Background: Ophiocordyceps sinensis (syn. Cordyceps sinensis), endemic to alpine regions on the Tibetan plateau, is one of the most valuable medicinal fungi in the world. Huge commercial demand has led to excessive harvest and a dramatic decline in its numbers. The diversity of terrains and climates on the Tibetan Plateau and the broad insect host range (more than 50 species in the family Hepialidae) may have resulted in substantial intraspecific genetic diversity for this fungus. The objective of this study was to evaluate the population distribution of O. sinensis from geographically diverse regions of the Tibetan Plateau based on nrDNA ITS and MAT1-2-1 gene sequences. Understanding of the genetic diversity and genesis of O. sinensis will provide important information for the evolution and conservation of this fungus. Results: Significant sequence variations in the ITS and MAT1-2-1 genes (27 and 23 informative sites, eight and seven haplotypes, respectively) were observed. Phylogenetic analysis based on ITS sequences, MAT1-2-1 sequences, or their combined data set, clustered isolates from northern regions in one clade (clade I), whereas isolates from southern regions were dispersed in all four clades (clade I-IV). Single-strand conformation polymorphism (SSCP) analyses of 2639 ITS clones from seven samples revealed 91 different SSCP patterns that were subsequently sequenced. ITS heterogeneity was found in XZ-LZ07-H1 (Nyingchi population), and 17 informative sites and five haplotypes were detected from 15 clones. The five haplotypes clustered into three clades (clade I, II, and IV). Conclusions: Significant genetic divergence in O. sinensis was observed and the genetic diversification was greater among southern isolates than that among northern isolates. The polymorphism of nrDNA ITS sequences suggested that O. sinensis spread from a center of origin (the Nyingchi District) to southern regions and subsequently to northern areas. These results suggest that southern populations are important reservoirs of genetic diversity and should be taken into account in conservation programs.

Fungal mitochondrial genes are often invaded by group I or II introns, which represent an ideal marker for understanding fungal evolution. A standard nomenclature of mitochondrial introns is needed to avoid confusion when comparing different fungal mitogenomes. Currently, there has been a standard nomenclature for introns present in rRNA genes, but there is a lack of a standard nomenclature for introns present in protein-coding genes. In this study, we propose a new nomenclature system for introns in fungal mitochondrial protein-coding genes based on (1) three-letter abbreviation of host scientific name, (2) host gene name, (3), one capital letter P (for group I introns), S (for group II introns), or U (for introns with unknown types), and (4) intron insertion site in the host gene according to the cyclosporin-producing fungus Tolypocladium inflatum. The suggested nomenclature was proved feasible by naming introns present in mitogenomes of 16 fungi of different phyla, including both basal and higher fungal lineages although minor adjustment of the nomenclature is needed to fit certain special conditions. The nomenclature also had the potential to name plant/protist/animal mitochondrial introns. We hope future studies follow the proposed nomenclature to ensure direct comparison across different studies.