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.     

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.     

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.     

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.     

Morphogenesis beyond Cytokinetic Arrest in Saccharomyces cerevisiae

The budding yeast lyt1 mutation causes cell lysis. We report here that lyt1 is an allele of cdc15, a gene which encodes a protein kinase that functions late in the cell cycle. Neither cdc15-1 nor cdc15-lyt1 strains are able to septate at 37°C, even though they may manage to rebud. Cells lyse after a shmoo-like projection appears at the distal pole of the daughter cell. Actin polarizes towards the distal pole but the septins remain at the mother–daughter neck. This morphogenetic response reflects entry into a new round of the cell cycle: the preference for polarization from the distal pole was lost in bud1 cdc15 double mutants; double cdc15-lyt1 cdc28-4 mutants, defective for START, did not develop apical projections and apical polarization was accompanied by DNA replication. The same phenomena were caused by mutations in the genes CDC14, DBF2, and TEM1, which are functionally related to CDC15. Apical polarization was delayed in cdc15 mutants as compared with budding in control cells and this delay was abolished in a septin mutant. Our results suggest that the delayed M/G1 transition in cdc15 mutants is due to a septin-dependent checkpoint that couples initiation of the cell cycle to the completion of cytokinesis.     

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.     

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.     

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.     

Genome-Wide Fine-Scale Recombination Rate Variation in Drosophila melanogaster

Estimating fine-scale recombination maps of Drosophila from population genomic data is a challenging problem, in particular because of the high background recombination rate. In this paper, a new computational method is developed to address this challenge. Through an extensive simulation study, it is demonstrated that the method allows more accurate inference, and exhibits greater robustness to the effects of natural selection and noise, compared to a well-used previous method developed for studying fine-scale recombination rate variation in the human genome. As an application, a genome-wide analysis of genetic variation data is performed for two Drosophila melanogaster populations, one from North America (Raleigh, USA) and the other from Africa (Gikongoro, Rwanda). It is shown that fine-scale recombination rate variation is widespread throughout the D. melanogaster genome, across all chromosomes and in both populations. At the fine-scale, a conservative, systematic search for evidence of recombination hotspots suggests the existence of a handful of putative hotspots each with at least a tenfold increase in intensity over the background rate. A wavelet analysis is carried out to compare the estimated recombination maps in the two populations and to quantify the extent to which recombination rates are conserved. In general, similarity is observed at very broad scales, but substantial differences are seen at fine scales. The average recombination rate of the X chromosome appears to be higher than that of the autosomes in both populations, and this pattern is much more pronounced in the African population than the North American population. The correlation between various genomic features—including recombination rates, diversity, divergence, GC content, gene content, and sequence quality—is examined using the wavelet analysis, and it is shown that the most notable difference between D. melanogaster and humans is in the correlation between recombination and diversity.     

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.