Palindromes and genomic stress fractures: bracing and repairing the damage

A reference catalog of DNA palindromes in the human genome and their variations in 1000 Genomes

A palindrome in DNA is like a palindrome in language, but when read backwards, it is a complement of the forward sequence; effectively, the two halves of a sequence complement each other from its midpoint like in a double strand of DNA. Palindromes are distributed throughout the human genome and play significant roles in gene expression and regulation. Palindromic mutations are linked to many human diseases, such as neuronal disorders, mental retardation, and various cancers. In this work, we computed and analyzed the palindromic sequences in the human genome and studied their conservation in personal genomes using 1000 Genomes data. We found that ~30% of the palindromes exhibit variation, some of which are caused by rare variants. The analysis of disease/trait-associated single-nucleotide polymorphisms in palindromic regions showed that disease-associated risk variants are 14 times more likely to be present in palindromic regions than in other regions. The catalog of palindromes in the reference genome and 1000 Genomes is being made available here with details on their variations in each individual genome to serve as a resource for future and retrospective whole-genome studies identifying statistically significant palindrome variations associated with diseases or traits and their roles in disease mechanisms.

A curious new role for MRN in Schizosaccharomyces pombe non-homologous end-joining

Chromosomal breaks can be healed by several repair processes, including one called non-homologous end-joining (NHEJ) where the two broken ends are ligated together with a loss of 0-5 bp of DNA. The protein requirements for NHEJ of cut DNA ends in the budding yeast Saccharomyces cerevisiae include its version of the Mre11-Rad50-Nbs1 (MRN) complex. In contrast, the fission yeast Schizosaccharomyces pombe and mammalian cells do not require MRN for this process. Recent work in S. pombe used transposon excision to generate breaks that were capped by DNA hairpins, which must be opened to produce ligatable ends. Repair in S. pombe was through an NHEJ reaction that now requires MRN. Surprisingly, wild type cells and MRN mutants that lack nuclease activity showed the same levels of excision. These genetic results suggest that MRN recruits an unknown hairpin-opening nuclease for this unusual NHEJ reaction.

Nonhomologous End-Joining with Minimal Sequence Loss Is Promoted by the Mre11-Rad50-Nbs1-Ctp1 Complex in Schizosaccharomyces pombe

While the Mre11-Rad50-Nbs1 (MRN) complex has known roles in repair processes like homologous recombination and microhomology-mediated end-joining, its role in nonhomologous end-joining (NHEJ) is unclear as Saccharomyces cerevisiae, Schizosaccharomyces pombe, and mammals have different requirements for repairing cut DNA ends. Most double-strand breaks (DSBs) require nucleolytic processing prior to DNA ligation. Therefore, we studied repair using the Hermes transposon, whose excision leaves a DSB capped by hairpin ends similar to structures generated by palindromes and trinucleotide repeats. We generated single Hermes insertions using a novel S. pombe transient transfection system, and used Hermes excision to show a requirement for MRN in the NHEJ of nonligatable ends. NHEJ repair was indicated by the >1000-fold decrease in excision in cells lacking Ku or DNA ligase 4. Most repaired excision sites had <5 bp of sequence loss or mutation, characteristic for NHEJ and similar excision events in metazoans, and in contrast to the more extensive loss seen in S. cerevisiaeS. pombe NHEJ was reduced >1000-fold in cells lacking each MRN subunit, and loss of MRN-associated Ctp1 caused a 30-fold reduction. An Mre11 dimer is thought to hold DNA ends together for repair, and Mre11 dimerization domain mutations reduced repair 300-fold. In contrast, a mre11 mutant defective in endonucleolytic activity, the same mutant lacking Ctp1, or the triple mutant also lacking the putative hairpin nuclease Pso2 showed wild-type levels of repair. Thus, MRN may act to recruit the hairpin opening activity that allows subsequent repair.

Translocation and deletion breakpoints in cancer genomes are associated with potential non-B DNA-forming sequences

Gross chromosomal rearrangements (including translocations, deletions, insertions and duplications) are a hallmark of cancer genomes and often create oncogenic fusion genes. An obligate step in the generation of such gross rearrangements is the formation of DNA double-strand breaks (DSBs). Since the genomic distribution of rearrangement breakpoints is non-random, intrinsic cellular factors may predispose certain genomic regions to breakage. Notably, certain DNA sequences with the potential to fold into secondary structures [potential non-B DNA structures (PONDS); e.g. triplexes, quadruplexes, hairpin/cruciforms, Z-DNA and single-stranded looped-out structures with implications in DNA replication and transcription] can stimulate the formation of DNA DSBs. Here, we tested the postulate that these DNA sequences might be found at, or in close proximity to, rearrangement breakpoints. By analyzing the distribution of PONDS-forming sequences within ±500 bases of 19 947 translocation and 46 365 sequence-characterized deletion breakpoints in cancer genomes, we find significant association between PONDS-forming repeats and cancer breakpoints. Specifically, (AT)_n, (GAA)_n and (GAAA)_n constitute the most frequent repeats at translocation breakpoints, whereas A-tracts occur preferentially at deletion breakpoints. Translocation breakpoints near PONDS-forming repeats also recur in different individuals and patient tumor samples. Hence, PONDS-forming sequences represent an intrinsic risk factor for genomic rearrangements in cancer genomes.

Short Inverted Repeats Are Hotspots for Genetic Instability: Relevance to Cancer Genomes

Analyses of chromosomal aberrations in human genetic disorders have revealed that inverted repeat sequences (IRs) often co-localize with endogenous chromosomal instability and breakage hotspots. Approximately 80% of all IRs in the human genome are short (<100 bp), yet the mutagenic potential of such short cruciform-forming sequences has not been characterized. Here, we find that short IRs are enriched at translocation breakpoints in human cancer and stimulate the formation of DNA double-strand breaks (DSBs) and deletions in mammalian and yeast cells. We provide evidence for replication-related mechanisms of IR-induced genetic instability and a novel XPF cleavage-based mechanism independent of DNA replication. These discoveries implicate short IRs as endogenous sources of DNA breakage involved in disease etiology and suggest that these repeats represent a feature of genome plasticity that may contribute to the evolution of the human genome by providing a means for diversity within the population.

GAP-Seq: a method for identification of DNA palindromes

BACKGROUND

Closely spaced long inverted repeats, also known as DNA palindromes, can undergo intrastrand annealing to form DNA hairpins. The ability to form these hairpins results in genome instability, difficulties in maintaining clones in Escherichia coli and major problems for most DNA sequencing approaches. Because of their role in genomic instability and gene amplification in some human cancers, it is important to develop systematic approaches to detect and characterize DNA palindromes.

RESULTS

We developed a new protocol to identify palindromes that couples the S1 nuclease treated Cot0 DNA (GAPF) with high-throughput sequencing (GAP-Seq). Unlike earlier protocols, it does not involve restriction enzymatic digestion prior to DNA snap-back thereby preserving longer DNA sequences. It also indicates the location of the novel junction, which can then be recovered. Using MCF-7 breast cancer cell line as the proof-of-principle analysis, we have identified 35 palindrome candidates and physically characterized the top 5 candidates and their junctions. Because this protocol eliminates many of the false positives that plague earlier techniques, we have improved palindrome identification.

CONCLUSIONS

The GAP-Seq approach underscores the importance of developing new tools for identifying and characterizing palindromes, and provides a new strategy to systematically assess palindromes in genomes. It will be useful for studying human cancers and other diseases associated with palindromes.

Palindrome-mediated and replication-dependent pathogenic structural rearrangements within the NF1 gene

Palindromic sequences can form hairpin structures or cruciform extrusions, which render them susceptible to genomic rearrangements. A 197-bp long palindromic AT-rich repeat (PATRR17) is located within intron 40 of the neurofibromatosis type 1 (NF1) gene (17q11.2). Through comprehensive NF1 analysis, we identified six unrelated patients with a rearrangement involving intron 40 (five deletions and one reciprocal translocation t(14;17)(q32;q11.2)). We hypothesized that PATRR17 may be involved in these rearrangements thereby causing NF1. Breakpoint cloning revealed that PATRR17 was indeed involved in all of the rearrangements. As microhomology was present at all breakpoint junctions of the deletions identified, and PATRR17 partner breakpoints were located within 7.1 kb upstream of PATRR17, fork stalling and template switching/microhomology-mediated break-induced replication was the most likely rearrangement mechanism. For the reciprocal translocation case, a 51 bp insertion at the translocation breakpoints mapped to a short sequence within PATRR17, proximal to the breakpoint, suggesting a multiple stalling and rereplication process, in contrast to previous studies indicating a purely replication-independent mechanism for PATRR-mediated translocations. In conclusion, we show evidence that PATRR17 is a hotspot for pathogenic intragenic deletions within the NF1 gene and suggest a novel replication-dependent mechanism for PATRR-mediated translocation.

Analysis of the t(3;8) of hereditary renal cell carcinoma: a palindrome-mediated translocation

It has emerged that palindrome-mediated genomic instability generates DNA-based rearrangements. The presence of palindromic AT-rich repeats (PATRRs) at the translocation breakpoints suggested a palindrome-mediated mechanism in the generation of several recurrent constitutional rearrangements: the t(11;22), t(17;22), and t(8;22). To date, all reported PATRR-mediated translocations include the PATRR on chromosome 22 (PATRR22) as a translocation partner. Here, the constitutional rearrangement, t(3;8)(p14.2;q24.1), segregating with renal cell carcinoma in two families, is examined. The chromosome 8 breakpoint lies in PATRR8 in the first intron of the RNF139 (TRC8) gene, whereas the chromosome 3 breakpoint is located in an AT-rich palindromic sequence in intron 3 of the FHIT gene (PATRR3). Thus, the t(3;8) is the first PATRR-mediated, recurrent, constitutional translocation that does not involve PATRR22. Furthermore, we detect de novo translocations similar to the t(11;22) and t(8;22), involving PATRR3 in normal sperm. The breakpoint on chromosome 3 is in proximity to FRA3B, the most common fragile site in the human genome and a site of frequent deletions in tumor cells. However, the lack of involvement of PATRR3 sequence in numerous FRA3B-related deletions suggests that there are several different DNA sequence-based etiologies responsible for chromosome 3p14.2 genomic rearrangements.

Sak1 kinase interacts with Pso2 nuclease in response to DNA damage induced by interstrand crosslink-inducing agents in Saccharomyces cerevisiae

By isolating putative binding partners through the two-hybrid system (THS) we further extended the characterization of the specific interstrand cross-link (ICL) repair gene PSO2 of Saccharomyces cerevisiae. Nine fusion protein products were isolated for Pso2p using THS, among them the Sak1 kinase, which interacted with the C-terminal β-CASP domain of Pso2p. Comparison of mutagen-sensitivity phenotypes of pso2Δ, sak1Δ and pso2Δsak1Δ disruptants revealed that SAK1 is necessary for complete WT-like repair. The epistatic interaction of both mutant alleles suggests that Sak1p and Pso2p act in the same pathway of controlling sensitivity to DNA-damaging agents. We also observed that Pso2p is phosphorylated by Sak1 kinase in vitro and co-immunoprecipitates with Sak1p after 8-MOP+UVA treatment. Survival data after treatment of pso2Δ, yku70Δ and yku70Δpso2Δ with nitrogen mustard, PSO2 and SAK1 with YKU70 or DNL4 single-, double- and triple mutants with 8-MOP+UVA indicated that ICL repair is independent of YKu70p and DNL4p in S. cerevisiae. Furthermore, a non-epistatic interaction was observed between MRE11, PSO2 and SAK1 genes after ICL induction, indicating that their encoded proteins act on the same substrate, but in distinct repair pathways. In contrast, an epistatic interaction was observed for PSO2 and RAD52, PSO2 and RAD50, PSO2 and XRS2 genes in 8-MOP+UVA treated exponentially growing cells.

De novo-generated small palindromes are characteristic of amplicon boundary junction of double minutes

Double minutes (DMs) are hallmarks of gene amplification. However, their molecular structure and the mechanisms of formation are largely unknown. To elucidate the structure and underlying molecular mechanism of DMs, we obtained and cloned DMs using microdissection; and degenerated oligonucleotide primed polymerase chain reaction (DOP-PCR) from the ovarian cancer cell line UACC-1598. Two large amplicons, the 284 kb AmpMYCN, originating from locus 2p24.3 and the 391 kb AmpEIF5A2, from locus 3q26.2, were found co-amplified on the same DMs. The two amplicons are joined through a complex 7 kb junction DNA sequence. Analysis of the junction has revealed three de novo created small palindromes surrounding the six breakpoints. Consistent with these observations, we further found that 70% of the 57 reported DM junction sequences have de novo creation of small palindromic sequences surrounding the breakpoints. Together, our findings indicate that de novo-generated small palindromic sequences are characteristic of amplicon boundary junctions on DMs. It is possible that the de novo-generated small palindromic sequences, which may be generated through non-homologous end joining in concert with a novel DNA repair machinery, play a common role in amplicon rejoining and gene amplification.

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.

Nonspaced inverted DNA repeats are preferential targets for homology-directed gene repair in mammalian cells

DNA repeats constitute potential sites for the nucleation of secondary structures such as hairpins and cruciforms. Studies performed mostly in bacteria and yeast showed that these noncanonical DNA structures are breakage-prone, making them candidate targets for cellular DNA repair pathways. Possible culprits for fragility at repetitive DNA sequences include replication and transcription as well as the action of structure–specific nucleases. Despite their patent biological relevance, the parameters governing DNA repeat-associated chromosomal transactions remain ill-defined. Here, we established an episomal recombination system based on donor and acceptor complementary DNA templates to investigate the role of direct and inverted DNA repeats in homologous recombination (HR) in mammalian cells. This system allowed us also to ascertain in a stringent manner the impact of repetitive sequence replication on homology-directed gene repair. We found that nonspaced DNA repeats can, per se, engage the HR pathway of the cell and that this process is primarily dependent on their spacing and relative arrangement (i.e. parallel or antiparallel) rather than on their sequence. Indeed, our data demonstrate that contrary to direct and spaced inverted repeats, nonspaced inverted repeats are intrinsically recombinogenic motifs in mammalian cells lending experimental support to their role in genome dynamics in higher eukaryotes.

Evolution of linear chromosomes and multipartite genomes in yeast mitochondria

Mitochondrial genome diversity in closely related species provides an excellent platform for investigation of chromosome architecture and its evolution by means of comparative genomics. In this study, we determined the complete mitochondrial DNA sequences of eight Candida species and analyzed their molecular architectures. Our survey revealed a puzzling variability of genome architecture, including circular- and linear-mapping and multipartite linear forms. We propose that the arrangement of large inverted repeats identified in these genomes plays a crucial role in alterations of their molecular architectures. In specific arrangements, the inverted repeats appear to function as resolution elements, allowing genome conversion among different topologies, eventually leading to genome fragmentation into multiple linear DNA molecules. We suggest that molecular transactions generating linear mitochondrial DNA molecules with defined telomeric structures may parallel the evolutionary emergence of linear chromosomes and multipartite genomes in general and may provide clues for the origin of telomeres and pathways implicated in their maintenance.

A palindrome-mediated recurrent translocation with 3:1 meiotic nondisjunction: the t(8;22)(q24.13;q11.21)

Palindrome-mediated genomic instability has been associated with chromosomal translocations, including the recurrent t(11;22)(q23;q11). We report a syndrome characterized by extremity anomalies, mild dysmorphia, and intellectual impairment caused by 3:1 meiotic segregation of a previously unrecognized recurrent palindrome-mediated rearrangement, the t(8;22)(q24.13;q11.21). There are at least ten prior reports of this translocation, and nearly identical PATRR8 and PATRR22 breakpoints were validated in several of these published cases. PCR analysis of sperm DNA from healthy males indicates that the t(8;22) arises de novo during gametogenesis in some, but not all, individuals. Furthermore, demonstration that de novo PATRR8-to-PATRR11 translocations occur in sperm suggests that palindrome-mediated translocation is a universal mechanism producing chromosomal rearrangements.

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.

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.

Nearby inverted repeats fuse to generate acentric and dicentric palindromic chromosomes by a replication template exchange mechanism

Gene amplification plays important roles in the progression of cancer and contributes to acquired drug resistance during treatment. Amplification can initiate via dicentric palindromic chromosome production and subsequent breakage-fusion-bridge cycles. Here we show that, in fission yeast, acentric and dicentric palindromic chromosomes form by homologous recombination protein-dependent fusion of nearby inverted repeats, and that these fusions occur frequently when replication forks arrest within the inverted repeats. Genetic and molecular analyses suggest that these acentric and dicentric palindromic chromosomes arise not by previously described mechanisms, but by a replication template exchange mechanism that does not involve a DNA double-strand break. We thus propose an alternative mechanism for the generation of palindromic chromosomes dependent on replication fork arrest at closely spaced inverted repeats.

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.

Chromosome aberrations resulting from double-strand DNA breaks at a naturally occurring yeast fragile site composed of inverted ty elements are independent of Mre11p and Sae2p

Genetic instability at palindromes and spaced inverted repeats (IRs) leads to chromosome rearrangements. Perfect palindromes and IRs with short spacers can extrude as cruciforms or fold into hairpins on the lagging strand during replication. Cruciform resolution produces double-strand breaks (DSBs) with hairpin-capped ends, and Mre11p and Sae2p are required to cleave the hairpin tips to facilitate homologous recombination. Fragile site 2 (FS2) is a naturally occurring IR in Saccharomyces cerevisiae composed of a pair of Ty1 elements separated by approximately 280 bp. Our results suggest that FS2 forms a hairpin, rather than a cruciform, during replication in cells with low levels of DNA polymerase. Cleavage of this hairpin results in a recombinogenic DSB. We show that DSB formation at FS2 does not require Mre11p, Sae2p, Rad1p, Slx4p, Pso2p, Exo1p, Mus81p, Yen1p, or Rad27p. Also, repair of DSBs by homologous recombination is efficient in mre11 and sae2 mutants. Homologous recombination is impaired at FS2 in rad52 mutants and most aberrations reflect either joining of two broken chromosomes in a "half crossover" or telomere capping of the break. In support of hairpin formation precipitating DSBs at FS2, two telomere-capped deletions had a breakpoint near the center of the IR. In summary, Mre11p and Sae2p are not required for DSB formation at FS2 or the subsequent repair of these DSBs.

Human SLX4 is a Holliday junction resolvase subunit that binds multiple DNA repair/recombination endonucleases

Structure-specific endonucleases resolve DNA secondary structures generated during DNA repair and recombination. The yeast 5' flap endonuclease Slx1-Slx4 has received particular attention with the finding that Slx4 has Slx1-independent key functions in genome maintenance. Although Slx1 is a highly conserved protein in eukaryotes, no orthologs of Slx4 were reported other than in fungi. Here we report the identification of Slx4 orthologs in metazoa, including fly MUS312, essential for meiotic recombination, and human BTBD12, an ATM/ATR checkpoint kinase substrate. Human SLX1-SLX4 displays robust Holliday junction resolvase activity in addition to 5' flap endonuclease activity. Depletion of SLX1 and SLX4 results in 53BP1 foci accumulation and H2AX phosphorylation as well as cellular hypersensitivity to MMS. Furthermore, we show that SLX4 binds the XPF(ERCC4) and MUS81 subunits of the XPF-ERCC1 and MUS81-EME1 endonucleases and is required for DNA interstrand crosslink repair. We propose that SLX4 acts as a docking platform for multiple structure-specific endonucleases.

Linkage disequilibrium between two high-frequency deletion polymorphisms: implications for association studies involving the glutathione-S transferase (GST) genes

Copy number variations (CNVs) represent a large source of genetic variation in humans and have been increasingly studied for disease association. A deletion polymorphism of the gene encoding the cytosolic detoxification enzyme glutathione S-transferase theta 1 (GSTT1) has been extensively studied for cancer susceptibility (919 studies, from HuGE navigator, http://www.hugenavigator.net/). However, clear conclusions have not been reached. Since the GSTT1 gene is located within a genomic region of segmental duplications (SD), there may be a confounding effect from another, yet-uncharacterized CNV at the same locus. Here we describe a previously uncharacterized 38-kilo-base (kb) long deletion polymorphism of GSTT2B located within a 61-kb DNA inverted repeat. GSTT2B is a duplicated copy of GSTT2, the only paralogue of GSTT1 in humans. A newly developed PCR assay revealed that a microhomology-mediated breakpoint appears to be shared among individuals at high frequency. The GSTT2B deletion polymorphism was in strong linkage disequilibrium (LD) (D′ = 0.841) with the neighboring GSTT1 deletion polymorphism in the Caucasian population. Alleles harboring a single deletion were significantly overrepresented (p = 2.22×10⁻¹⁶), suggesting a selection against alleles with both deletions. The deletion alleles are almost certainly the derived ones, because the GSTT2B-GSTT2-GSTT1 genes were strictly retained in chimpanzees. Extremely low GSTT2 mRNA expression was associated with the GSTT2B deletion, suggesting an influence of the deletion on the flanking region and loss of GSTT2 function. Genome-wide LD analysis between deletion polymorphisms further points to the uniqueness of two deletions, because strong LD between deletion polymorphisms might be very rare in humans. These results show a complex genomic organization and unexpected biological functions of CNVs within segmental duplications and emphasize the importance of detailed structural characterization for disease association studies.

Common diseases such as cancer are caused by interactions between multiple genetic and environmental factors. Glutathione S-transferases (GST) are key enzymes in eliminating carcinogens and harmful macromolecules from cells. Based on the assumption that individuals who do not have a particular type of GST genes are susceptible to cancers, a number of studies have been conducted to find a link between GST genotypes and cancer. However such associations remain inconclusive to date. Because GST genes are clustered in repetitive, complex regions in the genome, other previously uncharacterized variations/polymorphisms may have had an impact on the data. We describe here such a genotype, a 37-kb deletion of GSTT2B gene that is found very frequently among humans. The neighboring GSTT2 gene expression is greatly impaired by the GSTT2B deletion, conferring a potentially null allele at GSTT2. The GSTT2B deletion is non-randomly associated with another high frequency deletion of the GSTT1 gene. Therefore, a detailed characterization of this complex region of the genome revealed unexpected genetic and biological interactions of large deletion polymorphisms; this is essential to consider in future disease association studies.

Checkpoint responses to unusual structures formed by DNA repeats

DNA sequences that are prone to adopting non-B DNA secondary structures are associated with hotspots of genomic instability. The fine mechanisms by which alternative DNA structures induce phenomena such as repeat expansions, chromosomal fragility, or gross chromosomal rearrangements are under intensive studies. It is well established that DNA damage checkpoint responses play a crucial role in maintaining a stable genome. It is far less clear, however, whether and how the checkpoint machinery responds to alternative DNA structures. This review discusses the role of the interplay between DNA damage checkpoints and alternative DNA structures in genome maintenance.

Size-dependent palindrome-induced intrachromosomal recombination in yeast

Palindromic and quasi-palindromic sequences are important DNA motifs found in various cis-acting genetic elements, but are also known to provoke different types of genetic alterations. The instability of such motifs is clearly size-related and depends on their potential to adopt secondary structures known as hairpins and cruciforms. Here we studied the influence of palindrome size on recombination between two directly repeated copies of the yeast CYC1 gene leading to the loss of the intervening sequence ("pop-out" recombination). We show that palindromes inserted either within one copy or between the two copies of the CYC1 gene become recombinogenic only when they attain a certain critical size and we estimate this critical size to be about 70 bp. With the longest palindrome used in this study (150 bp) we observed a more than 20-fold increase in the pop-out recombination. In the sae2/com1 mutant the palindrome-stimulated recombination was completely abolished. Suppression of palindrome recombinogenicity may be crucial for the maintenance of genetic stability in organisms containing a significant number of large palindromes in their genomes, like humans.

Mus81-dependent double-strand DNA breaks at in vivo-generated cruciform structures in S. cerevisiae

Long DNA palindromes are implicated in chromosomal rearrangement, but their roles in the underlying molecular events remain a matter of conjecture. One notion is that palindromes induce DNA breaks after assuming a cruciform structure, the four-way DNA junction providing a target for cleavage by Holliday junction (HJ)-specific enzymes. Though compelling, few components of the "cruciform resolution" proposal are established. Here we address fundamental properties and genetic dependencies of palindromic DNA metabolism in eukaryotes. Plasmid-borne palindromes introduced into S. cerevisiae are site-specifically broken in vivo, and the breaks exhibit unique hallmarks of an HJ resolvase mechanism. In vivo resolution requires Mus81, for which the bacterial HJ resolvase RusA will substitute. These results provide confirmation of cruciform extrusion and resolution in the context of eukaryotic chromatin. Related observations are that, unchecked by a nuclease function provided by Mre11, episomal palindromes launch a self-perpetuating breakage-fusion-bridge-independent copy number increase termed "escape."

Evolution of a complex minisatellite DNA sequence

Minisatellites are tandem repeats of short DNA units widely distributed in genomes. However, the information on their dynamics in a phylogenetic context is very limited. Here we have studied the organization of the MsH43 locus in several species of primates and from these data we have reconstructed the evolutionary history of this complex minisatellite. Overall, with the exception of gibbon, MsH43 has an organization that is asymmetric, since the distribution of repeats is distinct between the 5' and 3' halves, and heterogeneous since there are many different repeats, some of them characteristic of each species. Inspection of the MsH43 arrays showed the existence of many duplications and deletions, suggesting the implication of slippage processes in the generation of polymorphism. Concerning the evolutionary history of this minisatellite, we propose that the birth of MsH43 may be situated before the divergence of Old World Monkeys since we found the existence of some MsH43 repeat motifs in prosimians and New World Monkeys. The analysis of MsH43 in apes revealed the existence of an evolutionary breakpoint in the pathway that originated African great apes and humans. Remarkably, human MsH43 is more homologous to orang-utan than to the corresponding sequence in gorilla and chimpanzee. This finding does not comply with the evolutionary paradigm that continuous alterations occur during the course of genome evolution. To adjust our results to the standard phylogeny of primates, we propose the existence of a wandering allele that was maintained almost unaltered during the period that extends between orang-utan and humans.

The MDR1/ABCB1 regional amplification in large inverted repeats with asymmetric sequences and microhomologies at the junction sites

A multidrug-resistant lung cancer cell line PTX250, established by treatment with the anti-cancer drug paclitaxel, has been demonstrated to have an increased copy number in the 7q21.12 region including the MDR1/ABCB1 gene. The amplicon is 2.7 megabases in size, and the copy number increase is 11-fold compared with the parental cell line. Here, we examined the amplicon structure and determined nucleotide sequences at both junctions of the amplicon. Fluorescence in situ hybridization analysis using an MDR1 probe demonstrated a cluster of fluorescent signals at the chromosomal end, suggesting an intra-chromosomal amplification. DNA fragments of both junctions were cloned and sequenced. The distal junction was a head-to-head fusion with a 4-base pair (bp) overlap separated by an asymmetric sequence of 1,265 bp, and the proximal junction was a tail-to-tail fusion with a 2-bp overlap intervened by an asymmetric sequence of 2,203 bp. These results suggest that the amplicon has a large palindromic structure with an asymmetric sequence and has been amplified through the breakage-fusion-bridge cycle. Specific sequences, which might be related to the occurrence of double-strand-breakages, were found at or near the junctions of the amplicon -- an inverted repeat in the distal junction and a highly AT-rich region near the proximal junction.

DNA palindromes with a modest arm length of greater, similar 20 base pairs are a significant target for recombinant adeno-associated virus vector integration in the liver, muscles, and heart in mice

Our previous study has shown that recombinant adeno-associated virus (rAAV) vector integrates preferentially in genes, near transcription start sites and CpG islands in mouse liver (H. Nakai, X. Wu, S. Fuess, T. A. Storm, D. Munroe, E. Montini, S. M. Burgess, M. Grompe, and M. A. Kay, J. Virol. 79:3606-3614, 2005). However, the previous method relied on in vivo selection of rAAV integrants and could be employed for the liver but not for other tissues. Here, we describe a novel method for high-throughput rAAV integration site analysis that does not rely on marker gene expression, selection, or cell division, and therefore it can identify rAAV integration sites in nondividing cells without cell manipulations. Using this new method, we identified and characterized a total of 997 rAAV integration sites in mouse liver, skeletal muscle, and heart, transduced with rAAV2 or rAAV8 vector. The results support our previous observations, but notably they have revealed that DNA palindromes with an arm length of greater, similar 20 bp (total length, greater, similar 40 bp) are a significant target for rAAV integration. Up to approximately 30% of total integration events occurred in the vicinity of DNA palindromes with an arm length of greater, similar 20 bp. Considering that DNA palindromes may constitute fragile genomic sites, our results support the notion that rAAV integrates at chromosomal sites susceptible to breakage or preexisting breakage sites. The use of rAAV to label fragile genomic sites may provide an important new tool for probing the intrinsic source of ongoing genomic instability in various tissues in animals, studying DNA palindrome metabolism in vivo, and understanding their possible contributions to carcinogenesis and aging.

The human genome-wide distribution of DNA palindromes

In this work, we performed a systematic study of perfect and nonspacer palindromes present in human genomic DNA, and we investigated palindrome distribution over the entire human genome and over the functional regions such as the exon, intron, intergenic, and upstream regions (2,000 bp upstream from translational start site). We found that 24 palindrome-abundant intervals are mostly located on G-bands, which condense early, replicate late, and are relatively A+T rich. In general, palindromes are overrepresented in introns but underrepresented in exons. Upstream region has enriched palindrome distribution, where palindromes can serve as transcription factor binding sites. We created a Human DNA Palindrome Database (HPALDB) which is accessible at http://vhp.ntu.edu.sg/hpaldb . It contains 12,556,994 entries covering all palindromes in the human genome longer than 6 bp. Queries can be performed in different ways. Each entry in the database is linked to its location on NCBI's human chromosome Map Viewer.