Variation in SSRs at different genomic regions and implications for the evolution and identification of Armillaria gallica

Armillaria spp. are devastating forest pathogens. Due to its low pathogenicity and abundant genetic variation, Armillaria gallica exhibited a unique and beneficial symbiosis with Gastrodia elata, which was used as a traditional Chinese medicine. However, the variation and population structure of A. gallica populations have rarely been investigated. Hence, we analyzed the evolution and variation in simple sequence repeats (SSRs) in three Armillaria genomes: A. gallica, A. cepistipes, and A. ostoyae to assess the genetic diversity and population structure of 14 A. gallica strains. Genome analysis revealed that SSRs were more abundant in the intergenic region than the intron and exon region, as was the SSR density. Compared with other two genomes, SSR density was the lowest in exon region and largest in the intron region of A. gallica, with significant variation in genic region. There were 17 polymorphic markers in A. gallica genome was identified, with 26.7% in genic region, which is higher than that of 18.8% in the intergenic region. Moreover, a total of 50 alleles and 42 polymorphic loci were detected among these A. gallica strains. The averaged polymorphism information content (PIC) was 0.4487, ranged from 0.2577 to 0.6786. Both principal coordinate analysis (PCoA) and population structure analyses based on the genotype data of SSRs divided the strains into two clusters. The cluster I included all the strains from high-altitude G. elata producing areas and some low-altitude areas, while the strains in Cluster II originated from low-altitude G. elata producing areas. These results indicated that substantial genome-specific variation in SSRs within the genic region of A. gallica and provide new insights for further studies on the evolution and breeding of A. gallica.