At the confluence of vicariance and dispersal: Phylogeography of cavernicolous springtails (Collembola: Arrhopalitidae, Tomoceridae) codistributed across a geologically complex karst landscape in Illinois and Missouri

A Study of Oribatid Mites as Potential Intermediate Hosts of Anoplocephalid Tapeworms of Tatra chamois and Tatra marmots from the Tatra Mountains, Central Europe, and Report of a New Intermediate Host for Andrya cuniculi, the Parasite of Leporidae

Tatra chamois (Rupicapra rupicapra tatrica (Blahout 1972)) and Tatra marmot (Marmota marmota latirostris (Kratochvíl 1961)) are significant endemic subspecies of the subalpine and alpine ranges of the Tatra Mountains in Central Europe. In four studied localities in the range of their typical biotopes in Slovakia and Poland, we investigated intestinal parasites of Tatra chamois and Tatra marmots, with an emphasis on anoplocephalid tapeworms. We also studied the occurrence, species diversity, and abundance of oribatid mites as intermediate hosts thereof, and the prevalence of cysticercoid larval stages of anoplocephalid tapeworms in collected oribatids using morphological and molecular methods. Coprological analyses revealed the average positivity of Moniezia spp. in chamois faeces at 23.5% and Ctenotaenia marmotae in marmot samples at 71.1%, with significant differences between the localities under study. Morphological analyses determined the presence of cysticercoids in five oribatid species: Ceratozetes gracilis, Edwardzetes edwardsi, Scheloribates laevigatus, Trichoribates novus, and Tectocepheus velatus sarekensis. This is the first record of T. v. sarekensis as an intermediate host of anoplocephalid tapeworms, as well as the first report of Andrya cuniculi occurrence in the territory of the Tatra Mountains, confirmed also by molecular methods.

New Insights into the Variation and Admixture of the Cave-Dwelling Spider Trogloneta yunnanensis in South China Karst

Subterranean karst caves can contain unexpected biodiversity, but few studies related to spider population genetics have been conducted in the karst area of Southern China. In this study, we investigated the population genetic structure of Trogloneta yunnanensis (Song & Zhu, 1994) based on 73 spider samples from six underground populations in South China Karst. Population genetic structure analysis showed a clear divergence (FST > 0.9 and Nm < 0.05) among populations according to mitochondrial genes. The phylogenetic gene tree constructed by BI and ML methods recovered six geographic clades. Divergence time estimation indicated that the divergence of these six populations can be traced back to the late Pleistocene. We supposed that the geographic isolation led to the extreme population structure. According to this study and previous studies about troglobites living in this region, the subterranean habitats of the Yunnan-Guizhou Plateau may contain many organisms with similar genetic structures. The subterranean biodiversity in the karst area of Southern China needs to be re-evaluated and protected.

Life-history traits drive spatial genetic structuring in Dinaric cave spiders

The subterranean ecosystem exerts strong selection pressures on the organisms that thrive in it. In response, obligate cave-dwellers have developed a series of morphological, physiological, and behavioral adaptations, such as eye reduction, appendage elongation, low metabolic rates or intermittent activity patterns, collectively referred to as troglomorphism. Traditionally, studies on cave organisms have been hampered by the difficulty of sampling (i.e., small population sizes, temporal heterogeneity in specimen occurrence, challenges imposed by the difficult-to-access nature of caves). Here, we circumvent this limitation by implementing a museomics approach. Specifically, we aim at comparing the genetic population structures of five cave spider species demonstrating contrasting life histories and levels of troglomorphism across different caves in the northern Dinarides (Balkans, Europe). We applied a genome-wide hybridization-capture approach (i.e., HyRAD) to capture DNA from 117 historical samples. By comparing the population genetic structures among five species and by studying isolation by distance, we identified deeper population structuring and more pronounced patterns of isolation by distance in the highly troglomorphic Parastalita stygia and Stalita pretneri ground dwellers, while the three web-building Troglohyphantes species, two of which can occasionally be found in surface habitats, showed less structured populations compatible with higher dispersal ability. The spatial distribution of genetic groups revealed common phylogeographic breaks among lineages across the studied species, which hint at the importance of environmental features in driving dispersal potential and shaping underground diversity.

Specific and Intraspecific Diversity of Symphypleona and Neelipleona (Hexapoda: Collembola) in Southern High Appalachia (USA)

Collembola, commonly known as springtails, are important detritivores, abundant in leaf litter and soil globally. Springtails are wingless hexapods with many North American species having wide distributions ranging from as far as Alaska to Mexico. Here, we analyze the occurrence and intraspecific diversity of springtails with a globular body shape (Symphypleona and Neelipleona), in southern high Appalachia, a significant biodiversity hotspot. The peaks of high Appalachia represent ‘sky islands’ due to their physical isolation, and they host numerous endemic species in other taxa. We surveyed globular Collembola through COI metabarcoding, assessing geographic and genetic diversity across localities and species. Intraspecific diversity in globular Collembola was extremely high, suggesting that considerable cryptic speciation has occurred. While we were able to associate morphospecies with described species in most of the major families in the region (Dicyrtomidae, Katiannidae, Sminthuridae, and Sminthurididae), other families (Neelidae, and Arrhopalitidae) are in more pressing need of taxonomic revision before species identities can be confirmed. Due to poor representation in databases, and high intraspecific variability, no identifications were accomplished through comparison with available DNA barcodes.

Hotspot in ferruginous rock may have serious implications in Brazilian conservation policy

A hotspot of subterranean Collembola in ferruginous rock caves and Mesovoid Shallow Substratum is revealed by the analysis of pseudocryptic diversity. The diversity is accessed by detailed description of chaetotaxy and slight variation in morphology of 11 new species of Trogolaphysa Mills, 1938 (Collembola, Paronellidae, Paronellinae) and the 50 previously recorded species of springtails from caves, using optical and electronic microscopy. When combined with recent subterranean surveys, our results show an important reservoir of cave diversity in the Mesovoid Shallow Substratum. Contrastingly the conservation policy for subterranean fauna in metallogenic areas in Brazil prioritizes the caves instead the cave species, which may be extremely detrimental to the fauna in the shallow subterranean habitats not accessible to humans.

Phylogeography and genetic diversity of the widespread katydid Ducetia japonica (Thunberg, 1815) across China

Habitat fragmentation can lower migration rates and genetic connectivity among remaining populations of native species. Ducetia japonica is one of the most widespread katydids in China, but little is known about its genetic structure and phylogeographic distribution. We combined the five-prime region of cytochrome c oxidase subunit I (COI-5P), 11 newly developed microsatellite loci coupled with an ecological niche model (ENM) to examine the genetic diversity and population structure of D. japonica in China and beyond to Laos and Singapore. Both Bayesian inference (BI) and haplotype network methods revealed six mitochondrial COI-5P lineages. The distribution of COI-5P haplotypes may not demonstrate significant phylogeographic structure (N ST > G ST, p > .05). The STRUCTURE analysis based on microsatellite data also revealed six genetic clusters, but discordant with those obtained from COI-5P haplotypes. For both COI-5P and microsatellite data, Mantel tests revealed a significant positive correlation between geographic and genetic distances in mainland China. Bayesian skyline plot (BSP) analyses indicated that the population size of D. japonica's three major mitochondrial COI-5P lineages were seemingly not affected by last glacial maximum (LGM, 0.015-0.025 Mya). The ecological niche models showed that the current distribution of D. japonica was similar to the species' distribution during the LGM period and only slightly extended in northern China. Further phylogeographic studies based on more extensive sampling are needed to identify specific locations of glacial refugia in northern China.

Inferring Evolutionary Timescales without Independent Timing Information: An Assessment of "Universal" Insect Rates to Calibrate a Collembola (Hexapoda) Molecular Clock

Previous estimates of nucleotide substitution rates are routinely applied as secondary or “universal” molecular clock calibrations for estimating evolutionary timescales in groups that lack independent timing information. A major limitation of this approach is that rates can vary considerably among taxonomic groups, but the assumption of rate constancy is rarely evaluated prior to using secondary rate calibrations. Here I evaluate whether an insect mitochondrial DNA clock is appropriate for estimating timescales in Collembola—a group of insect-like arthropods characterized by high levels of cryptic diversity. Relative rates of substitution in cytochrome oxidase subunit 1 (COI) were inferred via Bayesian analysis across a topologically constrained Hexapod phylogeny using a relaxed molecular clock model. Rates for Collembola did not differ significantly from the average rate or from the rates estimated for most other groups (25 of 30), suggesting that (1) their apparent cryptic diversity cannot be explained by accelerated rates of molecular evolution and (2) clocks calibrated using “universal” insect rates may be appropriate for estimating evolutionary timescales in this group. However, of the 31 groups investigated, 10 had rates that deviated significantly from the average (6 higher, 4 lower), underscoring the need for caution and careful consideration when applying secondary insect rate calibrations. Lastly, this study exemplifies a relatively simple approach for evaluating rate constancy within a taxonomic group to determine whether the use of secondary rates are appropriate for molecular clock calibrations.

Genetic diversity of soil invertebrates corroborates timing estimates for past collapses of the West Antarctic Ice Sheet

Changes in the extent of ice sheets through evolutionary timescales have influenced the connectivity of soil invertebrate populations across the Antarctic landscape. We use genetic divergences to estimate isolation times for soil invertebrates along the Transantarctic Mountains. Four species of Collembola (Arthropoda) each showed genetically distinct populations at locations likely isolated for millions of years. Two further species were less genetically diverse although also range restricted. Our genetic data corroborate climate reconstructions and estimates of past warm periods of reduced ice and absent ice shelf in the Ross Sea region, during which time open seaways would have facilitated dispersal of Collembola, and possibly other taxa.

During austral summer field seasons between 1999 and 2018, we sampled at 91 locations throughout southern Victoria Land and along the Transantarctic Mountains for six species of endemic microarthropods (Collembola), covering a latitudinal range from 76.0°S to 87.3°S. We assembled individual mitochondrial cytochrome c oxidase subunit 1 (COI) sequences (n = 866) and found high levels of sequence divergence at both small (<10 km) and large (>600 km) spatial scales for four of the six Collembola species. We applied molecular clock estimates and assessed genetic divergences relative to the timing of past glacial cycles, including collapses of the West Antarctic Ice Sheet (WAIS). We found that genetically distinct lineages within three species have likely been isolated for at least 5.54 My to 3.52 My, while the other three species diverged more recently (<2 My). We suggest that Collembola had greater dispersal opportunities under past warmer climates, via flotation along coastal margins. Similarly increased opportunities for dispersal may occur under contemporary climate warming scenarios, which could influence the genetic structure of extant populations. As Collembola are a living record of past landscape evolution within Antarctica, these findings provide biological evidence to support geological and glaciological estimates of historical WAIS dynamics over the last ca. 5 My.

Geographically structured genetic diversity in the cave beetle Darlingtonea kentuckensis Valentine, 1952 (Coleoptera, Carabidae, Trechini, Trechina)

Cave beetles of the eastern USA are one of many poorly studied groups of insects and nearly all previous work delimiting species is based solely on morphology. This study assesses genetic diversity in the monotypic cave carabid beetle genus DarlingtoneaValentine 1952, to test the relationship between putative geographical barriers to subterranean dispersal and the boundaries of genetically distinct groups. Approximately 400bp of the mitochondrial cytochrome oxidase I (COI) gene was sequenced from up to four individuals from each of 27 populations, sampled from caves along the escarpments of the Mississippian and Cumberland plateaus in eastern Kentucky, USA. The 81 individuals sequenced yielded 28 unique haplotypes. Hierarchical analyses of molecular variance (AMOVA) within and among geographically defined groups tested two a priori hypotheses of structure based on major and minor river drainages, as well as genetic distance clusters defined a posteriori from an unrooted analysis. High genetic differentiation (FST) between populations was found across analyses. The influence of isolation by distance could potentially account for much but not all of the variation found among geographically defined groups at both levels. High variability among the three northernmost genetic clusters (FCT), low variability among populations within clusters (FSC), and low within-cluster Mantel correlations indicate the importance of unidentified likely intra-karst barriers to gene flow separating closely grouped cave populations. Overall phylogeographic patterns are consistent with previous evidence of population isolation among cave systems in the region, revealing geographically structured cryptic diversity in Darlingtonea over its distribution. The landscape features considered a priori in this study were not predictive of the genetic breaks among the three northern clusters, which are genetically distinct despite their close geographic proximity.

Biogeographic barriers, Pleistocene refugia, and climatic gradients in the southeastern Nearctic drive diversification in cornsnakes (Pantherophis guttatus complex)

The southeastern Nearctic is a biodiversity hotspot that is also rich in cryptic species. Numerous hypotheses (e.g., vicariance, local adaptation, and Pleistocene speciation in glacial refugia) have been tested in an attempt to explain diversification and the observed pattern of extant biodiversity. However, previous phylogeographic studies have both supported and refuted these hypotheses. Therefore, while data support one or more of these diversification hypotheses, it is likely that taxa are forming within this region in species-specific ways. Here, we generate a genomic data set for the cornsnakes (Pantherophis guttatus complex), which are widespread across this region, spanning both biogeographic barriers and climatic gradients. We use phylogeographic model selection combined with hindcast ecological niche models to determine regions of habitat stability through time. This combined approach suggests that numerous drivers of population differentiation explain the current diversity of this group of snakes. The Mississippi River caused initial speciation in this species complex, with more recent divergence events linked to adaptations to ecological heterogeneity and allopatric Pleistocene refugia. Lastly, we discuss the taxonomy of this group and suggest there may be additional cryptic species in need of formal recognition.

Cave millipede diversity with the description of six new species from Guangxi, China

We synthesized the current knowledge of cave-dwelling millipede diversity from Guangxi Zhuang Autonomous Region (Guangxi), South China Karst, China and described six new millipede species from four caves from the Guilin area, northeastern Guangxi. Fifty-two cave-dwelling millipedes are known for the region consisting of 38 troglobionts and 14 troglophiles. Of the troglobionts, 24 are presently considered single-cave endemics. New species described here include Hyleoglomerisrukouqusp. nov. and Hyleoglomerisxuxiakeisp. nov. (Family Glomeridae), Hylomusyuanisp. nov. (Family Paradoxosomatidae), Eutrichodesmusjianjiasp. nov. (Family Haplodesmidae), Trichopeltisliangfengdongsp. nov. (Family Cryptodesmidae), and Glyphiulusmaocunsp. nov. (Family Cambalopsidae). Our work also resulted in range expansions of Pacidesmustrifidus Golovatch &amp; Geoffroy, 2014, Blingulussinicus Zhang &amp; Li, 1981 and Glyphiulusmelanoporus Mauriès &amp; Nguyen Duy-Jacquemin, 1997. As with many hypogean animals in Southeast Asia, intensive human activities threaten the persistence of both cave habitats and species. We provide both assessments on the newly described species’ distributions and recommendations for future research and conservation efforts.