- Research article
- Open Access
Comparative proteomic analysis of different developmental stages of the edible mushroom Termitomyces heimii
Biological Researchvolume 47, Article number: 30 (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
To compare the termite mushroom proteome at each of its three developmental stages, total proteins were extracted from the mycelia, primordia and fruiting bodies. For 2D-DIGE, protein extracts were labelled with fluorochromes and mixed prior to separation via 2D electrophoresis. Figures (1A) and (1B) show representative images of the gels from experiments 1 and 2, respectively, in which the separated proteomes can be observed. We normalised and detected the protein spots with Advanced PDQuest 2-D analysis software (Bio-Rad, Hercules, CA, USA). The optical density of each protein spot was evaluated using Student’s t-test. The results from the two experiments reveal approximately 271 protein spots with molecular masses ranging from 10 kDa to 250 kDa and pIs ranging from 3–10. Each stage has a different and unique proteome map (data not shown). For this preliminary study, and to avoid generating a massive amount of data, we selected 19 protein spots with an average ratio value > 1 and p value <0.01 that were present in both experiments for identification with MALDI-TOF.
The 19 proteins were subjected to MS and subsequent database searches; of these, 12 proteins (represented by 14 spots) had a known function, while five proteins had unknown functions (Table 1). Using functional information from the KEGG and Swiss-Prot pathway databases, we assigned biological functions to the identified proteins (Figure 2).
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 .
Aspartic peptidase A1 belongs to the protease family. We found that it was highly expressed during primordial development (spot #131), but its expression dramatically decreased during fruiting body formation. The enzyme was previously isolated from the basidiomycete fungus Piptoporus soloniensis, but it did not have a clear role in mushroom development. It is possible that aspartic peptidase A1 is down-regulated because the termite mushroom is a non-pathogenic fungus. The same expression pattern was observed for DNA primase (spot #240), although it was highly expressed after the mycelium stage, which may be due to the decrease in meiosis of fruiting body cells after a high level of meiosis during the primordial stage . ppGpp guanosine pentaphosphate (spot #244) was also highly expressed during fruiting body development. The enzyme is an alarmone, which is involved in the stringent bacterial response , resulting in the inhibition of RNA synthesis when there is a shortage of amino acids. To the best of our knowledge, this is the first report describing the up-regulation of ppGpp guanosine pentaphosphate in basidiomycetes. To further understand the protein-protein interactions that occur during mushroom development, we created a protein network (Figure 3) using the STRING protein interaction database (http://string-db.org/) using the identified proteins with known functions as input with a medium level of confidence . The two proteins involved in the glycolytic pathway, phosphoglucomutase and fructose 1,6-bisphosphate aldolase, had the highest number of interactions (7 and 6, respectively). The interaction network revealed the importance of these two proteins in glycolytic regulation, as these two proteins were connected to other enzymes involved in the pathway, such as glucose-6-phosphate isomerase and triose phosphate isomerase. The proteins with the highest number of connections were those that were up-regulated during fruiting body formation. The expression of these proteins may be necessary to develop the rigid cell wall and rapid cell expansion that are characteristic of fruiting body formation . Cysteine peroxiredoxin was highly expressed during fruiting body formation, and we identified a single connection to triose phosphate isomerase. Triose phosphate isomerase is related to the glycolytic pathway, and this interaction may explain the role of cysteine peroxiredoxin in plasma membrane integrity. Together with cell wall formation, the interaction between triose phosphate isomerase and cysteine peroxiredoxin may be necessary for fruiting body formation [14, 24].
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.
Termitomyces heimii were collected during the rainy and wet season at Semenyih (2.9500° N, 101.8500° E), Selangor, Malaysia. Mycelia of T. heimii were cultured on modified Hagem Modess (HM) media and incubated for 3 weeks (Figure 4). Samples at different stages of development (mycelia, primordia and fruiting bodies) were immediately frozen by submerging the samples in liquid nitrogen. The samples were then stored at -80°C for further use.
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.
Natarajan K: South Indian Agaricales V: Termitomyces heimii . Mycologia 1979, 71(4):853-855. 10.2307/3759201
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.2005
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.
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.
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.
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.073
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.200401251
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/pr700617s
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/pr070369o
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-195
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/jf3042159
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-3
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.
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-9
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.033
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-2
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.
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.145
Parker J: Elongation Factors; Translation. In Encyclopedia of Genetics. Edited by: Brenner S, Miller JH. New York: Academic Press; 2001:610-611.
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/jf301198b
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.
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-8
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.
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-663
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.x
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/jf2021495
Anderson E, Burns C, Zolan ME: Global gene expression in Coprinopsis cinerea meiotic mutants reflects checkpoint arrest. G3 (Bethesda) 2012, 2(10):1213-1221. 2012
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.001
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/gkq973
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-065
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.
The authors declare that they have no competing interests.
NR, JRA, NMR and NMS designed the experiment, MHY, NMJ and NSS were involved in sample collection and culture propagation, NR and NMR performed the experiment, JRA, NR, NMR and NBZ analyse the data, JRA and NR wrote the manuscript. All authors read and approved the final manuscript.