This research communication describes the profile of gene expression related to the synthesis of yak milk as determined via quantitative reverse transcription polymerase chain reaction (RT-qPCR). Significant up-regulation during lactation were observed in genes related to fatty acid (FA) uptake from blood (LPL, CD36), intracellular FA transport (FABP3), intracellular FA activation of long- and short-chain FAs (ACSS1, ACSS2, ACSL1), de novo synthesis (ACACA), desaturation (SCD), triacyglycerol (TAG) synthesis (AGPAT6, GPAM, LPIN1), lipid droplet formation (PLIN2, BTN1A1, XDH), ketone body utilisation (BDH1, OXCT1), and transcription regulation (THRSP, PPARGC1A). In particular, intracellular de novo FA synthesis (ACSS2, ACACA, and FABP3) and TAG synthesis (GPAM, AGPAT6, and LPIN1), whose regulation might be orchestrated as part of the gene network under the control of SERBF1 in the milk fat synthesis process, were more activated compared to levels in dairy cows. However, the genes involved in lipid droplet formation (PLIN2, XDH, and BTN1A1) were expressed at lower levels compared to those in dairy cows, where these genes are mainly controlled by the PPARG regulator.

The yak is one of the most important domestic animals in Tibetan life for providing basic resources such as milk, meat and transportation. Although yak milk production is not elevated, yak milk is superior to dairy cow milk in nutrient composition (protein and fat). However, the understanding of the metabolic mechanisms of yak mammary gland tissue during the lactation cycle remains elusive. In this study, GC-MS-based metabolomics was employed to study the metabolic variations in the yak mammary gland during the lactation cycle (pregnancy, lactation and dry period). Twenty-nine metabolites were up or downregulated during the lactation period. Compared to the dry period, during the lactation period the levels of oxalic acid were upregulated, while glycine and uridine were downregulated. Thirty-seven pathways were obtained when the 29 differential metabolites were imported into the KEGG pathway analysis. The most impacted pathways during the lactation cycle were glycine, serine and threonine metabolism; alanine, aspartate and glutamate metabolism; TCA cycle; glyoxylate and dicarboxylate metabolism; and pyrimidine metabolism. Our results provide important insights into the metabolic events involved in yak mammary gland development, lactogenesis and lactation, which can guide further research to improve milk yield and enhance the constituents of yak milk.

This research assessed the gene expression patterns related to the synthesis of milk in yak, which is characterized by high fat and protein content but low yield. The yak (Bos grunniens) is one of the most crucial domestic animals in Tibetan life; however, the genetic and molecular factors underlying yak milk protein synthesis remain understudied. Yak mammary biopsies harvested during late-pregnancy (d -15) through the end of subsequent lactation (d 1, 15, 30, 60, 180, and 240) were used to evaluate gene expression via real-time quantitative PCR. The expression pattern of 41 genes encompassing multiple pathways integral to milk protein synthesis including insulin, mammalian target of rapamycin (mTOR), 5' AMP-activated protein kinase, Jak2-Stat5 signaling, and the expression of glucose and AA transporters was evaluated. Our results confirmed that most upregulated genes increased from d -15 and peaked at d 30 or 60 and then remained relatively highly expressed. Specifically, there was an increased expression of mTOR-related amino acid transporters (SLC1A5, SLC7A5, and SLC36A1), glucose transporters (SLC2A1, SLC2A3, and SLC2A8), Jak2-Stat5 pathway (ELF5), and insulin signaling pathway components (IRS1, PDPK1, and AKT1). For activation of proteins synthesis, MTOR was significantly increased only at d 1. Among inhibitors of mTOR signaling, TSC1 and PRKAA2 were significantly upregulated during lactation. The RPL23 was downregulated among ribosomal components. In conclusion, a critical role for AA and glucose transporters and insulin signaling through mTOR for regulation of yak milk protein synthesis was revealed in this study of the yak mammary gland.

Background: The yak is a species of livestock which is crucial for local communities of the Qinghai-Tibet Plateau and adjacent regions and naturally owns one more thoracic vertebra than cattle. Recently, a sub-population of yak termed as the Jinchuan yak has been identified with over half its members own a thoracolumbar vertebral formula of T15L5 instead of the natural T14L5 arrangement. The novel T15L5 positioning is a preferred genetic trait leading to enhanced meat and milk production. Selective breeding of this trait would have great agricultural value and exploration of the molecular mechanisms underlying this trait would both accelerate this process and provide us insight into the development and regulation of somitogenesis. Results: Here we investigated the genetic background of the Jinchuan yak through resequencing fifteen individuals, comprising five T15L5 individuals and ten T14L5 individuals with an average sequencing depth of > 10X, whose thoracolumbar vertebral formulae were confirmed by anatomical observation. Principal component analysis, linkage disequilibrium analysis, phylogenetic analysis, and selective sweep analysis were carried out to explore Jinchuan yak’s genetic background. Three hundred and thirty candidate markers were identified as associated with the additional thoracic vertebrae and target sequencing was used to validate seven carefully selected markers in an additional 51 Jinchuan yaks. The accuracies of predicting 15 thoracic vertebrae and 20 thoracolumbar vertebrae with these 7 markers were 100.00 and 33.33% despite they both could only represent 20% of all possible genetic diversity. Two genes, PPP2R2B and TBLR1, were found to harbour the most candidate markers associated with the trait and likely contribute to the unique somitic number and identity according to their reported roles in the mechanism of somitogenesis. Conclusions: Our findings provide a clear depiction of the Jinchuan yak’s genetic background and a solid foundation for marker-assistant selection. Further exploitation of this unique population and trait could be promoted with the aid of our genomic resource.

The yak is primarily found throughout the Tibetan high plateau and the surrounding mountainous area of south central Asia; among its others attributes, its milk is very important for the local population. A key concern in the field of yak research is the better understanding of which genes control the production and composition of milk. The most accurate and sensitive method for gene expression analysis is quantitative reverse transcription polymerase chain reaction (RT-qPCR). It is essential for reliable RT-qPCR to be able to the normalize the data using internal control genes (ICGs). However, it is critical to assess the reliability of the normalization by testing multiple ICGs. Our objective was to uncover a reliable normalization for RT-qPCR data obtained from yak mammary tissue during the lactation cycle. We assessed the reliability of 10 ICGs (ACTB, EIF6, GAPDH, LRP10, MRPL39, MRPS15, MTG1, RPS8, RPS23, and UXT) using geNorm. The analysis revealed that all of the tested ICGs can be considered to be reliable, but the use of the 6 most stable ICGs should be applied to yield a reliable normalization factor (NF). We compared the results of 3 target genes (CSN1S1, ESR1, and MYC) normalized using 6, 3, or 1 of the best ICGs. We did not observe overall differences between the 3 normalization strategies with the exception of 1 time point in MYC. The use of only a single ICG is not recommended; thus, we concluded that the calculation of the NF using the 3 best ICGs, MRPS15, RPS23, and UXT, is a reliable normalization strategy for RT-qPCR data obtained from yak mammary tissue during pregnancy and lactation. A dilution effect of the ICGs due to a large increase in the mRNA of abundantly expressed genes in bovine and porcine mammary tissue during the lactation cycle was previously observed. To test for the presence of a dilution effect in our study, we evaluated the pattern of non-normalized RT-qPCR data of ICGs from pregnancy to lactation and compared them with the total RNA concentration, milk yield, and non-normalized RT-qPCR data of 3 target genes. With a few exceptions, the non- normalized RT-qPCR data for the tested ICGs was significantly increased by lactation and had a positive correlation with total RNA and the non-normalized RT-qPCR data of CSN1S1. These data clearly indicated the presence of a “concentration effect” of single mRNA that remains unexplained but needs to be accounted for during the normalization of RT-qPCR data. Based on our findings, we recommend that the NF of the MRPS15, RPS23, and UXT genes should be used in the normalization of RT-qPCR data obtained from mammary tissue of lactating yaks during pregnancy and lactation.

We used the newly-developed Dynamic Impact Approach (DIA) and gene network analysis to study the sow mammary transcriptome at 80, 100, and 110 days of pregnancy. A swine oligoarray with 13,290 inserts was used for transcriptome profiling. An ANOVA with false discovery rate (FDR < 0.15) correction resulted in 1,409 genes with a significant time effect across time comparisons. The DIA uncovered that Fatty acid biosynthesis, Interleukin-4 receptor binding, Galactose metabolism, and mTOR signaling were among the most-impacted pathways. IL-4 receptor binding, ABC transporters, cytokine-cytokine receptor interaction, and Jak-STAT signaling were markedly activated at 110 days compared with 80 and 100 days. Epigenetic and transcription factor regulatory mechanisms appear important in coordinating the final stages of mammary development during pregnancy. Network analysis revealed a crucial role for TP53, ARNT2, E2F4, and PPARG. The bioinformatics analyses revealed a number of pathways and functions that perform an irreplaceable role during late gestation to farrowing.