Deep RNA sequencing (RNAseq) has opened a new horizon for understanding global gene expression. The functional annotation of non-model mammalian genomes including bovines is still poor compared to that of human and mouse. This particularly applies to tissues without direct significance for milk and meat production, like skin, in spite of its multifunctional relevance for the individual. Thus, applying an RNAseq approach, we performed a whole transcriptome analysis of pigmented and nonpigmented bovine skin to describe the comprehensive transcript catalogue of this tissue.

Several studies in a variety of breeds have reported at least two QTL for milk production traits, including milk fat synthesis on bovine chromosome 6 (BTA6), comprising a region that comparatively has been mapped to equivalent syntenic chromosome intervals in human, pig, and mouse harboring loci associated with type II diabetes and obesity-related traits. We identified the bovine peroxysome proliferator-activated receptor-γ coactivator-1α gene (PPARGC1A) as a plausible positional and functional candidate gene for a previously described QTL for milk fat yield on BTA6 because of its chromosomal position and its key role in energy, fat, and glucose metabolism. To analyze the role of the bovine PPARGC1A gene in regulation of milk fat synthesis in dairy cattle, we determined its cDNA sequence, genomic organization, chromosomal localization, and expression pattern. The bovine PPARGC1A gene is organized in 13 exons comprising 6,261 bp and is expressed at different levels in a large number of tissues. Bovine PPARGC1A cDNA and protein sequences showed substantial similarity (92–95%) to its respective orthologs from human, rat, and mouse. Screening for polymorphisms in the coding sequence, exon/intron boundaries, 5′- and 3′-untranslated regions, and promoter region of the PPARGC1A gene in sires with a different genotype at the QTL for milk fat yield as well as in a multibreed panel revealed a total of 11 polymorphic loci. A significant association between an SNP in intron 9 of the PPARGC1A gene and milk fat yield was observed in a major dairy cattle population, indicating that the PPARGC1A gene could be involved in genetic variation underlying the QTL for milk fat synthesis on BTA6.

Background: Systems biology enables the identification of gene networks that modulate complex traits. Comprehensive metabolomic analyses provide innovative phenotypes that are intermediate between the initiator of genetic variability, the genome, and raw phenotypes that are influenced by a large number of environmental effects. The present study combines two concepts, systems biology and metabolic analyses, in an approach without prior functional hypothesis in order to dissect genes and molecular pathways that modulate differential growth at the onset of puberty in male cattle. Furthermore, this integrative strategy was applied to specifically explore distinctive gene interactions of non-SMC condensin I complex, subunit G (NCAPG) and myostatin (GDF8), known modulators of pre- and postnatal growth that are only partially understood for their molecular pathways affecting differential body weight. Results: Our study successfully established gene networks and interacting partners affecting growth at the onset of puberty in cattle. We demonstrated the biological relevance of the created networks by comparison to randomly created networks. Our data showed that GnRH (Gonadotropin-releasing hormone) signaling is associated with divergent growth at the onset of puberty and revealed two highly connected hubs, BTC and DGKH, within the network. Both genes are known to directly interact with the GnRH signaling pathway. Furthermore, a gene interaction network for NCAPG containing 14 densely connected genes revealed novel information concerning the functional role of NCAPG in divergent growth. Conclusions: Merging both concepts, systems biology and metabolomic analyses, successfully yielded new insights into gene networks and interacting partners affecting growth at the onset of puberty in cattle. Genetic modulation in GnRH signaling was identified as key modifier of differential cattle growth at the onset of puberty. In addition, the benefit of our innovative concept without prior functional hypothesis was demonstrated by data suggesting that NCAPG might contribute to vascular smooth muscle contraction by indirect effects on the NO pathway via modulation of arginine metabolism. Our study shows for the first time in cattle that integration of genetic, physiological and metabolomics data in a systems biology approach will enable (or contribute to) an improved understanding of metabolic and gene networks and genotype-phenotype relationships.

The recently constructed river buffalo whole-genome radiation hybrid panel (BBURH₅₀₀₀) has already been used to generate preliminary radiation hybrid (RH) maps for several chromosomes, and buffalo-bovine comparative chromosome maps have been constructed. Here, we present the first-generation whole genome RH map (WG-RH) of the river buffalo generated from cattle-derived markers. The RH maps aligned to bovine genome sequence assembly Btau_4.0, providing valuable comparative mapping information for both species. A total of 3990 markers were typed on the BBURH₅₀₀₀ panel, of which 3072 were cattle derived SNPs. The remaining 918 were classified as cattle sequence tagged site (STS), including coding genes, ESTs, and microsatellites. Average retention frequency per chromosome was 27.3% calculated with 3093 scorable markers distributed in 43 linkage groups covering all autosomes (24) and the X chromosomes at a LOD ≥ 8. The estimated total length of the WG-RH map is 36,933 cR₅₀₀₀. Fewer than 15% of the markers (472) could not be placed within any linkage group at a LOD score ≥ 8. Linkage group order for each chromosome was determined by incorporation of markers previously assigned by FISH and by alignment with the bovine genome sequence assembly (Btau_4.0). We obtained radiation hybrid chromosome maps for the entire river buffalo genome based on cattle-derived markers. The alignments of our RH maps to the current bovine genome sequence assembly (Btau_4.0) indicate regions of possible rearrangements between the chromosomes of both species. The river buffalo represents an important agricultural species whose genetic improvement has lagged behind other species due to limited prior genomic characterization. We present the first-generation RH map which provides a more extensive resource for positional candidate cloning of genes associated with complex traits and also for large-scale physical mapping of the river buffalo genome.

The molecular mechanisms regulating the physiological adaptation of tissues important for nutrient partitioning and metabolism in lactating cows are still not completely understood. The aim of our study was to identify tissue-specific regulatory mechanisms necessary to accommodate metabolic changes associated with different genetic potential for milk performance. For this purpose, we analyzed mRNA expression of genes involved in energy metabolism of segregating F₂beef type cows with a combined genetic dairy and beef background (Charolais × German Holstein cross, CH×GH) in contrast to purebred German Holstein (GH) dairy cows. Three groups of cows differing in milk performance were examined using quantitative real-time PCR in liver, mammary gland, and skeletal muscle. Our results describe substantial tissue-specific differences in mRNA transcription profiles between cow groups in relation to their genetic potential for milk performance and highlight genes exhibiting specific, partially yet-unknown functions in dairy and beef type cows, e.g., upregulation of PCK2 transcripts in the mammary gland and FBP2 transcripts in skeletal muscle of dairy cows. Noticeably, PCCA and PPARGC1A mRNA abundance varied significantly across experimental groups in all three tissues, pointing to potential key gene functions in the metabolic adaptation relative to divergent milk production performance. Correlations of mRNA expression levels to milk performance traits indicate that gene transcriptional processes may play a regulatory role in liver, mammary gland, and skeletal muscle to enable cows with different genetic potential for milk performance to cope with metabolic lactation-associated challenges.

Background: The intramuscular fat deposition and the fatty acid profiles of beef affect meat quality. High proportions of unsaturated fatty acids are related to beef flavor and are beneficial for the nutritional value of meat. Moreover, a variety of clinical and epidemiologic studies showed that particularly long-chain omega-3 fatty acids from animal sources have a positive impact on human health and disease. Results: To screen for genetic factors affecting fatty acid profiles in beef, we initially performed a microsatellite-based genome scan in a F₂ Charolais × German Holstein resource population and identified a quantitative trait locus (QTL) for fatty acid composition in a region on bovine chromosome 27 where previously QTL affecting marbling score had been detected in beef cattle populations. The long-chain acyl-CoA synthetase 1 (ACSL1) gene was identified as the most plausible functional and positional candidate gene in the QTL interval due to its direct impact on fatty acid metabolism and its position in the QTL interval. ACSL1 is necessary for synthesis of long-chain acyl-CoA esters, fatty acid degradation and phospholipid remodeling. We validated the genomic annotation of the bovine ACSL1 gene by in silico comparative sequence analysis and experimental verification. Re-sequencing of the complete coding, exon-flanking intronic sequences, 3' untranslated region (3'UTR) and partial promoter region of the ACSL1 gene revealed three synonymous mutations in exons 6, 7, and 20, six noncoding intronic gene variants, six polymorphisms in the promoter region, and four variants in the 3' UTR region. The association analysis identified the gene variant in intron 5 of the ACSL1 gene (c.481-233A>G) to be significantly associated with the relative content of distinct fractions and ratios of fatty acids (e.g., n-3 fatty acids, polyunsaturated, n-3 long-chain polyunsaturated fatty acids, trans vaccenic acid) in skeletal muscle. A tentative association of the ACSL1 gene variant with intramuscular fat content indicated that an indirect effect on fatty acid composition via modulation of total fat content of skeletal muscle cannot be excluded. Conclusions: The initial QTL analysis suggested the ACSL1 gene as a positional and functional candidate gene for fatty acid composition in bovine skeletal muscle. The findings of subsequent association analyses indicate that ACSL1 or a separate gene in close proximity might play a functional role in mediating the lipid composition of beef.