Concurrently mapping quantitative trait loci associations from multiple subspecies within hybrid populations.

Many of the world's agriculturally important plant and animal populations consist of hybrids of subspecies. Cattle in tropical and sub-tropical regions for example, originate from two subspecies, Bos taurus indicus (Bos indicus) and Bos taurus taurus (Bos taurus). Methods to derive the underlying genetic architecture for these two subspecies are essential to develop accurate genomic predictions in these hybrid populations. We propose a novel method to achieve this. First, we use haplotypes to assign SNP alleles to ancestral subspecies of origin in a multi-breed and multi-subspecies population. Then we use a BayesR framework to allow SNP alleles originating from the different subspecies differing effects. Applying this method in a composite population of B. indicus and B. taurus hybrids, our results show that there are underlying genomic differences between the two subspecies, and these effects are not identified in multi-breed genomic evaluations that do not account for subspecies of origin effects. The method slightly improved the accuracy of genomic prediction. More significantly, by allocating SNP alleles to ancestral subspecies of origin, we were able to identify four SNP with high posterior probabilities of inclusion that have not been previously associated with cattle fertility and were close to genes associated with fertility in other species. These results show that haplotypes can be used to trace subspecies of origin through the genome of this hybrid population and, in conjunction with our novel Bayesian analysis, subspecies SNP allele allocation can be used to increase the accuracy of QTL association mapping in genetically diverse populations.

Use of whole-genome sequence data and novel genomic selection strategies to improve selection for age at puberty in tropically-adapted beef heifers.

BACKGROUND: In tropically-adapted beef heifers, application of genomic prediction for age at puberty has been limited due to low prediction accuracies. Our aim was to investigate novel methods of pre-selecting whole-genome sequence (WGS) variants and alternative analysis methodologies; including genomic best linear unbiased prediction (GBLUP) with multiple genomic relationship matrices (MGRM) and Bayesian (BayesR) analyses, to determine if prediction accuracy for age at puberty can be improved. METHODS: Genotypes and phenotypes were obtained from two research herds. In total, 868 Brahman and 960 Tropical Composite heifers were recorded in the first population and 3695 Brahman, Santa Gertrudis and Droughtmaster heifers were recorded in the second population. Genotypes were imputed to 23 million whole-genome sequence variants. Eight strategies were used to pre-select variants from genome-wide association study (GWAS) results using conditional or joint (COJO) analyses. Pre-selected variants were included in three models, GBLUP with a single genomic relationship matrix (SGRM), GBLUP MGRM and BayesR. Five-way cross-validation was used to test the effect of marker panel density (6 K, 50 K and 800 K), analysis model, and inclusion of pre-selected WGS variants on prediction accuracy. RESULTS: In all tested scenarios, prediction accuracies for age at puberty were highest in BayesR analyses. The addition of pre-selected WGS variants had little effect on the accuracy of prediction when BayesR was used. The inclusion of WGS variants that were pre-selected using a meta-analysis with COJO analyses by chromosome, fitted in a MGRM model, had the highest prediction accuracies in the GBLUP analyses, regardless of marker density. When the low-density (6 K) panel was used, the prediction accuracy of GBLUP was equal (0.42) to that with the high-density panel when only six additional sequence variants (identified using meta-analysis COJO by chromosome) were included. CONCLUSIONS: While BayesR consistently outperforms other methods in terms of prediction accuracies, reasonable improvements in accuracy can be achieved when using GBLUP and low-density panels with the inclusion of a relatively small number of highly relevant WGS variants.

Genetic control of temperament traits across species: association of autism spectrum disorder risk genes with cattle temperament.

BACKGROUND: Temperament traits are of high importance across species. In humans, temperament or personality traits correlate with psychological traits and psychiatric disorders. In cattle, they impact animal welfare, product quality and human safety, and are therefore of direct commercial importance. We hypothesized that genetic factors that contribute to variation in temperament among individuals within a species will be shared between humans and cattle. Using imputed whole-genome sequence data from 9223 beef cattle from three cohorts, a series of genome-wide association studies was undertaken on cattle flight time, a temperament phenotype measured as the time taken for an animal to cover a short-fixed distance after release from an enclosure. We also investigated the association of cattle temperament with polymorphisms in bovine orthologs of risk genes for neuroticism, schizophrenia, autism spectrum disorders (ASD), and developmental delay disorders in humans. RESULTS: Variants with the strongest associations were located in the bovine orthologous region that is involved in several behavioural and cognitive disorders in humans. These variants were also partially validated in independent cattle cohorts. Genes in these regions (BARHL2, NDN, SNRPN, MAGEL2, ABCA12, KIFAP3, TOPAZ1, FZD3, UBE3A, and GABRA5) were enriched for the GO term neuron migration and were differentially expressed in brain and pituitary tissues in humans. Moreover, variants within 100 kb of ASD susceptibility genes were associated with cattle temperament and explained 6.5% of the total additive genetic variance in the largest cattle cohort. The ASD genes with the most significant associations were GABRB3 and CUL3. Using the same 100 kb window, a weak association was found with polymorphisms in schizophrenia risk genes and no association with polymorphisms in neuroticism and developmental delay disorders risk genes. CONCLUSIONS: Our analysis showed that genes identified in a meta-analysis of cattle temperament contribute to neuron development functions and are differentially expressed in human brain tissues. Furthermore, some ASD susceptibility genes are associated with cattle temperament. These findings provide evidence that genetic control of temperament might be shared between humans and cattle and highlight the potential for future analyses to leverage results between species.

Adequate supply of dietary taurine stimulates expression of molecular markers of growth and protein turnover in juvenile barramundi (Lates calcarifer).

A trial was conducted to investigate the effect of dietary taurine (Tau) supply on the plasma amino acid composition and hepatic expression of several genes in juvenile barramundi (Lates calcarifer) after feeding. Triplicate tanks of fish (average weight, 89.3 g) were fed diets containing either a deficient (1 g kg-1), adequate (8 g kg-1) or excessive (19 g kg-1) level of dietary Tau. Liver tissues collected before feeding, and at 2- and 4-h post-feeding, were analysed for expression of genes involved in pathways of sulphur amino acid turnover, Tau biosynthesis and transport, target of rapamycin (TOR) signalling, the somatotropic axis and protein turnover. The treatment had no significant effect on the profiles of any amino acid in plasma collected over time after feeding, other than Tau and glycine. The expression profile of cystine and Tau synthetic genes suggested an effect of Tau excess on the metabolism of cystine. Markers of two pathways of Tau biosynthesis appear to be active in this species, providing proof that this species possesses the ability to synthesise Tau from SAA precursors. A marker for the regulation of Tau transport and homeostasis was shown to be directly regulated by Tau availability, whilst a link between adequate supply of Tau and TOR pathway-mediated growth stimulation was also apparent. An observed depression in expression of genes of the somatotropic axis, coupled with upregulation of the proteolytic and TOR-suppressing genes, in response to excessive Tau supply in the diet, signalled that excessive Tau may not be conducive to optimal growth in this species.

Genome-wide association for milk production and female fertility traits in Canadian dairy Holstein cattle.

BACKGROUND: Genome-wide association studies (GWAS) are a powerful tool for detecting genomic regions explaining variation in phenotype. The objectives of the present study were to identify or refine the positions of genomic regions affecting milk production, milk components and fertility traits in Canadian Holstein cattle, and to use these positions to identify genes and pathways that may influence these traits. RESULT: Several QTL regions were detected for milk production (MILK), fat production (FAT), protein production (PROT) and fat and protein deviation (FATD, PROTD respectively). The identified QTL regions for production traits (including milk production) support previous findings and some overlap with genes with known relevant biological functions identified in earlier studies such as DGAT1 and CPSF1. A significant region on chromosome 21 overlapping with the gene FAM181A and not previous linked to fertility in dairy cattle was identified for the calving to first service interval and days open. A functional enrichment analysis of the GWAS results yielded GO terms consistent with the specific phenotypes tested, for example GO terms GO:0007595 (lactation) and GO:0043627 (response to estrogen) for milk production (MILK), GO:0051057 (positive regulation of small GTPase mediated signal transduction) for fat production (FAT), GO:0040019 (positive regulation of embryonic development) for first service to calving interval (CTFS) and GO:0043268 (positive regulation of potassium ion transport) for days open (DO). In other cases the connection between the enriched GO terms and the traits were less clear, for example GO:0003279 (cardiac septum development) for FAT and GO:0030903 (notochord development) for DO trait. CONCLUSION: The chromosomal regions and enriched pathways identified in this study confirm several previous findings and highlight new regions and pathways that may contribute to variation in production or fertility traits in dairy cattle.

Harnessing cross-species alignment to discover SNPs and generate a draft genome sequence of a bighorn sheep (Ovis canadensis).

BACKGROUND: Whole genome sequences (WGS) have proliferated as sequencing technology continues to improve and costs decline. While many WGS of model or domestic organisms have been produced, a growing number of non-model species are also being sequenced. In the absence of a reference, construction of a genome sequence necessitates de novo assembly which may be beyond the ability of many labs due to the large volumes of raw sequence data and extensive bioinformatics required. In contrast, the presence of a reference WGS allows for alignment which is more tractable than assembly. Recent work has highlighted that the reference need not come from the same species, potentially enabling a wide array of species WGS to be constructed using cross-species alignment. Here we report on the creation a draft WGS from a single bighorn sheep (Ovis canadensis) using alignment to the closely related domestic sheep (Ovis aries). RESULTS: Two sequencing libraries on SOLiD platforms yielded over 865 million reads, and combined alignment to the domestic sheep reference resulted in a nearly complete sequence (95% coverage of the reference) at an average of 12x read depth (104 SD). From this we discovered over 15 million variants and annotated them relative to the domestic sheep reference. We then conducted an enrichment analysis of those SNPs showing fixed differences between the reference and sequenced individual and found significant differences in a number of gene ontology (GO) terms, including those associated with reproduction, muscle properties, and bone deposition. CONCLUSION: Our results demonstrate that cross-species alignment enables the creation of novel WGS for non-model organisms. The bighorn sheep WGS will provide a resource for future resequencing studies or comparative genomics.

Non-synonymous mutations mapped to chromosome X associated with andrological and growth traits in beef cattle.

BACKGROUND: Previous genome-wide association analyses identified QTL regions in the X chromosome for percentage of normal sperm and scrotal circumference in Brahman and Tropical Composite cattle. These traits are important to be studied because they are indicators of male fertility and are correlated with female sexual precocity and reproductive longevity. The aim was to investigate candidate genes in these regions and to identify putative causative mutations that influence these traits. In addition, we tested the identified mutations for female fertility and growth traits. RESULTS: Using a combination of bioinformatics and molecular assay technology, twelve non-synonymous SNPs in eleven genes were genotyped in a cattle population. Three and nine SNPs explained more than 1% of the additive genetic variance for percentage of normal sperm and scrotal circumference, respectively. The SNPs that had a major influence in percentage of normal sperm were mapped to LOC100138021 and TAF7L genes; and in TEX11 and AR genes for scrotal circumference. One SNP in TEX11 was explained ~13% of the additive genetic variance for scrotal circumference at 12 months. The tested SNP were also associated with weight measurements, but not with female fertility traits. CONCLUSIONS: The strong association of SNPs located in X chromosome genes with male fertility traits validates the QTL. The implicated genes became good candidates to be used for genetic evaluation, without detrimentally influencing female fertility traits.

Evidence for positive selection of taurine genes within a QTL region on chromosome X associated with testicular size in Australian Brahman cattle.

BACKGROUND: Previous genome-wide association studies have identified significant regions of the X chromosome associated with reproductive traits in two Bos indicus-influenced breeds: Brahman cattle and Tropical Composites. Two QTL regions on this chromosome were identified in both breeds as strongly associated with scrotal circumference measurements, a reproductive trait previously shown to be useful for selection of young bulls. Scrotal circumference is genetically correlated with early age at puberty in both male and female offspring. These QTL were located at positions 69-77 and 81-92 Mb respectively, large areas each to which a significant number of potential candidate genes were mapped. RESULTS: To further characterise these regions, a bioinformatic approach was undertaken to identify novel non-synonymous SNP within the QTL regions of interest in Brahman cattle. After SNP discovery, we used conventional molecular assay technologies to perform studies of two candidate genes in both breeds. Non-synonymous SNP mapped to Testis-expressed gene 11 (Tex11) were associated (P < 0.001) with scrotal circumference in both breeds, and associations with percentage of normal sperm cells were also observed (P < 0.05). Evidence for recent selection was found as Tex11 SNP form a haplotype segment of Bos taurus origin that is retained within Brahman and Tropical Composite cattle with greatest reproductive potential. CONCLUSIONS: Association of non-synonymous SNP presented here are a first step to functional genetic studies. Bovine species may serve as a model for studying the role of Tex11 in male fertility, warranting further in-depth molecular characterisation.

Whole-genome resequencing of Hanwoo (Korean cattle) and insight into regions of homozygosity.

BACKGROUND: Hanwoo (Korean cattle), which originated from natural crossbreeding between taurine and zebu cattle, migrated to the Korean peninsula through North China. Hanwoo were raised as draft animals until the 1970s without the introduction of foreign germplasm. Since 1979, Hanwoo has been bred as beef cattle. Genetic variation was analyzed by whole-genome deep resequencing of a Hanwoo bull. The Hanwoo genome was compared to that of two other breeds, Black Angus and Holstein, and genes within regions of homozygosity were investigated to elucidate the genetic and genomic characteristics of Hanwoo. RESULTS: The Hanwoo bull genome was sequenced to 45.6-fold coverage using the ABI SOLiD system. In total, 4.7 million single-nucleotide polymorphisms and 0.4 million small indels were identified by comparison with the Btau4.0 reference assembly. Of the total number of SNPs and indels, 58% and 87%, respectively, were novel. The overall genotype concordance between the SNPs and BovineSNP50 BeadChip data was 96.4%. Of 1.6 million genetic differences in Hanwoo, approximately 25,000 non-synonymous SNPs, splice-site variants, and coding indels (NS/SS/Is) were detected in 8,360 genes. Among 1,045 genes containing reliable specific NS/SS/Is in Hanwoo, 109 genes contained more than one novel damaging NS/SS/I. Of the genes containing NS/SS/Is, 610 genes were assigned as trait-associated genes. Moreover, 16, 78, and 51 regions of homozygosity (ROHs) were detected in Hanwoo, Black Angus, and Holstein, respectively. 'Regulation of actin filament length' was revealed as a significant gene ontology term and 25 trait-associated genes for meat quality and disease resistance were found in 753 genes that resided in the ROHs of Hanwoo. In Hanwoo, 43 genes were located in common ROHs between whole-genome resequencing and SNP chips in BTA2, 10, and 13 coincided with quantitative trait loci for meat fat traits. In addition, the common ROHs in BTA2 and 16 were in agreement between Hanwoo and Black Angus. CONCLUSIONS: We identified 4.7 million SNPs and 0.4 million small indels by whole-genome resequencing of a Hanwoo bull. Approximately 25,000 non-synonymous SNPs, splice-site variants, and coding indels (NS/SS/Is) were detected in 8,360 genes. Additionally, we found 25 trait-associated genes for meat quality and disease resistance among 753 genes that resided in the ROHs of Hanwoo. These findings will provide useful genomic information for identifying genes or casual mutations associated with economically important traits in cattle.

Genome-wide copy number variation in Hanwoo, Black Angus, and Holstein cattle.

Hanwoo, Korean native cattle, is indigenous to the Korean peninsula. They have been used mainly as draft animals for about 5,000 years; however, in the last 30 years, their main role has been changed to meat production by selective breeding which has led to substantial increases in their productivity. Massively parallel sequencing technology has recently made possible the systematic identification of structural variations in cattle genomes. In particular, copy number variation (CNV) has been recognized as an important genetic variation complementary to single-nucleotide polymorphisms that can be used to account for variations of economically important traits in cattle. Here we report genome-wide copy number variation regions (CNVRs) in Hanwoo cattle obtained by comparing the whole genome sequence of Hanwoo with Black Angus and Holstein sequence datasets. We identified 1,173 and 963 putative CNVRs representing 16.7 and 7.8 Mbp from comparisons between Black Angus and Hanwoo and between Holstein and Hanwoo, respectively. The potential functional roles of the CNVRs were assessed by Gene Ontology enrichment analysis. The results showed that response to stimulus, immune system process, and cellular component organization were highly enriched in the genic-CNVRs that overlapped with annotated cattle genes. Of the 11 CNVRs that were selected for validation by quantitative real-time PCR, 9 exhibited the expected copy number differences. The results reported in this study show that genome-wide CNVs were detected successfully using massively parallel sequencing technology. The CNVs may be a valuable resource for further studies to correlate CNVs and economically important traits in cattle.

Impact of feed efficiency and diet on adaptive variations in the bacterial community in the rumen fluid of cattle.

Limited knowledge of the structure and activities of the ruminal bacterial community prevents the understanding of the effect of population dynamics on functional bacterial groups and on host productivity. This study aimed to identify particular bacteria associated with host feed efficiency in steers with differing diets and residual feed intake (RFI) using culture-independent methods: PCR-denaturing gradient gel electrophoresis (DGGE) and quantitative real-time PCR analysis. PCR-DGGE profiles were generated from the ruminal fluid of 55 steers fed a low-energy-density diet and then switched to a high-energy-density diet. Bacterial profile comparisons by multivariate statistical analysis showed a trend only for RFI-related clusters on the high-energy diet. When steers (n = 19) belonging to the same RFI group under both diets were used to identify specific bacterial phylotypes related to feed efficiency traits, correlations were detected between dry matter intake, average daily gain, and copy numbers of the 16S rRNA gene of Succinivibrio sp. in low-RFI (efficient) steers, whereas correlations between Robinsoniella sp. and RFI (P < 0.05) were observed for high-RFI (inefficient) animals. Eubacterium sp. differed significantly (P < 0.05) between RFI groups that were only on the high-energy diet. Our work provides a comprehensive framework to understand how particular bacterial phylotypes contribute to differences in feed efficiency and ultimately influence host productivity, which may either depend on or be independent from diet factors.

Genome-wide analysis of zygotic linkage disequilibrium and its components in crossbred cattle.

BACKGROUND: Linkage disequilibrium (LD) between genes at linked or independent loci can occur at gametic and zygotic levels known asgametic LD and zygotic LD, respectively. Gametic LD is well known for its roles in fine-scale mapping of quantitative trait loci, genomic selection and evolutionary inference. The less-well studied is the zygotic LD and its components that can be also estimated directly from the unphased SNPs. RESULTS: This study was set up to investigate the genome-wide extent and patterns of zygotic LD and its components in a crossbred cattle population using the genomic data from the Illumina BovineSNP50 beadchip. The animal population arose from repeated crossbreeding of multiple breeds and selection for growth and cow reproduction. The study showed that similar genomic structures in gametic and zygotic LD were observed, with zygotic LD decaying faster than gametic LD over marker distance. The trigenic and quadrigenic disequilibria were generally two- to three-fold smaller than the usual digenic disequilibria (gametic or composite LD). There was less power of testing for these high-order genic disequilibria than for the digenic disequilibria. The power estimates decreased with the marker distance between markers though the decay trend is more obvious for the digenic disequilibria than for high-order disequilibria. CONCLUSIONS: This study is the first major genome-wide survey of all non-allelic associations between pairs of SNPs in a cattle population. Such analysis allows us to assess the relative importance of gametic LD vs. all other non-allelic genic LDs regardless of whether or not the population is in HWE. The observed predominance of digenic LD (gametic or composite LD) coupled with insignificant high-order trigenic and quadrigenic disequilibria supports the current intensive focus on the use of high-density SNP markers for genome-wide association studies and genomic selection activities in the cattle population.

Short reads, circular genome: skimming solid sequence to construct the bighorn sheep mitochondrial genome.

As sequencing technology improves, an increasing number of projects aim to generate full genome sequence, even for nonmodel taxa. These projects may be feasibly conducted at lower read depths if the alignment can be aided by previously developed genomic resources from a closely related species. We investigated the feasibility of constructing a complete mitochondrial (mt) genome without preamplification or other targeting of the sequence. Here we present a full mt genome sequence (16,463 nucleotides) for the bighorn sheep (Ovis canadensis) generated though alignment of SOLiD short-read sequences to a reference genome. Average read depth was 1240, and each base was covered by at least 36 reads. We then conducted a phylogenomic analysis with 27 other bovid mitogenomes, which placed bighorn sheep firmly in the Ovis clade. These results show that it is possible to generate a complete mitogenome by skimming a low-coverage genomic sequencing library. This technique will become increasingly applicable as the number of taxa with some level of genome sequence rises.

Functional genomics approach for identification of molecular processes underlying neurodegenerative disorders in prion diseases.

Prion diseases or transmissible spongiform encephalopathies (TSEs) are infectious neurodegenerative disorders leading to death. These include Cresutzfeldt-Jakob disease (CJD), familial, sporadic and variant CJD and kuru in humans; and animal TSEs include scrapie in sheep, bovine spongiform encephalopathy (BSE) in cattle, chronic wasting disease (CWD) of mule deer and elk, and transmissible mink encephalopathy. All these TSEs share common pathological features such as accumulation of mis-folded prion proteins in the central nervous system leading to cellular dysfunction and cell death. It is important to characterize the molecular pathways and events leading to prion induced neurodegeneration. Here we discuss the impact of the functional genomics approaches including microarrays, subtractive hybridization and microRNA profiling in elucidating transcriptional cascades at different stages of disease. Many of these transcriptional changes have been observed in multiple neurodegenerative diseases which may aid in identification of biomarkers for disease. A comprehensive characterization of expression profiles implicated in neurodegenerative disorders will undoubtedly advance our understanding on neuropathology and dysfunction during prion disease and other neurodegenerative disorders. We also present an outlook on the future work which may focus on analysis of structural genetic variation, genome and transcriptome sequencing using next generation sequencing with an integrated approach on animal and human TSE related studies.

Extensive Variation in Gene Expression is Revealed in 13 Fertility-Related Genes Using RNA-Seq, ISO-Seq, and CAGE-Seq From Brahman Cattle.

Fertility is a key driver of economic profitability in cattle production. A number of studies have identified genes associated with fertility using genome wide association studies and differential gene expression analysis; however, the genes themselves are poorly characterized in cattle. Here, we selected 13 genes from the literature which have previously been shown to have strong evidence for an association with fertility in Brahman cattle (Bos taurus indicus) or closely related breeds. We examine the expression variation of the 13 genes that are associated with cattle fertility using RNA-seq, CAGE-seq, and ISO-seq data from 11 different tissue samples from an adult Brahman cow and a Brahman fetus. Tissues examined include blood, liver, lung, kidney, muscle, spleen, ovary, and uterus from the cow and liver and lung from the fetus. The analysis revealed several novel isoforms, including seven from SERPINA7. The use of three expression characterization methodologies (5' cap selected ISO-seq, CAGE-seq, and RNA-seq) allowed the identification of isoforms that varied in their length of 5' and 3' untranslated regions, variation otherwise undetectable (collapsed as degraded RNA) in generic isoform identification pipelines. The combinations of different sequencing technologies allowed us to overcome the limitations of relatively low sequence depth in the ISO-seq data. The lower sequence depth of the ISO-seq data was also reflected in the lack of observed expression of some genes that were observed in the CAGE-seq and RNA-seq data from the same tissue. We identified allele specific expression that was tissue-specific in AR, IGF1, SOX9, STAT3, and TAF9B. Finally, we characterized an exon of TAF9B as partially nested within the neighboring gene phosphoglycerate kinase 1. As this study only examined two animals, even more transcriptional variation may be present in a genetically diverse population. This analysis reveals the large amount of transcriptional variation within mammalian fertility genes and illuminates the fact that the transcriptional landscape cannot be fully characterized using a single technology alone.

Whole genome resequencing of black Angus and Holstein cattle for SNP and CNV discovery.

BACKGROUND: One of the goals of livestock genomics research is to identify the genetic differences responsible for variation in phenotypic traits, particularly those of economic importance. Characterizing the genetic variation in livestock species is an important step towards linking genes or genomic regions with phenotypes. The completion of the bovine genome sequence and recent advances in DNA sequencing technology allow for in-depth characterization of the genetic variations present in cattle. Here we describe the whole-genome resequencing of two Bos taurus bulls from distinct breeds for the purpose of identifying and annotating novel forms of genetic variation in cattle. RESULTS: The genomes of a Black Angus bull and a Holstein bull were sequenced to 22-fold and 19-fold coverage, respectively, using the ABI SOLiD system. Comparisons of the sequences with the Btau4.0 reference assembly yielded 7 million single nucleotide polymorphisms (SNPs), 24% of which were identified in both animals. Of the total SNPs found in Holstein, Black Angus, and in both animals, 81%, 81%, and 75% respectively are novel. In-depth annotations of the data identified more than 16 thousand distinct non-synonymous SNPs (85% novel) between the two datasets. Alignments between the SNP-altered proteins and orthologues from numerous species indicate that many of the SNPs alter well-conserved amino acids. Several SNPs predicted to create or remove stop codons were also found. A comparison between the sequencing SNPs and genotyping results from the BovineHD high-density genotyping chip indicates a detection rate of 91% for homozygous SNPs and 81% for heterozygous SNPs. The false positive rate is estimated to be about 2% for both the Black Angus and Holstein SNP sets, based on follow-up genotyping of 422 and 427 SNPs, respectively. Comparisons of read depth between the two bulls along the reference assembly identified 790 putative copy-number variations (CNVs). Ten randomly selected CNVs, five genic and five non-genic, were successfully validated using quantitative real-time PCR. The CNVs are enriched for immune system genes and include genes that may contribute to lactation capacity. The majority of the CNVs (69%) were detected as regions with higher abundance in the Holstein bull. CONCLUSIONS: Substantial genetic differences exist between the Black Angus and Holstein animals sequenced in this work and the Hereford reference sequence, and some of this variation is predicted to affect evolutionarily conserved amino acids or gene copy number. The deeply annotated SNPs and CNVs identified in this resequencing study can serve as useful genetic tools, and as candidates in searches for phenotype-altering DNA differences.

SNP discovery and allele frequency estimation by deep sequencing of reduced representation libraries.

High-density single-nucleotide polymorphism (SNP) arrays have revolutionized the ability of genome-wide association studies to detect genomic regions harboring sequence variants that affect complex traits. Extensive numbers of validated SNPs with known allele frequencies are essential to construct genotyping assays with broad utility. We describe an economical, efficient, single-step method for SNP discovery, validation and characterization that uses deep sequencing of reduced representation libraries (RRLs) from specified target populations. Using nearly 50 million sequences generated on an Illumina Genome Analyzer from DNA of 66 cattle representing three populations, we identified 62,042 putative SNPs and predicted their allele frequencies. Genotype data for these 66 individuals validated 92% of 23,357 selected genome-wide SNPs, with a genotypic and sequence allele frequency correlation of r = 0.67. This approach for simultaneous de novo discovery of high-quality SNPs and population characterization of allele frequencies may be applied to any species with at least a partially sequenced genome.

Gene expression in the medulla following oral infection of cattle with bovine spongiform encephalopathy.

The identification of variations in gene expression in response to bovine spongiform encephalopathy (BSE) may help to elucidate the mechanisms of neuropathology and prion replication and discover biomarkers for disease. In this study, genes that are differentially expressed in the caudal medulla tissues of animals infected with different doses of PrP(BSE) at 12 and 45 mo post infection were compared using array containing 24,000 oligonucleotide probes. Data analysis identified 966 differentially expressed (DE) genes between control and infected animals. Genes identified in at least two of four experiments (control versus 1-g infected animals at 12 and 45-mo; control versus 100-g infected animals at 12 and 45 mo) were considered to be the genes that may be associated with BSE disease. From the 176 DE genes associated with BSE, 84 had functions described in the Gene Ontology (GO) database. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of 14 genes revealed that prion infection may cause dysfunction of several different networks, including extracellular matrix (ECM), cell adhesion, neuroactive ligand-receptor interaction, complement and coagulation cascades, MAPK signaling, neurodegenerative disorder, SNARE interactions in vesicular transport, and the transforming growth factor (TGF) beta signaling pathways. The identification of DE genes will contribute to a better understanding of the molecular mechanisms of neuropathology in bovine species. Additional studies on larger number of animals are in progress in our laboratory to investigate the roles of these DE genes in pathogenesis of BSE.

Transcriptome analysis of the medulla tissue from cattle in response to bovine spongiform encephalopathy using digital gene expression tag profiling.

Bovine spongiform encephalopathy (BSE) is a transmissible, fatal neurodegenerative disorder of cattle produced by prions. The use of excessive parallel sequencing for comparison of gene expression in bovine control and infected tissues may help to elucidate the molecular mechanisms associated with this disease. In this study, tag profiling Solexa sequencing was used for transcriptome analysis of bovine brain tissues. Replicate libraries were prepared from mRNA isolated from control and infected (challenged with 100 g of BSE-infected brain) medulla tissues 45 mo after infection. For each library, 5-6 million sequence reads were generated and approximately 67-70% of the reads were mapped against the Bovine Genome database to approximately 13,700-14,120 transcripts (each having at least one read). About 42-47% of the total reads mapped uniquely. Using the GeneSifter software package, 190 differentially expressed (DE) genes were identified (>2.0-fold change, p < .01): 73 upregulated and 117 downregulated. Seventy-nine DE genes had functions described in the Gene Ontology (GO) database and 16 DE genes were involved in 38 different pathways described in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Digital analysis expression by tag profiling may be a powerful approach to comprehensive transcriptome analysis to identify changes associated with disease progression, leading to a better understanding of the underlying mechanism of pathogenesis of BSE.

A principal component regression based genome wide analysis approach reveals the presence of a novel QTL on BTA7 for MAP resistance in holstein cattle.

Bovine Johne's disease (JD), caused by Mycobacterium avium spp. paratuberculosis (MAP), causes significant losses to the dairy and beef cattle industries. Effective vaccination or therapeutic strategies against this disease are currently unavailable and infected animals either get culled or die due to clinical disease. An alternative strategy to manage the disease is to selectively breed animals with enhanced resistance to MAP infection. Therefore, the objective of this study was to identify genetic loci putatively associated with MAP infection in a resource population consisting of Holstein cattle using a genome-wide association approach. The BovineSNP50 BeadChip, containing 54,001 single nucleotide polymorphisms (SNPs), was used to genotype 232 animals with known MAP infection status. Since, traditional case-control analytical techniques are based on single-marker analysis and do not account for the existence of linkage disequilibrium (LD) between markers, we used a novel principal component regression approach, where each SNP was fit in a logistic regression model, along with principal components of other SNPs on the same chromosome showing association with the trait, as covariates. Such an approach allowed us to account for the LD that exists between multiple markers showing an association on the same chromosome. Our analysis revealed the presence of at least 12 genomic regions on BTA1, 5, 6, 7, 10, 11 and 14 that were associated with the MAP infection status of our resource population. A brief description of these genomic regions, and a discussion of the analysis used in this study, have been presented.