Comparative proteomic analysis of different developmental stages of the edible mushroom Termitomyces heimii
© Rahmad et al.; licensee BioMed Central Ltd. 2014
Received: 18 March 2014
Accepted: 25 June 2014
Published: 3 July 2014
Termitomyces heimii is a basidiomycete fungus that has a symbiotic relationship with termites, and it is an edible mushroom with a unique flavour and texture. T. heimii is also one of the most difficult mushrooms to cultivate throughout the world. Little is known about the growth and development of these mushrooms, and the available information is insufficient or poor. The purpose of this study was to provide a base of knowledge regarding the biological processes involved in the development of T. heimii. The proteomic method of 2 dimensional difference gel electrophoresis 2D-DIGE was used to determine and examine the protein profiles of each developmental stage (mycelium, primordium and fruiting body). Total proteins were extracted by TCA-acetone precipitation.
A total of 271 protein spots were detected by electrophoresis covering pH 3–10 and 10–250 kDa. Selected protein spots were subjected to mass spectrometric analyses with matrix-assisted laser desorption/ionisation (MALDI TOF/TOF). Nineteen protein spots were identified based on peptide mass fingerprinting by matching peptide fragments to the NCBI non-redundant database using MASCOT software. The 19 protein spots were categorised into four major groups through KEGG pathway analysis, as follows: carbohydrate metabolism, energy metabolism, amino acid metabolism and response to environmental stress.
The results from our study show that there is a clear correlation between the changes in protein expression that occur during different developmental stages. Enzymes related to cell wall synthesis were most highly expressed during fruiting body formation compared to the mycelium and primordial stages. Moreover, enzymes involved in cell wall component degradation were up-regulated in the earlier stages of mushroom development.
The edible mushroom Termitomyces heimii is an agaric-type fungus in the family Lyophyllaceae. It was first discovered in India during the late 1970s . The fungus is also known as “termite mushroom” due to its symbiotic relationship with termites, insects belonging to the subfamily Macrotermitinae. These basidiomycete mushrooms are distributed throughout tropical and subtropical areas of Africa and Asia. Termitomyces fungi are cultivated on a special substrate in the termite nest called the fungus comb . In most Macrotermes termites, the fungus comb is composed of their primary faeces, which is mostly undigested dead plant matter . After the termites grow, they feed on the mature parts of the fungus comb and aggregated asexual spores, called fungal nodules. Dead plant materials are thought to be efficiently and strongly degraded in the symbiotic system because fungus-growing termites, unlike many other termites, do not produce faeces rich in organic matter . Thus, the symbiotic fungi are responsible for the efficient decomposition of plant materials by fungus-growing termites. Termitomyces heimii is valuable because it has a delicious taste. In addition, the fungus has potential therapeutic uses in blood pressure and blood lipid regulation, the immune response, apoptosis and infection . Its cultivation remains challenging, and numerous worldwide attempts have been made with varying degrees of success. However, all attempts have fallen short as economical cultivation methods . There are limited reports on the biochemical composition of the fungus. This may be due to the difficulty in obtaining wild fruit bodies, as these are seasonal mushrooms. Studies to determine the mechanisms involved in the mushroom’s development are necessary to provide insight into potential fungus cultivation processes. Proteomic approaches, such as two-dimensional difference gel electrophoresis (2D-DIGE) in combination with mass spectrometry (MS), have been widely applied to identify and profile proteins expressed in pathogenic [7, 8] and edible basidiomycetes . A number of studies have convincingly demonstrated that DIGE and MS are complementary proteomics approaches, and the use of both approaches in combination provides deeper insight into the proteome. To better understand the mushroom’s systematic response to environmental stress and protein-protein interactions during different developmental stages, we employed proteomic approaches in this study.
2D-DIGE PAGE at different stages of development
Proteins identified from the differential 2D-DIGE analysis obtained from 3 different developmental stages (mycelium, primordium and fruiting bodies) of the edible mushroom Termitomyces heimii
Protein score C.I%
Total ion score C.I%
Average ratio value
Predicted protein [Laccaria bicolor S238N-H82]
NADP-glutamate dehydrogenase [Laccaria bicolor]
NADP-glutamate dehydrogenase [Laccaria bicolor]
Translation elongation factor 1-alpha [Uromyces polygoni-avicularis]
Cysteine peroxiredoxin [Laccaria bicolor S238N-H82]
Carbohydrate esterase family 4 protein [Laccaria bicolor S238N-H82]
Cysteine peroxiredoxin [Laccaria bicolor S238N-H82]
Phosphoglucomutase [Micromonaspusilla CCMP1545]
Fructose 1,6-bisphosphate aldolase [Laccaria bicolor S238N-H82]
Predicted protein [Laccaria bicolor S238N-H82]
Putative sulfite reductase [Janibacter sp. HTCC2649]
Aspartic peptidase A1 [Laccaria bicolor S238N-H82]
Hypothetical protein MPER_03344 [Moniliophthoraperniciosa FA553]
DNA primase [Oceanicolagranulosus HTCC2516]
(p)ppGppsynthetase I SpoT/RelA [Actinobacillussuccinogenes 130Z]
Hypothetical protein MPER_10776 [Moniliophthoraperniciosa FA553]
Metallophosphoesterase [Spirosomalinguale DSM 74]
Hypothetical protein Phum_PHUM233470 [Pediculushumanuscorporis]
NADH dehydrogenase (quinone) [Opitutus terrae PB90-1]
During the development of all basidiomycete fungi, the vegetative mycelium stage transforms into a fruiting body via an intermediate primordium stage, which is a complicated and critical process . Understanding the molecular mechanisms underlying this process has long been a goal in basidiomycete research. By defining the proteome at each developmental stage, we can obtain a broader understanding of the mushroom fruiting process .
This study is one of the first reports of the termite mushroom proteome. Unlike classical proteomic studies that produce large quantities of information, we screened different developmental stages of the termite mushroom for proteins with significantly but consistently altered expression levels during development. Therefore, we focused on 19 spots that passed our screening criteria. NADP-dependent glutamate dehydrogenase (NADP-GDH) is an intermediate enzyme in fungal ammonium absorption and converts 2-oxoglutarate to glutamate . NADP-GDH is positioned at the border of carbon and nitrogen metabolism and is potentially involved in the regulation of these processes . Although the cofactor enzyme NADH dehydrogenase (spot #173) was down-regulated during fruiting body development, we found that NADP-GDH (spot #109 and 106) was slightly up-regulated in the primordia stage compared to the mycelium. However, when we compared its expression between fruiting bodies and primordia, it was similarly, but not significantly, up-regulated. Metallophosphoesterase (spot #242) had a similar expression profile. It was not significantly down-regulated during the primordial stage; however, it was highly up-regulated in fruiting bodies. This was expected, as the levels of cell wall enzymes likely increase after fruiting bodies are harvested in basidiomycetes . Carbohydrate esterase 4 (Spot #175) hydrolyses proteins, especially pectins. We found that it was significantly up-regulated during primordial formation, but its expression significantly decreased during fruiting body formation. It was previously reported that carbohydrate esterase 4 was up-regulated in fungal cells during their interactions with their host cells [15, 16]. In addition, Yang et al.  found that carbohydrate esterase was expressed during fruiting body formation of the closely related edible mushroom Termitomyces albuminosus, which suggests that termite mushrooms can degrade lignin. The mushroom cooperates with the enzymes produced by the fungus-growing termite to degrade the cellulose and hemicellulose of plant cells, which supports the well-coordinated cooperation hypothesis and potentially explains the symbiotic relationship between the termite and the fungus .
We found that translation elongation factor 1-alpha (TEF α1) was up-regulated in the primordial stage compared to the mycelium. While the expression of the protein in the fruiting stage was not significantly down-regulated compared to its expression in the mycelium stage, additional translation factors are likely involved. TEF α1 is thought to play an important and major role in the regulation of cell wall morphogenesis . Similarly, TEF α1 was up-regulated in the basidiomycete Pleurotus tuberregium during the transition from mycelium to primordium .
Two spots (#199 and #200) represented cysteine peroxiredoxin. Peroxiredoxins are a class of fungal antioxidants that reduce hydroperoxides to alcohols . We found that the protein was highly expressed in the primordia and was significantly higher in fruiting bodies compared to the mycelium. Cysteine peroxiredoxin usually localises to the cytosol and mitochondria. In eukaryotic cells, cysteine peroxiredoxin is highly expressed in various isoforms and, despite being less efficient than catalase, cysteine peroxiredoxin likely contributes to the protection of the plasma membrane against lipid peroxide oxidation . This may be why cysteine peroxiredoxin levels increase during fungal development. Two different proteins involved in the glycolytic pathway, phosphoglucomutase (spot #61) and fructose 1,6 bisphosphate aldolase (spot #146), did not significantly change during the primordia stage; however, both proteins were significantly up-regulated during fruiting body formation. Insignificant changes in the levels of proteins related to the glycolytic pathway may lead to changes in glycolysis, gluconeogenesis, and mitochondrial function, which may be beneficial for nascent fungal cells adapting to an anaerobic environment . Increased levels of these proteins may be necessary for the development of the rigid cell wall during fruiting body formation . Sulphite reductase (spot #126) is another enzyme we detected that slightly changes during primordia formation and significantly changes during fruiting body formation. Although studies have been performed on bacterial cells, there is little information available on the functional role of sulphite reductase in mushroom development, except that it plays a role in sulphur assimilation in yeast during fermentation .
The identification of proteins with potential roles in T. heimii mushroom development could provide useful information regarding fruiting body formation, which is possibly one of the most important developmental events in the life cycle of this mushroom, as it allows the formation of propagative structures. While altered protein expression may be physiologically relevant to the fungus and can facilitate changes in metabolism via different pathways, it may also be directly involved in the fruiting body development of T. heimii, a rare, highly desirable, edible seasonal fungus. We are currently working to obtain more information on the proteome of each developmental stage of the fungus by utilising various proteomic tools, such as liquid chromatography-based separation and label-free quantitative proteomics. These tools can be used to ultimately improve the cultivation methods of T. heimii.
Separate samples were ground using a mortar and pestle in liquid nitrogen. Protein extraction was carried out using TCA acetone precipitation . Pellets were resuspended in lysis buffer containing 30 mM Tris–HCl, pH 8.8, 7 M urea, 2 M thiourea and 4% CHAPS. Protein quantitation was conducted with a Bradford kit (Bio-Rad, USA) with bovine serum albumin as a standard.
Two experiments with 3 biological replicates each were performed. Each proteome was labelled with an appropriate fluorochrome (Cy2 for the internal standard, Cy3 for the mycelial protein and Cy5 for the primordial protein). In the second experiment, the proteins were labelled similarly to those in the 1st experiment, except that Cy5 was used to label the fruiting body protein as shown in Figure 1B. For each sample, 30 μg protein was mixed with 1.0 μl diluted CyDye, and kept in the dark on ice for 30 min. Samples from each pair were labelled with Cy3 and Cy5, respectively, and an internal standard with Cy2 was run on each gel. The labelling reaction was terminated by adding 1.0 μl 10 mM lysine to each sample and incubating in the dark on ice for an additional 15 min. The labelled samples were then mixed together. The 2× 2-D sample buffer (8 M urea, 4% CHAPS, 20 mg/ml DTT, 2% pharmalytes and a trace amount of bromophenol blue), 100 μl destreak solution and rehydration buffer (7 M urea, 2 M thiourea, 4% CHAPS, 20 mg/ml DTT, 1% pharmalytes and a trace amount of bromophenol blue) were added to the labelling mixture to bring the total volume to 250 μl. The solution was mixed well and spun down before loading the labelled samples into a strip holder.
IEF and SDS-PAGE
The isoelectric focusing step (IEF) was carried out on a PROTEAN® i12™ IEF cell (BioRad) by applying the samples to 13 cm Immobiline™ DryStrip gels (IPG), with a linear pH range of 3–10 (GE Healthcare). Rehydration was performed at 250 V for 2 hours in the linear mode followed by 250 V constant for 2 hours, 250–5000 V for 4 hours in the linear mode and 5000 V constant for a total of 8 hours and 35 kvh. Upon completion of IEF, the IPG strips were incubated in freshly made equilibration buffer-1 (50 mM Tris–HCl, pH 8.8, 6 M urea, 30% glycerol, 2% SDS, a trace amount of bromophenol blue and 10 mg/ml DTT) for 15 minutes with gentle shaking. The strips were then rinsed in freshly made equilibration buffer-2 (50 mM Tris–HCl, pH 8.8, 6 M urea, 30% glycerol, 2% SDS, a trace amount of bromophenol blue and 45 mg/ml DTT) for 10 minutes with gentle shaking. Next, the IPG strips were rinsed in SDS gel running buffer and placed on the top of the SDS gels (12% gels). The SDS gels were run at 15°C until the dye front ran out of the gels.
Image scan and data analysis
Gel images were immediately scanned following SDS-PAGE using Typhoon TRIO (GE Healthcare). The scanned images were analysed with PDQuest 2-D analysis software (Bio-Rad), and multiplex gel images were generated by cross-gel analysis using DeCyder software (version 6.5, GE Healthcare). Quantitative and statistical analyses of changes in the spots were performed with DeCyder software (GE Healthcare).
Protein identification by mass spectrometry
Spot picking and trypsin digestion
Spots of interest were picked up using the Ettan Spot Picker (Amersham BioSciences) based on the in-gel analysis and spot picking design by DeCyder software. The gel spots were washed several times and digested in-gel with modified porcine trypsin protease (Trypsin Gold, Promega). The digested tryptic peptides were desalted by Zip-tip C18 (Millipore). Peptides were eluted from the Zip-tip with 0.5 μl matrix solution (α-cyano-4-hydroxycinnamic acid, 5 mg/ml in 50% acetonitrile, 0.1% trifluoroacetic acid, 25 mM ammonium bicarbonate) and spotted on the AB SCIEX MALDI plate (Opti-TOF™ 384 well insert).
MALDI-TOF MS and TOF/TOF tandem MS/MS were performed on an AB SCIEX TOF/TOF™ 5800 System (AB SCIEX, Framingham, MA). MALDI-TOF mass spectra were acquired in reflectron positive ion mode, averaging 4000 laser shots per spectrum. TOF/TOF tandem MS fragmentation spectra were acquired for each sample, averaging 4000 laser shots per fragmentation spectrum on each of the 7–10 most abundant ions present in each sample (excluding trypsin autolytic peptides and other known background ions).
Both the resulting peptide mass spectra and associated fragmentation spectra were submitted to a GPS Explorer workstation equipped with the MASCOT search engine (Matrix science) to search the non-redundant database of the National Center for Biotechnology Information (NCBInr). Searches were performed without constraining protein molecular weight or isoelectric point, with variable carbamidomethylation of cysteine and oxidation of methionine residues and with one allowed missed cleavage site. Candidates with either a protein score C.I.% or ion C.I.% more than 95 were considered significant.
2 Dimensional difference gel electrophoresis
- MALDI TOF/TOF:
Matrix-assisted laser desorption/ionization/ time-of-flight
Kyoto encyclopedia of genes and genomes
Immobilized pH gradient
Search Tool for the retrieval of interacting genes.
This project was funded by the grant from Ministry of Science Technology and Innovation MOSTI (Grant No.08-05-ABI PB030); the authors would like to thank Dr. Mohd Nazir Basiran (ABI Executive Director) for his kind assistant and help to facilitate this work.
- Natarajan K: South Indian Agaricales V: Termitomyces heimii . Mycologia 1979, 71(4):853-855. 10.2307/3759201View ArticleGoogle Scholar
- Taprab Y, Johjima T, Maeda Y, Moriya S, Trakulnaleamsai S, Noparatnaraporn N, Ohkuma M, Kudo T: Symbiotic fungi produce laccases potentially involved in phenol degradation in fungus combs of fungus-growing termites in Thailand. Appl Environ Microbiol 2005, 71(12):7696-7704. 10.1128/AEM.71.12.7696-7704.2005PubMed CentralView ArticlePubMedGoogle Scholar
- Wood TG, Thomas RJ: The mutualistic association between Macrotermitinae and Termitomyces. In Insect-Fungus Interactions. Edited by: Wilding N, Collins NM, Hammond PM, Webber JF. London: Academic Press; 1998:69-92.Google Scholar
- Rouland-Lefèvre C: Termites as Pests of Agriculture. In Biology of Termites: a Modern Synthesis. Edited by: Bignell DE, Roisin Y, Lo N. Netherlands: Springer; 2011:499-517.Google Scholar
- Abd Malek SN, Kanagasabapathy G, Sabaratnam V, Abdullah N, Yaacob H: Lipid Components of a Malaysian Edible Mushroom, Termitomyces heimii Natarajan. Int J Food Prop 2011, 15(4):809-814.View ArticleGoogle Scholar
- Tsai SY, Wu TP, Huang SJ, Mau JL: Nonvolatile taste components of Agaricus bisporus harvested at different stages of maturity. Food Chem 2007, 103: 1457-1464. 10.1016/j.foodchem.2006.10.073View ArticleGoogle Scholar
- Shimizu M, Yuda N, Nakamura T, Tanaka H, Wariishi H: Metabolic regulation at the tricarboxylic acid and glyoxylate Cycles of the lignin-degrading basidiomycetes Phanerochaete chrysosporium against exogenous addition of vanillin. Proteomics 2005, 5: 3919-3931. 10.1002/pmic.200401251View ArticlePubMedGoogle Scholar
- Matsuzaki F, Shimizu M, Wariishi H: Proteomic and metabolomic analyses of the white-rot fungus Phanerochaete chrysosporium exposed to exogenous benzoic acid. J Proteome Res 2008, 7: 2342-2350. 10.1021/pr700617sView ArticlePubMedGoogle Scholar
- Horie K, Rakwal R, Hirano M, Shibato J, Nam HW, Kim YS, Kouzuma Y, Agrawal GK, Masuo Y, Yonekura M: Proteomics of two cultivated mushrooms sparassis crispa and hericium erinaceum provides insight into their numerous functional protein components and diversity. J Proteome Res 2008, 7(5):1819-1835. 10.1021/pr070369oView ArticlePubMedGoogle Scholar
- Cheng CK, Au CH, Wilke S, Stajich J, Zolan M, Pukkila P, Kwan HS: 5'-Serial Analysis of Gene Expression studies reveal a transcriptomic switch during fruiting body development in Coprinopsis cinerea . BMC Genomics 2013, 14(1):195. 10.1186/1471-2164-14-195PubMed CentralView ArticlePubMedGoogle Scholar
- Lau CC, Abdullah N, Shuib AS, Aminudin N: Proteomic analysis of antihypertensive proteins in edible mushrooms. J Agric Food Chem 2012, 60(50):12341-12348. 10.1021/jf3042159View ArticlePubMedGoogle Scholar
- Baars JJ, Op den Camp HJ, van der Drift C, Joordens JJ, Wijmenga SS, Van Griensven LJ, Vogels GD: 15 N-NMR study of ammonium assimilation in Agaricus bisporus . Biochim Biophys Acta 1996, 1310(1):74-80. 10.1016/0167-4889(95)00157-3View ArticlePubMedGoogle Scholar
- Noor S, Punekar NS: Allosteric NADP-glutamate dehydrogenase from aspergilli: purification, characterization and implications for metabolic regulation at the carbon-nitrogen interface. Microbiology 2005, 151(Pt 5):1409-1419.View ArticlePubMedGoogle Scholar
- Sakamoto Y, Nakade K, Sato T: Characterization of the post-harvest changes in gene transcription in the gill of the Lentinula edodes fruiting body. Curr Genet 2009, 55(4):409-423. 10.1007/s00294-009-0255-9View ArticlePubMedGoogle Scholar
- Spanikova S, Biely P: Glucuronoyl esterase–novel carbohydrate esterase produced by Schizophyllum commune. FEBS Lett 2006, 580(19):4597-4601. 10.1016/j.febslet.2006.07.033View ArticlePubMedGoogle Scholar
- van den Brink J, De Vries RP: Fungal enzyme sets for plant polysaccharide degradation. Appl Microbiol Biotechnol 2011, 91(6):1477-1492. 10.1007/s00253-011-3473-2PubMed CentralView ArticlePubMedGoogle Scholar
- Yang F, Xu B, Li J, Huang Z: Transcriptome analysis of Termitomyces albuminosus reveals the biodegradation of lignocellulose. Wei Sheng Wu Xue Bao 2012, 52(4):466-477.PubMedGoogle Scholar
- Ohkuma M, Shimizu H, Thongaram T, Kosono S, Moriya K, Trakulnaleamsai S, Noparatnaraporn N, Kudo T: An alkaliphilic and xylanolytic Paenibacillus species isolated from the gut of a soil-feeding termite. Microbes Environ 2003, 18(3):145-151. 10.1264/jsme2.18.145View ArticleGoogle Scholar
- Parker J: Elongation Factors; Translation. In Encyclopedia of Genetics. Edited by: Brenner S, Miller JH. New York: Academic Press; 2001:610-611.View ArticleGoogle Scholar
- Chen L, Zhang B-B, Cheung PCK: Comparative proteomic analysis of mushroom cell wall proteins among the different developmental stages of pleurotus tuber-regium. J Agric Food Chem 2012, 60(24):6173-6182. 10.1021/jf301198bView ArticlePubMedGoogle Scholar
- Belozerskaia TA, Gessler NN: [Reactive oxygen species and the strategy of the antioxidant defense in fungi: a review]. Prikl Biokhim Mikrobiol 2007, 43(5):565-575.PubMedGoogle Scholar
- Wood ZA, Schroder E, Robin Harris J, Poole LB: Structure, mechanism and regulation of peroxiredoxins. Trends Biochem Sci 2003, 28(1):32-40. 10.1016/S0968-0004(02)00003-8View ArticlePubMedGoogle Scholar
- Cobos R, Barreiro C, Mateos RM, Coque JJ: Cytoplasmic- and extracellular-proteome analysis of Diplodia seriata: a phytopathogenic fungus involved in grapevine decline. Proc Natl Acad Sci U S A 2010, 8: 46.Google Scholar
- Patyshakuliyeva A, Jurak E, Kohler A, Baker A, Battaglia E, De Bruijn W, Burton KS, Challen MP, Coutinho PM, Eastwood DC, Gruben BS, Mäkelä MR, Martin F, Nadal M, van den Brink J, Wiebenga A, Zhou M, Henrissat B, Kabel M, Gruppen H, de Vries RP: Carbohydrate utilization and metabolism is highly differentiated in Agaricus bisporus . BMC Genomics 2013, 14: 663. 10.1186/1471-2164-14-663PubMed CentralView ArticlePubMedGoogle Scholar
- Cordente AG, Heinrich A, Pretorius IS, Swiegers JH: Isolation of sulfite reductase variants of a commercial wine yeast with significantly reduced hydrogen sulfide production. FEMS Yeast Res 2009, 9(3):446-459. 10.1111/j.1567-1364.2009.00489.xView ArticlePubMedGoogle Scholar
- El-Baky HA, Linke D, Nimtz M, Berger RG: PsoP1, a milk-clotting aspartic peptidase from the basidiomycete fungus Piptoporus soloniensis . J Agric Food Chem 2011, 59(18):10311-10316. 10.1021/jf2021495View ArticlePubMedGoogle Scholar
- Anderson E, Burns C, Zolan ME: Global gene expression in Coprinopsis cinerea meiotic mutants reflects checkpoint arrest. G3 (Bethesda) 2012, 2(10):1213-1221. 2012View ArticleGoogle Scholar
- Srivatsan A, Wang JD: Control of bacterial transcription, translation and replication by (p)ppGpp. Curr Opin Microbiol 2008, 11(2):100-105. 10.1016/j.mib.2008.02.001View ArticlePubMedGoogle Scholar
- Szklarczyk D, Franceschini A, Kuhn M, Simonovic M, Roth A, Minguez P, Doerks T, Stark M, Muller J, Bork P, Jensen LJ, Von Mering C: The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucleic Acids Res 2011, 39: D561-D568. 10.1093/nar/gkq973PubMed CentralView ArticlePubMedGoogle Scholar
- Lakshman DK, Natarajah SS, Garett WM, Dhar AK: Optimized protein extraction methods for proteomic analysis of Rhizoctimia solani. Mycologia 2008, 100(6):867-875. 10.3852/08-065View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. 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.