Next Article in Journal
Distribution Characteristics and Ecological Risk Assessment of Nitrogen, Phosphorus, and Some Heavy Metals in the Sediments of Yueliang Lake in Western Jilin Province, Northeast China
Previous Article in Journal
Monitoring and Analysis of Water Surface Slope of the Yarlung Zangbo River Based on Remote Sensing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Differential Response of Phaeodactylum tricornutum and Cylindrotheca fusiformis to High Concentrations of Cu2+ and Zn2+

1
State Key Laboratory of Marine Resource Utilization in the South China Sea, Hainan University, Haikou 570228, China
2
Laboratory of Development and Utilization of Marine Microbial Resource, Hainan University, Haikou 570228, China
3
Key Laboratory of Tropical Hydrobiology and Biotechnology of Hainan Province, Haikou 570228, China
4
College of Marine Sciences, Hainan University, Haikou 570228, China
5
School of Life Sciences, Hainan University, Haikou 570228, China
*
Author to whom correspondence should be addressed.
These authors contribute equally to this work.
Water 2022, 14(20), 3305; https://doi.org/10.3390/w14203305
Submission received: 18 August 2022 / Revised: 12 October 2022 / Accepted: 13 October 2022 / Published: 19 October 2022
(This article belongs to the Section Biodiversity and Functionality of Aquatic Ecosystems)

Abstract

:
Diatoms can be used as biosensors to assess aquatic environment quality, because they are widely distributed in almost all aquatic environments and show varied sensitivities toward heavy metal ions. The marine planktonic diatoms Phaeodactylum tricornutum (P. tricornutum) and Cylindrotheca fusiformis (C. fusiformis) are typical representatives of planktonic diatoms and benthic diatoms, respectively. C. fusiformis is very sensitive to changes in the concentration of heavy metal ions, and can be used as an indicator of the quality of the sedimental environment, while P. tricornutum can tolerate higher concentrations of heavy metal ions. To explore the potential difference in responses to heavy metal ions between planktonic and benthic diatoms, we compared the transcriptome of P. tricornutum and C. fusiformis under Cu2+ and Zn2+ treatment. The results indicated that P. tricornutum has several genes involved in ion transmembrane transport and ion homeostasis, which are significantly downregulated under Cu2+ and Zn2+ treatment. However, this enrichment of ion transmembrane transport- and ion homeostasis-related genes was not observed in C. fusiformis under Cu2+ and Zn2+ treatment. Additionally, genes related to heavy metal ion stress response such as peroxiredoxin, peroxidase, catalase, glutathione metabolism, phytochelatin, oxidative stress and disulfide reductase, were upregulated in P. tricornutum under Cu2+ and Zn2+ treatment, whereas most of them were downregulated in C. fusiformis under Cu2+ and Zn2+ treatment. This difference in gene expression may be responsible for the difference in sensitivity to heavy metals between P. tricornutum and C. fusiformis.

1. Introduction

The distribution and composition of biological communities are controlled or influenced by environmental variations such as disturbances, stressors, and biotic interactions and change in resources and hydraulic conditions [1]; therefore, such biological communities can be used as indicators of environmental conditions. Diatoms, for example, can be used as biosensors to assess aquatic environment quality, because diatoms are widely distributed in almost all aquatic environments [2], and different species of diatoms show varying sensitivities toward heavy metal ions [3]. Therefore, their species and distribution can be used as an indicator of the degree of heavy metal pollution in aquatic environments [4,5,6,7].
Since the relationship between diatoms and river pollution was revealed 70 years ago, the suitability of diatoms as bioassessment indicators for monitoring river quality has been demonstrated [8]. The sensitivity of diatoms to heavy metal ions is closely related to their response mechanisms [9]. Under high concentrations of heavy metal ions, diatoms tend to increase the synthesis of antioxidants or/and metal chelators, maintain ion balance through transporters, and increase extracellular carbohydrate production [9]. Moreover, it is reported that motile diatoms can tolerate higher concentrations of heavy metal ions than non-motile diatoms [2], indicating that there might be differences in response mechanisms between planktonic and benthic diatoms.
The marine planktonic diatom Phaeodactylum tricornutum (P. tricornutum) is rich in polyunsaturated fatty acids, lipids, and fucoxanthin [10]. Therefore, it can be used as a food for aquaculture animals and as raw materials for biodiesel and health products [11,12]. Additionally, due to its clear genomic background [13], universal molecular toolbox [14], and stable transgene expression system [15,16], P. tricornutum is also considered as a model single-cell organism for studying physiology, evolution, and biochemistry in microalgae. Cylindrotheca fusiformis (C. fusiformis) is a benthic diatom with a weakly silicificated cell wall, and is rich in nutrients which can induce the attachment and metamorphosis of benthic animal seedlings; thus, it can be used as open bait for sea cucumbers, abalones, sea urchins, and other marine treasure seedlings [17,18]. C. fusiformis grows rapidly under aerated conditions, and sinks to the bottom quickly after stopping aerating, making it very easy to be collected. In addition, the suitable temperature for most diatoms ranges from 10 to 25 °C, whereas the optimum temperature for C. fusiformis is approximately 30 °C. This can ensure the supply of seedling bait in the high-temperature season.
Therefore, P. tricornutum and C. fusiformis are typical representatives of planktonic and benthic diatoms, respectively. A comparative analysis of P. tricornutum and C. fusiformis will help to understand the different response mechanisms of planktonic and benthic diatoms. It is reported that C. fusiformis is very sensitive to changes in the concentration of heavy metals, and can be used as an indicator of the quality of the sedimental environment, while P. tricornutum can tolerate higher concentrations of heavy metal ions [4,6,19]. We propose that this may be related to their varying response mechanisms.
In this study, we aimed to explore the potential differential responses to heavy metal ions between planktonic and benthic diatoms. We compared the growth of P. tricornutum and C. fusiformis under different Cu2+ and Zn2+ concentrations, and transcriptome analyses were conducted. Moreover, we explored the mechanisms by which P. tricornutum responds to heavy metal ions, and why C. fusiformis is more sensitive to heavy metal ions.

2. Materials and Methods

2.1. Cell Culture and Treatments

P. tricornutum and C. fusiformis were obtained from the Microalgae Culture Center at the Ocean University of China. For P. tricornutum and C. fusiformis, algal cells were cultured using sterilized artificial seawater supplemented with f/2 nutrients at 20 °C and with four times of f/2 nutrients (2f) at 25 °C, respectively [20]. All cultures were grown under a 12:12 dark:light cycle under cool white fluorescent light (approximately 100 μmol m−2 s−1). Cell growth was detected by measuring the absorbance at 730 nm using a UV/visible spectrophotometer (UV-1800, Shimadzu, Tokyo, Japan).
For treatment with high concentrations of Cu2+ and Zn2+, P. tricornutum and C. fusiformis cells were treated with Cu2+ or Zn2+ at final concentrations of 30 μM and 60 μM. Control cells were cultured in f/2 (for P. tricornutum) or 2f medium (for C. fusiformis). Each treatment was performed in triplicate in 250 mL flasks. Cell growth was detected on days 0, 1, 3, 5 and 7.

2.2. Scanning Electron Microscope (SEM) and Energy Dispersive Spectroscopy (EDS) Analysis

For SEM-EDS analysis, P. tricornutum and C. fusiformis cells were treated with Cu2+ at a final concentration of 5 μM (PTCu and CFCu) and Zn2+ at a final concentration of 30 μM (PTZn and CFZn). Control cells were cultured in f/2 (for P. tricornutum, PTC) or 2f medium (for C. fusiformis, CFC). Each treatment was performed in triplicate in 2 L flasks. After 48 h, cell pellets were collected, further washed using distilled sea water, and centrifuged at 5000× g for 4 min. The pellets were fixed with 2.5% glutaraldehyde (4 °C) overnight and sequentially dehydrated for 15 min each in 30%, 50%, 70%, 80%, 90%, 100% and 100% EtOH, followed by CO2 critical point drying. Dried cells were placed on a conductive silicone rubber plate and treated with Gold sputtering, then viewed under the SEM (Hitachi’s TM4000 Plus, Hitachi Limited, Tokyo, Japan). EDS was performed with IXRF’s TM4-EDS.

2.3. Transcriptomic Analysis

For transcriptomic analysis, P. tricornutum and C. fusiformis cells were treated with Cu2+ at a final concentration of 5 μM (PTCu and CFCu) and Zn2+ at a final concentration of 30 μM (PTZn and CFZn). Control cells were cultured in f/2 (for P. tricornutum, PTC) or 2f medium (for C. fusiformis, CFC). Each treatment was performed in triplicate in 2 L flasks. After 48 h, cell pellets were collected, further washed using distilled sea water, and centrifuged at 5000× g for 4 min. The pellets were frozen in liquid nitrogen and stored at −80 °C.
Total RNA was extracted using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. High-quality total RNA (OD260/280 = 1.8–2.2, OD260/230 ≥ 2.0, RIN ≥ 6.5, 28S:18S ≥ 1.0, >2 μg) was used to construct cDNA libraries for high-throughput RNA sequencing. Overall, 1 μg of total RNA was used to construct an RNA-seq transcriptome library, using the TruSeqTM RNA sample preparation Kit from Illumina (Illumina, San Diego, CA, USA) as per the manufacturer’s instructions. Furthermore, cDNA libraries were selected for cDNA target fragments of 200–300 base pairs in 2% low-range ultra-agarose, and further amplified using Phusion DNA polymerase (New England Biolabs (Beijing), Beijing, China). The amplified cDNA libraries were loaded into a NovaSeq 6000 sequencing system Illumina (Illumina, San Diego, CA, USA).
To generate clean reads, raw sequence reads were trimmed using SeqPrep (https://github.com/jstjohn/SeqPrep accessed on 5 October 2016), and the quality of the raw reads was controlled using Sickle (https://github.com/najoshi/sickle accessed on 15 March 2015) with default parameters. The clean reads were annotated according to Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Groups of proteins (COG), the NCBI non-redundant protein sequences database (NR), Swiss-Prot, and Pfam databases. The mapped reads were further normalized using the reads per kb per million methods for the identification of differentially expressed genes (DEGs). Abundant genes were quantified using RSEM (http://deweylab.biostat.wisc.edu/rsem/ accessed on 14 February 2020) [21]. Differential gene expression was determined using the “edgeR” package in R (http://www.bioconductor.org/packages/2.12/bioc/html/edgeR.html accessed on 9 May 2018), based on the following threshold parameters: log2 fold-change > 2 and p-value < 0.05 [22]. Functional annotation and enrichment analyses were performed and classified using the GO and KEGG databases.

3. Results

3.1. Effects of Cu2+ and Zn2+ on Growth of P. tricornutum and C. fusiformis

To investigate the effect of Cu2+, Zn2+ on P. tricornutum and C. fusiformis, we compared the growth rates of P. tricornutum and C. fusiformis at different concentrations (0, 30 and 60 μM Cu2+ and 0, 30 and 60 μM Zn2+). The results showed that 60 μM Cu2+ significantly decreased the growth of P. tricornutum (p-value < 0.05), while the influence of 30 μM Cu2+ was not significant (p-value > 0.05), although there was also a tendency to decrease compared with that of the control group (Figure 1a). Both 30 μM and 60 μM Cu2+ significantly decreased the growth of C. fusiformis (Figure 1b). Neither 30 nor 60 μM Zn2+ significantly influenced the growth of P. tricornutum (Figure 1a), while 60 μM Zn2+ decreased the growth of C. fusiformis on day 5 (Figure 1b). As both 30 μM and 60 μM Cu2+ significantly decreased the growth of C. fusiformis, the concentration of Cu2+ for the transcriptomic analysis was set to a lower level (5 μM). As 30 μM Zn2+ did not significantly influence the growth of both P. tricornutum and C. fusiformis, the concentration of Zn2+ for the transcriptomic analysis was set to 30 μM.

3.2. Effects of Cu2+ and Zn2+ on Cell Morphology of P. tricornutum and C. fusiformis

To investigate the effect of Cu2+ and Zn2+ on cell morphology of P. tricornutum and C. fusiformis, we observed the cells with SEM. The results showed that both Cu2+ and Zn2+ did not significantly change the cell morphology of P. tricornutum (Figure 2a–c), while both Cu2+ and Zn2+ significantly changed the cell morphology of C. fusiformis (Figure 2d–f). This indicated that P. tricornutum was tolerant to Cu2+ and Zn2+, while C. fusiformis was more sensitive to Cu2+ and Zn2+.

3.3. Accumulation of Cu2+ and Zn2+ on Biosilica Shell of P. tricornutum and C. fusiformis

To investigate the accumulation of Cu2+ and Zn2+ on the biosilica shell of P. tricornutum and C. fusiformis, EDS was conducted to analyze the concentration of Si, Cu and Zn on the cell surface of P. tricornutum and C. fusiformis. The results showed that in the control group (containing 0.04 μM Cu2+ and 0.08 μM Zn2+ in the medium) of P. tricornutum, the content of Cu and Zn was 16.72% and 13.76% (Figure 3a), respectively. In the Cu2+ group (containing 5 μM Cu2+ and 0.08 μM Zn2+ in the medium) of P. tricornutum, the content of Cu and Zn was 18.70% and 13.75% (Figure 3b), respectively. In the Zn2+ group (containing 0.04 μM Cu2+ and 30 μM Zn2+ in the medium) of P. tricornutum (Figure 3c), the content of Cu and Zn was 6.53% and 17.76%, respectively. While in the control group (containing 0.16 μM Cu2+ and 0.32 μM Zn2+ in the medium) of C. fusiformis, the content of Cu and Zn was 0% (Figure 3d). In the Cu2+ group (containing 5 μM Cu2+ and 0.32 μM Zn2+ in the medium) of C. fusiformis, the content of Cu and Zn was 23.53% and 14.58% (Figure 3e), respectively. In the Zn2+ group (containing 0.16 μM Cu2+ and 30 μM Zn2+ in the medium) of C. fusiformis, the content of Cu and Zn was 9.92% and 20.81% (Figure 3f), respectively. These results indicated that both P. tricornutum and C. fusiformis accumulated Cu and Zn on the cell surface.

3.4. Effects of Cu2+ and Zn2+ on Gene Transcription in P. tricornutum

3.4.1. Annotation of P. tricornutum Transcriptome

To investigate the potential effect of Cu2+ and Zn2+ on gene transcription in P. tricornutum, we analyzed the transcriptome of P. tricornutum exposed to 5 μM Cu2+ (PTCu) and 30 μM Zn2+ (PTZn) for 48 h, with control (PTC) with no treatment of heavy metals. An average of 46,069,016 raw reads and 45,693,261 clean reads were generated from the total RNA extracted from P. tricornutum. A total of 98.39% of the clean read bases had a Q-value ≥ 20, and 94.99% of the clean read bases had a Q-value ≥ 30 (Table S1). De novo assembly generated 10,754 unigenes, including 10,167 known genes and 587 new genes. Figure 4 shows the length distribution of unigenes in P. tricornutum.
The acquired unigenes were annotated according to the GO, KEGG, COG, NR, Swiss-Prot, and Pfam databases. Of all the assembled unigenes, 82.28%, 46.02%, 72.75%, 99.19%, 56.1%, and 74.18% were annotated by GO, KEGG, COG, NR, Swiss-Prot, and Pfam, respectively (Figure 5, Table S2).

3.4.2. Identification and Functional Enrichment Analysis of Different Express Genes (DEGs) in P. tricornutum upon Cu2+ Treatment

Transcriptome analysis of DEGs in P. tricornutum exposed to 5 μM Cu2+ was performed, using high-throughput RNA sequencing. A total of 2006 genes, including 1119 up- and 887 downregulated genes were detected to be significantly regulated (p < 0.05) under Cu2+ treatment, with a 2-fold change in abundance considered as the criterion of biologically significant difference (Table S3). DEGs were classified into three main functional categories of GO terms: molecular function (MF), biological process (BP), and cellular component (CC; Figure 6). The GO enrichment analysis for upregulated genes is shown in Figure 6a, in which only 20 annotation categories with the most significantly enriched DEPs are shown. For BP, DEGs were assigned to 13 subcategories involved in photosynthesis, carbon metabolism, and energy metabolism, with the three most abundant clusters being ’protein-chromophore linkage’, ‘photosynthesis, light harvesting in photosystem I’, and ‘photosynthesis, light harvesting’. For CC, DEGs were classified into five subcategories involved in photosynthesis, ’thylakoid membrane’, ‘chloroplast thylakoid membrane’, ‘plastid thylakoid membrane’, ‘photosynthetic membrane’, and ‘light-harvesting complex’. In the MF category, DEGs were divided into the two subcategories ‘chlorophyll-binding’ and ‘tetrapyrrole binding’. The GO enrichment analysis for downregulated genes is shown in Figure 6b, in which only 20 annotation categories with the most significantly enriched DEPs are shown. For BP, DEGs were assigned to 11 subcategories involved in metal ion homeostasis, cation homeostasis, and ion transport. For CC, DEGs were classified into four subcategories involved in the integral component of (plasma) membrane and intrinsic component of (plasma) membrane. In the MF category, the DEGs were divided into four subcategories involved in (inorganic) cation and inorganic molecular entity transmembrane transporter activity.
Overall, 19 DEGs involved in heavy metal ion stress response are listed in Table 1. These genes were mainly related to antioxidants such as peroxiredoxin, peroxidase, catalase, glutathione metabolism, phytochelatin, oxidative stress, and disulfide reductase. Most (14 out of 19) of these genes were upregulated, indicating their important roles in response to the high concentration of heavy metal ions.

3.4.3. Effects of Zn2+ on Gene Transcription in P. tricornutum

Transcriptome analysis of differential gene expression in P. tricornutum exposed to 30 μM Zn2+ was performed, using high-throughput RNA sequencing. A total of 4043 genes, including 2184 up- and 1859 downregulated genes were detected to be significantly regulated (p < 0.05) under Zn2+ treatment (Table S4). The GO enrichment analysis for DEGs in P. tricornutum under Zn2+ treatment is shown in Figure 7, in which only 20 annotation categories with the most significantly enriched DEGs are shown. The GO enrichment for DEGs in P. tricornutum under Zn2+ treatment was similar to that under Cu2+ treatment, in which the upregulated genes were mainly involved in photosynthesis (Figure 7a), whereas the downregulated genes were mainly involved in ion homeostasis, cation homeostasis, and ion transport (Figure 7b).
Overall, 24 DEGs involved in heavy metal ion stress response are listed in Table 2. These genes were mainly related to antioxidants such as peroxidase, catalase, peroxiredoxin, glutathione metabolism, phytochelatin biosynthetic process, oxidative stress, mutase superoxide dismutase, and disulfide reductase. Most (17 out of 24) of these genes were upregulated, indicating their important roles in response to the high concentration of heavy metal ions.

3.5. Effects of Cu2+ and Zn2+ on Gene Transcription in C. fusiformis

3.5.1. Annotation of C. fusiformis Transcriptome

To investigate the potential effect of Cu2+ and Zn2+ treatment on transcription in C. fusiformis, we analyzed the transcriptome of C. fusiformis exposed to 5 μM Cu2+ (CFCu) and 30 μM Zn2+ (CFZn) for 48 h, with no addition of heavy metal ions as the control (CFC). An average of 43,832,802 raw reads and 43,323,647 clean reads were generated from total RNA extracted from C. fusiformis. A total of 98.25% of the clean read bases had a Q-value ≥ 20, and 94.64% of the clean read bases had a Q-value ≥ 30 (Table S5). De novo assembly generated 26,146 unigenes. Figure 8 shows the length distribution of unigenes.
The acquired unigenes were annotated according to GO, KEGG, COG, NR, Swiss-Prot, and Pfam databases. Of all the assembled unigenes, 36.72%, 35.9%, 56.07%, 38.88%, 39.87%, and 55.78% were annotated by GO, KEGG, COG, NR, Swiss-Prot, and Pfam, respectively (Figure 9, Table S6).

3.5.2. Effects of Cu2+ on Gene Transcription in C. fusiformis

Transcriptome analysis of differential gene expression in C. fusiformis exposed to 5 μM Cu2+ was performed, using high-throughput RNA sequencing. A total of 1133 genes, including 315 up- and 818 downregulated genes were detected to be significantly regulated (p < 0.05) under Cu2+ treatment (Table S7). The GO enrichment analysis for upregulated genes is shown in Figure 10a, in which only 20 annotation categories with the most significantly enriched DEPs are shown. For BP, DEGs were assigned to 17 subcategories involved in signal transduction, nucleotide biosynthetic, organophosphate biosynthetic, etc. For CC, DEGs were classified into 1 subcategory, the plasma membrane. In the MF category, the unigenes were divided into 2 subcategories, 3′,5′-cyclic-nucleotide phosphodiesterase activity, and cyclic-nucleotide phosphodiesterase activity. The GO enrichment analysis for downregulated genes is shown in Figure 10b. For BP, no DEGs were enriched. For CC, DEGs were classified into 2 subcategories as an intrinsic component of the membrane and integral component of the membrane. In the MF category, the unigenes were divided into 3 subcategories, phospholipid transporter, glutamyl-tRNA reductase, and lipase activities.
Overall, 8 DEGs involved in antioxidants are listed in Table 3, including 1 peroxiredoxin, 1 glutathione synthetase, 1 glutathione S-transferase, 1 glutathione peroxidase, 1 hydroxyacylglutathione hydrolase, 1 deaminated glutathione amidase, and 1 peroxinectin. In total, 7 out of 8 genes were downregulated, indicating considerable differences between C. fusiformis and P. tricornutum in response to the high concentration of heavy metal ions.

3.5.3. Effects of Zn2+ on gene transcription in C. fusiformis

Transcriptome analysis of differential gene expression in C. fusiformis exposed to 30 μM Zn2+ was performed using high-throughput RNA sequencing. A total of 1900 genes, including 854 up- and 1046 downregulated genes were detected to be significantly regulated (p < 0.05) under Zn2+ treatment (Table S8). The GO enrichment analysis for upregulated genes is shown in Figure 11a. For BP, DEGs were enriched in 1 subcategory of cellular modified amino acid metabolic process. For CC, DEGs were classified into 3 subcategories, 3-oxoacyl-[acyl-carrier-protein] synthase activity, arginase activity, and cullin family protein binding. The GO enrichment analysis for downregulated genes is shown in Figure 11b, in which only 20 annotation categories with the most significantly enriched DEPs are shown. For BP, DEGs were assigned to 8 subcategories involved in the regulation of biological quality, homeostasis, posttranslational modification (amino acid modification), organelle assembly, response to topologically incorrect protein, etc. For CC, DEGs were classified into 3 subcategories as an intrinsic component of membrane, an integral component of membrane, and endoplasmic reticulum lumen. In the MF category, the unigenes were divided into 9 subcategories involved in catalytic activity, ATPase activity, tubulin (cytoskeletal protein, calcium ion, microtubule) binding, primary active transmembrane transporter activity, and protein kinase activity.
Overall, 16 DEGs involved in antioxidants are listed in Table 4, including 1 thioredoxin-like protein, 1 glutathione synthetase, 4 glutathione S-transferase, peroxiredoxin, 1 glutathionyl-hydroquinone reductase, 1 glutathione peroxidase, 2 phytochelatin biosynthesis-related genes, 1 light-harvesting complex stress-related protein, 1 thyroid peroxidase, 1 oxidative stress-related Abc1-like protein, 1 catalase-peroxidase, 1 methionine sulfoxide reductase, and 1 peroxinectin. Half of them were downregulated, and the rest were upregulated, which is different from the result under Cu2+ treatment.

4. Discussion

Cu2+ and Zn2+ are crucial micronutrients for diatoms. When Cu2+ and Zn2+ are present in an adequate amount, diatoms exhibit a stronger fitness and grow faster. Cu2+ and Zn2+ are components of many enzymes in algae cells. Cu2+ is involved in the electron transport of photosynthesis by serving as a ligand of cytochrome oxidase and plastocyanin [23,24]. In addition, Cu2+ is a component of Cu-tyroninase and multicopper oxidase, which is involved in Fe-deficiency response [25,26]. Zn2+ is an important component for carbonic anhydrases, which are involved in CO2 fixation, and zinc finger transcription factors, which are involved in gene transcription [27,28]. In addition, both Cu2+ and Zn2+ are important for Cu/Zn-SOD (superoxide dismutase) which is involved in anti-oxidation [29]. However, excess Cu2+ or Zn2+ will interfere with cellular physiology and biological processes, resulting in decreased cell growth and even death. Different types of cells have different types and amounts of enzymes, thus their demands for Cu2+ and Zn2+ are various. Meanwhile, as the shielding and permeation properties of cell membranes for heavy metal ions are different in various species, their tolerances to heavy metal ions are also various. In this study, the growth of both P. tricornutum and C. fusiformis was inhibited at 60 μM Cu2+, while 30 μM Cu2+ decreased the growth of C. fusiformis, yet did not have significant effect on the growth of P. tricornutum (Figure 1). Neither 30 nor 60 μM Zn2+ significantly influenced the growth of P. tricornutum (Figure 1a), while 60 μM Zn2+ decreased the growth of C. fusiformis on day five (Figure 1b). This indicated that P. tricornutum and C. fusiformis show different sensitivities to Cu2+ and Zn2+.
To explore the mechanism underlying the difference in susceptibility to heavy metals between P. tricornutum and C. fusiformis, transcriptomic analysis was conducted. Ion transport is reported to be a response mechanism to the high concentration of heavy metal ions. In this study, it has been shown that under high concentrations of both Cu2+ and Zn2+, most DEGs involved in photosynthesis were upregulated, indicating the effect of both Cu2+ and Zn2+ on photosynthesis in P. tricornutum. Meanwhile, most genes downregulated in P. tricornutum under Cu2+ treatment were involved in metal ion homeostasis and transmembrane ion transport. This indicated that ion homeostasis and transmembrane transport might be the main mechanisms for P. tricornutum to respond to high Cu2+ concentrations. Moreover, this enrichment of downregulated genes in metal ion transport was observed in P. tricornutum under Zn2+ treatment. However, the enrichment of downregulated genes in metal ion homeostasis-related genes did not occur under Zn2+ treatment, indicating a different response mechanism for Zn2+ to that for Cu2+.
Besides genes related to metal ion homeostasis and transmembrane ion transport, some other genes were previously reported to be involved in heavy metal stress response, including genes related to catalase, antioxidation, ascorbate metabolism, glutathione metabolism, phytochelatin, and oxidative stress [9]. These genes are listed in Table 1, Table 2, Table 3 and Table 4. Most of these genes were upregulated in P. tricornutum under both Cu2+ and Zn2+ treatments; however, only a few were upregulated in C. fusiformis, indicating that the response of C. fusiformis to heavy metal ion stress is different from that of P. tricornutum. Moreover, the enrichment of DEGs in ion homeostasis and transmembrane transport-related genes was not observed in C. fusiformis either. It is reported that C. fusiformis is sensitive to heavy metal ions [4,6], whereas P. tricornutum is more tolerant to Cu2+ stress [19]. We suspect that difference in gene expression might be one of the mechanism’s responses to the difference in susceptibility to heavy metals between P. tricornutum and C. fusiformis.
In addition, since the metal toxicity for cells is more related to intracellular metal bioaccumulation than to the metal concentration in water, and the fact that both P. tricornutum and C. fusiformis are widely considered biofilm-producing mixed diatoms [30,31,32,33], the role of metal ion management of the biofilm should be considered when exploring the mechanism underlying the difference in susceptibility to heavy metals between P. tricornutum and C. fusiformis. Using SEM and EDS analysis, we found that both P. tricornutum and C. fusiformis accumulated Cu2+ and Zn2+ onto the biosilica shell. In future work it would be informative to determine the intracellular concentrations of Cu2+ and Zn2+.

5. Conclusions

Transcriptome analysis of P. tricornutum and C. fusiformis under Cu2+ and Zn2+ treatments indicated that genes involved in metal ion homeostasis and transmembrane ion transport, and those related to catalase, antioxidation, ascorbate metabolism, glutathione metabolism, phytochelatin, and oxidative stress, might play important roles in the response of diatoms to heavy metal stress.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w14203305/s1, Table S1: Clean data statistics for Phaeodactylum tricornutum (P. tricornutum) transcriptome; Table S2: Annotation statistics for P. tricornutum transcriptome; Table S3: Annotations for differentially expressed genes (DEGs) in P. tricornutum under Cu2+ treatment; Table S4: Annotations for DEGs in P. tricornutum under Zn2+ treatment; Table S5: Clean data statistics for Cylindrotheca fusiformis (C. fusiformis) transcriptome; Table S6: Annotation statistics for C. fusiformis transcriptome; Table S7: Annotations for DEGs in C. fusiformis under Cu2+ treatment; Table S8: Annotations for DEGs in C. fusiformis under Zn2+ treatment.

Author Contributions

Conceptualization, A.H.; methodology, A.H.; software, Y.W.; validation, Y.W., J.D. and S.G.; formal analysis, Y.W.; investigation, Y.W.; resources, J.D.; data curation, Y.W.; writing—original draft preparation, A.H. and Y.W.; writing—review and editing, A.H. and Z.X.; visualization, J.D.; supervision, A.H.; project administration, A.H.; funding acquisition, Z.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (41876158), the Natural Science Foundation of Hainan Province (420QN219), and the Scientific Research Foundation of Hainan University (KYQD(ZR)20060).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Materials; further inquiries can be directed to the corresponding authors.

Acknowledgments

We would like to express our thanks to Xinchun Zhang (Chinese Academy of Tropical Agricultural Sciences Environment and Plant Protection Institute) for technical service of SEM and EDS.

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Larned, S.T. A prospectus for periphyton: Recent and future ecological research. J. N. Am. Benthol. Soc. 2010, 29, 182–206. [Google Scholar] [CrossRef] [Green Version]
  2. Butcher, R.W. Studies in the Ecology of Rivers: VII. The Algae of Organically Enriched Waters. J. Ecol. 1947, 35, 186–191. [Google Scholar] [CrossRef]
  3. Tudesque, L.; Grenouillet, G.; Gevrey, M.; Khazraie, K.; Brosse, S. Influence of small-scale gold mining on French Guiana streams: Are diatom assemblages valid disturbance sensors? Ecol. Indic. 2012, 14, 100–106. [Google Scholar] [CrossRef]
  4. Sbihi, K.; Cherifi, O.; Bertrand, M. Toxicity and biosorption of chromium from aqueous solutions by the diatom Planothidium lanceolatum (Brébisson) Lange-Bertalot. Am. J. Sci. 2012, 3, 27–38. [Google Scholar] [CrossRef]
  5. De Stefano, L.; Rotiroti, L.; De Stefano, M.; Lamberti, A.; Lettieri, S.; Setaro, A.; Maddalena, P. Marine diatoms as optical biosensors. Biosens. Bioelectron. 2009, 24, 1580–1584. [Google Scholar] [CrossRef]
  6. Marie, M.; Kirsten, H.; Pamela, Q.; Negri, A.P. Additive toxicity of herbicide mixtures and comparative sensitivity of tropical benthic microalgae. Mar. Pollut. Bull. 2010, 60, 1978–1987. [Google Scholar]
  7. Satoh, A.; Vudikaria, L.Q.; Kurano, N.; Miyachi, S. Evaluation of the sensitivity of marine microalgal strains to the heavy metals, Cu, As, Sb, Pb and Cd. Environ. Int. 2005, 31, 713–722. [Google Scholar] [CrossRef]
  8. Rimet, F. Recent views on river pollution and diatoms. Hydrobiologia 2012, 683, 1–24. [Google Scholar] [CrossRef]
  9. Masmoudi, S.; Nguyen-Deroche, N.; Caruso, A.; Ayadi, H.; Morant-Manceau, A.; Tremblin, G.; Bertrand, M.; Schoefs, B. Cadmium, copper, sodium and zinc effects on diatoms: From heaven to hell—A review. Cryptogam. Algol. 2013, 34, 185–225. [Google Scholar] [CrossRef]
  10. Owens, T.G.; Wold, E.R. Light-Harvesting Function in the Diatom Phaeodactylum tricornutum: I. Isolation and Characterization of Pigment-Protein Complexes. Plant Physiol. 1986, 80, 732. [Google Scholar] [CrossRef] [Green Version]
  11. Patil, V.; Reitan, K.I.; Knutsen, G.; Mortensen, L.M.; Källqvist, T.; Olsen, E.; Vogt, G.; Gislerød, H.R. Microalgae as source of polyunsaturated fatty acids for aquaculture. Curr. Top. Plant Biol. 2005, 6, 57–65. [Google Scholar]
  12. Alipanah, L.; Rohloff, J.; Winge, P.; Bones, A.M.; Brembu, T. Whole-cell response to nitrogen deprivation in the diatom Phaeodactylum tricornutum. J. Exp. Bot. 2015, 66, 6281–6296. [Google Scholar] [CrossRef] [PubMed]
  13. Bowler, C.; Allen, A.E.; Badger, J.H.; Grimwood, J.; Jabbari, K.; Kuo, A.; Maheswari, U.; Martens, C.; Maumus, F.; Otillar, R.P. The Phaeodactylum genome reveals the evolutionary history of diatom genomes. Nature 2008, 456, 239–244. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Siaut, M.; Heijde, M.; Mangogna, M.; Montsant, A.; Coesel, S.; Allen, A.; Manfredonia, A.; Falciatore, A.; Bowler, C. Molecular toolbox for studying diatom biology in Phaeodactylum tricornutum. Gene 2007, 406, 23–35. [Google Scholar] [CrossRef]
  15. De Risco, V.; Raniello, R.; Maumus, F.; Rogato, A.; Bowler, C.; Falciatore, A. Gene silencing in the marine diatom Phaeodactylum tricornutum. Nucleic Acids Res. 2009, 37, 96. [Google Scholar] [CrossRef] [Green Version]
  16. Stukenberg, D.; Zauner, S.; Dell’Aquila, G.; Maier, U.G. Optimizing CRISPR/Cas9 for the Diatom Phaeodactylum tricornutum. Front. Plant Sci. 2018, 9, 740. [Google Scholar] [CrossRef] [Green Version]
  17. Kawamura, T.; Roberts, R.D.; Nicholson, C.M. Factors affecting the food value of diatom strains for post-larval abalone Haliotis iris. Aquaculture 1998, 160, 81–88. [Google Scholar] [CrossRef]
  18. Gallardo, W.G.; Buen, S.M.A. Evaluation of mucus, Navicula, and mixed diatoms as larval settlement inducers for the tropical abalone Haliotis asinina. Aquaculture 2003, 221, 357–364. [Google Scholar] [CrossRef]
  19. Cid, A.; Torres, E.; Herrero, C.; Abalde, J.E. Disorders provoked by copper in the marine diatom Phaeodactylum tricornutum in short-time exposure assays. Cah. Biol. Mar. 1997, 38, 201–206. [Google Scholar]
  20. Guillard, R. Culture of Marine Invertebrate Animals, 1st ed.; Springer: New York, NY, USA, 1975; pp. 29–60. [Google Scholar]
  21. Li, B.; Dewey, C.N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011, 12, 323. [Google Scholar] [CrossRef] [Green Version]
  22. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Hänsch, R.; Mende, R.R. Physiological functions of mineral micronutrients (Cu, Zn, Mn, Fe, Ni, Mo, B, Cl). Curr. Opin. Plant Biol. 2009, 12, 259–266. [Google Scholar] [CrossRef] [PubMed]
  24. Peers, G.; Price, N.M. Copper-containing plastocyanin used for electron transport by an oceanic diatom. Nature 2006, 441, 341–344. [Google Scholar] [CrossRef] [PubMed]
  25. Peers, G.; Quesnel, S.A.; Price, N.M. Copper requirements for iron acquisition and growth of coastal and oceanic diatoms. Limnol. Oceanogr. 2005, 50, 1149–1158. [Google Scholar] [CrossRef]
  26. Maldonado, M.T.; Allen, A.E.; Chong, J.S.; Lin, K.; Leus, D.; Karpenko, N.; Harris, S.L. Copper-dependent iron transport in coastal and oceanic diatoms. Limnol. Oceanogr. 2006, 51, 1729–1743. [Google Scholar] [CrossRef] [Green Version]
  27. Cox, E.H.; McLendon, G.L.; Morel, F.M.M.; Lane, T.W.; Prince, R.C.; Pickering, I.J.; George, G.N. The Active Site Structure of Thalassiosira weissflogii Carbonic Anhydrase 1. Biochemistry 2000, 39, 12128–12130. [Google Scholar] [CrossRef]
  28. Rayko, E.; Maumus, F.; Maheswari, U.; Jabbari, K.; Bowler, C. Transcription factor families inferred from genome sequences of photosynthetic stramenopiles. New Phytol. 2010, 188, 52–66. [Google Scholar] [CrossRef]
  29. Rijstenbil, J.W.; Derksen, J.W.M.; Gerringa, L.J.A.; Poortvliet, T.C.W.; Sandee, A.; Berg, M.; Drie, J.; Wijnholds, J.A. Oxidative stress induced by copper: Defense and damage in the marine planktonic diatom Ditylum brightwellii, grown in continuous cultures with high and low zinc levels. Mar. Biol. 1994, 119, 583–590. [Google Scholar] [CrossRef] [Green Version]
  30. Buhmann, M.T.; Schulze, B.; Foerderer, A.; Schleheck, D.; Kroth, P.G. Bacteria may induce the secretion of mucin-like proteins by the diatom Phaeodactylum tricornutum. J. Phycol. 2016, 52, 463–474. [Google Scholar] [CrossRef] [Green Version]
  31. Willis, A.; Chiovitti, A.; Dugdale, T.M.; Wetherbee, R. Characterization of the extracellular matrix of Phaeodactylum tricornutum (Bacillariophyceae): Structure, composition, and adhesive characteristics. J. Phycol. 2013, 49, 937–949. [Google Scholar] [CrossRef]
  32. Tong, C.Y.; Derek, C.J.C. The role of substrates towards marine diatom Cylindrotheca fusiformis adhesion and biofilm development. J. Appl. Phycol 2021, 33, 2845–2862. [Google Scholar] [CrossRef]
  33. Tong, C.Y.; Derek, C.J.C. Biofilm formation of benthic diatoms on commercial polyvinylidene fluoride membrane. Algal Res. 2021, 55, 102260. [Google Scholar] [CrossRef]
Figure 1. Growth of P. tricornutum (a) and C. fusiformis (b) under different Cu2+ or Zn2+ concentrations (0, 30 and 60 μM). Data points are the means of triplicates, and error bars represent the standard deviation.
Figure 1. Growth of P. tricornutum (a) and C. fusiformis (b) under different Cu2+ or Zn2+ concentrations (0, 30 and 60 μM). Data points are the means of triplicates, and error bars represent the standard deviation.
Water 14 03305 g001
Figure 2. Cell morphology of P. tricornutum (ac) and C. fusiformis (df) under different Cu2+ or Zn2+ concentrations (control, 5 μM Cu2+, and 30 μM Zn2+) recorded on TM4000 Plus SEM. The scale represents 10 μm.
Figure 2. Cell morphology of P. tricornutum (ac) and C. fusiformis (df) under different Cu2+ or Zn2+ concentrations (control, 5 μM Cu2+, and 30 μM Zn2+) recorded on TM4000 Plus SEM. The scale represents 10 μm.
Water 14 03305 g002
Figure 3. EDS analysis of P. tricornutum (ac) and C. fusiformis (df) under different Cu2+ or Zn2+ concentrations (control, 5 μM Cu2+, and 30 μM Zn2+) recorded on TM4-EDS. The box in the figure indicates the area scanned by EDS. The data in the figure reflects the percentage of Si, Cu, and Zn elements.
Figure 3. EDS analysis of P. tricornutum (ac) and C. fusiformis (df) under different Cu2+ or Zn2+ concentrations (control, 5 μM Cu2+, and 30 μM Zn2+) recorded on TM4-EDS. The box in the figure indicates the area scanned by EDS. The data in the figure reflects the percentage of Si, Cu, and Zn elements.
Water 14 03305 g003
Figure 4. Length distribution of transcripts in P. tricornutum.
Figure 4. Length distribution of transcripts in P. tricornutum.
Water 14 03305 g004
Figure 5. Functional annotation of unigenes in P. tricornutum.
Figure 5. Functional annotation of unigenes in P. tricornutum.
Water 14 03305 g005
Figure 6. Gene Ontology (GO) enrichment analysis of differentially expressed genes (DEGs) in P. tricornutum under Cu2+ treatment. (a) Upregulated genes in PTCu/PTC. (b) Downregulated genes in PTCu/PTC.
Figure 6. Gene Ontology (GO) enrichment analysis of differentially expressed genes (DEGs) in P. tricornutum under Cu2+ treatment. (a) Upregulated genes in PTCu/PTC. (b) Downregulated genes in PTCu/PTC.
Water 14 03305 g006
Figure 7. Gene Ontology (GO) enrichment analysis of differentially expressed genes (DEGs) in P. tricornutum under Zn2+ treatment. (a) Upregulated genes in PTZn/PTC. (b) Downregulated genes in PTZn/PTC.
Figure 7. Gene Ontology (GO) enrichment analysis of differentially expressed genes (DEGs) in P. tricornutum under Zn2+ treatment. (a) Upregulated genes in PTZn/PTC. (b) Downregulated genes in PTZn/PTC.
Water 14 03305 g007
Figure 8. Length distribution of transcripts in C. fusiformis.
Figure 8. Length distribution of transcripts in C. fusiformis.
Water 14 03305 g008
Figure 9. Functional annotation of unigenes in C. fusiformis.
Figure 9. Functional annotation of unigenes in C. fusiformis.
Water 14 03305 g009
Figure 10. Gene Ontology (GO) enrichment analysis of differentially expressed genes (DEGs) in C. fusiformis under Cu2+ treatment. (a) Upregulated genes in CFCu/CFC. (b) Downregulated genes in CFCu/CFC.
Figure 10. Gene Ontology (GO) enrichment analysis of differentially expressed genes (DEGs) in C. fusiformis under Cu2+ treatment. (a) Upregulated genes in CFCu/CFC. (b) Downregulated genes in CFCu/CFC.
Water 14 03305 g010
Figure 11. Gene Ontology (GO) enrichment analysis of differentially expressed genes (DEGs) in C. fusiformis under Zn2+ treatment. (a) Upregulated genes in CFZn/CFC. (b) Downregulated genes in CFZn/CFC.
Figure 11. Gene Ontology (GO) enrichment analysis of differentially expressed genes (DEGs) in C. fusiformis under Zn2+ treatment. (a) Upregulated genes in CFZn/CFC. (b) Downregulated genes in CFZn/CFC.
Water 14 03305 g011
Table 1. DEGs involved in heavy metal stress response in P. tricornutum under Cu2+.
Table 1. DEGs involved in heavy metal stress response in P. tricornutum under Cu2+.
Gene_idfcRegulatenrPathsSwissprot
Pt04g035505.7upXP_002181744.1 (predicted protein)map00480 (Glutathione metabolism); map00053 (Ascorbate and aldarate metabolism)Probable L-ascorbate peroxidase 8
Pt05g022602.1upXP_002186090.1 (catalase-peroxidase)map00360
(Phenylalanine metabolism); map00380 (Tryptophan metabolism)
Catalase-peroxidase
Pt08g0213030.8upXP_002179007.1 (predicted protein)map00480 (Glutathione metabolism)Probable cytosol aminopeptidase
Pt14g009800.3downXP_002181057.1 (predicted protein)
Pt02g055503.3upXP_002177701.1 (predicted protein)map00480 (Glutathione metabolism)
Pt03g031500.4downXP_002185216.1 (glyoxalase)map00620 (Pyruvate metabolism)Hydroxyacylglutathione hydrolase
Pt02g039602.8upXP_002177790.1 (predicted protein)
Pt07g010500.4downXP_002185856.1 (predicted protein)map00620 (Pyruvate metabolism)
Pt15g026909.6upXP_002182163.1 (predicted protein)
Pt12g009302.7upXP_002180005.1 (predicted protein)map00480 (Glutathione metabolism)Glutathione S-transferase DHAR2
Pt01g092000.5downXP_002177254.1 (predicted protein, partial) Glutathione gamma-glutamylcysteinyltransferase
Pt11g0190010.3upXP_002182079.1 (predicted protein)
Pt14g036502.1upXP_002180739.1 (glutathione peroxidase, partial)map00590 (Arachidonic acid metabolism); map00480 (Glutathione metabolism)Phospholipid hydroperoxide glutathione peroxidase
Pt05g024702.2upXP_002186390.1 (predicted protein)
Pt11g010905.9upXP_002182079.1 (predicted protein)
Pt07g041705.1upXP_002176312.1 (peroxidase domain-containing protein) Putative heme-binding peroxidase
Pt11g0313017.5upXP_002181851.1 (predicted protein)
Pt14g036502.1upXP_002180739.1 (glutathione peroxidase, partial)map00590 (Arachidonic acid metabolism); map00480 (Glutathione metabolism)Phospholipid hydroperoxide glutathione peroxidase
Pt21g012200.3downXP_002183862.1 (predicted protein)
Table 2. DEGs involved in heavy metal stress response in P. tricornutum under Zn2+.
Table 2. DEGs involved in heavy metal stress response in P. tricornutum under Zn2+.
Gene_idfcRegulatenrPathsSwissprot
Pt04g035506.91upXP_002181744.1 (predicted protein)map00480 (Glutathione metabolism); map00053 (Ascorbate and aldarate metabolism)Probable L-ascorbate peroxidase 8
Pt20g016502.55upXP_002182954.1 (catalase)map00630 (Glyoxylate and dicarboxylate metabolism); map00380 (Tryptophan metabolism); map04146 (Peroxisome)Catalase
Pt13g019100.50downXP_002180671.1 (predicted protein)
Pt23g002202.22upXP_002184868.1 (predicted protein) Peroxiredoxin-6
Pt10g010300.03downXP_002179508.1 (predicted protein)map00480 (Glutathione metabolism)Glutathione S-transferase
Pt05g042809.78upXP_002186195.1 (UDP-glucose 6-dehydrogenase)map00520 (Amino sugar and nucleotide sugar metabolism); map00040 (Pentose and glucuronate interconversions); map00053 (Ascorbate and aldarate metabolism)UDP-glucose 6-dehydrogenase 1
Pt14g012700.24downXP_002180872.1 (l-ascorbate peroxidase, partial)map00480 (Glutathione metabolism); map00053 (Ascorbate and aldarate metabolism)Putative heme-binding peroxidase
Pt16g008806.11upXP_002179589.1 (nad-dependent epimerase/dehydratase)map00520 (Amino sugar and nucleotide sugar metabolism); map00053 (Ascorbate and aldarate metabolism)GDP-mannose 3,5-epimerase
Pt08g031905.92upXP_002178726.1 (predicted protein)map00460 (Cyanoamino acid metabolism); map00480 (Glutathione metabolism); map00430 (Taurine and hypotaurine metabolism)Glutathione hydrolase-like YwrD proenzyme
Pt04g015100.39downXP_002183098.1 (glutathione peroxidase domain-containing protein)map00590 (Arachidonic acid metabolism); map00480 (Glutathione metabolism)Probable phospholipid hydroperoxide glutathione peroxidase
Pt02g039602.85upXP_002177790.1 (predicted protein)
Pt12g009304.58upXP_002180005.1 (predicted protein)map00480 (Glutathione metabolism)Glutathione S-transferase DHAR2
Pt21g0220025.23upXP_002183815.1 (predicted protein)map02010 (ABC transporters)Glutathione-binding protein GsiB
Pt18g021906.06upXP_002185391.1 (predicted protein)
Pt05g024702.74upXP_002186390.1 (predicted protein)
Pt23g011508.11upXP_002184892.1 (predicted protein) Glutathione gamma-glutamylcysteinyltransferase 2
Pt08g027300.47downGAX19067.1 (hypothetical protein FisN_8Hh293 [Fistulifera solaris]) ABC transporter G family member 1
Pt12g031600.48downXP_002180322.1 (glutathione reductase)map00480 (Glutathione metabolism)Glutathione reductase
PtUn01s1135.37upXP_002177253.1 (mutase superoxide dismutase)map04146 (Peroxisome)Superoxide dismutase
Pt13g0293012.70upXP_002180497.1 (precursor of mutase superoxide dismutase [Fe/Mn], partial)map04146 (Peroxisome)Superoxide dismutase
Pt01g091907.47upXP_002177253.1 (mutase superoxide dismutasemap04146 (Peroxisome)Superoxide dismutase
Pt05g044700.40downXP_002186201.1 (5′-Nucleotidase or metallophosphoesterase)
Pt07g041702.73upXP_002176312.1 (peroxidase domain-containing protein) Putative heme-binding peroxidase
Pt20g012202.27upXP_002182845.1 (predicted protein)map04146 (Peroxisome)Peroxiredoxin-2C
Table 3. DEGs involved in heavy metal stress response in C. fusiformis under Cu2+.
Table 3. DEGs involved in heavy metal stress response in C. fusiformis under Cu2+.
Gene_idnr_DescriptionfcRegulatePathsSwissprot
TRINITY_DN14518_c0_g1thioredoxin-like protein2.68upmap00940 (Phenylpropanoid biosynthesis)1-Cys peroxiredoxin A
TRINITY_DN1479_c0_g1glutathione synthetase0.28downmap00270 (Cysteine and methionine metabolism); map00480 (Glutathione metabolism)Glutathione synthetase
TRINITY_DN495_c0_g2hypothetical protein0.19downmap00480 (Glutathione metabolism)Glutathione S-transferase
TRINITY_DN7366_c1_g1glutathione peroxidase0.43downmap00590 (Arachidonic acid metabolism); map00480 (Glutathione metabolism)Hydroperoxy fatty acid reductase gpx1
TRINITY_DN1758_c0_g1hydroxyacylglutathione hydrolase0.45downmap00620 (Pyruvate metabolism); map00790 (Folate biosynthesis)Hydroxyacylglutathione hydrolase
TRINITY_DN3215_c0_g1hypothetical protein0.46down
TRINITY_DN2680_c0_g1hypothetical protein0.49downmap00270 (Cysteine and methionine metabolism); map00480 (Glutathione metabolism)Glutamate--cysteine ligase catalytic subunit
TRINITY_DN6304_c0_g1hypothetical protein0.38down
Table 4. DEGs involved in heavy metal stress response in C. fusiformis under Zn2+.
Table 4. DEGs involved in heavy metal stress response in C. fusiformis under Zn2+.
Gene_idnr_DescriptionfcSignificantRegulatePathsSwissprot
TRINITY_DN14518_c0_g1thioredoxin-like protein3.17yesupmap00940 (Phenylpropanoid biosynthesis)1-Cys peroxiredoxin A
TRINITY_DN1479_c0_g1glutathione synthetase0.33yesdownmap00270 (Cysteine and methionine metabolism); map00480 (Glutathione metabolism)Glutathione synthetase
TRINITY_DN1711_c0_g1hypothetical protein2.81yesupmap00590 (Arachidonic acid metabolism); map00480 (Glutathione metabolism)Glutathione S-transferase
TRINITY_DN1711_c0_g2hypothetical protein6.12yesupmap00590 (Arachidonic acid metabolism); map00480 (Glutathione metabolism)Glutathione S-transferase 1
TRINITY_DN2013_c0_g1glutathione-S-transferase0.44yesdownmap00590 (Arachidonic acid metabolism); map00480 (Glutathione metabolism)Glutathione S-transferase
TRINITY_DN327_c0_g2glutathione S-transferase2.73yesupmap00590 (Arachidonic acid metabolism); map00480 (Glutathione metabolism)Glutathione S-transferase 1
TRINITY_DN6449_c0_g1hypothetical protein2.06yesup Glutathionyl-hydroquinone reductase
TRINITY_DN7366_c1_g1glutathione peroxidase0.37yesdownmap00590 (Arachidonic acid metabolism); map00480 (Glutathione metabolism)Hydroperoxy fatty acid reductase
TRINITY_DN2338_c0_g3hypothetical protein0.40yesdown DEP domain-containing mTOR-interacting protein
TRINITY_DN3_c0_g4mercuric reductase0.43yesdown General L-amino acid-binding periplasmic protein Aap
TRINITY_DN17353_c0_g1LhcSR2.98yesupmap00196 (Photosynthesis—Antenna proteins)Light-harvesting complex stress-related protein
TRINITY_DN1775_c0_g1hypothetical protein0.30yesdown Thyroid peroxidase
TRINITY_DN319_c0_g1oxidative stress-related Abc1-like protein2.10yesup Protein ACTIVITY OF BC1 COMPLEX KINASE 8
TRINITY_DN3894_c1_g1catalase peroxidase2.40yesupmap00940 (Phenylpropanoid biosynthesis); map00380 (Tryptophan metabolism); map00360 (Phenylalanine metabolism)Catalase-peroxidase
TRINITY_DN5279_c0_g2methionine sulfoxide reductase B0.45yesdown Peptide methionine sulfoxide reductase
TRINITY_DN6304_c0_g1hypothetical protein0.39yesdown Peroxinectin A
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Huang, A.; Wang, Y.; Duan, J.; Guo, S.; Xie, Z. Differential Response of Phaeodactylum tricornutum and Cylindrotheca fusiformis to High Concentrations of Cu2+ and Zn2+. Water 2022, 14, 3305. https://doi.org/10.3390/w14203305

AMA Style

Huang A, Wang Y, Duan J, Guo S, Xie Z. Differential Response of Phaeodactylum tricornutum and Cylindrotheca fusiformis to High Concentrations of Cu2+ and Zn2+. Water. 2022; 14(20):3305. https://doi.org/10.3390/w14203305

Chicago/Turabian Style

Huang, Aiyou, Yujue Wang, Jiawen Duan, Shiyi Guo, and Zhenyu Xie. 2022. "Differential Response of Phaeodactylum tricornutum and Cylindrotheca fusiformis to High Concentrations of Cu2+ and Zn2+" Water 14, no. 20: 3305. https://doi.org/10.3390/w14203305

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop