Ecological thresholds under atmospheric nitrogen deposition for 1200 herbaceous species and 24 communities across the United States.

Nitrogen (N) emissions and atmospheric deposition have increased significantly during the last century and become a stressor for many N-sensitive plant species. Understanding individual and community herbaceous plant species thresholds to atmospheric N deposition can inform emissions reduction policy. Here, we present results using Threshold Indicator Taxa Analysis (TITAN) applied to more than 1200 unique plant species and 24 vegetation communities (i.e., alliances) across the United States (US) to assess vulnerability to N deposition. Alliance-level thresholds (change points) for species decreasing in abundance along the gradient ranged from 1.8 to 14.3 kg N ha─1  year─1 and tended to be lower in the west than the east, which suggests that eastern communities, where N deposition has been historically higher, may have already lost many sensitive species. For the species that were present in more than one alliance, over half had a variable response to the N deposition gradient, suggesting that local factors affect vulnerability. Significant progress has been made during the past 30 years to reduce N emissions, which has reduced the percentage of plots at risk to N deposition from 72% to 35%. Nevertheless, over a third of plots remain at risk, and an average reduction of N deposition of 20% would protect half of the plots where N deposition exceeds community thresholds. Furthermore, the alliance- and species-level change points determined in this study may be used to inform N critical loads.

Arthropod dark taxa provide new insights into diversity responses to bark beetle infestations.

Natural disturbances are increasing around the globe, also impacting protected areas. Although previous studies have indicated that natural disturbances result in mainly positive effects on biodiversity, these analyses mostly focused on a few well established taxonomic groups, and thus uncertainty remains regarding the comprehensive impact of natural disturbances on biodiversity. Using Malaise traps and meta-barcoding, we studied a broad range of arthropod taxa, including dark and cryptic taxa, along a gradient of bark beetle disturbance severities in five European national parks. We identified order-level community thresholds of disturbance severity and classified barcode index numbers (BINs; a cluster system for DNA sequences, where each cluster corresponds to a species) as negative or positive disturbance indicators. Negative indicator BINs decreased above thresholds of low to medium disturbance severity (20%-30% of trees killed), whereas positive indicator BINs benefited from high disturbance severity (76%-98%). BINs allocated to a species name contained nearly as many positive as negative disturbance indicators, but dark and cryptic taxa, particularly Diptera and Hymenoptera in our data, contained higher numbers of negative disturbance indicator BINs. Analyses of changes in the richness of BINs showed variable responses of arthropods to disturbance severity at lower taxonomic levels, whereas no significant signal was detected at the order level due to the compensatory responses of the underlying taxa. We conclude that the analyses of dark taxa can offer new insights into biodiversity responses to disturbances. Our results suggest considerable potential for forest management to foster arthropod diversity, for example by maintaining both closed-canopy forests (>70% cover) and open forests (<30% cover) on the landscape.

Long-term abundance trends of insect taxa are only weakly correlated.

Changes in the abundances of animals, such as with the ongoing concern about insect declines, are often assumed to be general across taxa. However, this assumption is largely untested. Here, we used a database of assemblage-wide long-term insect and arachnid monitoring to compare abundance trends among co-occurring pairs of taxa. We show that 60% of co-occurring taxa qualitatively showed long-term trends in the same direction-either both increasing or both decreasing. However, in terms of magnitude, temporal trends were only weakly correlated (mean freshwater r = 0.05 (±0.03), mean terrestrial r = 0.12 (±0.09)). The strongest correlation was between trends of beetles and those of moths/butterflies (r = 0.26). Overall, even though there is some support for directional similarity in temporal trends, we find that changes in the abundance of one taxon provide little information on the changes of other taxa. No clear candidate for umbrella or indicator taxa emerged from our analysis. We conclude that obtaining a better picture of changes in insect abundances will require monitoring of multiple taxa, which remains uncommon, especially in the terrestrial realm.

Response of macroinvertebrate community to water quality factors and aquatic ecosystem health assessment in a typical river in Beijing, China.

Healthy aquatic ecosystems can offer basic ecological services for the sustainable development of humans and society. Water quality greatly influences the macroinvertebrate community in aquatic ecosystems and can alter the aquatic ecosystem's health status. However, the quantitative relationship between macroinvertebrate community and water quality factors in rivers remains unclear, particularly in urban rivers, which are strongly affected by human activities. Therefore, a new framework for the quantitative analysis between macroinvertebrate community and key water quality driving factors was developed in the study, meanwhile, the aquatic ecosystem health conditions were evaluated and validated by different methods. The framework was applied to a typical urban river, the North Canal River, which is regarded as the "mother river" of Beijing. Combined with the redundancy analysis (RDA) and the threshold indicator taxa analysis (TITAN), the water quality driving factors and their indicator species were identified and the corresponding response threshold was determined. Based on the benthic index of biotic integrity (B-IBI), the multi-metric rapid bioassessment method, and the biological monitoring working party (BMWP) score, the aquatic ecosystem health condition in the basin was comprehensively evaluated. The results show that fluoride, biochemical oxygen demand, ammonia-nitrogen and total phosphorus were the key water quality driving factors influencing the community structure of macroinvertebrates. Four indicator species of ammonia-nitrogen were identified by the TITAN method with a threshold range of 1.09-6.94 mg L-1, and three indicator species of total phosphorus were identified with a threshold range of 0.48-1.27 mg L-1. According to the results of the aquatic ecosystem health assessment, the river ecosystem was generally unhealthy and the upstream was better than downstream; the health condition in the mountainous areas of Changping district was the best, while that in Chaoyang district and the central city area was the worst. The framework could provide a strong basis for ecological restoration and pollution control of the urban rivers and become an important tool for the rehabilitation of aquatic ecosystems.

DNA metabarcoding effectively quantifies diatom responses to nutrients in streams.

Nutrient pollution from human activities remains a common problem facing stream ecosystems. Identifying ecological responses to phosphorus and nitrogen can inform decisions affecting the protection and management of streams and their watersheds. Diatoms are particularly useful because they are a highly diverse group of unicellular algae found in nearly all aquatic environments and are sensitive responders to increased nutrient concentrations. Here, we used DNA metabarcoding of stream diatoms as an approach to quantifying effects of total phosphorus (TP) and total nitrogen (TN). Threshold indicator taxa analysis (TITAN) identified operational taxonomic units (OTUs) that increased or decreased along TP and TN gradients along with nutrient concentrations at which assemblages had substantial changes in the occurrences and relative abundances of OTUs. Boosted regression trees showed that relative abundances of gene sequence reads for OTUs identified by TITAN as low P, high P, low N, or high N diatoms had strong relationships with nutrient concentrations, which provided support for potentially using these groups of diatoms as metrics in monitoring programs. Gradient forest analysis provided complementary information by characterizing multi-taxa assemblage change using multiple predictors and results from random forest models for each OTU. Collectively, these analyses showed that notable changes in diatom assemblage structure and OTUs began around 20 µg TP/L, low P diatoms decreased substantially and community change points occurred from 75 to 150 µg/L, and high P diatoms became increasingly dominant from 150 to 300 µg/L. Diatoms also responded to TN with large decreases in low N diatoms occurring from 280 to 525 µg TN/L and a transition to dominance by high N diatoms from 525-850 µg/L. These diatom responses to TP and TN could be used to inform protection efforts (i.e., anti-degradation) and management goals (i.e., nutrient reduction) in streams and watersheds. Our results add to the growing support for using diatom metabarcoding in monitoring programs.

Water Beetles as Models in Ecology and Evolution.

Beetles have colonized water many times during their history, with some of these events involving extensive evolutionary radiations and multiple transitions between land and water. With over 13,000 described species, they are one of the most diverse macroinvertebrate groups in most nonmarine aquatic habitats and occur on all continents except Antarctica. A combination of wide geographical and ecological range and relatively accessible taxonomy makes these insects an excellent model system for addressing a variety of questions in ecology and evolution. Work on water beetles has recently made important contributions to fields as diverse as DNA taxonomy, macroecology, historical biogeography, sexual selection, and conservation biology, as well as predicting organismal responses to global change. Aquatic beetles have some of the best resolved phylogenies of any comparably diverse insect group, and this, coupled with recent advances in taxonomic and ecological knowledge, is likely to drive an expansion of studies in the future.

Dung beetle-mammal associations: methods, research trends and future directions.

Dung beetles are increasingly used as a study taxon-both as bioindicators of environmental change, and as a model system for exploring ecosystem functioning. The advantages of this focal taxon approach are many; dung beetles are abundant in a wide range of terrestrial ecosystems, speciose, straightforward to sample, respond to environmental gradients and can be easily manipulated to explore species-functioning relationships. However, there remain large gaps in our understanding of the relationship between dung beetles and the mammals they rely on for dung. Here we review the literature, showing that despite an increase in the study of dung beetles linked to ecosystem functioning and to habitat and land use change, there has been little research into their associations with mammals. We summarize the methods and findings from dung beetle-mammal association studies to date, revealing that although empirical field studies of dung beetles rarely include mammal data, those that do, indicate mammal species presence and composition has a large impact on dung beetle species richness and abundance. We then review the methods used to carry out diet preference and ecosystem functioning studies, finding that despite the assumption that dung beetles are generalist feeders, there are few quantitative studies that directly address this. Together this suggests that conclusions about the effects of habitat change on dung beetles are based on incomplete knowledge. We provide recommendations for future work to identify the importance of considering mammal data for dung beetle distributions, composition and their contributions to ecosystem functioning; a critical step if dung beetles are to be used as a reliable bioindicator taxon.

Anthropogenic disturbance equalizes diversity levels in arbuscular mycorrhizal fungal communities.

The arbuscular mycorrhizal (AM) symbiosis is a key plant-microbe interaction in sustainable functioning ecosystems. Increasing anthropogenic disturbance poses a threat to AM fungal communities worldwide, but there is little empirical evidence about its potential negative consequences. In this global study, we sequenced AM fungal DNA in soil samples collected from pairs of natural (undisturbed) and anthropogenic (disturbed) plots in two ecosystem types (10 naturally wooded and six naturally unwooded ecosystems). We found that ecosystem type had stronger directional effects than anthropogenic disturbance on AM fungal alpha and beta diversity. However, disturbance increased alpha and beta diversity at sites where natural diversity was low and decreased diversity at sites where natural diversity was high. Cultured AM fungal taxa were more prevalent in anthropogenic than natural plots, probably due to their efficient colonization strategies and ability to recover from disturbance. We conclude that anthropogenic disturbance does not have a consistent directional effect on AM fungal diversity; rather, disturbance equalizes levels of diversity at large scales and causes changes in community functional structure.

Bacterial communities under long-term conventional and transgenic cotton farming systems using V3-V5 and V5-V9 of 16s rDNA.

Understanding the community structure of soil microbes is required to evaluate the potential effects of genetically modified (GM) plants on ecological environments. Bacterial communities in soil planted with conventional cotton (CC) and transgenic cultivar (TC) in a natural ecosystem for three years were characterized by 454 pyrosequencing of the V3-V5 and V5-V9 regions of 16S rDNA from June to September 2013. V3-V5 and V5-V9 regions yielded a total of 12,848 and 10,541 OTUs, respectively. The V5-V9 amplicon was additionally used to detect phyla that were poorly sequenced by V3-V5 (such as Chlamydiae, Crenarchaeota and Archaea). Among the species detected by each primer pair, 46% of the species identified from V3-V5 and 60% of those identified from V5-V9 were detected by both primer pairs. Although distinct bacterial compositions existed between the two amplified regions, statistical analysis revealed no significant difference in the diversity indexes or phylogenetic patterns in TC versus compared to those in the CC control. Further, clustering analysis in both regions indicated that there was no unambiguous aggregation in TC compared to that in CC control. Of all 26 phyla detected by both regions, each region detected 2 distinct phyla exhibiting significant variations in abundance. The species unique to each treatment field accounted for less than 27% of all species and were rare taxa (abundance < 0.15%). However, a small fraction of diagnostic taxa with specific ecological functions differed significantly between TC and CC. These differences were not driven by any obvious environmental factors. The results established a comprehensive inventory of the bacterial communities associated with GM plants and indicated that transgenic cotton may not significantly affect soil microorganisms compared with conventional cotton over a three-year period. Furthermore, diagnostic taxa were provided for monitoring the perturbation in soil, but further verification in future studies is required.

Temperature threshold models for benthic macroinvertebrates in Idaho wadeable streams and neighboring ecoregions.

Water temperatures are warming throughout the world including the Pacific Northwest, USA. Benthic macroinvertebrates are one of the most important and widely used indicators of freshwater impairment; however, their response to increased water temperatures and their use for monitoring water temperature impairment has been hindered by lack of knowledge of temperature occurrences, threshold change points, or indicator taxa. We present new analysis of a large macroinvertebrate database provided by Idaho Department of Environmental Quality from wadeable streams in Idaho that is to be used in conjunction with our previous analyses. This new analysis provides threshold change points for over 400 taxa along an increasing temperature gradient and provides a list of statistically important indicator taxa. The macroinvertebrate assemblage temperature change point for the taxa that decreased with increased temperatures was determined to be about 20.5 °C and for the taxa assemblage that increased with increased temperatures was about 11.5 °C. Results of this new analysis combined with our previous analysis will also be useful for others in neighboring regions where these taxa occur.

Influences of environmental factors on biomass of phytoplankton in the northern part of Tai Lake, China, from 2000 to 2012.

Long-term (2000 to 2012) monthly data on communities of phytoplankton, and environmental variables were measured in water collected from Meiliang Bay and Wuli Lake of Tai Lake, China. Redundancy analysis (RDA) was conducted to explore relationships between the phytoplankton communities and environmental variables. Change points for concentrations of nutrients, which serve as early warnings of state shifts in lacustrine ecosystems, were identified using the Threshold Indicator Taxa Analysis (TITAN). The biomass of phytoplankton was positively correlated with the concentrations of total phosphorus (TP), suspended solids (SS), water temperature (WT), and pH but negatively correlated with the N/P ratio (by mass) and Secchi disk depth (SD). Furthermore, TP, rather than other factors, was a controlling factor limiting the primary production of phytoplankton in most of this region. The change points for concentrations of TP controlling the occurrences of sensitive and tolerant taxa were 56.1 and 103.5 µg TP/L, respectively. These results imply that an abrupt change in this lacustrine ecosystem has occurred in most parts of the study area, and the turbid state of this lake can be altered by reducing TP loading. This study provides an alternative ecological method for exploring the production of algal blooms and could advance the understanding of HABs.

Metropolitan pressures: Significant biodiversity declines and strong filtering of functional traits in fish assemblages.

Accelerating global urbanization is leading to drastic losses and restructuring of biodiversity. Although it is crucial to understand urban impacts on biodiversity to develop mitigation strategies, there is a dearth of knowledge on the functional structure of fish assemblages spanning the entire city-scale spectrum of urbanization intensity. Here, using environmental DNA sampled from 109 water sites in Beijing, we investigated the taxonomic and functional diversity patterns of fish assemblages across the city and uncovered community-, trait-, and species-level responses to various environmental stressors. By ranking sampling sites into three disturbance levels according to water physiochemical and landcover conditions, we found that both native and non-native fish taxonomic and functional α-diversity decreased significantly with elevating disturbance, as strong disturbance led to the disappearance of many species. However, the quantitative taxonomic and functional ß-diversity components of native and non-native fish showed distinct patterns; assemblage turnover dominated native fish ß-diversity and decreased with increasing disturbance, whereas species/trait richness differences dominated non-native fish ß-diversity and increased with disturbance intensity particularly in lotic waters. RLQ and fourth-corner analyses revealed that fish size, fecundity, diet, and reproductive behaviors were significantly correlated with water quality, with pollution-tolerant, larger-sized native and omnivorous non-native fishes being urban winners, which indicates strong trait-dependent environmental filtering. Potential ecological indicator species were identified based on the sensitivity of fish responses to pollution loads; these were mostly small native species, and many have bivalve-dependent reproduction. Our results demonstrate that, along with native fish assemblage simplification and homogenization, urban stressors exert profound impacts on community trait composition, highlighting the need to consider both biodiversity loss and functional reorganization in combating disturbance of aquatic ecosystems under global urbanization. Furthermore, correlations between cropland cover and water nutrient level suggested that the management of agricultural runoff might be critically important for safeguarding urban water quality.

Maximising signal-to-noise ratios in environmental DNA-based monitoring.

Man's impacts on global ecosystems are increasing and there is a growing demand that these activities be appropriately monitored. Monitoring requires measurement of a response metric ('signal') that changes maximally and consistently in response to the monitored activity irrespective of other factors ('noise'), thus maximising the signal-to-noise ratio. Indices derived from time-consuming morphology-based taxonomic identification of organisms are a core part of many monitoring programmes. Metabarcoding is an alternative to morphology-based identification and involves the sequencing of short fragments of DNA ('markers') from multiple taxa simultaneously. DNA suitable for metabarcoding includes that extracted from environmental samples (eDNA). Metabarcoding outputs DNA sequences that can be identified (annotated) by matching them against archived annotated sequences. However, sequences from most organisms are not archived - preventing annotation and potentially limiting metabarcoding in monitoring applications. Consequently, there is growing interest in using unannotated sequences as response metrics in monitoring programmes. We compared the sequences from three commonly used markers (16S (V3/V4 regions), 18S (V1/V2 regions) and COI) and, sampling along steep impact gradients, showed that the 16S and COI sequences were associated with the largest and smallest signal-to-noise ratio respectively. We trialled four separate, intuitive, noise-reduction approaches and demonstrated that removing less frequent sequences improved the signal-to-noise ratio, partitioning an additional 25 % from noise to explanatory factors in non-parametric ANOVA (NPA) and reducing dispersion in the data. For the 16S marker, retaining only the most frequently observed sequence, per sample, resulting in nine sequences across 150 samples, generated a near-maximal signal-to-noise ratio (95 % of the variance explained in NPA). We recommend that NPA, combined with rigorous elimination of less frequent sequences, be used to pre-filter sequences/taxa being used in monitoring applications. Our approach will simplify downstream analysis, for example the identification of key taxa and functional associations.

Karst spring microbiome: Diversity, core taxa, and community response to pathogens and antibiotic resistance gene contamination.

Karst aquifers are important water resources for drinking water supplies worldwide. Although they are susceptible to anthropogenic contamination due to their high permeability, there is a lack of detailed knowledge on the stable core microbiome and how contamination may affect these communities. In this study, eight karst springs (distributed across three different regions in Romania) were sampled seasonally for one year. The core microbiota was analysed by 16S rRNA gene amplicon sequencing. To identify bacteria carrying antibiotic resistance genes and mobile genetic elements, an innovative method was applied, consisting of high-throughput antibiotic resistance gene quantification performed on potential pathogen colonies cultivated on Compact Dry™ plates. A taxonomically stable bacterial community consisting of Pseudomonadota, Bacteroidota, and Actinomycetota was revealed. Core analysis reaffirmed these results and revealed primarily freshwater-dwelling, psychrophilic/psychrotolerant species affiliated to Rhodoferax, Flavobacterium, and Pseudomonas genera. Both sequencing and cultivation methods indicated that more than half of the springs were contaminated with faecal bacteria and pathogens. These samples contained high levels of sulfonamide, macrolide, lincosamide and streptogramins B, and trimethoprim resistance genes spread primarily by transposase and insertion sequences. Differential abundance analysis found Synergistota, Mycoplasmatota, and Chlamydiota as suitable candidates for pollution monitoring in karst springs. This is the first study highlighting the applicability of a combined approach based on high-throughput SmartChip™ antibiotic resistance gene quantification and Compact Dry™ pathogen cultivation for estimating microbial contaminants in karst springs and other challenging low biomass environments.

Towards a scientific community consensus on designating Vulnerable Marine Ecosystems from imagery.

Management of deep-sea fisheries in areas beyond national jurisdiction by Regional Fisheries Management Organizations/Arrangements (RFMO/As) requires identification of areas with Vulnerable Marine Ecosystems (VMEs). Currently, fisheries data, including trawl and longline bycatch data, are used by many RFMO/As to inform the identification of VMEs. However, the collection of such data creates impacts and there is a need to collect non-invasive data for VME identification and monitoring purposes. Imagery data from scientific surveys satisfies this requirement, but there currently is no established framework for identifying VMEs from images. Thus, the goal of this study was to bring together a large international team to determine current VME assessment protocols and establish preliminary global consensus guidelines for identifying VMEs from images. An initial assessment showed a lack of consistency among RFMO/A regions regarding what is considered a VME indicator taxon, and hence variability in how VMEs might be defined. In certain cases, experts agreed that a VME could be identified from a single image, most often in areas of scleractinian reefs, dense octocoral gardens, multiple VME species' co-occurrence, and chemosynthetic ecosystems. A decision flow chart is presented that gives practical interpretation of the FAO criteria for single images. To further evaluate steps of the flow chart related to density, data were compiled to assess whether scientists perceived similar density thresholds across regions. The range of observed densities and the density values considered to be VMEs varied considerably by taxon, but in many cases, there was a statistical difference in what experts considered to be a VME compared to images not considered a VME. Further work is required to develop an areal extent index, to include a measure of confidence, and to increase our understanding of what levels of density and diversity correspond to key ecosystem functions for VME indicator taxa. Based on our results, the following recommendations are made: 1. There is a need to establish a global consensus on which taxa are VME indicators. 2. RFMO/As should consider adopting guidelines that use imagery surveys as an alternative (or complement) to using bycatch and trawl surveys for designating VMEs. 3. Imagery surveys should also be included in Impact Assessments. And 4. All industries that impact the seafloor, not just fisheries, should use imagery surveys to detect and identify VMEs.

Soil microbial communities and their co-occurrence networks in response to long-term Pb-Zn contaminated soil in southern China.

Mining causes extreme heavy metal (HM) contamination to surrounding environments and poses threats to soil microbial community. The effects of HMs on soil microbial communities are not only related to their total amounts but also associated with the distribution of chemical fractions. However, the effects of chemical fractions on soil microbes and their interactions remain largely unclear. Here we investigated soil physicochemical properties and bacterial and fungal communities of soil samples from the control area and lightly (L), moderately (M), and heavily (H) contaminated areas, respectively, which were collected from long-term Pb-Zn slag contamination area in the southern China. The results showed that bacterial and fungal community composition and structure were significantly affected by HMs, while community diversity was not significantly affected by HMs. The critical environmental factor affecting bacterial and fungal communities was pH, and the impacts of chemical fractions on their changes were more significant than the total amounts of HMs. Variance partitioning analysis (VPA) revealed fungal community changes were mostly driven by HM total amounts, but bacterial community changes were mostly driven by soil chemical properties. Co-occurrence network indicated that interactions among species of fungal network were sparser than that of bacterial network, but fungal network was more stable, due to a more significant number of keystone taxa and a lower percentage of positive associations. These illustrated that the fungal community might serve as indicator taxa for HM-contaminated status, and specific HM-responsive fungal species such as Triangularia mangenotii, Saitozyma podzolica, and Cladosporium endophytica, and genus Rhizophagus can be considered relevant bioindicators due to their less relative abundance in contaminated areas. Additionally, HM-responsive bacterial OTUs representing five genera within Sulfurifustis, Thiobacillus, Sphingomonas, Qipengyuania, and Sulfurirhabdus were found to be tolerant to HM stress due to their high relative abundance in contaminated levels.

Effects of Tree Composition and Soil Depth on Structure and Functionality of Belowground Microbial Communities in Temperate European Forests.

Depending on their tree species composition, forests recruit different soil microbial communities. Likewise, the vertical nutrient gradient along soil profiles impacts these communities and their activities. In forest soils, bacteria and fungi commonly compete, coexist, and interact, which is challenging for understanding the complex mechanisms behind microbial structuring. Using amplicon sequencing, we analyzed bacterial and fungal diversity in relation to forest composition and soil depth. Moreover, employing random forest models, we identified microbial indicator taxa of forest plots composed of either deciduous or evergreen trees, or their mixtures, as well as of three soil depths. We expected that forest composition and soil depth affect bacterial and fungal diversity and community structure differently. Indeed, relative abundances of microbial communities changed more across soil depths than in relation to forest composition. The microbial Shannon diversity was particularly affected by soil depth and by the proportion of evergreen trees. Our results also reflected that bacterial communities are primarily shaped by soil depth, while fungi were influenced by forest tree species composition. An increasing proportion of evergreen trees did not provoke differences in main bacterial metabolic functions, e.g., carbon fixation, degradation, or photosynthesis. However, significant responses related to specialized bacterial metabolisms were detected. Saprotrophic, arbuscular mycorrhizal, and plant pathogenic fungi were related to the proportion of evergreen trees, particularly in topsoil. Prominent microbial indicator taxa in the deciduous forests were characterized to be r-strategists, whereas K-strategists dominated evergreen plots. Considering simultaneously forest composition and soil depth to unravel differences in microbial communities, metabolic pathways and functional guilds have the potential to enlighten mechanisms that maintain forest soil functionality and provide resistance against disturbances.

[Quantitative Analysis of the Correlation Between Macrobenthos Community and Water Environmental Factors and Aquatic Ecosystem Health Assessment in the North Canal River Basin of Beijing].

Macrobenthos can reflect the cumulative effect of various ecological threats on the water environment and are closely related to the health of river ecosystems. In this study, taking the North Canal River basin, a typical basin in Beijing, as an example, ecological data from 34 stations were investigated in the summer of 2015. Characteristics of the macrobenthos communities were analyzed, and driving environmental factors were identified using typical correspondence analysis. Thresholds and response species of those driving environmental factors were conducted using the thresholds indicator taxa analysis (TITAN). In this study, the health status of the river ecosystem was evaluated by the multi-metrics method and benthic index of biotic integrity (B-IBI). The benthic community was dominated by pollution-tolerant aquatic insects and mollusks, with a low-level Shannon-wiener diversity index between 0-1.01; fluoride, biochemical oxygen demand, ammonia-nitrogen, and total phosphorus were driving environmental factors influencing the community structure of macrobenthos. Indicator species of ammonia-nitrogen were identified by the TITAN in the North Canal River basin with a threshold range of 1.09-6.94 mg·L-1; three indicator species of total phosphorus were identified with a threshold range of 0.48-1.27 mg·L-1, which were all positive response species. According to the health assessment, the river ecosystem in the North Canal River basin was generally unhealthy, and the upstream ecosystem was better than that downstream; the health conditions in the mountainous areas of Changping district were the best, whereas those in Chaoyang and central city districts were the worst. This study can provide a basis for ecological restoration and pollution control of rivers and also provide a reference for the water ecological civilization construction in other cities.

Taxon-specific sensitivities to flow intermittence reveal macroinvertebrates as potential bioindicators of intermittent rivers and streams.

As complex mosaics of lotic, lentic, and terrestrial habitats, intermittent rivers and ephemeral streams (IRES) support high biodiversity. Despite their ecological importance, IRES are poorly represented in routine monitoring programs, but recent recognition of their considerable-and increasing-spatiotemporal extent is motivating efforts to better represent IRES in ecological status assessments. We examine response patterns of aquatic macroinvertebrate communities and taxa to flow intermittence (FI) across three European climatic regions. We used self-organizing map (SOM) to ordinate and classify sampling sites based on community structure in regions with continental, Mediterranean and oceanic climates. The SOM passively introduced FI, quantified as the mean annual % flow, and visualized its variability across classified communities, revealing a clear association between community structure and FI in all regions. Indicator species analysis identified taxa indicative of low, intermediate and high FI. In the continental region, the amphipod Niphargus was indicative of high FI and was associated with groundwater-fed IRES, whereas indicators of Mediterranean IRES comprised Odonata, Coleoptera and Heteroptera taxa, which favor lentic conditions. In the oceanic region, taxa indicative of relatively high FI included leuctrid stoneflies and a limnephilid caddisfly, likely reflecting the colonization of IRES by aerial adults from nearby perennial reaches. The Diptera families Chironomidae and Simuliidae showed contrasting FI preferences among regions, reflecting environmental heterogeneity between regions and the coarse taxonomic resolution to which these organisms were identified. These region-specific community and taxon responses of aquatic biota to FI highlight the need to adapt standard biotic indices to enable effective ecological status assessments in IRES.

Developing Indicators of Nutrient Pollution in Streams Using 16S rRNA Gene Metabarcoding of Periphyton-Associated Bacteria.

Indicators based on nutrient-biota relationships in streams can inform water quality restoration and protection programs. Bacterial assemblages could be particularly useful indicators of nutrient effects because they are species-rich, important contributors to ecosystem processes in streams, and responsive to rapidly changing conditions. Here, we sampled 25 streams weekly (12-14 times each) and used 16S rRNA gene metabarcoding of periphyton-associated bacteria to quantify the effects of total phosphorus (TP) and total nitrogen (TN). Threshold indicator taxa analysis identified assemblage-level changes and amplicon sequence variants (ASVs) that increased or decreased with increasing TP and TN concentrations (i.e., low P, high P, low N, and high N ASVs). Boosted regression trees confirmed that relative abundances of gene sequence reads for these four indicator groups were associated with nutrient concentrations. Gradient forest analysis complemented these results by using multiple predictors and random forest models for each ASV to identify portions of TP and TN gradients at which the greatest changes in assemblage structure occurred. Synthesized statistical results showed bacterial assemblage structure began changing at 24 µg TP/L with the greatest changes occurring from 110 to 195 µg/L. Changes in the bacterial assemblages associated with TN gradually occurred from 275 to 855 µg/L. Taxonomic and phylogenetic analyses showed that low nutrient ASVs were commonly Firmicutes, Verrucomicrobiota, Flavobacteriales, and Caulobacterales, Pseudomonadales, and Rhodobacterales of Proteobacteria, whereas other groups, such as Chitinophagales of Bacteroidota, and Burkholderiales, Rhizobiales, Sphingomonadales, and Steroidobacterales of Proteobacteria comprised the high nutrient ASVs. Overall, the responses of bacterial ASV indicators in this study highlight the utility of metabarcoding periphyton-associated bacteria for quantifying biotic responses to nutrient inputs in streams.