School of Biology, University of Leeds, Leeds, United Kingdom
Correspondence: E-mail: ejp20{at}le.ac.uk.
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Key Words: concerted evolution Chorthippus parallelus gene conversion ITS1 ribosomal DNA homogenization
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Ribosomal DNA units do not evolve independently but, instead, appear to evolve in a concerted fashion such that copies are similar within a species and different between species. Their evolution is dependent on the interactions among natural selection, genetic drift, and molecular turnover mechanisms (Dover 1982). The relative contributions of these phenomena are poorly understood and are likely to vary among species. Because of a lack of empirical data, the process of concerted evolution has primarily been investigated using mathematical models and computer simulations (Nagylaki and Petes 1982; Ohta and Dover 1983, 1984; Walsh 1987; Stephan 1989). The results of these studies suggest that molecular turnover mechanisms such as gene conversion and unequal crossing-over could play a significant role in the concerted evolution of repeated sequences but their impact depends heavily on the rates of these processes relative to recombination and mutation (reviewed by Elder and Turner 1995).
Although it has been shown that the ITS1 region is not freely evolving (Schlötterer et al. 1994), the level of selective constraint acting on this region is expected to be minimal. Consequently, the pattern of ITS1 sequence variation at the intra-individual level will primarily depend on the spread of base substitutions or sequence rearrangements caused by molecular turnover mechanisms. Among individuals, the level of variation observed will depend on the rate at which random chromosomal assortment (meiosis) disperses these rearrangements among the individuals of the next generation and on genetic drift. It has been proposed that the interaction between molecular turnover and sexual reproduction will maintain all individuals of a population in a similar state at any one point in time (Dover 1982). This theory has important implications because an absence of significant variation among individuals will remove the potential for selection to act, meaning that molecular turnover mechanisms could potentially become the driving force behind the evolution of regions such as ITS1. By quantifying sequence variation at the individual, population, and subspecies level, this study investigated whether the populations of ITS1 copies within individuals are in a similar evolutionary state among individuals of the meadow grasshopper, Chorthippus parallelus (Orthoptera: Acrididae). In addition, these data provide an insight into the way in which the ITS1 region is evolving and indicate which molecular turnover mechanisms are significantly contributing to this process.
C. parallelus represents an ideal model system, as the recent evolutionary history of the species has been intensively studied and, in contrast to Drosophila, the genome of this species contains multiple NORs. During the last ice age(s) this species was repeatedly confined to southerly refugia (Hewitt 1993). When C. parallelus expanded its range northwards again at the end of the last ice age (10,000 years before present (YBP)), hybrid zones formed where the different refugial populations met. The level of divergence that had accrued was sufficient for the French population (originating from the Balkan refuge) and the Spanish population (derived from a refuge in southern Spain) to be regarded as distinct subspecies, namely Chorthippus parallelus parallelus (Cpp) and Chorthippus parallelus erythropus (Cpe), respectively (Reynolds 1980). Based on mitochondrial sequence data (COI), the estimated time of divergence for these two species is 500,000 years (Lunt, Ibrahim, and Hewitt 1998), which is equivalent to 500,000 generations.
A diagnostic difference between Cpp and Cpe involves their haploid complement of NORs. Both subspecies have two autosomal NORs (on chromosomes L2 and L3), but Cpp has an additional NOR located on its X chromosome (Gosalvez et al. 1990). It is currently unclear whether the common ancestor had two or three NORs, but the discovery that Italian C. parallelus have two NORs (Flanagan et al. 1999) suggests that a NOR gain may have occurred recently in the ancestor of Cpp.
Sequencing multiple ITS1 copies from several C. parallelus individuals (Parkin 2002) revealed the presence of polymorphic nucleotide sites that will be referred to in this paper as "within-individual polymorphisms" (WIPs). Four of these WIPs were targeted in this study (identified as WIP1, 2, 3, and 4), chosen because they are present in more than one individual with the coexisting base types occurring at intermediate frequencies. These polymorphisms represent the transition period when a new mutation has arisen in one rDNA copy and begun to spread to other copies in the genome. Depending on the rate of spread, these substitutions may now be species, subspecies, or population-specific or may be shared across species.
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
ITS1 PCR
The ITS1 region was amplified using the primers ITS1A (5'-CCTTTGTACACACCGCCCGT) and ITS1D (5'-GTTCATGTGTCCTGCAGTTCAC). This primer pair was originally designed for use on mosquitoes (Sharpe, Harbach, and Butlin 2000) and targets the flanking rRNA gene regions (18S and 5.8S) to produce a PCR product of 602 bp in C. parallelus (GenBank accession number AY 585651). All PCR amplifications were performed on a Hybaid OmniGene thermal cycler and carried out in 25 µl volumes overlaid with one drop of mineral oil. Each 25 µl volume contained: approximately 5 ng genomic DNA template, 3.0 mM MgCl2, 12.5 pmol each primer, 0.12 mM dNTPs, 0.25 units Thermoprime Plus DNA polymerase (ABgene) with the associated ReddyMix PCR buffer at 1x concentration. After an initial denaturation of 5 min at 95°C there were 35 cycles of amplification (1 min at 94°C, 1 min at 67°C, and 2 min at 72°C) before a final incubation of 10 min at 72°C. The products were purified on columns (Promega) prior to cloning or sequencing.
Cloning
Purified PCR products were cloned using the pGEM-T Easy Vector system (Promega) and transformed into JM109 competent cells. Negative controls (water) were used to check for contamination. For each transformation, which represents an individual grasshopper, ITS1 PCR products were generated from 20 randomly chosen clones. A second ITS1 PCR (mismatch PCR) was performed using the same colonies, but the ITS1A primer was replaced by the primer ITS1-4 (5'-GTCTCAAAACTAGCAAACTCGGTATCAT) and the annealing temperature lowered to 55°C. ITS1-4 is a C. parallelus-specific primer designed specifically for studying WIP4. This primer anneals to a region immediately upstream of the WIP4 site but deliberately differs from the DNA template by two mutational changes in order to create a restriction site (Nde I) that encompasses the WIP4 site. Without this modification, the WIP4 site could not be investigated, as no commercial restriction endonuclease is available that could be used to determine the base type at that site.
Sequencing
For each population, a single female grasshopper was chosen from among the sample of 10 individuals. ITS1 PCR products from between five and 11 clones per individual were sequenced in both directions, using the ITS1A and ITS1D primers. Sequence reactions were performed with TaqFS dye-terminator fluorescent chemistry (Applied Biosystems) and run on an Applied Biosystems Model 373 automated sequencer.
Restriction Digests
Digests were performed to determine the base type at the four target ITS1 nucleotide sites (WIP1, 2, 3, and 4). For each clone, three digests were performed using the restriction endonucleases: Bsr I (NEB), Tai I, and Nde I (MBI Fermentas). Following the manufacturer's instructions, 5 µl of PCR product was digested in a total reaction volume of 15 µl. The Tai I and Bsr I digests were performed using PCR products consisting of the entire ITS1 region, whereas the Nde I digest could only be performed on the modified products generated by the mismatch PCR. For all digests there were additional restriction sites within the ITS1 region that were expected to consistently cut, providing a positive control for the reaction. The digested PCR products were visualized on ethidium bromide-stained 2.5% agarose gels. All 15 µl of the digest reaction was loaded, and to ensure accuracy in interpreting the RFLPs, a selection of standards was run in every gel. These standards were composed of digested PCR products whose full ITS1 sequence, and therefore the restriction endonuclease cutting sites, were known. Any fragments whose digest patterns could not be interpreted were sequenced. WIP patterns determined from clones were compared with restriction digests of the original PCR products from selected individuals as an additional check for contamination.
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
In total, nine WIPs were identified among the 30 sequences. Four of these WIPsidentified as WIP1, 2, 3, and 4were subsequently investigated using the PCR-RFLP approach and occur at nucleotide sites 58, 312, 15, and 124 of the C. parallelus' ITS1 sequence, respectively. Intra-individual variation was consistently found, with several haplotypes coexisting in the genome of each individual (fig. 1). Only one WIP (WIP1) was consistently found in all four individuals. From these sequence data the remaining eight WIPs appear to be subspecies-specific. However, the RFLP results show that the nucleotides at the WIP2, 3, and 4 sites are actually polymorphic in both subspecies, although in each case the second base type is relatively rare in one of the subspecies (table 1). This suggests that each mutational change arose prior to the divergence of the two subspecies and that homogeneity has yet to be restored at any of these nucleotide sites after at least 500,000 generations.
|
|
|
|
|
|
|
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
In C. parallelus the populations are evolving cohesively, in that individuals typically share the same ITS1 haplotypes (fig. 3). However, contrary to some predictions (Ohta and Dover 1983), the presence of significant levels of variation among individuals reveals that individuals are not in a similar state at any one point in time (fig. 4). Thus, molecular turnover mechanisms in highly repeated gene families need not preclude the operation of selection. The phenotypic effects of substitutions in ITS1 are unlikely to be sufficient to generate substantial fitness difference between individuals with varying mixes of haplotypes. However, for other repeated loci, including the rDNA coding regions, such differences could well be important.
The high level of variation among C. parallelus individuals indicates that there is substantial heterogeneity in haplotype composition within at least one population of NOR-carrying chromosomes. This means that significant differences typically exist between the maternal and paternal NORs of an individual, with respect to their haplotype distributions. Reciprocal recombination would tend to break down such differences, reducing the heterogeneity in the population. This must be outweighed by the spread of new variants within chromosomes, most likely by intrachromosomal gene conversion. Although gene conversion may only involve a small number of rDNA copies, it has been shown that it can be very effective at altering variant frequencies, especially if it is biased in nature (Hillis et al. 1991).
Ohta and Dover (1984) investigated the cohesive population genetics of molecular drive within multigene families such as rDNA. Their models express the outcome in terms of the variance in copy number relative to the mean copy number in the population of chromosomes. They find that this relative variance (RV) is dependent on the relative rates of recombination and intrachromosomal gene conversion (unequal crossovers were not considered). To compare our data with the predictions of these models, sample estimates of the base type frequencies associated with each WIP (table 1) were used to derive measures of sample relative variance (among individuals). By simulating the sampling strategy employed to generate these data from Cpe populations, the relationship between sample relative variance at the individual level and the underlying relative variance in the population of chromosomes (RV) was established (see Parkin [2002] for details). Our data suggest a population relative variance of approximately 0.1 in terms of variant frequency, with a lower 95% confidence limit of approximately 0.01.
To make comparisons to the Ohta and Dover (1984) models, the RV estimate must be converted to a measure based on numbers of variants, rather than variant frequency, by multiplying the estimate by the number of rDNA copies within the genome. For C. parallelus there is evidence that the rDNA copy number may exceed 1,000 per genome (Parkin 2002), which suggests an extremely high RV value of 100 (lower limit 10). Even with the rDNA copy number reduced to 200, to match Drosophila (Long and Dawid 1980), the resultant RV estimate of 20 (lower limit 2) still exceeds the range of predicted values derived from the models of Ohta and Dover (1984). RV is expected to increase as the rate of intrachromosomal gene conversion increases relative to the rate of recombination. Crossing-over between homologous chromosomes tends to equalize the frequencies of haplotypes among the chromosomes, reducing the variance among individuals. The high estimates of RV in C. parallelus indicate that gene conversion is regularly altering variant frequencies at the intrachromosomal level, and this process more than counteracts the homogenizing effects of interchromosomal recombination. The lowest value of relative variance consistent with our data (RV = 2, n = 200, lower confidence limit of RV) requires the per copy rate of intrachromosomal gene conversion ( in Ohta and Dover [1984]) to be of the same order as the per family rate of interchromosomal recombination per generation ([n 1]ß in Ohta and Dover [1984], where n is the copy number and ß is the per copy recombination rate).
Our approach did not enable us to compare the haplotypes present in the different NORs, except to the extent that this would be revealed from differences between the sexes of Cpp. Variation among NORs is expected because gene conversion is likely to be more common between homologous than nonhomologous chromosomes and has been observed in other species (Copenhaver and Pikaard 1996). However, variation among NORs cannot generate the differences in haplotype distributions that we observe among individuals. This is because each individual contains a sample of each of the NORs (two copies of the L2 NOR and two copies of the L3 NOR in Cpe). Our sample of 20 ITS1 copies from each individual therefore contains a mix of copies from the two NORs, but this mix should not vary among individuals, except due to sampling effects.
Linkage disequilibrium among the ITS1 haplotypes within individuals of C. parallelus also suggests that intrachromosomal gene conversion occurs more frequently than recombination between chromosomes. This is because gene conversion can spread intact haplotypes within repeat arrays, whereas recombination tends to break the association between WIPs. This conclusion was also reached by Schlötterer and Tautz (1994) and Polanco et al. (1998), based on their discovery that the rDNA arrays on the X and Y chromosomes of Drosophila melanogaster are homogenized for different polymorphic alleles. A similar pattern emerged between the two rDNA arrays of Arabidopsis thaliana, with specific polymorphic alleles largely homogeneous in each array. The coexisting polymorphic alleles within these arrays are organized into homogeneous adjacent units (Copenhaver and Pikaard 1996), suggesting that intrachromosomal events spread variants rapidly along an array. Based on the results of this study, the expectation is that the chromosomes of C. parallelus are predominantly populated by different ITS1 haplotypes. This would mean that the ITS1 region is evolving along independent chromosome lineages because the chromosome is, in effect, acting as a barrier to the spread of new mutations.
In C. parallelus significant differences were observed between the ITS1 haplotype distributions of the two subspecies. These differences have evolved in a period of 500,000 generations. Molecular turnover mechanisms have not resulted in rapid divergence because a common haplotype is still shared between subspecies, but new haplotypes have spread gradually over this long time period. Although populations of C. parallelus are now large, the species has a history of range expansion and contraction (Hewitt 1993), so genetic drift may have played a significant role in generating the observed divergence.
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
Richard Thomas, Associate Editor
![]() |
Literature Cited |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Copenhaver, G. P., and C. S. Pikaard. 1996. Two-dimensional RFLP analyses reveal megabase-sized clusters of rRNA gene variants in Arabidopsis thaliana, suggesting local spreading of variants as the mode for gene homogenization during concerted evolution. Plant J. 9:273-282.[CrossRef][ISI][Medline]
Dover, G. 1982. Molecular drivea cohesive mode of species evolution. Nature 299:111-117.[ISI][Medline]
Elder, J. F., and B. J. Turner. 1995. Concerted evolution of repetitive DNA sequences in eukaryotes. Quart. Rev. Biol. 70:297-320.[ISI][Medline]
Flanagan, N. S., P. L. Mason, J. Gosalvez, and G. M. Hewitt. 1999. Chromosomal differentiation through an alpine hybrid zone in the grasshopper Chorthippus parallelus. J. Evol. Biol. 12:577-585.[CrossRef][ISI]
Gosalvez, J., C. Lopez-Fernandez, J. Serrano, J. M. Rubio, and G. M. Hewitt. 1990. Chromosome structure in hybrid zones of Chorthippus parallelus parallelus and Chorthippus parallelus erythropus. Bol. San. Veg. 20:341-348.
Hewitt, G. 1993. After the ice ageparallelus meets erythropus in the Pyrenees. Pp. 140164 in R. G. Harrison, ed. Hybrid zones and the evolutionary process. Oxford University Press, Oxford.
Hillis, D. M., and M. T. Dixon. 1991. Ribosomal DNAmolecular evolution and phylogenetic inference. Quart. Rev. Biol. 66:410-453.[ISI]
Hillis, D. M., C. Moritz, C. A. Porter, and R. J. Baker. 1991. Evidence for biased gene conversion in concerted evolution of ribosomal DNA. Science 251:308-310.[ISI][Medline]
Long, E. O., and I. B. Dawid. 1980. Repeated genes in eukaryotes. Ann. Rev. Biochem. 49:727-764.[CrossRef][ISI][Medline]
Lunt, D. H., K. M. Ibrahim, and G. M. Hewitt. 1998. MtDNA phylogeography and postglacial patterns of subdivision in the meadow grasshopper Chorthippus parallelus. Heredity 80:633-641.[CrossRef][ISI][Medline]
Nagylaki, T., and T. D. Petes. 1982. Intrachromosomal gene conversion and the maintenance of sequence homogeneity among repeated genes. Genetics 100:315-337.
Ohta, T., and G. A. Dover. 1983. Population genetics of multigene families that are dispersed into 2 or more chromosomes. Proc. Natl. Acad. Sci. USA 80:4079-4083.[Abstract]
Ohta, T., and G. A. Dover. 1984. The cohesive population genetics of molecular drive. Genetics 108:501-521.
Parkin, E. J. 2002. Concerted evolution of ribosomal DNA loci in Chorthippus parallelus. Doctoral thesis, University of Leeds, Leeds, United Kingdom.
Polanco, C., A. I. Gonzalez, A. de la Fuente, and G. A. Dover. 1998. Multigene family of ribosomal DNA in Drosophila melanogaster reveals contrasting patterns of homogenization for IGS and ITS spacer regions: A possible mechanism to resolve this paradox. Genetics 149:243-256.
Reeder, R. H. 1990. Ribosomal RNA synthesis in the nucleolus. Trends Genet. 6:390-395.[CrossRef][ISI][Medline]
Reynolds, W. J. 1980. A re-examination of the characters separating Chorthippus montanus and C. parallelus (Orthoptera: Acrididae). J. Nat. Hist. 14:283-303.[ISI]
Schlötterer, C., M. T. Hauser, A. Vonhaeseler, and D. Tautz. 1994. Comparative evolutionary analysis of rDNA ITS regions in Drosophila. Mol. Biol. Evol. 11:513-522.[Abstract]
Schlötterer, C., and D. Tautz. 1994. Chromosomal homogeneity of Drosophila ribosomal DNA arrays suggests intrachromosomal exchanges drive concerted evolution. Curr. Biol. 4:777-783.[ISI][Medline]
Schneider, S., D. Roessli, and L. Excoffier. 2000. Arlequin v. 2.000: a software for population genetics data analysis.
Sharpe, R. G., R. E. Harbach, and R. K. Butlin. 2000. Molecular variation and phylogeny of members of the Minimus group of Anopheles subgenus Cellia (Diptera : Culicidae). Syst. Entomol. 25:263-272.[CrossRef][ISI]
Stephan, W. 1989. Tandem-repetitive non-coding DNA: forms and forces. Mol. Biol. Evol. 6:198-212.[Abstract]
Sunnucks, P, and D. F. Hales. 1996. Numerous transposed sequences of mitochondrial cytochrome oxidase I-II in aphids of the genus Sitobion (Hemiptera: Aphididae). Mol. Biol. Evol. 13:510-524.[Abstract]
Tautz, D., J. M. Hancock, D. A. Webb, C. Tautz, and G. A. Dover. 1988. Complete sequences of the ribosomal RNA genes of Drosophila melanogaster. Mol. Biol. Evol. 5:366-376.[Abstract]
Thompson, J. D., T. J. Gibson, F. Plewniak, F. Jeanmougin, and D. G. Higgins. 1997. The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 25:4876-4882.
Walsh, J. B. 1987. Persistence of tandem arrays: Implications for satellite and simple sequence DNAs. Genetics 115:553-567.
Weir, B. S. 1996. Genetic data analysis II. Sinauer Associates, Sunderland, Mass.