The Evolution of Ty1-copia Group Retrotransposons in Gymnosperms

Caroline Stuart-Rogers1, and Andrew J. Flavell

Department of Biochemistry, University of Dundee, Scotland


    Abstract
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 literature cited
 
A diverse collection of Ty1-copia group retrotransposons has been characterized from the genome of Picea abies (Norway spruce) by degenerate PCR amplification of a region of the reverse transcriptase gene. The occurrence of these retrotransposable elements in the gymnosperms was investigated by Southern blot hybridization analysis. The distribution of the different retrotransposons across the gymnosperms varies greatly. All of the retrotransposon clones isolated are highly conserved within the Picea (spruce) genus, many are also present in Pinus (pine) and/or Abies (fir) genera, and some share strongly homologous sequences with one or more of cedar, larch, Sequoia, cypress, and Ginkgo. Further subclones of one of the most strongly conserved retrotransposon sequences, Tpa28, were obtained from Ginkgo and P. abies. Comparisons of individual sequence pairs between the two species show nucleotide cross-homologies of around 80%–85%, corresponding to nucleotide substitution rates similar to those of nuclear protein-coding genes. Analysis of Tpa28 consensus sequences reveals that strong purifying selection has acted on this retrotransposon in the lineages connecting Ginkgo and Picea. Collectively, these data suggest, first, that the evolution of the Ty1-copia retrotransposon group in the gymnosperms is dominated by germ line vertical transmission, with strong selection for reverse transcriptase sequence, and, second, that extinction of individual retrotransposon types has been comparatively rare in gymnosperm species lineages compared with angiosperms. If this very high level of sequence conservation is a general property of the retrotransposons, then their extreme sequence diversity implies that they are extremely ancient, and the major element lineages seen today may have arisen early in eukaryote evolution. The data are also consistent with horizontal transmission of particular retrotransposons between species, but such a mechanism is unnecessary to explain the results.


    Introduction
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 literature cited
 
Retrotransposons are the dominant class of dispersed repetitious DNA in many eukaryote nuclear genomes. These mobile genetic elements transpose replicatively through RNA intermediates, and this can result in their accumulation to very high copy numbers (Pearce et al. 1996Citation ; SanMiguel et al. 1996Citation ; Kumar and Bennetzen 1999Citation ). For at least one eukaryote, maize, retrotransposons have expanded to occupy around one half of the entire genome in a few million years, causing large-scale changes in genome size (SanMiguel et al. 1996, 1998Citation ).

Ty1-copia group retrotransposons are present as large, highly heterogeneous populations within all higher plant divisions (Flavell et al. 1992Citation ; Flavell, Smith, and Kumar 1992Citation ; Voytas et al. 1992Citation ). Many phylogenetic comparisons among these populations of mobile elements have been carried out, both within and between species. These studies have shown the existence of homologous retrotransposon pairs in divergent species and discontinuous distributions of many retrotransposons among species groups (Flavell et al. 1992Citation ; Flavell, Smith, and Kumar 1992Citation ; Voytas et al. 1992Citation ). This led to suggestions that retrotransposons may have been transmitted horizontally between plant species. However, these early studies were inconclusive, because the relatively modest levels of sequence similarity detected are also consistent with the alternative hypothesis of vertical transmission, combined with low germ line transposition rates and strong selection for retrotransposon sequence conservation (Vanderwiel, Voytas, and Wendel 1993Citation ; Capy, Anxolabéhere, and Langin 1994Citation ; Flavell et al. 1997Citation ). Similarly, the discontinuities in element distributions do not prove horizontal transfer and can be explained by extinctions in multiple species lineages (Capy, Anxolabéhere, and Langin 1994Citation ). Subsequent phylogenetic studies have failed to provide convincing evidence for horizontal transfer of plant retrotransposons (Vanderwiel, Voytas, and Wendel 1993Citation ; Gribbon et al. 1999Citation ; Matsuoka and Tsunewaki 1999Citation ), and it is now believed to be either rare or nonexistent. The prevailing view of plant retrotransposon evolution is that these genetic elements already existed early in plant evolution and diverged into heterogeneous subgroups before modern plant orders arose, and the descendants of these subgroups have been transmitted vertically to the descendant species.

Recently, however, strong evidence for horizontal transmission of animal retrotransposons has been reported (Gonzales and Lessios 1999Citation ; Jordan, Matyunina, and McDonald 1999Citation ), and this suggests a reappraisal of the above model for the evolution of plant retrotransposons. In this study, the heterogeneous population of Ty1-copia group retrotransposons in Norway spruce (Picea abies) was characterized and the distribution of the different retrotransposons within the branch Gymnospermae was investigated. The picture that emerged was of a very diverse family of mobile elements within P. abies which shows a complex distribution pattern within the gymnosperms, varying from genus-specific retrotransposons to extremely widely distributed and very well conserved elements. These data support the hypothesis that two counteracting processes have dominated the evolution of this transposon group in the gymnosperms. The first of these is vertical transmission down lineages, with strong selection for the sequence of at least the reverse transcriptase, and the second is extinction of retrotransposons in particular lineages. The data do not exclude the possibility of horizontal transmission of retrotransposons between species, but they do not require it.


    Materials and Methods
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 literature cited
 
Isolation and Characterization of Ty1-copia Group Retrotransposons
Ty1-copia group retrotransposon reverse transcriptase gene fragments were amplified by degenerate PCR of P. abies DNA as previously described (Flavell et al. 1992Citation ). The PCR products were subcloned into HincII-digested M13mp19 vector and sequenced using Sequenase. Tpa9Tpa19 correspond to EMBL database accession numbers AJ288018–AJ288028, respectively, Tpa21Tpa44 to numbers AJ288030—AJ288053, respectively, and Tpa45Tpa55 to numbers AJ288171—AJ288181, respectively. Sequences Tgb1Tgb19 correspond to EMBL database accession numbers AJ288182—AJ288200, respectively.

Sequence Comparisons and Phylogenetic Trees
Computing was carried out at the SEQNET facility of the Daresbury Laboratory and the HGMP facility at Hinxton, CB10 1SB, United Kingdom (http://www.hgmp.mrc.ac.uk/). Nucleotide sequences were conceptually translated in three reading frames, which were compared visually with a reference set of plant retrotransposon peptide sequences to identify frameshifts. The corrected peptide sequences were aligned, and a phylogenetic tree was generated using CLUSTAL W (Thompson, Higgins, and Gibson 1994Citation ). Ratios of synonymous to nonsynonymous nucleotide substitutions were calculated using DIVERGE, version 10.0 (Wisconsin Package, Genetics Computer Group [GCG], Madison, Wis.). Nucleotide substitution rates were calculated by the method of Jukes and Cantor (1969)Citation . Consensus nucleotide sequences were obtained using CONSENSUS (part of the GCG package).

Southern Analysis
DNAs were isolated from fresh needles or leaves using Nucleon Phytopure kits. In some cases, DNA samples which could not be restriction-digested after purification were banded by CsCl density gradient ultracentrifugation. Southern blotting to Biodyne A nylon membranes (Pall) used the manufacturer's conditions. Blot hybridization to radioactively labeled retrotransposon probes prepared by the random oligo-primed labeling method was at Tm - 20°C in 0.75 M Na+. Final washing of the blots was in 2 x SSC (0.3 M Na+) at 64°C (fig. 3 ) or in 2 x SSC at 64°C or 0.1 x SSC (0.015 M Na+) at 65°C (fig. 4 ). Blots were stripped by the manufacturer's recommended method and re-exposed to check stripping efficiency before rehybridization.



View larger version (67K):
[in this window]
[in a new window]
 
Fig. 3.—Conservation of Tpa28 in the Gymnospermae at high wash stringency. a, A Southern blot probed with Tpa28 was washed at moderate stringency (2 x SSC, 64°C) and exposed to autoradiography for 1 day. b, The same blot was then washed at high stringency (0.1 x SSC, 65°C) and exposed for 3 days

 


View larger version (26K):
[in this window]
[in a new window]
 
Fig. 4.—Phylogenetic tree of Tpa28-homologous retrotransposons from Picea abies and Ginkgo biloba. Tpa retrotransposons are derived from P. abies and Tgb retrotransposons are from G. biloba. Divergences in distance units between sequences (Thompson, Higgins, and Gibson 1994Citation ) are indicated by horizontal branch lengths. Bootstrap values of 50% and higher are shown (1,000 trials). The tree is unrooted, so the relative lengths of the two horizontal sections of the deepest branch have no significance, but the total of the two lengths is equal to the distance separating Tpa45/Tpa28 from the other sequences. The groups of sequences (A–E) from which consensuses were derived for nucleotide substitution analysis (table 1 ) are bracketed

 
PCR Analysis of Tpa28 Reverse Transcriptase Fragments from Picea and Ginkgo
DNA sequence comparisons of Tpa28 and its closest homologs identified in P. abies identified a primer just upstream of the conserved YVDDM motif of the RT gene (CCACAATGAAGTATAGGTTGG) which is highly specific for Tpa28 alone. PCR with this primer and the conserved upstream primer (encoding TAFLHG; Flavell, Smith, and Kumar 1992Citation ) generated bands of the expected size (~220 bp) from both P. abies and Ginkgo biloba. These bands were subcloned and sequenced as above. The possibility of cross-contamination artifacts between the Picea and the Ginkgo DNA preparations was excluded by PCR using primer pairs specific for P. abies and G. biloba phytochrome genes. PCR primers selective for the P. abies phytochrome gene (EMBL accession number X80298) were ACTCCCGTGCAGGTTATCCAGTC and AGCATAGCGGAGAGGAAACGGA (band size = 254 bp), and the corresponding primers for G. biloba (U08162) were TTGCTGTGCAACTCCTGTT and CCTTCTTCATCAGTTCCATTG (band size = 172 bp). These primer pairs produced strong bands of the predicted sizes from their respective parent species, but no bands were detected in cross-species PCRs (data not shown).


    Results
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 literature cited
 
Characterization of Ty1-copia Group Retrotransposons of P. abies and Their Occurrence in Other Gymnosperms
Ty1-copia group retrotransposon sequences were amplified from P. abies by degenerate PCR of a part of the reverse transcriptase gene using primers which have been successfully used for a variety of higher plant species (Flavell et al. 1992Citation ; Flavell, Smith, and Kumar 1992Citation ; Pearce et al. 1996Citation ; Gribbon et al. 1999Citation ). The PCR products were subcloned, inserts were sequenced, and the inferred amino acid sequences were aligned. A phylogenetic tree for these sequences is shown in figure 1a. A highly heterogeneous collection of Ty1-copia group retrotransposons was evident, as is typical for plants with large genomes (Flavell et al. 1992Citation ; Flavell, Smith, and Kumar 1992Citation ; Vanderwiel, Voytas, and Wendel 1993Citation ; Pearce et al. 1996Citation ; Gribbon et al. 1999Citation ).



View larger version (53K):
[in this window]
[in a new window]
 
Fig. 1.—a, Phylogenetic tree of Picea abies retrotransposons. Divergences in distance units (Thompson, Higgins, and Gibson 1994Citation ) between retrotransposon sequences are indicated by horizontal branch lengths. Percentage bootstrap values are shown (1,000 trials). The tree is unrooted, so the relative lengths of the two horizontal sections of the deepest branch have no significance, but the total of the two lengths is equal to the distance separating Tpa42 from the other sequences. b, Summary of Southern data showing P. abies retrotransposon distributions in gymnosperms. Clusters of probes generating similar Southern hybridization patterns are bracketed. The relative hybridization signals are from subjective assessment and are not quantitative (see fig. 2 for actual data for six representative probes)

 
To determine the distribution of these retrotransposable elements across the Gymnospermae, each of the 35 retrotransposon clones was hybridized to Southern blots (figs. 1b and 2 ) of DNAs from 18 species which represent the diversity of the gymnosperm branch but are skewed in favor of Picea (six species) and its sister genus Pinus (five species). For both of these genera, species from the Old and New Worlds were represented. The other seven species comprise one representative each from the Abies (fir), Cedrus (cedar), Larix (larch), Ginkgo (maidenhair tree), Sequoiadendron (redwood), Taxodium (cypress), and Welwitschia (welwitschia) genera.

The different retrotransposon probes produced a wide variety of Southern patterns. Some of these are shown in figure 2 , and the overall results are summarized in figure 1b. The distribution of the different retrotransposons across the gymnosperms varies greatly. For example, Tpa43 is strictly confined to the genus Picea (fig. 2a ), where it is present in high copy numbers (300 copies per genome would appear as a faint band under the conditions used here), suggesting that most of these sequences are present in very high copy numbers in these large genomes (ca. 2 x 1010 bp). Tpa12 shows additional very weak hybridization in Pinus (fig. 2b ), despite the much lower signal in Picea. Tpa17 and Tpa35 are also strongly represented across all tested Picea species (fig. 2c and d ), where the Southern patterns are strongly conserved, and faint signals are discernible in Pinus. The within-genus conservation of Southern pattern is shared by all of the Ty1-copia group retrotransposons studied here (fig. 2 and data not shown), in agreement with previous results for the Tpe1 element in Pinus (Kamm et al. 1996Citation ).



View larger version (51K):
[in this window]
[in a new window]
 
Fig. 2.—Distribution of particular Picea retrotransposon homologs in the Gymnospermae shown by Southern analysis at moderate wash stringency (2 x SSC, 64°C). Lane 1, marker; lanes 2–7, Picea samples; lane 2, P. abies; lane 3, P. glauca; lane 4, P. obovata (partially degraded input DNA sample); lane 5, P. orientalis; lane 6, P. rubens; lane 7, P. sitchensis; lanes 8–12, Pinus samples; lane 8, P. cembra; lane 9, P. coulterii; lane 10, P. elliottii; lane 11, P. sylvestris var. Scotica; lane 12, P. thunbergii; lane 13, Abies alba (failed digestion in i); lane 14, Cedrus libani breviliensis; lane 15, Larix europaea; lane 16, Ginkgo biloba; lane 17, Sequoiadendron gigantum; lane 18, Taxodium distichum; lane 19, Welwitschia mirabilis; lane 20, marker

 
Other retrotransposons show higher conservation across larger evolutionary gaps. For example, Picea Tpa29 homologs are detected strongly in the genera Picea and Abies and weakly in Ginkgo (fig. 2e ). Tpa27 homologs also share this species distribution, with the addition that strongly homologous sequences are also seen in Pinus species and weaker signals are seen in Cedrus and Larix (fig. 2f ). Tpa22 shows similarities in species distribution to Tpa29, except for the absence of a signal in Abies alba and Ginkgo (fig. 2g ). The most extreme examples of strong conservation of Ty1-copia group retrotransposons across the range of gymnosperms tested are found for Tpa28 and Tpa13 (fig. 2h and i ). Tpa28-homologous sequences are present in every gymnosperm tested except Taxodium and Welwitschia, and Tpa13-homologous sequences are found in all of the species tested. There are many cases in which particular retrotransposon restriction fragments are conserved across large evolutionary gaps (fig. 2 and data not shown), and this is particularly apparent for Tpa28 and Tpa13 (fig. 2h and i ), showing that two six-base sequences and the size gap between them are conserved between plants which are separated by at least 250 Myr (see below).

The Southern blots in figure 2 had been washed at moderate stringency (2 x SSC, 64°C). To get a better idea of the level of homology shared between one of the best conserved retrotransposons, Tpa28, in the gymnosperms, high-stringency Southern wash conditions (0.1 x SSC, 65°C) were used (fig. 3b ). These conditions are 2°C below the calculated Tm for a perfect Tpa28 hybrid. This treatment reduced the Picea hybridization signals roughly 10-fold relative to a moderate-stringency wash (fig. 3a ; the exposure time for fig. 3b is threefold longer than that used for fig. 3a ). The strongest Ginkgo-specific bands were still just visible after this treatment, showing that close homologs of Tpa28 were present in the Ginkgo DNA sample.

Characterization and Comparison of Ty1-copia Group Retrotransposons of P. abies and G. biloba
The phylogenetic relationship between Picea and Ginkgo homologs of one of the best conserved retrotransposons, Tpa28, was studied further. PCR primers, which specifically amplify Tpa28 alone out of all 35 P. abies elements cloned, were designed (see Materials and Methods). PCR with either of these primers and the degenerate upstream primer used in the original cloning of all of these elements (Flavell et al. 1992Citation ) was carried out to isolate Tpa28 homologs from both P. abies and G. biloba. Multiple subclones of these were sequenced as before, and the phylogenetic relationships among these sequences are shown in figure 4 . The Ginkgo Tpa28 homologs form a well-defined clade with a bootstrap value of 97 which is related to a clade of three P. abies retrotransposons, Tpa52, Tpa53, and Tpa55.

The strong similarity among Tpa28 homologs from widely diverged species and that which is evident from the Southern data in figure 2 for other retrotransposons can be explained by two contrasting models. The first proposes that retrotransposons have been transmitted vertically down the species lineages, with very strong selection acting on at least the reverse transcriptase gene. The second model proposes that some retrotransposons have been transmitted horizontally between gymnosperm species in the relatively recent evolutionary past. Vertical transmission is the normal way in which eukaryotes acquire their retroelements, but extreme sequence conservation across large evolutionary distances has shown that horizontal transmission of retrotransposons can occur (Gonzales and Lessios 1999Citation ; Jordan, Matyunina, and McDonald 1999Citation ).

The vertical-transmission model rests on the hypothesis that strong selective pressure has acted on the retrotransposon sequence in the species lineages. We therefore looked more closely at the nucleotide substitution rates of the Tpa28 homologs in figure 4 , both within and between P. abies and Ginkgo, to see whether this hypothesis was valid. Tpa55 and Tgb12 share 85.6% nucleotide identity (fig. 5 ), corresponding to an overall nucleotide substitution rate for this region of 6.1 x 10-10 per site per year (assuming a divergence time of 260 MYA for these two species; see Discussion). The ratio of synonymous substitutions to nonsynonymous substitutions (Ks/Ka) between Tgb12 andTpa55 is 6.4, implying quite strong selection for peptide sequence in this region. For comparison, we determined the conservation between a region of the phytochrome gene over this time gap. The two phytochrome gene homologs X80298 and gb08162, of Pinus and Ginkgo, respectively, share 80% nucleotide identity and a Ks/Ka ratio of 10.3. Therefore, these species lineages appear to have been under lower selective pressure for Tpa28 reverse transcriptase sequence than they have for phytochrome, but their sequence has been better preserved; i.e., this selective pressure might appear to be insufficient to have conserved these two sequences to the degree evident here. However, retrotransposons exhibit a "pseudogene effect" (McAllister and Werren 1997Citation ), whereby newly transposed copies are unlikely to represent sources of new elements and therefore accumulate mutations at the same rate as a pseudogene. A pseudogene effect might distort the Ks/Ka ratio for Tgb12 and Tpa55. Consistent with this hypothesis, both Tpa55 and Tgb12 contain stop codons, and the high frequency of stops and frameshifts among the Tpa28 homologs sequenced in this study (data not shown) suggest that the majority are pseudogenes.



View larger version (43K):
[in this window]
[in a new window]
 
Fig. 5.—Alignment of Ginkgo biloba (Tgb12) and Picea abies (Tpa55) homologs of Tpa28.

 
To correct for the pseudogene effect, it is necessary to identify the particular Tpa28 retrotransposons which are being selected for transposition, i.e., those copies which are transposing replicatively. This is a very difficult task for a high-copy-number retrotransposon, but the sequence(s) of transposing elements can be inferred from the consensus of the genomic copies, which derive from the transposing elements. Therefore, consensuses were generated for the five groups of sequences shown in figure 4 , and substitution rates were calculated between these consensus sequences (table 1 ). The Ks/Ka ratio between the consensus sequence for the clade containing Tgb12 (clade D) and the consensus of all Picea Tpa28 homologs (group A) is 15.9, and the average Ks/Ka ratio between clade D and the other three groupings (A–C) is 12.8 (SD = 2.5), showing that the pseudogene effect has indeed partially masked the strong conservation of these two retrotransposon lineages and there has been strong selection acting on the reverse transcriptase of transposing homologs of Tpa28. This result is consistent with the vertical-transmission model, and, although horizontal transmission could also explain the high level of sequence conservation between the Tpa28 homologs of Picea and Ginkgo, such a mechanism is unnecessary and less parsimonious.


View this table:
[in this window]
[in a new window]
 
Table 1 Nucleotide Divergences Between Consensus Sequences of Tpa28 Homologs in Picea abies and Ginkgo biloba

 
Strong conservation for the encoded protein sequence of clade C, the other Ginkgo Tpa28 variant, with the Picea group A and clade D consensuses was also seen (Ks/Ka values of 10.2 and 12.9, respectively). Therefore, both Ginkgo Tpa28 variants have been strongly conserved against both each other and their homologs in Picea. Paradoxically, there was little pseudogene effect apparent for clade B, the group of three sequences containing Tpa55 itself and its two closest homologs. Comparisons between clade B and the Ginkgo clades gave an average Ks/Ka ratio of 7.8. This may be an artifact; Tpa52 and Tpa53 differ by only 2 nt; therefore, the consensus between these two sequences and Tpa55 is almost identical to Tpa52 or Tpa53.


    Discussion
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 literature cited
 
The results presented here show that in some ways the Ty1-copia group retrotransposons of gymnosperms are typical of other plants with large genome sizes. They are a complicated heterogeneous population of interrelated sequences, most of which are present in very high overall copy numbers, and all individual retrotransposons transcend species boundaries. However, the degree of conservation of some of these retrotransposons is the highest reported to date for plants (Flavell et al. 1992, 1997Citation ; Flavell, Smith, and Kumar 1992Citation ; Voytas et al. 1992Citation ; Vanderwiel, Voytas, and Wendel 1993Citation ; Matsuoka and Tsunewaki 1997Citation ; Kumar and Bennetzen 1999Citation ). In plants, it is common to see Southern blot cross-hybridization for retrotransposons between related genera (Matsuoka and Tsunewaki 1997Citation ; Pearce et al. 1997Citation ; Chavanne et al. 1998Citation ), and such signals occasionally persist between angiosperm tribes such as the Triticeae and the Poaceae (Matsuoka and Tsunewaki 1997Citation ), which are separated by ca. 60 Myr, but the data presented here comprise, to our knowledge, the first reported case in which strong cross-hybridization between retroelements extends to the taxonomic level of order, across hundreds of millions of years.

Two striking examples of high sequence conservation presented here are the Tpa28 and Tpa13 elements (fig. 2h and i ). These conifer-derived retrotransposons detect strongly homologous sequences across evolutionary gaps up to and including Ginkgo, the only remaining representative of the order Ginkgoales, which is first recorded definitively in the lower Permian, approximately 260 MYA (Stewart and Rothwell 1993Citation ). Modern conifers arose in the Triassic (>200 MYA). Virtually all phylogenetic analyses place Ginkgo basal to the conifer lineage (Crane 1985Citation ; Nixon et al. 1994Citation ; Qiu et al. 1999Citation ; Soltis, Soltis, and Chase 1999Citation ). There is therefore a minimum separation of 260 Myr between the modern representatives of the two orders. Further investigation of the Ginkgo and Picea homologs of Tpa28 shows a pair of nucleotide sequences from these species which are 85.6% identical, corresponding to an overall nucleotide substitution rate of 6.1 x 10-10 per site per year. The comparable rate for the phytochrome gene of these plants is 9.9 x 10-10 per site per year (EMBL accession numbers GB08162 and X80298). A more comprehensive comparison between four genes (from Pinus and Ginkgo; Qiu et al. 1999Citation ) shows substitution rates between 2.0 x 10-10 and 3.6 x 10-10 per site per year (data not shown). The conservation of the RT gene fragment described here is thus broadly similar to that seen for nuclear protein-coding genes.

For animals, nucleotide substitution rates similar to those seen here have been interpreted to be possibly due to horizontal transfer (ca. 3 x 10-10 to 8 x 10-10 per site per year for gypsy group retrotransposons of echinoids; Gonzales and Lessios 1999Citation ), and 100% conservation of retrotransposons in Drosophila and echinoid species has provided conclusive evidence for horizontal transmission of these retrotransposons between insect species (Gonzales and Lessios 1999Citation ; Jordan, Matyunina, and McDonald 1999Citation ). It is tempting to suggest that horizontal transmission has been responsible for the closely related Ty1-copia group retrotransposons in the gymnosperms. By such a model, Tpa28 and the other strongly conserved elements were introduced horizontally into one or both of the Ginkgo or Picea lineages long after these two lineages split but still millions of years ago, and in the ensuing period, they have become amplified by retrotransposition, with subsequent mutation yielding the many defective copies seen today.

To accept the horizontal-transmission model, it is necessary to exclude beyond reasonable doubt the alternative model that vertical transmission has been responsible for the data presented here. We are not yet convinced that this is the case, for the following reasons. First, studies on clades of echinoid retrotransposons, for which the retrotransposon phylogeny resembles the phylogenies of the corresponding organism lineages, show that the retrotransposon sequences can be conserved at rates between 1.3 x 10-9 per site per year and 2.6 x 10-8 per site per year (Gonzales and Lessios 1999Citation ). The lowest of these rates is only double that seen for Tpa28. Second, detailed study of the ratios of synonymous substitutions to nonsynonymous substitutions (Ks/Ka) between the Ginkgo and Picea Tpa28 homologs in figure 5 shows that they have been under strong purifying selection (see Results and table 1 ). The third reason for favoring the vertical-transmission model at present is the distribution of most of the retrotransposons across the gymnosperms (fig. 1b ). For the majority of the Ty1-copia group retrotransposons studied here, the hybridization signal is very roughly proportional to the evolutionary distance of a given species from Picea, and there are few examples of major signal discontinuities across the species tested here. Such a result is entirely consistent with vertical transmission, but horizontal transmission would be expected to lead to departures from this situation. Furthermore, most of the retrotransposons in this study are conserved between conifers and Ginkgo (fig 1b ), and if horizontal transmission has been responsible for the wide distribution of one of these (Tpa28), then presumably it has been involved in the wide distributions of the others; i.e., it has dominated the evolution of the majority of gymnosperm Ty1-copia group retrotransposons. This seems unlikely, because vertical transmission is far better established in all other eukaryote branches. In conclusion, our data do not rule out the possibility of horizontal transmission of Tpa28 or other retrotransposons between gymnosperm species, but they do not require it either, and the most parsimonious explanation is that vertical transmission, with strong sequence conservation, has dominated the evolution of this retrotransposon group in the gymnosperms.

There is no obvious reason why such strong sequence conservation should not act on the Ty1-copia group retrotransposons of all eukaryotes. The basic structures and life cycles of long-terminal-repeat (LTR) retrotransposons are shared between fungi, plants, and animals, and, as mentioned above, previous studies of animal LTR retrotransposons have revealed strong conservation for retrotransposon reverse transcriptase sequence (McAllister and Werren 1997Citation ; Gonzales and Lessios 1999Citation ). However, the large majority of the plant retrotransposons studied previously show much lower conservation across large evolutionary divides (Flavell et al. 1992Citation ; Flavell, Smith, and Kumar 1992Citation ; Voytas et al. 1992Citation ; Matsuoka and Tsunewaki 1997, 1999Citation ; Gribbon et al. 1999Citation ). We propose that this apparent contradiction derives from the fact that the large majority of the plant sequences studied previously are from angiosperms.

Why should retrotransposons show lower rates of sequence conservation in angiosperms than they do in gymnosperms? There does not seem to be any significant difference between the overall rates of nucleotide change of nuclear genes in angiosperms and gymnosperms (Qiu et al. 1999Citation ; unpublished data). One possible explanation is that true angiosperm retrotransposon homologs, descended from a unique element in the nearest common ancestor, are rarely found, and paralogous sequences have usually been compared instead. A possible gymnosperm example of this is Tpa29 (fig. 2e ). The ancestral element of Tpa29 of Picea, which shows weak homology in Ginkgo, may have arisen long before the Ginkgoales-Coniferales split, and at the time of the split, it had diverged considerably into two variants (let us call them variants 1 and 2). If variant 1 would eventually become Tpa29 in P. abies and variant 2 would become its Ginkgo homolog, then the sequence divergence between them would date back to the last common ancestor, long before the Ginkgoales-Coniferales split. Paralogy has previously been invoked (Capy, Anxolabéhere, and Langin 1994Citation ) to explain phylogenetic inconsistencies in transposon occurrence between species.

In the above example, if variant 1 survived in both lineages we would see close homology between the Ginkgo and the Picea Tpa29 elements; i.e., variant 1 must have become extinct in the Ginkgo lineage. In general, if true retrotransposon homologs between widely diverged species are much rarer in angiosperms than in gymnosperms, and this is not due to a difference in substitution rates, then extinction of retrotransposon lineages must have been much more common in angiosperms than in gymnosperms. This is an attractive model, because a potential mechanism for extinction is very active in angiosperms but is less so in gymnosperms. Drastic differences in genome sizes within a genus are common in angiosperms (http://www.rbgkew.org.uk/cval/database1.html), and such changes can occur over a few million years (SanMiguel et al. 1998Citation ), but the genomes of gymnosperms appear to be more stable. For example, the Pinus genus is approximately 200 Myr old, all species are diploid, containing the same chromosome number, and the genome sizes are all huge (between 21 and 30 pg per haploid genome; Wakamiya et al. 1993Citation ). A chronically large genome can tolerate multiple different retrotransposon types, but a small genome size places strong constraints on retrotransposon copy number (Charlesworth, Sniegowski, and Stephan 1994Citation ), and competition between different low-copy-number elements should lead to frequent extinctions. If most of the angiosperm lineages extant today had ancestors with small genomes, retrotransposon extinction in those ancestors may have eliminated many of the different ancestral retrotransposons which appear to have survived so well in the gymnosperms.

Our findings also suggest a reevaluation of the well-known extreme heterogeneity seen among retrotransposons (Flavell et al. 1992Citation ; Flavell, Smith, and Kumar 1992Citation ; Voytas et al. 1992Citation ). For example, comparison of the reverse transcriptase regions (studied here) between the recently transposed (and presumably active) white-apricot copia retrotransposon insertion of Drosophila and clade D consensus gives a Ka value (rate of nonsynonymous substitutions per nonsynonymous site) of 0.51, compared with 0.050 for clade D versus all Picea homologs. If these two retrotransposon lineages have been subject to levels of selection similar to that of Tpa28 in the gymnosperms, then the divergence time between copia and clade D is on the order of 2.5 billion years, which is around the beginning of eukaryote life. Similar comparisons between clade D and a consensus of eight BARE-1 retrotransposon sequences of barley (Gribbon et al. 1999Citation ) or the yeast Saccharomyces cerevisiae Ty1 element give Ka values of around 0.60, which is close to saturation. An alignment illustrating the relative sequence homologies between these retrotransposons is shown in figure 6 .



View larger version (53K):
[in this window]
[in a new window]
 
Fig. 6.—Alignment of inferred peptide sequences of Tpa28 homologs (group A, clade D) with other Ty1-copia group retrotransposons from angiosperm plants (BARE-1 consensus), metazoa (copia), and fungi (Ty1). Amino acid identities are indicated by black boxes, and conservative substitutions are indicated by gray boxes

 
We propose that the ancestors of clade D, copia, BARE-1, Ty1, and the many elements with similar diversity levels (Flavell et al. 1992Citation ; Flavell, Smith, and Kumar 1992Citation ; Voytas et al. 1992Citation ) split in the early phase of eukaryote life. The Ty1-copia group is only one of the three retrotransposon groups, and the oldest of these was most likely the non-LTR retrotransposons (Xiong and Eickbush 1991Citation ). Malik, Burke, and Eickbush (1999)Citation deduced from phylogenetic data that the non-LTR retrotransposons (the probable ancestral retroransposon group) arose and diversified early in eukaryote evolution. If our proposition for the Ty1-copia group is confirmed, it may well turn out to be the case that the antiquity of non-LTR retrotransposons has been underestimated. Finally, all of these propositions depend on the assumption that retrotransposons have evolved at roughly similar rates in multiple diverse lineages. This is a very big assumption, but it should be testable.


    Acknowledgements
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 literature cited
 
C.S.-R. was funded by a BBSRC earmark studentship. This work was supported by grant 960508 (TEBIODIV) from the European Commission under the Framework IV Biotechnology Program. We thank Steve Pearce for much helpful advice with the work; Alan Schulman for constructive criticism of the manuscript; Steve Pearce, Amar Kumar, Noel Ellis, Pat Heslop-Harrison, and Alan Schulman for many interesting discussions; and Michele Morgante for P. abies DNA samples.


    Footnotes
 
Pierre Capy, Reviewing Editor

1 Present address: University of Helsinki, Institute of Biotechnology, Helsinki, Finland. Back

1 Keywords: retrotransposon gymnosperm Picea Ginkgo evolution Back

2 Address for correspondence and reprints: Andrew J. Flavell, Department of Biochemistry, University of Dundee, DD1 4HN Scotland, United Kingdom. E-mail: a.j.flavell{at}dundee.ac.uk Back


    literature cited
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 literature cited
 

    Capy, P., D. Anxolabéhere, and T. Langin. 1994. The strange phylogenies of transposable elements: are horizontal transfers the only explanation? Trends Genetics 10:7–12.

    Charlesworth B., P. Sniegowski, and W. Stephan. 1994. The evolutionary dynamics of repetitive DNA in eukaryotes. Nature 371:215–220.

    Chavanne, F., D.-X. Zhang, M.-F. Liaud, and R. Cerff. 1998. Structure and evolution of Cyclops: a novel giant retrotransposon of the Ty3/gypsy family highly amplified in pea and other legume species. Plant Mol. Biol. 37:363–375.[ISI][Medline]

    Crane, P. R. 1985. Phylogenetic relationships in seed plants. Cladistics 1:329–348.

    Flavell, A. J., E. Dunbar, R. Anderson, S. R. Pearce, R. Hartley, and A. Kumar. 1992. Ty1-copia group retrotransposons are ubiquitous and heterogeneous in higher plants. Nucleic Acids Res. 20:3639–3644.[Abstract]

    Flavell A. J., S. R. Pearce, J. S. Heslop-Harrison, and A. Kumar. 1997. The evolution of Ty1-copia group retrotransposons in eukaryote genomes. Genetica 100:185–195.

    Flavell, A. J., D. B. Smith, and A. Kumar. 1992. Extreme heterogeneity of TyTy1-copia group retrotransposons in plants. Mol. Gen. Genet. 231:233–242.[ISI][Medline]

    Gonzalez, P., and H. A. Lessios. 1999. Evolution of sea urchin retroviral-like (SURL) elements: evidence from 40 Echinoid species. Mol. Biol. Evol. 16:938–952.[Abstract]

    Gribbon, B. M., S. R. Pearce, R. Kalendar, A. H. Schulman, L. Paulin, P. Jack, and A. J. Flavell. 1999. Phylogeny and transpositional activity of Ty1-copia group retrotransposons in cereal genomes. Mol. Gen. Genet. 261:883–891.[Medline]

    Jordan, I. K., L. V. Matyunina, and J. F. McDonald. 1999. Evidence for the recent horizontal transfer of long terminal repeat retrotransposon. Proc. Natl. Acad. Sci. USA 96:12621–12625.

    Jukes, T. H., and C. R. Cantor. 1969. Evolution of protein molecules. Pp. 21–132 in H. N. Munro, ed. Mammalian protein metabolism. Academic Press, New York.

    Kamm, A., R. L. Doudrick, J. S. Heslop-Harrison, and T. Schmidt. 1996. The genomic and physical organisation of Ty1-copia-like sequences as a component of large genomes in Pinus elliottii and other gymnosperms. Proc. Natl. Acad. Sci. USA 93:2708–2713.

    Konieczny, A., D. F. Voytas, M. P. Cummings, and F. M. Ausubel. 1991. A superfamily of Arabidopsis thaliana retrotransposons. Genetics 127:801–809.

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

    McAllister, B. F., and J. H. Werren. 1997. Evolutionary analysis of a retrotransposon with implications for strong evolutionary constraints on reverse transcriptase. Mol. Biol. Evol. 14:69–80.[Abstract]

    Malik, H. S., W. D. Burke, and T. H. Eickbush. 1999. The age and evolution of non-LTR retrotransposable elements. Mol. Biol. Evol. 16:793–805.[Abstract]

    Matsuoka, Y., and K. Tsunewaki. 1997. Presence of wheat retrotransposons in Gramineae species and the origin of wheat retrotransposon families. Genes Genet. Syst. 72:335–343.[ISI][Medline]

    ———. 1999. Evolutionary dynamics of Ty1-copia group retrotransposons in grass shown by reverse transcriptase domain analysis. Mol. Biol. Evol. 16:208–217.[Abstract]

    Nixon, K. C., W. L. Crepet, D. Stevenson, and E. M. Friis. 1994. A reevaluation of seed plant phylogeny. Ann. Mo. Bot. Gard. 81:484–533.

    Pearce, S. R., G. Harrison, J. S. Heslop-Harrison, A. J. Flavell, and A. Kumar. 1997. Characterisation and genomic organisation of Ty1-copia group retrotransposons in rye (Secale cereale). Genome 40:617–625.

    Pearce, S. R., G. Harrison, D. Li, J. S. Heslop-Harrison, A. Kumar, and A. J. Flavell. 1996. The Ty1-copia group retrotransposons in Vicia species: copy number, sequence heterogeneity and chromosomal localisation. Mol. Gen. Genet. 250:305–315.[ISI][Medline]

    Qiu, Y. L., J. Lee, F. Bernasconi-Quadroni, D. E. Soltis, P. S. Soltis, M. Zanis, E. A. Zimmer, Z. Chen, V. Savolainen, and M. W. Chase. 1999. The earliest angiosperms: evidence from mitochondrial, plastid and nuclear genomes. Nature 402:404–407.

    SanMiguel, P., B. S. Gaut, A. Tikhonov, Y. Nakijama, and J. L. Bennetzen. 1998. The paleontology of intergene retrotransposons of maize. Nat. Genet. 20:43–45.[ISI][Medline]

    SanMiguel, 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.

    Soltis, P. S., D. E. Soltis, and M. W. Chase. 1999. Nature 402:402–404.

    Stewart, W. N., and G. W. Rothwell. 1993. Paleobotany and the evolution of plants. 2nd edition. Cambridge University Press, Cambridge, England.

    Thompson, J. D., D. G. Higgins, and T. J. Gibson. 1994. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22:4673–4680.[Abstract]

    Vanderwiel, P. L., D. Voytas, and J. F. Wendel. 1993. Copia-like retrotransposable element evolution in diploid and polyploid cotton (Gossypium L.). J. Mol. Evol. 36:429–447.[ISI][Medline]

    Voytas, D. F., M. P. Cummings, A. Konieczny, F. M. Ausubel, and S. R. Rodermel. 1992 Copia-like retrotransposons are ubiquitous among plants. Proc. Natl. Acad. Sci. USA 89:7124–7128.

    Wakamiya, I., R. J. Newton, J. S. Johnston, and H. J. Price. 1993. Genome size and environmental factors in the genus Pinus. Am. J. Bot. 80:1235–1241.

    Xiong, Y., and T. H. Eickbush. 1990. Origin and evolution of retroelements based upon their reverse transcriptase sequences. EMBO J. 9:3353–3362.[Abstract]

Accepted for publication October 4, 2000.