Differentially expressed genes and gene networks involved in pig ovarian follicular atresia

ECSIT: Biological function and involvement in diseases

Evolutionary conserved signaling intermediate in Toll pathways (ECSIT), a multi-functional protein, was first identified as a cytosolic adaptor protein in Toll-like receptors (TLRs) signaling-mediated innate immune responses. In the past two decades, studies have expanded the understanding of ECSIT. Nevertheless, there are still large knowledge gaps due to the inadequate number of studies regarding ECSIT, especially an overall review of ECSIT is lacking. Here, we first comprehensively summarize the biological functions of ECSIT with particular focus on innate immune responses and mitochondrial homeostasis. Cumulative studies have reinforced that ECSIT is involved in the regulation of innate immune responses through activating NF-κB signaling and potentiating the Retinoic acid-induced gene Ⅰ (RIG-Ⅰ)/ mitochondrial antiviral- signaling protein (MAVS) pathway-mediated innate antiviral immunity. In addition, ECSIT determines the mitochondrial morphology and function including mitochondrial complex Ⅰ (CⅠ) assembly, mitochondrial reactive oxygen species (mROS) production, mitochondrial membrane potential (MMP) maintenance and mitochondrial quality control. Owing to these distinct functions, ECSIT is involved in the etiology and pathology of human diseases including Alzheimer's disease (AD), cardiac hypertrophy, musculoskeletal disintegration, cancer, extranodal natural killer/T cell lymphoma (ENKTL) and ischemic stroke. Collectively, the roles and mechanisms of ECSIT under physiological and pathological conditions are critically discussed to provide a clearer view of the therapeutic potential of ECSIT.

Comparative transcriptome profile analysis of granulosa cells from porcine ovarian follicles during early atresia

At various stages of ovarian follicular development, more than 99% of follicles will be eliminated through a degenerative process called atresia. The regulatory mechanisms of atresia have been elucidated to some extent, involving hormones, growth factors, cytokines, and other factors. However, the stimuli initiating atresia in follicular granulosa cells remain unknown. In this study, we isolated the granulosa cells from porcine ovarian follicles (3-5 mm diameter) divided into healthy follicles (HFs) and early atretic follicles (EAFs). We applied high-throughput RNA sequencing to identify and compare differentially expressed genes (DEGs) between HFs and EAFs. A total of 31,694 genes were detected, of which 21,806 were co-expressed in six samples, and 243 genes (p < 0.05; FDR < 0.05) were differentially expressed (DEGs), including 123 downregulated and 120 upregulated in EAFs. GO analysis highlighted hormone metabolism, plasma membrane localization, and transporter activity. The pathway analysis indicated that 51 DEGs, involved in steroidogenesis, cell adhesion molecules, and TGF-beta signaling pathways, were highly related to atresia. Additionally, the interaction network of DEGs (p < 0.01; FDR < 0.05) using STRING highlighted LHR, ACACB, and CXCR4 as central nodes. In summary, this transcriptome analysis enriched our knowledge of the shifted mechanisms in granulosa cells during early atresia and provided novel perspectives into the atresia initiation.

Unraveling the mysteries of chicken proteomics: Insights into follicle development and reproduction

Chicken proteomics is a valuable method for comprehending the many mechanisms involved in follicle growth and reproduction in birds. This study offers a thorough summary of the latest progress in chicken proteomics research, specifically highlighting the knowledge obtained regarding follicle development and reproductive physiology. Proteomic studies have revealed essential proteins and pathways that play a role in follicle development, including those that control oocyte size, maturation, and ovulation. Proteomic investigations have provided insight into the molecular pathways that govern reproductive processes. By utilizing advanced proteomic technologies, including mass spectrometry and protein microarray analysis, we have been able to identify and measure many proteins in chicken follicles at their different developmental stages. The utilization of proteomic methods has enabled the identification of previously unknown biomarkers for reproductive efficiency that expedited the creation of innovative diagnostic instruments for monitoring reproductive health in chicken. Chicken proteomics not only offers insights into follicle growth and reproduction but also uncovers the effects of environmental influences on reproductive function. This provides new opportunities for exploring the molecular pathways that cause these effects. The integration of current data with upcoming proteomic technologies offers the potential for innovative strategies to enhance chicken reproduction.

In-depth analysis of transcriptomes in ovarian cortical follicles from children and adults reveals interfollicular heterogeneity

The ovarian cortical reserve of follicles is vital for fertility. Some medical treatments are toxic to follicles, leading to premature ovarian insufficiency. Ovarian tissue cryopreservation is an established method to preserve fertility in adults and even applied in prepuberty despite unproven efficacy. Here, we analyze transcriptomes of 120 cortical follicles from children and adults for detailed comparison. We discover heterogeneity with two main types of follicles in both age groups: one with expected oocyte-granulosa profiles and another with predicted role in signaling. Transcriptional changes during growth to the secondary stage are similar overall in children and adults, but variations related to extracellular matrix, theca cells, and miRNA profiles are found. Notably, cyclophosphamide dose correlates with interferon signaling in child follicles. Additionally, morphology alone is insufficient for follicle categorization suggesting a need for additional markers. Marker genes for early follicle activation are determined. These findings will help refine follicular classification and fertility preservation techniques across critical ages.

A matter of new life and cell death: programmed cell death in the mammalian ovary

BACKGROUND

The mammalian ovary is a unique organ that displays a distinctive feature of cyclic changes throughout the entire reproductive period. The estrous/menstrual cycles are associated with drastic functional and morphological rearrangements of ovarian tissue, including follicular development and degeneration, and the formation and subsequent atrophy of the corpus luteum. The flawless execution of these reiterative processes is impossible without the involvement of programmed cell death (PCD).

MAIN TEXT

PCD is crucial for efficient and careful clearance of excessive, depleted, or obsolete ovarian structures for ovarian cycling. Moreover, PCD facilitates selection of high-quality oocytes and formation of the ovarian reserve during embryonic and juvenile development. Disruption of PCD regulation can heavily impact the ovarian functions and is associated with various pathologies, from a moderate decrease in fertility to severe hormonal disturbance, complete loss of reproductive function, and tumorigenesis. This comprehensive review aims to provide updated information on the role of PCD in various processes occurring in normal and pathologic ovaries. Three major events of PCD in the ovary-progenitor germ cell depletion, follicular atresia, and corpus luteum degradation-are described, alongside the detailed information on molecular regulation of these processes, highlighting the contribution of apoptosis, autophagy, necroptosis, and ferroptosis. Ultimately, the current knowledge of PCD aberrations associated with pathologies, such as polycystic ovarian syndrome, premature ovarian insufficiency, and tumors of ovarian origin, is outlined.

CONCLUSION

PCD is an essential element in ovarian development, functions and pathologies. A thorough understanding of molecular mechanisms regulating PCD events is required for future advances in the diagnosis and management of various disorders of the ovary and the female reproductive system in general.

Differential age-related transcriptomic analysis of ovarian granulosa cells in Kazakh horses

Introduction

The Kazakh horse, renowned for its excellence as a breed, exhibits distinctive reproductive traits characterized by early maturity and seasonal estrus. While normal reproductive function is crucial for ensuring the breeding and expansion of the Kazakh horse population, a noteworthy decline in reproductive capabilities is observed after reaching 14 years of age.

Methods

In this study, ovarian granulosa cells (GCs) were meticulously collected from Kazakh horses aged 1, 2, 7, and above 15 years old (excluding 15 years old) for whole transcriptome sequencing.

Results

The analysis identified and selected differentially expressed mRNAs, lncRNAs, miRNAs, and circRNAs for each age group, followed by a thorough examination through GO enrichment analysis. The study uncovered significant variations in the expression profiles of mRNAs, lncRNAs, miRNAs, and circRNAs within GCs at different stages of maturity. Notably, eca-miR-486-3p and miR-486-y exhibited the highest degree of connectivity. Subsequent GO, KEGG, PPI, and ceRNA network analyses elucidated that the differentially expressed target genes actively participate in signaling pathways associated with cell proliferation, apoptosis, and hormonal regulation. These pathways include but are not limited to the MAPK signaling pathway, Hippo signaling pathway, Wnt signaling pathway, Calcium signaling pathway, Aldosterone synthesis and secretion, Cellular senescence, and NF-kappa B signaling pathway-essentially encompassing signal transduction pathways crucial to reproductive processes.

Discussion

This research significantly contributes to unraveling the molecular mechanisms governing follicular development in Kazakh horses. It establishes and preliminarily validates a differential regulatory network involving lncRNA-miRNA-mRNA, intricately associated with processes such as cell proliferation, differentiation, and apoptosis and integral to the developmental intricacies of stromal follicles. The findings of this study provide a solid theoretical foundation for delving deeper into the realm of reproductive aging in Kazakh mares, presenting itself as a pivotal regulatory pathway in the context of horse ovarian development.

A spatiotemporal molecular atlas of the ovulating mouse ovary

Ovulation is essential for reproductive success, yet the underlying cellular and molecular mechanisms are far from clear. Here, we applied high-resolution spatiotemporal transcriptomics to map out cell type- and ovulation stage-specific molecular programs as function of time during follicle maturation and ovulation in mice. Our analysis revealed dynamic molecular transitions within granulosa cell types that occur in tight coordination with mesenchymal cell proliferation. We identified molecular markers for the emerging cumulus cell fate during the preantral-to-antral transition. We describe transcriptional programs that respond rapidly to ovulation stimulation and those associated with follicle rupture, highlighting the prominent roles of apoptotic and metabolic pathways during the final stages of follicle maturation. We further report stage-specific oocyte-cumulus cell interactions and diverging molecular differentiation in follicles approaching ovulation. Collectively, this study provides insights into the cellular and molecular processes that regulate mouse ovarian follicle maturation and ovulation with important implications for advancing therapeutic strategies in reproductive medicine.

Comparative miRNA expression profile analysis of porcine ovarian follicles: new insights into the initiation mechanism of follicular atresia

Follicular atresia occurs in every stage of ovarian development, which is relevant to female fertility. In the past decade, increasing studies have confirmed that miRNAs, a class of short non-coding RNAs, play an important role in follicular atresia by post-transcription regulation of their target genes. However, the function of miRNAs on follicular atresia initiation is unknown. In the present study, high-throughput small RNA sequencing was performed to analyze differential miRNA expression profiles between healthy (HF) follicles and early atretic (EAF) follicles. A total of 237 conserved miRNA were detected, and the miR-143 is the highest expressed in follicles. Meanwhile, we also found wide sequence variations (isomiRs) in porcine ovarian miRNA, including in 5'un-translation region, core seed sequences and 3'untranslation region. Furthermore, we identified 22 differentially expressed miRNAs in EAF groups compared to HF group, of which 3 miRNAs were upregulated, as well as 19 miRNAs were downregulated, and then the RT-PCR was performed to validate these profiles. The target genes of these differentially expressed miRNAs were predicted by using miRwalk, miRDB, and Targetscan database, respectively. Moreover, the gene ontology and KEGG pathway enrichment established that the regulating functions and signaling pathways of these miRNAs contribute to follicular atresia initiation and cell fate. In conclusion, this study provides new insights into the changes of miRNAs in early atretic follicles to demonstrate their molecular regulation in ovarian follicular atretic initiation.

Deciphering Cellular Heterogeneity and Communication Patterns in Porcine Antral Follicles by Single-Cell RNA Sequencing

The antral follicle stage is a critical period in mammalian oocyte maturation, marked by complex interactions between oocyte development and neighboring granulosa cells. Understanding the heterogeneity and communication patterns of granulosa cells within antral follicles is crucial for deciphering their roles in follicle development and oocyte maturation. Here, we employed single-cell RNA-sequencing to explore the molecular and cellular characteristics of porcine antral follicles. Our analysis revealed distinct subpopulations within mural and cumulus granulosa cells, indicating diverse cellular states and functions within the follicles. Functional enrichment analysis unveiled the involvement of specific subpopulations in steroid biosynthesis, cumulus expansion, and cellular communication. Moreover, comparing mature and less mature follicles highlighted differences in cell distribution and functions, indicating developmental-specific variations. Our findings shed light on the intricate cellular heterogeneity and communication network within porcine antral follicles, providing valuable insights into the regulation of follicle development and oocyte maturation in pigs. These results hold promise for improving pig reproductive efficiency and advancing human reproductive medicine.

MiR-202-5p Regulates Geese Follicular Selection by Targeting BTBD10 to Regulate Granulosa Cell Proliferation and Apoptosis

The regulation of granulosa cells (GCs) proliferation and apoptosis is the key step in follicular selection which determines the egg production performance of poultry. miR-202-5p has been reported to be involved in regulating the proliferation and apoptosis of mammalian ovarian GCs. However, its role in regulating the proliferation and apoptosis of goose GCs is still unknown. In the present study, the GCs of pre-hierarchical follicles (phGCs, 8-10 mm) and those of hierarchical follicles (hGCs, F2-F4) were used to investigate the role of miR-202-5p in cell proliferation and apoptosis during follicle selection. In phGCs and hGCs cultured in vitro, miR-202-5p was found to negatively regulate cell proliferation and positively regulate cell apoptosis. The results of RNA-seq showed that BTB Domain Containing 10 (BTBD10) is predicted to be a key target gene for miR-202-5p to regulate the proliferation and apoptosis of GCs. Furthermore, it is confirmed that miR-202-5p can inhibit BTBD10 expression by targeting its 3'UTR region, and BTBD10 was revealed to promote the proliferation and inhibit the apoptosis of phGCs and hGCs. Additionally, co-transfection with BTBD10 effectively prevented miR-202-5p mimic-induced cell apoptosis and the inhibition of cell proliferation. Meanwhile, miR-202-5p also remarkably inhibited the expression of Phosphatidylinositol-4,5-Bisphosphate 3-Kinase Catalytic Subunit Beta (PIK3CB) and AKT Serine/Threonine Kinase 1 (AKT1), while it was significantly restored by BTBD10. Overall, miR-202-5p suppresses the proliferation and promotes the apoptosis of GCs through the downregulation of PIK3CB/AKT1 signaling by targeting BTBD10 during follicular selection. Our study provides a theoretical reference for understanding the molecular mechanism of goose follicular selection, as well as a candidate gene for molecular marker-assisted breeding to improve the geese' egg production performance.

Regulatory mechanism and research progress of ferroptosis in obstetrical and gynecological diseases

Ferroptosis is a novel type of regulated cell death driven by iron-dependent lipid peroxidation, which is distinguished from traditional types of programmed cell death, such as apoptosis, proptosis and necrosis et al. Impaired iron homeostasis, lipid peroxidation and antioxidants depletion are three hallmarks of ferroptosis. Over the past years, emerging studies support the notion that ferroptosis might be involved in the pathology of obstetrical and gynecological diseases, including preeclampsia (PE), endometriosis (EMs) and polycystic ovarian syndrome (PCOS). In the PE condition, the high sensitivity of trophoblasts towards ferroptosis has been found to potentially link to inflammation, suboptimal vascular remodeling and aberrant hemodynamics, which are three prominent pathophysiological features of PE. As for EMs, compromised ferroptosis of endometrial cells was associated with the formation ectopic lesions, whereas in the nearby lesions, the presence of ferroptosis was suggested to promote the progression of EMs, contributing to the relative clinical manifestations. Ferroptosis has been implicated a crucial role in the initiation of ovarian follicular atresia, which might help to manage ovulation in PCOS patients. Taken together, this review explored the basis of ferroptosis mechanisms and comprehensively summarized the latest discovery of roles of ferroptosis on PE, EMs and PCOS, gaining a deeper insight into the pathogenesis of these obstetrical and gynecological diseases and investigation of novel therapeutic interventions.

Exploration the role of INHBA in Hu sheep granulosa cells using RNA-Seq

Activin/inhibin is an important factor for the fecundity of Hu sheep, and it is involved in follicular development in ovaries. Inhibin subunit beta A (INHBA) participates in the synthesis of activin A and inhibin A. In this study, we also noted a positive correlation between INHBA level and the secretion of both activin A and inhibin A in culture medium. Nevertheless, both knockdown and overexpression of INHBA downregulated the expression of Inhibin Subunit Alpha (INHA). Based on RNA-Sequencing, we further examined the effect and molecular mechanism of INHBA knockdown in GCs on mRNA expression. A total of 1,687 differentially expressed genes (DEGs) were identified (Fold change ≥ 2; False-discovory-rates (FDR) ≤ 0.01), of which 602 genes were upregulated and 1,087 genes were downregulated in the INHBA interference group compared with the control groups. Gene Ontology (GO) enrichment indicated that these DEGs were mainly involved in the regulation of cell cycle, protein serine/threonine kinase activity, and actin cytoskeleton reorganization. Moreover, DEGs were significantly enriched in 40 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including P53, progesterone-mediated oocyte maturation, and PI3K-AKT signaling pathways. We also noted a positive correlation between INHBA level and many PI3K/Akt/mTOR pathway-related genes at the gene or/and protein expression. Overall, this study may contribute to a better understanding of the roles of INHBA on GCs of prolific sheep, as well as the molecular effect of low INHBA expression on GCs, clarifying some reproductive failures.

Integrative analyses of the mRNA expression profile reveal the involvement of STC1 in chicken folliculogenesis

Efficient ovarian follicle development, maturation, and ovulation are critical for egg production performance. Previous research has underscored the importance of messenger RNAs (mRNAs) in regulating development and folliculogenesis in chicken ovarians. However, the molecular mechanism is not fully understood, especially in the late period of the laying cycle. In the present study, ovarian tissues from 80-week-old Hy-Line Brown layers (three with high and three with low rates of egg laying) were collected for transcriptome sequencing. A total of 306 differentially expressed genes (DEGs) were identified in this study, at a false discovery rate (FDR)-corrected P-value < 0.05 and a log2|fold change| (log2|FC|) ≥1.5. Among these DEGs, stanniocalcin 1 (STC1) was mainly related to cellular processes, single-organism processes, biological regulation, metabolic processes, developmental processes, and reproductive processes. Then, we further investigated the regulation of STC1 during chicken follicle development and found that STC1 inhibited the proliferation and stimulated the apoptosis of follicular granulosa cells (GCs), and decreased the expression of progesterone (P4) and estradiol (E2). Collectively, these results suggest that STC1 plays an important role in chicken follicle development by decreasing GC proliferation and steroidogenesis and stimulating GC apoptosis. This study contributes to the understanding of the reproductive biology of laying hens in the late period of the laying cycle and further lays a foundation for the improvement of egg production in poultry breeding.

Association Analysis of PRKAA2 and MSMB Polymorphisms and Growth Traits of Xiangsu Hybrid Pigs

In this study, Xiangsu hybrid pig growth traits were evaluated via PRKAA2 and MSMB as candidate genes. Sanger sequencing revealed three mutation sites in PRKAA2, namely, g.42101G>T, g.60146A>T, and g.61455G>A, and all these sites were intronic mutations. Moreover, six mutation sites were identified in MSMB: intronic g.4374G>T, exonic g.4564T>C, exonic g.6378G>A, exonic g.6386C>T, intronic g.8643G>A, and intronic g.8857A>G. Association analysis revealed that g.42101G>T, g.60146A>T, g.61455G>A, g.4374G>T, g.4564T>C, g.6378G>A, g.6386C>T, g.8643G>A, and g.8857A>G showed different relationship patterns among body weight, body length, body height, chest circumference, abdominal circumference, tube circumference, and chest depth. Real-time polymerase chain reaction results revealed that the expression of PRKAA2 was highest in the longissimus dorsi muscle, followed by that in the heart, kidney, liver, lung, and spleen. The expression of MSMB was highest in the spleen, followed by that in the liver, kidney, lung, heart, and longissimus dorsi muscle. These results suggest that PRKAA2 and MSMB can be used in marker-assisted selection to improve growth related traits in Xiangsu hybrid pigs, providing new candidate genes for Pig molecular breeding.

Enhancement of recombinant human IL-24 (rhIL-24) protein production from site-specific integrated engineered CHO cells by sodium butyrate treatment

Interleukin-24 (IL-24) has specific inhibitory effects on the proliferation of various tumor cells with almost no toxicity to normal cells. The antitumor activity of recombinant human IL-24 protein produced in mammalian cells is much higher than that of bacteria, but its expression level is extremely low. Sodium butyrate (NaBu) was utilized as a media additive to increase protein expression in Chinese hamster ovary cells. The site-specific integrated engineered cells FCHO/IL-24 were treated with NaBu under different culture conditions (10% and 0.5% serum adherent culture, 0.5% serum suspension culture). First, 3 days of 1 mmol/L NaBu treatment significantly increased rhIL-24 expression level in FCHO/IL-24 cells by 119.94 ± 1.5% (**p < 0.01), 57.49 ± 2.4% (**p < 0.01), and 20.17 ± 3.03% (*p < 0.05) under the above culture conditions. Second, NaBu has a time- and dose-dependent inhibitory effect on FCHO/IL-24 proliferation and induces G0/G1 phase arrest. Under 10% and 0.5% serum adherent culture, G0/G1 phase cells were increased by 11.3 ± 0.5% (**p < 0.01) and 15.0 ± 2.6% (**p < 0.01), respectively. No induction of apoptosis was observed under a high dosage of NaBu treatment. These results suggest that NaBu increases rhIL-24 secretion via inhibiting cell cycle progression, thereby trapping cells in the highly productive G0/G1 phase. Finally, with increasing NaBu dose, glucose concentration increased (**p < 0.01) while lactic acid and ammonia concentrations reduced significantly (**p < 0.01) in 10% and 0.5% serum adherent culture supernatant. RNA-seq showed that NaBu treatment affected multiple tumor and immune-related pathways. In conclusion, NaBu treatment dramatically promoted rhIL-24 production in engineered FCHO/IL-24 cells by altering downstream pathways and inducing G0/G1 cell arrest with little effect on apoptosis.

A single-cell atlas of the cycling murine ovary

The estrous cycle is regulated by rhythmic endocrine interactions of the nervous and reproductive systems, which coordinate the hormonal and ovulatory functions of the ovary. Folliculogenesis and follicle progression require the orchestrated response of a variety of cell types to allow the maturation of the follicle and its sequela, ovulation, corpus luteum formation, and ovulatory wound repair. Little is known about the cell state dynamics of the ovary during the estrous cycle and the paracrine factors that help coordinate this process. Herein, we used single-cell RNA sequencing to evaluate the transcriptome of >34,000 cells of the adult mouse ovary and describe the transcriptional changes that occur across the normal estrous cycle and other reproductive states to build a comprehensive dynamic atlas of murine ovarian cell types and states.

The TGF-β/SMAD Signaling Pathway Prevents Follicular Atresia by Upregulating MORC2

In mammals, female fertility is determined by the outcome of follicular development (ovulation or atresia). The TGF-β/SMAD signaling pathway is an important regulator of this outcome. However, the molecular mechanism by which the TGF-β/SMAD signaling pathway regulates porcine follicular atresia has not been fully elucidated. Microrchidia family CW-type zinc finger 2 (MORC2) is anovel epigenetic regulatory protein widely expressed in plants, nematodes, and mammals. Our previous studies showed that MORC2 is a potential downstream target gene of the TGF-β/SMAD signaling pathway. However, the role of MORC2 in porcine follicular atresia is unknown. To investigate this, qRT-PCR, western blotting, and TdT-mediated dUTP nick-end labeling were performed. Additionally, the luciferase activity assay was conductedto confirm that the TGF-β/SMAD signaling pathway regulates MORC2. Our results demonstrate that MORC2 is animportant anti-apoptotic molecule that prevents porcine follicular atresia via a pathway involving mitochondrial apoptosis, not DNA repair. Notably, this studyrevealsthat the TGF-β/SMAD signaling pathway inhibits porcine granulosa cell apoptosis by up-regulating MORC2. The transcription factor SMAD4 regulated the expression of MORC2 by binding to its promoter. Our results will help to reveal the mechanism underlying porcine follicular atresia and improve the reproductive efficiency of sows.

Expression profile analysis of sheep ovary after superovulation and estrus synchronisation treatment

Superovulation is a widely used reproductive technique in livestock production, but the mechanism of sheep's superovulation is not yet clear. Here, a method of superovulation and estrus synchronisation was used to treat female Duolang sheep. After treatment, there were significant differences in serum FSH and LH levels and the number of dominant follicles between the two groups of sheep. We identified a total of 5021 differentially expressed genes (11, 13 and 15 days after treatment) and performed RT‐qPCR analysis to identify several mRNA expression levels. GO and KEGG enrichment analysis revealed that differentially expressed genes were involved in the regulation of signalling pathways of follicular development, cell cycle, material synthesis, energy metabolism, such as COL3A1, RPS8, ACTA2, RPL7 RPS6 and TNFAIP6 may play a key role in regulating the development of follicles. Our results show a comprehensive expression profile after superovulation and estrus synchronisation treatment. We provide the basis for further research on breeding techniques to improve the ovulation rate and birth rate of livestock.

Discovery of Differentially Expressed MicroRNAs in Porcine Ovaries With Smaller and Larger Litter Size

The number of live births in a litter is an important reproductive trait, and is one of the main indicators which reflect the production level and economic benefit of a pig farm. The ovary is an important reproductive organ of the sow, and it undergoes a series of biological processes during each estrous cycle. A complex transcriptional network containing coding and non-coding RNAs in the ovary closely regulates the reproductive capability of sows. However, the molecular regulation mechanisms affecting sow litter size are still unclear. We investigated the expression profiles of microRNAs (miRNAs) in porcine ovaries from sows with smaller than average litter sizes (SLS) and those with larger litter sizes (LLS). In total, 411 miRNAs were identified, and of these 17 were significantly down-regulated and 16 miRNAs were up-regulated when comparing sows with LLS and SLS, respectively. We further characterized the role of miR-183 which was one of the most up-regulated miRNAs. CCK-8, EdU incorporation and western blotting assays demonstrated that miR-183 promoted the proliferation of granulosa cells (GCs) in pig ovaries. Moreover, miR-183 inhibited the synthesis of estradiol in GCs and promoted the synthesis of progesterone. These results will help in gaining understanding of the role of miRNAs in regulating porcine litter size.

Effect of liver receptor homolog-1 on cell apoptosis and steroid hormone secretion on granulosa cells from Hu sheep

The apoptosis of granulosa cells can result in follicular atresia, a key process in follicle selection and development. The related molecular mechanisms were explored to study the effects of liver receptor homolog-1 (Lrh-1) on apoptosis, sodium-titanium system, and steroid synthesis of granulosa cells from Hu sheep in vitro. After constructing, designing, and synthesizing the Lrh-1 overexpression vector, liposomes were used to transfect Hu sheep granulosa cells. Thereafter, qRT-PCR and Western blot were used to detect the mechanism of Lrh-1 on the titanium system and steroid synthesis of Hu sheep granulosa cells. The overexpression efficiency was determined by fluorescence 48 h after liposome transfection into the Hu sheep granulosa cells in vitro culture. The transfection efficiency was higher, hence providing a basis for subsequent experiments. However, the protein level and transcriptional level of Lrh-1 significantly increased after transfection with Lrh-1 overexpression vector in Hu sheep granulosa cells. We further revealed that after Lrh-1 overexpression in Hu sheep granulosa cells, the expression of the pro-apoptotic gene Bax decreased significantly, while that of the anti-apoptotic gene Bcl-2 increased significantly, as well as the sodium peptide system. Moreover, the expression levels of natriuretic peptide precursor A, B, and C (NPPA, NPPB, NPPC) all showed an upward trend, while the expression levels of steroid synthesis-related genes (P450arom and P450scc) decreased. The above results showed that Lrh-1 overexpression influenced the apoptosis, sodium peptide system, and steroid synthesis.

Exploration of molecular features of PCOS with different androgen levels and immune-related prognostic biomarkers associated with implantation failure

BACKGROUND

Polycystic ovary syndrome (PCOS), the most common heterogeneous reproductive disease afflicting women of childbearing age, has been recognized as a chronic inflammatory disease recently. Most PCOS patients have hyperandrogenism, indicating a poor prognosis and poor pregnancy outcomes. The molecular mechanism underlying PCOS development is still unknown. In the present study, we investigated the gene expression profiling characteristics of PCOS with hyperandrogenism (HA) or without hyperandrogenism (NHA) and identified immune-related factors that correlated with embryo implantation failure.

METHODS

PCOS and recurrent implantation failure (RIF) microarray datasets were obtained from the Gene Expression Omnibus (GEO) database. ClueGO software was used to perform enrichment analysis of differentially expressed genes (DEGs) in PCOS with varying androgen levels. The Weighted Co-Expression Network Analysis (WGCNA) was used to identify co-expressed modules and shared gene signatures between HA PCOS and RIF. Moreover, the upregulated DEGs of HA PCOS and RIF were intersected with shared gene signatures screening by WGCNA to excavate further key prognostic biomarkers related to implantation failure of HA PCOS. The selected biomarker was verified by qRT-PCR.

RESULTS

A total of 271 DEGs were found in HA PCOS granulosa cell samples, and 720 DEGs were found in NHA PCOS. According to CuleGO enrichment analysis, DEGs in HA PCOS are enriched in immune activation and inflammatory response. In contrast, DEGs in NHA PCOS are enriched in mesenchymal cell development and extracellular space. Using WGCNA analysis, we discovered 26 shared gene signatures between HA PCOS and RIF, which were involved in corticosteroid metabolism, bone maturation and immune regulation. DAPK2 was furtherly screened out and verified to be closely related with the development of HA PCOS, acting as an independent predictor biomarker of the embryo implantation failure. DAPK2 expression was negatively correlated to the embryo implantation rate (r=-0.474, P=0.003). The immune infiltration results suggested that upregulated DAPK2 expression was closely related with NK cell infiltration and macrophage M2, playing an essential role in the pathogenesis of implantation failure in HA PCOS.

CONCLUSION

Our research revealed the expression profiling of PCOS with different androgen levels and identified DAPK2 as a critical prognostic biomarker for implantation failure in PCOS.

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.

NORHA, a novel follicular atresia-related lncRNA, promotes porcine granulosa cell apoptosis via the miR-183-96-182 cluster and FoxO1 axis

BACKGROUND

Follicular atresia has been shown to be strongly associated with a low follicle utilization rate and female infertility, which are regulated by many factors such as microRNAs (miRNAs), which constitute a class of noncoding RNAs (ncRNAs). However, little is known about long noncoding RNAs (lncRNAs), which constitute another ncRNA family that regulate follicular atresia.

RESULTS

A total of 77 differentially expressed lncRNAs, including 67 upregulated and 10 downregulated lncRNAs, were identified in early atretic follicles compared to healthy follicles by RNA-Sequencing. We characterized a noncoding RNA that was highly expressed in atretic follicles (NORHA). As an intergenic lncRNA, NORHA was one of the upregulated lncRNAs identified in the atretic follicles. To determine NORHA function, RT-PCR, flow cytometry and western blotting were performed, and the results showed that NORHA was involved in follicular atresia by influencing GC apoptosis with or without oxidative stress. To determine the mechanism of action, bioinformatics analysis, luciferase reporter assay and RNA immunoprecipitation assay were performed, and the results showed that NORHA acted as a 'sponge', that directly bound to the miR-183-96-182 cluster, and thus prevented its targeted inhibition of FoxO1, a major sensor and effector of oxidative stress.

CONCLUSIONS

We provide a comprehensive perspective of lncRNA regulation of follicular atresia, and demonstrate that NORHA, a novel lncRNA related to follicular atresia, induces GC apoptosis by influencing the activities of the miR-183-96-182 cluster and FoxO1 axis.

Follicular fluid steroid profile in sows: relationship to follicle size and oocyte quality†

Identification of reliable characteristics of follicle quality and developmental competence has been pursued in numerous studies, but with inconsistent outcomes. Here, we aimed to identify these characteristics by analysis of the follicular fluid (FF) steroid profile in relation to cumulus-oocyte complex (COC) morphology and follicle size, followed by molecular substantiation. Multiparous sows at weaning were used to facilitate analysis at the start of the follicular phase of the oestrus cycle. Sows with a higher average follicle size (≥5 mm vs. < 5 mm) had a higher follicular fluid β-estradiol concentration, but did not differ in other measured steroids. Sows with high compared to low percentage high-quality COCs (<70% vs. ≥70% high-quality) had follicular fluid with a higher concentration of β-estradiol, 19-norandrostenedione, progesterone, and α-testosterone, while the concentration of cortisol was lower. Transcriptome analysis of granulosa cells of healthy follicles of sows with a high percentage high-quality COCs showed higher abundance of transcripts involved in ovarian steroidogenesis (e.g., CYP19A2 and 3, POR, VEGFA) and growth (IGF1) and differential abundance of transcripts involved in granulosa cell apoptosis (e.g., GADD45A, INHBB). Differences in aromatase transcript abundance (CYP19A1, 2 and 3) were confirmed at the protein level. In addition, sows with a high percentage high-quality COCs lost less weight during lactation and had higher plasma IGF1 concentration at weaning, which may have affected COC quality. To the best of our knowledge, this study is also the first to report the relation between FF steroid profile and COC quality.

Integrative Proteomic and Phosphoproteomic Analyses of Granulosa Cells During Follicular Atresia in Porcine

Ovarian follicular atresia is a natural physiological process; however, the mechanism is not fully understood. In this study, quantitative proteomic and phosphoproteomic analyses of granulosa cells (GCs) in healthy (H), slightly atretic (SA), and atretic follicles (A) of porcine were performed by TMT labeling, enrichment of phosphopeptides, and liquid chromatography with tandem mass spectrometry (LC-MS/MS) analysis. In total, 6,201 proteins were quantified, and 4,723 phosphorylation sites of 1,760 proteins were quantified. In total, 24 (11 up, 13 down) and 50 (29 up, 21 down) proteins with a fold change (FC) > 5 were identified in H/SA and H/A, respectively. In addition, there were 20 (H/SA, up) and 39 (H/A, up) phosphosites with an FC > 7 that could serve as potential biomarkers for distinguishing different quality categories of follicles. Western blotting and immunofluorescence confirmed the reliability of the proteomic analysis. Some key proteins (e.g., MIF, beta catenin, integrin β2), phosphosites (e.g., S76 of caspase6, S22 and S636 of lamin A/C), pathways (e.g., apoptosis, regulation of actin cytoskeleton pathway), transcription factors (e.g., STAT5A, FOXO1, and BCLAF1), and kinases (e.g., PBK, CDK5, CDK12, and AKT3) involved in the atresia process were revealed via further analysis of the differentially expressed proteins (DEPs) and phosphorylated proteins (DEPPs). Further study showed that mutant caspase6 Ser76 to Ala increased the ratios of cleaved caspase6/caspase6 and cleaved caspase3/caspase3 and dephosphorylation of caspase6 at Ser76 increased cell apoptotic rate, a new potential pathway of follicular atresia. Collectively, the proteomic and phosphoproteomic profiling and functional research in the current study comprehensively analyzed the dynamic changes in protein expression and phosphorylation during follicular atresia and provided some new explanations regarding the regulation of this process.

Transcriptome Analysis of Porcine Granulosa Cells in Healthy and Atretic Follicles: Role of Steroidogenesis and Oxidative Stress

One of the main causes of female infertility is a deregulated antral follicular atresia, a process of which the underlying molecular mechanisms are largely unknown. Our objective was therefore to characterize the complex transcriptome changes in porcine granulosa cells of healthy antral (HA) and advanced antral atretic (AA) follicles, using ELISA and RNA-Seq followed by qRT-PCR and immunohistochemistry. Granulosa cell RNA-Seq data revealed 2160 differentially expressed genes, 1483 with higher and 677 with lower mRNA concentrations in AA follicles. Bioinformatic analysis showed that the upregulated genes in AA follicles were highly enriched in inflammation and apoptosis processes, while the downregulated transcripts were mainly highlighted in the steroid biosynthesis pathway and response to oxidative stress processes including antioxidant genes (e.g., GSTA1, GCLC, GCLM, IDH1, GPX8) involved in the glutathione metabolism pathway and other redox-related genes (e.g., RRM2B, NDUFS4). These observations were confirmed by RT-qPCR and immunohistochemistry. Additionally, the granulosa cells of AA follicles express significantly stronger 8-OHdG immunostaining, a marker of oxidative DNA damage, implicating that oxidative stress may participate in follicular atresia. We hypothesize that the decrease in anti-apoptotic factors and steroid hormones coincides with increased oxidative stress markers and the expression of pro-apoptotic factors, all contributing to antral follicular atresia.

Gene analysis of major signaling pathways regulated by gonadotropins in human ovarian granulosa tumor cells (KGN)†

The female reproductive function largely depends on timing and coordination between follicle-stimulating hormone (FSH) and luteinizing hormone. Even though it was suggested that these hormones act on granulosa cells via shared signaling pathways, mainly protein kinases A, B, and C (PKA, PKB, and PKC), there is still very little information available on how these signaling pathways are regulated by each hormone to provide such differences in gene expression throughout folliculogenesis. To obtain a global picture of the principal upstream factors involved in PKA, PKB, and PKC signaling in granulosa cells, human granulosa-like tumor cells (KGN) were treated with FSH or specific activators (forskolin, SC79, and phorbol 12-myristate 13-acetate) for each pathway to analyze gene expression with RNA-seq technology. Normalization and cutoffs (FC 1.5, P ≤ 0.05) revealed 3864 differentially expressed genes between treatments. Analysis of major upstream regulators showed that PKA is a master kinase of early cell differentiation as its activation resulted in the gene expression profile that accompanies granulosa cell differentiation. Our data also revealed that the activation of PKC in granulosa cells is also a strong differentiation signal that could control "advanced" differentiation in granulosa cells and the inflammatory cascade that occurs in the dominant follicle. According to our results, PKB activation provides support for PKA-stimulated gene expression and is also involved in granulosa cell survival throughout follicular development. Taken together, our results provide new information on PKA, PKB, and PKC signaling pathways and their roles in stimulating a follicle at the crossroad between maturation/ovulation and atresia.

Cortical Granule Distribution and Expression Pattern of Genes Regulating Cellular Component Size, Morphogenesis, and Potential to Differentiation are Related to Oocyte Developmental Competence and Maturational Capacity In Vivo and In Vitro

Polyspermia is an adverse phenomenon during mammalian fertilization when more than one sperm fuses with a single oocyte. The egg cell is prepared to prevent polyspermia by, among other ways, producing cortical granules (CGs), which are specialized intracellular structures containing enzymes that aim to harden the zona pellucida and block the fusion of subsequent sperm. This work focused on exploring the expression profile of genes that may be associated with cortical reactions, and evaluated the distribution of CGs in immature oocytes and the peripheral density of CGs in mature oocytes. Oocytes were isolated and then processed for in vitro maturation (IVM). Transcriptomic analysis of genes belonging to five ontological groups has been conducted. Six genes showed increased expression after IVM (ARHGEF2, MAP1B, CXCL12, FN1, DAB2, and SOX9), while the majority of genes decreased expression after IVM. Using CG distribution analysis in immature oocytes, movement towards the cortical zone of the oocyte during meiotic competence acquisition was observed. CGs peripheral density decreased with the rise in meiotic competence during the IVM process. The current results reveal important new insights into the in vitro maturation of oocytes. Our results may serve as a basis for further studies to investigate the cortical reaction of oocytes.

The effects of LH inhibition with cetrorelix on cumulus cell gene expression during the luteal phase under ovarian coasting stimulation in cattle

Cumulus cells have an important role to play in the final preparation of the oocyte before ovulation. During the final phase of follicular differentiation, FSH levels are low and LH maintains follicular growth; however, it is not known if at that time LH has an influence on cumulus cells inside the follicle. In humans, LH is often inhibited to avoid a premature ovulatory LH surge. This procedure provides a tool to investigate the role of LH in follicular development. In this study, we investigated the impact of suppressing LH using the GnRH antagonist cetrorelix during an ovarian coasting stimulation protocol on the transcriptome of bovine cumulus cells (CC). Oocytes were collected twice from 6 dairy cows. For the first collection, the cows received FSH twice daily for 3 d, followed by FSH withdrawal for 68 h as a control protocol. For the second collection, the same stimulation protocol was used, but the cows also received, starting on day 2 of FSH stimulation, a GnRH antagonist once a day until recovery of the cumulus-oocyte complexes (COC). Half of the COC were subjected to in vitro maturation, fertilization, and culture to assess blastocyst rates. The other half of the COC underwent microarray analysis (n = 3 cows, 2 treatments, 6 oocyte collections) and qRT-PCR (n = 6 cows: 3 microarray cows +3 other cows, 2 treatments, 12 oocyte collections). The differential expression of specific genes was confirmed by RT-qPCR: decrease of ATP6AP2, SC4MOL, and OSTC and increase of PTGDS in the LH-inhibited condition. The global transcriptomic analysis of cumulus cells demonstrated that the inhibition of LH secretion may decrease survival and growth of the follicle. Moreover, the results suggested that LH may be important to cumulus for the maintenance of cellular mechanisms such as global RNA expression, protein and nucleic acid metabolism, and energy production. These results support the hypothesis that LH support is important during the final part of follicle maturation through its influence on the cumulus cells.

Transcriptome analyses of ovarian stroma: tunica albuginea, interstitium and theca interna

The ovary has specialised stromal compartments, including the tunica albuginea, interstitial stroma and theca interna, which develops concurrently with the follicular antrum. To characterise the molecular determinants of these compartments, stroma adjacent to preantral follicles (pre-theca), interstitium and tunica albuginea were laser microdissected (n = 4 per group) and theca interna was dissected from bovine antral follicles (n = 6). RNA microarray analysis showed minimal differences between interstitial stroma and pre-theca, and these were combined for some analyses and referred to as stroma. Genes significantly upregulated in theca interna compared to stroma included INSL3, LHCGR, HSD3B1, CYP17A1, ALDH1A1, OGN, POSTN and ASPN. Quantitative RT-PCR showed significantly greater expression of OGN and LGALS1 in interstitial stroma and theca interna versus tunica and greater expression of ACD in tunica compared to theca interna. PLN was significantly higher in interstitial stroma compared to tunica and theca. Ingenuity pathway, network and upstream regulator analyses were undertaken. Cell survival was also upregulated in theca interna. The tunica albuginea was associated with GPCR and cAMP signalling, suggesting tunica contractility. It was also associated with TGF-β signalling and increased fibrous matrix. Western immunoblotting was positive for OGN, LGALS1, ALDH1A1, ACD and PLN with PLN and OGN highly expressed in tunica and interstitial stroma (each n = 6), but not in theca interna from antral follicles (n = 24). Immunohistochemistry localised LGALS1 and POSTN to extracellular matrix and PLN to smooth muscle cells. These results have identified novel differences between the ovarian stromal compartments.

Activation of Steroidogenesis, Anti-Apoptotic Activity, and Proliferation in Porcine Granulosa Cells by RUNX1 Is Negatively Regulated by H3K27me3 Transcriptional Repression

H3K27me3 is an epigenetic modification that results in the repression of gene transcription. The transcription factor RUNX1 (the runt-related transcription factor 1) influences granulosa cells' growth and ovulation. This research uses ELISA, flow cytometry, EDU, ChIP-PCR, WB and qPCR to investigate steroidogenesis, cell apoptosis, and the proliferation effect of RUNX1 in porcine granulosa cells (pGCs) as regulated by H3K27me3. Decreased H3K27me3 stimulates the expression of steroidogenesis-related genes, including CYP11A1, PTGS2, and STAR, as well as prostaglandin. H3K27me3 transcriptionally represses RUNX1 here, whereas RUNX1 acts as an activator of FSHR, CYP11A1, and CYP19A1, promoting the production of androgen, estrogen, and prostaglandin, as well as increasing anti-apoptotic and cell proliferation activity, but decreasing progesterone. Both the complementary recovery of the H3K27me3 antagonist with the siRUNX1 signal, and the H3K27me3 agonist with the RUNX1 signal to maintain RUNX1 lead to the activation of CYP19A1, ER1, HSD17β4, and STAR here. Androgen and prostaglandin are significantly repressed but progesterone is markedly increased with the antagonist and siRUNX1. Prostaglandin is significantly promoted with the agonist and RUNX1. Furthermore, H3K27me3-RUNX1 affects the anti-apoptotic activity and stimulation of proliferation in pGCs. The present work verifies the transcriptional suppression of RUNX1 by H3K27me3 during antral follicular development and maturation, which determines the levels of hormone synthesis and cell apoptosis and proliferation in the pGC microenvironment.

Transcriptomic analysis of expression of genes regulating cell cycle progression in porcine ovarian granulosa cells during short-term in vitro primary culture

The primary function of ovarian granulosa cells (GCs) is the support of oocytes during maturation and development. Molecular analyses of granulosa cell-associated processes, leading to improvement of understanding of the cell cycle events during the formation of ovarian follicles (folliculogenesis), may be key to improve the in vitro fertilization procedures. Primary in vitro culture of porcine GCs was employed to examine the changes in the transcriptomic profile of genes belonging to “cell cycle”, “cell division”, “cell cycle process”, “cell cycle phase transition”, “cell cycle G1/S phase transition”, “cell cycle G2/M phase transition” and “cell cycle checkpoint” ontology groups. During the analysis, microarrays were employed to study the transcriptome of GCs, analyzing the total RNA of cells from specific periods of in vitro cultures. This research was based on material obtained from 40 landrace gilts of similar weight, age and the same living conditions. RNA was isolated at specific timeframes: before the culture was established (0 h) and after 48 h, 96 h and 144 h in vitro. Out of 133 differentially expressed genes, we chose the 10 most up-regulated (SFRP2, PDPN, PDE3A, FGFR2, PLK2, THBS1, ETS1, LIF, ANXA1, TGFB1) and the 10 most downregulated (IGF1, NCAPD2, CABLES1, H1FOO, NEK2, PPAT, TXNIP, NUP210, RGS2 and CCNE2). Some of these genes known to play key roles in the regulation of correct cell cycle passage (up-regulated SFRP2, PDE3A, PLK2, LIF and down-regulated CCNE2, TXNIP, NEK2). The data obtained provide a potential reference for studies on the process of mammalian folliculogenesis, as well as suggests possible new genetic markers for cell cycle progress in in vitro cultured porcine granulosa cells.

Transactivation of miR-202-5p by Steroidogenic Factor 1 (SF1) Induces Apoptosis in Goat Granulosa Cells by Targeting TGFβR2

MicroRNAs play key roles during ovary development, with emerging evidence suggesting that miR-202-5p is specifically expressed in female animal gonads. Granulosa cells (GCs) are somatic cells that are closely related to the development of female gametes in mammalian ovaries. However, the biological roles of miR-202-5p in GCs remain unknown. Here, we show that miR-202-5p is specifically expressed in GCs and accumulates in extracellular vesicles (EVs) from large growth follicles in goat ovaries. In vitro assays showed that miR-202-5p induced apoptosis and suppressed the proliferation of goat GCs. We further revealed that miR-202-5p is a functional miRNA that targets the transforming growth factor-beta type II receptor (TGFβR2). MiR-202-5p attenuated TGF-β/SMAD signaling through the degradation of TGFβR2 at both the mRNA and protein level, decreasing p-SMAD3 levels in GCs. Moreover, we verified that steroidogenic factor 1 (SF1) is a transcriptional factor that binds to the promoters of miR-202 and cytochrome P450 family 19 subfamily A member 1 (CYP19A1) through luciferase reporter and chromatin immunoprecipitation (ChIP) assays. That contributed to positive correlation between miR-202-5p and CYP19A1 expression and estradiol (E2) release. Furthermore, SF1 repressed TGFβR2 and p-SMAD3 levels in GCs through the transactivation of miR-202-5p. Taken together, these results suggest a mechanism by which miR-202-5p regulates canonical TGF-β/SMAD signaling through targeting TGFβR2 in GCs. This provides insight into the transcriptional regulation of miR-202 and CYP19A1 during goat ovarian follicular development.

Dynamics and Regulations of BimEL Ser65 and Thr112 Phosphorylation in Porcine Granulosa Cells during Follicular Atresia

BimEL protein is involved in follicular atresia by regulating granulosa cell apoptosis, but the dynamic changes of BimEL phosphorylation during follicular atresia are poorly understood. The aim of this study was to explore the changes of key BimEL phosphorylation sites and their upstream regulatory pathways. First, the levels of BimEL-Ser65 and BimEL-Thr112 phosphorylation (p-BimEL-S65, p-BimEL-T112) in granulosa cells (GC) from healthy (H), slightly-atretic (SA), and atretic (A) follicles and in cultured GC after different treatments were detected by Western blotting. Next, the effects of the corresponding site mutations of BIM on apoptosis of GC were investigated. Finally, the pathways of two phosphorylation sites were investigated by kinase inhibitors. The results revealed that p-BimEL-S65 levels were higher in GC from H than SA and A, whereas p-BimEL-T112 was reversed. The prosurvival factors like FSH and IGF-1 upregulated the level of p-BimEL-S65, while the proapoptotic factor, heat stress, increased the level of p-BimEL-T112 in cultured GC. Compared with the overexpression of wild BimEL, the apoptotic rate of the GC overexpressed BimEL-S65A (replace Ser65 with Ala) mutant was significantly higher, but the apoptotic rate of the cells overexpressing BimEL-T112A did not differ. In addition, inhibition of the ERK1/2 or JNK pathway by specific inhibitors reduced the levels of p-BimEL-S65 and p-BimEL-T112. In conclusion, the levels of p-BimEL-S65 and p-BimEL-T112 were reversed during follicular atresia. Prosurvival factors promote p-BimEL-S65 levels via ERK1/2 to inhibit GC apoptosis, whereas proapoptotic factor upregulates the level of p-BimEL-T112 via JNK to induce GC apoptosis.

‘Cell cycle’ and ‘cell death’- related genes are differentially expressed during long – term in vitro real-time cultivation of porcine oviductal epithelial cells

Alterations in cells depend on their genetic material, its activation and translation of the products. The genes responsible for the cell cycle processes and apoptosis of porcine oviductal cells have been presented in our study. The processes occurring in the reproductive system of females are extremely complex and require in-depth knowledge. Thanks to in vitro studies on the fallopian tube epithelium cells, we can get closer to understanding the biochemical and morphological changes occurring in mammalian organisms. Our research was conducted on fallopian tubes obtained from commercially bred pigs and its aim was to assess the expression profile of genes responsible for the most important processes of cellular life. Cell cultures were carried out for 30 days, with the obtained cells subjected to molecular analysis. We have shown significant regulation of “cell death” and “cell cycle” genes, some of which are related to the reproductive system. The alterations in transcriptomic profile and mutual relations between the genes were analyzed and related to the literature findings. The knowledge gained could help in identifying new potential markers of the in vitro occurrence of processes described by the ontology groups of interest.

Running title: pig, oocytes, microarray assays, in vitro maturation (IVM)

Detection of insertions/deletions (InDels) within the goat Runx2 gene and their association with litter size and growth traits

Runt-related transcription factor 2 (Runx2) is characterized by its critical functions in osteoblastic and ovulatory processes. The goal of this study was to explore the insertion/deletion (indel) variants of this gene and to evaluate their association with productive traits. Herein, a 12 bp and 6 bp insertion within the Runx2 gene was uncovered in Shaanbei white cashmere goats (SBWC; n = 1200). Chi-square analysis revealed that the 12 bp insertion was related to litter size (p < 0.01). Further association analysis also found this insertion was significantly associated with litter size (p = 1.1E-5). Interestingly, this insertion was also significantly associated with chest circumference (p = 0.018). Additionally, the 6 bp insertion was associated with body length (p = 0.003), chest width (p = 0.011), and chest circumference (p = 0.005). Furthermore, diplotype associations also uncovered that the combined genotypes of these two indels also significantly affected litter size and growth traits (p < 0.05). These findings suggested that these two insertions within the Runx2 gene were significantly associated with reproduction and growth traits, which would make them beneficial functional DNA markers that can be used in goat breeding.

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.

Analysis of protein-protein interaction network based on transcriptome profiling of ovine granulosa cells identifies candidate genes in cyclic recruitment of ovarian follicles

After pubertal, cohort of small antral follicles enters to gonadotrophin-sensitive development, called recruited follicles. This study was aimed to identify candidate genes in follicular cyclic recruitment via analysis of protein-protein interaction (PPI) network. Differentially expressed genes (DEGs) in ovine granulosa cells of small antral follicles between follicular and luteal phases were accumulated among gene/protein symbols of the Ensembl annotation. Following directed graphs, PTPN6 and FYN have the highest indegree and outdegree, respectively. Since, these hubs being up-regulated in ovine granulosa cells of small antral follicles during the follicular phase, it represents an accumulation of blood immune cells in follicular phase in comparison with luteal phase. By contrast, the up-regulated hubs in the luteal phase including CDK1, INSRR and TOP2A which stimulated DNA replication and proliferation of granulosa cells, they known as candidate genes of the cyclic recruitment.

Initiation of follicular atresia: gene networks during early atresia in pig ovaries

In mammals, more than 99% of ovarian follicles undergo a degenerative process known as atresia. The molecular events involved in atresia initiation remain incompletely understood. The objective of this study was to analyze differential gene expression profiles of medium antral ovarian follicles during early atresia in pig. The transcriptome evaluation was performed on cDNA microarrays using healthy and early atretic follicle samples and was validated by quantitative PCR. Annotation analysis applying current database (Sus scrofa 11.1) revealed 450 significantly differential expressed genes between healthy and early atretic follicles. Among them, 142 were significantly upregulated in early atretic with respect to healthy group and 308 were downregulated. Similar expression trends were observed between microarray data and quantitative RT-PCR confirmation, which indicated the reliability of the microarray analysis. Further analysis of the differential expressed genes revealed the most significantly affected biological functions during early atresia including blood vessel development, regulation of DNA-templated transcription in response to stress and negative regulation of cell adhesion. The pathway and interaction analysis suggested that atresia initiation associates with (1) a crosstalk of cell apoptosis, autophagy and ferroptosis rather than change of typical apoptosis markers, (2) dramatic shift of steroidogenic enzymes, (3) deficient glutathione metabolism and (4) vascular degeneration. The novel gene candidates and pathways identified in the current study will lead to a comprehensive view of the molecular regulation of ovarian follicular atresia and a new understanding of atresia initiation.

Genome-Wide Association Analyses Highlight the Potential for Different Genetic Mechanisms for Litter Size Among Sheep Breeds

Reproduction is an important trait in sheep breeding as well as in other livestock. However, despite its importance the genetic mechanisms of litter size in domestic sheep (Ovis aries) are still poorly understood. To explore genetic mechanisms underlying the variation in litter size, we conducted multiple independent genome-wide association studies in five sheep breeds of high prolificacy (Wadi, Hu, Icelandic, Finnsheep, and Romanov) and one low prolificacy (Texel) using the Ovine Infinium HD BeadChip, respectively. We identified different sets of candidate genes associated with litter size in different breeds: BMPR1B, FBN1, and MMP2 in Wadi; GRIA2, SMAD1, and CTNNB1 in Hu; NCOA1 in Icelandic; INHBB, NF1, FLT1, PTGS2, and PLCB3 in Finnsheep; ESR2 in Romanov and ESR1, GHR, ETS1, MMP15, FLI1, and SPP1 in Texel. Further annotation of genes and bioinformatics analyses revealed that different biological pathways could be involved in the variation in litter size of females: hormone secretion (FSH and LH) in Wadi and Hu, placenta and embryonic lethality in Icelandic, folliculogenesis and LH signaling in Finnsheep, ovulation and preovulatory follicle maturation in Romanov, and estrogen and follicular growth in Texel. Taken together, our results provide new insights into the genetic mechanisms underlying the prolificacy trait in sheep and other mammals, suggesting targets for selection where the aim is to increase prolificacy in breeding projects.

Identification of Differentially Expressed Genes in Porcine Ovaries at Proestrus and Estrus Stages Using RNA-Seq Technique

Estrus is an important factor for the fecundity of sows, and it is involved in ovulation and hormone secretion in ovaries. To better understand the molecular mechanisms of porcine estrus, the expression patterns of ovarian mRNA at proestrus and estrus stages were analyzed using RNA sequencing technology. A total of 2,167 differentially expressed genes (DEGs) were identified (P ≤ 0.05, |log2  Ratio| ≥ 1), of which 784 were upregulated and 1,383 were downregulated in the estrus compared with the proestrus group. Gene Ontology (GO) enrichment indicated that these DEGs were mainly involved in the cellular process, single-organism process, cell and cell part, and binding and metabolic process. In addition, a pathway analysis showed that these DEGs were significantly enriched in 33 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including cell adhesion molecules, ECM-receptor interaction, and cytokine-cytokine receptor interaction. Quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) confirmed the differential expression of 10 selected DEGs. Many of the novel candidate genes identified in this study will be valuable for understanding the molecular mechanisms of the sow estrous cycle.

miR-1275 controls granulosa cell apoptosis and estradiol synthesis by impairing LRH-1/CYP19A1 axis

miR-1275 is one of the microRNAs (miRNAs) that are differentially expressed during follicular atresia in pig ovaries, as identified by a miRNA microarray assay in our previous study [1]. However, its functions in follicular atresia remain unknown. In this study, we showed that miR-1275 promotes early apoptosis of porcine granulosa cells (pGCs) and the initiation of follicular atresia, and inhibits E2 release and expression of CYP19A1, the key gene in E2 production. Bioinformatics and luciferase reporter assays revealed that liver receptor homolog (LRH)-1, not CYP19A1, is a direct functional target of miR-1275. In vitro overexpression and knockdown experiments showed that LRH-1 significantly repressed apoptosis and induced E2 secretion and CYP19A1 expression in pGCs. LRH-1, whose expression was regulated by miR-1275, prevented apoptosis in pGCs. Furthermore, luciferase and chromatin immunoprecipitation assays demonstrated that LRH-1 protein bound to the CYP19A1 promoter and increased its activity. Our findings suggest that miR-1275 attenuates LRH-1 expression by directly binding to its 3'UTR. This prevents the interaction of LRH-1 protein with the CYP19A1 promoter, represses E2 synthesis, promotes pGC apoptosis, and initiates follicular atresia in porcine ovaries.

Transcriptome analysis of ovine granulosa cells reveals differences between small antral follicles collected during the follicular and luteal phases

Ovarian follicular growth occurs in both the follicular and luteal phases of the estrous cycle but in very different endocrine contexts. In both phases, many small antral follicles with similar morphologic and histologic characteristics are present within the ovaries as a reserve for the terminal folliculogenesis. However, there are several gaps in our molecular knowledge of the gene expression profiles of small antral follicles in the follicular and luteal phases. The aim of the present study was to use RNA sequencing to compare and analyze the global transcriptional profile of ovine granulosa cells collected from small antral follicles (1-3 mm) either during the follicular or the luteal phase of the estrous cycle, with the hypothesis that they should be differential. We identified 663 genes whose mRNA was differentially expressed or accumulated in the granulosa cell layer of small antral follicles in the two phases. A comprehensive interpretation of these data was performed through integrative analyses (Gene Ontology, Ingenuity Pathway Analysis) and the exploitation of already available transcriptomic data on follicular growth and atresia. In particular, we observed that the contrasted endocrine context between follicular and luteal phases may have an impact on estradiol, follicle-stimulating hormone (FSH), and on the activin/inhibin signaling pathways. Furthermore, we reveal the possible initiation of early follicular atresia in small antral follicles during the follicular phase in interaction with the presence of immune cells. This study provides new insights into the gene expression profile in ovine granulosa cells, and we suggest that these molecular changes may have an implication at the time of follicle selection.