Profiling the transcriptomic signatures and identifying the patterns of zygotic genome activation - a comparative analysis between early porcine embryos and their counterparts in other three mammalian species

The first lineage determination in mammals

This review describes in detail the morphological, cytoskeletal and gene expression events leading to the gene regulatory network bifurcation point of trophoblast and inner cell mass cells in a variety of mammalian preimplantation embryos. The interrelated processes of compaction and polarity establishment are discussed in terms of how they affect YAP/WWTR activity and the location and fate of cells. Comparisons between mouse, human, cattle, pig and rabbit embryos suggest a conserved role for YAP/WWTR signalling in trophoblast induction in eutherian animals though the mechanisms for, and timing of, YAP/WWTR activation differs among species. Downstream targets show further differences, with the trophoblast marker GATA3 being a direct target in all examined mammals, while CDX2-positive and SOX2-negative regulation varies.

Longitudinal profiling of human androgenotes through single-cell analysis unveils paternal gene expression dynamics in early embryo development

STUDY QUESTION

How do transcriptomics vary in haploid human androgenote embryos at single cell level in the first four cell cycles of embryo development?

SUMMARY ANSWER

Gene expression peaks at the fourth cell cycle, however some androcytes exhibit unique transcriptional behaviors.

WHAT IS KNOWN ALREADY

The developmental potential of an embryo is determined by the competence of the oocyte and the sperm. However, studies of the contribution of the paternal genome using pure haploid androgenotes are very scarce.

STUDY DESIGN, SIZE, DURATION

This study was performed analyzing the single-cell transcriptomic sequencing of 38 androcytes obtained from 10 androgenote bioconstructs previously produced in vitro (de Castro et al., 2023). These results were analyzed through different bioinformatics software such as g: Profiler, GSEA, Cytoscape, and Reactome.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Single cell sequencing was used to obtain the transcriptomic profiles of the different androcytes. The results obtained were compared between the different cycles studied using the DESeq2 program and functional enrichment pathways using g: Profiler, Cytoscape, and Reactome.

MAIN RESULTS AND THE ROLE OF CHANCE

A wave of paternally driven transcriptomic activation was found during the third-cell cycle, with 1128 upregulated and 225 downregulated genes and the fourth-cell cycle, with 1373 upregulated and 286 downregulated genes, compared to first-cell cycle androcytes. Differentially expressed routes related to cell differentiation, DNA-binding transcription, RNA biosynthesis and RNA polymerase II transcription regulatory complex, and cell death were found in the third and fourth with respect to the first-cell cycle. Conversely, in the fourth cell cycle, 153 downregulated and 332 upregulated genes were found compared with third cell cycle, associated with differentially expressed processes related to E-box binding and zinc finger protein 652 (ZNF652) transcription factor. Further, significant overexpression of LEUTX, PRAMEF1, DUXA, RFPL4A, TRIM43, and ZNF675 found in androgenotes, compared to biparental embryos, highlights the paternal contributions to zygote genome activation.

LARGE SCALE DATA

All raw sequencing data are available through the Gene Expression Omnibus (GEO) under accessions number: GSE216501.

LIMITATIONS, REASONS FOR CAUTION

Extrapolation of biological events from uniparental constructs to biparental embryos should be done with caution. Maternal and paternal genomes do not act independently of each other in a natural condition. The absence of one genome may affect gene transcription of the other. In this sense, the haploid condition of the bioconstructs could mask the transcriptomic patterns of the single cells.

WIDER IMPLICATIONS OF THE FINDINGS

The results obtained demonstrated the level of involvement of the human paternal haploid genome in the early stages of embryo development as well as its evolution at the transcriptomic level, laying the groundwork for the use of these bioconstructs as reliable models to dispel doubts about the genetic role played by the paternal genome in the early cycles of embryo development.

STUDY FUNDING/COMPETING INTEREST(S)

This study was funded by Instituto de Salud Carlos III (ISCIII) through the project 'PI22/00924', co-funded by European Regional Development Fund (ERDF); 'A way to make Europe'. F.D. was supported by the Spanish Ministry of Economy and Competitiveness through the Miguel Servet program (CPII018/00002). M.J.E. was supported by Instituto de Salud Carlos III (PI19/00577 [M.J.E.]) and FI20/00086. P.dC. was supported by a predoctoral grant for training in research into health (PFIS PI19/00577) from the Instituto de Salud Carlos III. All authors declare having no conflict of interest with regard to this trial.

DNA methylation regulates RNA m6A modification through transcription factor SP1 during the development of porcine somatic cell nuclear transfer embryos

Epigenetic modifications play critical roles during somatic cell nuclear transfer (SCNT) embryo development. Whether RNA N6-methyladenosine (m6A) affects the developmental competency of SCNT embryos remains unclear. Here, we showed that porcine bone marrow mesenchymal stem cells (pBMSCs) presented higher RNA m6A levels than those of porcine embryonic fibroblasts (pEFs). SCNT embryos derived from pBMSCs had higher RNA m6A levels, cleavage, and blastocyst rates than those from pEFs. Compared with pEFs, the promoter region of METTL14 presented a hypomethylation status in pBMSCs. Mechanistically, DNA methylation regulated METTL14 expression by affecting the accessibility of transcription factor SP1 binding, highlighting the role of the DNA methylation/SP1/METTL14 pathway in donor cells. Inhibiting the DNA methylation level in donor cells increased the RNA m6A level and improved the development efficiency of SCNT embryos. Overexpression of METTL14 significantly increased the RNA m6A level in donor cells and the development efficiency of SCNT embryos, whereas knockdown of METTL14 suggested the opposite result. Moreover, we revealed that RNA m6A-regulated TOP2B mRNA stability, translation level, and DNA damage during SCNT embryo development. Collectively, our results highlight the crosstalk between RNA m6A and DNA methylation, and the crucial role of RNA m6A during nuclear reprogramming in SCNT embryo development.

Single-Cell RNA Sequencing Reveals Differences in Chromatin Remodeling and Energy Metabolism among In Vivo-Developed, In Vitro-Fertilized, and Parthenogenetically Activated Embryos from the Oocyte to 8-Cell Stages in Pigs

In vitro-fertilized (IVF) and parthenogenetically activated (PA) embryos, key to genetic engineering, face more developmental challenges than in vivo-developed embryos (IVV). We analyzed single-cell RNA-seq data from the oocyte to eight-cell stages in IVV, IVF, and PA porcine embryos, focusing on developmental differences during early zygotic genome activation (ZGA), a vital stage for embryonic development. (1) Our findings reveal that in vitro embryos (IVF and PA) exhibit more similar developmental trajectories compared to IVV embryos, with PA embryos showing the least gene diversity at each stage. (2) Significant differences in maternal mRNA, particularly affecting mRNA splicing, energy metabolism, and chromatin remodeling, were observed. Key genes like SMARCB1 (in vivo) and SIRT1 (in vitro) played major roles, with HDAC1 (in vivo) and EZH2 (in vitro) likely central in their complexes. (3) Across different types of embryos, there was minimal overlap in gene upregulation during ZGA, with IVV embryos demonstrating more pronounced upregulation. During minor ZGA, global epigenetic modification patterns diverged and expanded further. Specifically, in IVV, genes, especially those linked to H4 acetylation and H2 ubiquitination, were more actively regulated compared to PA embryos, which showed an increase in H3 methylation. Additionally, both types displayed a distinction in DNA methylation. During major ZGA, IVV distinctively upregulated genes related to mitochondrial regulation, ATP synthesis, and oxidative phosphorylation. (4) Furthermore, disparities in mRNA degradation-related genes between in vivo and in vitro embryos were more pronounced during major ZGA. In IVV, there was significant maternal mRNA degradation. Maternal genes regulating phosphatase activity and cell junctions, highly expressed in both in vivo and in vitro embryos, were degraded in IVV in a timely manner but not in in vitro embryos. (5) Our analysis also highlighted a higher expression of many mitochondrially encoded genes in in vitro embryos, yet their nucleosome occupancy and the ATP8 expression were notably higher in IVV.

TET Family Members Are Integral to Porcine Oocyte Maturation and Parthenogenetic Pre-Implantation Embryogenesis

The ten-eleven translocation (TET) enzyme family, which includes TET1/2/3, participates in active DNA demethylation in the eukaryotic genome; moreover, TET1/2/3 are functionally redundant in mice embryos. However, the combined effect of TET1/2/3 triple-gene knockdown or knockout on the porcine oocytes or embryos is still unclear. In this study, using Bobcat339, a specific small-molecule inhibitor of the TET family, we explored the effects of TET enzymes on oocyte maturation and early embryogenesis in pigs. Our results revealed that Bobcat339 treatment blocked porcine oocyte maturation and triggered early apoptosis. Furthermore, in the Bobcat339-treated oocytes, spindle architecture and chromosome alignment were disrupted, probably due to the huge loss of 5-hydroxymethylcytosine (5hmC)and concurrent increase in 5-methylcytosine (5mC). After Bobcat339 treatment, early parthenogenetic embryos exhibited abnormal 5mC and 5hmC levels, which resulted in compromised cleavage and blastocyst rate. The mRNA levels of EIF1A and DPPA2 (ZGA marker genes) were significantly decreased, which may explain why the embryos were arrested at the 4-cell stage after Bobcat339 treatment. In addition, the mRNA levels of pluripotency-related genes OCT4 and NANOG were declined after Bobcat339 treatment. RNA sequencing analysis revealed differentially expressed genes in Bobcat339-treated embryos at the 4-cell stage, which were significantly enriched in cell proliferation, cell component related to mitochondrion, and cell adhesion molecule binding. Our results indicated that TET proteins are essential for porcine oocyte maturation and early embryogenesis, and they act by mediating 5mC/5hmC levels and gene transcription.

Whole-genome transcriptome and DNA methylation dynamics of pre-implantation embryos reveal progression of embryonic genome activation in buffaloes

BACKGROUND

During mammalian pre-implantation embryonic development (PED), the process of maternal-to-zygote transition (MZT) is well orchestrated by epigenetic modification and gene sequential expression, and it is related to the embryonic genome activation (EGA). During MZT, the embryos are sensitive to the environment and easy to arrest at this stage in vitro. However, the timing and regulation mechanism of EGA in buffaloes remain obscure.

RESULTS

Buffalo pre-implantation embryos were subjected to trace cell based RNA-seq and whole-genome bisulfite sequencing (WGBS) to draw landscapes of transcription and DNA-methylation. Four typical developmental steps were classified during buffalo PED. Buffalo major EGA was identified at the 16-cell stage by the comprehensive analysis of gene expression and DNA methylation dynamics. By weighted gene co-expression network analysis, stage-specific modules were identified during buffalo maternal-to-zygotic transition, and key signaling pathways and biological process events were further revealed. Programmed and continuous activation of these pathways was necessary for success of buffalo EGA. In addition, the hub gene, CDK1, was identified to play a critical role in buffalo EGA.

CONCLUSIONS

Our study provides a landscape of transcription and DNA methylation in buffalo PED and reveals deeply the molecular mechanism of the buffalo EGA and genetic programming during buffalo MZT. It will lay a foundation for improving the in vitro development of buffalo embryos.

Transcription factor ELK1 regulates the expression of histone 3 lysine 9 to affect developmental potential of porcine preimplantation embryos

A series of changes occur in the early embryo that are critical for subsequent development, and the pig is an excellent animal model of human disease, so understanding the regulatory mechanisms of early embryonic development in the pig is of very importance. To find key transcription factors regulating pig early embryonic development, we first profiled the transcriptome of pig early embryos, and confirmed that zygotic gene activation (ZGA) in porcine embryos starts from 4 cell stage. Subsequent enrichment analysis of up-regulated gene motifs during ZGA revealed that the transcription factor ELK1 ranked first. The expression pattern of ELK1 in porcine early embryos was analyzed by immunofluorescence staining and qPCR, and the results showed that the transcript level of ELK1 reached the highest at the 8 cell stage, while the protein level reached the highest at 4 cell stage. To further investigate the effect of ELK1 on early embryo development in pigs, we silenced ELK1 in zygotes and showed that ELK1 silencing significantly reduced cleavage rate, blastocyst rate as well as blastocyst quality. A significant decrease in the expression of the pluripotency gene Oct4 was also observed in blastocysts from the ELK1 silenced group by immunofluorescence staining. Silencing of ELK1 also resulted in decreased H3K9Ac modification and increased H3K9me3 modification at 4 cell stage. To investigate the effect of ELK1 on ZGA, we analyzed transcriptome changes in 4 cell embryos after ELK1 silencing by RNA seq, which revealed that ELK1 silencing resulted in significant differences in the expression of a total of 1953 genes at the 4 cell stage compared with their normal counterparts, including 1106 genes that were significantly upregulated and 847 genes that were significantly downregulated. Through GO and KEGG enrichment, we found that the functions and pathways of down-regulated genes were concentrated in protein synthesis, processing, cell cycle regulation, etc., while the functions of up-regulated genes were focused on aerobic respiration process. In conclusion, this study demonstrates that the transcription factor ELK1 plays an important role in regulation of preimplantation embryo development of pigs and deficiency of ELK1 leads to abnormal epigenetic reprogramming as well as zygotic genome activation, thus adversely affecting embryonic development. This study will provide important reference for the regulation of transcription factors in porcine embryo development.

Species and embryo genome origin affect lipid droplets in preimplantation embryos

Mammalian embryo development is affected by multiple metabolism processes, among which energy metabolism seems to be crucial. Therefore the ability and the scale of lipids storage in different preimplantation stages might affect embryos quality. The aim of the present studies was to show a complex characterization of lipid droplets (LD) during subsequent embryo developmental stages. It was performed on two species (bovine and porcine) as well as on embryos with different embryo origin [after in vitro fertilization (IVF) and after parthenogenetic activation (PA)]. Embryos after IVF/PA were collected at precise time points of development at the following stages: zygote, 2-cell, 4-cell, 8/16-cell, morula, early blastocyst, expanded blastocyst. LD were stained with BODIPY 493/503 dye, embryos were visualized under a confocal microscope and images were analyzed with the ImageJ Fiji software. The following parameters were analyzed: lipid content, LD number, LD size and LD area within the total embryo. The most important results show that lipid parameters in the IVF vs. PA bovine embryos differ at the most crucial moments of embryonic development (zygote, 8-16-cell, blastocyst), indicating possible dysregulations of lipid metabolism in PA embryos. When bovine vs. porcine species are compared, we observe higher lipid content around EGA stage and lower lipid content at the blastocyst stage for bovine embryos, which indicates different demand for energy depending on the species. We conclude that lipid droplets parameters significantly differ among developmental stages and between species but also can be affected by the genome origin.

Treatments of Porcine Nuclear Recipient Oocytes and Somatic Cell Nuclear Transfer-Generated Embryos with Various Reactive Oxygen Species Scavengers Lead to Improvements of Their Quality Parameters and Developmental Competences by Mitigating Oxidative Stress-Related Impacts

This study investigated the antioxidant effects of β-cryptoxanthin (BCX), hesperetin (HES), and icariin (ICA), and their effects on in vitro maturation of porcine oocytes and subsequent embryonic development of somatic cell nuclear transfer (SCNT). Treatment with 1 μM BCX (BCX-1) increased the developmental rate of porcine oocytes more than treatment with 100 μM HES (HES-100) or 5 μM ICA (ICA-5). The glutathione level and mRNA expression of antioxidant genes (NFE2L2, SOD1, and SOD2) were more increased in the BCX-1 group than in the HES-100 and ICA-5 groups, while the reactive oxygen species level was more decreased. Moreover, BCX improved the developmental capacity and quality of SCNT embryos. The total cell number, apoptotic cell rate, and development-related gene expression were modulated in the BCX-1 group to enhance embryonic development of SCNT. These results show that the antioxidant effects of BCX enhance in vitro maturation of porcine oocytes and subsequent embryonic development of SCNT.