Fine-scale comparative phylogeography of a sympatric sister species triplet of subterranean diving beetles from a single calcrete aquifer in Western Australia

Micro-endemic species of snails and amphipods show population genetic structure across very small geographic ranges

Understanding variation in population genetic structure, even across small distances and for species with extremely limited ranges, is critical for conservation planning and the development of effective management strategies for imperiled species. Organisms that occupy the same geographic extent can maintain different population structures, ranging from highly diverged to panmictic. Such differences can result from differences in biological characteristics such as dispersal ability or demographic history. We used microsatellite loci to evaluate population genetic structure and variation of four desert spring invertebrates having high to low dispersal ability: the lung snail Physa acuta, two species of gilled snails (Juturnia kosteri and Pyrgulopsis roswellensis; family Hydrobiidae) and the amphipod Gammarus desperatus. The study location represents entire species ranges for the micro-endemic hydrobiids and G. desperatus, while P. acuta is ubiquitous throughout much of North America. We found little evidence of significant population genetic structure for P. acuta and J. kosteri, but much more for P. roswellensis and G. desperatus. Our results demonstrate differences in habitat preference and/or dispersal ability between the species. This information provides insight into how gene flow shapes varying population genetic structure between species across small spatial scales (<100 km2). Most importantly, our results suggest that conservation agencies should not consider these micro-endemic species to be composed of single populations, but rather, that management plans for such species should account for population genetic variation across the species' ranges.

Connectivity, not short-range endemism, characterises the groundwater biota of a northern Australian karst system

Groundwater ecosystems have a diverse and unique fauna, often dominated by Crustacea and generally characterised by short range endemics confined to single aquifers. Much of this knowledge has come from studies conducted either in fractured rock aquifers or alluvial aquifers. Karstic subterranean environments are present in the Cambrian Limestone Aquifer (CLA) in the Northern Territory, Australia, a freshwater aquifer which spans an area of ~28,000 km². The presence of underground caverns and channels potentially allows extensive connectivity within this groundwater system. The emerging shale gas industry in the Beetaloo region, which underlies the CLA, provided the impetus to undertake the first survey of the potential existence of a stygofaunal community. Twenty-six groundwater wells (bores) and two springs were sampled in August and October 2019, across a distance of ~500 km, from the sub-tropical Mataranka region in the north to the semi-arid Barkly Tablelands in the south. Plankton nets and motorised pumps were used to collect water samples and conventional microscope-based morphological examinations in conjunction with environmental DNA (eDNA) were used to determine the presence of stygofauna. COI barcoding and 16S rRNA regions were also used for phylogenetic analysis. All stygofaunal communities were dominated by crustaceans, namely shrimps, amphipods, ostracods, copepods and syncarids. This fauna showed little affinity with the stygofauna recorded from more extensively sampled aquifers in north-western Australia, with new genera and species present in the CLA. eDNA analysis showed the presence of diverse biota at sites where direct water sampling for intact animals was difficult. COI and 16S analysis confirmed that a species of blind shrimp, Parisia unguis, occurred extensively throughout the aquifer, over a distance of at least ~300 km. The presence of Pa. unguis at widely separated sites across the CLA is consistent with substantial connectivity within the aquifer. This connectivity indicates that the risk of groundwater contamination from fracking chemicals needs to be adequately mitigated to prevent widespread effects.

The critical thermal maximum of diving beetles (Coleoptera: Dytiscidae): a comparison of subterranean and surface-dwelling species

Thermal tolerance limits in animals are often thought to be related to temperature and thermal variation in their environment. Recently, there has been a focus on studying upper thermal limits due to the likelihood for climate change to expose more animals to higher temperatures and potentially extinction. Organisms living in underground environments experience reduced temperatures and thermal variation in comparison to species living in surface habitats, but how these impact their thermal tolerance limits are unclear. In this study, we compare the thermal critical maximum (CTmax) of two subterranean diving beetles (Dytiscidae) to that of three related surface-dwelling species. Our results show that subterranean species have a lower CTmax (38.3-39.0°C) than surface species (42.0-44.5°C). The CTmax of subterranean species is ∼10°C higher than the highest temperature recorded within the aquifer. Groundwater temperature varied between 18.4°C and 28.8°C, and changes with time, depth and distance across the aquifer. Seasonal temperature fluctuations were 0.5°C at a single point, with the maximum heating rate being ∼1000x lower (0.008°C/hour) than that recorded in surface habitats (7.98°C/hour). For surface species, CTmax was 7-10°C higher than the maximum temperature in their habitats, with daily fluctuations from ∼1°C to 16°C and extremes of 6.9°C and 34.9°C. These findings suggest that subterranean dytiscid beetles are unlikely to reach their CTmax with a predicted warming of 1.3-5.1°C in the region by 2090. However, the impacts of long-term elevated temperatures on fitness, different life stages and other species in the beetle's trophic food web are unknown.

Evidence for speciation underground in diving beetles (Dytiscidae) from a subterranean archipelago

Most subterranean animals are assumed to have evolved from surface ancestors following colonization of a cave system; however, very few studies have raised the possibility of "subterranean speciation" in underground habitats (i.e., obligate cave-dwelling organisms [troglobionts] descended from troglobiotic ancestors). Numerous endemic subterranean diving beetle species from spatially discrete calcrete aquifers in Western Australia (stygobionts) have evolved independently from surface ancestors; however, several cases of sympatric sister species raise the possibility of subterranean speciation. We tested this hypothesis using vision (phototransduction) genes that are evolving under neutral processes in subterranean species and purifying selection in surface species. Using sequence data from 32 subterranean and five surface species in the genus Paroster (Dytiscidae), we identified deleterious mutations in long wavelength opsin (lwop), arrestin 1 (arr1), and arrestin 2 (arr2) shared by a sympatric sister-species triplet, arr1 shared by a sympatric sister-species pair, and lwop and arr2 shared among closely related species in adjacent calcrete aquifers. In all cases, a common ancestor possessed the function-altering mutations, implying they were already adapted to aphotic environments. Our study represents one of the first confirmed cases of subterranean speciation in cave insects. The assessment of genes undergoing pseudogenization provides a novel way of testing modes of speciation and the history of diversification in blind cave animals.

Extreme genetic diversity among springtails (Collembola) in subterranean calcretes of arid Australia

The subterranean islands hypothesis for calcretes of the Yilgarn region in Western Australia applies to many stygobitic (subterranean-aquatic) species that are "trapped" evolutionarily within isolated aquifers due to their aquatic lifestyles. In contrast, little is known about the distribution of terrestrial-subterranean invertebrates associated with the calcretes. We used subterranean Collembola from the Yilgarn calcretes to test the hypothesis that troglobitic species, those inhabiting the subterranean unsaturated (non-aquatic) zone of calcretes, are also restricted in their distribution and represent reciprocally monophyletic and endemic lineages. We used the barcoding fragment of the mtDNA cytochrome c oxidase subunit 1 (COI) gene from 183 individuals to reconstruct the phylogenetic history of the genus Pseudosinella Schäffer (Collembola, Lepidocyrtidae) from 10 calcretes in the Yilgarn. These calcretes represent less than 5% of the total possible calcretes in this region, yet we show that their diversity for subterranean Collembola comprises a minimum of 25 new species. Regionally, multiple levels of diversity exist in Pseudosinella, indicative of a complex evolutionary history for this genus in the Yilgarn. These species have probably been impacted by climatic oscillations, facilitating their dispersal across the landscape. The results represent a small proportion of the undiscovered diversity in Australia's arid zone.

Cutaneous respiration by diving beetles from underground aquifers of Western Australia (Coleoptera: Dytiscidae)

Insects have a gas-filled respiratory system, which provides a challenge for those that have become aquatic secondarily. Diving beetles (Dytiscidae) use bubbles on the surface of their bodies to supply O2 for their dives and passively gain O2 from the water. However, these bubbles usually require replenishment at the water's surface. A highly diverse assemblage of subterranean dytiscids has evolved in isolated calcrete aquifers of Western Australia with limited/no access to an air–water interface, raising the question of how they are able to respire. We explored the hypothesis that they use cutaneous respiration by studying the mode of respiration in three subterranean dytiscid species from two isolated aquifers. The three beetle species consume O2 directly from the water, but they lack structures on their bodies that could have respiratory function. They also have a lower metabolic rate than other insects. O2 boundary layers surrounding the beetles are present, indicating that O2 diffuses into the surface of their bodies via cutaneous respiration. Cuticle thickness measurements and other experimental results were incorporated into a mathematical model to understand whether cutaneous respiration limits beetle size. The model indicates that the cuticle contributes considerably to resistance in the O2 cascade. As the beetles become larger, their metabolic scope narrows, potentially limiting their ability to allocate energy to mating, foraging and development at sizes above approximately 5 mg. However, the ability of these beetles to utilise cutaneous respiration has enabled the evolution of the largest assemblage of subterranean dytiscids in the world.

Blind scolopendrid centipedes of the genus Cormocephalus from subterranean habitats in Western Australia (Myriapoda: Scolopendromorpha: Scolopendridae)

Only a single blind species is known in the centipede family Scolopendridae, representing the monotypic genus Tonkinodentus Schileyko, 1992, from Vietnam. All of more than 400 other species have four ocelli on each side of the cephalic plate. A complex of three new blind species of the genus Cormocephalus Newport, 1844, is described from the subterranean fauna of the central Pilbara region of Western Australia. Phylogenies based on sequence data for the barcode region of COI and a concatenated matrix that also includes 12S rRNA, 28S rRNA and ITS2 unite the blind Pilbara species as a monophyletic group, albeit with moderate bootstrap support, informally named the C. sagmus species group. Cormocephalus sagmus, C. pyropygus and C. delta spp. nov. supplement 17 epigean congeners previously described from Australia. The new species are all morphologically similar, but can be distinguished using the shape and spinulation of the ultimate leg prefemur. Two additional genetically distinct lineages were recovered that are not described, owing to the specimens being immature or lacking diagnostic morphological characters. The subterranean radiation in the Pilbara is more closely related to species from forests in the south-west of Western Australia than to congeners from the arid zone. http://zoobank.org/urn:lsid:zoobank.org:pub:6F67FD31-A373-4DC5-A5FD-374D32DEE02C

Understanding subterranean variability: the first genus of Bathynellidae (Bathynellacea, Crustacea) from Western Australia described through a morphological and multigene approach

The number of subterranean taxa discovered in the north of Western Australia has substantially increased due to the requirements for environmental surveys related to mining development. Challenges in estimating subterranean biodiversity and distributions are related to lack of knowledge of taxa with convergent morphological characters in a largely unobservable ecosystem setting. An integrated approach is warranted to understand such complexity. Bathynellidae occur in most Australian aquifers, but only one species has been described so far, and the group lacks a reliable taxonomic framework. A new genus and one new species from the Pilbara region of Western Australia, Pilbaranella ethelensis, gen. et sp. nov., is described using both morphological and molecular data. Three additional species of Pilbaranella are defined through mitochondrial and nuclear genes, using Automatic Barcode Gap Discovery and Poisson Tree Processes species delimitation methods. A comparison of morphology and 18S rRNA sequences between Pilbaranella, gen. nov. and known lineages provides the evidentiary basis for the decision to establish a new genus. This study provides a morphological and molecular framework to work with Bathynellidae, especially in Australia where a highly diverse fauna remains still undescribed.

The first complete mitochondrial genomes of subterranean dytiscid diving beetles (Limbodessus and Paroster) from calcrete aquifers of Western Australia

Comparative analyses of mitochondrial (mt) genomes may provide insights into the genetic changes, associated with metabolism, that occur when surface species adapt to living in underground habitats. Such analyses require comparisons among multiple independently evolved subterranean species, with the dytiscid beetle fauna from the calcrete archipelago of central Western Australia providing an outstanding model system to do this. Here, we present the first whole mt genomes from four subterranean dytiscid beetle species of the genera Limbodessus (L. palmulaoides) and Paroster (P. macrosturtensis, P. mesosturtensis and P. microsturtensis) and compare genome sequences with those from surface dytiscid species. The mt genomes were sequenced using a next-generation sequencing approach employing the Illumina Miseq system and assembled de novo. All four mt genomes are circular, ranging in size from 16 504 to 16 868 bp, and encode 37 genes and a control region. The overall structure (gene number, orientation and order) of the mt genomes is the same as that found in eight sequenced surface species, but with genome size variation resulting from length variation of intergenic regions and the control region . Our results provide a basis for future investigations of adaptive evolutionary changes that may occur in mt genes when species move underground.

Comparative phylogeography in the Atlantic forest and Brazilian savannas: pleistocene fluctuations and dispersal shape spatial patterns in two bumblebees

Background

Bombus morio and B. pauloensis are sympatric widespread bumblebee species that occupy two major Brazilian biomes, the Atlantic forest and the savannas of the Cerrado. Differences in dispersion capacity, which is greater in B. morio, likely influence their phylogeographic patterns. This study asks which processes best explain the patterns of genetic variation observed in B. morio and B. pauloensis, shedding light on the phenomena that shaped the range of local populations and the spatial distribution of intra-specific lineages.

Results

Results suggest that Pleistocene climatic oscillations directly influenced the population structure of both species. Correlative species distribution models predict that the warmer conditions of the Last Interglacial contributed to population contraction, while demographic expansion happened during the Last Glacial Maximum. These results are consistent with physiological data suggesting that bumblebees are well adapted to colder conditions. Intra-specific mitochondrial genealogies are not congruent between the two species, which may be explained by their documented differences in dispersal ability.

Conclusions

While populations of the high-dispersal B. morio are morphologically and genetically homogeneous across the species range, B. pauloensis encompasses multiple (three) mitochondrial lineages, and show clear genetic, geographic, and morphological differences. Because the lineages of B. pauloensis are currently exposed to distinct climatic conditions (and elevations), parapatric diversification may occur within this taxon. The eastern portion of the state of São Paulo, the most urbanized area in Brazil, represents the center of genetic diversity for B. pauloensis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0803-0) contains supplementary material, which is available to authorized users.

Stygofauna enhance prokaryotic transport in groundwater ecosystems

More than 97% of the world’s freshwater reserves are found in aquifers, making groundwater one of the most important resources on the planet. Prokaryotic communities in groundwater underpin the turnover of energy and matter while also maintaining groundwater purity. Thus, knowledge of microbial transport in the subsurface is crucial for maintaining groundwater health. Here, we describe for the first time the importance of stygofauna as vectors for prokaryotes. The “hitch-hiking” prokaryotes associated with stygofauna may be up to 5 orders of magnitude higher in abundance and transported up to 34× faster than bulk groundwater flow. We also demonstrate that prokaryotic diversity associated with stygofauna may be higher than that of the surrounding groundwater. Stygofauna are a newly recognized prokaryotic niche in groundwater ecosystems that have the potential to transport remediating, water purifying and pathogenic prokaryotes. Therefore, stygofauna may influence ecosystem dynamics and health at a microbial level, and at a larger scale could be a new source of prokaryotic diversity in groundwater ecosystems.

Go West: A One Way Stepping-Stone Dispersion Model for the Cavefish Lucifuga dentata in Western Cuba

Consistent with the limited dispersal capacity of most troglobitic animals, almost all Lucifuga cavefish species have very narrow geographic distribution in Cuba. However, one species, L. dentata, has a wide but disjointed distribution over 300 km in the west of the island. In order to estimate the relative role of vicariance and dispersal in the unexpected L. dentata distribution, we obtained partial sequences of the mitochondrial DNA (mtDNA) cytochrome b (cytb) gene and control region (CR), and then applied Approximate Bayesian Computation (ABC), based on the identification of five genetic and geographic congruent groups of populations. The process that best explains the distribution of genetic diversity in this species is sequential range expansion from east Matanzas to the western Pinar del Río provinces, followed by isolation of groups of populations. We found relative high haplotype diversity and low nucleotide diversity in all but the Havana group, which has high values for both diversity parameters, suggesting that this group has been demographically stable over time. For two groups of populations (Cayuco and Bolondrón), the mismatch distribution analyses suggests past demographic expansion. In the case of the Cayuco region, the star like relationships of haplotypes in the network suggests a recent founding event, congruent with other evidence indicating that this is the most recently colonized region. Over all, the results suggest that a combination of habitat availability, temporal interconnections, and possibly the biological properties of this species, may have enabled its dispersal and range expansion compared to other species of the genus, which are more geographically restricted.

Molecular phylogenetic analyses reveal a new southern hemisphere oniscidean family (Crustacea : Isopoda) with a unique water transport system

A significant diversity of terrestrial oniscidean isopods was recently discovered in the subterranean ‘calcrete islands’ of Western Australia, but the species and higher-level systematic status of much of the fauna are currently uncertain. Here we focus on one group of species that was initially assigned to the genus Trichorhina (Platyarthridae), based on several shared characters, and investigate the phylogenetic relationships of these species to 21 oniscidean genera, including 13 known families, using 18S rDNA sequence data. We then present phylogenetic analyses using 28S-only and combined 18S, 28S rDNA and mitochondrial cytochrome c oxidase subunit I (COI) data for a more restricted sampling of taxa, and present results for a detailed morphological study of the antennae and other cephalic structures of exemplar taxa. Bayesian and maximum likelihood analyses of the extended 18S-only, the 28S-only and multi-gene datasets provide strong evidence for a distinct well-supported monophyletic group comprising the new Western Australian and one South American taxon. This clade is unrelated to all included members of Platyarthridae, which appears to be polyphyletic, and it forms a distinct group relative to other oniscidean families. Given these findings and the results of the morphological study, a new southern hemisphere oniscidean family, Paraplatyarthridae Javidkar & King, fam. nov. is erected based on Paraplatyarthrus subterraneus Javidkar & King, gen. & sp. nov. (type genus and species), and several undescribed taxa which occur in the arid (terrestrial and subterranean) regions of Western Australia and subtropical South America. Paraplatyarthridae is distinguishable from all other oniscidian families on a combination of character states including, among others, the presence of fan-like scale setae on the dorsal body, and the ventral second antenna with leaf-like scale setae and a furrow containing elongated hair-like capillary setae that form part of a water conducting system unique within Oniscidea. This study has important implications for the higher-level classification of oniscidean crustaceans and points to the need for a more detailed molecular phylogeny that includes a comprehensive sampling of southern hemisphere taxa.

Opsin transcripts of predatory diving beetles: a comparison of surface and subterranean photic niches

The regressive evolution of eyes has long intrigued biologists yet the genetic underpinnings remain opaque. A system of discrete aquifers in arid Australia provides a powerful comparative means to explore trait regression at the genomic level. Multiple surface ancestors from two tribes of diving beetles (Dytiscidae) repeatedly invaded these calcrete aquifers and convergently evolved eye-less phenotypes. We use this system to assess transcription of opsin photoreceptor genes among the transcriptomes of two surface and three subterranean dytiscid species and test whether these genes have evolved under neutral predictions. Transcripts for UV, long-wavelength and ciliary-type opsins were identified from the surface beetle transcriptomes. Two subterranean beetles showed parallel loss of all opsin transcription, as expected under 'neutral' regressive evolution. The third species Limbodessus palmulaoides retained transcription of a long-wavelength opsin (lwop) orthologue, albeit in an aphotic environment. Tests of selection on lwop indicated no significant differences between transcripts derived from surface and subterranean habitats, with strong evidence for purifying selection acting on L. palmulaoides lwop. Retention of sequence integrity and the lack of evidence for neutral evolution raise the question of whether we have identified a novel pleiotropic role for lwop, or an incipient phase of pseudogene development.

Opsin transcripts of predatory diving beetles: a comparison of surface and subterranean photic niches

The regressive evolution of eyes has long intrigued biologists yet the genetic underpinnings remain opaque. A system of discrete aquifers in arid Australia provides a powerful comparative means to explore trait regression at the genomic level. Multiple surface ancestors from two tribes of diving beetles (Dytiscidae) repeatedly invaded these calcrete aquifers and convergently evolved eye-less phenotypes. We use this system to assess transcription of opsin photoreceptor genes among the transcriptomes of two surface and three subterranean dytiscid species and test whether these genes have evolved under neutral predictions. Transcripts for UV, long-wavelength and ciliary-type opsins were identified from the surface beetle transcriptomes. Two subterranean beetles showed parallel loss of all opsin transcription, as expected under 'neutral' regressive evolution. The third species Limbodessus palmulaoides retained transcription of a long-wavelength opsin (lwop) orthologue, albeit in an aphotic environment. Tests of selection on lwop indicated no significant differences between transcripts derived from surface and subterranean habitats, with strong evidence for purifying selection acting on L. palmulaoides lwop. Retention of sequence integrity and the lack of evidence for neutral evolution raise the question of whether we have identified a novel pleiotropic role for lwop, or an incipient phase of pseudogene development.

Can environment predict cryptic diversity? The case of Niphargus inhabiting Western Carpathian groundwater

In the last decade, several studies have shown that subterranean aquatic habitats harbor cryptic species with restricted geographic ranges, frequently occurring as isolated populations. Previous studies on aquatic subterranean species have implied that habitat heterogeneity can promote speciation and that speciation events can be predicted from species’ distributions. We tested the prediction that species distributed across different drainage systems and karst sectors comprise sets of distinct species. Amphipods from the genus Niphargus from 11 caves distributed along the Western Carpathians (Romania) were investigated using three independent molecular markers (COI, H3 and 28S). The results showed that: 1) the studied populations belong to eight different species that derive from two phylogenetically unrelated Niphargus clades; 2) narrow endemic species in fact comprise complexes of morphologically similar species that are indistinguishable without using a molecular approach. The concept of monophyly, concordance between mitochondrial and nuclear DNA, and the value of patristic distances were used as species delimitation criteria. The concept of cryptic species is discussed within the framework of the present work and the contribution of these species to regional biodiversity is also addressed.

Microsatellite development and first population size estimates for the groundwater isopod Proasellus walteri

Effective population size (N e) is one of the most important parameters in, ecology, evolutionary and conservation biology; however, few studies of N e in surface freshwater organisms have been published to date. Even fewer studies have been carried out in groundwater organisms, although their evolution has long been considered to be particularly constrained by small N e. In this study, we estimated the contemporary effective population size of the obligate groundwater isopod: Proaselluswalteri (Chappuis, 1948). To this end, a genomic library was enriched for microsatellite motifs and sequenced using 454 GS-FLX technology. A total of 54,593 reads were assembled in 10,346 contigs or singlets, of which 245 contained candidate microsatellite sequences with suitable priming sites. Ninety-six loci were tested for amplification, polymorphism and multiplexing properties, of which seven were finally selected for N e estimation. Linkage disequilibrium and approximate Bayesian computation methods revealed that N e in this small interstitial groundwater isopod could reach large sizes (> 585 individuals). Our results suggest that environmental conditions in groundwater, while often referred to as extreme, are not necessarily associated with small N e.

Evolutionary refugia and ecological refuges: key concepts for conserving Australian arid zone freshwater biodiversity under climate change

Refugia have been suggested as priority sites for conservation under climate change because of their ability to facilitate survival of biota under adverse conditions. Here, we review the likely role of refugial habitats in conserving freshwater biota in arid Australian aquatic systems where the major long-term climatic influence has been aridification. We introduce a conceptual model that characterizes evolutionary refugia and ecological refugees based on our review of the attributes of aquatic habitats and freshwater taxa (fishes and aquatic invertebrates) in arid Australia. We also identify methods of recognizing likely future refugia and approaches to assessing the vulnerability of arid-adapted freshwater biota to a warming and drying climate. Evolutionary refugia in arid areas are characterized as permanent, groundwater-dependent habitats (subterranean aquifers and springs) supporting vicariant relicts and short-range endemics. Ecological refugees can vary across space and time, depending on the dispersal abilities of aquatic taxa and the geographical proximity and hydrological connectivity of aquatic habitats. The most important are the perennial waterbodies (both groundwater and surface water fed) that support obligate aquatic organisms. These species will persist where suitable habitats are available and dispersal pathways are maintained. For very mobile species (invertebrates with an aerial dispersal phase) evolutionary refugia may also act as ecological refugees. Evolutionary refugia are likely future refugia because their water source (groundwater) is decoupled from local precipitation. However, their biota is extremely vulnerable to changes in local conditions because population extinction risks cannot be abated by the dispersal of individuals from other sites. Conservation planning must incorporate a high level of protection for aquifers that support refugial sites. Ecological refuges are vulnerable to changes in regional climate because they have little thermal or hydrological buffering. Accordingly, conservation planning must focus on maintaining meta-population processes, especially through dynamic connectivity between aquatic habitats at a landscape scale.

Patterns of population genetic variation in sympatric chiltoniid amphipods within a calcrete aquifer reveal a dynamic subterranean environment

Calcrete aquifers from the Yilgarn region of arid central Western Australia contain an assemblage of obligate groundwater invertebrate species that are each endemic to single aquifers. Fine-scale phylogeographic and population genetic analyses of three sympatric and independently derived species of amphipod (Chiltoniidae) were carried out to determine whether there were common patterns of population genetic structure or evidence for past geographic isolation of populations within a single calcrete aquifer. Genetic diversity in amphipod mitochondrial DNA (cytochrome c oxidase subunit I gene) and allozymes were examined across a 3.5 km(2) region of the Sturt Meadows calcrete, which contains a grid of 115 bore holes (=wells). Stygobiont amphipods were found to have high levels of mitochondrial haplotype diversity coupled with low nucleotide diversity. Mitochondrial phylogeographic structuring was found between haplogroups for one of the chiltoniid species, which also showed population structuring for nuclear markers. Signatures of population expansion in two of the three species, match previous findings for diving beetles at the same site, indicating that the system is dynamic. We propose isolation of populations in refugia within the calcrete, followed by expansion events, as the most likely source of intraspecific genetic diversity, due to changes in water level influencing gene flow across the calcrete.

Evolution of blind beetles in isolated aquifers: a test of alternative modes of speciation

Evidence is growing that not only allopatric but also sympatric speciation can be important in the evolution of species. Sympatric speciation has most convincingly been demonstrated in laboratory experiments with bacteria, but field-based evidence is limited to a few cases. The recently discovered plethora of subterranean diving beetle species in isolated aquifers in the arid interior of Australia offers a unique opportunity to evaluate alternative modes of speciation. This naturally replicated evolutionary experiment started 10-5 million years ago, when climate change forced the surface species to occupy geographically isolated subterranean aquifers. Using phylogenetic analysis, we determine the frequency of aquifers containing closely related sister species. By comparing observed frequencies with predictions from different statistical models, we show that it is very unlikely that the high number of sympatrically occurring sister species can be explained by a combination of allopatric evolution and repeated colonisations alone. Thus, diversification has occurred within the aquifers and likely involved sympatric, parapatric and/or microallopatric speciation.

Species diversity and genetic differentiation of stygofauna (Syncarida:Bathynellacea) across an alluvial aquifer in north-eastern Australia

Recent research suggests that alluvial aquifers in southern and eastern Australia may contain a diverse subterranean aquatic fauna (i.e. stygofauna). However, to date only a limited number of alluvial aquifers have been studied and little molecular data are available to assess species-level diversity and spatial patterns of genetic variation within stygofaunal species. In this paper, we present the initial results of a stygofaunal survey of the Burdekin River alluvial aquifer in Queensland, extending the northern range of alluvial aquifers along the east coast of Australia that have been investigated. The survey resulted in the collection of bathynellid stygofauna (Syncarida: Bathynellacea) and genetic analyses were conducted to determine species level diversity using the mitochondrial cytochrome oxidase subunit I (COI) gene. We further investigated the phylogenetic relationships of the species with bathynellids from western and southern Australia to assess the generic status of species. Four highly divergent COI lineages within the Parabathynellidae and one lineage within the Bathynellidae were found. These lineages did not group within any described genera, and phylogenetic analyses indicated that both local radiations and the retention of a lineage that was more apical in the genealogy account for the diversity within the Parabathynellidae in the Burdekin River alluvial aquifer. Most COI lineages were sampled from only a single bore, although one taxon within the Parabathynellidae was found to be more widespread in the aquifer. Haplotypes within this taxon were not shared among bores (ΦST = 0.603, P < 0.001). Overall, the high species diversity for bathynellaceans from an alluvial aquifer reported here, and surveys of bathynellaceans in several other alluvial systems in south-eastern Australia, suggests that groundwater ecosystems of eastern Australia may contain high stygofaunal diversity by Australian and world standards, particularly at the generic level for parabathynellids.

Explosive radiation of the genus Schizopera on a small subterranean island in Western Australia (Copepoda:Harpacticoida): unravelling the cases of cryptic speciation, size differentiation and multiple invasions

A previously unsurveyed calcrete aquifer in the Yilgarn region of Western Australia revealed an unprecedented diversity of copepods, representing 67% of that previously recorded in this whole region. Especially diverse was the genus Schizopera, with up to four morphospecies per bore and a significant size difference between them. Aims of this study were to: (1) survey the extent of this diversity using morphological and molecular tools; (2) derive a molecular phylogeny based on COI; and (3) investigate whether high diversity is a result of an explosive radiation, repeated colonisations, or both, size differentiation is a result of parallel evolution or different phylogeny, and whether Schizopera is a recent invasion in inland waters. More than 300 samples were analysed and the COI fragment successfully amplified by PCR from 43 specimens. Seven species and one subspecies are described as new, and three possible cryptic species were detected. Reconstructed phylogenies reveal that both explosive radiation and multiple colonisations are responsible for this richness, and that Schizopera is probably a recent invasion in these habitats. No evidence for parallel evolution was found, interspecific size differentiation being a result of different phylogeny. Sister species have parapatric distributions and show niche partitioning in the area of overlap.

Molecular phylogeny and phylogeography of the Cuban cave-fishes of the genus Lucifuga: evidence for cryptic allopatric diversity

Underground environments are increasingly recognized as reservoirs of faunal diversity. Extreme environmental conditions and limited dispersal ability of underground organisms have been acknowledged as important factors promoting divergence between species and conspecific populations. However, in many instances, there is no correlation between genetic divergence and morphological differentiation. Lucifuga Poey is a stygobiotic fish genus that lives in Cuban and Bahamian caves. In Cuba, it offers a unique opportunity to study the influence of habitat fragmentation on the genetic divergence of stygobiotic species and populations. The genus includes four species and one morphological variant that have contrasting geographical distributions. In this study, we first performed a molecular phylogenetic analysis of the Lucifuga Cuban species using mitochondrial and nuclear markers. The mitochondrial phylogeny revealed three deeply divergent clades that were supported by nuclear and morphological characters. Within two of these main clades, we identified five lineages that are candidate cryptic species and a taxonomical synonymy between Lucifuga subterranea and Lucifuga teresinarum. Secondly, phylogeographic analysis using a fragment of the cytochrome b gene was performed for Lucifuga dentata, the most widely distributed species. We found strong geographical organization of the haplotype clades at different geographic scales that can be explained by episodes of dispersal and population expansion followed by population fragmentation and restricted gene flow. At a larger temporal scale, these processes could also explain the diversification and the distribution of the different species.

Evidence for population fragmentation within a subterranean aquatic habitat in the Western Australian desert

The evolution of subterranean animals following multiple colonisation events from the surface has been well documented, but few studies have investigated the potential for species diversification within cavernicolous habitats. Isolated calcrete (carbonate) aquifers in central Western Australia have been shown to contain diverse assemblages of aquatic subterranean invertebrate species (stygofauna) and to offer a unique model system for exploring the mechanisms of speciation in subterranean ecosystems. In this paper, we investigated the hypothesis that microallopatric speciation processes (fragmentation and isolation by distance (IBD)) occur within calcretes using a comparative phylogeographic study of three stygobiontic diving beetle species, one amphipod species and a lineage of isopods. Specimens were sequenced for the mitochondrial cytochrome c oxidase 1 gene from three main sites: Quandong Well, Shady Well (SW) and Mt. Windarra (MW), spanning a 15 km region of the Laverton Downs Calcrete. Phylogenetic and haplotype network analyses revealed that each species possessed a single divergent clade of haplotypes that were present only at the southern MW site, despite the existence of other haplotypes at MW that were shared with SW. IBD between MW and SW was evident, but the common phylogeographic pattern most likely resulted from fragmentation, possibly by a salt lake adjacent to MW. These findings suggest that microallopatric speciation within calcretes may be a significant diversifying force, although the proportion of stygofauna species that may have resulted from in situ speciation in this system remains to be determined.

Evolution underground: shedding light on the diversification of subterranean insects

A recent study in BMC Evolutionary Biology has reconstructed the molecular phylogeny of a large Mediterranean cave-dwelling beetle clade, revealing an ancient origin and strong geographic structuring. It seems likely that diversification of this clade in the Oligocene was seeded by an ancestor already adapted to subterranean life.

VIEWPOINT. Is the Australian subterranean fauna uniquely diverse?

Australia was historically considered a poor prospect for subterranean fauna but, in reality, the continent holds a great variety of subterranean habitats, with associated faunas, found both in karst and non-karst environments. This paper critically examines the diversity of subterranean fauna in several key regions for the mostly arid western half of Australia. We aimed to document levels of species richness for major taxon groups and examine the degree of uniqueness of the fauna. We also wanted to compare the composition of these ecosystems, and their origins, with other regions of subterranean diversity world-wide. Using information on the number of ‘described’ and ‘known’ invertebrate species (recognised based on morphological and/or molecular data), we predict that the total subterranean fauna for the western half of the continent is 4140 species, of which ~10% is described and 9% is ‘known’ but not yet described. The stygofauna, water beetles, ostracods and copepods have the largest number of described species, while arachnids dominate the described troglofauna. Conversely, copepods, water beetles and isopods are the poorest known groups with less than 20% described species, while hexapods (comprising mostly Collembola, Coleoptera, Blattodea and Hemiptera) are the least known of the troglofauna. Compared with other regions of the world, we consider the Australian subterranean fauna to be unique in its diversity compared with the northern hemisphere for three key reasons: the range and diversity of subterranean habitats is both extensive and novel; direct faunal links to ancient Pangaea and Gondwana are evident, emphasising their early biogeographic history; and Miocene aridification, rather than Pleistocene post-ice age driven diversification events (as is predicted in the northern hemisphere), are likely to have dominated Australia’s subterranean speciation explosion. Finally, we predict that the geologically younger, although more poorly studied, eastern half of the Australian continent is unlikely to be as diverse as the western half, except for stygofauna in porous media. Furthermore, based on similar geology, palaeogeography and tectonic history to that seen in the western parts of Australia, southern Africa, parts of South America and India may also yield similar subterranean biodiversity to that described here.