Instability of repetitive DNA sequences: The role of replication in multiple mechanisms

Rearrangements between tandem sequence homologies of various lengths are a major source of genomic change and can be deleterious to the organism. These rearrangements can result in either deletion or duplication of genetic material flanked by direct sequence repeats. Molecular genetic analysis of repetitive sequence instability in Escherichia coli has provided several clues to the underlying mechanisms of these rearrangements. We present evidence for three mechanisms of RecA-independent sequence rearrangements: simple replication slippage, sister-chromosome exchange-associated slippage, and single-strand annealing. We discuss the constraints of these mechanisms and contrast their properties with RecA-dependent homologous recombination. Replication plays a critical role in the two slipped misalignment mechanisms, and difficulties in replication appear to trigger rearrangements via all these mechanisms.

Stabilization of diverged tandem repeats by mismatch repair: evidence for deletion formation via a misaligned replication intermediate.

A functional methyl-directed mismatch repair pathway in Escherichia coli prevents the formation of deletions between 101-bp tandem repeats with 4% sequence divergence. Deletions between perfectly homologous repeats are unaffected. Deletion in both cases occurs independently of the homologous recombination gene, recA. Because the methyl-directed mismatch repair pathway detects and excises one strand of a mispaired duplex, an intermediate for RecA-independent deletion of tandem repeats must therefore be a heteroduplex formed between strands of each repeat. We find that MutH endonuclease, which in vivo incises specifically the newly replicated strand of DNA, and the Dam methylase, the source of this strand-discrimination, are required absolutely for the exclusion of "homeologous" (imperfectly homologous) tandem deletion. This supports the idea that the heteroduplex intermediate for deletion occurs during or shortly after DNA replication in the context of hemi-methylation. Our findings confirm a "replication slippage" model for deletion formation whereby the displacement and misalignment of the nascent strand relative to the repeated sequence in the template strand accomplishes the deletion.

Evolution of the secondary structure of the rRNA internal transcribed spacer 2 (ITS2) in hard ticks (Ixodidae, Arthropoda)

ITS2 sequences are used extensively in molecular taxonomy and population genetics of arthropods and other animals yet little is known about the molecular evolution of ITS2. We studied the secondary structure of ITS2 in species from each of the six main lineages of hard ticks (family Ixodidae). The ITS2 of these ticks varied in length from 679 bp in Ixodes scapularis to 1547 bp in Aponomma concolor. Nucleotide content varied also: the ITS2 of ticks from the Prostriata lineage (Ixodes spp.) had 46-49% GC whereas ITS2 sequences of ticks from the Metastriata lineage (all other hard ticks) had 61-62% GC. Despite variation in nucleotide sequence, the secondary structure of the ITS2 of all of these ticks apparently has five domains. Stems 1, 3, 4 and 5 of this secondary structure were obvious in all of the species studied. However, stem 2 was not always obvious despite the fact that it is flanked by highly conserved sequence motifs in the adjacent stems, stems 1 and 3. The ITS2 of hard ticks has apparently evolved mostly by increases and decreases in length of the nucleotide sequences, which caused increases, and decreases in the length of stems of the secondary structure. This is most obvious when stems of the secondary structures of the Prostriata (Ixodes spp.) are compared to those of the Metastriata (all other hard ticks). Increases in the size of the ITS2 may have been caused by replication slippage which generated large repeats, like those seen in Haemaphysalis humerosa and species from the Rhipicepalinae lineage, and the small repeats found in species from the other lineages of ticks.

Origin and evolution of homologous repeated sequences in the mitochondrial DNA control region of shrews.

The complete mitochondrial DNA (mtDNA) control region was amplified and directly sequenced in two species of shrew, Crocidura russula and Sorex araneus (Insectivora, Mammalia). The general organization is similar to that found in other mammals: a central conserved region surrounded by two more variable domains. However, we have found in shrews the simultaneous presence of arrays of tandem repeats in potential locations where repeats tend to occur separately in other mammalian species. These locations correspond to regions which are associated with a possible interruption of the replication processes, either at the end of the three-stranded D-loop structure or toward the end of the heavy-strand replication. In the left domain the repeated sequences (R1 repeats) are 78 bp long, whereas in the right domain the repeats are 12 bp long in C. russula and 14 bp long in S. araneus (R2 repeats). Variation in the copy number of these repeated sequences results in mtDNA control region length differences. Southern blot analysis indicates that level of heteroplasmy (more than one mtDNA form within an individual) differs between species. A comparative study of the R2 repeats in 12 additional species representing three shrew subfamilies provides useful indications for the understanding of the origin and the evolution of these homologous tandemly repeated sequences. An asymmetry in the distribution of variants within the arrays, as well as the constant occurrence of shorter repeated sequences flanking only one side of the R2 arrays, could be related to asymmetry in the replication of each strand of the mtDNA molecule. The pattern of sequence and length variation within and between species, together with the capability of the arrays to form stable secondary structures, suggests that the dominant mechanism involved in the evolution of these arrays in unidirectional replication slippage.

Dissecting the role of low-complexity regions in the evolution of vertebrate proteins

Low-complexity regions (LCRs) in proteins are tracts that are highly enriched in one or a few aminoacids. Given their high abundance, and their capacity to expand in relatively short periods of time through replication slippage, they can greatly contribute to increase protein sequence space and generate novel protein functions. However, little is known about the global impact of LCRs on protein evolution. We have traced back the evolutionary history of 2,802 LCRs from a large set of homologous protein families from H.sapiens, M.musculus, G.gallus, D.rerio and C.intestinalis. Transcriptional factors and other regulatory functions are overrepresented in proteins containing LCRs. We have found that the gain of novel LCRs is frequently associated with repeat expansion whereas the loss of LCRs is more often due to accumulation of amino acid substitutions as opposed to deletions. This dichotomy results in net protein sequence gain over time. We have detected a significant increase in the rate of accumulation of novel LCRs in the ancestral Amniota and mammalian branches, and a reduction in the chicken branch. Alanine and/or glycine-rich LCRs are overrepresented in recently emerged LCR sets from all branches, suggesting that their expansion is better tolerated than for other LCR types. LCRs enriched in positively charged amino acids show the contrary pattern, indicating an important effect of purifying selection in their maintenance. We have performed the first large-scale study on the evolutionary dynamics of LCRs in protein families. The study has shown that the composition of an LCR is an important determinant of its evolutionary pattern.

The structure of interrupted human AC microsatellites

Microsatellite lengths change over evolutionary time through a process of replication slippage. A recently proposed model of this process holds that the expansionary tendencies of slippage mutation are balanced by point mutations breaking longer microsatellites into smaller units and that this process gives rise to the observed frequency distributions of uninterrupted microsatellite lengths. We refer to this as the slippage/point-mutation theory. Here we derive the theory's predictions for interrupted microsatellites comprising regions of perfect repeats, labeled segments, separated by dinucleotide interruptions containing point mutations. These predictions are tested by reference to the frequency distributions of segments of AC microsatellite in the human genome, and several predictions are shown not to be supported by the data, as follows. The estimated slippage rates are relatively low for the first four repeats, and then rise initially linearly with length, in accordance with previous work. However, contrary to expectation and the experimental evidence, the inferred slippage rates decline in segments above 10 repeats. Point mutation rates are also found to be higher within microsatellites than elsewhere. The theory provides an excellent fit to the frequency distribution of peripheral segment lengths but fails to explain why internal segments are shorter. Furthermore, there are fewer microsatellites with many segments than predicted. The frequencies of interrupted microsatellites decline geometrically with microsatellite size measured in number of segments, so that for each additional segment, the number of microsatellites is 33.6% less. Overall we conclude that the detailed structure of interrupted microsatellites cannot be reconciled with the existing slippage/point-mutation theory of microsatellite evolution, and we suggest that microsatellites are stabilized by processes acting on interior rather than on peripheral segments.

Complex rearrangements in patients with duplications of MECP2 can occur by fork stalling and template switching

Duplication at the Xq28 band including the MECP2 gene is one of the most common genomic rearrangements identified in neurodevelopmentally delayed males. Such duplications are non-recurrent and can be generated by a non-homologous end joining (NHEJ) mechanism. We investigated the potential mechanisms for MECP2 duplication and examined whether genomic architectural features may play a role in their origin using a custom designed 4-Mb tiling-path oligonucleotide array CGH assay. Each of the 30 patients analyzed showed a unique duplication varying in size from similar to 250 kb to similar to 2.6 Mb. Interestingly, in 77% of these non-recurrent duplications, the distal breakpoints grouped within a 215 kb genomic interval, located 47 kb telomeric to the MECP2 gene. The genomic architecture of this region contains both direct and inverted low-copy repeat (LCR) sequences; this same region undergoes polymorphic structural variation in the general population. Array CGH revealed complex rearrangements in eight patients; in six patients the duplication contained an embedded triplicated segment, and in the other two, stretches of non-duplicated sequences occurred within the duplicated region. Breakpoint junction sequencing was achieved in four duplications and identified an inversion in one patient, demonstrating further complexity. We propose that the presence of LCRs in the vicinity of the MECP2 gene may generate an unstable DNA structure that can induce DNA strand lesions, such as a collapsed fork, and facilitate a Fork Stalling and Template Switching event producing the complex rearrangements involving MECP2.

Origins of sequence diversity in the malaria vaccine candidate merozoite surface protein-2 (MSP-2) in Amazonian isolates of Plasmodium falciparum

The recent evolution of Plasmodium falciparum is at odds with the extensive polymorphism found in most genes coding for antigens. Here, we examined the patterns and putative mechanisms of sequence diversification in the merozoite surface protein-2 (MSP-2), a major malarial repetitive surface antigen. We compared the msp-2 gene sequences from closely related clones derived from sympatric parasite isolates from Brazilian Amazonia and used microsatellite typing to examine, in these same clones, the haplotype background of chromosome 2, where msp-2 is located. We found examples of msp-2 sequence rearrangements putatively created by nonreciprocal recombinational events, such as replication slippage and gene conversion, while maintaining the chromosome haplotype. We conclude that these nonreciprocal recombination events may represent a major source of antigenic diversity in MSP-2 in P falciparum populations with low rates of classical meiotic recombination. (c) 2006 Elsevier B.V. All rights reserved.

In-frame triplet deletions in RHD alter the D antigen phenotype

BACKGROUND: The deletion of three adjacent nucleotides in an exon may cause the lack of a single amino acid, while the protein sequence remains otherwise unchanged. Only one such in-frame deletion is known in the two RH genes, represented by the RHCE allele ceBP expressing a "very weak e antigen." STUDY DESIGN AND METHODS: Blood donor samples were recognized because of discrepant results of D phenotyping. Six samples came from Switzerland and one from Northern Germany. The molecular structures were determined by genomic DNA nucleotide sequencing of RHD. RESULTS: Two different variant D antigens were explained by RHD alleles harboring one in-frame triplet deletion each. Both single-amino-acid deletions led to partial D phenotypes with weak D antigen expression. Because of their D category V-like phenotypes, the RHD(Arg229del) allele was dubbed DVL-1 and the RHD(Lys235del) allele DVL-2. These in-frame triplet deletions are located in GAGAA or GAAGA repeats of the RHD exon 5. CONCLUSION: Partial D may be caused by a single-amino-acid deletion in RhD. The altered RhD protein segments in DVL types are adjacent to the extracellular loop 4, which constitutes one of the most immunogenic parts of the D antigen. These RhD protein segments are also altered in all DV, which may explain the similarity in phenotype. At the nucleotide level, the triplet deletions may have resulted from replication slippage. A total of nine amino acid positions in an Rhesus protein may be affected by this mechanism.

Unexpected Diversity of Cellular Immune Responses against Nef and Vif in HIV-1-Infected Patients Who Spontaneously Control Viral Replication

Background: HIV-1-infected individuals who spontaneously control viral replication represent an example of successful containment of the AIDS virus. Understanding the anti-viral immune responses in these individuals may help in vaccine design. However, immune responses against HIV-1 are normally analyzed using HIV-1 consensus B 15-mers that overlap by 11 amino acids. Unfortunately, this method may underestimate the real breadth of the cellular immune responses against the autologous sequence of the infecting virus. Methodology and Principal Findings: Here we compared cellular immune responses against nef and vif-encoded consensus B 15-mer peptides to responses against HLA class I-predicted minimal optimal epitopes from consensus B and autologous sequences in six patients who have controlled HIV-1 replication. Interestingly, our analysis revealed that three of our patients had broader cellular immune responses against HLA class I-predicted minimal optimal epitopes from either autologous viruses or from the HIV-1 consensus B sequence, when compared to responses against the 15-mer HIV-1 type B consensus peptides. Conclusion and Significance: This suggests that the cellular immune responses against HIV-1 in controller patients may be broader than we had previously anticipated.

Wolbachia replication and host cell division in Aedes albopictus

Wolbachia pipientis is an obligate intracellular endosymbiont of a range of arthropod species. The microbe is best known for its manipulations of host reproduction that include inducing cytoplasmic incompatibility, parthenogenesis, feminization, and male-killing. Like other vertically transmitted intracellular symbionts, Wolbachiarsquos replication rate must not outpace that of its host cells if it is to remain benign. The mosquito Aedes albopictus is naturally infected both singly and doubly with different strains of Wolbachia pipientis. During diapause in mosquito eggs, no host cell division is believed to occur. Further development is triggered only by subsequent exposure of the egg to water. This study uses diapause in Wolbachia-infected Aedes albopictus eggs to determine whether symbiont replication slows or stops when host cell division ceases or whether it continues at a low but constant rate. We have shown that Wolbachia densities in eggs are greatest during embryonation and then decline throughout diapause, suggesting that Wolbachia replication is dependent on host cell replication.

Site-directed mutants of RTP of Bacillus subtilis and the mechanism of replication fork arrest

DNA replication fork arrest during the termination phase of chromosome replication in Bacillus subtilis is brought about by the replication terminator protein (RTP) bound to specific DNA terminator sequences (Tev sites) distributed throughout the terminus region. An attractive suggestion by others was that crucial to the functioning of the RTP-Ter complex is a specific interaction between RTP positioned on the DNA and the helicase associated with the approaching replication fork. Ln support of this was the behaviour of two site-directed mutants of RTP. They appeared to bind Ter DNA normally but were ineffective in fork arrest as ascertained by in vitro Escherichia coli DnaB helicase and replication assays. We describe here a system for assessing the fork-arrest behaviour of RTP mutants in a bona fide in vivo assay in B. subtilis. One of the previously studied mutants, RTP.Y33N, was non-functional in fork arrest in vivo, as predicted. But through extensive analyses, this RTP mutant was shown to be severely defective in binding to Ter DNA, contrary to expectation. Taken in conjunction with recent findings on the other mutant (RTP.E30K), it is concluded that there is as yet no substantive evidence from the behaviour of RTP mutants to support the Rm-helicase interaction model for fork arrest. In an extension of the present work on RTP.Y33N, we determined the dissociation rates of complexes formed by wild-type (wt) RTP and another RTP mutant with various terminator sequences. The functional wtRTP-TerI complex was quite stable (half-life of 182 minutes), reminiscent of the great stability of the E. coli Tus-Ter complex. More significant were the exceptional stabilities of complexes comprising wtRTP and an RTP double-mutant (E39K.R42Q) bound to some particular terminator sequences. From the measurement of in vivo fork-arrest activities of the various complexes, it is concluded that the stability (half-life) of the whole RTP-Ter complex is not the overriding determinant of arrest, and that the RTP-Ter complex must be actively disrupted, or RTP removed, by the action of the approaching replication fork. (C) 1999 Academic Press.

A tale of two terminators: Crystal structures sharpen the debate on DNA replication fork arrest mechanisms

The structure of the Tus-Ter DNA replication fork arrest complex of Escherichia coli reveals a novel architecture for the bound Tus protein and a new type of DNA-binding motif, The structure of the complex may explain how Tus can block movement of a replication fork approaching from one direction and not the other.

Reorganization of terminator DNA upon binding replication terminator protein: Implications for the functional replication fork arrest complex

Termination of DNA replication in Bacillus subtilis involves the polar arrest of replication forks by a specific complex formed between the replication terminator protein (RTP) and DNA terminator sites. While determination of the crystal structure of RTP has facilitated our understanding of how a single RTP dimer interacts with terminator DNA, additional information is required in order to understand the assembly of a functional fork arrest complex, which requires an interaction between two RTP dimers and the terminator site. In this study, we show that the conformation of the major B. subtilis DNA terminator, Terl, becomes considerably distorted upon binding RTP. Binding of the first dimer of RTP to the B site of Terl causes the DNA to become slightly unwound and bent by similar to 40 degrees. Binding of a second dimer of RTP to the A site causes the bend angle to increase to similar to 60 degrees. We have used this new data to construct two plausible models that might explain how the ternary terminator complex can block DNA replication in a polar manner, in the first model, polarity of action is a consequence of the two RTP-DNA half-sites having different conformations. These different conformations result from different RTP-DNA contacts at each half-site (due to the intrinsic asymmetry at the terminator DNA), as well as interactions (direct or indirect) between the RTP dimers on the DNA. In the second model, polar fork arrest activity is a consequence of the different affinities of RTP for the A and B sites of the terminator DNA, modulated significantly by direct or indirect interactions between the RTP dimers.