Nucleotide Sequencing and Identification of Some Wild Mushrooms

Das, Sudip Kumar; Mandal, Aninda; Datta, Animesh K.; Gupta, Sudha; Paul, Rita; Saha, Aditi; Sengupta, Sonali; Dubey, Priyanka Kumari

doi:https://doi.org/10.1155/2013/403191

The Scientific World Journal

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2013 | Article ID 403191 | https://doi.org/10.1155/2013/403191

Nucleotide Sequencing and Identification of Some Wild Mushrooms

Sudip Kumar Das,¹Aninda Mandal,¹Animesh K. Datta,¹Sudha Gupta,²Rita Paul,³Aditi Saha,⁴Sonali Sengupta,⁵and Priyanka Kumari Dubey¹

Academic Editor: X. Liang, J. Yoon

Received26 Aug 2013

Accepted14 Oct 2013

Published30 Dec 2013

Abstract

The rDNA-ITS (Ribosomal DNA Internal Transcribed Spacers) fragment of the genomic DNA of 8 wild edible mushrooms (collected from Eastern Chota Nagpur Plateau of West Bengal, India) was amplified using ITS1 (Internal Transcribed Spacers 1) and ITS2 primers and subjected to nucleotide sequence determination for identification of mushrooms as mentioned. The sequences were aligned using ClustalW software program. The aligned sequences revealed identity (homology percentage from GenBank data base) of Amanita hemibapha [CN (Chota Nagpur) 1, % identity 99 (JX844716.1)], Amanita sp. [CN 2, % identity 98 (JX844763.1)], Astraeus hygrometricus [CN 3, % identity 87 (FJ536664.1)], Termitomyces sp. [CN 4, % identity 90 (JF746992.1)], Termitomyces sp. [CN 5, % identity 99 (GU001667.1)], T. microcarpus [CN 6, % identity 82 (EF421077.1)], Termitomyces sp. [CN 7, % identity 76 (JF746993.1)], and Volvariella volvacea [CN 8, % identity 100 (JN086680.1)]. Although out of 8 mushrooms 4 could be identified up to species level, the nucleotide sequences of the rest may be relevant to further characterization. A phylogenetic tree is constructed using Neighbor-Joining method showing interrelationship between/among the mushrooms. The determined nucleotide sequences of the mushrooms may provide additional information enriching GenBank database aiding to molecular taxonomy and facilitating its domestication and characterization for human benefits.

1. Introduction

Mushrooms are defined as “a macro fungi with a distinctive fruiting body” [1] and are traditionally used worldwide as nutritious food and as medicinal sources [2–7] including antioxidant activity [8, 9]. Further, mushrooms are great recyclers and decomposers [10] and therefore play a significant role in the ecosystem. To date, about 3000 species are regarded as prime edible mushrooms [11]. Proper identification knowledge of edible mushrooms is essential for effective exploration in human benefits.

Molecular markers, PCR (Polymerase Chain Reaction) [12] and non-PCR based [13], are widely used for mushroom identification and characterization. However, direct sequencing of PCR product of repetitive nuclear DNA [14–16] of mushrooms is powerful tool for identification and phylogenetic studies [17]. The present study describes the nucleotide sequencing of 8 wild edible mushrooms collected from Eastern Chota Nagpur Plateau of West Bengal, India, using genomic DNA from fruit bodies. The rDNA-ITS (Ribosomal DNA Internal Transcribed Spacers) fragments of the genomic DNA were amplified using ITS1 (Internal Transcribed Spacers 1) and ITS2 primers. The nucleotide sequences of 8 mushrooms were matched from the available known sequences of GenBank database. The ITS region is rather useful for molecular characterization in fungi at the species level and within the species [11]. The objective of the work is to gain proper identification knowledge of the mushrooms, which may provide direction towards domestication and commercialization of the wild species for economic benefits apart from aiding molecular taxonomy.

2. Materials and Methods

2.1. Germplasms

Eight morphologically different types of wild edible mushrooms (CN (Chota Nagpur) [1–8]) were collected from their natural habitat in Eastern Chota Nagpur Plateau of West Bengal (Paschim Medinipur having sandy loam to loamy soil of reddish brown color cover; latitude 22°N to 25°30′N, longitude 83°47′E to 87°50′E, altitude 610 m) during the rainy (June-July) and postrainy (August to October) seasons of 3 consecutive years (2008 to 2010). The mushroom samples were characterized morphologically following the methodology suggested by Largent and Stuntz [18]. The samples verified up to genus level [19] were as follows: CN 1—Amanita sp.; CN 2—Amanita sp.; CN 3—Astraeus sp.; CN 4—Termitomyces sp.; CN 5—Termitomyces sp.; CN 6—Termitomyces sp.; CN 7—Termitomyces sp.; and CN 8—Volvariella sp.

2.2. DNA Isolation

Field isolated samples were dried (45°C ± 1°C) and the dried samples were used for fungal DNA isolation from fruit bodies using the Fungal Genomic DNA MiniSpin Kit (Chromous Biotech, Bangalore) as per methodology suggested by Graham et al. [20]. Isolated DNA samples were checked in agarose gel (1% agarose gel prepared in TAE (Tris-acetate-ethylenediaminetetraacetic acid) buffer to which 2 μL ethidium bromide was added) to confirm the yield of DNA, purity, and concentration. A ladder (1 kb plus) was loaded in the gel in order to compare the size of the isolated DNA samples.

2.3. Primer Designing

Based on morphological identification of the samples, available nucleotide sequences ranging from 2 to 5 kb were downloaded from GenBank data base for primer designing. The FASTA format of all those sequences was arranged in a file and then the sequences were aligned according to their homology of nucleotide using ClustalW software program (http://www.ebi.ac.uk/Tools/msa/clustalw2/). The stretches of conserved sequences (ITS1 and ITS2) were considered for primer designing. The length of the primers was kept within 18 to 25 bp and the GC% was within 40 to 55% so that the melting temperature of the primers will not be high during PCR. It was also considered that the amplicon size should not cross more than 500 bp. Forward and reverse primers (Table 1) for 8 mushroom samples were synthesized commercially (IDT, USA) for amplification of fungal DNA.

2.4. PCR Reactions

Polymerase chain reaction was performed using the following: 10x PCR buffer 2.5 μL, dNTPs (25 mM of each) 2.0 μL, MgCl₂ (25 mM) 2.0 μL, forward primer (10 pm/μL) 1.0 μL, reverse primer (10 pm/μL) 1.0 μL, DNA template 10 ng for each sample, Taq. polymerase (2 U/μL) 0.5 μL, and DNase- and RNase-free water where volume was adjusted to 25 μL. The PCR conditions were determined according to Taq. polymerase, primer pair, and amplicon size. The amplification reactions were performed in a DNA Thermal Cycler (Eppendorf AG, Hamburg, Germany) programed as follows: 1st cycle of 5 min at 95°C (initial denaturation) followed by 30 cycles of 45 sec at 95°C (denaturation), 30 sec at 50°C (annealing), 1 min at 72°C (extension), and 1 cycle of 10 min at 72°C (final extension). The final step was held at 4°C. The PCR products were stored at −20°C and subsequently run in agarose (Hi-media, USA; in TAE buffer) gel (1%). PCR product (10 μL) was loaded in the gel by adding 6x DNA loading dye along with 2-log DNA ladder. The gel was run at 50 volts for 1 hour and the amplified products were visualized in the UV trans-illuminator and photographed in Gel Doc system (Bio-Rad, USA).

2.5. Sequencing

The PCR products were loaded in the coupled lane of low melting point agarose gel (Sigma, USA). The run was at 40 volts for 1 hour and the bands were cut using sterile blade. The DNA was extracted from the agarose using columns (QIAGEN, Germany). This DNA was quantified in 1% agarose gel and then sequenced using the respective forward and reverse primers. The kit which was used for DNA sequencing was Big Dye Terminator Kits, Applied Biosystems, USA and the sequencing was done in Capillary Electrophoresis DNA Sequencer, Applied Biosystems, USA.

2.6. Evolutionary Relationships Analysis

Following alignment of the nucleotide sequences, a phylogenetic tree was constructed using the Neighbor-Joining method [21]. The optimal tree is drawn with the sum of branch length equal to 1.4246. The tree is drawn to scale, with branch lengths in the same unit as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the Maximum Composite Likelihood method [22] and are in the unit of the number of base substitutions per site. The analysis involved 8 nucleotide sequences. Codon positions included were 1st + 2nd + 3rd + noncoding. All positions containing gaps and missing data were eliminated. There were a total of 102 positions in the final data set. Evolutionary analyses were conducted in Molecular Evolutionary Genetics Analysis 5 (MEGA5) [23].

3. Results

The size of the DNA of 8 mushroom samples was around 15 kb. The PCR amplification products (Figure 1) showed that CN 1 and 2 gave around 200 to 250 bp amplified band, while CN 3 and 6 had 300 bp amplified band. CN 4, 5, and 7 documented a band around 500 bp, while CN 8 showed 400 bp band. These PCR products were gel purified, run in 1% agarose gel, and processed for nucleotide sequencing.

The partial nucleotide sequences (Table 2) of 8 samples were obtained and analyzed for Basic Local Alignment Search Tool (BLAST) search program (National Center for Biotechnology Information (NCBI) site) against the whole GenBank data base of nucleotide sequences (using ClustalW software) for identification. Each PCR product was sequenced using forward and reverse primers and both were used for BLAST. The BLAST result is presented in Table 3. Based on % identity, mushrooms may be identified as follows: CN 1—Amanita hemibapha [family: Amanitaceae; % identity 99 (JX844716.1)]; CN 2—Amanita sp. [% identity 98 (JX844763.1)]; CN 3—Astraeus hygrometricus [Family: Diplocystaceae; % identity 87 (FJ536664.1)]; CN 4—Termitomyces sp. [Family: Lyophyllaceae; % identity 90 (JF746992.1)]; CN 5—Termitomyces sp. [% identity 99 (GU001667.1)]; CN 6—T. microcarpus [% identity 82 (EF421077.1)]; CN 7—Termitomyces sp. [% identity 76 (JF746993.1)]; and CN 8—Volvariella volvacea [Family: Pluteaceae; % identity 100 (JN086680.1)]. Out of 8 samples, 4 belong to Termitomyces spp. of which only 1 could be identified up to species level from the available GenBank data base. Molecular identification of the mushroom samples corroborated morphological identification of the samples up to genus level (Figures 2(a)–2(h)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Phylogenetic tree (Figure 3) constructed revealed a close relationship between CN 1 and 2 (A. hemibapha and Amanita sp.), 4 and 5 (Termitomyces sp.), and 3 and 6 (Astraeus hygrometricus and T. microcarpus). CN 8 (V. volvacea) showed close affinity with CN 1, 2, 4, and 5; however, CN 7 (Termitomyces sp.) showed distant relationship with other samples.

4. Discussion

The nucleotide sequences of 8 mushroom samples were blasted against available sequences from GenBank data base for identification. One sample, CN 8 (Volvariella sp.), matched 100% Indian sample V. volvacea, accession number JN086680.1, reported from Solan, Himachal Pradesh. On the other hand, the rest of the samples matched the samples that were reported from outside of India. However, CN 1 identified as A. hemibapha (pileus color: yellow) from GenBank data base matched morphologically the reported mushroom A. hemibapha (pileus color: white) by Roy and Samajpati [24] from Bankura district, West Bengal, adjoining to the current location of study. On the contrary, CN 2 was also identified as Amanita sp. that possessed white color pileus. Therefore, CN 1 and 2 could be different Amanita species as evinced from nucleotide sequences. CN 6 was identified to be as T. microcarpus from the available sequence of GenBank data base (82% homology). CN 4 and 5 belonging to Termitomyces spp. seem to be of genetic lineage different from that of the other Termitomyces spp. (CN 6 and 7) as noted from the phylogenetic tree; however, all Termitomyces species were collected from the same locality. From molecular point of view CN 4 and 5 were identified as Termitomyces spp., and these samples matched morphologically T. clypeatus and T. eurhizus [25], respectively. Tang et al. [26] also characterized Termitomyces species (T. longiradicatus, T. quilonensis, and T. poonensis) different morphologically from India. From phylogenetic tree, V. volvacea showed distinct relationship to A. hemibapha and Amanita sp. although the fungi belong to different families and therefore the finding is rather interesting from taxonomic point of view and may further be explored. The present study possibly suggests that GenBank data base for mushroom is not sufficiently rich in India. However, molecular characterization of 4 mushroom samples up to species level is performed and is essential and the provided nucleotide sequences of the rest of the samples may be relevant to GenBank data base for further exploration in the field.

Rajaratnam and Thiagarajan [11] extracted genomic DNA from the fruit body of wild mushroom and subjected it to nucleotide sequencing using ITS1 and ITS4 conserved primer stretches. The sequence was aligned using Jukes-Cantor Correlated Distance model and the aligned sequence (559 bp) revealed 88% matched score with Perenniporia sp. (GQ982890.1). Dung et al. [27] characterized 6 Oyster mushroom samples based on both morphological and molecular data and found that both were corroborating to each other. Lee et al. [16] performed identification of 3 medicinal mushroom (Ganoderma lucidum, Coriolus versicolor, and Fomes fomentarius) species from Korea based on nuclear large subunit rDNA sequences. Nucleotide sequencing of Termitomyces albuminosus [28], Ganoderma lucidum [29], and Agaricus bisporus [30], among others was also performed.

5. Conclusion

Molecular identification of the fungal samples may enrich and provide additional information to mushroom biodiversity and GenBank data base resource aiding to molecular phylogenetic analysis. Further, identification knowledge may also be significant for human benefits.

Acknowledgments

The research is grant-aided by Science & Society Division, Department of Science and Technology, Government of India. The authors sincerely acknowledge the kind help of Dr. Soma Chattopadhyay, Scientist, Institute of Life Sciences (Government of India), during work and preparation of the paper.

References

S. T. Chang and P. G. Miles, “Mushroom biology—a new discipline,” The Mycologist, vol. 6, pp. 64–65, 1992.
View at: Google Scholar
T. N. Kaul, Introduction to Mushroom Science: Systematics, Oxford & IBH Publishing, New Delhi, India, 1997.
P. Manzi, L. Gambelli, S. Marconi, V. Vivanti, and L. Pizzoferrato, “Nutrients in edible mushrooms: an inter-species comparative study,” Food Chemistry, vol. 65, no. 4, pp. 477–482, 1999.
View at: Publisher Site | Google Scholar
S. G. Jonathan and I. O. Fasidi, “Antimicrobial activities of two Nigerian edible macrofungi-Lycoperdonpusilum (Bat. Ex) and Lycoperdongigantum(Pers),” African Journal of Biomedical Research, vol. 6, pp. 85–90, 2003.
View at: Google Scholar
O. O. Kuforiji, I. O. Fasidi, and S. A. Odunfa, “Nutritive value of Pleurotus tuberregium cultivated on different agro-industrial wastes,” Nigerian Journal of Microbiology, vol. 17, pp. 68–70, 2003.
View at: Google Scholar
J. Gbolagade and I. Fasidi, “Antimicrobial activities of some selected Nigerian mushrooms,” African Journal of Biomedical Sciences, vol. 8, no. 2, pp. 83–87, 2005.
View at: Google Scholar
T. Yaltirak, B. Aslim, S. Ozturk, and H. Alli, “Antimicrobial and antioxidant activities of Russula delica Fr.,” Food and Chemical Toxicology, vol. 47, no. 8, pp. 2052–2056, 2009.
View at: Publisher Site | Google Scholar
L. Barros, M. Ferreira, B. Queirós, I. C. F. R. Ferreira, and P. Baptista, “Total phenols, ascorbic acid, β-carotene and lycopene in Portuguese wild edible mushrooms and their antioxidant activities,” Food Chemistry, vol. 103, no. 2, pp. 413–419, 2007.
View at: Publisher Site | Google Scholar
D. Kumari, M. S. Reddy, and R. C. Upadhyay, “Nutritional composition and antioxidant activities of 18 different wild Cantharellus mushrooms of Northwestern Himalayas,” Food Science and Technology International, vol. 17, no. 6, pp. 557–567, 2011.
View at: Publisher Site | Google Scholar
C. Manoharachary, K. Sridhar, R. Singh et al., “Fungal biodiversity: distribution, conservation and prospecting of fungi from India,” Current Science, vol. 89, no. 1, pp. 58–71, 2005.
View at: Google Scholar
S. Rajaratnam and T. Thiagarajan, “Molecular characterization of wild mushroom,” European Journal of Experimental Biology, vol. 2, no. 2, pp. 369–373, 2012.
View at: Google Scholar
G. G. Fonseca, E. A. Gandra, L. F. Sclowitz, A. P. A. Correa, J. A. J. Costa, and J. A. Levy, “Oyster mushrooms species differentiation through molecular markers RAPD,” International Journal of Plant Breeding and Genetics, vol. 2, pp. 13–18, 2008.
View at: Google Scholar
A. Pawlik, G. Janusz, J. Koszerny, W. Małek, and J. Rogalski, “Genetic diversity of the edible mushroom Pleurotus sp. by amplified fragment length polymorphism,” Current Microbiology, vol. 65, no. 4, pp. 438–445, 2012.
View at: Google Scholar
D. J. Royse, B. A. Bunyard, and M. S. Nicholson, “Molecular genetic analysis of diversity in populations of edible mushrooms,” in Mushroom Biology and Mushroom Products, S. T. Chang, J. A. Buswell, and S. Chiu, Eds., pp. 49–54, The Chinese University Press, Hong Kong, 1993.
View at: Google Scholar
L. Savard, P. Li, S. H. Strauss, M. W. Chase, M. Michaud, and J. Bousquet, “Chloroplast and nuclear gene sequences indicate Late Pennsylvanian time for the last common ancestor of extant seed plants,” Proceedings of the National Academy of Sciences of the United States of America, vol. 91, no. 11, pp. 5163–5167, 1994.
View at: Publisher Site | Google Scholar
J. S. Lee, M. O. Lim, K. Y. Cho, J. H. Cho, S. Y. Chang, and D. H. Nam, “Identification of medicinal mushroom species based on nuclear large subunit rDNA sequences,” Journal of Microbiology, vol. 44, no. 1, pp. 29–34, 2006.
View at: Google Scholar
M. Nei, Molecular Evolutionary Genetics, Columbia University Press, New York, NY, USA, 1987.
D. L. Largent and D. E. Stuntz, How to Identify Mushrooms to Genus I: Macroscopic Features, Mad River Press, Revised edition, 1986.
S. K. Das, Morphological, biochemical and molecular characterization of some wild edible mushroom of Eastern Chotonagpur part of West Bengal [Ph.D. thesis], University of Kalyani, West Bengal, India, 2013.
G. C. Graham, P. Mayers, and R. J. Henry, “A simplified method for the preparation of fungal genomic DNA for PCR and RAPD analysis,” BioTechniques, vol. 16, no. 1, pp. 48–50, 1994.
View at: Google Scholar
N. Saitou and M. Nei, “The neighbor-joining method: a new method for reconstructing phylogenetic trees,” Molecular Biology and Evolution, vol. 4, no. 4, pp. 406–425, 1987.
View at: Google Scholar
K. Tamura, M. Nei, and S. Kumar, “Prospects for inferring very large phylogenies by using the neighbor-joining method,” Proceedings of the National Academy of Sciences of the United States of America, vol. 101, no. 30, pp. 11030–11035, 2004.
View at: Publisher Site | Google Scholar
K. Tamura, D. Peterson, N. Peterson, G. Stecher, M. Nei, and S. Kumar, “MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods,” Molecular Biology and Evolution, vol. 28, no. 10, pp. 2731–2739, 2011.
View at: Publisher Site | Google Scholar
A. Roy and N. Samajpati, “Agaricales of West Bengal. III,” Indian Journal of Mycological Research, vol. 18, pp. 13–24, 1980.
View at: Google Scholar
R. P. Purkayastha and A. Chandra, Manual of Indian Edible Mushrooms, pp. 12–225, Today and Tomorrow's Publication, New Delhi, India, 1985.
B.-H. Tang, T.-Z. Wei, and Y.-J. Yao, “Type revision of three Termitomyces species from India,” Mycotaxon, vol. 94, pp. 93–102, 2005.
View at: Google Scholar
N. T. P. Dung, D. B. Tuyen, and P. H. Quang, “Morphological and genetic characteristics of oyster mushrooms and conditions effecting on its spawn growing,” International Food Research Journal, vol. 19, no. 1, pp. 347–352, 2012.
View at: Google Scholar
F. Yang, B. Xu, S. Zhao et al., “De novo sequencing and analysis of the termite mushroom (Termitomyces albuminosus) transcriptome to discover putative genes involved in bioactive component biosynthesis,” Journal of Bioscience and Bioengineering, vol. 114, no. 2, pp. 228–231, 2012.
View at: Google Scholar
S. Chen, J. Xu, C. Liu et al., “Genome sequence of the model medicinal mushroom Ganoderma lucidum,” Nature Communications, vol. 3, article 913, 2012.
View at: Google Scholar
E. Morin, A. Kohler, A. R. Baker et al., “Genome sequence of the button mushroom Agaricus bisporus reveals mechanisms governing adaptation to a humic-rich ecological niche,” Proceedings of the National Academy of Sciences, vol. 109, no. 43, pp. 17501–17506, 2012.
View at: Google Scholar

Copyright

Copyright © 2013 Sudip Kumar Das et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

6693

Downloads

3227

Citations