PacBio high-throughput multi-locus sequencing reveals high genetic diversity in mushroom-forming fungi

Multi-locus sequence data are widely used in fungal systematic and taxonomic studies to delimit species and infer evolutionary relationships. We developed and assessed the efficacy of a multi-locus pooled sequencing method using PacBio long-read high-throughput sequencing. Samples included fresh and dried voucher specimens, cultures and archival DNA extracts of Agaricomycetes with an emphasis on the order Cantharellales. Of the 283 specimens sequenced, 93.6% successfully amplified at one or more loci with a mean of 3.3 loci amplified. Our method recovered multiple sequence variants representing alleles of rDNA loci and single copy protein-coding genes rpb1, rpb2 and tef1. Within-sample genetic variation differed by locus and taxonomic group, with the greatest genetic divergence observed among sequence variants of rpb2 and tef1 from corticioid Cantharellales. Our method is a cost-effective approach for generating accurate multi-locus sequence data coupled with recovery of alleles from polymorphic samples and multi-organism specimens. These results have important implications for understanding intra-individual genomic variation among genetic loci commonly used in species delimitation of fungi.