Retrotransposons and Genomic Stability in Populations of the Young Allopolyploid Species Spartina anglica C.E. Hubbard (Poaceae)

Alex Baumel*, Malika Ainouche*, Ruslan Kalendar{dagger} and Alan H. Schulman{dagger}{ddagger}

*Population, Species and Evolution Lab., UMR CNRS 6553, University of Rennes 1. Campus Scientifique de Beaulieu, Rennes Cedex, France;
{dagger}Plant Genomics Laboratory, Institute of Biotechnology, University of Helsinki;
{ddagger}MTT Agrifood Research Finland, Jokioinen, Finland


    Abstract
 TOP
 Abstract
 Introduction
 Plant Material
 Methods
 Results
 Discussion
 Acknowledgements
 References
 
Spartina x townsendii arose during the end of the 19th century in England by hybridization between the indigenous Spartina maritima and the introduced Spartina alterniflora, native to the eastern seaboard of North America. Duplication of the hybrid genome gave rise to Spartina anglica, a vigorous allopolyploid involved in natural and artificial invasions on several continents. This system allows investigation of the early evolutionary changes that accompany stabilization of new allopolyploid species. Because allopolyploidy may be a genomic shock, eliciting retroelement insertional activity, we examined whether retrotransposons present in the parental species have been activated in the genome of S. anglica. For this purpose we used inter-retrotransposon amplified polymorphism (IRAP) and retrotransposons-microsatellite amplified polymorphism (REMAP) markers, which are multilocus PCR-based methods detecting retrotransposon integration events in the genome. IRAP and REMAP allowed the screening of insertional polymorphisms in populations of S. anglica. The populations are composed mainly of one major multilocus genotype, identical to the first-generation hybrid S. x townsendii. Few new integration sites were encountered in the young allopolyploid genome. We also found strict additivity of the parental subgenomes in the allopolyploid. Both these findings indicate that the genome of S. anglica has not undergone extensive changes since its formation. This contrasts with previous results from the literature, which report rapid structural changes in experimentally resynthesized allopolyploids.


    Introduction
 TOP
 Abstract
 Introduction
 Plant Material
 Methods
 Results
 Discussion
 Acknowledgements
 References
 
Allopolyploidy is a well-known process in plant speciation involving the duplication of an interspecific hybrid genome. However, evolution of the two duplicated parental genomes (i.e., the homoeologous subgenomes) in the allopolyploid remains obscure in most cases because of the lack of information on the initial conditions (Wendel 2000Citation ). In most natural polyploids, the parental species have evolved or become extinct since the hybridization event. Spartina anglica (Poaceae) is among the few well known recent allopolyploids. The history of this species has been well documented on the basis of morphological (Marchant 1967Citation ; Hubbard 1968Citation , pp. 358–359), cytogenetic (Marchant 1968Citation ), and isozyme (Guénégou, Citharel, and Levasseur 1988Citation ; Gray, Benham, and Raybould 1990Citation ; Raybould et al. 1991a,Citation 1991bCitation ) data. Spartina anglica arose around 1890 in southern England (Southampton Water), after hybridization between the indigenous S. maritima Fernald (2n = 60) and the introduced Eastern North American S. alterniflora Loiseleur (2n = 62). Chromosome doubling in the sterile hybrid (Spartina x townsendii) gave rise to the new allopolyploid species S. anglica (2n = 120, 122, and 124, according to Marchant 1968Citation ). This species is a fertile, perennial plant that has been actively colonizing the salt marshes and estuaries of Britain and France since 1890 and 1906, respectively. Both significant seed production and extensive lateral clonal growth characterize this rhizomatous species. As a nascent species (less than 150 years old), S. anglica is a particularly suitable model to trace allopolyploid genome evolution because its parental species are unambiguously identified and are both still extant.

Previous studies have shown no variation at isozyme loci in British (Raybould et al. 1991aCitation ) and French populations (M. T. Misset, personal communication). A recent study of several French populations (Baumel, Ainouche, and Levasseur 2001Citation ), using RAPD markers, also revealed low molecular variation, thus supporting a unique origin of the allopolyploid (i.e., unique hybridization event or multiple events involving identical parental genotypes). Both parental species (S. alterniflora and S. maritima) lack genetic variation in western Europe (Raybould et al. 1991a,Citation 1991b;Citation Yannic, Baumel, and Ainouche, unpublished data).

In this study, we focus on potentially rapidly evolving parts of the S. anglica genome. Because changes can occur at both single and multiple loci, dispersed components of the genome must be investigated to test the null hypothesis of strict parental genome additivity in the allopolyploid. It is well known that retroelements constitute an important fraction of the DNA content of many plants, particularly in grasses (Vicient et al. 2001Citation ). Their abundance, dispersion across the nuclear genome, and their insertional activity indicate that they play a significant role in plant genome evolution (Bennetzen 2000Citation ). Retrotransposons propagate by a replicational mechanism that has potential to remodel the genome (Kumar and Bennetzen 1999Citation ). In wild barley, Hordeum spontaneum, the retrotransposon BARE-1 generates genetic polymorphism and increases genome size by its activity (Kalendar et al. 2000Citation ). Studies on BARE-1 in Hordeum indicate that it is an active component of genome size variation in Hordeum species, which appears prone to adaptation to environmental stresses such as drought (Vicient et al. 1999Citation ; Kalendar et al. 2000Citation ; Wendel and Wessler 2000Citation ).

Retrotransposon activity is regulated by various factors, and biotic and abiotic stresses can enhance transposition in natural populations (Wessler 1996Citation ; Kalendar et al. 2000Citation ). Although it is known that transposable elements can be involved in speciation, as suggested by the example of hybrid dysgenesis in Drosophila (Kidwell and Lisch 2000Citation ), little is known about their role during the first stages of homoploid hybrid and allopolyploid speciation in plants. Interspecific crosses could result in genome stress (McClintock 1984Citation ; O'Neill, O'Neill, and Marshall Graves 1998Citation ; Liu and Wendel 2000Citation ), increasing retroelement insertional activity. Because polyploid genomes contain duplicates of all genes, they are thought to be buffered from the deleterious consequences of transposition (Matzke, Mittelsten Scheid, and Matzke 1999Citation ).

The high copy number and general dispersion of retrotransposons throughout the genome, as well as the large local change caused by their insertion, provide an excellent basis for the development of DNA-based markers. Such markers allow investigation of genetic variation related to the distribution and activity of retrotransposons, both at the intraspecific and the interspecific levels. Retrotransposon marker systems rely on PCR to generate fingerprints. The fingerprints are multilocus profiles revealing polymorphisms in the insertion of members of given families of retrotransposons. The inter-retrotransposon amplified polymorphism (IRAP) technique (Kalendar et al. 1999Citation ) generates PCR products from retrotransposons inserted near enough to each other so as to allow efficient amplification and generally uses primers matching the outer segments of the long terminal repeats (LTRs). Retrotransposons-microsatellite amplified polymorphism (REMAP) (Kalendar et al. 1999Citation ) uses one LTR primer, together with a primer designed for annealing to the 3' end of a stretch of simple sequence repeats (SSR), and detects retrotransposons inserted near to SSRs. Therefore, REMAP can amplify a pool of three different sequences: fragments situated between an LTR and a microsatellite locus, sequences situated between two microsatellite loci, and sequences situated between retroelements. Because many retroelements have the tendency to insert into regions rich in tandemly repeated DNA (Kalendar et al. 1999Citation ; Ramsay et al. 1999Citation ), as well as near to, or nested within, each other (San Miguel et al. 1996Citation ; Shirasu et al. 2000Citation ), REMAP and IRAP are efficient at detecting insertion events. In this study, fragments situated between microsatellite loci (Inter Simple Sequence Repeat or ISSR markers, Wolfe and Liston 1998Citation ) were also amplified in order to detect eventual variation affecting the repetitive fraction of the allopolyploid genome.

The goal of this study was to determine whether polyploid formation has led to retrotransposon activation, generating molecular variation in S. anglica populations. Because this clonal plant is characterized by vigorous vegetative propagation, mutations that have occurred since its appearance are likely to be fixed in populations. Our previous studies revealed a striking lack of genetic variation (Baumel, Ainouche, and Levasseur 2001Citation ), indicating that almost all the French populations of S. anglica are composed of the same major multilocus genotype. In this study, several samples from the center of origin of the species (Southampton salt marshes, England) have been analyzed by IRAP, REMAP, and ISSR markers. The plants examined include both the parental species (S. maritima and S. anglica) and the first-generation hybrid S. x townsendii.


    Plant Material
 TOP
 Abstract
 Introduction
 Plant Material
 Methods
 Results
 Discussion
 Acknowledgements
 References
 
A list of the Spartina samples analyzed in this study is presented in table 1 . The samples of S. x townsendii are from the Hythe (UK) salt marshes, where this species was recorded for the first time in 1870. These samples were identified as S. x townsendii by morphological observation (sterile stamens) and by genome size estimated with flow cytometry (unpublished data). Four populations of S. anglica have been sampled near this site, in Falley and Keyhaven, separated by approximately 5 and 15 km, respectively, from Hythe. Studland is in the southeastern part of Pool Harbour, about 50 km from Hythe. The fourth population is at Bull Island in Ireland. Individual plants were collected in the field with roots, leaves, and rhizomes to allow both for DNA extraction (from leaves) and for transplanting into the greenhouse. Because Spartina species are present at each site as clonal plants, we tried to maximize the distance between samples to avoid collecting from the same individual. Ten individuals per population were then collected for S. anglica and S. x townsendii. The parental species, S. maritima and S. alterniflora, were included for comparison. Three samples of S. maritima have been collected in different French salt marshes, and S. alterniflora has been collected in Marchwood (near Hythe, UK). As mentioned previously, these species lack genetic variation in western Europe; then, we might reasonably expect this sampling to be representative from the actual progenitors of S. anglica.


View this table:
[in this window]
[in a new window]
 
Table 1 List and Origin of the Plant Material (* = Chromosome Numbers from Marchant 1968)

 

    Methods
 TOP
 Abstract
 Introduction
 Plant Material
 Methods
 Results
 Discussion
 Acknowledgements
 References
 
DNA Extraction
Samples of DNA were prepared by the Cetyltrimethylammonium bromide (CTAB) method (Ausubel et al. 1995Citation , pp. 4.7.1–4.7.8). Approximately 30 mg of fresh leaves were ground in 2% CTAB solution (1.4 M NaCl, 20 mM EDTA, 100 mM Tris-HCl pH 8.0, 2% CTAB). The extract was incubated for 1 h at 65°C, then 500 µl of chloroform were added. After mixing, 500 µl of 100% isopropanol was added to the supernatant for precipitation. After microcentrifugation, the DNA pellet was washed with 70% ethanol, then resuspended in the Tris-EDTA buffer.

IRAP and REMAP Procedures
The IRAP and REMAP procedures were based on earlier published methods (Kalendar et al. 1999Citation ; Schulman et al. 2002). The IRAP reactions of 20 µl contained: 0.075 M Tris-HCl pH 8.8, 0.02 M (NH4)2SO4, 1.5 mM MgCl2, 0.01% Tween-20, 0.2 mM dNTPs, 0.2 µM each primer (presented in table 2 ), 1.5 U Taq polymerase, and 10 ng DNA template. The amplifications were performed on a Mastercycler Gradient (Eppendorf) using the following program: 2 min at 94°C; 33 cycles of 30 s at 94°C, 30 s at 58°C, 2 min at 72°C; termination by 10 min of final extension at 72°C.


View this table:
[in this window]
[in a new window]
 
Table 2 Primers

 
The reaction conditions and amplification program for REMAP were similar to those of IRAP, except that we used less Taq1 polymerase (1 U per reaction) and fewer PCR cycles (29 instead of 33).

Inter Simple Sequence Repeat Procedure
ISSRs are PCR markers commonly displaying a high degree of polymorphism (Wolfe and Liston 1998Citation ). In order to determine the specific contribution of retrotransposons to the observed polymorphism, we performed ISSR amplifications on the same Spartina material used for IRAP and REMAP. We used the following SSR primers: (GTG)7C, (CTC)6G, (TGC)6C, (AGC)6G. Amplification conditions were the same as for IRAP but with only 30 cycles.

Electrophoresis and Data Analysis
The three fingerprinting techniques produce numerous bands. Therefore, high-resolution electrophoresis is required for their reaction products. PCR products were run on 2% agarose gels (RESULT LE, BIOzym, Landgraaf, The Netherlands) for 7–8 h at 85 volts, stained with ethidium bromide, and photographed with a 35-mm film-based reflex camera.

As some primer combinations have one primer in common, these combinations might amplify the same polymorphic band. This eventual bias in polymorphism estimates was checked by band size comparisons and by recording the absence or presence of these bands among the samples.

IRAP and REMAP data were analyzed in two ways. First, a table scoring the presence or absence of the bands observed on the gels was created in order to record the variation in S. anglica populations. Invariable, ambiguous and nonreproducible bands were not included. Second, reaction products for the allopolyploid and its parents were run on the same gel, and interspecific comparisons were performed. Under the assumption that bands of the same size are derived from the same loci in each of these species, seven possible outcomes could be scored, with a given band being present or absent in one, two, or three species.


    Results
 TOP
 Abstract
 Introduction
 Plant Material
 Methods
 Results
 Discussion
 Acknowledgements
 References
 
Adaptation of IRAP and REMAP to Spartina
In the absence of characterized native retrotransposons, primers matching retrotransposons from other grass genera were tested for use in Spartina. The approach seemed reasonable because active retrotransposon families appear to cross generic boundaries in the grasses (Vicient et al. 2001Citation ). Several LTR primers designed for barley were tested in single-primer amplifications. Only those producing numerous, strong bands and a weak background, as well as polymorphism between the parental species, were selected for IRAP and REMAP. Because the annealing conditions of IRAP, REMAP, and ISSR were stringent, reproducible band patterns were expected. Reproducibility was tested and polymorphism confirmed by multiple reactions and electrophoretic separations.

Retrotransposon primers used in this study are listed in table 2 and match the LTRs of different families of elements. Two primers were designed to anneal to sequences of the miniature transposable elements Sire and Stowaway. Retrotransposon primers were also tested as pairs for IRAP and, in combination with SSR primers, for REMAP (table 2 ). We selected the following LTR/LTR primers for IRAP: C0795/C0945, C0795/Sire, C0795/Stowaway, C0699/C0945, C0947/C0905, 9900/C0905, 9900/C0699. For REMAP, the following LTR/ISSR primers were used: C0795/(GAG)6C, C0795/(GTG)7A, 9900/(GTG)7A, 9900/(GTG)7C.

Genetic Variation in S. anglica Populations Examined by IRAP, REMAP, and ISSR
Forty samples of S. anglica and 10 samples of S. x townsendii were analyzed with eight IRAP and four REMAP primer combinations and with four SSR primers. Over 296 bands were generated by IRAP and REMAP markers in S. anglica. IRAP produced an average of 26 bands for each primer combination, whereas REMAP yielded an average of 22 bands per primer combination. The ISSR procedure generated 92 bands with an average of 23 bands per primer.

In S. anglica, we found seven polymorphic bands with IRAP, three with REMAP, and three with ISSR analysis. The primer combinations that revealed polymorphism are presented in table 3 . These 13 polymorphic bands constituted 14 multilocus genotypes (table 3 ). Most samples (40%) belonged to one genotype. The other 13 genotypes were both of low frequency and very similar to the major genotype because the differences involve no more than absence or presence of 1-4 bands. Two IRAP bands were polymorphic in S. x townsendii (Hythe population), two samples each lacking one band. Except for these two samples, genotypes of the other eight accessions of S. x townsendii are similar to the major genotype of S. anglica. This major genotype is additive with respect to the parental bands. We assume that the major genotype is the founder of both hybrid S. x townsendii and S. anglica and that therefore the variation observed in these species is the result of mutations.


View this table:
[in this window]
[in a new window]
 
Table 3 IRAP, REMAP and ISSR Multilocus Genotypes and Their Distribution in Spartina anglica Populations. The Last Column Indicates the Parental Origin of the Marker When it was Possible to Identify it: Mar for S. maritima, Alt for S. alterniflora and Mar/Alt for both; Ang Indicates that the Marker is Specific to the Allopolyploid. N = Number of Samples Recorded for Each Genotype, M = Number of Individuals Recorded for Each Mutation Relative to the Major Genotype

 
Most polymorphisms (8/13) were restricted to one population, whereas five were recorded in more than one population. The two most frequent (C0947/C0905-C; C0795/[GAG]6A-A) occurred in 17.5% and 20% of the samples, respectively, and are each found in two or three populations. Furthermore, excepting the major genotype, only two multilocus genotypes (6 and 13) were recorded in more than one population. Of these, genotype 6 was found at Bull Island and Studland, and genotype 13 was recorded at Keyhaven and Fawley. Of the polymorphic bands in S. anglica (table 3 ), 9 out of 13 represent band absence from the major genotype. Six belong to the S. maritima genome, two are specific to the S. alterniflora genome, and one is present in both parents. Only three new IRAP bands and one new ISSR band were observed in the allopolyploid S. anglica. Two out of three new IRAP bands occurred only in one individual from Keyhaven. These two IRAP bands could be amplified by most of the combinations that involve the BARE-1 and Sukkula LTR primers (C0795 and 9900, respectively). Their reproducibility was confirmed by independent DNA preparations from another leaf. These DNA fragments were robustly amplified (fig. 1 ).



View larger version (113K):
[in this window]
[in a new window]
 
Fig. 1.—Electrophoretic pattern generated by IRAP, using primers 9900 and C0699. Lanes 1–10: S. anglica from Keyhaven, lane 11: S. maritima, lane 12: S. alterniflora, lanes 13–22: S. anglica from Studland, lanes 23–32: S. anglica from Bull Island

 
We cloned and sequenced these two polymorphic IRAP bands and determined their identity by BLAST searching. Band K2 (accession AF 450249), amplified by the primer 9900, contains a region (nt 325–959) comprising about 40% of a the LTR of a BARE-1 retrotransposon (accession Z17327), as well as a region (nt 26–324) which may be protein-encoding on the basis of EST and protein database matches. The structure of this novel band indicates that at least two, nested BARE-1 integration events, as seen earlier in the barley and maize genomes (San Miguel et al. 1996Citation ; Shirasu et al. 2000Citation ), occurred to generate the amplifiable domain. The other sequence, B2 (accession AF 450250), was generated with Sukkula and Stowaway primers and contains an internal domain which may encode part of a proline-rich protein.

Comparison Between S. anglica and its Parents
IRAP, REMAP, and ISSR patterns were compared among the species by recording the presence or absence of PCR bands (table 4 ). The comparison was based on the major genotype of S. anglica, which is identical to that of S. x townsendii. The three data sets were congruent and produced the same result. Most markers of S. maritima and S. alterniflora were species-specific, and only a minority (15%) was found in both species. The three techniques, but particularly IRAP and REMAP, indicated that S. maritima contributed slightly more markers (41%) in the allopolyploid than S. alterniflora (34.5%). The first three combinations of markers, representing the highest share (90.5%) of the total number of bands, result from the additivity of the parental markers in the allopolyploid (table 4 ).


View this table:
[in this window]
[in a new window]
 
Table 4 Presence-Absence of Markers in the Three Spartina Species, Alt = S. alterniflora, Mar = S. maritima, Ang = S. anglica. Comparison is Based on 224 IRAP Markers, 100 REMAP Markers, and 108 ISSR Markers

 
A minority of the marker bands (9.5%) was found only in one of the parental species. The major genotype does not exhibit any marker specific to S. anglica (i.e., not present in the parental species).


    Discussion
 TOP
 Abstract
 Introduction
 Plant Material
 Methods
 Results
 Discussion
 Acknowledgements
 References
 
LTR Primer Specificity
The IRAP and REMAP methods were first implemented in Hordeum (Kalendar et al. 1999Citation ). They were subsequently applied in Hordeum vulgare for fingerprinting and biodiversity studies (Kalendar et al. 2000Citation ; Vicient et al. 2001Citation ) and for mapping of traits in H. vulgare (Manninen et al. 2000Citation ) and Aegilops tauschii (Boyko et al. 2002Citation ). Within the Poaceae, Spartina belongs to the subfamily Chlorodoideae, and is not closely related to Hordeum, which is in the Pooideae subfamily (Hsiao et al. 1999Citation ; Kellog 2001Citation ). Nevertheless, most of the primers designed in Hordeum amplify bands in Spartina. This result is consistent with the evidence that the closely related retrotransposons are transcribed and translated in diverse grass species belonging to different subfamilies (Vicient et al. 2001Citation ). The primers produced reproducible patterns (fig. 1 ) and a greater number of bands than seen with RAPDs (Baumel, Ainouche, and Levasseur 2001Citation ) or ISSRs (this study). The abundance of marker bands indicates that retroelements closely related to many of those in barley are abundant in Spartina and display a similar genome organization (clustering, association with SSRs), similar to the pattern identified in barley (Shirasu et al. 2000Citation ).

IRAP and REMAP in S. anglica Populations
Fourteen genotypes were identified in the populations of S. anglica on the basis of 13 polymorphic bands. One genotype is most frequent in all populations, representing 40% of all samples. Most populations are composed of the same major genotype, identical to that of the first-generation hybrid S. x townsendii. The other genotypes encountered in this study differ only by the absence of 1-4 markers or, very rarely, by the appearance of a new band. Previous investigations in populations of S. anglica using multilocus RAPD markers (Baumel, Ainouche, and Levasseur 2001Citation ) also revealed a general lack of genomic variation. Taken together, these results support a dramatic genetic bottleneck at the time of formation of the species in Southampton, perhaps as a consequence of the unique origin of the species.

Evolution of the Parental Repetitive Sequences in the Allopolyploid
Both ISSR and the retrotransposon-based IRAP and REMAP reveal numerous markers that differentiate the parental genomes S. alterniflora and S. maritima as do the RAPD data. (Baumel, Ainouche, and Levasseur 2001Citation ). Our recent phylogenetic analysis (Baumel et al. 2002Citation ) indicates that although S. alterniflora and S. maritima belong to the same lineage, they contain divergent nuclear and chloroplast DNA sequences. This ability to differentiate the two parental genomes on the basis of several loci provides the opportunity to compare them with their homologues duplicated in the allopolyploid. A comparison of the three data sets (IRAP, REMAP, and ISSR) between S. anglica and its parents shows that the parental genomes remain largely unchanged in the allopolyploid, 90.5% of the markers showing additivity for the parental genomes. The remaining 9.5% of parental markers, found neither in the allopolyploid nor in S. x townsendii, might result from polymorphism in the parental species (i.e., markers not present in the genotypes involved in the parentage of S. x townsendii), or, alternatively, might result from parental marker loss in the hybrid and the allopolyploid, as reported in Brassica (Song et al. 1995Citation ) and wheat (Liu et al. 1998a,Citation 1998b;Citation Ozkan, Levy, and Feldman 2001Citation ), though a particular mechanism for these losses has to be elucidated. Aside from cryptic parental polymorphism and PCR artifacts, marker loss could occur from recombination, mutation, or deletion events, that separate or destroy PCR amplification sites. Because our analyses reveal very few mutations affecting the initial (major) genotype of S. anglica, we can conclude that the young allopolyploid genome has not undergone extensive change since its formation. Spartina anglica populations therefore are composed of the initial genotype that, in some samples, displays a few mutations that do not extensively affect the allopolyploid genome as a whole.

Nondeleterious mutations may be maintained in natural populations as a result of the vigorous clonal reproduction of S. anglica. Mutations may progressively lead to divergence between the parental (homoeologous) subgenomes present in the allopolyploid and their homologues in the parental species. Currently, however, such divergence is limited. Our multilocus investigations reject the hypothesis of rapid structural changes in S. anglica subgenomes. The results presented here and those published previously (Baumel, Ainouche, and Levasseur 2001Citation ) in the 150 year old S. anglica contrast with the rapid (1-3 generations old) and extensive changes reported in synthetic Brassica and Triticum allopolyploids (Song et al. 1995Citation ; Feldman et al. 1997Citation ; Liu et al. 1998a,Citation 1998b;Citation Ozkan, Levy, and Feldman 2001Citation ; Shaked et al. 2001Citation ). These authors demonstrate that rapid, non-Mendelian changes, involving preferential sequence elimination or modification of methylation patterns, may occur in the earliest stages of polyploid stabilization. As a consequence, young polyploids generally are not assumed to be completely additive with respect to their progenitors (Wendel 2000Citation ). However, analysis of the young allopolyploid B. juncea revealed no structural changes when compared with its diploid progenitors (Axelson et al. 2000Citation ), in parallel with our data for Spartina. Furthermore, both genomic additivity and epigenetic stasis were found in artificial (resynthesized) Gossypium polyploids (Liu et al. 2001Citation ).

In S. anglica, 9 of 13 polymorphisms result from the lack of parental bands from one of the parents. Among these, six are from S. maritima and two are from S. alterniflora. One band is present in S. maritima, S. alterniflora, the hybrid S. x townsendii, and present or absent in S. anglica. The multiplex REMAP and IRAP methods, unlike the monoplex RBIP technique (Flavell et al. 1998Citation ), result in dominant markers, which means that the allele lacking the retrotransposon insertion cannot be directly scored. The assumption that a missing band means that an empty allele has been inherited, rather than that a mutation has eliminated PCR primer binding, has not been systematically examined and requires isolation and sequencing of the loci from the genome of the polyploid. However, previous extensive sequencing of the LTR termini (Suoniemi, Schmidt, and Schulman 1997Citation ), in which REMAP and IRAP primers lie, demonstrated that mutations are rare in these regions critical to retrotransposon function and hence to band polymorphism. However, the one band in S. anglica that is present in both parents but missing in a few allopolyploid individuals may be caused by such a mutation. In a previous study (Baumel, Ainouche, and Levasseur 2001Citation ), six of seven RAPD polymorphic bands also originated from S. maritima. On the basis of a limited number of polymorphic loci (16), this asymmetrical parental subgenome variation may indicate that the mutation rate is higher in the S. maritima subgenome of the allopolyploid S. anglica than in the S. alterniflora subgenome. Because S. alterniflora has been shown to be the maternal parent of S. anglica (Ferris, King, and Gray 1997Citation ), this would indicate that the paternal genome is preferentially modified. Similar findings have been reported in the young synthetic allopolyploids involving reciprocal crosses between the diploids Brassica rapa and Brassica nigra and between B. rapa and Brassica oleracea (Song et al. 1995Citation ).

Genome duplication in polyploids may buffer the potentially deleterious effects of mutations that might result from bursts of retroelement activity (Matzke and Matzke 1998Citation ). The genome stress of interspecific hybridization, in turn, may activate such bursts (McClintock 1984Citation ). We therefore expected that new retroelement insertions might be detectable in the young allopolyploid genome of S. anglica. In the allotetraploid G. hirsutum, Zhao et al. (1998)Citation found that colonization of one subgenome by repetitive elements from the other subgenome followed the polyploidization event of about a million years ago. Overall, however, this genome showed remarkably few changes after the formation of the artificial Gossypium polyploids (Liu et al. 2001Citation ). The data of Liu and Wendel (2000)Citation indicate that, after genetic introgression in Zizania latifolia, retroelement activity is ephemeral and rapidly repressed.

Of particular interest are the two new IRAP bands detected in the population from Keyhaven (sample K2), each generated by Sukkula or BARE-1 LTR primers. These new bands correspond to integrations of a retrotransposons, which likely have occurred since the origin of the allopolyploid. The integration of new retrotransposon copies indicates that the BARE-1 family of retrotransposons remains active in S. anglica. The two sites are 0.68% of the total number of retrotransposon insertion sites (296) scored in S. anglica. Liu and Wendel (2000)Citation demonstrated retrotransposon activation in the rice-Z. latifolia introgression lines that they had constructed. To our knowledge, though, the new IRAP bands found in S. anglica represent the first new retrotransposon integration events documented to have occurred in a natural population during a defined period of time, in this case less than 150 years. New integration events are, however, often observed in plants regenerated from tissue culture (Hirochika et al. 1996Citation ), which is known to activate retrotransposons and which served in the isolation of the first such active element from plants (Grandbastien, Spielmann, and Caboche 1989Citation ).

On the basis of these data, one can attempt to estimate the number of integration events represented by these two new insertions. Spartina anglica has a base chromosome number of 10, as do other species in the Chloridoideae (Marchant 1968Citation ); given a genome size of 5.341 x 109 bp (2N = 122; http://www.rbgkew.org.uk/cval/homepage.html), the basic haploid genome is 4.378 x 108 bp or about twice that of rice (2.45 x 108). The Arabidopsis thaliana genome, the only plant genome for which complete information exists (Bevan et al. 1998Citation ), contains 3 x 104 genes of 4.8 kb each. Assuming a similar density in other plants (data on 1,000 genes indicates a slightly lower density in rice), genes in the S. anglica genome occupy about 1.44 x 108 bp or roughly one-third of the genome, compared with perhaps 3% in barley (Vicient et al. 2001Citation ). Hence, the repetitive fraction would constitute 4.36 x 109 bp in barley but 2.94 x 108 bp in the S. anglica basic genome. Barley contains about 8.7 x 104 BARE-1 LTRs (Vicient et al. 1999Citation ), the priming sites for IRAP and REMAP based on BARE-1 and about the same number of Sukkula elements (unpublished data). If the number of BARE-1 and Sukkula elements are reduced proportionately according to the fraction of repetitive DNA in the genome, then in S. anglica, the basic (haploid) genome might contain about 6,000 BARE-1 or Sukkula LTRs and the polyploid genome about 7.1 x 104. All marker combinations involving BARE-1 or Sukkula primers amplified a total of 154 bands from S. anglica. If all the BARE-1 and Sukkula LTRs could generate PCR products, one would expect 6,000 products from each haploid genome. The 154 bands observed represent 2.6% (154/6,000) of the total expected if the each haploid (x) chromosome set were to give products of identical size (because of the retroelements being identically organized in each genome). If, however, the constituent genomes are fully distinct regarding their retrotransposon organization, 154 bands would represent only 0.02% of the total possible (154/7.1 x 104).

Taking these calculations on retrotransposon sampling efficiency, repetitive DNA fraction, and retrotransposon copy number into account, two insertion events would represent between 76 (2/0.026) and 929 (2/0.0022) insertion events for these two retrotransposon families. This implies that the two novel insertions detected may represent between 75 and 1,000 retrotransposon integration events in the less than 150 years since the formation of the S. anglica genome. Although this suggests that the complement of BARE-1 and Sukkula elements grew by at most 1% in the S. anglica genome during the speciation event, it is nevertheless consistent with the hypothesis of activation by polyploidization, but such ploidy-specific variation, if any, has changed the genome of the allopolyploid very little. These calculations, in addition to being subject to sampling error, may be underestimates of retrotransposon activity in the S. anglica genome for two reasons: first, ancient element families, conserved between relatively distant genera, may be less active than more recently evolving, native elements; furthermore, nested integration events such as seen in other grasses (San Miguel et al. 1996Citation ; Shirasu et al. 2000Citation ) might not have been detected.

Recent progress in molecular analysis of allopolyploid systems resulted in significant insights into polyploid genome evolution (Wendel 2000Citation ). One of the best-investigated models in this respect is genus Gossypium. Contrasting results have been derived from different portions of this genome, revealing independent evolution of the parental subgenomes for low- and single-copy sequences (Cronn, Small, and Wendel 1999Citation ; Small and Wendel 2000Citation ) but nonindependent and concerted evolution of repetitive sequences (Wendel, Schnabel, and Seelanan 1995Citation ; Zhao et al. 1998Citation ). As mentioned previously, insights from experimentally resynthesized allopolyploids also reveal conflicting results: either genomic or phenotypic instability after polyploidization (Song et al. 1995Citation ; Feldman et al. 1997Citation ; Liu et al. 1998a,Citation 1998b;Citation Comai et al. 2000Citation ) or genomic stasis (Axelson et al. 2000Citation ; Liu et al. 2001Citation ). All natural alloplyploid models investigated to date are more or less stabilized species. Insights from early stages of the polyploid species formation have been exclusively recorded from studies of experimentally created polyploids. Spartina anglica is to our knowledge the first natural nascent species investigated in this respect. Our results indicate overall additivity and overall structural stasis of the parental subgenomes in this allopolyploid, concerning both low-copy (Baumel, Ainouche, and Levasseur 2001Citation ) and repetitive (this study) sequences, although some changes might be encountered in local populations, as noticed in Keyhaven. Comparisons among different natural young systems are needed for a better understanding of the short-term evolutionary process after duplication of a hybrid genome under natural selection pressures.


    Acknowledgements
 TOP
 Abstract
 Introduction
 Plant Material
 Methods
 Results
 Discussion
 Acknowledgements
 References
 
This work has been supported by grants from the European Molecular Biology Organization (to A.B.), from the CNRS of France (UMR Ecobio 6553, University of Rennes 1 to M.A.), by contract QLK5-CT-1999-0149 to the Commission of the European Communities (to R.K.), and by the Academy of Finland (to A.S.). The authors are grateful to A. J. Gray and A. F. Raybould (Centre for Ecology and Hydrology, Dorset, UK) for their kind assistance and helpful advice on Spartina sampling in Southern England. M. Mc Corry (University College, Dublin) is thanked for providing Spartina material from Ireland. The authors deeply appreciated the constructive comments of two anonymous reviewers who helped to improve the manuscript.


    Footnotes
 
Brandon Gaut, Reviewing Editor

Keywords: allopolyploidy retrotransposon IRAP REMAP Spartina anglica genome stress genomic stability Back

Address for correspondence and reprints: Malika Ainouche, Population, Species and Evolution Lab., UMR CNRS 6553, University of Rennes 1. Campus Scientifique de Beaulieu, 35042 Rennes Cedex, France. malika.ainouche{at}univ-rennes1.fr Back


    References
 TOP
 Abstract
 Introduction
 Plant Material
 Methods
 Results
 Discussion
 Acknowledgements
 References
 

    Ausubel F. M., R. Brent, R. E. Kingston, D. D. Moore, J. G. Seidman, J. A. Smith, K. Struhl, L. M. Albright, D. M. Coen, A. Varki, 1995 Current protocols in molecular biology John Wiley & Sons, New York

    Axelson T., C. W. Bowman, A. G. Sharpe, D. J. Lydiate, U. Lagercranz, 2000 Amphidiploid Brassica juncea contains conserved progenitor genomes Genome 43:679-688[ISI][Medline]

    Baumel A., M. L. Ainouche, R. J. Bayer, A. K. Ainouche, M. T. Misset, 2002 Molecular phylogeny of hybridizing species from the genus Spartina Schreb. (Poaceae) Mol. Phylogenet. Evol 22:303-314[ISI][Medline]

    Baumel A., M. L. Ainouche, J. E. Levasseur, 2001 Molecular investigations in populations of Spartina anglica C. E. Hubbard (Poaceae) invading coastal Brittany (France) Mol. Ecol 10:1689-1701[ISI][Medline]

    Bennetzen J. L., 2000 Transposable element contributions to plant gene and genome evolution Plant Mol. Biol 42:251-269[ISI][Medline]

    Bevan M., I. Bancroft, E. Bent, et al. (65 co-authors) 1998 Analysis of 1.9 Mb of contiguous sequence from chromosome 4 of Arabidopsis thaliana Nature 391:485-488[ISI][Medline]

    Boyko E., R. Kalendar, V. Korzun, B. Gill, A. H. Schulman, 2002 Combined mapping of Aegilops tauschii by retrotransposon, microsatellite, and gene markers Plant Mol. Biol. (in press)

    Comai L., A. P. Tyagi, K. Winter, R. Holmes-Davis, S. H. Reynolds, Y. Stevens, B. Byers, 2000 Phenotypic instability and rapid gene silencing in newly formed Arabidopsis allotetraploids Plant cell 12:1551-1567[Abstract/Free Full Text]

    Cronn R., R. L. Small, J. F. Wendel, 1999 Duplicated genes evolve independently following polyploid formation in cotton Proc. Natl. Acad. Sci. USA 96:14406-14411[Abstract/Free Full Text]

    Feldman M., B. Liu, G. Segal, S. Abbo, A. A. Levy, J. M. Vega, 1997 Rapid elimination of low-copy DNA sequences in polyploid wheat: a possible mechanism for differentiation of homoeologous chromosomes Genetics 147:1381-1387[Abstract/Free Full Text]

    Ferris C., R. A. King, A. J. Gray, 1997 Molecular evidence for the maternal parentage in the hybrid origin of Spartina anglica C. E. Hubbard Mol. Ecol 6:185-187[ISI]

    Flavell A. J., M. R. Knox, S. R. Pearce, T. H. N. Ellis, 1998 Retrotransposon-based insertion polymorphisms (RBIP) for high throughput marker analysis Plant J 16:643-650[ISI][Medline]

    Grandbastien M. A., A. Spielmann, M. Caboche, 1989 Tnt1, a mobile retroviral-like transposable element of tobacco isolated by plant cell genetics Nature 337:376-380[ISI][Medline]

    Gray A. J., P. E. M. Benham, A. F. Raybould, 1990 Spartina anglica—the evolutionary and ecological background Pp. 5–10 in A. J. Gray and P. E. M. Bonham, eds. Spartina anglica—a research review. HMSO, London

    Guénégou M. C., J. Citharel, J. E. Levasseur, 1988 The hybrid status of Spartina anglica (Poaceae). Enzymatic analysis of the species and the presumed parents Can. J. Bot 66:1830-1833[ISI]

    Hirochika H., K. Sugimoto, Y. Otsuki, H. Tsugawa, M. Kanda, 1996 Retrotransposons of rice involved in mutations induced by tissue culture Proc. Natl. Acad. Sci. USA 93:7783-7788[Abstract/Free Full Text]

    Hsiao C., S. W. L. Jacobs, N. J. Chatterton, K. H. Asay, 1999 A molecular phylogeny of grass family (Poaceae) based on the sequences of nuclear ribosomal DNA (ITS) Aust. Syst. Bot 11:667-668[ISI]

    Hubbard C. E., 1968 Grasses. 2nd edition Penguin Books, London

    Kalendar R., T. Grob, M. T. Regina, A. Suoniemi, A. H. Schulman, 1999 IRAP and REMAP: two new retrotransposon-based DNA fingerprinting techniques Theor. Appl. Genet 98:704-711[ISI]

    Kalendar R., J. Tanskanen, S. Immonen, E. Nevo, A. H. Schulman, 2000 Genome evolution of wild barley (Hordeum spontaneum) by BARE-1 retrotransposon dynamics in response to sharp microclimatic divergence Proc. Natl. Acad. Sci. USA 97:6603-6607[Abstract/Free Full Text]

    Kellog E. A., 2001 Evolutionary history of grasses Plant Phys 125:1198-1205[Free Full Text]

    Kidwell M. G., D. R. Lisch, 2000 Transposable elements and host genome evolution Trends Evol. Ecol 15:95-99

    Kumar A., J. L. Bennetzen, 1999 Plant Retrotransposons Annu. Rev. Genet 33:479-532[ISI][Medline]

    Liu B., C. L. Brubaker, G. Mergeai, R. C. Cronn, J. F. Wendel, 2001 Polyploid formation in cotton is not accompanied by rapid genomic changes Genome 44:321-330[ISI][Medline]

    Liu B., J. M. Vega, G. Segal, S. Abbo, M. Rodova, M. Feldman, 1998a. Rapid genomic changes in newly synthesized amphiploids of Triticum and Aegilops. I. Changes in low-copy noncoding DNA sequences Genome 41:272-277[ISI]

    ———. 1998b. Rapid genomic changes in newly synthesized amphiploids of Triticum and Aegilops. II. Changes in low-copy coding DNA sequences Genome 41:535-542[ISI][Medline]

    Liu B., J. F. Wendel, 2000 Retroelement activation followed by rapid repression in interspecific hybrid plants Genome 43:874-880[ISI][Medline]

    Manninen O., R. Kalendar, J. Robinson, A. H. Schulman, 2000 Application of BARE-1 retrotransposon markers to map a major resistance gene for net blotch in barley Mol. Gen. Genet 264:325-334[ISI][Medline]

    Marchant C. J., 1967 Corrected chromosomes numbers for Spartina x townsendii and its parents species Nature 31:929.

    ———. 1968 Evolution in Spartina (Gramineae). II. Chromosomes, basic relationships and the problem of Spartina x townsendii agg Bot. J. Linn. Soc 60:381-409

    Matzke M. A., A. J. M. Matzke, 1998 Polyploidy and transposons Trends Ecol. Evol 13:241.

    Matzke M. A., O. Mittelsten Scheid, A. J. M. Matzke, 1999 Rapid structural and epigenetic changes in polyploid and aneuploid genomes Bioessay 21:761-767[ISI][Medline]

    McClintock B., 1984 The significance of responses of the genome to challenge Science 256:792-801

    O'Neill R. J. W., M. J. O'Neill, J. A. Marshall Graves, 1998 Undermethylation associated with retroelement activation and chromosome remodelling in an interspecific mammalian hybrid Nature 393:68-72[ISI][Medline]

    Ozkan H., A. A. Levy, M. Feldman, 2001 Allopolyploidy-induced rapid genome evolution in the wheat (Aegilops-Triticum) group Plant Cell 13:1735-1747[Abstract/Free Full Text]

    Ramsay L., M. Macaulay, L. Cardle, M. Morgante, S. Degli-Ivanissevich, E. Maestri, W. Powell, R. Waugh, 1999 Intimate association of microsatellite repeats with retrotransposons and other dispersed repetitive elements in Barley Plant J 17:415-423[ISI][Medline]

    Raybould A. F., A. J. Gray, M. J. Lawrence, D. F. Marshall, 1991a. The evolution of Spartina anglica C. E. Hubbard (Gramineae): origin and genetic variability. Bot J. Linn. Soc 43:111-126

    ———. 1991b. The evolution of Spartina anglica C. E. Hubbard (Gramineae): genetic variation and status of the parental species in Britain Bot. J. Linn. Soc 44:369-380

    San Miguel P., A. Tikhonov, Y. K. Jin, et al. (11 co-authors) 1996 Nested retrotransposons in the intergenic regions of the maize genome Science 274:765-768[Abstract/Free Full Text]

    Schulman A. H., A. J. Flavell, T. H. N. Ellis, 2002 The application of LTR retrotransposons as molecular markers in plants In W. Miller and J. P. Capy, eds. Mobile genetic elements and their application in genomics. Humana Press, Totawa, NJ. (in press)

    Shaked H., K. Kashkush, H. Ozkan, M. Feldman, A. A. Levy, 2001 Sequence elimination and cytosine methylation are rapid and reproducible responses of the genome to wide hybridisation and allopolyploidy in wheat Plant Cell 13:1749-1759[Abstract/Free Full Text]

    Shirasu K., A. H. Schulman, T. Lahaye, P. Schulze-Lefert, 2000 A contiguous 66 kb barley DNA sequence provides evidence for reversible genome expansion Genome Res 10:908-915[Abstract/Free Full Text]

    Small R. L., J. F. Wendel, 2000 Copy number lability and evolutionary dynamics of the Adh gene family in diploid and tetraploid cotton (Gossypium) Genetics 155:1913-1926[Abstract/Free Full Text]

    Song K., P. Lu, K. Tang, T. C. Osborn, 1995 Rapid genome change in synthetic polyploids of Brassica and its implications for polyploid evolution Proc. Natl. Acad. Sci. USA 92:7719-7723[Abstract]

    Suoniemi A., D. Schmidt, A. H. Schulman, 1997 BARE-1 insertion site preferences and evolutionary conservation of RNA and cDNA processing sites Genetica 100:219-230[ISI][Medline]

    Vicient C. M., M. Jääskeläinen, R. Kalendar, A. H. Schulman, 2001 Active retrotransposons are a common feature of grass genomes Plant Physiol 125:1283-1292[Abstract/Free Full Text]

    Vicient C. M., A. Suoniemi, K. Anamthawat-Jónsson, J. Tanskanen, A. Beharav, E. Nevo, A. H. Schulman, 1999 Retrotransposon BARE-1 and its role in genome evolution in the genus Hordeum Plant Cell 11:1769-1784[Abstract/Free Full Text]

    Wendel J. F., 2000 Genome evolution in polyploids Plant Mol. Biol 42:225-249[ISI][Medline]

    Wendel J. F., A. Schnabel, T. Seelanan, 1995 Bidirectional interlocus concerted evolution following allopolyploid speciation in cotton (Gossypium) Proc. Natl. Acad. Sci. USA 92:280-284[Abstract]

    Wendel J. F., S. R. Wessler, 2000 Retrotransposon-mediated genome evolution on a local ecological scale Proc. Natl. Acad. Sci. USA 97:6250-6252[Free Full Text]

    Wessler S. R., 1996 Plant retrotransposons: turned on by stress Curr. Biol 6:959-961[ISI][Medline]

    Wolfe A. D., A. Liston, 1998 Contributions of the polymerase chain reaction to plant systematics Pp. 43–86 in D. E. Soltis, P. S. Soltis, and J. J. Doyle, eds. Molecular systematics of plants II: DNA sequencing. Kluwer, New York

    Zhao X. P., Y. Si, R. E. Hanson, C. R. Crane, H. J. Price, D. M. Stelly, J. F. Wendel, A. H. Paterson, 1998 Dispersed repetitive DNA has colonized new genomes since polyploid formation in cotton Genome Res 8:479-492[Abstract/Free Full Text]

Accepted for publication February 7, 2002.