Identification and verification of differentially expressed genes in yak mammary tissue during the lactation cycle

Linkages between rumen microbiome, host, and environment in yaks, and their implications for understanding animal production and management

Livestock on the Qinghai-Tibetan Plateau is of great importance for the livelihood of the local inhabitants and the ecosystem of the plateau. The natural, harsh environment has shaped the adaptations of local livestock while providing them with requisite eco-services. Over time, unique genes and metabolic mechanisms (nitrogen and energy) have evolved which enabled the yaks to adapt morphologically and physiologically to the Qinghai-Tibetan Plateau. The rumen microbiota has also co-evolved with the host and contributed to the host's adaptation to the environment. Understanding the complex linkages between the rumen microbiota, the host, and the environment is essential to optimizing the rumen function to meet the growing demands for animal products while minimizing the environmental impact of ruminant production. However, little is known about the mechanisms of host-rumen microbiome-environment linkages and how they ultimately benefit the animal in adapting to the environment. In this review, we pieced together the yak's adaptation to the Qinghai-Tibetan Plateau ecosystem by summarizing the natural selection and nutritional features of yaks and integrating the key aspects of its rumen microbiome with the host metabolic efficiency and homeostasis. We found that this homeostasis results in higher feed digestibility, higher rumen microbial protein production, higher short-chain fatty acid (SCFA) concentrations, and lower methane emissions in yaks when compared with other low-altitude ruminants. The rumen microbiome forms a multi-synergistic relationship among the rumen microbiota services, their communities, genes, and enzymes. The rumen microbial proteins and SCFAs act as precursors that directly impact the milk composition or adipose accumulation, improving the milk or meat quality, resulting in a higher protein and fat content in yak milk and a higher percentage of protein and abundant fatty acids in yak meat when compared to dairy cow or cattle. The hierarchical interactions between the climate, forage, rumen microorganisms, and host genes have reshaped the animal's survival and performance. In this review, an integrating and interactive understanding of the host-rumen microbiome environment was established. The understanding of these concepts is valuable for agriculture and our environment. It also contributes to a better understanding of microbial ecology and evolution in anaerobic ecosystems and the host-environment linkages to improve animal production.

Dynamic Profile of the Yak Mammary Transcriptome during the Lactation Cycle

The objective of this study was to assess the transcriptome of the mammary tissue of four yaks during the whole lactation cycle. For this purpose, biopsies of the mammary gland were performed at -30, -15, 1, 15, 30, 60, 120, 180, and 240 days relative to parturition (d). The transcriptome analysis was performed using a commercial bovine microarray platform and the results were analyzed using several bioinformatic tools. The statistical analysis using an overall false discovery rate ≤ 0.05 for the effect of whole lactation and p < 0.05 for each comparison identified >6000 differentially expressed genes (DEGs) throughout lactation, with a large number of DEGs observed at the onset (1 d vs. -15 d) and at the end of lactation (240 d vs. 180 d). Bioinformatics analysis revealed a major role of genes associated with BTA3, BTA4, BTA6, BTA9, BTA14, and BTA28 in lactation. Functional analysis of DEG underlined an overall induction of lipid metabolism, suggesting an increase in triglycerides synthesis, likely regulated by PPAR signaling. The same analysis revealed an induction of amino acid metabolism and secretion of protein, with a concomitant decrease in proteasome, indicating a major role of amino acid handling and reduced protein degradation in the synthesis and secretion of milk proteins. Glycan biosynthesis was induced for both N-glycan and O-glycan, suggesting increased glycan content in the milk. The cell cycle and immune response, especially antigen processing and presentation, were strongly inhibited during lactation, suggesting that morphological changes are minimized during lactation, while the mammary gland prevents immune hyper-response. Transcripts associated with response to radiation and low oxygen were enriched in the down-regulated DEG affected by the stage of lactation. Except for this last finding, the functions affected by the transcriptomic adaptation to lactation in mammary tissue of yak are very similar to those observed in dairy cows.

Integrating genome-wide methylation and transcriptome-wide analyses to reveal the genetic mechanism of milk traits in Kazakh horses

Horse Milk has important quantitative characteristics and high economic value. However, the DNA methylation regulators involved in horse milk traits have not been clarified. To explore the important role of genome-wide DNA methylation in regulating equine milk yield, this study systematically investigated the genome-wide DNA methylation profiles of Kazakh horse blood by comparing a high-production group (HP, average daily milk yield of 7.5 kg) and low-production group (LP, average daily milk yield of 3.2 kg) using deep whole-genome bisulfite sequencing. First, both groups showed similar proportions of methylation at CpG sites. Subsequently, we identified 26,677 differential methylated regions (DMRs) of CG, 15 DMRs of CHG, 480 DMRs of CHH and 8268 DMR-related genes (DMGs). GO and KEGG analyses revealed that some DMGs were involved in regulating milk and milk component formation. By combining the WGBS-seq and the previous RNA-seq data, a total of 94 overlapping genes were obtained. Finally, we found that 9 DMGs are likely involved in milk production by Kazakh horses.

CEBPβ binding directly to the promoter region drives CEBPɑ transcription and improves FABP4 transcriptional activity in adipose tissue of yak (Bos grunniens)

Fatty acid binding protein 4 (FABP4) was crucial to fatty acid uptake and intracellular transport. However, the mechanisms regulating yak (Bos grunniens) FABP4 transcription were not determined. In the current study, predominant expression levels of yak FABP4 were identified in subcutaneous fat and longissimus dorsi muscles by quantitative real-time polymerase chain reactions (qPCR). The CCAAT/enhancer binding protein alpha (CEBPα) and myocyte enhancer factor 2A (MEF2A), as transcriptional activator or repressor in the promoter region of FABP4, were confirmed by both site-directed mutagenesis experiment and chromatin immunoprecipitation assay. Additionally, molecular mechanisms of CEBPɑ regulation were analyzed to explore the transcriptional regulatory property of FABP4, which indicated that transcriptional activity of CEBPɑ depended on CCAAT/ enhancer binding protein beta (CEBPβ) transcription factor. Our results demonstrated that CEBPβ binding directly to the promoter region drove CEBPɑ transcription, improving yak FABP4 transcriptional activity in adipocytes. This mechanism expanded the information on the transcriptional regulatory network of adipogenesis.

First DNA Sequencing in Beninese Indigenous Cattle Breeds Captures New Milk Protein Variants

This study investigated polymorphisms in the milk protein genes CSN1S1, CSN2, CSN1S2, CSN3, LALBA, and LGB, and casein haplotypes in Beninese indigenous cattle. Considering 67 animals, DNA sequencing of the genes' exons, flanking regions and parts of the 5'-upstream regions identified 1058 genetic variants including 731 previously unknown. In addition, four novel milk protein variants were detected, including CSN3K (p.Ala66Val), LALBAF (p.Arg58Trp), LGBB1 (p.Ala134Val) and LGBK (p.Thr92Asnfs*13). CSN3K is caused by a novel SNP (BTA6:85656526C>T, exon 4) whereas LALBAF and LGBB1 are due to rs714688595C>T (exon 1) and rs109625649C>T (exon 4), respectively. Regarding LGBK, a frameshift insertion of one adenine residue at BTA11:103257980 (exon 3) induces a premature translation termination resulting in a 46% reduction of the reference protein sequence. The casein polymorphisms formed five main CSN1S1-CSN2-CSN1S2-CSN3 haplotypes including B-A1-A-B, B-A1-A-A and C-A2-A-B which are predominant in the investigated cattle breeds. Moreover, in silico analyses of polymorphisms within the 5'- and 3'- untranslated regions of all six milk proteins revealed effects on microRNA and transcription factor binding sites. This study suggests a large genetic variation of milk protein genes in Beninese cattle, which should be investigated in further studies for their effects on milk production, including quality and yield traits.