Evolution of endoglucanase genes in subterranean and surface isopod crustaceans from Central Western Australia

Recent studies have identified a significant number of endogenous cellulase genes in various arthropods, including isopods, allowing them to process hydrocarbons efficiently as a food source. While this research has provided insight into underlying gene-level processes in cellulose decomposition by arthropods, little is known about the existence and expression of cellulase genes in species from cave environments where carbohydrates are sparse. To investigate whether endogenous cellulase genes are maintained in subterranean species, we sequenced the transcriptomes of two subterranean paraplatyarthrid isopod species from calcrete (carbonate) aquifers of central Western Australia and a related surface isopod species. Seven protein-coding open-reading frames associated with endoglucanase genes were identified in all species. Orthology inference analyses, using a wide range of cellulase sequences from available databases, supported the endogenous origin of the putative endoglucanase genes. Selection analyses revealed that these genes are primarily subject to purifying selection in most of the sites for both surface and subterranean isopod species, indicating that they are likely to encode functional peptides. Furthermore, evolutionary branch models supported the hypothesis of an adaptive shift in selective pressure acting on the subterranean lineages compared with the ancestral lineage and surface species. Branch-site models also revealed a few amino acid sites on the subterranean branches to be under positive selection, suggesting the acquisition of novel adaptations to the subterranean environments. These findings also imply that hydrocarbons exist in subsurface aquifers, albeit at reduced levels, and have been utilized by subterranean isopods as a source of energy for millions of years.

Phylogenomics and biogeography of arid-adapted Chlamydogobius goby fishes

The progressive aridification of the Australian continent from ∼ 20 million years ago posed severe challenges for the persistence of its resident biota. A key question involves the role of refugial habitats - specifically, their ability to mediate the effects of habitat loss and fragmentation, and their potential to shape opportunities for allopatric speciation. With freshwater species, for example, the patchiness, or absence, of water will constrain distributions. However, aridity may not necessarily isolate populations if disjunct refugia experience frequent hydrological connections. To investigate this potential dichotomy, we explored the evolutionary history of the Chlamydogobius gobies (Gobiiformes: Gobiidae), an arid-adapted genus of six small, benthic fish species that exploit all types of waterbodies (i.e. desert springs, waterholes and bore-fed wetlands, coastal estuarine creeks and mangroves) across parts of central and northern Australia. We used Anchored Phylogenomics to generate a highly resolved phylogeny of the group from sequence data for 260 nuclear loci. Buttressed by companion allozyme and mtDNA datasets, our molecular findings infer the diversification of Chlamydogobius in arid Australia, and provide a phylogenetic structure that cannot be simply explained by invoking allopatric speciation events reflecting current geographic proximity. Our findings are generally consistent with the existing morphological delimitation of species, with one exception: at the shallowest nodes of phylogenetic reconstruction, the molecular data do not fully support the current dichotomous delineation of C. japalpa from C. eremius in Kati Thanda-Lake Eyre-associated waterbodies. Together these findings illustrate the ability of structural (hydrological) connections to generate patterns of connectivity and isolation for an ecologically moderate disperser in response to ongoing habitat aridification. Finally, we explore the implications of these results for the immediate management of threatened (C. gloveri) and critically endangered (C. micropterus, C. squamigenus) congeners.

Differential transcriptomic responses to heat stress in surface and subterranean diving beetles

Subterranean habitats are generally very stable environments, and as such evolutionary transitions of organisms from surface to subterranean lifestyles may cause considerable shifts in physiology, particularly with respect to thermal tolerance. In this study we compared responses to heat shock at the molecular level in a geographically widespread, surface-dwelling water beetle to a congeneric subterranean species restricted to a single aquifer (Dytiscidae: Hydroporinae). The obligate subterranean beetle Paroster macrosturtensis is known to have a lower thermal tolerance compared to surface lineages (CT_max 38 °C cf. 42–46 °C), but the genetic basis of this physiological difference has not been characterized. We experimentally manipulated the thermal environment of 24 individuals to demonstrate that both species can mount a heat shock response at high temperatures (35 °C), as determined by comparative transcriptomics. However, genes involved in these responses differ between species and a far greater number were differentially expressed in the surface taxon, suggesting it can mount a more robust heat shock response; these data may underpin its higher thermal tolerance compared to subterranean relatives. In contrast, the subterranean species examined not only differentially expressed fewer genes in response to increasing temperatures, but also in the presence of the experimental setup employed here alone. Our results suggest P. macrosturtensis may be comparatively poorly equipped to respond to both thermally induced stress and environmental disturbances more broadly. The molecular findings presented here have conservation implications for P. macrosturtensis and contribute to a growing narrative concerning weakened thermal tolerances in obligate subterranean organisms at the molecular level.

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.

Development and evaluation of a custom bait design based on 469 single-copy protein-coding genes for exon capture of isopods (Philosciidae: Haloniscus)

Transcriptome-based exon capture approaches, along with next-generation sequencing, are allowing for the rapid and cost-effective production of extensive and informative phylogenomic datasets from non-model organisms for phylogenetics and population genetics research. These approaches generally employ a reference genome to infer the intron-exon structure of targeted loci and preferentially select longer exons. However, in the absence of an existing and well-annotated genome, we applied this exon capture method directly, without initially identifying intron-exon boundaries for bait design, to a group of highly diverse Haloniscus (Philosciidae), paraplatyarthrid and armadillid isopods, and examined the performance of our methods and bait design for phylogenetic inference. Here, we identified an isopod-specific set of single-copy protein-coding loci, and a custom bait design to capture targeted regions from 469 genes, and analysed the resulting sequence data with a mapping approach and newly-created post-processing scripts. We effectively recovered a large and informative dataset comprising both short (<100 bp) and longer (>300 bp) exons, with high uniformity in sequencing depth. We were also able to successfully capture exon data from up to 16-year-old museum specimens along with more distantly related outgroup taxa, and efficiently pool multiple samples prior to capture. Our well-resolved phylogenies highlight the overall utility of this methodological approach and custom bait design, which offer enormous potential for application to future isopod, as well as broader crustacean, molecular studies.

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.

Genomic data reveal similar genetic differentiation in aquifer species with different dispersal capabilities and life histories

Little is known about the life histories, genetic structure and population connectivity of shallow groundwater organisms. We used next-generation sequencing (RAD-seq) to analyse population genomic structure in two aquifer species: Paraperla frontalis (Banks, 1902), a stonefly with groundwater larvae and aerial (winged) adults; and Stygobromus sp., a groundwater-obligate amphipod. We found similar genetic differentiation in each species between floodplains separated by ~70 river km in the Flathead River basin of north-west Montana, USA. Given that Stygobromus lacks the above-ground life stage of P. frontalis, our findings suggest that connectivity and the magnitude of genetic structure cannot be definitively assumed from life history differences.

Molecular and taxonomic analyses in troglobiotic Alpioniscus (Illyrionethes) species from the Dinaric Karst (Isopoda: Trichoniscidae)

Species richness of terrestrial isopods is high in caves of the Dinaric Karst, which hosts ~10% of the world’s nominal oniscidean troglobionts. The most widespread taxon is the southern European genus Alpioniscus, which consists of two subgenera: Alpioniscus s.s. and Illyrionethes. Before this study, 14 nominal troglobiotic Illyrionethes taxa were recorded from the Dinaric Karst. Our molecular analyses using two mitochnodrial DNA (16S rRNA and COI) fragments and a nuclear gene (H3) fragment on all known Dinaric taxa identified three distinct lineages: strasseri-, heroldi- and magnus-lineage. Our results confirmed the validity of most nominal species. The exceptions are Alpioniscus balthasari, which consists of two different species including Alpioniscus iapodicus, and Alpioniscus heroldi, which is paraphyletic with respect to Alpioniscus bosniensis. The strasseri-lineage was highly supported by all phylogenetic methods used; therefore, we performed a detailed morphological analysis to distinguish and characterize the species of this group. New morphological characters, such as body part ratios, are proposed for future species identification. In addition, we redescribe three known species (Alpioniscus strasseri, Alpioniscus christiani and Alpioniscus balthasari) and describe two new ones (Alpioniscus hirci sp. nov. and Alpioniscus velebiticus sp. nov.). As a result, 15 nominal species of Illyrionethes are currently known from the Dinaric Karst.

Too hot to handle: Cenozoic aridification drives multiple independent incursions of Schizomida (Hubbardiidae) into hypogean environments

The formation of the Australian arid zone, Australia's largest and youngest major biome, has been recognized as a major driver of rapid evolutionary radiations in terrestrial plants and animals. Here, we investigate the phylogenetic diversity and evolutionary history of subterranean short-tailed whip scorpions (Schizomida: Hubbardiidae), which are a significant faunal component of Western Australian hypogean ecosystems. We sequenced two mitochondrial (12S, COI) and three nuclear DNA markers (18S, 28S, ITS2) from ∼600 specimens, largely from the genera Draculoides and Paradraculoides, including 20 previously named species and an additional 56 newly identified operational taxonomic units (OTUs). Phylogenetic analyses revealed a large and rapid species radiation congruent with Cenozoic aridification of the continent, in addition to the identification of a new genus in Western Australia and the first epigean schizomid from the Pilbara. Here, we also synonymise Paradraculoides with Draculoides (new synonymy), due to paraphyly and a lack of reliable characters to define the two genera. Our results are consistent with multiple colonisations of the subterranean realm from epigean ancestors as their forest habitat fragmented and retracted, with ongoing fragmentation and diversification of lineages underground. These findings illustrate the remarkable diversity and high incidence of short-range endemism of Western Australia's subterranean fauna, which has important implications for identifying and managing short-range endemic subterranean fauna. They also highlight the advantages of including molecular data in subterranean fauna surveys as all specimens can be utilized, regardless of sex and life stage. Additionally, we have provided the first multi-gene phylogenetic framework for Australian schizomids, which will enable researchers and environmental consultants to identify new taxa or align them to existing lineages.

The genus Alpioniscus Racovitza, 1908 in Sardinia: taxonomy and natural history (Isopoda, Oniscidea, Trichoniscidae)

The genus Alpioniscus Racovitza, 1908 (Trichoniscidae) from Sardinia is revised. Three new cave-dwelling species are described: A.onnisi Taiti & Argano, sp. n., A.stochi Taiti & Argano, sp. n., and A.sideralis Taiti & Argano, sp. n.. The genus Utopioniscus Schmalfuss, 2005 is considered to be a junior synonym of Alpioniscus, after morphological and molecular analyses. Alpioniscusfragilis (Budde-Lund, 1909) and A.kuehni from Grotta del Bue Marino are illustrated. With the new species, the genus Alpioniscus in Sardinia comprises six species: two troglobionts (A.fragilis and A.onnisi), one endogean and troglobiont (A.thanit Taiti & Argano, 2009), and three stygobionts (A.kuehni, A.stochi, and A.sideralis). All the species occur in karstic areas in the central-eastern and south-eastern part of the island. A key to all the Sardinian species of Alpioniscus is provided.

Integrative taxonomy reveals a new genus and new species of Philosciidae (Crustacea: Isopoda: Oniscidea) from the Neotropical region

Although new methods and data are conquering space in the field of taxonomy, such as integrative taxonomy, most terrestrial isopod species are still described based only on morphology. Species of the genus Atlantoscia Ferrara and Taiti, 1981 were the first and are the unique terrestrial isopods from the Neotropics for which a molecular phylogeny was already conducted. Previous results indicated that this genus could be paraphyletic, and a more detailed analysis would be required. Our aim was to reconstruct the phylogeny of Atlantoscia using mitochondrial and nuclear markers and test its monophyly by integrating molecular and morphological data. We observed that, indeed, Atlantoscia is paraphyletic. Atlantoscia ituberasensis Campos-Filho, Lisboa and Araujo, 2013 and Atlantoscia rubromarginata Araujo and Leistikow, 1999 were placed in a new genus of terrestrial isopods, Paratlantoscia gen. nov., together with a new species described in the present study, Paratlantoscia robusta sp. nov. The new genus is defined by the presence of specialized respiratory areas in the pleopod exopods and its validity is highly corroborated by molecular analyses and by biogeographic information. This study highlights the importance of multiple and complementary perspectives as a way to improve the quality of species hypothesis and associated descriptions.

Multiple molecular markers reinforce the systematic framework of unique Australian cave fishes (Milyeringa : Gobioidei)

Australia was once thought to be a biodiversity desert when considering the subterranean world; however, recent work has revealed a fascinating collection of cave creatures, many with surprising biogeographic histories. This has especially been so in the karstic regions of north-western Australia (Cape Range peninsula, Barrow Island, Pilbara), which is home not only to a diverse collection of subterranean invertebrates, but also to the continent’s only known underworld-adapted vertebrates, which includes the cave fish in the genus Milyeringa. These cave gudgeons have recently been in a state of taxonomic flux, with species being both split and lumped, but this was done with limited data (incomplete geographic sampling and no nuclear DNA sequence data). Therefore, we have revisited the systematic status of Milyeringa in a total-evidence molecular approach by integrating all existing data (mitochondrial, allozymes) with new DNA sequences from nuclear and mitochondrial loci and new multilocus allozyme data. Our conclusion, that there are two species, matches the most recent taxonomic treatment, with Milyeringa veritas present on both the eastern and western sides of the Cape Range peninsula, and Milyeringa justitia on Barrow Island. This has implications for future research in the linked fields of biogeography and conservation.

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.

Groundwater oligochaetes show complex genetic patterns of distribution in the Pilbara region of Western Australia

Patterns of genetic diversity in the groundwater fauna of Australia have largely focused on obligate stygobites of relatively large size, namely, crustaceans. Oligochaete worms, with their smaller size and broader ecological niches, provide a contrasting model in which to examine such patterns. Genetic diversity in subterranean oligochaetes in the Pilbara region of Western Australia were examined using one nuclear (18S) and two mitochondrial (COI, 12S) regions. The observed variation was assessed at three levels of hydrology – river basin, creek catchment, and individual bore or site – to document geographic patterns. Most species appeared to be restricted to an individual catchment; however, five species, representing three families, were widespread, with some haplotypes being shared between bores, catchments and even basins. General patterns suggest that while hydrology plays a role in the distribution of oligochaete species, it does not always confine them to catchments, in contrast to patterns observed in groundwater isopods and amphipods in the region. We suggest that intrinsic characteristics of oligochaetes, such as body size, shape, reproductive strategy and ecological requirements, may have allowed them greater dispersal within the subterranean biome of the Pilbara. In particular, oligochaetes may occupy subterranean and surface waters, increasing their opportunities for dispersal.

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.

Molecular evidence for mid-Pleistocene divergence of populations of three freshwater amphipod species (Talitroidea : Chiltoniidae) on Kangaroo Island, South Australia, with a new spring-associated genus and species

Recent molecular and morphological analyses have shown that chiltoniid amphipods, once thought to be a relictual group, are a diverse and speciose family of Australian freshwater amphipods. As part of a larger examination of the family, chiltoniids from Kangaroo Island in South Australia were collected and analysed using molecular (COI and 28S) and morphological methods in order to understand species distributional patterns and relationships. Kartachiltonia moodyi gen. nov., sp. nov., a spring-associated species endemic to the island, was discovered and populations of three additional mainland species (Austrochiltonia australis, A. dalhousiensis and A. subtenuis) were examined. The island populations of A. australis, A. dalhousiensis and A. subtenuis were found to form natural groups with differing haplotype coalescence times dating from the Early to Mid-Pleistocene. Numerous cycles of regional climate change throughout the Pleistocene are likely to have driven speciation in chiltoniid amphipods in southern Australia and the presence of multiple chiltoniid species at Kangaroo Island indicates that it exists at a likely convergence of species distribution patterns. Three possible hypotheses to explain the evolution and diversity of chiltoniids in southern Australia are discussed as are evidence for potential introduction and long-distance dispersal events.

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.

Molecular species boundaries in the phreatoicidean genus Amphisopus (Isopoda : Amphisopidae) and evidence for a new freshwater isopod species from Western Australia

The freshwater isopod genus Amphisopus is one of only two phreatoicidean genera in Western Australia with wide distributions and containing multiple described species. Two species (Amphisopus annectans and A. lintoni) are known from the south-western part of the state. With recent sampling extending the known range of Amphisopus and the recorded possibility of an undescribed species, this study aimed to examine genetic species boundaries and to detect additional species diversity. Isopods were sampled from across the range, and genetic structure was examined using mtDNA sequence data from a COI fragment and data from ten polymorphic allozyme loci. While allozyme data supported the clear separation of the known species, phylogenetic analyses presented three divergent monophyletic, geographically restricted clades occurring in the western, eastern and central parts of the Amphisopus distribution. These corresponded to A. annectans, A. lintoni and an undescribed species, respectively. The latter’s status was supported by its phylogenetic position relative to the known species and by the extent of differentiation observed among other phreatoicidean taxa. As with other taxa from the region, divergence times among these three species supported aridity in the Miocene–Pliocene as a driver of diversification, and suggested the existence of wetter refugial areas in south-western Australia.

Two new subterranean ameirids (Crustacea : Copepoda : Harpacticoida) expose weaknesses in the conservation of short-range endemics threatened by mining developments in Western Australia

The discovery of two new non-marine ameirids from the southern Yilgarn region significantly extends the geographic range for this group in Australia and exposes weaknesses in the conservation and environmental impact assessment (EIA) of subterranean species potentially threatened by mining developments. Megastygonitocrella embe, sp. nov. differs from seven previously described Australian congeners by the armature of the second leg endopod and absence of spinules on the somites. A key to world species of Megastygonitocrella is presented. Phylogenetic analysis based on 57 morphological characters and 30 species belonging to the Stygonitocrella s.l. group suggests that Antistygonitocrella pardalotos, gen. et sp. nov. has no close relatives anywhere in the world. Superficial similarities between the two new species are either plesiomorphies or homoplasies. The habitats of these new short-range endemic species are fractured-rock aquifers developed in Archaean greenstone, where the groundwater is characterised by acid conditions, high salinity and low dissolved oxygen. The population of A. pardalotos is threatened by a mining development. Despite the advanced level of environmental protection policy in Western Australia, our taxonomic study highlights limitations in EIA practices and discusses potential improvements which have global relevance in regions where short-range endemics coincide with extraction of mineral resources.

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.

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.

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.

Multiple biogeographical barriers identified across the monsoon tropics of northern Australia: phylogeographic analysis of the brachyotis group of rock-wallabies

The monsoon tropics of northern Australia are a globally significant biodiversity hotspot, but its phylogeography is poorly known. A major challenge for this region is to understand the biogeographical processes that have shaped the distribution and diversity of taxa, without detailed knowledge of past climatic and environmental fluctuations. Although molecular data have great potential to address these questions, only a few species have been examined phylogeographically. Here, we use the widely distributed and abundant short-eared rock-wallaby (Petrogale brachyotis; n = 101), together with the sympatric monjon (P. burbidgei; n = 11) and nabarlek (P. concinna; n = 1), to assess historical evolutionary and biogeographical processes in northern Australia. We sequenced ∼1000 bp of mitochondrial DNA (control region, ND2) and ∼3000 bp of nDNA (BRCA1, ω-globin and two anonymous loci) to investigate phylogeographic structuring and delineate the time-scale of diversification within the region. Our results indicate multiple barriers between the Top End (Northern Territory) and Kimberley (Western Australia), which have caused divergence throughout the Plio-Pleistocene. Eight geographically discrete and genetically distinct lineages within the brachyotis group were identified, five of which are separated by major river valleys (Ord, Victoria, Daly), arid lowlands and discontinuous sandstone ranges. It is likely that these barriers have similarly influenced genetic structure in other monsoonal biota.

Difficulties barcoding in the dark: the case of crustacean stygofauna from eastern Australia

The eastern Australian aquifers remain mostly unexplored; however, recent surveys suggest that there could be substantial levels of subterranean biodiversity hidden in these aquifers. Groundwater fauna (stygofauna) is often characterised by short-range endemism. Furthermore, high levels of cryptic species, and lack of formal taxonomic descriptions and taxonomic expertise for many of the groups demand innovative approaches for assessing subterranean biodiversity. Here we evaluate the potential of using DNA barcoding as a rapid biodiversity assessment tool for the subterranean groundwater fauna of New South Wales, Australia. We experienced low amplification success using universal and more taxon-specific primers for PCR amplification of the barcoding gene (COI) in a range of crustacean stygofauna. Sequence comparisons of the most commonly used COI universal primers in selected crustacean taxa revealed high levels of variability. Our results suggest that successful amplification of the COI region from crustacean stygofauna is not straightforward using the standard ‘universal’ primers. We propose that the development of a multiprimer (taxon specific) and multigene approach for DNA barcode analyses, using next-generation sequencing methodologies, will help to overcome many of the technical problems reported here and provide a basis for using DNA barcoding for rapid biodiversity assessments of subterranean aquatic ecosystems.

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.

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.

The Australian freshwater amphipods Austrochiltonia australis and Austrochiltonia subtenuis (Amphipoda:Talitroidea:Chiltoniidae) confirmed and two new cryptic Tasmanian species revealed using a combined molecular and morphological approach

Given the complex nature of freshwater catchment divides and emerging evidence of high levels of genetic diversity, there is great potential for cryptic species to exist among Australian freshwater amphipod groups. Among the chiltoniid amphipods, two congeneric species, Austrochiltonia australis (Sayce, 1901) and A. subtenuis (Sayce, 1902), have been widely recorded across southern Australia yet are poorly known and contentiously defined. A large fragment of the mitochondrial DNA cytochrome c oxidase I (COI) gene was examined and morphological diversity among populations assessed across the reported geographic range of the two putative species. The results confirmed A. australis and A. subtenuis as morphological and molecular species. In addition, two previously undetected and cryptic species from Tasmania are recognised – sister species to A. subtenuis and A. australis. Working conclusions provide evidence towards a more comprehensive systematic revision of the Chiltoniidae and present species information relevant to conservation and management efforts of Australian river systems. A key is presented to the chiltoniid amphipods of southern Australia.

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.

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

Calcrete aquifers in the arid Yilgarn region of central Western Australia are a biodiversity hotspot for stygofauna. A distinct pattern of interspecific size class variation among subterranean dytiscid beetle species has been observed in 29 of these aquifers where either two or three small, medium and/or large sympatric species are found that are in some cases sister species. We used a 3.5 km(2) grid of bores to sample dytiscids on a fine-scale and employed a comparative phylogeographical and population genetic approach to investigate the origins of a sympatric sister species triplet of diving beetles from a single aquifer. Mitochondrial DNA sequence data from the Cytochrome oxidase c subunit I gene revealed that all three species have high levels of haplotype diversity with ancient (approximately 1 million years ago) intra-specific coalescence of haplotypes, but low levels of nucleotide diversity. Population analyses provide evidence for multiple expansion events within each species. There was spatial heterogeneity in the distribution of genetic variation and abundance both within and among the three taxa. Population analyses revealed significant fine-scale differentiation with isolation by distance for Paroster macrosturtensis and P. mesosturtensis, but not the smallest species P. microsturtensis. Haplotype network analyses provided limited or no evidence for past population fragmentation within the large and small species, but substantial historical divergence was observed in P. mesosturtensis that was not spatially structured. A patchy population structure with contemporaneous and historical isolation by distance in the three species is likely to have been a significant isolating and diversifying force, preventing us from ruling out a potential role for allopatric divergence during speciation of this beetle sister triplet.