Skip to main content

Addressing widespread misidentifications of traditional medicinal mushrooms in Sanghuangporus (Basidiomycota) through ITS barcoding and designation of reference sequences

Abstract

“Sanghuang” refers to a group of important traditionally-used medicinal mushrooms belonging to the genus Sanghuangporus. In practice, species of Sanghuangporus referred to in medicinal studies and industry are now differentiated mainly by a BLAST search of GenBank with the ITS barcoding region as a query. However, inappropriately labeled ITS sequences of “Sanghuang” in GenBank restrict accurate species identification and, to some extent, the utilization of these species as medicinal resources. We examined all available 271 ITS sequences related to “Sanghuang” in GenBank including 31 newly submitted sequences from this study. Of these sequences, more than half were mislabeled so we have now corrected the corresponding species names. The mislabeled sequences mainly came from strains utilized by non-taxonomists. Based on the analyses of ITS sequences submitted by taxonomists as well as morphological characters, we separate the newly described Sanghuangporus subbaumii from S. baumii and treat S. toxicodendri as a later synonym of S. quercicola. Fourteen species of Sanghuangporus are accepted, with intraspecific distances up to 1.30% (except in S. vaninii, S. weirianus and S. zonatus) and interspecific distances above 1.30% (except between S. alpinus and S. lonicerinus, and S. baumii and S. subbaumii). To stabilize the concept of these 14 species of Sanghuangporus, their taxonomic information and reliable ITS reference sequences are provided. Moreover, ten potential diagnostic sequences are provided for Hyperbranched Rolling Circle Amplification to rapidly confirm three common commercial species, viz. S. baumii, S. sanghuang, and S. vaninii. Our results provide a practical method for ITS barcoding-based species identification of Sanghuangporus and will promote medicinal studies and commercial development from taxonomically correct material.

INTRODUCTION

Many macrofungi are established in traditional medicine and possess diverse properties (Wu et al. 2019a). “Sanghuang” comprises an important group of wood-inhabiting mushrooms that have been utilized in traditional medicine in China and adjacent countries for 2000 years (Zhou et al. 2020). Modern scientific studies have revealed several medicinal attributes of “Sanghuang”, including antitumor, antioxidant, anti-inflammation, and immunomodulation activities (Zhou et al. 2020). This fungal resource has also attracted the attentions of fungal chemists and pharmacologists outside Asia (Chepkirui et al. 2018; Cheng et al. 2019). Natural products, such as polysaccharides, polyphenols, pyrones and terpenes are the bioactive compounds responsible for the medicinal properties of “Sanghuang” (Zhou et al. 2020). Today, “Sanghuang” is mainly consumed in a brewed tea made from small pieces of cultivated basidiomes or occasionally powdered mycelia.

Like other wood-inhabiting traditional medicinal mushrooms, such as “Lingzhi” (Cao et al. 2012; Wang et al. 2012; Yao et al. 2013, 2020; Dai et al. 2017), “Niuchangchih” (Wu et al. 2012b, 2012c) and “Fuhling” (Redhead and Ginns 2006), there has been much debate about the taxonomic identity of “Sanghuang”. Most fungal taxonomists now agree that “Sanghuang” is represented by species of Sanghuangporus (Zhou et al. 2020). Fourteen species have been described and accepted as members of Sanghuangporus: 11 species in Asia, and one in each of Africa, Europe, and North America (Zhou et al. 2020). In addition, more new species await to be described from Africa (Chepkirui et al. 2018; Cheng et al. 2019) and perhaps other parts of the world. Besides morphological and ecological (host preference) characters, the ITS barcoding region provides the most powerful tool for differentiating species of the genus. For example, more than half of the known species of Sanghuangporus were discovered with the aid of the ITS region alone (Wu et al. 2012a, 2019b; Tian et al. 2013; Ghobad-Nejhad 2015; Tomšovský 2015; Zhu et al. 2017). Moreover, the reliability of the ITS region for species differentiation in the genus has been substantiated by a multilocus-based phylogenetic analysis (Zhu et al. 2019). Consequently, Zhou et al. (2020) reported ITS sequences from reliably identified voucher collections of the known species in the genus.

Transdisciplinary studies on Sanghuangporus have been performed to promote the utilization of this medicinal resource (Zhou et al. 2016; Cai et al. 2019; Zhu et al. 2019; Shao et al. 2020). Most of these studies aimed to identify their materials via a BLAST search of GenBank (https://www.ncbi.nlm.nih.gov/genbank/) using the ITS barcoding region as the query. However, even though each of the 14 species of Sanghuangporus has a reliable ITS sequence accession number (Zhou et al. 2020), it is not always easy to determine material in hand by a simple ITS-based BLAST search. This is a consequence of redundant and even incorrectly labeled ITS sequences in GenBank (Nilsson et al. 2006; Hofstetter et al. 2019). With inaccurately identified sequences emerging as potential matches, more collections will inevitably be inaccurately identified and the ITS sequences generated from the inaccurately identified collections will be submitted to GenBank compounding the issue and presenting new obstacles for later accurate identification. This means that there is high likelihood of medicinal and other attributes being attributed to incorrectly named species of “Sanghuang”. Meanwhile, before the erection of the genus Sanghuangporus (Zhou et al. 2016), ITS sequences generated from “Sanghuang” were labeled under other generic names, such as Inonotus and Phellinus, even though with the correct epithets. This phenomenon confuses researchers who lack taxonomic knowledge, and results in a misapplication of species names to medicinal properties, which then has a negative effect on obtaining permissions from regulatory authorities for commercial development (Zhou 2020).

As stated by Zhou (2020), the use of correct scientific names for fungal species is crucial to studies of traditional Chinese medicine and their commercial exploitation. To facilitate the rational medicinal utilization of Sanghuangporus, all ITS sequences related to “Sanghuang” in GenBank should be re-examined to assist species identification. The aim of the current study is therefore to assess the utility of the ITS region for species discrimination in Sanghuangporus, and reset the species circumscriptions on the basis of the ITS barcoding region, in order to facilitate the correction of previously mislabeled ITS sequences in GenBank, and to provide candidate diagnostic ITS sequences for use in rapid species identification of Sanghuangporus using Hyperbranched Rolling Circle Amplification (HRCA).

MATERIALS AND METHODS

Morphological examination

The newly sequenced specimens and strains are deposited in HMAS, IFP and BJFC. The specimens were observed with an Olympus BX43 light microscope (Tokyo, Japan) at magnifications up to 1000×. Microscopic procedure followed Zhou et al. (2016). Specimen sections were prepared in Cotton blue (CB), Melzer’s reagent (IKI), and 5% potassium hydroxide (KOH). All measurements were made from material mounted in heated CB. When presenting the variation of basidiospore sizes, 5% of the measurements were excluded from each end of the range and are given in parentheses. Drawings were made with the aid of a drawing tube. In the text, L = mean basidiospore length (arithmetic average of all measured basidiospores), W = mean basidiospore width (arithmetic average of all measured basidiospores), Q = variation in the L/W ratios between the studied specimens, and (a/b) = number of basidiospores (a) measured from given number (b) of specimens.

Molecular sequencing

A small piece of the basidiome or culture was taken for DNA extraction, which was performed using a CTAB rapid plant genome extraction kit-DN14 (Aidlab Biotechnologies, Beijing). The crude DNA was used as templates for the PCR amplifications of the ITS region. The primer pairs ITS1F/ITS4 and ITS5/ITS4 (White et al. 1990; Gardes and Bruns 1993) were selected for amplification and subsequent sequencing at the Beijing Genomics Institute. The PCR procedure was as follows: initial denaturation at 95 °C for 3 min, followed by 34 cycles at 94 °C for 40 s, 57.2 °C for 45 s and 72 °C for 1 min, and a final extension at 72 °C for 10 min. All newly generated sequences are deposited in GenBank (Table 1).

Table 1 Information of analyzed ITS sequences of Sanghuangporus

Downloading sequences from GenBank

The genus name Sanghuangporus and the epithets of 14 Sanghuangporus species were used first as queries to search GenBank. Meanwhile, the reliable sequences of 14 Sanghuangporus species (Zhou et al. 2020) were used as queries to perform BLAST searches in GenBank. The cut-off value of similarity for the resulting sequences was set as 95%. All the ITS sequences matching these queries that had been deposited until 30 April 2020 were retrieved from GenBank (Table 1). In addition, recently published papers related to the taxonomy of Sanghuangporus were checked for supplementary information on collections generating these sequences (Wu et al. 2012a, 2019b; Zhou and Qin 2012; Tian et al. 2013; Ghobad-Nejhad 2015; Tomšovský 2015; Han et al. 2016; Zhou et al. 2016; Zhu et al. 2019; Huo et al. 2020; Shao et al. 2020).

Phylogenetic analyses

Two datasets of ITS sequences were assembled, one consisting of all sequences recovered from searches of GenBank and newly generated sequences, and the other consisting of the subset of sequences originating from material identified by taxonomists. The datasets were separately aligned using MAFFT 7.110 (Katoh and Standley 2013) under the G-INS-i option (Katoh et al. 2005). All resulting alignments are deposited in TreeBASE (http://www.treebase.org; accession number S26272). jModelTest (Guindon and Gascuel 2003; Posada 2008) was used to estimate the best-fit evolutionary model for each alignment with calculations made under the corrected Akaike information criterion. Following the estimated models, Maximum Likelihood (ML) and Bayesian Inference (BI) algorithms were used to construct midpoint-rooted trees for the alignments. The ML algorithm was performed using raxmlGUI 2.0 (Stamatakis 2014; Edler et al. 2021), and the bootstrap (BS) replicates were calculated under the auto FC option (Pattengale et al. 2010). The BI algorithm was performed using MrBayes 3.2 (Ronquist et al. 2012), which employed two independent runs each with four chains and starting from random trees. Trees were sampled every 1000th generation, of which the first 25% were removed as burn-in and the other 75% were retained for constructing a 50% majority consensus tree and calculating Bayesian posterior probabilities (BPPs). Tracer 1.5 (http://tree.bio.ed.ac.uk/software/tracer/) was used to judge the convergence of the chains.

Evaluation of molecular species delimitation

Molecular species delimitation was estimated using multi-rate Poisson Tree Processes (mPTP) method (Kapli et al. 2017). The Newick tree file generated from the ML algorithm was directly uploaded to the web-service version (https://mptp.h-its.org/#/tree) with no outgroup taxon.

Evaluation of genetic distances of ITS sequences

The genetic distances of an alignment of ITS sequences were estimated using MEGA X (Kumar et al. 2018; Stecher et al. 2020). For genetic distances between and within species of Sanghuangporus, the parameters were set as follows: a BS method of variance estimation with 1000 BS replications, a p-distance substitution model including transitions and transversions, uniform rates among sites, and a pairwise deletion treatment of gaps and missing data.

Identification of diagnostic ITS sequences

Identification of diagnostic ITS sequences was according to the alignment of the ITS sequences generated using MAFFT 7.110 (Katoh and Standley 2013) under the G-INS-i option (Katoh et al. 2005); if a fragment was more than one nucleotide long and was unique for one species and not variant within this species then this fragment was identified as a potential diagnostic sequence for this species.

RESULTS

A total of 13 specimens and 18 strains were newly sequenced, and the resulting ITS sequences were submitted to GenBank (Table 1). According to our criteria, 240 ITS sequences were downloaded from GenBank, but two sequences (HQ845057 and KP974834, originally identified as Inonotus vaninii and Sanghuangporus baumii, respectively) showed unexpectedly large differences from other sequences of Sanghuangporus by BLAST search, and thus were considered not to belong to the genus and were excluded from subsequent phylogenetic analyses (Table 1). Eventually, a dataset of all available 269 ITS sequences (31 newly sequenced and 238 downloaded from GenBank) from Sanghuangporus species was used to construct a preliminary phylogenetic framework for this genus. An alignment of 941 characters resulted from this dataset, and HKY + G was estimated as the best-fit evolutionary model for phylogenetic analysis. The ML search stopped after 850 bootstrap replicates. All chains in BI converged after ten million generations, which is indicated by the estimated sample sizes (ESSs) of all parameters above 500 and the potential scale reduction factors (PSRFs) close to 1.000. The ML and BI algorithms generated nearly congruent topologies in the main lineages (Additional file 1: Tree S1, Additional file 2: Tree S2). Therefore, only the topology from the ML algorithm is visualized in a circle form here; the midpoint-rooted tree recovered 13 species and four undescribed lineages of Sanghuangporus (Fig. 1). The one species gap compared with the 14 accepted species is a result of collections previously identified as S. quercicola and S. toxicodendri (this species is represented by collections Wu 1805–2, Wu 1805–3, Wu 1805–5, Wu 1807–2, Wu 1807–3 and Wu 1807–4) nesting within a single clade (Fig. 1). Of the 13 recovered species of Sanghuangporus, the clades of S. lonicericola and S. sanghuang did not receive good statistical support, the clade of S. alpinus was strongly supported just by the BI algorithm, and the other species were all strongly supported by both the ML and the BI algorithms (Additional file 1: Tree S1, Additional file 2: Tree S2). Sanghuangporus microcystideus merged with S. sp. 1 in the tree inferred from the ML algorithm (Fig. 1, Additional file 1: Tree S1), but was separated from S. sp. 1 in the BI tree (Additional file 2: Tree S2). The relationship between S. microcystideus and S. sp. 1 is still not clear, so we tentatively treat the specimen O 915609 as the single representative of S. microcystideus. One undescribed lineage including seven collections BZ-2029, BZ-2030, Dai 13360, LWZ 20190722–18, SFC 970527–1, Wu 0910–54 and Yuan 2444 showed a close relationship with S. baumii (Fig. 1).

Fig. 1
figure 1

The phylogenetic tree inferred from 269 ITS sequences. The topology was generated from the maximum likelihood algorithm. The tips in green represent mislabeled specimens, while those in red represent mislabeled strains

In GenBank, species names from 10 out of 77 phylogenetically analyzed specimens were misapplied (tips labeled in green in Fig. 1), while those from 134 out of 192 phylogenetically analyzed strains were wrongly identified to species level (tips labeled in red in Fig. 1). Furthermore, two ITS sequences (HQ845057 and KP974834) of strains labeled as species of Sanghuangporus were extremely deviant and did not belong to the genus (Table 1). Most of these errors came from submissions by non-taxonomists. Therefore, to circumscribe species in Sanghuangporus, we selected the ITS sequences submitted to GenBank by taxonomists for a new round of phylogenetic analysis (Table 1). The new dataset included 122 ITS sequences and resulted in an alignment of 871 characters with HKY + I + G as the best-fit evolutionary model. The ML search stopped after 450 bootstrap replicates. All chains in BI converged after four million generations, which is indicated by the ESSs of all parameters above 1000 and the PSRFs close to 1.000. The ML and BI algorithms generated nearly congruent topologies in the main lineages, and so only the midpoint-rooted ML tree is presented along with the BPPs at the nodes (Fig. 2). As in Fig. 1, this tree also recovered 13 species of Sanghuangporus with S. quercicola and S. toxicodendri nested within a single clade (Fig. 2). Among these 13 species, the clade of S. lonicericola was still not strongly supported, and the clades of S. alpinus and S. sanghuang were moderately supported from the ML algorithm and fully supported from the BI algorithm, while the clades of all other species received strong statistical support from both the ML and the BI algorithms (Fig. 2). Moreover, in the seven collections of the undescribed lineage close to S. baumii in Fig. 1, four were sampled in the new dataset, and the independence of these four collections and their affinity to S. baumii were also strongly supported (Fig. 2). Therefore, this undescribed lineage is described as a new species, S. subbaumii, below.

Fig. 2
figure 2

The phylogenetic tree inferred from ITS sequences submitted by taxonomists. The topology was generated from the maximum likelihood algorithm, and bootstrap values and Bayesian posterior probabilities simultaneously above 50% and 0.8, respectively, are presented at the nodes

Molecular species delimitation was estimated on the tree generated from the new dataset with 122 selected ITS sequences. The mPTP method supported the independence of 11 species, while Sanghuangporus alpinus, S. lonicerinus and S. weigelae were recovered as a single species (Additional file 3: Fig. S1).

To further explore the species relationships among Sanghuangporus, the alignment with 122 selected ITS sequences underwent a genetic distance analysis. The ranges of the within and between species genetic distances are mostly non-overlapping (Additional file 4: Table S1). Sanghuangporus microcystideus and S. pilatii, each represented by a single collection, were excluded from the within species analysis. Regarding other species of Sanghuangporus, the genetic distances within S. vaninii, S. weirianus and S. zonatus were 0–1.72%, 2.68% and 0–1.71%, respectively, whereas those within other species were no more than 1.30% and as low as 0.00% within S. ligneus (Additional file 4: Table S1). Regarding the genetic distances between species, all were above 1.30% except that those between S. alpinus and S. lonicerinus, and S. baumii and S. subbaumii were 1.03–2.86% and 1.19–3.07%, respectively. Across all pairwise comparisons between species, most (84 of 91) had distances above the maximum within species distance of 2.68% (Additional file 4: Table S1). Furthermore, distances between S. microcystideus and all other species were more than 8.90% and those between S. pilatii and all other species were more than 2.69% (Additional file 4: Table S1).

Based on an integrative taxonomic approach, 14 species of Sanghuangporus are accepted here. Their taxonomic information and reliable ITS sequences (from holotypes where possible) are provided below. Regarding S. baumii, S. lonicericola, S. lonicerinus, S. microcystideus, S. pilatii, S. vaninii, and S. weirianus, their holotypes were too old (50 years old or more) and so were unlikely to be successfully sequenced. Moreover, certain institutions did not make holotypes available for sequencing. Therefore, we use ITS sequences from other reference collections as reliable ITS sequences for those species.

Fifty-four ITS sequences of S. baumii, S. sanghuang and S. vaninii, the most common species in medicinal studies and products (Zhou et al. 2020), were further retrieved from the dataset with 122 selected sequences. These 54 sequences were realigned and the alignment is presented with shaded background (Additional file 5: Fig. S2). From this alignment, ten potential diagnostic sequences with two to six nucleotide differences were identified for HRCA to differentiate species: two for S. baumii, two for S. sanghuang and six for S. vaninii (Additional file 5: Fig. S2, Table 2).

Table 2 Diagnostic sequences with potential for discriminating Sanghuangporus baumii, S. sanghuang, and S. vaninii using Hyperbranched Rolling Circle Amplification. Label and position in alignment are as in Additional file 5: Fig. S2

TAXONOMY

Sanghuangporus alpinus (Y.C. Dai & X.M. Tian) L.W. Zhou & Y.C. Dai, Fungal Diversity 77: 340 (2016).

Basionym: Inonotus alpinus Y.C. Dai & X.M. Tian, Fungal Diversity 58: 162 (2013).

Type: China: Tibet: Linzhi County, Lulang, on living angiosperm tree, 24 Sept. 2010, B.K. Cui, Cui 9658 (BJFC – holotype).

ITS barcoding sequence: JQ860310 (from holotype).

Sanghuangporus baumii (Pilát) L.W. Zhou & Y.C. Dai, Fungal Diversity 77: 340 (2016).

Basionym: Phellinus baumii Pilát, Bull. trimest. Soc. mycol. Fr. 48: 25 (1932).

Synonym: Inonotus baumii (Pilát) T. Wagner & M. Fisch., Mycologia 94: 1009 (2002).

Type: Russia: Primorsky Krai: Vladivostok, on trunk of Syringae, 5 June 1928, M.K. Ziling 267 (PRM 189012 – holotype).

Reference collection: China: Heilongjiang: Yichun, Fenglin nature reserve, on living trunk of Syringa, 8 Sept. 2002, Y.C. Dai, Dai 3683 (IFP)

ITS barcoding sequence: JN642567 (from the reference collection cited above, proposed by Zhou et al. (2020) and accepted here).

Sanghuangporus ligneus Ghob.-Nejh., Mycol. Progr. 14(90): 2 (2015).

Type: Iran: East Azerbaijan: Khoda-Afarin, Kalaleh-Eslami, Darana, deciduous forest with Quercus macranthera, Lonicera, Cornus mas, and Crataegus, on stem of living Lonicera caucasica, 10 May 2008, M. Ghobad-Nejhad, Ghobad-Nejhad 1152 (ICH – holotype).

ITS barcoding sequence: KR073081 (from holotype).

Sanghuangporus lonicericola (Parmasto) L.W. Zhou & Y.C. Dai, Fungal Diversity 77: 340 (2016).

Basionym: Phellinus lonicericola Parmasto, Folia cryptog. Estonica 38: 59 (2001).

Synonym: Inonotus lonicericola (Parmasto) Y.C. Dai, Fungal Diversity 45: 276 (2010).

Type: Russia: Primorsky Krai: Lazovsky Nature Reserve, Petrov island, on trunk of Lonicera ruprechtiana in Taxus mixed forest, 2 Sept. 1961, E. Parmasto (TAA-M 013933 – holotype).

Reference collection: China: Heilongjiang: Ningan County, Jingpohu National Scenic Area, on living trunk of Lonicera, 8 Sept. 2007, Y.C. Dai, Dai 8376 (IFP)

ITS barcoding sequence: JQ860308 (from the reference collection cited above, proposed by Zhou et al. (2020) and accepted here).

Sanghuangporus lonicerinus (Bondartsev) Sheng H. Wu et al., Fungal Diversity 77: 340 (2016).

Basionym: Fomes lonicerinus Bondartsev, Acta Inst. Bot. Acad. Sci. USSR Plant. Crypt., Ser. II: no. 500 (1935).

Synonyms: Phellinus lonicerinus (Bondartsev) Bondartsev & Singer, Annls mycol. 39: 56 (1941).

Cryptoderma lonicerinum (Bondartsev) Imazeki, Bull. Tokyo Sci. Mus. 6: 107 (1943).

Porodaedalea lonicerina (Bondartsev) Imazeki, Col. Ill. Mushrooms Japan, 2: 191 (1989).

Inonotus lonicerinus (Bondartsev) Sheng H. Wu et al., Bot. Studies (Taipei) 53: 140 (2012).

Type: Uzbekistan: Samarkand: Sarymat, on trunk of Lonicera tatarica, 1926, E. Czerniakowsk (LE 22512 – lectotype designated by Bondartsev 1953).

Reference collection: Turkmenistan: Bakharden: Bakharden, Arvaz, Montes Kopet-dagh, on Lonicera, 17 Oct. 1971, E. Parmasto (TAA 55428)

ITS barcoding sequence: JN642575 (from the reference collection cited above, proposed by Zhou et al. (2020) and accepted here).

Sanghuangporus microcystideus (Har. & Pat.) L.W. Zhou & Y.C. Dai, Fungal Diversity 77: 340 (2016).

Basionym: Phellinus microcystideus Har. & Pat., Bull. Mus. natn. Hist. nat., Paris 15: 90 (1909).

Synonym: Fomes microcystideus (Har. & Pat.) Sacc. & Trotter, Syll. Fung. 21: 286 (1912).

Type: Congo: Moyen Oubangui: Grande Forêt, M.A. Chevalier 11431 (FH – holotype).

Reference collection: Tanzania: Arusha: Arusha National Park, Mount Meru, on trunk of Olea africana, 18 Feb. 1976, R. Harjula (O 915609)

ITS barcoding sequence: KP030787 (from the reference collection cited above, proposed by Zhou et al. (2020) and accepted here).

Sanghuangporus pilatii (Černý) Tomšovský, Phytotaxa 239: 84 (2015).

Basionym: Phellinus pilatii Černý, Česká Mykol. 22(1): 2 (1968).

Synonym: Porodaedalea pilatii (Černý) Fiasson & Niemelä, Karstenia 24(1): 26 (1984).

Type: Czech Republic: Břeclav: Tvrdonice, 8 Oct. 1955, A. Černý (PRM 628393 – holotype).

Reference collection: Czech Republic: Břeclav: Nové Mlýny, Křivé jezero National Nature Reserve, on Populus alba, 22 Oct. 2011, M. Tomšovský 41/2011 (BRNM 771989)

ITS barcoding sequence: KT428764 (from the reference collection cited above, proposed by Zhou et al. (2020) and accepted here).

Sanghuangporus quercicola Lin Zhu & B.K. Cui, Phytotaxa 311: 271 (2017).

Synonym: Sanghuangporus toxicodendri Sheng H. Wu et al., MycoKeys 57: 106 (2019).

Type: China: Henan: Neixiang County, Baotianman Nature Reserve, on dead tree of Quercus, 25 Aug. 2006, J. Li, Li 1149 (BJFC – holotype).

ITS barcoding sequence: KY328312 (from holotype).

Sanghuangporus sanghuang (Sheng H. Wu et al.) Sheng H. Wu et al., Fungal Diversity 77: 340 (2016).

Basionym: Inonotus sanghuang Sheng H. Wu et al., Bot. Studies (Taipei) 53: 140 (2012).

Type: China: Jilin: Baishan City, on Morus sp., Mar. 2009, S.H. Wu, Wu 0903–1 (TNM – holotype).

ITS barcoding sequence: JN794061 (from holotype).

Sanghuangporus subbaumii Shan Shen, Y.C. Dai & L.W. Zhou, sp. nov. (Figs. 3 and 4).

Fig. 3
figure 3

Basidiomes of Sanghuangporus subbaumii in situ. a Dai 13360 (holotype). b LWZ 20190722–18 (paratype). Bars: 2 cm

Fig. 4
figure 4

Microscopic structures of Sanghuangporus subbaumii (drawn from Dai 13360, holotype). a Basidiospores. b basidia and basidioles. c cystidioles. d hymenial setae. e hyphae from trama. f hyphae from context. Bars: a = 5 μm, b–e = 10 μm

MycoBank MB838235.

Etymology: subbaumii (Lat.), refers to the similarity to Sanghuangporus baumii.

Diagnosis: Differing from S. baumii in having resupinate, effused-reflexed to pileate basidiomes, acute pileal margin and longer hymenial setae (> 20 μm in length).

Type: China: Shanxi: Jiaocheng County, Pangquangou Nature Reserve, on fallen trunk of Prunus sp., 10 Aug. 2013, Y.C. Dai, Dai 13360 (BJFC – holotype; HMAS 281653 – isotype).

Description: Basidiomes perennial, resupinate, effused-reflexed to pileate, without odor or taste and hard corky when fresh, woody hard when dry; to 20 cm long and 5 cm wide when resupinate. Pilei dimidiate, ungulate in section, projecting to 3.5 cm wide, 6 cm long and 4 cm thick at base. Pileal surface dark brown and velutinate when juvenile, mouse-grey to black, glabrous and cracked with age, concentrically zonate and narrowly sulcate; margin yellow brown, acute. Pore surface yellowish brown, glancing; sterile margin distinct, yellowish; pores angular to circular, 5–7 per mm; dissepiments thin, entire. Context yellowish brown to dark brown, woody hard, to 3.5 cm thick. Tubes yellowish brown, darker than pore surface, woody hard, to 0.5 cm long.

Hyphal system monomitic in context, dimitic in trama; generative hyphae simple septate; tissue darkening but otherwise unchanged in KOH. Context generative hyphae occasionally slightly thick-walled with a wide lumen and yellowish, mostly thick-walled with a narrow lumen and yellowish brown, unbranched, frequently septate, more or less regularly arranged, 3.5–4 μm diam. Tubes generative hyphae thin to slightly thick-walled, hyaline, occasionally branched, frequently septate, 3–4.5 μm diam; skeletal hyphae dominant, thick-walled with a narrow lumen, yellowish brown, unbranched, rarely septate, subparallel along the tubes, 2.2–3.7 μm diam. Hymenial setae frequent in the mature hymenium, subulate to ventricose, dark brown, thick-walled, 20–35 × 7–12 μm. Cystidioles subulate, with narrow and tapering apex, hyaline, 15–20 × 4–6 μm. Basidia barrel-shaped to broadly clavate, with four sterigmata and a simple septum at the base, hyaline, 20–25 × 7–9 μm; basidioles in shape similar to basidia, but slightly smaller. Basidiospores broadly ellipsoid to subglobose, yellowish, slightly thick-walled, smooth, non-amyloid, non-dextrinoid, moderately cyanophilous, (3.8–)4–4.9(− 5.2) × 3.1–3.8(− 3.9) μm, L = 4.35 μm, W = 3.41 μm, Q = 1.24–1.31 (n = 60/2).

Notes: Sanghuangporus subbaumii mostly resembles S. baumii, but the latter species differs in having pileate basidiomes always, obtuse pileal margin and shorter hymenial setae (< 20 μm in length; Dai 2010). The resupinate to pileate basidiomes make S. subbaumii similar to S. vaninii, but S. vaninii lacks cystidioles and has a thin black zone separating heterogeneous context (Dai 2010).

ITS barcoding sequence: MT348580 (from holotype).

Additional specimen examined: China: Beijing: Shangfangshan Forest Park, on fallen angiosperm trunk, 22 July 2019, L.W. Zhou, LWZ 20190722–18 (HMAS 281654).

Sanghuangporus vaninii (Ljub.) L.W. Zhou & Y.C. Dai, Fungal Diversity 77: 340 (2016).

Basionym: Phellinus vaninii Ljub., Bot. Mater. 15: 115 (1962).

Synonym: Inonotus vaninii (Ljub.) T. Wagner & M. Fisch., Mycologia 94: 1009 (2002).

Type: Russia: Primorsky Krai: Shkotovsky District, watershed of the Maykhe river, Maykhinsky forestry, Verkhne-Maykhinskaya forest area, Peyshula, quarter 119, in valley of pine-broadleaved forest, on dried aspen tree, 14 Aug. 1951, L.V. Lyubarskiy (LE 22523 – holotype).

Reference collection: China: Jilin: Antu County, Changbaishan, on fallen trunk of Populus davidiana, 26 Aug. 2005, Y.C. Dai, Dai 7011 (IFP)

ITS barcoding sequence: JN642591 (from the reference collection cited above, proposed by Zhou et al. (2020) and accepted here).

Sanghuangporus weigelae (T. Hatt. & Sheng H. Wu) Sheng H. Wu et al., Fungal Diversity 77: 340 (2016).

Basionym: Inonotus weigelae T. Hatt. & Sheng H. Wu, Bot. Studies (Taipei) 53: 143 (2012).

Synonym: Inonotus tenuicontextus L.W. Zhou & W.M. Qin, Mycol. Progr. 11: 793 (2012).

Type: Japan: Nagano: Chino, Minoto, on Weigela coraeensis, 19 Sept. 1993, T. Hattori, F16899 (TFM – holotype).

ITS barcoding sequence: JN642596 (from holotype).

Sanghuangporus weirianus (Bres.) L.W. Zhou & Y.C. Dai, Fungal Diversity 77: 340 (2016).

Basionym: Fomes weirianus Bres., Stud. Trent., Classe II, Sci. Nat. Econ. 7(1): 5 (1926).

Synonyms: Phellinus weirianus (Bres.) Gilb., J. Ariz. Acad. Sci. 7: 137 (1972).

Inonotus weirianus (Bres.) T. Wagner & M. Fisch., Mycologia 94: 1009 (2002).

Type: USA: New Mexico: on trunk of Juglans rupestris, 25 Oct. 1911, G.G. Hedgcock & W.H. Long (BPI 235278 – holotype).

Reference collection: USA: Arizona: on Juglans major, 27 Aug. 1967, R.L. Gilbertson 6975-S (IMSNU 32021)

ITS barcoding sequence: AF110989 (from the reference collection cited above, proposed by Zhou et al. (2020) and accepted here).

Sanghuangporus zonatus (Y.C. Dai & X.M. Tian) L.W. Zhou & Y.C. Dai, Fungal Diversity 77: 341 (2016).

Basionym: Inonotus zonatus Y.C. Dai & X.M. Tian, Fungal Diversity 58: 165 (2013).

Type: China: Hainan: Jianfengling Nature Reserve, on living angiosperm tree, 11 May 2009, B.K. Cui, Cui 6631 (BJFC – holotype).

ITS barcoding sequence: JQ860305 (from holotype).

DISCUSSION

In this study, we summarized all available ITS barcoding sequences bearing the name “Sanghuang” in GenBank. A total of 271 ITS sequences related to “Sanghuang”, including 31 newly generated sequences from this study, were analyzed. In association with previous information of morphology, hosts, and multilocus-based phylogeny, 14 species are accepted as members of Sanghuangporus including the new species S. subbaumii described herein. We also synonymize S. toxicodendri under S. quercicola.

Sanghuangporus subbaumii has a phylogenetically close relationship to S. baumii; however, these two species form two distinct lineages with strong support (Additional file 1: Tree S1, Additional file 2: Tree S2, Fig. 2). Moreover, S. subbaumii and S. baumii were also estimated as two independent species using the mPTP method (Additional file 3: Fig. S1), and for ITS the interspecific distance is 1.19–3.07%, generally above the cut-off value of interspecific distances (1.30%) within Sanghuangporus (Additional file 4: Table S1). Besides molecular evidence, morphological differences between these two species are also clear. Geographically, S. subbaumii is only known from North China, whereas Chinese collections of S. baumii are distributed in north-east China (Table 1).

Sanghuangporus toxicodendri was recently described from specimens collected from Toxicodendron sp. in Hubei, central China (Wu et al. 2019b) and resembles S. quercicola, another species originally described from central China (Zhu et al. 2017). However, in the publication introducing S. toxicodendri (Wu et al. 2019b) the separation from S. quercicola was not well-supported phylogenetically. Moreover, the morphological differences between these two species are slight (such as for basidiospore length) or involve variable characters that do not have taxonomic signal (such as the surface color of the pileal margin) (Zhu et al. 2017; Wu et al. 2019b). In the current phylogenetic analyses, the six specimens of S. toxicodendri, three specimens of S. quercicola and four additional collections merged in a fully supported clade (Additional file 1: Tree S1, Additional file 2: Tree S2, Fig. 2). The mPTP-based estimation of species delimitation also treated S. toxicodendri and S. quercicola as a single species (Additional file 3: Fig. S1) and the intraspecific distances among ITS sequences under both names were 0–1.11%, well below the threshold of 1.30% (Additional file 4: Table S1). Therefore, S. toxicodendri and S. quercicola are considered conspecific, and S. quercicola has priority by publication date over S. toxicodendri.

The clade of S. lonicericola was present but not well-supported in our phylogenetic analyses (Additional file 1: Tree S1, Additional file 2: Tree S2, Fig. 2). Similarly, the clades of S. alpinus and S. sanghuang were not strongly supported by the ML algorithm (Fig. 2). For S. lonicericola and S. alpinus, despite the lack of support in one or both analyses, each formed a distinct clade, and for both species distances to other species were above the threshold of 1.30% (S. lonicericola minimum 2.19% and S. sanghuang minimum 2.90%; Additional file 4: Table S1). In addition, S. alpinus, S. lonicerinus, and S. weigelae, even though forming three independent lineages, were considered conspecific by the mPTP method (Additional file 3: Fig. S1). However, the interspecific distances for ITS between S. weigelae and each of S. alpinus and S. lonicerinus are above the cut-off value of interspecific distances (1.30%) within Sanghuangporus (Additional file 4: Table S1). Regarding the pair of S. alpinus and S. lonicerinus, for ITS the between species distance (1.03–2.86%) was generally above the intraspecific distances within either species (0–1.08% and 0–1.18%, respectively; Additional file 4: Table S1). Moreover, the monophyly of S. alpinus was strongly supported by the BI algorithm and that of S. lonicerinus was strongly supported by both the ML and the BI algorithms (Fig. 2). Besides, morphological delimitations among these five species are stable (Wu et al. 2012a; Tian et al. 2013; Zhou et al. 2016). Taking all this into account, we accept S. alpinus, S. lonicericola, S. lonicerinus, S. sanghuang, and S. weigelae as five independent species.

Sanghuangporus vaninii, S. weirianus, and S. zonatus are the only three species with intraspecific ITS distances of more than 1.30% (0–1.72%, 2.68% and 0–1.71%, respectively; Additional file 4: Table S1). However, they all received strong support as independent species (Additional file 1: Tree S1, Additional file 2: Tree S2, Fig. 2, Additional file 3: Fig. S1). As one of the most commonly cultivated species, several cultivars of S. vaninii were included in the evaluation of genetic distances of ITS sequences (Zhou et al. 2020; Table 1). The procedure of cultivation with continuous passage culture can dramatically accelerate the accumulation of genetic variation, which may result in the higher intraspecific ITS difference in S. vaninii. Noteworthily, branch lengths of the only two available collections of S. weirianus were markedly different even though the two strains were from the same original isolate (Fig. 2). Regarding S. zonatus, two collections from Hainan, South China grouped together with full statistical support, and then formed a fully supported clade with a collection from Yunnan, Southwest China (Table 1, Figs. 1 and 2). Both S. weirianus and S. zonatus are poorly collected species, and a more comprehensive sampling of these two species in phylogenetic analyses will further clarify their intraspecific relationships. For now, we tentatively accept them as monophyletic species.

A study by Nilsson et al. (2006) revealed that about 10–21% of 51,000 fungal ITS sequences available at that time in the International Nucleotide Sequence Databases were annotated with incorrect taxonomic information. More recently, this proportion has increased to almost 30% (Hofstetter et al. 2019). Regarding “Sanghuang”, more than half (or say 146) of the ITS sequences labeled as such, were found to be mislabeled, implying that the proportion of incorrectly labeled ITS sequences for “Sanghuang” is much higher than the average proportion for all fungal groups. This phenomenon may be attributable to the medicinal properties of “Sanghuang”, which attracts much more attention from non-taxonomists who submit ITS sequences to GenBank. Consequently, the numerous errors result in chaos with BLAST searches, especially for non-taxonomists. Although the RefSeq Targeted Loci (RTL) database has been initiated for fungal ITS sequences from type collections (Schoch et al. 2014), only two species of Sanghuangporus, viz. S. alpinus and S. zonatus were reannotated and deposited under accession numbers of NR_158887 and NR_166366. Actually, ITS sequences from six holotypes of accepted Sanghuangporus species are available in GenBank. This number increases to eight, if two synonyms of other species of Sanghuangporus, viz. Inonotus tenuicontextus and S. toxicodendri are considered. In UNITE (Nilsson et al. 2019), tens of species hypotheses belonging to Sanghuangporus are available under various threshold values at species level; however, not all accepted species of Sanghuangporus (such as S. ligneus, S. pilatii, and S. quercicola) are referred to and the reference sequences for some species hypotheses are not always those from holotypes. Moreover, both RTL and UNITE are not familiar to mycologists working on medicinal studies and government officers in charge of the policy of medicinal fungi, who normally take the first hit of a BLAST search in GenBank as the species name. Therefore, the accuracy of ITS sequences of “Sanghuang” in GenBank is crucial for medicinal studies and commercial development of this fungal genus.

Compared with specimens, many more mislabeled ITS sequences of Sanghuangporus came from cultured strains, and most of those sequences were submitted by non-taxonomists. A typical case is the recent paper on genome sequencing of “Sanghuang” that also submitted six ITS sequences to GenBank (Shao et al. 2020). In GenBank, all these six sequences were labeled as Inonotus sp. rather than species of Sanghuangporus (MN242716–MN242721), while the six strains generating these sequences were named as S. sanghuang (Shao et al. 2020). However, five of the six strains, including the one (labeled as KangNeng) subjected to genome sequencing, are actually S. vaninii (Fig. 1, Zhou et al. 2020); i.e. five out of six strains were wrongly identified to species level. Therefore, this species misidentification means that the whole genome sequence of “Sanghuang” may be misapplied in future studies. Shao et al. (2020) also stated that these six strains are commercially cultivated, which further results in the name chaos for commercial products of “Sanghuang”. Another publication on genome sequencing identified the genome sequenced strain S12 as Phellinus gilvus according to ITS barcoding region (Huo et al. 2020). However, the corresponding ITS sequence (MT275660) annotated as Fuscoporia gilva in GenBank represents S. vaninii (Fig. 1, Zhou et al. 2020). Another case is a paper devoted to the species identity of “Sanghuang” strains (Han et al. 2016). Thirty strains deposited in the Agricultural Sciences Institute culture collection (Mushroom Research Division, Rural Development Administration, Republic of Korea) were correctly identified as S. vaninii and S. sanghuang according to an ITS-based phylogenetic analysis; however, unfortunately, most of these ITS sequences were mislabeled when being submitted to GenBank.

Ten mislabeled ITS sequences found in the current study came from basidiomes. These errors were caused mainly by taxonomic revisions of certain species. Six sequences of specimens Wu 1805–2, Wu 1805–3, Wu 1805–5, Wu 1807–2, Wu 1807–3 and Wu 1807–4 that were originally labeled as Sanghuangporus sp. but later cited under S. toxicodendri by Wu et al. (2019b) are accepted to represent S. quercicola. Yuan 2444, previously considered as S. baumii, was nested within the lineage segregated from S. baumii as a new species S. subbaumii (Figs. 1 and 2, Additional file 3: Fig. S1). Consequently, the ITS sequence of Yuan 2444 (JX069836) is corrected to S. subbaumii (Table 1). Another mislabeled sequence was generated from a specimen originally described as Inonotus tenuicontextus (Zhou and Qin 2012). Although this species was published online earlier than Inonotus weigelae (basionym of S. weigelae; Wu et al. 2012a; Tian et al. 2013), its online date is before 1 January 2012 and thus the name was not effectively published online according to Art. 29.1 of the ICNafp (Turland et al. 2018). Inonotus tenuicontextus was then treated as a later synonym of I. weigelae (Tian et al. 2013). Therefore, this mislabeled sequence is accepted to represent S. weigelae (Table 1).

Although intact mature basidiomes of “Sanghuang” are not difficult to identify to species level morphologically and in a short time by taxonomists working on this group, most of the commercial products are small pieces or even powders. Normally, it is impossible to rapidly determine which species those commercial products represent. As for other traditional medicinal mushrooms (Raja et al. 2017), species names of Sanghuangporus are sometimes misapplied to certain products of “Sanghuang” (Shao et al. 2020). This confused situation to some extent restricts the commercial development of “Sanghuang” (Zhou 2020). Therefore, to standardize the “Sanghuang” industry, ten reference sequences are provided for HRCA based on the accurate boundaries among three commonly studied and cultivated species, viz. S. baumii, S. sanghuang, and S. vaninii (Lin et al. 2017; Zhou et al. 2020). HRCA is an isothermal amplification approach and thus provides a rapid, simple and low-cost detection of specific nucleic acid sequences (Nilsson et al. 1994; Lizardi et al. 1998) even for single nucleotide differences (Nilsson et al. 1997). This approach has been widely used for the clinical detection of human pathogenic microfungi (Zhou et al. 2008; Trilles et al. 2014; Rodrigues et al. 2015) and, recently, was also reported for the rapid detection of poisonous macrofungi (He et al. 2019a, 2019b). Regarding lethal Amanita species, nucleotide differences greater than two allowed species identification using the α-amanitin gene (He et al. 2019a). Here, for Sanghuangporus a set of candidates for future testing is provided that have diagnostic sequences containing between two and six nucleotide differences.

CONCLUSION

In order to promote medicinal studies and industrial development, the ITS barcoding region of Sanghuangporus species is here comprehensively analyzed to enable accurate species identification. Firstly, the ITS region is confirmed as an effective barcode in Sanghuangporus. Secondly, the names of all available ITS sequences in GenBank related to “Sanghuang” are carefully revised and where necessary corrected. Thirdly, the intraspecific ITS difference for each species of Sanghuangporus is evaluated to be up to 1.30% (except S. vaninii, S. weirianus, and S. zonatus), while the interspecific ITS difference is above 1.30% (except between S. alpinus and S. lonicerinus, and S. baumii and S. subbaumii). This provides a practical cut-off value for BLAST search-based species identification. Finally, ten potential diagnostic sequences are provided for HRCA assay to rapidly differentiate the three commonly studied and cultivated species, viz. S. baumii, S. sanghuang, and S. vaninii. As a follow up, we will suggest reannotation of ITS sequences related to “Sanghuang” to the GenBank administrators, especially to ensure that sequences from holotypes and reference collections for each species of Sanghuangporus are designated as such. Further, we will liaise with UNITE to ensure that appropriate reference sequences are designated for UNITE species hypotheses within Sanghuangporus.

Availability of data and materials

The materials are available as Additional files 1, 2, 3, 4 and 5. All sequence data generated for this study can be accessed via GenBank: https://www.ncbi.nlm.nih.gov/genbank/. Alignments are available at TreeBase (ID: 26272).

Abbreviations

BI:

Bayesian inference

BPP:

Bayesian posterior probability

CB:

Cotton Blue

CTAB:

Cetyl-trimethyl-ammonium bromide

IKI:

Melzer’s reagent

ITS:

Nuclear ribosomal internal transcribed spacer

KOH:

5% potassium hydroxide

ML:

Maximum likelihood

mPTP:

Multi-rate Poisson Tree Processes

PCR:

Polymerase chain reaction

RTL:

RefSeq Targeted Loci

References

Download references

Acknowledgements

We appreciate Genevieve Gates (Tasmanian Institute of Agriculture, Australia) for help in improving this manuscript, editors and anonymous reviewers for providing constructive comments on the earlier versions of the manuscript, and curators from FH, LE and TAAM for forwarding information of holotypes. Yan Yang and He-Nan Zhang (Institute of Edible Fungi, Shanghai Academy of Agricultural Sciences, China) are thanked for kindly providing cultivated strains for sequencing. Yu-Cheng Dai (BJFC) and Hai-Sheng Yuan (IFP) are thanked for kindly forwarding specimens and strains as loans for sequencing.

Adherence to national and international regulations

Not applicable.

Funding

The research was financed by the National Natural Science Foundation of China (No. 31970012), Youth Innovation Promotion Association of the Chinese Academy of Sciences (No. 2017240), and Biological Resources Programme, Chinese Academy of Sciences (KFJ-BRP-017-12).

Author information

Authors and Affiliations

Authors

Contributions

SS, S-LL and L-WZ retrieved and analyzed all data. J-HJ prepared fungal samples and performed molecular sequencing. L-WZ conceived the work and wrote the manuscript. All authors approved the manuscript.

Corresponding author

Correspondence to Li-Wei Zhou.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

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

Additional file 1: Tree S1.

The phylogenetic tree inferred from 269 ITS sequences. The topology was generated from the maximum likelihood algorithm and bootstrap values are presented at the nodes.

Additional file 2: Tree S2.

The phylogenetic tree inferred from 269 ITS sequences. The topology was generated from the Bayesian inference algorithm and Bayesian posterior probabilities are presented at the nodes.

Additional file 3: Figure S1.

Molecular species delimitation estimated from the Newick tree file of Fig. 2 using multi-rate Poisson Tree Processes method. The continuous red branches represent a single species.

Additional file 4: Table S1.

Genetic distances of ITS sequences between and within species of Sanghuangporus.

Additional file 5: Figure S2.

The alignment of Sanghuangporus baumii, S. sanghuang and S. vaninii generated from ITS sequences submitted by taxonomists. Ten potential diagnostic sequences for Hyperbranched Rolling Circle Amplification are labeled in capital letters.

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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Shen, S., Liu, SL., Jiang, JH. et al. Addressing widespread misidentifications of traditional medicinal mushrooms in Sanghuangporus (Basidiomycota) through ITS barcoding and designation of reference sequences. IMA Fungus 12, 10 (2021). https://doi.org/10.1186/s43008-021-00059-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s43008-021-00059-x

Keywords