Next Article in Journal
Prediction of Rainfall Time Series Using the Hybrid DWT-SVR-Prophet Model
Previous Article in Journal
Performance Analysis of Anaerobic Digestion Coupled with Simultaneous Nitrification and Denitrification Process for Treating Alcohol Precipitation Wastewater of Chinese Patent Medicine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of the Effects of Different Phosphorus Sources on Microcystis aeruginosa Growth and Microcystin Production via Transcriptomic Surveys

Department of Environmental Science and Engineering, Fudan University, Shanghai 200433, China
*
Author to whom correspondence should be addressed.
Water 2023, 15(10), 1938; https://doi.org/10.3390/w15101938
Submission received: 13 April 2023 / Revised: 17 May 2023 / Accepted: 17 May 2023 / Published: 20 May 2023
(This article belongs to the Section Water Quality and Contamination)

Abstract

:
Microcystis aeruginosa (M. aeruginosa) is a dominant species among frequent cyanobacterial blooms and is well known for its toxin-producing ability. Phosphorus (P) is a typical growth-limiting element for M. aeruginosa. Although alterations in physiological reactions have been documented in response to various P sources, the underlying molecular processes and transcriptional patterns remain poorly understood. This study evaluated the physiological and molecular responses of M. aeruginosa to different P sources. The growth of M. aeruginosa was promoted by both dissolved inorganic phosphorus (DIP) and dissolved organic phosphorus (DOP) at a concentration of 0.4 mg/L with an initial cell density of 1.50 (±0.05) × 106 cells/mL. The cell density reached 1.38 (±0.05) × 107 cells/mL in the DIP group on day 14, a value which was higher than that in the DOP group. Most photosynthesis genes had higher levels of upregulated expression in the DIP group. For instance, gene psbA was upregulated by 0.45 Log2Fold Change (Log2FC). In the DOP group, it is interesting that the Pi (PO4-P) concentration increased to 0.09 mg/L on day 14. Meanwhile, the expression of the gene encoding alkaline phosphatase-like protein was significantly upregulated, with a value of 1.56 Log2FC, and the alkaline phosphatase concentration increased in the DOP group. The intracellular microcystin (IMC) concentration decreased with time in both groups. However, the concentration of extracellular microcystins (EMCs) increased with incubation time in both groups. Phosphorus participates in the regulation of microcystin synthesis, mainly by regulating ATP synthesis. Based on the physiological and molecular investigations in this study, the results provide crucial insights into the physiological adaptations and the role of P in modulating harmful algal bloom formation, microcystin synthesis, and potential molecular responses of M. aeruginosa.

1. Introduction

Eutrophication is induced by nitrogen and phosphorus loadings in aquatic systems. High levels of eutrophication in aquatic systems are typically associated with harmful algal blooms (HABs). HABs are rapid growths of cyanobacteria that can cause remarkable damage to aquatic ecosystems [1,2]. Accordingly, their proliferation has become a global environmental issue and has gained considerable attention [3]. Microcystis aeruginosa (M. aeruginosa), one of the most widespread bloom-forming cyanobacteria, mainly develops in organic-rich freshwater bodies such as lakes and ponds. Microcystins (MCs), which are cyclic heptapeptide compounds generated by M. aeruginosa, are extremely hepatotoxic [4,5]. Controlling HABs and toxin generation have been priorities of research globally, with critical implications for human health.
Prior studies on the effects of nitrogen on M. aeruginosa have been relatively comprehensive and have covered various forms of nitrogen, nitrogen starvation, and nitrogen sources at different concentrations [6,7]. Additionally, the effects of these on an organism’s physiological state and transcriptome level have been deemed relatively sufficient [7,8]. Moreover, a previous study reported that M. aeruginosa exhibited the highest growth rate and maximal optical density when cultured in an N:P ratio of 0.1 [9], indicating that the growth of M. aeruginosa may be stimulated by an excessive amount of nitrogen and a consistent level of phosphorus in the culture medium. Phosphorus (P) is an essential nutrient for the development and metabolism of cyanobacteria as it is required for the formation of ribosomes, DNA, and cell membranes [10,11]. Furthermore, the availability of phosphorus has an important impact on cyanobacterial growth and metabolism and is regarded as a limiting factor for the occurrence of cyanobacterial blooms [12,13]. To date, most investigations, ranging from field to gene expression studies, have focused on the role of phosphorus in promoting cyanobacterial blooms and toxicity. Accordingly, physiological changes and MC production based on molecular processes through transcriptomic analysis in response to different phosphorus sources have not been well studied.
Among all the phosphorus compounds in natural waters, dissolved inorganic phosphorus (DIP) and dissolved organic phosphorus (DOP) are the dominant forms [14,15,16]. Zhang et al. discovered that the growth of M. aeruginosa was enhanced by elevated concentrations of DIP, but inhibited at high concentrations (0.6 and 1.0 mg P/L) of DOP after P starvation. The high concentration of DIP was found to promote the production of MCs in M. aeruginosa, while the high concentration of DOP induced the release of intracellular MCs without affecting MC production [17]. DIP was inferred to be a more bioavailable phosphorus source because M. aeruginosa growth was promoted more significantly by DIP than by DOP. Cyanobacteria always bloom during the summer with a large uptake of phosphorus, thereby leading to DIP depletion in freshwater [18,19,20]. The DIP content in surface water is extremely low, causing widespread phosphorus stress or even starvation [10,20]. To overcome P stress or starvation during this period, M. aeruginosa may shift to assimilate DOP to survive and remain the dominant species in natural waters. However, these conditions may trigger physiological variations and affect toxin production mechanisms. Studies regarding the cellular responses in M. aeruginosa under the condition that DOP serves as the sole P source may help evaluate toxin production and toxin persistence in cyanobacterial blooms. Moreover, further assessments can provide basic knowledge that will improve the phosphorus control of cyanobacteria in natural waters.
Numerous studies have determined the effects of various phosphorus sources on cyanobacteria. Li et al. evaluated the effects of three different phosphorus compounds, sodium-β-glycerophosphate (β-gly), adenosine triphosphate (ATP), and orthophosphate (K2HPO4), on M. aeruginosa [21]. The outcomes of their study demonstrated that M. aeruginosa can utilize both DIP and DOP. However, the intracellular responses of M. aeruginosa to different forms of P remain unknown. Fluorescence quantitative PCR was performed to identify those transcripts of the target genes related to P transportation (sphX, pstS, and phoX) and MC synthesis (mcyA and mcyD) and in which the pstS gene was found to be expressed with the highest abundance in the lowest concentration (0.02 mg/L) in the DOP group [17]. The effects of DIP and DOP on phosphorus uptake and toxin synthesis by M. aeruginosa were examined at the molecular level. However, understanding the expressions of several of these genes is insufficient because cyanobacterial metabolism is regulated by multiple genes. To fully comprehend how M. aeruginosa responds to DIP and DOP, the genes that are differentially expressed must be identified. Transcriptomics enables the investigation of gene transcription and gene expression at the RNA level. The response of cyanobacteria to different phosphorus sources can be more accurately assessed if genetic changes are combined with physiological changes. To develop innovative strategies to control cyanobacterial blooms, the growth, P metabolism, and toxicity of M. aeruginosa under DIP and DOP must be understood.
This study investigated the response of M. aeruginosa to different phosphorus sources through transcriptome analysis. Incubations were performed to evaluate the development of M. aeruginosa on DIP and DOP substrates. Furthermore, the changes in P concentrations were monitored, and the variations in P utilization were compared. The transcriptome data were compared with the physiological evaluation data of cell density, chlorophyll-a (chl-a), photosynthesis, MC production, and alkaline phosphatase (AP) concentration. This study aimed to (1) determine the effects of two different phosphorus sources on the growth of M. aeruginosa, (2) reveal the molecular mechanism of P utilization, MC production, and AP activity in M. aeruginosa, and (3) improve our understanding of the relationship between MC production and P. This study elucidates the physiological responses of M. aeruginosa to different phosphorus sources at the molecular level.

2. Materials and Methods

2.1. Microcystis aeruginosa and Cultures

M. aeruginosa (FACHB-905), isolated from cyanobacterial blooms in Dianchi Lake in China, was obtained from the Institute of Freshwater Algae Culture Collection at the Institute of Hydrobiology, Chinese Academy of Sciences (Wuhan, China). M. aeruginosa was cultured and grown in BG11 medium [22]. All solutions were prepared with deionized water. M. aeruginosa was cultivated in a climatic chamber at 25 °C with a light intensity of 27 µmol m−2 s−1 on a 12:12 light/dark cycle (BIC-300, Boxun, Shanghai, China). All flasks were gently shaken by hand three times daily to avoid cell aggregation or settlement [7].

2.2. Experimental Design

M. aeruginosa cells in the exponential growth phase were collected via centrifugation (3800× g for 10 min), washed three times with sterilized deionized water, and inoculated in the sterilized phosphorus-free BG-11 medium for 7 days to exhaust intracellular phosphate. Subsequently, the P-free medium was divided into two different groups: (1) one that used dipotassium phosphate medium (K2HPO4) as DIP and (2) one that used disodium beta-glycerophosphate pentahydrate medium (β-gly, containing C-O-P bonds) as DOP, with both containing 0.4 mg/L P, the maximum total phosphorus concentration required for Class V water according to China’s surface water environmental quality standards [23]. The control group was M. aeruginosa medium on Day 0. The initial cell biomass was adjusted to 1.5 (±0.05) × 106 cells/mL, and M. aeruginosa was inoculated in 1 L Erlenmeyer flasks. All experiments were done in triplicates.
The experimental period lasted for 14 days. The cell density, chl-a concentration photosynthetic activity, and concentration of total phosphorus (TP) and PO4-P (Pi) were measured every day at the same time until day 8 and then measured every other day, whereas MCs and AP were analyzed every two days.

2.3. Physiological Measurements

2.3.1. Measurement of Cell Density and Growth Rate

A UV-vis spectrophotometer (L6S, Leng-Guang, Shanghai, China) was used to measure the optical density at 680 nm (OD680) [24], which was found to be positively associated with cell numbers. The specific growth rate (μ) was calculated as per Equation (1):
μ = (lnN2 − lnN1)/(t2 − t1)
where N1 and N2 are the cell densities at times t1 and t2, respectively.

2.3.2. Extraction and Analysis of Chlorophyll a

10 mL cell culture was collected and centrifuged at 3600× g for 10 min, and the supernatant was discarded. The cells were rinsed with 5 mL of 90% acetone and left overnight at 4 °C in darkness. Then, the extracts were filtered with a 0.45 μm filter membrane, and absorbances at 630, 645, 663, and 750 nm were measured using a UV–vis spectrophotometer (L6S, Leng-Guang, Shanghai, China) to calculate the content of chl-a according to Equations (2) and (3), with a 90% acetone solution used as the blank [25].
chl-a (mg/L) = [11.64 × (OD663 − OD750) − 2.16 × (OD645 − OD750) +
0.10 × (OD630 − OD750)]·V1V−1 δ−1
chl-a (pg/cell) = (2) × 109/cell density
Here, V1 is the constant volume of the extracted cultures (mL), V denotes the sample volume, and δ is the optical path of the cuvette.

2.3.3. Analysis of Photosynthetic Activity

To evaluate photosynthetic activity, the Fv/Fm ratio was calculated. A 4-mL portion of the cell culture was dark-incubated for 15 min. The Fluorometer (AquaPen AP 100, Aqualabo, Brno, Czech Republic) was used to monitor the cells’ photosynthetic activity by measuring the variable-to-maximum chlorophyll fluorescence ratio (Fv/Fm ratio).

2.3.4. Measurement of Total Phosphorus and PO4-P

Subsamples of the DIP and DOP were filtered with a 0.45 μm filter. According to the instructions for the use of the Hach total phosphorus reagent, 5 mL of supernatants were added in the reagent tube, then potassium persulfate was added and digested at 150 °C for 30 min. After the mix was cooled to room temperature, 2 mL of NaOH and a packet of phosver3 phosphate reagent were added to the tube. The sample was used to determine the levels of TP by using the ammonium molybdate spectrophotometric method, and the undigested supernatants were used to determine the levels of Pi by using the Mo-Sb Anti-Spectrophotometer in the culture [26]. TP and Pi concentration were measured with an ultraviolet-visible spectrophotometer (HACH, DR6000).

2.3.5. Extraction and Measurement of Microcystins

The samples were centrifugated at 3800× g for 10 min. After centrifugation, the supernatant was collected to assess extracellular microcystin (EMC) levels, and the cyanobacterial cells were utilized in order to identify the intracellular microcystins (IMCs). The cell pellet was frozen at −20 °C for 24 h [7]. The IMCs were extracted by using the freeze-thaw method, wherein the cyanobacterial cells were freeze-thawed three times [17]. MC concentrations were evaluated using an enzyme-linked immunosorbent assay (ELISA) kit (Product No.062964, Shanghai MLBIO Biotechnology Co., Ltd., Shanghai, China) and detected through Microplate Reader–Spectra Max M2 (Molecular Devices, 3860 N First Street, San Jose, CA, Canada). The required operating methods were followed exactly as instructed by the kit [27].

2.3.6. Extraction and Analysis of Alkaline Phosphatase

The supernatant was collected to determine AP concentration via an enzyme-linked immunosorbent assay (ELISA) kit (Product No. 062906, Shanghai MLBIO Biotechnology Co., Ltd., Shanghai, China) and also detected through Microplate Reader–Spectra Max M2 (Molecular Devices, 3860 N First Street San Jose, CA, Canada). The absorbance was measured at a 450-nm wavelength, and the content of AP in the sample was calculated by using a standard curve.

2.4. RNA Extraction and Sequencing

For transcriptomic analysis, the samples were collected in 30 mL centrifuge tubes and rapidly harvested via centrifugation at 2500× g for 10 min at 4 °C. The samples were immediately frozen in liquid N after harvesting and stored at −80 °C until sequencing. All sample procedures were conducted within 30 min. To minimize temporal variations in gene expression and cell physiology, cultures were collected at the same time every day for all groups. All groups were analyzed in triplicate [28]. Total RNA was extracted using TRIzol reagent (Invitrogen Life Technologies, Carlsbad, CA, USA). The mRNA was isolated using a Ribo-Zero rRNA Removal Kit (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions. A NanoDrop spectrophotometer (NC2000, Thermo Scientific, Waltham, MA, USA) and a Bioanalyzer 2100 system (Agilent, Santa Clara, CA, USA) were used to assess RNA quality and integrity. The abundances of the transcripts were measured in terms of RPKM (Reads Per Kilobase per Million reads) for the analysis of the transcriptome data. DESeq (version 1.18.0) tools were used to perform a differential expression analysis for comparisons between samples with three biological replicates. The metabolic pathways were mapped using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The Microcystis DEGs were mapped to the NCBI (National Center for Biotechnology Information) genome database for KEGG analysis.

3. Results and Discussion

3.1. Growth and Photosynthetic Activity of Microcystis aeruginosa

3.1.1. Physiological Changes in Growth and Photosynthesis

Cell density and chl-a concentration act as indicators of growth, while the Fv/Fm ratio is used to evaluate the photosynthetic activity of M. aeruginosa. As shown in Figure 1A,B, M. aeruginosa growth was promoted under both phosphorus conditions, but the cell density was higher in the DIP group than in the DOP group. In the DIP group, the specific growth rate remained high in the first 10 days and then steadily decreased over the next several days. As shown in Figure 1B, the specific growth rate in the DIP group was higher than that in the DOP group after day 4, consistent with the higher cell density in the DIP group. The maximum cell density reached 1.38 (±0.05) × 107 cells/mL by the end of the experiment. In the DOP group, M. aeruginosa maintained a stable specific growth rate, and the cell density on the last day was 0.74-fold that in the DIP group. This result is consistent with that reported by Yuan et al. [29], who showed that M. aeruginosa grows better in the presence of K2HPO4 as the phosphorus source than in the presence of β-gly and ATP. However, a previous study found that the maximal cell density of M. aeruginosa in DOP groups was only 15% of that in the DIP groups under an initial concentration of 0.6 mg P/L [17]. This apparent difference may be attributed to the different initial phosphorus concentration settings and experimental periods. DIP uptake stimulates the growth and cell division of M. aeruginosa, leading to the formation of large blooms. In contrast, DOP uptake stimulates the production of extracellular enzymes, which can cleave DOP into inorganic forms. This can increase the availability of DIP and promote the further growth of M. aeruginosa [30].
Chl-a concentrations are a key marker of photosynthetic activity in M. aeruginosa [6]. As shown in Figure 1C, the chl-a content increased significantly from day 0 to day 10 and then showed a slight increase in both groups, peaking at 0.62 pg/cell on day 10 in the DIP group and 0.57 pg/cell on day 10 in the DOP group. The change in chl-a content displayed a similar trend to that of cell density, aligning with previous findings [7]. These results indicate that both DIP and DOP can stimulate the photosynthetic activity of M. aeruginosa. It is worth noting that the increase in chl-a content was slightly higher in the DIP group compared with that in the DOP group; this may be attributed to the higher bioavailability of DIP.
Fv/Fm values are positively correlated with photosystem (PS) II performance, which acts as a photosynthetic efficiency parameter [31]. As illustrated in Figure 1D, compared to the initial Fv/Fm value (0.12), the maximum values corresponded to 0.40 (DIP) and 0.36 (DOP), showing a significant increase in photosynthetic efficiency in both groups. Previous studies reported that the Fv/Fm values in DIP (K2HPO4) groups were higher than those in DOP (β-gly) groups under the condition of 0.2 mg P/L [17], in alignment with our findings, indicating higher cell density in DIP groups. The Fv/Fm values of both groups increased and showed gradual declines until the end of cultivation, thereby aiding cyanobacterial growth and proliferation.
Our results showed that the growth changes in cell density displayed similar patterns in both DIP and DOP groups, and that the growth rate first increased and then decreased. Besides, when DOP was used as the sole P source, the chl-a contents and the Fv/Fm ratio were marginally lower than they were in the DIP group, indicating that DOP was bioavailable to M. aeruginosa and could help M. aeruginosa alleviate DIP depletion.

3.1.2. Changes in Photosynthetic Genes’ Response to Different Phosphorus Sources

The transcriptional abundance of most genes involved in photosynthesis was variably expressed after the addition of different phosphorus sources. According to the literature, both PS I and PS II take part in the conversion of light energy into the energy of chemical bonds [32]. As shown in Figure 2A, compared to what was observed in the control, most PS I genes in the DIP group showed significant changes, and the expression of psaC was markedly elevated by 1.39 Log2Fold Change (Log2FC). PS I gene expression in the DOP group exhibited a similar pattern, with an elevated expression range of 0.13 to 0.38 Log2FC. Most of the detected genes involved in PS II, such as psbC—which encodes the PS II chlorophyll-binding protein CP43—and psbU—which encodes the PS II 12 kDa extrinsic protein [33]—had upregulated transcript abundance in the DIP group. However, other genes, such as psbY, psbZ, psb27, and psb28, exhibited a reduction in relative abundance. For example, psb28, one of the assembly factors of PS II [34], was downregulated by 0.19 Log2FC on day 5. Moreover, approximately half of the PS II-related genes were upregulated in the DOP group: a number fewer than that of the PS II-related genes in the DIP group. This result was consistent with the value of Fv/Fm in the DOP group (Figure 1D), which was lower than that in the DIP group.
All genes encoding ATP synthase subunits showed reductions in expression, as shown in Figure 2A. In the DIP group, the expressions of atpA (encoding ATP synthase F1 and alpha subunit), atpB (encoding ATP synthase F1 and beta subunit), and atpI (encoding ATP synthase I family protein) were downregulated by 1.17, 1.49, and 1.56 Log2FC, respectively. Similarly, in the DOP group, the transcript abundances of these three genes (atpA, atpB, and atpI) were downregulated by 1.18, 1.64, and 1.63 Log2FC, respectively. These results indicate that the energies of photosynthesis in both DIP and DOP group were lower than that in the control. This might be due to the phosphorus supply in the DIP and DOP groups being severely depleted, resulting in decreased ATPase activity in both groups. The cyanobacteria in the DIP and DOP groups had both decreased their energy consumption by reducing their photosynthesis activity.
Cyanobacteria are the principal producers in freshwater environments and are considered ideal model organisms for photosynthesis studies owing to their high photosynthetic efficiency and rapid development rate. The photosynthetic apparatus of cyanobacteria mainly consists of PS II, cytochrome b6f, PS I, phycobilisomes, and ATP synthase [35], as shown in Figure 2B. In PS II, the gene psbA, which encodes the photosynthetic reaction center protein [36], was upregulated by 0.45 Log2FC in the DIP group. However, this gene was downregulated by 0.85 Log2FC in the DOP group, as shown in Figure 2A. A previous report revealed that the dinoflagellate Peridinium bipes, a typical component of the phytoplankton group, can form blooms in freshwater ecosystems, and that the expression of psbA in the DOP (β-gly, 52 µmol/L) groups was downregulated compared to that in the DIP (K2HPO4, 52 µmol/L) groups, ultimately aligning with the findings of our investigation [37]. In addition, the expressions of psbC, which encodes light-harvesting proteins in PS II [38], were similar between the DIP and DOP groups. Accordingly, rather than being connected with the activity of the light-harvesting protein complex, the PS II photoreactive central protein of cyanobacteria was more active in the DIP group, and may have resulted in greater cell development in the DIP group.
As shown in Figure 2B, PS I is the second PS involved in the photo-cooperative photoreaction of cyanobacteria. A full membrane protein complex called PS I uses light energy to produce the high-energy carriers ATP and NADPH [39]. The large subunits PsaA and PsaB, which are the core subunits in PS I, are located in the center of the PS I monomer, and encode the PS I P700 chlorophyll-a apoproteins A1 and A2, respectively. Moreover, most chlorophyll and carotenoids in the antenna system are found in subunits PsaA and PsaB [33]. Therefore, the photosynthetic activity of cyanobacteria was markedly influenced by the expression of psaA and psaB. In this experiment, both the DIP and DOP groups displayed considerably high expression of psaA, which was upregulated by 0.88 Log2FC in the DIP group and 0.36 Log2FC in the DOP group. Although both genes were upregulated, the fold change in the DIP group was greater than that in the DOP group. In addition, the upregulated expression of psaB in the DIP group was markedly higher than that in the DOP group. Overall, the cells in the DIP group grew substantially more quickly than those in the DOP group owing to variations in the expression of genes involved in photosynthesis.
Our results indicated that both phosphorus sources promoted the upregulated expression of photosynthetic genes in M. aeruginosa. Additionally, we found that the gene encoding the expression of the PS II photoreactive central protein was significantly differentially expressed, compared to other genes, between the DIP and the DOP groups, leading to the concentration of chl-a in the DIP group being higher than that in the DOP group.

3.2. Variations in Related Genes’ Responses to Dissolved Inorganic Phosphorus and Dissolved Organic Phosphorus

3.2.1. Utilization of Different Phosphorus Sources and Alkaline Phosphatase Activity

The assimilation of P by M. aeruginosa was indicated by the concentrations of TP and Pi in the culture medium. The initial TP concentration in both groups was 0.4 mg/L. As shown in Figure 3A, after six days of inoculation, the concentration declined to 0.08 mg/L and fluctuated slightly on subsequent days in the DIP group. In contrast to the rapid reduction in TP concentrations in the DIP group, TP concentrations declined slowly and subsequently remained stable over the study period when β-gly was employed as a P source. The decreasing TP concentrations in both DIP and DOP groups indicate that DIP and DOP can be utilized by M. aeruginosa in natural water to facilitate its growth [29]. Li et al. reported that after P starvation, when M. aeruginosa utilized DIP (K2HPO4) and other DOPs (β-gly, ATP), TP in the DIP group decreased sharply during the first 3 days of incubation and remained low regardless of the initial phosphorus concentration (0.1, 1.0, 2.0, or 4.0 mg/L) [21]. The results demonstrate that DIP was more rapidly assimilated than DOP and was used directly by M. aeruginosa.
Figure 3B shows the variations in Pi concentrations in both groups. Interestingly, a minor quantity of Pi was detected in the DOP group. The concentration of Pi grew steadily, reaching 0.09 mg/L at the end of cultivation. This result agreed with that of a previous study [17], which found that Pi concentrations increased to 0.03 and 0.11 mg P/L at 0.2 and 0.6 mg P/L in the DOP (β-gly) groups. The Pi concentration of the DOP (β-gly) groups was reported to peak at 9.217 mg/L on day 8 with Nostoc sp., in comparison to the concentration of 2 mg/L on the first day [40]. A possible reason is that AP, a critical enzyme of M. aeruginosa, converts complex organophosphates to phosphate; this may explain the observed increase in Pi concentration in DOP group over time [41].
Figure 4 shows the variations in AP concentrations. AP is crucial in order for cyanobacteria to utilize phosphate esters of all organic phosphorus sources, such as phosphate monoesters and phosphate diesters [42]. M. aeruginosa grown in the absence of P has been found to have slower growth rates and higher AP concentration [43]. As shown in Figure 4, AP was detected in the supernatant of the cyanobacterial culture after phosphorus deprivation for one week (0.8 ng/mL, day 0). When M. aeruginosa was inoculated in the DIP medium, the AP levels declined quickly to a low level on day 2 and remained stable for the following days at 0.15 ± 0.05 ng/mL. In contrast, the AP concentration increased at the beginning of cultivation in the DOP group. On day 6, the maximal concentration of AP in the β-gly groups reached 1.25 ng/mL, a value which was higher than that of the initial AP. This result suggests that M. aeruginosa may have a higher capacity to hydrolyze DOP compounds into Pi under conditions of DIP limitation. Therefore, DOP may represent an important phosphorus source for cyanobacterial growth in aquatic ecosystems, especially under DIP limitation.
The results demonstrated that DIP was a preferred P source for M. aeruginosa because the P was depleted faster in the DIP group. Moreover, we discovered that the AP concentration in the DIP group was much lower than that in the DOP group. The AP concentration increased slightly in the DOP group, possibly because M. aeruginosa required more AP, which might hydrolyze extracellular DOP and convert it into Pi for absorption and utilization by M. aeruginosa.

3.2.2. Variations in Related Genes Encoding Phosphorus Metabolism under Different Phosphorus Sources

Cyanobacteria thrive in a permanently changing aquatic environment, undergoing seasonal challenges and possibly overcoming inorganic phosphate deficiency. These conditions trigger physiological variations, through methods such as boosting the hydrolysis of organophosphorus by adjusting intracellular storage or using exogenous dissolved organophosphorus to activate multiple phosphorus usage pathways. AP is an enzyme that is essential for phosphorus metabolism because it releases phosphate from complex organophosphate compounds [44]. Figure 5A shows a heatmap of the genes related to P sensing, transportation, and AP. An unsymbolized gene that regulates an AP-like protein was discovered in this study, with its expression abundance exceeding 1.56 Log2FC in the DOP group. This result may indicate that this gene regulates the synthesis of substances involved in the hydrolysis of β-gly in order to be assimilated by M. aeruginosa. When DOP accounts for a significant proportion of the total dissolved P, the bioavailability of DOP to M. aeruginosa may serve as a potential reason that leads M. aeruginosa to become a dominant species [21,30].
Pst is an ATP-dependent ABC transporter-type system that includes phosphate-binding and transmembrane protein subunits [45]. As shown in Figure 5A, gene pstS is involved in phosphate ABC transporter substrate-binding protein, the expression of which was upregulated by 0.55 Log2FC in the DIP group. However, pstS expression was downregulated by 0.08 Log2FC in the DOP group. The transcript abundance of pstS has also been found to be lower in the β-gly groups than in the DIP groups in Raphidiopsis raciborskii [46]. PstA and PstC are transmembrane subunits that form a channel in the inner membrane. Both genes had similar expression levels in the DIP and DOP groups, indicating that the expression of these genes was not significantly affected by different phosphorus sources. Furthermore, the expression of pstB, which is involved in the phosphate transport system ATP-binding protein, was upregulated by 0.27 and 0.05 Log2FC in the DIP and DOP groups. In this study, the expression of related genes was generally upregulated when the Pst system responded to different phosphorus sources. The expression difference between pstS and pstB was more significant, indicating that phosphate-binding protein activity (PstS and PstB) was more affected than transmembrane protein activity (PstA and PstC) in the DIP group.
An overview of the phosphate transport mechanism is presented in Figure 5B, where four proteins (encoded by PstS, PstC, PstA, and PstB) constitute an ABC transporter in the Pst system required for high-affinity phosphate capture in the cytoplasm [47]. phoU, phoR, and phoB were found to play key roles in phosphate sensing and regulation, as shown in Figure 5B. PhoU, a regulatory protein that interacts with PstB and PhoB, regulates Pi import via the Pst system. The expression of phoU was upregulated by 0.19 and 0.07 Log2FC in the DIP and DOP groups, respectively. Previous studies revealed that the PhoU protein can sense a conformational change in the Pst system when Pi is imported [48], indicating that both phosphorus sources can be absorbed by cyanobacteria as Pi. PhoR is the Pi-sensing site, and its expression was downregulated by 1.50 and 1.09 Log2FC in the DIP and DOP groups, respectively. Moreover, the expression of the phosphate regulon response regulator phoB was downregulated by 0.43 Log2FC in the DIP group. The PhoB/R two-component system was induced by Pi starvation [49]. Therefore, the expression of these two genes in the DIP group was lower than that in the control group and appeared to be reasonable. However, the expression of phoB was upregulated by 0.04 Log2FC in the DOP group, and this might be due to the regulatory effect of phoU. In addition to coordinating Pi input, phoU may affect the expression of phoB [50], thereby affecting the expressions of part of the AP genes in PhoB regulators. The AP genes receive information on Pi stress from the Pi sensor PhoR or PhoB regulators and then hydrolyze the organophosphorus to Pi. Therefore, genes related to AP were obviously upregulated in the DOP group.
The corresponding hydrolysis and utilization pathways in cyanobacteria are shown in Figure 5B. In the presence of adequate phosphorus, cyanobacteria absorb extra phosphorus and store it in their cells. According to previous research on the cyanobacterial proteome and transcriptome, cyanobacterial cells store phosphorus as polyphosphate in phosphorus-rich environments [51]. Phosphorus is stored as polyphosphate by polyphosphatases (PPP), allowing cells to simultaneously employ phosphorus from polyphosphate and convert ATP to ADP.

3.3. Effects of Different Phosphorus Sources on Microcystins Production

The variations in IMC contents under different phosphorus sources are shown in Figure 6A. The concentrations of IMCs followed a similar trend in both the DIP and DOP groups, initially increasing and then decreasing. The greatest level of IMC was achieved on the second day, with 51.1 and 53.8 μg/L in the DIP and DOP groups, respectively. The initial increase in IMC content could be attributed to the stimulation of microcystin production by the availability of phosphorus. After day 8, IMC content in the two groups declined and remained stable. By day 10, the IMC content decreased to 23.9 μg/L in the DIP group. This result is consistent with that of Zhang et al. [17], who demonstrated that IMC concentrations in the DIP and DOP groups significantly decreased over time (on day 18 compared to day 8).
As shown in Figure 6B, the EMC content displayed an upward trend. The EMC contents of the DIP and DOP groups increased until they reached their peak values on day 14. The peak concentrations of EMC in the DIP and DOP groups were 10.34 and 9.68 µg/L, respectively. These values indicate 5.59- and 3.40-fold higher concentrations than the initial concentrations. Considering Figure 1A, it is clear that the changes in EMC concentration aligned with the cyanobacterial biomass variations, which increased with time in both groups. This conclusion aligns with that of earlier studies, which showed that EMC levels were compatible with the growth rates of cyanobacterial cells [17]. In addition, the lysis and the death of cyanobacterial cells cause the release of IMCs into the surrounding environment, which increases EMCs [52].

3.4. Relation between Microcystins Biosynthesis and Phosphorus

Previous research linked the synthesis of MCs to nitrogen in the environment [53] because nitrogen is a component of MCs. Although phosphorus is not a component of MCs, most studies have discovered that phosphorus considerably impacts M. aeruginosa proliferation and MC production [54,55].
Phosphorus is essential for trapping, storing, transferring, and distributing energy via high-energy phosphate bonds [56]. Phosphorus potentially participates in metabolic processes to control MC production, as shown in Figure 7. DOP must first be hydrolyzed by AP secreted by M. aeruginosa, and then Pi can be taken up into the cell through the phosphate transport system [57]. DIP, on the other hand, can be taken up directly by M. aeruginosa through membrane-bound phosphate transport proteins [58]. Pi was then assimilated to support M. aeruginosa metabolism. M. aeruginosa has a sophisticated regulatory system for the assimilation of phosphorus. Thylakoid membranes can mediate fully functional photosynthetic and respiratory electron transport within the photosynthetic lamellae of cyanobacteria [59]. Phosphorus is absorbed to create ATP, which provides energy for MC synthesis and hence regulates enzyme performance and the MC production pathway. Additionally, polyphosphate (polyP), which serves as a nutrient store for cyanobacterial cells, may be formed from ATP by polyphosphate synthase [60]. In summary, phosphorus as an energy source plays a major role in MC synthesis.
In this study, most key genes encoding F-type ATP synthase were significantly downregulated in the DIP and DOP groups (Figure 2A), possibly decreasing ATP levels and influencing intracellular energy metabolism and hence affecting MC synthesis. These results align with the observation that IMC concentrations decreased with incubation time. Overall, these findings imply that phosphorus may influence the metabolism of M. aeruginosa and regulate MC production by affecting adenylate energy charge. Further, the MC contents and energy charge had a direct relationship, and MC synthesis was determined by energy charge. Our results, at the molecular level, align with those of the previous study [61]. Phosphorus is a crucial component in the production of ATP, which is required in order to provide energy for the synthesis of MCs. Compared to the effects of different phosphorus sources on genes involved in MC synthesis, this might be the primary mechanism by which phosphorus governs the production of MCs [17].
Improving our understanding of MCs’ gene expression and regulation will help us identify their potential cellular functions. Enzymes encoded by the mcy gene cluster play a crucial role in the biosynthesis of MCs. The MCs’ synthase gene cluster comprises 10 genes and is transcribed as two polycistronic transcripts (mcyABC and mcyDEFGHIJ) [62]. The relative expressions of the detected genes involved in MC production (mcyABCDEGJ) are shown in Figure 8. In this study, mcyD expression was downregulated by 0.12 Log2FC in the DIP group and 0.20 Log2FC in the DOP group. mcyD is reported to be essential for MC synthesis because it is involved in the creation of the amino acid Adda (2S,3S,8S,9S)-3-amino-9-methoxy-2,6,8-trimethyl-10-phenyldeca-4,6-dienoicacid, which is one of the key components of MCs [63]. Based on our findings, the downregulated expression of mcyD may be the primary cause of the final decline in IMC contents in cyanobacterial cells, supporting the result of Anam et al. [64], who discovered that mcyD gene expression was consistent with MC contents. Increased transcription levels of the microcystin-related gene mcyA in M. aeruginosa could lead to an elevated production of MCs [65]. The expression of mcyA was upregulated by 0.19 Log2FC in the DIP group, consistent with the findings of a previous study [66]. Meanwhile, the mcyA expression was downregulated by 0.29 Log2FC, possibly having led to lower IMC contents in the DOP group compared to in the DIP group from day 6. Furthermore, mcyJ is known to be a tailoring gene, putatively involved in O-methylation [62], and its expression was slightly upregulated by 0.19 and 0.09 Log2FC in the DIP and DOP groups, respectively, indicating that the synthesis activity of MCs was higher in the DIP group. The differential expression of mcyA and mcyJ may have resulted in higher concentrations of both EMC and IMC in the DIP group. Additionally, in both DIP and DOP groups, no appreciable differences were found in the expressions of some genes involved in peptide synthesis, such as mcyBCEG. Although DIP is considered a more bioavailable P source, DOP often accounts for most of the ambient P pool in aquatic ecosystems. Further, an outbreak of cyanobacterial blooms of M. aeruginosa is more likely to occur and be sustained in lakes and reservoirs with DOP being predominant [67]. In the future, managing eutrophication and lowering MC production from cyanobacterial blooms will require consideration of DOP in aquatic environments.

4. Conclusions

DIP and DOP both promoted the growth of M. aeruginosa, and the values of the chl-a content and Fv/Fm in the DIP group were higher than those in the DOP group, indicating that DIP has a more evident influence on the proliferation of cyanobacteria. The transcriptomic results revealed that the variations in most PS I and PS II genes were slightly different. The PS II central protein PsbA was more frequently active in the DIP group; the difference in PsbA expression might affect photosynthesis, which would, in turn, affect cell density. Notably, the AP synthesis-related genes were significantly more abundant in the DOP group than that in the DIP group, indicating that M. aeruginosa could absorb phosphorus by hydrolyzing extracellular DOP by AP at the molecular level. IMC concentrations first increased and then decreased with time. The expressions of most mcy genes showed a little difference between the DOP and DIP group. The expression of mcyD, the key gene encoding MC synthesis in the mcy cluster, was downregulated in both DIP and DOP groups. This indicates that P controls ATP production, which indirectly controls MC synthesis. These results emphasize the complexities of cyanobacterial cell physiology and related gene expression in the presence of different P sources and provide vital insights into the role of P in controlling HAB generation and toxin production.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w15101938/s1, Table S1: NPvsKD5.DESeq; Table S2: NPvsBGPD5.DESeq.

Author Contributions

Z.L.: investigation, data curation, formal analysis, visualization, writing—original draft, writing—review & editing. L.A.: investigation, formal analysis. F.Y.: investigation, data curation. W.S.: investigation, validation. W.D.: investigation. R.D.: funding acquisition, writing—review & editing, conceptualization. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (No. 51678159).

Data Availability Statement

Acknowledgments

The author is extremely appreciative of the teacher’s advice and the assistance of the research group’s pupils.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Paerl, H.W.; Huisman, J. Blooms like it hot. Science 2008, 320, 57–58. [Google Scholar] [CrossRef]
  2. Guo, L. Doing battle with the green monster of Taihu Lake. Science 2007, 317, 1166. [Google Scholar] [CrossRef]
  3. Bullerjahn, G.S.; McKay, R.M.; Davis, T.W.; Baker, D.B.; Boyer, G.L.; D’Anglada, L.V.; Doucette, G.J.; Ho, J.C.; Irwin, E.G.; Kling, C.L. Global solutions to regional problems: Collecting global expertise to address the problem of harmful cyanobacterial blooms. A Lake Erie case study. Harmful Algae 2016, 54, 223–238. [Google Scholar] [CrossRef] [PubMed]
  4. Krishnan, A.; Mou, X. A brief review of the structure, cytotoxicity, synthesis, and biodegradation of microcystins. Water 2021, 13, 2147. [Google Scholar] [CrossRef]
  5. Dittmann, E.; Wiegand, C. Cyanobacterial toxins–occurrence, biosynthesis and impact on human affairs. Mol. Nutr. Food Res. 2006, 50, 7–17. [Google Scholar] [CrossRef]
  6. Zhou, Y.; Zhang, X.; Li, X.; Jia, P.; Dai, R. Evaluation of changes in Microcystis aeruginosa growth and microcystin production by urea via transcriptomic surveys. Sci. Total Environ. 2019, 655, 181–187. [Google Scholar] [CrossRef] [PubMed]
  7. Zhou, Y.; Li, X.; Xia, Q.; Dai, R. Transcriptomic survey on the microcystins production and growth of Microcystis aeruginosa under nitrogen starvation. Sci. Total Environ. 2020, 700, 134501. [Google Scholar] [CrossRef] [PubMed]
  8. Yang, X.; Bi, Y.; Ma, X.; Dong, W.; Wang, X.; Wang, S. Transcriptomic analysis dissects the regulatory strategy of toxic cyanobacterium Microcystis aeruginosa under differential nitrogen forms. J. Hazard. Mater. 2022, 428, 128276. [Google Scholar] [CrossRef] [PubMed]
  9. Zhou, B.; Wang, Z. Effects of nitrogen and phosphorus on Microcystis aeruginosa growth and microcystin production. Green Process. Synth. 2022, 11, 64–70. [Google Scholar] [CrossRef]
  10. Smith, V.H. Low nitrogen to phosphorus ratios favor dominance by blue-green algae in lake phytoplankton. Science 1983, 221, 669–671. [Google Scholar] [CrossRef]
  11. Bradford, M.E.; Peters, R.H. The relationship between chemically analyzed phosphorus fractions and bioavailable phosphorus 1. Limnol. Oceanogr. 1987, 32, 1124–1137. [Google Scholar] [CrossRef]
  12. Yuan, Y.; Bi, Y.; Hu, Z. Phytoplankton communities determine the spatio-temporal heterogeneity of alkaline phosphatase activity: Evidence from a tributary of the Three Gorges Reservoir. Sci. Rep. 2017, 7, 16404. [Google Scholar] [CrossRef]
  13. Paerl, H.W.; Hall, N.S.; Calandrino, E.S. Controlling harmful cyanobacterial blooms in a world experiencing anthropogenic and climatic-induced change. Sci. Total Environ. 2011, 409, 1739–1745. [Google Scholar] [CrossRef] [PubMed]
  14. Wan, L.; Chen, X.; Deng, Q.; Yang, L.; Li, X.; Zhang, J.; Song, C.; Zhou, Y.; Cao, X. Phosphorus strategy in bloom-forming cyanobacteria (Dolichospermum and Microcystis) and its role in their succession. Harmful Algae 2019, 84, 46–55. [Google Scholar] [CrossRef]
  15. Ren, L.; Ren, L.; Wang, P.; Wang, C.; Chen, J.; Hou, J.; Qian, J. Algal growth and utilization of phosphorus studied by combined mono-culture and co-culture experiments. Environ. Pollut. 2017, 220, 274–285. [Google Scholar] [CrossRef] [PubMed]
  16. Qin, C.; Liu, H.; Liu, L.; Smith, S.; Sedlak, D.L.; Gu, A.Z. Bioavailability and characterization of dissolved organic nitrogen and dissolved organic phosphorus in wastewater effluents. Sci. Total Environ. 2015, 511, 47–53. [Google Scholar] [CrossRef] [PubMed]
  17. Zhang, Q.; Chen, Y.; Wang, M.; Zhang, J.; Chen, Q.; Liu, D. Molecular responses to inorganic and organic phosphorus sources in the growth and toxin formation of Microcystis aeruginosa. Water Res. 2021, 196, 117048. [Google Scholar] [CrossRef] [PubMed]
  18. Wang, C.; Yang, Y.; Yang, B.; Lin, H.; Miller, T.R.; Newton, R.J.; Guo, L. Causal relationship between alkaline phosphatase activities and phosphorus dynamics in a eutrophic coastal lagoon in Lake Michigan. Sci. Total Environ. 2021, 787, 147681. [Google Scholar] [CrossRef]
  19. Likens, G.; Bartsch, A.; Lauff, G.; Hobbie, J. Nutrients and eutrophication. Science 1971, 172, 873–874. [Google Scholar] [CrossRef]
  20. Hecky, R.; Kilham, P. Nutrient limitation of phytoplankton in freshwater and marine environments: A review of recent evidence on the effects of enrichment 1. Limnol. Oceanogr. 1988, 33, 796–822. [Google Scholar] [CrossRef]
  21. Li, J.; Wang, Z.; Cao, X.; Wang, Z.; Zheng, Z. Effect of orthophosphate and bioavailability of dissolved organic phosphorous compounds to typically harmful cyanobacterium Microcystis aeruginosa. Mar. Pollut. Bull. 2015, 92, 52–58. [Google Scholar] [CrossRef] [PubMed]
  22. Rippka, R.; Deruelles, J.; Waterbury, J.B.; Herdman, M.; Stanier, R.Y. Generic assignments, strain histories and properties of pure cultures of cyanobacteria. Microbiology 1979, 111, 1–61. [Google Scholar] [CrossRef]
  23. Zhang, F.; Lin, L.; Li, W.; Fang, D.; Lv, Z.; Li, M.; Ma, G.; Wang, Y.; Wang, L.; He, L. Long-Term Study of Monitoring History and Change Trends in Surface Water Quality in China. Water 2022, 14, 2134. [Google Scholar] [CrossRef]
  24. Dai, R.; Liu, H.; Qu, J.; Zhao, X.; Hou, Y. Effects of amino acids on microcystin production of the Microcystis aeruginosa. J. Hazard. Mater. 2009, 161, 730–736. [Google Scholar] [CrossRef]
  25. Lorenzen, C.J. Determination of Chlorophyll and pheopigments: Spectrophotometric equations. Limnol. Ocean. 1967, 12, 343–346. [Google Scholar] [CrossRef]
  26. Sun, Y.; Lin, S.; Wang, D.; Kestel, S.; Lee, G. The coagulation–ultrafiltration process from pilot testing to full-scale design for the achievement of ultra-low phosphorus limit. Water Pract. Technol. 2016, 11, 556–568. [Google Scholar] [CrossRef]
  27. Metcalf, J.; Bell, S.; Codd, G. Production of novel polyclonal antibodies against the cyanobacterial toxin microcystin-LR and their application for the detection and quantification of microcystins and nodularin. Water Res. 2000, 34, 2761–2769. [Google Scholar] [CrossRef]
  28. Ma, Y.; Yan, F.; An, L.; Shen, W.; Tang, T.; Li, Z.; Dai, R. Transcriptome analysis of changes in M. aeruginosa growth and microcystin production under low concentrations of ethinyl estradiol. Sci. Total Environ. 2023, 859, 160226. [Google Scholar] [CrossRef]
  29. Yuan, R.; Li, J.; Li, Y.; Ren, L.; Wang, S.; Kong, F. Formation mechanism of the Microcystis aeruginosa bloom in the water with low dissolved phosphorus. Mar. Pollut. Bull. 2019, 148, 194–201. [Google Scholar] [CrossRef]
  30. Zhang, T.; Lu, X.; Yu, R.; Qin, M.; Wei, C.; Hong, S. Response of extracellular and intracellular alkaline phosphatase in Microcystis aeruginosa to organic phosphorus. Environ. Sci. Pollut. Res. 2020, 27, 42304–42312. [Google Scholar] [CrossRef]
  31. Murchie, E.H.; Lawson, T. Chlorophyll fluorescence analysis: A guide to good practice and understanding some new applications. J. Exp. Bot. 2013, 64, 3983–3998. [Google Scholar] [CrossRef] [PubMed]
  32. Dobrogojski, J.; Adamiec, M.; Luciński, R. The chloroplast genome: A review. Acta Physiol. Plant. 2020, 42, 98. [Google Scholar] [CrossRef]
  33. Li, M.; Hong, J.; Wei, Z.; Yang, J.; Qiao, Z.; Liu, T.; Wang, H.; Bai, F.; Liu, L. Research Progress of Photosynthetic Proteins in Photosynthesis. Hans J. Comput. Biol. 2018, 8, 1–7. [Google Scholar] [CrossRef]
  34. Xiao, Y.; Huang, G.; You, X.; Zhu, Q.; Wang, W.; Kuang, T.; Han, G.; Sui, S.-F.; Shen, J.-R. Structural insights into cyanobacterial photosystem II intermediates associated with Psb28 and Tsl0063. Nat. Plants 2021, 7, 1132–1142. [Google Scholar] [CrossRef]
  35. Gan, F.; Zhang, S.; Rockwell, N.C.; Martin, S.S.; Lagarias, J.C.; Bryant, D.A. Extensive remodeling of a cyanobacterial photosynthetic apparatus in far-red light. Science 2014, 345, 1312–1317. [Google Scholar] [CrossRef] [PubMed]
  36. Debus, R.; Sevrinck, J.; Saetaert, P.; McIntosh, L. Nucleotide sequence of a second psbA gene from the unicellular cyanobacterium Synechocystis 6803. Nucleic Acids Res. 1989, 17, 3991. [Google Scholar]
  37. Yang, Y.; Shi, J.; Jia, Y.; Bai, F.; Yang, S.; Mi, W.; He, S.; Wu, Z. Unveiling the impact of glycerol phosphate (DOP) in the dinoflagellate Peridinium bipes by physiological and transcriptomic analysis. Environ. Sci. Eur. 2020, 32, 38. [Google Scholar] [CrossRef]
  38. Nugent, J.H. Oxygenic photosynthesis: Electron transfer in photosystem I and photosystem II. Eur. J. Biochem. 1996, 237, 519–531. [Google Scholar] [CrossRef]
  39. Yamori, W.; Shikanai, T. Physiological functions of cyclic electron transport around photosystem I in sustaining photosynthesis and plant growth. Annu. Rev. Plant Biol. 2016, 67, 81–106. [Google Scholar] [CrossRef]
  40. Dong, C.; Zhang, H.; Yang, Y.; He, X.; Liu, L.; Fu, J.; Shi, J.; Wu, Z. Physiological and transcriptomic analyses to determine the responses to phosphorus utilization in Nostocs. Harmful Algae 2019, 84, 10–18. [Google Scholar] [CrossRef]
  41. Zhang, Q.; Song, Q.; Wang, C.; Zhou, C.; Lu, C.; Zhao, M. Effects of glufosinate on the growth of and microcystin production by Microcystis aeruginosa at environmentally relevant concentrations. Sci. Total Environ. 2017, 575, 513–518. [Google Scholar] [CrossRef]
  42. Kathuria, S.; Martiny, A.C. Prevalence of a calcium-based alkaline phosphatase associated with the marine cyanobacterium Prochlorococcus and other ocean bacteria. Environ. Microbiol. 2011, 13, 74–83. [Google Scholar] [CrossRef] [PubMed]
  43. Ma, J.; Wang, P.; Ren, L.; Wang, X.; Paerl, H.W. Using alkaline phosphatase activity as a supplemental index to optimize predicting algal blooms in phosphorus-deficient lakes: A case study of Lake Taihu, China. Ecol. Indic. 2019, 103, 698–712. [Google Scholar] [CrossRef]
  44. Xie, E.; Su, Y.; Deng, S.; Kontopyrgou, M.; Zhang, D. Significant influence of phosphorus resources on the growth and alkaline phosphatase activities of Microcystis aeruginosa. Environ. Pollut. 2021, 268, 115807. [Google Scholar] [CrossRef] [PubMed]
  45. Lamarche, M.G.; Wanner, B.L.; Crépin, S.; Harel, J. The phosphate regulon and bacterial virulence: A regulatory network connecting phosphate homeostasis and pathogenesis. FEMS Microbiol. Rev. 2008, 32, 461–473. [Google Scholar] [CrossRef]
  46. Shi, J.-Q.; Ou-Yang, T.; Yang, S.-Q.; Zhao, L.; Ji, L.-L.; Wu, Z.-X. Transcriptomic responses to phosphorus in an invasive cyanobacterium, Raphidiopsis raciborskii: Implications for nutrient management. Harmful Algae 2022, 111, 102150. [Google Scholar] [CrossRef] [PubMed]
  47. Harke, M.J.; Berry, D.L.; Ammerman, J.W.; Gobler, C.J. Molecular response of the bloom-forming cyanobacterium, Microcystis aeruginosa, to phosphorus limitation. Microb. Ecol. 2012, 63, 188–198. [Google Scholar] [CrossRef]
  48. Wanner, B. Gene regulation by phosphate in enteric bacteria. J. Cell. Biochem. 1993, 51, 47–54. [Google Scholar] [CrossRef]
  49. Torriani, A. From cell membrane to nucleotides: The phosphate regulon in Escherichia coli. Bioessays 1990, 12, 371–376. [Google Scholar] [CrossRef]
  50. Lubin, E.A.; Henry, J.T.; Fiebig, A.; Crosson, S.; Laub, M.T. Identification of the PhoB regulon and role of PhoU in the phosphate starvation response of Caulobacter crescentus. J. Bacteriol. 2016, 198, 187–200. [Google Scholar] [CrossRef]
  51. Dyhrman, S.T.; Jenkins, B.D.; Rynearson, T.A.; Saito, M.A.; Mercier, M.L.; Alexander, H.; Whitney, L.P.; Drzewianowski, A.; Bulygin, V.V.; Bertrand, E.M. The transcriptome and proteome of the diatom Thalassiosira pseudonana reveal a diverse phosphorus stress response. PLoS ONE 2012, 7, e33768. [Google Scholar] [CrossRef] [PubMed]
  52. Pietsch, J.; Bornmann, K.; Schmidt, W. Relevance of intra-and extracellular cyanotoxins for drinking water treatment. Acta Hydrochim. Hydrobiol. 2002, 30, 7–15. [Google Scholar] [CrossRef]
  53. Tilahun, S.; Kifle, D.; Zewde, T.W.; Johansen, J.A.; Demissie, T.B.; Hansen, J.H. Temporal dynamics of intra-and extra-cellular microcystins concentrations in Koka reservoir (Ethiopia): Implications for public health risk. Toxicon 2019, 168, 83–92. [Google Scholar] [CrossRef] [PubMed]
  54. Wei, N.; Song, L.; Gan, N. Quantitative proteomic and microcystin production response of Microcystis aeruginosa to phosphorus depletion. Microorganisms 2021, 9, 1183. [Google Scholar] [CrossRef] [PubMed]
  55. Peng, G.; Lin, S.; Fan, Z.; Wang, X. Transcriptional and physiological responses to nutrient loading on toxin formation and photosynthesis in Microcystis aeruginosa FACHB-905. Toxins 2017, 9, 168. [Google Scholar] [CrossRef]
  56. Rowan, K. Phosphorus metabolism in plants. Int. Rev. Cytol. 1966, 19, 301–390. [Google Scholar]
  57. Bai, F.; Shi, J.; Yang, S.; Yang, Y.; Wu, Z. Interspecific competition between Cylindrospermopsis raciborskii and Microcystis aeruginosa on different phosphorus substrates. Environ. Sci. Pollut. Res. 2020, 27, 42264–42275. [Google Scholar] [CrossRef]
  58. Qiu, H.; Geng, J.; Ren, H.; Xia, X.; Wang, X.; Yu, Y. Physiological and biochemical responses of Microcystis aeruginosa to glyphosate and its Roundup® formulation. J. Hazard. Mater. 2013, 248, 172–176. [Google Scholar] [CrossRef]
  59. Singh, V.P.; Passarini, M.R.Z.; Kumar, A. Cyanobacterial biology in twenty-first century. Front. Microbiol. 2023, 14, 1184669. [Google Scholar] [CrossRef]
  60. Sanz-Luque, E.; Bhaya, D.; Grossman, A.R. Polyphosphate: A multifunctional metabolite in cyanobacteria and algae. Front. Plant Sci. 2020, 11, 938. [Google Scholar] [CrossRef]
  61. Dai, R.; Liu, H.; Qu, J.; Zhao, X.; Ru, J.; Hou, Y. Relationship of energy charge and toxin content of Microcystis aeruginosa in nitrogen-limited or phosphorous-limited cultures. Toxicon 2008, 51, 649–658. [Google Scholar] [CrossRef]
  62. Tillett, D.; Dittmann, E.; Erhard, M.; Von Döhren, H.; Börner, T.; Neilan, B.A. Structural organization of microcystin biosynthesis in Microcystis aeruginosa PCC7806: An integrated peptide–polyketide synthetase system. Chem. Biol. 2000, 7, 753–764. [Google Scholar] [CrossRef] [PubMed]
  63. Ko, S.-R.; Lee, Y.-K.; Srivastava, A.; Park, S.-H.; Ahn, C.-Y.; Oh, H.-M. The selective inhibitory activity of a fusaricidin derivative on a bloom-forming cyanobacterium Microcystis sp. J. Microbiol. Biotechnol. 2019, 29, 59–65. [Google Scholar] [CrossRef] [PubMed]
  64. Anam, G.B.; Guda, D.R.; Ahn, Y.-H. Hormones induce the metabolic growth and cytotoxin production of Microcystis aeruginosa under terpinolene stress. Sci. Total Environ. 2021, 769, 145083. [Google Scholar] [CrossRef]
  65. Zhou, J.; Gao, L.; Lin, Y.; Pan, B.; Li, M. Micrometer scale polystyrene plastics of varying concentrations and particle sizes inhibit growth and upregulate microcystin-related gene expression in Microcystis aeruginosa. J. Hazard. Mater. 2021, 420, 126591. [Google Scholar] [CrossRef]
  66. Srivastava, A.; Ko, S.-R.; Ahn, C.-Y.; Oh, H.-M.; Ravi, A.K.; Asthana, R.K. Microcystin biosynthesis and mcyA expression in geographically distinct Microcystis strains under different nitrogen, phosphorus, and boron regimes. BioMed Res. Int. 2016, 2016, 5985987. [Google Scholar] [CrossRef]
  67. Zhang, X.; Wang, Z.; Luo, Z.; Chen, Y.; Huang, X. Insights into the conversion of dissolved organic phosphorus favors algal bloom, arsenate biotransformation and microcystins release of Microcystis aeruginosa. J. Environ. Sci. 2023, 125, 205–214. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Changes in cell density (A), specific growth rate (B), and chlorophyll-a (C) and Fv/Fm (D) of the M. aeruginosa with time under DIP and DOP conditions. The error bars represent the standard errors.
Figure 1. Changes in cell density (A), specific growth rate (B), and chlorophyll-a (C) and Fv/Fm (D) of the M. aeruginosa with time under DIP and DOP conditions. The error bars represent the standard errors.
Water 15 01938 g001
Figure 2. Circle heat map of gene expression of genes involved in photosynthesis (A) and photosynthesis on the thylakoid membrane of M. aeruginosa with key genes involved in PS I and PS II (B). PQ: plastoquinone; PQH2: plastid hydroquinone; PC: plastocyanin. A reduction in transcript abundance corresponds to a blue hue, whereas an increase in transcript abundance corresponds to a red color. White indicates that there is no difference between the experimental and control conditions.
Figure 2. Circle heat map of gene expression of genes involved in photosynthesis (A) and photosynthesis on the thylakoid membrane of M. aeruginosa with key genes involved in PS I and PS II (B). PQ: plastoquinone; PQH2: plastid hydroquinone; PC: plastocyanin. A reduction in transcript abundance corresponds to a blue hue, whereas an increase in transcript abundance corresponds to a red color. White indicates that there is no difference between the experimental and control conditions.
Water 15 01938 g002
Figure 3. Changes in phosphorus concentration under different phosphorus source media with time. (A) TP concentrations of both groups; (B) Pi concentrations of both groups.
Figure 3. Changes in phosphorus concentration under different phosphorus source media with time. (A) TP concentrations of both groups; (B) Pi concentrations of both groups.
Water 15 01938 g003
Figure 4. Variations of Alkaline phosphatase concentration.
Figure 4. Variations of Alkaline phosphatase concentration.
Water 15 01938 g004
Figure 5. The heatmap of genes involved in P transportation and AP activity (A) and the schematic diagram of phosphorus transportation under DIP and DOP and metabolic pathways potentially involved in the utilization of different P sources (B). DIP and DOP are represented by circles and squares, respectively. Reductions in transcript abundance correspond to blue hues, whereas increases in transcript abundance correspond to red colors. White indicates that there is no difference between the experimental and control conditions. PPP: polyphosphatase.
Figure 5. The heatmap of genes involved in P transportation and AP activity (A) and the schematic diagram of phosphorus transportation under DIP and DOP and metabolic pathways potentially involved in the utilization of different P sources (B). DIP and DOP are represented by circles and squares, respectively. Reductions in transcript abundance correspond to blue hues, whereas increases in transcript abundance correspond to red colors. White indicates that there is no difference between the experimental and control conditions. PPP: polyphosphatase.
Water 15 01938 g005aWater 15 01938 g005b
Figure 6. Variations in MCs under a different form of P groups. (A) IMC (µg/L) and (B) EMC (µg/L). The error bars represent the standard errors.
Figure 6. Variations in MCs under a different form of P groups. (A) IMC (µg/L) and (B) EMC (µg/L). The error bars represent the standard errors.
Water 15 01938 g006
Figure 7. Phosphorus potentially participates in metabolic processes that control the production of MCs.
Figure 7. Phosphorus potentially participates in metabolic processes that control the production of MCs.
Water 15 01938 g007
Figure 8. Variations of the mcy gene expressions in the clusters of the DIP and DOP groups.
Figure 8. Variations of the mcy gene expressions in the clusters of the DIP and DOP groups.
Water 15 01938 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, Z.; An, L.; Yan, F.; Shen, W.; Du, W.; Dai, R. Evaluation of the Effects of Different Phosphorus Sources on Microcystis aeruginosa Growth and Microcystin Production via Transcriptomic Surveys. Water 2023, 15, 1938. https://doi.org/10.3390/w15101938

AMA Style

Li Z, An L, Yan F, Shen W, Du W, Dai R. Evaluation of the Effects of Different Phosphorus Sources on Microcystis aeruginosa Growth and Microcystin Production via Transcriptomic Surveys. Water. 2023; 15(10):1938. https://doi.org/10.3390/w15101938

Chicago/Turabian Style

Li, Zihao, Lili An, Feng Yan, Wendi Shen, Wenjun Du, and Ruihua Dai. 2023. "Evaluation of the Effects of Different Phosphorus Sources on Microcystis aeruginosa Growth and Microcystin Production via Transcriptomic Surveys" Water 15, no. 10: 1938. https://doi.org/10.3390/w15101938

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