Search PubMedSearch

SEARCH · Search PubMed

Results for “transcription factor evolution”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2Linked to original sources

Host-driven evolution shapes the polysaccharide utilization profiles of alga-associated Flavobacteriaceae.

BACKGROUND: Marine algae represent major producers of complex polysaccharides and serve as hosts for diverse microbial communities in the phycosphere. Flavobacteriaceae are among the key bacterial taxa involved in polysaccharide degradation and carbon remineralization in this environment. However, the extent to which algal hosts drive the divergence of polysaccharide utilization profiles in these bacteria remains unclear. RESULTS: We conducted a genome-resolved analysis of 103 cultured Flavobacteriaceae strains isolated from red, green, and brown macroalgae, as well as from diatoms and dinoflagellates. We found that macroalga-associated strains generally harbored more abundant and diverse CAZyme-encoding genes than their microalga-associated counterparts. Moreover, strains associated with different algal phyla showed distinct metabolic specializations that aligned with the typical polysaccharides of their respective hosts, strongly supporting host-specific adaptation. In four widely distributed genera (Maribacter, Flagellimonas, Polaribacter, Winogradskyella), CAZyme profile dissimilarity and key glycoside hydrolase gene divergence exhibited phylogenetic congruence with algal host phylogeny (Mantel r up to 0.76 and 0.85, respectively), indicative of host-associated functional adaptation. Using Maribacter as a model, cultivation experiments and transcriptome characterization demonstrated that polysaccharide utilization efficiency is not solely linked to the organization of genes into polysaccharide utilization loci (PULs), but also associated with the expression dynamics of key transcription factors (TFs), particularly those from AraC and DeoR families, whose expression patterns were coordinated with laminarin degradation. Notably, these two TF families also exhibited host-associated divergence patterns similar to those of CAZyme-encoding genes. Furthermore, analysis of the Tara Oceans metagenomic data indicated that, within the AraC and DeoR families, a higher proportion of genes were positively correlated with chlorophyll a content compared to other TF families, reinforcing their specialized roles in alga-associated bacterial lifestyles. CONCLUSIONS: Our integrative genomic and transcriptomic analyses reveal evolutionary and regulatory adaptation of marine Flavobacteriaceae to distinct algal hosts. These findings highlight algae-derived habitats as specialized niches that shape microbial metabolic potential, and suggest that carbohydrate metabolism plays a key role in host-driven bacterial evolution across global oceans. Video Abstract.

Flavobacteriaceae

Polyploidy-mediated variations in glutamate receptor proteins linked to Fusarium wilt resistance in upland cotton.

Cotton production in the US faces a serious threat from Fusarium oxysporum f. sp. vasinfectum race 4 (FOV4), a soil-borne fungus causing Fusarium wilt by infecting the roots and vascular system of susceptible cotton, leading to rapid wilting and death. Here, we investigate genetic mechanisms of resistance to FOV4 in the highly resistant upland cotton genotype "U1" using an early-generation segregating biparental population ("U1" × "CSX8308") with comprehensive genomic resources. Reference-grade genomic assemblies of the parents revealed minor structural variations between "U1" haplotypes, a high degree of collinearity at chromosome synteny and micro-synteny levels, and significant divergence from "CSX8308" with 8.9 million SNPs. QTL analysis identified significant markers on chromosomes D03 and A02 linked to reduced Fusarium wilt severity. Within these regions, two glutamate-receptor-like (GLR) genes showed structural variation and overlapped between translocated segments on A02 and D03, suggesting a rare but important reinforcing effect of parallel evolution between susceptible and resistant genotypes. Transcriptome profiles of "U1" under FOV4 infection reveal activation of calcium-binding proteins and transcription factors regulating plant hormones (ethylene, abscisic acid, jasmonic acid, and salicylic acid), along with enzymes involved in cell wall remodeling and phytoalexin production. Advancing cotton improvement depends on incorporating durable genetic disease resistance into high-yielding, high-quality cultivars.

Fusarium

DNA sequence arrangement and preliminary evidence on its evolution.

Some recent measurements of the sequence arrangement and evolution of the eukaryotic genome are reviewed. The range of genome sizes and extent of sequence transcribed into nuclear and messenger RNA indicate that the majority of the single copy DNA is not made up of structural genes. The rate of base substitution in the single copy DNA among the primates is similar to that of the codons for certain rapidly changing amino acid residues. This leads to the hypothesis that there is a "basal" rate of change in the genome not strongly affected by selection. The DNA of most higher animals shows a large amount of short period interspersion of repetitive and single copy DNA sequences and a smaller amount of long repetitive regions. The sequence divergence among the short interspersed repetitive sequences is greater than that of the sequences in long repetitive regions. The long repetitive regions are most probably recent additions to the genome and the short interspersed repetitive sequences result from a history of base substitution and translocation. The process of sequence rearrangement appears to be a significant part of the evolution of the genome and may have a much greater effect on the evolution of the phenotype than sequence alteration by base substitution.

Alleles

Genome-Wide Identification of NLP Family Genes in Cultivated Strawberry (Fragaria × ananassa Duch.) and Analysis of Their Expression Under Heat and Botrytis cinerea Stresses.

Nodule inception (NIN)-like proteins (NLPs) are plant-specific transcription factors regulating nutrient absorption, growth, and stress tolerance; however, their roles in stress responses remain largely uncharacterized. Cultivated strawberry (Fragaria × ananassa 'Camarosa') serves as an ideal model for dissecting the evolution and function of NLP genes. In this study, 37 FaNLP genes were identified genome-wide. Phylogenetic analysis classified them into three subfamilies, which are evenly distributed across seven chromosomes. Divergent exon-intron structures and conserved motif compositions suggest functional differentiation among FaNLPs. Quantitative real-time PCR (qRT-PCR) revealed distinct expression profiles under heat stress and Botrytis cinerea infection. Notably, FaNLPs were significantly more upregulated in the cultivar 'Shuxing' than in 'Benihoppe'. Heat stress inhibited photosynthesis and altered catalase activity (CAT), superoxide dismutase activities (SOD) and peroxidase activities (POD), whereas fungal infection enhanced chitinase activity in both cultivars. Comparative genomics with Arabidopsis and rice revealed strawberry-specific evolutionary patterns of NLPs. Subcellular localization prediction indicates that FaNLP proteins primarily localize to the nucleus, implying their potential roles as transcriptional regulators. This study links FaNLP sequence characteristics with stress response phenotypes, providing a foundation for elucidating NLP-mediated regulatory networks in strawberry. Future functional assays, including overexpression and knockout analyses, will further clarify the biological roles of FaNLPs.

Fragaria

Immune-mediated indirect interaction between gut microbiota and bacterial pathogens.

BACKGROUND: In many animals, survival during infection depends on the ability to coordinate interactions between the host immune system and gut microbiota. These tripartite interactions, in turn, potentially shape pathogen virulence evolution. A key regulator of the immune system and, hence, bipartite interactions in insects is the immune deficiency (Imd) pathway, which modulates gut microbiota and pathogens by synthesizing antimicrobial peptides (AMPs) through the NF-κB transcription factor Relish. However, whether Imd-dependent AMPs mediate indirect interactions between gut microbiota and pathogens in a tripartite context remains unclear. Using RNAi-mediated knockdown of Tenebrio molitor Relish (TmRelish), we hypothesized that Imd-dependent AMPs influence indirect interaction between Providencia burhodogranariea_B (P. b_B) infection and the gut microbiota. RESULTS: TmRelish knockdown altered bipartite interactions by disrupting gut microbiota load and composition, increasing pathogen load, and ultimately leading to higher host mortality during infection. However, we did not find support for our tripartite hypothesis that Imd-dependent AMPs mediate indirect interactions between the gut microbiota and P. b_B infection, suggesting the involvement of alternative regulatory pathways or Imd-independent mechanisms. Nevertheless, our investigations of tripartite interactions showed a positive effect of P. b_B infection on gut microbiota load, which in turn stimulated the expression of a subset of AMPs. However, this upregulation of AMPs did not result in reduced P. b_B load. Notably, the gut microbiota did not affect pathogen load but promoted host survival during P. b_B infection, indicating a role in increasing host tolerance rather than resistance. CONCLUSIONS: These findings suggest that while Imd-dependent AMPs may not mediate tripartite interactions in our system, microbiota-host interactions, such as microbiota-mediated immune priming and changes in microbiota load, can shape infection outcomes. These effects on infection outcomes almost certainly exert important selective pressures on the evolution of bacterial virulence.

Animals

Phytoplasma-plant interactions: effector-mediated host reprogramming, hormonal crosstalk, metabolic alterations and plant-mediated vector manipulation.

Phytoplasmas are wall-less, phloem-restricted bacterial pathogens that infect over 1,000 plant species, causing substantial losses in agriculture, horticulture, and forestry worldwide. Despite their reduced genomes and limited metabolic autonomy, these obligate parasites colonize diverse hosts through secreted effector proteins that extensively reprogram plant development, metabolism, immune signalling, and vector interactions. Advances in genomics, transcriptomics, proteomics, metabolomics, and functional studies have substantially clarified the molecular basis of phytoplasma pathogenicity and symptom development. This review synthesizes current understanding of phytoplasma-plant interactions, covering phytoplasma biology, genome evolution, and the infection cycle across plant and insect vector hosts. We examine the molecular functions of key effectors, SAP11, SAP54/PHYL1, SAP05, TENGU, SWP1, and recently identified virulence factors, focusing on how they target host transcription factors, phytohormone networks, protein degradation pathways, and immune responses to promote colonization and disease progression. We further discuss how phytoplasma infection disrupts phytohormone signalling, primary and secondary metabolism, and developmental programs to produce characteristic disease symptoms, with particular attention to pathogen-induced changes in host volatiles and nutritional quality that alter vector behaviour and enhance transmission. Finally, we summarize insights from multi-omics studies and emerging management strategies, including CRISPR-based genome editing, RNAi, rapid molecular diagnostics, resistant cultivars, microbiome-based approaches, and sustainable vector control, and highlight key knowledge gaps and priorities for developing effective, environmentally sustainable phytoplasma disease management.

Phytoplasma

Mirror worlds: The shared regulatory architecture of cell fate in development and cancer.

Lineage plasticity has emerged as a central mechanism through which cancer cells adapt to therapeutic pressure, evade immune surveillance, and acquire aggressive phenotypes. Although recognized across tumor types, the regulatory principles governing how cancer cells reprogram cellular identity remain incompletely understood. In this review, we propose that lineage plasticity in cancer reflects the redeployment of regulatory frameworks established during normal development. Rather than representing a stochastic byproduct of genomic instability, cancer plasticity frequently unfolds within gene regulatory architectures that also govern cell fate specification, lineage commitment, and controlled state transitions during embryogenesis and tissue homeostasis. Developmental transcription factors, including members of the SOX family, FOXA1, ASCL1, NKX2-1, and epithelial-mesenchymal transition regulators, function as lineage gatekeepers during development but are repurposed in cancer to destabilize lineage commitment and enable phenotypic switching. Similarly, epigenetic regulators that guide developmental trajectories, including chromatin remodeling complexes, Polycomb group proteins, and DNA methylation machinery, are frequently dysregulated or redistributed in tumors, altering the repression of lineage-stabilizing and alternative lineage programs and thereby weakening epigenetic barriers to lineage transitions. Together, these observations support a model in which development and cancer operate as mirror regulatory systems: one establishing and stabilizing cellular identity, the other exploiting the same regulatory architecture to permit adaptive reprogramming under selective pressure. We further discuss how emerging single-cell and spatial multi-omics technologies, integrated with artificial intelligence-based modeling, enable mapping of cell state landscapes and transitional trajectories, transforming lineage plasticity from a descriptive phenomenon into a measurable and predictable property of tumor evolution.

Humans

A mammalian tripartite enhancer cluster controls hypothalamic Pomc expression, food intake, and body weight.

Food intake and energy balance are tightly regulated by a group of hypothalamic arcuate neurons expressing the proopiomelanocortin (POMC) gene. In mammals, arcuate-specific POMC expression is driven by two cis-acting transcriptional enhancers known as nPE1 and nPE2. Because mutant mice lacking these two enhancers still showed hypothalamic Pomc mRNA, we searched for additional elements contributing to arcuate Pomc expression. By combining molecular evolution with reporter gene expression in transgenic zebrafish and mice, here, we identified a mammalian arcuate-specific Pomc enhancer that we named nPE3, carrying several binding sites also present in nPE1 and nPE2 for transcription factors known to activate neuronal Pomc expression, such as ISL1, NKX2.1, and ER&#x3b1;. We found that nPE3 originated in the lineage leading to placental mammals and remained under purifying selection in all mammalian orders, although it was lost in Simiiformes (monkeys, apes, and humans) following a unique segmental deletion event. Interestingly, ablation of nPE3 from the mouse genome led to a drastic reduction (>70%) in hypothalamic Pomc mRNA during development and only moderate (<33%) in adult mice. Comparison between double (nPE1 and nPE2) and triple (nPE1, nPE2, and nPE3) enhancer mutants revealed the relative contribution of nPE3 to hypothalamic Pomc expression and its importance in the control of food intake and adiposity in male and female mice. Altogether, these results demonstrate that nPE3 integrates a tripartite cluster of partially redundant enhancers that originated upon a triple convergent evolutionary process in mammals and that is critical for hypothalamic Pomc expression and body weight homeostasis.

Animals

Missense variants in human forkhead transcription factors reveal determinants of forkhead DNA bispecificity.

Recognition of specific DNA sequences by transcription factors (TFs) is a key step in transcriptional control of gene expression. While most forkhead (FH) TFs bind either an FKH (RYAAAYA) or an FHL (GACGC) recognition motif, some FHs can bind both motifs. Mechanisms that control whether an FH is monospecific vs. bispecific have remained unknown. Screening a library of 12 reference FH proteins, 61 naturally occurring missense variants including clinical variants, and 22 designed mutant FHs for DNA-binding activity using universal ("all 10-mer") protein-binding microarrays revealed non-DNA-contacting residues that control mono- vs. bispecificity. Variation in non-DNA-contacting amino acid residues of TFs is associated with human traits and may play a role in the evolution of TF DNA-binding activities and gene regulatory networks.

Humans

Comparative analysis of DREB gene family in buckwheat: the role of FtDREB02 in the delphinidin biosynthesis and drought stress response.

Dehydration response element binding (DREB) transcription factors play a pivotal role in plant abiotic stress responses, but its evolutionary and functional characterization in buckwheat remains unexplored. Here, we conducted a comprehensive analysis of the DREB gene family across three buckwheat species, revealing segmental duplication as the primary driver of family expansion and potential purifying selection during evolution. A FtDREB02 gene, classified as group A2, was identified through genome-wide association analysis (GWAS) on drought tolerance and delphinidin content. Functional validation in Arabidopsis thaliana and the hairy root of Tartary buckwheat (Fagopyrum tataricum) demonstrated that overexpression of this gene promotes delphinidin biosynthesis and enhances plant resistance to water scarcity. Through the integration of DAP-seq and PEG transcriptome cluster analysis, a FtANS candidate was screened. Functional studies showed that FtDREB02 regulates delphinidin content by binding directly to DRE elements of the FtANS promoter. This research identifies and comprehensively analyzes the DREB family within buckwheat species, elucidating the regulatory mechanisms of FtDREB02 in controlling flavonoid biosynthesis and drought resistance, providing potential genetic resources for breeding buckwheat varieties with excellent agronomic traits.

Anthocyanins

Functional convergence of regulatory regions provides vital insights into mammalian gliding adaptation.

Uncovering the key genetic basis of complex phenotypic convergence in distantly related species has been a long-standing focus in evolutionary biology and genetics, and the convergent evolution of gliding in mammals offers a valuable opportunity to address this question. Here, we investigated the genomic basis of convergent evolution of gliding in mammals by analyzing both protein-coding genes and conserved non-coding elements (CNEs). We first de novo assembled and annotated two chromosome-level genomes of gliding mammals, the red and white giant flying squirrel (Petaurista alborufus) and sugar gliders (Petaurus breviceps), and conducted comprehensive comparative genomic analysis combined with another gliding mammal, the Sunda flying lemur (Galeopterus variegatus) and 14 background species. We found that the convergent evolution of protein-coding genes provided relatively limited but functionally relevant evidence linked to gliding phenotypes. By contrast, we found that gliding-accelerated CNEs (GACNEs) cluster near functionally equivalent genes and frequently aggregate into highly diverged yet functionally convergent hotspot regions. Across the three gliding lineages, both GACNEs and hotspot GACNEs show strong convergence in their functional enrichment profiles, suggesting a broad genetic basis underlying the convergent gliding phenotype. Furthermore, we identified 72 core transcription factors underpinning the genetic basis of gliding convergence, including EMX2 and ZFHX3, potentially involved in multiple aspects of gliding adaptation. Our study highlights the role of functional convergence in regulatory regions as a key mechanism in mammalian gliding convergence, offering valuable insights and strategies for uncovering the genetic basis of complex convergent traits, thereby advancing understanding of the molecular basis of convergent traits.

Petaurista alborufus

The genome of Lespedeza potaninii reveals biased subgenome evolution and drought adaptation.

Lespedeza potaninii, a xerophytic subshrub belonging to the legume family, is native to the Tengger Desert and is highly adapted to drought. It has important ecological value due to its drought adaptability, but the underlying molecular mechanisms remain largely unknown. Here, we report a 1.24 Gb chromosome-scale assembly of the L. potaninii genome (contig N50 = 15.75&#x2005;Mb). Our results indicate that L. potaninii underwent an allopolyploid event with 2 subgenomes, A and B, presenting asymmetric evolution and B subgenome dominance. We estimate that the 2 diploid progenitors of L. potaninii diverged around 3.6 million years ago (MYA) and merged around 1.0 MYA. We revealed that the expansion of hub genes associated with drought responses, such as the binding partner 1 of accelerated cell death 11 (ACD11) (BPA1), facilitated environmental adaptations of L. potaninii to desert habitats. We found a novel function of the BPA1 family in abiotic stress tolerance in addition to the known role in regulating the plant immune response, which could improve drought tolerance by positively regulating reactive oxygen species homeostasis in plants. We revealed that bZIP transcription factors could bind to the BPA1 promoter and activate its transcription. Our work fills the genomic data gap in the Lespedeza genus and the tribe Desmodieae, which should provide theoretical support both in the study of drought tolerance and in the molecular breeding of legume crops.

Genome, Plant

Cell fate specification modes shape transcriptome evolution in the highly conserved spiral cleavage.

Early animal development can be remarkably variable, influenced by lineage-specific reproductive strategies and adaptations. Yet, early embryogenesis is also strikingly conserved in certain groups, such as Spiralia. In this clade, a shared cleavage program (i.e., spiral cleavage) and similar cell lineages are ancestral to at least seven phyla. Why early development is so conserved in specific groups and plastic in others is not fully understood. Here, we investigated two annelid species (Owenia fusiformis and Capitella teleta) with spiral cleavage but different modes of specifying their primary progenitor cells. By generating high-resolution transcriptomic time courses from the oocyte to gastrulation, we demonstrate that transcriptional dynamics differ markedly between these species during spiral cleavage and instead reflect their distinct timings of embryonic organiser specification. However, the end of cleavage and gastrulation exhibit high transcriptomic similarity, when orthologous transcription factors share gene expression domains, suggesting this period is a previously overlooked mid-developmental transition in annelid embryogenesis. Together, our data reveal hidden transcriptomic plasticity during spiral cleavage, indicating an evolutionary decoupling of morphological and transcriptomic conservation during early embryogenesis.

Animals

The ETS domain transcription factor Elk-1 contains a novel class of repression domain.

The ETS domain transcription factor Elk-1 serves as an integration point for different mitogen-activated protein (MAP) kinase pathways. Phosphorylation of Elk-1 by MAP kinases triggers its activation. However, while the activation process is well understood, its downregulation-inactivation is less well characterized. The ETS DNA-binding domain plays a role in the downregulation of Elk-dependent promoter activity following mitogenic activation by recruiting the mSin3A-HDAC complex. Here we have identified a novel evolutionarily conserved repression domain in Elk-1, termed the R motif, which serves to reduce the basal transcriptional activity of Elk-1 and dampen its response to mitogenic signals. This domain is highly potent and portable and can repress transcription in trans. The R motif is related to the CRD1 repression domain in p300 and can functionally replace this domain and confer p21(waf1/cip1) inducibility on p300. However, the R motif acts in a context-dependent manner and is not p21(waf1/cip1) responsive in Elk-1. Thus, the Elk-1 R motif and the p300 CRD1 motif represent a new class of repression domains that are regulated in a context-dependent manner.

Amino Acid Motifs

Function-informed transcriptome analysis of Drosophila renal tubule.

BACKGROUND: Comprehensive, tissue-specific, microarray analysis is a potent tool for the identification of tightly defined expression patterns that might be missed in whole-organism scans. We applied such an analysis to Drosophila melanogaster Malpighian (renal) tubule, a defined differentiated tissue. RESULTS: The transcriptome of the D. melanogaster Malpighian tubule is highly reproducible and significantly different from that obtained from whole-organism arrays. More than 200 genes are more than 10-fold enriched and over 1,000 are significantly enriched. Of the top 200 genes, only 18 have previously been named, and only 45% have even estimates of function. In addition, 30 transcription factors, not previously implicated in tubule development, are shown to be enriched in adult tubule, and their expression patterns respect precisely the domains and cell types previously identified by enhancer trapping. Of Drosophila genes with close human disease homologs, 50 are enriched threefold or more, and eight enriched 10-fold or more, in tubule. Intriguingly, several of these diseases have human renal phenotypes, implying close conservation of renal function across 400 million years of divergent evolution. CONCLUSIONS: From those genes that are identifiable, a radically new view of the function of the tubule, emphasizing solute transport rather than fluid secretion, can be obtained. The results illustrate the phenotype gap: historically, the effort expended on a model organism has tended to concentrate on a relatively small set of processes, rather than on the spread of genes in the genome.

Animals

Conserved HSFA1-dependent chromatin dynamics drive heat stress responses in plants.

Eukaryotic organisms remodel chromatin landscapes to regulate gene expression in response to environmental stress. In plants, heat stress (HS) induces widespread chromatin changes, yet the role of heat shock transcription factors (HSFs) in chromatin remodeling and their evolutionary conservation remains unclear. Using Marchantia polymorpha Mphsf mutants and Arabidopsis thaliana Athsfa1s mutants, we identify HSFA1 as a key regulator of HS-induced cis-regulatory element (CRE) accessibility, a mechanism conserved across land plants, mice, and humans. Gene regulatory network modeling reveals parallel transcription factor subnetworks, with MpWRKY10 and MpABI5B acting as indirect and negative HS regulators. We further showed that ABA modulates gene expression in an HSFA1-dependent manner without inducing chromatin remodeling. Finally, we develop a machine learning framework integrating chromatin accessibility and CRE information to predict gene expression across species, revealing stress-responsive regulatory logic at the transcriptional level. These findings provide insights into how TFs coordinate chromatin architecture to drive stress adaptation.

Heat-Shock Response

Loss of neurogenesis in Hydra leads to compensatory regulation of neurogenic and neurotransmission genes in epithelial cells.

Hydra continuously differentiates a sophisticated nervous system made of mechanosensory cells (nematocytes) and sensory-motor and ganglionic neurons from interstitial stem cells. However, this dynamic adult neurogenesis is dispensable for morphogenesis. Indeed animals depleted of their interstitial stem cells and interstitial progenitors lose their active behaviours but maintain their developmental fitness, and regenerate and bud when force-fed. To characterize the impact of the loss of neurogenesis in Hydra, we first performed transcriptomic profiling at five positions along the body axis. We found neurogenic genes predominantly expressed along the central body column, which contains stem cells and progenitors, and neurotransmission genes predominantly expressed at the extremities, where the nervous system is dense. Next, we performed transcriptomics on animals depleted of their interstitial cells by hydroxyurea, colchicine or heat-shock treatment. By crossing these results with cell-type-specific transcriptomics, we identified epithelial genes up-regulated upon loss of neurogenesis: transcription factors (Dlx, Dlx1, DMBX1/Manacle, Ets1, Gli3, KLF11, LMX1A, ZNF436, Shox1), epitheliopeptides (Arminins, PW peptide), neurosignalling components (CAMK1D, DDCl2, Inx1), ligand-ion channel receptors (CHRNA1, NaC7), G-Protein Coupled Receptors and FMRFRL. Hence epitheliomuscular cells seemingly enhance their sensing ability when neurogenesis is compromised. This unsuspected plasticity might reflect the extended multifunctionality of epithelial-like cells in early eumetazoan evolution.

Animals

A time-resolved single-cell roadmap of the logic driving anterior neural crest diversification from neural border to migration stages.

Neural crest cells exemplify cellular diversification from a multipotent progenitor population. However, the full sequence of early molecular choices orchestrating the emergence of neural crest heterogeneity from the embryonic ectoderm remains elusive. Gene-regulatory-networks (GRN) govern early development and cell specification toward definitive neural crest. Here, we combine ultradense single-cell transcriptomes with machine-learning and large-scale transcriptomic and epigenomic experimental validation of selected trajectories, to provide the general principles and highlight specific features of the GRN underlying neural crest fate diversification from induction to early migration stages using Xenopus frog embryos as a model. During gastrulation, a transient neural border zone state precedes the choice between neural crest and placodes which includes multiple converging gene programs. During neurulation, transcription factor connectome, and bifurcation analyses demonstrate the early emergence of neural crest fates at the neural plate stage, alongside an unbiased multipotent-like lineage persisting until epithelial-mesenchymal transition stage. We also decipher circuits driving cranial and vagal neural crest formation and provide a broadly applicable high-throughput validation strategy for investigating single-cell transcriptomes in vertebrate GRNs in development, evolution, and disease.

Animals