Elucidating stygofaunal trophic web interactions via isotopic ecology

Groundwater is a hidden global keystone ecosystem

Groundwater is a vital ecosystem of the global water cycle, hosting unique biodiversity and providing essential services to societies. Despite being the largest unfrozen freshwater resource, in a period of depletion by extraction and pollution, groundwater environments have been repeatedly overlooked in global biodiversity conservation agendas. Disregarding the importance of groundwater as an ecosystem ignores its critical role in preserving surface biomes. To foster timely global conservation of groundwater, we propose elevating the concept of keystone species into the realm of ecosystems, claiming groundwater as a keystone ecosystem that influences the integrity of many dependent ecosystems. Our global analysis shows that over half of land surface areas (52.6%) has a medium-to-high interaction with groundwater, reaching up to 74.9% when deserts and high mountains are excluded. We postulate that the intrinsic transboundary features of groundwater are critical for shifting perspectives towards more holistic approaches in aquatic ecology and beyond. Furthermore, we propose eight key themes to develop a science-policy integrated groundwater conservation agenda. Given ecosystems above and below the ground intersect at many levels, considering groundwater as an essential component of planetary health is pivotal to reduce biodiversity loss and buffer against climate change.

Microbial diversity and geochemistry of groundwater impacted by steel slag leachates

To understand long-term impacts of steel slag material on aquifer geochemistry and microbial communities, we conducted four sampling campaigns in the Gier alluvial groundwater (Loire, France). In its northern part, the aquifer flows under a 200,000 m³ steel slag exhibiting high levels of chromium and molybdenum. Geochemical analyses of the water table revealed the existence of water masses with different chemical signatures. They allowed us to identify an area particularly contaminated by leachates from the slag heap, whatever the sampling period. Water samples from this area were compared to non-contaminated samples, with geochemical characteristics similar to the river samples. To follow changes in microbial communities, the V3-V4 region of 16 s rRNA gene was sequenced. Overall, we observed lower diversity indices in contaminated areas, with higher relative abundances of Verrucomicrobiota and Myxococcota phyla, while several Proteobacteria orders exhibited lower relative abundances. In particular, one single genus among the Verrucomicrobiota, Candidatus Omnitrophus, represented up to 36 % of total taxon abundance in areas affected by steel slag leachates. A large proportion of taxa identified in groundwater were also detected in the upstream river, indicating strong river-groundwater interactions. Our findings pave the way for future research work on C. Omnitrophus remediation capacities.

The biogeochemical responses of hyporheic groundwater to the long-run managed aquifer recharge: Linking microbial communities to hydrochemistry and micropollutants

Interactions of surface water and groundwater (SW-GW) in hyporheic zones produce biogeochemical hotspots. However, response patterns of hyporheic groundwater to external influences remain unclear. In this study, three datasets (hydrochemistry, antibiotics, and microbiome) were collected over a hydrological year to explore the influence of a 12-year managed aquifer recharge (MAR) project. We observed that the long-term MAR practice elevated nutrient and antibiotic levels while reduced redox potential in hyporheic groundwater, and these impacts depended on decreasing SW-GW interaction intensity with aquifer depth. In contrast, the long-term MAR practice increased community dissimilarity of 30-m groundwater but had little impact on 50-m or 80-m groundwater. Moreover, hyporheic community assembly was dominated by dispersal limitation, and thereby co-varied hydrochemistry and antibiotics only attributed to small community variability. The long-term MAR practice decreased species-interaction intensity and changed the abundance of metabolic functions in hyporheic groundwater. Furthermore, predicted community functions involving carbon, nitrogen, sulfur, and manganese cycles for 30-m groundwater showed higher abundances than those for 50- and 80-m groundwater. Collectively, we showed that hyporheic groundwater was sensitive to the SW-GW interaction and human activities, with the interactions of hydrochemistry, contaminants, and microbiome linking to hyporheic groundwater quality and ecosystem functioning.

Subterranean carbon flows from source to stygofauna: a case study on the atyid shrimp Stygiocaris stylifera (Holthuis, 1960) from Barrow Island (WA)

Groundwater biota are crucial for the ecological functioning of subterranean ecosystems. However, while knowledge of the taxonomic diversity of groundwater invertebrates (stygofauna) is increasing, functional ecological information is still limited. Here, we investigate seldom empirically tested assumptions around stygofaunal trophic plasticity in coping with oligotrophic habitats. We focus on Barrow Island (Western Australia), an ideal natural laboratory due to the occurrence of natural oil seeps in association with aquifers. The trophic position and food source use of the endemic atyid shrimp Stygiocaris stylifera (Holthuis, 1960) were assessed via δ¹³C and δ¹⁵N stable isotope analysis (SIA). Background information on the environmental conditions was gathered through hydrochemical data and δ¹³C SIA combined with ¹⁴C data from dissolved inorganic/organic carbon and particulate organic carbon from groundwater samples. Our results indicate carbon enrichment in proximity to the natural oil seepage coupled with changes in trophic positions of S. stylifera from higher consumers/predators to biofilm grazers/decomposers. These results are consistent with an increased involvement of hydrocarbon seeps and associated microbial communities in the carbon flows and confirm potential for the trophic flexibility in stygofauna. Further investigations involving other trophic groups will help elucidate the functioning of the ecosystems at a community level.

eDNA in subterranean ecosystems: Applications, technical aspects, and future prospects

Monitoring of biota is pivotal for the assessment and conservation of ecosystems. Environments worldwide are being continuously and increasingly exposed to multiple adverse impacts, and the accuracy and reliability of the biomonitoring tools that can be employed shape not only the present, but more importantly, the future of entire habitats. The analysis of environmental DNA (eDNA) metabarcoding data provides a quick, affordable, and reliable molecular approach for biodiversity assessments. However, while extensively employed in aquatic and terrestrial surface environments, eDNA-based studies targeting subterranean ecosystems are still uncommon due to the lack of accessibility and the cryptic nature of these environments and their species. Recent advances in genetic and genomic analyses have established a promising framework for shedding new light on subterranean biodiversity and ecology. To address current knowledge and the future use of eDNA methods in groundwaters and caves, this review explores conceptual and technical aspects of the application and its potential in subterranean systems. We briefly introduce subterranean biota and describe the most used traditional sampling techniques. Next, eDNA characteristics, application, and limitations in the subsurface environment are outlined. Last, we provide suggestions on how to overcome caveats and delineate some of the research avenues that will likely shape this field in the near future. We advocate that eDNA analyses, when carefully conducted and ideally combined with conventional sampling techniques, will substantially increase understanding and enable crucial expansion of subterranean community characterisation. Given the importance of groundwater and cave ecosystems for nature and humans, eDNA can bring to the surface essential insights, such as study of ecosystem assemblages and rare species detection, which are critical for the preservation of life below, as well as above, the ground.

Getting to the Root of Organic Inputs in Groundwaters: Stygofaunal Plant Consumption in a Calcrete Aquifer

Groundwater environments interact with and support subterranean biota as well as superficial aquatic and terrestrial ecosystems. However, knowledge of subterranean energy flows remains incomplete. Cross-boundary investigations are needed to better understand the trophic structures of groundwater ecosystems and their reliance on carbon inputs from aboveground. In this study we used carbon and nitrogen stable isotope analyses combined with radiocarbon fingerprints to characterise organic flows in groundwater ecosystems. We coupled these data with DNA metabarcoding of the gut contents of stygofauna to further elucidate organic matter (OM) sources and shifts in diet preferences. Samples were collected from the arid zone Sturt Meadows calcrete aquifer under low rainfall (LR) and high rainfall (HR) conditions. Bayesian modelling of Δ14C, δ13C, and δ15N data indicated that primary consumers (copepods) incorporated mainly particulate organic carbon (POC) under LR but during HR shifted to root derived material (either exudates or direct root grazing). By contrast, diets of secondary consumers (amphipods) were dominated by root material under both LR and HR. Our DNA metabarcoding-based results indicate that amphipods relied primarily on root inputs from perennial trees (likely Eucalyptus and Callitris) during the dry season (LR). Under HR, diets of both amphipods and copepods also included organic material derived from a broad range of more shallow rooted shrubs, and ephemeral herbs and grasses. Our findings illustrate the complexity of functional linkages between groundwater biota and surface terrestrial ecosystems in environments where aboveground productivity, diversity and OM flux to groundwater are intimately linked to often episodic rainfall.

Inferring predator-prey interaction in the subterranean environment: a case study from Dinaric caves

Predator–prey interactions are among the most important biotic interactions shaping ecological communities and driving the evolution of defensive traits. These interactions and their effects on species received little attention in extreme and remote environments, where possibilities for direct observations and experimental manipulation of the animals are limited. In this paper, we study such type of environment, namely caves of the Dinarides (Europe), combining spatial and phylogenetic methods. We focused on several species of Niphargus amphipods living in phreatic lakes, as some of them use the dorsal spines as putative morphological defensive traits. We predicted that these spines represent a defense strategy against the olm (Proteus anguinus), a top predator species in the subterranean waters. We tested for spatial overlap of the olm and Niphargus species and showed that spined species live in closer proximity to and co-occur more frequently with the olm than non-spined species. Modeling of the evolution of the spines onto Niphargus phylogeny implies coevolution of this trait in the presence of olm. We conclude that these spines likely evolved as defensive traits in a predator–prey arms race. Combining multiple analyses, we provide an example for a methodological framework to assess predator–prey interactions when in-situ or laboratory observations are not possible.

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.

Linking Hydrogeology and Ecology in Karst Landscapes: The Response of Epigean and Obligate Groundwater Copepods (Crustacea: Copepoda)

Groundwater invertebrate communities in karst landscapes are known to vary in response to multiple environmental factors. This study aims to explore the invertebrate assemblages’ composition of an Apennine karst system in Italy mainly described by the Rio Gamberale surface stream and the Stiffe Cave. The stream sinks into the carbonate rock and predominantly feeds the saturated karst into the cave. For a minor portion, groundwater flows from the epikarst and the perched aquifer within it. The spatial distribution of the species belonging to the selected target group of the Crustacea Copepoda between the surface stream and the groundwater habitats inside the cave highlighted a different response of surface-water species and obligate groundwater dwellers to the hydrogeological traits of the karst unit. Our results suggest that fast endorheic infiltration routes promoted the drift of epigean species from the surface to groundwater via the sinking stream while most of the obligate groundwater dwellers come from the perched aquifer in the epikarst from diffuse infiltration pathways.

Rainfall as a trigger of ecological cascade effects in an Australian groundwater ecosystem

Groundwaters host vital resources playing a key role in the near future. Subterranean fauna and microbes are crucial in regulating organic cycles in environments characterized by low energy and scarce carbon availability. However, our knowledge about the functioning of groundwater ecosystems is limited, despite being increasingly exposed to anthropic impacts and climate change-related processes. In this work we apply novel biochemical and genetic techniques to investigate the ecological dynamics of an Australian calcrete under two contrasting rainfall periods (LR—low rainfall and HR—high rainfall). Our results indicate that the microbial gut community of copepods and amphipods experienced a shift in taxonomic diversity and predicted organic functional metabolic pathways during HR. The HR regime triggered a cascade effect driven by microbes (OM processors) and exploited by copepods and amphipods (primary and secondary consumers), which was finally transferred to the aquatic beetles (top predators). Our findings highlight that rainfall triggers ecological shifts towards more deterministic dynamics, revealing a complex web of interactions in seemingly simple environmental settings. Here we show how a combined isotopic-molecular approach can untangle the mechanisms shaping a calcrete community. This design will help manage and preserve one of the most vital but underrated ecosystems worldwide.

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.

Refining trophic dynamics through multi-factor Bayesian mixing models: A case study of subterranean beetles

Food web dynamics are vital in shaping the functional ecology of ecosystems. However, trophic ecology is still in its infancy in groundwater ecosystems due to the cryptic nature of these environments. To unravel trophic interactions between subterranean biota, we applied an interdisciplinary Bayesian mixing model design (multi-factor BMM) based on the integration of faunal C and N bulk tissue stable isotope data (δ13C and δ15N) with radiocarbon data (Δ14C), and prior information from metagenomic analyses. We further compared outcomes from multi-factor BMM with a conventional isotope double proxy mixing model (SIA BMM), triple proxy (δ13C, δ15N, and Δ14C, multi-proxy BMM), and double proxy combined with DNA prior information (SIA + DNA BMM) designs. Three species of subterranean beetles (Paroster macrosturtensis, Paroster mesosturtensis, and Paroster microsturtensis) and their main prey items Chiltoniidae amphipods (AM1: Scutachiltonia axfordi and AM2: Yilgarniella sturtensis), cyclopoids and harpacticoids from a calcrete in Western Australia were targeted. Diet estimations from stable isotope only models (SIA BMM) indicated homogeneous patterns with modest preferences for amphipods as prey items. Multi-proxy BMM suggested increased-and species-specific-predatory pressures on amphipods coupled with high rates of scavenging/predation on sister species. SIA + DNA BMM showed marked preferences for amphipods AM1 and AM2, and reduced interspecific scavenging/predation on Paroster species. Multi-factorial BMM revealed the most precise estimations (lower overall SD and very marginal beetles' interspecific interactions), indicating consistent preferences for amphipods AM1 in all the beetles' diets. Incorporation of genetic priors allowed crucial refining of the feeding preferences, while integration of more expensive radiocarbon data as a third proxy (when combined with genetic data) produced more precise outcomes but close dietary reconstruction to that from SIA + DNA BMM. Further multidisciplinary modeling from other groundwater environments will help elucidate the potential behind these designs and bring light to the feeding ecology of one the most vital ecosystems worldwide.

Distribution patterns, carbon sources and niche partitioning in cave shrimps (Atyidae: Typhlatya)

Cave shrimps of the Typhlatya genus are common and widespread in fresh, brackish and marine groundwater throughout the Yucatan Peninsula (Mexico). These species are ideal models to test niche partitioning within sympatric species in oligotrophic systems. Nevertheless, their food sources remain unidentified, and despite their frequency and functional importance, distribution and abundance patterns of these species within caves have not been fully recognized. Here, we describe the abundance of three Typhlatya species in different temporal and spatial scales, investigate changes in water conditions, and potential sources of carbon as an indication of food origin. Species composition and abundance varied markedly in space and time revealing patterns that differed from one system to another and in relation to environmental parameters. Isotope analysis showed that each species reflects a particular δ¹³C and Δ¹⁴C fingerprint, suggesting they feed in different proportions from the available carbon sources. Overall, our findings suggest a niche partitioning of habitat and feeding sources amongst the three Typhlatya species investigated, where environmental characteristics and physiological differences could play an important role governing their distribution patterns.

Ecotoxicological effects of anthropogenic stressors in subterranean organisms: A review

How anthropogenic stressors affect biodiversity is a central question in a changing world. Subterranean ecosystems and their biodiversity are particularly vulnerable to change, yet, these species are frequently neglected in analyses of global biodiversity and assessments of ecological status and risk. Are these hidden species affected by anthropogenic stressors? Do they survive outside of the current thermal limits of their ecosystems? These and other important questions can be addressed with ecotoxicological testing, relating contaminants and temperature resistance of species with measured environmental concentrations and climatic data. Ecotoxicological knowledge specific to subterranean ecosystems is crucial for establishing thresholds for their protection, but such data are both scarce and scattered. Here, we review the existing ecotoxicological studies of these impacts to subterranean-adapted species. An effort that includes 167 measured endpoints and presents a database containing experimentally derived species' tolerance data for 28 contaminants and temperature, for 46 terrestrial and groundwater species, including fungi and animals. The lack of standard data among the studies is currently the major impediment to evaluate how stressors affect subterranean-adapted species and how differently they respond from their relatives at surface. Improving understanding of ecotoxicological effects on subterranean-adapted species will require extensive analysis of physiological responses to a wide range of untested stressors, standardization of testing protocols and evaluation of exposures under realistic scenarios.