Plant Pathol J > Volume 35(6); 2019 > Article
Kim, Do, Cho, Jeong, and Kwak: Comparison of Microbial Community of Rhizosphere and Endosphere in Kiwifruit


Understanding the microbial community and function are crucial knowledge for crop management. In this study, bacterial and fungal community structures both rhizosphere and endosphere in kiwifruit were analyzed to gain our knowledge in kiwifruit microbiome. Microbial community in rhizosphere was less variation than endosphere community. Functional prediction results demonstrated that abundance of saprotrophic fungi was similar in both rhizosphere and endosphere, but potential pathogenic fungi was more abundance in endosphere than in rhizosphere. This finding suggested that maintain healthy soil is the first priority to protect the host plant against biotic stresses.

Plants have association with enormous microbial organisms and concept of holobiont comprises the host plant and its associated microbiota including bacteria, archaea, and fungi (Vandenkoornhuyse et al., 2015). Therefore, conventional ‘the disease triangle’ in plant pathology discipline has include microbiota as a factor in the disease management strategy. Metagenomic approach has made possible to unraveling both culturable and unculturable microbes as well as define microbial community structure. Phytobiome can be divided into rhizospheric, phyllospheric, and endophytic microbes depending on the host plant tissue parts (Quiza et al., 2015). Ottesen et al. (2013) provided a basic description and several distinctive of tomato microbiota. Their results indicated the distinct microbial groups with different plant tissue compartments and similar correlation to the distance from the soil. Cregger et al. (2018) also provided holistic microbial structures in eastern cottonwood (Populus deltoides) as the model perennial plant. Community structure of fungi, bacteria, and archaea were primarily varied to plant compartments regardless of the cottonwood genotypes. In comparison microbial structure in domestic and imported kiwifruit pollens, the imported pollens showed the lower diversity than the domestic pollens but a bacterial canker pathogen, Pseudomonas syringae pv. actinidiae, was only detected in the imported pollen (Cho et al., 2018). Up to now, only limited information in kiwifruit microbiome is available, therefore, we compared rhizosphere and endosphere microbial communities in kiwifruit.
Experimental kiwifruit trees (cv. Deliwoong, 10 years olds) were designated in Namhae Fruit Tree Genetic Resource Center in Republic of Korea (34°48′56.8″N, 127°55′42.7″E). Respective samples of rhizosphere soil and sap were collected from the three independent trees as replications. Rhizosphere soil was collected from main root and rootlet (depths between 25 to 40 cm) by a trowel and a brush into 50 ml conical tubes on March 24, 2017 (daily temperature 9.3 ± 4°C and a relative humidity of 68.8%). Sap was taken by drilling trunks and securing the taps into 50 ml conical tubes on April 24, 2017 (temp. 15.9 ± 6°C and relative humidity of 28.8%). DNA of Rhizosphere soil (10 g) and sap (10 ml) were extracted using Power-Max soil DNA isolation kit (Qiagen, Hilden, Germany) by three repetitions. Internal transcribed spacer (ITS) and 16S rRNA libraries were sequenced by Macrogen Co. (Seoul, Korea) with Illumina Miseq 250 paired end platform. Primers ITS3 (TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGGCA TCGATGAAGAACGCAGC) and ITS4 (GTCTCGTGGGCTCGGAGATGTGTATAA GAGACAGTCCTCCGCTTATTGATATGC) for fungal ITS library, 341F (CCTACGGGNGGCWGCAG) and 805R (GACTACHVGGGTATCTAATCC) for bacterial 16S rRNA libraries were employed. Raw sequences were processed quality control and clustering to operational taxonomic units (OTUs) by R version 3.5.2 package DADA2 version 1.10 (Callahan et al., 2016). Base on DADA2 guideline (, 16S rRNA libraries without reverse reads were tunicate to 220 bp according to Phred quality score of 30 or higher and the ITS libraries were not truncated in consideration of length variation depending on fungal diversity. The truncated and non-truncated each sequence were clustered to OTUs with divisive amplicon denoising algorithm (DADA) (Callahan et al., 2016) and removed chimera. The OTUs were assigned taxonomy by Silva version 132 (Quast et al., 2012) for bacteria and UNITE version 8.0 (Kõljalg et al., 2013) for fungi. The bacterial OTUs reads were normalized based on known 16S rRNA copy and predicted enzyme activity and metabolic pathway using PiCRUSt2 version 2.0.4-b (Langille et al., 2013) which re-developed version of PiCRUSt. Fungal OTUs were predicted functional groups by FUNGuild version 1.1 (Nguyen et al., 2016). The metagenomic analyses were performed with R package phyloseq version 1.22.3 (McMurdie and Holmes et al., 2013) and DEseq2 version 1.22.2 (Love et al., 2014). Each graph of metagenomic analysis was visualized by ggplot2 version 3.1.0 in R (Wickham, 2016).
Illumina sequencing data revealed 467,735 and 1,008,735 filtered read numbers for fungi and bacteria, respectively (Supplementary Table 1). Rarefaction curve of each sample in both bacteria and fungi was satisfied to analyze the microbial communities (Supplementary Fig. 1). Principal coordinates analysis and α-diversity showed that rhizosphere community had a relatively similar structure than those of endosphere (Figs. 1 and 2). Both bacteria and fungi diversities were also observed more variation in the endosphere than the rhizosphere. Among top fungal genera, all fungi were assigned into 3 phyla (Table 1) and Ascomycota (63%), Mortierellomycota (21%), and Basidiomycota (11%) were detected in both rhizosphere and endosphere of kiwifruit. Abundance bacterial family also were classified into 3 phyla, including Proteobacteria (51%), Actinobacteria (16%), and Firmicutes (13%) (Table 2). As expected diversity of kiwifruit microbiota, the rhizosphere community showed more microbial diversity than endosphere in both bacteria and fungi (Supplementary Fig. 2). Additionally, shared OTUs between rhizosphere and endosphere presented that only small group of microbes were existed in both rhizospheric soil and sap. This finding suggested that among numerous microbes in rhizosphere, certain microbes may enter into the endosphere as either entophytic or invasion pathogen. List of common fungal OTUs in both rhizosphere and endosphere (Supplementary Table 2), Didymella bellidis and Cladosporium cladosporioides were detected in rhizosphere as 1.8% and 0.4%, respectively. However, abundancy of D. bellidis and C. cladosporioides in sap were 16% and 8.1%, respectively. D. bellidis was previously classified as Phoma beliidis, which is well known the plant pathogenic fungus with wide host range as well as it is one of ubiquitous fungus worldwide (Chen et al., 2015; Montel et al., 1991). C. cladosporioides has been studied intensively because the fungus is an endophyte in Taxus brevifolia, taxol producing tree, a famous anti-cancer drug (Miao et al., 2018). C. cladosporioides is an endophytic fungus not only tree but also marine alga, and an recent discovering reported that the fungus biogenic synthesized even gold nanoparticles (Hulikere et al., 2017). Qin et al. (2016) reported that C. cladosporioides can promote plant seed germination and growth. The fungus was detected at all part of American sweetgum plant tissues, including rhizosphere, phyllosphere, and endosphere. Torres et al. (2017) observed and verified that C. cladosporioides suppressed rust disease occurrence due to antagonistic and hyperparasitic effects of the fungus against the pathogen. Assignment of functional guild using FUNGuild revealed difference between rhizosphere and endosphere fungal population (Fig. 3, Supplementary Fig. 3). Abundance of saprophytic fungi was no difference between rhizosphere and endosphere. However, abundance of pathotroph mode fungi were significantly higher in the sap than the rhizosphere. The finding suggested that putative pathogenic fungi may actively invade the plant through root system from rhizosphere.
In common bacterial OTUs in rhizosphere and endosphere, two strains were abundant in the sap than the rhizosphere. Curtobacterium sp. and Sphingomonas sp. were existed as 0.04% and 0.29% in rhizosphere, 1.22% and 0.76% in sap, respectively (Supplementary Table 3). Curtobacterium sp. is Gram-positive corynebacteria, soil inhabitant, and isolated from plant tissues including endosphere (Lumactud et al., 2016). The bacterium is also known as a potential biocontrol agent for sweet orange (Garrido et al., 2016). Sphingomonas sp. LK11, endophytic strain in Tephrosia apollinea, contributed plant growth through producing gibberellins and indole acetic acid (Khan et al., 2014). The LK11 stain also significantly enhanced tomato growth. KEGG BRITE mapping results presented enzymatic activities of the bacterial communities in rhizosphere and endosphere of kiwifruit (Supplementary Fig. 4). Among 41 categories of enzyme activity have been assigned and endospheric community showed higher enzymes activities than rhizospheric microbiome community. This funding suggested that although bacterial diversity and population density were less in endosphere than rhizosphere, the endosphere bacterial community much more active than rhizosphere community and the activity may relate plant growth promotion effect or protection against both biotic and abiotic stresses.
The microbial communities of kiwifruit are probably shaped by a broad variation of environmental and host-microbe related factors (Coleman-Derr et al., 2015). The rhizosphere has relatively high microbial diversity, while sap has low diversity. The fact that rhizosphere microbiome is richer than endophytes microbiome was well established (Daniel, 2005). Our findings supplemented to these general aspects in both fungal and bacterial community compositions and relationships between migration and function in the same kiwifruit tree. Phylum Ascomycota of fungi and Proteobacteria of bacteria were the most dominant groups in kiwifruit microbiome. Interestingly enough, putative pathogenic fungi and beneficial bacteria were mostly abundant groups in the shared OTUs between rhizosphere and endosphere. Bacteria and fungi competed with the effective endophytic rhizobia in the limited space and nutrients (Miao et al., 2017). The fungal endophytes have been recognized as two major groups clavicipitaceous endophytes (C-endophytes) class I and nonclavicipitaceous endophytes (NC-endophytes) class II, III, IV to differences in evolutionary relevance, taxonomy, host range and ecological functions (Rodriguez et al., 2008). Ascomycetous fungi have broad and diverse range in all the classes. Among Proteobacteria classes, α, β, and γ-Proteobacteria were reported as the most abundant bacteria in plant endophytic communities (Santoyo et al., 2016). The most commonly found endophyte genera are also rhizosphere inhabitants. Ascending migration of endophytes from roots to leaves had been demonstrated inside rice plants (Chi et al., 2005). Even though endophytic migration is adapted to specific plant genotypes, various endophytic communities are influenced by abiotic and biotic factors, such as environmental conditions, microbe-microbe interactions and plant-microbe interactions (Miliute et al., 2015). Through this kiwifruit microbiota study, results may provide a foundation for selection and development of kiwifruit growth promoting microbes or biological control agents.

Supplemental Materials


This research was supported by Next-Generation BioGreen 21 program (PJ013250) and by the 2019 R&D fund of the Animal and Plant Quarantine Agency.


Electronic Supplementary Material

Supplementary materials are available at The Plant Pathology Journal website (

Fig. 1
Principal coordinates analysis (PCoA) based on Bray-Curtis distance. Sap and rhizosphere samples (A) and sap samples (C) only at the operational taxonomic units (OTUs) level show theirs distance of fungal and bacterial composition similarity. Dots representing OTUs for sap and rhizosphere (B) and sap only sample (D) show that each OTU closer to the sample position in PCoA contributes more. The red, green, and blue OTU dots mean fungal genus Didymella, bacterial genus Pseudomonas and Streptomyces.
Fig. 2
Relative abundance at family level. Unitedness of bacterial abundance (A) and fungal abundance (B) are calculated to local contribution to beta diversity (LCBD).
Fig. 3
Predicted fungal trophic mode using FUNGuild. Bar heights mean average abundance (%). Their standard deviations are described by error bars.
Table 1
Each top 5 fungal genera in top 3 phyla component (%)
Taxa Combined Rhizosphere Sap
Ascomycota 63.29 ± 17.88 42.82 ± 7.76 79.75 ± 6.04
Didymella 8.96 ± 12.09 1.90 ± 0.83 16.02 ± 13.85
Cladosporium 7.04 ± 6.99 0.42 ± 0.19 13.65 ± 3.17
Paraphoma 6.05 ± 6.92 7.05 ± 9.63 5.05 ± 1.03
Not assigned family Microascaceae at genus level (OTU_F2) 5.85 ± 6.75 11.70 ± 4.78 0.00 ± 0.00
Aureobasidium 3.87 ± 7.32 0.02 ± 0.02 7.71 ± 8.81
Mortierellomycota 21.49 ± 18.50 39.13 ± 7.52 3.85 ± 2.44
Mortierella 21.44 ± 18.48 39.05 ± 7.52 3.82 ± 2.42
Basidiomycota 11.31 ± 7.66 15.29 ± 8.96 7.34 ± 2.34
Filobasidium 2.95 ± 6.21 0.03 ± 0.02 5.88 ± 7.75
Soilcocozyma 1.09 ± 0.52 1.49 ± 0.23 0.69 ± 0.41
Cystofilobasidium 0.74 ± 0.58 0.63 ± 0.52 0.85 ± 0.61
Cryptococcus 0.73 ± 1.00 0.00 ± 0.00 1.45 ± 0.97
Tausonia 0.53 ± 0.66 1.04 ± 0.63 0.04 ± 0.04

OTU, operational taxonomic unit.

Table 2
Each top 5 bacterial families in top 3 phyla component (%)
Taxa Combined Rhizosphere Sap
Proteobacteria 51.42 ± 28.37 27.60 ± 1.89 75.24 ± 21.71
Pseudomonadaceae 27.20 ± 38.18 1.20 ± 0.82 53.20 ± 39.55
Beijerinckiaceae 3.90 ± 6.04 0.63 ± 0.00 7.16 ± 7.19
Sphingomonadaceae 3.21 ± 3.10 2.43 ± 0.70 3.99 ± 4.19
Burkholderiaceae 2.92 ± 2.56 2.40 ± 0.70 3.43 ± 3.47
SC-I-84 1.58 ± 1.60 3.15 ± 0.47 0.01 ± 0.01
Actinobacteria 15.95 ± 9.78 16.73 ± 0.64 15.17 ± 13.77
Streptomycetaceae 2.92 ± 2.55 3.57 ± 0.53 0.18 ± 0.14
Geodermatophilaceae 1.60 ± 2.80 0.26 ± 0.02 2.95 ± 3.48
Nocardiaceae 1.53 ± 1.72 0.19 ± 0.04 2.86 ± 1.53
Microbacteriaceae 1.52 ± 1.90 0.41 ± 0.04 2.63 ± 2.19
Micromonosporaceae 0.90 ± 0.77 1.64 ± 0.25 0.15 ± 0.13
Firmicutes 13.09 ± 12.00 25.73 ± 0.57 1.15 ± 1.53
Bacillaceae 6.73 ± 6.46 13.18 ± 0.20 0.28 ± 0.36
Paenibacillaceae 3.36 ± 3.34 6.69 ± 0.30 0.03 ± 0.04
Planococcaceae 2.21 ± 2.16 4.33 ± 0.58 0.08 ± 0.12
Alicyclobacillaceae 0.30 ± 0.33 0.59 ± 0.24 0.01 ± 0.02
Thermoactinomycetaceae 0.21 ± 0.39 0.06 ± 0.02 0.36 ± 0.51


Callahan, BJ, McMurdie, PJ, Rosen, MJ, Han, AW, Johnson, AJA and Holmes, SP 2016. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 13:581-583.
crossref pmid pmc
Chen, Q, Jiang, JR, Zhang, GZ, Cai, L and Crous, PW 2015. Resolving the Phoma enigma. Stud Mycol. 82:137-217.
crossref pmid pmc
Chi, F, Shen, S-H, Cheng, H-P, Jing, Y-X, Yanni, YG and Dazzo, FB 2005. Ascending migration of endophytic rhizobia, from roots to leaves, inside rice plants and assessment of benefits to rice growth physiology. Appl Environ Microbiol. 71:7271-7278.
crossref pmid pmc
Cho, G, Kim, M-J, Kwon, Y and Kwak, Y-S 2018. Comparison of endophytic microbial community in kiwifruit plant cultivars. Plant Pathol J. 34:341-346.
crossref pmid pmc
Coleman-Derr, D, Desgarennes, D, Fonseca-Garcia, C, Gross, S, Clingenpeel, S, Woyke, T, North, G, Visel, A, Partida-Martinez, LP and Tringe, SG 2015. Plant compartment and biogeography affect microbiome composition in cultivated and native Agave species. New Phytol. 209:798-811.
crossref pmid pmc
Cregger, MA, Veach, AM, Yang, ZK, Crouch, MJ, Vilgalys, R, Tuskan, GA and Schadt, CW 2018. The Populus holobiont: dissecting the effects of plant niches and genotype on the microbiome. Microbiome. 6:31
crossref pmid pmc
Daniel, R 2005. The metagenomics of soil. Nat Rev Microbiol. 3:470-478.
crossref pmid
Garrido, LM, Alves, JMP, Oliveria, LS, Gruber, A, Padilla, G and Araújo, WL 2016. Draft genome sequence of Curtobacterium sp. strains ER1/6, and endophytic strain isolated from Citrus sinensis with potential to be used as a biocontrol agent. Genome Announc. 4:e01264-16.
pmid pmc
Hulikere, MM, Joshi, CG, Danagoudar, A, Poyya, J, Kudva, AK and Dhanajaya, BL 2017. Biogenic synthesis of gold nanoparticles by marine endophytic fungus-Cladosporium cladosporioides isolated from seaweed and evaluation of their antioxidant and antimicrobial properties. Process Biochem. 63:137-144.
Khan, AL, Waqas, M, Kang, S-M, Al-Harrasi, A, Hussain, J, Al-Rawahi, A, Al-Khiziri, S, Ullah, I, Ali, L, Jung, H-Y and Lee, I-J 2014. Bacterial endophyte Sphingomonas sp. LK11 produces gibberellins and IAA and promotes tomato plant growth. J Microbiol. 52:689-695.
crossref pmid
Kõljalg, U, Nilsson, RH, Abarenkov, K, Tedersoo, L, Taylor, AFS, Bahram, M, Bates, ST, Bruns, TD, Bengtsson-Palme, J, Callaghan, TM, Douglas, B, Drenkhan, T, Eberhardt, U, Dueñas, M, Grebenc, T, Griffith, GW, Hartmann, M, Kirk, PM, Kohout, P, Larsson, E, Lindahl, BD, Lücking, R, Martín, MP, Matheny, PB, Nguyen, NH, Niskanen, T, Oja, J, Peay, KG, Peintner, U, Peterson, M, Põldmaa, K, Saag, L, Saar, I, Schüßler, A, Scott, JA, Senés, C, Smith, ME, Suija, A, Taylor, DL, Telleria, MT, Weiss, M and Larsson, KH 2013. Towards a unified paradigm for sequence-based identification of fungi. Mol Ecol. 22:5271-5277.
crossref pmid
Langille, MGI, Zaneveld, J, Caporaso, JG, McDonald, D, Knights, D, Reyes, JA, Clemente, JC, Burkepile, DE, Vega Thurber, RL, Knight, R, Beiko, RG and Huttenhower, C 2013. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 31:814-821.
crossref pmid pmc
Love, MI, Huber, W and Anders, S 2014. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550
crossref pmid pmc
Lumactud, R, Shen, SY, Lau, M and Fulthorpe, R 2016. Bacterial endophytes isolated from plants in natural oil seep soils with chronic hydrocarbon contamination. Front Microbiol. 7:755
crossref pmid pmc
McMurdie, PJ and Holmes, S 2013. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 8:e61217
crossref pmid pmc
Miao, L-Y, Mo, X-C, Xi, X-Y, Zhou, L, De, G, Ke, Y-S, Liu, P, Song, F-J, Jin, W-W and Zhang, P 2018. Transcriptome analysis of a taxol-producing endophytic fungus Cladosporium cladosporioides MD2. AMB Expr. 8:21
Miao, Y-Y, Shi, S-L, Zhang, J-G and Mohamad, OA 2017. Migration, colonization and seedling growth of rhizobia with matrine treatment in alfalfa (Medicago sativa L.). Acta Agric Scand B Soil Plant Sci. 68:26-38.
Miliute, I, Buzaite, O, Baniulis, D and Stanys, V 2015. Bacterial endophytes in agricultural crops and their role in stress tolerance: a review. Zemdirbyste. 102:465-478.
Montel, E, Bridge, PD and Sutton, BC 1991. An integrated approach to Phoma systematics. Mycopathologia. 115:89-103.
crossref pmid
Nguyen, NH, Song, Z, Bates, ST, Branco, S, Tedersoo, L, Menke, J, Schilling, JS and Kennedy, PG 2016. FUNGuild: an open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 20:241-248.
Qin, Y, Pan, X and Yuan, Z 2016. Seed endophytic microbiota in a coastal plant and phytobeneficial properties of the fungus Cladosporium cladosporioides. Fungal Ecol. 24:53-60.
Ottesen, AR, Peña, AG, White, JR, Pettengill, JB, Li, C, Allard, S, Rideout, S, Allard, M, Hill, T, Evans, P, Strain, E, Musser, S, Knight, R and Brown, E 2013. Baseline survey of the anatomical microbial ecology of an important food plant: Solanum lycopersicum (tomato). BMC Microbiol. 13:114
crossref pmid pmc
Quast, C, Pruesse, E, Yilmaz, P, Gerken, J, Schweer, T, Yarza, P, Peplies, J and Glöckner, FO 2012. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41:D590-D596.
crossref pmid pmc
Quiza, L, St-Arnaud, M and Yergeau, E 2015. Harnessing phytomicrobiome signaling for rhizosphere microbiome engineering. Front Plant Sci. 6:507
crossref pmid pmc
Rodriguez, RJ, White, JF Jr, Arnold, AE and Redman, RS 2008. Fungal endophytes: diversity and functional roles. New Phytol. 182:314-330.
Santoyo, G, Moreno-Hagelsieb, G, del Carmen Orozco-Mosqueda, M and Glick, BR 2016. Plant growth-promoting bacterial endophytes. Microbiol Res. 183:92-99.
crossref pmid
Torres, DE, Rojas-Martínez, RI, Zavaleta-Mejía, E, Guevara-Fefer, P, Márquez-Guzmán, GJ and Pérez-Martínez, C 2017. Cladosporium cladosporioides and Cladosporium pseudocladosporioides as potential new fungal antagonists of Puccinia horiana Henn., the causal agent of chrysanthemum white rust. PLoS ONE. 12:e0170782
crossref pmid pmc
Vandenkoornhuyse, P, Quaiser, A, Duhamel, M, Le Van, A and Dufresne, A 2015. The importance of the microbiome of the plant holobiont. New Phytol. 206:1196-1206.
crossref pmid
Wickham, H 2016. ggplot2: elegant graphics for data analysis. 2nd ed. Springer International Publishing, Cham, Switzerland. 260.

Editorial Office
Rm,904 (New Bldg.) The Korean Science & Technology Center 22,
Teheran-ro 7-Gil, Gangnamgu, Seoul 06130, Korea
Tel: +82-2-557-9360    Fax: +82-2-557-9361    E-mail:                

Copyright © 2023 by Korean Society of Plant Pathology.

Developed in M2PI

Close layer
prev next