“Pitfalls in relative abundance estimation using eDNA metabarcoding”

Metabarcoding provides unparalleled information of overall ecosystem ecology and diversity, from occurrence to relative abundances, from micro‐ to macro‐organisms, including cryptic species and larval stages that would otherwise be impossible to identify. PCR is still the most cost‐effective and routinely used method to extract information from environmental samples coupled with high‐throughput sequencing (HTS). Nichols et al. ( 2018 ) metabarcoded nine soil environmental samples from Alaska and six mixtures of synthetic oligonucleotides, using a short fragment of the chloroplast trnL gene. The authors tested the effect of PCR replicates, different polymerase mixes and amplicon GC content on diversity estimates. Results show that occurrence data was not biased for any of the parameters tested and only insufficient PCR replicates could impair the detection of rare taxa. For abundance data, polymerases were found to differentially bias estimates of the relative amplicon abundances and this depended on GC content.     

High‐throughput sequencing provides insight into manipulated soil fungal community structure and diversity during temperate forest restoration

The process of community assembly in fungal communities is poorly understood and may have important implications for restoration. However, there is a shortage of data describing fungal community composition at various stages of restoration. This study describes how microbial inoculation with field‐collected soils or a commercial inoculum influenced fungal communities during temperate tree restoration. We utilized Illumina Mi‐Seq sequencing technology to examine fungal community structure in the rhizosphere soils of trees at the conclusion of one growing season. Inoculation treatment was found to be a significant determinant of fungal community structure in one of our three experimental tree species (Liriodendron tulipifera). We also found a marginally significant influence of inoculation method on fungal community structure in the rhizosphere soils of Quercus rubra, an ectomycorrhizal tree species. Importantly, within these taxa, the use of commercial inocula, while failing to lead to detectable abundances of the inoculated taxa, strongly influenced the resulting fungal community structure after 4 months in the field, relative to control trees that received no such inoculation. We observed lower abundances of Hebeloma, a potentially important ectomycorrhizal genera, in Quercus trees receiving the commercial inoculum compared with control trees; thus, the commercial inoculum might have unexpected consequences for fungal community assembly. Such unintended legacy effects of soil inoculation should be considered in ecological restoration. Furthermore, by taking a time series approach to sampling, high‐throughput sequencing approaches could be used to test the principles of ecological assembly theory, including legacy effects of taxa no longer detectable in the community.     

Multimodality in the individual size distributions of bird communities

Aim Body size often plays a significant role in community assembly through its impacts on the life history and ecological attributes of species. Insight into the importance of size in structuring communities can be gained by examining the distribution of sizes of individuals [i.e. the individual size distribution (ISD) or size spectrum] in a community. ISDs have been studied extensively in aquatic and tree communities, but have received little attention in terrestrial animal communities. Here, we conduct the first macroecological analysis of ISDs in terrestrial animal communities to determine whether they show broad‐scale consistency in shape. Location North America, north of Mexico. Methods Using likelihood‐based methods and Gaussian mixture modelling, coupled with data from the Breeding Bird Survey and Christmas Bird Count, we determine whether the ISDs for thousands of breeding and wintering North American bird communities are: (1) monotonically decreasing, (2) unimodal or (3) multimodal. Results We find that avian ISDs are consistently multimodal, with most characterized by more than five modes in both breeding and wintering communities from local to continental scales. In addition, the positions of these modes along the size axis are remarkably consistent. Main conclusions The striking consistency in the ISD within bird communities, as with tree and aquatic communities, indicates that the ISD is an important and informative characterization of resource utilization within an ecological assemblage. The differences in shape of the ISD among these groups also suggest that differences in body size‐related constraints affect interactions within a group and with the environment. Our results confirm that avian assemblages do exhibit structure along the body size axis, and therefore it will be fruitful to explore this pattern in greater detail.     

The effect of excluding taxa with low abundances or taxa with small distribution ranges on ecological assessment

The present study aims to investigate whether taxa with a small distribution range or taxa with low abundances indicate specific habitats or a high ecological quality and what the effect is if these taxa are excluded from ecological assessment. We compared autecological features between stream dwelling taxa with a mean abundance >5 individuals per sample and a mean abundance ≤5 individuals per sample as well as between taxa with a small distribution range and taxa with a large distribution range. The number of rare taxa (either with a small distribution range or with low abundances) in a sample was related to the ecological quality classes. To test the effect of exclusion of rare taxa we constructed 8 data sets all including 142 samples of Dutch lowland streams. From each data set we stepwise excluded taxa that had low abundances or taxa that were known to be restricted in their distribution range. With help of the AQEM assessment software we calculated the final ecological quality classes and the metrics that were included in the multimetric for the original data and the 8 selected data sets. Autecological features of the taxa within the different selections showed that taxa with small distribution ranges were often running water taxa, living on stones and gravel and indicating oligosaprobic water conditions in contrast to taxa that had a large distribution range. There were only small differences between taxa with low and high abundances. However, current velocity preference was lower for taxa with abundance ≤5 individuals per sample, saprobic values were higher and scores for typical stream habitats, such as lithal, psammal and akal were lower compared to high abundant taxa. If taxa with low abundances were excluded a higher ecological quality class was achieved in most cases, while excluding taxa with a small distribution range resulted in lower ecological quality classes. In conclusion, excluding taxa with a small distribution range led to worse ecological quality classes because these taxa have special autecological features that often indicate natural streams. On the other hand, excluding taxa with low abundances resulted in higher ecological quality classes because these taxa indicate more disturbed situations and because the number of taxa per sample was strongly reduced. Although the documentation of rare taxa (either with low abundances or with small distribution ranges) is often time and cost-intensive regarding field work, laboratory work, data processing, and analyses, the indicative power of these taxa for natural circumstances is essential and therefore rare taxa should be included in ecological assessment studies.     

Spatial patterns of bacterial taxa in nature reflect ecological traits of deep branches of the 16S rRNA bacterial tree

Whether bacteria display spatial patterns of distribution and at which level of taxonomic organization such patterns can be observed are central questions in microbial ecology. Here we investigated how the total and relative abundances of eight bacterial taxa at the phylum or class level were spatially distributed in a pasture by using quantitative PCR and geostatistical modelling. The distributions of the relative abundance of most taxa varied by a factor of 2.5–6.5 and displayed strong spatial patterns at the field scale. These spatial patterns were taxon‐specific and correlated to soil properties, which indicates that members of a bacterial clade defined at high taxonomical levels shared specific ecological traits in the pasture. Ecologically meaningful assemblages of bacteria at the phylum or class level in the environment provides evidence that deep branching patterns of the 16S rRNA bacterial tree are actually mirrored in nature.     

The accuracy of absolute differential abundance analysis from relative count data

Concerns have been raised about the use of relative abundance data derived from next generation sequencing as a proxy for absolute abundances. For example, in the differential abundance setting, compositional effects in relative abundance data may give rise to spurious differences (false positives) when considered from the absolute perspective. In practice however, relative abundances are often transformed by renormalization strategies intended to compensate for these effects and the scope of the practical problem remains unclear. We used simulated data to explore the consistency of differential abundance calling on renormalized relative abundances versus absolute abundances and find that, while overall consistency is high, with a median sensitivity (true positive rates) of 0.91 and specificity (1—false positive rates) of 0.89, consistency can be much lower where there is widespread change in the abundance of features across conditions. We confirm these findings on a large number of real data sets drawn from 16S metabarcoding, expression array, bulk RNA-seq, and single-cell RNA-seq experiments, where data sets with the greatest change between experimental conditions are also those with the highest false positive rates. Finally, we evaluate the predictive utility of summary features of relative abundance data themselves. Estimates of sparsity and the prevalence of feature-level change in relative abundance data give reasonable predictions of discrepancy in differential abundance calling in simulated data and can provide useful bounds for worst-case outcomes in real data.     

Structure of small tropical island freshwater fish and crustacean communities: A niche‐or dispersal‐based process?

Determining the relative importance of niche‐ and dispersal‐based processes in the structuring of animal communities is central in ecology. Freshwater fish and crustacean communities of small tropical islands can bring new insights for understanding these processes as all their species present a pelagic larval stage which gives them important dispersal capacities. Consequently, we hypothesized that dispersal‐based process may be preponderant for structuring these communities from the regional to the local, that is, survey site, scales. Gradient forest analyses allowed us to model the effect of 12 environmental variables on freshwater taxa abundances in two southwestern Indian Ocean islands: Mayotte (26 taxa) and Reunion (21). A total of 153 surveys in Mayotte and 266 in Reunion were used for building the models. Despite the strong heterogeneity of environmental conditions between the two islands, the main factors structuring freshwater fish and crustacean communities in both islands were the elevation and the slope of the sites. The observed structure appeared more pronounced for predatory species than for primary consumers and omnivores. As predators generally have limited locomotor capacities, it is concluded that dispersal‐based process structures these communities not only at large geographical scales but also at the intra‐watershed scale, by limiting the inland penetration (or dispersion) of species in relation to their locomotor capacities. However, more knowledge concerning ecological traits and taxonomic status of many species is needed to confirm this assumption. Abstract in French is available with online material     

Distinct ecological niches of marine symbiotic N2‐fixing cyanobacterium Candidatus Atelocyanobacterium thalassa sublineages

A recently described symbiosis between the metabolically streamlined nitrogen‐fixing cyanobacterium UCYN‐A and a single‐celled eukaryote prymnesiophyte alga is widely distributed throughout tropical and subtropical marine waters, and is thought to contribute significantly to nitrogen fixation in these regions. Several UCYN‐A sublineages have been defined based on UCYN‐A nitrogenase (nifH) sequences. Due to the low abundances of UCYN‐A in the global oceans, currently existing molecular techniques are limited for detecting and quantifying these organisms. A targeted approach is needed to adequately characterize the diversity of this important marine cyanobacterium, and to advance understanding of its ecological importance. We present findings on the distribution of UCYN‐A sublineages based on high throughput sequencing of UCYN‐A nifH PCR amplicons from 78 samples distributed throughout many major oceanic provinces. These UCYN‐A nifH fragments were used to define oligotypes, alternative taxonomic units defined by nucleotide positions with high variability. The data set was dominated by a single oligotype associated with the UCYN‐A1 sublineage, consistent with previous observations of relatively high abundances in tropical and subtropical regions. However, this analysis also revealed for the first time the widespread distribution of the UCYN‐A3 sublineage in oligotrophic waters. Furthermore, distinct assemblages of UCYN‐A oligotypes were found in oligotrophic and coastally influenced waters. This unique data set provides a framework for determining the environmental controls on UCYN‐A distributions and the ecological importance of the different sublineages.     

Higher contribution of globally rare bacterial taxa reflects environmental transitions across the surface ocean

Microbial taxa range from being ubiquitous and abundant across space to extremely rare and endemic, depending on their ecophysiology and on different processes acting locally or regionally. However, little is known about how cosmopolitan or rare taxa combine to constitute communities and whether environmental variations promote changes in their relative abundances. Here we identified the Spatial Abundance Distribution (SpAD) of individual prokaryotic taxa (16S rDNA‐defined Operational Taxonomic Units, OTUs) across 108 globally‐distributed surface ocean stations. We grouped taxa based on their SpAD shape (“normal‐like”‐ abundant and ubiquitous; “logistic”‐ globally rare, present in few sites; and “bimodal”‐ abundant only in certain oceanic regions), and investigated how the abundance of these three categories relates to environmental gradients. Most surface assemblages were numerically dominated by a few cosmopolitan “normal‐like” OTUs, yet there was a gradual shift towards assemblages dominated by “logistic” taxa in specific areas with productivity and temperature differing the most from the average conditions in the sampled stations. When we performed the SpAD categorization including additional habitats (deeper layers and particles of varying sizes), the SpAD of many OTUs changed towards fewer “normal‐like” shapes, and OTUs categorized as globally rare in the surface ocean became abundant. This suggests that understanding the mechanisms behind microbial rarity and dominance requires expanding the context of study beyond local communities and single habitats. We show that marine bacterial communities comprise taxa displaying a continuum of SpADs, and that variations in their abundances can be linked to habitat transitions or barriers that delimit the distribution of community members.     

Sampling,sequencing and the SAD

The Species Abundance Distribution (SAD) is a common metric for characterizing macroscopic ecological communities. Recently, this metric has been applied to analysis of microbial communities as well. However, as compared to macroscopic communities, sampling of microscopic communities is different. In particular, most microbial communities are studied using sequencing techniques. These techniques have known biases that result in certain taxa being detected more often than others, even if the taxa are present in the sample at equivalent abundances. There are, for example, amplification biases that result in some sequences being amplified more than others. Likewise, differences in genome size across organisms can result in different numbers of reads from different taxa, again resulting in biased detection. A number of bioinformatics methods have been devised to account for biases in sequencing data, allowing for more accurate estimates of relative taxon abundances. However, because the sampling process itself is affected by biased detection, and because sampling (and under-sampling in particular) can influence the shape of the SAD, it is possible that, even when corrected for through re-scaling, detection biases can affect SAD predictions from sequencing data. To test this hypothesis, we construct a simulation model of the sampling process, focusing on biased detection in shotgun sequencing that arises from genome size differences across microbial taxa. Interestingly, we find that, although genome size itself does not impact SAD predictions, predictions can vary depending on the range of genome sizes that are represented in a community, as well as how genome size is distributed (i.e., whether the majority of species have small versus large genomes). Our results suggest that care should be taken when comparing SADs across environments, particularly when those environments might have taxa with different genome size distributions. Furthermore, our results indicate that relatively deep sequencing might be required to avoid drawing spurious inferences about ecological differences across microbial communities.     

The ecology of the planktonic diatom Cyclotella and its implications for global environmental change studies

The fossil record of diatoms in lake sediments can be used to assess the effects of climate variability on lake ecosystems if ecological relationships between diatom community structure and environmental parameters are well understood. Cyclotella sensu lato taxa are a key group of diatoms that are frequently dominant members of phytoplankton communities in low‐ to moderate‐productivity lakes. Their relative abundances have fluctuated significantly in palaeolimnological records spanning over a century in arctic, alpine, boreal and temperate lakes. This suggests that these species are sensitive to environmental change and may serve as early indicators of ecosystem effects of global change. Yet patterns of change in Cyclotella species are not synchronous or unidirectional across, or even within, regions, raising the question of how to interpret these widespread changes in diatom community structure. We suggest that the path forward in resolving seemingly disparate records is to identify clearly the autecology of Cyclotella species, notably the role of nutrients, dissolved organic carbon and light, coupled with better consideration of both the mechanisms controlling lake thermal stratification processes and the resulting effects of changing lake thermal regimes on light and nutrients. Here we begin by reviewing the literature on the resource requirements of common Cyclotella taxa, illustrating that many studies reveal the importance of light, nitrogen, phosphorus, and interactions among these resources in controlling relative abundances. We then discuss how these resource requirements can be linked to shifts in limnological processes driven by environmental change, including climate‐driven change in lakewater temperature, thermal stratification and nutrient loading, as well as acidification‐driven shifts in nutrients and water clarity. We examine three case studies, each involving two lakes from the same region that have disparate trends in the relative abundances of the same species, and illustrate how the mechanisms by which these species abundances are changing can be deciphered. Ultimately, changes in resource availability and water clarity are key factors leading to shifts in Cyclotella abundances. Tighter integration of the autecology of this important group of diatoms with environmental change and subsequent alterations in limnological processes will improve interpretations of palaeolimnological records, and clarify the drivers of seemingly disparate patterns in fossil records showing widespread and rapid changes across the northern hemisphere.     

Projected effects of Climate‐change‐induced flow alterations on stream macroinvertebrate abundances

Global change has the potential to affect river flow conditions which are fundamental determinants of physical habitats. Predictions of the effects of flow alterations on aquatic biota have mostly been assessed based on species ecological traits (e.g., current preferences), which are difficult to link to quantitative discharge data. Alternatively, we used empirically derived predictive relationships for species’ response to flow to assess the effect of flow alterations due to climate change in two contrasting central European river catchments. Predictive relationships were set up for 294 individual species based on (1) abundance data from 223 sampling sites in the Kinzig lower‐mountainous catchment and 67 sites in the Treene lowland catchment, and (2) flow conditions at these sites described by five flow metrics quantifying the duration, frequency, magnitude, timing and rate of flow events using present‐day gauging data. Species’ abundances were predicted for three periods: (1) baseline (1998–2017), (2) horizon 2050 (2046–2065) and (3) horizon 2090 (2080–2099) based on these empirical relationships and using high‐resolution modeled discharge data for the present and future climate conditions. We compared the differences in predicted abundances among periods for individual species at each site, where the percent change served as a proxy to assess the potential species responses to flow alterations. Climate change was predicted to most strongly affect the low‐flow conditions, leading to decreased abundances of species up to ?42%. Finally combining the response of all species over all metrics indicated increasing overall species assemblage responses in 98% of the studied river reaches in both projected horizons and were significantly larger in the lower‐mountainous Kinzig compared to the lowland Treene catchment. Such quantitative analyses of freshwater taxa responses to flow alterations provide valuable tools for predicting potential climate‐change impacts on species abundances and can be applied to any stressor, species, or region.     

Evaluating multiplexed next‐generation sequencing as a method in palynology for mixed pollen samples

The identification of pollen plays an important role in ecology, palaeo‐climatology, honey quality control and other areas. Currently, expert knowledge and reference collections are essential to identify pollen origin through light microscopy. Pollen identification through molecular sequencing and DNA barcoding has been proposed as an alternative approach, but the assessment of mixed pollen samples originating from multiple plant species is still a tedious and error‐prone task. Next‐generation sequencing has been proposed to avoid this hindrance. In this study we assessed mixed pollen probes through next‐generation sequencing of amplicons from the highly variable, species‐specific internal transcribed spacer 2 region of nuclear ribosomal DNA. Further, we developed a bioinformatic workflow to analyse these high‐throughput data with a newly created reference database. To evaluate the feasibility, we compared results from classical identification based on light microscopy from the same samples with our sequencing results. We assessed in total 16 mixed pollen samples, 14 originated from honeybee colonies and two from solitary bee nests. The sequencing technique resulted in higher taxon richness (deeper assignments and more identified taxa) compared to light microscopy. Abundance estimations from sequencing data were significantly correlated with counted abundances through light microscopy. Simulation analyses of taxon specificity and sensitivity indicate that 96% of taxa present in the database are correctly identifiable at the genus level and 70% at the species level. Next‐generation sequencing thus presents a useful and efficient workflow to identify pollen at the genus and species level without requiring specialised palynological expert knowledge.     

Intertidal assemblages on seawalls and vertical rocky shores in Sydney Harbour,Australia

Abstract Understanding the ecological role of artificial structures, such as seawalls, in shallow coastal waters is necessary in order to plan sound strategies of conservation and management of natural habitats. In Sydney Harbour (NSW, Australia), about 50% of the foreshore is made of retaining seawalls This study evaluates the changes caused to natural assemblages of organisms by these structures, by comparing intertidal assemblages between seawalls and vertical rocky shores. The following hypotheses were tested: that assemblages on seawalls would differ from those on rocky shores at mid‐, but not at low‐shore levels; where assemblages differ between habitats, there would be differences in cover/abundances of widespread species; patterns would be consistent among locations and through time; the variability of assemblages at the scales of 10s of cm and metres would differ between seawalls and rocky shores at mid‐ and low‐shore levels. To test these hypotheses, assemblages on seawalls and rocky shores were sampled at three locations, at roughly 4‐monthly intervals, over a period of about 18 months. Results indicated that mid‐shore assemblages on seawalls were different from those on rocky shores, but this was not the case at low‐shore levels. Few taxa were unique to either habitat. Cover of common species of algae and sessile animals and abundances of mobile grazers were variable with few consistent patterns. Variability at the scales sampled differed between habitats and heights on the shore. Seawalls and rocky shores, in general, supported a similar suite of species, but patterns of abundance and variation differed among locations and from height to height in each habitat. The implications of these findings for the future management of seawalls are briefly considered.     

Not all are free‐living: high‐throughput DNA metabarcoding reveals a diverse community of protists parasitizing soil metazoa

Protists, the most diverse eukaryotes, are largely considered to be free‐living bacterivores, but vast numbers of taxa are known to parasitize plants or animals. High‐throughput sequencing (HTS) approaches now commonly replace cultivation‐based approaches in studying soil protists, but insights into common biases associated with this method are limited to aquatic taxa and samples. We created a mock community of common free‐living soil protists (amoebae, flagellates, ciliates), extracted DNA and amplified it in the presence of metazoan DNA using 454 HTS. We aimed at evaluating whether HTS quantitatively reveals true relative abundances of soil protists and at investigating whether the expected protist community structure is altered by the co‐amplification of metazoan‐associated protist taxa. Indeed, HTS revealed fundamentally different protist communities from those expected. Ciliate sequences were highly over‐represented, while those of most amoebae and flagellates were under‐represented or totally absent. These results underpin the biases introduced by HTS that prevent reliable quantitative estimations of free‐living protist communities. Furthermore, we detected a wide range of nonadded protist taxa probably introduced along with metazoan DNA, which altered the protist community structure. Among those, 20 taxa most closely resembled parasitic, often pathogenic taxa. Therewith, we provide the first HTS data in support of classical observational studies that showed that potential protist parasites are hosted by soil metazoa. Taken together, profound differences in amplification success between protist taxa and an inevitable co‐extraction of protist taxa parasitizing soil metazoa obscure the true diversity of free‐living soil protist communities.     

COASTAL DIATOM‐ENVIRONMENT RELATIONSHIPS IN THE BRACKISH BALTIC SEA1

High‐quality calibration data sets are required when diatom assemblages are used for monitoring ecological change or reconstructing palaeo‐environments. The quality of such data sets can be validated, in addition to other criteria, by the percentage of significant unimodal species responses as a measure of the length of an environmental gradient. This study presents diatom‐environment relationships analyzed from a robust data set of diatom communities living on submerged stones along a 2,000 km long coastline in the Baltic Sea area, including 524 samples taken at 135 sites and covering a salinity gradient from 0.4 to 11.4. Altogether, 487 diatom taxa belonging to 102 genera were recorded. Detrended canonical correspondence analysis showed that salinity was the overriding environmental factor regulating diatom community composition, while exposure to wave action and nutrient concentrations were of secondary importance. Modeling the abundances of the 58 most common diatom taxa yielded significant relationships with salinity for 57 taxa. Twenty‐three taxa showing monotonic responses were species with optimum distributions in freshwater or marine waters. Thirty‐four taxa showing unimodal responses were brackish‐water species with maximum distributions at different salinities. Separate analyses for small (cell biovolume <1,000 μm³) and large (≥1,000 μm³) taxa yielded similar results. In previous studies along shorter salinity gradients, large and small epilithic diatom taxa responded differently. From our large data, we conclude that counts of large diatom taxa alone seem sufficient for indicating salinity changes in coastal environments with high precision.     

Comparing High‐throughput Platforms for Sequencing the V4 Region of SSU‐rDNA in Environmental Microbial Eukaryotic Diversity Surveys

High‐throughput sequencing platforms are continuing to increase resulting read lengths, which is allowing for a deeper and more accurate depiction of environmental microbial diversity. With the nascent Reagent Kit v3, Illumina MiSeq now has the ability to sequence the eukaryotic hyper‐variable V4 region of the SSU‐rDNA locus with paired‐end reads. Using DNA collected from soils with analyses of strictly‐ and nearly identical amplicons, here we ask how the new Illumina MiSeq data compares with what we can obtain with Roche/454 GS FLX with regard to quantity and quality, presence and absence, and abundance perspectives. We show that there is an easy qualitative transition from the Roche/454 to the Illumina MiSeq platforms. The ease of this transition is more nuanced quantitatively for low‐abundant amplicons, although estimates of abundances are known to also vary within platforms.     

Remembrance of things past: modelling the relationship between species' abundances in living communities and death assemblages

Accumulations of dead skeletal material are a valuable archive of past ecological conditions. However, such assemblages are not equivalent to living communities because they mix the remains of multiple generations and are altered by post-mortem processes. The abundance of a species in a death assemblage can be quantitatively modelled by successively integrating the product of an influx time series and a post-mortem loss function (a decay function with a constant half-life). In such a model, temporal mixing increases expected absolute dead abundance relative to average influx as a linear function of half-life and increases variation in absolute dead abundance values as a square-root function of half-life. Because typical abundance distributions of ecological communities are logarithmically distributed, species' differences in preservational half-life would have to be very large to substantially alter species' abundance ranks (i.e. make rare species common or vice-versa). In addition, expected dead abundances increase at a faster rate than their range of variation with increased time averaging, predicting greater consistency in the relative abundance structure of death assemblages than their parent living community.     

Climate‐driven substitution of habitat‐forming species leads to reduced biodiversity within a temperate marine community

Aim

In marine ecosystems, habitat‐forming species (HFS) such as reef‐building corals and canopy‐forming macroalgae alter local environmental conditions and can promote biodiversity by providing biogenic living space for a vast array of associated organisms. We examined community‐level impacts of observed climate‐driven shifts in the relative abundances of two superficially similar HFS, the warm‐water kelp Laminaria ochroleuca and the cool‐water kelp Laminaria hyperborea.

Location

Western English Channel, north‐east Atlantic

Methods

We compared algal and invertebrate assemblages associated with kelp stipes and holdfasts, across multiple sites and sampling events. Significant differences were recorded in the structure of assemblages between the host kelp species at each site and event.

Results

Assemblages associated with stipes of the cool‐water HFS were, on average, >12 times more diverse and supported >3600 times more biomass compared with the warm‐water HFS. Holdfast assemblages also differed significantly between species, although to a lesser extent than those associated with stipes. Overall, assemblages associated with the warm‐water HFS were markedly impoverished and comprised far fewer rare or unique taxa.

Main conclusions

While previous research has shown how climate‐driven loss of HFS can cause biodiversity loss, our study demonstrates that climate‐driven substitutions of HFS can also lead to impoverished assemblages. The indirect effects of climate change remain poorly resolved, but shifts in the distributions and abundances of HFS may invoke widespread ecological change, especially in marine ecosystems where facilitative interactions are particularly strong.