The male infertility of cattleyak resulted from spermatogenic arrest has greatly restricted the effective utilization of the heterosis from crossbreeding of cattle and yak. Based on our previous studies, the significant divergences of the transcriptomic and proteomic sequencing between yak and cattleyak prompt us to investigate the critical roles of microRNAs in post-transcriptional regulation of gene expression during spermatogenesis. TUNEL-POD analysis presented sharply decreased spermatogenic cell types and the increased apoptotic spermatogonia in cattleyak. The STA-PUT velocity sedimentation was employed to obtain spermatogonia and spermatocytes from cattle, yak and cattleyak and these spermatogenic cells were verified by the morphological and phenotypic identification. MicroRNA microarray showed that 27 differentially expressed miRNAs were simultaneously identified both in cattleyak vs cattle and in cattleyak vs yak comparisons. Further analysis revealed that the down-regulation of bta-let-7 families, bta-miR-125 and bta-miR-23a might impair the RA-induced differentiation of spermatogonia. Target gene analysis for differentially expressed miRNAs revealed that miRNAs targeted major players involved in vesicle-mediated transport, regulation of protein kinase activity and Pathways in cancer. In addition, spermatogonia transfection analysis revealed that the down-regulation of bta-miR-449a in the cattleyak might block the transition of male germ cells from the mitotic cycle to the meiotic program. The present study provided valuable information for future elucidating the regulatory roles of miRNAs involved in spermatogenic arrest of cattleyak.

Milk fatty acids secreted by the mammary gland are one of the most important determinants of the nutritional value of goat milk. Unlike cow milk, limited data are available on the transcriptome-wide changes across stages of lactation in dairy goats. In this study, goat mammary gland tissue collected at peak lactation, cessation of milking, and involution were analyzed with digital gene expression (DGE) sequencing to generate longitudinal transcript profiles. A total of 51,299 unigenes were identified and further annotated to 12,763 genes, of which 9,131 were differentially expressed across various stages of lactation. Most abundant genes and differentially expressed genes (DEGs) were functionally classified through clusters of euKaryotic Orthologous Groups (KOG), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. A total of 16 possible expression patterns were uncovered, and 13 genes were deemed novel candidates for regulation of lactation in the goat: POLG, SPTA1, KLC, GIT2, COPS3, PDP, CD31, USP16/29/37, TLL1, NCAPH, ABI2, DNAJC4, and MAPK8IP3. In addition, PLA2, CPT1, PLD, GGA, SRPRB, and AP4S1 are proposed as novel and promising candidates regulating mammary fatty acid metabolism. “Butirosin and neomycin biosynthesis” and “Glyoxylate and dicarboxylate metabolism” were the most impacted pathways, and revealed novel metabolic alterations in lipid metabolism as lactation progressed. Overall, the present study provides new insights into the synthesis and metabolism of fatty acids and lipid species in the mammary gland along with more detailed information on molecular regulation of lactogenesis. The major findings will benefit efforts to further improve milk quality in dairy goats.

A prior study in bovine mammary (MACT) cells indicated that long-chain fatty acids (LCFA) C16:0 and C18:0, but not unsaturated LCFA, control transcription of milk fat-related genes partly via the activation of peroxisome proliferator-activated receptor gamma (PPARγ). However, in that study, the activation of PPARγ by LCFA was not demonstrated but only inferred. Prior data support a lower response of PPARγ to agonists in goat mammary cells compared to bovine mammary cells. The present study aimed to examine the hypothesis that LCFA alter the mRNA abundance of lipogenic genes in goat mammary epithelial cells (GMEC) at least in part via PPARγ. Triplicate cultures of GMEC were treated with a PPARγ agonist (rosiglitazone), a PPARγ inhibitor (GW9662), several LCFA (C16:0, C18:0, t10,c12-CLA, DHA, and EPA), or a combination of GW9662 with each LCFA. Transcription of 28 genes involved in milk fat synthesis was measured using RT-qPCR. The data indicated that a few measured genes were targets of PPARγ in GMEC (SCD1, FASN, and NR1H3) while more genes required a basal activation of PPARγ to be transcribed (e.g., LPIN1, FABP3, LPL, and PPARG). Among the tested LCFA, C16:0 had the strongest effect on upregulating transcription of measured genes followed by C18:0; however, for the latter most of the effect was via the activation of PPARγ. Unsaturated LCFA downregulated transcription of measured genes, with a lesser effect by t10,c12-CLA and a stronger effect by DHA and EPA; however, a basal activation of PPARγ was essential for the effect of t10,c12-CLA while the activation of PPARγ blocked the effect of DHA. The transcriptomic effect of EPA was independent from the activation of PPARγ. Data from the present study suggest that saturated LCFA, especially C18:0, can modulate milk fat synthesis partly via PPARγ in goats. The nutrigenomic effect of C16:0 is not via PPARγ but likely via unknown transcription factor(s) while PPARγ plays an indirect role on the nutrigenomic effect of polyunsaturated LCFA (PUFA) on milk fat related genes, particularly for CLA (permitting effect) and DHA (blocking effect).

The interactions between the rumen microbiota and the host are crucial for the digestive and absorptive processes of ruminants, and they are heavily influenced by the climatic conditions of their habitat. Owing to the harsh conditions of the high-altitude habitat, little is known about how ruminants regulate the host transcriptome and the composition of their rumen microbiota. Using the model species of goats, we examined the variations in the rumen microbiota, transcriptome regulation, and climate of the environment between high altitude (Lhasa, Xizang; 3650 m) and low altitude (Chengdu, Sichuan, China; 500 m) goats. The results of 16 S rRNA sequencing revealed variations in the abundance, diversity, and composition of rumen microbiota. Papillibacter, Quinella, and Saccharofermentans were chosen as potential microbes for the adaptation of Xizang goats to the harsh climate of the plateau by the Spearman correlation study of climate and microbiota. Based on rumen transcriptome sequencing analysis, 244 genes were found to be differentially expressed between Xizang goats and low-altitude goats, with 127 genes showing up-regulation and 117 genes showing down-regulation. SLC26A9, GPX3, ARRDC4, and COX1 were identified as potential candidates for plateau adaptation in Xizang goats. Moreover, the metabolism of fatty acids, arachidonic acids, pathway involving cytokines and their receptors could be essential for adaptation to plateau hypoxia and cold endurance. The expression of GPX3, a gene linked to plateau acclimatization in Xizang goats, was linked to the abundance of Anaerovibrio, and the expression of SLC26A9 was linked to the quantity of Selenomonas, according to ruminal microbiota and host Spearman correlation analysis. Our findings imply that in order to adapt harsh plateau conditions, Xizang goats have evolved to maximize digestion and absorption as well as to have a rumen microbiota suitable for the composition of their diet.