Comparative analyses of Flammulina filiformis mitochondrial genomes reveal high length polymorphism in intergenic regions and multiple intron gain/loss in cox1

The golden-needle mushroom Flammulina filiformis is one of the bulk mushroom products in the world. This study obtained complete mitogenomes of 44 wild isolates collected from nine provinces and two artificially bred cultivars of F. filiformis, together with three Flammulina rossica isolates and one Flammulina fennae isolate for comparison. The mitogenome of F. filiformis ranged from 83,540 bp to 90,938 bp, consisting of 14 conserved protein-coding genes (PCGs), two rRNA genes, and 25 tRNA genes. To the best of our knowledge, it contained the highest proportion of intergenic regions compared to the other known Basidiomycota mitogenomes. Introns and intergenic regions were two major contributing factors to the total size of the F. filiformis mitogenome. The conserved PCG cox3 is located in an intron of another conserved PCG, nad5. This is a unique phenomenon in all known fungal mitogenomes. Gain/loss of introns was observed in cox1, nad5, and rnl. Length polymorphism was widely observed in intergenic regions. Accordingly, primers were designed as useful markers for rapid identi-fication of F. filiformis isolates with differentiated mitogenomes. Our findings provide a basis for further studies related to variety identification and population genetics of this economically important mushroom.