Previous Article | Next Article ![]()
Eukaryotic Cell, April 2008, p. 639-646, Vol. 7, No. 4
1535-9778/08/$08.00+0 doi:10.1128/EC.00447-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Ethan Dawson-Baglien, and
Barbara B. Sears*
Department of Plant Biology, Michigan State University, East Lansing, Michigan 48824-1312
Received 9 December 2007/ Accepted 30 January 2008
|
|
|---|
|
|
|---|
Because of their abundance and variability, microsatellite loci have been used extensively as genetic markers in evolutionary and ecological studies of natural populations in eukaryotes (14, 37) and also as highly polymorphic markers for forensics and genotyping of animals (39). Well-saturated microsatellite maps have been developed for nuclear genomes from a number of plants, including rice, maize, barley, Arabidopsis, and soybean (14, 23). In some cases, nuclear microsatellite instability is remarkably high: in chickpea, indels were found at an average rate of 2.9 x 10–3 to 10 x 10–3 per (TAA)n locus per generation (36). An even higher rate of somatic instability was found for mononucleotide microsatellites in reporter genes in Arabidopsis thaliana, with individual leaves having many sectors (2, 3).
Microsatellites and larger tandem repeats have also been identified in all chloroplast genomes from which sequence data are available (reviewed in reference 25), including the unicellular green alga Chlamydomonas reinhardtii (20). In barley, rice, and pine, chloroplast microsatellites have been used to reveal much higher levels of diversity than can be observed through traditional chloroplast restriction fragment length polymorphism analysis (reviewed in reference 25). From their observations of sequence data from Pinus species, Provan et al. (24) concluded that SSR length polymorphisms occur in chloroplast DNA (cpDNA) at frequencies of 3 x 10–5 to 8 x 10–5 per site per year.
The aforementioned studies provide a useful starting point for assessing microsatellite variability, but interpretation of SSR data based solely on an evolutionary context is problematic, because single changes can easily involve multiple repeats. Because the repeats are identical, independent deletions or duplications can result in the same DNA variations. Although such "synapomorphies" occur in parallel, they would be scored as arising from a single event. With a goal of providing an accurate assessment of the rate of chloroplast microsatellite variability, we created an experimental system that would allow us to monitor change in a microsatellite reporter in the cpDNA of the green alga C. reinhardtii. Since such mutational analyses require data from a large number of individuals, the algal system was chosen over a higher plant due to the ease in handling large populations in a limited space. In addition, the easy transformability and selection procedures for Chlamydomonas cpDNA have made it a uniquely accessible system for such a study. To create the microsatellite slippage reporter, a stretch of repeats was introduced into a cpDNA gene essential for photosynthesis (rbcL) so that it created a disruption in the reading frame. Restoration of photosynthetic competence could be achieved only by duplication or deletion events.
|
|
|---|
200 rpm under 245-mol per m2 per s photosynthetic photon flux (245-PPF) continuous "high" light or at about 10 PPF for "low-light" conditions. Growth curves were plotted from hemocytometer counts (13) of aliquots removed from 50-ml cultures grown in TAP medium lacking acetate (TMP). When testing for photosynthetic competence on plates, either TMP or HS (10) medium was used, with incubation under aquarium lights at
40 PPF. Crosses were performed as previously described (28). Zygospores were induced to germinate, and meiotic progeny were separated as described previously (13). Construction of the microsatellite reporter. The plasmids for chloroplast transformation, pIF24 and pOF28, were constructed by using a 3.9-kb BamHI-EcoRI subclone of p67 (9) named p698, which contains the wild-type rbcL gene, in the pUC8 vector. Plasmid p547, which contains the aadA gene with 5' and 3' control regions from psbA (kindly provided by Heriberto Cerutti), was used as the source of an aadA chloroplast expression cassette, which was inserted into the MfeI cut site of p698, giving plasmid p699 (Fig. 1). A 1.2-kb EcoRV-EcoRI fragment at the right end of p699 was subcloned in pBlueScript+ for inserting the slippage substrates at the BspMI cut site in the sixth codon position of rbcL to avoid other BspMI sites present in the larger clone. The engineered EcoRV-EcoRI fragment was then used to replace the equivalent portion of p699. The pIF24 construct contained six 4-bp repeats, giving a 24-bp in-frame insertion in the rbcL gene, while pOF28 contained seven repeats, resulting in a 28-bp insertion at the same site. These constructs were transformed into the dam mutant and dcm mutant GM2163 bacterial strain. The oligonucleotide pairs used to generate the 24-bp (IF24) and 28-bp (OF28) insertions were as follows:
|
View larger version (8K): [in a new window] |
FIG. 1. Diagram of the cpDNA segment in the transforming plasmids pIF24 and pOF28. The shaded boxes represent the chloroplast-borne genes of C. reinhardtii, and the arrows indicate the directions of their transcription. Two triangles represent the insertions in the transforming plasmid construct: the aadA cassette is the selectable marker that was inserted between two chloroplast genes; the IF24 (24-bp) and OF28 (28-bp) SSRs were inserted at the sixth codon position of the rbcL gene in plasmids pIF24 and pOF28, respectively. The enzymes used for cloning are indicated. The size of the reporter construct was 8.1 kb, including the pUC8 vector.
|
IF24-R, 3'-TTTC TTTC TTTC TTTC TTTC TTTG-5'
OF28-F, 5'-AAAC AAAG AAAG AAAG AAAG AAAG AAAG-3'
OF28-R, 3'-TTTC TTTC TTTC TTTC TTTC TTTC TTTG-5'
Small-scale plasmid isolations followed the alkaline lysis procedure of Sambrook et al. (26). Larger-scale preparations used a kit produced by Qiagen (Valencia, CA).
The biolistic method (5, 15) was used for chloroplast transformation, with CC-125 as the recipient strain and the plasmids pIF24 and pOF28 as the transforming DNAs. After the biolistic procedure, the cells were resuspended in TAP and immediately transferred to TAP agar supplemented with 100 µg/ml spectinomycin (Sigma Chemical Co., St. Louis, MO) for direct selection of the transformants.
Assay for changes in the microsatellite reporter in cpDNA of C. reinhardtii.
The cells with the out-of-frame microsatellite (OF28) were grown in liquid TAP culture to late log phase (2 x 106 to 4 x 106 cells/ml) on a rotary shaker for about 7 days at room temperature in low light. Cell counts were determined with a hemocytometer. The cells were concentrated about 10-fold by centrifugation. An aliquot of the concentrated cells was used for serial dilution onto agar-solidified TAP and was kept on low-light shelves to determine the number of viable cells plated. The remaining cells were plated as 500-µl aliquots on medium lacking acetate (either HS agar [10] or TMP agar) to test for photosynthetic competence. The plates were placed under lights (
40 PPF) on a shelf at room temperature.
DNA isolation, PCR amplification, and sequencing.
Chlamydomonas total genomic DNA was prepared as described previously (11). The rbcL5 (5'-GGCCCTTTCTATGCTCGACTG-3') and rbcLmid (5'-CCGAATACGTTACCTAC-3') primers were used in the PCR amplification reactions as the forward and reverse primers, respectively, to amplify a 560-bp segment (including the microsatellite insertion site) from the rbcL genes of the wild-type, OF, IF, and revertant colonies. In some experiments, an
180-bp PCR product was generated using a different set of primers, rbcLfor2 (5'-CTACGTAATCAGGTGTGTAG-3') and rbcLrev (5'-CCGGACAGATTAATTTTAGGA-3'), in order to readily visualize changes in the repeat number (Fig. 2). DNA (1 to 5 ng) was added to individual 100-µl reaction mixtures containing 200 µmol of each dideoxynucleotide, 2.5 units of Taq polymerase (Gibco BRL), 10 pmol of each primer, and 1x PCR buffer (Gibco BRL). Conditions for the 30 cycles of PCR were as follows: denaturation at 95°C for 2 min and 85°C for 5 min, followed by 94°C for 40 s, 30 s at 48°C to anneal the primers with the template, and 30 s at 72°C for DNA extension. A final extension was done for 10 min at 72°C. All of the PCRs were performed with a minicycler model PTC-150-16 (MR Research Inc., Watertown, MA). The PCR-amplified products were separated on 1.5 or 2% agarose gels in Tris-borate-EDTA buffer for the 560-bp and 180-bp products, respectively.
![]() View larger version (32K): [in a new window] |
FIG. 2. Visualization of DNA changes at the chloroplast microsatellite reporter. (A) Diagram of microsatellite indels with PCR primers. The original out-of-frame construct (OF28) contained seven 4-bp repeats (a). Deletion of one repeat (b) or four repeats (c) or insertion of two repeats (d) made the insertion in frame. (B) Diagram of differential mobility of PCR products.
|
N recA gene, the recA
N1 oligonucleotide (7) was used as the forward primer and an oligonucleotide called recA3end (5'-CTGGCATGCTTAAAAATCTTCGTTAGTTTC-3') was used as the reverse primer to generate a 1-kb PCR product. After the initial denaturation at 95°C for 2 min and 85°C for 5 min, the conditions for the PCR program were repeated for 30 cycles: 94°C for 1 min, 2 min at 36°C to anneal the primers and template, and 2 min at 68°C for DNA extension. The final extension was as described above. The PCR products were run on 1% agarose gels. PCR products from the rbcL5-rbcLmid amplification were submitted for sequencing by the Michigan State University Genomics Technology and Sequence Facility (http://www.genomics.msu.edu) using the internal primer rbcLfor2. Sequence analyses used Megalign Software from the DNAStar Program (Madison, WI).
Mutation rate. The microsatellite mutation rate was determined by averaging the frequencies of reversion to photosynthetic competence on parallel plates in two separate experiments and dividing the number of revertants by the number of viable cells calculated to have been spread on each plate.
Reconstruction experiment. Logarithmically growing cultures of two nonphotosynthetic strains, the OF28 reporter construct and the 18-7G rbcL point mutant, were diluted to 1 x 106 cells/ml, as were two photosynthetically competent strains, CC-124 and the IF24 control line. All four strains were taken through a 10-fold dilution series for individual platings on TAP and TMP media. From the dilutions of CC-124 and IF24 that should have had 100, 1,000, and 10,000 cells/ml, 20 µl was added to 2-ml aliquots of the nonphotosynthetic cell lines (1 x 106 cells/ml), and two 500-µl aliquots were plated on TMP medium and placed in the light.
|
|
|---|
![]() View larger version (20K): [in a new window] |
FIG. 3. Sequences and amino acid predictions of the rbcL replication slippage constructs. The wild-type rbcL gene is shown at the top, with boldface characters indicating the 4 bases of the 3' overhang at the BspMI cut site. The OF28 line has a 28-bp insertion in its rbcL gene, creating a stop codon that renders the cells nonphotosynthetic. The IF24 line has a six-codon insertion relative to the wild type. IF12 and IF36 are two photosynthetically competent lines recovered from the initial OF28 line.
|
![]() View larger version (23K): [in a new window] |
FIG. 4. Growth curves. Cultures were inoculated at 104 cells/ml, and cell densities were determined over a 6-day period. The out-of-frame (OF28) and in-frame (IF24) constructs were present in each of two independent transformants. CC-125 is a wild-type control, and 18-7G is an rbcL mutant.
|
|
View this table: [in a new window] |
TABLE 1. Frequency of mutation of microsatellite reporter
|
|
View this table: [in a new window] |
TABLE 2. Mixing experimenta
|
![]() View larger version (57K): [in a new window] |
FIG. 5. Slippage events at the cpDNA microsatellite reporter. Shown are examples of the PCR-amplified rbcL segment from the original out-of-frame insertion strain (OF28 [OF]), the wild-type C. reinhardtii strain (WT), and several photosynthetically competent colonies derived from the OF28 strain (1 to 3), using primers rbcLfor2 and rbcLrev. Lane M contains a 123-bp DNA ladder (Invitrogen).
|
![]() View larger version (35K): [in a new window] |
FIG. 6. Microsatellite changes in six photosynthetic colonies (1 to 6) isolated from the OF28 cell line. See the legend to Fig. 5 for abbreviations. The PCR products in lane 3 are 4 bp larger and 4 bp smaller than the original OF28 insert.
|
N recA construct (7) by classical crossing procedures. The E. coli
N recA gene has been shown to inhibit chloroplast recombination when inserted into cpDNA (7). Our first assessment tested the presence and activity of the
N recA gene in the original strain, CC-3455, by PCR amplification using appropriate primers to visualize the gene (described in Materials and Methods) and then by confirming the strain's hypersensitivity to 5 mM 5-fluorodeoxyuridine relative to wild-type cells, as described previously (7). This line was crossed to an mt– cell line containing the slippage construct with the linked aadA cassette. The spectinomycin resistance conferred by the aadA gene was used to select biparental progeny from which the desired recombinant could be isolated. Initially, very few biparental progeny were recovered, so the mt+ gametes were exposed to a handheld UV lamp for 2 min prior to mating, and we were able to obtain 5.8% biparental zygotes.
From dissected meiotic progeny, spectinomycin-resistant cells were selected and checked for the presence of both the
N recA gene and the slippage construct in the cpDNAs by PCR amplification using the recA
N1/recA3end and rbcLfor2/rbcLrev primer pairs, respectively. As represented by the meiotic progeny analyzed in Fig. 7B, about half of the biparental meiotic progeny (12 of 27) carried the
N recA gene. In contrast, all the progeny that showed resistance to spectinomycin also contained the 28-bp out-of-frame insertion in the closely linked rbcL gene (Fig. 7A). Two meiotic progeny that contained both the slippage construct and the
N recA gene were chosen for further experiments. Since the presence of
N recA should reduce chloroplast recombination, we tested the cell lines for the frequency of restoration of photosynthetic competence due to insertion/deletion events at the reporter gene. As shown in Table 3, the recovery of photosynthetically competent colonies from the slippage reporter in a recombination-deficient background occurred at a frequency of 5.3 x 10–5, a rate that is approximately three times higher than that of the recombination-proficient lines.
![]() View larger version (98K): [in a new window] |
FIG. 7. PCR analysis to identify meiotic progeny carrying both the cpDNA microsatellite reporter and the E. coli N42RecA gene. (A) PCR screening of spectinomycin-resistant progeny for the microsatellite reporter. PCR products from 13 progeny (1 to13) and the wild type (WT) were amplified with the rbcL5 and rbcLmid primers to allow visualization of the region containing the microsatellite. (B) PCR products indicating the presence of the N42RecA gene. The lanes are labeled as in panel A. The faint lower band is derived from excess oligonucleotides in the reaction mixture.
|
|
View this table: [in a new window] |
TABLE 3. Rates of spontaneous replication slippage events in the cpDNA of C. reinhardtii containing the OF28 construct in a recombination-deficient background
|
![]() View larger version (34K): [in a new window] |
FIG. 8. PCR products from E. coli subclones carrying the pOF28 plasmid. The subclones were randomly selected; DNA was extracted and amplified with oligonucleotides rbcLfor2 and rbcLrev. Lane M contains the 123-bp DNA ladder.
|
|
|
|---|
In order to achieve those goals, we used the green alga C. reinhardtii, which frequently serves as a model system for chloroplast studies. Since most cpDNA microsatellites are present in the noncoding regions of the plastid genome, a reporter construct was created to monitor replication slippage events in vivo using an essential gene of the cpDNA, rbcL, in order to provide a phenotypic selection for the slippage event. The rbcL gene product is essential for the fixation of carbon from CO2; an out-of-frame insertion in the gene would disrupt its reading frame, resulting in the death of the algal cells on media without an organic carbon source (acetate). Changes in the number of repeats within the microsatellite would be able to be observed easily in these cells: a deletion of one or four repeat units or an insertion of two repeats would result in photosynthetic competence. Since this assay required that the rbcL gene product be functional with the addition of extra codons, initial experiments were done to test whether an in-frame insertion of additional codons in the rbcL gene of C. reinhardtii would allow the cells to remain photosynthetically competent. As shown in Fig. 4, an in-frame insertion of 24 bp containing AAAG repeats (IF24) in the sixth codon position did not affect the photosynthetic competence or growth properties of the cells. However, an out-of-frame insertion of 28 bp (OF28) caused the cells to be sensitive to high light and to require acetate for growth, traits that are characteristic of rbcL– mutations (33). The OF28 construct thus served as a reporter to monitor microsatellite deletion and duplication events.
In experiments analogous to the Luria-Delbruck fluctuation tests (19), cells from large logarithmically growing cultures were subdivided among many selective plates. As shown in Table 1, an average of 150 revertants arose for every 10 million viable cells plated (a frequency of 1.5 x 10–5). A mixing experiment (Table 2) showed that if photosynthetic cells were present in the large culture, they would have produced green colonies within a week, whereas in our tracking of the microsatellite reporter, 4 weeks elapsed before the first green, photosynthetically competent colonies were observed. Because cell division stops soon after the cells are plated on the selective medium, the cpDNA microsatellite variability must occur within the nondividing cells. A change in the microsatellite that reestablishes the rbcL reading frame should be dominant; hence, the reversion would not have to reach fixation to produce a photosynthetically competent cell. However, if the microsatellite variation occurs through replication slippage, cpDNA replication must occur in spite of an arrest of cell division. This idea is not as unorthodox at it might at first seem, since the replication of Chlamydomonas cpDNA is not synchronized with its nuclear replication, but rather, occurs steadily during the growth of liquid cultures (reviewed in reference 27). Furthermore, evidence exists for turnover of the Chlamydomonas cpDNA (29).
Precautions were taken to exclude any preexisting event that might have been present in the starter cells of the assay: these included plating of the early preculture cells on media without acetate to test for the presence of early revertants and checking the sizes of the microsatellite insert with PCR amplification at the beginning of the culture. If anything, our calculation of the frequency is an underestimate, since only a third of the possible deletion or insertion events (those that reestablished the reading frame) would be detected by our assay. Multiplying our quantifications by 3 produces an estimated indel mutation rate for the 4-bp microsatellites of 4.5 x 10–5. This value is equivalent to estimates from evolutionary comparisons (24) for rates of cpDNA replication slippage of mononucleotide repeats 10 or more bases long in intergenic spacers. However, in our experimental system, mononucleotide repeats (a stretch of 13 As) varied at such a high frequency that as soon as we recovered a biolistically transformed cell line through the linked selectable marker, some cells in the colony had become photosynthetically competent due to indels within the microsatellite. Hence, we were unable to quantify the rate of mononucleotide repeat variation, but it must be several orders of magnitude higher than for the 4-bp microsatellite.
In our PCR screens of a subset of the photosynthetic revertants of the AAAG microsatellite, only one colony appeared to contain reciprocal-recombination products (Fig. 5). Hence, replication slippage seemed to be responsible for the production of the microsatellite indels within the chloroplast, as is true for other biological systems. However, to more thoroughly assess the contribution of recombination to microsatellite variability, the microsatellite reporter was combined with a chloroplast-located dominant-negative allele of the E. coli recA gene. To our surprise, the frequency of variation resulting in photosynthetically competent colonies in the recombination-deficient strain was about three times higher than in the recombination-proficient lines (Table 3). Since slippage that restores the rbcL reading frame should produce a dominant allele, we believe that this result indicates that the majority of these mutations are eliminated through recombination (gene conversion). Reducing recombination among the highly polyploid cpDNAs would enhance the survival and recovery of the revertant alleles.
With one exception, the slipped colonies that we recovered experienced deletion of one or four repeat units; the two deletion types occurred with equal frequencies. In one later experiment, we recovered a two-repeat-unit duplication that restored photosynthetic competence (IF36). Although photosynthetically competent, IF36 grows somewhat more slowly than the wild type, and hence, there may be some bias against the recovery of longer forms. Nonetheless, the fact that we were able to recognize and recover a duplication variant leads us to conclude that replication slippage in the chloroplast DNA of Chlamydomonas has a bias toward deletions rather than insertions. This is in marked contrast to what occurs to this same substrate when it is maintained in a plasmid in E. coli. As shown in Fig. 8, expansions rather than deletions typified the changes in the bacteria.
The observation that deletions are the prevalent form of replication slippage in the Chlamydomonas chloroplast suggests that the template strand of the cpDNA is more prone to slippage than is the daughter strand. A deletion bias during replication slippage has been reported for other organisms, including cpDNA of a number of plant species, such as petunia and alfalfa (1), mitochondrial microsatellites of yeast (31), and nuclear genomes of animals, such as snails (39). In the Chlamydomonas chloroplast genome, "short dispersed repeats" have been reported to have proliferated to a great abundance (20), and therefore, we were surprised to find a strong deletion bias. However, the short dispersed repeats described earlier (20) are significantly longer (30 to 36 bp) than those we analyzed and are rarely tandem direct repeats, so it is unlikely that they are substrates for replication slippage.
We have previously shown that base substitution at a particular target site in the cpDNA occurs at a rate of 0.15 x 10–9 to 11 x 10–9 cells in C. reinhardtii (11). This means that cpDNA replication slippage at a 4-base repeat occurs 1,000 to 100,000 times more frequently than base substitution. Although our results have the advantage of providing base substitution and replication slippage rates in the same organism, the values may not be directly comparable, since the base substitution rates were observed for the 16S rRNA gene within the cpDNA inverted repeat, while the slippage events occurred in a reporter construct in the single-copy region. Others have found that copy correction in the inverted repeats reduces the mutation rate by half (reviewed in reference 21).
In bacteria, the mismatch repair system defined by the MutHLS repair system is known to be the key surveillance system for correcting mismatches generated by replication slippage events (4). During replication through microsatellites, strand slippage results in the formation of a looped structure with unpaired bases on either the template or the daughter strand. Such a looped structure is recognized by the MutS protein of the mismatch repair pathway and targeted for excision (18). It has been hypothesized that in prokaryotes, a looped structure containing one to three unpaired bases in the loop is more readily recognized and repaired than is a repeat of four or more nucleotides (18, 22). Our chloroplast data paint a different picture; variations due to slippage of the AAAG repeats were definitely recovered less frequently than changes in a mononucleotide repeat, which were so frequent that we could not quantify them. The observation that the same 4-bp-repeat microsatellite in E. coli experiences daughter strand slippage at an extremely high rate suggests that endosymbiosis may have benefited the chloroplast by providing its genetic system with repair enzymes that minimize replication slippage of the daughter strand.
This work was supported by grant MCB 9982600 from the National Science Foundation to B.B.S. and was carried out in compliance with current laws governing genetic experimentation in the United States.
Published ahead of print on 8 February 2008. ![]()
Present address: Department of Natural Sciences, Notre Dame de Namur University, 1500 Ralston Avenue, Belmont, CA 94002. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»