Multilocus microsatellite genotyping array for investigation of genetic epidemiology of Pneumocystis jirovecii

Pneumocystis jirovecii is a symbiotic respiratory fungus that causes pneumonia (PcP) in immunosuppressed patients. Because P. jirovecii cannot be reliably cultured in vitro, it has proven difficult to study and gaps in our understanding of the organism persist. The release of a draft genome for the organism opens the door for the development of new genotyping approaches for studying its molecular epidemiology and global population structure. We identified and validated 8 putatively neutral microsatellite markers and 1 microsatellite marker linked to the dihydropteroate synthase gene (dhps), the enzymatic target of sulfa drugs used for PcP prevention and treatment. Using these tools, we analyzed P. jirovecii isolates from HIV-infected patients from three geographically distant populations: Uganda, the United States, and Spain. Among the 8 neutral markers, we observed high levels of allelic heterozygosity (average He, 0.586 to 0.842). Consistent with past reports, we observed limited global population structuring, with only the Ugandan isolates showing minor differentiation from the other two populations. In Ugandan isolates that harbored mutations in dhps, the microsatellite locus linked to dhps demonstrated a depressed He, consistent with positive directional selection for sulfa resistance mutations. Using a subset of these microsatellites, analyses of individual and paired samples from infections in San Francisco, CA, showed reliable typeability within a single infection and high discriminatory power between infections. These features suggest that this novel microsatellite typing approach will be an effective tool for molecular-epidemiological investigations into P. jirovecii population structure, transmission, and drug resistance.