Skip to content

Advertisement

You're viewing the new version of our site. Please leave us feedback.

Learn more

BMC Genomics

Open Access

Use of Both Cumulus Cells’ Transcriptomic Markers and Zona Pellucida Birefringence to Select Developmentally Competent Oocytes in Human Assisted Reproductive Technologies

BMC Genomics201516(Suppl 1):S9

https://doi.org/10.1186/1471-2164-16-S1-S9

Published: 15 January 2015

ABSTRACT

Background

Selection of the best oocyte for subsequent steps of fertilization and embryo transfer was shown to be the crucial step in human infertility treatment procedure. Oocyte selection using morphological criteria mainly Zona pellucida (ZP) has been the gold standard method in assisted reproductive technologies (ART) clinics, but this selection approach has limitations in terms of accuracy, objectivity and constancy. Recent studies using OMICs-based approaches have allowed the identification of key molecular markers that quantitatively and non-invasively predict the oocyte quality for higher pregnancy rates and efficient infertility treatment. These biomarkers are a valuable reinforcement of the morphological selection criteria widely used in in vitro fertilization (IVF) clinics. In this context, this study was designed to investigate the relationship between transcriptomic predictors of oocyte quality found by our group and the conventional morphological parameters of oocyte quality mainly the ZP birefringence.

Results

Microarray data revealed that 48 and 27 differentially expressed candidate genes in cumulus cells (CCs) were respectively overexpressed and underexpressed in the ZGP (Zona Good Pregnant) versus ZBNP (Zona Bad Non Pregnant) groups. More than 70% of previously reported transcriptomic biomarkers of oocyte developmental competence were confirmed in this study. The analysis of possible association between ZP birefringence versus molecular markers approach showed an absence of correlation between them using the current set of markers.

Conclusions

This study suggested a new integrative approach that matches morphological and molecular approaches used to select developmentally competent oocytes able to lead to successful pregnancy and the delivery of healthy baby. For each ZP birefringence score, oocytes displayed a particular CCs' gene expression pattern. However, no correlations were found between the 7 gene biomarkers of oocyte developmental potential and the ZP birefringence score. Further studies using larger lists of candidate markers are required to identify suitable genes that are highly correlated with the morphological criteria, and therefore able to reinforce the accuracy of oocyte selection and the effectiveness of infertility treatment.

BACKGROUND

Infertility is “a disease” that affects the reproductive apparatus making it unable to carry out a successful clinical pregnancy after at least one year of regular and unprotected intercourse [1]. The worldwide prevalence of this disease keeps rising and is estimated to 1 out of 6 couples of reproductive age. Consequently, an increase need for assisted reproductive technologies (ART) to help couples achieving their parenthood project has been recorded. It is important to highlight that regardless of the infertility cause, the oocyte quality is recognized as the main cornerstone amongst others of any successful ART treatment [2]. Consequently, numerous studies focused on the criteria allowing an efficient selection of high-quality oocytes, named also competent oocytes, for assisted reproduction procedure using both human and animal models. Such oocyte competence is defined as its ability to properly achieve full maturation at nuclear, cytoplasmic and molecular levels, to achieve successful steps of fertilization and early embryo development, and to ultimately lead to a viable and healthy offspring [35]. The cumulus-oocyte complex (COC) morphology and cytology are amongst the main parameters that have been used in ART clinics to select the best oocyte for subsequent steps of fertilization and early embryo development in human and livestock species [6, 7]. These parameters include the oocyte diameter, its cytoplasm granulation and the first polar body integrity [7, 8]. For its bordering cumulus cells (CCs), their appearance, compaction and number of layers have been also correlated to the oocyte developmental potential [912]. It is established that the oocyte communication with CCs are a prerequisite to acquire its full developmental competence. In fact, CCs were shown to enhance the oocyte quality by their physical presence (cell-cell adhesion structures), their morphology [12, 13], their metabolic and signal transduction functions [1417].

Zona pellucida (ZP), which is a filamentous matrix of well-structured and glycosylated glycoproteins surrounding the oocyte, is another morphological criterion of oocyte selection. This matrix is made of three proteins encoded by three different genes ZP1, ZP2 and ZP3 [18, 19]. An additional protein expressed by a ZP1-paralogue gene, the ZP4, was also reported in human [2022]. ZP genes’ transcription and glycoprotein synthesis are reported in both the oocyte and CCs of most mammals. This ZP acts as an oocyte’s coat since the oocyte growth and maturation, remains during the early embryo development until the blastocyst hatching [23, 24]. Its thickness increases with the oocyte maturation progression to reach around 17 µm at the human MII stage [25]. ZP is also marked by layers of long filaments with repetitive protein heterodimers ZP2-ZP3 cross-linked by ZP1 homodimers [26, 27]. ZP proteins’ glycosylation is maximal in the outer layers and decreases progressively to a more compact and less porous inner region [24, 28]. In human and mice, ZP3 is crucial for sperm recognition and fertilization [27, 2931]. A positive correlation between ZP Thickness and uniformity, and the human oocyte quality was reported in several studies [9, 3234]. That’s why many human IVF clinics used the ZP thickness and birefringence as a morphological criterion to select oocytes with higher ability to achieve successful pregnancy [32, 3436].

Despite their contribution to a relative improvement of IVF outcomes, these morphological and microscopic-based selection approaches have limitations in terms of accuracy, objectivity and constancy. They rely more on the embryologist/clinician experience since the molecular events underlying the oocyte competence acquisition and their manifestation on its morphology are still poorly understood [37, 38]. The alternative of transferring more than one embryo to balance such weakness has led to a significant increase of multiple pregnancies rates, premature births, mother health problems and associated costly health treatments [3941]. Therefore, the qualitative and quantitative increase of ART outcomes and the reduction of associated health care costs for both the mother and the offspring are still challenging [42, 43].

The rapid development of high-throughput OMICs technologies has been a precious adjunct to ART in order to develop more effective strategies to improve the selection of high-quality oocytes using non-invasive and quantitative approaches, allowing therefore more effective elective single embryo transfer (eSET) procedure [44, 45]. In this context, some promising strategies based on indirect prediction of oocyte competence/quality using the follicular cells in the oocyte neighborhood mainly CCs were recently suggested. The focus on CCs is mainly due to their continuous bidirectional molecular interplay with the oocyte shown to be crucial for her subsequent developmental competence. Therefore, CCs are thought to reflect the oocyte competence level. Several molecular biomarkers using follicular cells and mainly CCs have been suggested recently by our group [4, 4651] and elsewhere [5255]. We assume that these molecular biomarkers are a valuable reinforcement to the morphological selection criteria accumulated by the clinicians’ cumulative experience and used in IVF clinics worldwide. In this context, this study was designed to study possible correlations between some molecular predictors of oocyte quality and the conventional morphological parameters of oocyte quality mainly ZP birefringence. We expect to report a positive correlation between the ZP birefringence score and transcriptomic markers of oocyte competence expressed in CCs. Such pioneering approach based on combining the OMICs tools with established oocyte quality parameters mainly the ZP birefringence should lead to an integrative strategy that will accurately predict the oocyte developmental potential. A clinical tool that use both molecular biomarkers and ZP birefringence to efficiently select high quality oocytes in IVF clinics will allow the selection of the best oocyte(s), a successful elective single embryo transfer (eSET) and consequently more effective infertility treatment.

RESULTS

Gene expression analysis

RNA hybridizations were performed on a custom-made microarray platform enriched in CCs, granulosa cells (GCs), oocytes and early embryos ESTs using suppressive subtractive hybridization (SSH) of cDNA sequences as described elsewhere [56, 57]. Results of microarrays hybridization of CCs of ZGP against those of ZBNP showed that 48 and 27 candidate genes were respectively overexpressed and underexpressed in the ZGP compared to the ZBNP (ZGP/ZBNP ratio ≥ 2; FDR=5%) (see tables S1 and S2, additional files 1 and 2). Interestingly, at least three (3) positive oocyte biomarkers from the (ZGP vs ZGNP) reported in our previous study [49] appeared also in the overexpressed gene list of the (ZGP vs ZBNP). These biomarkers of oocyte quality are PSMD6, CALM1 and NRP1.

Cumulus cells-induced molecular, cellular and physiological functions

The analysis of the CCs differentially expressed genes (48 overexpressed; 27 underexpressed) (see tables S1 and S2, additional files 1 and 2) showed that several candidates are associated to key functions related to the reproductive function at both molecular and cellular levels (Table 1). Around 23% and 19% of the differentially expressed genes were respectively associated to the reproductive function and inflammatory-like response. These results confirmed the major role of CCs in the final maturation of the oocyte and the ovulation process, known by its inflammatory-like aspect. At the molecular and cellular levels, the time of CCs removal prior to intracytoplasmic sperm injection (ICSI) is marked by an important expansion of extracellular matrix (ECM), active molecular and signaling events to prepare ovulation, steroidogenesis and cell differentiation. That’s why important functions associated to cell morphology (≈ 30% of candidate genes), cell signaling (≈ 17 % of candidate genes) and cellular assembly (≈ 12 % of candidate genes) were activated (Table 1). Conversely, the cell cycle function (≈ 7% of candidate genes) was under-represented since the CCs are in a differentiation stage characterized by more apoptosis rather than proliferation (Table 1).
Table 1

Main molecular, cellular and physiological functions triggered by CCs in ZGP versus ZGNP. These functions are ranked by p-value and their proportion of differentially expressed genes (both overexpressed and underexpressed) according to the IPA software.

Molecular, cellular and physiological functions

p-value

No. genes

Proportion of differentially expressed genes involved (%)*

Physiological functions / conditions

Reproductive and embryonic development

4.13E-05 - 2.95E-02

17

22.67

 

Organ morphology

5.92E-06 - 2.95E-02

16

21.33

 

Inflammatory-like response

2.86E-04 - 2.95E-02

14

18.67

Molecular & cellular functions

Cell Morphology

5.92E-06 - 2.95E-02

22

29.33

 

Small molecule biochemistry & signaling

1.09E-04 - 2.59E-02

13

17.33

 

Cellular assembly and organization

4.13E-05 - 2.95E-02

9

12.00

 

Cell cycle

1.37E-04 - 2.59E-02

5

6.67

*: The proportion of differentially expressed genes included both over- and under-expressed genes

Real time PCR analysis

In addition to the validation of the microarray data, quantitative real time PCR was performed on the initial biological samples to assess the expression levels of the 7 gene biomarkers found in our previous study to be overexpressed in ZGP group compared to its ZGNP counterpart [49]. This QPCR validation was achieved on CCs of oocytes with ZGP versus those with ZBNP (successful pregnancy and HZB versus pregnancy failure and LZB). Among the 7 biomarkers associated to both HZB and successful pregnancy reported previously, 5 genes (out of 7) were also significantly different between the two groups studied (ZGP versus the ZBNP) (i.e. 71.4 % of biomarkers confirmed), which is a further quantitative confirmation of the validity of our biomarkers on separate biological samples. These candidates are NRP1 (p = 0.003), CALM1 (p = 0.005), DPP8 (p = 0.006), UBQLN1 (p = 0.025) and PSMD6 (p = 0.048) (Figure 1). However, the association analysis between these biomarkers and ZP birefringence revealed that they were not significantly different between ZGNP and ZBNP (p > 0.05; Table 2).
Figure 1

Real-time PCR analysis of differentially expressed genes in individual CCs of ZGP group versus ZBNP. Gene candidates were ranked according to their p-values, which were determined following a T-test analysis achieved on normalized data at α = 0.05

Table 2

Quantitative real-time PCR comparative analysis of differentially expressed genes in individual CCs among the three ZP quality groups (ZGP versus ZGNP versus ZBNP). Candidates with different letters are significantly different following a F-test analysis at α = 0.05

Gene

CALM1

DPP8

HIST1H4C

NRP1

PSMD6

TOM1

UBQLN1

ZP quality groups

ZGP

a

a

a

a

a

b

a

 

ZGNP

b

b

b

b

b

a

b

 

ZBNP

b

b

a

b

b

a,b

b

Gene networks’ analysis

Both overexpressed and underexpressed candidate genes of ZGP vs ZBNP groups were mapped to their potential molecular, cellular and biological functions using the Ingenuity Pathways Analysis (IPA) software [58]. Four (4) main networks were selected for further discussion and a comprehensive analysis of differentially expressed genes in CCs collected prior to ICSI and their association with ZP and follicular function. These top-ranked networks were made of differentially over- and under-expressed genes in CCs and focused on: i) cell morphology, ii) intracellular signaling, iii) immune-like response, iv) steroidogenesis and v) apoptosis (Figures 2, 3, 4).
Figure 2

Main candidate genes differentially expressed in CCs and involved in both cell and ZP morphologies as revealed by IPA software. (Red): overexpressed genes. (Green): underexpressed genes. (Gray): ZP genes. Genes are CASP9 (caspase 9), CCNA2 (cyclin A2), RAN (RAN, member RAS oncogene family), CDKN1A (cyclin-dependent kinase inhibitor 1A), CYP11A1 (cytochrome P450, family 11, subfamily A, polypeptide 1), GLB1 (galactosidase, beta 1), HBA1/HBA2 (hemoglobin, alpha 1/ hemoglobin, alpha 2), HBB (hemoglobin, beta), HDAC2 (histone deacetylase 2), HSP90B1 (heat shock protein 90, beta, member 1), LAMB1 (laminin, beta 1), PRDX2 (peroxiredoxin 2), SGCB (sarcoglycan, beta), TIE1 (tyrosine kinase with immunoglobulin-like and EGF-like domains 1), ZP1 (zona pellucida glycoprotein 1 (sperm receptor)), ZP2 (zona pellucida glycoprotein 2 (sperm receptor)), ZP3 (zona pellucida glycoprotein 3 (sperm receptor)), and ZPBP (zona pellucida binding protein).

Figure 3

Gene network of differentially expressed candidates in CCs involved in both cell morphology and intracellular signalling pathways cross-talk. (Red): overexpressed genes. (Green): underexpressed genes. (White): other genes. Overexpressed genes are PRDX2 (peroxiredoxin 2), HBA1/HBA2 (hemoglobin, alpha 1/ hemoglobin, alpha 2), HBB (hemoglobin, beta), RAN (RAN, member RAS oncogene family), CASP9 (caspase 9), HDAC2 (histone deacetylase 2), HSP90B1 (heat shock protein 90, beta, member 1), TIE1 (tyrosine kinase with immunoglobulin-like and EGF-like domains 1), GLB1 (galactosidase, beta 1), CDKN1A (cyclin-dependent kinase inhibitor 1A) and RGS3 (regulator of G-protein signaling 3). Underexpressed genes are CCNA2 (cyclin A2) and CYP11A1 (cytochrome P450, family 11, subfamily A, polypeptide 1). Other genes are ADRB (Adrenoreceptor, Beta), Ck2 (casein Kinase II), PI3K (phosphatidylinositol-4,5-bisphosphate 3-kinase), PKC (protein Kinase C), MAPK and ERK1/2 (mitogen-activated protein kinase family), HSP70 (heat shock protein 70), RNA polymerase II, Insulin, AP1 (activator protein-1), CD3 and LH (luteinizing hormone).

Figure 4

Gene network of several differentially expressed CCs genes involved in both apotosis and anti-inflammatory-like response. (Red): overexpressed genes. (Green): underexpressed genes. (White): other genes. Overexpressed genes are CCAR1 (cell division cycle and apoptosis regulator 1), RGS3 (regulator of G-protein signaling 3), PRDX2 (peroxiredoxin 2), HSP90 / HSP90B1 (heat shock protein 90, beta, member 1), PANX1 (pannexin 1), CASP9 (caspase 9), CDKN1A (cyclin-dependent kinase inhibitor 1A), PSIP1(PC4 and SFRS1 interacting protein 1), RAN (RAN, member RAS oncogene family) and HDAC / HDAC2 (histone deacetylase 2). Underexpressed genes are FOSB (FBJ murine osteosarcoma viral oncogene homolog B), FDPS (farnesyl diphosphate synthetase), TAX1BP1 (Tax1 (human T-cell leukemia virus type I) binding protein 1), HSPA8/ HSP70 (heat shock protein 70 kDa, protein 8), CCNA2 (cyclin A2) and ENO1 (Enolase 1, alpha). The other genes of the network are ERK1/2 (Extracellular signal-regulated kinase 1 /2, mitogen-activated protein kinase family), HSP27 (heat shock 27kDa protein), BCR (complex) (B-cell receptor complex), PP2A (protein phosphatase type 2a), Histone H1, PKG (protein kinase G), LDL (low density lipoprotein) and IL1 (interleukin 1).

DISCUSSION

Study approach

This study was designed to investigate potential relationships/association between the oocyte quality/competence reflected by its developmental potential and the ZP morphology. For the oocyte competence, it was predicted using transcriptomic markers differentially expressed in its CCs through a comparative analysis of gene expression patterns between ZGP group versus ZGNP in our previous study [49]. Concerning the ZP morphology, it was analyzed using a polarizing microscope (Polscope) linked to a robot-like micromanipulation system (to ensure zona 3D visualisation and scoring) and a suitable software that automatically and objectively (user-independent) measure both ZP density and uniformity. The birefringence score is determined based on approximately 180 measurements that cover all the ZP and was shown to be correlated to the oocyte quality [32, 35, 59, 60]. Taken together, the interesting progress in both ZP birefringence scoring accuracy and the non-invasive and quantitative biomarkers of oocyte quality were behind this original study. The goal is to investigate possible correlations between oocyte developmental competence markers and ZP birefringence in order to assess if these two approaches are additive, and hence able to strengthen each other and enhance the selection of highly competent oocytes procedure. Such positive association once confirmed is expected to improve the infertility treatment approaches through the suggestion of an integrative method based on confirmed molecular predictors that reinforce the morphological selection based on ZP birefringence and the embryologist experience.

Cumulus cells’ molecular and cellular gene networks

Microarray data analysis revealed a non exhaustive list of 75 differentially expressed genes between ZGP versus ZBNP treatment groups (see tables S1 and S2, additional files 1 and 2) that reflected the behavior and roles of both CCs and ZP prior to ovulation. The analysis of the molecular, cellular and physiological processes triggered by both overexpressed and underexpressed genes revealed key gene networks. The molecular and cellular functions/ networks having the highest scores were selected for further analysis.

Cell morphology

It was interesting to report that around 30% of differentially expressed genes in CCs were involved in the regulation of cell morphology (Table 1; Figure 2). These findings showed that important ultrastructural changes occurred in cumulus cells during final maturation and ECM expansion through microtubule and microfilaments organization (RAN, PRDX2), cell differentiation (CDKNA1, HDAC2, GLB1, CCNA2, PRDX2) and intracellular signaling (CDKNA1, TIE1, RGS3) (Figure 3). Similar morphological and metabolic changes in CCs due to the remodeling of their cytoskeleton [61, 62], the transzonal projections number and signaling [29, 63] and the disruption of CCs shapes and cell-cell contacts [64, 65] were also reported. In several mammalian species, it was documented that these CCs intracellular events are essential to support oocyte final maturation, to convey COC some signaling factors and viscoelastic properties that ease its release during ovulation as well as to contribute in sperm recognition and fertilization [6669]. ZP is also an important structure involved in the organization of its own matrix as well as the oocyte and CCs morphology through dynamic arrangements of its structure and microfilament networks (Figure 2; [25, 70]). In this study, it was used as a crucial distinctive parameter between the two treatment groups (ZGP vs ZBNP) which may further explain the upregulation of the cell and organ morphology functions (Table 1).

Intracellular signaling

In addition to cell morphology, CCs are characterized by a variety of signaling pathways that finely regulated their functions. Such signaling events are governed by gonadotropins [7173], the oocyte-secreted factors [7476] as well as the intra-ovarian endocrine and paracrine environment [77, 78]. These signaling pathways involve several players including second messengers’ receptors, transcription factors and kinases as illustrated by figure 3. It is important to note that oocyte competence acquisition and ovulation are still poorly understood and this high cross-talk between key signaling pathways mainly kinases (figure 3) is expected. In fact, CCs is considered as an interface between the oocyte and its surrounding environment through functional signaling cascades involving the pathways of these kinases (ERK, p38MAPK, PKC, PKA and PI3K) known to be crucial for mammalian oocyte quality, CCs expansion, ovulation and even early embryo development [17, 7883].

Inflammatory-like response

CCs were retrieved prior to ovulation which is a complex mechanism that allows the rupture of the follicle and the ovarian epithelium to release the COC in the fallopian tube. It is a crucial differentiation step in the reproductive function marked by a local inflammatory reaction that involves different signals from the blood supply in the theca cells as well as steroidogenic and proteases factors produced within the follicle. These signals triggered the CCs expression of many genes associated to an inflammatory-like condition (around 19% of total differentially expressed genes; Table 1) including mainly the underexpression of anti-inflammatory regulators (TAX1BP1, MT-CYB, ENO1) and heat shock proteins (HSPA8, HSP90B1, HSP70); and the overexpression of immune (HLA-DRA, SCARA5, PRDX2) and cell signalling (RGS3, RAN, PSPI1) factors ( (Figure 4, see tables S1 and S2, additional files 1 and 2).

These inflammatory-related genes were suggested as protectors of the oocyte and CCs in an inflammatory and pro-apoptotic environment during ovulation [84, 85]. Moreover, this inflammatory-like response is considered as prerequisite for ovulation and subsequent fertilization [84, 8688]. Recent genome-wide association studies (GWAS) showed that CCs gained some immune and neuronal functions required during ovulation in many mammal species including human [3, 89, 90], mouse [85, 91], rat [92] and bovine [56]. Overall, these findings support a significant role of CCs in ovulation.

Apoptosis

At this preovulatory period, CCs have reached full expansion marked by weaker cell-cell connections and increased level of apoptosis [93, 94]. This apoptotic phenotype was clearly reflected in our study by several differentially expressed genes in CCs and involved in pro-apoptotic pathways as heat shock proteins (HSPA8, HSP90B1, HSP70), caspases (CASP9), cell cycle under-regulation and cell differentiation (FOSB, CCAR1, CCNA2, CDKNA1, HDAC2, PRDX2, PANX1) as well as cell death-related intracellular signaling and transport (FDPS, CDKNA1, RGS3, TAX1BP1, RAN, PSPI1) (Figure 4). Such pro-apoptotic gene pattern was associated to an obvious underexpression of cell cycle genes (less than 7% of differentially expressed candidate genes; Table 1). This apoptosis increase was also reported in follicular cells (cumulus and granulosa) in primates and cattle [6, 95, 96]. Interestingly, a slight apoptosis in the outer layers of CCs is considered as a morphological criterion that is positively correlated with the oocyte developmental competence [79, 93, 95, 97, 98]. Additionally, this apoptotic behavior is expected at this stage since CCs have achieved most of their main roles required to both ovarian function and oocyte quality. We assume that these apoptotic signs are an intrafollicular message about the oocyte preparedness for ovulation or atresia. This apoptotic condition may also amplify the immune and inflammatory-like response that prepares the ovulation process.

Steroidogenesis

It is established that follicular growth is governed mainly by gonadotropins and steroids. The nearness onset of ovulation is marked by a fine regulation of steroidogenesis to increase estradiol levels and initiate progesterone synthesis [99, 100]. Several genes involved in steroidogenesis (around 23%) were differentially expressed as shown by the top position of the reproductive function (Table 1). This CCs active hormonal activity and early luteinisation was also clearly reflected by several differentially expressed genes in CCs (see tables S1 and S2, additional files 1 and 2) and confirmed in previous reports [95, 101104].

Zona pellucida and ovarian function

According to CCs gene expression data, we confirmed herein the complexity of molecular and signalling events that drives the final steps of oocyte maturation, its interplay with CCs and ovulation. Since this study is comparing two treatment groups with different ZP birefringences, it is important to highlight that this oocyte’s coat remain in the oocyte immediate vicinity since the oocyte growth and maturation and remains during the early embryo development until hatching [23, 24]. ZP structure is based on filaments of ZP2-ZP3 heterodimers proteins, linked in some binding sites by ZP1[26]. It is involved in vital functions related to ovarian function and oocyte quality as summarized in figures 2 and 5 following gene function analysis using IPA software. In fact, ZP is a filamentous matrix located in the perivitelline space around the oocyte and composed of well-structured and highly- glycosylated glycoproteins [24, 28]. ZP genes transcription and their glycoprotein synthesis are reported in both the oocyte and CCs of most mammals and involved in maintaining both COCs and ZP morphologies (Figures 2 and 4), prevent polyspermy and protect early embryo development until the blastocyst hatching [25].
Figure 5

Main ZP genes’ functions related to ovarian function and oocyte quality following a functional analysis using the IPA software. ZP1 (zona pellucida glycoprotein 1 (sperm receptor)), ZP2 (zona pellucida glycoprotein 2 (sperm receptor)), ZP3 (zona pellucida glycoprotein 3 (sperm receptor)), and ZPBP (zona pellucida binding protein).

It is important to highlight that according to our gene expression data of this study and elsewhere [49], it appears that each ZP group has its particular gene expression pattern that is dependent on both ZP birefringence and oocyte quality.

Association ZP birefringence versus oocyte transcriptomic markers

Our data confirmed that CCs compartment is the site of specific signalling and gene expression events which support and reflect the oocyte’s quality. ZP was also showed to be crucial in key functions associated to the oocyte maturation, fertilization and early embryo development. Taken together, CCs and ZP look to be key sites of a well space and time-coordinated sequence of molecular events that determine the oocyte quality and its subsequent developmental competence. The main objective of this study is to study the possibility of combination of these two non-invasive parameters (molecular markers and ZP birefringence) to enhance the effectiveness of high-quality oocyte selection and reinforce the already used morphological criteria. To do this, we have analyzed a list of CCs genomic markers of oocytes having high ZP birefringence (HZB) and were able to achieve successful pregnancy as described in our previous study [49]. Correlation analysis of QPCR results of each one of these transcriptomic markers and both ZP birefringence and pregnancy outcomes revealed that 71.4 % (5 out of 7) of these biomarkers were confirmed using the ZBNP group. This finding is a further quantitative confirmation of the validity of our biomarkers identified on the (ZGP versus ZGNP) study [49] and using separate biological samples. These candidates are NRP1 (p = 0.003), CALM1 (p = 0.005), DPP8 (p = 0.006), UBQLN1 (p = 0.025) and PSMD6 (p = 0.048) (Figure 1). Such results are in line with previous reports confirming the usefulness of CCs transcriptomic markers for selection of the best human oocytes for subsequent fertilization and embryo transfer [42, 47, 49, 98, 105], increasing therefore the fertility treatment success rates.

Regarding potential additive effects between these biomarkers and the ZP birefringence, this correlation, if confirmed, could offer an interesting integrative approach based on the combination of morphological (ZP) and molecular (biomarkers) criteria allowing an accurate and non invasive selection of high competent oocytes. Thus, more reliable clinical tools/kits using both parameters could be developed for IVF clinics. However, our data revealed that most of the gene biomarkers were not significantly expressed between the two groups with failed pregnancy (ZGNP vs ZBNP; p> 0.05) except for HIST1H4C (Table 2). Therefore and despite having two opposite ZP birefringence scores (HZB versus LZB), few or no gene biomarkers differences were found between ZGNP versus ZBNP. According to these data, it looks that the ZP morphological phenotype is not directly associated to the selected transcriptomic markers of oocyte competence used in this study. Although the ZP birefringence was positively correlated to oocyte developmental competence, this selection criterion did not revealed significant differences at the gene expression levels once validated by QPCR. As a result, it was difficult to establish a positive correlation/association between the ZP morphology and the molecular marker of oocyte competence expressed in CCs.

Although these data are preliminary and require more investigations, several reasons could explain the lack of correlation between these two parameters. In fact, only 7 gene biomarkers have been used to assess such correlation which might be too small to cover potential genes directly affected by the ZP morphology and/or functions. Larger studies using whole genome microarrays and increased number of selected biomarkers are required to profoundly explore this relationship. Moreover and as mentioned before, the ZP formation started at the secondary follicular stage and is achieved at the preovulatory stage. Therefore at the time of COCs collection for ICSI, low correlation between the ZP birefringence and CCs gene expression may be expected despite their both involvement in subsequent fertilization. Another study by Dr. Montag group [106], reported a differential gene expression gradient within CCs when moving from the very differentiated corona radiata cells (CRCs) at the ZP vicinity to the outer layers. Therefore and due to the high similarities of the compared tissues, only corona CRCs may show a transcriptomic difference associated to ZP birefringence. This assumption is further supported by possible involvement of the CRCs in ZP protein biosynthesis. Therefore, the analysis of the whole CCs in our case may dilute possible transcriptional difference in CRCs with opposite ZP birefringence scores. The possibility that ZP birefringence is an independent parameter that does not unlikely correlate with the genomic predictors of oocyte competence might be also considered.

It is important to highlight that other promising technologies as time-lapse embryo monitoring (TLEM) and preimplantation genetic screening (PGS) can be also combined to our approach [107, 108]. Such complementary technologies will even strengthen more our approach through the combination of high quality oocyte selection (suggested herein) and morphologically good quality embryos (TLEM and/or PGS) would offer a very accurate and integrative selection approach that will allow efficient infertility treatments and higher pregnancy successful rates.

CONCLUSION

Our results confirmed that ZP properties’ variation is associated to a CCs gene expression difference. We also confirmed 5 gene biomarkers of oocyte developmental competence using different biological samples. Surprisingly, no correlations between our 7 CCs gene biomarkers of oocyte developmental potential and the ZP birefringence score. It looks that the ZP morphology is associated to a transcriptomic gene pattern that is not directly related to these biomarkers. These findings highlighted the complexity of the molecular events underlying the developmental competence pathway. Hence, further studies using larger lists of candidate markers are required to identify suitable genes that are highly correlated with the morphological criteria, and therefore able to reinforce the accuracy of oocyte selection. Ultimately, the idea of finding ZP-related biomarkers of oocyte quality would reinforce the accuracy of oocyte selection. Additive or correlated oocyte selection criteria should strengthen each other offering therefore a reliable prognostic tool of the pregnancy outcome. The improvement of oocyte selection procedures would increase the efficiency of infertility treatment as well as the pregnancy outcomes. Additional promising embryo selection technologies as TLEM and PGS are expected to reinforce our current approach through a later selection at the early embryo stage. Therefore, efficient single embryo transfer could be envisioned.

MATERIALS AND METHODS

Ethics statement

Patients involved in this study were recruited from the IVF clinic of the University of Bonn Medical School after signing a written informed consent. All procedures and protocols used in this study were approved by the institutional review boards of both the University of Bonn Medical School, Germany; and Laval University, Canada.

Patients

Seven consenting patients (n=7) were meticulously selected for ovarian stimulation and ICSI (intracytoplasmic sperm injection). Patient selection was based on diagnosis of male factor infertility induced due to low quality sperm (low concentration, morphology or motility) according to World Health Organization (WHO) guidelines [109]. To proper infertility diagnosis and prevent any potential bias, all patients were first subjected to extensive andrological, gynaecological and cytogenetic examinations. ICSI was the recommended procedure for all the patients included in the study. Based on ZP birefringence, CCs of each individual oocyte were divided into two classes: i) HZB: includes 8 CCs of oocytes with High Zona Birefringence; and ii) LZB: includes 3 CCs of oocytes with Low Zona Birefringence.

IVF and culture media

All media used for oocyte retrieval, denuding, gametes’ washing, ICSI treatment and subsequent embryo culture (Fertilization, Cleavage, Gamete, PVP, Hyaluronidase, Culture oil) were of pharmaceutical, embryo-tested grade and free of phenol red. Unless mentioned, all embryo culture media were provided by Cook Company, Brisbane, Australia.

Ovarian stimulation and cumulus-oocyte-complexes recovery

A gonadotropin releasing hormone agonist (GnRHa)-based stimulation regimen using triptorelin acetate (Decapeptyl (0.1 mg/day), Ferring, Germany) and HMG (HMG; Menogon, Organon) / FSH (Gonal-F, Serono, Germany) was given to patients as described before [49]. FSH/HMG doses (average of 225 IU) were individually adjusted for each patient using ultrasound monitoring mainly follicular diameter and estradiol. A dose of 10,000 IU of hCG (human chorionic gonadotrophin) was used to induce ovulation and COCs were transvaginally punctured 36 to 38 h later.

Cumulus cells removal and ZP birefringence assessment

Collected COCs were collected, washed and individually incubated at warm temperature using pre-mixed gas with low oxygen (6% CO2, 5% O2, 89% N2) at 37°C. CCs were manually collected using a sterile scalpel in HEPES-buffered medium under oil and put at -80°C. Cumulus-free MII oocytes were selected following a hyaluronidase treatment and thereafter individually incubated in 5-μL droplets of fertilization medium covered with mineral oil at 37°C for subsequent ZP imaging. The analysis of ZP birefringence was automatically done using an automatic module Octax polairAide™ (Octax ICSI Guard™, OCTAX Microscience GmbH, Altdorf, Germany) connected to a polarization imaging software (OCTAX Eyeware™) as reported elsewhere [32]. Zona birefringence visualization and scoring were automatically and non-invasively done at 180 points under a heating platform allowing 37.0 ± 0.5°C in the medium droplet during micromanipulation and microscopic observation. ZP images were captured and used to determine the ZP scores which reflected the uniformity and the intensity of the ZP around the oocyte [32]. While LZB MII oocytes are characterized by an irregular and/or low birefringence distribution, the HZB MII ones had uniform and high intensity birefringence ZP (Figure 6).
Figure 6

ZP birefringence classes as shown by the polscope. (A): is a low zona birefringence MII oocyte (LZB); (B): is a high zona birefringence MII oocyte (HZB).

Intracytoplasmic sperm injection (ICSI) and embryo culture

Prior to ICSI, patients were subjected to andrological, gynaecological and cytogenetic examination. The spermatozoa ejaculate was washed, centrifuged, suspended in fertilization media and stored in a CO2 incubator. Following ZP imaging, all MII oocytes with know ZP birefringence and good morphology from patients’ group selected for their male factor-induced infertility received ICSI procedure according to routine standard protocols and using a spermatozoon with good morphology and motility as described elsewhere [49]. Following ICSI, selected zygotes with successful fertilization (presence of two pronuclei (2PN) of equal size in the center of the ooplasm) were individually cultured in a Minc benchtop incubator at 5% O2, 6% CO2, 98% N2 until day 3.

Embryo transfer (ET) and pregnancy assessment

For each patient and as per law, the two fertilized oocytes with the highest ZP birefringence scores were selected for transfer and the supernumerary oocytes were cryopreserved. The priority for embryo transfer was given to early embryos produced by HZB MII oocytes whenever available; otherwise those coming from LZB ovules were transferred. ET was performed using a transvaginal intrauterine Sydney IVF catheter as described previously [32]. Pregnancy was checked first using a hCG test at day 14 post transfer and confirmed by ultrasounds to detect the presence of gestational sacs and a positive heart beat (viable embryo) at 5 weeks of pregnancy.

Experimental groups

Based on ZP birefringence and subsequent pregnancy results, CCs coming from each individual oocyte were considered as biological replicates and were split into two main experimental groups (Figure 7):
Figure 7

Study’s experimental design

  • ZGP (Zona Good Pregnant) group: CCs of HZB oocytes that lead to successful pregnancy.

  • ZBNP (Zona Bad Non Pregnant) group: CCs of LZB oocytes and pregnancy failure.

While the ZPG includes 8 CCs and is considered as the positive treatment, the ZBNP one acted as a negative control for ZP relationship with pregnancy.

RNA extraction and amplification

Total RNA from each CC sample were performed using the PicoPure RNA Isolation Kit (Arcturus, Molecular Devices Analytical Technologies, Sunnyvale, CA) according to the user manual guidelines and eluted in 30 μL of Elution buffer (EB) provided in the kit. Good RNA quality and concentration were assessed and confirmed using the Agilent 2100 bioanalyzer (Agilent Technologies, Waldbronn, Germany) according to the manufacturer’s protocol. For each experimental group, two pools (biological replicates) of 10 ng of total RNA was used for linear amplification of messenger RNA (mRNA) using the 2-round in vitro transcription (IVT) following the instructions of the RiboAmpplus RNA Amplification kit (Arcturus, Molecular Devices Analytical Technologies) as described in the manufacturer's manual. Amplified mRNA yield was eluted in 30 μL of RNA eluted buffer (RE) and quantified by spectrophotometry at 260 nm using the NanoDrop ND-1000 (NanoDrop Technologies, Wilmington, DE).

Microarray hybridizations

Amplified messenger RNA pools of each group (ZGP vs ZGNP) were labelled using the Universal Linkage System (ULS™) aRNA Fluorescent Labelling Kit (KREATECH Biotechnology, Amsterdam, The Netherlands) according to the manufacturer’s instructions and used in a dye-swap design (Figure 7). Prior to hybridization, an equimolar mixture of the two labelled probes of each ZGP group’s pool and its counterpart in ZGNP group was prepared based on labelling dosage. Hybridization was achieved using a partial custom-made array enriched with transcripts associated to good quality oocytes as described elsewhere [49, 57]. Hybridizations were done in the ArrayBooster using the Advacard AC3C (The Gel Company, San Francisco, CA) for 18 h at 50°C using Slide Hyb#1 (Ambion, Austin, TX). The slides were washed successively in (2X SSC/0.5% SDS), (0.5X SSC/0.5% SDS) and 1X SSC buffers, then spin-dried and immediately used for subsequent scanning and analysis.

Microarray data analysis

Hybridized microarray slides were scanned using the VersArray ChipReader 3.1 System (Bio-Rad, Mississauga, Canada) and analyzed using the ArrayPro Analyzer software (Media Cybernetics, Bethesda, MD). Raw microarray data were first Loess-normalized and corrected for background as described before [56, 57]. Ratio of net fluorescence intensities of our dye-swap experiments between positive (ZGP; pregnant) and negative (ZBNP; non-pregnant) groups was assessed using Array Analysis Tool (Baltimore, MD) developed at the National Institute on Aging (NIA) [110] at FDR=5% and a minimum cut-off limit of 2.25. Given that each clone was printed twice on each microarray slide, these two additional technical sub-replicates were considered during data statistical analysis. Two lists of respectively over-expressed and under-expressed gene candidates with more than two-fold change ZGP / ZBNP were generated for subsequent analysis and data validation.

Gene networks’ analysis

The analysis of gene networks and the molecular pathways triggered by the CCs differentially expressed genes was performed using the QIAGEN’s Ingenuity Pathway Analysis software (IPA) [58]. Briefly, the candidate genes were uploaded with their official name and fold change into the IPA. Such software allowed potential connections and downstream signalling pathways between candidate genes and previous studies to generate and rank causal and/or hierarchical network analysis that may explain the observed gene expression changes in CCs of ZGP vs. ZBNP groups. Such approach was helpful to demystify the molecular pathways and cellular functions associated to ZP birefringence and their potential relationship with the oocyte developmental competence.

Quantitative real-time PCR of oocyte quality markers

Selected gene candidates mainly the molecular markers of high developmental capacity oocytes reported in a previous study [49] were validated on the original samples of ZGP, ZGNP and ZBNP groups using quantitative real-time PCR. Briefly, equal amounts of total RNA were taken from each replicate on individual CCs of each patient group, denatured and reversed transcribed using the SensiScript reverse transcriptase kit (Qiagen, Mississauga, ON, Canada). Real-time PCR was performed on the selected candidates in LightCycler capillaries (Roche Applied Science, Mannheim, Germany) using the LightCycler FastStart DNA Master SYBR Green I kit (Roche) as well as custom sets of primers reported elsewhere [49] according to the manufacturer’s guidelines. Three housekeeping/control genes ACTB (β-actin), GAPDH, and PPIA were quantified and used in GeNorm Normalization [111, 112]. The two housekeeping genes (ACTB and PPIA; P > 0.05) were the most stable and therefore used for QPCR data normalization.

Association between transcriptomic markers versus ZP birefringence

In order to assess the relationship between the 7 CCs transcriptomic biomarkers of high developmental potential oocytes associated to successful pregnancy and the zona pellucida birefringence, association study between the levels of expression of such biomarkers in 3 patient groups (ZGP vs ZBNP vs ZGNP) having different ZP birefringences and pregnancy outcomes was performed. Therefore and for each gene biomarker, its gene expression levels revealed by QPCR in CCs of the three patients’ groups were compared and analyzed based both the pregnancy output and ZP birefringence score.

Statistical analysis

Statistical analysis of QPCR results was done by ANOVA followed by Fisher’s protected least significant difference (LSD) test at α=0.05 and using the GraphPad Prism 5 software (GraphPad Software, San Diego, CA). Similarly, association analysis between oocyte quality biomarkers and ZP birefringence was performed for each gene biomarker using the QPCR results of CCs gene expression in the three patients’ groups: ZGP, ZBNP and ZGNP and ANOVA/ F-test as described before.

ABBREVIATIONS

ZP: 

Zona pellucida

ART: 

assisted reproductive technologies

IVF: 

in vitro fertilization

CCs: 

cumulus cells

COC: 

cumulus-oocyte complex

MII: 

metaphase II

eSET: 

elective single embryo transfer

GCs: 

granulosa cells

SSH: 

suppressive subtractive hybridization

ICSI: 

intracytoplasmic sperm injection

ECM: 

extracellular matrix

GWAS: 

genome-wide association studies

CRCs: 

corona radiata cells

IVT: 

in vitro transcription

ANOVA: 

analysis of variance

Declarations

ACKNOWLEDGMENT

This study was supported by The Natural Sciences and Engineering Research Council of Canada (NSERC), Canada ; and King AbdulAziz City for Science and Technology (KACST Large grant no. ARP-34-210).

DECLARATIONS

The publication charges for this article have been funded by the Center of Excellence in Genomic Medicine Research, King Abdulaziz University, Jeddah, Saudi Arabia.

This article has been published as part of BMC Genomics Volume 16 Supplement 1, 2015: Selected articles from the 2nd International Genomic Medical Conference (IGMC 2013): Genomics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/16/S1

Authors’ Affiliations

(1)
Center of Excellence in Genomic Medicine Research, King AbdulAziz University
(2)
KACST Technology Innovation Center in Personalized Medicine, King AbdulAziz University
(3)
Centre de Recherche en Biologie de la Reproduction, Laval University
(4)
Department of Gynecological Endocrinology and Reproductive Medicine, Bonn University
(5)
ilabcomm Gm bH

References

  1. Zegers-Hochschild F, Adamson GD, de Mouzon J, Ishihara O, Mansour R, Nygren K, Sullivan E, Vanderpoel S: International Committee for Monitoring Assisted Reproductive Technology (ICMART) and the World Health Organization (WHO) revised glossary of ART terminology, 2009. Fertility and sterility. 2009, 92 (5): 1520-1524. 10.1016/j.fertnstert.2009.09.009.PubMedGoogle Scholar
  2. Krisher RL: The effect of oocyte quality on development. J Anim Sci. 2004, 82 (E-Suppl): E14-23.PubMedGoogle Scholar
  3. Feuerstein P, Cadoret V, Dalbies-Tran R, Guerif F, Bidault R, Royere D: Gene expression in human cumulus cells: one approach to oocyte competence. Hum Reprod. 2007, 22 (12): 3069-3077. 10.1093/humrep/dem336.PubMedGoogle Scholar
  4. Sirard M-A, Assidi M: Screening for Oocyte Competence. Oocyte Physiology and Development in Domestic Animals. Edited by: Krisher RL. 2013, Wiley-Blackwell, 191-206.Google Scholar
  5. Sirard MA, Richard F, Blondin P, Robert C: Contribution of the oocyte to embryo quality. Theriogenology. 2006, 65 (1): 126-136. 10.1016/j.theriogenology.2005.09.020.PubMedGoogle Scholar
  6. Blondin P, Sirard MA: Oocyte and follicular morphology as determining characteristics for developmental competence in bovine oocytes. Mol Reprod Dev. 1995, 41 (1): 54-62. 10.1002/mrd.1080410109.PubMedGoogle Scholar
  7. Ebner T, Moser M, Sommergruber M, Tews G: Selection based on morphological assessment of oocytes and embryos at different stages of preimplantation development: a review. Human reproduction update. 2003, 9 (3): 251-262. 10.1093/humupd/dmg021.PubMedGoogle Scholar
  8. Yu Y, Mai Q, Chen X, Wang L, Gao L, Zhou C, Zhou Q: Assessment of the developmental competence of human somatic cell nuclear transfer embryos by oocyte morphology classification. Hum Reprod. 2009, 24 (3): 649-657.PubMedGoogle Scholar
  9. Chan PJ: Developmental potential of human oocytes according to zona pellucida thickness. J In Vitro Fert Embryo Transf. 1987, 4 (4): 237-241. 10.1007/BF01533763.PubMedGoogle Scholar
  10. Borini A, Lagalla C, Cattoli M, Sereni E, Sciajno R, Flamigni C, Coticchio G: Predictive factors for embryo implantation potential. Reproductive biomedicine online. 2005, 10 (5): 653-668. 10.1016/S1472-6483(10)61675-6.PubMedGoogle Scholar
  11. Tanghe S, Van Soom A, Nauwynck H, Coryn M, de Kruif A: Minireview: Functions of the cumulus oophorus during oocyte maturation, ovulation, and fertilization. Mol Reprod Dev. 2002, 61 (3): 414-424. 10.1002/mrd.10102.PubMedGoogle Scholar
  12. Coticchio G, Sereni E, Serrao L, Mazzone S, Iadarola I, Borini A: What criteria for the definition of oocyte quality?. Annals of the New York Academy of Sciences. 2004, 1034: 132-144. 10.1196/annals.1335.016.PubMedGoogle Scholar
  13. Mehri S, Levi Setti PE, Greco K, Sakkas D, Martinez G, Patrizio P: Correlation between follicular diameters and flushing versus no flushing on oocyte maturity, fertilization rate and embryo quality. Journal of assisted reproduction and genetics. 2014, 31 (1): 73-77. 10.1007/s10815-013-0124-9.PubMed CentralPubMedGoogle Scholar
  14. Thompson JG, Lane M, Gilchrist RB: Metabolism of the bovine cumulus-oocyte complex and influence on subsequent developmental competence. Society of Reproduction and Fertility supplement. 2007, 64: 179-190.PubMedGoogle Scholar
  15. Sugiura K, Pendola FL, Eppig JJ: Oocyte control of metabolic cooperativity between oocytes and companion granulosa cells: energy metabolism. Dev Biol. 2005, 279 (1): 20-30. 10.1016/j.ydbio.2004.11.027.PubMedGoogle Scholar
  16. Sutton-McDowall M, Gilchrist R, Thompson J: The pivotal role of glucose metabolism in determining oocyte developmental competence. Reproduction. 2010Google Scholar
  17. Assidi M, Richard FJ, Sirard MA: FSH in vitro versus LH in vivo: similar genomic effects on the cumulus. Journal of ovarian research. 2013, 6 (1): 68-10.1186/1757-2215-6-68.PubMed CentralPubMedGoogle Scholar
  18. Wassarman PM: Zona pellucida glycoproteins. Annu Rev Biochem. 1988, 57: 415-442. 10.1146/annurev.bi.57.070188.002215.PubMedGoogle Scholar
  19. Pokkyla RM, Lakkakorpi JT, Nuojua-Huttunen SH, Tapanainen JS: Sequence variations in human ZP genes as potential modifiers of zona pellucida architecture. Fertility and sterility. 2011, 95 (8): 2669-2672. 10.1016/j.fertnstert.2011.01.168.PubMedGoogle Scholar
  20. Lefièvre L, Conner SJ, Salpekar A, Olufowobi O, Ashton P, Pavlovic B, Lenton W, Afnan M, Brewis IA, Monk M, et al: Four zona pellucida glycoproteins are expressed in the human. Hum Reprod. 2004, 19 (7): 1580-1586. 10.1093/humrep/deh301.PubMedGoogle Scholar
  21. Hughes DC, Barratt CL: Identification of the true human orthologue of the mouse Zp1 gene: evidence for greater complexity in the mammalian zona pellucida?. Biochimica et biophysica acta. 1999, 1447 (2-3): 303-306. 10.1016/S0167-4781(99)00181-5.PubMedGoogle Scholar
  22. Chirinos M, Carino C, Gonzalez-Gonzalez ME, Arreola E, Reveles R, Larrea F: Characterization of human sperm binding to homologous recombinant zona pellucida proteins. Reprod Sci. 2011, 18 (9): 876-885. 10.1177/1933719111398146.PubMedGoogle Scholar
  23. Dunbar BS, Avery S, Lee V, Prasad S, Schwahn D, Schwoebel E, Skinner S, Wilkins B: The mammalian zona pellucida: its biochemistry, immunochemistry, molecular biology, and developmental expression. Reproduction, fertility, and development. 1994, 6 (3): 331-347. 10.1071/RD9940331.PubMedGoogle Scholar
  24. Lunn MO, Wright SJ: Analysis of the Ultrastructure of the Canine Zona Pellucida. Microscopy and microanalysis : the official journal of Microscopy Society of America, Microbeam Analysis Society, Microscopical Society of Canada. 2006, 12 (2): 270-271.Google Scholar
  25. Gook DA, Edgar DH, Borg J, Martic M: Detection of zona pellucida proteins during human folliculogenesis. Hum Reprod. 2008, 23 (2): 394-402.PubMedGoogle Scholar
  26. Green DP: Three-dimensional structure of the zona pellucida. Rev Reprod. 1997, 2 (3): 147-156. 10.1530/ror.0.0020147.PubMedGoogle Scholar
  27. Wassarman PM, Litscher ES: Mammalian fertilization: the egg's multifunctional zona pellucida. Int J Dev Biol. 2008, 52 (5-6): 665-676. 10.1387/ijdb.072524pw.PubMedGoogle Scholar
  28. Jimenez-Movilla M, Aviles M, Gomez-Torres MJ, Fernandez-Colom PJ, Castells MT, de Juan J, Romeu A, Ballesta J: Carbohydrate analysis of the zona pellucida and cortical granules of human oocytes by means of ultrastructural cytochemistry. Hum Reprod. 2004, 19 (8): 1842-1855. 10.1093/humrep/deh311.PubMedGoogle Scholar
  29. Familiari G, Heyn R, Relucenti M, Nottola SA, Sathananthan AH: Ultrastructural dynamics of human reproduction, from ovulation to fertilization and early embryo development. Int Rev Cytol. 2006, 249: 53-141.PubMedGoogle Scholar
  30. Ni Y, Li K, Xu W, Song L, Yao K, Zhang X, Huang H, Zhang Y, Shi QX: Acrosome reaction induced by recombinant human zona pellucida 3 peptides rhuZP3a22 approximately 176 and rhuZP3b177 approximately 348 and their mechanism. Journal of andrology. 2007, 28 (3): 381-388.PubMedGoogle Scholar
  31. Brewis IA, Clayton R, Barratt CL, Hornby DP, Moore HD: Recombinant human zona pellucida glycoprotein 3 induces calcium influx and acrosome reaction in human spermatozoa. Molecular human reproduction. 1996, 2 (8): 583-589. 10.1093/molehr/2.8.583.PubMedGoogle Scholar
  32. Montag M, Schimming T, Koster M, Zhou C, Dorn C, Rosing B, van der Ven H, Ven der Ven K: Oocyte zona birefringence intensity is associated with embryonic implantation potential in ICSI cycles. Reproductive biomedicine online. 2008, 16 (2): 239-244. 10.1016/S1472-6483(10)60580-9.PubMedGoogle Scholar
  33. Shi W, Xu B, Wu LM, Jin RT, Luan HB, Luo LH, Zhu Q, Johansson L, Liu YS, Tong XH: Oocytes with a Dark Zona Pellucida Demonstrate Lower Fertilization, Implantation and Clinical Pregnancy Rates in IVF/ICSI Cycles. PloS one. 2014, 9 (2): e89409-10.1371/journal.pone.0089409.PubMed CentralPubMedGoogle Scholar
  34. Omidi M, Khalili MA, Ashourzadeh S, Rahimipour M: Zona pellucida birefringence and meiotic spindle visualisation of human oocytes are not influenced by IVM technology. Reproduction, fertility, and development. 2014, 26 (3): 407-413. 10.1071/RD13001.PubMedGoogle Scholar
  35. Madaschi C, Aoki T, de Almeida Ferreira Braga DP, de Cassia Savio Figueira R, Semiao Francisco L, Iaconelli A, Borges E: Zona pellucida birefringence score and meiotic spindle visualization in relation to embryo development and ICSI outcomes. Reproductive biomedicine online. 2009, 18 (5): 681-686. 10.1016/S1472-6483(10)60014-4.PubMedGoogle Scholar
  36. Ebner T, Balaban B, Moser M, Shebl O, Urman B, Ata B, Tews G: Automatic user-independent zona pellucida imaging at the oocyte stage allows for the prediction of preimplantation development. Fertility and sterility. 2009Google Scholar
  37. Labrecque R, Sirard MA: The study of mammalian oocyte competence by transcriptome analysis: progress and challenges. Molecular human reproduction. 2014, 20 (2): 103-116. 10.1093/molehr/gat082.PubMedGoogle Scholar
  38. Guerif F, Le Gouge A, Giraudeau B, Poindron J, Bidault R, Gasnier O, Royere D: Limited value of morphological assessment at days 1 and 2 to predict blastocyst development potential: a prospective study based on 4042 embryos. Hum Reprod. 2007, 22 (7): 1973-1981. 10.1093/humrep/dem100.PubMedGoogle Scholar
  39. Templeton A: Joseph Price oration. The multiple gestation epidemic: the role of the assisted reproductive technologies. American journal of obstetrics and gynecology. 2004, 190 (4): 894-898. 10.1016/j.ajog.2004.02.017.PubMedGoogle Scholar
  40. Adashi EY, Barri PN, Berkowitz R, Braude P, Bryan E, Carr J, Cohen J, Collins J, Devroey P, Frydman R, et al: Infertility therapy-associated multiple pregnancies (births): an ongoing epidemic. Reproductive biomedicine online. 2003, 7 (5): 515-542. 10.1016/S1472-6483(10)62069-X.PubMedGoogle Scholar
  41. Bromer JG, Seli E: Assessment of embryo viability in assisted reproductive technology: shortcomings of current approaches and the emerging role of metabolomics. Current opinion in obstetrics & gynecology. 2008, 20 (3): 234-241. 10.1097/GCO.0b013e3282fe723d.Google Scholar
  42. Okun N, Sierra S, Wilson RD, Audibert F, Brock JA, Campagnolo C, Carroll J, Cartier L, Chitayat D, Gagnon A, et al: Pregnancy outcomes after assisted human reproduction. Journal of obstetrics and gynaecology Canada : JOGC = Journal d'obstetrique et gynecologie du Canada : JOGC. 2014, 36 (1): 64-83.PubMedGoogle Scholar
  43. Gerris J: Single-embryo transfer versus multiple-embryo transfer. Reproductive biomedicine online. 2009, 18 (Suppl 2): 63-70.PubMedGoogle Scholar
  44. Forman EJ, Tao X, Ferry KM, Taylor D, Treff NR, Scott RT: Single embryo transfer with comprehensive chromosome screening results in improved ongoing pregnancy rates and decreased miscarriage rates. Hum Reprod. 2012, 27 (4): 1217-1222. 10.1093/humrep/des020.PubMed CentralPubMedGoogle Scholar
  45. Straughen JK, Salihu HM, Keith L, Petrozzino J, Jones C: Obligatory versus elective single embryo transfer in in vitro fertilization. A population-based analysis of data from the U.K. Human Fertilisation and Embryology Authority. The Journal of reproductive medicine. 2013, 58 (3-4): 95-100.PubMedGoogle Scholar
  46. Hamel M, Dufort I, Robert C, Leveille MC, Leader A, Sirard MA: Identification of follicular marker genes as pregnancy predictors for human IVF: new evidence for the involvement of luteinization process. Molecular human reproduction. 2010a, 16 (8): 548-556. 10.1093/molehr/gaq051.Google Scholar
  47. Hamel M, Dufort I, Robert C, Leveille MC, Leader A, Sirard MA: Genomic assessment of follicular marker genes as pregnancy predictors for human IVF. Molecular human reproduction. 2010, 16 (2): 87-96. 10.1093/molehr/gap079.PubMedGoogle Scholar
  48. Assidi M, Sirard M-A: Cumulus Cell Gene Expression as a Marker of Oocyte Quality. Oogenesis. Edited by: Coticchio G, Albertini DF, De Santis L. 2013, London: Springer, 231-252.Google Scholar
  49. Assidi M, Montag M, Van der Ven K, Sirard MA: Biomarkers of human oocyte developmental competence expressed in cumulus cells before ICSI: a preliminary study. Journal of assisted reproduction and genetics. 2011, 28 (2): 173-188. 10.1007/s10815-010-9491-7.PubMed CentralPubMedGoogle Scholar
  50. Seli E, Robert C, Sirard MA: OMICS in assisted reproduction: possibilities and pitfalls. Molecular human reproduction. 2010Google Scholar
  51. Hamel M, Dufort I, Robert C, Gravel C, Léveillé M-C, Leader A, Sirard M-A: Identification of differentially expressed markers in human follicular cells associated with competent oocytes. Human Reproduction. 2008Google Scholar
  52. Tesfaye D, Hossain MM, Schellander K: MicroRNAs in Oocyte Physiology and Development. Oocyte Physiology and Development in Domestic Animals. 2013, Wiley-Blackwell, 71-83.Google Scholar
  53. Geraedts J, Montag M, Magli MC, Repping S, Handyside A, Staessen C, Harper J, Schmutzler A, Collins J, Goossens V, et al: Polar body array CGH for prediction of the status of the corresponding oocyte. Part I: clinical results. Hum Reprod. 2011, 26 (11): 3173-3180. 10.1093/humrep/der294.PubMed CentralPubMedGoogle Scholar
  54. Ursula E-R: Essential Features of the Mature Oocyte. Preservation of Human Oocytes. 2009, CRC Press, 95-103.Google Scholar
  55. Lasiene K, Vitkus A, Valanciute A, Lasys V: Morphological criteria of oocyte quality. Medicina (Kaunas). 2009, 45 (7): 509-515.Google Scholar
  56. Assidi M, Dufort I, Ali A, Hamel M, Algriany O, Dielemann S, Sirard MA: Identification of Potential Markers of Oocyte Competence Expressed in Bovine Cumulus Cells Matured with Follicle-Stimulating Hormone and/or Phorbol Myristate Acetate In Vitro. Biology of reproduction. 2008, 79 (2): 209-222. 10.1095/biolreprod.108.067686.PubMedGoogle Scholar
  57. Hamel M, Dufort I, Robert C, Gravel C, Leveille MC, Leader A, Sirard MA: Identification of differentially expressed markers in human follicular cells associated with competent oocytes. Hum Reprod. 2008, 23 (5): 1118-1127. 10.1093/humrep/den048.PubMedGoogle Scholar
  58. IPA: Ingenuity Pathways Analysis. Ingenuity® Systems, wwwingenuitycom. 2014Google Scholar
  59. Pelletier C, Keefe DL, Trimarchi JR: Noninvasive polarized light microscopy quantitatively distinguishes the multilaminar structure of the zona pellucida of living human eggs and embryos. Fertility and sterility. 2004, 81 (Suppl 1): 850-856.PubMedGoogle Scholar
  60. Rama Raju GA, Prakash GJ, Krishna KM, Madan K: Meiotic spindle and zona pellucida characteristics as predictors of embryonic development: a preliminary study using PolScope imaging. Reproductive biomedicine online. 2007, 14 (2): 166-174. 10.1016/S1472-6483(10)60784-5.PubMedGoogle Scholar
  61. Allworth AE, Albertini DF: Meiotic maturation in cultured bovine oocytes is accompanied by remodeling of the cumulus cell cytoskeleton. Dev Biol. 1993, 158 (1): 101-112. 10.1006/dbio.1993.1171.PubMedGoogle Scholar
  62. Sutovsky P, Flechon JE, Flechon B, Motlik J, Peynot N, Chesne P, Heyman Y: Dynamic changes of gap junctions and cytoskeleton during in vitro culture of cattle oocyte cumulus complexes. Biology of reproduction. 1993, 49 (6): 1277-1287. 10.1095/biolreprod49.6.1277.PubMedGoogle Scholar
  63. Suzuki H, Jeong BS, Yang X: Dynamic changes of cumulus-oocyte cell communication during in vitro maturation of porcine oocytes. Biology of reproduction. 2000, 63 (3): 723-729. 10.1095/biolreprod63.3.723.PubMedGoogle Scholar
  64. Sato E, Yokoo M: Morphological and biochemical dynamics of porcine cumulus-oocyte complexes: role of cumulus expansion in oocyte maturation. Ital J Anat Embryol. 2005, 110 (2 Suppl 1): 205-217.PubMedGoogle Scholar
  65. Motta PM, Nottola SA, Pereda J, Croxatto HB, Familiari G: Ultrastructure of human cumulus oophorus: a transmission electron microscopic study on oviductal oocytes and fertilized eggs. Hum Reprod. 1995, 10 (9): 2361-2367. 10.1093/oxfordjournals.humrep.a136299.PubMedGoogle Scholar
  66. Russell DL, Salustri A: Extracellular matrix of the cumulus-oocyte complex. Semin Reprod Med. 2006, 24 (4): 217-227. 10.1055/s-2006-948551.PubMedGoogle Scholar
  67. Gilchrist RB, Lane M, Thompson JG: Oocyte-secreted factors: regulators of cumulus cell function and oocyte quality. Human reproduction update. 2008, 14 (2): 159-177. 10.1093/humupd/dmm040.PubMedGoogle Scholar
  68. Van Soom A, Tanghe S, De Pauw I, Maes D, de Kruif A: Function of the cumulus oophorus before and during mammalian fertilization. Reproduction in domestic animals = Zuchthygiene. 2002, 37 (3): 144-151. 10.1046/j.1439-0531.2002.00345.x.PubMedGoogle Scholar
  69. Shimada M, Richards JS: Cumulus Cells are an Essential Mediator of Ovulation Stimuli from Granulosa Cells to Oocyte. Journal of Mammalian Ova Research. 2010, 27 (1): 2-10. 10.1274/jmor.27.2.Google Scholar
  70. Familiari G, Relucenti M, Heyn R, Micara G, Correr S: Three-dimensional structure of the zona pellucida at ovulation. Microsc Res Tech. 2006, 69 (6): 415-426. 10.1002/jemt.20301.PubMedGoogle Scholar
  71. Choi JH, Choi KC, Auersperg N, Leung PC: Gonadotropins upregulate the epidermal growth factor receptor through activation of mitogen-activated protein kinases and phosphatidyl-inositol-3-kinase in human ovarian surface epithelial cells. Endocrine-related cancer. 2005, 12 (2): 407-421. 10.1677/erc.1.00896.PubMedGoogle Scholar
  72. Freimann S, Ben-Ami I, Dantes A, Ron-El R, Amsterdam A: EGF-like factor epiregulin and amphiregulin expression is regulated by gonadotropins/cAMP in human ovarian follicular cells. Biochemical and biophysical research communications. 2004, 324 (2): 829-834. 10.1016/j.bbrc.2004.09.129.PubMedGoogle Scholar
  73. Dickey RP, Brinsden PR: Ovulation induction for intrauterine insemination: gonadotropins and oral drug-gonadotropin combinations. Manual of Intrauterine Insemination and Ovulation Induction. Edited by: Dickey RP, Brinsden PR, Pyrzak R. 2010, Cambridge University Press, 80-92.Google Scholar
  74. Pangas SA, Matzuk MM: The art and artifact of GDF9 activity: cumulus expansion and the cumulus expansion-enabling factor. Biology of reproduction. 2005, 73 (4): 582-585. 10.1095/biolreprod.105.042127.PubMedGoogle Scholar
  75. Su YQ, Wu X, O'Brien MJ, Pendola FL, Denegre JN, Matzuk MM, Eppig JJ: Synergistic roles of BMP15 and GDF9 in the development and function of the oocyte-cumulus cell complex in mice: genetic evidence for an oocyte-granulosa cell regulatory loop. Dev Biol. 2004, 276 (1): 64-73. 10.1016/j.ydbio.2004.08.020.PubMedGoogle Scholar
  76. Su YQ, Sugiura K, Wigglesworth K, O'Brien MJ, Affourtit JP, Pangas SA, Matzuk MM, Eppig JJ: Oocyte regulation of metabolic cooperativity between mouse cumulus cells and oocytes: BMP15 and GDF9 control cholesterol biosynthesis in cumulus cells. Development. 2008, 135 (1): 111-121.PubMedGoogle Scholar
  77. Gregory SJ, Kaiser UB: Regulation of gonadotropins by inhibin and activin. Semin Reprod Med. 2004, 22 (3): 253-267. 10.1055/s-2004-831901.PubMedGoogle Scholar
  78. Yamashita Y, Shimada M: The release of EGF domain from EGF-like factors by a specific cleavage enzyme activates the EGFR-MAPK3/1 pathway in both granulosa cells and cumulus cells during the ovulation process. J Reprod Dev. 2012, 58 (5): 510-514. 10.1262/jrd.2012-056.PubMedGoogle Scholar
  79. Shimada M, Hernandez-Gonzalez I, Gonzalez-Robanya I, Richards JS: Induced expression of pattern recognition receptors in cumulus oocyte complexes: novel evidence for innate immune-like functions during ovulation. Mol Endocrinol. 2006, 20 (12): 3228-3239. 10.1210/me.2006-0194.PubMedGoogle Scholar
  80. Li M, Liang CG, Xiong B, Xu BZ, Lin SL, Hou Y, Chen DY, Schatten H, Sun QY: PI3-kinase and mitogen-activated protein kinase in cumulus cells mediate EGF-induced meiotic resumption of porcine oocyte. Domestic animal endocrinology. 2008a, 34 (4): 360-371. 10.1016/j.domaniend.2007.10.001.Google Scholar
  81. Chen X, Zhou B, Yan J, Xu B, Tai P, Li J, Peng S, Zhang M, Xia G: Epidermal growth factor receptor activation by protein kinase C is necessary for FSH-induced meiotic resumption in porcine cumulus-oocyte complexes. J Endocrinol. 2008, 197 (2): 409-419. 10.1677/JOE-07-0592.PubMedGoogle Scholar
  82. Sriraman V, Rudd MD, Lohmann SM, Mulders SM, Richards JS: Cyclic guanosine 5'-monophosphate-dependent protein kinase II is induced by luteinizing hormone and progesterone receptor-dependent mechanisms in granulosa cells and cumulus oocyte complexes of ovulating follicles. Mol Endocrinol. 2006, 20 (2): 348-361. 10.1210/me.2005-0317.PubMedGoogle Scholar
  83. Sirard MA, Desrosier S, Assidi M: In vivo and in vitro effects of FSH on oocyte maturation and developmental competence. Theriogenology. 2007, 68 (Suppl 1): S71-76.PubMedGoogle Scholar
  84. Espey LL: Current status of the hypothesis that mammalian ovulation is comparable to an inflammatory reaction. Biology of reproduction. 1994, 50 (2): 233-238. 10.1095/biolreprod50.2.233.PubMedGoogle Scholar
  85. Richards JS: Genetics of ovulation. Semin Reprod Med. 2007, 25 (4): 235-242. 10.1055/s-2007-980217.PubMedGoogle Scholar
  86. Richards JS: Ovulation: new factors that prepare the oocyte for fertilization. Molecular and cellular endocrinology. 2005, 234 (1-2): 75-79. 10.1016/j.mce.2005.01.004.PubMedGoogle Scholar
  87. Hernandez-Gonzalez I, Gonzalez-Robayna I, Shimada M, Wayne CM, Ochsner SA, White L, Richards JS: Gene expression profiles of cumulus cell oocyte complexes during ovulation reveal cumulus cells express neuronal and immune-related genes: does this expand their role in the ovulation process?. Mol Endocrinol. 2006, 20 (6): 1300-1321. 10.1210/me.2005-0420.PubMedGoogle Scholar
  88. Russell DL, Robker RL: Molecular mechanisms of ovulation: co-ordination through the cumulus complex. Human reproduction update. 2007, 13 (3): 289-312. 10.1093/humupd/dml062.PubMedGoogle Scholar
  89. McKenzie LJ, Pangas SA, Carson SA, Kovanci E, Cisneros P, Buster JE, Amato P, Matzuk MM: Human cumulus granulosa cell gene expression: a predictor of fertilization and embryo selection in women undergoing IVF. Hum Reprod. 2004, 19 (12): 2869-2874. 10.1093/humrep/deh535.PubMedGoogle Scholar
  90. Zhang X, Jafari N, Barnes RB, Confino E, Milad M, Kazer RR: Studies of gene expression in human cumulus cells indicate pentraxin 3 as a possible marker for oocyte quality. Fertility and sterility. 2005, 83 (Suppl 1): 1169-1179.PubMedGoogle Scholar
  91. Shimada M, Hernandez-Gonzalez I, Gonzalez-Robanya I, Richards JS: Induced expression of pattern recognition receptors in cumulus oocyte complexes: novel evidence for innate immune-like functions during ovulation. Mol Endocrinol. 2006a, 20 (12): 3228-3239. 10.1210/me.2006-0194.Google Scholar
  92. Espey LL, Yoshioka S, Russell DL, Robker RL, Fujii S, Richards JS: Ovarian expression of a disintegrin and metalloproteinase with thrombospondin motifs during ovulation in the gonadotropin-primed immature rat. Biology of reproduction. 2000, 62 (4): 1090-1095. 10.1095/biolreprod62.4.1090.PubMedGoogle Scholar
  93. Lourenco B, Sousa AP, Almeida-Santos T, Ramalho-Santos J: Relation of cumulus cell status with single oocyte maturity, fertilization capability and patient age. Journal of reproduction & infertility. 2014, 15 (1): 15-21.Google Scholar
  94. Host E, Gabrielsen A, Lindenberg S, Smidt-Jensen S: Apoptosis in human cumulus cells in relation to zona pellucida thickness variation, maturation stage, and cleavage of the corresponding oocyte after intracytoplasmic sperm injection. Fertility and sterility. 2002, 77 (3): 511-515. 10.1016/S0015-0282(01)03006-0.PubMedGoogle Scholar
  95. Stouffer RL, Xu F, Duffy DM: Molecular control of ovulation and luteinization in the primate follicle. Front Biosci. 2007, 12: 297-307. 10.2741/2065.PubMedGoogle Scholar
  96. van Montfoort AP, Geraedts JP, Dumoulin JC, Stassen AP, Evers JL, Ayoubi TA: Differential gene expression in cumulus cells as a prognostic indicator of embryo viability: a microarray analysis. Molecular human reproduction. 2008, 14 (3): 157-168. 10.1093/molehr/gam088.PubMedGoogle Scholar
  97. Lee KS, Joo BS, Na YJ, Yoon MS, Choi OH, Kim WW: Cumulus cells apoptosis as an indicator to predict the quality of oocytes and the outcome of IVF-ET. Journal of assisted reproduction and genetics. 2001, 18 (9): 490-498. 10.1023/A:1016649026353.PubMed CentralPubMedGoogle Scholar
  98. Assou S, Haouzi D, De Vos J, Hamamah S: Human cumulus cells as biomarkers for embryo and pregnancy outcomes. Molecular human reproduction. 2010Google Scholar
  99. Akison LK, Robker RL: The critical roles of progesterone receptor (PGR) in ovulation, oocyte developmental competence and oviductal transport in mammalian reproduction. Reproduction in domestic animals = Zuchthygiene. 2012, 47 (Suppl 4): 288-296.PubMedGoogle Scholar
  100. Zafeiriou S, Loutradis D, Michalas S: The role of gonadotropins in follicular development and their use in ovulation induction protocols for assisted reproduction. Eur J Contracept Reprod Health Care. 2000, 5 (2): 157-167. 10.1080/13625180008500389.PubMedGoogle Scholar
  101. Hamel M, Dufort I, Robert C, Leveille MC, Leader A, Sirard MA: Identification of follicular marker genes as pregnancy predictors for human IVF: new evidence for the involvement of luteinization process. Molecular human reproduction. 2010, 16 (8): 548-556. 10.1093/molehr/gaq051.PubMedGoogle Scholar
  102. Chian RC, Ao A, Clarke HJ, Tulandi T, Tan SL: Production of steroids from human cumulus cells treated with different concentrations of gonadotropins during culture in vitro. Fertility and sterility. 1999, 71 (1): 61-66. 10.1016/S0015-0282(98)00416-6.PubMedGoogle Scholar
  103. Grosse J, Bulling A, Brucker C, Berg U, Amsterdam A, Mayerhofer A, Gratzl M: Synaptosome-associated protein of 25 kilodaltons in oocytes and steroid-producing cells of rat and human ovary: molecular analysis and regulation by gonadotropins. Biology of reproduction. 2000, 63 (2): 643-650. 10.1095/biolreprod63.2.643.PubMedGoogle Scholar
  104. Sirard MA, Desrosier S, Assidi M: In vivo and in vitro effects of FSH on oocyte maturation and developmental competence. Theriogenology. 2007Google Scholar
  105. Uyar A, Torrealday S, Seli E: Cumulus and granulosa cell markers of oocyte and embryo quality. Fertility and sterility. 2013, 99 (4): 979-997. 10.1016/j.fertnstert.2013.01.129.PubMedGoogle Scholar
  106. Van Der Ven K, Montag M, Drengner C, Köster M, Van Der Ven H: Differential gene expression profiles in cells of the corona radiata and outer cumulus oophorus in relation to oocyte competence. Abstracts of the 25th Annual Meeting of ESHRE, Amsterdam, The Netherlands, 28 June - 1 July, 2009: 28 June - 1 July, 2009 2009; Amsterdam, The Netherlands. 2009, i32-Google Scholar
  107. Kaser DJ, Racowsky C: Clinical outcomes following selection of human preimplantation embryos with time-lapse monitoring: a systematic review. Human reproduction update. 2014Google Scholar
  108. Polanski LT, Coelho Neto MA, Nastri CO, Navarro PA, Ferriani RA, Raine-Fenning N, Martins WP: Time-lapse embryo imaging for improving reproductive outcomes: a systematic review and meta-analysis. Ultrasound in obstetrics & gynecology : the official journal of the International Society of Ultrasound in Obstetrics and Gynecology. 2014Google Scholar
  109. Organization WH: WHO laboratory manual for the examination and processing of human semen. 2010, Geneva, Switzerland: WHOGoogle Scholar
  110. NIA: National Institute on Aging (NIA) Array Analysis Tool http://lgsun.grc.nia.nih.gov/ANOVA/. 2013, National Institutes of Health(Laboratory of Genetics, Baltimore MD , USA)Google Scholar
  111. GeNorm: http://medgen.ugent.be/~jvdesomp/genorm/. Biogazelle's qbase+ software. 2013,Google Scholar
  112. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome biology. 2002, 3 (7): RESEARCH0034-PubMed CentralPubMedGoogle Scholar

Copyright

© ASSIDI et al; licensee BioMed Central Ltd. 2015

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Advertisement