Sean McWilliam

A Genome Wide Survey of SNP Variation Reveals the Genetic Structure of Sheep Breeds

The genetic structure of sheep reflects their domestication and subsequent formation into discrete breeds. Understanding genetic structure is essential for achieving genetic improvement through genome-wide association studies, genomic selection and the dissection of quantitative traits. After identifying the first genome-wide set of SNP for sheep, we report on levels of genetic variability both within and between a diverse sample of ovine populations. Then, using cluster analysis and the partitioning of genetic variation, we demonstrate sheep are characterised by weak phylogeographic structure, overlapping genetic similarity and generally low differentiation which is consistent with their short evolutionary history. The degree of population substructure was, however, sufficient to cluster individuals based on geographic origin and known breed history. Specifically, African and Asian populations clustered separately from breeds of European origin sampled from Australia, New Zealand, Europe and North America. Furthermore, we demonstrate the presence of stratification within some, but not all, ovine breeds. The results emphasize that careful documentation of genetic structure will be an essential prerequisite when mapping the genetic basis of complex traits. Furthermore, the identification of a subset of SNP able to assign individuals into broad groupings demonstrates even a small panel of markers may be suitable for applications such as traceability.

The sheep genome illuminates biology of the rumen and lipid metabolism

A genome for ewe and ewe Sheep-specific genetic changes underlie differences in lipid metabolism between sheep and other mammals, and may have contributed to the production of wool. Jiang et al. sequenced the genome of two Texel sheep, a breed that produces high-value meat, milk, and wool. The genome information will provide an important resource for livestock production and aid in the understanding of mammalian evolution. Science, this issue p. 1168 A genomic analysis of sheep explains specializations in digestive system physiology and wool production. Sheep (Ovis aries) are a major source of meat, milk, and fiber in the form of wool and represent a distinct class of animals that have a specialized digestive organ, the rumen, that carries out the initial digestion of plant material. We have developed and analyzed a high-quality reference sheep genome and transcriptomes from 40 different tissues. We identified highly expressed genes encoding keratin cross-linking proteins associated with rumen evolution. We also identified genes involved in lipid metabolism that had been amplified and/or had altered tissue expression patterns. This may be in response to changes in the barrier lipids of the skin, an interaction between lipid metabolism and wool synthesis, and an increased role of volatile fatty acids in ruminants compared with nonruminant animals.

Gene expression profiling of muscle tissue in Brahman steers during nutritional restriction

Expression profiling using microarrays allows for the detailed characterization of the gene networks that regulate an animals response to environmental stresses. During nutritional restriction, processes such as protein turnover, connective tissue remodeling, and muscle atrophy take place in the skeletal muscle of the animal. These processes and their regulation are of interest in the context of managing livestock for optimal production efficiency and product quality. Here we expand on recent research applying complementary DNA (cDNA) microarray technology to the study of the effect of nutritional restriction on bovine skeletal muscle. Using a custom cDNA microarray of 9,274 probes from cattle muscle and s.c. fat libraries, we examined the differential gene expression profile of the LM from 10 Brahman steers under three different dietary treatments. The statistical approach was based on mixed-model ANOVA and model-based clustering of the BLUP solutions for the gene x diet interaction effect. From the results, we defined a transcript profile of 156 differentially expressed array elements between the weight loss and weight gain diet substrates. After sequence and annotation analyses, the 57 upregulated elements represented 29 unique genes, and the 99 downregulated elements represented 28 unique genes. Most of these co-regulated genes cluster into groups with distinct biological function related to protein turnover and cytoskeletal metabolism and contribute to our mechanistic understanding of the processes associated with remodeling of muscle tissue in response to nutritional stress.

Origin, Evolution, and Biological Role of miRNA Cluster in DLK-DIO3 Genomic Region in Placental Mammals

MicroRNAs (miRNAs) are a rapidly growing family of small regulatory RNAs modulating gene expression in plants and animals. In animals, most of the miRNAs discovered in early studies were found to be evolutionarily conserved across the whole kingdom. More recent studies, however, have identified many miRNAs that are specific to a particular group of organisms or even a single species. These present a question about evolution of the individual miRNAs and their role in establishing and maintaining lineage-specific functions and characteristics. In this study, we describe a detailed analysis of the miRNA cluster (hereafter mir-379/mir-656 cluster) located within the imprinted DLK-DIO3 region on human chromosome 14. We show that orthologous miRNA clusters are present in all sequenced genomes of the placental (eutherian) mammals but not in the marsupial (metatherian), monotreme (prototherian), or any other vertebrate genomes. We provide evidence that the locus encompassing this cluster emerged in an early eutherian ancestor prior to the radiation of modern placental mammals by tandem duplication of the ancient precursor sequence. The original amplified cluster may have contained in excess of 250 miRNA precursor sequences, most of which now appear to be inactive. Examination of the eutherian genomes showed that the cluster has been maintained in evolution for approximately 100 Myr. Analysis of genes that contain predicted evolutionarily conserved targets for miRNAs from this cluster revealed significant overrepresentation of the Gene Ontology terms associated with biological processes such as neurogenesis, embryonic development, transcriptional regulation, and RNA metabolism. Consistent with these findings, a survey of the miRNA expression data within the cluster demonstrates a strong bias toward brain and placenta samples from adult organisms and some embryonic tissues. Our results suggest that emergence of the mir-379/mir-656 miRNA cluster was one of the factors that facilitated evolution of the placental mammals. Overrepresentation of genes involved in regulation of neurogenesis among predicted miRNAs targets indicates an important role of the mir-379/mir-656 cluster in this biological process in the placental mammals.

Using comparative genomics to reorder the human genome sequence into a virtual sheep genome

BackgroundIs it possible to construct an accurate and detailed subgene-level map of a genome using bacterial artificial chromosome (BAC) end sequences, a sparse marker map, and the sequences of other genomes?ResultsA sheep BAC library, CHORI-243, was constructed and the BAC end sequences were determined and mapped with high sensitivity and low specificity onto the frameworks of the human, dog, and cow genomes. To maximize genome coverage, the coordinates of all BAC end sequence hits to the cow and dog genomes were also converted to the equivalent human genome coordinates. The 84,624 sheep BACs (about 5.4-fold genome coverage) with paired ends in the correct orientation (tail-to-tail) and spacing, combined with information from sheep BAC comparative genome contigs (CGCs) built separately on the dog and cow genomes, were used to construct 1,172 sheep BAC-CGCs, covering 91.2% of the human genome. Clustered non-tail-to-tail and outsize BACs located close to the ends of many BAC-CGCs linked BAC-CGCs covering about 70% of the genome to at least one other BAC-CGC on the same chromosome. Using the BAC-CGCs, the intrachromosomal and interchromosomal BAC-CGC linkage information, human/cow and vertebrate synteny, and the sheep marker map, a virtual sheep genome was constructed. To identify BACs potentially located in gaps between BAC-CGCs, an additional set of 55,668 sheep BACs were positioned on the sheep genome with lower confidence. A coordinate conversion process allowed us to transfer human genes and other genome features to the virtual sheep genome to display on a sheep genome browser.ConclusionWe demonstrate that limited sequencing of BACs combined with positioning on a well assembled genome and integrating locations from other less well assembled genomes can yield extensive, detailed subgene-level maps of mammalian genomes, for which genomic resources are currently limited.

An interactive bovine in silico SNP database (IBISS)

An interactive bovine in silico SNP (IBISS) database has been created through the clustering and aligning of bovine EST and mRNA sequences. Approximately 324,000 EST and mRNA sequences were clustered to produce 29,965 clusters (producing 48,679 consensus sequences) and 48,565 singletons. A SNP screening regime was placed on variations detected in the multiple sequence alignment files to determine which SNPs are more likely to be real rather than sequencing errors. A small subset of predicted SNPs was validated on a diverse set of bovine DNA samples using PCR amplification and sequencing. Fifty percent of the predicted SNPs in the “putative >1” category were polymorphic in the population sampled. The IBISS database represents more than just a SNP database; it is also a genomic database containing uniformly annotated predicted gene mRNA and protein sequences, gene structure, and genomic organization information.

The genome of bovine ephemeral fever rhabdovirus contains two related glycoprotein genes.

A 3789 nucleotide region of the bovine ephemeral fever virus (BEFV) genome, located 1.65 kb downstream of the N gene, has been cloned and sequenced. The region contains two long open reading frames (ORFs) which are bounded by putative consensus (AACAGG) and polyadenylation (CATG[A]7) sequences and are separated by an intergenic region of 53 nucleotides. Discrete mRNAs corresponding to each ORF have been identified. The first ORF encodes a polypeptide comprising 623 residues which was identified by peptide sequencing as the virion G protein. The deduced amino acid sequence of the G protein includes putative signal and transmembrane domains and five potential glycosylation sites. The second ORF encodes a polypeptide of 586 amino acids which also has characteristics of a rhabdovirus glycoprotein, including putative signal and transmembrane domains and eight potential glycosylation sites, and appears to correspond to a 90-kDa nonstructural glycoprotein (GNS) identified in BEFV-infected cells (Walker et al. [1991] J. Gen. Virol. 72, 67-74). A database search indicated that both the G and GNS proteins share significant amino acid sequence homology with other rhabdovirus G proteins and with each other. Highest homology scores for each protein were with sigma virus and vesicular stomatitis virus serotypes.

Gene expression profiling of bovine in vitro adipogenesis using a cDNA microarray

The gene expression profile of bovine bone marrow stromal cells undergoing adipogenesis was established using a custom cDNA microarray. Cells that were treated with adipogenic stimulants and those that were not were collected at each of the six time points, and gene expression differences between the treated and untreated samples within each time point were compared using a microarray. Statistical analyses revealed that 158 genes showed a minimum fold change of 2 in at least one of the five post-differentiation time points. These genes are involved in various cellular pathways and functions, including lipogenesis, glycolysis, cytoskeleton remodelling, extracellular matrix, transcription as well as various signalling pathways such as insulin, calcium and wingless signalling. The experiment also identified 17 differentially expressed (DE) microarray elements with no assigned function. Quantitative real-time PCR was employed to validate eight DE genes, and the PCR data were found to reproduce the microarray data for these eight genes. Subsequent gene ontology annotation was able to provide a global overview of the molecular function of DE genes during adipogenesis. This analysis was able to indicate the importance of different gene categories at various stages of adipogenic conversion, thereby providing further insights into the molecular changes during bovine adipogenesis.

GENOME ORGANIZATION AND TRANSCRIPTION STRATEGY IN THE COMPLEX GNS-L INTERGENIC REGION OF BOVINE EPHEMERAL FEVER RHABDOVIRUS

A 1622 nucleotide region of the bovine ephemeral fever virus (BEFV) genome, located between the second glycoprotein (GNS) gene and the polymerase (L) gene, has been cloned and sequenced in Australian (BB7721) and Chinese (Beijing-1) isolates of the virus. In the Australian isolate, the region contains five long open reading frames (ORFs) organized into three coding regions (alpha, beta and gamma), each of which are bound by a consensus transcription initiation and transcription termination-polyadenylation-like sequences. The alpha coding region contains three long ORFs (alpha 1, alpha 2 and alpha 3). The alpha 1 ORF encodes a 10.6 kDa polypeptide which contains hydrophobic and highly basic regions characteristic of a viroporin. The alpha 2 ORF encodes a 13.7 kDa polypeptide and overlaps the alpha 3 ORF which encodes a 5.7 kDa polypeptide. The beta coding region contains a single long ORF encoding a polypeptide of 12.2 kDa. The gamma coding region, which does not occur in Adelaide River virus (ARV), contains a single long ORF encoding a polypeptide of 13.4 kDa. The Chinese isolate shares 91% nucleotide sequence identity with the Australian isolate. The organization of the alpha, beta and gamma coding regions is preserved and the sequences of the encoded polypeptides are similar to those of BB7721. The major transcription products of the region were identified in BB7721 as polycistronic alpha (alpha 1-alpha 2-alpha 3) and beta-gamma mRNAs. Sequence similarities in the BEFV alpha-beta and beta-gamma gene junctions, and the gamma-L and beta-L gene junctions of BEFV and ARV, suggest that the gamma gene may have evolved from the beta-gene by sequence duplication.

Validation of alternative methods of data normalization in gene co-expression studies

MOTIVATION Clusters of genes encoding proteins with related functions, or in the same regulatory network, often exhibit expression patterns that are correlated over a large number of conditions. Protein associations and gene regulatory networks can be modelled from expression data. We address the question of which of several normalization methods is optimal prior to computing the correlation of the expression profiles between every pair of genes. RESULTS We use gene expression data from five experiments with a total of 78 hybridizations and 23 diverse conditions. Nine methods of data normalization are explored based on all possible combinations of normalization techniques according to between and within gene and experiment variation. We compare the resulting empirical distribution of gene x gene correlations with the expectations and apply cross-validation to test the performance of each method in predicting accurate functional annotation. We conclude that normalization methods based on mixed-model equations are optimal.