Retrotransposon BARE-1: expression of encoded proteins and formation of virus-like particles in barley cells

Arabidopsis retrotransposon virus-like particles and their regulation by epigenetically activated small RNA

In Arabidopsis, LTR retrotransposons are activated by mutations in the chromatin gene DECREASE in DNA METHYLATION 1 (DDM1), giving rise to 21- to 22-nt epigenetically activated siRNA (easiRNA) that depend on RNA DEPENDENT RNA POLYMERASE 6 (RDR6). We purified virus-like particles (VLPs) from ddm1 and ddm1rdr6 mutants in which genomic RNA is reverse transcribed into complementary DNA. High-throughput short-read and long-read sequencing of VLP DNA (VLP DNA-seq) revealed a comprehensive catalog of active LTR retrotransposons without the need for mapping transposition, as well as independent of genomic copy number. Linear replication intermediates of the functionally intact COPIA element EVADE revealed multiple central polypurine tracts (cPPTs), a feature shared with HIV in which cPPTs promote nuclear localization. For one member of the ATCOPIA52 subfamily (SISYPHUS), cPPT intermediates were not observed, but abundant circular DNA indicated transposon “suicide” by auto-integration within the VLP. easiRNA targeted EVADE genomic RNA, polysome association of GYPSY (ATHILA) subgenomic RNA, and transcription via histone H3 lysine-9 dimethylation. VLP DNA-seq provides a comprehensive landscape of LTR retrotransposons and their control at transcriptional, post-transcriptional, and reverse transcriptional levels.

Interspecific hybridisation and LTR-retrotransposon mobilisation-related structural variation in plants: A case study

The dynamics of long-terminal-repeat retrotransposons in two poplar species (Populus deltoides and P. nigra) and in an interspecific hybrid, recently synthesized, were investigated by analyzing the genomic abundance and transcription levels of a collection of 828 full-length retroelements identified in the genome sequence of P. trichocarpa, all occurring also in the genomes of P. deltoides and P. nigra. Overall, genomic abundance and transcription levels of many retrotransposons in the hybrid resulted higher or lower than expected by calculating the mean of the parental values. A bioinformatics procedure was established to ascertain the occurrence of the same retrotransposon loci in the three genotypes. The results indicated that retrotransposon abundance variations between the hybrid and the mean value of the parents were due to i) co-segregation of retrotransposon high- or low-abundant haplotypes; ii) new retroelement insertions; iii) retrotransposon loss. Concerning retrotransposon expression, this was generally low, with only 14/828 elements over- or under-expressed in the hybrid than expected by calculating the mean of the parents. It is concluded that interspecific hybridisation between the two poplar species determine quantitative variation and differential expression of some retrotransposons, with possible consequences for the genetic differentiation of the hybrid.

Retrotransposon BARE displays strong tissue-specific differences in expression

The BARE retrotransposon comprises c. 10% of the barley (Hordeum vulgare) genome. It is actively transcribed, translated and forms virus-like particles (VLPs). For retrotransposons, the inheritance of new copies depends critically on where in the plant replication occurs. In order to shed light on the replication strategy of BARE in the plant, we have used immunolocalization and in situ hybridization to examine expression of the BARE capsid protein, Gag, at a tissue-specific level. Gag is expressed in provascular tissues and highly localized in companion cells surrounding the phloem sieve tubes in mature vascular tissues. BARE Gag and RNA was not seen in the shoot apical meristem of young seedlings, but appeared, following transition to flowering, in the developing floral spike. Moreover, Gag has a highly specific localization in pre-fertilization ovaries. The strong presence of Gag in the floral meristems suggests that newly replicated copies there will be passed to the next generation. BARE expression patterns are consistent with transcriptional regulation by predicted response elements in the BARE promoter, and in the ovary with release from epigenetic transcriptional silencing. To our knowledge, this is the first analysis of the expression of native retrotransposon proteins within a plant to be reported.

Retrotransposon replication in plants

Retrotransposons comprise the bulk of large plant genomes, replicating via an RNA intermediate whereby the original, integrated element remains in place. Of the two main orders, the LTR retrotransposons considerably outnumber the LINEs. LINEs integrate into target sites simultaneously with the RNA transcript being copied into cDNA by target-primed reverse transcription. LTR retrotransposon replication is basically equivalent to the intracellular phase of retroviral life cycles. The envelope gene giving extracellular mobility to retroviruses is in fact widespread in plants and their retrotransposons. Evolutionary analyses of the retrotransposons and retroviruses suggest that both form an ancient monophyletic group. The particular adaptations of LTR retrotransposons to plant life cycles enabling their success remain to be clarified.

BARE retrotransposons are translated and replicated via distinct RNA pools

The replication of Long Terminal Repeat (LTR) retrotransposons, which can constitute over 80% of higher plant genomes, resembles that of retroviruses. A major question for retrotransposons and retroviruses is how the two conflicting roles of their transcripts, in translation and reverse transcription, are balanced. Here, we show that the BARE retrotransposon, despite its organization into just one open reading frame, produces three distinct classes of transcripts. One is capped, polyadenylated, and translated, but cannot be copied into cDNA. The second is not capped or polyadenylated, but is destined for packaging and ultimate reverse transcription. The third class is capped, polyadenylated, and spliced to favor production of a subgenomic RNA encoding only Gag, the protein forming virus-like particles. Moreover, the BARE2 subfamily, which cannot synthesize Gag and is parasitic on BARE1, does not produce the spliced sub-genomic RNA for translation but does make the replication competent transcripts, which are packaged into BARE1 particles. To our knowledge, this is first demonstration of distinct RNA pools for translation and transcription for any retrotransposon.

The application of LTR retrotransposons as molecular markers in plants

Retrotransposons are a major agent of genome evolution. Various molecular marker systems have been developed that exploit the ubiquitous nature of these genetic elements and their property of stable integration into dispersed chromosomal loci that are polymorphic within species. The key methods, SSAP, IRAP, REMAP, RBIP, and ISBP, all detect the sites at which the retrotransposon DNA, which is conserved between families of elements, is integrated into the genome. Marker systems exploiting these methods can be easily developed and inexpensively deployed in the absence of extensive genome sequence data. They offer access to the dynamic and polymorphic, nongenic portion of the genome and thereby complement methods, such as gene-derived SNPs, that target primarily the genic fraction.

BARE retrotransposons produce multiple groups of rarely polyadenylated transcripts from two differentially regulated promoters

The BARE retrotransposon family comprises more than 10(4) copies in the barley (Hordeum vulgare) genome. The element is bounded by long terminal repeats (LTRs, 1829 bp) containing promoters and RNA-processing motifs required for retrotransposon replication. Members of the BARE1 subfamily are transcribed, translated, and form virus-like particles. Very similar retrotransposons are expressed as RNA and protein in other cereals and grasses. The BARE2 subfamily is, however, non-autonomous because it cannot produce the GAG capsid protein. The pattern of plant development implies that inheritance of integrated copies should critically depend, in the first instance, on cell-specific and tissue-specific expression patterns. We examined transcription of BARE within different barley tissues and analyzed the promoter function of the BARE LTR. The two promoters of the LTR vary independently in activity by tissue. In embryos TATA1 was almost inactive, whereas transcription in callus appears to be less tightly regulated than in other tissues. Deletion analyses of the LTR uncovered strong positive and negative regulatory elements. The promoters produce multiple groups of transcripts that are distinct by their start and stop points, by their sequences, and by whether they are polyadenylated. Some of these groups do not share the common end structures needed for template switching during replication. Only about 15% of BARE transcripts are polyadenylated. The data suggest that distinct subfamilies of transcripts may play independent roles in providing the proteins and replication templates for the BARE retrotransposon life cycle.

Virus-like particle formation and translational start site choice of the plant retrotransposon Tto1

Ty1/copia group retrotransposon Tto1 from tobacco was put under control of an inducible promoter for expression in Arabidopsis thaliana. The system was used to analyze intermediates of the transposition process. The Tto1 RNA 5' region has a complex structure and contains several AUG codons. We therefore sought to experimentally define the translation initiation site. Constructs starting at various positions within the structural gag region were expressed in planta and functionally characterized. We found that gag proteins starting at the first ATG of the gag-pol ORF (ATG1), but also those starting at the third ATG of the gag-pol ORF (ATG3), can form virus-like particles (VLPs). However, gag protein expressed by the inducible Tto1 element had a size similar to gag starting at ATG1, and mutation of ATG1 in the inducible element abolished reverse transcription. This suggested that translation initiation at ATG1 is essential for the Tto1 life cycle. To support this conjecture, gag protein starting at ATG1, or gag protein shortened amino-terminally by nine amino acids (starting at the second ATG of the gag region, ATG2), was co-expressed with Tto1 carrying mutations at ATG1 and ATG2. Trans-complementation of the defective Tto element by gag starting at ATG1, but not by gag starting at ATG2, defines ATG1 as the functional translation initiation site.

Parasitism and the retrotransposon life cycle in plants: a hitchhiker's guide to the genome

LTR (long terminal repeat) retrotransposons are the main components of higher plant genomic DNA. They have shaped their host genomes through insertional mutagenesis and by effects on genome size, gene expression and recombination. These Class I transposable elements are closely related to retroviruses such as the HIV by their structure and presumptive life cycle. However, the retrotransposon life cycle has been closely investigated in few systems. For retroviruses and retrotransposons, individual defective copies can parasitize the activity of functional ones. However, some LTR retrotransposon groups as a whole, such as large retrotransposon derivatives and terminal repeats in miniature, are non-autonomous even though their genomic insertion patterns remain polymorphic between organismal accessions. Here, we examine what is known of the retrotransposon life cycle in plants, and in that context discuss the role of parasitism and complementation between and within retrotransposon groups.

Quantification of the retrotransposon BARE-1 reveals the dynamic nature of the barley genome

We used quantitative real-time PCR analysis to measure the copy number of the BARE-1 retrotransposon in 5 cultivars of barley (Hordeum vulgare), as well as in samples from its wild relative, Hordeum spontaneum. Two sets of PCR primers were used to amplify regions within the long terminal repeat (LTR) and the reverse transcriptase (RT) gene of BARE-1 (GenBank accession Z17327). The LTR primers detected an average of 2.148 x 105 +/- 0.012 x 105 copies per haploid genome among barley samples, whereas the RT primers detected an average of 1.588 x 104 +/- 0.085 x 104 copies. The average ratio of LTR:RT was estimated to be 13.5:1. This finding indicates that more than 7% of the barley genome is occupied by BARE-1 elements in the form of solo LTRs and another 2.6% of the genome is occupied by the full-length element. Taken together, BARE-1 sequences represent approximately 9.6% of the barley genome among the barley plants used in this study. For the above estimation, a genome size of 5.44 x 103 Mb for H. vulgare and 5.39 x 103 Mb for H. spontaneum were assumed. Our study on quantification results of the BARE-1 for a small group of barley cultivars showed that there are significant differences among cultivars in terms of BARE-1 copy number, providing further evidence that BARE-1 is active and has a major role in shaping the barley genome as a result of breeding and selection. Quantification results also showed that most of the elements (> 90%) are present as truncated copies (solo LTRs). These results show that there is a high level of recombination leading to the formation of truncated elements and a subsequent DNA loss from the genome. Taken together, our study provides a glimpse into a dynamic micro-evolutionary process that is the by-product of genome reshuffling and directional selection in barley breeding programs.

Retrotransposons: Metaparasites and Agents of Genome Evolution

Transposable elements comprise the bulk of higher plant genomes. The majority of these elements are the Class I LTR retrotransposons, which transpose via an RNA intermediate in a "Copy-and-Paste" mechanism. Because retrotransposons use cellular resources and their own enzymes to replicate independently of the genome as a whole, and have thereby become in many cases more predominant than the cellular genes, they have been considered "selfish DNA" and nuclear parasites. They are thought to share many features of the internal life cycle of retroviruses such as HIV (lentiviruses). However, whereas at least some of the retroviruses arriving in an organism during an infection must be functional in order for the infection to proceed, some LTR retrotransposon families appear to completely lack active members even though they remain mobile. Furthermore, the process of retrotransposition is inherently error-prone and mutagenic, giving rise to "pseudospecies," or clusters of imperfect copies. The non-autonomous retrotransposons are able to cis- and trans-parasitize host retrotransposons to gain mobility, much as do defective interfering particles of RNA viruses. Hence, a complex dynamic is set up, whereby the impact of retrotransposons on genomes can be under selection on the organismal level; the impact of non-autonomous retrotransposons on autonomous ones can likewise be under selection if there is selection on the autonomous elements themselves. We are exploring the retrotransposon life cycle and the causes and possible consequences of non-autonomy at each stage regarding genome evolution.

A movable feast: diverse retrotransposons and their contribution to barley genome dynamics

Cellular genes comprise at most 5% of the barley genome; the rest is occupied primarily by retrotransposons. Retrotransposons move intracellularly by a replicative mechanism similar to that of retroviruses. We describe the major classes of retrotransposons in barley, including the two nonautonomous groups that were recently identified, and detail the evidence supporting our current understanding of their life cycle. Data from analyses of long contiguous segments of the barley genome, as well as surveys of the prevalence of full-length retrotransposons and their solo LTR derivatives in the genus Hordeum, indicate that integration and recombinational loss of retrotransposons are major factors shaping the genome. The sequence conservation and integrative capacity of barley retrotransposons have made them excellent sources for development of molecular marker systems.

Genetic diversity among barley cultivars assessed by sequence-specific amplification polymorphism

We analyzed the genetic structure and relationships among barley cultivars (Hordeum vulgare L.) with sequence-specific amplification polymorphisms (S-SAPs). Polymorphisms were identified in 824 individual barley plants representing 103 cultivars (eight plants per cultivar) widely grown in Canada and the United States, using PCR primers designed from the long terminal repeat of the barley retrotransposon BARE-1 and a subset of four selective MseI primers. From the 404 bands scored, 150 were polymorphic either within or between cultivars. Genetic structure assessed with analysis of molecular variance attributed the largest component of variation to the within groups of cultivars (69-86%). Within-cultivar genetic variation was estimated as average gene diversity over loci and ranged from 0 (completely homogenous) to 0.076 (most heterogeneous cultivar). Only 17 out of 103 cultivars (16%) were judged to be homogenous by this criterion. Relationships among cultivars were analyzed by cluster analysis using unweighted pair-groups using arithmetic averages and found groups similar to those determined by agriculturally significant phenotypic traits such as spike morphology (two-rowed or six-rowed), cultivar type (malting or feed), seed characteristic (hull-less or hulled), and growth habit (winter or spring), with minor overlaps. Discriminant analysis of groups determined by these phenotypic traits fully supported the different groups with minor overlaps between the malting/feed. S-SAP markers generated from retrotransposons such as BARE-1 are invaluable tools for the study of genetic diversity in organisms with a narrow genetic base such as barley. In this study, S-SAP analysis revealed significant amounts of cryptic variation in closely related cultivars including somaclonal variation, which could not be inferred by the pedigree analysis.

The DNA sequence of a gypsy element from Gossypium hirsutum L. and characterization of gypsy elements in three Gossypium species

A strategy was developed to isolate a complete gypsy-element from Gossypium hirsutum L. based on the sequence of a RAPD that was polymorphic in near isoganic lines of cotton that varied in leaf shape. A 5998 nt clone was isolated and its gene order and sequence confirmed it was a gypsy-type retroelement. Sequences homologous to the gag portion of this clone were also found in Gossypium herbaceum L. and Gossypium raimondii L. A portion of the open reading frame of the integrase gene was amplified from three Gossypium species under investigation. PCR products of the expected size were amplified from all three Gossypium species. G. raimondii sequences were statistically more degenerate than sequences from either G. hirsutum or G. herbaceum. Finally, analysis of the open reading frames of selected integrase clones from the three Gossypium species revealed a second gypsy element based on similarity of the deduced amino acid sequence.

Comparison of the utility of barley retrotransposon families for genetic analysis by molecular marker techniques

The Sequence-Specific Amplification Polymorphism (S-SAP) method, and the related molecular marker techniques IRAP (inter-retrotransposon amplified polymorphism) and REMAP (retrotransposon-microsatellite amplified polymorphism), are based on retrotransposon activity, and are increasingly widely used. However, there have been no systematic analyses of the parameters of these methods or of the utility of different retrotransposon families in producing polymorphic, scorable fingerprints. We have generated S-SAP, IRAP, and REMAP data for three barley (Hordeum vulgare L.) varieties using primers based on sequences from six retrotransposon families (BARE-1, BAGY-1, BAGY-2, Sabrina, Nikita and Sukkula). The effect of the number of selective bases on the S-SAP profiles has been examined and the profiles obtained with eight MseI+3 selective primers compared for all the elements. Polymorphisms detected in the insertion pattern of all the families show that each can be used for S-SAP. The uniqueness of each transposition event and differences in the historic activity of each family suggest that the use of multiple retrotransposon families for genetic analysis will find applications in mapping, fingerprinting, and marker-assisted selection and evolutionary studies, not only in barley and other Hordeum species and related taxa, but also more generally.

Plant LTR-retrotransposons and MITEs: control of transposition and impact on the evolution of plant genes and genomes

Transposons are genetic elements that can move, and sometimes spread, within genomes, and that constitute an important fraction of eukaryote genomes. Two types of transposons, long terminal repeat (LTR)-retrotransposons and miniature inverted-repeat transposable elements (MITEs), are highly represented in plant genomes, and can account for as much as 50-80% of the total DNA content. In the last few years it has been shown that, in spite of their mutagenic capacity, both LTR-retrotransposons and MITEs can be found associated to genes, suggesting that their activity has influenced the evolution of plant genes. In this review we will summarise recent data on the control of the activity and the impact of both LTR-retrotransposons and MITEs on the evolution of plant genes and genomes.

Dasheng: a recently amplified nonautonomous long terminal repeat element that is a major component of pericentromeric regions in rice

A new and unusual family of LTR elements, Dasheng, has been discovered in the genome of Oryza sativa following database searches of approximately 100 Mb of rice genomic sequence and 78 Mb of BAC-end sequence information. With all of the cis-elements but none of the coding domains normally associated with retrotransposons (e.g., gag, pol), Dasheng is a novel nonautonomous LTR element with high copy number. Over half of the approximately 1000 Dasheng elements in the rice genome are full length (5.6-8.6 kb), and 60% are estimated to have amplified in the past 500,000 years. Using a modified AFLP technique called transposon display, 215 elements were mapped to all 12 rice chromosomes. Interestingly, more than half of the mapped elements are clustered in the heterochromatic regions around centromeres. The distribution pattern was further confirmed by FISH analysis. Despite clustering in heterochromatin, Dasheng elements are not nested, suggesting their potential value as molecular markers for these marker-poor regions. Taken together, Dasheng is one of the highest-copy-number LTR elements and one of the most recent elements to amplify in the rice genome.

Plant transposable elements: where genetics meets genomics

Transposable elements are the single largest component of the genetic material of most eukaryotes. The recent availability of large quantities of genomic sequence has led to a shift from the genetic characterization of single elements to genome-wide analysis of enormous transposable-element populations. Nowhere is this shift more evident than in plants, in which transposable elements were first discovered and where they are still actively reshaping genomes.

Envelope-class retrovirus-like elements are widespread, transcribed and spliced, and insertionally polymorphic in plants

Retrotransposons and retroviruses share similar intracellular life cycles and major encoded proteins, but retrotransposons lack the envelope (env) critical for infectivity. Retrotransposons are ubiquitous and abundant in plants and active retroviruses are known in animals. Although a few env-containing retroelements, gypsy-like Athila, Cyclops, and Calypso and copia-like SIRE-1, have been identified in plants, the general presence and functionality of the domain remains unclear. We show here that env-class elements are present throughout the flowering plants and are widely transcribed. Within the grasses, we show the transcription of the env domain itself for Bagy-2 and related retrotransposons, all members of the Athila group. Furthermore, Bagy-2 transcripts undergo splicing to generate a subgenomic env product as do those of retroviruses. Transcription and the polymorphism of their insertion sites in closely related barley cultivars suggests that at least some are propagationally active. The putative ENV polypeptides of Bagy-2 and rice Rigy-2 contain predicted leucine zipper and transmembrane domains typical of retroviral ENVs. These findings raise the prospect of active retroviral agents among the plants.

Linear DNA intermediates of the Tto1 retrotransposon in Gag particles accumulated in stressed tobacco and Arabidopsis thaliana

The active transcription of some plant retrotransposons under diverse stress conditions suggests active transposition. However, transposition has been demonstrated only during tissue/cell culture. To examine whether transposition is activated under conditions other than tissue/cell culture, DNA intermediates for retrotransposition of the tobacco retrotransposon Tto1 were analysed. Using transgenic Arabidopsis callus expressing high levels of Tto1 RNA in a ddm1 hypomethylation mutant background, the existence of extrachromosomal Tto1 linear DNA molecules in a Gag-particle fraction was demonstrated. By combination with ligation-mediated PCR amplification, we detected Tto1 linear DNA molecules in particle fractions from callus and methyl jasmonate-treated leaves of tobacco, but not from non-stressed leaves. Tto1 DNA intermediates could not be detected in the tobacco corolla where Tto1 is expressed. These results indicate that the transcriptional activation of Tto1 by defence-related stresses leads to the synthesis of DNA intermediates, whereas post-transcriptional suppression of Tto1 activity is suggested in the corolla.

Active retrotransposons are a common feature of grass genomes

A large fraction of the genomes of grasses, members of the family Graminae, is composed of retrotransposons. These elements resemble animal retroviruses in their structure and possess a life cycle similar to theirs that includes transcription, translation, and integration of daughter copies. We have investigated if retrotransposons are generally transcribed in the grasses and other plants, and whether the various families of elements are translationally and integrationally active in multiple grass species. A systematic search of 7.8 x 10(5) publicly available expressed sequence tags from plants revealed widespread retrotransposon transcripts at a frequency of one in 1,000. Monocot retrotransposons found relatively more expressed sequence tags from non-source species than did those of dicots. Antibodies were raised to the capsid protein, GAG, of BARE-1, a transcribed and translated copia-like retrotransposon of barley (Hordeum vulgare). These detected immunoreactive proteins of sizes identical to those of the BARE-1 GAG and polyprotein, respectively, in other species of the tribe Triticeae as well as in oats (Avena sativa) and rice (Oryza sativa). Retrotransposon-based markers showed integrational polymorphisms for BARE-1 in different subfamilies of the Graminae. The results suggest that grasses share families of transcriptionally, translationally, and integrationally active retrotransposons, enabling a comparative and integrative approach to understanding the life cycle of retrotransposons and their impact on the genome.

A contiguous 66-kb barley DNA sequence provides evidence for reversible genome expansion

Organisms with large genomes contain vast amounts of repetitive DNA sequences, much of which is composed of retrotransposons. Amplification of retrotransposons has been postulated to be a major mechanism increasing genome size and leading to "genomic obesity." To gain insights into the relation between retrotransposons and genome expansion in a large genome, we have studied a 66-kb contiguous sequence at the Rar1 locus of barley in detail. Three genes were identified in the 66-kb contig, clustered within an interval of 18 kb. Inspection of sequences flanking the gene space unveiled four novel retroelements, designated Nikita, Sukkula, Sabrina, and BAGY-2 and several units of the known BARE-1 element. The retroelements identified are responsible for at least 15 integration events, predominantly arranged as multiple nested insertions. Strikingly, most of the retroelements exist as solo LTRs (Long Terminal Repeats), indicating that unequal crossing over and/or intrachromosomal recombination between LTRs is a common feature in barley. Our data suggest that intraelement recombination events deleted most of the original retrotransposon sequences, thereby providing a possible mechanism to counteract retroelement-driven genome expansion.

Genome evolution of wild barley (Hordeum spontaneum) by BARE-1 retrotransposon dynamics in response to sharp microclimatic divergence

The replicative spread of retrotransposons in the genome creates new insertional polymorphisms, increasing retrotransposon numbers and potentially both their share of the genome and genome size. The BARE-1 retrotransposon constitutes a major, dispersed, active component of Hordeum genomes, and BARE-1 number is positively correlated with genome size. We have examined genome size and BARE-1 insertion patterns and number in wild barley, Hordeum spontaneum, in Evolution Canyon, Lower Nahal Oren, Mount Carmel, Israel, along a transect presenting sharply differing microclimates. BARE-1 has been sufficiently active for its insertional pattern to resolve individuals in a way consonant with their ecogeographical distribution in the canyon and to distinguish them from provenances outside the canyon. On both slopes, but especially on the drier south-facing slope, a simultaneous increase in the BARE-1 copy number and a decrease in the relative number lost through recombination, as measured by the abundance of solo long terminal repeats, appear to have driven the BARE-1 share of the genome upward with the height and dryness of the slope. The lower recombinational loss would favor maintenance of more full-length copies, enhancing the ability of the BARE-1 family to contribute to genome size growth. These local data are consistent with regional trends for BARE-1 in H. spontaneum across Israel and therefore may reflect adaptive selection for increasing genome size through retrotransposon activity.