Abstract
Previous studies for enumerating fungal communities on pine needles relied entirely on phenotypic diversity (microscopy) or identification based on DNA sequence data from those taxa recovered via cultural studies. To bypass limitations of the culturing methods and provide a more realistic diversity estimate, we employed and assessed a PCR-DGGE based method coupled with rDNA phylogenetic sequence analyses to characterize fungal taxa associated with pine needles. Fresh (living) and decayed needles from three hosts of the Pinaceae (Keteleeria fortunei, Pinus elliottii and P. massoniana) were examined. Morphological studies reveal that the most abundant species associated with decayed needles were Cladosporium cladosporioides and an unidentified Trichoderma species followed by Gliocephalotrichum sp., Gliocladium sp., Lophodermium pinastri, Paecilomyces varioti, Phaeostalagmus cyclosporus and a Phoma sp, which are commonly occurring fungi. Community genomic data from freshly collected and decayed pine needles recovered 40 operational taxonomic units, which appear to be mostly undetected members of the natural fungal consortium. Sequence analyses revealed a number of phylotypes or “species” that were not recovered using traditional morphological and cultural approaches previously used. Phylogenetic data from partial 18S rDNA sequence data reveal that most phylotypes represent potential novel phylogenetic fungal lineages with affinities to the Dothideomycetes, Leotiomycetes, Lecanoromycetes and Sordariomycetes and were not identical to previously known endophytes or saprobes. Although the major ecological roles of these phylotypes in pine needles are still enigmatic, this study provides new insights in hidden fungal diversity that mycologists are possibly ignoring given the discrepancies associated with available methods. To what extent do previously recovered identified species (either as saprobes or endophytes) from morphological or culturing studies act as pioneer decomposers or constitute an integral part of endophytic community warrants further investigation.
Similar content being viewed by others
Introduction
Fungi, especially from pines, are important contributors to primary productivity and nutrient cycling. They are mostly isolated as endophytes, pathogens or saprobes and are critical components in plant ecosystems1,2,3. Several studies have focused on the characterization of fungi colonising pine needles from ecological and pathological perspectives (e.g.4,5,6,7,8,9,10). The pioneering work on mycota of pine litter by Kendrick and Burges4 prompted other researchers to focus on this host, which is an important component of many forests. Traditional methods documenting fungal diversity, however, rely mostly on taxonomic identification of species based on spore morphology and DNA sequence data of endophytes from cultural dependent methods. This has a major drawback given that not all fungi grow readily in highly selective medium and sporulation is dependent on species, host, seasonality, growth conditions, and other environmental factors which can largely underestimate diversity.
In many cases, fungi exist predominantly as vegetative hyphae, which make microscopic examination difficult3,7. In addition, the use of artificial media does not reflect natural conditions and only those microorganisms that thrive well under these conditions will grow. It is thought that only a small fraction (0.1 to 10%) of microorganisms existing in the nature can be cultured artificially11. The use of previous classical methods to enumerate fungal taxa or characterise fungal communities from pine needles has therefore underestimated diversity. Generally, most studies were limited to temperate regions and have consistently isolated only a few dominant fungi, the most common of which are Lophodermium pinastri and L. piceae12. Recently, several studies based on DNA sequence analyses, especially those targeting endophytes, have also revealed different groups of fungal species not recovered via traditional methods (e.g.8,9,10).
To circumvent the challenges associated with spore production and cultivation, a number of molecular based approaches have been widely used (as detailed in Anderson et al.13). Among these, PCR-DGGE has proven to be a powerful technique to detect and characterize fungal communities in soil samples, plant material (roots and leaves), rocks, water, alga and glass13,14,15,16,17. Studies have revealed that dominant microorganisms isolated from environmental samples using traditional and cultivation-dependent methodology are genetically different from those identified via DGGE analysis of rDNA gene (e.g.18). Diversity of freshwater fungi in decaying leaves of white birch (Betula papayrifera), several maple species (Acer rubrum, A. saccharum, and A. spicatum) and white spruce (Picea glauca) has been reported using a traditional and a molecular approach19. They found that species richness estimates based on DGGE are higher than those from T-RFLP analysis and much higher than those based on spore identification. Subsequent DGGE analyses of aquatic hyphomycetes have also revealed that presence or absence of fungal taxa and biomass from community DNA is regulated primarily by season and type of substrates20. Based on sequence analyses of DGGE bands, Vainio & Hantula15 successfully detected a higher fungal diversity from environmental samples of Picea abies (Norway spruce stumps) than when using culture techniques. Cullings et al.21 used DGGE to evaluate changes in soil communities in response to the effects of artificial defoliation of pines in a pine-spruce forest. From a pathological perspective, DGGE has also been valuable in profiling fungal taxa that are associated with healthy and infected black spruce seedlings (Picea marina)22.
Our previous studies on pine needles from Keteleeria fortunei, Pinus elliottii and Pinus massoniana (Family: Pinaceae) in Hong Kong revealed a few taxa new to science23,24 and a large number of unidentified taxa, especially those that exist predominantly in their asexual stages and sterile hyphae. To the best of our knowledge, there is no known studies on pine needles of these 3 members of the Pinaceae that have analysed the fungal species (especially those colonizing internal tissues) based on PCR-DGGE followed by 18S rDNA sequence analyses. Mahajan et al.9 have attempted to characterize microbial communities from Pinus roxburghii of the Himalayan subtropical region (India) based on PCR-DGGE using ITS-GC clamped primers but all PCR amplifications failed and no fungal taxa could be sequenced. Aneja et al.25 used DGGE to characterize bacteria and fungi but their focus was detection of potential saprobic fungi that are potential decomposers in spruce planted in vitro. In this study we attempt to recover fungal species found in pine needles using a polyphasic approach (direct morphological comparison, cultural based methods and PCR-DGGE analyses of 18S rDNA coupled with phylogenetic analyses). The objectives of this study were to (1) profile the main fungal species associated with different decay stages and (2) find out whether fungal species from pine needles are more diverse than expected and, (3) whether “taxa” recovered from the molecular approach are related to previously identified endophytes or whether they constitute unknown phylogenetic lineages.
Results
Morphological-based methodology
A total of 39 taxa were identified from the three hosts and 2 needle types (Table 1). A large number of anamorphic (asexual) fungi (77%) were recovered, while the average occurrence of sexual morph ascomycetes was quite low (6.25%). The most abundant taxa were Cladosporium cladosporioides, Gliocephalotrichum sp., Gliocladium sp., Paecilomyces varioti, Phaeostalagmus cyclosporus, Phoma sp., and a Trichoderma sp. The percentage of non-sporulating taxa was 14.58%. A proportion of fungi (14.5%) existed as cultures that could not be identified as no sporulating structures were observed. These were termed as “Mycelia sterilia” (sensu3,26) and have the highest frequency of occurrence. Lophodermium pinastri has the highest frequency of occurrence among the three ascomycetes present. Anamorphic fungi were dominant on both, L and F needle types (with 100% dominance on Keteleeria fortunei (L-type needles), 94.74% on K. fortunei (F-type needles), 96.77% on Pinus elliottii and 93.75% on P. massoniana).
The highest species diversity as calculated by Shannon diversity index was recorded from L-type needles on Pinus elliottii (1.21), followed by P. massoniana (L-type) (1.20), Keteleeria fortunei (L-type) (1.08) and the lowest was on K. fortunei (F-type) (1.04). Cluster analyses (Bray-Curtis similarity) and multi-dimensional scaling (MDS) were applied to show the similarities and differences of the fungal communities on Keteleeria fortunei, Pinus elliottii and P. massoniana. The L-type needles and F-type needles of Keteleeria fortunei were similar in fungal diversity. Samples were interspersed in three different clades rather than forming distinct clades based on needle types. Dendograms generated also indicated that species composition from Pinus elliottii and P. massoniana were different from those of Keteleeria fortunei. Samples from Pinus elliottii and P. massoniana (both F and L needle types) are assembled in one cluster whereas samples from Keteleeria fortunei (both F and L needle types) were found to form two different clusters without any relationships to needle types. Results also reveal that species composition from L-type needles of Pinus elliottii and P. massoniana was very similar and have similar fungal diversity but clades corresponding to either host were not observed. In addition, we also noted that there was no clear relationship between fungal diversity and seasonal variation among Pinus elliottii and P. massoniana (results not shown). For instance, samples collected from winter (Nov-Mar) and summer (May-Sep) seasons were found in all the clades. However, in the case of Keteleeria fortunei, samples collected in winter appeared in distinct clades and were distant from those collected in summer months.
Molecular-Based Methodology
DGGE analyses
At least two sets of DGGE profiles were generated in order to ensure the results obtained were reproducible. All bands that migrated at different positions were sequenced to ensure that the maximum number of phylotypes was identified. This also enabled us to compare whether bands migrating at different positions are phylogenetically similar or dissimilar. The most appropriate denaturant gradient, which gave better resolution, was 5–55% urea/formamide. The first DGGE analysis revealed 35 bands (Fig. 1) while the second one revealed 36 bands (Fig. 2). Most of the bands found in Fig. 1 are also present in Fig. 2 but there were a few exceptions: bands 1, 2, 3, 10, 17 and 35, present in the first gel (Fig. 1) were not observed in the second gel (Fig. 2). In addition, bands G, H, V, W, Z, A3, A8 and A9 were present in the second gel but not in the first gel. Bands were selected and excised based on their positions and mobility from both gels. Forty-four bands were sequenced and 39 were successful (87% success).
Some bands were weak but this did not affect the quality of sequences as these bands were purified and reamplified again prior to sequencing. Fungal diversity of Keteleeria fortunei (especially on L-type needles) was the highest among the three Pinaceae species as more taxa were detected. Eight, 16 and 2 bands were recovered from fresh (living), L-type needles and F1-type needles of Keteleeria fortunei respectively. DGGE of DNA samples from Pinus massoniana resulted in 9 bands (4 and 5 from living and L-type needles respectively). The lowest phylotypic diversity was recorded from Pinus elliottii (3 and 1 from living and L-type needles respectively).
Sequence analysis
Based on blast search results, 30 phylotypes (75%) were similar to known bitunicate fungi (Dothideomycetes); 2 were similar to the Eurotiales; 5 similar to the Letiomycetes and 1 to the Lecanoromycetes; 1 was a unitunicate taxon and similar to the Hypocreales. All sequences display reasonable similarity (94% or more) to fungal sequences in GenBank. Those having highest sequence similarity were included in the phylogenetic analyses. Analyses based on different criteria yielded essentially similar trees and the phylogram generated from Bayesian analyses with bootstrap support generated under Maximum Likelihood is presented here. Some bands present in one needle type were also found in the other needle types. For example, bands B, D and F from L-type needles of Keteleeria fortunei share similar mobility gradient as bands A5 and A6, recovered from fresh needles of Pinus elliottii. However, blast search results and phylogenies indicate that co-migrating bands are not necessarily similar or related. For instance, bands A or 33 (Figs 1 and 2 respectively; from the living needles of Keteleeria fortunei), M or 31 (from L-type needles of K. fortunei) and 3 (Fig. 1; from L-type needles of Pinus massoniana) were found to be at the same position but sequence analyses revealed that they were totally different taxa. Band M or 31 was found to be highly similar and related to the Rhytismatales (Letiomycetes), while band A or 33 was similar and related to Cochliobolus sativus (Pleosporales) and band 3 was similar and related to an uncultured clone, whose classification within the Dothideomycetes, is unknown. Similar results were obtained with bands 17, L (or 32) and A4 (or 20), which migrated to the same position on the gel. Co-migrating bands were also found within one needle type but results also indicate that they are not be similar or related. For instance, bands O (or 13) and 10 were from L-type needles of Keteleeria fortunei but sequence comparison revealed that band O (or 13) and 10 were dissimilar and phylogenetically unrelated. Band O is related to an uncultured clone (Dothideomycetes) whereas band 10 was a sister taxon to the Rhytismatales. Similar results were obtained with bands A (or 33) and C (or 15) recovered from living needles of Keteleeria fortunei and L (or 32) and Q (or 12) from L-type needles of Keteleeria fortunei. This phenomenon clearly indicates that bands (phylotypes) co-migrating at the same position have distinct evolutionary linages. This is not surprising as it has already been demonstrated in previous studies that phylogenetically distant taxa can have co-migrating bands and that one band does not necessarily mean one unique phylotype27,28,29.
Sequences obtained from respective DGGE bands were analysed under different optimality criteria as stated before. All trees generated yielded same topologies with regards to the position of our phylotypes (except in one ML analysis where the position of the Eurotiomycetes and Dothideomycetes were slightly different but this relationship was unstable with no support). The Bayesian tree, rooted with Saccharomyces is selected to infer phylogenetic position of our phylotypes (Fig. 3). Phylogeny shows that most of the phylotypes sequenced form a distinct and moderately supported monophyletic group (78% BS) basal to other known bitunicate ascomycetes, especially to the family Sympoventuriaceae (Fig. 3). In all phylogenies, phylotype 2 PM clusters with a few cultured Scolecobasidium clones with strong support while phylotype Q12 KF is nested in between Verruconis and Ochroconis species. Phylotypes A (33), E (16), and 1 KF belong to the family Pleosporaceae and are closely related to Cochlobolius sativus and related species (Pleosporaceae). Phylogram also depicts that our two phylotypes (G KF and 17 KF) constitute an independent and strongly supported subclade (100% BS) within the genus Penicillium (Eurotiales, Eurotiomycetes) while uncultured isolate H KF is related to Metus conglomeratus and Lecanora hybocarpa and is best accommodated within the Lecanorales. The unitunicate phylotype P KF clusters with Geosmithia putterilli with moderate support and also shares close affinities to Gliocladium viride (Hypocreales). Phylotypes 10 KF, C15 KF, M31 KF, J29 KF and L32 KF have close phylogenetic affinities to other anamorphic Leotiomycetes (especially Lophodermium pinastri and other relatives of the Rhytismataceae, Rhytismatales).
Discussion
Morphology based diversity studies
Anamorphic fungi were dominant and accounted for approximately 80% of the total fungal taxa. This is largely consistent with previous morphological studies where anamorphic fungi constituted the bulk of fungal communities within pine needles25. This abundance of anamorphic taxa indicates that the cultural techniques tend to favour faster growing anamorphic fungi, which outcompete others. Fungal taxa with the highest frequency of occurrence were those that failed to produce any spores. They were categorised as mycelia sterila and not sequenced herein as these endophytes are rather well documented (e.g.5,6,30) and the main aim in this study was to use PCR-DGGE to investigate if pine needles harbour further potential unknown fungal taxa. Common saprobes, such as Cladosporium cladosporioides and Trichoderma sp., were also high in occurrence. These two taxa have previously been reported to be dominant and pioneer decomposers of Pinus4,30,31,32,33.
In this study, Pinus elliottii was found to have the highest fungal species richness among the three-tree species. This figure is much lower than other plant substrates such as grasses and palms but results corroborate other diversity studies on pine needles (e.g.7). Eighty fungal taxa were identified from two Pandanus species (Pandanus furcatus and P. tectorius) of the family Pandanaceae34 and 205 fungal species were identified from Gramineae and Cyperaceae35. Tokumasu36 reported that the number of fungal species recorded from the L-type needles of Pinus densiflora, P. sylvestris, P. banksiana, P. rigida and P. strobus were 26, 12, 20, 26 and 22 respectively. Pine needle morphology and structure may account for this lower diversity. The presence of thick waxy cuticles covering the needles is a major characteristic of the Pinaceae. Penetration of cuticle is the crucial step for fungal invasion of the leaves and this waxy cuticle is a primary defense against microbial invasion37.
Keteleeria fortunei, Pinus elliottii and P. massoniana yielded 26, 31 and 30 taxa, respectively. Twenty taxa were common on all tree species. For Keteleeria fortunei, species richness in F-type needle was lower than L-type needles (19 as compared to 25) and this supports other studies that fungal diversity is usually higher in L-type needles36. One explanation to account for this difference is that fungi in L-type needles are mostly those occurring at leaf senescence (possibly endophytes). Another explanation might be related to the amount of sugars present within pine needles. Millar38 suggested that there are inadequate sugar residues within F-type needles as compared to L-type and therefore less nutrient availability. This phenomenon has already been demonstrated in numerous other studies (e.g.33,36,39). A similar species richness between Pinus elliotti and P. massoniana can be explained by the fact that both hosts are related, have similar needle morphology and are from the same location. Fungal composition of needle litter within Pinus species have also been reported to be similar32,33,36. Lophodermium pinastri (Leotiomycetes) was only found in Pinus elliottii and P. massoniana but not in Keteleeria fortunei. It is possible that this fungal species, as stated in12, is restricted to Pinus but no other tree species in the family Pinaceae. Our phylogenetic data based on uncultured sequences, however, seem to contradict what Minter12 postulated. Based on sequence analyses in this study, it cannot be confirmed that we recovered L. pinastri but at least four phylotypes from K. fortunei from fresh and L type needles that were related to the Leotiomycetes. Our molecular data clearly demonstrate that diversity of Letiomycetous fungi is higher than expected. Phylogenies also indicated that phylotypes from the same host or needle types (either fresh or decayed) were not necessarily related.
Molecular Based Diversity Studies
Phylotype diversity and community shift
DGGE fingerprints suggest that fungal diversity was high in the living (fresh) and during early needle decay (L type) whereas only 1–2 bands could be detected in the F type needles. This suggests that only 1–2 fungal phylotypes are dominant during that particular decay stage (F type). A higher phylotypic diversity was also observed in Keteleeria fortunei (Figs 1 and 2). In most cases, we found that more bands (with clear signals and more conspicuous) were recovered from the L-type needles, except in Pinus ellioti. DGGE banding profiles also revealed a diversity shift of fungal 18S rDNA from fresh to L and F decay stages. This obviously indicates that the fungal community was initially represented by more fungi in the living and L type needles and then diversity of the fungal community gradually decreases as the needles start to decay. As further decomposition of pine needles occurs (F-type needles), only 1–2 abundant fungi dominate and this indicates a community shift from one that is highly diverse to one with a reduced number of abundant fungi (e.g.25). The latter are presumably those that have greater role in the decomposition process but such an assumption needs to be verified. It should be noted, however, that such a diversity shift and distinct phylotypes associated with different decay stages have not been found in the other two Pinus species investigated. Based on available information, it is very difficult to account for this difference. It might be that (i) phenolic compounds found in large amounts in Pinus needles could have been inhibitory to DNA and PCR amplication40 or (ii) limitations in our molecular approach (e.g. DGGE resolution) that failed to recover other taxa41,42.
Phylogeny of phylotypes and comparison with phenotypic diversity
To improve our understanding of the relationship between diversity and taxonomy, we generated phylogenies and made an attempt to identify and classify the uncultured phylotypes. Phylotypes 17 and G isolated from fresh and L type needles respectively cluster together with high support, and are fungi related and belong to Penicillium. The latter are common asexual fungi that have also been obtained from cultural studies in high frequency (Table 1) and similar findings have also been reported9,43. Phylotype P belongs to the order Hypocreales and could possibly be a Geosmithia species given their close phylogenetic affinity. Jankowiak et al.44 reported nine unnamed species of Geosmithia that are associated with bark beetle species infecting pine branches while Kolařík et al.45 reported Geosmithia as either saprobes or pathogens. Our phylotype herein is presumably a saprobe as it was collected from F type needles but the possibly that it is an endophyte and behaves as a pathogen cannot be ruled out as endophytes are known to adopt different ecological strategies (e.g.2). Taxa from the Eurotiales and Hypocreales have been recovered herein given their high abundance, ubiquitous nature and are also commonly recovered from microbial communities associated with healthy and diseased spruce seedlings22.
A common endophytic genus from pine needles is Lophodermium, a leotiomycetous fungus of the order Rhytismatales6,32,36,46. Based on previous studies, we have hypothesised that Lophodermium would be the predominant taxon but such was not the case. Surprisingly only 5 phylotypes (10 KF, C15 KF, M31 KF, J29 KF and L32 KF) constitute an independent lineage and share a close relatedness to Lophodermium pinastri, Meria laricis, Botryotinia fuckelina and Mollisia cinerea, which are members of the Rhytismatales and Helotiales (Leotiomycetes, Fig. 3). Even whether these phylotypes are Lophodermium species is doubtful and our results corroborates those of Ono et al.47 who reported few OTU’s isolates of the Leotiomycetes from pines. Numerous endophytic studies on pine needles based on cultural methods have recovered Lophodermium species (e.g5). It has generally been assumed that Lophodermium species, given their abundance, play an important and primary role in needle decomposition. However, are Lophodermium species pioneer decomposers? It could be possible that the large number of phylotypes in this study could also be a significant contributor to decomposition. However, function of most of them is still elusive and it is impossible to draw real conclusions as to their role. It is highly possible that they perish soon after needle senescence, as some of them were not recovered from cultural and morphological studies. For those similar taxa that were recovered from fresh and decayed needles, the assumption that they change their lifestyle (i.e. change from endophytic mode to saprophytic mode) to enable them to make maximum use of available nutrients or to flourish in other ecological niches cannot be overlooked. Such a competitive strategy has already been demonstrated in leaves of Magnolia lillifera2.
Contrary to our expectation, DNA sequence data analysed demonstrate that the majority of phylotypes from living leaves and L type needles (A10 PM-X25 PM) constitute two major clades close to other bitunicate taxa (Dothideomycetes) especially to Ochroconis, Scolecobasidium, and Ascorhombispora, which belongs to the family Sympoventuriaceae (Venturiales). Could these phylotypes represent isolated lineages that can provide clues to unaccounted fungal diversity and potential ecological roles? Phylogeny indicates that this novel monophyletic lineage has accumulated more nucleotide substitution than other bitunicate taxa and it might be possible that these uncultured taxa diverged recently. It is, for the time being, difficult to propose any taxonomic rank (Class, Order, or Family) to accommodate this new lineage, given that the systematics of its sister taxa are unclear and unresolved.
Another interesting finding is the detection of a particular phylotype (H KF) that has a close phylogenetic affiliation to Lecanora hybocarpa, Flavocetraria nivalis, Metus conglomeratus and Pilophorus species (order Lecanorales). To date, there are no reports of any lichenised fungi closely related to those genera from pine needles. This clearly shows that DGGE analyses allow (i) recovery of other leotiomycetous and lichenised fungi that possibly previously escaped morphological detection and (ii) suggest that diversity among these groups of fungi within pine needles is higher than expected. Many taxa of the Pleosporales appear to be quite common from environmental samples6,21,22. DGGE analyses on community DNA from rhizosphere of diseased seedling of Picea mariana (Black spruce) recovered three environmental clones that were related to Paraphaeosphaeria and Leptosphaeria (Pleosporales)22. Another uncultured clone related to Phoma herbarum (Pleosporales) was isolated from cankers of Spruce seedling. Pleosporales fungi are economically important plant pathogens and there is a high probability therefore that these phylotypes are predominantly endophytes and colonise pine needles without causing any apparent disease. Our molecular data corroborate with those obtained from6 who recovered eight environmental clones from needles of Pinus taeda that were related to Alternaria (Pleosporales). Our molecular approach recovered three phylotypes that fits within the Pleosporaceae and related to Cochlobolius species which are closely related to Pleospora herbarum. It is worth pointing out that our morphological studies also yielded identified taxa such as Curvularia species (anamorphs of Cochliobolus), Embellisia and Phoma (anamorphic Alternaria) but it is difficult to ascertain whether the DGGE phylotypes in this study were any of these taxa (but it cannot be ruled out). Another major discrepancy noted was the absence of Cladosporium from our molecular approach but this was abundant from our morphological survey as well as from environmental clones6. In addition, no basidiomycetes and xylariaceous taxa were recovered here but these were common from Pinus taeda6.
Significance and limitations of the study
Our molecular approach therefore, indicates that the majority of taxa are not necessarily related to those isolated from conventional techniques. Arnold et al.6 also reported that a number of operational taxonomic units recovered were distinct from those obtained from the cultural studies. We presume that those “DNA based taxa” might play a more active ecological and functional role in the community, whereas those that are easily identified from the cultural and morphological approach represent only selective endophytes- those that colonise aerial parts or inner surfaces and those that are fast growing on selective media but not necessarily abundant in the tissues. Apart from inadequate knowledge about the precise ecological nature and function of these phylotypes, this study raises a few questions about the reliability of DGGE and only one molecular approach in documenting the diversity of natural fungal consortia. Although primers targeting partial ribosomal DNA (this and other studies) appear to detect fungal DNA from environmental samples, they usually amplify only a short fragment that is either too conserved (e.g. 18S rDNA) or too variable (e.g ITS rDNA). Only one primer pair which recovers about 350 bp amplicons was tested and it appears to be taxonomically biased and too specific as no basidiomycetes or other fungal groups (apart from Ascomycetes) were recovered. This study could have therefore reflected only a fraction of those cryptic taxa within pine needles. The small number of bands recovered could possibly be linked to (i) the limited resolution of DGGE (ii) individual bands that were masked and possibly contain more than one phylotypes, and (iii) specificity of primers and size of DNA fragment amplified. In bacterial diversity studies, it has been shown that different DNA protocols and purification methods yield different DGGE profiles48. Similar research are required in fungal studies in addition to the use of a combination of different primer pairs that would target a wider group of fungi. Another intriguing finding was the abundance of taxa that constitutes new phylogenetic lineages. Given that only DNA sequences or other genomic data are available for uncultured “taxa”, an appropriate nomenclatural system and “species” definition based on DNA sequence characters are also critical issues that need to be dealt with. Potential bias of rDNA primers and application of DGGE in estimating fungal diversity have already been detailed49 and are not discussed here.
Future avenues with regards to environmental sequences
Phylogenetic information obtained here cannot be linked to function but it does provide insights into the relative abundance of uncultured phylotypes or “taxa” that could be more functional and metabolically active than those obtained from morphological and cultural methods. It is obvious that fungal diversity knowledge is still rudimentary. Our study reveals that even small substrates, such as pine needles harbour a diversity of enigmatic and cryptic operational taxonomic units (potential novel species) that can be recovered only by a molecular approach. We suspect that there could be a myriad of more cryptic taxa. Given that the roles of these phylotypes are still enigmatic, we should not only rely on findings based on taxonomic characterization but the way forward would be to proceed with metagenomic sequencing to enable analyses of functional gene content and properly understand community structure. Future work should also target to look into possible avenues about how to induce those taxa to sporulate in culture, describe their morphology, assess their physiological characteristics and exploit their potential ecological and functional roles. While our molecular approach herein permits recovery of potential new phylotypes, we are also fully cognizant of the intricacies associated with reliability of the 18S rDNA region and possibly more gene regions with greater sequence variability (e.g. ITS) and their drawbacks should be investigated to provide a more realistic picture of mycobiota. Should we consider those “loner” sequences/environmental phylotypes with distinct lineages as neglected properties or orphans simply because we cannot give them scientific names? There is a need to translate those molecular sequences into formal taxonomic names and consider them nomenclaturally valid to facilitate taxonomy, as postulated by Hibbett50 and Jeewon et al.51 but with caution52.
Experimental Procedures
Morphological based methods
Collection sites and substrates
Two study sites in Hong Kong (Cape D’ Aguilar and Sai Kung Yung Shue O) and senescent needles were collected randomly from the ground from three hosts (Keteleeria fortunei, and Pinus elliottii and P. massoniana) at 2 monthly intervals over a year. Decaying needles were categorised as L-types (senescent needles recently fallen light brown and occurring in an uncompacted layer at the surface) and F-types (grey to black needles, with softened tissues, low tensile strength and a higher moisture content than L-type needles forming a more closely packed stratum) as described by Kendrick and Burges4.
Plating technique
To document fungi, a modified washing and plating technique of Tokumasu36,53 was used for the treatment of needles. Five needles were placed into centrifuge tubes with the addition of sterilized 0.005% Aerozol OT solution (Di-iso-octyl sodium sulfosuccinate). The capped tubes were shaken vigorously in a vortical shaker at 300 oscillations per minute for one minute. Solutions were then discarded and the procedure repeated five times. Needles were then washed with sterile water three times and dried on sterilised filter paper for one day (to prevent bacterial growth). Needles were then placed onto the surface of half-strength cornmeal agar plates (only one needle in one agar plate) and incubated at room temperature (~24 °C). The samples were examined weekly over 1–2 months to check for any developing sporulating structures under the microscope. Fungal species diversity one each tree species was calculated using the Shannon Diversity index (H1).
Molecular Based Methods
DNA extraction and amplification
Living needles (green needles still attached to branches), L-type needles and F-type needles of Keteleeria fortunei and Living needles and L-type needles of Pinus elliottii and Pinus massoniana were also collected from Cape D’ Aguilar. All the samples (at least 3 samples from each needle type) were kept in paper bags and stored in the fridge at 4 °C and worked on the following days. Every single needle was cut into segments and ground with liquid nitrogen (with 1% PVPP) and DNA extracted following Jeewon et al.54 but using hot phenol which yielded better DNA suitable for amplification. The 5′end of the 18S ribosomal DNA (rDNA) was amplified with the primer pairs NS1 (5′-GTA GTC ATA TGC TTG TCT C-3′) and fungal-specific oligonucleotide GCfung (5′-CGC CCG CCG CGC CCC GCG CCC GGC CCG CCG CCC CCG CCC CAT TCC CCG TTA CCC GTT G-3′) as suggested by May et al.55. One PCR reaction mixture included 5 μl of 10x PCR buffer, 4 μl of 10 μM dNTPs, 0.5 μl of 0.1% bovine serum albumin (BSA), 0.5 μl of 1% PVP, 0.3 μl Taq polymerase, 31.7 μl nanopure water, 1.5 μl of 10 μM of each primer, 3 μl of 30 ng of each DNA sample. PCR profile was as follows: Initial denaturation for 3 min at 95 °C followed by 35 amplification cycles at 95 °C for 1 min, 50 °C for 1 min and 72 °C for 3 min, followed by a final extension of 72 °C for 10 min. PCR products were used for DGGE analysis.
DGGE analysis
An aliquot of 25 μl of each PCR product was mixed with and equal amount of dye and loaded in urea/formamide (100% denaturing solution with 7M urea and 40% formamide) denaturing gradient in a 7% acrylamide gel and run at 150 V in 1x TAE buffer (pH 8) at 60 °C (DGGE-2001, CBS Scientific Co.). The denaturing gradient was first optimised using a range of denaturant concentrations and electrophoresis with time. Gradient of 5–55% urea/formamide was selected based on a series of optimisation experiments. The gel was stained in 1x TAE containing ethidium bromide (EtBr) for 30 min and de-stained in 1x TAE for 20 min. The gel images were captured under UV trans-illumination and analysed by the computer program Gel Doc. Bands were excised and each fragment was placed into a 1.5 ml marked microcentrifuge tubes. They were reamplified with the original primers and purified using a DNA purification kit (Amersham Biosciences GFXTM PCR DNA and Gel Band Purification Kit) prior to automated sequencing using NS1 primer (API Prism 3730, Applied Biosystems). DGGE analyses were performed twice under the same conditions.
Phylogenetic analyses
Partial 18S rDNA gene sequences recovered from DGGE were aligned in Bioedit56, in MAFFT v. 7 with the web server (http://mafft.cbrc.jp/alignment/server) with similar sequences from different fungal orders57,58,59,60. Alignments were also checked and edited manually to optimise alignment. Phylogenetic analyses were conducted in PAUP version 4.0b10 as described by61,62,63. In addition to Maximum Parsimony, alternative analyses based on the maximum likelihood (ML) and Bayesian analyses (BI) were also conducted (as described in64,65,66). Phylograms were visualised with FigTree v. 1.4.0 and annotated in Microsoft PowerPoint (2007) as described in65.
Change history
15 July 2020
An amendment to this paper has been published and can be accessed via a link at the top of the paper.
References
Promputtha, I. et al. Ribosomal DNA fingerprinting in the identification of non-sporulating endophytes from Magnolia liliifera (Magnoliaceae). Fungal Diversity 20, 167–186 (2005).
Promputtha, I. et al. A phylogenetic evaluation of whether endophytes become saprotrophs at host senescence. Microbial Ecology 53, 579–590 (2007).
Jeewon, R. et al. DNA based identification and phylogenetic characterisation of endophytic and saprobic fungi from Antidesma madagascariense, a medicinal plant in Mauritius. Journal of Mycology 2013, 1–13 (2013).
Kendrick, W. B. & Burges, A. Biological aspects of the decay of Pinus sylvestris leaf litter. Nova Hedwigia 4, 313–342 (1962).
Ganley, R. J., Brunsfeld, S. J. & Newcombe, G. A community of unknown, endophytic fungi in western white pine. Proceedings of the National Academy of Sciences of the United States of America 101, 10107–10112 (2004).
Arnold, A. E. et al. Diversity and phylogenetic affinities of foliar fungal endophytes in loblolly pine inferred by culturing and environmental PCR. Mycologia 99, 185–206 (2007).
Zamora, P., Martínez-Ruiz, C. & Diez, J. J. Fungi in needles and twigs of pine plantations from northern Spain. Fungal Diversity 30, 171–184 (2008).
Terhonen, E. et al. The effect of latitude, season and needle-age on the mycota of Scots pine (Pinus sylvestris) in Finland. Silva Fennica 45, 301–317 (2011).
Mahajan, R. et al. Autochthonous microbial community associated with pine needle forest litterfall influences its degradation under natural environmental conditions. Environmental monitoring and assessment 188, 1–10 (2016).
Nguyen, D. et al. Do foliar fungal communities of Norway spruce shift along a tree species diversity gradient in mature European forests? Fungal Ecology 23, 97–108 (2016).
Muyzer, G. et al. Profiling of complex microbial populations by denaturing gradient gel electrophoresis analyses of polymerase chain reaction-amplified genes coding for 16S rDNA. Applied and Environmental Microbiology 59, 695–700 (1993).
Minter, D. W. Lophodermium on pines. Mycological Papers 147, 1–54 (1981).
Anderson, I. C., Campbell, C. D. & Prosser, J. I. Potential bias of fungal 18S rDNA and internal transcribed spacer polymerase chain reaction primers for estimating fungal biodiversity in soil. Environmental Microbiology 5, 36–47 (2003).
Smit, E. et al. Analyses of fungal diversity in the wheat rhizosphere by sequencing of cloned PCR-amplified genes encoding 18S rRNA and temperature gradient gel electrophoresis. Applied and Environmental Microbiology 65, 2614–2621 (1999).
Vainio, E. J. & Hantula, J. Direct analysis of wood-inhabiting fungi using denaturing gradient gel electrophoresis of amplified ribosomal DNA. Mycological Research 104, 927–936 (2000).
Duong, L. M. et al. DGGE coupled with ribosomal DNA phylogenies reveal uncharacterized fungal phylotypes on living leaves of Magnolia liliifera. Fungal Diversity 23, 121–138 (2006).
Xiao, C. et al. Effects of cultivation ages and modes on microbial diversity in the rhizosphere soil of Panax ginseng. Journal of ginseng research 40, 28–37 (2016).
Jie, W., Liu, X. & Cai, B. Diversity of Rhizosphere Soil Arbuscular Mycorrhizal Fungi in Various Soybean Cultivars under Different Continuous Cropping Regimes. PLoS One 8, e72898 (2013).
Nikolcheva, L. G., Cockshutt, A. M. & Bärlocher, F. Determining diversity of freshwater fungi on decaying leaves: comparison of traditional and molecular approaches. Applied and Environmental Microbiology 69, 2548–2554 (2003).
Nikolcheva, L. G. & Barlocher, F. Seasonal and substrate preferences of fungi colonizing leaves in streams: traditional vs. molecular evidence. Environmental Microbiology 7, 270–280 (2005).
Cullings, K., Raleigh, C., New, M. H. & Henson, J. Effects of artificial defoliation of pines on the structure and physiology of the soil fungal community of a mixed pine-spruce forest. Applied and environmental microbiology. 71, 1996–2000 (2005).
Filion, M., Hamelin, R. C., Bernier, L. & St-Arnaud, M. Molecular profiling of rhizosphere microbial communities associated with healthy and diseased black spruce (Picea mariana) seedlings grown in a nursery. Applied and environmental microbiology 70, 3541–3551 (2004).
Yeung, Q. S. Y., Jeewon, R. & Hyde, K. D. Cancellidium pinicola sp. nov. from Pinus massoniana and its phylogeny. Cryptogamie Mycologie 27, 295–304 (2006).
Jeewon, R., Yeung, Q. S. Y. & Hyde, K. D. Novel phylogenetic group within Thozetella (Chaetosphaeriaceae): a new taxon based on morphology and DNA sequence analyses. Canadian Journal of Microbiology 55, 680–7 (2009).
Aneja, M. K. et al. Microbial Colonization of Beech and Spruce Litter—Influence of Decomposition Site and Plant Litter Species on the Diversity of Microbial Community. Microbial Ecology 52, 127–135 (2006).
Lacap, D. C. et al. An evaluation of the fungal ‘morphotype’ concept based on ribosomal DNA sequences. Fungal Diversity 12, 53–66 (2003).
MacNaughton, S. J. et al. Microbial population changes during bioremediation of an experimental oil spill. Applied and Environmental Microbiology 65, 3566–3574 (1999).
Boon, N. et al. Evaluation of nested PCR-DGGE (denaturing gradient gel electrophoresis) with group-specific 16S rRNA primers for the analysis of bacterial communities from diferent wastewater treatment plants. FEMS Microbiology Ecology 39, 101–112 (2001).
Sekiguchi, H. et al. A single band does not always represent single bacterial strains in denaturing gradient gel electrophoresis analysis. Biotechnology Letters 23, 1205–1208 (2001).
Osono, T. & Takeda, H. Comparison of litter decomposing ability among diverse fungi in a cool temperate deciduous forest in Japan. Mycologia 94, 421–427 (2002).
Paul, N. C. & Yu, S. H. Two Species of Endophytic Cladosporium in Pine Trees in Korea. Mycobiology 36, 211–216 (2008).
Legault, D. & Dessureault, M. Mycoflora of Pinus banksiana and Pinus resinosa needles. II. Epiphytic fungi. Canadian Journal of Botany 67, 2061–2065 (1988).
Manoch, L., Tokumasu, S. & Tubaki, K. A. preliminary survey of microfungal flora of pine leaf litter in Thailand. Transactions of the Mycological Society of Japan 27, 159–165 (1986).
Whitton, S. R. Microfungi on the Pandanaceae. PhD thesis. University of Hong Kong, Hong Kong (1999).
Wong, M. K. M. & Hyde, K. D. Diversity of fungi on six species of Gramineae and one species of Cyperaceae in Hong Kong. Mycological Research 105, 1485–1491 (2001).
Tokumasu, S. Leaf litter fungi of the forests of Pinus densiflora and four introduced pines at Sugadaira, central Japan. Transactions of the Mycological Society of Japan 19, 383–390 (1978).
Gutschick, V. P. Research reviews biotic and abiotic consequences of differences in leaf structure. New Phytologist 143, 3–18 (1999).
Millar, C. S. Decomposition of coniferous leaf litter. In: Biology of Plant Litter Decomposition 1 (eds C. H. Dickinson and G. J. F., Pugh). Academic Press, London, England: 105–128 (1974).
Tokumasu, S. Microfungal flora on decaying pine needles collected from Kitaiojima, a subtropical island in the Pacific. Transactions of the Mycological Society of Japan 26, 481–486 (1985).
Smart, M. & Roden, L. C. A small-scale RNA isolation protocol useful for high-throughput extractions from recalcitrant plants. South African Journal of Botany 76, 375–379 (2010).
Jeewon, R. & Hyde, K. D. Diversity and detection of fungi from environmental samples: traditional versus molecular approaches. Advanced techniques in soil microbiology. Soil Biology Series. Springer, Berlin, 1–15 (2007).
Al-Mailem, D. M., Mayada, M. K. & Radwan, S. S. Capabilities and limitations of DGGE for the analysis of hydrocarbonoclastic prokaryotic communities directly in environmental samples. MicrobiologyOpen 6(5) (2017).
Romeralo, C., Diez, J. J. & Santiago, N. F. Presence of fungi in Scots pine needles found to correlate with air quality as measured by bioindicators in northern. Forest Pathology 42, 443–453 (2012).
Jankowiak, R., Kolařík, M. & Bilański, P. Association of Geosmithia fungi (Ascomycota: Hypocreales) with pine- and spruce-infesting bark beetles in Poland. Fungal Ecology 11, 71–79 (2014).
Kolařík, M., Freeland, M., Utlet, C. & Tisserat, N. Geosmithia morbida s p. n ov., a new phytopathogenic species living in symbiosis with the walnut twig beetle (Pityophthorus juglandis) on Juglans in USA. Mycologia 103, 325–332 (2011).
Tokumasu, S., Aoki, T. & Oberwinker, F. Fungal succession on pine needles in Germany. Mycoscience 35, 29–37 (1994).
Ono, R. et al. A comparison of the community diversity of foliar fungal endophytes between seedling and adult loblolly pines (Pinus taeda). Fungal. Biology 119, 917–928 (2015).
Maarit-Niemi, R., Heiskanen, I., Wallenius, K. & Lindstrom, K. Extraction and purification of DNA in rhizosphere soil samples for PCR-DGGE analysis of bacterial consortia. Journal of microbiological methods 45, 155–165 (2001).
Jeewon R. & Hyde K. D. Diversity and detection of Fungi from environmental samples: Traditional versus Molecular approaches. In (Eds Varma, A. & Oelmuller, R.) Advanced Techniques in Soil Microbiology pp 1–11, (2007) Soil Biology Series. Springer-Verlag Press.
Hibbett, D. The invisible dimension of fungal diversity. Science 351, 1150–1151 (2016).
Jeewon, R. et al. Nomenclatural and identification pitfalls of endophytic mycota based on DNA sequence analyses of ribosomal and protein genes phylogenetic markers: A taxonomic dead end? Mycosphere 8, 1802–1817 (2017).
Hongsanang, S. et al. Can we use environmental DNA as holotypes? Fungal Diversity 91, 1–30 (2018).
Tokumasu, S. Fungal succession on pine needles fallen at different seasons: the succession of surface colonisers. Mycoscience 39, 417–423 (1998).
Jeewon, R., Liew, E. C. Y. & Hyde, K. D. Phylogenetic relationships of Pestalotiopsis and allied genera inferred from ribosomal DNA sequences and morphological characters. Molecular Phylogenetics and Evolution 25, 378–392 (2002).
May, L. A., Smiley, B. & Schmidt, M. G. Comparative denaturing gradient gel electrophoresis analyses of fungal communities associated with whole plant corn silage. Canadian Journal of Botany 47, 829–841 (2001).
Hall, T. A. BioEdit: a user-friendly biological sequence alignment editor and analyses program for Windows 95/98/NT. Nucleic Acids Symposium Series 41, 95–98 (1999).
Tibpromma, S. et al. Fungal diversity notes 491–603: taxonomic and phylogenetic contributions to fungal taxa. Fungal Diversity 83, 1–261 (2017).
Wang, H. K. et al. The polyphyletic nature of Pleosporales: An example from Massariosphaeria based on ribosomal DNA and RBP2 gene phylogenies. Fungal. Biology 111, 1268–1276 (2007).
Zhang, Y. et al. Multi-gene phylogeny and morphotaxonomy of Amniculicola lignicola: novel freshwater fungus from France and its relationships to the Pleosporales. Fungal. Biology 112, 1186–1194 (2008).
Thongkantha, S. et al. Molecular phylogeny of Magnaporthaceae (Sordariomycetes) with a new species, Ophioceras chiangdaoensis from Dracaena loureiri in Thailand. Fungal Diversity 34, 157–173 (2009).
Jeewon, R. et al. Phylogenetic significance of morphological characters in the taxonomy of Pestalotiopsis species. Molecular Phylogenetics and Evolution 27, 372–383 (2003).
Jeewon, R., Liew, E. C. Y. & Hyde, K. D. Phylogenetic evaluation of species nomenclature of Pestalotiopsis in relation to host association. Fungal Diversity 17, 39–55 (2004).
Shenoy, B. D. et al. Ribosomal and RPB2 DNA sequence analyses suggest that Sporidesmium and morphologically similar genera are polyphyletic. Mycological Research 110, 916–928 (2006).
Tang, A. M. C., Jeewon, R. & Hyde, K. D. Phylogenetic utility of protein (RPB2, b-tubulin) and ribosomal (18S, 28S) gene sequences in the systematics of Sordariomycetes (Ascomycota, Fungi). Antonie van Leeuwenhoek 91, 327–349 (2007).
Wanasinghe, D. N. et al. Fenestellaceae with descriptions of new Fenestella species and Neocucurbitaria gen. nov. Mycosphere 8, 397–414 (2017).
Liu, J. K. et al. Ranking higher taxa using divergence times: a case study in Dothideomycetes. Fungal Diversity 84, 75–99 (2017).
Acknowledgements
This study was funded by the Hong Kong Research Grants Council (HKU 7322/04 M and CRCG 10205773). The University of Hong Kong is acknowledged for providing Yeung with a postgraduate scholarship to pursue her Master’s degree and University of Mauritius and Zheijiang University are acknowledged for research support.
Author information
Authors and Affiliations
Contributions
R.J. designed the project, contributed to lab work, DNA analyses, paper writing and supervised the project. Q.S.Y. Yeung was involved in collection of samples as well as performed lab work icw morphology and DNA samples. D.N.W. contributed to DNA analyses and D.N.W., S.R., D.P., H.K.W. and K.D.H. reviewed and assisted in manuscript preparation and review.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Jeewon, R., Yeung, Q.S.Y., Wannasinghe, D.N. et al. Hidden mycota of pine needles: Molecular signatures from PCR-DGGE and Ribosomal DNA phylogenetic characterization of novel phylotypes. Sci Rep 8, 18053 (2018). https://doi.org/10.1038/s41598-018-36573-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-018-36573-z
This article is cited by
-
Harmoniella junipericola (Strobiloscyphaceae, Pezizales), a new hyphomycete inhabiting on leaf litter of Juniperus chinensis in Japan, its molecular phylogenetic classification, and brief notes on other Harmoniella species
Mycological Progress (2023)
-
Foliar endophytic fungi: diversity in species and functions in forest ecosystems
Symbiosis (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.