Chromatin remodeling in bovine embryos indicates species-specific regulation of genome activation

Functional COPA is indispensable for early embryo development beyond major genome activation in bovines

Embryonic genome activation is divided into a minor and a major wave of transition to endogenous transcription. In bovines, minor genome activation begins early in the 2-cell stage and is completed by the 8-cell stage when major genome activation becomes dominant. While the activation of genes known to regulate early development have been studied extensively, genes involved in more central cellular functions have not been examined. Taking advantage of the CRISPR Cas9 system, the present study investigated the effect of knocking out the Golgi retrograde protein transporter COPA on early bovine development. After the electroporation of presumptive zygotes with Cas9 ribonucleoproteins targeting COPA exon 6, sequences of 2 (11 %) and 4-cell (16 %) embryos showed knockouts of COPA whereas 8-cell embryos and blastocysts did not, demonstrating that COPA is necessary for development to the 8-cell stage and beyond. Using a repair template containing silent mutations along the target site, COPA loss of wildtype was observed in 5 blastocysts, with successful knock-in of the template on at least one allele. This shows that an edited yet functional copy of COPA can save the developmental capacity of the embryo and demonstrates that Cas9 activity at the target region itself is not responsible for the loss of function. Together, the present study revealed that COPA is necessary for embryonic development, and that the timing of this necessity is before major genome activation onset. More generally, this study further demonstrates the utility of genome editing within reproductive biotechnology for the interrogation of gene function and early embryonic development.

Multiple Retrotransposon-mediated NF-YA Gene Duplication Events Recurred in Diverse Groups of Mammals at Different Ancestry Levels

NF-Y is a transcription factor trimer formed by Histone Fold subunits NF-YB/NF-YC and NF-YA, which confers sequence-specificity for the CCAAT box, an important cis regulatory element. The subunits are extremely conserved in all eukaryotes and in mammals they are typically encoded by single copy genes. We describe here the presence of a second NF-YA termed NF-YAr for retrogene in diverse groups of mammals (Cetacea, Ruminantia, Ursidae, Sciuridae, hippopotamus, and greater horseshoe bat). NF-YAr retrogenes are located on different chromosomes with respect to the parental gene; they are compact and intronless, or with few annotated introns. Phylogenetic and synteny analyses indicate multiple independent retrotransposition events in the different orders. Analysis of RNA-seq data of Bos taurus suggests expression confined to spermatozoa. Conservation of translation initiation signals around predicted start codons, and of 5'UTR sequences, are consistent with protein expression, suggesting that NF-YAr is a translated, retroposed NF-YA. 3D-informed structural considerations of the predicted protein sequences point at deleterious changes for CCAAT-binding and, potentially, for trimer formation. These findings indicate that multiple independent NF-YA retrotransposition events were fixed in selected orders of mammals, generating a second NF-YA with a strict tissue distribution.

Super RNA Pol II domains enhance minor ZGA through 3D interaction to ensure the integrity of major transcriptional waves in late-ZGA mammals

Zygotic genome activation (ZGA) occurs at distinct stages across mammals, with mice initiating ZGA at the 2-cell stage and bovines and humans activating the process in the 4- to 8-cell stages. RNA polymerase II (RNA Pol II) gradually initiates ZGA in mice, but regulation in late-ZGA species remains unclear. Here, RNA Pol II profiling in bovine embryos identified strong intergenic clusters that boost minor ZGA gene expression via chromatin interactions and are named super RNA Pol II domains (SPDs). CRISPRi perturbation of SPDs in bovine embryos decreases the expression of minor ZGA genes, whereas the knockdown of these genes disrupts major ZGA and embryogenesis. Rapid enhancement of minor ZGA genes also occurs in human embryos. Alternatively, mouse and porcine oocytes precociously express these minor ZGA genes without SPDs. Thus, SPDs appear to be an adaptation in bovine embryos, promoting minor ZGA gene expression to comparable levels as early-ZGA species, illuminating species-specific regulation of ZGA timing.

Genome-wide association analyses identify single-nucleotide polymorphisms associated with in vitro embryo cleavage and blastocyst rates in Holstein bulls

Reproductive success is an essential component of profitable and sustainable dairy operations. Although selection for production traits such as milk yield has led to a decline in fertility in dairy cattle, strategies including assisted reproductive technologies such as in vitro fertilization and embryo transfer, as well as genomic selection for fertility traits, are being implemented to help mitigate this loss. Previous studies have identified genetic markers associated with fertility traits such as daughter pregnancy rate, conception rate, and interval to first conception. However, genetic markers associated with in vitro embryo development have yet to be explored. Therefore, this study aimed to identify genetic markers associated with in vitro embryo cleavage and blastocyst rates with the future goal of creating an index to select bulls with superior abilities to produce embryos in vitro. Rigorous in vitro fertilization trials identified Holstein bulls with divergent abilities to produce embryos that successfully cleave and develop into blastocysts in vitro. A total of 40 bulls with embryo cleavage and blastocyst rate phenotypes were whole genome sequenced. Sequencing was performed on an Illumina platform with PE150 reads followed by cleaning, mapping to ARS-UCD1.2, variant calling with GATK HaplotypeCaller, and quality control following the 1000 Bulls Genomes Project best practices. Genome-wide association analyses identified 819 and 442 SNPs associated with embryo cleavage and blastocyst rates, respectively. Significant regions on chromosomes 15, 18, 21, 22, and 23 were in linkage disequilibrium with QTL previously associated with reproductive traits in cattle. Further, a region on chromosome 28 with the most significant variant associated with cleavage rate contained several synonymous variants in EGLN1, which is a critical component of the hypoxia inducible pathway. An additional region on chromosome 18 associated with cleavage rate contained a missense variant in SMG9, a gene involved in the nonsense-mediated mRNA decay pathway. A region on chromosome 21 associated with blastocyst rate also contained variants in regulatory regions downstream from AEN, a gene required for efficient DNA fragmentation during the p53-dependent apoptosis pathway. In summary, this study identified a preliminary collection of genomic regions associated with embryo cleavage and blastocyst rates in vitro that contain relevant genes and other QTL associated with reproductive traits in cattle. These regions, after validation, could contribute to a selection index for identification of bulls with genetic predispositions to produce a greater number of embryos in vitro for use in assisted reproduction techniques.

Epigenome dynamics in early mammalian embryogenesis

During early embryonic development in mammals, the totipotency of the zygote - which is reprogrammed from the differentiated gametes - transitions to pluripotency by the blastocyst stage, coincident with the first cell fate decision. These changes in cellular potency are accompanied by large-scale alterations in the nucleus, including major transcriptional, epigenetic and architectural remodelling, and the establishment of the DNA replication programme. Advances in low-input genomics and loss-of-function methodologies tailored to the pre-implantation embryo now enable these processes to be studied at an unprecedented level of molecular detail in vivo. Such studies have provided new insights into the genome-wide landscape of epigenetic reprogramming and chromatin dynamics during this fundamental period of pre-implantation development.

The Complexities of Interspecies Somatic Cell Nuclear Transfer: From Biological and Molecular Insights to Future Perspectives

For nearly three decades, interspecies somatic cell nuclear transfer (iSCNT) has been explored as a potential tool for cloning, regenerative medicine, and wildlife conservation. However, developmental inefficiencies remain a major challenge, largely due to persistent barriers in nucleocytoplasmic transport, mitonuclear communication, and epigenome crosstalk. This review synthesized peer-reviewed English articles from PubMed, Web of Science, and Scopus, spanning nearly three decades, using relevant keywords to explore the molecular mechanisms underlying iSCNT inefficiencies and potential improvement strategies. We highlight recent findings deepening the understanding of interspecies barriers in iSCNT, emphasizing their interconnected complexities, including the following: (1) nucleocytoplasmic incompatibility may disrupt nuclear pore complex (NPC) assembly and maturation, impairing the nuclear transport of essential transcription factors (TFs), embryonic genome activation (EGA), and nuclear reprogramming; (2) mitonuclear incompatibility could lead to nuclear and mitochondrial DNA (nDNA-mtDNA) mismatches, affecting electron transport chain (ETC) assembly, oxidative phosphorylation, and energy metabolism; (3) these interrelated incompatibilities can further influence epigenetic regulation, potentially leading to incomplete epigenetic reprogramming in iSCNT embryos. Addressing these challenges requires a multifaceted, species-specific approach that balances multiple incompatibilities rather than isolating a single factor. Gaining insight into the molecular interactions between the donor nucleus and recipient cytoplast, coupled with optimizing strategies tailored to specific pairings, could significantly enhance iSCNT efficiency, ultimately transforming experimental breakthroughs into real-world applications in reproductive biotechnology, regenerative medicine, and species conservation.

The maternal-to-zygotic transition: reprogramming of the cytoplasm and nucleus

A fertilized egg is initially transcriptionally silent and relies on maternally provided factors to initiate development. For embryonic development to proceed, the oocyte-inherited cytoplasm and the nuclear chromatin need to be reprogrammed to create a permissive environment for zygotic genome activation (ZGA). During this maternal-to-zygotic transition (MZT), which is conserved in metazoans, transient totipotency is induced and zygotic transcription is initiated to form the blueprint for future development. Recent technological advances have enhanced our understanding of MZT regulation, revealing common themes across species and leading to new fundamental insights about transcription, mRNA decay and translation.

An atlas of transcription initiation reveals regulatory principles of gene and transposable element expression in early mammalian development

Transcriptional activation of the embryonic genome (EGA) is a major developmental landmark enabling the embryo to become independent from maternal control. The magnitude and control of transcriptional reprogramming during this event across mammals remains poorly understood. Here, we developed Smart-seq+5' for high sensitivity, full-length transcript coverage and simultaneous capture of 5' transcript information from single cells and single embryos. Using Smart-seq+5', we profiled 34 developmental stages in 5 mammalian species and provide an extensive characterization of the transcriptional repertoire of early development before, during, and after EGA. We demonstrate widespread transposable element (TE)-driven transcription across species, including, remarkably, of DNA transposons. We identify 19,657 TE-driven genic transcripts, suggesting extensive TE co-option in early development over evolutionary timescales. TEs display similar expression dynamics across species and species-specific patterns, suggesting shared and divergent regulation. Our work provides a powerful resource for understanding transcriptional regulation of mammalian development.

Transcriptomic insights into developmental arrest in fluorescent labeling transgenic Asian elephant (Elephas maximus) embryos via inter-order cloning

Introduction

Asian elephants (Elephas maximus) provide a unique model for studying cloning in large mammals. As an endangered species with declining populations and limited oocyte availability, interspecies somatic cell nuclear transfer (iSCNT) combined with transcriptomic analysis holds promise for advancing iSCNT embryonic arrest development and further facilitating applications in conservation efforts, therapeutic cloning, and regenerative medicine.

Methods

This study conducted low-input RNA sequencing analyses on transgenic Asian elephant-pig (AE-P) inter-order cloned embryos expressing enhanced green fluorescent protein (EGFP) at the 2- and 4-cell stages. Differential gene expressions, pathway enrichment, and hub gene analyses were performed to identify the molecular mechanisms and core genes influencing normal and arrest development.

Results and Discussion

Approximately 25% of clean reads successfully aligned with the Asian elephant genome. The transcriptomic analysis revealed that inter-order cloned embryos with earlier cleavage at the 2- and 4-cell stages exhibited signs of residual transcriptomic memory and incomplete epigenetic reprogramming, while arrested embryos showed indications of nucleocytoplasmic incompatibility and nDNA-mtDNA mismatch. Hub gene analyses indicated core genes such as NDUFC2, NDUFS3, NDUFAB1, SDHC, SDHB, NUP54, NUP43, NUP37, NDC1, CDK1, and CCNB1 linked to energy production, nucleocytoplasmic transport, and cell cycle regulation highlighting the overall challenges in cloning Asian elephant inter-order embryos. Altogether, the analysis of high-throughput sequencing enhances the reliability of iSCNT production in this study, advancing our understanding of cellular reprogramming and molecular roadblocks in AE-P inter-order cloned embryos. Transcriptomic analyses have identified key factors contributing to developmental barriers in iSCNT, offering valuable insights into the complexities of these challenges.

Beta-hydroxybutyrate alters bovine preimplantation embryo development through transcriptional and epigenetic mechanisms†

Developing embryos are susceptible to fluctuations in the nutrients and metabolites concentrations within the reproductive tract, which can lead to alterations in their developmental trajectory. Ketotic dairy cows have diminished fertility, and elevated levels of the ketone body beta-hydroxybutyrate (BHB) have been associated with poor embryonic development. We used an in vitro model based on either in vitro fertilization (IVF) or parthenogenesis to investigate the effects of BHB on the preimplantation bovine embryo development, epigenome, and transcriptome. Embryo culture medium was supplemented with BHB at a similar concentration to that present in the blood of cows suffering with severe ketosis, followed by analysis of blastocysts formation rate, diameter, total number of cells, levels of H3K9 beta-hydroxybutyrylation (H3K9bhb), apoptosis, and transcriptional alterations. As a result, we observed that BHB reduced the blastocysts rates, the diameter and the total number of cells in both parthenotes and IVF embryos. Exposure to BHB for either 3 or 7 days greatly increased the H3K9bhb levels in parthenotes at the 8-cells and blastocyst stages, and affected the expression of HDAC1, TET1, DNMT1, KDM6B, NANOG, and MTHFD2 genes. Additionally, culture of IVF embryos with BHB for 7 days dramatically increased H3K9bhb and reduced NANOG in blastocysts. RNA-seq analysis of IVF blastocysts revealed that BHB modulated the expression of 118 genes, which were involved with biological processes such as embryonic development, implantation, reproduction, proliferation, and metabolism. These findings provided valuable insights into the mechanisms through which BHB disrupts preimplantation embryonic development and affects the fertility in dairy cows.

The role(s) of NF-Y in development and differentiation

NF-Y is a conserved sequence-specific trimeric Transcription Factor -TF- binding to the CCAAT element. We review here the role(s) in development, from pre-implantation embryo to terminally differentiated tissues, by rationalizing and commenting on genetic, genomic, epigenetic and biochemical studies. This effort brings to light the impact of NF-YA isoforms on stemness and differentiation, as well as binding to distal vs promoter proximal sites and connections with selected TFs.

METTL3 improves the development of somatic cell nuclear transfer embryos through AURKB and H3S10ph in goats

Developmental abnormalities are more common in somatic cell nuclear transfer (SCNT) embryos due to epigenetic barriers that occur during the maternal-to-zygotic transition (MZT). N6-methyladenosine (m6A) is an RNA epigenetic modification that plays a significant role in numerous biological processes. However, the relationship between m6A and SCNT embryonic development is largely unexplored. In the present study, we found that the low expression of m6A methyltransferase-like 3 (METTL3) was associated with developmental arrest before zygotic genome activation (ZGA) in goat SCNT embryos and that karyokinesis defects were evident during their development. Notably, we demonstrated that METTL3 overexpression rescued the karyokinesis abnormalities, enhanced embryonic development and elevated the blastocyst formation rate. Further experiments revealed that METTL3 could mitigate the defects of maternal mRNA degradation, enhance the translation of Aurora kinase B (AURKB) and increase the phosphorylation of serine 10 on histone H3 (H3S10ph) to ensure the normal karyokinesis in SCNT embryos before ZGA in goats. Overall, our study highlights the essential role of METTL3 in enhancing the development of goat SCNT embryos. These findings indicate that METTL3 is critical for optimal SCNT efficiency and advance our understanding of m6A's role in embryonic development.

Molecular cloning of PRD-like homeobox genes expressed in bovine oocytes and early IVF embryos

BACKGROUND

Embryonic genome activation (EGA) is a critical step in early embryonic development, as it marks the transition from relying on maternal factors to the initiation of transcription from embryo's own genome. The factors associated with EGA are not well understood and need further investigation. PRD-like (PRDL) homeodomain transcription factors (TFs) are considered to play crucial roles in this early event during development but these TFs have evolved differently, even within mammalian lineages. Different numbers of PRDL TFs have been predicted in bovine (Bos taurus); however, their divergent evolution requires species-specific confirmation and functional investigations.

RESULTS

In this study, we conducted molecular cloning of mRNAs for the PRDL TFs ARGFX, DUXA, LEUTX, NOBOX, TPRX1, TPRX2, and TPRX3 in bovine oocytes or in vitro fertilized (IVF) preimplantation embryos. Our results confirmed the expression of PRDL TF genes in early bovine development at the cDNA level and uncovered their structures. For each investigated PRDL TF gene, we isolated at least one homeodomain-encoding cDNA fragment, indicative of DNA binding and thus potential role in transcriptional regulation in developing bovine embryos. Additionally, our cDNA cloning approach allowed us to reveal breed-related differences in bovine, as evidenced by the identification of a high number of single nucleotide variants (SNVs) across the PRDL class homeobox genes. Subsequently, we observed the prediction of the 9aa transactivation domain (9aaTAD) motif in the putative protein sequence of TPRX3 leading us to conduct functional analysis of this gene. We demonstrated that the TPRX3 overexpression in bovine fibroblast induces not only protein-coding genes but also short noncoding RNAs involved in splicing and RNA editing. We supported this finding by identifying a shared set of genes between our and published bovine early embryo development datasets.

CONCLUSIONS

Providing full-length cDNA evidence for previously predicted homeobox genes that belong to PRDL class improves the annotation of the bovine genome. Updating the annotation with seven developmentally-important genes will enhance the accuracy of RNAseq analysis with datasets derived from bovine preimplantation embryos. In addition, the absence of TPRX3 in humans highlights the species-specific and TF-specific regulation of biological processes during early embryo development.

Characterization of bovine long non-coding RNAs, OOSNCR1, OOSNCR2 and OOSNCR3, and their roles in oocyte maturation and early embryonic development

In mammals, early embryogenesis relies heavily on the regulation of maternal transcripts including protein-coding and non-coding RNAs stored in oocytes. In this study, the expression of three bovine oocyte expressed long non-coding RNAs (lncRNAs), OOSNCR1, OOSNCR2, and OOSNCR3, was characterized in somatic tissues, the ovarian follicle, and throughout early embryonic development. Moreover, the functional requirement of each transcript during oocyte maturation and early embryonic development was investigated using a siRNA-mediated knockdown approach. Tissue distribution analysis revealed that OOSNCR1, OOSNCR2 and OOSNCR3 are predominantly expressed in fetal ovaries. Follicular cell expression analysis revealed that these lncRNAs are highly expressed in the oocytes, with minor expression detected in the cumulus cells (CCs) and mural granulosa cells (mGCs). The expression for all three genes was highest during oocyte maturation, decreased at fertilization, and ceased altogether by the 16-cell stage. Knockdown of OOSNCR1, OOSNCR2 and OOSNCR3 in immature oocytes was achieved by microinjection of the cumulus-enclosed germinal vesicle (GV) oocytes with siRNAs targeting these lncRNAs. Knockdown of OOSNCR1, OOSNCR2 and OOSNCR3 did not affect cumulus expansion, but oocyte survival at 12 h post-insemination was significantly reduced. In addition, knockdown of OOSNCR1, OOSNCR2 and OOSNCR3 in immature oocytes resulted in a decreased rate of blastocyst development, and reduced expression of genes associated with oocyte competency such as nucleoplasmin 2 (NPM2), growth differentiation factor 9 (GDF9), bone morphogenetic protein 15 (BMP15), and JY-1 in MII oocytes. The data herein suggest a functional requirement of OOSNCR1, OOSNCR2, and OOSNCR3 during bovine oocyte maturation and early embryogenesis.

Transcription of Endogenous Retroviruses: Broad and Precise Mechanisms of Control

Endogenous retroviruses (ERVs) are the remnants of retroviral germline infections and are highly abundant in the genomes of vertebrates. At one time considered to be nothing more than inert 'junk' within genomes, ERVs have been tolerated within host genomes over vast timescales, and their study continues to reveal complex co-evolutionary histories within their respective host species. For example, multiple instances have been characterized of ERVs having been 'borrowed' for normal physiology, from single copies to ones involved in various regulatory networks such as innate immunity and during early development. Within the cell, the accessibility of ERVs is normally tightly controlled by epigenetic mechanisms such as DNA methylation or histone modifications. However, these silencing mechanisms of ERVs are reversible, and epigenetic alterations to the chromatin landscape can thus lead to their aberrant expression, as is observed in abnormal cellular environments such as in tumors. In this review, we focus on ERV transcriptional control and draw parallels and distinctions concerning the loss of regulation in disease, as well as their precise regulation in early development.

Identification of novel cattle (Bos taurus) genes and biological insights of their function in pre-implantation embryo development

BACKGROUND

Appropriate regulation of genes expressed in oocytes and embryos is essential for acquisition of developmental competence in mammals. Here, we hypothesized that several genes expressed in oocytes and pre-implantation embryos remain unknown. Our goal was to reconstruct the transcriptome of oocytes (germinal vesicle and metaphase II) and pre-implantation cattle embryos (blastocysts) using short-read and long-read sequences to identify putative new genes.

RESULTS

We identified 274,342 transcript sequences and 3,033 of those loci do not match a gene present in official annotations and thus are potential new genes. Notably, 63.67% (1,931/3,033) of potential novel genes exhibited coding potential. Also noteworthy, 97.92% of the putative novel genes overlapped annotation with transposable elements. Comparative analysis of transcript abundance identified that 1,840 novel genes (recently added to the annotation) or potential new genes were differentially expressed between developmental stages (FDR < 0.01). We also determined that 522 novel or potential new genes (448 and 34, respectively) were upregulated at eight-cell embryos compared to oocytes (FDR < 0.01). In eight-cell embryos, 102 novel or putative new genes were co-expressed (|r|> 0.85, P < 1 × 10-8) with several genes annotated with gene ontology biological processes related to pluripotency maintenance and embryo development. CRISPR-Cas9 genome editing confirmed that the disruption of one of the novel genes highly expressed in eight-cell embryos reduced blastocyst development (ENSBTAG00000068261, P = 1.55 × 10-7).

CONCLUSIONS

Our results revealed several putative new genes that need careful annotation. Many of the putative new genes have dynamic regulation during pre-implantation development and are important components of gene regulatory networks involved in pluripotency and blastocyst formation.

The developmental and evolutionary characteristics of transcription factor binding site clustered regions based on an explainable machine learning model

Gene expression is temporally and spatially regulated by the interaction of transcription factors (TFs) and cis-regulatory elements (CREs). The uneven distribution of TF binding sites across the genome poses challenges in understanding how this distribution evolves to regulate spatio-temporal gene expression and consequent heritable phenotypic variation. In this study, chromatin accessibility profiles and gene expression profiles were collected from several species including mammals (human, mouse, bovine), fish (zebrafish and medaka), and chicken. Transcription factor binding sites clustered regions (TFCRs) at different embryonic stages were characterized to investigate regulatory evolution. The study revealed dynamic changes in TFCR distribution during embryonic development and species evolution. The synchronization between TFCR complexity and gene expression was assessed across species using RegulatoryScore. Additionally, an explainable machine learning model highlighted the importance of the distance between TFCR and promoter in the coordinated regulation of TFCRs on gene expression. Our results revealed the developmental and evolutionary dynamics of TFCRs during embryonic development from fish, chicken to mammals. These data provide valuable resources for exploring the relationship between transcriptional regulation and phenotypic differences during embryonic development.

Maternal nutrition and fetal imprinting of the male progeny

The global population as well as the demand for human food is rapidly growing worldwide, which necessitates improvement of efficiency in livestock operations. In this context, environmental factors during fetal and/or neonatal life have been observed to influence normal physical and physiological function of an individual during adulthood, and this phenomenon is called fetal or developmental programming. While numerous studies have reported the impact of maternal factors on development of the female progeny, limited information is available on the potential effects of fetal programming on reproductive function of the male offspring. Therefore, the objective for this review article was to focus on available literature regarding the impact of maternal factors, particularly maternal nutrition, on reproductive system of the male offspring. To this end, we highlighted developmental programming of the male offspring in domestic species (i.e., pig, cow and sheep) as well as laboratory species (i.e., mice and rat) during pregnancy and lactation. In this sense, we pointed out the effects of maternal nutrition on various functions of the male offspring including hypothalamic-pituitary axis, hormonal levels, testicular tissue and semen parameters.

Efficiency of embryo complementation and pluripotency maintenance following multiple passaging of in vitro-derived bovine embryos

Context Current methods to obtain bovine embryos of high genetic merit include approaches that require skilled techniques for low-efficiency cloning strategies. Aims The overall goal herein was to identify the efficacy of alternative methods for producing multiple embryos through blastomere complementation while determining maintenance of cell pluripotency. Methods Bovine oocytes were fertilised in vitro to produce 4-cell embryos from which blastomeres were isolated and cultured as 2-cell aggregates using a well-of-the-well system. Aggregates were returned to incubation up to 7days (Passage 1). A second passage of complement embryos was achieved by splitting 4-cell Passage 1 embryos. Passaged embryos reaching the blastocyst stage were characterised for cell number and cell lineage specification in replicate with non-reconstructed zona-intact embryos. Key results Passage 1 and 2 embryo complements yielded 29% and 25% blastocyst development, respectively. Passage 1 embryos formed blastocysts, but with a reduction in expression of SOX2 and decreased size compared to non-reconstructed zona-intact embryos. Passage 2 embryos had a complete lack of SOX2 expression and a reduction in transcript abundance of SOX2 and SOX17, suggesting loss of pluripotency markers that primarily affected inner cell mass (ICM) and hypoblast formation. Conclusions In vitro fertilised bovine embryos can be reconstructed with multiple passaging to generate genetically identical embryos. Increased passaging drives trophectoderm cell lineage specification while compromising ICM formation. Implications These results may provide an alternative strategy for producing genetically identical bovine embryos through blastomere complementation with applications towards the development of trophoblast and placental models of early development.

Single-cell transcriptional landscapes of bovine peri-implantation development

Supporting healthy pregnancy outcomes requires a comprehensive understanding of the molecular and cellular programs of peri-implantation development, when most pregnancy failure occurs. Here, we present single-cell transcriptomes of bovine peri-implantation embryo development at day 12, 14, 16, and 18 post-fertilization. We defined the cellular composition and gene expression of embryonic disc, hypoblast, and trophoblast lineages in bovine peri-implantation embryos, and identified markers and pathway signaling that represent distinct stages of bovine peri-implantation lineages; the expression of selected markers was validated in peri-implantation embryos. Using detailed time-course transcriptomic analyses, we revealed a previously unrecognized primitive trophoblast cell lineage. We also characterized conserved and divergence peri-implantation lineage programs between bovine and other mammalian species. Finally, we established cell-cell communication signaling underlies embryonic and extraembryonic cell interaction to ensure proper early development. These data provide foundational information to discover essential biological signaling underpinning bovine peri-implantation development.

Jump-starting life: balancing transposable element co-option and genome integrity in the developing mammalian embryo

Remnants of transposable elements (TEs) are widely expressed throughout mammalian embryo development. Originally infesting our genomes as selfish elements and acting as a source of genome instability, several of these elements have been co-opted as part of a complex system of genome regulation. Many TEs have lost transposition ability and their transcriptional potential has been tampered as a result of interactions with the host throughout evolutionary time. It has been proposed that TEs have been ultimately repurposed to function as gene regulatory hubs scattered throughout our genomes. In the early embryo in particular, TEs find a perfect environment of naïve chromatin to escape transcriptional repression by the host. As a consequence, it is thought that hosts found ways to co-opt TE sequences to regulate large-scale changes in chromatin and transcription state of their genomes. In this review, we discuss several examples of TEs expressed during embryo development, their potential for co-option in genome regulation and the evolutionary pressures on TEs and on our genomes.

Expression and function of NF-Y subunits in cancer

NF-Y is a Transcription Factor (TF) targeting the CCAAT box regulatory element. It consists of the NF-YB/NF-YC heterodimer, each containing an Histone Fold Domain (HFD), and the sequence-specific subunit NF-YA. NF-YA expression is associated with cell proliferation and absent in some post-mitotic cells. The review summarizes recent findings impacting on cancer development. The logic of the NF-Y regulome points to pro-growth, oncogenic genes in the cell-cycle, metabolism and transcriptional regulation routes. NF-YA is involved in growth/differentiation decisions upon cell-cycle re-entry after mitosis and it is widely overexpressed in tumors, the HFD subunits in some tumor types or subtypes. Overexpression of NF-Y -mostly NF-YA- is oncogenic and decreases sensitivity to anti-neoplastic drugs. The specific roles of NF-YA and NF-YC isoforms generated by alternative splicing -AS- are discussed, including the prognostic value of their levels, although the specific molecular mechanisms of activity are still to be deciphered.

Species-Specific Transcription Factors Associated with Long Terminal Repeat Promoters of Endogenous Retroviruses: A Comprehensive Review

Endogenous retroviruses (ERVs) became a part of the eukaryotic genome through endogenization millions of years ago. Moreover, they have lost their innate capability of virulence or replication. Nevertheless, in eukaryotic cells, they actively engage in various activities that may be advantageous or disadvantageous to the cells. The mechanisms by which transcription is triggered and implicated in cellular processes are complex. Owing to the diversity in the expression of transcription factors (TFs) in cells and the TF-binding motifs of viruses, the comprehensibility of ERV initiation and its impact on cellular functions are unclear. Currently, several factors are known to be related to their initiation. TFs that bind to the viral long-terminal repeat (LTR) are critical initiators. This review discusses the TFs shown to actively associate with ERV stimulation across species such as humans, mice, pigs, monkeys, zebrafish, Drosophila, and yeast. A comprehensive summary of the expression of previously reported TFs may aid in identifying similarities between animal species and endogenous viruses. Moreover, an in-depth understanding of ERV expression will assist in elucidating their physiological roles in eukaryotic cell development and in clarifying their relationship with endogenous retrovirus-associated diseases.

Evaluating histone modification analysis of individual preimplantation embryos

BACKGROUND

We previously reported a modification of the CUT&Tag method (NTU-CAT) that allows genome-wide histone modification analysis in individual preimplantation embryos. In the present study, NTU-CAT was further simplified by taking advantage of the Well-of-the-Well (WOW) system, which enables the processing of multiple embryos in a shorter time with less reagent and cell loss during the procedure (WOW-CUT&Tag, WOW-CAT).

RESULTS

WOW-CAT allowed histone modification profiling from not only a single blastocyst but also from a portion of it. WOW-CAT generated similar H3K4me3 profiles as NTU-CAT, but they were closer to the profiles produced by chromatin immunoprecipitation-sequencing, such as a valley-like trend and relatively lower false positive rates, indicating that WOW-CAT may attenuate the bias of Tn5 transposase to cut open chromatin regions. Simultaneous WOW-CAT of two halves of single blastocysts was conducted to analyze two different histone modifications (H3K4me3 and H3K27ac) within the same embryo. Furthermore, trophectoderm cells were biopsied and subjected to WOW-CAT in anticipation of preimplantation diagnosis of histone modifications. WOW-CAT allowed the monitoring of epigenetic modifications in the main body of the embryo. For example, analysis of H3K4me3 modifications of XIST and DDX3Y in trophectoderm biopsies could be used to sex embryos in combination with quantitative PCR, but without the need for deep sequencing.

CONCLUSIONS

These results suggest the applicability of WOW-CAT for flexible epigenetic analysis of individual embryos in preimplantation epigenetic diagnosis.

Cyclophosphamide reduces gene transcriptional activity and embryo in vitro development by inhibiting NF-κB expression through decreasing AcH4K12

Cyclophosphamide (CTX), a widely used chemotherapeutic agent for cancer treatment, has been associated with long-term toxicity and detrimental effects on oocytes and ovaries, resulting in female reproductive dysfunction. This study aimed to investigate the potential impact of CTX on in vitro maturation (IVM) injury of porcine oocytes and subsequent embryonic development, as well as its effects on epigenetic modification and gene activation during early embryonic development. The results demonstrated that CTX treatment caused aberrant spindle structure and mitochondrial dysfunction during oocyte maturation, inducing DNA damage and early apoptosis, which consequently disrupted meiotic maturation. Indeed, CTX significantly reduced the in vitro developmental capacity of porcine embryos, and induced DNA damage and apoptosis in in vitro fertilization (IVF) blastocysts. Importantly, CTX induced abnormal histone modification of AcH4K12 in early porcine embryos. Moreover, addition of LBH589 before zygotic genome activation (ZGA) effectively increased AcH4K12 levels and restored the protein expression of NF-κB, which can effectively enhance the in vitro developmental potential of IVF embryos. The DNA damage and apoptosis induced by CTX compromised the quality of the blastocysts, which were recovered by supplementation with LBH589. This restoration was accompanied by down-regulation of BAX mRNA expression and up-regulation of BCL2, POU5F1, SOX2 and SOD1 mRNA expression. These findings indicated that CTX caused abnormal histone modification of AcH4K12 in early porcine embryos and reduced the protein expression of NF-κB, a key regulator of early embryo development, which may block subsequent ZGA processes.

Pleomorphic adenoma gene1 in reproduction and implication for embryonic survival in cattle: a review

The pleomorphic adenoma gene1 (PLAG1) encodes a DNA-binding, C2H2 zinc-finger protein which acts as a transcription factor that regulates the expression of diverse genes across different organs and tissues; hence, the name pleomorphic. Rearrangements of the PLAG1 gene, and/or overexpression, are associated with benign tumors and cancers in a variety of tissues. This is best described for pleomorphic adenoma of the salivary glands in humans. The most notable expression of PLAG1 occurs during embryonic and fetal development, with lesser expression after birth. Evidence has accumulated of a role for PLAG1 protein in normal early embryonic development and placentation in mammals. PLAG1 protein influences the expression of the ike growth factor 2 (IGF2) gene and production of IGF2 protein. IGF2 is an important mitogen in ovarian follicles/oocytes, embryos, and fetuses. The PLAG1-IGF2 axis, therefore, provides one pathway whereby PLAG1 protein can influence embryonic survival and pregnancy. PLAG1 also influences over 1,000 other genes in embryos including those associated with ribosomal assembly and proteins. Brahman (Bos indicus) heifers homozygous for the PLAG1 variant, rs109815800 (G > T), show greater fertility than contemporary heifers with either one, or no copy, of the variant. Greater fertility in heifers homozygous for rs109815800 could be the result of early puberty and/or greater embryonic survival. The present review first looks at the broader roles of the PLAG1 gene and PLAG1 protein and then focuses on the emerging role of PLAG1/PLAG1 in embryonic development and pregnancy. A deeper understanding of factors which influence embryonic development is required for the next transformational increase in embryonic survival and successful pregnancy for both in vivo and in vitro derived embryos in cattle.

Single-cell RNA-seq and single-cell bisulfite-sequencing reveal insights into yak preimplantation embryogenesis

Extensive epigenetic reprogramming occurs during preimplantation embryonic development. However, the impact of DNA methylation in plateau yak preimplantation embryos and how epigenetic reprogramming contributes to transcriptional regulatory networks are unclear. In this study, we quantified gene expression and DNA methylation in oocytes and a series of yak embryos at different developmental stages and at single-cell resolution using single-cell bisulfite-sequencing and RNA-seq. We characterized embryonic genome activation and maternal transcript degradation and mapped epigenetic reprogramming events critical for embryonic development. Through cross-species transcriptome analysis, we identified 31 conserved maternal hub genes and 39 conserved zygotic hub genes, including SIN3A, PRC1, HDAC1/2, and HSPD1. Notably, by combining single-cell DNA methylation and transcriptome analysis, we identified 43 candidate methylation driver genes, such as AURKA, NUSAP1, CENPF, and PLK1, that may be associated with embryonic development. Finally, using functional approaches, we further determined that the epigenetic modifications associated with the histone deacetylases HDAC1/2 are essential for embryonic development and that the deubiquitinating enzyme USP7 may affect embryonic development by regulating DNA methylation. Our data represent an extensive resource on the transcriptional dynamics of yak embryonic development and DNA methylation remodeling, and provide new insights into strategies for the conservation of germplasm resources, as well as a better understanding of mammalian early embryonic development that can be applied to investigate the causes of early developmental disorders.

PDGFRβ Activation Induced the Bovine Embryonic Genome Activation via Enhanced NFYA Nuclear Localization

Embryonic genome activation (EGA) is a critical step during embryonic development. Several transcription factors have been identified that play major roles in initiating EGA; however, this gradual and complex mechanism still needs to be explored. In this study, we investigated the role of nuclear transcription factor Y subunit A (NFYA) in bovine EGA and bovine embryonic development and its relationship with the platelet-derived growth factor receptor-β (PDGFRβ) by using a potent selective activator (PDGF-BB) and inhibitor (CP-673451) of PDGF receptors. Activation and inhibition of PDGFRβ using PDGF-BB and CP-673451 revealed that NFYA expression is significantly (p < 0.05) affected by the PDGFRβ. In addition, PDGFRβ mRNA expression was significantly increased (p < 0.05) in the activator group and significantly decreased (p < 0.05) in the inhibitor group when compared with PDGFRα. Downregulation of NFYA following PDGFRβ inhibition was associated with the expression of critical EGA-related genes, bovine embryo development rate, and implantation potential. Moreover, ROS and mitochondrial apoptosis levels and expression of pluripotency-related markers necessary for inner cell mass development were also significantly (p < 0.05) affected by the downregulation of NFYA while interrupting trophoblast cell (CDX2) differentiation. In conclusion, the PDGFRβ-NFYA axis is critical for bovine embryonic genome activation and embryonic development.

Molecular and cellular programs underlying the development of bovine pre-implantation embryos

Early embryonic mortality is a major cause of infertility in cattle, yet the underlying molecular causes remain a mystery. Over the past half century, assisted reproductive technologies such as in vitro fertilisation and somatic cell nuclear transfer have been used to improve cattle reproductive efficiency; however, reduced embryo developmental potential is seen compared to their in vivo counterparts. Recent years have seen exciting progress across bovine embryo research, including genomic profiling of embryogenesis, new methods for improving embryo competence, and experimenting on building bovine embryos from stem cell cultures. These advances are beginning to define bovine embryo molecular and cellular programs and could potentially lead to improved embryo health. Here, I highlight the current status of molecular determinants and cellular programs of bovine embryo development and new opportunities to improve the bovine embryo health.

An organism-wide ATAC-seq peak catalog for the bovine and its use to identify regulatory variants

We report the generation of an organism-wide catalog of 976,813 cis-acting regulatory elements for the bovine detected by the assay for transposase accessible chromatin using sequencing (ATAC-seq). We regroup these regulatory elements in 16 components by nonnegative matrix factorization. Correlation between the genome-wide density of peaks and transcription start sites, correlation between peak accessibility and expression of neighboring genes, and enrichment in transcription factor binding motifs support their regulatory potential. Using a previously established catalog of 12,736,643 variants, we show that the proportion of single-nucleotide polymorphisms mapping to ATAC-seq peaks is higher than expected and that this is owing to an approximately 1.3-fold higher mutation rate within peaks. Their site frequency spectrum indicates that variants in ATAC-seq peaks are subject to purifying selection. We generate eQTL data sets for liver and blood and show that variants that drive eQTL fall into liver- and blood-specific ATAC-seq peaks more often than expected by chance. We combine ATAC-seq and eQTL data to estimate that the proportion of regulatory variants mapping to ATAC-seq peaks is approximately one in three and that the proportion of variants mapping to ATAC-seq peaks that are regulatory is approximately one in 25. We discuss the implication of these findings on the utility of ATAC-seq information to improve the accuracy of genomic selection.

Chromatin accessibility memory of donor cells disrupts bovine somatic cell nuclear transfer blastocysts development

The post-transfer developmental capacity of bovine somatic cell nuclear transfer (SCNT) blastocysts is reduced, implying that abnormalities in gene expression regulation are present at blastocyst stage. Chromatin accessibility, as an indicator for transcriptional regulatory elements mediating gene transcription activity, has heretofore been largely unexplored in SCNT embryos, especially at blastocyst stage. In the present study, single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) of in vivo and SCNT blastocysts were conducted to segregate lineages and demonstrate the aberrant chromatin accessibility of transcription factors (TFs) related to inner cell mass (ICM) development in SCNT blastocysts. Pseudotime analysis of lineage segregation further reflected dysregulated chromatin accessibility dynamics of TFs in the ICM of SCNT blastocysts compared to their in vivo counterparts. ATAC- and ChIP-seq results of SCNT donor cells revealed that the aberrant chromatin accessibility in the ICM of SCNT blastocysts was due to the persistence of chromatin accessibility memory at corresponding loci in the donor cells, with strong enrichment of trimethylation of histone H3 at lysine 4 (H3K4me3) at these loci. Correction of the aberrant chromatin accessibility through demethylation of H3K4me3 by KDM5B diminished the expression of related genes (e.g., BCL11B) and significantly improved the ICM proliferation in SCNT blastocysts. This effect was confirmed by knocking down BCL11B in SCNT embryos to down-regulate p21 and alleviate the inhibition of ICM proliferation. These findings expand our understanding of the chromatin accessibility abnormalities in SCNT blastocysts and BCL11B may be a potential target to improve SCNT efficiency.

β-hydroxybutyrate impairs bovine oocyte maturation via pyruvate dehydrogenase (PDH) associated energy metabolism abnormality

Background: Ketosis is one of the most frequent and costly metabolic disorders in high-producing dairy cows, and negatively associated with the health and reproductive performance of bovine. Ketosis is mainly caused by the accumulation of ketone body β-hydroxybutyric acid and its diagnosis is based on β-hydroxybutyrate (βHB) concentration in blood. Methods: In this study, we investigated the effects of βHB on bovine oocyte maturation in the concentration of subclinical (1.2 mM) βHB and clinical (3.6 mM). Results: The results showed βHB disrupted bovine oocyte maturation and development capacity. Further analysis showed that βHB induced oxidative stress and mitochondrial dysfunction, as indicated by the increased level of reactive oxygen species (ROS), disrupted mitochondrial structure and distribution, and depolarized membrane potential. Furthermore, oxidative stress triggered early apoptosis, as shown by the enhanced levels of Caspase-3 and Annexin-V. Moreover, 3.6 mM βHB induced the disruption of the pyruvate dehydrogenase (PDH) activity, showing with the decrease of the global acetylation modification and the increase of the abnormal spindle rate. Conclusion: Our study showed that βHB in subclinical/clinical concentration had toxic effects on mitochondrial function and PDH activity, which might affect energy metabolism and epigenetic modification of bovine oocytes and embryos.

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.

Single-cell transcriptional landscapes of bovine peri-implantation development

Supporting healthy pregnancy outcomes requires a comprehensive understanding of the cellular hierarchy and underlying molecular mechanisms during peri-implantation development. Here, we present a single-cell transcriptome-wide view of the bovine peri-implantation embryo development at day 12, 14, 16 and 18, when most of the pregnancy failure occurs in cattle. We defined the development and dynamic progression of cellular composition and gene expression of embryonic disc, hypoblast, and trophoblast lineages during bovine peri-implantation development. Notably, the comprehensive transcriptomic mapping of trophoblast development revealed a previously unrecognized primitive trophoblast cell lineage that is responsible for pregnancy maintenance in bovine prior to the time when binucleate cells emerge. We analyzed novel markers for the cell lineage development during bovine early development. We also identified cell-cell communication signaling underling embryonic and extraembryonic cell interaction to ensure proper early development. Collectively, our work provides foundational information to discover essential biological pathways underpinning bovine peri-implantation development and the molecular causes of the early pregnancy failure during this critical period.

SIGNIFICANCE STATEMENT

Peri-implantation development is essential for successful reproduction in mammalian species, and cattle have a unique process of elongation that proceeds for two weeks prior to implantation and represents a period when many pregnancies fail. Although the bovine embryo elongation has been studied histologically, the essential cellular and molecular factors governing lineage differentiation remain unexplored. This study profiled the transcriptome of single cells in the bovine peri-implantation development throughout day 12, 14, 16, and 18, and identified peri-implantation stage-related features of cell lineages. The candidate regulatory genes, factors, pathways and embryonic and extraembryonic cell interactions were also prioritized to ensure proper embryo elongation in cattle.

Chromatin Reprogramming of In Vitro Fertilized and Somatic Cell Nuclear Transfer Bovine Embryos During Embryonic Genome Activation

Reprogramming of the gamete into a developmentally competent embryo identity is a fundamental aspect of preimplantation development. One of the most important processes of this reprogramming is the transcriptional awakening during embryonic genome activation (EGA), which robustly occurs in fertilized embryos but is defective in most somatic cell nuclear transfer (SCNT) embryos. However, little is known about the genome-wide underlying chromatin landscape during EGA in SCNT embryos and how it differs from a fertilized embryo. By profiling open chromatin genome-wide in both types of bovine embryos, we find that SCNT embryos fail to reprogram a subset of the EGA gene targets that are normally activated in fertilized embryos. Importantly, a small number of transcription factor (TF) motifs explain most chromatin regions that fail to open in SCNT embryos suggesting that over-expression of a limited number of TFs may provide more robust reprogramming. One such TF, the zygotically-expressed bovine gene DUXC which is a homologue of EGA factors DUX/DUX4 in mouse/human, is alone capable of activating ∼84% of all EGA transcripts that fail to activate normally in SCNT embryos. Additionally, single-cell chromatin profiling revealed low intra-embryo heterogeneity but high inter-embryo heterogeneity in SCNT embryos and an uncoupling of cell division and open chromatin reprogramming during EGA. Surprisingly, our data also indicate that transcriptional defects may arise downstream of promoter chromatin opening in SCNT embryos, suggesting additional mechanistic insights into how and why transcription at EGA is dysregulated. We anticipate that our work will lead to altered SCNT protocols to increase the developmental competency of bovine SCNT embryos.

Histone remodeling reflects conserved mechanisms of bovine and human preimplantation development

How histone modifications regulate changes in gene expression during preimplantation development in any species remains poorly understood. Using CUT&Tag to overcome limiting amounts of biological material, we profiled two activating (H3K4me3 and H3K27ac) and two repressive (H3K9me3 and H3K27me3) marks in bovine oocytes, 2-, 4-, and 8-cell embryos, morula, blastocysts, inner cell mass, and trophectoderm. In oocytes, broad bivalent domains mark developmental genes, and prior to embryonic genome activation (EGA), H3K9me3 and H3K27me3 co-occupy gene bodies, suggesting a global mechanism for transcription repression. During EGA, chromatin accessibility is established before canonical H3K4me3 and H3K27ac signatures. Embryonic transcription is required for this remodeling, indicating that maternally provided products alone are insufficient for reprogramming. Last, H3K27me3 plays a major role in restriction of cellular potency, as blastocyst lineages are defined by differential polycomb repression and transcription factor activity. Notably, inferred regulators of EGA and blastocyst formation strongly resemble those described in humans, as opposed to mice. These similarities suggest that cattle are a better model than rodents to investigate the molecular basis of human preimplantation development.

Regulation, functions and transmission of bivalent chromatin during mammalian development

Cells differentiate and progress through development guided by a dynamic chromatin landscape that mediates gene expression programmes. During development, mammalian cells display a paradoxical chromatin state: histone modifications associated with gene activation (trimethylated histone H3 Lys4 (H3K4me3)) and with gene repression (trimethylated H3 Lys27 (H3K27me3)) co-occur at promoters of developmental genes. This bivalent chromatin modification state is thought to poise important regulatory genes for expression or repression during cell-lineage specification. In this Review, we discuss recent work that has expanded our understanding of the molecular basis of bivalent chromatin and its contributions to mammalian development. We describe the factors that establish bivalency, especially histone-lysine N-methyltransferase 2B (KMT2B) and Polycomb repressive complex 2 (PRC2), and consider evidence indicating that PRC1 shapes bivalency and may contribute to its transmission between generations. We posit that bivalency is a key feature of germline and embryonic stem cells, as well as other types of stem and progenitor cells. Finally, we discuss the relevance of bivalent chromtin to human development and cancer, and outline avenues of future research.

Gene Expression of Aquaporins (AQPs) in Cumulus Oocytes Complex and Embryo of Cattle

Aquaporins (AQPs) are proteins with various functions related to proper cell function and early development in mammals. The aim of this study was to evaluate the presence of AQPs and determine their mRNA levels in the cumulus oocyte complex (COC) of four bovine breeds and in blastocysts of five bovine crosses. Grade I, II and III COCs were collected by ovum pick up from non-lactating heifers of the Brahaman, Holstein, Gir and Romosinuano breeds. Embryos were produced in vitro up to the blastocyst stage of the bovine ♀Gir × ♂Holstein, ♀Holstein × ♂Gir, ♀Brahman × ♂Holstein, ♀Holstein × ♂Brahman, and ♀Romosinuano × ♂Holstein crosses. mRNA expression of AQP1-AQP12b was estimated in COC and embryos by real-time-PCR. The presence of the twelve AQPs in the COCs and bovine embryos was established. Additionally, significant differences were determined in the expression of AQP6 and AQP12b in COCs, as well as in transcripts levels of AQP4, AQP8 and AQP9 from bovine embryos. Gene expression of AQPs in COCs and bovine embryos is consistent with the previously described biological functions. This is the first report of AQPs in COC of Gir, Brahman, Holstein and Romosinuano and embryos of five crossbreeds between Bos indicus and B. taurus.

Dynamic Change of R-Loop Implicates in the Regulation of Zygotic Genome Activation in Mouse

In mice, zygotic genome activation (ZGA) occurs in two steps: minor ZGA at the one-cell stage and major ZGA at the two-cell stage. Regarding the regulation of gene transcription, minor ZGA is known to have unique features, including a transcriptionally permissive state of chromatin and insufficient splicing processes. The molecular characteristics may originate from extremely open chromatin states in the one-cell stage zygotes, yet the precise underlying mechanism has not been well studied. Recently, the R-loop, a triple-stranded nucleic acid structure of the DNA/RNA hybrid, has been implicated in gene transcription and DNA replication. Therefore, in the present study, we examined the changes in R-loop dynamics during mouse zygotic development, and its roles in zygotic transcription or DNA replication. Our analysis revealed that R-loops persist in the genome of metaphase II oocytes and preimplantation embryos from the zygote to the blastocyst stage. In particular, zygotic R-loop levels dynamically change as development proceeds, showing that R-loop levels decrease as pronucleus maturation occurs. Mechanistically, R-loop dynamics are likely linked to ZGA, as inhibition of either DNA replication or transcription at the time of minor ZGA decreases R-loop levels in the pronuclei of zygotes. However, the induction of DNA damage by treatment with anticancer agents, including cisplatin or doxorubicin, does not elicit genome-wide changes in zygotic R-loop levels. Therefore, our study suggests that R-loop formation is mechanistically associated with the regulation of mouse ZGA, especially minor ZGA, by modulating gene transcription and DNA replication.

The single-cell expression profile of transposable elements and transcription factors in human early biparental and uniparental embryonic development

Transposable elements (TEs) and transcription factors (TFs) are involved in the precise regulation of gene expression during the preimplantation stage. Activation of TEs is a key event for mammalian embryonic genome activation and preimplantation early embryonic development. TFs are involved in the regulation of drastic changes in gene expression patterns, but an inventory of the interplay between TEs and TFs during normal/abnormal human embryonic development is still lacking. Here we used single-cell RNA sequencing data generated from biparental and uniparental embryos to perform an integrative analysis of TE and TF expression. Our results showed that endogenous retroviruses (ERVs) are mainly expressed during the minor embryonic genome activation (EGA) process of early embryos, while Alu is gradually expressed in the middle and later stages. Some important ERVs (e.g., LTR5_Hs, MLT2A1) and Alu TEs are expressed at significantly lower levels in androgenic embryos. Integrative analysis revealed that the expression of the transcription factors CTCF and POU5F1 is correlated with the differential expression of ERV TEs. Comparative coexpression network analysis further showed distinct expression levels of important TFs (e.g., LEUTX and ZSCAN5A) in dizygotic embryos vs. parthenogenetic and androgenic embryos. This systematic investigation of TE and TF expression in human early embryonic development by single-cell RNA sequencing provides valuable insights into mammalian embryonic development.

Genome-wide chromatin accessibility analysis unveils open chromatin convergent evolution during polyploidization in cotton

Allopolyploidization, resulting in divergent genomes in the same cell, is believed to trigger a “genome shock”, leading to broad genetic and epigenetic changes. However, little is understood about chromatin and gene-expression dynamics as underlying driving forces during allopolyploidization. Here, we examined the genome-wide DNase I-hypersensitive site (DHS) and its variations in domesticated allotetraploid cotton ( Gossypium hirsutum and Gossypium barbadense , AADD) and its extant AA ( Gossypium arboreum ) and DD ( Gossypium raimondii ) progenitors. We observed distinct DHS distributions between G. arboreum and G. raimondii . In contrast, the DHSs of the two subgenomes of G. hirsutum and G. barbadense showed a convergent distribution. This convergent distribution of DHS was also present in the wild allotetraploids Gossypium darwinii and G. hirsutum var. yucatanense , but absent from a resynthesized hybrid of G. arboreum and G. raimondii , suggesting that it may be a common feature in polyploids, and not a consequence of domestication after polyploidization. We revealed that putative cis -regulatory elements (CREs) derived from polyploidization-related DHSs were dominated by several families, including Dof, ERF48, and BPC1. Strikingly, 56.6% of polyploidization-related DHSs were derived from transposable elements (TEs). Moreover, we observed positive correlations between DHS accessibility and the histone marks H3K4me3, H3K27me3, H3K36me3, H3K27ac, and H3K9ac, indicating that coordinated interplay among histone modifications, TEs, and CREs drives the DHS landscape dynamics under polyploidization. Collectively, these findings advance our understanding of the regulatory architecture in plants and underscore the complexity of regulome evolution during polyploidization.

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.

CRISPR/Cas gene editing in the human germline

The ease and efficacy of CRISPR/Cas9 germline gene editing in animal models paved the way to human germline gene editing (HGGE), by which permanent changes can be introduced into the embryo. Distinct genes can be knocked out to examine their function during embryonic development. Alternatively, specific sequences can be introduced which can be applied to correct disease-causing mutations. To date, it has been shown that the success of HGGE is dependent on various experimental parameters and that various hurdles (i.e. loss-of-heterozygosity and mosaicism) need to be overcome before clinical applications should be considered. Due to the shortage of human germline material and the ethical constraints concerning HGGE, alternative models such as stem cells have been evaluated as well, in terms of their predictive value on the genetic outcome for HGGE approaches. This review will give an overview of the state of the art of HGGE in oocytes and embryos, and its accompanying challenges.

Construction of a transposase accessible chromatin landscape reveals chromatin state of repeat elements and potential causal variant for complex traits in pigs

Background

A comprehensive landscape of chromatin states for multiple mammalian tissues is essential for elucidating the molecular mechanism underlying regulatory variants on complex traits. However, the genome-wide chromatin accessibility has been only reported in limited tissue types in pigs.

Results

Here we report a genome-wide landscape of chromatin accessibility of 20 tissues in two female pigs at ages of 6 months using ATAC-seq, and identified 557,273 merged peaks, which greatly expanded the pig regulatory element repository. We revealed tissue-specific regulatory elements which were associated with tissue-relevant biological functions. We identified both positive and negative significant correlations between the regulatory elements and gene transcripts, which showed distinct distributions in terms of their strength and distances from corresponding genes. We investigated the presence of transposable elements (TEs) in open chromatin regions across all tissues, these included identifications of porcine endogenous retroviruses (PERVs) exhibiting high accessibility in liver and homology of porcine specific virus sequences to universally accessible transposable elements. Furthermore, we prioritized a potential causal variant for polyunsaturated fatty acid in the muscle.

Conclusions

Our data provides a novel multi-tissues accessible chromatin landscape that serve as an important resource for interpreting regulatory sequences in tissue-specific and conserved biological functions, as well as regulatory variants of loci associated with complex traits in pigs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00767-3.

Low Expression of Mitofusin 1 Gene Leads to Mitochondrial Dysfunction and Embryonic Genome Activation Failure in Ovine-Bovine Inter-Species Cloned Embryos

Inter-species somatic cell nuclear transfer (iSCNT) is significant in the study of biological problems such as embryonic genome activation and the mitochondrial function of embryos. Here, we used iSCNT as a model to determine whether abnormal embryo genome activation was caused by mitochondrial dysfunction. First, we found the ovine-bovine iSCNT embryos were developmentally blocked at the 8-cell stage. The reactive oxygen species level, mitochondrial membrane potential, and ATP level in ovine-bovine cloned embryos were significantly different from both bovine-bovine and IVF 8-cell stage embryos. RNA sequencing and q-PCR analysis revealed that mitochondrial transport, mitochondrial translational initiation, mitochondrial large ribosomal subunit, and mitochondrial outer membrane genes were abnormally expressed in the ovine-bovine embryos, and the mitochondrial outer membrane and mitochondrial ribosome large subunit genes, mitochondrial fusion gene 1, and ATPase Na+/K+ transporting subunit beta 3 gene were expressed at lower levels in the ovine-bovine cloned embryos. Furthermore, we found that overexpression and knockdown of Mfn1 significantly affected mitochondrial fusion and subsequent biological functions such as production of ATP, mitochondrial membrane potential, reactive oxygen species and gene expressions in cloned embryos. These findings enhance our understanding of the mechanism by which the Mfn1 gene regulates embryonic development and embryonic genome activation events.

Ruminant-specific retrotransposons shape regulatory evolution of bovine immunity

Cattle are an important livestock species, and mapping the genomic architecture of agriculturally relevant traits such as disease susceptibility is a major challenge in the bovine research community. Lineage-specific transposable elements (TEs) are increasingly recognized to contribute to gene regulatory evolution and variation, but this possibility has been largely unexplored in ruminant genomes. We conducted epigenomic profiling of the type II interferon (IFN) response in bovine cells and found thousands of ruminant-specific TEs including MER41_BT and Bov-A2 elements predicted to act as IFN-inducible enhancer elements. CRISPR knockout experiments in bovine cells established that critical immune factors including IFNAR2 and IL2RB are transcriptionally regulated by TE-derived enhancers. Finally, population genomic analysis of 38 individuals revealed that a subset of polymorphic TE insertions may function as enhancers in modern cattle. Our study reveals that lineage-specific TEs have shaped the evolution of ruminant IFN responses and potentially continue to contribute to immune gene regulatory differences across modern breeds and individuals. Together with previous work in human cells, our findings demonstrate that lineage-specific TEs have been independently co-opted to regulate IFN-inducible gene expression in multiple species, supporting TE co-option as a recurrent mechanism driving the evolution of IFN-inducible transcriptional networks.

MiR-202-3p determines embryo viability during mid-blastula transition

Developmental growth is an intricate process involving the coordinated regulation of the expression of various genes, and microRNAs (miRNAs) play crucial roles in diverse processes throughout animal development. The mid-blastula transition (MBT) is a developmental milestone when maternal RNAs are cleared and the zygotic genome programmed asynchronous cell division begins to drive embryogenesis. While mechanisms underlying MBT have been intensively revealed, factors regulating cell proliferation at the transition remain largely unknown. We report here a microRNA, miR-202-3p to be a key factor that determines embryonic fate during MBT in zebrafish. A miR-202-3p antagomir specifically terminated embryo development at the mid-blastula stage. In vivo deletion of the miR-202 locus recapitulated the fatal phenotypes, which were rescued only by miR-202-3p or its precursor. Transcriptome comparison revealed >250 RNAs including both maternal and zygotic origins were dysregulated at MBT in the miR-202^−/− embryos, corresponding with arrays of homeostatic disorders leading to massive apoptosis. A trio of genes: nfkbiaa, perp and mgll, known to be intimately involved with cell proliferation and survival, were identified as direct targets of miR-202-3p. Importantly, over- or under-expression of any of the trio led to developmental delay or termination at the blastula or gastrula stages. Furthermore, nfkbiaa and perp were shown to inter-regulate each other. Thus, miR-202-3p mediates a regulatory network whose components interact closely during MBT to determine embryonic viability and development.

Identification of the DNA binding element of ZNFO, an oocyte-specific zinc finger transcription factor in cattle

The maternal effect genes are essential components of oocyte competence, which orchestrate the early developmental events before zygotic genome activation (ZGA). The Krüppel-associated box (KRAB) domain-containing zinc finger proteins (KRAB-ZFPs) constitute the largest transcription factor family in mammals. As a novel maternal effect gene, ZNFO was identified previously in our laboratory. The gene codes for a KRAB-ZFP specifically expressed in bovine oocytes and early embryos and gene silencing experiments have demonstrated that ZNFO is required for early embryonic development in cattle. In the present study, we identified a consensus sequence, ATATCCTGTTTAAACCCC, as the DNA binding element of ZNFO (ZNFOBE) using a library of random oligonucleotides by cyclic amplification of sequence target (CAST) analysis. Sequence-specific binding of ZNFO to the DNA binding element was confirmed by an electrophoretic mobility shift assay (EMSA), and the key nucleotides in the ZNFOBE that are required for specific binding by ZNFO were further determined by a competitive EMSA using mutant competitors. Through a luciferase-based reporter assay, it was confirmed that the interaction between ZNFO and ZNFOBE is required for the repressive function of ZNFO. These results provide an essential step towards the identification of ZNFO regulated genes that play important roles during early embryonic development.

Genome-wide profiling of histone H3K4me3 and H3K27me3 modifications in individual blastocysts by CUT&Tag without a solid support (NON-TiE-UP CUT&Tag)

Individual analysis of the epigenome of preimplantation embryos is useful for characterizing each embryo and for investigating the effects of environmental factors on their epigenome. However, it is difficult to analyze genome-wide epigenetic modifications, especially histone modifications, in a large number of single embryos due to the small number of cells and the complexity of the analysis methods. To solve this problem, we further modified the CUT&Tag method, which can analyze histone modifications in a small number of cells, such that the embryo is handled as a cell mass in the reaction solutions in the absence of the solid-phase magnetic beads that are used for antibody and enzyme reactions in the conventional method (NON-TiE-UP CUT&Tag; NTU-CAT). By using bovine blastocysts as a model, we showed that genome-wide profiles of representative histone modifications, H3K4me3 and H3K27me3, could be obtained by NTU-CAT that are in overall agreement with the conventional chromatin immunoprecipitation-sequencing (ChIP-seq) method, even from single embryos. However, this new approach has limitations that require attention, including false positive and negative peaks and lower resolution for broad modifications. Despite these limitations, we consider NTU-CAT a promising replacement for ChIP-seq with the great advantage of being able to analyze individual embryos.