High-resolution annotation of the mouse preimplantation embryo transcriptome using long-read sequencing

Alternative splicing of CARM1 regulated by LincGET-guided paraspeckles biases the first cell fate in mammalian early embryos

The heterogeneity of CARM1 controls first cell fate bias during early mouse development. However, how this heterogeneity is established is unknown. Here, we show that Carm1 mRNA is of a variety of specific exon-skipping splicing (ESS) isoforms in mouse two-cell to four-cell embryos that contribute to CARM1 heterogeneity. Disruption of paraspeckles promotes the ESS of Carm1 precursor mRNAs (pre-mRNAs). LincGET, but not Neat1, is required for paraspeckle assembly and inhibits the ESS of Carm1 pre-mRNAs in mouse two-cell to four-cell embryos. We further find that LincGET recruits paraspeckles to the Carm1 gene locus through HNRNPU. Interestingly, PCBP1 binds the Carm1 pre-mRNAs and promotes its ESS in the absence of LincGET. Finally, we find that the ESS seen in mouse two-cell to four-cell embryos decreases CARM1 protein levels and leads to trophectoderm fate bias. Our findings demonstrate that alternative splicing of CARM1 has an important role in first cell fate determination.

Self-organization of embryonic stem cells into a reproducible embryo model through epigenome editing

Embryonic stem cells (ESCs) can self-organize in vitro into developmental patterns with spatial organization and molecular similarity to that of early embryonic stages. This self-organization of ESCs requires transmission of signaling cues, via addition of small molecule chemicals or recombinant proteins, to induce distinct embryonic cellular fates and subsequent assembly into structures that can mimic aspects of early embryonic development. During natural embryonic development, different embryonic cell types co-develop together, where each cell type expresses specific fate-inducing transcription factors through activation of non-coding regulatory elements and interactions with neighboring cells. However, previous studies have not fully explored the possibility of engineering endogenous regulatory elements to shape self-organization of ESCs into spatially-ordered embryo models. Here, we hypothesized that cell-intrinsic activation of a minimum number of such endogenous regulatory elements is sufficient to self-organize ESCs into early embryonic models. Our results show that CRISPR-based activation (CRISPRa) of only two endogenous regulatory elements in the genome of pluripotent stem cells is sufficient to generate embryonic patterns that show spatial and molecular resemblance to that of pre-gastrulation mouse embryonic development. Quantitative single-cell live fluorescent imaging showed that the emergence of spatially-ordered embryonic patterns happens through the intrinsic induction of cell fate that leads to an orchestrated collective cellular motion. Based on these results, we propose a straightforward approach to efficiently form 3D embryo models through intrinsic CRISPRa-based epigenome editing and independent of external signaling cues. CRISPRa-Programmed Embryo Models (CPEMs) show highly consistent composition of major embryonic cell types that are spatially-organized, with nearly 80% of the structures forming an embryonic cavity. Single cell transcriptomics confirmed the presence of main embryonic cell types in CPEMs with transcriptional similarity to pre-gastrulation mouse embryos and revealed novel signaling communication links between different embryonic cell types. Our findings offer a programmable embryo model and demonstrate that minimum intrinsic epigenome editing is sufficient to self-organize ESCs into highly consistent pre-gastrulation embryo models.

Nuclear RNA catabolism controls endogenous retroviruses, gene expression asymmetry, and dedifferentiation

Endogenous retroviruses (ERVs) are remnants of ancient parasitic infections and comprise sizable portions of most genomes. Although epigenetic mechanisms silence most ERVs by generating a repressive environment that prevents their expression (heterochromatin), little is known about mechanisms silencing ERVs residing in open regions of the genome (euchromatin). This is particularly important during embryonic development, where induction and repression of distinct classes of ERVs occur in short temporal windows. Here, we demonstrate that transcription-associated RNA degradation by the nuclear RNA exosome and Integrator is a regulatory mechanism that controls the productive transcription of most genes and many ERVs involved in preimplantation development. Disrupting nuclear RNA catabolism promotes dedifferentiation to a totipotent-like state characterized by defects in RNAPII elongation and decreased expression of long genes (gene-length asymmetry). Our results indicate that RNA catabolism is a core regulatory module of gene networks that safeguards RNAPII activity, ERV expression, cell identity, and developmental potency.

Relationship of quantitative reverse transcription polymerase chain reaction (RT-PCR) to RNA Sequencing (RNAseq) transcriptome identifies mouse preimplantation embryo reference genes†

Numerous reference genes for use with quantitative reverse transcription polymerase chain reaction (RT-qPCR) have been used for oocytes, eggs, and preimplantation embryos. However, none are actually suitable because of their large variations in expression between developmental stages. To address this, we produced a standardized and merged RNA sequencing (RNAseq) data set by combining multiple publicly available RNAseq data sets that spanned mouse GV oocytes, MII eggs, and 1-cell, 2-cell, 4-cell, 8-cell, morula, and blastocyst stage embryos to identify transcripts with essentially constant expression across all stages. Their expression was then measured using RT-qPCR, with which they did not exhibit constant expression but instead revealed a fixed quantitative relationship between measurements by the two techniques. From this, the relative amounts of total messenger RNA at each stage from the GV oocyte through blastocyst stages were calculated. The quantitative relationship between measurements by RNAseq and RT-qPCR was then used to find genes predicted to have constant expression across stages in RT-qPCR. Candidates were assessed by RT-qPCR to confirm constant expression, identifying Hmgb3 and Rb1cc1 or the geometric mean of those plus either Taf1d or Cd320 as suitable reference genes. This work not only identified transcripts with constant expression from mouse GV oocytes to blastocysts, but also determined a general quantitative relationship between expression measured by RNAseq and RT-qPCR across stages that revealed the relative levels of total mRNA at each stage. The standardized and merged RNA data set should also prove useful in determining transcript expression in mouse oocytes, eggs, and embryos.

Isoform-resolved transcriptome of the human preimplantation embryo

Human preimplantation development involves extensive remodeling of RNA expression and splicing. However, its transcriptome has been compiled using short-read sequencing data, which fails to capture most full-length mRNAs. Here, we generate an isoform-resolved transcriptome of early human development by performing long- and short-read RNA sequencing on 73 embryos spanning the zygote to blastocyst stages. We identify 110,212 unannotated isoforms transcribed from known genes, including highly conserved protein-coding loci and key developmental regulators. We further identify 17,964 isoforms from 5,239 unannotated genes, which are largely non-coding, primate-specific, and highly associated with transposable elements. These isoforms are widely supported by the integration of published multi-omics datasets, including single-cell 8CLC and blastoid studies. Alternative splicing and gene co-expression network analyses further reveal that embryonic genome activation is associated with splicing disruption and transient upregulation of gene modules. Together, these findings show that the human embryo transcriptome is far more complex than currently known, and will act as a valuable resource to empower future studies exploring development.

Dot1l cooperates with Npm1 to repress endogenous retrovirus MERVL in embryonic stem cells

Dot1l is a histone methyltransferase without a SET domain and is responsible for H3K79 methylation, which marks active transcription. In contradiction, Dot1l also participates in silencing gene expression. The target regions and mechanism of Dot1l in repressing transcription remain enigmatic. Here, we show that Dot1l represses endogenous retroviruses in embryonic stem cells (ESCs). Specifically, the absence of Dot1l led to the activation of MERVL, which is a marker of 2-cell-like cells. In addition, Dot1l deletion activated the 2-cell-like state and predisposed ESCs to differentiate into trophectoderm lineage. Transcriptome analysis revealed activation of 2-cell genes and meiotic genes by Dot1l deletion. Mechanistically, Dot1l interacted with and co-localized with Npm1 on MERVL, and depletion of Npm1 similarly augmented MERVL expression. The catalytic activity and AT-hook domain of Dot1l are important to suppress MERVL. Notably, Dot1l-Npm1 restricts MERVL by regulating protein level and deposition of histone H1. Furthermore, Dot1l is critical for Npm1 to efficiently interact with histone H1 and inhibit ubiquitination of H1 whereas Npm1 is essential for Dot1l to interact with MERVL. Altogether, we discover that Dot1l represses MERVL through chaperoning H1 by collaborating with Npm1. Importantly, our findings shed light on the non-canonical transcriptional repressive role of Dot1l in ESCs.

Deep annotation of long noncoding RNAs by assembling RNA-seq and small RNA-seq data

Long noncoding RNAs (lncRNAs) are increasingly being recognized as modulators in various biological processes. However, due to their low expression, their systematic characterization is difficult to determine. Here, we performed transcript annotation by a newly developed computational pipeline, termed RNA-seq and small RNA-seq combined strategy (RSCS), in a wide variety of cellular contexts. Thousands of high-confidence potential novel transcripts were identified by the RSCS, and the reliability of the transcriptome was verified by analysis of transcript structure, base composition, and sequence complexity. Evidenced by the length comparison, the frequency of the core promoter and the polyadenylation signal motifs, and the locations of transcription start and end sites, the transcripts appear to be full length. Furthermore, taking advantage of our strategy, we identified a large number of endogenous retrovirus-associated lncRNAs, and a novel endogenous retrovirus-lncRNA that was functionally involved in control of Yap1 expression and essential for early embryogenesis was identified. In summary, the RSCS can generate a more complete and precise transcriptome, and our findings greatly expanded the transcriptome annotation for the mammalian community.

A DddA ortholog-based and transactivator-assisted nuclear and mitochondrial cytosine base editors with expanded target compatibility

Bacterial double-stranded DNA (dsDNA) cytosine deaminase DddA_tox-derived cytosine base editor (DdCBE) and its evolved variant, DddA11, guided by transcription-activator-like effector (TALE) proteins, enable mitochondrial DNA (mtDNA) editing at TC or HC (H = A, C, or T) sequence contexts, while it remains relatively unattainable for GC targets. Here, we identified a dsDNA deaminase originated from a Roseburia intestinalis interbacterial toxin (riDddA_tox) and generated CRISPR-mediated nuclear DdCBEs (crDdCBEs) and mitochondrial CBEs (mitoCBEs) using split riDddA_tox, which catalyzed C-to-T editing at both HC and GC targets in nuclear and mitochondrial genes. Moreover, transactivator (VP64, P65, or Rta) fusion to the tail of DddA_tox- or riDddA_tox-mediated crDdCBEs and mitoCBEs substantially improved nuclear and mtDNA editing efficiencies by up to 3.5- and 1.7-fold, respectively. We also used riDddA_tox-based and Rta-assisted mitoCBE to efficiently stimulate disease-associated mtDNA mutations in cultured cells and in mouse embryos with conversion frequencies of up to 58% at non-TC targets.

Exploring differential exon usage via short- and long-read RNA sequencing strategies

Alternative splicing produces various mRNAs, and thereby various protein products, from one gene, impacting a wide range of cellular activities. However, accurate reconstruction and quantification of full-length transcripts using short-reads is limited, due to their length. Long-reads sequencing technologies may provide a solution by sequencing full-length transcripts. We explored the use of both Illumina short-reads and two long Oxford Nanopore Technology (cDNA and Direct RNA) RNA-Seq reads for detecting global differential splicing during mouse embryonic stem cell differentiation, applying several bioinformatics strategies: gene-based, isoform-based and exon-based. We detected the strongest similarity among the sequencing platforms at the gene level compared to exon-based and isoform-based. Furthermore, the exon-based strategy discovered many differential exon usage (DEU) events, mostly in a platform-dependent manner and in non-differentially expressed genes. Thus, the platforms complemented each other in the ability to detect DEUs (i.e. long-reads exhibited an advantage in detecting DEUs at the UTRs, and short-reads detected more DEUs). Exons within 20 genes, detected in one or more platforms, were here validated by PCR, including key differentiation genes, such as Mdb3 and Aplp1. We provide an important analysis resource for discovering transcriptome changes during stem cell differentiation and insights for analysing such data.

Improved zebra finch brain transcriptome identifies novel proteins with sex differences

The zebra finch (Taeniopygia guttata), a representative oscine songbird species, has been widely studied to investigate behavioral neuroscience, most notably the neurobiological basis of vocal learning, a rare trait shared in only a few animal groups including humans. In 2019, an updated zebra finch genome annotation (bTaeGut1_v1.p) was released from the Ensembl database and is substantially more comprehensive than the first version published in 2010. In this study, we utilized the publicly available RNA-seq data generated from Illumina-based short-reads and PacBio single-molecule real-time (SMRT) long-reads to assess the bird transcriptome. To analyze the high-throughput RNA-seq data, we adopted a hybrid bioinformatic approach combining short and long-read pipelines. From our analysis, we added 220 novel genes and 8,134 transcript variants to the Ensembl annotation, and predicted a new proteome based on the refined annotation. We further validated 18 different novel proteins by using mass-spectrometry data generated from zebra finch caudal telencephalon tissue. Our results provide additional resources for future studies of zebra finches utilizing this improved bird genome annotation and proteome.

A developmentally programmed splicing failure contributes to DNA damage response attenuation during mammalian zygotic genome activation

Transition from maternal to embryonic transcriptional control is crucial for embryogenesis. However, alternative splicing regulation during this process remains understudied. Using transcriptomic data from human, mouse, and cow preimplantation development, we show that the stage of zygotic genome activation (ZGA) exhibits the highest levels of exon skipping diversity reported for any cell or tissue type. Much of this exon skipping is temporary, leads to disruptive noncanonical isoforms, and occurs in genes enriched for DNA damage response in the three species. Two core spliceosomal components, Snrpb and Snrpd2, regulate these patterns. These genes have low maternal expression at ZGA and increase sharply thereafter. Microinjection of Snrpb/d2 messenger RNA into mouse zygotes reduces the levels of exon skipping at ZGA and leads to increased p53-mediated DNA damage response. We propose that mammalian embryos undergo an evolutionarily conserved, developmentally programmed splicing failure at ZGA that contributes to the attenuation of cellular responses to DNA damage.

Two Compact Cas9 Ortholog-Based Cytosine Base Editors Expand the DNA Targeting Scope and Applications In Vitro and In Vivo

CRISPR/Cas9-based base editing tools enable precise genomic installation and hold great promise for gene therapy, whereas the big size of Cas9 nucleases and its reliability on specific protospacer adjacent motif (PAM) sequences as well as target site preferences restrict the extensive applications of base editing tools. Here, we generate two cytosine base editors (CBEs) by fusing cytidine deaminases with two compact codon-optimized Cas9 orthologs from Streptococcus_gordonii_str._Challis_substr._CH1 (ancSgo-BE4) and Streptococcus_thermophilus_LMG_18311 (ancSth1a-BE4), which are much smaller than Streptococcus pyogenes (SpCas9) and recognize NNAAAG and NHGYRAA PAM sequences, respectively. Both CBEs display high activity, high fidelity, a different editing window, and low by-products for cytosine base editing with minimal DNA and RNA off-targeting activities in mammalian cells. Moreover, both editors show comparable or higher editing efficiencies than two engineered SpCas9 variant (SpCas9-NG and SpRY)-based CBEs in our tested target sites, which perfectly match the PAM sequences for ancSgo-BE4 or ancSth1a-BE4. In addition, we successfully generate two mouse models harboring clinically relevant mutations at the Ar gene via ancSgo-BE4 and ancSth1a-BE4, which display androgen insensitivity syndrome and/or developmental lethality in founder mice. Thus, the two novel CBEs broaden the base editing tool kits with expanded targeting scope and window for efficient gene modification and applications, respectively.

Maternal gestational mercury exposure in relation to cord blood T cell alterations and placental gene expression signatures

The immunotoxic impacts of mercury during early life is poorly understood. We investigated the associations between gestational mercury exposure and frequency of cord blood T cells as well as placental gene expression. Frequency of natural Treg cells was positively associated with prenatal and postpartum mercury toenail concentrations. Frequency of NKT and activated naïve Th cells was positively associated with prenatal toenail mercury concentrations and number of maternal silver-mercury dental amalgams, respectively. Placental gene expression analyses revealed distinct gene signatures associated with mercury exposure. Decreased placental expression of a histone demethylase, KDM4DL, was associated with both higher prenatal and postpartum maternal toenail mercury levels among male infants and remained statistically significant after adjustment for fish and seafood consumption. The results suggest that gestational exposure to mercury concentrations contribute to alterations in both T cells and gene expression in placenta at birth. These alterations may inform mechanisms of mercury immunotoxicity.

Base editing-mediated perturbation of endogenous PKM1/2 splicing facilitates isoform-specific functional analysis in vitro and in vivo

OBJECTIVES

PKM1 and PKM2, which are generated from the alternative splicing of PKM gene, play important roles in tumourigenesis and embryonic development as rate-limiting enzymes in glycolytic pathway. However, because of the lack of appropriate techniques, the specific functions of the 2 PKM splicing isoforms have not been clarified endogenously yet.

MATERIALS AND METHODS

In this study, we used CRISPR-based base editors to perturbate the endogenous alternative splicing of PKM by introducing mutations into the splicing junction sites in HCT116 cells and zebrafish embryos. Sanger sequencing, agarose gel electrophoresis and targeted deep sequencing assays were utilized for identifying mutation efficiencies and detecting PKM1/2 splicing isoforms. Cell proliferation assays and RNA-seq analysis were performed to describe the effects of perturbation of PKM1/2 splicing in tumour cell growth and zebrafish embryo development.

RESULTS

The splicing sites of PKM, a 5' donor site of GT and a 3' acceptor site of AG, were efficiently mutated by cytosine base editor (CBE; BE4max) and adenine base editor (ABE; ABEmax-NG) with guide RNAs (gRNAs) targeting the splicing sites flanking exons 9 and 10 in HCT116 cells and/or zebrafish embryos. The mutations of the 5' donor sites of GT flanking exons 9 or 10 into GC resulted in specific loss of PKM1 or PKM2 expression as well as the increase in PKM2 or PKM1 respectively. Specific loss of PKM1 promoted cell proliferation of HCT116 cells and upregulated the expression of cell cycle regulators related to DNA replication and cell cycle phase transition. In contrast, specific loss of PKM2 suppressed cell growth of HCT116 cells and resulted in growth retardation of zebrafish. Meanwhile, we found that mutation of PKM1/2 splicing sites also perturbated the expression of non-canonical PKM isoforms and produced some novel splicing isoforms.

CONCLUSIONS

This work proved that CRISPR-based base editing strategy can be used to disrupt the endogenous alternative splicing of genes of interest to study the function of specific splicing isoforms in vitro and in vivo. It also reminded us to notice some novel or undesirable splicing isoforms by targeting the splicing junction sites using base editors. In sum, we establish a platform to perturbate endogenous RNA splicing for functional investigation or genetic correction of abnormal splicing events in human diseases.

DevOmics: an integrated multi-omics database of human and mouse early embryo

Transcriptomic and epigenetic alterations during early embryo development have been proven to play essential roles in regulating the cell fate. Nowadays, advances in single-cell transcriptomics and epigenomics profiling techniques provide large volumes of data for understanding the molecular regulatory mechanisms in early embryos and facilitate the investigation of assisted reproductive technology as well as preimplantation genetic testing. However, the lack of integrated data collection and unified analytic procedures greatly limits their usage in scientific research and clinical application. Hence, it is necessary to establish a database integrating the regulatory information of human and mouse early embryos with unified analytic procedures. Here, we introduce DevOmics (http://devomics.cn/), which contains normalized gene expression, DNA methylation, histone modifications (H3K4me3, H3K9me3, H3K27me3, H3K27ac), chromatin accessibility and 3D chromatin architecture profiles of human and mouse early embryos spanning six developmental stages (zygote, 2cell, 4cell, 8cell, morula and blastocyst (ICM, TE)). The current version of DevOmics provides Search and Advanced Search for retrieving genes a researcher is interested in, Analysis Tools including the differentially expressed genes (DEGs) analysis for acquiring DEGs between different types of samples, allelic explorer for displaying allele-specific gene expression as well as epigenetic modifications and correlation analysis for showing the dynamic changes in different layers of data across developmental stages, as well as Genome Browser and Ortholog for visualization. DevOmics offers a user-friendly website for biologists and clinicians to decipher molecular regulatory mechanisms of human and mouse early embryos.

Long read, isoform aware sequencing of mouse nucleus accumbens after chronic cocaine treatment

To better understand the full-length transcriptome of the nucleus accumbens (NAc)-a key brain reward region-in chronic cocaine treatment, we perform the first single molecule, long-read sequencing analysis using the Iso-seq method to detect 42,114 unique transcripts from mouse NAc polyadenylated RNA. Using GENCODE annotation as a reference, we find that over half of the Iso-seq derived transcripts are annotated, while 46% of them harbor novel splicing events in known genes; around 1% of them correspond to other types of novel transcripts, such as fusion, antisense and intergenic. Approximately 34% of the novel transcripts are matched with a compiled transcriptome assembled from published short-read data from various tissues, with the remaining 69% being unique to NAc. These data provide a more complete picture of the NAc transcriptome than existing annotations and can serve as a comprehensive reference for future transcriptomic analyses of this important brain reward region.