Next Article in Journal
New Trends in Meat Packaging
Previous Article in Journal
Optimizing a Suspension Culture Method with a Decreased Cost to Detect Enteroviruses in Water to Increase Surveillance Access
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Divergent Microbiota Dynamics along the Coastal Marine Ecosystem of Puerto Rico

by
Clifford Jaylen Louime
1,*,
Frances Vazquez-Sanchez
2,
Dieunel Derilus
1 and
Filipa Godoy-Vitorino
2,*
1
Department of Environmental Sciences, University of Puerto Rico, San Juan, PR 00931, USA
2
Department of Microbiology & Medical Zoology, School of Medicine, University of Puerto Rico, San Juan, PR 00936, USA
*
Authors to whom correspondence should be addressed.
Submission received: 24 June 2020 / Accepted: 30 August 2020 / Published: 21 December 2020

Abstract

:
Understanding the different factors shaping the spatial and temporal distribution of marine microorganisms is fundamental in predicting their responses to future environmental disturbances. There has been, however, little effort to characterize the microbial diversity including the microbiome dynamics among regions in the Caribbean Sea. Toward this end, this study was designed to gain some critical insights into microbial diversity within the coastal marine ecosystem off the coast of Puerto Rico. Using Illumina MiSeq, the V4 region of the 16S rRNA gene was sequenced with the goal of characterizing the microbial diversity representative of different coastal sites around the island of Puerto Rico. This study provided valuable insights in terms of the local bacterial taxonomic abundance, α and β diversity, and the environmental factors shaping microbial community composition and structure. The most dominant phyla across all 11 sampling sites were the Proteobacteria, Bacteroidetes, and Planctomycetes, while the least dominant taxonomic groups were the NKB19, Tenericutes, OP3, Lentisphaerae, and SAR406. The geographical area (Caribbean and Atlantic seas) and salinity gradients were the main drivers shaping the marine microbial community around the island. Despite stable physical and chemical features of the different sites, a highly dynamic microbiome was observed. This highlights Caribbean waters as one of the richest marine sources for a microbial biodiversity hotspot. The data presented here provide a basis for further temporal evaluations aiming at deciphering microbial taxonomic diversity around the island, while determining how microbes adapt to changes in the climate.

1. Introduction

A key objective in microbial ecology is to understand the main factors affecting the spatial and temporal distributions of microbial taxa in the environment. This is fundamental in predicting the ecosystem response to upcoming environmental changes. Marine microbiota is highly dilute and differs according to temperature, geography, and other environmental factors [1], changing in composition with symbiotic association to marine animal systems, such as corals [2,3,4,5,6] or sponges [7,8,9] or even anthropogenic impact [10].The communities conformed by these microorganisms, respond to changes in the environment, such as shifts in the temperature as a result of climate change and nutrient cycling [11,12]. Indeed, by controlling the cycling of elements essential for life on earth, microbes are the backbone of all our planet’s ecosystems. From nitrogen/carbon fixation to sulfur metabolism, the role of microbes on Earth is so immensely significant, that life on Earth as we know it, would not exist, if all microbes were to disappear [13].
Next Generation Sequencing (NGS) projects, such as the Global Ocean Sampling (GOS) expedition in the Sargasso Sea—one of the largest metagenomic studies so far, comes with twice the data size of the human genome, found 1.2 million genes, and inferred the presence of at least 1800 bacterial species [14,15,16]. The Tara Expeditions have also revealed that single-celled eukaryotes (protists) are critical players in global biogeochemical cycling of nutrients and energy in the oceans [17]. This confirms the role of microbes in the global biogeochemical cycles of carbon and nitrogen as already discussed and the underlying evidence that genomic and metabolic plasticity are the basis of microbial adaptation in marine ecosystems [15]. However, the marine environment, with its dispersion due to water currents has also an inherent complexity, in contrast to the static terrestrial environment. This limits the understanding of microbiome dynamics in oceans, and therefore, highlights the importance of localized studies that define core microbiota analyses in different geographical locations so that in the coming years, a better picture of the global microbial composition and structure is available. Furthermore, understanding the general mechanisms that shape the spatial distribution of bacterial taxa in specific marine microclimates, is a prerequisite to predict these ecosystems responses to eventual and future environmental changes.
The tropical island of Puerto Rico is part of a volcanic island platform, in the northernmost Atlantic side marked by higher wave energy and terrigenous sediments, while the Caribbean side is dominated by carbonate platforms with coral reefs and dominantly carbonate sediments [18]. These inherent coastal differences are also expected to result in different microbial community structures. Understanding the different environmental parameters that shape the distribution of costal microbiomes is crucial to predict and anticipate their response to the ongoing and accelerated environmental changes. In the present study, we aimed to compare the microbial community structure across coastal sites, across littoral Puerto Rico, and samples were grouped in the areas of the Caribbean and Atlantic seas, studied according to the dissolved oxygen, pH, temperature, and salinity gradients. The results of this work pointed out the relative importance of considering some key environmental variables in understanding tropical marine ecosystem biomarkers.

2. Materials and Methods

2.1. Sample Collection and Processing

The collection and processing of seawater samples was performed according to a previous procedure [19]. During the month of May 2016, right at the transition of the dry and the rainy season on the island, a total of 11 samples, six seawater and five mostly from sediment, were taken at different sites across the geography of Puerto Rico including, the north, south east, and western coasts (Figure 1 and Table 1). For the meta-taxonomic analysis, replicate samples were drawn from the water at a depth between 3 and 10 m, using sterile 1-L polypropylene bottles and kept on ice until processing in the laboratory.
During the sampling expedition, various in situ physical measurements were performed on the water samples using the CTD probe (Fisher Scientific—Meter #6534251, Pittsburgh, PA, USA) equipped with water temperature, pH, salinity, and dissolved oxygen sensor. For genomic characterization, 1 L of each aggregated sample (water, sediment, and biofilms) was spun down at a moderate speed of 5000 rpm in order to avoid shearing of the DNA and the resulting pellet was stored frozen at −80 °C for later processing. Metagenomic DNA was isolated with the FastDNA™ Soil Kit (MPBiomedicals, Inc., Irvine, CA, USA), according to the manufacturer’s instructions. The quality and quantity of the DNA was first verified with the Nanodrop 1000 then run on a 1% DNA gel. The resulting purified DNA samples were stored at −80 °C for PCR analysis.

2.2. 16S rRNA Gene Amplicons and Sequence Analysis

The V4 hypervariable region of the 16S ribosomal RNA (~291 bp length) was amplified by PCR using the universal bacterial and archaeal primers: 515F (5′ GTGCCAGCMGCCGCGGTAA 3′) and 806R (5′ GGACTACHVGGGTWTCTAAT 3′) as used in the Earth Microbiome Project [20]. Amplicons of ~300 bp were multiplex paired-end sequenced in the Illumina Miseq platform. The resulting 16S rRNA gene sequences per sample were deposited in NCBI’s Bioproject PRJNA392167, ID: 392167 (NCBI 2018).
The software package QIIME v.1.9 (Quantitative Insights into Microbial Ecology) was used to analyze the generated sequences [21,22]. Processing was done as detailed [23]. Sequenced reads were demultiplexed and processed with a default Phred score of 20 (99% accuracy of the base calls), size (>200 bp), and chimera filtering with the Usearch Hierarchical Clustering Method [24].
Sequences were binned (97% homology) via the open reference-based method, implementing UCLUST [24] as the default OTU clustering tool used in QIIME’s workflows. Taxonomic classification was performed using the Ribosomal Database Project RDP v 2.2 [25] and the OTU table was built using the RDP previously assigned taxonomy and OTU map, which contains the OTU assignment and all of its sequences. The generated OTU table served as the basis for alpha diversity and beta diversity across sites. The same OTU table was used to statistically estimate the average abundance-based coverage of the samples. Finally, the taxa unique to each sample or shared among samples were analyzed and separated statistically using the R software package [26] to perform the regression analysis based on the metadata collected. Beta diversity analyses using Non-Metric Multidimensional Scaling (nMDS) ordinations of Bray Curtis were done using distance metrics computed from the rarefied OTU table and the metadata. We used nMDS ordination, achieved by the metaMDS wrapper function as well as boxplots, which were generated using the vegan package [27]. The ordination was applied such that the data were scaled down to two dimensions. The heatmap showing the taxa that significantly differed in abundance between depths (p < 0.05) was built using the heatmap.3 function in R [28].

3. Results and Discussion

Puerto Rico is a small island of 100 by 35 miles. Therefore, by collecting samples along all four sides of the quasi rectangle-shaped island (Figure 1), we were able to get a representative coverage of the bacterioplankton community diversity around the island. The 10 collected samples resulted in an average of 108,095 ± 51,855 raw reads. Due to the variability inherent to water sampling, and the fact that current NGS technologies such as Illumina, have higher error rates than traditional Sanger sequencing (Ratan 2013), we decided to perform a conservative data analyses using only OTUs corresponding to the 100% core microbiome at different categories. The number of sequences used in the analyses was an average of 29,324 ± 10,149 reads (Table 1). These sequences yielded ~2309 ± 221 OTUs.
A total of 32 phyla were present in the 11 sampling sites including one unclassified bacterial and two archaeal phyla. The most dominant phyla across all 11 samples were the Proteobacteria, Bacteroidetes, and Planctomycetes, while the least dominant taxonomic groups were the NKB19, Tenericutes, OP3, Lentisphaerae, and SAR406 (Table S1). A total of 293 Family and/or Genus level-OTUs were found, and the top dominant taxa were the Piscirickettsiaceae, Rhodobacteraceae (Rhodovulum), and Vibrio (Table S2), which correspond to the core microbiome since they were overrepresented in almost all of the samples. The least dominant taxa at the family and genus levels were the Chthoniobacteraceae, Pelagibacteraceae (SAR11), Roseovarius, or Reinekea (Table S2). The low detection of ubiquitous Pelagibacteraceae (SAR11) (<0.5%) is notable and was not expected. This is probably associated to the PCR bias of the V4 primer used in this survey. In the GOS dataset, 30.9% of the total 16S rDNA reads retrieved from the shotgun metagenomic dataset belonged to the SAR1 [29], whereas only 0.5–1.45% of the 16S rDNA amplicons reads in the BR (Benham Rise, Philippines, Western Pacific Ocean) dataset belonged to this group [30]. The shotgun metagenomics of the same marine samples revealed that members of SAR11 represent more than 28% of the total prokaryote abundance, meanwhile 16S rDNA amplicons revealed that these ubiquitous group represented less than 2.5% of the total prokaryotic community [31]. The primer set 515F-806R [32] used in this study, has been reported for low bending efficiency of the reverse primer for the SAR11 cluster [33]. Here, we can assume that the underrepresentation of the SAR11 taxa in these pelagic samples is likely due to the inherent PCR primer bias.
We found no significant differences in the environmental parameters between sampling sites (dissolved oxygen, pH, temperature), nor between the four cardinal areas, as well as between the sediment and water (p > 0.05). However, significant differences in the relative abundance of certain taxa were found between the Atlantic and Caribbean sides of the island. Indeed, we found a clear separation on the beta diversity plot, between the Caribbean and the Atlantic areas (p-value = 0.004) (Figure 2A). Interestingly, communities were also separated according to salinity (Figure 2B) (p-value = 0.004). Alpha diversity (species diversity) based on the Shannon index, was found to be significantly higher in the Caribbean sampling sites compared to the Atlantic Areas (Figure 3A; p-value = 0.007). A significant value was taking into account that the Caribbean area is only composed of four samples, while the Atlantic area is composed of six samples. Similarly, diversity was found to be significantly higher in lower salinity samples 34–36 ppt, compared to those areas with a higher salinity 37–39 ppt (p-value = 0.012—Figure 3B). The 34–36 ppt range is made up of four samples and the 37–39 ppt range is composed of six samples. This result suggests that, salinity has a very high impact on microbial community along the littoral coast of Puerto Rico. This result is consistent with previous studies that showed that the species richness and evenness declined as salinity increased along an estuarine salinity gradient [34].
The main differences between the microbial communities of the Atlantic and Caribbean include a higher abundance of Verrucomicrobia and Bacteroidetes in the Caribbean, while Preoteobacteria, Actinobacteria, and Crenarchaeota were more abundant in the Atlantic sites (Figure S1). At the family and genus levels, we found the Atlantic areas to be significantly more abundant in Vibrio, Rhodobacteraceae, Rhodovulum, and Pisciricketsiaceae, while Chromatiales, Desulfococcus, Desulfobacteraceae (sulfate reducing bacteria), Flameovirgaceae or Alteromonadales, Pirellulaceae, and other Planctomycetes were more abundant in the Caribbean (Figure 4A,B and Figure 5). Sulfate reducing bacterial consortia such as Desulfococcus and Desulfobulbus have also been found to be abundant in the Caribbean sites (Figure 5), and were found mostly associated with coral systems [2,35,36]. The increase in sulfate-reducing microorganisms in the cycling of elements in Caribbean sites, is thermodynamically favored in relation to, for example, methanogenesis, and could be related to a decrease in the methane flux in the Caribbean sites. An increase in halophilic taxa such as genus Thiohalorhabdus is also shown in the Caribbean sites.
Some of the bacterial taxa which were found in our survey include a diversity of aerobic generalistic hydrocarbon degraders, which can utilize hydrocarbons and non-hydrocarbon substrates as a source of carbon and energy such as Pisciricketsiaceae (Figure 4B and Figure 5). These findings supported previous ones on microbial taxa found in the Deepwater Horizon oil spill [37]. Planctomycetes, abundant in the Caribbean, are reported worldwide in marine water samples and are especially known to be associated with microalgae [38]. This association with marine bacterioplankton and algae is of considerable interest as the bacteria may play an active role in controlling harmful algal blooms in oceans. Moreover, the increase of this taxa in the Caribbean, which maybe associated to ammonia is oxidized by nitrate to nitrogen gas, yielding energy [39].
Therefore, the Caribbean is a good model for these expected changes due to lower salinity levels and elevated silica [40]. In addition, understanding the microbial ecology of the different marine ecosystems is essential for our ability to assess the importance of biogeochemical cycles-climate feedbacks in the climate change era.

4. Conclusions

Microbes regulate the health and dynamics of most ecosystems. Therefore, it is vital to continuously monitor not only terrestrial but also aquatic ecosystems in order to develop quick responses to disturbances. Our 16S rRNA survey of samples collected across the island of Puerto Rico revealed that the geographical location and salinity were the main drivers of microbial structure in pelagic water. Our findings suggest that microbial community diversity across the island of Puerto Rico changes depending on the geographic location of either the Atlantic and the Caribbean with an increase in sulfate reducers and ammonia oxidizers. Remarkably, we found that water salinity is the other main driver of the bacterial assemblages in the pelagic water across this Caribbean Island. Sea surface salinity is an essential variable of climate change, which plays a fundamental role in ocean circulation, water circle, and climate. This result suggests that marine microbial communities are very sensitive to very small changes in sea surface salinity. It is important to highlight the limitation of this survey and some other microbiome analysis, as these rely on PCR amplification of 16S rRNA and some key taxa may fail to be detected due to the inherent primer bias and sensitivity. Furthermore, a high diversity of the metabolic capacities of the microbial communities could not be assessed using the 16S rRNA approach. Further studies with more extensive sampling around the coastal water of this Caribbean region with the use of shotgun metagenomics will provide more insight regarding the response of the microbial communities and functions to environmental changes in this local marine ecosystem.
As demonstrated before, the role played by anthropogenic and natural impacts can be significant as well [41]. We found the Caribbean Sea samples to be an interesting model to study climate change due to the relative abundance of sulfate reducing bacteria, previously associated to the demise of coral ecosystems [42]. Nonetheless, the few numbers of samples suggest that the microbiota around the island of Puerto Rico warrant further time-scaled investigations. The island is climatically characterized by two different seasons, the rainy season that goes from July through November and the dry season from December until June. Heavy precipitations or drought are known to alter the salinity or pH of water bodies, thereby significantly changing the local microbial community [43]. Understanding how the natural water microbiota assemblages change in the Caribbean, may shed light on how global phytoplankton will cope with a changing climate. Although they have been reported to contribute very little to global warming, island nations seem to be especially vulnerable to the impacts of climate change [44].
One of the most alarming effects is the sea level rise directly caused by the expansion of seawater as it warms. These entail among others, marine species dying out or changing their metabolism and range. Such loss in biodiversity is proven to be detrimental to the biogeochemical cycling of nutrients including atmospheric gas production and carbon sequestration by microbes in the oceans. The balance of all these cycles and compounds controls the dynamics of all ocean biomes. Therefore, it is paramount to continue looking at microorganisms in order to get a better understanding of global change in the quest of finding solutions to unprecedented and undesirable changes.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2036-7481/11/2/9/s1.

Author Contributions

All authors contributed equally to the manuscript: Conceptualization, C.J.L.; methodology, C.J.L., D.D., and F.G.-V.; funding, C.J.L. and F.G.-V.; software (analysis, validation), F.G.-V., F.V.-S., and D.D.; writing—original draft preparation, C.J.L. and F.G.-V.; writing—review and editing, all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the US Department of Education through their HSI-STEM Program under grant number DE-AD13-076ET34652. This work was also partially supported by an INBRE (NIH/NIGMS- award number P20 GM103475). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Acknowledgments

This research originated from the “Center for Renewable Energy and Sustainability (CRES-Center)” at the University of Puerto Rico Rio Piedras.

Conflicts of Interest

The authors certify that they have no affiliations with or involvement in any organization or entity with any financial interest or non-financial interest in the subject matter or materials discussed in this manuscript.

References

  1. Sunagawa, S.; Coelho, L.P.; Chaffron, S.; Kultima, J.R.; Labadie, K.; Salazar, G.; Djahanschiri, B.; Zeller, G.; Mende, D.R.; Alberti, A.; et al. Ocean plankton. Structure and function of the global ocean microbiome. Science 2015, 348, 1261359. [Google Scholar] [CrossRef] [Green Version]
  2. Godoy-Vitorino, F.; Ruiz-Diaz, C.P.; Rivera-Seda, A.; Ramirez-Lugo, J.S.; Toledo-Hernandez, C. The microbial biosphere of the coral Acropora cervicornis in Northeastern Puerto Rico. PeerJ 2017, 5, e3717. [Google Scholar] [CrossRef]
  3. Lawler, S.N.; Kellogg, C.A.; France, S.C.; Clostio, R.W.; Brooke, S.D.; Ross, S.W. Coral-Associated Bacterial Diversity Is Conserved across Two Deep-Sea Anthothela Species. Front. Microbiol. 2016, 7, 458. [Google Scholar] [CrossRef]
  4. Sun, W.; Anbuchezhian, R.; Li, Z. (Eds.) Association of Coral-Microbes, and the Ecological Roles of Microbial Symbionts in Corals. In The Cnidaria, Past, Present and Future; Springer: Cham, Switzerland, 2016; pp. 347–357. [Google Scholar]
  5. Ainsworth, T.D.; Fordyce, A.J.; Camp, E.F. The Other Microeukaryotes of the Coral Reef Microbiome. Trends Microbiol. 2017, 25, 980–991. [Google Scholar] [CrossRef]
  6. Hernandez-Agreda, A.; Gates, R.D.; Ainsworth, T.D. Defining the Core Microbiome in Corals’ Microbial Soup. Trends Microbiol. 2017, 25, 125–140. [Google Scholar] [CrossRef]
  7. Fan, L.; Liu, M.; Simister, R.; Webster, N.S.; Thomas, T. Marine microbial symbiosis heats up: The phylogenetic and functional response of a sponge holobiont to thermal stress. ISME J. 2013, 7, 991–1002. [Google Scholar] [CrossRef]
  8. Karlinska-Batres, K.; Worheide, G. Phylogenetic diversity and community structure of the symbionts associated with the coralline sponge Astrosclera willeyana of the Great Barrier Reef. Microb. Ecol. 2013, 65, 740–752. [Google Scholar] [CrossRef]
  9. Steinert, G.; Taylor, M.W.; Deines, P.; Simister, R.L.; de Voogd, N.J.; Hoggard, M.; Schupp, P.J. In four shallow and mesophotic tropical reef sponges from Guam the microbial community largely depends on host identity. PeerJ 2016, 4, e1936. [Google Scholar] [CrossRef]
  10. Gilbert, J.A.; Steele, J.A.; Caporaso, J.G.; Steinbruck, L.; Reeder, J.; Temperton, B.; Huse, S.; McHardy, A.C.; Knight, R.; Joint, I.; et al. Defining seasonal marine microbial community dynamics. ISME J. 2012, 6, 298–308. [Google Scholar] [CrossRef] [Green Version]
  11. Harley, C.D.; Anderson, K.M.; Demes, K.W.; Jorve, J.P.; Kordas, R.L.; Coyle, T.A.; Graham, M.H. Effects of Climate Change on Global Seaweed Communities. J. Phycol. 2012, 48, 1064–1078. [Google Scholar] [CrossRef]
  12. Harley, C.D.; Randall Hughes, A.; Hultgren, K.M.; Miner, B.G.; Sorte, C.J.; Thornber, C.S.; Rodriguez, L.F.; Tomanek, L.; Williams, S.L. The impacts of climate change in coastal marine systems. Ecol. Lett. 2006, 9, 228–241. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Tortell, P.D.; Maldonado, M.T.; Granger, J.; Price, N.M. Marine bacteria and biogeochemical cycling of iron in the oceans. FEMS Microbiol. Ecol. 2006, 29. [Google Scholar] [CrossRef]
  14. Piganeau, G.; Desdevises, Y.; Derelle, E.; Moreau, H. Picoeukaryotic sequences in the Sargasso sea metagenome. Genome Biol. 2008, 9, R5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Yooseph, S.; Nealson, K.H.; Rusch, D.B.; McCrow, J.P.; Dupont, C.L.; Kim, M.; Johnson, J.; Montgomery, R.; Ferriera, S.; Beeson, K.; et al. Genomic and functional adaptation in surface ocean planktonic prokaryotes. Nature 2010, 468, 60–66. [Google Scholar] [CrossRef] [Green Version]
  16. Williamson, S.J.; Rusch, D.B.; Yooseph, S.; Halpern, A.L.; Heidelberg, K.B.; Glass, J.I.; Andrews-Pfannkoch, C.; Fadrosh, D.; Miller, C.S.; Sutton, G.; et al. The Sorcerer II Global Ocean Sampling Expedition: Metagenomic characterization of viruses within aquatic microbial samples. PLoS ONE 2008, 3, e1456. [Google Scholar] [CrossRef]
  17. Seeleuthner, Y.; Mondy, S.; Lombard, V.; Carradec, Q.; Pelletier, E.; Wessner, M.; Leconte, J.; Mangot, J.F.; Poulain, J.; Labadie, K.; et al. Single-cell genomics of multiple uncultured stramenopiles reveals underestimated functional diversity across oceans. Nat. Commun. 2018, 9, 310. [Google Scholar] [CrossRef]
  18. Hubbard, D.; Burke, R.; Gill, P.; Ramirez, W.; Sherman, C. Coral-reef Geology: Puerto Rico and the US Virgin Islands. In Coral Reefs of the USA; Springer Science: New York, NY, USA, 1970. [Google Scholar]
  19. Suh, S.S.; Park, M.; Hwang, J.; Kil, E.J.; Jung, S.W.; Lee, S.; Lee, T.K. Seasonal Dynamics of Marine Microbial Community in the South Sea of Korea. PLoS ONE 2015, 10, e0131633. [Google Scholar] [CrossRef] [Green Version]
  20. Gilbert, J.A.; Meyer, F.; Jansson, J.; Gordon, J.; Pace, N.; Tiedje, J.; Ley, R.; Fierer, N.; Field, D.; Kyrpides, N.; et al. The Earth Microbiome Project: Meeting report of the “1 EMP meeting on sample selection and acquisition” at Argonne National Laboratory October 6 2010. Stand. Genom. Sci. 2010, 3, 249–253. [Google Scholar] [CrossRef] [Green Version]
  21. Kuczynski, J.; Stombaugh, J.; Walters, W.A.; Gonzalez, A.; Caporaso, J.G.; Knight, R. Using QIIME to analyze 16S rRNA gene sequences from microbial communities. Curr. Protoc. Bioinform. 2011, 21, 1E.5.1–1E.5.20. [Google Scholar] [CrossRef] [Green Version]
  22. Caporaso, J.G.; Lauber, C.L.; Walters, W.A.; Berg-Lyons, D.; Huntley, J.; Fierer, N.; Owens, S.M.; Betley, J.; Fraser, L.; Bauer, M.; et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012, 6, 1621–1624. [Google Scholar] [CrossRef] [Green Version]
  23. Derilus, D.; Forestil, A.; Fortuné, J.; Polyanska, O.; Louime, C.; Gervais, G.; Massey, S.E. Functional Metagenomics Characterization of an Anaerobic Saltwater Bioreactor. J. Renew. Energy 2019, 1–15. [Google Scholar] [CrossRef]
  24. Edgar, R.C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010, 26, 2460–2461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Cole, J.R.; Chai, B.; Farris, R.J.; Wang, Q.; Kulam-Syed-Mohideen, A.S.; McGarrell, D.M.; Bandela, A.M.; Cardenas, E.; Garrity, G.M.; Tiedje, J.M. The ribosomal database project (RDP-II): Introducing myRDP space and quality controlled public data. Nucleic Acids Res. 2007, 35, D169–D172. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. R Development Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2008; Available online: http://www.R-project.org (accessed on 31 May 2020).
  27. Oksanen, J.; Blanchet, F.G.; Kindt, R.; Legendre, P.; O’hara, R.; Simpson, G.L.; Solymos, P.; Stevens, M.H.H.; Wagner, H. Vegan: Community Ecology Package. R Package Version 2.4-4; R Development Core Team: Vienna, Austria, 2010. [Google Scholar]
  28. Zhao, S.; Guo, Y.; Sheng, Q.; Shyr, Y. Advanced heat map and clustering analysis using heatmap3. Biomed. Res. Int. 2014, 986048. [Google Scholar] [CrossRef] [Green Version]
  29. Biers, E.J.; Sun, S.; Howard, E.C. Prokaryotic genomes and diversity in surface ocean waters: Interrogating the global ocean sampling metagenome. Appl. Environ. Microbiol. 2009, 75, 2221–2229. [Google Scholar] [CrossRef] [Green Version]
  30. Gajigan, A.P.; Yniguez, A.T.; Villanoy, C.L.; San Diego-McGlone, M.L.; Jacinto, G.S.; Conaco, C. Diversity and community structure of marine microbes around the Benham Rise underwater plateau, northeastern Philippines. PeerJ 2018, 6, e4781. [Google Scholar] [CrossRef] [Green Version]
  31. Wilson, B.; Müller, O.; Nordmann, E.-L.; Seuthe, L.; Bratbak, G.; Øvreås, L. Changes in Marine Prokaryote Composition with Season and Depth Over an Arctic Polar Year. Front. Mar. Sci. 2017, 4. [Google Scholar] [CrossRef] [Green Version]
  32. Caporaso, J.G.; Lauber, C.L.; Walters, W.A.; Berg-Lyons, D.; Lozupone, C.A.; Turnbaugh, P.J.; Fierer, N.; Knight, R. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl. Acad. Sci. USA 2011, 108 (Suppl. 1), 4516–4522. [Google Scholar] [CrossRef] [Green Version]
  33. Apprill, A.; McNally, S.; Parsons, R.; Weber, L. Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton. Aquat. Microb. Ecol. 2015, 75, 129–137. [Google Scholar] [CrossRef] [Green Version]
  34. Campbell, B.J.; Kirchman, D.L. Bacterial diversity, community structure and potential growth rates along an estuarine salinity gradient. ISME J. 2013, 7, 210–220. [Google Scholar] [CrossRef] [Green Version]
  35. Sekar, R.; Mills, D.K.; Remily, E.R.; Voss, J.D.; Richardson, L.L. Microbial communities in the surface mucopolysaccharide layer and the black band microbial mat of black band-diseased Siderastrea siderea. Appl. Environ. Microbiol. 2006, 72, 5963–5973. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Godoy-Vitorino, F.; Toledo-Hernandez, C. Reef-Building Corals as a Tool for Climate Change Research in the Genomics Era. Results Probl. Cell Differ. 2018, 65, 529–546. [Google Scholar] [CrossRef] [PubMed]
  37. Berry, D.; Gutierrez, T. Evaluating the Detection of Hydrocarbon-Degrading Bacteria in 16S rRNA Gene Sequencing Surveys. Front. Microbiol. 2017, 8, 896. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Bondoso, J.; Godoy-Vitorino, F.; Balague, V.; Gasol, J.M.; Harder, J.; Lage, O.M. Epiphytic Planctomycetes communities associated with three main groups of macroalgae. FEMS Microbiol. Ecol. 2017, 93. [Google Scholar] [CrossRef] [Green Version]
  39. Liu, C.; Yu, D.; Wang, Y.; Chen, G.; Tang, P.; Huang, S. A novel control strategy for the partial nitrification and anammox process (PN/A) of immobilized particles: Using salinity as a factor. Bioresour. Technol. 2020, 302, 122864. [Google Scholar] [CrossRef] [PubMed]
  40. Corredor, J.; Morell, J. Seasonal variation of physical and biogeochemical features in eastern Caribbean surface water. J. Geophys. Res. 2001, 106, 4517–4525. [Google Scholar] [CrossRef]
  41. Kirchman, D.L.; Cottrell, M.T.; Lovejoy, C. The structure of bacterial communities in the western Arctic Ocean as revealed by pyrosequencing of 16S rRNA genes. Environ. Microbiol. 2010, 12, 1132–1143. [Google Scholar] [CrossRef] [PubMed]
  42. Rudiak-Gould, P. Promiscuous corroboration and climate change translation: A case study from the Marshall Islands. Glob. Environ. Chang. 2012, 22, 46–54. [Google Scholar] [CrossRef]
  43. Gifford, S.; Sharma, S.; Booth, M.; Moran, M. Expression patterns reveal niche diversification in a marine microbial assemblage. ISME J. 2013, 7, 281–298. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Patel, S. Climate science: A sinking feeling. Nature 2006, 440, 734–736. [Google Scholar] [CrossRef]
Figure 1. Diagram of the 11 sampling locations across the geography of Puerto Rico.
Figure 1. Diagram of the 11 sampling locations across the geography of Puerto Rico.
Microbiolres 11 00009 g001
Figure 2. Non-metric multidimensional scaling (NMS) ordination showing the level of dissimilarity of microbial community composition between the geographical area (A) and salinity (B), as explained by NMS axes 1 (NMDS1) and 2 (NMDS2).
Figure 2. Non-metric multidimensional scaling (NMS) ordination showing the level of dissimilarity of microbial community composition between the geographical area (A) and salinity (B), as explained by NMS axes 1 (NMDS1) and 2 (NMDS2).
Microbiolres 11 00009 g002
Figure 3. Box plots showing the alpha diversity measures (Shannon index) along the geographical area (A) and salinity (B).
Figure 3. Box plots showing the alpha diversity measures (Shannon index) along the geographical area (A) and salinity (B).
Microbiolres 11 00009 g003
Figure 4. Taxonomic assignment according to the geographical areas. Bar plot showing the relative abundance of different taxa (mostly at the genus and family levels) in the 11 samples (A). Heatmap showing the most significant OTUs by area (p-value < 0.05) (B).
Figure 4. Taxonomic assignment according to the geographical areas. Bar plot showing the relative abundance of different taxa (mostly at the genus and family levels) in the 11 samples (A). Heatmap showing the most significant OTUs by area (p-value < 0.05) (B).
Microbiolres 11 00009 g004
Figure 5. Boxplots of significantly different genus-level OTUs in the Atlantic and Caribbean sites.
Figure 5. Boxplots of significantly different genus-level OTUs in the Atlantic and Caribbean sites.
Microbiolres 11 00009 g005
Table 1. Description of the 11 samples based on the sampling area, environment, and number of sequences and operational taxonomic units (OTUs).
Table 1. Description of the 11 samples based on the sampling area, environment, and number of sequences and operational taxonomic units (OTUs).
Sample IDSiteCoordinatesSampleRegionNumber of SequencesNumber of OTUs
C6Cabo RojoLat: 18° 05′ 15″ N, Lon: 067° 08′ 49″ WWaterCaribbean21,4932268
F9YabucoaLat: 18° 03′ 01″ N, Lon: 066° 52′ 45″ W25,0222716
A10HumacaoLat: 18° 09′ 18″ N, Lon: 065° 49′ 13″ WAtlantic25,3662162
L1FajardoLat: 18° 19′ 32″ N, Lon: 065° 39′ 08″ W22,7882461
J2San JuanLat: 18° 27′ 55″ N, Lon: 066° 06′ 20″ W25,5832618
Y03AreciboLat: 18° 20′ 47″ N, Lon: 066° 45′ 10″ W57,1102268
L8SalinasLat: 17° 58′ 39″ N, Lon: 066° 17′ 452″ WSedimentCaribbean29,6052027
J22PonceLat: 18° 00′ 39″ N, Lon: 066° 36′ 50″ W27,6702239
N6MayagüezLat: 18° 12′ 04″ N, Lon: 067° 08′ 42″ WAtlantic29,0662231
N41AguadillaLat: 18° 25′ 38″ N, Lon: 067° 09′ 14″ W29,5372102
Y31AreciboLat: 18° 20′ 47″ N, Lon: 066° 45′ 10″ W31,1341052
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Louime, C.J.; Vazquez-Sanchez, F.; Derilus, D.; Godoy-Vitorino, F. Divergent Microbiota Dynamics along the Coastal Marine Ecosystem of Puerto Rico. Microbiol. Res. 2020, 11, 45-55. https://0-doi-org.brum.beds.ac.uk/10.3390/microbiolres11020009

AMA Style

Louime CJ, Vazquez-Sanchez F, Derilus D, Godoy-Vitorino F. Divergent Microbiota Dynamics along the Coastal Marine Ecosystem of Puerto Rico. Microbiology Research. 2020; 11(2):45-55. https://0-doi-org.brum.beds.ac.uk/10.3390/microbiolres11020009

Chicago/Turabian Style

Louime, Clifford Jaylen, Frances Vazquez-Sanchez, Dieunel Derilus, and Filipa Godoy-Vitorino. 2020. "Divergent Microbiota Dynamics along the Coastal Marine Ecosystem of Puerto Rico" Microbiology Research 11, no. 2: 45-55. https://0-doi-org.brum.beds.ac.uk/10.3390/microbiolres11020009

Article Metrics

Back to TopTop