Quantitative food web analysis supports the energy-limitation hypothesis in cave stream ecosystems

Comparative mitogenomic analysis of subterranean and surface amphipods (Crustacea, Amphipoda) with special reference to the family Crangonyctidae

Mitochondrial genomes play important roles in studying genome evolution, phylogenetic analyses, and species identification. Amphipods (Class Malacostraca, Order Amphipoda) are one of the most ecologically diverse crustacean groups occurring in a diverse array of aquatic and terrestrial environments globally, from freshwater streams and lakes to groundwater aquifers and the deep sea, but we have a limited understanding of how habitat influences the molecular evolution of mitochondrial energy metabolism. Subterranean amphipods likely experience different evolutionary pressures on energy management compared to surface-dwelling taxa that generally encounter higher levels of predation and energy resources and live in more variable environments. In this study, we compared the mitogenomes, including the 13 protein-coding genes involved in the oxidative phosphorylation (OXPHOS) pathway, of surface and subterranean amphipods to uncover potentially different molecular signals of energy metabolism between surface and subterranean environments in this diverse crustacean group. We compared base composition, codon usage, gene order rearrangement, conducted comparative mitogenomic and phylogenomic analyses, and examined evolutionary signals of 35 amphipod mitogenomes representing 13 families, with an emphasis on Crangonyctidae. Mitogenome size, AT content, GC-skew, gene order, uncommon start codons, location of putative control region (CR), length of rrnL and intergenic spacers differed between surface and subterranean amphipods. Among crangonyctid amphipods, the spring-dwelling Crangonyx forbesi exhibited a unique gene order, a long nad5 locus, longer rrnL and rrnS loci, and unconventional start codons. Evidence of directional selection was detected in several protein-encoding genes of the OXPHOS pathway in the mitogenomes of surface amphipods, while a signal of purifying selection was more prominent in subterranean species, which is consistent with the hypothesis that the mitogenome of surface-adapted species has evolved in response to a more energy demanding environment compared to subterranean amphipods. Overall, gene order, locations of non-coding regions, and base-substitution rates points to habitat as an important factor influencing the evolution of amphipod mitogenomes.

Environmental DNA of insects and springtails from caves reveals complex processes of eDNA transfer in soils

Subterranean environments host a substantial amount of biodiversity, however assessing the distribution of species living underground is still extremely challenging. Environmental DNA (eDNA) metabarcoding is a powerful tool to estimate biodiversity in poorly known environments and has excellent performance for soil organisms. Here, we tested 1) whether eDNA metabarcoding from cave soils/sediments allows to successfully detect springtails (Hexapoda: Collembola) and insects (Hexapoda: Insecta); 2) whether eDNA mostly represents autochthonous (cave-dwelling) organisms or it also incorporates information from species living in surface environments; 3) whether eDNA detection probability changes across taxa with different ecology. Environmental DNA metabarcoding analyses detected a large number of Molecular Operational Taxonomic Units (MOTUs) for both insects and springtails. For springtails, detection probability was high, with a substantial proportion of hypogean species, suggesting that eDNA provides good information on the distribution of these organisms in caves. Conversely, for insects most of MOTUs represented taxa living outside caves, and the majority of them represented taxa/organisms living in freshwater environments (Ephemeroptera, Plecoptera and Trichoptera). The eDNA of freshwater insects was particularly abundant in deep sectors of caves, far from the entrance. Furthermore, average detection probability of insects was significantly lower than the one of springtails. This suggests that cave soils/sediments act as "conveyer belts of biodiversity information", possibly because percolating water lead to the accumulation of eDNA of organisms living in nearby areas. Cave soils hold a complex mix of autochthonous and allochthonous eDNA. eDNA provided unprecedented information on the understudied subterranean cave organisms; analyses of detection probability and occupancy can help teasing apart local eDNA from the eDNA representing spatially-integrated biodiversity for whole landscape.

Longevity in Cave Animals

An extraordinary longevity has been observed in some cave species, and this raised the hypothesis that a longer lifespan may be considered one of the characteristic traits of these animals. However, only a few cave species have been studied thus far, and a firm conclusion remains to be drawn. Here we review the available knowledge on the longevity of subterranean species, point out the limitations of previous studies, and provide suggestions for future studies to answer important questions regarding the longevity in cave animals, its adaptive value and the related promoting factors. We also argue that studying the longevity in cave animals will contribute to the field of aging, especially to understanding the evolution of this phenomenon.

Postojna-Planina Cave System in Slovenia, a Hotspot of Subterranean Biodiversity and a Cradle of Speleobiology

The Postojna-Planina Cave System (PPCS) in central Slovenia is a globally exceptional site of subterranean biodiversity, comprised of many interconnected caves with cumulative passage length exceeding 34 km. Two rivers sink into the caves of the PPCS, called the Pivka and Rak, and join underground into Unica River, which emerges to the surface. The studies of fauna of PPCS began in the 19th century with the first scientific descriptions of specialized cave animals in the world, making it “the cradle of speleobiology”. Currently, the species list of PPCS contains 116 troglobiotic animal species belonging to eight phyla, confirming its status as the richest in the world. Of these, 47 species have been scientifically described from the PPCS, and more than 10 await formal taxonomic descriptions. We expect that further sampling, detailed analyses of less studied taxa, and the use of molecular methods may reveal more species. To keep the cave animals’ checklist in PPCS up-to-date, we have supplemented the printed checklist with an online interface. As the revised checklist is a necessary first step for further activities, we discuss the importance of PPCS in terms of future research and conservation.

Abiotic Community Constraints in Extreme Environments: Epikarst Copepods as a Model System

The general hypothesis that the overall presence or absence of one or more species in an extreme habitat is determined by physico-chemical factors was investigated using epikarst copepod communities as a model system, an example of an extreme environment with specialized, often rare species. The relationship between the presence or absence of epikarst copepods from drips in six Slovenian caves and 12 physico-chemical factors (temperature, conductivity, pH, Ca2+, Na+, K+, Mg2+, NH4+, and Cl−, NO2−, NO3−, and SO42−) was explored. Statistical analyses included principal components analysis, logistic mixed models, stepwise logistic multivariate regression, classification trees, and random forests. Parametric statistical analyses demonstrated the overall importance of two variables—temperature and conductivity. The more flexible statistical approaches, namely categorical trees and random forests, indicate that temperature and concentrations of Ca2+ and Mg2+ were important. This may be because they are essential nutrients or, at least in the case of Ca2+, its importance in molting. The correlation of Cl− and NO3− with copepod abundance may be due to unmeasured variables that vary at the scale of individual cave, but in any case, the values have an anthropogenic component. This contrasts with factors important in individual species’ niche separation, which overlap with the community parameters only for NO3−.

Cave-adapted evolution in the North American amblyopsid fishes inferred using phylogenomics and geometric morphometrics

Cave adaptation has evolved repeatedly across the Tree of Life, famously leading to pigmentation and eye degeneration and loss, yet its macroevolutionary implications remain poorly understood. We use the North American amblyopsid fishes, a family spanning a wide degree of cave adaptation, to examine the impact of cave specialization on the modes and tempo of evolution. We reconstruct evolutionary relationships using ultraconserved element loci, estimate the ancestral histories of eye-state, and examine the impact of cave adaptation on body shape evolution. Our phylogenomic analyses provide a well-supported hypothesis for amblyopsid evolutionary relationships. The obligate blind cavefishes form a clade and the cave-facultative eyed spring cavefishes are nested within the obligate cavefishes. Using ancestral state reconstruction, we find support for at least two independent subterranean colonization events within the Amblyopsidae. Eyed and blind fishes have different body shapes, but not different rates of body shape evolution. North American amblyopsids highlight the complex nature of cave-adaptive evolution and the necessity to include multiple lines of evidence to uncover the underlying processes involved in the loss of complex traits.

Ecophysiological and life-history adaptations of Gammarus balcanicus (Schäferna, 1922) in a sinking-cave stream from Western Carpathians (Romania)

Freshwater gammarids are known to colonise occasionally sinking-cave streams, providing contrasting morphological, life-history and ecophysiological adaptations compared to their surface conspecifics. In this study, a subterranean and a surface population of the species Gammarus balcanicus was surveyed for one year in a sinking-cave stream from the Western Carpathians (Romania). The results showed that the cave-dwelling population comprised individuals that were significantly larger compared to their surface conspecifics, had larger body-size at sexual maturity and that the females produced fewer, but larger eggs, compared to the population situated outside the cave. The trophic position and the omnivory were significantly higher for the cave-dwelling compared to surface population and the elemental imbalance for C:P molar ratios lower, but similar for C:N. However, the subterranean population did not present troglomorphic characters or longer lifespan as known for other cave-surface paired crustaceans. This, together with the rather extensive hydrological connection of the habitats, suggests active gene-flow between populations and similar response to seasonality for body-size distributions, indicating that the observed ecophysiological and life-history differences are rather the consequence of phenotypic plasticity than the result of genetic adaptation.

Metabolic responses to long-term food deprivation in subterranean and surface amphipods

A long-standing hypothesis in subterranean biology posits that organisms living in poor resource subsurface habitats can withstand long periods of bioenergetic shortages due to an innate reduced metabolic rate when compared to their epigean counterparts. However, previous studies have proposed that caves with ample energy resources may not evolve organisms with reduced metabolic rate. The equivocal nature of previous findings suggests that there is a need to compare food deprivation responses of subterranean and surface species in order to elucidate whether there are widespread adaptations to low energy systems in subterranean taxa. The purpose of the study was to examine patterns in basal metabolism and the effects of food deprivation in closely related subterranean- and epigean- amphipods,Stygobromus peckiandSynurellasp. from central and east Texas, USA, respectively. Basal metabolic rates (measured as O2consumption) differed between species, withS. peckihaving substantially lower rates thanSynurella. Individuals of both species were food deprived for a pre-determined time interval and changes in total body protein, lipids, and carbohydrates were measured throughout food deprivation experiments.Stygobromus peckihad larger initial energy stores thanSynurellaand were more conservative in the use of energetic reserves over a prolonged period of food deprivation. Thus, it appears that althoughS. peckiare currently found in shallow phreatic and spring opening environments, they have maintained more efficient metabolic adaptations to deal with prolonged periods of food deprivation.

Predation risk and microhabitat selection by cave salamanders, Eurycea lucifuga (Rafinesque, 1822)

Habitat selection is driven by many factors, but no one location is likely to be best for all factors; thus, individuals are subject to trade-offs when selecting habitat. Caves provide a clear example of such trade-offs because these habitats are energy deprived. Cave salamanders (Eurycea lucifuga) commonly inhabit caves at least in part because this habitat is cool and wet. We tested the hypothesis that caves also provide cave salamanders with a reduction in predation risk. We used clay models to test for differences in predation risk in caves vs. forests and at low (e.g., ground) vs. elevated (e.g., cave wall) positions, and recorded locations of cave salamanders to assess vertical (i.e., wall) vs. non-vertical (e.g., ground) substrate selection in a cave. Overall, a mean of 3.2 models in caves were damaged and a mean of 8.2 were damaged in forests. Cave salamanders selected vertical substrate more often than non-vertical substrate (, ), and in caves, low-positioned models were more likely to be damaged than models on walls, although there was no effect of model height in forests. This study suggests that caves provide salamanders a refuge from increased predation pressure, a benefit that likely compensates for the costs of moving between caves and the richer foraging grounds outside them.

Less effective selection leads to larger genomes

The evolutionary origin of the striking genome size variations found in eukaryotes remains enigmatic. The effective size of populations, by controlling selection efficacy, is expected to be a key parameter underlying genome size evolution. However, this hypothesis has proved difficult to investigate using empirical data sets. Here, we tested this hypothesis using 22 de novo transcriptomes and low-coverage genomes of asellid isopods, which represent 11 independent habitat shifts from surface water to resource-poor groundwater. We show that these habitat shifts are associated with higher transcriptome-wide [Formula: see text] After ruling out the role of positive selection and pseudogenization, we show that these transcriptome-wide [Formula: see text] increases are the consequence of a reduction in selection efficacy imposed by the smaller effective population size of subterranean species. This reduction is paralleled by an important increase in genome size (25% increase on average), an increase also confirmed in subterranean decapods and mollusks. We also control for an adaptive impact of genome size on life history traits but find no correlation between body size, or growth rate, and genome size. We show instead that the independent increases in genome size measured in subterranean isopods are the direct consequence of increasing invasion rates by repeat elements, which are less efficiently purged out by purifying selection. Contrary to selection efficacy, polymorphism is not correlated to genome size. We propose that recent demographic fluctuations and the difficulty of observing polymorphism variation in polymorphism-poor species can obfuscate the link between effective population size and genome size when polymorphism data are used alone.

No Evidence That Nitrogen Limitation Influences the Elemental Composition of Isopod Transcriptomes and Proteomes

The field of stoichiogenomics aims at understanding the influence of nutrient limitations on the elemental composition of the genome, transcriptome, and proteome. The 20 amino acids and the 4 nt differ in the number of nutrients they contain, such as nitrogen (N). Thus, N limitation shall theoretically select for changes in the composition of proteins or RNAs through preferential use of N-poor amino acids or nucleotides, which will decrease the N-budget of an organism. While these N-saving mechanisms have been evidenced in microorganisms, they remain controversial in multicellular eukaryotes. In this study, we used 13 surface and subterranean isopod species pairs that face strongly contrasted N limitations, either in terms of quantity or quality. We combined in situ nutrient quantification and transcriptome sequencing to test if N limitation selected for N-savings through changes in the expression and composition of the transcriptome and proteome. No evidence of N-savings was found in the total N-budget of transcriptomes or proteomes or in the average protein N-cost. Nevertheless, subterranean species evolving in N-depleted habitats displayed lower N-usage at their third codon positions. To test if this convergent compositional change was driven by natural selection, we developed a method to detect the strand-asymmetric signature that stoichiogenomic selection should leave in the substitution pattern. No such signature was evidenced, indicating that the observed stoichiogenomic-like patterns were attributable to nonadaptive processes. The absence of stoichiogenomic signal despite strong N limitation within a powerful phylogenetic framework casts doubt on the existence of stoichiogenomic mechanisms in metazoans.

More than a corridor: use of a main stem stream as supplemental foraging habitat by a brook trout metapopulation

Coldwater fishes in streams, such as brook trout (Salvelinus fontinalis), typically are headwater specialists that occasionally expand distributions downstream to larger water bodies. It is unclear, however, whether larger streams function simply as dispersal corridors connecting headwater subpopulations, or as critical foraging habitat needed to sustain large mobile brook trout. Stable isotopes (δ(13)C and δ(15)N) and a hierarchical Bayesian mixing model analysis was used to identify brook trout that foraged in main stem versus headwater streams of the Shavers Fork watershed, West Virginia. Headwater subpopulations were composed of headwater and to a lesser extent main stem foraging individuals. However, there was a strong relationship between brook trout size and main stem prey contributions. The average brook trout foraging on headwater prey were limited to 126 mm standard length. This size was identified by mixing models as a point where productivity support switched from headwater to main stem dependency. These results, similar to other studies conducted in this watershed, support the hypothesis that productive main stem habitat maintain large brook trout and potentially facilitates dispersal among headwater subpopulations. Consequently, loss of supplementary main stem foraging habitats may explain loss of large, mobile fish and subsequent isolation of headwater subpopulations in other central Appalachian watersheds.