Cumulus cell gene expression associated with pre-ovulatory acquisition of developmental competence in bovine oocytes

Phthalate exposure and risk of ovarian dysfunction in endometriosis: human and animal data

Objective: We aimed to explore the correlations between and possible mechanisms of common environmental endocrine disruptors, phthalates, and ovarian dysfunction in endometriosis. Methods: Subjects were included in the case group (n = 107) who were diagnosed with endometriosis by postoperative pathology in Fujian Maternal and Child Hospital from February 2018 to February 2021, and the women who were excluded from endometriosis by surgery were as the control group (n = 70). The demographic information of the subjects were evaluated by questionnaire, and the clinical characteristics were evaluated by medical records and 3-year follow-up results. Gas chromatography‒mass spectrometry was used to quantify 10 metabolites of phthalates, including dimethyl ortho-phthalate (DMP), mono-n-methyl phthalate (MMP), dioctyl ortho-phthalate (DEP), mono-ethyl phthalate (MEP), di-n-butyl ortho-phthalate (DBP), mono-butyl phthalate (MBP), benzylbutyl phthalate (BBzP), mono-benzyl; phthalate (MBzP), diethylhexyl phthalate (DEHP) and mono-ethylhexyl phthalate (MEHP), in the urine samples of the subjects. Furthermore, a total of 54 SD rats were exposed to DEHP 0, 5, 50, 100, 250, 500, 1,000, 2000, and 3,000 mg/kg/day for 2 weeks. The SD rats' body weight, oestrus cycle changes, and serum anti-mullerian hormone (AMH) levels were evaluated. After sacrifice, the mass index of the rat uterus and bilateral ovaries were calculated. Finally, bioinformatics analysis of rat ovarian tissues was performed to explore the possible mechanism. SPSS 24.0 (IBM, United States) was used for data analysis. p-value <0.05 was considered statistically significant. Results: The human urinary levels of DMP (p < 0.001), MMP (p = 0.001), DEP (p = 0.003), MEP (p = 0.002), DBP (p = 0.041), MBP (p < 0.001), BBzP (p = 0.009), DEHP (p < 0.001), and MEHP (p < 0.001) were significantly higher in women with endometriosis than in controls. Notably, DEHP was a significant risk factor for endometriosis (OR: 11.0, 95% CI: 5.4-22.6). The area under the ROC curve increased when multiple phthalates were diagnosed jointly, reaching 0.974 as the highest value, which was helpful for the diagnosis of endometriosis. In vivo experiments showed that after DEHP exposure in rats, the mass index of the ovary and uterus decreased in a dose-dependent manner; the oestrus cycle of SD rats was irregularly prolonged and disordered; and the serum AMH level was negatively correlated with the DEHP exposure dose (Rho = -0.8, p < 0.001). Bioinformatics analysis of rat ovarian tissues showed that seven genes involved in the steroid biosynthesis pathway were upregulated and may play a negative role in ovarian function. Conclusion: Exposure to phthalates, especially DEHP, is associated with the occurrence of endometriosis and affects women's reproductive prognosis and ovarian function. The steroid biosynthesis pathway may be related to ovarian dysfunction. The detection of phthalate in urine may become a new biological target for the diagnosis of endometriosis.

TBX2 affects proliferation, apoptosis and cholesterol generation by regulating mitochondrial function and autophagy in bovine cumulus cell

BACKGROUND

T-box transcription factor 2 (TBX2) is a member of T-box gene family whose members are highly conserved in evolution and encoding genes and are involved in the regulation of developmental processes. The encoding genes play an important role in growth and development. Although TBX2 has been widely studied in cancer cell growth and development, its biological functions in bovine cumulus cells remain unclear.

OBJECTIVES

This study aimed to investigate the regulatory effects of TBX2 in bovine cumulus cells.

METHODS

TBX2 gene was knockdown with siRNA to clarify the function in cellular physiological processes. Cell proliferation and cycle changes were determined by xCELLigence cell function analyzer and flow cytometry. Mitochondrial membrane potential and autophagy were detected by fluorescent dye staining and immunofluorescence techniques. Western blot and quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) were used to detect the expression changes of proliferation and autophagy-related proteins. Aadenosine triphosphate (ATP) production, glucose metabolism, and cholesterol synthesis of cumulus cells were measured by optical density and chemiluminescence analysis.

RESULTS

After inhibition of TBX2, the cell cycle was disrupted. The levels of apoptosis, ratio of light chain 3 beta II/I, and reactive oxygen species were increased. The proliferation, expansion ability, ATP production, and the amount of cholesterol secreted by cumulus cells were significantly decreased.

CONCLUSIONS

TBX2 plays important roles in regulating the cells' proliferation, expansion, apoptosis, and autophagy; maintaining the mitochondrial function and cholesterol generation of bovine cumulus cells.

RNA-sequencing reveals genes linked with oocyte developmental potential in bovine cumulus cells

Cumulus cells provide an interesting biological material to perform analyses to understand the molecular clues determining oocyte competence. The objective of this study was to analyze the transcriptional differences between cumulus cells from oocytes exhibiting different developmental potentials following individual in vitro embryo production by RNA-seq. Cumulus cells were allocated into three groups according to the developmental potential of the oocyte following fertilization: (1) oocytes developing to blastocysts (Bl+), (2) oocytes cleaving but arresting development before the blastocyst stage (Bl-), and (3) oocytes not cleaving (Cl-). RNAseq was performed on 4 (Cl-) or 5 samples (Bl+ and Bl-) of cumulus cells pooled from 10 cumulus-oocyte complexes per group. A total of 49, 50, and 18 differentially expressed genes (DEGs) were detected in the comparisons Bl+ versus Bl-, Bl+ versus Cl- and Bl- versus Cl-, respectively, showing a fold change greater than 1.5 at an adjusted p value <0.05. Focussing on DEGs in cumulus cells from Bl+ group, 10 DEGs were common to both comparisons (10/49 from Bl+ vs. Bl-, 10/50 from Bl+ vs. Cl-). These DEGs correspond to 6 upregulated genes (HBE1, ITGA1, PAPPA, AKAP12, ITGA5, and SLC1A4), and 4 downregulated genes (GSTA1, PSMB8, FMOD, and SFRP4) in Bl+ compared to the other groups, from which 7 were validated by quantitative PCR (HBE1, ITGA1, PAPPA, AKAP12, ITGA5, PSMB8 and SFRP4). These genes are involved in critical biological functions such as integrin-mediated cell adhesion, oxygen availability, IGF and Wnt signaling or PKA pathway, highlighting specific biological processes altered in incompetent in vitro maturation oocytes.

Follicular-fluid proteomics during equine follicle development

The complex composition of the follicular fluid (FF), the intimate proximity to the oocyte, and the continual changes in their composition have a major effect on folliculogenesis and oogenesis. To date, the profiling of FF proteomes during follicle selection, development, and ovulation has not been comprehensively investigated. Therefore, a shotgun proteomics approach and bioinformatics analyses were used to profile the proteomes of equine FF harvested in vivo from follicles at the following development stages: predeviation (18-20 mm), deviation (22-25 mm), postdeviation (26-29 mm), preovulatory (30-35 mm), and impending ovulation. A total of 294 proteins were detected in FF (FDR <1%), corresponding to 65 common proteins and 124, 142, 167, 132, and 142 proteins in the predeviation, deviation, postdeviation, preovulatory, and impending ovulation groups, respectively. The higher expression of properdin and several other proteins belonging to the complement system during the deviation time and ovulation suggested their contribution in the selection of the future dominant follicle and ovulation. Apolipoprotein A-1 and antithrombin-III appeared to be important throughout folliculogenesis. The "complement and coagulation cascades" was the major KEGG pathway across all stages of follicle development. The significant expression of several proteins belonging to the serine-type endopeptidase indicated their likely contribution to follicle and oocyte development. Our data provide an extensive description and functional analyses of the equine FF proteome during follicle selection, development, and ovulation. This information will help improve understanding of the ovarian function and ovulatory dysfunctions and might serve as a reference for future biomarker discovery for oocyte quality assessment.

Nurturing the egg: the essential connection between cumulus cells and the oocyte

The determinants of oocyte quality remain uncertain. Under suitable conditions, which have yet to be defined, the gamete grows and acquires the competence to resume meiosis, be fertilised and undergo embryonic development at least beyond genome activation, after which the blastomere is autonomous enough to adapt to the specificity of its environment. This review describes the central role played by the oocyte in reproductive success and how communication between cumulus cells and the oocyte are essential to proper oogenesis and the quality of the resulting gamete. While most attempts to improve oocyte quality have been directed at gonadotrophin-based systemic endocrine signalling, it is proposed that parallel control of fertility may act locally within ovarian follicles through intimate cooperation between somatic cells and the oocyte via the network of transzonal projections. This intercellular communication may prove to be more sensitive to environmental conditions than systemic endocrine signalling, which is essential for many non-reproductive tissues.

Transcriptome profile of goat folliculogenesis reveals the interaction of oocyte and granulosa cell in correlation with different fertility population

To understand the molecular and genetic mechanisms related to the litter size in one species of two different populations (high litter size and low litter size), we performed RNA-seq for the oocytes and granulosa cells (GCs) at different developmental stages of follicle, and identified the interaction of genes from both sides of follicle (oocyte and GCs) and the ligand-receptor pairs from these two sides. Our data were very comprehensive to uncover the difference between these two populations regarding the folliculogenesis. First, we identified a set of potential genes in oocyte and GCs as the marker genes which can be used to determine the goat fertility capability and ovarian reserve ability. The data showed that GRHPR, GPR84, CYB5A and ERAL1 were highly expressed in oocyte while JUNB, SCN2A, MEGE8, ZEB2, EGR1and PRRC2A were highly expressed in GCs. We found more functional genes were expressed in oocytes and GCs in high fertility group (HL) than that in low fertility group (LL). We uncovered that ligand-receptor pairs in Notch signaling pathway and transforming growth factor-β (TGF-β) superfamily pathways played important roles in goat folliculogenesis for the different fertility population. Moreover, we discovered that the correlations of the gene expression in oocytes and GCs at different stages in the two populations HL and LL were different, too. All the data reflected the gene expression landscape in oocytes and GCs which was correlated well with the fertility capability.

Oocyte Selection for In Vitro Embryo Production in Bovine Species: Noninvasive Approaches for New Challenges of Oocyte Competence

Simple Summary

The efficiency of producing embryos using in vitro technologies in cattle species remains lower when compared to mice, indicating that the proportion of female gametes that fail to develop after in vitro manipulation is considerably large. Considering that the intrinsic quality of the oocyte is one of the main factors affecting embryo production, the precise identification of noninvasive markers that predict oocyte competence is of major interest. The aim of this review was to explore the current literature on different noninvasive markers associated with oocyte quality in the bovine model. Apart from some controversial findings, the presence of cycle-related structures in ovaries, a follicle size between 6 and 10 mm, a large slightly expanded investment without dark areas, large oocyte diameter (>120 microns), dark cytoplasm, and the presence of a round and smooth first polar body have been associated with better embryonic development. In addition, the combination of oocyte and zygote selection, spindle imaging, and the anti-Stokes Raman scattering microscopy together with studies decoding molecular cues in oocyte maturation have the potential to further optimize the identification of oocytes with better developmental competence for in vitro technologies in livestock species.

Abstract

The efficiency of producing embryos using in vitro technologies in livestock species rarely exceeds the 30–40% threshold, indicating that the proportion of oocytes that fail to develop after in vitro fertilization and culture is considerably large. Considering that the intrinsic quality of the oocyte is one of the main factors affecting blastocyst yield, the precise identification of noninvasive cellular or molecular markers that predict oocyte competence is of major interest to research and practical applications. The aim of this review was to explore the current literature on different noninvasive markers associated with oocyte quality in the bovine model. Apart from some controversial findings, the presence of cycle-related structures in ovaries, a follicle size between 6 and 10 mm, large number of surrounding cumulus cells, slightly expanded investment without dark areas, large oocyte diameter (>120 microns), dark cytoplasm, and the presence of a round and smooth first polar body have been associated with better competence. In addition, the combination of oocyte and zygote selection via brilliant cresyl blue (BCB) test, spindle imaging, and the anti-Stokes Raman scattering microscopy together with studies decoding molecular cues in oocyte maturation have the potential to further optimize the identification of oocytes with better developmental competence for in-vitro-derived technologies in livestock species.

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

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

Combinational approach of retrospective clinical evidence and transcriptomics highlight AMH superiority to FSH, as successful ICSI outcome predictor

Objective

Combination of transcriptomic and retrospective clinical data, to assess anti-Mullerian hormone (AMH) functionality at a cumulus cell level and evaluate AMH potential as a suitable marker for IVF outcomes (oocytes retrieved, number of day 3 embryos, gestation outcomes).

Design

Raw RNA-sequencing data of cumulus cells sourced from younger (n = 10) patient group (group A) (age 29 (1 year of age), baseline FSH 7.4 (0.5 mIU/ml), AMH 4.67 (1.56 ng/ml)) and older (n = 10) patient group (group B) (age 43 (± 0.55 years of age), baseline FSH 8 (0.8 mIU/ml), AMH 1.07 (0.44 ng/ml)) were employed to derive transcriptomic differences among high vs. low AMH groups. We collected retrospectively patient data from 80 infertile patients selected according to pre-specified inclusion criteria.

Setting

Publicly available raw RNA-sequencing data were retrieved from the SRA database of NCBI resource GEO Accession (GSM21575/35-44; GEO Accession: GSM21575/45-55). Retrospective data were collected from referrals to the Institute of Reproductive Medicine, Lito Hospital of Athens and the Institute of Life, Iaso Hospital of Athens, between the periods of March 2015 and April 2018.

Intervention(s)

A fixed human menopausal gonadotropin (hMG) antagonist protocol was used for all patients. All patients had serum AMH levels measured within a 3-month period prior to stimulation and serum levels of FSH and estradiol (day 2 of menstrual cycle; E2) (Clinical Trial code NV24042014).

Main outcome measure(s)

The primary outcomes were identification of transcriptomic variations among high (group A) vs. low (group B) AMH patients. Retrospective data primary outcomes were number of oocytes retrieved, fertilized successfully (grades A and B, day 2 embryos), and total number of day 3 embryos. Secondary outcome was live birth rate. Finally, we compared primary outcomes with AMH and FSH level as well as their genetic pathways (interacting genes) to demonstrate the predictive accuracy.

Results

Essential players of the AMH signaling cascade, namely, SMAD1, SMAD4, SMAD5, ALK1, and LEF1, were significantly upregulated in group A (n 10) transcriptome. This biological clue was further supported by retrospective clinical data (n 80 participants), where AMH was positively correlated with both oocytes retrieved and fertilized as well as number of day 3 (grades A and B) embryos from patients undergoing IVF, in a statistically significant manner. AMH was further positive trend of association with successful pregnancy outcomes.

Conclusion

Overall, this study offers new insight on AMH effects upon cumulus cells and new aspects on how AMH might promote oocyte integrity and embryo viability at a biochemical level as well as add to the current body of evidence supporting AMH clinical potential as a more sensitive marker of IVF outcomes in comparison with FSH, regarding numbers of oocytes received and high-quality day 2 and day 3 embryos.

Electronic supplementary material

The online version of this article (10.1007/s10815-020-01802-w) contains supplementary material, which is available to authorized users.

Chromosome Missegregation in Single Human Oocytes Is Related to the Age and Gene Expression Profile

The growing trend for women to postpone childbearing has resulted in a dramatic increase in the incidence of aneuploid pregnancies. Despite the importance to human reproductive health, the events precipitating female age-related meiotic errors are poorly understood. To gain new insight into the molecular basis of age-related chromosome missegregation in human oocytes, we combined the transcriptome profiles of twenty single oocytes (derived from females divided into two groups according to age <35 and ≥35 years) with their chromosome status obtained by array comparative genomic hybridization (aCGH). Furthermore, we compared the transcription profile of the single oocyte with the surrounding cumulus cells (CCs). RNA-seq data showed differences in gene expression between young and old oocytes. Dysregulated genes play a role in important biological processes such as gene transcription regulation, cytoskeleton organization, pathways related to RNA maturation and translation. The comparison of the transcription profile of the oocyte and the corresponding CCs highlighted the differential expression of genes belonging to the G protein-coupled receptor superfamily. Finally, we detected the loss of a X chromosome in two oocytes derived from women belonging to the ≥35 years age group. These aneuploidies may be caused by the detriment of REEP4, an endoplasmic reticulum protein, in women aged ≥35 years. Here we gained new insight into the complex regulatory circuit between the oocyte and the surrounding CCs and uncovered a new putative molecular basis of age-related chromosome missegregation in human oocytes.

Dual-specificity phosphatase 1 regulates cell cycle progression and apoptosis in cumulus cells by affecting mitochondrial function, oxidative stress, and autophagy

Dual-specificity phosphatase 1 ( DUSP1) is differentially expressed in cumulus cells of different physiological states, but its specific function and mechanism of action remain unclear. In this study, we explored the effects of DUSP1 expression inhibition on cell cycle progression, proliferation, apoptosis, and lactate and cholesterol levels in cumulus cells and examined reactive oxygen species levels, mitochondrial function, autophagy, and the expression of key cytokine genes. The results showed that inhibition of DUSP1 in cumulus cells caused abnormal cell cycle progression, increased cell proliferation, decreased apoptosis rates, increased cholesterol synthesis and lactic acid content, and increased cell expansion. The main reason for these effects was that inhibition of DUSP1 reduced ROS accumulation, increased glutathione level and mitochondrial membrane potential, and reduced autophagy levels in cells. These results indicate that DUSP1 limits the biological function of bovine cumulus cells under normal physiological conditions and will greatly contribute to further explorations of the physiological functions of cumulus cells and the interactions of the cumulus-oocyte complex.

The use of adenosine to inhibit oocyte meiotic resumption in Bos taurus during pre-IVM and its potential to improve oocyte competence

One of the major challenges of artificial reproductive technologies is to develop new methods for producing greater numbers of embryos. An oocyte fosters the ability to develop into an embryo before oocyte meiotic resumption. The aim of the present study was to assess the effect of adenosine (ADO), a purine nucleoside found in follicular fluid, on the inhibition of oocyte meiotic resumption and the production of blastocysts. The results showed the efficacy of ADO to inhibit oocyte meiotic resumption. The use of ADO (3 mM) during a pre-in vitro maturation (pre-IVM) culture period of 6 h resulted in a significant increase (p < 0.05) of blastocysts compared to control conditions with no pre-IVM culture period. No effect on the percentage of cleavage was observed. The effect of adenosine on blastocyst yield was time- and concentration-dependent with an optimum effect at 3 mM for 6 h. Supplementing the ADO pre-IVM culture medium with estradiol, follicle-stimulating hormone, progesterone, epidermal growth factor, insulin-like growth factor-2 or reelin did not improve the blastocyst yield. Transcriptional analyses of ADO-treated cumulus cells revealed that NRP1, RELN, MAN1A1, THRA and GATM were up-regulated. Finally, bioinformatic analysis identified mitochondrial function as the top canonical pathway affected by ADO. This opens up new opportunities for further investigations.

Trio, a novel high fecundity allele: I. Transcriptome analysis of granulosa cells from carriers and noncarriers of a major gene for bovine ovulation rate

A major gene for bovine ovulation rate has been mapped to a 1.2 Mb region of chromosome 10. Screening of coding regions of positional candidate genes within this region failed to reveal a causative polymorphism, leading to the hypothesis that the phenotype results from differences in candidate gene expression rather than alteration of gene structure. This study tested differences in expression of positional candidate genes in granulosa cells between carriers and noncarriers of the high fecundity allele, as well as characterizing differences in the transcriptomic profile between genotypes. Five carriers and five noncarriers, female descendants of "Trio," a carrier of the high fecundity allele were initially used in an RNA-seq analysis of gene expression. Four of ten samples were contaminated with theca cells, so that six samples were used in the final analysis (three of each genotype). Of 14 973 genes expressed, 143 were differentially expressed (false discovery rate P < 0.05) in carriers versus noncarriers. Among the positional candidate genes, SMAD6 was 6.6-fold overexpressed in the carriers compared to noncarriers (P < 5 × 10-5). This result was replicated in an independent group of 12 females (7 carriers and 5 noncarriers) using quantitative real-time PCR; SMAD6 was 9.3-fold overexpressed in carriers versus noncarriers (P = 1.17 × 10-6). Association of overexpression of SMAD6, an inhibitor of the BMP/SMAD signaling pathway, with high ovulation rate corresponds well with disabling mutations in ligands (BMP15 and GDF9) and a receptor (BMPR1B) of this pathway that cause increased ovulation rate in sheep.

Transcriptomics-genomics data integration and expression quantitative trait loci analyses in oocyte donors and embryo recipients for improving invitro production of dairy cattle embryos

In this paper we first provide a brief review of main results from our previously published studies on genome-wide gene expression (transcriptomics) in donor and recipient cattle used in invitro production (IVP) of embryos and embryo transfer (ET). Then, we present novel results from applying integrative systems genomics and biological analyses where transcriptomics data are combined with genomic data in both donor and recipient cattle to map expression quantitative trait loci (eQTLs). The eQTLs are genetic markers that can regulate or control the expression of genes in the entire genome, via complex molecular mechanisms, and thus can act as a powerful tool for genomic and gene-assisted selection. We identified significant eQTLs potentially controlling the expression of 13 candidate genes for donor cow quality (IVP parameters; e.g. cyclin B1 (CCNB1), outer dense fiber of sperm tails 2 like (ODF2L)) and 19 candidate genes for recipient cows quality (endometrial receptivity; e.g. ER membrane protein complex subunit 9 (EMC9), mannosidase beta (MANBA), peptidase inhibitor 16 (PI16)). Annotation and colocation of detected eQTLs show that some of the eQTLs are in the same genomic regions previously reported as QTLs for reproduction-related traits. However, eQTLs and the candidate genes identified should be further validated in larger populations before implementation as genetic markers or used in genomic selection for improving IVP and ET performance.

RNA Sequencing-Based Whole-Transcriptome Analysis of Friesian Cattle Fed with Grape Pomace-Supplemented Diet

Simple Summary

Grape pomace (GPO) is an important source of polyphenols which are known to have antioxidant properties. In the past decade, GPO has received some attention as a bioactive dietary component in farm animals’ diet. In this study, we have analyzed the whole-transcriptome of Friesian calves fed with a GPO-supplemented diet using RNA-sequencing. We noted that the most affected pathway was the cholesterol lipid biosynthesis and this effect was consistent with a reduction in both serum cholesterol and lipid oxidation in the carcasses. This study provides evidence on the antioxidant property of GPO-supplemented diet, from a molecular biology standpoint.

Abstract

Grape pomace (GPO), the main by-product of the wine making process, is a rich source of polyphenols with potent antioxidant properties. Recently, GPO has emerged as a potential feed additive in livestock nutrition, with several reports describing its beneficial effects on animals’ overall health status or production traits. However, little is known about it from a molecular biology standpoint. In the present study, we report the first RNA sequencing-based whole-transcriptome profiling of Friesian calves fed with a GPO-supplemented diet. We identified 367 differentially expressed genes (p < 0.05) in the GPO-supplemented calves (n = 5), when compared with unsupplemented control group (n = 5). The pathway analysis showed that ‘cholesterol lipid biosynthesis’ was the most negatively-enriched (p < 0.001) pathway in the GPO-supplemented animals. In specific terms, five important genes coding for cholesterol biosynthesis enzymes, namely the Farnesyl-diphosphate Farnesyltransferase 1 (FDFT-1), Squalene Epoxidase (SQLE), NAD(P)-dependent Steroid Dehydrogenase-like (NSDHL), Methylsterol Monooxygenase (MSMO)-1, and Sterol-C5-desaturase (SC5D), two major transcription factors (the Sterol Regulatory Element-binding Transcription Factor 1 and 2), as well as the Low-Density Lipoprotein Receptor (LDLR), were all downregulated following GPO supplementation. Such an effect was mirrored by a reduction of blood cholesterol levels (p = 0.07) and a lowered (p < 0.001) Malondialdehyde (lipid oxidation marker) level in carcasses. We provide evidence on the effects of GPO-supplemented diets on the whole-transcriptome signature in veal calves, which mainly reflects an antioxidant activity.

Lipid Identification and Transcriptional Analysis of Controlling Enzymes in Bovine Ovarian Follicle

Ovarian follicle provides a favorable environment for enclosed oocytes, which acquire their competence in supporting embryo development in tight communications with somatic follicular cells and follicular fluid (FF). Although steroidogenesis in theca (TH) and granulosa cells (GC) is largely studied, and the molecular mechanisms of fatty acid (FA) metabolism in cumulus cells (CC) and oocytes are emerging, little data is available regarding lipid metabolism regulation within ovarian follicles. In this study, we investigated lipid composition and the transcriptional regulation of FA metabolism in 3–8 mm ovarian follicles in bovine. Using liquid chromatography and mass spectrometry (MS), 438 and 439 lipids were identified in FF and follicular cells, respectively. From the MALDI-TOF MS lipid fingerprints of FF, TH, GC, CC, and oocytes, and the MS imaging of ovarian sections, we identified 197 peaks and determined more abundant lipids in each compartment. Transcriptomics revealed lipid metabolism-related genes, which were expressed constitutively or more specifically in TH, GC, CC, or oocytes. Coupled with differential lipid composition, these data suggest that the ovarian follicle contains the metabolic machinery that is potentially capable of metabolizing FA from nutrient uptake, degrading and producing lipoproteins, performing de novo lipogenesis, and accumulating lipid reserves, thus assuring oocyte energy supply, membrane synthesis, and lipid-mediated signaling to maintain follicular homeostasis.

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.

Effects of aging on gene expression and mitochondrial DNA in the equine oocyte and follicle cells

We hypothesised that advanced mare age is associated with follicle and oocyte gene alterations. The aims of the study were to examine quantitative and temporal differences in mRNA for LH receptor (LHR), amphiregulin (AREG) and epiregulin (EREG) in granulosa cells, phosphodiesterase (PDE) 4D in cumulus cells and PDE3A, G-protein-coupled receptor 3 (GPR3), growth differentiation factor 9 (GDF9), bone morphogenetic protein 15 (BMP15) and mitochondrial (mt) DNA in oocytes. Samples were collected from dominant follicles of Young (3-12 years) and Old (≥20 years) mares at 0, 6, 9 and 12h after administration of equine recombinant LH. LHR mRNA declined after 0h in Young mares, with no time effect in Old mares. For both ages, gene expression of AREG was elevated at 6 and 9h and EREG was expression was elevated at 9h, with higher expression in Old than Young mares. Cumulus cell PDE4D expression increased by 6h (Old) and 12h (Young). Oocyte GPR3 expression peaked at 9 and 12h in Young and Old mares, respectively. Expression of PDE3A increased at 6h, with the increase greater in oocytes from Old than Young mares at 6 and 9h. Mean GDF9 and BMP15 transcripts were higher in Young than Old, with a peak at 6h. Copy numbers of mtDNA did not vary over time in oocytes from Young mares, but a temporal decrease was observed in oocytes from Old mares. The results support an age-associated asynchrony in the expression of genes that are essential for follicular and oocyte maturation before ovulation.

Comparative analysis of granulosa cell gene expression in association with oocyte competence in FSH-stimulated Holstein cows

Ovarian stimulation with exogenous FSH followed by FSH withdrawal or ‘coasting’ is an effective means of increasing the number of oocytes obtainable for the in vitro production of cattle embryos. However, the quality of the oocytes thus obtained varies considerably from one cow to the next. The aim of the present study was to gain a better understanding of the follicular conditions associated with low oocyte developmental competence. Granulosa cells from 94 Holstein cows in a commercial embryo production facility were collected following ovarian stimulation and coasting. Microarray analysis showed 120 genes expressed with a differential of at least 1.5 when comparing donors of mostly competent with donors of mostly incompetent oocytes. Using ingenuity pathway analysis, we revealed the main biological functions and potential upstream regulators that distinguish donors of mostly incompetent oocytes. These are involved in cell proliferation, apoptosis, lipid metabolism, retinol availability and insulin signalling. In summary, we demonstrated that differences in follicle maturity at collection could explain differences in oocyte competence associated with individual animals. We also revealed deficiencies in lipid metabolism and retinol signalling in granulosa cells from donors of mostly incompetent oocytes.

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.

Differences in cumulus cell gene expression indicate the benefit of a pre-maturation step to improve in-vitro bovine embryo production

STUDY QUESTION

Does the gene expression profile of cumulus cells (CC) accompanying oocytes with different degrees of chromatin compaction within the germinal vesicle (GV) reflect the oocyte's quality and response in culture during in-vitro embryo production (IVP).

SUMMARY ANSWER

The transcriptomic profile of the CC is related to oocyte competence, setting the stage for the development of customized pre-maturation strategies to improve IVP.

WHAT IS KNOWN ALREADY

Oocytes complete the acquisition of their competence during antral follicle development. During this period, the chromatin configuration within the GV changes dynamically and is indicative of oocyte's developmental potential. The interactions between somatic and germ cells modulate chromatin morphology and function and are critical for acquisition of oocyte competence.

STUDY DESIGN, SIZE, DURATION

Bovine cumulus-oocyte complexes (COC) were isolated from 0.5 to 6 mm antral follicles. Surrounding CC were separated from the oocyte and classified as GV0, GV1, GV2 and GV3 according to the degree of the oocyte's chromatin compaction.

PARTICIPANTS/MATERIALS, SETTING, METHOD

RNA extracted from CC of each group was amplified and hybridized on a bovine embryo-specific 44 K Agilent slide. The CC_GV1, CC_GV2 and CC_GV3 classes were each hybridized against the CC_GV0 class, representing an early oocyte differentiation stage with poor development competence. The data were normalized and fold changes of the differentially expressed genes were determined. Microarray data were validated using quantitative RT-PCR on selected targets. Microarray data were further analyzed through: (i) between-group analysis (BGA), which classifies the samples according to their transcriptomic profiles; (ii) cluster analysis according to the expression profile of each gene; and (iii) Ingenuity Pathway Analysis (IPA) to study gene regulation patterns and predicted functions. Furthermore, CC of each GV group were cultured and apoptotic cells were assessed after 3 h by caspase analysis. Finally, based on the analysis of CC transcriptomic profiles and the relationship between morphological features of the COC and the oocyte chromatin configuration, a customized, stage-dependent oocyte pre-maturation (pre-IVM) system was used to improve oocyte developmental potential before IVM. For this, the blastocyst rate and quality were assessed after in-vitro maturation and fertilization of pre-matured oocytes.

MAIN RESULTS AND THE ROLE OF CHANCE

Overall, quantitative RT-PCR results of a subset of five selected genes were consistent with the microarray data. Clustering analysis generated 16 clusters representing the main profiles of transcription modulation. Of the 5571 significantly differentially expressed probes, the majority (25.49%) best fitted with cluster #6 (downregulation between CC_GV0 and CC_GV1 and stable low levels in successive groups). IPA identified the most relevant functions associated with each cluster. Genes included in cluster #1 were mostly related to biological processes such as 'cell cycle' and 'cell death and survival', whereas genes included in cluster #5 were mostly related to 'gene expression'. Interestingly, 'lipid metabolism' was the most significant function identified in clusters #6, #9 and #12. IPA of gene lists obtained from each contrast (i.e., CC_GV0 vs. CC_GV1; CC_GV0 vs. CC_GV2; CC_GV0 vs. CC_GV3) revealed that the main affected function in each contrast was 'cell death and survival'. Importantly, apoptosis was predicted to be inhibited in CC_GV1 and CC_GV2, but activated in CC_GV3. Caspase analysis indicated that a low percentage of CC_GV0 was prone to undergo apoptosis but apoptosis increased significantly in CC from oocytes with condensed chromatin, reaching a peak in CC_GV3 (P < 0.05). Finally, the tailored oocyte pre-maturation strategy, based on morphological features of the COC and the oocyte chromatin configuration, demonstrated that pre-IVM improved the developmental capability of oocytes at early stages of differentiation (GV1-enriched COC) but was detrimental for oocytes at more advanced stages of development (GV2 and GV3-enriched COC).

LARGE SCALE DATA

The data are available through the GEO series accession number GSE79886.

LIMITATIONS, REASONS FOR CAUTION

This study was conducted with bovine samples. Whether or not the results are applicable to human oocytes requests further elucidation. Embryo transfer experiments are required to determine whether the improvement in blastocyst rates in the tailored system leads to increased live birth rates.

WIDER IMPLICATIONS OF THE FINDINGS

The identification of multiple non-invasive biomarkers predictive of oocyte quality can greatly strengthen the pre-IVM approach aimed to improve IVM outcomes. These results have potentially important implications in treating human infertility and in developing breeding schemes for domestic mammals.

STUDY FUNDING/COMPETING INTERESTS

This work was supported in part by NSERC Strategic Network EmbryoGENE, Canada and in part by CIG-Marie Curie Actions-Reintegration Grants within the EU 7FP (n. 303640, 'Pro-Ovum'). The authors declare no potential conflict of interest.

FSH stimulation of anestrous cats improves oocyte quality and development of parthenogenetic embryos

In the domestic cat, the efficiency of in vitro embryo production systems is negatively affected during the nonbreeding season. The objective of this research was to evaluate the effect of FSH stimulation in anestrous cats, on quality of cumulus-oocyte complexes (COCs) and in vitro developmental competence after parthenogenetic activation. To accomplish this purpose, anestrous cats were grouped into: (1) FSH treated (serial doses of 5 mg of porcine FSH each, every 24 hours, for 4 days) and (2) untreated control. The COCs were classified morphologically and a proportion of grade I and II COCs was used for expression analysis of FSHR, LHCGR, EGFR, PTGS2, EGR1, GDF9, and GATM by RT-qPCR. In addition, another proportion of grade I and II COCs was matured in vitro and used for parthenogenetic activation. After 8 days in culture, blastocyst and hatching blastocyst rates were assessed, and the expression of OCT4, SOX2, NANOG, CDX2, and GATA6 was evaluated. The COCs in the FSH group had an enhanced quality, a higher expression of LHCGR and a lower expression of GATM than did COCs from the control group (P < 0.05). Furthermore, embryos in the FSH group had increased blastocyst and hatching blastocyst rates, and those embryos had a higher expression of OCT4 and GATA than their counterparts from the control group (P < 0.05). In conclusion, ovarian stimulation of anestrous cats with FSH improved quality and increased the expression of LHCGR in COCs. The enhanced in vitro developmental competence, after parthenogenetic activation of oocytes from FSH-treated cats, coincided with an increased expression of OCT4 and GATA6 in blastocysts and hatching blastocysts.

Transcriptome meta-analysis of three follicular compartments and its correlation with ovarian follicle maturity and oocyte developmental competence in cows

Oocyte developmental competence in superstimulated cows is dependent in part on the duration of the FSH coasting. FSH coasting refers to superstimulation with FSH (2 days of endogenous FSH following follicle ablation and 3 days of FSH injections) followed by no FSH for a specific duration. The optimal duration varies among individuals. FSH coasting appears to modulate the transcriptome of different follicular compartments, which cooperate as a single functional unit. However, the integrative effects of FSH coasting on different follicular compartments remain ambiguous. Meta-analysis of three independent transcriptome studies, each focused on a single cell type (granulosa, cumulus, and oocyte) during FSH coasting, allowed the identification of 12 gene clusters with similar time-course expression patterns in all three compartments. Network analysis identified HNF4A (involved in metabolic functions) and ELAVL1 (an RNA-binding protein) as hub genes regulated respectively upward and downward in the clusters enriched at the optimal coasting time, and APP (involved in mitochondrial functions) and COPS5 (a member of the COP9 signalosome) as hub genes regulated respectively upwards and downwards in the clusters enriched progressively throughout the coasting period. We confirmed the effects on HNF4A downstream targets (TTR, PPL) and other hub genes (ELAVL1, APP, MYC, and PGR) in 30 cows with RT-quantitative PCR. The correlation of hub gene expression levels with FSH coasting indicated that a combination of these genes could predict oocyte competence with 83% sensitivity, suggesting that they are potential biomarkers of follicle differentiation. These findings could be used to optimize FSH coasting on an individual basis.

Molecular markers for oocyte competence in bovine cumulus cells

This study aimed to quantify the expression of candidate genes in cumulus cells (CCs) from cumulus-oocyte complexes (COCs) with high and low potential for in vitro development up to the blastocyst stage. First, the effects of individual culture and biopsy on embryo development were evaluated. Individuals cultured using the well of the well system were compared with individuals cultured in 20 μL droplets (microdroplets) and those cultured in groups (control). Blastocyst rates were lower for the individual culture systems (P < 0.05; well of the well = 17.9%, n = 95; microdrop = 26.3%, n = 95) than for the control group (45.0%, n = 209). Second, the effects of biopsy on embryo production were compared between the control and microdroplet cultures, and no effects (P > 0.05) were observed for either group. Finally, the expression profiles of glypican 4 (GPC4), IGF4-binding protein, follicle-stimulating hormonereceptor, growth hormone receptor, epidermal growth factor receptor, fibroblast growth factor 11, solute carrier family 2 member 1, solute carrier family 2 member 3,sprouty homolog 1, versican, and keratin protein 8 in CCs obtained by biopsy were quantified by real-time polymerase chain reaction. Cumulus cells were categorized on the basis of the fates of the COCs: expanded blastocyst, cleaved and arrested, and uncleaved. The GPC4 gene was overexpressed (P = 0.007) in CCs from oocytes that formed embryos compared with those that produced cleaved and arrested embryos. We concluded that individual culture reduced blastocyst production; however, biopsy did not affect embryo development. The profile of GPC4 expression can be used as a marker to distinguish COCs with potential for embryo development from those with limited developmental potential.

The effects of human recombinant granulocyte-colony stimulating factor treatment during in vitro maturation of porcine oocyte on subsequent embryonic development

Granulocyte colony-stimulating factor (G-CSF) is required for proliferation, differentiation, and survival of cells. It is also a biomarker of human oocyte developmental competence for embryo implantation. In humans, the G-CSF concentration peaks during the ovulatory phase of the ovarian cycle. In this study, the expressions of G-CSF and its receptor were analyzed by polymerase chain reaction in granulosa cells (GCs), CL, cumulus cells (CCs), and oocytes. Cumulus-oocyte complexes were aspirated from antral follicles of 1 to 3 mm (small follicles) and 4 to 6 mm (medium follicles). Cumulus-oocyte complexes from two kinds of follicles were matured in protein-free maturation medium supplemented with various concentrations of G-CSF (0, 10, and 100 ng/mL). By real-time polymerase chain reaction, the expressions of G-CSF and its receptor were detected in GCs, CL, CCs, and oocytes. Interestingly, the G-CSF transcript levels were significantly lower in oocytes than in the other cell types, whereas the G-CSF receptor transcript levels in oocytes were similar to those in GCs. After 44 hours of IVM, no differences in the rate of nuclear maturation were detected; however, the intracellular reactive oxygen species levels in oocytes from both groups of follicles matured with 10 ng/mL of human recombinant G-CSF (hrG-CSF) groups were significantly lower (P < 0.05). After parthenogenetic activation, the cleavage rates were significantly (P < 0.05) higher in 100 ng/mL hrG-CSF-treated small (63.3%) follicles than in 0, 10 ng/mL hrG-CSF-treated small (38.6% and 49.0%, respectively) follicles and 0 ng/mL hrG-CSF-treated medium (52.1%) follicles, and the cleavage rates were significantly (P < 0.05) higher in 10 ng/mL hrG-CSF-treated medium (76.3%) follicles than in all other groups. The blastocyst formation rates were significantly (P < 0.05) higher in 100 ng/mL hrG-CSF-treated small (31.2%) follicles than in 0 and 10 ng/mL hrG-CSF small (10.4% and 15.6%, respectively) follicles, and the 10 ng/mL hrG-CSF medium (45.7%) follicle was significantly (P < 0.05) higher than in all other groups. The total cell number in blastocysts from the 10 ng/mL hrG-CSF medium (106.5) follicles was significantly (P < 0.05) increased compared to 0, 10, 100 ng/mL hrG-CSF small (55.0, 73.7 and 59.5, respectively) follicles and 0, 100 ng/mL hrG-CSF-treated medium (82.5 and 93.5, respectively) follicles. After IVF, the blastocysts stage was significantly (P < 0.05) increased in 10 ng/mL hrG-CSF-treated medium (36.4%) follicles. Fertilization efficiency was significantly high in 100 ng/mL of small (29.1%) and 10 ng/mL of medium (44.0%) follicles. We also examined the Bcl2 and ERK2 transcript levels and found that they were significantly higher in the small and medium follicle treatment groups. In conclusion, these results indicate that hrG-CSF improve the viability of porcine embryos.

Global gene expression in granulosa cells of growing, plateau and atretic dominant follicles in cattle

BACKGROUND

The physiological state of the dominant follicle is important as it may be linked to the competence of the oocyte within. The objective of this study was to analyze, by transcriptomic analysis, the changes occurring in granulosa cells from dominant follicles at different phases of follicular growth.

METHODS

Granulosa cells were collected from slaughterhouse dairy cattle follicles with a diameter greater than 9 mm, and were classified at different phases of follicle growth based on flow cytometry profiles of DNA content after staining with propidium iodide. Three phases were identified based on the proportion of cells in -G1 (less than 2n DNA), G0-G1 (2n DNA) or S-M (more than 2n DNA) and follicles were thus allocated to the growing, plateau or atresia group. Between group analysis (BGA) showed clear segregation of the three groups, and the groups were contrasted against each other in a loop design to identify differently expressed genes. Ingenuity Pathway Analysis (IPA) was used to identify the functions and upstream regulators associated with the observed differently expressed genes.

RESULTS

Major differences were observed between the growth phases. Granulosa cells from follicles in the plateau phase had increased expression of TYRO3 and downregulation of JAM2 compared to growing follicles, supporting the idea of a shift from proliferation to differentiation. On the other hand, genes regulating the response to oxidative stress (VNN1) and angiogenesis (ANGPT2) were upregulated in granulosa cells from atretic follicles. While the predicted activated functions in cells at the plateau stage compared to cells at the growing stage included synthesis and transport of molecules, the predictions for atretic follicles relative to plateau ones included an increase in apoptosis and cell death.

CONCLUSION

Consistent with previous studies, these observations allowed us to match the presence of specific gene transcripts to a particular physiological status and consequently to classify follicles. The results also demonstrated that the plateau phase is not a simple 'in between' status between growth and atresia, as several characteristics are unique to this stage.