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.

Impacts of climate change on herpetofauna diversity in the Qinghai-Tibetan Plateau.

Although numerous studies on the impacts of climate change on biodiversity have been published, only a handful are focused on the intraspecific level or consider population-level models (separate models per population). We endeavored to fill this knowledge gap relative to the Qinghai-Tibetan plateau (QTP) by combining species distribution modeling (SDMs) with population genetics (i.e., population-level models) and phylogenetic methods (i.e., phylogenetic tree reconstruction and phylogenetic diversity analyses). We applied our models to 11 endemic and widely distributed herpetofauna species inhabiting high elevations in the QTP. We aimed to determine the influence of environmental heterogeneity on species' responses to climate change, the magnitude of climate-change impacts on intraspecific diversity, and the relationship between species range loss and intraspecific diversity losses under 2 shared socioeconomic pathways (SSP245 and SSP585) and 3 future periods (2050s, 2070s, and 2090s). The effects of global climatic change were more pronounced at the intraspecific level (22% of haplotypes lost and 36% of populations lost) than the morphospecies level in the SSP585 climate change scenario. Maintenance of genetic diversity was in general determined by a combination of factors including range changes, species genetic structure, and the part of the range predicted to be lost. This is owing to the fact that the loss and survival of populations were observed in species irrespective of the predicted range changes (contraction or expansion). In the southeast (mountainous regions), climate change had less of an effect on range size (>100% in 3 species) than in central and northern QTP plateau regions (range size <100% in all species). This may be attributed to environmental heterogeneity, which provided pockets of suitable climate in the southeast, whereas ecosystems in the north and central regions were homogeneous. Generally, our results imply that mountainous regions with high environmental heterogeneity and high genetic diversity may buffer the adverse impacts of climate change on species distribution and intraspecific diversity. Therefore, genetic structure and characteristics of the ecosystem may be crucial for conservation under climate change.

Impactos del cambio climático sobre la diversidad de herpetofauna en la meseta Qinghai-Tíbet Región Aunque se han publicado numerosos estudios sobre los impactos del cambio climática en la biodiversidad, son muy pocos los que se enfocan en el nivel intraespecífico o que consideran modelos a nivel poblacional (modelos separados por población). Intentamos cerrar este vacío de conocimiento en relación con la meseta Qinghai-Tíbet (MQT) con la combinación entre modelos de distribución de especies (MDE) y genética poblacional (modelos a nivel poblacional) y métodos filogenéticos (reconstrucción de árboles filogenéticos y análisis de diversidad filogenética). Aplicamos nuestros modelos a once especies endémicas de herpetofauna con distribución amplia en las elevaciones más altas de la MQT. Nos planteamos determinar la influencia de la heterogeneidad de las especies sobre la respuesta de las especies al cambio climático, la magnitud de los impactos del cambio climático sobre la diversidad intraespecífica y la relación entre la pérdida de distribución de la especie y las pérdidas de diversidad intraespecífica bajo dos vías socioeconómicas (SSP245 y SSP585) y tres periodos del futuro (2050s, 2070s y 2090s). Los efectos del cambio climático global fueron más pronunciados a nivel intraespecífico (22% de pérdida en los haplotipos y 36% en las poblaciones) que al nivel morfoespecie en el escenario de cambio climático SSP585. El mantenimiento de la diversidad genética casi siempre estuvo determinado por una combinación de factores que incluyen cambios en la distribución, estructura genética de las especies y la parte de la distribución que se pronosticó se perdería. Esto se debe a que observamos la pérdida y supervivencia de las poblaciones sin importar los cambios pronosticados en la distribución (contracción o expansión). En las regiones montañosas del sureste, el cambio climático tuvo un efecto menor sobre la distribución (>100% en tres especies) comparado con las regiones de la meseta central y del norte de la MQT (distribución <100% en todas las especies). Esto puede atribuirse a la heterogeneidad ambiental, la cual proporciona recovecos de clima adecuado en el sureste, mientras que los ecosistemas en las regiones central y norte fueron homogéneos. De manera general, nuestros resultados implican que las regiones montañosas con una elevada heterogeneidad ambiental y una gran diversidad genética podrían reducir los impactos adversos del cambio climático sobre la distribución de las especies y la diversidad intraespecífica. Por lo tanto, la estructura genética y las características del ecosistema pueden ser cruciales para conservar bajo el cambio climático.

DNA barcoding of Chinese snakes reveals hidden diversity and conservation needs.

DNA barcoding has greatly facilitated studies of taxonomy, biodiversity, biological conservation, and ecology. Here, we establish a reliable DNA barcoding library for Chinese snakes, unveiling hidden diversity with implications for taxonomy, and provide a standardized tool for conservation management. Our comprehensive study includes 1638 cytochrome c oxidase subunit I (COI) sequences from Chinese snakes that correspond to 17 families, 65 genera, 228 named species (80.6% of named species) and 36 candidate species. A barcode gap analysis reveals gaps, where all nearest neighbour distances exceed maximum intraspecific distances, in 217 named species and all candidate species. Three species-delimitation methods (ABGD, sGMYC, and sPTP) recover 320 operational taxonomic units (OTUs), of which 192 OTUs correspond to named and candidate species. Twenty-eight other named species share OTUs, such as Azemiops feae and A. kharini, Gloydius halys, G. shedaoensis, and G. intermedius, and Bungarus multicinctus and B. candidus, representing inconsistencies most probably caused by imperfect taxonomy, recent and rapid speciation, weak taxonomic signal, introgressive hybridization, and/or inadequate phylogenetic signal. In contrast, 43 species and candidate species assign to two or more OTUs due to having large intraspecific distances. If most OTUs detected in this study reflect valid species, including the 36 candidate species, then 30% more species would exist than are currently recognized. Several OTU divergences associate with known biogeographic barriers, such as the Taiwan Strait. In addition to facilitating future studies, this reliable and relatively comprehensive reference database will play an important role in the future monitoring, conservation, and management of Chinese snakes.

First national record of Microhylahmongorum Hoang, Nguyen, Phan, Pham, Ninh, Wang, Jiang, Ziegler and Nguyen, 2022 (Anura, Microhylidae, Microhyla) in China.

Background: To date, 10 species of the genus Microhyla have been recorded in China, of which six were distributed in Yunnan Province. Microhylahmongorum Hoang, Nguyen, Phan, Pham, Ninh, Wang, Jiang, Ziegler, and Nguyen, 2022 was also speculated to be distributed in Xishuangbana, Yunnan Province, China. However, there is no evidence of documentation of M.hmongorum. New information: We report the first country record of Microhylahmongorum, based on specimens collected from Yunnan border region. Morphologically, the specimen was consistent with the original descriptions of M.hmongorum. Phylogenetically, the sequences of the specimens from China clustered with the sequence of type specimens of M.hmongorum from Vietnam, with uncorrected pairwise distances of 0.9% at the 16S gene fragment analysed. Therefore, we report M.hmongorum as a new record species in China.

The highest-elevation frog provides insights into mechanisms and evolution of defenses against high UV radiation.

Defense against ultraviolet (UV) radiation exposure is essential for survival, especially in high-elevation species. Although some specific genes involved in UV response have been reported, the full view of UV defense mechanisms remains largely unexplored. Herein, we used integrated approaches to analyze UV responses in the highest-elevation frog, Nanorana parkeri. We show less damage and more efficient antioxidant activity in skin of this frog than those of its lower-elevation relatives after UV exposure. We also reveal genes related to UV defense and a corresponding temporal expression pattern in N. parkeri. Genomic and metabolomic analysis along with large-scale transcriptomic profiling revealed a time-dependent coordinated defense mechanism in N. parkeri. We also identified several microRNAs that play important regulatory roles, especially in decreasing the expression levels of cell cycle genes. Moreover, multiple defense genes (i.e., TYR for melanogenesis) exhibit positive selection with function-enhancing substitutions. Thus, both expression shifts and gene mutations contribute to UV adaptation in N. parkeri. Our work demonstrates a genetic framework for evolution of UV defense in a natural environment.

Discovery of a wild, genetically pure Chinese giant salamander creates new conservation opportunities.

Effective conservation of threatened biota relies on accurate assessments and scientific guidance. As an unfortunate example, Chinese giant salamanders ( Andrias, CGS) remain critically endangered in nature. Misguided conservation efforts, e.g., commercial propagation and releasing of millions of likely non-indigenous or interspecific hybrids, have further compromised conservation initiatives. Limited information on wild populations of CGS poses a significant conservation challenge. Following 18-month long field monitoring, we now report the discovery of a wild population of CGS in a closed nature reserve in Jiangxi Province, China. Genomic assessments reveal its genetic distinctiveness and do not detect genetic admixture with other species. Based on morphological and molecular evidences, we describe this CGS as a new species Andrias jiangxiensis sp. nov. This is the only known species of CGS today with a genetically pure, reproducing, in situ population. This discovery emphasizes the important role that closed nature reserves play in protecting species, and the necessity of integrating long-term field monitoring and genetic assessments. It sets a new pathway for discovering and conserving endangered species, especially for those biotas that are similarly being extirpated by anthropogenic translocations and overexploitation.

Species Persistence with Hybridization in Toad-Headed Lizards Driven by Divergent Selection and Low Recombination.

Speciation plays a central role in evolutionary studies, and particularly how reproductive isolation (RI) evolves. The origins and persistence of RI are distinct processes that require separate evaluations. Treating them separately clarifies the drivers of speciation and then it is possible to link the processes to understand large-scale patterns of diversity. Recent genomic studies have focused predominantly on how species or RI originate. However, we know little about how species persist in face of gene flow. Here, we evaluate a contact zone of two closely related toad-headed lizards (Phrynocephalus) using a chromosome-level genome assembly and population genomics. To some extent, recent asymmetric introgression from Phrynocephalus putjatai to P. vlangalii reduces their genomic differences. However, their highly divergent regions (HDRs) have heterogeneous distributions across the genomes. Functional gene annotation indicates that many genes within HDRs are involved in reproduction and RI. Compared with allopatric populations, contact areas exhibit recent divergent selection on the HDRs and a lower population recombination rate. Taken together, this implies that divergent selection and low genetic recombination help maintain RI. This study provides insights into the genomic mechanisms that drive RI and two species persistence in the face of gene flow during the late stage of speciation.

Herpetological phylogeographic analyses support a Miocene focal point of Himalayan uplift and biological diversification.

The Himalaya are among the youngest and highest mountains in the world, but the exact timing of their uplift and origins of their biodiversity are still in debate. The Himalayan region is a relatively small area but with exceptional diversity and endemism. One common hypothesis to explain the rich montane diversity is uplift-driven diversification-that orogeny creates conditions favoring rapid in situ speciation of resident lineages. We test this hypothesis in the Himalayan region using amphibians and reptiles, two environmentally sensitive vertebrate groups. In addition, analysis of diversification of the herpetofauna provides an independent source of information to test competing geological hypotheses of Himalayan orogenesis. We conclude that the origins of the Himalayan herpetofauna date to the early Paleocene, but that diversification of most groups was concentrated in the Miocene. There was an increase in both rates and modes of diversification during the early to middle Miocene, together with regional interchange (dispersal) between the Himalaya and adjacent regions. Our analyses support a recently proposed stepwise geological model of Himalayan uplift beginning in the Paleocene, with a subsequent rapid increase of uplifting during the Miocene, finally giving rise to the intensification of the modern South Asian Monsoon.

Multi-locus genetic analyses of Quasipaa from throughout its distribution.

Montane frogs of the genus Quasipaa Dubois, 1992 occur from southern China to Southeast Asia (Frost 2021). Analyses of mtDNA (Cytb) and nuDNA data (Rag1, Rag2, Rhod, Tyr) for samples from 93 localities throughout its distribution yield a phylogeny. Clades A and B occur in Southeast Asia, clade C in northern Yangtze River, China, clade D in southwestern China, and clades E and F in southeastern China. Results place Q. yei within monophyletic Quasipaa and identify two new species. Based on nuDNA data, the basal split of clade A and B indicates an Indochinese origin of Quasipaa. The west-east diversification of five species across South China (Q. spinosa, Q. exilispinosa, Q. jiulongensis, Q. shini, Q. boulengeri) corresponds to topographic terrains II and III of China. Divergence of species from southeastern China (Q. shini, Q. jiulongensis, Q. spinosa, Q. exilispinosa) and southwestern China (Q. boulengeri) dates to 15.30-16.56 Ma (million years ago). A principal component analysis (PCA) and t-test involving 19 bioclimatic variables identifies significantly different environmental conditions between the two regions. Species' distribution models (SDM) for Q. spinosa and Q. boulengeri identify the best areas to be eastern and western South China, respectively. Thus, environmental variation appears to have influenced the genetic divergence and distributions of Quasipaa in South China. Mito-nuclear discordance indicates that some individuals of Q. exilispinosa and Q. spinosa hybridized historically.

Author Correction: 863 genomes reveal the origin and domestication of chicken.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

863 genomes reveal the origin and domestication of chicken.

Despite the substantial role that chickens have played in human societies across the world, both the geographic and temporal origins of their domestication remain controversial. To address this issue, we analyzed 863 genomes from a worldwide sampling of chickens and representatives of all four species of wild jungle fowl and each of the five subspecies of red jungle fowl (RJF). Our study suggests that domestic chickens were initially derived from the RJF subspecies Gallus gallus spadiceus whose present-day distribution is predominantly in southwestern China, northern Thailand and Myanmar. Following their domestication, chickens were translocated across Southeast and South Asia where they interbred locally with both RJF subspecies and other jungle fowl species. In addition, our results show that the White Leghorn chicken breed possesses a mosaic of divergent ancestries inherited from other subspecies of RJF. Despite the strong episodic gene flow from geographically divergent lineages of jungle fowls, our analyses show that domestic chickens undergo genetic adaptations that underlie their unique behavioral, morphological and reproductive traits. Our study provides novel insights into the evolutionary history of domestic chickens and a valuable resource to facilitate ongoing genetic and functional investigations of the world's most numerous domestic animal.

A revised taxonomy of Asian snail-eating snakes Pareas (Squamata, Pareidae): evidence from morphological comparison and molecular phylogeny.

The Asian snail-eating snakes Pareas is the largest genus of the family Pareidae (formerly Pareatidae), and widely distributed in Southeast Asia. However, potential diversity remains poorly explored due to their highly conserved morphology and incomplete samples. Here, on basis of more extensive sampling, interspecific phylogenetic relationships of the genus Pareas were reconstructed using two mitochondrial fragments (cyt b and ND4) and two nuclear genes (c-mos and Rag1), and multivariate morphometrics conducted for external morphological data. Both Bayesian Inference and Maximum Likelihood analyses consistently showed that the genus Pareas was comprised of two distinct, monophyletic lineages with moderate to low support values. Based on evidences from molecular phylogeny and morphological data, cryptic diversity of this genus was uncovered and two new species were described. In additional, the validity of P. macularius is confirmed.

An integrative phylogenomic approach illuminates the evolutionary history of Old World tree frogs (Anura: Rhacophoridae).

Rhacophoridae are one of the most speciose and ecologically diverse families of amphibians. Resolution of their evolutionary relationships is key to understanding the accumulation of biodiversity, yet previous hypotheses based on Sanger sequencing exhibit much discordance amongst generic relationships. This conflict precludes the making of sound macroevolutionary conclusions. Herein, we conduct the first phylogenomic study using broad-scale sampling and sequences of 352 nuclear DNA loci obtained using anchored hybrid enrichment targeted sequencing. The robust time-calibrated phylogenetic hypothesis clarifies several long-disputed relationships and facilitates the testing of evolutionary hypotheses on spatiotemporal diversification and reproductive modes. The major extant lineages of Rhacophoridae appear to have radiated in mainland Asia, and the spatiotemporal process corresponds with several common accumulations of biodiversity in Asia. Analyses do not detect any case of "Out of Himalaya" in Rhacophoridae. All transitions of reproductive modes appear to have evolved in an ordered, gradual sequence associated with gaining independence of standing water for larval development. The different reproductive modes are phylogenetically conserved and the completion of their transitions appear to have occurred over a period of ~30 Ma, which does not fit a pattern of a rapid burst of diversification. Innovations in reproductive modes associate statistically with the uneven distribution of species-richness between clades, where higher diversification is linked to increased terrestrial modes of reproduction. These results strengthen the hypothesis that breeding innovations drive diversification by providing new opportunities for ecological release and dispersion.

Genetic variation and cryptic lineage diversity of the Nigerian red-headed rock agama Agama agama associate with eco-geographic zones.

Nigeria is an Afrotropical region with considerable ecological heterogeneity and levels of biotic endemism. Among its vertebrate fauna, reptiles have broad distributions, thus, they constitute a compelling system for assessing the impact of ecological variation and geographic isolation on species diversification. The red-headed rock agama, Agama agama, lives in a wide range of habitats and, thus, it may show genetic structuring and diversification. Herein, we tested the hypothesis that ecology affects its genetic structure and population divergence. Bayesian inference phylogenetic analysis of a mitochondrial DNA (mtDNA) gene recovered four well-supported matrilines with strong evidence of genetic structuring consistent with eco-geographic regions. Genetic differences among populations based on the mtDNA also correlated with geographic distance. The ecological niche model for the matrilines had a good fit and robust performance. Population divergence along the environmental axes was associated with climatic conditions, and temperature ranked highest among all environmental variables for forest specialists, while precipitation ranked highest for the forest/derived savanna, and savanna specialists. Our results cannot reject the hypothesis that niche conservatism promotes geographic isolation of the western populations of Nigerian A. agama. Thus, ecological gradients and geographic isolation impact the genetic structure and population divergence of the lizards. This species might be facing threats due to recent habitat fragmentation, especially in western Nigeria. Conservation actions appear necessary.

A new species of Acanthosaura Gray 1831 (Reptilia: Agamidae) from central Vietnam.

Herein, we describe the new agamid species Acanthosaura phongdienensis sp. nov. from central Vietnam based on morphological and genetic data. Males of the new medium-sized species have a snout-vent length of up to 77.4 mm and females up to 64.7 mm. In both genders, the tail is longer than the snout-vent length. A spine occurs on each side of the neck, the rostral scute is entire, and the nuchal and dorsal crests are continuous. A postorbital spine is about half as long as the diameter of the eye-ball, and small lateral scales intermix with large, keeled scales whose tips point backwards and upwards. The new species differs genetically from the closely related species A. brachypoda, A. coronata, and A. lepidogaster by uncorrected p-distances of 13.6%, 21.7%, and 12.8%, respectively. This is the 13th species of Acanthosaura and the eighth species in the genus from Vietnam.

A new species of the genus Raorchestes (Anura: Rhacophoridae) from Yunnan Province, China.

A new bush frog species, Raorchestes cangyuanensis sp. nov., from Cangyuan, Yunnan Province, China, is described based on morphological and molecular analyses. It differs from all known congeners by a combination of the following characters: body size small, adult snout-vent length (SVL) 16.1-20.0 mm in males (n=3); tympanum indistinct; tips of all fingers and toes expanded into discs with circummarginal grooves; rudimentary webbing between toes; fingers and toes with lateral dermal fringes; inner and outer metacarpal tubercles present; heels meeting when limbs held at right angles to body; crotch with a distinct black patch; discs of fingers and toes orange; male with external single subgular vocal sac and reddish nuptial pad at the base of first finger.

A new species of the genus Xenophrys Anura Megophryidae from northern Thailand.

Species of Xenophrys are conserved morphologically and live primarily in forests. In Thailand, the genus harbors many cryptic species. Herein we report the collection of specimens from Doi Inthanon, Chiang Mai Province, northern Thailand, which were identified previously as X. minor. Molecular and morphological analyses find that these specimens differ significantly from other known congeners, and therefore we describe a new species. Further, our phylogenetic analyses indicate that X. latidactyla is a junior synonym of X. palpebralespinosa.

Genomic and transcriptomic investigations of the evolutionary transition from oviparity to viviparity.

Viviparous (live-bearing) vertebrates have evolved repeatedly within otherwise oviparous (egg-laying) clades. Over two-thirds of these changes in vertebrate reproductive parity mode happened in squamate reptiles, where the transition has happened between 98 and 129 times. The transition from oviparity to viviparity requires numerous physiological, morphological, and immunological changes to the female reproductive tract, including eggshell reduction, delayed oviposition, placental development for supply of water and nutrition to the embryo by the mother, enhanced gas exchange, and suppression of maternal immune rejection of the embryo. We performed genomic and transcriptomic analyses of a closely related oviparous-viviparous pair of lizards (Phrynocephalus przewalskii and Phrynocephalus vlangalii) to examine these transitions. Expression patterns of maternal oviduct through reproductive development of the egg and embryo differ markedly between the two species. We found changes in expression patterns of appropriate genes that account for each of the major aspects of the oviparity to viviparity transition. In addition, we compared the gene sequences in transcriptomes of four oviparous-viviparous pairs of lizards in different genera (Phrynocephalus, Eremias, Scincella, and Sphenomorphus) to look for possible gene convergence at the sequence level. We discovered low levels of convergence in both amino acid replacement and evolutionary rate shift. This suggests that most of the changes that produce the oviparity-viviparity transition are changes in gene expression, so occasional reversals to oviparity from viviparity may not be as difficult to achieve as has been previously suggested.

Revalidation and resurrection of Panaspis massaiensis (Angel, 1924) and the description of a new species of Panaspis Cope (Squamata: Scincidae) from south-eastern Kenya.

Historically Panaspis wahlbergi (Smith, 1849) has been the only assignable species present in Kenya. Recent studies have shown that it comprises multiple cryptic species and the nominal species is now restricted to southern Africa. Newly collected mitochondrial data (16S rRNA) helped to resolve the status of the Kenyan populations, which revealed the presence of two distant related species. Pairwise distances show average 5.87% differences between the two Kenyan species, and 3.58-5.27% and 8.62-9.15% to nominal P. wahlbergi and P. maculicollis Jacobsen Broadley, 2000 respectively. Ablepharus massaiensis Angel, 1924 was described from the Maasai plains near Nairobi, but has long been considered a junior synonym of P. wahlbergi. We herein resurrect Panaspis massaiensis comb. nov. as a valid species and describe a new species, Panaspis tsavoensis sp. nov. from the Tsavo Conservation Area in south-eastern Kenya. Morphological examinations of specimens reveal minor differences from each other as well as nominal forms of P. wahlbergi and P. maculicollis. Panaspis massaiensis comb. nov. shares with the P. wahlbergi group a white ventrolateral stripe but can be distinguished by the presence of 26 midbody scale rows versus 24. Panaspis tsavoensis sp. nov. on the other hand, lacks the white ventrolateral stripe, most similar to the P. maculicollis group but differs in that P. maculicollis breeding males have a black patch on the neck with diagonal rows of white spots. Panaspis massaiensis comb. nov. is widespread in the Kenyan and northern Tanzanian highlands, isolated dryland montane forests and rocky hills, while P. tsavoensis sp. nov. occur in the expansive arid lowlands of Tsavo Conservation Area and should be present in similar arid lowlands in northern Kenya as well as in adjacent Tanzania.