Search

Microbial volatile organic compounds (mVCs) are formed in the metabolism of microorganisms and widely distributed in nature and pose threats to human health. However, the air pollution by microorganisms is a situation which is poorly understood. In this study, the cytotoxicity of E. aerogenes VCs was evaluated in the model organism Saccharomyces cerevisiae. E. aerogenes VCs inhibited the survival of yeast and triggered the formation of intracellular reactive oxygen species (ROS). The hypersensitive of MAP kinase mpk1/slt2 and 19S regulatory assembly chaperone adc17 mutants to the E. aerogenes VCs indicated cell wall integrity (CWI) pathway together with stress-inducible proteasome assembly regulation are essentially involved in mVCs tolerance mechanism. Furthermore, exposure to the mVCs resulted in the transcriptional upregulation of the CWI pathway, the regulatory particle assembly chaperones, and genes involved in proteasome regulations. Our research suggested that the ROS/MAPK signaling and proteasome regulatory pathway play pivotal roles in the integration and fine-tuning of the mVCs stress response. This study provides a molecular framework for future study of the effects of mVCs on more complex organisms, such as humans.

Phenotypic variation can arise from differences in the protein coding sequence and in the regulatory elements. However, little is known about the contribution of regulatory difference to the expression divergence, especially the cis and trans regulatory variation to the expression divergence in intraspecific populations. In this study, we used two different yeast strains, BY4743 and RM11-1a/α, to study the regulatory variation to the expression divergence between BY and RM under oxidative stress condition. Our results indicated that the expression divergence of BY and RM is mainly due to trans regulatory variations under both normal and oxidative stress conditions. However, cis regulatory variation seems to play a very important role in oxidative stress response in yeast because 36% of genes showed an increase in cis regulatory variation effect compared with 13% of genes that showed an increase in trans regulatory variation effect after oxidative stress. Our data also indicated that genes located on the longer arm of the chromosomes are more susceptible to cis variation effect under oxidative stress than genes on the shorter arm of the chromosomes.

Primers for eight microsatellites were developed; they successfully amplified DNA from 20 domesticated Formosan Sambar deer (Cervus unicolor swinhoei). All loci were polymorphic, with 10-19 alleles per locus. The average observed heterozygosity across loci and samples was 0.310, ranging from 0 to 0.750 at each locus. All loci but one, CU18, deviated from Hardy-Weinberg equilibrium due to excessive homozygosity in these domesticated broodstocks, reflecting inbreeding. These microsatellite loci will be useful, not only for assessment of population structure and genetic variability, but also for conservation of wild deer populations in Taiwan.

Gene regulation change has long been recognized as an important mechanism for phenotypic evolution. We used the evolution of yeast aerobic fermentation as a model to explore how gene regulation has evolved and how this process has contributed to phenotypic evolution and adaptation. Most eukaryotes fully oxidize glucose to CO2 and H2O in mitochondria to maximize energy yield, whereas some yeasts, such as Saccharomyces cerevisiae and its relatives, predominantly ferment glucose into ethanol even in the presence of oxygen, a phenomenon known as aerobic fermentation. We examined the genome-wide gene expression levels among 12 different yeasts and found that a group of genes involved in the mitochondrial respiration process showed the largest reduction in gene expression level during the evolution of aerobic fermentation. Our analysis revealed that the downregulation of these genes was significantly associated with massive loss of binding motifs of Cbf1p in the fermentative yeasts. Our experimental assays confirmed the binding of Cbf1p to the predicted motif and the activator role of Cbf1p. In summary, our study laid a foundation to unravel the long-time mystery about the genetic basis of evolution of aerobic fermentation, providing new insights into understanding the role of cis-regulatory changes in phenotypic evolution.

Bostrycin, a red antibacterial agent produced by Nigrospora sp. no. 407, is considered for meat processing. To optimize production, the culture conditions of submerged fermentation (SmF) and solid-state fermentation (SSF) were investigated. The optimal SmF conditions were a medium containing 1.0% cane molasses and incubation at 30 °C and 150 rpm for 6 days. In SSF, other than bostrycin, less pigment was produced and the optimal ratio of bagasse to water was 1:2 for 10 days. The production and recovery rate of bostrycin by SmF were 120 mg/L and 40%, respectively. Bostrycin exhibited thermostable, pH-dependent color change and dose-dependent antibacterial activity against Clostridium botulinum. Bostrycin-modified meat turned strong red for at least 24 h and could not be removed by washing; bostrycin maintained its antibacterial activity with a bacteriostasis rate of 91% on Staphylcoccus aureus. This is an easy and inexpensive means of acquiring bostrycin from molasses and sugarcane.

Gene expression is regulated both by cis elements, which are DNA segments closely linked to the genes they regulate, and by trans factors, which are usually proteins capable of diffusing to unlinked genes. Understanding the patterns and sources of regulatory variation is crucial for understanding phenotypic and genome evolution. Here, we measure genome-wide allele-specific expression by deep sequencing to investigate the patterns of cis and trans expression variation between two strains of Saccharomyces cerevisiae. We propose a statistical modeling framework based on the binomial distribution that simultaneously addresses normalization of read counts derived from different parents and estimating the cis and trans expression variation parameters. We find that expression polymorphism in yeast is common for both cis and trans, though trans variation is more common. Constraint in expression evolution is correlated with other hallmarks of constraint, including gene essentiality, number of protein interaction partners, and constraint in amino acid substitution, indicating that both cis and trans polymorphism are clearly under purifying selection, though trans variation appears to be more sensitive to selective constraint. Comparing interspecific expression divergence between S. cerevisiae and S. paradoxus to our intraspecific variation suggests a significant departure from a neutral model of molecular evolution. A further examination of correlation between polymorphism and divergence within each category suggests that cis divergence is more frequently mediated by positive Darwinian selection than is trans divergence.

Both cis and trans mutations contribute to gene expression divergence within and between species. We used Saccharomyces cerevisiae as a model organism to estimate the relative contributions of cis and trans variations to the expression divergence between a laboratory (BY) and a wild (RM) strain of yeast. We examined whether genes regulated by a single transcription factor (TF; single input module, SIM genes) or genes regulated by multiple TFs (multiple input module, MIM genes) are more susceptible to trans variation. Because a SIM gene is regulated by a single immediate upstream TF, the chance for a change to occur in its trans-acting factors would, on average, be smaller than that for a MIM gene. We chose 232 genes that exhibited expression divergence between BY and RM to test this hypothesis. We examined the expression patterns of these genes in a BY–RM coculture system and in a BY–RM diploid hybrid. We found that trans variation is far more important than cis variation for expression divergence between the two strains. However, because in 75% of the genes studied, cis variation has significantly contributed to expression divergence, cis change also plays a significant role in intraspecific expression evolution. Interestingly, we found that the proportion of genes with diverged expression between BY and RM is larger for MIM genes than for SIM genes; in fact, the proportion tends to increase with the number of transcription factors that regulate the gene. Moreover, MIM genes are, on average, subject to stronger trans effects than SIM genes, though the difference between the two types of genes is not conspicuous.

The yeast Saccharomyces cerevisiae proliferates rapidly in glucose-containing media. As glucose is getting depleted, yeast cells enter the transition from fermentative to nonfermentative metabolism, known as the diauxic shift, which is associated with major changes in gene expression. To understand the expression evolution of genes involved in the diauxic shift and in nonfermentative metabolism within species, a laboratory strain (BY), a wild strain (RM), and a clinical isolate (YJM) were used in this study. Our data showed that the RM strain enters into the diauxic shift approximately 1 h earlier than the BY strain with an earlier, higher induction of many key transcription factors (TFs) involved in the diauxic shift. Our sequence data revealed sequence variations between BY and RM in both coding and promoter regions of the majority of these TFs. The key TF Cat8p, a zinc-finger cluster protein, is required for the expression of many genes in gluconeogenesis under nonfermentative growth, and its derepression is mediated by deactivation of Mig1p. Our kinetic study of CAT8 expression revealed that CAT8 induction corresponded to the timing of glucose depletion in both BY and RM and CAT8 was induced up to 50- to 90-folds in RM, whereas only 20- to 30-folds in BY. In order to decipher the relative importance of cis- and trans-variations in expression divergence in the gluconeogenic pathway during the diauxic shift, we studied the expression levels of MIG1, CAT8, and their downstream target genes in the cocultures and in the hybrid diploids of BY-RM, BY-YJM, and RM-YJM and in strains with swapped promoters. Our data showed that the differences between BY and RM in the expression of MIG1, the upstream regulator of CAT8, were affected mainly by changes in cis-elements, though also by changes in trans-acting factors, whereas those of CAT8 and its downstream target genes were predominantly affected by changes in trans-acting factors.

Internal duplication can enhance the function of a gene or provide raw material for the emergence of a new function in a gene. Therefore, it is interesting to see whether the frequency of internal duplication has increased during metazoan evolution. The growing number of sequenced eukaryotic genomes provides an excellent opportunity to study the change in the pattern of internal duplication in the course of metazoan evolution. We studied repeated segments in proteins in the proteomes of 11 eukaryotes. We found that the frequency of internal duplication in Caenorhabditis elegans and Drosophila melanogaster (two protostomes) is higher than that in fungi but lower than that in chordates. Moreover, the frequencies of internal duplication for the chordates studied are largely similar. We classified orthologous proteins of chordates into three antiquity groups and found that more recently derived proteins in the metazoan lineage have higher repetitiveness than older ones. Our analysis suggests that lineage-specific internal duplication in protein evolution increases with organismal complexity before the emergence of chordates but not so afterward. Proteins with repeated regions might have been preferred before the protostome-chordate split. This finding supports the suggestion that exon-shuffling occurred more frequently after the first multicellular organism appeared and might have contributed to the metazoan radiation.

Background Comparative genomics provides insights into the diversification of bacterial species. Bacterial speciation usually takes place with lasting homologous recombination, which not only acts as a cohering force between diverging lineages but brings advantageous alleles favored by natural selection, and results in ecologically distinct species, e.g., frequent host shift in Xanthomonas pathogenic to various plants. Results Using whole-genome sequences, we examined the genetic divergence in Xanthomonas campestris that infected Brassicaceae, and X. citri, pathogenic to a wider host range. Genetic differentiation between two incipient races of X. citri pv. mangiferaeindicae was attributable to a DNA fragment introduced by phages. In contrast to most portions of the genome that had nearly equivalent levels of genetic divergence between subspecies as a result of the accumulation of point mutations, 10% of the core genome involving with homologous recombination contributed to the diversification in Xanthomonas, as revealed by the correlation between homologous recombination and genomic divergence. Interestingly, 179 genes were under positive selection; 98 (54.7%) of these genes were involved in homologous recombination, indicating that foreign genetic fragments may have caused the adaptive diversification, especially in lineages with nutritional transitions. Homologous recombination may have provided genetic materials for the natural selection, and host shifts likely triggered ecological adaptation in Xanthomonas. To a certain extent, we observed positive selection nevertheless contributed to ecological divergence beyond host shifting. Conclusion Altogether, mediated with lasting gene flow, species formation in Xanthomonas was likely governed by natural selection that played a key role in helping the deviating populations to explore novel niches (hosts) or respond to environmental cues, subsequently triggering species diversification. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-1369-8) contains supplementary material, which is available to authorized users.

YqjH and YqjW are Bacillus subtilis homologs of the UmuC/DinB or Y superfamily of DNA polymerases that are involved in SOS-induced mutagenesis in Escherichia coli . While the functions of YqjH and YqjW in B . subtilis are still unclear, the comparisons of protein structures demonstrate that YqjH has 36% identity to E . coli DNA polymerase IV (DinB protein), and YqjW has 26% identity to E . coli DNA polymerase V (UmuC protein). In this report, we demonstrate that both YqjH and the products of the yqjW operon are involved in UV-induced mutagenesis in this bacterium. Furthermore, resistance to UV-induced damage is significantly reduced in cells lacking a functional YqjH protein. Analysis of stationary-phase mutagenesis indicates that absences of YqjH, but not that of YqjW, decreases the ability of B . subtilis to generate revertants at the hisC952 allele via this system. These data suggest a role for YqjH in the generation of at least some types of stationary-phase-induced mutagenesis.

To begin characterizing DNA repair capability among Xiphophorus species, we adapted oligonucleotide-based DNA repair assays to extracts of fish tissues. Here, we report the initial results of relative base excision repair (BER) capability among 3 inbred Xiphophorus fish lines representing 2 species (X. maculatus and X. couchianus), and interspecies hybrid F(1) animals produced by crossing them. Overall, data from uracil- N-glycosylase (UNG)-initiated BER assay (UNG-BER) indicate that brain tissue extracts generally possess higher BER activity than do gill and liver extracts. UNG-BER activities in gill and liver extracts were similar. The BER activity in the tissues of F(1) interspecies hybrids followed the activity pattern of the X. couchianus parent in gill and liver extracts, was similar to the X. maculatus parent in brain extracts, but exhibited a reduced repair capacity in gill tissue extracts compared with either parent. We discuss the use of oligonucleotide-based DNA repair assays to elucidate the role that inheritance of DNA repair potential may play in susceptibility to disease and tumorigenesis in the intact organism.

To study the evolutionary relationships among the four living gymnosperm orders and the interfamilial relationships in each order, a set of 65 nuclear 18S rRNA sequences from ferns, gymnosperms, and angiosperms was analyzed using the neighbor-joining and maximum-parsimony methods. With Selaginella as the outgroup, the analysis strongly indicates that the seed plants form a monophyletic group with the ferns as a sister group. Within the seed plants the angiosperms are clearly a monophyletic group. Although the bootstrap support for the monophyly of the gymnosperm clade is moderate, the monophyly is further supported by its lack of angiosperm-specific indels. Within the gymnosperms there appear to be three monophyletic clades: Cycadales-Ginkgoales, Gnetales, and Coniferales. The cycad-ginkgo clade is the earliest gymnosperm lineage. Given the strong support for the sister group relationship between Gnetales and Coniferales, it is unlikely that Gnetales is a sister group of the angiosperms, contrary to the view of many plant taxonomists. Within Coniferales, Pinaceae is monophyletic and basal to the remaining conifer families, among which there are three monophyletic clades: Phyllocladaceae-Podocarpaceae, Araucariaceae, and Sciadopityaceae-Taxaceae-Cephalotaxaceae-Taxodiacea e-Cupressaceae. Within the latter clade, Sciadopityaceae may be an outgroup to the other four families. Among the angiosperms, no significant cluster at the level of subclass was found, but there was evidence that Nymphaeaceae branched off first. Within the remaining angiosperms, the monocots included in this study are nested and form a monophyletic group. This study attests to the utility of nuclear 18S rRNA sequences in addressing relationships among living gymnosperms. Considerable variation in substitution rates was observed among the ferns and seed plants.

Background As a strong fermentator, Saccharomyces cerevisiae has the potential to be an excellent host for ethanol production by consolidated bioprocessing. For this purpose, it is necessary to transform cellulose genes into the yeast genome because it contains no cellulose genes. However, heterologous protein expression in S. cerevisiae often suffers from hyper-glycosylation and/or poor secretion. Thus, there is a need to genetically engineer the yeast to reduce its glycosylation strength and to increase its secretion ability. Results Saccharomyces cerevisiae gene-knockout strains were screened for improved extracellular activity of a recombinant exocellulase (PCX) from the cellulose digesting fungus Phanerochaete chrysosporium. Knockout mutants of 47 glycosylation-related genes and 10 protein-trafficking-related genes were transformed with a PCX expression construct and screened for extracellular cellulase activity. Twelve of the screened mutants were found to have a more than 2-fold increase in extracellular PCX activity in comparison with the wild type. The extracellular PCX activities in the glycosylation-related mnn10 and pmt5 null mutants were, respectively, 6 and 4 times higher than that of the wild type; and the extracellular PCX activities in 9 protein-trafficking-related mutants, especially in the chc1, clc1 and vps21 null mutants, were at least 1.5 times higher than the parental strains. Site-directed mutagenesis studies further revealed that the degree of N-glycosylation also plays an important role in heterologous cellulase activity in S. cerevisiae. Conclusions Systematic screening of knockout mutants of glycosylation- and protein trafficking-associated genes in S. cerevisiae revealed that: (1) blocking Golgi-to-endosome transport may force S. cerevisiae to export cellulases; and (2) both over- and under-glycosylation may alter the enzyme activity of cellulases. This systematic gene-knockout screening approach may serve as a convenient means for increasing the extracellular activities of recombinant proteins expressed in S. cerevisiae.

Background Neocallimastix patriciarum is one of the common anaerobic fungi in the digestive tracts of ruminants that can actively digest cellulosic materials, and its cellulases have great potential for hydrolyzing cellulosic feedstocks. Due to the difficulty in culture and lack of a genome database, it is not easy to gain a global understanding of the glycosyl hydrolases (GHs) produced by this anaerobic fungus. Results We have developed an efficient platform that uses a combination of transcriptomic and proteomic approaches to N. patriciarum to accelerate gene identification, enzyme classification and application in rice straw degradation. By conducting complementary studies of transcriptome (Roche 454 GS and Illumina GA IIx) and secretome (ESI-Trap LC-MS/MS), we identified 219 putative GH contigs and classified them into 25 GH families. The secretome analysis identified four major enzymes involved in rice straw degradation: β-glucosidase, endo-1,4-β-xylanase, xylanase B and Cel48A exoglucanase. From the sequences of assembled contigs, we cloned 19 putative cellulase genes, including the GH1, GH3, GH5, GH6, GH9, GH18, GH43 and GH48 gene families, which were highly expressed in N. patriciarum cultures grown on different feedstocks. Conclusions These GH genes were expressed in Pichia pastoris and/or Saccharomyces cerevisiae for functional characterization. At least five novel cellulases displayed cellulytic activity for glucose production. One β-glucosidase (W5-16143) and one exocellulase (W5-CAT26) showed strong activities and could potentially be developed into commercial enzymes.

Agricultural Biotechnology Research Center, Academia Sinica, Taipei, Taiwan18.1 Introduction18.2 Next generation: renewable energy biomass program18.3 Main feedstocks for next-generation biofuels18.4 Identification of cellulase genes by genomic approaches18.5 How can transcriptomic study help identify cellulase genes in a microbe?18.6 How can proteomic study help identify cellulase genes in a microbe?18.7 Conversion of cellulose to ethanol18.8 Concluding remarksReferences

Previously, using a chromosomal reversion assay system, we established that an adaptive mutagenic process occurs in nongrowing Bacillus subtilis cells under stress, and we demonstrated that multiple mechanisms are involved in generating these mutations (41, 43). In an attempt to delineate how these mutations are generated, we began an investigation into whether or not transcription and transcription-associated proteins influence adaptive mutagenesis. In B. subtilis , the Mfd protein (transcription repair coupling factor) facilitates removal of RNA polymerase stalled at transcriptional blockages and recruitment of repair proteins to DNA lesions on the transcribed strand. Here we demonstrate that the loss of Mfd has a depressive effect on stationary-phase mutagenesis. An association between Mfd mutagenesis and aspects of transcription is discussed.

Adaptive (stationary-phase) mutagenesis occurs in the gram-positive bacterium Bacillus subtilis . Furthermore, taking advantage of B. subtilis as a paradigm for the study of prokaryotic differentiation and development, we have shown that this type of mutagenesis is subject to regulation involving at least two of the genes that are involved in the regulation of post-exponential phase prokaryotic differentiation, i.e., comA and comK . On the other hand, a functional RecA protein was not required for this type of mutagenesis. The results seem to suggest that a small subpopulation(s) of the culture is involved in adaptive mutagenesis and that this subpopulation(s) is hypermutable. The existence of such a hypermutable subpopulation(s) raises important considerations with respect to evolution, the development of specific mutations, the nature of bacterial populations, and the level of communication among bacteria in an ecological niche.

We utilized a custom-synthesized double-strand oligonucleotide containing a single O 6 -methylguanine ( O 6 -MG) residue within a restriction endonuclease recognition site to determine O 6 -methylguanine-DNA-methyltransferase ( O 6 -MGMT) activity in various tissue extracts prepared from Xiphophorus fish. The results suggest Xiphophorus fish O 6 -MGMT activity has many of the same characteristics as Escherichia coli and mammalian O 6 -MGMT’s including rapid reaction kinetics consistent with stoichiometric removal of methyl groups, but exhibits a temperature optimum of 23°C. Results from protein extract activity assays indicate O 6 -MGMT activity patterns among four Xiphophorus tissues followed the order: brain≥testes>gill≥liver. In mammals, O 6 -MGMT activity is high in liver, while activity in brain is minimal (i.e. ≈9% of liver); however, we report that in the Xiphophorus fishes examined, brain tissue extracts exhibited much higher (≈six-fold) O 6 -MGMT activity levels than liver. Comparison of O 6 -MGMT activity between Xiphophorus species employed in tumor induction experiments did not indicate significant differences in ability to clear the pre-mutagenic O 6 -MG from the oligonucleotide substrate.

During an investigation of the parameters controlling mutations in Bacillus subtilis we observed that this bacterium exhibits a transient growth requirement for two nonessential amino acids (glutamic acid and isoleucine) during a type of postexponential growth on a minimal medium.

To determine the evolutionary positions of the conifer genera Amentotaxus, Phyllocladus, and Nageia, we obtained 18S rRNA sequences from 11 new taxa representing the major living orders and families of gymnosperms. With the published Chlamydomonas as an outgroup, phylogenetic analyses of our new data and available sequences indicate that (1) the Gnetales form a monophyletic group, which is an outgroup to the conifers, (2) the conifers are monophyletic, (3) Taxaceae, Cephalotaxaceae, Cupressaceae, and Taxodiaceae form a monophyletic group, (4) Amentotaxus is closer to Torreya than to Cephalotaxus, suggesting that Amentotaxus is better to be classified as a member of Taxaceae, (5) Phyllocladus, Dacrycarpus, Podocarpus, and Nageia form a monophyletic group, and (6) Pinaceae is an out-group to the other families of conifers. Our finding that Phyllocladus is a sister group of the Podocarpaceae disagrees with the suggestion that the phylloclade of the genus is an ancient structure and that the genus is a terminal taxon within the Podocarpaceae. The genus Nageia is more closely related to Podocarpus than to Dacrycarpus and was derived from within the Podocarpaceae. In conclusion, our data indicate that in conifers, the uniovulate cone occurred independently in Taxacaeae and Cephalotaxaceae, and in Podocarpaceae after the three families separated from Pinaceae, and support the hypothesis that the uniovulate cone is derived from reduction of a multiovulate cone.

Phenotypic variation among individuals in a population can be due to DNA sequence variation in protein coding regions or in regulatory elements. Recently, many studies have indicated that mutations in regulatory elements may be the major cause of phenotypic evolution. However, the mechanisms for evolutionary changes in gene expression are still not well understood. Here, we studied the relative roles of cis and trans regulatory changes in Saccharomyces cerevisiae cells to cope with heat stress. It has been found that the expression level of ~300 genes was induced at least two fold and that of ~500 genes was repressed at least two fold in response to heat shock. From the former set of genes, we randomly selected 65 genes that showed polymorphism(s) between the BY and RM strains for pyrosequencing analysis to explore the relative contributions of cis and trans regulatory variations to the expression divergence between BY and RM. Our data indicated that the expression divergence between BY and RM was mainly due to trans regulatory variations under either the normal condition or the heat stress condition. However, the relative contribution of trans regulatory variation was decreased from 76.9% to 61.5% after the heat shock stress. These results indicated that the cis regulatory variation may play an important role in the adaption to heat stress. In our data, 43.1% (28 genes) of the 65 genes showed the same trend of cis or trans variation effect after the heat shock stress, 35.4% (23 genes) showed an increased cis variation effect and 21.5% (14 genes) showed an increased trans variation effect after the heat shock stress. Thus, our data give insights into the relative roles of cis and trans variations in response to heat shock in yeast.

Background 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. Results 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. Conclusions 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. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-1422-7) contains supplementary material, which is available to authorized users.