Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2020 Sep 16;15(9):e0238405.
doi: 10.1371/journal.pone.0238405. eCollection 2020.

Persistent El Niño driven shifts in marine cyanobacteria populations

Affiliations

Persistent El Niño driven shifts in marine cyanobacteria populations

Alyse A Larkin et al. PLoS One. .

Abstract

In the California Current Ecosystem, El Niño acts as a natural phenomenon that is partially representative of climate change impacts on marine bacteria at timescales relevant to microbial communities. Between 2014-2016, the North Pacific warm anomaly (a.k.a., the "blob") and an El Niño event resulted in prolonged ocean warming in the Southern California Bight (SCB). To determine whether this "marine heatwave" resulted in shifts in microbial populations, we sequenced the rpoC1 gene from the biogeochemically important picocyanobacteria Prochlorococcus and Synechococcus at 434 time points from 2009-2018 in the MICRO time series at Newport Beach, CA. Across the time series, we observed an increase in the abundance of Prochlorococcus relative to Synechococcus as well as elevated frequencies of ecotypes commonly associated with low-nutrient and high-temperature conditions. The relationships between environmental and ecotype trends appeared to operate on differing temporal scales. In contrast to ecotype trends, most microdiverse populations were static and possibly reflect local habitat conditions. The only exceptions were microdiversity from Prochlorococcous HLI and Synechococcus Clade II that shifted in response to the 2015 El Niño event. Overall, Prochlorococcus and Synechococcus populations did not return to their pre-heatwave composition by the end of this study. This research demonstrates that extended warming in the SCB can result in persistent changes in key microbial populations.

PubMed Disclaimer

Conflict of interest statement

The authors have declared that no competing interests exist.

Figures

Fig 1
Fig 1. Seasonal and interannual changes in environmental conditions and ecotype frequencies.
(A) Temperature (°C), nitrate and phosphate concentrations (μM), and the log10-transformed ratio of Prochlorococcus sequence counts relative to Synechococcus sequence counts in the MICRO time series. (B) Interpolated relative abundance of rpoC1 reads taxonomically assigned to the six most abundant Prochlorococcus (Pro.) and Synechococcus (Syn.) ecotypes. The sampling dates are marked along the top of the time series. (C) Linear model regression coefficients for environmental variables and relative ecotype abundance as a function month and year (see Methods). A trend of 0 is marked with horizontal dashed lines.
Fig 2
Fig 2. Pearson correlation coefficients (r) between relative abundance trends and environmental trends.
The significance of each correlation is marked as p-value: * < 0.05, ** < 0.01, *** <0.001. Bars without a mark had a non-significant correlation.
Fig 3
Fig 3. Microdiversity clusters within the most abundant Prochlorococcus and Synechococcus ecotypes were temporally and statistically stable.
Dendrograms represent the average-linkage hierarchical clustering of single nucleotide polymorphisms (SNP) differences between a rpoC1 reference sequence and the majority consensus sequence at each sampling time point. Each row of the dendrogram represents the consensus rpoC1 sequence at a time point, where only SNPs (in comparison to the reference sequence) are depicted. Only time points where the ecotype was present are shown. The reference rpoC1 sequence is depicted at the top of each dendrogram. Reference genomes used to compare rpoC1 sequences included: HLI–MIT9515; LLI–NATL2A; HLII–MIT9301; Clade I–WH8016; Clade IV–CC9902; and, Clade II–CC9605. Within-cluster sum of squares was used to determine the optimum number of SNP-based clusters (i.e., haplotypes). Cluster labels (red) indicate bootstrapped Jaccard similarity values for the most stable clusters (values < 0.5 are unstable, 0.6–0.75 weak stability, 0.75–0.85 stable, and > 0.85 highly stable). Stable clusters are alternately shaded in grey for emphasis.
Fig 4
Fig 4. Prochlorococcus HLI haplotypes (Pro.HLI.1 and Pro.HLI.2) and Synechococcus II haplotypes (Syn.II.1 and Syn.II.2) showed significant shifts in composition across the time series.
(A) Interpolated relative abundance of the haplotypes. Samples with less than 10 reads were removed from the analysis. (B) Linear regression coefficients for monthly and yearly trends in the Pro.HLI.2 and Syn.II.2 haplotypes. (C) Pearson’s correlation coefficients (ρ) comparing haplotype to environmental monthly and yearly trends. The significance of each correlation is marked as p-value: * < 0.05. Bars without a mark had a non-significant correlation.
Fig 5
Fig 5. Comparison of Prochlorococcus HLII and Synechococcus Clade II rpoC1 highly conserved SNPs detected at the Newport MICRO times series and on Pacific Ocean cruise NH1418 (September-October 2014; 19.00°N, -158.00°W to -3.00°N, -149.67°W; 5m depth).
Newport and NH1418 samples largely clustered separately but had some overlap in SNPs. Samples with less than 10 sequences were removed from the analysis. Clade labels (red) indicate bootstrapped Jaccard similarity values for the most stable clusters with more than two samples (values < 0.5 are unstable, 0.6–0.75 weak stability, 0.75–0.85 stable, and > 0.85 highly stable).
Fig 6
Fig 6. Conceptual diagram of picocyanobacterial community shifts at the MICRO time series from 2009–2018.
Synechococcus cells are outlined in purple and Prochlorococcus cells are outlined in green. Blue cell shading indicates cold temperature-adapted ecotypes and yellow cell shading indicates warm temperature-adapted ecotypes. The blue line indicates the nitrate trend (N) and the red line indicates the temperature trend (T).

Similar articles

Cited by

References

    1. Polovina JJ, Dunne JP, Woodworth PA, Howell EA. Projected expansion of the subtropical biome and contraction of the temperate and equatorial upwelling biomes in the North Pacific under global warming. ICES Journal of Marine Science. 2011;68(6):986–95.
    1. Moore JK, Lindsay K, Doney SC, Long MC, Misumi K. Marine Ecosystem Dynamics and Biogeochemical Cycling in the Community Earth System Model CESM1(BGC): Comparison of the 1990s with the 2090s under the RCP4.5 and RCP8.5 Scenarios. Journal of Climate. 2013;26(23):9291–312.
    1. Fuhrman JA, Hewson I, Schwalbach MS, Steele JA, Brown MV, Naeem S. Annually reoccurring bacterial communities are predictable from ocean conditions. Proceedings of the National Academy of Sciences of the United States of America. 2006;103(35):13104–9. 10.1073/pnas.0602399103 - DOI - PMC - PubMed
    1. Gilbert JA, Steele JA, Caporaso JG, Steinbrueck L, Reeder J, Temperton B, et al. Defining seasonal marine microbial community dynamics. ISME Journal. 2012;6(2):298–308. 10.1038/ismej.2011.107 - DOI - PMC - PubMed
    1. Martiny AC, Talarmin A, Mouginot C, Lee JA, Huang JS, Gellene AG, et al. Biogeochemical interactions control a temporal succession in the elemental composition of marine communities. Limnology and Oceanography. 2016;61(2):531–42.

Publication types

MeSH terms

Grants and funding

Financial support for this work was provided by the UCI Undergraduate Research Opportunities Program (to AJF and AF), NSF Graduate Research Fellowship Program and UCI Chancellor’s Club Fellowship (to ARM), and NSF Biological Oceanography (OCE-1848576 and OCE-1948842 to ACM). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.