Fpls 14 1193668
Fpls 14 1193668
Fpls 14 1193668
*CORRESPONDENCE Ruohan Huang 1†, Wan Liu 1†, Jinghua Su 2, Shihao Li 1,3,
Wei Zhang
[email protected] Liqing Wang 1, Erik Jeppesen 4,5,6,7 and Wei Zhang 1*
† 1
Key laboratory of Exploration and Utilization of Aquatic Genetic Resources of the Ministry of
These authors have contributed equally to
this work Education, Engineering Research Center of Environmental DNA and Ecological Water Health
Assessment, Shanghai Ocean University, Shanghai, China, 2 Research Institute of Natural Ecology
RECEIVED 25 March 2023 Conservation, Shanghai Academy of Environmental Sciences, Shanghai, China, 3 Shanghai Aquatic
ACCEPTED 19 June 2023 Technology Co., Ltd, Shanghai, China, 4 Department of Ecoscience, Aarhus University,
PUBLISHED 05 July 2023 Aarhus, Denmark, 5 Sino-Danish Centre for Education and Research (SDC), University of Chinese
Academy of Sciences, Beijing, China, 6 Limnology Laboratory and EKOSAM, Department of Biological
CITATION
Sciences, Middle East Technical University, Ankara, Türkiye, 7 Institute of Marine Sciences, Middle East
Huang R, Liu W, Su J, Li S, Wang L,
Technical University, Mersin, Türkiye
Jeppesen E and Zhang W (2023) Keystone
microalgae species determine the removal
efficiency of sulfamethoxazole: a case
study of Chlorella pyrenoidosa and
microalgae consortia.
In recent years, antibiotics pollution has caused serious harm to the aquatic
Front. Plant Sci. 14:1193668.
doi: 10.3389/fpls.2023.1193668 environment, and microalgae mediated degradation of antibiotics has attracted
COPYRIGHT
increasing attention. However, the potential toxicity of antibiotics to keystone
© 2023 Huang, Liu, Su, Li, Wang, Jeppesen microalgae species or their microalgae consortia, and the impact of microalgal
and Zhang. This is an open-access article
diversity on antibiotic removal need to be further studied. In this study, we
distributed under the terms of the Creative
Commons Attribution License (CC BY). The investigated the removal efficiency and tolerance of five freshwater microalgae
use, distribution or reproduction in other (Chlorella pyrenoidosa, Scenedesmus quadricauda, Dictyosphaerium sp.,
forums is permitted, provided the original
author(s) and the copyright owner(s) are
Haematoccocus pluvialis, and Botryococcus braunii) and their microalgae
credited and that the original publication in consortia to sulfamethoxazole (SMX). We found that the removal efficiency of
this journal is cited, in accordance with SMX by C. pyrenoidosa reached 49%, while the other four microalgae ranged
accepted academic practice. No use,
distribution or reproduction is permitted between 9% and 16%. In addition, C. pyrenoidosa, S. quadricauda, and
which does not comply with these terms. Dictyosphaerium sp. had better tolerance to SMX than H. pluvialis, and their
growth and photosynthesis were less affected. At 10 and 50 mg/L SMX, the
removal capacity of SMX by mixed microalgae consortia was lower than that of C.
pyrenoidos except for the consortium with C. pyrenoidos and S. quadricauda.
The consortia generally showed higher sensitivity towards SMX than the
individual species, and the biochemical characteristics (photosynthetic
pigment, chlorophyll fluorescence parameters, superoxide anion (O 2 - ),
superoxide dismutase activity (SOD), malondialdehyde (MDA) and extracellular
enzymes) were significantly influenced by SMX stress. Therefore, the removal of
antibiotics by microalgae consortia did not increase with the number
of microalgae species. Our study provides a new perspective for the selection
of microalgal consortia to degrade antibiotics.
KEYWORDS
1 Introduction not enhance but potentially reduce the removal efficiency of SMX.
Our results contribute to a better understanding of the ecotoxicity
The relationship between biodiversity and ecosystem functioning of SMX and its removal by microalgae and microalgae consortia.
is a subject of widespread concern and has been studied in both
terrestrial and aquatic ecosystems (Diaz and Cabido, 2001). Increased
species diversity typically leads to higher functional differences in the 2 Article types
community, which may result in more effective access to and use of
resources in turn (Loreau and Hector, 2001; Fargione et al., 2007). Original Research.
Several studies have shown that high species diversity of plants,
animals, and microorganisms can resist the negative effects of a
changing external environment (Ghorbannezhad et al., 2018; 3 Materials and methods
Norberg, 2000; Vacek et al., 2021). In addition, species diversity
may have effect on the degradation of pollutants in the 3.1 Chemicals
environment. For example, inoculating microbial communities
isolated from in-situ contaminated soil can effectively improve the SMX was purchased from Sangon Biotech Co., Ltd (Shanghai,
degradation of polycyclic aromatic hydrocarbons (Li et al., 2008). And China) with >98% purity. We used 0.1% mol/L NaOH as a solvent
microalgae consortia have proven to significantly improve the to fully dissolve 0.1 g SMX. Then, we added 0.1 mol/L HCl to adjust
degradation of cefradine, p-chlorophenol, and other pollutants in the initial pH 7 of the test medium (Chen et al., 2020b). The
the water environment in comparison with monocultures (Lima et al., chemicals were used of analytical grade.
2004; Xiao et al., 2021). However, Choudhury et al. (2018) and Brisson
et al. (2020) pointed out that mixed cultivation of plants is not
necessarily better than single cultivation in resource utilization and 3.2 Algal cultures
pollutant removal.
Antibiotics, widely used in aquaculture andanimal husbandry, Five microalgae: Chlorella pyrenoidosa (FACHB-10),
have become an important class of contaminants of emerging Scenedesmus quadricauda (FACHB-44), Dictyosphaerium sp.
concern (Kümmerer, 2009; Liu et al., 2018; Li et al., 2020). (FACHB-2072), Haematoccocus pluvialis (FACHB-1928) and
Previous studies have shown that the removal of antibiotics Botryococcus braunii (FACHB-357) were purchased from the
depends on physical adsorption, advanced oxidation technology Freshwater Algae Culture Collection of the Institute of
and microbial degradation (Dutta and Mala, 2020; Huang et al., Hydrobiology (FACHB-Collection), Wuhan City, China. The
2021). Recently, microalgae-mediated biodegradation has attracted microalgae were cultivated in 3000 mL Erlenmeyer flasks
attention and proven to be an environmentally effective and safe containing BG-11 medium and incubated in a homoeothermic
antibiotic removal technology (Xiong et al., 2021; Yu et al., 2022). incubator at 25 ± 1°C under 3000 lux illumination with a light-
Microalgae such as Chlorella sp., Dictyostelium sp., and dark period of 12:12 h. To maintain the logarithmic growth of the
Scenedesmus obliquus have been shown to effectively degrade microalgae and reduce any effect caused by minor differences in
antibiotics (Chen et al., 2020a; Chen et al., 2020b; Yang et al., photo irradiance, the flasks were arranged randomly and gently
2020; Aydin et al., 2022). Nevertheless, microalgae have specific shaken three times a day (Wan et al., 2015). All experimental
preferences for antibiotic types, and different microalgae have devices used for the algal culture/for algal cultivation were
different cell hydrophobicity, which affect the degradation of autoclaved at 121°C for 30 min before use.
antibiotics (Bai and Acharya, 2016; Xiong et al., 2017c).
Therefore, recent studies have focused on the consortia of
microalgae and bacteria or other microalgae to improve the 3.3 Experimental design
removal efficiency of antibiotics (Rodrigues et al., 2020; Rodrigues
et al., 2021; Wang et al., 2022). However, whether the determinant To explore differences in toxicity to and removal efficiency of
of antibiotic removal efficiency of microalgae consortia is higher antibiotics by keystone microalgae species in monocultures or when
than for keystone species needs further studies. combined, we designed two experiments: microalgae screening and
We investigated the removal and tolerance capacity to microalgae consortia.
sulfamethoxazole (SMX) of five green microalgae (C. pyrenoidosa,
S. quadricaudas, Dictyostelium sp., H. pluvialis and Botryococcus 3.3.1 Microalgae screening
braunii) and then screened several combinations of these According to previous research and our pre-experiment results
microalgae for consortia effects. To explore the difference in the our preliminary experimental results on nine microalgae species
efficiency of removing antibiotics by different microalgae consortia, (Supplementary Figure 1), we selected five microalgae, C.
we measured a series of biochemical indicators. We hypothesized pyrenoidosa, S. quadricauda, Dictyosphaerium sp., H. pluvialis and
that the keystone microalgae species had a strong removal capacity B. braunii, which have shown tolerance and removal capacity to
for SMX. Besides, consortia of individual keystone microalgae with organic pollutants in the environment (Cheirsilp et al., 2016; Chen
a high removal efficiency of SMX and other microalgae species may et al., 2020a; Yang et al., 2020; Aydin et al., 2022), as experimental
objects and explored their tolerance to and removal capacity of sample was centrifuged at 6500 rpm for 15 min. After removing the
SMX. The SMX solution was prepared with BG-11 medium with supernatant, 5 mL of 95% methanol was added to the sample, which
the experimental concentration of 5 mg/L and added to the five was then placed in a 60 °C constant temperature water bath for
microalgae cultures. The experiment was then conducted in 500-mL 10 min. The mixed samples were centrifuged at 6500 rpm for
Erlenmeyer flasks containing 400 mL of BG-11 medium inoculated 15 min, and absorbance of the supernatant at 665, 652, and 470 nm
with 15 mg/L microalgal cell suspension. Besides, the control was measured with a spectrophotometer. Chlorophyll a, chlorophyll
comprised an algal medium without added SMX, and an abiotic b and carotenoids were calculated as follows:
control group contained only SMX ([CK]). Three replicates were set
up for each treatment and control group, and all flasks were Chlorophyll a (Ca, mg=L) = 16:82 A665 − 9:28 A652
incubated under the same conditions used for the inoculum
culture for 16 days. Samples were collected on days 2, 4, 6, 8, 10, Chlorophyll b (Cb, mg=L) = 36:92 A652 − 16:54 A665
12, 14, and 16 to determine each index.
Carotenoids (mg=L) = (1000 A470 − 1: 91Ca − 95: 15Cb)=225
3.3.2 Microalgae consortia Moreover, chlorophyll fluorescence parameters, including
Based on the results of our first experiment, we selected three of photosynthetic system II (PSII) maximum photosynthetic rate (Fv/
the five species: C. pyrenoidosa , S. quadricauda and Fm) and actual photosynthetic efficiency (Yield), were measured
Dictyosphaerium sp. C. pyrenoidosa, which had the strongest using a pulse amplitude-modulated fluorometer (Phyto-PAM,
tolerance and removal capacity to SMX, was selected as control Walz, Effeltrich, Germany) equipped with an emitter-detector-
(C), and groups where C. pyrenoidosa was combined with the other fiberoptic unit with an irradiance of 16 mmol photons/m2/sec PAR.
two microalgae respectively acted as treatments (Table 1). All
microalgae in the treatment group were inoculated at the same
fresh weight ratio, and the total content of microalgae in the final 3.6 Antioxidant responses
system was 15mg/L. Microalgae consortia experiments were carried
out with 10 and 50 mg/L SMX, respectively. According to the Samples of microalgae suspension (5 mL) were collected on day
microalgae screening experimental results, SMX can be mostly 3 and 7 to determine superoxide anion. The samples were
removed by microalgae in the first 7 days. Therefore, a 7-days centrifuged at 7000 rpm for 20 min, and then the supernatant
exposure experiment was conducted where photosynthetic was removed. Then 1 ml of extract was added, crushed with a cell
pigments, photosynthesis, and the activities of intracellular and crusher for 18 min, and centrifuged at 10000 rpm at 4 °C for 20 min.
extracellular enzymes of microalgae consortia were measured. The concentration of O2- was determined according to the assay kit
(Beijing Solabao Technology Co., Ltd., China).
In addition, the concentrations of SOD and MDA were
3.4 Measurement of cell growth determined using assay kits (Nanjing Jiancheng Bioengineering
Institute, China) following the manufacturer’s protocol. The
Algal cell growth in each treatment was measured according to extraction steps were as follows: 1) 10 mL samples were collected
Peng et al. (2014). The algal biomass of five microalgae was and centrifuged at 6500 rpm for 10 min at 4°C (Eppendorf,
calculated according to its linear relationship with an optical Centrifuge 5810R), 2) the biomass pellet was washed with 0.01 M
density of 680nm (OD680) (Supplementary Table 1). sodium phosphate buffer (pH 7.4), crushed and centrifuged, and 3)
the enzymes obtained from the disrupted microalgae were extracted
using 500mL of sodium phosphate buffer.
3.5 Measurement of photosynthetic
pigment and chlorophyll
fluorescence parameters 3.7 Determination of extracellular
enzyme content
The modified hot methanol method was used to determine the
cytochrome of microalgae (Xiong et al., 2016). Briefly, 5 mL of the The change characteristics of extracellular enzyme activity of
microalgae were determined by API ZYM (Merier, France), which
TABLE 1 Experimental setup of the combined experiment. is a laboratory kit for semi-quantitative analysis of the production of
hydrolytic enzymes. 3mL samples were collected and centrifuged at
Group name Microalgae species 7000 rpm for 15 min to obtain the supernatant containing
CK Abiotic control group extracellular enzymes. Supernatants (65 mL) were dispensed into
the 20 microcupules of API ZYM strips and incubated at 37°C for
C C. pyrenoidosa
4 h. After that, 30 mL of the commercial reagents ZYM A and ZYM
C+S C. pyrenoidosa + S. quadricauda B was added to all the microcupules to develop chromogenic
C+D C. pyrenoidosa + Dictyosphaerium sp. substrates. The color reactions were read after 5 min, and the
results were recorded at four grades from 0 to 4 according to the
C+S+D C. pyrenoidosa + S. quadricauda + Dictyosphaerium sp.
color rendering degree in the microcupules (Martinez et al., 2016).
3.8 Determination of antibiotic pyrenoidosa exposed to SMX was significantly higher than that of
concentration the control group (p<0.05) (Figure 2A), Fv/Fm of Dictyostelium sp.
and S. quadricauda were not higher (Figures 2B, C), for H. pluvialis
The samples were filtered through a 0.22 mm polytetrafluoroethylene and B. braunaii it was significantly lower (p<0.05) (Figures 2D, E).
(PTFE) filter, and the residual SMX concentration was determined by SMX did not exert obvious stress on C. pyrenoidosa, S. quadricauda,
high performance liquid chromatography (HPLC, Alliance 2695 system, and Dictyostelium sp. but promoted the synthesis of chlorophyll a
USA). A total of 2 mL of the samples was separated using an ACQUITY (Supplementary Figure 2). In addition, during the whole
UPLC BEH C18 chromatography Column [250 × 4.6 mm, 5 μm, waters, experimental period, the Yield of the microalgae treatment group
USA) at a column temperature of 30°C. The mobile phase consisted of a was different from that of the control group, especially that of the B.
mixture of solution A (acetonitrile and LCMS-grade water (30: 70 v/v)] braunaii treatment group showed a very significant decrease
and solution B (0.1% formic acid in LCMS-grade water) at a flow rate of (p<0.01) (Supplementary Figure 3).
1.0 mL/min.
The removal rate (%) of SMX was calculated using the following 4.1.2 SMX removal of five microalgae
formula: C. pyrenoidosa had the highest removal rate of SMX (49%)
during the 16 days exposure. The removal rates of S. quadricauda
Removal rate ( % ) = (1 − Ct =C0 ) 100 % and Dictyostelium sp. were 16% and 13%, respectively, and 10% and
where C0 (mg/L) is the initial concentration of antibiotics at 9% for H. pluvialis and B. braunii, respectively. (Figure 3).
time 0, and Ct (mg/L) is the concentration of antibiotics at time t. Compared with CK group, the removal rate of SMX by C.
pyrenoidosa was significantly increased (p<0.01), and higher than
that of the other four microalgae. Besides, the removal rate of SMX
3.9 Statistical analysis by S. quadricauda and Dictyostelium sp. was also significantly
higher than that of CK group (p<0.05).
To determine significant differences between the control and
treatment groups, a one-way analysis of variance (ANOVA) was
conducted using SPSS statistics software (version 22, Chicago, IL, 4.2 Microalgae consortia
USA). Values with p< 0.05 and p< 0.01 were considered significant
and very significant, respectively. To reveal the relationship between 4.2.1 Effects of SMX on pigment and
the variation of various indicators of microalgae and the photosynthetic parameters
experimental period, principal component analysis (PCA) was At different concentrations of SMX treatment, the content of
applied using R 4.1.1 statistical computation software. In addition, three photosynthetic pigments in each treatment group showed
hierarchy cluster analysis of the extracellular enzyme was upward trends (Figure 4). However, when the concentration of
performed using “pheatmap” package in R 4.1.1 statistical SMX increased to 50 mg/L, the total amount of photosynthetic
computation software. All figures were generated using Origin pigment synthesis decreased (Figure 4). When treated with 10 mg/L
2021 software (OriginLab, Northampton, MA, USA). SMX, there was no significant difference in chlorophyll a content
between the treatment and control groups at the end of the
experiment. However, when SMX reached 50 mg/L, the
4 Results chlorophyll a content of [C+S] group and [C+S+D] group was
significantly lower than in the control group (p<0.05) (Figures 4A,
4.1 Microalgal screening B). at the two SMX concentrations, the content of chlorophyll b
in [C+S] group fluctuated greatly, but the content of chlorophyll b
4.1.1 Effects of SMX on algal growth in [C+D] group and [C+S+D] group was significantly lower than in
and photosynthesis the control group (p< 0.05) (Figures 4C, D). In addition, the change
The growth ability of the five microalgae differed, and they also of carotenoids in the [C+S] group fluctuated more at the two SMX
had different responses to SMX (Figure 1). Compared with the concentrations. (Figures 4E, F). When exposed to SMX on day 2,
control groups, C. pyrenoidosa, S. quadricauda and Dictyostelium Fv/Fm of all groups was significantly inhibited (p<0.05) and then
sp. grew better in the medium containing SMX, especially the gradually recovered on day 3 (Supplementary Figures 4A, B). Yield
biomass of C. pyrenoidosa increased significantly compared with was significantly higher for all groups than in the control group
the control group (p<0.05) (Figures 1A, C, D). By contrast, the (p<0.05) (Supplementary Figures 4C, D).
growth of H. pluvialis and B. braunaii was significantly inhibited
(p<0.05) (Figures 1B, E). During the course of the experiment, the m 4.2.2 Effect of SMX on the intracellular
of five microalgae showed a decline to a varying degree, most antioxidant system of microalgae
pronounced for H. pluvialis. In the 10 mg/L SMX treatment, the concentrations of O2-, SOD,
After exposure to SMX, F v /F m of the microalgae was and MDA in the [C] group, [C+S] group, and [C+D] group showed
significantly reduced (p<0.05) except for S. quadricauda early in no significant change on day 3 and day 7 (Figures 5A, C, E). However,
the experiment (Figure 2). On day 14 and 16, Fv/Fm of C. the contents of O2-, SOD, and MDA in [C+S+D] group were
A B
C D
FIGURE 1
Differences in the algal biomass of five microalgae (A: C. pyrenoidosa, B: S. quadricauda, C: Dictyosphaerium sp., D: H pluvialis and E: B braunii)
during 16 days of cultivation. Error bars represent standard deviation (n=3). Columns with different letters indicate significant differences (p<0.05)
between the control and treatment.
significantly higher than in the other groups, and the content on day 3 treatment (p<0.05), and the content on day 3 was 1.7-fold higher than
was significantly higher than on day 7 (p<0.05). When SMX reached on day 7. The content of MDA in [C+S+D] group was not
50 mg/L, the contents of O2- and SOD in the [C+S+D] group significantly different from that in the 10 mg/L SMX treatment,
increased significantly compared with those in the 10 mg/L SMX and the MDA content exhibited no significant change on day 3 and
treatment, and the contents on day 3 were 1-fold and 1.6-fold higher day 7 as the [C] group and the [C+D] group (Figure 5F).
than on day 7, respectively. In addition, the contents of O2- and SOD
in [C] group and [C+S] group were significantly higher on day 3 than 4.2.3 Effects of SMX on extracellular enzyme
on day 7 (p<0.05), but there was no significant change in [C+D] The heatmap revealed changes in the activity of extracellular
group (Figures 5B, D). However, the content of MDA in [C+S] group enzymes throughout the experiment. Overall, 7 of the 20 enzymes
was significantly higher than in the other groups below 50 mg/L SMX could be detected on day 3 and day 7 of the experiment, respectively,
A B
C D
FIGURE 2
Effects of the SMX on the maximum photosynthetic efficiency (Fv/Fm) of five microalgae (A: C. pyrenoidosa, B: S. quadricauda, C: Dictyosphaerium
sp., D: H pluvialis and E: B braunii) during 16 days of cultivation. Error bars represent standard deviation (n=3). Asterisks indicate significant
differences between the control and treatment groups (p< 0.05*).
namely acid phosphatase, naphthol-AS-BI-phosphate, esterase (C4), phosphatase gradually decreased in all treatment groups but the
lipase esterase (C8), leucine aminopeptidase, valine aminopeptidase activity on day 7 of the experiment was higher than that on day 3.
and cystine aminopeptidase (Figure 6). At both SMX concentrations,
C4 and C8 did not appear on day 3 of the experiment but were 4.2.4 Removal of SMX by different
recorded on day 7, and the number and activity decreased with the microalgae consortia
increasing SMX concentration. The activities of the three The antibiotic content of treatment groups decreased
aminopeptidases were greater on day 3 than on day 7, being higher significantly during the first two days of the experiment (p<0.05).
in the [C+D] and [C+S+D] groups and lower in the [C+S] group. In With enhanced SMX concentration, the removal capacity of all
addition, with the increased SMX concentration, the activity of treatment groups decreased. [C] group and [C+S] group had a
FIGURE 3
The removal rate (%) of SMX for five microalgae during 16 days of exposure. CK was an abiotic control containing 5 mg/L SMX. Error bars represent
standard deviation (n=3). Asterisks indicate significant differences between the control and treatment groups (p< 0.05*; p< 0.01**).
stronger removal capacity than the other groups whose removal even decompose and absorb antibiotics as organic substances (Kiki
rates reached 11.1% and 10.9%, respectively, at 10 mg/L SMX. et al., 2020). Some studies have shown that Dictyosphaerium sp. can
However, although [C+D] and [C+S+D] had the capacity to remove combine with the photocatalyst bismuth vanadate (BiVO4) to
the different concentrations of SMX, SMX was significantly lower improve the removal of sulfamethazine (SM2), while H. pluvialis
than for the [C] and [C+S] groups (p<0.05) (Figure 7). can efficiently degrade nearly 10 kinds of antibiotics through
membrane photobioreactors (Chen et al., 2020c; Claude et al.,
4.2.5 PCA 2022). Due to the different modes of action and drug affinity of
On day 3 and 7 of the SMX treatment, PC1 constituted 50.4% of microalgae, these microalgae did not show high degradation rates
the total variance in all groups, and PC2 accounted for an additional for SMX in this study. Therefore, C. Pyrenoidosa is a keystone
30.3% (Figure 8). The distribution of each microalgae consortium in microalgae species for degrading SMX. However, when we
PCA was significantly different on day 3 and day 7 after exposure to combined C. pyrenoidosa with two other microalgae, the
SMX. On day 3, the contents of O2-, SOD, and MDA of the microalgae consortia showed a higher photosynthetic pigment
microalgae consortia had changed most significantly, while on content and stronger photosynthesis than in the C. pyrenoidosa
day 7 mainly photosynthetic pigment and chlorophyll monoculture after the SMX treatment, and the activity of O2- and
fluorescence parameters of microalgae consortia had changed. antioxidant enzymes in the antioxidant system of the microalgae
consortia was higher than for the individual microalgae. The PCA
results also revealed that the effects of SMX on antioxidant enzymes
5 Discussion and photosynthesis of microalgae consortia were stronger than that
of individual microalgae species C. pyrenoidosa. Chlorophyll is an
We found that the tolerance and removal efficiency of SMX by important component of light capture and energy transduction in
the five different microalgae varied. C. pyrenoidosa showed the photosynthesis and photosynthetic reaction, including light
highest removal capacity, followed by S. quadricauda and reaction, Calvin cycle, and starch synthesis. Some emerging
Dictyosphaerium sp. Previous studies have shown that Chlorella pollutants can affect the accumulation of chlorophyll content in
sp. has excellent dissipation to many antibiotics, so it is often the microalgae (Telfer, 2002; Del Campo et al., 2007). And the increase
preferred microalgae to remove antibiotics (Song et al., 2019; Chen in chlorophyll content in cells can serve as a protective mechanism
et al., 2020b; Xiong et al., 2020). Besides, C. pyrenoidosa had better to eliminate the accumulation of reactive oxygen species in
growth and higher photosynthesis in the culture system containing chloroplasts (Kasahara et al., 2002; Tsiaka et al., 2013; Xiong
SMX. Other studies have shown increases in the biomass of C. et al., 2016). In addition, the antioxidant properties of carotenoids
pyrenoidosa treated with both chloramphenicol and florfenicol (Lai can prevent lipid peroxidation by inhibiting the production of
et al., 2009). Stimulated growth might reflect that C. pyrenoidosa singlet oxygen and free radicals in cells, thereby promoting the
can resist the persecution of antibiotics at low concentrations and stability of photosynthesis and protecting cells (Jahns and
A B
C D
E F
FIGURE 4
Effects of two concentrations of SMX on the cytochromes of different microalgae consortia during 7 days. (A, C, E) 10 mg/L; (B, D, F) 50 mg/L. Error
bars represent standard deviation (n=3). Asterisks indicate significant differences between the control and treatment groups (p< 0.05*).
Holzwarth, 2012). When the antioxidant defense system is When C. pyrenoidosa was combined with other microalgae, the
stimulated under stress conditions, its antioxidant enzyme activity removal efficiency of SMX was never higher than that of the
will increase (Qian et al., 2008; Singh et al., 2018). Throughout the individual keystone microalgae C. pyrenoidosa, and it even declined
entire experimental cycle, the changes in photosynthetic pigments for some of the consortia: the removal efficiency of [C+S] to SMX was
and antioxidant enzyme activity of microalgae consortia were close to that of keystone microalgae C. pyrenoidosa, but the removal
significantly stronger than those of individual microalgae. Our capacity of [C+S+D] group was weaker than that of the keystone
results, therefore, indicate that the microalgae consortia were microalgae and the consortia containing two microalgae species ([C
more sensitive to SMX than the individual microalgae species, +S] and [C+D]). Thus, the increase in the richness of the different
which is consistent with the findings of Xiong et al. (2017b) in microalgae taxa did not improve the removal rate of SMX. Lima et al.
their study of the effects of enrofloxacin on microalgae consortia. (2003) found that the higher the proportion of the core microalgae
A B
C D
E F
FIGURE 5
Effects of different concentrations of SMX on the intracellular antioxidant system of microalgae consortia during 7 days of exposure. Superoxide
anion (O2-) content (A, B) superoxide dismutase (SOD) activity (C, D) malondialdehyde (MDA) content (E, F). Error bars represent standard deviation
(n=3). Columns with different letters indicate significant differences (p<0.05) between the control and treatment groups.
Coenochloris pyrenoidosa in the microalgae consortia, the higher the cultivation of multiple microalgae may lead to competition for
degradation rate of the microalgae consortia to p-nitrophenol. organic matter in the substrate, thereby affecting the metabolism of
Therefore, the lower removal rate of SMX in the [C+S+D] group SMX by the microalgae itself. In co-culture systems, there are
may likely be due to the decreased proportion of the keystone interactions between microalgae species, such as competition or
microalgae C. pyrenoidosa in the consortia because of resource symbiosis, so their capacity to degrade pollutants may also be
competition. Besides, In the study on the co- metabolism of SMX affected by such interactions (White et al., 2014; Goncalves et al.,
by C. pyrenoidosa, it was found that the addition of organic substrate 2017). The reported mechanisms of antibiotic removal by microalgae
or other nutrients can accelerate the potential of microbial include photolysis, biosorption, bioaccumulation, and intracellular
biodegradation of antibiotics, which is related to the characteristics and extracellular biodegradation (Liu et al., 2021; Yu et al., 2022). In
of additional substrates that can maintain biomass production and our study, photodegradation to remove SMX was negligible, and
act as an electron donor for co-metabolism of non-growth substrates other studies (Xiong et al., 2016; Xiong et al., 2017a) have revealed
(Xiong et al., 2020; Gao et al., 2022). Therefore, we speculate that co- that the biosorption and bioaccumulation of emerging pollutants
A B
FIGURE 6
Heat map analysis of the extracellular enzyme of different microalgae consortia on day 3 and day 7 of exposure to SMX at three concentrations: 10
mg/L (A) and 50 mg/L (B). The names of the different types of extracellular enzymes are distinguished by letters (C4, esterase; C8, lipase esterase; C,
Cystine aminopeptidase; A, Acid phosphatase; N, Naphthyl-AS-BI-phosphate; L, Leucine aminopeptidase; V, Valine aminopeptidase). The color of
each section is proportional to the significance of change in the biochemical parameters (red, relatively high; blue, relatively low).
contribute little to microalgae-mediated bioremediation. Therefore, microorganisms. Therefore, there may be a symbiotic relationship
the main mechanism of SMX removal by microalgae is the between C. pyrenoidosa and S. quadricauda (Johnson and Admassu,
biodegradation induced by the enzyme system of microalgae. The 2013; Koreiviene et al., 2014), which can maintain a stable growth
results of the API-ZYM reagent strip demonstrated that the state under SMX stress. The symbiotic relationship between them
extracellular enzyme activity of [C+S] group was always low even may have been destroyed when adding Dictyosphaerium sp., because
with the increase of SMX concentration, while the extracellular they need to mobilize more extracellular enzymes to decompose and
aminopeptidase and phosphatase of [C+D] group and [C+S+D] cycle nutrients. SMX can promote gene expression by affecting the
group were active. Nitrogen compounds are the products of the metabolic process of non-coding RNA, leading to the destruction of
nitrogen cycle, which begins with the hydrolysis of protein by the ultrastructure of microalgae cells, affecting the permeability of cell
aminopeptidase (Mayer, 1989), and phosphatase participates in the membranes and the activity of antioxidant enzymes (Xu et al., 2022).
phosphorus cycle (Martinez et al., 2016; Hu et al., 2022). The In this study, the intracellular antioxidant enzymes of the [C+S+D]
substances produced by these cycles can maintain the growth of group were always in a more active state, and the content of O2- was
A B
FIGURE 7
The removal rate (%) at different concentrations of SMX in three different microalgae consortia treatment groups (A) 10 mg/L; (B) 50 mg/L. Error bars
represent standard deviation (n=3). Asterisks indicate significant differences between the control and treatment groups (p< 0.05*; p< 0.01**).
FIGURE 8
Biplot from PCA integrating all the measured variables (Chl a, Chl b, carotenoids, Fv/Fm, Yield, MDA, O2-, SOD) for the two sampling times (days 3
and 7) and eight different treatments (shapes and colors interact in pairs: ◼-10 mg/L, ▲-50 mg/L; green-[C] group, red-[C+S] group, blue-[C+D]
group, orange-[C+S+D] group).
always at a higher level. This shows that under the stress of SMX, the review and editing. JS: Supervision, Writing-review and editing. SL:
balance between the generation and elimination of O2- was destroyed, Experiment, Analysis and Writing-review. LW: Supervision, Project
and the antioxidant enzyme could not effectively remove excessive administration. EJ: Writing-review and editing. WZ: Supervision,
O2-, so the physiological function and growth of the [C+S+D] group Project administration, Writing-review and editing. All authors
were under more severe stress. Additionally, when the SMX reached contributed to the article and approved the submitted version.
50 mg/L, the MDA content of [C+S] group significantly increased
and exceeded that of the other microalgae consortia. MDA is an
aldehyde product produced by lipid peroxidation caused by ROS,
which reflects the degree of lipid peroxidation of cell membranes and Funding
the tolerance of plants to stress (Celekli et al., 2013). MDA can also be
used as an indicator of NADPH-dependent oxidase and peroxidase, This work was supported by the Science and Technology
Commission of Shanghai Municipality (Grant No.19DZ1203404
which can effectively transform the accumulated pollutants in cells
through various enzyme systems (Apel and Hirt, 2013). Therefore, and No. 19DZ1204504). EJ was supported by the TÜBITAK
differently from the control group [C], the increase of MDA may be a program BIDEB2232 (project 118C250).
way to improve the antibiotic removal efficiency in the [C+S] group.
Acknowledgments
Data availability statement
We thank Anne Mette Poulsen for English edition and the
The original contributions presented in the study are included Shanghai Municipal Bureau of Ecology and Environment for its
in the article/Supplementary Material. Further inquiries can be support and assistance in this project.
directed to the corresponding author.
The remaining authors declare that the research was conducted reviewers. Any product that may be evaluated in this article, or
in the absence of any commercial or financial relationships that claim that may be made by its manufacturer, is not guaranteed or
could be construed as a potential conflict of interest. endorsed by the publisher.
References
Apel, K., and Hirt, H. (2013). Reactive oxygen species: metabolism, oxidative stress, Goncalves, A. L., Pires, J. C. M., and Simoes, M. (2017). A review on the use of
and signal transduction. Annu. Rev. Plant Biol. 55, 373–399. doi: 10.1146/ microalgal consortia for wastewater treatment. Algal. Res. 24, 403–415. doi: 10.1016/
annurev.arplant.55.031903.141701 j.algal.2016.11.008
Aydin, S., Unlu, I. D., Arabaci, D. N., and Duru, O. A. (2022). Evaluating the effect of Hu, X., Arif, M., Ding, D., Li, J., He, X., and Li, C. (2022). Invasive plants and species
microalga Haematococcus pluvialis bioaugmentation on aerobic membrane bioreactor richness impact litter decomposition in riparian zones. Front. Plant Sci. 13.
in terms of performance, membrane fouling and microbial community structure. Sci. doi: 10.3389/fpls.2022.955656
Total Environ. 807, 149908. doi: 10.1016/j.scitotenv.2021.149908 Huang, A., Yan, M., Lin, J., Xu, L., and Gong, H. (2021). A review of processes for
Bai, X., and Acharya, K. (2016). Removal of trimethoprim, sulfamethoxazole, and removing antibiotics from breeding wastewater. Int. J. Env. Res. Pub. He. 18, 4909.
triclosan by the green alga nannochloris sp. J. Hazard. Mater. 315, 70–75. doi: 10.1016/ doi: 10.3390/ijerph18094909
j.jhazmat.2016.04.067 Jahns, P., and Holzwarth, A. R. (2012). The role of the xanthophyll cycle and of
Brisson, J., Rodriguez, M., Martin, C. A., and Proulx, R. (2020). Plant diversity effect lutein in photoprotection of photosystem II. Biochimica. Et. Biophysica. Acta 1817,
on water quality in wetlands: a meta-analysis based on experimental systems. Ecol. 182–193. doi: 10.1016/j.bbabio.2011.04.012
Appl. 30, 904–915. doi: 10.1002/eap.2074 Johnson, K. R., and Admassu, W. (2013). Mixed algae cultures for low cost
Celekli, A., Kapi, M., and Bozkurt, H. (2013). Effect of cadmium on biomass, environmental compensation in cultures grown for lipid production and wastewater
pigmentation, malondialdehyde, and proline of Scenedesmus quadricauda var. remediation. J. Chem. Technol. Biotechnol. 88, 992–998. doi: 10.1002/jctb.3943
longispina. Bull. Environ. Contam. Toxicol. 91, 571–576. doi: 10.1007/s00128-013- Kasahara, M., Kagawa, T., Oikawa, K., Suetsugu, N., Miyao, M., and Wada, M.
1100-x (2002). Chloroplast avoidance movement reduces photodamage in plants. Nature 420,
Cheirsilp, B., Tippayut, J., Romprom, P., and Prasertsan, P. (2016). 829. doi: 10.1038/nature01213
Phytoremediation of secondary effluent from palm oil mill by using oleaginous Kiki, C., Rashid, A., Wang, Y., Li, Y., Zeng, Q., Yu, C. P., et al. (2020). Dissipation of
microalgae for integrated lipid production and pollutant removal. Waste Biomass antibiotics by microalgae: kinetics, identification of transformation products and
Valori. 8, 2889–2897. doi: 10.1007/s12649-016-9751-6 pathways. J. Hazard. Mater. 387, 121985. doi: 10.1016/j.jhazmat.2019.121985
Chen, S., Li, J., Feng, W., Yuan, M., Zhang, W., Xu, H., et al. (2020a). Biochemical Koreiviene, J., Valciukas, R., Karosiene, J., and Baltrenas, P. (2014). Testing of
responses of the freshwater microalga dictyosphaerium sp. upon exposure to three Chlorella/ scenedesmus microalgae consortia for remediation of wastewater, CO2
sulfonamides. J. Environ. Sci. 97, 141–148. doi: 10.1016/j.jes.2020.05.018 mitigation and algae biomass feasibility for lipid production. J. Environ. Eng. Landsc.
Chen, S., Yuan, M., Feng, W., Liu, W., Zhang, W., Xu, H., et al. (2020c). Catalytic 22, 105–114. doi: 10.3846/16486897.2013.911182
degradation mechanism of sulfamethazine via photosynergy of monoclinic BiVO4 and Kümmerer, K. (2009). Antibiotics in the aquatic environment – a review – part I.
microalgae under visible-light irradiation. Water Res. 185, 116220. doi: 10.1016/ Chemosphere 75, 417–434. doi: 10.1016/j.chemosphere.2008.11.086
j.watres.2020.116220
Lai, H., Hou, J., Su, C., and Chen, C. (2009). Effects of chloramphenicol, florfenicol,
Chen, S., Zhang, W., Li, J., Yuan, M., Zhang, J., Xu, F., et al. (2020b). Ecotoxicological and thiamphenicol on growth of algae Chlorella pyrenoidosa, Isochrysis galbana, and
effects of sulfonamides and fluoroquinolones and their removal by a green alga Tetraselmis chui. Ecotoxicol. Environ. Saf. 72, 329–334. doi: 10.1016/
(Chlorella vulgaris) and a cyanobacterium (Chrysosporum ovalisporum). Environ. j.ecoenv.2008.03.005
pollut. 263, 114554. doi: 10.1016/j.envpol.2020.114554
Li, X., Li, P., Lin, X., Zhang, C., Li, Q., and Gong, Z. (2008). Biodegradation of aged
Choudhury, M. I., McKie, B. G., Hallin, S., and Ecke, F. (2018). Mixtures of polycyclic aromatic hydrocarbons (PAHs) by microbial consortia in soil and slurry
macrophyte growth forms promote nitrogen cycling in wetlands. Sci. Total Environ. phases. J. Hazard. Mater. 150, 21–16. doi: 10.1016/j.jhazmat.2007.04.040
635, 1436–1443. doi: 10.1016/j.envpol.2020.114554
Li, Z., Li, M., Zhang, Z., Li, P., Zang, Y., and Liu, X. (2020). Antibiotics in aquatic
Claude, K., Ye, X., Li, X., Adyari, B., Hu, A., Qin, D., et al. (2022). Continuous environments of China: a review and meta-analysis. Ecotoxicol. Environ. Saf. 199,
antibiotic attenuation in algal membrane photobioreactor: performance and kinetics. J. 110668. doi: 10.1016/j.ecoenv.2020.110668
Haz. Mat. 434, 128910. doi: 10.1016/j.jhazmat.2022.128910
Del Campo, J. A., Mercedes, G. G., and Guerrero, M. G. (2007). Outdoor cultivation Lima, S. A. C., Castro, P. M. L., and Morais, R. M. (2003). Biodegradation of p-
of microalgae for carotenoid production: current state and perspectives. Appl. nitrophenol by microalgae. J. Appl. Phycol. 15, 137–142. doi: 10.1023/A:1023877420364
Microbiol. Biot. 74 (6), 1163–1174. doi: 10.1007/s00253-007-0844-9 Lima, S. A. C., Raposo, M. F. J., Castro, P. M. L., and Morais, R. M. (2004).
Diaz, S., and Cabido, M. (2001). Vive la difference: plant functional diversity matters to Biodegradation of p-chlorophenol by a microalgae consortium. Water Res. 38, 97–102.
ecosystem processes. Trends Ecol. Evol. 16, 646–655. doi: 10.1016/S0169-5347(01)02283-2 doi: 10.1016/j.watres.2003.09.005
Dutta, J., and Mala, A. A. (2020). Removal of antibiotic from the water environment Liu, R., Li, S., Tu, Y., and Hao, X. (2021). Capabilities and mechanisms of microalgae
by the adsorption technologies: a review. Water Sci. Technol. 82, 401–426. doi: 10.2166/ on removing micropollutants from wastewater: a review. J. Environ. Manage. 285,
wst.2020.335 112149. doi: 10.1016/j.jenvman.2021.112149
Fargione, J., Tilman, D., Dybzinski, R., Lambers, J., Clark, C., Harpole, W. S., et al. Liu, X., Lu, S. Y., Guo, W., Xi, B., and Wang, W. (2018). Antibiotics in the aquatic
(2007). From selection to complementarity: shifts in the causes of biodiversity- environments: a review of lakes, China. Sci. Total Environ. 627, 1195–1208.
productivity relationships in a long-term biodiversity experiment. Proc. R. Soc B. doi: 10.1016/j.scitotenv.2018.01.271
274, 871–876. doi: 10.1098/rspb.2006.0351
Loreau, M., and Hector, A. (2001). Partitioning selection and complementarity in
Gao, F., Yang, L., Chen, A., Zhou, W., Chen, D., and Chen, J. (2022). Promoting biodiversity experiments. Nature 413, 548–548. doi: 10.1038/35083573
effect of plant hormone gibberellin on co-metabolism of sulfamethoxazole by
microalgae Chlorella pyrenoidosa. Bioresour. Technol. 351, 126900. doi: 10.1016/ Martinez, D., Molina, M. J., Sanchez, J., Moscatelli, M. C., and Marinari, S. (2016).
j.biortech.2022.126900 API ZYM assay to evaluate enzyme fingerprinting and microbial functional diversity in
relation to soil processes. Biol. Fertil. Soils 52, 77–89. doi: 10.1007/s00374-015-1055-7
Ghorbannezhad, H., Moghimi, H., and Dastgheib, S. M. M. (2018). Evaluation of
heavy petroleum degradation using bacterial-fungal mixed cultures. Ecotoxicol. Mayer, L. M. (1989). Extracellular proteolytic enzyme activity in sediments of an
Environ. Saf. 164, 434–439. doi: 10.1016/j.ecoenv.2018.08.046 intertidal mudflat. Limnol. Oceanogr. 34, 973–981. doi: 10.4319/lo.1989.34.6.0973
Norberg, J. (2000). Resource-niche complementarity and autotrophic compensation Wang, Y., He, Y., Li, X., Nagarajan, D., and Chang, J. S. (2022). Enhanced
determines ecosystem-level responses to increased cladoceran species richness. biodegradation of chlortetracycline via a microalgae-bacteria consortium. Bioresour.
Oecologia 122, 264–272. doi: 10.1007/PL00008855 Technol. 343, 126149. doi: 10.2139/ssrn.3932832
Peng, F., Ying, G., Yang, B., Liu, S., Lai, H., Liu, Y., et al. (2014). Biotransformation of White, L. H., Martin, D. W., Witt, K. K., and Vogt, F. (2014). Impacts of nutrient
progesterone and norgestrel by two freshwater microalgae (Scenedesmus obliquus and competition on microalgae biomass production. J. Chemom. 28, 448–461. doi: 10.1016/
Chlorella pyrenoidosa): transformation kinetics and products identification. j.biortech.2021.126149
Chemosphere 95, 581–588. doi: 10.1016/j.chemosphere.2013.10.013 Xiao, G., Chen, J., Show, P., Yang, Q., Ke, J., Zhao, Q., et al. (2021). Evaluating the
Qian, H., Chen, W., Sheng, G. D., Xu, X., Liu, W., and Fu, Z. (2008). Effects of application of antibiotic treatment using algae-algae/activated sludge system.
glufosinate on antioxidant enzymes, subcellular structure, and gene expression in the Chemosphere 282, 130966. doi: 10.1016/j.chemosphere.2021.130966
unicellular green alga Chlorella vulgaris. Aquat. Toxicol. 88, 301–307. doi: 10.1016/ Xiong, Q., Hu, L. X., Liu, Y. S., Zhao, J., He, L., and Ying, G. (2021). Microalgae-
j.aquatox.2008.05.009 based technology for antibiotics removal: from mechanisms to application of
Rodrigues, D. A. D., da Cunha, C., Freitas, M. G., de Barros, A. L. C., Castro, P., innovational hybrid systems. Environ. Int. 155, 106594. doi: 10.1016/
Pereira, A. R., et al. (2020). Biodegradation of sulfamethoxazole by microalgae-bacteria j.envint.2021.106594
consortium in wastewater treatment plant effluents. Sci. Total Environ. 749, 141441. Xiong, J., Kurade, M. B., Abou-Shanab, R. A. I., Ji, M. K., Choi, J., Kim, J. O., et al.
doi: 10.1016/j.scitotenv.2020.141441 (2016). Biodegradation of carbamazepine using freshwater microalgae Chlamydomonas
Rodrigues, D. A. D., da Cunha, C., Santo, D. R. D., de Barros, A. L. C., Pereira, A. R., Silva, mexicana and Scenedesmus obliquus and the determination of its metabolic fate.
S. D., et al. (2021). Removal of cephalexin and erythromycin antibiotics, and their resistance Bioresour. Technol. 205, 183–190. doi: 10.1016/j.biortech.2016.01.038
genes, by microalgae-bacteria consortium from wastewater treatment plant secondary Xiong, J., Kurade, M. B., and Jeon, B. H. (2017a). Biodegradation of levofloxacin by
effluents. Environ. Sci. pollut. Res. Int. 28, 67822–67832. doi: 10.1007/s11356-021-15351-x an acclimated freshwater microalga, Chlorella vulgaris. Chem. Eng. J. 313, 1251–1257.
Singh, R., Upadhyay, A. K., Chandra, P., and Singh, D. P. (2018). Sodium chloride doi: 10.1016/j.cej.2016.11.017
incites reactive oxygen species in green algae Chlorococcum humicola and Chlorella Xiong, J., Kurade, M. B., and Jeon, B. H. (2017b). Ecotoxicological effects of
vulgaris: implication on lipid synthesis, mineral nutrients and antioxidant system. enrofloxacin and its removal by monoculture of microalgal species and their
Bioresour. Technol. 270, 489–497. doi: 10.1016/j.biortech.2018.09.065 consortium. Environ. pollut. 226, 486–493. doi: 10.1016/j.envpol.2017.04.044
Song, C., Wei, Y., Qiu, Y., Qi, Y., Li, Y., and Kitamura, Y. (2019). Biodegradability Xiong, J., Kurade, M. B., Kim, J. R., Roh, H.-S., and Jeon, B. H. (2017c). Ciprofloxacin
and mechanism of florfenicol via chlorella sp. UTEX1602 and L38: experimental study. toxicity and its co-metabolic removal by a freshwater microalga Chlamydomonas
Bioresour. Technol. 272, 529–534. doi: 10.1016/j.biortech.2018.10.080 mexicana. J. Hazard. Mater 323, 212–219. doi: 10.1016/j.jhazmat.2016.04.073
Telfer, A. (2002). What is b–carotene doing in the photosystem II reaction centre? Xiong, Q., Liu, Y. S., Hu, L., Shi, Z., Cai, W., He, L., et al. (2020). Co-Metabolism of
Phil. Trans. R. Soc Lond. B. 357, 1431–1440. doi: 10.1098/rstb.2002.1139 sulfamethoxazole by a freshwater microalga Chlorella pyrenoidosa. Water Res. 175,
Tsiaka, P., Tsarpali, V., Ntaikou, I., Kostopoulou, M. N., Lyberatos, G., and Dailianis, 115656. doi: 10.1016/j.watres.2020.115656
S. (2013). Carbamazepine-mediated pro-oxidant effects on the unicellular marine algal
Xu, D., Xie, Y., and Li, J. (2022). Toxic effects and molecular mechanisms of
species Dunaliella tertiolecta and the hemocytes of mussel Mytilus galloprovincialis.
sulfamethoxazole on Scenedesmus obliquus. Ecotox. Environ. Safe. 232, 113258.
Ecotoxicology 22 (8), 1208–1220. doi: 10.1007/s10646-013-1108-3
doi: 10.1016/j.ecoenv.2022.113258
Vacek, Z., Prokupkova, A., Vacek, S., Bulusek, D., Simunek, V., Hajek, V., et al.
(2021). Mixed vs. monospecific mountain forests in response to climate change: Yang, L., Ren, L., Tan, X., Chu, H., Chen, J., Zhang, Y., et al. (2020). Removal of
structural and growth perspectives of Norway spruce and European beech. For. Ecol. ofloxacin with biofuel production by oleaginous microalgae Scenedesmus obliquus.
Manage. 488, 119019. doi: 10.1016/j.foreco.2021.119019 Bioresour. Technol. 315, 123738. doi: 10.1016/j.biortech.2020.123738
Wan, J., Guo, P., Peng, X., and Wen, K. (2015). Effect of erythromycin exposure on Yu, C., Pang, H., Wang, J., Chi, Z., Zhang, Q., Kong, F., et al. (2022). Occurrence of
the growth, antioxidant system and photosynthesis of Microcystis flos-aquae. J. Hazard. antibiotics in waters, removal by microalgae-based systems, and their toxicological
Mater. 283, 778–786. doi: 10.1016/j.jhazmat.2014.10.026 effects: a review. Sci. Total Environ. 813, 151891. doi: 10.1016/j.scitotenv.2021.151891