Tomato microbiome under long‐term organic and conventional farming

Abstract

The compartment niche is the main reason behind the shifts in endophytic bacterial communities. Long-term organic greenhouse exerted limited influence on the variations of endophytic bacterial communities. Organic greenhouse and root had more complex co-occurrence networks than conventional greenhouse and stem, respectively. Cultivable method results found that Protecbacteria, Bacteriodes, and Actinobacteria are the dominant phyla in the endophytes. Greenhouse vegetable production in China has exploded in 2014, with a yield of 260 million tons, representing 35% of the overall vegetable production accounts, and tomato (Solanum lycopersicum L.) as an important economic vegetable reached 55 million tons in 2019 [1]. Organic farming production excludes inorganic fertilizers and operates on the basis of green resource cycling in an ecological manner by applying bio-fertilizer instead of pesticide or fertilizer, which delivered greater ecological services, reduced adverse impacts on the environment, and increased sustainability of agricultural systems compared to conventional farming systems [2, 3]. Intriguingly, a lot of evidence has shown that a positive link between the microbiome and plant health by plant growth promoting rhizosphere (PGPR), which could enhance the nutrition metabolism and resistance to biological stress as plant's second genome [4]. Various compartment niches (such as roots, stems, leaves, flowers, and fruits) are inhabited by a large number of highly diverse microorganisms (i.e., the plant microbiome), including bacteria, fungi, archaea, protists, and so on, together with the host plant form a “symbiotic totality,” interacting and coevolving together [5, 6]. The plant microbiome plays an important role in plant growth and development, nutrient absorption, biotic and abiotic stress resistance, and so on [7]. Endophytes normally completed their whole life cycle within host plants without occurring plant diseases, and their multiproduction might be influenced by host plants and the growing environment [8]. In addition, plant bacteria can provide protection against fungal pathogen diseases by directly promoting plant growth or antibiotic production, nutrient solubilization and transportation, and nitrogen fixation [9, 10]. Existing studies have shown that soil and plant microbiome are determined by many factors, including plant genetic diversity, growth period, soil physical and chemical properties, and soil nutrient status [11, 12]. The majority of plant-associated microbiomes acted as biofertilizers [13], illustrating that the bacterial communities associated with crops are meaningful to unearth potentially beneficial candidates for biologically controlling crop diseases [14]. In recent years, cultivation-independent approaches, especially metagenomics approaches, have allowed a comprehensive analysis of bacterial diversity from various types of plants, including agronomic crops, dicotyledon, and herb plants [7, 9, 14, 15]. Nonetheless, previous studies on microorganisms mostly focused on rhizosphere soil bacterial quantity and community diversity, but the root and stem bacteria under long-term greenhouse trials are largely unknown. Exploring the endophytic bacterial communities will benefit exploring the mechanisms of selectivity adaption in different niches and identifying potential biologic control or promoting growth microbiome candidates in agricultural production [16, 17]. The principal objective of this study was (1) to assess how host (compartment niche) and environmental factors (organic farming or conventional farming) interactively shape tomato endophytic microbiome assemblies and concurrence patterns based on long-term field trials since 2002, and (2) to isolate the potential beneficial tomato microbiomes as a candidate for agricultural inoculants. This study will help to provide comprehensive insight into the tomato endophytic bacterial communities cultivated in a greenhouse agroecosystem and provide useful information for using the beneficial microbiome in agricultural production. To explore the variation of microbiota in long-term organic (ORG) and conventional (CON) systems, the experiments were procedure in two semi-round arched greenhouses (Figure 1A). The change of physicochemical parameters during the key crop season with different agricultural practices is shown in Figure 1B–D. The total nitrogen (TN) concentration increased both in the organic and conventional greenhouses, but the TN concentration was higher in the organic greenhouse than that of the conventional greenhouse (Student's t test, p < 0.001) NH4+–N showed a similar trend with TN, and the concentration of ammonia nitrogen increased with the plant development (Student's t test, p < 0.001), and the ammonium concentration in the organic greenhouse was higher than that of conventional greenhouse (Student's t test, p < 0.001). The soil organic matter content decreased from 54.31 to 33.85 g/kg in the organic greenhouse and from 50.44 to 29.65 g/kg in the conventional greenhouse (Student's t test, p < 0.001). See Supporting Information: Table S1 for the results of the two-way analysis of variance (ANOVA) for soil properties as affected by agricultural manipulation (treat), growing phase (time), and the interaction (treat × time). A total of 648 operational taxonomy units (OTUs) and 2,014,992 high-quality, nonplastid sequences were acquired (Supporting Information: Figure S1). The taxonomic analysis found that the tomato bacteria mainly comprised four phyla (Figure 1F), among which Proteobacteria was the most abundant (40.63%), followed by Firmicutes (21.76%), Actinobacteria (20.3%), and Bacteroidetes (8.7%). In addition, the relative abundance of the top 30 most abundant genera is shown in Supporting Information: Figure S2. The dominant genera belonged to Weissella, Bacillus, Mesorhizobium, and Chryseobacterium.sp in the endophytes. The abundance of Firmicutes in the bacterial community was decreased and the abundance of Proteobacteria increased significantly, especially after the fruiting phase in the root endophytes (Supporting Information: Figure S3), and a similar phenomenon occurred in stem endophytic bacterial communities: the relative abundance of Arthrobacter.sp, Rhizobium.sp gradually increased, while the Bacillus.sp decreased dramatically during the late crop season (Figure 1F and Supporting Information: Figure S2). Results suggested that compartment niche had a strong effect on bacterial diversity (Shannon index), and roots have higher Shannon diversity than stems. Shannon index was higher in the late growth phase than in the early growth phase; however, it had no difference between different manipulation (Figure 2A–C). See Supporting Information: Figure S4 for Chao1 index. Results based on Bray–Curtis distance and ADONIS analysis of the complete data set suggested that the variation of the bacterial community derived from compartment-niche in the organic greenhouse (R2 = 0.076, p = 0.001) (Figure 2D). Moreover, bacterial communities in the conventional greenhouses were similar, of which the variations in microbial communities were mainly explained by compartment niche (R2 = 0.056, p = 0.021) (Figure 2E,F). Nonmetric multidimensional scaling (NMDS) showed that bacteria in four crop development phases separated clusters in the first two coordinate axes (Figure 2G, R2 = 0.106, p = 0.002). Root and stem samples grouped with organic, conventional agricultural practice did not form a distinct cluster (R2 = 0.04, p = 0.135 and R2 = 0.02, p = 0.59, for root and stem, respectively, Figure 2H,I). The LDA effect size (LEfSe) analysis demonstrated that Flavobacterium, Sphingomonadales, xanthomnadaccae, and Rhizobiom were enriched in the organic tomato root, while Bacteriodes were contained in the conventional tomato root. Enterobacteriaceae are enriched in the organic stem, and corynbacterium is enriched in the conventional stem (Figure 3A,B). Bacterial function prediction analyses for these dominant taxa showed that the top four abundant functional groups (accounting for >75% of total groups) included lipid metabolism and signal transduction. Little function discrepancy was found between the organic and conventional tomato endophytic bacteria; however, a greater relative abundance of energy metabolism functional group was observed in root samples (3.8%) than that in stem (0.8%) samples (Figure 3C,D). To explore the effect of organic greenhouse and niche selection, the crop phase is not the main factor of our concern. Thus, four co-occurrence networks were constructed (conventional system root endophytes, conventional system stem endophytes, organic system root endophytes, and organic system stem endophytes) merged by four plant crop phases (Figure 4A–D). The topological features of the co-occurrence network are shown in Figure 4E and Supporting Information: Table S2, and details concerning the bacteria classified as nodes in the connector and module hub area can be found in Supporting Information: Table S3. The taxonomic composition of the networks showed no difference among the four treatments, with the majority of nodes belonging to Proteobacteria, Firmicutes, Bacteroidetes, and Actinobacteria, respectively. Network complexity gradually decreased from organic greenhouse to conventional greenhouse. The average connectivity degree of roots and stems in the organic greenhouse was 30.172 and 6.215, while that in the conventional greenhouse decreased by 19.415 and 5.893. All the nodes with Zi ≥ 2.5 or Pi ≥ 0.62 were determined as the keystone species, that is, nodes in the area of connectors (0.24%), module hubs (0.65%), and network hubs (0%) played a crucial role in the co-occurrence networks. Among them, the connectors and module hubs of phyla Proteobacteria (e.g., Rhizobiaceae, Enterobacteriaceae, Sphingomonadaceae, and Pseudomonadaceae), Bacteroidetes (e.g., Sphingobacteriaceae), and Actinobacteria (e.g., Anaerolinea) were accounted for 85%, 1%, and 0.5%, respectively (see Supporting Information: Table S3). RDA analysis was performed to reveal that soil TN was considered the major driver, which led to a tomato endophytic bacterial community shift (Figure 4F). To further characterize the compartment niche selection effect on tomato microbiota, and the key functional bacteria for TN identified based on the network analysis, microbial networks were generated (Figure 5A,B). More root bacteria were beneficial for soil organic formation than the stem. Network complexity gradually decreased from root to stem (with an average degree of 5.59 in the root and 4.30 in the stem); (see Supporting Information: Table S4). Results also demonstrated that most of the nodes (99.11%) were in the peripheral region, and their connections were primarily connected to the nodes in the module, at the same time, nodes in the area of connectors (1.8%), module hubs (2.2%), and network hubs (0%) played a crucial role in the co-occurrence networks. Among them, the connectors and module hubs of phyla Proteobacteria (e.g., Rhizobiaceae, Enterobacteriaceae, Sphingomonadaceae, and Pseudomonadaceae), Bacteroidetes (e.g., Sphingobacteriaceae), and Bacteroides (e.g., Sphingobacteriaceae) were accounted for 77% and 21%, respectively (see Supporting Information: Table S3). We isolated a total of 663 bacteria by the cultural-dependent method, but after 16S sequence alignment to remove redundancy, we obtained 60 bacteria from a different genus. Results showed that the number of endophytic bacterial increased with a plant growing, ranging from 4.21 to 7.45 (log10 CFU/g.FW), reaching a peak in the fruiting phase (Figure 6A). Organic agricultural practice and conventional practice nearly almost have no influence on the number of endophytic bacteria, and the isolated number of bacteria was more in roots than that in stems (p < 0.01, Figure 6A). Bacteroidetes and Proteobacteria were dominant endophytic bacterial phylum, specifically, Bacillus is more prevalent in the organic farming greenhouse than that in the conventional farming greenhouse (Figure 6B). In this study, a greenhouse survey demonstrated that tomato endophytic microbiomes were assembled largely driven at the compartment niche perspective, but had insignificant influence from fertilization environmental practice. Numerous studies indicated that the plant microbiome is shaped largely by host selection, such as genetics and organ compartment niche, while environmental factors (i.e., fertilizer regimes) exerted limited influence [18]. The lower diversity in the endosphere than that in the rhizosphere and soil is likely due to the host selective effect [19]. Endophytic bacterial showed environmental insensitivity, while rhizosphere soil microbial communities fluctuated over the course of plant development. This may be due to changes in root exudate profiles having a more significant influence on rhizosphere soil microbiome over growth phases [20]. Han et al. [20, 21] have shown that the tomato bacterial communities in different compartment niches are significantly different, with host selection sequentially increased from soils to epiphytes to endophytes, and the successful colonization of microorganisms in different parts of plants requires the ability to overcome host immunity and abiotic pressures in different microhabitats. Similar results were observed in sugarcane [22], rice [23], and corn [24] associated microbiome, that is, compartment niche played a key role either in the composition of endophyte or epiphyte [25]. Xia et al. [26] revealed the relative contributions of host and environmental factors to maize [27], wheat [28], and barley [29] microbiome community construction, confirming that crop microbiome community construction was mainly shaped by site niche and host species, while location (and site-related soil and climate factors) or fertilization practices had little effect. Co-occurrence patterns of bacterial communities of different agricultural practices were revealed in this study, and organic regimes established a more complicated network than conventional practices, and root networks are more complex than stem networks in both agriculture practices. Furthermore, the analysis of network correlation to soil chemical properties indicated that the positive and negative effects of organic on bacterial community up to 50% of the time, whereas TN accounted for 40.4% and 59.6% in root and stem, respectively. Similar to the soil microbiome study findings by Owens, S. M et al. [30]. The TN accounted for the positive effect of almost 70% (Figure 5), highlighting the importance of TN in driving bacterial community composition and function [31, 32]. In addition, the KEGG categories (L2 level) indicated that the amino acid metabolism of the microbial community in the stem niche was significantly higher than that in the root, while the signal transduction, energy metabolism, and metabolism of cofactors and vitamins were lower than root, and these results may cause less impact of stem bacterial community in soil organic content compared with root. Moreover, the regulation of soil TN may be more important for remodeling the root microbial community [33]. Several biomarkers were identified at the genus level in different treatments based on Lefse analysis. Results showed that Bacillales were enriched in the organic tomato, while Enterobacteriales and Actinobacteria were enriched in the conventional tomato (LDA > 3.5, p < 0.05, Figure 3A,B). Surprisingly, more taxonomy discrepancies occurred between root and stem in this study (Supporting Information: Figure S4). Zhao et al. [34] found that plant roots recruited a higher proportion of Proteobacteria and Bacteroides, which played a key role in modulating nutrition intake [35], pathogen suppression [36], and plant tolerance to stress [37]. In addition, the KEGG metabolism pathway (L2 level) of compartment niche is insignificantly different among root and stem (Figure 3C). The tomato root microbiome responds to environmental stresses by regulating specific metabolic pathways, such as amino acid metabolism and carbohydrate metabolism functions, and the activity of signal transduction was significantly improved to increase the resistance to adverse environments [38] (Figure 3D). This study revealed that the organ compartment niche was a key driver shaping endophytic bacteria and root microbiome but not that in the stem. Cultivable isolation further demonstrated that most of the dominant taxa belonged to Proteobacteria and Firmicutes, and the majority of biocontrol microbiome belonged to this phylum, having antagonism against soil-borne diseases [39]. Comprehensive analysis of the relative contribution of compartment niche and environmental factors to the tomato endophytic microbiome assembly in this study. The findings revealed that compartment niche can drive community-scale shifts in tomato endophytic bacterial composition; however, the agricultural practice had little influence on endophytic bacterial communities. Organic practices remodeled a more complicated network than conventional practices, and root endophytic bacterial networks are more complex than that stem networks. The amino acid metabolism of the microbial community in the stem was significantly higher than that in the root, while the signal transduction, energy metabolism, and metabolism of cofactors and vitamins were lower than root. Proteobacteria and Firmicutes are dominant taxonomies in the cultivable endophytic bacteria, and Bacillus is more highly abundant in organic farming. In conclusion, this study provides a comprehensive view of compartment niche and agricultural practice factors shaping the tomato endophytic bacterial communities under greenhouse conditions. These efforts will provide important information for the incorporation of beneficial bacteria into tomato agricultural production. Long-term trials have been conducted in Quzhou county, Hebei province (36°52′N, 115°01′E) since 2002. Both the ORG and CON systems were part of this experiment. There were two semi-round arched greenhouses (52 m in length and 7 m in width for each greenhouse) as shown in Figure 1A. Crop cultivars, irrigation, and tillage schemes were the same in all two systems in the same growing season, and a previous study had documented the agricultural management in each greenhouse in detail [40]. ORG manipulation was carried out in accordance with the International Federation of Organic Agriculture Movements' (IFOAM) guidelines, employing compost obtained from chicken and cow dung, as well as biocontrol for plant protection (i.e., Sticky yellow paper traps and insect net to control pests, mechanical removal of diseased plants and ridging planting to control plant diseases, mechanical/mulch to suppress weeds, etc.) for plant protection. The CON system used fertilizers, insecticides, and poultry to produce greenhouse vegetables locally [41, 42]. Tomato samples (Solanum lycopersicum cultivar “Money Marker”) were collected from organic and conventional greenhouses at four-time points (i.e., seedling phase, flowering phase, fruiting phase, and harvesting phase), and both roots and stems parts were sampled. Each treatment had five replications; thus 40 root and 40 stem samples were sampled in total (Figure 1A). R2A agar plate (Reasoner's 2A agar) is a culture medium, which is a nutritionally rich media supporting the growth of fast-growing bacteria [43, 44]. Soil pH was measured using a pH meter and a soil-to-water ratio of 1:5 (PHS-3C). The elemental analyzer was used to measure soil chemical properties, such as TN and organic content (Vario EL III-Elementer), NH4+–N was analyzed by Flow Analysis Instrument (AutoAnalyzer3). Soil chemical properties were recorded by Excel (v.2019) and significance was tested by ANOVA in SPSS (v.19). Plant samples were syndicated in Diagenode Bioruptor at a low for physical removal of soil and and all samples were on further at was performed as bacterial growth was after the roots on the R2A agar plate at for 7 as well as no of the 16S was observed as and was microbial community from samples was using a for soil The bacterial 16S were with the and with were with using the of All sequence were to each based on and and the were for the The and sequences > without and average were used for further of the taxonomic was performed as of were performed with by the All sequences under the number in this study, and the were in the The microbiome composition was by using and the effect of different factors, including time, and plant on microbial community was tested by ADONIS using the in and correlation analyses were performed to the between microbial communities and environmental using the in more than of samples were for the network analysis to and the network topological parameters were in and analysis and a significance were used to explore the most using and respectively were analysis with a of relative to diseased change relative to diseased significance with p < was used to the tomato endophytic functional and KEGG functional was used to information on the of at each KEGG function was used for and using t to the between the two were with pH with for at The roots and stems were to be all were after at at and bacteria as as the roots and stems were up into and a of root of the plant including roots as well as root was collected and for bacteria All bacteria were isolated by and each isolate was by analysis, and sequence were using and the was identified and the the data and and performed the data analysis. and and and and All to the and the This study was by the of China We the and on the of this The no of All data analyzed during this study are included in this and Supporting The data in this study can be found in The of the and can be found and in the Supporting and may be found in the or The is not for the content or of supporting information by the than be to the for the

Topics & Concepts