Search

Choosing an aquatic invertebrate sampling method for biomonitoring depends upon study goals, resources, and ecosystem conditions. In this study, we compared two methods that are widely used in stream ecology, but have not been directly compared: Hester–Dendy (HD) and Hess sampling. Hester–Dendy sampling uses artificial substrate that invertebrates colonize over a specific period of time. In contrast, Hess samplers surround a fixed area of natural substrate with a net. To compare approaches, we combined 5 years of simultaneous HD and Hess data collection (2010–2014) from the Niobrara River with a 14-year (1996–2009) historical HD data set for the same study sites. We used this full 19-year data set to assess how ecosystem health has changed in the Niobrara River over time, while also testing the influence of HD versus Hess data (2010–2014) on historical trends (1996–2009). Our results showed that HD samples are more taxonomically variable and bias bioassessment metrics because they collect more sensitive taxa versus Hess sampling. However, when combined with the 1996–2009 HD data set, both recent HD and Hess data sets recovered the same trend of declining ecosystem health in the Niobrara River. These results provide empirical evidence that even when historical HD data are combined with recent Hess data, long-term bioassessment trends remain unchanged despite more accurate perspectives of invertebrate assemblages being collected.

Pteronarcys californica (Newport 1848) is commonly referred to as the giant salmonfly and is the largest species of stonefly (Insecta: Plecoptera) in the western United States. Historically, it was widespread and abundant in western rivers, but populations have experienced a substantial decline in the past few decades, becoming locally extirpated in numerous rivers in Utah, Colorado, and Montana. Although previous research has explored the ecological variables conducive to the survivability of populations of the giant salmonfly, a lack of genomic resources hampers exploration of how genetic variation is spread across extant populations. To accelerate research on this imperiled species, we present a de novo chromosomal-length genome assembly of P. californica generated from PacBio HiFi sequencing and Hi-C chromosome conformation capture. Our assembly includes 14 predicted pseudo chromosomes and 98.8% of Insecta universal core orthologs. At 2.40 gigabases, the P. californica assembly is the largest of available stonefly assemblies, highlighting at least 9.5-fold variation in assembly size across the order. Repetitive elements (REs) account for much of the genome size increase in P. californica relative to other stonefly species, with the content of Class I retroelements alone exceeding the entire assembly size of all but two other species studied. We also observed preliminary suborder-specific trends in genome size that merit testing with more robust taxon sampling., (© The American Genetic Association. 2024.)

Caddisflies (Trichoptera) are among the most diverse groups of freshwater animals with more than 16 000 described species. They play a fundamental role in freshwater ecology and environmental engineering in streams, rivers and lakes. Because of this, they are frequently used as indicator organisms in biomonitoring programmes. Despite their importance, key questions concerning the evolutionary history of caddisflies, such as the timing and origin of larval case making, remain unanswered owing to the lack of a well-resolved phylogeny. Here, we estimated a phylogenetic tree using a combination of transcriptomes and targeted enrichment data for 207 species, representing 48 of 52 extant families and 174 genera. We calibrated and dated the tree with 33 carefully selected fossils. The first caddisflies originated approximately 295 million years ago in the Permian, and major suborders began to diversify in the Triassic. Furthermore, we show that portable case making evolved in three separate lineages, and shifts in diversification occurred in concert with key evolutionary innovations beyond case making.

We present a genome assembly from one male Meconema thalassinum (the oak bush-cricket; Arthropoda; Insecta; Orthoptera; Tettigoniidae). The genome sequence is 9,039.1 megabases in span. Most of the assembly is scaffolded into 15 chromosomal pseudomolecules, including the X sex chromosome. The mitochondrial genome has also been assembled and is 15.63 kilobases in length., Competing Interests: No competing interests were disclosed., (Copyright: © 2024 Crowley LM et al.)

Repetitive elements (REs) are integral to the composition, structure, and function of eukaryotic genomes, yet remain understudied in most taxonomic groups. We investigated REs across 601 insect species and report wide variation in RE dynamics across groups. Analysis of associations between REs and protein-coding genes revealed dynamic evolution at the interface between REs and coding regions across insects, including notably elevated RE-gene associations in lineages with abundant long interspersed nuclear elements (LINEs). We leveraged this large, empirical data set to quantify impacts of long-read technology on RE detection and investigate fundamental challenges to RE annotation in diverse groups. In long-read assemblies, we detected ∼36% more REs than short-read assemblies, with long terminal repeats (LTRs) showing 162% increased detection, whereas DNA transposons and LINEs showed less respective technology-related bias. In most insect lineages, 25%-85% of repetitive sequences were "unclassified" following automated annotation, compared with only ∼13% in Drosophila species. Although the diversity of available insect genomes has rapidly expanded, we show the rate of community contributions to RE databases has not kept pace, preventing efficient annotation and high-resolution study of REs in most groups. We highlight the tremendous opportunity and need for the biodiversity genomics field to embrace REs and suggest collective steps for making progress toward this goal., (© 2023 Sproul et al.; Published by Cold Spring Harbor Laboratory Press.)

Arthropod silk is vital to the evolutionary success of hundreds of thousands of species. The primary proteins in silks are often encoded by long, repetitive gene sequences. Until recently, sequencing and assembling these complex gene sequences has proven intractable given their repetitive structure. Here, using high-quality long-read sequencing, we show that there is extensive variation-both in terms of length and repeat motif order-between alleles of silk genes within individual arthropods. Further, this variation exists across two deep, independent origins of silk which diverged more than 500 Mya: the insect clade containing caddisflies and butterflies and spiders. This remarkable convergence in previously overlooked patterns of allelic variation across multiple origins of silk suggests common mechanisms for the generation and maintenance of structural protein-coding genes. Future genomic efforts to connect genotypes to phenotypes should account for such allelic variation.

Background: Generating the most contiguous, accurate genome assemblies given available sequencing technologies is a long-standing challenge in genome science. With the rise of long-read sequencing, assembly challenges have shifted from merely increasing contiguity to correctly assembling complex, repetitive regions of interest, ideally in a phased manner. At present, researchers largely choose between two types of long read data: longer, but less accurate sequences, or highly accurate, but shorter reads (i.e., >Q20 or 99% accurate). To better understand how these types of long-read data as well as scale of data (i.e., mean length and sequencing depth) influence genome assembly outcomes, we compared genome assemblies for a caddisfly, Hesperophylax magnus, generated with longer, but less accurate, Oxford Nanopore (ONT) R9.4.1 and highly accurate PacBio HiFi (HiFi) data. Next, we expanded this comparison to consider the influence of highly accurate long-read sequence data on genome assemblies across 6750 plant and animal genomes. For this broader comparison, we used HiFi data as a surrogate for highly accurate long-reads broadly as we could identify when they were used from GenBank metadata., Results: HiFi reads outperformed ONT reads in all assembly metrics tested for the caddisfly data set and allowed for accurate assembly of the repetitive ~ 20 Kb H-fibroin gene. Across plants and animals, genome assemblies that incorporated HiFi reads were also more contiguous. For plants, the average HiFi assembly was 501% more contiguous (mean contig N50 = 20.5 Mb) than those generated with any other long-read data (mean contig N50 = 4.1 Mb). For animals, HiFi assemblies were 226% more contiguous (mean contig N50 = 20.9 Mb) versus other long-read assemblies (mean contig N50 = 9.3 Mb). In plants, we also found limited evidence that HiFi may offer a unique solution for overcoming genomic complexity that scales with assembly size., Conclusions: Highly accurate long-reads generated with HiFi or analogous technologies represent a key tool for maximizing genome assembly quality for a wide swath of plants and animals. This finding is particularly important when resources only allow for one type of sequencing data to be generated. Ultimately, to realize the promise of biodiversity genomics, we call for greater uptake of highly accurate long-reads in future studies., (© 2023. The Author(s).)

Long-read sequencing is driving a new reality for genome science in which highly contiguous assemblies can be produced efficiently with modest resources. Genome assemblies from long-read sequences are particularly exciting for understanding the evolution of complex genomic regions that are often difficult to assemble. In this study, we utilized long-read sequencing data to generate a high-quality genome assembly for an Antarctic eelpout, Ophthalmolycus amberensis, the first for the globally distributed family Zoarcidae. We used this assembly to understand how O. amberensis has adapted to the harsh Southern Ocean and compared it to another group of Antarctic fishes: the notothenioids. We showed that selection has largely acted on different targets in eelpouts relative to notothenioids. However, we did find some overlap; in both groups, genes involved in membrane structure, thermal tolerance and vision have evidence of positive selection. We found evidence for historical shifts of transposable element activity in O. amberensis and other polar fishes, perhaps reflecting a response to environmental change. We were specifically interested in the evolution of two complex genomic loci known to underlie key adaptations to polar seas: haemoglobin and antifreeze proteins (AFPs). We observed unique evolution of the haemoglobin MN cluster in eelpouts and related fishes in the suborder Zoarcoidei relative to other Perciformes. For AFPs, we identified the first species in the suborder with no evidence of afpIII sequences (Cebidichthys violaceus) in the genomic region where they are found in all other Zoarcoidei, potentially reflecting a lineage-specific loss of this cluster. Beyond polar fishes, our results highlight the power of long-read sequencing to understand genome evolution., (© 2022 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Glaciers are rapidly receding under climate change. A melting cryosphere will dramatically alter global sea levels, carbon cycling, and water resource availability. Glaciers host rich biotic communities that are dominated by microbial diversity, and this biodiversity can impact surface albedo, thereby driving a feedback loop between biodiversity and cryosphere melt. However, the microbial diversity of glacier ecosystems remains largely unknown outside of major ice sheets, particularly from a temporal perspective. Here, we characterized temporal dynamics of bacteria, eukaryotes, and algae on the Paradise Glacier, Mount Rainier, USA, over nine time points spanning the summer melt season. During our study, the glacier surface steadily darkened as seasonal snow melted and darkening agents accumulated until new snow fell in late September. From a community-wide perspective, the bacterial community remained generally constant while eukaryotes and algae exhibited temporal progression and community turnover. Patterns of individual taxonomic groups, however, were highly stochastic. We found little support for our a priori prediction that autotroph abundance would peak before heterotrophs. Notably, two different trends in snow algae emerged-an abundant early- and late-season operational taxonomic unit (OTU) with a different midsummer OTU that peaked in August. Overall, our results highlight the need for temporal sampling to clarify microbial diversity on glaciers and that caution should be exercised when interpreting results from single or few time points. IMPORTANCE Microbial diversity on mountain glaciers is an underexplored component of global biodiversity. Microbial presence and activity can also reduce the surface albedo or reflectiveness of glaciers, causing them to absorb more solar radiation and melt faster, which in turn drives more microbial activity. To date, most explorations of microbial diversity in the mountain cryosphere have only included single time points or focused on one microbial community (e.g., bacteria). Here, we performed temporal sampling over a summer melt season for the full microbial community, including bacteria, eukaryotes, and fungi, on the Paradise Glacier, Washington, USA. Over the summer, the bacterial community remained generally constant, whereas eukaryote and algal communities temporally changed through the melt season. Individual taxonomic groups, however, exhibited considerable stochasticity. Overall, our results highlight the need for temporal sampling on glaciers and that caution should be exercised when interpreting results from single or few time points.

Insect silk is a versatile biomaterial. Lepidoptera and Trichoptera display some of the most diverse uses of silk, with varying strength, adhesive qualities, and elastic properties. Silk fibroin genes are long (>20 Kbp), with many repetitive motifs that make them challenging to sequence. Most research thus far has focused on conserved N- and C-terminal regions of fibroin genes because a full comparison of repetitive regions across taxa has not been possible. Using the PacBio Sequel II system and SMRT sequencing, we generated high fidelity (HiFi) long-read genomic and transcriptomic sequences for the Indianmeal moth ( Plodia interpunctella ) and genomic sequences for the caddisfly Eubasilissa regina . Both genomes were highly contiguous (N50  = 9.7 Mbp/32.4 Mbp, L50  = 13/11) and complete (BUSCO complete  = 99.3%/95.2%), with complete and contiguous recovery of silk heavy fibroin gene sequences. We show that HiFi long-read sequencing is helpful for understanding genes with long, repetitive regions., Competing Interests: The authors declare that they have no competing interests., (© The Author(s) 2022.)

In less than 25 y, the field of animal genome science has transformed from a discipline seeking its first glimpses into genome sequences across the Tree of Life to a global enterprise with ambitions to sequence genomes for all of Earth's eukaryotic diversity [H. A. Lewin et al. , Proc. Natl. Acad. Sci. U.S.A. 115, 4325-4333 (2018)]. As the field rapidly moves forward, it is important to take stock of the progress that has been made to best inform the discipline's future. In this Perspective, we provide a contemporary, quantitative overview of animal genome sequencing. We identified the best available genome assemblies in GenBank, the world's most extensive genetic database, for 3,278 unique animal species across 24 phyla. We assessed taxonomic representation, assembly quality, and annotation status for major clades. We show that while tremendous taxonomic progress has occurred, stark disparities in genomic representation exist, highlighted by a systemic overrepresentation of vertebrates and underrepresentation of arthropods. In terms of assembly quality, long-read sequencing has dramatically improved contiguity, whereas gene annotations are available for just 34.3% of taxa. Furthermore, we show that animal genome science has diversified in recent years with an ever-expanding pool of researchers participating. However, the field still appears to be dominated by institutions in the Global North, which have been listed as the submitting institution for 77% of all assemblies. We conclude by offering recommendations for improving genomic resource availability and research value while also broadening global representation., Competing Interests: The authors declare no competing interest., (Copyright © 2021 the Author(s). Published by PNAS.)

The field of plant genome sequencing has grown rapidly in the past 20 years, leading to increases in the quantity and quality of publicly available genomic resources. The growing wealth of genomic data from an increasingly diverse set of taxa provides unprecedented potential to better understand the genome biology and evolution of land plants. Here we provide a contemporary view of land plant genomics, including analyses on assembly quality, taxonomic distribution of sequenced species and national participation. We show that assembly quality has increased dramatically in recent years, that substantial taxonomic gaps exist and that the field has been dominated by affluent nations in the Global North and China, despite a wide geographic distribution of study species. We identify numerous disconnects between the native range of focal species and the national affiliation of the researchers studying them, which we argue are rooted in colonialism-both past and present. Luckily, falling sequencing costs, widening availability of analytical tools and an increasingly connected scientific community provide key opportunities to improve existing assemblies, fill sampling gaps and empower a more global plant genomics community., (© 2021. The Author(s).)

Modern genetic data sets present unprecedented opportunities to understand the evolutionary origins of diverse taxonomic groups. When the timing of key events is known, it is possible to investigate biogeographic history in the context of major phenomena (e.g., cooling of a major ocean). In this study, we investigated the biogeographic history of the suborder Zoarcoidei, a globally distributed fish group that includes species inhabiting both poles that produce antifreeze proteins to survive chronic subfreezing temperatures. We first generated a multi-locus, time-calibrated phylogeny for the group. We then used biogeographic modeling to reconstruct ancestral ranges across the tree and to quantify the type and frequency of biogeographic events (e.g., founder, dispersal). With these results, we considered how the cooling of the Southern and Arctic Oceans, which reached their present-day subfreezing temperatures 10-15 million years ago (Mya) and 2-3 Mya, respectively, may have shaped the group's evolutionary history, with an emphasis on the most speciose and widely distributed family, eelpouts (family Zoarcidae). Our phylogenetic results clarified the Zoarcoidei taxonomy and showed that the group began to diversify in the Oligocene ~31-32 Mya, with the center of origin for all families in north temperate waters. Within-area speciation was the most common biogeographic event in the group's history (80% of all events) followed by dispersal (20%). Finally, we only found evidence, albeit limited, for ocean cooling underpinning diversification of eelpouts living in the high Antarctic over the last 10 million years., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Peer-review and subject-matter editing is the backbone of scientific publishing. However, early-career researchers (ECRs) are given few opportunities to participate in the editorial process beyond reviewing articles. Thus, a disconnect exists: science needs high-quality editorial talent to conduct, oversee and improve the publishing process, yet we dedicate few resources to building editorial talent nor giving ECRs formal opportunities to influence publishing from within. ECRs can contribute to the publishing landscape in unique ways given their insight into new and rapidly developing publishing trends (e.g. open science). Here, we describe a two-way fellowship model that gives ECRs a "seat" at the editorial table of a field-leading journal. We describe both the necessary framework and benefits that can stem from editorial fellowships for ECRs, editors, journals, societies, and the ​broader scientific community., (© 2021 John Wiley & Sons Ltd. This article has been contributed to by US Government employees and their work is in the public domain in the USA.)

The North American tiger salamander species complex, including its best-known species, the Mexican axolotl, has long been a source of biological fascination. The complex exhibits a wide range of variation in developmental life history strategies, including populations and individuals that undergo metamorphosis; those able to forego metamorphosis and retain a larval, aquatic lifestyle (i.e., paedomorphosis); and those that do both. The evolution of a paedomorphic life history state is thought to lead to increased population genetic differentiation and ultimately reproductive isolation and speciation, but the degree to which it has shaped population- and species-level divergence is poorly understood. Using a large multilocus dataset from hundreds of samples across North America, we identified genetic clusters across the geographic range of the tiger salamander complex. These clusters often contain a mixture of paedomorphic and metamorphic taxa, indicating that geographic isolation has played a larger role in lineage divergence than paedomorphosis in this system. This conclusion is bolstered by geography-informed analyses indicating no effect of life history strategy on population genetic differentiation and by model-based population genetic analyses demonstrating gene flow between adjacent metamorphic and paedomorphic populations. This fine-scale genetic perspective on life history variation establishes a framework for understanding how plasticity, local adaptation, and gene flow contribute to lineage divergence. Many members of the tiger salamander complex are endangered, and the Mexican axolotl is an important model system in regenerative and biomedical research. Our results chart a course for more informed use of these taxa in experimental, ecological, and conservation research., Competing Interests: The authors declare no competing interest.

Mountains are global biodiversity hotspots where cold environments and their associated ecological communities are threatened by climate warming. Considerable research attention has been devoted to understanding the ecological effects of alpine glacier and snowfield recession. However, much less attention has been given to identifying climate refugia in mountain ecosystems where present-day environmental conditions will be maintained, at least in the near-term, as other habitats change. Around the world, montane communities of microbes, animals, and plants live on, adjacent to, and downstream of rock glaciers and related cold rocky landforms (CRL). These geomorphological features have been overlooked in the ecological literature despite being extremely common in mountain ranges worldwide with a propensity to support cold and stable habitats for aquatic and terrestrial biodiversity. CRLs are less responsive to atmospheric warming than alpine glaciers and snowfields due to the insulating nature and thermal inertia of their debris cover paired with their internal ventilation patterns. Thus, CRLs are likely to remain on the landscape after adjacent glaciers and snowfields have melted, thereby providing longer-term cold habitat for biodiversity living on and downstream of them. Here, we show that CRLs will likely act as key climate refugia for terrestrial and aquatic biodiversity in mountain ecosystems, offer guidelines for incorporating CRLs into conservation practices, and identify areas for future research., (© 2021 John Wiley & Sons Ltd.)

Alpine regions are changing rapidly due to loss of snow and ice in response to ongoing climate change. While studies have documented ecological responses in alpine lakes and streams to these changes, our ability to predict such outcomes is limited. We propose that the application of fundamental rules of life can help develop necessary predictive frameworks. We focus on four key rules of life and their interactions: the temperature dependence of biotic processes from enzymes to evolution; the wavelength dependence of the effects of solar radiation on biological and ecological processes; the ramifications of the non-arbitrary elemental stoichiometry of life; and maximization of limiting resource use efficiency across scales. As the cryosphere melts and thaws, alpine lakes and streams will experience major changes in temperature regimes, absolute and relative inputs of solar radiation in ultraviolet and photosynthetically active radiation, and relative supplies of resources (e.g., carbon, nitrogen, and phosphorus), leading to nonlinear and interactive effects on particular biota, as well as on community and ecosystem properties. We propose that applying these key rules of life to cryosphere-influenced ecosystems will reduce uncertainties about the impacts of global change and help develop an integrated global view of rapidly changing alpine environments. However, doing so will require intensive interdisciplinary collaboration and international cooperation. More broadly, the alpine cryosphere is an example of a system where improving our understanding of mechanistic underpinnings of living systems might transform our ability to predict and mitigate the impacts of ongoing global change across the daunting scope of diversity in Earth's biota and environments., (© 2020 John Wiley & Sons Ltd.)

Climate change is altering conditions in high-elevation streams worldwide, with largely unknown effects on resident communities of aquatic insects. Here, we review the challenges of climate change for high-elevation aquatic insects and how they may respond, focusing on current gaps in knowledge. Understanding current effects and predicting future impacts will depend on progress in three areas. First, we need better descriptions of the multivariate physical challenges and interactions among challenges in high-elevation streams, which include low but rising temperatures, low oxygen supply and increasing oxygen demand, high and rising exposure to ultraviolet radiation, low ionic strength, and variable but shifting flow regimes. These factors are often studied in isolation even though they covary in nature and interact in space and time. Second, we need a better mechanistic understanding of how physical conditions in streams drive the performance of individual insects. Environment-performance links are mediated by physiology and behavior, which are poorly known in high-elevation taxa. Third, we need to define the scope and importance of potential responses across levels of biological organization. Short-term responses are defined by the tolerances of individuals, their capacities to perform adequately across a range of conditions, and behaviors used to exploit local, fine-scale variation in abiotic factors. Longer term responses to climate change, however, may include individual plasticity and evolution of populations. Whether high-elevation aquatic insects can mitigate climatic risks via these pathways is largely unknown., (© 2020 John Wiley & Sons Ltd.)

Rapid glacier recession is altering the physical conditions of headwater streams. Stream temperatures are predicted to rise and become increasingly variable, putting entire meltwater-associated biological communities at risk of extinction. Thus, there is a pressing need to understand how thermal stress affects mountain stream insects, particularly where glaciers are likely to vanish on contemporary timescales. In this study, we measured the critical thermal maximum (CT MAX ) of stonefly nymphs representing multiple species and a range of thermal regimes in the high Rocky Mountains, USA. We then collected RNA-sequencing data to assess how organismal thermal stress translated to the cellular level. Our focal species included the meltwater stonefly, Lednia tumana, which was recently listed under the U.S. Endangered Species Act due to climate-induced habitat loss. For all study species, critical thermal maxima (CT MAX  > 20°C) far exceeded the stream temperatures mountain stoneflies experience (<10°C). Moreover, while evidence for a cellular stress response was present, we also observed constitutive expression of genes encoding proteins known to underlie thermal stress (i.e., heat shock proteins) even at low temperatures that reflected natural conditions. We show that high-elevation aquatic insects may not be physiologically threatened by short-term exposure to warm temperatures and that longer-term physiological responses or biotic factors (e.g., competition) may better explain their extreme distributions., (© 2020 John Wiley & Sons Ltd.)

Climate change is proceeding rapidly in high mountain regions worldwide. Rising temperatures will impact insect physiology and associated fitness and will shift populations in space and time, thereby altering community interactions and composition. Shifts in space are expected as insects move upslope to escape warming temperatures and shifts in time will occur with changes in phenology of resident high-elevation insects. Clearly, spatiotemporal shifts will not affect all species equally. Terrestrial insects may have more opportunities than aquatic insects to exploit microhabitats, potentially buffering them from warming. Such responses of insects to warming may also fuel evolutionary change, including hitchhiking of maladaptive alleles and genetic rescue. Together, these considerations suggest a striking restructuring of high-elevation insect communities that remains largely unstudied., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Aquatic insects comprise 10% of all insect diversity, can be found on every continent except Antarctica, and are key components of freshwater ecosystems. However, aquatic insect genome biology lags dramatically behind that of terrestrial insects. If genomic effort was spread evenly, one aquatic insect genome would be sequenced for every ~9 terrestrial insect genomes. Instead, ~24 terrestrial insect genomes have been sequenced for every aquatic insect genome. This discrepancy is even more dramatic if the quality of genomic resources is considered; for instance, while no aquatic insect genome has been assembled to the chromosome level, 29 terrestrial insect genomes spanning four orders have. We argue that a lack of aquatic insect genomes is not due to any underlying difficulty (e.g., small body sizes or unusually large genomes), yet it is severely hampering aquatic insect research at both fundamental and applied scales. By expanding the availability of aquatic insect genomes, we will gain key insight into insect diversification and empower future research for a globally important taxonomic group.

Aquatic invertebrates are excellent indicators of ecosystem quality; however, choosing a sampling method can be difficult. Each method and associated protocol has advantages and disadvantages, and finding the approach that minimizes biases yet fulfills management objectives is crucial. To test the effects of both sampling methods and sample handling-i.e., to composite samples or leave them as replicates-we collected aquatic invertebrates from the Niobrara River at Agate Fossil Beds National Monument, Nebraska, using three methods and two sample handling protocols. We compared aquatic invertebrate assemblages collected with a Hester-Dendy multi-plate sampler, Hess sampler, and a D-frame dipnet. We calculated six common bioassessment metrics from composite (combined) and replicate (separate) samples. Hess samples contained the highest taxonomic richness (capturing 77% of all taxa observed) and dipnet samples the least (47%). Hester-Dendy samples had the greatest proportion of Ephemeroptera, and Ephemeroptera, Plecoptera, and Trichoptera (EPT). Dipnet samples had the lowest evenness values. In terms of sample handling, composite samples had inflated richness, diversity, and evenness compared with replicate samples, but bioassessment metrics calculated from proportions or averages (i.e., Hilsenhoff's Biotic Index and the proportion of EPT taxa) did not differ between them. The proportion of invertebrate groups from composite samples were not statistically different among sampling methods, but several groups differed between replicate samples collected by different methods. Ultimately, we recommend collecting replicate samples with a Hess sampler when the goal of the study is to detect ecosystem change, among locations or differences in variables of interest.

Glaciers are important drivers of environmental heterogeneity and biological diversity across mountain landscapes. Worldwide, glaciers are receding rapidly due to climate change, with important consequences for biodiversity in mountain ecosystems. However, the effects of glacier loss on biodiversity have never been quantified across a mountainous region, primarily due to a lack of adequate data at large spatial and temporal scales. Here, we combine high-resolution biological and glacier change (ca. 1850-2015) datasets for Glacier National Park, USA, to test the prediction that glacier retreat reduces biodiversity in mountain ecosystems through the loss of uniquely adapted meltwater stream species. We identified a specialized cold-water invertebrate community restricted to the highest elevation streams primarily below glaciers, but also snowfields and groundwater springs. We show that this community and endemic species have unexpectedly persisted in cold, high-elevation sites, even in catchments that have not been glaciated in ∼170 y. Future projections suggest substantial declines in suitable habitat, but not necessarily loss of this community with the complete disappearance of glaciers. Our findings demonstrate that high-elevation streams fed by snow and other cold-water sources continue to serve as critical climate refugia for mountain biodiversity even after glaciers disappear., Competing Interests: The authors declare no competing interest., (Copyright © 2020 the Author(s). Published by PNAS.)

Alpine streams are dynamic habitats harboring substantial biodiversity across small spatial extents. The diversity of alpine stream biota is largely reflective of environmental heterogeneity stemming from varying hydrological sources. Globally, alpine stream diversity is under threat as meltwater sources recede and stream conditions become increasingly homogeneous. Much attention has been devoted to macroinvertebrate diversity in alpine headwaters, yet to fully understand the breadth of climate change threats, a more thorough accounting of microbial diversity is needed. We characterized microbial diversity (specifically Bacteria and Archaea) of 13 streams in two disjunct Rocky Mountain subranges through 16S rRNA gene sequencing. Our study encompassed the spectrum of alpine stream sources (glaciers, snowfields, subterranean ice, and groundwater) and three microhabitats (ice, biofilms, and streamwater). We observed no difference in regional (γ) diversity between subranges but substantial differences in diversity among (β) stream types and microhabitats. Within-stream (α) diversity was highest in groundwater-fed springs, lowest in glacier-fed streams, and positively correlated with water temperature for both streamwater and biofilm assemblages. We identified an underappreciated alpine stream type-the icy seep-that are fed by subterranean ice, exhibit cold temperatures (summer mean <2°C), moderate bed stability, and relatively high conductivity. Icy seeps will likely be important for combatting biodiversity losses as they contain similar microbial assemblages to streams fed by surface ice yet may be buffered against climate change by insulating debris cover. Our results show that the patterns of microbial diversity support an ominous trend for alpine stream biodiversity; as meltwater sources decline, stream communities will become more diverse locally, but regional diversity will be lost. Icy seeps, however, represent a source of optimism for the future of biodiversity in these imperiled ecosystems., (© 2019 John Wiley & Sons Ltd.)

Disentangling the contemporary and historical factors underlying the spatial distributions of species is a central goal of biogeography. For species with broad distributions but little capacity to actively disperse, disconnected geographical distributions highlight the potential influence of passive, long-distance dispersal (LDD) on their evolutionary histories. However, dispersal alone cannot completely account for the biogeography of any species, and other factors-e.g. habitat suitability, life history-must also be considered. North American ice worms ( Mesenchytraeus solifugus) are ice-obligate annelids that inhabit coastal glaciers from Oregon to Alaska. Previous studies identified a complex biogeographic history for ice worms, with evidence for genetic isolation, unexpectedly close relationships among geographically disjunct lineages, and contemporary migration across large (e.g. greater than 1500 km) areas of unsuitable habitat. In this study, we analysed genome-scale sequence data for individuals from most of the known ice worm range. We found clear support for divergence between populations along the Pacific Coast and the inland flanks of the Coast Mountains (mean F ST = 0.60), likely precipitated by episodic ice sheet expansion and contraction during the Pleistocene. We also found support for LDD of ice worms from Alaska to Vancouver Island, perhaps mediated by migrating birds. Our results highlight the power of genomic data for disentangling complex biogeographic patterns, including the presence of LDD.

Extreme environments typically require costly adaptations for survival, an attribute that often translates to an elevated influence of habitat conditions on biotic communities. Microbes, primarily bacteria, are successful colonizers of extreme environments worldwide, yet in many instances, the interplay between harsh conditions, dispersal, and microbial biogeography remains unclear. This lack of clarity is particularly true for habitats where extreme temperature is not the overarching stressor, highlighting a need for studies that focus on the role other primary stressors (e.g., toxicants) play in shaping biogeographic patterns. In this study, we leveraged a naturally paired stream system in southern Mexico to explore how elevated hydrogen sulfide (H 2 S) influences microbial diversity. We sequenced a portion of the 16S rRNA gene using bacterial primers for water sampled from three geographically proximate pairings of streams with high (> 20 μM) or low (~ 0 μM) H 2 S concentrations. After exploring bacterial diversity within and among sites, we compared our results to a previous study of macroinvertebrates and fish for the same sites. By spanning multiple organismal groups, we were able to illuminate how H 2 S may differentially affect biodiversity. The presence of elevated H 2 S had no effect on overall bacterial diversity (p = 0.21), a large effect on community composition (25.8% of variation explained, p < 0.0001), and variable influence depending upon the group-whether fish, macroinvertebrates, or bacteria-being considered. For bacterial diversity, we recovered nine abundant operational taxonomic units (OTUs) that comprised a core H 2 S-rich stream microbiome in the region. Many H 2 S-associated OTUs were members of the Epsilonproteobacteria and Gammaproteobacteria, which both have been implicated in endosymbiotic relationships between sulfur-oxidizing bacteria and eukaryotes, suggesting the potential for symbioses that remain to be discovered in these habitats.

In alpine regions worldwide, climate change is dramatically altering ecosystems and affecting biodiversity in many ways. For streams, receding alpine glaciers and snowfields, paired with altered precipitation regimes, are driving shifts in hydrology, species distributions, basal resources, and threatening the very existence of some habitats and biota. Alpine streams harbour substantial species and genetic diversity due to significant habitat insularity and environmental heterogeneity. Climate change is expected to affect alpine stream biodiversity across many levels of biological resolution from micro- to macroscopic organisms and genes to communities. Herein, we describe the current state of alpine stream biology from an organism-focused perspective. We begin by reviewing seven standard and emerging approaches that combine to form the current state of the discipline. We follow with a call for increased synthesis across existing approaches to improve understanding of how these imperiled ecosystems are responding to rapid environmental change. We then take a forward-looking viewpoint on how alpine stream biologists can make better use of existing data sets through temporal comparisons, integrate remote sensing and geographic information system (GIS) technologies, and apply genomic tools to refine knowledge of underlying evolutionary processes. We conclude with comments about the future of biodiversity conservation in alpine streams to confront the daunting challenge of mitigating the effects of rapid environmental change in these sentinel ecosystems., (© 2017 Cambridge Philosophical Society.)

Glacier ecosystems are teeming with life on, beneath, and to a lesser degree, within their icy masses. This conclusion largely stems from polar research, with less attention paid to mountain glaciers that overlap environmentally and ecologically with their polar counterparts in some ways, but diverge in others. One difference lies in the susceptibility of mountain glaciers to the near-term threat of climate change, as they tend to be much smaller in both area and volume. Moreover, mountain glaciers are typically steeper, more dependent upon basal sliding for movement, and experience higher seasonal precipitation. Here, we provide a modern synthesis of the microbial ecology of mountain glacier ecosystems, and particularly those at low- to mid-latitudes. We focus on five ecological zones: the supraglacial surface, englacial interior, subglacial bedrock-ice interface, proglacial streams and glacier forefields. For each, we discuss the role of microbiota in biogeochemical cycling and outline ecological and hydrological connections among zones, underscoring the interconnected nature of these ecosystems. Collectively, we highlight the need to: better document the biodiversity and functional roles of mountain glacier microbiota; describe the ecological implications of rapid glacial retreat under climate change and resolve the relative contributions of ecological zones to broader ecosystem function., (© 2017 The Authors. Environmental Microbiology published by Society for Applied Microbiology and John Wiley & Sons Ltd.)

Climate warming is causing rapid loss of glaciers and snowpack in mountainous regions worldwide. These changes are predicted to negatively impact the habitats of many range-restricted species, particularly endemic, mountaintop species dependent on the unique thermal and hydrologic conditions found only in glacier-fed and snow melt-driven alpine streams. Although progress has been made, existing understanding of the status, distribution, and ecology of alpine aquatic species, particularly in North America, is lacking, thereby hindering conservation and management programs. Two aquatic insects - the meltwater stonefly (Lednia tumana) and the glacier stonefly (Zapada glacier) - were recently proposed for listing under the U.S. Endangered Species Act due to climate-change-induced habitat loss. Using a large dataset (272 streams, 482 total sites) with high-resolution climate and habitat information, we describe the distribution, status, and key environmental features that limit L. tumana and Z. glacier across the northern Rocky Mountains. Lednia tumana was detected in 113 streams (175 sites) within Glacier National Park (GNP) and surrounding areas. The probability of L. tumana occurrence increased with cold stream temperatures and close proximity to glaciers and permanent snowfields. Similarly, densities of L. tumana declined with increasing distance from stream source. Zapada glacier was only detected in 10 streams (24 sites), six in GNP and four in mountain ranges up to ~600 km southwest. Our results show that both L. tumana and Z. glacier inhabit an extremely narrow distribution, restricted to short sections of cold, alpine streams often below glaciers predicted to disappear over the next two decades. Climate warming-induced glacier and snow loss clearly imperils the persistence of L. tumana and Z. glacier throughout their ranges, highlighting the role of mountaintop aquatic invertebrates as sentinels of climate change in mid-latitude regions., (© 2016 Published by John Wiley & Sons, Ltd. This article has been contributed to by US Government employees and their work is in the public domain in the USA.)

Perhaps the most important recent advance in species delimitation has been the development of model-based approaches to objectively diagnose species diversity from genetic data. Additionally, the growing accessibility of next-generation sequence data sets provides powerful insights into genome-wide patterns of divergence during speciation. However, applying complex models to large data sets is time-consuming and computationally costly, requiring careful consideration of the influence of both individual and population sampling, as well as the number and informativeness of loci on species delimitation conclusions. Here, we investigated how locus number and information content affect species delimitation results for an endangered Mexican salamander species, Ambystoma ordinarium. We compared results for an eight-locus, 137-individual data set and an 89-locus, seven-individual data set. For both data sets, we used species discovery methods to define delimitation models and species validation methods to rigorously test these hypotheses. We also used integrated demographic model selection tools to choose among delimitation models, while accounting for gene flow. Our results indicate that while cryptic lineages may be delimited with relatively few loci, sampling larger numbers of loci may be required to ensure that enough informative loci are available to accurately identify and validate shallow-scale divergences. These analyses highlight the importance of striking a balance between dense sampling of loci and individuals, particularly in shallowly diverged lineages. They also suggest the presence of a currently unrecognized, endangered species in the western part of A. ordinarium's range., (© 2016 John Wiley & Sons Ltd.)

Much remains unknown about the genetic status and population connectivity of high-elevation and high-latitude freshwater invertebrates, which often persist near snow and ice masses that are disappearing due to climate change. Here we report on the conservation genetics of the meltwater stonefly Lednia tumana (Ricker) of Montana, USA, a cold-water obligate species. We sequenced 1530 bp of mtDNA from 116 L. tumana individuals representing "historic" (>10 yr old) and 2010 populations. The dominant haplotype was common in both time periods, while the second-most-common haplotype was found only in historic samples, having been lost in the interim. The 2010 populations also showed reduced gene and nucleotide diversity and increased genetic isolation. We found lower genetic diversity in L. tumana compared to two other North American stonefly species, Amphinemura linda (Ricker) and Pteronarcys californica Newport. Our results imply small effective sizes, increased fragmentation, limited gene flow, and loss of genetic variation among contemporary L. tumana populations, which can lead to reduced adaptive capacity and increased extinction risk. This study reinforces concerns that ongoing glacier loss threatens the persistence of L. tumana, and provides baseline data and analysis of how future environmental change could impact populations of similar organisms.

Implementation of the coalescent model in a Bayesian framework is an emerging strength in genetically based species delimitation studies. By providing an objective measure of species diagnosis, these methods represent a quantitative enhancement to the analysis of multilocus data, and complement more traditional methods based on phenotypic and ecological characteristics. Recognized as two species 20 years ago, mouse lemurs (genus Microcebus) now comprise more than 20 species, largely diagnosed from mtDNA sequence data. With each new species description, enthusiasm has been tempered with scientific scepticism. Here, we present a statistically justified and unbiased Bayesian approach towards mouse lemur species delimitation. We perform validation tests using multilocus sequence data and two methodologies: (i) reverse-jump Markov chain Monte Carlo sampling to assess the likelihood of different models defined a priori by a guide tree, and (ii) a Bayes factor delimitation test that compares different species-tree models without a guide tree. We assess the sensitivity of these methods using randomized individual assignments, which has been used in bpp studies, but not with Bayes factor delimitation tests. Our results validate previously diagnosed taxa, as well as new species hypotheses, resulting in support for three new mouse lemur species. As the challenge of multiple researchers using differing criteria to describe diversity is not unique to Microcebus, the methods used here have significant potential for clarifying diversity in other taxonomic groups. We echo previous studies in advocating that multiple lines of evidence, including use of the coalescent model, should be trusted to delimit new species., (© 2016 John Wiley & Sons Ltd.)