Whole-transcriptome splicing profiling of E7.5 mouse primary germ layers reveals frequent alternative promoter usage during mouse early embryogenesis

Alternative splicing dynamics during gastrulation in mouse embryo

Alternative splicing (AS) plays an essential role in development, differentiation and carcinogenesis. However, the mechanisms underlying splicing regulation during mouse embryo gastrulation remain unclear. Based on spatial-temporal transcriptome and epigenome data, we detected the dynamics of AS and revealed its regulatory mechanisms across primary germ layers during mouse gastrulation, spanning developmental stages from E6.5 to E7.5. Subsequently, the dynamic expression of splicing factors (SFs) during gastrulation was characterized, while the expression patterns and functions of germ layer-specific SFs were identified. The results indicate that AS and differential alternative splicing events (DASEs) exhibit dynamic changes and are significantly abundant during the late stage of gastrulation. Similarly, SFs demonstrate stage-specific expression, with elevated levels observed during the middle and late stages of gastrulation. Epigenetic signals associated with SFs and AS sites demonstrate significant enrichment and undergo dynamic changes throughout gastrulation. Overall, this study offers a systematic analysis of AS during mouse gastrulation, identifies primary germ layer-specific AS events, and characterizes the expression patterns of SFs and the associated epigenetic signals. These findings enhance the understanding of the mechanisms underlying the formation of the three germ layers during mammalian gastrulation, with a focus on pre-mRNA AS.

Deciphering lineage specification during early embryogenesis in mouse gastruloids using multilayered proteomics

Gastrulation is a critical stage in embryonic development during which the germ layers are established. Advances in sequencing technologies led to the identification of gene regulatory programs that control the emergence of the germ layers and their derivatives. However, proteome-based studies of early mammalian development are scarce. To overcome this, we utilized gastruloids and a multilayered mass spectrometry-based proteomics approach to investigate the global dynamics of (phospho) protein expression during gastruloid differentiation. Our findings revealed many proteins with temporal expression and unique expression profiles for each germ layer, which we also validated using single-cell proteomics technology. Additionally, we profiled enhancer interaction landscapes using P300 proximity labeling, which revealed numerous gastruloid-specific transcription factors and chromatin remodelers. Subsequent degron-based perturbations combined with single-cell RNA sequencing (scRNA-seq) identified a critical role for ZEB2 in mouse and human somitogenesis. Overall, this study provides a rich resource for developmental and synthetic biology communities endeavoring to understand mammalian embryogenesis.

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.

Ethiopian indigenous goats offer insights into past and recent demographic dynamics and local adaptation in sub-Saharan African goats

Knowledge on how adaptive evolution and human socio-cultural and economic interests shaped livestock genomes particularly in sub-Saharan Africa remains limited. Ethiopia is in a geographic region that has been critical in the history of African agriculture with ancient and diverse human ethnicity and bio-climatic conditions. Using 52K genome-wide data analysed in 646 individuals from 13 Ethiopian indigenous goat populations, we observed high levels of genetic variation. Although runs of homozygosity (ROH) were ubiquitous genome-wide, there were clear differences in patterns of ROH length and abundance and in effective population sizes illustrating differences in genome homozygosity, evolutionary history, and management. Phylogenetic analysis incorporating patterns of genetic differentiation and gene flow with ancestry modelling highlighted past and recent intermixing and possible two deep ancient genetic ancestries that could have been brought by humans with the first introduction of goats in Africa. We observed four strong selection signatures that were specific to Arsi-Bale and Nubian goats. These signatures overlapped genomic regions with genes associated with morphological, adaptation, reproduction and production traits due possibly to selection under environmental constraints and/or human preferences. The regions also overlapped uncharacterized genes, calling for a comprehensive annotation of the goat genome. Our results provide insights into mechanisms leading to genome variation and differentiation in sub-Saharan Africa indigenous goats.

Retinoic acid promotes metabolic maturation of human Embryonic Stem Cell-derived Cardiomyocytes

Cardiomyocytes differentiated from human embryonic stem cells (hESCs) represent a promising cell source for heart repair, disease modeling and drug testing. However, improving the differentiation efficiency and maturation of hESC-derived cardiomyocytes (hESC-CMs) is still a major concern. Retinoic acid (RA) signaling plays multiple roles in heart development. However, the effects of RA on cardiomyocyte differentiation efficiency and maturation are still unknown.

Methods: RA was added at different time intervals to identify the best treatment windows for cardiomyocyte differentiation and maturation. The efficiency of cardiomyocyte differentiation was detected by quantitative real-time PCR and flow cytometry. Cardiomyocytes maturation was detected by immunofluorescence staining, metabolic assays and patch clamp to verify structural, metabolic and electrophysiological maturation, respectively. RNA sequencing was used for splicing analysis.

Results: We found that RA treatment at the lateral mesoderm stage (days 2-4) significantly improved cardiomyocyte differentiation, as evidenced by the upregulation of TNNT2, NKX2.5 and MYH6 on day 10 of differentiation. In addition, flow cytometry showed that the proportion of differentiated cardiomyocytes in the RA-treated group was significantly higher than that in control group. RA treatment on days 15-20 increased cardiomyocyte area, sarcomere length, multinucleation and mitochondrial copy number. RNA sequencing revealed RA promoted RNA isoform switch to the maturation-related form. Meanwhile, RA promoted electrophysiological maturation and calcium handling of hESC-CMs. Importantly, RA-treated cardiomyocytes showed decreased glycolysis and enhanced mitochondrial oxidative phosphorylation, with the increased utilization of fatty acid and exogenous pyruvate but not glutamine.

Conclusion: Our data indicated that RA treatment at an early time window (days 2-4) promotes the efficiency of cardiomyocyte differentiation and that RA treatment post beating (days 15-20) promotes cardiomyocyte maturation. The biphasic effects of RA provide new insights for improving cardiomyocyte differentiation and quality.

Targeted deletion of an NRL- and CRX-regulated alternative promoter specifically silences FERM and PDZ domain containing 1 (Frmpd1) in rod photoreceptors

Regulation of cell type-specific gene expression is critical for generating neuronal diversity. Transcriptome analyses have unraveled extensive heterogeneity of transcribed sequences in retinal photoreceptors because of alternate splicing and/or promoter usage. Here we show that Frmpd1 (FERM and PDZ domain containing 1) is transcribed from an alternative promoter specifically in the retina. Electroporation of Frmpd1 promoter region, -505 to +382 bp, activated reporter gene expression in mouse retina in vivo. A proximal promoter sequence (-8 to +33 bp) of Frmpd1 binds to neural retina leucine zipper (NRL) and cone-rod homeobox protein (CRX), two rod-specific differentiation factors, and is necessary for activating reporter gene expression in vitro and in vivo. Clustered regularly interspaced short palindromic repeats/Cas9-mediated deletion of the genomic region, including NRL and CRX binding sites, in vivo completely eliminated Frmpd1 expression in rods and dramatically reduced expression in rod bipolar cells, thereby overcoming embryonic lethality caused by germline Frmpd1 deletion. Our studies demonstrate that a cell type-specific regulatory control region is a credible target for creating loss-of-function alleles of widely expressed genes.

Expression profile and potential functional differentiation of the Speedy/RINGO family in mice

As novel cyclin-dependent kinase (CDK) activators, Speedy/RINGO (hereafter named Speedy) proteins can directly regulate the cell cycle of vertebrates by binding to and activating various CDKs. Previous studies have shown that Speedy genes are highly associated with different types of cancer and other diseases. However, Speedy genes have not been systematically identified in mice, and their function and expression profiles remain elusive, which greatly hinders the functional and mechanistic study of Speedy genes in vivo. Here, we comprehensively identified Speedy genes in the mouse genome. Phylogenetic analysis showed that the Speedy gene family should be divided into three subfamilies, rather than the previously reported two subfamilies. Mice have two of the three subfamilies of Speedy genes, namely, subfamilies A and E. Speedy subfamily C genes have been lost from the mouse genome. By combining experimental and bioinformatics approaches, we found that the genes from subfamilies A and E have different expression profiles, indicating their functional divergence, which was also consistent with the phylogenetic results. The genes belonging to subfamily E showed only slightly different expression profiles, indicating their similar functions. Coexpression network analysis showed that the genes coexpressed with mouse Speedy genes were primarily enriched in reproduction-related mechanisms and there were significant functional differences between genes from subfamilies A and E, further demonstrating functional differentiation. In summary, we provide a comprehensive landscape (from evolution to expression and function) of the Speedy family in mice; we also demonstrate that Speedy genes mainly participate in reproduction-related mechanisms and that they have undergone functional differentiation in mice.

WNT signaling memory is required for ACTIVIN to function as a morphogen in human gastruloids

Self-organization of discrete fates in human gastruloids is mediated by a hierarchy of signaling pathways. How these pathways are integrated in time, and whether cells maintain a memory of their signaling history remains obscure. Here, we dissect the temporal integration of two key pathways, WNT and ACTIVIN, which along with BMP control gastrulation. CRISPR/Cas9-engineered live reporters of SMAD1, 2 and 4 demonstrate that in contrast to the stable signaling by SMAD1, signaling and transcriptional response by SMAD2 is transient, and while necessary for pluripotency, it is insufficient for differentiation. Pre-exposure to WNT, however, endows cells with the competence to respond to graded levels of ACTIVIN, which induces differentiation without changing SMAD2 dynamics. This cellular memory of WNT signaling is necessary for ACTIVIN morphogen activity. A re-evaluation of the evidence gathered over decades in model systems, re-enforces our conclusions and points to an evolutionarily conserved mechanism.