Diverse ETS transcription factors mediate FGF signaling in the Ciona anterior neural plate

A microfluidic chip for immobilization and imaging of Ciona intestinalis larvae

Sea squirts (Tunicata) are chordates and develop a swimming larva with a small and defined number of individually identifiable cells. This offers the prospect of connecting specific stimuli to behavioral output and characterizing the neural activity that links these together. Here, we describe the development of a microfluidic chip that allows live larvae of the sea squirt Ciona intestinalis to be immobilized and recorded. By generating transgenic larvae expressing GCaAMP6m in defined cells, we show that calcium ion levels can be recorded from immobilized larvae, while microfluidic control allows larvae to be exposed to specific waterborne stimuli. We trial this on sea water carrying increased levels of carbon dioxide, providing evidence that larvae can sense this gas.

Early transcriptional similarities between two distinct neural lineages during ascidian embryogenesis

In chordates, the central nervous system arises from precursors that have distinct developmental and transcriptional trajectories. Anterior nervous systems are ontogenically associated with ectodermal lineages while posterior nervous systems are associated with mesoderm. Taking advantage of the well-documented cell lineage of ascidian embryos, we asked to what extent the transcriptional states of the different neural lineages become similar during the course of progressive lineage restriction. We performed single-cell RNA sequencing (scRNA-seq) analyses on hand-dissected neural precursor cells of the two distinct lineages, together with those of their sister cell lineages, with a high temporal resolution covering five successive cell cycles from the 16-cell to neural plate stages. A transcription factor binding site enrichment analysis of neural specific genes at the neural plate stage revealed limited evidence for shared transcriptional control between the two neural lineages, consistent with their different ontogenies. Nevertheless, PCA analysis and hierarchical clustering showed that, by neural plate stages, the two neural lineages cluster together. Consistent with this, we identified a set of genes enriched in both neural lineages at the neural plate stage, including miR-124, Celf3.a, Zic.r-b, and Ets1/2. Altogether, the current study has revealed genome-wide transcriptional dynamics of neural progenitor cells of two distinct developmental origins. Our scRNA-seq dataset is unique and provides a valuable resource for future analyses, enabling a precise temporal resolution of cell types not previously described from dissociated embryos.

Lhx3/4 initiates a cardiopharyngeal-specific transcriptional program in response to widespread FGF signaling

Individual signaling pathways, such as fibroblast growth factors (FGFs), can regulate a plethora of inductive events. According to current paradigms, signal-dependent transcription factors (TFs), such as FGF/MapK-activated Ets family factors, partner with lineage-determining factors to achieve regulatory specificity. However, many aspects of this model have not been rigorously investigated. One key question relates to whether lineage-determining factors dictate lineage-specific responses to inductive signals or facilitate these responses in collaboration with other inputs. We utilize the chordate model Ciona robusta to investigate mechanisms generating lineage-specific induction. Previous studies in C. robusta have shown that cardiopharyngeal progenitor cells are specified through the combined activity of FGF-activated Ets1/2.b and an inferred ATTA-binding transcriptional cofactor. Here, we show that the homeobox TF Lhx3/4 serves as the lineage-determining TF that dictates cardiopharyngeal-specific transcription in response to pleiotropic FGF signaling. Targeted knockdown of Lhx3/4 leads to loss of cardiopharyngeal gene expression. Strikingly, ectopic expression of Lhx3/4 in a neuroectodermal lineage subject to FGF-dependent specification leads to ectopic cardiopharyngeal gene expression in this lineage. Furthermore, ectopic Lhx3/4 expression disrupts neural plate morphogenesis, generating aberrant cell behaviors associated with execution of incompatible morphogenetic programs. Based on these findings, we propose that combinatorial regulation by signal-dependent and lineage-determinant factors represents a generalizable, previously uncategorized regulatory subcircuit we term "cofactor-dependent induction." Integration of this subcircuit into theoretical models will facilitate accurate predictions regarding the impact of gene regulatory network rewiring on evolutionary diversification and disease ontogeny.

Comparative analysis of PacBio and ONT RNA sequencing methods for Nemopilema Nomurai venom identification

Recent studies on marine organisms have made use of third-generation sequencing technologies such as Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT). While these specialized bioinformatics tools have different algorithmic designs and performance capabilities, they offer scalability and can be applied to various datasets. We investigated the effectiveness of PacBio and ONT RNA sequencing methods in identifying the venom of the jellyfish species Nemopilema nomurai. We conducted a detailed analysis of the sequencing data from both methods, focusing on key characteristics such as CD, alternative splicing, long-chain noncoding RNA, simple sequence repeat, transcription factor, and functional transcript annotation. Our findings indicate that ONT generally produced higher raw data quality in the transcriptome analysis, while PacBio generated longer read lengths. PacBio was found to be superior in identifying CDs and long-chain noncoding RNA, whereas ONT was more cost-effective for predicting alternative splicing events, simple sequence repeats, and transcription factors. Based on these results, we conclude that PacBio is the most specific and sensitive method for identifying venom components, while ONT is the most cost-effective method for studying venogenesis, cnidocyst (venom gland) development, and transcription of virulence genes in jellyfish. Our study has implications for future sequencing technologies in marine jellyfish, and highlights the power of full-length transcriptome analysis in discovering potential therapeutic targets for jellyfish dermatitis.

Nervous system-related gene regulatory networks and functional evolution of ETS proteins across species

The ETS domain transcription factor family is one of the major transcription factor superfamilies that play regulatory roles in development, cell growth, and cancer progression. Although different functions of ETS member proteins in the nervous system have been demonstrated in various studies, their role in neuronal cell differentiation and the evolutionary conservation of its target genes have not yet been extensively studied. In this study, we focused on the regulatory role of ETS transcription factors in neuronal differentiation and their functional evolution by comparative transcriptomics. In order to investigate the regulatory role of ETS transcription factors in neuronal differentiation across species, transcriptional profiles of ETS members and their target genes were investigated by comparing differentially expressed genes and gene regulatory networks, which were analyzed using human, gorilla, mouse, fruit fly and worm transcriptomics datasets. Bioinformatics approaches to examine the evolutionary conservation of ETS transcription factors during neuronal differentiation have shown that ETS member proteins regulate genes associated with neuronal differentiation, nervous system development, axon, and synaptic regulation in different organisms. This study is a comparative transcriptomic study of ETS transcription factors in terms of neuronal differentiation using a gene regulatory network inference algorithm. Overall, a comparison of gene regulation networks revealed that ETS members are indeed evolutionarily conserved in the regulation of neuronal differentiation. Nonetheless, ETS, PEA3, and ELF subfamilies were found to be relatively more active transcription factors in the transcriptional regulation of neuronal differentiation.

New developments in the biology of fibroblast growth factors

The fibroblast growth factor (FGF) family is composed of 18 secreted signaling proteins consisting of canonical FGFs and endocrine FGFs that activate four receptor tyrosine kinases (FGFRs 1-4) and four intracellular proteins (intracellular FGFs or iFGFs) that primarily function to regulate the activity of voltage-gated sodium channels and other molecules. The canonical FGFs, endocrine FGFs, and iFGFs have been reviewed extensively by us and others. In this review, we briefly summarize past reviews and then focus on new developments in the FGF field since our last review in 2015. Some of the highlights in the past 6 years include the use of optogenetic tools, viral vectors, and inducible transgenes to experimentally modulate FGF signaling, the clinical use of small molecule FGFR inhibitors, an expanded understanding of endocrine FGF signaling, functions for FGF signaling in stem cell pluripotency and differentiation, roles for FGF signaling in tissue homeostasis and regeneration, a continuing elaboration of mechanisms of FGF signaling in development, and an expanding appreciation of roles for FGF signaling in neuropsychiatric diseases. This article is categorized under: Cardiovascular Diseases > Molecular and Cellular Physiology Neurological Diseases > Molecular and Cellular Physiology Congenital Diseases > Stem Cells and Development Cancer > Stem Cells and Development.

ERK signaling dissolves ERF repression condensates in living embryos

Phase separation underlies the organization of the nucleus, including the biogenesis of nucleoli and the packaging of heterochromatin. Here we explore the regulation of transcription factor condensates involved in gene repression by ERK signaling in gastrulating embryos of a simple proto-vertebrate (Ciona). ERK signaling induces nuclear export of the transcriptional repressor Ets-2 repressive factor (ERF), which has been linked to various human developmental disorders. Using high-resolution imaging, we show that ERF is localized within discrete nuclear condensates that dissolve upon ERK activation. Interestingly, we observe dynamic pulses of assembly and dissociation during interphase, providing visualization of a nuclear phase separation process regulated by cell signaling. We discuss the implications of these observations for producing sharp on/off switches in gene activity and suppressing noise in cell-cell signaling events.

Comprehensive multi-omics integration identifies differentially active enhancers during human brain development with clinical relevance

BACKGROUND

Non-coding regulatory elements (NCREs), such as enhancers, play a crucial role in gene regulation, and genetic aberrations in NCREs can lead to human disease, including brain disorders. The human brain is a complex organ that is susceptible to numerous disorders; many of these are caused by genetic changes, but a multitude remain currently unexplained. Understanding NCREs acting during brain development has the potential to shed light on previously unrecognized genetic causes of human brain disease. Despite immense community-wide efforts to understand the role of the non-coding genome and NCREs, annotating functional NCREs remains challenging.

METHODS

Here we performed an integrative computational analysis of virtually all currently available epigenome data sets related to human fetal brain.

RESULTS

Our in-depth analysis unravels 39,709 differentially active enhancers (DAEs) that show dynamic epigenomic rearrangement during early stages of human brain development, indicating likely biological function. Many of these DAEs are linked to clinically relevant genes, and functional validation of selected DAEs in cell models and zebrafish confirms their role in gene regulation. Compared to enhancers without dynamic epigenomic rearrangement, DAEs are subjected to higher sequence constraints in humans, have distinct sequence characteristics and are bound by a distinct transcription factor landscape. DAEs are enriched for GWAS loci for brain-related traits and for genetic variation found in individuals with neurodevelopmental disorders, including autism.

CONCLUSION

This compendium of high-confidence enhancers will assist in deciphering the mechanism behind developmental genetics of human brain and will be relevant to uncover missing heritability in human genetic brain disorders.

Brain Sensory Organs of the Ascidian Ciona robusta: Structure, Function and Developmental Mechanisms

During evolution, new characters are designed by modifying pre-existing structures already present in ancient organisms. In this perspective, the Central Nervous System (CNS) of ascidian larva offers a good opportunity to analyze a complex phenomenon with a simplified approach. As sister group of vertebrates, ascidian tadpole larva exhibits a dorsal CNS, made up of only about 330 cells distributed into the anterior sensory brain vesicle (BV), connected to the motor ganglion (MG) and a caudal nerve cord (CNC) in the tail. Low number of cells does not mean, however, low complexity. The larval brain contains 177 neurons, for which a documented synaptic connectome is now available, and two pigmented organs, the otolith and the ocellus, controlling larval swimming behavior. The otolith is involved in gravity perception and the ocellus in light perception. Here, we specifically review the studies focused on the development of the building blocks of ascidians pigmented sensory organs, namely pigment cells and photoreceptor cells. We focus on what it is known, up to now, on the molecular bases of specification and differentiation of both lineages, on the function of these organs after larval hatching during pre-settlement period, and on the most cutting-edge technologies, like single cell RNAseq and genome editing CRISPR/CAS9, that, adapted and applied to Ciona embryos, are increasingly enhancing the tractability of Ciona for developmental studies, including pigmented organs formation.

Single-cell analysis of cell fate bifurcation in the chordate Ciona

BACKGROUND

Inductive signaling interactions between different cell types are a major mechanism for the further diversification of embryonic cell fates. Most blastomeres in the model chordate Ciona robusta become restricted to a single predominant fate between the 64-cell and mid-gastrula stages. The deeply stereotyped and well-characterized Ciona embryonic cell lineages allow the transcriptomic analysis of newly established cell types very early in their divergence from sibling cell states without the pseudotime inference needed in the analysis of less synchronized cell populations. This is the first ascidian study to use droplet scRNAseq with large numbers of analyzed cells as early as the 64-cell stage when major lineages such as primary notochord first become fate restricted.

RESULTS AND CONCLUSIONS

We identify 59 distinct cell states, including new subregions of the b-line neural lineage and the early induction of the tail tip epidermis. We find that 34 of these cell states are directly or indirectly dependent on MAPK-mediated signaling critical to early Ciona patterning. Most of the MAPK-dependent bifurcations are canalized with the signal-induced cell fate lost upon MAPK inhibition, but the posterior endoderm is unique in being transformed into a novel state expressing some but not all markers of both endoderm and muscle. Divergent gene expression between newly bifurcated sibling cell types is dominated by upregulation in the induced cell type. The Ets family transcription factor Elk1/3/4 is uniquely upregulated in nearly all the putatively direct inductions. Elk1/3/4 upregulation together with Ets transcription factor binding site enrichment analysis enables inferences about which bifurcations are directly versus indirectly controlled by MAPK signaling. We examine notochord induction in detail and find that the transition between a Zic/Ets-mediated regulatory state and a Brachyury/FoxA-mediated regulatory state is unexpectedly late. This supports a "broad-hourglass" model of cell fate specification in which many early tissue-specific genes are induced in parallel to key tissue-specific transcriptional regulators via the same set of transcriptional inputs.

Genome-Wide Identification, Comparison, and Expression Analysis of Transcription Factors in Ascidian Styela clava

Tunicates include diverse species, as they are model animals for evolutionary developmental biology study. The embryonic development of tunicates is known to be extensively regulated by transcription factors (TFs). Styela clava, the globally distributed invasive tunicate, exhibits a strong capacity for environmental adaptation. However, the TFs were not systematically identified and analyzed. In this study, we reported 553 TFs categorized into 60 families from S. clava, based on the whole genome data. Comparison of TFs analysis among the tunicate species revealed that the gene number in the zinc finger superfamily displayed the most significant discrepancy, indicating this family was under the highly evolutionary selection and might be related to species differentiation and environmental adaptation. The greatest number of TFs was discovered in the Cys2His2-type zinc finger protein (zf-C2H2) family in S. clava. From the point of temporal view, more than half the TFs were expressed at the early embryonic stage. The expression correlation analysis revealed the existence of a transition for TFs expression from early embryogenesis to the later larval development in S. clava. Eight Hox genes were identified to be located on one chromosome, exhibiting different arrangement and expression patterns, compared to Ciona robusta (C. intestinalis type A). In addition, a total of 23 forkhead box (fox) genes were identified in S. clava, and their expression profiles referred to their potential roles in neurodevelopment and sensory organ development. Our data, thus, provides crucial clues to the potential functions of TFs in development and environmental adaptation in the leathery sea squirt.

Single-cell transcriptome profiling of the Ciona larval brain

The tadpole-type larva of Ciona has emerged as an intriguing model system for the study of neurodevelopment. The Ciona intestinalis connectome has been recently mapped, revealing the smallest central nervous system (CNS) known in any chordate, with only 177 neurons. This minimal CNS is highly reminiscent of larger CNS of vertebrates, sharing many conserved developmental processes, anatomical compartments, neuron subtypes, and even specific neural circuits. Thus, the Ciona tadpole offers a unique opportunity to understand the development and wiring of a chordate CNS at single-cell resolution. Here we report the use of single-cell RNAseq to profile the transcriptomes of single cells isolated by fluorescence-activated cell sorting (FACS) from the whole brain of Ciona robusta (formerly intestinalis Type A) larvae. We have also compared these profiles to bulk RNAseq data from specific subsets of brain cells isolated by FACS using cell type-specific reporter plasmid expression. Taken together, these datasets have begun to reveal the compartment- and cell-specific gene expression patterns that define the organization of the Ciona larval brain.

Transcriptional regulation of the Ciona Gsx gene in the neural plate

The ascidian neural plate consists of a defined number of identifiable cells organized in a grid of rows and columns, representing a useful model to investigate the molecular mechanisms controlling neural patterning in chordates. Distinct anterior brain lineages are specified via unique combinatorial inputs of signalling pathways with Nodal and Delta-Notch signals patterning along the medial-lateral axis and FGF/MEK/ERK signals patterning along the anterior-posterior axis of the neural plate. The Ciona Gsx gene is specifically expressed in the a9.33 cells in the row III/column 2 position of anterior brain lineages, characterised by a combinatorial input of Nodal-OFF, Notch-ON and FGF-ON. Here, we identify the minimal cis-regulatory element (CRE) of 376 bp, which can recapitulate the early activation of Gsx. We show that this minimal CRE responds in the same way as the endogenous Gsx gene to manipulation of FGF- and Notch-signalling pathways and to overexpression of Snail, a mediator of Nodal signals, and Six3/6, which is required to demarcate the anterior boundary of Gsx expression at the late neurula stage. We reveal that sequences proximal to the transcription start site include a temporal regulatory element required for the precise transcriptional onset of gene expression. We conclude that sufficient spatial and temporal information for Gsx expression is integrated in 376 bp of non-coding cis-regulatory sequences.

Multiple inputs into a posterior-specific regulatory network in the Ciona notochord

The gene regulatory networks underlying Ciona notochord fate specification and differentiation have been extensively investigated, but the regulatory basis for regionalized expression within the notochord is not understood. Here we identify three notochord-expressed genes, C11.331, C12.115 and C8.891, with strongly enriched expression in the secondary notochord cells at the posterior tip of the tail. C11.331 and C12.115 share a distinctive expression pattern that is highly enriched in the secondary notochord lineage but also graded within that lineage with the strongest expression at the posterior tip. Both genes show similar responses to pharmacological perturbations of Wnt and FGF signaling, consistent with an important role for Wnt and FGF ligands expressed at the tail tip. Reporter analysis indicates that the C11.331 cis-regulatory regions are extensively distributed, with multiple non-overlapping regions conferring posterior notochord-enriched expression. Fine-scale analysis of a minimal cis-regulatory module identifies discrete positive and negative elements including a strong silencer. Truncation of the silencer region leads to increased expression in the primary notochord, indicating that C11.331 expression is influenced by putative regulators of primary versus secondary notochord fate. The minimal CRM contains predicted ETS, GATA, LMX and Myb sites, all of which lead to reduced expression in secondary notochord when mutated. These results show that the posterior-enriched notochord expression of C11.331 depends on multiple inputs, including Wnt and FGF signals from the tip of the tail, multiple notochord-specific regulators, and yet-to-be identified regulators of regional identity within the notochord.

Open chromatin dynamics reveals stage-specific transcriptional networks in hiPSC-based neurodevelopmental model

Chromatin accessibility to transcription factors (TFs) strongly influences gene transcription and cell differentiation. However, a mechanistic understanding of the transcriptional control during the neuronal differentiation of human induced pluripotent stem cells (hiPSCs), a promising cellular model for mental disorders, remains elusive. Here, we carried out additional analyses on our recently published open chromatin regions (OCRs) profiling at different stages of hiPSC neuronal differentiation. We found that the dynamic changes of OCR during neuronal differentiation highlighted cell stage-specific gene networks, and the chromatin accessibility at the core promoter region of a gene correlates with the corresponding transcript abundance. Within the cell stage-specific OCRs, we identified the binding of cell stage-specific TFs and observed a lag of a neuronal TF binding behind the mRNA expression of the corresponding TF. Interestingly, binding footprints of NEUROD1 and NEUROG2, both of which induce high efficient conversion of hiPSCs to glutamatergic neurons, were among those most enriched in the relatively mature neurons. Furthermore, TF network analysis showed that both NEUROD1 and NEUROG2 were present in the same core TF network specific to more mature neurons, suggesting a pivotal mechanism of epigenetic control of neuronal differentiation and maturation. Our study provides novel insights into the epigenetic control of glutamatergic neurogenesis in the context of TF networks, which may be instrumental to improving hiPSC modeling of neuropsychiatric disorders.

An FGF-driven feed-forward circuit patterns the cardiopharyngeal mesoderm in space and time

In embryos, multipotent progenitors divide to produce distinct progeny and express their full potential. In vertebrates, multipotent cardiopharyngeal progenitors produce second-heart-field-derived cardiomyocytes, and branchiomeric skeletal head muscles. However, the mechanisms underlying these early fate choices remain largely elusive. The tunicate Ciona emerged as an attractive model to study early cardiopharyngeal development at high resolution: through two asymmetric and oriented divisions, defined cardiopharyngeal progenitors produce distinct first and second heart precursors, and pharyngeal muscle (aka atrial siphon muscle, ASM) precursors. Here, we demonstrate that differential FGF-MAPK signaling distinguishes between heart and ASM precursors. We characterize a feed-forward circuit that promotes the successive activations of essential ASM determinants, Hand-related, Tbx1/10 and Ebf. Finally, we show that coupling FGF-MAPK restriction and cardiopharyngeal network deployment with cell divisions defines the timing of gene expression and permits the emergence of diverse cell types from multipotent progenitors.

Nodal and FGF coordinate ascidian neural tube morphogenesis

Formation of the vertebrate neural tube represents one of the premier examples of morphogenesis in animal development. Here, we investigate this process in the simple chordate Ciona intestinalis Previous studies have implicated Nodal and FGF signals in the specification of lateral and ventral neural progenitors. We show that these signals also control the detailed cellular behaviors underlying morphogenesis of the neural tube. Live-imaging experiments show that FGF controls the intercalary movements of ventral neural progenitors, whereas Nodal is essential for the characteristic stacking behavior of lateral cells. Ectopic activation of FGF signaling is sufficient to induce intercalary behaviors in cells that have not received Nodal. In the absence of FGF and Nodal, neural progenitors exhibit a default behavior of sequential cell divisions, and fail to undergo the intercalary and stacking behaviors essential for normal morphogenesis. Thus, cell specification events occurring prior to completion of gastrulation coordinate the morphogenetic movements underlying the organization of the neural tube.

Co-expression of Foxa.a, Foxd and Fgf9/16/20 defines a transient mesendoderm regulatory state in ascidian embryos

In many bilaterian embryos, nuclear β-catenin (nβ-catenin) promotes mesendoderm over ectoderm lineages. Although this is likely to represent an evolutionary ancient developmental process, the regulatory architecture of nβ-catenin-induced mesendoderm remains elusive in the majority of animals. Here, we show that, in ascidian embryos, three nβ-catenin transcriptional targets, Foxa.a, Foxd and Fgf9/16/20, are each required for the correct initiation of both the mesoderm and endoderm gene regulatory networks. Conversely, these three factors are sufficient, in combination, to produce a mesendoderm ground state that can be further programmed into mesoderm or endoderm lineages. Importantly, we show that the combinatorial activity of these three factors is sufficient to reprogramme developing ectoderm cells to mesendoderm. We conclude that in ascidian embryos, the transient mesendoderm regulatory state is defined by co-expression of Foxa.a, Foxd and Fgf9/16/20.

DOI:http://dx.doi.org/10.7554/eLife.14692.001