Modulation of tooth regeneration through opposing responses to Wnt and BMP signals in teleosts.

Most vertebrate species undergo tooth replacement throughout adult life. This process is marked by the shedding of existing teeth and the regeneration of tooth organs. However, little is known about the genetic circuitry regulating tooth replacement. Here, we tested whether fish orthologs of genes known to regulate mammalian hair regeneration have effects on tooth replacement. Using two fish species that demonstrate distinct modes of tooth regeneration, threespine stickleback (Gasterosteus aculeatus) and zebrafish (Danio rerio), we found that transgenic overexpression of four different genes changed tooth replacement rates in the direction predicted by a hair regeneration model: Wnt10a and Grem2a increased tooth replacement rate, whereas Bmp6 and Dkk2 strongly inhibited tooth formation. Thus, similar to known roles in hair regeneration, Wnt and BMP signals promote and inhibit regeneration, respectively. Regulation of total tooth number was separable from regulation of replacement rates. RNA sequencing of stickleback dental tissue showed that Bmp6 overexpression resulted in an upregulation of Wnt inhibitors. Together, these data support a model in which different epithelial organs, such as teeth and hair, share genetic circuitry driving organ regeneration.

Evolution of Spatial and Temporal cis-Regulatory Divergence in Sticklebacks.

Cis-regulatory changes are thought to play a major role in adaptation. Threespine sticklebacks have repeatedly colonized freshwater habitats in the Northern Hemisphere, where they have evolved a suite of phenotypes that distinguish them from marine populations, including changes in physiology, behavior, and morphology. To understand the role of gene regulatory evolution in adaptive divergence, here we investigate cis-regulatory changes in gene expression between marine and freshwater ecotypes through allele-specific expression (ASE) in F1 hybrids. Surveying seven ecologically relevant tissues, including three sampled across two developmental stages, we identified cis-regulatory divergence affecting a third of genes, nearly half of which were tissue-specific. Next, we compared allele-specific expression in dental tissues at two timepoints to characterize cis-regulatory changes during development between marine and freshwater fish. Applying a genome-wide test for selection on cis-regulatory changes, we find evidence for lineage-specific selection on several processes between ecotypes, including the Wnt signaling pathway in dental tissues. Finally, we show that genes with ASE, particularly those that are tissue-specific, are strongly enriched in genomic regions of repeated marine-freshwater divergence, supporting an important role for these cis-regulatory differences in parallel adaptive evolution of sticklebacks to freshwater habitats. Altogether, our results provide insight into the cis-regulatory landscape of divergence between stickleback ecotypes across tissues and during development, and support a fundamental role for tissue-specific cis-regulatory changes in rapid adaptation to new environments.

Site pleiotropy of a stickleback Bmp6 enhancer.

Development and regeneration are orchestrated by gene regulatory networks that operate in part through transcriptional enhancers. Although many enhancers are pleiotropic and are active in multiple tissues, little is known about whether enhancer pleiotropy is due to 1) site pleiotropy, in which individual transcription factor binding sites (TFBS) are required for activity in multiple tissues, or 2) multiple distinct sites that regulate expression in different tissues. Here, we investigated the pleiotropy of an intronic enhancer of the stickleback Bone morphogenetic protein 6 (Bmp6) gene. This enhancer was previously shown to regulate evolved changes in tooth number and tooth regeneration, and is highly pleiotropic, with robust activity in both fins and teeth throughout embryonic, larval, and adult life, and in the heart and kidney in adult fish. We tested the hypothesis that the pleiotropy of this enhancer is due to site pleiotropy of an evolutionarily conserved predicted Foxc1 TFBS. Transgenic analysis and site-directed mutagenesis experiments both deleting and scrambling this predicted Foxc1 TFBS revealed that the binding site is required for enhancer activity in both teeth and fins throughout embryonic, larval, and adult development, and in the heart and kidney in adult fish. Collectively these data support a model where the pleiotropy of this Bmp6 enhancer is due to site pleiotropy and this putative binding site is required for enhancer activity in multiple anatomical sites from the embryo to the adult.

Population Genomics of Variegated Toad-Headed Lizard Phrynocephalus versicolor and Its Adaptation to the Colorful Sand of the Gobi Desert.

The variegated toad-headed agama, Phrynocephalus versicolor, lives in the arid landscape of the Chinese Gobi Desert. We analyzed populations from three different locations which vary in substrate color and altitude: Heishankou (HSK), Guazhou County (GZ), and Ejin Banner (EJN). The substrate color is either light-yellow (GZ-y), yellow (EJN-y), or black (HSK-b); the corresponding lizard population colors largely match their substrate in the degree of melanism. We assembled the P. versicolor genome and sequenced over 90 individuals from the three different populations. Genetic divergence between populations corresponds to their geographic distribution. We inferred the genetic relationships among these populations and used selection scans and differential expression to identify genes that show signatures of selection. Slc2a11 and akap12, among other genes, are highly differentiated and may be responsible for pigment adaptation to substrate color in P. versicolor.

Evolved Bmp6 enhancer alleles drive spatial shifts in gene expression during tooth development in sticklebacks.

Mutations in enhancers have been shown to often underlie natural variation but the evolved differences in enhancer activity can be difficult to identify in vivo. Threespine sticklebacks (Gasterosteus aculeatus) are a robust system for studying enhancer evolution due to abundant natural genetic variation, a diversity of evolved phenotypes between ancestral marine and derived freshwater forms, and the tractability of transgenic techniques. Previous work identified a series of polymorphisms within an intronic enhancer of the Bone morphogenetic protein 6 (Bmp6) gene that are associated with evolved tooth gain, a derived increase in freshwater tooth number that arises late in development. Here, we use a bicistronic reporter construct containing a genetic insulator and a pair of reciprocal two-color transgenic reporter lines to compare enhancer activity of marine and freshwater alleles of this enhancer. In older fish, the two alleles drive partially overlapping expression in both mesenchyme and epithelium of developing teeth, but the freshwater enhancer drives a reduced mesenchymal domain and a larger epithelial domain relative to the marine enhancer. In younger fish, these spatial shifts in enhancer activity are less pronounced. Comparing Bmp6 expression by in situ hybridization in developing teeth of marine and freshwater fish reveals similar evolved spatial shifts in gene expression. Together, these data support a model in which the polymorphisms within this enhancer underlie evolved tooth gain by shifting the spatial expression of Bmp6 during tooth development, and provide a general strategy to identify spatial differences in enhancer activity in vivo.

Lamprey lecticans link new vertebrate genes to the origin and elaboration of vertebrate tissues.

The evolution of vertebrates from an invertebrate chordate ancestor involved the evolution of new organs, tissues, and cell types. It was also marked by the origin and duplication of new gene families. If, and how, these morphological and genetic innovations are related is an unresolved question in vertebrate evolution. Hyaluronan is an extracellular matrix (ECM) polysaccharide important for water homeostasis and tissue structure. Vertebrates possess a novel family of hyaluronan binding proteins called Lecticans, and studies in jawed vertebrates (gnathostomes) have shown they function in many of the cells and tissues that are unique to vertebrates. This raises the possibility that the origin and/or expansion of this gene family helped drive the evolution of these vertebrate novelties. In order to better understand the evolution of the lectican gene family, and its role in the evolution of vertebrate morphological novelties, we investigated the phylogeny, genomic arrangement, and expression patterns of all lecticans in the sea lamprey (Petromyzon marinus), a jawless vertebrate. Though both P. marinus and gnathostomes each have four lecticans, our phylogenetic and syntenic analyses are most consistent with the independent duplication of one of more lecticans in the lamprey lineage. Despite the likely independent expansion of the lamprey and gnathostome lectican families, we find highly conserved expression of lecticans in vertebrate-specific and mesenchyme-derived tissues. We also find that, unlike gnathostomes, lamprey expresses its lectican paralogs in distinct subpopulations of head skeleton precursors, potentially reflecting an ancestral diversity of skeletal tissue types. Together, these observations suggest that the ancestral pre-duplication lectican had a complex expression pattern, functioned to support mesenchymal histology, and likely played a role in the evolution of vertebrate-specific cell and tissue types.

Distinct tooth regeneration systems deploy a conserved battery of genes.

BACKGROUND: Vertebrate teeth exhibit a wide range of regenerative systems. Many species, including most mammals, reptiles, and amphibians, form replacement teeth at a histologically distinct location called the successional dental lamina, while other species do not employ such a system. Notably, a 'lamina-less' tooth replacement condition is found in a paraphyletic array of ray-finned fishes, such as stickleback, trout, cod, medaka, and bichir. Furthermore, the position, renewal potential, and latency times appear to vary drastically across different vertebrate tooth regeneration systems. The progenitor cells underlying tooth regeneration thus present highly divergent arrangements and potentials. Given the spectrum of regeneration systems present in vertebrates, it is unclear if morphologically divergent tooth regeneration systems deploy an overlapping battery of genes in their naïve dental tissues. RESULTS: In the present work, we aimed to determine whether or not tooth progenitor epithelia could be composed of a conserved cell type between vertebrate dentitions with divergent regeneration systems. To address this question, we compared the pharyngeal tooth regeneration processes in two ray-finned fishes: zebrafish (Danio rerio) and threespine stickleback (Gasterosteus aculeatus). These two teleost species diverged approximately 250 million years ago and demonstrate some stark differences in dental morphology and regeneration. Here, we find that the naïve successional dental lamina in zebrafish expresses a battery of nine genes (bmpr1aa, bmp6, cd34, gli1, igfbp5a, lgr4, lgr6, nfatc1, and pitx2), while active Wnt signaling and Lef1 expression occur during early morphogenesis stages of tooth development. We also find that, despite the absence of a histologically distinct successional dental lamina in stickleback tooth fields, the same battery of nine genes (Bmpr1a, Bmp6, CD34, Gli1, Igfbp5a, Lgr4, Lgr6, Nfatc1, and Pitx2) are expressed in the basalmost endodermal cell layer, which is the region most closely associated with replacement tooth germs. Like zebrafish, stickleback replacement tooth germs additionally express Lef1 and exhibit active Wnt signaling. Thus, two fish systems that either have an organized successional dental lamina (zebrafish) or lack a morphologically distinct successional dental lamina (sticklebacks) deploy similar genetic programs during tooth regeneration. CONCLUSIONS: We propose that the expression domains described here delineate a highly conserved "successional dental epithelium" (SDE). Furthermore, a set of orthologous genes is known to mark hair follicle epithelial stem cells in mice, suggesting that regenerative systems in other epithelial appendages may utilize a related epithelial progenitor cell type, despite the highly derived nature of the resulting functional organs.

Evolution of the endothelin pathway drove neural crest cell diversification.

Neural crest cells (NCCs) are migratory, multipotent embryonic cells that are unique to vertebrates and form an array of clade-defining adult features. The evolution of NCCs has been linked to various genomic events, including the evolution of new gene-regulatory networks1,2, the de novo evolution of genes3 and the proliferation of paralogous genes during genome-wide duplication events4. However, conclusive functional evidence linking new and/or duplicated genes to NCC evolution is lacking. Endothelin ligands (Edns) and endothelin receptors (Ednrs) are unique to vertebrates3,5,6, and regulate multiple aspects of NCC development in jawed vertebrates7-10. Here, to test whether the evolution of Edn signalling was a driver of NCC evolution, we used CRISPR-Cas9 mutagenesis11 to disrupt edn, ednr and dlx genes in the sea lamprey, Petromyzon marinus. Lampreys are jawless fishes that last shared a common ancestor with modern jawed vertebrates around 500 million years ago12. Thus, comparisons between lampreys and gnathostomes can identify deeply conserved and evolutionarily flexible features of vertebrate development. Using the frog Xenopus laevis to expand gnathostome phylogenetic representation and facilitate side-by-side analyses, we identify ancient and lineage-specific roles for Edn signalling. These findings suggest that Edn signalling was activated in NCCs before duplication of the vertebrate genome. Then, after one or more genome-wide duplications in the vertebrate stem, paralogous Edn pathways functionally diverged, resulting in NCC subpopulations with different Edn signalling requirements. We posit that this new developmental modularity facilitated the independent evolution of NCC derivatives in stem vertebrates. Consistent with this, differences in Edn pathway targets are associated with differences in the oropharyngeal skeleton and autonomic nervous system of lampreys and modern gnathostomes. In summary, our work provides functional genetic evidence linking the origin and duplication of new vertebrate genes with the stepwise evolution of a defining vertebrate novelty.

Evolution of vertebrate gill covers via shifts in an ancient Pou3f3 enhancer.

Whereas the gill chambers of jawless vertebrates open directly into the environment, jawed vertebrates evolved skeletal appendages that drive oxygenated water unidirectionally over the gills. A major anatomical difference between the two jawed vertebrate lineages is the presence of a single large gill cover in bony fishes versus separate covers for each gill chamber in cartilaginous fishes. Here, we find that these divergent patterns correlate with the pharyngeal arch expression of Pou3f3 orthologs. We identify a deeply conserved Pou3f3 arch enhancer present in humans through sharks but undetectable in jawless fish. Minor differences between the bony and cartilaginous fish enhancers account for their restricted versus pan-arch expression patterns. In zebrafish, mutation of Pou3f3 or the conserved enhancer disrupts gill cover formation, whereas ectopic pan-arch Pou3f3b expression generates ectopic skeletal elements resembling the multimeric covers of cartilaginous fishes. Emergence of this Pou3f3 arch enhancer >430 Mya and subsequent modifications may thus have contributed to the acquisition and diversification of gill covers and respiratory strategies during gnathostome evolution.

The origin and diversification of the developmental mechanisms that pattern the vertebrate head skeleton.

The apparent evolvability of the vertebrate head skeleton has allowed a diverse array of shapes, sizes, and compositions of the head in order to better adapt species to their environments. This encompasses feeding, breathing, sensing, and communicating: the head skeleton somehow participated in the evolution of all these critical processes for the last 500 million years. Through evolution, present head diversity was made possible via developmental modifications to the first head skeletal genetic program. Understanding the development of the vertebrate common ancestor's head skeleton is thus an important step in identifying how different lineages have respectively achieved their many innovations in the head. To this end, cyclostomes (jawless vertebrates) are extremely useful, having diverged from jawed vertebrates approximately 400 million years ago, at the deepest node within living vertebrates. From this ancestral vantage point (that is, the node connecting cyclostomes and gnathostomes) we can best identify the earliest major differences in development between vertebrate classes, and start to address how these might translate onto morphology. In this review we survey what is currently known about the cell biology and gene expression during head development in modern vertebrates, allowing us to better characterize the developmental genetics driving head skeleton formation in the most recent common ancestor of all living vertebrates. By pairing this vertebrate composite with information from fossil chordates, we can also deduce how gene regulatory modules might have been arranged in the ancestral vertebrate head. Together, we can immediately begin to understand which aspects of head skeletal development are the most conserved, and which are divergent, informing us as to when the first differences appear during development, and thus which pathways or cell types might be involved in generating lineage specific shape and structure.

Embryonic expression of endothelins and their receptors in lamprey and frog reveals stem vertebrate origins of complex Endothelin signaling.

Neural crest cells (NCCs) are highly patterned embryonic cells that migrate along stereotyped routes to give rise to a diverse array of adult tissues and cell types. Modern NCCs are thought to have evolved from migratory neural precursors with limited developmental potential and patterning. How this occurred is poorly understood. Endothelin signaling regulates several aspects of NCC development, including their migration, differentiation, and patterning. In jawed vertebrates, Endothelin signaling involves multiple functionally distinct ligands (Edns) and receptors (Ednrs) expressed in various NCC subpopulations. To test the potential role of endothelin signaling diversification in the evolution of modern, highly patterned NCC, we analyzed the expression of the complete set of endothelin ligands and receptors in the jawless vertebrate, the sea lamprey (Petromyzon marinus). To better understand ancestral features of gnathostome edn and ednr expression, we also analyzed all known Endothelin signaling components in the African clawed frog (Xenopus laevis). We found that the sea lamprey has a gnathsotome-like complement of edn and ednr duplicates, and these genes are expressed in patterns highly reminiscent of their gnathostome counterparts. Our results suggest that the duplication and specialization of vertebrate Endothelin signaling coincided with the appearance of highly patterned and multipotent NCCs in stem vertebrates.

CRISPR/Cas9-mediated mutagenesis in the sea lamprey Petromyzon marinus: a powerful tool for understanding ancestral gene functions in vertebrates.

Lamprey is one of only two living jawless vertebrates, a group that includes the first vertebrates. Comparisons between lamprey and jawed vertebrates have yielded important insights into the origin and evolution of vertebrate physiology, morphology and development. Despite its key phylogenetic position, studies of lamprey have been limited by their complex life history, which makes traditional genetic approaches impossible. The CRISPR/Cas9 system is a bacterial defense mechanism that was recently adapted to achieve high-efficiency targeted mutagenesis in eukaryotes. Here we report CRISPR/Cas9-mediated disruption of the genes Tyrosinase and FGF8/17/18 in the sea lamprey Petromyzon marinus, and detail optimized parameters for producing mutant F0 embryos. Using phenotype and genotype analyses, we show that CRISPR/Cas9 is highly effective in the sea lamprey, with a majority of injected embryos developing into complete or partial mutants. The ability to create large numbers of mutant embryos without inbred lines opens exciting new possibilities for studying development in lamprey and other non-traditional model organisms with life histories that prohibit the generation of mutant lines.

A gene expression map of the larval Xenopus laevis head reveals developmental changes underlying the evolution of new skeletal elements.

The morphology of the vertebrate head skeleton is highly plastic, with the number, size, shape, and position of its components varying dramatically between groups. While this evolutionary flexibility has been key to vertebrate success, its developmental and genetic bases are poorly understood. The larval head skeleton of the frog Xenopus laevis possesses a unique combination of ancestral tetrapod features and anuran-specific novelties. We built a detailed gene expression map of the head mesenchyme in X. laevis during early larval development, focusing on transcription factor families with known functions in vertebrate head skeleton development. This map was then compared to homologous gene expression in zebrafish, mouse, and shark embryos to identify conserved and evolutionarily flexible aspects of vertebrate head skeleton development. While we observed broad conservation of gene expression between X. laevis and other gnathostomes, we also identified several divergent features that correlate to lineage-specific novelties. We noted a conspicuous change in dlx1/2 and emx2 expression in the second pharyngeal arch, presaging the differentiation of the reduced dorsal hyoid arch skeletal element typical of modern anamniote tetrapods. In the first pharyngeal arch we observed a shift in the expression of the joint inhibitor barx1, and new expression of the joint marker gdf5, shortly before skeletal differentiation. This suggests that the anuran-specific infrarostral cartilage evolved by partitioning of Meckel's cartilage with a new paired joint. Taken together, these comparisons support a model in which early patterning mechanisms divide the vertebrate head mesenchyme into a highly conserved set of skeletal precursor populations. While subtle changes in this early patterning system can affect skeletal element size, they do not appear to underlie the evolution of new joints or cartilages. In contrast, later expression of the genes that regulate skeletal element differentiation can be clearly linked to the evolution of novel skeletal elements. We posit that changes in the expression of downstream regulators of skeletal differentiation, like barx1 and gdf5, is one mechanism by which head skeletal element number and articulation are altered during evolution.

Evolution of the new vertebrate head by co-option of an ancient chordate skeletal tissue.

A defining feature of vertebrates (craniates) is a pronounced head that is supported and protected by a robust cellular endoskeleton. In the first vertebrates, this skeleton probably consisted of collagenous cellular cartilage, which forms the embryonic skeleton of all vertebrates and the adult skeleton of modern jawless and cartilaginous fish. In the head, most cellular cartilage is derived from a migratory cell population called the neural crest, which arises from the edges of the central nervous system. Because collagenous cellular cartilage and neural crest cells have not been described in invertebrates, the appearance of cellular cartilage derived from neural crest cells is considered a turning point in vertebrate evolution. Here we show that a tissue with many of the defining features of vertebrate cellular cartilage transiently forms in the larvae of the invertebrate chordate Branchiostoma floridae (Florida amphioxus). We also present evidence that during evolution, a key regulator of vertebrate cartilage development, SoxE, gained new cis-regulatory sequences that subsequently directed its novel expression in neural crest cells. Together, these results suggest that the origin of the vertebrate head skeleton did not depend on the evolution of a new skeletal tissue, as is commonly thought, but on the spread of this tissue throughout the head. We further propose that the evolution of cis-regulatory elements near an ancient regulator of cartilage differentiation was a major factor in the evolution of the vertebrate head skeleton.

Roles for FGF in lamprey pharyngeal pouch formation and skeletogenesis highlight ancestral functions in the vertebrate head.

A defining feature of vertebrates (craniates) is a pronounced head supported and protected by a cellularized endoskeleton. In jawed vertebrates (gnathostomes), the head skeleton is made of rigid three-dimensional elements connected by joints. By contrast, the head skeleton of modern jawless vertebrates (agnathans) consists of thin rods of flexible cellular cartilage, a condition thought to reflect the ancestral vertebrate state. To better understand the origin and evolution of the gnathostome head skeleton, we have been analyzing head skeleton development in the agnathan, lamprey. The fibroblast growth factors FGF3 and FGF8 have various roles during head development in jawed vertebrates, including pharyngeal pouch morphogenesis, patterning of the oral skeleton and chondrogenesis. We isolated lamprey homologs of FGF3, FGF8 and FGF receptors and asked whether these functions are ancestral features of vertebrate development or gnathostome novelties. Using gene expression and pharmacological agents, we found that proper formation of the lamprey head skeleton requires two phases of FGF signaling: an early phase during which FGFs drive pharyngeal pouch formation, and a later phase when they directly regulate skeletal differentiation and patterning. In the context of gene expression and functional studies in gnathostomes, our results suggest that these roles for FGFs arose in the first vertebrates and that the evolution of the jaw and gnathostome cellular cartilage was driven by changes developmentally downstream from pharyngeal FGF signaling.

BMP, Wnt and FGF signals are integrated through evolutionarily conserved enhancers to achieve robust expression of Pax3 and Zic genes at the zebrafish neural plate border.

Neural crest cells generate a range of cells and tissues in the vertebrate head and trunk, including peripheral neurons, pigment cells, and cartilage. Neural crest cells arise from the edges of the nascent central nervous system, a domain called the neural plate border (NPB). NPB induction is known to involve the BMP, Wnt and FGF signaling pathways. However, little is known about how these signals are integrated to achieve temporally and spatially specific expression of genes in NPB cells. Furthermore, the timing and relative importance of these signals in NPB formation appears to differ between vertebrate species. Here, we use heat-shock overexpression and chemical inhibitors to determine whether, and when, BMP, Wnt and FGF signaling are needed for expression of the NPB specifiers pax3a and zic3 in zebrafish. We then identify four evolutionarily conserved enhancers from the pax3a and zic3 loci and test their response to BMP, Wnt and FGF perturbations. We find that all three signaling pathways are required during gastrulation for the proper expression of pax3a and zic3 in the zebrafish NPB. We also find that, although the expression patterns driven by the pax3a and zic3 enhancers largely overlap, they respond to different combinations of BMP, Wnt and FGF signals. Finally, we show that the combination of the two pax3a enhancers is less susceptible to signaling perturbations than either enhancer alone. Taken together, our results reveal how BMPs, FGFs and Wnts act cooperatively and redundantly through partially redundant enhancers to achieve robust, specific gene expression in the zebrafish NPB.