Analyses of the large-scale population structure of pathogens enable the identification of migration patterns, diversity reservoirs or longevity of populations, the understanding of current evolutionary trajectories and the anticipation of future ones. A detailed analysis of populations in centre of diversity should enable to infer the adaptive capacity of the pathogen and identify potential sources for new invasions. Puccinia striiformis f.sp. tritici (PST) is the causal agent of wheat yellow/stripe rust, and despite a worldwide distribution, this fungus remains a model species for invasion studies, due to its long-distance migration capacity and recurrent local emergence of new strains. Little is known about the ancestral relationship of the worldwide PST population with unknown center of origin. We used multilocus microsatellite genotyping to infer the worldwide population structure of PST and the origin of new invasions, analysing a set of isolates representative of sampling performed over six continents. Bayesian and multivariate clustering methods partitioned the isolates into six distinct genetic groups, corresponding to distinct geographic areas. The assignment analysis confirmed the Middle East-Red Sea Area as the most likely source of newly spreading, high-temperature-adapted strains; Europe as the source of South American, North American and Australian populations; and Mediterranean-Central Asian populations as the origin of South African populations. The existence of strong population subdivision at worldwide level shows that major genetic groups are not markedly affected by recent dispersal events. However, the sources for recent invasions and the migration routes identified emphasize the importance of human activities on the recent long-distance spread of the disease. The analyses of linkage disequilibrium and genotypic diversity indicated a strong regional heterogeneity in levels of recombination, with clear signatures of recombination in the Himalayan (Nepal and Pakistan) and near-Himalayan (China) regions and a predominant clonal population structure in other regions. To explain the variability in diversity and recombination of worldwide PST populations, we assessed their sex ability in terms of telial production, the sex-specific structures that are obligatory for PST sexual cycle, in a set of 56 isolates representative of these worldwide geographical origins. We confirmed that the variability in genotypic diversity/ recombination was linked with the sex ability, pinpointing the Himalayan region as the possible center of origin of PST, from where it then spread worldwide. The reduced sex ability in clonal populations certainly reflects a loss of sexual function, associated to migration in areas where sexual alternate host is lacking, or not necessary for the completion of epidemic cycle. Approximate Bayesian computation analyses confirmed an out of Himalaya spread of PST, with Pakistan and China being the most ancestral population. A detailed analysis of Pakistani population at regional level revealed the existence of a strong population subdivision, a high genotypic diversity and the existence of recombination signature at each location reflecting the role of sexual recombination in the temporal maintenance at local level. A time spaced sampling of PST in the valley of Tianshui (China) inspired the development of a new estimator, allowing to quantify the relative contribution of sexual reproduction and effective population size on the basis of clonal resampling within and between years. A sexual reproduction rate of 74% (95% confidence interval [CI]: 38-95%) and effective population size of 1735 (95% CI: 675-2800) was quantified in Chinese PST population. The description of the origin and migration routes of PST populations worldwide and at its centre of diversity contributes to our understanding of PST evolutionary potential, and is helpful to build disease management strategies.
Population biology and invasion history of puccinia striformis F.SP. tritici at worldwide and local scale