Concerted evolution of the tandem array encoding primate U2 snRNA (the RNU2 locus) is accompanied by dramatic remodeling of the junctions with flanking chromosomal sequences.

The genes encoding primate U2 snRNA are organized as a nearly perfect tandem array (the RNU2 locus) that has been evolving concertedly for >35 Myr since the divergence of baboons and humans. Thus the repeat units of the tandem array are essentially identical within each species, but differ between species. Homogeneity is maintained because any change in one repeat unit is purged from the array or fixed in all other repeats. Intriguingly, the cytological location of RNU2 has remained unchanged despite concerted evolution of the tandem array. We had found previously that junction sequences between the U2 tandem array and flanking DNA were subject to remodeling over a region of 200-300 bp during the past 5 Myr in the hominid lineage. Here we show that the junctions between the U2 tandem array and flanking DNA have undergone dramatic rearrangements over a region of 1 to >10 kbp in the 35 Myr since divergence of the Old World Monkey and hominid lineages. We argue that these rearrangements reflect the high level of genetic activity required to sustain concerted evolution, and propose a model to explain why maintenance of homogeneity within a tandemly repeated multigene family would lead to junctional diversity.     

Concerted Evolution of the Tandemly Repeated Genes Encoding Primate U2 Small Nuclear RNA (the RNU2 Locus) Does Not Prevent Rapid Diversification of the (CT)n · (GA)n Microsatellite Embedded within the U2 Repeat Unit

The RNU2 locus encoding human U2 small nuclear RNA (snRNA) is organized as a nearly perfect tandem array containing 5 to 22 copies of a 5.8-kb repeat unit. Just downstream of the U2 snRNA gene in each 5.8-kb repeat unit lies a large (CT)_n · (GA)_n dinucleotide repeat (n ≈ 70). This form of genomic organization, in which one repeat is embedded within another, provides an unusual opportunity to study the balance of forces maintaining the homogeneity of both kinds of repeats. Using a combination of field inversion gel electrophoresis and polymerase chain reaction, we have been able to study the CT microsatellites within individual U2 tandem arrays. We find that the CT microsatellites within an RNU2 allele exhibit significant length polymorphism, despite the remarkable homogeneity of the surrounding U2 repeat units. Length polymorphism is due primarily to loss or gain of CT dinucleotide repeats, but other types of deletions, insertions, and substitutions are also frequent. Polymorphism is greatly reduced in regions where pure (CT)_n tracts are interrupted by occasional G residues, suggesting that irregularities stabilize both the length and the sequence of the dinucleotide repeat. We further show that the RNU2 loci of other catarrhine primates (gorilla, chimpanzee, orangutan, and baboon) contain orthologous CT microsatellites; these also exhibit length polymorphism, but are highly divergent from each other. Thus, although the CT microsatellite is evolving far more rapidly than the rest of the U2 repeat unit, it has persisted through multiple speciation events spanning >35 Myr. The persistence of the CT microsatellite, despite polymorphism and rapid evolution, suggests that it might play a functional role in concerted evolution of the RNU2 loci, perhaps as an initiation site for recombination and/or gene conversion.     

Concerted evolution of the tandem array encoding primate U2 snRNA occurs in situ, without changing the cytological context of the RNU2 locus.

In primates, the tandemly repeated genes encoding U2 small nuclear RNA evolve concertedly, i.e. the sequence of the U2 repeat unit is essentially homogeneous within each species but differs somewhat between species. Using chromosome painting and the NGFR gene as an outside marker, we show that the U2 tandem array (RNU2) has remained at the same chromosomal locus (equivalent to human 17q21) through multiple speciation events over > 35 million years leading to the Old World monkey and hominoid lineages. The data suggest that the U2 tandem repeat, once established in the primate lineage, contained sequence elements favoring perpetuation and concerted evolution of the array in situ, despite a pericentric inversion in chimpanzee, a reciprocal translocation in gorilla and a paracentric inversion in orang utan. Comparison of the 11 kb U2 repeat unit found in baboon and other Old World monkeys with the 6 kb U2 repeat unit in humans and other hominids revealed that an ancestral U2 repeat unit was expanded by insertion of a 5 kb retrovirus bearing 1 kb long terminal repeats (LTRs). Subsequent excision of the provirus by homologous recombination between the LTRs generated a 6 kb U2 repeat unit containing a solo LTR. Remarkably, both junctions between the human U2 tandem array and flanking chromosomal DNA at 17q21 fall within the solo LTR sequence, suggesting a role for the LTR in the origin or maintenance of the primate U2 array.     

Structure and evolution of the U2 small nuclear RNA multigene family in primates: gene amplification under natural selection?

The organization of U2 genes was compared in apes, Old World monkeys, and the prosimian galago. In humans and all apes (gibbon, orangutan, gorilla, and chimpanzee), the U2 genes were organized as a tandem repeat of a 6-kb element; however, the restriction maps of the 6-kb elements in these divergent species differed slightly, demonstrating that mechanisms must exist for maintaining sequence homogeneity within this tandem array. In Old World monkeys, the U2 genes were organized as a tandem repeat of an 11-kb element; the restriction maps of the 11-kb elements in baboon and two closely related macaques, bonnet and rhesus monkeys, also differed slightly, confirming that efficient sequence homogenization is an intrinsic property of the U2 tandem array. Interestingly, the 11-kb monkey repeat unit differed from the 6-kb hominid repeat unit by a 5-kb block of monkey-specific sequence. Finally, we found that the U2 genes of the prosimian galago were dispersed rather than tandemly repeated, suggesting that the hominid and Old World monkey U2 tandem arrays resulted from independent amplifications of a common ancestral U2 gene. Alternatively, the 5-kb monkey-specific sequence could have been inserted into the 6-kb array or deleted from the 11-kb array soon after divergence of the hominid and Old World monkey lineages.     

Repeat instability at human minisatellites arising from meiotic recombination.

Little is known about the role of meiotic recombination processes such as unequal crossover in driving instability at tandem repeat DNA. Methods have therefore been developed to detect meiotic crossovers within two different GC-rich minisatellite repeat arrays in humans, both in families and in sperm DNA. Both loci normally mutate in the germline by complex conversion-like transfer of repeats between alleles. Analysis shows that inter-allelic unequal crossovers also occur at both loci, although at low frequency, to yield simple recombinant repeat arrays with exchange of flanking markers. Equal crossovers between aligned alleles, resulting in recombinant alleles but without change in repeat copy number, also occur in sperm at a similar frequency to unequal crossovers. Both crossover and conversion show polarity in the repeat array and are co-suppressed in an allele showing unusual germline stability. This provides evidence that minisatellite conversion and crossover arise by a common mechanism, perhaps by alternative processing of a meiotic recombination initiation complex, and implies that minisatellite instability is a by-product of meiotic recombination in repeat DNA. While minisatellite recombination is infrequent, crossover rates indicate that the unstable end of a human minisatellite can act as a recombination warm-spot, even between sequence-heterologous alleles.     

Meiotic recombination and flanking marker exchange at the highly unstable human minisatellite CEB1 (D2S90)

Unequal crossover has long been suspected to play a role in the germline-specific instability of tandem-repeat DNA, but little information exists on the dynamics and processes of unequal exchange. We have therefore characterized new length alleles associated with flanking-marker exchange at the highly unstable human minisatellite CEB1, which mutates in the male germline by a complex process often resulting in the gene conversion-like transfer of repeats between alleles. DNA flanking CEB1 is rich in single-nucleotide polymorphisms (SNPs) and shows extensive haplotype diversity, consistent with elevated recombinational activity near the minisatellite. These SNPs were used to recover mutant CEB1 molecules associated with flanking-marker exchange, directly from sperm DNA. Mutants with both proximal and distal flanking-marker exchange were shown to contribute significantly to CEB1 turnover and suggest that the 5' end of the array is very active in meiotic unequal crossover. Coconversions involving the interallelic transfer of repeats plus immediate flanking DNA were also common, were also polarized at the 5' end of CEB1, and appeared to define a conversion gradient extending from the repeat array into adjacent DNA. Whereas many mutants associated with complete exchange resulted in simple recombinant-repeat arrays that show reciprocity, coconversions were highly gain-biased and were, on average, more complex, with allele rearrangements similar to those seen in the bulk of sperm mutants. This suggests distinct recombination-processing pathways producing, on the one hand, simple crossovers in CEB1 and, on the other hand, complex conversions that sometimes extend into flanking DNA.     

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.     

Genotyping HapSTR loci: phase determination from direct sequencing of PCR products

HapSTRs combine information from a microsatellite (or simple tandem repeat, STR) with one or more single nucleotide polymorphisms in the DNA sequence immediately flanking the STR. These loci may offer increased power for the estimation of demographic parameters, but also present some challenges for data collection and analysis. We describe a process for inferring HapSTR alleles, including the flanking haplotypes, STR alleles and their phase relative to each other, directly from DNA sequence electropherograms of PCR products from heterozygous individuals. Our approach eliminates the need for more costly and time-consuming processes, such as cloning or acrylamide gel electrophoresis to separate alleles prior to sequencing.     

Tandem repeat copy-number variation in protein-coding regions of human genes

Background

Tandem repeat variation in protein-coding regions will alter protein length and may introduce frameshifts. Tandem repeat variants are associated with variation in pathogenicity in bacteria and with human disease. We characterized tandem repeat polymorphism in human proteins, using the UniGene database, and tested whether these were associated with host defense roles.

Results

Protein-coding tandem repeat copy-number polymorphisms were detected in 249 tandem repeats found in 218 UniGene clusters; observed length differences ranged from 2 to 144 nucleotides, with unit copy lengths ranging from 2 to 57. This corresponded to 1.59% (218/13,749) of proteins investigated carrying detectable polymorphisms in the copy-number of protein-coding tandem repeats. We found no evidence that tandem repeat copy-number polymorphism was significantly elevated in defense-response proteins (p = 0.882). An association with the Gene Ontology term 'protein-binding' remained significant after covariate adjustment and correction for multiple testing. Combining this analysis with previous experimental evaluations of tandem repeat polymorphism, we estimate the approximate mean frequency of tandem repeat polymorphisms in human proteins to be 6%. Because 13.9% of the polymorphisms were not a multiple of three nucleotides, up to 1% of proteins may contain frameshifting tandem repeat polymorphisms.

Conclusion

Around 1 in 20 human proteins are likely to contain tandem repeat copy-number polymorphisms within coding regions. Such polymorphisms are not more frequent among defense-response proteins; their prevalence among protein-binding proteins may reflect lower selective constraints on their structural modification. The impact of frameshifting and longer copy-number variants on protein function and disease merits further investigation.     

Intragenic duplication and divergence in the spectrin superfamily of proteins

Many structural, signaling, and adhesion molecules contain tandemly repeated amino acid motifs. The alpha-actinin/spectrin/dystrophin superfamily of F-actin-crosslinking proteins contains an array of triple alpha-helical motifs (spectrin repeats). We present here the complete sequence of the novel beta-spectrin isoform beta(Heavy)- spectrin (beta H). The sequence of beta H supports the origin of alpha- and beta-spectrins from a common ancestor, and we present a novel model for the origin of the spectrins from a homodimeric actin-crosslinking precursor. The pattern of similarity between the spectrin repeat units indicates that they have evolved by a series of nested, nonuniform duplications. Furthermore, the spectrins and dystrophins clearly have common ancestry, yet the repeat unit is of a different length in each family. Together, these observations suggest a dynamic period of increase in repeat number accompanied by homogenization within each array by concerted evolution. However, today, there is greater similarity of homologous repeats between species than there is across repeats within species, suggesting that concerted evolution ceased some time before the arthropod/vertebrate split. We propose a two-phase model for the evolution of the spectrin repeat arrays in which an initial phase of concerted evolution is subsequently retarded as each new protein becomes constrained to a specific length and the repeats diverge at the DNA level. This evolutionary model has general applicability to the origins of the many other proteins that have tandemly repeated motifs.      

Molecular analysis of a polymorphic domain of alpha satellite from the human X chromosome.

Alpha satellite DNA, a diverse family of tandemly repeated DNA sequences located at the centromeric region of each human chromosome, is organized in a highly chromosome-specific manner and is characterized by a high frequency of restriction-fragment-length polymorphism. To examine events underlying the formation and spread of these polymorphisms within a tandem array, we have cloned and sequenced a representative copy of a polymorphic array from the X chromosome and compared this polymorphic copy with the predominant higher-order repeat form of X-linked alpha satellite. Sequence data indicate that the polymorphism arose by a single base mutation that created a new restriction site (for HindIII) in the sequence of the predominant repeat unit. This variant repeat unit, marked by the new HindIII site, was subsequently amplified in copy number to create a polymorphic domain consisting of approximately 500 copies of the variant repeat unit within the X-linked array of alpha satellite. We propose that a series of intrachromosomal recombination events between misaligned tandem arrays, involving multiple rounds of either unequal crossing-over or sequence conversion, facilitated the spread and fixation of this variant HindIII repeat unit.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Unusual evolution of interspersed repeat sequences in the Drosophila ananassae subgroup

New repeat sequences were found in the Drosophila ananassae genome sequence. They accounted for approximately 1.2% of the D. ananassae genome and were estimated to be more abundant in genomes of its closely related species belonging to the Drosophila bipectinata complex, whereas it was entirely absent in the Drosophila melanogaster genome. They were interspersed throughout euchromatic regions of the genome, usually as short tandem arrays of unit sequences, which were mostly 175-200 bp long with two distinct peaks at 180 and 189 bp in the length distribution. The nucleotide differences among unit sequences within the same array (locus) were much smaller than those between separate loci, suggesting within-locus concerted evolution. The phylogenetic tree of the repeat sequences from different loci showed that divergences between sequences from different chromosome arms occurred only at earlier stages of evolution, while those within the same chromosome arm occurred thereafter, resulting in the increase in copy number. We found RNA polymerase III promoter sequences (A box and B box), which play a critical role in retroposition of short interspersed elements. We also found conserved stem-loop structures, which are possibly associated with certain DNA rearrangements responsible for the increase in copy number within a chromosome arm. Such an atypical combination of characteristics (i.e., wide dispersal and tandem repetition) may have been generated by these different transposition mechanisms during the course of evolution.     

Preferential homogenization between adjacent and alternate subrepeats in wheat rDNA.

DNA from the "non-transcribed spacer" (NTS) of two wheat ribosomal RNA gene (rDNA) clones was sequenced. The regions flanking the internal subrepeat arrays are highly conserved between the two clones; the nucleotide sequence differ by less than one-half percent. In contrast, the consensus sequences of the subrepeats in the two arrays differ by three percent. Mutations unique to each array, yet found in more than one subrepeat of the array, are preferentially found in adjacent and alternate subrepeats. The similarity of the DNA sequences of the flanking regions is consistent with a model of homogenization among rDNA gene units by intergenic conversion. We propose that a different mechanism, preferential conversion between neighboring subrepeats, is largely responsible for the homogenization of subrepeats within an array.     

Evolution of Tandemly Repeated Sequences: What Happens at the End of an Array?

Tandemly repeated sequences are a major component of the eukaryotic genome. Although the general characteristics of tandem repeats have been well documented, the processes involved in their origin and maintenance remain unknown. In this study, a region on the paternal sex ratio (PSR) chromosome was analyzed to investigate the mechanisms of tandem repeat evolution. The region contains a junction between a tandem array of PSR2 repeats and a copy of the retrotransposon NATE, with other dispersed repeats (putative mobile elements) on the other side of the element. Little similarity was detected between the sequence of PSR2 and the region of NATE flanking the array, indicating that the PSR2 repeat did not originate from the underlying NATE sequence. However, a short region of sequence similarity (11/15 bp) and an inverted region of sequence identity (8 bp) are present on either side of the junction. These short sequences may have facilitated nonhomologous recombination between NATE and PSR2, resulting in the formation of the junction. Adjacent to the junction, the three most terminal repeats in the PSR2 array exhibited a higher sequence divergence relative to internal repeats, which is consistent with a theoretical prediction of the unequal exchange model for tandem repeat evolution. Other NATE insertion sites were characterized which show proximity to both tandem repeats and complex DNAs containing additional dispersed repeats. An ``accretion model' is proposed to account for this association by the accumulation of mobile elements at the ends of tandem arrays and into ``islands' within arrays. Mobile elements inserting into arrays will tend to migrate into islands and to array ends, due to the turnover in the number of intervening repeats. Received: 18 August 1997 / Accepted: 18 September 1998     

Origins of polymorphism at a polypurine hypervariable locus.

We present characterisation of a hypervariable locus, D8S210, mapped to the telomeric region of the short arm of chromosome 8. The locus is highly polymorphic with alleles varying in size from 1.8 kb to 24 kb. Sequence data from 7 alleles shows that the variable region is entirely polypurine on one strand with a tetranucleotide repeating unit GGAA at the margins and diverged versions of this motif internally. The margins are conserved between alleles; polymorphism occurring in the internal regions of the repeat. Alleles are inherited in a Mendelian manner and one new mutation has been observed in analysis of 51 meioses. Use of single copy flanking sequences to elaborate the polymorphism revealed loss of single copy DNA in 3 unrelated families and in 2 other unrelated individuals. Restriction mapping shows that this loss is similar for different sized alleles in all three families suggesting that it was an early event that may have involved a flanking Alu sequence. We present evidence that the polypurine region can adopt triplex conformations in vitro. Such structures may facilitate loss or gain of unique sequences in the genome, contribute to mutation at conformation transition points and drive the hypervariability (> 99% heterozygosity) of this locus.     

The interleukin-1 receptor antagonist gene and the inhibitor of kappa B-like protein gene polymorphisms are not associated with myocardial infarction in Slovene population with type 2 diabetes

In this study we investigated the association of the interleukin-1 receptor antagonist gene variable number tandem repeat (IL1RN VNTR) polymorphism and of the inhibitor of kappa B-like protein (IKBL) gene polymorphism with myocardial infarction (MI) in a group of patients with type 2 diabetes. The IL1RN VNTR and the IKBL+ 738T > C gene polymorphisms were tested in 374 Caucasians: 151 cases with MI and 223 subjects with no history of coronary artery disease. The IL1RN VNTR polymorphism was not a risk factor for MI in Caucasians with type 2 diabetes (genotype 22 vs. the rest: odds ratio (OR) 1.6; 95% confidence interval (CI) = 0.8-3.5; p = 0.2). We also failed to demonstrate that IKBL+ 738T > C gene polymorphism was associated with MI in patients with type 2 diabetes (OR = 0.9; 95% CI = 0.3-2.6; p = 0.9). We provide evidence that the IL1RN VNTR and the IKBL + 738T > C gene polymorphisms are not risk factors for MI in Caucasians with type 2 diabetes.