Transcriptional Analysis of Microcystis aeruginosa Co-Cultured with Algicidal Bacteria Brevibacillus laterosporus

Harmful algal blooms caused huge ecological damage and economic losses around the world. Controlling algal blooms by algicidal bacteria is expected to be an effective biological control method. The current study investigated the molecular mechanism of harmful cyanobacteria disrupted by algicidal bacteria. Microcystis aeruginosa was co-cultured with Brevibacillus laterosporus Bl-zj, and RNA-seq based transcriptomic analysis was performed compared to M. aeruginosa, which was cultivated separately. A total of 1706 differentially expressed genes were identified, which were mainly involved in carbohydrate metabolism, energy metabolism and amino acid metabolism. In the co-cultured group, the expression of genes mainly enriched in photosynthesis and oxidative phosphorylation were significantly inhibited. However, the expression of the genes related to fatty acid synthesis increased. In addition, the expression of the antioxidant enzymes, such as 2-Cys peroxiredoxin, was increased. These results suggested that B. laterosporus could block the electron transport by attacking the PSI system and complex I of M. aeruginosa, affecting the energy acquisition and causing oxidative damage. This further led to the lipid peroxidation of the microalgal cell membrane, resulting in algal death. The transcriptional analysis of algicidal bacteria in the interaction process can be combined to explain the algicidal mechanism in the future.


Introduction
With the increase in water eutrophication, harmful algal blooms (HABs) have frequently occurred worldwide. In HABs, cyanobacteria are the most common group of microalgae, which can produce toxins [1][2][3]. Cyanobacterial blooms can cause serious water environmental problems, such as water oxygen deprivation, odorous substances and toxins, which are harmful to drinking water, aquatic animals, even threatening to human health [4][5][6]. Various methods have been proposed for removing or inhibiting cyanobacterial blooms, which include physical, chemical and biological methods. However, physical and chemical methods are expensive and cause secondary pollution, which is potentially harmful to aquatic ecosystems [7]. In contrast, biological methods are effective strategies to control HABs [8].
Algicidal bacteria are a group of bacteria that can kill or inhibit microalgae. At present, many algicidal bacteria, such as Bacillus [9], Chryseobacterium [10], Sulfitobacter [11], Hahella [12] and Stenotrophomonas [13], have been previously reported. There are two main mechanisms of algicidal bacteria: (1) direct attack: algicidal bacteria directly invade the microalgal cells; (2) indirect attack: algicidal bacteria inhibit microalgal growth by secreting extracellular active substances to lyse them or by competing with them for nutrients. The 2 of 13 first study of algicidal bacteria was reported in Myxobacter that inhibited the growth of blue-green algae by secreting some cellulolytic enzymes and digesting the cell wall of the host [14]. Similarly, Furusawa et al. [15] found that Saprospira sp. SS98-5 exhibited an algicidal effect on Chaetoceros ceratosporum, using its filamentous cells to trap the diatom cells and dissolve their cytoderm, which were in contact with the bacteria. The specific or non-specific extracellular substances secreted by algicidal bacteria during their metabolism can destroy the structure of algal cells and lyse them. Extracellular active substances often act through various mechanisms to dissolve the cyanobacteria, including blocking electron transfer, reducing algal photosynthetic system activity, causing oxidative damage, inhibiting intracellular protein and carbohydrate synthesis and disrupting normal cell metabolism and cellular structure [16][17][18]. Yang et al. found that prodigiosin (PG) from Hahella sp. KA22 inhibited Microcystis aeruginosa by generating reactive oxygen species (ROS), which induce lipid peroxidation, disrupt the membrane system and obliterate the function of the photosystem [12]. In addition, many proteins of M. aeruginosa involved in important metabolic processes in response to PG stress were differentially expressed, which may trigger necrotic-like or apoptotic-like cell death with features similar to those in eukaryotes [12].
The algicidal mechanism has been extensively studied by measuring the physiological and biochemical characteristics of microalgae. A few studies have been conducted on the mechanism of algal-lysing at the molecular and protein levels. RNAseq and microarray analysis were determined to study the interactions between Roseobacters and dinoflagellates. The data suggested that at the early mutualistic phase of the symbiosis, polyhydroxyalkanoate (PHA) degradation might be the main carbon and energy source of Dinoroseobacter shibae, supplemented in the light by degradation of dimethylsulfoniopropionate (DMSP) and aerobic anoxygenic photosynthesis [19]. Proteomic analysis was applied to investigate the algicidal process of Trametes versicolor F21a on M. aeruginosa, and 30 fungal enzymes with endo-or exoglycosidase activities were significantly up-regulated, suggesting that these enzymes may degrade lipopolysaccharides, peptidoglycans and alginic acid of algal cells [20]. All these studies were analyzed from the algicidal bacteria or fungi, while the analysis of lysed microalgae at a molecular level may provide a new understanding of the algicidal mechanism.
Brevibacillus laterosporus is a Gram-positive bacterium that can produce antibacterial molecules and other secondary metabolites that have a wide range of biological activities [21]. In a previous study, the B. laterosporus strain Bl-zj was isolated from the intertidal soil of a park in Zhanjiang, which has been proved to have an algae-lytic effect on cyanobacteria, and the potential algicidal factors were screened by genome analysis [22,23].
The research on the mechanism of algicidal bacteria has been mainly conducted using various physiological and biochemical analyses; however, there is a scarcity of reports at the gene level, especially from the lysed microalgae. In this study, the algicidal bacteria B. laterosporus Bl-zj and M. aeruginosa were co-cultured, and transcriptome analyses were performed to investigate the response of M. aeruginosa to algicidal bacteria.

Cyanobacteria and Algicidal Bacterium
M. aeruginosa FACHB 905 was purchased from the Freshwater Algae Culture Collection at the Institute of Hydrobiology, Wuhan, China. For inoculum preparation, it was cultured for 7 d to reach the exponential growth phase. It was cultivated in BG11 [24] medium at 28 ± 1 • C under 50 µmol·m −2 ·s −1 with a photoperiod of 12:12 h light: dark cycle. The concentration of M. aeruginosa was adjusted to 1 × 10 7 cells/mL for further experiments.
B. laterosporus Bl-zj was isolated and preserved by the Laboratory of Algae Resource Development and Culture Environment Ecological Restoration of Guangdong Ocean University. It was cultured using Beef extract-peptone (BP) medium and incubated for 24 h at 28 • C and 150 rpm to reach the logarithmic growth phase. The concentration of B. laterosporus was adjusted to 1 × 10 7 cells/mL for further experiments.

Measurement of the Algicidal Efficiency
Equivalent volumes of microalga and bacterium (10 mL) were mixed, and total culture volume was added to 100 mL with BG11 medium (MB group). The control group was just 10 mL of M. aeruginosa added in BG11 medium to the total volume of 100 mL (BG11 group). Each group contained triplicate samples. The experiment lasted 4 days under the same condition as 2.1 when the MB group culture was turned yellow, and the microalgal biomass was at a low proportion.
Spectrophotometer with a 1 cm light path cuvette was used to measure and calculate chlorophyll a concentration, according to the methods of Jeffrey and Humphrey [25]. In brief, the cultures of each group were harvested at 0, 1, 2, 3, 4 d, and chlorophyll a was extracted with 90% (v/v) acetone for 24 h in the dark. The absorbance of the supernatant at 630 nm, 647 nm and 664 nm was determined. The chlorophyll a concentration was calculated according to the formula: Chl a (mg/L) = (11.85 × A 664 − 1.54 × A 647 − 0.08 × A 630 ) × V 1 /V 2 , where A 630 , A 647 and A 664 represent absorbance, V 1 and V 2 are the volume of 90% (v/v) acetone (mL) and cultures (mL), respectively. The removal rate of microalgal cells (R, %) was represented as the algicidal effect and calculated by the formula: R = (C 0 − C 1 )/C 0 × 100%, where C 0 and C 1 represent the chlorophyll a concentration of BG11 group and MB group, respectively.

Transcriptomic Samples Preparation
The M. aeruginosa and B. laterosporus were co-cultured at a volume ratio of 1:1 (50 mL M. aeruginosa, 50 mL B. laterosporus and 400 mL BG11 medium). M. aeruginosa was separately cultured as control (50 mL M. aeruginosa and 450 mL BG11 medium). Samples from the control group (named "CMA") and experimental group were collected on the second and fourth day (named "MB2" and "MB4", respectively), as the algicidal effect was obviously on the fourth day. The cocultures were centrifuged at 5000 rpm, 4 • C for 10 min to collect the pellets and then quickly frozen in liquid nitrogen and stored at −80 • C. Each group contained triplicate samples.
The differences in gene expression were analyzed by HTSeq (version0.6.1p2) for the screening of differentially expressed genes (DEGs) and principal components analysis. The read count mapped to each gene was calculated as the original expression of the gene and the FPKM (Fragments Per Kilo bases per Million fragments) was used to normalize the expression levels. The DEGs screening conditions were as follows: log 2 |fold change| > 1, p-value < 0.05. DEGs were further annotated to Gene Ontology (GO) database and Kyoto Encyclopedia of Genes and Genomes (KEGG) database.

Quantitative Real-Time PCR Validation
Quantitative PCR (qPCR) was performed to verify the results of mRNA-seq analysis. The cDNA template was synthesized by reverse transcription using a kit (HiScript R III RT SuperMix for qPCR (+gDNA wiper), Vazyme Biotechnology Co., LTD., Nanjing, China). The qPCR was performed using SYBR green (AceQ Universal SYBR qPCR Master Mix, Vazyme Biotechnology Co., LTD., Nanjing, China) on a 7500 Fast Real-time PCR System (Applied Biosystems, Foster, CA, USA). The primers were designed according to the gene sequences from M. aeruginosa PCC 7806SL complete genome, with 16S rRNA [26] as the reference gene and synthesized by GENEWIZ (Suzhou, China) (the primers were shown in Supplementary Table S1).

The Algicidal Efficiency
Compared with the control group (BG11), the chlorophyll a concentration of M. aeruginosa significantly decreased with the addition of B. laterosporus from the second day ( Figure 1). The removal rates of M. aeruginosa in the MB group were 34.39%, 72.36% and 92.30% on the second, third and fourth days, respectively, showing an increasing algicidal effect over time.
Vazyme Biotechnology Co., LTD., Nanjing, China) on a 7500 Fast Real-tim (Applied Biosystems, Foster, CA, USA). The primers were designed accordi sequences from M. aeruginosa PCC 7806SL complete genome, with 16S rRN reference gene and synthesized by GENEWIZ (Suzhou, China) (the primer in Supplementary Table S1).

The Algicidal Efficiency
Compared with the control group (BG11), the chlorophyll a concentrat ginosa significantly decreased with the addition of B. laterosporus from th ( Figure 1). The removal rates of M. aeruginosa in the MB group were 34.39% 92.30% on the second, third and fourth days, respectively, showing an incre effect over time.

Illumina Sequencing Assembly Data Quality Analysis
Nine cDNA libraries from three groups were sequenced to study the t of M. aeruginosa under the effect of B. laterosporus (Table 1). An average of bp raw data and 64,858,225 raw reads were obtained from each sample. Aft adaptors and low-quality sequences from raw data, an average of 7,800,58 data and 52,003,890 clean reads were screened. After mapping to the refer an average of 37,785,780 mapped reads were obtained. All Q20 and Q30 valu sequences in the samples exceeded 96.94% and 92.9%, respectively.

Illumina Sequencing Assembly Data Quality Analysis
Nine cDNA libraries from three groups were sequenced to study the transcriptomes of M. aeruginosa under the effect of B. laterosporus (Table 1). An average of 9,728,733,733 bp raw data and 64,858,225 raw reads were obtained from each sample. After filtering the adaptors and low-quality sequences from raw data, an average of 7,800,583,533 bp clean data and 52,003,890 clean reads were screened. After mapping to the reference genome, an average of 37,785,780 mapped reads were obtained. All Q20 and Q30 values of the read sequences in the samples exceeded 96.94% and 92.9%, respectively. Principal component analysis (PCA) was performed for each sample according to the expression to observe similarities between sample groups ( Figure 2). The three replicates of each group (CMA, MB2 and MB4) were gathered with close distance, but different treatment groups separated from others at a long distance, indicating good repeatability data. Principal component analysis (PCA) was performed for each sample accordi expression to observe similarities between sample groups ( Figure 2). The three r of each group (CMA, MB2 and MB4) were gathered with close distance, but treatment groups separated from others at a long distance, indicating good repe data.

Identification of DEGs
In this study, transcriptomic analysis was used to investigate the impact of bacteria B. laterosporus Bl-zj on M. aeruginosa metabolism at the molecular level showed that a total of 1324 DEGs were found in the MB2 group compared to which 612 were up-regulated and 712 were down-regulated ( Figure 3). There w DEGs in the MB4 group compared to CMA, of which 450 were up-regulated and down-regulated. A total of 1706 DEGs were identified; among them, there were regulated and 528 down-regulated DEGs both in MB2 vs. CMA comparison and CMA comparison.

Identification of DEGs
In this study, transcriptomic analysis was used to investigate the impact of algicidal bacteria B. laterosporus Bl-zj on M. aeruginosa metabolism at the molecular level. Results showed that a total of 1324 DEGs were found in the MB2 group compared to CMA, of which 612 were up-regulated and 712 were down-regulated ( Figure 3). There were 1289 DEGs in the MB4 group compared to CMA, of which 450 were up-regulated and 839 were down-regulated. A total of 1706 DEGs were identified; among them, there were 377 up-regulated and 528 down-regulated DEGs both in MB2 vs. CMA comparison and MB4 vs. CMA comparison. Principal component analysis (PCA) was performed for each sample accordin expression to observe similarities between sample groups ( Figure 2). The three re of each group (CMA, MB2 and MB4) were gathered with close distance, but d treatment groups separated from others at a long distance, indicating good repea data.

Identification of DEGs
In this study, transcriptomic analysis was used to investigate the impact of a bacteria B. laterosporus Bl-zj on M. aeruginosa metabolism at the molecular level. showed that a total of 1324 DEGs were found in the MB2 group compared to C which 612 were up-regulated and 712 were down-regulated ( Figure 3). There we DEGs in the MB4 group compared to CMA, of which 450 were up-regulated and 8 down-regulated. A total of 1706 DEGs were identified; among them, there were regulated and 528 down-regulated DEGs both in MB2 vs. CMA comparison and CMA comparison.

Functional Classification of the DEGs by GO and KEGG Pathway Analysis
The GO and KEGG pathways were used to analyze the DEGs and determine their main biological functions to explore the molecular changes in M. aeruginosa under the effect of B. laterosporus Bl-zj.
The top 10 most significant enriched genes in three GO categories were summarized (Figure 4). In the cellular component category, the GO terms mainly enriched in the MB2 vs. CMA comparison included "photosystem I", "cell part" and "cytoplasm". The "intracellular part", "intracellular" and "cell part" were mainly enriched in the MB4 vs. CMA comparison. Similarly, in the molecular function category, "ion binding", "quinone binding" and "small molecule binding" were significantly enriched in the MB2 vs. CMA comparison. The "ion binding", "small molecule binding" and "oxidoreductase activity" were significantly enriched in the MB4 vs. CMA comparison. Further, in the biological process, "prophyrin-containing compound biosynthetic process", "oxidative phosphorylation" and "prophyrin-containing compound metabolic process" were significantly enriched in the MB2 vs. CMA comparison. The "photosynthesis", "pigment biosynthetic process" and "oxidation-reduction process" were significantly enriched in the MB4 vs. CMA comparison.
The biological functions associated with DEGs were further analyzed using the KEGG database. The up-regulated and down-regulated genes were involved in carbohydrate metabolism, amino acid metabolism, metabolism of cofactors and vitamins and energy metabolism in MB2 vs. CMA comparison and MB4 vs. CMA comparison ( Table 2). In addition, the genes of lipid metabolism and translation were mainly up-regulated and the genes of energy metabolism, replication and repair and membrane transport were mainly down-regulated. The GO and KEGG pathways were used to analyze the DEGs and determine their main biological functions to explore the molecular changes in M. aeruginosa under the effect of B. laterosporus Bl-zj.
The top 10 most significant enriched genes in three GO categories were summarized (Figure 4). In the cellular component category, the GO terms mainly enriched in the MB2 vs. CMA comparison included "photosystem I", "cell part" and "cytoplasm". The "intracellular part", "intracellular" and "cell part" were mainly enriched in the MB4 vs. CMA comparison. Similarly, in the molecular function category, "ion binding", "quinone binding" and "small molecule binding" were significantly enriched in the MB2 vs. CMA comparison. The "ion binding", "small molecule binding" and "oxidoreductase activity" were significantly enriched in the MB4 vs. CMA comparison. Further, in the biological process, "prophyrin-containing compound biosynthetic process", "oxidative phosphorylation" and "prophyrin-containing compound metabolic process" were significantly enriched in the MB2 vs. CMA comparison. The "photosynthesis", "pigment biosynthetic process" and "oxidation-reduction process" were significantly enriched in the MB4 vs. CMA comparison.   The biological functions associated with DEGs were further analyzed using the KEGG database. The up-regulated and down-regulated genes were involved in carbohydrate metabolism, amino acid metabolism, metabolism of cofactors and vitamins and energy metabolism in MB2 vs. CMA comparison and MB4 vs. CMA comparison ( Table 2). In addition, the genes of lipid metabolism and translation were mainly up-regulated and the genes of energy metabolism, replication and repair and membrane transport were mainly down-regulated.

Algicidal-Related Gene and Pathway Analysis
The expression of genes related to photosynthesis, oxidative phosphorylation and fatty acid synthesis in M. aeruginosa were significantly changed under the effect of B. laterosporus Bl-zj (Table 3). There were nine photosynthesis-related genes that were down-regulated on the second and fourth days, compared to the CMA group, among which seven genes were related to photosystem I (PSI). Nine genes related to oxidative phosphorylation were downregulated both on the second and fourth days, compared to the control. The expression of ND5 and Ndufs8 were most significantly decreased with log 2 (fold changes) of 3.38 and 3.43 on the fourth day compared to the control, respectively. In addition, the expression of two genes (FabF and FabG1) involved in fatty acid synthesis increased on the second day compared to the control, while YOXD was increased on the fourth day compared to the control. Similarly, three genes (FabH, FabG2 and FabZ) were up-regulated on the second and fourth days compared to the control.

Quantitative Real-Time PCR Validation
To confirm the RNA-seq results, seven genes were selected for further qRT-PCR analysis. Figure 5a,b shows the expression levels of each gene on the second and fourth days, respectively. Although the fold changes were different, the expression trends of these seven genes were consistent in both RNA-seq and qRT-PCR results.

Discussion
In the present study, second-and fourth-day transcriptomes of M. aeruginosa in coculture with B. laterosporus Bl-zj were obtained and compared using next-generation sequencing technology. We found that several pathways significantly changed, which were probably the critical factors as B. laterosporus attacked the cells of M. aeruginosa.
Photosynthesis is a critical cellular process that determines the growth of microalgae. In this process, light energy is captured and used to synthesize carbohydrates and produce oxygen while consuming carbon dioxide. Transcriptomic analysis showed that PSI was seriously affected by B. laterosporus Bl-zj. The principal subunits of the reaction center (PsaA and PsaB) and other subunits (PsaC, PsaE, PsaD, PsaF and PsaL) were significantly down-regulated on the second day and continuously inhibited on the fourth day (Table  3). PSI is a supercomplex of a reaction center and light-harvesting complexes [27]. PsaA and PsaB are the central core proteins that bonded all cofactors of the electron transport chain except the cluster of ferredoxins and harbor the electron transport chain [28]. PsaC can bind to the terminal electron acceptor and participate in the transfer of electrons to ferredoxin [29]. PsaE was responsible for the dissociation of ferredoxin from the PSI complex and involved in the cyclic electron flow that occurred around the PSI complex [30]. PsaD has been regarded as a key subunit in the assembly, stability and functionality of PSI [31]. PsaF has been found to be related to efficient electron transfer from both plastocyanin and cytochrome c6 to PSI [32]. It has been reported that PsaL was critical for energy transfer from phosphorylated LHCII to the PSI reaction center [33]. In this study, it was observed that Psb28 significantly declined after 2 days compared to the control. Psb28 was the photosystem II reaction center protein, which plays an important role in PSII repair [34]. In addition, plastocyanin and ferredoxin significantly decreased after the second and fourth days compared to the control. Plastocyanin and ferredoxin are small soluble copper protein and iron-sulfur protein, respectively, which are important PSI turnover products induced by light energy involved in electron transfer from plastocyanin to ferredoxin [28,35]. These results showed that under the action of B. laterosporus Bl-zj photosynthesis, especially the PSI of M. aeruginosa, is significantly inhibited. The reduction of photosynthesis-related gene transcripts impeded electron transport, resulting in the loss of equivalent products necessary for the process of carbon assimilation.
Oxidative phosphorylation is the main energy source of aerobic cells and the main pathway of ATP production [36]. It is an enzymatic process in which ATP is formed when

Discussion
In the present study, second-and fourth-day transcriptomes of M. aeruginosa in coculture with B. laterosporus Bl-zj were obtained and compared using next-generation sequencing technology. We found that several pathways significantly changed, which were probably the critical factors as B. laterosporus attacked the cells of M. aeruginosa.
Photosynthesis is a critical cellular process that determines the growth of microalgae. In this process, light energy is captured and used to synthesize carbohydrates and produce oxygen while consuming carbon dioxide. Transcriptomic analysis showed that PSI was seriously affected by B. laterosporus Bl-zj. The principal subunits of the reaction center (PsaA and PsaB) and other subunits (PsaC, PsaE, PsaD, PsaF and PsaL) were significantly downregulated on the second day and continuously inhibited on the fourth day (Table 3). PSI is a supercomplex of a reaction center and light-harvesting complexes [27]. PsaA and PsaB are the central core proteins that bonded all cofactors of the electron transport chain except the cluster of ferredoxins and harbor the electron transport chain [28]. PsaC can bind to the terminal electron acceptor and participate in the transfer of electrons to ferredoxin [29]. PsaE was responsible for the dissociation of ferredoxin from the PSI complex and involved in the cyclic electron flow that occurred around the PSI complex [30]. PsaD has been regarded as a key subunit in the assembly, stability and functionality of PSI [31]. PsaF has been found to be related to efficient electron transfer from both plastocyanin and cytochrome c6 to PSI [32]. It has been reported that PsaL was critical for energy transfer from phosphorylated LHCII to the PSI reaction center [33]. In this study, it was observed that Psb28 significantly declined after 2 days compared to the control. Psb28 was the photosystem II reaction center protein, which plays an important role in PSII repair [34]. In addition, plastocyanin and ferredoxin significantly decreased after the second and fourth days compared to the control. Plastocyanin and ferredoxin are small soluble copper protein and iron-sulfur protein, respectively, which are important PSI turnover products induced by light energy involved in electron transfer from plastocyanin to ferredoxin [28,35]. These results showed that under the action of B. laterosporus Bl-zj photosynthesis, especially the PSI of M. aeruginosa, is significantly inhibited. The reduction of photosynthesis-related gene transcripts impeded electron transport, resulting in the loss of equivalent products necessary for the process of carbon assimilation.
Oxidative phosphorylation is the main energy source of aerobic cells and the main pathway of ATP production [36]. It is an enzymatic process in which ATP is formed when electrons are transferred from substrate to oxygen and ADP is phosphorylated. Transcriptomic analysis showed that the genes related to oxidative phosphorylation significantly decreased after the second and fourth days, relative to the control (Table 3). NADH: ubiquinone oxidoreductase (Complex I) is the first and most complex enzyme required in this process. It connected NADH and ubiquinone electron transfer to proton transmembrane transport, which helped to generate proton power, which was necessary for ATP synthesis. Further, ND1, ND2, ND4, ND5, Ndufs1, Ndufs8 were the core subunits of Complex I [37]. In Chlamydomonas, the loss of ND1 or ND5 subunit prevented the assembly of whole mitochondrial Complex I, whereas the loss of ND4 led to the formation of a subcomplex of 650 kDa present in a reduced amount [38,39]. In addition, COX10 related to cytochrome c oxidase (Complex IV) and ATPF 0 B, ATPF 1 D related ATP synthase significantly decreased after the second and fourth days compared to the control group. Complex IV catalyzes the last step of the mitochondrial electron transfer chain and is considered one of the main regulatory sites of oxidative phosphorylation [40]. COX10 participates in the biosynthesis of heme a, a key cofactor of cytochrome c oxidase [41]. ATP synthase can synthesize cellular ATP from ADP and inorganic phosphate. It is comprised of two parts, an ATP-driven motor F 1 and a proton-driven F 0 motor [42].
The down-regulation of oxidative phosphorylation expression led to a reduction in the transmembrane gradient and ATP synthesis, suggesting that B. laterosporus Bl-zj decreased the amount of cellular ATP in M. aeruginosa cells. Moreover, the absolute value of log 2 (fold changes) of MB4 vs. CMA comparison were almost larger than that of MB2 vs. CMA comparison. This indicated that oxidative phosphorylation of M. aeruginosa was continuously inhibited. Further, the ATP synthesis decreased, and the cyanobacteria cells could not obtain enough energy to survive. In addition, the electron transport chain and oxidative phosphorylation are critical cellular processes that sustain life, and the absence of either can affect cellular respiration and even lead to cell death [43]. The decrease of electron transport chain and oxidative phosphorylation genes expression could inhibit the algal growth and even cause death.
The response and adaptability of microalgae to environmental changes were closely related to numerous changes in the composition of lipids in cells and the ability to synthesize a series of special lipids [44]. Under stress, some algal cells cannot reproduce but accumulate lipids instead [45,46]. Many deep-sea bacteria adapted to the deep-sea environment by synthesizing more unsaturated lipids [47]. In this study, the expression of fatty acid biosynthesis-related genes FabH and FabZ were up-regulated in both MB2 and MB4 groups compared to control, while FabF was up-regulated only in MB2 vs. CMA comparison (Table 3). FabF catalyzed the elongation of fatty acyl chains and NADPH-specific reduction of long-chain β-ketoacyl derivatives [48]. FabH was considered to catalyze the first elongation reaction (Claisen condensation) of type II fatty acid synthesis, resulting in the production of short-chain fatty acid primers [49]. FabZ showed ubiquitous distribution in type II fatty acid synthase, and the dehydratase efficiently catalyzed the dehydration of short and long-chain saturated and unsaturated beta-hydroxyacyl-ACPs [50]. The upregulated expression of fatty acid synthesis genes revealed that M. aeruginosa was under stress. When M. aeruginosa was under the stress of pyrogallol, a compound with strong inhibition on M. aeruginosa, the expression of FabZ was also increased [51]. Moreover, the expression of gene Y755 was up-regulated. 2-Cys peroxiredoxin was a subclass of peroxiredoxin proteins (Prx-s), and Prx-s were a ubiquitous family of antioxidant enzymes and important components of the cellular antioxidant defense system as well as redox homeostasis [52]. These suggested that membranes may undergo oxidative damage. In order to maintain the normal function of membranes, more unsaturated fatty acids need to be integrated into membranes.

Conclusions
In the present study, we provide the molecular responses of M. aeruginosa to B. laterosporus treatment. The B. laterosporus treatment led to a suppression of photosynthesis and oxidative phosphorylation of M. aeruginosa. We also found that the expression of the genes related to fatty acid synthesis increased. In addition, the expression of the antioxidant enzymes, such as 2-Cys peroxiredoxin, was increased. The inhibition of photosynthesis and oxidative phosphorylation of M. aeruginosa could block electron transport and affect energy acquisition. The increase in fatty acid synthesis-related genes and 2-Cys peroxiredoxin indicated that cyanobacteria exhibited oxidative damage, causing algal cell membrane lipid peroxidation. These changes inhibited cyanobacteria growth, leading to algae death.