The human fungal pathogen Mucor circinelloides develops spontaneous resistance to an antifungal drug both through mutation and through a newly identified epigenetic RNA-mediated pathway; RNA interference is spontaneously triggered to silence the fkbA gene, giving rise to drug-resistant epimutants that revert to being drug-sensitive once again when grown in the absence of drug. RNA interference (RNAi) is a mechanism conserved across eukaryotes that controls multiple cellular functions. This study reports that the opportunistic human pathogen Mucor circinelloides can develop spontaneous resistance to the antifungal drug FK506 (tacrolimus) via two distinct mechanisms. One is through conventional Mendelian mutation, whereas the other, surprisingly, is via a newly identified epigenetic RNAi-mediated pathway. Joseph Heitman and colleagues show that RNAi is spontaneously triggered to silence a gene, fkbA, that encodes the peptidylprolyl isomerase FKBP12. This enzyme interacts with the drug to form a complex that inhibits calcineurin, blocking the transition to hyphae. The resulting drug-resistant 'epimutants' revert to drug sensitivity when grown in the absence of drug. Microorganisms evolve via a range of mechanisms that may include or involve sexual/parasexual reproduction, mutators, aneuploidy, Hsp90 and even prions. Mechanisms that may seem detrimental can be repurposed to generate diversity. Here we show that 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 peptidylprolyl 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 and restore hyphal growth. In parallel, RNAi is spontaneously triggered to silence the fkbA gene, giving rise to drug-resistant epimutants. FK506-resistant epimutants readily reverted to the drug-sensitive wild-type phenotype when grown without exposure to the drug. The establishment of these epimutants is accompanied by generation of abundant fkbA small RNAs and requires the RNAi pathway as well as other factors that constrain or reverse the epimutant state. Silencing involves the generation of a double-stranded RNA trigger intermediate using the fkbA mature mRNA as a template 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.

Transformation of Mucor circinelloides with self‐replicative plasmids containing a wild‐type copy of the carotenogenic gene carB causes silencing of the carB function in 3% of transformants. Genomic analyses revealed a relationship between silenced phenotype and number of copies of plasmids. This phenotype results from a reduction of the steady‐state levels of carB mRNA, a reduction that is not due to differences in the level of transcription, indicating that silencing is post‐transcriptional. Small sense and antisense RNAs have been found to be associated with gene silencing in M.circinelloides . Two size classes of small antisense RNAs, differentially accumulated during the vegetative growth of silenced transformants, have been detected: a long 25‐nucleotide RNA and a short 21‐nucleotide RNA. Secondary sense and antisense RNAs corresponding to sequences of the endogenous gene downstream of the initial triggering molecule have also been detected, revealing the existence of spreading of RNA targeting in fungi. These findings, together with the self‐replicative nature of the triggering molecules, make M.circinelloides a suitable organism for investigating some unresolved questions in RNA silencing.

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.

MRNA profiling is routinely used to identify microRNA targets, however, this high-throughput technology is not suitable for identifying targets regulated only at protein level. Here, we have developed and validated a novel methodology based on computational analysis of promoter sequences combined with mRNA microarray experiments to reveal transcription factors that are direct microRNA targets at the protein level. Using this approach we identified Smad3, a key transcription factor in the TGFβ signaling pathway, as a direct miR-140 target. We showed that miR-140 suppressed the TGFβ pathway through repression of Smad3 and that TGFβ suppressed the accumulation of miR-140 forming a double negative feedback loop. Our findings establish a valid strategy for the discovery of microRNA targets regulated only at protein level, and we propose that additional targets could be identified by re-analysis of existing microarray datasets.

MicroRNAs (miRNAs) are short noncoding RNA molecules regulating the expression of mRNAs. Target identification of miRNAs is computationally difficult due to the relatively low homology between miRNAs and their targets. We present here an experimental approach to target identification where the cartilage-specific miR-140 was overexpressed and silenced in cells it is normally expressed in separate experiments. Expression of mRNAs was profiled in both experiments and the intersection of mRNAs repressed by miR-140 overexpression and derepressed by silencing of miR-140 was identified. The intersection contained only 49 genes, although both treatments affected the accumulation of hundreds of mRNAs. These 49 genes showed a very strong enrichment for the miR-140 seed sequence implying that the approach is efficient and specific. Twenty-one of these 49 genes were predicted to be direct targets based on the presence of the seed sequence. Interestingly, none of these were predicted by the published target prediction methods we used. One of the potential target mRNAs, Cxcl12, was experimentally validated by Northern blot analysis and a luciferase reporter assay.