Factors affecting inverted repeat stimulation of recombination and deletion in Saccharomyces cerevisiae

RPA and Rad27 limit templated and inverted insertions at DNA breaks

Formation of templated insertions at DNA double-strand breaks (DSBs) is very common in cancer cells. The mechanisms and enzymes regulating these events are largely unknown. Here, we investigated templated insertions in yeast at DSBs using amplicon sequencing across a repaired locus. We document very short (most ∼5-34 bp), templated inverted duplications at DSBs. They are generated through a foldback mechanism that utilizes microhomologies adjacent to the DSB. Enzymatic requirements suggest a hybrid mechanism wherein one end requires Polδ-mediated synthesis while the other end is captured by nonhomologous end joining (NHEJ). This process is exacerbated in mutants with low levels or mutated RPA ( rtt105 Δ; rfa1 -t33) or extensive resection mutant ( sgs1 Δ exo1 Δ). Templated insertions from various distant genomic locations also increase in these mutants as well as in rad27 Δ and originate from fragile regions of the genome. Among complex insertions, common events are insertions of two sequences, originating from the same locus and with inverted orientation. We propose that these inversions are also formed by microhomology-mediated template switching. Taken together, we propose that a shortage of RPA typical in cancer cells is one possible factor stimulating the formation of templated insertions.

A unifying model that explains the origins of human inverted copy number variants

With the release of the telomere-to-telomere human genome sequence and the availability of both long-read sequencing and optical genome mapping techniques, the identification of copy number variants (CNVs) and other structural variants is providing new insights into human genetic disease. Different mechanisms have been proposed to account for the novel junctions in these complex architectures, including aberrant forms of DNA replication, non-allelic homologous recombination, and various pathways that repair DNA breaks. Here, we have focused on a set of structural variants that include an inverted segment and propose that they share a common initiating event: an inverted triplication with long, unstable palindromic junctions. The secondary rearrangement of these palindromes gives rise to the various forms of inverted structural variants. We postulate that this same mechanism (ODIRA: origin-dependent inverted-repeat amplification) that creates the inverted CNVs in inherited syndromes also generates the palindromes found in cancers.

The evolutionary loss of the Eh1 motif in FoxE1 in the lineage of placental mammals

Forkhead box E1 (FoxE1) protein is a transcriptional regulator known to play a major role in the development of the thyroid gland. By performing sequence alignments, we detected a deletion in FoxE1, which occurred in the evolution of mammals, near the point of divergence of placental mammals. This deletion led to the loss of the majority of the Eh1 motif, which was important for interactions with transcriptional corepressors. To investigate a potential mechanism for this deletion, we analyzed replication through the deletion area in mammalian cells with two-dimensional gel electrophoresis, and in vitro, using a primer extension reaction. We demonstrated that the area of the deletion presented an obstacle for replication in both assays. The exact position of polymerization arrest in primer extension indicated that it was most likely caused by a quadruplex DNA structure. The quadruplex structure hypothesis is also consistent with the exact borders of the deletion. The exact roles of these evolutionary changes in FoxE1 family proteins are still to be determined.

Widely spaced and divergent inverted repeats become a potent source of chromosomal rearrangements in long single-stranded DNA regions

DNA inverted repeats (IRs) are widespread across many eukaryotic genomes. Their ability to form stable hairpin/cruciform secondary structures is causative in triggering chromosome instability leading to several human diseases. Distance and sequence divergence between IRs are inversely correlated with their ability to induce gross chromosomal rearrangements (GCRs) because of a lesser probability of secondary structure formation and chromosomal breakage. In this study, we demonstrate that structural parameters that normally constrain the instability of IRs are overcome when the repeats interact in single-stranded DNA (ssDNA). We established a system in budding yeast whereby >73 kb of ssDNA can be formed in cdc13-707fs mutants. We found that in ssDNA, 12 bp or 30 kb spaced Alu-IRs show similarly high levels of GCRs, while heterology only beyond 25% suppresses IR-induced instability. Mechanistically, rearrangements arise after cis-interaction of IRs leading to a DNA fold-back and the formation of a dicentric chromosome, which requires Rad52/Rad59 for IR annealing as well as Rad1-Rad10, Slx4, Msh2/Msh3 and Saw1 proteins for nonhomologous tail removal. Importantly, using structural characteristics rendering IRs permissive to DNA fold-back in yeast, we found that ssDNA regions mapped in cancer genomes contain a substantial number of potentially interacting and unstable IRs.

The Chloroplast Genome of the Lichen Photobiont Trebouxiophyceae sp. DW1 and Its Phylogenetic Implications

Lichens are symbiotic associations of algae and fungi. The genetic mechanism of the symbiosis of lichens and the influence of symbiosis on the size and composition of the genomes of symbiotic algae have always been intriguing scientific questions explored by lichenologists. However, there were limited data on lichen genomes. Therefore, we isolated and purified a lichen symbiotic alga to obtain a single strain (Trebouxiophyceae sp. DW1), and then obtained its chloroplast genome information by next-generation sequencing (NGS). The chloroplast genome is 129,447 bp in length, and the GC content is 35.2%. Repetitive sequences with the length of 30–35 bp account for 1.27% of the total chloroplast genome. The simple sequence repeats are all mononucleotide repeats. Codon usage analysis showed that the genome tended to use codon ending in A/U. By comparing the length of different regions of Trebouxiophyceae genomes, we found that the changes in the length of exons, introns, and intergenic sequences affect the size of genomes. Trebouxiophyceae had an unstable chloroplast genome structure, with IRs repeatedly losing during evolution. Phylogenetic analysis showed that Trebouxiophyceae is paraphyletic, and Trebouxiophyceae sp. DW1 is sister to the clade of Koliella longiseta and Pabia signiensis.

High-throughput techniques enable advances in the roles of DNA and RNA secondary structures in transcriptional and post-transcriptional gene regulation

The most stable structure of DNA is the canonical right-handed double helix termed B DNA. However, certain environments and sequence motifs favor alternative conformations, termed non-canonical secondary structures. The roles of DNA and RNA secondary structures in transcriptional regulation remain incompletely understood. However, advances in high-throughput assays have enabled genome wide characterization of some secondary structures. Here, we describe their regulatory functions in promoters and 3’UTRs, providing insights into key mechanisms through which they regulate gene expression. We discuss their implication in human disease, and how advances in molecular technologies and emerging high-throughput experimental methods could provide additional insights.

Migrating bubble synthesis promotes mutagenesis through lesions in its template

Break-induced replication (BIR) proceeds via a migrating D-loop for hundreds of kilobases and is highly mutagenic. Previous studies identified long single-stranded (ss) nascent DNA that accumulates during leading strand synthesis to be a target for DNA damage and a primary source of BIR-induced mutagenesis. Here, we describe a new important source of mutagenic ssDNA formed during BIR: the ssDNA template for leading strand BIR synthesis formed during D-loop migration. Specifically, we demonstrate that this D-loop bottom template strand (D-BTS) is susceptible to APOBEC3A (A3A)-induced DNA lesions leading to mutations associated with BIR. Also, we demonstrate that BIR-associated ssDNA promotes an additional type of genetic instability: replication slippage between microhomologies stimulated by inverted DNA repeats. Based on our results we propose that these events are stimulated by both known sources of ssDNA formed during BIR, nascent DNA formed by leading strand synthesis, and the D-BTS that we describe here. Together we report a new source of mutagenesis during BIR that may also be shared by other homologous recombination pathways driven by D-loop repair synthesis.

Nucleic acid folding simulations using a physics-based atomistic free energy model

Performing full-resolution atomistic simulations of nucleic acid folding has remained a challenge for biomolecular modeling. Understanding how nucleic acids fold and how they transition between different folded structures as they unfold and refold has important implications for biology. This paper reports a theoretical model and computer simulation of the ab initio folding of DNA inverted repeat sequences. The formulation is based on an all-atom conformational model of the sugar-phosphate backbone via chain closure, and it incorporates three major molecular-level driving forces-base stacking, counterion-induced backbone self-interactions, and base pairing-via separate analytical theories designed to capture and reproduce the effects of the solvent without requiring explicit water and ions in the simulation. To accelerate computational throughput, a mixed numerical/analytical algorithm for the calculation of the backbone conformational volume is incorporated into the Monte Carlo simulation, and special stochastic sampling techniques were employed to achieve the computational efficiency needed to fold nucleic acids from scratch. This paper describes implementation details, benchmark results, and the advantages and technical challenges with this approach.

High-throughput characterization of the role of non-B DNA motifs on promoter function

lternative DNA conformations, termed non-B DNA structures, can affect transcription, but the underlying mechanisms and their functional impact have not been systematically characterized. Here, we used computational genomic analyses coupled with massively parallel reporter assays (MPRAs) to show that certain non-B DNA structures have a substantial effect on gene expression. Genomic analyses found that non-B DNA structures at promoters harbor an excess of germline variants. Analysis of multiple MPRAs, including a promoter library specifically designed to perturb non-B DNA structures, functionally validated that Z-DNA can significantly affect promoter activity. We also observed that biophysical properties of non-B DNA motifs, such as the length of Z-DNA motifs and the orientation of G-quadruplex structures relative to transcriptional direction, have a significant effect on promoter activity. Combined, their higher mutation rate and functional effect on transcription implicate a subset of non-B DNA motifs as major drivers of human gene-expression-associated phenotypes.

Protein innovation through template switching in the Saccharomyces cerevisiae lineage

DNA polymerase template switching between short, non-identical inverted repeats (IRs) is a genetic mechanism that leads to the homogenization of IR arms and to IR spacer inversion, which cause multinucleotide mutations (MNMs). It is unknown if and how template switching affects gene evolution. In this study, we performed a phylogenetic analysis to determine the effect of template switching between IR arms on coding DNA of Saccharomyces cerevisiae. To achieve this, perfect IRs that co-occurred with MNMs between a strain and its parental node were identified in S. cerevisiae strains. We determined that template switching introduced MNMs into 39 protein-coding genes through S. cerevisiae evolution, resulting in both arm homogenization and inversion of the IR spacer. These events in turn resulted in nonsynonymous substitutions and up to five neighboring amino acid replacements in a single gene. The study demonstrates that template switching is a powerful generator of multiple substitutions within codons. Additionally, some template switching events occurred more than once during S. cerevisiae evolution. Our findings suggest that template switching constitutes a general mutagenic mechanism that results in both nonsynonymous substitutions and parallel evolution, which are traditionally considered as evidence for positive selection, without the need for adaptive explanations.

Mitochondrial Genomics of Six Cacao Pathogens From the Basidiomycete Family Marasmiaceae

Thread blight disease has recently been described as an emerging disease on cacao (Theobroma cacao) in Ghana. In Ghana, thread blight disease is caused by multiple species of the Marasmiaceae family: Marasmius tenuissimus, M. crinis-equi, M. palmivorus, and Marasmiellus scandens. Interestingly, two additional members of the Marasmiaceae; Moniliophthora roreri (frosty pod rot) and Moniliophthora perniciosa (witches’ broom disease), are major pathogens of cacao in the Western hemisphere. It is important to accurately characterize the genetic relationships among these economically important species in support of their disease management. We used data from Illumina NGS-based genome sequencing efforts to study the mitochondrial genomes (mitogenomes) of the four cacao thread blight associated pathogens from Ghana and compared them with published mitogenomes of Mon. roreri and Mon. perniciosa. There is a remarkable interspecies variation in mitogenome size within the six cacao-associated Marasmiaceae species, ranging from 43,121 to 109,103 bp. The differences in genome lengths are primarily due to the number and lengths of introns, differences in intergenic space, and differences in the size and numbers of unidentified ORFs (uORF). Among seven M. tenuissimus mitogenomes sequenced, there is variation in size and sequence pointing to divergent evolution patterns within the species. The intronic regions show a high degree of sequence variation compared to the conserved sequences of the 14 core genes. The intronic ORFs identified, regardless of species, encode GIY-YIG or LAGLIDADG domain-containing homing endonuclease genes. Phylogenetic relationships using the 14 core proteins largely mimic the phylogenetic relationships observed in gene order patterns, grouping M. tenuissimus with M. crinis-equi, and M. palmivorus with Mon. roreri and Mon. perniciosa, leaving Mar. scandens as an outlier. The results from this study provide evidence of independent expansion/contraction events and sequence diversification in each species and establish a foundation for further exploration of the evolutionary trajectory of the fungi in Marasmiaceae family.

Nanopore Sequencing Indicates That Tandem Amplification of Chromosome 20q11.21 in Human Pluripotent Stem Cells Is Driven by Break-Induced Replication

Copy number variants (CNVs) are genomic rearrangements implicated in numerous congenital and acquired diseases, including cancer. The appearance of culture-acquired CNVs in human pluripotent stem cells (PSCs) has prompted concerns for their use in regenerative medicine. A particular problem in PSC is the frequent occurrence of CNVs in the q11.21 region of chromosome 20. However, the exact mechanism of origin of this amplicon remains elusive due to the difficulty in delineating its sequence and breakpoints. Here, we have addressed this problem using long-read Nanopore sequencing of two examples of this CNV, present as duplication and as triplication. In both cases, the CNVs were arranged in a head-to-tail orientation, with microhomology sequences flanking or overlapping the proximal and distal breakpoints. These breakpoint signatures point to a mechanism of microhomology-mediated break-induced replication in CNV formation, with surrounding Alu sequences likely contributing to the instability of this genomic region.

Structural parameters of palindromic repeats determine the specificity of nuclease attack of secondary structures

Palindromic sequences are a potent source of chromosomal instability in many organisms and are implicated in the pathogenesis of human diseases. In this study, we investigate which nucleases are responsible for cleavage of the hairpin and cruciform structures and generation of double-strand breaks at inverted repeats in Saccharomyces cerevisiae. We demonstrate that the involvement of structure-specific nucleases in palindrome fragility depends on the distance between inverted repeats and their transcriptional status. The attack by the Mre11 complex is constrained to hairpins with loops <9 nucleotides. This restriction is alleviated upon RPA depletion, indicating that RPA controls the stability and/or formation of secondary structures otherwise responsible for replication fork stalling and DSB formation. Mus81-Mms4 cleavage of cruciforms occurs at divergently but not convergently transcribed or nontranscribed repeats. Our study also reveals the third pathway for fragility at perfect and quasi-palindromes, which involves cruciform resolution during the G2 phase of the cell cycle.

Molluscan mitochondrial genomes break the rules

The first animal mitochondrial genomes to be sequenced were of several vertebrates and model organisms, and the consistency of genomic features found has led to a 'textbook description'. However, a more broad phylogenetic sampling of complete animal mitochondrial genomes has found many cases where these features do not exist, and the phylum Mollusca is especially replete with these exceptions. The characterization of full mollusc mitogenomes required considerable effort involving challenging molecular biology, but has created an enormous catalogue of surprising deviations from that textbook description, including wide variation in size, radical genome rearrangements, gene duplications and losses, the introduction of novel genes, and a complex system of inheritance dubbed 'doubly uniparental inheritance'. Here, we review the extraordinary variation in architecture, molecular functioning and intergenerational transmission of molluscan mitochondrial genomes. Such features represent a great potential for the discovery of biological history, processes and functions that are novel for animal mitochondrial genomes. This provides a model system for studying the evolution and the manifold roles that mitochondria play in organismal physiology, and many ways that the study of mitochondrial genomes are useful for phylogeny and population biology. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Palindromes in DNA-A Risk for Genome Stability and Implications in Cancer

A palindrome in DNA consists of two closely spaced or adjacent inverted repeats. Certain palindromes have important biological functions as parts of various cis-acting elements and protein binding sites. However, many palindromes are known as fragile sites in the genome, sites prone to chromosome breakage which can lead to various genetic rearrangements or even cell death. The ability of certain palindromes to initiate genetic recombination lies in their ability to form secondary structures in DNA which can cause replication stalling and double-strand breaks. Given their recombinogenic nature, it is not surprising that palindromes in the human genome are involved in genetic rearrangements in cancer cells as well as other known recurrent translocations and deletions associated with certain syndromes in humans. Here, we bring an overview of current understanding and knowledge on molecular mechanisms of palindrome recombinogenicity and discuss possible implications of DNA palindromes in carcinogenesis. Furthermore, we overview the data on known palindromic sequences in the human genome and efforts to estimate their number and distribution, as well as underlying mechanisms of genetic rearrangements specific palindromic sequences cause.

Variable Number Tandem Repeats in the Mitochondrial DNA of Lentinula edodes

Variable number tandem repeats (VNTRs) in mitochondrial DNA (mtDNA) of Lentinula edodes are of interest for their role in mtDNA variation and their application as genetic marker. Sequence analysis of three L. edodes mtDNAs revealed the presence of VNTRs of two categories. Type I VNTRs consist of two types of repeat units in a symmetric distribution, whereas Type II VNTRs contain tandemly arrayed repeats of 7- or 17-bp DNA sequences. The number of repeat units was variable depending on the mtDNA of different strains. Using the variations in VNTRs as a mitochondrial marker and the A mating type as a nuclear type marker, we demonstrated that one of the two nuclei in the donor dikaryon preferentially enters into the monokaryotic cytoplasm to establish a new dikaryon which still retains the mitochondria of the monokaryon in the individual mating. Interestingly, we found 6 VNTRs with newly added repeat units from the 22 mates, indicating that elongation of VNTRs occurs during replication of mtDNA. This, together with comparative analysis of the repeating pattern, enables us to propose a mechanistic model that explains the elongation of Type I VNTRs through reciprocal incorporation of basic repeat units, 5’-TCCCTTTAGGG-3’ and its complementary sequence (5’-CCCTAAAGGGA-3’).

Genome plasticity in Candida albicans is driven by long repeat sequences

Genome rearrangements resulting in copy number variation (CNV) and loss of heterozygosity (LOH) are frequently observed during the somatic evolution of cancer and promote rapid adaptation of fungi to novel environments. In the human fungal pathogen Candida albicans, CNV and LOH confer increased virulence and antifungal drug resistance, yet the mechanisms driving these rearrangements are not completely understood. Here, we unveil an extensive array of long repeat sequences (65-6499 bp) that are associated with CNV, LOH, and chromosomal inversions. Many of these long repeat sequences are uncharacterized and encompass one or more coding sequences that are actively transcribed. Repeats associated with genome rearrangements are predominantly inverted and separated by up to ~1.6 Mb, an extraordinary distance for homology-based DNA repair/recombination in yeast. These repeat sequences are a significant source of genome plasticity across diverse strain backgrounds including clinical, environmental, and experimentally evolved isolates, and represent previously uncharacterized variation in the reference genome.

The role of fork stalling and DNA structures in causing chromosome fragility

Alternative non-B form DNA structures, also called secondary structures, can form in certain DNA sequences under conditions that produce single-stranded DNA, such as during replication, transcription, and repair. Direct links between secondary structure formation, replication fork stalling, and genomic instability have been found for many repeated DNA sequences that cause disease when they expand. Common fragile sites (CFSs) are known to be AT-rich and break under replication stress, yet the molecular basis for their fragility is still being investigated. Over the past several years, new evidence has linked both the formation of secondary structures and transcription to fork stalling and fragility of CFSs. How these two events may synergize to cause fragility and the role of nuclease cleavage at secondary structures in rare and CFSs are discussed here. We also highlight evidence for a new hypothesis that secondary structures at CFSs not only initiate fragility but also inhibit healing, resulting in their characteristic appearance.

Noncanonical secondary structures arising from non-B DNA motifs are determinants of mutagenesis

Somatic mutations show variation in density across cancer genomes. Previous studies have shown that chromatin organization and replication time domains are correlated with, and thus predictive of, this variation. Here, we analyze 1809 whole-genome sequences from 10 cancer types to show that a subset of repetitive DNA sequences, called non-B motifs that predict noncanonical secondary structure formation can independently account for variation in mutation density. Combined with epigenetic factors and replication timing, the variance explained can be improved to 43%-76%. Approximately twofold mutation enrichment is observed directly within non-B motifs, is focused on exposed structural components, and is dependent on physical properties that are optimal for secondary structure formation. Therefore, there is mounting evidence that secondary structures arising from non-B motifs are not simply associated with increased mutation density-they are possibly causally implicated. Our results suggest that they are determinants of mutagenesis and increase the likelihood of recurrent mutations in the genome. This analysis calls for caution in the interpretation of recurrent mutations and highlights the importance of taking non-B motifs that can simply be inferred from the reference sequence into consideration in background models of mutability henceforth.

Single-strand annealing between inverted DNA repeats: Pathway choice, participating proteins, and genome destabilizing consequences

Double strand DNA breaks (DSBs) are dangerous events that can result from various causes including environmental assaults or the collapse of DNA replication. While the efficient and precise repair of DSBs is essential for cell survival, faulty repair can lead to genetic instability, making the choice of DSB repair an important step. Here we report that inverted DNA repeats (IRs) placed near a DSB can channel its repair from an accurate pathway that leads to gene conversion to instead a break-induced replication (BIR) pathway that leads to genetic instabilities. The effect of IRs is explained by their ability to form unusual DNA structures when present in ssDNA that is formed by DSB resection. We demonstrate that IRs can form two types of unusual DNA structures, and the choice between these structures depends on the length of the spacer separating IRs. In particular, IRs separated by a long (1-kb) spacer are predominantly involved in inter-molecular single-strand annealing (SSA) leading to the formation of inverted dimers; IRs separated by a short (12-bp) spacer participate in intra-molecular SSA, leading to the formation of fold-back (FB) structures. Both of these structures interfere with an accurate DSB repair by gene conversion and channel DSB repair into BIR, which promotes genomic destabilization. We also report that different protein complexes participate in the processing of FBs containing short (12-bp) versus long (1-kb) ssDNA loops. Specifically, FBs with short loops are processed by the MRX-Sae2 complex, whereas the Rad1-Rad10 complex is responsible for the processing of long loops. Overall, our studies uncover the mechanisms of genomic destabilization resulting from re-routing DSB repair into unusual pathways by IRs. Given the high abundance of IRs in the human genome, our findings may contribute to the understanding of IR-mediated genomic destabilization associated with human disease.

Efficient and accurate repair of double-strand DNA breaks (DSBs), resulting from the exposure of cells to ionizing radiation or various chemicals, is crucial for cell survival. Conversely, faulty DSB repair can generate genomic instability that can lead to birth defects or cancer in humans. Here we demonstrate that inverted DNA repeats (IRs) placed in the vicinity of a DSB, interfere with the accurate repair of DSBs and promote genomic rearrangements and chromosome loss. This results from annealing between inverted repeats, located either in different DNA molecules or in the same molecule. In addition, we describe a new role for the Rad1-Rad10 protein complex in processing fold-back (FB) structures formed by intra-molecular annealing involving IRs separated by long spacers. In contrast, FBs with short spacers are processed by the Mre11-Rad50-Xrs2/-Sae2 complex. Overall, we describe several pathways of DSB promoted interaction between IRs that can lead to genomic instability. Given the large number of IRs in the human genome, our findings are relevant to the mechanisms driving genomic destabilization in humans contributing to the development of cancer and other diseases.

Saccharomyces cerevisiae Red1 protein exhibits nonhomologous DNA end-joining activity and potentiates Hop1-promoted pairing of double-stranded DNA

Elucidation of the function of synaptonemal complex (SC) in Saccharomyces cerevisiae has mainly focused on in vivo analysis of recombination-defective meiotic mutants. Consequently, significant gaps remain in the mechanistic understanding of the activities of various SC proteins and the functional relationships among them. S. cerevisiae Hop1 and Red1 are essential structural components of the SC axial/lateral elements. Previous studies have demonstrated that Hop1 is a structure-selective DNA-binding protein exhibiting high affinity for the Holliday junction and promoting DNA bridging, condensation, and pairing between double-stranded DNA molecules. However, the exact mode of action of Red1 remains unclear, although it is known to interact with Hop1 and to suppress the spore viability defects of hop1 mutant alleles. Here, we report the purification and functional characterization of the full-length Red1 protein. Our results revealed that Red1 forms a stable complex with Hop1 in vitro and provided quantitative insights into their physical interactions. Mechanistically, Red1 preferentially associated with the Holliday junction and 3-way junction rather than with single- or double-stranded DNA with overhangs. Although Hop1 and Red1 exhibited similar binding affinities toward several DNA substrates, the two proteins displayed some significant differences. Notably, Red1, by itself, lacked DNA-pairing ability; however, it potentiated Hop1-promoted intermolecular pairing between double-stranded DNA molecules. Moreover, Red1 exhibited nonhomologous DNA end-joining activity, thus revealing an unexpected role for Red1 in recombination-based DNA repair. Collectively, this study presents the first direct insights into Red1's mode of action and into the mechanism underlying its role in chromosome synapsis and recombination.

Lirex: A Package for Identification of Long Inverted Repeats in Genomes

Long inverted repeats (LIRs) are evolutionarily and functionally important structures in genomes because of their involvement in RNA interference, DNA recombination, and gene duplication. Identification of LIRs is highly complicated when mismatches and indels between the repeats are permitted. Long inverted repeat explorer (Lirex) was developed and introduced in this report. Written in Java, Lirex provides a user-friendly interface and allows users to specify LIR searching criteria, such as length of the region, as well as pattern and size of the repeats. Recombinogenic LIRs can be selected on the basis of mismatch rate and internal spacer size from identified LIRs. Lirex, as a cross-platform tool to identify LIRs in a genome, may assist in designing following experiments to explore the function of LIRs. Our tool can identify more LIRs than other LIR searching tools. Lirex is publicly available at http://124.16.219.129/Lirex.

Characterization and distribution of mating-type genes of the turfgrass pathogen Sclerotinia homoeocarpa on a global scale

Sclerotinia homoeocarpa F.T. Bennett is a filamentous member of Ascomycota that causes dollar spot, the most economically important disease of turfgrass worldwide. We sequenced and characterized the mating-type (MAT) locus of four recently-collected contemporary strains causing dollar spot, four historical type strains used to describe the fungus, and three species of Rutstroemiaceae. Moreover, we developed a multiplex PCR assay to screen 1019 contemporary isolates for mating-type. The organization of the MAT loci of all strains examined could be classified into one of four categories: (1) putatively heterothallic, as exemplified by all contemporary strains and three of four historical type strains; (2) putatively heterothallic with a deleted putative gene in the MAT1-2 idiomorph, as detected in strains from two recently-collected populations in the United Kingdom that show more similarity to historical strains; (3) putatively homothallic with close physical linkage between MAT1-1-1 and MAT1-2-1, as found in one historical type strain of S. homoeocarpa and two strains of Rutstroemia cuniculi; and (4) an unresolved but apparently homothallic organization in which strains contained both MAT1-1-1 and MAT1-2-1 but linkage between these genes and between the two flanking genes could not be confirmed, as identified in R. paludosa and Poculum henningsianum. In contemporary S. homoeocarpa populations there was no significant difference in the frequency of the two mating types in clone-corrected samples when analyzed on regional and local scales, suggesting sex may be possible in this pathogen. However, two isolates from Italy and twenty from California were heterokaryotic for both complete heterothallic MAT idiomorphs. Results from this study contribute to knowledge about mating systems in filamentous fungi and enhance our understanding of the evolution and biology of an important plant pathogen.

Correlations between long inverted repeat (LIR) features, deletion size and distance from breakpoint in human gross gene deletions

Long inverted repeats (LIRs) have been shown to induce genomic deletions in yeast. In this study, LIRs were investigated within ±10 kb spanning each breakpoint from 109 human gross deletions, using Inverted Repeat Finder (IRF) software. LIR number was significantly higher at the breakpoint regions, than in control segments (P < 0.001). In addition, it was found that strong correlation between 5' and 3' LIR numbers, suggesting contribution to DNA sequence evolution (r = 0.85, P < 0.001). 138 LIR features at ±3 kb breakpoints in 89 (81%) of 109 gross deletions were evaluated. Significant correlations were found between distance from breakpoint and loop length (r = -0.18, P < 0.05) and stem length (r = -0.18, P < 0.05), suggesting DNA strands are potentially broken in locations closer to bigger LIRs. In addition, bigger loops cause larger deletions (r = 0.19, P < 0.05). Moreover, loop length (r = 0.29, P < 0.02) and identity between stem copies (r = 0.30, P < 0.05) of 3' LIRs were more important in larger deletions. Consequently, DNA breaks may form via LIR-induced cruciform structure during replication. DNA ends may be later repaired by non-homologous end-joining (NHEJ), with following deletion.

High variability of mitochondrial gene order among fungi

From their origin as an early alpha proteobacterial endosymbiont to their current state as cellular organelles, large-scale genomic reorganization has taken place in the mitochondria of all main eukaryotic lineages. So far, most studies have focused on plant and animal mitochondrial (mt) genomes (mtDNA), but fungi provide new opportunities to study highly differentiated mtDNAs. Here, we analyzed 38 complete fungal mt genomes to investigate the evolution of mtDNA gene order among fungi. In particular, we looked for evidence of nonhomologous intrachromosomal recombination and investigated the dynamics of gene rearrangements. We investigated the effect that introns, intronic open reading frames (ORFs), and repeats may have on gene order. Additionally, we asked whether the distribution of transfer RNAs (tRNAs) evolves independently to that of mt protein-coding genes. We found that fungal mt genomes display remarkable variation between and within the major fungal phyla in terms of gene order, genome size, composition of intergenic regions, and presence of repeats, introns, and associated ORFs. Our results support previous evidence for the presence of mt recombination in all fungal phyla, a process conspicuously lacking in most Metazoa. Overall, the patterns of rearrangements may be explained by the combined influences of recombination (i.e., most likely nonhomologous and intrachromosomal), accumulated repeats, especially at intergenic regions, and to a lesser extent, mobile element dynamics.

Excision of unstable artificial gene-specific inverted repeats mediates scar-free gene deletions in Escherichia coli

Inverted repeat and palindromic sequences have the propensity to form non-beta cruciform structures during DNA replication, leading to perturbations within the genome or plasmid replicon. In this study, the tolerance of the Escherichia coli genome to inverted repeat sequences from 25 to 1200 bp was investigated. Genomic inverted repeats were readily created via the homologous insertion of an overlap extension PCR product containing a gene-specific region of the genome together with thyA coding sequence, creating inverted repeat sequences of various lengths flanking the thyA selection marker in the resulting genome. Inverted repeat sequences below 100 bp were stably propagated, while those above and up to 1200 bp were found to be transiently unstable under auxotrophic thymine selection. Excision efficiency improves with increases of the inverted repeat until 600-800 bp, indicating that the genomic stability of inverted repeat sequences is due to secondary structure formation. Its effectiveness of creating precise and scar-free gene deletions was further demonstrated by deleting a number of genes in E. coli. The procedure can be readily adapted for sequence integration and point mutations in E. coli genome. It also has the potential for applications on other bacteria for efficient gene deletions.

Genome-wide high-resolution mapping of chromosome fragile sites in Saccharomyces cerevisiae

In mammalian cells, perturbations in DNA replication result in chromosome breaks in regions termed "fragile sites." Using DNA microarrays, we mapped recombination events and chromosome rearrangements induced by reduced levels of the replicative DNA polymerase-α in the yeast Saccharomyces cerevisiae. We found that the recombination events were nonrandomly associated with a number of structural/sequence motifs that correlate with paused DNA replication forks, including replication-termination sites (TER sites) and binding sites for the helicase Rrm3p. The pattern of gene-conversion events associated with cross-overs suggests that most of the DNA lesions that initiate recombination between homologs are double-stranded DNA breaks induced during S or G2 of the cell cycle, in contrast to spontaneous recombination events that are initiated by double-stranded DNA breaks formed prior to replication. Low levels of DNA polymerase-α also induced very high rates of aneuploidy, as well as chromosome deletions and duplications. Most of the deletions and duplications had Ty retrotransposons at their breakpoints.

Genome-wide screen reveals replication pathway for quasi-palindrome fragility dependent on homologous recombination

Inverted repeats capable of forming hairpin and cruciform structures present a threat to chromosomal integrity. They induce double strand breaks, which lead to gross chromosomal rearrangements, the hallmarks of cancers and hereditary diseases. Secondary structure formation at this motif has been proposed to be the driving force for the instability, albeit the mechanisms leading to the fragility are not well-understood. We carried out a genome-wide screen to uncover the genetic players that govern fragility of homologous and homeologous Alu quasi-palindromes in the yeast Saccharomyces cerevisiae. We found that depletion or lack of components of the DNA replication machinery, proteins involved in Fe-S cluster biogenesis, the replication-pausing checkpoint pathway, the telomere maintenance complex or the Sgs1-Top3-Rmi1 dissolvasome augment fragility at Alu-IRs. Rad51, a component of the homologous recombination pathway, was found to be required for replication arrest and breakage at the repeats specifically in replication-deficient strains. These data demonstrate that Rad51 is required for the formation of breakage-prone secondary structures in situations when replication is compromised while another mechanism operates in DSB formation in replication-proficient strains.

Inverted repeats are found in many eukaryotic genomes including humans. They have a potential to cause chromosomal breakage and rearrangements that contribute to genome polymorphism and the development of diseases. Instability of inverted repeats is accounted for by their propensity to adopt DNA secondary structures that is negatively affected by the distance between the repeats and level of sequence divergence. However, the genetic factors that promote the abnormal structure formation or affect the ability of the repeats to break are largely unknown. Here, using a genome-wide screen we identified 38 mutants that destabilize imperfect human inverted Alu repeats and predispose them to breakage. The proteins that are required to maintain repeat stability belong to the core of the DNA replication machinery and to the accessory proteins that help replication fork to move through the difficult templates. Remarkably, when replication machinery is compromised, the proteins involved in homologous recombination promote the formation of secondary structures and replication block thereby triggering breakage at the inverted repeats. These results reveal a powerful pathway for the destabilization of chromosomes containing inverted repeats that requires the activity of homologous recombination.

Comparative analysis of the inverted repeat of a chalcone synthase pseudogene between yellow soybean and seed coat pigmented mutants

In soybean, the I gene inhibits pigmentation over the entire seed coat, resulting in yellow seeds. It is thought that this suppression of seed coat pigmentation is due to naturally occurring RNA silencing of chalcone synthase genes (CHS silencing). Fully pigmented seeds can be found among harvested yellow seeds at a very low percentage. These seed coat pigmented (scp) mutants are generated from yellow soybeans by spontaneous recessive mutation of the I gene. A candidate for the I gene, GmIRCHS, contains a perfect inverted repeat (IR) of a CHS pseudogene (pseudoCHS3) and transcripts of GmIRCHS form a double-stranded CHS RNA that potentially triggers CHS silencing. One CHS gene, ICHS1, is located 680 bp downstream of GmIRCHS. Here, the GmIRCHS-ICHS1 cluster was compared in scp mutants of various origins. In these mutants, sequence divergence in the cluster resulted in complete or partial loss of GmIRCHS in at least the pseudoCHS3 region. This result is consistent with the notion that the IR of pseudoCHS3 is sufficient to induce CHS silencing, and further supports that GmIRCHS is the I gene.

Mitochondrial inverted repeats strongly correlate with lifespan: mtDNA inversions and aging

Mitochondrial defects are implicated in aging and in a multitude of age-related diseases, such as cancer, heart failure, Parkinson's disease, and Huntington's disease. However, it is still unclear how mitochondrial defects arise under normal physiological conditions. Mitochondrial DNA (mtDNA) deletions caused by direct repeats (DRs) are implicated in the formation of mitochondrial defects, however, mitochondrial DRs show relatively weak (Pearson's r = -0.22, p<0.002; Spearman's ρ = -0.12, p = 0.1) correlation with maximum lifespan (MLS). Here we report a stronger correlation (Pearson's r = -0.55, p<10(-16); Spearman's ρ = -0.52, p<10(-14)) between mitochondrial inverted repeats (IRs) and lifespan across 202 species of mammals. We show that, in wild type mice under normal conditions, IRs cause inversions, which arise by replication-dependent mechanism. The inversions accumulate with age in the brain and heart. Our data suggest that IR-mediated inversions are more mutagenic than DR-mediated deletions in mtDNA, and impose stronger constraint on lifespan. Our study identifies IR-induced mitochondrial genome instability during mtDNA replication as a potential cause for mitochondrial defects.

Size of gene specific inverted repeat--dependent gene deletion In Saccharomyces cerevisiae

We describe here an approach for rapidly producing scar-free and precise gene deletions in S. cerevisiae with high efficiency. Preparation of the disruption gene cassette in this approach was simply performed by overlap extension-PCR of an invert repeat of a partial or complete sequence of the targeted gene with URA3. Integration of the prepared disruption gene cassette to the designated position of a target gene leads to the formation of a mutagenesis cassette within the yeast genome, which consists of a URA3 gene flanked by the targeted gene and its inverted repeat between two short identical direct repeats. The inherent instability of the inverted sequences in close proximity facilitates the self-excision of the entire mutagenesis cassette deposited in the genome and promotes homologous recombination resulting in a seamless deletion via a single transformation. This rapid assembly circumvents the difficulty during preparation of disruption gene cassettes composed of two inverted repeats of the URA3, which requires the engineering of unique restriction sites for subsequent digestion and T4 DNA ligation in vitro. We further identified that the excision of the entire mutagenesis cassette flanked by two DRs in the transformed S. cerevisiae is dependent on the length of the inverted repeat of which a minimum of 800 bp is required for effective gene deletion. The deletion efficiency improves with the increase of the inverted repeat till 1.2 kb. Finally, the use of gene-specific inverted repeats of target genes enables simultaneous gene deletions. The procedure has the potential for application on other yeast strains to achieve precise and efficient removal of gene sequences.

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.

Double-strand breaks (DSBs) are DNA lesions that can be fatal to a cell if left unrepaired. They can be caused by exogenous sources, such as gamma radiation, or endogenous stresses, such as high levels of transcription. Yeast cells primarily repair DSBs that are initiated outside of meiosis by mitotic recombination, which can result in physical exchanges between chromosomes, known as crossovers. We created a mitotic recombination map of one chromosome arm, representing 10% of the genome. This recombination map allows us to determine which regions of the chromosome arm are more susceptible to DNA damage than other regions. We were able to determine that most DSBs that result in detectable genomic changes were initiated prior to DNA replication and that some secondary DNA structures can be recombination hotspots. Recombination can also occur during meiosis, as a method of ensuring proper chromosome segregation. However, previously reported meiotic recombination maps have no correlation with our mitotic recombination map.

Compositional bias is a major determinant of the distribution pattern and abundance of palindromes in Drosophila melanogaster

Palindromic sequences are important DNA motifs related to gene regulation, DNA replication and recombination, and thus, investigating the evolutionary forces shaping the distribution pattern and abundance of palindromes in the genome is substantially important. In this article, we analyzed the abundance of palindromes in the genome, and then explored the possible effects of several genomic factors on the palindrome distribution and abundance in Drosophila melanogaster. Our results show that the palindrome abundance in D. melanogaster deviates from random expectation and the uneven distribution of palindromes across the genome is associated with local GC content, recombination rate, and coding exon density. Our data suggest that base composition is the major determinant of the distribution pattern and abundance of palindromes and the correlation between palindrome density and recombination is a side-product of the effect of compositional bias on the palindrome abundance.

DNA secondary structures: stability and function of G-quadruplex structures

In addition to the canonical double helix, DNA can fold into various other inter- and intramolecular secondary structures. Although many such structures were long thought to be in vitro artefacts, bioinformatics demonstrates that DNA sequences capable of forming these structures are conserved throughout evolution, suggesting the existence of non-B-form DNA in vivo. In addition, genes whose products promote formation or resolution of these structures are found in diverse organisms, and a growing body of work suggests that the resolution of DNA secondary structures is critical for genome integrity. This Review focuses on emerging evidence relating to the characteristics of G-quadruplex structures and the possible influence of such structures on genomic stability and cellular processes, such as transcription.

Localized DNA cleavage secondary to genotoxic exposure adjacent to an Alu inverted repeat

Radiation is a potent inducer of DNA damage leading to both random DNA loss and mutation. As part of a study focused on the mechanism whereby cells undergo loss of heterozygosity (LOH), a region of common LOH telomeric termination at 11q24 was observed in clones of H292 mucoepidermoid cells established after irradiation (IR). A 10-kbp region including the telomeric extent of LOH termination was analyzed after IR using six sets of ligation-mediated polymerase chain reaction (PCR) primers to detect the presence of DNA breaks. A cluster of DNA breaks was detected that closely mapped to the telomeric extent of LOH and which were observed up to 8 hr after IR. Repeating the experiment in the presence of the inhibitor of apoptosis, zVAD.fmk, did not change the location or amount of cleavage. A similar distribution of breaks was also seen in the MCF-10A breast cancer cell line after IR. Further inspection of the involved region showed that 22/32 and 7/7 DNA breaks found in H292 and MCF-10A cells, respectively, were located either in or immediately adjacent to an AluSx1 sequence, itself ≈ 1 kbp 5' to an AluSq2 that was in an inverted orientation to the AluSx1. The region between the inverted Alu repeats was notable for both DNAse hypersensitivity and an open chromatin conformation inferred from histone modification data. These factors may contribute to genomic instability at this location.

FLNA genomic rearrangements cause periventricular nodular heterotopia

OBJECTIVE

To identify copy number variant (CNV) causes of periventricular nodular heterotopia (PNH) in patients for whom FLNA sequencing is negative.

METHODS

Screening of 35 patients from 33 pedigrees on an Affymetrix 6.0 microarray led to the identification of one individual bearing a CNV that disrupted FLNA. FLNA-disrupting CNVs were also isolated in 2 other individuals by multiplex ligation probe amplification. These 3 cases were further characterized by high-resolution oligo array comparative genomic hybridization (CGH), and the precise junctional breakpoints of the rearrangements were identified by PCR amplification and sequencing.

RESULTS

We report 3 cases of PNH caused by nonrecurrent genomic rearrangements that disrupt one copy of FLNA. The first individual carried a 113-kb deletion that removes all but the first exon of FLNA. A second patient harbored a complex rearrangement including a deletion of the 3' end of FLNA accompanied by a partial duplication event. A third patient bore a 39-kb deletion encompassing all of FLNA and the neighboring gene EMD. High-resolution oligo array CGH of the FLNA locus suggests distinct molecular mechanisms for each of these rearrangements, and implicates nearby low copy repeats in their pathogenesis.

CONCLUSIONS

These results demonstrate that FLNA is prone to pathogenic rearrangements, and highlight the importance of screening for CNVs in individuals with PNH lacking FLNA point mutations.

Mutagenic Inverted Repeats Assisted Genome Engineering (MIRAGE) in Saccharomyces cerevisiae: deletion of gal7

MIRAGE is a unique in vivo genome editing technique that exploits the inherent instability of inverted repeats (palindromes) in the Saccharomyces cerevisiae chromosome. As a technique able to quickly create deletions as well as precise point mutations, it is valuable in applications that require creation of designer strains of this yeast. In particular, it has various potential applications in metabolic engineering, systems biology, synthetic biology, and molecular genetics.

Real-time detection of cruciform extrusion by single-molecule DNA nanomanipulation

During cruciform extrusion, a DNA inverted repeat unwinds and forms a four-way junction in which two of the branches consist of hairpin structures obtained by self-pairing of the inverted repeats. Here, we use single-molecule DNA nanomanipulation to monitor in real-time cruciform extrusion and rewinding. This allows us to determine the size of the cruciform to nearly base pair accuracy and its kinetics with second-scale time resolution. We present data obtained with two different inverted repeats, one perfect and one imperfect, and extend single-molecule force spectroscopy to measure the torque dependence of cruciform extrusion and rewinding kinetics. Using mutational analysis and a simple two-state model, we find that in the transition state intermediate only the B-DNA located between the inverted repeats (and corresponding to the unpaired apical loop) is unwound, implying that initial stabilization of the four-way (or Holliday) junction is rate-limiting. We thus find that cruciform extrusion is kinetically regulated by features of the hairpin loop, while rewinding is kinetically regulated by features of the stem. These results provide mechanistic insight into cruciform extrusion and help understand the structural features that determine the relative stability of the cruciform and B-form states.

Friedreich's ataxia (GAA)n•(TTC)n repeats strongly stimulate mitotic crossovers in Saccharomyces cerevisae

Expansions of trinucleotide GAA•TTC tracts are associated with the human disease Friedreich's ataxia, and long GAA•TTC tracts elevate genome instability in yeast. We show that tracts of (GAA)₂₃₀•(TTC)₂₃₀ stimulate mitotic crossovers in yeast about 10,000-fold relative to a “normal” DNA sequence; (GAA)_n•(TTC)_n tracts, however, do not significantly elevate meiotic recombination. Most of the mitotic crossovers are associated with a region of non-reciprocal transfer of information (gene conversion). The major class of recombination events stimulated by (GAA)_n•(TTC)_n tracts is a tract-associated double-strand break (DSB) that occurs in unreplicated chromosomes, likely in G1 of the cell cycle. These findings indicate that (GAA)_n•(TTC)_n tracts can be a potent source of loss of heterozygosity in yeast.

Although meiotic recombination has been much more studied than mitotic recombination, mitotic recombination is a universal property. Meiotic recombination rates are quite variable within the genome, with some chromosomal regions (hotspots) having much higher levels of exchange than other regions (coldspots). For mitotic recombination, although some types of DNA sequences are known to be associated with elevated recombination rates (highly-transcribed genes, inverted repeated sequences), relatively few hotspots have been described. In this report, we show that a 690 base pair region consisting of 230 copies of the (GAA)_n•(TTC)_n trinucleotide repeat stimulates mitotic crossovers in yeast 10,000-fold more strongly than an “average” yeast sequence. This sequence is a preferred site for chromosome breakage in stationary phase yeast cells. Our findings may be relevant to understanding the expansions of the (GAA)_n•(TTC)_n trinucleotide repeat tracts that are associated with the human disease Friedreich's ataxia.

The distribution of inverted repeat sequences in the Saccharomyces cerevisiae genome

Although a variety of possible functions have been proposed for inverted repeat sequences (IRs), it is not known which of them might occur in vivo. We investigate this question by assessing the distributions and properties of IRs in the Saccharomyces cerevisiae (SC) genome. Using the IRFinder algorithm we detect 100,514 IRs having copy length greater than 6 bp and spacer length less than 77 bp. To assess statistical significance we also determine the IR distributions in two types of randomization of the S. cerevisiae genome. We find that the S. cerevisiae genome is significantly enriched in IRs relative to random. The S. cerevisiae IRs are significantly longer and contain fewer imperfections than those from the randomized genomes, suggesting that processes to lengthen and/or correct errors in IRs may be operative in vivo. The S. cerevisiae IRs are highly clustered in intergenic regions, while their occurrence in coding sequences is consistent with random. Clustering is stronger in the 3' flanks of genes than in their 5' flanks. However, the S. cerevisiae genome is not enriched in those IRs that would extrude cruciforms, suggesting that this is not a common event. Various explanations for these results are considered.

Second-generation HIF-activated oncolytic adenoviruses with improved replication, oncolytic, and antitumor efficacy

There is a need to develop more potent oncolytic adenoviruses that exhibit increased anti-tumor activity in patients. The HYPR-Ads are targeted oncolytic adenoviruses that specifically kill tumor cells which express active hypoxia-inducible factor (HIF). While therapeutically efficacious, the HYPR-Ads exhibited attenuated replication and oncolytic activity. To overcome these deficiencies and improve anti-tumor efficacy, we created new HIF-activated oncolytic Ads, HIF-Ad and HIF-Ad-IL4, which have two key changes: (i) a modified HIF-responsive promoter to regulate the E1A replication gene and (ii) insertion of the E3 gene region. The HIF-Ads demonstrated conditional activation of E1A expression under hypoxia. Importantly, the HIF-Ads exhibit hypoxia-dependent replication, oncolytic, and cellular release activities and potent anti-tumor efficacy, all of which are significantly greater than the HYPR-Ads. Notably, HIF-Ad-IL4 treatment led to regressions in tumor size by 70% and extensive tumor infiltration by leukocytes resulting in an anti-tumor efficacy that is up to 6-fold greater than the HYPR-Ads, HIF-Ad, and wild-type adenovirus treatment. These studies demonstrate that treatment with a HIF-activated oncolytic adenovirus leads to a measurable therapeutic response. The novel design of the HIF-Ads represents a significant improvement compared to first-generation oncolytic Ads and has great potential to increase the efficacy of this cancer therapy.

Tandem repeat coupled with endonuclease cleavage (TREC): a seamless modification tool for genome engineering in yeast

The complete synthetic Mycoplasma genitalium genome (∼583 kb) has been assembled and cloned as a circular plasmid in the yeast Saccharomyces cerevisiae. Attempts to engineer the cloned genome by standard genetic methods involving the URA3/5-fluoroorotic acid (5-FOA) counter-selection have shown a high background of 5-FOA resistant clones derived from spontaneous deletions of the bacterial genome maintained in yeast. Here, we report a method that can seamlessly modify the bacterial genome in yeast with high efficiency. This method requires two sequential homologous recombination events. First, the target region is replaced with a mutagenesis cassette that consists of a knock-out CORE (an18-bp I-SceI recognition site, the SCEI gene under the control of the GAL1 promoter, and the URA3 marker) and a DNA fragment homologous to the sequence upstream of the target site. The replacement generates tandem repeat sequences flanking the CORE. Second, galactose induces the expression of I-SceI, which generates a double-strand break (DSB) at the recognition site. This DSB promotes intra-molecular homologous recombination between the repeat sequences, and leads to an excision of the CORE. As a result, a seamless modification is generated. This method can be adapted for a variety of genomic modifications and may provide an important tool to modify and design natural or synthetic genomes propagated in yeast.

The transposon Galileo generates natural chromosomal inversions in Drosophila by ectopic recombination

Background

Transposable elements (TEs) are responsible for the generation of chromosomal inversions in several groups of organisms. However, in Drosophila and other Dipterans, where inversions are abundant both as intraspecific polymorphisms and interspecific fixed differences, the evidence for a role of TEs is scarce. Previous work revealed that the transposon Galileo was involved in the generation of two polymorphic inversions of Drosophila buzzatii.

Methodology/Principal Findings

To assess the impact of TEs in Drosophila chromosomal evolution and shed light on the mechanism involved, we isolated and sequenced the two breakpoints of another widespread polymorphic inversion from D. buzzatii, 2z³. In the non inverted chromosome, the 2z³ distal breakpoint was located between genes CG2046 and CG10326 whereas the proximal breakpoint lies between two novel genes that we have named Dlh and Mdp. In the inverted chromosome, the analysis of the breakpoint sequences revealed relatively large insertions (2,870-bp and 4,786-bp long) including two copies of the transposon Galileo (subfamily Newton), one at each breakpoint, plus several other TEs. The two Galileo copies: (i) are inserted in opposite orientation; (ii) present exchanged target site duplications; and (iii) are both chimeric.

Conclusions/Significance

Our observations provide the best evidence gathered so far for the role of TEs in the generation of Drosophila inversions. In addition, they show unequivocally that ectopic recombination is the causative mechanism. The fact that the three polymorphic D. buzzatii inversions investigated so far were generated by the same transposon family is remarkable and is conceivably due to Galileo's unusual structure and current (or recent) transpositional activity.