Highly diversified fungi are associated with the achlorophyllous orchid Gastrodia flavilabella
- Tsunglin Liu†1,
- Ching-Min Li†2,
- Yue-Lun Han†1,
- Tzen-Yuh Chiang2,
- Yu-Chung Chiang3 and
- Huang-Mo Sung2Email author
© Liu et al.; licensee BioMed Central. 2015
Received: 27 August 2014
Accepted: 28 February 2015
Published: 14 March 2015
Mycoheterotrophic orchids are achlorophyllous plants that obtain carbon and nutrients from their mycorrhizal fungi. They often show strong preferential association with certain fungi and may obtain nutrients from surrounding photosynthetic plants through ectomycorrhizal fungi. Gastrodia is a large genus of mycoheterotrophic orchids in Asia, but Gastrodia species’ association with fungi has not been well studied. We asked two questions: (1) whether certain fungi were preferentially associated with G. flavilabella, which is an orchid in Taiwan and (2) whether fungal associations of G. flavilabella were affected by the composition of fungi in the environment.
Using next-generation sequencing, we studied the fungal communities in the tubers of Gastrodia flavilabella and the surrounding soil. We found (1) highly diversified fungi in the G. flavilabella tubers, (2) that Mycena species were the predominant fungi in the tubers but minor in the surrounding soil, and (3) the fungal communities in the G. flavilabella tubers were clearly distinct from those in the surrounding soil. We also found that the fungal composition in soil can change quickly with distance.
G. flavilabella was associated with many more fungi than previously thought. Among the fungi in the tuber of G. flavilabella, Mycena species were predominant, different from the previous finding that adult G. elata depends on Armillaria species for nutritional supply. Moreover, the preferential fungus association of G. flavilabella was not significantly influenced by the composition of fungi in the environment.
Mycorrhizal association between plants and fungi is a common phenomenon in plants. In fact, most orchid species depend on interactions with mycorrhizal fungi for completing their life cycle, particularly during their early developmental stages because orchids lack endosperm or seed-based nutrient reserves . Some orchid species are achlorophyllous through their entire life and must obtain carbon and nutrients for their growth and survival through mycorrhizal fungi . It has been suggested that fully mycoheterotrophic orchids often show high specificity toward their mycorrhizal fungi [3-5]. Conversely, mycorrhizal association may vary in photosynthetic orchids [6,7]. Thus, identifying mycorrhizal fungi that are important for the survival and growth of mycoheterotrophic orchids may provide insights into the evolutionary dynamics between orchids and their fungal associates.
Gastrodia species form one of the largest achlorophyllous and mycoheterotrophic orchid genera, and are distributed throughout Oceania to Asia, including Australia, Vietnam, China, Taiwan, Japan and South Korea . The rhizomes of Gastrodia elata are a prominent herbal medicine for human diseases such as vertigo, blackout and headache . It has been suggested that G. elata depends on Mycena fungi for germination [10-12] and relies on Armillaria fungi for carbon and nutritional supply after germination [13,14]. In addition, some Gastrodia species are associated with litter- or wood-decaying fungi [4,15,16]. Ogura-Tsujita et al. first demonstrated that the adult G. confusa gain carbon through several wood- or litter-decaying Mycena species . In contrast, G. similis was found to be associated largely with the saprotrophic Resinicium fungi . Moreover, G. sesamoides, a common obligate mycoheterotrophic orchid species, relies on the saprotrophic fungi Campanella and Marasmius for its carbon supply instead of on an ectomycorrhizal fungal partner of a photosynthetic plant . The above four Gastrodia species all depend on saprotrophic fungi for carbon and nutrition regardless of their geographic regions. In this study, we examined the fungal composition in the tubers of G. flavilabella, an endemic species of Taiwan , to identify the predominant fungi species associated with this plant. The identification is important because Gastrodia species do not always associate with the same fungal species for their growth.
Concurrent association of several fungi with an individual orchid is common, especially in photosynthetic orchids [17-20], suggesting the importance of studying the whole community of mycorrhizal fungi rather than the presence of individual fungal species. However, current mycorrhizal fungi identification requires microscopic identification and culturing or DNA segment cloning/sequencing, which may not provide sufficient resolution for the whole fungal community in the plant [14,16,17]. In this study, we used a deep sequencing approach to investigate the fungal community structure in the tubers of G. flavilabella and in the surrounding soil. We also investigated if G. flavilabella is preferentially associated with certain fungi.
Environmental and ecological factors may affect the extent of mycorrhizal preference of plant and the degree of dependence of plant on fungi for carbon and nutrients [7,21]. For example, the habitat characteristics and the presence of mycorrhizal fungi are important factors for the distribution of the orchid Cypripedium californicum because C. californicum is associated with multiple fungal families and may switch among different mycorrhizal fungi relatively easily . In addition, McCormick et al. showed that the distribution of a mycoheterotrophic orchid, Corallorhiza odontorhiza, and three green orchids correlated strongly with the abundance of the required mycorrhizal fungi . Thus, we also asked if the preferential fungal association of G. flavilabella is influenced by the fungal composition of the environment and/or by other environmental factors.
Fungal samples and the 28S rDNA sequences
Information of the Gastrodia flavilabella (Fla) tuber and soil samples used in this study.
No. of 28S raw PE reads in tuber; soil
Length of read1; 2 (bp)
No. and % of high-quality merged reads in tuber; soil
567041 (61.1%); 598310 (63.4%)
494672 (40.4%); 294985 (27.2%)
432956 (60.5%); 473080 (60.9%)
359376 (58.6%); 416882 (66.5%)
412503 (57.9%); 487634 (61.9%)
Taxonomic classification of the amplicon data
Statistics of the tuber and soil PE read data from different taxonomy domains
read number (%)
read number (%)
read number (%)
read number (%)
Comparison of fungal communities by amplicon size distribution
Highly diversified fungal communities in tubers and the surrounding soil
Number of OTUs in the reads sampled from the tuber and soil data
No. of OTUs in 100,000 reads sampled from tuber data
No. of OTUs in 100,000 reads sampled from soil data
210.1 ± 4.4
945.6 ± 9.1
129.2 ± 3.8
990.5 ± 7.5
353.7 ± 1.6
1017.9 ± 5.2
88.8 ± 3.3
1191.8 ± 4.5
82.8 ± 3.5
1018.4 ± 8.7
Fungal communities in the tubers and the soil are different
Fungal communities and compositions in the tubers and the soil
At all taxonomic levels, especially the lower ones, the fungal composition between a tuber sample and the corresponding soil sample was apparently different (Figure 5). This and the following observation were consistent with Figure 4. First, the fungal communities of the Fla2, Fla4, and Fla5 tuber samples were more similar to each other than to the Fla1 and Fla3 tuber samples. Second, among the soil samples, Fla3 and Fla4 showed more similar fungal compositions compared to other samples. Third, the fungal communities of soil were more complex than those observed in the tubers. Compared to the tuber samples, the higher percentages of unclassified amplicons in the soil samples might indicate a large fraction of novel fungal species in the soil.
Mycena species are the predominant fungi in G. flavilabella
The species closely related to Mycena cf. quiniaultensis OSC 67121 constituted 84.6-97.0% of the fungi communities in the Fla2, Fla4, and Fla5 tuber samples (Figure 5). In the Fla1 tuber sample, 47.4% of fungi communities were Mycena species, of which 27.1% and 20.1% were closely related to Mycena cf. quiniaultensis OSC 67121 and Mycena chlorophos, respectively. Thus, Mycena species were the predominant fungi in the G. flavilabella, except that only 8.1% of Mycena species were in the Fla3 tuber sample (Figure 5). Note that the taxonomy of some abundant OTUs did not appear in Figure 5 when they were unclassified at the taxonomic level (Additional file 1: Table S1). Take Fla3 tuber sample for example, the most abundant OTU (48.9% of fungal reads) was classified as the Agaricales; however, it was perfectly aligned to 13 different Mycena species and 7 non-Mycena species. That is, the most abundant OTU in the Fla3 tuber might still be Mycena species, but they were surely distinct from Mycena cf. quiniaultensis OSC 67121. Our data also suggested that the predominance of Mycena fungi in the G. flavilabella tubers was not influenced by the fungal composition of the environment because the percentages (0.3-2.5%) of Mycena fungi in all soil samples were much lower than those in the tuber samples (Figure 5).
Preferential association between Mycena fungi and G. flavilabella
Based on the clear distinction between the fungal communities in the G. flavilabella tubers and those in the surrounding soil (Figure 2, 3, and 4), we proposed that specific fungal communities had developed in the tubers of G. flavilabella. Moreover, the species closely related to Mycena cf. quiniaultensis OSC 67121 (identity ≥ 98%) was the predominant fungus in G. flavilabella tubers because it constituted 84.6-97.0% of the fungi communities in the Fla2, Fla4, and Fla5 tuber samples, but only a small fraction in all the soil samples (Figure 4). Interestingly, some differences were observed among the tubers. For example, Mycena species closely related to Mycena cf. quiniaultensis OSC 67121 or Mycena chlorophos were the top abundant fungi in the Fla1 tuber sample, while other Mycena and non-Mycena species were abundant in the Fla3 tuber sample. However, since they are all Mycena fungi, we concluded that adult G. flavilabella has a strong preference for association with Mycena fungi.
The highly diversified fungi in the tubers of G. flavilabella suggested that G. flavilabella does not have a strong defense mechanism to prevent fungal invasion and this perhaps reflects a fully mycoheterotrophic orchid’s nature to exploit any possible nutrient source. But why did the preference of G. flavilabella for Mycena fungi occur? It is possible that Mycena fungi were the first few fungi that invaded the tubers of G. flavilabella and G. flavilabella had since became dependent on Mycena fungi. It is also possible that the preference developed simply because Mycena fungi are saprotrophic and can provide G. flavilabella suitable nutrients. On the other hand, G. flavilabella might have changed its physiology to provide a stable environment for Mycena species or G. flavilabella might be able to stimulate the growth of Mycena species as has been observed in some mycorrhizal fungi .
Association of saprotrophic fungi with G. flavilabella
To date, achlorophyllous orchids were the only plants reported to solely rely on the saprotrophic fungi for carbon and nutrition [4,16]. Most of the mycoheterotrophic plant associated mycorrhizal fungi also form symbiotic relationship with photosynthetic plants, so that many mycoheterotrophic plants may also obtain carbon and nutrients from surrounding photosynthetic plants through ectomycorrhizal fungi . Ogura-Tsujita et al. reported that G. confusa is associated with several Mycena species and obtains carbon from these saprotrophic fungi but not from the surrounding photosynthetic plants . It has been suggested that adult G. elata depends on Armillaria species for nutritional supply and can be cultivated without photosynthetic plants [13,14]. Now, there is more evidence that Gastrodia species is associated with litter- or wood-decaying fungi such as Resinicium, Campanella and Marasmius species for carbon and nutrients [15,16]. We found that G. flavilabella has preferential association with Mycena fungi, and possibly also with Dictyopanus, Mycenoporella, Favolaschia, Panellus, Cruentomycena, and Resinomycena fungi. Association of these litter- or wood-decaying saprotrophic fungi suggested that G. flavilabella is a fully mycoheterotrophic orchid that does not rely on photosynthetic plants. However, we cannot rule out that G. flavilabella also obtains nutrients from surrounding photosynthetic plants.
Environmental factors might affect the fungal communities of tubers and the surrounding soil
The environmental factors of soil
Soil moisture a (%)
Soil organic matter b (%)
Soil total- nitrogen content c (%)
51.9 ± 0.3
35.0 ± 0.9
1.32 ± 0.01
32.4 ± 0.2
12.6 ± 0.4
0.45 ± 0.01
57.3 ± 0.3
40.7 ± 0.5
1.35 ± 0.04
36.6 ± 0
13.9 ± 0.2
0.46 ± 0.02
36.3 ± 0.3
14.0 ± 0.3
0.47 ± 0.01
Metazoa in the orchid samples
Metazoa can infect plants with no mutual benefits, leading to plant parasitism. For example, the nematode Meloidogyne incognita can invade the roots of almost all cultivated plants [26,27] and induce the re-differentiation of root cells . Many herbivore insects are also plant parasites and some arthropods can live within plants; e.g., phylloxera is an endo-parasite of grapes . It is thus not peculiar to observe metazoan sequences in our tuber samples. Among the metazoan sequences in our Fla tuber samples, nematodes and arthropods were the two major phyla in most cases (Additional file 1: Figure S4b).
The metazoan sequences might be in part due to contamination from the surrounding soil. In this scenario, we expected the metazoan community in the Fla tubers to be, to some extent, similar to that in the surrounding soil. However, the Fla1, Fla4, and Fla5 tubers did not cluster with the surrounding soil in UniFrac analysis (Additional file 1 Figure S5a). When focusing on the identities of metazoa, i.e., unweighted analysis, only the Fla3 tuber and soil samples clustered together (Additional file 1: Figure S5b). Thus, contamination should not be a major concern. Interestingly, we did not observe a common pattern of metazoan compositions in the Fla tubers, indicating no host specificity for metazoans.
LSU sequences have been collected in the databases RDP  and SILVA . However, we decided to curate our own LSU sequences because many LSU sequences in NCBI were not included in these databases. The RDP database (release 11) provides 62,860 fungal specific 28S rDNA sequences, so that they can not be used to classify non-fungi amplicons. In contrast, SILVA contains many non-fungal sequences but relatively few fungal sequences. Among the 39,412 LSU sequences in SILVA (release 115), 1,959, 5,600, and 4,145 belong to fungi, metazoa and viridiplantae, respectively. In contrast, we collected 268,614 LSU sequences, of which 120,617, 95,458, and 14,540 were fungi, metazoan and viridiplantae, respectively. Our more comprehensive LSU references allowed us to infer taxonomy with a higher accuracy. Note that some of our LSU sequences were not classified at all the seven taxonomic levels and some were only partial 28S rDNA sequences. This is perhaps one reason why other databases do not collect such sequences. Nevertheless, incomplete information was better than missing information and the partial 28S rDNA sequences were useful for our analysis.
Several computational tools, e.g., RDP classifier , Greengenes classifier  based on NAST alignments  and the BLAST-based method , have been proposed for taxonomic classification. The BLAST-based method was shown comparable to RDP classifier and Greengenes classifier in terms of internal consistency of taxonomic assignments, but was less sensitive than Greengenes . At the time of this study, NAST was not available. Moreover, Greengenes was only available as a web-server and thus could not accommodate our large datasets. We had tried the RDP classifier, but found some practical challenges. First, the RDP classifier requires references to have complete taxonomy classifications from domain to genus. This immediately dismissed those perfectly aligned or highly identical references that had missing taxonomy information. Second, the reference sequences of RDP classifier were not comparable to our amplicons in length. Our BLAST-based nearest neighbor approach classified most (92.9-98.4%) of the good-quality merged tuber amplicons to fungi, metazoa or viridiplantae (Table 2). Although 32.9-56.5% of the soil amplicons could not be classified to fungi, metazoa or viridiplantae, those unclassified amplicons in the soil samples might correspond to the novel species in the soil.
Confidence of taxonomy classification
Our own LSU references were still not sufficient for giving all the sequences the right classification. For example, most of the viridiplantae amplicons in all tuber samples were classified as “Viridiplantae; Streptophyta; Liliopsida; Asparagales; Orchidaceae; Campylocentrum; Campylocentrum micranthum” instead of being G. flavilabella. This was the result of no Gastrodia LSU sequence in the NCBI database. The classification, however, was reasonable because Gastrodia belongs to the Orchidaceae family. Moreover, the alignment identity of the corresponding OTU representative to the Campylocentrum micranthum reference was only 92.4% (data not shown).
The OTU representative classified as Mycena cf. quiniaultensis OSC 67121 was 99.4% identical to the corresponding reference sequence (Additional file 1: Table S1). A similar observation was made for the OTU classified as Mycena chlorophos, the third most abundant OTU in the Fla1 tuber sample. In the Fla3 tuber sample, the three most abundant OTUs were classified as Agaricales, Agaricomycetes and Ascomycota fungi, respectively, but were aligned to multiple taxonomic groups at the family or higher level. For example, the most abundant OTU was classified as Agaricales fungi; however, the OTU representative was perfectly aligned to 13 different Mycena species and 7 non-Mycena species. Moreover, none of the best alignments were to the Mycena cf. quiniaultensis OSC 67121, indicating that the dominating Mycena species in the Fla3 tuber sample was different from those in other Fla tuber samples. The taxonomic classification for the second and third most abundant OTUs were Agaricomycetes and Ascomycota fungi; however, none of the best alignments was to the Tricholomataceae family, indicating non-Mycena species of these two OTUs. We emphasized again that our BLAST-based nearest neighbor approach to species classifications might not be 100 percent correct if the species sequence information in the database was missing or not correct. However, it is certain that the sequences classified differently, e.g., as Mycena cf. quiniaultensis OSC 67121 and Mycena chlorophos, were distinct, thus truly representing different species.
This is the first study using a NGS deep sequencing approach for identifying the fungal communities in the tubers of G. flavilabella and in the surrounding soil. We found highly diversified fungal communities in the tubers, and Mycena species were the predominant fungi in the tubers of G. flavilabella, indicating that the mycoheterotrophic G. flavilabella has a preferential association with saprotrophic Mycena fungi. So far, it has been shown that Gastrodia species such as G. elata, G. confuse, G. similis and G. sesamoides are associated with litter- or wood-decaying fungi such as Armillaria, Mycena, Resinicium, Campanella and Marasmius species for carbon and nutrients [4,13-16]. Finally, the environmental factors such as soil water and nutrient contents might have a greater impact on the fungal communities in the tubers than the fungal communities in the soil.
In this study, we collected five tuber samples (Fla1-5) and five soil samples from their surrounding soil in the Hsitou area of Taiwan (Figure 1). Briefly, these samples were collected from an area with a 75 meter radius. Fla4 and Fla5 were located only about 0.05 meters apart from each other. Fla3 was ~150 meters away from Fla4 and Fla5. Fla1 and Fla2 were located between Fla3 and Fla4/Fla5 and the distance between Fla1 and Fla2 was ~12 meters. We individually collected one tuber of G. flavilabella and one surrounding soil sample (3 cm depth; 3 cm diameter) at each of five localities. The orchid tubers were surface-sterilized with 1% sodium hypochlorite for 30s and rinsed three times in sterile water for 30s to avoid contamination. The sterilized samples were then sagittal sectioned to about 1 g per fragment. Each plant surrounding soil was sieved to remove rocks and roots, and to break up the large soil aggregates. All of the orchid tuber and the surrounding soil samples were kept at −80°C for further analysis.
Genomic DNA extraction
The CTAB extraction method [35,36] was used to extract genomic DNAs from the plant tuber samples (~1 g). Genomic DNAs from the surrounding soil samples were extracted by a commercial DNA extraction kit (PowerSoil DNA Isolation Kit, MoBio).
PCR amplification of 28S rDNA markers
For each tuber or soil sample, the PCR reaction was carried out in a 50 μl reaction with 1ul template DNA, 0.5 unit of GoTaq polymerase (Promega, US), 5 μl of 10X PCR buffer, 5 μl of 1X dNTP mix (2.5 μM each dNTP), 3 μl of MgCl2 (25 mM), and 4 μl of each fungal-specific primer (3 μM). The PCRs were run with 35 cycles of 95°C for 30s, 56.2°C for 30s, and 72°C for 30s, and a final extension of 72°C for 10 min and the PCR products were purified with QiAquick PCR purification kit (Qiagen) following manufacture’s instruction. We designed a set of fungal specific LSU rDNA primers, AACACGGACCAAGGAGTC (forward) and CAGGCATAGTTCACCATCTT (reverse) which target a LSU region that is conserved in fungi but not in other organisms and amplify a variable region of the LSU rDNA. The amplified LSU regions are expected to be 167–218 bp in size, which can be covered by the Illumina paired-end reads.
Illumina paired-end (PE) sequencing
DNA sequencing was carried out on Illumina HiSeq 2500 platform (the Fla2 sample was on HiSeq 2000) at Yourgene Biosciences, Taiwan, following the manufacturer’s protocols. Briefly, amplicon DNA was A-tailed using the polymerase activity of Klenow fragment. Indexed adapters were then ligated to the DNA fragments by DNA ligase followed by PCR reaction of 10 to 18 cycles to enrich the adapter-modified DNA fragments. Before sequencing, the libraries were validated by QPCR, Expersion and Qubit.
For each Illumina PE library, we first aligned the 28S primer sequences to all reads using BLAST (v2.2.29+, options: −word_size 5 -evalue 0.001). A read was considered qualified if a whole primer was aligned on the positive strand of the read and the primer was aligned only once to the read. It was possible that both forward and reverse primers appeared on a read. In that case, the mate-read was also required contained both primers if the mate-read was long enough. For each qualified read, we trimmed the segment outside primer from the read if there was any. A PE was qualified if both reads were qualified and the corresponding primers formed a pair. The paired reads of all qualified PE were further merged by FLASH  (v1.2.7) with default parameters (requiring ≥10 bp overlap) and short merged reads (<100 bp) were discarded. High-quality merged reads were then selected for analysis. Specifically, we scanned each merged read using a window size of 5 bp and required the average quality to be at least 30 throughout the read. For OTU analysis by UPARSE, read orientation was adjusted so that the forward primer was at the 5′ end of all merged reads.
Large subunit rDNA references
We obtained from NCBI  non-redundant nucleotides (nt, last modified on 2014.03.11) and parsed out the large subunit (LSU) rDNA sequences. Specifically, we kept the sequences whose description ($d) contained the keywords “large”, “LSU”, “28S”, and “rRNA” (perl script: ($d = ~/[Ll]arge/ || $d = ~/[(nr)\s\(\-)LSU[\s\],]/ || $d = ~/[\s\(\/\-)2[3-8]S[\s\-]/] && ($d = ~/r[DR]NA/ || ($d = ~/ribosomal/ && $d = ~/[DR]NA/)). The sequences whose description contained the keyword “mRNA” or “spacer” were further excluded. The taxonomies at seven levels (domain, phylum, class, order, family, genus, and species) of the retained LSU rDNA sequences were then assign based on the gene ID. The mappings between gene ID and taxonomy ID were obtained from the file “gi_taxid_nucl.dmp” in the NCBI Taxonomy database (last modified on 2014.03.10). We parsed the two files “name.dmp” and “nodes.dmp” in the same database for the full lineages of all taxonomy IDs. Note that the parsed LSU rDNAs came from several domains, including fungi, metazoa, viridiplantae, etc. Also note that the taxonomy classifications of many sequences were null at some taxonomy levels, which were considered as unclassified at those levels.
To determine the taxonomy of the high-quality merged reads, we aligned the reads to the LSU rRNA references using MegaBLAST [39,40] (v2.2.29+, option: −word_size 16 -evalue 1e-10). For each read, the top alignment(s) with the lowest E-value, i.e., a nearest neighbor approach, was used for annotation. The corresponding taxonomy lineages, which could be more than one, were modified as follows. First, the lineages containing the keyword “uncultured” or “fungal_sp” were treated as unclassified lineages. Unclassified lineages were discarded if classified one(s) existed. Second, the lineages with the most complete classifications were selected. Third, the classification at a taxonomy level was considered unclassified if there was more than one classification. Based on the taxonomy, we split the reads by domain.
For each sample, we used the UPARSE pipeline  (v7.0.1090) to cluster the fungal reads into operational taxonomic units (OTUs) as follows. First, we took only the unique reads and recorded the read counts (usearch –derep_fulllength –sizeout). Second, the unique reads were sorted by read count and the singletons were removed (−sortbysize –minsize 2). Third, OTU clustering was performed (−cluster_otus). Fourth, chimeric OTUs were further filtered using NCBI LSU sequences as references (−uchime_ref –db ncbilsu.fa –strand plus –nonchimeras). Lastly, the unique reads were re-assigned to the OTUs with an identity cutoff 0.97 (−usearch_global –starnd plus –id 0.97) [42,43]. The statistics of OTUs were then obtained and the OTU representative sequences were collected for comparing fungal communities.
Comparison of fungal communities
To compare fungal communities, the OTU representatives of all tuber and soil samples were first collected for constructing a phylogenetic tree using the mothur package  (v1.32, commands: align.seqs, dist.seqs, clearcut). Mothur requires aligned reference sequences for tree construction. To prepare the references, we extracted the fungal LSU rRNA references and ran e-PCR  (v2.3.11, fahash option: −w 3, re-PCR option: −n 2 -g 2 1–300) to locate binding sites for our 28S primers. The fungal LSU rDNAs containing the binding sites of both primers were kept and the corresponding segments between the two binding sites were extracted. We then used Clustal Omega  (v1.2, default parameters) to obtain the multiple sequence alignment of the extracted segments, which served as the curated references for tree construction. Based on the tree, the distances between samples were calculated by Fast UniFrac  (v1.5.3). We did both weighted (i.e., considering the numbers of reads represented by OTU representatives) and un-weighted principle component and clustering analyses. The distances between fungal communities were visualized using the R package  (v3.0.1).
To show the fungal compositions of a sample, we re-assigned taxonomy to the OTU representatives as above and calculated their fractions. Note that we used only the fungal LSU sequences as references. The results of all samples at different taxonomy levels were shown in stacked histograms.
Availability of supporting data
The merged PE reads of the ten tuber and soil samples were deposited in NCBI SRA database under the SRA ID SRP054374. Detail information of DNA sequences along with the supporting Tables and Figures were included in the Additional file 1.
This study was supported by NSC grant (NSC 100- 2621- B-006- 004- MY3) to HMS.
- Rasmussen HN. Terrestrial Orchids, from Seed To Mycotrophic Plant. Cambridge: Cambridge University Press; 1995.View ArticleGoogle Scholar
- Leake JR. Myco-heterotroph/epiparasitic plant interactions with ectomycorrhizal and arbuscular mycorrhizal fungi. Curr Opin Plant Biol. 2004;7(4):422–8.View ArticlePubMedGoogle Scholar
- Barrett CF, Freudenstein JV, Taylor DL, Koljalg U. Rangewide analysis of fungal associations in the fully mycoheterotrophic Corallorhiza striata complex (Orchidaceae) reveals extreme specificity on ectomycorrhizal Tomentella (Thelephoraceae) across North America. Am J Bot. 2010;97(4):628–43.View ArticlePubMedGoogle Scholar
- Ogura-Tsujita Y, Gebauer G, Hashimoto T, Umata H, Yukawa T. Evidence for novel and specialized mycorrhizal parasitism: the orchid Gastrodia confusa gains carbon from saprotrophic Mycena. Proc Biol Sci. 2009;276(1657):761–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Ogura-Tsujita Y, Yukawa T. High mycorrhizal specificity in a widespread mycoheterotrophic plant, Eulophia zollingeri (Orchidaceae). Am J Bot. 2008;95(1):93–7.View ArticlePubMedGoogle Scholar
- Jacquemyn H, Honnay O, Cammue BP, Brys R, Lievens B. Low specificity and nested subset structure characterize mycorrhizal associations in five closely related species of the genus Orchis. Mol Ecol. 2010;19(18):4086–95.View ArticlePubMedGoogle Scholar
- Shefferson RP, Weiss M, Kull T, Taylor DL. High specificity generally characterizes mycorrhizal association in rare lady’s slipper orchids, genus Cypripedium. Mol Ecol. 2005;14(2):613–26.View ArticlePubMedGoogle Scholar
- Hsu T-C. Master, Taxonomy of Gastrodia (Orchidaceae) in Taiwan. Taipei: National Taiwan University; 2008.Google Scholar
- Xu J-T, Guo S-X. Retrospect on the research of the cultivation of Gastrodia elata Bl, a rare traditional Chinese medicine. Chin Med J (Engl). 2000;113(8):686–92.Google Scholar
- Fan L, Guo S-X. Interaction between protocorms of Gastrodia elata (orchidaceae) and Mycena dendrobii in symbiotic germination. Mycosystema. 1999;18(2):219–25.Google Scholar
- Fan L, Guo S-X, Xiao P-G. Interaction between protocorms of Gastrodia elata (orchidaceae) and Mycena anoectochila during symbiotic germination. Mycosystema. 2001;20(4):539–46.Google Scholar
- Kim YI, Chang KJ, Ka KH, Hur H, Hong IP, Shim JO, et al. Seed Germination of Gastrodia elata Using Symbiotic Fungi, Mycena osmundicola. Mycobiology. 2006;34(2):79–82.View ArticlePubMed CentralPubMedGoogle Scholar
- Cha JY, Igarashi T. Armillaria species associated with Gastrodia elata in Japan. Eur J Forest Pathol. 1995;25(6–7):319–26.View ArticleGoogle Scholar
- Sekizaki H, Kuninaga S, Yamamoto M, Asazu SN, Sawa S, Kojoma M, et al. Identification of Armillaria nabsnona in gastrodia tubers. Biol Pharm Bull. 2008;31(7):1410–4.View ArticlePubMedGoogle Scholar
- Dearnaley JDW, Bougoure JJ. Isotopic and molecular evidence for saprotrophic Marasmiaceae mycobionts in rhizomes of Gastrodia sesamoides. Fungal Ecol. 2010;3(4):288–94.View ArticleGoogle Scholar
- Martos F, Dulormne M, Pailler T, Bonfante P, Faccio A, Fournel J, et al. Independent recruitment of saprotrophic fungi as mycorrhizal partners by tropical achlorophyllous orchids. New Phytol. 2009;184(3):668–81.View ArticlePubMedGoogle Scholar
- Lievens B, van Kerckhove S, Juste A, Cammue BP, Honnay O, Jacquemyn H. From extensive clone libraries to comprehensive DNA arrays for the efficient and simultaneous detection and identification of orchid mycorrhizal fungi. J Microbiol Methods. 2010;80(1):76–85.View ArticlePubMedGoogle Scholar
- Martos F, Munoz F, Pailler T, Kottke I, Gonneau C, Selosse MA. The role of epiphytism in architecture and evolutionary constraint within mycorrhizal networks of tropical orchids. Mol Ecol. 2012;21(20):5098–109.View ArticlePubMedGoogle Scholar
- Roy M, Yagame T, Yamato M, Iwase K, Heinz C, Faccio A, et al. Ectomycorrhizal Inocybe species associate with the mycoheterotrophic orchid Epipogium aphyllum but not its asexual propagules. Ann Bot. 2009;104(3):595–610.View ArticlePubMed CentralPubMedGoogle Scholar
- Xing X, Ma X, Deng Z, Chen J, Wu F, Guo S. Specificity and preference of mycorrhizal associations in two species of the genus Dendrobium (Orchidaceae). Mycorrhiza. 2013;23(4):317–24.View ArticlePubMedGoogle Scholar
- Taylor DL, Bruns TD. Community structure of ectomycorrhizal fungi in a Pinus muricata forest: minimal overlap between the mature forest and resistant propagule communities. Mol Ecol. 1999;8(11):1837–50.View ArticlePubMedGoogle Scholar
- McCormick MK, Whigham DF, O’Neill JP, Becker JJ, Werner S, Rasmussen HN, et al. Abundance and distribution of Corallorhiza odontorhiza reflect variations in climate and ectomycorrhizae. Ecol Monogr. 2009;79(4):619–35.View ArticleGoogle Scholar
- Bidartondo MI. The evolutionary ecology of myco-heterotrophy. New Phytol. 2005;167(2):335–52.View ArticlePubMedGoogle Scholar
- Leake JR. Plants parasitic on fungi: unearthing the fungi in myco-heterotrophs and debunking the ‘saprophytic’ plant myth. Mycologist. 2005;19(3):113–22.Google Scholar
- Koorem K, Gazol A, Opik M, Moora M, Saks U, Uibopuu A, et al. Soil nutrient content influences the abundance of soil microbes but not plant biomass at the small-scale. PLoS One. 2014;9(3):e91998.View ArticlePubMed CentralPubMedGoogle Scholar
- Blaxter ML. Nematoda: genes, genomes and the evolution of parasitism. Adv Parasitol. 2003;54:101–95.View ArticlePubMedGoogle Scholar
- Trudgill DL, Blok VC. Apomictic, polyphagous root-knot nematodes: exceptionally successful and damaging biotrophic root pathogens. Annu Rev Phytopathol. 2001;39:53–77.View ArticlePubMedGoogle Scholar
- Caillaud MC, Lecomte P, Jammes F, Quentin M, Pagnotta S, Andrio E, et al. MAP65-3 microtubule-associated protein is essential for nematode-induced giant cell ontogenesis in Arabidopsis. Plant Cell. 2008;20(2):423–37.View ArticlePubMed CentralPubMedGoogle Scholar
- Granett J, Walker MA, Kocsis L, Omer AD. Biology and management of grape phylloxera. Annu Rev Entomol. 2001;46:387–412.View ArticlePubMedGoogle Scholar
- Cole JR, Wang Q, Fish JA, Chai B, McGarrell DM, Sun Y, et al. Ribosomal Database Project: data and tools for high throughput rRNA analysis. Nucleic Acids Res. 2014;42(Database issue):D633–42.View ArticlePubMed CentralPubMedGoogle Scholar
- Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41(Database issue):D590–6.View ArticlePubMed CentralPubMedGoogle Scholar
- DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72(7):5069–72.View ArticlePubMed CentralPubMedGoogle Scholar
- DeSantis Jr TZ, Hugenholtz P, Keller K, Brodie EL, Larsen N, Piceno YM, et al. a multiple sequence alignment server for comparative analysis of 16S rRNA genes. Nucleic Acids Res. 2006;34(Web Server issue):W394–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Liu Z, DeSantis TZ, Andersen GL, Knight R. Accurate taxonomy assignments from 16S rRNA sequences produced by highly parallel pyrosequencers. Nucleic Acids Res. 2008;36(18):e120.View ArticlePubMed CentralPubMedGoogle Scholar
- Gardes M, Bruns TD. ITS primers with enhanced specificity for basidiomycetes–application to the identification of mycorrhizae and rusts. Mol Ecol. 1993;2(2):113–8.View ArticlePubMedGoogle Scholar
- Murray MG, Thompson WF. Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 1980;8(19):4321–5.View ArticlePubMed CentralPubMedGoogle Scholar
- Magoc T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27(21):2957–63.View ArticlePubMed CentralPubMedGoogle Scholar
- Coordinators NR. Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2014;42(Database issue):D7–17.View ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.View ArticlePubMedGoogle Scholar
- Morgulis A, Coulouris G, Raytselis Y, Madden TL, Agarwala R, Schaffer AA. Database indexing for production MegaBLAST searches. Bioinformatics. 2008;24(16):1757–64.View ArticlePubMed CentralPubMedGoogle Scholar
- Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10(10):996–8.View ArticlePubMedGoogle Scholar
- Brown SP, Rigdon-Huss AR, Jumpponen A. Analyses of ITS and LSU gene regions provide congruent results on fungal community responses. Fungal Ecol. 2014;9:65–8.View ArticleGoogle Scholar
- Weber CF, Vilgalys R, Kuske CR. Changes in Fungal Community Composition in Response to Elevated Atmospheric CO2 and Nitrogen Fertilization Varies with Soil Horizon. Front Microbiol. 2013;4:78.View ArticlePubMed CentralPubMedGoogle Scholar
- Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75(23):7537–41.View ArticlePubMed CentralPubMedGoogle Scholar
- Rotmistrovsky K, Jang W, Schuler GD. A web server for performing electronic PCR. Nucleic Acids Res. 2004;32(Web Server issue):W108–12.View ArticlePubMed CentralPubMedGoogle Scholar
- Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011;7:539.View ArticlePubMed CentralPubMedGoogle Scholar
- Hamady M, Lozupone C, Knight R. Fast UniFrac: facilitating high-throughput phylogenetic analyses of microbial communities including analysis of pyrosequencing and PhyloChip data. ISME J. 2010;4(1):17–27.View ArticlePubMed CentralPubMedGoogle Scholar
- Team RC. R: A language and environment for statistical computing. 2013.Google Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.