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.

Analysis of the DNA-binding properties of Alx1, an evolutionarily conserved regulator of skeletogenesis in echinoderms

Alx1, a homeodomain-containing transcription factor, is a highly conserved regulator of skeletogenesis in echinoderms. In sea urchins, Alx1 plays a central role in the differentiation of embryonic primary mesenchyme cells (PMCs) and positively regulates the transcription of most biomineralization genes expressed by these cells. The alx1 gene arose via duplication and acquired a skeletogenic function distinct from its paralog (alx4) through the exonization of a 41–amino acid motif (the D2 domain). Alx1 and Alx4 contain glutamine-50 paired-type homeodomains, which interact preferentially with palindromic binding sites in vitro. Chromatin immunoprecipitation sequencing (ChIP-seq) studies have shown, however, that Alx1 binds both to palindromic and half sites in vivo. To address this apparent discrepancy and explore the function of the D2 domain, we used an endogenous cis-regulatory module associated with Sp-mtmmpb, a gene that encodes a PMC-specific metalloprotease, to analyze the DNA-binding properties of Alx1. We find that Alx1 forms dimeric complexes on TAAT-containing half sites by a mechanism distinct from the well-known mechanism of dimerization on palindromic sites. We used transgenic reporter assays to analyze the functional roles of half sites in vivo and demonstrate that two sites with partially redundant functions are essential for the PMC-specific activity of the Sp-mtmmpb cis-regulatory module. Finally, we show that the D2 domain influences the DNA-binding properties of Alx1 in vitro, suggesting that the exonization of this motif may have facilitated the acquisition of new transcriptional targets and consequently a novel developmental function.

CRISPR-Cas9 editing of non-coding genomic loci as a means of controlling gene expression in the sea urchin

We seek to manipulate gene function here through CRISPR-Cas9 editing of cis-regulatory sequences, rather than the more typical mutation of coding regions. This approach would minimize secondary effects of cellular responses to nonsense mediated decay pathways or to mutant protein products by premature stops. This strategy also allows for reducing gene activity in cases where a complete gene knockout would result in lethality, and it can be applied to the rapid identification of key regulatory sites essential for gene expression. We tested this strategy here with genes of known function as a proof of concept, and then applied it to examine the upstream genomic region of the germline gene Nanos2 in the sea urchin, Strongylocentrotus purpuratus. We first used CRISPR-Cas9 to target established genomic cis-regulatory regions of the skeletogenic cell transcription factor, Alx1, and the TGF-β signaling ligand, Nodal, which produce obvious developmental defects when altered in sea urchin embryos. Importantly, mutation of cis-activator sites (Alx1) and cis-repressor sites (Nodal) result in the predicted decreased and increased transcriptional output, respectively. Upon identification of efficient gRNAs by genomic mutations, we then used the same validated gRNAs to target a deadCas9-VP64 transcriptional activator to increase Nodal transcription directly. Finally, we paired these new methodologies with a more traditional, GFP reporter construct approach to further our understanding of the transcriptional regulation of Nanos2, a key gene required for germ cell identity in S. purpuratus. With a series of reporter assays, upstream Cas9-promoter targeted mutagenesis, coupled with qPCR and in situ RNA hybridization, we concluded that the promoter of Nanos2 drives strong mRNA expression in the sea urchin embryo, indicating that its primordial germ cell (PGC)-specific restriction may rely instead on post-transcriptional regulation. Overall, we present a proof-of-principle tool-kit of Cas9-mediated manipulations of promoter regions that should be applicable in most cells and embryos for which CRISPR-Cas9 is employed.

Crucial and Overlapping Roles of Six1 and Six2 in Craniofacial Development

SIX1 and SIX2 encode closely related transcription factors of which disruptions have been associated with distinct craniofacial syndromes, with mutations in SIX1 associated with branchiootic syndrome 3 (BOS3) and heterozygous deletions of SIX2 associated with frontonasal dysplasia defects. Whereas mice deficient in Six1 recapitulated most of the developmental defects associated with BOS3, mice lacking Six2 function had no obvious frontonasal defects. We show that Six1 and Six2 exhibit partly overlapping patterns of expression in the developing mouse embryonic frontonasal, maxillary, and mandibular processes. We found that Six1 ^-/- Six2 ^-/- double-mutant mice were born with severe craniofacial deformity not seen in the Six1 ^-/- or Six2 ^-/- single mutants, including skull bone agenesis, midline facial cleft, and syngnathia. Moreover, whereas Six1 ^-/- mice exhibited partial transformation of maxillary zygomatic bone into a mandibular condyle-like structure, Six1 ^-/-Six2 ^+/- mice exhibit significantly increased penetrance of the maxillary malformation. In addition to ectopic Dlx5 expression at the maxillary-mandibular junction as recently reported in E10.5 Six1 ^-/- embryos, the E10.5 Six1 ^-/- Six2 ^+/- embryos showed ectopic expression of Bmp4, Msx1, and Msx2 messenger RNAs in the maxillary-mandibular junction. Genetically inactivating 1 allele of either Ednra or Bmp4 significantly reduced the penetrance of maxillary malformation in both Six1 ^-/- and Six1 ^-/- Six2 ^+/- embryos, indicating that Six1 and Six2 regulate both endothelin and bone morphogenetic protein-4 signaling pathways to pattern the facial structures. Furthermore, we show that neural crest-specific inactivation of Six1 in Six2 ^-/- embryos resulted in midline facial cleft and frontal bone agenesis. We show that Six1 ^-/- Six2 ^-/- embryos exhibit significantly reduced expression of key frontonasal development genes Alx1 and Alx3 as well as increased apoptosis in the developing frontonasal mesenchyme. Together, these results indicate that Six1 and Six2 function partly redundantly to control multiple craniofacial developmental processes and play a crucial neural crest cell-autonomous role in frontonasal morphogenesis.

Hidden genetic history of the Japanese sand dollar Peronella (Echinoidea: Laganidae) revealed by nuclear intron sequences

The marine environment around Japan experienced significant changes during the Cenozoic Era. In this study, we report findings suggesting that this dynamic history left behind traces in the genome of the Japanese sand dollar species Peronella japonica and P. rubra. Although mitochondrial Cytochrome C Oxidase I sequences did not indicate fragmentation of the current local populations of P. japonica around Japan, two different types of intron sequence were found in the Alx1 locus. We inferred that past fragmentation of the populations account for the presence of two types of nuclear sequences as alleles in the Alx1 intron of P. japonica. It is likely that the split populations have intermixed in recent times; hence, we did not detect polymorphisms in the sequences reflecting the current localization of the species. In addition, we found two allelic sequences of theAlx1 intron in the sister species P. rubra. The divergence times of the two types of Alx1 intron sequences were estimated at approximately 14.9 and 4.0 million years ago for P. japonica and P. rubra, respectively. Our study indicates that information from the intron sequences of nuclear genes can enhance our understanding of past genetic events in organisms.

Zebrafish zic2 controls formation of periocular neural crest and choroid fissure morphogenesis

The vertebrate retina develops in close proximity to the forebrain and neural crest-derived cartilages of the face and jaw. Coloboma, a congenital eye malformation, is associated with aberrant forebrain development (holoprosencephaly) and with craniofacial defects (frontonasal dysplasia) in humans, suggesting a critical role for cross-lineage interactions during retinal morphogenesis. ZIC2, a zinc-finger transcription factor, is linked to human holoprosencephaly. We have previously used morpholino assays to show zebrafish zic2 functions in the developing forebrain, retina and craniofacial cartilage. We now report that zebrafish with genetic lesions in zebrafish zic2 orthologs, zic2a and zic2b, develop with retinal coloboma and craniofacial anomalies. We demonstrate a requirement for zic2 in restricting pax2a expression and show evidence that zic2 function limits Hh signaling. RNA-seq transcriptome analysis identified an early requirement for zic2 in periocular neural crest as an activator of alx1, a transcription factor with essential roles in craniofacial and ocular morphogenesis in human and zebrafish. Collectively, these data establish zic2 mutant zebrafish as a powerful new genetic model for in-depth dissection of cell interactions and genetic controls during craniofacial complex development.

microRNA-31 modulates skeletal patterning in the sea urchin embryo

MicroRNAs (miRNAs) are small non-coding RNAs that repress the translation and reduce the stability of target mRNAs in animal cells. microRNA-31 (miR-31) is known to play a role in cancer, bone formation and lymphatic development. However, studies to understand the function of miR-31 in embryogenesis have been limited. We examined the regulatory role of miR-31 in early development using the sea urchin as a model. miR-31 is expressed at all stages of development and its knockdown (KD) disrupts the patterning and function of primary mesenchyme cells (PMCs), which form the embryonic skeleton spicules. We identified that miR-31 directly represses Pmar1, Alx1, Snail and VegfR7 within the PMC gene regulatory network using reporter constructs. Further, blocking the miR-31-mediated repression of Alx1 and/or VegfR7 in the developing embryo resulted in defects in PMC patterning and skeletogenesis. The majority of the mislocalized PMCs in miR-31 KD embryos did not express VegfR10, indicating that miR-31 regulates VegfR gene expression within PMCs. In addition, miR-31 indirectly suppresses Vegf3 expression in the ectoderm. These results indicate that miR-31 coordinately suppresses genes within the PMCs and in the ectoderm to impact PMC patterning and skeletogenesis. This study identifies the novel function and molecular mechanism of miR-31-mediated regulation in the developing embryo.