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

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.

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.

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.

Sperm DNA methylation patterns at discrete CpGs and genes involved in embryonic development are related to bull fertility

Background

Despite a multifactorial approach being taken for the evaluation of bull semen quality in many animal breeding centres worldwide, reliable prediction of bull fertility is still a challenge. Recently, attention has turned to molecular mechanisms, which could uncover potential biomarkers of fertility. One of these mechanisms is DNA methylation, which together with other epigenetic mechanisms is essential for the fertilising sperm to drive normal embryo development and establish a viable pregnancy. In this study, we hypothesised that bull sperm DNA methylation patterns are related to bull fertility. We therefore investigated DNA methylation patterns from bulls used in artificial insemination with contrasting fertility scores.

Results

The DNA methylation patterns were obtained by reduced representative bisulphite sequencing from 10 high-fertility bulls and 10 low-fertility bulls, having average fertility scores of − 6.6 and + 6.5%, respectively (mean of the population was zero). Hierarchical clustering analysis did not distinguish bulls based on fertility but did highlight individual differences. Despite this, using stringent criteria (DNA methylation difference ≥ 35% and a q-value < 0.001), we identified 661 differently methylated cytosines (DMCs). DMCs were preferentially located in intergenic regions, introns, gene downstream regions, repetitive elements, open sea, shores and shelves of CpG islands. We also identified 10 differently methylated regions, covered by 7 unique genes (SFRP1, STXBP4, BCR, PSMG4, ARSG, ATP11A, RXRA), which are involved in spermatogenesis and early embryonic development.

Conclusion

This study demonstrated that at specific CpG sites, sperm DNA methylation status is related to bull fertility, and identified seven differently methylated genes in sperm of subfertile bulls that may lead to altered gene expression and potentially influence embryo development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08614-5.

KDM5-mediated redistribution of H3K4me3 is required for oocyte-to-embryonic transition in cattle

Reprogramming of histone modifications is critical to safeguard correct gene expression profile during preimplantation development. Of interest, trimethylation of lysine 4 on histone 3 (H3K4me3) exhibits a unique and dynamic landscape with a potential species-specific feature. Here, we address how it is reprogrammed and its functional significance during oocyte maturation and early embryonic development in cows. Notably, the overall signal of H3K4me3 decreased sharply during embryonic genome activation (EGA). By using low input ChIP-seq, we find widespread broad H3K4me3 domains in oocytes and early cleaved embryos. The broad domains are gradually removed after fertilization, which is obviously seen during EGA. Meanwhile, H3K4me3 become enriched at promoter regions after the removal of broad H3K4me3. Interestingly, the gene expression level displays a positive correlation with the relative H3K4me3 signal of their promoters when embryos reach 16-cell stage. Importantly, disruption of H3K4me3 demethylases KDM5 increases H3K4me3 level, decreases the embryonic developmental rate and results in dysregulation of over a thousand genes. Meanwhile, KDM5 deficiency causes a redistribution of H3K4me3 across genome. In particular, H3K4me3 in gene body or intergenic regions can't be removed and H3K4me3 in promoter regions is aberrantly reduced. Besides, the positive correlation between promoter H3K4me3 enrichment and gene expression level disappear. Overall, we describe the genomic reprogramming of H3K4me3 with a greater resolution during bovine preimplantation development and propose that KDM5-mediated redistribution of H3K4me3 plays an important role in modulating oocyte-to-embryonic transition.

Profiling of open chromatin in developing pig (Sus scrofa) muscle to identify regulatory regions

There is very little information about how the genome is regulated in domestic pigs (Sus scrofa). This lack of knowledge hinders efforts to define and predict the effects of genetic variants in pig breeding programs. To address this knowledge gap, we need to identify regulatory sequences in the pig genome starting with regions of open chromatin. We used the "Improved Protocol for the Assay for Transposase-Accessible Chromatin (Omni-ATAC-Seq)" to identify putative regulatory regions in flash-frozen semitendinosus muscle from 24 male piglets. We collected samples from the smallest-, average-, and largest-sized male piglets from each litter through five developmental time points. Of the 4661 ATAC-Seq peaks identified that represent regions of open chromatin, >50% were within 1 kb of known transcription start sites. Differential read count analysis revealed 377 ATAC-Seq defined genomic regions where chromatin accessibility differed significantly across developmental time points. We found regions of open chromatin associated with downregulation of genes involved in muscle development that were present in small-sized fetal piglets but absent in large-sized fetal piglets at day 90 of gestation. The dataset that we have generated provides a resource for studies of genome regulation in pigs and contributes valuable functional annotation information to filter genetic variants for use in genomic selection in pig breeding programs.

The USR domain of USF1 mediates NF-Y interactions and cooperative DNA binding

The trimeric CCAAT-binding NF-Y is a "pioneer" Transcription Factor -TF- known to cooperate with neighboring TFs to regulate gene expression. Genome-wide analyses detected a precise stereo-alignment -10/12 bp- of CCAAT with E-box elements and corresponding colocalization of NF-Y with basic-Helix-Loop-Helix (bHLH) TFs. We dissected here NF-Y interactions with USF1 and MAX. USF1, but not MAX, cooperates in DNA binding with NF-Y. NF-Y and USF1 synergize to activate target promoters. Reconstruction of complexes by structural means shows independent DNA binding of MAX, whereas USF1 has extended contacts with NF-Y, involving the USR, a USF-specific amino acid sequence stretch required for trans-activation. The USR is an intrinsically disordered domain and adopts different conformations based on E-box-CCAAT distances. Deletion of the USR abolishes cooperative DNA binding with NF-Y. Our data indicate that the functionality of certain unstructured domains involves adapting to small variation in stereo-alignments of the multimeric TFs sites.

Evolutionary epigenomic analyses in mammalian early embryos reveal species-specific innovations and conserved principles of imprinting

While mouse remains the most popular model, the conservation of parental-to-embryonic epigenetic transition across mammals is poorly defined. Through analysis of oocytes and early embryos in human, bovine, porcine, rat, and mouse, we revealed remarkable species-specific innovations as no single animal model fully recapitulates the human epigenetic transition. In rodent oocytes, transcription-dependent DNA methylation allows methylation of maternal imprints but not intergenic paternal imprints. Unexpectedly, prevalent DNA hypermethylation, paralleled by H3K36me2/3, also occurs in nontranscribed regions in porcine and bovine oocytes, except for megabase-long “CpG continents (CGCs)” where imprinting control regions preferentially reside. Broad H3K4me3 and H3K27me3 domains exist in nonhuman oocytes, yet only rodent H3K27me3 survives beyond genome activation. Coincidently, regulatory elements preferentially evade H3K27me3 in rodent oocytes, and failure to do so causes aberrant embryonic gene repression. Hence, the diverse mammalian innovations of parental-to-embryonic transition center on a delicate “to-methylate-or-not” balance in establishing imprints while protecting other regulatory regions.

TDG is a pig-specific epigenetic regulator with insensitivity to H3K9 and H3K27 demethylation in nuclear transfer embryos

Pig cloning by somatic cell nuclear transfer (SCNT) frequently undergoes incomplete epigenetic remodeling during the maternal-to-zygotic transition, which leads to a significant embryonic loss before implantation. Here, we generated the first genome-wide landscapes of histone methylation in pig SCNT embryos. Excessive H3K9me3 and H3K27me3, but not H3K4me3, were observed in the genomic regions with unfaithful embryonic genome activation and donor-cell-specific gene silencing. A combination of H3K9 demethylase KDM4A and GSK126, an inhibitor of H3K27me3 writer, were able to remove these epigenetic barriers and restore the global transcriptome in SCNT embryos. More importantly, thymine DNA glycosylase (TDG) was defined as a pig-specific epigenetic regulator for nuclear reprogramming, which was not reactivated by H3K9me3 and H3K27me3 removal. Both combined treatment and transient TDG overexpression promoted DNA demethylation and enhanced the blastocyst-forming rates of SCNT embryos, thus offering valuable methods to increase the cloning efficiency of genome-edited pigs for agricultural and biomedical purposes.

•

Identification of reprogramming-resistant genes and regions in porcine SCNT embryos

•

H3K9me3 and H3K27me3 are enriched in reprogramming-resistant genes and regions

•

Removing H3K9me3 and H3K27me3 by KDM4A and GSK126 facilitates nuclear reprogramming

•

Transient TDG overexpression promotes DNA demethylation and improves reprogramming

Conserved and species-specific chromatin remodeling and regulatory dynamics during mouse and chicken limb bud development

Chromatin remodeling and genomic alterations impact spatio-temporal regulation of gene expression, which is central to embryonic development. The analysis of mouse and chicken limb development provides important insights into the morphoregulatory mechanisms, however little is known about the regulatory differences underlying their morphological divergence. Here, we identify the underlying shared and species-specific epigenomic and genomic variations. In mouse forelimb buds, we observe striking synchrony between the temporal dynamics of chromatin accessibility and gene expression, while their divergence in chicken wing buds uncovers species-specific regulatory heterochrony. In silico mapping of transcription factor binding sites and computational footprinting establishes the developmental time-restricted transcription factor-DNA interactions. Finally, the construction of target gene networks for HAND2 and GLI3 transcriptional regulators reveals both conserved and species-specific interactions. Our analysis reveals the impact of genome evolution on the regulatory interactions orchestrating vertebrate limb bud morphogenesis and provides a molecular framework for comparative Evo-Devo studies.

The vertebrate limb bud is a paradigm to uncover the fundamental mechanisms that govern embryogenesis and evolutionary diversification. Here the authors compare mouse and chicken limb bud development to study the impact of genome evolution on conserved and divergent gene regulatory interactions.

Genome activation in equine in vitro-produced embryos

Embryonic genome activation is a critical event in embryo development, in which the transcriptional program of the embryo is initiated. The timing and regulation of this process are species-specific. In vitro embryo production is becoming an important clinical and research tool in the horse; however, very little is known about genome activation in this species. The objective of this work was to identify the timing of genome activation, and the transcriptional networks involved, in in vitro-produced horse embryos. RNA-Seq was performed on oocytes and embryos at eight stages of development (MII, zygote, 2-cell, 4-cell, 8-cell, 16-cell, morula, blastocyst; n = 6 per stage, 2 from each of 3 mares). Transcription of seven genes was initiated at the 2-cell stage. The first substantial increase in gene expression occurred at the 4-cell stage (minor activation), followed by massive gene upregulation and downregulation at the 8-cell stage (major activation). An increase in intronic nucleotides, indicative of transcription initiation, was also observed at the 4-cell stage. Co-expression network analyses identified groups of genes that appeared to be regulated by common mechanisms. Investigation of hub genes and binding motifs enriched in the promoters of co-expressed genes implicated several transcription factors. This work represents, to the best of our knowledge, the first genomic evaluation of embryonic genome activation in horse embryos.

Impact of Global Transcriptional Silencing on Cell Cycle Regulation and Chromosome Segregation in Early Mammalian Embryos

The onset of an early development is, in mammals, characterized by profound changes of multiple aspects of cellular morphology and behavior. These are including, but not limited to, fertilization and the merging of parental genomes with a subsequent transition from the meiotic into the mitotic cycle, followed by global changes of chromatin epigenetic modifications, a gradual decrease in cell size and the initiation of gene expression from the newly formed embryonic genome. Some of these important, and sometimes also dramatic, changes are executed within the period during which the gene transcription is globally silenced or not progressed, and the regulation of most cellular activities, including those mentioned above, relies on controlled translation. It is known that the blastomeres within an early embryo are prone to chromosome segregation errors, which might, when affecting a significant proportion of a cell within the embryo, compromise its further development. In this review, we discuss how the absence of transcription affects the transition from the oocyte to the embryo and what impact global transcriptional silencing might have on the basic cell cycle and chromosome segregation controlling mechanisms.

Genome-wide association analyses identify genotype-by-environment interactions of growth traits in Simmental cattle

Understanding genotype-by-environment interactions (G × E) is crucial to understand environmental adaptation in mammals and improve the sustainability of agricultural production. Here, we present an extensive study investigating the interaction of genome-wide SNP markers with a vast assortment of environmental variables and searching for SNPs controlling phenotypic variance (vQTL) using a large beef cattle dataset. We showed that G × E contribute 10.1%, 3.8%, and 2.8% of the phenotypic variance of birth weight, weaning weight, and yearling weight, respectively. G × E genome-wide association analysis (GWAA) detected a large number of G × E loci affecting growth traits, which the traditional GWAA did not detect, showing that functional loci may have non-additive genetic effects regardless of differences in genotypic means. Further, variance-heterogeneity GWAA detected loci enriched with G × E effects without requiring prior knowledge of the interacting environmental factors. Functional annotation and pathway analysis of G × E genes revealed biological mechanisms by which cattle respond to changes in their environment, such as neurotransmitter activity, hypoxia-induced processes, keratinization, hormone, thermogenic and immune pathways. We unraveled the relevance and complexity of the genetic basis of G × E underlying growth traits, providing new insights into how different environmental conditions interact with specific genes influencing adaptation and productivity in beef cattle and potentially across mammals.

Functional roles of the chromatin remodeler SMARCA5 in mouse and bovine preimplantation embryos†

Upon fertilization, extensive chromatin reprogramming occurs during preimplantation development. Growing evidence reveals species-dependent regulations of this process in mammals. ATP-dependent chromatin remodeling factor SMARCA5 (also known as SNF2H) is required for peri-implantation development in mice. However, the specific functional role of SMARCA5 in preimplantation development and if it is conserved among species remain unclear. Herein, comparative analysis of public RNA-seq datasets reveals that SMARCA5 is universally expressed during oocyte maturation and preimplantation development in mice, cattle, humans, and pigs with species-specific patterns. Immunostaining analysis further describes the temporal and spatial changes of SMARCA5 in both mouse and bovine models. siRNA-mediated SMARCA5 depletion reduces the developmental capability and compromises the specification and differentiation of inner cell mass in mouse preimplantation embryos. Indeed, OCT4 is not restricted into the inner cell mass and the formation of epiblast and primitive endoderm disturbed with reduced NANOG and SOX17 in SMARCA5-deficient blastocysts. RNA-seq analysis shows SMARCA5 depletion causes limited effects on the transcriptomics at the morula stage, however, dysregulates 402 genes, including genes involved in transcription regulation and cell proliferation at the blastocyst stage in mice. By comparison, SMARCA5 depletion does not affect the development through the blastocyst stage but significantly compromises the blastocyst quality in cattle. Primitive endoderm formation is greatly disrupted with reduced GATA6 in bovine blastocysts. Overall, our studies demonstrate the importance of SMARCA5 in fostering the preimplantation development in mice and cattle while there are species-specific effects.