Expansions and contractions in 36-bp minisatellites by gene conversion in yeast

Genome-wide characterization of human minisatellite VNTRs: population-specific alleles and gene expression differences

Variable Number Tandem Repeats (VNTRs) are tandem repeat (TR) loci that vary in copy number across a population. Using our program, VNTRseek, we analyzed human whole genome sequencing datasets from 2770 individuals in order to detect minisatellite VNTRs, i.e., those with pattern sizes ≥7 bp. We detected 35 638 VNTR loci and classified 5676 as commonly polymorphic (i.e. with non-reference alleles occurring in >5% of the population). Commonly polymorphic VNTR loci were found to be enriched in genomic regions with regulatory function, i.e. transcription start sites and enhancers. Investigation of the commonly polymorphic VNTRs in the context of population ancestry revealed that 1096 loci contained population-specific alleles and that those could be used to classify individuals into super-populations with near-perfect accuracy. Search for quantitative trait loci (eQTLs), among the VNTRs proximal to genes, indicated that in 187 genes expression differences correlated with VNTR genotype. We validated our predictions in several ways, including experimentally, through the identification of predicted alleles in long reads, and by comparisons showing consistency between sequencing platforms. This study is the most comprehensive analysis of minisatellite VNTRs in the human population to date.

The role of the Mre11-Rad50-Nbs1 complex in double-strand break repair-facts and myths

Homologous recombination (HR) initiates double-strand break (DSB) repair by digesting 5'-termini at DSBs, the biochemical reaction called DSB resection, during which DSBs are processed by nucleases to generate 3' single-strand DNA. Rad51 recombinase polymerizes along resected DNA, and the resulting Rad51-DNA complex undergoes homology search. Although DSB resection by the Mre11 nuclease plays a critical role in HR in Saccharomyces cerevisiae, it remains elusive whether DSB resection by Mre11 significantly contributes to HR-dependent DSB repair in mammalian cells. Depletion of Mre11 decreases the efficiency of DSB resection only by 2- to 3-fold in mammalian cells. We show that although Mre11 is required for efficient HR-dependent repair of ionizing-radiation-induced DSBs, Mre11 is largely dispensable for DSB resection in both chicken DT40 and human TK6 B cell lines. Moreover, a 2- to 3-fold decrease in DSB resection has virtually no impact on the efficiency of HR. Thus, although a large number of researchers have reported the vital role of Mre11-mediated DSB resection in HR, the role may not explain the very severe defect in HR in Mre11-deficient cells, including their lethality. We here show experimental evidence for the additional roles of Mre11 in (i) elimination of chemical adducts from DSB ends for subsequent DSB repair, and (ii) maintaining HR intermediates for their proper resolution.

The role of variable DNA tandem repeats in bacterial adaptation

DNA tandem repeats (TRs), also designated as satellite DNA, are inter- or intragenic nucleotide sequences that are repeated two or more times in a head-to-tail manner. Because TR tracts are prone to strand-slippage replication and recombination events that cause the TR copy number to increase or decrease, loci containing TRs are hypermutable. An increasing number of examples illustrate that bacteria can exploit this instability of TRs to reversibly shut down or modulate the function of specific genes, allowing them to adapt to changing environments on short evolutionary time scales without an increased overall mutation rate. In this review, we discuss the prevalence and distribution of inter- and intragenic TRs in bacteria and the mechanisms of their instability. In addition, we review evidence demonstrating a role of TR variations in bacterial adaptation strategies, ranging from immune evasion and tissue tropism to the modulation of environmental stress tolerance. Nevertheless, while bioinformatic analysis reveals that most bacterial genomes contain a few up to several dozens of intra- and intergenic TRs, only a small fraction of these have been functionally studied to date.

Break-induced replication is highly inaccurate

DNA replication initiated by one-ended homologous recombination at a double-strand break is highly inaccurate, as it greatly stimulates frameshift mutations over the entire path of the replication fork.

DNA must be synthesized for purposes of genome duplication and DNA repair. While the former is a highly accurate process, short-patch synthesis associated with repair of DNA damage is often error-prone. Break-induced replication (BIR) is a unique cellular process that mimics normal DNA replication in its processivity, rate, and capacity to duplicate hundreds of kilobases, but is initiated at double-strand breaks (DSBs) rather than at replication origins. Here we employed a series of frameshift reporters to measure mutagenesis associated with BIR in Saccharomyces cerevisiae. We demonstrate that BIR DNA synthesis is intrinsically inaccurate over the entire path of the replication fork, as the rate of frameshift mutagenesis during BIR is up to 2,800-fold higher than during normal replication. Importantly, this high rate of mutagenesis was observed not only close to the DSB where BIR is less stable, but also far from the DSB where the BIR replication fork is fast and stabilized. We established that polymerase proofreading and mismatch repair correct BIR errors. Also, dNTP levels were elevated during BIR, and this contributed to BIR-related mutagenesis. We propose that a high level of DNA polymerase errors that is not fully compensated by error-correction mechanisms is largely responsible for mutagenesis during BIR, with Pol δ generating many of the mutagenic errors. We further postulate that activation of BIR in eukaryotic cells may significantly contribute to accumulation of mutations that fuel cancer and evolution.

Accurate transmission of genetic information requires the precise replication of parental DNA. Mutations (which can be beneficial or deleterious) arise from errors that remain uncorrected. DNA replication occurs during S-phase of the cell cycle and is extremely accurate due to highly selective DNA polymerases coupled with effective error-correction mechanisms. In contrast, DNA synthesis associated with short-patch DNA repair is often error-prone. Break-induced replication (BIR) presents an interesting case of large-scale DNA duplication that occurs in the context of DNA repair. In this study we employed a yeast-based system to investigate the level of mutagenesis associated with BIR compared to mutagenesis during normal DNA replication. We report that frameshifts, which are the most deleterious kind of point mutation, are much more frequent during BIR than during normal DNA replication. Surprisingly, we observed that the majority of mutations associated with BIR were created by polymerases responsible for normal DNA replication, which are assumed to be highly precise. Overall, we propose that BIR is a novel source of mutagenesis that may contribute to disease genesis and evolution.

Energetic coupling between clustered lesions modulated by intervening triplet repeat bulge loops: allosteric implications for DNA repair and triplet repeat expansion

Clusters of closely spaced oxidative DNA lesions present challenges to the cellular repair machinery. When located in opposing strands, base excision repair (BER) of such lesions can lead to double strand DNA breaks (DSB). Activation of BER and DSB repair pathways has been implicated in inducing enhanced expansion of triplet repeat sequences. We show here that energy coupling between distal lesions (8oxodG and/or abasic sites) in opposing DNA strands can be modulated by a triplet repeat bulge loop located between the lesion sites. We find this modulation to be dependent on the identity of the lesions (8oxodG vs. abasic site) and the positions of the lesions (upstream vs. downstream) relative to the intervening bulge loop domain. We discuss how such bulge loop-mediated lesion crosstalk might influence repair processes, while favoring DNA expansion, the genotype of triplet repeat diseases.

Dynamic evolution of megasatellites in yeasts

Megasatellites are a new family of long tandem repeats, recently discovered in the yeast Candida glabrata. Compared to shorter tandem repeats, such as minisatellites, megasatellite motifs range in size from 135 to more than 300 bp, and allow calculation of evolutionary distances between individual motifs. Using divergence based on nucleotide substitutions among similar motifs, we determined the smallest distance between two motifs, allowing their subsequent clustering. Motifs belonging to the same cluster are recurrently found in different megasatellites located on different chromosomes, showing transfer of genetic information between megasatellites. In comparison, evolution of the few similar tandem repeats in Saccharomyces cerevisiae FLO genes mainly involves subtelomeric homologous recombination. We estimated selective constraints acting on megasatellite motifs and their host genes, and found that motifs are under strong purifying selection. Surprisingly, motifs inserted within pseudogenes are also under purifying selection, whereas the pseudogenes themselves evolve neutrally. We propose that megasatellite motifs propagate by a combination of three different molecular mechanisms: (i) gene duplication, (ii) ectopic homologous recombination and (iii) transfer of motifs from one megasatellite to another one. These mechanisms actively cooperate to create new megasatellites, that may play an important role in the adaptation of Candida glabrata to its human host.

Comparative genomics and molecular dynamics of DNA repeats in eukaryotes

Repeated elements can be widely abundant in eukaryotic genomes, composing more than 50% of the human genome, for example. It is possible to classify repeated sequences into two large families, "tandem repeats" and "dispersed repeats." Each of these two families can be itself divided into subfamilies. Dispersed repeats contain transposons, tRNA genes, and gene paralogues, whereas tandem repeats contain gene tandems, ribosomal DNA repeat arrays, and satellite DNA, itself subdivided into satellites, minisatellites, and microsatellites. Remarkably, the molecular mechanisms that create and propagate dispersed and tandem repeats are specific to each class and usually do not overlap. In the present review, we have chosen in the first section to describe the nature and distribution of dispersed and tandem repeats in eukaryotic genomes in the light of complete (or nearly complete) available genome sequences. In the second part, we focus on the molecular mechanisms responsible for the fast evolution of two specific classes of tandem repeats: minisatellites and microsatellites. Given that a growing number of human neurological disorders involve the expansion of a particular class of microsatellites, called trinucleotide repeats, a large part of the recent experimental work on microsatellites has focused on these particular repeats, and thus we also review the current knowledge in this area. Finally, we propose a unified definition for mini- and microsatellites that takes into account their biological properties and try to point out new directions that should be explored in a near future on our road to understanding the genetics of repeated sequences.

Coding tandem repeats generate diversity in Aspergillus fumigatus genes

Genes containing multiple coding mini- and microsatellite repeats are highly dynamic components of genomes. Frequent recombination events within these tandem repeats lead to changes in repeat numbers, which in turn alters the amino acid sequence of the corresponding protein. In bacteria and yeasts, the expansion of such coding repeats in cell wall proteins is associated with alterations in immunogenicity, adhesion, and pathogenesis. We hypothesized that identification of repeat-containing putative cell wall proteins in the human pathogen Aspergillus fumigatus may reveal novel pathogenesis-related elements. Here, we report that the genome of A. fumigatus contains as many as 292 genes with internal repeats. Fourteen of 30 selected genes showed size variation of their repeat-containing regions among 11 clinical A. fumigatus isolates. Four of these genes, Afu3g08990, Afu2g05150 (MP-2), Afu4g09600, and Afu6g14090, encode putative cell wall proteins containing a leader sequence and a glycosylphosphatidylinositol anchor motif. All four genes are expressed and produce variable-size mRNA encoding a discrete number of repeat amino acid units. Their expression was altered during development and in response to cell wall-disrupting agents. Deletion of one of these genes, Afu3g08990, resulted in a phenotype characterized by rapid conidial germination and reduced adherence to extracellular matrix suggestive of an alteration in cell wall characteristics. The Afu3g08990 protein was localized to the cell walls of dormant and germinating conidia. Our findings suggest that a subset of the A. fumigatus cell surface proteins may be hypervariable due to recombination events in their internal tandem repeats. This variation may provide the functional diversity in cell surface antigens which allows rapid adaptation to the environment and/or elusion of the host immune system.

Complex minisatellite rearrangements generated in the total or partial absence of Rad27/hFEN1 activity occur in a single generation and are Rad51 and Rad52 dependent

Genomes contain tandem repeat blocks that are at risk of expansion or contraction. The mechanisms of destabilization of the human minisatellite CEB1 (arrays of 36- to 43-bp repeats) were investigated in a previously developed model system, in which CEB1-0.6 (14 repeats) and CEB1-1.8 (42 repeats) alleles were inserted into the genome of Saccharomyces cerevisiae. As in human cells, CEB1 is stable in mitotically growing yeast cells but is frequently rearranged in the absence of the Rad27/hFEN1 protein involved in Okazaki fragments maturation. To gain insight into this mode of destabilization, the CEB1-1.8 and CEB1-0.6 human alleles and 47 rearrangements derived from a CEB1-1.8 progenitor in rad27Delta cells were sequenced. A high degree of polymorphism of CEB1 internal repeats was observed, attesting to a large variety of homology-driven rearrangements. Simple deletion, double deletion, and highly complex events were observed. Pedigree analysis showed that all rearrangements, even the most complex, occurred in a single generation and were inherited equally by mother and daughter cells. Finally, the rearrangement frequency was found to increase with array size, and partial complementation of the rad27Delta mutation by hFEN1 demonstrated that the production of novel CEB1 alleles is Rad52 and Rad51 dependent. Instability can be explained by an accumulation of unresolved flap structures during replication, leading to the formation of recombinogenic lesions and faulty repair, best understood by homology-dependent synthesis-strand displacement and annealing.

Evidence that Protein Length Expansion and Contraction Is Partly Due to Mutational Events in Premeiotic Cells

Studies on the rate of evolution of proteins typically concentrate on rates of change of orthologous amino acids rather than on changes in size (i.e., generation of nonorthologous domains). Recent work has focused attention on Ser/Thr-rich regions in yeast as these tend to undergo size changes rapidly, with size polymorphisms commonly being found, especially in proteins with cell-surface localization. The underlying mechanism generating the indels is presently unclear though, due to a lack of correlation with the location of meiotic double-strand breaks, it has, by exclusion, been conjectured to be replication slippage. Here we provide new evidence to support this possibility. Notably, we show that Ser/Thr-rich repeat regions are more generally associated with the location of Mre11p in premeiotic cells. This is to be expected if the repeats were produced by mutational events in mitotic cells possibly through replication slippage.

Molecular evolution of minisatellites in hemiascomycetous yeasts

Minisatellites are DNA tandem repeats exhibiting size polymorphism among individuals of a population. This polymorphism is generated by two different mechanisms, both in human and yeast cells, "replication slippage" during S-phase DNA synthesis and "repair slippage" associated to meiotic gene conversion. The Saccharomyces cerevisiae genome contains numerous natural minisatellites. They are located on all chromosomes without any obvious distribution bias. Minisatellites found in protein-coding genes have longer repeat units and on the average more repeat units than minisatellites in noncoding regions. They show an excess of cytosines on the coding strand, as compared to guanines (negative GC skew). They are always multiples of three, encode serine- and threonine-rich amino acid repeats, and are found preferably within genes encoding cell wall proteins, suggesting that they are positively selected in this particular class of genes. Genome-wide, there is no statistically significant association between minisatellites and meiotic recombination hot spots. In addition, minisatellites that are located in the vicinity of a meiotic hot spot are not more polymorphic than minisatellites located far from any hot spot. This suggests that minisatellites, in S. cerevisiae, evolve probably by strand slippage during replication or mitotic recombination. Finally, evolution of minisatellites among hemiascomycetous yeasts shows that even though many minisatellite-containing genes are conserved, most of the time the minisatellite itself is not conserved. The diversity of minisatellite sequences found in orthologous genes of different species suggests that minisatellites are differentially acquired and lost during evolution of hemiascomycetous yeasts at a pace faster than the genes containing them.

Enhanced oligonucleotide-directed gene targeting in mammalian cells following treatment with DNA damaging agents

Targeted gene repair, a form of oligonucleotide-directed mutagenesis, employs end-modified single-stranded DNA oligonucleotides to mediate single-base changes in chromosomal DNA. In this work, we use a specific 72-mer to direct the repair of a mutated eGFP gene stably integrated in the genome of DLD-1 cells. Corrected cells express eGFP that can be identified and quantitated by FACS. The repair of this mutant gene is dependent on the presence of a specifically designed oligonucleotide and the frequency with which the mutation is reversed is affected by the induction of DNA damage. We used hydroxyurea, VP16 (etoposide), and thymidine to modulate the rate of DNA replication through the stalling of the replication forks or the introduction of lesions. Addition of hydroxyurea or VP16 before the electroporation of the oligonucleotide, results in an accumulation of double-strand breaks (DSB) whose repair is facilitated by either nonhomologous end joining (NHEJ) or homologous recombination (HR). The addition of thymidine results in DNA damage within replication forks, damage that is repaired through the process of homologous recombination. Our data suggest that gene repair activity is elevated when DNA damage induces or activates the homologous recombination pathway.

Evolution of the regular zone of histone H1 in fabaceae plants

An analysis of the histone H1 subtype, H1-1, in eight legumes belonging to four genera of the tribe Vicieae (Pisum, Lathyrus, Lens, and Vicia), revealed an extended region consisting of the tandemly repeated AKPAAK motifs. We named this region the Regular zone (RZ). The AKPAAK motifs are organized into two blocks separated by a short (two or six amino acids) intervening sequence (IS). The distal block contains six AKPAAK motifs, while the number of repeats in the proximal block varies from six in V. faba to seven in the other species. In V. hirsuta, the first two repeated units of the proximal block are octapeptides AKAKPAAK. The apparent rate of synonymous substitutions in the blocks of RZ is much higher than in the rest of the gene. This can be explained by repeat shuffling within each block. In the C-domain of the orthologous H1 subtype from Medicago truncatula (tribe Trifolieae), a region corresponding to the RZ of Vicieae species was found. It also consists of two blocks of AKPAAK motifs (four and three repeats in the proximal and distal blocks, respectively). These blocks are separated by a 20-amino acid IS. The first 20 amino acids of the Medicago RZ are not part of AKPAAK repeats. We hypothesise that the RZ has most probably evolved as a result of an expansion of AKPAAK repeats from two separate sites in the C-domain. This process started tens of millions of years ago and was most likely directed by positive selection.

Microhomology-dependent end joining and repair of transposon-induced DNA hairpins by host factors in Saccharomyces cerevisiae

The maize, cut-and-paste transposon Ac/Ds is mobile in Saccharomyces cerevisiae, and DNA sequences of repair products provide strong genetic evidence that hairpin intermediates form in host DNA during this transposition, similar to those formed for V(D)J coding joints in vertebrates. Both DNA strands must be broken for Ac/Ds to excise, suggesting that double-strand break (DSB) repair pathways should be involved in repair of excision sites. In the absence of homologous template, as expected, Ac excisions are repaired by nonhomologous end joining (NHEJ) that can involve microhomologies close to the broken ends. However, unlike repair of endonuclease-induced DSBs, repair of Ac excisions in the presence of homologous template occurs by gene conversion only about half the time, the remainder being NHEJ events. Analysis of transposition in mutant yeast suggests roles for the Mre11/Rad50 complex, SAE2, NEJ1, and the Ku complex in repair of excision sites. Separation-of-function alleles of MRE11 suggest that its endonuclease function is more important in this repair than either its exonuclease or Rad50-binding properties. In addition, the interstrand cross-link repair gene PSO2 plays a role in end joining hairpin ends that is not seen in repair of linearized plasmids and may be involved in positioning transposase cleavage at the transposon ends.

Prion protein gene polymorphisms in Saccharomyces cerevisiae

The yeast Saccharomyces cerevisiae genome encodes several proteins that, in laboratory strains, can take up a stable, transmissible prion form. In each case, this requires the Asn/Gln-rich prion-forming domain (PrD) of the protein to be intact. In order to further understand the evolutionary significance of this unusual property, we have examined four different prion genes and their corresponding PrDs, from a number of naturally occurring strains of S. cerevisiae. In 4 of the 16 strains studied we identified a new allele of the SUP35 gene (SUP35delta19) that contains a 19-amino-acid deletion within the N-terminal PrD, a deletion that eliminates the prion property of Sup35p. In these strains a second prion gene, RNQ1, was found to be highly polymorphic, with eight different RNQ1 alleles detected in the six diploid strains studied. In contrast, for one other prion gene (URE2) and the sequence of the NEW1 gene encoding a PrD, no significant degree of DNA polymorphism was detected. Analysis of the naturally occurring alleles of RNQ1 and SUP35 indicated that the various polymorphisms identified were associated with DNA tandem repeats (6, 12, 33, 42 or 57 bp) within the coding sequences. The expansion and contraction of DNA repeats within the RNQ1 gene may provide an evolutionary mechanism that can ensure rapid change between the [PRION+] and [prion-] states.

Hypermutable minisatellites, a human affair?

Minisatellites are a class of highly polymorphic GC-rich tandem repeats. They include some of the most variable loci in the human genome, with mutation rates ranging from 0.5% to >20% per generation. Structurally, they consist of 10- to 100-bp intermingled variant repeats, making them ideal tools for dissecting mechanisms of instability at tandem repeats. Distinct mutation processes generate rare intra-allelic somatic events and frequent complex conversion-like germline mutations in these repeats. Furthermore, turnover of repeats at human minisatellites is controlled by intense recombinational activity in DNA flanking the repeat array. Surprisingly, whereas other mammalian genomes possess minisatellite-like sequences, hypermutable loci have not been identified that suggest human-specific turnover processes at minisatellite arrays. Attempts to transfer minisatellite germline instability to the mouse have failed. However, yeast models are now revealing valuable information regarding the mechanisms regulating instability at these tandem repeats. Finally, minisatellites and tandem repeats provide exquisitely sensitive molecular tools to detect genomic insults such as ionizing radiation exposure. Surprisingly, by a mechanism that remains elusive, there are transgenerational increases in minisatellite instability.

Strand invasion and DNA synthesis from the two 3' ends of a double-strand break in Mammalian cells

Analysis of the crossover products recovered following transformation of mammalian cells with a sequence insertion ("ends-in") gene-targeting vector revealed a novel class of recombinant. In this class of recombinants, a single vector copy has integrated into an ectopic genomic position, leaving the structure of the cognate chromosomal locus unaltered. Thus, in this respect, the recombinants resemble simple cases of random vector integration. However, the important difference is that the two paired 3' vector ends have acquired endogenous, chromosomal sequences flanking both sides of the vector-borne double-strand break (DSB). In some cases, copying was extensive, extending >16 kb into nonhomologous flanking DNA. The results suggest that mammalian homologous recombination events can involve strand invasion and DNA synthesis by both 3' ends of the DSB. These DNA interactions are a central, predicted feature of the DSBR model of recombination.

RAD1 controls the meiotic expansion of the human HRAS1 minisatellite in Saccharomyces cerevisiae

Minisatellite DNA is repetitive DNA with a repeat unit length from 15 to 100 bp. While stable during mitosis, it destabilizes during meiosis, altering both in length and in sequence composition. The basis for this instability is unknown. To investigate the factors controlling minisatellite stability, a minisatellite sequence 3' of the human HRAS1 gene was introduced into the Saccharomyces cerevisiae genome, replacing the wild-type HIS4 promoter. The minisatellite tract exhibited the same phenotypes in yeast that it exhibited in mammalian systems. The insertion stimulated transcription of the HIS4 gene; mRNA production was detected at levels above those seen with the wild-type promoter. The insertion stimulated meiotic recombination and created a hot spot for initiation of double-strand breaks during meiosis in the regions immediately flanking the repetitive DNA. The tract length altered at a high frequency during meiosis, and both expansions and contractions in length were detected. Tract expansion, but not contraction, was controlled by the product of the RAD1 gene. RAD1 is the first gene identified that controls specifically the expansion of minisatellite tracts. A model for tract length alteration based on these results is presented.

Break-induced replication: a review and an example in budding yeast

Break-induced replication (BIR) is a nonreciprocal recombination-dependent replication process that is an effective mechanism to repair a broken chromosome. We review key roles played by BIR in maintaining genome integrity, including restarting DNA replication at broken replication forks and maintaining telomeres in the absence of telomerase. Previous studies suggested that gene targeting does not occur by simple crossings-over between ends of the linearized transforming fragment and the target chromosome, but involves extensive new DNA synthesis resembling BIR. We examined gene targeting in Saccharomyces cerevisiae where only one end of the transformed DNA has homology to chromosomal sequences. Linearized, centromere-containing plasmid DNA with the 5' end of the LEU2 gene at one end was transformed into a strain in which the 5' end of LEU2 was replaced by ADE1, preventing simple homologous gene replacement to become Leu2(+). Ade1(+) Leu2(+) transformants were recovered in which the entire LEU2 gene and as much as 7 kb of additional sequences were found on the plasmid, joined by microhomologies characteristic of nonhomologous end-joining (NHEJ). In other experiments, cells were transformed with DNA fragments lacking an ARS and homologous to only 50 bp of ADE2 added to the ends of a URA3 gene. Autonomously replicating circles were recovered, containing URA3 and as much as 8 kb of ADE2-adjacent sequences, including a nearby ARS, copied from chromosomal DNA. Thus, the end of a linearized DNA fragment can initiate new DNA synthesis by BIR in which the newly synthesized DNA is displaced and subsequently forms circles by NHEJ.