The genetic diversity of the soybean rust pathogen Phakopsora pachyrhizi has been driven by two major evolutionary lineages

Phakopsora pachyrhizi, an obligate biotrophic rust fungus, is the causal agent of Asian Soybean Rust (ASR) disease. Here, we utilized whole-genome data to explore the evolutionary patterns and population structure across 45 P. pachyrhizi isolates collected from 1972 to 2017 from diverse geographic regions worldwide. We also characterized in-silico mating-type (MAT) genes of P. pachyrhizi, in the predicted proteome of three isolates, to investigate the sexual compatibility system. Our molecular phylogenetic analysis in P. pachyrhizi inferred two distinct evolutionary lineages structured on a temporal scale, with lineage Pp1 grouping isolates obtained from 1972 to 1994, while more recently collected isolates formed a second lineage, Pp2. We found high levels of genetic diversity in lineage Pp1, whereas lineage Pp2 exhibited a strong clonal genetic structure, with a significant lower diversity. The widespread propagation of P. pachyrhizi clonal spores across soybean-growing regions likely explains the absence of a large-scale spatial genetic structure within each lineage. Two independent isolates (TW72-1 and AU79-1) showed moderate levels of genetic admixture, suggesting potential somatic hybridization between the two P. pachyrhizi lineages. We observed no clear congruence between virulence levels of P. pachyrhizi isolates and their phylogenetic patterns. Our findings support a probable tetrapolar mating system in P. pachyrhizi. Taken together, our study offers new insights into the evolutionary history of P. pachyrhizi and demonstrates that multiple MAT genes are highly expressed during the later stages of soybean infection, suggesting their potential role in the formation of urediniospores within the life cycle of P. pachyrhizi. AUTHOR SUMMARYThe Asian Soybean Rust (ASR) disease, caused by basidiomycetes fungus P. pachyrhizi, represents a critical threat to soybean crops worldwide. With the recent availability of high-quality genome assemblies for P. pachyrhizi, we are committed to exploring the genetic diversity of this destructive pathogen. This study analyzed whole-genome resequencing data from P. pachyrhizi isolates collected over several decades in various geographic regions. We identified recent diversification patterns in P. pachyrhizi, with two major lineages. The origin of the two lineages is likely due to temporal shifts in the genetic structure of P. pachyrhizi. Additionally, we investigated the genes responsible for sexual compatibility, known as mating-type (MAT) genes. Transcriptome data indicated that MAT genes are actively expressed during the later stages of the P. pachyrhizi-soybean infection. We conclude that the evolutionary dynamic of P. pachyrhizi is shaped by divergent lineages, which exhibit varying levels of virulence on differential soybean lines containing ASR resistance genes.