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

Genetic Regulation of Liver Metabolites and Transcripts Linking to Biochemical-Clinical Parameters

Given the central metabolic role of the liver, hepatic metabolites and transcripts reflect the organismal physiological state. Biochemical-clinical plasma biomarkers, hepatic metabolites, transcripts, and single nucleotide polymorphism (SNP) genotypes of some 300 pigs were integrated by weighted correlation networks and genome-wide association analyses. Network-based approaches of transcriptomic and metabolomics data revealed linked of transcripts and metabolites of the pentose phosphate pathway (PPP). This finding was evidenced by using a NADP/NADPH assay and HDAC4 and G6PD transcript quantification with the latter coding for first limiting enzyme of this pathway and by RNAi knockdown experiments of HDAC4. Other transcripts including ARG2 and SLC22A7 showed link to amino acids and biomarkers. The amino acid metabolites were linked with transcripts of immune or acute phase response signaling, whereas the carbohydrate metabolites were highly enrich in cholesterol biosynthesis transcripts. Genome-wide association analyses revealed 180 metabolic quantitative trait loci (mQTL) (p < 10-4). Trans-4-hydroxy-L-proline (p = 6 × 10-9), being strongly correlated with plasma creatinine (CREA), showed strongest association with SNPs on chromosome 6 that had pleiotropic effects on PRODH2 expression as revealed by multivariate analysis. Consideration of shared marker association with biomarkers, metabolites, and transcripts revealed 144 SNPs associated with 44 metabolites and 69 transcripts that are correlated with each other, representing 176 mQTL and expression quantitative trait loci (eQTL). This is the first work to report genetic variants associated with liver metabolite and transcript levels as well as blood biochemical-clinical parameters in a healthy porcine model. The identified associations provide links between variation at the genome, transcriptome, and metabolome level molecules with clinically relevant phenotypes. This approach has the potential to detect novel biomarkers displaying individual variation and promoting predictive biology in medicine and animal breeding

Identification of novel transcripts and noncoding RNAs in bovine skin by deep next generation sequencing

BACKGROUND: 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. RESULTS: A total of 39,577 unique primary skin transcripts were mapped to the bovine reference genome assembly. The majority of the transcripts were mapped to known transcriptional units (65%). In addition to the reannotation of known genes, a substantial number (10,884) of unknown transcripts (UTs) were discovered, which had not previously been annotated. The classification of UTs was based on the prediction of their coding potential and comparative sequence analysis, subsequently followed by meticulous manual curation. The classification analysis and experimental validation of selected UTs confirmed that RNAseq data can be used to amend the annotation of known genes by providing evidence for additional exons, untranslated regions or splice variants, by approving genes predicted in silico and by identifying novel bovine loci. A large group of UTs (4,848) was predicted to potentially represent long noncoding RNA (lncRNA). Predominantly, potential lncRNAs mapped in intergenic chromosome regions (4,365) and therefore, were classified as potential intergenic lncRNA. Our analysis revealed that only about 6% of all UTs displayed interspecies conservation and discovered a variety of unknown transcripts without interspecies homology but specific expression in bovine skin. CONCLUSIONS: The results of our study demonstrate a complex transcript pattern for bovine skin and suggest a possible functional relevance of novel transcripts, including lncRNA, in the modulation of pigmentation processes. The results also indicate that the comprehensive identification and annotation of unknown transcripts from whole transcriptome analysis using RNAseq data remains a tremendous future challenge

Different Blood Cell-Derived Transcriptome Signatures in Cows Exposed to Vaccination Pre- or Postpartum

Periparturient cows have been found to reveal immunosuppression, frequently associated with increased susceptibility to uterine and mammary infections. To improve understanding of the causes and molecular regulatory mechanisms accounting for this phenomenon around calving, we examined the effect of an antigen challenge on gene expression modulation on cows prior to (BC) or after calving (AC) using whole transcriptome sequencing (RNAseq). The transcriptome analysis of the cows’ blood identified a substantially higher number of loci affected in BC cows (2,235) in response to vaccination compared to AC cows (208) and revealed a divergent transcriptional profile specific for each group. In BC cows, a variety of loci involved in immune defense and cellular signaling processes were transcriptionally activated, whereas protein biosynthesis and posttranslational processes were tremendously impaired in response to vaccination. Furthermore, energy metabolism in the blood cells of BC cows was shifted from oxidative phosphorylation to the glycolytic system. In AC cows, the number and variety of regulated pathways involved in immunomodulation and maintenance of immnunocompetence are considerably lower after vaccination, and upregulation of arginine degradation was suggested as an immunosuppressive mechanism. Elevated transcript levels of erythrocyte-specific genes involved in gas exchange processes were a specific transcriptional signature in AC cows pointing to hematopoiesis activation. The divergent and substantially lower magnitude of transcriptional modulation in response to vaccination in AC cows provides evidence for a suppressed immune capacity of early lactating cows on the molecular level and demonstrates that an efficient immune response of cows is related to their physiological and metabolic status

Unmapped reads from cattle RNAseq data: A source for missing and misassembled sequences in the reference assemblies and for detection of pathogens in the host

Usually, reads from transcriptome sequencing data unmapped to the target species' reference genome are disregarded. A recent RNAseq project on the new fatal disease Bovine Neonatal Pancytopenia had indicated an unexplained immune response signature to a double-stranded RNA virus. To unravel its background, contigs were de novo assembled from unmapped RNAseq reads and aligned against the bovine genome assemblies and multispecies NCBI databases. Lack of genuine virus sequence contigs rejected the hypothesis of a live virus being causal for the unexplained immune response. Alignment data also demonstrated incomplete bovine reference genome assemblies. In addition, we found that several parasite and virus genome reference assemblies in NCBI were contaminated with bovine DNA and confirmed recombination of bovine DNA into BVD virus strains. Exploring unmapped reads can extract useful biological information regarding the presence of microorganisms and can highlight issues with reference genome assemblies of host and pathogen species

A novel RNAseq–assisted method for MHC class I genotyping in a non-model species applied to a lethal vaccination-induced alloimmune disease

BACKGROUND: MHC class I genotyping is essential for a wide range of biomedical, immunological and biodiversity applications. Whereas in human a comprehensive MHC class I allele catalogue is available, respective data in non-model species is scarce in spite of decades of research. RESULTS: Taking advantage of the new high-throughput RNA sequencing technology (RNAseq), we developed a novel RNAseq-assisted method (RAMHCIT) for MHC class I typing at nucleotide level. RAMHCIT is performed on white blood cells, which highly express MHC class I molecules enabling reliable discovery of new alleles and discrimination of closely related alleles due to the high coverage of alleles with reads. RAMHCIT is more comprehensive than previous methods, because no targeted PCR pre-amplification of MHC loci is necessary, which avoids preselection of alleles as usually encountered, when amplification with MHC class I primers is performed prior to sequencing. In addition to allele identification, RAMHCIT also enables quantification of MHC class I expression at allele level, which was remarkably consistent across individuals. CONCLUSIONS: Successful application of RAMHCIT is demonstrated on a data set from cattle with different phenotype regarding a lethal, vaccination-induced alloimmune disease (bovine neonatal pancytopenia), for which MHC class I alleles had been postulated as causal agents. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s12864-016-2688-0) contains supplementary material, which is available to authorized users

Epigenome-wide skeletal muscle DNA methylation profiles at the background of distinct metabolic types and ryanodine receptor variation in pigs

BACKGROUND: Epigenetic variation may result from selection for complex traits related to metabolic processes or appear in the course of adaptation to mediate responses to exogenous stressors. Moreover epigenetic marks, in particular the DNA methylation state, of specific loci are driven by genetic variation. In this sense, polymorphism with major gene effects on metabolic and cell signaling processes, like the variation of the ryanodine receptors in skeletal muscle, may affect DNA methylation. METHODS: DNA-Methylation profiles were generated applying Reduced Representation Bisulfite Sequencing (RRBS) on 17 Musculus longissimus dorsi samples. We examined DNA methylation in skeletal muscle of pig breeds differing in metabolic type, Duroc and Pietrain. We also included F2 crosses of these breeds to get a first clue to DNA methylation sites that may contribute to breed differences. Moreover, we compared DNA methylation in muscle tissue of Pietrain pigs differing in genotypes at the gene encoding the Ca2+ release channel (RYR1) that largely affects muscle physiology. RESULTS: More than 2000 differently methylated sites were found between breeds including changes in methylation profiles of METRNL, IDH3B, COMMD6, and SLC22A18, genes involved in lipid metabolism. Depending on RYR1 genotype there were 1060 differently methylated sites including some functionally related genes, such as CABP2 and EHD, which play a role in buffering free cytosolic Ca2+ or interact with the Na+/Ca2+ exchanger. CONCLUSIONS: The change in the level of methylation between the breeds is probably the result of the long-term selection process for quantitative traits involving an infinite number of genes, or it may be the result of a major gene mutation that plays an important role in muscle metabolism and triggers extensive compensatory processes

Genetically regulated hepatic transcripts and pathways orchestrate haematological, biochemical and body composition traits

The liver is the central metabolic organ and exhibits fundamental functions in haematological traits. Hepatic expression, haematological, plasma biochemical, and body composition traits were assessed in a porcine model (n = 297) to establish tissue-specific genetic variations that influence the function of immune-metabolism-correlated expression networks. At FDR (false discovery rate) <1%, more than 3,600 transcripts were jointly correlated (r = |0.22–0.48|) with the traits. Functional enrichment analysis demonstrated common links of metabolic and immune traits. To understand how immune and metabolic traits are affected via genetic regulation of gene expression, eQTLs were assessed. 20517 significant (FDR < 5%) eQTLs for 1401 transcripts were identified, among which 443 transcripts were associated with at least one of the examined traits and had cis-eQTL (such as ACO1 (6.52 × 10−7) and SOD1 (6.41 × 10−30). The present study establishes a comprehensive view of hepatic gene activity which links together metabolic and immune traits in a porcine model for medical research

Mitochondrial-nuclear crosstalk, haplotype and copy number variation distinct in muscle fiber type, mitochondrial respiratory and metabolic enzyme activities

Genes expressed in mitochondria work in concert with those expressed in the nucleus to mediate oxidative phosphorylation (OXPHOS), a process that is relevant for muscle metabolism and meat quality. Mitochondrial genome activity can be efficiently studied and compared in Duroc and Pietrain pigs, which harbor different mitochondrial haplotypes and distinct muscle fiber types, mitochondrial respiratory activities, and fat content. Pietrain pigs homozygous-positive for malignant hyperthermia susceptibility (PiPP) carried only haplotype 8 and showed the lowest absolute mtDNA copy number accompanied by a decrease transcript abundance of mitochondrial-encoded subunits ND1, ND6, and ATP6 and nuclear-encoded subunits NDUFA11 and NDUFB8. In contrast, we found that haplotype 4 of Duroc pigs had significantly higher mitochondrial DNA (mtDNA) copy numbers and an increase transcript abundance of mitochondrial-encoded subunits ND1, ND6, and ATP6. These results suggest that the variation in mitochondrial and nuclear genetic background among these animals has an effect on mitochondrial content and OXPHOS system subunit expression. We observed the co-expression pattern of mitochondrial and nuclear encoded OXPHOS subunits suggesting that the mitochondrial-nuclear crosstalk functionally involves in muscle metabolism. The findings provide valuable information for understanding muscle biology processes and energy metabolism, and may direct use for breeding strategies to improve meat quality and animal health