Phylogenetic assessment of length variation at a microsatellite locus

Sixty-six haplotypes at a locus containing a simple dinucleotide (CA)n microsatellite repeat were isolated by PCR–single-strand conformational polymorphism from populations of the horseshoe crab Limulus polyphemus. These haplotypes were sequenced to assess nucleotide variation directly. Thirty-four distinct sequences (alleles) were identified in a region 570 bp long that included the microsatellite motif. In the repeat region itself, CA-number varied in integer values from 5 to 11 across alleles, except that a (CA)8 class was not observed. Differences among alleles were due also to polymorphisms at 22 sites in regions immediately flanking the microsatellite repeats. Nucleotide substitutions in these regions were used to estimate phylogenetic relationships among alleles, and the gene phylogeny was used to trace the evolution of length variation and CA repeat numbers. A low correlation between size variation and genealogical relationships among alleles suggests that absolute fragment size (as normally scored in microsatellite assays) is an unreliable indicator of historical affinities among alleles. This finding on the molecular fine structure of microsatellite variation suggests the need for caution in the use of repeat counts at microsatellite loci as secure indicators of allelic relationships.

Equilibrium distributions of microsatellite repeat length resulting from a balance between slippage events and point mutations

We describe and test a Markov chain model of microsatellite evolution that can explain the different distributions of microsatellite lengths across different organisms and repeat motifs. Two key features of this model are the dependence of mutation rates on microsatellite length and a mutation process that includes both strand slippage and point mutation events. We compute the stationary distribution of allele lengths under this model and use it to fit DNA data for di-, tri-, and tetranucleotide repeats in humans, mice, fruit flies, and yeast. The best fit results lead to slippage rate estimates that are highest in mice, followed by humans, then yeast, and then fruit flies. Within each organism, the estimates are highest in di-, then tri-, and then tetranucleotide repeats. Our estimates are consistent with experimentally determined mutation rates from other studies. The results suggest that the different length distributions among organisms and repeat motifs can be explained by a simple difference in slippage rates and that selective constraints on length need not be imposed.

Microsatellite spreading in the human genome: evolutionary mechanisms and structural implications.

Microsatellites are tandem repeat sequences abundant in the genomes of higher eukaryotes and hitherto considered as "junk DNA." Analysis of a human genome representative data base (2.84 Mb) reveals a distinct juxtaposition of A-rich microsatellites and retroposons and suggests their coevolution. The analysis implies that most microsatellites were generated by a 3'-extension of retrotranscripts, similar to mRNA polyadenylylation, and that they serve in turn as "retroposition navigators," directing the retroposons via homology-driven integration into defined sites. Thus, they became instrumental in the preservation and extension of primordial genomic patterns. A role is assigned to these reiterating A-rich loci in the higher-order organization of the chromatin. The disease-associated triplet repeats are mostly found in coding regions and do not show an association with retroposons, constituting a unique set within the family of microsatellite sequences.

Microsatellite variability and genetic distances.

We analyze the within- and between-population dynamics of the distribution of the number of repeats at multiple microsatellite DNA loci subject to stepwise mutation. Analytical expressions for moments up to the fourth order within a locus and the variance of between-locus variance at mutation-drift equilibrium have been obtained. These statistics may be used to test the appropriateness of the one-step mutation model and to detect between-locus variation in the mutation rate. Published data are compatible with the one-step mutation model, although they do not reject the two-step model. Using both multinomial sampling and diffusion approximations for the analysis of the genetic distance introduced by Goldstein et al. [Goldstein, D. B., Linares, A. R., Cavalli-Sforza, L. L. & Feldman, M. W. (1995) Proc. Natl. Acad. Sci. USA 92, 6723-6727], we show that this distance follows a chi 2 distribution with degrees of freedom equal to the number of loci when there is no variation in mutation rates among the loci. In the presence of such variation, the variance of the distance is obtained. We conclude that the number of microsatellite loci required for the construction of phylogenetic trees with reliable branch lengths may be several hundred. Also, mutations that change repeat scores by several units, even though extremely rare, may dramatically influence estimates of population parameters.

Orientation-dependent and sequence-specific expansions of CTG/CAG trinucleotide repeats in Saccharomyces cerevisiae

A quantitative and selective genetic assay was developed to monitor expansions of trinucleotide repeats (TNRs) in yeast. A promoter containing 25 repeats allows expression of a URA3 reporter gene and yields sensitivity to the drug 5-fluoroorotic acid. Expansion of the TNR to 30 or more repeats turns off URA3 and provides drug resistance. When integrated at either of two chromosomal loci, expansion rates were 1 × 10−5 to 4 × 10−5 per generation if CTG repeats were replicated on the lagging daughter strand. PCR analysis indicated that 5–28 additional repeats were present in 95% of the expanded alleles. No significant changes in CTG expansion rates occurred in strains deficient in the mismatch repair gene MSH2 or the recombination gene RAD52. The frequent nature of CTG expansions suggests that the threshold number for this repeat is below 25 in this system. In contrast, expansions of the complementary repeat CAG occurred at 500- to 1,000-fold lower rates, similar to a randomized (C,A,G) control sequence. When the reporter plasmid was inverted within the chromosome, switching the leading and lagging strands of replication, frequent expansions were observed only when CTG repeats resided on the lagging daughter strand. Among the rare CAG expansions, the largest gain in tract size was 38 repeats. The control repeats CTA and TAG showed no detectable rate of expansions. The orientation-dependence and sequence-specificity data support the model that expansions of CTG and CAG tracts result from aberrant DNA replication via hairpin-containing Okazaki fragments.

Induction of microsatellite instability by oxidative DNA damage

Instability of repetitive sequences, both in intronic sequences and within coding regions, has been demonstrated to be a hallmark of genomic instability in human cancer. Understanding how these mutational events arise may provide an opportunity for prevention or early intervention in cancer development. To study the source of this instability, we have identified a region of the β-lactamase gene that is tolerant to the insertion of fragments of exogenous DNA as large as 1,614 bp with minimal loss of enzyme activity, as determined by antibiotic resistance. Fragments inserted out-of-frame render Escherichia coli sensitive to antibiotic, and compensatory frameshift mutations that restore the reading frame of β-lactamase can be selected on the basis of antibiotic resistance. We have utilized this site to insert a synthetic microsatellite sequence within the β-lactamase gene and selected for mutations yielding frameshifts. This assay provides for detection of one frameshift mutation in a background of 106 wild-type sequences. Mismatch repair deficiency increased the observed frameshift frequency ≈300-fold. Exposure of plasmid containing microsatellite sequences to hydrogen peroxide resulted in frameshift mutations that were localized exclusively to the microsatellite sequences, whereas DNA damage by UV or N-methyl-N′-nitro-N-nitrosoguanidine did not result in enhanced mutagenesis. We postulate that in tumor cells, endogenous production of oxygen free radicals may be a major factor in promoting instability of microsatellite sequences. This β-lactamase assay may provide a sensitive methodology for the detection and quantitation of mutations associated with the development of cancer.

Epigenetic phenotypes distinguish microsatellite-stable and -unstable colorectal cancers

Aberrant DNA methylation is a common phenomenon in human cancer, but its patterns, causes, and consequences are poorly defined. Promoter methylation of the DNA mismatch repair gene MutL homologue (MLH1) has been implicated in the subset of colorectal cancers that shows microsatellite instability (MSI). The present analysis of four MspI/HpaII sites at the MLH1 promoter region in a series of 89 sporadic colorectal cancers revealed two main methylation patterns that closely correlated with the MSI status of the tumors. These sites were hypermethylated in tumor tissue relative to normal mucosa in most MSI(+) cases (31/51, 61%). By contrast, in the majority of MSI(−) cases (20/38, 53%) the same sites showed methylation in normal mucosa and hypomethylation in tumor tissue. Hypermethylation displayed a direct correlation with increasing age and proximal location in the bowel and was accompanied by immunohistochemically documented loss of MLH1 protein both in tumors and in normal tissue. Similar patterns of methylation were observed in the promoter region of the calcitonin gene that does not have a known functional role in tumorigenesis. We propose a model of carcinogenesis where different epigenetic phenotypes distinguish the colonic mucosa in individuals who develop MSI(+) and MSI(−) tumors. These phenotypes may underlie the different developmental pathways that are known to occur in these tumors.

Telomeres in the mouse have large inter-chromosomal variations in the number of T2AG3 repeats

The ultra-long telomeres that have been observed in mice are not in accordance with the concept that critical telomere shortening is related to aging and immortalization. Here, we have used quantitative fluorescence in situ hybridization to estimate (T2AG3)n lengths of individual telomeres in various mouse strains. Telomere lengths were very heterogeneous, but specific chromosomes of bone marrow cells and skin fibroblasts from individual mice had similar telomere lengths. We estimate that the shortest telomeres are around 10 kb in length, indicating that each mouse cell has a few telomeres with (T2AG3)n lengths within the range of human telomeres. These short telomeres may be critical in limiting the replicative potential of murine cells.

Complexities in the genetic structure of Anopheles gambiae populations in west Africa as revealed by microsatellite DNA analysis

Chromosomal forms of Anopheles gambiae, given the informal designations Bamako, Mopti, and Savannah, have been recognized by the presence or absence of four paracentric inversions on chromosome 2. Studies of karyotype frequencies at sites where the forms occur in sympatry have led to the suggestion that these forms represent species. We conducted a study of the genetic structure of populations of An. gambiae from two villages in Mali, west Africa. Populations at each site were composed of the Bamako and Mopti forms and the sibling species, Anopheles arabiensis. Karyotypes were determined for each individual mosquito and genotypes at 21 microsatellite loci determined. A number of the microsatellites have been physically mapped to polytene chromosomes, making it possible to select loci based on their position relative to the inversions used to define forms. We found that the chromosomal forms differ at all loci on chromosome 2, but there were few differences for loci on other chromosomes. Geographic variation was small. Gene flow appears to vary among different regions within the genome, being lowest on chromosome 2, probably due to hitchhiking with the inversions. We conclude that the majority of observed genetic divergence between chromosomal forms can be explained by forces that need not involve reproductive isolation, although reproductive isolation is not ruled out. We found low levels of gene flow between the sibling species Anopheles gambiae and Anopheles arabiensis, similar to estimates based on observed frequencies of hybrid karyotypes in natural populations.

Chromosomal arm replacement generates a high level of intraspecific polymorphism in the terminal inverted repeats of the linear chromosomal DNA of Streptomyces ambofaciens

The chromosomal DNA of the bacteria Streptomyces ambofaciens DSM40697 is an 8-Mb linear molecule that ends in terminal inverted repeats (TIRs) of 210 kb. The sequences of the TIRs are highly variable between the different linear replicons of Streptomyces (plasmids or chromosomes). Two spontaneous mutant strains harboring TIRs of 480 and 850 kb were isolated. The TIR polymorphism seen is a result of the deletion of one chromosomal end and its replacement by 480 or 850 kb of sequence identical to the end of the undeleted chromosomal arm. Analysis of the wild-type sequences involved in these rearrangements revealed that a recombination event took place between the two copies of a duplicated DNA sequence. Each copy was mapped to one chromosomal arm, outside of the TIR, and encoded a putative alternative sigma factor. The two ORFs, designated hasR and hasL, were found to be 99% similar at the nucleotide level. The sequence of the chimeric regions generated by the recombination showed that the chromosomal structure of the mutant strains resulted from homologous recombination events between the two copies. We suggest that this mechanism of chromosomal arm replacement contributes to the rapid evolutionary diversification of the sequences of the TIR in Streptomyces.

Characteristic enrichment of DNA repeats in different genomes

Using computer programs developed for this purpose, we searched for various repeated sequences including inverted, direct tandem, and homopurine–homopyrimidine mirror repeats in various prokaryotes, eukaryotes, and an archaebacterium. Comparison of observed frequencies with expectations revealed that in bacterial genomes and organelles the frequency of different repeats is either random or enriched for inverted and/or direct tandem repeats. By contrast, in all eukaryotic genomes studied, we observed an overrepresentation of all repeats, especially homopurine–homopyrimidine mirror repeats. Analysis of the genomic distribution of all abundant repeats showed that they are virtually excluded from coding sequences. Unexpectedly, the frequencies of abundant repeats normalized for their expectations were almost perfect exponential functions of their size, and for a given repeat this function was indistinguishable between different genomes.

Control of alternative pre-mRNA splicing by distributed pentameric repeats

Multiple copies of the hexamer TGCATG have been shown to regulate fibronectin pre-mRNA alternative splicing. GCATG repeats also are clustered near the regulated calcitonin-specific 3′ splice site in the rat calcitonin/CGRP gene. Specific mutagenesis of these repeats in calcitonin/CGRP pre-mRNA resulted in the loss of calcitonin-specific splicing, suggesting that the native repeats act to enhance alternative exon inclusion. Mutation of subsets of these elements implies that alternative splicing requires a minimum of two repeats, and that the combination of one intronic and one exonic repeat is necessary for optimal cell-specific splicing. However, multimerized intronic repeats inhibited calcitonin-specific splicing in both the wild-type context and in a transcript lacking endogenous repeats. These results suggest that both the number and distribution of repeats may be important features for the regulation of tissue-specific alternative splicing. Further, RNA containing a single repeat bound cell-specific protein complexes, but tissue-specific differences in protein binding were not detected by using multimerized repeats. Together, these data support a novel model for alternative splicing regulation that requires the cell-specific recognition of multiple, distributed sequence elements.

Multiple imprinted sense and antisense transcripts, differential methylation and tandem repeats in a putative imprinting control region upstream of mouse Igf2

The mouse insulin-like growth factor 2 (Igf2) locus is a complex genomic region that produces multiple transcripts from alternative promoters. Expression at this locus is regulated by parental imprinting. However, despite the existence of putative imprinting control elements in the Igf2 upstream region, imprinted transcriptional repression is abolished by null mutations at the linked H19 locus. To clarify the extent to which the Igf2 upstream region contains autonomous imprinting control elements we have performed functional and comparative analyses of the region in the mouse and human. Here we report the existence of multiple, overlapping imprinted (maternally repressed) sense and antisense transcripts that are associated with a tandem repeat in the mouse Igf2 upstream region. Regions flanking the repeat exhibit tissue-specific parental allelic methylation patterns, suggesting the existence of tissue-specific control elements in the upstream region. Studies in H19 null mice indicate that both parental allelic methylation and monoallelic expression of the upstream transcripts depends on an intact H19 gene acting in cis. The homologous region in human IGF2 is structurally conserved, with the significant exception that it does not contain a tandem repeat. Our results support the proposal that tandem repeats act to target methylation to imprinted genetic loci.

Identification of Darlin, a Dictyostelium Protein with Armadillo-like Repeats that Binds to Small GTPases and Is Important for the Proper Aggregation of Developing Cells

We purified from Dictyostelium lysates an 88-kDa protein that bound to a subset of small GTPases, including racE, racC, cdc42Hs, and TC4ran, but did not bind to R-ras or rabB. Cloning of the gene encoding this 88-kDa protein revealed that it contained multiple armadillo-like repeats most closely related to the mammalian GTP exchange factor smgGDS. We named this protein darlin (Dictyostelium armadillo-like protein). Disruption of the gene encoding darlin demonstrated that this protein is not essential for cytokinesis, pinocytosis, phagocytosis, or development. However, the ability of darlin null cells to aggregate in response to starvation is severely affected. When starved under liquid medium, the mutant cells were unable to form aggregation centers and streams, possibly because of a defect in cAMP relay signaling. This defect was not due to an inability of the darlin mutants to activate adenylate cyclase in response to G protein stimulation. These results suggest that the darlin protein is involved in a signaling pathway that may modulate the chemotactic response during early development.

Xenopus Kielin: A dorsalizing factor containing multiple chordin-type repeats secreted from the embryonic midline

The midline tissues are important inductive centers of early vertebrate embryos. By signal peptide selection screening, we isolated a secreted factor, Kielin, which contains multiple cys-rich repeats similar to those in chordin (Chd). Expression of Kielin starts at midgastrula stages in the notochord and is detected in the floor plate of neurula embryos. Kielin is induced in mesoderm and in ectoderm by nodal-related genes. Chd is sufficient to activate Kielin expression in mesoderm whereas Shh or HNF-3β in addition to Chd is required for induction in ectoderm. Kielin has a distinct biological activity from that of Chd. Injection of Kielin mRNA causes dorsalization of ventral marginal zone explants and expansion of MyoD expression in neurula embryos. Unlike Chd, Kielin does not efficiently induce neural differentiation of animal cap ectoderm, suggesting that the activity of Kielin is not simply caused by BMP4 blockade. Kielin is a signaling molecule that mediates inductive activities of the embryonic midline.