Toward Understanding the Genetic Basis of Yak Ovary Reproduction: A Characterization and Comparative Analyses of Estrus Ovary Transcriptiome in Yak and Cattle

Hypoxia-Inducible Factor 1α Affects Yak Oocyte Maturation and Early Embryonic Development by Regulating Autophagy

In animal assisted reproductive technology, the production of high-quality oocytes is crucial. The yak, having lived in the Qinghai-Tibet Plateau for an extended period, has reproductive cells that are regulated by hypoxia-inducible factor 1α (HIF-1α). This study aimed to investigate the impact of HIF-1α on yak oocyte maturation and early embryonic development in vitro through the regulation of autophagy. The in vitro maturation process of yak oocytes involved the addition of the HIF-1α inducer DFOM and the inhibitor LW6 to examine their effects on yak oocyte maturation, early embryonic development, cell autophagy, cytochrome P450s (CYP450s) enzyme expression, and cumulus diffusion factors. The findings revealed that DFOM significantly upregulated the expression of HIF-1α, resulting in increased the cumulus diffusion area, elevated first polar body expulsion rate of oocytes, enhanced mitochondrial and actin levels, decreased ROS production, and reduced early apoptosis levels of oocytes. Moreover, DFOM promoted the expression of autophagy-related proteins, CYP450s enzymes, and cumulus diffusion factors, thereby enhancing oocyte maturation and early embryonic development. Conversely, LW6 exhibited opposite effects. The inhibition of autophagy levels with 3-MA during DFOM treatment yielded similar outcomes. Furthermore, reducing autophagy led to increased apoptosis levels at all stages of early embryonic development, as well as a significant decrease in total cell number and ICM/TE ratio of blastocysts. Studies have shown that during the in vitro maturation of yak oocytes, HIF-1α can affect the cumulus expansion area of oocytes by regulating autophagy, the first polar body excretion rate, mitochondrial level, actin level, ROS and early apoptosis level, the CYP450s enzyme, and the expression of cumulus expansion factors, thereby improving the in vitro maturation and early embryonic development of yak oocytes. These findings offer valuable insights into the reproductive regulation mechanism of yaks in hypoxic environments and suggest potential strategies for the advancement of yak assisted reproductive technology.

Long noncoding RNAs and mRNAs profiling in ovary during laying and broodiness in Taihe Black-Bone Silky Fowls (Gallus gallus Domesticus Brisson)

BACKGROUND

Broodiness significantly impacts poultry egg production, particularly notable in specific breeds such as the black-boned Silky, characterized by pronounced broodiness. An understanding of the alterations in ovarian signaling is essential for elucidating the mechanisms that influence broodiness. However, comparative research on the characteristics of long non-coding RNAs (lncRNAs) in the ovaries of broody chickens (BC) and high egg-laying chickens (GC) remains scant. In this investigation, we employed RNA sequencing to assess the ovarian transcriptomes, which include both lncRNAs and mRNAs, in eight Taihe Black-Bone Silky Fowls (TBsf), categorized into broody and high egg-laying groups. This study aims to provide a clearer understanding of the genetic underpinnings associated with broodiness and egg production.

RESULTS

We have identified a total of 16,444 mRNAs and 18,756 lncRNAs, of which 349 mRNAs and 651 lncRNAs exhibited significantly different expression (DE) between the BC and GC groups. Furthermore, we have identified the cis-regulated and trans-regulated target genes of differentially abundant lncRNA transcripts and have constructed an lncRNA-mRNA trans-regulated interaction network linked to ovarian follicle development. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation analyses have revealed that DE mRNAs and the target genes of DE lncRNAs are associated with pathways including neuroactive ligand-receptor interaction, CCR6 chemokine receptor binding, G-protein coupled receptor binding, cytokine-cytokine receptor interaction, and ECM-receptor interaction.

CONCLUSION

Our research presents a comprehensive compilation of lncRNAs and mRNAs linked to ovarian development. Additionally, it establishes a predictive interaction network involving differentially abundant lncRNAs and differentially expressed genes (DEGs) within TBsf. This significantly contributes to our understanding of the intricate interactions between lncRNAs and genes governing brooding behavior.

Exploring the induction of endometrial epithelial cell apoptosis in clinical-type endometritis in yaks through the cyt-c/caspase-3 signaling axis

Endometritis is a significant contributor to reduced productivity in yaks in Tibet, China. The Cyt-c/Caspase-3 signaling axis plays a crucial role in the mitochondrial pathway that triggers cell apoptosis due to endogenous factors. In this study, we examined the endometrial epithelial tissue of yaks with endometritis using pathological examination, immunohistochemical analysis, TUNEL staining, qRT-PCR, and Western blot. The results indicated significant changes in the apoptotic factors of the Cyt-c/Caspase-3 signaling axis. The expression levels of Bak1, Bax, Cyt-c, Apaf-1, Caspase-9, and Caspase-3 were significantly increased (P < 0.05), while the expression level of Bcl-2 was significantly decreased. Immunohistochemistry results revealed significant increase in Bak1, Bax, Cyt-c, Apaf-1, Caspase-9, and Caspase-3 expression in the cytoplasm compared to the healthy group, except for Bcl-2, which showed a significant decrease. Pathological section analysis demonstrated that clinical endometritis in yaks led to structural damage, bleeding, congestion, and inflammatory cell infiltration in the endometrial epithelium. Our study findings indicated that clinical endometritis in yaks can modulate apoptosis of endometrial epithelial cells via the Cyt-c/Caspase-3 signaling pathway, resulting in different levels of damage. This research is pioneering in exploring cell apoptosis induced by clinical endometritis in yaks, offering novel insights and potential strategies for the future prevention and treatment of endometritis in yaks.

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.

Differentially expressed transcripts study during pregnancy and postpartum anestrus of yak (Bos grunniens)

Background: Yak is the main livestock species in the plateau area, and its reproductive performance is low, usually two years or three years. A very few of yaks recover within a certain period of time after delivery and smoothly enter the next estrous cycle, while most of them enter the postpartum anestrus and show no estrus performance. However, the key biological factors and influencing mechanisms that cause postpartum anestrus in yaks are not clear. Objective: To study the expression of differential transcripts in ovaries of yak during pregnancy and postpartum anestrus. Methods: Each three yaks in pregnancy and anestrus under natural grazing conditions in Haiyan County, Qinghai Province were selected and slaughtered, and their ovaries were collected and sent to Biomarker Technologies. Oxford Nanopore Technologies single-molecule real-time electrical signal sequencing technology was used to perform full-length transcriptome sequencing. Astalavista software was used to identify the types of alternative splicing events in yak estrus and pregnancy, and TAPIS pipeline was used to identify alternative polyadenylation. Results: The results showed that there were 1751 differentially expressed transcripts (DETs) between pregnancy and anestrus in yak, of which 808 were upregulated and 943 were downregulated. GO analysis showed that the biological processes of DETs were mainly reproductive, reproductive and rhythmic processes. KEGG analysis showed that the DET cell junction-related adhesion junction protein (β-catenin) and amino terminal kinase (JNK) were involved in FAs (local adhesion). Phosphatidylinositol-3-kinase (PI3K) is involved in the PI3K/AKT/mTOR signaling pathway. Circadian rhythm output cycle failure (Clock) and brain and muscle tissue aromatic hydrocarbon receptor nuclear transporter-like protein 1 (Bmal1) are involved in circadian rhythm signaling pathway. Conclusion: This study found that β-catenin, JNK, PI3K, Clock and Bmal1 were closely related to postpartum anestrus in yak.

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.

Effects of Two Feeding Patterns on Growth Performance, Rumen Fermentation Parameters, and Bacterial Community Composition in Yak Calves

The health of young ruminants is highly dependent on early rumen microbial colonization. In this study, the effects of milk replacer on growth performance, rumen fermentation, and the rumen microflora in yak calves were evaluated. Sixty yak calves (body weight = 22.5 ± 0.95 kg, age = 30 ± 1 d) were assigned to the CON group (breastfed) or TRT group (milk replacer fed) and evaluated over 120 d. At 120 d, ruminal fluid samples were collected from 14 calves and then conducted for rumen fermentation and microbiota analyses. There was no difference in growth performance; however, calf survival was higher in the TRT group than in the CON group. The concentration of total volatile fatty acids and the molar proportion of butyric acid and lactic acid were increased with milk replacer feed in the TRT group (p < 0.05), but iso-valeric acid concentration was highest in the CON group (p < 0.05). Firmicutes and Bacteroidetes were the most dominant phyla in the CON and TRT groups, respectively. In the TRT group, Bacteroidetes, Prevotellaceae, Bacteroidia, Bacteroidetes, and Prevotella_1 were the dominant flora in the rumen of calves. The relative abundances of various taxa were correlated with rumen fermentation parameters; the relative abundance of Quinella and iso-butyrate levels were positively correlated (r = 0.57). The relative abundances of the Christensenellaceae_R-7_group and A/P were positively correlated (r = 0.57). In summary, milk replacer is conducive to the development of the rumen microflora, the establishment of rumen fermentation function, and the implementation of early weaning in yaks.

Sirtuin 7 is essential for the survival and synthesis of oestrogen in yak (Bos grunniens) cumulus granulosa cells

Cumulus granulosa cells (CGCs) are a type of important ovarian somatic cells that carries out various functions related to oogenesis, follicular development and embryogenesis. The study on the development and function of CGCs facilitates the understanding of reproductive regulation in female animals. Sirtuin 7 (SIRT7) is a member of the sirtuin family of NAD+-dependent deacetylases mediating numerous biological processes. In this study, we detected the localization of SIRT7 in yak ovaries as well as explored the function of SIRT7 in yak CGCs. The results revealed that the SIRT7 protein was mainly localized in the cytoplasm of oocytes, granulosa cells and theca cells. The knockdown of SIRT7 in yak CGCs repressed cell proliferation and impacted the expressions of several apoptosis-related genes. Furthermore, oestrogen synthesis was also inhibited in SIRT7-deficient yak CGCs. The expressions of several sterogenesis-related genes decreased significantly following SIRT7 knockdown. In addition, the lack of SIRT7 in yak CGCs resulted in decreased levels of the TGFB/SMAD family members TGFB1, TGFBR1 SMAD2 and SMAD3. Moreover, the activation of the TGFB/SMAD pathway by adding TGFB/SMAD pathway activator SRI-011381 partially rescued the level of oestrogen secreted by SIRT7-deficient yak CGCs, as well as the expressions of steroidogenesis-related genes NR5A1 and CYP19A1. This research is the first to focus on the role of SIRT7 in yak ovary, and the outcomes offer new insights into the mechanism governing yak reproduction.

Characterization of pig skeletal muscle transcriptomes in response to low temperature

OBJECTIVES

The response of mammals to cold environment is a complex physiological activity, and its underlying mechanism must be analyzed from multiple perspectives. Skeletal muscle is an important thermogenic tissue that maintains body temperature in mammals. We dissected the molecular mechanism of pig skeletal muscle response to a cold environment by performing comparative transcriptome analysis in the Enshi black pig.

METHODS

Three pigs were subjected to acute cold stress (3 days), three pigs were subjected to cold acclimation (58 days), and three pigs were used as controls. RNA-seq was used to screen the differentially expressed genes (DEGs) of skeletal muscle.

RESULTS

Using RNA-seq methods, we identified 1241 DEGs within the acute cold stress group and 1886 DEGs within the cold acclimation group. Prolonged cold exposure induced more gene expression changes. A total of 540 key cold-responsive DEGs were found, and their trends were consistent within the acute cold stress group and cold acclimation group. Gene expression pattern analysis showed that there were significant differences between the low-temperature treatment groups and the control group, and there were also differences between individuals after long-term low-temperature treatment. Analysis of DEGs revealed that 134 pathways were significantly enriched in the cold adaptation group, 98 pathways were significantly enriched in the acute cold stress group, and 71 pathways were shared between the two groups. The 71 shared pathways were mainly related to lipid, amino acid, and carbohydrate metabolism; signal transduction; endocrine, immune, and nervous system; cardiovascular disease; infectious diseases caused by bacteria or viruses; and neurodegenerative disease.

CONCLUSIONS

In conclusion, this study provides insights into the molecular mechanism of porcine skeletal muscle response under low-temperature environment. The data may assist further research on the mechanism of pig response to cold exposure.

Exosomes Derived from Yak Follicular Fluid Increase 2-Hydroxyestradiol Secretion by Activating Autophagy in Cumulus Cells

Exosomes in the follicular fluid can carry and transfer regulatory molecules to recipient cells, thus influencing their biological functions. However, the specific effects of yak follicular fluid exosomes on 2-hydroxyestrodiol (2-OHE₂) secretion remain unknown. Here, we investigated whether yak follicular fluid exosomes can increase 2-OHE₂ secretion through the activation of autophagy in cumulus cells (YCCs). In vitro cultured YCCs were treated with yak follicular fluid exosomes for 6, 12, and 24 h. The effects of yak follicular fluid exosomes on autophagy and 2-OHE₂ secretion were evaluated through real-time quantitative fluorescence PCR (RT-qPCR), Western blotting (WB), transfected with RFP-GFP-LC3, immunohistochemistry, and ELISA. To further investigate whether 2-OHE₂ secretion was related to autophagy, YCCs were administered with yak follicular fluid exosomes, 3-methyladenine (3-MA), and rapamycin (RAPA). The results revealed that treatment with yak follicular fluid exosomes activated autophagy in YCCs and increased 2-OHE₂ secretion. Conversely, the inhibition of autophagy with 3-MA blocked these effects, suggesting that autophagy has an important role in 2-OHE₂ secretion in YCCs. Treatment of YCCs with rapamycin showed similar results with yak follicular fluid exosomes as there was an increase in 2-OHE2 secretion due to the activation of autophagy in the treated cumulus cells. Our results demonstrate that autophagy is enhanced by yak follicular fluid exosomes, and this is associated with an increase in 2-OHE₂ secretion in YCCs.

Transcriptomic analysis of the thyroid and ovarian stroma reveals key pathways and potential candidate genes associated with egg production in ducks

The importance of thyroid-related genes has been repeatedly mentioned in the transcriptome studies of poultry with different laying performance, yet there are few systematic studies to unravel the regulatory mechanisms of the thyroid-ovary axis in the poultry egg production process. In this study, we compared the transcriptome profiles in the thyroid and ovarian stroma between high egg production (GP) and low egg production (DP) ducks, and then revealed the pathways and candidate genes involved in the process. We identified 1,114 and 733 differentially expressed genes (DEGs) in the thyroid and ovarian stroma, separately. The Gene Ontology (GO) analysis showed that a total of 504 and 189 GO terms were identified in the thyroid and ovarian stroma (P < 0.05). Three common GO terms were identified from the top 5 GO terms with the highest significant level in two tissues, including extracellular space, calcium ion binding, and integral component of plasma membrane. The enrichment analysis of the Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that 15 and 14 KEGG pathways were significantly (P < 0.05) enriched in the thyroid and ovarian stroma, respectively. And, there were 8 common pathways, including neuroactive ligand-receptor interaction, calcium signaling pathway, ECM-receptor interaction, PPAR signaling pathway, melanogenesis, wnt signaling pathway, vascular smooth muscle contraction, and cytokine-cytokine receptor interaction. Notably, the neuroactive ligand-receptor interaction pathway was the most significantly enriched by the DEGs both in the thyroid and ovarian stroma. The interaction among DEGs enriched in the neuroactive ligand-receptor interaction and ECM-receptor interaction suggested that the thyroid may regulate ovarian development by these genes. Through integrated analysis of the protein-protein interaction (PPI) network and KEGG pathway maps, 9 key DEGs (PTH, THBS2, THBS4, CD36, ADIPOQ, ACSL6, PRKAA2, CRH, and PCK1) were identified, which could play crucial roles in the thyroid to regulate ovarian function and then affect egg-laying performance between GP and DP. This study serves as a basis to explore the molecular mechanism of the thyroid affecting ovarian function and egg production in female ducks and may help to identify molecular markers that can be used for duck genetic selection.

House feeding system improves the estrus rate in yaks (Bos grunniens) by increasing specific fecal microbiota and myo-inositol content in serum

Grazing (G) yaks (Bos grunniens) are generally of low fertility, which severely limits the income of local pastoralists. However, we recently found that yaks had a 52% higher estrus rate in house feeding (HF) than in G. Gas chromatography-mass spectrometry (GC-MS) and 16S rRNA gene sequencing were used to analyze serum metabolites and fecal microbiota of 20 rutting yaks in the G and HF systems, respectively, to explain this phenomenon. The results showed that 73 total metabolites differed significantly (p &lt; 0.05 and VIP &gt; 1) between the G and HF systems. In the HF system, 53 were upregulated and 20 were downregulated compared with the G system. Organic oxygen compounds, organic acids and their derivatives, and lipids and lipid-like molecules were the most common differential metabolites. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway mapper revealed that 25 metabolic signaling pathways differed significantly between the two systems. The top three enriched pathways included central carbon metabolism in cancer, aminoacyl–tRNA biosynthesis, and ABC transporters. The 16S rRNA gene sequencing data showed no significant differences in Chao 1 index between the two systems. According to principal component analysis (PCA), the HF and G systems were distinctly and separately clustered in terms of fecal microbiota distribution. The G system showed significantly higher abundances of Firmicutes. The HF system showed significantly higher abundances of Alistipes, Treponema, and Rikenellaceae_ RC9_ gut_ group. Pearson's correlation analysis and core network analysis revealed that Rikenellaceae_RC9_ gut_ group, Alistipes, and Treponema were positively correlated with myo-inositol and formed the core bacteria. In summary, the HF system promoted the estrus rate and changed the composition of yak fecal microbiota and serum metabolites. Increased estrus rate might be obtained due to enhanced myo-inositol content in yak serum via the HF system. Correlation analysis suggested that myo-inositol content might also be partly increased via yak-specific fecal microbiota, contributing to the estrus rate. These findings could lead to a novel therapeutic strategy for G yaks due to their low estrus rate.

First Isolation and Characteristics of Bovine Parainfluenza Virus Type 3 from Yaks

The yaks belong to the genus Bos within the family Bovidae that live in the Tibet Plateau and is an indispensable economic resource for the local herders. Respiratory tract infections are common diseases in yaks caused by various pathogens; however, there have been no reports of bovine parainfluenza virus type 3 (BPIV3) infection. This study was conducted to investigate the pathogens and analyze their characteristics from the four yak lung samples with severe respiratory tract infection symptoms in the yak farm. Results showed that out of four lung samples, three were identified as BPIV3-positive by RT-PCR. A BPIV3 strain (10^6.5 TCID₅₀/mL) was successfully isolated from the BPIV3-positive lung samples using Madin–Darby bovine kidney cells. The isolate caused systemic infection in the BALB/c mice and induced pathological changes in the lungs. Moreover, three complete BPIV3 genomes were amplified from the clinical samples. Phylogenetic trees based on the complete genomes, hemagglutinin-neuraminidase protein (HN), phosphoprotein (P), and large polymerase subunit protein (L) amino acid sequences showed that the complete BPIV3 genomes belonged to BPIV3 genotype C, and clustered into a large branch with the Chinese strains, although the three yak BPIV3 strains were clustered into a small branch. Compared to known BPIV3 genotype C strains in GenBank, the three genomes of yak BPIV3 showed four identical amino acid mutations in the HN, P and L proteins, suggesting a unique genetic evolution of BPIV3 in yaks. This study first isolated and characterized the BPIV3 from yaks, which contributed to the understanding of the infection and evolution of BPIV3 in yaks in the Tibet Plateau.

CircRNAs in Xiang pig ovaries among diestrus and estrus stages

BACKGROUND

The fecundity of sows is a trait of major economic in pig industry. The molecular regulation of estrus cycles can affect the fecundity of female animals. Compared with the other pig breeds, Xiang pig exhibits the special estrus behaviors. CircRNAs are thought to involve in regulation of multiple biological processes. However, the potential roles of circRNAs in ovary regulation on Xiang pig estrus are largely unknown.

RESULTS

8,937 circRNAs were identified from eight libraries constructed from the ovarian samples of Xiang pig at estrus and diestrus stages by RNA sequencing method. Of which, 1,995 were high confidence circRNAs detected at least two junction reads in each ovary sample and seven circRNAs were validated by RT-PCR method. Furthermore, we identified 290 upregulated and 15 downregulated circRNAs in estrus ovaries. These differentially expressed circRNAs (DECs) derived from 273 host genes. And 207 miRNAs were identified to be targets sponged by 156 DECs with 432 binding sites, containing more than one miRNA binding site in each circRNA. Function enrichment analysis revealed that the host genes and the targets of miRNAs sponged by DECs were enriched in several reproduction-related signaling pathways, such as ovarian steroidogenesis, oocyte maturation, circadian rhythm, estrogen signaling pathway, GnRH signaling pathway, circadian entrainment, and oocyte meiosis. The circRNA-miRNA-mRNA networks revealed that 153 miRNAs interacting with 122 DECs and 86 miRNAs interacting with 84 DECs were involved in ovarian functions and ovarian circadian entrainment and circadian rhythm respectively. The DEC-miRNA-DEG (differentially expressed gene, DEG) networks associated with reproduction-related signaling pathways contained 22 DECs,18 miRNAs and 7 DEGs. 22 DECs were recognized as hub circRNAs during the estrus phase of Xiang pigs.

CONCLUSIONS

The circRNAs that function as miRNA sponges could play a key role in post-transcriptional regulation of gene expression during Xiang pig's estrus cycle.

Feminization of channel catfish with 17β-oestradiol involves methylation and expression of a specific set of genes independent of the sex determination region

Exogenous oestrogen 17β-oestradiol (E2) has been shown to effectively induce feminization in teleosts. However, the molecular mechanisms underlying the process remain unclear. Here, we determined global DNA methylation and gene expression profiles of channel catfish (Ictalurus punctatus) during early sex differentiation after E2 treatment. Overall, the levels of global DNA methylation after E2 treatment were not significantly different from those of controls. However, a specific set of genes were differentially methylated, which included many sex differentiation-related pathways, such as MARK signalling, adrenergic signalling, Wnt signalling, GnRH signalling, ErbB signalling, and ECM-receptor interactions. Many genes involved in these pathways were also differentially expressed after E2 treatment. Specifically, E2 treatments resulted in upregulation of female-related genes and downregulation of male-related genes in genetic males during sex reversal. However, E2-induced sex reversal did not cause sex-specific changes in methylation profiles or gene expression within the sex determination region (SDR) on chromosome 4, suggesting that E2-induced sex reversal was a downstream process independent of the sex determination process that was regulated by sex-specific methylation within the SDR.

Ovarian Transcriptomic Analysis of Ninghai Indigenous Chickens at Different Egg-Laying Periods

Egg production is an essential indicator of poultry fertility. The ovary is a crucial organ involved in egg production; however, little is known about the key genes and signaling pathways involved in the whole egg-laying cycle of hens. In order to explore the mechanism of egg production at different stages of the egg-laying process, ovarian tissues from four chickens were randomly selected for transcriptome analysis at each of the three ages (145 d, 204 d, and 300 d in the early, peak, and late stages of egg laying). A total of 12 gene libraries were constructed, and a total of 8433 differential genes were identified from NH145d vs. NH204d, NH145d vs. NH300d and NH300d vs. NH204d (Ninghai 145-day-old, Ninghai 204-day-old, and Ninghai 300-day-old), with 1176, 1653 and 1868 up-regulated genes, and 621, 1955 and 1160 down-regulated genes, respectively. In each of the two comparison groups, 73, 1004, and 1030 differentially expressed genes were found to be co-expressed. We analyzed the differentially expressed genes and predicted nine genes involved in egg production regulation, including LRP8, BMP6, ZP4, COL4A1, VCAN, INHBA, LOX, PTX3, and IHH, as well as several essential egg production pathways, such as regulation adhesion molecules (CAMs), calcium signaling pathways, neuroactive ligand–receptor interaction, and cytokine–cytokine receptor interaction. Transcriptional analysis of the chicken ovary during different phases of egg-lay will provide a useful molecular basis for study of the development of the egg-laying ovary.

Changes in Transcriptomic Profiles in Different Reproductive Periods in Yaks

Simple Summary

The molecular regulation mechanism of yak ovarian activity has attracted extensive attention. This study investigated the global gene expression profiles in different reproductive stages (anestrus, estrus, and pregnancy) by RNA-seq technology. Enrichment analysis revealed that DEGs were involved in the process of follicular growth, ovulation, and hormone metabolism. This study explored the regulation mechanism of the yak ovary in the reproductive cycle and laid a theoretical foundation for further understanding the reproductive characteristics of yak.

Abstract

Yak reproductive characteristics have received extensive attention, though the molecular regulation mechanism of its ovarian activity remains to be explored. Therefore, this study initially conducted a comparative analysis of yak ovarian activities in anestrus, estrus, and pregnancy regarding their morphology and histology, followed by implementing RNA sequencing (RNA-seq) technology to detect the overall gene expression and biological mechanism in different reproductive stages. H&E staining showed that there were more growing follicles and mature follicles in ovarian tissue sections during estrus than ovarian tissues during non-estrus. The RNA-seq analysis of yak ovary tissues in three periods showed that DEGs related to follicular development and hormone metabolism were screened in the three comparison groups, such as COL1A2, NR4A1, THBS2, PTGS2, SCARB1, STAR, and WNT2B. Bioinformatics analysis showed that these DEGs are involved in ion binding, cell development, metabolic processes, enriched in ECM–receptor interactions, steroid biosynthesis, together with aldosterone generation/discharge and Wnt/PI3K-Akt signaling pathways. In addition, we speculate alternate splice development events to have important role/s in regulating ovarian functional genomic expression profiles. These results provide essential knowledge aimed at scrutinizing pivotal biomarkers for yak ovarian activity, together with paving the way for enhancing researchers’ focus on improving yak reproductive performance.

Adaptation Mechanisms of Yak (Bos grunniens) to High-Altitude Environmental Stress

Living at a high altitude involves many environmental challenges. The combined effects of hypoxia and cold stress impose severe physiological challenges on endothermic animals. The yak is integral to the livelihood of the people occupying the vast, inhospitable Qinghai-Tibetan plateau and the surrounding mountainous region. Due to long-term selection, the yak exhibits stable and unique genetic characteristics which enable physiological, biochemical, and morphological adaptations to a high altitude. Thus, the yak is a representative model for mammalian plateau-adaptability studies. Understanding coping mechanisms provides unique insights into adaptive evolution, thus informing the breeding of domestic yaks. This review provides an overview of genetic adaptations in Bos grunniens to high-altitude environmental stress. Combined genomics and theoretical advances have informed the genetic basis of high-altitude adaptations.

Molecular cloning, sequence, and expression patterns of DNA damage induced transcript 3 (DDIT3) gene in female yaks (Bos grunniens)

Endoplasmic reticulum stress (ERS) plays an important role in regulating the reproductive process of female mammals, mainly involved in follicular atresia and corpus luteum regression. DNA damage induced transcript 3 (DDIT3) is a marker gene of ERS. The objectives of the present study were to clone and analyze the sequence and tissue expression characteristics of DDIT3 gene in female yaks. By reverse transcriptase-polymerase chain reaction (RT-PCR) strategy, we obtained full-length 507-bp DDIT3-cDNA, encoding for 168-aa protein. Yak DDIT3 exhibited highest and least identity with that of bison and horse, respectively. Real-time PCR analyses revealed that the expression level of DDIT3 gene in ovary was higher than that in heart, liver, kidney, spleen, lung, uterus and oviduct (p < 0.05). DDIT3 expression level in ovary and uterus during pregnancy was higher than that in follicular phase, luteal phase and fetus stage. DDIT3 was highly expressed in metaphase II oocytes and granulosa cells than that in germinal vesicle and metaphase I oocytes (p < 0.05), respectively. This is the first molecular characterization and expression patterns of DDIT3 gene in female yaks. These results indicated that the DDIT3 gene possibly plays an important role in regulating ovary function and pregnancy maintenance in yaks.

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.

Effects of Premating Calcium and Phosphorus Supplementation on Reproduction Efficiency of Grazing Yak Heifers

Simple Summary

This study was realized to explore the effects of calcium chloride (CaCl) and monocalcium phosphate (MCP) supplementation on the reproductive efficiency of grazing yak heifers. The body weight, serum markers of bone metabolism, and conception and calving rate of grazing yaks in control group and supplementary feeding groups were compared. The results revealed that supplementation with MCP but not CaCl could significantly improve the reproductive performance, possibly due to the improvement in body weight and bone phosphorus storage providing better estrous physiological conditions for grazing yak heifers. The findings of this study may be helpful and instructional to improve the reproductive efficiency of yaks on the Qinghai Tibet Plateau.

Abstract

Reproductive efficiency is the main factor limiting yak production on the Tibet Plateau. The purpose of this study was to investigate the effect of supplementation with calcium chloride (CaCl) and monocalcium phosphate (MCP) for 30 days before breeding on body weight (BW) change, serum bone metabolism biomarkers, conception rate, and calving rate of grazing yaks. Ninety 3 year old yak heifers (153.05 ± 6.56 kg BW) were assigned to three treatments (n = 30 per treatment): grazing without supplementation (CONT), grazing plus calcium chloride supplementation (CaCl), and grazing plus monocalcium phosphate supplementation (MCP). Compared with the CONT group, supplementation with CaCl increased the serum concentrations of osteocalcin and decreased the alkaline phosphatase (ALP) levels (p < 0.05); supplementation with MCP increased the average daily gain (ADG), serum concentrations of phosphorus (P) and osteocalcin, conception rate, and calving rate (p < 0.05), whereas it decreased the serum concentrations of hydroxyproline, ALP, and calcitonin (p < 0.05). Both CaCl and MCP supplementation had no effect on serum calcium (Ca) concentration. The ADG, conception rate, and calving rate were higher in the MCP group than in the CaCl group (p < 0.05), while the serum concentrations of hydroxyproline and calcitonin were lower (p < 0.05). It could be concluded that premating supplementation with MCP increased the body weight gain and subsequent conception and calving rate of grazing yaks. Supplementation with MCP had a positive effect on body condition and bone metabolism, thus providing a better estrous condition for grazing yak heifers, which could contribute to enhancing reproduction efficiency.

Ovarian transcriptomic analysis and follicular development of Leizhou black duck

This study investigated the factors that caused the differences in egg production during the development of ovarian follicles in Leizhou black ducks. Leizhou black ducks population was divided into 2 groups as high-yield group (HG) and low-yield group (LG). The number of eggs (NE), age at first egg (AFE), weight at first egg, and egg weight (EW) of both groups were recorded, and differences were analyzed using the t test. The logistic model was used to simulate the egg production curves to analyze the production rules. The ovarian follicles of both duck groups were collected to count the number of different grades sized follicles, weigh the ovaries, and observe follicular sections to analyze the developmental differences. Ovarian transcriptomic sequencing was performed to investigate differentially expressed genes and signal pathways in both duck groups. The results revealed a significant difference (P < 0.01) in the NE laid, AFE, and EW between both groups. Comparatively, HG had significantly more (P < 0.01) large yellow follicles (LYF) than LG. The density of medullary layer cells of the follicle section was greater in HG than LG ducks. Transcriptome sequencing revealed a total of 1,027 differentially expressed genes between the HG and LG ducks of which 495 genes were upregulated, and 532 genes were downregulated. Fifty genes were related to reproduction and reproductive processes. Kyoto Encyclopedia of Genes and Genomes–enriched signaling pathways revealed 274 signal pathways enriched in these differentially expressed genes of which the steroid biosynthesis pathway was significantly enriched. Analysis (Q < 0.05) showed that HSD3β → gonadotropin-releasing hormone (GnRH) and estrogen receptor (ESR) → LHβ/ERK1/2 were enriched in the steroid biosynthesis signal pathway. Follicle-stimulating hormone signal pathway mediated by HSD3β → GnRH and ESR → LHβ/ERK1/2 may be involved in ovarian follicle development to regulate LYF reserve process and affect its ovulation cycle, which in turn influence the egg production of Leizhou black ducks.

The expression of hypoxia-inducible factor-1 alpha in primary reproductive organs of the female yak (Bos grunniens) at different reproductive stages

The yak (Bos grunniens) is the most important livestock animal in high-altitude regions owing to its prominent adaptability to cold conditions, nutritional deficiencies and hypoxia. The reproductive organs exhibit different histological appearances and physiological processes at different reproductive stages. Hypoxia-inducible factor-1 alpha (HIF-1α) is the regulatory subunit of HIF-1 that crucially regulates the response to hypoxia in mammalian organisms. The goal of our study was to investigate the expression and distribution of HIF-1α in the primary yak reproductive organs at different reproductive stages. Samples of the ovary, oviduct and uterus of 15 adult female yaks were collected and used in the experiment. The expression and localization of HIF-1α proteins and mRNA were investigated using quantitative real-time polymerase chain reaction (qRT-PCR), Western blot (WB) and immunohistochemistry (IHC). The results indicated that the expression of HIF-1α protein in the ovary was higher during the luteal phase than during the follicular phase and gestation period (p < .05). In the oviduct, HIF-1α protein was also more highly expressed during the luteal phase than during the follicular phase and gestation period (p < .01). However, in the uterus, the HIF-1α protein had stronger expression during the gestation period than during the follicular phase (p < .01) and luteal phase (p < .05). The expression of HIF-1α mRNA was similar to that of its protein. Immunohistochemical analysis revealed intense immunostaining of HIF-1α proteins in the follicular granulosa cells, granular luteal cells, villous epithelial cells of the oviduct, endometrial glandular epithelium and luminal epithelium, foetal villous trophoblast, and epithelia of caruncular crypts. This study showed that the expression of HIF-1α in the ovary, oviduct and uterus varies according to the stage of the reproductive cycle. This implies that HIF-1α plays an important role in regulating the stage-specific physiological function of yak reproductive organs under hypoxic environments.

The Expression of ERK1/2 in Female Yak (Bos grunniens) Reproductive Organs

The main reproductive organs undergo different histological appearances and physiological processes under different reproductive statuses. The variation of these organs depends on a delicate regulation of cell proliferation, differentiation, and apoptosis. Extracellular signal-regulated kinases1/2 (ERK1/2) are members of the mitogen-activated protein kinase (MAPK) super family. They have important roles in regulating various biological processes of different cells, tissues, and organ types. Activated ERK1/2 generally promotes cell survival, but under certain conditions, ERK1/2 also have the function of inducing apoptosis. It is widely believed that ERK1/2 play a significant role in regulating the reproductive processes of mammals. The goal of our research is to investigate the expression and distribution of ERK1/2 in the yak's main reproductive organs during different stages. In the present study, samples of the ovary, oviduct, and uterus of 15 adult female yak were collected and used in the experiment. The ERK1/2 proteins, localization, and quantitative expression of their mRNA were investigated using immunohistochemistry (IHC), western blot (WB) and relative quantitative real-time polymerase chain reaction (RT-PCR). The results indicated that ERK1/2 proteins and their mRNA were highly expressed in the ovary of the luteal phase and gestation period, in the oviduct of the luteal phase, and in the uterus of the luteal phase and gestation period. Immunohistochemical analysis revealed a strong distribution of ERK1/2 proteins in follicular granulosa cells, granular luteal cells, villous epithelial cells of the oviduct, endometrial glandular epithelium, and luminal epithelium. These results demonstrated that the expression of ERK1 and ERK2 proteins and their mRNA in the yak's ovary, oviduct, and uterus varies with the stage of the reproductive cycle. The variation character of ERK1 and ERK 2 expression in the yak's main reproductive organs during different stages implies that they play an important role in regulating the reproductive function under different physiological statuses.

Transcriptome analysis revealed key signaling networks regulating ovarian activities in the domestic yak

Domestic yaks are the most important livestock species on the Qinghai-Tibetan Plateau. Adult female yaks normally breed in warm season and enter anestrous in cold season. Currently, how the ovarian activity is regulated at the molecular level remains to be determined. This study was conducted to investigate follicular development and gene expression patterns of yak ovarian tissues in the warm and cold seasons. Dynamics of follicular development was evaluated based on histological analyses and global gene expression was examined by using RNA-sequencing (RNA-seq) technology. Firstly, we found that follicle development of yak cows in cold season was different from that in warm season. Interestingly, ovaries collected from yaks in cold season contained a significant higher number of antral follicles and some of these follicles showed signs of polycystic structure, indicating abnormal granulosa cell function. RNA-seq analyses of ovarian tissues from non-pregnant adult yaks in cold and warm season revealed that a list of 320 transcripts were differentially expressed, specifically, 79 were up-regulated and 241 were down-regulated in the ovaries from yaks during the cold season. Further analysis demonstrated that transcripts associated with estrogen secretion and metabolism signaling pathway were altered, including FST, CYP1A1, PIK3R1 and PIK3R2. This study showed histological features of follicle development and revealed candidate genes that may have important roles in regulating ovarian activities in the yak seasonal reproduction.

Identification and characterization of microRNAs in the gonads of Crassostrea hongkongensis using high-throughput sequencing

Crassostrea hongkongensis is one of the three most-commonly cultivated oyster species in China. Although microRNAs (miRNAs) expression in the gonads have been widely investigated, few studies of miRNAs in mollusk gonads are available, particularly in oyster. In the present study, we analyzed the miRNAs expressed in the ovaries and testes of C. hongkongensis. We obtained 14,166,409 and 15,133,900 raw reads from the ovaries and testes, respectively, yielding 13,634,997 (ovarian) and 14,494,149 (testicular) 18-35-nt sequences. We mapped these sequences to the C. hongkongensis genome reference sequence, and identified 8,771,717 (ovarian) and 9,926,014 (testicular) sequences corresponding to miRNAs in the Rfam database. After blasting the miRNA sequences against the miRBase database, we identified 50 known mature miRNAs and 53 novel miRNAs. Of these, 27 miRNAs were significantly upregulated in ovaries as compared to the testes, and 43 miRNAs were significantly upregulated in the testes as compared to the ovaries. To validate the differential expression results generated by Illumina sequencing, we used RT-real-time quantitative PCR (RT-qPCR) to characterize the expression patterns of the six most differently expressed miRNAs (lgi-miR-1990, lgi-miR-1986, lgi-miR-263b, lgi-miR-279, lgi-miR-1992, and novel_98) as well as two miRNAs associated with gonad development (lgi-miR-29 and lgi-miR-8). Most of the RT-qPCR miRNA expression patterns were similar to those recovered by high-throughput sequencing with the exceptions of novel_98 and lgi-miR-1992. Gene Ontology (GO) annotations indicated that the multi-organism cellular process GO category was enriched with the target genes of the differentially expressed miRNAs. Additionally, the target genes were enriched in several KEGG pathways, including the ECM-receptor interaction, galactose metabolism, phagosome, and notch signaling pathway. These pathways are involved in gonadal differentiation and the maintenance of gonad function. This identification and characterization of the miRNAs differentially expressed between the ovaries and testes of C. hongkongensis will increase our understanding of the role of miRNAs in gonad differentiation in the oyster.

The genetic mechanism of high prolificacy in small tail han sheep by comparative proteomics of ovaries in the follicular and luteal stages

To investigate the genetic mechanism of sheep prolificacy, protein profiling of ovaries in the follicular and luteal phases was conducted. The tandem mass tag technique was used to analyze the proteomes of ovaries from STH sheep that did not have the FecB mutation in the bone morphogenetic protein receptor 1B gene. Parallel Reaction Monitoring (PRM) was operated to validate the target differentially abundant proteins (DAPs). The result showed, a total of 34,037 peptides were found, and 5074 proteins were identified. The screened DAPs strictly related to energy metabolism, hormone synthesis, ovarian function were significantly enriched in oxidative phosphorylation(COX7A, ND5, and UQCR10), ovarian steroidogenesis(StAR and HSD3B), taurine and hypotaurine metabolism(CSAD), glycosaminoglycan biosynthesis-heparin sulfate/heparin(GLCE), necroptosis(H2AX, AIFM1, and FTH1), protein digestion and absorption(COL4A1 and COL4A5) and glycosaminoglycan degradation(HYAL2 and HEXB) pathways. These analyses indicated that the reproductive performance of sheep is regulated through different pathways. In consequence, these findings are an important resource that can be used in future studies of the genetic mechanism of high fecundity traits in sheep, and these DAPs can be further investigated as candidate markers to predict prolificacy of sheep. SIGNIFICANCE: Litter size is an important quantitative trait, but the genetic mechanism of high-prolificacy is still unclear in sheep. Our study identified potential signaling pathways and differentially abundant proteins related to reproductive performance. These findings will facilitate a better revealing the mechanism and provide possible targets for molecular design breeding for the formation of polytocous traits in sheep.

Comparative transcriptomic analysis of high and low egg-producing duck ovaries

The egg-laying rate is an important indicator of egg production of laying ducks. Egg production directly impacts the economic benefits of the duck industry. In order to obtain better insight into the molecular mechanisms associated with the process of egg production, comparative transcriptomic analysis of the ovaries of Jinding ducks with high and low egg production was performed using the Illumina HiSeq 2500 system. A total of 843 differentially expressed genes (DEGs) was identified, 367 that were down-regulated and 476 that were up-regulated in high egg production (HEP) ovaries, as compared with low egg production (LEP) ovaries. Some genes, such as MC5R, APOD, ORAI1, and DYRK4, were more active in HEP ovaries, indicating that these genes may play important roles in regulation of egg production. Among these 843 DEGs, 685 were assigned to gene ontology (GO) categories. Of these, 25 genes were related to reproduction, and 30 were related to the reproductive process, including some associated with ovarian follicle development, circadian regulation of gene expression, circadian rhythm, and estrogen receptor binding. Furthermore, some important functional pathways were revealed, such as the steroid biosynthesis pathway, the endocrine and other factor-regulated calcium reabsorption pathways, circadian rhythm, the neuroactive ligand-receptor interaction pathway, fatty acid biosynthesis, and the calcium-signaling pathway, which appear to be much more active in the HEP group, as compared to those of the LEP group. The results of this study provide very useful information that may contribute to future functional studies of genes involved in bird reproduction.

Analysis of alternative splicing events by RNA sequencing in the ovaries of Xiang pig at estrous and diestrous

Alternative splicing (AS) is an important mechanism for controlling gene expression, which regulates multiple biological processes in higher organisms. Chinese indigenous Xiang pig has distinctive biological characteristics, such as small size, early sexual maturity, lower litter size and not very clear exhibition of estrous behaviors. To further understand how AS responds to estrous cycles in Xiang pig, the genome-wide analysis of AS events was performed by RNA-seq method in Xiang pig ovaries at diestrous and estrous. Using ASprofile program, we analyzed twelve basic AS events in Xiang pig ovaries and identified 68,775 AS events in 15,142 genes from diestrous ovaries and 69,493 AS events in 15,291 genes from the estrous ovaries with average 4.54 splicing events. 94.4-95.5% of expressed genes underwent alternative splicing in this tissues. The frequencies of AS events were similar to each other at diestrous and estrous. Transcription start site (TSS) was the predominant type of AS events, followed by transcription terminal site (TTS), and skipped exon (SKIP). The remaining type of AS events, e.g., intron retention (IR) and alternative exon ends (AE), showed the lower frequencies. Further comparison analysis of gene expression indicated that 4,433 genes had at least one splice variant differentially expressed during estrous, whereas only 2,382 of them were differentially expressed at gene level. Numerous genes involved in gonad development and hormone metabolism were differentially regulated through AS. Twelve genes with different types of alternatively splicing were validated by using RT-PCR method. The GO annotation and KEGG pathway analysis clearly revealed that a lot of DEGs (differentially expressed genes) and DSGs (differentially spliced genes) were involved in follicular development and ovarian steroid biosynthesis. A large number of DSGs, although not differentially expressed, were enriched in circadian rhythm and several signaling pathways. These pathways potentially regulated the female animal reproductive function in gene and/or AS level. Our results suggested that alternative splicing play an essential role in regulation of gene expression in female pigs during estrous. Numerous genes involved in gonad development, hormone metabolism, circadian rhythm were differentially regulated through alternative splicing.

Genetic Diversity, Molecular Phylogeny, and Selection Evidence of Jinchuan Yak Revealed by Whole-Genome Resequencing

Jinchuan yak, a newly discovered yak breed, not only possesses a large proportion of multi-ribs but also exhibits many good characteristics, such as high meat production, milk yield, and reproductive performance. However, there is limited information about its overall genetic structure, relationship with yaks in other areas, and possible origins and evolutionary processes. In this study, 7,693,689 high-quality single-nucleotide polymorphisms were identified by resequencing the genome of Jinchuan yak. Principal component and population genetic structure analyses showed that Jinchuan yak could be distinguished as an independent population among the domestic yak population. Linkage disequilibrium analysis showed that the decay rate of Jinchuan yak was the lowest of the domestic yak breeds, indicating that the degree of domestication and selection intensity of Jinchuan yak were higher than those of other yak breeds. Combined with archaeological data, we speculated that the origin of domestication of Jinchuan yak was ∼6000 yr ago (4000–10,000 yr ago). The quantitative dynamics of population growth history in Jinchuan yak was similar to that of other breeds of domestic and wild yaks, but was closer to that of the wild yak. No significant gene exchange between Jinchuan and other domestic yaks occurred. Compared with other domestic yaks, Jinchuan yak possessed 339 significantly and positively selected genes, several of which relate to physiological rhythm, histones, and the breed’s excellent production characteristics. Our results provide a basis for the discovery of the evolution, molecular origin, and unique traits of Jinchuan yak.