Yunyan Deng
Kui Wang
Zhangxi Hu
Qiang Hu6
Ying Zhong Tang- 1CAS Key Laboratory of Marine Ecology and Environmental Sciences, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China
- 2Laboratory of Marine Ecology and Environmental Science, Qingdao National Laboratory for Marine Science and Technology, Qingdao, China
- 3Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao, China
- 4Institute for Advanced Study, Shenzhen University, Shenzhen, China
- 5College of Fisheries, Guangdong Ocean University, Zhanjiang, China
- 6Faculty of Synthetic Biology, CAS Key Laboratory of Quantitative Engineering Biology, Shenzhen Institute of Synthetic Biology, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China
Identification of a core microbiome (a group of taxa commonly present and consistently abundant in most samples of host populations) is important to capture the key microbes closely associated with a host population, as this process may potentially contribute to further revealing their spatial distribution, temporal stability, ecological influence, and even impacts on their host’s functions and fitness. The naked dinoflagellate Karlodinium veneficum is a cosmopolitan and toxic species, which is also notorious in forming harmful algal blooms (HABs) and causing massive fish-kills. Here we reported the core microbiome tightly associated with 19 strains of K. veneficum that were originally isolated from 6 geographic locations along the coast of China and from an estuary of Chesapeake Bay, United States, and have been maintained in the laboratory for several months to over 14 years. Using high-throughput metabarcoding of the partial 16S rRNA gene amplicons, a total of 1,417 prokaryotic features were detected in the entire bacterial microbiome, which were assigned to 17 phyla, 35 classes, 90 orders, 273 families, and 716 genera. Although the bacterial communities associated with K. veneficum cultures displayed heterogeneity in feature (sequences clustered at 100% sequence similarity) composition among strains, a core set of 6 genera were found persistent in their phycospheres, which could contribute up to 74.54% of the whole bacterial microbiome. Three γ-proteobacteria members of the “core,” namely, Alteromonas, Marinobacter, and Methylophaga, were the predominant core genera and made up 83.25% of the core bacterial microbiome. The other 3 core genera, Alcanivorax, Thalassospira, and Ponticoccus, are reported to preferably utilize hydrocarbons as sole or major source of carbon and energy, and two of which (Alcanivorax and Ponticoccus) are recognized as obligate hydrocarbonoclastic bacteria (OHCB). Since OHCB generally present in extremely low abundance in marine water and elevate their abundance mostly in petroleum-impacted water, our detection in K. veneficum cultures suggests that the occurrence of obligate and generalist hydrocarbon-degrading bacteria living with dinoflagellates may be more frequent in nature. Our work identified a core microbiome with stable association with the harmful alga K. veneficum and opened a window for further characterization of the physiological mechanisms and ecological implications for the dinoflagellate-bacteria association.
Introduction
In aquatic habitats, phytoplankton communities within the photic zones are characterized by highly heterogeneous assemblages of unicellular microalgae and a wide variety of associated microbes. The phototrophic microalgae capture solar energy, transform inorganic matter into organic matter, and release organic compounds to the region immediately surrounding themselves, known as the phycosphere, which is enriched in organic molecules exuded by the algal cell into the surrounding water, and is the aquatic analogue of the rhizosphere (Seymour et al., 2017; Durán et al., 2022). The phycosphere provides a niche for colonization of heterotrophic bacteria, which utilize a major fraction of fixed carbon through catabolism of microalgal metabolites. It could be considered as a marketplace where cross-kingdom communications are mediated by the release and uptake of organic compounds (Kouzuma and Watanabe, 2015; Seymour et al., 2017). The exchange of metabolites and infochemicals at this interface governs multifarious and sophisticated microalgae-bacteria associations, which span from cooperative to competitive relationships (e.g., mutualism, commensalism, parasitism and competition; Amin et al., 2009; Kouzuma and Watanabe, 2015; Aharonovich and Sher, 2016; Hennon et al., 2018). Since modulating nutrient cycling and biomass production at the base of the marine food web, interactions between microalgae and bacteria represent one of the fundamental ecological relationships in aquatic environments, which are considered to have multilevel influences and substantial impacts on ecosystem-scale processes (e. g. biogeochemical cycling, ecosystem productivity, harmful algal blooms; Buchan et al., 2014; Seymour et al., 2017).
The concept of “core microbiome” originally emerged for human microbiome and refers to common groups of microbes or genes that were likely to be particularly important for host biological function (Shade and Handelsman, 2012 and the references therein), and further expanded to other host-associated microbiomes. Broadly, the core microbiome has been referred to any set of microbial taxa, or the genomic and functional attributes associated with those taxa, that are characteristic of a host or environment of interest (Astudillo-Garcia et al., 2017; Risely, 2020; Neu et al., 2021). Most commonly, core microbiomes are measured as the suite of taxa shared among two or more samples from a particular host or environment (Shade and Handelsman, 2012; Astudillo-Garcia et al., 2017; Neu et al., 2021). Identifying the core members is the first step to understand the stable, consistent components across complex microbial assemblages, because those commonly occurring taxa that appear in all assemblages are proposed to be associated with a particular habitat and are critical to the function of that type of community (Shade and Handelsman, 2012; Astudillo-Garcia et al., 2017; Risely, 2020; Neu et al., 2021). For terrestrial plants, the composition and function of core microbiomes within a rhizosphere or phyllosphere have been explored in some model plants, such as Arabidopsis, maize, rice, barley, and soybean (Lundberg et al., 2012; Peiffer et al., 2013; Mendes et al., 2014; Edwards et al., 2015). Several works conveyed a notion that the soil types and host plant genotypes were the main factors affecting the associated bacterial consortia (Lundberg et al., 2012; Peiffer et al., 2013; Edwards et al., 2015; Grady et al., 2019), though a few of them addressed the effects of microbial assemblages on host metabolism (Risely, 2020 and the references therein).
Dinoflagellates are the second most abundant phytoplankton in coastal marine ecosystems and infamous for causing harmful algal blooms (HABs). Approximately 200 species in this lineage have been recorded as HABs-causing agents, accounting for ca. 40% of the total recorded species forming HABs globally (Jeong et al., 2021). With the ever-increasing number of dinoflagellate HAB events worldwide in the past decades, the importance of dinoflagellate-bacteria interactions have been increasingly emphasized due to amassed evidence of linkages between bacterioplankton and the population dynamics of dinoflagellates in the blooming field (Mayali et al., 2011; Yang et al., 2015; Zhou et al., 2018; Hattenrath-Lehmann et al., 2019; Cui et al., 2020; Kang et al., 2021; Maire et al., 2021). Meanwhile, bacteria within certain lineages (e.g., Roseobacter clade in α-proteobacteria, Marinobacter and Alteromonas clades in γ-proteobacteria) have been commonly found both in the phycospheres of laboratory-cultured dinoflagellates (Bolch et al., 2017; Park et al., 2018; Shin et al., 2018; Li et al., 2019; Guo et al., 2020; Tarazona-Janampa et al., 2020; Deng et al., 2022; Wu et al., 2022) and field samples during dinoflagellate blooms (Mayali et al., 2011; Yang et al., 2015; Zhou et al., 2018; Hattenrath-Lehmann et al., 2019; Cui et al., 2020; Kang et al., 2021), implying that specific interactions could be established between dinoflagellates and particular groups of bacteria. However, the relationships between phytoplankton and surrounding bacteria are intimate, complex, and dynamic (Amin et al., 2012; Behringer et al., 2018; Tarazona-Janampa et al., 2020; Maire et al., 2021; Wu et al., 2022). Compared with microbiota of sessile plants, the bacterial consortia co-existed with phytoplankton in aquatic habitats are particularly difficult to be accurately defined, owing to the majority of microalgae being unicellular organisms living in constantly changing and dynamic environments (Behringer et al., 2018). Furthermore, some case studies recorded that the community composition of the associated bacterial flora was correlated with the growth stage of host phytoplankton (Grossart et al., 2005; Lupette et al., 2016; Bolch et al., 2017; Park et al., 2018; Li et al., 2019), suggesting that their interactions are more complex than expected. Although previous works have highlighted some common and/or important bacterial taxa interacting with many dinoflagellate species, the association of specific bacterial taxa with multiple strains of a specific dinoflagellate species, which would suggest “intimate” associations between the hosting species and its “microbiome” (Behringer et al., 2018), has become a new research perspective and hotspot. Recent research efforts have been undertaken to investigate Symbiodiniaceae members, the crucial endosymbionts of coral reefs (Lawson et al., 2018; Maire et al., 2021). Lawson et al. (2018) characterized the bacterial communities association with a wide phylogenetic diversity of Symbiodinium cultures (18 types spanning 5 clades), showing the existence of a consistent core bacterial microbiome among different Symbiodinium clades. From 11 cultures of Symbiodiniaceae covering 6 genera and 9 species, Maire et al. (2021) identified a highly conserved suite of intracellular bacteria across all these species of Symbiodiniaceae, suggesting these bacteria are involved in Symbiodiniaceae physiology. The area of core bacterial microbiome associated with dinoflagellates is still in the mist and underexplored.
The naked Karlodinium veneficum (Ballantine) J. Larsen is a notorious harmful algal bloom-forming dinoflagellate, which has a cosmopolitan distribution covering coastal ecosystems of all continents except Antarctica (Yang et al., 2020 and the references therein). Since first been isolated form United Kingdom, numerous bloom events of K. veneficum have been reported worldwide, e.g., South Africa, Europe, North America and Australia (as reviewed in Yang et al., 2019). In China, K. veneficum blooms have increasingly occurred since 2000. K. veneficum is a toxic species that produces karlotoxins, which have lethal effect to predators (Deeds et al., 2002) and may also lead to human skin irritation, swelling, rashes and itchiness (Lim et al., 2014). In addition, the species is an omnivorous phagotroph that could obtain nutrients and energy via feeding on aquatic organisms (Yang et al., 2020), and also exhibits potent allelopathic effects on other co-occurring algae (Yang et al., 2019). Due to these characteristics, its blooms usually cause massive mortality of fishes, mussels and zooplankton (as reviewed in Yang et al., 2021). For instance, a long-lasting K. veneficum bloom along the Yellow Sea coastal waters (Sanggou Bay, Shandong Province) in 2012 caused severe impacts on the local aquaculture (Xu et al., 2012). Massive fish-kills caused by K. veneficum were frequently reported in the Chesapeake Bay, United States from 2015 to 2017, killing nearly half million fish (Maryland Department of the Environment, 2017). K. veneficum was also responsible for the fish kill in coastal Finnish waters (Ersöströmmen, Finland) in 2016, which is the first case shown to be related to this species in the Baltic Sea (Finnish Environment Institute, 2016). Given the frequent bloom events of K. veneficum and their severe impacts on fishery and aquaculture, it has drawn great attention from the scientific community and thus has been intensively studied from multiple aspects (toxicology, toxins, trophic mode, blooming dynamics, allelopathy, and life cycle). As an important facet of its ecology, the bacterial consortia harbored by this species, however, have not been adequately examined. During the past years, we have established 19 clonal cultures K. veneficum that have been isolated from 6 geographic regions along the coast of China and from an estuary of Chesapeake Bay, United States, and maintained for months to more than 14 years, which provided us a highly valuable asset in terms of the abovementioned research objective. Therefore, we characterized the bacterial communities associated with these laboratory cultures of K. veneficum with the aim to define the core (constitutive) bacterial microbiome that stably co-exists with the species. Our findings provided new insights into the possible mechanisms underpinning the possible intimate associations and their potential contributions to physiological functioning of K. veneficum.
Materials and methods
Algal cultures of Karlodinium veneficum
A total of 19 strains of K. veneficum were employed in this study. The strain KV7 + 8, was isolated from the Chesapeake Bay estuary, United States in 2006 (Table 1 and Supplementary Figure 1). The other 18 strains were isolated from 5 different geographic origins of China (Fujian, Guangdong, Guangxi, Hebei, and Shandong Provinces; Table 1 and Supplementary Figure 1), which were initially established via single-cell pipetting either from coastal waters or via resting cyst germination from 2014 to 2020 (see Table 1 for details). The samples harvested for genomic DNA isolation (see below) were performed on December 20, 2020, thus these cultures have been maintained in the laboratory for several months to over 14 years. All cultures were routinely maintained in f/2-Si medium (i.e., f/2 medium without NaSiO3; Guillard, 1975) prepared with sterilized and filtered (0.22 mm; Millipore Corp.) seawater (salinity 32–33), in an incubator that was set at a temperature of 20 ± 1°C, an light/dark cycle of 12:12 h L/D, and an irradiance intensity of 100 μmol photons m–2 s–1.
Table 1. Information on 19 clonal cultures of Karlodinium veneficum used in this study.
Sample collection and genomic deoxyribonucleic acid extraction
Cultures at exponential growth stage were inoculated into six-well culture plates (Corning, United States) containing 10 mL of the seawater-based f/2-Si medium and then incubated for 15 days under the same temperature and illumination as the routine culture conditions. Samples were then collected when all cultures were at their stationary growth stage as pre-determined. All cells in each sample (approximately 104∼105 cells) were pelleted in a 1.5 mL centrifuge tube and immediately used for genomic DNA isolation. Genomic DNA was extracted using the Plant DNA Extraction Kit (Tiangen, Beijing, China) referring to the manufacturer’s protocols. The total DNA of each sample was eluted with 50 μl TE buffer. Nuclear-free water processed through DNA extraction were used as sample blanks. The DNA quality and purity were estimated by NanoDrop™ 1000 spectrophotometer (Thermo Fisher Scientific, United States), then stored at –80°C for further PCR amplification.
28S and 16S rRNA genes amplicon sequencing
The highly variable D2 domain and parts of the more conservative D1 and D3 domains of eukaryotic 28S rRNA gene was PCR amplified using an eukaryotic universal primer set of LSU335 (5’-ACCGATAGCA(G)AACAAGTA-3’) and LSU714 (5’-TCCTTGGTCCGTGTTTCA-3’), following the PCR protocol that described in Chai et al. (2018). The V3–V4 variable region of the bacterial 16S rRNA gene was amplified using the primer set of 341F (5’-CCTACGGGNGGCWGCAG-3’) and 805R (5’-GACTACHVGGGTATCTAATCC-3’; Logue et al., 2016). The PCR protocol were performed as that described in Deng et al. (2022). All PCR reactions were conducted in 25 μl volume including 12.5 μl of 2 × Phusion® Hot Start Flex Master Mix, 2.5 μl of each primer (1 μM), and 50 ng of template DNA. Nuclease-free water served as negative control. The resulting amplicons were checked on an agarose gel electrophoresis and purified with the Gel Extraction Kit (Axygen Biosciences, United States). The size and quantity of the purified amplicon libraries were assessed on Agilent 2100 Bioanalyzer (Agilent, United States) and then pyrosequenced on the NovaSeq PE250 platform (LC-Bio Technology Company, Hangzhou, China).
Sequencing data processing and bioinformatic analyses
The quality control of raw data was performed in Fqtrim software (version 0.9.7), and the chimeric reads were further excluded using Vsearch tool (version 2.0.3; Rognes et al., 2016). Generally, raw paired-end reads with 10 bp of minimal overlapping, the length less than 200 bp, average quality score less than 20, and 20% of maximum mismatch rate were discarded. Paired-end reads was sorted into samples based on their unique barcode and merged using FLASH (version 1.2.6; Magoè and Salzberg, 2011). The amplicon sequence variants (ASVs, sequences clustered at 100% sequence similarity) were generated with DADA2 package (version 3.6.1; Callahan et al., 2016). The eukaryotic and bacterial features aligned to NCBI GenBank and SILVA databases (retrieved on April 23, 2021), respectively, were denominated at the domain, phylum, class, order, family, and genus levels. Relative abundance of each feature was estimated based on its read counts normalized to the total number of good quality reads. Alpha diversity indices (Shannon diversity, Simpson evenness, Chao1 richness, Observed species richness, and Goods coverage) and beta diversity of PCoA (Principal coordinated analysis) based on the weighted-uniFrac distance were conducted with QIIME (version 2) plugin (Bolyen et al., 2019). Functional annotations of the presented common bacterial taxa were predicted from 16S rRNA gene-based microbial feature compositions using the PICRUSt algorithm (version 2.3.0-b) to make inferences from KEGG database (Douglas et al., 2020).
Identification of the core bacterial microbiome of Karlodinium veneficum across different strains
In this study, bacterial taxa that were commonly and abundantly present across different samples were regarded as the core or stable constituents of the dinoflagellate-associated microbiome, a definition similarly adopted in the previous works on dinoflagellates (Lawson et al., 2018; Sörenson et al., 2019; Maire et al., 2021; Wu et al., 2022). Hence, the core microbiome was defined hereafter as those taxa present in 100% of all samples with a minimum relative abundance of 0.1%. The core bacterial microbiome was delineated by using the “get.coremicrobiome” function of the Mothur software platform with an abundance threshold of 0.1% (Astudillo-Garcia et al., 2017).
Results
General descriptions of pyrosequencing the eukaryotic 28S rRNA and prokaryotic 16S rRNA gene amplicons
The eukaryotic 28S rRNA gene pyrosequencing yielded 1,598,796 raw reads, with an average of 84,147 sequences per sample (Supplementary Table 1). The raw reads were deposited into the NCBI SRA database (BioProject: PRJNA824339). Good coverage value, an indicator of sample completeness (Chao and Jost, 2012), was 1.00 in all 19 samples (Supplementary Table 2), indicating that the sequencing depth sufficed to capture the feature diversity of eukaryotic taxa. Upon removing short sequences, poor quality sequences, and chimeric sequences, 1,546,779 clean sequences remained. After dereplication using DADA2 plugin within the QIIME 2 tool, the 28S rDNA sequencing data set finally generated 729 effective features (sequences clustered at 100% sequence similarity; Supplementary Table 3). Based on the annotations in NCBI database, features assigned to “Karlodinium veneficum” accounted for 97.13–99.93% of the total features in all 19 samples (Supplementary Table 3), manifesting the identity of all cultures used in our study as K. veneficum.
The prokaryotic 16S rRNA gene sequencing generated 1,601,282 raw reads (BioProject: PRJNA771505), with sequences per sample ranging from 80,045 to 87,348 (Supplementary Table 1). Goods coverage of all samples were 1.00 (Supplementary Table 2), indicating sufficient sequences were harvested to reveal the majority of taxa in the prokaryotic assemblages. A total of 1,375,395 clean reads were obtained from the 19 samples (Supplementary Table 1), with a quality control efficiency of 85.89%. The dereplication using DADA2 plugin yielded 1,458 effective features (Supplementary Table 3). Based on the assignments in Silva database, the 1,417 features were used for subsequent analyses, whereas the remained 41 features annotated as “unclassified” were excluded from further analysis.
Feature composition of bacterial communities across strains of Karlodinium veneficum
The 1,417 features identified as bacteria covered 17 phyla, 35 classes, 90 orders, 273 families, and 716 genera. The number of features per sample varied from 74 to 592 (mean = 157; Supplementary Table 3). Generally, Proteobacteria was the absolutely predominant phylum, account for 94.45% of all features, followed by Bacteroidetes (4.30%), Planctomycetes (0.47%), and Cyanobacteria (0.44%; Figure 1A). At the class level, the bacterial community was dominated by γ-proteobacteria (73.50%), α-proteobacteria (17.82%), δ-proteobacteria (3.03%), Bacteroidia (2.87%), and Rhodothermia (0.98%). At the genus level, Alteromonas (22.37%), Marinobacter (20.46%), Methylophaga (19.22%), Alcanivorax (5.74%), Thalassospira (3.52%), UC (unclassified) γ-proteobacteria (3.29%), Ponticoccus (3.22%), Marivita (2.48%), UC δ-proteobacteria (2.27%), UC α-proteobacteria (1.41%), and Roseibacterium (1.03%) were the most dominant genera for all samples as a whole (relative abundance > 1%; Figure 1B).
Figure 1. Relative abundance of bacteria at phylum (A) and genus [top 30; (B)] levels. The abundance is presented in terms of percentage in total effective features in a sample. UC, unclassified.
The PCoA plot was performed to visualize similarity and difference of bacterial composition among different samples at the feature level (sequences clustered at 100% sequence similarity). Notably, the samples did not cluster according to their geographic origins (Figure 2). Six (KV8, KV10, KV12, KV13, KV14, KV15) of the 8 strains from Jiaozhou Bay (Shandong Province, China) and one (KV9) from Ningde (Fujian Province, China) formed cluster A. Two strains (KV1 and KV2) from Ningde (Fujian Province, China) formed cluster C. All the 3 strains from Beibu Gulf (Guangxi Province, China), one strain (KV16) from Pearl River Estuary (Guangdong Province, China), and one member (KV11) from Jiaozhou Bay (Shandong Province, China) formed cluster B (Figure 2). The remained 5 samples exhibited indiscernible affiliation with one another (Figure 2).
Figure 2. Principal coordinate analysis (PCoA) of bacterial community based on weighted-unifrac distances. The color of the symbols indicate the original isolation sources of samples.
Functional predictions of the whole bacterial microbiome of Karlodinium veneficum
To primarily disentangle the putative ecological roles played by the individual bacterial microbiomes in the associated phycospheres, the functional potentials of bacterial consortia were mapped to the reference pathways in the KEGG database. Generally, the whole bacterial microbiome showed high activity in the 3 terms representing major functions of metabolism (47.27%), genetic information processing (19.17%), and environmental information processing (15.16%) at the KEGG level 1, followed by terms of unclassified (7.95%), cellular processes (5.59%), and organismal systems (4.85%; Figure 3A). For the KEGG level 2, the terms membrane transport, amino acid metabolism, and carbohydrate metabolism exhibited relative high abundances to be the top three predominant functional modules (Figure 3B).
Figure 3. Functional prediction of the whole bacterial microbiome associated with Karlodinium veneficum in KEGG categories at level 1 (A) and level 2 (B). Gene functions were predicted from 16S rRNA gene-based microbial compositions using the PICRUSt algorithm to make inferences from KEGG annotated databases. Relative signal intensity was normalized by the number of the genes for each indicated metabolic pathway. a: amino acid metabolism; b: carbohydrate metabolism; c: energy metabolism; d: metabolism of cofactors and vitamins; e: lipid metabolism; f: nucleotide metabolism; g: xenobiotics biodegradation and metabolism; h: glycan biosynthesis and metabolism; m: metabolism of terpenoids and polyketides; n: biosynthesis of other secondary metabolites.
Determination of the core (constitutive) microbiome of Karlodinium veneficum
The 6 genera, Alteromonas, Marinobacter, Methylophaga, Alcanivorax, Thalassospira, and Ponticoccus, were identified as the set of core bacterial microbiome living with K. veneficum, which were referred as “core genera” hereafter. These core genera account for 74.55% of the whole bacterial microbiome, in which Alteromonas was the most predominant, making up 30.01% of the core bacterial microbiome. The next two prevailing genera were Marinobacter (27.45%) and Methylophaga (25.79%; Figure 4). While Alcanivorax made up 5.74%, the genera Thalassospira and Ponticoccus both exceeded 4% in the core microbiome (Figure 4).
Figure 4. Taxonomic composition of core bacterial microbiome associated with Karlodinium veneficum across different strains at the genus level. Half circle diagram represents relative abundances of bacterial genera for the core bacterial microbiome.
Discussion
Bacterial communities of high diversity associated with laboratory-cultured Karlodinium veneficum displayed heterogeneity in features composition among strains from different original sources (geographic origins)
Compared with the relatively high biodiversity observed in terrestrial plant microbiomes (e.g., rhizospheres typically harbor > 103 bacterial taxa; Mendes et al., 2013 and references therein; Ling et al., 2022), some case-studies have reported less diversity in the bacterial communities of microalgae (Krohn-Molt et al., 2017). Especially for works employed laboratory-cultured microalgae, in which presumably well-adapted microbial communities have established during long-term maintenance (Buchan et al., 2014; Seymour et al., 2017). Generally, fewer than 30 bacterial isolates on a species level were affiliated with these microalgae-bacteria biofilms (Krohn-Molt et al., 2017 and the references therein). Here, the associated whole bacterial microbiome of dinoflagellate K. veneficum consisting of 1,417 features, with the number of recovered bacteria at the feature level per sample varying from 37 to 397 (mean = 76), reflecting high diversity and heterogeneity of bacterial flora harbored by laboratory cultures of K. veneficum. Metabolic inferences predicted broad and diverse functional repertoire of the whole bacterial communities in phycospheres, which still required further experimental validation. The PCoA analysis suggested that the geographic origins could influence the composition of the bacterial assemblages to some extent, as evident by the feature composition similarity shared by some strains from the same geographic origin. Yet strains with same and/or close sources could also displayed little or no affiliation with each other, indicating that the original sources was not necessarily the sole and critical determinant for their associated bacterial microbiome structures. We speculated that the excluded algal metabolites and chemical signals as well as the given cultivation conditions might shape and balance the associated bacterial communities of K. veneficum across different strains in our study.
Karlodinium veneficum strains from different geographic origins harbor a core set of constitutive bacterial taxa in their phycospheres, which included three genera commonly reported as dinoflagellate-associated taxa
The whole bacterial microbiome of K. veneficum recovered in our study was generally congruent with the majority of previously characterized bacterial assemblages associated with laboratory cultures of dinoflagellates, in which α- and γ-proteobacteria are the predominant players (Li et al., 2019; Guo et al., 2020; Tarazona-Janampa et al., 2020; Deng et al., 2022; Wu et al., 2022; and the references therein). Despite the heterogeneity of feature composition, bacterial communities of strains from different geographic origins consisted a core set of constitutive taxa. Among them, the Alteromonas and Marinobacter (Alteromonadaceae) are the dominant phylotypes ubiquitously associated with a wide variety of dinoflagellates in both laboratory cultures and field samples (Mayali et al., 2011; Yang et al., 2015; Bolch et al., 2017; Park et al., 2018; Shin et al., 2018; Zhou et al., 2018; Hattenrath-Lehmann et al., 2019; Li et al., 2019; Cui et al., 2020; Guo et al., 2020; Tarazona-Janampa et al., 2020; Kang et al., 2021; Deng et al., 2022; Wu et al., 2022). Metabolically active Marinobacter may be essential for trace metal uptake of dinoflagellates via a beneficial interaction process well known as Iron-Carbon mutualism (Amin et al., 2009). While certain Alteromonas members exhibit algicidal/supportive effect on co-cultured phytoplankton depending on strain and relative abundance. They could inhibit phytoplankton growth by producing algicidal compounds or/and promote growth by removing hydrogen peroxide to relieve oxidative stress on phytoplankton (Aharonovich and Sher, 2016; Hennon et al., 2018). The genus Methylophaga is a unique group of aerobic, halophilic, non-methane-utilizing methylotrophs (Neufeld et al., 2007). Few documents have recorded their associations with dinoflagellates in field samples (Neufeld et al., 2008; Hattenrath-Lehmann et al., 2019), likely due to the carbon sources that methylotrophic bacteria depending on generally present at very low concentrations marine environments. Our recent study has revealed that Methylophaga could be the most significantly enriched bacterial group associated with laboratory-cultured dinoflagellates, relative to that associated with non-dinoflagellates (Deng et al., 2022). In the marine ecosystem, dinoflagellates are the most prolific producers of dimethylsulfoniopropionate (DMSP; Caruana and Malin, 2014). Case study on K. veneficum showed that its intracellular concentrations of DMSP could reach ∼11 mM (Caruana et al., 2012). DMSP could be subsequently degraded and released dimethylsulfide (DMS). The genus Methylophaga included the most majority of marine bacteria knowingly capable of using DMS as a sole carbon source (Boden et al., 2010). Plausibly, the tight association between Methylophaga and K. veneficum observed in our study support such notion.
Highly specialized hydrocarbonoclastic bacteria stably living with Karlodinium veneficum implying that their association with dinoflagellates may be more common in nature than previously perceived
The hydrocarbonoclastic bacteria are prokaryotic hydrocarbon degraders capable of growing in pure culture using hydrocarbons as sole, or major source of carbon for growth (Prince et al., 2010; Thompson et al., 2017). Among them, obligate hydrocarbonoclastic bacteria (OHCB) are a unique group with narrow nutritional spectrum of hydrocarbons almost exclusively as a sole source of carbon and energy (Head et al., 2006; Yakimov et al., 2007). Of the 3 core genera reported here, Alcanivorax, Thalassospira and Ponticoccus, are among those hydrocarbon bacterial utilizers, while two of them (Alcanivorax and Ponticoccus) are recognized as obligate hydrocarbonoclastic taxa (Prince et al., 2010 and the references therein). Alcanivorax strains belonging to the γ-subclass of the Proteobacteria are merely able to grow on a highly restricted spectrum of substrates, predominantly linear and branched n-alkanes and their derivate, with carbon chain length ranging between 9 and 32 carbon molecules (Yakimov et al., 2007). Bacteria affiliated with the genus Thalassospira are marine α-proteobacteria that could degrade several polycyclic aromatic hydrocarbons (PAHs), including naphthalene, dibenzothiophene, phenanthrene and fluorene (Kodama et al., 2008). The γ-proteobacteria Ponticoccus are members that specialize in the PAHs degradation (Gutierrez et al., 2012).
In natural environment, despite the narrow nutritional requirement, hydrocarbonoclastic bacteria are nonetheless cosmopolitan and found throughout the global ocean. Interestingly, their global distribution appears to be solely confined to the marine environments, largely owing to oil-impacted seawater provide sources of carbon and energy to these hydrocarbon-degrading bacteria (see below). Petroleum is a natural, heterogeneous mixture of hydrocarbons, consisting mainly of alkanes with different chain lengths and branch points, cycloalkanes, mono-aromatic and polycyclic aromatic hydrocarbons (McGenity et al., 2012). Hundreds of millions of liters of petroleum enter the marine environment from both natural and anthropogenic sources every year (Head et al., 2006; McGenity et al., 2012). The ecological importance of hydrocarbonoclastic bacteria in the removal of hydrocarbon pollutants was evidenced in the reports documenting their strong enrichment in petrochemical input marine environments (Head et al., 2006; Yakimov et al., 2007 and references therein). In nature, they are usually present at low or undetectable levels before the polluting event and become predominant in petrochemical contamination seawater (Head et al., 2006; Yakimov et al., 2007). Until now, our knowledge about these bacteria (especially OHCB) is largely from few studies on seawater samples collected at petroleum-contaminated sites either by cultivation methods or sequencing surveys (Kodama et al., 2008; McGenity et al., 2012; Thompson et al., 2017). The ecology and physiology of hydrocarbon-degrading bacteria and particularly, their interactions with algae and other microbes during algal blooms remain to be explicitly explored.
Generally, the success of an algae-bacteria relationship has much to do with access, by the bacteria, to an available source of carbon and energy in the form of algal exudates (Bell and Mitchell, 1972) and reciprocal benefits to the phytoplankton through bacteria-mediated bioavailability (McGenity et al., 2012). Here the long-term (from >14 years to several months) and stable co-existence of 3 genera of hydrocarbon-degrading bacteria in the laboratory cultures of K. veneficum could be assumed to at least partly attribute to successful access of hydrocarbonoclastic bacteria to available source of hydrocarbon or hydrocarbon-like substrates (e.g., fatty acids, fatty alcohols, pristine, phytane, and PAHs). Dinoflagellates can be a biogenic source of hydrocarbons (Blumer et al., 1971; McKay et al., 1996; Shaw et al., 2010), and the excluded organic carbon could be subsequently utilized by bacteria such as OHCB. Moreover, either through biogenic synthesis or adsorption, the phytoplankton cell surface could accumulate PAHs (Andelman and Suess, 1970; Gutierrez et al., 2012, 2013; Thompson et al., 2017). Thus, the available hydrocarbon or hydrocarbon-like compounds in the phycosphere could play significant roles in supporting lives of the surrounding hydrocarbon-degrading bacteria, especially in the case for OHCB. However, it is not clear what the exact roles these prokaryotic hydrocarbon degraders played during their association with dinoflagellate hosts. Heretofore, only limited documents have reported the obligate and generalist hydrocarbonoclastic bacteria, either from laboratory cultures of dinoflagellates that have been maintained in cultivation for years, or through sequencing field samples of dinoflagellates freshly collected from places without recognized source of petrochemical input (Gutierrez et al., 2012, 2013; Shin et al., 2018; Guo et al., 2020; Tarazona-Janampa et al., 2020). The bacterial community associated with any non-axenic laboratory culture of phytoplankton is logically a close representative of the community that existed in the phycosperes of the cells in the field from which the culture strain was initially isolated and established (Jasti et al., 2005; Gutierrez et al., 2012, 2013). Hence, the observed stable hydrocarbonoclastic bacteria are very likely to be members of the microbiota associated with the 19 strains of K. veneficum at the time when they were originally established from different water sources. These findings collectively highlighted the occurrence of hydrocarbonoclastic bacteria, especially OHCB, with dinoflagellates may be more frequent in nature than what we have perceived. Our results also raised the possibility that the phycospheres of certain lineages of dinoflagellates (e.g., Kareniaceae) may represent a previously unrecognized niche attractive to bacterial taxa possessing specialized nutritional preference for hydrocarbons. This association brings important questions with respect to the roles that these prokaryotic hydrocarbon degraders played in their relationships with dinoflagellate hosts. Collectively, our results provided a stepwise gain of fundamental insights into the respective and interactive ecology of the key harmful K. veneficum and OHCB. The work also opened a window to further explore the ecological functions exerted by hydrocarbon degraders living with phytoplankton in natural seawaters.
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 in the article/Supplementary material.
Author contributions
YD and KW performed the experiments, provided feedback on the experiments and results, and wrote the article with contributions of all authors. ZH maintained the algal cultures and prepared the samples. QH designed the research. YT designed the research and edited the manuscript. All authors read and approved the final manuscript.
Funding
This research was financially supported by the National Science Foundation of China (Grant No. 42176207), the Key Deployment Project of Centre for Ocean Mega-Research of Science, Chinese Academy of Sciences (Grant No. COMS2019Q09), and the Natural Science Foundation of Shenzhen, China (Grant No. JCYJ20210324094013037).
Acknowledgments
We sincerely thank Fengting Li for her help with part of the experiments.
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.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.967610/full#supplementary-material
References
Aharonovich, D., and Sher, D. (2016). Transcriptional response of Prochlorococcus to co-culture with a marine Alteromonas: Differences between strains and the involvement of putative infochemicals. ISME J. 10, 2892–2906. doi: 10.1038/ismej.2016.70
Amin, S. A., Green, D. H., Hart, M. C., Küpper, F. C., Sunda, W. G., and Carrano, C. J. (2009). Photolysis of iron-siderophore chelates promotes bacterial-algal mutualism. Proc. Natl. Acad. Sci. U.S.A. 106, 17071–17076. doi: 10.1073/pnas.0905512106
Amin, S. A., Parker, M. S., and Armbrust, E. V. (2012). Interactions between diatoms and bacteria. Microbiol. Mol. Biol. Rev. 76, 667–684. doi: 10.1128/MMBR.00007-12
Andelman, J. B., and Suess, M. J. (1970). Polynuclear aromatic hydrocarbons in the water environment. Bull. World Health Organ. 43, 479–508.
Astudillo-Garcia, C., Bell, J. J., Webster, N. S., Glasl, B., Jompa, J., Montoya, J. M., et al. (2017). Evaluating the core microbiota in complex communities: A systematic investigation. Environ. Microbiol. 19, 1450–1462. doi: 10.1111/1462-2920.13647
Behringer, G., Ochsenkühn, M. A., Fei, C., Fanning, J., Koester, J. A., and Amin, S. A. (2018). Bacterial communities of diatoms display strong conservation across strains and time. Front. Microbiol. 9:659. doi: 10.3389/fmicb.2018.00659
Bell, W., and Mitchell, R. (1972). Chemotactic and growth responses of marine bacteria to algal extracellular products. Biol. Bull. 143, 265–277. doi: 10.2307/1540052
Blumer, M., Guillard, R. R. L., and Chase, T. (1971). Hydrocarbons of marine phytoplankton. Mar. Biol. 8, 183–189. doi: 10.1007/BF00355214
Boden, R., Kelly, D. P., Murrell, J. C., and Schäfer, H. (2010). Oxidation of dimethylsulfide to tetrathionate by Methylophaga thiooxidans sp. nov.: A new link in the sulfur cycle. Environ. Microbiol. 12, 2688–2699. doi: 10.1111/j.1462-2920.2010.02238.x
Bolch, C. J. S., Bejoy, T. A., and Green, D. H. (2017). Bacterial associates modify growth dynamics of the dinoflagellate Gymnodinium catenatum. Front. Microbiol. 8:670. doi: 10.3389/fmicb.2017.00670
Bolyen, E., Rideout, J. R., Dillon, M. R., Bokulich, N. A., Abnet, C. C., Al-Ghalith, G. A., et al. (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37, 852–857. doi: 10.1038/s41587-019-0209-9
Buchan, A., LeCleir, G., Gulvik, C., and González, J. M. (2014). Master recyclers: Features and functions of bacteria associated with phytoplankton blooms. Nat. Rev. Microbiol. 12, 686–698. doi: 10.1038/nrmicro3326
Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869
Caruana, A., Steinke, M., Turner, S., and Malin, G. (2012). Concentrations of dimethylsulphoniopropionate and activities of dimethylsulphide-producing enzymes in batch cultures of nine dinoflagellate species. Biogeochemistry 110, 87–107. doi: 10.1007/s10533-012-9705-4
Caruana, A. M., and Malin, G. (2014). The variability in DMSP content and DMSP lyase activity in marine dinoflagellates. Prog. Oceanogr. 120, 410–424. doi: 10.1016/j.pocean.2013.10.014
Chai, Z. Y., He, Z. L., Deng, Y. Y., Yang, Y. F., and Tang, Y. Z. (2018). Cultivation of seaweed Gracilaria lemaneiformis enhanced biodiversity in an eukaryotic plankton community as revealed via metagenomic analyses. Mol. Ecol. 27, 1081–1093. doi: 10.1111/mec.14496
Chao, A., and Jost, L. (2012). Coverage-based rarefaction and extrapolation: Standardizing samples by completeness rather than size. Ecology 93, 2533–2547. doi: 10.1890/11-1952.1
Cui, Y., Chun, S. J., Baek, S. S., Baek, S. H., Kim, P. J., Son, M., et al. (2020). Unique microbial module regulates the harmful algal bloom (Cochlodinium polykrikoides) and shifts the microbial community along the Southern Coast of Korea. Sci. Total Environ. 721:137725. doi: 10.1016/j.scitotenv.2020.137725
Deeds, J. R., Terlizzi, D. E., Adolf, J. E., Stoecker, D. K., and Place, A. R. (2002). Toxic activity from cultures of Karlodinium micrum (=Gyrodinium galatheanum) (Dinophyceae)-a dinoflagellate associated with fish mortalities in an estuarine aquaculture facility. Harmful Algae 1, 169–189. doi: 10.1016/S1568-9883(02)00027-6
Deng, Y. Y., Wang, K., Hu, Z. X., and Tang, Y. Z. (2022). Abundant species diversity and essential functions of bacterial communities associated with dinoflagellates as revealed from metabarcoding sequencing for laboratory-raised clonal cultures. Int. J. Environ. Res. Public Health 19:4446. doi: 10.3390/ijerph19084446
Douglas, G. M., Maffei, V. J., Zaneveld, J. R., Yurgel, S. N., Brown, J. R., Taylor, C. M., et al. (2020). PICRUSt2 for prediction of metagenome functions. Nat. Biotechnol. 38, 685–688. doi: 10.1038/s41587-020-0548-6
Durán, P., Flores-Uribe, J., Wippel, K., Zhang, P., Guan, R., Melkonian, M., et al. (2022). Shared features and reciprocal complementation of the Chlamydomonas and Arabidopsis microbiota. Nat. Commun. 13:406. doi: 10.1038/s41467-022-28055-8
Edwards, J., Johnson, C., Santos-Medellin, C., Lurie, E., Podishetty, N. K., Bhatnagar, S., et al. (2015). Structure, variation, and assembly of the root-associated microbiomes of rice. Proc. Natl. Acad. Sci. U.S.A. 112, E911–E920. doi: 10.1073/pnas.1414592112
Finnish Environment Institute (2016). Last summer’s fish kill was caused by a toxic dinoflagellate: Emerging algal toxins in coastal Finnish waters. Available Online at: https://www.ymparisto.fi/en-US/Sea/Last_summers_fish_kill_was_caused_by_a_t(38306) [accessed July 13, 2022].
Grady, K. L., Sorensen, J. W., Stopnisek, N., Guittar, J., Shade, A., et al. (2019). Assembly and seasonality of core phyllosphere microbiota on perennial biofuel crops. Nat. Commun. 10:4135. doi: 10.1038/s41467-019-11974-4
Grossart, H., Levold, F., Allgaier, M., Simon, M., and Brinkhoff, T. (2005). Marine diatom species harbour distinct bacterial communities. Environ. Microbiol. 7, 860–873. doi: 10.1111/j.1462-2920.2005.00759.x
Guillard, R. R. L. (1975). “Culture of phytoplankton for feeding marine invertebrates,” in Culture of marine invertebrate animals, eds W. L. Smith and M. H. Chanley (New York, NY: Plenum Press). doi: 10.1007/978-1-4615-8714-9_3
Guo, R. Y., Wang, P. B., Lu, D. D., and Dai, X. F. (2020). Comparison of bacterial communities associated with Prorocentrum donghaiense and Karenia mikimotoi strains from Chinese coastal waters. Mar. Freshw. Res. 71, 1662–1671. doi: 10.1071/MF20035
Gutierrez, T., Green, D. H., Nichols, P. D., Whitman, W. B., Semple, K. T., and Aitken, M. D. (2013). Polycyclovorans algicola gen. nov., sp. nov., an aromatic-hydrocarbon-degrading marine bacterium found associated with laboratory cultures of marine phytoplankton. Appl. Environ. Microbiol. 79, 205–214. doi: 10.1128/AEM.02833-12
Gutierrez, T., Nichols, P. D., Whitman, W. B., and Aitken, M. D. (2012). Porticoccus hydrocarbonoclasticus sp. nov., an aromatic hydrocarbon-degrading bacterium identified in laboratory cultures of marine phytoplankton. Appl. Environ. Microbiol. 78, 628–637. doi: 10.1128/AEM.06398-11
Hattenrath-Lehmann, T. K., Jankowiak, J., Koch, F., and Gobler, C. J. (2019). Prokaryotic and eukaryotic microbiomes associated with blooms of the ichthyotoxic dinoflagellate Cochlodinium=(Margalefidinium) polykrikoides in New York, USA, estuaries. PLoS One 14:e0223067. doi: 10.1371/journal.pone.0223067
Head, I. M., Jones, D. M., and Röling, W. F. (2006). Marine microorganisms make a meal of oil. Nat. Rev. Microbiol. 4, 173–182. doi: 10.1038/nrmicro1348
Hennon, G. M., Morris, J. J., Haley, S. T., Zinser, E. R., Durrant, A. R., Entwistle, E., et al. (2018). The impact of elevated CO2 on Prochlorococcus and microbial interactions with ‘helper’ bacterium Alteromonas. ISME J. 12, 520–531. doi: 10.1038/ismej.2017.189
Jasti, S., Sieracki, M. E., Poulton, N. J., Giewat, M. W., and Rooney-Varga, J. N. (2005). Phylogenetic diversity and specificity of bacteria closely associated with Alexandrium spp. and other phytoplankton. Appl. Environ. Microbiol. 71, 3483–3494. doi: 10.1128/AEM.71.7.3483-3494.2005
Jeong, H. J., Kang, H. C., Lim, A. S., Jang, S. H., Lee, K., Lee, S. Y., et al. (2021). Feeding diverse prey as an excellent strategy of mixotrophic dinoflagellates for global dominance. Sci. Adv. 7:eabe4214. doi: 10.1126/sciadv.abe4214
Kang, J., Park, J. S., Jung, S. W., Kim, H. J., Joo, H. M., Kang, D., et al. (2021). Zooming on dynamics of marine microbial communities in the phycosphere of Akashiwo sanguinea (Dinophyta) blooms. Mol. Ecol. 30, 207–221. doi: 10.1111/mec.15714
Kodama, Y., Stiknowati, L., Ueki, A., Ueki, K., and Watanabe, K. (2008). Thalassospira tepidiphila sp. nov., a polycyclic aromatic hydrocarbon degrading bacterium isolated from seawater. Int. J. Syst. Evol. Microbiol. 58, 711–715. doi: 10.1099/ijs.0.65476-0
Kouzuma, A., and Watanabe, K. (2015). Exploring the potential of algae/bacteria interactions. Curr. Opin. Biotechnol. 33, 125–129. doi: 10.1016/j.copbio.2015.02.007
Krohn-Molt, I., Alawi, M., Förstner, K. U., Wiegandt, A., Burkhardt, L., Indenbirken, D., et al. (2017). Insights into microalga and bacteria interactions of selected phycosphere biofilms using metagenomic, transcriptomic, and proteomic approaches. Front. Microbiol. 8:1941. doi: 10.3389/fmicb.2017.01941
Lawson, C. A., Raina, J. B., Kahlke, T., Seymour, J. R., and Suggett, D. J. (2018). Defining the core microbiome of the symbiotic dinoflagellate, Symbiodinium. Environ. Microbiol. Rep. 10, 7–11. doi: 10.1111/1758-2229.12599
Li, S. F., Chen, M. C., Chen, Y. F., Tong, J., Wang, L. Y., Xu, Y., et al. (2019). Epibiotic bacterial community composition in red-tide dinoflagellate Akashiwo sanguinea culture under various growth conditions. FEMS Microbiol. Ecol. 95:fiz057. doi: 10.1093/femsec/fiz057
Lim, H. C., Leaw, C. P., Tan, T. H., Kon, N. F., Yek, L. H., Hii, K. S., et al. (2014). A bloom of Karlodinium australe (Gymnodiniales, Dinophyceae) associated with mass mortality of cage-cultured fishes in West Johor Strait, Malaysia. Harmful Algae 40, 51–62. doi: 10.1016/j.hal.2014.10.005
Ling, N., Wang, T., and Kuzyakov, Y. (2022). Rhizosphere bacteriome structure and functions. Nat. Commun. 13:836. doi: 10.1038/s41467-022-28448-9
Logue, J. B., Stedmon, C. A., Kellerman, A. M., Nielsen, N. J., Andersson, A. F., Laudon, H., et al. (2016). Experimental insights into the importance of aquatic bacterial community composition to the degradation of dissolved organic matter. ISME J. 10, 533–545. doi: 10.1038/ismej.2015.131
Lundberg, D. S., Lebeis, S. L., Paredes, S. H., Yourstone, S., Gehring, J., Malfatti, S., et al. (2012). Defining the core Arabidopsis thaliana root microbiome. Nature 488, 86–90. doi: 10.1038/nature11237
Lupette, J., Lami, R., Krasovec, M., Grimsley, N., Moreau, H., Piganeau, G., et al. (2016). Marinobacter dominates the bacterial community of the Ostreococcus tauri phycosphere in culture. Front. Microbiol. 7:1414. doi: 10.3389/fmicb.2016.01414
Magoè, T., and Salzberg, S. L. (2011). FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507
Maire, J., Girvan, S. K., Barkla, S. E., Perez-Gonzalez, A., Suggett, D. J., Blackall, L. L., et al. (2021). Intracellular bacteria are common and taxonomically diverse in cultured and in hospite algal endosymbionts of coral reefs. ISME J. 15, 2028–2042. doi: 10.1038/s41396-021-00902-4
Maryland Department of the Environment (2017). Karlodinium veneficum and fish kills what you need to know. Available Online at: https://mde.maryland.gov/programs/Water/FishandShellfish/Documents/karlofactsheetfinal032017.pdf [accessed July 13, 2022].
Mayali, X., Franks, P. J., and Burton, R. S. (2011). Temporal attachment dynamics by distinct bacterial taxa during a dinoflagellate bloom. Aquat. Microb. Ecol. 63, 111–122. doi: 10.3354/ame01483
McGenity, T. J., Folwell, B. D., McKew, B. A., and Sanni, G. O. (2012). Marine crude-oil biodegradation: A central role for interspecies interactions. Aquat. Biosyst. 8:10. doi: 10.1186/2046-9063-8-10
McKay, W. A., Turner, M. F., Jones, B. M. R., and Halliwell, C. M. (1996). Emissions of hydrocarbons from marine phytoplankton-some results from controlled laboratory experiments. Atmos. Environ. 30, 2583–2593. doi: 10.1016/1352-2310(95)00433-5
Mendes, L. W., Kuramae, E. E., Navarrete, A. A., van Veen, J. A., and Tsai, S. M. (2014). Taxonomical and functional microbial community selection in soybean rhizosphere. ISME J. 8, 1577–1587. doi: 10.1038/ismej.2014.17
Mendes, R., Garbeva, P., and Raaijmakers, J. M. (2013). The rhizosphere microbiome: Significance of plant beneficial, plant pathogenic, and human pathogenic microorganisms. FEMS Microbiol. Rev. 37, 634–663. doi: 10.1111/1574-6976.12028
Neu, A. T., Allen, E. E., and Roy, K. (2021). Defining and quantifying the core microbiome: Challenges and prospects. Proc. Natl. Acad. Sci. U.S.A. 118:e2104429118. doi: 10.1073/pnas.2104429118
Neufeld, J. D., Boden, R., Moussard, H., Schäfer, H., and Murrell, J. C. (2008). Substrate-specific clades of active marine methylotrophs associated with a phytoplankton bloom in a temperate coastal environment. Appl. Environ. Microbol. 74, 7321–7328. doi: 10.1128/AEM.01266-08
Neufeld, J. D., Schäfer, H., Cox, M. J., Boden, R., McDonald, I. R., and Murrell, J. C. (2007). Stable-isotope probing implicates Methylophaga spp. and novel Gammaproteobacteria in marine methanol and methylamine metabolism. ISME J. 1, 480–491. doi: 10.1038/ismej.2007.65
Park, B. S., Guo, R., Lim, W. A., and Ki, J. S. (2018). Importance of free living and particle-associated bacteria for the growth of the harmful dinoflagellate Prorocentrum minimum: Evidence in culture stages. Mar. Freshw. Res. 69, 290–299. doi: 10.1071/MF17102
Peiffer, J. A., Spor, A., Koren, O., Jin, Z., Tringe, S. G., Dangl, J. L., et al. (2013). Diversity and heritability of the maize rhizosphere microbiome under field conditions. Proc. Natl. Acad. Sci. U.S.A. 110, 6548–6553. doi: 10.1073/pnas.1302837110
Prince, R. C., Gramain, A., and McGenity, T. J. (2010). “Prokaryotic hydrocarbon degraders,” in Handbook of Hydrocarbon and Lipid Microbiology, ed. K. N. Timmis (Berlin: Springer-Verlag), 1669–1692. doi: 10.1007/978-3-540-77587-4_118
Risely, A. (2020). Applying the core microbiome to understand host-microbe systems. J. Anim. Ecol. 89, 1549–1558. doi: 10.1111/1365-2656.13229
Rognes, T., Flouri, T., Nichols, B., Quince, C., and Mahé, F. (2016). VSEARCH: A versatile open source tool for metagenomics. Peer J. 4:e2584. doi: 10.7717/peerj.2584
Seymour, J. R., Amin, S. A., Raina, J. B., and Stocker, R. (2017). Zooming in on the phycosphere: The ecological interface for phytoplankton-bacteria relationships. Nat. Microbiol. 2, 1–12. doi: 10.1038/nmicrobiol.2017.65
Shade, A., and Handelsman, J. (2012). Beyond the Venn diagram: The hunt for a core microbiome. Environ. Microbiol. 14, 4–12. doi: 10.1111/j.1462-2920.2011.02585.x
Shaw, S. L., Gantt, B., and Meskhidze, N. (2010). Production and emissions of marine isoprene and monoterpenes: A review. Adv. Meteorol. 24:408696. doi: 10.1155/2010/408696
Shin, H., Lee, E., Shin, J., Ko, S. R., Oh, H. S., Ahn, C. Y., et al. (2018). Elucidation of the bacterial communities associated with the harmful microalgae Alexandrium tamarense and Cochlodinium polykrikoides using nanopore sequencing. Sci. Rep. 8:5323. doi: 10.1038/s41598-018-23634-6
Sörenson, E., Bertos-Fortis, M., Farnelid, H., Kremp, A., Krüger, K., Lindehoff, E., et al. (2019). Consistency in microbiomes in cultures of Alexandrium species isolated from brackish and marine waters. Environ. Microbiol. Rep. 11, 425–433. doi: 10.1111/1758-2229.12736
Tarazona-Janampa, U. I., Cembella, A. D., Pelayo-Zárate, M. C., Pajares, S., Márquez-Valdelamar, L. M., Okolodkov, Y. B., et al. (2020). Associated bacteria and their effects on growth and toxigenicity of the dinoflagellate Prorocentrum lima species complex from epibenthic substrates along Mexican coasts. Front. Marine Sci. 7:569. doi: 10.3389/fmars.2020.00569
Thompson, H., Angelova, A., Bowler, B., Jones, M., and Gutierrez, T. (2017). Enhanced crude oil biodegradative potential of natural phytoplankton-associated hydrocarbonoclastic bacteria. Environ. Microbiol. 19, 2843–2861. doi: 10.1111/1462-2920.13811
Wu, Z., Lee, W. H., Liu, Z., Lin, S., and Lam, P. K. S. (2022). Microbiome associated with Gambierdiscus balechii cultures under different toxicity conditions. Front. Mar. Sci. 9:760553. doi: 10.3389/fmars.2022.760553
Xu, N., Pang, S. J., and Liu, F. (2012). Molecular identification of a bloom-forming species isolated from Sanggou Bay in Shandong Province. Mar. Sci. 36, 13–18.
Yakimov, M. M., Timmis, K. N., and Golyshin, P. N. (2007). Obligate oil-degrading marine bacteria. Curr. Opin. Biotechnol. 18, 257–266. doi: 10.1016/j.copbio.2007.04.006
Yang, C., Yi, L., Zhou, B., Zhou, Y., Wei, Z., Yun, T., et al. (2015). Illumina sequencing-based analysis of free-living bacterial community dynamics during an Akashiwo sanguinea bloom in Xiamen sea, China. Sci. Rep. 5:8476. doi: 10.1038/srep08476
Yang, H. J., Hu, Z. X., Shang, L. X., Deng, Y. Y., and Tang, Y. Z. (2020). A strain of the toxic dinoflagellate Karlodinium veneficum isolated from the East China Sea is an omnivorous phagotroph. Harmful Algae 93:101775. doi: 10.1016/j.hal.2020.101775
Yang, H. J., Hu, Z. X., and Tang, Y. Z. (2021). Plasticity and multiplicity of trophic modes in the dinoflagellate Karlodinium and their pertinence to population maintenance and bloom dynamics. J. Mar. Sci. Eng. 9:51. doi: 10.3390/jmse9010051
Yang, H. J., Hu, Z. X., Xu, N., and Tang, Y. Z. (2019). A comparative study on the allelopathy and toxicity of four strains of Karlodinium veneficum with different culturing histories. J. Plankton Res. 41, 17–29. doi: 10.1093/plankt/fby047
Keywords: Alcanivorax, core bacterial microbiome, hydrocarbon-degrading bacteria, Karlodinium veneficum, Ponticoccus, Thalassospira
Citation: Deng Y, Wang K, Hu Z, Hu Q and Tang YZ (2022) Identification and implications of a core bacterial microbiome in 19 clonal cultures laboratory-reared for months to years of the cosmopolitan dinoflagellate Karlodinium veneficum. Front. Microbiol. 13:967610. doi: 10.3389/fmicb.2022.967610
Received: 14 June 2022; Accepted: 18 July 2022;
Published: 04 August 2022.
Edited by:
Jin Zhou, Tsinghua University, ChinaReviewed by:
Yang Chen, Hebei Normal University of Science and Technology, ChinaKyoko Yarimizu, National Research Institute of Fisheries Science, Japan
Copyright © 2022 Deng, Wang, Hu, Hu and Tang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ying Zhong Tang, yingzhong.tang@qdio.ac.cn
†These authors have contributed equally to this work and share first authorship