Luhua Yang1,2
Shun Li
Yunzeng Zhang
Cheng Gao
Likai Hao- 1Key Laboratory of Marine Environmental Corrosion and Biofouling, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China
- 2Key Laboratory of Urban Environment and Health, Ningbo Urban Environment Observation and Research Station, Institute of Urban Environment, Chinese Academy of Sciences, Xiamen, China
- 3State Key Laboratory of Environmental Geochemistry, Institute of Geochemistry, Chinese Academy of Sciences, Guiyang, China
- 4Guizhou Academy of Tobacco Science, Guiyang, China
- 5University of Chinese Academy of Sciences, Beijing, China
- 6Zhejiang Key Laboratory of Urban Environmental Processes and Pollution Control, CAS Haixi Industrial Technology Innovation Center in Beilun, Ningbo, China
- 7College of Bioscience and Biotechnology, Yangzhou University, Yangzhou, China
- 8State Key Laboratory of Microbial Resources, Institute of Microbiology, Chinese Academy of Sciences, Beijing, China
- 9Bei Bu Zhan Qu CDC, Shenyang, China
- 10CAS Center for Excellence in Quaternary Science and Global Change, Xi’an, China
The plant microbiota is believed to be an accessory genome that extends plant functions, forming holobionts together with the host plant. Plant disease resistance, therefore, is inextricably linked with plant microbiota, which play important roles in plant growth and health. To explore the relationship between plant microbiota and disease resistance, we investigated the tobacco microbiome of two varieties with contrasting disease-resistance levels to bacterial wilt and black shank diseases. Comparative microbiome analysis indicated that the resistant variety assembled a distinct microbiota with higher network complexity and diversity. While Pseudomonas and Ensifer, which contain biocontrol and beneficial members, were enriched in the rhizosphere of the resistant variety, Ralstonia, a genus including the known causative pathogen, was enriched in the susceptible variety. Metagenome sequencing revealed that biocontrol functions, such as hydrogen cyanide synthase, pyochelin biosynthesis, and arthrofactin-type cyclic lipopeptide synthetase, were more abundant in the resistant variety. Further analysis indicated that contigs encoding the corresponding genes were mostly assigned to Pseudomonas. Among all the metagenome-assembled genomes, positive selection was suggested in the genome assigned to Pseudomonas only in the rhizosphere of the resistant variety. The search of biosynthetic gene clusters in the Pseudomonas genome revealed a non-ribosomal peptide synthetase, the compound of which was brabantamide A, with known antimicrobial activity. Collectively, our study suggests that the plant microbiota might be involved in microbe-mediated disease resistance. Particularly, our results highlight Pseudomonas in the rhizosphere of the disease-resistant variety as a promising biocontrol candidate. Our study may facilitate further screening of bacterial isolates and the targeted design of microbial communities.
Introduction
Recent studies have highlighted that plants cannot be viewed as standalone entities (Vandenkoornhuyse et al., 2015; Wei and Jousset, 2017) but rather as a part of the functional unit interdependent with their associated microbiota, which is commonly referred to as the “plant holobiont” (Sanchez-Canizares et al., 2017; Trivedi et al., 2020; Lyu et al., 2021). The plant microbiota is believed to extend plant functions as an accessory genome and play critical roles in plant growth and health (Chialva et al., 2022). Consequently, plant microbiota are inextricably linked to plant disease resistance (Vannier et al., 2019), with growing evidence indicating that they can confer disease resistance to host plants at both individual and community levels (Carrion et al., 2019; Choi et al., 2020; Gu et al., 2020; Matsumoto et al., 2021; Li et al., 2022).
Recently, studies have pinpointed the specific plant genes within barley and tomato that actively shape the plant microbiome (Escudero-Martinez et al., 2022; Oyserman et al., 2022), which could mediate disease resistance. This discovery suggests a profound interconnection between the inherent plant disease resistance and reciprocal traits exhibited by rhizobacteria. Therefore, gaining insights into the plant microbiome of varieties with varying disease resistance becomes pivotal for unraveling resistance mechanisms, particularly in plants lacking known host-resistance genes. While most previous studies emphasized the state with the presence of pathogens (Gao et al., 2021; Yin et al., 2021; Gu et al., 2022), very few have evaluated the healthy states under natural conditions. Accordingly, the interplay between the natural microbiome and the resulting disease performance remains elusive.
Tobacco (Nicotiana tabacum) is one of the most economically significant agricultural crops worldwide. Bacterial wilt (caused by Ralstonia solanacearum) and black shank (caused by Phytophthora nicotianae) are two major diseases affecting tobacco plants (Gallup et al., 2018; Ahmed et al., 2022). To explore the relationship between disease resistance and plant microbiota, we carefully selected two tobacco varieties with contrasting disease-resistance levels. We hypothesize that the disease-resistant variety hosts a more diverse and stable natural microbiome, enriched with a higher abundance of beneficial microorganisms. In contrast, the disease-susceptible variety exhibits a less diverse microbiome. These findings aim to contribute to a comprehensive understanding of plant health, potentially enhancing agricultural sustainability.
Materials and methods
Sample collection and DNA extraction
Two varieties of flue-cured tobacco (Nicotiana tabacum) were chosen for investigation: HD and GZ36. HD, a conventional variety widely cultivated, exhibits high susceptibility to black shank (caused by Phytophthora nicotianae) and bacterial wilt (caused by Ralstonia solanacearum). In contrast, GZ36, a recently developed variety, showcases remarkable resistance to prevalent tobacco diseases.
The experimental cultivation took place in fields within the tobacco-planting farmlands of the typical karst region in Tianma, Anshun, Guizhou, China (26° 24′ 19.548” N, 106°15′24.588″ E). Destructive sampling was performed on 26 July 2020, precisely during the topping stage, approximately 10 days subsequent to the elimination of terminal buds. The selection of this sampling point was underpinned by two principal considerations. First, the topping stage signifies a pivotal juncture in the growth and production of tobacco plants, marking the transition from the reproductive phase to the vegetative phase. Second, disease resistance typically experiences a decline after the topping stage.
The bulk soil, roots with soil attached, and leaves were sampled and promptly transferred to the laboratory on ice. Triplicates were taken for each variety. The soil cores (at 0–20 cm depth) from each field were collected and pooled as one bulk soil replication, containing approximately 80 g. The plants were further divided into four compartments, as described below. First, the loosely attached soil was removed by shaking. The roots were then placed into a 50-ml Falcon tube with PBS buffer and were washed on a shaking platform for 20 min at 180 rpm. The washing buffer was then subjected to centrifugation (1,500 g, 20 min), and the pellet obtained was defined as the rhizosphere compartment. Roots were transferred to a new 50-ml Falcon tube. After a second washing step (20 min, 180 rpm) with PBS buffer and surface sterilization using 75% ethanol (5 min, 180 rpm), roots were flash-frozen with liquid nitrogen and were defined as the root compartment (root endosphere). The leaves were divided into leaf epiphytes and leaf endophytes in the same way as the roots. All compartments were flash-frozen and were stored at −80 °C. The DNA was extracted using the Mag-MK Soil Genome DNA extraction kit (Sangon, China). In total, 30 samples, including two varieties and five compartments, were obtained.
Amplicon sequencing
The primer pairs, 799F and 1193R (Chelius and Triplett, 2001) were used for the bacterial 16S rRNA gene amplification to avoid the co-amplification of the chloroplast. The PCR reaction systems contained 25 μL of 2x Premix Taq (Takara Biotechnology, Japan), 1 μL of each primer (10 μM), and 3 μL of DNA (20 ng/μl) template in a volume of 50 μL. The thermocycling program is as follows: 5 min at 94°C for initialization; 30 cycles of 30-s denaturation at 94°C, 30 s annealing at 52°C, and 30-s extension at 72°C, followed by a 10-min final elongation at 72°C. The sequencing libraries were generated using the NEBNext® Ultra™ II DNA Library Prep Kit for Illumina® (New England Biolabs, United States) following the manufacturer’s instructions. The quality was assessed on the Qubit@ 2.0 Fluorometer (Thermo Fisher Scientific, United States). In the end, the library was sequenced on an Illumina Nova6000 platform (Illumina, United States).
Metagenomic sequencing
After quantifying genomic concentrations using the Qubit™ dsDNA HS Assay Kit (Thermo Fisher Scientific, United States), libraries with an insert length of approximately 500 bp were prepared with a total of approximately 500 ng of DNA. The DNA of each sample was mechanically sheared to approximately 500 bp fragments using Covaris S220 (Covaris, United States). The sheared DNA was processed using the NEBNext Ultra DNA Library Prep Kit for Illumina (New England Biolabs, United States) for end repair and adaptor ligation. The fragmented DNA was recovered using 1x Hieff NGS™ DNA Selection Beads (Yeasen Biotechnology, China). The purified PCR products were assayed using a Qubit 4 Fluorometer (Thermo Fisher Scientific, United States), were pooled with equal concentrations, and were sequenced on the Illumina HiSeq platform (Illumina, United States).
Bioinformatics
The amplicon sequencing analysis was performed using USEARCH (v11.0.667). After the removal of primers and PhiX contamination, the reads were merged and filtered. The clean reads were denoised using the UNOISE3 algorithm (Edgar, 2016a). The taxonomy assignment of the 16S rRNA gene was performed using the RDP training set (v18) (Cole et al., 2014) with the SINTAX taxonomy prediction algorithm (Edgar, 2016b). The reads assigned to chloroplasts and mitochondria were filtered out.
For metagenome sequencing, the quality control and the removal of plant-derived reads were performed by KneadData (v0.10.0).1 The taxonomy was assigned using MetaPhlAn (v3.0.13) (Beghini et al., 2021). The metagenomic reads were assembled by megahit (v1.2.9) (Li et al., 2015). The reads assembled for eukaryotes, viruses, and archaea were filtered out using MMseqs2 (v13.45) (Mirdita et al., 2021). Gene annotation, clustering, and quantification were performed using Prodigal (v1.14.6) (Hyatt et al., 2010), MMseqs2 (v13.45) (Mirdita et al., 2021), and CoverM (v0.3.2),2 respectively. Functional annotation of the reassembled bins was conducted using eggNOG-Mapper (v2.0.1) (Huerta-Cepas et al., 2017, 2019).
The genome reconstruction, including binning and refinement, was conducted using MetaWRAP (v1.3.2) (Uritskiy et al., 2018). The taxonomy of bins was checked by GTDK-Tk (v2.1.1) (Chaumeil et al., 2022) with the database release207_v2 (Parks et al., 2022). Bins with a completeness of more than 80% and contamination of less than 5% were considered to be of high quality and were kept for further analysis. The prediction of biosynthetic gene clusters in the metagenome-assembled genomes was performed using antiSMASH (v6.1.1) (Blin et al., 2021) with the default settings.
The metagenomic reads of each bin were recruited using anvi’o v7.1 (Eren et al., 2021) and the described workflow (Shaiber et al., 2020). The variants in metagenomic read recruitment results were calculated using the microbial population genetics framework implemented in anvi’o. The polymorphism rates of individual codon sites from allele frequencies were calculated to infer the synonymous (pS(site)) and non-synonymous (pN(site)) changes as defined previously (Shenhav and Zeevi, 2020; Kiefl et al., 2022). The ratio of non-synonymous substitution (dN) rates and synonymous substitution (dS) rates (dN/dS) is widely used for studying the selection, which could quantify substitution rates between diverged species (Kryazhimskiy and Plotkin, 2008). In this study, we calculated the per-site rates of synonymous (pS) and non-synonymous polymorphism (pN), as pS and pN could resolve shorter evolutionary timescales than the typical fixation rate. Moreover, pS and pN can be calculated using metagenomic data without complete haplotypes and could specify rates on a per-sample basis, allowing for inter-sample comparisons (Kiefl et al., 2022). In the calculation of pN/pS, the minimum length of contigs was set at 200 bp. The minimum coverage was tenfold, and the minimum occurrence in samples was specified as two. While pN/pS >1 indicates positive selection owing to increased rates of non-synonymous substitutions, pN/pS <1 indicates purifying selection, and pN/pS =1 implies neutral evolution (Rocha et al., 2006).
Statistics
In the study of alpha diversity, we utilized Hill numbers as the diversity index and used the R package “hilldiv” (Alberdi and Gilbert, 2019) for the calculation. Hill numbers, also called effective numbers, were proposed as a unified diversity concept using the scaling parameter q, also known as the diversity order (Chao et al., 2014). The parameter q determines the weight given to the relative abundance of an operational taxonomic unit (OTU) in a community. For instance, when q is 0, the relative abundance is not considered; when q is 1, the OTUs/amplicon sequence variants (ASVs) are weighted precisely according to their relative abundance; and when q is 2, more weight is given to OTUs/ASVs with high relative abundance. The R package “phyloseq” was used for ordination studies (McMurdie and Holmes, 2013). The community dissimilarity in PCoA was tested using adonis in the R package “vegan” (Oksanen et al., 2013).
The rhizosphere and root endosphere samples were combined as root compartments to construct the molecular ecological networks. Only OTUs with more than 20 reads present in no less than half of the samples were included for correlation calculation. The Spearman correlations were performed using the R package “psych” (Revelle, 2019). The cutoff threshold of the correlation and the p-value were set at 0.7 and 0.001, respectively. The same approach was adopted for the phyllosphere samples. In total, four co-occurrence networks, each representing one variety (GZ36/HD) of one compartment (root/leaf), were constructed. The nodes were divided into four subcategories according to the values of within-module connectivity (Zi) and among-module connectivity (Pi): (i) highly linked connector nodes (Zi ≤ 2.5, Pi >0.62), (ii) module hubs (Zi > 2.5, Pi ≤0.62), (iii) network hubs (Zi > 2.5, Pi >0.62), and (iv) peripheral nodes (Zi ≤ 2.5, Pi ≤0.62) (Olesen et al., 2007). The nodes that belonged to connectors, module hubs, or network hubs were considered as keystone OTUs. We then calculated cohesion (Herren and McMahon, 2017) as a property of network complexity using the ‘taxa shuffle’ null model. Negative cohesion indicates disparate niches and/or negative interactions between taxa, whereas positive cohesion reflects niche overlap and/or positive interactions (Hernandez et al., 2021; Yuan et al., 2021).
The taxa specifically inhabiting different microhabitats were identified using the linear discriminant analysis effect size (LEfSe) method (Segata et al., 2011) using the visualization web server ImageGP (Chen et al., 2022). The cutoff of the Wilcoxon and Kruskal–Wallis (KW) tests were both set at 0.01. The cutoff LAD value was set at 3. Differences in functional potentials in the two varieties were determined using the R package “DESeq2” (Love et al., 2014). The results with an adjusted p-value less than 0.05 and a log2 fold change of more than two were considered significant.
Results
Variety effects were strong in the root microbiota but were less prominent in the phyllosphere
In this study, we investigated the root and leaf microbiota of the disease-resistant tobacco variety GZ36 and the disease-susceptible variety HD. In the root samples, contrasting trends of alpha diversity were observed in the rhizosphere and endosphere (Figure 1A). While the diversity of the rhizosphere was higher in GZ36, the diversity of the endosphere was higher in HD. However, statistical differences were only found in the endosphere at q = 0 (equivalent to richness) and q = 1 (equivalent to Shannon index), despite the strong trends observed. In the phyllosphere, both the epiphytic and endophytic compartments exhibited higher diversity in GZ36, although the results were not statistically significant.
Figure 1. A comparative study of the tobacco microbiome between two varieties with contrasting disease resistance levels. The disease-resistant variety GZ36 was represented by the green color, while the disease-susceptible variety HD was represented by the red color. (A) The alpha diversity and (B) beta diversity of the root and leaf microbiome associated with the disease-resistant variety GZ36 and the disease-susceptible variety HD. The asterisk symbol in panel (A) indicates statistical significance. The different symbols in panel (B) represent plant compartments.
The community structure was further compared using principal coordinate analysis (PCoA) (Figure 1B). Varieties appeared to be the strong drivers of root microbiota, which were corroborated by the permutational multivariate analysis (adonis, p < 0.01). However, no statistical significance was observed in the phyllosphere microbiota between the two varieties.
Molecular ecological networks were constructed to investigate the microbial patterns of the two varieties in each compartment (Figure 2A). The complexity, i.e., the connectivities of the microbial communities, was quantified using the index of cohesion. The absolute value of negative cohesion was found to be higher in the disease-resistant variety GZ36 in both root and leaf compartments. The positive cohesion, however, only differed in the root compartment (higher in GZ36), whereas no statistical differences were observed in the phyllosphere (Figure 2B).
Figure 2. Microbial co-occurrence network of the disease-resistant tobacco variety GZ36 exhibited higher complexity than that of the disease-susceptible variety HD. (A) The co-occurrence networks of root and leaf microbiota associated with the two tobacco varieties. (B) The corresponding complexity of each co-occurrence network is quantified by the cohesion index. The absolute value of the cohesion index indicates the degree of complexity. The asterisk denotes statistical significance, while ns represents no statistical significance.
The keystone OTUs of each network were identified according to their topological features (Supplementary material S1; Supplementary Figure S1). In the root compartment, 30 OTUs were detected as key stone OTUs in GZ36 (Actinobacteria, Bacteroidetes, and Proteobacteria), while 40 OTUs were detected as key OTUs in HD (Acidobacteria, Actinobacteria, Candidatus_Saccharibacteria, Gemmatimonadetes, and Proteobacteria) (Table 1). Although the keystone OTUs appeared more diverse in HD, the phylum Bacteroidetes was found exclusively in GZ36 (Table 1). In the phyllosphere, 29 keystone OTUs were observed in GZ36 (including 3 unassigned, 2 Actinobacteria, 15 Bacteroidetes, 2 Firmicutes, and 7 Proteobacteria), whereas only 8 OTUs (2 unassigned, 2 Bacteroidetes, 2 Firmicutes, and 2 Proteobacteria) were found in HD. A much higher proportion of Bacteroidetes were revealed in the resistant variety GZ36.
Table 1. Topological features of the molecular ecological networks, including the number and the taxonomic distribution of key stone OTUs.
Distinct microbiota and KEGG orthologs were enriched in the two tobacco varieties
Based on the results of 16S rRNA gene sequencing, 22 genera were found to be more abundant in the rhizosphere of GZ36, among which Sphingopyxis, Ensifer, Pseudomonas, Glycomyces, and Chryseobacterium were the dominant groups (Figure 3A). In contrast, nine genera displayed higher abundance in the rhizosphere of HD, including Ralstonia, Humibacter, Povalibacter, Novimethylophilus, Ramlibacter, and Methylotenera. In the root endosphere, equal numbers of genera (11 genera) were enriched in each variety. Notably, Agrobacterium and Sphingomonas appeared prominent in the enriched groups of the two varieties. However, in the phyllosphere, contrasting trends were observed. The genera, namely Stenotrophomonas and Rosenbergiella, in the endosphere were enriched exclusively in the variety HD. In the epiphytic leaf compartment, 8 genera were more abundantly present in GZ36, while 10 genera showed higher abundance in HD.
Figure 3. The taxa enriched in each tobacco variety. (A) The genera with statistically different abundances between the two varieties based on amplicon sequencing. (B) The taxa with statistically different abundances between the two varieties based on metagenome sequencing. The green color denotes the disease-resistant variety GZ36, whereas the red color denotes the disease-susceptible variety HD.
A taxonomy comparison based on metagenomic sequencing was also carried out. However, differences between the two varieties were only observed in the root samples but not in the leaf samples (Figure 3B). Similar to the results based on amplicon sequencing, Pseudomonas, Ensifer, and Sphingopyxis were enriched in the rhizosphere of GZ36, whereas two species of Ralstonia (Ralstonia_pickettii and Ralstonia_sp_NFACC01) were enriched in HD.
We further combined the epiphytic and endophytic compartments to form root samples and leaf samples. While similar results were obtained in the root samples as in the rhizosphere, Stenotrophomonas and Kosakonia were found to be enriched in the phyllosphere of HD (Supplementary material S1; Supplementary Figure S2).
We further investigated the functional differences in root and leaf microbiota between the two varieties. In the rhizosphere, 66 KEGG Orthologs (KOs) showed higher abundance in GZ36, whereas 16 KOs were more abundant in HD (Figure 4). Among the KOs enriched in the rhizosphere of GZ36, two KOs were found to be related to hydrogen cyanide synthase (K10814 hydrogen cyanide synthase HcnA and K10816 hydrogen cyanide synthase HcnC), and one KO was related to siderophore production (K12241 pyochelin biosynthesis protein). Notably, the KO with the highest abundance was arthrofactin-type cyclic lipopeptide synthetase (K15659). We then retrieved the contigs assigned to the corresponding KOs. A total of 5 entries were found for hydrogen cyanide synthase, 3 entries for pyochelin biosynthesis protein, and 47 entries for arthrofactin-type cyclic lipopeptide synthetase. A further NCBI BLAST search found that all five contigs containing hydrogen cyanide synthase were assigned to Pseudomonas. Among the three contigs containing genes encoding the pyochelin biosynthesis protein, one contig showed high similarity with Streptomyces, while the other two had no significant results. As to the contigs encoding arthrofactin-type cyclic lipopeptide synthetase, 38 out of the 47 contigs belonged to Pseudomonas, and the rest were identified as Streptomyces, Burkholderia, and Xanthomonas. In the root endosphere, only one KO (K21003) was shown to be more abundant in GZ36, annotated as the polysaccharide biosynthesis protein PslJ.
Figure 4. KEGG orthologs (KOs) enriched in (A) the disease-resistant variety GZ36 and (B) the disease-susceptible variety HD. The size of the dots was in correspondence with the absolute reads number of the KO in the enriched variety. The color of the KOs corresponds to the negative value of the base 10 logarihm of the adjusted p-value.
In the phyllosphere, 226 KOs were more abundant in the leaf epiphytes of GZ36, whereas only 26 KOs showed higher abundance in HD (Supplementary material S2). Among the KOs enriched in GZ36, two KOs were related to polyketide synthase (K12434 polyketide synthase 7 and K12437 polyketide synthase) (Supplementary material S1; Supplementary Figure S3). However, in the endophytic leaf compartment, no statistical differences were observed in KOs between the two varieties.
Positive selection and biosynthetic gene clusters in the root microbiota were associated with the disease-resistant variety
By metagenomic binning and further refinement, we obtained three high-quality metagenome-assembled genomes (MAGs) and eight high-quality MAGs from tobacco root and leaf microbiota, respectively (Supplementary material S1; Supplementary Table S1). The MAGs from the root microbiota were assigned to Rhodocyclaceae, Pseudomonas, and Streptomyces, whereas MAGs from the leaf microbiota were assigned to Agrobacterium, Frateuria, Kosakonia, Neokomagataea, Pantoea, Pseudomonas, and Stenotrophomonas.
To investigate the selective pressures driving protein evolution within microorganisms inhabiting the two varieties, we analyzed the polymorphism of codons in each MAG. By splitting each single codon variant (SCV) into its synonymous and non-synonymous proportions, we estimated per-site rates of synonymous polymorphism and non-synonymous polymorphism as pS and pN. The ratio between pN and pS (pN/pS) was calculated as a proxy for evolutionary pressure.
In the phyllosphere, we found the polymorphism to be overwhelmingly synonymous. None of the pN/pS exceeded 1 (Supplementary material S1; Supplementary Figure S4). In the root microbiota, the occurrence of pN/pS higher than 1 was observed in two MAGs, namely bin 1 (assigned to Streptomyces) and bin 10 (assigned to Pseudomonas) (Figure 5). A closer investigation indicated that 4 proteins in bin 10 showed a pN/pS value higher than 1 in GZ36 but not in HD, suggesting a possible positive selection of the 4 proteins in GZ36. The four proteins were annotated with functions of iron-containing redox enzyme (PF14518.9), inhibitor of vertebrate lysozyme (Ivy) (PF08816.14), carbohydrate kinase (PF01256.20), and YjeF-related protein N-terminus (PF03853.18) (Supplementary material S1; Supplementary Table S2).
Figure 5. The pN/pS in the metagenome-assembled genomes retrieved from the root microbiota. The green dots represent the Pfams from the disease-resistant variety GZ36, while the red dots represent the Pfams from the disease-susceptible variety HD.
We further searched the biosynthetic gene clusters (BGCs) of the MAG (bin 10) assigned to Pseudomonas. A total of six BGCs, including NRPS (non-ribosomal peptide synthetase), NRPS-like, RiPP-like (RiPP: ribosomally synthesized and post-translationally modified peptide product), redox-cofactor, hserlactone (homoserine lactone cluster), and aryl polyene (aryl polyene cluster), were found to be present in this MAG (Supplementary material S1; Supplementary Table S3). Notably, the compound of the NRPS BGC was predicted to be brabantamide A (with 80% gene similarity), and the compound of the hserlactone BGC was likely to be corpeptin A/corpeptin B (72% gene similarity), according to the results of the KnownClusterBlast analysis searched against the Minimum Information about a Biosynthetic Gene (MIBiG) cluster repository.
Discussion
In this study, our comparative study revealed that the disease-resistant tobacco variety GZ36 and the disease-susceptible variety HD assembled distinct bacterial microbiota. Notably, we found that the resistant variety exhibited higher complexity of the microbial network as well as trends of higher diversity. Recent studies have revealed that microbiome network connectivity and composition are linked to disease resistance in strawberry plants (Hassani et al., 2023). Consistent with previous studies, our results suggest that higher complexity and diversity may contribute to the disease resistance of the tobacco variety GZ36.
Potential biocontrol and growth-promoting bacteria were enriched in the rhizosphere of the resistant variety, whereas Ralstonia was enriched in the susceptible variety
Three genera, including Pseudomonas, Ensifer, and Sphingopyxis, appeared with higher abundance in the variety GZ36 based on both amplicon and metagenomic sequencing. The metagenomic analysis suggested that the enriched Pseudomonas belonged to species SHC52. The strain Pseudomonas SHC52 was reported to be a biocontrol agent with antimicrobial activity, which could inhibit the hyphal growth of Rhizoctonia solani and various other fungal, oomycete, and bacterial pathogens (Schmidt et al., 2014). Ensifer is a genus of nitrogen-fixing bacteria, which could promote the growth and health of host plants (Rogel et al., 2001). Moreover, Ensifer adhaerens, as suggested by the metagenome sequencing, could attach to other bacteria and may cause the lysis of these bacteria, acting as a bacterial predator (Germida and Casida, 1983). In contrast to the known beneficial or biocontrol functions of the above two strains, Sphingopyxis was less studied and is most likely a commensal bacteria (Sharma et al., 2021).
The highly disease-susceptible tobacco variety, however, exhibited a higher abundance of Ralstonia, which could be further classified as Ralstonia pickettii and Ralstonia sp. NFACC01. Although the two strains were non-pathogenic (Ryan et al., 2006), they might largely overlap with the pathogen Ralstonia solanacearum in niches and resources (Wei et al., 2018). Thus, it is highly possible that Ralstonia solanacearum is more likely to survive in the rhizosphere of this variety, which further leads to disease susceptibility.
Biocontrol-related KOs were enriched in the rhizosphere of the disease-resistant variety
In the rhizosphere, genes encoding hydrogen cyanide synthase were found to be more abundant in the variety GZ36. Hydrogen cyanide (HCN) is a well-known poison. Some HCN-producing bacteria have important effects against plant diseases (Sehrawat et al., 2022). Historically, this effect has been attributed to a direct poisonous effect (Michelsen and Stougaard, 2012). However, recent studies demonstrated that HCN plays a positive role in the plant response to pathogens independently of its toxicity, and the bacteria stimulate plant growth depending on their capacity to emit HCN (Díaz-Rueda et al., 2023).
In addition, the pyochelin biosynthesis protein showed higher abundance in the disease-resistant variety GZ36. Pyochelin is a siderophore to assimilate iron. It has been shown that siderophores serve a role beyond merely transporting iron. They are crucial mediators, facilitating interactions among microbial communities and the eukaryotic hosts they reside in (Kramer et al., 2019). Competition for iron through the secretion of siderophores has long been considered a major mechanism in the biological control of plant diseases (Miethke and Marahiel, 2007). Particularly, it has been experimentally confirmed that pyochelin-producing bacteria as well as purified pyoverdines could inhibit the growth of Ralstonia solanacearum (Gu et al., 2020).
Notably, arthrofactin-type cyclic lipopeptide synthetase appeared to be the most abundant ortholog in the enriched KOs. Cyclic lipodepsipeptides (CLPs) are secondary metabolites produced by non-ribosomal peptide synthetases (NRPSs) with antimicrobial functions (Raaijmakers et al., 2006). Among them, arthrofactin is one of the most potent biosurfactants (Geudens and Martins, 2018). Several CLP-producing bacteria have been inscribed for the biocontrol of plant diseases associated with CLP metabolite production (Arbeit et al., 2004).
Thus the above-mentioned KOs might be promising candidate traits potentially related to pathogen suppression. Interestingly, contigs containing corresponding genes were mostly assigned to Pseudomonas (and a small part to Streptomyces), which was consistent with the taxonomic results that Pseudomonas was enriched in the resistant variety.
It is suggested that the plant microbiome can determine the outcome of pathogen infections. Plants may establish interactions with protective microbes before pathogens arrive, preventing the onset of diseases (Pereira et al., 2023). While the microbiome is suggested to be linked with the disease resistance of strawberries (Hassani et al., 2023), our results provide possible explanations in terms of the functions provided by the plant microbiota. Recent studies have shown that plant genetic traits could determine the assembly of plant microbiota, which in turn prevent pathogen invasion (Escudero-Martinez et al., 2022; Oyserman et al., 2022). The plant and the associated microbiota thus form a self-reinforcing immunity and a recruitment loop, coining the concept “holobiont.” However, we acknowledge that the results presented in our study represent a plausible correlation rather than concrete causality. Further studies are still needed to corroborate the benefits of the microbiota.
The positive selection indicated by pN/pS
In this study, four proteins were observed with pN/pS ratios >1 in GZ36 but not in HD in the MAG assigned to Pseudomonas, indicating that the four proteins might have been under positive selection in the resistant variety. The protein with the highest pN/pS was assigned to the function of an iron-containing redox enzyme. Iron is not only an essential nutrient for the growth of microorganisms but is also critically involved in the redox response (Cornelis et al., 2011). Competition for iron was shown to drive phytopathogen control in the natural rhizosphere (Gu et al., 2020). The second protein, an inhibitor of vertebrate lysozyme (Ivy), is reported as a stress response to cell wall damage (Abergel et al., 2007). It is possible that the two proteins of Pseudomonas were under positive selection.
NRPS identified in the MAG
The search of BGCs in the MAG assigned to Pseudomonas suggested the presence of NRPS, which was likely to be brabantamide A. Brabantamide A has been isolated from Pseudomonas fluorescens DSM 11579 (Thirkettle et al., 2000) and Pseudomonas brassicacearum (Andersson et al., 2012). The genetic origin and biosynthetic pathways of brabantamide A were elucidated clearly in Pseudomonas sp., SHC52 (Schmidt et al., 2014). Significant antibacterial and antifungal activity was reported with brabantamide A. This is also consistent with the taxonomic results showing that Pseudomonas sp., SHC52 was enriched in the resistant variety. With the high similarity of brabantamide A in our MAG, it is highly possible that the Pseudomonas obtained by metagenome binning is a promising biocontrol agent.
Additional investigations are warranted to isolate the bacterial strain Pseudomonas and substantiate its protective efficacy against plant pathogens. Subsequent research endeavors should focus on elucidating the mechanisms underlying this protective function, potentially through genomic and metabolomic analyses. Moreover, experimentation aimed at determining the transferability and reproducibility of the protective effect through microbial transplantation would provide valuable insights. These efforts could pave the way for innovative applications in agriculture, such as developing biocontrol agents or engineering plant microbiomes for enhanced disease resistance.
Conclusion
In this study, we carried out a comparative microbiome study of two tobacco varieties, which are highly resistant and susceptible to the common tobacco diseases, respectively. Our results revealed that the resistant variety assembled microbiota with higher network complexity and diversity. Notably, Pseudomonas SHC52, a known biocontrol agent, was enriched in the rhizosphere of the resistant variety, whereas Ralstonia, a genus containing the known causative agent, was enriched in the susceptible variety. Moreover, biocontrol functions such as hydrogen cyanide synthase, pyochelin biosynthesis, and arthrofactin-type cyclic lipopeptide synthetase were enriched in the resistant variety, which were largely derived from Pseudomonas, as indicated by sequence retrieval and BLAST analysis. Metagenome-assembled genomes were further obtained by binning. In the assembled genome assigned to Pseudomonas, positive selection was suggested for the resistant variety. In addition, a non-ribosomal peptide synthetase encoding brabantamide A was identified, which was known to have antimicrobial activity. Collectively, our study suggests Pseudomonas as a promising biocontrol agent that might contribute to the disease resistance of the host plant.
Data availability statement
The sequencing datasets presented in this study can be found in GSA database (https://bigd.big.ac.cn/gsa/) under the accession number CRA008960.
Author contributions
LY: Data curation, Formal analysis, Visualization, Writing – original draft. YG: Data curation, Investigation, Methodology, Writing – review & editing. HY: Conceptualization, Data curation, Writing – review & editing. SL: Writing – review & editing. YZ: Writing – review & editing. CG: Writing – review & editing. TW: Writing – review & editing. LH: Conceptualization, Funding acquisition, Project administration, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by grants from the Guizhou Provincial Department of Science and Technology (E2DF028), the Startup Funding of the Chinese Academy of Sciences (2017-020), the National Natural Science Foundation of China (41877400), the National Key Research and Development Project of China (2018YFC1802601), the Strategic Priority Research Program of Chinese Academy of Sciences (XDB40020300), and the State Key Laboratory of Environmental Geochemistry (SKLEG2018911).
Acknowledgments
The authors thank colleagues Dr. Y. Wang, Dr. Q. Huang, Dr. H. Yang, L. Liu, and Y. Fu from the Institute of Geochemistry of the Chinese Academy of Sciences in Guiyang for sample preparation and handling.
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.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
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.2024.1361883/full#supplementary-material
Footnotes
References
Abergel, C., Monchois, V., Byrne, D., Chenivesse, S., Lembo, F., Lazzaroni, J. C., et al. (2007). Structure and evolution of the ivy protein family, unexpected lysozyme inhibitors in gram-negative bacteria. Proc. Natl. Acad. Sci. USA 104, 6394–6399. doi: 10.1073/pnas.0611019104
Alberdi, A., and Gilbert, M. T. P. (2019). A guide to the application of hill numbers to DNA-based diversity analyses. Mol. Ecol. Resour. 19, 804–817. doi: 10.1111/1755-0998.13014
Andersson, P. F., Levenfors, J., and Broberg, A. (2012). Metabolites from Pseudomonas brassicacearum with activity against the pink snow mould causing pathogen Microdochium nivale. BioControl 57, 463–469. doi: 10.1007/s10526-011-9411-5
Ahmed, W., Yang, J., Tan, Y., Munir, S., Liu, Q., Zhang, J., et al. (2022). Ralstonia solanacearum, a deadly pathogen: Revisiting the bacterial wilt biocontrol practices in tobacco and other Solanaceae. Rhizosphere 21: 100479. doi: 10.1016/j.rhisph.2022.100479
Arbeit, R. D., Maki, D., Tally, F. P., Campanaro, E., and Eisenstein, B. I.Daptomycin and Investigators (2004). The safety and efficacy of daptomycin for the treatment of complicated skin and skin-structure infections. Clin. Infect. Dis. 38, 1673–1681. doi: 10.1086/420818
Beghini, F., McIver, L. J., Blanco-Miguez, A., Dubois, L., Asnicar, F., Maharjan, S., et al. (2021). Integrating taxonomic, functional, and strain-level profiling of diverse microbial communities with bioBakery 3. elife 10:65088. doi: 10.7554/eLife.65088
Blin, K., Shaw, S., Kautsar, S. A., Medema, M. H., and Weber, T. (2021). The antiSMASH database version 3: increased taxonomic coverage and new query features for modular enzymes. Nucleic Acids Res. 49, D639–D643. doi: 10.1093/nar/gkaa978
Carrion, V. J., Perez-Jaramillo, J., Cordovez, V., Tracanna, V., de Hollander, M., Ruiz-Buck, D., et al. (2019). Pathogen-induced activation of disease-suppressive functions in the endophytic root microbiome. Science 366, 606–612. doi: 10.1126/science.aaw9285
Chao, A., Chiu, C.-H., and Jost, L. (2014). Unifying species diversity, phylogenetic diversity, functional diversity, and related similarity and differentiation measures through hill numbers. Annu. Rev. Ecol. Evol. Syst. 45, 297–324. doi: 10.1146/annurev-ecolsys-120213-091540
Chaumeil, P. A., Mussig, A. J., Hugenholtz, P., and Parks, D. H. (2022). GTDB-Tk v2: memory friendly classification with the genome taxonomy database. Bioinformatics 38, 5315–5316. doi: 10.1093/bioinformatics/btac672
Chelius, M. K., and Triplett, E. W. (2001). The diversity of Archaea and Bacteria in association with the roots of Zea mays L. Microb. Ecol. 41, 252–263. doi: 10.1007/s002480000087
Chen, T., Liu, Y. X., and Huang, L. (2022). ImageGP: an easy-to-use data visualization web serverfor scientific researchers. iMeta 1:e5. doi: 10.1002/imt2.5
Chialva, M., Lanfranco, L., and Bonfante, P. (2022). The plant microbiota: composition, functions, and engineering. Curr. Opin. Biotechnol. 73, 135–142. doi: 10.1016/j.copbio.2021.07.003
Choi, K., Choi, J., Lee, P. A., Roy, N., Khan, R., Lee, H. J., et al. (2020). Alteration of bacterial wilt resistance in tomato plant by microbiota transplant. Front. Plant Sci. 11:1186. doi: 10.3389/fpls.2020.01186
Cole, J. R., Wang, Q., Fish, J. A., Chai, B., McGarrell, D. M., Sun, Y., et al. (2014). Ribosomal database project: data and tools for high throughput rRNA analysis. Nucleic Acids Res. 42, D633–D642. doi: 10.1093/nar/gkt1244
Cornelis, P., Wei, Q., Andrews, S. C., and Vinckx, T. (2011). Iron homeostasis and management of oxidative stress response in bacteria. Metallomics 3, 540–549. doi: 10.1039/c1mt00022e
Díaz-Rueda, P., de los Ríos, L. M., Romero, L. C., and García, I. (2023). Old poisons, new signaling molecules: the case of hydrogen cyanide. J. Exp. Bot. 74, 6040–6051. doi: 10.1093/jxb/erad317
Edgar, R. C. (2016a). UNOISE2: improved error-correction for Illumina 16S and ITS amplicon sequencing. bioRxiv :081257. doi: 10.1101/081257
Edgar, R. C. (2016b). SINTAX: a simple non-Bayesian taxonomy classifier for 16S and ITS sequences. bioRxiv :074161. doi: 10.1101/074161
Eren, A. M., Kiefl, E., Shaiber, A., Veseli, I., Miller, S. E., Schechter, M. S., et al. (2021). Community-led, integrated, reproducible multi-omics with anvi'o. Nat. Microbiol. 6, 3–6. doi: 10.1038/s41564-020-00834-3
Escudero-Martinez, C., Coulter, M., Alegria Terrazas, R., Foito, A., Kapadia, R., Pietrangelo, L., et al. (2022). Identifying plant genes shaping microbiota composition in the barley rhizosphere. Nat. Commun. 13:3443. doi: 10.1038/s41467-022-31022-y
Gallup, C. A., McCorkle, K. L., Ivors, K. L., and Shew, D. (2018). Characterization of the Black Shank Pathogen, Phytophthora nicotianae, Across North Carolina Tobacco Production Areas. Plant dis. 102, 1108–1114. doi: 10.1094/PDIS-02-17-0295-RE
Gao, M., Xiong, C., Gao, C., Tsui, C. K. M., Wang, M. M., Zhou, X., et al. (2021). Disease-induced changes in plant microbiome assembly and functional adaptation. Microbiome 9:187. doi: 10.1186/s40168-021-01138-2
Germida, J. J., and Casida, L. E. (1983). Ensifer adhaerens predatory activity against other Bacteria in soil, as monitored by indirect phage analysis. Appl. Environ. Microbiol. 45, 1380–1388. doi: 10.1128/aem.45.4.1380-1388.1983
Geudens, N., and Martins, J. C. (2018). Cyclic Lipodepsipeptides from Pseudomonas spp.—biological Swiss-Army knives. Front. Microbiol. 9:1867. doi: 10.3389/fmicb.2018.01867
Gu, Y., Banerjee, S., Dini-Andreote, F., Xu, Y., Shen, Q., Jousset, A., et al. (2022). Small changes in rhizosphere microbiome composition predict disease outcomes earlier than pathogen density variations. ISME J. 16, 2448–2456. doi: 10.1038/s41396-022-01290-z
Gu, S., Wei, Z., Shao, Z., Friman, V. P., Cao, K., Yang, T., et al. (2020). Competition for iron drives phytopathogen control by natural rhizosphere microbiomes. Nat. Microbiol. 5, 1002–1010. doi: 10.1038/s41564-020-0719-8
Hassani, M. A., Gonzalez, O., Hunter, S. S., Holmes, G. J., Hewavitharana, S. S., Ivors, K., et al. (2023). Microbiome network connectivity and composition linked to disease resistance in strawberry plants. Phytobio. J. 7, 298–311. doi: 10.1094/PBIOMES-10-22-0069-R
Hernandez, D. J., David, A. S., Menges, E. S., Searcy, C. A., and Afkhami, M. E. (2021). Environmental stress destabilizes microbial networks. ISME J. 15, 1722–1734. doi: 10.1038/s41396-020-00882-x
Herren, C. M., and McMahon, K. D. (2017). Cohesion: a method for quantifying the connectivity of microbial communities. ISME J. 11, 2426–2438. doi: 10.1038/ismej.2017.91
Huerta-Cepas, J., Forslund, K., Coelho, L. P., Szklarczyk, D., Jensen, L. J., von Mering, C., et al. (2017). Fast genome-wide functional annotation through Orthology assignment by eggNOG-mapper. Mol. Biol. Evol. 34, 2115–2122. doi: 10.1093/molbev/msx148
Huerta-Cepas, J., Szklarczyk, D., Heller, D., Hernandez-Plaza, A., Forslund, S. K., Cook, H., et al. (2019). eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 47, D309–D314. doi: 10.1093/nar/gky1085
Hyatt, D., Chen, G. L., Locascio, P. F., Land, M. L., Larimer, F. W., and Hauser, L. J. (2010). Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinform. 11:119. doi: 10.1186/1471-2105-11-119
Kiefl, E., Esen, O. C., Miller, S. E., Kroll, K. L., Willis, A. D., Rappé, M. S., et al. (2022). Structure-informed microbial population genetics elucidate selective pressures that shape protein evolution. Sci. Adv. 9:eabq4632. doi: 10.1126/sciadv.abq4632
Kramer, J., Özkaya, Ö., and Kümmerli, R. (2019). Bacterial siderophores in community and host interactions. Nat. Rev. Microbiol. 18, 152–163. doi: 10.1038/s41579-019-0284-4
Kryazhimskiy, S., and Plotkin, J. B. (2008). The population genetics of dN/dS. PLoS Genet. 4:e1000304. doi: 10.1371/journal.pgen.1000304
Li, D., Liu, C. M., Luo, R., Sadakane, K., and Lam, T. W. (2015). MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 31, 1674–1676. doi: 10.1093/bioinformatics/btv033
Li, M., Pommier, T., Yin, Y., Wang, J., Gu, S., Jousset, A., et al. (2022). Indirect reduction of Ralstonia solanacearum via pathogen helper inhibition. ISME J. 16, 868–875. doi: 10.1038/s41396-021-01126-2
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8
Lyu, D., Zajonc, J., Page, A., Tanney, C. A. S., Shah, A., Monjezi, N., et al. (2021). Plant Holobiont theory: the Phytomicrobiome plays a central role in evolution and success. Microorganisms 9:675. doi: 10.3390/microorganisms9040675
Matsumoto, H., Fan, X. Y., Wang, Y., Kusstatscher, P., Duan, J., Wu, S. L., et al. (2021). Bacterial seed endophyte shapes disease resistance in rice. Nat. Plants 7, 60–72. doi: 10.1038/s41477-020-00826-5
McMurdie, P. J., and Holmes, S. (2013). Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8:e61217. doi: 10.1371/journal.pone.0061217
Michelsen, C. F., and Stougaard, P. (2012). Hydrogen cyanide synthesis and antifungal activity of the biocontrol strain Pseudomonas fluorescens In5 from Greenland is highly dependent on growth medium. Can. J. Microbiol. 58, 381–390. doi: 10.1139/w2012-004
Miethke, M., and Marahiel, M. A. (2007). Siderophore-based iron acquisition and pathogen control. Microbiol. Mol. Biol. Rev. 71, 413–451. doi: 10.1128/MMBR.00012-07
Mirdita, M., Steinegger, M., Breitwieser, F., Soding, J., and Levy Karin, E. (2021). Fast and sensitive taxonomic assignment to metagenomic contigs. Bioinformatics 37, 3029–3031. doi: 10.1093/bioinformatics/btab184
Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O’hara, R. B., et al. (2013). vegan: Community Ecology Package. Software. Available at: http://CRAN.R-project.org/package=vegan
Olesen, J. M., Bascompte, J., Dupont, Y. L., and Jordano, P. (2007). The modularity of pollination networks. Proc. Natl. Acad. Sci. USA 104, 19891–19896. doi: 10.1073/pnas.0706375104
Oyserman, B. O., Flores, S. S., Griffioen, T., Pan, X., van der Wijk, E., Pronk, L., et al. (2022). Disentangling the genetic basis of rhizosphere microbiome assembly in tomato. Nat. Commun. 13:3228. doi: 10.1038/s41467-022-30849-9
Parks, D. H., Chuvochina, M., Rinke, C., Mussig, A. J., Chaumeil, P. A., and Hugenholtz, P. (2022). GTDB: an ongoing census of bacterial and archaeal diversity through a phylogenetically consistent, rank normalized and complete genome-based taxonomy. Nucleic Acids Res. 50, D785–D794. doi: 10.1093/nar/gkab776
Pereira, L. B., Thomazella, D. P. T., and Teixeira, P. J. P. L. (2023). Plant-microbiome crosstalk and disease development. Curr. Opin. Plant Biol. 72:102351. doi: 10.1016/j.pbi.2023.102351
Raaijmakers, J. M., de Bruijn, I., and de Kock, M. J. (2006). Cyclic lipopeptide production by plant-associated Pseudomonas spp.: diversity, activity, biosynthesis, and regulation. Mol. Plant-Microbe Interact. 19, 699–710. doi: 10.1094/MPMI-19-0699
Revelle, W. (2019). Psych: procedures for psychological. R Package Version 1.9.12. Northwestern University, Evanston: Psychometric, and Personality Research. Available at: https://CRAN.R-project.org/package=psych.
Rocha, E. P., Smith, J. M., Hurst, L. D., Holden, M. T., Cooper, J. E., Smith, N. H., et al. (2006). Comparisons of dN/dS are time dependent for closely related bacterial genomes. J. Theor. Biol. 239, 226–235. doi: 10.1016/j.jtbi.2005.08.037
Rogel, M. A., Hernandez-Lucas, I., Kuykendall, L. D., Balkwill, D. L., and Martinez-Romero, E. (2001). Nitrogen-fixing nodules with Ensifer adhaerens harboring Rhizobium tropici symbiotic plasmids. Appl. Environ. Microbiol. 67, 3264–3268. doi: 10.1128/AEM.67.7.3264-3268.2001
Ryan, M. P., Pembroke, J. T., and Adley, C. C. (2006). Ralstonia pickettii: a persistent gram-negative nosocomial infectious organism. J. Hosp. Infect. 62, 278–284. doi: 10.1016/j.jhin.2005.08.015
Sanchez-Canizares, C., Jorrin, B., Poole, P. S., and Tkacz, A. (2017). Understanding the holobiont: the interdependence of plants and their microbiome. Curr. Opin. Microbiol. 38, 188–196. doi: 10.1016/j.mib.2017.07.001
Schmidt, Y., van der Voort, M., Crusemann, M., Piel, J., Josten, M., Sahl, H. G., et al. (2014). Biosynthetic origin of the antibiotic cyclocarbamate brabantamide A (SB-253514) in plant-associated Pseudomonas. Chembiochem 15, 259–266. doi: 10.1002/cbic.201300527
Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:R60. doi: 10.1186/gb-2011-12-6-r60
Sehrawat, A., Sindhu, S., and Glick, B. R. (2022). Hydrogen cyanide production by soil bacteria: biological control of pests and promotion of plant growth in sustainable agriculture. Pedosphere 32, 15–38. doi: 10.1016/S1002-0160(21)60058-9
Shaiber, A., Willis, A. D., Delmont, T. O., Roux, S., Chen, L. X., Schmid, A. C., et al. (2020). Functional and genetic markers of niche partitioning among enigmatic members of the human oral microbiome. Genome Biol. 21:292. doi: 10.1186/s13059-020-02195-w
Sharma, M., Khurana, H., Singh, D. N., and Negi, R. K. (2021). The genus Sphingopyxis: systematics, ecology, and bioremediation potential—A review. J. Environ. Manag. 280:111744. doi: 10.1016/j.jenvman.2020.111744
Shenhav, L., and Zeevi, D. (2020). Resource conservation manifests in the genetic code. Science 370, 683–687. doi: 10.1126/science.aaz9642
Thirkettle, J., Alvarez, E., Boyd, H., Brown, M., Diez, E., Hueso, J., et al. (2000). SB-253514 and analogues; novel inhibitors of lipoprotein-associated phospholipase A2 produced by Pseudomonas fluorescens DSM 11579. I. Fermentation of producing strain, isolation and biological activity. J. Antibiot. (Tokyo) 53, 664–669. doi: 10.7164/antibiotics.53.664
Trivedi, P., Leach, J. E., Tringe, S. G., Sa, T., and Singh, B. K. (2020). Plant-microbiome interactions: from community assembly to plant health. Nat. Rev. Microbiol. 18, 607–621. doi: 10.1038/s41579-020-0412-1
Uritskiy, G. V., DiRuggiero, J., and Taylor, J. (2018). MetaWRAP-a flexible pipeline for genome-resolved metagenomic data analysis. Microbiome 6:158. doi: 10.1186/s40168-018-0541-1
Vandenkoornhuyse, P., Quaiser, A., Duhamel, M., Le Van, A., and Dufresne, A. (2015). The importance of the microbiome of the plant holobiont. New Phytol. 206, 1196–1206. doi: 10.1111/nph.13312
Vannier, N., Agler, M., and Hacquard, S. (2019). Microbiota-mediated disease resistance in plants. PLoS Pathog. 15:e1007740. doi: 10.1371/journal.ppat.1007740
Wei, Z., Hu, J., Gu, Y. A., Yin, S., Xu, Y., Jousset, A., et al. (2018). Ralstonia solanacearum pathogen disrupts bacterial rhizosphere microbiome during an invasion. Soil Biol. Biochem. 118, 8–17. doi: 10.1016/j.soilbio.2017.11.012
Wei, Z., and Jousset, A. (2017). Plant breeding Goes microbial. Trends Plant Sci. 22, 555–558. doi: 10.1016/j.tplants.2017.05.009
Yin, C., Casa Vargas, J. M., Schlatter, D. C., Hagerty, C. H., Hulbert, S. H., and Paulitz, T. C. (2021). Rhizosphere community selection reveals bacteria associated with reduced root disease. Microbiome 9:86. doi: 10.1186/s40168-020-00997-5
Keywords: plant microbiota, tobacco microbiome, disease resistance, positive selection, biocontrol
Citation: Yang L, Guo Y, Yang H, Li S, Zhang Y, Gao C, Wei T and Hao L (2024) Distinct microbiota assembly and functional patterns in disease-resistant and susceptible varieties of tobacco. Front. Microbiol. 15:1361883. doi: 10.3389/fmicb.2024.1361883
Edited by:
Marika Pellegrini, University of L’Aquila, ItalyReviewed by:
Sara Fareed Mohamed Wahdan, Suez Canal University, EgyptYong-Xin Liu, Chinese Academy of Agricultural Sciences, China
Copyright © 2024 Yang, Guo, Yang, Li, Zhang, Gao, Wei and Hao. 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: Likai Hao, haolikai@mail.gyig.ac.cn