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.
(χ
2
= 1.025;
P
= 0.3113),
P. venezuelensis
(χ
2
= 0.692;
P
= 0.4054), and
P. lutzii
(χ
2
= 0.027;
P
= 0.8694), supporting random mating within each species. In contrast, skewed distributions were found ...