Tandem duplication. A novel type of triplet repeat instability

SVSR: A Program to Simulate Structural Variations and Generate Sequencing Reads for Multiple Platforms

Structural variation accounts for a major fraction of mutations in the human genome and confers susceptibility to complex diseases. Next generation sequencing along with the rapid development of computational methods provides a cost-effective procedure to detect such variations. Simulation of structural variations and sequencing reads with real characteristics is essential for benchmarking the computational methods. Here, we develop a new program, SVSR, to simulate five types of structural variations (indels, tandem duplication, CNVs, inversions, and translocations) and SNPs for the human genome and to generate sequencing reads with features from popular platforms (Illumina, SOLiD, 454, and Ion Torrent). We adopt a selection model trained from real data to predict copy number states, starting from the first site of a particular genome to the end. Furthermore, we utilize references of microbial genomes to produce insertion fragments and design probabilistic models to imitate inversions and translocations. Moreover, we create platform-specific errors and base quality profiles to generate normal, tumor, or normal-tumor mixture reads. Experimental results show that SVSR could capture more features that are realistic and generate datasets with satisfactory quality scores. SVSR is able to evaluate the performance of structural variation detection methods and guide the development of new computational methods.

Structural basis for the identification of an i-motif tetraplex core with a parallel-duplex junction as a structural motif in CCG triplet repeats

CCG triplet repeats can fold into tetraplex structures, which are associated with the expansion of (CCG)n trinucleotide sequences in certain neurological diseases. These structures are stabilized by intertwining i-motifs. However, the structural basis for tetraplex i-motif formation in CCG triplet repeats remains largely unknown. We report the first crystal structure of a CCG-repeat sequence, which shows that two dT(CCG)3 A strands can associate to form a tetraplex structure with an i-motif core containing four C:C(+) pairs flanked by two G:G homopurine base pairs as a structural motif. The tetraplex core is attached to a short parallel-stranded duplex. Each hairpin itself contains a central CCG loop in which the nucleotides are flipped out and stabilized by stacking interactions. The helical twists between adjacent cytosine residues of this structure in the i-motif core have an average value of 30°, which is greater than those previously reported for i-motif structures.

The structural basis of actinomycin D-binding induces nucleotide flipping out, a sharp bend and a left-handed twist in CGG triplet repeats

The potent anticancer drug actinomycin D (ActD) functions by intercalating into DNA at GpC sites, thereby interrupting essential biological processes including replication and transcription. Certain neurological diseases are correlated with the expansion of (CGG)_n trinucleotide sequences, which contain many contiguous GpC sites separated by a single G:G mispair. To characterize the binding of ActD to CGG triplet repeat sequences, the structural basis for the strong binding of ActD to neighbouring GpC sites flanking a G:G mismatch has been determined based on the crystal structure of ActD bound to ATGCGGCAT, which contains a CGG triplet sequence. The binding of ActD molecules to GCGGC causes many unexpected conformational changes including nucleotide flipping out, a sharp bend and a left-handed twist in the DNA helix via a two site-binding model. Heat denaturation, circular dichroism and surface plasmon resonance analyses showed that adjacent GpC sequences flanking a G:G mismatch are preferred ActD-binding sites. In addition, ActD was shown to bind the hairpin conformation of (CGG)₁₆ in a pairwise combination and with greater stability than that of other DNA intercalators. Our results provide evidence of a possible biological consequence of ActD binding to CGG triplet repeat sequences.

The binding of the Co(II) complex of dimeric chromomycin A3 to GC sites with flanking G:G mismatches

Some neurological diseases are correlated with expansion of (CXG)n trinucleotide repeats, which contain many contiguous GpC flanked by mismatched X/X base pair. This study focused on the binding of the Co(II) complex of dimeric chromomycin A3(Chro), Co(II)(Chro)2, to DNA with CXG trinucleotide repeats. The present study showed that GC sites with flanking G:G mismatches provide an excellent binding site for Co(II)(Chro)2 as shown by surface plasmon resonance and fluorescence analysis, compared to GC sites with flanking A:A, T:T, or C:C mismatches. In addition, we measured the ability of Co(II)(Chro)2 to act on the hairpin DNA of (CGG)16. We observed that Co(II)(Chro)2 could stabilize and trap the cruciform conformation of (CGG)16. Furthermore, two Co(II)(Chro)2 molecules may bind at the two GpC sites separated by at least one GC site in the hairpin structure of (CGG)16. In a synthetic self-priming DNA model, 5'-(CGG)16(CCG)6-3', Co(II)(Chro)2 can interfere with the expansion process of CGG triplet repeats, as shown by a gel electrophoretic expansion assay. Here, we first report the acting of Co(II)(Chro)2, the groove-binding drug, to trinucleotide repeats. Our results provide the possible biological consequence of Co(II)(Chro)2 bound to CGG triplet repeat sequences.

Chromosomal model for analysis of a long CTG/CAG tract stability in wild-type Escherichia coli and its nucleotide excision repair mutants

Many human hereditary neurological diseases, including fragile X syndrome, myotonic dystrophy, and Friedreich's ataxia, are associated with expansions of the triplet repeat sequences (TRS) (CGG/CCG, CTG/CAG, and GAA/TTC) within or near specific genes. Mechanisms that mediate mutations of TRS include DNA replication, repair, and gene conversion and (or) recombination. The involvement of the repair systems in TRS instability was investigated in Escherichia coli on plasmid models, and the results showed that the deficiency of some nucleotide excision repair (NER) functions dramatically affects the stability of long CTG inserts. In such models in which there are tens or hundreds of plasmid molecules in each bacterial cell, repetitive sequences may interact between themselves and according to a recombination hypothesis, which may lead to expansions and deletions within such repeated tracts. Since one cannot control interaction between plasmids, it is also sometimes difficult to give precise interpretation of the results. Therefore, using modified lambda phage (lambdaInCh), we have constructed a chromosomal model to study the instability of trinucleotide repeat sequences in E. coli. We have shown that the stability of (CTG/CAG)68 tracts in the bacterial chromosome is influenced by mutations in NER genes in E. coli. The absence of the uvrC or uvrD gene products greatly enhances the instability of the TRS in the chromosome, whereas the lack of the functional UvrA or UvrB proteins causes substantial stabilization of (CTG/CAG) tracts.

Mutations causing Friedreich ataxia

Friedreich ataxia (FRDA) is an autosomal recessive degenerative disease primarily affecting the nervous system and the heart. The identification of the FRDA gene (FXN) encoding the mitochondrial protein frataxin allowed to better define the clinical and pathological spectrum of the disease, to investigate pathogenesis and to eventually propose novel treatments for this so far incurable disease. Frataxin has homologs in all eukaryotes and in gram-negative bacteria. Studies in yeast and mouse models and biochemical investigations indicate a role in the assembly of iron–sulfur clusters in the mitochondrion. Frataxin deficiency leads to abnormal mitochondrial iron metabolism, decreased activities of iron–sulfur cluster-containing enzymes, reduced oxidative phosphorylation and possibly increased oxidative stress. Complete absence of frataxin is lethal at an embryonic level in all investigated multicellular organisms, in yeast it leads to loss of mitochondrial function and mitochondrial DNA. Patients with FRDA have a unique mutation that results in low levels of structurally and functionally normal frataxin. The mutation is an unstable hyperexpansion of a GAA triplet-repeat sequence in the first intron of the gene. Pathological lengths of the GAA repeat adopt in vitro a non-B, triple helical DNA structure that inhibits transcription. In the nucleus, chromatin condenses at and around the GAA repeat into a heterochromatic state that is associated with repressed transcription. It is not yet known whether the triple helical structure forms in the nucleus and has a role in triggering chromatin condensation. Rare mutations, always occurring in heterozygosity with the GAA repeat expansion in FRDA patients, alter the coding frataxin sequence and result in non-functional or absent protein. Restoring frataxin transcription is an exciting perspective for treatment. Preliminary results in cell models suggest that some inhibitors of histone deacetylases may have this effect.

Novel repressor of the human FMR1 gene - identification of p56 human (GCC)(n)-binding protein as a Krüppel-like transcription factor ZF5

A series of relatively short (GCC)(n) triplet repeats (n = 3-30) located within regulatory regions of many mammalian genes may be considered as putative cis-acting transcriptional elements (GCC-elements). Fragile X-mental retardation syndrome is caused by an expansion of (GCC)(n) triplet repeats within the 5'-untranslated region of the human fragile X-mental retardation 1 (FMR1) gene. The present study aimed to characterize a novel human (GCC)(n)-binding protein and investigate its possible role in the regulation of the FMR1 gene. A novel human (GCC)(n)-binding protein, p56, was isolated and identified as a Krüppel-like transcription factor, ZF5, by MALDI-TOF analysis. The capacity of ZF5 to specifically interact with (GCC)(n) triplet repeats was confirmed by the electrophoretic mobility shift assay with purified recombinant ZF5 protein. In cotransfection experiments, ZF5 overexpression repressed activity of the GCC-element containing mouse ribosomal protein L32 gene promoter. Moreover, RNA interference assay results showed that endogenous ZF5 acts as a repressor of the human FMR1 gene. Thus, these data identify a new class of ZF5 targets, a subset of genes containing GCC-elements in their regulatory regions, and raise the question of whether transcription factor ZF5 is implicated in the pathogenesis of fragile X syndrome.

The Myotonic Dystrophy Type 1 Triplet Repeat Sequence Induces Gross Deletions and Inversions

The capacity of (CTG.CAG)n and (GAA.TTC)n repeat tracts in plasmids to induce mutations in DNA flanking regions was evaluated in Escherichia coli. Long repeats of these sequences are involved in the etiology of myotonic dystrophy type 1 and Friedreich's ataxia, respectively. Long (CTG.CAG)n (where n = 98 and 175) caused the deletion of most, or all, of the repeats and the flanking GFP gene. Deletions of 0.6-1.8 kbp were found as well as inversions. Shorter repeat tracts (where n = 0 or 17) were essentially inert, as observed for the (GAA.TTC)176-containing plasmid. The orientation of the triplet repeat sequence (TRS) relative to the unidirectional origin of replication had a pronounced effect, signaling the participation of replication and/or repair systems. Also, when the TRS was transcribed, the level of deletions was greatly elevated. Under certain conditions, 30-50% of the products contained gross deletions. DNA sequence analyses of the breakpoint junctions in 47 deletions revealed the presence of 1-8-bp direct or inverted homologies in all cases. Also, the presence of non-B folded conformations (i.e. slipped structures, cruciforms, or triplexes) at or near the breakpoints was predicted in all cases. This genetic behavior, which was previously unrecognized for a TRS, may provide the basis for a new type of instability of the myotonic dystrophy protein kinase (DMPK) gene in patients with a full mutation.

Small de novo duplication in the repeat region of the TATA-box-binding protein gene manifest with a phenotype similar to variant Creutzfeldt-Jakob disease

A 20-year-old North American patient developed rapidly progressive cognitive decline and pronounced ataxia, a phenotype compatible with prion disease. No structural changes were found in the PRNP gene, which excludes genetic prion disease, but the patient's PRNP codon 129 Met/Met genotype is known to predispose to variant Creutzfeldt-Jakob disease (vCJD). Further studies identified an expanded allele with 55 CAG/CAA repeats in the TBP gene. The increase of trinucleotide repeat number in the coding region of the TBP gene has previously been associated with spinocerebellar ataxia type 17 (SCA17). The patient's unaffected parents and siblings show normal-size TBP alleles with 37-38 repeats. Haplotype and nucleotide sequence analyses clearly indicate that the mutation has occurred de novo on a paternal chromosome by insertion/duplication of a (CAA)(CAG)(CAA)(CAG)(15) sequence. This report presents a second fully investigated sporadic case of SCA17 occurring as a result of a DNA rearrangement within the polymorphic TBP trinucleotide repeat region. Our findings suggest that patients suspected of vCJD should undergo testing for SCA17, Huntington's disease and other neurodegenerative disorders having phenotypic similarities with vCJD.

DNA Double-strand Breaks Induce Deletion of CTG·CAG Repeats in an Orientation-dependent Manner in Escherichia coli

The influences of double-strand breaks (DSBs) within a triplet repeat sequence on its genetic instabilities (expansions and deletions) related to hereditary neurological diseases was investigated. Plasmids containing 43 or 70 CTG.CAG repeats or 43 CGG.CCG repeats were linearized in vitro near the center of the repeats and were transformed into parental, RecA-dependent homologous recombination-deficient, or RecBC exonuclease-deficient Escherichia coli. The resulting repair process considerably increased deletion of the repeating sequence compared to the circular DNA controls. Unexpectedly, the orientation of the insert relative to the unidirectional ColE1 origin of replication affected the amount of instability generated during the repair of the DSB. When the CTG strand was the template for lagging-strand synthesis, instability was increased, most markedly in the recA- strain. Results indicated that RecA and/or RecBC might play a role in DSB repair within the triplet repeat. Altering the length, orientation, and sequence composition of the triplet repeat suggested an important role of DNA secondary structures during repair intermediates. Hence, we hypothesize that ColE1 origin-dependent replication was involved during the repair of the DSB. A model is presented to explain the mechanisms of the observed genetic instabilities.

Structure-dependent Recombination Hot Spot Activity of GAA·TTC Sequences from Intron 1 of the Friedreich's Ataxia Gene

The recombinational properties of long GAA.TTC repeating sequences were analyzed in Escherichia coli to gain further insights into the molecular mechanisms of the genetic instability of this tract as possibly related to the etiology of Friedreich's ataxia. Intramolecular and intermolecular recombination studies showed that the frequency of recombination between the GAA.TTC tracts was as much as 15 times higher than the non-repeating control sequences. Homologous, intramolecular recombination between GAA.TTC tracts and GAAGGA.TCCTTC repeats also occurred with a very high frequency (approximately 0.8%). Biochemical analyses of the recombination products demonstrated the expansions and deletions of the GAA.TTC repeats. These results, together with our previous studies on the CTG.CAG sequences, suggest that the recombinational hot spot characteristics may be a common feature of all triplet repeat sequences. Unexpectedly, we found that the recombination properties of the GAA.TTC tracts were unique, compared with CTG.CAG repeats, because they depended on the DNA secondary structure polymorphism. Increasing the length of the GAA.TTC repeats decreased the intramolecular recombination frequency between these tracts. Also, a correlation was found between the propensity of the GAA.TTC tracts to adopt the sticky DNA conformation and the inhibition of intramolecular recombination. The use of novobiocin to modulate the intracellular DNA topology, i.e. the lowering of the negative superhelical density, repressed the formation of the sticky DNA structure, thereby restoring the expected positive correlation between the length of the GAA.TTC tracts and the frequency of intramolecular recombination. Hence, our results demonstrate that sticky DNA exists and functions in E. coli.

Long CTG tracts from the myotonic dystrophy gene induce deletions and rearrangements during recombination at the APRT locus in CHO cells

Expansion of CTG triplet repeats in the 3' untranslated region of the DMPK gene causes the autosomal dominant disorder myotonic dystrophy. Instability of CTG repeats is thought to arise from their capacity to form hairpin DNA structures. How these structures interact with various aspects of DNA metabolism has been studied intensely for Escherichia coli and Saccharomyces cerevisiae but is relatively uncharacterized in mammalian cells. To examine the stability of (CTG)(17), (CTG)(98), and (CTG)(183) repeats during homologous recombination, we placed them in the second intron of one copy of a tandemly duplicated pair of APRT genes. Cells selected for homologous recombination between the two copies of the APRT gene displayed distinctive patterns of change. Among recombinants from cells with (CTG)(98) and (CTG)(183), 5% had lost large numbers of repeats and 10% had suffered rearrangements, a frequency more than 50-fold above normal levels. Analysis of individual rearrangements confirmed the involvement of the CTG repeats. Similar changes were not observed in proliferating (CTG)(98) and (CTG)(183) cells that were not recombinant at APRT. Instead, they displayed high frequencies of small changes in repeat number. The (CTG)(17) repeats were stable in all assays. These studies indicate that homologous recombination strongly destabilizes long tracts of CTG repeats.

Molecular mechanisms of TRS instability

To date several neurodegenerative disorders including myotonic dystrophy, Huntington's disease, Kennedy's disease, fragile X syndrome, spinocerebellar ataxias or Friedreich's ataxia have been linked to the expanding trinucleotide sequences. Although phenotypic features vary among these debilitating diseases, the structural abnormalities of the triplet repeat containing DNA sequences is the primary cause for all of these disorders. Expansions of the CAG repeat within coding regions of miscellaneous genes result in the synthesis of aberrant proteins containing enormously long polyglutamine stretches. Such proteins acquire toxic functions and/or may direct cells into the apoptotic cycle. On the other hand, massive expansions of various triplet repeats (i.e., CTG/CAG, CGG/CCG/, GAA/TTC) inside the noncoding regions lead to the silencing of transcription and therefore affect expression of the adjacent genes. The repetitive character of TRS allows stretches of such tracts to form slipped-stranded structures, self-complementary hairpins, triplexes or more complex configurations called "sticky DNA", which are not equally processed by some cellular mechanisms, as compared to random DNA. It is likely that the instability of the short TRS (below the threshold level) occurs due to the SILC pathway, which is driven by the DNA slippage. Accumulation of the short expansions leads to the disease premutation state where the MLC pathway becomes predominant. Independent of which mechanism is involved in the MLC pathway (replication, transcription, repair or recombination) the process of complementary strand synthesis is crucial for the TRS instability. Generally, dependent on the location of the tract which has higher potential to form secondary DNA structure, further processing of such tract may result in expansions (secondary structure formed at the newly synthesized strand) or deletions (structure present on the template strand). Analyses of molecular mechanisms of the TRS genetic instability using bacteria, yeast, cell lines and transgenic animals as models allowed the scientists to better understand the role of some major cellular processes in the development of neurodegenerative disorders in humans. However, it is necessary to remember that most of these investigations were focused on the involvement of each particular process separately. Much less of this work though was dedicated to the search for the interactions between such cellular systems that in effect could result in different rate of TRS expansions. Thus, more intensive studies are necessary in order to fully understand the phenomenon ofthe dynamic mutations leading to the human hereditary neurodegenerative diseases.

Crystal structure of actinomycin D bound to the CTG triplet repeat sequences linked to neurological diseases

The potent anticancer drug actinomycin D (ActD) acts by binding to DNA GpC sequences, thereby interfering with essential biological processes including replication, transcription and topoisomerase. Certain neurological diseases are correlated with expansion of (CTG)n trinucleotide sequences, which contain many contiguous GpC sites separated by a single base pair. In order to characterize the binding of ActD to CTG triplet repeat sequences, we carried out heat denaturation and CD analyses, which showed that adjacent GpC sequences flanking a T:T mismatch are preferred ActD-binding sites, and that ActD binding results in a conformational transition to A-type structure. The structural basis of the strong binding of ActD to neighboring GpC sites flanking a T:T mismatch was provided by the crystal structure of ActD bound to ATGCTGCAT, which contains a CTG triplet sequence. Binding of two ActD molecules to GCTGC causes a kink in the DNA helix. In addition, using a synthetic self-priming DNA model, 5'-(CAG)4(CTG)(16)-3', we observed that ActD can trap the cruciform or duplexes of (CTG)n and interfere with the expansion process of CTG triplet repeats as shown by gel electrophoretic expansion assay. Our results may provide the possible biological consequence of ActD bound to CTG triplet repeat sequences.

Specific interaction of p53 with target binding sites is determined by DNA conformation and is regulated by the C-terminal domain

Transcriptional activation of p53-regulated genes is initiated by sequence-specific DNA binding of p53 to target binding sites. Regulation of sequence-specific DNA binding is complex and occurs at various levels. We demonstrate that DNA topology is an important parameter for regulating the selective and highly specific interaction of p53 with its target binding sites. Specific binding of wild-type p53 is greatly enhanced when cognate binding sites are present in a non-linear stem-loop conformation. The C-terminal domain plays a key role in regulating the specific interactions of p53 with target binding sites in a DNA conformation-dependent manner. The C-terminal domain is required for binding to target sites in a non-linear DNA conformation in contrast to the strong inhibitory effects of the C terminus on p53 interaction with linear DNA. We propose that selective binding of p53 to various promoters may be determined by the DNA conformation within p53 cognate sites.

Long CTG.CAG repeats from myotonic dystrophy are preferred sites for intermolecular recombination

Homologous recombination was shown to enable the expansion of CTG.CAG repeat sequences. Other prior investigations revealed the involvement of replication and DNA repair in these genetic instabilities. Here we used a genetic assay to measure the frequency of homologous intermolecular recombination between two CTG.CAG tracts. When compared with non-repeating sequences of similar lengths, long (CTG.CAG)(n) repeats apparently recombine with an approximately 60-fold higher frequency. Sequence polymorphisms that interrupt the homogeneity of the CTG.CAG repeat tracts reduce the apparent recombination frequency as compared with the pure uninterrupted repeats. The orientation of the repeats relative to the origin of replication strongly influenced the apparent frequency of recombination. This suggests the involvement of DNA replication in the recombination process of triplet repeats. We propose that DNA polymerases stall within the CTG.CAG repeat tracts causing nicks or double-strand breaks that stimulate homologous recombination. The recombination process is RecA-dependent.

Long CTG.CAG repeat sequences markedly stimulate intramolecular recombination

Previous studies have shown that homologous recombination is a powerful mechanism for generation of massive instabilities of the myotonic dystrophy CTG.CAG sequences. However, the frequency of recombination between the CTG.CAG tracts has not been studied. Here we performed a systematic study on the frequency of recombination between these sequences using a genetic assay based on an intramolecular plasmid system in Escherichia coli. The rate of intramolecular recombination between long CTG.CAG tracts oriented as direct repeats was extraordinarily high; recombinants were found with a frequency exceeding 12%. Recombination occurred in both RecA(+) and RecA(-) cells but was approximately 2-11 times higher in the recombination proficient strain. Long CTG.CAG tracts recombined approximately 10 times more efficiently than non-repeating control sequences of similar length. The recombination frequency was 60-fold higher for a pair of (CTG.CAG)(165) tracts compared with a pair of (CTG.CAG)(17) sequences. The CTG.CAG sequences in orientation II (CTG repeats present on a lagging strand template) recombine approximately 2-4 times more efficiently than tracts of identical length in the opposite orientation relative to the origin of replication. This orientation effect implies the involvement of DNA replication in the intramolecular recombination between CTG.CAG sequences. Thus, long CTG.CAG tracts are hot spots for genetic recombination.

Expansion/contraction of mammalian mitochondrial DNA repeats in Escherichia coli mimics the mitochondrial heteroplasmy

Length polymorphism due to tandem repeats is a common feature in animal mitochondrial DNA. The rabbit mitochondrial genome contains a 20 bp repeat domain, which generates a general heteroplasmic state. The observed polymorphic patterns suggest a dynamic equilibrium between gain and loss of units that maintains the copy number in the range 3-19 repeat units. In the apparent absence of recombination, slipped-strand mispairing during replication appears to be the primary cause of additions and deletions. To investigate this hypothesis we have set up a plasmid assay in Escherichia coli. A variable number of repeat units was inserted into a plasmid in both orientations relative to the colE1 origin of replication. Our data show that (i) a minimum unit number (>3) is necessary to generate length polymorphs, (ii) the number of events increases with the length tract, (iii) an excess of additions over deletions is found when the copy number is less than 10 and the trend is reversed when it is over 10, (iv) the frequency of deletions-additions is dependent on the orientation, (v) the polymorphism patterns are different according to the orientation. The length polymorphic pattern generated in the bacteria, in one orientation, mimics that observed in the mitochondria, suggesting that slipped mispairing between repeated sequences during DNA replication is responsible for the mitochondrial heteroplasmic state.