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

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.

Scratching the surface of subterranean biodiversity: Molecular analysis reveals a diverse and previously unknown fauna of Parabathynellidae (Crustacea: Bathynellacea) from the Pilbara, Western Australia

Like other crustacean families, the Parabathynellidae is a poorly studied subterranean and aquatic (stygobiontic) group in Australia, with many regions of available habitat having not yet been surveyed. Here we used a combined approach of molecular species delimitation methods, applied to mitochondrial and nuclear genetic data, to identify putative new species from material obtained from remote subterranean habitats in the Pilbara region of Western Australia. Based on collections from these new localities, we delineated a minimum of eight and up to 24 putative new species using a consensus from a range of molecular delineation methods and additional evidence. When we placed our new putative species into the broader phylogenetic framework of Australian Parabathynellidae, they grouped with two known genera and also within one new and distinct Pilbara-only clade. These new species significantly expand the known diversity of Parabathynellidae in that they represent a 22% increase to the 109 currently recognised species globally. Our investigations showed that sampling at new localities can yield extraordinary levels of new species diversity, with the majority of species showing likely restricted endemic geographical ranges. These findings represent only a small sample from a region comprising less than 2.5% of the Australian continent.

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.

Water Beetles as Models in Ecology and Evolution

Beetles have colonized water many times during their history, with some of these events involving extensive evolutionary radiations and multiple transitions between land and water. With over 13,000 described species, they are one of the most diverse macroinvertebrate groups in most nonmarine aquatic habitats and occur on all continents except Antarctica. A combination of wide geographical and ecological range and relatively accessible taxonomy makes these insects an excellent model system for addressing a variety of questions in ecology and evolution. Work on water beetles has recently made important contributions to fields as diverse as DNA taxonomy, macroecology, historical biogeography, sexual selection, and conservation biology, as well as predicting organismal responses to global change. Aquatic beetles have some of the best resolved phylogenies of any comparably diverse insect group, and this, coupled with recent advances in taxonomic and ecological knowledge, is likely to drive an expansion of studies in the future.

Geographic isolation facilitates the evolution of reproductive isolation and morphological divergence

Geographic isolation is known to contribute to divergent evolution, resulting in unique phenotypes. Oftentimes morphologically distinct populations are found to be interfertile while reproductive isolation is found to exist within nominal morphological species revealing the existence of cryptic species. These disparities can be difficult to predict or explain especially when they do not reflect an inferred history of common ancestry which suggests that environmental factors affect the nature of ecological divergence. A series of laboratory experiments and observational studies were used to address what role biogeographic factors may play in the ecological divergence of Hyalella amphipods. It was found that geographic isolation plays a key role in the evolution of reproductive isolation and divergent morphology and that divergence cannot be explained by molecular genetic variation.

Environmental DNA in subterranean biology: range extension and taxonomic implications for Proteus

Europe’s obligate cave-dwelling amphibian Proteus anguinus inhabits subterranean waters of the north-western Balkan Peninsula. Because only fragments of its habitat are accessible to humans, this endangered salamander’s exact distribution has been difficult to establish. Here we introduce a quantitative real time polymerase chain reaction-based environmental DNA (eDNA) approach to detect the presence of Proteus using water samples collected from karst springs, wells or caves. In a survey conducted along the southern limit of its known range, we established a likely presence of Proteus at seven new sites, extending its range to Montenegro. Next, using specific molecular probes to discriminate the rare black morph of Proteus from the closely related white morph, we detected its eDNA at five new sites, thus more than doubling the known number of sites. In one of these we found both black and white Proteus eDNA together. This finding suggests that the two morphs may live in contact with each other in the same body of groundwater and that they may be reproductively isolated species. Our results show that the eDNA approach is suitable and efficient in addressing questions in biogeography, evolution, taxonomy and conservation of the cryptic subterranean fauna.

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.

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.

Desert springs: deep phylogeographic structure in an ancient endemic crustacean (Phreatomerus latipes)

Desert mound springs of the Great Artesian Basin in central Australia maintain an endemic fauna that have historically been considered ubiquitous throughout all of the springs. Recent studies, however, have shown that several endemic invertebrate species are genetically highly structured and contain previously unrecognised species, suggesting that individuals may be geographically 'stranded in desert islands'. Here we further tested the generality of this hypothesis by conducting genetic analyses of the obligate aquatic phreatoicid isopod Phreatomerus latipes. Phylogenetic and phylogeographic relationships amongst P. latipes individuals were examined using a multilocus approach comprising allozymes and mtDNA sequence data. From the Lake Eyre region in South Australia we collected data for 476 individuals from 69 springs for the mtDNA gene COI; in addition, allozyme electrophoresis was conducted on 331 individuals from 19 sites for 25 putative loci. Phylogenetic and population genetic analyses showed three major clades in both allozyme and mtDNA data, with a further nine mtDNA sub-clades, largely supported by the allozymes. Generally, each of these sub-clades was concordant with a traditional geographic grouping known as spring complexes. We observed a coalescent time between ∼2-15 million years ago for haplotypes within each of the nine mtDNA sub-clades, whilst an older total time to coalescence (>15 mya) was observed for the three major clades. Overall we observed that multiple layers of phylogeographic history are exemplified by Phreatomerus, suggesting that major climate events and their impact on the landscape have shaped the observed high levels of diversity and endemism. Our results show that this genus reflects a diverse fauna that existed during the early Miocene and appears to have been regionally restricted. Subsequent aridification events have led to substantial contraction of the original habitat, possibly over repeated Pleistocene ice age cycles, with P. latipes populations becoming restricted in the distribution to desert springs.

What lies beneath: molecular phylogenetics and ancestral state reconstruction of the ancient subterranean Australian Parabathynellidae (Syncarida, Crustacea)

The crustacean family Parabathynellidae is an ancient and significant faunal component of subterranean ecosystems. Molecular data were generated in order to examine phylogenetic relationships amongst Australian genera and assess the species diversity of this group within Australia. We also used the resultant phylogenetic framework, in combination with an ancestral state reconstruction (ASR) analysis, to explore the evolution of two key morphological characters (number of segments of the first and second antennae), previously used to define genera, and assess the oligomerization principle (i.e. serial appendage reduction over time), which is commonly invoked in crustacean systematics. The ASR approach also allowed an assessment of whether there has been convergent evolution of appendage numbers during the evolution of Australian parabathynellids. Sequence data from the mtDNA COI and nDNA 18S rRNA genes were obtained from 32 parabathynellid species (100% of described genera and ~25% of described species) from key groundwater regions across Australia. Phylogenetic analyses revealed that species of each known genus, defined by traditional morphological methods, were monophyletic, suggesting that the commonly used generic characters are robust for defining distinct evolutionary lineages. Additionally, ancestral state reconstruction analysis provided evidence for multiple cases of convergent evolution for the two morphological characters evaluated, suggesting that caution needs to be shown when using these characters for elucidating phylogenetic relationships, particularly when there are few morphological characters available for reconstructing relationships. The ancestral state analysis contradicted the conventional view of parabathynellid evolution, which assumes that more simplified taxa (i.e. those with fewer-segmented appendages and setae) are derived and more complex taxa are primitive.

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.