Abstract
Pyropia is an economically important edible red alga worldwide. The aquaculture industry and Pyropia production have grown considerably in recent decades. Microbial communities inhabit the algal surface and produce a variety of compounds that can influence host adaptation. Previous studies on the Pyropia microbiome were focused on the microbial components or the function of specific microbial lineages, which frequently exclude metabolic information and contained only a small fraction of the overall community. Here, we performed a genome-centric analysis to study the metabolic potential of the Pyropia haitanensis phycosphere bacteria. We reconstructed 202 unique metagenome-assembled genomes (MAGs) comprising all major taxa present within the P. haitanensis microbiome. The addition of MAGs to the genome tree containing all publicly available Pyropia-associated microorganisms increased the phylogenetic diversity by 50% within the bacteria. Metabolic reconstruction of the MAGs showed functional redundancy across taxa for pathways including nitrate reduction, taurine metabolism, organophosphorus, and 1-aminocyclopropane-1-carboxylate degradation, auxin, and vitamin B12 synthesis. Some microbial functions, such as auxin and vitamin B12 synthesis, that were previously assigned to a few Pyropia-associated microorganisms were distributed across the diverse epiphytic taxa. Other metabolic pathways, such as ammonia oxidation, denitrification, and sulfide oxidation, were confined to specific keystone taxa.
Introduction
Marine algae are an ancient and functionally important component of aquatic ecosystems (Raven et al., 2002; ; ), providing valuable sources of food and habitat for a variety of marine microorganisms and animals (Seymour et al., 2017). Similar to the rhizosphere, microbes inhabit the “phycosphere” () surrounding seaweed-based cells, which extend from the cell surface through the dispersal boundary layer of seaweed-derived dissolved organic compounds, secondary metabolites, and exopolymeric substances (). Interactions between algae and microbes are thought to occur in this specific environment (). In recent years, many phycologists and microbiologists have explored potential interactions between algae and microbes (Wichard, 2015; ; ). The currently accepted hypothesis is that microorganisms have beneficial, neutral, and detrimental effects on algae, with beneficial effects being the most widely studied (; van Tol et al., 2017). Algae-associated bacteria are involved in the exchange of diverse chemical currencies, including vitamins (; Xie et al., 2013), hormones (), quorum sensing signals (Geng and Belas, 2010; Huang et al., 2018), dissolved organic carbon (Teeling et al., 2012), and nutrient remineralization (Klawonn et al., 2019). These studies have mainly been conducted on some ecologically important phytoplankton (Geng and Belas, 2010; Xie et al., 2013; ; ; Huang et al., 2018; Foster and Zehr, 2019; Nowinski et al., 2019; Rambo et al., 2020) and macroalgae species (Twigg et al., 2014; Vollmers et al., 2017; Roth-Schulze et al., 2018), and only few studies on farmed seaweed, such as Pyropia (Wang et al., 2021).
Pyropia is an edible red alga (phylum: Rhodophyta; class: Bangiophyceae), which is mainly consumed as processed food products (). Because Pyropia contains high levels of vitamins (e.g., vitamin B12), minerals (e.g., iron), and protein contents (), it is usually regarded as a health-promoting food. In East Asian countries such as China, Japan, and South Korea, Pyropia is widely cultivated () and its production was over 28.7 million tons in 2018 according to FAO statistics (Food and Agriculture Organization [FAO], 2020). Under economic incentives, intensive cultivation breeding and practices increase cultivar susceptibility to diseases and induce colonization and infection of harmful microbes (Yan et al., 2019), leading to the reduced production of Pyropia. In contrast, some microorganisms promote algal growth () and morphological development. Previous studies have shown that marine macroalgae, including Ulva (Provasoli and Pintner, 1980; Nakanishi et al., 1996; Marshall et al., 2006), Monostroma (Matsuo et al., 2003), and Pyropia (Fukui et al., 2014), cannot develop normal morphology under axenic conditions.
Recent efforts continue to explore potential Pyropia-microbial interactions in different environments and statuses. For instance, analysis of the epiphytic microorganisms of P. haitanensis using amplicon sequencing revealed that the microbial communities differed significantly during various life stages of P. haitanensis and hypothesized that microbial taxa could produce some plant growth regulators to promote the growth of P. haitanensis (Wang et al., 2020). Moreover, a study on the microbial community of P. yezoensis infected with red rot disease revealed the interactions between the disease and the epiphytic and planktonic microbial communities, and the potential of using community differentiation to forecast disease occurrence (Yan et al., 2019). The above studies have made important contributions to the community structure and dynamics of microorganisms in Pyropia and have also shown a complex relationship between the host and microbiome. However, these features can only be assigned at taxon level and do not provide information on the functional characteristics of a community. Our previous study on the thalli of laboratory-cultured P. haitanensis (PH40) revealed microbial gene functions (Wang et al., 2021). Gene-centric analysis of the whole community has provided critical functional insights into the Pyropia microbiome, but could not link function with phylogeny. In this study, we performed a genome-centric analysis to characterize 202 MAGs from the phycosphere of P. haitanensis cultivars in China. Our sequenced data represented more than 95% of the entire microbiome. These MAGs could comprise approximately 65% of the sequenced microbes, providing an unprecedented opportunity to determine the metabolic potential of all major taxa within a laver microbiome.
Materials and Methods
Sample Collection and DNA Extraction
P. haitanensis samples were collected from farming regions in southern and northern China in November 2019. Two farms near coastal water were located in Rizhao City (35°42′N, 119°50′E), Shandong Province, China; and Ningde City (26°85′N, 120°03′E), Fujian Province, China. Three sampling points were designated at each farm; each sampling point was approximately half nautical mile apart. A total of 18 samples were obtained. Each sample contained approximately 100 P. haitanensis thalli. The thalli were placed in sterile bags containing seawater from the same sampling location and transported to the laboratory within 8 h in a low-temperature (4°C) box. Parts of the thalli were cleaned to extract metagenomic DNA, and the rest were stored in a freezer at −80°C.
To remove loosely associated microbes and other attachments, 20 thalli from each sample were randomly selected and washed three times with sterilized seawater (). Then, the clean thalli were cut into fragments of approximately 3 mm × 3 mm using sterilized blades and cleaned again with sterilized seawater. These thalli were used as material for metagenomic DNA extraction. All steps were performed on a clean bench. Metagenomic DNA was extracted using the DNeasy PowerSoil Kit (QIAGEN, Germany), according to the manufacturer’s instructions (Yan et al., 2019; ). The concentration of the extracted DNA was determined using the NanoDrop ND-2000 spectrophotometer (NanoDrop, Wilmington, DE, United States), and the mass was evaluated using a 1% agarose gel. To minimize extraction bias, the three replicates from each sampling point were pooled, resulting in six DNA samples that were denoted as RZ-A, RZ-B, RZ-C, ND-A, ND-B, and ND-C.
Library Preparation, Sequence Preprocessing, and Assembly
Construction of a metagenomic DNA library and sequencing of paired ends (2 × 150 bp) were performed by Shanghai Persenor according to the standard protocol.1 Each sample contained approximately 30 GB of raw reads. The 3′ and 5′ end adapters were removed using Trim Galore v0.5.0,2 and reads with length < 20 bp were removed. Paired reads were mapped to the P. haitanensis genome (GCA_008729055.1) using Bowtie2 v2.3.2.3 To minimize the impact of host contamination, we filtered out read sequences aligned to the P. haitanensis genome. The PCR replicates were removed using FastUniq v1.1 (Xu et al., 2012). The contamination rate (%) of clean reads was assessed by remapping the reads to the host genome. Non-pareil v3.30 was used to assess the average coverage of each metagenome (Rodriguez-R et al., 2018). The clean reads of each sample were assembled using SPAdes v3.13.0 (metaspades model, kmers were 21, 33, and 55) (). Then, we merged all of the unassembled reads and co-assembled them. Finally, all contigs were moved to a single file for subsequent analysis. The assembly results were statistically analyzed using QUAST v4.5 (Gurevich et al., 2013). Contigs with lengths of less than 1,000 bp were filtered. Finally, to check the assembly effect and quality of the data, clean reads were realigned to contigs.
Taxonomic Assignment and Metagenomic Binning
To observe the microbial community, Kaiju v1.5.0 (Menzel et al., 2016) was used for taxonomic assignment based on the clean reads of each sample. Non-redundant NCBI BLAST was selected as the taxonomic annotation database. The rarefaction curves for each metagenome were plotted using R language v4.0.1 [vegan package, rarecurve(), step = 2,000] based on the Kaiju results.
Binning was performed using MetaBAT2 v2.12.1 (Kang et al., 2019), MaxBin2 v2.2.4 (-max_iteration 5) (Wu et al., 2016), and CONCOCT (). The “bin_refinement” function of MetaWRAP v1.1.8 (Uritskiy et al., 2018) was used to optimize the draft genomes (complete > 50%, contamination < 10%). The dRep v2.5.4 (Olm et al., 2017) was used to remove the redundant MAGs.
Using the Genome Taxonomy Database R89 (GTDB) (Parks et al., 2018), taxonomy was assigned to MAGs using GTDB-Tk v1.1.0 (). The “pipe” and “appraise” functions of SingleM v0.13.24 were used to scan for single copy marker genes in both the MAGs and the sequenced raw reads to estimate the proportion of marker genes recovered in the MAGs. The “genome” function of coverM v0.5.05 was used to evaluate the relative abundance of MAGs in each sample based on read alignments. MAGs with a cumulative abundance > 1% (sum of each MAG abundance in each sample) were included in the heatmap visualized using R language (packages including pheatmap and RColorBrewer), as shown in Supplementary Figure 1.
Phylogenetic Tree Building and Metabolic Reconstruction of the Metagenome-Assembled Genomes
We filtered out MAGs with quality < 40% (calculated as “Quality = Genome completeness (%) - [5 × Contamination (%)]”). A concatenated marker gene tree was inferred using the 190 MAGs after quality control, together with 231 available IMG/M6Pyropia isolated genomes (by 29 October 2020) that were ≥ 50% complete and contained ≤ 10% contamination. For comparison, taxonomy was reassigned to the previously Pyropia-associated genomes using the GTDB-R89 database (Supplementary Table 1). Thereafter, the bacterial tree was inferred using GTDB-Tk and embellished using ITOL (Letunic and Bork, 2021). As no archaeal genome was isolated and submitted to public databases before this study, phylogenetic distance (PD) and phylogenetic gain (PG) for the bacterial tree were calculated using GenomeTreeTk v0.0.417 to determine the degree to which the phylogenetic diversity from the current study was added to the tree.
The annotation and annotation functions of EnrichM v0.5.08 were used to annotate MAGs with the Kyoto Encyclopedia of Genes and Genome (KEGG) Orthology database (Kanehisa et al., 2021). Genes that were not annotated were assigned as a “hypothetical protein”. Considering that we obtained incomplete microbial genomes, to assign a pathway/module to MAG, the appropriate KEGG pathway/module needed to be ≥ 60% complete, and all key enzymes must be present (as defined in the literature) (Engelberts et al., 2020). Because the KEGG modules had no auxin biosynthesis module and only a partial cobalamin synthesis module, we created specific modules (Supplementary Note 1). The results of the feature annotation are shown in the phylogenetic tree with KO numbers (Supplementary Figures 2–7).
Results
Microbial Community Structures of the Phycosphere
After the quality control step, 209 GB of data was obtained. Approximately 100% of the reads had a base quality greater than 99% (quality score > 20). The rate of clean reads mapped to the host genome was less than 0.1%.
On average, 47% of the reads were assigned to the taxonomic database. A total of 174 phyla were detected, in which five phyla with relative abundances of over 1% (Figure 1A). A total of 3,433 genera were identified, including 15 genera with relative abundances of more than 1% (Figure 1B). The rarefaction curves (Figure 1C) and non-pareil curves (Figure 1D) indicated that the sequencing data obtained in this study covered more than 95% of the whole microbiome.
FIGURE 1
Dereplication of the MAGs reduced the total number to 202 (75.05 ± 24.95% completeness, 4.74 ± 4.74% contamination), including 1 archaeal MAG-3 belonging to the phylum Thaumarchaeota (Crenarchaeota as per GTDB) and 201 bacterial MAGs belonging to eight phyla, including Actinobacteriota (12 MAGs), Bacteroidota (100 MAGs), Bdellovibrionota (5 MAGs), Chloroflexota (6 MAGs), Myxococcota (4 MAGs), Patescibacteria (15 MAGs), Proteobacteria (54 MAGs), and Verrucomicrobiota (5 MAGs) (Supplementary Table 2). All dominant taxa with a relative abundance of > 1% were detected in the MAGs. The MAGs contained an average of 3,466 ± 2,747 genes, of which 30.38 ± 20.24% were hypothetical proteins (Supplementary Table 2). Moreover, the analysis of single-copy marker genes presented in the microbial community and MAGs using SingleM showed that 202 MAGs represented 65.2% of the sequenced microbiome. The above results suggested that we obtained a representative dataset of the Pyropia microbiome and were able to obtain the metabolic potential of all dominant taxa.
A concatenated genome tree was constructed from 231 Porphyra/Pyropia-isolated microbial genomes from the IMG/M database and 190 MAGs after quality control in this study (Figure 2). Although previous genomes were widely distributed across the phylogenetic tree, the 190 MAGs from this study still provided a phylogenetic gain (the additional branch length contributed by a set of taxa) of 50% for the bacteria (Table 1).
FIGURE 2
TABLE 1
| No. taxa | PD | Percent PD (%) | |
| Full tree | 421 | 75.71 | 100.00 |
| Outgroup taxa (PD)a | 231 | 37.54 | 49.58 |
| Ingroup taxa (PD)b | 190 | 51.07 | 67.46 |
| Ingroup taxa (PG) | 190 | 38.18 | 50.42 |
Calculated phylogenetic distance (PD) and phylogenetic gain (PG) for the bacterial genome trees, per taxon.
aOutgroup taxa included previous Pyropia symbiont genomes in the IMG/M database and bingroup taxa included MAGs retrieved from this study.
Autotrophic Carbon Fixation Pathways in the Pyropia haitanensis Microbiome
We searched the MAGs for the six known prokaryotic carbon fixation modules in map00720 (Kanehisa et al., 2021), including the Wood-Ljungdahl (WL) pathway (M00377), reductive citric acid (rTCA) cycle (M00173), 3-hydroxypropionate/4-hydroxybutyrate (HP-HB) cycle (M00375), 3-hydroxypropionate (3-HP) cycle (M00376), dicarboxylate/4-hydroxybutyrate (DC-HB) cycle (M00374), and the Calvin Benson Bassham (CBB) cycle (M00165) (Supplementary Table 3). Only the CBB cycle was observed in MAG-99 and MAG-191 (family Rhodobacteraceae) (Supplementary Figure 2). Although the rTCA and the 3-HP cycles had more than 60% module integrity in some MAGs, both lacked the key enzymes ATP-citrate lyases (aclAB) and malonyl-CoA reductase (mcr), respectively (Zarzycki et al., 2009; Hugler and Sievert, 2011). The other three carbon fixation modules (HP-HB, DC-HB, and WL) were < 60% complete in the MAGs.
Nitrogen-Transforming Network in the Pyropia haitanensis Microbiome
Nitrate can be reduced to ammonium via assimilatory nitrate reduction (ANRA) and dissimilatory nitrate reduction (DNRA) (). ANRA includes assimilatory nitrate reductases (nasA or narB) and ferredoxin-nitrite reductase (nirA), while DNRA includes nitrate reductases (narGHI or napAB) and nitrite reductases (nirBD or nrfAH) (Fiore et al., 2010). The nasA and narB were identified in Acidimicrobiia (3 MAGs), Alphaproteobacteria (12 MAGs), Gammaproteobacteria (8 MAGs), Anaerolineae (1 MAG), and Bacteroidia (9 MAGs) (Figure 3). Meanwhile, nirA was identified in the family DEV007 (2 MAGs). For DNRA, narGHI or napAB were identified in the genera Halocynthiibacter (MAG-90), Arenitalea (MAG-201), and Colwellia (MAG-91). Four phyla, Proteobacteria (20 MAGs), Bacteroidetes (8 MAGs), Actinobacteria (3 MAGs), and Chloroflexota (1 MAG), contained nirBD. On the other hand, denitrification comprises four steps: nitrate reduction (narGHI or napAB), nitrite reduction (nirK or nirS), nitric oxide reduction (norBC), and nitrous oxide reduction (nosZ) (Fiore et al., 2010; Figure 4). The nirK was identified in MAG-3 (family Nitrosopumilaceae) and MAG-201. Furthermore, MAG-201 had norBC and nosZ, whereas nosZ was also identified in MAG-91. The nxrAB encode nitrite oxidoreductase, which catalyzes the nitrite oxidation (Sohaskey and Wayne, 2003). Both genes were identified in MAG-90. Except for the two genes (amoBC) of ammonia oxidation identified in MAG-3, KEGG Orthologys (Kos) associated with nitrification (Kanehisa et al., 2021) were absent from other MAGs, as were the KOs involved in nitrogen fixation (Supplementary Figure 3).
FIGURE 3
FIGURE 4
Dissolved Organic Phosphorus Hydrolyzation and Organosulfur Metabolism
One of the most important mechanisms that allow marine organisms to cope with low orthophosphate (Pi) stress is the utilization of Dissolved Organic Phosphorus (DOP) by alkaline phosphatase (AP, phoABD) (Lin et al., 2016). The AP genes were identified in 136 MAGs, including Actinobacteria (8 MAGs), Bacteroidetes (86 MAGs), Bdellovibrionota (5 MAGs), Chloroflexota (2 MAGs), Chloroflexota (1 MAG), Myxococcota (1 MAG), Patescibacteria (3 MAGs), Proteobacteria (26 MAGs), and Verrucomicrobiota (5 MAGs) (Figure 3 and Supplementary Figure 4).
Microorganisms can metabolize the host taurine via ATP-binding cassette (ABC) transporters (tauACB), and taurine is subsequently converted to sulfite via taurine dioxygenases (tauD) (Rambo et al., 2020; Figure 5). The tauD was identified in MAGs belonging to Actinobacteria (8 MAGs), Bacteroidota (1 MAG), Chloroflexota (2 MAGs), Alphaproteobacteria (9 MAGs), and Gammaproteobacteria (4 MAGs), of which tauACB transporters were also identified in Actinobacteriota (6 MAGs), Chloroflexota (2 MAGs), Alphaproteobacteria (5 MAGs), and Gammaproteobacteria (3 MAGs) (Figure 3). In addition to the cleavage of taurine, sulfite may also originate from sulfide by dissimilatory sulfite reductase (dsrAB) (Engelberts et al., 2020; Figure 5). Both genes were identified in MAG-68 (Gammaproteobacteria) (Supplementary Figure 5). In turn, sulfite can be oxidized to sulfate via sulfite dehydrogenase/oxidase (sorAB/SUOX/soeABC) (Rambo et al., 2020). The SUOX/soeABC were identified in all phyla except Bdellovibrionota, Patescibacteria, Verrucomicrobiota, and Crenarchaeota (Figure 3). The assimilatory sulfate reduction module (M00176), which can reduce sulfate to sulfide (Kanehisa et al., 2021), was found in the phyla Actinobacteria (4 MAGs), Bacteroidetes (47 MAGs), Chloroflexota (5 MAGs), Thaumarchaeota (1 MAG), Myxococcota (3 MAGs), Proteobacteria (9 MAGs), and Verrucomicrobiota (4 MAGs) (Figure 3). The sulfur oxidation system (SOX) is a pathway that includes soxAB and soxXYZ (Rambo et al., 2020). These genes were found in five Alphaproteobacteria and six Gammaproteobacteria MAGs.
FIGURE 5
Metabolic Potential for the Synthesis of Cobalamin
In the present study, potential cobalamin producers in MAGs were identified in different lineages, including the family Cellvibrionaceae (1 MAG), Flavobacteriaceae (1 MAG), Rhizobiaceae (2 MAGs), Rhodobacteraceae (13 MAGs), SZUA-35 (1 MAG), UBA10066 (1 MAG), and UBA6668 (3 MAGs) (Figure 3). Except for MAG-116, MAG-184, and MAG-48, 19 MAGs had the potential to synthesize cobalamin via both aerobic and anaerobic routes (Supplementary Figure 6). In addition to de novo biosynthesis, cobalamin may also originate from the lower axial ligand conversion of pseudocobalamin (Hoffmann et al., 2000). The bluB, which encodes for 5,6-dimethylbenzimidazole synthase (DMB), was present in the families Rhodobacteraceae (6 MAGs), Rhizobiaceae (2 MAGs), Nitrosopumilaceae (1 MAG), and Flavobacteriaceae (1 MAG). Amidohydrolase (CbiZ), and cobinamide-phosphate synthase (CbiB) can remodel pseudocobalamin into cobalamin by replacing adenine with DMB (Yi et al., 2012). These two genes were only identified in MAG-179 (family UBA11606).
The btuB, which encodes for TonB-dependent outer membrane cobalamin receptor and transporter (Noinaj et al., 2010), was identified in the phyla Bacteroidetes (85 MAGs), Bdellovibrionota (3 MAGs), Proteobacteria (19 MAGs), Myxococcota (2 MAGs), and Verrucomicrobiota (4 MAGs) (Supplementary Figure 6). However, few cobalamin biosynthesis genes were found in these MAGs. Moreover, 157 MAGs lacked metE (B12-independent methionine synthase), while 115 MAGs only included metH (B12-dependent methionine synthase) (Supplementary Table 3).
Potential of Metagenome-Assembled Genomes to Produce Auxin and Influence Ethylene Biosynthesis
We searched all the MAGs for the four known 3-indoleacetic acid (IAA) biosynthesis pathways: the indole-3-acetamide (IAM), indole-3-pyruvate (IPyA), tryptamine (TAM), and indole-3-acetonitrile (IAN) pathways (Supplementary Note 1). The metabolic potential of the IAM pathway was present in the families Granulosicoccaceae (3 MAGs) and UBA11606 (2 MAGs) (Figure 3). The TAM pathway was present in 104 MAGs, including the phyla Actinobacteria (8 MAGs), Bacteroidetes (47 MAGs), Bdellovibrionota (1 MAG), Chloroflexota (4 MAGs), Myxococcota (2 MAGs), Patescibacteria (8 MAGs), Proteobacteria (33 MAGs), and Verrucomicrobiota (1 MAG) (Figure 3). However, for the IPyA pathway, the gene encoding the key enzyme (indolepyruvate decarboxylase, ipdC) (Duca et al., 2014) was absent in all MAGs (Supplementary Table 3). Furthermore, although some MAGs (7 MAGs in order Chitinophagales) were contained YUCCA (encoding the homologs of indole-3-pyruvate monooxygenase) (Mano and Nemoto, 2012), while Tam1 and IL4I1 (encoding aminotransferases and L-amino-acid oxidase, respectively) were absent (Supplementary Table 3). The enzymes encoded by Tam1 and IL4I1 can transform tryptophan into IPyA (Duca et al., 2014). More than 60% of the genes necessary for carrying out the IAN pathway were found in 11 MAGs; however, indoleacetaldoxime dehydratase (CYP71A13) and myrosinase (E3.2.1.147) were absent (Supplementary Figure 7). These two genes are related to the conversion of indole-3-acetaldoxime (IAOX) and glucobrassicin to IAN (Zhang et al., 2019).
Microorganisms can promote plant growth indirectly by converting the ethylene precursor 1-aminocyclopropane-1-carboxylate (ACC) to α-ketobutyrate and ammonia by ACC deaminase (acdS) (). The acdS was identified in the phyla Actinobacteria (4 MAGs), Bacteroidetes (77 MAGs), Proteobacteria (7 MAGs), and Chloroflexota (1 MAG) (Figure 3), of which 46 MAGs are also potential IAA producers. MAG-144 may contain both TAM and IAM pathways (Supplementary Figure 7).
Discussion
Heterotrophic Microorganisms Dominate Pyropia haitanensis Microbiota
The products of photosynthesis by plants and algae, such as rhamnose, xylose, and glucose, can provide a carbon source for epiphytic microorganisms (; Rigonato et al., 2012; Thapa and Prasanna, 2018). Microorganisms that fix carbon make epiphytic microorganisms less dependent on photosynthetic products and allow them to set a niche for the colonization of other microbes (Rigonato et al., 2012; Thapa and Prasanna, 2018; Flores-Nunez et al., 2020). However, it is still unclear whether all epiphytic microbes are dependent on the carbon provided by Pyropia. Here, only two members of Rhodobacteraceae (MAG-99 and MAG-191) fix carbon using the CBB cycle (Supplementary Figure 2), suggesting that epiphytic microorganisms were largely heterotrophic need exogenous carbon to support their life activities. In addition, host-acquired autotrophic carbon from microbial fixation could also contribute to overall high environmental tolerance, particularly to ocean warming (). However, the relative contribution and mechanisms of microbial carbon fixation to marine algal nutrition should be studied further.
Metabolic Potential of Nutrients
Except for Cyanobacteria, plant cannot fix nitrogen on their own, and the process of nitrogen fixation usually requires the assistance of plant symbiotic microorganisms, such as rhizobia (; Kuypers et al., 2018). In our study, no genes associated with nitrogen fixation were found in MAGs, suggesting that the P. haitanensis epiphytic microbiomes are unlikely to catalyze the conversion of airborne N2 to NH3. Within the nitrate reduction reaction, microorganisms can convert nitrate and nitrite to ammonium through ANRA and DNRA (). A recent study showed the presence of ammonium transporters in P. yezoensis (Kakinuma et al., 2017), implying that microbial ammonium may also be trade with the host. Excess nitrite has been shown to have toxic effects on organisms, such as affecting the proton permeability of cell membranes (Vadivelu et al., 2006; Zhou et al., 2010) and inhibiting photosynthesis (Singh et al., 2007). The microbial consumption of nitrite also reduces the risk of excessive nitrite (Vaucheret et al., 1992; Duncanson et al., 1993). Nitrate or nitrite can also be converted to nitrogen oxides (N2O and NO) via denitrification (Fiore et al., 2010). In addition to being released into the environment, NO is an effective signaling molecule in plant-rhizobacteria interactions, promoting plant growth and development in an auxin-dependent manner (Lamattina et al., 2003; Molina-Favero et al., 2008). N2O is a potent greenhouse gas and ozone-depleting substance, and N2O emissions can cause atmospheric pollution (Ravishankara et al., 2009). Hence, MAG-91 and MAG-201 may also contribute to the reduction of N2O emissions to the environment.
Phosphorus is a central component of nucleic acids and phospholipids and plays a central role in the production of chemical energy (Forlani et al., 2008). However, algae grown in coastal regions may experience Pi starvation because of the rapid consumption of Pi for photosynthesis (Lin et al., 2016). Moreover, anthropogenic N also increases the N:P ratio in receiving waters (Zhang et al., 2007; Zhou et al., 2017), thus limiting the Pi of surrounding organisms. Although there are other forms of phosphorus sources in the ocean, organisms preferentially utilize Pi as they can be taken up directly (Huang et al., 2005; Lin et al., 2016). The major mechanism by which marine phytoplankton and bacteria can convert DOP to bioavailable Pi is the induction of alkaline phosphatase (phoABD) (Suzumura et al., 2012; Lin et al., 2019). Here, phoABD were identified in 136 MAGs (Supplementary Figure 4). DOP hydrolysis has costly energetic metabolic activity because of the stable C-P bonds (Lin et al., 2016). In terrestrial plants, while both plant and microbial phosphatases are efficient in releasing Pi from soil organic phosphorus, there is evidence that microbial enzymes show a higher efficiency of phosphorus release (Tarafdar et al., 2001). Therefore, it is possible that the “extra” microbial Pi can be assimilated to support host metabolism and growth (Jones, 1972; ; ; Hoppe, 2003).
Taurine, dietary fibers, polyunsaturated fatty acids, and sulfated polysaccharides are enriched in Pyropia (; ). Our results showed the potential of 16 MAGs to obtain foreign taurine via ABC transporters and convert it to sulfite via dioxygenases (Supplementary Figure 5). Sulfate produced during dissimilatory sulfate oxidation and dissimilatory thiosulfate oxidation may provide the raw material for assimilatory sulfate reduction for cysteine (Takahashi et al., 2011) and other cell material biosynthesis. Similar to the nitrogen cycle, genes required for certain reactions in the sulfur cycle have been identified only in specific microbial lineages, implies that they may play prominent roles in nitrogen/sulfur cycling within the Pyropia phycosphere.
Cobalamin and Phytohormone Metabolism
In many organisms with metH, vitamin B12 usually acts as an enzyme cofactor participating in catalyzing primary biochemical reactions (Sanudo-Wilhelmy et al., 2014), such as amino acid () and DNA synthesis (). Currently, cobalamin is thought to be synthesized by a relatively small set of prokaryotes, and algae may obtain cobalamin from these associated microbes (Giovannoni, 2012; Lynch and Neufeld, 2015). Our previous study revealed six potential cobalamin producers, all of which belong to Alphaproteobacteria (Wang et al., 2021). Here, 22 MAGs in five phyla were considered potential producers of cobalamin through aerobic or anaerobic pathway (Supplementary Figure 6). This result suggested that vitamin B12 producers were not restricted to a particular microbial lineage. Another source of cobalamin is the substitution of the adenine lower axial ligand of pseudocobalamin by DMB (Helliwell et al., 2016). Pseudocobalamin is a vitamer that can be synthesized in cyanobacteria (Hoffmann et al., 2000). Furthermore, amidohydrolase (cbiZ) and cobinamide-phosphate synthase (cbiB) can replace the adenine lower axial ligand with DMB (Yi et al., 2012). In this study, 10 MAGs were identified to have DMB synthesis genes, and one MAG (MAG-179) was identified to have remodeled pseudocobalamin genes. These MAGs may play a key role in maximizing the effect of cobalamin production (Lu et al., 2020). In terms of potential cobalamin consumption, we found 113 MAGs with the TonB-dependent cobalamin transporter gene (btuB) and 115 MAGs with the B12-dependent methionine synthase gene (metH). The MAGs with btuB had low complete of the cobalamin synthesis pathway (Supplementary Figure 6). This phenomenon was also presented in soil cobalamin research (Lu et al., 2020); our study was more specific by presenting at the individual level. Both btuB and metH were identified in more than half of the MAGs, suggesting that cobalamin producers play important roles in the microbial community (; Lu et al., 2020).
Auxin is one of the important plant hormones that influence plant growth and development (Vanneste and Friml, 2009; ). In the rhizosphere, the potential to synthesize auxin is considered the main characteristic of plant growth-promoting rhizobacteria (PGPR) (Spaepen et al., 2007). A previous study showed that co-culture with IAA producers can positively affect the growth of P. yezoensis (Matsuda et al., 2018). Furthermore, bacterial IAA is also reported to regulate microbial physiology, including overcoming stress, regulating cellular processes, and providing a carbon and nitrogen source (; Tromas and Perrot-Rechenmann, 2010; Van Puyvelde et al., 2011). Here, the IAM and TAM pathways were present in 109 MAGs, including diverse microbial lineages (Supplementary Figure 7). Both pathways are also widely distributed in rhizobacteria (Zhang et al., 2019). The IPyA pathway was usually discovered in plant growth-promoting rhizobacteria (Fu et al., 2012); however, the key gene (ipdC) was absent in the MAGs (Supplementary Figure 7). Notably, although the IPyA and IAN pathways lack genes encoding key enzymes, they both have a high degree of integrity in some MAGs (Supplementary Table 3). Therefore, it does not exclude the possibility that bacteria can synthesize IAA from intermediates of other microbial lineages (Zhang et al., 2019), even of host ().
Ethylene is associated with plant senescence, abscission, flower development, and pathogen defense signaling (Lin et al., 2009; Van de Poel et al., 2015). Ethylene biosynthesis is induced under diverse stress conditions, thereby inhibiting plant growth (). In this process, ACC is an important intermediate (Lin et al., 2009). A previous report showed that exogenous addition of ACC can increase ethylene production and promote the formation of spermatia and zygotospores in P. yezoensis (Uji et al., 2020). In the rhizosphere, microbial ACC deaminase can metabolize ACC exuded by plants, resulting in a reduction or stalling of ethylene biosynthesis inside the host owing to a decrease in the precursor (Fu et al., 2012), and it is also another feature of PGPR (Wang et al., 2007; ). In this study, the 89 MAGs possessed the ACC deaminase gene, of which 46 MAGs also possessed IAA synthesis potential (Supplementary Figure 7). These MAGs with “dual identity” may be advantageous for future screening of beneficial bacteria (Lugtenberg and Kamilova, 2009).
Metabolic Potential of Archaea
Although archaea are extremely widespread in the biosphere (Vaucheret et al., 1992), their potential role in macroalgae is rarely mentioned (Trias et al., 2012; Karray et al., 2017; Ghosh et al., 2019), especially in Pyropia. In this study, we obtained an archaeal draft genome (MAG-3) belonging to the Nitrosopumilaceae family. Our result filled a gap in the understanding of the archaeal genome in the epiphytic microorganisms of Pyropia. Further functional analysis of MAG-3 revealed its potential role in MAG-3 ammonia oxidation (Supplementary Figure 3), nitrite reduction, DMB biosynthesis, and assimilatory sulfate reduction (Supplementary Figures 5, 6). However, there may be some other functional features that have not been identified owing to the incomplete genome.
Different Metabolic Potentials Between Two Sampling Sites
We also evaluated the differences between microbial communities though the abundance of MAGs (Supplementary Figure 9). Regarding inorganic salt metabolism, the relative abundance of MAGs with DOP metabolic potential was higher in the ND group (RZ vs. ND: 56.6% vs. 78.31%); in contrast, the relative abundance of MAGs with sulfur and nitrogen metabolic potential was higher in the RZ group (S: 61.58% vs. 52.93%; N: 17.16% vs. 15.52%). In terms of phytohormone anabolism, the relative abundance of MAGs with IAA synthesis potential exceeded 60% in both groups (RZ vs. ND: 77.58% vs. 61.76%); the relative abundance of MAGs with metabolic ACC potential was higher in the ND group (35.78% vs. 53.05%). Both environment and algae are capable of influencing the community structure of epiphytic microorganisms, so we speculate that the reason for these differences may be the result of a combination of environmental and algal effects (Yan et al., 2019; Wang et al., 2020; ).
It should be emphasized that our results and conclusions were based on bioinformatics methods. Although we recovered 65.2% of the sequencing community, the remaining 35.8% of unbinned sequences may contain functions that we have not characterized. In addition, bioinformatics-based predictions may include incomplete metabolic pathways because of incomplete genome analysis. We believe that future studies aimed at understanding Pyropia as a holobiont should also move toward the retrieval of host genes and genomes so that the metabolic potential of Pyropia and its epiphytic microorganisms can be assessed for metabolic complementarity.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, PRJNA726403.
Author contributions
JW, ZM, and YM designed the experiments and prepared the manuscript. JW and XT collected the materials. JW did the experiments and bioinformatics analysis. All authors were involved in revision of the manuscript and approved its final version.
Funding
This work was supported by the National Natural Science Foundation of China (Grant No. 32060829), the Fundamental Research Funds for the Central Universities (202064006), the National Key R&D Program of China (Grant No. 2020YFD0901101), the 2020 Research Program of Sanya Yazhou Bay Science and Technology City (Grant No. SKJC-2020-02-009), Special Project of Central Government Guiding Local Science and Technology Development (Grant No. ZY2020HN02), and the Innovation Platform for Academicians of Hainan Province.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.857901/full#supplementary-material
Supplementary Figure 1Relative abundance of the dominant MAGs (>1% abundance across the data). “Unmapped” represents reads that could not be mapped to any of the 202 MAGs. RZ, Rizhao; ND, Ningde.
Supplementary Figure 2KO absence/presence in the MAGs retrieved from Pyropia haitanensis that are involved in six autotrophic carbon fixation pathways found in marine prokaryotes, from left to right: Wood-Ljungdahl pathway, reductive citric acid cycle, hydroxypropionate-hydroxybutyrate cycle, 3-hydroxypropionate bi-cycle, dicarboxylate-hydroxyburyrate cycle, and the Calvin-Benson-Bassham cycle. Corresponding KEGG modules are displayed on the left of the figure. Similar colored squares indicate KOs involved in the same step of the module as defined by KEGG. Blue colored lineages have metabolic potential for the Calvin-Benson-Bassham cycle. Branch labels display taxonomy at the lowest inferred level. Bacterial tree was rooted to GCA_002453875.1. The KEGG module was retrieved from https://www.genome.jp/on December 7, 2020.
Supplementary Figure 3KO absence/presence in the MAGs retrieved from Pyropia haitanensis that are involved in nitrogen fixation, nitrification, denitrification, assimilatory nitrate reduction, and dissimilatory nitrate reduction. Corresponding KEGG modules are displayed on the left of the figure. Similar colored squares indicate KOs involved in the same step of the module as defined by KEGG. Branch labels display taxonomy at the lowest inferred level. Bacterial tree was rooted to GCA_002453875.1. KEGG modules were retrieved from https://www.genome.jp/on December 7, 2020.
Supplementary Figure 4KO absence/presence in the MAGs retrieved from Pyropia haitanensis that are involved in taurine transport, sulfide and sulfite oxidation, and sulfate reduction. Corresponding KEGG modules are displayed on the left of the figure. Similar colored squares indicate KOs involved in the same step of the module as defined by KEGG. Branch labels display taxonomy at the lowest inferred level. Bacterial tree was rooted to GCA_002453875.1. KEGG modules were retrieved from https://www.genome.jp/on December 7, 2020.
Supplementary Figure 5KO absence/presence in the MAGs retrieved from Pyropia haitanensis that are involved in hydrolyzing dissolved organic phosphorus (DOP). Similar colored squares indicate KOs involved in the same step of the module as defined by KEGG. Blue colored lineages have the metabolic potential to hydrolyze DOP. Branch labels display taxonomy at the lowest inferred level. Bacterial tree was rooted to GCA_002453875.1.
Supplementary Figure 6KO absence/presence in MAGs retrieved from Pyropia haitanensis that are involved in indole-3-acetic acid synthesis and ACC deaminase (K01505). Corresponding KEGG modules are displayed on the left of the figure. Similar colored squares indicate KOs involved in the same step of the module as defined by KEGG. Blue colored lineages have the metabolic potential to produce IAA by the TAM pathway. Bold lineages have the metabolic potential to produce IAA by the IAM pathway. Italic lineages have the metabolic potential to produce ACC deaminase. Branch labels display taxonomy at the lowest inferred level. Bacterial tree was rooted to GCA_002453875.1. KEGG modules were retrieved from https://www.genome.jp/on December 7, 2020.
Supplementary Figure 7KO absence/presence in MAGs retrieved from Pyropia haitanensis that are involved in cobalamin (vitamin B12) metabolism. Blue colored lineages can metabolize cobalamin through either one of the two pathways. Bold lineages have the metabolic potential to produce 5,6-dimethylbenzimidazole (DMB). Red colored lineages can remodel cobalamin-like compounds (such as pseudocobalamin) into cobalamin. Branch labels display taxonomy at the lowest inferred level. Bacterial tree was rooted to GCA_002453875.1.
Supplementary Figure 8The ongrowing technique and morphology of Pyropia haitanensis. The fixed nets are hung between poles (A,B). (C) Is morphology of Pyropia haitanensis.
Supplementary Figure 9Evaluation of the differences between microbial communities. The X-axis shows the relative abundance of MAGs with metabolic potential. ACC, 1-aminocyclopropane-1-carboxylate; DOP, dissolved organic phosphorus; IAM, Indole-3-acetamide pathway; N, Nitrogen metabolism; S, Sulfur metabolism; TAM, Tryptamine pathway; RZ, Rizhao; ND, Ningde.
Supplementary Table 1Taxonomy and statistics of the IMG/M database with Pyropia symbiont genomes included in the concatenated marker gene tree. The GTDB R89 taxonomy is presented.
Supplementary Table 2Taxonomy and statistics of 202 MAGs retrieved from this study. The GTDB R89 taxonomy is presented.
Supplementary Table 3KO count for each of the 202 MAGs retrieved in this study. The GTDB R89 taxonomy is presented. KOs are involved in the described metabolic pathways (i.e., carbon fixation, nitrogen, sulfur, phytohormone, and vitamin B12).
Supplementary Note 1Auxin biosynthesis modules and cobalamin synthesis modules.
Footnotes
2.^https://github.com/FelixKrueger/TrimGalore
3.^https://github.com/BenLangmead/bowtie2
4.^https://github.com/wwood/singlem
5.^https://github.com/wwood/CoverM
6.^https://img.jgi.doe.gov/cgi-bin/m/main.cgi
References
1
AbelesF. B.MorganP. W.SaltveitM. E.Jr. (2012). Ethylene in Plant Biology.London: Academic Press.
2
AhmedA.KhurshidA.TangX. H.WangJ. H.KhanT. U.MaoY. X. (2021). Structural and functional impacts of microbiota on pyropia yezoensis and surrounding seawater in cultivation farms along coastal areas of the yellow sea.Microorganisms9:1291. 10.3390/microorganisms9061291
3
AlnebergJ.BjarnasonB. S.de BruijnI.SchirmerM.QuickJ.IjazU. Z.et al (2014). Binning metagenomic contigs by coverage and composition.Nat. Methods111144–1146. 10.1038/Nmeth.3103
4
AminS. A.HmeloL. R.van TolH. M.DurhamB. P.CarlsonL. T.HealK. R.et al (2015). Interaction and signalling between a cosmopolitan phytoplankton and associated bacteria.Nature52298–101. 10.1038/nature14488
5
AntunesJ.LeaoP.VasconcelosV. (2019). Marine biofilms: diversity of communities and of chemical cues.Environ. Microbiol. Rep.11287–305. 10.1111/1758-2229.12694
6
ArandaM.LiY.LiewY. J.BaumgartenS.SimakovO.WilsonM. C.et al (2016). Genomes of coral dinoflagellate symbionts highlight evolutionary adaptations conducive to a symbiotic lifestyle.Sci. Rep.6:39734. 10.1038/srep39734
7
ArmbrustE. V.BergesJ. A.BowlerC.GreenB. R.MartinezD.PutnamN. H.et al (2004). The genome of the diatom thalassiosira pseudonana: ecology, evolution, and metabolism.Science30679–86. 10.1126/science.1101156
8
BanerjeeR. V.MatthewsR. G. (1990). Cobalamin-dependent methionine synthase.FASEB J.41450–1459.
9
BankevichA.NurkS.AntipovD.GurevichA. A.DvorkinM.KulikovA. S.et al (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing.J. Comput. Biol.19455–477. 10.1089/cmb.2012.0021
10
BellJ. J.BennettH. M.RovelliniA.WebsterN. S. (2018). Sponges to be winners under near-future climate scenarios.Bioscience68955–968. 10.1093/biosci/biy142
11
BellW.MitchellR. (1972). Chemotactic and growth responses of marine bacteria to algal extracellular products.Biol. Bull.143265–277.
12
BiancoC.ImperliniE.CalogeroR.SenatoreB.AmoresanoA.CarpentieriA.et al (2006). Indole-3-acetic acid improves escherichia coli’s defences to stress.Arch. Microbiol.185373–382. 10.1007/s00203-006-0103-y
13
BidonB.KabbaraS.CourdavaultV.GlevarecG.OudinA.HericourtF.et al (2020). Cytokinin and ethylene cell signaling pathways from prokaryotes to eukaryotes.Cells9:2526. 10.3390/cells9112526
14
BlakleyR. L.BarkerH. A. (1964). Cobamide stimulation of the reduction of ribotides to deoxyribotides in lactobacillus leichmannii.Biochem. Biophys. Res. Commun.16391–397. 10.1016/0006-291x(64)90363-8
15
BlouinN. A.BrodieJ. A.GrossmanA. C.XuP.BrawleyS. H. (2011). Porphyra: a marine crop shaped by stress.Trends Plant Sci.1629–37. 10.1016/j.tplants.2010.10.004
16
BrodieJ.ChanC. X.De ClerckO.CockJ. M.CoelhoS. M.GachonC.et al (2017b). The algal revolution.Trends Plant Sci.22726–738. 10.1016/j.tplants.2017.05.005
17
BrodieJ.BallS. G.BougetF. Y.ChanC. X.De ClerckO.CockJ. M.et al (2017a). Biotic interactions as drivers of algal origin and evolution.New Phytol.216670–681. 10.1111/nph.14760
18
BrucknerC. G.RehmC.GrossartH. P.KrothP. G. (2011). Growth and release of extracellular organic compounds by benthic diatoms depend on interactions with bacteria.Environ. Microbiol.131052–1063. 10.1111/j.1462-2920.2010.02411.x
19
BuchanA.LeCleirG. R.GulvikC. A.GonzalezJ. M. (2014). Master recyclers: features and functions of bacteria associated with phytoplankton blooms.Nat. Rev. Microbiol.12686–698. 10.1038/nrmicro3326
20
BulgarelliD.SchlaeppiK.SpaepenS.van ThemaatE. V. L.Schulze-LefertP. (2013). Structure and functions of the bacterial microbiota of plants.Annu. Rev. Plant Biol.64807–838. 10.1146/annurev-arplant-050312-120106
21
BurkeC.KjellebergS.ThomasT. (2009). Selective extraction of bacterial DNA from the surfaces of macroalgae.Appl. Environ. Microbiol.75252–256. 10.1128/Aem.01630-08
22
CaoJ.WangJ.WangS.XuX. (2016). Porphyra species: a mini-review of its pharmacological and nutritional properties.J. Med. Food19111–119. 10.1089/jmf.2015.3426
23
ChaumeilP. A.MussigA. J.HugenholtzP.ParksD. H. (2020). GTDB-Tk: a toolkit to classify genomes with the genome taxonomy database.Bioinformatics361925–1927. 10.1093/bioinformatics/btz848
24
ChengQ. (2008). Perspectives in biological nitrogen fixation research.J. Integrat. Plant Biol.50786–798. 10.1111/j.1744-7909.2008.00700.x
25
ChoT. J.RheeM. S. (2020). Health functionality and quality control of laver (porphyra, pyropia): current issues and future perspectives as an edible seaweed.Mar. Drugs18:14. 10.3390/md18010014
26
ColeJ. J. (1982). Interactions between bacteria and algae in aquatic ecosystems.Annu. Rev. Ecol. Syst.13291–314. 10.1146/annurev.es.13.110182.001451
27
CroftM. T.LawrenceA. D.Raux-DeeryE.WarrenM. J.SmithA. G. (2005). Algae acquire vitamin B12 through a symbiotic relationship with bacteria.Nature43890–93. 10.1038/nature04056
28
CurrieD. J. (1990). Large-scale variability and interactions among phytoplankton, bacterioplankton, and phosphorus.Limnol. Oceanogr.351437–1455. 10.4319/lo.1990.35.7.1437
29
del Carmen Orozco-MosquedaM.GlickB. R.SantoyoG. (2020). ACC deaminase in plant growth-promoting bacteria (PGPB): an efficient mechanism to counter salt stress in crops.Microbiol. Res.235:126439. 10.1016/j.micres.2020.126439
30
DevolA. H. (2015). Denitrification, anammox, and N(2) production in marine sediments.Ann. Rev. Mar. Sci.7403–423. 10.1146/annurev-marine-010213-135040
31
DucaD.LorvJ.PattenC. L.RoseD.GlickB. R. (2014). Indole-3-acetic acid in plant–microbe interactions.Antonie Van Leeuwenhoek10685–125. 10.1007/s10482-013-0095-y
32
DuncansonE.GilkesA. F.KirkD. W.ShermanA.WrayJ. L. (1993). Nir1, a conditional-lethal mutation in barley causing a defect in nitrite reduction.Mol. Gen. Genet.236275–282. 10.1007/Bf00277123
33
EngelbertsJ. P.RobbinsS. J.de GoeijJ. M.ArandaM.BellS. C.WebsterN. S. (2020). Characterization of a sponge microbiome using an integrative genome-centric approach.ISME J.141100–1110. 10.1038/s41396-020-0591-9
34
FioreC. L.JarettJ. K.OlsonN. D.LesserM. P. (2010). Nitrogen fixation and nitrogen transformations in marine symbioses.Trends Microbiol.18455–463. 10.1016/j.tim.2010.07.001
35
Flores-NunezV. M.Fonseca-GarciaC.DesgarennesD.Eloe-FadroshE.WoykeT.Partida-MartinezL. P. (2020). Functional signatures of the epiphytic prokaryotic microbiome of agaves and cacti.Front. Microbiol.10:3044. 10.3389/fmicb.2019.03044
36
Food and Agriculture Organization [FAO] (2020). FAO Yearbook of Fishery and Aquaculture Statistics.Rome: FAO.
37
ForlaniG.PavanM.GramekM.KafarskiP.LipokJ. (2008). Biochemical bases for a widespread tolerance of cyanobacteria to the phosphonate herbicide glyphosate.Plant Cell Physiol.49443–456. 10.1093/pcp/pcn021
38
FosterR. A.ZehrJ. P. (2019). Diversity, genomics, and distribution of phytoplankton- cyanobacterium single-cell symbiotic associations.Annu. Rev. Microbiol.73435–456. 10.1146/annurev-micro-090817-062650
39
FuM. Z.WangZ. L.PuX. M.XuZ. J.ZhuM. Y. (2012). Changes of nutrient concentrations and N:P:Si ratios and their possible impacts on the huanghai sea ecosystem.Acta Oceanol. Sin.31101–112. 10.1007/s13131-012-0224-x
40
FukuiY.AbeM.KobayashiM.YanoY.SatomiM. (2014). Isolation of hyphomonas strains that induce normal morphogenesis in protoplasts of the marine red alga pyropia yezoensis.Microb. Ecol.68556–566. 10.1007/s00248-014-0423-4
41
GengH. F.BelasR. (2010). Molecular mechanisms underlying roseobacter-phytoplankton symbioses.Curr. Opin. Biotechnol.21332–338. 10.1016/j.copbio.2010.03.013
42
GhoshS.GnaimR.GreisermanS.FadeevL.GozinM.GolbergA. (2019). Macroalgal biomass subcritical hydrolysates for the production of polyhydroxyalkanoate (PHA) by Haloferax mediterranei.Bioresour. Technol.271166–173. 10.1016/j.biortech.2018.09.108
43
GiovannoniS. J. (2012). Vitamins in the sea.Proc. Natl. Acad. Sci. U.S.A.10913888–13889. 10.1073/pnas.1211722109
44
GurevichA.SavelievV.VyahhiN.TeslerG. (2013). QUAST: quality assessment tool for genome assemblies.Bioinformatics291072–1075. 10.1093/bioinformatics/btt086
45
HelliwellK. E.LawrenceA. D.HolzerA.KudahlU. J.SassoS.KräutlerB.et al (2016). Cyanobacteria and eukaryotic algae use different chemical variants of vitamin B12.Curr. Biol.26999–1008. 10.1016/j.cub.2016.02.041
46
HoffmannB.OberhuberM.StupperichE.BotheH.BuckelW.KonratR.et al (2000). Native corrinoids from clostridium cochlearium are adeninylcobamides: spectroscopic analysis and identification of pseudovitamin B(12) and factor A.J. Bacteriol.1824773–4782. 10.1128/JB.182.17.4773-4782.2000
47
HoppeH. G. (2003). Phosphatase activity in the sea.Hydrobiologia493187–200. 10.1023/A:1025453918247
48
HuangB. Q.OuL. J.HongH. S.LuoH. W.WangD. Z. (2005). Bioavailability of dissolved organic phosphorus compounds to typical harmful dinoflagellate Prorocentrum donghaiense Lu.Mar. Pollut. Bull.51838–844. 10.1016/j.marpolbul.2005.02.035
49
HuangX. Q.ZhuJ. M.CaiZ. H.LaoY. M.JinH.YuK.et al (2018). Profiles of quorum sensing (QS)-related sequences in phycospheric microorganisms during a marine dinoflagellate bloom, as determined by a metagenomic approach.Microbiol. Res.2171–13. 10.1016/j.micres.2018.08.015
50
HuglerM.SievertS. M. (2011). Beyond the calvin cycle: autotrophic carbon fixation in the ocean.Annu. Rev. Mar. Sci.3261–289. 10.1146/annurev-marine-120709-142712
51
JonesJ. (1972). Studies on freshwater bacteria: association with algae and alkaline phosphatase activity.J. Ecol.6059–75. 10.1007/s11356-021-16255-6
52
KakinumaM.NakamotoC.KishiK.CouryD. A.AmanoH. (2017). Isolation and functional characterization of an ammonium transporter gene, PyAMT1, related to nitrogen assimilation in the marine macroalga Pyropia yezoensis (Rhodophyta).Mar. Environ. Res.12876–87. 10.1016/j.marenvres.2016.08.007
53
KanehisaM.FurumichiM.SatoY.Ishiguro-WatanabeM.TanabeM. (2021). KEGG: integrating viruses and cellular organisms.Nucleic Acids Res.49D545–D551. 10.1093/nar/gkaa970
54
KangD. W. D.LiF.KirtonE.ThomasA.EganR.AnH.et al (2019). MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies.PeerJ7:e7359. 10.7717/peerj.7359
55
KarrayR.KarrayF.LoukilS.MhiriN.SayadiS. (2017). Anaerobic co-digestion of tunisian green macroalgae ulva rigida with sugar industry wastewater for biogas and methane production enhancement.Waste Manag.61171–178. 10.1016/j.wasman.2016.11.042
56
KlawonnI.BonagliaS.WhitehouseM. J.LittmannS.TienkenD.KuypersM. M. M.et al (2019). Untangling hidden nutrient dynamics: rapid ammonium cycling and single-cell ammonium assimilation in marine plankton communities.ISME J.131960–1974. 10.1038/s41396-019-0386-z
57
KuypersM. M. M.MarchantH. K.KartalB. (2018). The microbial nitrogen-cycling network.Nat. Rev. Microbiol.16263–276. 10.1038/nrmicro.2018.9
58
LamattinaL.Garcia-MataC.GrazianoM.PagnussatG. (2003). Nitric oxide: the versatility of an extensive signal molecule.Annu. Rev. Plant Biol.54109–136. 10.1146/annurev.arplant.54.031902.134752
59
LetunicI.BorkP. (2021). Interactive tree of life (iTOL) v5: an online tool for phylogenetic tree display and annotation.Nucleic Acids Res.49W293–W296. 10.1093/nar/gkab301
60
LinS. J.LitakerR. W.SundaW. G. (2016). Phosphorus physiological ecology and molecular mechanisms in marine phytoplankton.J. Phycol.5210–36. 10.1111/jpy.12365
61
LinX.GuoC. T.LiL.LiT. C.LinS. J. (2019). Non-conventional metal Ion cofactor requirement of dinoflagellate alkaline phosphatase and translational regulation by phosphorus limitation.Microorganisms7:232. 10.3390/microorganisms7080232
62
LinZ. F.ZhongS. L.GriersonD. (2009). Recent advances in ethylene research.J. Exp. Bot.603311–3336. 10.1093/jxb/erp204
63
LuX. D.HealK. R.IngallsA. E.DoxeyA. C.NeufeldJ. D. (2020). Metagenomic and chemical characterization of soil cobalamin production.ISME J.1453–66. 10.1038/s41396-019-0502-0
64
LugtenbergB.KamilovaF. (2009). Plant-growth-promoting rhizobacteria.Annu. Rev. Microbiol.63541–556. 10.1146/annurev.micro.62.081307.162918
65
LynchM. D. J.NeufeldJ. D. (2015). Ecology and exploration of the rare biosphere.Nat. Rev. Microbiol.13217–229. 10.1038/nrmicro3400
66
ManoY.NemotoK. (2012). The pathway of auxin biosynthesis in plants.J. Exp. Bot.632853–2872. 10.1093/jxb/ers091
67
MarshallK.JointI.CallowM. E.CallowJ. A. (2006). Effect of marine bacterial isolates on the growth and morphology of axenic plantlets of the green alga Ulva linza.Microb. Ecol.52302–310. 10.1007/s00248-006-9060-x
68
MatsudaR.HandayaniM. L.SasakiH.TakechiK.TakanoH.TakioS. (2018). Production of indoleacetic acid by strains of the epiphytic bacteria Neptunomonas spp. isolated from the red alga pyropia yezoensis and the seagrass zostera marina.Arch. Microbiol.200255–265. 10.1007/s00203-017-1439-1
69
MatsuoY.SuzukiM.KasaiH.ShizuriY.HarayamaS. (2003). Isolation and phylogenetic characterization of bacteria capable of inducing differentiation in the green alga monostroma oxyspermum.Environ. Microbiol.525–35. 10.1046/j.1462-2920.2003.00382.x
70
MenzelP.NgK. L.KroghA. (2016). Fast and sensitive taxonomic classification for metagenomics with Kaiju.Nat. Commun.7:11257. 10.1038/ncomms11257
71
Molina-FaveroC.CreusC. M.SimontacchiM.PuntaruloS.LamattinaL. (2008). Aerobic nitric oxide production by azospirillum brasilense Sp245 and its influence on root architecture in tomato.Mol. Plant-Microbe Interact.211001–1009. 10.1094/Mpmi-21-7-1001
72
NakanishiK.NishijimaM.NishimuraM.KuwanoK.SagaN. (1996). Bacteria that induce morphogenesis in ulva pertusa (chlorophyta) grown unde axenic conditions1.J. Phycol.32479–482. 10.1111/j.0022-3646.1996.00479.x
73
NoinajN.GuillierM.BarnardT. J.BuchananS. K. (2010). TonB-dependent transporters: regulation, structure, and function.Ann. Rev. Microbiol.6443–60. 10.1146/annurev.micro.112408.134247
74
NowinskiB.SmithC. B.ThomasC. M.EssonK.MarinR.PrestonC. M.et al (2019). Microbial metagenomes and metatranscriptomes during a coastal phytoplankton bloom.Sci. Data6:129. 10.1038/s41597-019-0132-4
75
OlmM. R.BrownC. T.BrooksB.BanfieldJ. F. (2017). dRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication.ISME J.112864–2868. 10.1038/ismej.2017.126
76
ParksD. H.ChuvochinaM.WaiteD. W.RinkeC.SkarshewskiA.ChaumeilP. A.et al (2018). A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life.Nat. Biotechnol36996–1004. 10.1038/nbt.4229
77
ProvasoliL.PintnerI. J. (1980). Bacteria induced polymorphism in an axenic laboratory strain of Ulva lactuca (Chlorophyceae) 1.J. Phycol.16196–201. 10.1111/j.0022-3646.1980.00196.x
78
RamboI. M.DombrowskiN.ConstantL.ErdnerD.BakerB. J. (2020). Metabolic relationships of uncultured bacteria associated with the microalgae Gambierdiscus.Environ. Microbiol.221764–1783. 10.1111/1462-2920.14878
79
RavenJ. A.JohnstonA. M.KublerJ. E.KorbR.McInroyS. G.HandleyL. L.et al (2002). Seaweeds in cold seas: evolution and carbon acquisition.Ann. Bot.90525–536. 10.1093/aob/mcf171
80
RavishankaraA. R.DanielJ. S.PortmannR. W. (2009). Nitrous Oxide (N2O): the dominant ozone-depleting substance emitted in the 21st century.Science326123–125. 10.1126/science.1176985
81
RigonatoJ.AlvarengaD. O.AndreoteF. D.DiasA. C. F.MeloI. S.KentA.et al (2012). Cyanobacterial diversity in the phyllosphere of a mangrove forest.Fems Microbiol. Ecol.80312–322. 10.1111/j.1574-6941.2012.01299.x
82
Rodriguez-RL. M.GunturuS.TiedjeJ. M.ColeJ. R.KonstantinidisK. T. (2018). Nonpareil 3: fast estimation of metagenomic coverage and sequence diversity.Msystems3e39–e18. 10.1128/mSystems.00039-18
83
Roth-SchulzeA. J.PintadoJ.Zozaya-ValdesE.CremadesJ.RuizP.KjellebergS.et al (2018). Functional biogeography and host specificity of bacterial communities associated with the marine green alga ulva spp.Mol. Ecol.271952–1965. 10.1111/mec.14529
84
Sanudo-WilhelmyS. A.Gomez-ConsarnauL.SuffridgeC.WebbE. A. (2014). The role of B vitamins in marine biogeochemistry.Ann. Rev. Mar. Sci.6339–367. 10.1146/annurev-marine-120710-100912
85
SeymourJ. R.AminS. A.RainaJ. B.StockerR. (2017). Zooming in on the phycosphere: the ecological interface for phytoplankton-bacteria relationships.Nat. Microbiol.2:17065. 10.1038/nmicrobiol.2017.65
86
SinghP.JajooA.SahayA.BhartiS. (2007). Relation between the mode of binding of nitrite and the energy distribution between the two photosystems.Physiol. Plantarum129447–454. 10.1111/j.1399-3054.2006.00808.x
87
SohaskeyC. D.WayneL. G. (2003). Role of narK2X and narGHJI in hypoxic upregulation of nitrate reduction by mycobacterium tuberculosis.J. Bacteriol.1857247–7256. 10.1128/Jb.185.24.7247-7256.2003
88
SpaepenS.VanderleydenJ.RemansR. (2007). Indole-3-acetic acid in microbial and microorganism-plant signaling.Fems Microbiol. Rev.31425–448. 10.1111/j.1574-6976.2007.00072.x
89
SuzumuraM.HashihamaF.YamadaN.KinouchiS. (2012). Dissolved phosphorus pools and alkaline phosphatase activity in the euphotic zone of the western North Pacific Ocean.Front. Microbiol.3:99. 10.3389/fmicb.2012.00099
90
TakahashiH.KoprivaS.GiordanoM.SaitoK.HellR. (2011). Sulfur assimilation in photosynthetic organisms: molecular functions and regulations of transporters and assimilatory Enzymes.Ann. Rev. Plant Biol.62157–184. 10.1146/annurev-arplant-042110-103921
91
TarafdarJ. C.YadavR. S.MeenaS. C. (2001). Comparative efficiency of acid phosphatase originated from plant and fungal sources.J. Plant Nutr. Soil Sci.164279–282. 10.1002/1522-2624(200106)164:3<279::aid-jpln279>3.0.co;2-l
92
TeelingH.FuchsB. M.BecherD.KlockowC.GardebrechtA.BennkeC. M.et al (2012). Substrate-controlled succession of marine bacterioplankton populations induced by a phytoplankton bloom.Science336608–611. 10.1126/science.1218344
93
ThapaS.PrasannaR. (2018). Prospecting the characteristics and significance of the phyllosphere microbiome.Ann. Microbiol.68229–245. 10.1007/s13213-018-1331-5
94
TriasR.Garcia-LledoA.SanchezN.Lopez-JuradoJ. L.HallinS.BanerasL. (2012). Abundance and composition of epiphytic bacterial and archaeal ammonia oxidizers of marine red and brown macroalgae.Appl. Environ. Microbiol.78318–325. 10.1128/Aem.05904-11
95
TromasA.Perrot-RechenmannC. (2010). Recent progress in auxin biology.Comptes Rendus Biol.333297–306. 10.1016/j.crvi.2010.01.005
96
TwiggM. S.TaitK.WilliamsP.AtkinsonS.CamaraM. (2014). Interference with the germination and growth of Ulva zoospores by quorum-sensing molecules from Ulva-associated epiphytic bacteria.Environ. Microbiol.16445–453. 10.1111/1462-2920.12203
97
UjiT.EndoH.MizutaH. (2020). Sexual reproduction via a 1-aminocyclopropane-1-carboxylic acid-dependent pathway through redox modulation in the marine red alga pyropia yezoensis (rhodophyta).Front. Plant Sci.11:609686. 10.3389/fpls.2020.609686
98
UritskiyG. V.DiRuggieroJ.TaylorJ. (2018). MetaWRAP-a flexible pipeline for genome-resolved metagenomic data analysis.Microbiome6:158. 10.1186/s40168-018-0541-1
99
VadiveluV. M.YuanZ. G.FuxC.KellerJ. (2006). The inhibitory effects of free nitrous acid on the energy generation and growth processes of an enriched nitrobacter culture.Environ. Sci. Technol.404442–4448. 10.1021/es051694k
100
Van de PoelB.SmetD.Van Der StraetenD. (2015). Ethylene and hormonal cross talk in vegetative growth and development.Plant Physiol.16961–72. 10.1104/pp.15.00724
101
Van PuyveldeS.ClootsL.EngelenK.DasF.MarchalK.VanderleydenJ.et al (2011). Transcriptome analysis of the rhizosphere bacterium azospirillum brasilense reveals an extensive auxin response.Microb. Ecol.61723–728. 10.1007/s00248-011-9819-6
102
van TolH. M.AminS. A.ArmbrustE. V. (2017). Ubiquitous marine bacterium inhibits diatom cell division.ISME J.1131–42. 10.1038/ismej.2016.112
103
VannesteS.FrimlJ. (2009). Auxin: a trigger for change in plant development.Cell1361005–1016. 10.1016/j.cell.2009.03.001
104
VaucheretH.KronenbergerJ.LepingleA.VilaineF.BoutinJ. P.CabocheM. (1992). Inhibition of tobacco nitrite reductase activity by expression of antisense RNA.Plant J.2559–569.
105
VollmersJ.FrentrupM.RastP.JoglerC.KasterA. K. (2017). Untangling genomes of novel planctomycetal and verrucomicrobial species from monterey bay kelp forest metagenomes by refined binning.Front. Microbiol.8:472. 10.3389/fmicb.2017.00472
106
WangJ.MaoY.DuG.LiX.TangX. (2021). On microbial community of pyropia haitanensis by metagenomic analysis.J. Oceanol. Limnol.391091–1102. 10.1007/s00343-020-0189-0
107
WangR. C.XingX. J.CrawfordN. (2007). Nitrite acts as a transcriptome signal at micromolar concentrations in arabidopsis roots.Plant Physiol.1451735–1745. 10.1104/pp.107.108944
108
WangW.WuL.XuK.XuY.JiD.ChenC.et al (2020). The cultivation of pyropia haitanensis has important impacts on the seawater microbial community.J. Appl. Phycol.322561–2573. 10.1007/s10811-020-02068-6
109
WichardT. (2015). Exploring bacteria-induced growth and morphogenesis in the green macroalga order ulvales (chlorophyta).Front. Plant Sci.6:86. 10.3389/fpls.2015.00086
110
WuY. W.SimmonsB. A.SingerS. W. (2016). MaxBin 2.0: an automated binning algorithm to recover genomes from multiple metagenomic datasets.Bioinformatics32605–607. 10.1093/bioinformatics/btv638
111
XieB.BishopS.StessmanD.WrightD.SpaldingM. H.HalversonL. J. (2013). Chlamydomonas reinhardtii thermal tolerance enhancement mediated by a mutualistic interaction with vitamin B12-producing bacteria.ISME J.71544–1555. 10.1038/ismej.2013.43
112
XuH.LuoX.QianJ.PangX.SongJ.QianG.et al (2012). FastUniq: a fast de novo duplicates removal tool for paired short reads.PLoS One7:e52249. 10.1371/journal.pone.0052249
113
YanY. W.YangH. C.TangL.LiJ.MaoY. X.MoZ. L. (2019). Compositional shifts of bacterial communities associated with pyropia yezoensis and surrounding seawater co-occurring with red rot disease.Front. Microbiol.10:1666. 10.3389/fmicb.2019.01666
114
YiS.SethE. C.MenY. J.StablerS. P.AllenR. H.Alvarez-CohenL.et al (2012). Versatility in corrinoid salvaging and remodeling pathways supports corrinoid-dependent metabolism in dehalococcoides mccartyi.Appl. Environ. Microbiol.787745–7752. 10.1128/Aem.02150-12
115
ZarzyckiJ.BrechtV.MullerM.FuchsG. (2009). Identifying the missing steps of the autotrophic 3-hydroxypropionate CO2 fixation cycle in chloroflexus aurantiacus.Proc. Natl. Acad. Sci. U.S.A.10621317–21322. 10.1073/pnas.0908356106
116
ZhangJ.LiuS.RenJ.WuY.ZhangG. (2007). Nutrient gradients from the eutrophic changjiang (yangtze River) estuary to the oligotrophic kuroshio waters and re-evaluation of budgets for the east china sea shelf.Progr. Oceanogr.74449–478.
117
ZhangP. F.JinT.SahuS. K.XuJ.ShiQ.LiuH.et al (2019). The distribution of tryptophan-dependent indole-3-acetic acid synthesis pathways in bacteria unraveled by large-scale genomic analysis.Molecules24:1411. 10.3390/molecules24071411
118
ZhouY.GandaL.LimM.YuanZ. G.KjellebergS.NgW. J. (2010). Free nitrous acid (FNA) inhibition on denitrifying poly-phosphate accumulating organisms (DPAOs).Appl. Microbiol. Biotechnol.88359–369. 10.1007/s00253-010-2780-3
119
ZhouY. P.ZhangY. M.LiF. F.TanL. J.WangJ. T. (2017). Nutrients structure changes impact the competition and succession between diatom and dinoflagellate in the east china sea.Sci. Total Environ.574499–508. 10.1016/j.scitotenv.2016.09.092
Summary
Keywords
Pyropia haitanensis, microbial community, metagenome, binning, microbial metabolic function
Citation
Wang J, Tang X, Mo Z and Mao Y (2022) Metagenome-Assembled Genomes From Pyropia haitanensis Microbiome Provide Insights Into the Potential Metabolic Functions to the Seaweed. Front. Microbiol. 13:857901. doi: 10.3389/fmicb.2022.857901
Received
19 January 2022
Accepted
28 February 2022
Published
23 March 2022
Volume
13 - 2022
Edited by
Zhiyong Li, Shanghai Jiao Tong University, China
Reviewed by
Nianjun Xu, Ningbo University, China; Guangce Wang, Institute of Oceanology (CAS), China
Updates
Copyright
© 2022 Wang, Tang, Mo and Mao.
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: Yunxiang Mao, yxmao@hntou.edu.cn
This article was submitted to Aquatic Microbiology, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.