MITOCHONDRIAL INTROGRESSION AND INCOMPLETE LINEAGE SORTING THROUGH SPACE AND TIME: PHYLOGENETICS OF CROTAPHYTID LIZARDS

Exploring Mitonuclear Discordance: Ghost Introgression From an Ancient Extinction Lineage in the Odorrana swinhoana Complex

Mitonuclear discordance, the incongruence between mitochondrial DNA (mtDNA) and nuclear DNA (nuDNA), is a well-documented phenomenon with various potential explanations. One emerging hypothesis, ghost introgression, refers to the genetic contribution of an ancient, extinct or unsampled lineage and can now be tested using modern genomic data and demographic models. In this study, we investigated the evolutionary history of the Odorrana swinhoana complex (Anura: Ranidae), which includes O. swinhoana, O. utsunomiyaorum and an unidentified population with highly divergent mtDNA. While mitochondrial phylogeny suggested this population as a basal lineage, nuclear data from ddRADseq revealed it as a mixture of the most derived O. swinhoana nuclear sequences combined with ancient mtDNA. Demographic modelling further supported ghost introgression, as all models incorporating a ghost population outperformed those without it. These findings suggest that an eastward expansion of western O. swinhoana replaced an ancient Odorrana lineage, leaving only its mtDNA and fragments of its nuclear genome in the hybrid population. Our results provide one of the first documented cases of ghost introgression in amphibians and highlight its potential as a widespread evolutionary process. This study also underscores the risks of relying solely on mtDNA for phylogenetic reconstruction and species delimitation.

The Limits of the Metapopulation: Lineage Fragmentation in a Widespread Terrestrial Salamander (Plethodon cinereus)

In vicariant species formation, divergence results primarily from periods of allopatry and restricted gene flow. Widespread species harboring differentiated, geographically distinct sublineages offer a window into what may be a common mode of species formation, whereby a species originates, spreads across the landscape, then fragments into multiple units. However, incipient lineages usually lack reproductive barriers that prevent their fusion upon secondary contact, blurring the boundaries between a single, large metapopulation-level lineage and multiple independent species. Here, we explore this model of species formation in the Eastern Red-backed Salamander (Plethodon cinereus), a widespread terrestrial vertebrate with at least 6 divergent mitochondrial clades throughout its range. Using anchored hybrid enrichment data, we applied phylogenomic and population genomic approaches to investigate patterns of divergence, gene flow, and secondary contact. Genomic data broadly match most mitochondrial groups but reveal mitochondrial introgression and extensive admixture at several contact zones. While species delimitation analyses in Bayesian Phylogenetics and Phylogeography supported 5 lineages of P. cinereus, genealogical divergence indices (gdi) were highly sensitive to the inclusion of admixed samples and the geographic representation of candidate species, with increasing support for multiple species when removing admixed samples or limiting sampling to a single locality per group. An analysis of morphometric data revealed differences in body size and limb proportions among groups, with a reduction of forelimb length among warmer and drier localities consistent with increased fossoriality. We conclude that P. cinereus is a single species, but one with highly structured component lineages of various degrees of independence.

Phylogenomic Discordance is Driven by Wide-Spread Introgression and Incomplete Lineage Sorting During Rapid Species Diversification Within Rattlesnakes (Viperidae: Crotalus and Sistrurus)

-Phylogenomics allows us to uncover the historical signal of evolutionary processes through time and estimate phylogenetic networks accounting for these signals. Insight from genome-wide data further allows us to pinpoint the contributions to phylogenetic signal from hybridization, introgression, and ancestral polymorphism across the genome. Here, we focus on how these processes have contributed to phylogenetic discordance among rattlesnakes (genera Crotalus and Sistrurus), a group for which there are numerous conflicting phylogenetic hypotheses based on a diverse array of molecular datasets and analytical methods. We address the instability of the rattlesnake phylogeny using genomic data generated from transcriptomes sampled from nearly all known species. These genomic data, analyzed with coalescent and network-based approaches, reveal numerous instances of rapid speciation where individual gene trees conflict with the species tree. Moreover, the evolutionary history of rattlesnakes is dominated by incomplete speciation and frequent hybridization, both of which have likely influenced past interpretations of phylogeny. We present a new framework in which the evolutionary relationships of this group can only be understood in light of genome-wide data and network-based analytical methods. Our data suggest that network radiations, like those seen within the rattlesnakes, can only be understood in a phylogenomic context, necessitating similar approaches in our attempts to understand evolutionary history in other rapidly radiating species.

Whole-genome analysis reveals phylogenetic and demographic history of Eurasian perch

The contemporary diversity and distribution of species are shaped by their evolutionary and ecological history. This can be deciphered with the help of phylogenetic and demographic analysis methods, ideally combining and supplementing information from mitochondrial and nuclear genomes. In this study, we investigated the demographic history of Eurasian perch (Perca fluviatilis), a highly adaptable teleost with a distribution range across Eurasia. We combined whole-genome resequencing data with available genomic resources to analyse the phylogeny, phylogeography, and demographic history of P. fluviatilis populations from Europe and Siberia. We identified five highly diverged evolutionary mtDNA lineages, three of which show a strong signal of admixture in the Baltic Sea region. The estimated mean divergence time between these lineages ranged from 0.24 to 1.42 million years. Based on nuclear genomes, two distinct demographic trajectories were observed in European and Siberian samples reflecting contrasting demographic histories ca. 30,000-100,000 years before the present. A comparison of mtDNA and nuclear DNA evolutionary trees and AMOVA revealed concordances, as well as incongruences, between the two types of data, most likely reflecting recent postglacial colonization and hybridization events. Overall, our findings demonstrate the power and usefulness of genome-wide information for delineating historical processes that have shaped the genome of P. fluviatilis. We also highlight the added value of data-mining existing transcriptomic resources to complement novel sequence data, helping to shed light on putative glacial refugia and postglacial recolonization routes.

Cryptic Taxa Revealed through Combined Analysis of Chromosomes and DNA Barcodes: The Polyommatus ripartii Species Complex in Armenia and NW Iran

The detection of cryptic species in complexes that have undergone recent speciation is often difficult, since many standard nuclear markers have not yet accumulated differences between closely related taxa, and differences in mitochondrial markers can be leveled out due to mitochondrial introgressions. In these cases, the use of derived chromosomal characters such as non-ancestral chromosomal numbers and/or unusual karyotype features may be a solution to the species delimitation problem. However, non-ancestral but similar karyotypes may arise secondarily as a result of homoplastic evolution, and their interpretation as homologies may lead to incorrect taxonomic conclusions. In our study, we show that the combined use of mitochondrial DNA barcodes and karyotypes helps to solve this problem and identifies cryptic species in situations where each of these markers does not work individually. Using this approach, we show that the fauna of Armenia and adjacent Iran includes the following cryptic taxa of the Polyommatus ripartii species complex (haploid chromosome number, n in parentheses): P. ripartii paralcestis (n = 90), P. ripartii kalashiani, subsp. nov (n close to 90), P. emmeli, sp. nov. (n = 77-79), P. keleybaricus, sp. nov. (n = 86), P. demavendi belovi (n = 73-75), P. demavendi antonius, subsp. nov. (n = 71-73), P. admetus anatoliensis (n = 79) and P. eriwanensis (n = 29-34). Polyommatus admetus yeranyani is synonymized with P. admetus anatoliensis.

A Palaeogene stem crotaphytid (Aciprion formosum) and the phylogenetic affinities of early fossil pleurodontan iguanians

Pleurodonta is an ancient, diverse clade of iguanian lizard distributed primarily in the Western Hemisphere. Although the clade is a frequent subject of systematic research, phylogenetic resolution among the major pleurodontan clades is elusive. That uncertainty has complicated the interpretations of many fossil pleurodontans. I describe a fossil skull of a pleurodontan lizard from the Palaeogene of Wyoming that was previously allocated to the puzzling taxon Aciprion formosum, and provide an updated morphological matrix for iguanian lizards. Phylogenetic analyses using Bayesian inference demonstrate that the fossil skull is the oldest and first definitive stem member of Crotaphytidae (collared and leopard lizards), establishing the presence of that clade in North America during the Palaeogene. I also discuss new or revised hypotheses for the relationships of several early pleurodontans. In particular, I examine potential evidence for crown-Pleurodonta in the Cretaceous of Mongolia (Polrussia), stem Pleurodonta in the Cretaceous of North America (Magnuviator) and a stem anole in the Eocene of North America (Afairiguana). I suggest that the placement of the fossil crotaphytid is stable to the uncertain phylogeny of Pleurodonta, but recognize the dynamic nature of fossil diagnosis and the potential for updated systematic hypotheses for the other fossils analysed here.

Revision of the African Neoperla Needham, 1905 (Plecoptera: Perlidae: Perlinae) based on morphological and molecular data

Adults of the African species of the genus Neoperla Needham, 1905 (Plecoptera: Perlidae: Perlinae) are revised, and 82 valid species are recognised. Of the original 29 named species, 14 valid ones are redescribed from types, nine lectotypes are designated, and nine new synonymies are proposed. Sixty-two new species are named, and several additional new species only known from material insufficient for formal description are listed. There are one unavailable nomen nudum and three doubtful names. Diagnostic morphological traits of the recognised species are described and illustrated in detail where available of males, females, and eggs. Dichotomous keys to all species are provided. Mitochondrial DNA-data of 71 species-group taxa are available, which permits reliable association of sexes and minimises the risk of synonymies between species known from only one sex. Importantly, this includes DNA sequences from 50 holotypes. The DNA-data will also permit future monitoring with eDNA sequencing and identification of nymphs, which are probably important but so far unidentifiable bioindicators in streams. The genus Neoperla occurs largely in the northern hemisphere and in Africa, where it is restricted to the Ethiopian region south of the Sahara but unknown from Madagascar. A single specimen not clearly distinct from a widespread mainland species was taken on Comoro Island. All African species are endemic, but one African species group has a few outlying members in Asia. The Asian N. montivaga-group is not known from Africa. The following new species are described and named: N. aethiopica n. sp., N. amoena n. sp., N. angolana n. sp., N. bareensis n. sp., N. bella n. sp., N. benti n. sp., N. beta n. sp., N. bipolaris n. sp., N. biserrata n. sp., N. brachyphallus n. sp., N. caeleps n. sp., N. cataractae n. sp., N. claviger n. sp., N. coffea n. sp., N. costata n. sp., N. crenulata n. sp., N. crustata n. sp., N. decorata n. sp., N. dianae n. sp., N. dolium n. sp., N. dundoana n. sp., N. duodeviginti n. sp., N. erinaceus n. sp., N. excavata n. sp., N. filamentosa n. sp., N. funiculata n. sp., N. gibbosa n. sp., N. gordius n. sp., N. heideae n. sp., N. ivanloebli n. sp., N. juxtadidita n. sp., N. kalengonis n. sp., N. larvata n. sp., N. lineata n. sp., N. luhohonis n. sp., N. massevensis n. sp., N. multiserrata n. sp., N. muyukae n. sp., N. nichollsi n. sp., N. occulta n. sp., N. orthonema n. sp., N. pallidogigas n. sp., N. panafricana n. sp., N. pickeri n. sp., N. pilulifera n. sp., N. pirus n. sp., N. planidorsum n. sp., N. plicata n. sp., N. proxima n. sp., N. pusilla n. sp., N. rostrata n. sp., N. sambarua n. sp., N. sassandrae n. sp., N. schuelei n. sp., N. serrula n. sp., N. simplex n. sp., N. socia n. sp., N. sorella n. sp., N. spaghetti n. sp., N. spectabilis n. sp., N. spironema n. sp., N. tangana n. sp., N. tansanica n. sp., N. usambara n. sp., and N. vicina n. sp..

Diversification of the African legless skinks in the subfamily Acontinae (Family Scincidae)

Cladogenic diversification is often explained by referring to climatic oscillations and geomorphic shifts that cause allopatric speciation. In this regard, southern Africa retains a high level of landscape heterogeneity in vegetation, geology, and rainfall patterns. The legless skink subfamily Acontinae occurs broadly across the southern African subcontinent and therefore provides an ideal model group for investigating biogeographic patterns associated with the region. A robust phylogenetic study of the Acontinae with comprehensive coverage and adequate sampling of each taxon has been lacking up until now, resulting in unresolved questions regarding the subfamily's biogeography and evolution. In this study, we used multi-locus genetic markers (three mitochondrial and two nuclear) with comprehensive taxon coverage (all currently recognized Acontinae species) and adequate sampling (multiple specimens for most taxa) of each taxon to infer a phylogeny for the subfamily. The phylogeny retrieved four well-supported clades in Acontias and supported the monophyly of Typhlosaurus. Following the General Lineage Concept (GLC), many long-standing phylogenetic enigmas within Acontias occidentalis and the A. kgalagadi, A. lineatus and A. meleagris species complexes, and within Typhlosaurus were resolved. Our species delimitation analyses suggest the existence of hidden taxa in the A. occidentalis, A. cregoi and A. meleagris species groups, but also suggest that some currently recognized species in the A. lineatus and A. meleagris species groups, and within Typhlosaurus, should be synonymised. We also possibly encountered "ghost introgression" in A. occidentalis. Our inferred species tree revealed a signal of gene flow, which implies possible cross-over in some groups. Fossil evidence calibration dating results showed that the divergence between Typhlosaurus and Acontias was likely influenced by cooling and increasing aridity along the southwest coast in the mid-Oligocene caused by the opening of the Drake Passage. Further cladogenesis observed in Typhlosaurus and Acontias was likely influenced by Miocene cooling, expansion of open habitat, uplifting of the eastern Great Escarpment (GE), and variation in rainfall patterns, together with the effect of the warm Agulhas Current since the early Miocene, the development of the cold Benguela Current since the late Miocene, and their co-effects. The biogeographic pattern of the Acontinae bears close resemblance to that of other herpetofauna (e.g., rain frogs and African vipers) in southern Africa.

Complex cycles of divergence and migration shape lineage structure in the common kingsnake species complex

Aim

The Nearctic is a complex patchwork of habitats and geologic features that form barriers to gene flow resulting in phylogeographic structure and speciation in many lineages. Habitats are rarely stable over geologic time, and the Nearctic has undergone major climatic changes in the past few million years. We use the common kingsnake species complex to study how climate, geography, and history influence lineage formation over a large, complex landscape.

Location

Nearctic/North America.

Taxon

Common kingsnake, Lampropeltis getula, species complex.

Methods

We analyzed genome-wide sequence data from 51 snakes spanning the majority of the species complex's range. We used population clustering, generalized dissimilarity modeling, and coalescent methods to identify the number of genetic clusters within the L. getula complex, infer the environmental correlates of genetic differentiation, and estimate models of divergence and gene flow among lineages.

Results

We identified three major lineages within the L. getula complex and further continuous spatial structure within lineages. The most important ecological correlates of genetic distance in the complex are related to aridity and precipitation, consistent with lineage breaks at the Great Plains/Desert ecotone and the Cochise Filter Barrier. Lineages are estimated to have undergone multiple rounds of isolation and secondary contact, with highly asymmetric migration occurring at present.

Main conclusions

Changing climates combined with a large and geologically complex landscape have resulted in a mosaic of discrete and spatially continuous genetic structure. Multiple rounds of isolation and secondary contact as climate fluctuated over the past ~4.4 My have likely driven the evolution of discrete lineages that maintain high levels of gene flow. Continuous structure is strongly shaped by aridity and precipitation, suggesting roles for major precipitation gradients in helping to maintain lineage identity in the face of gene flow when lineages are in geographic contact.

Testing assertions of widespread introgressive hybridization in a clade of neotropical toads with low mate selectivity (Rhinella granulosa species group)

Discordance between different genomic regions, often identified through multilocus sequencing of selected markers, presents particular difficulties in identifying historical processes which drive species diversity and boundaries. Mechanisms causing discordance, such as incomplete lineage sorting or introgression due to interspecific hybridization, are better identified based on population-level genomic datasets. In the toads of the Rhinella granulosa species group, patterns of mito-nuclear discordance and potential hybridization have been reported by several studies. However, these patterns were proposed based on few loci, such that alternative mechanisms behind gene-tree heterogeneity cannot be ruled out. Using genome-wide ddRADseq loci from a subset of species within this clade, we found only partial concordance between currently recognized species-level taxon boundaries and patterns of genetic structure. While most taxa within the R. granulosa group correspond to clades, genetic clustering analyses sometimes grouped distinct taxonomic units into a single cluster. Moreover, levels of admixture between inferred clusters were limited and restricted to a single taxon pair which is best explained by incomplete lineage sorting as opposed to introgressive hybridization, according to D-statistics results. These findings contradict previous assertions of widespread cryptic diversity and gene flow within the R. granulosa clade. Lastly, our analyses suggest that diversification events within the Rhinella granulosa group mostly dated back to the early Pliocene, being generally younger than species divergences in other closely related clades that present high levels of cross-species gene flow. This finding uniquely contradicts common assertions that this young clade of toads exhibits interspecific hybridization.

A tale about vipers’ tails: phylogeography of black-tailed rattlesnakes

The phylogenetic relationships among black-tailed rattlesnakes remain poorly understood and some authors indicated that the diversity of this group has been underestimated and additional analyses are required to clarify the biogeographic patterns throughout its distribution in Mexico. Therefore, the aim of this study was to elucidate the phylogenetic relationships among black-tailed rattlesnakes across their range, identifying relative divergence times among the main clades and reconstructing the biogeographical history of the group. Three partial mitochondrial genes (ND4, cytb and ATPase6) and one nuclear gene (RAG1) were sequenced to infer the phylogenetic relationships, through the maximum likelihood and Bayesian inference-based methods; demographic history reconstruction was investigated through Bayesian Skyline plot analysis and the ancestral area reconstruction was carried out considering a Bayesian framework. We found strong evidence that the black-tailed rattlesnakes’ group is composed of six clades, which is in agreement with subspecies previously reported. Divergence time estimation indicated that the origin of the C. molossus group could be traced to the middle of the Miocene (~7.71 Mya). Ancestral area reconstruction indicated that early divergence events occurred in Central Mexico, probably related to the geological dynamics of the Trans-Mexican Volcanic Belt. The lineage C. m. oaxacus is the basal member of the C. molossus group. Furthermore, the combination of geological events and changes in Quaternary vegetation may have contributed to the divergence of C. molossus clades. Our results suggest several clades within C. molossus complex could be potentially recognized as separate species.

Ecological speciation by sympatric host shifts in a clade of herbivorous sea slugs, with introgression and localized mitochondrial capture between species

Host shifting in insect-plant systems was historically important to the development of ecological speciation theory, yet surprisingly few studies have examined whether host shifting drives diversification of marine herbivores. When small-bodied consumers feed and also mate on a preferred host, disruptive selection can split a population into host races despite gene flow. Support for host shifts is notably lacking for invertebrates associated with macroalgae, where the scale of dispersal by planktonic larvae often far exceeds the grain of host patchiness, and adults are typically less specialized than terrestrial herbivores. Here, we present a candidate example of ecological speciation in a clade of sea slugs that primarily consume green algae in the genus Caulerpa, including highly invasive species. Ancestral character state reconstructions supported 'sea grapes' (C. racemosa, C. lentillifera) as the ancestral host for a tropical radiation of 12 Elysia spp., with one shift onto alternative Caulerpa spp. in the Indo-Pacific. A Caribbean radiation of three species included symaptric host shifts to Rhipocephalus brevicaulis in the ancestor of E. pratensis Ortea & Espinosa, 1996, and to C. prolifera in E. hamanni Krug, Vendetti & Valdes 2016, plus a niche expansion to a range of Caulerpa spp. in E. subornata Verrill, 1901. All three species are broadly sympatric across the Caribbean but are host-partitioned at a fine grain, and distinct by morphology and at nuclear loci. However, non-recombining mtDNA revealed a history of gene flow between E. pratensis and E. subornata: COI haplotypes from E. subornata were 10.4% divergent from E. pratensis haplotypes from four sites, but closely related to all E. pratensis haplotypes sampled from six Bahamian islands, indicating historical introgression and localized "mitochondrial capture." Disruptive selective likely fueled divergence and adaptation to distinct host environments, indicating ecological speciation may be an under-appreciated driver of diversification for marine herbivores as well as epibionts and other resource specialists.

Diversity and conservation of terrestrial vertebrates (birds, mammals, and reptiles) of Sierra Cucapá, Mexicali, Baja California, Mexico

Knowledge about the biodiversity of Baja California has been obtained mainly from natural protected areas (NPAs), while some unprotected natural areas have been poorly studied. The Sierra Cucapá in the northeast of the peninsula is one example. The objectives of this study are 1) to integrate existing knowledge of bird, mammal, and reptile diversity in Cucapá from public databases, citizen science platforms, and information generated from fieldwork, 2) to identify the spatial distribution of records in the study area, 3) to compare the composition of vertebrate species of Cucapá with that of NPAs of northern part of the peninsula, and 4) to assess the biological conservation value of Cucapá. We obtained records of 150 species of native vertebrates (102 birds, 34 mammals, and 14 reptiles) of which 10 species of birds, four mammals, and seven reptiles are included in a risk extinction category. The different sources of information contributed in a complementary way to the species inventories. Large areas in western and northern Cucapá lack records. The total difference in species composition between Cucapá and nearby NPAs ranged between 58 and 69% for birds, 61 and 79% for mammals, and 69 and 87% for reptiles. The species richness of Cucapá, its particular species composition, the presence of species in risk extinction categories, and the number and size of unexplored areas indicate that this area represents an opportunity for biological conservation in the northern part of the Peninsula. This work provides compelling data for the protection of Cucapá.

An integrative approach reveals high species diversity in the primitively segmented spider genus

Accurate species delimitation is crucial for our understanding of evolution, biodiversity and conservation. However, morphology-based species delimitation alone appears to be prone to taxonomic errors and ineffective for taxa with high interspecific morphological homogeneity or intraspecific morphological variations, as is the case for mesothele and mygalomorph spiders. Combined molecular–morphology species delimitation has shown great potential to delimit species boundaries in such ancient lineages. In the present study, molecular and morphological evidence were integrated to delimit species of the primitively segmented spider genus Songthela Ono, 2000. The cytochrome c oxidase subunit I gene (COI) was sequenced for 192 novel specimens belonging to 12 putative morphospecies. The evolutionary relationships within Songthela and the 12-morphospecies hypothesis were tested in two steps – species discovery and species validation – using four single-locus species delimitation approaches. All species delimitation analyses supported the 12-species hypothesis. Phylogenetic analyses yielded three major clades in Songthela, which are consistent with morphology. Accordingly, we assigned 19 known and 11 new species (S. aokoulong, sp. nov., S. bispina, sp. nov., S. dapo, sp. nov., S. huayanxi, sp. nov., S. lianhe, sp. nov., S. lingshang, sp. nov., S. multidentata, sp. nov., S. tianmen, sp. nov., S. unispina, sp. nov., S. xiujian, sp. nov., S. zizhu, sp. nov.) of Songthela to three species-groups: the bispina-group, the multidentata-group and the unispina-group. Another new species, S. zimugang, sp. nov., is not included in any species groups, but forms a sister lineage to the bispina- and unispina-groups. These results elucidate a high species diversity of Songthela in a small area and demonstrate that integrating morphology with COI-based species delimitation is fast and cost-effective in delimiting species boundaries. http://zoobank.org/urn:lsid:zoobank.org:pub:AF0F5B31-AFAF-4861-9844-445AE8678B67

Biogeographic consequences of shifting climate for the western massasauga (Sistrurus tergeminus)

The western massasauga (Sistrurus tergeminus) is a small pit viper with an extensive geographic range, yet observations of this species are relatively rare. They persist in patchy and isolated populations, threatened by habitat destruction and fragmentation, mortality from vehicle collisions, and deliberate extermination. Changing climates may pose an additional stressor on the survival of isolated populations. Here, we evaluate historic, modern, and future geographic projections of suitable climate for S. tergeminus to outline shifts in their potential geographic distribution and inform current and future management. We used maximum entropy modeling to build multiple models of the potential geographic distribution of S. tergeminus. We evaluated the influence of five key decisions made during the modeling process on the resulting geographic projections of the potential distribution, allowing us to identify areas of model robustness and uncertainty. We evaluated models with the area under the receiver operating curve and true skill statistic. We retained 16 models to project both in the past and future multiple general circulation models. At the last glacial maximum, the potential geographic distribution associated with S. tergeminus occurrences had a stronghold in the southern part of its current range and extended further south into Mexico, but by the mid-Holocene, its modeled potential distribution was similar to its present-day potential distribution. Under future model projections, the potential distribution of S. tergeminus moves north, with the strongest northward trends predicted under a climate scenario increase of 8.5 W/m2. Some southern populations of S. tergeminus have likely already been extirpated and will continue to be threatened by shifting availability of suitable climate, as they are already under threat from desertification of grasslands. Land use and habitat loss at the northern edge of the species range are likely to make it challenging for this species to track suitable climates northward over time.

San Josecito Cave and Its Paleoecological Contributions for Quaternary Studies in Mexico

San Josecito Cave (2250 m elevation) is located nearby Aramberri, Nuevo León, northeastern Mexico, with excavations occurring in 1935–1941 and 1990. It is a paleontological cave and the significance of its faunal data rests in the understanding of the Quaternary ecosystems of the Mexican Plateau and the Southern Plains. This significance is underpinned by a consideration of associated stratigraphic and geochronological data. The fauna is composed of mollusks, amphibians, reptiles, birds, and mammals. More than 30 extinct vertebrate species have been identified, constituting one of the most important Quaternary localities in the Americas. Radiocarbon dates and faunal correlations indicate the excavated deposits represent an interval of time between 45,000 and 11,000 14C years BP. The current synthesis demonstrates that the previous view of the assemblage as a single local fauna is erroneous and that, instead, several successive local faunas are present within a stratigraphic framework. This finding underscores the need for detailed studies of single localities in building paleoenvironmental models. As a corollary, results point to the necessity of including all vertebrate classes represented from a locality in building those models. In addition, the field and analytical methodologies demonstrate the importance of very detailed paleontological excavations, with precise spatial and temporal controls, to assess the taphonomic history of a locality, construct a stratigraphic and geochronological framework, and infer the paleoecological conditions during the time span considered based on the number of local faunas represented. The recognition of San Josecito Cave as an important Late Pleistocene vertebrate paleontological locality is enhanced with the consideration of its faunal data for paleoenvironment reconstruction and possible contribution to Quaternary paleoclimatic modeling.

Population Genetics for Inferring Introduction Sources of the Oriental Fruit Fly, Bactrocera dorsalis: A Test for Quarantine Use in Korea

To infer the introduction sources of the oriental fruit fly, Bactrocera dorsalis, we used a mitochondrial marker to reconstruct the haplotype network and 15 microsatellite loci to reveal genetic structure and relationships between the geographically or temporally different collections from Asia. We performed Approximate Bayesian computations to infer a global origin and a source of the quarantine collections found in Korea. As a result, the 40 populations were divided into three groups, of which genetic similarity is not related to the geographic vicinity. Korean samples had a similar genetic structure to Taiwan and Thailand ones. Our results suggest that the place of origin of the B. dorsalis specimens found in Korea's border quarantine is likely to be Taiwan or Thailand. As the global origin of B. dorsalis, we estimated that Taiwan and Thailand were most likely the global origins of Southeast Asian populations by testing hypothetical scenarios by the approximate Bayesian computation analyses. Our results will allow easier identification of the source region of the forthcoming invasion of quarantined B. dorsalis specimens.

Sympatric occurrence of deeply diverged mitochondrial DNA lineages in Siberian geometrid moths (Lepidoptera: Geometridae): cryptic speciation, mitochondrial introgression, secondary admixture or effect of Wolbachia?

The divergent sympatric mitochondrial lineages within traditionally recognized species present a challenge regularly faced by taxonomists and evolutionary biologists. We encountered this problem when studying the Siberian geometrid moths, Alcis deversata and Thalera chlorosaria. Within each of these species we found two deeply diverged mitochondrial lineages that demonstrated a level of genetic differentiation exceeding the standard interspecific DNA barcode threshold. Using analyses of nuclear genes, morphology, ecological preferences and Wolbachia endosymbionts, we tested five hypotheses that might explain the mitochondrial pattern observed: cryptic speciation, ancestral polymorphism, interspecific mitochondrial introgression, secondary admixture of allopatrically evolved populations and an effect of intracellular Wolbachia endosymbionts. We demonstrate that in A. deversata and Th. chlorosaria the mitochondrial differences are not correlated with differences in nuclear genes, morphology, ecology and Wolbachia infection status, thus not supporting the hypothesis of cryptic species and an effect of Wolbachia. Mitochondrial introgression can lead to a situation in which one species has both its own mitochondrial lineage and the lineage obtained from another species. We found this situation in the species pair Alcis repandata and Alcis extinctaria. We conclude that the mitochondrial heterogeneity in A. deversata and Th. chlorosaria is most likely to be attributable to the secondary admixture of allopatrically evolved populations.

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.

Cytonuclear discordance in the crowned-sparrows, Zonotrichia atricapilla and Zonotrichia leucophrys

The golden-crowned (Zonotrichia atricapilla) and white-crowned (Z. leucophrys) sparrows have been presented as a compelling case for rapid speciation. They display divergence in song and plumage with overlap in their breeding ranges implying reproductive isolation, but have almost identical mitochondrial genomes. Previous research proposed hybridization and subsequent mitochondrial introgression as an alternate explanation, but lacked robust nuclear gene trees to distinguish between introgression and incomplete lineage sorting. We test for signatures of these processes between Z. atricapilla and Z. leucophrys, and investigate the relationships among Z. leucophrys subspecies, using mitochondrial sequencing and a reduced representation nuclear genomic dataset. Contrary to the paraphyly evident in mitochondrial gene trees, we confirmed the reciprocal monophyly of Z. atricapilla and Z. leucophrys using large panels of single nucleotide polymorphisms (SNPs). The pattern of cytonuclear discordance is consistent with limited, historical hybridization and mitochondrial introgression, rather than a recent origin and incomplete lineage sorting between recent sister species. We found evidence of nuclear phylogeographic structure within Z. leucophrys with two distinct clades. Altogether, our results indicate deeper divergences between Z. atricapilla and Z. leucophrys than inferred using mitochondrial markers. Our results demonstrate the limitations of relying solely on mitochondrial DNA for taxonomy, and raise questions about the possibility of selection on the mitochondrial genome during temperature oscillations (e.g. during the Pleistocene). Historical mitochondrial introgression facilitated by past environmental changes could cause erroneous dating of lineage splitting in other taxa when based on mitochondrial DNA alone.

Integrative species delimitation based on COI, ITS, and morphological evidence illustrates a unique evolutionary history of the genus Paracercion (Odonata: Coenagrionidae)

Paracercion are common ‘blue and black’ colored damselflies. We explore the species boundaries of Paracercion (Odonata: Coenagrionidae) using ABGD, bPTP, GMYC and Distance-based clustering. We finally got the molecular data of all nine species of Paracercion. P. hieroglyphicum and P. melanotum were combined into one putative species based on cytochrome c oxidase I (COI). However, they were separated into two putative species based on the nuclear segment including ITS1-5.8S-ITS2 (ITS). This suggests the introgression of mtDNA in Paracercion. Paracercion barbatum and Paracercion melanotum can be separated into two species based on COI, whereas they were combined into one putative species based on ITS, which suggests a hybridization event between them. The lower interspecific divergence (COI: 0.49%) between P. barbatum and Paracercion v-nigrum indicates a recent speciation event in Paracercion. Paracercion sieboldii and P. v-nigrum can be separated into two putative species based on COI, while they were frequently merged into the same putative species based on ITS. This can be explained by incomplete lineage sorting in nDNA. Besides, P. pendulum and P. malayanum were synonymized as junior synonyms of P. melanotum. P. luzonicum was confirmed not to belong to Paracercion. The possibility of introgression, hybridization, recent speciation and incomplete lineage sorting makes species delimitation, based on molecular data, difficult and complicates understanding of the evolutionary history of Paracercion. The discordance in COI and ITS also indicates the value of using markers from different sources in species delimitation studies.

Four Species Linked by Three Hybrid Zones: Two Instances of Repeated Hybridization in One Species Group (Genus Liolaemus)

Hybridization is an evolutionary process that can generate diverse outcomes, such as reinforcing species boundaries, generating new species, or facilitating the introgression of locally-adapted alleles into new genomic backgrounds. Liolaemus is a highly diverse clade of South American lizards with ~260 species and as many as ten new species are described each year. Previous Liolaemus studies have detected gene flow and introgression among species using phylogenetic network methods and/or through comparisons of nuclear and mitochondrial DNA patterns, yet no study has systematically studied hybrid zones between Liolaemus species. Here, we compared three hybrid zones between four species in the Liolaemus fitzingerii group of lizards in Central Argentina where two species, L. melanops and L. xanthoviridis, each hybridize with two other species (L. shehuen and L. fitzingerii). We sampled three transects that were each ~120 km in length and sequenced both mitochondrial and genome-wide SNP data for 267 individuals. In our analyses of nuclear DNA, we also compared bi-allelic SNPs to phased alleles (50 bp RAD loci). Population structure analyses confirmed that boundaries separating species are sharp, and all clines are <65 km wide. Cline center estimates were consistent between SNPs and phased alleles, but cline width estimates were significantly different with the SNPs producing wider estimates. The mitochondrial clines are narrower and shifted 4–20 km southward in comparison to the nuclear clines in all three hybrid zones, indicating that either each of the species has sex-biased dispersal (males northward or females southward), the population densities are unequal, or that the hybrid zones are moving north over time. These comparisons indicate that some patterns of hybridization are similar across hybrid zones (mtDNA clines all narrower and shifted to the south), whereas cline width is variable. Hybridization in the L. fitzingerii group is common and geographically localized; further studies are needed to investigate whether hybrid zones act as hard species boundaries or promoters of speciation through processes such as reinforcement. Nonetheless, this study provides insights into both biotic and abiotic mechanisms helping to maintain species boundaries within the speciose Liolaemus system.

How mitonuclear discordance and geographic variation have confounded species boundaries in a widely studied snake

As DNA sequencing technologies and methods for delimiting species with genomic data become more accessible and numerous, researchers have more tools than ever to investigate questions in systematics and phylogeography. However, easy access to sophisticated computational tools is not without its drawbacks. Choosing the right approach for one's question can be challenging when presented with multitudinous options, some of which fail to distinguish between species and intraspecific population structure. Here, we employ a methodology that emphasizes intensive geographic sampling, particularly at contact zones between populations, with a focus on differentiating intraspecific genetic clusters from species in the Pantherophis guttatus complex, a group of North American ratsnakes. Using a mitochondrial marker as well as ddRADseq data, we find evidence of mitonuclear discordance which has contributed to historical confusion about the relationships within this group. Additionally, we identify geographically and genetically structured populations within the species Pantherophis emoryi that are congruent with previously described morphological variation. Importantly, we find that these structured populations within P. emoryi are highly admixed throughout the range of the species and show no evidence of any reproductive isolation. Our data support a revision of the taxonomy of this group, and we recognize two species within the complex and three subspecies within P. emoryi. This study illustrates the importance of thorough sampling of contact zones and consideration of gene flow when delimiting species in widespread complexes containing parapatric lineages.

A paradigm shift in our view of species drives current trends in biological classification

Discontent about changes in species classifications has grown in recent years. Many of these changes are seen as arbitrary, stemming from unjustified conceptual and methodological grounds, or leading to species that are less distinct than those recognised in the past. We argue that current trends in species classification are the result of a paradigm shift toward which systematics and population genetics have converged and that regards species as the phylogenetic lineages that form the branches of the Tree of Life. Species delimitation now consists of determining which populations belong to which individual phylogenetic lineage. This requires inferences on the process of lineage splitting and divergence, a process to which we have only partial access through incidental evidence and assumptions that are themselves subject to refutation. This approach is not free of problems, as horizontal gene transfer, introgression, hybridisation, incorrect assumptions, sampling and methodological biases can mislead inferences of phylogenetic lineages. Increasing precision is demanded through the identification of both sister relationships and processes blurring or mimicking phylogeny, which has triggered, on the one hand, the development of methods that explicitly address such processes and, on the other hand, an increase in geographical and character data sampling necessary to infer/test such processes. Although our resolving power has increased, our knowledge of sister relationships - what we designate as species resolution - remains poor for many taxa and areas, which biases species limits and perceptions about how divergent species are or ought to be. We attribute to this conceptual shift the demise of trinominal nomenclature we are witnessing with the rise of subspecies to species or their rejection altogether; subspecies are raised to species if they are found to correspond to phylogenetic lineages, while they are rejected as fabricated taxa if they reflect arbitrary partitions of continuous or non-hereditary variation. Conservation strategies, if based on taxa, should emphasise species and reduce the use of subspecies to avoid preserving arbitrary partitions of continuous variation; local variation is best preserved by focusing on biological processes generating ecosystem resilience and diversity rather than by formally naming diagnosable units of any kind. Since many binomials still designate complexes of species rather than individual species, many species have been discovered but not named, geographical sampling is sparse, gene lineages have been mistaken for species, plenty of species limits remain untested, and many groups and areas lack adequate species resolution, we cannot avoid frequent changes to classifications as we address these problems. Changes will not only affect neglected taxa or areas, but also popular ones and regions where taxonomic research remained dormant for decades and old classifications were taken for granted.

Maximum Likelihood Estimation of Species Trees from Gene Trees in the Presence of Ancestral Population Structure

Though large multilocus genomic data sets have led to overall improvements in phylogenetic inference, they have posed the new challenge of addressing conflicting signals across the genome. In particular, ancestral population structure, which has been uncovered in a number of diverse species, can skew gene tree frequencies, thereby hindering the performance of species tree estimators. Here we develop a novel maximum likelihood method, termed TASTI (Taxa with Ancestral structure Species Tree Inference), that can infer phylogenies under such scenarios, and find that it has increasing accuracy with increasing numbers of input gene trees, contrasting with the relatively poor performances of methods not tailored for ancestral structure. Moreover, we propose a supertree approach that allows TASTI to scale computationally with increasing numbers of input taxa. We use genetic simulations to assess TASTI's performance in the three- and four-taxon settings and demonstrate the application of TASTI on a six-species Afrotropical mosquito data set. Finally, we have implemented TASTI in an open-source software package for ease of use by the scientific community.

Rapid colour shift by reproductive character displacement in Cupido butterflies

Reproductive character displacement occurs when competition for successful breeding imposes a divergent selection on the interacting species, causing a divergence of reproductive traits. Here, we show that a disputed butterfly taxon is actually a case of male wing colour shift, apparently produced by reproductive character displacement. Using double digest restriction-site associated DNA sequencing and mitochondrial DNA sequencing we studied four butterfly taxa of the subgenus Cupido (Lepidoptera: Lycaenidae): Cupido minimus and the taxon carswelli, both characterized by brown males and females, plus C. lorquinii and C. osiris, both with blue males and brown females. Unexpectedly, taxa carswelli and C. lorquinii were close to indistinguishable based on our genomic and mitochondrial data, despite displaying strikingly different male coloration. In addition, we report and analysed a brown male within the C. lorquinii range, which demonstrates that the brown morph occurs at very low frequency in C. lorquinii. Such evidence strongly suggests that carswelli is conspecific with C. lorquinii and represents populations with a fixed male brown colour morph. Considering that these brown populations occur in sympatry with or very close to the blue C. osiris, and that the blue C. lorquinii populations never do, we propose that the taxon carswelli could have lost the blue colour due to reproductive character displacement with C. osiris. Since male colour is important for conspecific recognition during courtship, we hypothesize that the observed colour shift may eventually trigger incipient speciation between blue and brown populations. Male colour seems to be an evolutionarily labile character in the Polyommatinae, and the mechanism described here might be at work in the wide diversification of this subfamily of butterflies.

Phylogeography of Schizopygopsis malacanthus Herzenstein (Cypriniformes, Cyprinidae) in relation to the tectonic events and Quaternary climatic oscillations in the Shaluli Mountains Region

The tectonic events and Quaternary climatic oscillations in the Shaluli Mountains (the margin of the southeastern Tibetan Plateau) had an extensive effect on the genetic structure and distribution patterns of this region's terrestrial fauna and flora. It is not yet clear whether similar mechanisms influence this region's fish fauna. Schizopygopsis malacanthus is limited to high-elevation rivers and lakes, and exhibits distinct adaptations to the mountains of, and near to, the Qinghai-Tibet Plateau (QTP). Therefore, this species is a good candidate for investigating patterns of genetic variation resulting from palaeoenvironmental fluctuations in the Shaluli Mountains (China). Here, we used microsatellite and mitochondrial DNA control region sequences to analyze six populations of S. malacanthus collected from the Jinsha and Yalong River drainages. Genealogical analyses identified four maternal lineages and perhaps even four putative species, of which the Ouqu River lineage played a pivotal role during the course of the species' evolution. Two lineages from the Yalong River drainage did not cluster together, whereas those from different drainages grouped together, suggesting tectonic event impacts that possibly altered regional river drainages have been highly influential in population connectivity and gene flow. Population genetic analysis indicated that the geographic barriers and this species preference for higher elevations both played key roles in the divergence of S. malacanthus populations. Demographic tests suggested large-scale spatial synchrony in population fluctuations of S. malacanthus, accompanying dramatic Pleistocene climatic oscillations. It appears that palaeoenvironmental changes in the Shaluli Mountains influenced the distribution and evolution of studied S. malacanthus populations, which provide important information about the factors that influenced the phylogeographic history of this region's fish fauna. Additionally, our study also has implications for ongoing conservation of this vulnerable fish.

Introgression obscures lineage boundaries and phylogeographic history in the western banded gecko, Coleonyx variegatus (Squamata: Eublepharidae)

The geomorphological formation of the Baja California peninsula and the Gulf of California is a principal driver of diversification for the reptiles of North America’s warm deserts. The western banded gecko, Coleonyx variegatus, is distributed throughout the Mojave, Sonoran and Peninsular deserts. In this study we use multilocus sequence data to address deep phylogeographic structure within C. variegatus. Analyses of mtDNA data recover six divergent clades throughout the range of C. variegatus. Topology of the mtDNA gene tree suggests separate origins of peninsular populations with an older lineage in the south and a younger one in the north. In contrast, analyses of multilocus nuclear data provide support for four lineages, corresponding to the subspecies C. v. abbotti, C. v. peninsularis, C. v. sonoriensis and C. v. variegatus. Phylogenetic analyses of the nuclear data recover C. v. abbotti and C. v. peninsularis as a clade, indicating a single origin of the peninsular populations. Discordance between the nuclear and mtDNA data is largely the result of repeated episodes of mtDNA introgression that have obscured both lineage boundaries and biogeographic history. Dating analyses of the combined nuclear and mtDNA data suggest that the peninsular clade diverged from the continental group in the Late Miocene.

Molecular species delimitation in the primitively segmented spider genus Heptathela endemic to Japanese islands

Determining species boundaries forms an important foundation for biological research. However, the results of molecular species delimitation can vary with the data sets and methods that are used. Here we use a two-step approach to delimit species in the genus Heptathela, a group of primitively segmented trapdoor spiders that are endemic to Japanese islands. Morphological evidence suggests the existence of 19 species in the genus. We tested this initial species hypothesis by using six molecular species-delimitation methods to analyse 180 mitochondrial COI sequences of Heptathela sampled from across the known range of the genus. We then conducted a set of more focused analyses by sampling additional genetic markers from the subset of taxa that were inconsistently delimited by the single-locus analyses of mitochondrial DNA. Multilocus species delimitation was performed using two Bayesian approaches based on the multispecies coalescent. Our approach identified 20 putative species among the 180 sampled individuals of Heptathela. We suggest that our two-step approach provides an efficient strategy for delimiting species while minimizing costs and computational time.

Unraveling the Chloroplast Genomes of Two Prosopis Species to Identify Its Genomic Information, Comparative Analyses and Phylogenetic Relationship

Genus Prosopis (family Fabaceae) are shrubby trees, native to arid and semi-arid regions of Asia, Africa, and America and known for nitrogen fixation. Here, we have sequenced the complete chloroplast (cp) genomes of two Prosopis species (P. juliflora and P. cineraria) and compared them with previously sequenced P. glandulosa, Adenanthera microsperma, and Parkia javanica belonging to the same family. The complete genome sequences of Prosopis species and related species ranged from 159,389 bp (A. microsperma) to 163,677 bp (P. cineraria). The overall GC contents of the genomes were almost the similar (35.9–36.6%). The P. juliflora and P. cineraria genomes encoded 132 and 131 genes, respectively, whereas both the species comprised of 85 protein-coding genes higher than other compared species. About 140, 134, and 129 repeats were identified in P. juliflora, P. cineraria and P. glandulosa cp genomes, respectively. Similarly, the maximum number of simple sequence repeats were determined in P. juliflora (88), P. cineraria (84), and P. glandulosa (78). Moreover, complete cp genome comparison determined a high degree of sequence similarity among P. juliflora, P. cineraria, and P. glandulosa, however some divergence in the intergenic spacers of A. microsperma and Parkia javanica were observed. The phylogenetic analysis showed that P. juliflora is closer to P. cineraria than P. glandulosa.

Genetic variation across trophic levels: A test of the correlation between population size and genetic diversity in sympatric desert lizards

Understanding the causes of genetic variation in real populations has been elusive. Competing theories claim that neutral vs. selective processes have a greater influence on the genetic variation within a population. A key difference among theories is the relationship between population size and genetic diversity. Our study tests this empirically by sampling two species of herbivorous lizards (Dipsosaurus dorsalis and Sauromalus ater) and two species of carnivorous lizards (Crotaphytus bicinctores and Gambelia wislizenii) that vary in population size at the same locality, and comparing metrics of genetic diversity. Contrary to neutral expectations, results from four independent loci showed levels of diversity were usually higher for species with smaller population sizes. This suggests that selective processes may be having an important impact on intraspecific diversity in this reptile community, although tests showed little evidence for selection on the loci sequenced for this study. It is also possible that idiosyncratic histories of the focal species may be overriding predictions from simple neutral models. If future studies show that lack of correlation between population size and genetic diversity is common, methods using genetic diversity to estimate population parameters like population size or time to common ancestor should be used with caution, as these estimates are based on neutral theory predictions.

Mitochondrial genomes of the bird genus Piranga: rates of sequence evolution, and discordance between mitochondrial and nuclear markers

We report the characteristics of the mitochondrial genomes of 22 individuals in the bird genus Piranga, including all currently recognized species in the genus (n = 11). Elements follow the standard avian mitogenome series, including two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, 13 protein coding genes, and the mitochondrial control region. Excluding tRNA sequences, sequence divergence rate was lowest in rRNA genes and highest in genes encoding NADH (specifically ND1, ND2, ND3) and the control region. Gene trees assembled from 16 elements (non-tRNAs) varied greatly in topological concordance compared to the recognized species tree (based on thousands of nuclear loci), with no one gene tree precisely recovering all relationships in the genus. We also investigated patterns of concordance between the mitogenome tree and the nuclear species tree and found some discrepancies. Across non-tRNA gene trees (n = 16), the species tree topology was recovered by as few as three elements at a particular node and complete concordance (i.e. 16/16 gene trees matched the species tree topology) was recovered at only one node. We found mitochondrial gene regions that are often used in vertebrate systematics (e.g. CytB, ND2) recovered nearly the exact same topology as the nuclear species tree topology.

Local ancestry analysis reveals genomic convergence in extremophile fishes

The molecular basis of convergent phenotypes is often unknown. However, convergence at a genomic level is predicted when there are large population sizes, gene flow among diverging lineages or strong genetic constraints. We used whole-genome resequencing to investigate genomic convergence in fishes ( Poecilia spp.) that have repeatedly colonized hydrogen sulfide (H₂S)-rich environments in Mexico. We identified genomic similarities in both single nucleotide polymorphisms (SNPs) and structural variants (SVs) among independently derived sulfide spring populations, with approximately 1.2% of the genome being shared among sulfidic ecotypes. We compared these convergent genomic regions to candidate genes for H₂S adaptation identified from transcriptomic analyses and found that a significant proportion of these candidate genes (8%) were also in regions where sulfidic individuals had similar SNPs, while only 1.7% were in regions where sulfidic individuals had similar SVs. Those candidate genes included genes involved in sulfide detoxification, the electron transport chain (the main toxicity target of H₂S) and other processes putatively important for adaptation to sulfidic environments. Regional genomic similarity across independent populations exposed to the same source of selection is consistent with selection on standing variation or introgression of adaptive alleles across divergent lineages. However, combined with previous analyses, our data also support that adaptive changes in mitochondrially encoded subunits arose independently via selection on de novo mutations. Pressing questions remain on what conditions ultimately facilitate the independent rise of adaptive alleles at the same loci in separate populations, and thus, the degree to which evolution is repeatable or predictable. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Evaluating methodologies for species delimitation: the mismatch between phenotypes and genotypes in lichenized fungi (Bryoria sect. Implexae, Parmeliaceae)

In many lichen-forming fungi, molecular phylogenetic analyses lead to the discovery of cryptic species within traditional morphospecies. However, in some cases, molecular sequence data also questions the separation of phenotypically characterised species. Here we apply an integrative taxonomy approach - including morphological, chemical, molecular, and distributional characters - to re-assess species boundaries in a traditionally speciose group of hair lichens, Bryoria sect. Implexae. We sampled multilocus sequence and microsatellite data from 142 specimens from a broad intercontinental distribution. Molecular data included DNA sequences of the standard fungal markers ITS, IGS, GAPDH, two newly tested loci (FRBi15 and FRBi16), and SSR frequencies from 18 microsatellite markers. Datasets were analysed with Bayesian and maximum likelihood phylogenetic reconstruction, phenogram reconstruction, STRUCTURE Bayesian clustering, principal coordinate analysis, haplotype network, and several different species delimitation analyses (ABGD, PTP, GMYC, and DISSECT). Additionally, past population demography and divergence times are estimated. The different approaches to species recognition do not support the monophyly of the 11 currently accepted morphospecies, and rather suggest the reduction of these to four phylogenetic species. Moreover, three of these are relatively recent in origin and cryptic, including phenotypically and chemically variable specimens. Issues regarding the integration of an evolutionary perspective into taxonomic conclusions in species complexes, which have undergone recent diversification, are discussed. The four accepted species, all epitypified by sequenced material, are Bryoria fuscescens, B. glabra, B. kockiana, and B. pseudofuscescens. Ten species rank names are reduced to synonymy. In the absence of molecular data, they can be recorded as the B. fuscescens complex. Intraspecific phenotype plasticity and factors affecting the speciation of different morphospecies in this group of Bryoria are outlined.

Morphometrics and phylogeography of the cave-obligate land snail Helicodiscus barri (Gastropoda, Stylommatophora, Helicodiscidae)

Molecular studies have recently led to the detection of many cryptic species complexes within morphologically ambiguous species formerly undescribed by the scientific community. Organisms such as land snails are at a particularly higher risk of species misidentification and misinterpretation, in that gastropod systematics are based almost entirely on external shell morphology. Subterranean ecosystems are associated with especially high degrees of cryptic speciation, largely owing to the abiotic similarities of these systems. In this study, we attempt to diagnose the potential cryptic diversity in the troglobitic land snail Helicodiscusbarri. Land snails are generally associated with having low vagility, and as such this species’ broad, mosaic distribution indicates the misdiagnosis of this organism as a single species. We analyze both mitochondrial (16S, CO1) and nuclear (28S, H3) genetic data for 23 populations. Phylogeny for H.barri was reconstructed using both maximum-likelihood and Bayesian approaches to assess relationships among populations, and two species delimitation methods (mPTP and ABGD) were used to detect the presence of unique molecular operational taxonomic units (MOTUs). Species delimitation results revealed seven and sixteen MOTUs respectively, suggesting the presence of several cryptic lineages within H.barri. To assess how external shell morphology corresponds with patterns of genetic and environmental variation, two morphometric approaches were used incorporating 115 shells from 31 populations. Both morphometric approaches reveal a significant environmental influence on shell morphology, and one approach showed the significance of MOTU groups. We discuss the delimitation and morphometric results and additionally provide discussion on the taxonomic and conservation implications of this study.

Introgressive hybridization in a Spiny-Tailed Iguana, Ctenosaura pectinata, and its implications for taxonomy and conservation

Introgression, the transmission of genetic material of one taxon into another through hybridization, can have various evolutionary outcomes. Previous studies have detected signs of introgression between western populations of the Mexican endemic and threatened spiny-tailed iguana, Ctenosaura pectinata. However, the extent of this phenomenon along the geographic distribution of the species is unknown. Here, we use multilocus data together with detailed geographic sampling to (1) define genotypic clusters within C. pectinata; (2) evaluate geographic concordance between maternally and biparentally inherited markers; (3) examine levels of introgression between genotypic clusters, and (4) suggest taxonomic modifications in light of this information. Applying clustering methods to genotypes of 341 individuals from 49 localities of C. pectinata and the closely related C. acanthura, we inferred the existence of five genotypic clusters. Contact zones between genotypic clusters with signatures of interbreeding were detected, showing different levels of geographic discordance with mtDNA lineages. In northern localities, mtDNA and microsatellites exhibit concordant distributions, supporting the resurrection of C. brachylopha. Similar concordance is observed along the distribution of C. acanthura, confirming its unique taxonomic identity. Genetic and geographic concordance is also observed for populations within southwestern Mexico, where the recognition of a new species awaits in depth taxonomic revision. In contrast, in western localities a striking pattern of discordance was detected where up to six mtDNA lineages co-occur with only two genotypic clusters. Given that the type specimen originated from this area, we suggest that individuals from western Mexico keep the name C. pectinata. Our results have profound implications for conservation, management, and forensics of Mexican iguanas.

Genomic and mitochondrial evidence of ancient isolations and extreme introgression in the four-lined snake

Comparing mitochondrial and genomic phylogenies is an essential tool for investigating speciation processes, because each genome carries different inheritance properties and evolutionary characteristics. Furthermore, mitonuclear discordance may arise from ecological adaptation, historic isolation, population size changes, and sex-biased dispersal. Closely related taxa are expected to experience gene flow; however, this may not be true for insular populations or populations isolated in refugia. The four-lined snake Elaphe quatuorlineata has a fragmented distribution, separating populations of the Italian and Balkan Peninsulas, whereas several insular Aegean populations of significantly smaller body size (Cyclades island group and Skyros Island, Greece) are currently considered distinct subspecies. We constructed the species-tree phylogeny of this species utilizing genome-wide single nucleotide polymorphisms and a gene-tree based on complete cytochrome b sequences, aiming to detect convergence and discrepancies between biparentally and maternally inherited genomes. Population structuring, phylogenetic patterns and migration events among geographically defined lineages supported our hypothesis of isolation in multiple sub-refugia. Where biogeographical barriers did not restrict subsequent dispersal, extensive genetic exchange occurred between mainland Balkan populations. This process has led to the mitochondrial sweep of an ancestral mitolineage that survived only in peripheral (East Greece) and insular populations (North Cyclades and Skyros). The Central Cyclades represent an ancient lineage for both molecular markers that emerged almost 3.3 Mya. Considering their distinct morphology, insular E. quatuorlineata populations should be the future focus of an extensive sampling, especially since the mitonuclear discordance observed in this species could be related to ecological adaptations, such as the island-dwarfism phenomenon.

Chloroplast genome analysis of box-ironbark Eucalyptus

Eucalyptus L'Hérit. (Myrtaceae) is a taxonomically complex and highly speciose genus that dominates much of Australia's woody vegetation. However, very little information is available about the molecular biology and chloroplast diversity of certain groups, such as Eucalyptus section Adnataria, which is found in many woodland habitats of eastern Australia. We report four new complete chloroplast genomes of Eucalyptus, including three genomes from species previously lacking any chloroplast reference sequences. Plastomes of E. albens, E. conica, E. crebra and E. melliodora assembled using a de novo approach were shown to be largely identical to each other, and similar in size and structure to previously published chloroplast genomes from Eucalyptus. A total of 132 genes (114 single-copy genes and 18 duplicated genes in the IR regions) were identified, and shown to be highly conserved in terms of gene order, content and organization. Slightly higher divergence in the intergenic spacers was identified through comparative genomic analyses. Chloroplast sequences of 35 additional individuals representing 12 species were assembled using a reference guided approach. Rates of nucleotide substitution varied among the protein coding genes, with 17 genes under possible positive selection, and 29 invariant genes. Phylogenetic analysis of either the whole reconstructed plastome sequences or the individual genes revealed extreme discordance with expected species boundaries or higher-level relationships. Plastome relationships were better predicted by geography than by nuclear DNA or taxonomic relationships, suggesting a substantial influence of gene flow over and above the effects of incomplete lineage sorting. These results provide resources for future research and valuable insights into the prevalence of interspecific gene flow among Eucalyptus species.

Population Connectivity and Traces of Mitochondrial Introgression in New Zealand Black-Billed Gulls (Larus bulleri)

Black-billed gulls (Larus bulleri) are endemic to New Zealand and are suspected to be undergoing substantial population declines. They primarily breed on open gravel beds in braided rivers of the South Island—a habitat that is diminishing and becoming increasingly modified. Although management of this species is increasing, little has been published on their movements and demographics. In this study, both mitochondrial DNA (mtDNA) control region domain I and nuclear single nucleotide polymorphisms (SNPs) were examined to help understand the connectivity and population structure of black-billed gulls across the country and to help inform management decisions. Mitochondrial DNA showed no population structure, with high haplotype and low nucleotide diversity, and analyses highlighted mitochondrial introgression with the closely related red-billed gulls (Larus novaehollandiae scopulinus). Nuclear DNA analyses, however, identified two groups, with Rotorua birds in the North Island being distinct from the rest of New Zealand, and isolation-by-distance evident across the South Island populations. Gene flow primarily occurs between nearby colonies with a stepwise movement across the landscape. The importance from a genetic perspective of the more isolated North Island birds (1.6% of total population) needs to be further evaluated. From our results, we infer that the South Island black-billed gull management should focus on maintaining several populations within each region rather than focusing on single specific colonies or river catchments. Future study is needed to investigate the genetic structure of populations at the northern limit of the species’ range, and identify the mechanisms behind, and extent of, the hybridisation between red-billed and black-billed gulls.

The Biogeography of Deep Time Phylogenetic Reticulation

Most phylogenies are typically represented as purely bifurcating. However, as genomic data have become more common in phylogenetic studies, it is not unusual to find reticulation among terminal lineages or among internal nodes (deep time reticulation; DTR). In these situations, gene flow must have happened in the same or adjacent geographic areas for these DTRs to have occurred and therefore biogeographic reconstruction should provide similar area estimates for parental nodes, provided extinction or dispersal has not eroded these patterns. We examine the phylogeny of the widely distributed New World kingsnakes (Lampropeltis), determine if DTR is present in this group, and estimate the ancestral area for reticulation. Importantly, we develop a new method that uses coalescent simulations in a machine learning framework to show conclusively that this phylogeny is best represented as reticulating at deeper time. Using joint probabilities of ancestral area reconstructions on the bifurcating parental lineages from the reticulating node, we show that this reticulation likely occurred in northwestern Mexico/southwestern US, and subsequently, led to the diversification of the Mexican kingsnakes. This region has been previously identified as an area important for understanding speciation and secondary contact with gene flow in snakes and other squamates. This research shows that phylogenetic reticulation is common, even in well-studied groups, and that the geographic scope of ancient hybridization is recoverable.

Genetic analyses reveal cryptic introgression in secretive marsh bird populations

Hybridization is common in bird populations but can be challenging for management, especially if one of the two parent species is of greater conservation concern than the other. King rails (Rallus elegans) and clapper rails (R. crepitans) are two marsh bird species with similar morphologies, behaviors, and overlapping distributions. The two species are found along a salinity gradient with the king rail in freshwater marshes and the clapper in estuarine marshes. However, this separation is not absolute; they are occasionally sympatric, and there are reports of interbreeding. In Virginia, USA, both king and clapper rails are identified by the state as Species of Greater Conservation Need, although clappers are thought to be more abundant and king rails have a higher priority ranking. We used a mitochondrial DNA marker and 13 diagnostic nuclear single nucleotide polymorphisms (SNPs) to identify species, classify the degree of introgression, and explore the evolutionary history of introgression in two putative clapper rail focal populations along a salinity gradient in coastal Virginia. Genetic analyses revealed cryptic introgression with site-specific rates of admixture. We identified a pattern of introgression where clapper rail alleles predominate in brackish marshes. These results suggest clapper rails may be displacing king rails in Virginia coastal waterways, most likely as a result of ecological selection. As introgression can result in various outcomes from outbreeding depression to local adaptation, continued monitoring of these populations would allow further exploration of hybrid fitness and inform conservation management.

The effect of Hengduan Mountains Region (HMR) uplift to environmental changes in the HMR and its eastern adjacent area: Tracing the evolutionary history of Allium section Sikkimensia (Amaryllidaceae)

Exploring the effects of orographic events and climatic shifts on geographic distribution of organism in the Hengduan Mountains Region (HMR) and its eastern adjacent area is crucial to the understanding of the environmental changes to organismal evolution. To gain further insight into these processes, we reconstruct evolutionary history of ten species in Allium section Sikkimensia, distributed across regions abovementioned. Using chloroplast and nuclear sequence variation of 79 populations of these ten Allium species with known morphological preferences, we elucidate the phylogenetic relationship, divergence time, ancestral area and genetic structures. Climatic variables analysis, Isolation by distance (IBD) and environment (IBE) and Species distribution modeling (SDM) were analyzed along different genetic clades. These analyses indicated that the initial split of Sikkimensia was triggered by climate changes following Qinghai-Tibet Plateau sensu lato (QTPsl) uplift during the late Miocene. Subsequently, divergences within lineage (lineage A)/among lineages (lineage C and D) in Sikkimensia may be induced by the intense uplift of the HMR around 3-4 Ma and abrupt intensifying of the Asian monsoon regimes. Furthermore, Sikkimensia populations exhibited lopsided demographic history in the Last Glacial Maximum (LGM), as was indicated by the expansion of their range in the QDM and contraction in the HMR. Our findings appear to suggest that the HMR uplift could have strengthened the orographic difference between the HMR and its eastern adjacent area and led to a colder climate in the HMR, while geological topography also played an important role for taxa to respond the climate change that had taken place in the HMR and its eastern adjacent area during the Pleistocene.