Exploring Paleogene Tibet's warm temperate environments through target enrichment and phylogenetic niche modelling of Himalayan spiny frogs (Paini, Dicroglossidae).

The Cenozoic topographic development of the Himalaya-Tibet orogen (HTO) substantially affected the paleoenvironment and biodiversity patterns of High Asia. However, concepts on the evolution and paleoenvironmental history of the HTO differ massively in timing, elevational increase and sequence of surface uplift of the different elements of the orogen. Using target enrichment of a large set of transcriptome-derived markers, ancestral range estimation and paleoclimatic niche modelling, we assess a recently proposed concept of a warm temperate paleo-Tibet in Asian spiny frogs of the tribe Paini and reconstruct their historical biogeography. That concept was previously developed in invertebrates. Because of their early evolutionary origin, low dispersal capacity, high degree of local endemism, and strict dependence on temperature and humidity, the cladogenesis of spiny frogs may echo the evolution of the HTO paleoenvironment. We show that diversification of main lineages occurred during the early to Mid-Miocene, while the evolution of alpine taxa started during the late Miocene/early Pliocene. Our distribution and niche modelling results indicate range shifts and niche stability that may explain the modern disjunct distributions of spiny frogs. They probably maintained their (sub)tropical or (warm)temperate preferences and moved out of the ancestral paleo-Tibetan area into the Himalaya as the climate shifted, as opposed to adapting in situ. Based on ancestral range estimation, we assume the existence of low-elevation, climatically suitable corridors across paleo-Tibet during the Miocene along the Kunlun, Qiangtang and/or Gangdese Shan. Our results contribute to a deeper understanding of the mechanisms and processes of faunal evolution in the HTO.

Hidden hotspots of amphibian biodiversity in China.

Identifying and protecting hotspots of endemism and species richness is crucial for mitigating the global biodiversity crisis. However, our understanding of spatial diversity patterns is far from complete, which severely limits our ability to conserve biodiversity hotspots. Here, we report a comprehensive analysis of amphibian species diversity in China, one of the most species-rich countries on Earth. Our study combines 20 y of field surveys with new molecular analyses of 521 described species and also identifies 100 potential cryptic species. We identify 10 hotspots of amphibian diversity in China, each with exceptional species richness and endemism and with exceptional phylogenetic diversity and phylogenetic endemism (based on a new time-calibrated, species-level phylogeny for Chinese amphibians). These 10 hotspots encompass 59.6% of China's described amphibian species, 49.0% of cryptic species, and 55.6% of species endemic to China. Only four of these 10 hotspots correspond to previously recognized biodiversity hotspots. The six new hotspots include the Nanling Mountains and other mountain ranges in South China. Among the 186 species in the six new hotspots, only 9.7% are well covered by protected areas and most (88.2%) are exposed to high human impacts. Five of the six new hotspots are under very high human pressure and are in urgent need of protection. We also find that patterns of richness in cryptic species are significantly related to those in described species but are not identical.

Male spiny frogs enter the underwater battlefield with loose skin exhibiting enhanced penetration of capillaries into the epidermis.

The skin has multiple functions, and capillaries can penetrate the epidermis to shorten the diffusion path while allowing maintenance of overall epidermal thickness for nonrespiratory roles. However, a method for quantifying the capillary penetration extent is lacking. Such a method may facilitate making comparisons and detecting associations, potentially making the extent a useful variable in biological studies. We quantified the extent as the ratio of the average minimum thickness of epidermis overlying each capillary to the average epidermal thickness along a skin section and then explored its performance in the Emei mustache toad, Leptobrachium boringii, a species in which breeding males with loose skin call and fight each other with maxillary spines underwater. The ratio showed informative associations with other variables, such as perfused capillary density. It displayed small intragroup variation and could be more sensitive than other variables in revealing structural differences in the skin. The ratio estimates were lowest and were correlated with epidermal and stratum compactum thicknesses in breeding males, i.e., a covariation but not reinforcement against stabbing, constituting early evidence consistent with the increased extensibility of loose skin conferring a defensive advantage during combat in amphibians. In addition, our results lead to the hypothesis that high hemoglobin density along subepidermal capillaries favors the maintenance of low blood partial oxygen pressure and hence increases cutaneous oxygen uptake. We also provide evidence supporting the new idea that the cooccurrence of loose skin and underwater calling found in some frogs can be explained by the latter benefiting from a large functional respiratory surface area. Awareness of the usefulness of the ratio may promote its application and the quantification of the penetration. Regarding exchange surface design, these findings for L. boringii imply a case in which looseness increases surface area as well as prevents damage.

Sex-Linked Loci on the W Chromosome in the Multi-Ocellated Racerunner (Eremias multiocellata) Confirm Genetic Sex-Determination Stability in Lacertid Lizards.

The multi-ocellated racerunner, Eremias multiocellata, was considered to have temperature-dependent sex determination (TSD), as its sex ratio can be influenced at different temperatures. However, such an observation contrasts with recent findings that suggest TSD is less common than previously thought. Here, a genotyping-by-sequencing (GBS) approach was employed to identify sex-linked markers in the E. multiocellata, for which the mechanism choice of TSD or GSD is still controversial. We preliminarily identified 119 sex-linked markers based on sex-associated sex-specific sequences, 97% of which indicated female heterogamety. After eliminating the false positives, 38 sex-linked markers were recognized, all of which showed the ZW/ZZ system. Then, eight of the novel markers were verified by PCR amplification from 15 populations of E. multiocellata, which support the GSD in E. multiocellata without geographic variation. To test the conservation of sex chromosome in Eremias, the eight markers were further cross-tested by PCR amplification in 10 individuals of the Mongolian racerunner (Eremias argus), two of which exhibited cross-utility. The novel sex-linked markers could be mapped on the W chromosome of the sand lizard (Lacerta agilis). Our finding that the sex-linked markers are shared in closely related species, along with a conserved synteny of the W chromosome, further supports the homology and conservation of sex chromosomes in the lacertid lizards.

Genomic differentiation with isolation by distance along a latitudinal gradient in the spotted-leg treefrog Polypedates megacephalus.

The patterns of isolation by distance (IBD) entailing increased genetic differentiation among populations have aroused extensive concerns for evolutionary biologists. Although the IBD may act on spatial processes contributing to the genetic differentiation among populations in anuran species, the factors shaping the IBD of frogs among populations in natural systems are largely unknown. Here, we studied the genetic differentiation among six populations with 24 individuals of the spotted-leg treefrog along a latitudinal gradient (1860.31 km) based on 1020 single nucleotide polymorphisms from restriction site-associated DNA sequencing. The results showed that the genetic diversity differed significantly among populations and that the insular populations had higher genetic diversity than the mainland populations. Furthermore, we also found a significant genetic differentiation among populations (FST = 0.277) and no sign of inbreeding (FUNI = -0.145). The IBD was detected for all populations, and a higher degree of the IBD was indicated when controlling for the effects of the isolation between Hainan and mainland populations caused by the Qiongzhou Strait. Our findings suggest that the form of the Qiongzhou Strait plays a key role in shaping the genetic diversity and population differentiation in treefrogs.

UCE Phylogenomics, detection of a putative hybrid population, and one older mitogenomic node age of Batrachuperus salamanders.

The prevalence of incomplete lineage sorting complicates the examination of hybridization and species-level paraphyly with gene trees of a small number of loci. In Asian mountain salamanders of the genus Batrachuperus, possible hybridization and species paraphyly had been identified by utilizing mitochondrial genealogy and fixed allozyme differences. Here we sampled 2909 UCEs in 44 local populations from all six Batrachuperus species, inferred gene and species trees, compared them with mitochondrial and allozyme results, and examined the potential hybridization and species paraphyly. The clustering pattern of single-locus trees, increased proportion of heterozygous SNPs, allele frequency-based migration edge estimation, and intrapopulation long branches (as expected from an increase of genetic lineage and nucleotide diversity) support that an eastern B. karlschmidti population has experienced admixture with B. tibetanus. On the 2909-UCE concatenated and species trees, lower nodal supports were observed when similar proportions of loci agreed with alternative topologies, i.e., a reciprocal monophyly between a Pengxian lineage and the remainder of B. pinchonii (0.379) or a paraphyly of the latter with respect to Pengxian (0.362). The UCE phylogenomics agreed with the relatively recent groupings in the allozyme dendrogram. Despite incomplete lineage sorting, the mitochondrial trees were similar to the UCE trees for deeper relationships of the genus. However, one significant branch-length level discordance was identified. The branch between the common ancestor of B. daochengensis and B. yenyuanensis and common ancestor of the genus was approximately three times shorter on the mitochondrial tree than on the UCE tree, suggesting that the split of the mitochondrial lineages was likely a few million years earlier than the split of species. This finding supports considering possible ancestral polymorphism when interpreting different divergence dates estimated from mitochondrial and genome-wide data.

The complete mitochondrial genome of Batrachuperus sp. 2 (Caudata: Hynobiidae).

We describe the mitochondrial genome sequence of a mountain salamander, Batrachuperus sp.2 (Hynobiidae). The mitogenome was extracted and assembled from a large genome by next-generation sequencing. The mitochondrial genome sequences were 16,394 bp in size with 37 genes including 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and a non-coding control region (D-loop). The overall base composition of Batrachuperus sp.2 is 34.1% for A, 19.3% for C, 13.5% for G, 33.1% for T, and 32.8% for GC content. The phylogenetic tree was constructed based on 13 protein-coding genes of Batrachuperus sp.2 and 6 closely related species by RAxML.

Dynamics behind disjunct distribution, hotspot-edge refugia, and discordant RADseq/mtDNA variability: insights from the Emei mustache toad.

BACKGROUND: The distribution of genetic diversity and the underlying processes are important for conservation planning but are unknown for most species and have not been well studied in many regions. In East Asia, the Sichuan Basin and surrounding mountains constitute an understudied region that exhibits a "ring" of high species richness overlapping the eastern edge of the global biodiversity hotspot Mountains of Southwest China. We examine the distributional history and genetic diversification of the Emei mustache toad Leptobrachium boringii, a typical "ring" element characterized by disjunct ranges in the mountains, by integrating time-calibrated gene tree, genetic variability, individual-level clustering, inference of population splitting and mixing from allele frequencies, and paleoclimatic suitability modeling. RESULTS: The results reveal extensive range dynamics, including secondary contact after long-term isolation via westward dispersal accompanied by variability loss. They allow the proposal of a model that combines recurrent contractions caused by Quaternary climatic changes and some failed expansions under suitable conditions for explaining the shared disjunct distribution pattern. Providing exceptional low-elevation habitats in the hotspot area, the eastern edge harbors both long-term refugial and young immigrant populations. This finding and a synthesis of evidence from other taxa demonstrate that a certain contributor to biodiversity, one that preserves and receives low-elevation elements of the east in this case, can be significant for only a particular part of a hotspot. By clarifying the low variability of these refugial populations, we show that discordant mitochondrial estimates of diversity can be obtained for populations that experienced admixture, which would have unlikely left proportional immigrant alleles for each locus. CONCLUSIONS: Dispersal after long-term isolation can explain much of the spatial distribution of genetic diversity in this species, while secondary contact and long-term persistence do not guarantee a large variation. The model for the formation of disjunct ranges may apply to many other taxa isolated in the mountains surrounding the Sichuan Basin. Furthermore, this study provides insights into the heterogeneous nature of hotspots and discordant variability obtained from genome-wide and mitochondrial data.

Microsatellite development from genome skimming and transcriptome sequencing: comparison of strategies and lessons from frog species.

BACKGROUND: Even though microsatellite loci frequently have been isolated using recently developed next-generation sequencing (NGS) techniques, this task is still difficult because of the subsequent polymorphism screening requires a substantial amount of time. Selecting appropriate polymorphic microsatellites is a critical issue for ecological and evolutionary studies. However, the extent to which assembly strategy, read length, sequencing depth, and library layout produce a measurable effect on microsatellite marker development remains unclear. Here, we use six frog species for genome skimming and two frog species for transcriptome sequencing to develop microsatellite markers, and investigate the effect of different isolation strategies on the yield of microsatellites. RESULTS: The results revealed that the number of isolated microsatellites increases with increased data quantity and read length. Assembly strategy could influence the yield and the polymorphism of microsatellite development. Larger k-mer sizes produced fewer total number of microsatellite loci, but these loci had a longer repeat length, suggesting greater polymorphism. However, the proportion of each type of nucleotide repeats was not affected; dinucleotide repeats were always the dominant type. Finally, the transcriptomic microsatellites displayed lower levels of polymorphisms and were less abundant than genomic microsatellites, but more likely to be functionally linked loci. CONCLUSIONS: These observations provide deep insight into the evolution and distribution of microsatellites and how different isolation strategies affect microsatellite development using NGS.

Examining the interglacial high-elevation refugia scenario in East Asian subtropical mountain systems with the frog species Leptobrachium liui.

The effects of Quaternary climatic oscillations on the distributions of organisms in different parts of the world are not equally well understood, limiting the ability to understand the determinants of biodiversity. Compared with the mountain regions in southern Europe and southwestern North America, such effects on high-elevation species in the East Asian subtropical mountain systems located in southern and southeastern China have seldom been addressed. In this study, using Leptobrachium liui (Megophryidae), we made one of the earliest attempts to examine the interglacial high-elevation refugia scenario in these Asian mountains. Based on our current understanding of the study system, we formulated a hypothesis that these frogs of western origin were distributed more widely and continuously during glacial phases, allowing eastward dispersal, and that they are currently isolated in interglacial refugia at higher elevations. Microsatellite data and mitochondrial and nuclear sequence data were obtained with extensive sampling followed by the synthesis of phylogeographic and population genetic analyses and modeling of the species distribution. The analyses revealed a sequential eastward divergence of microsatellite clusters and gene lineages accompanied by a decline in genetic diversity. Molecular dating estimates revealed divergence events in the Pleistocene, and a reduction in local populations was inferred to have occurred at a time comparable to the end of the last glacial. Strong genetic isolation by distance reflecting a more continuous historical distribution was detected. Furthermore, environmental niche models inferred a wide planar distribution during the last glacial maximum, providing further support for the hypothesis.

Molecular Phylogenies indicate a Paleo-Tibetan Origin of Himalayan Lazy Toads (Scutiger).

The Himalaya presents an outstanding geologically active orogen and biodiversity hotspot. However, our understanding of the historical biogeography of its fauna is far from comprehensive. Many taxa are commonly assumed to have originated from China-Indochina and dispersed westward along the Himalayan chain. Alternatively, the "Tibetan-origin hypothesis" suggests primary diversification of lineages in Paleo-Tibet, and secondary diversification along the slopes of the later uplifted Greater Himalaya. We test these hypotheses in high-mountain megophryid anurans (Scutiger). Extensive sampling from High Asia, and analyses of mitochondrial (2839 bp) and nuclear DNA (2208 bp), using Bayesian and Maximum likelihood phylogenetics, suggest that the Himalayan species form a distinct clade, possibly older than those from the eastern Himalaya-Tibet orogen. While immigration from China-Indochina cannot be excluded, our data may indicate that Himalayan Scutiger originated to the north of the Himalaya by colonization from Paleo-Tibet and then date back to the Oligocene. High intraspecific diversity of Scutiger implies limited migration across mountains and drainages along the Himalaya. While our study strengthens support for a "Tibetan-origin hypothesis", current sampling (10/22 species; 1 revalidated: S. occidentalis) remains insufficient to draw final conclusions on Scutiger but urges comparative phylogeographers to test alternative, geologically supported hypotheses for a true future understanding of Himalayan biogeography.

Next-generation sequencing of mixed genomic DNA allows efficient assembly of rearranged mitochondrial genomes in Amolops chunganensis and Quasipaa boulengeri.

Recent improvements in next-generation sequencing (NGS) technologies can facilitate the obtainment of mitochondrial genomes. However, it is not clear whether NGS could be effectively used to reconstruct the mitogenome with high gene rearrangement. These high rearrangements would cause amplification failure, and/or assembly and alignment errors. Here, we choose two frogs with rearranged gene order, Amolops chunganensis and Quasipaa boulengeri, to test whether gene rearrangements affect the mitogenome assembly and alignment by using NGS. The mitogenomes with gene rearrangements are sequenced through Illumina MiSeq genomic sequencing and assembled effectively by Trinity v2.1.0 and SOAPdenovo2. Gene order and contents in the mitogenome of A. chunganensis and Q. boulengeri are typical neobatrachian pattern except for rearrangements at the position of "WANCY" tRNA genes cluster. Further, the mitogenome of Q. boulengeri is characterized with a tandem duplication of trnM. Moreover, we utilize 13 protein-coding genes of A. chunganensis, Q. boulengeri and other neobatrachians to reconstruct the phylogenetic tree for evaluating mitochondrial sequence authenticity of A. chunganensis and Q. boulengeri. In this work, we provide nearly complete mitochondrial genomes of A. chunganensis and Q. boulengeri.

Intraspecific rearrangement of mitochondrial genome suggests the prevalence of the tandem duplication-random loss (TDLR) mechanism in Quasipaa boulengeri.

BACKGROUND: Tandem duplication followed by random loss (TDRL) is the most frequently invoked model to explain the diversity of gene rearrangements in metazoan mitogenomes. The initial stages of gene rearrangement are difficult to observe in nature, which limits our understanding of incipient duplication events and the subsequent process of random loss. Intraspecific gene reorganizations may represent intermediate states, and if so they potentially shed light on the evolutionary dynamics of TDRL. RESULTS: Nucleotide sequences in a hotspot of gene-rearrangement in 28 populations of a single species of frog, Quasipaa boulengeri, provide such predicted intermediate states. Gene order and phylogenetic analyses support a single tandem duplication event and a step-by-step process of random loss. Intraspecific gene rearrangements are not commonly found through comparison of all mitochondrial DNA records of amphibians and squamate reptiles in GenBank. CONCLUSIONS: The intraspecific variation in Q. boulengeri provides insights into the rate of partial duplications and deletions within a mitogenome, and reveals that fixation and gene-distribution in mitogenomic reorganization is likely non-adaptive.

Combining phylogenomic and supermatrix approaches, and a time-calibrated phylogeny for squamate reptiles (lizards and snakes) based on 52 genes and 4162 species.

Two common approaches for estimating phylogenies in species-rich groups are to: (i) sample many loci for few species (e.g. phylogenomic approach), or (ii) sample many species for fewer loci (e.g. supermatrix approach). In theory, these approaches can be combined to simultaneously resolve both higher-level relationships (with many genes) and species-level relationships (with many taxa). However, fundamental questions remain unanswered about this combined approach. First, will higher-level relationships more closely resemble those estimated from many genes or those from many taxa? Second, will branch support increase for higher-level relationships (relative to the estimate from many taxa)? Here, we address these questions in squamate reptiles. We combined two recently published datasets, one based on 44 genes for 161 species, and one based on 12 genes for 4161 species. The likelihood-based tree from the combined matrix (52 genes, 4162 species) shared more higher-level clades with the 44-gene tree (90% vs. 77% shared). Branch support for higher level-relationships was marginally higher than in the 12-gene tree, but lower than in the 44-gene tree. Relationships were apparently not obscured by the abundant missing data (92% overall). We provide a time-calibrated phylogeny based on extensive sampling of genes and taxa as a resource for comparative studies.

Do missing data influence the accuracy of divergence-time estimation with BEAST?

Time-calibrated phylogenies have become essential to evolutionary biology. A recurrent and unresolved question for dating analyses is whether genes with missing data cells should be included or excluded. This issue is particularly unclear for the most widely used dating method, the uncorrelated lognormal approach implemented in BEAST. Here, we test the robustness of this method to missing data. We compare divergence-time estimates from a nearly complete dataset (20 nuclear genes for 32 species of squamate reptiles) to those from subsampled matrices, including those with 5 or 2 complete loci only and those with 5 or 8 incomplete loci added. In general, missing data had little impact on estimated dates (mean error of ∼5Myr per node or less, given an overall age of ∼220Myr in squamates), even when 80% of sampled genes had 75% missing data. Mean errors were somewhat higher when all genes were 75% incomplete (∼17Myr). However, errors increased dramatically when only 2 of 9 fossil calibration points were included (∼40Myr), regardless of missing data. Overall, missing data (and even numbers of genes sampled) may have only minor impacts on the accuracy of divergence dating with BEAST, relative to the dramatic effects of fossil calibrations.

The evolution of mitochondrial genomes in modern frogs (Neobatrachia): nonadaptive evolution of mitochondrial genome reorganization.

BACKGROUND: Although mitochondrial (mt) gene order is highly conserved among vertebrates, widespread gene rearrangements occur in anurans, especially in neobatrachians. Protein coding genes in the mitogenome experience adaptive or purifying selection, yet the role that selection plays on genomic reorganization remains unclear. We sequence the mitogenomes of three species of Glandirana and hot spots of gene rearrangements of 20 frog species to investigate the diversity of mitogenomic reorganization in the Neobatrachia. By combing these data with other mitogenomes in GenBank, we evaluate if selective pressures or functional constraints act on mitogenomic reorganization in the Neobatrachia. We also look for correlations between tRNA positions and codon usage. RESULTS: Gene organization in Glandirana was typical of neobatrachian mitogenomes except for the presence of pseudogene trnS (AGY). Surveyed ranids largely exhibited gene arrangements typical of neobatrachian mtDNA although some gene rearrangements occurred. The correlation between codon usage and tRNA positions in neobatrachians was weak, and did not increase after identifying recurrent rearrangements as revealed by basal neobatrachians. Codon usage and tRNA positions were not significantly correlated when considering tRNA gene duplications or losses. Change in number of tRNA gene copies, which was driven by genomic reorganization, did not influence codon usage bias. Nucleotide substitution rates and dN/dS ratios were higher in neobatrachian mitogenomes than in archaeobatrachians, but the rates of mitogenomic reorganization and mt nucleotide diversity were not significantly correlated. CONCLUSIONS: No evidence suggests that adaptive selection drove the reorganization of neobatrachian mitogenomes. In contrast, protein-coding genes that function in metabolism showed evidence for purifying selection, and some functional constraints appear to act on the organization of rRNA and tRNA genes. As important nonadaptive forces, genetic drift and mutation pressure may drive the fixation and evolution of mitogenomic reorganizations.

The complete mitochondrial genome of Quasipaa boulengeri (Anura: Dicroglossidae).

Quasipaa boulengeri, an Asian spiny frog that belongs to Neobatrachia, has been reported as a de novo case of chromosomal rearrangements. In this paper, we determined the whole mitochondrial DNA (mtDNA) sequence of Q. boulengeri. The complete mtDNA sequence of Q. boulengeri is 17,741 bp in length, which contains 13 protein-coding genes (PCGs), two rRNAs, 23 tRNAs and one control region. The overall nucleotides base composition of the complete mtDNA is A (28.86%), G (14.43%), C (26.09%), T (30.62%), with a much higher A + T content. Compared to other Neobatrachia species, the mtDNA genome of Q. boulengeri shows a novel gene order at WANCY region.

A de novo case of floating chromosomal polymorphisms by translocation in Quasipaa boulengeri (Anura, Dicroglossidae).

Very few natural polymorphisms involving interchromosomal reciprocal translocations are known in amphibians even in vertebrates. In this study, thirty three populations, including 471 individuals of the spiny frog Quasipaa boulengeri, were karyotypically examined using Giemsa stain or FISH. Five different karyomorphs were observed. The observed heteromorphism was autosomal but not sex-related, as the same heteromorphic chromosomes were found both in males and females. Our results indicated that the variant karyotypes resulted from a mutual interchange occurring between chromosomes 1 and 6. The occurrence of a nearly whole-arm translocation between chromosome no. 1 and no. 6 gave rise to a high frequency of alternate segregation and probably resulted in the maintenance of the translocation polymorphisms in a few populations. The translocation polymorphism is explained by different frequencies of segregation modes of the translocation heterozygote during meiosis. Theoretically, nine karyomorphs should be investigated, however, four expected karyotypes were not found. The absent karyomorphs may result from recessive lethal mutations, position effects, duplications and deficiencies. The phylogenetic inference proved that all populations of Q. boulengeri grouped into a monophyletic clade. The mutual translocation likely evolved just once in this species and the dispersal of the one karyomorph (type IV) can explain the chromosomal variations among populations.

Coalescence patterns of endemic Tibetan species of stream salamanders (Hynobiidae: Batrachuperus).

Orogenesis of topographically diverse montane regions often drives complex evolutionary histories of species. The extensive biodiversity of the eastern edge of the Tibetan Plateau, which gradually decreases eastwardly, facilitates a comparison of historical patterns. We use coalescence methods to compare species of stream salamanders (Batrachuperus) that occur at high and low elevations. Coalescent simulations reveal that closely related species are likely to have been influenced by different drivers of diversification. Species living in the western high-elevation region with its northsouth extending mountains appear to have experienced colonization via dispersal followed by isolation and divergence. In contrast, species on the eastern low-elevation region, which has many discontinuous mountain ranges, appear to have experienced fragmentation, sometimes staged, of wide-ranging ancestral populations. The two groups of species appear to have been affected differently by glaciation. High-elevation species, which are more resistant to cooler temperatures, appear to have experienced population declines as recently as the last glaciation (0.016-0.032Ma). In contrast, salamanders dwelling in the warmer and wetter habitats at low-elevation environs appear to have been affected less by the relatively recent, milder glaciation, and more so by harsher, extensive glaciations (0.5-0.175 Ma). Thus, elevation, topography and cold tolerance appear to drive evolutionary patterns of diversification and demography even among closely related taxa. The comparison of multiple species in genealogical analyses can lead to an understanding of the evolutionary drivers.

COI is better than 16S rRNA for DNA barcoding Asiatic salamanders (Amphibia: Caudata: Hynobiidae).

The 5' region of the mitochondrial DNA (mtDNA) gene cytochrome c oxidase I (COI) is the standard marker for DNA barcoding. However, because COI tends to be highly variable in amphibians, sequencing is often challenging. Consequently, another mtDNA gene, 16S rRNA gene, is often advocated for amphibian barcoding. Herein, we directly compare the usefulness of COI and 16S in discriminating species of hynobiid salamanders using 130 individuals. Species identification and classification of these animals, which are endemic to Asia, are often based on morphology only. Analysis of Kimura 2-parameter genetic distances (K2P) documents the mean intraspecific variation for COI and 16S rRNA genes to be 1.4% and 0.3%, respectively. Whereas COI can always identify species, sometimes 16S cannot. Intra- and interspecific genetic divergences occasionally overlap in both markers, thus reducing the value of a barcoding gap to identify genera. Regardless, COI is the better DNA barcoding marker for hynobiids. In addition to the comparison of two potential markers, high levels of intraspecific divergence in COI (>5%) suggest that both Onychodactylus fischeri and Salamandrella keyserlingii might be composites of cryptic species.