Providing a stable methodological basis for comparing transcript abundance of developing embryos using microarrays

Comprehensive transcriptomic analysis of long non-coding RNAs in bovine ovarian follicles and early embryos

Long non-coding RNAs (lncRNAs) have been the subject of numerous studies over the past decade. First thought to come from aberrant transcriptional events, lncRNAs are now considered a crucial component of the genome with roles in multiple cellular functions. However, the functional annotation and characterization of bovine lncRNAs during early development remain limited. In this comprehensive analysis, we review lncRNAs expression in bovine ovarian follicles and early embryos, based on a unique database comprising 468 microarray hybridizations from a single platform designed to target 7,724 lncRNA transcripts, of which 5,272 are intergenic (lincRNA), 958 are intronic, and 1,524 are antisense (lncNAT). Compared to translated mRNA, lncRNAs have been shown to be more tissue-specific and expressed in low copy numbers. This analysis revealed that protein-coding genes and lncRNAs are both expressed more in oocytes. Differences between the oocyte and the 2-cell embryo are also more apparent in terms of lncRNAs than mRNAs. Co-expression network analysis using WGCNA generated 25 modules with differing proportions of lncRNAs. The modules exhibiting a higher proportion of lncRNAs were found to be associated with fewer annotated mRNAs and housekeeping functions. Functional annotation of co-expressed mRNAs allowed attribution of lncRNAs to a wide array of key cellular events such as meiosis, translation initiation, immune response, and mitochondrial related functions. We thus provide evidence that lncRNAs play diverse physiological roles that are tissue-specific and associated with key cellular functions alongside mRNAs in bovine ovarian follicles and early embryos. This contributes to add lncRNAs as active molecules in the complex regulatory networks driving folliculogenesis, oogenesis and early embryogenesis all of which are necessary for reproductive success.

Folliculogenesis and acquisition of oocyte competence in cows

IVF success depends on hundreds of factors and details but the oocyte quality remains the most important and problematic issue. All antral follicles contain oocytes and all of them have that have reached their full size, can be aspirated, can mature and can be fertilized in vitro. But only a few will make it to embryo unless harvested at a very specific time/status. The conditions impacting the oocyte competence are essentially dependant on the follicular status. Growing follicles contains oocytes that have not completed their preparation, as they are still writing information (RNA), later, dominant follicles or follicles at the plateau phase, stop transcription and become candidates for development. Once in transcriptional arrest, the oocytes, if not ovulated in a short amount of time, do not always make good embryos. This window is affected by time and follicle size and looks like a bell curve. The following review further explain the physiological and molecular evidences that we have to illustrate the competence window and provides clues on how to optimize ovarian stimulation to maximise oocyte quality.

Transcriptome Profiling of In-Vivo Produced Bovine Pre-implantation Embryos Using Two-color Microarray Platform

Early embryonic loss is a large contributor to infertility in cattle. Moreover, bovine becomes an interesting model to study human preimplantation embryo development due to their similar developmental process. Although genetic factors are known to affect early embryonic development, the discovery of such factors has been a serious challenge. Microarray technology allows quantitative measurement and gene expression profiling of transcript levels on a genome-wide basis. One of the main decisions that have to be made when planning a microarray experiment is whether to use a one- or two-color approach. Two-color design increases technical replication, minimizes variability, improves sensitivity and accuracy as well as allows having loop designs, defining the common reference samples. Although microarray is a powerful biological tool, there are potential pitfalls that can attenuate its power. Hence, in this technical paper we demonstrate an optimized protocol for RNA extraction, amplification, labeling, hybridization of the labeled amplified RNA to the array, array scanning and data analysis using the two-color analysis strategy.

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.

Insulin exposure during in vitro bovine oocyte maturation changes blastocyst gene expression and developmental potential

Metabolic imbalance impairs fertility, because changes in concentrations of metabolites and hormones in the blood and follicular fluid create an unfavourable environment for early embryonic development. Insulin is a key metabolic hormone known for its effects on fertility: insulin concentrations are increased during energy balance disturbances in diabetes or metabolic syndrome. Still, insulin is frequently used at supraphysiological concentrations for embryo in vitro culture with unknown consequences for the developmental potential of the offspring. In the present study we investigated the effects of insulin exposure during in vitro bovine oocyte maturation on developmental rates, embryo quality and gene expression. Supplementation of the maturation media with insulin at 10 or 0.1 µg mL-1 decreased blastocyst rates compared with an insulin-free control (19.8 ± 1.3% and 20.4 ± 1.3% vs 23.8 ± 1.3%, respectively; P < 0.05) and led to increased cell numbers (nearly 10% more cells on Day 8 compared with control; P < 0.05). Transcriptome analysis revealed significant upregulation of genes involved in lipid metabolism, nuclear factor (erythroid-derived 2)-like 2 (NRF2) stress response and cell differentiation, validated by quantitative polymerase chain reaction. To conclude, the results of the present study demonstrate that insulin exposure during in vitro oocyte maturation has a lasting effect on the embryo until the blastocyst stage, with a potential negative effect in the form of specific gene expression perturbations.

Breed-specific factors influence embryonic lipid composition: comparison between Jersey and Holstein

Some embryos exhibit better survival potential to cryopreservation than others. The cause of such a phenotype is still unclear and may be due to cell damage during cryopreservation, resulting from overaccumulation and composition of lipids. In cattle embryos, in vitro culture conditions have been shown to impact the number of lipid droplets within blastomeres. Thus far, the impact of breed on embryonic lipid content has not been studied. In the present study were compared the colour, lipid droplet abundance, lipid composition, mitochondrial activity and gene expression of in vivo-collected Jersey breed embryos, which are known to display poor performance post-freezing, with those of in vivo Holstein embryos, which have good cryotolerance. Even when housed and fed under the same conditions, Jersey embryos were found to be darker and contain more lipid droplets than Holstein embryos, and this was correlated with lower mitochondrial activity. Differential expression of genes associated with lipid metabolism and differences in lipid composition were found. These results show genetic background can impact embryonic lipid metabolism and storage.

Studying bovine early embryo transcriptome by microarray

Microarrays represent a significant advantage when studying gene expression in early embryo because they allow for a speedy study of a large number of genes even if the sample of interest contains small quantities of genetic material. Here we describe the protocols developed by the EmbryoGENE Network to study the bovine transcriptome in early embryo using a microarray experimental design.

Ribosomal RNA depletion or exclusion has negligible effect on the detection of viruses in a pan viral microarray

•

Ribosomal RNA depletion protocols were assessed to improve microarray performance.

•

The outcome was compared with random amplification protocol.

•

Ribosomal RNA depletion had little effect on the microarray performance.

Pan viral DNA microarrays, which can detect known, novel and multiple viral infections, are major laboratory assets contributing to the control of infectious diseases. The large quantity of ribosomal RNA (rRNA) found in tissue samples is thought to be a major factor contributing to the comparatively lower sensitivity of detecting RNA viruses, as a sequence-independent PCR is used to amplify unknown samples for microarray analysis. This study aimed to determine whether depletion or exclusion of rRNA can improve microarray detection and simplify its analysis. Therefore, two different rRNA depletion and exclusion protocols, RiboMinus™ technology and non-rRNA binding hexanucleotides, were applied to the microarray sample processing and the outcome was compared with those of the sequence-independent amplification protocol. This study concludes that the two procedures, described to deplete or exclude rRNA, have negligible effect on the microarrays detection and analysis and might only in combination with further techniques result in a significant enhancement of sensitivity. Currently, existing protocols of random amplification and background adjustment are pertinent for the purpose of sample processing for microarray analysis.

Toxicogenomic approaches in developmental toxicology testing

The emergence of toxicogenomic applications provides new tools to characterize, classify, and potentially predict teratogens. However, due to the vast number of experimental and statistical procedural steps, toxicogenomic studies are challenging. Here, we guide researchers through the basic framework of conducting toxicogenomic investigations in the field of developmental toxicology, providing examples of biological and technical factors that may influence response and interpretation. Furthermore, we review current, diverse applications of toxicogenomic-based approaches in teratology testing, including exposure-response characterization (dose and duration), chemical classification studies, and cross-model comparisons study designs. This review is intended to guide scientists through the challenging and complex structure of conducting toxicogenomic analyses, while considering the many applications of using toxicogenomics in study designs and the future of these types of "omics" approaches in developmental toxicology.

Genomic assessment of human cumulus cell marker genes as predictors of oocyte developmental competence: impact of various experimental factors

Background

Single embryo transfer (SET) is the most successful way to reduce the frequency of multiple pregnancies following in vitro fertilisation. However, selecting the embryo for SET with the highest chances of pregnancy remains a difficult challenge since morphological and kinetics criteria provide poor prediction of both developmental and implantation ability. Partly through the expression of specific genes, the oocyte-cumulus interaction helps the oocyte to acquire its developmental competence. Our aim was therefore to identify at the level of cumulus cells (CCs) genes related to oocyte developmental competence.

Methodology/Principal Findings

197 individual CCs were collected from 106 patients undergoing an intra-cytoplasmic sperm injection procedure. Gene expression of CCs was studied using microarray according to the nuclear maturity of the oocyte (immature vs. mature oocyte) and to the developmental competence of the oocyte (ability to reach the blastocyst stage after fertilisation). Microarray study was followed by a meta-analysis of the behaviour of these genes in other datasets available in Gene Expression Omnibus which showed the consistency of this list of genes. Finally, 8 genes were selected according to oocyte developmental competence from the 308 differentially expressed genes (p<0.0001) for further validation by quantitative PCR (qPCR). Three of these 8 selected genes were validated as potential biomarkers (PLIN2, RGS2 and ANG). Experimental factors such as inter-patient and qPCR series variability were then assessed using the Generalised Linear Mixed Model procedure, and only the expression level of RGS2 was confirmed to be related to oocyte developmental competence. The link between biomarkers and pregnancy was finally evaluated and level of RGS2 expression was also correlated with clinical pregnancy.

Conclusion/Significance

RGS2, known as a regulator of G protein signalling, was the only gene among our 8 selected candidates biomarkers of oocyte competence to cover many factors of variability, including inter-patient factors and experimental conditions.

A review of toxicogenomic approaches in developmental toxicology

Over the past decade, the use of gene expression profiling (i.e., toxicogenomics or transcriptomics) has been established as the vanguard "omics" technology to investigate exposure-induced molecular changes that underlie the development of disease. As this technology quickly advances, researchers are striving to keep pace in grasping the complexity of toxicogenomic response while at the same time determine its applicability for the field of developmental toxicology. Initial studies suggest toxicogenomics to be a promising tool for multiple types of study designs, including exposure-response investigations (dose and duration), chemical classification, and model comparisons. In this review, we examine the use of toxicogenomics in developmental toxicology, discussing biological and technical factors that influence response and interpretation. Additionally, we provide a framework to guide toxicogenomic investigations in the field of developmental toxicology.

Combining resources to obtain a comprehensive survey of the bovine embryo transcriptome through deep sequencing and microarrays

While most assisted reproductive technologies (ART) are considered routine for the reproduction of species of economical importance, such as the bovine, the impact of these manipulations on the developing embryo remains largely unknown. In an effort to obtain a comprehensive survey of the bovine embryo transcriptome and how it is modified by ART, resources were combined to design an embryo-specific microarray. Close to one million high-quality reads were produced from subtracted bovine embryo libraries using Roche 454 Titanium deep sequencing technology, which enabled the creation of an augmented bovine genome catalog. This catalog was enriched with bovine embryo transcripts, and included newly discovered indel type and 3'UTR variants. Using this augmented bovine genome catalog, the EmbryoGENE Bovine Microarray was designed and is composed of a total of 42,242 probes, including 21,139 known reference genes; 9,322 probes for novel transcribed regions (NTRs); 3,677 alternatively spliced exons; 3,353 3'-tiling probes; and 3,723 controls. A suite of bioinformatics tools was also developed to facilitate microrarray data analysis and database creation; it includes a quality control module, a Laboratory Information Management System (LIMS) and microarray analysis software. Results obtained during this study have already led to the identification of differentially expressed blastocyst targets, NTRs, splice variants of the indel type, and 3'UTR variants. We were able to confirm microarray results by real-time PCR, indicating that the EmbryoGENE bovine microarray has the power to detect physiologically relevant changes in gene expression.

Comprehensive cross production system assessment of the impact of in vitro microenvironment on the expression of messengers and long non-coding RNAs in the bovine blastocyst

In vitroproduction (IVP) of cattle embryos over the past two decades has revealed several negative impacts that have been attributed to the artificial microenvironment. Studies on embryos producedin vitroclearly point to aberrant gene expression levels. So far, the causal association between phenotype and measured gene expression has not led to substantial improvement of IVP systems. The aim of this study was to generate a unique dataset composed of microarray-derived relative transcript abundance values for blastocysts produced in tenin vitrosystems differing primarily in culture medium formulation. Between-group comparisons determine the level of overall similarity among systems relative toin vivoreference embryos. The use of the dataset to contrast allin vitrotreatments with thein vivoblastocysts pointed to a single common gene network. The ‘boutique’ array contained a panel of novel uncharacterized transcripts that were variably expressed depending on the medium in which the blastocysts were produced. These novel transcripts were differentially expressed in blastocysts even as carryover from conditions encountered 7 days earlier during oocyte maturation. All of the selected novel candidates thus expressed were from intergenic regions. The function of this long non-coding RNA remains unknown but clearly points to an additional level of complexity in early embryo development.

Microarray analysis of gene expression during early development: a cautionary overview

The rise of the 'omics' technologies started nearly a decade ago and, among them, transcriptomics has been used successfully to contrast gene expression in mammalian oocytes and early embryos. The scarcity of biological material that early developmental stages provide is the prime reason why the field of transcriptomics is becoming more and more popular with reproductive biologists. The potential to amplify scarce mRNA samples and generate the necessary amounts of starting material enables the relative measurement of RNA abundance of thousands of candidates simultaneously. So far, microarrays have been the most commonly used high-throughput method in this field. Microarray platforms can be found in a wide variety of formats, from cDNA collections to long or short oligo probe sets. These platforms generate large amounts of data that require the integration of comparative RNA abundance values in the physiological context of early development for their full benefit to be appreciated. Unfortunately, significant discrepancies between datasets suggest that direct comparison between studies is difficult and often not possible. We have investigated the sample-handling steps leading to the generation of microarray data produced from prehatching embryo samples and have identified key steps that significantly impact the downstream results. This review provides a discussion on the best methods for the preparation of samples from early embryos for microarray analysis and focuses on the challenges that impede dataset comparisons from different platforms and the reasons why methodological benchmarking performed using somatic cells may not apply to the atypical nature of prehatching development.