Comparison of mitochondrial genomes provides insights into intron dynamics and evolution in the caterpillar fungus Cordyceps militaris

Intra-specific comparison of mitochondrial genomes can help elucidate the evolution of a species, however it has not been performed for hypocrealean fungi that form diverse symbiotic associations with other organisms. In this study, comparative analyses of three completely sequenced mitochondrial genomes of a hypocrealean fungus, Cordyceps militaris, the type species of Cordyceps genus, revealed that the introns were the main contributors to mitochondrial genome size variations among strains. Mitochondrial genes in C. militaris have been invaded by group I introns in at least eight positions. PCR assays of various C. militaris isolates showed abundant variations of intron presence/absence among strains at seven of the eight intronic loci. Although the ancestral intron pattern was inferred to contain all eight introns, loss and/or gain events occurred for seven of the eight introns. These introns invaded the C. militaris mitochondrial genome probably by horizontal transfer from other fungi, and intron insertions into intronless genes in C. militaris were accompanied by co-conversions of upstream exon sequences especially for those introns targeting protein-coding genes. We also detected phylogenetic congruence between the intron and exon trees at each individual locus, consistent with the ancestral mitochondria of C. militaris as having all eight introns. This study helps to explain the evolution of C. militaris mitochondrial genomes and will facilitate population genetic studies of this medicinally important fungus.