Inverted repeats as genetic elements for promoting DNA inverted duplication: implications in gene amplification

Inverted repeats in the monkeypox virus genome are hot spots for mutation

The current monkeypox virus (MPXV) strain differs from the strain arising in 2018 by 50+ single nucleotide polymorphisms (SNPs) and is mutating much faster than expected. The cytidine deaminase apolipoprotein B messenger RNA editing enzyme, catalytic subunit B (APOBEC3) was hypothesized to be driving this increased mutation. APOBEC has recently been identified to preferentially mutate cruciform DNA secondary structures formed by inverted repeats (IRs). IRs were recently identified as hot spots for mutation in severe acute respiratory syndrome coronavirus 2, and we aimed to identify whether IRs were also hot spots for mutation within MPXV genomes. We found that MPXV genomes were replete with IR sequences. Of the 50+ SNPs identified in the 2022 outbreak strain, 63.9% of these were found to have arisen within IR regions in the 2018 reference strain (MT903344.1). Notably, IR sequences found in the 2018 reference strain were significantly lost over time, with an average of 32.5% of these sequences being conserved in the 2022 MPXV genomes. This evidence was highly indicative that mutations were arising within IRs. This data provides further support to the hypothesis that APOBEC may be driving MPXV mutation and highlights the necessity for greater surveillance of IRs of MPXV genomes to detect new mutations.

The complexity landscape of viral genomes

BACKGROUND

Viruses are among the shortest yet highly abundant species that harbor minimal instructions to infect cells, adapt, multiply, and exist. However, with the current substantial availability of viral genome sequences, the scientific repertory lacks a complexity landscape that automatically enlights viral genomes' organization, relation, and fundamental characteristics.

RESULTS

This work provides a comprehensive landscape of the viral genome's complexity (or quantity of information), identifying the most redundant and complex groups regarding their genome sequence while providing their distribution and characteristics at a large and local scale. Moreover, we identify and quantify inverted repeats abundance in viral genomes. For this purpose, we measure the sequence complexity of each available viral genome using data compression, demonstrating that adequate data compressors can efficiently quantify the complexity of viral genome sequences, including subsequences better represented by algorithmic sources (e.g., inverted repeats). Using a state-of-the-art genomic compressor on an extensive viral genomes database, we show that double-stranded DNA viruses are, on average, the most redundant viruses while single-stranded DNA viruses are the least. Contrarily, double-stranded RNA viruses show a lower redundancy relative to single-stranded RNA. Furthermore, we extend the ability of data compressors to quantify local complexity (or information content) in viral genomes using complexity profiles, unprecedently providing a direct complexity analysis of human herpesviruses. We also conceive a features-based classification methodology that can accurately distinguish viral genomes at different taxonomic levels without direct comparisons between sequences. This methodology combines data compression with simple measures such as GC-content percentage and sequence length, followed by machine learning classifiers.

CONCLUSIONS

This article presents methodologies and findings that are highly relevant for understanding the patterns of similarity and singularity between viral groups, opening new frontiers for studying viral genomes' organization while depicting the complexity trends and classification components of these genomes at different taxonomic levels. The whole study is supported by an extensive website (https://asilab.github.io/canvas/) for comprehending the viral genome characterization using dynamic and interactive approaches.

Two New M23 Peptidoglycan Hydrolases With Distinct Net Charge

Bacterial peptidoglycan hydrolases play an essential role in cell wall metabolism during bacterial growth, division, and elongation (autolysins) or in the elimination of closely related species from the same ecological niche (bacteriocins). Most studies concerning the peptidoglycan hydrolases present in Gram-positive bacteria have focused on clinically relevant Staphylococcus aureus or the model organism Bacillus subtilis, while knowledge relating to other species remains limited. Here, we report two new peptidoglycan hydrolases from the M23 family of metallopeptidases derived from the same staphylococcal species, Staphylococcus pettenkoferi. They share modular architecture, significant sequence identity (60%), catalytic and binding residue conservation, and similar modes of activation, but differ in gene distribution, putative biological role, and, strikingly, in their isoelectric points (pIs). One of the peptides has a high pI, similar to that reported for all M23 peptidases evaluated to date, whereas the other displays a low pI, a unique feature among M23 peptidases. Consequently, we named them SpM23_B (Staphylococcus pettenkoferi M23 "Basic") and SpM23_A (Staphylococcus pettenkoferi M23 "Acidic"). Using genetic and biochemical approaches, we have characterized these two novel lytic enzymes, both in vitro and in their physiological context. Our study presents a detailed characterization of two novel and clearly distinct peptidoglycan hydrolases to understand their role in bacterial physiology.

Morphology and Mitochondrial Genome of Fischoederius sp. 1 in Thailand

A rumen fluke Fischoederius elongatus is assigned to the type species of genus Fischoederius, family Gastrothylacidae. However, the mitochondrial sequences recently published are thought to be of inconsistent species, suggesting that several morphologically similar but genetically distinct species might be classified as Fischoederius elongatus. Thus, mentions of F. elongatus from South, Southeast, and East Asia might unintentionally refer to different species. The present work describes morphology and a full mitochondrial genome sequence of one of these species. The fluke specimens were collected from 2 infected cattle in Thailand. An interesting finding was the presence of a second tRNA-Asp gene next to a partial ND1 gene. It is suggested that these duplicated sequences are the remnants of non-reciprocal recombination events caused by inverted repeats located between ND2 and ND1 mitochondrial genes.

DNA Functionality with Photoswitchable Hydrazone Cytidine*

A new family of hydrazone modified cytidine phosphoramidite building block was synthesized and incorporated into oligodeoxynucleotides to construct photoswitchable DNA strands. The E-Z isomerization triggered by the irradiation of blue light with a wavelength of 450 nm was investigated and confirmed by ¹ H NMR spectroscopy and HPLC in the contexts of both nucleoside and oligodeoxynucleotide. The light activated Z form isomer of this hydrazone-cytidine with a six-member intramolecular hydrogen bond was found to inhibit DNA synthesis in the primer extension model by using Bst DNA polymerase. In addition, the hydrazone modification caused the misincorporation of dATP together with dGTP into the growing DNA strand with similar selectivity, highlighting a potential G to A mutation. This work provides a novel functional DNA building block and an additional molecular tool that has potential chemical biology and biomedicinal applications to control DNA synthesis and DNA-enzyme interactions using the cell friendly blue light irradiation.

Mitochondrial genomes of the green macroalga Ulva pertusa (Ulvophyceae, Chlorophyta): novel insights into the evolution of mitogenomes in the Ulvophyceae

To further understand the trends in the evolution of mitochondrial genomes (mitogenomes or mtDNAs) in the Ulvophyceae, the mitogenomes of two separate thalli of Ulva pertusa were sequenced. Two U. pertusa mitogenomes (Up1 and Up2) were 69,333 bp and 64,602 bp in length. These mitogenomes shared two ribosomal RNAs (rRNAs), 28 transfer RNAs (tRNAs), 29 protein-coding genes, and 12 open reading frames. The 4.7 kb difference in size was attributed to variation in intron content and tandem repeat regions. A total of six introns were present in the smaller U. pertusa mtDNA (Up2), while the larger mtDNA (Up1) had eight. The larger mtDNA had two additional group II introns in two genes (cox1 and cox2) and tandem duplication mutations in noncoding regions. Our results showed the first case of intraspecific variation in chlorophytan mitogenomes and provided further genomic data for the undersampled Ulvophyceae.

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.

Francisella tularensis Subtype A.II Genomic Plasticity in Comparison with Subtype A.I

Although Francisella tularensis is considered a monomorphic intracellular pathogen, molecular genotyping and virulence studies have demonstrated important differences within the tularensis subspecies (type A). To evaluate genetic variation within type A strains, sequencing and assembly of a new subtype A.II genome was achieved for comparison to other completed F. tularensis type A genomes. In contrast with the F. tularensis A.I strains (SCHU S4, FSC198, NE061598, and TI0902), substantial genomic variation was observed between the newly sequenced F. tularensis A.II strain (WY-00W4114) and the only other publically available A.II strain (WY96-3418). Genome differences between WY-00W4114 and WY96-3418 included three major chromosomal translocations, 1580 indels, and 286 nucleotide substitutions of which 159 were observed in predicted open reading frames and 127 were located in intergenic regions. The majority of WY-00W4114 nucleotide deletions occurred in intergenic regions, whereas most of the insertions and substitutions occurred in predicted genes. Of the nucleotide substitutions, 48 (30%) were synonymous and 111 (70%) were nonsynonymous. WY-00W4114 and WY96-3418 nucleotide polymorphisms were predominantly G/C to A/T allelic mutations, with WY-00W4114 having more A+T enrichment. In addition, the A.II genomes contained a considerably higher number of intact genes and longer repetitive sequences, including transposon remnants than the A.I genomes. Together these findings support the premise that F. tularensis A.II may have a fitness advantage compared to the A.I subtype due to the higher abundance of functional genes and repeated chromosomal sequences. A better understanding of the selective forces driving F. tularensis genetic diversity and plasticity is needed.

Formation of linear amplicons with inverted duplications in Leishmania requires the MRE11 nuclease

Extrachromosomal DNA amplification is frequent in the protozoan parasite Leishmania selected for drug resistance. The extrachromosomal amplified DNA is either circular or linear, and is formed at the level of direct or inverted homologous repeated sequences that abound in the Leishmania genome. The RAD51 recombinase plays an important role in circular amplicons formation, but the mechanism by which linear amplicons are formed is unknown. We hypothesized that the Leishmania infantum DNA repair protein MRE11 is required for linear amplicons following rearrangements at the level of inverted repeats. The purified LiMRE11 protein showed both DNA binding and exonuclease activities. Inactivation of the LiMRE11 gene led to parasites with enhanced sensitivity to DNA damaging agents. The MRE11^−/− parasites had a reduced capacity to form linear amplicons after drug selection, and the reintroduction of an MRE11 allele led to parasites regaining their capacity to generate linear amplicons, but only when MRE11 had an active nuclease activity. These results highlight a novel MRE11-dependent pathway used by Leishmania to amplify portions of its genome to respond to a changing environment.

Extrachromosomal DNA amplification is frequent in the human protozoan parasite Leishmania when challenged with drug or other stressful conditions. DNA amplicons, either circular or linear, are formed by recombination between direct or inverted repeats spread throughout the genome of the parasite. The recombinase RAD51 is involved in the formation of circular amplicons, but the mechanism by which linear amplicons are formed is still unknown in this parasite. Studies in other organisms have provided some evidence that a DNA break is required for linear amplifications, and that the DNA repair protein MRE11 can be involved in this process. In this work, we present our biochemical, cellular and molecular characterization of the Leishmania infantum MRE11 orthologue and provide evidence that this nuclease is involved in the formation of linear amplicons in Leishmania. Our results highlight a novel MRE11-dependent pathway used by Leishmania to amplify portions of its genome to respond to a changing environment.

Large inverted duplications in the human genome form via a fold-back mechanism

Inverted duplications are a common type of copy number variation (CNV) in germline and somatic genomes. Large duplications that include many genes can lead to both neurodevelopmental phenotypes in children and gene amplifications in tumors. There are several models for inverted duplication formation, most of which include a dicentric chromosome intermediate followed by breakage-fusion-bridge (BFB) cycles, but the mechanisms that give rise to the inverted dicentric chromosome in most inverted duplications remain unknown. Here we have combined high-resolution array CGH, custom sequence capture, next-generation sequencing, and long-range PCR to analyze the breakpoints of 50 nonrecurrent inverted duplications in patients with intellectual disability, autism, and congenital anomalies. For half of the rearrangements in our study, we sequenced at least one breakpoint junction. Sequence analysis of breakpoint junctions reveals a normal-copy disomic spacer between inverted and non-inverted copies of the duplication. Further, short inverted sequences are present at the boundary of the disomic spacer and the inverted duplication. These data support a mechanism of inverted duplication formation whereby a chromosome with a double-strand break intrastrand pairs with itself to form a “fold-back” intermediate that, after DNA replication, produces a dicentric inverted chromosome with a disomic spacer corresponding to the site of the fold-back loop. This process can lead to inverted duplications adjacent to terminal deletions, inverted duplications juxtaposed to translocations, and inverted duplication ring chromosomes.

Chromosomes with large inverted duplications and terminal deletions cause neurodevelopmental disorders in children. These chromosome rearrangements typically involve hundreds of genes, leading to significant changes in gene dosage. Though inverted duplications adjacent to terminal deletions are a relatively common type of chromosomal imbalance, the DNA repair mechanism responsible for their formation is not known. In this study, we analyze the genomic organization of the largest collection of human inverted duplications. We find a common inverted duplication structure, consistent with a model that requires DNA to fold back and form a dicentric chromosome intermediate. These data provide insight into the formation of nonrecurrent inverted duplications in the human genome.

Probing mercury(II)-DNA interactions by nanopore stochastic sensing

In this work, DNA-Hg(II) interactions were investigated by monitoring the translocation of DNA hairpins in a protein ion channel in the absence and presence of metal ions. Our experiments demonstrate that target-specific hairpin structures could be stabilized much more significantly by mercuric ions than by the stem length and the loop size of the hairpin due to the formation of Thymine-Hg(II)-Thymine complexes. In addition, the designed DNA probe allows the development of a highly sensitive nanopore sensor for Hg(2+) with a detection limit of 25 nM. Further, the sensor is specific, and other tested metal ions including Pb(2+), Cu(2+), Cd(2+), and so on with concentrations of up to 2 orders of magnitude greater than that of Hg(2+) would not interfere with the mercury detection.

Instability of the octarepeat region of the human prion protein gene

Prion diseases are a family of unique fatal transmissible neurodegenerative diseases that affect humans and many animals. Sporadic Creutzfeldt-Jakob disease (sCJD) is the most common prion disease in humans, accounting for 85–90% of all human prion cases, and exhibits a high degree of diversity in phenotypes. The etiology of sCJD remains to be elucidated. The human prion protein gene has an octapeptide repeat region (octarepeats) that normally contains 5 repeats of 24–27 bp (1 nonapeptide and 4 octapeptide coding sequences). An increase of the octarepeat numbers to six or more or a decrease of the octarepeat number to three is linked to genetic prion diseases with heterogeneous phenotypes in humans. Here we report that the human octarepeat region is prone to either contraction or expansion when subjected to PCR amplification in vitro using Taq or Pwo polymerase and when replicated in wild type E. coli cells. Octarepeat insertion mutants were even less stable, and the mutation rate for the wild type octarepeats was much higher when replicated in DNA mismatch repair-deficient E.coli cells. All observed octarepeat mutants resulting from DNA replication in E.coli were contained in head-to-head plasmid dimers and DNA mfold analysis (http://mfold.rna.albany.edu/?q=mfold/DNA-Folding-Form) indicates that both DNA strands of the octarepeat region would likely form multiple stable hairpin structures, suggesting that the octarepeat sequence may form stable hairpin structures during DNA replication or repair to cause octarepeat instability. These results provide the first evidence supporting a somatic octarepeat mutation-based model for human sCJD etiology: 1) the instability of the octarepeat region leads to accumulation of somatic octarepeat mutations in brain cells during development and aging, 2) this instability is augmented by compromised DNA mismatch repair in aged cells, and 3) eventually some of the octarepeat mutation-containing brain cells start spontaneous de novo prion formation and replication to initiate sCJD.

A framework for automated enrichment of functionally significant inverted repeats in whole genomes

BACKGROUND

RNA transcripts from genomic sequences showing dyad symmetry typically adopt hairpin-like, cloverleaf, or similar structures that act as recognition sites for proteins. Such structures often are the precursors of non-coding RNA (ncRNA) sequences like microRNA (miRNA) and small-interfering RNA (siRNA) that have recently garnered more functional significance than in the past. Genomic DNA contains hundreds of thousands of such inverted repeats (IRs) with varying degrees of symmetry. But by collecting statistically significant information from a known set of ncRNA, we can sort these IRs into those that are likely to be functional.

RESULTS

A novel method was developed to scan genomic DNA for partially symmetric inverted repeats and the resulting set was further refined to match miRNA precursors (pre-miRNA) with respect to their density of symmetry, statistical probability of the symmetry, length of stems in the predicted hairpin secondary structure, and the GC content of the stems. This method was applied on the Arabidopsis thaliana genome and validated against the set of 190 known Arabidopsis pre-miRNA in the miRBase database. A preliminary scan for IRs identified 186 of the known pre-miRNA but with 714700 pre-miRNA candidates. This large number of IRs was further refined to 483908 candidates with 183 pre-miRNA identified and further still to 165371 candidates with 171 pre-miRNA identified (i.e. with 90% of the known pre-miRNA retained).

CONCLUSIONS

165371 candidates for potentially functional miRNA is still too large a set to warrant wet lab analyses, such as northern blotting, on all of them. Hence additional filters are needed to further refine the number of candidates while still retaining most of the known miRNA. These include detection of promoters and terminators, homology analyses, location of candidate relative to coding regions, and better secondary structure prediction algorithms. The software developed is designed to easily accommodate such additional filters with a minimal experience in Perl.

E. coli SbcCD and RecA control chromosomal rearrangement induced by an interrupted palindrome

Survival and genome stability are critical characteristics of healthy cells. DNA palindromes pose a threat to genome stability and have been shown to participate in a reaction leading to the formation of inverted chromosome duplications centered around themselves. There is considerable interest in the mechanism of this rearrangement given its likely contribution to genome instability in cancer cells. This study shows that formation of large inverted chromosome duplications can be observed in the chromosome of Escherichia coli. They are formed at the site of a 246 bp interrupted DNA palindrome in the absence of the hairpin nuclease SbcCD and the recombination protein RecA. The genetic requirements for this spontaneous rearrangement are consistent with a pathway involving DNA degradation and hairpin formation, as opposed to a cruciform cleavage pathway. Accordingly, the formation of palindrome-dependent hairpin intermediates can be induced by an adjacent DNA double-stand break.

Discovery of a long inverted repeat in human POTE genes

POTE gene family is tightly related to prostate, ovary, testis and placenta cancers. We recently identified an intronic long inverted repeat (LIR) in some members of the POTE gene family. Due to the capacity of inducing gene amplification, the POTE intronic LIRs may be involved in over-expression of the POTE genes. Our study aimed to understand the origin of the LIR in primates. We collected the LIR and its flanking sequences within rhesus monkey, chimpanzee and human genomes. The rhesus monkey genome only has half-sized LIRs (lack one repeat copy), whereas the human and chimpanzee genomes contain both full-sized and half-sized LIRs. Phylogenetic tree indicates that the LIR is formed after divergence of rhesus monkey and the common ancestor of human and chimpanzee. The POTE genes containing a full-sized LIR were amplified in the human genome.

Acquisition of inverted GSTM exons by an intron of primate GSTM5 gene

The human GSTM gene family is composed of five gene members, GSTM1-5, and plays an important role in detoxification. In this study, the human GSTM5 gene was found to have a long inverted repeat (LIR) in intron 5. The LIR is able to form a stem-loop structure with a 31-bp stem and a 9-nt loop. The intronic LIR was also identified in other primates but not in non-primates. The human and chimpanzee LIRs had undergone compensating mutations that make the stem loop more stable, suggesting a functional role for the LIR. Sequence homology showed that the LIR was actually a part of inverted exons acquired by the intron. Results of phylogenetic analysis indicate that the inverted exons were derived from exon 5 of GSTM4 and exon 5 of GSTM1. The intronic LIR and inverted GSTM exons can probably introduce complexity in the expression of GSTM gene family.

A study on genomic distribution and sequence features of human long inverted repeats reveals species-specific intronic inverted repeats

The inverted repeats present in a genome play dual roles. They can induce genomic instability and, on the other hand, regulate gene expression. In the present study, we report the distribution and sequence features of recombinogenic long inverted repeats (LIRs) that are capable of forming stable stem-loops or palindromes within the human genome. A total of 2551 LIRs were identified, and 37% of them were located in long introns (largely > 10 kb) of genes. Their distribution appears to be random in introns and is not restrictive, even for regions near intron-exon boundaries. Almost half of them comprise TG/CA-rich repeats, inversely arranged Alu repeats and MADE1 mariners. The remaining LIRs are mostly unique in their sequence features. Comparative studies of human, chimpanzee, rhesus monkey and mouse orthologous genes reveal that human genes have more recombinogenic LIRs than other orthologs, and over 80% are human-specific. The human genes associated with the human-specific LIRs are involved in the pathways of cell communication, development and the nervous system, as based on significantly over-represented Gene Ontology terms. The functional pathways related to the development and functions of the nervous system are not enriched in chimpanzee and mouse orthologs. The findings of the present study provide insight into the role of intronic LIRs in gene regulation and primate speciation.

The presence of distal and proximal promoters for rat mitochondrial glycerol-3-phosphate acyltransferase

Sequence analysis using the Promoser program predicted two promoter-like regions for rat mtGPAT: a distal promoter approximately 30kb upstream and a proximal promoter near the first translational codon. Rat liver cells transfected with pGL3-basic vector containing the distal and proximal promoter resulted in 10.8- and 4.8-fold increase in the luciferase activity, respectively. Results of electromobility shift assay and chromatin immunoprecipitation suggested binding of transcription factors to the distal and proximal promoter regions. 5' RACE PCR showed two transcripts with different transcriptional start sites. When transfected rat liver cells were starved and refed, there was about 2.7-fold increase in the luciferase activity with cells transfected with the distal promoter while the proximal promoter showed no change. Thus, the two promoters could be functionally distinguished. Taken together, the results suggest that there are two promoters for rat mtGPAT gene and that the transcriptional regulation is mediated through the distal promoter.

Double-strand breaks in the myotonic dystrophy type 1 and the fragile X syndrome triplet repeat sequences induce different types of mutations in DNA flanking sequences in Escherichia coli

The putative role of double-strand breaks (DSBs) created in vitro by restriction enzyme cleavage in or near CGG*CCG or CTG*CAG repeat tracts on their genetic instabilities, both within the repeats and in their flanking sequences, was investigated in an Escherichia coli plasmid system. DSBs at TRS junctions with the vector generated a large number of mutagenic events in flanking sequences whereas DSBs within the repeats elicited no similar products. A substantial enhancement in the number of mutants was caused by transcription of the repeats and by the absence of recombination functions (recA-, recBC-). Surprisingly, DNA sequence analyses on mutant clones revealed the presence of only single deletions of 0.4-1.6 kb including the TRS and the flanking sequence from plasmids originally containing (CGG*CCG)43 but single, double and multiple deletions as well as insertions were found for plasmids originally containing (CTG*CAG)n (where n = 43 or 70). Non-B DNA structures (slipped structures with loops, cruciforms, triplexes and tetraplexes) as well as microhomologies are postulated to participate in the recombination and/or repair processes.

On the Paucity of Duplicated Genes in Caenorhabditis elegans Operons

Spliced leader trans-splicing is an mRNA maturation process used by a small set of eukaryotes, including the nematode C. elegans, to cap the downstream genes of operons. We analyzed the frequency of duplication of operonic genes in C. elegans and confirmed that they are duplicated less often in the genome than monocistronic genes. Because operons account for about 15% of the genes in C. elegans, this lower duplication frequency might place a large constraint on the plasticity of the genome. Further analyses suggest that this paucity of duplicated genes results from operon organization hindering specific types of gene duplication.

Long inverted repeats in eukaryotic genomes: recombinogenic motifs determine genomic plasticity

Inverted repeats are unstable motifs in a genome, having a causal relation to fragment rearrangements and recombination events. We have investigated long inverted repeats (LIR) of > 30 bp in length in eukaryotic genomes to assess their contribution to genome stability. An algorithm was first designed for searching for LIRs with < 2 kb internal spacers and >85% identity (degree of homology between repeat copies of a LIR). There are much fewer LIRs in yeast, fruitfly, pufferfish and chicken than in Caenorhabditis elegans, zebrafish, frog and human. However, the high LIR frequencies do not necessarily imply high genome instability because of variant internal spacers and stem lengths and identities. From the collection of identified LIRs, we selected recombinogenic LIRs that had a short internal spacer and a high copy identity and were prone to induce high instability. We found that a relatively high proportion (5-9.8%) of the LIRs in C. elegans, zebrafish and frog were recombinogenic LIRs. In contrast, the proportions in human and mouse LIRs were quite low (0.4-1.1%) basically accounting for long internal spacers. We suggest that C. elegans, zebrafish and frog genomes are unstable in terms of the LIR frequency and the proportion of recombinogenic LIRs. For the other genomes, LIRs most likely have a minor impact.

A mechanism of palindromic gene amplification in Saccharomyces cerevisiae

Selective gene amplification is associated with normal development, neoplasia, and drug resistance. One class of amplification events results in large arrays of inverted repeats that are often complex in structure, thus providing little information about their genesis. We made a recombination substrate in Saccharomyces cerevisiae that frequently generates palindromic duplications to repair a site-specific double-strand break in strains deleted for the SAE2 gene. The resulting palindromes are stable in sae2Delta cells, but unstable in wild-type cells. We previously proposed that the palindromes are formed by invasion and break-induced replication, followed by an unknown end joining mechanism. Here we demonstrate that palindrome formation can occur in the absence of RAD50, YKU70, and LIG4, indicating that palindrome formation defines a new class of nonhomologous end joining events. Sequence data from 24 independent palindromic duplication junctions suggest that the duplication mechanism utilizes extremely short (4-6 bp), closely spaced (2-9 bp), inverted repeats to prime DNA synthesis via an intramolecular foldback of a 3' end. In view of our data, we present a foldback priming model for how a single copy sequence is duplicated to generate a palindrome.

The Structure and Early Evolution of Recently Arisen Gene Duplicates in theCaenorhabditis elegansGenome

The significance of gene duplication in provisioning raw materials for the evolution of genomic diversity is widely recognized, but the early evolutionary dynamics of duplicate genes remain obscure. To elucidate the structural characteristics of newly arisen gene duplicates at infancy and their subsequent evolutionary properties, we analyzed gene pairs with ≤10% divergence at synonymous sites within the genome of Caenorhabditis elegans. Structural heterogeneity between duplicate copies is present very early in their evolutionary history and is maintained over longer evolutionary timescales, suggesting that duplications across gene boundaries in conjunction with shuffling events have at least as much potential to contribute to long-term evolution as do fully redundant (complete) duplicates. The median duplication span of 1.4 kb falls short of the average gene length in C. elegans (2.5 kb), suggesting that partial gene duplications are frequent. Most gene duplicates reside close to the parent copy at inception, often as tandem inverted loci, and appear to disperse in the genome as they age, as a result of reduced survivorship of duplicates located in proximity to the ancestral copy. We propose that illegitimate recombination events leading to inverted duplications play a disproportionately large role in gene duplication within this genome in comparison with other mechanisms.

Diverse sequences within Tlr elements target programmed DNA elimination in Tetrahymena thermophila

Tlr elements are a novel family of approximately 30 putative mobile genetic elements that are confined to the germ line micronuclear genome in Tetrahymena thermophila. Thousands of diverse germ line-limited sequences, including the Tlr elements, are specifically eliminated from the differentiating somatic macronucleus. Macronucleus-retained sequences flanking deleted regions are known to contain cis-acting signals that delineate elimination boundaries. It is unclear whether sequences within deleted DNA also play a regulatory role in the elimination process. In the current study, an in vivo DNA rearrangement assay was used to identify internal sequences required in cis for the elimination of Tlr elements. Multiple, nonoverlapping regions from the approximately 23-kb Tlr elements were independently sufficient to stimulate developmentally regulated DNA elimination when placed within the context of flanking sequences from the most thoroughly characterized family member, Tlr1. Replacement of element DNA with macronuclear or foreign DNA abolished elimination activity. Thus, diverse sequences dispersed throughout Tlr DNA contain cis-acting signals that target these elements for programmed elimination. Surprisingly, Tlr DNA was also efficiently deleted when Tlr1 flanking sequences were replaced with DNA from a region of the genome that is not normally associated with rearrangement, suggesting that specific flanking sequences are not required for the elimination of Tlr element DNA.

The Mre11 complex is required for repair of hairpin-capped double-strand breaks and prevention of chromosome rearrangements

Inverted repeats (IRs) that can form a hairpin or cruciform structure are common in the human genome and may be sources of instability. An IR involving the human Alu sequence (Alu-IR) has been studied as a model of such structures in yeast. We found that an Alu-IR is a mitotic recombination hotspot requiring MRE11/RAD50/XRS2 and SAE2. Using a newly developed approach for mapping rare double-strand breaks (DSBs), we established that induction of recombination results from breaks that are terminated by hairpins. Failure of the mre11, rad50, xrs2, and sae2 mutants to process the hairpins blocks recombinational repair of the DSBs and leads to generation of chromosome inverted duplications. Our results suggest an additional role for the Mre11 complex in maintaining genome stability.

Suppression of gene amplification and chromosomal DNA integration by the DNA mismatch repair system

Mismatch repair (MMR)-deficient cells are shown to produce >15-fold more methotrexate-resistant colonies than MMR normal cells. The increased resistance to methotrexate is primarily due to gene amplification since all the resistant clones contain double-minute chromosomes and increased copy numbers of the DHFR gene. In addition, integration of linearized or retroviral DNAs into chromosomes is also significantly elevated in MMR-deficient cells. These results suggest that in addition to microsatellite instability and homeologous recombination, MMR is also involved in suppression of other genome instabilities such as gene amplification and chromosomal DNA integration.