Different mitochondrial DNA copy number in liver and mammary gland of lactating cows with divergent genetic background for milk production.

Adequate metabolic adaptation of key tissues playing an essential role for bioenergetic homeostasis and lactogenesis is critical in cows to adapt to changes in energy requirements and physiological processes during the lactation period. Mitochondria are recognized as central to meet energy needs and maintaining of metabolic homeostasis because mitochondrial DNA (mtDNA) is template for several polypeptides of the respiratory chain complexes essential for ATP generation. The quantity of mtDNA in a cell has been widely used as a surrogate marker for the capacity of cells for energy generation. In our study we analyzed the mtDNA copy number and the mRNA expression of important nuclear encoded genes controlling mitochondrial biogenesis in liver and mammary gland. We compared cows with a nuclear genome dairy × beef crossbred make-up to purebred German Holstein dairy cows. The study revealed tissue-specific variations of mtDNA copy number and expression levels of nuclear genes involved in mitochondrial biogenesis when comparing lactating cows with different genetic predisposition regarding milk performance. This may reflect nuclear genome-determined genetic differences between the cow groups in coping with metabolic demands and physiological changes during lactation. The results indicate that mitochondrial biogenesis processes in the liver and mammary gland appear to be impaired in high lactating dairy cows, which consequently, would point to a disturbed energy adaptation. The results provide a basis to further elucidate the adaptive and regulatory modulation of the mitochondrial biogenesis in response to lactation-associated metabolic challenges in lactating cows.

Epistatic interactions between at least three loci determine the "rat-tail" phenotype in cattle.

BACKGROUND: The "rat-tail" syndrome (RTS) is an inherited hypotrichosis in cattle, which is exclusively expressed in diluted coloured hair. The affected animals also suffer from disturbed thermoregulation, which impairs their health and growth performance. Phenotypic features that are similar to RTS are observed in dogs with black hair follicle dysplasia. RESULTS: We used a resource cross population between German Holstein and Charolais cattle breeds to prove that epistatic interactions between at least three independent genetic loci are required for the expression of the RTS phenotype. In this population, the RTS is exclusively expressed in animals with a eumelanic background that is due to the dominant E (D) allele at the melanocortin 1 receptor gene located on Bos taurus autosome (BTA) 18. In addition, only the individuals that are heterozygous at the dilution locus on BTA5 that corresponds to the premelanosome protein or silver gene variant c.64G>A were classified as displaying a RTS phenotype. Linkage and whole-genome association analyses using different models and different pedigrees allowed us to map a third locus (hereafter referred to as the RTS locus) that is essential for the expression of the RTS phenotype to the chromosomal region between 14 and 22 Mb on BTA5. Our findings clearly demonstrate that the RTS and dilution loci are distinct loci on BTA5. CONCLUSIONS: Our study provides evidence that the RTS locus has effects on hair conformation and coat colour dilution and that the effect on coat colour dilution is clearly independent from that of the dilution locus. Finally, our results excluded several other loci that were previously reported to be associated with or to underlie hair conformation or pigmentation traits as the causal mutations of RTS and also several major functional candidate genes that are associated with hypotrichosis in humans. Our finding on the identification of a three-locus interaction that underlies RTS provides a prime example of epistatic interaction between several independent loci that is required for the expression of a distinct phenotype.

Systems biology analysis merging phenotype, metabolomic and genomic data identifies Non-SMC Condensin I Complex, Subunit G (NCAPG) and cellular maintenance processes as major contributors to genetic variability in bovine feed efficiency.

Feed efficiency is a paramount factor for livestock economy. Previous studies had indicated a substantial heritability of several feed efficiency traits. In our study, we investigated the genetic background of residual feed intake, a commonly used parameter of feed efficiency, in a cattle resource population generated from crossing dairy and beef cattle. Starting from a whole genome association analysis, we subsequently performed combined phenotype-metabolome-genome analysis taking a systems biology approach by inferring gene networks based on partial correlation and information theory approaches. Our data about biological processes enriched with genes from the feed efficiency network suggest that genetic variation in feed efficiency is driven by genetic modulation of basic processes relevant to general cellular functions. When looking at the predicted upstream regulators from the feed efficiency network, the Tumor Protein P53 (TP53) and Transforming Growth Factor beta 1 (TGFB1) genes stood out regarding significance of overlap and number of target molecules in the data set. These results further support the hypothesis that TP53 is a major upstream regulator for genetic variation of feed efficiency. Furthermore, our data revealed a significant effect of both, the Non-SMC Condensin I Complex, Subunit G (NCAPG) I442M (rs109570900) and the Growth /differentiation factor 8 (GDF8) Q204X (rs110344317) loci, on residual feed intake and feed conversion. For both loci, the growth promoting allele at the onset of puberty was associated with a negative, but favorable effect on residual feed intake. The elevated energy demand for increased growth triggered by the NCAPG 442M allele is obviously not fully compensated for by an increased efficiency in converting feed into body tissue. As a consequence, the individuals carrying the NCAPG 442M allele had an additional demand for energy uptake that is reflected by the association of the allele with increased daily energy intake as observed in our study.

A systems biology approach using metabolomic data reveals genes and pathways interacting to modulate divergent growth in cattle.

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.

Multiple splice variants within the bovine silver homologue (SILV) gene affecting coat color in cattle indicate a function additional to fibril formation in melanophores.

BACKGROUND: The silver homologue(SILV) gene plays a major role in melanosome development. SILV is a target for studies concerning melanoma diagnostics and therapy in humans as well as on skin and coat color pigmentation in many species ranging from zebra fish to mammals. However, the precise functional cellular mechanisms, in which SILV is involved, are still not completely understood. While there are many studies addressing SILV function upon a eumelaneic pigment background, there is a substantial lack of information regarding the further relevance of SILV, e.g. for phaeomelanosome development. RESULTS: In contrast to previous results in other species reporting SILV expression exclusively in pigmented tissues, our experiments provide evidence that the bovine SILV gene is expressed in a variety of tissues independent of pigmentation. Our data show that the bovine SILV gene generates an unexpectedly large number of different transcripts occurring in skin as well as in non-pigmented tissues, e.g. liver or mammary gland. The alternative splice sites are generated by internal splicing and primarily remove complete exons. Alternative splicing predominantly affects the repeat domain of the protein, which has a functional key role in fibril formation during eumelanosome development. CONCLUSION: The expression of the bovine SILV gene independent of pigmentation suggests SILV functions exceeding melanosome development in cattle. This hypothesis is further supported by transcript variants lacking functional key elements of the SILV protein relevant for eumelanosome development. Thus, the bovine SILV gene can serve as a model for the investigation of the putative additional functions of SILV. Furthermore, the splice variants of the bovine SILV gene represent a comprehensive natural model to refine the knowledge about functional domains in the SILV protein. Our study exemplifies that the extent of alternative splicing is presumably much higher than previously estimated and that alternatively spliced transcripts presumably can generate molecules of deviating function compared to their constitutive counterpart.

A gene-based high-resolution comparative radiation hybrid map as a framework for genome sequence assembly of a bovine chromosome 6 region associated with QTL for growth, body composition, and milk performance traits.

BACKGROUND: A number of different quantitative trait loci (QTL) for various phenotypic traits, including milk production, functional, and conformation traits in dairy cattle as well as growth and body composition traits in meat cattle, have been mapped consistently in the middle region of bovine chromosome 6 (BTA6). Dense genetic and physical maps and, ultimately, a fully annotated genome sequence as well as their mutual connections are required to efficiently identify genes and gene variants responsible for genetic variation of phenotypic traits. A comprehensive high-resolution gene-rich map linking densely spaced bovine markers and genes to the annotated human genome sequence is required as a framework to facilitate this approach for the region on BTA6 carrying the QTL. RESULTS: Therefore, we constructed a high-resolution radiation hybrid (RH) map for the QTL containing chromosomal region of BTA6. This new RH map with a total of 234 loci including 115 genes and ESTs displays a substantial increase in loci density compared to existing physical BTA6 maps. Screening the available bovine genome sequence resources, a total of 73 loci could be assigned to sequence contigs, which were already identified as specific for BTA6. For 43 loci, corresponding sequence contigs, which were not yet placed on the bovine genome assembly, were identified. In addition, the improved potential of this high-resolution RH map for BTA6 with respect to comparative mapping was demonstrated. Mapping a large number of genes on BTA6 and cross-referencing them with map locations in corresponding syntenic multi-species chromosome segments (human, mouse, rat, dog, chicken) achieved a refined accurate alignment of conserved segments and evolutionary breakpoints across the species included. CONCLUSION: The gene-anchored high-resolution RH map (1 locus/300 kb) for the targeted region of BTA6 presented here will provide a valuable platform to guide high-quality assembling and annotation of the currently existing bovine genome sequence draft to establish the final architecture of BTA6. Hence, a sequence-based map will provide a key resource to facilitate prospective continued efforts for the selection and validation of relevant positional and functional candidates underlying QTL for milk production and growth-related traits mapped on BTA6 and on similar chromosomal regions from evolutionary closely related species like sheep and goat. Furthermore, the high-resolution sequence-referenced BTA6 map will enable precise identification of multi-species conserved chromosome segments and evolutionary breakpoints in mammalian phylogenetic studies.