Department of Biochemistry, University of Dundee, Scotland
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Ty1-copia group retrotransposons are present as large, highly heterogeneous populations within all higher plant divisions (Flavell et al. 1992
; Flavell, Smith, and Kumar 1992
; Voytas et al. 1992
). 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. 1992
; Flavell, Smith, and Kumar 1992
; Voytas et al. 1992
). 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 1993
; Capy, Anxolabéhere, and Langin 1994
; Flavell et al. 1997
). 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 1994
). Subsequent phylogenetic studies have failed to provide convincing evidence for horizontal transfer of plant retrotransposons (Vanderwiel, Voytas, and Wendel 1993
; Gribbon et al. 1999
; Matsuoka and Tsunewaki 1999
), 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 1999
; Jordan, Matyunina, and McDonald 1999
), 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 |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
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 1994
). 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)
. 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.
|
|
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
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. 1996
).
|
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. 1992
) 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 1999
; Jordan, Matyunina, and McDonald 1999
).
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 1997
), 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.
|
|
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
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 1993
). Modern conifers arose in the Triassic (>200 MYA). Virtually all phylogenetic analyses place Ginkgo basal to the conifer lineage (Crane 1985
; Nixon et al. 1994
; Qiu et al. 1999
; Soltis, Soltis, and Chase 1999
). 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. 1999
) 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 1999
), 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 1999
; Jordan, Matyunina, and McDonald 1999
). 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 1999
). 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 1997
; Gonzales and Lessios 1999
). However, the large majority of the plant retrotransposons studied previously show much lower conservation across large evolutionary divides (Flavell et al. 1992
; Flavell, Smith, and Kumar 1992
; Voytas et al. 1992
; Matsuoka and Tsunewaki 1997, 1999
; Gribbon et al. 1999
). 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. 1999
; 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 1994
) 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. 1998
), 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. 1993
). 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 1994
), 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. 1992
; Flavell, Smith, and Kumar 1992
; Voytas et al. 1992
). 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. 1999
) 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
.
|
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
1 Present address: University of Helsinki, Institute of Biotechnology, Helsinki, Finland.
1 Keywords: retrotransposon
gymnosperm
Picea
Ginkgo
evolution
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
![]() |
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:712.
Charlesworth B., P. Sniegowski, and W. Stephan. 1994. The evolutionary dynamics of repetitive DNA in eukaryotes. Nature 371:215220.
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:363375.[ISI][Medline]
Crane, P. R. 1985. Phylogenetic relationships in seed plants. Cladistics 1:329348.
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:36393644.[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:185195.
Flavell, A. J., D. B. Smith, and A. Kumar. 1992. Extreme heterogeneity of TyTy1-copia group retrotransposons in plants. Mol. Gen. Genet. 231:233242.[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:938952.[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:883891.[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:1262112625.
Jukes, T. H., and C. R. Cantor. 1969. Evolution of protein molecules. Pp. 21132 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:27082713.
Konieczny, A., D. F. Voytas, M. P. Cummings, and F. M. Ausubel. 1991. A superfamily of Arabidopsis thaliana retrotransposons. Genetics 127:801809.
Kumar, A., and J. L. Bennetzen. 1999. Plant retrotransposons. Annu. Rev. Genet. 33:479532.[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:6980.[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:793805.[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:335343.[ISI][Medline]
. 1999. Evolutionary dynamics of Ty1-copia group retrotransposons in grass shown by reverse transcriptase domain analysis. Mol. Biol. Evol. 16:208217.[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:484533.
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:617625.
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:305315.[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:404407.
SanMiguel, P., B. S. Gaut, A. Tikhonov, Y. Nakijama, and J. L. Bennetzen. 1998. The paleontology of intergene retrotransposons of maize. Nat. Genet. 20:4345.[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:765768.
Soltis, P. S., D. E. Soltis, and M. W. Chase. 1999. Nature 402:402404.
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:46734680.[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:429447.[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:71247128.
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:12351241.
Xiong, Y., and T. H. Eickbush. 1990. Origin and evolution of retroelements based upon their reverse transcriptase sequences. EMBO J. 9:33533362.[Abstract]