The Expression and Epigenetic Characteristics of the HSF2 Gene in Cattle-Yak and the Correlation with Its Male Sterility.

Aberrant expression of the heat shock proteins and factors was revealed to be closely associated with male reproduction. Heat shock factor 2 (HSF2) is a transcription factor that is involved in the regulation of diverse developmental pathways. However, the role and the corresponding molecular mechanism of HSF2 in male cattle-yak sterility are still poorly understood. Therefore, the aim of this study was to obtain the sequence and the biological information of the cattle-yak HSF2 gene and to investigate the spatiotemporal expression profiles of the locus during the development of cattle-yak testes. Additionally, the differential expression was analyzed between the cattle-yak and the yak, and the methylation of corresponding promoter regions was compared. Our results showed an additional 54 bp fragment and a missense mutation (lysine to glutamic acid) were presented in the cattle-yak HSF2 gene, which correlated with enriched expression in testicular tissue. In addition, the expression of the HSF2 gene showed dynamic changes during the growth of the testes, reaching a peak in adulthood. The IHC indicated that HSF2 protein was primarily located in spermatocytes (PS), spermatogonia (SP), and Sertoli cells (SC) in cattle-yak testes, compared with the corresponding cells of cattle and the yak. Furthermore, bisulfite-sequencing PCR (BSP) revealed that the methylated CpG sites in the promoter region of the cattle-yak HSF2 were more numerous than in the yak counterpart, which suggests hypermethylation of this region in the cattle-yak. Taken together, the low expression abundance and hypermethylation of HSF2 may underpin the obstruction of spermatogenesis, which leads to male cattle-yak infertility. Our study provided a basic guideline for the HSF2 gene in male reproduction and a new insight into the mechanisms of male cattle-yak sterility.

Molecular characterization and expression profile of the ALDH1A1 gene and its functions in yak luteal cells.

The ALDH1A1 gene encodes a cytoplasmic member of the aldehyde dehydrogenase 1 family, which plays an important role in regulating animal reproductive performance, including estrus cycle and embryonic development. The aim of this study was to characterize ALDH1A1 activity in ovaries of 3-5 year-old yaks and to determine its effects on cell proliferation, apoptosis, and progesterone secretion in luteal cells (LCs). The coding sequence (CDS) of the ALDH1A1 gene was cloned by reverse transcription-PCR and immunohistochemical analysis was used to confirm localization of the ALDH1A1 protein in the ovary. To assess the activity of ALDH1A1 in regulating progesterone secretion, si-ALDH1A1 was transfected into LCs in vitro and progesterone levels in LC supernatants were measured by ELISA. The interference efficiency was assessed by real-time quantitative PCR (RT-qPCR) and immunofluorescence staining, and cell proliferation and apoptosis were evaluated by EdU and TUNEL staining, respectively. The cloned ALDH1A1 sequence contained 1462 bp, encoding 487 amino acids. Immunohistochemical analysis showed that ALDH1A1 protein expression, which was significantly higher in LCs, was mainly found in antral follicles and the corpus luteum (CL). The expression of ALDH1A1 mRNA in LCs was effectively inhibited by si-ALDH1A1transfection, and progesterone secretion was markedly decreased along with the significant down-regulation of progesterone pathway-related genes, STAR, CYP11A1, CYP19A1, CYP17A1, 3ß-HSD, and HSD17B1. Knockdown of ALDH1A1 mRNA expression decreased cell proliferation and increased apoptosis in LCs. The mRNA expression of the proliferation-related genes, PCNA, CCND1, CCNB1 and CDC25A, was significantly down-regulated, while expression of the apoptosis-promoting CASP3 gene was significantly increased. In summary, we characterized the yak ALDH1A1 gene and revealed that ALDH1A1 knockdown promoted apoptosis, repressed cell proliferation, and decreased progesterone secretion by yak LCs, potentially by regulating the mRNA expression of genes related to proliferation, apoptosis, and progesterone synthesis and secretion.

Pangenome and multi-tissue gene atlas provide new insights into the domestication and highland adaptation of yaks.

BACKGROUND: The genetic diversity of yak, a key domestic animal on the Qinghai-Tibetan Plateau (QTP), is a vital resource for domestication and breeding efforts. This study presents the first yak pangenome obtained through the de novo assembly of 16 yak genomes. RESULTS: We discovered 290 Mb of nonreference sequences and 504 new genes. Our pangenome-wide presence and absence variation (PAV) analysis revealed 5,120 PAV-related genes, highlighting a wide range of variety-specific genes and genes with varying frequencies across yak populations. Principal component analysis (PCA) based on binary gene PAV data classified yaks into three new groups: wild, domestic, and Jinchuan. Moreover, we proposed a 'two-haplotype genomic hybridization model' for understanding the hybridization patterns among breeds by integrating gene frequency, heterozygosity, and gene PAV data. A gene PAV-GWAS identified a novel gene (BosGru3G009179) that may be associated with the multirib trait in Jinchuan yaks. Furthermore, an integrated transcriptome and pangenome analysis highlighted the significant differences in the expression of core genes and the mutational burden of differentially expressed genes between yaks from high and low altitudes. Transcriptome analysis across multiple species revealed that yaks have the most unique differentially expressed mRNAs and lncRNAs (between high- and low-altitude regions), especially in the heart and lungs, when comparing high- and low-altitude adaptations. CONCLUSIONS: The yak pangenome offers a comprehensive resource and new insights for functional genomic studies, supporting future biological research and breeding strategies.

MicroRNA-542-3p targets Pten to inhibit the myoblasts proliferation but suppresses myogenic differentiation independent of targeted Pten.

BACKGROUND: Non-coding RNA is a key epigenetic regulation factor during skeletal muscle development and postnatal growth, and miR-542-3p was reported to be conserved and highly expressed in the skeletal muscle among different species. However, its exact functions in the proliferation of muscle stem cells and myogenesis remain to be determined. METHODS: Transfection of proliferative and differentiated C2C12 cells used miR-542-3p mimic and inhibitor. RT-qPCR, EdU staining, immunofluorescence staining, cell counting kit 8 (CCK-8), and Western blot were used to evaluate the proliferation and myogenic differentiation caused by miR-542-3p. The dual luciferase reporter analysis and rescued experiment of the target gene were used to reveal the molecular mechanism. RESULTS: The data shows overexpression of miR-542-3p downregulation of mRNA and protein levels of proliferation marker genes, reduction of EdU+ cells, and cellular vitality. Additionally, knocking it down promoted the aforementioned phenotypes. For differentiation, the miR-542-3p gain-of-function reduced both mRNA and protein levels of myogenic genes, including MYOG, MYOD1, et al. Furthermore, immunofluorescence staining immunized by MYHC antibody showed that the myotube number, fluorescence intensity, differentiation index, and myotube fusion index all decreased in the miR-542-3p mimic group, compared with the control group. Conversely, these phenotypes exhibited an increased trend in the miR-542-3p inhibitor group. Mechanistically, phosphatase and tensin homolog (Pten) was identified as the bona fide target gene of miR-542-3p by dual luciferase reporter gene assay, si-Pten combined with miR-542-3p inhibitor treatments totally rescued the promotion of proliferation by loss-function of miR-542-3p. CONCLUSIONS: This study indicates that miR-542-3p inhibits the proliferation and differentiation of myoblast and Pten is a dependent target gene of miR-542-3p in myoblast proliferation, but not in differentiation.

Differential metabolites screening in yak (Bos grunniens) seminal plasma after cryopreservation and the evaluation of the effect of galactose on post-thaw sperm motility.

Sperm survival and activity depend on the provision of energy and nutrients from seminal plasma (SP). This study aimed to investigate the variations of metabolites within SP before and after freezing and subsequently explore the potential regulatory mechanisms affecting yak sperm cryodamage due to changes in metabolites in the SP. Untargeted metabolomics analysis was performed to screen for differential metabolites, followed by KEGG analysis to identify enriched signaling pathways. The combinatorial analysis of metabolomics and sperm proteomics revealed the influence of key SP metabolites on sperm proteins. Subsequently, the relevant differentially expressed proteins were verified by Western blot analysis. Finally, the mechanism underlying the positive effect of galactose on sperm motility was determined by assessing the change in ATP content in sperm before and after freezing and thawing. The data showed that a total of 425 and 269 metabolites were identified in the positive and negative ion modes, respectively. Freezing and thawing resulted in the up-regulation of 70 metabolites and the down-regulation of 29 metabolites in SP. The primary impact of freezing and thawing was observed in carbohydrate metabolism, including pyruvate metabolism, pentose phosphate pathway, galactose metabolism, the TCA cycle, and butanoate metabolism. In the combined analysis and Western blot results, a significant positive correlation was observed between galactose and Aldo-keto reductase family 1 member B1 (AKR1B1) (P < 0.05), which has the ability to convert galactose into galactol. Furthermore, the addition of galactose to thawed semen improved sperm motility by increasing AKR1B1 protein in sperm and was associated with the content of ATP. These data identify differential metabolites between fresh and frozen-thawed SP and suggest that galactose is a valuable additive for cryopreserved sperm, providing a theoretical basis for further exploration of the refrigerant formula for yak sperm cryopreservation.

Cloning and Characterization of Yak DHODH Gene and Its Functional Studies in a Bisphenol S-Induced Ferroptosis Model of Fetal Fibroblasts.

Dihydroorotate dehydrogenase (DHODH) is a rate-limiting enzyme of de novo biosynthesis of pyrimidine. Although the involvement of DHODH in resisting ferroptosis has been successively reported in recent years, which greatly advanced the understanding of the mechanism of programmed cell death (PCD), the genetic sequence of the yak DHODH gene and its roles in ferroptosis are still unknown. For this purpose, we firstly cloned the coding region sequence of DHODH (1188 bp) from yak liver and conducted a characterization analysis of its predictive protein that consists of 395 amino acids. We found that the coding region of the yak DHODH gene presented high conservation among species. Second, the expression profile of the DHODH gene in various yak tissues was investigated using RT-qPCR. The results demonstrated that DHODH was widely expressed in different yak tissues, with particularly high levels in the spleen, heart, and liver. Third, to investigate the involvement of DHODH in regulating ferroptosis in cells, yak skin fibroblasts (YSFs) were isolated from fetuses. And then, bisphenol S (BPS) was used to induce the in vitro ferroptosis model of YSFs. We observed that BPS decreased the cell viability (CCK8) and membrane potential (JC-1) of YSFs in a dose-dependent manner and induced oxidative stress by elevating reactive oxygen species (ROS). Simultaneously, it was evident that BPS effectively augmented the indicators associated with ferroptosis (MDA and BODIPY staining) and reduced GSH levels. Importantly, the co-administration of Ferrostatin-1 (Fer), a potent inhibitor of ferroptosis, significantly alleviated the aforementioned markers, thereby confirming the successful induction of ferroptosis in YSFs by BPS. Finally, overexpression plasmids and siRNAs of the yak DHODH gene were designed and transfected respectively into BPS-cultured YSFs to modulate DHODH expression. The findings revealed that DHODH overexpression alleviated the occurrence of BPS-induced ferroptosis, while interference of DHODH intensified the ferroptosis process in YSFs. In summary, we successfully cloned the coding region of the yak DHODH gene, demonstrating its remarkable conservation across species. Moreover, using BPS-induced ferroptosis in YSFs as the model, the study confirmed the role of the DHODH gene in resisting ferroptosis in yaks. These results offer valuable theoretical foundations for future investigations into the functionality of the yak DHODH gene and the underlying mechanisms of ferroptosis in this species.

Factors Influencing the Maturation and Developmental Competence of Yak (Bos grunniens) Oocytes In Vitro.

The yak (Bos grunniens) is a unique breed living on the Qinghai-Tibet Plateau and its surrounding areas, providing locals with a variety of vital means of living and production. However, the yak has poor sexual maturity and low fertility. High-quality mature oocytes are the basis of animal breeding technology. Recently, in vitro culturing of oocytes and embryo engineering technology have been applied to yak breeding. However, compared to those observed in vivo, the maturation rate and developmental capacity of in vitro oocytes are still low, which severely limits the application of in vitro fertilization and embryo production in yaks. This review summarizes the endogenous and exogenous factors affecting the in vitro maturation (IVM) and developmental ability of yak oocytes reported in recent years and provides a theoretical basis for obtaining high-quality oocytes for in vitro fertilization and embryo production in yaks.

Testis-specific knockout of Kdm2a reveals nonessential roles in male fertility but partially compromises spermatogenesis.

Lysine-specific histone demethylase 2 (Kdm2a) is a regulatory factor of histone modifications that participates in gametogenesis and embryonic development. The mis-regulation of Kdm2a can lead to aberrant gene expression, thereby contributing to abnormal cell proliferation, differentiation, apoptosis, and tumorigenesis. However, due to the potential confounding effects that are secondary to the loss of Kdm2a function from the soma in existing whole-animal mutants, the in vivo function of Kdm2a in spermatogenesis for male fertility remains unknown. Herein, we focus on exploring the spatiotemporal expression profile and biological functions of Kdm2a in the spermatogenesis and fertility of male mice. A testis-specific knockout Kdm2a model (Kdm2a cKO) was established by using the Stra8-Cre/loxP recombinase system to explore the roles of Kdm2a in male fertility. Our results showed that Kdm2a was ubiquitously expressed and dynamically distributed in multiple tissues and cell types in the testis of mice. Surprisingly, Kdm2a-deficient adult males were completely fertile and comparable with their control (Kdm2aflox/flox) counterparts. Despite the significantly reduced total number of sperm and density of seminiferous tubules in Kdm2a cKO testis accompanied by the degeneration of spermatogenesis, the fertilization ability and embryonic developmental competence of the Kdm2a cKO were comparable with those of their control littermates, suggesting that Kdm2a disruption did not markedly affect male fertility, at least during younger ages. Furthermore, Kdm2a homozygous mutants exhibited a lower total number and motility of sperm than the control group and showed notably affected serum 17ß-estradiol concentration. Interestingly, the transcriptome sequencing revealed that the loss of Kdm2a remarkably upregulated the expression level of Kdm2b. This effect, in turn, may induce compensative effects in the case of Kdm2a deficiency to maintain normal male reproduction. Together, our results reveal that Kdm2a shows spatiotemporal expression during testicular development and that its loss is insufficient to compromise the production of spermatozoa completely. The homologous Kdm2b gene might compensate for the loss of Kdm2a. Our work provides a novel Kdm2a cKO mouse allowing for the efficient deletion of Kdm2a in a testis-specific manner, and further investigated the biological function of Kdm2a and the compensatory effects of Kdm2b. Our study will advance our understanding of underlying mechanisms in spermatogenesis and male fertility.

Molecular cloning and characterization of STC1 gene and its functional analyses in yak (Bos grunniens) cumulus granulosa cells.

Cumulus granulosa cells (CGCs), an important type of ovarian somatic cells, carries out various functions related to oogenesis, follicular development, and steroidogenesis. Studying the biological mechanisms involved in the development and function of CGCs makes a great contribution to understanding the reproductive regulation in female animals. Stanniocalcin-1 (STC1) is an important Ca2+-regulated glycoprotein hormone that exhibits high expression levels in ovaries. In this study, we cloned the coding sequence of the yak STC1, predicted the structure of STC1 protein, detected the expression and localization of STC1 in yak ovaries, and analyzed the functions of STC1 in yak CGCs. The CDS (coding sequence) region of yak STC1 gene was found to be 744 bp and encoded 247 amino acids. Homology comparison revealed that STC1 protein was highly conserved among mammals. The STC1 mRNA displayed dynamic expression profiles in different stages of yak ovaries, and the highest expression was found in the follicular phase. Regarding localization, STC1 protein was widely distributed in various kinds of yak ovarian cells, including oocytes, mural granulosa cells, CGCs, and thecal cells. Repressing the expression of STC1 resulted in defective proliferation and survival of yak CGCs. In addition, knockdown the expression of STC1 repressed the secretion of progesterone and promoted the secretion of estrogen. Overexpression of STC1 partially rescued the proliferation of CGCs and resulted in opposite effects on the secretion of progesterone and estrogen. Several apoptosis and steroidogenesis-related genes, including BAX, BCL2, HSD3B1, HSD17B1, CYP11A1 and CYP17A1 showed altered expressions after repressing or increasing the expression of STC1 in yak CGCs. To the best of our knowledge, this study is the first to focus on the role of STC1 in yak CGCs, and the outcomes offer fresh insights into the mechanism governing yak reproduction.

Maternal Kdm2a-mediated PI3K/Akt signaling and E-cadherin stimulate the morula-to-blastocyst transition revealing crucial roles in early embryonic development.

Histone methylation plays an essential role in oocyte growth and preimplantation embryonic development. The modification relies on histone methyl-transferases and demethylases, and one of these, lysine-specific demethylase 2a (Kdm2a), is responsible for modulating histone methylation during oocyte and early embryonic development. The mechanism of how Kdm2a deficiency disrupts early embryonic development and fertility remains elusive. To determine if maternally deposited Kdm2a is required for preimplantation embryonic development, the expression profile of Kdm2a during early embryos was detected via immunofluorescence staining and RT-qPCR. The Kdm2a gene in oocytes was specifically deleted with the Zp3-Cre/LoxP system and the effects of maternal Kdm2a loss were studied through a comprehensive range of female reproductive parameters including fertilization, embryo development, and the number of births. RNA transcriptome sequencing was performed to determine differential mRNA expression, and the interaction between Kdm2a and the PI3K/Akt pathway was studied with a specific inhibitor and activator. Our results revealed that Kdm2a was continuously expressed in preimplantation embryos and loss of maternal Kdm2a suppressed the morula-to-blastocyst transition, which may have been responsible for female subfertility. After the deletion of Kdm2a, the global H3K36me2 methylation in mutant embryos was markedly increased, but the expression of E-cadherin decreased significantly in morula embryos compared to controls. Mechanistically, RNA-seq analysis revealed that deficiency of maternal Kdm2a altered the mRNA expression profile, especially in the PI3K/Akt signaling pathway. Interestingly, the addition of a PI3K/Akt inhibitor (LY294002) to the culture medium blocked embryo development at the stage of morula; however, the developmental block caused by maternal Kdm2a loss was partially rescued with a PI3K/Akt activator (SC79). In summary, our results indicate that loss of Kdm2a influences the transcriptome profile and disrupts the PI3K/Akt signaling pathway during the development of preimplantation embryo. This can result in embryo block at the morula stage and female subfertility, which suggests that maternal Kdm2a is a potential partial redundancy with other genes encoding enzymes in the dynamics of early embryonic development. Our results provide further insight into the role of histone modification, especially on Kdm2a, in preimplantation embryonic development in mice.

RFRP-3 Influences Apoptosis and Steroidogenesis of Yak Cumulus Cells and Compromises Oocyte Meiotic Maturation and Subsequent Developmental Competence.

RF amide-related peptide 3 (RFRP-3), a mammalian ortholog of gonadotropin-inhibitory hormone (GnIH), is identified to be a novel inhibitory endogenous neurohormonal peptide that regulates mammalian reproduction by binding with specific G protein-coupled receptors (GPRs) in various species. Herein, our objectives were to explore the biological functions of exogenous RFRP-3 on the apoptosis and steroidogenesis of yak cumulus cells (CCs) and the developmental potential of yak oocytes. The spatiotemporal expression pattern and localization of GnIH/RFRP-3 and its receptor GPR147 were determined in follicles and CCs. The effects of RFRP-3 on the proliferation and apoptosis of yak CCs were initially estimated by EdU assay and TUNEL staining. We confirmed that high-dose (10-6 mol/L) RFRP-3 suppressed viability and increased the apoptotic rates, implying that RFRP-3 could repress proliferation and induce apoptosis. Subsequently, the concentrations of E2 and P4 were significantly lower with 10-6 mol/L RFRP-3 treatment than that of the control counterparts, which indicated that the steroidogenesis of CCs was impaired after RFRP-3 treatment. Compared with the control group, 10-6 mol/L RFRP-3 treatment decreased the maturation of yak oocytes efficiently and subsequent developmental potential. We sought to explore the potential mechanism of RFRP-3-induced apoptosis and steroidogenesis, so we observed the levels of apoptotic regulatory factors and hormone synthesis-related factors in yak CCs after RFRP-3 treatment. Our results indicated that RFRP-3 dose-dependently elevated the expression of apoptosis markers (Caspase and Bax), whereas the expression levels of steroidogenesis-related factors (LHR, StAR, 3ß-HSD) were downregulated in a dose-dependent manner. However, all these effects were moderated by cotreatment with inhibitory RF9 of GPR147. These results demonstrated that RFRP-3 adjusted the expression of apoptotic and steroidogenic regulatory factors to induce apoptosis of CCs, probably through binding with its receptor GPR147, as well as compromised oocyte maturation and developmental potential. This research revealed the expression profiles of GnIH/RFRP-3 and GPR147 in yak CCs and supported a conserved inhibitory action on oocyte developmental competence.

TMT-based quantitative proteomics analysis reveals the differential proteins between fresh and frozen-thawed sperm of yak (Bos grunniens).

Sperm cryopreservation is one of the most effective methods for the conservation of germplasm resources and used of superior sires widely. However, the motility of yak (Bos grunniens) sperm was low after thawing and the proteomics changes in sperm cryopreservation remain unknown. Therefore, the aim of this study was to explore the differences between fresh sperm and frozen sperm of yak through the proteomic analysis and thus improve the understanding of sperm cryodamage. The Tandem Mass Tags (TMT) technology was used to screen differentially expressed proteins (DEPs) before and after freezing. Then, GO and KEGG analysis were conducted to analyze the DEPs enriched signaling pathways. Finally, the DEPs, including superoxide dismutase 1 (SOD1) and NADH ubiquinone oxidoreductase core subunit S8 (NDUFS8) were verified by the immunofluorescence technique. The results showed that there were 229 DEPs between fresh and frozen-thawed yak sperm. Compared with the fresh sperm, 120 proteins were up-regulated and 109 proteins were down-regulated in frozen-thawed sperm. The GO annotation showed that the up-regulated proteins enriched in metabolic and cytoskeleton-related processes, including lipoprotein metabolic process, lipid transport, extracellular region and intermediate filament cytoskeleton organization. In contrast, the down-regulated proteins enriched in biological processes including single fertilization, sperm capacitation and response to unfolded protein. KEGG pathway analysis indicated that freezing and thawing affected the oxidative phosphorylation pathway, the fructose and mannose metabolic pathway and the glycerolipid metabolic pathway of yak sperm. Immunofluorescence results showed that the protein expression level of SOD1 protein in the frozen group was significantly lower than that in the fresh group (P < 0.01), and the protein expression level of NDUFS8 protein was significantly higher in frozen group (P < 0.01). This study revealed the DEPs between fresh and frozen-thawed sperm and provides a theoretical basis to further explore the exertion of normal biological functions of yak sperm after freezing and thawing.

The Biological Characteristics and Differential Expression Patterns of TSSK1B Gene in Yak and Its Infertile Hybrid Offspring.

This study aimed to investigate the spatially and temporally expressed patterns and biological characteristics of TSSK1B in male yaks and explore the potential correlation between TSSK1B and male sterility of the yak hybrid offspring (termed cattle-yak). First, the coding sequence (CDS) of TSSK1B was cloned by RT-PCR, and bioinformatics analysis was conducted with relevant software. Quantitative real-time PCR (RT-qPCR) was employed to detect the expression profile of TSSK1B in various tissues of male adult yaks, the spatiotemporal expression of TSSK1B in different stages of yak testes, and the differential expression of TSSK1B between yak and cattle-yak testes. The cellular localization of TSSK1B was determined by immunohistochemistry (IHC). Furthermore, the methylation status of the TSSK1B promoter region was analyzed by bisulfite-sequencing PCR (BSP). The results showed that TSSK1B was 1235 bp long, including 1104 bp of the CDS region, which encoded 367 amino acids. It was a conserved gene sharing the highest homology with Bos mutus (99.67%). In addition, the bioinformatics analysis revealed that TSSK1B was an unstable hydrophilic protein mainly containing the alpha helix of 34.06% and a random coil of 44.41%, with a transmembrane structure of 29 amino acids long. The RT-qPCR results demonstrated that TSSK1B was specifically expressed in yak testes compared with that in other tissues and especially highly expressed in adult yak testes. On the contrary, TSSK1B was hardly expressed in the testis of adult cattle-yak. IHC confirmed that TSSK1B protein was more strongly expressed in the testes of adult yaks than in their fetal and juvenile counterparts. Interestingly, nearly no expression was observed in the testes of cattle-yak compared with the corresponding testes of yak. Bisulfite-sequencing PCR (BSP) revealed that the methylated CpG sites in the TSSK1B promoter region of cattle-yak was significantly higher than that in the yak. Taken together, this study revealed that TSSK1B was specifically expressed in yak testes and highly expressed upon sexual maturity. Moreover, the rare expression in cattle-yak may be related to the hypermethylation of the promoter region, thereby providing a basis for further studies on the regulatory mechanism of TSSK1B in male cattle-yak sterility.

Identification of microRNA transcriptome throughout the lifespan of yak (Bos grunniens) corpus luteum.

The corpus luteum (CL) is a temporary organ that plays a critical role for female fertility by maintaining the estrous cycle. MicroRNA (miRNA) is a class of non-coding RNAs involved in various biological processes. However, there exists limited knowledge of the role of miRNA in yak CL. In this study, we used high-throughput sequencing to study the transcriptome dynamics of miRNA in yak early (eCL), middle (mCL) and late-stage CL (lCL). A total of 6,730 miRNAs were identified, including 5,766 known and 964 novels miRNAs. Three miRNAs, including bta-miR-126-3p, bta-miR-143 and bta-miR-148a, exhibited the highest expressions in yak CLs of all the three stages. Most of the miRNAs were 20-24 nt in length and the peak was at 22 nt. Besides, most miRNAs with different lengths displayed significant uracil preference at the 5'-end. Furthermore, 1,067, 280 and 112 differentially expressed (DE) miRNAs were found in eCL vs. mCL, mCL vs. lCL, and eCL vs. lCL, respectively. Most of the DE miRNAs were down-regulated in the eCL vs. mCL and eCL vs. lCL groups, and up-regulated in the mCL vs. lCL group. A total of 18,904 target genes were identified, with 18,843 annotated. Pathway enrichment analysis of the DE miRNAs target genes illustrated that the most enriched cellular process in each group included pathways in cancer, PI3K-Akt pathway, endocytosis, and focal adhesion. A total of 20 putative target genes in 47 DE miRNAs were identified to be closely associated with the formation, function or regression of CL. Three DE miRNAs, including bta-miR-11972, novel-miR-619 and novel-miR-153, were proved to directly bind to the 3'-UTR of their predicated target mRNAs, including CDK4, HSD17B1 and MAP1LC3C, respectively. Both of these DE miRNAs and their target mRNAs exhibited dynamic expression profiles across the lifespan of yak CL. This study presents a general basis for understanding of the regulation of miRNA on yak CL and also provides a novel genetic resource for future analysis of the gene network during the estrous cycle in the yak.

The molecular characteristic analysis of TRIB2 gene and its expressional patterns in Bos grunniens tissue and granulosa cells.

Tribbles homolog 2 (TRIB2) plays an important role in the follicular development of female mammals. However, its expression and function in the yak (Bos grunniens) are still unclear. In this study, we predicted the molecular characteristics of TRIB2, and revealed its expression pattern in yak (Bos grunniens) tissues and ovarian granulosa cells. We cloned the full length of the yak TRIB2 gene obtained by RT-PCR was 1368 bp and the coding sequence (CDS) was 624 bp, encoding 207 amino acids (AA). Homology analysis showed that the yak TRIB2 is highly conserved among species. TRIB2 was detected to be extensively expressed in seven tissues of the yak liver, spleen, lung, kidney, ovary, oviduct and uterus by qPCR. The expression of TRIB2 mRNA in the ovary during gestation was significantly lower than that in the non-pregnant (p < 0.05). At each stage of follicle development, the TRIB2 mRNA in granulosa cells showed a significant upward trend with the development of follicles. The expression of TRIB2 gradually decreased with the increase of the culture time of the granulosa cells in vitro. In conclusion, these results suggest that TRIB2 may play an important role in the follicular development of yaks.

Analysis of PPP1R11 expression in granulosa cells during developmental follicles of yak and its effects on cell function.

The aims of this study were to analyse the protein phosphatase 1 regulatory subunit 11 (PPP1R11) expression and cellular localization in yak follicles and investigate its effects on cell proliferation, apoptosis and oestrogen secretion in granulosa cells (GCs). Ten healthy and non-pregnant female yaks (4-year-old) were used as experimental animals. The mRNA relative expression level of PPP1R11 in GCs from small (<3.0 mm), medium (3.0-5.9 mm) and large (6.0-9.0 mm) follicles was detected by RT-qPCR, and the cellular localization of PPP1R11 protein was detected by immunohistochemistry staining (IHC). After isolation, culture and identification of yak GCs in vitro, si-PPP1R11 and si-NC (negative control) were transfected into GCs. RT-qPCR and immunofluorescence staining were used to evaluate the interference efficiency, and ELISA was performed to detect oestrogen concentration. Then, EdU staining and TUNEL staining were conducted to analyse cell proliferation and apoptosis. In addition, the oestrogen synthesis, proliferation- and apoptosis-related genes were detected by RT-qPCR after knockdown PPP1R11. The results showed that PPP1R11 is mainly located in ovarian GCs, and the expression levels of PPP1R11 in GCs from large follicles were significantly higher than that from medium and small follicles. Transfection of si-PPP1R11 into GCs could significantly inhibit the expression of PPP1R11. Interestingly, the oestrogen secretion ability and the expression level of oestrogen pathway-related genes (STAR, CYP11A1, CYP19A1 and HSD17B1) were also significantly downregulated. Moreover, the proportion of positive cells was decreased, and cellular proliferation-related genes (PCNA, CCNB1 and CDC25A) were significantly downregulated after knockdown PPP1R11. However, the proportion of apoptotic cells was increased, and apoptosis-related genes (BAX, CASP3 and P53) were significantly upregulated. Taken together, this study was the first revealed the expression and cellular localization of PPP1R11 in yak follicles. Interference PPP1R11 could reduce oestrogen secretion, inhibit proliferation and promote apoptosis in GCs, which provided a basis for further studies on the regulatory mechanism of PPP1R11 in follicle development.

Oocyte-Specific Knockout of Histone Lysine Demethylase KDM2a Compromises Fertility by Blocking the Development of Follicles and Oocytes.

The methylation status of histones plays a crucial role in many cellular processes, including follicular and oocyte development. Lysine-specific demethylase 2a (KDM2a) has been reported to be closely associated with gametogenesis and reproductive performance, but the specific function and regulatory mechanism have been poorly characterized in vivo. We found KDM2a to be highly expressed in growing follicles and oocytes of mice in this study. To elucidate the physiological role of Kdm2a, the zona pellucida 3-Cre (Zp3-Cre)/LoxP system was used to generate an oocyte Kdm2a conditional knockout (Zp3-Cre; Kdm2aflox/flox, termed Kdm2a cKO) model. Our results showed that the number of pups was reduced by approximately 50% in adult Kdm2a cKO female mice mating with wildtype males than that of the control (Kdm2aflox/flox) group. To analyze the potential causes, the ovaries of Kdm2a cKO mice were subjected to histological examination, and results indicated an obvious difference in follicular development between Kdm2a cKO and control female mice and partial arrest at the primary antral follicle stage. The GVBD and matured rates of oocytes were also compromised after conditional knockout Kdm2a, and the morphological abnormal oocytes increased. Furthermore, the level of 17ß-estradiol of Kdm2a cKO mice was only 60% of that in the counterparts, and hormone sensitivity decreased as the total number of ovulated and matured oocytes decreased after superovulation. After deletion of Kdm2a, the patterns of H3K36me2/3 in GVBD-stage oocytes were remarkedly changed. Transcriptome sequencing showed that the mRNA expression profiles in Kdm2a cKO oocytes were significantly different, and numerous differentially expressed genes were involved in pathways regulating follicular and oocyte development. Taken together, these results indicated that the oocyte-specific knockout Kdm2a gene led to female subfertility, suggesting the crucial role of Kdm2a in epigenetic modification and follicular and oocyte development.

MicroRNA-34b-5p targets PPP1R11 to inhibit proliferation and promote apoptosis in cattleyak Sertoli cells by regulating specific signaling pathways.

Cattleyaks, a hybrid of (♂) and yak (♀), exhibit the marked productivity and adaptability of plateau, but suffer from male infertility. Small non-coding RNAs, especially miRNAs, play crucial roles in spermatogenesis and affect the growth of Sertoli cells (SCs). The objective of the present study was to explore the interaction between miR-34b-5p and protein phosphatase 1 regulatory inhibitor subunit 11 (PPP1R11) and its effect on cattleyak SCs. RT-qPCR was used to determine the expression pattern of miR-34b-5p and PPP1R11, while the cellular and subcellular localization of PPP1R11 was determined by immunohistochemistry and immunocytochemistry. The interaction between MiR-34b-5p and PPP1R11 was evaluated by immunofluorescence, proliferation, apoptosis, and western blotting assays. The potential binding sites between miR-34b-5p and PPP1R11 were uncovered through targeted search of an online database, and verified using a dual luciferase reporter system. Our data show that miR-34b-5p is differentially expressed in the testes and SCs of cattleyaks compared to yaks. Overexpression of miR-34b-5p in SCs suppressed proliferation and induced apoptosis, while the effects of miR-34b-5p knockdown were the reverse. The 3'UTR of PPP1R11 was identified as a potential target site of miR-34b-5p, and this was validated by online database searches and our data from the dual-luciferase reporter assay, and it displayed an inverse expression pattern to miR-34b-5p in SCs. The effects of silencing PPP1R11 by siRNA were similar to the results of miR-34b-5p upregulation, but significantly different from miR-34b-5p downregulation in cattleyak SCs. The effects with PPP1R11 overexpression were opposite, suggesting a novel biofunctional role of PPP1R11 inactivation in depressing cattleyak SCs growth. Lastly, we confirmed that miR-34b-5p inhibited PPP1R11 expression and induced apoptosis by regulating proliferation- and apoptosis-related genes in SCs. Thus, miR-34b-5p regulates the apoptosis and proliferation of cattleyak SCs via targeting PPP1R11, which can provide an innovative direction for exploring the mechanism of cattleyak male sterility.

Molecular characterization of TRIB1 gene and its role in regulation of steroidogenesis in bos grunniens granulosa cells.

To explore the expression pattern of the TRIB1 gene in yak follicles and its effect on the steroidogenesis of granulosa cells (GCs). Here, 4-5 years old female yaks were treated as the subjects. Immunohistochemically assay found that TRIB1 protein was expressed in different developmental follicles. Among different cell types of follicles, including cumulus cells (CCs), granulosa cells (GCs) and theca cells (TCs), the TRIB1 protein was most abundant in GCs (P < 0.0001). In addition, we cloned the coding sequence (CDS) of the yak TRIB1 gene, which is 1119 bp, encoding 372 amino acids (AA). The amino acid sequence homology of TRIB1 is >80% to those of other species, except for zebrafish. To further explore the function of TRIB1 in steroidogenesis, the pcDNA3.1(+)-TRIB1 eukaryotic expression vector was constructed and then transfected into GCs. The data showed that overexpression of TRIB1 significantly reduced the progesterone (P4) secretion of granulosa cells measured by ELISA assay (P < 0.05), but not Estradiol (E2) secretion. Consistently, TRIB1 gain-of-function downregulated the mRNA levels of steroidogenesis related genes steroidogenic acute regulatory protein (StAR), cytochrome P450 family 11 subfamily A member 1 (CYP11A1) and 3ß-hydroxysteroid dehydrogenase (3ß-HSD) (P < 0.01), while cytochrome P450 family 17 subfamily A member 1 (CYP17A1) and cytochrome P450 family 19 subfamily A member 1 (CYP19A1) had no significant difference (P > 0.05). Interestingly, mito-tracker staining showed that mitochondrial number significantly decreased in TRIB1 overexpressed GCs (P < 0.01). Further, overexpression of TRIB1 inhibited mRNA levels of mitochondrial biogenesis related genes, including Mitochondrial transcription factor (TFAM) and Peroxisome proliferator-activated receptor alpha co-activator (PPARGC1A) (P < 0.05). Conclusively, this work indicates that TRIB1 inhibited progesterone synthesis of GCs might be involved in the reduction of the mitochondria number.

MicroRNA-342-3p regulates yak oocyte meiotic maturation by targeting DNA methyltransferase 1.

MicroRNAs (miRNAs) play vital roles in the development of oocytes and ovarian follicles. We have previously shown differential expression of miR-342-3p during yak oocyte maturation. In this study, we investigated the role of miR-342-3p in meiotic maturation of yak oocytes and the underlying mechanism. The profile of ovarian DNA methyltransferase 1 (DNMT1) expression was investigated in yak by RT-qPCR and western blot analyses. The pattern of Dnmt1 expression in various meiotic stages (GV stage, MI stage and MII stage) of yak oocyte maturation was then measured by immunofluorescence staining. The interaction between Dnmt1 and miR-342-3p was verified by dual-luciferase reporter assay. Finally, miR-342-3p inhibitors were microinjected into yak cumulus-oocyte complex to evaluate the effects on oocyte maturation. MiR-342-3p expression was upregulated in oocytes during meiotic maturation, with significantly higher levels in the MII stage compared with the GV- and MI stages (p < .05), whereas the opposite pattern of Dnmt1 expression was detected. In the period to sexual maturity (3-year-old), DNMT1 showed an age-related pattern of ovarian expression at both the gene and protein levels. Immunohistochemistry analysis also indicated maturation-stage-related differences in DNMT1 expression in the ovarian follicles and corpus luteum, with expression predominantly detected in cumulus cells and oocytes. MiR-342-3p inhibitors effectively upregulated Dnmt1 expression and significantly inhibited oocyte meiotic maturation. Taken together, our results indicate that miR-342-3p plays a vital role in the meiotic maturation of yak oocytes by targeting the 3'-untranslated regions (UTR) of Dnmt1 and provide a new perspective on the mechanism of this process.