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.

Minisatellite instability in severe combined immunodeficiency mouse cells

We have recently found that okadaic acid, which shows strong inhibitory activity on protein serine/threonine phosphatases and tumor-promoting activity in vivo and in vitro, induces minisatellite mutation (MSM). Human tumors and chemically induced counterparts in experimental animals are also sometimes associated with MSM. In the present study, we demonstrated minisatellite (MS) instability in severe combined immunodeficiency (SCID) cells in which the DNA-dependent protein kinase catalytic subunit (DNA-PKcs) is impaired. Cells from a SCID fibroblast cell line transformed by simian virus 40 large tumor antigen, SC3VA2, and from an embryonal SCID fibroblast cell line, SC1K, were cloned and propagated to 107 to 108 cells, and then subjected to subcloning. After propagation of each subclone to 107 to 108 cells, DNA samples were digested with HinfI and analyzed by Southern blotting using the Pc-1 MS sequence as a probe. Under low-stringency conditions, about 40 MS bands were detected, with 45% ± 6% and 37% ± 3% of SC3VA2 and SC1K cells, respectively, having MSM. In contrast, cells from the RD13B2 cell line, which was established from SCVA2 by introducing human chromosome 8q fragments, on which DNA-PKcs is known to reside, to complement the SCID phenotype, showed a very low frequency of MSM (3% ± 3%). The high frequencies of MSM in SC3VA2 and SC1K were significant, with no difference between the two. The present study clearly demonstrates that MS instability exists in SCID fibroblasts, suggesting that DNA-PKcs might be involved in the stable maintenance of MS sequences in the genome.

Alternative genetic pathways in colorectal carcinogenesis

The comparative typing of matched tumor and blood DNAs at dinucleotide repeat (microsatellite) loci has revealed in tumor DNA the presence of alleles that are not observed in normal DNA. The occurrence of these additional alleles is possibly due to replication errors (RERs). Although this observation has led to the recognition of a subtype of colorectal cancer with a high incidence of RERs (caused by a deficiency in DNA mismatch repair), a thorough analysis of the RER frequency in a consecutive series of colorectal cancers had not been reported. It is shown here that the extensive typing of 88 colorectal tumors reveals a bimodal distribution for the frequency of RER at microsatellite loci. Within the major mode (75 tumors, RER− subtype), the probability that a locus exhibited instability did not differ significantly among loci and tumors, being 0.02. The subsequent development of a statistical test for an operational discrimination between the RER− and RER+ subtypes indicated that the probability of misclassification did not exceed 0.001 in this series. The frequency of K-ras mutation was found to be equivalent in the two subtypes. However, in the RER+ tumors, the p53 gene mutation was less frequently detected, the adenomatous polyposis coli (APC) mutation was rare, and the biallelic inactivation of either of these genes was not observed. Furthermore, the concomitant occurrence of APC and tumor growth factor β receptor type II gene alterations was found only once. These data suggest that the repertoires of genes that are frequently altered in RER+ and RER− tumors may be more different than previously thought.

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.

The onset and extent of genomic instability in sporadic colorectal tumor progression

Cancer cell genomes contain alterations beyond known etiologic events, but their total number has been unknown at even the order of magnitude level. By sampling colorectal premalignant polyp and carcinoma cell genomes through use of the technique inter-(simple sequence repeat) PCR, we have found genomic alterations to be considerably more abundant than expected, with the mean number of genomic events per carcinoma cell totaling approximately 11,000. Colonic polyps early in the tumor progression pathway showed similar numbers of events. These results indicate that, as with certain hereditary cancer syndromes, genomic destabilization is an early step in sporadic tumor development. Together these results support the model of genomic instability being a cause rather than an effect of malignancy, facilitating vastly accelerated somatic cell evolution, with the observed orderly steps of the colon cancer progression pathway reflecting the consequences of natural selection.

Genetic instability of cancer cells is proportional to their degree of aneuploidy

Genetic and phenotypic instability are hallmarks of cancer cells, but their cause is not clear. The leading hypothesis suggests that a poorly defined gene mutation generates genetic instability and that some of many subsequent mutations then cause cancer. Here we investigate the hypothesis that genetic instability of cancer cells is caused by aneuploidy, an abnormal balance of chromosomes. Because symmetrical segregation of chromosomes depends on exactly two copies of mitosis genes, aneuploidy involving chromosomes with mitosis genes will destabilize the karyotype. The hypothesis predicts that the degree of genetic instability should be proportional to the degree of aneuploidy. Thus it should be difficult, if not impossible, to maintain the particular karyotype of a highly aneuploid cancer cell on clonal propagation. This prediction was confirmed with clonal cultures of chemically transformed, aneuploid Chinese hamster embryo cells. It was found that the higher the ploidy factor of a clone, the more unstable was its karyotype. The ploidy factor is the quotient of the modal chromosome number divided by the normal number of the species. Transformed Chinese hamster embryo cells with a ploidy factor of 1.7 were estimated to change their karyotype at a rate of about 3% per generation, compared with 1.8% for cells with a ploidy factor of 0.95. Because the background noise of karyotyping is relatively high, the cells with low ploidy factor may be more stable than our method suggests. The karyotype instability of human colon cancer cell lines, recently analyzed by Lengnauer et al. [Lengnauer, C., Kinzler, K. W. & Vogelstein, B. (1997) Nature (London) 386, 623–627], also corresponds exactly to their degree of aneuploidy. We conclude that aneuploidy is sufficient to explain genetic instability and the resulting karyotypic and phenotypic heterogeneity of cancer cells, independent of gene mutation. Because aneuploidy has also been proposed to cause cancer, our hypothesis offers a common, unique mechanism of altering and simultaneously destabilizing normal cellular phenotypes.

Evolutionary instability of the major histocompatibility complex class I loci in New World primates

Homologues of the human major histocompatibility complex (MHC) HLA-A, -B, -E, -F, and -G loci are present in all the Catarrhini (Old World primates, apes, and humans), and some of their allelic lineages have survived several speciation events. Analysis of 26 MHC class I cDNAs from seven different genera of New World primates revealed that the Callitrichinae (tamarins and marmosets) are an exception to these rules of MHC stability. In gene trees of primate MHC class I genes, sequences from the Callitrichinae cluster in a genus-specific fashion, whereas in the other genera of New World primates, as in the Catarrhini, they cluster in a transgeneric way. The genus-specific clustering of the Callitrichinae cDNAs indicates that there is no orthology between MHC class I loci in genera of this phyletic group. Additionally, the Callitrichinae genera exhibit limited variability of their MHC class I genes, in contrast to the high variability displayed by all other primates. Each Callitrichinae genus, therefore, expresses its own set of MHC class I genes, suggesting that an unusually high rate of turnover of loci occurs in this subfamily. The limited variability of MHC class I genes in the Callitrichinae is likely the result of the recent origin of these loci.

The DNA Helicase Activity of BLM Is Necessary for the Correction of the Genomic Instability of Bloom Syndrome Cells

Bloom syndrome (BS) is a rare autosomal recessive disorder characterized by growth deficiency, immunodeficiency, genomic instability, and the early development of cancers of many types. BLM, the protein encoded by BLM, the gene mutated in BS, is localized in nuclear foci and absent from BS cells. BLM encodes a DNA helicase, and proteins from three missense alleles lack displacement activity. BLM transfected into BS cells reduces the frequency of sister chromatid exchanges and restores BLM in the nucleus. Missense alleles fail to reduce the sister chromatid exchanges in transfected BS cells or restore the normal nuclear pattern. BLM complements a phenotype of a Saccharomyces cerevisiae sgs1 top3 strain, and the missense alleles do not. This work demonstrates the importance of the enzymatic activity of BLM for its function and nuclear localization pattern.

Mutation rate varies among alleles at a microsatellite locus: Phylogenetic evidence

The understanding of the mutational mechanism that generates high levels of variation at microsatellite loci lags far behind the application of these genetic markers. A phylogenetic approach was developed to study the pattern and rate of mutations at a dinucleotide microsatellite locus tightly linked to HLA-DQB1 (DQCAR). A random Japanese population (n = 129) and a collection of multiethnic samples (n = 941) were typed at the DQB1 and DQCAR loci. The phylogeny of DQB1 alleles was then reconstructed and DQCAR alleles were superimposed onto the phylogeny. This approach allowed us to group DQCAR alleles that share a common ancestor. The results indicated that the DQCAR mutation rate varies drastically among alleles within this single microsatellite locus. Some DQCAR alleles never mutated during a long period of evolutionary time. Sequencing of representative DQCAR alleles showed that these alleles lost their ability to mutate because of nucleotide substitutions that shorten the length of uninterrupted CA repeat arrays; in contrast, all mutating alleles had relatively longer perfect CA repeat sequences.

Analysis of microsatellite mutations in the mitochondrial DNA of Saccharomyces cerevisiae

In the nuclear genome of Saccharomyces cerevisiae, simple, repetitive DNA sequences (microsatellites) mutate at rates much higher than nonrepetitive sequences. Most of these mutations are deletions or additions of repeat units. The yeast mitochondrial genome also contains many microsatellites. To examine the stability of these sequences, we constructed a reporter gene (arg8m) containing out-of-frame insertions of either poly(AT) or poly(GT) tracts within the coding sequence. Yeast strains with this reporter gene inserted within the mitochondrial genome were constructed. Using these strains, we showed that poly(GT) tracts were considerably less stable than poly(AT) tracts and that alterations usually involved deletions rather than additions of repeat units. In contrast, in the nuclear genome, poly(GT) and poly(AT) tracts had similar stabilities, and alterations usually involved additions rather than deletions. Poly(GT) tracts were more stable in the mitochondria of diploid cells than in haploids. In addition, an msh1 mutation destabilized poly(GT) tracts in the mitochondrial genome.

Chromosomal instability and cytoskeletal defects in oral cancer cells

Oral squamous cell carcinomas are characterized by complex, often near-triploid karyotypes with structural and numerical variations superimposed on the initial clonal chromosomal alterations. We used immunohistochemistry combined with classical cytogenetic analysis and spectral karyotyping to investigate the chromosomal segregation defects in cultured oral squamous cell carcinoma cells. During division, these cells frequently exhibit lagging chromosomes at both metaphase and anaphase, suggesting defects in the mitotic apparatus or kinetochore. Dicentric anaphase chromatin bridges and structurally altered chromosomes with consistent long arms and variable short arms, as well as the presence of gene amplification, suggested the occurrence of breakage–fusion–bridge cycles. Some anaphase bridges were observed to persist into telophase, resulting in chromosomal exclusion from the reforming nucleus and micronucleus formation. Multipolar spindles were found to various degrees in the oral squamous cell carcinoma lines. In the multipolar spindles, the poles demonstrated different levels of chromosomal capture and alignment, indicating functional differences between the poles. Some spindle poles showed premature splitting of centrosomal material, a precursor to full separation of the microtubule organizing centers. These results indicate that some of the chromosomal instability observed within these cancer cells might be the result of cytoskeletal defects and breakage–fusion–bridge cycles.

Mutational inactivation of the proapoptotic gene BAX confers selective advantage during tumor clonal evolution

A remarkable instability at simple repeated sequences characterizes gastrointestinal cancer of the microsatellite mutator phenotype (MMP). Mutations in the DNA mismatch repair gene family underlie the MMP, a landmark for hereditary nonpolyposis colorectal cancer. These tumors define a distinctive pathway for carcinogenesis because they display a particular spectrum of mutated cancer genes containing target repeats for mismatch repair deficiency. One such gene is BAX, a proapoptotic member of the Bcl-2 family of proteins, which plays a key role in programmed cell death. More than half of colon and gastric cancers of the MMP contain BAX frameshifts in a (G)8 mononucleotide tract. However, the functional significance of these mutations in tumor progression has not been established. Here we show that inactivation of the wild-type BAX allele by de novo frameshift mutations confers a strong advantage during tumor clonal evolution. Tumor subclones with only mutant alleles frequently appeared after inoculation into nude mice of single-cell clones of colon tumor cell lines with normal alleles. In contrast, no clones of BAX-expressing cells were found after inoculation of homozygous cell clones without wild-type BAX. These results support the interpretation that BAX inactivation contributes to tumor progression by providing a survival advantage. In this context, survival analyses show that BAX mutations are indicators of poor prognosis for both colon and gastric cancer of the MMP.

The human papillomavirus type 16 E6 and E7 oncoproteins cooperate to induce mitotic defects and genomic instability by uncoupling centrosome duplication from the cell division cycle

Loss of genomic integrity is a defining feature of many human malignancies, including human papillomavirus (HPV)-associated preinvasive and invasive genital squamous lesions. Here we show that aberrant mitotic spindle pole formation caused by abnormal centrosome numbers represents an important mechanism in accounting for numeric chromosomal alterations in HPV-associated carcinogenesis. Similar to what we found in histopathological specimens, HPV-16 E6 and E7 oncoproteins cooperate to induce abnormal centrosome numbers, aberrant mitotic spindle pole formation, and genomic instability. The low-risk HPV-6 E6 and E7 proteins did not induce such abnormalities. Whereas the HPV-16 E6 oncoprotein has no immediate effects on centrosome numbers, HPV-16 E7 rapidly induces abnormal centrosome duplication. Thus our results suggest a model whereby HPV-16 E7 induces centrosome-related mitotic disturbances that are potentiated by HPV-16 E6.

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.

Sequence-tagged microsatellite profiling (STMP): a rapid technique for developing SSR markers

We describe a technique, sequence-tagged microsatellite profiling (STMP), to rapidly generate large numbers of simple sequence repeat (SSR) markers from genomic or cDNA. This technique eliminates the need for library screening to identify SSR-containing clones and provides an ∼25-fold increase in sequencing throughput compared to traditional methods. STMP generates short but characteristic nucleotide sequence tags for fragments that are present within a pool of SSR amplicons. These tags are then ligated together to form concatemers for cloning and sequencing. The analysis of thousands of tags gives rise to a representational profile of the abundance and frequency of SSRs within the DNA pool, from which low copy sequences can be identified. As each tag contains sufficient nucleotide sequence for primer design, their conversion into PCR primers allows the amplification of corresponding full-length fragments from the pool of SSR amplicons. These fragments permit the full characterisation of a SSR locus and provide flanking sequence for the development of a microsatellite marker. Alternatively, sequence tag primers can be used to directly amplify corresponding SSR loci from genomic DNA, thereby reducing the cost of developing a microsatellite marker to the synthesis of just one sequence-specific primer. We demonstrate the utility of STMP by the development of SSR markers in bread wheat.