Effect of Didymosphenia geminata coverage on the phytobenthic community in an Andean basin of Chile

Abstract

Background

The invasive diatom Didymosphenia geminata (Lyngbye) M. Schmidt has received considerable attention in recent years due to its rapid geographic expansion and massive proliferation, which have altered habitat availability for benthic species and triggered negative effects on stream ecosystems. We evaluated the changes in the community structure of phytobenthos caused by the temporal variation in D. geminata coverage, in addition to the environmental variables correlated with the temporal variation of this invasive microalga in the Andean sector of the Biobio River basin, Chile.

Methods

Environmental parameters were measured during the austral summer of 2014–2015, when phytobenthos samples were collected and used to develop a relative abundance matrix of taxa by calculating species richness and Shannon diversity. Multivariate techniques were used to establish the relationships among environmental variables, including D. geminata coverage, and the phytobenthic community.

Results

Massive proliferation of D. geminata occurs during summer (December–January). According to multiple regression analysis, electrical conductivity, temperature and total phosphorus were the variables that best explained the variation in D. geminata coverage. When D. geminata coverage was over 50%, phytobenthic species richness was significantly higher than at the uninvaded site, without a significant change in Shannon diversity. In addition, the % coverage of this invasive microalga and total phosphorus concentration were variables that differentiated phytobenthic communities among the study sites.

Conclusions

Environmental factors such as conductivity, temperature and total phosphorus concentration influenced the temporal variability of D. geminata mats. In addition, the massive growth of this invasive diatom caused a higher species richness without altering Shannon diversity. Our results suggest that the spatio-temporal variability of D. geminata correlated with environmental variables will help predict the habitat suitability of this alga in other Andean rivers and allow a better understanding of ecological habitat alterations.

Background

Invasive organisms are species that can settle outside of their natural distribution areas, triggering negative ecological and economic consequences [1]. During the last decade, the benthic diatom D. geminata (Lyngbye) M. Schmidt has been recognized as an invasive species that is capable of displacing benthic species [2, 3] given its significant proliferation and rapid dispersion in freshwater ecosystems [4]. This diatom has an apical pore through which it secretes a mucopolysaccharide stalk, allowing it to attach to various substrates [5], forming large macroscopic conglomerates that are resistant to bacterial degradation and may remain in the environment for weeks and months, even after cell death [6]. D. geminata is native to Europe, New Zealand and North America and inhabits lakes and rivers with low nutrient concentrations [7], although in recent decades it has been observed that it is expanding its geographical range, forming massive growths in rivers and lakes at high latitudes in South America [8,9,10].

Observations have shown that D. geminata has a high tolerance to changes in water quality [11,12,13] and its massive proliferation is mainly associated with low phosphorus concentrations [6, 14,15,16], low current velocity and depth [17,18,19]. In addition, its ability to completely cover various riverbed substrates, make this benthic diatom not only a threat to other aquatic organisms, but also one of the most damaging invasive species in river ecosystems [7, 16, 20, 21], as it alters the availability of habitat for benthic species, causing changes in the food chain [2, 22]. Much research in this area has focused on assessing the environmental variables that favor the settlement and proliferation of D. geminata [15, 19], but few studies have evaluated the effect of D. geminata mats on phytobenthic communities [23, 24].

In South America, studies on the proliferation of D. geminata have focused on taxonomic and phylogenetic aspects [5, 25] and the possibility of this microalga spreading to new water bodies [8, 9]. Bhatt et al. (2008) suggested the need for studies focused on the monitoring of D. geminata population growth and possible interactions with other biotic communities. In this regard, changes to the taxonomic composition of phytobenthic communities can act as a roadmap to estimate the alteration of water bodies that have undergone ecological impacts [26]. Thus, communities, apart from providing more knowledge about microalgal flora [27], could also be considered an indicator of the impact of the massive proliferation of D. geminata.

The main objective of our study was to determine the effect of the variation in D. geminata coverage on the structure of phytobenthic communities. We also analyzed the environmental variables that correlate with the temporal variation of this invasive microalga in the Andean sector of the Biobio River basin, Chile.

Materials and methods

Study area

The Biobío River basin has a drainage area of 24,260 km² and is considered a multiple-use river system (e.g., hydroelectricity generation, domestic and industrial effluent receiving, drinking water supply and fishing) [28]. Fieldwork was done between December 2014 and March 2015, in the austral summer. Sampling sites were selected in the main course of the Biobio River. Three of the analyzed sites were invaded by D. geminata (Marimenuco, MAR; Tallón, TAL; Lolén, LOL); a tributary free of Didymo blooms was included as non-invaded site (Liucura, LIU) (Fig. 1). The topographic composition of the study area results in highly variable and seasonal climate patterns, with cold-wet winters and hot-dry summers [29]. Average daily precipitation during the study period varied between 0 and 0.66 mm, while average daily air temperature fluctuated between 7.11 and 20.63 °C (Dirección General de Aguas, DGA). The vegetation cover of this area consists of an Araucaria forest (Araucaria araucana, (Molina) K. Koch) with understory of Antarctic beech (Nothofagus antarctica, (G. Forster) Oerst). The species Festuca scabriuscula (Phil.), Acaena sericea (Phil.), Baccharis magellanica (Lam.) and Rumex acetosella (L.) dominate the herbaceous stratum [29, 30].

Fig. 1

Map of Andean section of the Biobio river basin (Chile). White circles indicate sites invaded by D. geminata and the black circle the site that was not invaded. The arrows indicate the flow direction

River habitat characteristics

At each study site, a section of approximately 50 m in length was selected considering the following criteria: similar substrate size, current velocity, depth and open canopy. Water quality variables were established with in situ measurements of temperature (°C), conductivity (μS/cm), pH and dissolved oxygen (mg/L), using a Quanta multiparameter probe. Total nitrogen (mg/L), total phosphorus (μg/L) and soluble phosphorus (μg/L) were analyzed in the laboratory following APHA methodology [31]. Hydrological variables included flow velocity (m/s) and depth (m), which were measured at substrate level using an FP111 Global Flow Probe. At each site, a total of six replicates for each variable measured in situ and three replicates for nutrients were considered. The River Habitat Index (IHF) described by Pardo et al. [32] was used to characterize the river habitat at each site. This index includes seven categories, which assess hydrogeomorphological aspects of the channel, including the presence and frequency of rapids, substrate composition, velocity and depth regimes, channel exposure, habitat heterogeneity and aquatic vegetation coverage.

D. geminata coverage and phytobenthos sampling

The presence of D. geminata was characterized according to the criteria described by [33], which are composed of five categories: Absent, no mucilaginous mats observed. Initial growth, less than 20% coverage, small brownish mats were observed. Medium growth, visual coverage ≥20% and no more than 50%, some coalesced patches completely cover some rocks. High growth, coverage is ≥50% and less than 80%, with the mats increasing in thickness and in dry conditions taking on the appearance of wet paper. Finally, growth was deemed massive when the thickness of the mats increased such that the bedrock and riparian vegetation were in large part covered (coverage > 80%).

To determine D. geminata coverage at each invaded site, three similarly-sized submerged cobbles (greater than 12 cm) were randomly chosen [33]. Then, a 16 cm² area was defined using a plastic grid and photographed to determine the percent coverage using the IMAGE J program. This % coverage of each cobble was included as an environmental variable in the data analysis.

Subsequently, epilithic samples were taken from the cobbles by scraping their upper surfaces on the aforementioned area and brushing and washing them in a plastic bowl until the selected area was completely clear. The collected material was deposited in 15-mL plastic bottles and preserved in Lugol’s solution at a temperature of 4 ± 2 °C [25]. In the laboratory, samples were analyzed and up to 400 individuals were identified in each sample with a light microscope at a magnification of 1000 [34] using an Olympus CX31 equipped with a digital camera for microphotography. For specific diatom identification the samples were cleaned and mounted [35]. Taxonomic identification was carried out based on specialized literature [36,37,38,39,40]. After counting the cells, a taxa abundance matrix was created, expressed in relative abundance (%). In addition, the community parameters – species richness (number of species) and the Shannon diversity index for each site and month of sampling (between December 2014 and March 2015) – were calculated.

Statistical analysis

We used multiple linear regressions to identify environmental variables that best explained the temporal variability of D. geminata. Prior to this analysis, the independent variables that did not meet the assumptions of normality were transformed to ln (x + 1). This analysis only considered sites invaded by D. geminata. A monthly comparison of the variables selected in the linear model, between the non-invaded site and each impacted site, was performed using nonparametric Mann-Whitney U comparisons. Differences in species richness and Shannon-Weaver diversity between the non-invaded and invaded sites were tested using nonparametric Mann-Whitney U comparisons. This analysis was conducted in R [41].

Principal coordinate analysis (PCoA) was performed to estimate similarity between the phytobenthos communities at different sampling sites. We used two-way permutational multivariate analysis of variance (PERMANOVA, [42]) to check whether there were differences in the phytobenthic community composition between the invaded and non-invaded sample sites and coverage categories. PCoA analysis methods were conducted using a Chord distance matrix [43] based on relative abundance of the phytobenthos community transformed to ln (x + 1). In addition, a SIMPER analysis was performed to establish the similarity percentage in order to estimate the contribution of different phytobenthos species in the previously established groups (SIMPER, [44]). Finally, a distance-based linear model (DISTLM) was executed in order to determine the relationship between environmental variables and the phytobenthic community. Correlated environmental variables were eliminated from the analysis (r > 0.70). DISTLM analysis was carried out using forward selection under the Akaike information criterion (AIC) and coefficient of determination R² to obtain a more detailed explanation of the behavior of the variables [45]. Significant differences were assumed when the p value was ≤0.05. PCoA, PERMANOVA and DISTLM analysis methods were performed with PRIMER 6.0 [46].

Results

River habitat characterization and environmental variables

D. geminata mats covered the majority of available substrates, particularly pebbles, stones, gravel, exposed roots, trunks and branches. Coverage of the invaded sites ranged between < 20% (initial growth) and > 80% (massive growth). The greatest D. geminata coverage was recorded at MAR during December–January, with values > 80% (massive growth). Conversely, the lowest coverage values were recorded at TAL, peaking at coverage ≥20% and no more than 50% (medium growth) during December (Table 1).

Table 1 Mean values and standard deviation of the most significant environmental variables according to the multivariate linear model. Values of the River habitat index (IHF), D. geminata coverage, coverage categories described by Diaz et al. (2013) and the significant Mann-Whitney U test (***: p < 0.001, **: p < 0.01; *: p < 0.05; ns: non-significance)

The river habitat characterization of the non-invaded site presented IHF values that fluctuated between 74 and 81, indicating high habitat heterogeneity with low anthropogenic intervention and high colonization potential. Meanwhile, at the invaded sites IHF fluctuated between 58 and 69, indicating significant modifications in the stream habitat. The MAR, TAL and LOL sites presented lower values compared to LIU. The low values reported in the index are associated with a reduction in components of rapid frequency, velocity/depth regime and aquatic vegetation composition. The lowest IHF values were recorded during the months of greatest D. geminata coverage (December, January and February) (Table 1).

Multiple linear regression analysis suggested that conductivity, temperature and total phosphorus were the variables that best represented the variation in D. geminata coverage (R² = 0.80; adjusted R² = 0.72; p < 0.01). Throughout the study period, significant fluctuations in these variables between the invaded sites and non-invaded site were evident (Table 1). Conductivity presented lower values at MAR and higher values at downstream sites TAL and LOL. The temperature was higher at the invaded sites (MAR, TAL and LOL) than at the not-invaded site. In contrast, water depth presented minimal fluctuations, with no differences among the months. Finally, the three invaded sites presented lower total phosphorus concentrations compared to LIU. Paired comparisons between non-invaded and invaded sites during the study period showed that the major significant differences (p < 0.05) occurred in conductivity, temperature and total phosphorus (Table 1).

Phytobenthos community

A total of 51 taxa were observed in the phytobenthic community structure during the sampling period (December–March), the majority of which belong to the class Bacillarophyceae. When D. geminata coverage was greater than 50%, the specific richness increased significantly with respect to the non-invaded site (p < 0.05), while when the coverage was less than 30%, the richness was similar to that of the non-invaded site (Fig. 2). D. geminata was the taxon with the highest relative abundance at the MAR, TAL and LOL sites; however, it was observed that the abundance of the diatoms Rhopalodia gibba and Melosira varians present at the non-invaded site (LIU) also increased at the invaded sites. Meanwhile, the Shannon index showed no significant differences in species diversity between the non-invaded site and invaded sites.

Fig. 2

a Richness of phytobenthic taxa at the different study sites. b D. geminata coverage (%) at the different study sites. The dashed line indicates 30% D. geminata coverage and continuous line 50% coverage

The PCoA analysis based on the phytobenthos community differentiated the sites into two groups: one for the site that was not invaded by D. geminata (LIU) and the other for the invaded sites, explaining a total variation of 43.4% (Fig. 3). A two-way PERMANOVA test presented a significant effect of the sites (Pseudo-F = 3.29, p = 0.0021) and D. geminata coverage categories (Pseudo-F = 2.46, p = 0.004) on the phytobenthos community. There was also significant interaction between both factors (Pseudo-F = 1.83, p = 0.047) (Table 2).

Fig. 3

Principal Coordinate Analysis (PCoA) based on relative abundance of the phytobenthos community for the studied months and sites (white circle = site not invaded by D. geminata; black circle = sites invaded by D. geminata)

Table 2 Details of the two-way PERMANOVA test conducted on ln(x + 1)-transformed phytobenthos community abundance data recorded at the sampling sites

The relationship between environmental variables and the phytobenthos community determined by the DISTLM analysis demonstrated that under the Akaike criteria, D. geminata coverage (AIC = 56.33; Pseudo-F = 3.53; p = 0.001) and total phosphorus (AIC = 55.71; Pseudo-F = 2.32; p = 0.013) were the only variables included in the final model. Meanwhile, the coefficient of determination R² included the total number of variables in the final model (R² = 0.74); the only significant variables in the study site distribution were coverage (R² = 0.20; Pseudo-F = 3.53; p = 0.001) and total phosphorus (R² = 0.32; Pseudo-F = 2.32; p = 0.014). With both criteria, D. geminata coverage and total phosphorus concentration are variables that differentiate the phytobenthos communities among the sites.

The SIMPER analysis identified the species that contributed to differentiation between the non-invaded and invaded sites, with 81.17% dissimilarity between groups (Table 3). The top three dominant species at the invaded sites that most contributed to this dissimilarity percentage were Didymopshenia geminata (13.91%), Aulacoseira granulata (Ehrenberg) Simonsen (4.99%) and Rhopalodia gibba (3.66). At the non-invaded site (LIU) the species that most contributed to this dissimilarity were Fragilaria sp. (4.61%), Cocconeis placentula Ehrenberg (4.52%) and Gomphoneis minuta (Stone) Kociolek & Stoermer (3.71%). These results suggest that the change in species abundance was responsible for the high level of dissimilarity between sites.

Table 3 Similarity percentage analysis (SIMPER). Average abundance values are shown for phytobenthos species at sites non-invaded (LIU) and invaded (MAR, TAL, and LOL) by D. geminata. The comparison between both sites is shown, along with average dissimilarity values and contribution percentages (%). Bacillariophyceae (B), Cyanophyta (C), Chlorophyta (Ch)

The taxa that were positively related to the first dbRDA axis were those that were present at LIU, but absent or present at very low relative abundances at MAR, TAL and LOL. These include the species Gomphoneis minuta, Cymbella sp., Fragilaria sp., Cymbella affinis Kützing, Cocconeis placentula, Mougeotia sp. and Melosira hustedtii Krasske (Table 3, Fig. 4). Associated with massive D. geminata growth, an increase in relative abundance was observed for Rhopalodia gibba (Ehrenberg) Otto Müller, Melosira varians Agardh and Dolichospermum sp. Furthermore, Aulacoseira granulata, Fragilaria acus (Kützing) Lange-Bertalot and Scenedesmus quadricauda (Turpin) Brébisson were present only in invaded areas.

Fig. 4

Distance-based redundancy analysis of phytobenthos community and environmental variables. dbRDA analysis resulting from DISTLM for sampling sites (white circle = site without D. geminata, black circle = sites with D. geminata). The environmental variables shown are those determined by the coefficient of determination R², with the most significant highlighted. The lower edge shows the taxa that contribute to 60% of the dissimilarity given by SIMPER between the sites that were invaded and not invaded by this microalga

Discussion

In the last decade, massive proliferation of D. geminata has been observed in geographical areas where it had not been previously reported. This migration is probably due to human activity, which is the main vector for microalga propagation. However, [47] recently proposed that the expansion of the distribution range of D. geminata blooms is due to environmental changes. For example, the first cell registry of D. geminata in Chile dates back to the 1960s [48], although the first blooms in several rivers in southern Chile were reported 50 years later [49]. These observations show that D. geminata was not recently introduced to this region and that its cellular presence does not always cause massive blooms [50]. Thus, policies and management should focus on understanding the environmental factors that promote D. geminata blooms [47] and consider the variables related to the habitat window for the growth and persistence of D. geminata in river ecosystems. Our results indicate that the environmental parameters analyzed in this investigation showed little variability in the study period, in contrast to the great variability observed in D. geminata coverage. Thus, based on these results and in accordance with [17], the spatial variability in the coverage of D. geminata mats is explained by the habitat characteristics, and it is necessary to acquire data for long time series (annual, inter-annual) to highlight the importance of the environmental and hydrological variables that drive these changes. In our study, we emphasize that the application of the multivariate linear model indicated that conductivity, temperature and total phosphorus concentration were highly correlated with coverage of this microalga. Our results are consistent with previous studies [15, 51], which pointed out that such variables condition the establishment and temporal variation of D. geminata in rivers in Canada and the United States. Light availability was not measured in this study, but the reported low depths appeared to be directly related to light availability and the temperature in the water column, and, thus, to directly influence D. geminata growth. In this context, the combined effects of these factors (light-temperature) are determinants in the photosynthetic production of the invasive diatom D. geminata, favoring its success as a colonizer [52]. In the present work it was observed that during the months with the highest D. geminata coverage the temperature was between 15 and 21 °C, which is similar to the range found in Spanish rivers [53].

Low concentrations of total phosphorus during the study period were related to the % coverage of D. geminata mats. Similarly, [15] showed the existence of a negative relationship between the concentration of total phosphorus and D. geminata mat coverage. Low phosphorus concentrations promote cell division of this microalga and production of stalks that raise the cells above the substrate for better nutrient absorption, creating large mats that can remain in the environment for a long time [54, 55]. In contrast, during March at MAR an inverse relationship, with declining D. geminata coverage and increasing total phosphorus concentration, was observed, implying that other periphyton species have a competitive advantage over D. geminata [19]. Nonetheless, the low total phosphorus concentration recorded during March did not cause an increase in coverage. This was due to an increase in the flow rate as a result of the first rain of the month, causing D. geminata shedding. Previous reports [56, 57] indicated that alterations in the flow regime are conditions for the establishment and continued presence of this invasive microalga. During the study period, the invaded sites had a decrease in the index compared to the non-invaded site (LIU) as a result of the decreasing values of the substrate composition and the velocity regime. Such physical alteration is associated with the ability of D. geminata to colonize all available benthic substrates [57].

The result of total substrate occupation by D. geminata generated a shift in phytobenthos community structure (Table 3, Fig. 4). A negative effect observed in this study was the decrease in relative abundance of the species G. minuta, C. affinis, and Cymbella sp. at the invaded sites (MAR, TAL, LOL) compared with the non-invaded site (LIU). This change could be attributed to the competitive pressure exerted by D. geminata on these phytobenthic species with similar ecological requirements, which were ultimately displaced by the massive proliferation of D. geminata [16, 58]. Indeed, the total occupation of the substrate by this invasive microalga reduces the nutrient availability for the different strata of phytobenthos, which remain isolated after being covered by the filamentous mat formed by D. geminata. This is consistent with what was described by [59], who indicated that nutrient availability is lower at the base of the periphytic community than on its surface, implying that some species can monopolize space and other resources according to their level of morphological organization [60]. Similar conclusions were reached by [61], who suggested that the thickness of the periphytic layer affects nutrient availability, especially the lower layers, due to the speed of nutrient diffusion. In addition, the decrease in the abundance of G. minuta and Cymbella sp., dominant species in late successional stages, increases the availability of substrate, which makes a site more susceptible to the establishment of D. geminata [58]. Conversely, the presence of A. granulata, S. quadricauda, and F. acus only at sites invaded by D. geminata is due to the fact that these species come from Galletué Lake; when moving downstream, probably during periods of low river flow, they settle on the substrate, occupying the microhabitat formed by D. geminata stalks. Our results are similar to those reported by [16], who described a positive relationship between the biomass of D. geminata and the relative abundance of other small diatoms. In addition, we observed the presence of a species of Chlorophyta after the growth of D. geminata, which is consistent with the findings of [53].

Although in this study no significant differences in species diversity were found (with respect to the non-invaded site), others have reported that the new microhabitats created by D. geminata filaments decrease or increase species diversity [16, 23, 53].

Different studies have shown the effect of mats formed by D. geminata on the basis of the food web and higher trophic levels [2, 24], making it necessary to identify a threshold (% of D. geminata coverage) that allows a determination of when the taxa change is generated. An interesting observation in this study is that richness is significantly higher when the coverage percentage of D. geminata is greater than 50%. Subsequently, this richness decreases with the progressive reduction in the coverage of D. geminata, reaching values similar to those found at the non-invaded site. This alteration in the phytobenthic community was similar to that reported by [53], who found that the massive growth of D. geminata changes the taxonomic composition of phytobenthos, directly affecting its species richness.

Conclusions

We confirmed that conductivity, temperature and total phosphorus had a significant influence on variation in % coverage of D. geminata. In addition, both the % coverage of this invasive microalga and total phosphorus concentration were factors that influenced phytobenthos community structure, increasing phytobenthic species richness without altering taxa diversity. Meanwhile, a decrease in the abundance of other benthic diatoms at the invaded sites generated greater substrate availability, making them more susceptible to the establishment of D.geminata and resulting in an increase in the coverage of this invasive microalga. Nevertheless, the lack of studies regarding the presence of D. geminata in phytobenthic communities makes it difficult to understand the competition for available resources. Therefore, further studies that include competitive availability for resources, as well as ecophysiological and molecular aspects, should be carried out in order to better understand the interaction between this invasive microalga and other aquatic communities.