Historical mitochondrial genome introgression confounds species delimitation: evidence from phylogenetic inference in the Odorrana grahami species complex.

Species delimitation is essential to informing conservation policy and understanding ecological and evolutionary processes. Most of our recent gains in knowledge on animal diversity rely on morphological characteristics and mitochondrial (mt) DNA variation. Concordant results based on both have led to an unprecedented acceleration in the identification of new species and enriched the field of taxonomy. However, discordances are also found commonly between morphological and mtDNA evidence. This confounds species delimitation, especially when gene flow or mt genome introgression has occurred. Here, we illustrate how mt genome introgression among species of the Odorrana grahami complex confounds species delimitation using the combined evidence of morphological characters, mt variation, and thousands of nuclear single-nucleotide polymorphisms (SNPs) from genotyping-by-sequencing (GBS). Fifty-eight samples across the distribution of the O. grahami complex were included. The mtDNA matrilineal genealogy indicated 2 clades, with O. grahami and Odorrana junlianensis clustered together. In contrast, all nuclear evidence including gene trees, species trees, and genetic structure analyses based on GBS data support 3 species with distinct genetic clusters. These 3 distinct genetic clusters also correspond to distinct morphological characters. They affirm the distinct taxonomic entities of both O. grahami and O. junlianensis, as well as a third clade distinct from either. Which species the third clade belongs to remains unclear and will require further testing. The nuclear genomic loci contradict the COI evidence, with indications of rampant historical mt genome introgression among the species of the O. grahami complex. These discordant signals previously confused species delimitation efforts in this group. Based on these findings, we recommend the integration of independent data, especially nuclear genomic evidence, in species delimitation so as to be robust against the pitfalls of mt introgression.

Benchmarking ultra-high molecular weight DNA preservation methods for long-read and long-range sequencing.

BACKGROUND: Studies in vertebrate genomics require sampling from a broad range of tissue types, taxa, and localities. Recent advancements in long-read and long-range genome sequencing have made it possible to produce high-quality chromosome-level genome assemblies for almost any organism. However, adequate tissue preservation for the requisite ultra-high molecular weight DNA (uHMW DNA) remains a major challenge. Here we present a comparative study of preservation methods for field and laboratory tissue sampling, across vertebrate classes and different tissue types. RESULTS: We find that storage temperature was the strongest predictor of uHMW fragment lengths. While immediate flash-freezing remains the sample preservation gold standard, samples preserved in 95% EtOH or 20-25% DMSO-EDTA showed little fragment length degradation when stored at 4°C for 6 hours. Samples in 95% EtOH or 20-25% DMSO-EDTA kept at 4°C for 1 week after dissection still yielded adequate amounts of uHMW DNA for most applications. Tissue type was a significant predictor of total DNA yield but not fragment length. Preservation solution had a smaller but significant influence on both fragment length and DNA yield. CONCLUSION: We provide sample preservation guidelines that ensure sufficient DNA integrity and amount required for use with long-read and long-range sequencing technologies across vertebrates. Our best practices generated the uHMW DNA needed for the high-quality reference genomes for phase 1 of the Vertebrate Genomes Project, whose ultimate mission is to generate chromosome-level reference genome assemblies of all ∼70,000 extant vertebrate species.

Invasive Bullfrogs Maintain MHC Polymorphism Including Alleles Associated with Chytrid Fungal Infection.

Maintenance of genetic diversity at adaptive loci may facilitate invasions by non-native species by allowing populations to adapt to novel environments, despite the loss of diversity at neutral loci that typically occurs during founder events. To evaluate this prediction, we compared genetic diversity at major histocompatibility complex (MHC) and cytochrome b (cytb) loci from 20 populations of the American bullfrog (Rana catesbeiana) across theinvasive and native ranges in North America and quantified the presence of the pathogen Batrachochytrium dendrobatidis (Bd). Compared to native populations, invasive populations had significantly higher Bd prevalence and intensity, significantly higher pairwise MHC and cytb FST, and significantly lower cytb diversity, but maintained similar levels of MHC diversity. The two most common MHC alleles (LiCA_B and Rapi_33) were associated with a significant decreased risk of Bd infection, and we detected positive selection acting on four peptide binding residues. Phylogenetic analysis suggested invasive populations likely arose from a single founding population in the American Midwest with a possible subsequent invasion in the northwest. Overall, our study suggests that the maintenance of diversity at adaptive loci may contribute to invasion success and highlights the importance of quantifying diversity at functional loci to assess the evolutionary potential of invasive populations.

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.

From Gondwana to the Yellow Sea, evolutionary diversifications of true toads Bufo sp. in the Eastern Palearctic and a revisit of species boundaries for Asian lineages.

Taxa with vast distribution ranges often display unresolved phylogeographic structures and unclear taxonomic boundaries resulting in hidden diversity. This hypothesis-driven study reveals the evolutionary history of Bufonidae, covering the phylogeographic patterns found in Holarctic bufonids from the West Gondwana to the phylogenetic taxonomy of Asiatic true toads in the Eastern Palearctic. We used an integrative approach relying on fossilized birth-death calibrations, population dynamics, gene-flow, species distribution, and species delimitation modeling to resolve the biogeography of the clade and highlight cryptic lineages. We verified the near-simultaneous Miocene radiations within Western and Eastern Palearctic Bufo, c. 14.49-10.00 Mya, temporally matching with the maximum dust outflows in Central Asian deserts. Contrary to earlier studies, we demonstrated that the combined impacts of long dispersal and ice-age refugia equally contributed to the current genetic structure of Bufo in East Asia. Our findings reveal a climate-driven adaptation in septentrional Eastern Asian Bufo, explaining its range shifts toward northern latitudes. We resolve species boundaries within the Eastern Palearctic Bufo, and redefine the taxonomic and conservation units of the northeastern species: B. sachalinensis and its subspecies.

The east Asian Asiatic toad (also known by the latin name Bufo gargarizans) lives in a wide range of habitats across East Asia including forests, meadows and cultivated land. However, it remains unclear how these toads evolved and became so widespread ­ partly because it has proved difficult for researchers to clearly define the species and what distinguishes it from other closely-related species of toads (collectively known as Bufo toads). Othman et al. combined several bioinformatics techniques to study Asiatic toads and 38 other species of bufonid toads from across the globe. This approach found that Bufo toads first emerged in eastern Asia between 14 and 10 million years ago. This coincides with a point in time when large swathes of land in central Asia turned from adequate to sustain toad populations into desert, suggesting the change in climate prompted the toads to migrate eastwards from central Asia. The Bufo toads then divided into two groups of species: one in mainland East Asia and the other in Japan. Furthermore, the study revealed there is more genetic diversity ­ that is, more variety in the DNA of individuals ­ in Asiatic toads than previously thought. The findings also help to define several other species of Bufo toads more clearly and describe a new toad species restricted to the Korean Peninsula, northeastern China and the Russian Primorye region: the Sakhalin toad (Bufo sachalinensis). This work demonstrates that a large-scale study of many species across the globe can be used to understand how the species evolved and more clearly distinguish one species from another. The findings of Othman et al. will be of interest to both professional and citizen scientists interested in the natural history of Asia. Furthermore, as several species of Bufo toads are in decline in the wild, they provide evidence that may aid future efforts to conserve them.

A combined approach of mitochondrial DNA and anchored nuclear phylogenomics sheds light on unrecognized diversity, phylogeny, and historical biogeography of the torrent frogs, genus Amolops (Anura: Ranidae).

The genus Amolops ("torrent frogs") is one of the most species-rich genera in Ranidae, with 59 recognized species. This genus currently includes six species groups diagnosed mainly by morphology. Several recent molecular studies indicated that the classification of species groups within Amolops remains controversial, and key nodes in the phylogeny have been inadequately resolved. In addition, the diversity of Amolops remains poorly understood, especially for those from incompletely sampled regions. Herein, we investigate species-level diversity within the genus Amolops throughout southern China and Southeast Asia, and infer evolutionary relationships among the species using mtDNA data (16S, COI, and ND2). Molecular analyses indicate nine unnamed species, mostly distributed in the Himalayas. We then utilized anchored hybrid enrichment to generate a dataset representing the major mitochondrial lineages to resolve phylogenetic relationships, biogeography, and pattern of species diversification. Our resulting phylogeny strongly supports the monophyly of four previously identified species groups (the A. ricketti, A. daiyunensis, A. hainanensis, and A. monticola groups), but paraphyly for the A. mantzorum and A. marmoratus groups, as previously defined. We erect one new species group, the A. viridimaculatus group, and recognize Dubois' (1992) subgenus Amo as the A. larutensis species group. Biogeographic analysis suggests that Amolops originated on the Indo-Burma/Thai-Malay Peninsula at the Eocene/Oligocene boundary, and dispersed outward, exemplifying a common pattern observed for the origin of Asian biodiversity. The early divergence within Amolops coincides with the Himalayan uplift and the lateral extrusion of Indochina at the Oligocene/Miocene boundary. Our results show that paleoclimatic and geomorphological events have profoundly influenced the patterns of lineage diversification within Amolops.

A new species of the genus Megophrys Gunther, 1864 (Amphibia: Anura: Megophryidae) from Mount Wuyi, China.

A new species of horned toad, Megophrys ombrophila sp. nov., is described based on specimens found from Guadun village from Mount Wuyi in northwestern Fujian Province, China. The species is distinguished from other described Megophrys by morphology, bioacoustics, and molecular data of the 12S and 16S mitochondrial loci. The species is characterized by its rotund stature, head length approximately equal to head width, canthus rostralis well developed, tympanum large and distinct, vomerine teeth absent, margin of tongue smooth, not notched from behind, heels of the feet not meeting when femurs are held at 90° to the axis of the body and tibias are depressed against the femur, toes weakly webbed at base, dorsal skin mostly smooth with scattered granules and ridges, usually the presence of two discontinuous dorsolateral ridges, supratympanic fold distinct and well-developed, and with females ranging from 32.8-35 mm snout-to-vent length, and males ranging from 27.4-34.5 mm SVL. In its type locality, the species is sympatric with M. boettgeri and M. kuatunensis; and phylogenetically, the most closely related described species is M. obesa from southwest Guangdong province. As more cryptic species are described within the subfamily Megophryinae, the more apparent becomes the need for fine-scale molecular and phenotypic assessment to capture the many forms of this hyper-diverse group.

Hidden diversity within the depauperate genera of the snake tribe Lampropeltini (Serpentes, Colubridae).

Accurate representation of lineage diversity through complete taxon sampling is crucial to understanding the evolution of biodiversity, particularly when using molecular phylogenetics to estimate evolutionary relationships. In this interest, taxonomic diversity is often used as a proxy for lineage diversity even though the two concepts are not synonymous. We explore this within the snake tribe Lampropeltini which includes some of the most conspicuous and heavily studied snakes in North America. Both the taxonomy and hypothesized relationships within this tribe have been in flux. The number of species has increased from 23 to 51 over the last thirty years, predominately within three of the nine genera (Lampropeltis, Pantherophis, Pituophis). The remaining six depauperate genera (Arizona, Bogertophis, Cemophora, Pseudelaphe, Rhinocheilus, and Senticolis) have been poorly represented in phylogenetic studies. To estimate evolutionary relationships and determine if the dichotomy in depauperate and speciose genera within Lampropeltini is a function of taxon sampling or truly represents the lineage diversity, we estimated the phylogeny of this group using nuclear and mitochondrial loci in a concatenated and coalescent framework with the largest sampling of the six depauperate genera to date. In addition, we estimated the divergence dates among the genera to assess whether the instability of Lampropeltini phylogenetic relationships is due to an adaptive radiation. While some nodes still remain unresolved, the generic-level relationships we recovered agree with those of a recent next-generation study that used a much larger set of loci for fewer individuals. We also tested two putative species, Arizona pacata and Pseudelaphe phaescens, for the first time phylogenetically and find evidence that they are distinct lineages. Overall, we find that the taxonomic and genetic diversity are not correlated in Lampropeltini and that representing putative diversity in phylogenies will lead to a better estimate of evolutionary histories, especially in groups with complex radiations.