Shedding subspecies: The influence of genetics on reptile subspecies taxonomy

Range-Wide Phylogeography and Ecological Niche Modeling Provide Insights into the Evolutionary History of the Mongolian Racerunner (Eremias argus) in Northeast Asia

The Mongolian racerunner, Eremias argus, is a small lizard endemic to Northeast Asia that can serve as an excellent model for investigating how geography and past climate change have jointly influenced the evolution of biodiversity in this region. To elucidate the processes underlying its diversification and demography, we reconstructed the range-wide phylogeographic pattern and evolutionary trajectory, using phylogenetic, population genetic, landscape genetic, Bayesian phylogeographic reconstruction and ecological niche modeling approaches. Phylogenetic analyses of the mtDNA cyt b gene revealed eight lineages that were unbounded by geographic region. The genetic structure of E. argus was mainly determined by geographic distance. Divergence dating indicated that E. argus and E. brenchleyi diverged during the Mid-Pliocene Warm Period. E. argus was estimated to have coalesced at~0.4351 Ma (Marine Isotope Stage 19). Bayesian phylogeographic diffusion analysis revealed out-of-Inner Mongolia and rapid colonization events from the end of the Last Interglacial to the Last Glacial Maximum, which is consistent with the expanded suitable range of the Last Glacial Maximum. Pre-Last Glacial Maximum growth of population is presented for most lineages of E. argus. The Glacial Maximum contraction model and the previous multiple glacial refugia hypotheses are rejected. This may be due to an increase in the amount of climatically favorable habitats in Northeast Asia. Furthermore, E. argus barbouri most likely represents an invalid taxon. The present study is the first to report a range-wide phylogeography of reptiles over such a large region in Northeast Asia. Our results make a significant contribution towards understanding the biogeography of the entire Northeast Asia.

Acknowledging more biodiversity without more species

Delimiting and naming biodiversity is a vital step toward wildlife conservation and research. However, species delimitation must be consistent across biota so that the limited resources available for nature protection can be spent effectively and objectively. To date, newly discovered lineages typically are either left undescribed and thus remain unprotected or are being erroneously proposed as new species despite mixed evidence for completed speciation, in turn contributing to the emerging problem of taxonomic inflation. Inspired by recent conceptual and methodological progress, we propose a standardized workflow for species delimitation that combines phylogenetic and hybrid zone analyses of genomic datasets ("genomic taxonomy"), in which phylogeographic lineages that do not freely admix are ranked as species, while those that have remained fully genetically compatible are ranked as subspecies. In both cases, we encourage their formal taxonomic naming, diagnosis, and description to promote social awareness toward biodiversity. The use of loci throughout the genome overcomes the unreliability of widely used barcoding genes when phylogeographic patterns are complex, while the evaluation of divergence and reproductive isolation unifies the long-opposed concepts of lineage species and biological species. We suggest that a shift in conservation assessments from a single level (species) toward a two-level hierarchy (species and subspecies) will lead to a more balanced perception of biodiversity in which both intraspecific and interspecific diversity are valued and more adequately protected.

Super-high procoagulant activity of gecko thrombin: A gift from sky dragon

AIMS

Gecko, the "sky dragon" named by Traditional Chinese Medicine, undergoes rapid coagulation and scarless regeneration following tail amputation in the natural ecology, providing a perfect opportunity to develop the efficient and safe drug for blood clotting. Here, gecko thrombin (gthrombin) was recombinantly prepared and comparatively studied on its procoagulant activity.

METHODS

The 3D structure of gthrombin was constructed using the homology modeling method of I-TASSER. The active gthrombin was prepared by the expression of gecko prethrombin-2 in 293 T cells, followed by purification with Ni2+ -chelating column chromatography prior to activation by snake venom-derived Ecarin. The enzymatic activities of gthrombin were assayed by hydrolysis of synthetic substrate S-2238 and the fibrinogen clotting. The vulnerable nerve cells were used to evaluate the toxicity of gthrombin at molecular and cellular levels.

RESULTS

The active recombinant gthrombin showed super-high catalytic and fibrinogenolytic efficiency than those of human under different temperatures and pH conditions. In addition, gthrombin made nontoxic effects on the central nerve cells including neurons, contrary to those of mammalian counterparts, which contribute to neuronal damage, astrogliosis, and demyelination.

CONCLUSIONS

A super-high activity but safe procoagulant candidate drug was identified from reptiles, which provided a promising perspective for clinical application in rapid blood clotting.

Geography and past climate changes have shaped the evolution of a widespread lizard in arid Central Asia

The complex orogenic history and structure of Central Asia, coupled with Pleistocene glacial cycles have generated its stepwise aridification. Such events would have significantly influenced the evolution of many mid-latitude species in arid Central Asia (ACA). In this study, we employed two mitochondrial genes (CO1 and ND2) and genome-wide SNPs, coupled with ecological niche modeling, to investigate the lineage diversification and historical demography within a widespread lizard Phrynocepahlus helioscopus, and their associations with geography and past climate change. We obtained themtDNA dataset for 300 individuals from 96 localities within the known range of the lizard, among which 51 individuals from 27 localities were selected for generating the SNP dataset via genotyping-by-sequencing approach. Phylogenetic analyses of the concatenated mtDNA dataset revealed eight geographically correlated lineages that diverged by 4.21-10.41% for the CO1 gene, which were estimated to have coalesced ∼4.47 million years ago. However, we observed mito-nuclear discordance pattern regarding the population of Clade V (P. helioscopus sergeevi) from Tajikistan. Ancestral area estimations suggested that P. helioscopus originated from the Fergana Valley and then dispersed into the adjacent areas in ACA along with a history of multiple allopatric divergence processes, suggesting that Fergana may have been the cradle of diversification of P. helioscopus. The intensification of aridification across Central Asia during the Late Pliocene may have facilitated the rapid radiation of this arid-adapted lizard throughout this vast territory. Subsequently, the geological events (e.g., uplift of the Hissar-Alay, transgressions of the Caspian Sea) and geographic barriers (e.g., Amu Darya River, Zerarshan River) during the Pleistocene triggered the progressive diversification of P. helioscopus. Interestingly, Clade VIII (P. helioscopus varius) experienced rapid population growth coupled with range expansion while Clade IV (P. helioscopus cameranoi) underwent drastic population expansion associated with range contraction during the Last Glacial Maximum. In Clade IV, but not in Clade VIII, environmental turnover contributes more to mitochondrial genetic distinctiveness than geographic distance. Overall, the SNP dataset demonstrates that geographic distance plays a greater role than environmental distance. Both the mtDNA dataset and the SNP dataset suggest local-scale genetic differentiation in Clade IV and Clade VIII, revealing potential geographic barriers in the Ili River Valley and the Junggar Basin, respectively. Twenty-seven outlier SNPs associated with environmental factors (precipitation and temperature) were identified, which supports the signature of local adaptation to the arid desert environment. Finally, our finding suggests taxonomic implications, such as support for full species status for P. saidalievi (Clade II) and P. meridionalis (Clade I). Future analyses based on further evidence and increased taxon and geographic sampling should be carried out to corroborate our findings.

NA. The taxonomic status of the kukri snake Oligodon arenarius Vassilieva, 2015 with a redescription of Oligodon macrurus (Angel, 1927) (Squamata, Serpentes, Colubridae)

We investigated the taxonomic status of the recently described kukri snake Oligodon arenarius Vassilieva, 2015 and the morphologically similar Oligodon macrurus (Angel, 1927), two species endemic to the southern coast of Vietnam. Based on phylogenetic analyses using three mitochondrial genes (12S–16S rRNA, cytochrome b), we recovered O. arenarius and O. macrurus in a clade within the O. cyclurus-taeniatus species group, agreeing with previous intrageneric classifications. Genetic distances between O. arenarius and O. macrurus are extremely low (less than 0.5% based on 12S–16S) and render O. arenarius paraphyletic. All preserved specimens of O. arenarius and O. macrurus convey little to no differences in color pattern, hemipenial morphology and osteological features; the latter of which is based on three dimensional micro computer tomography (µCT) scans of one specimen per species. Contrasting these results, univariate and multivariate analyses revealed significant differences in relative tail length, and the number of ventral and subcaudal scales between both species. Although the molecular and morphological datasets present conflicting results, integrating the evidence leads us to synonymize O. arenarius with O. macrurus. We provide a formal redescription of O. macrurus, designate a neotype specimen to avoid future taxonomic confusion, and provide the first detailed osteological description of this species. Oligodon macrurus sensu stricto is endemic to coastal dunefields and adjacent forest habitats in southern Vietnam, where ongoing human development, tourism and road mortality pose significant threats to its conservation. Consequently, we suggest that O. macrurus should be listed as “Vulnerable” based on the assessment criteria of the International Union for Conservation of Nature (IUCN).

Empirical and philosophical problems with the subspecies rank

Species-level taxonomy derives from empirical sources (data and techniques) that assess the existence of spatiotemporal evolutionary lineages via various species "concepts." These concepts determine if observed lineages are independent given a particular methodology and ontology, which relates the metaphysical species concept to what "kind" of thing a species is in reality. Often, species concepts fail to link epistemology back to ontology. This lack of coherence is in part responsible for the persistence of the subspecies rank, which in modern usage often functions as a placeholder between the evolutionary events of divergence or collapse of incipient species. Thus, prospective events like lineages merging or diverging require information from unknowable future information. This is also conditioned on evidence that the lineage already has a detectably distinct evolutionary history. Ranking these lineages as subspecies can seem attractive given that many lineages do not exhibit intrinsic reproductive isolation. We argue that using subspecies is indefensible on philosophical and empirical grounds. Ontologically, the rank of subspecies is either identical to that of species or undefined in the context of evolutionary lineages representing spatiotemporally defined individuals. Some species concepts more inclined to consider subspecies, like the Biological Species Concept, are disconnected from evolutionary ontology and do not consider genealogy. Even if ontology is ignored, methods addressing reproductive isolation are often indirect and fail to capture the range of scenarios linking gene flow to species identity over space and time. The use of subspecies and reliance on reproductive isolation as a basis for an operational species concept can also conflict with ethical issues governing the protection of species. We provide a way forward for recognizing and naming species that links theoretical and operational species concepts regardless of the magnitude of reproductive isolation.

Split it up and see: using proxies to highlight divergent inter-populational performances in aquaculture standardised conditions

Background

Considering wild inter-populational phenotypic differentiation can facilitate domestication and subsequent production of new species. However, comparing all populations across a species range to identify those exhibiting suitable key traits for aquaculture (KTA; i.e. important for domestication and subsequent production) expressions is not feasible. Therefore, proxies highlighting inter-populational divergences in KTA are needed. The use of such proxies would allow to identify, prior to bioassays, the wild population pairs which are likely to present differentiations in KTA expressions in aquaculture conditions. Here, we assessed the relevance of three alternative proxies: (i) genetic distance, (ii) habitat divergence, and (iii) geographic/hydrologic distances. We performed this evaluation on seven allopatric populations of Perca fluviatilis for which divergences in KTA had already been shown.

Results

We showed differences in the correlation degree between the alternative proxy-based and KTA-based distance matrices, with the genetic proxy being correlated to the highest number of KTA. However, no proxy was correlated to all inter-populational divergences in KTA.

Conclusion

For future domestication trials, we suggest using a multi-proxy assessment along with a prioritisation strategy to identify population pairs which are of interest for further evaluation in bioassays.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01937-z.

Cryptic diversity and phylogeography of the Rhabdophis nuchalis group (Squamata: Colubridae)

Previous studies, have found that the rapid uplift of the Tibetan plateau accelerated the diversification of species. However, there are few relevant biogeographic data for the Colubridae in this region. We conducted a comprehensive study of the Rhabdophis nuchalis Group, which presently contains four nominal species, R. nuchalis, R. pentasupralabialis, R. leonardi, and R. chiwen. Building upon previous studies with specimens we have recently examined, greater interspecific and intraspecific diversity has been revealed. Here we address three questions: (1) Do the intraspecific differences represent only geographic variation within lineages, or are there cryptic species? (2) What are the interspecific relationships among members of the R. nuchalis Group? (3) What has been the biogeographic history of this species group? To resolve these questions we used four mitochondrial gene sequences and one nuclear sequence to investigate the molecular phylogenetic and geographic relationships among populations. Our molecular analysis reveals cryptic species diversity within the R. nuchalis Group, and seven clades were identified in the analysis. Ancestral area estimation suggests that the R. nuchalis Group originated in the Hengduan Mountains approximately 6.24 Mya and expanded its range northward to the Qinling-Daba Mountains. The Sichuan Basin appears to have been a barrier to migration. Species divergence seems to have been related to the rapid uplift of the Qinghai-Tibet Plateau.

Quaternary range dynamics and taxonomy of the Mediterranean collared dwarf racer, Platyceps collaris (Squamata: Colubridae)

The geological and geographical settings of the Eastern Mediterranean have resulted in complex patterns of intraspecific diversifications and phylogeographical histories that can be observed in squamates. In this study, we examined genetic differentiation of the Collared dwarf racer (Platyceps collaris) using a multilocus genetic dataset with a sampling that covered the entire range of the species. We developed distribution models in current and past climatic conditions to assess the dynamics of the species distribution through time. We sequenced a fragment of the cytochrome b mitochondrial gene of the holotype and eight paratypes of Coluber rubriceps thracius, which is considered a synonym of Platyceps collaris. Our results show that there are two distinct clades within P. collaris, one occupying the Balkans and western and southern Anatolia (termed the Balkan–Anatolian clade), the other in the Levant (termed the Levantine clade). All type specimens of C. r. thracius are genetically identical and cluster within the Balkan–Anatolian clade. Distribution models indicate the presence of two refugia during climatically challenging periods. One was in western Anatolia and served as a source for the colonization of the Balkans and southern Anatolia, and the other was in the northern Levant, from where P. collaris dispersed further south. According to our results, we revise the subspecific taxonomy of P. collaris.

Painted black: Iguana melanoderma (Reptilia, Squamata, Iguanidae) a new melanistic endemic species from Saba and Montserrat islands (Lesser Antilles)

The Lesser Antilles, in the Eastern Caribbean, is inhabited by three Iguana species: the Lesser Antillean iguanaIguana delicatissima, which is endemic to the northernmost islands of the Lesser Antilles, the introduced common iguana from South America, Iguana iguana iguana, represented also by the two newly described endemic subspecies Iguana iguana sanctaluciae from Saint Lucia and Iguana iguana insularis from Saint Vincent and the Grenadines, and Grenada, and the introduced Iguana rhinolopha from Central America. Drawing on both morphological and genetic data, this paper describes the Iguana populations from Saba and Montserrat as a new species, Iguana melanoderma. This species is recognized on the basis of the following combination of characteristics: private microsatellite alleles, unique mitochondrial ND4 haplotypes, a distinctive black spot between the eye and tympanum, a dorsal carpet pattern on juveniles and young adults, a darkening of body coloration with aging (except for the anterior part of the snout), a black dewlap, pink on the jowl, the high number of large tubercular nape scales, fewer than ten medium sized-triangular dewlap spikes, high dorsal spikes, and lack of horns on the snout. This new melanistic taxon is threatened by unsustainable harvesting (including for the pet trade) and both competition and hybridization from escaped or released invasive alien iguanas (I. iguana iguana and I. rhinolopha) from South and Central America, respectively. The authors call for action to conserve Iguana melanoderma in Saba and Montserrat and for further research to investigate its relationship to other melanistic iguanas from the Virgin Islands and coastal islands of Venezuela.

Genomic analyses reveal two species of the matamata (Testudines: Chelidae: Chelus spp.) and clarify their phylogeography

The matamata is one of the most charismatic turtles on earth, widely distributed in northern South America. Debates have occurred over whether or not there should be two subspecies or species recognized due to its geographic variation in morphology. Even though the matamata is universally known, its natural history, conservation status and biogeography are largely unexplored. In this study we examined the phylogeographic differentiation of the matamata based on three mitochondrial DNA fragments (2168 bp of the control region, cytochrome oxidase subunit I, and the cytochrome b gene), one nuclear genomic DNA fragment (1068 bp of the R35 intron) and 1661 Single Nucleotide Polymorphisms (SNPs). Our molecular and morphological analyses revealed the existence of two distinct, genetically deeply divergent evolutionary lineages of matamatas that separated in the late Miocene (approximately 12.7 million years ago), corresponding well to the time when the Orinoco Basin was established. As a result of our analyses, we describe the genetically and morphologically highly distinct matamata from the Orinoco and Río Negro Basins and the Essequibo drainage as a species new to science (Chelus orinocensis sp. nov.). Chelus fimbriata sensu stricto is distributed in the Amazon Basin and the Mahury drainage. Additionally, the analyses revealed that each species displays phylogeographic differentiation. For C. orinocensis, there is moderate mitochondrial differentiation between the Orinoco and the Río Negro. For C. fimbriata, there is more pronounced differentiation matching different river systems. One mitochondrial clade was identified from the Amazon, Ucayali, and Mahury Rivers, and another one from the Madeira and Jaci Paraná Rivers. The C. orinocensis in the Essequibo and Branco Rivers have haplotypes that constitute a third clade clustering with C. fimbriata. Phylogenetic analyses of the R35 intron and SNP data link the matamatas from the Essequibo and Branco with the new species, suggesting past gene flow and old mitochondrial introgression. Chelus orinocensis is collected for the pet trade in Colombia and Venezuela. However, neither the extent of the harvest nor its impact are known. Hence, it is crucial to gather more information and to assess its exploitation throughout its distribution range to obtain a better understanding of its conservation status and to design appropriate conservation and management procedures. RESUMEN: La matamata es una de las tortugas más carismáticas del mundo, ampliamente distribuida en el norte de Sudamérica. Debido a su variación morfológica geográfica, se debate sobre el reconocimiento de dos subespecies o especies. A pesar de que la matamata es universalmente conocida, su historia natural, estado de conservación y biogeografía han sido muy poco estudiados. En este estudio examinamos la diferenciación filogeográfica de las matamatas en base ​​a tres fragmentos de ADN mitocondrial (2168 pb de la región de control, la subunidad I del citocromo oxidasa y el gen del citocromo b), un fragmento de ADN genómico nuclear (1068 pb del intrón R35) y 1661 polimorfismos de nucleótido único (SNPs). Nuestros análisis moleculares y morfológicos revelaron la existencia de dos linajes evolutivos distintos de matamatas, genéticamente divergentes que se separaron en el Mioceno tardio (hace aproximadamente 12.7 millones de años), correspondiendo al tiempo en que se estableció la cuenca del Orinoco. Como resultado de nuestros análisis, describimos las genéticamente y morfológicamente distintas matamatas de las cuencas del Orinoco, Río Negro y Essequibo como una especie nueva para la ciencia (Chelus orinocensis sp. nov.). Chelus fimbriata sensu stricto se distribuye en la cuenca del Amazonas y en el drenaje del Mahury. Adicionalmente, los análisis revelaron que cada especie muestra diferenciación filogeográfica. Para C. orinocensis, hay una moderada diferenciación mitocondrial entre el Orinoco y el Río Negro. Para C. fimbriata, hay una diferenciación más pronunciada, concordando con los diferentes sistemas fluviales. Se identificó un clado de los ríos Amazonas, Ucayali y Mahury y otro de los ríos Madeira y Jaci Paraná. Las C. orinocensis de los ríos Essequibo y Branco tienen haplotipos que constituyen un tercer clado que se agrupa con C. fimbriata. Los análisis filogenéticos del intrón R35 y los datos de SNP asocian las matamatas de Essequibo y Branco con la nueva especie, sugiriendo flujo de genes pasado ​​e introgresión mitocondrial antigua. Chelus orinocensis se colecta para el comercio de mascotas en Colombia y Venezuela. Sin embargo, ni se conoce el alcance de las colectas ni su impacto. Por lo tanto, es crucial recopilar más información y evaluar su explotación en todo su rango de distribución, comprender mejor su estado de conservación y para diseñar acciones apropiadas de conservación y manejo.

Conservation Below the Species Level: Suitable Evolutionarily Significant Units among Mountain Vipers (the Montivipera raddei complex) in Iran

Northern and western mountains of Iran are among the most important biodiversity and endemism hot spots for reptiles in the Middle East. Among herpetofauna, the montivipers represent an emblematic and fragmented endemic group for which estimating their level of genetic differentiation and defining conservation priorities is urgently needed. Here, we present the most comprehensive phylogenetic study on the Montivipera raddei species group comprising all 5 known taxa, among which 3 are endemic to Iran. Based on 2 mitochondrial genes, phylogenetic and phylogeographic analyses revealed 3 major lineages each presenting very contrasting distribution areas. The Iranian montivipers are highly structured in clades showing low genetic diversity and corresponding to high altitude summits. Molecular dating revealed the role of Quaternary paleo-climatic oscillations and altitudinal movements of montivipers in shaping genetic diversity and differentiation of these sky-island taxa. In addition, the best scenario of historical biogeography allowed identifying 3 possible refugial areas in Iran most likely arising by vicariance. Based on our mitochondrial results and pending additional data, we recognize 3 candidate species among the M. raddei complex: M. raddei, Montivipera latifii, and Montivipera kuhrangica that are coherent with their geographical distribution. We propose that the most appropriate evolutionary significant units for conservation of the montivipers are represented by 13 units among which 6 are recognized as high priority. Finally, we suggest some recommendations to the IUCN as well as to the Iranian conservation policies with respect to conservation prioritization.

Subspecies definitions and legislation: from eastern wallaroo (Osphranter robustus robustus) to euro (Osphranter robustus erubescens)

As a consequence of genetic studies of population structuring, the usefulness of subspecies has been questioned, with opinions divided. The situation is further confused by the use of varying species and subspecies concepts. Most alternatives require each taxon to be following an independent evolutionary trajectory. These include traditional approaches and the more recent phylogenetic species concept. The latter has led to large increases in the apparent number of taxa in some groups, though strong objections have been raised to this approach. An alternative, the ecological species concept, has been opposed by phylogeneticists. These two approaches are compared using morphological and genetic data from common wallaroo (Osphranter robustus) populations. The different taxonomies that might result (many species, one species with two subspecies, one species with no subspecies) can have significant consequences for legislative and management decisions. The ecological approach is considered preferable for subspecies and the present taxonomy is maintained. A potential location of the boundary between the wallaroo subspecies is proposed. How the use of the different subspecies definitions would affect legislative decisions is explored.

Mitochondrial Divergence between Western and Eastern Great Bustards: Implications for Conservation and Species Status

The great bustard is the heaviest bird capable of flight and an iconic species of the Eurasian steppe. Populations of both currently recognized subspecies are highly fragmented and critically small in Asia. We used DNA sequence data from the mitochondrial cytochrome b gene and the mitochondrial control region to estimate the degree of mitochondrial differentiation and rates of female gene flow between the subspecies. We obtained genetic samples from 51 individuals of Otis tarda dybowskii representing multiple populations, including the first samples from Kazakhstan and Mongolia and samples from near the Altai Mountains, the proposed geographic divide between the subspecies, allowing for better characterization of the boundary between the 2 subspecies. We compared these with existing sequence data (n = 66) from Otis tarda tarda. Our results suggest, though do not conclusively prove, that O. t. dybowskii and O. t. tarda may be distinct species. The geographic distribution of haplotypes, phylogenetic analysis, analyses of molecular variance, and coalescent estimation of divergence time and female migration rates indicate that O. t. tarda and O. t. dybowskii are highly differentiated in the mitochondrial genome, have been isolated for approximately 1.4 million years, and exchange much less than 1 female migrant per generation. Our findings indicate that the 2 forms should at least be recognized and managed as separate evolutionary units. Populations in Xinjiang, China and Khövsgöl and Bulgan, Mongolia exhibited the highest levels of genetic diversity and should be prioritized in conservation planning.

Phylogeography and evolution of the Red Salamander (Pseudotriton ruber)

Phylogeographic studies frequently result in the elevation of subspecific taxa to species given monophyly, or the synonymy of subspecies that are not monophyletic. However, given limited or incongruent datasets, retention of subspecies can be useful to describe hypothesized incipient species or to illustrate interesting biological phenomena driving morphological diversity. Four subspecific taxa have been used to describe largely allopatric geographic variation within the species Pseudotriton ruber, a plethodontid salamander occupying stream and spring habitats across eastern North America: P. r. vioscai occurs in lowland Coastal Plain habitats, while P. r. ruber, P. r. nitidus, and P. r. schencki occupy upland regions in and around the Appalachian Mountains. Pseudotriton ruber co-occurs through its distribution with the aposematic newt Notophthalmus viridescens, and both species are hypothesized to be part of a Müllerian mimicry complex. In this study, we sequenced regions of two mitochondrial (cytochrome b, NADH dehydrogenase subunit 2) and one single copy nuclear protein-coding gene (pro-opiomelanocortin) from individuals sampled across much of the distribution of P. ruber and then used maximum-likelihood and Bayesian phylogenetic inference to test the monophyly of subspecies, reconstruct biogeographic history, and make inferences about morphological evolution. Phylogeographic hypotheses from mitochondrial and nuclear datasets described structure among populations of P. ruber which separated Coastal Plain and upland Appalachian populations, but subspecies were not monophyletic. Biogeographic reconstruction estimated the ancestor of all populations to have occupied and initially diverged in the Coastal Plain during the Pliocene (∼3.6mya), before one lineage subsequently invaded upland areas of Appalachia. Bold bright coloration of high elevation subspecies P. r. nitidus and P. r. schencki appears to have evolved twice. We hypothesize that the Müllerian mimicry complex with N. viridescens and P. ruber may provide a selective mechanism driving the co-evolution of striking bright and dull morphological variation among populations of both species. While P. ruber subspecies were not consistent with our criteria for diagnosing species (monophyly) and therefore could not be elevated to species, we advocate for the retention of subspecies because they describe hypotheses about an incipient species (P. r. vioscai) and how Müllerian mimicry may shape morphological diversity of species.

Gekko japonicus genome reveals evolution of adhesive toe pads and tail regeneration

Reptiles are the most morphologically and physiologically diverse tetrapods, and have undergone 300 million years of adaptive evolution. Within the reptilian tetrapods, geckos possess several interesting features, including the ability to regenerate autotomized tails and to climb on smooth surfaces. Here we sequence the genome of Gekko japonicus (Schlegel's Japanese Gecko) and investigate genetic elements related to its physiology. We obtain a draft G. japonicus genome sequence of 2.55 Gb and annotated 22,487 genes. Comparative genomic analysis reveals specific gene family expansions or reductions that are associated with the formation of adhesive setae, nocturnal vision and tail regeneration, as well as the diversification of olfactory sensation. The obtained genomic data provide robust genetic evidence of adaptive evolution in reptiles.

Geckos are small, agile reptiles with nocturnal habits. Here, the authors sequence the genome of the Schlegel's Japanese Gecko and reveal gene family expansions and reductions associated with formation of adhesive setae, nocturnal vision, tail regeneration, and diversification of olfactory sensation.