Freshwater connectivity transforms spatially integrated signals of biodiversity.

Aquatic ecosystems offer a continuum of water flow from headwater streams to inland lakes and coastal marine systems. This spatial connectivity influences the structure, function and dynamics of aquatic communities, which are among the most threatened and degraded on the Earth. Here, we determine the spatial resolution of environmental DNA (eDNA) in dendritic freshwater networks, which we use as a model for connected metacommunities. Our intensive sampling campaign comprised over 420 eDNA samples across 21 connected lakes, allowing us to analyse detections at a variety of scales, from different habitats within a lake to entire lake networks. We found strong signals of within-lake variation in eDNA distribution reflective of typical habitat use by both fish and zooplankton. Most importantly, we also found that connecting channels between lakes resulted in an accumulation of downstream eDNA detections in lakes with a higher number of inflows, and as networks increased in length. Environmental DNA achieves biodiversity surveys in these habitats in a high-throughput, spatially integrated way. These findings have profound implications for the interpretation of eDNA detections in aquatic ecosystems in global-scale biodiversity monitoring observations.

Environmental transcriptomics under heat stress: Can environmental RNA reveal changes in gene expression of aquatic organisms?

To safeguard biodiversity in a changing climate, taxonomic information about species turnover and insights into the health of organisms are required. Environmental DNA approaches are increasingly used for species identification, but cannot provide functional insights. Transcriptomic methods reveal the physiological states of macroorganisms, but are currently species-specific and require tissue sampling or animal sacrifice, making community-wide assessments challenging. Here, we test whether broad functional information (expression level of the transcribed genes) can be harnessed from environmental RNA (eRNA), which includes extra-organismal RNA from macroorganisms along with whole microorganisms. We exposed Daphnia pulex as well as phytoplankton prey and microorganism colonizers to control (20°C) and heat stress (28°C) conditions for 7 days. We sequenced eRNA from tank water (after complete removal of Daphnia) as well as RNA from Daphnia tissue, enabling comparisons of extra-organismal and organismal RNA-based gene expression profiles. Both RNA types detected similar heat stress responses of Daphnia. Using eRNA, we identified 32 Daphnia genes to be differentially expressed following heat stress. Of these, 17 were also differentially expressed and exhibited similar levels of relative expression in organismal RNA. In addition to the extra-organismal Daphnia response, eRNA detected community-wide heat stress responses consisting of distinct functional profiles and 121 differentially expressed genes across eight taxa. Our study demonstrates that environmental transcriptomics based on extra-organismal eRNA can noninvasively reveal gene expression responses of macroorganisms following environmental changes, with broad potential implications for the biomonitoring of health across the trophic chain.

Accelerated rates of large-scale mutations in the presence of copper and nickel.

Mutation rate variation has been under intense investigation for decades. Despite these efforts, little is known about the extent to which environmental stressors accelerate mutation rates and influence the genetic load of populations. Moreover, most studies on stressors have focused on unicellular organisms and point mutations rather than large-scale deletions and duplications (copy number variations [CNVs]). We estimated mutation rates in Daphnia pulex exposed to low levels of environmental stressors as well as the effect of selection on de novo mutations. We conducted a mutation accumulation (MA) experiment in which selection was minimized, coupled with an experiment in which a population was propagated under competitive conditions in a benign environment. After an average of 103 generations of MA propagation, we sequenced 60 genomes and found significantly accelerated rates of deletions and duplications in MA lines exposed to ecologically relevant concentrations of metals. Whereas control lines had gene deletion and duplication rates comparable to other multicellular eukaryotes (1.8 × 10-6 per gene per generation), the presence of nickel and copper increased these rates fourfold. The realized mutation rate under selection was reduced to 0.4× that of control MA lines, providing evidence that CNVs contribute to mutational load. Our CNV breakpoint analysis revealed that nonhomologous recombination associated with regions of DNA fragility is the primary source of CNVs, plausibly linking metal-induced DNA strand breaks with higher CNV rates. Our findings suggest that environmental stress, in particular multiple stressors, can have profound effects on large-scale mutation rates and mutational load of multicellular organisms.

Genotype diversity promotes the persistence of Daphnia populations exposed to severe copper stress.

When environmental stressors of high intensity are sustained for long periods of time, populations face high probabilities of being extirpated. However, depending on the intensity of the stressor, large populations with sufficient genetic diversity may persist. We report the results of an experiment that tracked the persistence of Daphnia populations exposed to copper contamination. We assessed whether genotypic diversity reduced the risk of extinction. We created monoclonal and multiclonal populations and monitored their population sizes during a 32-week experiment. Cu was applied at a sub-lethal concentration and then increased every week until the population sizes dropped to about 10% of the carrying capacity (Cu at 180 µg/L). The concentration was then increased up to 186 µg/L and held stable until the end of the experiment. A survival analysis showed that clonal diversity extended the persistence of Daphnia populations, but copper contamination caused a substantial genetic erosion followed by population extirpation. However, some Cu-treated populations, mostly multiclonal, showed U-shaped patterns of growth consistent with evolutionary rescue but these did not lead to lasting population recovery. These results highlight the importance of genetic variation for population persistence, but they also show how quickly it can be lost in contaminated environments.

Landscape modification and nutrient-driven instability at a distance.

Almost 50 years ago, Michael Rosenzweig pointed out that nutrient addition can destabilise food webs, leading to loss of species and reduced ecosystem function through the paradox of enrichment. Around the same time, David Tilman demonstrated that increased nutrient loading would also be expected to cause competitive exclusion leading to deleterious changes in food web diversity. While both concepts have greatly illuminated general diversity-stability theory, we currently lack a coherent framework to predict how nutrients influence food web stability across a landscape. This is a vitally important gap in our understanding, given mounting evidence of serious ecological disruption arising from anthropogenic displacement of resources and organisms. Here, we combine contemporary theory on food webs and meta-ecosystems to show that nutrient additions are indeed expected to drive loss in stability and function in human-impacted regions. Our models suggest that destabilisation is more likely to be caused by the complete loss of an equilibrium due to edible plant species being competitively excluded. In highly modified landscapes, spatial nutrient transport theory suggests that such instabilities can be amplified over vast distances from the sites of nutrient addition. Consistent with this theoretical synthesis, the empirical frequency of these distant propagating ecosystem imbalances appears to be growing. This synthesis of theory and empirical data suggests that human modification of the Earth is strongly connecting distantly separated ecosystems, causing rapid, expansive and costly nutrient-driven instabilities over vast areas of the planet. Similar to existing food web theory, the corollary to this spatial nutrient theory is that slowing down spatial nutrient pathways can be a potent means of stabilising degraded ecosystems.

Speciation in Daphnia.

The microcrustacean Daphnia is arguably one of the most studied zooplankton species, having a well understood ecology, life history, and a relatively well studied evolutionary history. Despite this wealth of knowledge, species boundaries within closely related species in this genus often remain elusive and the major evolutionary forces driving the diversity of daphniids remain controversial. This genus contains more than 80 species with multiple cryptic species complexes, with many closely related species able to hybridize. Here, we review speciation research in Daphnia within the framework of current speciation theory. We evaluate the role of geography, ecology, and biology in restricting gene flow and promoting diversification. Of the 253 speciation studies on Daphnia, the majority of studies examine geographic barriers (55%). While evidence shows that geographic barriers play a role in species divergence, ecological barriers are also probably prominent in Daphnia speciation. We assess the contribution of ecological and nonecological reproductive isolating barriers between closely related species of Daphnia and found that none of the reproductive isolating barriers are restricting gene flow completely. Research on reproductive isolating barriers has disproportionally focused on two species complexes, the Daphnia pulex and Daphnia longispina species complexes. Finally, we identify areas of research that remain relatively unexplored and discuss future research directions that build our understanding of speciation in daphniids.

Integrating physiology and environmental dynamics to operationalize environmental DNA (eDNA) as a means to monitor freshwater macro-organism abundance.

Research has demonstrated consistent positive correlations between organism abundance and absolute environmental DNA (eDNA) concentrations. Robust correlations in laboratory experiments indicate strong functional links, suggesting the potential for eDNA to monitor organism abundance in nature. However, correlations between absolute eDNA concentrations and organism abundance in nature tend to be weaker because myriad biotic and abiotic factors influence steady-state eDNA concentrations, decoupling its direct functional link with abundance. Additional technical challenges can also weaken correlations between relative organism abundance and relative eDNA data derived from metabarcoding. Future research must account for these factors to improve the inference of organism abundance from eDNA, including integrating the effects of organism physiology on eDNA production, eDNA dynamics in lentic/lotic systems, and key environmental parameters that impact estimated steady-state concentrations. Additionally, it is critical to manage expectations surrounding the accuracy and precision that eDNA can provide - eDNA, for example, cannot provide abundance estimates comparable to intensively managed freshwater fisheries that enumerate every individual fish. Recent developments, however, are encouraging. Current methods could provide meaningful information regarding qualitative conservation thresholds and emergent research has demonstrated that eDNA concentrations in natural ecosystems can provide rough quantitative estimates of abundance, particularly when models integrate physiology and/or eDNA dynamics. Operationalizing eDNA to infer abundance will probably require more than simple correlations with organism biomass/density. Nevertheless, the future is promising - models that integrate eDNA dynamics in nature could represent an effective means to infer abundance, particularly when traditional methods are considered too "costly" or difficult to obtain.

Thermal stratification and fish thermal preference explain vertical eDNA distributions in lakes.

Significant advances have been made towards surveying animal and plant communities using DNA isolated from environmental samples. Despite rapid progress, we lack a comprehensive understanding of the "ecology" of environmental DNA (eDNA), particularly its temporal and spatial distribution and how this is shaped by abiotic and biotic processes. Here, we tested how seasonal variation in thermal stratification and animal habitat preferences influences the distribution of eDNA in lakes. We sampled eDNA depth profiles of five dimictic lakes during both summer stratification and autumn turnover, each containing warm- and cool-water fishes as well as the cold-water stenotherm, lake trout (Salvelinus namaycush). Habitat use by S. namaycush was validated by acoustic telemetry and was significantly related to eDNA distribution during stratification. Fish eDNA became "stratified" into layers during summer months, reflecting lake stratification and the thermal niches of the species. During summer months, S. namaycush, which rarely ventured into shallow waters, could only be detected at the deepest layers of the lakes, whereas the eDNA of warm-water fishes was much more abundant above the thermocline. By contrast, during autumn lake turnover, the fish species assemblage as detected by eDNA was homogenous throughout the water column. These findings contribute to our overall understanding of the "ecology" of eDNA within lake ecosystems, illustrating how the strong interaction between seasonal thermal structure in lakes and thermal niches of species on very localized spatial scales influences our ability to detect species.

Alternative splicing is highly variable among Daphnia pulex lineages in response to acute copper exposure.

BACKGROUND: Despite being one of the primary mechanisms of gene expression regulation in eukaryotes, alternative splicing is often overlooked in ecotoxicogenomic studies. The process of alternative splicing facilitates the production of multiple mRNA isoforms from a single gene thereby greatly increasing the diversity of the transcriptome and proteome. This process can be important in enabling the organism to cope with stressful conditions. Accurate identification of splice sites using RNA sequencing requires alignment to independent exonic positions within the genome, presenting bioinformatic challenges, particularly when using short read data. Although technological advances allow for the detection of splicing patterns on a genome-wide scale, very little is known about the extent of intraspecies variation in splicing patterns, particularly in response to environmental stressors. In this study, we used RNA-sequencing to study the molecular responses to acute copper exposure in three lineages of Daphnia pulex by focusing on the contribution of alternative splicing in addition to gene expression responses. RESULTS: By comparing the overall gene expression and splicing patterns among all 15 copper-exposed samples and 6 controls, we identified 588 differentially expressed (DE) genes and 16 differentially spliced (DS) genes. Most of the DS genes (13) were not found to be DE, suggesting unique transcriptional regulation in response to copper that went unnoticed with conventional DE analysis. To understand the influence of genetic background on gene expression and alternative splicing responses to Cu, each of the three lineages was analyzed separately. In contrast to the overall analysis, each lineage had a higher proportion of unique DS genes than DE genes suggesting that genetic background has a larger influence on DS than on DE. Gene Ontology analysis revealed that some pathways involved in stress response were jointly regulated by DS and DE genes while others were regulated by only transcription or only splicing. CONCLUSIONS: Our findings suggest an important role for alternative splicing in shaping transcriptome diversity in response to metal exposure in Daphnia, highlighting the importance of integrating splicing analyses with gene expression surveys to characterize molecular pathways in evolutionary and environmental studies.

The evolution of reproductive isolation in Daphnia.

BACKGROUND: The process by which populations evolve to become new species involves the emergence of various reproductive isolating barriers (RIB). Despite major advancements in understanding this complex process, very little is known about the order in which RIBs evolve or their relative contribution to the total restriction of gene flow during various stages of speciation. This is mainly due to the difficulties of studying reproductive isolation during the early stages of species formation. This study examines ecological and non-ecological RIB within and between Daphnia pulex and Daphnia pulicaria, two recently diverged species that inhabit distinct habitats and exhibit an unusual level of intraspecific genetic subdivision. RESULTS: We find that while ecological prezygotic barriers are close to completion, none of the non-ecological barriers can restrict gene flow between D. pulex and D. pulicaria completely when acting alone. Surprisingly, we also identified high levels of postzygotic reproductive isolation in 'conspecific' interpopulation crosses of D. pulex. CONCLUSIONS: While the ecological prezygotic barriers are prevalent during the mature stages of speciation, non-ecological barriers likely dominated the early stages of speciation. This finding indicates the importance of studying the very early stages of speciation and suggests the contribution of postzygotic isolation in initiating the process of speciation.

Trophic structure modulates community rescue following acidification.

Community rescue occurs when a community that experiences lethal stress persists only through the spread of rare types, either genotypes or species, resistant to the stress. Rescue interacts with trophic structure because physical stress experienced by a focal assemblage within the community may also be experienced by its predators and prey. In general, trophic structure will facilitate rescue only when a stress has a less severe effect on a focal assemblage than on its predators. In other circumstances, when stress affects prey or has only a weak effect on predators, trophic structure is likely to hamper rescue. We exposed a community of phytoplankton and zooplankton derived from a natural lake to acidification in outdoor mesocosms large enough to support trophically complex communities. Rescue of the phytoplankton from severe acidification was facilitated by prior exposure to sublethal stress, confirming previous results from microcosm experiments. Even communities that have previously been less highly stressed were eventually rescued, however, because their zooplankton predators were more sensitive to acidification and became extinct. Our experiment shows how community rescue following severe stress is modulated by the differential effect of the stress relative to trophic level.

Proteomic Profile of Daphnia pulex using Data-Independent Acquisition Mass Spectrometry and Ion Mobility Separation.

Daphnia pulex is a keystone species for aquatic habitats and an ecological/evolution model organism. Although significant progress has been made on characterizing its genome, the D. pulex proteome remains largely uncharacterized partially due to abnormally high protein degradation during homogenization and emphasis on genomic analysis. In this study, various sample preparation and mass spectrometry acquisition methods are performed for the purpose of improving D. pulex proteome exploration. Benefits for employing both in-gel and in-solution methods of trypsin digestion are observed. Furthermore, acquisition methods employing ion mobility separation greatly increase peptide identification and more than doubled the proteome coverage. Bioinformatic analysis suggests that mitochondrial and hydrolytic activities are enriched in D. pulex compared to closely related invertebrates or Homo sapiens. Also, novel D. pulex proteins possessing putative genome modifying functional domains are identified. Data are available via ProteomeXchange with identifier PXD008455.

Spontaneous Mutation Accumulation in Daphnia pulex in Selection-Free vs. Competitive Environments.

Understanding the rates, spectra, and fitness effects of spontaneous mutations is fundamental to answering key questions in evolution, molecular biology, disease genetics, and conservation biology. To estimate mutation rates and evaluate the effect of selection on new mutations, we propagated mutation accumulation (MA) lines of Daphnia pulex for more than 82 generations and maintained a non-MA population under conditions where selection could act. Both experiments were started with the same obligate asexual progenitor clone. By sequencing 30 genomes and implementing a series of validation steps that informed the bioinformatic analyses, we identified a total of 477 single nucleotide mutations (SNMs) in the MA lines, corresponding to a mutation rate of 2.30 × 10-9 (95% CI 1.90-2.70 × 10-9) per nucleotide per generation. The high overall rate of loss of heterozygosity (LOH) of 4.82 × 10-5 per site per generation was due to a large ameiotic recombination event spanning an entire arm of a chromosome (∼6 Mb) and several hemizygous deletion events spanning ∼2 kb each. In the non-MA population, we found significantly fewer mutations than expected based on the rate derived from the MA experiment, indicating purifying selection was likely acting to remove new deleterious mutations. We observed a surprisingly high level of genetic variability in the non-MA population, which we propose to be driven by balancing selection. Our findings suggest that both positive and negative selection on new mutations is powerful and effective in a strictly clonal population.

The effects of parameter choice on defining molecular operational taxonomic units and resulting ecological analyses of metabarcoding data.

The combination of DNA barcoding and high-throughput (next-generation) sequencing (metabarcoding) provides many promises but also serious challenges. Generating a reliable comparable estimate of biodiversity remains a central challenge to the application of the technology. Many approaches have been used to turn millions of sequences into distinct taxonomic units. However, the extent to which these methods impact the outcome of simple ecological analyses is not well understood. Here we performed a simple analysis of dietary overlap by skinks and shrews on Ile Aux Aigrettes, Mauritius. We used a combination of filtering thresholds and clustering algorithms on a COI metabarcoding dataset and demonstrate that all bioinformatics parameters will have interacting effects on molecular operational taxonomic unit (MOTU) recovery rates. These effects generated estimates covering two orders of magnitude. However, the effect on a simple ecological analysis was not large and, despite the wide variation in estimates of niche overlap, the same ecological conclusion was drawn in most cases. We advise that a conservative clustering programme coupled with larger sequence divergences to define a cluster, the removal of singletons, rigorous length filtering, and stringent match criteria for Molecular Identifier tags are preferable to avoid MOTU inflation and that the same parameters be used in all comparative analyses.

Genetic reconstructions of invasion history.

A diverse array of molecular markers and constantly evolving analytical approaches have been employed to reconstruct the invasion histories of the most notorious invasions. Detailed information on the source(s) of introduction, invasion route, type of vectors, number of independent introductions and pathways of secondary spread has been corroborated for a large number of biological invasions. In this review, I present the promises and limitations of current techniques while discussing future directions. Broad phylogeographic surveys of native and introduced populations have traced back invasion routes with surprising precision. These approaches often further clarify species boundaries and reveal complex patterns of genetic relationships with noninvasive relatives. Moreover, fine-scale analyses of population genetics or genomics allow deep inferences on the colonization dynamics across invaded ranges and can reveal the extent of gene flow among populations across various geographical scales, major demographic events such as genetic bottlenecks as well as other important evolutionary events such as hybridization with native taxa, inbreeding and selective sweeps. Genetic data have been often corroborated successfully with historical, geographical and ecological data to enable a comprehensive reconstruction of the invasion process. The advent of next-generation sequencing, along with the availability of extensive databases of repository sequences generated by barcoding projects opens the opportunity to broadly monitor biodiversity, to identify early invasions and to quantify failed invasions that would otherwise remain inconspicuous to the human eye.

Synergistic interactions of biotic and abiotic environmental stressors on gene expression.

Understanding the response of organisms to multiple stressors is critical for predicting if populations can adapt to rapid environmental change. Natural and anthropogenic stressors often interact, complicating general predictions. In this study, we examined the interactive and cumulative effects of two common environmental stressors, lowered calcium concentration, an anthropogenic stressor, and predator presence, a natural stressor, on the water flea Daphnia pulex. We analyzed expression changes of five genes involved in calcium homeostasis - cuticle proteins (Cutie, Icp2), calbindin (Calb), and calcium pump and channel (Serca and Ip3R) - using real-time quantitative PCR (RT-qPCR) in a full factorial experiment. We observed strong synergistic interactions between low calcium concentration and predator presence. While the Ip3R gene was not affected by the stressors, the other four genes were affected in their transcriptional levels by the combination of the stressors. Transcriptional patterns of genes that code for cuticle proteins (Cutie and Icp2) and a sarcoplasmic calcium pump (Serca) only responded to the combination of stressors, changing their relative expression levels in a synergistic response, while a calcium-binding protein (Calb) responded to low calcium stress and the combination of both stressors. The expression pattern of these genes (Cutie, Icp2, and Serca) were nonlinear, yet they were dose dependent across the calcium gradient. Multiple stressors can have complex, often unexpected effects on ecosystems. This study demonstrates that the dominant interaction for the set of tested genes appears to be synergism. We argue that gene expression patterns can be used to understand and predict the type of interaction expected when organisms are exposed simultaneously to natural and anthropogenic stressors.

Revealing population demographics with environmental RNA.

The analyses of environmental DNA (eDNA) and environmental RNA (eRNA) released by organisms into their surrounding environment (water, soil and air) have emerged as powerful tools for monitoring biodiversity. While eDNA has been widely adopted for the non-invasive detection of species and characterization of community composition, the utilization of eRNA is still in its infancy. Due to its functional nature, eRNA holds intriguing potential for biodiversity monitoring offering new avenues of research beyond species detection. For example, conspecifics that are almost genetically identical can exhibit distinct transcriptomic differences depending on their life stage. In this issue of Molecular Ecology Resources, Parsley and Goldberg (2024) demonstrate, through a lab-validated field study, that eRNA can be used to detect distinct life stages of amphibians. This study elegantly demonstrates that eRNA can be used not only to detect invasive or endangered species but also to reveal population demographic information important for guiding effective conservation strategies.