A comprehensive transcriptome survey, or gene atlas, provides information essential for a complete understanding of the genomic biology of an organism. We present an atlas of RNA abundance for 92 adult, juvenile and fetal cattle tissues and three cattle cell lines. The Bovine Gene Atlas was generated from 7.2 million unique digital gene expression tag sequences (300.2 million total raw tag sequences), from which 1.59 million unique tag sequences were identified that mapped to the draft bovine genome accounting for 85% of the total raw tag abundance. Filtering these tags yielded 87,764 unique tag sequences that unambiguously mapped to 16,517 annotated protein-coding loci in the draft genome accounting for 45% of the total raw tag abundance. Clustering of tissues based on tag abundance profiles generally confirmed ontology classification based on anatomy. There were 5,429 constitutively expressed loci and 3,445 constitutively expressed unique tag sequences mapping outside annotated gene boundaries that represent a resource for enhancing current gene models. Physical measures such as inferred transcript length or antisense tag abundance identified tissues with atypical transcriptional tag profiles. We report for the first time the tissue-specific variation in the proportion of mitochondrial transcriptional tag abundance. The Bovine Gene Atlas is the deepest and broadest transcriptome survey of any livestock genome to date. Commonalities and variation in sense and antisense transcript tag profiles identified in different tissues facilitate the examination of the relationship between gene expression, tissue, and gene function.

Background: WC1 co-receptors belong to the scavenger receptor cysteine-rich (SRCR) superfamily and are encoded by a multi-gene family. Expression of particular WC1 genes defines functional subpopulations of WC1⁺γδ T cells. We have previously identified partial or complete genomic sequences for thirteen different WC1 genes through annotation of the bovine genome Btau_3.1 build. We also identified two WC1 cDNA sequences from other cattle that did not correspond to sequences in the Btau_3.1 build. Their absence in the Btau_3.1 build may have reflected gaps in the genome assembly or polymorphisms among animals. Since the response of γδ T cells to bacterial challenge is determined by WC1 gene expression, it was critical to understand whether individual cattle or breeds differ in the number of WC1 genes or display polymorphisms.Results: Real-time quantitative PCR using DNA from the animal whose genome was sequenced (“Dominette”) and sixteen other animals representing ten breeds of cattle, showed that the number of genes coding for WC1 co-receptors is thirteen. The complete coding sequences of those thirteen WC1 genes is presented, including the correction of an error in theWC1-2gene due to mis-assembly in the Btau_3.1 build. All other cDNA sequences were found to agree with the previous annotation of complete or partial WC1 genes. PCR amplification and sequencing of the most variable N-terminal SRCR domain (domain 1 which has the SRCR “a” pattern) of each of the thirteen WC1 genes showed that the sequences are highly conserved among individuals and breeds. Of 160 sequences of domain 1 from three breeds of cattle, no additional sequences beyond the thirteen described WC1 genes were found. Analysis of the complete WC1 cDNA sequences indicated that the thirteen WC1 genes code for three distinct WC1 molecular forms.Conclusion: The bovine WC1 multi-gene family is composed of thirteen genes coding for three structural forms whose sequences are highly conserved among individual cattle and breeds. The sequence diversity necessary for WC1 genes to function as a multi-genic pattern recognition receptor array is encoded in the genome, rather than generated by recombinatorial diversity or hypermutation.

A microsatellite-based genome scan of a Wagyu × Limousin F₂ cross population previously demonstrated QTL affecting LM area and fatty acid composition were present in regions near the centromere of BTA2. In this study, we used 70 SNP markers to examine the centromeric 24 megabases (Mb) of BTA2, including the Limousin-specific F94L myostatin allele (AB076403.1; 415C > A) located at approximately 6 Mb on the draft genome sequence of BTA2. A significant effect of the F94L marker was observed (F = 60.17) for LM area, which indicated that myostatin is most likely responsible for the effect. This is consistent with previous reports that the substitution of Leu for Phe at AA 94 of myostatin (caused by the 415C > A transversion) is associated with increased muscle growth. Surprisingly, several fatty acid trait QTL, which affected the amount of unsaturated fats, also mapped to or very near the myostatin marker, including the ratio of C16:1 MUFA to C16:0 saturated fat (F = 16.72), C18:1 to C18:0 (F = 18.88), and total content of MUFA (F = 17.12). In addition, QTL for extent of marbling (F = 14.73) approached significance (P = 0.05), and CLA concentration (F = 9.22) was marginally significant (P = 0.18). We also observed associations of SNP located at 16.3 Mb with KPH (F = 15.00) and for the amount of SFA (F = 12.01). These results provide insight into genetic differences between the Wagyu and Limousin breeds and may lead to a better tasting and healthier product for consumers through improved selection for lipid content of beef. Copyright © 2009. . Copyright 2009 Journal of Animal Science