Previous Article | Next Article ![]()
Applied and Environmental Microbiology, June 2008, p. 3481-3489, Vol. 74, No. 11
0099-2240/08/$08.00+0 doi:10.1128/AEM.02893-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

School of Natural Resources & Environment, University of Michigan, Ann Arbor, Michigan 48109,1 Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, Michigan 481092
Received 21 December 2007/ Accepted 1 April 2008
|
|
|---|
|
|
|---|
Although fungi are key mediators of plant litter decay in terrestrial ecosystems, we have an incomplete knowledge of how the composition and function of fungal communities change during the process of litter decay. For example, litter is thought to be initially colonized by "sugar fungi" (sensu Garrett, 1963) which have neither cellulolytic or ligninolytic ability. These species are replaced over time by cellulolytic fungi, and finally the litter is colonized by ligninolytic fungi (12); however, the general validity of this scheme has been questioned (11, 20, 29, 43). Moreover, testing this hypothesis has been difficult, because determining the function of fungal species typically requires isolating organisms and examining their effects on defined substrates in pure culture (7, 15, 30-32), which has well-documented limitations (2). The recent development of degenerate primer sets for both basidiomycete (22) and ascomycete laccase genes (25) provides cultivation-independent tools for assessing the genetic diversity and activity of a "lignin-degrading guild" within fungal communities. To our knowledge, a complementary tool to examine the "cellulolytic guild" is still lacking, despite an increasing knowledge of cellulase gene structure (4, 6, 14, 33).
Because cellulose is a large, insoluble molecule with both microcrystalline and amorphous regions, it must be broken down into smaller oligosaccharides before microbial uptake, and this is achieved through the action of extracellular cellulase enzymes (24). There are two types of cellulase, endoglucanases (EC 3.2.1.4), which randomly cleave cellulose molecules, and cellobiohydrolases (EC 3.2.1.91 and EC 3.2.1.-), which remove cellobiose or glucose from the reducing or nonreducing ends of the cellulose polymer (24). These enzymes act in a coordinated and synergistic manner to hydrolyze cellulose into small oligosaccharides and ultimately glucose (24). Cellobiohydrolases are classified into three glycoside hydrolase families (GH6, GH7, and GH48) according to amino acid sequence similarity (16). Of these families, only GH7 is thought to be exclusively fungal, and this family contains CBHI cellobiohydrolases and EG1 endoglucanases from both ascomycete and basidiomycete fungi. Our goal was to develop a general primer set that would permit the recovery of GH7 cellobiohydrolase (cbhI) genes and to use this tool to examine the diversity and distribution of these genes in the forest floor and surface soil of two contrasting forest ecosystems.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Cellobiohydrolase protein sequences used to develop fungcbhI primer set
|
![]() View larger version (27K): [in a new window] |
FIG. 1. Consensus sequence of 40 fungal (CBHI) partial cellobiohydrolase amino acid sequences. The mean pairwise similarity of the alignment was 60.3%, and 83 sites were identical in all sequences. The amino acids shown are those conserved in at least 75% of the sequences. The amino acid motifs in bold type within boxes were used to design degenerate primers fungcbhIF and fungcbhIR. The mean pairwise similarity of the region defined by these two conserved motifs is 67.9%.
|
In addition, forest floor and surface mineral soil were collected from a white oak-black oak (Quercus alba-Q. velutina) forest stand (MST58A; 44°18.60'N, 85°53.83'W) and a sugar maple-basswood (Acer saccharum-Tilia americana) forest stand (MST24; 44°13.40'N, 85°45.31'W) in Wexford County, MI. The soil in the oak forest is an Entic Haplorthod with a pH of 3.9 and 4.4% soil organic C, whereas the maple soil is a Typic Haplorthod with a pH of 5.5 and a carbon content of 5.5% (45). The forest floor of the oak forest has a lignin content of 34% and a cellulose content of 44%. In the maple forest, the forest floor has lower lignin (26%) and slightly higher cellulose (50%) (1). In each stand, six forest floor and six soil samples were collected at approximately 5-m intervals along each of two parallel transects, providing 12 samples of each horizon. Forest floor was collected as a grab sample, and soil was collected to a depth of 5 cm with a 2.5-cm-diameter corer. Forest floor and soil samples were separately combined to create one composite forest floor and one composite soil sample for each stand; these samples were stored on ice prior to DNA extraction.
DNA extraction.
DNA was extracted from sporocarp tissue by a modified cetyltrimethylammonium bromide method after the tissue was ground in Carlson lysis buffer (3). Soil samples were homogenized by hand and passed through a 2-mm sieve to remove roots and coarse material. Ten grams of sieved, homogenized soil was shaken with 100 ml of 0.1 M sodium pyrophosphate for 1 h, and the resulting slurry was then washed through sieves of 0.25-mm and 0.053-mm mesh with
4 liters of deionized water in order to remove spores and maximize the hyphal yield (23). DNA was extracted from the 0.053- to 0.25-mm fraction with PowerMax Soil DNA Isolation kits (MO BIO Laboratories, Inc., Carlsbad, CA). Forest floor samples were chopped (Hamilton Beach R 10-speed blender) to facilitate homogenization, and DNA was extracted from 2.5 g of this material as described above; DNA was stored at –80°C until use.
PCR amplification and cloning.
Approximately 20 ng of sporocarp DNA, or 50 ng of DNA extracted from the two soil horizons, was used as the starting template for PCR amplification. In addition to the template, the reaction cocktail contained 1.5 mM MgCl2, PCR buffer, each deoxynucleoside triphosphate at 0.2 µM, each primer at 0.5 µM, 50 µg of bovine serum albumin, and 1 U of Expand High Fidelity Taq polymerase (Roche); negative controls with no DNA template were included in each batch. PCR was performed by the following protocol: initial denaturation at 94°C for 3 min, followed by 30 (sporocarps) or 35 (soil and forest floor) cycles of 94°C for 30 s, 48°C for 45 s, and 72°C for 90 s, with a final extension at 72°C for 15 min. Success of PCR was assessed by ethidium bromide visualization after agarose gel electrophoresis. We considered any lane showing a band or bands of at least 500 bp to contain a potential cbhI gene fragment.
PCR products were cleaned (PCR Cleanup; MO BIO Laboratories, Inc.) and cloned into pCR2.1-TOPO with the TOPO TA Cloning kit (Invitrogen). For the sporocarp libraries, 12 clones per species were picked and grown overnight in Luria broth before plasmid extraction (Wizard Plus Miniprep kit; Promega, Madison, WI) and sequencing at the University of Michigan Core Sequencing Facility with primers M13F and M13R. For the forest floor and soil samples, 96 clones per ecosystem (48 soil and 48 litter) were grown overnight in Luria broth supplemented with 10% glycerol and sent to the University of Georgia for bidirectional sequencing with primers M13F and M13R.
Sequence analysis.
Sequences were edited with BioEdit, and contigs were constructed in Geneious 3.0.5. For the sequences derived from each sporocarp, a multiple-sequence alignment was used to determine sequence similarity, and for sequences with similarities of >98%, a consensus sequence was derived. The coding region of each DNA sequence was deduced after alignment against known cbhI mRNA sequences (Agaricus bisporus [Z50094], Athelia rolfsii [AB103461], Lentinula edodes [AF411250], Volvariella volvacea [AY559102, AY559103], Aspergillus aculeatus [AB002821], Aspergillus fumigatus [XM745951], and Thermoascus aurantiacus [AY840982]), and partial protein sequences were predicted by translation after removal of introns. The similarity of the putative CBHI protein fragments to known CBHI proteins was assessed with BLAST P. Deduced protein sequences were aligned along with the protein sequences of the reference species (Table 1), as well as the best BLAST P matches, with the ClustalX module in Geneious Pro 3.0.5. The alignment was checked by eye, exported to PAUP*4.0 (38), and analyzed through maximum parsimony. For the maximum-parsimony analysis, all 178 characters in the alignment were included at equal weight, gaps were scored as missing data, and a heuristic search was performed with starting trees obtained by random stepwise addition and TBR branch swapping. Bootstrapping was performed with 100 replications. Trees were rooted with an outgroup composed of the GH7 endoglucanase sequences of Aspergillus oryzae BAE66197 and Aspergillus nidulans AM54071. The diversity of cellulolytic protein in each sample was estimated with p-distance (27). The p-distance values and bootstrapped standard errors were calculated with MEGA 4 (39).
Nucleotide sequence accession numbers.
Newly determined sequences were deposited in GenBank under accession numbers EU345437 to EU345472 and EU359569 to EU359616.
|
|
|---|
![]() View larger version (72K): [in a new window] |
FIG. 2. fungcbhIF-fungcbhIR PCR amplicons recovered from sporocarp and environmental genomic DNA. The markers in the extreme left and right lanes are X174 (Promega). Lanes: 1 and 5, Sarcoscypha occidentalis; 2, Xylaria polymorpha; 3, Clitopilus prunulus; 4, Clitocybe gibba (ALASKA); 6, Clitocybe gibba (MICH); 7, Clitocybe nuda; 8, Chlorophyllum molybdites; 9, oak soil horizon; 10, oak forest floor; 11, sugar maple soil horizon; 12, sugar maple forest floor; 13, no DNA.
|
|
View this table: [in a new window] |
TABLE 2. Cellobiohydrolase (cbhI) gene fragments recovered from sporocarp DNAs of six fungal species and best BLAST P matches of putative CBHI protein fragments recovered after removal of introns and translationa
|
Two sequences 568 and 569 bp in length and 84.1% identical were recovered from the Chlorophyllum molybdites library and termed cbhI568 and cbhI569 (Table 2). Both sequence types had 55-bp introns that differed in splice site (Table 2) and in sequence. Translation recovered 170-amino-acid-long sequences with similarity to Agaricus bisporus exocellobiohydrolase and V. volvacea CBHI-II (Table 2).
For the remaining basidiomycete, Clitopilus prunulus, a 564-bp sequence was recovered (Table 1). The presence of an intron was inferred, and translation recovered a 170-residue protein fragment with similarity to V. volvacea CBHI-II (Table 2).
A 615-bp sequence was recovered in the library of the ascomycete Sarcoscypha occidentalis, and the presence of two short introns was inferred (Table 2). Removal of these introns and subsequent translation produced a 173-amino-acid protein fragment with similarity to the 1,4-β-cellobiosidase of the basidiomycete Pleurotus ostreatus sp. "Florida" (Table 2).
Finally, for the ascomycete Xylaria polymorpha, a 519-bp sequence was recovered. There was no evidence of an intron, and translation recovered a 172-residue putative protein sequence with similarity to the GH7 protein of Sclerotinia sclerotiorum (Table 2).
Environmental clone library.
A total of 160 clones containing inserts of 521 to 692 bp obtained from the four sublibraries (black oak-white oak litter, black oak-white oak soil, sugar maple-basswood litter, and sugar maple-basswood soil) were sequenced, and 50 unique DNA sequences were recovered. Alignments against mRNA revealed that the majority of the sequences contained a single intron; however, introns were absent in seven sequences and five contained two introns. Translation of the 50 sequences uncovered 45 putative amino acid sequences (21 from soil, 23 from litter, and 1 from both soil and litter) ranging in length from 165 to 174 residues, and all 45 of these putative amino acid sequences gave BLAST P matches (65.5 to 90.6% identity) to known GH7 cellobiohydrolase sequences (Table 3). Pairwise similarity of the deduced protein fragments ranged from 54.2 to 99.4%. The p-distance values were 0.257 ± 0.023, 0.294 ± 0.022, 0.322 ± 0.023, and 0.323 ± 0.023 for the oak litter, maple litter, oak soil, and maple soil, respectively.
|
View this table: [in a new window] |
TABLE 3. Cellobiohydrolase genotypes recovered with fungcbh primer set from soil and litter samples of two northern hardwood forests in Michigan
|
![]() View larger version (32K): [in a new window] |
FIG. 3. Phylogenetic relationships among the reference cbhI-encoded protein sequences, environmental clones, and best BLAST P matches. The data set included 178 characters. Shown is 1 of 411 most-parsimonious trees. Bootstrap frequencies of >50% are shown at supported branches. Asterisks indicate branches that collapse in strict consensus. The tree is rooted with an outgroup of GH7 endocellulases.
|
|
|
|---|
Multiple forms of cbhI genes were recovered from nearly all of the fungal species that we examined. The same is true of fungal laccase genes or basidiomycete actin genes recovered from individual species, and this observation has been attributed to a combination of allelic variation and the presence in the fungal genome of distinct gene families derived from different ancestral genes (13, 40). In three of the four basidiomycete species that we examined, we found multiple forms of the cbhI gene that differ in length and DNA sequence, as well as intron size and sequence (Table 1), suggesting that they are distinct cellobiohydrolase genes. These distinct forms encode different CBHI proteins (Fig. 3), further suggesting that they may represent distinct gene families. This result extends previous studies which have documented the presence of cbhI gene families in several fungal species (6, 19, 28), as well as from recently sequenced fungal genomes which have revealed two distinct GH7 genes in Aspergillus terreus, A. fumigatus, and A. flavus, three in Sclerotinia sclerotiorum, six in Coprinopsis cinerea, and eight in Puccinia graminis. Ecologically, having multiple different forms of the cbhI gene may be advantageous, if the proteins encoded by these genes have different environmental optima (42). For example, in the basidiomycete Phanerochaete chrysosporium, which has six forms of the cbhI gene, different patterns of gene expression are observed, depending on whether the fungus is metabolizing powdered cellulose or wood chips (36, 42). In all of the fungal species that we examined, many of the cbhI genes were represented by sequences exhibiting minor coding variation; for example, the Clitocybe nuda cbhI570 gene was represented by three sequence types with a mean pairwise similarity of 98.6%. Although we used a proofreading polymerase, we cannot entirely exclude the possibility of sequencing error. Nevertheless, most of the observed sequence variation was concentrated in third codon positions and putative protein sequences rarely varied by more than a single amino acid. In combination, these observations suggest that allelic forms of the cbhI genes may exist in these species. Further studies with single spore isolates are required to more conclusively determine the allelic status of these cbhI genes.
Cellobiohydrolase genes were recovered from the soil and forest floor of both ecosystems that we sampled. The majority of these environmental cbhI genotypes could not be linked to a known species, which was not unexpected given that fewer than 100 known or presumed cbhI genes reside in public databases and the mean BLAST P score obtained was 79.8% (Table 3). By means of comparison, the similarity of the cbhI genes of Neurospora crassa (Ascomycota, Sordariales) and Phanerochaete chrysosporium (Basidiomycota, Corticiales) across the region defined by the fungcbhI primer set is 73.5%. The best matches we obtained were actually to cbhI genes isolated as part of this study from Clitocybe nuda (100% match to oak clone A01) and Xylaria polymorpha (98.8% to maple clone G12), showing that increased sampling from culture and herbarium collections is required to improve cbhI identification in the future.
Phylogenetic analysis (Fig. 3) provided some insight into the taxonomic affinity of the environmental clones; however, the small number of known sequences relative to the taxonomic breadth of cellulolytic fungi made identification of deep relationships difficult or even impossible. Within the two ecosystems that we studied, more than 70% of the CBHI protein fragments recovered from the forest floor libraries clustered in just three groups (Fig. 3). One group, with good bootstrap support, clustered with Xylaria polymorpha (Xylaria group, Fig. 3). The other two groups, Unknown Forest Floor I and II (Fig. 3), are poorly supported and show the highest similarities (65 to 75% similarity at best) to the CBHI proteins of the basidiomycete Athelia rolfsii and the ascomycete Sclerotinia sclerotiorum. Given the example above of Neurospora crassa and Phanerochaete chrysosporium, the true taxonomic affiliations of these environmental sequences remain unknown. Of the remaining forest floor protein fragments, only one (oak A01, Fig. 3) can be unequivocally identified as a basidiomycete, and only one clustered in a well-supported clade (Alternaria group, Fig. 3). Interestingly, despite the differences in litter biochemistry between oak and maple forests, the cellulolytic communities seem to show a high degree of overlap. Soil-derived CBHI protein fragments from both the maple and oak ecosystems were also placed in the Alternaria group, and this soil horizon of both ecosystems was also represented in a well-supported Trichoderma clade (Trichoderma group, Fig. 3). CBHI protein fragments potentially derived from other microfungal ascomycete species within the Eurotiales, Helotiales, and Sordariales, as well as from Basidiomycetes, were also recovered from the soils of both ecosystems (Fig. 3).
We used p-distance values as a measurement of cellobiohydrolase protein diversity within our ecosystems, and on the basis of these p-distance values, the soil horizons appear to harbor a greater diversity of cellulolytic proteins than the forest floor. Because it is possible that cbhI genes from groups not represented during primer design are not being recovered, comparisons of diversity must be considered speculative. Our primers were designed with 45 of the 62 full-length or nearly full-length cellobiohydrolase sequences curated at the CAZY database (http://www.cazy.org/fam/GH7.html), including all of the sequences for which cellobiohydrolase activity has been confirmed. However, as revealed by inspection of Table 1, this selection is biased toward Ascomycetes, especially from the Sordariomycetes, Eurotiomycetes, and Dothidiomycetes, as well as toward Agaricales within the Basidiomycetes. Nevertheless, our results demonstrate that the primers can amplify cbhI genes from species, such as Sarcoscypha occidentalis (Pezizomycetes), not represented during primer design.
Fungi are key agents of plant litter decay, and we have successfully developed a culture-independent method for studying the community of cellulolytic fungi in soil. We selected the gene for cellobiohydrolase as the target gene because it encodes a key enzyme mediating cellulose decomposition; it also likely plays a critical role in the initial phases of litter colonization and exploitation, when crystalline cellulose is present in the plant cell wall (9). Moreover, by focusing on a single cellulase type, our primers also offer the potential to examine relationships between gene expression and enzyme activity. In addition to cbhI genes, many fungi also possess cbhII genes which encode the cellobiohydrolase enzymes of family GH6. However, coordinated expression of both cbhI and cbhII genes has been observed in the basidiomycete Polyporus arcularius (28), and similar results have been obtained for Trichoderma reesei (17), V. volvacea (19), and Agaricus bisporus (4, 44); thus, focusing solely on the cbhI gene should provide a reasonable indication of cellobiohydrolase potential. The application of the fungcbhI primers in gene expression studies has the potential to provide an indication of fungal cellobiohydrolase activity, as well as to allow the factors that affect cellobiohydrolase activity to be examined in complex communities.
Published ahead of print on 11 April 2008. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»