Differences in cumulus cells gene expression between modified natural and stimulated in vitro fertilization cycles

Transcriptomic integrity of human oocytes used in ARTs: technical and intrinsic factor effects

BACKGROUND

Millions of children have been born throughout the world thanks to ARTs, the harmlessness of which has not yet been fully demonstrated. For years, efforts to evaluate the specific effects of ART have focused on the embryo; however, it is the oocyte quality that mainly dictates first and foremost the developmental potential of the future embryo. Ovarian stimulation, cryopreservation, and IVM are sometimes necessary steps to obtain a mature oocyte, but they could alter the appropriate expression of the oocyte genome. Additionally, it is likely that female infertility, environmental factors, and lifestyle have a significant influence on oocyte transcriptomic quality, which may interfere with the outcome of an ART attempt.

OBJECTIVE AND RATIONALE

The objective of this review is to identify transcriptomic changes in the human oocyte caused by interventions specific to ART but also intrinsic factors such as age, reproductive health issues, and lifestyle. We also provide recommendations for future good practices to be conducted when attempting ART.

SEARCH METHODS

An in-depth literature search was performed on PubMed to identify studies assessing the human oocyte transcriptome following ART interventions, or in the context of maternal aging, suboptimal lifestyle, or reproductive health issues.

OUTCOMES

ART success is susceptible to external factors, maternal aging, lifestyle factors (smoking, BMI), and infertility due to endometriosis or polycystic ovary syndrome. Indeed, all of these are likely to increase oxidative stress and alter mitochondrial processes in the foreground. Concerning ART techniques themselves, there is evidence that different ovarian stimulation regimens shape the oocyte transcriptome. The perturbation of processes related to the mitochondrion, oxidative phosphorylation, and metabolism is observed with IVM. Cryopreservation might dysregulate genes belonging to transcriptional regulation, ubiquitination, cell cycle, and oocyte growth pathways. For other ART laboratory factors such as temperature, oxygen tension, air pollution, and light, the evidence remains scarce. Focusing on genes involved in chromatin-based processes such as DNA methylation, heterochromatin modulation, histone modification, and chromatin remodeling complexes, but also genomic imprinting, we observed systematic dysregulation of such genes either after ART intervention or lifestyle exposure, as well as due to internal factors such as maternal aging and reproductive diseases. Alteration in the expression of such epigenetic regulators may be a common mechanism linked to adverse oocyte environments, explaining global transcriptomic modifications.

WIDER IMPLICATIONS

Many IVF factors and additional external factors have the potential to impair oocyte transcriptomic integrity, which might not be innocuous for the developing embryo. Fortunately, it is likely that such dysregulations can be minimized by adapting ART protocols or reducing adverse exposure.

Protein Expression and Bioinformatics Study of Granulosa Cells of Polycystic Ovary Syndrome Expressed Under the Influence of DHEA

Background

The reproductive system is heavily dependent on ovarian follicles, which are made up of germ cells (oocytes) and granulosa cells (GCs), including cumulus granulosa cells (CGCs) and mural granulosa cells (MGCs). Understanding their normal and steroid-induced functions is the key to understanding the pathophysiology of endocrinal diseases in women.

Objective

This study investigated the differentially expressed proteins by CGCs and MGCs of patients with polycystic ovarian syndrome (PCOS) and without subsequent exposure to dehydroepiandrosterone sulfate (DHEAS) and functional differentiation.

Design

The present study was observational and experimental study carried out in hospital involving 80 female patients undergoing IVF for infertility.

Methods

In this study, we isolated CGCs and MGCs from the follicular fluid of both PCOS and non-PCOS patients undergoing in vitro fertilization (IVF). The cells were cultured and treated with DHEAS for 48 hours, and these cells were extracted, digested, and analyzed by tandem mass spectrometry followed by processing of the results using open-source bioinformatics tools.

Results

The present investigation discovered 276 and 341 proteins in CGCs and MGCs, respectively. DHEAS reduced the number of proteins expressed by CGCs and MGCs to 34 and 57 from 91 and 94, respectively. Venn results of CGCs revealed 49, 53, 36, and 21 proteins in normal CGCs, PCOS-CGCs, post-DHEAS, and PCOS-CGCs, respectively. Venn analysis of MGCs showed 51 proteins specific to PCOS and 29 shared by normal and PCOS samples after DHEAS therapy. MGCs express the most binding and catalytic proteins, whereas CGCs express transporter-related proteins. A protein pathway study demonstrated considerable differences between normal and PCOS samples, while DHEAS-treated samples of both cell lines showed distinct pathways. String findings identified important network route components such as albumin, actin, apolipoprotein, complement component C3, and heat shock protein.

Conclusion

This is the first study to show how DHEAS-induced stress affects the expression of proteins by MGCs and CGCs isolated from normal and PCOS patients. Further studies are recommended to identify PCOS biomarkers from CGCs and MGCs expressed under the influence of DHEAS.

The Effect of Clomiphene Citrate and Letrozole in Apoptotic Pathways and Cell Cycle in Human Primary Cumulus Cells and the Protective Effect of Estradiol

Clomiphene citrate (CC) and letrozole are ovulatory stimulants that, despite high ovulation rates, achieve low pregnancy rates. This study aimed to investigate the in vitro effects of CC and letrozole, alone or in combination with estradiol, on apoptosis in human cumulus cells. We performed a controlled prospective study using primary cumulus cell cultures from patients undergoing in vitro fertilization (n=22). Alpha-inhibin immunocytochemistry was used to assess cell culture purity and morphology. Cell viability was evaluated by MTT assay, cell cycle status by flow cytometry, and Caspase-3, Bax and SOD-2, and S26 gene expression by qPCR. Cells were treated for 24 hours in 5 conditioned media: CC, CC + estradiol, letrozole, letrozole + estradiol and control. None of the treatments affected cell viability, but letrozole reduced the mean percentage of cells in the S phase compared to control (24.79 versus 21.70, p=0.0014). Clomiphene treatment increased mRNA expression of Bax (4 fold) and SOD-2 (2 fold), which was reversed by co-treatment with estradiol. SOD-2 expression increased in cells treated with letrozole compared to control (4 fold), which was also reversed by estradiol. These findings suggest that clomiphene citrate and letrozole do not significantly affect the viability of human cumulus cells. Still, the expression of genes involved in apoptosis was modulated by these drugs alone and in association with estradiol, suggesting that CC and letrozole may have direct effects on cumulus cells beyond their known mechanisms of action.

Comparative analyses in transcriptome of human granulosa cells and follicular fluid micro-environment between poor ovarian responders with conventional controlled ovarian or mild ovarian stimulations

BACKGROUND

Both mild and conventional controlled ovarian stimulation are the frequently used protocols for poor ovarian responders. However, there are some debates about which treatment is better. Moreover, little is known about the follicular physiology after the two ovarian stimulation protocols. This study was intended to investigate the features in granulosa cells and follicular fluid micro-environment after the two different ovarian stimulation protocols in poor responders.

METHODS

Granulosa cells RNA were sequenced using Illumina Hiseq technology. Specific differently expressed genes and proteins were verified by real-time quantitative PCR and Western blot analysis. Moreover, hormone and cytokine concentrations in the follicular fluid were measured by electrochemiluminescence immunoassay and enzyme-linked immunoabsorbent assay. The correlation between the results of molecular experiments and the laboratory outcomes were analyzed by Spearman correlation analysis.

RESULTS

The differentially expressed genes between the two groups were involved in 4 signaling pathways related to the follicular development; three proteins pertinent to the TGF-β signaling pathway were expressed differently in granulosa cells between the two, and the constituents in the follicular fluid were also different. Further, a correlation between the TGF-β signaling pathway and the good-quality embryo was observed.

CONCLUSIONS

The present study made a comparison for the first time in the transcriptome of human granulosa cells and the follicular fluid micro-environment between poor responders with the conventional controlled ovarian stimulation or the mild ovarian stimulation, showing that the TGF-β signaling pathway may correlate with the good-quality of embryos in the mild group, which may be instrumental to the choice of optimal management for IVF patients.

Stratifying Cumulus Cell Samples Based on Molecular Profiling to Help Resolve Biomarker Discrepancies and to Predict Oocyte Developmental Competence

To increase the efficiency of assisted reproductive techniques (ART), molecular studies have been performed to identify the best predictive biomarkers for selecting the most suitable germ cells for fertilization and the best embryo for intra-uterine transfer. However, across different studies, no universal markers have been found. In this study, we addressed this issue by generating gene expression and CpG methylation profiles of outer cumulus cells obtained during intra-cytoplasmic sperm injection (ICSI). We also studied the association of the generated genomic data with the clinical parameters (spindle presence, zona pellucida birefringence, pronuclear pattern, estrogen level, endometrium size and lead follicle size) and the pregnancy result. Our data highlighted the presence of several parameters that affect analysis, such as inter-individual differences, inter-treatment differences, and, above all, specific treatment protocol differences. When comparing the pregnancy outcome following the long protocol (GnRH agonist) of ovarian stimulation, we identified the single gene markers (NME6 and ASAP1, FDR < 5%) which were also correlated with endometrium size, upstream regulators (e.g., EIF2AK3, FSH, ATF4, MKNK1, and TP53) and several bio-functions related to cell death (apoptosis) and cellular growth and proliferation. In conclusion, our study highlighted the need to stratify samples that are very heterogeneous and to use pathway analysis as a more reliable and universal method for identifying markers that can predict oocyte development potential.

Ovarian stimulated cycles reduce protection of follicular fluid against free radicals

Controlled ovarian hyperstimulation cycle with exogenous gonadotropins (COH) is associated with clinical complications. The aim of this work was to determine whether COH alters the physiological antioxidant status of follicular fluid in women with no reproductive dysfunction, compared to the natural cycle (NC). In this longitudinal study, forty-one women (oocyte donors) consecutively underwent NC and COH. Follicular fluid was collected at oocyte retrieval and different redox biomarkers were determined: total antioxidant activity (TAA), oxygen radical absorbance capacity (ORAC), nitric oxide, α- and γ-tocopherol, the fatty acid composition, activities of superoxide dismutase, catalase, total and Se-dependent glutathione peroxidases, and the antioxidant paraoxonase (PON) family. Results showed that TAA (1.70 ± 0.03 mM versus 1.86 ± 0.03 mM, p < 0.05), α-tocopherol (4.37 ± 0.26 μM versus 5.74 ± 0.30 μM, p < 0.05), PON1 paraoxonase (245 ± 24 nmol/min/ml versus 272 ± 27 nmol/min/ml, p < 0.05), PON1 arylesterase (87.2 ± 4.6 μmol/min/ml versus 99.3 ± 4.8 μmol/min/ml, p < 0.05), and PON3 simvastatinase (13.48 ± 0.52 nmol/min/ml versus 16.29 ± 0.72 nmol/min/ml, p < 0.001) were significantly lower in COH versus NC. Fatty acids from COH were more saturated, increasing palmitate and decreasing the n-6 and total polyunsaturated fatty acids (PUFAs). Docosahexaenoic acid also increased (p < 0.05). Results suggest that COH could lead to premature ovarian aging and provide new insights into the possible prevention of the adverse effects of ovarian hyperstimulation by directing therapeutic applications to the maintenance of the redox balance and fatty acid status, with special attention to paraoxonase proteins and docosahexaenoic acid.

Genes for anti-Müllerian hormone and androgen receptor are underexpressed in human cumulus cells surrounding morphologically highly graded oocytes

Objectives:

The aim of this study was to investigate the expression of genes crucial for the quality of the oocyte and whether expression levels of these genes in cumulus cells can be biological markers for the quality of the oocyte, zygote or embryo, or even for achievement of pregnancy after the assisted reproductive technology procedure. We examined the expression profile of the anti-Müllerian hormone (AMH) gene and its respective receptors: anti-Müllerian hormone receptor type 2 (AMHR2), follicle-stimulating hormone receptor (FSHR) and androgen receptor (AR) in cumulus cells (CCs) surrounding the oocyte, as well as AMH concentrations in follicular fluid of the associated follicle. The obtained gene expression levels were correlated with the morphological quality of the associated oocyte, zygote and embryo as well as with assisted reproductive technology outcome following the intracytoplasmic sperm injection procedure.

Methods:

This study involved 129 cumulus cells and 35 follicular fluid samples, taken from 58 patients undergoing the intracytoplasmic sperm injection procedure. Oocytes, zygotes and embryos were assessed for morphological quality. The relative gene expression of AMH, AMHR2, FSHR and AR was calculated using the delta–delta Ct method. Anti-Müllerian hormone concentrations in follicular fluids were measured by enzyme-linked immunosorbent assay.

Results:

The results yielded suggest a relationship between AMH, AR and oocyte morphology: AMH and AR gene expression levels in CCs surrounding morphologically optimal oocytes were significantly lower than in CCs surrounding oocytes with suboptimal morphology (p = 0.011 and p = 0.008, respectively). Statistically significant positive correlation was found between mRNA expression levels of AMH and FSHR (p < 0.001), AMH and AR (p = 0.001), AMHR2 and FSHR (p < 0.001), AMHR2 and AR (p < 0.001), as well as between FSHR and AR (p < 0.001).

Conclusion:

Assessed results point to AMH and AR relation with oocyte maturity, but not with its fertilization potential, or with embryo quality.

Transcriptome analysis of porcine immature oocytes and surrounding cumulus cells after vitrification and in vitro maturation

Cryopreservation impairs oocyte quality, which may be associated with abnormal gene expression. Currently, alteration of mRNA levels in vitrified porcine oocytes has not been well characterized. The aim of this study was to analyze transcriptome profiles with RNA sequencing (RNA-seq) in porcine immature oocytes and their surrounding cumulus cells (CCs) after vitrification and in vitro maturation (IVM). There were 19 upregulated and 18 downregulated genes differentially expressed in vitrified oocytes, with no significant GO enrichment or KEGG pathway identified for these genes. In addition, CCs derived from vitrified oocytes had 40 significantly upregulated and 100 significantly downregulated genes. In total, 7 GO terms were significantly enriched in molecular function and biological process, and only MAPK signaling pathway reached significant enrichment based on KEGG analysis. Moreover, selected differentially expressed genes had similar expression patterns through comparison between results from qRT-PCR and RNA-Seq. In conclusion, our data provided detailed information on mRNA transcriptomes in porcine immature oocytes and CCs after vitrification and IVM, which offered now insights regarding reduced developmental potential of the vitrified oocytes.

Functional signaling and gene regulatory networks between the oocyte and the surrounding cumulus cells

Background

The maturation and successful acquisition of developmental competence by an oocyte, the female gamete, during folliculogenesis is highly dependent on molecular interactions with somatic cells. Most of the cellular interactions identified, thus far, are modulated by growth factors, ions or metabolites. We hypothesized that this interaction is also modulated at the transcriptional level, which leads to the formation of gene regulatory networks between the oocyte and cumulus cells. We tested this hypothesis by analyzing transcriptome data from single oocytes and the surrounding cumulus cells collected from antral follicles employing an analytical framework to determine interdependencies at the transcript level.

Results

We overlapped our transcriptome data with putative protein-protein interactions and identified hundreds of ligand-receptor pairs that can transduce paracrine signaling between an oocyte and cumulus cells. We determined that 499 ligand-encoding genes expressed in oocytes and cumulus cells are functionally associated with transcription regulation (FDR < 0.05). Ligand-encoding genes with specific expression in oocytes or cumulus cells were enriched for biological functions that are likely associated with the coordinated formation of transzonal projections from cumulus cells that reach the oocyte’s membrane. Thousands of gene pairs exhibit significant linear co-expression (absolute correlation > 0.85, FDR < 1.8 × 10^− 5) patterns between oocytes and cumulus cells. Hundreds of co-expressing genes showed clustering patterns associated with biological functions (FDR < 0.5) necessary for a coordinated function between the oocyte and cumulus cells during folliculogenesis (i.e. regulation of transcription, translation, apoptosis, cell differentiation and transport).

Conclusion

Our analyses revealed a complex and functional gene regulatory circuit between the oocyte and surrounding cumulus cells. The regulatory profile of each cumulus-oocyte complex is likely associated with the oocytes’ developmental potential to derive an embryo.

Electronic supplementary material

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

Cumulus cell transcriptome profiling is not predictive of live birth after in vitro fertilization: a paired analysis of euploid sibling blastocysts

OBJECTIVE

To compare the transcriptome of cumulus cells associated with a euploid embryo that resulted in live birth with that of a sibling euploid embryo without sustained implantation.

DESIGN

Paired analysis.

SETTING

Academic institution.

PATIENT(S)

Couples undergoing in vitro fertilization (IVF)/intracytoplasmic sperm injection with preimplantation genetic screening with female age ≤42 years and normal ovarian reserve.

INTERVENTION(S)

Transcriptome profiling of cumulus cells from sibling oocytes for correlation with live birth after euploid blastocyst transfer. Embryos were individually cultured to facilitate association with clinical outcomes. The cumulus cell transcriptome from the embryo resulting in live birth was compared with that of its sibling embryo without sustained implantation to investigate potential biomarkers that may aid in embryo selection.

MAIN OUTCOME MEASURE(S)

Differential gene expression in cumulus cells associated with a euploid embryo resulting in live birth and its sibling euploid embryo without sustained implantation using next-generation RNA sequencing (RNAseq).

RESULT(S)

Cumulus cell RNAseq of 34 samples (from 17 patients) generated an average of 10.4 ± 4 × 10⁶ reads per sample. A total of 132 differentially expressed genes between sibling embryos that resulted in a live birth and those that did not were identified (P<.05). However, after correcting for multiple testing none of the genes remained significantly differentially expressed (false discovery rate <.05).

CONCLUSION(S)

The RNAseq profiles were similar between cumulus cells associated with a euploid embryo resulting in live birth and its sibling embryo that did not sustain implantation. The cumulus cell transcriptome is not predictive of live birth within an individual patient's cohort of euploid embryos.

Time-lapse observation and transcriptome analysis of a case with repeated multiple pronuclei after IVF/ICSI

PURPOSE

The purpose of this study was to investigate the cause of repeated multipronucleus (MPN) formation in zygotes in a patient after both in vitro fertilization (IVF) and intracytoplasmic sperm injection (ICSI).

METHOD

This is a case study. A patient had unexplained primary infertility with recurring total MPN zygotes after IVF and ICSI cycles. Time-lapse monitoring of pronucleus formation was carried out. Embryos developed from MPN zygotes were analyzed by fluorescence in situ hybridization (FISH). Single-cell RNA-seq analysis was used to identify gene expression profiles of the patient's oocyte and zygote, and these were compared to the data from oocytes and zygotes from donors with normal fertilization (patient, n = 1; donors, n = 4). Oocyte-specific genes with differential expression were selected by the Amazonia!

DATABASE

RESULTS

From time-lapse analysis, we observed the formation of multiple micronuclei near the site of the second polar body extrusion. These micronuclei migrated, expanded, and juxtaposed with the male pronucleus leading to a multipronucleus. None of these MPN zygotes could develop to the blastocyst stage, and FISH analysis revealed a chaotic chromosomal complement in the arrested embryos. RNA-seq analysis showed 113 differentially expressed genes (DEGs) between the patient and the donor oocytes and zygotes. Moreover, 25 of the 113 DEGs were unique or highly expressed in oocytes and early embryos. From 25 DEGs, three genes, DYNC2LI1, NEK2, and CCNH, which are involved in meiosis and the chromosome separation process, were further validated by real-time PCR.

CONCLUSION

We identified several candidate genes affecting pronucleus formation as a new cause of infertility.

A curated transcriptomic dataset collection relevant to embryonic development associated with in vitro fertilization in healthy individuals and patients with polycystic ovary syndrome

The collection of large-scale datasets available in public repositories is rapidly growing and providing opportunities to identify and fill gaps in different fields of biomedical research. However, users of these datasets should be able to selectively browse datasets related to their field of interest. Here we made available a collection of transcriptome datasets related to human follicular cells from normal individuals or patients with polycystic ovary syndrome, in the process of their development, during in vitro fertilization. After RNA-seq dataset exclusion and careful selection based on study description and sample information, 12 datasets, encompassing a total of 85 unique transcriptome profiles, were identified in NCBI Gene Expression Omnibus and uploaded to the Gene Expression Browser (GXB), a web application specifically designed for interactive query and visualization of integrated large-scale data. Once annotated in GXB, multiple sample grouping has been made in order to create rank lists to allow easy data interpretation and comparison. The GXB tool also allows the users to browse a single gene across multiple projects to evaluate its expression profiles in multiple biological systems/conditions in a web-based customized graphical views. The curated dataset is accessible at the following link: http://ivf.gxbsidra.org/dm3/landing.gsp.

Inside the granulosa transcriptome

The somatic component of follicular structure is a mixture of different cell types, represented by Granulosa cells (GCs) that are the paracrine regulators of the oocyte growth. GCs finely support this process by a continuous bidirectional talk with oocyte, which ensure oocyte quality and competence. Specific pathways are involved in the cross-talk and in both GCs and oocyte development. This review summarizes data from GCs gene expression analysis concerning both their physiological role and their interaction with oocyte. We also explore the CGs transcriptome modifications induced by controlled ovarian stimulation (COS) or pathological conditions and their impact in reproduction. The transcriptome analysis of GCs could be a powerful tool to improve our knowledge about the pathways involved in oocyte development. This approach, associated with new technologies as RNA-seq could allow the identifications of new noninvasive biological markers of oocyte quality to increase the efficiency of clinical IVF. Moreover, GCs expression analysis could be useful to shed light on new therapeutic targets by providing new options for the treatment of infertility.

High-throughput analysis of ovarian granulosa cell transcriptome

The quality of follicular oocytes depends on interactions with surrounding granulosa cells. Development of molecular techniques and methods enables better understanding of processes underlying mammalian reproduction on cellular level. The success in reproductive biology and medicine in different species depends on reliable assessment of oocyte and embryo viability which presently mainly bases on embryo morphology. Although successful pregnancies have been achieved using this approach, its precision still should be improved and completed with other, more objective, and accurate assessment strategies. Global profiling of gene expression in follicular cumulus cells using microarrays is continuously leading to the establishment of new biomarkers which can be used to select oocytes with highest developmental potential. Even more potential applications and greater precision could be achieved using next generation sequencing (NGS) of granulosa and cumulus cell RNA (RNA-seq). However, due to the high cost, this method is not used as frequently as microarrays at the moment. In any case, high-throughput technologies offer the possibilities and advantages in ovarian somatic cell analysis on scale that has not been noted so far. The aim of this work is to present current directions and examples of global molecular profiling of granulosa cells and underline its impact on reproductive biology and medicine.