Intraspecific mitochondrial variability was studied in ten strains of A. bisporus var. bisporus, in a strain representative of A. bisporus var. eurotetrasporus and in a strain of the closely related species Agaricus devoniensis. In A. bisporus, the cox1 gene is the richest in group I introns harboring homing endonuclease genes (heg). This study led to identify group I introns as the main source of cox1 gene polymorphism. Among the studied introns, two groups were distinguished according to the heg they contained. One group harbored heg maintained putatively functional. The other group was composed of eroded heg sequences that appeared to evolve toward their elimination. Low nucleotide substitution rates were found in both types of intronic sequences. This feature was also shared by all types of studied mitochondrial sequences, not only intronic but also genic and intergenic ones, when compared with nuclear sequences. Hence, the intraspecific evolution of A. bisporus mitochondrial genome appears characterized by both an important mobility (presence/absence) of large group I introns and by low nt substitution rates. This stringent conservation of mitochondrial sequences, when compared with their nuclear counterparts, appears irrespective of their apparent functionality and contrasts to what is widely accepted in fungal sequence evolution. This strengthens the usefulness of mtDNA sequences to get clues on intraspecific evolution.

Keywords