Allele-biased expression of the bovine APOB gene associated with the cholesterol deficiency defect suggests cis-regulatory enhancer effects of the LTR retrotransposon insertion.

The insertion of an endogenous retroviral long terminal repeat (LTR) sequence into the bovine apolipoprotein B (APOB) gene is causal to the inherited genetic defect cholesterol deficiency (CD) observed in neonatal and young calves. Affected calves suffer from developmental abnormalities, symptoms of incurable diarrhoea and often die within weeks to a few months after birth. Neither the detailed effects of the LTR insertion on APOB expression profile nor the specific mode of inheritance nor detailed phenotypic consequences of the mutation are undisputed. In our study, we analysed German Holstein dairy heifers at the peak of hepatic metabolic load and exposed to an additional pathogen challenge for clinical, metabolic and hepatic transcriptome differences between wild type (CDF) and heterozygote carriers of the mutation (CDC). Our data revealed that a divergent allele-biased expression pattern of the APOB gene in heterozygous CDC animals leads to a tenfold higher expression of exons upstream and a decreased expression of exons downstream of the LTR insertion compared to expression levels of CDF animals. This expression pattern could be a result of enhancer activity induced by the LTR insertion, in addition to a previously reported artificial polyadenylation signal. Thus, our data support a regulatory potential of mobile element insertions. With regard to the phenotype generated by the LTR insertion, heterozygote CDC carriers display significantly differential hepatic expression of genes involved in cholesterol biosynthesis and lipid metabolism. Phenotypically, CDC carriers show a significantly affected lipomobilization compared to wild type animals. These results reject a completely recessive mode of inheritance for the CD defect, which should be considered for selection decisions in the affected population. Exemplarily, our results illustrate the regulatory impact of mobile element insertions not only on specific host target gene expression but also on global transcriptome profiles with subsequent biological, functional and phenotypic consequences in a natural in-vivo model of a non-model mammalian organism.

Metabogenomic analysis to functionally annotate the regulatory role of long non-coding RNAs in the liver of cows with different nutrient partitioning phenotype.

Long non-coding RNAs (lncRNAs) hold gene regulatory potential, but require substantial further functional annotation in livestock. Applying two metabogenomic approaches by combining transcriptomic and metabolomic analyses, we aimed to identify lncRNAs with potential regulatory function for divergent nutrient partitioning of lactating crossbred cows and to establish metabogenomic interaction networks comprising metabolites, genes and lncRNAs. Through correlation analysis of lncRNA expression with transcriptomic and metabolomic data, we unraveled lncRNAs that have a putative regulatory role in energy and lipid metabolism, the urea and tricarboxylic acid cycles, and gluconeogenesis. Especially FGF21, which correlated with a plentitude of differentially expressed genes, differentially abundant metabolites, as well as lncRNAs, suggested itself as a key metabolic regulator. Notably, lncRNAs in close physical proximity to coding-genes as well as lncRNAs with natural antisense transcripts appear to perform a fine-tuning function in gene expression involved in metabolic pathways associated with different nutrient partitioning phenotypes.

Single-cell RNA sequencing of freshly isolated bovine milk cells and cultured primary mammary epithelial cells.

Bovine mammary function at molecular level is often studied using mammary tissue or primary bovine mammary epithelial cells (pbMECs). However, bulk tissue and primary cells are heterogeneous with respect to cell populations, adding further transcriptional variation in addition to genetic background. Thus, understanding of the variation in gene expression profiles of cell populations and their effect on function are limited. To investigate the mononuclear cell composition in bovine milk, we analyzed a single-cell suspension from a milk sample. Additionally, we harvested cultured pbMECs to characterize gene expression in a homogeneous cell population. Using the Drop-seq technology, we generated single-cell RNA datasets of somatic milk cells and pbMECs. The final datasets after quality control filtering contained 7,119 and 10,549 cells, respectively. The pbMECs formed 14 indefinite clusters displaying intrapopulation heterogeneity, whereas the milk cells formed 14 more distinct clusters. Our datasets constitute a molecular cell atlas that provides a basis for future studies of milk cell composition and gene expression, and could serve as reference datasets for milk cell analysis.

Identification of Regulatory Functions of LncRNAs Associated With T. circumcincta Infection in Adult Sheep.

Several recent studies have demonstrated the role of long non-coding RNAs (lncRNAs) in regulating the defense mechanism against parasite infections, but no studies are available that investigated their relevance for immune response to nematode infection in sheep. Thus, the aim of the current study was to (i) detect putative lncRNAs that are expressed in the abomasal lymph node of adult sheep after an experimental infection with the gastrointestinal nematode (GIN) Teladorsagia circumcincta and (ii) to elucidate their potential functional role associated with the differential host immune response. We hypothesized that putative lncRNAs differentially expressed (DE) between samples from animals that differ in resistance to infection may play a significant regulatory role in response to nematode infection in adult sheep. To obtain further support for our hypothesis, we performed co-expression and functional gene enrichment analyses with the differentially expressed lncRNAs (DE lncRNAs). In a conservative approach, we included for this predictive analysis only those lncRNAs that are confirmed and supported by documentation of expression in gastrointestinal tissues in the current sheep gene atlas. We identified 9,105 putative lncRNA transcripts corresponding to 7,124 gene loci. Of these, 457 were differentially expressed lncRNA loci (DELs) with 683 lncRNA transcripts. Based on a gene co-expression analysis via weighted gene co-expression network analysis, 12 gene network modules (GNMs) were found significantly correlated with at least one of 10 selected target DE lncRNAs. Based on the principle of "guilt-by-association," the DE genes from each of the three most significantly correlated GNMs were subjected to a gene enrichment analysis. The significant pathways associated with DE lncRNAs included ERK5 Signaling, SAPK/JNK Signaling, RhoGDI Signaling, EIF2 Signaling, Regulation of eIF4 and p70S6K Signaling and Oxidative Phosphorylation pathways. They belong to signaling pathway categories like Cellular Growth, Proliferation and Development, Cellular Stress and Injury, Intracellular and Second Messenger Signaling and Apoptosis. Overall, this lncRNA study conducted in adult sheep after GIN infection provided first insights into the potential functional role of lncRNAs in the differential host response to nematode infection.

A 50-kb deletion disrupting the RSPO2 gene is associated with tetradysmelia in Holstein Friesian cattle.

BACKGROUND: Tetradysmelia is a rare genetic disorder that is characterized by an extremely severe reduction of all limb parts distal of the scapula and pelvic girdle. We studied a Holstein Friesian backcross family with 24 offspring, among which six calves displayed autosomal recessive tetradysmelia. In order to identify the genetic basis of the disorder, we genotyped three affected calves, five dams and nine unaffected siblings using a Bovine Illumina 50 k BeadChip and sequenced the whole genome of the sire. RESULTS: Pathological examination of four tetradysmelia cases revealed a uniform and severe dysmelia of all limbs. Applying a homozygosity mapping approach, we identified a homozygous region of 10.54 Mb on chromosome 14 (Bos taurus BTA14). Only calves that were diagnosed with tetradysmelia shared a distinct homozygous haplotype for this region. We sequenced the whole genome of the cases' sire and searched for heterozygous single nucleotide polymorphisms (SNPs) and small variants on BTA14 that were uniquely present in the sire and absent from 3102 control whole-genome sequences of the 1000 Bull Genomes Project, but none were identified in the 10.54-Mb candidate region on BTA14. Therefore, we subsequently performed a more comprehensive analysis by also considering structural variants and detected a 50-kb deletion in the targeted chromosomal region that was in the heterozygous state in the cases' sire. Using PCR, we confirmed that this detected deletion segregated perfectly within the family with tetradysmelia. The deletion spanned three exons of the bovine R-spondin 2 (RSPO2) gene, which encode three domains of the respective protein. R-spondin 2 is a secreted ligand of leucine-rich repeats containing G protein-coupled receptors that enhance Wnt signalling and is involved in a broad range of developmental processes during embryogenesis. CONCLUSIONS: We identified a 50-kb deletion on BTA14 that disrupts the coding sequence of the RSPO2 gene and is associated with bovine tetradysmelia. To our knowledge, this is the first reported candidate causal mutation for tetradysmelia in a large animal model. Since signalling pathways involved in limb development are conserved across species, the observed inherited defect may serve as a model to further elucidate fundamental pathways of limb development.

Indication of Premelanosome Protein (PMEL) Expression Outside of Pigmented Bovine Skin Suggests Functions Beyond Eumelanogenesis.

The premelanosome protein (PMEL) is important for fibril formation within melanosomes during vertebrate melanogenesis. Fibrils form a matrix for pigment deposition within pigmented tissues such as skin and hair. PMEL mutations are known to modulate eumelanic pigmentation in vertebrates. However, in bovines, PMEL mutations were also found to alter pheomelanic pigmentation resulting in coat color dilution. Furthermore, epistatic effects of a mutated PMEL allele were detected in the phenotypic expression of the bovine hair defect "rat-tail syndrome" (RTS) characterized by charcoal coat color and hair deformation. Reports about PMEL gene expression in non-pigmented tissues raised the hypothesis that there may be unknown functions of the PMEL protein beyond eumelanin deposition to PMEL fibrils. In our study, we analysed the PMEL protein expression in pigmented skin and non-pigmented bovine tissues (non-pigmented skin, thyroid gland, rumen, liver, kidney, and adrenal gland cortex). We found that a processed form of the bovine PMEL protein is expressed in pigmented as well as in non-pigmented tissues, which is in line with gene expression data from targeted RT-PCR and whole transcriptome RNAseq analysis. The PMEL protein is located in membranes and within the cytosol of epithelial cells. Based on our data from bovine tissues, we concluded that at least in cattle PMEL potentially has additional, yet unexplored functions, which might contribute to effects of PMEL mutations on pheomelanin coat color dilution and charcoal coat color in RTS animals. However, indication of PMEL protein in unpigmented cells and tissues will require further confirmation in the future, because there have been no confirmed reports before, which had detected bovine PMEL protein with specific antibodies either in pigmented or unpigmented tissue.

Identification and Annotation of Potential Function of Regulatory Antisense Long Non-Coding RNAs Related to Feed Efficiency in Bos taurus Bulls.

Long non-coding RNAs (lncRNAs) can influence transcriptional and translational processes in mammalian cells and are associated with various developmental, physiological and phenotypic conditions. However, they remain poorly understood and annotated in livestock species. We combined phenotypic, metabolomics and liver transcriptomic data of bulls divergent for residual feed intake (RFI) and fat accretion. Based on a project-specific transcriptome annotation for the bovine reference genome ARS-UCD.1.2 and multiple-tissue total RNA sequencing data, we predicted 3590 loci to be lncRNAs. To identify lncRNAs with potential regulatory influence on phenotype and gene expression, we applied the regulatory impact factor algorithm on a functionally prioritized set of loci (n = 4666). Applying the algorithm of partial correlation and information theory, significant and independent pairwise correlations were calculated and co-expression networks were established, including plasma metabolites correlated with lncRNAs. The network hub lncRNAs were assessed for potential cis-actions and subjected to biological pathway enrichment analyses. Our results reveal a prevalence of antisense lncRNAs positively correlated with adjacent protein-coding genes and suggest their participation in mitochondrial function, acute phase response signalling, TCA-cycle, fatty acid ß-oxidation and presumably gluconeogenesis. These antisense lncRNAs indicate a stabilizing function for their cis-correlated genes and a putative regulatory role in gene expression.

Hepatic Transcriptome Analysis Identifies Divergent Pathogen-Specific Targeting-Strategies to Modulate the Innate Immune System in Response to Intramammary Infection.

Mastitis is one of the major risks for public health and animal welfare in the dairy industry. Two of the most important pathogens to cause mastitis in dairy cattle are Staphylococcus aureus (S. aureus) and Escherichia coli (E. coli). While S. aureus generally induces a chronic and subclinical mastitis, E. coli is an important etiological pathogen resulting in an acute and clinical mastitis. The liver plays a central role in both, the metabolic and inflammatory physiology of the dairy cow, which is particularly challenged in the early lactation due to high metabolic and immunological demands. In the current study, we challenged the mammary glands of Holstein cows with S. aureus or E. coli, respectively, mimicking an early lactation infection. We compared the animals' liver transcriptomes with those of untreated controls to investigate the hepatic response of the individuals. Both, S. aureus and E. coli elicited systemic effects on the host after intramammary challenge and seemed to use pathogen-specific targeting strategies to bypass the innate immune system. The most striking result of our study is that we demonstrate for the first time that S. aureus intramammary challenge causes an immune response beyond the original local site of the mastitis. We found that in the peripheral liver tissue defined biological pathways are switched on in a coordinated manner to balance the immune response in the entire organism. TGFB1 signaling plays a crucial role in this context. Important pathways involving actin and integrin, key components of the cytoskeleton, were downregulated in the liver of S. aureus infected cows. In the hepatic transcriptome of E. coli infected cows, important components of the complement system were significantly lower expressed compared to the control cows. Notably, while S. aureus inhibits the cell signaling by Rho GTPases in the liver, E. coli switches the complement system off. Also, metabolic hepatic pathways (e.g., lipid metabolism) are affected after mammary gland challenge, demonstrating that the liver restricts metabolic tasks in favor of the predominant immune response after infection. Our results provide new insights for the infection-induced modifications of the dairy cow's hepatic transcriptome following mastitis.

Biological Network Approach for the Identification of Regulatory Long Non-Coding RNAs Associated With Metabolic Efficiency in Cattle.

Background: Genomic regions associated with divergent livestock feed efficiency have been found predominantly outside protein coding sequences. Long non-coding RNAs (lncRNA) can modulate chromatin accessibility, gene expression and act as important metabolic regulators in mammals. By integrating phenotypic, transcriptomic, and metabolomic data with quantitative trait locus data in prioritizing co-expression network analyses, we aimed to identify and functionally characterize lncRNAs with a potential key regulatory role in metabolic efficiency in cattle. Materials and Methods: Crossbred animals (n = 48) of a Charolais x Holstein F2-population were allocated to groups of high or low metabolic efficiency based on residual feed intake in bulls, energy corrected milk in cows and intramuscular fat content in both genders. Tissue samples from jejunum, liver, skeletal muscle and rumen were subjected to global transcriptomic analysis via stranded total RNA sequencing (RNAseq) and blood plasma samples were used for profiling of 640 metabolites. To identify lncRNAs within the indicated tissues, a project-specific transcriptome annotation was established. Subsequently, novel transcripts were categorized for potential lncRNA status, yielding a total of 7,646 predicted lncRNA transcripts belonging to 3,287 loci. A regulatory impact factor approach highlighted 92, 55, 35, and 73 lncRNAs in jejunum, liver, muscle, and rumen, respectively. Their ensuing high regulatory impact factor scores indicated a potential regulatory key function in a gene set comprising loci displaying differential expression, tissue specificity and loci overlapping with quantitative trait locus regions for residual feed intake or milk production. These were subjected to a partial correlation and information theory analysis with the prioritized gene set. Results and Conclusions: Independent, significant and group-specific correlations (|r| > 0.8) were used to build a network for the high and the low metabolic efficiency group resulting in 1,522 and 1,732 nodes, respectively. Eight lncRNAs displayed a particularly high connectivity (>100 nodes). Metabolites and genes from the partial correlation and information theory networks, which each correlated significantly with the respective lncRNA, were included in an enrichment analysis indicating distinct affected pathways for the eight lncRNAs. LncRNAs associated with metabolic efficiency were classified to be functionally involved in hepatic amino acid metabolism and protein synthesis and in calcium signaling and neuronal nitric oxide synthase signaling in skeletal muscle cells.

Heat stress directly impairs gut integrity and recruits distinct immune cell populations into the bovine intestine.

High ambient temperature has multiple potential effects on the organism such as hyperthermia, endotoxemia, and/or systemic inflammation. However, it is often difficult to discriminate between cause and consequence of phenotypic effects, such as the indirect influence of heat stress via reduced food intake. Lactating dairy cows are a particularly sensitive model to examine the effects of heat stress due to their intensive metabolic heat production and small surface:volume ratio. Results from this model show heat stress directly induced a so-far unknown infiltration of yet uncategorized cells into the mucosa and submucosa of the jejunum. Due to a pair-feeding design, we can exclude this effect being a consequence of the concurrent heat-induced reduction in feed intake. Isolation and characterization of the infiltrating cells using laser capture microdissection and RNA sequencing indicated a myeloic origin and macrophage-like phenotype. Furthermore, targeted transcriptome analyses provided evidence of activated immune- and phagocytosis-related pathways with LPS and cytokines as upstream regulators directly associated with heat stress. Finally, we obtained indication that heat stress may directly alter jejunal tight junction proteins suggesting an impaired intestinal barrier. The penetration of toxic and bacterial compounds during heat stress may have triggered a modulated immune repertoire and induced an antioxidative defense mechanism to maintain homeostasis between commensal bacteria and the jejunal immune system. Our bovine model indicates direct effects of heat stress on the jejunum of mammals already at moderately elevated ambient temperature. These results need to be considered when developing concepts to combat the negative consequences of heat stress.

Development and evaluation of a milk protein transcript depletion method for differential transcriptome analysis in mammary gland tissue.

BACKGROUND: In the mammary gland transcriptome of lactating dairy cows genes encoding milk proteins are highly abundant, which can impair the detection of lowly expressed transcripts and can bias the outcome in global transcriptome analyses. Therefore, the aim of this study was to develop and evaluate a method to deplete extremely highly expressed transcripts in mRNA from lactating mammary gland tissue. RESULTS: Selective RNA depletion was performed by hybridization of antisense oligonucleotides targeting genes encoding the caseins (CSN1S1, CSN1S2, CSN2 and CSN3) and whey proteins (LALBA and PAEP) within total RNA followed by RNase H-mediated elimination of the respective transcripts. The effect of the RNA depletion procedure was monitored by RNA sequencing analysis comparing depleted and non-depleted RNA samples from Escherichia coli (E. coli) challenged and non-challenged udder tissue of lactating cows in a proof of principle experiment. Using RNase H-mediated RNA depletion, the ratio of highly abundant milk protein gene transcripts was reduced in all depleted samples by an average of more than 50% compared to the non-depleted samples. Furthermore, the sensitivity for discovering transcripts with marginal expression levels and transcripts not yet annotated was improved. Finally, the sensitivity to detect significantly differentially expressed transcripts between non-challenged and challenged udder tissue was increased without leading to an inadvertent bias in the pathogen challenge-associated biological signaling pathway patterns. CONCLUSIONS: The implementation of selective RNase H-mediated RNA depletion of milk protein gene transcripts from the mammary gland transcriptome of lactating cows will be highly beneficial to establish comprehensive transcript catalogues of the tissue that better reflects its transcriptome complexity.

Analysis of pig transcriptomes suggests a global regulation mechanism enabling temporary bursts of circular RNAs.

To investigate the dynamics of circRNA expression in pig testes, we designed specific strategies to individually study circRNA production from intron lariats and circRNAs originating from back-splicing of two exons. By applying these methods on seven Total-RNA-seq datasets sampled during the testicular puberty, we detected 126 introns in 114 genes able to produce circRNAs and 5,236 exonic circRNAs produced by 2,516 genes. Comparing our RNA-seq datasets to datasets from the literature (embryonic cortex and postnatal muscle stages) revealed highly abundant intronic and exonic circRNAs in one sample each in pubertal testis and embryonic cortex, respectively. This abundance was due to higher production of circRNA by the same genes in comparison to other testis samples, rather than to the recruitment of new genes. No global relationship between circRNA and mRNA production was found. We propose ExoCirc-9244 (SMARCA5) as a marker of a particular stage in testis, which is characterized by a very low plasma estradiol level and a high abundance of circRNA in testis. We hypothesize that the abundance of testicular circRNA is associated with an abrupt switch of the cellular process to overcome a particular challenge that may have arisen in the early stages of steroid production. We also hypothesize that, in certain circumstances, isoforms and circular transcripts from different genes share functions and that a global regulation of circRNA production is established. Our data indicate that this massive production of circRNAs is much more related to the structure of the genes generating circRNAs than to their function. Abbreviations: PE: Paired Ends; CR: chimeric Read; SR: Split Read; circRNA: circular RNA; NC: non conventional; ExoCirc-RNA: exonic circular RNA; IntroLCirc-: name of a porcine intronic lariat circRNA; ExoCirc-: name of a porcine exonic circRNA; IntronCircle-: name of a porcine intron circle; sisRNA: stable intronic sequence RNA; P: porcine breed Pietrain; LW: porcine breed Large White; RT: reverse transcription/reverse transcriptase; Total-RNA-seq: RNA-seq obtained from total RNA after ribosomal depletion; mRNA-seq: RNA-seq of poly(A) transcripts; TPM: transcripts per million; CR-PM: chimeric reads per million; RBP: RNA binding protein; miRNA: micro RNA; E2: estradiol; DHT: dihydrotestesterone.

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.

Long noncoding RNAs are associated with metabolic and cellular processes in the jejunum mucosa of pre-weaning calves in response to different diets.

Long noncoding RNAs (lncRNAs) emerged as important regulatory component of mechanisms involved in gene expression, chromatin modification and epigenetic processes, but they are rarely annotated in the bovine genome. Our study monitored the jejunum transcriptome of German Holstein calves fed two different milk diets using transcriptome sequencing (RNA-seq). To identify potential lncRNAs within the pool of unknown transcripts, four bioinformatic lncRNA prediction tools were applied. The intersection of the alignment-free lncRNA prediction tools (CNCI, PLEK and FEELnc) predicted 1,812 lncRNA transcripts concordantly comprising a catalogue of 1,042 putative lncRNA loci expressed in the calves' intestinal mucosa. Nine lncRNA loci were differentially expressed (DE lncRNAs) between both calf groups. To elucidate their biological function, we applied a systems biology approach that combines weighted gene co-expression network analysis with functional enrichment and biological pathway analysis. Four DE lncRNAs were found to be strongly correlated with a gene network module (GNM) enriched for genes from canonical pathways of remodeling of epithelial adherens junction, tight junction and integrin signaling. Another DE lncRNA was strongly correlated with a GNM enriched for genes associated with energy metabolism and maintaining of cellular homeostasis with a focus on mitochondrial processes. Our data suggest that these DE lncRNAs may play potential regulatory roles in modulating biological processes associated with energy metabolism pathways and cellular signaling processes affecting the barrier function of intestinal epithelial cells of calves in response to different feeding regimens in the pre-weaning period.

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.

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.

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.

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.

The mammalian cervical vertebrae blueprint depends on the T (brachyury) gene.

A key common feature all but three known mammalian genera is the strict seven cervical vertebrae blueprint, suggesting the involvement of strong conserving selection forces during mammalian radiation. This is further supported by reports indicating that children with cervical ribs die before they reach reproductive age. Hypotheses were put up, associating cervical ribs (homeotic transformations) to embryonal cancer (e.g., neuroblastoma) or ascribing the constraint in cervical vertebral count to the development of the mammalian diaphragm. Here, we describe a spontaneous mutation c.196A > G in the Bos taurus T gene (also known as brachyury) associated with a cervical vertebral homeotic transformation that violates the fundamental mammalian cervical blueprint, but does not preclude reproduction of the affected individual. Genome-wide mapping, haplotype tracking within a large pedigree, resequencing of target genome regions, and bioinformatic analyses unambiguously confirmed the mutant c.196G allele as causal for this previously unknown defect termed vertebral and spinal dysplasia (VSD) by providing evidence for the mutation event. The nonsynonymous VSD mutation is located within the highly conserved T box of the T gene, which plays a fundamental role in eumetazoan body organization and vertebral development. To our knowledge, VSD is the first unequivocally approved spontaneous mutation decreasing cervical vertebrae number in a large mammal. The spontaneous VSD mutation in the bovine T gene is the first in vivo evidence for the hypothesis that the T protein is directly involved in the maintenance of the mammalian seven-cervical vertebra blueprint. It therefore furthers our knowledge of the T-protein function and early mammalian notochord development.