Pei Wang1,2†
Xiao Kong3,4†
Hongsong Chen5†
Youlun Xiao2
Huijun Liu6
Xiaojuan Li2
Zhuo Zhang1,2
Xinqiu Tan2
Diandong Wang7
Decai Jin4,8*
Ye Deng4,8
Tomislav Cernava9- 1College of Plant Protection, Hunan Agricultural University, Changsha, China
- 2Hunan Plant Protection Institute, Hunan Academy of Agricultural Sciences, Changsha, China
- 3School of Public Health, Qingdao University, Qingdao, China
- 4Key Laboratory of Environmental Biotechnology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing, China
- 5Guangxi Key Laboratory for Biology of Crop Diseases and Insect Pests, Institute of Plant Protection, Guangxi Academy of Agricultural Sciences, Nanning, China
- 6Beijing Key Laboratory of Detection and Control of Spoilage Organisms and Pesticide Residues in Agricultural Products, Beijing University of Agriculture, Beijing, China
- 7College of Life Science and Technology, Yangtze Normal University, Chongqing, China
- 8College of Resources and Environment, University of Chinese Academy of Sciences, Beijing, China
- 9Institute of Environmental Biotechnology, Graz University of Technology, Graz, Austria
Microbial communities associated with the plant phyllosphere and endosphere can have both beneficial as well as detrimental effects on their hosts. There is an ongoing debate to which extend the phyllosphere and endosphere microbiome assembly is controlled by the host plant how pronounced cultivar effects are. We investigated the bacterial and fungal communities from the phyllosphere and endosphere of 10 different rice cultivars grown under identical environmental conditions in the frame of a targeted approach to identify drivers of community assembly. The results indicated that the endophytic bacterial communities were clearly separated into two groups. The α-diversity and microbial network complexity within Group I were significantly lower than in Group II. Moreover, the genera Nocardioides, Microvirga, and Gaiella were significantly more abundant in Group II and only present in the interaction networks of this group. These three genera were significantly correlated with α- and β-diversity of the endophytic bacterial community and thus identified as major drivers of the endosphere community. We have identified keystone taxa that shape endophytic bacterial communities of different rice cultivars. Our overall findings provide new insights into plant-microbe interactions, and may contribute to targeted improvements of rice varieties in the future.
Introduction
With only a few rare exceptions, plants are generally associated with highly diversified microbial communities (Heil, 2015; Manzotti et al., 2020). These microorganisms are closely associated with ecosystem functions, which can be either beneficial or harmful (Graham et al., 2016). The phyllosphere of plants is likely the largest microbial habitat on earth (Redford et al., 2010) with a global leaf surface area estimated to be ∼6.4 × 108 km2 and contains up to 1026 bacteria (Lindow and Brandl, 2003; Vorholt, 2012). This habitat constitutes an oligotrophic environment and microorganisms that colonize it are affected by many factors, such as temperature, UV radiation, growing season, CO2 concentration, which all affect microbial communities to a certain extend (Remus-Emsermann et al., 2012; Copeland et al., 2015; Ren et al., 2015; Pan et al., 2019). In addition to microbes that colonize plant surfaces, many microorganisms, collectively referred to as endophytes, are also found within plant tissues.
Plant endophytes are associated with almost every plant on earth and fulfill various essential role therein (Strobel et al., 2004). These specific microorganism usually colonize root hairs, leaves, vascular tissues, epidermal cells, intercellular spaces, and the cytoplasm (Strobel et al., 2004). Endophytes use roots, flowers and natural openings on the leaves and stem lenticels as entry points to inner plant tissues (Lodewyckx et al., 2002; Kandel et al., 2017). In general, endophytes that grow within healthy plants often contribute with certain functions to their host’s growth and health for a mutual benefit. They can increase host disease resistance through the secretion of bioactive compounds (Busby et al., 2016; Matsumoto et al., 2021), augment nitrogen transfer from soil to plant (Behie and Bidochka, 2014), as well as improve crop performance (Liu et al., 2017).
Plant endophytes provide the basis for a broad range of applications in agriculture. They can be used to increase plant biomass, improve plant health, or degrade environmental contaminants (James et al., 2002; Sandhu et al., 2007; Puri et al., 2016; Feng et al., 2017a). In contrast, pathogenic microorganisms that can occur in the same niches in plants can cause various plant as well as human diseases (De Kempeneer et al., 2004; Teplitski et al., 2011). Therefore, a better understanding of indigenous plant-associated microbial communities will likely provide the means to counteract harmful aspects that are accompanied by pathogen invasion.
Rice (Oryza sativa L.) is the most commonly cultivated crop in China and also serves as a staple food in many other countries of the world. There are thousands of different rice varieties in the world, with new cultivars continuously arising from extensive breeding approaches. While phenotypic traits can substantially differ between cultivars, there is contradicting evidence to which extend the plant-associated microbiome is determined by the plant genotype. For example, Hunter et al. demonstrated that lettuce genotypes significantly affect phyllosphere bacterial and fungal communities (Hunter et al., 2010; Hunter et al., 2015). However, Lundberg et al. (2012) showed that the host genotypes had a weak effect on the root-associated microbiota of Arabidopsis thaliana. In terms of rice plants, there is currently no consensus on whether the plant genotype affects the associated microbial communities and to which extend they are affected. Hardoim et al. (2011) demonstrated that the formation of bacterial communities associated with rice roots was affected by the genotype. Different plant genotypes can have diverging compositions of their root exudates, which influences microorganisms in the rhizosphere. Interestingly, Wu et al. (2009) demonstrated that variations in the rice cultivar and nitrogen fertilizers only influenced methanogens that are tightly associated with rice roots. Due to the high variations in the plant phyllosphere, there is an ongoing debate on whether the microbiome assembly in this microenvironment is controlled by the host genotype. The intrinsic potential of microbial communities in the phyllosphere to modulate the community structure and composition was so far mostly ignored.
In this study, 10 rice cultivars were grown within the same plot and under the same management conditions to explore the relationship between plant cultivar and leaf-associated microbial communities. The bacterial and fungal communities from the phyllosphere and endosphere of these cultivars were analyzed by high-throughput sequencing of marker genes. The aims of the study were as follows: (1) To decipher the differences in the cultivars’ phyllosphere communities; (2) to investigate the potential of each cultivar to shape the phyllosphere microbiome; (3) to reveal potential keystone species that influence the rice leaf-associated microbiome.
Materials and Methods
Site Description and Sample Collection
All samples for this study were collected from an experimental field of the Rice Blast Identification Center in Taojiang, Hunan, China (112°06′34″E, 28°38′55″W). The 10 rice cultivars that were implemented in this study are described in Table 1. All cultivars were planted in one experimental field, separated by ditches, and received the same management program (including irrigation, fertilization, and pesticide treatements). The soil is a fertile, heavy loam in which rice can be grown during the whole year. The samples were collected at the tillering stage with no visible occurrence of pests or plant diseases.
Table 1. Rice cultivars implemented in the present study.
Five leaves from each plant and six replicates for each cultivar were randomly collected and placed in an isolated box with cooling packs for transport to the laboratory. The total community DNA from phyllosphere epiphytes and endophytes was collected as follows. The leaves were immersed in sterile PBS buffer (0.02 mM, pH 7.0, 0.1% Tween 80), and shaken on a thermostatic oscillator (BSD-YF2600, Shanghai Boxun Medical Biological Instrument Crop., Shanghai, China) for 30 min at 30°C, followed by 10 min of sonication (Ningbo Xingzhi Biological Technology Co., Ltd., Ningbo, China). To enrich epiphytic phyllosphere microorganisms the suspension was passed through a 0.22 μm membrane using a vacuum pump, the membrane was then stored at 4°C until DNA extraction.
Once phyllosphere microorganisms were collected the leaves were surface sterilized with 75% ethyl alcohol for 3 min, 2.5% sodium hypochlorite (NaClO) for 5 min, and washed with sterile water five times to remove residual reagents. The leaves were then ground with a sterile mortar and pestle in PBS, silica sand was used to increase the friction. The resulting supernatant was stored at 4°C until DNA extraction.
Total Community DNA Extraction and PCR Amplification of Marker Genes
Genomic DNA was extracted with the Fast DNA spin kit for soil (MP Biomedicals LLC, United States) according to the manufacturer’s instructions. The 16S rRNA gene was amplified with primer pair 799F (AACMGGATTAGATACCCKG)/1115R (AGGGTTGCGCTCGTTG) (Kong et al., 2018), while the ITS region of fungi was amplified with primer pair gITS7F (GTGARTCATCGARTCTTTG)/ITS4 (TCCTCCGCTTATTGATATGC) (Wei et al., 2018). Each primer pair contained a unique barcode (12 bp) to distinguish different samples after sequencing. The PCR reaction was as follows: 5 μL 10 × PCR buffer, 4 μL dNTP, 0.5 μL DNA polymerase (TaKaRa Biotech, Beijing, China), 1.5 μL each primer (10 μM) and 1 μL DNA template. The PCR thermocycling conditions were as follows: Initial denaturation at 94°C for 1 min, denaturation at 94°C for 20 s, annealing at 57°C for 25 s, and extension at 72°C for 45 s, the denaturation, annealing and extension processes were repeated 30 and 35 times for the amplification of 16S rDNA and ITS target fragment, respectively. A final extension at 72°C for 10 min was conducted for both reactions. The PCR products were visualized by electrophoresis in a 2% agarose gel and purified with E.Z.N.A. Gel Extraction Kit (Omega Bio-tek, Norcross, GA, United States) following the manufacturer’s instructions. Sequencing was conducted by the commercial provider Annoroad Gene Technology (Beijing, China), and the Illumina Miseq platform for 2 × 250 bp paired-end sequencing. A total of 240 samples was sequenced in this study (phyllosphere and endosphere samples, bacterial and fungal samples, ten cultivars, six replicates; 2 × 2 × 10 × 6 = 240).
Analysis of High-Throughput Sequence Data and Statistical Analysis
The raw data was subjected to initial processing with the online pipeline1 (Feng et al., 2017b). In the first step, the forward and reverse sequence files were merged into full length sequences by using FLASH (Magoc and Salzberg, 2011). Reads with average quality score <20 or containing ambiguous nucleotides (N) were removed from the analysis. Then, sequences were trimmed based on their length. All bacterial 16S rRNA gene fragment sequences with a length <295 or >305 bp were removed. For the fungal ITS region reads with <300 or >400 bp were removed. An OTU (operational taxonomic unit) table was generated at 97% similarity level with the UPARSE algorithm (Edgar, 2013). For the fungal analysis, the ITS sequences were verified with the ITSx tool in order to remove other eukaryotic sequences (Bengtsson-Palme et al., 2013). Taxonomic assignments were conducted with the RDP database and the Warcup database (Wang et al., 2007). Finally, to account for variations in sequencing depth, the OTU tables were resampled randomly and the resulting tables were used for downstream analysis. The sequencing data is publicly available at the NCBI Sequence Read Archive under accession no. SRP230471 (PRJNA586741).
Before deepening analyses were conducted, it was confirmed that the rarefaction curves were saturated which suggested that the sequencing depth was sufficient to reflect the microbial community (Supplementary Figure 1). The alpha diversity of microbial communities was assessed with the Shannon index. The community composition was analyzed at different phylogenetic levels (phylum and genus) and the beta diversity was visualized with a principal coordinate analysis (PCoA) based on a weighted UniFrac matrix. Dissimilarity tests based on the multi-response permutation procedure (MRPP), adonis (PERMANOVA), and analysis of similarities (ANOSIM) were used to evaluate differences among groups. Functional groups were identified with the FAPROTAX (Functional Annotation of Prokaryotic Taxa) (Louca et al., 2016) and FunGuild bioinformatic tools (Nguyen et al., 2016) for the bacterial and fungal communities, respectively. A Mantel test was used to calculate the correlation between environmental factors and the observed microbial community.
Construction of Microbial Interaction Networks
Interaction networks were constructed to assess the interaction between different taxa in the rice phyllosphere. For fungi, only an epiphyte network was constructed due to the limited number of replicates for fungal endophytes after the removal of samples with insufficient read numbers. The networks were constructed with an online pipeline available at http://ieg4.rccc.ou.edu/mena (Deng et al., 2012). First, the resampled OTU tables of each rice cultivar were used to construct individual networks. The filter was set to 80%, so that only operational taxonomic units (OTUs) present in a minimum of five of the six replicates were considered. Logarithmic data transformation was omitted and correlation calculations were performed with Spearman’s rho. The networks were constructed based on the RMP method with a threshold (cutoff) of 0.91, 0.95, and 0.86 for phyllosphere bacteria, endophytic bacteria, and phyllosphere fungi, respectively. Network properties, such as average connectivity, average path length, and average clustering coefficient (avgCC) were separately calculated for each dataset. Finally, the networks were visualized with Cytoscape 3.3.0 (Shannon et al., 2003).
Statistical Analyses
The statistical analyses were carried out with IBM SPSS statistics 21. Significance tests were performed with one-way ANOVA and the correlation analysis was carried out by assessing Pearson’s correlation coefficient.
Results
General Data Statistics and Microbial Community Diversity
Following initial quality filtering, several samples within the fungal endophyte dataset contained <10,000 reads after quality control, therefore only three replicates were available for the downstream analysis of most cultivars. In total 7,166,642 bacterial and 7,780,134 fungal high-quality sequences were retained after quality control. Moreover, ranges of 76 – 4,540 and 39 – 386 OTUs per sample were obtained after taxonomic classification for bacteria and fungi, respectively. Following resampling for data normalization of 16S rDNA gene fragment and ITS region samples, 10,635 and 11,711 reads per sample were retained, respectively.
The Shannon index indicated differences in microbial diversity of different samples (Figure 1). Among the endophytic bacterial community samples, cultivars E, F, G, H, and J (Group II) harbored more diverse communities than Group I comprising cultivars A, B, C, D, and K (P < 0.05). Among all samples, cultivar F was associated with the highest diversity (H’ = 6.08) while cultivar C (H’ = 1.44) had the lowest alpha diversity. For phyllosphere and endophytic fungal communities, no significant difference was found among the cultivars. Within the Group II cultivars, the order of diversity was: endophytic bacteria ≥ phyllosphere bacteria > endophytic fungi ≥ phyllosphere fungi, while for Group I cultivars, the diversity was: phyllosphere bacteria ≥ endophytic bacteria > endophytic fungi ≥ phyllosphere fungi.
Figure 1. Shannon indices of phyllosphere and endosphere microbial communities from ten different rice cultivars. BE, endophytic bacteria; BP, phyllosphere bacteria; FE, endophytic fungi; FP, phyllosphere fungi. The capital letters designate significant differences (P < 0.05) between different cultivars and the lowercase letters represent significant differences (P < 0.05) among different samples of the same cultivar (BE, BP, FE, and FP).
Microbial Community Composition of Rice
The bacterial community was primarily composed of the phyla Proteobacteria (53.01–97.21%), Actinobacteria (0.16–27.56%), and Firmicutes (0.68–13.19%) while the fungal community was dominated by Basidiomycota (77.11–100%) and Zygomycota (0–18.33%) (Supplementary Figure 2). The relative abundance of Actinobacteria and Firmicutes in the endosphere of Group II was higher than in Group I. Zygomycota were only present in endosphere of cultivars A, B, G, H, J, and K. The fungal community of phylloshere samples was almost entirely composed of the phylum Basidiomycota (99.88–100%).
In a deepening assessment the microbial community composition was assessed at genus level for all datasets. The endophytic bacterial community compositions of Group II cultivars were more similar to each other while the same trend was observed within Group I (Figure 2A). Within Group I, there were six bacterial genera in common for all cultivars, including Acinetobacter, Buttiauxella, and Serratia, while a total of 93 genera were in common for Group II, including the highly abundant Bacillus, Sphingomonas, Pantoea, and Nocardioides. Six genera, Acinetobacter, Buttiauxella, Cedecea, Pantoea, Serratia, and Sphingomonas, were present in all 10 rice cultivars. Significance analyses for differential abundance indicated that members of the genera Buttiauxella were significantly more abundant in the endosphere of Group I while Citrobacter, Arthrobacter, Nocardioides, Microvirga, Delftia, Gaiella, Planomiicrobium, Undibacterium, and Azohydromonas were significantly more abundant in the endosphere of Group II (P < 0.05). Overall, the species richness of Group II was higher than Group I, while the abundance was lower. In the bacterial phyllosphere community, a total of 14 genera, including Acinetobacter, Buttiauxella, Acidovorax, Sphingomonas, Cedecea, and Methylobacterium, were detected in all 10 rice cultivars. In addition, members of the genus Herbaspirillum were more abundant in cultivars A, H, J, and K. Buttiauxella was more abundant in cultivars B–D while Ralstonia was more abundant in cultivars E and F. The relative abundance of other bacteria showed no significant difference among the different cultivars.
Figure 2. Community composition within the bacterial (A) and fungal (B) community at genus level. For both communities the phyllosphere and the endosphere are separately shown. The fraction labeled with “other” represents the sum of all microorganisms with a relative abundance <2% in all community types.
The fungal genera Malassezia, Moesziomyces, and Tilletiaria were present in the endosphere and phyllosphere (epiphytes) of all ten rice cultivars. The endophytic fungal community was dominated by Pseudozyma (6.08–73.11%), Malassezia (9.76–70.03%) and Mortierella (0.11–18.20%) while the epiphytic phyllosphere fungal community was dominated by Pseudozyma (64.99–3.98%) and Tilletiopsis (11.46–30.12%) (Figure 2B). The relative abundance of Malassezia in the endosphere of cultivars F and K was lower than in the other cultivars. The relative abundance of Sterigmatomyces in the endosphere of cultivar C and Rhizophydium in the endosphere of cultivar A was higher than in the other cultivars. Within the phyllosphere samples, the relative abundance of Malassezia was significantly higher in cultivar E and Tilletiopsis was higher in cultivars D, G, and K when compared to other cultivars.
Microbial Community Structures Among Different Cultivars and Sample Types
A PCoA plot based on a weighted UniFrac matrix was used to visualize differences in community structures (Figure 3). The first axis (PC1) accounted for 39.37 and 78.57% difference while the second axis (PC2) accounted for 21.96 and 8.47% of the difference for bacterial and fungal communities, respectively. The phyllosphere and endophytic communities, both bacterial and fungal, were significantly different from each other. Also, the endophytic bacterial communities of the ten cultivars clearly clustered into two groups; Group I: cultivars A, B, C, D, and K and Group II: cultivars E, F, G, H, and J. The community structure of the phyllosphere bacteria, phyllosphere fungi, and endophytic fungi showed no significant difference among the cultivars.
Figure 3. Principal coordinate analysis (PCoA) plots based on weighted UniFrac matrix of bacterial (A) and fungal communities (B). The initial letters B and F stand for bacterial and fungal, respectively. The middle letters E and P stand for endosphere and phyllosphere, respectively. The last letter A–K designate the implemented cultivars.
Identification of Community Modulators via Microbial Interactions Networks
Different networks were constructed in order to infer intra-community interactions. Following initial assessments, only networks of the endophytic bacterial communities were further analyzed, because no significant differences were found in the networks from phyllosphere samples. The overall network topologies are shown in Table 2. The number of total nodes, total links, average degree (avgK), and number of modules in Group II were significantly higher than in Group I, indicating that cultivars E, F, G, H, and J hosted more complex networks. The degree of distribution in these networks fitted well with the power law (R2 of power law: 0.309 – 1), suggesting that these networks have scale-free properties, which indicates that most nodes have few neighbors while a few nodes have many neighbors. The avgCC and modularity in these empirical networks were significantly different from their corresponding random networks, which confirmed that these networks were not random. The average path lengths (GD) were within the range of 1.091 to 9.237, which were close to the logarithmic values of the total number of network nodes, an indication that these nodes were closely related to each other, and that the networks showed small-world properties. Modularity ranges from 0.504 to 0.857, indicated that these networks could be divided into modules and that a module was a functional unit.
Table 2. Topological properties of the empirical molecular ecological networks (MENs) and their associated random MENs.
Complementary Z-P plot analyses were implemented to further define the roles of individual nodes within the networks (Supplementary Figure 3). The network nodes were grouped according to four possible roles: peripherals, connectors, module hubs, and network hubs. In this study, the majority of nodes were assigned to peripherals (99.06%), while connectors and module hubs, were only found in cultivars E, F, G, and J, where they made up 0.39 and 0.55% of the nodes, respectively. No network hubs were found at all. The connectors and module hubs were assigned to Proteobacteria (Sphingomonas sp., Filomicrobium sp., Microvirga sp., and Shinella sp.), Actinobacteria (Blastococcus sp., Arthrobacter sp., Gaiella sp., and Bryobacter sp.), and Firmicutes (Cohnella sp.).
The ten networks mainly consisted of Proteobacteria (46.3 – 100%), which was also the only identified phylum present in cultivar D (Figure 4). Actinobacteria were the second most abundant phylum in the obtained networks (16.22 – 30.46%), except for cultivars C and D. The relative abundances of Proteobacteria and Actinobacteria were higher and lower, respectively, in Group II as compared to Group I. The ratio of positive interactions in Group II was also higher than in Group I. At the genus level, Nocardioides was the most abundant in the networks of Group II while it was not found in the networks of Group I. In addition, the genera Bacillus, Gaiella, Microvirga, Arthrobacter, Massilia, and Blastococcus were also predominant in the networks of Group II. On the contrary, the composition of Group I networks showed substantially less complexity with Herbaspirillum, Sphingomonas, Cedecea, Acinetobacter, and Methylobacterium as the dominant genera in different networks.
Figure 4. Bacterial endophyte interaction networks of ten rice cultivars. The different sections (A–J) represent the cultivars (A–H,J,K), respectively. The node colors represent different phyla and the node size represent the number of links with other nodes. Blue and red links indicate positive and negative interactions, respectively.
Identification of Keystone Taxa That Affect the Endophytic Bacterial Communities
Genera with relative abundance higher that 1% were implemented in correlation analyses with α-diversity of endophytic bacterial communities among different cultivars. The results indicated total 24 genera were significantly associated with α-diversity. These bacteria were further used to analyze the correlation with endophyte community structures. This resulted in ten genera, whereof three genera were also shown to be part of the microbial networks indicating intra-community interactions. These three genera were assigned to Nocardioides, Gaiella, and Microvirga, and their relative abundances in Group II (0.18–4.38, 0.85–2.98, and 0.86–2.82%) were substantially higher than that in Group I (0–0.03%, 0.01–0.13%, and 0–0.08%, respectively). The three genera showed a significantly positive correlation with α- and β-diversity (P < 0.05, Table 3). In a complementary assessment that included the plant phenotype, the leaf length was found to be negatively correlated with the relative abundance of these three genera (P < 0.05).
Table 3. Identified keystone bacterial genera that affect the endophytic bacterial community across different rice cultivars.
Discussion
Microbial communities within the plant phyllosphere and endosphere play an important role in promoting the growth and health of their hosts (Ren et al., 2015). In this study, analyses of microbial community compositions and network interactions within the phyllosphere and endosphere of 10 different rice cultivars were conducted.
While the phyllosphere epiphytic communities showed no clear trends, the endophytic bacterial communities of the analyzed rice cultivars were clearly separable into two groups. The diversity of the endophytic bacterial community in Group I (cultivars A, B, C, D, and K) was significantly lower than that of Group II (cultivars E, F, G, H, and J). This finding is similar to an observation by Jiang et al. (2017) in which the rhizosphere bacterial communities of blueberry cultivars clearly clustered into three groups. In general, higher diversity in a biological system correlates with higher resistance to environmental changes (Loreau and de Mazancourt, 2013), which suggests that the microbial ecosystem stability of Group I is lower than that of Group II. In addition, we observed that the diversity of the phyllsophere bacterial community in Group I was higher than the endophytic bacterial community, while in Group II the opposite trend was found. The observed differences between bacterial communities in the phyllosphere and the endosphere might be due to differing permeability of the cultivars for colonization by endophytes. However, no significant difference was found in the diversity of fungal communities located in the phyllosphere and endosphere among the ten cultivars.
The relative abundance of the bacterial phyla Actinobacteria and Firmicutes in the endosphere was significantly higher in Group II then compared to Group I. Although Group I contained fewer species, the genera Herbaspirillum, Buttiauxella, and Acinetobacter were more abundant therein. For Group II, a higher number of species was detected, but no dominant genus for the entire group was found. Herbaspirillum, is a Gram-negative bacterial genus that is commonly found in the endosphere of rice. On one hand, it was found to be a plant-growth promoting bacterium, while on the other hand it can also act as a mild pathogen that impairs rice growth by augmenting ethylene levels (Valdameri et al., 2017). In addition, nitrogen fixation is characteristic for this genus (Chubatsu et al., 2011). The genus Buttiauxella belongs to the Enterobacteriaceae family and was previously found in the endosphere of rice seeds (Wang et al., 2020). It was shown to improve plant growth and to enhance the remediation efficiency for cadmium (Cd) in the host plant (Wu et al., 2018). Various Acinetobacter species are known to exert plant growth-promoting effects and decrease adverse effects of PAHs in rice (Gao et al., 2006).
When fungal communities were analyzed, the endosphere of the ten cultivars was dominated by Malassezia and Pseudozyma while the phyllosphere fungal community was dominated by Pseudozyma and Tilletiopsis. Pseudozyma sp. is commonly found on leaves of rice and other plants (Kitamoto et al., 2011; Nasanit et al., 2015). It can provide protection against plant-pathogenic fungi by producing distinct extracellular fatty acids (Avis and Bélanger, 2002). Tilletiopsis spp. were previously found in intercellular fluid from leaf blades and sheaths of rice, and they were shown to be promising biocontrol agents for the control of powdery mildew (Haggag et al., 2015). In the present study, the relative abundance of potentially beneficial microorganisms differed among rice cultivars, indicating that their ecological implications might be redundant, i.e., replaceable by other community members.
In the last decade, various studies reported that plant cultivars shape their associated microbiomes. For example, Manter et al. (2010) demonstrated that the endophyte bacterial community in potato roots is cultivar-dependent. Analogous to this finding, the study of Lopez-Velasco et al. (2011) indicated that the phyllosphere bacterial community of spinach is also determined by the cultivar. Furthermore, Materatski et al. (2019) demonstrated that the plant cultivar is one of the main factors that shapes the endophytic fungal community. However, there is still some controversy regarding the relative importance of the environment in terms of shaping the microbiome when compared to cultivar effects (Morella et al., 2020). For example, Kong et al. (2020) indicated that bacterial and fungal communities located in the phyllosphere of three different maize cultivars showed no significant difference. Singh et al. (2019) demonstrated that the phyllosphere microbial community of grapes was affected by the grape species and the growing year, with the latter having a more pronounced effect. In the present study, in order to exclude the impact of environmental factors, all 10 rice cultivars were grown under identical environmental conditions in the same field. Our results indicated that the differences among phyllosphere bacterial, phyllosphere fungal, and endophytic fungal communities were not statistically significant. However, the endophytic bacterial community was separable into two groups, with cultivars A, B, C, D, and K (Group I) and cultivars E, F, G, H, and J (Group II) clustering closer together and showing the same trends in their α-diversity. The genotypic relatedness of the different cultivars could not explain the group formation. For example, cultivar H and cultivar K, cultivar A and cultivar E, were derived from the same female parent, but their endophytic bacterial communities were significantly different. This is in contrast to cultivars of other plants where endophytic communities were shown to correspond to the genotype (Abdullaeva et al., 2020; Chen et al., 2020). Our observation therefore indicated that the analyzed rice cultivars were not the main drivers of the structure of leaf-associated microbial communities. In order to explore intrinsic factors within microbial communities, interaction networks of endophytic bacterial communities were screened for keystone taxa.
Networks analyses are suitable to infer interactions among different microbial species, such as competition and mutualism (Zhou et al., 2011; Deng et al., 2012). In the present study, the interaction networks of bacterial endosphere communities showed that the network complexity of cultivar Group I was significantly lower than Group II, indicating that the microbial inter-connectivity within Group II was higher than within Group I. In addition, Nocardioides, Microvirga, and Gaiella were found as the most important members within interaction networks of cultivar Group II while they were not present in networks of Group I. Nocardioides sp. is a common endophyte and has been isolated from the endosphere of many plants, such as Jatropha curcas L. (Qin et al., 2012), maize (Glaeser et al., 2014; Kampfer et al., 2016), mugwort and horse-weed (Han et al., 2013), Carex scabrifolia Steud (Song et al., 2011), and Perilla frutescens (Du et al., 2013). Nocardioides sp. contains dozens of species and some of them exhibit plant growth-promoting activities through nitrogen fixation (Nafis et al., 2019). Other species, such as N. alkalitolerans, N. jensenii, N. rotundus, and N. mesophilus are capable of reducing nitrate (Toth et al., 2008; Dastager et al., 2010; Wang et al., 2016). Representatives of Microvirga were also isolated from some plant tissues, such as the root endosphere of rapeseed (Jimenez-Gomez et al., 2019), roots endosphere of Abulitonindicum (Bhadania et al., 2019), as well as the rice rhizosphere (Li et al., 2020). Moreover, Microvirga spp. were also found to be nitrogen-fixing and plant growth-promoting bacteria (Jimenez-Gomez et al., 2019; Zhao et al., 2019). Gaiella spp. were previously detected in the endosphere of Lolium perenne L. and found to play an important role in the plant’s interaction network (Wemheuer et al., 2017). Hu et al. (2020) detected Gaiella spp. in bulk soil and rhizosphere soil of tobacco. In addition, Gaiella can reduce nitrate to nitrite (Albuquerque et al., 2011) which is essential for nitrogen cycling. These three genera were significantly correlated with both α- and β-diversity of the endophytic bacterial communities indicating a key role in community modulation. While direct cultivar effects on microbial communities can be low in distinct plants, we hypothesize that a low number of host-selected microbes (e.g., Nocardioides, Microvirga, and Gaiella found in in Group II of the present study) can exert substantial intra-community modulating effects. Host-microbe interplay was further confirmed in the present study by correlation analyses between the relative abundance of the three identified community-shaping genera and their host plant’s leaf length. The results indicated distinct phenotypic adaptions of the host in response to these taxa. Further research will be required to decipher the underlying mechanisms and to identify potential effects on the overall host performance as a response to these microbes.
Conclusion
Our study provides important clues related to the assembly of microbial communities in rice. It was found that 10 rice cultivars can be either assigned to a group with microbial communities of relatively low complexity (Group I) or to a group with relatively high complexity (Group II). Correlation analyses revealed that the genera Nocardioides, Microvirga, and Gaiella significantly correlated with the α-diversity as well as the β-diversity of the endophytic bacterial communities. The relative abundance of the three genera was significantly higher in Group II where they were found to be embedded in complex intra-community interaction networks, which is indicative for an important role in community modulation. In summary, our results indicate that the plant cultivar can enrich distinct endophytes in certain cases, which then exert further modulating effects on the remaining microbial community. To the best of our knowledge, this is the first study to comprehensively investigate plant microbial communities among different rice cultivars and mechanistically uncover potential causes for community variance.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/genbank/, SRP230471.
Author Contributions
All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.
Funding
This study received financial support by the National Key Research and Development Program of China (No. 2017YFD0800102), the National Natural Science Foundation of China (Nos. 31471831 and 41977122), and the opening funding of Guangxi Key Laboratory of Biology for Crop Diseases and Insect Pests (2019-KF-03).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.629852/full#supplementary-material
Footnotes
References
Abdullaeva, Y., Ambika Manirajan, B., Honermeier, B., Schnell, S., and Cardinale, M. (2020). Domestication affects the composition, diversity, and co-occurrence of the cereal seed microbiota. J. Adv. Res. 2020:8. doi: 10.1016/j.jare.2020.12.008
Albuquerque, L., França, L., Rainey, F. A., Schumann, P., Nobre, M. F., and da Costa, M. S. (2011). Gaiella occulta gen. nov., sp. nov., a novel representative of a deep branching phylogenetic lineage within the class Actinobacteria and proposal of Gaiellaceae fam. nov. and Gaiellales ord. nov. Syst. Appl. Microbiol. 34, 595–599. doi: 10.1016/j.syapm.2011.07.001
Avis, T. J., and Bélanger, R. R. (2002). Mechanisms and means of detection of biocontrol activity of Pseudozyma yeasts against plant-pathogenic fungi. FEMS Yeast Res. 2, 5–8. doi: 10.1016/S1567-1356(01)00058-7
Behie, S. W., and Bidochka, M. J. (2014). Nutrient transfer in plant-fungal symbioses. Trends Plant Sci. 19, 734–740. doi: 10.1016/j.tplants.2014.06.007
Bengtsson-Palme, J., Ryberg, M., Hartmann, M., Branco, S., Wang, Z., Godhe, A., et al. (2013). Improved software detection and extraction of ITS1 and ITS2 from ribosomal ITS sequences of fungi and other eukaryotes for analysis of environmental sequencing data. Methods Ecol. Evol. 4, 914–919. doi: 10.1111/2041-210x.12073
Bhadania, R. A., Parkhiya, M. V., Akbari, D. L., Akbari, L. F., Bhalani, H. N., and Golakiya, B. A. (2019). Microscopic colonization of salt tolerant plant–bacterium Microvirga zambiensis strain WSM 3693 in cotton tissue using auto fluorescent proteins GFP. Int. J. Curr. Microbiol. App. Sci. 8, 3449–3453.
Busby, P. E., Ridout, M., and Newcombe, G. (2016). Fungal endophytes: modifiers of plant disease. Plant Mol. Biol. 90, 645–655. doi: 10.1007/s11103-015-0412-0
Chen, X. Y. L., Krug, L., Yang, H., Li, H. X., Yang, M. F., Berg, G., et al. (2020). Nicotiana tabacum seed endophytic communities share a common core structure and genotype-specific signatures in diverging cultivars. Comput. Struct. Biotechnol. J. 18, 287–295. doi: 10.1016/j.csbj.2020.01.004
Chubatsu, L. S., Monteiro, R. A., de Souza, E. M., de Oliveira, M. A. S., Yates, M. G., Wassem, R., et al. (2011). Nitrogen fixation control in Herbaspirillum seropedicae. Plant Soil 356, 197–207. doi: 10.1007/s11104-011-0819-6
Copeland, J. K., Yuan, L., Layeghifard, M., Wang, P. W., and Guttman, D. S. (2015). Seasonal community succession of the phyllosphere microbiome. Mol. Plant Microbe Interact. 28, 274–285. doi: 10.1094/MPMI-10-14-0331-FI
Dastager, S. G., Lee, J. C., Pandey, A., and Kim, C. J. (2010). Nocardioides mesophilus sp. nov., isolated from soil. Int. J. Syst. Evol. Microbiol. 60(Pt 10), 2288–2292. doi: 10.1099/ijs.0.019059-0
De Kempeneer, L., Sercu, B., Vanbrabant, W., Van Langenhove, H., and Verstraete, W. (2004). Bioaugmentation of the phyllosphere for the removal of toluene from indoor air. Appl. Microbiol. Biotechnol. 64, 284–288. doi: 10.1007/s00253-003-1415-3
Deng, Y., Jiang, Y. H., Yang, Y. F., He, Z. L., Luo, F., and Zhou, J. Z. (2012). Molecular ecological network analyses. BMC Bioinformatics 13:113. doi: 10.1186/1471-2105-13-113
Du, H. J., Wei, Y. Z., Su, J., Liu, H. Y., Ma, B. P., Guo, B. L., et al. (2013). Nocardioides perillae sp. nov., isolated from surface-sterilized roots of Perilla frutescens. Int. J. Syst. Evol. Microbiol. 63(Pt 3), 1068–1072. doi: 10.1099/ijs.0.044982-0
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604
Feng, F., Ge, J., Li, Y., He, S., Zhong, J., Liu, X., et al. (2017a). Enhanced degradation of chlorpyrifos in rice (Oryza sativa L.) by five strains of endophytic bacteria and their plant growth promotional ability. Chemosphere 184, 505–513. doi: 10.1016/j.chemosphere.2017.05.178
Feng, K., Zhang, Z. J., Cai, W. W., Liu, W. Z., Xu, M. Y., Yin, H. Q., et al. (2017b). Biodiversity and species competition regulate the resilience of microbial biofilm community. Mol. Ecol. 26, 6170–6182. doi: 10.1111/mec.14356
Gao, Y., Yu, X. Z., Wu, S. C., Cheung, K. C., Tam, N. F. Y., Qian, P. Y., et al. (2006). Interactions of rice (Oryza sativa L.) and PAH-degrading bacteria (Acinetobacter sp.) on enhanced dissipation of spiked phenanthrene and pyrene in waterlogged soil. Sci. Total Environ. 372, 1–11. doi: 10.1016/j.scitotenv.2006.09.029
Glaeser, S. P., McInroy, J. A., Busse, H. J., and Kampfer, P. (2014). Nocardioides zeae sp. nov., isolated from the stem of Zea mays. Int. J. Syst. Evol. Microbiol. 64(Pt 7), 2491–2496. doi: 10.1099/ijs.0.061481-0
Graham, E. B., Knelman, J. E., Schindlbacher, A, Siciliano, S., Breulmann, M., Yannarell, A., et al. (2016). Microbes as engines of ecosystem function: When does community structure enhance predictions of ecosystem processes? Front. Microbiol. 7:214. doi: 10.3389/fmicb.2016.00214
Haggag, W. M., Abouziena, H. F., Hoballah, E., and Saber, M. (2015). Tilletiopsis pallescens: Optimization, fermentation and manufacture for managing mango powdery mildew. J. Chemical Pharmaceut. Res. 7, 341–352.
Han, J. H., Kim, T. S., Joung, Y., Kim, M. N., Shin, K. S., Bae, T., et al. (2013). Nocardioides endophyticus sp. nov. and Nocardioides conyzicola sp. nov., isolated from herbaceous plant roots. Int. J. Syst. Evol. Microbiol. 63(Pt 12), 4730–4734. doi: 10.1099/ijs.0.054619-0
Hardoim, P. R., Andreote, F. D., Reinhold-Hurek, B., Sessitsch, A., van Overbeek, L. S., and van Elsas, J. D. (2011). Rice root-associated bacteria: insights into community structures across 10 cultivars. FEMS Microbiol. Ecol. 77, 154–164. doi: 10.1111/j.1574-6941.2011.01092.x
Heil, M. (2015). “The Microbe-Free Plant: Fact or Artifact,” in Biological Nitrogen Fixation, ed. F. J. D. Bruijn (Hoboken, NJ: John Wiley & Sons, Inc), 1163–1174. doi: 10.1002/9781119053095.ch114
Hu, Q. L., Tan, L., Gu, S. S., Xiao, Y., Xiong, X. Y., Zeng, W. A., et al. (2020). Network analysis infers the wilt pathogen invasion associated with non-detrimental bacteria. NPJ Biofilms Microbi. 6:8. doi: 10.1038/s41522-020-0117-2
Hunter, P. J., Hand, P., Pink, D., Whipps, J. M., and Bending, G. D. (2010). Both leaf properties and microbe-microbe interactions influence within-species variation in bacterial population diversity and structure in the lettuce (Lactuca Species) phyllosphere. Appl. Environ. Microbiol. 76, 8117–8125. doi: 10.1128/AEM.01321-10
Hunter, P. J., Pink, D. A. C., and Bending, G. D. (2015). Cultivar-level genotype differences influence diversity and composition of lettuce (Lactuca sp.) phyllosphere fungal communities. Fungal Ecol. 17, 183–186. doi: 10.1016/j.funeco.2015.05.007
James, E. K., Gyaneshwar, P., Mathan, N., Barraquio, W. L., Reddy, P. M., Iannetta, P. P., et al. (2002). Infection and colonization of rice seedlings by the plant growth-promoting bacterium Herbaspirillum seropedicae Z67. Mol. Plant Microbe Interact. 15, 894–906. doi: 10.1094/MPMI.2002.15.9.894
Jiang, Y. J., Li, S. Z., Li, R. P., Zhang, J., Liu, Y. H., Lv, L. F., et al. (2017). Plant cultivars imprint the rhizosphere bacterial community composition and association networks. Soil Biol. Biochem. 109, 145–155. doi: 10.1016/j.soilbio.2017.02.010
Jimenez-Gomez, A., Saati-Santamaria, Z., Igual, J. M., Rivas, R., Mateos, P. F., and Garcia-Fraile, P. (2019). Genome Insights into the Novel Species Microvirga brassicacearum, a Rapeseed Endophyte with Biotechnological Potential. Microorganisms 7:7090354. doi: 10.3390/microorganisms7090354
Kampfer, P., Glaeser, S. P., McInroy, J. A., and Busse, H. J. (2016). Nocardioides zeicaulis sp. nov., an endophyte actinobacterium of maize. Int. J. Syst. Evol. Microbiol. 66, 1869–1874. doi: 10.1099/ijsem.0.000959
Kandel, S. L., Joubert, P. M., and Doty, S. L. (2017). Bacterial endophyte colonization and distribution within plants. Microorganisms 5:77. doi: 10.3390/microorganisms5040077
Kitamoto, H. K., Shinozaki, Y., Cao, X. H., Morita, T., Konishi, M., Tago, K., et al. (2011). Phyllosphere yeasts rapidly break down biodegradable plastics. AMB Express 1:44. doi: 10.1186/2191-0855-1-44
Kong, X., Han, Z. F., Tai, X., Jin, D. C., Ai, S., Zheng, X. X., et al. (2020). Maize (Zea mays L. Sp.) varieties significantly influence bacterial and fungal community in bulk soil, rhizosphere soil and phyllosphere. FEMS Microbiol. Ecol. 96:fiaa020. doi: 10.1093/femsec/fiaa020
Kong, X., Jin, D. C., Jin, S. L., Wang, Z. G., Yin, H. Q., Xu, M. Y., et al. (2018). Responses of bacterial community to dibutyl phthalate pollution in a soil-vegetable ecosystem. J. Hazard. Mater. 353, 142–150. doi: 10.1016/j.jhazmat.2018.04.015
Li, J., Gao, R. Y., Chen, Y., Xue, D., Han, J. H., Wang, J., et al. (2020). Isolation and Identification of Microvirga thermotolerans HR1, a Novel Thermo-Tolerant Bacterium, and Comparative Genomics among Microvirga Species. Microorganisms 8:8010101. doi: 10.3390/microorganisms8010101
Lindow, S. E., and Brandl, M. T. (2003). Microbiology of the Phyllosphere. Appl. Environ. Microbiol. 69, 1875–1875. doi: 10.1128/aem.69.4.1875-1883.2003
Liu, H., Carvalhais, L. C., Crawford, M., Singh, E., Dennis, P. G., Pieterse, C. M. J., et al. (2017). Inner Plant Values: Diversity, Colonization and Benefits from Endophytic Bacteria. Front. Microbiol. 8:2552. doi: 10.3389/fmicb.2017.02552
Lodewyckx, C., Vangronsveld, J., Porteous, F., Moore, E. R. B., Taghavi, S., Mezgeay, M., et al. (2002). Endophytic Bacteria and Their Potential Applications. Crit. Rev. Plant Sci. 21, 583–606. doi: 10.1080/0735-260291044377
Lopez-Velasco, G., Welbaum, G. E., Falkinham, J. O. III, and Ponder, M. A. (2011). Phyllopshere bacterial community structure of Spinach (Spinacia oleracea) as affected by cultivar and environmental conditions at time of harvest. Diversity 3, 721–738. doi: 10.3390/d30x000x
Loreau, M., and de Mazancourt, C. (2013). Biodiversity and ecosystem stability: a synthesis of underlying mechanisms. Ecol. Lett. 16, 106–115. doi: 10.1111/ele.12073
Louca, S., Parfrey, L. W., and Doebeli, M. (2016). Decoupling function and taxonomy in the global ocean microbiome. Science 353:1272. doi: 10.1126/science.aaf4507
Lundberg, D. S., Lebeis, S. L., Paredes, S. H., Yourstone, S., Gehring, J., Malfatti, S., et al. (2012). Defining the core Arabidopsis thaliana root microbiome. Nature 488, 86–90. doi: 10.1038/nature11237
Magoc, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507
Manter, D. K., Delgado, J. A., Holm, D. G., and Stong, R. A. (2010). Pyrosequencing reveals a highly diverse and cultivar-specific bacterial endophyte community in potato roots. Microb. Ecol. 60, 157–166. doi: 10.1007/s00248-010-9658-x
Manzotti, A., Bergna, A., Burow, M., Jorgensen, H. J. L., Cernava, T., Berg, G., et al. (2020). Insights into the community structure and lifestyle of the fungal root endophytes of tomato by combining amplicon sequencing and isolation approaches with phytohormone profiling. FEMS Microbiol. Ecol. 96:fiaa052. doi: 10.1093/femsec/fiaa052
Materatski, P., Varanda, C., Carvalho, T., Dias, A. B., Campos, M. D., Rei, F., et al. (2019). Spatial and temporal variation of fungal endophytic richness and diversity associated to the phyllosphere of olive cultivars. Fungal Biol. 123, 66–76. doi: 10.1016/j.funbio.2018.11.004
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
Morella, N. M., Weng, F. C., Joubert, P. M., Metcalf, C. J. E., Lindow, S., and Koskella, B. (2020). Successive passaging of a plant-associated microbiome reveals robust habitat and host genotype-dependent selection. Proc. Natl. Acad. Sci. U S A 117, 1148–1159. doi: 10.1073/pnas.1908600116
Nafis, A., Raklami, A., Bechtaoui, N., El Khalloufi, F., El Alaoui, A., Glick, B. R., et al. (2019). Actinobacteria from Extreme Niches in Morocco and Their Plant Growth-Promoting Potentials. Diversity 11:139. doi: 10.3390/d11080139
Nasanit, R., Krataithong, K., Tantirungkij, M., and Limtong, S. (2015). Assessment of epiphytic yeast diversity in rice (Oryza sativa) phyllosphere in Thailand by a culture-independent approach. Antonie Van Leeuwenhoek 107, 1475–1490. doi: 10.1007/s10482-015-0442-2
Nguyen, N. H., Song, Z., Bates, S. T., Branco, S., Tedersoo, L., Menke, J., et al. (2016). FUNGuild: An open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 20, 241–248. doi: 10.1016/j.funeco.2015.06.006
Pan, J. G., Jin, D. C., Jiang, H. M., Leng, X. Y., Zhang, A. A., and Bai, Z. H. (2019). The phyllosphere indigenous microbiota of Brassica campestris L. change its diversity in responding to di-n-butyl phthalate pollution. World J. Microbiol. Biotechnol. 35:21. doi: 10.1007/s11274-019-2589-x
Puri, A., Padda, K. P., and Chanway, C. P. (2016). Seedling growth promotion and nitrogen fixation by a bacterial endophyte Paenibacillus polymyxa P2b-2R and its GFP derivative in corn in a long-term trial. Symbiosis 69, 123–129. doi: 10.1007/s13199-016-0385-z
Qin, S., Yuan, B., Zhang, Y. J., Bian, G. K., Tamura, T., Sun, B. Z., et al. (2012). Nocardioides panzhihuaensis sp. nov., a novel endophytic actinomycete isolated from medicinal plant Jatropha curcas L. Antonie van Leeuwenhoek 102, 353–360. doi: 10.1007/s10482-012-9745-8
Redford, A. J., Bowers, R. M., Knight, R., Linhart, Y., and Fierer, N. (2010). The ecology of the phyllosphere: geographic and phylogenetic variability in the distribution of bacteria on tree leaves. Environ. Microbiol. 12, 2885–2893. doi: 10.1111/j.1462-2920.2010.02258.x
Remus-Emsermann, M. N., Tecon, R., Kowalchuk, G. A., and Leveau, J. H. (2012). Variation in local carrying capacity and the individual fate of bacterial colonizers in the phyllosphere. ISME J. 6, 756–765. doi: 10.1038/ismej.2011.209
Ren, G., Zhu, C., Alam, M. S., Tokida, T., Sakai, H., Nakamura, H., et al. (2015). Response of soil, leaf endosphere and phyllosphere bacterial communities to elevated CO2 and soil temperature in a rice paddy. Plant Soil 392, 27–44. doi: 10.1007/s11104-015-2503-8
Sandhu, A., Halverson, L. J., and Beattie, G. A. (2007). Bacterial degradation of airborne phenol in the phyllosphere. Environ. Microbiol. 9, 383–392. doi: 10.1111/j.1462-2920.2006.01149.x
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303
Singh, P., Santoni, S., Weber, A., This, P., and Peros, J. P. (2019). Understanding the phyllosphere microbiome assemblage in grape species (Vitaceae) with amplicon sequence data structures. Sci. Rep. 9:14294. doi: 10.1038/s41598-019-50839-0
Song, G. C., Yasir, M., Bibi, F., Chung, E. J., Jeon, C. O., and Chung, Y. R. (2011). Nocardioides caricicola sp. nov., an endophytic bacterium isolated from a halophyte, Carex scabrifolia Steud. Int. J. Syst. Evol. Microbiol. 61(Pt 1), 105–109. doi: 10.1099/ijs.0.019919-0
Strobel, G., Daisy, B., Castillo, U., and Harper, J. (2004). Natural Products from Endophytic Microorganisms. J. Nat. Prod. 67, 257–268. doi: 10.1021/np030397v
Teplitski, M., Warriner, K., Bartz, J., and Schneider, K. R. (2011). Untangling metabolic and communication networks: interactions of enterics with phytobacteria and their implications in produce safety. Trends Microbiol. 19, 121–127. doi: 10.1016/j.tim.2010.11.007
Toth, E. M., Keki, Z., Homonnay, Z. G., Borsodi, A. K., Marialigeti, K., and Schumann, P. (2008). Nocardioides daphniae sp. nov., isolated from Daphnia cucullata (Crustacea: Cladocera). Int. J. Syst. Evol. Microbiol. 58, 78–83. doi: 10.1099/ijs.0.65305-0
Valdameri, G., Alberton, D., Moure, V. R., Kokot, T. B., Kukolj, C., Brusamarello-Santos, L. C. C., et al. (2017). Herbaspirillum rubrisubalbicans, a mild pathogen impairs growth of rice by augmenting ethylene levels. Plant Mol. Biol. 94, 625–640. doi: 10.1007/s11103-017-0629-1
Vorholt, J. A. (2012). Microbial life in the phyllosphere. Nat. Rev. Microbiol. 10, 828–840. doi: 10.1038/nrmicro2910
Wang, L. N., Li, J. L., and Zhang, G. Y. (2016). Nocardioides rotundus sp. nov., isolated from deep seawater. Int. J. Syst. Evol. Microbiol. 66, 1932–1936. doi: 10.1099/ijsem.0.000966
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/aem.00062-07
Wang, Z. S., Zhu, Y. Q., Jing, R. X., Wu, X. Y., Li, N., Liu, H., et al. (2020). High-throughput sequencing-based analysis of the composition and diversity of endophytic bacterial community in seeds of upland rice. Arch. Microbiol. 2020, 1–12. doi: 10.21203/rs.3.rs-23885/v1
Wei, Z. Y., Liu, Y. Y., Feng, K., Li, S. Z., Wang, S., Jin, D. C., et al. (2018). The divergence between fungal and bacterial communities in seasonal and spatial variations of wastewater treatment plants. Sci. Total Environ. 628-629, 969–978. doi: 10.1016/j.scitotenv.2018.02.003
Wemheuer, F., Kaiser, K., Karlovsky, P., Daniel, R., Vidal, S., and Wemheuer, B. (2017). Bacterial endophyte communities of three agricultural important grass species differ in their response towards management regimes. Sci. Rep. 7:40914. doi: 10.1038/srep40914
Wu, K., Luo, J., Li, J., An, Q., Yang, X., Liang, Y., et al. (2018). Endophytic bacterium Buttiauxella sp. SaSR13 improves plant growth and cadmium accumulation of hyperaccumulator Sedum alfredii. Environ. Sci. Pollut. Res. 25, 21844–21854. doi: 10.1007/s11356-018-2322-6
Wu, L. Q., Ma, K., Li, Q., Ke, X. B., and Lu, Y. H. (2009). Composition of archaeal community in a paddy field as affected by rice cultivar and N fertilizer. Microb. Ecol. 58, 819–826. doi: 10.1007/s00248-009-9554-4
Zhao, Y. X., Jiang, H. Y., Cheng, X., Zhu, Y. X., Fan, Z. X., Dai, Z. L., et al. (2019). Neonicotinoid thiacloprid transformation by the N2-fixing bacterium Microvirga flocculans CGMCC 1.16731 and toxicity of the amide metabolite. Int. Biodeterior. Biodegrad. 145:104806. doi: 10.1016/j.ibiod.2019.104806
Keywords: Oryza sativa, phyllosphere, endosphere, microbial community, plant-microbe interactions
Citation: Wang P, Kong X, Chen H, Xiao Y, Liu H, Li X, Zhang Z, Tan X, Wang D, Jin D, Deng Y and Cernava T (2021) Exploration of Intrinsic Microbial Community Modulators in the Rice Endosphere Indicates a Key Role of Distinct Bacterial Taxa Across Different Cultivars. Front. Microbiol. 12:629852. doi: 10.3389/fmicb.2021.629852
Received: 16 November 2020; Accepted: 20 January 2021;
Published: 16 February 2021.
Edited by:
Prem Lal Kashyap, Indian Institute of Wheat and Barley Research (ICAR), IndiaReviewed by:
Clarisse Brígido, University of Evora, PortugalKulandaivelu Velmourougane, Central Institute for Cotton Research (ICAR), India
Copyright © 2021 Wang, Kong, Chen, Xiao, Liu, Li, Zhang, Tan, Wang, Jin, Deng and Cernava. 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: Decai Jin, dcjin@rcees.ac.cn
†These authors have contributed equally to this work