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.

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...

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...

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.

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.

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.

Sequence motifs that are potentially recognized by DNA-binding proteins occur far more often in genomic DNA than do observed in vivo protein–DNA interactions. To determine how chromatin influences the utilization of particular DNA-binding sites, we compared the in vivo genome-wide binding location of the yeast transcription factor Leu3 to the binding location observed on the same genomic DNA in the absence of any protein cofactors. We found that the DNA-sequence motif recognized by Leu3 in vitro and in vivo was functionally indistinguishable, but Leu3 bound different genomic locations under the two conditions. Accounting for nucleosome occupancy in addition to DNA-sequence motifs significantly improved the prediction of protein–DNA interactions in vivo, but not the prediction of sites bound by purified Leu3 in vitro. Use of histone modification data does not further improve binding predictions, presumably because their effect is already manifest in the global histone distribution. Measurements of nucleosome occupancy in strains that differ in Leu3 genotype show that low nucleosome occupancy at loci bound by Leu3 is not a consequence of Leu3 binding. These results permit quantitation of the epigenetic influence that chromatin exerts on DNA binding-site selection, and provide evidence for an instructive, functionally important role for nucleosome occupancy in determining patterns of regulatory factor targeting genome-wide.