Abstract
Many macrofungal cryptic species remain unidentified. A possible solution is to increase the number of loci analyzed and use rigorous statistics for macrofungal species delimitation. To validate this assumption, cryptic species of the Hypholoma fasciculare complex, a group of common wood-decomposing fungi, were attempted to be delineated. Massively parallel sequencing of mitochondrial ribosomal RNA (mt_rRNA), nuclear ribosomal internal transcribed spacer (ITS) region, and 24 single-copy genes were performed for 96 specimens collected in Japan. Then, the species boundaries were inferred using comparative gene genealogies (mt_rRNA vs. ITS), Bayesian Poisson tree process (bPTP) model for the phylogeny of concatenated nuclear sequences, and analysis of molecular variance (AMOVA) for single nucleotide polymorphisms. In both the mt_rRNA and ITS phylogenies, the H. fasciculare complex was not divided into well-supported clades. Nevertheless, based on the bPTP, two mitochondrial haplotypes were inferred to represent distinct species (H. fasciculare and H. subviride). The results of AMOVA also indicated that the differentiation of nuclear loci can be explained mostly by differences between haplotype. These results suggest that it is necessary to increase the number of target loci to 20 or more and use both phylogeny-based and population genetics-based statistics for the accurate delimitation of macrofungal species.