Comparative proteomic analysis of tail regeneration in the green anole lizard, Anolis carolinensis.

As amniote vertebrates, lizards are the most closely related organisms to humans capable of appendage regeneration. Lizards can autotomize, or release their tails as a means of predator evasion, and subsequently regenerate a functional replacement. Green anoles (Anolis carolinensis) can regenerate their tails through a process that involves differential expression of hundreds of genes, which has previously been analyzed by transcriptomic and microRNA analysis. To investigate protein expression in regenerating tissue, we performed whole proteomic analysis of regenerating tail tip and base. This is the first proteomic data set available for any anole lizard. We identified a total of 2,646 proteins - 976 proteins only in the regenerating tail base, 796 only in the tail tip, and 874 in both tip and base. For over 90% of these proteins in these tissues, we were able to assign a clear orthology to gene models in either the Ensembl or NCBI databases. For 13 proteins in the tail base, 9 proteins in the tail tip, and 10 proteins in both regions, the gene model in Ensembl and NCBI matched an uncharacterized protein, confirming that these predictions are present in the proteome. Ontology and pathways analysis of proteins expressed in the regenerating tail base identified categories including actin filament-based process, ncRNA metabolism, regulation of phosphatase activity, small GTPase mediated signal transduction, and cellular component organization or biogenesis. Analysis of proteins expressed in the tail tip identified categories including regulation of organelle organization, regulation of protein localization, ubiquitin-dependent protein catabolism, small GTPase mediated signal transduction, morphogenesis of epithelium, and regulation of biological quality. These proteomic findings confirm pathways and gene families activated in tail regeneration in the green anole as well as identify uncharacterized proteins whose role in regrowth remains to be revealed. This study demonstrates the insights that are possible from the integration of proteomic and transcriptomic data in tail regrowth in the green anole, with potentially broader application to studies in other regenerative models.

A chromosome-level genome assembly for the eastern fence lizard (Sceloporus undulatus), a reptile model for physiological and evolutionary ecology.

BACKGROUND: High-quality genomic resources facilitate investigations into behavioral ecology, morphological and physiological adaptations, and the evolution of genomic architecture. Lizards in the genus Sceloporus have a long history as important ecological, evolutionary, and physiological models, making them a valuable target for the development of genomic resources. FINDINGS: We present a high-quality chromosome-level reference genome assembly, SceUnd1.0 (using 10X Genomics Chromium, HiC, and Pacific Biosciences data), and tissue/developmental stage transcriptomes for the eastern fence lizard, Sceloporus undulatus. We performed synteny analysis with other snake and lizard assemblies to identify broad patterns of chromosome evolution including the fusion of micro- and macrochromosomes. We also used this new assembly to provide improved reference-based genome assemblies for 34 additional Sceloporus species. Finally, we used RNAseq and whole-genome resequencing data to compare 3 assemblies, each representing an increased level of cost and effort: Supernova Assembly with data from 10X Genomics Chromium, HiRise Assembly that added data from HiC, and PBJelly Assembly that added data from Pacific Biosciences sequencing. We found that the Supernova Assembly contained the full genome and was a suitable reference for RNAseq and single-nucleotide polymorphism calling, but the chromosome-level scaffolds provided by the addition of HiC data allowed synteny and whole-genome association mapping analyses. The subsequent addition of PacBio data doubled the contig N50 but provided negligible gains in scaffold length. CONCLUSIONS: These new genomic resources provide valuable tools for advanced molecular analysis of an organism that has become a model in physiology and evolutionary ecology.

Brain transcriptomic responses of Yarrow's spiny lizard, Sceloporus jarrovii, to conspecific visual or chemical signals.

Species with multimodal communication integrate information from social cues in different modalities into behavioral responses that are mediated by changes in gene expression in the brain. Differences in patterns of gene expression between signal modalities may shed light on the neuromolecular mechanisms underlying multisensory processing. Here, we use RNA-Seq to analyze brain transcriptome responses to either chemical or visual social signals in a territorial lizard with multimodal communication. Using an intruder challenge paradigm, we exposed 18 wild-caught, adult, male Sceloporus jarrovii to either male conspecific scents (femoral gland secretions placed on a small pebble), the species-specific push-up display (a programmed robotic model), or a control (an unscented pebble). We conducted differential expression analysis with both a de novo S. jarrovii transcriptome assembly and the reference genome of a closely related species, Sceloporus undulatus. Despite some inter-individual variation, we found significant differences in gene expression in the brain across signal modalities and the control in both analyses. The most notable differences occurred between chemical and visual stimulus treatments, closely followed by visual stimulus versus the control. Altered expression profiles could explain documented aggression differences in the immediate behavioral response to conspecific signals from different sensory modalities. Shared differentially expressed genes between visually- or chemically-stimulated males are involved in neural activity and neurodevelopment and several other differentially expressed genes in stimulus-challenged males are involved in conserved signal-transduction pathways associated with the social stress response, aggression and the response to territory intruders across vertebrates.

Anatomical and histological analyses reveal that tail repair is coupled with regrowth in wild-caught, juvenile American alligators (Alligator mississippiensis).

Reptiles are the only amniotes that maintain the capacity to regenerate appendages. This study presents the first anatomical and histological evidence of tail repair with regrowth in an archosaur, the American alligator. The regrown alligator tails constituted approximately 6-18% of the total body length and were morphologically distinct from original tail segments. Gross dissection, radiographs, and magnetic resonance imaging revealed that caudal vertebrae were replaced by a ventrally-positioned, unsegmented endoskeleton. This contrasts with lepidosaurs, where the regenerated tail is radially organized around a central endoskeleton. Furthermore, the regrown alligator tail lacked skeletal muscle and instead consisted of fibrous connective tissue composed of type I and type III collagen fibers. The overproduction of connective tissue shares features with mammalian wound healing or fibrosis. The lack of skeletal muscle contrasts with lizards, but shares similarities with regenerated tails in the tuatara and regenerated limbs in Xenopus adult frogs, which have a cartilaginous endoskeleton surrounded by connective tissue, but lack skeletal muscle. Overall, this study of wild-caught, juvenile American alligator tails identifies a distinct pattern of wound repair in mammals while exhibiting features in common with regeneration in lepidosaurs and amphibia.

Immune and sex-biased gene expression in the threatened Mojave desert tortoise, Gopherus agassizii.

The immune system of ectotherms, particularly non-avian reptiles, remains poorly characterized regarding the genes involved in immune function, and their function in wild populations. We used RNA-Seq to explore the systemic response of Mojave desert tortoise (Gopherus agassizii) gene expression to three levels of Mycoplasma infection to better understand the host response to this bacterial pathogen. We found over an order of magnitude more genes differentially expressed between male and female tortoises (1,037 genes) than differentially expressed among immune groups (40 genes). There were 8 genes differentially expressed among both variables that can be considered sex-biased immune genes in this tortoise. Among experimental immune groups we find enriched GO biological processes for cysteine catabolism, regulation of type 1 interferon production, and regulation of cytokine production involved in immune response. Sex-biased transcription involves iron ion transport, iron ion homeostasis, and regulation of interferon-beta production to be enriched. More detailed work is needed to assess the seasonal response of the candidate genes found here. How seasonal fluctuation of testosterone and corticosterone modulate the immunosuppression of males and their susceptibility to Mycoplasma infection also warrants further investigation, as well as the importance of iron in the immune function and sex-biased differences of this species. Finally, future transcriptional studies should avoid drawing blood from tortoises via subcarapacial venipuncture as the variable aspiration of lymphatic fluid will confound the differential expression of genes.

Discovery of a New TLR Gene and Gene Expansion Event through Improved Desert Tortoise Genome Assembly with Chromosome-Scale Scaffolds.

Toll-like receptors (TLRs) are a complex family of innate immune genes that are well characterized in mammals and birds but less well understood in nonavian sauropsids (reptiles). The advent of highly contiguous draft genomes of nonmodel organisms enables study of such gene families through analysis of synteny and sequence identity. Here, we analyze TLR genes from the genomes of 22 tetrapod species. Findings reveal a TLR8 gene expansion in crocodilians and turtles (TLR8B), and a second duplication (TLR8C) specifically within turtles, followed by pseudogenization of that gene in the nonfreshwater species (desert tortoise and green sea turtle). Additionally, the Mojave desert tortoise (Gopherus agassizii) has a stop codon in TLR8B (TLR8-1) that is polymorphic among conspecifics. Revised orthology further reveals a new TLR homolog, TLR21-like, which is exclusive to lizards, snakes, turtles, and crocodilians. These analyses were made possible by a new draft genome assembly of the desert tortoise (gopAga2.0), which used chromatin-based assembly to yield draft chromosomal scaffolds (L50 = 26 scaffolds, N50 = 28.36 Mb, longest scaffold = 107 Mb) and an enhanced de novo genome annotation with 25,469 genes. Our three-step approach to orthology curation and comparative analysis of TLR genes shows what new insights are possible using genome assemblies with chromosome-scale scaffolds that permit integration of synteny conservation data.

Virus Discovery in Desert Tortoise Fecal Samples: Novel Circular Single-Stranded DNA Viruses.

The Sonoran Desert tortoise Gopherus morafkai is adapted to the desert, and plays an important ecological role in this environment. There is limited information on the viral diversity associated with tortoises (family Testudinidae), and to date no DNA virus has been identified associated with these animals. This study aimed to assess the diversity of DNA viruses associated with the Sonoran Desert tortoise by sampling their fecal matter. A viral metagenomics approach was used to identify the DNA viruses in fecal samples from wild Sonoran Desert tortoises in Arizona, USA. In total, 156 novel single-stranded DNA viruses were identified from 40 fecal samples. Those belonged to two known viral families, the Genomoviridae (n = 27) and Microviridae (n = 119). In addition, 10 genomes were recovered that belong to the unclassified group of circular-replication associated protein encoding single-stranded (CRESS) DNA virus and five circular molecules encoding viral-like proteins.

Anolis carolinensis as a model to understand the molecular and cellular basis of foveal development.

The fovea is an anatomical specialization of the central retina containing closely packed cone-photoreceptors providing an area of high acuity vision in humans and primates. Despite its key role in the clarity of vision, little is known about the molecular and cellular basis of foveal development, due to the absence of a foveal structure in commonly used laboratory animal models. Of the amniotes the retina in birds of prey and some reptiles do exhibit a typical foveal structure, but they have not been studied in the context of foveal development due to lack of availability of embryonic tissue, lack of captive breeding programs, and limited genomic information. However, the genome for the diurnal bifoveate reptile species Anolis carolinensis (green anole) was recently published and it is possible to collect embryos from this species in captivity. Here, we tested the feasibility of using the anole as a model to study foveal development. Eyes were collected at various stages of development for histological analysis, immunofluorescence, and apoptosis. We show that at embryonic stage (ES) 10 there is peak ganglion cell density at the incipient central foveal region and a single row of cone photoreceptor nuclei. At ES17 the foveal pit begins to form and at this stage there are 3-4 rows of cone nuclei. Post-hatching a further increase in cone density and lengthening of inner and outer segments is observed. A yellowish pigment was seen in the adult central foveal region, but not in the temporal fovea. At ES14 Pax6 was localized across the entire retina, but was more prominent in the ganglion cell layer (GCL) and the part of the inner nuclear layer (INL) containing amacrine cell bodies. However, at ES17 Pax6 expression in the ganglion cells of the central retina was markedly reduced. Bioinformatic analysis revealed that 86% of human candidate foveal hypoplasia genes had an orthologous gene or DNA sequence in the green anole. These findings provide the first insight into foveal morphogenesis in the green anole and suggest that it could be a very useful model for investigating the molecular signals driving foveal development, and thus inform on human foveal development and disease.

Comparative Genomics Reveals Accelerated Evolution in Conserved Pathways during the Diversification of Anole Lizards.

Squamates include all lizards and snakes, and display some of the most diverse and extreme morphological adaptations among vertebrates. However, compared with birds and mammals, relatively few resources exist for comparative genomic analyses of squamates, hampering efforts to understand the molecular bases of phenotypic diversification in such a speciose clade. In particular, the ∼400 species of anole lizard represent an extensive squamate radiation. Here, we sequence and assemble the draft genomes of three anole species-Anolis frenatus, Anolis auratus, and Anolis apletophallus-for comparison with the available reference genome of Anolis carolinensis. Comparative analyses reveal a rapid background rate of molecular evolution consistent with a model of punctuated equilibrium, and strong purifying selection on functional genomic elements in anoles. We find evidence for accelerated evolution in genes involved in behavior, sensory perception, and reproduction, as well as in genes regulating limb bud development and hindlimb specification. Morphometric analyses of anole fore and hindlimbs corroborated these findings. We detect signatures of positive selection across several genes related to the development and regulation of the forebrain, hormones, and the iguanian lizard dewlap, suggesting molecular changes underlying behavioral adaptations known to reinforce species boundaries were a key component in the diversification of anole lizards.

Developmental and adult-specific processes contribute to de novo neuromuscular regeneration in the lizard tail.

Peripheral nerves exhibit robust regenerative capabilities in response to selective injury among amniotes, but the regeneration of entire muscle groups following volumetric muscle loss is limited in birds and mammals. In contrast, lizards possess the remarkable ability to regenerate extensive de novo muscle after tail loss. However, the mechanisms underlying reformation of the entire neuromuscular system in the regenerating lizard tail are not completely understood. We have tested whether the regeneration of the peripheral nerve and neuromuscular junctions (NMJs) recapitulate processes observed during normal neuromuscular development in the green anole, Anolis carolinensis. Our data confirm robust axonal outgrowth during early stages of tail regeneration and subsequent NMJ formation within weeks of autotomy. Interestingly, NMJs are overproduced as evidenced by a persistent increase in NMJ density 120 and 250 days post autotomy (DPA). Substantial Myelin Basic Protein (MBP) expression could also be detected along regenerating nerves indicating that the ability of Schwann cells to myelinate newly formed axons remained intact. Overall, our data suggest that the mechanism of de novo nerve and NMJ reformation parallel, in part, those observed during neuromuscular development. However, the prolonged increase in NMJ number and aberrant muscle differentiation hint at processes specific to the adult response. An examination of the coordinated exchange between peripheral nerves, Schwann cells, and newly synthesized muscle of the regenerating neuromuscular system may assist in the identification of candidate molecules that promote neuromuscular recovery in organisms incapable of a robust regenerative response.

Identification of satellite cells from anole lizard skeletal muscle and demonstration of expanded musculoskeletal potential.

The lizards are evolutionarily the closest vertebrates to humans that demonstrate the ability to regenerate entire appendages containing cartilage, muscle, skin, and nervous tissue. We previously isolated PAX7-positive cells from muscle of the green anole lizard, Anolis carolinensis, that can differentiate into multinucleated myotubes and express the muscle structural protein, myosin heavy chain. Studying gene expression in these satellite/progenitor cell populations from A. carolinensis can provide insight into the mechanisms regulating tissue regeneration. We generated a transcriptome from proliferating lizard myoprogenitor cells and compared them to transcriptomes from the mouse and human tissues from the ENCODE project using XGSA, a statistical method for cross-species gene set analysis. These analyses determined that the lizard progenitor cell transcriptome was most similar to mammalian satellite cells. Further examination of specific GO categories of genes demonstrated that among genes with the highest level of expression in lizard satellite cells were an increased number of genetic regulators of chondrogenesis, as compared to mouse satellite cells. In micromass culture, lizard PAX7-positive cells formed Alcian blue and collagen 2a1 positive nodules, without the addition of exogenous morphogens, unlike their mouse counterparts. Subsequent quantitative RT-PCR confirmed up-regulation of expression of chondrogenic regulatory genes in lizard cells, including bmp2, sox9, runx2, and cartilage specific structural genes, aggrecan and collagen 2a1. Taken together, these data suggest that tail regeneration in lizards involves significant alterations in gene regulation with expanded musculoskeletal potency.

Summary of the first inaugural joint meeting of the International Consortium for scoliosis genetics and the International Consortium for vertebral anomalies and scoliosis, March 16-18, 2017, Dallas, Texas.

Scoliosis represents the most common musculoskeletal disorder in children and affects approximately 3% of the world population. Scoliosis is separated into two major phenotypic classifications: congenital and idiopathic. Idiopathic scoliosis is defined as a curvature of the spine of 10° or greater visualized on plane radiograph and does not have associated vertebral malformations (VM). "Congenital" scoliosis (CS) due to malformations in vertebrae is frequently associated with other birth defects. Recently, significant advances have been made in understanding the genetic basis of both conditions. There is evidence that both conditions are etiologically related. A 2-day conference entitled "Genomic Approaches to Understanding and Treating Scoliosis" was held at Scottish Rite Hospital for Children in Dallas, Texas, to synergize research in this field. This first combined, multidisciplinary conference featured international scoliosis researchers in basic and clinical sciences. A major outcome of the conference advancing scoliosis research was the proposal and subsequent vote in favor of merging the International Consortium for Vertebral Anomalies and Scoliosis (ICVAS) and International Consortium for Scoliosis Genetics (ICSG) into a single entity called International Consortium for Spinal Genetics, Development, and Disease (ICSGDD). The ICSGDD is proposed to meet annually as a forum to synergize multidisciplinary spine deformity research.

Somitogenesis and Axial Development in Reptiles.

Among amniote vertebrates, reptiles display the greatest variation in axial skeleton morphology. Only recently have they been used in gene expression studies of somitogenesis , challenging previous assumptions about the segmentation clock and axial patterning. An increasing number of reptile genomes and transcriptomes are becoming available as next-generation sequencing becomes more affordable. Information regarding gene sequence and structure can be used to design and synthesize labeled riboprobes by in vitro transcription for gene expression analysis by in situ hybridization, thus, enabling the characterization of spatial and temporal expression patterns of genes involved in somitogenesis, a topic of great interest within evolutionary developmental studies of vertebrates.

The Agassiz's desert tortoise genome provides a resource for the conservation of a threatened species.

Agassiz's desert tortoise (Gopherus agassizii) is a long-lived species native to the Mojave Desert and is listed as threatened under the US Endangered Species Act. To aid conservation efforts for preserving the genetic diversity of this species, we generated a whole genome reference sequence with an annotation based on deep transcriptome sequences of adult skeletal muscle, lung, brain, and blood. The draft genome assembly for G. agassizii has a scaffold N50 length of 252 kbp and a total length of 2.4 Gbp. Genome annotation reveals 20,172 protein-coding genes in the G. agassizii assembly, and that gene structure is more similar to chicken than other turtles. We provide a series of comparative analyses demonstrating (1) that turtles are among the slowest-evolving genome-enabled reptiles, (2) amino acid changes in genes controlling desert tortoise traits such as shell development, longevity and osmoregulation, and (3) fixed variants across the Gopherus species complex in genes related to desert adaptations, including circadian rhythm and innate immune response. This G. agassizii genome reference and annotation is the first such resource for any tortoise, and will serve as a foundation for future analysis of the genetic basis of adaptations to the desert environment, allow for investigation into genomic factors affecting tortoise health, disease and longevity, and serve as a valuable resource for additional studies in this species complex.

Evolution of Dosage Compensation in Anolis carolinensis, a Reptile with XX/XY Chromosomal Sex Determination.

In species with highly heteromorphic sex chromosomes, the degradation of one of the sex chromosomes will result in unequal gene expression between the sexes (e.g. between XX females and XY males) and between the sex chromosomes and the autosomes. Dosage compensation is a process whereby genes on the sex chromosomes achieve equal gene expression. We compared genome-wide levels of transcription between males and females, and between the X chromosome and the autosomes in the green anole, Anolis carolinensis. We present evidence for dosage compensation between the sexes, and between the sex chromosomes and the autosomes. When dividing the X chromosome into regions based on linkage groups, we discovered that genes in the first reported X-linked region, anole linkage group b (LGb), exhibit complete dosage compensation, although the rest of the X-linked genes exhibit incomplete dosage compensation. Our data further suggest that the mechanism of this dosage compensation is upregulation of the X chromosome in males. We report that approximately 10% of coding genes, most of which are on the autosomes, are differentially expressed between males and females. In addition, genes on the X chromosome exhibited higher ratios of nonsynonymous to synonymous substitution than autosomal genes, consistent with the fast-X effect. Our results from the green anole add an additional observation of dosage compensation in a species with XX/XY sex determination.

XGSA: A statistical method for cross-species gene set analysis.

MOTIVATION: Gene set analysis is a powerful tool for determining whether an experimentally derived set of genes is statistically significantly enriched for genes in other pre-defined gene sets, such as known pathways, gene ontology terms, or other experimentally derived gene sets. Current gene set analysis methods do not facilitate comparing gene sets across different organisms as they do not explicitly deal with homology mapping between species. There lacks a systematic investigation about the effect of complex gene homology on cross-species gene set analysis. RESULTS: In this study, we show that not accounting for the complex homology structure when comparing gene sets in two species can lead to false positive discoveries, especially when comparing gene sets that have complex gene homology relationships. To overcome this bias, we propose a straightforward statistical approach, called XGSA, that explicitly takes the cross-species homology mapping into consideration when doing gene set analysis. Simulation experiments confirm that XGSA can avoid false positive discoveries, while maintaining good statistical power compared to other ad hoc approaches for cross-species gene set analysis. We further demonstrate the effectiveness of XGSA with two real-life case studies that aim to discover conserved or species-specific molecular pathways involved in social challenge and vertebrate appendage regeneration. AVAILABILITY AND IMPLEMENTATION: The R source code for XGSA is available under a GNU General Public License at http://github.com/VCCRI/XGSA CONTACT: jho@victorchang.edu.au.

Differential expression of conserved and novel microRNAs during tail regeneration in the lizard Anolis carolinensis.

BACKGROUND: Lizards are evolutionarily the most closely related vertebrates to humans that can lose and regrow an entire appendage. Regeneration in lizards involves differential expression of hundreds of genes that regulate wound healing, musculoskeletal development, hormonal response, and embryonic morphogenesis. While microRNAs are able to regulate large groups of genes, their role in lizard regeneration has not been investigated. RESULTS: MicroRNA sequencing of green anole lizard (Anolis carolinensis) regenerating tail and associated tissues revealed 350 putative novel and 196 known microRNA precursors. Eleven microRNAs were differentially expressed between the regenerating tail tip and base during maximum outgrowth (25 days post autotomy), including miR-133a, miR-133b, and miR-206, which have been reported to regulate regeneration and stem cell proliferation in other model systems. Three putative novel differentially expressed microRNAs were identified in the regenerating tail tip. CONCLUSIONS: Differentially expressed microRNAs were identified in the regenerating lizard tail, including known regulators of stem cell proliferation. The identification of 3 putative novel microRNAs suggests that regulatory networks, either conserved in vertebrates and previously uncharacterized or specific to lizards, are involved in regeneration. These findings suggest that differential regulation of microRNAs may play a role in coordinating the timing and expression of hundreds of genes involved in regeneration.

Assessing models of speciation under different biogeographic scenarios; an empirical study using multi-locus and RNA-seq analyses.

Evolutionary biology often seeks to decipher the drivers of speciation, and much debate persists over the relative importance of isolation and gene flow in the formation of new species. Genetic studies of closely related species can assess if gene flow was present during speciation, because signatures of past introgression often persist in the genome. We test hypotheses on which mechanisms of speciation drove diversity among three distinct lineages of desert tortoise in the genus Gopherus. These lineages offer a powerful system to study speciation, because different biogeographic patterns (physical vs. ecological segregation) are observed at opposing ends of their distributions. We use 82 samples collected from 38 sites, representing the entire species' distribution and generate sequence data for mtDNA and four nuclear loci. A multilocus phylogenetic analysis in *BEAST estimates the species tree. RNA-seq data yield 20,126 synonymous variants from 7665 contigs from two individuals of each of the three lineages. Analyses of these data using the demographic inference package ∂a∂i serve to test the null hypothesis of no gene flow during divergence. The best-fit demographic model for the three taxa is concordant with the *BEAST species tree, and the ∂a∂i analysis does not indicate gene flow among any of the three lineages during their divergence. These analyses suggest that divergence among the lineages occurred in the absence of gene flow and in this scenario the genetic signature of ecological isolation (parapatric model) cannot be differentiated from geographic isolation (allopatric model).

Transcriptomic analysis of tail regeneration in the lizard Anolis carolinensis reveals activation of conserved vertebrate developmental and repair mechanisms.

Lizards, which are amniote vertebrates like humans, are able to lose and regenerate a functional tail. Understanding the molecular basis of this process would advance regenerative approaches in amniotes, including humans. We have carried out the first transcriptomic analysis of tail regeneration in a lizard, the green anole Anolis carolinensis, which revealed 326 differentially expressed genes activating multiple developmental and repair mechanisms. Specifically, genes involved in wound response, hormonal regulation, musculoskeletal development, and the Wnt and MAPK/FGF pathways were differentially expressed along the regenerating tail axis. Furthermore, we identified 2 microRNA precursor families, 22 unclassified non-coding RNAs, and 3 novel protein-coding genes significantly enriched in the regenerating tail. However, high levels of progenitor/stem cell markers were not observed in any region of the regenerating tail. Furthermore, we observed multiple tissue-type specific clusters of proliferating cells along the regenerating tail, not localized to the tail tip. These findings predict a different mechanism of regeneration in the lizard than the blastema model described in the salamander and the zebrafish, which are anamniote vertebrates. Thus, lizard tail regrowth involves the activation of conserved developmental and wound response pathways, which are potential targets for regenerative medical therapies.

Reptile genomes open the frontier for comparative analysis of amniote development and regeneration.

Developmental genetic studies of vertebrates have focused primarily on zebrafish, frog and mouse models, which have clear application to medicine and well-developed genomic resources. In contrast, reptiles represent the most diverse amniote group, but have only recently begun to gather the attention of genome sequencing efforts. Extant reptilian groups last shared a common ancestor ?280 million years ago and include lepidosaurs, turtles and crocodilians. This phylogenetic diversity is reflected in great morphological and behavioral diversity capturing the attention of biologists interested in mechanisms regulating developmental processes such as somitogenesis and spinal patterning, regeneration, the evolution of "snake-like" morphology, the formation of the unique turtle shell, and the convergent evolution of the four-chambered heart shared by mammals and archosaurs. The complete genome of the first non-avian reptile, the green anole lizard, was published in 2011 and has provided insights into the origin and evolution of amniotes. Since then, the genomes of multiple snakes, turtles, and crocodilians have also been completed. Here we will review the current diversity of available reptile genomes, with an emphasis on their evolutionary relationships, and will highlight how these genomes have and will continue to facilitate research in developmental and regenerative biology.