- Research article
- Open Access
Rice transposable elements are characterized by various methylation environments in the genome
BMC Genomics volume 8, Article number: 469 (2007)
Recent studies using high-throughput methods have revealed that transposable elements (TEs) are a comprehensive target for DNA methylation. However, the relationship between TEs and their genomic environment regarding methylation still remains unclear. The rice genome contains representatives of all known TE families with different characteristics of chromosomal distribution, structure, transposition, size, and copy number. Here we studied the DNA methylation state around 12 TEs in nine genomic DNAs from cultivated rice strains and their closely related wild strains.
We employed a transposon display (TD) method to analyze the methylation environments in the genomes. The 12 TE families, consisting of four class I elements, seven class II elements, and one element of a different class, were differentially distributed in the rice chromosomes: some elements were concentrated in the centromeric or pericentromeric regions, but others were located in euchromatic regions. The TD analyses revealed that the TE families were embedded in flanking sequences with different methylation degrees. Each TE had flanking sequences with similar degrees of methylation among the nine rice strains. The class I elements tended to be present in highly methylated regions, while those of the class II elements showed widely varying degrees of methylation. In some TE families, the degrees of methylation were markedly lower than the average methylation state of the genome. In two families, dramatic changes of the methylation state occurred depending on the distance from the TE.
Our results demonstrate that the TE families in the rice genomes can be characterized by the methylation states of their surroundings. The copy number and degree of conservation of the TE family are not likely to be correlated with the degree of methylation. We discuss possible relationships between the methylation state of TEs and their surroundings. This is the first report demonstrating that TEs in the genome are associated with a particular methylation environment that is a feature of a given TE.
Cytosine methylation is not only responsible for regulation of gene expression and immobilization of transposable elements (TEs) [1, 2], but also gives cues about genomic environments . In eukaryotic chromosomes, condensed regions such as the centromere and heterochromatin have intensively accumulated TEs, resulting in heavy DNA methylation in organisms that possess methylation machinery . Recent genome-wide analyses using high-throughput methods revealed that plant genomes show differential methylation states between genes and repetitive sequences: gene-rich regions are present in the hypomethylated fraction, while repetitive sequences are often found in the hypermethylated fraction of the genome [5–7]. A comparative study of methylation status between Arabidopsis ecotypes revealed that TEs are commonly methylated, while genic methylation is highly polymorphic . On the other hand, in mammalian genomes, genes and repeats did not show differential DNA methylation states . Every repetitive sequence has a unique set of features regarding structure, transpositional manner, size, copy number, etc. Tran et al.  reported that the major epigenetic systems (DNA methylation and histone H3 lysine-9) in the Arabidopsis genome target TEs in a genome-wide manner, irrespective of element type and position. However, plant genomes contain an enormous number of repetitive sequences that are widely distributed over chromosomal segments that include hypomethylated euchromatic regions [11, 12]. Little is known about the detailed methylation states of the regions around assorted repetitive sequences in plant genomes. We are interested in the methylation state surrounding the various types of plant repetitive sequences.
According to the decoded sequences , the rice genome comprises 389 mega base pairs (Mb) in total nucleotide length, and is one of the most compact genomes in monocotyledonous species. Despite its small genome size, TEs account for at least 35% of the total genome in rice. Representatives from all known TEs exist in the rice genome. The TEs in the rice genome consist mainly of two types of elements categorized into class I and class II, whose contents account for 20% and 13%, respectively, of the total length of the rice genome. Class I elements, including long terminal-repeat retrotransposons (LTR elements) and non-LTR retrotransposons (longand short-interspersed nucleotide elements: LINEs and SINEs), which transpose in the copy-and-paste manner. In particular, LTR elements are largely localized in gene-poor, heterochromatic regions such as the centromeric and pericentromeric regions . Class II elements, which transpose in the cut-and-paste manner, outnumber class I elements by more than two fold, and often appear in euchromatic regions. The average size of class II elements is less than 300 bp, while that of class I elements is greater than 1000 bp, indicating that a large number of class II elements are scattered in the rice genome . These differentially distributed TEs should reside in various differen genomic environments, particularly in terms of DNA methylation.
Previously, Takata et al.  detected DNA methylation in the flanking regions of two miniature inverted repeat transposable elements (MITEs) in AA-genome Oryza species using a transposon-display (TD) method. The proportion of methylated fragments in the total number of fragments was clearly different in the two MITEs. The methylation frequency estimated for one of the two MITEs was much lower than the average methylation frequency of the entire genomes of Oryza species. The results implied that MITEs show different preferences of the methylation state in their flanking sequences. This led us to pose the question of whether the sequences around TEs are hypermethylated or in specific methylation states.
In order to address the question, here we performed TD analyses targeting 12 rice TEs, which include class I, II, and other-class elements. The degrees of methylation of the TE flanking sequences were investigated by comparing the TD profiles between Msp I and Hpa II-digested samples. The 12 TEs tested here were found to show different features with regard to chromosomal distribution, copy number, and degree of conservation within a cognate family. We attempted to find relationships between the methylation state in the sequences adjacent to the TE and copy number or degree of conservation within the cognate family. The individual TE families reside in different genomic segments with varying degrees of methylation, but each family is present in an environment with a particular degree of methylation.
Results and Discussion
Selection of rice TEs
To examine the methylation states of the sequences flanking various repetitive elements, we selected 12 rice TEs that were classified into three groups: class I, class II, and a third group containing unknown type elements (Table 1). The representative sequences for the individual TEs examined here were those published in the papers listed in Table 1. Of the five elements in class I, noaCRR1 , RIRE5 (Taguchi and Otsubo unpublished result), and RIRE7  belong to the LTR elements, and p-SINE1  is a representative of rice SINE s. Two Ty3/gypsy type elements, noaCRR1 and RIRE7, were distributed around the centromeres. RIRE5, which is a member of the Ty1/copia type elements, was broadly distributed over the rice chromosomes, but not in the centromere. p-SINE1 is small sized (less than 500 bp) and was distributed randomly in euchromatic regions.
The Class II elements tested here consisted of Akan , Kiddo , Kiseru , mPing [20–22], Mashu , TabitoII , and Toya , all of which are miniature inverted repeat transposable elements (MITEs). Although mPing and TabitoII are affiliated with the Tourist family, the two structures are markedly different. Toya may constitute a large group of plant TEs, because elements that resemble Toya have been discovered in the tobacco  and wheat genomes . For Akan, Kiddo, and Mashu, no counterparts have been identified so far. The class II elements in rice are dominated by the remarkable range of MITEs. The sequences for Basho, belonging to the unknown group , are present in certain portions of the rice genome, for instance Turcotte et al.  estimated the proportion of the unknown TE-like sequences in 910 kb of rice genomic DNA sequences to be 1.3%. BLAST searches did not find any specific features of the distribution patterns of the six MITEs and Basho in rice chromosomes, and they appeared to be dispersed over all euchromatic regions (data not shown).
Copy number and conservation of TE sequences
Copy number is an important factor when one considers the cause of homology-dependent DNA methylation. However, it is difficult to estimate the precise numbers of repetitions even if the complete database is available. This is because most of the repetitive sequences have diversified during evolution. The older elements are more differentiated within the cognate sequences than the new elements. In order to survey the copy number and the degree of sequence conservation for the 12 TEs, we performed BLAST searches specifying the rice genome ( as of April, 2007), which takes account of gapped alignments between the query and subject sequences (Table 1). The outcomes of the BLAST searches using the representative sequences published for the 12 TEs as the queries revealed the total hit numbers, alignments of the hit sequences, and the bit scores. The total hit numbers of the 12 TEs were widely distributed from 11 to 835 (Table 1). We estimated the degrees of conservation in the individual TE families indicating the proportion of the conserved copy that was designated over 80% of the maximum bit score in the BLAST search using the most representative sequence, i.e., the 10% conservation degree does not mean 10% identity, but means that 10% of the total copy number in the TE family possesses more than 80% of the maximum bit score. The degree of conservation of the TEs analyzed here varied widely, ranging from 5% to 96%. Low conservation was seen in noaCRR1, RIRE5, and Basho, while the conserved copies of mPing, which is an active element in rice, accounted for 96% of the total copy number (Table 1). The information about the copy number and conservation degree was used for consideration of the relationships between DNA methylation and TEs.
Transposon-displays for the cultivated and wild rice strains
The detection of methylation in TE flanking sequences was performed by the TD method based on amplified fragment length polymorphisms (AFLP), which is basically similar to that described by Takata et al. . We designed the primers taking account of the conserved region within the 60-bp consensus sequences for the 12 TEs (see Additional file 1). All the 60-bp consensus sequences are located in the 5' end of the TEs. The detection of methylated fragments by TD requires a comparison between DNAs digested with methylation-sensitive and insensitive-restriction enzyme isoschizomers, e.g., Hpa II and Msp I, which recognize the same tetranucleotide sequence (5'-CCGG-3'). Hence, we surveyed the presence of cytosine methylation at 5'-CCGG-3' sites, which were located upstream of the TE copies. If a fragment from an Msp I-digested sample was located at the same position as that from the Hpa II-digested sample, the genomic segment was considered to be unmethylated. A fragment was considered to be methylated when the fragment derived from the Msp I-digested sample did not correspond to that from the Hpa II-digested sample. Nine rice genomes from Japonica, Indica, and O. rufipogon were tested in this study, all of which are closely related to each other, and whose hybrids produce fertile seeds. Evaluation of the nine strains was of benefit to verify the efficacy of the TD method for the detection of genome wide methylation. In the analyses, the same gel profile was reproduced in three independent trials of each primer pair.
Because of the ambiguity of fragments larger than 500 bp detected by the TD method, only fragments below 500 bp were considered valid for evaluation in this study. The TD analyses showed different numbers of bands among the 12 TEs (Figure 1), reflecting the copy number of the TE family and its degree of conservation. In all the TEs examined here except mPing, the numbers of fragments below 500 bp detected in the Msp I-digested samples were not significantly different among the nine strains examined (Figure 1). Although the TD for mPing produced 27 - 38 bands below 500 bp in Japonica strains, as expected from the rice genome database, limited numbers (11-4) of bands for mPing were detectable in Indica and O. rufipogon strains. Indeed, it has been reported that Indica strains possess lower copy numbers of mPing than Japonica strains [20, 21].
Detection of DNA methylation
In order to define the average degrees of cytosine methylation in the whole genomes, the AFLP method was employed using two sets of the digested samples (Hpa II/EcoR I and Msp I/EcoR I) for the nine genomes from the cultivated and wild rice strains, which were the same materials as used here. The proportion of methylated fragments in the total number of fragments detected by AFLP ranged between 23% and 28% among the nine genomes (Figure 2). The results implied that an average of 24% of the CCGG motifs in the examined genomes were methylated, as the distribution of Eco RI sites was not strongly influenced by methylation.
Comparison of the TD patterns for the Hpa II and Msp I-digested samples allowed us to estimate the proportions of the methylated fragments flanked by the 5' end of the TEs (see Figure 3). Figure 2 shows the average proportions of methylated fragments in the total MspI fragments of the nine strains. The small standard deviations found in noaCRR1, p-SINE1, and Mashu validate the efficacy of the TD method as capable of detecting methylated fragments equivalently in the different strains. The highest proportions of methylated sequences (over 50%) were observed in the flanking regions of noaCRR1, RIRE5, RIRE7, and Akan. These results were appropriate because the noaCRR1 family and RIRE7 elements are concentrated in the centromeric or pericentromeric regions in rice, where the genomic DNA is heavily methylated [15, 16, 28]. The other LTR element, RIRE5, which is not specifically localized in centromeric or pericentromeric regions, showed a similar tendency to the above LTR elements. These results for the LTR elements agreed with those obtained by high-throughput methods that revealed TEs account for a high proportion of the heavily methylated genomic fraction [5, 6, 9]. Akan is the only class II element that resides in heavily methylated fractions that might be related to heterochromatic regions, but it is localized neither in centromeric nor pericentromeric regions according to BLAST searches (data not shown).
In contrast, the flanking sequences of Mashu (Figure 3) and Basho were hypomethylated (under 20%) in all nine genomes and showed very small standard deviations. Moderately methylated sequences (around 30%) were found in the flanking sequences of Kiddo and Toya. Those of the copies of p-SINE1, KiseruI, mPing, and TabitoII had similar methylation levels to the average methylation state of the whole genome. Thus, individual TEs reside in chromosomal regions with varying levels of methylation, although stable methylation states of the TE flanking sequences were observed within the same TE family. In the class II elements, the proportions of the methylated fragments differed widely, ranging from that of Akan to that of Mashu. Like noaCRR1 and RIRE7, the other individual TE families were present in particular preferential methylation environments of the Oryza genomes.
Relationships with copy number
Mashu has more than 500 copies in the rice genome (Table 1 and Figure 3). Despite the large copy number, Mashu is not an old element, because the copy organization is relatively conserved and the insertion sites are highly polymorphic between Japonica and Indica strains . The hypomethylation environment around the copies of Mashu suggests that the copy number of this TE does not influence the methylation state of the adjacent sequences. The sequence of Mashu itself might be associated with the low methylation state, as are the surrounding sequences. It is not necessarily the case that a TE with an abundant copy number is correlated with a high methylation state. TE activities might not always be controlled by DNA methylation, and instead may be controlled in alternative manner(s) such as by gap repair [30, 31] or unknown host function(s) to prevent harmful behaviors by TEs. Akan, mPing, and Toya have different methylation states in the surrounding sequences among the nine strains, as shown by the large standard deviations (Figure 2). These are low copy number elements (Figure 1), for which it might be difficult to establish stable methylation states around the elements. Besides the low copy number, the instability of the methylation degree around mPing might be due to some activity of the element, which enables it to be inserted into various methylation environments in the genome. Unfortunately, the other active elements in the rice genome were unsuitable for TD analysis due to their very low copy numbers (Tos17  and Karma ) and the presence of a Msp I/Hpa II site at the terminal of the element (nDart [34, 35]). The other class II elements, which should be immobilized or very low active elements, were embedded in relatively stable methylation states of the genomic sequences.
Relationships with conservation degree
The degree of conservation is a rough indication of copy divergence in a TE family (Table 1). Old TE family members tend to have a low degree of conservation and divergent copy organization, while new family members might have better-conserved copy constituents. The conserved elements may be active, like mPing, which was the only active class II element in this study [20–22]. Hence, the degree of conservation is certainly related to the activity of the TE element. Here, we tested whether there is a correlation between TE conservation and methylation around the TE. As shown in Table 1 and Figure 2, we did not find a particular correlation between them. This implies an independence of the TE copy organization from the methylation environment in the genome. However, we do not know whether the methylation environment in the genome affect the activity of the TE elements.
Spread of methylation to the adjacent sequences
When the TD fragments of the Msp I-digested samples were sorted into those with distances from the element below 200 bp and above 200 bp, the degrees of methylation were found to be different in the proximal and distal regions of p-SINE1 and TabitoII (Figure 3 and 4). These proximal regions showed an even heavier methylation state than the distal regions. It seems likely that the methylation state of the TE has spread to the neighboring sequences after insertion. If so, the methylation state of the TE itself should be linked with that of the adjacent sequences. The higher degrees of methylation of sequences within 200-bp of p-SINE1 and TabitoII might reflect the hypermethylation state of the elements (Figure 4). In fact, p-SINE1 copies in the rice genome were strongly suppressed by methylation, whereas p-SINE1 could be transcribed when the methylated state was cancelled by anti-sense genes homologous to Arabidopsis DDM1 or by 5-azacytidine (Ohtsubo et al. unpublished data). In a transgenic rice line retaining a β-glucuronidase reporter gene, the methylation initially occurred in its promoter region derived from the rice pararetrovirus (rice tungro bacilliform virus: RTBV) and gradually expanded to the β-glucuronidase reporter coding sequence in subsequent generations of homozygous offspring . In fact, the rice genome contains about 30 methylated sequences homologous to RTBV , which might induce methylation from the promoter to the reporter coding sequence. In fission yeast, the spread of heterochromatin requires read-through transcription to generate RNA interference (RNAi) via small interfering RNAs (siRNAs) [38, 39]. Similarly, the spread of DNA methylation might require read-through transcription for the primary substrates. As in the case of the LTR elements where the flanking regions are methylated, the methylated sequence could force methylation of the nearby low-methylated sequences. On the other hand, most of the class II elements examined here were present in low-methylated regions. These elements produce few transcripts corresponding to their genomic sequences, because no well-matched transcripts or ESTs were deposited in the databases (data not shown). Therefore, the methylation states of these elements should be as low as their surroundings if methylation is dependent on siRNAs. Alternatively, it is also possible that the methylation states of the surroundings have influenced those of the elements.
TEs characterized by degree of methylation of the surroundings
Target site preferences consisting of only a few base pairs have been identified in certain types of MITEs [40–42]. However, it is still difficult to predict the insertion sequences of TEs in the genome even if the known flanking sequences are available, since these sequences show no consensus. However, as shown by Figure 2, TEs are nicely characterized by the degree of methylation of the TE flanking sequences. Although we do not know yet whether a TE can select an insertion site based on the methylation state, our results clearly show that each TE family is characterized by the methylation state of its surroundings. This is the first report demonstrating that TEs in the genome are associated with a particular methylation environment, which is a feature of a given TE.
TD analysis facilitated a comparison of the methylation degrees among different TEs. In the rice genome database, the 12 TEs analyzed here showed different features of chromosomal distribution, copy number, and conservation degree within a cognate family. The TD analyses revealed that the individual TE families reside in different genomic segments with varying degrees of methylation, but the methylation environments are mostly similar within a given family. The sequences in the vicinity of retrotransposons tend to have a heavier methylation state than those of transposons. A markedly lower methylation state was seen for the flanking sequences of Mashu and Basho. We failed to find a relationship between the copy number of a TE and the methylation state of the adjacent sequences. Our results also detected no correlation between TE conservation and the methylation state around the TE. In the cases of p-SINE1 and TabitoII, the degree of methylation decreased with increasing distance from the elements. If most TEs are methylated, as demonstrated by high-throughput analyses [5–7, 10], the methylation state of a TE itself might not be simply associated with that of its surroundings. Our results demonstrate that rice TEs are localized in genomic environments with different degrees of methylation characteristic of the particular TE.
BLAST and BLAT analyses
In order to investigate the copy numbers and degrees of conservation of the 12 TEs, we used the BLAST website  specified for the rice genome at NCBI. This BLAST analysis considered gaps between the query and subject. All the queries used were representative sequences that were given accession numbers or were selected from reference papers. The other online genome analysis called the Rice Annotation Project DataBase (RAP-DB) developed by the National Institute of Genetics, Japan, can find perfect sequence matches of 33 bases, and sometimes down to 20 bases. The chromosomal locations of the elements can also be searched for by using RAP-BLAT . This system does not consider gaps between the query and subject.
Plant materials and genomic DNA extraction
The nine rice strains used in these experiments were as follows: A58, T65, and Nipponbare (Oryza sativa, Japonica), IR36, Kasalath, and 108 (O. sativa, Indica), W107 and W630 (O. rufipogon, annual), and W593 (O. rufipogon, perennial). These strains are closely related species that can be cross-pollinated with each other. They were maintained as inbred strains by several selfings. In order to obtain genomic DNA with uniform genomic environmental conditions for methylation, all the seeds were sown at the same time and the plants were grown for 3 weeks at about 25°C. Total genomic DNA was isolated from the above-ground parts of the plants (excluding roots), as described by Nagano et al. .
Amplified fragment polymorphism (AFLP) and transposon display (TD) to analyze methylation sites
In order to investigate the average frequencies of the methylation-sensitive sites in the rice and related genomes, we used an AFLP method that was basically similar to that reported by Takata et al. . The isoschizomers Hpa II and Msp I were used instead of Mse I and Taq I as 'frequent cutter' enzymes. Both Hpa II and Msp I recognize the same tetranucleotide sequence (5'-CCGG-3'), but exhibit differential sensitivity to DNA methylation. Hpa II is sensitive to the methylation of any cytosine of 5'-CCGG-3', whereas Msp I cuts 5'-C5mCGG-3', but not 5'-5mCCGG-3'. Thus, 5'-5mCCGG-3' or 5'-5mC5mCGG-3' at the proximal site of TEs may not give band in this TD method. Eco RI (5'-GAATTC-3') was employed as the 'less-frequent cutter' enzyme. The oligonucleotides used for adaptors and primers in this study are listed in Additional file 2. All the adaptors were prepared by annealing the two adaptor primers at room temperature after boiling. DNA digestion, adaptor ligation, and PCR reactions were performed according to the method of Takagi et al. . Two sets of digestion reactions (Hpa II/Eco RI and Msp I/Eco RI) were carried out simultaneously. Approximately 0.25 μg of rice genomic DNA was digested with 5 units of each enzyme in 25 μl of a suitable digestion buffer. Each of the Hpa II/Msp I and Eco RI adaptors (50 μM) was ligated to the digested genomic DNA (200 ng) by the addition of 10 μl of a mix containing 1 × ligation buffer, 0.5 μM ATP, and 1 Weiss unit of T4 DNA ligase and incubation at 16°C overnight. The first PCR reaction was performed in 20 μl of a reaction mixture containing 1 × PCR buffer, 0.2 mM dNTPs (Takara Shuzo Co Ltd., Osaka, Japan), 15 ng of the adaptor-ligated DNA, 10 pmol of each first primer, and 1 unit of EX Taq polymerase (Takara Shuzo Co Ltd.). The first PCR reaction was carried out as follows: 72°C/2 min; 94°C/3 min; 24 cycles of 94°C/30 sec, 54°C/30 sec and 72°C/1 min, and a final cycle of 72°C/5 min. The PCR product was then diluted 20-fold with sterilized water. Nested primers were used for the second PCR: the 5' end of the second primer was labeled with FITC (lambda em = 495 nm), and the 3' end of the second adaptor primers was overhung by three selected nucleotides (ATC, CTC, GAT, GTT or TTA). These five different types of 3' ends of the second adaptor primers were used for each trial (see Additional file 2). The second PCR reaction was performed with the same reaction mixture as the first PCR plus 2 μl of 50% DMSO. For temperature cycling, we modified the conditions used in the "touchdown method": 94°C/5 min, then 10 cycles beginning with 94°C/30 sec, 68°C/30 sec, 72°C/1 min, followed by a reduction of the annealing temperature of the first cycle (68°C) with a 1°C increment per cycle, and then 27 cycles at an annealing temperature of 58°C, and a final incubation at 72°C/5 min. The reaction products were electrophoresed on a 3.5% denaturing polyacrylamide gel containing 50 ml of a mixture of 3.5 ml of Long Ranger (Takara Shuzo Co Ltd.) (equivalent to 50% polyacrylamide:bisacrylamide solution = 19:1) and 7.0 M urea. Internal Lane Standard 600 (Promega, Madison, Wisconsin, USA) was used as the size marker from 60 bp to 600 bp, which was labeled with carboxy-X-rhodamine (lambda em = 576 nm). After the gel was pre-run using 1 × TBE as running buffer, the samples were electrophoresed at a constant 1700 V for 80 min. The fragments and size markers were visualized by scanning the fluorescent signals in the gel plate using a Typhoon 8600 scanner (Molecular Dynamics Inc. Sunnyvale, California).
TD analyses were applied for detection of methylation in the flanking sequences of the 12 TEs. The TD basically followed the AFLP method described above with the following modifications. In order to design the primers in the 5' side of the individual TE sequences, we made consensuses for 60 bp of the 5' end sequences of the individual TEs after surveying sequences homologous to the representative TEs using RAP-BLAT  (see Additional file 1). Based on the consensus for each element, nested primers for the TEs were designed, and the 5' end of the second TE primer was labeled with FITC (see Additional file 2). No overhung nucleotide was added at the 3' end of the second adaptor primers in the TD method. In the cases of LTR elements, the PCR reactions elongate to both the outside and inside of the elements, because the primer sequences reside in the tandem repeats of the LTR elements. A very limited number of the TD bands were directed to the internal sequence of LTR elements due to its higher conservation among the copies compared with the flanking sequences. Thus, we did not take these bands away from the total band numbers when the proportions of methylated bands were calculated.
amplified fragment length polymorphism
miniature inverted repeat transposable element
long interspersed nucleotide element
- LTR element:
long terminal-repeat retrotransposon
short interspersed nucleotide element
Yoder JA, Walsh CP, Bestor TH: Cytosine methylation and the ecology of intragenomic parasites. Trends Genet. 1997, 13 (8): 335-340. 10.1016/S0168-9525(97)01181-5.
Bender J: DNA methylation and epigenetics. Annu Rev Plant Biol. 2004, 55: 41-68. 10.1146/annurev.arplant.55.031903.141641.
Rangwala SH, Richards EJ: The value-added genome: building and maintaining genomic cytosine methylation landscapes. Curr Opin Genet Dev. 2004, 14 (6): 686-691. 10.1016/j.gde.2004.09.009.
Wong LH, Choo KH: Evolutionary dynamics of transposable elements at the centromere. Trends Genet. 2004, 20 (12): 611-616. 10.1016/j.tig.2004.09.011.
Zhang X, Yazaki J, Sundaresan A, Cokus S, Chan SW, Chen H, Henderson IR, Shinn P, Pellegrini M, Jacobsen SE, Ecker JR: Genome-wide high-resolution mapping and functional analysis of DNA methylation in arabidopsis. Cell. 2006, 126 (6): 1189-1201. 10.1016/j.cell.2006.08.003.
Zilberman D, Gehring M, Tran RK, Ballinger T, Henikoff S: Genome-wide analysis of Arabidopsis thaliana DNA methylation uncovers an interdependence between methylation and transcription. Nat Genet. 2007, 39 (1): 61-69. 10.1038/ng1929.
Rabinowicz PD, Citek R, Budiman MA, Nunberg A, Bedell JA, Lakey N, O'Shaughnessy AL, Nascimento LU, McCombie WR, Martienssen RA: Differential methylation of genes and repeats in land plants. Genome Res. 2005, 15 (10): 1431-1440. 10.1101/gr.4100405.
Vaughn MW, Tanurd Ic M, Lippman Z, Jiang H, Carrasquillo R, Rabinowicz PD, Dedhia N, McCombie WR, Agier N, Bulski A, Colot V, Doerge RW, Martienssen RA: Epigenetic Natural Variation in Arabidopsis thaliana. PLoS Biol. 2007, 5 (7): e174-10.1371/journal.pbio.0050174.
Rabinowicz PD, Palmer LE, May BP, Hemann MT, Lowe SW, McCombie WR, Martienssen RA: Genes and transposons are differentially methylated in plants, but not in mammals. Genome Res. 2003, 13 (12): 2658-2664. 10.1101/gr.1784803.
Tran RK, Zilberman D, de Bustos C, Ditt RF, Henikoff JG, Lindroth AM, Delrow J, Boyle T, Kwong S, Bryson TD, Jacobsen SE, Henikoff S: Chromatin and siRNA pathways cooperate to maintain DNA methylation of small transposable elements in Arabidopsis. Genome Biol. 2005, 6 (11): R90-10.1186/gb-2005-6-11-r90.
Mao L, Wood TC, Yu Y, Budiman MA, Tomkins J, Woo S, Sasinowski M, Presting G, Frisch D, Goff S, Dean RA, Wing RA: Rice transposable elements: a survey of 73,000 sequence-tagged-connectors. Genome Res. 2000, 10 (7): 982-990. 10.1101/gr.10.7.982.
Arabidopsis_Genome_Initiative: Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature. 2000, 408 (6814): 796-815. 10.1038/35048692.
International_Rice_Genome_Sequencing_Project: The map-based sequence of the rice genome. Nature. 2005, 436 (7052): 793-800. 10.1038/nature03895.
Takata M, Kishima Y, Sano Y: DNA methylation polymorphisms in rice and wild rice strains: Detection of epigenetic markers. Breeding Science. 2005, 55 (1): 57-63. 10.1270/jsbbs.55.57.
Nagaki K, Neumann P, Zhang D, Ouyang S, Buell CR, Cheng Z, Jiang J: Structure, divergence, and distribution of the CRR centromeric retrotransposon family in rice. Mol Biol Evol. 2005, 22 (4): 845-855. 10.1093/molbev/msi069.
Kumekawa N, Ohmido N, Fukui K, Ohtsubo E, Ohtsubo H: A new gypsy-type retrotransposon, RIRE7: preferential insertion into the tandem repeat sequence TrsD in pericentromeric heterochromatin regions of rice chromosomes. Mol Genet Genomics. 2001, 265 (3): 480-488. 10.1007/s004380000436.
Mochizuki K, Umeda M, Ohtsubo H, Ohtsubo E: Characterization of a plant SINE, p-SINE1, in rice genomes. Jpn J Genet. 1992, 67 (2): 155-166. 10.1266/jjg.67.155.
Nagano H, Kunii M, Azuma T, Kishima Y, Sano Y: Characterization of the repetitive sequences in a 200-kb region around the rice waxy locus: diversity of transposable elements and presence of veiled repetitive sequences. Genes Genet Syst. 2002, 77 (2): 69-79. 10.1266/ggs.77.69.
Yang G, Dong J, Chandrasekharan MB, Hall TC: Kiddo, a new transposable element family closely associated with rice genes. Mol Genet Genomics. 2001, 266 (3): 417-424. 10.1007/s004380100530.
Jiang N, Bao Z, Zhang X, Hirochika H, Eddy SR, McCouch SR, Wessler SR: An active DNA transposon family in rice. Nature. 2003, 421 (6919): 163-167. 10.1038/nature01214.
Kikuchi K, Terauchi K, Wada M, Hirano HY: The plant MITE mPing is mobilized in anther culture. Nature. 2003, 421 (6919): 167-170. 10.1038/nature01218.
Nakazaki T, Okumoto Y, Horibata A, Yamahira S, Teraishi M, Nishida H, Inoue H, Tanisaka T: Mobilization of a transposon in the rice genome. Nature. 2003, 421 (6919): 170-172. 10.1038/nature01219.
Schenke D, Sasabe M, Toyoda K, Inagaki YS, Shiraishi T, Ichinose Y: Genomic structure of the NtPDR1 gene, harboring the two miniature inverted-repeat transposable elements, NtToya1 and NtStowaway101. Genes Genet Syst. 2003, 78 (6): 409-418. 10.1266/ggs.78.409.
Saito M, Yonemaru J, Ishikawa G, Nakamura T: A candidate autonomous version of the wheat MITE Hikkoshi is present in the rice genome. Mol Genet Genomics. 2005, 273 (5): 404-414. 10.1007/s00438-005-1144-7.
Inukai T, Sano Y: Sequence rearrangement in the AT-rich minisatellite of the novel rice transposable element Basho. Genome. 2002, 45 (3): 493-502. 10.1139/g02-010.
Turcotte K, Srinivasan S, Bureau T: Survey of transposable elements from rice genomic sequences. Plant J. 2001, 25 (2): 169-179. 10.1046/j.1365-313x.2001.00945.x.
BLAST Plant genomes. [http://www.ncbi.nlm.nih.gov/BLAST/Genome/PlantBlast.shtml?10]
Mizuno H, Ito K, Wu J, Tanaka T, Kanamori H, Katayose Y, Sasaki T, Matsumoto T: Identification and mapping of expressed genes, simple sequence repeats and transposable elements in centromeric regions of rice chromosomes. DNA Res. 2006, 13 (6): 267-274. 10.1093/dnares/dsm001.
Takagi K, Nagano H, Kishima Y, Sano Y: MITE-transposon display efficiently detects polymorphisms among the Oryza AA-genome species. Breeding Science. 2003, 53 (2): 125-132. 10.1270/jsbbs.53.125.
Rubin E, Levy AA: Abortive gap repair: underlying mechanism for Ds element formation. Mol Cell Biol. 1997, 17 (11): 6294-6302.
Yamashita S, Takano-Shimizu T, Kitamura K, Mikami T, Kishima Y: Resistance to gap repair of the transposon Tam3 in Antirrhinum majus: a role of the end regions. Genetics. 1999, 153 (4): 1899-1908.
Hirochika H, Sugimoto K, Otsuki Y, Tsugawa H, Kanda M: Retrotransposons of rice involved in mutations induced by tissue culture. Proc Natl Acad Sci USA. 1996, 93 (15): 7783-7788. 10.1073/pnas.93.15.7783.
Komatsu M, Shimamoto K, Kyozuka J: Two-step regulation and continuous retrotransposition of the rice LINE-type retrotransposon Karma. Plant Cell. 2003, 15 (8): 1934-1944. 10.1105/tpc.011809.
Tsugane K, Maekawa M, Takagi K, Takahara H, Qian Q, Eun CH, Iida S: An active DNA transposon nDart causing leaf variegation and mutable dwarfism and its related elements in rice. Plant J. 2006, 45 (1): 46-57. 10.1111/j.1365-313X.2005.02600.x.
Fujino K, Sekiguchi H, Kiguchi T: Identification of an active transposon in intact rice plants. Mol Genet Genomics. 2005, 273 (2): 150-157. 10.1007/s00438-005-1131-z.
Kloti A, He X, Potrykus I, Hohn T, Futterer J: Tissue-specific silencing of a transgene in rice. Proc Natl Acad Sci USA. 2002, 99 (16): 10881-10886. 10.1073/pnas.152330299.
Kunii M, Kanda M, Nagano H, Uyeda I, Kishima Y, Sano Y: Reconstruction of putative DNA virus from endogenous rice tungro bacilliform virus-like sequences in the rice genome: implications for integration and evolution. BMC Genomics. 2004, 5 (1): 80-10.1186/1471-2164-5-80.
Irvine DV, Zaratiegui M, Tolia NH, Goto DB, Chitwood DH, Vaughn MW, Joshua-Tor L, Martienssen RA: Argonaute slicing is required for heterochromatic silencing and spreading. Science. 2006, 313 (5790): 1134-1137. 10.1126/science.1128813.
Volpe TA, Kidner C, Hall IM, Teng G, Grewal SI, Martienssen RA: Regulation of heterochromatic silencing and histone H3 lysine-9 methylation by RNAi. Science. 2002, 297 (5588): 1833-1837. 10.1126/science.1074973.
Bureau TE, Wessler SR: Stowaway: a new family of inverted repeat elements associated with the genes of both monocotyledonous and dicotyledonous plants. Plant Cell. 1994, 6 (6): 907-916. 10.1105/tpc.6.6.907.
Bureau TE, Ronald PC, Wessler SR: A computer-based systematic survey reveals the predominance of small inverted-repeat elements in wild-type rice genes. Proc Natl Acad Sci USA. 1996, 93 (16): 8524-8529. 10.1073/pnas.93.16.8524.
Naito K, Cho E, Yang G, Campbell MA, Yano K, Okumoto Y, Tanisaka T, Wessler SR: Dramatic amplification of a rice transposable element during recent domestication. Proc Natl Acad Sci USA. 2006, 103 (47): 17620-17625. 10.1073/pnas.0605421103.
RAP-DB: rice annotation project database. [http://rapdb.dna.affrc.go.jp/tools/blat]
We thank Hironori Nagano, Hokkaido University, for technical suggestions and comments on the manuscript and Yu Taguchi, University of Tokyo, for providing unpublished data. This work was supported by grants from the Ministry of Education, Science, and Culture, Japan.
MK, AK, and AT conceived the experiments and carried out the molecular analyses and databases searches. YK conceived of the study and drafted the manuscript. HO participated in the sequence alignment, database analysis, and the design of the study. YK and YS participated in the basic design of this study and coordination. All authors read and approved the final manuscript.
Miwako Takata, Akihiro Kiyohara, Atsuko Takasu contributed equally to this work.
Electronic supplementary material
Additional file 1: The 60-bp consensus sequences of the 5' end of 12 TEs. The data provided represent the 60-bp consensus sequences determined from the rice database for the 5' end of 12 TEs. (PDF 31 KB)
About this article
Cite this article
Takata, M., Kiyohara, A., Takasu, A. et al. Rice transposable elements are characterized by various methylation environments in the genome. BMC Genomics 8, 469 (2007). https://doi.org/10.1186/1471-2164-8-469
- Amplify Fragment Length Polymorphism
- Methylation State
- Rice Genome
- Pericentromeric Region
- Amplify Fragment Length Polymorphism