Next Article in Journal
Bioinformatics and Experimental Analyses Reveal NFIC as an Upstream Transcriptional Regulator for Ischemic Cardiomyopathy
Next Article in Special Issue
Variability in the Pre-Analytical Stages Influences Microbiome Laboratory Analyses
Previous Article in Journal
MarZIC: A Marginal Mediation Model for Zero-Inflated Compositional Mediators with Applications to Microbiome Data
Previous Article in Special Issue
Gut Microbiota as a Potential Predictive Biomarker in Relapsing-Remitting Multiple Sclerosis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Microbial Richness of Marine Biofilms Revealed by Sequencing Full-Length 16S rRNA Genes

College of Marine Life Sciences, Ocean University of China, Qingdao 266005, China
*
Author to whom correspondence should be addressed.
Genes 2022, 13(6), 1050; https://doi.org/10.3390/genes13061050
Submission received: 4 May 2022 / Revised: 6 June 2022 / Accepted: 7 June 2022 / Published: 12 June 2022
(This article belongs to the Special Issue When Genes Meet Microbial Ecology and Evolution)

Abstract

:
Marine biofilms are a collective of microbes that can grow on many different surfaces immersed in marine environments. Estimating the microbial richness and specificity of a marine biofilm community is a challenging task due to the high complexity in comparison with seawater. Here, we compared the resolution of full-length 16S rRNA gene sequencing technique of a PacBio platform for microbe identification in marine biofilms with the results of partial 16S rRNA gene sequencing of traditional Illumina PE250 platform. At the same time, the microbial richness, diversity, and composition of adjacent seawater communities in the same batch of samples were analyzed. Both techniques revealed higher species richness, as reflected by the Chao1 index, in the biofilms than that in the seawater communities. Moreover, compared with Illumina sequencing, PacBio sequencing detected more specific species for biofilms and less specific species for seawater. Members of Vibrio, Arcobacter, Photobacterium, Pseudoalteromonas, and Thalassomonas were significantly enriched in the biofilms, which is consistent with the previous understanding of species adapted to a surface-associated lifestyle and validates the taxonomic analyses in the current study. To conclude, the full-length sequencing of 16S rRNA genes has probably a stronger ability to analyze more complex microbial communities, such as marine biofilms, the species richness of which has probably been under-estimated in previous studies.

Graphical Abstract

1. Introduction

Marine microbial communities in the ocean account for more than 80% of biomass on earth and have crucial ecological roles [1,2]. For example, marine microorganisms are the most important primary producers in the ocean; they convert carbon dioxide into organic matter through photosynthesis and chemosynthesis and support the whole marine ecosystem through food chains and food webs [2,3]. At the same time, microbes are considered to be the main decomposers in ecosystems, which reconvert complex organic matter into relatively simple organics [4,5]. For these reasons, there is a strong motivation to study the diversity of marine microorganisms. The methods originally used to study microbial diversity were mostly direct observation under the microscope and plate counting [6], but assessing marine microbial diversity is a challenging task due to their sheer variety and resistance to cultivation [7].
In the marine environment, a biofilm can be defined as a community of microbes attached to the surface of a solid material, such as the surface of rocks, floating artificial panels, and animal guts [8,9,10]. Zhang et al. collected biofilms from various substrates, including natural and artificial materials, in eight countries around the world and obtained 101 biofilm metagenomes using Illumina sequencing. The results indicated very high microbial richness in the biofilm communities, which contained microbes (more than 7300 species as revealed by the miTags extracted from the 101 biofilm metagenomes) with undetectable abundance in the seawater [11]. As such, biofilms hugely supplement marine microbial richness and represent bounteous reservoirs for exploring marine microbial resources.
The emergence of DNA-sequencing technologies at the end of last century has made it possible to directly detect microbial diversity in natural environments, which has greatly benefited the field of microbial ecology [6]. So far, three main techniques have been used to sequence microbial diversity. The Illumina platforms that are used to sequence environmental metagenomes are effective means of characterizing both taxonomic and functional diversity; 16S rRNA genes extracted from metagenomes are defined as miTags, which have the advantage of more accurately reflecting real community richness [12], and as mentioned above, this method has been used in analyzing marine biofilms [11]. Illumina platforms can also be employed to detect target variable regions in the highly conserved 16S rRNA genes after PCR using universal primers. The V3–V4 and V4–V5 variable regions are the most common sequences used in the study of natural environment samples [13,14,15]. The V2–V3 variable regions have higher resolution for lower taxonomic levels (i.e., genus and species level) [16], and the V6–V8 variable regions have better recognition ability for lower abundance taxa [17], whereas the V1–V3 [18,19], V3–V5 [20], V4–V6 [21], and V5–V9 [22] regions are also used for profiling microbial community compositions under specific environments. In third-generation sequencing techniques, i.e., PacBio sequencing, the read length completely covers the full-length 16S rRNA genes, with an average length of 1.5 kb, providing more information on native DNA bases [23,24]. Studies have shown that PacBio sequencing can detect more operational taxonomic units (OTUs) with higher resolution at the genus and species levels and demonstrated it to be a useful way to profile marine microbial communities [25,26].
The major purpose of this study is to analyze the microbial richness, diversity, and composition of marine biofilms through sequencing full-length 16S rRNA genes and the V3–V4 regions; we also test whether PacBio sequencing is more effective in dealing with complex microbial communities, through comparing the results of biofilms with the results of seawater communities. After extracting DNA from the biofilm and seawater samples, 16S rRNA gene amplicon sequencing was conducted using PacBio and Illumina techniques to sequence the full-length and the V3–V4 region of 16S rRNA genes, respectively. Then, the microbial richness, diversity, and community structure were comparatively analyzed.

2. Methods

2.1. Biofilm and Seawater Sampling

Following the steps described in our previous study [11], biofilms were developed on polystyrene Petri dishes with a diameter of 9 cm. The Petri dishes were placed in a 20 × 30 cm nylon net bag and then put into the subtidal zone along the coast of Qingdao (120.145, 39.915). After they were immersed in seawater for two weeks, the panels were transferred to the laboratory and subjected to sample preparation immediately. Briefly, the biofilms attached to the plate were scraped off with a cotton swab and put into 50 mL centrifuge tubes with artificial seawater to vibrate for 15 min. The seawater samples were filtered with 0.22-micron filter membranes, which were then cut, crushed, and put into 50 mL centrifuge tubes filled with artificial seawater to vibrate for 15 min. Finally, the centrifuge tubes were put into a centrifuge (6000 rpm) to obtain bacterial cells for DNA extraction.

2.2. DNA Extraction, PCR, and Sequencing

DNA from the biofilm and seawater samples was extracted using a One-Step-Lysis Bacteria DNA Kit (Nobelab Biotechnology, Beijing, China). For the PacBio sequencing, a 50 μL amplification system was prepared, which consisted of 25 μL of Prime Star Max premix (R045, Takara, Beijing, China), 2.5 μL of 27F primers (5′-AGAGTTTGATCMTGGCTCAG-3′) and 2.5 μL of 1492R primers (5′-GGTTACCTTGTTACGACTT-3′) with six different barcodes, 1 μL of DNA template, and distilled water to the remaining volume. The PCR amplification conditions included a denaturation step at 95 °C for 5 min, 30 cycles of amplification (95 °C for 30 s, 57 °C for 30 s, and 72 °C for 2 min), and a final extending step at 72 °C for 5 min. PCR products were examined using agarose gel electrophoresis and a Qubit assay before being processed on a PacBio RS II sequencing platform by Novogene (Beijing, China). The sequencing depth was 100,000 reads per sample. For Illumina sequencing, the same PCR procedures were conducted using the 341F primers (5′-CCTACGGGNGGCWGCAG-3′) and 805R primers (5′-GACTACHVGGGTATCTAATCC-3′) with six different pairs of barcodes. The PCR products were sequenced on an Illumina PE250 platform by Novogene (Beijing, China), with a sequencing depth of 100,000 pair-ended reads per sample.

2.3. Data Correction and Filtration

The offline data from PacBio sequencing were distinguished and segregated based on the different barcode sequences added in the PCR process. The raw sequences were saved in bam format without correction. Then, circular consensus sequencing (CCS) (SMRT link v7.0) was used to correct the sequences. The correction parameters were the following: CCS = 3, minimum accuracy = 0.99. The sequence lengths ranged between 1340 bp and 1640 bp. The files were stored in fastq and fasta formats. Sequences containing consecutive identical base numbers of more than 8 were filtered out using cutadapt [27] to generate the final valid data (clean reads) for analysis. Pair-ended Illumina reads were assigned to samples based on their unique barcodes and truncated by cutting off the barcode and primer sequences. FLASH (v1.2.7) software was used to splice the overlapping read pairs [28]. QIIME2 [29] was used to filter out sequences with undetected bases (represented by N) or low-quality bases (score < 75%) and, finally, the chimeric sequences generated by PCR errors. The reads obtained after the above processing steps were considered the final valid data (clean reads).

2.4. Data Analysis

Parallel-META3 (v3.5) [30] was employed to classify the PacBio and Illumina clean reads using a similarity cutoff of 97%. The mapped OTU IDs in the Parallel-META3 database file were used as the first column in the construction of an OTU table that integrated both PacBio and Illumina reads. OTUs that contained less than two reads (total number) were removed from further analysis. The relative abundance of a given OTU was defined as the read number it contained. The relative abundance of a given phylum or genus was indicated by summarizing the read numbers contained by all the OTUs that were classified into the same phylum or genus. The multiple_rarefactions.py, α_diversity.py, and collate_α.py script documented in QIIME2 [29] were used to draw rarefaction curves at intervals of 1000 sequences with 10 replicated calculations. α-diversities based on the Shannon diversity (H index) and Chao1 richness were calculated. β-diversity, as indicated by the principal coordinate analysis (PCoA), was performed using PAST (v4.03) [31] and checked using ASAP2 [32] after transforming the OTU compositions into a Bray–Curtis distance matrix. Finally, a statistical analysis was performed using a two-tailed Student’s t-test after confirming that the data were normally distributed.

3. Results

3.1. Richness, Diversity, and Rarefaction Curves

PacBio RS II was used to sequence the full-length 16S rRNA genes of six samples, including three biofilm and three seawater samples, which generated 140,859, 191,358, 178,835, 154,650, 130,362, and 126,329 clean reads for Biofilm 1, 2, and 3 and Seawater 1, 2, and 3, respectively (Table 1). Illumina PE250 was used to sequence the 16S rRNA genes (V3–V4 region) of the same samples, which generated 97,651, 97,487, 99,715, 97,732, 96,248, and 104,682 respective clean reads (Table 1). Hereafter, the PacBio data were designated PB-Biofilms and PB-Seawater, while the Illumina data were designated Illu-Biofilms and Illu-Seawater. The OTU table comprising all the datasets was constructed by mapping the sequences to a reference OTU database documented in Parallel-META3 [30]. The detected OTUs ranged from 1203 to 3937 in the 12 data files. Overall, more OTUs were detected in the PacBio-derived biofilm data (3851 OTUs on average), which detected 34.306% more OTUs in comparison with the Illumina-derived biofilm data (2867 OTUs on average) (Table 1). However, for the seawater community, Illumina sequencing detected more OTUs than the PacBio sequencing did (1711 vs. 1399) (Table 1).
Rarefaction curves were drawn based on the OTU table to examine whether the sequencing depth was sufficient to reflect species richness and to compare diversities. Chao1 richness as well as Shannon diversity (H index) were calculated based on the PacBio and Illumina datasets, followed by rarefaction curve construction. In all the examined samples, the curves reached saturation, suggesting that the sequencing depth was enough to cover most species (Figure 1). In most of the comparisons, α-diversity values were higher for the biofilms than for the seawater communities; interestingly, however, Illumina sequencing generated higher Shannon diversity for the seawater than for the biofilms (Figure 1), suggesting that microbial community evenness is influenced by sequencing type. On the other hand, the biofilms’ Chao1 and Shannon values generated by PacBio sequencing were higher than those generated by Illumina sequencing; whereas, for the seawater samples, Illumina sequencing revealed similar or even higher microbial diversity than PacBio sequencing did (Figure 1).

3.2. PCoA

To reveal similarities between the datasets generated by different techniques, as well as between the biofilm and seawater samples, a β-diversity analysis was performed using PCoA. Based on Bray–Curtis distances, which consider both the identity and abundance of the OTUs, the six biofilm datasets were separated from the six seawater datasets along the PC1 axis, which explained 34.026% of the difference (Figure 2). Moreover, dissimilarities between the PB-Biofilms and Illu-Biofilms were observed, as was the case for the PB-Seawater and Illu-Seawater communities (Figure 2).

3.3. Venn Analysis

To further visualize the dissimilarity between biofilm and seawater samples, a Venn analysis was conducted. OTUs from the three PB-Biofilms were pooled to compare with the three combined PB-Seawater samples. As shown by the results, 2003 (35.75%) OTUs were shared by the biofilm and seawater groups, while 2809 (50.13%) and 791 OTUs (14.12%) were specific to the biofilm group and seawater group, respectively (Figure 3A). In parallel, OTUs from the three Illu-Biofilms were pooled to compare with the combined Illu-Seawater samples, which revealed that 1855 (37.82%) OTUs were shared by the biofilm and seawater groups, while 2167 (44.18%) and 883 OTUs (18.00%) were specific to the biofilm group and seawater group, respectively (Figure 3B). Taken together, PacBio sequencing revealed the presence of more specific OTUs for biofilms and less specific species for seawater communities.

3.4. Taxonomic Composition

To further study the dissimilarity between the biofilm and seawater samples, the taxonomic structure was profiled at the phylum level using the OTU table. Based on the PacBio datasets, 52 phyla were identified in the biofilms, whereas 47 phyla were identified in the seawater samples. The biofilms were primarily composed of Proteobacteria (70.958–79.604%), Bacteroides (8.979–15.448%), and Cyanobacteria (3.526–8.301%), whereas the seawater samples were predominantly composed of Proteobacteria (59.787–60.258%), Cyanobacteria (10.993–16.527%), Bacterioidetes (12.263–15.171%), and Actinobacteria (7.238–9.204%) (Figure 4A). PacBio sequencing showed that 44 phyla were shared among the biofilms and seawater communities, while AncK6, FCPU426, WPS_2, Caldithrix, WS5, Hyd24_12, KSB3, and Armatimonadetes were only present in the biofilms, and Crenarchaeota, Deferribacteres, and WS1 were only present in the seawater communities (Figure 4A). In parallel, Illumina sequencing revealed the presence of 57 and 52 phyla in the biofilm and seawater samples, respectively. Proteobacteria (84.310–94.900%) and Bacteroidetes (2.230–6.807%) made up most of the read sequences for the biofilms, whereas Proteobacteria (45.451–60.317%), Actinobacteria (15.582–19.125%), Cyanobacteria (5.992–13.110%), and Bacteroidetes (9.567–12.814%) were responsible for most of the read sequences for the seawater communities (Figure 4B). Illumina sequencing showed that 48 phyla were present in both biofilms and seawater, while Anck6, LD1, WWE1, WS2, OP1, KSB3, Poribacteria, Fcpu426, and Caldiserica were specific to the biofilms, and TPD_58 sar406, AC1, and Thermotogae were specific to the seawater communities (Figure 4B). Taken together, the results of both sequencing techniques indicated that there were more phyla in the biofilms than the seawater communities, and the biofilms possessed more specific phyla than the seawater communities. For both the biofilms and seawater communities, Illumina sequencing indicated the presence of more phyla than PacBio sequencing.
OTUs were then classified to the genus level. Based on the PacBio datasets, 678 genera were identified in the biofilms and were dominated by Vibrio, which matched about half of the sequences (40.221–49.222%), followed by Arcobacter (4.110–7.097%) and Photobacterium (3.177–5.942%). According to the PacBio datasets, only 526 genera appeared in the seawater samples, including Candidatus_Portiera (24.113–29.628%), Candidatus_Pelagibacter (11.571–18.514%), Candidatus_Actinomarina (6.452–8.675%), and Synechococcus (3.304–5.345%) (Figure 5A). PacBio sequencing showed that there were 225 genera specific for the biofilms and 73 genera specific for the seawater communities. In parallel, the Illumina datasets revealed 605 genera in the biofilms, which were dominated by Vibrio (52.511–66.727%), followed by Catenoccus (7.761–8.412%), Pseudoalteromonas (4.337–7.389%), and Photobacterium (4.843–7.138%) (Figure 5B). According to the Illumina datasets, there were 537 genera in the seawater samples, including Candidatus_Portiera (14.403–18.685%), Candidatus_Actinomarina (12.379–15.769%), Roseovarius (8.350–14.460%), Synechococcus (3.284–11.054%), and NS4 (2.490–7.841%). Illumina sequencing showed that 174 genera were specific to the biofilms, while 105 genera were specific to seawater. Taken together, the results of both sequencing techniques indicated that the biofilms possessed more genera as well as more specific genera than the seawater communities. Interestingly, PacBio sequencing of the biofilms revealed more genera than Illumina sequencing did, whereas PacBio sequencing of the seawater communities revealed fewer genera than Illumina sequencing did.
We further performed a statistical analysis to identify genera that were significantly enriched in the biofilms in comparison with the seawater communities. According to the PacBio datasets, Vibrio, Thalassomonas, Pseudoalteromonas, Photobacterium, Arcobacter, and Aquibacter were significantly (p-value <0.01) enriched in the biofilms (Figure 6A). In particular, the PacBio datasets suggested that Vibrio species were enriched by more than 100-fold in the biofilms (53.64% in biofilms compared to 0.52% in the seawater communities) (Figure 6A). According to the Illumina datasets, Vibrio, Pseudoalteromonas, and Photobacterium were also significantly enriched in the biofilms (Figure 6B). In addition, the Illumina datasets revealed Catenococcus to be a biofilm-enriched genus instead of Arcobacter (not detected) and Aquibacter (enriched, but not significantly) (Figure 6B).

4. Discussion

In the present study, we performed comparative analyses of the microbial richness, diversity, and composition of biofilms and adjacent seawater samples, based on 16S rRNA gene sequences generated by PacBio and Illumina sequencing. The results revealed higher microbial richness in the biofilms than the seawater communities. For the biofilms, PacBio sequencing revealed more species than Illumina sequencing did. The taxonomy was further analyzed to indicate biofilm specificity as well as the advantages of the PacBio sequencing for profiling complex microbial compositions.
The higher richness of the marine biofilms compared with the seawater communities observed in the present study and several previous studies [11,33] can be attributed to biofilms facilitating the proliferation of many species that are very rare in seawater. Such proliferation is probably related to the special biofilm architecture, which provides suitable physiological and chemical conditions, as well as ecological advantages, for these microbes. It is well known that biofilm matrixes formed by the polymerization of a large number of metabolites can improve the utilization efficiency of internal resources and signal transduction [34,35] and also improve the tolerance of bacteria in the community to external stresses, including metal ions and oxidative stresses [33,36]. Moreover, the uneven distribution of polymers and microbes within a biofilm facilitates the formation of many microenvironments with different oxygen gradients. The formation of such oxygen gradients allows the coexistence of microbes with different respiration abilities, i.e., preferences for different electron acceptors [37,38]. Actually, microbes can cooperate with each other to promote electron transduction and energy generation, which are important driving forces of biofilm development in environments with poor carbon and energy sources [39,40]. In line with these views, both PacBio and Illumina sequencing approaches identified Vibrio, Catenococcus, and Pseudoalteromonas as dominant and specific biofilm inhabitants in the present study, suggesting that these genera have evolved surface-associated lifestyles. A biofilm-forming ability brings survival advantages, such as enabling Vibrio cholerae to survive protozoan grazing [41] and enhancing Vibrio virulence [42]. Pseudoalteromonas species also have a strong ability to form biofilms. Many members of Pseudoalteromonas secrete anti-metabolites and can form a protective film inside the biofilm to enhance the resistance of the whole community to antibiotics [43]. At the same time, several Pseudoalteromonas species have antibacterial abilities and can parasitize the hosts to form biofilms that help the host fight against pathogenic bacteria [43]. The detection of these microbes has also validated the taxonomic profiling results in the present study.
Compared with Illumina sequencing, PacBio sequencing detected more species as well as more specific species in the biofilms. These results are consistent with previous conclusions that PacBio has great advantages for categorization to the species level and is more accurate than techniques measuring some variable regions [25,26]. On the other hand, the PCoA showed substantial dissimilarity between the PacBio- and Illumina-sequenced biofilms, and only 18 of the top 30 genera in the PacBio-sequenced biofilms were recalled by Illumina sequencing. Such a substantial difference may have arisen because of the full-length sequencing of 16S RNA, which not only contains the sequences comprising the V3–V4 variable region, but also information on all other variable regions. This view is also supported by findings that show that 16S rRNA regions have different advantages for identifying distinct microbes, largely depending on their taxonomy. For example, the V1–V2 region can be used to accurately identify a large number of clinically important strains [44,45], whereas the V1–V3 variable region can reveal a large number of endemic microbial groups [46]. However, it is worth pointing out that the primers used for the full-length 16S rRNA amplification and PacBio sequencing (27F/1492R) fail to fully cover the conserved regions of archaea, whose identification often rely on the 23F/1492R primer pair [47,48], and this also contributes to sequencing technique-generated taxonomic differences. Moreover, considering the different microbial richness between biofilms and seawater communities, the advantages of different sequencing strategies are probably related to the complexity of the target microbial communities.

5. Conclusions

To conclude, this work has highlighted the complexity of marine biofilms with extremely high and specific microbial richness in comparison with free-living microbial communities in seawater. Considering that most published microbial community research has focused on seawater, future studies on biofilms should further change our understanding of the ecology and lifestyle of marine microbes. More importantly, we found that using PacBio sequencing can better reflect the profiles of complex microbial communities, such as biofilms. In addition, we believe that the development and application of novel DNA sequencing techniques will continuously refresh our understanding of the microbial richness in the global ocean, which might have been severely underestimated in global sampling projects performed in past years [49,50].

Author Contributions

Experiment design, W.Z., X.H., W.D. and S.W.; experiments, S.W. and X.S.; software, H.C. and M.W.; sampling resources, X.H. and W.D.; data analyses, S.W. and W.Z.; manuscript writing, S.W.; funding, W.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ocean University of China (Funding No. 842041010) and the Organization Department of China (Funding No. 862105020028).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available on the NCBI database (Bioproject No. PRJNA816818).

Conflicts of Interest

The authors declare that there are no conflict of interest.

References

  1. Azam, F.; Worden, A.Z. Microbes, molecules, and marine ecosystems. Science 2004, 303, 1622–1624. [Google Scholar] [CrossRef]
  2. Cavicchioli, R.; Ripple, W.J.; Timmis, K.N.; Azam, F.; Bakken, L.R.; Baylis, M.; Behrenfeld, M.J.; Boetius, A.; Boyd, P.W.; Classen, A.T.; et al. Scientists’ warning to humanity: Microorganisms and climate change. Nat. Rev. Microbiol. 2019, 17, 569–586. [Google Scholar] [CrossRef]
  3. Verde, C.; Giordano, D.; Bellas, C.M.; di Prisco, G.; Anesio, A.M. Polar marine microorganisms and climate change. Adv. Microb. Physiol. 2016, 69, 187–215. [Google Scholar] [CrossRef]
  4. Palit, K.; Rath, S.; Chatterjee, S.; Das, S. Microbial diversity and ecological interactions of microorganisms in the mangrove ecosystem: Threats, vulnerability, and adaptations. Environ. Sci. Pollut. Res. Int. 2022, 29, 32467–32512. [Google Scholar] [CrossRef]
  5. Hutchins, D.A.; Fu, F. Microorganisms and ocean global change. Nat. Microbiol. 2017, 2, 17058. [Google Scholar] [CrossRef]
  6. Salazar, G.; Sunagawa, S. Marine microbial diversity. Curr. Biol. 2017, 27, R489–R494. [Google Scholar] [CrossRef]
  7. Kaeberlein, T.; Lewis, K.; Epstein, S.S. Isolating “uncultivable” microorganisms in pure culture in a simulated natural environment. Science 2002, 296, 1127–1129. [Google Scholar] [CrossRef]
  8. Zhang, W.; Watanabe, H.K.; Ding, W.; Lan, Y.; Tian, R.-M.; Sun, J.; Chen, C.; Cai, L.; Li, Y.; Oguri, K.; et al. Gut microbial divergence between two populations of the hadal amphipod Hirondellea gigas. Appl. Environ. Microbiol. 2019, 85, e02032-18. [Google Scholar] [CrossRef]
  9. Hutchinson, N.; Nagarkar, S.; Aitchison, J.C.; Williams, G.A. Microspatial variation in marine biofilm abundance on intertidal rock surfaces. Aquat. Microb. Ecol. 2006, 42, 187–197. [Google Scholar] [CrossRef]
  10. Salta, M.; Wharton, J.A.; Blache, Y.; Stokes, K.R.; Briand, J.F. Marine biofilms on artificial surfaces: Structure and dynamics. Environ. Microbiol. 2013, 15, 2879–2893. [Google Scholar] [CrossRef]
  11. Zhang, W.; Ding, W.; Li, Y.X.; Tam, C.; Bougouffa, S.; Wang, R.; Pei, B.; Chiang, H.; Leung, P.; Lu, Y.; et al. Marine biofilms constitute a bank of hidden microbial diversity and functional potential. Nat. Commun. 2019, 10, 517. [Google Scholar] [CrossRef]
  12. Logares, R.; Sunagawa, S.; Salazar, G.; Cornejo-Castillo, F.M.; Ferrera, I.; Sarmento, H.; Hingamp, P.; Ogata, H.; de Vargas, C.; Lima-Mendez, G.; et al. Metagenomic 16S rDNA Illumina tags are a powerful alternative to amplicon sequencing to explore diversity and structure of microbial communities. Environ. Microbiol. 2014, 16, 2659–2671. [Google Scholar] [CrossRef]
  13. Fadeev, E.; Cardozo-Mino, M.G.; Rapp, J.Z.; Bienhold, C.; Salter, I.; Salman-Carvalho, V.; Molari, M.; Tegetmeyer, H.E.; Buttigieg, P.L.; Boetius, A. Comparison of two 16S rRNA primers (V3-V4 and V4-V5) for studies of Arctic microbial communities. Front. Microbiol. 2021, 12, 637526. [Google Scholar] [CrossRef]
  14. Hoshino, T.; Doi, H.; Uramoto, G.I.; Wormer, L.; Adhikari, R.R.; Xiao, N.; Morono, Y.; D’Hondt, S.; Hinrichs, K.U.; Inagaki, F. Global diversity of microbial communities in marine sediment. Proc. Natl. Acad. Sci. USA 2020, 117, 27587–27597. [Google Scholar] [CrossRef]
  15. Guo, K.; Wu, N.; Li, W.; Baattrup-Pedersen, A.; Riis, T. Microbial biofilm community dynamics in five lowland streams. Sci. Total Environ. 2021, 798, 149169. [Google Scholar] [CrossRef]
  16. Bukin, Y.S.; Galachyants, Y.P.; Morozov, I.V.; Bukin, S.V.; Zakharenko, A.S.; Zemskaya, T.I. The effect of 16S rRNA region choice on bacterial community metabarcoding results. Sci. Data 2019, 6, 190007. [Google Scholar] [CrossRef]
  17. Willis, C.; Desai, D.; LaRoche, J. Influence of 16S rRNA variable region on perceived diversity of marine microbial communities of the Northern North Atlantic. FEMS Microbiol. Lett. 2019, 366, fnz152. [Google Scholar] [CrossRef]
  18. Sun, Z.; Li, G.; Wang, C.; Jing, Y.; Zhu, Y.; Zhang, S.; Liu, Y. Community dynamics of prokaryotic and eukaryotic microbes in an estuary reservoir. Sci. Rep. 2014, 4, 6966. [Google Scholar] [CrossRef]
  19. Jordaan, K.; Bezuidenhout, C.C. Bacterial community composition of an urban river in the North West Province, South Africa, in relation to physico-chemical water quality. Environ. Sci. Pollut. Res. Int. 2016, 23, 5868–5880. [Google Scholar] [CrossRef]
  20. Cai, H.; Jiang, H.; Krumholz, L.R.; Yang, Z. Bacterial community composition of size-fractioned aggregates within the phycosphere of cyanobacterial blooms in a eutrophic freshwater lake. PLoS ONE 2014, 9, e102879. [Google Scholar] [CrossRef]
  21. Darwish, N.; Shao, J.; Schreier, L.L.; Proszkowiec-Weglarz, M. Choice of 16S ribosomal RNA primers affects the microbiome analysis in chicken ceca. Sci. Rep. 2021, 11, 11848. [Google Scholar] [CrossRef] [PubMed]
  22. Roesch, L.F.; Fulthorpe, R.R.; Riva, A.; Casella, G.; Hadwin, A.K.; Kent, A.D.; Daroub, S.H.; Camargo, F.A.; Farmerie, W.G.; Triplett, E.W. Pyrosequencing enumerates and contrasts soil microbial diversity. ISME J. 2007, 1, 283–290. [Google Scholar] [CrossRef]
  23. Hu, T.; Chitnis, N.; Monos, D.; Dinh, A. Next-generation sequencing technologies: An overview. Hum. Immunol. 2021, 82, 801–811. [Google Scholar] [CrossRef] [PubMed]
  24. Heather, J.M.; Chain, B. The sequence of sequencers: The history of sequencing DNA. Genomics 2016, 107, 1–8. [Google Scholar] [CrossRef] [PubMed]
  25. Zhang, J.; Su, L.; Wang, Y.; Deng, S. Improved high-throughput sequencing of the human oral microbiome: From Illumina to PacBio. Can. J. Infect. Dis. Med. Microbiol. 2020, 2020, 6678872. [Google Scholar] [CrossRef]
  26. Ji, B.; Zhang, X.; Zhang, S.; Song, H.; Kong, Z. Insights into the bacterial species and communities of a full-scale anaerobic/anoxic/oxic wastewater treatment plant by using third-generation sequencing. J. Biosci. Bioeng. 2019, 128, 744–750. [Google Scholar] [CrossRef]
  27. Marcel, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011, 17, 10–12. [Google Scholar] [CrossRef]
  28. Magoc, T.; Salzberg, S.L. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 2011, 27, 2957–2963. [Google Scholar] [CrossRef]
  29. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F.; et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 2019, 37, 852–857. [Google Scholar] [CrossRef]
  30. Jing, G.; Sun, Z.; Wang, H.; Gong, Y.; Huang, S.; Ning, K.; Xu, J.; Su, X. Parallel-META 3: Comprehensive taxonomical and functional analysis platform for efficient comparison of microbial communities. Sci. Rep. 2017, 7, 40371. [Google Scholar] [CrossRef]
  31. Hammer, Ø.; Harper, D.A.T.; Ryan, P.D. PAST: Paleontological statistics software package for education and data analysis. Palaeontol. Electron. 2001, 4, 4–9. [Google Scholar]
  32. Tian, R.; Imanian, B. ASAP 2: A pipeline and web server to analyze marker gene amplicon sequencing data automatically and consistently. BMC Bioinform. 2022, 23, 27. [Google Scholar] [CrossRef]
  33. Ding, W.; Zhang, W.; Alikunhi, N.M.; Batang, Z.; Pei, B.; Wang, R.; Chen, L.; Al-Suwailem, A.; Qian, P.Y. Metagenomic analysis of zinc surface-associated marine biofilms. Microb. Ecol. 2019, 77, 406–416. [Google Scholar] [CrossRef] [PubMed]
  34. Wang, R.; Ding, W.; Long, L.; Lan, Y.; Tong, H.; Saha, S.; Wong, Y.H.; Sun, J.; Li, Y.; Zhang, W.; et al. Exploring the influence of signal molecules on marine biofilms development. Front. Microbiol. 2020, 11, 571400. [Google Scholar] [CrossRef] [PubMed]
  35. Wang, R.; Zhang, W.; Ding, W.; Liang, Z.; Long, L.; Wong, W.C.; Qian, P.Y. Profiling signal transduction in global marine biofilms. Front. Microbiol. 2021, 12, 768926. [Google Scholar] [CrossRef] [PubMed]
  36. Zhang, W.; Wang, Y.; Lee, O.O.; Tian, R.; Cao, H.; Gao, Z.; Li, Y.; Yu, L.; Xu, Y.; Qian, P.Y. Adaptation of intertidal biofilm communities is driven by metal ion and oxidative stresses. Sci. Rep. 2013, 3, 3180. [Google Scholar] [CrossRef] [PubMed]
  37. Stewart, P.S.; Franklin, M.J. Physiological heterogeneity in biofilms. Nat. Rev. Microbiol. 2008, 6, 199–210. [Google Scholar] [CrossRef]
  38. Zhou, Y.; Li, R.; Guo, B.; Zhang, L.; Zhang, H.; Xia, S.; Liu, Y. Three-dimension oxygen gradient induced low energy input for grey water treatment in an oxygen-based membrane biofilm reactor. Environ. Res. 2020, 191, 110124. [Google Scholar] [CrossRef]
  39. Zhang, W.; Ding, W.; Yang, B.; Tian, R.; Gu, S.; Luo, H.; Qian, P.Y. Genomic and transcriptomic evidence for carbohydrate consumption among microorganisms in a cold seep brine pool. Front. Microbiol. 2016, 7, 1825. [Google Scholar] [CrossRef]
  40. Zhang, W.; Wang, Y.; Bougouffa, S.; Tian, R.; Cao, H.; Li, Y.; Cai, L.; Wong, Y.H.; Zhang, G.; Zhou, G.; et al. Synchronized dynamics of bacterial niche-specific functions during biofilm development in a cold seep brine pool. Environ. Microbiol. 2015, 17, 4089–4104. [Google Scholar] [CrossRef]
  41. Matz, C.; McDougald, D.; Moreno, A.M.; Yung, P.Y.; Yildiz, F.H.; Kjelleberg, S. Biofilm formation and phenotypic variation enhance predation-driven persistence of Vibrio cholerae. Proc. Natl. Acad. Sci. USA 2005, 102, 16819–16824. [Google Scholar] [CrossRef] [PubMed]
  42. Silva, A.J.; Benitez, J.A. Vibrio cholerae biofilms and cholera pathogenesis. PLoS Negl. Trop. Dis. 2016, 10, e0004330. [Google Scholar] [CrossRef]
  43. Offret, C.; Desriac, F.; Le Chevalier, P.; Mounier, J.; Jegou, C.; Fleury, Y. Spotlight on antimicrobial metabolites from the marine bacteria Pseudoalteromonas: Chemodiversity and ecological significance. Mar. Drugs 2016, 14, 129. [Google Scholar] [CrossRef] [PubMed]
  44. Watts, G.S.; Youens-Clark, K.; Slepian, M.J.; Wolk, D.M.; Oshiro, M.M.; Metzger, G.S.; Dhingra, D.; Cranmer, L.D.; Hurwitz, B.L. 16S rRNA gene sequencing on a benchtop sequencer: Accuracy for identification of clinically important bacteria. J. Appl. Microbiol. 2017, 123, 1584–1596. [Google Scholar] [CrossRef]
  45. Sune, D.; Rydberg, H.; Augustinsson, A.N.; Serrander, L.; Jungestrom, M.B. Optimization of 16S rRNA gene analysis for use in the diagnostic clinical microbiology service. J. Microbiol. Methods 2020, 170, 105854. [Google Scholar] [CrossRef] [PubMed]
  46. O’Dea, C.; Huerlimann, R.; Masters, N.; Kuballa, A.; Veal, C.; Fisher, P.; Stratton, H.; Katouli, M. Microbial diversity profiling of gut microbiota of Macropus giganteus using three hypervariable regions of the bacterial 16S rRNA. Microorganisms 2021, 9, 1721. [Google Scholar] [CrossRef]
  47. Vlasceanu, L.; Sarbu, S.M.; Engel, A.S.; Kinkle, B.K. Acidic cave-wall biofilms located in the Frasassi Gorge, Italy. Geomicrobiol. J. 2000, 17, 125–139. [Google Scholar] [CrossRef]
  48. Baker, B.J.; Comolli, L.R.; Dick, G.J.; Hauser, L.J.; Hyatt, D.; Dill, B.D.; Land, M.L.; Verberkmoes, N.C.; Hettich, R.L.; Banfield, J.F. Enigmatic, ultrasmall, uncultivated Archaea. Proc. Natl. Acad. Sci. USA 2010, 107, 8806–8811. [Google Scholar] [CrossRef]
  49. Cao, S.; Zhang, W.; Ding, W.; Wang, M.; Fan, S.; Yang, B.; McMinn, A.; Wang, M.; Xie, B.B.; Qin, Q.L.; et al. Structure and function of the Arctic and Antarctic marine microbiota as revealed by metagenomics. Microbiome 2020, 8, 47. [Google Scholar] [CrossRef]
  50. 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. Structure and function of the global ocean microbiome. Science 2015, 348, 1261359. [Google Scholar] [CrossRef]
Figure 1. Rarefaction curves showing sequence number and α-diversity correlations. The three biofilm and three seawater samples are indicated by different colors. Chao1 richness and Shannon (H index) diversity were calculated based on the PacBio (A) and the Illumina datasets (B). For each sample, 100,000 sequences were included for analysis, and the average value of 10 permutations is shown.
Figure 1. Rarefaction curves showing sequence number and α-diversity correlations. The three biofilm and three seawater samples are indicated by different colors. Chao1 richness and Shannon (H index) diversity were calculated based on the PacBio (A) and the Illumina datasets (B). For each sample, 100,000 sequences were included for analysis, and the average value of 10 permutations is shown.
Genes 13 01050 g001
Figure 2. Principal coordinate analysis (PCoA) based on OTU compositions of biofilms and seawater communities. Bray–Curtis distances were used in the calculation. PB-Biofilm and PB-Seawater represent the three biofilm and three seawater samples sequenced by PacBio, respectively. Illu-Biofilm and Illu-Seawater represent the same samples sequenced by Illumina.
Figure 2. Principal coordinate analysis (PCoA) based on OTU compositions of biofilms and seawater communities. Bray–Curtis distances were used in the calculation. PB-Biofilm and PB-Seawater represent the three biofilm and three seawater samples sequenced by PacBio, respectively. Illu-Biofilm and Illu-Seawater represent the same samples sequenced by Illumina.
Genes 13 01050 g002
Figure 3. Venn diagrams showing the distribution of OTUs across biofilm and seawater samples. The three replicates were pooled together for analysis. (A) PacBio-sequenced samples. (B) Illumina-sequenced samples.
Figure 3. Venn diagrams showing the distribution of OTUs across biofilm and seawater samples. The three replicates were pooled together for analysis. (A) PacBio-sequenced samples. (B) Illumina-sequenced samples.
Genes 13 01050 g003
Figure 4. Phylum-level microbial community composition of the six samples. Proteobacteria were classified down to the class level. The 30 most abundant genera in terms of maximum relative abundance are shown, and the rest are merged into ‘Others’. (A) Biofilm and seawater microbial communities by PacBio sequencing. (B) Biofilm and seawater microbial communities by Illumina sequencing.
Figure 4. Phylum-level microbial community composition of the six samples. Proteobacteria were classified down to the class level. The 30 most abundant genera in terms of maximum relative abundance are shown, and the rest are merged into ‘Others’. (A) Biofilm and seawater microbial communities by PacBio sequencing. (B) Biofilm and seawater microbial communities by Illumina sequencing.
Genes 13 01050 g004
Figure 5. Genus-level taxonomic structure of biofilm and seawater samples. The top 30 abundant genera in terms of maximum relative abundance are shown, and the rest are merged into ‘Others’. (A) Biofilm and seawater microbial communities by PacBio sequencing. (B) Biofilm and seawater microbial communities by Illumina sequencing.
Figure 5. Genus-level taxonomic structure of biofilm and seawater samples. The top 30 abundant genera in terms of maximum relative abundance are shown, and the rest are merged into ‘Others’. (A) Biofilm and seawater microbial communities by PacBio sequencing. (B) Biofilm and seawater microbial communities by Illumina sequencing.
Genes 13 01050 g005
Figure 6. Statistical analysis of the 30 most abundant genera in biofilm and seawater samples. There were significant differences in genera between the biofilm and the seawater communities. ‘**’ represents a p-value < 0.01, calculated by two-tailed Student’s t-test. (A) Biofilm and seawater microbial communities by PacBio sequencing. (B) Biofilm and seawater microbial communities by Illumina sequencing.
Figure 6. Statistical analysis of the 30 most abundant genera in biofilm and seawater samples. There were significant differences in genera between the biofilm and the seawater communities. ‘**’ represents a p-value < 0.01, calculated by two-tailed Student’s t-test. (A) Biofilm and seawater microbial communities by PacBio sequencing. (B) Biofilm and seawater microbial communities by Illumina sequencing.
Genes 13 01050 g006
Table 1. Reads and OTUs information.
Table 1. Reads and OTUs information.
SampleClean ReadsOTU NumberAverage OTU Number
PB-Biofilm1140,85937423851
PB-Biofilm2191,3583874
PB-Biofilm3178,8353937
PB-Seawater1154,65016901399
PB-Seawater2130,3621203
PB-Seawater3126,3291305
Illu-Biofilm197,65129532867
Illu-Biofilm297,4872462
Illu-Biofilm399,7153187
Illu-Seawater197,73216921711
Illu-Seawater296,2481837
Illu-Seawater3104,6821606
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, S.; Su, X.; Cui, H.; Wang, M.; Hu, X.; Ding, W.; Zhang, W. Microbial Richness of Marine Biofilms Revealed by Sequencing Full-Length 16S rRNA Genes. Genes 2022, 13, 1050. https://doi.org/10.3390/genes13061050

AMA Style

Wang S, Su X, Cui H, Wang M, Hu X, Ding W, Zhang W. Microbial Richness of Marine Biofilms Revealed by Sequencing Full-Length 16S rRNA Genes. Genes. 2022; 13(6):1050. https://doi.org/10.3390/genes13061050

Chicago/Turabian Style

Wang, Shougang, Xiaoyan Su, Han Cui, Meng Wang, Xiaoli Hu, Wei Ding, and Weipeng Zhang. 2022. "Microbial Richness of Marine Biofilms Revealed by Sequencing Full-Length 16S rRNA Genes" Genes 13, no. 6: 1050. https://doi.org/10.3390/genes13061050

APA Style

Wang, S., Su, X., Cui, H., Wang, M., Hu, X., Ding, W., & Zhang, W. (2022). Microbial Richness of Marine Biofilms Revealed by Sequencing Full-Length 16S rRNA Genes. Genes, 13(6), 1050. https://doi.org/10.3390/genes13061050

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop