A comparative phylogeographic study reveals discordant evolutionary histories of alpine ground beetles (Coleoptera, Carabidae)

Synergistic effects of Pleistocene geological and climatic events on complex phylogeographic history of widespread sympatric species of Megaloptera in East Asia

Unraveling the phylogeographic histories of species remains a key endeavor for comprehending the evolutionary processes contributing to the rich biodiversity and high endemism found in East Asia. In this study, we explored the phylogeographic patterns and demographic histories of three endemic fishfly and dobsonfly species ( Neochauliodes formosanus, Protohermes costalis, and Neoneuromus orientalis) belonging to the holometabolan order Megaloptera. These species, which share a broad and largely overlapping distribution, were analyzed using comprehensive mitogenomic data. Our findings revealed a consistent influence of vicariance on the population isolation of Neoc. formosanus and P. costalis between Hainan, Taiwan, and the East Asian mainland during the early Pleistocene, potentially hindering subsequent colonization of the later diverged Neon. orientalis to these islands. Additionally, we unveiled the dual function of the major mountain ranges in East Asia, serving both as barriers and conduits, in shaping the population structure of all three species. Notably, we demonstrated that these co-distributed species originated from Southwest, Southern, and eastern Central China, respectively, then subsequently migrated along multi-directional routes, leading to their sympatric distribution on the East Asian mainland. Furthermore, our results highlighted the significance of Pleistocene land bridges along the eastern coast of East Asia in facilitating the dispersal of mountain-dwelling insects with low dispersal ability. Overall, this study provides novel insight into the synergistic impact of Pleistocene geological and climatic events in shaping the diversity and distribution of aquatic insects in East Asia.

Evidence for Admixture and Rapid Evolution During Glacial Climate Change in an Alpine Specialist

The pace of current climate change is expected to be problematic for alpine flora and fauna, as their adaptive capacity may be limited by small population size. Yet, despite substantial genetic drift following post-glacial recolonization of alpine habitats, alpine species are notable for their success surviving in highly heterogeneous environments. Population genomic analyses demonstrating how alpine species have adapted to novel environments with limited genetic diversity remain rare, yet are important in understanding the potential for species to respond to contemporary climate change. In this study, we explored the evolutionary history of alpine ground beetles in the Nebria ingens complex, including the demographic and adaptive changes that followed the last glacier retreat. We first tested alternative models of evolutionary divergence in the species complex. Using millions of genome-wide SNP markers from hundreds of beetles, we found evidence that the N. ingens complex has been formed by past admixture of lineages responding to glacial cycles. Recolonization of alpine sites involved a distributional range shift to higher elevation, which was accompanied by a reduction in suitable habitat and the emergence of complex spatial genetic structure. We tested several possible genetic pathways involved in adaptation to heterogeneous local environments using genome scan and genotype-environment association approaches. From the identified genes, we found enriched functions associated with abiotic stress responses, with strong evidence for adaptation to hypoxia-related pathways. The results demonstrate that despite rapid demographic change, alpine beetles in the N. ingens complex underwent rapid physiological evolution.

Differentiation patterns of emperor moths (Lepidoptera: Saturniidae: Saturniinae) of a continental island: divergent evolutionary history driven by Pleistocene glaciations

Background

On the basis of molecular dating, Pleistocene glaciations have been proposed as the major driving force of biota speciation in the Palearctic and the pre-Quaternary origin of Amazonian taxa. However, the major driving factors in East Asia remain unclear. All 16 saturniine species inhabiting Taiwan with congeners of populations, subspecies, or species in East Asia constitute research objects for addressing the mode of speciation because of the repeated formation and disappearance of a landbridge from the Asian mainland to Taiwan during glacial cycles.

Methods

The genetic divergences of mitochondrial cytochrome c oxidase subunit I (COI) and 16S rDNA and the nuclear 28S rDNA of the saturniine species from Taiwan and the Asian mainland were assessed to determine the monophyly of each genus and species of Saturniinae. Moreover, 519 saturniine COI sequences of 114 taxa from adjacent East and Southeast Asian populations and closely related species were retrieved from GenBank and analyzed. The differentiation timing and possible origination of the insular saturniines were elucidated based on phylogenetic relationships, haplotype networks, and lineage calibrations.

Results

Approximately 90% of intraspecific COI divergence was <2%; all divergences exceeding 2% originated from comparisons between allopatric populations or subspecies. Relationship analyses revealed that multiple introductions likely occurred in insular saturniines and that some East Asian saturniines were paraphyletic as deduced by analyzing endemic insular species. Calibration dating revealed that Taiwanese endemic saturniines split from sibling Asian species 0.2-2.7 million years ago (Mya), whereas subspecific-level and population-level splitting events occurred 0.1-1.7 Mya and 0.2-1.2 Mya, respectively. Moreover, phylogenetic patterns combined with geographical distributions revealed that hill-distributed Taiwanese saturniines are closely related to those from southern China and Southeast Asia, whereas saturniines inhabiting altitudes higher than 1,500 m in Taiwan have siblings distributed in temperate Northeast Asia.

Discussion

The Global DNA Barcoding Initiative was successfully applied to study the population genetic structure in species. Most Formosan saturniines are distinct and monophyletic, reflecting the vicariant barrier of the Taiwan Strait; Pleistocene glacial cycles provided opportunities for insular saturniines to experience repeated isolation from and secondary contact with the continental mainland. Each insular saturniine may have evolved with a unique differentiation timing pattern that possibly emerged in the Early, Middle, or Late Pleistocene with these patterns differing from the consistent pattern that occurred in the temperate Palearctic and tropical Amazonian regions. Moreover, multiple migrations or artificial genetic admixtures may have also occurred, as suggested by the coexistence of two divergent lineages in a few Taiwanese saturniines.

Phylogeny of the supertribe Nebriitae (Coleoptera, Carabidae) based on analyses of DNA sequence data

The phylogeny of the carabid beetle supertribe Nebriitae is inferred from analyses of DNA sequence data from eight gene fragments including one nuclear ribosomal gene (28S), four nuclear-protein coding genes (CAD, topoisomerase 1, PEPCK, and wingless), and three mitochondrial gene fragments (16S + tRNA-Leu + ND1, COI ("barcode" region) and COI ("Pat/Jer" region)). Our taxon sample included 264 exemplars representing 241 species and subspecies (25% of the known nebriite fauna), 39 of 41 currently accepted genera and subgenera (all except Notiokasis and Archileistobrius), and eight outgroup taxa. Separate maximum likelihood (ML) analyses of individual genes, combined ML analyses of nuclear, nuclear protein-coding, and mitochondrial genes, and combined ML and Bayesian analyses of the eight-gene-fragment matrix resulted in a well-resolved phylogeny of the supertribe, with most nodes in the tree strongly supported. Within Nebriitae, 167 internal nodes of the tree (out of the maximum possible 255) are supported by maximum-likelihood bootstrap values of 90% or more. The tribes Notiophilini, Opisthiini, Pelophilini, and Nebriini are well supported as monophyletic but relationships among these are not well resolved. Nippononebria is a distinct genus more closely related to Leistus than Nebria. Archastes, Oreonebria, Spelaeonebria, and Eurynebria, previously treated as distinct genera by some authors, are all nested within a monophyletic genus Nebria. Within Nebria, four major clades are recognized: (1) the Oreonebria Series, including eight subgenera arrayed in two subgeneric complexes (the Eonebria and Oreonebria Complexes); (2) the Nebriola Series, including only subgenus Nebriola; (3) the Nebria Series, including ten subgenera arrayed in two subgeneric complexes, the Boreonebria and Nebria Complexes, with the latter further subdivided into three subgeneric subcomplexes (the Nebria, Epinebriola, and Eunebria Subcomplexes)); and (4) the Catonebria Series, including seven subgenera arrayed in two subgeneric complexes (the Reductonebria and Catonebria Complexes). A strong concordance of biogeography with the inferred phylogeny is noted and some evident vicariance patterns are highlighted. A revised classification, mainly within the Nebriini, is proposed to reflect the inferred phylogeny. Three genus-group taxa (Nippononebria, Vancouveria and Archastes) are given revised status and seven are recognized as new synonymies (Nebriorites Jeannel, 1941 and Marggia Huber, 2014 = Oreonebria Daniel, 1903; Pseudonebriola Ledoux & Roux, 1989 = Boreonebria Jeannel, 1937; Patrobonebria Bänninger, 1923, Paranebria Jeannel, 1937 and Barbonebriola Huber & Schmidt, 2017 = Epinebriola Daniel & Daniel, 1904; and Asionebria Shilenkov, 1982 = Psilonebria Andrewes, 1923). Six new subgenera are proposed and described for newly recognized clades: Parepinebriola Kavanaugh subgen. nov. (type species: Nebria delicata Huber & Schmidt, 2017), Insulanebria Kavanaugh subgen. nov. (type species: Nebria carbonaria Eschscholtz, 1829), Erwinebria Kavanaugh subgen. nov. (type species Nebria sahlbergii Fischer von Waldheim, 1828), Nivalonebria Kavanaugh subgen. nov. (type species: Nebria paradisi Darlington, 1931), Neaptenonebria Kavanaugh subgen. nov. (type species: Nebria ovipennis LeConte, 1878), and Palaptenonebria Kavanaugh subgen. nov. (type species: Nebria mellyi Gebler, 1847). Future efforts to better understand relationships within the supertribe should aim to expand the taxon sampling of DNA sequence data, particularly within subgenera Leistus and Evanoleistus of genus Leistus and the Nebria Complex of genus Nebria.

Ancestral Haplotype Retention and Population Expansion Determine the Complicated Population Genetic Structure of the Hilly Lineage of Neolucanus swinhoei Complex (Coleoptera, Lucanidae) on the Subtropical Taiwan Island

The present study demonstrates that the complicated genetic structure of the hilly lineage of the Neolucanus swinhoei complex was driven by its biological features and habitat requirements as well as hindrance by the CMR during periodical Pleistocene glaciations. The results revealed a tendency of geographical differentiation and major and sub- lineage divergences before and after the Riss glaciation, followed by stable population growth during Würm glaciation. At least four refugia were inferred for N. swinhoei during the Riss-Würm glaciations. The ancestral haplotype retention in the cytochrome oxidase subunit I (COI) gene and compensated substitution in 16S rRNA gene is a possible evolutionary scenario resulting in the inconsistent evolution pattern between COI and 16S rRNA gene coupled with the long-distance dispersal of N. swinhoei. Although the CMR did hinder the dispersal of N. swinhoei, its ancestors may have dispersed to eastern Taiwan through the northern and southern low mountains of the CMR before the Riss glaciation. Our finding suggests that the population growth in the Würm glaciation led a dispersal back to western Taiwan, which is contrast to the more common dispersal scenario from western Taiwan to eastern populations proposed in other studies.

Drainage basins serve as multiple glacial refugia for alpine habitats in the Sierra Nevada Mountains, California

The evolutionary histories of alpine species are often directly associated with responses to glaciation. Deep divergence among populations and complex patterns of genetic variation have been inferred as consequences of persistence within glacier boundaries (i.e., on nunataks), while shallow divergence and limited genetic variation are assumed to result from expansion from large refugia at the edge of ice shields (i.e., massifs de refuge). However, for some species, dependence on specific microhabitats could profoundly influence their spatial and demographic response to glaciation, and such a simple dichotomy may obscure the localization of actual refugia. In this study, we use the Nebria ingens complex (Coleoptera: Carabidae), a water-affiliated ground beetle lineage, to test how drainage basins are linked to their observed population structure. By analysing mitochondrial COI gene sequences and genome-wide single nucleotide polymorphisms, we find that the major drainage systems of the Sierra Nevada Mountains in California best explain the population structure of the N. ingens complex. In addition, we find that an intermediate morphotype within the N. ingens complex is the product of historical hybridization of N. riversi and N. ingens in the San Joaquin basin during glaciation. This study highlights the importance of considering ecological preferences in how species respond to climate fluctuations and provides an explanation for discordances that are often observed in comparative phylogeographical studies.

Rapid speciation and ecological divergence into North American alpine habitats: the Nippononebria (Coleoptera: Carabidae) species complex

The climate-driven species pump hypothesis has been supported in a number of phylogeographic studies of alpine species. Climate-driven shifts in distribution, coupled with rapid demographic change, have led to strong genetic drift and lineage diversification. Although the species pump has been linked to rapid speciation in a number of studies, few studies have demonstrated that ecological divergence accompanies rapid speciation. Here we examine genetic, morphological and physiological variation in members of the ground beetle taxon Nippononebria, to test three competing hypotheses of evolutionary diversification: isolation and incomplete lineage sorting (no speciation), recent speciation without ecological divergence, or recent speciation with ecological divergence into alpine habitats. Genetic data are consistent with recent divergence, with major lineages forming in the last million years. A species tree analysis, in conjunction with morphological divergence in male reproductive traits, support the formation of three recognized Nippononebria taxa. Furthermore, both morphological and physiological traits demonstrate ecological divergence in alpine lineages, with convergent shifts in body shape and thermal tolerance breadth. This provides strong evidence that the climate-driven species pump can generate ecological novelty, though it is argued that spatial scale may be a key determinant of broader patterns of macroevolution in alpine communities.

Comparative phylogeography of trans-Andean freshwater fishes based on genome-wide nuclear and mitochondrial markers

The Neotropical region represents one of the greatest biodiversity hot spots on earth. Despite its unparalleled biodiversity, regional comparative phylogeographic studies are still scarce, with most focusing on model clades (e.g. birds) and typically examining a handful of loci. Here, we apply a genome-wide comparative phylogeographic approach to test hypotheses of codiversification of freshwater fishes in the trans-Andean region. Using target capture methods, we examined exon data for over 1,000 loci combined with complete mitochondrial genomes to study the phylogeographic history of five primary fish species (>150 individuals) collected from eight major river basins in Northwestern South America and Lower Central America. To assess their patterns of genetic structure, we inferred genealogical concordance taking into account all major aspects of phylogeography (within loci, across multiple genes, across species and among biogeographic provinces). Based on phylogeographic concordance factors, we tested four a priori biogeographic hypotheses, finding support for three of them and uncovering a novel, unexpected pattern of codiversification. The four emerging inter-riverine patterns are as follows: (a) Tuira + Atrato, (b) Ranchería + Catatumbo, (c) Magdalena system and (d) Sinú + Atrato. These patterns are interpreted as shared responses to the complex uplifting and orogenic processes that modified or sundered watersheds, allowing codiversification and speciation over geological time. We also find evidence of cryptic speciation in one of the species examined and instances of mitochondrial introgression in others. These results help further our knowledge of the historical geographic factors shaping the outstanding biodiversity of the Neotropics.

Subspecific Differentiation Events of Montane Stag Beetles (Coleoptera, Lucanidae) Endemic to Formosa Island

Taxonomic debates have been carrying on for decades over Formosan stag beetles, which consist of a high proportion of endemic species and subspecies featuring morphological variations associated with local adaptation. With the influence of periodical Pleistocene glaciations and the presence of several mountain ranges, the genetic differentiation and taxonomic recognition, within this medium-size island, of two endemic subspecies for each of four montane stag beetles, i.e. Lucanus ogakii, L. kanoi, Prismognathus davidis, and Neolucanus doro, has been an appealing issue. Based on monophyletic lineages and population structure, possible divergent scenarios have been proposed to clarify the subspecific status for each of the above mentioned stag beetles. Phylogenetic inferences based on COI+16S rDNA+28S rDNA of 240 Formosan lucanids have confirmed most species are monophyletic groups; and the intraspecific (<2%) and interspecific (>2%) genetic distances of the two mitochondrial genes could be applied concordantly for taxonomic identification. On account of Bayesian-based species delimitation, geographic distribution, population structure, and sequence divergences, the subspecific status for L. ogakii, L. kanoi, and Pri. davidis are congruent with their geographic distribution in this island; and the calibration time based on the mitochondrial genes shows the subspecific split events occurred 0.7–1 million years ago. In addition, a more complicated scenario, i.e. genetic differentiation including introgression/hybridization events, might have occurred among L. ogakii, L. kanoi, and L. maculifemoratus. The geological effects of mountain hindrance accompanied by periodical glaciations could have been vital in leading to the geographical subspecific differentiation of these montane stag beetles.