Microorganisms evolve via mechanisms spanning sexual/parasexual reproduction, mutators, aneuploidy, Hsp90, and even prions. Mechanisms that may seem detrimental can be repurposed to generate diversity. Here we show the human fungal pathogen Mucor circinelloides develops spontaneous resistance to the antifungal drug FK506 (tacrolimus) via two distinct mechanisms. One involves Mendelian mutations that confer stable drug resistance; the other occurs via an epigenetic RNA interference (RNAi)-mediated pathway resulting in unstable drug resistance. The peptidyl-prolyl isomerase FKBP12 interacts with FK506 forming a complex that inhibits the protein phosphatase calcineurin1. Calcineurin inhibition by FK506 blocks M. circinelloides transition to hyphae and enforces yeast growth2. Mutations in the fkbA gene encoding FKBP12 or the calcineurin cnbR or cnaA genes confer FK506 resistance (FK506R) and restore hyphal growth. In parallel, RNAi is spontaneously triggered to silence the FKBP12 fkbA gene, giving rise to drug-resistant epimutants. FK506R epimutants readily reverted to the drug-sensitive wild-type (WT) phenotype when grown without drug. The establishment of these epimutants is accompanied by generation of abundant fkbA small RNA (sRNA) and requires the RNAi pathway as well as other factors that constrain or reverse the epimutant state. Silencing involves generation of a double-stranded RNA (dsRNA) trigger intermediate from the fkbA mature mRNA to produce antisense fkbA RNA. This study uncovers a novel epigenetic RNAi-based epimutation mechanism controlling phenotypic plasticity, with possible implications for antimicrobial drug resistance and RNAi-regulatory mechanisms in fungi and other eukaryotes.

Endogenous short RNAs (esRNAs) play diverse roles in eukaryotes and usually are produced from double-stranded RNA (dsRNA) by Dicer. esRNAs are grouped into different classes based on biogenesis and function but not all classes are present in all three eukaryotic kingdoms. The esRNA register of fungi is poorly described compared to other eukaryotes and it is not clear what esRNA classes are present in this kingdom and whether they regulate the expression of protein coding genes. However, evidence that some dicer mutant fungi display altered phenotypes suggests that esRNAs play an important role in fungi. Here, we show that the basal fungus Mucor circinelloides produces new classes of esRNAs that map to exons and regulate the expression of many protein coding genes. The largest class of these exonic-siRNAs (ex-siRNAs) are generated by RNA-dependent RNA Polymerase 1 (RdRP1) and dicer-like 2 (DCL2) and target the mRNAs of protein coding genes from which they were produced. Our results expand the range of esRNAs in eukaryotes and reveal a new role for esRNAs in fungi.

Mucor circinelloides is a zygomycete fungus and an emerging opportunistic pathogen in immunocompromised patients, especially transplant recipients and in some cases otherwise healthy individuals. We have discovered a novel example of size dimorphism linked to virulence. M. circinelloides is a heterothallic fungus: (+) sex allele encodes SexP and (−) sex allele SexM, both of which are HMG domain protein sex determinants. M. circinelloides f. lusitanicus (Mcl) (−) mating type isolates produce larger asexual sporangiospores that are more virulent in the wax moth host compared to (+) isolates that produce smaller less virulent sporangiospores. The larger sporangiospores germinate inside and lyse macrophages, whereas the smaller sporangiospores do not. sexMΔ mutants are sterile and still produce larger virulent sporangiospores, suggesting that either the sex locus is not involved in virulence/spore size or the sexP allele plays an inhibitory role. Phylogenetic analysis supports that at least three extant subspecies populate the M. circinelloides complex in nature: Mcl, M. circinelloides f. griseocyanus, and M. circinelloides f. circinelloides (Mcc). Mcc was found to be more prevalent among clinical Mucor isolates, and more virulent than Mcl in a diabetic murine model in contrast to the wax moth host. The M. circinelloides sex locus encodes an HMG domain protein (SexP for plus and SexM for minus mating types) flanked by genes encoding triose phosphate transporter (TPT) and RNA helicase homologs. The borders of the sex locus between the three subspecies differ: the Mcg sex locus includes the promoters of both the TPT and the RNA helicase genes, whereas the Mcl and Mcc sex locus includes only the TPT gene promoter. Mating between subspecies was restricted compared to mating within subspecies. These findings demonstrate that spore size dimorphism is linked to virulence of M. circinelloides species and that plasticity of the sex locus and adaptations in pathogenicity have occurred during speciation of the M. circinelloides complex.