ORIGINAL RESEARCH article

Front. Microbiol., 13 October 2022

Sec. Microbial Physiology and Metabolism

Volume 13 - 2022 | https://doi.org/10.3389/fmicb.2022.1038136

Acclimation and stress response of Prochlorococcus to low salinity

  • 1. Key Laboratory of Tropical Marine Bio-resources and Ecology, South China Sea Institute of Oceanology, Chinese Academy of Sciences, Guangzhou, China

  • 2. College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, Huairou, Beijing, China

  • 3. Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou), Guangzhou, China

Article metrics

View details

4

Citations

2,1k

Views

865

Downloads

Abstract

Prochlorococcus is an obligate marine microorganism and the dominant autotroph in tropical and subtropical open ocean. However, the salinity range for growing and response to low salinity exposure of Prochlorococcus are still unknown. In this study, we found that low-light adapted Prochlorococcus stain NATL1A and high-light adapted strain MED4 could be acclimated in the lowest salinity of 25 and 28 psu, respectively. Analysis of the effective quantum yield of PSII photochemistry (Fv/Fm) indicated that both strains were stressed when growing in salinity lower than 34 psu. We then compared the global transcriptome of low salinity (28 psu) acclimated cells and cells growing in normal seawater salinity (34 psu). The transcriptomic responses of NATL1A and MED4 were approximately different, with more differentially expressed genes in NATL1A (525 genes) than in MED4 (277 genes). To cope with low salinity, NATL1A down-regulated the transcript of genes involved in translation, ribosomal structure and biogenesis and ATP-production, and up-regulated photosynthesis-related genes, while MED4 regulated these genes in an opposite way. In addition, both strains up-regulated an iron ABC transporter gene, idiA, suggesting low salinity acclimated cells could be iron limited. This study demonstrated the growing salinity range of Prochlorococcus cells and their global gene expression changes due to low salinity stress.

Introduction

Cyanobacterium Prochlorococcus is the smallest and most abundant photosynthetic, oxygen-evolving organism on Earth, playing a significant role in carbon fixation and biogeochemical cycles in the ocean (Guillard et al., 1985; Goericke and Welschmeyer, 1993; Liu et al., 1997). The prokaryotic Prochlorococcus cells contain divinyl-chlorophyll a and both monovinyl and divinyl-chlorophyll b as their primary photosynthetic pigments, which are unique to other cyanobacteria that contain chlorophyll a and phycobiliprotein (Chisholm et al., 1992; Hess et al., 1996). Prochlorococcus is believed to be an obligate marine organism that is predominantly found in oligotrophic open oceans, as well as in some coastal waters, but barely seen in low salinity estuarine waters (Flombaum et al., 2013). Prochlorococcus thrives throughout the euphotic zone in the tropical and subtropical oceans from 45° N to 40° S (Scanlan et al., 2009). This genus of marine picocyanobacteria is divided into high-light (HL) adapted and low-light (LL) adapted ecotypes, which are also phylogenetically distinct (Ferris and Palenik, 1998; Moore and Chisholm, 1999). HL ecotypes are usually distributed in upper euphotic zone, while LL ecotypes are generally distributed in the lower to bottom euphotic zone (Johnson et al., 2006; Zinser et al., 2007). Besides the light-related niche partitioning of HL and LL ecotypes, two HL ecotypes, HLI and HLII, also display temperature-related niche partitioning that HLII ecotypes dominate the warmer oceans between 30° N and 30° S while HLI ecotypes dominate the higher latitude oceans (West et al., 2001; Rocap et al., 2003; Mühling, 2012; Voigt et al., 2014). Despite comprising diverse phylogenetic lineages, Prochlorococcus is monophyletic on the phylogenetic tree built on 16S rRNA sequences of cyanobacteria (Rocap et al., 2002). Synechococcus is the sister genus of Prochlorococcus. However, Synechococcus is a provisional genus containing polyphyletic clusters which are scattering on the phylogenetic tree of cyanobacteria (Robertson et al., 2001). Marine Synechococcus is affiliated with cluster 5, which comprises subclusters 5.1, 5.2 and 5.3. In contrast to Prochlorococcus, marine Synechococcus is much more widely distributed, existing from estuary to open ocean and from equatorial to polar regions (Partensky et al., 1999; Zwirglmaier et al., 2008).

Salinity is a crucial factor affecting the growth and biogeography of cyanobacteria (Scanlan et al., 2009). There were plenty of studies on cyanobacteria’s salt acclimation and salt stress response (Hagemann, 2011). However, most of those studies were conducted mainly on freshwater cyanobacteria such as the euryhaline Synechococcus strain PCC 7002 and moderately halotolerant Synechocystis strain PCC 6803 rather than typical marine cyanobacteria such as Prochlorococcus or marine Synechococcus (Hagemann, 2011). For example, when growing at high salinity, Synechococcus PCC 7002 had increased expression of genes involved in compatible solute biosynthesis and electron transport, while only minor changes were observed when cells were grown at low salinity (Ludwig and Bryant, 2012). It also has been revealed that 200–300 genes were up-regulated and a comparable number of genes were down-regulated after the addition of salt in Synechocystis PCC 6803 (Kanesaki et al., 2002; Marin et al., 2003). Secondly, very few studies focus on the acclimation and stress response of marine cyanobacteria to low salinity. Lastly, compared to Synechococcus, salinity-related physiological studies on Prochlorococcus are even more seldom. A recent study showed that Prochlorococcus strain AS9601 could be acclimated to a high salt concentration of 5% (w/v; Al-Hosani et al., 2015). The authors compared the growth rate and transcriptome of AS9601 at salinities 3.8% (w/v) and 5% (w/v), and found that, under high salt concentration, approximately one-third of the genome expressed differentially.

The strict biogeographic distribution of Prochlorococcus in oceanic waters suggests that this organism cannot be adapted to low salinity. However, what is the lowest salinity that Prochlorococcus can survive and what is the stress response of Prochlorococcus cells to low salinity are still unclear. In this study we first tested the salinity range of two Prochlorococcus strains, NATL1A and MED4, and then acclimatized the two strains under different salinities. We found that the lowest acclimation salinity is 25 psu for MED4 and 28 psu for NATL1A. Both NATL1A and MED4 cells were stressed when growing in salinities lower than 34 psu. We also found that the transcriptomic response of the two strains to low salinity stress were highly different.

Materials and methods

Strains and growth conditions

Prochlorococcus strains MED4 and NATL1A were obtained from Jiao Nianzhi Lab, Xiamen University. Cultures were maintained in Pro99 natural seawater medium with a salinity of 34 psu, at 21°C and under a constant light intensity of 10 μE m−2 s−1. We used canted neck polystyrene flasks (Corning Inc., Corning, NY, United States) of different volumes to culture the Prochlorococcus strains.

Experiment setup and growth rate calculation

Preparation of Pro99 medium followed the protocol from the Chisholm Laboratory.1 The seawater from the South China Sea basin was filtered through 0.22 μm polycarbonate membrane, and the salinity was pre-adjusted to 22 psu ~ 60 psu with a 2 psu interval using ddH2O or NaCl. Salinity was measured using an ATAGO PAL-06S refractometer (ATAGO, Japan). These seawaters were autoclaved at 121°C for 15 min. Macronutrient (NH4Cl and NaH2PO4) stocks and the trace metal stock were prepared in advance, and they were added into the above seawater base. Prochlorococcus cultures growing in the Pro99 medium of salinity 34 psu were inoculated into the salinity gradient mediums. The salinity was finally adjusted to 22 psu ~ 60 psu using the ddH2O with Pro99 nutrients. Prochlorococcus growth was monitored every day for 2 weeks by measuring the OD440 absorbance using a multimode plate reader (PerkinElmer, Waltham, MA, United States) and measuring the cell abundance using a flow cytometer (BD Accuri C6, BD Biosciences, CA, United States). Three biological replicates were set up for the experiment. Growth rate was calculated based on the two monitoring methods, respectively. Growth rate was calculated according to Mackey et al. (2013): Td = Ln (Ni + 1/Ni), Ni + 1 is the number of cells on day i + 1, Ni is the number of cells on day i, Td is the growth rate of cells. The average growth rate of cells was calculated during the logarithmic phase.

Low salinity acclimation

Prochlorococcus strains MED4 and NATL1A were acclimated to different salinities (24 psu, 25 psu, 26 psu, 27 psu, 28 psu, 30 psu, 32 psu, 34 psu) by consecutive transfers from exponential growing cultures to fresh media. Three biological replicates were set up for each salinity. Five rounds of transfer were conducted for each strain. Using flow cytometry, cell abundance was monitored at day 0, day 5 and day 10. To assess the stress to low salinity, each strain’s dark-adapted photochemical efficiency (Fv/Fm) was monitored on day 10 in each round, using a handheld fluorometer (AquaPen AP 110/C, Photon Systems Instruments). To measure Fv/Fm, 1 ml culture was dark-adapted in the sample cuvette for 15–30 min. The maximal fluorescence levels (Fm) were measured in the dark and under bright purple light (455 nm, 100 μEm−2 s−1), where F0 is the basal fluorescence level and Fv is the variable fluorescence. The PSII quantum yield was calculated as Fv/Fm = (Fm–F0)/Fm.

RNAseq analysis

To acclimate the Prochlorococcus strains, MED4 and NATL1A were growing in the Pro99 medium of salinity 28 psu and 34 psu for five rounds of inoculation. Then the acclimated cultures were inoculated in fresh medium of salinity 28 psu and 34 psu, with salinity 34 psu being the control. Three biological replicates were set up. During the exponential growth phase, 100 ml cultures were filtered onto 0.22 μm polycarbonate membrane to collect cells and the membranes were immediately flash frozen in RNAlater by liquid nitrogen and stored at −80°C until RNA extraction. Total RNA was extracted from the membrane using the MagZol Reagent (Magen Biotech, Guangzhou, China). Sequencing libraries were prepared using VAHTS™ Stranded mRNA-seq Library Prep Kit for Illumina® (Vazyme biotech co., Ltd., Nanjing, China) following the manufacturer’s instructions. Libraries were multiplexed and sequencing was carried out on an Illumina HiSeq system with the 2 × 150 paired-end (PE) configuration (GENEWIZ). Cutadapt (v1. 9. 1) was used to remove adapters, primers, and reads with a base quality <20 based on FASTQ files. Clean data were aligned to the MED4 and NATL1A genomes via Bowtie2 software (v2. 1. 0). HTSeq (v0. 6. 1p1) was used to estimate gene expression levels from clean data. Differential expression analysis was performed using the DESeq Bioconductor package, a model based on negative binomial distribution. After adjusting using Benjamini and Hochberg’s approach for controlling the false discovery rate, differentially expressed genes were considered significant at value of p < 0.05. Highly induced or suppressed genes were considered as meeting both false discovery rate p < 0.05 and magnitude of log2fold change with values greater than 1 (highly induced) or less than −1 (highly suppressed). These two different criterions were also used in a previous study (Al-Hosani et al., 2015). Transcriptomic data have been deposited in NCBI Gene Expression Omnibus (GEO) under the accession number GSE195946.

Results and discussion

Salinity range and acclimation to low salinity of Prochlorococcus

High-light adapted Prochlorococcus strain MED4 and low-light adapted strain NATL1A were tested for growth in different salinities ranging from 22 psu to 60 psu. Cell counting through flow cytometry (Figure 1A) and absorbance measurement at 440 nm (Figure 1B) were used to monitor the growth of Prochlorococcus cells. MED4 could grow in the salinity range from 22 psu to 50 psu and NATL1A could grow in the range from 26 psu to 50 psu. The optimal salinity ranges of MED4 and NATL1A were similar, from 30 psu to 40 psu. This result came from the first transfer of cultures from salinity 34 psu to other salinities. Under the same growing temperature (21°C) and light intensity (10 μE m−2 s−1), the LL strain NATL1A grew faster than the HL strain MED4.

Figure 1

Figure 1

Growth rate of Prochlorococcus strains NATL1A and MED4 growing in Pro99 medium with salinities from 22 psu to 60 psu. Flow cytometry (A) and absorbance (OD440) measurement of chlorophyll (B) were used to monitor growth.

During the acclimation experiment, NATL1A could grow in salinities higher than 26 psu in the first round, but could not survive in salinity lower than 28 psu in the last round (Figure 2B). Interestingly, MED4 showed a gradually changing growth rate in the salinity gradient from 25 psu to 28 psu, while NATL1A showed a sharp change between salinity 27 psu and salinity 28 psu (Figures 2C,D). The effective quantum yield of PSII photochemistry (Fv/Fm) was measured on the 10th day at the end of each incubation round (Figure 3). Both strains showed reduced yield when growing in low salinities from 24 psu to 32 psu, compared to the yield when growing in salinity 34 psu, and the lower the salinity resulted in lower yield. The yield of MED4 growing in salinity 24 psu was not detectable after round 4, while the yield of NATL1A growing in salinities 27 psu and below was not detectable after round 2. Together, these data showed that Prochlorococcus MED4 and NATL1A could be acclimated in salinities 25 psu and 28 psu, respectively. Interesting, the high-light adapted strain MED4 and low-light adapted strain NATL1A showed different tolerance to low salinity.

Figure 2

Figure 2

Acclimation of Prochlorococcus strains NATL1A and MED4 in different salinities. Five rounds of transfers were carried out and the growth (A,B) were monitored by flow cytometry. The pictures (C,D) showed the last round of cultures.

Figure 3

Figure 3

Dark-adapted photochemical efficiency (Fv/Fm) of Prochlorococcus strains MED4 (A) and NATL1A (B) growing in different salinities during acclimation.

It is well known that Prochlorococcus is an oceanic microorganism (Partensky et al., 1999), although a few studies claimed that Prochlorococcus-like populations existed in estuarine and even freshwater environments (Corzo et al., 1999; Shang et al., 2007; Mitbavkar et al., 2012; Zhang et al., 2013). However, these studies all only depended on flow cytometry investigation, and could not confirm that those “populations” on the flow cytometry diagram were indeed Prochlorococcus. Our acclimation study suggests that Prochlorococcus cannot live in salinity lower than 25 psu for a long time period (50 days in this study). This study provides evidence supporting that Prochlorococcus is an oceanic organism.

Differentially expressed genes in low salinity acclimated Prochlorococcus cells

RNA-seq was performed to assess the response of acclimated Prochlorococcus cells to low salinity (28 psu), with the salinity 34 psu being the control. When the filter criteria of significance meet the value of p < 0.05, there were 525 differentially expressed genes in the low salinity acclimated cells of NATL1A, with 286 genes being induced and 239 genes being repressed (Table 1). By contrast, MED4 appears to be less fluctuant under low salinity stress, with only 277 differentially expressed genes, among which 146 were induced and 131 were repressed (Table 1). A previous study compared the transcriptomes of Prochlorococcus AS9601 under high salt stress (5.0%) and under normal salt concentration (3.8%), and found 627 differentially expressed genes (Al-Hosani et al., 2015). Together, these results suggest that Prochlorococcus is sensitive to salinity changes.

Table 1

Function categoriesTotal no.Total no. differential expressedNo. inducedNo. repressedPrevalent expression profile
NATL1AMED4NATL1AMED4NATL1AMED4NATL1AMED4NATL1AMED4
Translation, ribosomal structure and biogenesis1281303825524331RepressedInduced
Transcription26271067531InducedInduced
Signal transduction mechanisms2624733241EqualEqual
Secondary metabolites biosynthesis, transport, and catabolism1822473215EqualRepressed
Replication, recombination, and repair44212100EqualEqual
Posttranslational modification, protein turnover, chaperones878223191351014InducedRepressed
Nucleotide transport and metabolism48501458461EqualInduced
Lipid metabolism29311063373RepressedEqual
Intracellular trafficking and secretion1412422121EqualEqual
Inorganic ion transport and metabolism5857865234EqualEqual
Energy production and conversion7778281419896InducedEqual
DNA replication, recombination, and repair68651277552EqualInduced
Defense mechanisms1716212001EqualEqual
Coenzyme metabolism105101201164147RepressedRepressed
Cell wall/ membrane/envelope biogenesis12111001EqualEqual
Cell envelope biogenesis, outer membrane93913010166144EqualEqual
Cell division and chromosome partitioning1416433112EqualEqual
Cell cycle control, cell division, chromosome partitioning22010100EqualEqual
Carbohydrate transport and metabolism464810104664EqualEqual
Amino acid transport and metabolism12811838171911196EqualInduced
Function unknown8688181310489
General function prediction only1371324413238215
Others65200020
Not in COGs1,02184119696125437153
Total2,2392042525277286146239131

Functional categorization of differentially expressed genes (p < 0.05) in low salinity acclimated cells of NATL1A and MED4.

Subsequently, the differentially expressed genes of these two strains were functionally classified according to Cyanobase definitions (Fujisawa et al., 2014). Firstly, the numbers of induced and repressed genes were equal for most functional modules (60–70%) in each of the two strains. Secondly, compared to the control group, the changed prevalent expression profiles between the two strains were different. In low salinity acclimated NATL1A, some genes involved in translation, ribosomal structure and biogenesis, lipid metabolism and coenzyme metabolism were down-regulated, while genes involved in transcription, posttranslational modification, protein turnover, chaperones and energy production and conversion were up-regulated (Table 1). However, in MED4, regulation profile of the functions mentioned above is different from NATL1A, except for coenzyme metabolism.

Contrasting regulation between NATL1A and MED4

Most interestingly, among the genes involved in translation, ribosomal structure and biogenesis, five were up-regulated and 33 were down-regulated in low salinity acclimated NATL1A cells compared to control (Tables 1, 2). However, in MED4, 24 genes of those genes were up-regulated and only one was down-regulated (Tables 1, 3). Strikingly, the regulation of genes involved in energy production and conversion were also in distinct patterns between NATL1A and MED4. In low salinity acclimated NATL1A cells, the ATP-producing genes were down-regulated (atpA, atpC, atpD, atpH and other ATP synthase genes), while many genes involved in photosynthesis (psaC, psb27, rbcS), cytochrome oxidation (cyoA, cyoB, ctaE), NADH dehydrogenase (ndhA, ndhH) were up-regulated (Table 2). However, in low salinity acclimated MED4 cells, genes for photosynthesis were down-regulated, such as photosystems II (psbA, psbB, psbD, psbN), cytochrome F (petA), and electron transport chain intermediate (ndhD), while most ATP-producing genes (acnB, atpG, atpF, atpH, atpD) were up-regulated. This striking contrasting transcriptional regulation indicated the two strains processed different response mechanisms to low salinity stress. It is likely that, to respond to low salinity stress, NATL1A enhanced photosynthesis but repressed ATP production and translation and biosynthesis. In contrast, MED4 repressed photosynthesis but enhanced ATP production, translation and biosynthesis. The reason is possible that NATL1A and MED4 were in different stress level under the salinity 28 psu, which appears to be slightly stressful for MED4, but extremely stressful for NATL1A. This is the reason why the differentially expressed genes of NATL1A were more than those of MED4.

Table 2

Gene IDGene nameProductp-Valuelog2FC
Energy production and conversion
gene-NATL1_05001cyoAputative cytochrome c oxidase, subunit 2< 0.0011.249
gene-NATL1_17081acoAPyruvate dehydrogenase E1 alpha subunit< 0.0010.979
gene-NATL1_04991cyoBCytochrome c oxidase, subunit I< 0.0010.897
gene-NATL1_06051rbcSRibulose bisphosphate carboxylase, small chain0.0100.780
gene-NATL1_02471ndhHputative NADH dehydrogenase subunit< 0.0010.770
gene-NATL1_20041NATL1_20041NADH dehydrogenase I subunit N< 0.0010.723
gene-NATL1_05651psb27possible Photosystem II reaction center Psb27 protein< 0.0010.6925
gene-NATL1_05981chlNLight-independent protochlorophyllide reductase subunit N< 0.0010.650
gene-NATL1_20591psaCPhotosystem I subunit PsaC0.0040.633
gene-NATL1_04561pdhCDihydrolipoamide acetyltransferase0.0010.607
gene-NATL1_04981ctaECytochrome c oxidase, subunit III0.0020.580
gene-NATL1_20451icdIsocitrate dehydrogenase0.0020.547
gene-NATL1_04171petBCytochrome b60.0060.532
gene-NATL1_02331ndhAputative respiratory-chain NADH dehydrogenase subunit0.0060.498
gene-NATL1_17231NATL1_17231FAD/FMN-containing dehydrogenases0.0240.484
gene-NATL1_03751rubprobable rubredoxin0.01580.453
gene-NATL1_21811acnBAconitate hydratase B0.01880.394
gene-NATL1_03311psbIphotosystem II reaction center PsbI protein0.0263−0.446
gene-NATL1_19381NATL1_19381Fe-S oxidoreductase0.021−0.482
gene-NATL1_18501atpHATP synthase, delta (OSCP) subunit0.008−0.497
gene-NATL1_18491atpAATP synthase F1, alpha subunit0.006−0.565
gene-NATL1_19601psaIphotosystem I subunit VIII (PsaI)0.008−0.574
gene-NATL1_00561NATL1_00561Flavoprotein, FldA0.009−0.647
gene-NATL1_18481NATL1_18481ATP synthase gamma subunit< 0.001−0.668
gene-NATL1_18511NATL1_18511ATP synthase B/B′ CF(0)< 0.001−0.671
gene-NATL1_18381atpDATP synthase F1, beta subunit< 0.001−0.723
gene-NATL1_14931gldAputative glycerol dehydrogenase0.004−0.755
gene-NATL1_18391atpCATP synthase, Epsilon subunit< 0.001−1.089
Inorganic ion transport and metabolism
gene-NATL1_16181afuAputative iron ABC transporter, substrate binding protein< 0.0011.208
gene-NATL1_19031NATL1_19031Ferric uptake regulator family< 0.0010.853
gene-NATL1_05281nhaPputative Na+/H+ antiporter, CPA1 family0.0050.597
gene-NATL1_20831mgtEMgtE family, putative magnesium transport protein0.0090.501
gene-NATL1_03411amtBAmmonium transporter family0.0220.414
gene-NATL1_03071met3ATP-sulfurylase0.002−0.582
gene-NATL1_15081petHferredoxin-NADP oxidoreductase (FNR)< 0.001−0.687
Molecular chaperone
gene-NATL1_09851NATL1_09851Molecular chaperone DnaK, heat shock protein hsp700.0010.668
gene-NATL1_21861NATL1_21861Molecular chaperone DnaK2, heat shock protein hsp70-20.0100.624
Translation, ribosomal structure and biogenesis
gene-NATL1_18631NATL1_18631FtsJ cell division protein: S4 domain:Hemolysin A0.0010.837
gene-NATL1_04781NATL1_04781tRNA/rRNA methyltransferase (SpoU)0.0060.573
gene-NATL1_00131NATL1_00131tRNA-dihydrouridine synthase0.0400.566
gene-NATL1_04521lrtAlight repressed protein A-like protein0.0020.563
gene-NATL1_03171ileSIsoleucyl-tRNA synthetase0.0160.425
gene-NATL1_04021rps1a30S ribosomal protein S1, protein A0.036−0.379
gene-NATL1_17711rplU50S ribosomal protein L210.025−0.445
gene-NATL1_16641rpsN30S Ribosomal protein S140.011−0.455
gene-NATL1_05781frrRibosome recycling factor0.022−0.457
gene-NATL1_19971rpsC30S ribosomal protein S30.016−0.480
gene-NATL1_07951glySGlycyl-tRNA synthetase beta subunit0.013−0.485
gene-NATL1_19921rplX50S ribosomal protein L240.022−0.490
gene-NATL1_07891rpsB30S ribosomal protein S20.015−0.497
gene-NATL1_02771rplL50S ribosomal protein L7/L120.004−0.564
gene-NATL1_19481rpsJ30S ribosomal protein S100.001−0.568
gene-NATL1_21621aspSAspartyl-tRNA synthetase0.005−0.570
gene-NATL1_19521rpsL30S ribosomal protein S120.001−0.582
gene-NATL1_19871rpsE30S ribosomal protein S50.001−0.605
gene-NATL1_02781rplJ50S ribosomal protein L100.002−0.610
gene-NATL1_19991rpsS30S Ribosomal protein S190.001−0.616
gene-NATL1_19891rplF50S ribosomal protein L6< 0.001−0.623
gene-NATL1_09331gatAGlutamyl-tRNA (Gln) amidotransferase A subunit0.001−0.631
gene-NATL1_05331mapputative methionine aminopeptidase< 0.001−0.635
gene-NATL1_16221glyQglycyl-tRNA synthetase, alpha subunit0.034−0.637
gene-NATL1_20021rplD50S ribosomal protein L4< 0.001−0.705
gene-NATL1_19881rplR50S ribosomal protein L18< 0.001−0.736
gene-NATL1_19951rpmC50S ribosomal protein L29< 0.001−0.752
gene-NATL1_00581alaSAlanyl-tRNA synthetase0.002−0.756
gene-NATL1_10481fmtputative Methionyl-tRNA formyltransferase0.008−0.777
gene-NATL1_17561tyrSTyrosyl-tRNA synthetase0.026−0.777
gene-NATL1_03281pthPeptidyl-tRNA hydrolase0.050−0.785
gene-NATL1_20011rplW50S ribosomal protein L23< 0.001−0.791
gene-NATL1_07901tsfputative Elongation factor Ts< 0.001−0.821
gene-NATL1_19511rpsG30S ribosomal protein S7< 0.001−0.868
gene-NATL1_21311rplT50S ribosomal protein L20< 0.001−0.883
gene-NATL1_06131tdcFPutative translation initiation inhibitor, yjgF family0.001−0.939
gene-NATL1_10131rpsR30S Ribosomal protein S180.001−1.140
gene-NATL1_10191cspRputative tRNA/rRNA methyltransferase (SpoU family)0.009−1.690

List of a part of differentially expressed genes (p < 0.05) in low salinity acclimated Prochlorococcus NATL1A.

Table 3

Gene IDGene nameProductp-Valuelog2FC
Energy production and conversion
gene-PMM0930pdhBPyruvate dehydrogenase E1 beta subunit< 0.0010.691
gene-PMM0317psbMpossible Photosystem II reaction center M protein (PsbM)0.0330.590
gene-PMM0544chlBLight-independent protochlorophyllide reductase subunit B0.0040.507
gene-PMM1452atpH,atpDATP synthase, delta (OSCP) subunit0.0420.476
gene-PMM0785prk,cbbPphosphoribulokinase0.0070.475
gene-PMM1700acnBAconitate hydratase B0.0090.460
gene-PMM1454atpGATP synthase B/B′ CF(0)0.0190.411
gene-PMM1453atpFATP synthase B/B′ CF(0)0.0380.392
gene-PMM0223psbAPhotosystem II PsbA protein (D1)0.018−0.427
gene-PMM1157psbDPhotosystem II PsbD protein (D2)0.008−0.470
gene-PMM0315psbBPhotosystem II PsbB protein (CP47)0.014−0.476
gene-PMM0461petACytochrome f0.021−0.477
gene-PMM1171isiBFlavodoxin0.045−0.522
gene-PMM1229PMM1229Dehydrogenase, E1 component0.009−0.577
gene-PMM0594ndhDputative NADH Dehydrogenase (complex I) subunit (chain 4)0.001−0.663
gene-PMM0252psbNPhotosystem II reaction center N protein (psbN)0.001−0.919
gene-PMM0366PMM0366Type-1 copper (blue) domain0.002−0.960
gene-PMM0316PMM0316possible ferredoxin< 0.001−1.358
gene-PMM0926psb28possible Photosystem II reaction center Psb28 protein0.041−1.707
Inorganic ion transport and metabolism
gene-PMM1032PMM1032ABC transporter, substrate binding protein, possibly Mn.0.0060.751
gene-PMM1164futA/afuA/idiAputative iron ABC transporter, substrate binding protein< 0.0010.950
gene-PMM0808PMM0808Rieske iron–sulfur protein 2Fe-2S subunit0.019−0.501
gene-PMM0227cysDATP-sulfurylase< 0.001−0.893
gene-PMM1701PMM1701putative chloride channel< 0.001−0.984
gene-PMM0504PMM0504CutA1 divalent ion tolerance protein0.004−2.817
Molecular chaperone
gene-PMM1704dnaK2Molecular chaperone DnaK2, heat shock protein hsp70-20.0470.360
Translation, ribosomal structure and biogenesis
gene-PMM1537rps13, rpsM30S ribosomal protein S130.0011.774
gene-PMM1538rpmJ, rpl3650S Ribosomal protein L36< 0.0011.081
gene-PMM1688aspSAspartyl-tRNA synthetase0.0020.952
gene-PMM1507rpsJ, rps1030S ribosomal protein S10< 0.0010.844
gene-PMM0068defputative formylmethionine deformylase0.0250.745
gene-PMM1661rpl35, rpmI50S ribosomal protein L350.0060.724
gene-PMM1534rpl17, rplQ50S ribosomal protein L170.0020.697
gene-PMM1550rpl29, rpmC50S ribosomal protein L290.0410.609
gene-PMM1545rps8, rpsH30S ribosomal protein S80.0010.599
gene-PMM1191pnppolyribonucleotide nucleotidyltransferase0.0030.592
gene-PMM0597thrSThreonyl-tRNA synthetase0.0160.561
gene-PMM0312rps1a, rpsA130S ribosomal protein S1, homolog A0.0020.560
gene-PMM1548rpl14, rplN50S Ribosomal protein L140.0470.539
gene-PMM1280PMM1280putative bifuntional enzyme: tRNA methyltransferase: 2-C-methyl-D-erythritol 2, 4-cyclodiphosphate synthase0.0290.551
gene-PMM0202rpl10, rplJ50S ribosomal protein L100.0020.507
gene-PMM1662rpl20, rplT50S ribosomal protein L200.0170.498
gene-PMM0870rpl33, rpmG50S Ribosomal protein L330.0300.467
gene-PMM1706rps6, rpsF30S ribosomal protein S60.0410.465
gene-PMM1508tufAElongation factor Tu0.0060.459
gene-PMM0238ileSIsoleucyl-tRNA synthetase0.0270.429
gene-PMM1532rpl13, rplM50S ribosomal protein L130.0440.425
gene-PMM0203rpl1, rplA50S ribosomal protein L10.0210.420
gene-PMM1546rpl5, rplE50S ribosomal protein L50.0280.384
gene-PMM1509fusAElongation factor G0.0300.370

List of a part of differentially expressed genes (p < 0.05) in low salinity acclimated Prochlorococcus MED4.

It has been pointed out that the response of photosystem gene expression to high salt stress might be dependent on the organism under study, based on the investigations on Prochlorococcus strain AS9601, Synechocystis PCC 6803 and Synechococcus PCC 7002 (Al-Hosani et al., 2015). In high salt acclimated AS9601 cells, many genes coding for components of Photosystem I, Photosystem II and chlorophyll were down-regulated. By contrast, in high salt acclimated PCC 7002, PSI genes were down-regulated but PSII genes were not changed significantly (Ludwig and Bryant, 2012). Similarly, in this study, NATL1A and MED4 also showed heterogeneity in response to low salinity stress.

Compatible solute and transporters

Cyanobacteria generally use the salt-out strategy for salt acclimation, in which cells maintain low intracellular ion concentration and accumulate compatible solutes to establish turgor (Hagemann, 2011). Compatible solutes are low-molecular-weight organic compounds, with sucrose, glucosylglycerol (GG), glucosylglycerate (GGA) and glycine betaine (GB) being the most common ones utilized by cyanobacteria (Klähn and Hagemann, 2011). Prochlorococcus cells probably use GGA and sucrose as their main compatible solutes (Scanlan et al., 2009). In NATL1A, we observed significant decrease in transcript abundance of the GGA synthesis genes (gpgP, encoding glucosyl-phosphoglycerate phosphatase, and gpgS, encoding glucosyl-phosphoglycerate synthase) in the low salinity acclimated cells compared to control cells (Table 4). However, we did not observe significant change on the sucrose synthesis gene spsA (encoding sucrose phosphate synthase). Moreover, in MED4, all the three genes did not show significant change in transcript abundance. These results suggest that, to cope with low salinity stress, NATL1A probably reduced the concentration of intracellular compatible solute GGA, while MED4 did not reduce the concentration of the compatible solutes. Again, this different observations may be due to that the two strains were at different stress level when growing in the medium with salinity 28 psu. In another study, high salt acclimated Prochlorococcus AS9601 cells up-regulated the gpgS gene and a sodium transporter, suggesting that active extrusion of sodium ions and accumulation of GGA are involved in AS9601 acclimation to high salt stress (Al-Hosani et al., 2015). Together, these results suggests that compatible solute GGA may play an important role in the adaptation of Prochlorococcus to salinity changes.

Table 4

NATL1AMED4
Locusp-Valuelog2FCLocusp-Valuelog2FC
gpgPNATL1_057210.001−0.809PMM05150.3670.214
gpgSNATL1_091310.012−0.431PMM09620.476−0.132
spsANATL1_219510.421−0.147PMM17110.143−0.371

Expression change on genes responsible for compatible solute biosynthesis.

Na+/H+ antiporter is closely related to plant salinity tolerance, and it is one of the critical factors of plant salt tolerance. To adapt to a high salt environment, plants will reduce the plasma membrane Na+ level through Na+/H+ antiporter (Apse et al., 1999; Hasegawa et al., 2000). Besides, cyanobacteria cells involved in salt stress tolerance was correlated with the activity of Na+/H+ antiporter (Allakhverdiev et al., 1999, 2000). However, in this study, the transcript level of Na+/H+ antiporter (nhaP) was increased under low salinity stress in NATL1A cells (Table 2). It is not clear what is the mechanism involved in this phenomenon. Perhaps the increasing expression of Na+/H+ antiporter would help to reduce the cytoplasm Na+ level which has already adapted to high salinity level of seawater.

Iron transporter and molecular chaperone

Interestingly, a periplasmic ABC-type Fe3+ transporter (afuA/idiA/futA) was up-regulated in low salinity acclimated cells of both NATL1A and MED4, compared to control cells. Moreover, NATL1A also up-regulated a ferric uptake regulator (NATL1_19031). It has been demonstrated that the transcript levels of idiA gene in Synechococcus PCC 6301 and Prochlorococcus MED4 were increased under iron deficiency conditions (Michel et al., 1999; Webb et al., 2001; Thompson et al., 2011). This suggests that cells may be iron-limited under low salt-stress. Previously, afuA was found to be down-regulated in high salt stressed cells of Prochlorococcus AS9601 (Al-Hosani et al., 2015). The authors attributed this to the reduced expression of iron required proteins under high salt condition. They also concluded that AS9601 was not iron limited because no difference in ferredoxin expression level was found between salt acclimated cells and control cells. It has been also revealed that iron requirement and siderophore production in cells is lower under high salinity (Boyle et al., 1977; Ruebsam et al., 2018). Together, these results indicates that there is a tight link between iron requirement and salt conditions in Prochlorococcus. However, the gene isiB (flavodoxin), which was induced in low iron stress (Erdner and Anderson, 1999; McKay et al., 1999), was not observed to be up-regulated in this study (Table 3). Hence, the specific relationship between low salinity stress and iron homeostasis remains to be investigated.

Up-regulated expression of dnaK was observed in both MED4 and NATL1A, which suggests that this gene could play a role in low salinity acclimation (Tables 2, 3). However, the molecular chaperone dnaK is one of the key factors for salt stress tolerance in halophiles, and over expression of dnaK can greatly reduce the growth lag period of the bacteria, allowing them to grow normally under salt stress (Sugimoto et al., 2003). Fukuda et al. (2001, 2002) cloned the dnaK gene from Tetragenococcus halophila JCM5888 and introduced it into E. coli, and found that the dnaK transcript abundance was increased approximately 3.5-fold under salt stress. Meanwhile, dnaK was also found to be present in the halotolerant cyanobacterium Aphanothece halophytica (Hibino et al., 1999). The gene product of dnaK, heat shock protein hsp70, likely plays an important role in stress resistance, no matter it is low salinity stress or high salt stress.

Highly differentially expressed genes

When the filter criterion was changed from only meeting the value of p (p < 0.05) to meeting both value of p and log2fold change with values greater than 1 (high induction) or less than-1 (high inhibition), there were 81 and 30 highly differentially expressed genes in NATL1A and MED4, respectively. These number are comparable to the previous study on Prochlorococcus AS9601, in which 69 highly differentially expressed genes were found in high salt acclimated cells compared to control cells (Al-Hosani et al., 2015). In NATL1A, 22 genes were down-regulated and 59 were up-regulated, while in MED4, 17 genes were down-regulated and 13 were up-regulated. There was no apparent gene enrichment pattern observed among these highly differentially expressed genes (Supplementary Tables S1, S2). For example, in low-salinity stress cells of NATL1A, many genes were highly inhibited, which were related to posttranslational modification (NATL1_02111 and NATL1_13731), signal transduction mechanisms (typA), cell envelope biogenesis, outer membrane (NATL1_08371 and NATL1_04491), translation, ribosomal structure and biogenesis (rpsR), coenzyme metabolism (folE), energy production and conversion (atpC), and amino acid transport and metabolism (proA). Nevertheless, in salinity acclimated cells of MED4, some other genes appear to be repressed, which were those involved in energy production and conversion (PMM0316), secondary metabolites biosynthesis, transport, and catabolism (PMM0280), DNA replication, recombination, and repair (ruvC), lipid metabolism (des, yocE) and posttranslational modification (PMM1006).

Conclusion

Prochlorococcus is the most abundant phototroph in the ocean. This organism has been adapted to open ocean areas with stable salt concentrations, and barely found in nearshore and estuarine waters with lower and variable salt concentrations. In this study, we showed that the lowest salinities for acclimation of high-light adapted Prochlorococcus strain MED4 and low-light adapted strain NATL1A were 25 psu and 28 psu, respectively. The optimal growing salinity of both MED4 and NATL1A were from 30 to 40 psu. Global transcriptome analysis showed that the two strains responded differently to low salinity stress. First, far more genes of NATL1A were impacted than those of MED4 in low salinity acclimated cells, suggesting NATL1A was more intensively stressed than MED4 under salinity 28 psu. Second, compared to control, low salinity acclimated cells of NATL1A repressed the expression of genes involved in translation, ribosomal structure and biogenesis and ATP production, but enhanced photosynthesis, while MED4 regulated these pathways in an opposite way. To cope with low salinity, NATL1A also reduced the transcript abundance of genes involved in compatible solute GGA, while MED4 did not. Interpreting from a previous study and this study, a tight link between iron transportation and salt condition was verified, with high salinity stressed cells coupling with up-regulation of iron transporters and low salinity stressed cells coupling with down-regulation of iron transporters. This study demonstrated the regulations of global transcriptome of Prochlorococcus under low salinity stress and the mechanisms within those regulations warrant further investigation.

Funding

This study was funded by the National Natural Science Foundation of China [grant number 42176116, 41576126]; Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou) [GML2019ZD0404]; the Natural Science Foundation of Guangdong Province [grant number 2017A030306020]; the Youth Innovation Promotion Association [grant number 2018377]; and the Rising Star Foundation of the South China Sea Institute of Oceanology [grant number NHXX2019ST0101].

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE195946.

Author contributions

SH, JD, and LL designed the experiments. XH and HL performed the experiments and analyzed the data. SH and XH wrote the manuscript. JD and LL provided resources and supervision. All authors contributed to the article and approved the submitted version.

Acknowledgments

We thank Rui Zhang and Xilin Xiao for sharing the Prochlorococcus strains.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.1038136/full#supplementary-material

References

  • 1

    Al-HosaniS.OudahM. M.HenschelA.YousefL. F. (2015). Global transcriptome analysis of salt acclimated Prochlorococcus AS9601. Microbiol. Res.176, 2128. doi: 10.1016/j.micres.2015.04.006

  • 2

    AllakhverdievS. I.NishiyamaY.SuzukiI.TasakaY.MurataN. (1999). Genetic engineering of the unsaturation of fatty acids in membrane lipids alters the tolerance of Synechocystis to salt stress. Proc. Natl. Acad. Sci.96, 58625867. doi: 10.1073/pnas.96.10.5862

  • 3

    AllakhverdievS. I.SakamotoA.NishiyamaY.InabaM.MurataN. (2000). Ionic and osmotic effects of NaCl-induced inactivation of photosystems I and II in Synechococcus sp. Plant Physiol.123, 10471056. doi: 10.1104/pp.123.3.1047

  • 4

    ApseM. P.AharonG. S.SneddenW. A.BlumwaldE. (1999). Salt tolerance conferred by overexpression of a vacuolar Na+/H+ antiport in Arabidopsis. Science285, 12561258. doi: 10.1126/science.285.5431.1256

  • 5

    BoyleE. A.EdmondJ. M.SholkovitzE. R. (1977). The mechanism of iron removal in estuaries. Geochim. Cosmochim. Acta41, 13131324. doi: 10.1016/0016-7037(77)90075-8

  • 6

    ChisholmS. W.FrankelS. L.GoerickeR.OlsonR. J.PalenikB.WaterburyJ. B.et al. (1992). Prochlorococcus marinus nov. gen. Nov. sp.: an oxyphototrophic marine prokaryote containing divinyl chlorophyll a and b. Arch. Microbiol.157, 297300. doi: 10.1007/BF00245165

  • 7

    CorzoA.Jimenez-GomezF.GordilloF. J. L.Garcia-RuizR.NiellF. X. (1999). Synechococcus and Prochlorococcus-like populations detected by flow cytometry in a eutrophic reservoir in summer. J. Plankton Res.21, 15751581. doi: 10.1093/plankt/21.8.1575

  • 8

    ErdnerD. L.AndersonD. M. (1999). Ferredoxin and flavodoxin as biochemical indicators of iron limitation during open-ocean iron enrichment. Limnol. Oceanogr.44, 16091615. doi: 10.4319/lo.1999.44.7.1609

  • 9

    FerrisM. J.PalenikB. (1998). Niche adaptation in ocean cyanobacteria. Nature396, 226228. doi: 10.1038/24297

  • 10

    FlombaumP.GallegosJ. L.GordilloR. A.RinconJ.ZabalaL. L.JiaoN. Z.et al. (2013). Present and future global distributions of the marine cyanobacteria Prochlorococcus and Synechococcus. Proc. Natl. Acad. Sci.110, 98249829. doi: 10.1073/pnas.1307701110

  • 11

    FujisawaT.OkamotoS.KatayamaT.NakaoM.YoshimuraH.Kajiya-KanegaeH.et al. (2014). CyanoBase and RhizoBase: databases of manually curated annotations for cyanobacterial and rhizobial genomes. Nucleic Acids Res.42, D666D670. doi: 10.1093/nar/gkt1145

  • 12

    FukudaD.WatanabeM.ShinoS.SonomotoK.IshizakiA. (2001). Cloning and characterization of dnaK operon of Tetragenococcus halophila. J. Fac. Agric. Kyushu Univ.46, 229241. doi: 10.5109/24435

  • 13

    FukudaD.WatanabeM.SonezakiS.SugimotoS.SonomotoK.IshizakiA. (2002). Molecular characterization and regulatory analysis of dnaK operon of halophilic lactic acid bacterium Tetragenococcus halophila. J. Biosci. Bioeng.93, 388394. doi: 10.1016/S1389-1723(02)80072-X

  • 14

    GoerickeR.WelschmeyerN. A. (1993). The marine prochlorophyte Prochlorococcus contributes significantly to phytoplankton biomass and primary production in the Sargasso Sea. Deep-Sea Res. I Oceanogr. Res. Pap.40, 22832294. doi: 10.1016/0967-0637(93)90104-B

  • 15

    GuillardR. R. L.MurphyL. S.FossP.Liaaen-JensenS. (1985). Synechococcus spp. as likely zeaxanthin-dominant ultraphytoplankton in the North Atlantic. Limnol. Oceanogr.30, 412414. doi: 10.4319/lo.1985.30.2.0412

  • 16

    HagemannM. (2011). Molecular biology of cyanobacterial salt acclimation. FEMS Microbiol. Rev.35, 87123. doi: 10.1111/j.1574-6976.2010.00234.x

  • 17

    HasegawaP. M.BressanR. A.ZhuJ. K.BohnertH. J. (2000). Plant cellular and molecular responses to high salinity. Annu. Rev. Plant Biol.51, 463499. doi: 10.1146/annurev.arplant.51.1.463

  • 18

    HessW. R.PartenskyF.Van der StaayG. W.Garcia-FernandezJ. M.BörnerT.VaulotD. (1996). Coexistence of phycoerythrin and a chlorophyll a/b antenna in a marine prokaryote. Proc. Natl. Acad. Sci.93, 1112611130. doi: 10.1073/pnas.93.20.11126

  • 19

    HibinoT.KakuN.YoshikawaH.TakabeT.TakabeT. (1999). Molecular characterization of DnaK from the halotolerant cyanobacterium Aphanothece halophytica for ATPase, protein folding, and copper binding under various salinity conditions. Plant Mol. Biol.40, 409418. doi: 10.1023/A:1006273124726

  • 20

    JohnsonZ. I.ZinserE. R.CoeA.McNultyN. P.WoodwardE. M. S.ChisholmS. W. (2006). Niche partitioning among Prochlorococcus ecotypes along ocean-scale environmental gradients. Science311, 17371740. doi: 10.1126/science.1118052

  • 21

    KanesakiY.SuzukiI.AllakhverdievS. I.MikamiK.MurataN. (2002). Salt stress and hyperosmotic stress regulate the expression of different sets of genes in Synechocystis sp. PCC 6803. Biochem. Biophys. Res. Commun.290, 339348. doi: 10.1006/bbrc.2001.6201

  • 22

    KlähnS.HagemannM. (2011). Compatible solute biosynthesis in cyanobacteria. Environ. Microbiol.13, 551562. doi: 10.1111/j.1462-2920.2010.02366.x

  • 23

    LiuH.NollaH. A.CampbellL. (1997). Prochlorococcus growth rate and contribution to primary production in the equatorial and subtropical North Pacific Ocean. Aquat. Microb. Ecol.12, 3947. doi: 10.3354/ame012039

  • 24

    LudwigM.BryantD. A. (2012). Synechococcus sp. strain PCC 7002 transcriptome: acclimation to temperature, salinity, oxidative stress, and mixotrophic growth conditions. Front. Microbiol.3:354. doi: 10.3389/fmicb.2012.00354

  • 25

    MackeyK. R.PaytanA.CaldeiraK.GrossmanA. R.MoranD.McIlvinM.et al. (2013). Effect of temperature on photosynthesis and growth in marine Synechococcus spp. Plant Physiol.163, 815829. doi: 10.1104/pp.113.221937

  • 26

    MarinK.SuzukiI.YamaguchiK.RibbeckK.YamamotoH.KanesakiY.et al. (2003). Identification of histidine kinases that act as sensors in the perception of salt stress in Synechocystis sp. PCC 6803. Proc. Natl. Acad. Sci.100, 90619066. doi: 10.1073/pnas.1532302100

  • 27

    McKayR. M. L.La RocheJ.YakuninA. F.DurnfordD. G.GeiderR. J. (1999). Accumulation of ferredoxin and flavodoxin in a marine diatom in response to Fe. J. Phycol.35, 510519. doi: 10.1046/j.1529-8817.1999.3530510.x

  • 28

    MichelK. P.KrügerF.PühlerA.PistoriusE. K. (1999). Molecular characterization of idiA and adjacent genes in the cyanobacteria Synechococcus sp. strains PCC 6301 and PCC 7942. Microbiology145, 14731484. doi: 10.1099/13500872-145-6-1473

  • 29

    MitbavkarS.RajaneeshK. M.AnilA. C.SundarD. (2012). Picophytoplankton community in a tropical estuary: detection of Prochlorococcus-like populations. Estuar. Coast. Shelf Sci.107, 159164. doi: 10.1016/j.ecss.2012.05.002

  • 30

    MooreL. R.ChisholmS. W. (1999). Photophysiology of the marine cyanobacterium Prochlorococcus: ecotypic differences among cultured isolates. Limnol. Oceanogr.44, 628638. doi: 10.4319/lo.1999.44.3.0628

  • 31

    MühlingM. (2012). On the culture-independent assessment of the diversity and distribution of Prochlorococcus. Environ. Microbiol.14, 567579. doi: 10.1111/j.1462-2920.2011.02589.x

  • 32

    PartenskyF.BlanchotJ.VaulotD. (1999). Differential distribution and ecology of Prochlorococcus and Synechococcus in oceanic waters: a review. Bull. Inst. Oceanogr. Monaco19, 457476.

  • 33

    RobertsonB. R.TezukaN.WatanabeM. M. (2001). Phylogenetic analyses of Synechococcus strains (cyanobacteria) using sequences of 16S rDNA and part of the phycocyanin operon reveal multiple evolutionary lines and reflect phycobilin content. Int. J. Syst. Evol. Microbiol.51, 861871. doi: 10.1099/00207713-51-3-861

  • 34

    RocapG.DistelD. L.WaterburyJ. B.ChisholmS. W. (2002). Resolution of Prochlorococcus and Synechococcus ecotypes by using 16S-23S ribosomal DNA internal transcribed spacer sequences. Appl. Environ. Microbiol.68, 11801191. doi: 10.1128/AEM.68.3.1180-1191.2002

  • 35

    RocapG.LarimerF. W.LamerdinJ.MalfattiS.ChainP.AhlgrenN. A.et al. (2003). Genome divergence in two Prochlorococcus ecotypes reflects oceanic niche differentiation. Nature424, 10421047. doi: 10.1038/nature01947

  • 36

    RuebsamH.KirschF.ReimannV.ErbanA.KopkaJ.HagemannM.et al. (2018). The iron-stress activated RNA 1 (IsaR1) coordinates osmotic acclimation and iron starvation responses in the cyanobacterium Synechocystis sp. PCC 6803. Environ. Microbiol.20, 27572768. doi: 10.1111/1462-2920.14079

  • 37

    ScanlanD. J.OstrowskiM.MazardS.DufresneA.GarczarekL.HessW. R.et al. (2009). Ecological genomics of marine picocyanobacteria. Microbiol. Mol. Biol. Rev.73, 249299. doi: 10.1128/MMBR.00035-08

  • 38

    ShangX.ZhangL. H.ZhangJ. (2007). Prochlorococcus-like populations detected by flow cytometry in the fresh and brackish waters of the Changjiang estuary. J. Mar. Biol. Assoc. U. K.87, 643648. doi: 10.1017/S0025315407055191

  • 39

    SugimotoS.NakayamaJ.FukudaD.SonezakiS.WatanabeM.TosukhowongA.et al. (2003). Effect of heterologous expression of molecular chaperone DnaK from Tetragenococcus halophilus on salinity adaptation of Escherichia coli. J. Biosci. Bioeng.96, 129133. doi: 10.1016/S1389-1723(03)90114-9

  • 40

    ThompsonA. W.HuangK.SaitoM. A.ChisholmS. W. (2011). Transcriptome response of high-and low-light-adapted Prochlorococcus strains to changing iron availability. ISME J.5, 15801594. doi: 10.1038/ismej.2011.49

  • 41

    VoigtK.SharmaC. M.MitschkeJ.LambrechtS. J.VoßB.HessW. R.et al. (2014). Comparative transcriptomics of two environmentally relevant cyanobacteria reveals unexpected transcriptome diversity. ISME J.8, 20562068. doi: 10.1038/ismej.2014.57

  • 42

    WebbE. A.MoffettJ. W.WaterburyJ. B. (2001). Iron stress in open-ocean cyanobacteria (Synechococcus, Trichodesmium, and Crocosphaera spp.): identification of the IdiA protein. Appl. Environ. Microbiol.67, 54445452. doi: 10.1128/AEM.67.12.5444-5452.2001

  • 43

    WestN. J.SchönhuberW. A.FullerN. J.AmannR. I.RippkaR.PostA. F.et al. (2001). Closely related Prochlorococcus genotypes show remarkably different depth distributions in two oceanic regions as revealed by in situ hybridization using 16S rRNA-targeted oligonucleotides. Microbiology147, 17311744. doi: 10.1099/00221287-147-7-1731

  • 44

    ZhangX.ShiZ.YeF.ZengY.HuangX. (2013). Picophytoplankton abundance and distribution in three contrasting periods in the Pearl River estuary, South China. Mar. Freshw. Res.64, 692705. doi: 10.1071/MF12303

  • 45

    ZinserE. R.JohnsonZ. I.CoeA.KaracaE.VenezianoD.ChisholmS. W. (2007). Influence of light and temperature on Prochlorococcus ecotype distributions in the Atlantic Ocean. Limnol. Oceanogr.52, 22052220. doi: 10.4319/lo.2007.52.5.2205

  • 46

    ZwirglmaierK.JardillierL.OstrowskiM.MazardS.GarczarekL.VaulotD.et al. (2008). Global phylogeography of marine Synechococcus and Prochlorococcus reveals a distinct partitioning of lineages among oceanic biomes. Environ. Microbiol.10, 147161. doi: 10.1111/j.1462-2920.2007.01440.x

Summary

Keywords

Prochlorococcus, transcriptome, low salinity acclimation, low salinity stress, RNAseq

Citation

He X, Liu H, Long L, Dong J and Huang S (2022) Acclimation and stress response of Prochlorococcus to low salinity. Front. Microbiol. 13:1038136. doi: 10.3389/fmicb.2022.1038136

Received

06 September 2022

Accepted

27 September 2022

Published

13 October 2022

Volume

13 - 2022

Edited by

Prashant Kumar Singh, Mizoram University, India

Reviewed by

Gulab Chand Arya, Agricultural Research Organization (ARO), Israel; Avinash Singh, New York University, United States; Deepanker Yadav, Guru Ghasidas Vishwavidyalaya, India

Updates

Copyright

*Correspondence: Sijun Huang,

This article was submitted to Microbial Physiology and Metabolism, a section of the journal Frontiers in Microbiology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics