Homogenization of Tandemly Repeated Nucleotide Sequences by Distance-Dependent Nucleotide Sequence Conversion

Previous work revealed that recurrent mutations (= mutation occurring more than once) in the tandemly repeated arrays present in nontranscribed spacers (NTS) of ribosomal RNA genes (rDNA) are clustered, i.e., they most frequently occur in repeats with adjacent or alternate distribution. A possible explanation is that the likelihood of heteroduplex formation, a prerequisite of gene conversion, decreases with the distance between repeats. To test this possibility, evolution of an array of 11 initially homogeneous repeats was computer simulated using three models, two assuming that the likelihood of heteroduplex formation decreases with increasing distance between the repeats and one assuming that it is constant. Patterns of mutation distribution obtained in computer simulations were compared with the distribution of mutations found in the repeated arrays in the NTS of seven rDNA clones. The patterns of mutations generated by the models assuming that the likelihood of heteroduplex formation decreases as distance between the repeats increases agreed with the patterns observed in rDNA; the patterns generated by the model assuming that the likelihood is independent of distance between repeats disagreed with the patterns observed in the rDNA clones. The topology of the heteroduplex formed between DNA in adjacent repeats predicts that the most frequently occurring conversions in the NTS repeated arrays will be shorter than the length of the repeat. The topology of the heteroduplex also predicts that if the heteroduplex leads to crossing over a circular repeat is excised. It is speculated that the circle can transpose or can be amplified via rolling circle replication and subsequently transpose. It is also shown that homogenization of the NTS repeated arrays proceeds at different rates in different species.     

Analysis of ALS5 and ALS6 allelic variability in a geographically diverse collection of Candida albicans isolates

The Candida albicans ALS (agglutinin-like sequence) gene family encodes eight cell-surface glycoproteins, some of which function in adhesion to host surfaces. ALS genes have a central tandem repeat-encoding domain comprised entirely of head-to-tail copies of a conserved 108-bp sequence. The number of copies of the tandemly repeated sequence varies between C. albicans strains and often between alleles within the same strain. Because ALS alleles can encode different-sized proteins that may have different functional characteristics, defining the range of allelic variability is important. Genomic DNA from C. albicans strains representing the major genetic clades was PCR amplified to determine the number of tandemly repeated sequence copies within the ALS5 and ALS6 central domain. ALS5 alleles had 2-10 tandem repeat sequence copies (mean=4.82 copies) while ALS6 alleles had 2-8 copies (mean=4.00 copies). Despite this variability, tandem repeat copy number was stable in C. albicans strains passaged for 3000 generations. Prevalent alleles and allelic distributions varied among the clades for ALS5 and ALS6. Overall, ALS6 exhibited less variability than ALS5. ALS5 deletions can occur naturally in C. albicans via direct repeats flanking the ALS5 locus. Deletion of both ALS5 alleles was associated particularly with clades III and SA. ALS5 exhibited allelic polymorphisms in the coding region 5' of the tandem repeats; some alleles resembled ALS1, suggesting recombination between these contiguous loci. Natural deletion of ALS5 and the sequence variation within its coding region suggest relaxed selective pressure on this locus, and that Als5p function may be dispensable in C. albicans or redundant within the Als family.     

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.      

Digital Genotyping of Macrosatellites and Multicopy Genes Reveals Novel Biological Functions Associated with Copy Number Variation of Large Tandem Repeats

Tandem repeats are common in eukaryotic genomes, but due to difficulties in assaying them remain poorly studied. Here, we demonstrate the utility of Nanostring technology as a targeted approach to perform accurate measurement of tandem repeats even at extremely high copy number, and apply this technology to genotype 165 HapMap samples from three different populations and five species of non-human primates. We observed extreme variability in copy number of tandemly repeated genes, with many loci showing 5–10 fold variation in copy number among humans. Many of these loci show hallmarks of genome assembly errors, and the true copy number of many large tandem repeats is significantly under-represented even in the high quality ‘finished’ human reference assembly. Importantly, we demonstrate that most large tandem repeat variations are not tagged by nearby SNPs, and are therefore essentially invisible to SNP-based GWAS approaches. Using association analysis we identify many cis correlations of large tandem repeat variants with nearby gene expression and DNA methylation levels, indicating that variations of tandem repeat length are associated with functional effects on the local genomic environment. This includes an example where expansion of a macrosatellite repeat is associated with increased DNA methylation and suppression of nearby gene expression, suggesting a mechanism termed “repeat induced gene silencing”, which has previously been observed only in transgenic organisms. We also observed multiple signatures consistent with altered selective pressures at tandemly repeated loci, suggesting important biological functions. Our studies show that tandemly repeated loci represent a highly variable fraction of the genome that have been systematically ignored by most previous studies, copy number variation of which can exert functionally significant effects. We suggest that future studies of tandem repeat loci will lead to many novel insights into their role in modulating both genomic and phenotypic diversity.     

A Gene with Specific and Global Effects on Recombination of Sequences from Tandemly Repeated Genes in Saccharomyces Cerevisiae

The preservation of sequence homogeneity and copy number of tandemly repeated genes may require specific mechanisms or regulation of recombination. We have identified mutations that specifically affect recombination among natural repetitions in the yeast Saccharomyces cerevisiae. The rrm3 mutation stimulates mitotic recombination in the naturally occurring tandem repeats of the rDNA and copper chelatin (CUP1) genes. This mutation does not affect recombination of several other types of repeated genes tested including Ty elements, mating type information and duplications created by transformation. In addition to stimulating exchange among the multiple CUP1 repeats at their natural chromosomal location, rrm3 also increases recombination of a duplication of CUP1 units present at his4. This suggests that the RRM3 gene may encode a sequence-specific factor that contributes to a global suppression of mitotic exchange in sequences that can be maintained as tandem arrays.     

Chromosomal locations of major tRNA gene clusters of Xenopus laevis

In Xenopus laevis eight tRNA genes are located in a 3.18 kb tandemly repeated unit. There are 150 copies of the unit at a single locus near the long arm telomere of one of the acrocentric chromosomes in the 14–17 group. Two additional classes of tRNA gene-containing repeats have been isolated (defined by clones p3.1 and p3.2) that have structures related to that of the 3.18 kb unit. Using in situ hybridization at the electron microscopic level, the p3.2 repeats are found clustered at a single locus in the subtelomeric region on one of the submetacentric chromosomes, whereas the p3.1 repeats are clustered at a locus indistinguishable from that containing the 3.18 kb repeats. This suggests that these tDNA tandem repeats can diverge in sequence from each other without being at distantly separated loci.     

Complex structure of knobs and centromeric regions in maize chromosomes

The recovery of maize (Zea mays L.) chromosome addition lines of oat (Avena sativa L.) from oat x maize crosses enables us to analyze the structure and composition of individual maize chromosomes via the isolation and characterization of chromosome-specific cosmid clones. Restriction fragment fingerprinting, sequencing, and in situ hybridization were applied to discover a new family of knob associated tandem repeats, the TR1, which are capable of forming fold-back DNA segments, as well as a new family of centromeric tandem repeats, CentC. Analysis of knob and centromeric DNA segments revealed a complex organization in which blocks of tandemly arranged repeating units are interrupted by insertions of other repeated DNA sequences, mostly represented by individual full size copies of retrotransposable elements. There is an obvious preference for the integration/association of certain retrotransposable elements into knobs or centromere regions as well as for integration of retrotransposable elements into certain sites (hot spots) of the 180-bp repeat. DNA hybridization to a blot panel of eight individual maize chromosome addition lines revealed that CentC, TR1, and 180-bp tandem repeats are found in each of these maize chromosomes, but the copy number of each can vary significantly from about 100 to 25,000. In situ hybridization revealed variation among the maize chromosomes in the size of centromeric tandem repeats as well as in the size and composition of knob regions. It was found that knobs may be composed of either 180-bp or TR1, or both repeats, and in addition to large knobs these repeated elements may form micro clusters which are detectable only with the help of in situ hybridization. The association of the fold-back elements with knobs, knob polymorphism and complex structure suggest that maize knob may be consider as megatransposable elements. The discovery of the interspersion of retrotransposable elements among blocks of tandem repeats in maize and some other organisms suggests that this pattern may be basic to heterochromatin organization for eukaryotes.     

The 1.688 repetitive DNA of Drosophila: concerted evolution at different genomic scales and association with genes

Concerted evolution leading to homogenization of tandemly repeated DNA arrays is widespread and important for genome evolution. We investigated the range and nature of the process at chromosomal and array levels using the 1.688 tandem repeats of Drosophila melanogaster where large arrays are present in the heterochromatin of chromosomes 2, 3, and X, and short arrays are found in the euchromatin of the same chromosomes. Analysis of 326 euchromatic and heterochromatic repeats from 52 arrays showed that the homogenization of 1.688 repeats occurred differentially for distinct genomic regions, from euchromatin to heterochromatin and from local arrays to chromosomes. We further found that most euchromatic arrays are either close to, or are within introns of, genes. The short size of euchromatic arrays (one to five repeats) could be selectively constrained by their role as gene regulators, a situation similar to the so-called "tuning knobs."     

Tandem-repetitive noncoding DNA: forms and forces

A model of sequence-dependent, unequal crossing-over and gene amplification (slippage replication) has been stimulated in order to account for various structural features of tandemly repeated DNA sequences. It is shown that DNA whose sequence is not maintained by natural selection will exhibit repetitive patterns over a wide range of recombination rates as a result of the interaction of unequal crossing-over and slippage replication, processes that depend on sequence similarity. At high crossing-over frequencies, the nucleotide patterns generated in the simulations are simple and highly regular, with short, nearly identical sequences repeated in tandem. Decreasing recombination rates increase the tendency to longer and more-complex repeat units. Periodicities have been observed down to very low recombination rates (one or more orders of magnitude lower than mutation rate). At such low rates, most of the sequences contain repeats which have an extensive substructure and a high degree of heterogeneity among each other; often higher-order structures are superimposed on a tandem array. These results are compared with various structural properties of tandemly repeated DNAs known from eukaryotes, the spectrum ranging from simple-sequence DNAs, particularly the hypervariable mini-satellites, to the classical satellite DNAs, located in chromosomal regions of low recombination, e.g., heterochromatin.     

Repetitive satellite-like sequences are present within or upstream from 3 avian protein-coding genes.

Peculiar DNA sequences made up by the tandem repetition of a 5 bp unit have been identified within or upstream from three avian protein-coding genes. One sequence is located within an intron of the chicken "ovalbumin-X" gene with 5'-TCTCC-3' as basic repeat unit (36 repeats). Another sequence made of 27 repeats of a 5'-GGAAG-3' basic unit is found 2500 base pairs upstream from the promoter of the chicken ovotransferrin (conalbumin) gene. A related but different sequence is present in the corresponding region of the ovotransferrin gene in the pheasant, with 5'-GGAAA-3' as the basic unit (55 repeats). These three satellite-like elements are thus characterized by a total assymetry in base distribution, with purines restricted to one strand, and pyrimidines to the other. Two of the basic repeat units can be derived from the third one (GGAAA) by a single base pair change. These related sequences are found repeated in three avian genomes, at degrees which vary both with the sequence type and the genome type. Evolution of tandemly repeated sequences (including satellites) is in general studied by analysing randomly picked elements. The presence of conserved protein-coding regions neighbouring satellite-like sequences allow to follow their evolution at a single locus, as exemplified by the striking comparison of the pheasant and chicken sequences upstream from the ovotransferrin gene.     

Arrays of repetitive DNA elements in the largest chromosomes of Toxoplasma gondii.

A novel tandemly repeated DNA structure of Toxoplasma gondii that meets the requirements assigned for satellital DNA was characterized. A DNA fragment of 1002 bp contains two different elements of repetitive DNA families named ABGTg7 and ABGTg8.2. Both repeats are members of a more complex tandem structure where ABGTg7-like monomers can be arranged either as direct tandems or flanked by other related or non-related repeats. Pulse-field gel electrophoresis analysis showed that these repeats hybridize with the largest T. gondii chromosomes. Bal31 sensitivity assays indicated that these elements are located near the telomeres and along other regions too. Five genomic lambda phages were isolated and two different completed clusters of the repeated structure were analyzed.     

Biased amino acid composition in repeat regions of Plasmodium antigens.

Many malarial antigens contain extensive arrays of tandemly repeated short amino acid sequences, and much of the antibody response induced by malaria infections is directed against these repeats. Indeed, it has been hypothesized that these repeats function to elicit a relatively ineffective T-cell-independent antibody response by the host. In order to test this hypothesis, tandem repeats of Plasmodium species were examined for a bias in composition favoring amino acids likely to form epitopes for the antibody. The genome of Plasmodium is very A+T-rich, and nucleotide compositional bias will, in itself, lead to a high proportion of hydrophilic amino acids. When this bias was controlled for, Plasmodium antigens did not show a higher proportion of hydrophilic amino acids than expected, but there was a significant reduction in the proportion of hydrophobic amino acids in the repeats of the antigens. The amino acid composition of the repeats was thus strikingly different from those seen both in the remainder of the antigens and in a sample of Plasmodium falciparum housekeeping genes.     

Expansion and contraction of the DUP240 multigene family in Saccharomyces cerevisiae populations

The influence of duplicated sequences on chromosomal stability is poorly understood. To characterize chromosomal rearrangements involving duplicated sequences, we compared the organization of tandem repeats of the DUP240 gene family in 15 Saccharomyces cerevisiae strains of various origins. The DUP240 gene family consists of 10 members of unknown function in the reference strain S288C. Five DUP240 paralogs on chromosome I and two on chromosome VII are arranged as tandem repeats that are highly polymorphic in copy number and sequence. We characterized DNA sequences that are likely involved in homologous or nonhomologous recombination events and are responsible for intra- and interchromosomal rearrangements that cause the creation and disappearance of DUP240 paralogs. The tandemly repeated DUP240 genes seem to be privileged sites of gene birth and death.     

An algorithm for approximate tandem repeats.

A perfect single tandem repeat is defined as a nonempty string that can be divided into two identical substrings, e.g., abcabc. An approximate single tandem repeat is one in which the substrings are similar, but not identical, e.g., abcdaacd. In this paper we consider two criterions of similarity: the Hamming distance (k mismatches) and the edit distance (k differences). For a string S of length n and an integer k our algorithm reports all locally optimal approximate repeats, r = umacro ?, for which the Hamming distance of umacro and ? is at most k, in O(nk log (n/k)) time, or all those for which the edit distance of umacro and ? is at most k, in O(nk log k log (n/k)) time. This paper concentrates on a more general type of repeat called multiple tandem repeats. A multiple tandem repeat in a sequence S is a (periodic) substring r of S of the form r = u(a)u', where u is a prefix of r and u' is a prefix of u. An approximate multiple tandem repeat is a multiple repeat with errors; the repeated subsequences are similar but not identical. We precisely define approximate multiple repeats, and present an algorithm that finds all repeats that concur with our definition. The time complexity of the algorithm, when searching for repeats with up to k errors in a string S of length n, is O(nka log (n/k)) where a is the maximum number of periods in any reported repeat. We present some experimental results concerning the performance and sensitivity of our algorithm. The problem of finding repeats within a string is a computational problem with important applications in the field of molecular biology. Both exact and inexact repeats occur frequently in the genome, and certain repeats occurring in the genome are known to be related to diseases in the human.     

Tandem repeat markers as novel diagnostic tools for high resolution fingerprinting of <Emphasis Type="Italic">Wolbachia</Emphasis>

BACKGROUND: Strains of the endosymbiotic bacterium Wolbachia pipientis are extremely diverse both genotypically and in terms of their induced phenotypes in invertebrate hosts. Despite extensive molecular characterisation of Wolbachia diversity, little is known about the actual genomic diversity within or between closely related strains that group tightly on the basis of existing gene marker systems, including Multiple Locus Sequence Typing (MLST). There is an urgent need for higher resolution fingerprinting markers of Wolbachia for studies of population genetics, horizontal transmission and experimental evolution. RESULTS: The genome of the wMel Wolbachia strain that infects Drosophila melanogaster contains inter- and intragenic tandem repeats that may evolve through expansion or contraction. We identified hypervariable regions in wMel, including intergenic Variable Number Tandem Repeats (VNTRs), and genes encoding ankyrin (ANK) repeat domains. We amplified these markers from 14 related Wolbachia strains belonging to supergroup A and were successful in differentiating size polymorphic alleles. Because of their tandemly repeated structure and length polymorphism, the markers can be used in a PCR-diagnostic multilocus typing approach, analogous to the Multiple Locus VNTR Analysis (MLVA) established for many other bacteria and organisms. The isolated markers are highly specific for supergroup A and not informative for other supergroups. However, in silico analysis of completed genomes from other supergroups revealed the presence of tandem repeats that are variable and could therefore be useful for typing target strains. CONCLUSIONS: Wolbachia genomes contain inter- and intragenic tandem repeats that evolve through expansion or contraction. A selection of polymorphic tandem repeats is a novel and useful PCR diagnostic extension to the existing MLST typing system of Wolbachia, as it allows rapid and inexpensive high-throughput fingerprinting of closely related strains for which polymorphic markers were previously lacking.     

Comparative analyses of human single- and multilocus tandem repeats

Using the compiled human genome sequence, we systematically cataloged all tandem repeats with periods between 20 and 2000 bp and defined two subsets whose consensus sequences were found at either single-locus tandem repeats (slTRs) or multilocus tandem repeats (mlTRs). Parameters compiled for these subsets provide insights into mechanisms underlying the creation and evolution of tandem repeats. Both subsets of tandem repeats are nonrandomly distributed in the genome, being found at higher frequency at many but not all chromosome ends and internal clusters of mlTRs were also observed. Despite the integral role of recombination in the biology of tandem repeats, recombination hotspots colocalized only with shorter microsatellites and not the longer repeats examined here. An increased frequency of slTRs was observed near imprinted genes, consistent with a functional role, while both slTRs and mlTRs were found more frequently near genes implicated in triplet expansion diseases, suggesting a general instability of these regions. Using our collated parameters, we identified 2230 slTRs as candidates for highly informative molecular markers.     

A B-cell-specific nuclear protein that binds to DNA sites 5'' to immunoglobulin S alpha tandem repeats is regulated during differentiation.

Immunoglobulin heavy-chain switching is effected by recombination events between sites associated with tandemly repeated switch sequences located 5' to immunoglobulin heavy-chain genes. Using the band mobility shift assay, we have identified two distinct sites 5' to the alpha heavy-chain switch sequence with affinity for a single B-cell-specific DNA-binding protein, S alpha-BP. S alpha-BP was present in nuclear extracts from pre-B and B cells but was not detected in extracts from plasmacytomas, B-cell hybridomas, T-cell lymphomas, or a macrophage cell line. It was also not detectable in other nonlymphoid cells tested. Evidence suggests there are S alpha-BP-binding sites near other immunoglobulin switch sequences. As with the S alpha sites, these sites appear to be distinct from the consensus tandem repeats characteristic of immunoglobulin switch sequences. The possible functions of S alpha-BP on contacting its binding sites are discussed in the context of immunoglobulin heavy-chain switch recombination.     

The immunoglobulin heavy chain switch: structural features of gamma 1 recombinant switch regions

The immunoglobulin heavy chain isotype switch is mediated by a DNA rearrangement involving specific genomic segments referred to as switch regions. Switch regions are composed of tandemly repeated simple sequences. The role of the tandemly repeated structure of switch regions in the switch recombination process is not understood. We mapped eight recombination sites--six in the gamma 1 and two in the gamma 3 tandem arrays. In addition, we obtained molecular clones representing three of the six gamma 1 rearrangements, and determined the nucleotide sequences of the recombination sites in each. In general, the rearrangements are confined to the tandem repeat units, and are not clustered in a particular portion of either the gamma 3 or gamma 1 switch region. Nucleotide sequence analysis of one of the recombinant clones, gamma M35, reveals evidence for a successive switch event wherein a recombination between S mu and S gamma 3 was followed by recombination 57 bp downstream with S gamma 1. gamma 1 sequence data from the molecular clones we obtained, together with similar data from other investigators regarding the gamma 1, gamma 2b, and gamma 2a switch regions, reveals that recombinations tend to occur at homologous positions of the respective gamma-unit repeats, adjacent to the elements AGCT and GGGG found in each. This finding suggests that the cutting and religation step of the recombination process is mediated by a recombinase common to the four gamma-isotypes.     

Recent and rapid amplification of the sperm basic nuclear protein genes in winter flounder

The high molecular weight basic nuclear proteins (HMrBNPs), which are tightly bound to sperm chromatin in winter flounder, are made up of imperfect reiterations of simple peptide sequences that contain phosphorylatable DNA-binding motifs. Genomic Southern blots hybridized with probes to the coding and non-coding regions of HMrBNP mRNA showed that HMrBNP sequences form a complex multi-gene family. Previously, one gene (2B) was used to establish an evolutionary link between histone H1 and the HMrBNPs. Further examination of this complex, multi-gene family has now revealed that the majority of the HMrBNP genes are linked as 4.5 kb direct tandem repeats that each contain a 2.8 kb coding region and a 1.7 kb intergenic region (IR). These findings, combined with the cloning of the IR, established that the tandemly repeated genes lack introns and code for the abundant 3 kb HMrBNP mRNAs that produce the prominent 110 kDa HMrBNP. Southern blotting of DNAs from other righteye flounder species showed that HMrBNP multi-gene families were present in closely related species, though with substantial differences in restriction patterns and band intensities, but were not detected in more distantly related flounders. These observations are consistent with recent and rapid elaboration of the HMrBNP gene family.     

Origin and evolution of genes and genomes. Crucial role of triplet expansions

A novel concept on mechanisms of evolution of genes and genomes is suggested: the sequences evolve largely by local events of triplet expansion and subsequent mutational changes in the repeats. The immediate memory about the earlier expansion events still resides in the sequences, in form of the frequently occurring segments of tandemly repeating codons. Other predicted fossils of the original repeats are: (I) the expanding triplets should be accompanied by their point mutation derivatives and (II) the remaining excess of codons formerly belonging to the tandem repeats should be reflected in overall codon usage biases. Both predictions are confirmed by analysis of largest available database of non-redundant protein coding sequences, of total size ～5?×?10(9) codons. One important conclusion also follows from the results. Life which, presumably, started with replication of expanding triplets and their subsequent mutational changes, is continuing to emerge within the genes and genomes, in form of new events of triplet expansion.