Microbiome and pathobiome analyses reveal changes in community structure by foliar pathogen infection in rice

Dastogeer, Khondoker M. G.; Yasuda, Michiko; Okazaki, Shin

doi:10.3389/fmicb.2022.949152

ORIGINAL RESEARCH article

Front. Microbiol., 02 August 2022

Sec. Microbe and Virus Interactions with Plants

Volume 13 - 2022 | https://doi.org/10.3389/fmicb.2022.949152

This article is part of the Research Topic Emerging Frontiers of Microbial Functions in Sustainable Agriculture View all 13 articles

Microbiome and pathobiome analyses reveal changes in community structure by foliar pathogen infection in rice

$\nKhondoker M. G. Dastogeer,$ Khondoker M. G. Dastogeer^1,2^*

Michiko Yasuda¹

Shin Okazaki¹^*

¹Plant Microbiology Laboratory, Tokyo University of Agriculture and Technology, Tokyo, Japan
²Department of Plant Pathology, Bangladesh Agricultural University, Mymensingh, Bangladesh

Increasing evidence suggests that the plant rhizosphere may recruit beneficial microbes to suppress soil-borne pathogens, but microbiome assembly due to foliar pathogen infection and ecological mechanisms that govern microbiome assembly and functions in the diseased host are not fully understood. To provide a comprehensive view of the rice-associated microbiome, we compared bacterial and fungal communities of healthy rice and those infected with Magnaporthe oryzae, the causal agent of blast disease. We found that the soil had a greater diversity of bacterial and fungal communities than plant endospheric communities. There was no significant dysbiosis of bacterial and fungal microbiome diversity due to disease, but it caused a substantial alteration of bacterial community structure in the root and rhizosphere compartments. The pathobiome analysis showed that the microbiome community structure of leaf and grain tissues was changed markedly at the pathogen infection site, although the alpha diversity did not change. Correspondingly, the relative abundances of some bacteria and fungi were clearly altered in symptomatic tissues. We noted an increase in Rhizobium bacteria and a decline of Tylospora, Clohesyomyces, and Penicillium fungi in the symptomatic leaf and grain tissues from both locations. According to the inferred microbial network, several direct interactions between M. oryzae and other microbes were identified. The majority of edges in the interaction network were positive in diseased samples; contrastingly, the number of edges was much lower in the healthy samples. With source tracking analysis, we observed a sharp contrast in the source of root endosphere bacteria due to Magnaporthe infection. Whereas the majority (71%) of healthy root bacteria could be tracked from the soil, only a very small portion (17%) could be tracked from the soil for diseased samples. These results advanced our understanding and provided potential ideas and a theoretical basis for studying pathobiome and exploiting the microbiome for sustainable agriculture.

Introduction

The world faces formidable challenges in achieving food security for the ever-increasing population (Mc Carthy et al., 2018; Barrett, 2021). Various biotic and abiotic stressors, such as drought, salinity, disease, and pests affect crop production. In particular, plant pathogens represent an insidious threat to agriculture, and they are accounted for a loss of ~16% of global crop yield (Oerke, 2006). Numerous studies have been conducted on the intricate interplay between plant and pathogen, including pathogenicity, disease progression, plant immunity, and disease management. However, in most of these studies, the main focus of the investigation has been the binary interactions between plant and pathogen and those under various environmental conditions. The enormous complexity of interactions among plant, pathogen, and other microorganisms and their outcome under diverse conditions have also received attention in the recent years (Bulgarelli et al., 2013; Santos and Olivares, 2021; Trivedi et al., 2022).

Similar to humans and other animals, plants harbor complex microbial communities called the “plant microbiome.” There is an increasing interest in understanding the composition and function of the microbiome for harnessing their potential, such as growth promotion and disease resistance of the host, as well as in understanding the basics of host–microbial symbioses (Busby et al., 2017; Song et al., 2020; Zhang et al., 2021). Pathogenic microbes cause changes in plant phenotypes through tissue damage and induction of plant defenses, which can alter the immunity of plants to colonization by microorganisms. Therefore, factors that influence the impact of the pathogen on hosts will likely affect the colonization and growth of plant-associated microorganisms. Various host and environmental factors influence microbiome structure and diversity in the plant (Compant et al., 2019; Dastogeer et al., 2020). Host immunity level is one of the major factors shaping plant microbiome community composition (Dastogeer et al., 2020). Recent studies suggest that similar to the gut microbiome, the plant microbiome can enhance the immune functions of the plant host (Vannier et al., 2019; Ma et al., 2021; Teixeira et al., 2021). It has been shown that plants can recruit selective microorganisms from the soil by exudating compounds in the rhizosphere to positively impact plant growth and health (Reinhold-Hurek et al., 2015; Sasse et al., 2018; Liu et al., 2021).

Accumulating data implies that there remains a constant battle between the host and its microbes to maintain microbiome homeostasis in the host (Paasch and He, 2021). In the human microbiome research, the concept of “healthy microbiome” has gained massive attention, although the definition of “healthy microbiome” is still not very clear-cut (Bäckhed et al., 2012; Shanahan et al., 2021). In addition to modifying host defense, the healthy microbiome maintains ecological stability in the host and thus prevents microbial intruders, such as those protected from pathogen attack. On the other hand, the disease causes a shift in the microbiome termed, “microbiome dysbiosis,” a situation in which microbiome homeostasis is disrupted and the organism becomes more vulnerable to potentially harmful microbial invaders. The pathobiome concept arose from human studies, which suggests a shift of “one pathogen–one disease paradigm” to a set of host-associated microorganisms with reduced (or possibly reduced) health status as a result of interactions between members of that set and the host (Defazio et al., 2014; Krezalek et al., 2016). To understand the pathobiome, it is important to clarify the nature of the interactions of the associated microbes among themselves and with the host, in addition to the identity of the community (Bass et al., 2019). It was shown that perturbations in the microbial community by external factors could significantly affect the susceptibility of humans and animals to several diseases (Sekirov et al., 2010; Ferreira et al., 2011; Willing et al., 2011; Vannier et al., 2019). For example, during the progression of mastitis, dysbiosis of the milk microbiome can occur with the increase of opportunistic pathogenic bacteria and reduction of healthy commensal bacteria (Patel et al., 2017). However, it is mostly unclear whether dysbiosis is a cause or consequence of disease (Bäckhed et al., 2012). An in-depth understanding of microbiome dysbiosis may be helpful in the effort to restore a microbial community so as to abate the host damage. There is a paucity of information in the plant system regarding if and to what extent dysbiosis of the microbiome occurs due to pathogen infection in the host. In a study, Kaushal et al. (2020) showed that when Fusarium oxysporum infects banana plants, the fungal and bacterial communities had a shift toward a less diverse community. Humphrey and Whiteman showed that bittercress affected by a leaf-mining fly, Scaptomyza nigrita, had overall higher bacterial densities than undamaged leaves, mainly due to the increased abundance of typical leaf bacteria, particularly Pseudomonas (Humphrey and Whiteman, 2020). Increased bacterial diversity in the damaged tissue is likely caused by increased released nutrients, increased jasmonic acid, and bacterial inoculation by insect secretions (Humphrey and Whiteman, 2020; Smets and Koskella, 2020). Increased understanding of microbiome alterations under stress and functions of microorganisms in improving plant fitness and defenses highlighted the need for the investigation of their role more elaborately across various plant-pathogen-environment studies.

Rice (Oryza sativa) is an important crop that constitutes the staple diet of over three billion people worldwide (Skamnioti and Gurr, 2009). Rice blast disease caused by Magnaporthe oryzae is a global problem that can cause 10–30% reduction in rice production each year, which could feed 60 million more people (Skamnioti and Gurr, 2009; Kirtphaiboon et al., 2021). Due to a broad host range of the pathogen and the evolution of new pathotypes, blast management is a daunting task (Valent, 2021; Devanna et al., 2022). Also, environmental sustainability necessitates the innovation of natural biocontrol agents in place of chemical fungicides. Therefore, rice–Magnaporthe interaction pathosystem has emerged as a model system to study host–pathogen interaction. The microbiome can play a significant role in host defense and pathogen infection, and the study needs to consider pathobiome for a better understanding of the roles of disease-associated microbes. This study compared the rice microbiomes of healthy and infected rice samples gathered from the same location to gain an insight into the potential role of the pathogens in shaping the microbiome composition. We investigated the microbiome (both bacteria and fungi) structure of non-symptomatic and symptomatic rice plants naturally infected by M. oryzae. Symptomatic rice plants were selected based on the morphological symptoms, including typical eye-shaped brown spots in leaves and yellowish lines surrounding the spots. In addition, we demonstrated the composition and assemblage of the naturally occurring microbiome in the bulk soil, rhizosphere, root, leaf, and grain samples. We believe our study is one of the pioneer investigations to describe an inventory of bacterial and fungal communities associated with the components of non-symptomatic and symptomatic rice plants infected by M. oryzae in rice.

Materials and methods

Headings collection and processing of soil, rhizosphere, and plant samples

Rice (Oryza sativa japonica) samples were collected from two locations, viz. Fukushima and Miyagi of Japan. Samples of Miyagi (cultivar: 54-3110) were collected from a field of Miyagi Prefectural Furukawa Agricultural Experiment Station (38.59652N, 140.91219E) on 9 September 2019, and the Fukushima samples (cultivar: Hitomebore) were collected from Nihonmatsu city (37.602373N, 140.586406E) grown in a farmer's field on 21 September 2019. The cultivars were susceptible to blast infection, and there were severe disease symptoms throughout the fields in both locations. In both places, samples were collected at the maturity stage of rice before the ripening phase. Rice plants were dug out with soil from around 20 cm using a shovel. We randomly chose and collected nine healthy-looking plants and nine blast-infected plants in different areas within a plot. Each plant was wrapped separately with a plastic bag and transported to the laboratory in a cooler box containing dry ice to maintain a low temperature to minimize potential microbial community disruptions. In the laboratory, the samples were further kept at 4°C until processed, and processing was completed within 48 h of collection. Sample fractionation into bulk soil, rhizosphere soil, the root, leaf, and grain compartments were performed within 48 h after sampling. At first, plants were cut above the root system. From the shoots, leaf and panicle were separated and washed with water to remove the adhering debris and named “leaf fraction” and “grain fraction,” respectively. The roots systems were shaken to collect soil manually, and the soil was labeled as “bulk soil fraction.” Large soil clods were broken with a hand-held tiller. All the debris was removed from the soil, placed in the zipper storage bag, and stored at −20°C until use. Using sterilized scissors, rice roots around 6–8 cm in length were cut. The excised roots (cut as needed to fit) were placed in a 50 ml falcon tube containing 35 ml of autoclaved phosphate buffer (6.33 g/L NaH₂PO₄, 8.5 g/L Na₂HPO₄ anhydrous, pH = 6.5, 200 μl/L silwet L-77 as surfactant). The tubes were shaken for 2 min to release the rhizosphere from the surface of the roots (rhizosphere fraction). The roots were taken out with sterilized forceps, blot dried on paper towels, and placed in a new 50 mL falcon tube (root fraction). The roots, leaves, and grains were surface sterilized by washing in 0.25% of NaOCl for 1 min, followed by 70% of EtOH for 40 s and subsequent washing in sterile water thrice. The efficacy of surface sterilization was evaluated by tissue imprint method (Greenfield et al., 2015). We used a cork borer to pierce and collect small tissues from the middle of the symptoms of leaves and named “symptomatic fraction” and from non-symptomatic tissues that did not show any apparent disease symptoms (non-symptomatic fraction). Sterilized tissues were blot dried in autoclaved paper towels, cut into small pieces, and stored at −80°C until further processed. Any plant parts or debris were removed from bulk and rhizosphere soil samples. After suspending in phosphate buffer, the soil was filtered through a sterile 100-μm-mesh cell strainer to remove any small plant parts and debris. Suspended soils were collected by centrifugation at 3000 × g for 5 min, and the pellets were stored at −20°C until DNA extraction. Rhizosphere soils and plant parts from three of the nine plants were bulked to make a biological replicate.

DNA extraction and amplicon analysis

Total genomic DNA from bulk rhizosphere and soil was extracted using NucleoSpin^® Soil (Macherey–Nagel, Duren, Germany) and from the plant tissues using DNeasy^® Plant Mini Kit (QIAGEN, Hilden, Germany) according to respective protocols in the manuals. The DNA samples were eluted in 50 μL of nuclease-free water and used for bacterial and fungal community profiling. The quantity and quality of DNA were measured using a Nanodrop 2000, diluted to 100 ng/ml, and stored at −20°C. Bacterial 16S rRNA using 515f/806r primer pair (515f: 5′- GTGCCAGCMGCCGCGGTAA-3′; 806r: 5′-GGACTACNVGGGTW TCTAAT-3 (Caporaso et al., 2011; Apprill et al., 2015) and the fungal ITS2 using ITS1f/ITS2 primer pair (ITS1F: 5′- GTGAATCATCGAATCTTTGAA-3′; ITS2R: 5′-TCCTCCGCTTATTGATATGC-3′) (White et al., 1990; Turenne et al., 1999) were PCR-amplified. Each sample was amplified in triplicate in a 20 μl reaction volume containing 0.2 μl of Ex-Taq DNA Polymerase, 2.0 μl of 10x Ex buffer (Thermo Fisher Scientific Inc., Waltham, MA, USA), 1 μl of each primer, 13.2 μl of Milli- Q water, and 1 ng of template DNA. PCR was performed with (94°C/3 min, 94°C/45 s, 50°C/60 s, 72°C/90 s, 72°C/10 min for 35 cycles) for 16S rRNA and (94°C/3 min, 94°C/40 s, 55°C/40 s, 72°C/60 s, 72°C/7 min for 30 cycles) for ITS2. PCR quality was controlled by loading 5 μl of each reaction on a 1% agarose gel and affirming that no band was detected within the negative control. The replicated reactions were combined to make one biological replication, and the barcoded Illumina libraries were sent for paired-end Illumina MiSeq sequencing (2 × 300 bp, Bioengineering Lab. Co., Sagamihara, Japan). The obtained 16S rRNA and ITS amplicon raw reads were deposited into the Sequence Read Archive (SRA) database of the National Center for Biotechnology Information (NCBI) under the project PRJNA824966.

Sequencing data processing and identification of amplicon sequence variants

Microbiome bioinformatics was performed with QIIME2 2020.8 (Bolyen et al., 2019). The adapter and primers were removed with Cutadapt v2.4 from the raw reads (Martin, 2011). The sequences were demultiplexed using the q2-demux plugin and followed by quality control, length trimming, denoising, chimera, and PhiX removal, and feature table construction by DADA2 with default settings except that “–ptrunc-len-f” and “–p-trunc-len-r” which were set at 250 and 200, respectively for 16S data and at 160 and 200 for ITS data, respectively (via q2-dada2) (Callahan et al., 2016). The resulting amplicon sequence variants (ASVs) were aligned with mafft (Katoh et al., 2002) and phylogenetic trees were constructed using fasttree2 (Price et al., 2010). Taxonomy was assigned to ASVs using the QIIME feature-classifier classify-skarn (Bokulich et al., 2018) with the pre-trained naïve Bayes SILVA classifier v132 trimmed to the V4 region of the 16S rDNA gene (Quast et al., 2012) for bacteria and pretrained UNITE ver8 99% database (UNITE Community, 2019), trained on the full reference sequences without any extraction for fungi. Non-bacterial and fungal reads were removed from the obtained ASV table. We normalized the library using scaling with ranked subsampling using the “SRS”-function in the “SRS” with “qiime srs SRS” (Beule and Karlovsky, 2020). The alpha diversity were evaluated with the chao1 estimator (richness) and the Shannon index and Faith PD in qiime2 and Kruskal–Wallis tests were used to compare the diversity. Box plots to display the alpha diversity indices were created using ggplot2 (Wickham et al., 2016) installed in R Core Team (2020). The Bray–Curtis dissimilarity matrix was calculated to assess differences in composition of bacterial and fungal communities and used for the principal coordinate analysis (PCoA) ordination plot in PAST4.04 (Hammer et al., 2001). We compared differences in community composition and structure among sampling sites (α- diversity) using the analysis of similarity (ANOSIM) (Clarke, 1993) and permutational multivariate analysis of variance (PERMANOVA) (Anderson, 2014). Least discriminant analysis (LDA) and random forest were performed using microbiomeanalyst (Chong et al., 2020).

Functional prediction, source tracking, and network analyses of fungal and bacterial community

The marker genes in the samples (16S rRNA and ITS) were inferred for the prediction of functional profiles by Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt2) algorithm software (Douglas et al., 2020). MetaCyc, the comprehensive reference gene function, was mapped to analyze gene functions for bacterial and fungal ASVs (Krieger et al., 2004). Pathways inferred with high-level functions were mapped against the MetaCyc database (Caspi et al., 2020). To predict the ecological function of fungi, we used the FUNGuild database (Nguyen et al., 2016), whereas, for bacteria, we used the FAPROTX (Louca et al., 2016). To estimate the association between bacterial and fungal communities in soil, the rhizosphere, roots, leaves, and grains tissues in healthy and diseased plants, we used SourceTracker2 in python (version 3.7.0) (Knights et al., 2011). We analyzed SourceTracker2 default parameters, and one by one, each environmental sample type was designated as a sink with all other environmental sample types as sources. For example, to investigate source-associated ASVs on leaves of healthy plants, root, rhizosphere soil, and bulk soil of healthy samples were specified as sources, and leaf was designated as the sink. Kruskal–Wallis tests were used to compare the mean percentages assigned to different sources. The percentage value was derived from the statistical average of the results of SourceTracker. We constructed networks of microbiome communities in healthy and infected samples using Pearson correlation with a significance of p < 0.05 and correlation coefficient R > 0.60. For this, we combined the ASV table of both bacterial and fungal communities. To simplify the network, we extracted only the ASVs that showed a significant correlation with M. oryzae and created a subnetwork of interactions between the pathogen and other microbes. Although different DNA extraction kits used for plant and soil samples might cause some variation, we considered the impact to be insignificant for our analysis.

Results

Soil microbial communities are more diverse and distinct from endospheric communities

We analyzed and compared the alpha and beta diversity of bacterial and fungal microbiome of different compartments (i.e., soil, the rhizosphere, root, and leaf) of healthy rice plants from both locations. As expected, the rhizosphere and bulk soil microbiome were more diverse, followed by the root and leaf endophytic communities in both locations (Figures 1A–H; Supplementary Figures 1A,B). Similar patterns were observed for other diversity indices, such as Chao1 and Faith's phylogenetic diversity (PD) (Supplementary Table 1). Samples from two different locations, i.e., Fukushima and Miyagi, showed variation in the diversity of bulk soil and leaf infosphere communities (Supplementary Tables 1, 2). The PCoA revealed distinct bacterial and fungal assemblages among different compartments in the healthy samples. It was noted that bulk soil and rhizosphere microbial communities overlapped but were separated from endosphere (the root and leaf) communities which were significantly different from each other (Supplementary Figures 1C,G; (Supplementary Table 3). We checked which bacterial groups differed among compartments using LDA. With few others, Delta and Gamma-proteobacteria, Bacteroidia, Verrucomicrobiae, Acidobacteria, and Clostridia were higher in both soil and rhizosphere samples but very low or absent in plant microbiomes (root and leaf). Oxyphotobacteria is the dominant bacteria in the leaf and root microbiome, which were less abundant in the belowground microbiome (Supplementary Figure 1D). We also checked which fungal groups were different among compartments by the illustration of taxa bar plots. Eurotiomycetes, Dothideomycetes, and Agaricomycetes were more abundant in the aboveground compartment (root and leaf), whereas Sordariomycetes and unknown fungal taxa were more in bulk soil and rhizosphere communities (Supplementary Figure 1H).

FIGURE 1

Figure 1. Bacterial and fungal microbiome community of rice as influenced by plant health status and location: Alpha diversity (Shannon Entropy) of bacterial (A–D) and fungal (E–H) community in bulk soil (A,E), rhizosphere soil (B,F), root endosphere (C,G), and leaf endosphere (D,H). Principal co-ordinate (PCoA), analysis of similarity (ANOSIM), and permutational multivariate analysis of variance (PERMANOVA)-based analyses of bacterial (I–L) and fungal (M–P) microbiome of healthy and diseased (non-symptomatic) samples of rice collected from two locations in Japan. Bray-Curtis coefficient of community similarity index was computed for comparing the bacterial community similarity among groups. The ANOSIM statistic R values (up to 1) and PERMANOVA pseudo-F indicate similarity/dissimilarity between groups. Ordination plots show the influence of plant health status, and location on microbial community assembly of bulk soil (I,M), the rhizosphere soil (J,N), root endosphere (K,O), and the leaf endosphere (L,P). Differential abundance of bacterial microbiome community of rhizosphere soil (Q) and root endosphere (R) rice as obtained by random forest analysis. Random forest analysis results of top 20 bacterial family with the highest discriminatory power between diseased and healthy samples of rhizosphere community are shown. Red fields show a high abundance and blue fields a low abundance of the particular bacterial family. FD, Fukushima diseased; FH, Fukushima healthy; MD, Miyagi diseased; MH, Miyagi healthy.

Disease-caused alteration of bacterial but not fungal community compositions in the rhizosphere and root samples

Since microbiome community composition and diversity were different in different compartments, we compared the bacterial and fungal microbiome of healthy and diseased samples separately from soil, the rhizosphere, root, and leaf. No significant difference was observed for the alpha diversity of bacteria (Shannon, Faith'sPD, and Chao1) between healthy and diseased samples regardless of the microbiome compartments (Figures 1A–D; Supplementary Table 1). PCoA, PERMANOVA, and ANOSIM analyses revealed that the rhizosphere and root compartment had distinctly different bacterial compositions in the diseased plants compared to the healthy samples (Figures 1J,K). The bacterial communities of bulk soil and leaf endophyte did not differ due to disease status. The bacterial community of the two locations was distinctly different regardless of the plant compartment. Several bacterial families were more abundant in the rhizosphere of healthy as compared to the diseased samples, such as Barnesiellaceae, Desulfobacteraceae, Thermomonosporaceae, Methanosaetaceae, Pirellulaceae, and Coriobacteriaceae in both locations. In contrast, bacteria of the families Brevinemataceae, Holophagaceae, Paracaedibacteraceae, and Solirubrobacteraceae were highly abundant in the rhizosphere communities of diseased samples (Figure 1Q). In the root communities, members of Streptomycetaceae, Chitinophagaceae, Burkholderiaceae, and Kineosporiaceae were present abundantly in the healthy samples, and Veillonellaceae were abundant in the diseased samples (Figure 1R). At the genus level, Turneriella, Pelolinea, Sulfuricurvum, Desulfocapsa, Cuspidothrix, and Arenimonas were more abundant in the healthy rhizosphere, and Desulfovirga, Roseiarcus, and Nakamurella in diseased samples (Supplementary Figure 2). In roots, Streptomyces, Methylosinus, Cephaloticoccus, Burkholderia, Paraburkholderia, Caulobacter, Acinetobacter, Geothrix, Ruminiclostridium, Hydrogenophaga, and Pleomorphomonas were more abundant in healthy plants in both locations (Supplementary Figure 2).

In the case of the fungal microbiome, no significant differences were observed for alpha diversity between healthy and diseased samples regardless of the microbiome compartments (Figures 1E–H; Supplementary Table 2). Fukushima samples had a higher diversity of fungal community compared to Miyagi samples which were very prominent in soil communities. However, the richness index (Chao1) did not exhibit any variation between locations except for bulk soil (Supplementary Table 2). From PCoA, we did not observe any influence of blast disease infection on the structure of the fungal community; there were locational variations in the soil as well as endospheric communities (Figures 1M–P; Supplementary Table 2).

Pathobiome analysis of symptomatic and non-symptomatic tissue of diseased plants

We were interested in the microbiome community in symptomatic and non-symptomatic tissues (leaf and grain) from the diseased plants. Shannon alpha diversity analysis indicated no substantial differences in the bacterial and fungal communities between symptomatic and non-symptomatic leaf and grain tissues (Figures 2A–D). PCoA analysis and similarity indices (PERMANOVA, ANOSIM) statistics suggested that there remained clear differences in the bacterial and fungal communities of symptomatic leaves and grains from non-symptomatic tissues and that sampling locations significantly influenced community structuring of leaf and grain communities (Figures 2E–H, Supplementary Tables 4, 5). A further inspection of bacterial group differences using random forest analysis revealed that bacteria belonging to certain families were variably present in non-symptomatic or symptomatic tissues. For instance, the symptomatic leaf tissues had higher abundances of Xanthomonadaceae, Sphingobacteriaceae, Beijerinckiaceae, Microbacteriaceae, and Rhizobiaceae, and the symptomatic grains were dominated by Sphingomonadaceae, Burkholderiaceae, and Caulobacteriaceae families. Enterobacteriaceae were found to be highly abundant in the non-symptomatic grains in both the location and the non-symptomatic leaf of Miyagi samples and of the Fukushima samples (Supplementary Figure 3). At the genus level, several bacteria were detected in significantly higher abundance in the non-symptomatic leaves in the Fukushima samples, but they were either lowly present or absent in the Miyagi sample. However, consistently higher abundances of Stenotrophomonas, Pedobacter, Methylobacterium, and Rhizobium were observed in the symptomatic leaves in both locations (Figure 2I). The bacterial genera Rhizobium, Novosphingobium, Darnella, Brevundimonas, Curtobacterium, Stenotrophomonas, and Acinetobacter had comparatively increased occurrences in the symptomatic grain tissues in both locations. Only Phenyracillin was present abundantly in non-symptomatic grain samples irrespective of the sample location (Figure 2J). As for fungi, members of certain families encountered higher in the non-symptomatic tissues. For instance, Thelephoraceae, Lindgomycetaceae, Atheliaceae, Cantharellaceae, and Aspergillaceae showed increased abundances in leaf and grain non-symptomatic samples (Supplementary Figures 3C,D). The fungal genera, such as Coniochaeta, Tylospora, Tomentella, Clohesyomyces, and Penicillium were found abundantly present in non-symptomatic tissues in both the leaf and grain samples regardless of the sampling location (Figures 2K,L).

FIGURE 2

Figure 2. Bacterial and fungal microbiome community of rice leaf and grain as influenced by fungal pathogen infection: Alpha diversity (Shannon Entropy) of bacterial (A,B) and fungal (C,D) community in leaf (A,C) and grain (B,D) tissues. PCoA, ANOSIM, and PERMANOVA-based analyses of bacterial (E,F) and fungal (G,H) microbiome of symptomatic and asymptomatic samples of leaf (E,G) and grain (F,H) samples of rice collected from two locations in Japan. Random forest analysis results of bacterial (top 20) (I,J) and fungal (K,L) genera with the highest discriminatory power between symptomatic and asymptomatic tissues of leaf (I,K) and grain (J,L). Red fields show a high abundance and blue fields a low abundance of the particular bacterial genera. FNS, Fukushima non-symptomatic; FS, Fukushima symptomatic; MNS, Miyagi non-symptomatic; MS, Miyagi symptomatic.

Functional prediction: Microbiota functional profiling of microbiome of rice

The potential metabolic functions of bacterial and fungal communities in symptomatic and non-symptomatic tissues were predicted by PICRUSt2. The partial least squares-discriminant analysis (PLS-DA) plots were determined to evaluate the similarity of the microbial functions between the symptomatic vs. non-symptomatic and diseased vs. healthy samples. The PLS-DA analyses revealed an evident clustering among the samples. There was a difference in the predicted functional pathways of bacterial microbiome communities in diseased rhizosphere and root endosphere (Supplementary Figure 4). The number of abundant functional metabolic pathways predicted from the rhizosphere and root microbiome was higher in diseased plants. For example, “ethylmalonyl-CoA pathway,” “Amine and Polyamine Degradation,” “Amino Acid Biosynthesis,” “isopropanol biosynthesis,” and several other bacterial metabolic pathways were predicted to be abundant in rhizospheres as well as the root endosphere community of diseased plants as compared to healthy plants (Figure 3A; Supplementary Figure 4). Higher abundance of certain other metabolites was predicted in the microbiome rhizosphere of healthy plants, such as “methyleketone biosynthesis,” “Glycan Biosynthesis,” Secondary Metabolite Biosynthesis,” “Cofactor, Prosthetic Group, Electron Carrier, and Vitamin Biosynthesis,” “Glycolysis,” but interestingly, these were comparatively more abundant in the roots of diseased plants (Figure 3A, Supplementary Figure 4). Significant differences in the functional pathways in both bacterial and fungal communities were observed in the symptomatic and non-symptomatic leaf and grain tissues (Figures 3B,C, Supplementary Figure 4). The predicted bacterial functional metabolic pathways, such as “TCA cycle,” “superpathway of glycolysis, pyruvate dehydrogenase, TCA, and glyoxylate bypass,” “fatty acid and lipid biosynthesis,” “fatty acid and lipid degradation,” “amino acid degradation,” and “amine and polyamine biosynthesis” were more abundant in symptomatic tissues in both grain as well as leaf communities (Figures 3B,C). When analyzed the fungal metabolic predictive pathways, certain pathways similar to those of bacterial pathways were found to be more abundant in the symptomatic tissues, such as “fatty acid and lipid degradation” and “amino acid degradation.” The pathways, such as “carbohydrate biosynthesis,” “pentose phosphate pathways,” “TCA cycle,” “inorganic nutrient metabolism,” and “degradation/utilization/ assimilation–other” were more abundant in the fungal metabolites of the non-symptomatic tissues as compared to symptomatic tissues (Supplementary Figure 4). However, some of these pathways predicted from the bacterial community seemed less abundant in the non-symptomatic tissues (Figures 3B,C).

FIGURE 3

Figure 3. Heatmap Clustering the relative abundance of functional signatures of bacterial microbiome of rice from (A) the rhizosphere, (B) leaf, and (C) grain samples as predicted by Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt2) within Meta-Cyc categories. S and NS represent symptomatic and non-symptomatic tissues, respectively.

Network of interactions of M. oryzae and other microbial members in the rhizosphere, leaf, and grain microbiome

To gain a deeper insight into the interactions of M. oryzae with other microorganisms, the networks in healthy and infected plants were visualized. The network structures of rhizosphere communities appeared to be significantly altered in M. oryzae-infected rice plants showing higher complexity and connectivity than those of healthy plants (Figures 4A,B). We noted more links (43 links) between Magaporthe and other microbial genera in the network of infected rhizosphere compartments than the healthy rhizosphere network (17 links), indicating a shift in the rhizosphere microbiome interactions due to foliar infection. The network of interactions in the infected plants revealed that M. oryzae were positively correlated with almost all bacterial genera belonging to Bacteroidetes, Proteobacteria, Firmicutes, Actinobacteria, Spirochaetes, and Fusobacteria, and negatively correlated with Desulfobacca (Proteobacteria) and a fungal genus, Penicillium. In contrast, most bacterial genera were negatively correlated with M. oryzae in diseased plant rhizo-microbiome (Figures 4A,B). When we compared the network structures of symptomatic and non-symptomatic tissues, we observed changes in the correlation pattern of M. oryzae and other microorganisms (Figures 4C–F). For instance, we noted that there were more links (22 links) between Magnaporthe and other microbial genera in the network of non-symptomatic leaf compartments than the symptomatic network (17 links), indicating a shift in the microbiome interactions due to M. oryzae infection. Similar was the case in the grain network, where there were more links (21 links) between Magaporthe and other microbial genera in the network of non-symptomatic than in the symptomatic network (11 links). Our observations indicate that pathogen invasion alters microbiome members who might assist in colonization through a mutualistic relationship or come to defend the plant against pathogen during the infection process of disease. In other words, microbial members in the non-symptomatic tissues might play roles in keeping the tissues decreased in the invasion of the pathogen that would otherwise cause disease symptoms.

FIGURE 4

Figure 4. Network of interactions between the pathogen (M. oryzae) and other species in the rhizosphere of healthy (A) and diseased plants (B), non–symptomatic (C) and symptomatic leaf (D), and non-symptomatic (E) and symptomatic grain (F). The size of the bubbles corresponds to their relative abundances. The red and blue colors of the edge represent positive and negative interaction with M. oryzae, respectively.

SourceTracker analysis of microbiome community from soil to other compartments

We utilized the SourceTracker program to study the proportion of bacterial and fungal communities derived from soils. According to the source apportionment results, there were no noticeable differences in the sources of rhizosphere bacterial communities between diseased and healthy plants, and around half of the total were from neighboring bulk soil communities (Figure 5). In healthy samples, the majority of root bacteria community members (71%) were derived from the soil communities (bulk and rhizosphere soils), but rare members (17%) of the root communities of the diseased plant were derived from the soil bacteria communities, indicating there is a clear boundary between the interior and exterior of healthy roots (Figure 5). Very rare members (<10%) of the aboveground (leaf and grain) bacteria were sourced from belowground communities, and the majority were from unknown sources irrespective of the plant health status. Regardless of the symptom expression in grains of diseased plants, most bacteria of gain communities were from leaves (~75%). Regarding the fungal communities, the rhizosphere soil communities were mainly (~80%) derived from the bulk soil in both healthy and diseased plants. More than 60% of root endophytic fungal communities were primarily derived from the soil communities (Bulk and rhizosphere). Notably, more than half (60%) of endophytic leaf communities were derived from the belowground in the soil in the healthy plants, which is significantly less (31%) in the non-symptomatic leaf of the diseased plant and very low (11%) in the symptomatic leaf tissues, indicating most of the aboveground fungal species in the healthy plants could be tracked back from the soils but not in the infected plants. Non-symptomatic grain communities derive 40% of fungi from the leaf and 15% from the belowground. In contrast, symptomatic tissues receive the majority (70%) of fungi from the leaf and are very rare (3%) from belowground communities (Figure 5).

FIGURE 5

Figure 5. Percentage of bacteria and fungi in each sample that was associated with source types as determined by SourceTracker. The asterisk (**) indicates significant higher percentage at p < 0.01 in healthy as compared to diseased samples whereas caret (∧) indicates higher percentage at p < 0.05 in non-symptomatic tissues as compared to symptomatic tissues of the diseased plants as obtained from Kruskal–Wallis test.

Discussion

Rice microbiome composition is shaped by compartments

We present a critical appraisal of bacterial and fungal communities sampled in the bulk soil, rhizosphere soil, the root, leaf, and grain in two different locations of non-symptomatic and symptomatic sample types to provide a comprehensive view of the rice-associated microbiome. Rice-associated bacterial and fungal communities have previously been studied by culture-dependent as well as culture-independent methods (Edwards et al., 2015; Bertani et al., 2016; Kanasugi et al., 2020; Kim and Lee, 2020; Sinong et al., 2020). We profiled the 16S and ITS2 regions to reveal the composition of bacterial and fungal communities. Our results support the concept that distinct plant components play a crucial role in engaging microbial communities, irrespective of the location (Supplementary Figure 1). We demonstrated that the soil (bulk and rhizosphere soil) had a greater diversity of bacterial and fungal communities in the rice plant than endospheric (root and leaf) communities. Bacterial and fungal communities of bulk soil and rhizosphere soil showed no significant difference; however, root endospheric microbial communities differed from leaf communities. The rhizosphere has a diverse and highly populated microbiome and is subject to chemical transformations brought about by the root exudates and metabolites of microbial degradation. The rhizosphere is distinct from the edaphosphere (bulk soil zone), with enhanced microbial activity mainly due to enhanced root exudates deposition (Chaudhary et al., 2012; Reinhold-Hurek et al., 2015; Vieira et al., 2020). However, it is interesting that in our analyses, we observed a slightly higher diversity of rhizosphere bacterial community over bulk soil community in Miyagi but not in Fukushima samples. Although several previous studies reported distinct and more diverse rhizosphere bacterial communities in different crops (Praeg et al., 2019; Hinsu et al., 2021) than in the surrounding bulk soil, we found no distinct separation of the bacterial community of these two zones. This could be due to the fact that the rhizosphere zone is dynamic, and the composition of the microbial community changes as a result of changing root exudation patterns that vary during the life cycle of plants (Edwards et al., 2015; Qu et al., 2020). We collected the samples at the grain maturity or near ripening stage, and microbial diversity is less than the active growth stage of the plant (Hinsu et al., 2021). Previous studies reported that there was no considerable difference in the bulk soil and rhizosphere soil microbial communities. Several others also reported reduced microbiome diversity in rhizosphere communities compared with the bulk soil community (Schlaeppi et al., 2014; Essel et al., 2019). Therefore, in addition to root exudation, the composition of the microbial community in the rhizosphere is linked with the plant developmental stage and several biological and environmental factors (Schreiter et al., 2014; Essel et al., 2019). Our analysis showed that fungal communities were not influenced by the rhizosphere environment. We assume that root exudate deposition might be selective, or bacteria might respond more strongly to some compound secreted by the rice plant. The lower richness and diversity of endophytic communities (root and leaf) compared to the soil (bulk and rhizosphere) samples are understandable. This finding is coherent with the reports from several plant species, including Arabidopsis, soybean, rice, Agave, and others (Knief et al., 2012; Dong et al., 2019). Intuitively, the root and leaf tissues had different microbial communities, and the diversity of bacteria in the leaf was lower than that in the root communities. Similar results were also reported in rice in another study by Wang et al. (2016) and other plant hosts, including Arabidopsis (Bodenhausen et al., 2013), tomato (Dong et al., 2019), and Agave species (Coleman-Derr et al., 2016). The higher bacterial richness in the root compartment may be attributed to the fact that it is the primary site of the interaction of plants with soil and represents one of the richest microbial ecosystems on earth (Hardoim et al., 2011; Tian and Zhang, 2017). It has been reported that the diversity of bacterial endophytes decreased as the distance from soil increased (Ottesen et al., 2013; Dong et al., 2019). However, the reverse pattern was found for the fungal diversity, i.e., higher diversity in the leaf than in the root tissues. The previous study also reported higher diversity of fungi in the shoots than in the roots and attributed to the assumption that most endophytic fungi are derived from aerial fungal spores (horizontal transmission), whereas most of the endophytic bacteria are derived from soils (Rodriguez et al., 2009; Wang et al., 2016).

Foliar fungal infection shifted bacterial but not fungal microbiome in the rhizosphere

It was found that the blast-infected samples had different bacterial community structures compared to the healthy plants in the root and rhizosphere, but there was no apparent difference in the fungal community (Figure 1). Our study indicated that bacteria of several families were enriched differentially in the rhizosphere and the roots of the healthy and diseased plants. In particular, members of Streptomycetaceae were found to be significantly more abundant in the root endosphere of healthy plants in both locations (Figure 3). The members of Streptomyces are famous for their ability to produce various bioactive compounds and play essential roles in the agricultural field through their biological control potential against phytopathogens, including M. oryzae fungi (Procópio et al., 2012; Ser et al., 2016; Law et al., 2017). Endophytic or N-fixing Burkholderia strains show immense potential as a biocontrol, as well as plant growth promotion and bioremediation agent (Compant et al., 2008; Elshafie and Camele, 2021). After risk assessment, these Burkholderia- products were withdrawn from the market, as some strains may pose a risk to human health. However, over the last several years, the number of new Burkholderia species that show plant-beneficial properties and are associated with clinical patients has increased enormously (Eberl and Vandamme, 2016). The Caulobacter genus is often isolated from the endosphere or rhizosphere of various plants (Fitzpatrick et al., 2018; Walters et al., 2018; Luo et al., 2019).

The role of plant-associated microbiota in protection against pathogenic fungi and oomycetes is well documented. Although we did not examine empirically, the differential bacterial assemblage may be linked to the differential compound exudation through the rhizosphere due to plant health status responsible for specific microbial group selection, thus regulating their community compositions. The “cry for help” hypothesis suggests that plants recruit microbial partners to maximize their survival and growth when affected by external stress and is likely a survival strategy conserved across the plant kingdom (Liu et al., 2019, 2020; Gao et al., 2021). However, we do not know if the differential microbiome assemblage is the effect of disease or the cause of disease manifestation. Either hypothesis (cause or effect) or even both simultaneously might be possible, and clarification is warranted with manipulative experiments. We detected the presence of Magnaporthe in the rhizosphere of healthy samples. Interestingly, the co-occurrence networks analysis of Magnaporthe with other taxa revealed that the rhizosphere of healthy plants had a robust and more complex network with a higher number of nodes and edges than the diseased network. Most edges of healthy rhizo-microbiome networks were positively correlated with Magnaporthe, whereas they were primarily negative in diseased plants (Figure 4). Again, there were no fungal taxa correlated with Magnaporthe in the diseased plant. Only one ascomycetous genera were positively correlated with it, suggesting the putative importance of bacterial taxa in plant health status under M. oryzae attack compared to the fungal microbiome. A previous study also highlighted the higher stability of bacterial than the fungal networks in healthy plants. Fusarium infection in Chili pepper decreased the complexity of bacterial networks but increased the complexity of fungal networks (Gao et al., 2021). Our functional prediction analysis indicated that the rhizosphere microbiome of diseased plants putatively carries out more diverse metabolic functions (Figure 3) than those associated with healthy individuals, which may be related to disease-induced microbiome recruitment to defend plants under stress. Metaproteomic and transcriptomic analyses, as well as microbiome manipulation of diseased vs. healthy plants, would provide an increased understanding of those microbial community and their functions.

Pathobiome showed clear separation of the microbiome from those of healthy tissues

The pathobiome concept has been defined as the totality of microbes interacting with a given pathogen species and their influence on pathogenesis (Vayssier-Taussat et al., 2014; Jakuschkin et al., 2016). Characterization of the components of the pathobiome is an important consideration for understanding the pathogenesis, persistence, transmission, and evolution of pathogenic agents (Vayssier-Taussat et al., 2014). Shifts in the composition, richness, and abundance of the microbiome have been shown to occur due to pathogen infection in plants (Jakuschkin et al., 2016; Musonerimana et al., 2020; Mannaa and Seo, 2021). Our results also revealed that the composition of the fungal and bacterial community of rice leaves changed markedly with the M. oryzae infection. Through the species classification, we found not only the ASVs assigned to potentially pathogenic Magnaporthe were rarely observed abundantly in symptomatic tissues but also from healthy leaf and grain, albeit with significantly low abundances (Figure 2). We could not technically affirm all those ASVs assigned to Magnaporthe are pathogenic by just using ITS sequences. The presence of Magnaporthe in the healthy tissue implies either there was a latent infection of the pathogen, or these were non-pathogenic strains of Magnaporthe. Correspondingly, the relative abundances of some bacteria and fungi were clearly altered in symptomatic tissues. There was a significant increase in the relative abundance of members, such as Rhizobiaceae, Microbacteriaceae, Beijerinkiacee, and Xanthomonadaceae in the symptomatic over the non-symptomatic tissues. As for fungi, there were higher abundances of Thelephoraceae, Lingomycetaceae, Atheliaceae, Catharellaceae, and Aspergillaceae families in the non-symptomatic leaves in both locations. These compositional changes could be a cause or consequence of the pathogen invasion. Species of Paenibacillus are well-known pathogen suppressors in planta. In many studies, several Paenibacillus have been identified and used as biocontrol agents with the ability to secrete antibiotics or other antimicrobial proteins and have been applied to prevent and control various plant pathogens, including Magnaporthe (Rybakova et al., 2016; Padda et al., 2017; Yu et al., 2019). In the plant endophytic compartments, the relative abundances of Paenibacillus, Enterobacter, and several other bacteria that are often considered plant-beneficial microbes showed a significant decrease in abundance compared to the healthy samples. These decreases imply that these endophytic taxa might be excluded due to compromised immune systems locally or outcompeted by more successful colonizers. The shift in the microbiome composition in the symptomatic tissue could be due to the degradation of plant tissue (necrotic/decomposed tissue) by the pathogen that results in colonization with different microbiomes or colonization of bacteria on dead fungal hyphae for their nutrition or bacterial endosymbiont of fungi or pathogen teaming up with commensals (Lundberg et al., 2012; Venturi and da Silva, 2012; Tláskal et al., 2016). On the other hand, the relative abundances of several bacteria, such as those belonging to Rhizobiaceae, Xanthomonadaceae, and Sphingobacteriaceae in the symptomatic leaf and grain and a few more only in the grain tissues samples showed a marked increase compared to non-symptomatic samples, suggesting they might be involved in the pathogenesis and have mutualistic relationships with the pathogen, or maybe they are opportunists and could take advantage of different ecological niches created by pathogen invasion (Lundberg et al., 2012; Hu et al., 2020). The relatively higher abundances of Rhizobium in the symptomatic tissues implied that this strain could either be involved in mutualism with the M. oryzae in the disease process or may remain as saprophyte by adopting an oligotrophic lifestyle (Poole et al., 2018). Legume-rhizobium symbiosis is well documented, and much of the rhizobia research has focused on specific parts of the rhizobia–legume symbioses. Our current understanding of rhizobia interaction with non-legumes and its ecology is not comprehensive. Several studies reported rhizobium from the rhizosphere and endosphere of non-legume crops, including rice, and described their roles in plant growth promotion (Peng et al., 2008; Zhang et al., 2014).

According to the inferred microbial network, several direct interactions between M.oryzae and other microbes may occur. In contrast to our expectation, the results suggested that mutualism, facilitation, and commensalism may dominate in the interactions of Magnaporthe with other microbes in leaves. Almost all the edges in the interaction network were indeed positive in the symptomatic tissues whereas only about half of the edges were positive in the network of healthy tissues. The positive interactions between microbes and the fungal pathogen may be accounted for by endosymbiosis, or they may also be accounted for by commensalism, where bacteria may use dead hyphae as a source of nutrients (Hoffman and Arnold, 2010; Venturi and da Silva, 2012). Bacteria may also facilitate fungal infection because they benefit from the changes in the plant metabolism induced by the fungal pathogen (Venturi and da Silva, 2012). Several Pantoea species have been isolated as endophytes without causing symptoms; some isolates are antibiotic producers and have been developed as biocontrol agents for plant diseases (Walterson and Stavrinides, 2015; Jiang et al., 2019). However, this bacterial genus is highly diverse (Walterson and Stavrinides, 2015) and has also been reported to cause disease in rice (Doni et al., 2019). As per the interaction network, there may be a direct, antagonistic (negative) ecological interaction between Magnaporthe and Pantoea (Figure 4). This finding suggests the putative roles of Pantoea to protect against the blast disease. Manipulative experiments are therefore required to validate the isolation and the antagonistic relationship between the two species and decipher the mechanism underlying the interaction.

Soil microbes have the potential to influence plant immune defense and pathogen invasion. Therefore, understanding how to plant endosphere microbes that interact with the soil communities may provide a “road map” to explain the pathogen infection and disease process. Previous studies used the SourceTracker program to estimate the proportions of the microbiome in a given community that comes from potential source environments and have been used to analyze the relationship between disease-associated endosphere microbial communities and rhizosphere communities (Knights et al., 2011; Hu et al., 2020). In this study, we utilized this program to track the source of plant rhizospheric and endospheric microbial communities in healthy vs. diseased plants. Previous studies showed that bacterial communities in the rhizosphere soils were mainly derived from the bulk soil (Mendes et al., 2013; Hu et al., 2020). We found a somewhat difference in that although most fungal communities of the rhizosphere were from bulk, only half of the total bacteria in the rhizosphere were derived from the nearby bulk soils (~80%). There is a sharp contrast between the source of root endosphere bacteria resulting from Magnaporthe infection. Whereas the majority (71%) of healthy root bacteria could be tracked from soil (rhizosphere and bulk combined), there was only a very small portion for diseased samples (Figure 5). However, it is not clear about the sources of most microbes, especially bacteria, in diseased plants. We assume that early in the pathogen infection, due to potential host immunity or “cry for help,” as mentioned before, the microbiome was localized prior to our sampling, and the diseased plant could not keep dynamic recruitment of bacteria from soil whereas the healthy individual could. However, controlled experiments with a defined microbiome (synthetic community) and tracking for microbiome sources could reveal important information.

Taken together, we infer that foliar infection by pathogenic M. oryzae causes a shift in the rhizosphere bacteria. M. oryzae in symptomatic leaf and grain tissues are associated with interactions with various bacteria, including Bacillus, Enterobacter, and several other bacteria, the role of which in the disease process needs to be clarified. The positively interacting bacteria putatively obtain benefits from invading pathogens, which might lead to the migration of many additional bacterial genera into the plant roots and leaves, eventually causing an outbreak of blast disease. These findings will provide potential ideas and a theoretical basis for isolating biocontrol agents controlling the blast disease in rice with further work.

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/, PRJNA824966.

Author contributions

KD: conceptualization, methodology, investigation, formal analysis, software, writing—original draft, and visualization. MY: conceptualization, methodology, writing—review and editing, and project administration. SO: conceptualization, resources, and writing—review and editing, supervision, project administration, and funding acquisition. All authors contributed to the interpretation of the data, critically revised the manuscript, read, and approved the final version before submission.

Funding

This work was supported by the Japan Society for the Promotion of Science (JSPS) Kakenhi (Grant Number 19H02860).

Acknowledgments

We thank Dr. Endo Takashi and Mr. Ohno Tatsuhiro for allowing us to collect rice samples from their fields and providing assistance in sampling.

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.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.949152/full#supplementary-material

References

Anderson, M. J. (2014). “Permutational multivariate analysis of variance (PERMANOVA),” in Wiley Statsref: Statistics Reference Online (Hoboken, NJ: Wiley), 1–15.

Google Scholar

Apprill, A., McNally, S., Parsons, R., and Weber, L. (2015). Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton. Aquatic Microb. Ecol. 75, 129–137. doi: 10.3354/ame01753