Apple endophytic microbiota of different rootstock/scion combinations suggests a genotype-specific influence
© The Author(s). 2018
Received: 1 October 2017
Accepted: 18 January 2018
Published: 27 January 2018
High-throughput amplicon sequencing spanning conserved portions of microbial genomes (16s rRNA and ITS) was used in the present study to describe the endophytic microbiota associated with three apple varieties, “Royal Gala,” “Golden Delicious,” and “Honey Crisp,” and two rootstocks, M.9 and M.M.111. The objectives were to (1) determine if the microbiota differs in different rootstocks and apple varieties and (2) determine if specific rootstock-scion combinations influence the microbiota composition of either component.
Results indicated that Ascomycota (47.8%), Zygomycota (31.1%), and Basidiomycota (11.6%) were the dominant fungal phyla across all samples. The majority of bacterial sequences were assigned to Proteobacteria (58.4%), Firmicutes (23.8%), Actinobacteria (7.7%), Bacteroidetes (2%), and Fusobacteria (0.4%). Rootstocks appeared to influence the microbiota of associated grafted scion, but the effect was not statistically significant. Pedigree also had an impact on the composition of the endophytic microbiota, where closely-related cultivars had a microbial community that was more similar to each other than it was to a scion cultivar that was more distantly-related by pedigree. The more vigorous rootstock (M.M.111) was observed to possess a greater number of growth-promoting bacterial taxa, relative to the dwarfing rootstock (M.9).
The mechanism by which an apple genotype, either rootstock or scion, has a determinant effect on the composition of a microbial community is not known. The similarity of the microbiota in samples with a similar pedigree suggests the possibility of some level of co-evolution or selection as proposed by the “holobiont” concept in which metaorganisms have co-evolved. Clearly, however, the present information is only suggestive, and a more comprehensive analysis is needed.
Apple (Malus x domestica) is considered one of the most important fruit crops in the world with an annual production exceeding 80 million tons (http://faostat.fao.org). Commercial production typically utilizes the use of a selected scion cultivar grafted on a specific rootstock. This horticultural practice has a wide variety of benefits in regard to size control and related management practices, disease resistance, and fruit quality traits (taste, size, firmness, storability, etc.). Malling 9 (M.9) is a commonly used rootstock. Although it has a good size control, which results in tree size of 45–50% of a standard tree size, it is susceptible to fire blight and wooly apple aphid, but fairly resistant to crown and other root rots. In contrast, Malling Merton 111 (M.M.111) rootstock has resistance to wooly apple aphid and is moderately resistant (tolerant) to fire blight, as well as crown and root rot. Nevertheless, it does not exert much size control, being classified as a semi-dwarf, with resulting grafted trees being about 85 to 100% the size of ungrafted trees. The various attributes that contribute to disease resistance and size control, however, have not always been clearly elucidated, though specific genetic markers for various important traits in apple have been established [1–4].
High-throughput amplicon sequencing has proven to be an excellent tool in studying the microbial diversity of several environments, including the soil , phyllosphere [6, 7], carposphere , and rhizosphere . In regard to apples, the most detailed studies have been conducted in relation to the rhizosphere and plant pathogens associated with apple replant disease [10–12]. The microbiota of apple flower and fruit has also been characterized. Shade et al.  reported the succession of changes in microbial community structure over the life of apple flowers and suggested that this information could be used as a basis for developing ecological approaches to disease management. Abdelfattah et al.  noted that different portions of an apple fruit (stem-end, peel, wound, and calyx-end), surveyed at the point-of-purchase (supermarket), possess a distinct microbiota and noted a predominance of yeast often associated with dermal disorders (rashes) in humans.
The relationship between the resident microbiota and their plant hosts is essential to understanding agricultural production systems. Bacterial endophyte communities have been shown to be impacted by management practices , and the soil microbiome has been reported to influence grapevine-associated microbiota [15, 16]. The stability and uniformity of the composition, as well as the determinants of plant-associated microbiota, will not be clarified in agricultural crops, however, until many additional studies are conducted. In the present study, we describe the endophytic microbiota associated with three apple varieties, “Royal Gala,” “Golden Delicious,” and “Honey Crisp,” and two rootstocks, “M.9” (dwarfing) and “M.M.111” (semi-dwarfing). The objectives of the study were to (1) determine if the microbiota differs in different rootstocks and apple varieties and (2) determine if specific rootstock-scion combinations influence the composition of the microbiota of either component. The compositions of both the fungal and bacterial communities were assessed.
High-throughput amplicon sequencing
Endophytic microbiota of shoots of apple rootstocks and scions
Members of the Ascomycota were the dominant fungal phylum across all samples, accounting for 47.8% of the total number of detected sequences (Additional file 2: Figure S1). This was followed by Zygomycota (31.1%), Basidiomycota (11.6%), Glomeromycota (3.3%), and unidentified fungi (5.8%). Rozellomycota and Chytridiomycota were also detected but at a very low frequency (0.1 and 0.3%, respectively). Ascomycota were largely identified as members of the classes Dothideomycetes (36.6%), Sordariomycetes (5%), and Eurotiomycetes (3.3%). The Zygomycota was almost exclusively composed of members of the order Entomophthorales of the class Incertae sedis (31.1%). The Basidiomycota were principally represented by members of the Agaricomycetes (3.7%) and Tremellomycetes (4.1%).
Bacterial OTUs were assigned to Proteobacteria (58.4%), Firmicutes (23.8%), Actinobacteria (7.7%), Bacteroidetes (2%), and Fusobacteria (0.4%). The Proteobacteria were mainly represented by Gammaproteobacteria (37.2%), Betaproteobacteria (13.3%), and Alphaproteobacteria (7.8%), and Firmicutes were represented only by the classes Bacilli (22.4%) and Clostridia (1.4%). Actinobacteria (7.7%) was the only bacterial class detected within the phylum Actinobacteria (7.7%). Whereas, Bacteroidetes (2%) was represented by Flavobacteriia (1.1%), Saprospirae (0.4%), Bacteroidia (0.2%), and Cytophagia (0.2%) (Additional file 2: Figure S1).
Differences in the composition of the microbiota in “M.9” vs. “M.M.111” rootstocks
Comparisons between scion cultivars
The comparison between the three scions using Student’s t test revealed that there were significant differences in the alpha diversity (based on the Shannon index) of the endophytic OTUs present in “Golden Delicious,” “Honey Crisp,” and “Royal Gala” (Additional file 4: Table S3). No significant differences were observed, however, between “Royal Gala” and “Golden Delicious,” except in the number of observed OTUs. The beta diversity of the endophytic fungal microbiota of “Royal Gala” vs. “Honey Crisp,” and “Golden Delicious” vs. “Honey Crisp” was also significantly different. No statistically significant differences were observed, however, between the cultivars in regard to the alpha or beta diversity of the bacterial OTUs (Additional file 4: Table S3).
Impact of different rootstocks on the microbiota of scion varieties
Other studies have reported mixed results in regard to the impact of genotype on the composition of the associated microbiota. Weinert et al.  reported that tuber-associated bacteria were only weakly affected by genotype. Leff et al.  in a study of 33 strains of sunflower (Helianthus annus) that differed in their degree of domestication found that plant-associated fungal communities in the rhizosphere were more strongly influenced by host genetic factors and plant breeding than bacterial communities. They also found that there was a minimal vertical transmission of fungi from seeds to adult plants. The findings of the present study are in agreement with the results of the latter study. Bálint et al.  reported that host genotype had a structuring effect on the composition of foliar fungal communities in balsam poplar (Populus balsamifera) trees planted in a common garden. They hypothesized that there was either a filtering mechanism that allowed trees to be selective in their recruitment of fungi from the environment or that the fungal organisms were present systemically, thus allowing them to be “passed on” during successive clonal propagations. Sun et al.  found that host genotype (species) explained 30.1% of the variability (as determined by redundancy analysis) in the endophytic composition of three forest species (Betula platyphylla, Quercus liaotungensis, and Ulmus macrocarpa) in a mixed temperate forest. Pérez-Izquierdo et al.  reported that tree genotype was pivotal in structuring fungal communities in Mediterranean pine forests. Collectively, these studies suggest a degree of co-evolution or co-inheritance between a genotype and its associated microbiota (holobiont concept); however, more detailed and comprehensive studies will need to be conducted in apple to support such an inference.
This is the first examination of the endophytic population of apple using high-throughput amplicon sequencing. Beneficial symbionts, such as Diversispora and Piloderma, are recognized as endophytes that have a beneficial influence on the growth and stress tolerance of plants [22, 23]. These two genera were more abundant in the “M.M.111” rootstock than in the “M.9,” where their presence was low or absent. It would be interesting to determine if the growth promoting taxa present in “M.M.111” and their absence in “M.9” plays a role in the growth-related attributes associated with these two different rootstocks. M.9 had a higher RA of genera such as Zoophthora, Acremonium, and species from the order Eurotiales that are either saprophytic or antagonists of insects. The high RA of Zoophthora was reflected in the dominant presence of Zygomycota in “Royal Gala” and “Golden Delicious.” While the dominant presence of the Zygomycota in these two cultivars contradicts the general belief that Ascomycota is the most dominant, plant-associated phylum, Ascomycota, was still the predominant phylum in the combined samples collectively.
The consistent high prevalence of Xanthomonadaceae in all rootstocks and scions was surprising. This family contains a long list of phytopathogenic bacteria including Xanthomonas species and Xylella fastidiosa which causes disease in a wide range of important agricultural crops. Paenibacillus has been detected in various environments but is commonly recognized as a predominantly endophytic bacterium in woody plants . It also functions in growth promotion and nitrogen fixation in plants . Interestingly, like the growth beneficial fungi, this growth-promoting bacterium was generally found to be more abundant in “M.M.111,” the vigorous rootstock. Members of the Propionibacterium had a similar abundance pattern as Paenibacillus. Propionibacteria are used in the production of vitamin B12 and vitamin B2 and the production of porphyrinogen and tetrapyrrole compounds and exopolysaccharides . As the name indicates, these bacteria are known to synthesize propionic acid, which can solubilize phosphate. Since more than 90% of phosphate present in the soil or absorbed by the plant is in an insoluble form, these bacterial taxa may play a role in making phosphorus available to the plant. Methylobacterium species are ubiquitous in the natural environment, both as free-living organisms in soil and water, but also on the phylloplane and in leaf, stem, and root tissues of plants. Some induce plant leaf and root nodule formation and can promote plant growth by the production of auxins . Some species of Burkholderiales are associated with plants and are known to have growth-promoting effects . Burkholderia species are common rhizosphere constituents of major agricultural crops including maize, rice, sugarcane, wheat, tomato, and potato. A valuable feature of some growth-promoting Burkholderia species is their capacity for biological nitrogen fixation, siderophore production, inorganic phosphates solubilization, indoleacetic acid production, and phytopathogen inhibition .
A greater diversity of fungi than bacteria was identified in the analysis of the microbiota of various apple rootstock/scion combinations. In regard to the fungi, a distinct influence of genotype (both rootstock and scion) on the composition of the fungal microbiota was observed as well as in the relative abundance of specific genera or species. The compositions of the fungal taxa in “Golden Delicious” and “Royal Gala,” two scions closely-related by pedigree, were determined to be more similar to each other than they were to the composition of the fungal taxa in “Honey Crisp,” a scion more distantly-related by pedigree. A similar relationship in bacterial taxa could not be established. Interestingly, “M.M.111” rootstock, which does not have strong growth control properties, had a great number of beneficial and growth promoting fungal and bacterial taxa than “M.9” rootstock which does exhibit significant size control. Also of interest was the presence of high relative abundance of the genus Zoophthora in “Golden Delicious” and “Royal Gala,” a fungal genus generally recognized as an insect pathogen. Evidence for a genotype-specific impact on the composition of the endophytic microbiome in apple shoots suggests that further studies should be conducted to determine if the “holobiont” concept is supported in apple. Only a few microbiome-related studies have been conducted on temperate fruit trees, and it is expected that future research will have important implications on many aspects of tree physiology and cultivar development.
High-throughput amplicon sequencing of the endophytic microbiota was conducted on two-year-old, clonally propagated and grafted apple trees grown in pots in a greenhouse located at the US Department of Agriculture-Agricultural Research Service (USDA-ARS), Appalachian Fruit Research Station, Kearneysville, WV, USA. The trees were part of a longer-term study on tree growth and root architecture and were maintained at ambient light conditions and temperatures, watered, and fertilized on a regular basis. Trees of self-rooted “M.9” and “M.M.111” rootstocks and three apple scion varieties (“Honey Crisp,” “Golden Delicious,” and “Royal Gala”) grafted on either rootstock were used in the study. The microbial analysis was conducted on a total of 24 stem samples (three varieties on two rootstocks and two rootstocks × three biological replicates), destructively sampled, frozen, and ground in liquid nitrogen, prior to DNA extraction.
Several shoots (20 cm in length) were collected from each sampled tree, surface sterilized with 5% sodium hypochlorite (v:v), and rinsed three times in sterile water. The bark of each shoot was then removed with a sterile razor. The remaining xylem tissues were then immersed in liquid nitrogen and maintained at − 20 °C until lyophilization. Lyophilized samples were homogenized in a 2010 Geno/Grinder (SPEX SamplePrep, Metuchen, NJ, USA) using autoclaved metal beads. Total DNA was extracted from the ground xylem tissues using a Wizard Genomic DNA Purification Kit according to the manufacturer’s protocol (Promega, Madison, WI, USA). The quantity of each of the DNA samples was determined using a spectrophotometer (Nanodrop; Thermo Fisher Scientific Inc.), and the total DNA concentration was adjusted to 5.0 ng μL−1. The fungal ITS2 region was amplified using the universal primers ITS3/KYO2 and ITS4 to amplify the ITS2 region of the ribosomal DNA . The bacterial 16S region was amplified using the protocol described by Lundberg et al. , including the use of a pair of peptide-nucleic-acids (PNA) that were incorporated into the PCR amplification in order to reduce the generation of non-target chloroplast and mitochondrial amplicons. The universal 16S primer pair 515F and 806R was used to generate bacterial-derived 16S amplicons . All primers were modified to include Illumina adaptors (www.illumina.com). PCR reactions were conducted in a total volume of 25 μL containing 12.5 μL of KAPA HiFi HotStart ReadyMix (Kapa Biosystems, Wilmington, MA, USA), 1.0 μL of each primer (10 μM), 2.5 μL of DNA template, and 8.0 μL nuclease-free water. The reactions were incubated in a T100 thermal cycler (Bio-Rad) for 3 min at 98 °C followed by 30 cycles of 30 s at 95 °C, 30 s at 50 °C, and 30 s at 72 °C. All reactions ended with a final extension of 1 min at 72 °C. Nuclease-free water (QIAGEN, Valencia, CA, USA) replaced template DNA in negative controls. All amplicons and amplification mixtures including negative controls were sequenced on a MiSeq platform using V2 chemistry (Illumina, San Diego, CA, USA).
Paired-end reads were merged using PEAR 0.9.6 Paired-End reAd mergeR  and default parameters. The CLC genomics workbench V8 (Qiagen) was used for primer and quality trimming with a minimum of Q20. Sequences without either primer were discarded. Chimeric sequences were identified and filtered using VSEARCH 1.4.0 . The UCLUST algorithm  of the software package QIIME 1.9.1  was used to cluster sequences at a similarity threshold of 97% against UNITE dynamic database  released on 31 January 2016 for ITS reads and Greengenes Database  for the 16S reads. Sequences that failed to cluster against the database were de novo clustered using the same algorithm. The most abundant sequences in each OTU were selected as representative sequences and used for the taxonomic assignment using the BLAST algorithm  as implemented in QIIME 1.9.1.
The OTU table was normalized by rarefaction to an even sequencing depth in order to remove sample heterogeneity. The rarefied OTU table was used to calculate alpha diversity indices including observed species (Sobs), Chao1, and Shannon metrics. MetagenomeSeq’s cumulative sum scaling (CSS) was used as a normalization method for other downstream analyses, including taxa relative abundance, β-diversity, and group significance . Alpha diversities were compared based on a two-sample t test using non-parametric (Monte Carlo) methods and 999 Monte Carlo permutations (999). The CSS-normalized OTU table was analyzed using Bray Curtis metrics  and utilized to evaluate the β-diversity and construct PCoA plots . Only OTUs that were found in at least 50% of the samples were used to calculate the differential OTU abundance using the t test and the Kruskal-Wallis test . In all tests, significance was determined using 999 Monte Carlo permutations, and the false discovery rate (FDR) was used to adjust the calculated P value, and when the FDR-P ≤ 0.05, it was considered significant. PC1 and PC2 scores of the PCoA results were subjected to ANOVA with the Fisher’s LSD post hoc test to determine the significance of the sample clustering using IBM SPSS Statistics version 22 (IBM Corp., Armonk, NY).
This study was supported by the National Key R&D Program of China (2017YFD0201100), a grant from the US-Israel Binational Agricultural Research and Development Fund (IS 2513-16), and Foundation for High-level Talents of Chongqing University of Arts and Sciences (R2016LX01).
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files.
MW conceived and designed the experiments. JL, AA, JN, and EB performed the experiments. JL, AA, LS, and SD analyzed the date. JL, AA, and MW were the major contributors in writing the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Chagné D, Dayatilake D, Diack R, Oliver M, Ireland H, Watson A, Gardiner SE, Johnston JW, Schaffer RJ, Tustin S. Genetic and environmental control of fruit maturation, dry matter and firmness in apple (Malus × domestica Borkh.). Hortic Res. 2014;1:14046.View ArticlePubMedPubMed CentralGoogle Scholar
- Costa F. MetaQTL analysis provides a compendium of genomic loci controlling fruit quality traits in apple. Tree Genet Genomes. 2015;11:819.View ArticleGoogle Scholar
- Muranty H, Troggio M, Sadok IB, Rifaï MA, Auwerkerken A, Banchi E, Velasco R, Stevanato P, van de Weg WE, Di Guardo M, Kumar S, Laurens F, Bink MC. Accuracy and responses of genomic selection on key traits in apple breeding. Hortic Res. 2015;2:15060.View ArticlePubMedPubMed CentralGoogle Scholar
- Duan N, Bai Y, Sun H, Wang N, Ma Y, Li M, Wang X, Jiao C, Legall N, Mao L, Wan S, Wang K, He T, Feng S, Zhang Z, Mao Z, Shen X, Chen X, Jiang Y, Wu S, Yin C, Ge S, Yang L, Jiang S, Xu H, Liu J, Wang D, Qu C, Wang Y, Zuo W, Xiang L, Liu C, Zhang D, Gao Y, Xu Y, Xu K, Chao T, Fazio G, Shu H, Zhong GY, Cheng L, Fei Z, Chen X. Genome re-sequencing reveals the history of apple and supports a two-stage model for fruit enlargement. Nat Commun. 2017;8:249.View ArticlePubMedPubMed CentralGoogle Scholar
- Buee M, Reich M, Murat C, Morin E, Nilsson RH, Uroz S, Martin F. 454 Pyrosequencing analyses of forest soils reveal an unexpectedly high fungal diversity. New Phytol. 2009;184:449–56.View ArticlePubMedGoogle Scholar
- Abdelfattah A, Wisniewski M, Li Destri Nicosia MG, Cacciola SO, Schena L. Metagenomic analysis of fungal diversity on strawberry plants and the effect of management practices on the fungal community structure of aerial organs. PLoS One. 2016;11:e0160470.View ArticlePubMedPubMed CentralGoogle Scholar
- Abdelfattah A, Cacciola SO, Mosca S, Zappia R, Schena L. Analysis of the fungal diversity in citrus leaves with greasy spot disease symptoms. Microb Ecol. 2017;73:739–49.View ArticlePubMedGoogle Scholar
- Abdelfattah A, Wisniewski M, Droby S, Schena L. Spatial and compositional variation in the fungal communities of organic and conventionally grown apple fruit at the consumer point-of-purchase. Hortic Res. 2016;3:6047.View ArticleGoogle Scholar
- Mendes R, Garbeva P, Raaijmakers JM. The rhizosphere microbiome: significance of plant beneficial, plant pathogenic, and human pathogenic microorganisms. FEMS Microbiol Rev. 2013;37:634–63.View ArticlePubMedGoogle Scholar
- Mazzola M, Cook RJ. Effects of fungal root pathogens on the population dynamics of biocontrol strains of fluorescent Pseudomonads in the wheat rhizosphere. Appl Environ Microbiol. 1991;57:2171–8.PubMedPubMed CentralGoogle Scholar
- Zhang Q, Sun J, Liu S, Wei Q. Manure refinement affects apple rhizosphere bacterial community structure: a study in sandy soil. PLoS One. 2013;10:e76937.View ArticleGoogle Scholar
- Mazzola M, Hewavitharana SS, Strauss SL. Brassica seed meal soil amendments transform the rhizoshpere microbiome and improve apple production through resistance to pathogen reinfestation. Phytopathology. 2015;105:460–9.View ArticlePubMedGoogle Scholar
- Shade A, McManus PS, Handelsman J. Unexpected diversity during community succession in the apple flower microbiome. MBio. 2013;4:e00602–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Wemheuer F, Kaiser K, Karlovsky P, Daniel R, Vidal S, Wemheuer B. Bacterial endophyte communities of three agricultural important grass species differ in their response towards management regimes. Sci Rep. 2017;7:40914.View ArticlePubMedPubMed CentralGoogle Scholar
- Campisano A, Antonielli L, Pancher M, Yousaf S, Pindo M, Pertot I. Bacterial endophytic communities in the grapevine depend on pest management. PLoS One. 2014;9:e112763.View ArticlePubMedPubMed CentralGoogle Scholar
- Zarraonaindia I, Owens SM, Weisenhorn P, West K, Hampton-Marcell J, Lax S, Bokulich NA, Mills DA, Martin G, Taghavi S, van der Lelie D, Gilbert JA. The soil microbiome influences grapevine-associated microbiota. MBio. 2015;6:e02527–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Weinert N, Meincke R, Gottwald C, Heuer H, Schloter M, Berg G, Smalla K. Bacterial diversity on the surface of potato tubers in soil and the influence of the plant genotype. FEMS Microbiol Ecol. 2010;74:114–23.View ArticlePubMedGoogle Scholar
- Leff JW, Lynch RC, Kane NC, Fierer N. Plant domestication and the assembly of bacterial and fungal communities associated with strains of the common sunflower, Helianthus annuus. New Phytol. 2017;214:412–23.View ArticlePubMedGoogle Scholar
- Bálint M, Tiffin P, Hallström B, O'Hara RB, Olson MS, Fankhauser JD, Piepenbring M, Schmitt I. Host genotype shapes the foliar fungal microbiome of balsam poplar (Populus balsamifera). PLoS One. 2013;8:e53987.View ArticlePubMedPubMed CentralGoogle Scholar
- Sun X, Ding Q, Hyde KD, Guo LD. Community structure and preference of endophytic fungi of three woody plants in a mixed forest. Fungal Ecol. 2012;5:624–32.View ArticleGoogle Scholar
- Pérez-Izquierdo L, Zabal-Aguirre M, Flores-Rentería D, González-Martínez SC, Buée M, Rincón A. Functional outcomes of fungal community shifts driven by tree genotype and spatial-temporal factors in Mediterranean pine forests. Environ Microbiol. 2017;19:1639–52.View ArticlePubMedGoogle Scholar
- Mäkipää R, Rajala T, Schigel D, Rinne KT, Pennanen T, Abrego N, Ovaskainen O. Interactions between soil- and dead wood-inhabiting fungal communities during the decay of Norway spruce logs. ISME J. 2017;11:1964–74.View ArticlePubMedPubMed CentralGoogle Scholar
- Zou YN, Wang P, Liu CY, Ni QD, Zhang DJ, Wu QS. Mycorrhizal trifoliate orange has greater root adaptation of morphology and phytohormones in response to drought stress. Sci Rep. 2017;7:41134.View ArticlePubMedPubMed CentralGoogle Scholar
- Ulrich K, Stauber T, Ewald D. Paenibacillus—a predominant endophytic bacterium colonising tissue cultures of woody plants. Platn Cell Tiss Org. 2008;93:347–51.View ArticleGoogle Scholar
- Puri A, Padda KP, Chanway CP. Seedling growth promotion and nitrogen fixation by a bacterial endophyte Paenibacillus polymyxa P2b-2R and its GFP derivative in corn in a long-term trial. Symbiosis. 2016;69:123–9.View ArticleGoogle Scholar
- Kiatpapan P, Murooka Y. Genetic manipulation system in propionibacteria. J Biosci Bioeng. 2002;93:1–8.View ArticlePubMedGoogle Scholar
- Roat C, Saraf M. Unravelling the interaction of plant and their phyllosphere microbiome. In: Singh R, Kothari R, Koringa P, Singh S, editors. Understanding host-microbiome interactions—an omics approach. Singapore: Springer; 2017. p. 157–72.View ArticleGoogle Scholar
- Castanheira N, Dourado AC, Kruz S, Alves PI, Delgado-Rodríguez AI, Pais I, Semedo J, Scotti-Campos P, Sánchez C, Borges N, Carvalho G, Barreto Crespo MT, Fareleira P. Plant growth-promoting Burkholderia species isolated from annual ryegrass in Portuguese soils. J Appl Microbiol. 2016;120:724–39.View ArticlePubMedGoogle Scholar
- Paungfoo-Lonhienne C, Lonhienne TG, Yeoh YK, Donose BC, Webb RI, Parsons J, Liao W, Sagulenko E, Lakshmanan P, Hugenholtz P, Schmidt S, Ragan MA. Crosstalk between sugarcane and a plant-growth promoting Burkholderia species. Sci Rep. 2016;6:37389.View ArticlePubMedPubMed CentralGoogle Scholar
- Toju H, Tanabe AS, Yamamoto S, Sato H. High-coverage ITS primers for the DNA-based identification of Ascomycetes and Basidiomycetes in environmental samples. PLoS One. 2012;7:e40863.View ArticlePubMedPubMed CentralGoogle Scholar
- Lundberg DS, Yourstone S, Mieczkowski P, Jones CD, Dangl JL. Practical innovations for high-throughput amplicon sequencing. Nat Methods. 2013;10:999–1002.View ArticlePubMedGoogle Scholar
- Zhang J, Kobert K, Flouri T, Stamatakis A. PEAR: a fast and accurate Illumina paired-end reAdmergeR. Bioinformatics. 2014;30:614–20.View ArticlePubMedGoogle Scholar
- Rognes T, Flouri T, Nichols B, Quince C, Mahé F. VSEARCH: a versatile open source tool for metagenomics. PeerJ. 2016;4:e2584.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.View ArticlePubMedGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Abarenkov K, Nilsson HR, Larsson KH, Alexander IJ, Eberhardt U, Erland S, Høiland K, Kjøller R, Larsson E, Pennanen T, Sen R, Taylor AF, Tedersoo L, Ursing BM, Vrålstad T, Liimatainen K, Peintner U, Kõljalg U. The UNITE database for molecular identification of fungi—recent updates and future perspectives. New Phytol. 2010;186:281–5.View ArticlePubMedGoogle Scholar
- DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevi D, Hu P, Andersen GL. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72:5069–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.View ArticlePubMedGoogle Scholar
- Paulson JN, Stine OC, Bravo HC, Pop M. Differential abundance analysis for microbial marker-gene surveys. Nat Methods. 2013;10:1200–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Bray JR, Curtis JT. An ordination of the upland forest communities of Southern Wisconsin. Ecol Monogr. 1957;27:325–49.View ArticleGoogle Scholar
- Vázquez-Baeza Y, Pirrung M, Gonzalez A, Knight R. EMPeror: a tool for visualizing high-throughput microbial community data. Gigascience. 2013;2:16.View ArticlePubMedPubMed CentralGoogle Scholar
- Kruskal WH, Wallis WA. Use of ranks in one-criterion variance analysis. J Am Stat Assoc. 1952;47:583–621.View ArticleGoogle Scholar