Single-cell transcriptomics reveals evolutionary reconfiguration of embryonic cell fate specification in the sea urchin Heliocidaris erythrogramma.

Altered regulatory interactions during development likely underlie a large fraction of phenotypic diversity within and between species, yet identifying specific evolutionary changes remains challenging. Analysis of single-cell developmental transcriptomes from multiple species provides a powerful framework for unbiased identification of evolutionary changes in developmental mechanisms. Here, we leverage a "natural experiment" in developmental evolution in sea urchins, where a major life history switch recently evolved in the lineage leading to Heliocidaris erythrogramma , precipitating extensive changes in early development. Comparative analyses of scRNA-seq developmental time courses from H. erythrogramma and Lytechinus variegatus (representing the derived and ancestral states respectively) reveals numerous evolutionary changes in embryonic patterning. The earliest cell fate specification events, and the primary signaling center are co-localized in the ancestral dGRN but remarkably, in H. erythrogramma they are spatially and temporally separate. Fate specification and differentiation are delayed in most embryonic cell lineages, although in some cases, these processes are conserved or even accelerated. Comparative analysis of regulator-target gene co-expression is consistent with many specific interactions being preserved but delayed in H. erythrogramma , while some otherwise widely conserved interactions have likely been lost. Finally, specific patterning events are directly correlated with evolutionary changes in larval morphology, suggesting that they are directly tied to the life history shift. Together, these findings demonstrate that comparative scRNA-seq developmental time courses can reveal a diverse set of evolutionary changes in embryonic patterning and provide an efficient way to identify likely candidate regulatory interactions for subsequent experimental validation.

A human-specific enhancer fine-tunes radial glia potency and corticogenesis.

Humans evolved an extraordinarily expanded and complex cerebral cortex, associated with developmental and gene regulatory modifications 1-3 . Human accelerated regions (HARs) are highly conserved genomic sequences with human-specific nucleotide substitutions. Although there are thousands of annotated HARs, their functional contribution to human-specific cortical development is largely unknown 4,5 . HARE5 is a HAR transcriptional enhancer of the WNT signaling receptor Frizzled8 (FZD8) active during brain development 6 . Here, using genome-edited mouse and primate models, we demonstrate that human (Hs) HARE5 fine-tunes cortical development and connectivity by controlling the proliferative and neurogenic capacity of neural progenitor cells (NPCs). Hs-HARE5 knock-in mice have significantly enlarged neocortices containing more neurons. By measuring neural dynamics in vivo we show these anatomical features correlate with increased functional independence between cortical regions. To understand the underlying developmental mechanisms, we assess progenitor fate using live imaging, lineage analysis, and single-cell RNA sequencing. This reveals Hs-HARE5 modifies radial glial progenitor behavior, with increased self-renewal at early developmental stages followed by expanded neurogenic potential. We use genome-edited human and chimpanzee (Pt) NPCs and cortical organoids to assess the relative enhancer activity and function of Hs-HARE5 and Pt-HARE5. Using these orthogonal strategies we show four human-specific variants in HARE5 drive increased enhancer activity which promotes progenitor proliferation. These findings illustrate how small changes in regulatory DNA can directly impact critical signaling pathways and brain development. Our study uncovers new functions for HARs as key regulatory elements crucial for the expansion and complexity of the human cerebral cortex.

Differential gene expression profiling implicates altered network development in rat postnatal day 4 cortex following 4-Methylimidazole (4-MeI) induced maternal seizures.

Compromised maternal health leading to maternal seizures can have adverse effects on the healthy development of offspring. This may be the result of inflammation, hypoxia-ischemia, and altered GABA signaling. The current study examined cortical tissue from F2b (2nd litter of the 2nd generation) postnatal day 4 (PND4) offspring of female Harlan SD rats chronically exposed to the seizuregenic compound, 4-Methylimidazole (0, 750, or 2500 ppm 4-MeI). Maternal seizures were evident only at 2500 ppm 4-MeI. GABA related gene expression as examined by qRT-PCR and whole genome microarray showed no indication of disrupted GABA or glutamatergic signaling. Canonical pathway hierarchical clustering and multi-omics combinatory genomic (CNet) plots of differentially expressed genes (DEG) showed alterations in genes associated with regulatory processes of cell development including neuronal differentiation and synaptogenesis. Functional enrichment analysis showed a similarity of cellular processes across the two exposure groups however, the genes comprising each cluster were primarily unique rather than shared and often showed different directionality. A dose-related induction of cytokine signaling was indicated however, pathways associated with individual cytokine signaling were not elevated, suggesting an alternative involvement of cytokine signaling. Pathways related to growth process and cell signaling showed a negative activation supporting an interpretation of disruption or delay in developmental processes at the 2500 ppm 4-MeI exposure level with maternal seizures. Thus, while GABA signaling was not altered as has been observed with maternal seizures, the pattern of DEG suggested a potential for alteration in neuronal network formation.

Recent reconfiguration of an ancient developmental gene regulatory network in Heliocidaris sea urchins.

Changes in developmental gene regulatory networks (dGRNs) underlie much of the diversity of life, but the evolutionary mechanisms that operate on regulatory interactions remain poorly understood. Closely related species with extreme phenotypic divergence provide a valuable window into the genetic and molecular basis for changes in dGRNs and their relationship to adaptive changes in organismal traits. Here we analyse genomes, epigenomes and transcriptomes during early development in two Heliocidaris sea urchin species that exhibit highly divergent life histories and in an outgroup species. Positive selection and chromatin accessibility modifications within putative regulatory elements are enriched on the branch leading to the derived life history, particularly near dGRN genes. Single-cell transcriptomes reveal a dramatic delay in cell fate specification in the derived state, which also has far fewer open chromatin regions, especially near conserved cell fate specification genes. Experimentally perturbing key transcription factors reveals profound evolutionary changes to early embryonic patterning events, disrupting regulatory interactions previously conserved for ~225 million years. These results demonstrate that natural selection can rapidly reshape developmental gene expression on a broad scale when selective regimes abruptly change. More broadly, even highly conserved dGRNs and patterning mechanisms in the early embryo remain evolvable under appropriate ecological circumstances.

Developmental single-cell transcriptomics in the Lytechinus variegatus sea urchin embryo.

Using scRNA-seq coupled with computational approaches, we studied transcriptional changes in cell states of sea urchin embryos during development to the larval stage. Eighteen closely spaced time points were taken during the first 24 h of development of Lytechinus variegatus (Lv). Developmental trajectories were constructed using Waddington-OT, a computational approach to 'stitch' together developmental time points. Skeletogenic and primordial germ cell trajectories diverged early in cleavage. Ectodermal progenitors were distinct from other lineages by the 6th cleavage, although a small percentage of ectoderm cells briefly co-expressed endoderm markers that indicated an early ecto-endoderm cell state, likely in cells originating from the equatorial region of the egg. Endomesoderm cells also originated at the 6th cleavage and this state persisted for more than two cleavages, then diverged into distinct endoderm and mesoderm fates asynchronously, with some cells retaining an intermediate specification status until gastrulation. Seventy-nine out of 80 genes (99%) examined, and included in published developmental gene regulatory networks (dGRNs), are present in the Lv-scRNA-seq dataset and are expressed in the correct lineages in which the dGRN circuits operate.

Methodologies for Following EMT In Vivo at Single Cell Resolution.

An epithelial-mesenchymal transition (EMT) occurs in almost every metazoan embryo at the time mesoderm begins to differentiate. Several embryos have a long record as models for studying an EMT given that a known population of cells enters the EMT at a known time thereby enabling a detailed study of the process. Often, however, it is difficult to learn the molecular details of these model EMT systems because the transitioning cells are a minority of the population of cells in the embryo and in most cases there is an inability to isolate that population. Here we provide a method that enables an examination of genes expressed before, during, and after the EMT with a focus on just the cells that undergo the transition. Single cell RNA-seq (scRNA-seq) has advanced as a technology making it feasible to study the trajectory of gene expression specifically in the cells of interest, in vivo, and without the background noise of other cell populations. The sea urchin skeletogenic cells constitute only 5% of the total number of cells in the embryo yet with scRNA-seq it is possible to study the genes expressed by these cells without background noise. This approach, though not perfect, adds a new tool for uncovering the mechanism of EMT in this cell type.

Myoepithelial Cells of Submucosal Glands Can Function as Reserve Stem Cells to Regenerate Airways after Injury.

Cells demonstrate plasticity following injury, but the extent of this phenomenon and the cellular mechanisms involved remain underexplored. Using single-cell RNA sequencing (scRNA-seq) and lineage tracing, we uncover that myoepithelial cells (MECs) of the submucosal glands (SMGs) proliferate and migrate to repopulate the airway surface epithelium (SE) in multiple injury models. Specifically, SMG-derived cells display multipotency and contribute to basal and luminal cell types of the SMGs and SE. Ex vivo expanded MECs have the potential to repopulate and differentiate into SE cells when grafted onto denuded airway scaffolds. Significantly, we find that SMG-like cells appear on the SE of both extra- and intra-lobular airways of large animal lungs following severe injury. We find that the transcription factor SOX9 is necessary for MEC plasticity in airway regeneration. Because SMGs are abundant and present deep within airways, they may serve as a reserve cell source for enhancing human airway regeneration.