Abstract
Microbiomes are important for crop performance. However, a deeper knowledge of crop-associated microbial communities is needed to harness beneficial host-microbe interactions. Here, by assessing the assembly and functions of maize microbiomes across soil types, climate zones, and genotypes, we found that the stem xylem selectively recruits highly conserved microbes dominated by Gammaproteobacteria. We showed that the proportion of bacterial taxa carrying the nitrogenase gene (nifH) was larger in stem xylem than in other organs such as root and leaf endosphere. Of the 25 core bacterial taxa identified in xylem sap, several isolated strains were confirmed to be active nitrogen-fixers or to assist with biological nitrogen fixation. On this basis, we established synthetic communities (SynComs) consisting of two core diazotrophs and two helpers. GFP-tagged strains and 15N isotopic dilution method demonstrated that these SynComs do thrive and contribute, through biological nitrogen fixation, 11.8% of the total N accumulated in maize stems. These core taxa in xylem sap represent an untapped resource that can be exploited to increase crop productivity.
Similar content being viewed by others
Introduction
Maize (Zea Mays L.) is grown worldwide and is of great importance for food production. In 2020, the global production of maize was estimated at 1.05 million thousand tons1, ranking first in terms of yield production2. Several studies have shown that microbiomes have essential functions in maize performance. For example, field experiments in nitrogen (N)-depleted soil in Mexico indicated that the mucilage of aerial roots of a maize landrace was enriched with diazotrophs, and their fixation of atmospheric N contributed 29–82% of the N nutrition of maize3. Similar to the human microbiome4, the plant microbiome is referred to as the second genome of the plant, and is crucial for plant health5. Thus, targeted manipulation of crop microbiomes represents a promising approach to maximise sustainable crop production in the future6.
It is now widely accepted that plants and microorganisms can form complex co-associations governed by specific assembly rules7. As such, the recruitment and the selection of host-adapted microorganisms is of importance for crop health and nutrition7,8. It has even been suggested that plants and their associated microbiome are collectively forming holobionts and therefore should no longer be considered as standalone entities9. Recent studies have reported that a “core microbiota” exists within the plant, i.e., a subset of the plant microbiota that is reproducibly associated with a particular crop species across a wide range of scales7,10. These core microorganisms have a persistent relationship with the host, irrespective of the environmental conditions11. Since the core microbiota has traits for efficient colonisation, nutrient acquisition, and stress tolerance12,13, it is likely to be particularly important for the biological functions of the host and may, therefore, affect crop health and nutrition. In recent years, the development and construction of synthetic communities (SynComs) has provided functional and mechanistic insights into microbe-host relationships and how these relationships influence plant fitness7. For example, inoculation of maize seedlings with a simplified and representative seven-strain SynCom resulted in a stronger biocontrol effect against the phytopathogenic fungus Fusarium verticillioides than inoculation with each strain separately, indicating a clear benefit to the host14.
Previous studies have mainly characterised the core microbiota in the belowground compartments of plants15,16, the phyllosphere17, and seeds18. For example, Lundberg et al. initially identified core taxa in Arabidopsis thaliana roots as Actinobacteria and other specific families15. Grady et al. identified seasonally persistent and dynamic core phyllosphere microbiotas in perennial biofuel crops17. However, most studies on the assembly of the core microbiome and its collective functions in terms of host health and nutrition have overlooked the other plant compartments, such as vascular tissues. The vascular tissue acts as an effective long-distance transport system, which is driven by hydrostatic pressure gradients between the root and the shoot19. This driving force ensures smooth transport of solutes and signals among plant organs20. Moreover, the sizes of the holes in the perforated plates between xylem elements are large enough to allow passage of bacteria21,22. A few studies have reported that some systemic bacterial colonisers can spread to aboveground plant compartments through transpiration-driven xylem flow22. Anguita-Maeso et al. characterized the bacterial communities inhabiting the olive xylem sap by culture-dependent and independent approaches and found that Sphingomonas was the most representative genera23. Nevertheless, little is known about the functional relationship between the xylem-inhabiting microbiota and plant growth and development.
Here, we determined whether maize has core microbiomes in different plant compartments that can perform functions increasing host fitness. For this purpose, we investigated the microbiome of seven plant compartments in maize genotypes grown across a range of environmental conditions (i.e., different fertilisation regimes, soil types, climate zones). We found that the microbiota was more conserved in the xylem across soils than in the other plant compartments. We identified a core microbiota of 25 operational taxonomic units (OTUs) that was consistently detected in all xylem sap samples. Strains isolated from this core microbiota were capable of biological N-fixation (BNF). Their endophytic behaviour and contribution to maize N nutrition was verified by analyses of GFP-tagged strains and a 15N isotopic dilution method. Our study reveals details of an overlooked core microbiome inhabiting the xylem sap and how it contributes to the N nutrition of maize.
Results
Effects of soil type and fertilisation on the maize microbiome
To detect differences in the maize microbiota among various plant compartments along the soil–maize continuum, we sampled plants from six long-term fertilisation experimental sites with three soil types (black soils, BSA; red soils, RSA; fluvo-aquic soils, FSA) across the middle temperate zone to the subtropical zone of China (Supplementary Data 1 and Supplementary Fig. 1). The soils in these areas have been fertilised for at least 29 years according to three regimes (no fertilisation, Control; chemical fertiliser N, phosphorus, and potassium, NPK; and organic manure plus chemical fertiliser, NPKM), which has led to drastic differences in soil fertility levels and bacterial communities (Supplementary Table 1 and Supplementary Fig. 2). Analyses of bacterial communities from the bulk soil (BS), rhizosphere soil (RS), root endosphere (RE), xylem sap (VE), stem endosphere (SE), leaf endosphere (LE), and phylloplane (P) by 16 S rRNA amplicon sequencing of the V5–V7 region suggested that plant compartments were the main driver of maize microbiome composition (plant compartments, R2 = 57.72%, P < 0.001; sites, R2 = 7.20%, P < 0.001; fertilisation treatments, R2 = 0.53%, P < 0.001, Fig. 1a). Moreover, an unconstrained principal coordinates analysis (PCoA) of data from the BS, RS, RE, VE, SE, LE, P samples showed that the beta-diversity of aboveground bacterial microbiomes differed from that of the root and soil microbiomes (Supplementary Fig. 2 and Supplementary Fig. 3). This trend was further captured by an independent permutational multivariate analysis of variance (PerMANOVA), which showed a gradually decreasing influence of site and fertilisation on bacterial communities from the soil to xylem sap (from 91.93% in BS to 28.42% in VE) (Fig. 1b; Supplementary Table 2). Alpha-diversity analyses also revealed significantly lower diversity of the microbiota in VE (PFDR < 0.05) compared with that in BS and other plant compartments (Fig. 1c and Supplementary Fig. 4).
To further explore the influence of the environment on the plant microbiome, distance-decay relationships (DDRs) were analysed. We observed stronger geographic, edaphic, and climatic distance-decay of community similarity for RS and BS (Fig. 1d–f; Supplementary Table 3) than for the other compartments. However, such relationships were not observed in VE across geographic distance (R2 < 0.001; P = 0.601) and climatic distance (R2 < 0.001; P = 0.237). Non-significant correlations between the microbiome in VE and most soil physicochemical properties (e.g., total carbon, available phosphorus, soil water content) were confirmed by the Mantel test (Supplementary Fig. 5 and Supplementary Table 4). Together, these results indicate that, despite large edaphic and climatic differences across sites, the specific recruitment of bacterial taxa inhabiting VE was robust across genotypes and different environmental conditions.
Changes in abundance and composition of maize-associated microbiomes along the soil-plant continuum
We investigated the bacterial community composition along the soil–plant continuum, in particular, the steady-state composition of bacterial microbiomes in VE. We found that maize microbiomes in BS and RS consisted mostly of Acidobacteria, Actinobacteria, Gammaproteobacteria, Alphaproteobacteria, and Chloroflexi, accounting for 69.42%–85.89% of total relative abundance depending on soil types (Fig. 2a and Supplementary Data 3). Gammaproteobacteria were more abundant in RE than in BS and RS, and were predominant in stem microbiota (VE and SE). The average relative abundance of Gammaproteobacteria varied from 12.5% (BS), 29.2% (RS), 59.3% (RE), to 93.2% (VE), indicating that the transition in relative abundance profiles of host-associated lineages among plant compartments was more gradual than discrete.
Next, we used fast expectation-maximisation microbial source tracking (FEAST) to track the potential sources of plant bacterial microbiomes (Fig. 2b). The host showed a strong selection effect on the compositional pattern of microbial communities in habitats from the surface to the interior (Fig. 2b and Supplementary Fig. 6a). For example, 60.61% of the microbiome in RE was derived from RS. The VE acted as a linker between aboveground and belowground parts, with 24.21% of its bacterial microbiome originating from RE.
We further observed that the Gammaproteobacterial OTUs enriched in VE were dominated by Enterobacteriaceae, Erwiniaceae, and Pseudomonadaceae, while Burkholderiaceae was predominant in the SE population (Supplementary Fig. 6b, c). We detected 1364 OTUs that were shared among all plant compartments as well as the bulk soil, suggesting that they had the capacity to adapt to very different niches and could thrive even in plant organs (Fig. 2c). Notably, this set of taxa accounted for 76% of the total OTU numbers in VE and 99% of the total relative abundance. Of these, 166 OTUs were enriched in VE compared with bulk soil, and they mainly belonged to the genera Klebsiella and Pantoea (Fig. 2d, e and Supplementary Data 4). A linear model analysis showed that the dominant Klebsiella in VE was more enriched in control soils (Fig. 2f), with a relative abundance up to 43.35%, than in fertilised soils (relative abundance of 22.95% in NPK and 21.83% in NPKM) (Fig. 2g). These results suggest that most taxa enriched in xylem likely originated from the soil, gradually moved into different plant compartments, and eventually thrived in the stem xylem, where they may play a role in maize nutrition.
Xylem-inhabiting bacterial communities enriched functional roles related to N cycling
Given the key role of the stem xylem as the channel for materials transport between the aboveground and belowground compartments, we performed a functional annotation of prokaryotic taxa (FAPROTAX) analysis to generate putative functional profiles of the microbiotas inhabiting xylem and other compartments based on their community composition (Fig. 3a). The microbiotas in the aboveground and belowground compartments were separated into two clear clusters in terms of their putative functional profiles. Most functions were more abundant in root-associated environments than in aboveground compartments, for example, cellulolysis, predatory/exoparasitic, chitinolysis, and ureolysis. However, potential functions related to fermentation and N-cycling were overrepresented in VE. Specifically, the number of OTUs involved in N-cycling was approximately 2-times higher in VE than in LE, and 1.5-times higher than in SE (Fig. 3b). The number of these OTUs in VE was also significantly (PFDR < 0.05) higher in control soils than in the two long-term fertilised soils (NPK and NPKM), suggesting that these plant-recruited bacteria were more likely to play a critical role in N availability in plants under soil nutrient stress (Fig. 3b). However, like other tools assigning ecological functions based on 16 S data, FAPROTAX has some limitations related to the number of reference strains and the affiliation of functions that are poorly conserved; therefore, our FAPROTAX prediction needs to be confirmed by other methods.
To test this hypothesis, we performed triple-qPCR (see Methods) to determine the ratio of the nitrogenase gene (nifH) to the 16 S rRNA gene, as a proxy of the relative N-fixation potential of the bacterial community. As expected, this ratio was significantly higher (PFDR < 0.05) in VE than in the other three aboveground compartments (SE, LE, and P) in control soil only (Fig. 3c).
Identification and isolation of xylem-inhabiting core bacterial communities
Abundance–occupancy analyses identified 25 OTUs (occupancy of 1) that were persistent in xylem in all samples (Fig. 4a and Supplementary Data 5). These taxa were dominated by Gammaproteobacteria (89.42%), Firmicutes (4.10%), Alphaproteobacteria (0.18%), and Actinobacteria (0.13%), with an overrepresented OTU belonging to the Klebsiella genus (29.2%) (Supplementary Fig. 7). Their ubiquitous presence and high abundance indicated that these taxa were core members of maize-associated bacterial communities. Next, we isolated 969 pure bacterial isolates from the xylem sap using three media, representing a total of 109 unique strains belonging to four phyla and 23 families (Supplementary Data 6). Further characterisation of these isolates revealed that nine out of 109 strains grew on N-free medium, harboured the nifH gene, and showed high nitrogenase activity (Fig. 4b). We then cross-referenced the V5–V7 region of the 16 S rRNA gene of the 109 isolated strains against the highly abundant OTUs (> 0.01%) in VE. We found that these isolated strains shared sequence identity with 14 of the 25 core OTUs in VE, which accounted for 81.33% of the total bacterial abundance in VE (Fig. 4c, Supplementary Fig. 7 and Supplementary Data 7). Further, two of the 14 strains displayed remarkably high nitrogenase activity (Klebsiella variicola MNAZ1050 and Citrobacter sp. MNAZ1397, their combination hereafter referred to KC), and harboured the entire nif gene cluster in their genomes (Fig. 4d).
To explore whether the non-N-fixers in the core VE microbiota can assist in biological N-fixation, we mixed 12 non-N-fixers with diazotrophs and determined the total nitrogenase activity of each combination. We found that, albeit to different extents, the nitrogenase activity of each combination was higher than that of the diazotrophs alone, especially for the combinations of KC + Acinetobacter sp. ACZLY512 and KC + Rosenbergiella epipactidis YCCK550 (Fig. 4e). By sequencing the genomes of these non-N-fixing core strains, we found that the total reads of genes encoding major enzymes involved in the tricarboxylic acid (TCA) cycle was higher in Acinetobacter sp. ACZLY512 and R. epipactidis YCCK550 than in other strains (Supplementary Fig. 8a and Supplementary Data 8). Therefore, these non-N-fixers may have strong respiratory metabolism, thereby creating a micro-oxygen environment conducive to biological N-fixation. This was further confirmed by acetylene reduction assays (ARA) under initial atmospheric oxygen levels (see Methods), which showed that the nitrogenase activity of the diazotrophs was no longer inhibited in the presence of air after adding Acinetobacter sp. ACZLY512 (Supplementary Fig. 8b). We also found that there were more functional gene reads involved in ammonia assimilation in the two non-N-fixers (Supplementary Fig. 8c and Supplementary Data 8). This implied that they might facilitate N-fixation, because the activation and transcription of nitrogenase genes is positively regulated by ammonia comsumption24,25.
Successful colonisation and N-fixation by synthetic communities in xylem vessels
We further verified the endophytic behaviour and N-fixation function of the core strains by constructing and analysing synthetic communities (SynComs). The SynComs consisted of two N-fixing bacteria (K. variicola MNAZ1050 and Citrobacter sp. MNAZ1397) and two non-N-fixing strains that may help to fix N (Acinetobacter sp. ACZLY512 and R. epipactidis YCCK550), which together accounted for 34.61% of the total OTU abundance in the VE bacterial microbiota (Fig. 4e). At 63 days after inoculation onto maize plants, the successful colonisation of soil-derived SynComs in the maize aboveground compartments was detected. As shown in the confocal laser scanning microscope (CLSM) images, only the inoculated maize plant had bright green fluorescent signals from the GFP-tagged SynComs within the stem and xylem sap (Fig. 5a and Supplementary Fig. 9).
The contribution of biological N-fixation of the inoculated SynComs to maize N-nutrition was measured by the 15N isotopic dilution method. The maize plants inoculated with the SynComs had a higher proportion of N from the air (lower δ15N values) than did the uninoculated plants. The contribution of BNF to total accumulated N in maize roots, stems, and leaves was 4.01%, 11.81%, and 0.46%, respectively (Table 1). Characterisation of maize plant phenotypes revealed that, compared with uninoculated control plants, those inoculated with the SynComs showed significantly higher root dry weight (Fig. 5b, c and Supplementary Fig. 9) and root-to-shoot ratio (Fig. 5d, independent samples T-test, P < 0.01).
Discussion
The conductive property of plant xylem vessels makes it a highly dynamic endosphere21. Xylem vessels with open perforation plates function as a free-flowing pipe, thereby allowing endophytes to spread longitudinally from the rhizosphere to other plant tissues22,26. However, the importance of the xylem in the plant microbiome assembly has been overlooked. In the present study, while the root- and leaf-associated bacterial communities varied significantly with environmental factors as previously reported27,28, the stem xylem recruited highly conserved taxa independent of geographic (P = 0.601) and climatic (P = 0.237) distances. A FEAST analysis showed that 60.61% of bacteria in the root endosphere were derived from rhizosphere soils; 24.21% of bacteria in the stem xylem were derived from the root endosphere and 9.15% of bacteria in the leaf endosphere were derived from xylem-inhabiting endophytes (Fig. 2b). It can therefore be assumed that these facultative endophytes positively selected by the host have the physiological plasticity to thrive in different plant compartments, with functionally important roles29,30. Accordingly, the proportion of bacterial taxa harbouring N-fixation genes was larger in the xylem than in other plant compartments. For example, the ratio of nifH to 16 S rRNA in xylem was about twice that in the root endosphere and four times that in the leaf endosphere, especially in control soils where this function is likely to be more important because plants are under nutrient stress. This is also confirmed by the fact that the average relative abundances of two core N-fixers in SynCom (K. variicola MNAZ1050 and Citrobacter sp. MNAZ1397) were significantly higher in the control soil than in two fertilised soils (NPK and NPKM) (Supplementary Fig. 10). Therefore, our results confirm those of previous studies, and suggest that the beneficial effects of endophytic bacteria could be as important as those of the microorganisms present in the rhizosphere29,31.
We further conducted an abundance–occupancy pattern analysis to identify the core subset of bacterial communities in xylem. This approach provides ecological insights into which taxa are most important for the plant32. We used a conservative threshold for the occupancy analysis, so that members entering the xylem due to chance events such as wounds and cracks would be excluded. Both the occupancy and abundance patterns indicated the presence of a core microbiome consisting of 25 OTUs with high abundance, consistent with previous studies that identified core microorganisms across a wide of sites10,33. The fact that the samples were collected from very diverse sources (six different long-term fertilisation experiments, from sites spanning three climate zones, and encompassing multiple maize genotypes) provides further evidence for the persistent association between the core subset and its host. The isolation of endophytic diazotrophic bacteria from root, stem, and leaf tissues of maize cultivars has been demonstrated in a previous study34. Here, concurrent culture-dependent methods yielded isolates that were cross-referenced against the core OTUs. Among these isolates, one strain mapping to the most abundant core OTU (29.2%) was identified as K. variicola and contained the entire nif gene cluster with nitrogenase activity of up to 1180 nmol C2H4 (mg protein hr)−1. A previous study reported that Klebsiella is present as an endophyte in maize stem and root tissue, and quantified its dinitrogenase reductase activity inside maize tissue after adding an exogenous carbon source35. Our results indicate that K. variicola is dominant within the maize core microbiota across different soil types, climate conditions, and maize genotypes, suggesting that its key beneficial role in maize might be more general. Furthermore, we found that the microbial communities in xylem include core non-N-fixers such as Acinetobacter sp. ACZLY512 and R. epipactidis YCCK550 with the potential to assist in N-fixation. The genome analysis of these strains provided a possible explanation for this mutually beneficial symbiosis, with both genomes showing a higher abundance of gene clusters potentially involved in the TCA cycle and ammonia assimilation. Bacteria with the ability to perform BNF are known to be sensitive to the concentrations of fixed N and external oxygen24. Thus, we speculate that these core non-N-fixers might assist in N-fixation by regulating N and oxygen concentrations in the xylem microenvironment where the diazotrophs are located.
Finally, we established SynComs consisting of two diazotrophs and two helper strains. Our results confirmed the endophytic lifestyle and N-fixation capacity of the SynComs in planta. The contribution of BNF to total N accumulation was higher in maize stems (11.81%) than in roots (4.01%) and leaves (0.46%), demonstrating that SynComs do thrive and contribute to plant N-nutrition in the stem xylem after being recruited from the soil. Interestingly, the SynComs also increased root dry biomass (Fig. 5b, c), suggesting that there is a multi-mechanism synergistic interaction between the core endophytes and their host. This is supported by the fact that the genomes of K. variicola MNAZ1050 and Acinetobacter sp. ACZLY512 contain genes involved in IAA biosynthesis and ethylene biosynthesis, respectively. Thus, hormones produced by these strains may promote root growth. Such synergistic effects of the bacterial consortium were previously observed in other plant-bacterial interactions. For example, Santhanam et al. showed that a five-strain SynCom was able to effectively protect Nicotiana attenuata against fungal sudden-wilt disease via complementary traits, while the individual strains or smaller combinations of them did not36. Overall, our results highlight the limitation of incorporating only the desirable traits into the design of SynComs since this method might overlook helper strains that could improve the overall beneficial activity of the SynComs.
In summary, our data indicate that maize specifically recruit a core microbiota in the xylem sap that is conserved across environmental conditions as well as genotypes. This core microbiota contributes to plant N-nutrition through BNF and promotes root development. Further research is needed to explore the mechanisms and functions involved in the association between the xylem sap microbiota and its host. This core microbiota could represent a promising resource to develop alternative microbial biotechnologies to enhance crop performance in sustainable agriculture.
Methods
Site description
The six long-term fertilisation field experiments were located across a latitudinal gradient in China from north to south, spanning three climate zones from the middle temperate zone to the subtropical zone. These sites were chosen to represent the three main agricultural production areas in China. The soils at the sites were black soils (Udolls or Typic Hapludoll according to USDA Soil Taxonomy, BSA) in Hailun (Heilongjiang Province) and Changchun (Jilin Province); fluvo-aquic soils (Aquic Inceptisol according to USDA Soil Taxonomy, FSA) in Yucheng (Shandong Province) and Yuanyang (Henan Province); and red soils (Ultisols according to USDA Soil Taxonomy, RSA) in Jinxian (Jiangxi Province) and Qiyang (Hunan Province). The two most distant sites (from Hailun to Qiyang) were more than 2,500 km apart, and the two closest test sites (from Yucheng to Yuanyang) were at least 300 km apart. Three treatments, i.e., no fertiliser (Control), chemical fertiliser N, P, and K (NPK), and organic manure plus chemical fertiliser (NPKM), have been applied in triplicate plots in each field for 29 years or more. Climate data corresponding to the sampling site coordinates were obtained from the China Meteorological Data Network (http://data.cma.cn/). Further details of experimental sites are provided in Supplementary Data 1.
Sample collection
Sampling was performed during the silking-maturity period of maize in 2019 and 2020, with the exact date of sampling depending on the developmental stage of plants at each location (Supplementary Data 2). In the FSA and RSA soils, three individual maize plants were selected from each subplot (a total of 27 maize plants per site). From each maize plant, we collected the following compartments in the field: mixed leaves, xylem sap, stem, roots, bulk soil. The mixed leaves sample consisted of the 2nd, 4th, and 6th leaves, which were removed from the plant stem using ethanol-sterilised scissors. To collect xylem sap, a proxy for xylem, we cut off the stem mid-way between the 2nd and 3rd node from the base of the plant, and sterilised absorbent cotton in sterilised bags was placed on the cut end of the shoot (see Supplementary Movie 1 for details of this procedure). Meanwhile, we inserted a steel stick (sterilised, 2 cm diameter, 30 cm length) into the soil at each subplot to simulate the collection of xylem sap and check for contaminants during the field operations (Supplementary Fig. 11). The stem sample consisted of the upper region between the 2nd and 3rd nodes, collected into sterilised bags. To collect root samples, we shook whole roots vigorously to remove all loose soil. The roots and root-adhered soil particles were collected for further separation of the roots and rhizosphere soil in the laboratory. The bulk soil sample was collected from between the rows of maize plants. At the sites with BSA soils, we only collected one individual maize plant from each subplot (a total of nine maize plants per site), because these two long-term experiments have strict requirements for sampling to avoid large-scale damage to the entire test field. Thus, for each plant, the 1st and 2nd, the 3rd and 4th, and the 5th and 6th leaves were collected as three replicates. Similarly, fine roots (< 2 mm) and thick roots (> 2 mm) were collected separately. Except for this difference, the other operations were the same as those at the other four experimental sites. All samples were placed on ice for transport and further processing within 48 h. The soil parameters of pH, total C (TC) and N (TN), ammonium (NH4+) and nitrate (NO3−), and soil available P (AP) and K (AK) are listed in Supplementary Table 1.
Sample processing
To recover the xylem sap absorbed in the cotton, each cotton ball was placed into a 50-mL sterile centrifuge tube with a filter and centrifuged at 6000 × g for 5 min. The collected sap was divided into two parts; one part was used for bacterial isolation, and the other part (stored at −80 °C) was used for culture-independent bacterial 16 S rRNA gene profiling.
Processing of root-associated samples
We used a modified protocol15 to separate the microbiome living on the plant surface (epiphytes) from the microbiome living within the plant (endophytes). Briefly, 5 g root tissue was weighed into a 100 mL conical flask containing 80 mL sterile PBS and 5 μL Tween 80. The mixture was vortexed, and the liquid was collected as the rhizosphere (root epiphyte) sample. To extract rhizosphere DNA, the sample was centrifuged at 10,000 × g for 5 min, and then 500 mg of the resulting tight pellet containing fine sediment and microorganisms was placed in a Lysing Matrix E tube (supplied in the FastDNA™ Spin Kit for Soil). To obtain endophytes from the root samples, the roots were washed with fresh PBS until the buffer was clear after vortexing. The roots were then sonicated using an ultrasonic cell disruptor (Scientz JY 88- IIN, Ningbo Scientz Biotechnology Co., Ltd., Zhejiang, China) at a low frequency for 10 min (30-s bursts followed by 30-s rests).
Processing of leaf-associated samples
The method for washing leaf samples was similar to that used to wash the roots, except for an additional step before collecting the phyllosphere. Each leaf sample (5 g) was added to a conical flask containing sterile PBS, which was subjected to two 5-min treatments in an ultrasound bath (25 °C, 40 KHZ), with vortexing for 30 s between the two ultrasonication treatments. This procedure released most of the phyllosphere microbes from the leaves. Then, the filtrate was collected and the leaves were washed again using the above procedure. Phyllosphere samples were collected by centrifugation (at 10,000 × g for 20 min) of the accumulated filtrate and were resuspended in 1 mL sodium phosphate buffer (FastDNA™ Spin Kit for Soil) before being transferred to Lysing Matrix E tubes (FastDNA™ Spin Kit for Soil). Finally, the leaf and stem samples were washed and sonicated in the same way as the roots. Sonicated root, leaf, and stem samples were snap-frozen in liquid N2 and stored at −80 °C until analysis.
DNA extraction, PCR amplification and sequencing
Total DNAs were extracted from the aforementioned samples with a FastDNA™ Spin Kit for Soil (MP Biomedicals, Solon, OH, USA) following the manufacturer’s instructions. The DNA concentration and purity were measured using a NanoDrop2000 spectrophotometer (NanoDrop2000, Thermo Fisher Scientific, Waltham, MA, USA). The DNAs extracted from soils (bulk and rhizosphere soil) were diluted 10-fold. For 16 S rRNA gene libraries, the V5–V7 region was amplified using the primers 799 F and 1193 R (Supplementary Table 5). Each DNA template was amplified in triplicate (together with a water control) in a 25-μL reaction volume. The PCR conditions were as follows: 12.5 µL 2× EasyTaq PCR SuperMix (TransGen Biotech, Beijing, China), 1.25 µL forward primers (10 µM), 1.25 µL barcoded reverse primers (10 µM), 1.25 µL template DNA, and 8.75 µL ddH2O. The PCR amplification program was as follows: 94 °C for 3 min; 28 cycles of 94 °C for 30 s, 55 °C for 30 s, 72 °C for 90 s; and 72 °C for 90 s. The products were stored at 4 °C until use. After mixing the triplicate PCR products of each sample, the bacterial 16 S rRNA gene amplicons were extracted from a 1% agarose gel using a Gel Extraction Kit (Omega Bio-tek Inc., Norcross, GA, USA). The DNAs were measured using a Quant-iT™ PicoGreen™ dsDNA Assay kit (Thermo Fisher Scientific) and pooled in equimolar concentrations. Sequencing libraries were generated using an Illumina TruSeq DNA PCR-Free Library Preparation Kit (Illumina, San Diego, CA, USA) and were sequenced on the NovaSeq-PE 250 platform (Illumina).
16 S rRNA gene amplicon sequence processing
Paired-end reads were checked by FastQC v.0.10.137 and merged using the USEARCH 11.0.66738 fastq_mergepairs script. Reads were assigned and demultiplexed to each sample according to the unique barcodes by QIIME 1.9.139. After removing barcodes and primers, low-quality reads were filtered and non-redundant reads were identified using VSEARCH 2.12.040. Unique reads with ≥ 97% similarity were assigned to the same OTU. Representative sequences were selected using UPARSE41 and the classify.seqs command in mothur42 was used to taxonomically classify each OTU with reference to the SILVA 138 database43. The OTUs classified as host plastids, cyanobacteria, and others not present in our samples were removed from the dataset.
Statistical analysis
Alpha and beta-diversity analyses were conducted with R script as described in previous studies44 and on the QIIME239 platform. PerMANOVA analyses were performed using the ‘Adonis’ function implemented in the vegan package45 of R. We used the microbial source-tracking method FEAST46 to determine the potential origin of the microbiota inhabiting the various compartments of maize plants. Bacterial OTUs were assigned into multiple functional groups using FAPROTAX v.1.2.147. All analyses were conducted in the R Environment48 except for beta-diversity analyses. All plots were generated with ggplot249 and GraphPad Prism 8.0.0 (GraphPad Software, San Diego, CA, USA, www.graphpad.com).
Distance-decay relationship analyses
We conducted distance-decay relationship analyses to assess the relationship between the similarity of communities in individual plant compartments and spatial distance, edaphic distance, and climatic distance. We used the Geosphere package to calculate the geographic distances in km from the latitude and longitude coordinates, and calculated edaphic and climatic distances separately as the Euclidian distance. We used the ‘vegdist’ function in the vegan package to calculate the Bray–Curtis similarity of microbial communities. Here, the variation in the slope of distance decay reflects the degree to which the similarity of microbial communities in plant compartments varies with environmental distances. The relationships between the Bray–Curtis similarity of each compartment and specific soil or climatic factors were determined by calculating Pearson’s correlation (r) values. The significance of r values was assessed with the Mantel test implemented in the vegan package.
Differential abundance testing
A negative binomial generalised linear model was implemented with the edgeR package50 to detect differences in OTU abundance among samples. We compared individual plant compartments (RS, RE, VE, SE, LE, and P) against bulk soil (BS) and conducted pairwise comparisons among fertilisation treatments. For each comparison, after constructing the DGEList object and filtering out low counts, the calcNormFactors function was used to obtain normalisation factors and the estimateDisp function was used to estimate tagwise, common, and trended dispersions. We then used the glmFit function to test the differential OTU abundance. The corresponding P values were corrected for multiple tests using FDR with α = 0.05.
Core taxa selection
We used the UpSetR package51 to visualise the OTUs that overlapped among all plant compartments and soils. The overlapping OTUs were defined as those detected in at least one sample from each compartment. We further identified the union of overlapped OTUs and enriched OTUs in the xylem using the EVenn online tool52. Importantly, abundance–occupancy analyses were conducted to identify the core OTUs across environmental gradients. We calculated occupancy with the most conservative approach, which restricted the core to only those OTUs that were detected in all xylem sap samples (i.e., occupancy=1).
Triple-qPCR to verify FAPROTAX results
To detect contamination with plant organelles, sample DNAs were amplified using the universal primer pair 799 F/1193R53, which amplifies both bacterial 16 S and mitochondrial 18 S rRNA but not chloroplast sequences; the mitochondrial-specific primer pair mito1345F/mito4130R54,55, which only amplifies mitochondrial 18 S rRNA, and the nifH gene primers PolF/PolR56 (Supplementary Table 5). To reflect the potential N-fixation of bacterial communities in each plant compartment, the following ratio was calculated:
All qPCR assays were run on a QuantStudio 6 Flex using SYBR Green Pro Taq HS Premix (Accurate Biotechnology, Changsha, China) in a 20 µL volume containing 200 nM of each primer and approximately 50 ng DNA per reaction. All three primer sets were amplified in a three-step qPCR57 run at 95 °C for 15 s, 55 °C for 30 s and 72 °C for 40 s for 40 cycles followed by a melting curve analysis. The amplification efficiency varied from 83 to 117%.
Isolation of bacteria from xylem sap
To isolate strains, four gradient dilutions (10−3, 10−4, 10−5, and 10−6) of xylem sap were incubated on TSB, R2A, and Ashby’s Nitrogen-Free Agar media for 5–7 days at 30 °C (Supplementary Data 9). After incubation, colonies were selected based on their character and colony morphology and were purified by triple serial colony isolation. The isolates were subjected to Sanger sequencing and identified on the basis of PCR analyses with 27 F and 1492 R primers, and alignment against reference. 16S rRNA gene sequences using the BLAST algorithm. Isolates belonging to the core taxa were identified by comparing the 16 S rRNA V5–V7 regions against the highly abundant OTUs (> 0.01%) using UCLUST with 98.65% similarity;58 this threshold has been reported to accurately distinguish two species. Cladograms were visualised by iTOL.v6.459. The isolated cultures were stored in 30% (v/v) glycerol.
Nitrogen-fixing capacity of bacterial isolates
We used three different methods to evaluate the N-fixing capacity of bacterial isolates. First, we observed growth on Ashby’s N-Free medium, and documented which strains grew well after streaking of diluted cultures. Then, these strains were analysed by PCR to detect nifH with the PolF/PolR primer set56. The positive control was the N-fixing strain, Azotobacter chroococcum ACCC10006 (Agricultural Culture Collection of China). Strains that did not yield a PCR product with this primer set were analysed using other nitrogenase gene primers including nifH-F/nifH-R60 primers and the nested PCR primers FGPH19/PolR (outer primers) and PolF/AQER (inner primers)56 (Supplementary Table 5). We also conducted acetylene reduction assays (ARA)61,62,63 to quantify the nitrogenase activity of putative N-fixing strains. Each tested strain was initially incubated overnight in TSB medium and then washed twice with sterile 0.9% NaCl solution. After centrifugation and re-suspension, the bacterial pellet was added to a 20 mL serum vial containing 5 mL Dobereiner’s N-free liquid medium (Supplementary Data 9), reaching a final OD600 of ~0.1. The vials were first flushed with argon to evacuate air, and then 1% and 10% of the headspace was replaced with pure and fresh O2 and C2H2, respectively. After incubation at 30 °C for 12 h, the gas phase was analysed with a gas chromatograph (Agilent Technologies 6890 N). Data are presented as mean values from five replicate cultures. To test the hypothesis that the core non-N-fixers might assist N-fixation by modifying the oxygen concentration, the nitrogenase activity was measured as described above except that the headspace atmosphere in the sealed vial was not adjusted to 1% O2 by flushing with argon gas so that the initial oxygen concentration was that found in ambient air.
Draft whole-genome sequencing of cross-referenced core strains
Isolated genomic DNA was extracted with a TIANamp Bacteria DNA Kit (Tiangen Biotech, Beijing, China). The purified genomic DNA was used to construct a sequencing library, which was generated using the NEB Next® Ultra™ DNA Library Prep Kit for Illumina (NEB, Beverly, MA, USA) following the manufacturer’s recommendations. Pooled libraries were sequenced on the NovaSeq-PE 150 platform. After trimming low-quality reads by fastq64, the clean reads were assembled into draft genomes (excluding contigs of < 300 bp) by SPAdes 3.13.165. Gene prediction and annotation were performed by NCBI PGAP66, and the putative genes were further annotated by searching against the eggNOG database67 by emaper. The functional mapping and analysis pipeline (FMAP 0.15)68 was used to align the filtered reads using BLAST against a KEGG Filtered UniProt69 reference cluster (e < 1e-5, identity > 70%) and to calculate the number of reads mapping to each KEGG Orthologous group (KO). Other data manipulation was performed using perl scripts developed in-house.
Potted plant experiment
Two potted plant experiments were performed to (1) reproduce the endophytic behaviour of isolates from xylem sap; and (2) verify their N-fixation potential in maize. We constructed SynComs consisting of two diazotrophs (K. variicola MNAZ1050 and Citrobacter sp. MNAZ1397) and two non-N-fixers (Acinetobacter sp. ACZLY512 and R. epipactidis YCCK550) based on a N-fixing capacity test. Each individual strain was cultured overnight in TSB medium at 30 °C and 180 rpm, then cells were collected by centrifugation and the pellet was suspended in sterile 0.9% NaCl solution. Four bacterial suspensions were mixed in equal amounts to a final OD600 of ~0.2.
The potted plants were grown in plastic pots filled with loose a soilless mixture consisting of perlite and vermiculite (sterilised by autoclaving). The nutrients needed for plant growth were added as base fertilisers (Supplementary Data 10). The seeds of maize “Zhengdan 958” were surface-disinfected for 15 min with sodium hypochlorite (approximately 2% active chlorine, with 200 μL Tween 80) and washed for 5 min, five times, with sterile water. The final rinse water (100 μL) was spread on TSA medium to check for other attached bacteria. Seeds were allowed to germinate, and then germinated seeds with similar primary root lengths were selected for inoculation with SynComs. Maize plants were grown in a greenhouse under a 16-h light/8-h photoperiod at 30 °C/25 °C (day/night).
Colonisation of maize tissues by GFP-tagged SynComs
To verify the endophytic behaviour, each of the four members of SynComs was tagged with green fluorescent protein (GFP) (vector pCPP6529-GFPuv). One GFP-tagged and three other wild bacterial suspensions were mixed as described above, giving a total of four different GFP-tagged combinations. The control was SynComs with no GFP tags. For inoculation, maize seedlings with the endosperm removed were soaked in four GFP-tagged combined solutions for 30 min. The same bacterial suspension was applied to the potted plants at days 10 and 30 after transplanting. After 63 days, stem samples from between the 2nd and 3rd nodes were surface-sterilised with 70% ethanol, collected, and then sectioned using a Leica VT 1000 S vibratome (Leica, Nussloch, Germany). Thin sections (60 μm) and xylem sap were observed under a confocal laser scanning microscope (CLSM, Zeiss LSM 880 confocal microscope, Jena, Germany).
15N isotope dilution method
To verify the N-fixation potential of endophytes in maize, the N fertiliser was replaced with 15N-labelled (NH4)2SO4 (30 % 15N atom, Shanghai Research Institute of Chemical Industry, China). Maize seedlings with endosperm removed were soaked in a bacterial suspension of SynComs (Treatment group) or autoclaved SynComs (Control group) for 30 min. The same SynComs, active or autoclaved, were re-applied to the potted plants at 10 and 30 days after transplanting, as described above. The roots, stems and leaves were harvested separately from each treatment (four replicates) on day 65. The roots were washed with deionised water to remove adhering isotope residues. The N content and 15N enrichment of plant tissue were determined using Elementar vario PYRO cube elemental analyser (Vario PYRO Cube, Elementar, Hanau, Germany) and Isoprime 100 isotope mass spectrometer (Isoprime, Cheadle, United Kingdom). The plants inoculated with autoclaved SynComs were used as the reference to calculate BNF with the following equations34:
where %Ndfa is the percentage of N derived from air, %15Na.e. (%15N atom excess) is the enrichment in plants inoculated with SynComs (I) and autoclaved SynComs (UI), N2-fixed is N derived from air, and N yield.I is the total N content of the whole inoculated plant.
Reporting summary
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
Data availability
Raw 16 S rRNA gene amplicon sequence data have been deposited in the Sequence Read Archive under BioProject PRJNA795590. The Whole Genome Shotgun data has been deposited at DDBJ/ENA/GenBank under the accession PRJNA798255. The SILVA database is available at https://www.arb-silva.de. The eggNOG database is available at http://eggnog5.embl.de/download/eggnog_5.0/. KEGG Filtered UniProt reference cluster is available at https://ftp.uniprot.org/pub/databases/uniprot/uniref/uniref90/. Source data are provided with this paper.
Code availability
All scripts are available in GitHub (https://github.com/PlantNutrition/Liyu.)
References
World Data Atlas. Available online: https://knoema.com/atlas/topics/Agriculture/Crops-Production-Quantity-tonnes/Maize-production (accessed on 12 January 2022).
FAOSTAT-Crops Data. Available online: http://www.fao.org/faostat/en/#data/QC (accessed on 11 January 2022).
Van Deynze, A. et al. Nitrogen fixation in a landrace of maize is supported by a mucilage-associated diazotrophic microbiota. PLoS Biol. 16, e2006352 (2018).
Grice, E. A. & Segre, J. A. The human microbiome: our second genome. Annu. Rev. Genomics Hum. Genet. 13, 151–170 (2012).
Berendsen, R. L., Pieterse, C. M. & Bakker, P. A. The rhizosphere microbiome and plant health. Trends Plant Sci. 17, 478–486 (2012).
Singh, B. K., Trivedi, P., Egidi, E., Macdonald, C. A. & Delgado-Baquerizo, M. Crop microbiome and sustainable agriculture. Nat. Rev. Microbiol. 18, 601–602 (2020).
Trivedi, P., Leach, J. E., Tringe, S. G., Sa, T. & Singh, B. K. Plant-microbiome interactions: from community assembly to plant health. Nat. Rev. Microbiol. 18, 607–621 (2020).
Afzal, I., Shinwari, Z. K., Sikandar, S. & Shahzad, S. Plant beneficial endophytic bacteria: Mechanisms, diversity, host range and genetic determinants. Microbiol. Res. 221, 36–49 (2019).
Vandenkoornhuyse, P. et al. The importance of the microbiome of the plant holobiont. N. Phytol. 206, 1196–1206 (2015).
Risely, A. Applying the core microbiome to understand host-microbe systems. J. Anim. Ecol. 89, 1549–1558 (2020).
Stopnisek, N. & Shade, A. Persistent microbiome members in the common bean rhizosphere: an integrated analysis of space, time, and plant genotype. ISME J. 15, 2708–2722 (2021).
Tian, B. et al. Beneficial traits of bacterial endophytes belonging to the core communities of the tomato root microbiome. Agric. Ecosyst. Environ. 247, 149–156 (2017).
Hamonts, K. et al. Field study reveals core plant microbiota and relative importance of their drivers. Environ. Microbiol. 20, 124–140 (2018).
Niu, B., Paulson, J. N., Zheng, X. Q. & Kolter, R. Simplified and representative bacterial community of maize roots. Proc. Natl Acad. Sci. USA 114, E2450–E2459 (2017).
Lundberg, D. S. et al. Defining the core Arabidopsis thaliana root microbiome. Nature 488, 86–93 (2012).
Xu, J. et al. The structure and function of the global citrus rhizosphere microbiome. Nat. Commun. 9, 4894 (2018).
Grady, K. L., Sorensen, J. W., Stopnisek, N., Guittar, J. & Shade, A. Assembly and seasonality of core phyllosphere microbiota on perennial biofuel crops. Nat. Commun. 10, 4135 (2019).
Links, M. G. et al. Simultaneous profiling of seed-associated bacteria and fungi reveals antagonistic interactions between microorganisms within a shared epiphytic microbiome on Triticum and Brassica seeds. N. Phytol. 202, 542–553 (2014).
Fricke, W. Xylem: differentiation, water transport and ecology. (eLS. John Wiley & Sons, Ltd: Chichester 2017).
Lucas, W. J. et al. The plant vascular system: evolution, development and functions. J. Integr. Plant Biol. 55, 294–388 (2013).
Bove, J. M. & Garnier, M. Phloem-and xylem-restricted plant pathogenic bacteria. Plant Sci. 163, 1083–1098 (2002).
Compant, S. et al. Endophytic colonization of Vitis vinifera L. by plant growth promoting bacterium Burkholderia sp. strain PsJN. Appl. Environ. Microbiol. 71, 1685–1693 (2005).
Anguita-Maeso, M. et al. Culture-dependent and culture-independent characterization of the olive xylem microbiota: Effect of sap extraction methods. Front. Plant Sci. 10, 1708 (2019).
Dixon, R. & Kahn, D. Genetic regulation of biological nitrogen fixation. Nat. Rev. Microbiol. 2, 621–631 (2004).
Bueno Batista, M. & Dixon, R. Manipulating nitrogen regulation in diazotrophic bacteria for agronomic benefit. Biochem. Soc. Trans. 47, 603–614 (2019).
James, E. K. et al. Further observations on the interaction between sugar cane and Gluconacetobacter diazotrophicus under laboratory and greenhouse conditions. J. Exp. Bot. 52, 747–760 (2001).
Edwards, J. et al. Structure, variation, and assembly of the root-associated microbiomes of rice. Proc. Natl Acad. Sci. USA 112, E911–E920 (2015).
Long, J. et al. Environmental factors influencing phyllosphere bacterial communities in giant pandas’ staple food bamboos. Front. Microbiol. 12, 748141 (2021).
Hardoim, P. R., van Overbeek, L. S. & van Elsas, J. D. Properties of bacterial endophytes and their proposed role in plant growth. Trends Microbiol. 16, 463–471 (2008).
Li, Y. et al. Diazotroph paenibacillus triticisoli BJ-18 drives the variation in bacterial, diazotrophic and fungal communities in the rhizosphere and root/shoot endosphere of maize. Int. J. Mol. Sci. 22, 1460 (2021).
Pillay, V. K. & Nowak, J. Inoculum density, temperature, and genotype effects on in vitro growth promotion and epiphytic and endophytic colonization of tomato (Lycopersicon esculentum L) seedlings inoculated with a pseudomonad bacterium. Can. J. Microbiol. 43, 354–361 (1997).
Shade, A. & Stopnisek, N. Abundance-occupancy distributions to prioritize plant core microbiome membership. Curr. Opin. Microbiol. 49, 50–58 (2019).
Neu, A. T., Allen, E. E. & Roy, K. Defining and quantifying the core microbiome: Challenges and prospects. Proc. Natl Acad. Sci. USA 118, 2104429118 (2021).
Montañez, A., Abreu, C., Gill, P. R., Hardarson, G. & Sicardi, M. Biological nitrogen fixation in maize (Zea mays L.) by 15N isotope-dilution and identification of associated culturable diazotrophs. Biol. Fertil. Soils 45, 253–263 (2008).
Chelius, M. K. & Triplett, E. W. Immunolocalization of dinitrogenase reductase produced by Klebsiella pneumoniae in association with Zea mays L. Appl. Environ. Microbiol. 66, 783–787 (2000).
Santhanam, R. et al. A suite of complementary biocontrol traits allows a native consortium of root-associated bacteria to protect their host plant from a fungal sudden-wilt disease. Mol. Ecol. 28, 1154–1169 (2019).
Andrews, S. FastQC: a quality control tool for high throughput sequence data. http://www.bioinformatics.babraham.ac.uk/projects/fastqc (2010).
Edgar, R. C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461 (2010).
Bolyen, E. et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37, 852–857 (2019).
Rognes, T., Flouri, T., Nichols, B., Quince, C. & Mahe, F. VSEARCH: a versatile open-source tool for metagenomics. PeerJ 4, 22 (2016).
Edgar, R. C. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998 (2013).
Schloss, P. D. et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541 (2009).
Quast, C. et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596 (2013).
Zhang, J. et al. NRT1.1B is associated with root microbiota composition and nitrogen use in field-grown rice. Nat. Biotechnol. 37, 676–684 (2019).
Oksanen, J. et al. vegan: Community Ecology Package (2018).
Shenhav, L. et al. FEAST: fast expectation-maximization for microbial source tracking. Nat. Methods 16, 627–632 (2019).
Louca, S., Parfrey, L. W. & Doebeli, M. Decoupling function and taxonomy in the global ocean microbiome. Science 353, 1272–1277 (2016).
R Core Team R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/ (2019).
Wickham, H. ggplot2: Elegant Graphics For Data Analysis (Springer, 2016).
Robinson, M. D., McCarthy, D. J. & Smyth, G. K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140 (2009).
Conway, J. R., Lex, A. & Gehlenborg, N. UpSetR: an R package for the visualization of intersecting sets and their properties. Bioinformatics 33, 2938–2940 (2017).
Chen, T., Zhang, H., Liu, Y., Liu, Y. X. & Huang, L. EVenn: Easy to create repeatable and editable Venn diagrams and Venn networks online. J. Genet. Genomics 48, 863–866 (2021).
Wippel, K. et al. Host preference and invasiveness of commensal bacteria in the Lotus and Arabidopsis root microbiota. Nat. Microbiol. 6, 1150–1162 (2021).
Rastogi, G., Tech, J. J., Coaker, G. L. & Leveau, J. H. A PCR-based toolbox for the culture-independent quantification of total bacterial abundances in plant environments. J. Microbiol. Methods 83, 127–132 (2010).
Nakano, M. 16S rRNA gene primer validation for bacterial diversity analysis of vegetable products. J. Food Prot. 81, 848–859 (2018).
Poly, F., Monrozier, L. J. & Bally, R. Improvement in the RFLP procedure for studying the diversity of nifH genes in communities of nitrogen fixers in soil. Res. Microbiol. 152, 95–103 (2001).
Pfeilmeier, S. et al. The plant NADPH oxidase RBOHD is required for microbiota homeostasis in leaves. Nat. Microbiol. 6, 852–864 (2021).
Kim, M., Oh, H. S., Park, S. C. & Chun, J. Towards a taxonomic coherence between average nucleotide identity and 16S rRNA gene sequence similarity for species demarcation of prokaryotes. Int. J. Syst. Evol. Microbiol. 64, 346–351 (2014).
Letunic, I. & Bork, P. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 49, W293–W296 (2021).
Rosch, C., Mergel, A. & Bothe, H. Biodiversity of denitrifying and dinitrogen-fixing bacteria in an acid forest soil. Appl. Environ. Microbiol. 68, 3818–3829 (2002).
Cannon, F. C., Dixon, R. A. & Postgate, J. R. Derivation and properties of F-prime factors in Escherichia coli carrying nitrogen fixation genes from Klebsiella pneumoniae. J. Gen. Microbiol. 93, 111–125 (1976).
Yang, J., Xie, X., Wang, X., Dixon, R. & Wang, Y. P. Reconstruction and minimal gene requirements for the alternative iron-only nitrogenase in Escherichia coli. Proc. Natl Acad. Sci. USA 111, E3718–E3725 (2014).
Haskett, T. L., Knights, H. E., Jorrin, B., Mendes, M. D. & Poole, P. S. A simple in situ assay to assess plant-associative bacterial nitrogenase activity. Front. Microbiol. 12, 690439 (2021).
Chen, S. F., Zhou, Y. Q., Chen, Y. R. & Gu, J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, 884–890 (2018).
Bankevich, A. et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477 (2012).
Tatusova, T. et al. NCBI prokaryotic genome annotation pipeline. Nucleic Acids Res. 44, 6614–6624 (2016).
Huerta-Cepas, J. et al. Fast genome-wide functional annotation through orthology assignment by eggNOG-Mapper. Mol. Biol. Evol. 34, 2115–2122 (2017).
Kim, J., Kim, M. S., Koh, A. Y., Xie, Y. & Zhan, X. FMAP: functional mapping and analysis pipeline for metagenomics and metatranscriptomics studies. BMC Bioinf. 17, 420 (2016).
UniProt, C. UniProt: a hub for protein information. Nucleic Acids Res. 43, D204–D212 (2015).
Acknowledgements
We thank Peng Xu for helpful advice regarding assembly, annotation and interpretation of the genome of core strains, and Yongliang Yan, PhD, and Xiaoxia Zhang, PhD, for assistance in laboratory work. We also thank Jennifer Smith, PhD, from Liwen Bianji (Edanz) (www.liwenbianji.cn/), for editing the English text of a draft of this manuscript. This work was financially supported by the Young Elite Scientists Sponsorship Program by CAST (2017QNRC001).
Author information
Authors and Affiliations
Contributions
C.A. and L.P. conceived the study and supervised the project. L.Z., C.A., M.Z., S.H. and X.X. performed the experiments. L.L., Q.G., Y.W., S.Z., S.H., L.Y., Y.W., K.L., X.Y., D.L. and L.Z. coordinated field experiments. H.W. provided the GFP-containing plasmid. W.Z. and P.H reviewed the manuscript. L.Z, C.A. and L.P. wrote the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Source data
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Zhang, L., Zhang, M., Huang, S. et al. A highly conserved core bacterial microbiota with nitrogen-fixation capacity inhabits the xylem sap in maize plants. Nat Commun 13, 3361 (2022). https://doi.org/10.1038/s41467-022-31113-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-022-31113-w
This article is cited by
-
Core microbiota drive multi-functionality of the soil microbiome in the Cinnamomum camphora coppice planting
BMC Microbiology (2024)
-
Selenium and Bacillus proteolyticus SES increased Cu-Cd-Cr uptake by ryegrass: highlighting the significance of key taxa and soil enzyme activity
Environmental Science and Pollution Research (2024)
-
Plant biotechnology research with single-cell transcriptome: recent advancements and prospects
Plant Cell Reports (2024)
-
Enhancement of beneficial microbiomes in plant–soil continuums through organic fertilization: Insights into the composition and multifunctionality
Soil Ecology Letters (2024)
-
Influences of plant organ, genotype, and cultivation site on the endophytic bacteriome of maize (Zea mays L.) in the semi-arid region of Pernambuco, Brazil
Brazilian Journal of Microbiology (2024)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.