A retrospective model of oocyte competence: global mRNA and housekeeping transcripts are not associated with in vitro developmental outcome

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.

Validation of reference genes for gene expression studies in bovine oocytes and cumulus cells derived from in vitro maturation

Quantitative real-time PCR (qPCR) is a valuable tool for gene expression studies and it is necessary to choose an ideal endogenous reference gene for data normalization. This work studied a set of reference genes in oocytes and cumulus cells of COCs (Cumulus-Oocyte Complexes) that are suitable for relative gene expression analyses after in vitro maturation (IVM) in bovine. Immature COCs were collected from ovaries of Nelore cattle (Bos indicus) and submitted to IVM. MII oocytes and cumulus cells were subjected to RNA extraction, reverse transcription and preamplification of cDNA. The expression level of eight reference genes (ACTB, GADPH, B2M, H2AFZ, GUSB, HPRT1, PPIA, and TBP) was measured by real time PCR and analyzed by geNorm software. The gene stability measure (M) was calculated and the ideal number of reference genes (RGs) was determined by the V value (pairwise variation). For oocyte samples, two RGs were the ideal number for relative quantification: HPRT1 and B2M and for bovine cumulus samples four were indicated: HPRT1, PPIA, B2M, and TBP genes. The normalization of a non-reference target gene (SOD1) by these reference genes was shown to be considerably different from normalization by less stable reference genes. Our results strengthen the importance of choosing good normalizing genes in order to analyze gene expression under specific experimental conditions and we suggest the use of these RGs in oocytes and cumulus cells of bovine cattle in in vitro matured COCs.

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.

Messenger RNAs in metaphase II oocytes correlate with successful embryo development to the blastocyst stage

The mRNAs accumulated in oocytes provide support for embryo development until embryo genomic activation. We hypothesized that the maternal mRNA stock present in bovine oocytes is associated with embryo development until the blastocyst stage. To test our hypothesis, we analyzed the transcriptome of the oocyte and correlated the results with the embryo development. Our goal was to identify genes expressed in the oocyte that correlate with its ability to develop to the blastocyst stage. A fraction of oocyte cytoplasm was biopsied using micro-aspiration and stored for further expression analysis. Oocytes were activated chemically, cultured individually and classified according to their capacity to develop in vitro to the blastocyst stage. Microarray analysis was performed on mRNA extracted from the oocyte cytoplasm fractions and correlated with its ability to develop to the blastocyst stage (good quality oocyte) or arrest at the 8-16-cell stage (bad quality oocyte). The expression of 4320 annotated genes was detected in the fractions of cytoplasm that had been collected from oocytes matured in vitro. Gene ontology classification revealed that enriched gene expression of genes was associated with certain biological processes: 'RNA processing', 'translation' and 'mRNA metabolic process'. Genes that are important to the molecular functions of 'RNA binding' and 'translation factor activity, RNA binding' were also enriched in oocytes. We identified 29 genes with differential expression between the two groups of oocytes compared (good versus bad quality). The content of mRNAs expressed in metaphase II oocytes influences the activation of the embryonic genome and enables further develop to the blastocyst stage.

The mechanism of oocyte activation influences the cell cycle-related genes expression during bovine preimplantation development

The first cleavage divisions and preimplantation embryonic development are supported by mRNA and proteins synthesized and stored during oogenesis. Thus, mRNA molecules of maternal origin decrease and embryonic development becomes gradually dependent on expression of genetic information derived from the embryonic genome. However, it is still unclear what the role of the sperm cell is during this phase and whether the absence of the sperm cell during the artificial oocyte activation affects subsequent embryonic development. The objective of this study was to determine, in bovine embryos, changes in cell cycle-associated transcript levels (cyclin A, cyclin B, cyclin E, CDC2, CDK2, and CDK4) after oocyte activation in the presence or absence of the sperm cell. To evaluate that, in vitro-produced (IVP) and parthenogenetically activated (PA) embryos (2-4 cells (2-4C), 8-16 cells (8-16C) and blastocysts) were evaluated by real-time PCR. There was no difference in cleavage and blastocyst rates between IVP and PA groups. Transcript level was higher in oocytes than in IVP and PA embryos. Cleaved PA embryos showed higher expression of cyclin A, cyclin B, cyclin E, and CDK2 and lower expression of CDC2 when compared with that from the IVP group. At the time of activation, all transcripts were expressed less in PA than in IVP embryos, whereas at the blastocyst stage, almost all genes were expressed at a higher level in the PA group. These results suggest that in both groups there is an initial consumption of these transcripts in the early stages of embryonic development. Furthermore, 8-16C embryos seem to synthesize more cell cycle-related genes than 2-4C embryos. However, in PA embryos, activation of the cell cycle genes seems to occur after the 8- to 16-cell stage, suggesting a failure in the activation process.