1. Introduction
Production and research on carbon-based nanomaterials (CNMs) undergo exponential growth over last decade [
1,
2,
3]. In 2021 the global CNMs market was valued at
$2.9 billion and projected to reach
$31.6 billion by 2031 with a compound annual growth rate (CAGR) of 27.7% according to the report of Allied Market Research (
https://www.alliedmarketresearch.com/carbon-nano-materials-market (accessed on 05 May 2023)). The main representatives of the family of CNMs are fullerene (C60 or C70), graphene (Gr), graphene oxide (GrO), carbon nanotubes (CNTs), carbon quantum dots, and other derivatives [
4]. Freixa et al. (2018) stated in their analytical work that fullerenes are the most studied group of CNMs, followed by multiwalled carbon nanotubes (MWCNTs), single walled carbon nanotubes (SWCNTs), Gr, and black carbon [
5]. The interest in CNMs based on its unique mechanical, electrical, thermal, optical, and chemical properties [
6]. Consequently, this type of materials found application in various industrial processes and consumer products, such as drug carriers [
2,
7], bio-sensors [
8], compounds for bionanocomposites [
9,
10], energy conversion and storage devices [
11,
12], environmental purification [
13,
14], recovery of valuable compounds such as rare earth elements and other metals [
15,
16], etc.
The growing production and application of CNMs lead to rise in the risks of environmental contamination. CNMs could enter aquatic systems during life cycle of all the products containing CNMs including production, transportation, application, and disposal [
17]. In this point of view, risk assessment of CNMs became a crucial problem for maintaining safe environment and minimize human health risks.
Among the main obstacles coupled with risk assessment of nanomaterials (NMs) in aquatic environment should be highlighted (1) dependence of toxic behavior on particle physical and chemical properties, and variety of sizes, shapes, surface area, functionalization, etc. [
18,
19,
20], (2) fate and transformation of NMs in aquatic environment and in the contact with different organisms, which could alter initial toxicity of materials [
21,
22], and (3) species specific toxicity [
23,
24]. The following paragraphs briefly discuss these three points.
Particle size, functional groups, oxygen content, surface charges, hydrophobicity, and defect sites can be emphasized among the properties of CNMs defining its toxicity [
24]. The generally accepted fact that shorter CNTs have higher toxicity [
25], and SWCNTs are more toxic than MWCNTs [
26,
27]. Moreover, well dispersed CNMs are more toxic than their aggregates [
28]. The most likely explanation of this size-dependent toxicity is the larger specific surface area of smaller particles and higher interaction with membranes of organisms [
29]. The dose-response relationship in the toxicity assessment of CNMs is generally considered as not linear, where test organisms may adapt to low concentrations and high concentrations cause strong negative response [
30,
31]. The presence of functional groups is a controversial factor which reported as inhibitor or stimulator of the toxicity compared to pristine unmodified NMs [
32,
33,
34,
35]. Surface functionalization of CNMs often change their main toxic effect between mechanical damage and oxidative stress [
32,
36].
The environmental transformation of CNMs could significantly change their toxicity. Colloidal behavior and biodegradation of carbon-based and graphene family nanomaterials in aquatic environment were reviewed in our previous work [
21]. It was previously shown that the addition of natural organic matter (NOM) caused concentration-dependent disaggregation of fullerene C60 crystals and increased potential toxic effect [
37]. It also has been reported that the stability of GrO increases in the presence of humic acid and decreases in the presence of polysaccharides [
38]. In general, laboratory toxicity comparison of CNMs complicated with variations of the used exposure protocols, including type of NMs, exposure time, sample preparation methods, etc.
Based on the literature, the sensitivity to CNMs varied between aquatic species, where the most sensitive group of organisms were algae, followed by crustaceans, fish and bacteria [
5]. Considering microalgae as one of the most sensitive organisms to CNMs, it should be mentioned that the toxic effect of CNMs to algal cells might be associated not only with direct exposure, but also with shading effect (light absorption and photosynthesis prevention) and nutrient depletion [
39]. It also worth noting that the toxicity of CNMs varies between different microalgae species [
40,
41]. Several studies reported overproduction of extracellular polymeric substances in algae exposed to high concentrations of CNMs as a defensive mechanism [
42,
43].
This paper focused on the evaluation of toxic levels and toxic effects caused by the impact of four typical representatives of CNMs in marine microalgae. The chosen CNMs were multiwalled carbon nanotubes (CNTs), fullerene (C60), graphene powder (Gr), and graphene oxide (GrO). Microalgae was chosen as a sensitive and useful modes for NMs environmental risk assessment. In this work we used a red algae
Porphyridium purpureum which is commonly used in ecotoxicology [
44,
45,
46].
2. Materials and Methods
2.1. Nanoparticles
In this work we used four types of CNMs (
Table 1), namely multiwalled carbon nanotubes (CNTs), fullerene (C60), graphene powder (Gr), and graphene oxide (GrO). These types of NPs were chosen to compare the toxic effects and the impact between different CNMs in marine microalgae.
2.2. Microalgae Cultures and Exposure Protocol
The culture of a red alga
Porphyridium purpureum (Bory de Saint-Vincent) Drew et Ross, 1965 (Rhodophyta) originally isolated from the Peter the Great Bay (Sea of Japan, Far Eastern Russia) were provided by the Resource Collection
Marine Biobank of the National Scientific Center of Marine Biology, Far Eastern Branch of the Russian Academy of Sciences (NSCMB FEB RAS). The
P. purpureum species (
Figure S1) was chosen based on their abundance among microalgae in the Sea of Japan [
47] and their suitability as test-organisms in ecotoxicology [
48,
49]. The morphology and physiology of
P. purpureum was carefully described before [
50,
51].
Culturing of microalgae and toxicity test conditions were maintained following the guidance of OECD No.201 [
63] with minor modifications as stated below. Microalgae were cultured with Guillard’s f/2 medium [
52]. Filtered (pore diameter of the filter was 0.22 µm) and sterilized seawater with salinity 33 ± 1‰, pH 8.0 ± 0.2 was used for the experiments. The cultivation was carried out at a temperature of 20 ± 2 °C with an illumination intensity of 300 µmol photons/m
2s, and a light:dark cycle of 12:12 h.
Before the experiment microalgae cells were cultivated in 250 mL Erlenmeyer’s flasks. Algal culture in the exponential growth phase was taken for bioassays. For the experiment, microalgae cells were transferred to 12-well plates where each well contained 2 mL of microalgae aliquots and 2 mL of the tested sample to facilitate. The initial cell density in each well was 5–6×104 cells/mL. The wells of a control group had only microalgae aliquots with addition of 2 mL of f/2 medium. The other wells had different concentration of the prepared NP suspension.
The stock suspensions of the four used CNMs were prepared in filtered seawater to eliminate additional negative impact on microalgae associated with salinity reduction. The stock concentration for all the CNMs was 250 mg/mL. To prevent the agglomeration of NPs, the stock suspensions were sonicated with ultrasound homogenizer Bandelin Sonopuls GM 3100 (Bandelin Electronic GmbH & Co. KG, Germany) with the high-frequency power of 100W for 30 minutes. The sonication was performed on ice in 40 mL Sonopuls Rosette cell RZ-2 to prevent sample heating. The prepared stock suspensions were used to obtain the final exposure concentrations of 1, 10, 25, 50, 75, 100, 125 mg/mL. Each concentration and control group were carried out in triplicate. The duration of microalgae exposure to the NPs was 96 h.
2.3. Flow Cytometry Measurement
The method of flow cytometry was used to evaluate the growth rate inhibition, size and biochemical changes in microalgae cells after the exposure to NPs. All the measurements were carried out with CytoFLEX flow cytometer (Beckman Coulter, Indianapolis, IN, USA) with the software package CytExpert v.2.5. The staining by specific fluorescent dyes was used to evaluate biochemical changes and distinguish live and dead microalgae cells. The used endpoints, biomarkers, and parameters of their registration are represented in
Table 2. The excitation source and emission channels were selected according to the maximum emission of the used fluorescent dyes, provided by the manufacturer (Molecular Probes, Eugene, OR, USA). In all the cases, the excitation source was a blue laser (488 nm) of the CytoFLEX flow cytometer.
The determination and count of microalgae cells in the analyzed samples was carried out using the parameters of microalgae cell size, granularity, and fluorescence of chlorophyll a (emission channel 690 nm). The dead cells were excluded from the count by the presence of intensive fluorescence in 610 nm emission channel (PI positive cells).
For the measurement of esterase activity, membrane potential, and ROS generation, after 24 h of exposure, the sample from each well of 12-well plates was gently pipetted and 500 µL of liquid was transferred in 48-well plate and stained. The staining was made with PI and one of the other dyes. In general, all the measurements were performed three times separately, namely PI and FDA, PI and FDA, PI and H2DCFDA. PI was used to exclude dead cells, FDA, DiOC6, and H2DCFDA were used to assess esterase activity, membrane potential, and ROS generation changes, respectively.
The data of each well was collected at a flow rate of 100 μl/min until 2000 cells were registered. The mean fluorescence intensity (MFI) of the registered cells in the emission channel of 525 nm was used for comparison. For all the types of the used CNMs and for each couple of the dyes, negative and positive controls was measured directly before the measurement of the wells with exposed microalgae cells. Negative control group was prepared by 98 °C heat treatment of not exposed cells in dry block heater for 15 minutes. The wells with addition of only f/2 media were used as positive control. The obtained MFI data were normalized using positive control as 100% and negative control as 0%.
The growth rate inhibition was measured after 96 h of exposure. The sample from each well of 12-well plates was gently pipetted and 100 µL of liquid was transferred in 96-well plate and stained with only PI. The data of each well was collected at a flow rate of 100 μl/min for 30 sec. The obtained data was collected as number of cells per mL and then compared with the number of cells in control group. The changes in forward scatter intensity and the used size calibration kit allowed to compare the changes in cell size distribution after the exposure to CNMs.
Based on the results of growth rate inhibition, esterase activity, membrane potential, and ROS generation change assays, we calculated the effective concentrations of the CNMs, which caused 10% (EC10) and 50% (EC50) inhibition of listed endpoints. The calculation of EC10 and EC50 values were performed by nonlinear regression fit in GraphPad Prism 8.0.2 (GraphPad Software, San Diego, CA, USA).
2.4. Microscopy
Morphological changes of microalgae cells were observed and captured by optical microscope Axio Observer A1 (Carl Zeiss, Oberkochen, Germany).
2.5. Statistical Analysis
Statistical analyzes were performed by GraphPad Prism 8.0.2 (GraphPad Software, San Diego, CA, USA). The statistical significance was tested by one-way ANOVA. Normality residuals was checked by the Anderson-Darling test. A value of p ≤ 0.05 was considered statistically significant.
3. Results
The no observable effect level (NOEL) and calculated effective concentrations of the four used CNMs which caused 10% (EC10) and 50% (EC50) inhibition of microalgal growth rate and corresponding changes in esterase activity, membrane potential, and ROS generation in microalgae cells are given in
Table 3.
Based on the obtained results, the tested CNMs can be listed depending on the level of toxic exposure in red microalgae P. purpureum. Therefore, the growth rate and esterase activity of microalgae reduce in the following order: CNTs > GrO > Gr > C60. Sample CNTs had the highest adverse effect on growth rate and esterase activity of the microalgae cells. Moreover, only this sample caused membrane depolarization (DiOC6 fluorescence decrease) and significantly increased the level of ROS generation (increase in H2DCFDA fluorescence) in microalgae cells. Samples C60 and Gr caused high membrane hyperpolarization in microalgae cells (DiOC6 fluorescence increase). At the same time, GrO demonstrated no significant effect on cell membrane polarization even at the highest used concentration. It should be noted that sample C60, which had the lowest toxic impact in the microalgae, was the only used type of CNMs demonstrated no significant effect on ROS generation in the cells of P. purpureum.
The changes in esterase activity, membrane potential, and ROS generation depending on the concentration of the CNs are visualized in
Figure 1. All the data of statistical significance calculations related to these parameters are listed in
Table S1.
In addition to the data of
Table 3,
Figure 1a demonstrated that esterase activity of
P. purpureum exposed to C60 increased at lower concentrations (1-25 mg/L), maintained normal level at middle concentrations (50-75 mg/L), and rapidly decreased at the higher concentrations (100-125 mg/L).
The changes in the size of
P. purpureum cells are visualized in
Figure 2. The data of statistical significance calculations related to the changes in cell size are listed in
Table S2.
The highest toxic effect caused by the exposure of
P. purpureum to CNTs (
Table 3,
Figure 1) correlates with decrease in the size of microalgae cells at concentrations of 25 and 50 mg/L. The low concentrations of samples C60 (1-10 mg/L) and Gr (1 mg/L) also caused moderate decrease in cell size. The higher concentrations of samples Gr (50-125 mg/L) and GrO (25-125 mg/L) caused noticeable increase in cell size. The most pronounced increase in cell size was observed after the exposure to GrO, and the concentrations of 25-125 mg/L evoked enlargement of the cells to the size more than 10 µm.
The microscopic observation demonstrated that the cells of
P. purpureum can agglomerate with the big clusters of CNMs. This effect occurred in the cases of CNTs (
Figure 3a), Gr (
Figure 3c), and GrO (
Figure 3d). C60 did not forms clusters and demonstrated lower affinity to
P. purpureum cells, but small particles of C60 was absorbed to the surface of microalgae cells (
Figure 3b).
4. Discussion
Although CNMs are assumed as a substances with relatively low toxicity [
17,
53], they has a great variety of allotropic forms which could demonstrate different toxic properties in different species and conditions [
54,
55]. The present study was designed to determine the differences in the effect of multiwalled carbon nanotubes, fullerene, graphene powder, and graphene oxide in marine microalgae
P. purpureum.
In our previous work, a red algae
P. purpureum was more sensitive to the exposure of CNTs compare to thee the other marine microalgae species because of a highly hydrophobic surface of
P. purpureum cells covered with exopolysaccharide coverage [
40]. It is known that the surface properties of CNMs is one of the determinant factors of their toxicity [
56,
57]. The properties of graphene family nanomaterials and related biological interactions were carefully described in the work of Sanchez et al. 2011 [
58]. Fu and Zhang (2018) in their work explained the relationship between adhesion and hydrophobicity of NPs [
59]. Consequently, hydrophobic NPs would have higher affinity to hydrophobic regions of the cell membrane and result in higher potential of accumulation and penetration across the cells [
59,
60].
In this work, all the used CNMs did not have any functionalization or surface coatings and initially had hydrophobic properties. Among the tested samples, the highest hydrophobicity demonstrated sample CNTs, which rapidly agglomerated in seawater after the sonication at high concentrations, and probably had higher adhesion with
P. purpureum cells. This assumption correlate with the fact that CNTs were found to cause significantly higher toxic effect (
Table 3) in the microalgae compare to the other used nanomaterials. Moreover, fullerene C60 revealed the lowest toxicity toward microalgae cells (
Table 3) and was the only used CNM did not formed clusters in seawater even at the highest used concentration (125 mg/L) and demonstrated lower adhesion with the cells (
Figure 3b).
It is known that graphene family NMs could directly penetrate into the cell membrane of algae through cell pores [
61,
62,
63]. It was reported that GrO enters into the cells of
Chlorella vulgaris and damage organelles, enhanced the generation of ROS, and disrupted antioxidant enzymes [
64]. On the contrary, our study revealed decrease of ROS generation in
P. purpureum under the exposure to GrO (
Figure 1c).
As demonstrated in
Figure 3, the agglomerated flakes of CNTs, Gr, and GrO were sedimented and covered microalgae cells. In this case, it is important to notice the role of shading effect in the toxicity of CNMs. In photosynthetic microorganisms like microalgae or cyanobacteria CNMs prevent the illumination of the cells due to light absorption CNMs [
41,
65]. This effect could limit photosynthetic activity and growth rate of microalgae cells, and cause metabolic disruption [
41]. This effect might be the reason of observed cell size change in
P. purpureum exposed to high concentrations of CNTs, Gr, and GrO (
Figure 2).
It is known that unlike graphene family nanomaterials, CNTs are grown catalytically and often contain residual metal catalysts [
58]. The presence of trace metal residuals is another factor caused higher toxicity of CNTs in microalgae [
66]. The used sample of CNTs contain residuals of Al and Co (
Table 1). Although Co is one of the essential metals for cell function, it may become toxic at high concentrations [
67,
68]. It was reported that Al induced oxidative stress, ultrastructural changes, changes in lipid metabolism, degradation of cellular organelles, and suppression of antioxidant enzymatic activity in microalgae [
69,
70]. These facts explain ROS generation increase in
P. purpureum cells after the exposure to CNTs (
Figure 1c). Intensive ROS formation causes destruction in proteins, lipids, and carbohydrates, and lead to oxidative stress in microalgae [
68].
Previous work demonstrated that oxidized multiwalled CNTs and GrO caused esterase enzyme inhibition but no oxidative stress in cyanobacteria
Microcystis aeruginosa, and EC50 level 7.38 and 11.1 mg/L for CNTs and GrO, respectively [
71]. These results are in agreement with the results of current study (
Table 3) in fact that CNTs had higher cytotoxicity in microalgae than GrO. Several works shown that Gr was more toxic than CNTs in microalgae
Scenedesmus obliquus and
Chlorella pyrenoidosa [
39,
72], The other study with
Scenedesmus obliquus stated a higher toxicity of Gr compared to GrO [
61]. The results of these studies did not agree with our work as the sensitivity of
P. purpureum to the tested CNMs had the following order: CNTs > GrO > Gr > C60. Current observations reveal the importance of multispecies toxicity assessment and the role of surface properties both of microalgae cells and CNMs.
Considering the surface properties of CNMs in aquatic environment, it should be highlighted that NPs inevitably undergo not only physical transformation (agglomeration, sedimentation, etc.) but also obtain surface modification as a result of interaction with NOM, absorption of proteins and “biomolecular corona” formation [
21,
73]. This issue attracts attention of scientific community [
74,
75] and required further study in order to extend the understanding of fate and toxicity of CNMs. The “biomolecular corona” formation was not assessed in current work, nevertheless, this phenomenon should be taken into account in further study.