Evolution Of Gene Regulatory Networks Research Articles

Evolution of the gene regulatory network of body axis by enhancer hijacking in amphioxus.

A central goal of evolutionary developmental biology is to decipher the evolutionary pattern of gene regulatory networks (GRNs) that control embryonic development, and the mechanism underlying GRNs evolution. The Nodal signaling that governs the body axes of deuterostomes exhibits a conserved GRN orchestrated principally by Nodal, Gdf1/3, and Lefty. Here we show that this GRN has been rewired in cephalochordate amphioxus. We found that while the amphioxus Gdf1/3 ortholog exhibited nearly no embryonic expression, its duplicate Gdf1/3-like, linked to Lefty, was zygotically expressed in a similar pattern as Lefty. Consistent with this, while Gdf1/3-like mutants showed defects in axial development, Gdf1/3 mutants did not. Further transgenic analyses showed that the intergenic region between Gdf1/3-like and Lefty could drive reporter gene expression as that of the two genes. These results indicated that Gdf1/3-like has taken over the axial development role of Gdf1/3 in amphioxus, possibly through hijacking Lefty enhancers. We finally demonstrated that, to compensate for the loss of maternal Gdf1/3 expression, Nodal has become an indispensable maternal factor in amphioxus and its maternal mutants caused axial defects as Gdf1/3-like mutants. We therefore demonstrated a case that the evolution of GRNs could be triggered by enhancer hijacking events. This pivotal event has allowed the emergence of a new GRN in extant amphioxus, presumably through a stepwise process. In addition, the co-expression of Gdf1/3-like and Lefty achieved by a shared regulatory region may have provided robustness during body axis formation, which provides a selection-based hypothesis for the phenomena called developmental system drift.

Evolution of the Spider Homeobox Gene Repertoire by Tandem and Whole Genome Duplication.

Gene duplication generates new genetic material that can contribute to the evolution of gene regulatory networks and phenotypes. Duplicated genes can undergo subfunctionalization to partition ancestral functions and/or neofunctionalization to assume a new function. We previously found there had been a whole genome duplication (WGD) in an ancestor of arachnopulmonates, the lineage including spiders and scorpions but excluding other arachnids like mites, ticks, and harvestmen. This WGD was evidenced by many duplicated homeobox genes, including two Hox clusters, in spiders. However, it was unclear which homeobox paralogues originated by WGD versus smaller-scale events such as tandem duplications. Understanding this is a key to determining the contribution of the WGD to arachnopulmonate genome evolution. Here we characterized the distribution of duplicated homeobox genes across eight chromosome-level spider genomes. We found that most duplicated homeobox genes in spiders are consistent with an origin by WGD. We also found two copies of conserved homeobox gene clusters, including the Hox, NK, HRO, Irx, and SINE clusters, in all eight species. Consistently, we observed one copy of each cluster was degenerated in terms of gene content and organization while the other remained more intact. Focussing on the NK cluster, we found evidence for regulatory subfunctionalization between the duplicated NK genes in the spider Parasteatoda tepidariorum compared to their single-copy orthologues in the harvestman Phalangium opilio. Our study provides new insights into the relative contributions of multiple modes of duplication to the homeobox gene repertoire during the evolution of spiders and the function of NK genes.

Evolutionary innovation through transcription factor rewiring in microbes is shaped by levels of transcription factor activity, expression, and existing connectivit.

The survival of a population during environmental shifts depends on whether the rate of phenotypic adaptation keeps up with the rate of changing conditions. A common way to achieve this is via change to gene regulatory network (GRN) connections-known as rewiring-that facilitate novel interactions and innovation of transcription factors. To understand the success of rapidly adapting organisms, we therefore need to determine the rules that create and constrain opportunities for GRN rewiring. Here, using an experimental microbial model system with the soil bacterium Pseudomonas fluorescens, we reveal a hierarchy among transcription factors that are rewired to rescue lost function, with alternative rewiring pathways only unmasked after the preferred pathway is eliminated. We identify 3 key properties-high activation, high expression, and preexisting low-level affinity for novel target genes-that facilitate transcription factor innovation. Ease of acquiring these properties is constrained by preexisting GRN architecture, which was overcome in our experimental system by both targeted and global network alterations. This work reveals the key properties that determine transcription factor evolvability, and as such, the evolution of GRNs.

Transcription factor expression levels and environmental signals constrain transcription factor innovation.

Evolutionary innovation of transcription factors frequently drives phenotypic diversification and adaptation to environmental change. Transcription factors can gain or lose connections to target genes, resulting in novel regulatory responses and phenotypes. However the frequency of functional adaptation varies between different regulators, even when they are closely related. To identify factors influencing propensity for innovation, we utilise a Pseudomonas fluorescens SBW25 strain rendered incapable of flagellar mediated motility in soft-agar plates via deletion of the flagellar master regulator (fleQ). This bacterium can evolve to rescue flagellar motility via gene regulatory network rewiring of an alternative transcription factor to rescue activity of FleQ. Previously, we have identified two members (out of 22) of the RpoN-dependent enhancer binding protein (RpoN-EBP) family of transcription factors (NtrC and PFLU1132) that are capable of innovating in this way. These two transcription factors rescue motility repeatably and reliably in a strict hierarchy – with NtrC the only route in a ∆fleQ background, and PFLU1132 the only route in a ∆fleQ∆ntrC background. However, why other members in the same transcription factor family have not been observed to rescue flagellar activity is unclear. Previous work shows that protein homology cannot explain this pattern within the protein family (RpoN-EBPs), and mutations in strains that rescued motility suggested high levels of transcription factor expression and activation drive innovation. We predict that mutations that increase expression of the transcription factor are vital to unlock evolutionary potential for innovation. Here, we construct titratable expression mutant lines for 11 of the RpoN-EBPs in P. fluorescens . We show that in five additional RpoN-EBPs (FleR, HbcR, GcsR, DctD, AauR and PFLU2209), high expression levels result in different mutations conferring motility rescue, suggesting alternative rewiring pathways. Our results indicate that expression levels (and not protein homology) of RpoN-EBPs are a key constraining factor in determining evolutionary potential for innovation. This suggests that transcription factors that can achieve high expression through few mutational changes, or transcription factors that are active in the selective environment, are more likely to innovate and contribute to adaptive gene regulatory network evolution.

Detection of TALEN-mediated genome cleavage during the early embryonic stage of the starfish Patiria pectinifera.

Echinoderms have long been utilized as experimental materials to study the genetic control of developmental processes and their evolution. Among echinoderms, the molecular study of starfish embryos has received considerable attention across research topics such as gene regulatory network evolution and larval regeneration. Recently, experimental techniques to manipulate gene functions have been gradually established in starfish as the feasibility of genome editing methods was reported. However, it is still unclear when these techniques cause genome cleavage during the development of starfish, which is critical to understand the timeframe and applicability of the experiment during early development of starfish. We herein reported that gene functions can be analyzed by the genome editing method TALEN in early embryos, such as the blastula of the starfish Patiria pectinifera. We injected the mRNA of TALEN targeting rar, which was previously constructed, into eggs of P. pectinifera and examined the efficiency of genome cleavage through developmental stages from 6 to 48 hours post fertilization. The results will be key knowledge not only when designing TALEN-based experiments but also when assessing the results.

Whole-genome Duplications and the Long-term Evolution of Gene Regulatory Networks in Angiosperms

Abstract Angiosperms have a complex history of whole-genome duplications (WGDs), with varying numbers and ages of WGD events across clades. These WGDs have greatly affected the composition of plant genomes due to the biased retention of genes belonging to certain functional categories following their duplication. In particular, regulatory genes and genes encoding proteins that act in multiprotein complexes have been retained in excess following WGD. Here, we inferred protein–protein interaction (PPI) networks and gene regulatory networks (GRNs) for seven well-characterized angiosperm species and explored the impact of both WGD and small-scale duplications (SSDs) in network topology by analyzing changes in frequency of network motifs. We found that PPI networks are enriched in WGD-derived genes associated with dosage-sensitive intricate systems, and strong selection pressures constrain the divergence of WGD-derived genes at the sequence and PPI levels. WGD-derived genes in network motifs are mostly associated with dosage-sensitive processes, such as regulation of transcription and cell cycle, translation, photosynthesis, and carbon metabolism, whereas SSD-derived genes in motifs are associated with response to biotic and abiotic stress. Recent polyploids have higher motif frequencies than ancient polyploids, whereas WGD-derived network motifs tend to be disrupted on the longer term. Our findings demonstrate that both WGD and SSD have contributed to the evolution of angiosperm GRNs, but in different ways, with WGD events likely having a more significant impact on the short-term evolution of polyploids.

Single-Cell Heterogeneity in Snake Venom Expression Is Hardwired by Co-Option of Regulators from Progressively Activated Pathways.

The ubiquitous cellular heterogeneity underlying many organism-level phenotypes raises questions about what factors drive this heterogeneity and how these complex heterogeneous systems evolve. Here, we use single-cell expression data from a Prairie rattlesnake (Crotalus viridis) venom gland to evaluate hypotheses for signaling networks underlying snake venom regulation and the degree to which different venom gene families have evolutionarily recruited distinct regulatory architectures. Our findings suggest that snake venom regulatory systems have evolutionarily co-opted trans-regulatory factors from extracellular signal-regulated kinase and unfolded protein response pathways that specifically coordinate expression of distinct venom toxins in a phased sequence across a single population of secretory cells. This pattern of co-option results in extensive cell-to-cell variation in venom gene expression, even between tandemly duplicated paralogs, suggesting this regulatory architecture has evolved to circumvent cellular constraints. While the exact nature of such constraints remains an open question, we propose that such regulatory heterogeneity may circumvent steric constraints on chromatin, cellular physiological constraints (e.g., endoplasmic reticulum stress or negative protein-protein interactions), or a combination of these. Regardless of the precise nature of these constraints, this example suggests that, in some cases, dynamic cellular constraints may impose previously unappreciated secondary constraints on the evolution of gene regulatory networks that favors heterogeneous expression.

Robustness and innovation in synthetic genotype networks

Genotype networks are sets of genotypes connected by small mutational changes that share the same phenotype. They facilitate evolutionary innovation by enabling the exploration of different neighborhoods in genotype space. Genotype networks, first suggested by theoretical models, have been empirically confirmed for proteins and RNAs. Comparative studies also support their existence for gene regulatory networks (GRNs), but direct experimental evidence is lacking. Here, we report the construction of three interconnected genotype networks of synthetic GRNs producing three distinct phenotypes in Escherichia coli. Our synthetic GRNs contain three nodes regulating each other by CRISPR interference and governing the expression of fluorescent reporters. The genotype networks, composed of over twenty different synthetic GRNs, provide robustness in face of mutations while enabling transitions to innovative phenotypes. Through realistic mathematical modeling, we quantify robustness and evolvability for the complete genotype-phenotype map and link these features mechanistically to GRN motifs. Our work thereby exemplifies how GRN evolution along genotype networks might be driving evolutionary innovation.

Development's Inaugural Outstanding Paper Prize.

Development's Inaugural Outstanding Paper Prize.

The ensemble of gene regulatory networks at mutation–selection balance

The evolution of diverse phenotypes both involves and is constrained by molecular interaction networks. When these networks influence patterns of expression, we refer to them as gene regulatory networks (GRNs). Here, we develop a model of GRN evolution analogous to work from quasi-species theory, which is itself essentially the mutation-selection balance model from classical population genetics extended to multiple loci. With this GRN model, we prove that-across a broad spectrum of selection pressures-the dynamics converge to a stationary distribution over GRNs. Next, we show from first principles how the frequency of GRNs at equilibrium is related to the topology of the genotype network, in particular, via a specific network centrality measure termed the eigenvector centrality. Finally, we determine the structural characteristics of GRNs that are favoured in response to a range of selective environments and mutational constraints. Our work connects GRN evolution to quasi-species theory-and thus to classical populations genetics-providing a mechanistic explanation for the observed distribution of GRNs evolving in response to various evolutionary forces, and shows how complex fitness landscapes can emerge from simple evolutionary rules.

Young transposable elements rewired gene regulatory networks in human and chimpanzee hippocampal intermediate progenitors.

The hippocampus is associated with essential brain functions, such as learning and memory. Human hippocampal volume is significantly greater than expected compared with that of non-human apes, suggesting a recent expansion. Intermediate progenitors, which are able to undergo multiple rounds of proliferative division before a final neurogenic division, may have played a role in evolutionary hippocampal expansion. To investigate the evolution of gene regulatory networks underpinning hippocampal neurogenesis in apes, we leveraged the differentiation of human and chimpanzee induced pluripotent stem cells into TBR2 (or EOMES)-positive hippocampal intermediate progenitor cells (hpIPCs). We found that the gene networks active in hpIPCs are significantly different between humans and chimpanzees, with ∼2500 genes being differentially expressed. We demonstrate that species-specific transposon-derived enhancers contribute to these transcriptomic differences. Young transposons, predominantly endogenous retroviruses and SINE-Vntr-Alus (SVAs), were co-opted as enhancers in a species-specific manner. Human-specific SVAs provided substrates for thousands of novel TBR2-binding sites, and CRISPR-mediated repression of these SVAs attenuated the expression of ∼25% of the genes that are upregulated in human intermediate progenitors relative to the same cell population in the chimpanzee.

Evolution of gene regulatory network of C4 photosynthesis in the genus Flaveria reveals the evolutionary status of C3-C4 intermediate species

C4 photosynthesis evolved from ancestral C3 photosynthesis by recruiting pre-existing genes to fulfill new functions. The enzymes and transporters required for the C4 metabolic pathway have been intensively studied and well documented; however, the transcription factors (TFs) that regulate these C4 metabolic genes are not yet well understood. In particular, how the TF regulatory network of C4 metabolic genes was rewired during the evolutionary process is unclear. Here, we constructed gene regulatory networks (GRNs) for four closely evolutionarily related species from the genus Flaveria, which represent four different evolutionary stages of C4 photosynthesis: C3 (F. robusta), type I C3-C4 (F. sonorensis), type II C3-C4 (F. ramosissima), and C4 (F. trinervia). Our results show that more than half of the co-regulatory relationships between TFs and core C4 metabolic genes are species specific. The counterparts of the C4 genes in C3 species were already co-regulated with photosynthesis-related genes, whereas the required TFs for C4 photosynthesis were recruited later. The TFs involved in C4 photosynthesis were widely recruited in the type I C3-C4 species; nevertheless, type II C3-C4 species showed a divergent GRN from C4 species. In line with these findings, a 13CO2 pulse-labeling experiment showed that the CO2 initially fixed into C4 acid was not directly released to the Calvin–Benson–Bassham cycle in the type II C3-C4 species. Therefore, our study uncovered dynamic changes in C4 genes and TF co-regulation during the evolutionary process; furthermore, we showed that the metabolic pathway of the type II C3-C4 species F. ramosissima represents an alternative evolutionary solution to the ammonia imbalance in C3-C4 intermediate species.

Heterosis of fitness and phenotypic variance in the evolution of a diploid gene regulatory network.

Heterosis describes the phenomenon, whereby a hybrid population has higher fitness than an inbred population, which has previously been explained by either Mendelian dominance or overdominance under the general assumption of a simple genotype-phenotype relationship. However, recent studies have demonstrated that genes interact through a complex gene regulatory network (GRN). Furthermore, phenotypic variance is reportedly lower for heterozygotes, and the origin of such variance-related heterosis remains elusive. Therefore, a theoretical analysis linking heterosis to GRN evolution and stochastic gene expression dynamics is required. Here, we investigated heterosis related to fitness and phenotypic variance in a system with interacting genes by numerically evolving diploid GRNs. According to the results, the heterozygote population exhibited higher fitness than the homozygote population, indicating fitness-related heterosis resulting from evolution. In addition, the heterozygote population exhibited lower noise-related phenotypic variance in expression levels than the homozygous population, implying that the heterozygote population is more robust to noise. Furthermore, the distribution of the ratio of heterozygote phenotypic variance to homozygote phenotypic variance exhibited quantitative similarity with previous experimental results. By applying dominance and differential gene expression rather than only a single gene expression model, we confirmed the correlation between heterosis and differential gene expression. We explain our results by proposing that the convex high-fitness region is evolutionarily shaped in the genetic space to gain noise robustness under genetic mixing through sexual reproduction. These results provide new insights into the effects of GRNs on variance-related heterosis and differential gene expression.

Evolution of miRNA-Binding Sites and Regulatory Networks in Cichlids.

The divergence of regulatory regions and gene regulatory network (GRN) rewiring is a key driver of cichlid phenotypic diversity. However, the contribution of miRNA-binding site turnover has yet to be linked to GRN evolution across cichlids. Here, we extend our previous studies by analyzing the selective constraints driving evolution of miRNA and transcription factor (TF)–binding sites of target genes, to infer instances of cichlid GRN rewiring associated with regulatory binding site turnover. Comparative analyses identified increased species-specific networks that are functionally associated to traits of cichlid phenotypic diversity. The evolutionary rewiring is associated with differential models of miRNA- and TF-binding site turnover, driven by a high proportion of fast-evolving polymorphic sites in adaptive trait genes compared with subsets of random genes. Positive selection acting upon discrete mutations in these regulatory regions is likely to be an important mechanism in rewiring GRNs in rapidly radiating cichlids. Regulatory variants of functionally associated miRNA- and TF-binding sites of visual opsin genes differentially segregate according to phylogeny and ecology of Lake Malawi species, identifying both rewired, for example, clade-specific and conserved network motifs of adaptive trait associated GRNs. Our approach revealed several novel candidate regulators, regulatory regions, and three-node motifs across cichlid genomes with previously reported associations to known adaptive evolutionary traits.

Transposable elements contribute to the spatiotemporal microRNA landscape in human brain development.

Transposable elements (TEs) contribute to the evolution of gene regulatory networks and are dynamically expressed throughout human brain development and disease. One gene regulatory mechanism influenced by TEs is the miRNA system of post-transcriptional control. miRNA sequences frequently overlap TE loci and this miRNA expression landscape is crucial for control of gene expression in adult brain and different cellular contexts. Despite this, a thorough investigation of the spatiotemporal expression of TE-embedded miRNAs in human brain development is lacking. Here, we identify a spatiotemporally dynamic TE-embedded miRNA expression landscape between childhood and adolescent stages of human brain development. These miRNAs sometimes arise from two apposed TEs of the same subfamily, such as for L2 or MIR elements, but in the majority of cases stem from solo TEs. They give rise to in silico predicted high-confidence pre-miRNA hairpin structures, likely represent functional miRNAs, and have predicted genic targets associated with neurogenesis. TE-embedded miRNA expression is distinct in the cerebellum when compared to other brain regions, as has previously been described for gene and TE expression. Furthermore, we detect expression of previously nonannotated TE-embedded miRNAs throughout human brain development, suggestive of a previously undetected miRNA control network. Together, as with non-TE-embedded miRNAs, TE-embedded sequences give rise to spatiotemporally dynamic miRNA expression networks, the implications of which for human brain development constitute extensive avenues of future experimental research. To facilitate interactive exploration of these spatiotemporal miRNA expression dynamics, we provide the “Brain miRTExplorer” web application freely accessible for the community.

Spatially and Temporally Distributed Complexity-A Refreshed Framework for the Study of GRN Evolution.

Irrespective of the heuristic value of interpretations of developmental processes in terms of gene regulatory networks (GRNs), larger-angle views often suffer from: (i) an inadequate understanding of the relationship between genotype and phenotype; (ii) a predominantly zoocentric vision; and (iii) overconfidence in a putatively hierarchical organization of animal body plans. Here, we constructively criticize these assumptions. First, developmental biology is pervaded by adultocentrism, but development is not necessarily egg to adult. Second, during development, many unicells undergo transcriptomic profile transitions that are comparable to those recorded in pluricellular organisms; thus, their study should not be neglected from the GRN perspective. Third, the putatively hierarchical nature of the animal body is mirrored in the GRN logic, but in relating genotype to phenotype, independent assessments of the dynamics of the regulatory machinery and the animal’s architecture are required, better served by a combinatorial than by a hierarchical approach. The trade-offs between spatial and temporal aspects of regulation, as well as their evolutionary consequences, are also discussed. Multicellularity may derive from a unicell’s sequential phenotypes turned into different but coexisting, spatially arranged cell types. In turn, polyphenism may have been a crucial mechanism involved in the origin of complex life cycles.

System Principles Governing the Organization, Architecture, Dynamics, and Evolution of Gene Regulatory Networks.

Synthetic biology aims to apply engineering principles for the rational, systematical design and construction of biological systems displaying functions that do not exist in nature or even building a cell from scratch. Understanding how molecular entities interconnect, work, and evolve in an organism is pivotal to this aim. Here, we summarize and discuss some historical organizing principles identified in bacterial gene regulatory networks. We propose a new layer, the concilion, which is the group of structural genes and their local regulators responsible for a single function that, organized hierarchically, coordinate a response in a way reminiscent of the deliberation and negotiation that take place in a council. We then highlight the importance that the network structure has, and discuss that the natural decomposition approach has unveiled the system-level elements shaping a common functional architecture governing bacterial regulatory networks. We discuss the incompleteness of gene regulatory networks and the need for network inference and benchmarking standardization. We point out the importance that using the network structural properties showed to improve network inference. We discuss the advances and controversies regarding the consistency between reconstructions of regulatory networks and expression data. We then discuss some perspectives on the necessity of studying regulatory networks, considering the interactions’ strength distribution, the challenges to studying these interactions’ strength, and the corresponding effects on network structure and dynamics. Finally, we explore the ability of evolutionary systems biology studies to provide insights into how evolution shapes functional architecture despite the high evolutionary plasticity of regulatory networks.

Architecture and evolution of the cis-regulatory system of the echinoderm kirrelL gene.

The gene regulatory network (GRN) that underlies echinoderm skeletogenesis is a prominent model of GRN architecture and evolution. KirrelL is an essential downstream effector gene in this network and encodes an Ig-superfamily protein required for the fusion of skeletogenic cells and the formation of the skeleton. In this study, we dissected the transcriptional control region of the kirrelL gene of the purple sea urchin, Strongylocentrotus purpuratus. Using plasmid- and bacterial artificial chromosome-based transgenic reporter assays, we identified key cis-regulatory elements (CREs) and transcription factor inputs that regulate Sp-kirrelL, including direct, positive inputs from two key transcription factors in the skeletogenic GRN, Alx1 and Ets1. We next identified kirrelL cis-regulatory regions from seven other echinoderm species that together represent all classes within the phylum. By introducing these heterologous regulatory regions into developing sea urchin embryos we provide evidence of their remarkable conservation across ~500 million years of evolution. We dissected in detail the kirrelL regulatory region of the sea star, Patiria miniata, and demonstrated that it also receives direct inputs from Alx1 and Ets1. Our findings identify kirrelL as a component of the ancestral echinoderm skeletogenic GRN. They support the view that GRN subcircuits, including specific transcription factor-CRE interactions, can remain stable over vast periods of evolutionary history. Lastly, our analysis of kirrelL establishes direct linkages between a developmental GRN and an effector gene that controls a key morphogenetic cell behavior, cell-cell fusion, providing a paradigm for extending the explanatory power of GRNs.

Common Themes and Future Challenges in Understanding Gene Regulatory Network Evolution.

A major driving force behind the evolution of species-specific traits and novel structures is alterations in gene regulatory networks (GRNs). Comprehending evolution therefore requires an understanding of the nature of changes in GRN structure and the responsible mechanisms. Here, we review two insect pigmentation GRNs in order to examine common themes in GRN evolution and to reveal some of the challenges associated with investigating changes in GRNs across different evolutionary distances at the molecular level. The pigmentation GRN in Drosophila melanogaster and other drosophilids is a well-defined network for which studies from closely related species illuminate the different ways co-option of regulators can occur. The pigmentation GRN for butterflies of the Heliconius species group is less fully detailed but it is emerging as a useful model for exploring important questions about redundancy and modularity in cis-regulatory systems. Both GRNs serve to highlight the ways in which redeployment of trans-acting factors can lead to GRN rewiring and network co-option. To gain insight into GRN evolution, we discuss the importance of defining GRN architecture at multiple levels both within and between species and of utilizing a range of complementary approaches.

Lessons from a transcription factor: Alx1 provides insights into gene regulatory networks, cellular reprogramming, and cell type evolution.

The skeleton-forming cells of sea urchins and other echinoderms have been studied by developmental biologists as models of cell specification and morphogenesis for many decades. The gene regulatory network (GRN) deployed in the embryonic skeletogenic cells of euechinoid sea urchins is one of the best understood in any developing animal. Recent comparative studies have leveraged the information contained in this GRN, bringing renewed attention to the diverse patterns of skeletogenesis within the phylum and the evolutionary basis for this diversity. The homeodomain-containing transcription factor, Alx1, was originally shown to be a core component of the skeletogenic GRN of the sea urchin embryo. Alx1 has since been found to be key regulator of skeletal cell identity throughout the phylum. As such, Alx1 is currently serving as a lens through which multiple developmental processes are being investigated. These include not only GRN organization and evolution, but also cell reprogramming, cell type evolution, and the gene regulatory control of morphogenesis. This review summarizes our current state of knowledge concerning Alx1 and highlights the insights it is yielding into these important developmental and evolutionary processes.