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.

We carried out a population genomic survey of Saccharomyces cerevisiae diploid isolates and find that many budding yeast strains have high levels of genomic heterozygosity, much of which is likely due to outcrossing. We demonstrate that variation in heterozygosity among strains is correlated with a life-history tradeoff that involves how readily yeast switch from asexual to sexual reproduction under nutrient stress. This trade-off is reflected in a negative relationship between sporulation efficiency and pseudohyphal development and correlates with variation in the expression of RME1, a transcription factor with pleiotropic effects on meiosis and filamentous growth. Selection for alternate lifehistory strategies in natural versus human-associated environments likely contributes to differential maintenance of genomic heterozygosity through its effect on the frequency that yeast lineages experience sexual cycles and hence the opportunity for inbreeding. In addition to elevated levels of heterozygosity, many strains exhibit large genomic regions of loss-of-heterozygosity (LOH), suggesting that mitotic recombination has a significant impact on genetic variation in this species. This study provides new insights into the roles that both outcrossing and mitotic recombination play in shaping the genome architecture of Saccharomyces cerevisiae. This study also provides a unique case where stark differences in the genomic distribution of genetic variation among individuals of the same species can be largely explained by a lifehistory trade-off.T he frequency of sex and the nature of breeding systems have a profound effect on genome variation and evolution. For example, inbred populations have an increased frequency of homozygous genotypes (1), lower effective rates of recombination (2), and smaller effective population sizes relative to outcrossed populations with the same number of individuals (3). Likewise, clonal populations are expected to exhibit high levels of heterozygosity coupled with increased allelic diversity but decreased genotypic diversity relative to sexual populations (4).The budding yeast Saccharomyces cerevisiae is one of the best studied model organisms, but relatively little is known about the importance of sexual versus asexual reproduction and inbreeding versus outcrossing in shaping genome evolution in this species. One recent study estimated that outcrossing occurs approximately once every 50,000 generations in S. cerevisiae (5), but low rates of outcrossing do not preclude the possibility that outcrossing has an important impact on genetic variation. Studies of the closely related yeast Saccharomyces paradoxus suggest that sexual cycles are rare relative to asexual cycles and that when sex does occur it primarily involves inbreeding (6, 7). However, S. paradoxus exhibits distinctly different intra-and interpopulation patterns of variation than does S. cerevisiae (8), and hence these findings may not be generalizable across the Saccharomyces genus.Patterns of heterozygosity are an important indica...

RNA polymerase II (Pol II) termination is triggered by sequences present in the nascent transcript. Termination of pre-mRNA transcription is coupled to recognition of cis-acting sequences that direct cleavage and polyadenylation of the pre-mRNA. Termination of nonpolyadenylated [non-poly(A)] Pol II transcripts in Saccharomyces cerevisiae requires the RNA-binding proteins Nrd1 and Nab3. We have used a mutational strategy to characterize non-poly(A) termination elements downstream of the SNR13 and SNR47 snoRNA genes. This approach detected two common RNA sequence motifs, GUA[AG] and UCUU. The first motif corresponds to the known Nrd1-binding site, which we have verified here by gel mobility shift assays. We also show that Nab3 protein binds specifically to RNA containing the UCUU motif. Taken together, our data suggest that Nrd1 and Nab3 binding sites play a significant role in defining non-poly(A) terminators. As is the case with poly(A) terminators, there is no strong consensus for non-poly(A) terminators, and the arrangement of Nrd1p and Nab3p binding sites varies considerably. In addition, the organization of these sequences is not strongly conserved among even closely related yeasts. This indicates a large degree of genetic variability. Despite this variability, we were able to use a computational model to show that the binding sites for Nrd1 and Nab3 can identify genes for which transcription termination is mediated by these proteins.RNA polymerases synthesize discrete transcripts by initiating and terminating transcription in response to specific sequence elements. To initiate transcription, DNA sequences direct the binding of transcription initiation factors and polymerase to appropriate sites upstream of genes. Termination signals, providing for the creation of functional 3Ј ends, may reside either in the DNA template or, alternatively, are found in the nascent RNA transcript.Eucaryotic RNA polymerases have evolved distinct mechanisms for termination. RNA polymerase III (Pol III) requires no protein factors but terminates efficiently after transcribing four to six consecutive U residues, presumably due to instability of the RNA-DNA hybrid in the enzyme active site (1, 7). RNA Pol I terminates in response to a protein factor, Reb1, which blocks further elongation by binding to a DNA sequence downstream of the termination site (28). The Reb1 site is situated in such a way that the paused polymerase contains an inherently unstable U-rich RNA-DNA hybrid in the active site.The RNA Pol II termination mechanism is more complex than those employed by the other eucaryotic RNA polymerases, requiring a large multiprotein complex that recognizes the poly(A) signal in the nascent transcript (4,17,20,21,27,43). Deletion or mutation of the poly(A) signal results in a failure to terminate messages at the appropriate site (9, 19). This observation suggests that processing of the nascent message is coupled to termination, but the mechanism of this coupling remains uncertain. In Saccharomyces cerevisiae, mutations in severa...

Nutrient stresses trigger a variety of developmental switches in the budding yeast Saccharomyces cerevisiae. One of the least understood of such responses is the development of complex colony morphology, characterized by intricate, organized, and strain-specific patterns of colony growth and architecture. The genetic bases of this phenotype and the key environmental signals involved in its induction have heretofore remained poorly understood. By surveying multiple strain backgrounds and a large number of growth conditions, we show that limitation for fermentable carbon sources coupled with a rich nitrogen source is the primary trigger for the colony morphology response in budding yeast. Using knockout mutants and transposon-mediated mutagenesis, we demonstrate that two key signaling networks regulating this response are the filamentous growth MAP kinase cascade and the Ras-cAMP-PKA pathway. We further show synergistic epistasis between Rim15, a kinase involved in integration of nutrient signals, and other genes in these pathways. Ploidy, mating-type, and genotype-by-environment interactions also appear to play a role in the controlling colony morphology. Our study highlights the high degree of network reuse in this model eukaryote; yeast use the same core signaling pathways in multiple contexts to integrate information about environmental and physiological states and generate diverse developmental outputs.

SUMMARY Gene inactivation by transposon insertion or allelic exchange is a powerful approach to probe gene function. Unfortunately, many microbes, including Chlamydia, are not amenable to routine molecular genetic manipulations. Here we describe an arrayed library of chemically-induced mutants of the genetically-intransigent pathogen Chlamydia trachomatis, in which all mutations have been identified by whole genome sequencing, providing a platform for reverse genetic applications. An analysis of possible loss-of-function mutations in the collection uncovered plasticity in the central metabolic properties of this obligate intracellular pathogen. We also describe the use of the library in a forward genetic screen that identified InaC as a bacterial factor that binds host ARF and 14-3-3 proteins to modulate F-actin assembly and Golgi redistribution around the pathogenic vacuole. This work provides a robust platform for reverse and forward genetic approaches in Chlamydia and should serve as a valuable resource to the community.

Species within the human pathogenic Cryptococcus species complex are major threats to public health, causing approximately 1 million annual infections globally. Cryptococcus amylolentus is the most closely known related species of the pathogenic Cryptococcus species complex, and it is non-pathogenic. Additionally, while pathogenic Cryptococcus species have bipolar mating systems with a single large mating type (MAT) locus that represents a derived state in Basidiomycetes, C. amylolentus has a tetrapolar mating system with 2 MAT loci (P/R and HD) located on different chromosomes. Thus, studying C. amylolentus will shed light on the transition from tetrapolar to bipolar mating systems in the pathogenic Cryptococcus species, as well as its possible link with the origin and evolution of pathogenesis. In this study, we sequenced, assembled, and annotated the genomes of 2 C. amylolentus isolates, CBS6039 and CBS6273, which are sexual and interfertile. Genome comparison between the 2 C. amylolentus isolates identified the boundaries and the complete gene contents of the P/R and HD MAT loci. Bioinformatic and chromatin immunoprecipitation sequencing (ChIP-seq) analyses revealed that, similar to those of the pathogenic Cryptococcus species, C. amylolentus has regional centromeres (CENs) that are enriched with species-specific transposable and repetitive DNA elements. Additionally, we found that while neither the P/R nor the HD locus is physically closely linked to its centromere in C. amylolentus, and the regions between the MAT loci and their respective centromeres show overall synteny between the 2 genomes, both MAT loci exhibit genetic linkage to their respective centromere during meiosis, suggesting the presence of recombinational suppressors and/or epistatic gene interactions in the MAT-CEN intervening regions. Furthermore, genomic comparisons between C. amylolentus and related pathogenic Cryptococcus species provide evidence that multiple chromosomal rearrangements mediated by intercentromeric recombination have occurred during descent of the 2 lineages from their common ancestor. Taken together, our findings support a model in which the evolution of the bipolar mating system was initiated by an ectopic recombination event mediated by similar repetitive centromeric DNA elements shared between chromosomes. This translocation brought the P/R and HD loci onto the same chromosome, and further chromosomal rearrangements then resulted in the 2 MAT loci becoming physically linked and eventually fusing to form the single contiguous MAT locus that is now extant in the pathogenic Cryptococcus species.

MutY is an adenine-DNA glycosylase with specificity for mismatches involving 8-oxoguanine (oG.A) or guanine (G.A). In addition to a 25 kDa catalytic domain common to all members of its DNA glycosylase superfamily, MutY has a 14 kDa C-terminal domain. Sequence analyses suggest that this C-terminal domain is distantly related to MutT, a pyrophosphohydrolase specific for 2'-deoxy-8-oxoguanosine triphosphate (doGTP). Here we present biochemical evidence that the MutT-like domain of MutY is the principal determinant of oG specificity. First, MutY dissociates approximately 1500-fold more slowly from oG-containing product DNA than from G-containing product, but a truncated protein lacking the C-terminal domain dissociates as rapidly from oG-DNA as the full-length protein dissociates from G-DNA. Second, MutY removes adenine from oG.A mismatches almost 30-fold faster than from G.A mismatches in a pre-steady-state assay, but deletion of the C-terminal domain reduces this specificity for oG.A to less than 4-fold. The kinetic data are consistent with a model in which binding of oG to the C-terminal domain of MutY accelerates the pre-steady-state glycosylase reaction by facilitating adenine base flipping. The observation that oG specificity derives almost exclusively from the C-terminal domain of MutY adds credence to the sequence analyses and suggests that specificity for oG.A mismatches was acquired by fusion of a MutT-like protein onto the core catalytic domain of an adenine-DNA glycosylase.

Calcineurin is a highly conserved Ca2+/calmodulin-dependent serine/threonine-specific protein phosphatase that orchestrates cellular Ca2+ signaling responses. In Cryptococcus neoformans, calcineurin is activated by multiple stresses including high temperature, and is essential for stress adaptation and virulence. The transcription factor Crz1 is a major calcineurin effector in Saccharomyces cerevisiae and other fungi. Calcineurin dephosphorylates Crz1, thereby enabling Crz1 nuclear translocation and transcription of target genes. Here we show that loss of Crz1 confers phenotypes intermediate between wild-type and calcineurin mutants, and demonstrate that deletion of the calcineurin docking domain results in the inability of Crz1 to translocate into the nucleus under thermal stress. RNA-sequencing revealed 102 genes that are regulated in a calcineurin-Crz1-dependent manner at 37°C. The majority of genes were down-regulated in cna1Δ and crz1Δ mutants, indicating these genes are normally activated by the calcineurin-Crz1 pathway at high temperature. About 58% of calcineurin-Crz1 target genes have unknown functions, while genes with known or predicted functions are involved in cell wall remodeling, calcium transport, and pheromone production. We identified three calcineurin-dependent response element motifs within the promoter regions of calcineurin-Crz1 target genes, and show that Crz1 binding to target gene promoters is increased upon thermal stress in a calcineurin-dependent fashion. Additionally, we found a large set of genes independently regulated by calcineurin, and Crz1 regulates 59 genes independently of calcineurin. Given the intermediate crz1Δ mutant phenotype, and our recent evidence for a calcineurin regulatory network impacting mRNA in P-bodies and stress granules independently of Crz1, calcineurin likely acts on factors beyond Crz1 that govern mRNA expression/stability to operate a branched transcriptional/post-transcriptional stress response network necessary for fungal virulence. Taken together, our findings reveal the core calcineurin-Crz1 stress response cascade is maintained from ascomycetes to a pathogenic basidiomycete fungus, but its output in C. neoformans appears to be adapted to promote fungal virulence.