Deregulation of oxidative phosphorylation pathways in embryos derived in vitro from prepubertal and pubertal heifers based on whole-transcriptome sequencing

BACKGROUND

Although, oocytes from prepubertal donors are known to be less developmentally competent than those from adult donors it does not restrain their ability to produce full-term pregnancies. The transcriptomic profile of embryos could be used as a predictor for embryo's individual developmental competence. The aim of the study was to compare transcriptomic profile of blastocysts derived from prepubertal and pubertal heifers oocytes. Bovine cumulus-oocyte complexes (COCs) were obtained by ovum pick- up method from prepubertal and pubertal heifers. After in vitro maturation COCs were fertilized and cultured to the blastocyst stage. Total RNA was isolated from both groups of blastocysts and RNA-seq was performed. Gene ontology analysis was performed by DAVID (Database for Annotation, Visualization and Integrated Discovery).

RESULTS

A higher average blastocyst rate was obtained in the pubertal than in the pre-pubertal group. There were no differences in the quality of blastocysts between the examined groups. We identified 436 differentially expressed genes (DEGs) between blastocysts derived from researched groups, of which 247 DEGs were downregulated in blastocysts derived from pubertal compared to prepubertal heifers oocytes, and 189 DEGs were upregulated. The genes involved in mitochondrial function, including oxidative phosphorylation (OXPHOS) were found to be different in studied groups using Kyoto Encyclopedia of Genes (KEGG) pathway analysis and 8 of those DEGs were upregulated and 1 was downregulated in blastocysts derived from pubertal compared to prepubertal heifers oocytes. DEGs associated with mitochondrial function were found: ATP synthases (ATP5MF-ATP synthase membrane subunit f, ATP5PD- ATP synthase peripheral stalk subunit d, ATP12A- ATPase H+/K + transporting non-gastric alpha2 subunit), NADH dehydrogenases (NDUFS3- NADH: ubiquinone oxidoreductase subunit core subunit S3, NDUFA13- NADH: ubiquinone oxidoreductase subunit A13, NDUFA3- NADH: ubiquinone oxidoreductase subunit A3), cytochrome c oxidase (COX17), cytochrome c somatic (CYCS) and ubiquinol cytochrome c reductase core protein 1 (UQCRC1). We found lower number of apoptotic cells in blastocysts derived from oocytes collected from prepubertal than those obtained from pubertal donors.

CONCLUSIONS

Despite decreased expression of genes associated with OXPHOS pathway in blastocysts from prepubertal heifers oocytes, the increased level of ATP12A together with the lower number of apoptotic cells in these blastocysts might support their survival after transfer.

RNA sequencing and gene co-expression network of in vitro matured oocytes and blastocysts of buffalo

In reproductive technologies, uncovering the molecular aspects of oocyte and embryo competence under different conditions is crucial for refining protocols and enhancing efficiency. RNA-seq generates high-throughput data and provides transcriptomes that can undergo additional computational analyses. This study presented the transcriptomic profiles of in vitro matured oocytes and blastocysts produced in vitro from buffalo crossbred (Bubalus bubalis), coupled with gene co-expression and module preservation analysis. Cumulus Oophorus Complexes, obtained from slaughterhouse-derived ovaries, were subjected to in vitro maturation to yield metaphase II oocytes (616) or followed in vitro fertilization and culture to yield blastocysts for sequencing (526). Oocyte maturation (72%, ±3.34 sd) and embryo development (21.3%, ±4.18 sd) rates were obtained from three in vitro embryo production routines following standard protocols. Sequencing of 410 metaphase II oocytes and 70 hatched blastocysts (grade 1 and 2) identified a total of 13,976 genes, with 62% being ubiquitously expressed (8,649). Among them, the differentially expressed genes (4,153) and the strongly variable genes with the higher expression (fold-change above 11) were highlighted in oocytes (BMP15, UCHL1, WEE1, NLRPs, KPNA7, ZP2, and ZP4) and blastocysts (APOA1, KRT18, ANXA2, S100A14, SLC34A2, PRSS8 and ANXA2) as representative indicators of molecular quality. Additionally, genes exclusively found in oocytes (224) and blastocysts (2,200) with specific biological functions were identified. Gene co-expression network and module preservation analysis revealed strong preservation of functional modules related to exosome components, steroid metabolism, cell proliferation, and morphogenesis. However, cell cycle and amino acid transport modules exhibited weak preservation, which may reflect differences in embryo development kinetics and the activation of cell signaling pathways between buffalo and bovine. This comprehensive transcriptomic profile serves as a valuable resource for assessing the molecular quality of buffalo oocytes and embryos in future in vitro embryo production assays.

Comparative transcriptome profile of embryos at different developmental stages derived from somatic cell nuclear transfer (SCNT) and in-vitro fertilization (IVF) in riverine buffalo (Bubalus bubalis)

Somatic cell nuclear transfer (SCNT) is a very important reproductive technology with many diverse applications, such as fast multiplication of elite animals, the production of transgenic animals and embryonic stem (ES) cells. However, low cloning efficiency, a low live birth rate and the abnormally high incidence of abnormalities in the offspring born are attributed to incomplete or aberrant nuclear reprogramming. In SCNT embryos, the aberrant expression pattern of the genes throughout embryonic development is responsible for the incomplete nuclear reprogramming. The present study was carried out to identify the differential gene expression (DEGs) profile and molecular pathways of the SCNT and IVF embryos at different developmental stages (2 cell, 8 cell and blastocyst stages). In the present study, 1164 (2 cell), 1004 (8 cell) and 530 (blastocyst stage) DEGs were identified in the SCNT embryos as compared to IVF embryos. In addition, several genes such as ZEB1, GDF1, HSF5, PDE3B, VIM, TNNC, HSD3B1, TAGLN, ITGA4 and AGMAT were affecting the development of SCNT embryos as compared to IVF embryos. Further, Gene Ontology (GO) and molecular pathways analysis suggested, SCNT embryos exhibit variations compared to their IVF counterparts and affected the development of embryos throughout the different developmental stages. Apart from this, q-PCR analysis of the GDF1, TMEM114, and IGSF22 genes were utilized to validate the RNA-seq data. These findings contribute valuable insights about the different genes and molecular pathways underlying SCNT embryo development and offer crucial information for improving SCNT efficiency.

Vitamin C down-regulates the H3K9me3-dependent heterochromatin in buffalo fibroblasts via PI3K/PDK1/SGK1/KDM4A signal axis

The success of reprogramming is dependent on the reprogramming factors enriched in the cytoplasm of recipient oocytes and the potential of donor nucleus to be reprogrammed. Histone 3 lysine 9 trimethylation (H3K9me3) was identified as a major epigenetic barrier impeding complete reprogramming. Treating donor cell with vitamin C (Vc) can enhance the developmental potential of cloned embryos, but the underlying mechanisms still need to be elucidated. In this study, we found that 20μg/mL Vc could promote proliferation and inhibit apoptosis of BFFs, as well as down-regulate the H3K9me3-dependent heterochromatin and increase chromatin accessibility. Inhibited the expression of KDM4A resulted in increasing apoptosis rate and the H3K9me3-dependent heterochromatin, which can be restored by Vc. Moreover, Vc up-regulated the expression of KDM4A through PI3K/PDK1/SGK1 pathway. Inhibiting any factor in the signal axis of this PI3K pathway not only suppressed the activity of KDM4A but also substantially increased the level of H3K9me3 modification and the expression of the HP1α protein. Finally, Vc can rescue those negative effects induced by the blocking the PI3K/PDK1/SGK1 pathway. Collectively, Vc can down-regulate the H3K9me3-dependent heterochromatin in BFFs via PI3K/PDK1/SGK1/KDM4A signal axis, suggesting that Vc can turn the chromatin status of donor cells to be reprogrammed more easily.

Early Cell Specification in Mammalian Fertilized and Somatic Cell Nuclear Transfer Embryos

Early cell specification in mammalian preimplantation embryos is an intricate cellular process that leads to coordinated spatial and temporal expression of specific genes. Proper segregation into the first two cell lineages, the inner cell mass (ICM) and the trophectoderm (TE), is imperative for developing the embryo proper and the placenta, respectively. Somatic cell nuclear transfer (SCNT) allows the formation of a blastocyst containing both ICM and TE from a differentiated cell nucleus, which means that this differentiated genome must be reprogrammed to a totipotent state. Although blastocysts can be generated efficiently through SCNT, the full-term development of SCNT embryos is impaired mostly due to placental defects. In this review, we examine the early cell fate decisions in fertilized embryos and compare them to observations in SCNT-derived embryos, in order to understand if these processes are affected by SCNT and could be responsible for the low success of reproductive cloning.

Global DNA methylation profiles of buffalo (Bubalus bubalis) preimplantation embryos produced by handmade cloning and in vitro fertilization

Somatic cell nuclear transfer technique (SCNT) has proved to be an outstanding method of multiplication of elite animals but accompanied with low efficiency and live birth rate of cloned animals. Epigenetic alterations of DNA has been one of the culprits behind this issue. Cloned embryos are found to deviate slightly from regular pattern of demethylation and re-methylation at the time of nuclear reprogramming and embryonic development when compared with embryos produced by in vitro fertilization (IVF). Thus, the present study was aimed at evaluating global DNA methylation profiles of cloned embryos at 2-cell, 8-cell and blastocyst stages and compare it with corresponding stages of embryos produced by IVF by using MeDIP-Sequencing on Illumina-based platform. We found out that cloned embryos exhibited significantly different DNA methylation pattern as compared to IVF embryos with respect to distribution of differentially methylated regions in different components of genome, CpG islands distribution and methylation status, gene ontological profiles and pathways affected throughout the developmental stages. The data generated from MeDIP-Seq was validated at blastocyst stage cloned and IVF embryos by bisulfite-sequencing PCR on five randomly selected gene regions.

Global MicroRNA Expression Profiling of Buffalo (Bubalus bubalis) Embryos at Different Developmental Stages Produced by Somatic Cell Nuclear Transfer and In-Vitro Fertilization Using RNA Sequencing

Despite the success of cloning technology in the production of offspring across several species, its application on a wide scale is severely limited by the very low offspring rate obtained with cloned embryos. The expression profile of microRNAs (miRNAs) in cloned embryos throughout embryonic development is reported to deviate from regular patterns. The present study is aimed at determining the dynamics of the global expression of miRNA profile in cloned and in-vitro fertilization (IVF) pre-implantation embryos at different developmental stages, i.e., the two-cell, eight-cell, and blastocyst stages, using next-generation sequencing. The results of this study suggest that there is a profound difference in global miRNA profile between cloned and IVF embryos. These differences are manifested throughout the course of embryonic development. The cloned embryos differ from their IVF counterparts in enriched Gene Ontology (GO) terms of biological process, molecular function, cellular component, and protein class categories in terms of the targets of differentially expressed miRNAs. The major pathways related to embryonic development, such as the Wnt signaling pathway, the apoptosis signaling pathway, the FGF signaling pathway, the p53 pathway, etc., were found to be affected in cloned relative to IVF embryos. Overall, these data reveal the distinct miRNA profile of cloned relative to IVF embryos, suggesting that the molecules or pathways affected may play an important role in cloned embryo development.

MeDIP-sequencing for profiling global DNA methylation in buffalo embryos produced by in vitro fertilization

Assisted reproductive technique like in vitro fertilization has contributed immensely in producing genetically improved livestock. Production of embryos under in vitro conditions can affect global DNA methylation pattern during the course of embryonic development. The present study is aimed at the generation and comparison of global DNA methylome of embryos at 2-cell, 8-cell and blastocyst stage of buffalo embryos produced by in vitro fertilization using MeDIP-Sequencing. It is observed that there is a profound difference in the global DNA methylation profile of IVF embryos at different developmental stages. These differences are manifested throughout the course of embryonic development. Pathways like Wnt signaling pathway, gonadotropin-releasing hormone receptor pathway and integrin signaling were found to be majorly affected by hypermethylation of DNA in IVF embryos throughout the development.

Effect of Dickkopf-1 and colony stimulating factor-2 on the developmental competence, quality, gene expression and live birth rate of buffalo (Bubalus bubalis) embryos produced by hand-made cloning

A functional canonical WNT signaling pathway exists in preimplantation embryos and inhibits embryonic development. Recent studies suggest that this pathway is over-expressed in nuclear transferred (NT), compared to IVF embryos. The present study investigated the effects of Dickkopf-1 (DKK1), an inhibitor of canonical WNT signaling pathway and colony stimulating factor-2 (CSF2), an embryokine, on the developmental competence, quality, gene expression and live birth rate of NT buffalo embryos produced by Hand-made cloning (HMC). Following supplementation of the in vitro culture medium on day 5 with DKK1 (100 ng/mL), CSF2 (10 ng/mL), DKK1+CSF2 or no supplementation (control), the blastocyst rate was higher (P < 0.05) with DKK1 and DKK1+CSF2 (42.6 ± 1.4% and 46.6 ± 0.9%, respectively) than with CSF2 or controls (40.6 ± 1.3% and 39.0 ± 1.3%, respectively). The apoptotic index of the blastocysts was lower (P < 0.05) for DKK1, CSF2 and DKK1+CSF2 groups (3.44 ± 0.14, 3.39 ± 0.11 and 3.11 ± 0.22, respectively) compared to controls (6.64 ± 0.25), and was similar to that of the IVF blastocysts (3.67 ± 0.18). Although the total cell number was similar for the DKK1, CSF2, DKK1+CSF2 and control groups (200.4 ± 3.05, 196.4 ± 3.73, 204.7 ± 3.71 and 205 ± 4.03, respectively), the inner cell mass:trophectoderm cell number ratio of DKK1, CSF2 and DKK1+CSF2 groups (0.21 ± 0.01, 0.17 ± 0.01 and 0.22 ± 0.02, respectively) was higher (P < 0.05) than controls (0.13 ± 0.01) and was similar to that of IVF blastocysts (0.19 ± 0.01). Treatment with DKK1 or CSF2 or both increased (P < 0.05) the expression level of OCT4, NANOG,SOX2, GATA6, BCL2, PTEN, P53, FGF4, GLUT1 and IFN-τ, and decreased that of C-MYC, CDX2, CASPASE, DNMT3a, TCF7 and LEF1 in blastocysts, compared to controls. Transfer of DKK1-treated embryos to 13 recipients resulted in 4 pregnancies (30.8%; 2 live births, one abortion and one currently at 9 months of pregnancy) whereas, transfer of DKK1+CSF2-treated embryos to 16 recipients, resulted in 4 pregnancies (25.0%), all of which resulted in live births. No pregnancy was obtained after transfer of control and CSF-treated embryos to 12 and 16 recipients, respectively. These results suggest that DKK1 treatment of NT embryos increases the blastocyst, conception and live birth rate, and improves their quality whereas, CSF2 treatment, does not affect the blastocyst, conception and live birth rate despite improvement in embryo quality.

Effect of season on the in-vitro maturation and developmental competence of buffalo oocytes after somatic cell nuclear transfer

Somatic cell nuclear transfer (SCNT) is a valuable technology tool with various uses in transgenic animals, regenerative medicine, and stem cell research. However, the efficiency of SCNT embryos appears to have poor developmental competency. Environmental issues may adversely affect SCNT embryos in buffalo. Thereafter, the present study aimed to explore the effect of season on the maturation of buffalo oocytes and subsequent developmental capability after parthenogenetic activation and SCNT in buffalo. Buffalo oocytes (n = 6353) were collected from local slaughterhouse at various seasons; spring (March-April), summer (May-August), autumn (September-November), and winter (December-January). A significant increase (p < 0.05) was recorded in the maturation rate (57.07%) at autumn compared with spring, summer, and winter (50.46, 50.93, and 50.66%, respectively). No significant differences were recorded in the fusion and the cleavage rates among all seasons. Blastocyst development rate was higher (p < 0.05) in autumn and winter (16.52 ± 8.45% and 15.98 ± 7.17%, respectively) than in spring and summer (9.47 ± 6.71% and 10.84 ± 6.58%, respectively) seasons. It could be concluded that the season had a significant effect on oocyte development competence which can be used for SCNT in buffalo.