Research Article |
Corresponding author: Georgios Zervakis ( zervakis@aua.gr ) Academic editor: María P. Martín
© 2020 Vassiliki Fryssouli, Georgios Zervakis, Elias Polemis, Milton A. Typas.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Fryssouli V, Zervakis GI, Polemis E, Typas MA (2020) A global meta-analysis of ITS rDNA sequences from material belonging to the genus Ganoderma (Basidiomycota, Polyporales) including new data from selected taxa. MycoKeys 75: 71-143. https://doi.org/10.3897/mycokeys.75.59872
|
Ganoderma P. Karst. is a cosmopolitan genus of white-rot fungi which comprises species with highly-prized pharmaceutical properties, valuable biotechnological applications and of significant phytopathological interest. However, the status of the taxonomy within the genus is still highly controversial and ambiguous despite the progress made through molecular approaches. A metadata analysis of 3908 nuclear ribosomal internal transcribed spacer (ITS) rDNA sequences obtained from GenBank/ENA/DDBJ and UNITE was performed by targeting sequences annotated as Ganoderma, but also sequences from environmental samples and from material examined for the first time. Ganoderma taxa segregated into five main lineages (Clades A to E). Clade A corresponds to the core of laccate species and includes G. shanxiense and three major well-supported clusters: Cluster A.1 (‘G. lucidum sensu lato’) consists of taxa from Eurasia and North America, Cluster A.2 of material with worldwide occurrence including G. resinaceum and Cluster A.3 is composed of species originating from all continents except Europe and comprises G. lingzhi. Clade B includes G. applanatum and allied species with a Holarctic distribution. Clade C comprises taxa from Asia and Africa only. Clade D consists of laccate taxa with tropical/subtropical occurrence, while clade E harbours the highest number of non-laccate species with a cosmopolitan distribution. The 92 Ganoderma-associated names, initially used for sequences labelling, correspond to at least 80 taxa. Amongst them, 21 constitute putatively new phylospecies after our application of criteria relevant to the robustness/support of the terminal clades, intra- and interspecific genetic divergence and available biogeographic data. Moreover, several other groups or individual sequences seem to represent distinct taxonomic entities and merit further investigation. A particularly large number of the public sequences was revealed to be insufficiently and/or incorrectly identified, for example, 87% and 78% of entries labelled as G. australe and G. lucidum, respectively. In general, ITS demonstrated high efficacy in resolving relationships amongst most of the Ganoderma taxa; however, it was not equally useful at elucidating species barriers across the entire genus and such cases are outlined. Furthermore, we draw conclusions on biogeography by evaluating species occurrence on a global scale in conjunction with phylogenetic structure/patterns. The sequence variability assessed in ITS spacers could be further exploited for diagnostic purposes.
Biogeography, fungal diversity, ITS, medicinal mushroom, phylogeny, taxonomy
The genus Ganoderma P. Karst. (Basidiomycota, Polyporales) is characterised by a cosmopolitan distribution and high species diversity, especially in the tropics and subtropics. It comprises white-rot fungi that possess an efficient ligninolytic mechanism which is exploited in various biotechnological applications (
Ganoderma was erected to include only one species, G. lucidum (Curtis: Fr.) P. Karst. (
Still, the species concept in Ganoderma is not universally accepted and remains inadequately established (
The internal transcribed spacer (ITS: ITS1-5.8S-ITS2) region of the nuclear ribosomal RNA was erected as the formal DNA barcode in Fungi since it demonstrates a clear barcoding gap for a wide range of lineages and is often in good agreement with morphological/biological species concepts and could, therefore, be exploited for identification purposes (
Although the ITS region has been used in more Ganoderma studies than any other marker for proposing new taxa (
On the basis of the discrepancies and shortcomings noted before, the objectives of this study were: (i) to perform a thorough metadata analysis on the basis of a global dataset of Ganoderma ITS sequences, (ii) to evaluate the accuracy of specimen identifications to species, (iii) to determine not fully assessed (i.e. “Ganoderma sp.”) or erroneously labelled sequences in GenBank (www.ncbi.nlm.nih.gov/genbank/) and other relevant databases, in order to associate taxonomic names with phylogenetic lineages, (iv) to expand the knowledge on distribution and biogeography of Ganoderma species, (v) to examine controversial boundaries amongst existing species and complexes and (vi) to contribute to the development of quick and efficient sequence-based tools suitable for identification of Ganoderma species through the large-scale assessment of molecular information existing in public databases.
5.8S: ribosomal DNA 5.8S gene; ASV: amplicon sequence variant; BI: Bayesian Inference phylogenetic analysis; BPP: Bayesian Posterior Probability; BS: Bootstrap Support; DDBJ: DNA Data Bank of Japan; DOI: digital object identifier; DS: dataset; ENA: European Nucleotide Archive, ITS: ribosomal DNA internal transcribed spacer region; ML: Maximum Likelihood phylogenetic analysis; NCBI: National Center for Biotechnology Information.
Dried specimens were obtained on loan from the fungaria of the Bulgarian Academy of Sciences (
We studied 54 specimens in the form of either dried material or pure cultures. They represented well-established Ganoderma taxa with European distribution, i.e. four laccate taxa [G. lucidum (11 specimens), G. carnosum (5), G. resinaceum (18) and G. pfeifferi (3)] and two non-laccate/dull taxa [G. adspersum (8) and G. applanatum (5)]. In addition, material not fully identified (Ganoderma sp.; 2), together with commercial strains labelled as G. lucidum (1) and G. tsugae (1), was examined. Initial species labelling was in accordance with the identification made by the respective collector; however, at the end of our study, some of them were re-assessed. Details on their identity appear in Suppl. material
Mycelia for DNA extraction were produced in static potato dextrose (Difco, USA) cultures. Following a 10–15 day incubation period at 25 °C, the mycelia were harvested by filtration and either directly processed for DNA extraction or stored at -20 °C. Mycelium or dried basidiome samples were pulverised by a micropestle in the presence of sterile sand and liquid nitrogen. Total genomic DNA was subsequently extracted through the silica Plant II DNA Extraction Miniprep System (Macherey and Nagel, Germany) by following the standard CTAB protocol provided by the manufacturer with minor modifications, i.e. the lysis step was extended to 1 h at 65 °C and the precipitation step to 1 h at room temperature, while the final elution step was performed at 65 °C for 1 h (
Sequences of the ITS region were generated for phylogenetic analyses. DNA samples were subjected to PCR amplification of the ITS region by using the primer pairs ITS1/ITS4 (
Amplified fragments were examined by electrophoresis on 1% agarose gels. PCR products of the expected size were purified by microcentrifugation using the PureLink PCR purification kit (Invitrogen/Thermo Fisher Scientific, USA), according to manufacturer’s protocol. Purified amplicons were processed for bidirectional Sanger sequencing at CEMIA (University of Thessaly, Greece; https://cemia.eu/) using the forward ITS1 and reverse ITS4 primers. The resulting chromatograms were proofread and assembled using Unipro UGENE v.31 (
An initial dataset was compiled by retrieving and examining all publicly available ITS sequences assigned to the genus Ganoderma (followed either by a species epithet or not fully identified and labelled as “Ganoderma sp.”), as well as associated environmental samples and misidentified entries, appearing in GenBank/ENA/DDBJ (The International Nucleotide Sequence Database Collaboration;
Summarised information on the Ganoderma ITS sequences used in this study. Species marked in bold appear as they are presented in Clades and/or Clusters according to the outcome of the phylogenetic analysis (Fig.
CLADES/Species (no. of sequences per taxon) | Sequences original labelling (no. of sequences per taxon) | Geographic origin of sequenced material | Host type | BS/BPP values |
---|---|---|---|---|
CLADE A | 69%/- | |||
G. shanxiense (2) | G. shanxiense (2) | China | AD | 100%/1.00 |
CLADE A, Cluster A.1 | 100%/1.00 | |||
G. tsugae (57) | G. tsugae (46)*, G. lucidum (3), uncultured Ganoderma (3), G. ahmadii (1)*, G. carnosum (1), G. valesiacum (1), Ganoderma sp. (1), Polyporus tsugae (1) | Canada, Germany, India, Pakistan, UK, USA | AD, GS | |
G. oregonense (27) | G. oregonense (15), G. carnosum (5), G. tsugae (4), G. lucidum (2), uncultured soil fungus (1) | Canada, Estonia, USA | GS | |
G. carnosum (26) | G. carnosum (22), G. lucidum (4) | Belgium, Czech Republic, France, Greece, Poland, Slovakia, Slovenia | GS | |
G. aff. carnosum (4) | uncultured soil fungus (2), G. carnosum (1), G. lucidum (1) | Estonia, UK | 93%/1.00 | |
G. lucidum (153) | G. lucidum (107), uncultured soil fungus (28), G. tsugae (12), G. oerstedii (3), G. oregonense (1), G. valesiacum (1), Ganoderma sp. (1) | Algeria, Argentina, Armenia, Belgium, Bulgaria, China, Czech Republic, Finland, France, Greece, India, Iran, Italy, Norway, Russia, Slovakia, Spain, South Korea, Sweden, Thailand, UK, USA, commercial | AD, GS | |
G. leucocontextum – G. weixiensis (33) | G. leucocontextum (24)*, G. lucidum (3), Ganoderma sp. (3), G. weixiensis (2)*, G. carnosum (1) | China (Tibet), Nepal, Pakistan | AD, GS | 81%/1.00 |
CLADE A, Cluster A.2 | 94%/1.00 | |||
G. austroafricanum (2) | G. austroafricanum (1)*, G. aff. austroafricanum (1) | South Africa | AD | 89%/- |
G. hoehnelianum (14) | G. hoehnelianum (9), Ganoderma sp. (5) | China, Gabon, Myanmar | 99%/1.00 | |
G. weberianum (12) | G. weberianum (9), G. microsporum (1)*, G. sichuanense (1), Ganoderma sp. (1) | Philippines, Taiwan | AD | 91%/0.99 |
G. sichuanense (19) | G. sichuanense (9)*, G. weberianum (4), G. lucidum (2), G. tenue (2), Ganoderma sp. (1), uncultured soil fungus (1) | Australia, China, India | AD | |
G. carocalcareum (13) | Ganoderma sp. (8), G. weberianum (3), G. carocalcareum (2)* | Gabon, Cameroon | AD | -/0.98 |
Ganoderma sp. A1 (17) | G. weberianum (17) | India | AD | 100%/1.00 |
G. aff. weberianum (5) | G. weberianum (2), G. cf. weberianum (1), G. subamboinense (1), Ganoderma sp. (1) | Brazil, China, India | AD | |
G. mexicanum (17) | G. mexicanum (6), G. subamboinense var. laevisporum (5)*, G. subamboinense (2), G. weberianum (2), G. sessiliforme (1), G. tuberculosum (1) | Argentina, Brazil, Martinique, Mexico, USA | AD | |
G. parvulum (23) | G. parvulum (10), G. weberianum (5), Ganoderma sp. (4), G. subamboinense var. laevisporum (2), G. stipitatum (1), G. subamboinense (1)* | Brazil, Colombia, Costa Rica, Cuba, French Guiana, Mexico, USA | AD | 77%/- |
Ganoderma sp. A2 (2) | G. resinaceum (1), Ganoderma sp. (1) | China | 100%/1.00 | |
G. resinaceum (131) | G. resinaceum (105), G. lucidum (8), Ganoderma sp. (8), G. pfeifferi (6), G. cf. resinaceum (2), Polyporales sp. (2) | Belgium, Bulgaria, China, Czech Republic, Egypt, France, Greece, India, Iran, Iraq, Italy, Netherlands, Poland, Slovakia, South Africa, South Korea, Tunisia, Turkey, UK | AD | |
Ganoderma sp. A3 (12) | G. resinaceum (5), G. cf. resinaceum (3), Ganoderma sp. (2), G. lucidum (1), uncultured Ganoderma (1) | Malaysia, Taiwan, commercial | AD | 99%/1.00 |
G. aff. sessile (4) | G. lucidum (4) | India, Turkey | AD | 68%/- |
G. aff. polychromum (10) | G. resinaceum (5), G. sessile (2), G. cf. sessile (1), G. platense (1), G. zonatum (1) | Argentina, USA | AD | 98%/1.00 |
G. polychromum (19) | G. polychromum (11), G. lucidum (6), G. sessile (2) | China, India, USA | AD | 93%/1.00 |
G. sessile (228) | G. sessile (134), G. resinaceum (60), Ganoderma sp. (15), G. lucidum (10), G. oregonense (2), G. boninense (1), G. lobatum (1), G. neojaponicum (1), G. polychromum (1), G. valesiacum (1), Hericium erinaceum (1), uncultured root-associated fungus (1) | Argentina, China, India, Japan, Russia, South Korea, Taiwan, USA, commercial | AD, AM | 78%/1.00 |
CLADE A, Cluster A.3 | 76%/1.00 | |||
G. concinnum (2) | G. chalceum (1), G. concinnum (1) | Brazil | ||
G. tuberculosum (37) | G. tuberculosum (28), Ganoderma sp. (6), Coriolopsis caperata (1), G. parvulum (1), G. resinaceum (1) | Brazil, Colombia, Cuba, Martinique, Mexico, Panama, USA | AD | 100%/1.00 |
Ganoderma sp. A4 (2) | G. lucidum (2) | Argentina | 99%/1.00 | |
G. wiiroense (14) | G. wiiroense (9)*, Ganoderma sp. (3), G. lucidum (2) | Ghana, India, Senegal | AD | 100%/1.00 |
G. flexipes (7) | G. flexipes (7) | China, Laos, Vietnam | AD, GS | 100%/1.00 |
Ganoderma sp. A5 (7) | G. multiplicatum (7) | China, Myanmar | AD | 100%/1.00 |
G. philippii (102) | G. pseudoferreum (75), G. philippii (15), Ganoderma sp. (9), G. australe (2), uncultured soil fungus (1) | China, Indonesia, Malaysia, Thailand | AD | 97%/1.00 |
G. lingzhi (615) | G. lingzhi (333)*, G. lucidum (206), Ganoderma sp. (37), G. sichuanense (27)*, G. tsugae (5), Amauroderma rugosum (1), G. boninense (1), G. calidophilum (1), G. cupreum (1), G. luteomarginatum (1), Haddowia longipes (1), Laccaria bicolor (1) | Bangladesh, China, India, Iran, Iraq, Japan, Laos, Malaysia, Myanmar, Nepal, South Korea, Thailand, commercial | AD, AM | 100%/1.00 |
G. curtisii (142) | G. curtisii (124), G. meredithae (11)*, G. lucidum (3), G. curtisii sp. meredithae (2), Ganoderma sp. (2) | Mexico, USA, commercial | AD | |
G. ravenelii (12) | G. ravenelii (6), G. curtisii (3), G. lucidum (2), uncultured fungus (1) | India, USA | AD, GS | 78%/1.00 |
G. multiplicatum (17) | G. multiplicatum (10), G. perzonatum (7) | Brazil, Colombia, Mexico | AD | 99%/1.00 |
G. destructans – G. dunense (43) | G. destructans (39)*, G. dunense (3)*, uncultured soil fungus (1) | Cameroon, South Africa | AD | |
G. mizoramense (3) | G. mizoramense (2)*, G. lucidum (1) | India | AD | |
G. steyaertanum (39) | G. steyaertanum (34), G. aff. steyaertanum (3), Ganoderma sp. (2) | Australia, Indonesia | AD | 89%/0.99 |
G. martinicense (49) | G. parvulum (24), G. martinicense (18)*, G. perzonatum (2), G. lucidum (1), G. oerstedii (1), G. tornatum (1), G. tuberculosum (1), Ganoderma sp. (1) | Argentina, Brazil, Colombia, Cuba, Martinique, Mexico, USA | AD | 93%/1.00 |
G. multipileum (243) | Ganoderma sp. (112), G. lucidum (105), G. multipileum (22), Agaricales sp. (1), G. leucocontextum (1), G. lingzhi (1), Polyporaceae sp. (1) | China, India, Nepal, Pakistan, Philippines, Taiwan, Thailand | AD, AM, GS | 74%/0.97 |
Ganoderma sp. A6 (15) | G. tropicum (15) | India | AD | 100%/1.00 |
G. tropicum (33) | G. tropicum (15)*, G. fornicatum (12), G. williamsianum (2), Vanderbylia fraxinea (2), Ganoderma sp. (1), uncultured soil fungus (1) | China, India, Laos, Taiwan, Thailand | AD | |
Ganoderma sp. A7 (3) | G. fornicatum (3) | Malaysia | 100%/1.00 | |
CLADE B | 96%/1.00 | |||
Ganoderma sp. B1 (4) | Ganoderma sp. (4) | China, USA | 100%/1.00 | |
Ganoderma sp. B2 (3) | G. applanatum (1), G. lingzhi (1), G. multipileum (1) | Nepal | 100%/1.00 | |
G. applanatum (424) | uncultured soil fungus (230), G. applanatum (119), G. lipsiense (21), Fungi (plant leaf) (15), uncultured Ganoderma (15), uncultured fungus (8), G. adspersum (5), G. applanatum cplx (3), Ganoderma sp. (2), fungal sp. (1), G. australe (1), G. cf. applanatum (1), G. lobatum (1), G. oregonense (1), Trametes sp. (1) | Antarctica, Armenia, Austria, Bulgaria, Canada, China, Czech Republic, Estonia, France, Germany, Greece, Hungary, India, Japan, Kyrgyzstan, Latvia, Lithuania, Netherlands, Poland, Russia, Slovakia, South Korea, Thailand, UK, USA, commercial | AD, AM, GS | 99%/1.00 |
CLADE C | 100%/1.00 | |||
CLADE C, Cluster C.1 | 99%/1.00 | |||
G. neojaponicum (10) | G. neojaponicum (7), G. calidophilum (2), Ganoderma sp. (1) | China, Laos, Myanmar, Taiwan | AD | 99%/1.00 |
CLADE C, Cluster C.2 | 100%/1.00 | |||
Ganoderma sp. C1 (2) | Ganoderma sp. (2) | Cameroon | AM | 82%/1.00 |
G. aridicola (7) | Ganoderma sp. (6), G. aridicola (1)* | Cameroon, South Africa | AD, AM, GS | |
Ganoderma sp. C2 (3) | Ganoderma sp. (3) | Cameroon | AD | 100%/1.00 |
G. enigmaticum – G. thailandicum (10) | G. enigmaticum (7)*, G. thailandicum (2)*, uncultured soil fungus (1) | Ghana, Ivory Coast, South Africa, Thailand | AD, GS | |
G. casuarinicola (63) | Ganoderma sp. (47), G. casuarinicola (6)*, G. enigmaticum (4), uncultured fungus (2), G. applanatum (1), G. carnosum (1), G. lucidum (1), uncultured Ganoderma (1) | China, India, Sri Lanka | AD, AM, GS | 71%/- |
CLADE D | ||||
CLADE D, Cluster D.1 | 98%/1.00 | |||
G. mbrekobenum (36) | Ganoderma sp. (21), G. mbrekobenum (11)*, G. applanatum (1), G. carnosum (1), G. lucidum (1), G. tsugae (1) | Ghana, India, Senegal, Sri Lanka | AD, AM, GS | 98%/1.00 |
CLADE D, Cluster D.2 | ||||
G. nasalanense (17) | G. australe (9), Ganoderma sp. (4), G. nasalanense (2)*, uncultured soil fungus (2) | India, Laos, Malaysia, Vietnam | AD | 98%/1.00 |
G. sinense (66) | G. sinense (45), Ganoderma sp. (8), G. lucidum (5), G. japonicum (4), G. subresinosum (2), G. atrum (1), G. formosanum (1) | China, Taiwan, Thailand | AD | 100%/1.00 |
CLADE D, Cluster D.3 | 75%/1.00 | |||
G. cupreum (8) | G. cupreum (4), G. australe (1), G. cf. cupreum (1), G. chalceum (1), uncultured fungus (1) | Cameroon, Gabon, Malaysia, South Africa, Tanzania | AD | 99%/1.00 |
G. orbiforme (6) | G. orbiforme (6) | Brazil | 89%/0.96 | |
G. subfornicatum (9) | G. ecuadoriense (5)*, Ganoderma sp. (2), G. subfornicatum (1)*, uncultured fungus (1) | Brazil, Ecuador, French Guiana, India, Peru | AD | 100%/0.96 |
G. mastoporum (123) | G. australe (60), G. orbiforme (19), G. mastoporum (13), Ganoderma sp. (11), G. cupreum (10), uncultured soil fungus (6), G. fornicatum (3), G. multicornum (1) | Australia, China, India, Indonesia, Laos, Malaysia, Myanmar, Taiwan, Thailand, Vietnam | AD, GS | 80%/1.00 |
CLADE D, Cluster D.4 | 93%/1.00 | |||
Group D.4.1 | 87%/1.00 | |||
G. angustisporum (16) | G. australe (8), Ganoderma sp. (5), G. angustisporum (3)* | Australia, China, India, Indonesia, Malaysia | AD, GS | |
Ganoderma sp. D1 (2) | G. applanatum (2) | Gabon | AD | 100%/1.00 |
Group D.4.2 | 91%/1.00 | |||
G. zonatum (84) | G. zonatum (84) | USA | AD, AM | 100%/1.00 |
Ganoderma sp. D2 (4) | Ganoderma sp. (4) | Colombia | AM | 100%/1.00 |
G. ryvardenii (22) | G. ryvardenii (15)*, Ganoderma sp. (6), G. wiiroense (1) | Cameroon, India | AD, AM | 100%/1.00 |
G. boninense (69) | G. boninense (32), Ganoderma sp. (29), G. miniatocinctum (3), G. zonatum (3), G. orbiforme (2) | China, Indonesia, Japan, Malaysia, Taiwan, Thailand, Vietnam | AM | |
Ganoderma sp. D3 (12) | Ganoderma sp. (12) | Indonesia | AM | 98%/1.00 |
CLADE E | 81%/1.00 | |||
CLADE E, Cluster E.1 | 99%/1.00 | |||
G. williamsianum (42) | G. australe (29), G. williamsianum (7), G. cf. australe (2), G. australe cplx (2), Ganoderma sp. (1), uncultured fungus (1) | China, Malaysia, Myanmar, Thailand | AD | 99%/1.00 |
CLADE E, Cluster E.2 | 84%/0.99 | |||
Ganoderma sp. E1 (23) | G. applanatum cplx (8), G. tornatum (7), Ganoderma sp. (4), G. lobatum (3), G. gibbosum (1) | Brazil, Colombia, Costa Rica, Ecuador, French Guyana, Mexico, Peru, USA | AD | 87%/0.99 |
Ganoderma sp. E2 (37) | G. gibbosum (12), G. tornatum (8), G. lobatum (7), Ganoderma sp. (6), G. applanatum cplx (2), G. australe (2) | Argentina, Brazil, Colombia, Cuba, Puerto Rico, USA | AD, AM | 96%/1.00 |
G. aff. gibbosum (51) | Ganoderma sp. (46), G. australe (3), G. gibbosum (1), G. ryvardenii (1) | India | AD, AM, GS | 69%/1.00 |
G. eickeri (4) | G. eickeri (2)*, Ganoderma sp. (2) | South Africa | AD | 92%/1.00 |
G. gibbosum (113) | G. gibbosum (61), G. applanatum (28), G. australe (10), Ganoderma sp. (3), G. australe cplx (2), G. australe IG1 (2), G. lingzhi (2), Agaricales sp. (1), G. fulvellum (1), G. lucidum (1), Fuscoporia viticola (1), uncultured Ganoderma (1) | China, Japan, Laos, Pakistan, South Korea, Taiwan, Thailand | AD, AM | |
G. ellipsoideum (80) | G. gibbosum (22), G. australe (10), G. australe cplx (10), Ganoderma sp. (10), G. adspersum (5), G. applanatum (5), uncultured soil fungus (5), G. ellipsoideum (5)*, G. applanatum cplx (3), G. tornatum (3), G. aff. steyaertanum (1), Tomophagus sp. (1) | Australia, Cambodia, China, India, Indonesia, Laos, Malaysia, Myanmar, Papua New Guinea, Sri Lanka, Thailand, USA, Vietnam | AD, AM, GS | -/0.99 |
Ganoderma sp. E3 (7) | G. australe (6), uncultured soil fungus (1) | Australia, Indonesia | 96%/1.00 | |
Ganoderma sp. E4 (13) | G. australe (12), G. tornatum (1) | Indonesia, Malaysia | AD, AM | 100%/1.00 |
CLADE E, Cluster E.3 | 71%/0.97 | |||
G. knysnamense (4) | G. knysnamense (4)* | South Africa | AD | 100%/1.00 |
G. mutabile (2) | G. mutabile (2)* | China | 100%/1.00 | |
G. cupreolaccatum (1) | G. cupreolaccatum (1) | |||
G. pfeifferi (17) | G. pfeifferi (17) | Czech Republic, Greece, Slovakia, UK | AD | 93%/- |
CLADE E, Cluster E.4 | -/0.99 | |||
G. chocoense (1) | G. chocoense (1)* | Ecuador | ||
G. podocarpense (2) | G. podocarpense (1)*, uncultured soil fungus (1) | Ecuador, Panama | 100%/1.00 | |
Ganoderma sp. E5 (8) | Ganoderma sp. (4), G. lobatum (2), G. tornatum (1), uncultured soil fungus (1) | Argentina | AD | 100%/1.00 |
Ganoderma sp. E6 (35) | G. australe (19), Ganoderma sp. (8), G. australe IG2 (2), G. applanatum (2), G. australe cplx (1), G. cf. australe (1), G. cf. philippii (1), uncultured Ganoderma (1), uncultured soil fungus (1) | China, India, Laos, New Zealand, Papua New Guinea, Taiwan, Thailand, Vietnam | AD | 95%/1.00 |
G. australe (80) | G. australe (27), G. australe cplx (14), Ganoderma sp. (13), uncultured soil fungus (10), G. annulare (2), G. applanatum cplx (2), G. brownii (2), G. lobatum (2), G. tornatum (2), fungal sp. (1), G. adspersum (1), G. applanatum (1), G. lipsiense (1), G. lucidum (1), uncultured Ganoderma (1) | Argentina, Australia, Brazil, Chile, Costa Rica, India, New Zealand, South Africa, UK, USA | AD, AM | |
CLADE E, Cluster E.5 | 97%/1.00 | |||
Ganoderma sp. E7 (17) | Ganoderma sp. (10), G. applanatum cplx (3), G. lobatum (3), uncultured soil fungus (1) | USA | 100%/1.00 | |
G. aff. adspersum (11) | G. adspersum (4), G. applanatum (4), G. australe cplx (2), uncultured soil fungus (1) | China, Japan, Korea | 72%/- | |
G. adspersum (144) | G. adspersum (113), G. australe (13), Ganoderma sp. (10), G. applanatum (4), basidiomycetes sp. (2), uncultured Ganoderma (1), uncultured fungus (1) | Armenia, Belgium, Croatia, France, Georgia, Germany, Greece, India, Iran, Italy, Slovakia, Spain, Tunisia, UK, USA | AD, GS | 89%/1.00 |
The principal phylogenetic analysis of ITS sequences, spanning the entire genus, was based on the construction of the main dataset (‘DS’), which was prepared by preferably using ASVs since they were considered more reliable. Selected singletons were also included in the DS when they corresponded to type material, when representing material of diverse origin or under various taxonomic names and in the absence of adequate number of ASVs (four was set as minimum and ten as maximum) for a particular clade. On the basis of the outcome of the phylogenetic analysis performed on the DS, six additional partial datasets (designated as ‘pDS’) were constructed in order to examine in more detail relationships within each of the main phylogenetic groups of the genus as they derived from the analyses performed by including all additional entries available (Table
Hence, multiple alignments of seven different matrices (Table
Phylogenetic analyses were based on Maximum Likelihood (ML) and Bayesian Inference (BI). The ML analyses were conducted with RaxML HPC BlackBox running on XSEDE (Stamatakis et al. 2014) under the general time-reversible (GTR) model of nucleotide substitution with gamma distributed rate heterogeneity (GTRGAMMA) for branch confidence with non-parametric bootstrap support (BS) according to MRE-based bootstrapping criteria assessed through the CIPRES Science Gateway-web portal/platform (
Datasets of Ganoderma sequences constructed and details of the phylogenetic analyses conducted in the frame of this study.
Datasets constructed and analysed | No. of sequences used/total | Represented entries/total entries | Alignment length | Constant characters* | Parsimony Informative characters* | No. of Rapid Bootstraps | ML Optimisation Likelihood | Model substitution AICc TS1/5.8S/ITS2 | No. of generations | Split frequency | 50% credible trees | ML trees presented |
---|---|---|---|---|---|---|---|---|---|---|---|---|
DS: entire genus and individual Clades/Clusters | 440/2119 | 2027/3908 | 713 | 328 | 307 | 504 | -10429.843785 | GTR+G/GTR/TVM+G | 37885000 | 0.004398 | 28415 | Figs |
pDS1a: Cluster A.1 (expanded) | 120/124 | 297/301 | 571 | 500 | 25 | 600 | -1701.266373 | JC/JC/JC | 13710000 | 0.004400 | 10285 | Suppl. material |
pDS1b: Cluster A.2 (expanded) | 263/274 | 517/528 | 608 | 436 | 86 | 600 | -2929.177334 | TIM2ef+G/GTR+G/TIM3+G | 63120000 | 0.004399 | 47341 | Suppl. material |
pDS1c: Cluster A.3 (expanded) | 341/641 | 694/1385 | 648 | 376 | 156 | 552 | -4501.116292 | GTR+G/TVM/TVM+G | 49780000 | 0.004400 | 37337 | Suppl. material |
pDS2/pDS3: Clades B & C (expanded) | 26/292 67/74 | 224/431 88/95 | 607 | 560/459 | 38/89 | 600 | -2652.570067 | K80/JC/JC | 7390000 | 0.004365 | 5545 | Suppl. material |
pDS4: Clade D (expanded) | 292/316 | 449/474 | 631 | 391 | 149 | 504 | -4424.884308 | TPM3uf+I+G/TPM2/TPM1uf+I+G | 40320000 | 0.004398 | 30241 | Suppl. material |
pDS5: Clade E (expanded) | 367/396 | 664/693 | 656 | 386 | 162 | 552 | -5123.351039 | TIM3+G/GTR/TPM3uf+G | 65585000 | 0.004398 | 49190 | Suppl. material |
Elaboration of ITS metadata made it possible to assess the phylogeny of Ganoderma species under study and the support it received; values of bootstrap support (BS) in ML and Bayesian Posterior Probability (BPP) in BI analyses were considered significant and retained when equal or higher than 65% and 0.95 in ML and BI analyses, respectively. Moreover, inter- and intra-specific pairwise genetic distances (on the basis of uncorrected p-values) within and between allied Ganoderma taxa were calculated in MEGA X as the proportion (p) of nucleotide sites at which two sequences was different and was obtained by dividing the number of nucleotide differences by the total number of nucleotides. In addition, ITS sequence similarities were computed in MAFFT through the EMBL-EBI portal. Indicative cases are illustrated by boxplot graphs depicting genetic distances and sequence similarities in the main clades, as well as within and amongst selected species of the genus Ganoderma.
Widely-adopted thresholds for separating amongst species in Basidiomycota are < 97% to 98% for ITS sequence similarity and > 0.010 to 0.020 for genetic distance uncorrected p-values (
Moreover, in order to provide additional information about the variation existing in the ITS spacers for the entire genus, as well as for each major clade/cluster derived from the phylogenetic analyses, the length and GC content of ITS1 and ITS2 were calculated in Geneious Prime version 11.1.4 (https://www.geneious.com) by examining all sequences used for the construction of the main phylogenetic tree (Fig.
In total, 3970 ITS entries were retrieved from the GenBank/ENA/DDBJ and UNITE databases; 62 sequences were removed from further analysis since they were either erroneously annotated as Ganoderma (58) or they could not be reliably identified (Suppl. material
a Initial labelling of 3908 Ganoderma sequences analysed in the present study: numbers in parentheses correspond to sequences deposited under the particular name in GenBank/ENA/DDBJ and UNITE, while species names appear underlined when ITS sequences derive from type material b final assigment of 3908 Ganoderma sequences to 80 species and six distinct groups as a result of the phylogenetic analyses performed in this study: numbers in parentheses correspond to the number of sequences grouped within each taxon (data deriving from Table
Almost half (45.3%) of all Ganoderma sequences deposited in GenBank/ENA/DDBJ and UNITE correspond to only eight species names, i.e. G. lucidum (12.3%), G. lingzhi (8.6%), G. australe (5.4%), G. resinaceum (4.5%), G. applanatum (4.3%), G. sessile (3.5%), G. adspersum (3.3%) and G. curtisii (3.3%) (Fig.
Basidiomes of Ganoderma spp. amongst those collected and analysed in this study (specimens codes appear in parantheses; Suppl. material
For inferring the phylogeny of the entire genus, the main dataset (DS) included 2027 entries (i.e. ca. 52% of the total number of Ganoderma entries analysed) corresponding to 161 singletons and 279 ASVs (Table
A comprehensive evaluation of all ITS sequences available permitted us to determine variation in the ITS1 and ITS2 spacers by examining their length and GC content; relevant data are presented below. In addition, the comparative assessment of ITS1 and ITS2 heterogeneity amongst Ganoderma species could subsequently contribute to species determination since they contain information of potential diagnostic value. The multiple sequence alignment revealed two polymorphic segments which could be potentially used for the identification of Ganoderma specimens at species level (Table
Summary of polymorphic regions in ITS1 and ITS2 spacers assessed in Ganoderma species/groups represented by ≥ 10 entries in the GenBank/ENA/DDBJ and UNITE databases. For each region, the length and position between conserved oligonucleotides (i.e. TGCAC to GAATG in ITS1 and AATCT to TAGCT in ITS2) are hereby provided. Additional pertinent information is included in Suppl. material
Species/Groups | ITS1 sequence of potential diagnostic value | Length (nt) | Position in the alignment of Suppl. material |
ITS2 sequence of potential diagnostic value | Length (nt) | Position in the alignment of Suppl. material |
---|---|---|---|---|---|---|
CLADE A | ||||||
Cluster A.1 | ||||||
G. oregonense | 13 | 25–37 | GCCTTTGCGGGTW | 26 | 17–42 | |
G. tsugae | TGTGAAGCGTGCT | 13 | 25/26–37/38 | TGYRGGCTTGGAC | 26 | 17–42 |
G. carnosum | 13 | 25–37 | AGCCTTGC | 8 | 16–23 | |
G. lucidum | TGAAGCGYNCCYY | 13 | 27–39 | nd | ||
G. leucocontextum – G. weixiensis | CGAAGCGTGC | 10 | 27–36 | nd | ||
Cluster A.2 | ||||||
G. hoehnelianum | CTTCAGTC | 8 | 16–23 | CTTGTGGGTT | 10 | 20–29 |
G. weberianum | nd | nd | ||||
G. sichuanense | nd | nd | ||||
G. carocalcareum | AACGTCGTKAAGCGGGC | 17 | 21–37 | nd | ||
Ganoderma sp. A1 | GGGTCTTTT | 9 | 34–42 | CGTCTTTC | 8 | 60–67 |
G. mexicanum | GCTCTTTACTGAGCC | 15 | 36–50 | CGGCCGGCTCCTCT | 21 | 65/67–85/87 |
G. parvulum | 15 | 36–50 | TAAATGC1 | 21 | 65/67–85/87 | |
G. resinaceum | AAGCGGCG | 8 | 55/56–62/63 | nd | ||
Ganoderma sp. A3 | GGATCGGCGT | 10 | 55–64 | ACAGATCT | 8 | 13–20 |
G. polychromum | ACACCTAT | 8 | 84–91 | nd | ||
G. sessile | CCACAAACTCTR | 12 | 78–89 | CTTACAAA | 8 | 10–17 |
Cluster A.3 | ||||||
G. tuberculosum | GATTGTCG | 8 | 21–28 | CCATGCCC | 8 | 58/59–65/66 |
G. wiiroense | GGCATTAT | 8 | 21–28 | TTCTCTTA | 8 | 71/72–78/79 |
G. philippii | TTGCTGGG | 8 | 39–46 | CTTTTGTGGYTTT | 13 | 18–30 |
G. lingzhi | CAGATTGC | 8 | 19–26 | 10 | 54–63 | |
G. curtisii | TGCGGAGCGCAT | 12 | 49–60 | CGGCCGTTAT | 10 | 54–63 |
G. ravenelii | GAGTGCAT | 8 | 53–60 | 10 | 54–63 | |
G. multiplicatum | CCCTTTAT | 8 | 35–42 | nd | ||
G. destructans – G. dunense | 9 | 22–30 | nd | |||
G. steyaertanum | ATCVTAAAA2 | 9 | 22–30 | CTCTTGGCC | 9 | 61–69 |
G. martinicense | 9 | 22–30 | CATTCTTG | 8 | 59–66 | |
G. multipileum | 9 | 22–30 | G(C)AAGCTTTTG | 10–11 | 13–22/23 | |
Ganoderma sp. A6 | TCCCAGGA | 8 | 50–56 | CTCCTCTCTT | 10 | 72–81 |
G. tropicum | ACCGGGCTTTGCA | 13 | 42–54 | nd | ||
CLADE B | ||||||
G. applanatum | GTGCTYTT | 8 | 32–39 | TAAGCTTKTGT | 11 | 14–24 |
CLADE C | ||||||
G. neojaponicum | ATGGATCGCG | 10 | 18–27 | AGGTGTTTG | 9 | 47–55 |
G. enigmaticum – G. thailandicum | CTTCTTGTC | 9 | 35–43 | TTGCAACC | 8 | 11–18 |
G. casuarinicola | GCTCTTGT | 8 | 34–41 | 8 | 11–18 | |
CLADE D | ||||||
G. mbrekobenum | TTWCAGASSGT | 11 | 16–26 | AGGCTATT | 8 | 48–55 |
G. nasalanense | CGTTTTCA | 8 | 70–77 | TCTTTAATA | 9 | 60/62–68/70 |
G. sinense | GGAGCTSGT | 9 | 41–49 | GTAAAGGC | 8 | 24–31 |
G. mastoporum | nd | TTTTTARYGRKTTTGTAGG | 19 | 19–37 | ||
G. angustisporum | GTGTAAAA | 8 | 27–34 | ATGGCTWGT | 8 | 24/28/29–32/36/37 |
G. zonatum | TCGCTCGC | 8 | 34–41 | TCTCTTCA | 8 | 3–10 |
G. ryvardenii | TCGTGCGG | 8 | 23–30 | CTTTAACT | 8 | 61–68 |
G. boninense | GTTTGACRAGTT | 12 | 40/44–51/55 | ATCTCTTTGY | 10 | 16–25 |
Ganoderma sp. D3 | GGCGTGGT | 8 | 24–31 | 10 | 16–25 | |
CLADE E | ||||||
G. williamsianum | CTTCAGGTC | 9 | 16–24 | CTTAATYGA | 9 | 21–29 |
Ganoderma sp. E1 | GTTTTACG | 8 | 15–22 | ATRAGCTTCT | 10 | 13–22 |
Ganoderma sp. E2 | 8 | 15–22 | TATGKGAG | 8 | 23–30 | |
G. aff. gibbosum | 13 | 27–39 | 10 | 60–69 | ||
G. gibbosum | TGARRSGGGCTYG3 | 13 | 27–39 | TCCYTTTACR3 | 10 | 60–69 |
G. ellipsoideum | 13 | 27–39 | 10 | 60–69 | ||
Ganoderma sp. E4 | RTTAAACG | 8 | 26–33 | GTCGGACTW4 | 9 | 59–67 |
G. pfeifferi | GGCCCGTTT5 | 9 | 34/35–42/43 | GCCTTTGTC6 | 9 | 57–65 |
Ganoderma sp. E6 | ACYGAGCYYGC | 11 | 41–51 | TCTTTGCGGGG | 11 | 19–29 |
G. australe | CGAAACGKGCTCG | 13 | 27–39 | 11 | 19–29 | |
Ganoderma sp. E7 | CCCCATGA | 8 | 83/84–90/91 | GTCTTTACA | 9 | 59–67 |
G. aff. adspersum | GGGCCCGTTC | 10 | 33–42 | CTTCTTGCGG | 10 | 18–27 |
G. adspersum | AGGCCCGTTC | 10 | 33–42 | AGGTTTGTAGGG | 12 | 27/28–38/39 |
The ITS analyses resulted in the formation of well-resolved/supported terminal subclades which led us to accept 80 Ganoderma taxa at species level in accordance with the criteria set and in conjunction with literature data available (Fig.
Summary tree of the genus Ganoderma inferred from ML analysis, based on ITS sequence data (main dataset, DS; Table
In all cases, ML and BI analyses provided almost identical tree topologies with minor differences and, thus, only the trees inferred from the ML analysis are presented. The genus Ganoderma exhibits a strongly-supported monophyly (BS: 90%, BPP: 1.00; Fig.
On the basis of ITS meta-analysis, Clade A is moderately supported only through ML analysis (69%; Fig.
Detail from Fig.
Clade A includes 1927 entries distributed across 881 unique ITS sequence types of which 240 appear as ASVs representing 1414 entries in GenBank/ENA/DDBJ and UNITE. Clade A could be further divided into three well-supported Clusters (A.1, A.2 and A.3) and to the recently-introduced G. shanxiense L. Fan & H. Liu (
Detail from Fig.
In the context of this work, Cluster A.1 corresponds to a well-supported clade (100%, 1.00; Fig.
The major group of Cluster A.1 corresponds to G. lucidum (Curtis) P. Karst. (G. lucidum sensu stricto), which is represented by 10 sequences (or 97 entries) in the main dataset (Fig.
G. carnosum Pat. was first described on Abies in southwest France (Pyrenees Mts.) and is distributed throughout Europe (
Sequences deriving from two U.K. specimens (deposited as G. carnosum and G. lucidum) and from two environmental samples (Estonia, labelled as “uncultured soil fungus”) formed a distinct group (93%, 1.00; Fig.
G. oregonense Murrill and G. tsugae Murrill are closely-associated taxa (
In the past, G. tsugae was occasionally reported to be conspecific with G. carnosum or G. lucidum on the basis of morphological observations (
Noteworthy cases pertain to sequences under the names of “G. valesiacum” and “G. ahmadii” which grouped within G. tsugae (Table
In the frame of this study, four ASVs representing 20 entries formed a well-supported group (81%, 1.00; Fig.
Cluster A.2 (94%, 1.00; Fig.
A major subclade is formed by specimens collected in southeast Asia and Australia growing mostly on angiosperms (93%, 1.00 and 72%, 1.00, Fig.
G. hoehnelianum Bres. forms a monophyletic group in Cluster A.2 (99%, 1.00; Table
A new phylogenetic species within Cluster A.2 is hereby proposed and is provisionally named “Ganoderma sp. A1” (corresponding to the UNITE DOIs SH1740420.08FU, SH1740444.08FU and SH1740445.08FU); its monophyly is strongly supported in both trees (100%, 1.00; Fig.
Other four entries of dubious identity derive from material originating from Asia (China and India) labelled as Ganoderma sp., “G. weberianum” (2) and “Ganoderma cf. weberianum”, as well as from Brazil under the name “G. subamboinense” (Table
Another major subclade consisting of material originating from the Neotropics is formed by a total of 17 entries. Four of them are singletons, while the other 13 are grouped in six ASVs (Suppl. material
Another well-supported terminal clade (100%, 1.00; Fig.
G. resinaceum Boud. is represented by 10 ASVs and 65 entries (Fig.
Another new phylogenetic species is hereby proposed, provisionally named “Ganoderma sp. A3” (corresponding to the UNITE DOI SH1723084.08FU). It is strongly supported in both trees (100%, 1.00 and 99%, 1.00 in Fig.
G. sessile Murrill is a well-supported (78%, 1.00; Fig.
A closely-related and well-supported (93%, 1.00; Fig.
A distinct phylogenetic group (98%, 1.00; Fig.
Another terminal group, although not adequately supported in any of the trees constructed, forms a sister clade to G. sessile/G. polychromum (Fig.
Cluster A.3 (76%, 1.00; Fig.
G. tuberculosum Murrill is strongly supported in the generated trees (100%, 1.00; Fig.
Furthermore, one sequence (JX310812) labelled as “G. chalceum” originating from Brazilian material and one sequence under the name G. concinnum Ryvarden (possibly of South American origin) form a terminal subclade which nested close to G. tuberculosum (0.97; Fig.
Two sister subclades (0.99, Suppl. material
Two new well-supported monophyletic species, provisionally named as “Ganoderma sp. A4” (no UNITE DOI available) and “Ganoderma sp. A5” (corresponding to the UNITE DOI SH1723120.08FU) are revealed in this study. Ganoderma sp. A4 is represented by two sequences deriving from Argentinian material which were originally identified as “G. lucidum” (100%, 1.00; Fig.
G. philippii (Bres. & Henn. ex Sacc.) Bres. constitutes a well-supported species (97%, 1.00; Fig.
Three species, i.e. G. lingzhi S.H. Wu, Y. Cao & Y.C. Dai, G. ravenelii Steyaert and G. curtisii (Berk.) Murrill, form a strongly-supported group (A.3.1; 98%, 1.00; Fig.
G. curtisii (70%, 1.00; Suppl. material
Six other closely-related species are found in Cluster A.3 (i.e. G. multiplicatum, G. destructans, G. steyaertanum, G. mizoramense, G. martinicense and G. multipileum) (Fig.
In this study, G. destructans M.P.A. Coetzee, Marinc., M.J. Wingf. represented by 39 entries (including the type material; Table
G. steyaertanum B.J. Smith & Sivasith. forms a well-supported group (89%, 0.99; Fig.
G. martinicense Welti & Courtec. is sister to G. multipileum (72%, 0.98; Fig.
G. multipileum Hou is phylogenetically supported (74%, 0.97; Fig.
A second group (A.3.2; 100%, 1.00; Fig.
A strongly-supported terminal subclade in Group A.3.2 (100%, 1.00; Fig.
Furthermore, three singletons from Malaysian material (in this particular case, geographic origin is inferred from the title of the study which appears on the respective GenBank records), initially identified as “G. fornicatum”, form another distinct well-supported group (100%, 1.00; Fig.
Clade B (96%, 1.00; Fig.
Detail from Fig.
Detail from Fig.
The other two species in Clade B form a well-supported sister clade (95%, 1.00; Fig.
Box plots of a ITS sequence similarity (%) and b genetic distances (p-values) within (intra) and between (inter) Ganoderma species for each one of the main lineages (Clades/Clusters) of the genus, as well as pairwise comparisons between selected species. The size of each box represents 50% of the values, the black horizontal line within each box indicates the median, the ‘x’ represents the average value, the error bars represent interquartile ranges and circles indicate outliers. The red-dotted horizontal line, transversing the plots, represents the value levels accepted in this study for proposing new phylogenetic species.
Clade C is strongly supported (100%, 1.00; Fig.
Cluster C.2 is strongly supported (100%, 1.00; Fig.
The other group (100%, 1.00; Fig.
Clade D includes 291 sequences representing 449 individual entries separated into four clusters (Suppl. material
Cluster D.1 is placed at the base of Clade D and it is strongly supported (98%, 1.00; Fig.
Cluster D.2 comprises two distantly-related species, i.e. G. sinense Zhao, Hsu & Zhang (
Cluster D.3 is well supported (75%, 1.00; Fig.
Another terminal subclade is composed of nine entries (98%, 1.00; Fig.
Five sequences, representing G. orbiforme, form a well-supported terminal group (89%, 0.96; Fig.
Finally, a well-represented and supported terminal subclade of Cluster D.3 (80%, 1.00; Fig.
Cluster D.4 is strongly supported (93%, 1.00; Fig.
Group D.4.1 (87%, 1.00; Fig.
Group D.4.2 (91%, 1.00; Fig.
The other two species comprise material from Asia only; the ‘core’ part corresponds to G. boninense Pat. and it is represented by 69 entries deposited as G. boninense (32), “G. miniatocinctum” (3), “G. orbiforme” (2), “G. zonatum” (3) and Ganoderma sp. (29) (Table
On the basis of the results presented above, it is apparent that G. zonatum is not related to G. sessile (the latter forms part of Cluster A.2), as previously reported by
Clade E is strongly supported (81%, 1.00; Fig.
Cluster E.1 corresponds to a single strongly-supported phylospecies (99%, 1.00; Fig.
Cluster E.2 (84%, 0.99; Fig.
A well-supported terminal subclade (87%, 0.99; Fig.
A strongly-supported (96%, 1.00; Fig.
Cluster E.2 also includes a large terminal subclade (65%, 1.00; Suppl. material
The other terminal clade (92%, 1.00; Fig.
G. ellipsoideum is a recently-described species from Hainan Island, China (
The other two phylogenetic species appearing on the terminal subclade of Cluster E.2 (96%, 1.00; Fig.
Cluster E.3 consists of material corresponding to the laccate taxa G. pfeifferi Bres. (17 sequences from Europe only; Table
The close phylogenetic position of G. pfeifferi and G. adspersum (Schulzer) Donk (Cluster E.5) is congruent with their similar pileus dark-brown context and the complex structure of the crust in contrast to the normal palisade appearance in laccate species of Clade A (
Cluster E.4 (0.99; Fig.
A sister group to the former (1.00; Fig.
A sister terminal subclade to the previous group (G. chocoense, G. podocarpense and Ganoderma sp. E5) consists of sequences originating from material collected in south and southeast Asia, while a single sequence indicates its presence also in Papua New Guinea (95%, 1.00; Fig.
Finally, a distinct group included 80 entries deriving from specimens mainly from southeast Asia, South America and Oceania, as well as from central/north America and South Africa. That said, one sequence originated in the vicinity of the Kew Botanical Gardens, UK, such that we suspect it to have been imported with plant material. These sequences were mostly deposited as “G. australe” (27), “G. australe complex” (14) and Ganoderma sp. (13). Despite the widespread use of the “australe” epithet to describe several entities (14) placed in other terminal clades of the present study, we believe that this particular terminal group coincides with G. australe (Fr.) Pat. after evaluating all available ITS sequence data, the geographic distribution of specimens analysed and pertinent publications (
Cluster E.5 (97%, 1.00; Fig.
A sister subclade to the aforementioned clade corresponds to G. adspersum sensu lato and is represented by 155 entries deriving from Europe, south and west Asia and North Africa (1.00; Suppl. material
In conclusion, the plasticity of morphological characters and substantial overlap of alleged diagnostic features is prevalent in the ‘dull’ taxa of this group; consequently, their taxonomic significance is dubious. Moreover, the massive, heavily-agglutinated matrix of the pileal crust in such non-laccate species often obstructs observation of discriminating features in pileal elements. The situation is further aggravated by the loss of type material of widely and commonly used species names (e.g. G. australe and G. gibbosum), the absence of (correct) neo-typification (as there are ambiguous synonymies regarding several taxa) and unclear species descriptions (e.g. G. tornatum and G. lobatum). All of the above could explain the obstacles which prevented the establishment of a stable classification system for species of Clade E. Still, the highly variable morphology of basidiomes and the ITS divergence are indicative of underestimated diversity and the presence of cryptic species is quite certain as is also indicated by the outcome of the present study.
As determined from the analysis performed in this study, the combined length of the two spacers (excluding the intercalary 5.8S gene) ranged from 378 (Ganoderma sp. D1) to 429 bases (Ganoderma sp. E4) with an average value of 400 (± 6.8) bases (Fig.
Box plots of a length (bases) and b GC (%) content of ITS1 and ITS2 sequences for each one of the main lineages (Clades/Clusters) of the genus Ganoderma. The size of each box represents 50% of the values, the black horizontal line within each box indicates the median, the ‘x’ represents the average value, the error bars represent interquartile ranges and circles indicate outliers.
In the ITS1 region (289 sites), 195 (67%) were variable and 159 (55%) were parsimony informative; the ITS2 region (262 sites) included 169 (65%) variable and 137 (52%) parsimony informative sites. This is in accordance with the outcome of previous reports indicating a larger variability for ITS1 in comparison to ITS2 in Fungi (
The GC content was almost identical in ITS1 and ITS2 spacers when calculated for the entire genus (49.1 ± 1.5% and 49.2 ± 1.7%, respectively) (Fig.
The application of criteria of wide applicability/suitability for delineating taxa in the genus Ganoderma is a particularly challenging task because species exhibit complex evolutionary backgrounds, widespread occurrence and/or problematic taxonomy as previously mentioned. Since this study was principally based on the analysis of ITS metadata, difficulties related to erroneous, fragmentary and/or incomplete information on the origin and the true identity of the sequenced material had to be addressed, together with labelling referring to the genus only (i.e. “Ganoderma sp.”) or with unidentified sequences (e.g. “Agaricales sp.”, “uncultured fungus” and “unidentified soil fungus”). The outcome demonstrated that the use of ITS rDNA could confer valuable data on the establishment of phylogenetic species within the genus since the majority of terminal clades were strongly supported and species boundaries were – in the majority of cases – elucidated, although relationships/affinities amongst particular sections or within certain species complexes were not adequately resolved.
Intraspecific ITS sequence similarity values were relatively high (i.e. overall average: 99.32 ± 0.59%) for the Ganoderma species examined. Moreover, genetic distance, based on uncorrected p-values, provided an additional effective tool for species delineation which was found to be congruent with ITS sequence similarities (overall average of intraspecies genetic distance within the genus: 0.005 ± 0.004). Hence, in the context of the criteria adopted, these parameters contributed significantly to delineating Ganoderma species and were generally in accordance with the concepts of already accepted morpho- and/or phylospecies. In addition, a species hypothesis using a threshold of interspecific values for ITS similarity (≤ 98%) and genetic distance (≥ 0.015) was effectively applied for the 21 putatively new phylospecies hereby proposed, while several other distinct entities of dubious status are revealed in the trees inferred. Finally, 59 terminal groups correspond to already established species of the genus which demonstrate a large variability in genetic distance and sequence similarity when selected pairs of phylogenetically-related taxa are compared within different clades (Fig.
Especially as regards taxa of Cluster A.1, average interspecific distances (0.008 ± 0.004) are close to the respective intraspecific values (0.004 ± 0.003); similarly, sequence similarity in interspecific comparisons is high (99.0 ± 0.6; Fig.
On the other hand, species in Clade B present a clear barcoding gap and distinct ‘sequence diameters’ (as defined by
As regards Clades C, D and E, barcoding gaps are quite pronounced and delimitation of many species could be made on the basis of the 98% sequence similarity and the 0.015 genetic distance values previously mentioned and used for the establishment of the phylospecies proposed herein (Fig.
In conclusion, ITS phylogeny, in conjunction with sequence similarity and genetic distance measurements, do not fully support the delineation of some well-established Ganoderma taxa; instead, they advocate their inclusion in monophyletic units representing species complexes. On the other hand, G. neojaponicum (intraspecific values: 0.019 ± 0.009; 97.51 ± 1.44%), Ganoderma sp. A6 (0.028 ± 0.005; 97.24 ± 1.02), Ganoderma sp. C1 (0.029, 95.50%) and Ganoderma sp. D1 (0.020, 97.76%) seem to harbour cryptic variation and might correspond to more than one phylospecies. Relatively high intraspecies genetic distance values were also detected in G. flexipes, G. mastoporum, G. mbrekobenum and Ganoderma sp. D3 (Fig.
Members of the genus Ganoderma exhibit relatively-complex microanatomy (
Three main lineages of the genus were identified: Clade A, Clade B and Clades C through E. Clades A and E include taxa with a cosmopolitan distribution, while species of Clade B are distributed across the Holarctic region; species of Clades C occur in the Paleotropics and taxa of Clade D exhibit a pantropical distribution. Their subsequent analysis results are consistent with the hypothesis of a Northern Hemisphere origin (tropical Asia) for Ganoderma species with subsequent range expansions to the Southern Hemisphere and by colonisation of the Neotropics through long distance dispersal (
Indicative cases of Ganoderma species distribution include:
Palearctic – Eurasian, Old World: G. lucidum, G. resinaceum and G. adspersum occur across Eurasia and share several common host plants, for example, the genera Quercus, Salix, Populus, Abies and Larix. A strictly European distribution is exhibited by G. carnosum, G. aff. carnosum and G. pfeifferi. Allopatric speciation seems to be under way between Eurasian G. resinaceum and Taiwanese collections corresponding to Ganoderma sp. A3. An Old Word distribution is also presented by species of the Cluster E.3, since each one of them has been reported exclusively from Africa (G. knysnamense), from Asia (G. mutabile) or from Europe (G. pfeifferi).
East Asia – Malay Archipelago – Oceania: G. weberianum complex (G. sichuanense – G. weberianum) and G. steyaertanum (Clade A), G. angustisporum and G. mastoporum (Clade D) and Ganoderma sp. E3 and Ganoderma sp. E6 (Clade E). Particularly as regards G. ellipsoideum, its distribution extends to USA on the basis of sequences deriving from environmental samples.
East Asia and South Africa (Paleotropic): All taxa of Clade C, as well as G. hoehnelianum (Cluster A.2), G. wiiroense (Cluster A.3), G. angustisporum and Ganoderma sp. D1 (Clade D). Moreover, G. carocalcareum and G. austroafricanum (Cluster A.2), G. destructans – G. dunense (Cluster A.3), G. aridicola, G. enigmaticum – G. thailandicum, Ganoderma sp. C1 and Ganoderma sp. C2 (Clade C), G. cupreum, Ganoderma sp. D1 and G. ryvardenii (Clade D), as well as G. eickeri (Clade E) have so far only been recorded in Africa.
Holoarctic/Nearctic – Palearctic: Within Clade B, G. applanatum presents an inter-continental distribution indicating gene exchange between the Palearctic and the Nearctic regions through land bridges which were widely accepted as corridors for such transfers (
East Asian – North American: Several biogeographic studies evidenced migration of fungi from east Asia to North America via the Bering Land Bridge route (Wu and Mueller 1997;
Neotropical: G. mexicanum – G. parvulum complex, G. aff. polychromum (Cluster A.2), G. tuberculosum and G. martinicense (Cluster A.3), Ganoderma sp. D2 (Clade D) and Ganoderma sp. E1 and sp. E2 (Clade E). As regards other species occurring in the Americas only, G. curtisii (Cluster A.3), G. zonatum (Clade D) and Ganoderma sp. E7 (Clade E) are confined to North America. The closely-related G. podocarpense and G. chocoense (Clade E) were reported only in Central America, whereas Ganoderma sp. A4, G. concinnum and G. multiplicatum (Cluster A.3), G. orbiforme (Clade D) and Ganoderma sp. E5 (Clade E) were recorded in South America.
Southern Hemisphere: The phylogenetic analysis of Southern Hemisphere species and complexes (species in Clusters A.2 and A.3, Clusters D.3 and D.4 and E.2 and E.4) indicated a restricted gene flow apparently due to geographic isolation, although episodic long-distance dispersal still occurs (
The majority of Ganoderma species were collected on angiosperm hosts, 41 species on eudicots and 19 species on monocots, while another 18 species were reported on gymnosperms (Table
As previously stated, one of the major obstacles for exploiting sequences present in public depositories is that many of them are inaccurate; errors in labelling of metadata were estimated to correspond to as much as 20% – or even 30% according to a recent report – of GenBank accessions including also recent deposits (
The study of a large dataset comprising almost four thousand ITS entries proved to be valuable in obtaining a significant amount of phylogenetic information which contributed to elucidating the status of Ganoderma species. In addition, it contributed to establishing robust relationships amongst the majority of them, while it also revealed limitations in the use of ITS (alone) to assess certain taxa which have to be addressed through a multigene approach. However, it is interesting that the outcome of recent publications employing more than one marker (by focusing on particular groups in the genus, for example,
At a more general level, this study evidences that significant – yet largely untapped – mycological explanatory power resides in the public DNA sequence corpus and we hope that other mycologists will start scrutinising the sequence data available for their fungal groups of expertise. Our results also demonstrated that the so-called environmental sequences – usually ignored in a taxonomic/phylogenetic context – should be included in such pursuits (cf.
We are indebted to the curators of the fungaria of the Bulgarian Academy of Sciences (
This research was performed in the frame of a project (THALIS – UOA – MIS 377062) cofinanced by the European Union (European Social Fund – ESF) and Greek National Funds through the Operational Programme ‘Education and Lifelong Learning’ of the National Strategic Reference Framework (NSRF).
Tables S1–S6
Data type: species data
Explanation note: Table S1. Information about the Ganoderma material/specimens analyzed for the first time in the frame of this work: species name (the initial name appears in parenthesis when different from the one determined in the present study), specimens code, plant host, geographic origin, collector and date of collection, type of material examined (H: herbarium specimen, C: pure culture) and GenBank accession number of the generated ITS sequences. Table S2. Detailed information on the Ganoderma ITS sequences used in this study. Table S3. Ganoderma sequences excluded from the analysis since they were either erroneously labelled as Ganoderma or whose identity could not be reliably resolved. Table S4. Ganoderma sequences (642) included in the data analysis but excluded from the trees inferred (Figs
Figure S1
Data type: molecular data
Explanation note: Hypervariable regions of potential diagnostic value in ITS1 and ITS2 spacers for the Ganoderma taxa appearing in Table
Figure S2a
Data type: molecular data
Explanation note: Phylogenetic reconstruction of the genus Ganoderma inferred from ML analysis based on ITS sequence data (pDS1a; Table
Figure S2b
Data type: molecular data
Explanation note: Phylogenetic reconstruction of the genus Ganoderma inferred from ML analysis based on ITS sequence data (pDS1b; Table
Figure S2c
Data type: molecular data
Explanation note: Phylogenetic reconstruction of the genus Ganoderma inferred from ML analysis based on ITS sequence data (pDS1c; Table
Figure S2d
Data type: molecular data
Explanation note: Phylogenetic reconstruction of the genus Ganoderma inferred from ML analysis based on ITS sequence data (pDS2 & pDS3; Table
Figure S2e
Data type: molecular data
Explanation note: Phylogenetic reconstruction of the genus Ganoderma inferred from ML analysis based on ITS sequence data (pDS4; Table
Figure S2f
Data type: molecular data
Explanation note: Phylogenetic reconstruction of the genus Ganoderma inferred from ML analysis based on ITS sequence data (pDS5; Table