Characterization of the complete mitochondrial genome of Hynobius bambusicolus Wang, Othman, Qiu and Borzée, 2023 (Amphibia, Caudata, Hynobiidae) and its phylogenetic implications

Hynobius bambusicolus (Caudata, Hynobiidae) is a recently described species, identified in 2022, and is thus not widely known. In this study, we sequenced and annotated the complete mitogenome of H. bambusicolus. The resulting mitochondrial genome is 16,406 bp in length and comprises 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNA), 22 transfer RNA genes, and a non-coding region. The base composition of the mitogenome is 33.3% A, 32.1% T, 20.8% C, and 13.8% G. The phylogenetic trees indicated that H. bambusicolus is the basal branch within the Southern Chinese Hynobius clade.

A time-calibrated salamander phylogeny including 765 species and 503 genes

Recent time-calibrated amphibian phylogenies agree on the family-level relationships among extant salamanders but had disparate sampling regimes and inferred very different divergence times. For example, a recent phylogenomic study based on 220 nuclear loci had limited taxon sampling (41 species) and estimated relatively young divergence dates, whereas a more extensive supermatrix study based on 15 genes and 481 species estimated dates that were 22-45 million years older for major clades. Here, we combined phylogenomic and supermatrix approaches to estimate the largest salamander phylogeny to date based on molecular markers. Our matrix contained 765 salamander species and 503 genes (with 92.3% missing data overall). We included 284 more species than the previous largest salamander phylogeny (59% increase) and sampled approximately 93% of all currently described salamander species. Our dating analyses incorporated more than twice as many fossil calibration points within salamanders as previous studies. Maximum-likelihood estimates of tree topology yielded family-level relationships that were consistent with earlier studies. Nearly all species were placed in the expected genera, despite extensive missing data in many species. Bootstrap support was generally high across the tree but was poor in some clades where sampling of genes was limited (e.g., among some bolitoglossine salamanders). The dating analyses yielded age estimates for major clades that were generally intermediate between those from the previous phylogenomic and supermatrix analyses. We also provide a set of 200 time-calibrated trees for use in comparative analyses.

The last of their kind: Is the genus Scutiger (Anura: Megophryidae) a relict element of the paleo-Transhimalaya biota?

The orographic evolution of the Himalaya-Tibet Mountain system continues to be a subject of controversy, leading to considerable uncertainty regarding the environment and surface elevation of the Tibetan Plateau during the Cenozoic era. As many geoscientific (but not paleontological) studies suggest, elevations close to modern heights exist in vast areas of Tibet since at least the late Paleogene, implicating the presence of large-scale alpine environments for more than 30 million years. To explore a recently proposed alternative model that assumes a warm temperate environment across paleo-Tibet, we carried out a phylogeographic survey using genomic analyses of samples covering the range of endemic lazy toads (Scutiger) across the Himalaya-Tibet orogen. We identified two main clades, with several, geographically distinct subclades. The long temporal gap between the stem and crown age of Scutiger may suggest high extinction rates. Diversification within the crown group, depending on the calibration, occurred either from the Mid-Miocene or Late-Miocene and continued until the Holocene. The present-day Himalayan Scutiger fauna could have evolved from lineages that existed on the southern edges of the paleo-Tibetan area (the Transhimalaya = Gangdese Shan), while extant species living on the eastern edge of the Plateau originated probably from the eastern edges of northern parts of the ancestral Tibetan area (Hoh Xil, Tanggula Shan). Based on the Mid-Miocene divergence time estimation and ancestral area reconstruction, we propose that uplift-associated aridification of a warm temperate Miocene-Tibet, coupled with high extirpation rates of ancestral populations, and species range shifts along drainage systems and epigenetic transverse valleys of the rising mountains, is a plausible scenario explaining the phylogenetic structure of Scutiger. This hypothesis aligns with the fossil record but conflicts with geoscientific concepts of high elevated Tibetan Plateau since the late Paleogene. Considering a Late-Miocene/Pliocene divergence time, an alternative scenario of dispersal from SE Asia into the East, Central, and West Himalaya cannot be excluded, although essential evolutionary and biogeographic aspects remain unresolved within this model.

Divergence and serial colonization shape genetic variation and define conservation units in Asian elephants

Asian elephants (Elephas maximus) are the largest extant terrestrial megaherbivores native to Asia, with 60% of their wild population found in India. Despite ecological and cultural importance, their population genetic structure and diversity, demographic history, and ensuing implications for management/conservation remain understudied. We analyzed 34 whole genomes (between 11× and 32×) from most known elephant landscapes in India and identified five management/conservation units corresponding to elephants in Northern (Northwestern/Northeastern), Central, and three in Southern India. Our data reveal signatures of divergence and serial colonization and a potential dilution of genetic diversity from north to south of India. The northern populations diverged from others more than 70,000 years ago, have higher genetic diversity, and have low inbreeding (pi = 0.0016 ± 0.0001; FROH > 1 MB = 0.09 ± 0.03). Two of three populations in Southern India have low diversity and are inbred, with very low effective population sizes compared with census sizes (pi = 0.0014 ± 0.00009 and 0.0015 ± 0.0001; FROH > 1 MB = 0.25 ± 0.09 and 0.17 ± 0.02). Analyses of genetic load reveal the purging of potentially high-effect insertion/deletion (indel) deleterious alleles in the southern populations and a decreasing number of deleterious alleles from north to south in India. However, despite dilution and purging for the damaging mutation load in Southern India, the load that remains is homozygous. High homozygosity of deleterious alleles, coupled with low neutral genetic diversity, make southernmost populations high priority for conservation attention. Most surprisingly, our study suggests that patterns of genetic diversity and genetic load can correspond to genomic signatures of serial founding events, even in large, highly mobile, endangered mammals.

Museum specimens shedding light on the evolutionary history and cryptic diversity of the hedgehog family Erinaceidae

The family Erinaceidae encompasses 27 extant species in two subfamilies: Erinaceinae, which includes spiny hedgehogs, and Galericinae, which comprises silky-furred gymnures and moonrats. Although they are commonly recognized by the general public, their phylogenetic history remains incompletely understood, and several species have never been included in any molecular analyses. Additionally, previous research suggested that the species diversity of Erinaceidae might be underestimated. In this study, we sequenced the mitochondrial genomes of 29 individuals representing 18 erinaceid species using 18 freshly collected tissue and 11 historical museum specimens. We also integrated previously published data for a concatenated analysis. We aimed to elucidate the evolutionary relationships within Erinaceidae, estimate divergence times, and uncover potential underestimated species diversity. Our data finely resolved intergeneric and interspecific relationships and presented the first molecular evidence for the phylogenetic position of Mesechinus wangi, Paraechinus micropus, and P. nudiventris. Our results revealed a sister relationship between Neotetracus and Neohylomys gymnures, as well as a sister relationship between Hemiechinus and Mesechinus, supporting previous hypotheses. Additionally, our findings provided a novel phylogenetic position for Paraechinus aethiopicus, placing it in a basal position within the genus. Furthermore, our study uncovered cryptic species diversity within Hylomys suillus as well as in Neotetracus sinensis, Atelerix albiventris, P. aethiopicus, and Hemiechinus auratus, most of which have been previously overlooked.

Gene Flow and Isolation in the Arid Nearctic Revealed by Genomic Analyses of Desert Spiny Lizards

The opposing forces of gene flow and isolation are two major processes shaping genetic diversity. Understanding how these vary across space and time is necessary to identify the environmental features that promote diversification. The detection of considerable geographic structure in taxa from the arid Nearctic has prompted research into the drivers of isolation in the region. Several geographic features have been proposed as barriers to gene flow, including the Colorado River, Western Continental Divide (WCD), and a hypothetical Mid-Peninsular Seaway in Baja California. However, recent studies suggest that the role of barriers in genetic differentiation may have been overestimated when compared to other mechanisms of divergence. In this study, we infer historical and spatial patterns of connectivity and isolation in Desert Spiny Lizards (Sceloporus magister) and Baja Spiny Lizards (Sceloporus zosteromus), which together form a species complex composed of parapatric lineages with wide distributions in arid western North America. Our analyses incorporate mitochondrial sequences, genomic-scale data, and past and present climatic data to evaluate the nature and strength of barriers to gene flow in the region. Our approach relies on estimates of migration under the multispecies coalescent to understand the history of lineage divergence in the face of gene flow. Results show that the S. magister complex is geographically structured, but we also detect instances of gene flow. The WCD is a strong barrier to gene flow, while the Colorado River is more permeable. Analyses yield conflicting results for the catalyst of differentiation of peninsular lineages in S. zosteromus. Our study shows how large-scale genomic data for thoroughly sampled species can shed new light on biogeography. Furthermore, our approach highlights the need for the combined analysis of multiple sources of evidence to adequately characterize the drivers of divergence.

Multiple transitions between realms shape relict lineages of Proteus cave salamanders

In comparison to biodiversity on Earth's surface, subterranean biodiversity has largely remained concealed. The olm (Proteus anguinus) is one of the most enigmatic extant cave inhabitants, and until now little was known regarding its genetic structure and evolutionary history. Olms inhabit subterranean waters throughout the Dinaric Karst of the western Balkans, with a seemingly uniform phenotypic appearance of cave-specialized traits: an elongate body, snout and limbs, degenerated eyes and loss of pigmentation ("white olm"). Only a single small region in southeastern Slovenia harbours olms with a phenotype typical of surface animals: pigmented skin, eyes, a blunt snout and short limbs ("black olm"). We used a combination of mitochondrial DNA and genome-wide single nucleotide polymorphism data to investigate the molecular diversity, evolutionary history and biogeography of olms along the Dinaric Karst. We found nine deeply divergent species-level lineages that separated between 17 and 4 million years ago, while molecular diversity within lineages was low. We detected no signal of recent admixture between lineages and only limited historical gene flow. Biogeographically, the contemporaneous distribution of lineages mostly mirrors hydrologically separated subterranean environments, while the historical separation of olm lineages follows microtectonic and climatic changes in the area. The reconstructed phylogeny suggests at least four independent transitions to the cave phenotype. Two of the species-level lineages have miniscule ranges and may represent Europe's rarest amphibians. Their rarity and the decline in other lineages call for protection of their subterranean habitats.

Phylogenetic relationships and genetic differentiation of two Salamandrella species as revealed via COI gene from Northeastern China

Due to traditional classification methods' limitations, some cryptic species remain undiscovered. To better explore the existence of the Schrenck salamander (Salamandrella tridactyla, a cryptic species of Siberian salamander S. keyserlingii) in China, we conducted a molecular phylogenetic analysis to confirm the taxonomic relationship among Salamandrella species and investigate genetic variation. We used complete sequences of the mitochondrial COI gene from 65 specimens collected across a wide range in Northeastern China. Thirty-five haplotypes were obtained from six populations. They showed medium-high haplotype diversity (Hd) and low nucleotide polymorphism (π). The phylogenetic tree and haplotype network analysis revealed that populations from Greater Khingan Ridge (Huma: HM) and Lesser Khingan Ridge (Tieli: TL) belong to S. keyserlingii, while populations from Changbai Mountain (Shangzhi-zhuziying: SZ, Shangzhi-cuijia: SC, Hailin: HL, and Baishan: BS) belong to S. tridactyla. This indicates the monophyly of Salamandrella and each of the two species. There was a substantial level of genetic differentiation between different species and within populations of the same species. This differentiation was significantly related to geographical distance. At last, the mismatch distribution and neutrality analyses indicated that the TL populations have undergone expansion of history. The study supplements the distributional range of Schrenck salamander. And it provides a theoretical basis for species conservation of Salamandrella species.

Crossing the Rift valley: using complete mitogenomes to infer the diversification and biogeographic history of ethiopian highlands Ptychadena (anura: Ptychadenidae)

The Ethiopian Highlands are considered a biodiversity hotspot, harboring a high number of endemic species. Some of the endemic species probably diversified in situ; this is, for example, the case of a monophyletic clade containing 12 known species of grass frogs of the genus Ptychadena. The different species occur at elevations ranging from 1,500 to above 3,400 m and constitute excellent models to study the process of diversification in the highlands as well as adaptations to high elevations. In this study, we sampled 294 specimens across the distribution of this clade and used complete mitogenomes and genome-wide SNP data to better understand how landscape features influenced the population structure and dispersal of these grass frogs across time and space. Using phylogenetic inference, population structure analyses, and biogeographic reconstructions, we found that the species complex probably first diversified on the south-east side of the Great Rift Valley. Later on, species dispersed to the north-west side, where more recent diversification occurred. We further demonstrate that Ptychadena species have dispersed across the Great Rift Valley at different times. Our analyses allowed for a more complete understanding of the contribution of geological events, biogeographic barriers and climatic changes as drivers of species diversification and adaptation in this important biogeographic region.

Hidden species diversity and mito-nuclear discordance within the Mediterranean cone snail, Lautoconus ventricosus

The Mediterranean cone snail, Lautoconus ventricosus, is currently considered a single species inhabiting the whole Mediterranean basin and the adjacent Atlantic coasts. Yet, no population genetic study has assessed its taxonomic status. Here, we collected 245 individuals from 75 localities throughout the Mediterranean Sea and used cox1 barcodes, complete mitochondrial genomes, and genome skims to test whether L. ventricosus represents a complex of cryptic species. The maximum likelihood phylogeny based on complete mitochondrial genomes recovered six main clades (hereby named blue, brown, green, orange, red, and violet) with sufficient sequence divergence to be considered putative species. On the other hand, phylogenomic analyses based on 437 nuclear genes only recovered four out of the six clades: blue and orange clades were thoroughly mixed and the brown one was not recovered. This mito-nuclear discordance revealed instances of incomplete lineage sorting and introgression, and may have caused important differences in the dating of main cladogenetic events. Species delimitation tests proposed the existence of at least three species: green, violet, and red+blue+orange (i.e., cyan). Green plus cyan (with sympatric distributions) and violet, had West and East Mediterranean distributions, respectively, mostly separated by the Siculo-Tunisian biogeographical barrier. Morphometric analyses of the shell using species hypotheses as factor and shell length as covariate showed that the discrimination power of the studied parameters was only 70.2%, reinforcing the cryptic nature of the uncovered species, and the importance of integrative taxonomic approaches considering morphology, ecology, biogeography, and mitochondrial and nuclear population genetic variation.

Phylogenomics of Anguis and Pseudopus (Squamata, Anguidae) indicates Balkan-Apennine mitochondrial capture associated with the Messinian event

A dated phylogenetic hypothesis on the evolutionary history of the extant taxa of the Western Palearctic lizards Anguis and Pseudopus is revised using genome-wide nuclear DNA and mitogenomes. We found overall concordance between nuclear and mitochondrial DNA phylogenies, with one significant exception - the Apennine A. veronensis. In mitochondrial DNA, this species forms a common clade with the earliest diverging lineage, the southern Balkan endemic A. cephallonica, while it clusters together with A. fragilis in nuclear DNA. The nuclear phylogeny conforms to the morphology, which is relatively similar between A. veronensis and A. fragilis. The most plausible explanation for the mitonuclear discordance is ancient mitochondrial capture from the Balkan ancestor of A. cephallonica to the Apennine population of the A. fragilis-veronensis ancestor. We hypothesize that this capture occurred only in a geographically restricted population. The dating of this presumed mitochondrial introgression and capture coincides with the Messinian event, when the Balkan and Apennine Peninsulas were presumably largely connected. The dated nuclear phylogenomic reconstruction estimated the divergence of A. cephallonica around 12 Mya, while the sister clade representing the A. fragilis species complex consisting of the sister species A. fragilis-A. veronensis and A. colchica-A. graeca further diversified around 7 Mya. The depth of nuclear divergence among the evolutionary lineages of Pseudopus (0.5-1.2 Mya) supports their subspecies status.

Performance of intron 7 of the β-fibrinogen gene for phylogenetic analysis: An example using gladiator frogs, Boana Gray, 1825 (Anura, Hylidae, Cophomantinae)

Boana, the third largest genus of Hylinae, has cryptic morphological species. The potential applicability of b-fibrinogen intron 7 – FGBI7 is explored to propose a robust phylogeny of Boana. The phylogenetic potential of FGBI7 was evaluated using maximum parsimony, MrBayes, and maximum likelihood analysis. Comparison of polymorphic sites and topologies obtained with concatenated analysis of FGBI7 and other nuclear genes (CXCR4, CXCR4, RHO, SIAH1, TYR, and 28S) allowed evaluation of the phylogenetic signal of FGBI7. Mean evolutionary rates were calculated using the sequences of the mitochondrial genes ND1 and CYTB available for Boana in GenBank. Dating of Boana and some of its groups was performed using the RelTime method with secondary calibration. FGBI7 analysis revealed high values at informative sites for parsimony. The absolute values of the mean evolutionary rate were higher for mitochondrial genes than for FGBI7. Dating of congruent Boana groups for ND1, CYTB, and FGBI7 revealed closer values between mitochondrial genes and slightly different values from those of FGBI7. Divergence times of basal groups tended to be overestimated when mtDNA was used and were more accurate when nDNA was used. Although there is evidence of phylogenetic potential arising from concatenation of specific genes, FGBI7 provides well-resolved independent gene trees. These results lead to a paradigm for linking data in phylogenomics that focuses on the uniqueness of species histories and ignores the multiplicities of individual gene histories.

Molecular phylogeny and taxonomy of the genus Nectogale (Mammalia: Eulipotyphla: Soricidae)

The elegant water shrew, Nectogale elegans, is one of the small mammal species most adapted to a semi-aquatic lifestyle. The taxonomy of the genus Nectogale has received little attention due to difficulties in specimen collection. In this study, we sequenced one mitochondrial and eight nuclear genes to infer the phylogenetic relationship of Nectogale. Phylogenetic analyses revealed two large clades within Nectogale. One clade represented N. elegans, and the other was regarded as N. sikhimensis. The split between N. elegans and N. sikhimensis dated back to the early Pleistocene (2.15 million years ago [Ma]), which might be relevant to the Qinghai-Tibet Plateau (QTP) uplift. The morphological comparison showed several distinguishing characters within Nectogale: the shape of the mastoids, the first lower unicuspid (a1), and the second upper molar (M2). Overall, the molecular and the morphological evidences supported that the genus Nectogale consists of two valid species: N. elegans and N. sikhimensis.

Comprehensive taxon sampling and vetted fossils help clarify the time tree of shorebirds (Aves, Charadriiformes)

Shorebirds (Charadriiformes) are a globally distributed clade of modern birds and, due to their ecological and morphological disparity, a frequent subject of comparative studies. While molecular phylogenies have been key to establishing the suprafamilial backbone of the charadriiform tree, a number of relationships at both deep and shallow taxonomic levels remain poorly resolved. The timescale of shorebird evolution also remains uncertain as a result of extensive disagreements among the published divergence dating studies, stemming largely from different choices of fossil calibrations. Here, we present the most comprehensive non-supertree phylogeny of shorebirds to date, based on a total-evidence dataset comprising 353 ingroup taxa (90% of all extant or recently extinct species), 27 loci (15 mitochondrial and 12 nuclear), and 69 morphological characters. We further clarify the timeline of charadriiform evolution by time-scaling this phylogeny using a set of 14 up-to-date and thoroughly vetted fossil calibrations. In addition, we assemble a taxonomically restricted 100-locus dataset specifically designed to resolve outstanding problems in higher-level charadriiform phylogeny. In terms of tree topology, our results are largely congruent with previous studies but indicate that some of the conflicts among earlier analyses reflect a genuine signal of pervasive gene tree discordance. Monophyly of the plovers (Charadriidae), the position of the ibisbill (Ibidorhyncha), and the relationships among the five subfamilies of the gulls (Laridae) could not be resolved even with greatly increased locus and taxon sampling. Moreover, several localized regions of uncertainty persist in shallower parts of the tree, including the interrelationships of the true auks (Alcinae) and anarhynchine plovers. Our node-dating and macroevolutionary rate analyses find support for a Paleocene origin of crown-group shorebirds, as well as exceptionally rapid recent radiations of Old World oystercatchers (Haematopodidae) and select genera of gulls. Our study underscores the challenges involved in estimating a comprehensively sampled and carefully calibrated time tree for a diverse avian clade, and highlights areas in need of further research.

First description of the mitogenome and phylogeny:Aedes vexansand Ochlerotatus caspius of the Tribe Aedini (Diptera: Culicidae)

Culicidae, the mosquito family, includes more than 3600 species subdivided into the subfamilies Anophelinae and Culicinae. One-third of mosquitoes belong to the Aedini tribe, which is subordinate to the subfamily Culicinae, which comprises common vectors of viral zoonoses. The tribe of Aedini is extremely diverse in morphology and geographical distribution and has high ecological and medical significance. However, knowledge about the systematics of the Aedini tribe is still limited owing to its large population and the similar morphological characteristics of its species. This study provides the first description of the complete mitochondrial (mt) genome sequence of Aedes vexans and Ochlerotatus caspius belonging to the Aedini tribe. The mt genomes of A. vexans and O. caspius are circular molecules that are 15,861 bp and 15,954 bp in size, with AT contents of 78.54% and 79.36%, respectively. Both the circular mt genomes comprise 37 functional subunits, including 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes (tRNAs), and a control region (also known as the AT-rich region). The most common start codons are ATT/ATG, apart from cox1 (TCG) and nad5 (GTG), while TAA is the termination codon for all PCGs. All tRNAs have a typical clover leaf structure, except tRNA Ser1. Phylogenetic analysis of the concatenated, aligned amino acid sequences of the 13 PCGs showed that A. vexans gathered with Aedes sp. in a sister taxon, and O. caspius gathered with Ochlerotatus sp. in a sister taxon. The findings from the present study support the concept of monophyly of all groups, ratify the current taxonomic classification, and provide vital molecular marker resources for further studies of the taxonomy, population genetics, and systematics of the Aedini tribe.

Middle Jurassic fossils document an early stage in salamander evolution

Salamanders are an important group of living amphibians and model organisms for understanding locomotion, development, regeneration, feeding, and toxicity in tetrapods. However, their origin and early radiation remain poorly understood, with early fossil stem-salamanders so far represented by larval or incompletely known taxa. This poor record also limits understanding of the origin of Lissamphibia (i.e., frogs, salamanders, and caecilians). We report fossils from the Middle Jurassic of Scotland representing almost the entire skeleton of the enigmatic stem-salamander Marmorerpeton. We use computed tomography to visualize high-resolution three-dimensional anatomy, describing morphologies that were poorly characterized in early salamanders, including the braincase, scapulocoracoid, and lower jaw. We use these data in the context of a phylogenetic analysis intended to resolve the relationships of early and stem-salamanders, including representation of important outgroups alongside data from high-resolution imaging of extant species. Marmorerpeton is united with Karaurus, Kokartus, and others from the Middle Jurassic-Lower Cretaceous of Asia, providing evidence for an early radiation of robustly built neotenous stem-salamanders. These taxa display morphological specializations similar to the extant cryptobranchid "giant" salamanders. Our analysis also demonstrates stem-group affinities for a larger sample of Jurassic species than previously recognized, highlighting an unappreciated diversity of stem-salamanders and cautioning against the use of single species (e.g., Karaurus) as exemplars for stem-salamander anatomy. These phylogenetic findings, combined with knowledge of the near-complete skeletal anatomy of Mamorerpeton, advance our understanding of evolutionary changes on the salamander stem-lineage and provide important data on early salamanders and the origins of Batrachia and Lissamphibia.

Mitochondrial characteristics of Pseudohynobius flavomaculatus a protected salamander in China, and biogeographical implications for the family Hynobiidae (Amphibia, Caudata)

Pseudohynobius flavomaculatus a provincially-protected salamander species, inhabits mountainous areas of Chongqing and surrounding provinces in China. In the present study, the complete mitochondrial genome of P. flavomaculatus was sequenced and analyzed. The mitogenome is 16,401 bp in length and consisted of 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a control region. We performed a novel phylogenetic analysis, which demonstrated a sister relationship between P. flavomaculatus and P. jinfo. The 95% confidence interval around our new divergence date estimate suggest that Hynobiidae originated at 101.62–119.84 (mean=110.87) Ma. Species within Hynobiidae diverged successively in the Cenozoic era, and hynobiid speciation coincides primarily with geologic events. Our biogeographical inference demonstrates that nearly all early hynobiids divergences correspond to geological estimates of orogeny, which may have contributed to the notably high dN/dS ratio in this clade. We conclude that orogeny is likely a primary, dynamic factor, which may have repeatedly initiated the process of speciation in the family Hynobiidae.

The distribution of the critically endangered salamander Paradactylodon (Afghanodon) mustersi (Smith, 1940) in Afghanistan

The Afghanistan Mountain Salamander, Paradactylodon (Afghanodon) mustersi (Smith, 1940), is an evolutionary old species, listed as Critically Endangered by IUCN Red List. The species is endemic to the Hindu Kush mountain range with records from only a few localities from four central-eastern Afghan provinces (Kabul, Ghazni, Parwan, and Wardak). Due to the long-term complicated political situation in the country which makes zoological research almost impossible, the current distribution and the presence of P. mustersi at previously known localities has remained unassessed for 40 years. We carried out recent, sporadic surveys between 2017 and 2021 to detect P. mustersi in three tributaries of the Paghman stream on the southern slopes of the Hindu Kush Mountains and the two nearby provinces (Panjsheer and Parwan), where the presence of the species was expected. We confirmed the occurrence of P. mustersi at all survey sites, and present the first record of the species for Panjsheer Province. We also confirmed that the species is currently endangered by human-mediated factors like habitat disturbance (increased visitors attendance, water pollution, construction activities), especially in the Paghman area.

Positive selection in multiple salivary gland proteins of Anophelinae reveals potential targets for vector control

Anopheles is a genus belonging to the Culicidae family, which has great medical importance due to its role as a vector of Plasmodium, the causative agent of malaria. Great focus has been given to the salivary gland proteins (SGPs) group from Anopheles' functional genomics. This class of proteins is essential to blood-feeding behavior as they have attributes such as vasodilators and anti-clotting properties. Recently, a comprehensive review on Anopheles SGPs was performed; however, the authors did not deeply explore the adaptive molecular evolution of these genes. In this context, this work aimed to perform a more detailed analysis of the adaptive molecular evolution of SGPs in Anopheles, carrying out positive selection and gene family evolution analysis on 824 SGPs. Our results show that most SGPs have positively selected codon sites that can be used as targets in developing new strategies for vector control and that younger SGPs evolve at a faster rate than older SGPs. Notably, we could not find any evidence of an accelerated shift in SGPs' rates of gene gain and loss compared with other proteins, as suggested in previous works.

Invasions of an obligate asexual daphnid species support the nearly neutral theory

To verify the "nearly neutral theory (NNT)," the ratio of nonsynonymous to synonymous substitutions (dN/dS) was compared among populations of different species. To determine the validity of NNT, however, populations that are genetically isolated from each other but share the same selection agents and differ in size should be compared. Genetically different lineages of obligate asexual Daphnia pulex invading Japan from North America are an ideal example as they satisfy these prerequisites. Therefore, we analyzed the whole-genome sequences of 18 genotypes, including those of the two independently invaded D. pulex lineages (JPN1 and JPN2) and compared the dN/dS ratio between the lineages. The base substitution rate of each genotype demonstrated that the JPN1 lineage having a larger distribution range diverged earlier and thus was older than the JPN2 lineage. Comparisons of the genotypes within lineages revealed that changes in dN/dS occurred after the divergence and were larger in the younger lineage, JPN2. These results imply that the JPN1 lineage has been more effectively subjected to purification selections, while slightly deteriorating mutations are less purged in JPN2 with smaller population size. Altogether, the lineage-specific difference in the dN/dS ratio for the obligate asexual D. pulex was well explained by the NNT.

Molecular dating of the blood pigment hemocyanin provides new insight into the origin of animals

The Neoproterozoic included changes in oceanic redox conditions, the configuration of continents and climate, extreme ice ages (Sturtian and Marinoan), and the rise of complex life forms. A much-debated topic in geobiology concerns the influence of atmospheric oxygenation on Earth and the origin and diversification of animal lineages, with the most widely popularized hypotheses relying on causal links between oxygen levels and the rise of animals. The vast majority of extant animals use aerobic metabolism for growth and homeostasis; hence, the binding and transportation of oxygen represent a vital physiological task. Considering the blood pigment hemocyanin (Hc) is present in sponges and ctenophores, and likely to be present in the common ancestor of animals, we investigated the evolution and date of Hc emergence using bioinformatics approaches on both transcriptomic and genomic data. Bayesian molecular dating suggested that the ancestral animal Hc gene arose approximately 881 Ma during the Tonian Period (1000-720 Ma), prior to the extreme glaciation events of the Cryogenian Period (720-635 Ma). This result is corroborated by a recently discovered fossil of a putative sponge ~890 Ma and modern molecular dating for the origin of metazoans of ~1,000-650 Ma (but does contradict previous inferences regarding the origin of Hc ~700-600 Ma). Our data reveal that crown-group animals already possessed hemocyanin-like blood pigments, which may have enhanced the oxygen-carrying capacity of these animals in hypoxic environments at that time or acted in the transport of hormones, detoxification of heavy metals, and immunity pathways.

Genome-scale data reveal deep lineage divergence and a complex demographic history in the Texas horned lizard (Phrynosoma cornutum) throughout the southwestern and central US

The southwestern and central United States serve as an ideal region to test alternative hypotheses regarding biotic diversification. Genomic data can now be combined with sophisticated computational models to quantify the impacts of paleoclimate change, geographic features, and habitat heterogeneity on spatial patterns of genetic diversity. In this study, we combine thousands of genotyping-by-sequencing (GBS) loci with mtDNA sequences (ND1) from the Texas horned lizard (Phrynosoma cornutum) to quantify relative support for different catalysts of diversification. Phylogenetic and clustering analyses of the GBS data indicate support for at least three primary populations. The spatial distribution of populations appears concordant with habitat type, with desert populations in AZ and NM showing the largest genetic divergence from the remaining populations. The mtDNA data also support a divergent desert population, but other relationships differ and suggest mtDNA introgression. Genotype–environment association with bioclimatic variables supports divergence along precipitation gradients more than along temperature gradients. Demographic analyses support a complex history, with introgression and gene flow playing an important role during diversification. Bayesian multispecies coalescent analyses with introgression (MSci) analyses also suggest that gene flow occurred between populations. Paleo-species distribution models support two southern refugia that geographically correspond to contemporary lineages. We find that divergence times are underestimated and population sizes are overestimated when introgression occurred and is ignored in coalescent analyses, and furthermore, inference of ancient introgression events and demographic history is sensitive to inclusion of a single recently admixed sample. Our analyses cannot refute the riverine barrier or glacial refugia hypotheses. Results also suggest that populations are continuing to diverge along habitat gradients. Finally, the strong evidence of admixture, gene flow, and mtDNA introgression among populations suggests that P. cornutum should be considered a single widespread species under the General Lineage Species Concept.

Sympatric and independently evolving lineages in the Thoropa miliaris - T. taophora species complex (Anura: Cycloramphidae)

Species delimitation can be challenging and affected by subjectivity. Sibling lineages that occur in sympatry constitute good candidates for species delimitation regardless of the adopted species concept. The Thoropa miliaris + T. taophora species complex exhibits high genetic diversity distributed in several lineages that occur sympatrically in the southeastern Atlantic Forest of Brazil. We used 414 loci obtained by anchored hybrid enrichment to characterize genetic variation in the Thoropa miliaris species group (T. saxatilis, T megatympanum, T. miliaris, and T. taophora), combining assignment analyses with traditional and coalescent phylogeny reconstruction. We also investigated evolutionary independence in co-occurring lineages by estimating gene flow, and validated lineages under the multispecies coalescent. We recovered most previously described lineages as unique populations in assignment analyses; exceptions include two lineages within T. miliaris that are further substructured, and the merging of all T. taophora lineages. We found very low probabilities of gene flow between sympatric lineages, suggesting independent evolution. Species tree inferences and species delimitation yielded resolved relationships and indicate that all lineages constitute putative species that diverged during the Pliocene and Pleistocene, later than previously estimated.

Evolution of a key enzyme of aerobic metabolism reveals Proterozoic functional subunit duplication events and an ancient origin of animals

The biological toolkits for aerobic respiration were critical for the rise and diversification of early animals. Aerobic life forms generate ATP through the oxidation of organic molecules in a process known as Krebs' Cycle, where the enzyme isocitrate dehydrogenase (IDH) regulates the cycle's turnover rate. Evolutionary reconstructions and molecular dating of proteins related to oxidative metabolism, such as IDH, can therefore provide an estimate of when the diversification of major taxa occurred, and their coevolution with the oxidative state of oceans and atmosphere. To establish the evolutionary history and divergence time of NAD-dependent IDH, we examined transcriptomic data from 195 eukaryotes (mostly animals). We demonstrate that two duplication events occurred in the evolutionary history of NAD-IDH, one in the ancestor of eukaryotes approximately at 1967 Ma, and another at 1629 Ma, both in the Paleoproterozoic Era. Moreover, NAD-IDH regulatory subunits β and γ are exclusive to metazoans, arising in the Mesoproterozoic. Our results therefore support the concept of an ''earlier-than-Tonian'' diversification of eukaryotes and the pre-Cryogenian emergence of a metazoan IDH enzyme.

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.

Global Diversity and Geographic Distributions of Padina Species (Dictyotales, Phaeophyceae): New Insights Based on Molecular and Morphological Analyses

The taxonomic status and species diversity of the brown algal genus Padina (Dictyotales, Phaeophyceae) was assessed based on DNA sequences and the morpho-anatomy of specimens collected worldwide, especially from tropical and subtropical western Pacific regions. Phylogenetic analyses using chloroplast rbcL and mitochondrial cox3 gene sequences demonstrated four distinct clades for newly collected samples with high bootstrap support. Each species clade possesses a suite of morphological features that are not shared by any known species of Padina. These are P. imbricata sp. nov., Padina lutea sp. nov., P. moffittianoides sp. nov., and P. nitida sp. nov. The occurrence of these and other species of Padina clearly points to an elevated diversity of the genus in tropical/subtropical waters of the western Pacific. Phylogenetic analyses provided new insights into biogeographic characteristics of the genus, with many species in the Pacific Ocean showing shared/overlapping distributions, whereas species from the Mediterranean/Atlantic and/or the Indian Ocean tend to be confined to particular regions. Consideration has also been given to the evolutionary time frame of the genus Padina based on molecular time trees: a time tree of the concatenated data set (rbcL + cox3) revealed the estimated divergence time in the mid-Cretaceous, whereas those of cox3 and rbcL showed older estimates pointing to the periods of mid-Jurassic and Early Cretaceous, respectively.

Impact of the Miocene orogenesis on Kaloula spp. radiation and implication of local refugia on genetic diversification

The phylogeography of the Kaloula genus in East Asia is still poorly understood. One of the difficulties is the absence of fossils to corroborate molecular dating estimates. Here, we examined the mitochondrial structure of Kaloula spp. in East Asia and focused on the impact of glaciations on the northernmost species: Kaloula borealis. We determined the phylogenetic relationships, molecular dating, and genetic connectivity assessments within the genus from 1211 bp of concatenated mitochondrial 12S and 16S. The relaxed clock analyses reveal the emergence of Kaloula spp. common ancestor in East and Southeast Asia between the Eocene and Oligocene, c. 38.47 Ma (24.69-53.65). The genetic diversification of lineages then increased on the East Asian Mainland during the Lower Miocene, c. 20.10 (8.73-30.65), most likely originating from the vicariance and radiation triggered by the orogeny of the Qinghai-Tibetan Plateau. Later, the dispersal towards the North East Asian Mainland during the Upper Miocene drove the population diversification of K. borealis c. 9.01 Ma (3.66-15.29). Finally, the central mainland population became isolated following orogenesis events and diverged into K. rugifera during the Pliocene, c. 3.06 Ma (0.02-10.90). The combination of population genetic and barrier analyses revealed a significant genetic isolation between populations of Kaloula spp. matching with the massive Qinling-Daba Mountain chain located in south-central China. Finally, we highlight a young divergence within the Eastern Mainland population of K. borealis, possibly attributed to refugia in south eastern China from which populations later expanded.

African lates perches (Teleostei, Latidae, Lates): Paraphyly of Nile perch and recent colonization of Lake Tanganyika

Lates perches of the genus Lates (Latidae) are large piscivorous fishes, with a strikingly disjunct distribution range in coastal areas and estuaries of the Indo-Pacific region and in some large African freshwater systems. Previous phylogenetic hypotheses based on osteological and ontogenetic data suggested paraphyly of the African representatives, or even the small Lake Tanganyika species assemblage, with respect to the remaining Lates species. Based on a multilocus phylogeny, however, we show that extant African lates perches are monophyletic. The Nile perch, L. niloticus, which is widely distributed in the Nilo-Sudan region and Central Africa, comprises three distinct lineages and is paraphyletic with respect to the four endemic Lake Tanganyika species. We find that diversification of extant African Lates happened only as recently as the Pliocene. With the extensive, in part much older fossil record, this suggests repeated extinction and (re-)colonization of hydrological systems. We further find that Lates started to diversify in Lake Tanganyika only in the Pleistocene, which is much more recent than other fish radiations endemic to Lake Tanganyika, implying that they radiated in the presence of other top predators already in this ecosystem.

Chronological Incongruences between Mitochondrial and Nuclear Phylogenies of Aedes Mosquitoes

One-third of all mosquitoes belong to the Aedini, a tribe comprising common vectors of viral zoonoses such as Aedes aegypti and Aedes albopictus. To improve our understanding of their evolution, we present an updated multigene estimate of Aedini phylogeny and divergence, focusing on the disentanglement between nuclear and mitochondrial phylogenetic signals. We first show that there are some phylogenetic discrepancies between nuclear and mitochondrial markers which may be caused by wrong taxa assignment in samples collections or by some stochastic effect due to small gene samples. We indeed show that the concatenated dataset is model and framework dependent, indicating a general paucity of signal. Our Bayesian calibrated divergence estimates point toward a mosquito radiation in the mid-Jurassic and an Aedes radiation from the mid-Cretaceous on. We observe, however a strong chronological incongruence between mitochondrial and nuclear data, the latter providing divergence times within the Aedini significantly younger than the former. We show that this incongruence is consistent over different datasets and taxon sampling and that may be explained by either peculiar evolutionary event such as different levels of saturation in certain lineages or a past history of hybridization throughout the genus. Overall, our updated picture of Aedini phylogeny, reveal a strong nuclear-mitochondrial incongruence which may be of help in setting the research agenda for future phylogenomic studies of Aedini mosquitoes.

Comparative osteology of the hynobiid complex Liua-Protohynobius-Pseudohynobius (Amphibia, Urodela): Ⅰ. Cranial anatomy of Pseudohynobius

Hynobiidae are a clade of salamanders that diverged early within the crown radiation and that retain a considerable number of features plesiomorphic for the group. Their evolutionary history is informed by a fossil record that extends to the Middle Jurassic Bathonian time. Our understanding of the evolution within the total group of Hynobiidae has benefited considerably from recent discoveries of stem hynobiids but is constrained by inadequate anatomical knowledge of some extant forms. Pseudohynobius is a derived hynobiid clade consisting of five to seven extant species living endemic to southwestern China. Although this clade has been recognized for over 37 years, osteological details of these extant hynobiids remain elusive, which undoubtedly has contributed to taxonomic controversies over the hynobiid complex Liua-Protohynobius-Pseudohynobius. Here we provide a bone-by-bone study of the cranium in the five extant species of Pseudohynobius (Ps. flavomaculatus, Ps. guizhouensis, Ps. jinfo, Ps. kuankuoshuiensis and Ps. shuichengensis) based on x-ray computer tomography data for 18 specimens. Our results indicate that the cranium in each of these species has a combination of differences in morphology, proportions and articulation patterns in both dermal and endochondral bones. Our study establishes a range of intraspecific differences that will serve as organizing hypotheses for future studies as more extensive collections of these species become available. Morphological features in the cranium for terrestrial ecological adaptation in Hynobiidae are summarized. Based on the results, we also discuss the evolution and development of several potential synapomorphies of Hynobiidae, including features of the orbitosphenoid and articular.

A new species of the genus Batrachuperus (Urodela: Hynobiidae) from Southwestern China

A new stream salamander species, Batrachuperus daochengensissp. nov., from southwestern China, is described herein based on morphological and molecular evidence. Molecular phylogeny derived from the mitochondrial gene together with previous nuclear data revealed that B. daochengensissp. nov. is sister to B. yenyuanensis. The new species differs from all other species of the genus by the following combination of characters: brown horny epidermis on tips of fingers and toes absent; tubercles on palms and soles absent; costal grooves 12; dorsal brown, mottled with blackish spots; fingers 2-3-4-1 in order of decreasing length; tips of longest digits of fore- and hindlimbs largely separated by one to two costal spaces when adpressed towards each other along sides of body. The new species is currently known in the central and southern Shaluli Mountains in southwestern China.

Complete mitochondrial genome of Hynobius dunni (Amphibia, Caudata, Hynobiidae) and its phylogenetic position

Hynobius dunni is a salamander species of the genus Hynobius endemically distributed in eastern Kyushu in southwestern Japan. In this study, we determined the complete mitochondrial genome sequence and clarified the phylogenetic position of this species. The mitochondrial genome was 16,47 bp in length and encoded 13 protein, 2 ribosomal RNA, and 22 transfer RNA genes. Phylogenetic tree based on 13 protein-coding genes revealed that H. nebulosus were the most closely related species within the Hynobius species. The data identified in this study will be useful for population and conservation genetic studies of Hynobius species.

Climatic and topographic changes since the Miocene influenced the diversification and biogeography of the tent tortoise (Psammobates tentorius) species complex in Southern Africa

BACKGROUND

Climatic and topographic changes function as key drivers in shaping genetic structure and cladogenic radiation in many organisms. Southern Africa has an exceptionally diverse tortoise fauna, harbouring one-third of the world's tortoise genera. The distribution of Psammobates tentorius (Kuhl, 1820) covers two of the 25 biodiversity hotspots in the world, the Succulent Karoo and Cape Floristic Region. The highly diverged P. tentorius represents an excellent model species for exploring biogeographic and radiation patterns of reptiles in Southern Africa.

RESULTS

We investigated genetic structure and radiation patterns against temporal and spatial dimensions since the Miocene in the Psammobates tentorius species complex, using multiple types of DNA markers and niche modelling analyses. Cladogenesis in P. tentorius started in the late Miocene (11.63-5.33 Ma) when populations dispersed from north to south to form two geographically isolated groups. The northern group diverged into a clade north of the Orange River (OR), followed by the splitting of the group south of the OR into a western and an interior clade. The latter divergence corresponded to the intensification of the cold Benguela current, which caused western aridification and rainfall seasonality. In the south, tectonic uplift and subsequent exhumation, together with climatic fluctuations seemed responsible for radiations among the four southern clades since the late Miocene. We found that each clade occurred in a habitat shaped by different climatic parameters, and that the niches differed substantially among the clades of the northern group but were similar among clades of the southern group.

CONCLUSION

Climatic shifts, and biome and geographic changes were possibly the three major driving forces shaping cladogenesis and genetic structure in Southern African tortoise species. Our results revealed that the cladogenesis of the P. tentorius species complex was probably shaped by environmental cooling, biome shifts and topographic uplift in Southern Africa since the late Miocene. The Last Glacial Maximum (LGM) may have impacted the distribution of P. tentorius substantially. We found the taxonomic diversify of the P. tentorius species complex to be highest in the Greater Cape Floristic Region. All seven clades discovered warrant conservation attention, particularly Ptt-B-Ptr, Ptt-A and Pv-A.

A large-scale systematic framework of Chinese snakes based on a unified multilocus marker system

Snakes are one of the most diverse groups of terrestrial vertebrates, with approximately 3500 extant species. A robust phylogeny and taxonomy of snakes is crucial for us to know, study and protect them. For a large group such as snakes, broad-scale phylogenetic reconstructions largely rely on data integration. Increasing the compatibility of the data from different researches is thus important, which can be facilitated by standardization of the loci used in systematic analyses. In this study, we proposed a unified multilocus marker system for snake systematics by conflating 5 mitochondrial markers, 19 vertebrate-universal nuclear protein coding (NPC) markers and 72 snake-specific noncoding intron markers. This marker system is an addition to the large squamate conserved locus set (SqCL) for studies preferring a medium-scale data set. We applied this marker system to over 440 snake samples and constructed the currently most comprehensive systematic framework of the snakes in China. Robust snake phylogenetic relationships were recovered at both deep and shallow evolutionary depths, demonstrating the usefulness of this multilocus marker system. Discordance was revealed by a parallel comparison between the snake tree based on the multilocus marker system and that based on only the mitochondrial loci, highlighting the necessity of using multiple types of markers to better understand the snake evolutionary histories. The divergence times of different snake groups were estimated with the nuclear data set. Our comprehensive snake tree not only confirms many important nodes inferred in previous studies but also contributes new insights into many snake phylogenetic relationships. Suggestions are made for the current Chinese snake taxonomy.

Mitogenomic Perspectives on the Adaptation to Extreme Alkaline Environment of Amur ide (Leuciscus waleckii)

Amur ide (Leuciscus waleckii, Family Cyprinidae) is widely distributed in Northeast Asia. L. waleckii usually inhabits freshwater environments but can also survive in the Lake Dali Nur, one of the most extreme aquatic environments on the earth, with an alkalinity up to 50 mmol/L (pH 9.6). To investigate mechanisms of mitogenomic evolution underlying adaptation to extreme environments, we determined 30 complete mitogenomes that included Lake Dali Nur (alkaline environment, AL) population and Amur basin (freshwater environment, FW) population. Through phylogenetic and divergence time analysis, we found that AL and FW populations forming distinct two groups which were consistent with geographic divergence (the formation of Lake Dali Nur). In addition, we found that almost of the windows exhibited higher nucleotide diversity in FW population (avg 0.0046) than AL population (avg 0.0012). This result indicated that severe environment selection had remarkably reduced the genetic diversity of mitogenome in AL population and suggested that severe environment selection had remarkably reduced the genetic diversity of mitogenome in the AL population. Compared with the FW population (ω = 0.064), the AL population (ω = 0.092) had a larger mean ω (dN/dS ratios) value for the 13 concatenated mitochondrial protein-coding genes, indicating that the high alkaline tolerated group had accumulated more nonsynonymous mutations. These nonsynonymous mutations had resulted in slightly beneficial amino acid changes that allowed adaption to the severe conditions. This study provides an additional view to decipher the adaptive mitogenome evolution of L. waleckii of the high alkaline environment.

Phylogeny of spiny frogs Nanorana (Anura: Dicroglossidae) supports a Tibetan origin of a Himalayan species group

Recent advances in the understanding of the evolution of the Asian continent challenge the long-held belief of a faunal immigration into the Himalaya. Spiny frogs of the genus Nanorana are a characteristic faunal group of the Himalaya-Tibet orogen (HTO). We examine the phylogeny of these frogs to explore alternative biogeographic scenarios for their origin in the Greater Himalaya, namely, immigration, South Tibetan origin, strict vicariance. We sequenced 150 Nanorana samples from 62 localities for three mitochondrial (1,524 bp) and three nuclear markers (2,043 bp) and complemented the data with sequence data available from GenBank. We reconstructed a gene tree, phylogenetic networks, and ancestral areas. Based on the nuDNA, we also generated a time-calibrated species tree. The results revealed two major clades (Nanorana and Quasipaa), which originated in the Lower Miocene from eastern China and subsequently spread into the HTO (Nanorana). Five well-supported subclades are found within Nanorana: from the East, Central, and Northwest Himalaya, the Tibetan Plateau, and the southeastern Plateau margin. The latter subclade represents the most basal group (subgenus Chaparana), the Plateau group (Nanorana) represents the sister clade to all species of the Greater Himalaya (Paa). We found no evidence for an east-west range expansion of Paa along the Himalaya, nor clear support for a strict vicariance model. Diversification in each of the three Himalayan subclades has probably occurred in distinct areas. Specimens from the NW Himalaya are placed basally relative to the highly diverse Central Himalayan group, while the lineage from the Tibetan Plateau is placed within a more terminal clade. Our data indicate a Tibetan origin of Himalayan Nanorana and support a previous hypothesis, which implies that a significant part of the Himalayan biodiversity results from primary diversification of the species groups in South Tibet before this part of the HTO was uplifted to its recent heights.

A High-Quality Assembly of the Nine-Spined Stickleback (Pungitius pungitius) Genome

The Gasterosteidae fish family hosts several species that are important models for eco-evolutionary, genetic, and genomic research. In particular, a wealth of genetic and genomic data has been generated for the three-spined stickleback (Gasterosteus aculeatus), the "ecology's supermodel," whereas the genomic resources for the nine-spined stickleback (Pungitius pungitius) have remained relatively scarce. Here, we report a high-quality chromosome-level genome assembly of P. pungitius consisting of 5,303 contigs (N50 = 1.2 Mbp) with a total size of 521 Mbp. These contigs were mapped to 21 linkage groups using a high-density linkage map, yielding a final assembly with 98.5% BUSCO completeness. A total of 25,062 protein-coding genes were annotated, and about 23% of the assembly was found to consist of repetitive elements. A comprehensive analysis of repetitive elements uncovered centromere-specific tandem repeats and provided insights into the evolution of retrotransposons. A multigene phylogenetic analysis inferred a divergence time of about 26 million years ago (Ma) between nine- and three-spined sticklebacks, which is far older than the commonly assumed estimate of 13 Ma. Compared with the three-spined stickleback, we identified an additional duplication of several genes in the hemoglobin cluster. Sequencing data from populations adapted to different environments indicated potential copy number variations in hemoglobin genes. Furthermore, genome-wide synteny comparisons between three- and nine-spined sticklebacks identified chromosomal rearrangements underlying the karyotypic differences between the two species. The high-quality chromosome-scale assembly of the nine-spined stickleback genome obtained with long-read sequencing technology provides a crucial resource for comparative and population genomic investigations of stickleback fishes and teleosts.

A High-Quality Assembly of the Nine-Spined Stickleback (Pungitius pungitius) Genome

The Gasterosteidae fish family hosts several species that are important models for eco-evolutionary, genetic, and genomic research. In particular, a wealth of genetic and genomic data has been generated for the three-spined stickleback (Gasterosteus aculeatus), the "ecology's supermodel," whereas the genomic resources for the nine-spined stickleback (Pungitius pungitius) have remained relatively scarce. Here, we report a high-quality chromosome-level genome assembly of P. pungitius consisting of 5,303 contigs (N50 = 1.2 Mbp) with a total size of 521 Mbp. These contigs were mapped to 21 linkage groups using a high-density linkage map, yielding a final assembly with 98.5% BUSCO completeness. A total of 25,062 protein-coding genes were annotated, and about 23% of the assembly was found to consist of repetitive elements. A comprehensive analysis of repetitive elements uncovered centromere-specific tandem repeats and provided insights into the evolution of retrotransposons. A multigene phylogenetic analysis inferred a divergence time of about 26 million years ago (Ma) between nine- and three-spined sticklebacks, which is far older than the commonly assumed estimate of 13 Ma. Compared with the three-spined stickleback, we identified an additional duplication of several genes in the hemoglobin cluster. Sequencing data from populations adapted to different environments indicated potential copy number variations in hemoglobin genes. Furthermore, genome-wide synteny comparisons between three- and nine-spined sticklebacks identified chromosomal rearrangements underlying the karyotypic differences between the two species. The high-quality chromosome-scale assembly of the nine-spined stickleback genome obtained with long-read sequencing technology provides a crucial resource for comparative and population genomic investigations of stickleback fishes and teleosts.

Historical museum collections clarify the evolutionary history of cryptic species radiation in the world's largest amphibians

Inaccurate taxonomic assessment of threatened populations can hinder conservation prioritization and management, with human-mediated population movements obscuring biogeographic patterns and confounding reconstructions of evolutionary history. Giant salamanders were formerly distributed widely across China, and are interpreted as a single species, Andrias davidianus. Previous phylogenetic studies have identified distinct Chinese giant salamander lineages but were unable to associate these consistently with different landscapes, probably because population structure has been modified by human-mediated translocations for recent commercial farming. We investigated the evolutionary history and relationships of allopatric Chinese giant salamander populations with Next-Generation Sequencing methods, using historical museum specimens and late 20th-century samples, and retrieved partial or near-complete mitogenomes for 17 individuals. Samples from populations unlikely to have been affected by translocations form three clades from separate regions of China, spatially congruent with isolation by either major river drainages or mountain ranges. Pliocene-Pleistocene divergences for these clades are consistent with topographic modification of southern China associated with uplift of the Qinghai-Tibet Plateau. General Mixed Yule Coalescent model analysis indicates that these clades represent separate species: Andrias davidianus (Blanchard, 1871) (northern Yangtze/Sichuan), Andrias sligoi (Boulenger, 1924) (Pearl/Nanling), and an undescribed species (Huangshan). Andrias sligoi is possibly the world's largest amphibian. Inclusion of additional reportedly wild samples from areas of known giant salamander exploitation and movement leads to increasing loss of biogeographic signal. Wild Chinese giant salamander populations are now critically depleted or extirpated, and conservation actions should be updated to recognize the existence of multiple species.

Biogeography, Systematics, and Ecomorphology of Pacific Island Anoles

: Anoles are regarded as important models for understanding dynamic processes in ecology and evolution. Most work on this group has focused on species in the Caribbean Sea, and recently in mainland South and Central America. However, the Eastern Tropical Pacific (ETP) is home to seven species of anoles from three unique islands (Islas Cocos, Gorgona, and Malpelo) that have been largely overlooked. Four of these species are endemic to single islands (Norops townsendi on Isla Cocos, Dactyloa agassizi on Isla Malpelo, D. gorgonae and N. medemi on Isla Gorgona). Herein, we present a phylogenetic analysis of anoles from these islands in light of the greater anole phylogeny to estimate the timing of divergence from mainland lineages for each species. We find that two species of solitary anoles (D. agassizi and N. townsendi) diverged from mainland ancestors prior to the emergence of their respective islands. We also present population-wide morphological data suggesting that both display sexual size dimorphism, similar to single-island endemics in the Caribbean. All lineages on Isla Gorgona likely arose during past connections with South America, and ecologically partition their habitat. Finally, we highlight the importance of conservation of these species and island fauna in general.

Speciation, gene flow, and seasonal migration in Catharus thrushes (Aves:Turdidae)

New World thrushes in the genus Catharus are small, insectivorous or omnivorous birds that have been used to explore several important questions in avian evolution, including the evolution of seasonal migration and plumage variation. Within Catharus, members of a clade of obligate long-distance migrants (C. fuscescens, C. minimus, and C. bicknelli) have also been used in the development of heteropatric speciation theory, a divergence process in which migratory lineages (which might occur in allopatry or sympatry during portions of their annual cycle) diverge despite low levels of gene flow. However, research on Catharus relationships has thus far been restricted to the use of small genetic datasets, which provide limited resolution of both phylogenetic and demographic histories. We used a large, multi-locus dataset from loci containing ultraconserved elements (UCEs) to study the demographic histories of the migratory C. fuscescens-minimus-bicknelli clade and to resolve the phylogeny of the migratory species of Catharus. Our dataset included more than 2000 loci and over 1700 variable genotyped sites, and analyses supported our prediction of divergence with gene flow in the fully migratory clade, with significant gene flow among all three species. Our phylogeny of the genus differs from past work in its placement of C. ustulatus, and further analyses suggest historic gene flow throughout the genus, producing genetically reticulate (or network) phylogenies. This raises questions about trait origins and suggests that seasonal migration and the resulting migratory condition of heteropatry is likely to promote hybridization not only during pairwise divergence and speciation, but also among non-sisters.

Reconstructing evolution at the community level: A case study on Mediterranean amphibians

Reconstructing reliable timescales for species evolution is an important and indispensable goal of modern biogeography. However, many factors influence the estimation of divergence times, and uncertainty in the inferred time trees remains a major issue that is often insufficiently acknowledged. We here focus on a fundamental problem of time tree analysis: the combination of slow-evolving (nuclear DNA) and fast-evolving (mitochondrial DNA) markers in a single time tree. Both markers differ in their suitability to infer divergences at different time scales (the 'genome-timescale-dilemma'). However, strategies to infer shallow and deep divergences in a single time tree have rarely been compared empirically. Using Mediterranean amphibians as model system that is exceptional in its geographic and taxonomic completeness of available genetic information, we analyze 202 lineages of western Palearctic amphibians across the entire Mediterranean region. We compiled data of four nuclear and five mitochondrial genes and used twelve fossil calibration points widely acknowledged for amphibian evolution. We reconstruct time trees for an extensive lineage-level data set and compare the performances of the different trees: the first tree is based on primary fossil calibration and mitochondrial DNA, while the second tree is based on a combination of primary fossil and on secondary calibrations taken from a nuclear tree using mitochondrial DNA (two-step protocol). Focusing on a set of nodes that are most likely explained by vicariance, we statistically compare the reconstructed alternative time trees by applying a biogeographical plausibility test. Our two-step protocol outperformed the alternative approach in terms of spatial and temporal plausibility. It allows us to infer scenarios for Mediterranean amphibian evolution in eight geographic provinces. We identified several tectonic and climatic events explaining the majority of Mediterranean amphibian divergences, with Plio-Pleistocene climatic fluctuations being the dominant driver for intrageneric evolution. However, often more than one event could be invoked for a specific split. We give recommendations for the use of secondary calibrations in future molecular clock analyses at the community level.

Osteology of Batrachuperus yenyuanensis (Urodela, Hynobiidae), a high-altitude mountain stream salamander from western China

Batrachuperus yenyuanensis, commonly known as Yenyuan Stream Salamander, is a hynobiid species inhabiting high-altitude (2440-4025 m above sea level) mountain stream and pond environments along the eastern fringe of the Qinghai-Tibetan Plateau in western Sichuan Province, China. Although the species has been known for almost 70 years since its initial discovery in 1950, a thorough osteological description has never been provided. Our study provides a detailed account of the bony anatomy of this species, based on micro computed tomography scanning of multiple specimens collected from the type locality Shuangertang at Bailinshan, Yanyuan County, and several other localities in Sichuan Province. Our revised species diagnosis utilizes both bony and soft anatomical features. Comparative study of the specimens from the type locality in Yanyuan area with those from the nearby Xichang and Mianning areas confirms that they all pertain to Batrachuperus yenyuanensis, thereby removing doubt on the occurrence of the species in the latter areas. With this confirmation, the distribution of the species is extended from the type locality northwards some 180 km to the Mianning area, on both the west and east sides of the Yalong River, which is a major tributary of the upper Yangtze River. This distribution pattern indicates that the biogeographic origin and historical evolution of the species are closely associated with the orogeny of the Hengduan Mountains and formation of the Yalong River. Given the basalmost position of Batrachuperus yenyuanensis in relation to other congeneric species based on molecular studies, however, early expansion of the species distribution by dispersal is expected following the origin of the genus in early-middle Miocene in western Sichuan Province. Thus, the species may well have achieved its current distribution in western Sichuan before the drastic uplift of the Qinghai-Tibetan Plateau in Pliocene.