Application of multi-omics data integration and machine learning approaches to identify epigenetic and transcriptomic differences between in vitro and in vivo produced bovine embryos

The metabolomic composition of the oviductal fluid is controlled by the periovulatory hormonal context in Bos indicus cows†

In cattle, oviductal function is controlled by the ovarian sex-steroids estradiol and progesterone. Here, we tested the hypothesis that the exposure to contrasting sex-steroid milieus differentially impacts the oviductal fluid composition. Estrous cycles of non-lactating, multiparous Nelore cows were pre-synchronized and then synchronized with a protocol designed two induce ovulation of large or small follciles. Larger preovulatory follicle (day 0) and corpora lutea (day 4) and greater estradiol (day 0) and progesterone (day 4) concentrations were observed in the large follciles group. Four days after induced ovulation, oviductal fluid was collected post-mortem. Quantitative mass spectrometry was used to determine the concentration of amino acids, biogenic amines, acylcarnitines, phosphatidylcholines, lysophosphatidylcholines, sphingomyelins, hexoses, prostaglandins, and related compounds. Multivariate analyses (orthogonal projections to latent structures discriminant analysis) were conducted to compare the metabolomic signatures of oviductal fluids. Correlation network analysis was conducted to measure the strength and hierarchy of associations among metabolites. Of the 205 metabolites quantified, 171 were detected in at least 50% of the samples and were included in further data analysis. After orthogonal projections to latent structures discriminant analysis analysis, samples of the large follciles and small follciles were divided clearly into two non-overlapping clusters. Twenty metabolites had different or tended to have different concentrations in the oviductal fluid when comparing groups. Seven of these 20 analytes had greater concentrations in large follciles cows. Moreover, total sum of biogenic amines, phosphatidylcholines, and prostaglandins were higher in the small follciles group. The correlation network showed that the large follciles group metabolites' concentrations were highly intercorrelated, which was not observed in the small follciles group. We concluded that the periovulatory endocrine milieu regulates the composition of the oviductal fluid.

In vitro culture alters cell lineage composition and cellular metabolism of bovine blastocyst†

Profiling bovine blastocyst transcriptome at the single-cell level has enabled us to reveal the first cell lineage segregation, during which the inner cell mass (ICM), trophectoderm (TE), and an undefined population of transitional cells were identified. By comparing the transcriptome of blastocysts derived in vivo (IVV), in vitro from a conventional culture medium (IVC), and in vitro from an optimized reduced nutrient culture medium (IVR), we found a delay of the cell fate commitment to ICM in the IVC and IVR embryos. Developmental potential differences between IVV, IVC, and IVR embryos were mainly contributed by ICM and transitional cells. Pathway analysis of these non-TE cells between groups revealed highly active metabolic and biosynthetic processes, reduced cellular signaling, and reduced transmembrane transport activities in IVC embryos that may lead to reduced developmental potential. IVR embryos had lower activities in metabolic and biosynthetic processes but increased cellular signaling and transmembrane transport, suggesting these cellular mechanisms may contribute to improved blastocyst development compared to IVC embryos. However, the IVR embryos had compromised development compared to IVV embryos with notably over-active transmembrane transport activities that impaired ion homeostasis.

Energy metabolism disorders during in vitro maturation of bovine cumulus-oocyte complexes interfere with blastocyst quality and metabolism

Glucose and fatty acids (FA) metabolism disturbances during oocyte in vitro maturation (IVM) affect their metabolism and surrounding cumulus cells, but only inhibition of glucose metabolism decreases embryo culture efficiency. Therefore, the present experiment aimed to reveal if glucose or FA metabolism inhibition leads to the disruption of embryo developmental potential, and to characterize the metabolic landscape of embryos reaching the blastocyst stage. Inhibitors of glucose (IO + DHEA) or FA (ETOMOXIR) metabolism were applied during IVM, and the control group was matured under standard conditions. Blastocysts obtained from experimental and control groups were analyzed with regard to lipidome and metabolome (mass spectrometry), transcriptome (RNA-Seq) and fluorescence lipid droplets staining (BODIPY). We showed that inhibition of glucose and fatty acid metabolism leads to cellular stress response compromising the quality of preimplantation embryos. The inhibition of energy metabolism affects membrane fluidity as well as downregulates fatty acids biosynthesis and gene expression of trophectoderm cell line markers. Therefore, we conclude that oocyte maturation environment exerts a substantial effect on preimplantation development programming at cellular and molecular levels.

Biomarker metabolite mating of viable frozen-thawed in vitro-produced bovine embryos with pregnancy-competent recipients leads to improved birth rates

Selection of competent recipients before embryo transfer (ET) is indispensable for improving pregnancy and birth rates in cattle. However, pregnancy prediction can fail when the competence of the embryo is ignored. We hypothesized that the pregnancy potential of biomarkers could improve with information on embryonic competence. In vitro-produced embryos cultured singly for 24 h (from d 6 to 7) were transferred to d 7 synchronized recipients as fresh or after freezing and thawing. Recipient blood was collected on d 0 (estrus; n = 108) and d 7 (4-6 h before ET; n = 107) and plasma was analyzed by nuclear magnetic resonance (1H+NMR). Spent embryo culture medium (CM) was collected and analyzed by ultra-high-performance liquid chromatography tandem mass spectrometry in a subset of n = 70 samples. Concentrations of metabolites quantified in plasma (n = 35) were statistically analyzed as a function of pregnancy diagnosed on d 40, d 62 and birth. Univariate analysis with plasma metabolites consisted of a block study with controllable fixed factors (i.e., embryo cryopreservation, recipient breed, and day of blood collection; Wilcoxon test and t-test). Metabolite concentrations in recipients and embryos were independently analyzed by iterations that reclassified embryos or recipients using the support vector machine. Iterations identified some competent embryos, but mostly competent recipients that had a pregnancy incompetent partner embryo. Misclassified recipients that could be classified as competent were reanalyzed in a new iteration to improve the predictive model. After subsequent iterations, the predictive potential of recipient biomarkers was recalculated. On d 0, creatine, acetone and l-phenylalanine were the most relevant biomarkers at d 40, d 62, and birth, and on d 7, l-glutamine, l-lysine, and ornithine. Creatine was the most representative biomarker within blocks (n = 20), with a uniform distribution over pregnancy endpoints and type of embryos. Biomarkers showed higher abundance on d 7 than d 0, were more predictive for d 40 and d 62 than at birth, and the pregnancy predictive ability was lower with frozen-thawed (F-T) embryos. Six metabolic pathways differed between d 40 pregnant recipients for fresh and F-T embryos. Within F-T embryos, more recipients were misclassified, probably due to pregnancy losses, but were accurately identified when combined with embryonic metabolite signals. After recalculation, 12 biomarkers increased receiver operator characteristic-area under the curve (>0.65) at birth, highlighting creatine (receiver operator characteristic-area under the curve = 0.851), and 5 new biomarkers were identified. Combining metabolic information of recipient and embryos improves the confidence and accuracy of single biomarkers.

Pre-Implantation Bovine Embryo Evaluation-From Optics to Omics and Beyond

Approximately 80% of the ~1.5 million bovine embryos transferred in 2021 were in vitro produced. However, only ~27% of the transferred IVP embryos will result in live births. The ~73% pregnancy failures are partly due to transferring poor-quality embryos, a result of erroneous stereomicroscopy-based morphological evaluation, the current method of choice for pre-transfer embryo evaluation. Numerous microscopic (e.g., differential interference contrast, electron, fluorescent, time-lapse, and artificial-intelligence-based microscopy) and non-microscopic (e.g., genomics, transcriptomics, epigenomics, proteomics, metabolomics, and nuclear magnetic resonance) methodologies have been tested to find an embryo evaluation technique that is superior to morphologic evaluation. Many of these research tools can accurately determine embryo quality/viability; however, most are invasive, expensive, laborious, technically sophisticated, and/or time-consuming, making them futile in the context of in-field embryo evaluation. However accurate they may be, using complex methods, such as RNA sequencing, SNP chips, mass spectrometry, and multiphoton microscopy, at thousands of embryo production/collection facilities is impractical. Therefore, future research is warranted to innovate field-friendly, simple benchtop tests using findings already available, particularly from omics-based research methodologies. Time-lapse monitoring and artificial-intelligence-based automated image analysis also have the potential for accurate embryo evaluation; however, further research is warranted to innovate economically feasible options for in-field applications.

Review: Overview of the transcriptomic landscape in bovine blastocysts and elongated conceptuses driving developmental competence

In cattle, pregnancy loss due to early embryonic mortality is a major concern that significantly impacts reproductive efficiency. Given the economic importance of cattle in livestock productivity, much research has been carried out to comprehend the regulatory mechanisms underlying this early embryo loss. Thus, understanding the molecular principles behind the reciprocal communication between the maternal uterus and the developing conceptus is paramount. Measurement of mRNA expression through a variety of techniques is widely used to unravel the complex and dynamic interaction between these two players. Development of high-throughput technologies, such as microarrays and RNA sequencing, have allowed global quantification of the full range of expressed mRNA, or transcriptome, of a biological sample. Therefore, numerous investigators have applied one or the other method to study the bovine embryo transcriptome at different developmental checkpoints and under different conditions. The goal of this article was to review studies involving the use of high-throughput techniques to study the transcriptome of the bovine embryo from the blastocyst (∼day 7) to the elongating conceptus stage (∼days 13-16) in terms of developmental capacity and the impact of procedures for in vitro embryo production. Furthermore, the differentially expressed genes reported by each study and enriched pathways were compared to determine common terms. The studies described here highlight differences in the transcriptome (i) between blastocysts with divergent ability to sustain a pregnancy, (ii) between age-matched elongated conceptuses with divergent developmental fates, and (iii) between blastocysts and elongated conceptuses produced in vitro or in vivo. Comparison between these works, supported by other studies involving transcriptomic data integration presented at the end of this review, highlights the involvement of pathways related to energy metabolism in embryonic competence, which may be altered because of the procedures involved in the in vitro production of embryos.

Embryo morphokinetics derived from fresh and vitrified bovine oocytes predict blastocyst development and nuclear abnormalities

Embryo development is a dynamic process and critical stages may go unnoticed with the use of traditional morphologic assessments, especially the timing of embryonic divisions and aberrant zygotic cleavage patterns. Bovine embryo development is impaired after oocyte vitrification, but little is known about the underlying morphokinetic behavior. Here, bovine zygotes from fresh (n = 708) and vitrified oocytes (n = 182) were monitored by time-lapse imaging and the timing and nature of early blastomere divisions were modeled to find associations with blastocyst development at day 8. The predictive potential of morphokinetic parameters was analyzed by logistic regression and receiver operating characteristic curve analysis to determine optimal cut-off values. Lag-phase was highly correlated with embryo development. Remarkably, 100% of zygotes that reached the blastocyst stage showed a lag-phase. Fast first cleavage increased the chance of blastocyst development to 30% with a cut-off of 32 h and 22 min. Aberrant zygotic cleavage events, including multipolar division, unequal blastomere sizes, and membrane ruffling resulted in decreased blastocyst development. Multipolar division leads to uneven blastomeres, which was associated with anuclear and multinuclear blastomeres, indicating genome segregation errors. Moreover, we described for the first time morphokinetics of embryos derived from vitrified bovine oocytes. Vitrification severely affected blastocyst development, although lower cryoprotectant concentration in equilibration solutions seems to be less detrimental for embryo yield. Impaired development was linked to slow cleavages, lower lag-phase incidence, and increased early embryonic arrest. Typically, less than 15% of the embryos produced from vitrified oocytes reached more than eight cells. Interestingly, the rate of abnormal first cleavage events was not affected by oocyte vitrification. In conclusion, time to first cleavage, the presence of a lag-phase, and the absence of aberrant zygotic cleavage were the best predictors of bovine blastocyst development for both fresh and vitrified oocytes.

Machine-learning methods applied to integrated transcriptomic data from bovine blastocysts and elongating conceptuses to identify genes predictive of embryonic competence

Early pregnancy loss markedly impacts reproductive efficiency in cattle. The objectives were to model a biologically relevant gene signature predicting embryonic competence for survival after integrating transcriptomic data from blastocysts and elongating conceptuses with different developmental capacities and to validate the potential biomarkers with independent embryonic data sets through the application of machine-learning algorithms. First, two data sets from in vivo-produced blastocysts competent or not to sustain a pregnancy were integrated with a data set from long and short day-15 conceptuses. A statistical contrast determined differentially expressed genes (DEG) increasing in expression from a competent blastocyst to a long conceptus and vice versa; these were enriched for KEGG pathways related to glycolysis/gluconeogenesis and RNA processing, respectively. Next, the most discriminative DEG between blastocysts that resulted or did not in pregnancy were selected by linear discriminant analysis. These eight putative biomarker genes were validated by modeling their expression in competent or noncompetent blastocysts through Bayesian logistic regression or neural networks and predicting embryo developmental fate in four external data sets consisting of in vitro-produced blastocysts (i) competent or not, or (ii) exposed or not to detrimental conditions during culture, and elongated conceptuses (iii) of different length, or (iv) developed in the uteri of high- or subfertile heifers. Predictions for each data set were more than 85% accurate, suggesting that these genes play a key role in embryo development and pregnancy establishment. In conclusion, this study integrated transcriptomic data from seven independent experiments to identify a small set of genes capable of predicting embryonic competence for survival.

High-resolution ribosome profiling reveals translational selectivity for transcripts in bovine preimplantation embryo development

High-resolution ribosome fractionation and low-input ribosome profiling of bovine oocytes and preimplantation embryos has enabled us to define the translational landscapes of early embryo development at an unprecedented level. We analyzed the transcriptome and the polysome- and non-polysome-bound RNA profiles of bovine oocytes (germinal vesicle and metaphase II stages) and early embryos at the two-cell, eight-cell, morula and blastocyst stages, and revealed four modes of translational selectivity: (1) selective translation of non-abundant mRNAs; (2) active, but modest translation of a selection of highly expressed mRNAs; (3) translationally suppressed abundant to moderately abundant mRNAs; and (4) mRNAs associated specifically with monosomes. A strong translational selection of low-abundance transcripts involved in metabolic pathways and lysosomes was found throughout bovine embryonic development. Notably, genes involved in mitochondrial function were prioritized for translation. We found that translation largely reflected transcription in oocytes and two-cell embryos, but observed a marked shift in the translational control in eight-cell embryos that was associated with the main phase of embryonic genome activation. Subsequently, transcription and translation become more synchronized in morulae and blastocysts. Taken together, these data reveal a unique spatiotemporal translational regulation that accompanies bovine preimplantation development.

In vitro and in vivo-produced male dairy calves show molecular differences in the hepatic and muscular energy regulation

In cattle, the in vitro production (IVP) of embryos is becoming more relevant than embryos produced in vivo, i.e. after multiple ovulation and embryo transfer (MOET). However, the effects of IVP on the developmental programming of specific organs in the postnatal calves are yet unknown. Previously, we reported an epigenomic and transcriptomic profile of the hypothalamus–pituitary–testicular axis compatible with its earlier activation in IVP calves compared to MOET animals. Here, we studied the hepatic and muscular epigenome and transcriptome of those same male dairy calves (n = 4 per group). Tissue samples from liver and semitendinosus muscle were obtained at 3 months of age, and the extracted gDNA and RNA were sequenced through whole-genome bisulfite sequencing and RNA-sequencing, respectively. Next, bioinformatic analyses determined differentially methylated cytosines or differentially expressed genes [false discovery rate (FDR) < 0.05] for each Omic dataset; and nonparametrically combined genes (NPCG) for both integrated omics (P < 0.05). KEGG pathways enrichment analysis showed that NPCG upregulated in the liver and the muscle of the IVP calves were involved in oxidative phosphorylation and the tricarboxylic acid cycle. In contrast, ribosome and translation were upregulated in the liver but downregulated in the muscle of the IVP calves compared to the MOET calves (FDR < 0.05). A model considering the effect of the methylation levels and the group on the expression of all the genes involved in these pathways confirmed these findings. In conclusion, the multiomics data integration approach indicated an altered hepatic and muscular energy regulation in phenotypically normal IVP calves compared to MOET calves.

Environmental Exposures around Conception: Developmental Pathways Leading to Lifetime Disease Risk

Environment around conception can influence the developmental programme with lasting effects on gestational and postnatal phenotype and with consequences for adult health and disease risk. Peri-conception exposure comprises a crucial part of the 'Developmental Origins of Health and Disease' (DOHaD) concept. In this review, we consider the effects of maternal undernutrition experienced during the peri-conception period in select human models and in a mouse experimental model of protein restriction. Human datasets indicate that macronutrient deprivation around conception affect the epigenome, with enduring effects on cardiometabolic and neurological health. The mouse model, comprising maternal low protein diet exclusively during the peri-conception period, has revealed a stepwise progression in altered developmental programming following induction through maternal metabolite deficiency. This progression includes differential effects in extra-embryonic and embryonic cell lineages and tissues, leading to maladaptation in the growth trajectory and increased chronic disease comorbidities. The timeline embraces an array of mechanisms across nutrient sensing and signalling, cellular, metabolic, epigenetic and physiological processes with a coordinating role for mTORC1 signalling proposed. Early embryos appear active participants in environmental sensing to optimise the developmental programme for survival but with the trade-off of later disease. Similar adverse health outcomes may derive from other peri-conception environmental experiences, including maternal overnutrition, micronutrient availability, pollutant exposure and assisted reproductive treatments (ART) and support the need for preconception health before pregnancy.

Epigenomic and transcriptomic analyses reveal early activation of the HPG axis in in vitro-produced male dairy calves

In cattle, several calves born after IVP ("in vitro" embryo production) present similar birthweight to those generated after MOET (multiple ovulation and embryo transfer). However, the underlying molecular patterns in organs involved in the developmental process are unknown and could indicate physiological programming. The objectives of this study were: (1) to compare epigenomic and transcriptomic modifications in the hypothalamus, pituitary, gonadal and adrenal organs between 3 months old ovum pick-up-IVP and MOET male calves (n = 4 per group) and (2) to use blood epigenomic data to proxy methylation of the inner organs. Extracted gDNA and RNA were sequenced through whole-genome bisulfite sequencing and RNA sequencing, respectively. Next, bioinformatic analyses determined differentially methylated cytosines (DMC) and differentially expressed genes (DEG) (FDR < 0.05) in IVP versus MOET samples and the KEGG pathways that were overrepresented by genes associated with DMC or DEG (FDR < 0.1). Pathways related to hypothalamus, pituitary, gonadal (HPG) axis activation (GnRH secretion in the hypothalamus, GnRH signaling in the pituitary, and steroidogenesis in the testicle) were enriched in IVP calves. Modeling the effect of the methylation levels and the group on the expression of all the genes involved in these pathways confirmed their upregulation in HPG organs in IVP calves. The application of the DIABLO method allowed the identification of 15 epigenetic and five transcriptomic biomarkers, which were able to predict the embryo origin using the epigenomic data from the blood. In conclusion, the use of an integrated epigenomic-transcriptomic approach suggested an early activation of the HPG axis in male IVP calves compared to MOET counterparts, and the identification of potential biomarkers allowed the use of blood samples to proxy methylation levels of the relevant internal organs.