High-Resolution Mapping of Two Types of Spontaneous Mitotic Gene Conversion Events in Saccharomyces cerevisiae

Gene conversions and crossovers are related products of the repair of double-stranded DNA breaks by homologous recombination. Most previous studies of mitotic gene conversion events have been restricted to measuring conversion tracts that are <5 kb. Using a genetic assay in which the lengths of very long gene conversion tracts can be measured, we detected two types of conversions: those with a median size of ∼6 kb and those with a median size of >50 kb. The unusually long tracts are initiated at a naturally occurring recombination hotspot formed by two inverted Ty elements. We suggest that these long gene conversion events may be generated by a mechanism (break-induced replication or repair of a double-stranded DNA gap) different from the short conversion tracts that likely reflect heteroduplex formation followed by DNA mismatch repair. Both the short and long mitotic conversion tracts are considerably longer than those observed in meiosis. Since mitotic crossovers in a diploid can result in a heterozygous recessive deleterious mutation becoming homozygous, it has been suggested that the repair of DNA breaks by mitotic recombination involves gene conversion events that are unassociated with crossing over. In contrast to this prediction, we found that ∼40% of the conversion tracts are associated with crossovers. Spontaneous mitotic crossover events in yeast are frequent enough to be an important factor in genome evolution.     

Multiple Pathways Suppress Non-Allelic Homologous Recombination during Meiosis in Saccharomyces cerevisiae

Recombination during meiosis in the form of crossover events promotes the segregation of homologous chromosomes by providing the only physical linkage between these chromosomes. Recombination occurs not only between allelic sites but also between non-allelic (ectopic) sites. Ectopic recombination is often suppressed to prevent non-productive linkages. In this study, we examined the effects of various mutations in genes involved in meiotic recombination on ectopic recombination during meiosis. RAD24, a DNA damage checkpoint clamp-loader gene, suppressed ectopic recombination in wild type in the same pathway as RAD51. In the absence of RAD24, a meiosis-specific recA homolog, DMC1, suppressed the recombination. Homology search and strand exchange in ectopic recombination occurred when either the RAD51 or the DMC1 recA homolog was absent, but was promoted by RAD52. Unexpectedly, the zip1 mutant, which is defective in chromosome synapsis, showed a decrease, rather than an increase, in ectopic recombination. Our results provide evidence for two types of ectopic recombination: one that occurs in wild-type cells and a second that occurs predominantly when the checkpoint pathway is inactivated.     

Meiotic Recombination Initiation in and around Retrotransposable Elements in Saccharomyces cerevisiae

Meiotic recombination is initiated by large numbers of developmentally programmed DNA double-strand breaks (DSBs), ranging from dozens to hundreds per cell depending on the organism. DSBs formed in single-copy sequences provoke recombination between allelic positions on homologous chromosomes, but DSBs can also form in and near repetitive elements such as retrotransposons. When they do, they create a risk for deleterious genome rearrangements in the germ line via recombination between non-allelic repeats. A prior study in budding yeast demonstrated that insertion of a Ty retrotransposon into a DSB hotspot can suppress meiotic break formation, but properties of Ty elements in their most common physiological contexts have not been addressed. Here we compile a comprehensive, high resolution map of all Ty elements in the rapidly and efficiently sporulating S. cerevisiae strain SK1 and examine DSB formation in and near these endogenous retrotransposable elements. SK1 has 30 Tys, all but one distinct from the 50 Tys in S288C, the source strain for the yeast reference genome. From whole-genome DSB maps and direct molecular assays, we find that DSB levels and chromatin structure within and near Tys vary widely between different elements and that local DSB suppression is not a universal feature of Ty presence. Surprisingly, deletion of two Ty elements weakened adjacent DSB hotspots, revealing that at least some Ty insertions promote rather than suppress nearby DSB formation. Given high strain-to-strain variability in Ty location and the high aggregate burden of Ty-proximal DSBs, we propose that meiotic recombination is an important component of host-Ty interactions and that Tys play critical roles in genome instability and evolution in both inbred and outcrossed sexual cycles.     

Rad51 Inhibits Translocation Formation by Non-Conservative Homologous Recombination in Saccharomyces cerevisiae

Chromosomal translocations are a primary biological response to ionizing radiation (IR) exposure, and are likely to result from the inappropriate repair of the DNA double-strand breaks (DSBs) that are created. An abundance of repetitive sequences in eukaryotic genomes provides ample opportunity for such breaks to be repaired by homologous recombination (HR) between non-allelic repeats. Interestingly, in the budding yeast, Saccharomyces cerevisiae the central strand exchange protein, Rad51 that is required for DSB repair by gene conversion between unlinked repeats that conserves genomic structure also suppresses translocation formation by several HR mechanisms. In particular, Rad51 suppresses translocation formation by single-strand annealing (SSA), perhaps the most efficient mechanism for translocation formation by HR in both yeast and mammalian cells. Further, the enhanced translocation formation that emerges in the absence of Rad51 displays a distinct pattern of genetic control, suggesting that this occurs by a separate mechanism. Since hypomorphic mutations in RAD51 in mammalian cells also reduce DSB repair by conservative gene conversion and stimulate non-conservative repair by SSA, this mechanism may also operate in humans and, perhaps contribute to the genome instability that propels the development of cancer.     

Synthesis of Morphinan Alkaloids in Saccharomyces cerevisiae

Morphinan alkaloids are the most powerful narcotic analgesics currently used to treat moderate to severe and chronic pain. The feasibility of morphinan synthesis in recombinant Saccharomyces cerevisiae starting from the precursor (R,S)-norlaudanosoline was investigated. Chiral analysis of the reticuline produced by the expression of opium poppy methyltransferases showed strict enantioselectivity for (S)-reticuline starting from (R,S)-norlaudanosoline. In addition, the P. somniferum enzymes salutaridine synthase (PsSAS), salutaridine reductase (PsSAR) and salutaridinol acetyltransferase (PsSAT) were functionally co-expressed in S. cerevisiae and optimization of the pH conditions allowed for productive spontaneous rearrangement of salutaridinol-7-O-acetate and synthesis of thebaine from (R)-reticuline. Finally, we reconstituted a 7-gene pathway for the production of codeine and morphine from (R)-reticuline. Yeast cell feeding assays using (R)-reticuline, salutaridine or codeine as substrates showed that all enzymes were functionally co-expressed in yeast and that activity of salutaridine reductase and codeine-O-demethylase likely limit flux to morphine synthesis. The results of this study describe a significant advance for the synthesis of morphinans in S. cerevisiae and pave the way for their complete synthesis in recombinant microbes.     

Comparative Genome-Wide Screening Identifies a Conserved Doxorubicin Repair Network That Is Diploid Specific in Saccharomyces cerevisiae

The chemotherapeutic doxorubicin (DOX) induces DNA double-strand break (DSB) damage. In order to identify conserved genes that mediate DOX resistance, we screened the Saccharomyces cerevisiae diploid deletion collection and identified 376 deletion strains in which exposure to DOX was lethal or severely reduced growth fitness. This diploid screen identified 5-fold more DOX resistance genes than a comparable screen using the isogenic haploid derivative. Since DSB damage is repaired primarily by homologous recombination in yeast, and haploid cells lack an available DNA homolog in G1 and early S phase, this suggests that our diploid screen may have detected the loss of repair functions in G1 or early S phase prior to complete DNA replication. To test this, we compared the relative DOX sensitivity of 30 diploid deletion mutants identified under our screening conditions to their isogenic haploid counterpart, most of which (n = 26) were not detected in the haploid screen. For six mutants (bem1Δ, ctf4Δ, ctk1Δ, hfi1Δ,nup133Δ, tho2Δ) DOX-induced lethality was absent or greatly reduced in the haploid as compared to the isogenic diploid derivative. Moreover, unlike WT, all six diploid mutants displayed severe G1/S phase cell cycle progression defects when exposed to DOX and some were significantly enhanced (ctk1Δ and hfi1Δ) or deficient (tho2Δ) for recombination. Using these and other “THO2-like” hypo-recombinogenic, diploid-specific DOX sensitive mutants (mft1Δ, thp1Δ, thp2Δ) we utilized known genetic/proteomic interactions to construct an interactive functional genomic network which predicted additional DOX resistance genes not detected in the primary screen. Most (76%) of the DOX resistance genes detected in this diploid yeast screen are evolutionarily conserved suggesting the human orthologs are candidates for mediating DOX resistance by impacting on checkpoint and recombination functions in G1 and/or early S phases.     

Towards the Systematic Mapping and Engineering of the Protein Prenylation Machinery in Saccharomyces cerevisiae

Protein prenylation is a widespread and highly conserved eukaryotic post-translational modification that endows proteins with the ability to reversibly attach to intracellular membranes. The dynamic interaction of prenylated proteins with intracellular membranes is essential for their signalling functions and is frequently deregulated in disease processes such as cancer. As a result, protein prenylation has been pharmacologically targeted by numerous drug discovery programs, albeit with limited success. To a large extent, this can be attributed to an insufficient understanding of the interplay of different protein prenyltransferases and the combinatorial diversity of the prenylatable sequence space. Here, we report a high-throughput, growth-based genetic selection assay in Saccharomyces cerevisiae based on the Ras Recruitment System which, for the first time, has allowed us to create a comprehensive map of prenylatable protein sequences in S. cerevisiae. We demonstrate that potential prenylatable space is sparsely (6.2%) occupied leaving room for creation of synthetic orthogonal prenylatable sequences. To experimentally demonstrate that, we used the developed platform to engineer mutant farnesyltransferases that efficiently prenylate substrate motives that are not recognised by endogenous protein prenyltransferases. These uncoupled mutants can now be used as starting points for the systematic engineering of the eukaryotic protein prenylation machinery.     

Suppression of Allelic Recombination and Aneuploidy by Cohesin Is Independent of Chk1 in Saccharomyces cerevisiae

Sister chromatid cohesion (SCC), which is established during DNA replication, ensures genome stability. Establishment of SCC is inhibited in G2. However, this inhibition is relived and SCC is established as a response to DNA damage, a process known as Damage Induced Cohesion (DIC). In yeast, Chk1, which is a kinase that functions in DNA damage signal transduction, is considered an activator of SCC through DIC. Nonetheless, here we show that, unlike SCC mutations, loss of CHK1 did not increase spontaneous or damage-induced allelic recombination or aneuploidy. We suggest that Chk1 has a redundant role in the control of DIC or that DIC is redundant for maintaining genome stability.     

Genealogy-Based Methods for Inference of Historical Recombination and Gene Flow and Their Application in Saccharomyces cerevisiae

Genetic exchange between isolated populations, or introgression between species, serves as a key source of novel genetic material on which natural selection can act. While detecting historical gene flow from DNA sequence data is of much interest, many existing methods can be limited by requirements for deep population genomic sampling. In this paper, we develop a scalable genealogy-based method to detect candidate signatures of gene flow into a given population when the source of the alleles is unknown. Our method does not require sequenced samples from the source population, provided that the alleles have not reached fixation in the sampled recipient population. The method utilizes recent advances in algorithms for the efficient reconstruction of ancestral recombination graphs, which encode genealogical histories of DNA sequence data at each site, and is capable of detecting the signatures of gene flow whose footprints are of length up to single genes. Further, we employ a theoretical framework based on coalescent theory to test for statistical significance of certain recombination patterns consistent with gene flow from divergent sources. Implementing these methods for application to whole-genome sequences of environmental yeast isolates, we illustrate the power of our approach to highlight loci with unusual recombination histories. By developing innovative theory and methods to analyze signatures of gene flow from population sequence data, our work establishes a foundation for the continued study of introgression and its evolutionary relevance.     

Kinetics of Phosphomevalonate Kinase from Saccharomyces cerevisiae

The mevalonate-based isoprenoid biosynthetic pathway is responsible for producing cholesterol in humans and is used commercially to produce drugs, chemicals, and fuels. Heterologous expression of this pathway in Escherichia coli has enabled high-level production of the antimalarial drug artemisinin and the proposed biofuel bisabolane. Understanding the kinetics of the enzymes in the biosynthetic pathway is critical to optimize the pathway for high flux. We have characterized the kinetic parameters of phosphomevalonate kinase (PMK, EC 2.7.4.2) from Saccharomyces cerevisiae, a previously unstudied enzyme. An E. coli codon-optimized version of the S. cerevisiae gene was cloned into pET-52b+, then the C-terminal 6X His-tagged protein was expressed in E. coli BL21(DE3) and purified on a Ni²⁺ column. The K_M of the ATP binding site was determined to be 98.3 µM at 30°C, the optimal growth temperature for S. cerevisiae, and 74.3 µM at 37°C, the optimal growth temperature for E. coli. The K_M of the mevalonate-5-phosphate binding site was determined to be 885 µM at 30°C and 880 µM at 37°C. The V_max was determined to be 4.51 µmol/min/mg enzyme at 30°C and 5.33 µmol/min/mg enzyme at 37°C. PMK is Mg²⁺ dependent, with maximal activity achieved at concentrations of 10 mM or greater. Maximum activity was observed at pH = 7.2. PMK was not found to be substrate inhibited, nor feedback inhibited by FPP at concentrations up to 10 µM FPP.     

Evolutionary Genomics of Transposable Elements in Saccharomyces cerevisiae

Saccharomyces cerevisiae is one of the premier model systems for studying the genomics and evolution of transposable elements. The availability of the S. cerevisiae genome led to unprecedented insights into its five known transposable element families (the LTR retrotransposons Ty1-Ty5) in the years shortly after its completion. However, subsequent advances in bioinformatics tools for analysing transposable elements and the recent availability of genome sequences for multiple strains and species of yeast motivates new investigations into Ty evolution in S. cerevisiae. Here we provide a comprehensive phylogenetic and population genetic analysis of all Ty families in S. cerevisiae based on a systematic re-annotation of Ty elements in the S288c reference genome. We show that previous annotation efforts have underestimated the total copy number of Ty elements for all known families. In addition, we identify a new family of Ty3-like elements related to the S. paradoxus Ty3p which is composed entirely of degenerate solo LTRs. Phylogenetic analyses of LTR sequences identified three families with short-branch, recently active clades nested among long branch, inactive insertions (Ty1, Ty3, Ty4), one family with essentially all recently active elements (Ty2) and two families with only inactive elements (Ty3p and Ty5). Population genomic data from 38 additional strains of S. cerevisiae show that the majority of Ty insertions in the S288c reference genome are fixed in the species, with insertions in active clades being predominantly polymorphic and insertions in inactive clades being predominantly fixed. Finally, we use comparative genomic data to provide evidence that the Ty2 and Ty3p families have arisen in the S. cerevisiae genome by horizontal transfer. Our results demonstrate that the genome of a single individual contains important information about the state of TE population dynamics within a species and suggest that horizontal transfer may play an important role in shaping the genomic diversity of transposable elements in unicellular eukaryotes.     

High-Resolution Mapping of Spontaneous Mitotic Recombination Hotspots on the 1.1 Mb Arm of Yeast Chromosome IV

Although homologous recombination is an important pathway for the repair of double-stranded DNA breaks in mitotically dividing eukaryotic cells, these events can also have negative consequences, such as loss of heterozygosity (LOH) of deleterious mutations. We mapped about 140 spontaneous reciprocal crossovers on the right arm of the yeast chromosome IV using single-nucleotide-polymorphism (SNP) microarrays. Our mapping and subsequent experiments demonstrate that inverted repeats of Ty retrotransposable elements are mitotic recombination hotspots. We found that the mitotic recombination maps on the two homologs were substantially different and were unrelated to meiotic recombination maps. Additionally, about 70% of the DNA lesions that result in LOH are likely generated during G1 of the cell cycle and repaired during S or G2. We also show that different genetic elements are associated with reciprocal crossover conversion tracts depending on the cell cycle timing of the initiating DSB.     

Distance-Independence of Mitotic Intrachromosomal Recombination in Saccharomyces Cerevisiae

Many genetic studies have shown that the frequency of homologous recombination depends largely on the distance in which recombination can occur. We have studied the effect of varying the length of duplicated sequences on the frequency of mitotic intrachromosomal recombination in Saccharomyces cerevisiae. We find that the frequency of recombination resulting in the loss of one of the repeats and the intervening sequences reaches a plateau when the repeats are short. In addition, the frequency of recombination to correct a point mutation contained in one of these repeats is not proportional to the size of the duplication but rather depends dramatically on the location of the mutation within the repeated sequences. However, the frequency of mitotic interchromosomal reciprocal recombination is dependent on the distance separating the markers. The difference in the response of intrachromosomal and interchromosomal mitotic recombination to increasing lengths of homology may indicate there are different rate-limiting steps for recombination in these two cases. These findings have important implications for the maintenance and evolution of duplicated sequences.     

High-Resolution Mapping of Crossover and Non-crossover Recombination Events by Whole-Genome Re-sequencing of an Avian Pedigree

Recombination is an engine of genetic diversity and therefore constitutes a key process in evolutionary biology and genetics. While the outcome of crossover recombination can readily be detected as shuffled alleles by following the inheritance of markers in pedigreed families, the more precise location of both crossover and non-crossover recombination events has been difficult to pinpoint. As a consequence, we lack a detailed portrait of the recombination landscape for most organisms and knowledge on how this landscape impacts on sequence evolution at a local scale. To localize recombination events with high resolution in an avian system, we performed whole-genome re-sequencing at high coverage of a complete three-generation collared flycatcher pedigree. We identified 325 crossovers at a median resolution of 1.4 kb, with 86% of the events localized to <10 kb intervals. Observed crossover rates were in excellent agreement with data from linkage mapping, were 52% higher in male (3.56 cM/Mb) than in female meiosis (2.28 cM/Mb), and increased towards chromosome ends in male but not female meiosis. Crossover events were non-randomly distributed in the genome with several distinct hot-spots and a concentration to genic regions, with the highest density in promoters and CpG islands. We further identified 267 non-crossovers, whose location was significantly associated with crossover locations. We detected a significant transmission bias (0.18) in favour of ‘strong’ (G, C) over ‘weak’ (A, T) alleles at non-crossover events, providing direct evidence for the process of GC-biased gene conversion in an avian system. The approach taken in this study should be applicable to any species and would thereby help to provide a more comprehensive portray of the recombination landscape across organism groups.     

Branched-Chain Aminotransferases Control TORC1 Signaling in Saccharomyces cerevisiae

The conserved target of rapamycin complex 1 (TORC1) integrates nutrient signals to orchestrate cell growth and proliferation. Leucine availability is conveyed to control TORC1 activity via the leu-tRNA synthetase/EGOC-GTPase module in yeast and mammals, but the mechanisms sensing leucine remain only partially understood. We show here that both leucine and its α-ketoacid metabolite, α-ketoisocaproate, effectively activate the yeast TORC1 kinase via both EGOC GTPase-dependent and -independent mechanisms. Leucine and α-ketoisocaproate are interconverted by ubiquitous branched-chain aminotransferases (BCAT), which in yeast are represented by the mitochondrial and cytosolic enzymes Bat1 and Bat2, respectively. BCAT yeast mutants exhibit severely compromised TORC1 activity, which is partially restored by expression of Bat1 active site mutants, implicating both catalytic and structural roles of BCATs in TORC1 control. We find that Bat1 interacts with branched-chain amino acid metabolic enzymes and, in a leucine-dependent fashion, with the tricarboxylic acid (TCA)-cycle enzyme aconitase. BCAT mutation perturbed TCA-cycle intermediate levels, consistent with a TCA-cycle block, and resulted in low ATP levels, activation of AMPK, and TORC1 inhibition. We propose the biosynthetic capacity of BCAT and its role in forming multicomplex metabolons connecting branched-chain amino acids and TCA-cycle metabolism governs TCA-cycle flux to activate TORC1 signaling. Because mammalian mitochondrial BCAT is known to form a supramolecular branched-chain α-keto acid dehydrogenase enzyme complex that links leucine metabolism to the TCA-cycle, these findings establish a precedent for understanding TORC1 signaling in mammals.     

Visualization of Atg3 during Autophagosome Formation in Saccharomyces cerevisiae

Macroautophagy (autophagy) is a highly conserved cellular recycling process involved in degradation of eukaryotic cellular components. During autophagy, macromolecules and organelles are sequestered into the double-membrane autophagosome and degraded in the vacuole/lysosome. Autophagy-related 8 (Atg8), a core Atg protein essential for autophagosome formation, is a marker of several autophagic structures: the pre-autophagosomal structure (PAS), isolation membrane (IM), and autophagosome. Atg8 is conjugated to phosphatidylethanolamine (PE) through a ubiquitin-like conjugation system to yield Atg8-PE; this reaction is called Atg8 lipidation. Although the mechanisms of Atg8 lipidation have been well studied in vitro, the cellular locale of Atg8 lipidation remains enigmatic. Atg3 is an E2-like enzyme that catalyzes the conjugation reaction between Atg8 and PE. Therefore, we hypothesized that the localization of Atg3 would provide insights about the site of the lipidation reaction. To explore this idea, we constructed functional GFP-tagged Atg3 (Atg3-GFP) by inserting the GFP portion immediately after the handle region of Atg3. During autophagy, Atg3-GFP transiently formed a single dot per cell on the vacuolar membrane. This Atg3-GFP dot colocalized with 2× mCherry-tagged Atg8, demonstrating that Atg3 is localized to autophagic structures. Furthermore, we found that Atg3-GFP is localized to the IM by fine-localization analysis. The localization of Atg3 suggests that Atg3 plays an important role in autophagosome formation at the IM.