A Computer Simulation Study of Vntr Population Genetics: Constrained Recombination Rules Out the Infinite Alleles Model

Extensive allelic diversity in variable numbers of tandem repeats (VNTRs) has been discovered in the human genome. For population genetic studies of VNTRs, such as forensic applications, it is important to know whether a neutral mutation-drift balance of VNTR polymorphism can be represented by the infinite alleles model. The assumption of the infinite alleles model that each new mutant is unique is very likely to be violated by unequal sister chromatid exchange (USCE), the primary process believed to generate VNTR mutants. We show that increasing both mutation rates and misalignment constraint for intrachromosomal recombination in a computer simulation model reduces simulated VNTR diversity below the expectations of the infinite alleles model. Maximal constraint, represented as slippage of single repeats, reduces simulated VNTR diversity to levels expected from the stepwise mutation model. Although misalignment rule is the more important variable, mutation rate also has an effect. At moderate rates of USCE, simulated VNTR diversity fluctuates around infinite alleles expectation. However, if rates of USCE are high, as for hypervariable VNTRs, simulated VNTR diversity is consistently lower than predicted by the infinite alleles model. This has been observed for many VNTRs and accounted for by technical problems in distinguishing alleles of neighboring size classes. We use sampling theory to confirm the intrinsically poor fit to the infinite alleles model of both simulated VNTR diversity and observed VNTR polymorphisms sampled from two Papua New Guinean populations.     

Genome-wide prediction of human VNTRs

Polymorphic minisatellites, also known as variable number of tandem repeats (VNTRs), are tandem repeat regions that show variation in the number of repeat units among chromosomes in a population. Currently, there are no general methods for predicting which minisatellites have a high probability of being polymorphic, given their sequence characteristics. An earlier approach has focused on potentially highly polymorphic and hypervariable minisatellites, which make up only a small fraction of all minisatellites in the human genome. We have developed a model, based on available minisatellite and VNTR sequence data, that predicts the probability that a minisatellite (unit size > or = 6 bp) identified by the computer program Tandem Repeats Finder is polymorphic (VNTR). According to the model, minisatellites with high copy number and high degree of sequence similarity are most likely to be VNTRs. This approach was used to scan the draft sequence of the human genome for VNTRs. A total of 157,549 minisatellite repeats were found, of which 29,224 are predicted to be VNTRs. Contrary to previous results, VNTRs appear to be widespread and abundant throughout the human genome, with an estimated density of 9.1 VNTRs/Mb.     

Slipped-strand mispairing at noncontiguous repeats in Poecilia reticulata: a model for minisatellite birth

The standard slipped-strand mispairing (SSM) model for the formation of variable number tandem repeats (VNTRs) proposes that a few tandem repeats, produced by chance mutations, provide the "raw material" for VNTR expansion. However, this model is unlikely to explain the formation of VNTRs with long motifs (e.g., minisatellites), because the likelihood of a tandem repeat forming by chance decreases rapidly as the length of the repeat motif increases. Phylogenetic reconstruction of the birth of a mitochondrial (mt) DNA minisatellite in guppies suggests that VNTRs with long motifs can form as a consequence of SSM at noncontiguous repeats. VNTRs formed in this manner have motifs longer than the noncontiguous repeat originally formed by chance and are flanked by one unit of the original, noncontiguous repeat. SSM at noncontiguous repeats can therefore explain the birth of VNTRs with long motifs and the "imperfect" or "short direct" repeats frequently observed adjacent to both mtDNA and nuclear VNTRs.     

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.     

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.     

Mechanisms of human minisatellite mutation in yeast

Minisatellites are tandem repeat loci, with repeat units ranging in size from 5 bp to 100 bp. The total lengths of repeat arrays vary from about 0.5 kb to 30 kb, and excessive variability in allele length at human minisatellite loci is the result of germline-specific complex recombination events generating new length alleles. Minisatellite alleles also mutate to new lengths in somatic cells, but this occurs at a much lower rate than in the germline. Since recombination is involved in minisatellite mutation, the yeast Saccharomyces cerevisiae is a suitable model organism that has been employed to further dissect the molecular basis of mutation events at human minisatellites. These studies have shown that the mutational behaviour of a minisatellite in meiosis is not determined by the intrinsic properties of the repeat array, but are highly dependent on the position of the minisatellite in the genome. The processes for minisatellite mutation in yeast and humans are identical in the sense that mutation is indeed driven by meiotic recombination, but differ with regard to the types of structural changes that are generated by the recombination events. Tetrad analyses showed that inter-allelic transfers of repeats occur by conversion and not crossing over, and that several chromatids can be involved in successive recombination events in one meiosis, resulting in mutant alleles in several spores. It has been demonstrated that the genes SPO11 and RAD50, involved in the initiation of recombination events, are required for human minisatellite mutation in yeast meiosis. Intrinsic properties of the repeat array appear to determine the stability of human minisatellites in yeast mitosis, since mitotic mutation rates in yeast are highly variable between minisatellites. The repair genes RAD27 and DNA2 stabilise human minisatellites in yeast mitosis, while RAD5 has no effect on mitotic stability. MSH2 depresses human minisatellite frequency in meiotic cells of yeast.     

In Depth Characterization of Repetitive DNA in 23 Plant Genomes Reveals Sources of Genome Size Variation in the Legume Tribe Fabeae

The differential accumulation and elimination of repetitive DNA are key drivers of genome size variation in flowering plants, yet there have been few studies which have analysed how different types of repeats in related species contribute to genome size evolution within a phylogenetic context. This question is addressed here by conducting large-scale comparative analysis of repeats in 23 species from four genera of the monophyletic legume tribe Fabeae, representing a 7.6-fold variation in genome size. Phylogenetic analysis and genome size reconstruction revealed that this diversity arose from genome size expansions and contractions in different lineages during the evolution of Fabeae. Employing a combination of low-pass genome sequencing with novel bioinformatic approaches resulted in identification and quantification of repeats making up 55–83% of the investigated genomes. In turn, this enabled an analysis of how each major repeat type contributed to the genome size variation encountered. Differential accumulation of repetitive DNA was found to account for 85% of the genome size differences between the species, and most (57%) of this variation was found to be driven by a single lineage of Ty3/gypsy LTR-retrotransposons, the Ogre elements. Although the amounts of several other lineages of LTR-retrotransposons and the total amount of satellite DNA were also positively correlated with genome size, their contributions to genome size variation were much smaller (up to 6%). Repeat analysis within a phylogenetic framework also revealed profound differences in the extent of sequence conservation between different repeat types across Fabeae. In addition to these findings, the study has provided a proof of concept for the approach combining recent developments in sequencing and bioinformatics to perform comparative analyses of repetitive DNAs in a large number of non-model species without the need to assemble their genomes.     

Analysis of VNTRs in the solute carrier family 6, member 18 (SLC6A18) and lack of association with hypertension

We report here the distribution of VNTRs (variable number of tandem repeats; minisatellites) and polymorphic analysis of SLC6A18, which is a member of the SLC6 Na(+)- and Cl(-)-dependent neurotransmitter transporter family. In this study, DNA was obtained from 300 unrelated individuals and 205 patients with essential hypertension (EH). We then analyzed the VNTRs in the genomic DNA by searching for minisatellites of SLC6A18 using the Tandem Repeat Finder program. Eight novel VNTRs were identified: five of which were polymorphic minisatellites (SLC6A18-MS1, -MS2, -MS4, -MS5, and -MS6) and three of which were monomorphic minisatellites (SLC6A18-MS3, -MS7, and -MS8). Next, we investigated the relationship between EH and four of the polymorphic minisatellites (SLC6A18-MS1, -MS2, -MS4, and -MS6). We excluded SLC6A18-MS5 from the common/rare allele analysis, because most individuals were heterozygous and hypervariable for this locus. There were no significant differences observed in the overall distribution of these minisatellites, which indicates that these polymorphisms are not responsible for EH susceptibility in the Korean population. A segregation analysis of the minisatellites in SLC6A18 was then conducted by analyzing genomic DNA obtained from two generations of five families and from three generations of two families. The five polymorphic minisatellites were transmitted through meiosis following Mendelian inheritance, which suggests that polymorphic minisatellites could be useful markers for paternity mapping and DNA fingerprinting. In summary, we discovered five novel VNTR polymorphisms in SLC6A18; however, these variations were not associated with EH.     

Large mitochondrial repeats multiplied during the polymerase chain reaction

The number of repeats in repetitive DNA like micro‐ and minisatellites is often determined by polymerase chain reaction (PCR). When we counted repeats in an array of mitochondrial repeats in the cattle tick (Boophilus microplus) we found that the number of repeats increased during PCR. Multiplication of the repeats was independent of the primers used to amplify the region, the PCR annealing temperature and the length of the PCR product. The use of PCR to determine the number of repeats in arrays needs to be reassessed. For long repeats, a subset of samples should always be analysed by Southern blot hybridization to confirm the PCR results.     

Influences of array size and homogeneity on minisatellite mutation.

Unstable minisatellites display high frequencies of spontaneous gain and loss of repeats in the human germline. Most length changes arise through complex recombination events including intra-allelic duplications/deletions and inter-allelic transfers of repeats. Definition of the factors modulating instability requires both measurement of mutation rate and detailed analysis of mutant structures at the level of individual alleles. We have measured mutation rates in sperm for a wide range of alleles of the highly unstable human minisatellite CEB1. Instability varies by three orders of magnitude between alleles and increases steadily with the size of the tandem array. Structural analysis of mutant molecules derived from six alleles revealed that it is the rate of intra-allelic rearrangements which increases with array size and that intra-allelic duplication events tend to cluster within homogeneous segments of alleles; both phenomena resemble features of trinucleotide repeat instability. In contrast, inter-allelic transfers occur at a fairly constant rate, irrespective of array length, and show a mild polarity towards one end of the minisatellite, suggesting the possible influence of flanking DNA on these conversion-like events.     

A census of protein repeats.

In this study, we analyzed all known protein sequences for repeating amino acid segments. Although duplicated sequence segments occur in 14 % of all proteins, eukaryotic proteins are three times more likely to have internal repeats than prokaryotic proteins. After clustering the repetitive sequence segments into families, we find repeats from eukaryotic proteins have little similarity with prokaryotic repeats, suggesting most repeats arose after the prokaryotic and eukaryotic lineages diverged. Consequently, protein classes with the highest incidence of repetitive sequences perform functions unique to eukaryotes. The frequency distribution of the repeating units shows only weak length dependence, implicating recombination rather than duplex melting or DNA hairpin formation as the limiting mechanism underlying repeat formation. The mechanism favors additional repeats once an initial duplication has been incorporated. Finally, we show that repetitive sequences are favored that contain small and relatively water-soluble residues. We propose that error-prone repeat expansion allows repetitive proteins to evolve more quickly than non-repeat-containing proteins.     

The influence of sequence divergence between alleles of the human MS205 minisatellite incorporated into the yeast genome on length-mutation rates and lethal recombination events during meiosis

Certain minisatellites exhibit hypervariability with respect to the number of repeat units and, thus, allele length. Such polymorphism is generated by germline-specific recombinational events that occur at high frequencies and lead to the gain or loss of repeat units. In order to elucidate the molecular details of mutagenesis in minisatellites, we have integrated human minisatellites into the yeast genome in the vicinity of a hotspot for meiotic double-strand breaks (DSBs). Here, we describe the results of tetrad analyses of mutations in the human MS205 minisatellite in yeast strains heterozygous for alleles composed of 51 and 31 repeat units, as well as in a strain homozygous for the same 51 repeat unit allele. The length-mutation rate was twice as high in the heterozygous strain as in the homozygous strain, suggesting that sequence divergence between alleles enhances the generation of length mutations. In the case of heterozygotes, the frequency of length mutants resulting from inter-allelic exchange was significantly higher in tetrads with three viable spores than in tetrads with four viable spores, indicating that there is a higher probability for spore mortality in tetrads originating from meioses during which inter-allelic exchange of repeat units occurs. In an attempt to explain these findings, we propose a model for minisatellite mutation involving recombination, in which sequence divergence between alleles results in a heteroduplex containing numerous mismatches. We suggest that convergent mismatch-repair tracts in this heteroduplex give rise to a DSB that may be repaired by an additional round of recombination resulting in mutation of a third allele, or be lethal if such recombination fails. It appears probable that the formation of such additional mutants is the major explanation for the difference in meiotic length-mutation rates between the heterozygous and homozygous yeast strains, and that this phenomenon contributes to high germline length-mutation frequencies at minisatellites in humans.     

How do chromosomes attach to synaptonemal complexes?

Fluorochrome-labeled oligonucleotides (n = 44) corresponding to mouse genome repetitive sequences were hybridized in situ with pachytene nuclei of mouse spermatocytes. Signals of the repetitive sequences MaLR, MER, and (GT)₂₂ were found to be dispersed through chromatin, and signals of B1 repeats and minisatellites were mostly attached to synaptonemal complexes immunostained with anti-SYCP3 antibodies. These results suggest that B1 repeats and minisatellites are candidates for sequences anchoring chromatin to synaptonemal complexes.     

Efficient repair of HO-induced chromosomal breaks in Saccharomyces cerevisiae by recombination between flanking homologous sequences.

Novel recombinational repair of a site-specific double-strand break (DSB) in a yeast chromosome was investigated. When the recognition site for the HO endonuclease enzyme is embedded in nonyeast sequences and placed between two regions of homology, expression of HO endonuclease stimulates recombination between the homologous flanking regions to yield a deletion, the apparent product of an intrachromosomal exchange between direct repeats. This deletion-repair event is very efficient, thus preventing essentially all the potential lethality due to the persistence of a DSB. Interestingly, unlike previous studies involving spontaneous recombination between chromosomal repeats, the recombination events stimulated by HO-induced DSBs are accompanied by loss of the sequences separating the homologous regions greater than 99.5% of the time. Repair is dependent on the RAD52 gene. The deletion-repair event provides an in vivo assay for the sensitivity of any particular recognition site to HO cleavage. By taking advantage of a galactose-inducible HO gene, it has been possible to follow the kinetics of this event at the DNA level and to search for intermediates in this reaction. Deletion-repair requires approximately 45 min and is inhibited when cycloheximide is added after HO endonuclease cleavage.     

Structure of the Y Chromosomal Su(ste) Locus in Drosophila Melanogaster and Evidence for Localized Recombination among Repeats

The structure of the Suppressor of Stellate [Su(Ste)] locus on the Drosophila melanogaster Y chromosome was examined by restriction analysis of both native and cloned genomic DNA. The locus consists of short subarrays of tandem repeats separated by members of other moderately repeated families. Both size variants and restriction variants proved to be common. Most repeats fell into two size classes--2.8 and 2.5 kb--but other size variants were also observed. Restriction variants showed a strong tendency to cluster, both at the gross level where some variants were present in only one of three subintervals of the locus, and at the fine level, where repeats from the same phage clone were significantly more similar than repeats from different clones. Restriction variants were shared freely among repeats of different size classes; however, size variants appeared to be randomly distributed among phage clones. These data indicate that recombination among tandem Su(Ste) repeats occurs at much higher frequencies between close neighbors than distant ones. In addition, they suggest that gene conversion rather than sister chromatid exchange may be the primary recombinational mechanism for spreading variation among repeats at the Su(Ste) locus.     

Mutations at the human minisatellite MS32 integrated in yeast occur with high frequency in meiosis and involve complex recombination events

Minisatellites are composed of tandem repetitive DNA sequences and are present at many positions in the human genome. They frequently mutate to new length alleles in the germline, by complex and incompletely understood recombination mechanisms which may operate during meiosis. In several minisatellites the mutation events are restricted to one end of the repeat array, indicating a possible association with elements that act in cis. Mutant alleles do not show exchange of flanking regions. To construct a model system suitable for further investigations of the mutation process, we have integrated the human minisatellite MS32, flanked by synthetic markers, in the vicinity of a meiotic recombination hot spot upstream of the LEU2 locus in the yeast Saccharomyces cerevisiae. Here we provide direct evidence for a meiotic origin of MS32 mutations. Mutation events were polarised towards both ends of the minisatellite and varied from simple duplications and deletions to complex intra- and interallelic events. Interallelic events were frequently accompanied by exchange of regions flanking the minisatellite. The results also support the notion that cis-acting elements are involved in the mutational process. The fact that MS32 mutant structures are similar in yeast and human shows that meiotic recombination plays a crucial role in both organisms and emphasises the usefulness of yeast strains harbouring minisatellites as a model system for the study of minisatellite mutation. Received: 1 March 1997 / Accepted: 16 May 1997     

Effects of recombination rate on human endogenous retrovirus fixation and persistence

Endogenous retroviruses (ERVs) result from germ line infections by exogenous retroviruses. They can proliferate within the genome of their host species until they are either inactivated by mutation or removed by recombinational deletion. ERVs belong to a diverse group of mobile genetic elements collectively termed transposable elements (TEs). Numerous studies have attempted to elucidate the factors determining the genomic distribution and persistence of TEs. Here we show that, within humans, gene density and not recombination rate correlates with fixation of endogenous retroviruses, whereas the local recombination rate determines their persistence in a full-length state. Recombination does not appear to influence fixation either via the ectopic exchange model or by indirect models based on the efficacy of selection. We propose a model linking rates of meiotic recombination to the probability of recombinational deletion to explain the effect of recombination rate on persistence. Chromosomes 19 and Y are exceptions, possessing more elements than other regions, and we suggest this is due to low gene density and elevated rates of human ERV integration in males for chromosome Y and segmental duplication for chromosome 19.     

Persistence of Tandem Arrays: Implications for Satellite and Simple-Sequence Dnas

Recombination processes acting on tandem arrays are suggested here to have probable intrinsic biases, producing an expected net decrease in array size following each event, in contrast to previous models which assume no net change in array size. We examine the implications of this by modeling copy number dynamics in a tandem array under the joint interactions of sister-strand unequal crossing over (rate gamma per generation per copy) and intrastrand recombination resulting in deletion (rate epsilon per generation per copy). Assuming no gene amplification or selection, the expected mean persistence time of an array starting with z excess copies (i.e., array size z + 1) is z(1 + gamma/epsilon) recombinational events. Nontrivial equilibrium distributions of array sizes exist when gene amplification or certain forms of selection are considered. We characterize the equilibrium distribution for both a simple model of gene amplification and under the assumption that selection imposes a minimal array size, n. For the latter case, n + 1/alpha is an upper bound for mean array size under fairly general conditions, where alpha(= 2 epsilon/gamma) is the scaled deletion rate. Further, the distribution of excess copies over n is bounded above by a geometric distribution with parameter alpha/(1 + alpha). Tandem arrays are unlikely to be greatly expanded by unequal crossing over unless alpha much less than 1, implying that other mechanisms, such as gene amplification, are likely important in the evolution of large arrays. Thus unequal crossing over, by itself, is likely insufficient to account for satellite DNA.     

Sequences flanking the repeat arrays of human minisatellites: association with tandem and dispersed repeat elements.

We present DNA sequences flanking cloned hypervariable human minisatellites. In addition to providing confirmatory evidence that minisatellites cluster with other tandem repeats, these flanking sequences contain a high frequency of interspersed repetitive elements. These elements include a retroviral LTR-like sequence, from which one of the minisatellites appears to have expanded, and a recently described short interspersed repeat. We present our own findings concerning this element, in particular that those examples studied do not show significant evolutionary conservation, despite suggestions that the element may have a cis-acting function.