Exploring genetic diversity, population structure, and phylogeography in Paracoccidioides species using AFLP markers

Paracoccidioidomycosis (PCM) is a life-threatening systemic fungal infection acquired after inhalation of Paracoccidioides propagules from the environment. The main agents include members of the P. brasiliensis complex (phylogenetically-defined species S1, PS2, PS3, and PS4) and P. lutzii. DNA-sequencing of protein-coding loci (e.g., GP43, ARF, and TUB1) is the reference method for recognizing Paracoccidioides species due to a lack of robust phenotypic markers. Thus, developing new molecular markers that are informative and cost-effective is key to providing quality information to explore genetic diversity within Paracoccidioides. We report using new amplified fragment length polymorphism (AFLP) markers and mating-type analysis for genotyping Paracoccidioides species. The bioinformatic analysis generated 144 in silico AFLP profiles, highlighting two discriminatory primer pairs combinations (#1 EcoRI-AC/MseI-CT and #2 EcoRI-AT/MseI-CT). The combinations #1 and #2 were used in vitro to genotype 165 Paracoccidioides isolates recovered from across a vast area of South America. Considering the overall scored AFLP markers in vitro (67-87 fragments), the values of polymorphism information content (PIC = 0.3345-0.3456), marker index (MI = 0.0018), effective multiplex ratio (E = 44.6788-60.3818), resolving power (Rp = 22.3152-34.3152), discriminating power (D = 0.5183-0.5553), expected heterozygosity (H = 0.4247-0.4443), and mean heterozygosity (H _avp  = 0.00002-0.00004), demonstrated the utility of AFLP markers to speciate Paracoccidioides and to dissect both deep and fine-scale genetic structures. Analysis of molecular variance (AMOVA) revealed that the total genetic variance (65-66 %) was due to variability among P. brasiliensis complex and P. lutzii (PhiPT = 0.651-0.658, P < 0.0001), supporting a highly structured population. Heterothallism was the exclusive mating strategy, and the distributions of MAT1-1 or MAT1-2 idiomorphs were not significantly skewed (1:1 ratio) for P. brasiliensis s. str. (χ² = 1.025; P = 0.3113), P. venezuelensis (χ² = 0.692; P = 0.4054), and P. lutzii (χ² = 0.027; P = 0.8694), supporting random mating within each species. In contrast, skewed distributions were found for P. americana (χ² = 8.909; P = 0.0028) and P. restrepiensis (χ² = 4.571; P = 0.0325) with a preponderance of MAT1-1. Geographical distributions confirmed that P. americana, P. restrepiensis, and P. lutzii are more widespread than previously thought. P. brasiliensis s. str. is by far the most widely occurring lineage in Latin America countries, occurring in all regions of Brazil. Our new DNA fingerprint assay proved to be rapid, reproducible, and highly discriminatory, to give insights into the taxonomy, ecology, and epidemiology of Paracoccidioides species, guiding disease-control strategies to mitigate PCM.