Metagenomic analyses of bacterial endophytes associated with the phyllosphere of a Bt maize cultivar and its isogenic parental line from South Africa

Structure and Diversity of Endophytic Bacteria in Maize Seeds and Germinating Roots

Seed endophytes in maize, which facilitate the transmission of microorganisms from one plant generation to the next, may play a crucial role in plant protection and growth promotion. This study aimed to investigate the effects of various maize varieties on the communities of endophytic bacteria in seeds and germinating roots. This study utilized Illumina high-throughput sequencing technology to examine the structural and diversity differences of endophytic bacterial communities within seed maize (BY1507), silage maize (QQ446), and wild maize (Teosinte) in both seeds and germinating roots. The results showed that 416 bacterial genera were detected, with Pantoea, Lachnospiraceae, Pararhizobium, Enterobacteriaceae, Stenotrophomonas, and Pseudonocardia being the most prevalent (relative abundance > 10%) at the genus level. No significant difference was observed in diversity indices (Chao1, ACE, Shannon, and Simpson) of seed endophytes among BY1507, QQ446, and Teosinte. The Shannon and Simpson indices for the germinating root endophyte from the wild variety (Teosinte) were significantly higher than the domesticated varieties (BY1507 and QQ446). PCoA revealed a notable overlap in the endophytic bacterial communities from the seeds of BY1507, QQ446, and Teosinte. Yet, clustering patterns were found. Co-occurrence network analysis showed that BY1507, QQ446, and Teosinte share a notable proportion of shared endophytic bacteria (>30%) between the seeds and germinating roots. This investigation elucidates the characteristics of endophytic microbial communities of seeds and germinating roots with seed maize, silage maize, and wild maize, offering data for future research on the physiological ecological adaptation of these endophytic microbial communities.

Seasonal Variation Imparts the Shift in Endophytic Bacterial Community Between Mango and its Hemiparasites

The study of community composition and community structure is important to know the ecological behaviour and community dynamics of the participating species and to understand the molecular interplay that lies between them. The community diversity greatly lies in the physiological status of the host and the environmental factors. The present study aims to explore the endophytic bacterial communities and their dynamics in the pre-flowering and post-flowering seasons in the horticulturally important Mango (Mangifera indica L.) and its hemiparasites: Loranthus parasiticus (L.) Marr. and Macrosolen colchinchinensis (Lour.) Tiegh. through a metagenomic approach using the sequence of V3 region of 16S rRNA gene. The genera Bacillus, Acinetobacter and Corynebacterium, under the phyla Firmicutes, Proteobacteria and Actinobacteria, respectively, were found to be the most abundant genera present in mango and its hemiparasites. It was found that during the post-flowering season, the twigs and leaves of mango had lesser endophytes than in other seasons while the alpha-diversity indices of the representative genera were the highest in L. parasiticus during the same seasons. However, in M. colchinchinensis, the alpha diversity was also higher in the post-flowering season similar to another hemiparasite plant L. parasiticus. The ecological, taxonomic and complex correlation studies unravelled that the hemiparasites act as the potent reservoirs of endophytic communities throughout the year and during favourable conditions, these bacterial communities disseminate to the mango plant.

Insights into endophytic bacterial diversity of rice grown across the different agro-ecological regions of West Bengal, India

Endophytes have recently garnered importance worldwide and multiple studies are being conducted to understand their important role and mechanism of interaction inside plants. But before we indulge in their functions it is necessary to dig into the microbiome. This will help to get a complete picture of the microbes intrinsic to their host and understand changes in community composition with respect to their habitats. To fulfil this requirement in our study we have attempted to dissect the endophytic diversity in roots of rice plant grown across the various agro-ecological zones of West Bengal by undergoing amplicon analysis of their 16S rRNA gene. Based on the measured environmental parameters agro-ecological zones can be divided into two groups: nutrient dense groups, representing zones like Gangetic, Northern hill and Terai-Teesta zone characterised by soil with higher levels of nitrogen (N) and total organic carbon and nutrient low groups representing Coastal saline, Red-laterite and Vindhyan zone mainly characterised by high electroconductivity and pH. Gammaproteobacteria, Alphaproteobacteria, Bacilli and Bacteroidetes were mostly abundant in nutrient dense sites whereas Clostridia and Planctomycetes were concentrated in nutrient low sites. Few genera (Aeromonas, Sulfurospirillum, Uliginosibacterium and Acidaminococcus) are present in samples cultivated in all the zones representing the core microbiome of rice in West Bengal, while some other genera like Lactococcus, Dickeya, Azonexus and Pectobacterium are unique to specific zone. Hence it can be concluded that this study has provided some insight in to the endophytic status of rice grown across the state of West Bengal.

Metagenomic analyses, isolation and characterization of endophytic bacteria associated with Eucalyptus urophylla BRS07-01 in vitro plants

Eucalyptus is the main species for the forestry industry in Brazil. Biotechnology and, more recently, gene editing offer significant opportunities for rapid improvements in Eucalyptus breeding programs. However, the recalcitrance of Eucalyptus species to in vitro culture is also a major limitation for commercial deployment of biotechnology techniques in Eucalyptus improvement. We evaluated various clones of Eucalyptus urophylla for their in vitro regeneration potential identified a clone, BRS07-01, with considerably higher regeneration rate (85%) in organogenesis, and significantly higher than most works described in literature. Endophytic bacteria are widely reported to improve in vitro plant growth and development. Hence, we believe that inclusion of endophytic plant growth promoting bacteria enhanced was responsible for the improved plantlets growth and development of this clone under in vitro culture. Metagenomic analysis was performed to isolate and characterize the prominent endophytic bacteria on BRS07-01 leaf tissue in vitro micro-cultures, and evaluate their impact on plant growth promotion. The analysis revealed the presence of the phyla Firmicutes (35%), Proteobacteria (30%) and much smaller quantities of Actinobacteria, Bacteroidetes, Gemmatimonadetes, Crenarchaeota, Euryarchaeota and Acidobacteria. Of the thirty endophytic bacterial strains isolated, eleven produced indole-3-acetic acid. Two of the isolates were identified as Enterobacter sp. and Paenibacillus polymyxa, which are nitrogen-fixing and capable of phosphate and produce ammonium. These isolates also showed similar positive effects on the germination of common beans (Phaseolus spp.). The isolates will now be tested as a growth promoter in Eucalyptus in vitro cultures. Graphical abstract for the methodology using cultivation independent and dependent methodologies to investigate the endophytic bacteria community from in vitro Eucalyptus urophylla BRS07-01.

Screening of Rhizosphere Bacteria and Nematode Populations Associated with Soybean Roots in the Mpumalanga Highveld of South Africa

Soybean is among South Africa’s top crops in terms of production figures. Over the past few years there has been increasingly more damage caused to local soybean by plant-parasitic nematode infections. The presence of Meloidogyne (root-knot nematodes) and Pratylenchus spp. (root lesion nematodes) in soybean fields can cripple the country’s production, however, little is known about the soil microbial communities associated with soybean in relation to different levels of Meloidogyne and Pratylenchus infestations, as well as the interaction(s) between them. Therefore, this study aimed to identify the nematode population assemblages and endemic rhizosphere bacteria associated with soybean using Next Generation Sequencing (NGS). The abundance of bacterial genera that were then identified as being significant using linear discriminant analysis (LDA) Effect Size (LEfSe) was compared to the abundance of the most prevalent plant-parasitic nematode genera found across all sampled sites, viz. Meloidogyne and Pratylenchus. While several bacterial genera were identified as significant using LEfSe, only two with increased abundance were associated with decreased abundance of Meloidogyne and Pratylenchus. However, six bacterial genera were associated with decreased Pratylenchus abundance. It is therefore possible that endemic bacterial strains can serve as an alternative method for reducing densities of plant-parasitic nematode genera and in this way reduce the damages caused to this economically important crop.

MYCOBACTERIA IN SKIN LESIONS AND THE HABITAT OF THE ENDANGERED HOUSTON TOAD (ANAXYRUS HOUSTONENSIS)

Head-starting of the federally endangered Houston toad (Anaxyrus houstonensis), that is, the release of egg strands, tadpoles, and metamorphic juveniles produced in captivity into the original breeding ponds, requires assessment of potential threats for the transmission of pathogens from captive to free-ranging toads. We used Illumina-based 16S rRNA V3 amplicon sequencing to investigate the community structure of bacteria from skin lesions of captive Houston toad and habitat (pond) samples. Proteobacteria, alone or together with Actinobacteria and, in some samples, Cyanobacteria represented virtually all reads in tissue lesion samples, whereas pond samples were much more diverse, with Acidobacteria, Actinobacteria, Bacteriodetes, Chloroflexi, Cyanobacteria, Firmicutes, Planctomycetes, Proteobacteria, and Verrucomicrobia present with little variation between samples. If present in lesions, Actinobacteria were largely represented by Mycobacteriaceae, and here mainly by one sequence identical to sequences of members of the Mycobacterium chelonae-abscessus complex. In pond samples, mycobacteria represented only a small portion of the actinobacteria, although at higher diversity with six distinct reads. Sequences for reads obtained from pond samples were identical to those representing the M. chelonae-abscessus complex, a group with Mycobacterium marinum, Mycobacterium kansasii, Mycobacterium avium, a group with Mycobacterium vaccae, Mycobacterium fortuitum, Mycobacterium poriferae, and a group with Mycobacterium elephantis and Mycobacterium celeriflavum, whereas sequences of high similarity were detected for reads related to those of Mycobacterium holsaticum, Mycobacterium pallens, and Mycobacterium obuense, and Mycobacterium goodii. Our results indicated that lesions observed on the Houston toad in captivity are not the result of mycobacteria in every case, and that the presence of mycobacteria in the captive colony does not represent a novel pathogen threat to the wild populations because such bacteria are also seen in the natural pond habitats for the Houston toad.

Genomic assessment of Stenotrophomonas indicatrix for improved sunflower plant

Diverse agriculturally important microbes have been studied with known potential in plant growth promotion. Providing several opportunities, Stenotrophomonas species are characterized as promising plant enhancers, inducers, and protectors against environmental stressors. The S. indicatrix BOVIS40 isolated from the sunflower root endosphere possessed unique features, as genome insights into the Stenotrophomonas species isolated from oilseed crops in Southern Africa have not been reported. Plant growth-promotion screening and genome analysis of S. indicatrix BOVIS40 were presented in this study. The genomic information reveals various genes underlining plant growth promotion and resistance to environmental stressors. The genome of S. indicatrix BOVIS40 harbors genes involved in the degradation and biotransformation of organic molecules. Also, other genes involved in biofilm production, chemotaxis, and flagellation that facilitate bacterial colonization in the root endosphere and phytohormone genes that modulate root development and stress response in plants were detected in strain BOVIS40. IAA activity of the bacterial strain may be a factor responsible for root formation. A measurable approach to the S. indicatrix BOVIS40 lifestyle can strategically provide several opportunities in their use as bioinoculants in developing environmentally friendly agriculture sustainably. The findings presented here provide insights into the genomic functions of S. indicatrix BOVIS40, which has set a foundation for future comparative studies for a better understanding of the synergism among microbes inhabiting plant endosphere. Hence, highlighting the potential of S. indicatrix BOVIS40 upon inoculation under greenhouse experiment, thus suggesting its application in enhancing plant and soil health sustainably.

Development of an inexpensive matrix-assisted laser desorption-time of flight mass spectrometry method for the identification of endophytes and rhizobacteria cultured from the microbiome associated with maize

Many endophytes and rhizobacteria associated with plants support the growth and health of their hosts. The vast majority of these potentially beneficial bacteria have yet to be characterized, in part because of the cost of identifying bacterial isolates. Matrix-assisted laser desorption-time of flight (MALDI-TOF) has enabled culturomic studies of host-associated microbiomes but analysis of mass spectra generated from plant-associated bacteria requires optimization. In this study, we aligned mass spectra generated from endophytes and rhizobacteria isolated from heritage and sweet varieties of Zea mays. Multiple iterations of alignment attempts identified a set of parameters that sorted 114 isolates into 60 coherent MALDI-TOF taxonomic units (MTUs). These MTUs corresponded to strains with practically identical (>99%) 16S rRNA gene sequences. Mass spectra were used to train a machine learning algorithm that classified 100% of the isolates into 60 MTUs. These MTUs provided >70% coverage of aerobic, heterotrophic bacteria readily cultured with nutrient rich media from the maize microbiome and allowed prediction of the total diversity recoverable with that particular cultivation method. Acidovorax sp., Pseudomonas sp. and Cellulosimicrobium sp. dominated the library generated from the rhizoplane. Relative to the sweet variety, the heritage variety c ontained a high number of MTUs. The ability to detect these differences in libraries, suggests a rapid and inexpensive method of describing the diversity of bacteria cultured from the endosphere and rhizosphere of maize.

Cultivable and metagenomic approach to study the combined impact of nanogypsum and Pseudomonas taiwanensis on maize plant health and its rhizospheric microbiome

In the present study we examined the effect of nanogypsum and Pseudomonas taiwanensis strain BCRC 17751on plant and soil health using conventional and metagenomics approaches. Soil physicochemical properties and agronomical parameters of maize plants were reported to be better when applied with nanogypsum and bacterial inoculum together. When compared to control a significant increase in total bacterial counts, nitrogen, phosphorus, potassium (NPK) solubilizing bacterial population and soil enzyme activities (fluorescein diacetate, alkaline phosphatase, dehydrogenase, β-glucosidase, arylesterase and amylase) was reported in treatments. The metagenomics studies revealed dominance of beneficial bacteria such as Proteobacteria, Bacteriodetes, Planctomycetes, Acidobacteria and Nitrospirae in treated soil. On the other hand some novel bacterial diversity was also reported in treated soil which was evident from presence of taxonomically unclassified sequences. Hence, it can be concluded that combined application of nanogypsum and Pseudomonas taiwanensis in maize help in improving the structure and function of soil which affects the plant health without causing any toxic effect. However, in situ validation of the prescribed treatment is required under field conditions on different crops in order to give maximum benefits to the farmers and the environment.

Insights into the structure and role of seed-borne bacteriome during maize germination

Seed germination events modulate microbial community composition, which ultimately influences seed-to-seedling growth performance. Here, we evaluate the germinated maize (variety SHS 5050) root bacterial community of disinfected seed (DS) and non-disinfected seed (NDS). Using a gnotobiotic system, sodium hypochlorite (1.25%; 30 min)-treated seeds showed a reduction of bacterial population size and an apparent increase of bacterial community diversity associated with a significant selective reduction of Burkholderia-related sequences. The shift in the bacterial community composition in DS negatively affects germination speed, seedling growth and reserve mobilization rates compared with NDS. A synthetic bacterial community (syncom) formed by 12 isolates (9 Burkholderia spp., 2 Bacillus spp., and 1 Staphylococcus sp.) obtained from natural microbiota maize seeds herein was capable of recovering germination and seedling growth when reintroduced in DS. Overall, results showed that changes in bacterial community composition and selective reduction of Burkholderia-related members' dominance interfere with germination events and the initial growth of the maize. By cultivation-dependent and -independent approaches, we deciphered seed-maize microbiome structure, bacterial niches location and bacterial taxa with relevant roles in seedling growth performance. A causal relationship between seed microbial community succession and germination performance opens opportunities in seed technologies to build-up microbial communities to boost plant growth and health.

Genes related to redox and cell curvature facilitate interactions between Caulobacter strains and Arabidopsis

Bacteria play an integral role in shaping plant growth and development. However, the genetic factors that facilitate plant-bacteria interactions remain largely unknown. Here, we demonstrated the importance of two bacterial genetic factors that facilitate the interactions between plant-growth-promoting (PGP) bacteria in the genus Caulobacter and the host plant Arabidopsis. Using homologous recombination, we disrupted the cytochrome ubiquinol oxidase (cyo) operon in both C. vibrioides CB13 and C. segnis TK0059 by knocking out the expression of cyoB (critical subunit of the cyo operon) and showed that the mutant strains were unable to enhance the growth of Arabidopsis. In addition, disruption of the cyo operon, metabolomic reconstructions, and pH measurements suggested that both elevated cyoB expression and acid production by strain CB13 contribute to the previously observed inhibition of Arabidopsis seed germination. We also showed that the crescent shape of the PGP bacterial strain C. crescentus CB15 contributes to its ability to enhance plant growth. Thus, we have identified specific genetic factors that explain how select Caulobacter strains interact with Arabidopsis plants.

A novel PCR-clamping assay reducing plant host DNA amplification significantly improves prokaryotic endo-microbiome community characterization

Due to the sequence homology between the bacterial 16S rRNA gene and plant chloroplast and mitochondrial DNA, the taxonomic characterization of plant microbiome using amplicon-based high throughput sequencing often results in the overwhelming presence of plant-affiliated reads, preventing the thorough description of plant-associated microbial communities. In this work we developed a PCR blocking primer assay targeting the taxonomically informative V5-V6 region of the 16S rRNA gene in order to reduce plant DNA co-amplification, and increase diversity coverage of associated prokaryotic communities. Evaluation of our assay on the characterization of the prokaryotic endophytic communities of Zea mays, Pinus taeda and Spartina alternifora leaves led to significantly reducing the proportion of plant reads, yielded 20 times more prokaryotic reads and tripled the number of detected OTUs compared to a commonly used V5-V6 PCR protocol. To expand the application of our PCR-clamping assay across a wider taxonomic spectrum of plant hosts, we additionally provide an alignment of chloroplast and mitochondrial DNA sequences encompassing more than 200 terrestrial plant families as a supporting tool for customizing our blocking primers.

An endophytic strain of the genus Bacillus isolated from the seeds of maize (Zea mays L.) has antagonistic activity against maize pathogenic strains

Plant endophytes are microbes that colonize plant internal tissues and are ubiquitously associated with plants. In this study, seven endophytic bacterial strains, 665L2, 725L2, 725R2, 92R2, 728R3, 728R4 and 2416T3, were isolated from seeds of five healthy maize varieties (Zea mays L.) and all identified as Bacillus velezensis by polyphasic taxonomy based on 16S rRNA and gyrA gene phylogenetic analysis. In addition, according to the genotyping results from random amplified polymorphic DNA (RAPD), 665L2, 725L2, 725R2 and 92R2 belonged to the same strain, while 728R3 and 2416T3 belonged to another strain. Pathogenic fungal strains 4, 5 and 6 were isolated from three diseased maize varieties (Zea mays L.), and they were identified as Talaromyces funiculosus, Penicillium oxalicum and Fusarium verticillioides, respectively, by polyphasic taxonomy based on morphological identification, ITS rDNA and bio-control gene phylogenetic analyses. Seven endophytic bacterial Bacillus velezensis strains had favourable antagonistic activity, and antagonistic testing was carried out against the three pathogenic strains, Talaromyces funiculosus 4, Penicillium oxalicum 5 and Fusarium verticillioides 6. Biological control lipopeptide antibiotic genes (bioA, bmyB, ituC, fenD, srfAA, srfAB, yngG and yndJ) were amplified using specific primers, and they were found in the seven endophytic bacterial Bacillus velezensis strains. This study provides a scientific basis for future research on the use of resistant endophytic bacterial resources to enhance crop production.

The plant organs and rhizosphere determine the common bean mycobiome

The plant microbiota diversity is often underestimated when approaches developed mainly for the identification of cultivable microorganisms are used. High-throughput sequencing allows a deeper understanding of the microbial diversity associated with plants. The amplification of ITS1 was used to analyze fungal diversity in several plant organs and rhizosphere of three common bean (Phaseolus vulgaris) varieties grown in a greenhouse. The fungal diversity diverged between those plant organs and the rhizosphere, with the highest found in the rhizosphere and the lowest in the stem. In each organ different numbers of genus, OTUs were identified, in a total of 283 OTUs evenly distributed among the varieties. In the co-occurrence network, a larger number of positive interactions were found in the organs of the aerial part in all varieties. We observed that the diversity of the endophytic microbiota differed more between plant organs than between common bean varieties. Our results show that the diversity of endophytic fungi can be efficiently accessed with the sequencing of ITS amplicons and that this diversity may vary among distinct plant organs and the rhizosphere of a single plant variety.

High-throughput sequence analysis of bacterial communities in commercial biofertiliser products marketed in South Africa: an independent snapshot quality assessment

The genetic and predicted functional diversity of bacterial communities in 12 commercial biofertiliser products were evaluated using high-throughput sequencing of the 16S rRNA gene. Proteobacteria, Firmicutes and Bacteroides dominated the bacterial communities, with the genera Pseudomonas, Lactobacillus, Bacillus, Bradyrhizobium and Rhizobium being prevalent. The manufacturer-specified species were detected in relatively high abundance in two of the products while a few or none of the specified species were detected in some products. A number of unspecified microbes were detected, including potential human and crop pathogens such as Alcaligenes, Clostridium, Escherichia-Shigella and Proteus. The functional prediction unravelled high prevalence of enzyme-coding genes such as nitrogenase, NifT, alkaline phosphatase and reductases of nitric oxide, nitrate and nitrite which contribute to nitrogen-fixation, phosphorus solubilisation and degradation of nitrates and nitrites. In addition, toxins such as leukocidin/hemolysin and colicin V protein that cause product quality damage were highly predicted in over 67% of the products. Overall, high-throughput sequence analysis of bacterial communities in biofertiliser products revealed that majority of the products were of poor quality. This result justifies the need for regular quality assessment and improvement in quality control systems during biofertiliser formulation.

High-Throughput Sequence Analyses of Bacterial Communities and Multi-Mycotoxin Profiling During Processing of Different Formulations of Kunu, a Traditional Fermented Beverage

Kunu is a traditional fermented single or mixed cereals-based beverage popularly consumed in many parts of West Africa. Presently, the bacterial community and mycotoxin contamination profiles during processing of various kunu formulations have never been comprehensively studied. This study, therefore, investigated the bacterial community and multi-mycotoxin dynamics during the processing of three kunu formulations using high-throughput sequence analysis of partial 16S rRNA gene (hypervariable V3-V4 region) and liquid chromatography tandem mass spectrometry (LC-MS/MS), respectively. A total of 2,303 operational taxonomic units (OTUs) were obtained across six processing stages in all three kunu formulations. Principal coordinate analysis biplots of the Bray-Curtis dissimilarity between bacterial communities revealed the combined influences of formulations and processing steps. Taxonomically, OTUs spanned 13 phyla and 486 genera. Firmicutes (phylum) dominated (relative abundance) most of the processing stages, while Proteobacteria dominated the rest of the stages. Lactobacillus (genus taxa level) dominated most processing stages and the final product (kunu) of two formulations, whereas Clostridium sensu stricto (cluster 1) dominated kunu of one formulation, constituting a novel observation. We further identified Acetobacter, Propionibacterium, Gluconacetobacter, and Gluconobacter previously not associated with kunu processing. Shared phylotypes between all communities were dominated by lactic acid bacteria including species of Lactobacillus, Lactococcus, Leuconostoc, Pediococcus, and Weissella. Other shared phylotypes included notable acetic acid bacteria and potential human enteric pathogens. Ten mycotoxins [3-Nitropropionic acid, aflatoxicol, aflatoxin B1 (AFB1), AFB2, AFM1, alternariol (AOH), alternariolmethylether (AME), beauvericin (BEAU), citrinin, and moniliformin] were quantified at varying concentrations in ingredients for kunu processing. Except for AOH, AME, and BEAU that were retained at minimal levels of < 2 μg/kg in the final product, most mycotoxins in the ingredients were not detectable after processing. In particular, mycotoxin levels were substantially reduced by fermentation, although simple dilution and sieving also contributed to mycotoxin reduction. This study reinforces the perception of kunu as a rich source of bacteria with beneficial attributes to consumer health, and provides in-depth understanding of the microbiology of kunu processing, as well as information on mycotoxin contamination and reduction during this process. These findings may aid the development of starter culture technology for safe and quality kunu production.

Identification of Major Rhizobacterial Taxa Affected by a Glyphosate-Tolerant Soybean Line via Shotgun Metagenomic Approach

The worldwide commercial cultivation of transgenic crops, including glyphosate-tolerant (GT) soybeans, has increased widely during the past 20 years. However, it is accompanied with a growing concern about potential effects of transgenic crops on the soil microbial communities, especially on rhizosphere bacterial communities. Our previous study found that the GT soybean line NZL06-698 (N698) significantly affected rhizosphere bacteria, including some unidentified taxa, through 16S rRNA gene (16S rDNA) V4 region amplicon deep sequencing via Illumina MiSeq. In this study, we performed 16S rDNA V5–V7 region amplicon deep sequencing via Illumina MiSeq and shotgun metagenomic approaches to identify those major taxa. Results of these processes revealed that the species richness and evenness increased in the rhizosphere bacterial communities of N698, the beta diversity of the rhizosphere bacterial communities of N698 was affected, and that certain dominant bacterial phyla and genera were related to N698 compared with its control cultivar Mengdou12. Consistent with our previous findings, this study showed that N698 affects the rhizosphere bacterial communities. In specific, N698 negatively affects Rahnella, Janthinobacterium, Stenotrophomonas, Sphingomonas and Luteibacter while positively affecting Arthrobacter, Bradyrhizobium, Ramlibacter and Nitrospira.