Characterization of the human cumulus cell transcriptome during final follicular maturation and ovulation

Single-cell RNA-seq reveals the effects of the FecB mutation on the transcriptome profile in ovine cumulus cells

Genetic variations in the ovine ovulation rate, which are associated with the FecB mutation, provide useful models by which to explore the mechanisms regulating the development of mammalian antral follicles. In order to study the effects of the FecB mutation on cumulus cell differentiation, preovulatory follicles were aspirated and cumulus cells were isolated from three FecB genotypes (homozygous, heterozygous and wild type) of Small Tail Han (STH) sheep superstimulated with FSH. Transcriptome information from tens of thousands of cumulus cells was determined with the 10 × Genomics single-cell RNA-seq technology. Under the superovulation treatment, the observed number of preovulatory follicles in the ovaries of FecB carriers was still significantly higher than that in the wild-type (P < 0.05). The expression patterns of cumulus cells differed between FecB carriers and wild-type ewes. The screened cumulus cells could also be further divided into different cell clusters, and the differentiation states and fates of each group of cumulus cells also remained different, which supports the notion that heterogeneity in gene expression is prevalent in single cells. The oxidative phosphorylation pathway was significantly enriched in differentially expressed genes among the cell differentiation branch nodes of cumulus cells and among the differentially expressed genes of cumulus cells from the three genotypes. Combined with the important role of oxidative phosphorylation in the maturation of COCs, we suggest that the oxidative phosphorylation pathway of cumulus cells plays a crucial role in the differentiation process of cumulus cells and the mutation effect of the FecB gene.

moSCminer: a cell subtype classification framework based on the attention neural network integrating the single-cell multi-omics dataset on the cloud

Single-cell omics sequencing has rapidly advanced, enabling the quantification of diverse omics profiles at a single-cell resolution. To facilitate comprehensive biological insights, such as cellular differentiation trajectories, precise annotation of cell subtypes is essential. Conventional methods involve clustering cells and manually assigning subtypes based on canonical markers, a labor-intensive and expert-dependent process. Hence, an automated computational prediction framework is crucial. While several classification frameworks for predicting cell subtypes from single-cell RNA sequencing datasets exist, these methods solely rely on single-omics data, offering insights at a single molecular level. They often miss inter-omic correlations and a holistic understanding of cellular processes. To address this, the integration of multi-omics datasets from individual cells is essential for accurate subtype annotation. This article introduces moSCminer, a novel framework for classifying cell subtypes that harnesses the power of single-cell multi-omics sequencing datasets through an attention-based neural network operating at the omics level. By integrating three distinct omics datasets-gene expression, DNA methylation, and DNA accessibility-while accounting for their biological relationships, moSCminer excels at learning the relative significance of each omics feature. It then transforms this knowledge into a novel representation for cell subtype classification. Comparative evaluations against standard machine learning-based classifiers demonstrate moSCminer's superior performance, consistently achieving the highest average performance on real datasets. The efficacy of multi-omics integration is further corroborated through an in-depth analysis of the omics-level attention module, which identifies potential markers for cell subtype annotation. To enhance accessibility and scalability, moSCminer is accessible as a user-friendly web-based platform seamlessly connected to a cloud system, publicly accessible at http://203.252.206.118:5568. Notably, this study marks the pioneering integration of three single-cell multi-omics datasets for cell subtype identification.

Long noncoding RNAs in human reproductive processes and diseases

Infertility has become a global disease burden. Although assisted reproductive technologies are widely used, the assisted reproduction birth rate is no more than 30% worldwide. Therefore, understanding the mechanisms of reproduction can provide new strategies to improve live birth rates and clinical outcomes of enhanced implantation. Long noncoding RNAs (lncRNAs) have been reported to exert regulatory roles in various biological processes and diseases in many species. In this review, we especially focus on the role of lncRNAs in human reproduction. We summarize the function and mechanisms of lncRNAs in processes vital to reproduction, such as spermatogenesis and maturation, sperm motility and morphology, follicle development and maturation, embryo development and implantation. Then, we highlight the importance and diverse potential of lncRNAs as good diagnostic molecular biomarkers and therapeutic targets for infertility treatment.

The expression levels of NOS2, HMOX1, and VEGFC in cumulus cells are markers of oocyte maturation and fertilization rate

Throughout the reproductive life of women, cumulus cells (CC) protect the dormant oocyte from damage, act as sensors of the follicular microenvironment, and act as a gatekeeper for oocyte developmental potential. One such mechanism relies on the hypoxia-tolerance response, which, with age, decreases systematically, including in the ovary. We aimed to evaluate the association between gene expression related to hypoxia and aging in CC and reproductive results in in vitro fertilization cycles. We recruited 94 women undergoing controlled ovarian stimulation. Total RNA was extracted from pooled CCs collected after oocyte pick-up (OPU) and reverse-transcribed to complementary DNA using random hexamers to test 14 genes related to hypoxia response via HIF1α activation, oxidative stress, and angiogenic responses. The expression of CLU, NOS2, and TXNIP had a positive correlation with age (rs  = 0.25, rs  = 0.24, and rs  = 0.35, respectively). Additionally, NOS2 and HMOX1 expression correlated positively with the retrieval of immature oocytes (rs  = 0.22 and rs  = 0.40, respectively). Moreover, VEGFC levels decreased overall with increasing fertilization rate, independently of age (rs  = -0.29). We found that the fertilization potential of a cohort of oocytes is related to the ability of CC to respond to oxidative stress and hypoxia with age, pointing at NOS2, HMOX1, and VEGFC expression as markers for oocyte maturation and fertilization success.

ABCC4 is a PGE2 efflux transporter in the ovarian follicle: A mediator of ovulation and a potential non-hormonal contraceptive target

The role of prostaglandins (PGs) in the ovulatory process is known. However, the role of the ATP binding cassette subfamily C member 4 (ABCC4), transmembrane PG carrier protein, in ovulation remains unknown. We report herein that ABCC4 expression is significantly upregulated in preovulatory human granulosa cells (GCs). We found that PGE2 efflux in cultured human GCs is mediated by ABCC4 thus regulating its extracellular concentration. The ABCC4 inhibitor probenecid demonstrated effective blocking of ovulation and affects key ovulatory genes in female mice in vivo. We postulate that the reduction in PGE2 efflux caused by the inhibition of ABCC4 activity in GCs decreases the extracellular concentration of PGE2 and its ovulatory effect. Treatment of female mice with low dose of probenecid as well as with the PTGS inhibitor indomethacin or Meloxicam synergistically blocks ovulation. These results support the hypothesis that ABCC4 has an important role in ovulation and might be a potential target for non-hormonal contraception, especially in combination with PGE2 synthesis inhibitors. These findings may fill the gap in understanding the role of ABCC4 in PGE2 signaling, enhance the understanding of ovulatory disorders, and facilitate the treatment and control of fertility.

RNA sequencing-based transcriptome analysis of granulosa cells from follicular fluid: Genes involved in embryo quality during in vitro fertilization and embryo transfer

BACKGROUND

Granulosa cells play an important role in folliculogenesis, however, the role of RNA transcripts of granulosa cells in assessing embryo quality remains unclear. Therefore, we aims to investigate that RNA transcripts of granulosa cells be used to assess the probability of the embryonic developmental capacity.

METHODS

This prospective cohort study was attempted to figure out the probability of the embryonic developmental capacity using RNA sequencing of granulosa cells. Granulosa cells were collected from 48 samples in good-quality embryo group and 79 in only poor- quality embryo group from women undergoing in vitro fertilization and embryo transfer treatment. Three samples from each group were used for RNA sequencing.

RESULTS

226 differentially expressed genes (DEGs) were related to high developmental competence of embryos. Gene Ontology enrichment analysis indicated that these DEGs were primarily involved in biological processes, molecular functions, and cellular components. Additionally, pathway analysis revealed that these DEGs were enriched in 13 Kyoto Encyclopedia of Genes and Genomes pathways. Reverse transcription quantitative polymerase chain reaction verified the differential expression of the 13 selected DEGs. Among them,10 genes were differently expressed in the poor-quality embryo group compared to good-quality embryo group, including CSF1R, CTSH, SERPINA1, CYP27A1, ITGB2, IL1β, TNF, TAB1, BCL2A1, and CCL4.

CONCLUSIONS

RNA sequencing data provide the support or confute granulosa expressed genes as non-invasive biomarkers for identifying the embryonic developmental capacity.

Long Non-Coding RNA GDAR Regulates Ovine Granulosa Cells Apoptosis by Affecting the Expression of Apoptosis-Related Genes

Short-term dietary supplementation of ewes during the luteal phase can increase fertility, most probably by stimulating glucose uptake by the follicles. However, the molecular mechanism of glucose regulation of follicular development has not yet been clarified, especially the further study of long non-coding RNA (lncRNA) in determining fertility during follicular development. We generated granulosa cell (GC) models of different doses of glucose (0, 2.1, 4.2, 8.4, 16.8 and 33.6 mM), and observed that the highest cell viability was recorded in the 8.4 mM group and the highest apoptosis rates were recorded in the 33.6 mM group. Therefore, a control group (n = 3, 0 mM glucose), a low glucose group (n = 3, add 8.4 mM glucose), and a high glucose group (n = 3, add 33.6 mM glucose) of GCs were created for next whole genomic RNA sequencing. In total, 18,172 novel lncRNAs and 510 annotated lncRNAs were identified in the GCs samples. Gene Ontology indicated that differentially expressed lncRNAs associated with cell apoptosis were highly enriched. Kyoto Encyclopedia of Genes and Genomes enrichment analysis of lncRNA target genes found that the apoptosis pathway and the p53 signaling pathway were both enriched. Furthermore, we focused on the function of a lncGDAR and verified that lncGDAR could influence cell apoptosis in GC development through affecting the mRNA and protein expression of apoptosis-related markers. These results provide the basis for further study of the lncRNA regulation mechanism in nutrition on female fertility.

Down-regulation of long non-coding RNAs in reproductive aging and analysis of the lncRNA-miRNA-mRNA networks in human cumulus cells

PURPOSE

Long non-coding RNAs (lncRNAs) control gene expression at multiple levels. By interacting with microRNAs (miRNAs), they regulate their mRNA targets creating dynamic regulatory networks involved in different cellular processes. Their role in follicle development and oocyte maturation has recently emerged. lncRNA deregulation has been found associated with different pathological conditions. In this study, we identified differentially expressed lncRNAs in cumulus cells (CCs) isolated from MII oocytes of advanced maternal age women and proposed ceRNA-networks involved in signaling pathways crucial in ovarian folliculogenesis and female germ cell maturation.

METHODS

We performed a high-throughput analysis of the expression profile of 68 lncRNAs from CCs of aged and young women by using NanoString technology. By miRNet, TarPmiR, miRTarBase, OKdb, and KEGG we predicted some ceRNA-networks involving the differentially expressed (DE) lncRNAs, miRNA interactors, and their mRNA target genes.

RESULTS

We identified 28 lncRNAs down-regulated in CC samples from aged women. The analysis revealed that the miRNAs binding 11 of the DE lncRNAs and their mRNA targets are included in ceRNA-networks involved in the regulation of the PI3K-Akt, FOXO, and p53 signaling pathways.

CONCLUSION

We proposed that the lncRNA down-regulation in CCs from aged women could influence the expression of genes encoding proteins deregulated in reproductive aging. A better understanding of the interplay of lncRNA-miRNA-mRNA networks in human CCs could increase our knowledge about the mechanisms of regulation of gene expression involved in aging, lead to the development of novel therapeutics, and improve reproductive outcomes in aged women.

Aging conundrum: A perspective for ovarian aging

Progressive loss of physiological integrity and accumulation of degenerative changes leading to functional impairment and increased susceptibility to diseases are the main features of aging. The ovary, the key organ that maintains female reproductive and endocrine function, enters aging earlier and faster than other organs and has attracted extensive attention from society. Ovarian aging is mainly characterized by the progressive decline in the number and quality of oocytes, the regulatory mechanisms of which have yet to be systematically elucidated. This review discusses the hallmarks of aging to further highlight the main characteristics of ovarian aging and attempt to explore its clinical symptoms and underlying mechanisms. Finally, the intervention strategies related to aging are elaborated, especially the potential role of stem cells and cryopreservation of embryos, oocytes, or ovarian tissue in the delay of ovarian aging.

Transcriptome analysis reveals that long noncoding RNAs contribute to developmental differences between medium-sized ovarian follicles of Meishan and Duroc sows

Ovulation rate is an extremely important factor affecting litter size in sows. It differs greatly among pig breeds with different genetic backgrounds. Long non-coding RNAs (lncRNAs) can regulate follicle development, granulosa cell growth, and hormone secretion, which in turn can affect sow litter size. In this study, we identified 3554 lncRNAs and 25,491 mRNAs in M2 follicles of Meishan and Duroc sows. The lncRNA sequence and open reading frame lengths were shorter than mRNAs, and lncRNAs had fewer exons, were less abundant, and more conserved than protein-coding RNAs. Furthermore, 201 lncRNAs were differentially expressed (DE) between breeds, and quantitative trait loci analysis of DE lncRNAs were performed. A total of 127 DE lncRNAs were identified in 119 reproduction trait-related loci. In addition, the potential target genes of lncRNAs in cis or trans configurations were predicted. Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed that some potential target genes were involved in follicular development and hormone secretion-related biological processes or pathways, such as progesterone biosynthetic process, estrogen metabolic process, ovarian steroidogenesis, and PI3K-Akt signaling pathway. Furthermore, we also screened 19 differentially expressed lncRNAs in the PI3K-Akt signaling pathway as candidates. This study provides new insights into the roles of lncRNAs in follicular growth and development in pigs.

Whole-Transcriptome Analysis of LncRNAs Mediated ceRNA Regulation in Granulosa Cells Isolated From Healthy and Atresia Follicles of Chinese Buffalo

Granulosa cells (GCs) are the main supporting cells in follicles and play an important role in the regulation of oocyte maturation and follicular atresia. Accumulating evidence indicates that non-coding RNAs participate in regulation of the physiological function of GCs. However, whole-transcriptome analysis for GCs of buffalo has yet to be reported. In this study, healthy follicles (HFs) and atretic follicles (AFs) were defined according to the apoptosis rate of GCs and the hormone level in follicular fluid. GCs were collected from HFs and AFs (n = 15, 5 < n < 8 mm) for whole-transcriptome analysis using second-generation high-throughput sequencing. A total of 1,861 and 1,075 mRNAs, 159 and 24 miRNAs, and 123 and 100 lncRNAs, were upregulated and downregulated between HFs and AFs, respectively. Enrichment of functions and signaling pathways of these differentially expressed (DE) genes showed that most of DEmRNAs and targets of DEmiRNAs were annotated to the categories of ECM–receptor interaction and focal adhesion, as well as PI3K-AKT, mTOR, TGF-beta, Rap1, and estrogen signaling pathways. The competing endogenous RNA (CeRNA) network was also constructed based on the ceRNA theory which further revealed regulatory roles of these DERNAs in GCs of buffalo follicles. Finally, we validated that lnc4040 regulated the expression of Hif1a as miR-709 sponge in a ceRNA mechanism, suggesting their critical functions in GCs of buffalo follicles. These results show that lncRNAs are dynamically expressed in GCs of HFs and AFs, and interacting with target genes in a ceRNA manner, suggesting their critical functions in buffalo follicular development and atresia.

The Involvement of Lumican in Human Ovulatory Processes

Based on a previous global transcriptome sequencing project, we hypothesized that Lumican (LUM) might play a role in ovulatory processes. We sought to determine LUM gene expression under various conditions in human preovulatory follicles. The in vitro expression of LUM mRNA in mural (MGCs) and cumulus (CGCs) granulosa cells was characterized using quantitative real-time polymerase chain reaction (qRT-PCR). Immunohistochemical staining was used to identify human LUM expression in follicles at different developmental stages. Cell signaling studies were performed by treating human MGCs with human chorionic gonadotropin (hCG) and both, different stimulators and inhibitors to determine their effect on LUM expression by using qRT-PCR. Cell confluence studies were carried out to study the correlation between LUM expression and follicle cell proliferation. Follicular MGCs and CGCs of women undergoing in vitro fertilization (IVF) procedures due to endometriosis were analyzed for differences in LUM expression patterns by qRT-PCR. LUM mRNA expression was significantly higher in MGCs as compared to CGCs. In CGCs, LUM mRNA was higher in mature metaphase II (MII) oocytes than in germinal vesicle (GV) and metaphase I (MI) oocytes. LUM expression was significantly upregulated in response to hCG in cultured MGCs. Immunohistochemistry of human ovaries revealed LUM was mostly present in MGCs of large preovulatory and postovulatory follicles and absent from primordial follicles. Using pharmacological activators and inhibitors, we demonstrated that LUM induction by luteinizing hormone (LH)/hCG is carried through the mitogen-activated protein kinase (MEK) pathway. LUM expression was induced in high-density cell cultures in a confluence-dependent manner. MGCs from follicles of subjects with endometriosis exhibited reduced mRNA transcription levels compared to control subjects. Our study confirms that LUM is a newly discovered ovulatory gene. LUM might play an important role during the preovulatory period up until ovulation as well as in endometriosis infertility. A better understanding of LUM's role might provide potential new treatment paradigms for some types of female infertility.

Characterization of Long Non-Coding RNA Profiles in Porcine Granulosa Cells of Healthy and Atretic Antral Follicles: Implications for a Potential Role in Apoptosis

Long non-coding RNAs (lncRNAs) play important roles in multiple biological processes including ovarian follicular development. Here we aimed to gain novel information regarding lncRNAs transcriptome profiles in porcine granulosa cells of advanced atretic antral (AA) and healthy antral (HA) follicles using RNA-seq. A total of 11,321 lncRNAs including 10,813 novel and 508 annotated lncRNAs were identified, of which 173 lncRNAs were differentially expressed (DE-lncRNAs); ten of these were confirmed by qRT-PCR. Gene Ontology indicated that DE-lncRNAs associated with developmental processes were highly enriched. Pathway analysis demonstrated predicted cis- and trans-targets of DE-lncRNAs. Potential mRNA targets of up-regulated DE-lncRNAs were mainly enriched in apoptosis related pathways, while targeted genes of downregulated DE-lncRNAs were primarily enriched in metabolism and ovarian steroidogenesis pathways. Linear regression analyses showed that expression of upregulated DE-lncRNAs was significantly associated with apoptosis related genes. NOVEL_00001850 is the most-downregulated DE-lncRNA (FDR = 0.04, FC = -6.53), of which miRNA binding sites were predicted. KEGG analysis of its downregulated target genes revealed that ovarian steroidogenesis was the second most highlighted pathway. qRT-PCR and linear regression analysis confirmed the expression and correlation of its potential targeted gene, CYP19A1, a key gene involved in estradiol synthesis. Our results indicate that lncRNAs may participate in granulosa cells apoptosis and thus antral follicular atresia.

Human bone morphogenetic protein 8A promotes expansion and prevents apoptosis of cumulus cells in vitro

Cumulus expansion is essential for ovulation and oocyte maturation in mammals. Previous studies suggest that this process requires certain cumulus expansion enabling factors, induced by LH surge, that activate SMAD signaling locally. However, their identities remain uncertain. Using a superovulated rat model, we showed that Bmp8 transcripts were abundant in cumulus cell-oocyte complexes (COCs) and their levels can be further induced during ovulation. By analyzing human COC-related transcriptomic datasets, BMP8 transcripts in cumulus cells were also found to be significantly elevated along with the maturation status and developmental competence of enclosed oocytes. In cultured rat COCs, treatment with recombinant BMP8A protein activated both SMAD1/5/8 and SMAD2/3 pathways; the resulting SMAD2/3 signaling induced COC expansion as well as the expression of COC expansion-related genes, whereas the resulting SMAD2/3 and SMAD1/5/8 activations were both required for protecting expanded cumulus cells from apoptosis. Taken together, our data demonstrated that addition of BMP8 protein in the in vitro rat COC cultures not only promotes cumulus expansion but also sustains survival of expanded cumulus cells via different SMAD downstreams. With these capabilities, BMP8 may have clinical applications to ameliorate the fertilizability and subsequent developmental competence of the enclosed oocytes when doing in vitro COC maturation.

Detection of selection signatures for response to Aleutian mink disease virus infection in American mink

Aleutian disease (AD) is the most significant health issue for farmed American mink. The objective of this study was to identify the genomic regions subjected to selection for response to infection with Aleutian mink disease virus (AMDV) in American mink using genotyping by sequencing (GBS) data. A total of 225 black mink were inoculated with AMDV and genotyped using a GBS assay based on the sequencing of ApeKI-digested libraries. Five AD-characterized phenotypes were used to assign animals to pairwise groups. Signatures of selection were detected using integrated measurement of fixation index (F_ST) and nucleotide diversity (θπ), that were validated by haplotype-based (hap-FLK) test. The total of 99 putatively selected regions harbouring 63 genes were detected in different groups. The gene ontology revealed numerous genes related to immune response (e.g. TRAF3IP2, WDR7, SWAP70, CBFB, and GPR65), liver development (e.g. SULF2, SRSF5) and reproduction process (e.g. FBXO5, CatSperβ, CATSPER4, and IGF2R). The hapFLK test supported two strongly selected regions that contained five candidate genes related to immune response, virus–host interaction, reproduction and liver regeneration. This study provided the first map of putative selection signals of response to AMDV infection in American mink, bringing new insights into genomic regions controlling the AD phenotypes.

Long non-coding RNA lnc-CCNL1-3:1 promotes granulosa cell apoptosis and suppresses glucose uptake in women with polycystic ovary syndrome

Polycystic ovary syndrome (PCOS) is a common endocrine and metabolic disease in premenopausal women. Long non-coding RNAs (lncRNAs) constitute important factors in numerous biological processes. However, their roles in PCOS pathogenesis require further clarification. Our study aims to elucidate the roles of lncRNA lnc-CCNL1-3:1 (CCNL) in PCOS. CCNL expression in human luteinized granulosa cells (hLGCs) derived from women with and without PCOS was detected. The full length of CCNL was obtained by 5' and 3' rapid amplification of cDNA ends. CCNL roles in granulosa cell apoptosis, mitochondrial function, and glucose uptake were evaluated. The binding relationship between CCNL and forkhead box O1 (FOXO1) was determined by RPISeq, RNA immunoprecipitation, subcellular fractionation, and immunofluorescence. In KGN cells and hLGCs, CCNL overexpression upregulated FOXO1 expression, promoted cell apoptosis, reduced glucose transport capability, and impaired mitochondrial function, and these effects were partially abolished by silencing FOXO1. The interaction of CCNL with FOXO1 might prevents FOXO1 exclusion from the nucleus and subsequent degradation in the cytosol. We determined that CCNL serve as a facilitator in the processes of PCOS. CCNL might participate in PCOS pathologies such as follicular atresia and insulin resistance.

Transcriptomics of cumulus cells - a window into oocyte maturation in humans

Background

Cumulus cells (CC) encapsulate growing oocytes and support their growth and development. Transcriptomic signatures of CC have the potential to serve as valuable non-invasive biomarkers for oocyte competency and potential. The present sibling cumulus-oocyte-complex (COC) cohort study aimed at defining functional variations between oocytes of different maturity exposed to the same stimulation conditions, by assessing the transcriptomic signatures of their corresponding CC. CC were collected from 18 patients with both germinal vesicle and metaphase II oocytes from the same cycle to keep the biological variability between samples to a minimum. RNA sequencing, differential expression, pathway analysis, and leading-edge were performed to highlight functional differences between CC encapsulating oocytes of different maturity.

Results

Transcriptomic signatures representing CC encapsulating oocytes of different maturity clustered separately on principal component analysis with 1818 genes differentially expressed. CCs encapsulating mature oocytes were more transcriptionally synchronized when compared with CCs encapsulating immature oocytes. Moreover, the transcriptional activity was lower, albeit not absent, in CC encapsulating mature oocytes, with 2407 fewer transcripts detected than in CC encapsulating immature (germinal vesicle - GV) oocytes. Hallmark pathways and ovarian processes that were affected by oocyte maturity included cell cycle regulation, steroid metabolism, apoptosis, extracellular matrix remodeling, and inflammation.

Conclusions

Herein we review our findings and discuss how they align with previous literature addressing transcriptomic signatures of oocyte maturation. Our findings support the available literature and enhance it with several genes and pathways, which have not been previously implicated in promoting human oocyte maturation. This study lays the ground for future functional studies that can enhance our understanding of human oocyte maturation.

Long non-coding RNAs in ovarian granulosa cells

Granulosa cells (GCs) are somatic cells surrounding oocytes within follicles and are essential for folliculogenesis. Pathological changes in GCs are found in several ovarian disorders. Recent reports have indicated that long non-coding RNAs (lncRNAs), which modulate gene expression via multiple mechanisms, are key regulators of the normal development of GCs, follicles, and ovaries. In addition, accumulating evidence has suggested that lncRNAs can be utilized as biomarkers for the diagnosis and prognosis of GC-related diseases, such as polycystic ovary syndrome (PCOS) and premature ovarian insufficiency (POI). Therefore, lncRNAs not only play a role in GCs that are involved in normal folliculogenesis, but they may also be considered as potential candidate biomarkers and therapeutic targets in GCs under pathological conditions. In the future, a detailed investigation of the in vivo delivery or targeting of lncRNAs and large-cohort-validation of the clinical applicability of lncRNAs is required.

Granulosa cell proliferation is inhibited by PGE2 in the primate ovulatory follicle

Prostaglandin E2 (PGE2) is a key paracrine mediator of ovulation. Few specific PGE2-regulated gene products have been identified, so we hypothesized that PGE2 may regulate the expression and/or activity of a network of proteins to promote ovulation. To test this concept, Ingenuity Pathway Analysis (IPA) was used to predict PGE2-regulated functionalities in the primate ovulatory follicle. Cynomolgus macaques underwent ovarian stimulation. Follicular granulosa cells were obtained before (0 h) or 36 h after an ovulatory dose of human chorionic gonadotropin (hCG), with ovulation anticipated 37–40 h after hCG. Granulosa cells were obtained from additional monkeys 36 h after treatment with hCG and the PTGS2 inhibitor celecoxib, which significantly reduced hCG-stimulated follicular prostaglandin synthesis. Granulosa cell RNA expression was determined by microarray and analyzed using IPA. No granulosa cell mRNAs were identified as being significantly up-regulated or down-regulated by hCG + celecoxib compared with hCG only. However, IPA predicted that prostaglandin depletion significantly regulated several functional pathways. Cell cycle/cell proliferation was selected for further study because decreased granulosa cell proliferation is known to be necessary for ovulation and formation of a fully-functional corpus luteum. Prospective in vivo and in vitro experiments confirmed the prediction that hCG-stimulated cessation of granulosa cell proliferation is mediated via PGE2. Our studies indicate that PGE2 provides critical regulation of granulosa cell proliferation through mechanisms that do not involve significant regulation of mRNA levels of key cell cycle regulators. Pathway analysis correctly predicted that PGE2 serves as a paracrine mediator of this important transition in ovarian structure and function.

The Interaction of lncRNA XLOC-2222497, AKR1C1, and Progesterone in Porcine Endometrium and Pregnancy

The endometrium is an important tissue for pregnancy and plays an important role in reproduction. In this study, high-throughput transcriptome sequencing was performed in endometrium samples of Meishan and Yorkshire pigs on days 18 and 32 of pregnancy. Aldo-keto reductase family 1 member C1 (AKR1C1) was found to be a differentially expressed gene, and was identified by quantitative real-time PCR (qRT-PCR) and Western blot. Immunohistochemistry results revealed the cellular localization of the AKR1C1 protein in the endometrium. Luciferase activity assay demonstrated that the AKR1C1 core promoter region was located in the region from -706 to -564, containing two nuclear factor erythroid 2-related factor 2 (NRF2) binding sites (antioxidant response elements, AREs). XLOC-2222497 was identified as a nuclear long non-coding RNA (lncRNA) highly expressed in the endometrium. XLOC-2222497 overexpression and knockdown have an effect on the expression of AKR1C1. Endocrinologic measurement showed the difference in progesterone levels between Meishan and Yorkshire pigs. Progesterone treatment upregulated AKR1C1 and XLOC-2222497 expression in porcine endometrial epithelial cells. In conclusion, transcriptome analysis revealed differentially expressed transcripts during the early pregnancy process. Further experiments demonstrated the interaction of XLOC-2222497/AKR1C1/progesterone in the endometrium and provided new potential targets for pregnancy maintenance and its control.

Comparative Transcriptome Analysis Suggests Key Roles for 5-Hydroxytryptamlne Receptors in Control of Goose Egg Production

To date, research on poultry egg production performance has only been conducted within inter or intra-breed groups, while those combining both inter- and intra-breed groups are lacking. Egg production performance is known to differ markedly between Sichuan white goose (Anser cygnoides) and Landes goose (Anser anser). In order to understand the mechanism of egg production performance in geese, we undertook this study. Here, 18 ovarian stromal samples from both Sichuan white goose and Landes goose at the age of 145 days (3 individuals before egg production initiation for each breed) and 730 days (3 high- and low egg production individuals during non-laying periods for each breed) were collected to reveal the genome-wide expression profiles of ovarian mRNAs and lncRNAs between these two geese breeds at different physiological stages. Briefly, 58, 347, 797, 777, and 881 differentially expressed genes (DEGs) and 56, 24, 154, 105, and 224 differentially expressed long non-coding RNAs (DElncRNAs) were found in LLD vs. HLD (low egg production Landes goose vs. high egg production Landes goose), LSC vs. HSC (low egg production Sichuan White goose vs. high egg production Sichuan white goose), YLD vs. YSC (young Landes goose vs. young Sichuan white goose), HLD vs. HSC (high egg production Landes goose vs. high egg production Sichuan white goose), and LLD vs. LSC (low egg production Landes goose vs. low egg production Sichuan white goose) groups, respectively. Functional enrichment analysis of these DEGs and DElncRNAs suggest that the "neuroactive ligand-receptor interaction pathway" is crucial for egg production, and particularly, members of the 5-hydroxytryptamine receptor (HTR) family affect egg production by regulating ovarian metabolic function. Furthermore, the big differences in the secondary structures among HTR1F and HTR1B, HTR2B, and HTR7 may lead to their different expression patterns in goose ovaries of both inter- and intra-breed groups. These results provide novel insights into the mechanisms regulating poultry egg production performance.

Comprehensive analysis of miRNAs, lncRNAs, and mRNAs reveals potential players of sexually dimorphic and left-right asymmetry in chicken gonad during gonadal differentiation

Despite thousands of sex-biased genes being found in chickens, the genetic control of sexually dimorphic and left-right asymmetry during gonadal differentiation is not yet completely understood. This study aimed to identify microRNAs (miRNAs), long noncoding RNAs (lncRNAs), messenger RNAs (mRNAs), and signaling pathways during gonadal differentiation in chick embryos (day 6/stage 29). The left and right gonads were collected for RNA sequencing. Sex-biased, side-biased miRNAs, lncRNAs, mRNAs, and shared differentially expressed miRNAs (DEmiRNA)–differentially expressed mRNAs (DEmRNA)–differentially expressed lncRNAs (DElncRNA) interaction networks were performed. A total of 8 DEmiRNAs, 183 DElncRNAs, and 123 DEmRNAs were identified for the sex-biased genes, and 7 DEmiRNAs, 189 DElncRNAs, and 183 DEmRNAs for the side-biased genes. The results of quantitative real-time PCR were generally consistent with the RNA-sequencing results. The study suggested that miRNAs and lncRNAs regulation were novel gene-specific dosage compensation mechanism and they could contribute to left-right asymmetry of chicken, but sex-biased and side-biased miRNAs, lncRNAs, and mRNAs were independent of each other. The competing endogenous RNA (ceRNA) networks showed that 17 target pairs including miR-7b (CYP19A1, FSHR, GREB1, STK31, CORIN, and TDRD9), miR-211 (FSHR, GREB1, STK31, CORIN, and TDRD9), miR-204 (FSHR, GREB1, CORIN, and TDRD9), and miR-302b-5p (CYP19A1 and TDRD9) may play crucial roles in ovarian development. These analyses provide new clues to uncover molecular mechanisms and signaling networks of ovarian development.

Elucidating Decorin's role in the preovulatory follicle

Background

DCN (decorin) is a proteoglycan known to be involved in regulating cell proliferation, collagen fibril organization and migration. In our global transcriptome RNA-sequencing approach to systematically identify new ovulation-associated genes, DCN was identified as one of the highly regulated genes. We therefore hypothesize that DCN may have a role in ovulatory processes such as cell migration and proliferation.

Aim

To characterize the expression, regulation and function of the proteoglycan DCN in the human ovarian follicles during the preovulatory period.

Methods

The in-vivo expression of DCN mRNA in mural (MGCs) and cumulus (CGCs) granulosa cells was characterized using quantitative RT-PCR and western blot.

A signaling study was performed by treating human MGCs cultures with gonadotropins and different stimulators and inhibitors to determine their effect on DCN expression by qRT- PCR and elucidate the pathways regulating these proteins.

In a functional study, KGN granulosa cell line was used to study cell migration with a scratch assay.

Results

DCN mRNA expression was significantly higher in MGCs compared to CGCs. DCN mRNA was significantly higher in CGCs surrounding mature metaphase II (MII) oocytes compared to CGCs of germinal vesicle (GV) and metaphase I (MI) oocytes.

hCG significantly increased DCN mRNA and protein expression levels in cultured MGCs. Using signal transduction activators and inhibitors, we demonstrated that DCN induction by LH/hCG is carried out via PKA, PKC, ERK/MEK, and PI3K pathways.

We showed that DCN expression is also induced in high-density cell cultures, in a dose-dependent pattern. In addition, progesterone induced a significant increase in DCN secretion to the media. MGCs from follicles of endometriosis patients exhibited reduced (about 20% of) mRNA transcriptions levels compared to MGCs follicles of control patients. More significantly, we found that DCN has an inhibiting effect on KGN cell migration.

Conclusions

Our study indicates that DCN is a unique ovulatory gene.

Our findings support the hypothesis that DCN plays an important new role during the preovulatory period and ovulation, and stress its involvement in endometriosis infertility. A better understanding of DCN role in ovulation and endometriosis may provide treatment for some types of infertility.

Genetic Parameters for Yolk Cholesterol and Transcriptional Evidence Indicate a Role of Lipoprotein Lipase in the Cholesterol Metabolism of the Chinese Wenchang Chicken

The yolk cholesterol has been reported to affect egg quality and breeding performance in chickens. However, the genetic parameters and molecular mechanisms regulating yolk cholesterol remain largely unknown. Here, we used the Wenchang chicken, a Chinese indigenous breed with a complete pedigree, as an experimental model, and we examined 24 sire families (24 males and 240 females) and their 362 daughters. First, egg quality and yolk cholesterol content were determined in 40-week-old chickens of two consecutive generations, and the heritability of these parameters was analyzed using the half-sib correlation method. Among first-generation individuals, the egg weight, egg shape index, shell strength, shell thickness, yolk weight, egg white height, Haugh unit, and cholesterol content were 45.36 ± 4.44 g, 0.81 ± 0.12, 3.07 ± 0.92 kg/cm², 0.340 ± 0.032 mm, 15.57 ± 1.64 g, 3.36 ± 1.15 mm, 58.70 ± 12.33, and 274.3 ± 36.73 mg/egg, respectively. When these indexes were compared to those of the following generation, no statistically significant difference was detected. Although yolk cholesterol content was not associated with egg quality in females, an increase in yolk cholesterol content was correlated with increased yolk weight and albumin height in sire families (p < 0.05). Moreover, the heritability estimates for the yolk cholesterol content were 0.328 and 0.530 in female and sire families, respectively. Therefore, the yolk cholesterol content was more strongly associated with the sire family. Next, chickens with low and high yolk cholesterol contents were selected for follicular membrane collection. Total RNA was extracted from these samples and used as a template for transcriptional sequencing. In total, 375 down- and 578 upregulated genes were identified by comparing the RNA sequencing data of chickens with high and low yolk cholesterol contents. Furthermore, Gene Ontology term and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses indicated the involvement of energy metabolism and immune-related pathways in yolk cholesterol deposition. Several genes participating in the regulation of the yolk cholesterol content were located on the sex chromosome Z, among which lipoprotein lipase (LPL) was associated with the peroxisome proliferator-activated receptor signaling pathway and the Gene Ontology term cellular component. Collectively, our data suggested that the ovarian steroidogenesis pathway and the downregulation of LPL played critical roles in the regulation of yolk cholesterol content.

Analysis of the equine "cumulome" reveals major metabolic aberrations after maturation in vitro

Background

Maturation of oocytes under in vitro conditions (IVM) results in impaired developmental competence compared to oocytes matured in vivo. As oocytes are closely coupled to their cumulus complex, elucidating aberrations in cumulus metabolism in vitro is important to bridge the gap towards more physiological maturation conditions. The aim of this study was to analyze the equine “cumulome” in a novel combination of proteomic (nano-HPLC MS/MS) and metabolomic (UPLC-nanoESI-MS) profiling of single cumulus complexes of metaphase II oocytes matured either in vivo (n = 8) or in vitro (n = 7).

Results

A total of 1811 quantifiable proteins and 906 metabolic compounds were identified. The proteome contained 216 differentially expressed proteins (p ≤ 0.05; FC ≥ 2; 95 decreased and 121 increased in vitro), and the metabolome contained 108 metabolites with significantly different abundance (p ≤ 0.05; FC ≥ 2; 24 decreased and 84 increased in vitro). The in vitro “cumulome” was summarized in the following 10 metabolic groups (containing 78 proteins and 21 metabolites): (1) oxygen supply, (2) glucose metabolism, (3) fatty acid metabolism, (4) oxidative phosphorylation, (5) amino acid metabolism, (6) purine and pyrimidine metabolism, (7) steroid metabolism, (8) extracellular matrix, (9) complement cascade and (10) coagulation cascade. The KEGG pathway “complement and coagulation cascades” (ID4610; n = 21) was significantly overrepresented after in vitro maturation. The findings indicate that the in vitro condition especially affects central metabolism and extracellular matrix composition. Important candidates for the metabolic group oxygen supply were underrepresented after maturation in vitro. Additionally, a shift towards glycolysis was detected in glucose metabolism. Therefore, under in vitro conditions, cumulus cells seem to preferentially consume excess available glucose to meet their energy requirements. Proteins involved in biosynthetic processes for fatty acids, cholesterol, amino acids, and purines exhibited higher abundances after maturation in vitro.

Conclusion

This study revealed the marked impact of maturation conditions on the “cumulome” of individual cumulus oocyte complexes. Under the studied in vitro milieu, cumulus cells seem to compensate for a lack of important substrates by shifting to aerobic glycolysis. These findings will help to adapt culture media towards more physiological conditions for oocyte maturation.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5836-5) contains supplementary material, which is available to authorized users.

HAS2-AS1 is a novel LH/hCG target gene regulating HAS2 expression and enhancing cumulus cells migration

Background

The cumulus expansion process is one of the LH mediated ovulatory processes. Hyaluronan synthase 2 (HAS2) regulates the synthesis of hyaluronic acid, the main component of the cumulus expansion process. Recently, the lncRNA HAS2 antisense RNA 1 (HAS2-AS1) was identified in our global transcriptome RNA-sequencing of novel ovulation associated genes. The role of HAS2-AS1 in HAS2 regulation w.as studied previously with contradictive results in different models but not in the ovary. Taken together the induction of HAS2-AS1 and the important role of HAS2 in the cumulus expansion process, we hypothesize that HAS2-AS1 regulate HAS2 expression and function in the ovary. Therefore we undertook to study the expression, regulation, and possible functional role of HAS2-AS1 in the human ovary.

Results

HAS2-AS1, located within the HAS2 gene that was highly regulated in our library. We found that HAS2-AS1 express mainly in cumulus cells (CCs). Furthermore, HAS2-AS1 showed low expression in immature CCs and a significant increase expression in mature CCs. Functional studies reveal that inhibition of HAS2-AS1 by siRNA caused decrease expression of HAS2. Furthermore, inhibition of HAS2-AS1 by siRNA results in decrease migration of granulosa cells.

Conclusions

Our results suggest that HAS2-AS1 is an LH/hCG target gene that plays a positive role in HAS2 expression and thus might play a role in regulating cumulus expansion and migration.

Long noncoding RNA and mRNA expression profiles following igf3 knockdown in common carp, Cyprinus carpio

As a novel IGF system member, igf3 plays an important role in gonadal development of teleost fish. Although studies have reported the unusual expression of igf3 in fish gonad, whether the igf3 affects the expression of long noncoding RNAs (lncRNAs) in gonad remains unknown. In this study, an igf3 knockdown common carp (Cyprinus carpio) model was established by RNA interference. Then RNA sequencing of C. carpio gonad after igf3 knockdown was performed. A total of 327,169,410 and 306,305,018 clean reads were identified from control and igf3-dsRNA interference group, respectively. After a stringent filtering, RNA-seq yielded 14199 lncRNA and 106932 mRNA transcripts with 124 and 353 differentially expressed lncRNAs and mRNAs. Our dataset provides an extensive resource for understanding the potential regulatory molecular mechanism of igf3 in early stage of gonadal development in C. carpio.

Filtered reproductive long non-coding RNAs by genome-wide analyses of goat ovary at different estrus periods

Background

The goat is an important farm animal. Reproduction is an important process of goat farming. The ovary is the most important reproductive organ for goats. In recent years, an increasing number of long non-coding RNAs (lncRNAs) have been implicated in the regulation of mammal reproduction. However, there are few studies on the function of lncRNAs in reproduction, particularly lncRNAs in the ovary.

Results

The sequencing of goat ovaries generated 1,122,014,112 clean reads, and 4926 lncRNAs and 1454 TUCPs (transcripts of uncertain coding potential) were identified for further analysis by using the coding potential analysis software, CNCI, CPC and Pfam-sca. There were 115 /22 differential lncRNAs /TUCPs transcripts between the ovaries of the luteal phase and the follicular phase. We predicted the related genes of lncRNA /TUCP based on co-expression and co-localization methods. In total, 2584 /904 genes were predicted by co-expression, and 326/73 genes were predicted by co-localization. The functions of these genes were further analyzed with GO and KEGG analysis. The results showed that lncRNAs /TUCPs, which are highly expressed in goat ovaries in the luteal phase, are mainly associated with the synthesis of progesterone, and we filtered the lncRNAs /TUCPs, such as XR_001918177.1 and TUCP_001362, which may regulate the synthesis of progesterone; lncRNAs /TUCPs, which are highly expressed in goat ovaries in the follicular phase, are mainly associated with oogenesis and the maturation of oocytes, and we filtered the lncRNAs /TUCPs that may regulate the oogenesis and maturation of oocyte, such as XR_001917388.1 and TUCP_000849.

Conclusion

The present study provided the genome expression profile of lncRNAs /TUCPs in goat ovaries at different estrus periods and filtered the potential lncRNAs /TUCPs associated with goat reproduction. These results are helpful to further study the molecular mechanisms of goat reproduction.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5268-7) contains supplementary material, which is available to authorized users.

Characterization of the miRNA regulators of the human ovulatory cascade

Ovarian follicular development and ovulation are complex and tightly regulated processes that involve regulation by microRNAs (miRNAs). We previously identified differentially expressed mRNAs between human cumulus granulosa cells (CGCs) from immature early antral follicles (germinal vesicle - GV) and mature preovulatory follicles (metaphase II - M2). In this study, we performed an integrated analysis of the transcriptome and miRNome in CGCs obtained from the GV cumulus-oocyte complex (COC) obtained from IVM and M2 COC obtained from IVF. A total of 43 differentially expressed miRNAs were identified. Using Ingenuity IPA analysis, we identified 7288 potential miRNA-regulated target genes. Two hundred thirty-four of these target genes were also found in our previously generated ovulatory gene library while exhibiting anti-correlated expression to the identified miRNAs. IPA pathway analysis suggested that miR-21 and FOXM1 cooperatively inhibit CDC25A, TOP2A and PRC1. We identified a mechanism for the temporary inhibition of VEGF during ovulation by TGFB1, miR-16-5p and miR-34a-5p. The linkage bioinformatics analysis between the libraries of the coding genes from our preliminary study with the newly generated library of regulatory miRNAs provides us a comprehensive, integrated overview of the miRNA-mRNA co-regulatory networks that may play a key role in controlling post-transcriptomic regulation of the ovulatory process.

Cumulus cells surrounding oocytes with high developmental competence exhibit down-regulation of phosphoinositol 1,3 kinase/protein kinase B (PI3K/AKT) signalling genes involved in proliferation and survival

STUDY QUESTION

Is the phosphoinositol 1,3-kinase/protein kinase B (PI3K/AKT) pathway expression profile in cumulus cells (CCs) a potential marker of oocyte competence and predictive of pregnancy outcome?

SUMMARY ANSWER

Eleven genes (AKT1, ARHGEF7, BCL2L1, CCND1, E2F1, HRAS, KCNH2, PIK3C2A, SHC1, SOS1 and SPP1) in the PI3K/AKT pathway were significantly down-regulated in CCs from oocytes that went on to produce a pregnancy compared to CCs associated with a negative outcome.

WHAT IS KNOWN ALREADY

The PI3K/AKT pathway plays a pivotal role in the interdependence and continuous feedback between the oocyte and CCs.

STUDY DESIGN SIZE, DURATION

The expression analysis of 92 transcripts in the PI3K/AKT pathway in CCs from patients with negative or positive pregnancy outcome, after single embryo transfer, was performed. Mouse CCs target gene expression was conducted to associate the expression profile of PI3K/AKT pathway to oocyte developmental profile.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Fifty-five good prognosis IVF patients who had been referred to IVF or intracytoplasmic sperm injection treatment for male-factor infertility or tubal disease were enroled. CCs from single cumulus-oocyte complexes (COCs) from 16 patients who underwent a single embryo transfer were analyzed. Twenty-five CD-1 mice were used to assess gene expression in CCs associated with oocytes with different competence in relation to hCG priming. A total 220 human COCs were collected. The RNA extracted from CCs of 16 selected patients was used to analyze PI3K/AKT pathway gene expression employing a 96-well custom TaqMan Array. Expression data of CCs associated to positive IVF outcome were compared to data from negative outcome samples. Mice were sacrificed after 9, 12, 15, 21 and 24 h post-hCG administration to obtain CCs from MII oocytes with different developmental competence. Akt1, Bcl2l2 and Shc1 expression were tested in the collected mouse CCs. In addition, the expression of upstream regulator ESR1, the gene encoding for the oestrogen receptor ERβ, and the downstream effectors of the pathway FOXO1, FOXO3 and FOXO4 was evaluated in human and mouse samples.

MAIN RESULTS AND THE ROLE OF CHANCE

Transcripts involved in the PI3K Signaling Pathway were selectively modulated according to the IVF/ICSI outcome of the oocyte. Eleven transcripts in this pathway were significantly down-regulated in all samples of CCs from oocytes with positive when compared those with a negative outcome. These outcomes were confirmed in mouse CCs associated with oocytes at different maturation stages. Expression data revealed that the down-regulation of ESR1 could be related to oocyte competence and is likely to be the driver of expression changes highlighted in the PI3K/AKT pathway.

LIMITATIONS REASONS FOR CAUTION

Small sample size and retrospective design.

WIDER IMPLICATIONS OF THE FINDINGS

The CCs expression profile of PI3K/AKT signaling genes, disclosed a specific CCs gene signature related to oocyte competence. It could be speculated that CCs associated with competent oocytes have completed their role in sustaining oocyte development and are influencing their fate in response to metabolic and hormonal changes by de-activating anti-apoptotic signals.

STUDY FUNDING/COMPETING INTEREST(S)

Supported by Merck Serono an affiliate of Merck KGaA, Darmstadt, Germany (research grant for the laboratory session; Merck KGaA reviewed the manuscript for medical accuracy only before journal submission. The authors are fully responsible for the content of this manuscript, and the views and opinions described in the publication reflect solely those of the authors). The authors declare no conflict of interest.

Differential long non-coding RNA expression profiles in human oocytes and cumulus cells

Progress in assisted reproductive technologies strongly relies on understanding the regulation of the dialogue between oocyte and cumulus cells (CCs). Little is known about the role of long non-coding RNAs (lncRNAs) in the human cumulus-oocyte complex (COC). To this aim, publicly available RNA-sequencing data were analyzed to identify lncRNAs that were abundant in metaphase II (MII) oocytes (BCAR4, C3orf56, TUNAR, OOEP-AS1, CASC18, and LINC01118) and CCs (NEAT1, MALAT1, ANXA2P2, MEG3, IL6STP1, and VIM-AS1). These data were validated by RT-qPCR analysis using independent oocytes and CC samples. The functions of the identified lncRNAs were then predicted by constructing lncRNA-mRNA co-expression networks. This analysis suggested that MII oocyte lncRNAs could be involved in chromatin remodeling, cell pluripotency and in driving early embryonic development. CC lncRNAs were co-expressed with genes involved in apoptosis and extracellular matrix-related functions. A bioinformatic analysis of RNA-sequencing data to identify CC lncRNAs that are affected by maternal age showed that lncRNAs with age-related altered expression in CCs are essential for oocyte growth. This comprehensive analysis of lncRNAs expressed in human MII oocytes and CCs could provide biomarkers of oocyte quality for the development of non-invasive tests to identify embryos with high developmental potential.

Advances of Long Noncoding RNAs-mediated Regulation in Reproduction

OBJECTIVE

Advances in genomics and molecular biology have led to the discovery of a large group of uncharacterized long noncoding RNAs (lncRNAs). Emerging evidence indicated that many lncRNAs function in multiple biological processes and its dysregulation often causes diseases. Recent studies suggested that almost all regulatory lncRNAs interact with biological macromolecules such as DNA, RNA, and protein. LncRNAs regulate gene expression mainly on three levels, including epigenetic modification, transcription, and posttranscription, through DNA methylation, histone modification, and chromatin remodeling. LncRNAs can also affect the development of diseases and therefore be used to diagnose and treat diseases. With new sequencing and microarray techniques, hundreds of lncRNAs involved in reproductive disorders have been identified, but their functions in these disorders are undefined.

DATA SOURCES

This review was based on articles published in PubMed databases up to July 10, 2017, with the following keywords: "long noncoding RNAs", "LncRNA", "placentation", and "reproductive diseases".

STUDY SELECTION

Original articles and reviews on the topics were selected.

RESULTS

LncRNAs widely participate in various physiological and pathological processes as a new class of important regulatory factors. In spermatogenesis, spermatocytes divide and differentiate into mature spermatozoa. The whole process is elaborately regulated by the expression of phase-specific genes that involve many strains of lncRNAs. Literature showed that lncRNA in reproductive cumulus cells may contribute to the regulation of oocyte maturation, fertilization, and embryo development.

CONCLUSIONS

LncRNA has been found to play a role in the development of reproduction. Meanwhile, we reviewed the studies on how lncRNAs participate in reproductive disorders, which provides a basis for the study of lncRNA in reproduction regulation.

Long Noncoding RNAs: Potential Regulators Involved in the Pathogenesis of Polycystic Ovary Syndrome

Polycystic ovary syndrome (PCOS) is the most common cause of anovulatory infertility in women of reproductive age, and its etiology remains poorly understood. Altered activities of long noncoding RNAs (lncRNAs) have been associated with human diseases and development. However, the roles of lncRNAs are unknown in reproductive medicine. We investigated the potential role of lncRNAs in the pathogenesis of PCOS, using human granulosa cells (GCs) and the KGN cell line. We used microarrays to compare lncRNA expression profiles in GCs from seven patients with PCOS and seven matched women. GC samples were collected during 2014 to 2016 from infertile women in Guangzhou, China. Quantitative real-time polymerase chain reaction was used to measure levels of the lncRNA HCG26 in GCs from 53 patients with PCOS and 50 controls. HCG26 was knocked down with locked nucleic acid GapmeRs in KGN cells to examine its role in cell proliferation, aromatase and follicle-stimulating hormone receptor gene expression, and estradiol production. A total of 862 lncRNA transcripts and 998 messenger RNA transcripts were differentially expressed (greater than or equal to twofold change; P < 0.05) in PCOS GCs compared with those of controls. HCG26 levels were upregulated in patients with PCOS and were associated with antral follicle count. HCG26 knockdown in KGN cells inhibited cell proliferation and cell-cycle progression and increased aromatase gene expression and estradiol production. Our study reports the lncRNA profiles in GCs from patients who have PCOS and those from healthy women and suggests that dysregulated lncRNAs may play vital roles in GC proliferation and steroidogenesis, providing insights into the pathogenesis of PCOS.

Obesity Modulates Inflammation and Lipid Metabolism Oocyte Gene Expression: A Single-Cell Transcriptome Perspective

CONTEXT

It is hypothesized that obesity adversely affects the ovarian environment, which can disrupt oocyte maturation and embryonic development.

OBJECTIVE

This study aimed to compare oocyte gene expression profiles and follicular fluid (FF) content from overweight/obese (OW) women and normal-weight (NW) women who were undergoing fertility treatments.

DESIGN

Using single-cell transcriptomic analyses, we investigated oocyte gene expression using RNA sequencing.

PATIENTS OR OTHER PARTICIPANTS

Eleven OW women and 13 NW women undergoing fertility treatments were enrolled.

MAIN OUTCOME MEASURES

Oocyte messenger RNA profiles as well as serum and FF hormone and lipid levels were assessed.

RESULTS

OW women had significantly higher body mass index, body fat percentage, and serum homeostatic model assessment-insulin resistance index compared with NW women (P < 0.01). Serum leptin and C-reactive protein (CRP) levels as well as FF leptin, CRP, and triglyceride levels were increased (P < 0.05) in OW compared with NW women. Oocytes from OW women had increased expression of proinflammatory (CXCL2; P = 0.071) and oxidative stress-related (DUSP1; P = 0.051) genes but had decreased expression of GAS7 (fat metabolism; P = 0.065), TXNIP (oxidative stress; P = 0.055), and transcription factors ID3 (P = 0.075) and TWIST1 (P = 0.099) compared with NW women.

CONCLUSIONS

These findings provide evidence for the significant influence of body composition on oocyte transcript abundance in women undergoing hormonal induction to retrieve oocytes. They further identify the potential for maternal diet to influence oocyte gene expression. The preconception period is, therefore, an important window of opportunity to consider for lifestyle interventions.

Karyomapping and how is it improving preimplantation genetics?

INTRODUCTION

Preimplantation genetic diagnosis and screening (PGD/PGS) has been applied clinically for >25 years however inherent drawbacks include the necessity to tailor each case to the trait in question, and that technology to detect monogenic and chromosomal disorders respectively is fundamentally different. Areas covered: The area of preimplantation genetics has evolved over the last 25 years, adapting to changes in technology and the need for more efficient, streamlined diagnoses. Karyomapping allows the determination of inheritance from the (grand)parental haplobocks through assembly of inherited chromosomal segments. The output displays homologous chromosomes, crossovers and the genetic status of the embryos by linkage comparison, as well as chromosomal disorders. It also allows for determination of heterozygous SNP calls, avoiding the risks of allele dropout, a common problem with other PGD techniques. Manuscripts documenting the evolution of preimplantation genetics, especially those investigating technologies that would simultaneously detect monogenic and chromosomal disorders, were selected for review. Expert commentary: Karyomapping is currently available for detection of single gene disorders; ~1000 clinics worldwide offer it (via ~20 diagnostic laboratories) and ~2500 cases have been performed. Due an inability to detect post-zygotic trisomy reliably however and confounding problems of embryo mosaicism, karyomapping has yet to be applied clinically for detection of chromosome disorders.

Non-coding RNAs in the Ovarian Follicle

The mammalian ovarian follicle is the complex reproductive unit comprising germ cell, somatic cells (Cumulus and Granulosa cells), and follicular fluid (FF): paracrine communication among the different cell types through FF ensures the development of a mature oocyte ready for fertilization. This paper is focused on non-coding RNAs in ovarian follicles and their predicted role in the pathways involved in oocyte growth and maturation. We determined the expression profiles of microRNAs in human oocytes and FF by high-throughput analysis and identified 267 microRNAs in FF and 176 in oocytes. Most of these were FF microRNAs, while 9 were oocyte specific. By bioinformatic analysis, independently performed on FF and oocyte microRNAs, we identified the most significant Biological Processes and the pathways regulated by their validated targets. We found many pathways shared between the two compartments and some specific for oocyte microRNAs. Moreover, we found 41 long non-coding RNAs able to interact with oocyte microRNAs and potentially involved in the regulation of folliculogenesis. These data are important in basic reproductive research and could also be useful for clinical applications. In fact, the characterization of non-coding RNAs in ovarian follicles could improve reproductive disease diagnosis, provide biomarkers of oocyte quality in Assisted Reproductive Treatment, and allow the development of therapies for infertility disorders.

Patent highlights June-July 2016

A snapshot of noteworthy recent developments in the patent literature of relevance to pharmaceutical and medical research and development.

Long non-coding RNAs in human early embryonic development and their potential in ART

BACKGROUND

Human long non-coding RNAs (lncRNAs) are an emerging category of transcripts with increasingly documented functional roles during development. LncRNAs and roles during human early embryo development have recently begun to be unravelled.

OBJECTIVE AND RATIONALE

This review summarizes the most recent knowledge on lncRNAs and focuses on their expression patterns and role during early human embryo development and in pluripotent stem cells (PSCs). Public mRNA sequencing (mRNA-seq) data were used to illustrate these expression signatures.

SEARCH METHODS

The PubMed and EMBASE databases were first interrogated using specific terms, such as 'lncRNAs', to get an extensive overview on lncRNAs up to February 2016, and then using 'human lncRNAs' and 'embryo', 'development', or 'PSCs' to focus on lncRNAs involved in human embryo development or in PSC.Recently published RNA-seq data from human oocytes and pre-implantation embryos (including single-cell data), PSC and a panel of normal and malignant adult tissues were used to describe the specific expression patterns of some lncRNAs in early human embryos.

OUTCOMES

The existence and the crucial role of lncRNAs in many important biological phenomena in each branch of the life tree are now well documented. The number of identified lncRNAs is rapidly increasing and has already outnumbered that of protein-coding genes. Unlike small non-coding RNAs, a variety of mechanisms of action have been proposed for lncRNAs. The functional role of lncRNAs has been demonstrated in many biological and developmental processes, including cell pluripotency induction, X-inactivation or gene imprinting. Analysis of RNA-seq data highlights that lncRNA abundance changes significantly during human early embryonic development. This suggests that lncRNAs could represent candidate biomarkers for developing non-invasive tests for oocyte or embryo quality. Finally, some of these lncRNAs are also expressed in human cancer tissues, suggesting that reactivation of an embryonic lncRNA program may contribute to human malignancies.

WIDER IMPLICATIONS

LncRNAs are emerging potential key players in gene expression regulation. Analysis of RNA-seq data from human pre-implantation embryos identified lncRNA signatures that are specific to this critical step. We anticipate that further studies will show that these new transcripts are major regulators of embryo development. These findings might also be used to develop new tests/treatments for improving the pregnancy success rate in IVF procedures or for regenerative medicine applications involving PSC.

Identification of molecular markers for oocyte competence in bovine cumulus cells

Cumulus cells (CCs) have an important role during oocyte growth, competence acquisition, maturation, ovulation and fertilization. In an attempt to isolate potential biomarkers for bovine in vitro fertilization, we identified genes differentially expressed in bovine CCs from oocytes with different competence statuses, through microarray analysis. The model of follicle size, in which competent cumulus-oocyte complexes (COCs) were recovered from bigger follicles (≥8.0 mm in diameter) and less competent ones from smaller follicles (1-3 mm), was used. We identified 4178 genes that were differentially expressed (P < 0.05) in the two categories of CCs. The list was further enriched, through the use of a 2.5-fold change in gene expression as a cutoff value, to include 143 up-regulated and 80 down-regulated genes in CCs of competent COCs compared to incompetent COCs. These genes were screened according to their cellular roles, most of which were related to cell cycle, DNA repair, energy metabolism, metabolism of amino acids, cell signaling, meiosis, ovulation and inflammation. Three candidate genes up-regulated (FGF11, IGFBP4, SPRY1) and three down-regulated (ARHGAP22, COL18A1 and GPC4) in CCs from COCs of big follicles (≥8.1 mm) were selected for qPCR analysis. The selected genes showed the same expression patterns by qPCR and microarray analysis. These genes may be potential genetic markers that predict oocyte competence in in vitro fertilization routines.

Competence Classification of Cumulus and Granulosa Cell Transcriptome in Embryos Matched by Morphology and Female Age

Objective

By focussing on differences in the mural granulosa cell (MGC) and cumulus cell (CC) transcriptomes from follicles resulting in competent (live birth) and non-competent (no pregnancy) oocytes the study aims on defining a competence classifier expression profile in the two cellular compartments. Design: A case-control study. Setting: University based facilities for clinical services and research. Patients: MGC and CC samples from 60 women undergoing IVF treatment following the long GnRH-agonist protocol were collected. Samples from 16 oocytes where live birth was achieved and 16 age- and embryo morphology matched incompetent oocytes were included in the study.

Methods

MGC and CC were isolated immediately after oocyte retrieval. From the 16 competent and non-competent follicles, mRNA was extracted and expression profile generated on the Human Gene 1.0 ST Affymetrix array. Live birth prediction analysis using machine learning algorithms (support vector machines) with performance estimation by leave-one-out cross validation and independent validation on an external data set.

Results

We defined a signature of 30 genes expressed in CC predictive of live birth. This live birth prediction model had an accuracy of 81%, a sensitivity of 0.83, a specificity of 0.80, a positive predictive value of 0.77, and a negative predictive value of 0.86. Receiver operating characteristic analysis found an area under the curve of 0.86, significantly greater than random chance. When applied on 3 external data sets with the end-point outcome measure of blastocyst formation, the signature resulted in 62%, 75% and 88% accuracy, respectively. The genes in the classifier are primarily connected to apoptosis and involvement in formation of extracellular matrix. We were not able to define a robust MGC classifier signature that could classify live birth with accuracy above random chance level.

Conclusion

We have developed a cumulus cell classifier, which showed a promising performance on external data. This suggests that the gene signature at least partly include genes that relates to competence in the developing blastocyst.

A maternal deletion upstream of the imprint control region 2 in 11p15 causes loss of methylation and familial Beckwith-Wiedemann syndrome

Beckwith-Wiedemann syndrome (BWS; OMIM #130650) is an overgrowth syndrome caused by different genetic or epigenetic alterations affecting imprinted regions on chromosome 11p15.5. Here we report a family with multiple offspring affected with BWS including giant omphalocoeles in which maternal transmission of a chromosomal rearrangement including an inversion and two deletions leads to hypomethylation of the imprint control region 2 (ICR2). As the deletion includes the promoter and 5' part of the KCNQ1 gene, we suggest that transcription of this gene may be involved in establishing the maternal methylation imprint of the ICR2, which is located in intron 10 of KCNQ1.

Long non-coding RNA regulation of reproduction and development

Noncoding RNAs (ncRNAs) have long been known to play vital roles in eukaryotic gene regulation. Studies conducted over a decade ago revealed that maturation of spliced, polyadenylated coding mRNA occurs by reactions involving small nuclear RNAs and small nucleolar RNAs; mRNA translation depends on activities mediated by transfer RNAs and ribosomal RNAs, subject to negative regulation by micro RNAs; transcriptional competence of sex chromosomes and some imprinted genes is regulated in cis by ncRNAs that vary by species; and both small-interfering RNAs and piwi-interacting RNAs bound to Argonaute-family proteins regulate post-translational modifications on chromatin and local gene expression states. More recently, gene-regulating noncoding RNAs have been identified, such as long intergenic and long noncoding RNAs (collectively referred to as lncRNAs)--a class totaling more than 100,000 transcripts in humans, which include some of the previously mentioned RNAs that regulate dosage compensation and imprinted gene expression. Here, we provide an overview of lncRNA activities, and then review the role of lncRNAs in processes vital to reproduction, such as germ cell specification, sex determination and gonadogenesis, sex hormone responses, meiosis, gametogenesis, placentation, non-genetic inheritance, and pathologies affecting reproductive tissues. Results from many species are presented to illustrate the evolutionarily conserved processes lncRNAs are involved in.

Transcriptomic Analysis and Meta-Analysis of Human Granulosa and Cumulus Cells

Specific gene expression in oocytes and its surrounding cumulus (CC) and granulosa (GC) cells is needed for successful folliculogenesis and oocyte maturation. The aim of the present study was to compare genome-wide gene expression and biological functions of human GC and CC. Individual GC and CC were derived from 37 women undergoing IVF procedures. Gene expression analysis was performed using microarrays, followed by a meta-analysis. Results were validated using quantitative real-time PCR. There were 6029 differentially expressed genes (q < 10⁻⁴); of which 650 genes had a log2 FC ≥ 2. After the meta-analysis there were 3156 genes differentially expressed. Among these there were genes that have previously not been reported in human somatic follicular cells, like prokineticin 2 (PROK2), higher expressed in GC, and pregnancy up-regulated nonubiquitous CaM kinase (PNCK), higher expressed in CC. Pathways like inflammatory response and angiogenesis were enriched in GC, whereas in CC, cell differentiation and multicellular organismal development were among enriched pathways. In conclusion, transcriptomes of GC and CC as well as biological functions, are distinctive for each cell subpopulation. By describing novel genes like PROK2 and PNCK, expressed in GC and CC, we upgraded the existing data on human follicular biology.

Novel contraceptive targets to inhibit ovulation: the prostaglandin E2 pathway

BACKGROUND

Prostaglandin E2 (PGE2) is an essential intrafollicular regulator of ovulation. In contrast with the one-gene, one-protein concept for synthesis of peptide signaling molecules, production and metabolism of bioactive PGE2 requires controlled expression of many proteins, correct subcellular localization of enzymes, coordinated PGE2 synthesis and metabolism, and prostaglandin transport in and out of cells to facilitate PGE2 action and degradation. Elevated intrafollicular PGE2 is required for successful ovulation, so disruption of PGE2 synthesis, metabolism or transport may yield effective contraceptive strategies.

METHODS

This review summarizes case reports and studies on ovulation inhibition in women and macaques treated with cyclooxygenase inhibitors published from 1987 to 2014. These findings are discussed in the context of studies describing levels of mRNA, protein, and activity of prostaglandin synthesis and metabolic enzymes as well as prostaglandin transporters in ovarian cells.

RESULTS

The ovulatory surge of LH regulates the expression of each component of the PGE2 synthesis-metabolism-transport pathway within the ovulatory follicle. Data from primary ovarian cells and cancer cell lines suggest that enzymes and transporters can cooperate to optimize bioactive PGE2 levels. Elevated intrafollicular PGE2 mediates key ovulatory events including cumulus expansion, follicle rupture and oocyte release. Inhibitors of the prostaglandin-endoperoxide synthase 2 (PTGS2) enzyme (also known as cyclooxygenase-2 or COX2) reduce ovulation rates in women. Studies in macaques show that PTGS2 inhibitors can reduce the rates of cumulus expansion, oocyte release, follicle rupture, oocyte nuclear maturation and fertilization. A PTGS2 inhibitor reduced pregnancy rates in breeding macaques when administered to simulate emergency contraception. However, PTGS2 inhibition did not prevent pregnancy in monkeys when administered to simulate monthly contraceptive use.

CONCLUSION

PTGS2 inhibitors alone may be suitable for use as emergency contraceptives. However, drugs of this class are unlikely to be effective as monthly contraceptives. Inhibitors of additional PGE2 synthesis enzymes or modulation of PGE2 metabolism or transport also hold potential for reducing follicular PGE2 and preventing ovulation. Approaches which target multiple components of the PGE2 synthesis-metabolism-transport pathway may be required to effectively block ovulation and lead to the development of novel contraceptive options for women. Therapies which target PGE2 may also impact disorders of the uterus and could also have benefits for women's health in addition to contraception.