Background
LTR Retrotransposons transpose through reverse transcription of an RNA intermediate and are ubiquitous components of all eukaryotic genomes thus far examined 1. Plant genomes, in particular, have been found to be comprised of a remarkably high number of LTR retrotransposons 234. For example, more than half of the maize 5 and over 90% of the wheat 6 genomes are comprised of LTR retrotransposons. There is a significant body of direct and indirect evidence that LTR retrotransposons have contributed to gene and genome evolution in both animals and plants 478910111213.
Rice (Oryza sativa L.) is a staple food for over half of the world's population. Unlike most other cereal grasses, rice has a relatively small genome estimated at 430 Mb 1415. Comparative genetic maps within the grass family indicate extensive regions of conserved gene content and orders 16. The development of many important tools for genetic analysis, including excellent genetic and physical maps 1718, efficient genetic transformation techniques 19, an ever increasing dataset of expressed sequence tags (ESTs) 20, and a large collection of diverse germplasm 21 has established rice as an ideal model for the study of mechanisms underlying plant evolution. The recent release of draft genome sequences of two rice subspecies, indica 22 and japonica 23, and ongoing efforts to compile a complete japonica rice genome sequence publicly available by International Rice Genome Sequencing Project (IRGSP) 15, promises to greatly facilitate the study of rice genome evolution.
Most TEs were earlier discovered by chance or by limited assays using conserved regions of retrotransposons in rice 242526. For example, based on the screening of genomic libraries with a conserved retrotransposon probe, Hirochika and his colleagues 24 made an estimate of ~1000 retrotransposons totally in the rice genome, which fall into 32 families. By hybridizing a rice BAC library with a variety of RT probes, however, Wang et al. 27 estimated that the genome contains only about 100 copia-like elements in the entire haploid genome. Computer based sequence similarity searches of the nucleic acid databases are considered to be a more accurate method to identify retrotransposon families and estimate TE content and distribution in rice genome 28293031. Using a new data-mining program (LTR retrotransposon structure program), we previously mined the GeneBank rice database (GBRD) as well as a more extensive Monsanto rice data set (MRD) (259 Mb) of LTR retrotransposons 32. Including numerous LTR retrotransposon families that have not been previously reported, our comprehensive survey indicated that there are at least 59 distinct O. sativa LTR retrotransposon families comprising ≈17% of the rice genome. Although the majority of rice retrotransposon families have been identified and their phylogenetic relationships have been inferred based on the sequence similarity of the RT domains, the evolutionary dynamics, substructure, and relative integration time or age of rice retrotransposon elements remain poorly understood.
Our laboratory has long been interested in the evolution of LTR retrotransposons and their contribution to genome evolution. The current availability of a considerable part of the rice sequences at NCBI database is providing an excellent and timely opportunity to analyze the evolutionary history of rice LTR-retrotransposons. The present study was undertaken (i) to assess retrotransposon diversity and timing of insertion events by identifying all distinguishable rice LTR-retrotransposon sequences in the published sequence database, and(ii) to further infer, based on the establishment of substructure of all characterized families using the sequence divergence among the more rapidly evolving LTRs, possible retrotransposon evolution. In order to gain insight into the chromosomal distribution of LTR-retrotransposons in the rice genome, we also selected chromosome 10 to conduct a detailed examination of the intra-chromosomal distribution of LTR-retrotransposons. Finally, we have made a preliminary attempt to assess the potential contribution of LTR-retrotransposons to the evolution of gene structure and function in rice by identifying elements located in or near putative genes.
Results and discussion
Characterization of O. sativa LTR retrotransposons
LTR retrotransposons displaying ≥ 90% reverse transcriptase (RTs) pairwise identity at the amino acid level were, by earlier convention, assigned to the same family 33. Using this criterion, McCarthy et al. 32 previously identified and characterized 55 families of LTR retrotransposons in the rice genome. Four additional groups of non-autonomous LTR retrotransposons identified in this previous study did not display a RT sequence homology but were, nevertheless, designated as families based on their distinct structures 32. Only 38 retrotransposon families for which both LTRs have been identified 32 were included in this study.
In this study, we have extended the earlier survey of full-length LTR retrotransposons in the rice genome to include all identifiable fragments of LTR retrotransposon sequences. In accord with previously established criteria 34, we classified rice LTR retrotransposon sequences into three major groups: 1) full-length elements: a)
autonomous full-length elements
contain all of the characteristic features of LTR retrotransposons including putative gag, pol and, in some cases, env genes flanked by LTRs; b) non-autonomous full-length elements
are a recently identified 3239 sub-class of full-length elements that have two LTRs flanking a series of repeating motifs of various lengths. Although numerous, relatively small and, in some cases, recently transposed, these non-autonomous elements encode none of the typical LTR retrotransposon ORFs and must acquire essential reverse transcription and intergration functions in trans; 2) solo LTRs are solitary LTR elements believed to be the products of recombination events between the flanking LTRs of full-length elements 1; 3) fragmented elements are defined as partially deleted or truncated LTR retrotransposon sequence. This third category is a "catch all" grouping that includes all LTR retrotransposon sequences that are neither full-length elements or solo LTRs. For information purposes, we report in this paper both the total number of fragmented elements, as well as, the number of these that are fragmented solo LTRs.
We identified a total of 1219 LTR retrotransposon sequences in the portion of the rice genome surveyed in this study distributed over 38 families (Tables 1, see Additional file: 1; Table 2). Gypsy-like elements are • 4 × more numerous than copia-like elements (990 vs. 229). The number of element sequences within families was found to vary considerably, ranging from 1(Osr19) to 255 (Osr34) (Table 2). We identified 217 full-length elements representing 17.8 % of all LTR retrotransposon sequences identified. Osr27 contains the largest number of full-length elements (61) in the genome, while 12 families (Osr5–7, 14, 19–20, 22–23, 28, 42–44) have only 1 full-length element. The most abundant copia-like family was Osr8 (59 sequences) while Osr34 is the most abundant gypsy-like family (255 sequences).
<p>Additional file 1</p>
Click here for file
<p>Table 2</p>
Number of full-length, fragmented, and solo LTRs in the sequenced O. sativa (Nipponbare) genome
LTR Retrotransposon Families
Full-length
Solo-LTRs
Fragmented retroelements (Fragmented solo-LTRs)
Subtotal
Osr1
3
4
11 (6)
18
Osr2
2
0
0
2
Osr3
0
0
3 (1)
3
Osr4
3
3
1 (1)
7
Osr5
1
0
2 (1)
3
Osr6
1
0
2 (0)
3
Osr7
1
0
3 (2)
4
Osr8
7
10
42 (37)
59
Osr9a
0
0
0
0
Osr10
2
7
21 (18)
30
Osr11a
0
0
0
0
Osr12
3
0
0
3
Osr13
11
5
19 (12)
35
Osr14
1
9
10 (10)
20
Osr15
5
1
12 (12)
18
Osr16
3
1
0
4
Osr17
3
1
10 (5)
14
Osr18a
0
0
0
0
Osr19
1
0
0
1
Osr20
1
0
1 (1)
2
Osr21a
0
0
0
0
Osr22
1
0
1 (0)
2
Osr23
1
0
1 (0)
2
Subtotal for Copia
50
41
138 (106)
229
Osr24
3
0
0
3
Osr25
34
4
40 (26)
78
Osr26
7
2
77 (50)
86
Osr27
61
39
40 (30)
132
Osr28
1
0
2 (2)
3
Osr29
6
5
21 (19)
32
Osr30
4
17
51 (43)
72
Osr31
3
2
2 (2)
7
Osr32
0
11
60 (53)
71
Osr33
4
5
84 (76)
93
Osr34
14
27
214 (176)
255
Osr35
2
0
0
2
Osr36
2
0
0
2
Osr37
8
13
28 (22)
49
Osr38a
0
0
0
0
Osr39a
0
0
0
0
Osr40
8
10
33 (27)
51
Osr41
7
4
16 (11)
27
Osr42
1
0
3 (2)
4
Osr43
1
0
7 (6)
8
Osr44
1
0
6 (6)
7
Subtotal for Gypsy
165
139
684 (551)
990
Overall
217
180
822 (660)
1219
a. No full-length element was identified for this family in the GenBank database of O. sativa. Location given is for an example RT in the database
Of the total LTR retrotransposon sequences identified, 822 or 67.4% were found to be fragmented or truncated elements. Fragmented elements were detected in 31 families varying in number from 1 (Osr4, 20, 22, and 23) to 214 (Osr34) per family. Seven families of LTR retrotransposons (Osr2, 12, 16, 19, 24, 35, and 36) were found to contain no fragmented elements. We identified 179 solo LTRs representing 14.8 % of all LTR retrotransposon sequences identified. The solo LTRs are distributed over 21 families (Osr1, 4, 8, 13–17, 25–27, 29–34, 37, and 40–41), and range in number from 1 (Osr15 and Osr16) to 39 (Osr27) per family. The Osr19 family consists of a single full-length element without any distinguishable fragmented elements or solo LTRs.
Other distinguishing characteristics of the LTR retrotransposon sequences identified in this study are presented in Table 1 (see Additional file: 1), including clone accession numbers, chromosomal location, sequence length, target site repeats (TSRs), LTR pairwise identities and estimated element age. TSRs result from a duplication of the unoccupied insertion site following element insertion 35. The TSRs of all rice LTR retrotransposon are five bp long except for Osr26.1 which has a TSR of seven bp. While members of O. sativa families share >90% RT sequence identity 32, sequence identity values among the more rapidly evolving LTRs are highly variable within families, ranging from 75.505% to 100%.
Phylogenetic substructure of O. sativa LTR retrotransposons
The slowly evolving RT encoding region of LTR retrotransposons is ideal for calculating evolutionary distances among even distantly related families of retroelements 636. However, sequence analyses of the more rapidly evolving LTRs are better suited for the characterization of phylogenetic substructure within families of LTR retrotransposons. We used LTR sequence divergence among elements within each family to identify sub-structure within the O. sativa LTR retrotransposon (Osr) phylogeny.
Significant sub-structure was detected in eleven families by neighbor-joining and parsimony criteria (Osr4, 10, 12, 25, 26, 27, 30, 34, 35, 43, and 44).
The Osr26 family, for example, is comprised of 32 elements falling into at least 5 distinct clades (Fig. 1). The Osr27 family, which is composed of 134 elements, displayed complicated substructure, consisting of at least 10 divergent clades with strong bootstrap support (Fig. 2). The remaining 27 families are closely related and displayed no significant intra-family substructure.
<p>Figure 1</p>
Phylogenetic trees of subfamily structure based on LTR nucleotide sequence data
Phylogenetic trees of subfamily structure based on LTR nucleotide sequence data. The 32 elements of Osr26 family fall into at least 5 clades, with Osr25 as the outgroup. To better exhibit tree structure, all Osr27 elements were removed. Insertions/deletions were ignored while performing phylogenetic analyses. Values on individual branches are bootstrap percentages using 1000 bootstrap repetitions. Each LTR in the tree is named by the genomic clone in which it was found. For elements with two LTRs, the 3' LTR is labeled by a lower case "b" while the 5' LTR is labeled by a lower case "a". Each tree is exhibited with a scale bar determined by the number of nucleotide substitutions per site between two sequences. The tight clustering seen in both families represents a high degree of nucleotide identity between elements within a subfamily.
![]()
<p>Figure 2</p>
Phylogenetic trees of subfamily structure based on LTR nucleotide sequence data
Phylogenetic trees of subfamily structure based on LTR nucleotide sequence data. The Osr27 family forms at least 10 divergent clades with strong bootstrap support, with Osr41 as the outgroup. To better exhibit tree structure, all Osr27 elements were removed. Insertions/deletions were ignored while performing phylogenetic analyses. Values on individual branches are bootstrap percentages using 1000 bootstrap repetitions. Each LTR in the tree is named by the genomic clone in which it was found. For elements with two LTRs, the 3' LTR is labeled by a lower case "b" while the 5' LTR is labeled by a lower case "a". Each tree is exhibited with a scale bar determined by the number of nucleotide substitutions per site between two sequences. The tight clustering seen in both families represents a high degree of nucleotide identity between elements within a subfamily.
![]()
Aging of the LTR-retrotransposons of O. sativa
LTR nucleotide identity can be used to estimate the time of LTR retrotransposon integration 373839. Of the 217 full-length elements that we have identified in O. sativa, 116 (53.5%) have >99% LTR similarity with 34 of these (15.7%) being identical. The remaining 101(46.5%)have relatively low levels of nucleotide divergence (<99%). Among them, 31(14%), 24(11%), 30(14%), and 15(7%) fell into different ranges of LTR similarities of 98–99%, 97–98%, 95–97%, and <95%, respectively. Identical LTRs imply that the elements have inserted recently and have not had time to accumulate mutations between LTRs. Using the average neutral substitution rate for grasses 373940 of 0.065 substitutions per site per million years and the intra-element LTR divergence calculated here, we have computed the estimated integration time of the 181 elements with LTR nucleotide divergence (Table 1, see Additional file: 1). Based on the full-length elements available in this study, we estimate over half of rice retroelements (53.5%) with >99% LTR identity have inserted within the last 770,000 years while the older elements (46.5%) have integrated over a span of 770,000–7,700,000 years ago. Figure 3, 4 and 5 show estimated ages of all the full-length elements found in our study.
<p>Figure 3</p>
O. sativa copia-like LTR-retrotransposon age calculated using intraelement LTR nucleotide similarities
O. sativa copia-like LTR-retrotransposon age calculated using intraelement LTR nucleotide similarities. Only are those elements that contain LTR nucleotide divergence values other than zero included.
![]()
<p>Figure 4</p>
O. sativa gypsy-like LTR-retrotransposon age calculated using intraelements LTR Nucleotide Similarities (A)
O. sativa gypsy-like LTR-retrotransposon age calculated using intraelements LTR Nucleotide Similarities (A). Only are those elements that contain LTR nucleotide divergence values other than zero included.
![]()
<p>Figure 5</p>
O. sativa gypsy-like LTR-retrotransposon age calculated using intraelements LTR Nucleotide Similarities (B)
O. sativa gypsy-like LTR-retrotransposon age calculated using intraelements LTR Nucleotide Similarities (B). Only are those elements that contain LTR nucleotide divergence values other than zero included.
![]()
The LTRs of one gypsy-like (Osr27.2) and one copia-like (Osr13.8) element displayed atypically high levels of sequence divergence indicating that these elements are exceptionally old or possibly that these elements are, in fact, hybrid elements generated by homologous recombination or some other recombination process. Indeed, such inter-element recombination events have been previously documented in yeast 4142. However, in this case, recombination is unlikely since the target site duplications of these (and indeed all full-length elements identified in this study) are identical. Thus, we have no direct evidence that any of the full-length rice LTR retrotransposons analyzed in this study were generated by recombination.
Distribution of Rice LTR-retrotransposons on chromosome 10
At the time of this study, about 48% rice genomic sequences were available in the public database, which included the almost entirely sequenced chromosome 1 and 10. In an initial effort to gain insight into the intra-chromosomal distribution of O. sativa LTR retrotransposons, we selected Chromosome 10 for a more detailed analysis. Tests were conducted to determine whether LTR retrotransposon sequences were randomly distributed on Chromosome 10. We found that the average density of elements on Chromosome 10 is 22.321/Mb (Table 3). There is a nonrandom clustering of both copia-like and gypsy-like LTR retrotransposons along the chromosome (Figure 6, 7 and 8).
<p>Table 3</p>
Number of full-length, fragmented, and solo LTRs in the completely sequenced chromosome 10 of O. sativa (Nipponbare) genome
LTR Retrotransposon Families
Full-length
Solo-LTRs
Fragmented retroelements (Fragmented solo-LTRs)
Subtotal
Osr1
4
1
15 (4)
20
Osr2
0
0
0
0
Osr3
0
0
0
0
Osr4
0
0
0
0
Osr5
0
0
1(1)
1
Osr6
0
0
1(0)
1
Osr7
0
0
0
0
Osr8
3
10
40 (36)
53
Osr9a
0
0
0
0
Osr10
2
5
9 (9)
16
Osr11a
0
0
0
0
Osr12
6
2
2 (2)
10
Osr13
7
5
16 (14)
28
Osr14
3
8
4 (4)
15
Osr15
3
0
6 (6)
9
Osr16
1
0
0
1
Osr17
3
1
2 (1)
6
Osr18a
0
0
0
0
Osr19
1
0
0
1
Osr20
0
0
0
0
Osr21a
0
0
0
0
Osr22
1
0
0
1
Osr23
0
0
1 (0)
1
Subtotal for Copia
34
32
97 (77)
163
Osr24
0
0
0
0
Osr25
19
1
25 (19)
45
Osr26
1
0
36 (18)
37
Osr27
9
12
8 (8)
29
Osr28
0
0
1 (1)
1
Osr29
1
3
19 (15)
23
Osr30
1
7
63 (59)
71
Osr31
3
0
0
3
Osr32
0
4
25 (24)
29
Osr33
2
7
138 (115)
147
Osr34
9
11
71 (60)
91
Osr35
2
0
0
2
Osr36
1
0
0
1
Osr37
6
4
18 (14)
28
Osr38a
0
0
0
0
Osr39a
0
0
0
0
Osr40
4
5
24 (23)
33
Osr41
6
0
15 (8)
21
Osr42
0
0
1 (1)
1
Osr43
0
0
1 (0)
1
Osr44
0
0
0
0
Subtotal for Gypsy
64
54
445 (365)
563
Overall
98
86
542 (442)
726
a. No full-length element was identified for this family in the GenBank database of O. sativa. Location given is for an example RT in the database
<p>Figure 6</p>
Distribution of full-length, fragmented and solo LTR element sequences on the Chromosome 10 in the O. sativa genome
Distribution of full-length, fragmented and solo LTR element sequences on the Chromosome 10 in the O. sativa genome. The figure shows the distribution of copia-like, gypsy-like and all elements on the chromosome 10. A genomic position for an element was estimated based on a PAC or BAC clone which contains the element. All elements were plotted to their respective chromosome location. Both copia-like and gypsy-like LTR retrotransposons exhibit nonrandom clustering along the chromosome. More LTR retrotransposons reside in the pericentric regions of the chromosome than in the arms, and short arm of Chromosome 10 also displays a higher density of the elements than long arm. The number of gypsy-like elements are more than that of copia-like elements on chromosome 10.
![]()
<p>Figure 7</p>
The distribution of full-length, solo and fragmented LTR elements along the chromosome 10
The distribution of full-length, solo and fragmented LTR elements along the chromosome 10. The borders of the pericentromeric regions of chromosome 10 were assigned as being 5 cM from the center of the centromere on each arm. The remainders of the chromosome were designated as arms. Both copia-like and gypsy-like LTR retrotransposons exhibit nonrandom clustering along the chromosome. More LTR retrotransposons reside in the pericentric regions of the chromosome than in the arms, and short arm of Chromosome 10 also displays a higher density of the elements than long arm. The number of gypsy-like elements are more than that of copia-like elements on chromosome 10.
![]()
<p>Figure 8</p>
Insertion frequencies of full-length, solo and fragmented LTR elements along the chromosome 10
Insertion frequencies of full-length, solo and fragmented LTR elements along the chromosome 10. DNA density for chromosomal arms and pericentromeric regions was calculated from the total DNA of three contigs (~10.4 cM and ~62.8 cM in arms, ~10.5 cM in pericentromeric regions) on Chromosome 10 (see Methods). Both copia-like and gypsy-like LTR retrotransposons exhibit nonrandom clustering along the chromosome. More LTR retrotransposons reside in the pericentric regions of the chromosome than in the arms, and short arm of Chromosome 10 also displays a higher density of the elements than long arm. The number of gypsy-like elements are more than that of copia-like elements on chromosome 10.
![]()
There are • 3 × more LTR retrotransposon sequences in the pericentric regions of the chromosome than in the arms (48.069 vs. 16.019 elements per Mb DNA, p < 0.005). Moreover, the short arm of Chromosome 10 displays a higher density of the elements than long arm (58.163 vs. 11.183 elements per Mb DNA). Consistent with what was observed on the whole genome level, there are proportionally more gypsy-like than copia-like elements on chromosome 10 (17 vs. 5 elements per Mb DNA). Comparing total number of retrotransposons in rice genome (Table 2) and retrotransposons in Chromosome 10), we noticed that two families (Osr12 and 33) have more elements in Chromosome 10 than in the whole sequenced genome. This is not unexpected because many sequences from Chromosome 10 were released after we stopped searching rice genomic sequences in GenBank as of January 1, 2002 (see methods).
O. sativa LTR retroelement sequences are associated with putative genes
In an initial effort to determine whether and how frequently LTR retrotransposon sequences are associated with putative genes in the rice genome, we examined elements from all 19 families of copia-like elements and 5 representative families gypsy-like elements present in O. sativa. Our results indicate that 111/510 or 22% of LTR retrotransposon sequences lie within putative or established rice genes over the region of the genome analyzed in this study (Table 4, Figure 9 and 10). Fragmented elements are more frequently associated with genes (16%), followed by solo LTRs (3%) and full-length elements (2%). While these numbers are likely to change somewhat as the rice genome is better annotated, these preliminary estimates indicate that the potential contribution of LTR retrotranspson sequences to the evolution of gene structure and function in rice may be significant.
<p>Table 4</p>
Numbers and frequencies of association between LTR retrotransposons and putative genes for all known copia and some gypsy representative families in the rice (Nipponbare) genome
LTR Retrotransposon Families
Full-length
Solo LTRs
Fragment
Subtotal
Percentages of Association Cases in Total Retroelements
Osr1
0
0
3
3
3/18 = 16.67%
Osr2
0
0
0
0
0
Osr3
0
0
2
2
2/3 = 66.67%
Osr4
0
0
0
0
0
Osr5
0
0
0
0
0
Osr6
0
0
0
0
0
Osr7
0
0
0
0
0
Osr8
1
0
8
9
10/59 = 16.95%
Osr9a
0
0
0
0
0
Osr10
0
1
2
3
3/30 = 10%
Osr11a
0
0
0
0
0
Osr12
0
0
0
0
0
Osr13
1
3
2
6
6/35 = 17.14%
Osr14
0
0
0
0
0
Osr15
0
0
5
5
5/18 = 27.58%
Osr16
1
0
0
1
1/4 = 25.00%
Osr17
2
0
1
3
3/13 = 23.08%
Osr18a
0
0
0
0
0
Osr19
0
0
0
0
0
Osr20
0
0
0
0
0
Osr21a
0
0
0
0
0
Osr22
0
0
0
0
0
Osr23
0
0
0
0
0
Subtotal for Copia
5
4
23
32
32/230 = 14.04%
Osr31
1
0
0
1
1/7 = 14.29%
Osr34
7
8
57
72
72/257 = 28.24%
Osr36
1
0
0
1
1/2 = 50.00%
Osr43
1
0
3
4
4/8 = 50.00%
Osr44
0
0
1
1
1/7 = 14.29%
Subtotal for Gypsy
10
8
61
79
79/280 = 28.42%
Overall
15
12
84
111
111/510 = 21.94%
a. No full-length element was identified for this family in the GenBank database of O. sativa. Location given is for an example RT in the database
<p>Figure 9</p>
Five cases of association between retroelements with O. sativa putative genes observed in this study
Five cases of association between retroelements with O. sativa putative genes observed in this study. Red arrows indicate positions of LTR retroelements with the direction of transcription. Green bars represent NCBI database-predicted gene regions with their orientation of transcription. All the five associations (from A to E) are located in the following genomic clones: AP003103, AP003105, H0423H10, AC018727 and H0806H05. (A) An entire full-length Osr43 retrotransposon is completely a part of a putative rice gene (B1085F09.5) of unknown function. (B) An Osr13 solo LTR is entirely a part of a putative rice gene (P0489A05) of unknown function. (C) An Osr15 fragmented LTR is completely a part of an unknown gene (H0423H10.8) and overlaps two introns (No. 1 and 2) and No. 2 extron of the gene. (D) An Osr17 full-length retrotransposon overlaps two putative genes nearby (OSJNBa0056G17.33 and OSJNBa0056G17.22) of unknown function. (E) An Osr34 fragmented LTR is associated with two putative genes nearby (H0806H05.13 and H0806H05.14) of unknown function. It constitutes extrons and introns of these two genes.
![]()
<p>Figure 10</p>
Four cases of association between retroelements with O. sativa putative genes observed in this study
Four cases of association between retroelements with O. sativa putative genes observed in this study. Red arrows indicate positions of LTR retroelements with the direction of transcription. Green bars represent NCBI database-predicted gene regions with their orientation of transcription. All the four associations (from F to I) are located in the following genomic clones: AP003627, AC079685, AP003021 and AP003144. (F) An Osr34 solo LTR is associated with two putative genes nearby (P0459B04.22 and P0459B04.23) of unknown function. (G) An Osr8 fragmented LTR is partly associated with a putative LeoOPT1 (oligopeptide transport) gene (OSJNBb0012A20.5). (H) An Osr13 solo LTR constitutes almost half of a putative gene (P0503E05.17) of unknown function; (I) An Osr13 full-length retrotransposon is part of a putative rice gene (P0507H05.22) of unknown function.
![]()