Skip to main content

Evolution of ribosomal DNA-derived satellite repeat in tomato genome

Abstract

Background

Tandemly repeated DNA, also called as satellite DNA, is a common feature of eukaryotic genomes. Satellite repeats can expand and contract dramatically, which may cause genome size variation among genetically-related species. However, the origin and expansion mechanism are not clear yet and needed to be elucidated.

Results

FISH analysis revealed that the satellite repeat showing homology with intergenic spacer (IGS) of rDNA present in the tomato genome. By comparing the sequences representing distinct stages in the divergence of rDNA repeat with those of canonical rDNA arrays, the molecular mechanism of the evolution of satellite repeat is described. Comprehensive sequence analysis and phylogenetic analysis demonstrated that a long terminal repeat retrotransposon was interrupted into each copy of the 18S rDNA and polymerized by recombination rather than transposition via an RNA intermediate. The repeat was expanded through doubling the number of IGS into the 25S rRNA gene, and also greatly increasing the copy number of type I subrepeat in the IGS of 25-18S rDNA by segmental duplication. Homogenization to a single type of subrepeat in the satellite repeat was achieved as the result of amplifying copy number of the type I subrepeat but eliminating neighboring sequences including the type II subrepeat and rRNA coding sequence from the array. FISH analysis revealed that the satellite repeats are commonly present in closely-related Solanum species, but vary in their distribution and abundance among species.

Conclusion

These results represent that the dynamic satellite repeats were originated from intergenic spacer of rDNA unit in the tomato genome. This result could serve as an example towards understanding the initiation and the expansion of the satellite repeats in complex eukaryotic genome.

Background

The large variety of genome sizes found throughout the plant kingdom is mainly attributed to species-specific differences in ploidy and repetitive DNA content [1]. Repetitive DNA can be divided into two categories: interspersed repeats, which are individual repeat units that are distributed around the genome in an apparently random fashion, and tandem repeated DNA, whose repeat units are placed next to each other in an array. Several previous studies have uncovered interspersed repeats, the retrotransposons, which are usually the most abundant form of repetitive DNA in plants with large genomes [2, 3]. The transposition mechanism of these repeats has been well characterized with respect to interspersed repeats. Multigene families including ribosomal RNAs (rRNA) as well as noncoding sequences such as satellite DNA, minisatellite sequences and microsatellite sequences are often arranged in tandem arrays [4, 5]. Tandemly repeated DNA is primarily found at centromeres, subtelomeric regions, and heterochromatin. Recently, a number of new satellite repeats have been described in higher plants using cytological techniques [6–8].

Ribosomal DNA is one of the most well-characterized tandem arrays and is made up of genes that are transcribed into the components of the ribosome [9]. The repeated unit consists of the 18S, 5.8S, and 25S rRNA genes, external transcribed spacers, internal transcribed spacers, and an intergenic spacer (IGS). The coding regions of rDNA are highly conserved among eukaryotic organisms, whereas the sequence of the noncoding IGS region varies broadly between even closely-related species. This observation has been explained by the model of horizontal or concerted evolution, originally proposed by Brown et al. [10]. The identity of coding sequences from different species can be explained to have occurred through the maintenance of sequences with strong purifying selection. However, in plants, such as legumes [11], potato [12], and tobacco [7, 13], highly amplified satellite repeats with sequence homologous to the IGS subrepeats of rDNA have been reported to exist in dispersed patterns over several chromosomes. In these genomes, the IGS subrepeat-homologous satellite sequences occur in blocks independent of the rRNA gene cluster [12, 14, 15]. The satellite sequences described in plants often have erratic distributions and large differences in abundance, between even closely related species [7, 15]. The discovery of satellite repeats homologous to the IGS of 45S rDNA induced speculation that satellite repeats might be originated from 45S rDNA. The mechanism of satellite repeat generation has been explained by several hypotheses. (1) Satellite repeats could have arisen through repeated and random unequal crossing over [16], (2) by replication slippage and unequal crossing over with subsequent expansion [17], and (3) by the products of rolling circle replication of extrachromosomal circular DNAs that became re-inserted into the genome [18, 19]. The segmental duplication of large arrays of satellite repeats has also been proposed to be the primary mechanism responsible for their amplification, contributing to the rapid reshuffling of CentO satellites in rice centromere [20, 21]. However, to date, there is no clear explanation on how the sequence of the rDNA unit escaped from the highly efficient concerted evolutionary mechanisms that keep it so well conserved. To better understand the origins of satellite repeats, it is necessary to find and compare sequences from genomes in different stages along the path of satellite repeat generation [22]. Here we report the discovery of a satellite repeat that is highly homologous to IGS of the 18S–25S rRNA genes in tomato. The comparison of sequences from several BAC clones containing rDNA in various stages of modification has provided a plausible explanation for how IGS homologous satellite repeats were developed from the well-conserved rDNA unit.

Results

Cytological localization of 45S rDNA and IGS-homologous repeats in the tomato genome

As part of the international tomato genome sequencing project, we employed fluorescence in situ hybridization (FISH) analysis for confirmation of the position of genetic marker-anchored BAC clones on chromosome 2 [23, 24]. LE-HBa0007F24, a clone anchored at genetic marker cLER-1-H17, produced very strong signals in the nucleolus organizing region (NOR) of chromosome 2 and on three other chromosomes (Figure 1A). However, hybridization with wheat 45S rDNA (pTa71, GeneBank accession number: X07841) produced only one signal focused on the short arm of tomato chromosome 2, indicating that this is the only location of rDNA repeats (Figure 1B). When the nucleotide sequences of tomato and wheat 45S rDNAs were compared, the 18S rRNA genes were 96% identical, whereas the IGS sequences had very low identity [25]. Therefore, we speculated that only the short arm of chromosome 2 contains the canonical 45S rDNA unit and that the foci on the other three chromosomes do not contain the coding sequences of 45S rDNA, but some IGS sequences.

Figure 1
figure 1

The distribution of IGS-homologous satellite repeats. (A) FISH signals (red) obtained with probe derived from Hba0007F24, containing tomato rDNA (red arrow head). Bar, 10 μm (B) FISH signal (red) obtained from heterologous rDNA probe, pTa71, for wheat 25-18S rDNA (red arrow head). Bar, 10 μm (C) FISH signals (green) obtained with probe, pIGS, for type I subrepeat of rDNA (arrow head). Bar, 10 μm (D) FISH signal on DNA fibers prepared from S. lycopersicum with pTa71 (green) and pIGS (red) probes. IGS homologous satellite repeat (arrow) and rDNA array were detected (arrow head).

In order to test our hypothesis, we performed FISH analysis with a tomato-specific IGS probe (pIGS) made from 483 bp of tomato sequence amplified from the type I subfamily IGS of 25-18S rDNA (Figure 1C). Like the tomato 45S rDNA probe and unlike the wheat 45S rDNA probe, the tomato type I IGS probe hybridized to loci on four chromosomes. FISH analysis on extended DNA fibers prepared from S. lycopersicum confirmed that there are two types of IGS organization in the tomato genome (Figure 1D). One is co-localized with the coding sequence of rRNA genes and the other is linearly stretched over a 300 kb region that lacks rRNA genes. Because FISH analysis had demonstrated that the HBa0007F24 clone is derived from the same chromosomes that hybridize to the tomato 45S rDNA probe, a probe was made from the partial 18S rDNA fragment of HBa0007F24 and tested in another FISH experiment (data not shown). This probe hybridized only to the NOR on the short arm of tomato chromosome 2. These data indicate that there is a single known canonical 45S rDNA block and three IGS homologous satellite repeats that are independent of rRNA genes in the S. lycopersicum genome.

Sequence analyses reveal that transition of rDNA is initiated in the NOR

To determine how the satellite repeats are generated from the 45S rDNA array, we determined full sequences of two BAC clones in distinct stage of divergence of rDNA repeat: HBa0007F24 (131,560 bp, GeneBank accession number: AC215351) of which end sequence was 89% identical to tomato 45S rDNA (GeneBank accession numbers: AY366528, AY366529) and Sle0089P21 (17,800 bp, GenebBank accession number: AC215459) which contains two copies of canonical rDNA (Figure 2A–B). Sequence comparison of two BAC sequences revealed dynamic changes in the HBa0007F24 sequence which consists of seven diverged truncated rDNA units with similar composition but different lengths of elements. Essentially every unit has 18S, 25S, 5.8S, and IGS. The modified 18S and 25S rDNAs showed 92–94% and 89–91% identity to typical rDNAs, respectively. Annotation of the sequence revealed that there are three significant modifications of the 45S rDNA in the HBa0007F24 clone (Figure 2). First, the 18S rDNAs were interrupted by LTR-type retrotransposons. Second, the 25S rDNAs were fragmented by the IGS-like sequences. Third, IGS sequences which have 3 to 4 times longer than normal rDNA unit were found between the 25-18S rDNAs.

Figure 2
figure 2

rDNA variants in S. lycopersicum genome. (A) Typical rDNA units in Sle0089P21 (18,122 bp). (B) Modified rDNA units found in BAC clone, HBa0007F24 (131,560 bp). There are seven variant repeats in the BAC clone. Retrotransposon sequences and extra-IGS regions are detected in 18S rRNA and 25S rRNA genes, respectively. Roman numerals in parentheses are the numbers of rDNA units in HBa0007F24. (C) More expanded IGS sequence in Sle0049A24 (28,040 bp).

To further characterize the IGS homologous repeat sequences in the tomato genome, we selected and sequenced the BAC clone, Sle0049A24 (28,040 bp), which has expanded IGS sequences at both ends, according to the results of a BLAST search of the BAC end sequence database (Figure 2C). The majority (85%, 23,858 bp) of the total 28,040 bp sequence of Sle0049A24 (GeneBank accession number: AC225927) is a long IGS stretch, which are 86% identical to subrepeat I. The remainder of the sequence (1,753 bp) is partial 25S rDNA and subrepeat II of IGS. The length of the repeated unit was well-conserved as 52–53 bp-long throughout the 21,607 bp of continuous subrepeat type I sequence.

LTR retrotransposons inserted in the 18S rDNAs were polymerized by uneven recombination

Sequence comparison of HBa0007F24 and Sle0089P21 showed that the tomato rDNA related retrotransposons (TRRTs) were inserted at the same location of the seven 18S rDNAs in HBa0007F24 (Figure 3). Following the accepted system for retrotransposon nomenclature [26], it was classified as a Ty3-Gypsy like LTR retrotransposon. TRRT2 through TRRT7 are orientated in the direction of rDNA transcription, and TRRT1 has the opposite orientation (Figure 3C). To determine the reason for the opposite orientation of TRRT1, we analyzed the flanking sequences of each TRRT (Figure 3D). We found that each TRRT has the same flanking sequences (CTAC), indicating that TRRTl or others were inverted by recombination after insertion rather than inserted in the opposite orientation. Therefore, we assume that this inversion was mediated by rearrangement among duplicated segments.

Figure 3
figure 3

Distribution of TRRT in the tomato genome. (A) Chromosomal localization of TRRT (red signal). FISH analysis of pachytene stage S. lycopersicum was probed with TRRT. TRRT was localized on the heterochromatin of all chromosome include NOR of chromosome 2 (arrow). Bar, 10 μm (B) Fiber FISH shows that the TRRTs (red signal) are localized in the rDNA array (pTA71, green signal). (C) Arrows indicate the direction of seven TRRTs present in the HBa007F24 clone. TRRT in unit IV has solo LTR sequence. Blue box (18S rDNA) and green box (25S rDNA) indicate fragmented rRNA genes. (D) Comparison of flanking sequences of TRRTs with those of typical 18S rRNA genes. Seven TRRTs are inserted in the same site and create 4 bp TDS (shadow boxes). (E) Organization of TRRT.

Phylogenetic analysis of 13 LTRs belonging to the seven TRRTs of HBa0007F24 demonstrated that segmental duplications were a major process for TRRT amplification (Figure 4C, see Additional file 1). The results from computing the proportion of nucleotide differences between each pair of LTR sequences showed that no LTR pair of a single LTR retrotransposon was clustered together implying that the retrotransposons were duplicated rather than re-inserted via a intermediate RNA. Two LTRs of TRRT1 were the most closely clustered, whereas LTRs of TRRT5 and 6 had the most diverged sequences between any two LTRs of a single LTR-retrotransposon. However, the same positions in TRRT5 and TRRT6 (RT5-5':RT6-5', RT5-3':RT6-3') were closely related, indicating they were the most recently duplicated by recombination. Phylogenetic analysis of 18S and 25S rDNAs of HBa0007F24 with typical rRNA genes shows similar results obtained in the analysis of the LTRs (Figure 4A, B).

Figure 4
figure 4

Phylogenetic analysis of variants in HBa0007F24 sequence. Neighbor-Joining tree obtained for 18S rDNA (A), 25S rDNA (B) and LTRs of TRRT(C). Fragmented 18S rDNA and 25S rDNA in HBa0007F24 put together deleting the inserts (18S joins, 25S joins). Comparison of phylogenetic distance of 13 LTR pairs of 7 retrotransposons shows that LTR pair of each retrotransposon is not clustered together implying that the retrotransposons was duplicated rather than transposition for LTR retrotransposon amplification. Opposite position LTRs (RT5 LTR 3': RT5 LTR 5' or RT6 LTR 3': RT6 LTR 5') of TRRT5 or TRRT6 have the most divergent sequence, but the same position LTR (RT5 LTR 3': RT6 LTR 3' or RT6 LTR 5': RT5 LTR 5') of TRRT5 and TRRT6 are closely clustered.

The retrotransposon encodes four proteins, retrotrans gag, reverstranscriptase, integrase core domain, chromatin organizing modifier, and has long terminal repeats (LTR) at both ends (Figure 3E). However, TRRT4 contains solo LTR, but does not encode the gag-pol gene, indicating that some of the sequence was lost through unequal recombination [2]. The retrotransposon insertions vary in length, from 5,645 to 6,028 bp, and share 93.61% sequence identity with each other. The 18S rDNAs associated with the TRRTs have a similar degree of sequence identity to the canonical 18S rRNA gene (92–94% identity).

Because all of the retrotransposons were found at the same position, 781 bp of the 18S rDNA, with the same flanking sequences, we examined whether the TRRT transposed site-specific manner. Using the retrotransposon sequence as a query to search GeneBank (BLASTN), we identified two tomato BAC clones, C02HBa0155E05 and C06HBa0169D11, which have the same retrotransposon, but are not associated with rDNA. The flanking sequences of these retrotransposons were different from those of HBa0007F24. FISH analysis also demonstrated that the retrotransposon sequence is present on other chromosomes as various sizes of blocks supporting that the retrotransposon was not integrated in a site-specific manner (Figure 3A).

Duplication of the IGS between 25-18S rDNA sequences into the middle of 25S rRNA

Sequence comparison of HBa0007F24 and Sle0089P21 showed that IGS-like sequences are inserted at the same position of 25S rDNA, between 1,388 bp and 1,641 bp of the canonical sequence, and each interrupted copy has lost 254 bp of 25S rDNA sequence where the IGS was inserted (Figure 5A and 8). Of seven 25S rDNAs, six had IGS-like sequences in the middle of the coding region doubling the number of IGS in the 45S rDNA unit. The remaining 25S rDNA (Figure 5A-II) was missing part of its 5' end.

Figure 5
figure 5

Duplicated IGS in 25S rDNA. (A) Deletion patterns of the 25S rDNA in the HBa007F24 clone. All IGS-insertion occurred at the same site, 1,388–1,641 bp of 25S rDNA sequence, losing 253 bp. Assuming that deletions were accumulated, the sequence of deletion events can be estimated to be IV→ V→ III from the deletion pattern. Numbers on shade box represent position of normal 25S rRNA gene. Position and length of deleted sites are indicated. (B) IGSs duplicated in 25S rDNA were compared with normal IGS of 25-18S rDNA. Two subrepeats were differently regulated: Type II subrepeats have been easily lost, but type I subrepeats were not.

The size and the number of deletion of each unit of 25S rDNA demonstrated that the deletions were accumulated (Fig 5A). For example, all units (II, III, IV, V, and VII) of 25S rDNAs in the BAC clone have a deletion of 81–82 bp at the same position 2,986–3,066 bp and the deletion frequency varies from 2 to 5. Unit IV has a single 82 bp deletion, while unit III has four such deletions. Following the number of deletions in each unit, the order of recombination can be deduced as: IV → V → III. These results indicate that 25S rDNAs harboring IGS were mainly multiplied by unequal recombination from a single variant rather than individual IGSs being recombined in parallel in their respective 25S rDNA sequences.

The duplicated IGS-like sequences in the 25S rDNA commonly contain conserved sub-family repeat sequences, namely type I subrepeat, AT rich regions, and type II subrepeat (Figure 5B). The length of IGS-like sequences inserted into 25S rDNAs varies from 1,476 bp to 2,074 bp and primarily depends on the length of the type II subrepeat, but not type I subrepeat. These data indicate that type I and type II subrepeats have been differentially regulated during molecular evolution.

Differential amplification between subfamily repeats in the rDNA intergenic spacer

Sequence comparison of HBa0007F24 and Sle0089P21 showed that the length of the IGSs of 25-18S rDNAs (8,400 – 11,408 bp) in HBa0007F24 was 3 to 4 times longer than the normal IGS sequences (3,395 bp) present in tomato genome (Figure 6A, B). In addition, as shown in Figure 6, the type II subrepeat present at the downstream of the transcription initiation site (TIS) has been replaced by the type I subrepeat. In the amplified IGS, the type I subrepeat in the upstream of the TIS was 10–17 times longer (4,771–7,921 bp) than conventional type I subrepeat (448 bp). Furthermore, another type I subrepeat found downstream of the TIS was 2–3 times longer (1,040–1,478 bp) than the conventional type I subrepeat (448 bp). Therefore, these results strongly indicated that duplication of the IGS into the 25S rDNA occurred before the type II subrepeat replacement by the type I subrepeat, and also before the type I subrepeat amplification of the IGS located between the 25S and 18S rDNAs. Even though the length of the type I subrepeat was expanded, the monomer length was well conserved as 53 bp.

Figure 6
figure 6

Comparison of typical IGS and modified IGS. (A) Dot blot analysis of normal IGS from Sle0089P21 and amplified IGS from unit IV of HBa0007F24. Red dot represents complimentary match between two sequences. Dot blot parameter: window = 15, mismatch = 0. (B) Diagram representing the structures of normal IGS and amplified IGS of 25-18S rDNA. Type II subrepeat in the modified IGS is replaced with type I subrepeat. Left-side box represents 25S rDNA and right-side box represents 18S rDNA. (C) Fragmental duplication revealed by the Neighbor-joining tree among 129 monomers of IGS unit IV. The closely related pairs of monomers are connected.

To figure out the molecular mechanisms of the repeat proliferation, we performed phylogenetic analysis of 129 repeat monomers that were identified in expanded IGS in the unit IV of HBa0007F24 sequence (Figure 6C). By analyzing the most related monomers revealed by the Neighbor-Joining tree obtained (see Additional file 2), we identified 25 pairs of monomers that are arranged in four duplicated clusters of monomers. The data obtained demonstrated that segmental duplication was occurred among the amplified type I subrepeat.

Distribution of IGS-homologous repeats among closely-related tomato species

We carried out FISH analysis to study the organization of the 45S rDNA locus and IGS-homologous repeats on the Eulycopersicon red fruited subgenera [27], including S. lycopersicum, S. lycopersicum var.cerasiforme, and S. pimpinellifolium, which are very closely-related species (Figure 7, see Additional file 3). FISH analysis was applied sequentially using the pTa71 probe for the 45S rDNA locus and pIGS probe for the IGS type I on the pachytene chromosome of S. lycopersicum, S. lycopersicum var. cerasiforme, and S. pimpinellifolium. When S. lycopersicum chromosomes were hybridized to pTa71, a single strong signal was detected on the short arm of chromosome 2; however, in S. lycopersicum var. cerasiforme and S. pimpinellifolium, four signals were detected in four separated heterochromatic regions (Figure 7A, C, E). The number of signals was in accordance with the number of signals detected with the pIGS probe in S. lycopersicum. The number of foci detected on the pachytene chromosomes with the pIGS probe varied, as follows: four signals in S. lycopersicum, seven signals in S. lycopersicum var.cerasiforme, and six signals in S. pimpinellifolium (Figure 7B, D, F). All foci were located in regions of the pericentromeric heterochromatin. Most of the signals corresponding to the IGS repeat were stronger and more numerous than the signals from the pTa71 probe. However, the foci detected on the short arm of chromosome 2 of S. lycopersicum and S. pimpinellifolium were of a similar intensity, whether detected as with pIGS or pTa71. Taken together, the numbers of the satellite repeat vary dramatically across closely-related species and they can divide into two groups, repeats with or without rDNA coding sequence (additional file 4).

Figure 7
figure 7

Distribution of 45S rDNA and IGS homologous satellite repeat in tomatoes. The number of foci is variable even among closely related species. (A) S. lycopersicum with pTa71 probe; (B) S. lycopersicum with pIGS probe; (C) S. pimpinellifolium with pTa71 probe; (D) S. pimpinellifolium with pIGS probe; (E) S. lycopersicum var.cerasiforme with pTa71 probe; (F) S. lycopersicum var.cerasiforme with pIGS probe. Bar, 10 μm

Discussion

FISH analysis and the sequences of three BAC clones described in this study provide a good explanation on the origin and developmental procedures involved in the evolution of IGS-homologous satellite repeats because they contain both the original form and early stages of the variants in a genome. They also allow us to compare transitional sequences that make it possible to compare with previously proposed models.

Multiple mechanisms have been postulated to explain the development of satellite repeats, including unequal crossover, gene conversion, satellite transposition, illegitimate recombination, and segmental duplication [2, 5, 16, 21, 28, 29]. However, the origin and the development of the early stages of satellite repeat have remained unclear because transition sequences have not been identified. Our results showed that rDNA is the origin of the satellite repeat, and repeated rearrangement and retrotransposon insertion were involved in satellite repeat initiation (Figure 8). The inserted retrotransposon in 18S rDNA might serve as sites of unequal or ectopic recombination [30]. Unequal crossover seems to be commonly employed to multiply modified rDNA units such as TRRT inserted 18S rDNAs and IGS inserted 25S rDNAs. Doubling the number of IGS into the 25S rDNA seems a very effective way of amplifying a repeated sequence. Duplication of IGS into the 25S rDNA sequence may occur prior to the amplification of the type I subrepeat of the IGS of 25-18S rDNA. Because duplicated IGS in the 25S rDNA sequence is similar with normal IGS in appearance feature while IGSs of 25-18S rDNA were highly amplified and reorganized. Segmental duplication of the repeat was also one of the major mechanisms of expanding satellite repeat [21]. The copy number of type I subrepeat in expanded IGS of 25-18S rDNA was also expanded by segmental duplication (Figure 6C)

Figure 8
figure 8

Model for IGS homologous satellite repeat development. Ty3-Gypsy type LTR retrotransposon (TRRT) inserted in 18S rDNA might be integrated in the array at the early stage of divergence. The TRRT inserted 18S rDNA persist and increase its copy number through recombination. Duplication of IGS into 25S rDNA seems to be occurred prior to replace type II subrepeat to type I subrepeat. Homogenized satellite repeat may be developed through amplifying type I subrepeat with segmental duplication, but removing neighboring sequences such as 25S rDNA and type II subrepeat.

Homogenization of the satellite repeat to single-type monomers seems to be actively conducted from the initial stage of the satellite repeat. Gene conversion and unequal crossover has been proposed as the mechanism for genome-wide homogenization of the satellite repeats [5, 16, 31]. Our results are consistent with the models in that both gene conversion and unequal crossover play a pivotal role in regional homogenization of the satellite repeats on initial step. Analysis of the duplicated IGS located in the 25S rDNA sequence indicated that the two subrepeats are differentially regulated in their copy number. While the type I subrepeat was conserved in copy number, the number of type II subrepeat sequence was decreased in different length between variants implying that repeated unequal crossover has occurred. Neighboring sequences of type I subrepeat such as 25S rDNA and the retrotransoposon may be deleted as a result of unequal crossover homogenizing the expanding repeat (Figure 2B). However, the replacement of the type II subrepeat to the type I subrepeat within the expanded IGS of 25-18S rDNA can be explained by gene conversion. Despite the amplification of copy number of type I subrepeats, the size of type I subrepeat monomers was relatively well conserved as 52–53 bp. It appears to be a common phenomenon for the size of satellite monomers to conserve nearly uniform within a genome [13, 32, 33], whereas the copy number of satellite monomers varies dramatically across species, within an organism, or on a specific chromosome between different subspecies or varieties [7, 21, 34].

In addition, we have postulated that satellite repeats sharing the same monomer can become separated from the original locus. The four signals detected with the pTa71 probe in S. lycopersicum var.cerasiforme and S. pimpinellifolium indicate that there are coding sequences of rRNA at those four loci. This situation can be explained if an rDNA block had moved to a new location outside of their original array, followed by coding region deletion and type I subrepeat amplification. However, it is more likely that the amplified IGS repeat sequence translocated outside of the rDNA array with partial coding regions attached, than that transposed sequences on multiple chromosomes independently developed satellite repeats. These assumptions are supported by our FISH results in that, with the exception of foci on the short arm of chromosome 2, very weak or no signal was detected with a partial 18S rDNA probe in S. lycopersicum var. cerasiforme and S. pimpinellifolium. Of course, rDNA coding sequence could be effectively eliminated in the satellite repeat of S. lycopersicum. In S. lycopersicum var.cerasiforme and S. pimpinellifolium, but not yet in S. lycopersicum, additional IGS-homologous satellite repeats may have been made from the original three IGS homologous satellite repeats lacking rRNA coding sequence.

Tandem repeat segments could be moved by recombinational excision of looped-out modified rDNA segments during meiotic unequal alignment of the rDNA repeating units and might re-integrate into a new location [17, 28]. Indeed, this mechanism of recombinational deletion of amplified repeats was discovered in the fourth unit of rDNA in the clone, HBa0007F24. Of course, unequal recombination or illegitimate recombination could also lead to the translocation of satellite repeats [16, 35].

Although the integration of mobile elements in rDNA loci has been reported in animal taxa, new insertions are rapidly eliminated from the rDNA locus by unequal crossover between sister-chromatids [35, 36] and that new insertions are subject to random crossovers. Strong selective pressure against inactive rDNA units eliminates these insertions from the loci [37, 38]. The LTR retrotransposon in plant genome appears to have adapted differently to survive in the rDNA locus. When the TRRT randomly inserted into the 18S rRNA gene, unequal recombination might have occurred to make the array uniform. However, instead of removing the TRRT from the rDNA, it increased the number of genes modified with the insert. Evidence of the effort to remove the TRRT is also present in the sequence: Solo LTR unaccompanied by the rest of the transposon genes may be a remnant of a retrotransposon that has been removed by unequal recombination [3]. Similar constitution of sequence containing a rDNA variant and retrotransposons was found on the heterochromatin of Brassica rapa. [39]. In addition, amplified satellite repeats with sequence homologous to the IGS subrepeat of rDNA have been reported in several plants [12–15], although the preferential amplification or elimination of a repeat is highly variable across species. For example, the 2D8 satellite repeat found in potatoes consists of ~3 kb monomers of AT-rich and GC-rich subrepeat clusters, showing high sequence similarity with type I and AT-rich regions of IGS in rDNA [12]. The A1/A2 satellite repeat of tobacco has sequence similarity to regions downstream of the TIS of IGS [7].

Conclusion

We identified IGS-derived satellite repeats in tomato genome. By analyses of multiple transitional sequences, we clearly showed the origin and the growing procedure of the satellite repeat in tomato genome. Our results also suggested the molecular mechanisms of proliferation and homogenization of the satellite repeat in tomato genome by showing multiplication procedure of modified rDNA units and amplification/deletion of different subrepeats within the 45S rDNA.

Methods

Plant materials

Tomato plants (S. lycopersicum (LA3911), S. lycopersicum var. cerasiforme (LA4352), and S. pimpinellifolium (LA0417)) were grown in a controlled environment room at 26°C with 16 h light/8 h dark cycles.

DNA probes

All BAC clones used for FISH mapping were identified by screening a Heinz 608 BAC library, kindly provided by Drs. S. Tanksley and J. Giovannoni at Cornell University, Ithaca, NY, USA. BAC probes were labeled with digoxigenin-11-dUTP or biotin-16-dUTP by nick translation, according to the manufacture's protocol (Roche, Switzerland). Other DNA probes used for FISH analysis were: an rDNA probe, pTa71, containing the coding sequences for the 18S–25S rRNA genes from wheat [25], pIGS, TRRT, partial 18S rRNA. Probes were amplified from the clone, HBa0007F24 with the following primer pairs: partial 18S rRNA primer (F-AATGATCCTTCCGCAG GTTTACCT, R-GCTCTGTATACATTAGCATGGGATA), 25S rRNA primer (F-CGCCCTCCTACTCTTCGGGGCCTGG, R-CAGGTTAGGCGGCATTACCCGCTGA), Retrotransposon primer (F-AATGTA TACAGAGCATAGTGTGATGTCC, R-AGTGCTCAAAGAAAGCCTACTGTCACGG), and pIGS primer (F-CGACGTACCATTT GTGCTT, R-TTACCTATGGGCAGCACACATGGTC). PCR products were cloned using pGEM-T easy Vector system (Promega, Madison, WI) and sequenced.

Fluorescence in situhybridization (FISH)

The FISH procedure applied to both mitotic and meiotic chromosomes was the same as previously reported by Koo et al. [40]. In brief, chromosomal DNA on the slides was denatured with 70% formamide at 70°C for 2.5 min, followed by dehydration in a 70%, 85%, 95%, and 100% ethanol series at -20°C for 3 min each. The probe mixture, containing 50% formamide (v/v), 10% dextran sulfate (w/v), 5 ng/μl salmon sperm DNA and 500 ng/μl of labeled probe DNA, was heated at 90°C for 10 min then kept on ice for 5 min. Twenty microliters of this mixture was applied to the denatured chromosomal DNA and covered with a glass cover-slip. Slides were then placed in a humid chamber at 37°C for 18 h. Probes were detected with avidin-FITC and anti-digoxigenin Cy3 (Roche, Switzerland). Chromosomes were counterstained with 1 μg/μl DAPI (Sigma). The signals were detected with a Cooled CCD Camera (CoolSNAP, Photometrics, Pleasanton, CA). Images were processed with software (Meta Imaging Series™ 4.6) using Leica epi-fluorescence microscope equipped with FITC-DAPI two-way or FITC-rhodamine-DAPI three-way filter sets (Leica, Japan). The final printed images were prepared with Adobe Photoshop 7.0.

Fiber-FISH

Leaf nuclei were prepared, as described by Jackson et al. [41]. A suspension of nuclei was deposited at one end of a poly-L-lysine coated slide and permitted to air dry for 10 min. STE lysis buffer (8 μl) was pipetted on top of the nuclei and the slide incubated at room temperature for 4 min. A clean cover-slip was used to slowly drag the contents along the slide. The preparation was air dried, fixed in ethanol:glacial acetic acid (3:1) for 2 min, and baked at 60°C for 30 min. The DNA fiber preparation was incubated with a probe mixture, covered with a 22 mm × 40 mm cover-slip, and sealed with rubber cement. The slide was placed in an 80°C oven in direct contact with a heated surface for 3 min, transferred to a wet chamber, which was pre-warmed in an 80°C oven for 2 min, and transferred to 37°C overnight. Post-hybridization washing stringency was the same as for FISH on chromosome spreads. Signal detection was performed according to Koo et al. [40].

BAC sequencing and sequence assembly

A shotgun sequencing library was constructed in the pUC118 vector for average insert sizes of 3–5 kb. BigDye Terminator chemistry v3.1 (ABI, Foster City, CA) was used for the sequencing reactions. The sequences were analyzed using an ABI3730XL automatic DNA sequencer (ABI, Foster City, CA). All of the initial sequence data obtained were analyzed with the Phred/Phrap/Consed processing [42]. Base-calling and assembling for the individual sequences were conducted through the Phred/Phrap software. The value of the Phred scores of the sequences was 30 or higher. The completely assembled sequence was edited using Consed. Sequence editing for consensus contig formation was generated by visual confirmation, using the Sequencher 4.1.5 program (Gene codes Corp., Ann Arbor, USA)

Sequence annotation

Putative genes were primarily identified using BLASTN and BLASTX searches of GenBank http://www.ncbi.nlm.nih.gov/blast/Blast.cgi. The consensus sequence of IGS monomers previously reported by Schmidt-Puchata et al. [43] was used to search against the BAC sequences, HBa007F24, SLe089P21, and Sle049A24 by BLAST. Tandem repeats were identified using Tandem Repeat Finder [44] and edited manually in DNAMAN (Lynnon Corporation, Quebec, Canada). The monomers were extracted from IGS unit IV of HBa007F24 and then aligned using ClustalX [45]. The alignments were edited manually if necessary. The Neighbor-Joining trees were built using the Kimura two-parameter method [46]. MEGA 4 was employed to calculate pairwise transition and transversion mutations [47].

References

  1. Flavell RB, Bennett MD, Smith JB, Smith DB: Genome size and the proportion of repeated nucleotide sequence DNA in plants. Biochem Genet. 1974, 12: 257-269. 10.1007/BF00485947.

    Article  PubMed  Google Scholar 

  2. Ma J, Bennetzen JL: Rapid recent growth and divergence of rice nuclear genomes. Proc Natl Acad Sci USA. 2004, 101: 12404-12410. 10.1073/pnas.0403715101.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Ma J, Devos KM, Bennetzen JL: Analyses of LTR-retrotransposon structures reveal recent and rapid genomic DNA loss in rice. Genome Res. 2004, 14: 860-869. 10.1101/gr.1466204.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Pavelitz T, Rusche L, Matera AG, Scharf JM, Weiner AM: Concerted evolution of the tandem array encoding primate U2 snRNA occurs in situ, without changing the cytological context of the RNU2 locus. EMBO J. 1995, 14: 169-177.

    PubMed  PubMed Central  Google Scholar 

  5. Charlesworth B, Sniegowski P, Stephan W: The evolutionary dynamics of repetitive DNA in eukaryotes. Nature. 1994, 371: 215-220. 10.1038/371215a0.

    Article  PubMed  Google Scholar 

  6. Macas J, Navratilova A, Meszaros T: Sequence subfamilies of satellite repeats related to rDNA intergenic spacer are differentially amplified on Vicia sativa chromosomes. Chromosoma. 2003, 112: 152-158. 10.1007/s00412-003-0255-3.

    Article  PubMed  Google Scholar 

  7. Lim KY, Skalicka K, Koukalova B, Volkov RA, Matyasek R, Hemleben V, Leitch AR, Kovarik A: Dynamic changes in the distribution of a satellite homologous to intergenic 26-18S rDNA spacer in the evolution of Nicotiana. Genetics. 2004, 166: 1935-1946. 10.1534/genetics.166.4.1935.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Tek AL, Song J, Macas J, Jiang J: Sobo, a recently amplified satellite repeat of potato, and its implications for the origin of tandemly repeated sequences. Genetics. 2005, 170: 1231-1238. 10.1534/genetics.105.041087.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Hillis DM, Dixon MT: Ribosomal DNA: molecular evolution and phylogenetic inference. Q Rev Biol. 1991, 66: 411-453. 10.1086/417338.

    Article  PubMed  Google Scholar 

  10. Brown DD, Wensink PC, Jordan E: A comparison of the ribosomal DNA's of Xenopus laevis and Xenopus mulleri: the evolution of tandem genes. J Mol Biol. 1972, 63: 57-73. 10.1016/0022-2836(72)90521-9.

    Article  PubMed  Google Scholar 

  11. Flaquet J, Creusot F, Drone M: Molecular analysis of Phaseolus vulgaris rDNA unit and characterization of a satellite DNA homologous to IGS subrepeats. Plant Physiol and Biochem. 1997, 35: 611-622.

    Google Scholar 

  12. Stupar RM, Song J, Tek AL, Cheng Z, Dong F, Jiang J: Highly condensed potato pericentromeric heterochromatin contains rDNA-related tandem repeats. Genetics. 2002, 162: 1435-1444.

    PubMed  PubMed Central  Google Scholar 

  13. Volkov RA, Borisjuk NV, Panchuk II, Schweizer D, Hemleben V: Elimination and rearrangement of parental rDNA in the allotetraploid Nicotiana tabacum. Mol Biol Evol. 1999, 16: 311-320.

    Article  PubMed  Google Scholar 

  14. Maggini F, Cremonini R, Zolfino C, Tucci GF, D'Ovidio R, Delre V, DePace C, Scarascia Mugnozza GT, Cionini PG: Structure and chromosomal localization of DNA sequences related to ribosomal subrepeats in Vicia faba. Chromosoma. 1991, 100: 229-234. 10.1007/BF00344156.

    Article  PubMed  Google Scholar 

  15. Unfried K, Schiebel K, Hemleben V: Subrepeats of rDNA intergenic spacer present as prominent independent satellite DNA in Vigna radiata but not in Vigna angularis. Gene. 1991, 99: 63-68. 10.1016/0378-1119(91)90034-9.

    Article  PubMed  Google Scholar 

  16. Smith GP: Evolution of repeated DNA sequences by unequal crossover. Science. 1976, 191: 528-535. 10.1126/science.1251186.

    Article  PubMed  Google Scholar 

  17. Walsh JB: Persistence of tandem arrays: implications for satellite and simple-sequence DNAs. Genetics. 1987, 115: 553-567.

    PubMed  PubMed Central  Google Scholar 

  18. Navratilova A, Koblizkova A, Macas J: Survey of extrachromosomal circular DNA derived from plant satellite repeats. BMC Plant Biol. 2008, 8: 90-10.1186/1471-2229-8-90.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Cohen S, Yacobi K, Segal D: Extrachromosomal circular DNA of tandemly repeated genomic sequences in Drosophila. Genome Res. 2003, 13: 1133-1145. 10.1101/gr.907603.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Henikoff S, Ahmad K, Malik HS: The centromere paradox: stable inheritance with rapidly evolving DNA. Science. 2001, 293: 1098-1102. 10.1126/science.1062939.

    Article  PubMed  Google Scholar 

  21. Ma J, Jackson SA: Retrotransposon accumulation and satellite amplification mediated by segmental duplication facilitate centromere expansion in rice. Genome Res. 2006, 16: 251-259. 10.1101/gr.4583106.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Strachan T, Webb D, Dover GA: Transition stages of molecular drive in multiple-copy DNA families in Drosophila. EMBO J. 1985, 4: 1701-1708.

    PubMed  PubMed Central  Google Scholar 

  23. Koo DH, Jo SH, Bang JW, Park HM, Lee S, Choi D: Integration of cytogenetic and genetic linkage maps unveils the physical architecture of tomato chromosome 2. Genetics. 2008, 179: 1211-1220. 10.1534/genetics.108.089532.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Mueller LA, Tanksley SD, Giovannoni JJ, van Eck J, Stack S, Choi D, Kim BD, Chen M, Cheng Z, Li C, et al: The Tomato Sequencing Project, the First Cornerstone of the International Solanaceae Project (SOL). Comp Funct Genomics. 2005, 6: 153-158. 10.1002/cfg.468.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Gerlach WL, Bedbrook JR: Cloning and characterization of ribosomal RNA genes from wheat and barley. Nucleic Acids Res. 1979, 7: 1869-1885. 10.1093/nar/7.7.1869.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Capy P: Classification and nomenclature of retrotransposable elements. Cytogenet Genome Res. 2005, 110: 457-461. 10.1159/000084978.

    Article  PubMed  Google Scholar 

  27. Peralta IE, Spooner DM: Classification of wild tomatoes: a review. Kurtziana. 2000, 28: 45-54.

    Google Scholar 

  28. Alexandrov IA, Mitkevich SP, Yurov YB: The phylogeny of human chromosome specific alpha satellites. Chromosoma. 1988, 96: 443-453. 10.1007/BF00303039.

    Article  PubMed  Google Scholar 

  29. Mestrovic N, Castagnone-Sereno P, Plohl M: Interplay of selective pressure and stochastic events directs evolution of the MEL172 satellite DNA library in root-knot nematodes. Mol Biol Evol. 2006, 23: 2316-2325. 10.1093/molbev/msl119.

    Article  PubMed  Google Scholar 

  30. Kumar A, Bennetzen JL: Plant Retrotransposons. Annu Rev Genet. 1999, 33: 479-532. 10.1146/annurev.genet.33.1.479.

    Article  PubMed  Google Scholar 

  31. Stephan W: Recombination and the evolution of satellite DNA. Genet Res. 1986, 47: 167-174.

    Article  PubMed  Google Scholar 

  32. Willard HF, Waye JS: Chromosome-specific subsets of human alpha satellite DNA: analysis of sequence divergence within and between chromosomal subsets and evidence for an ancestral pentameric repeat. J Mol Evol. 1987, 25: 207-214. 10.1007/BF02100014.

    Article  PubMed  Google Scholar 

  33. Zhang Y, Huang Y, Zhang L, Li Y, Lu T, Lu Y, Feng Q, Zhao Q, Cheng Z, Xue Y, et al: Structural features of the rice chromosome 4 centromere. Nucleic Acids Res. 2004, 32: 2023-2030. 10.1093/nar/gkh521.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Cheng Z, Dong F, Langdon T, Ouyang S, Buell CR, Gu M, Blattner FR, Jiang J: Functional rice centromeres are marked by a satellite repeat and a centromere-specific retrotransposon. Plant Cell. 2002, 14: 1691-1704. 10.1105/tpc.003079.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Eickbush TH, Eickbush DG: Finely orchestrated movements: evolution of the ribosomal RNA genes. Genetics. 2007, 175: 477-485. 10.1534/genetics.107.071399.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Kojima KK, Kuma K, Toh H, Fujiwara H: Identification of rDNA-specific non-LTR retrotransposons in Cnidaria. Mol Biol Evol. 2006, 23: 1984-1993. 10.1093/molbev/msl067.

    Article  PubMed  Google Scholar 

  37. Kojima KK, Fujiwara H: Evolution of target specificity in R1 clade non-LTR retrotransposons. Mol Biol Evol. 2003, 20: 351-361. 10.1093/molbev/msg031.

    Article  PubMed  Google Scholar 

  38. Long EO, Dawid IB: Repeated genes in eukaryotes. Annu Rev Biochem. 1980, 49: 727-764. 10.1146/annurev.bi.49.070180.003455.

    Article  PubMed  Google Scholar 

  39. Lim KB, Yang TJ, Hwang YJ, Kim JS, Park JY, Kwon SJ, Kim J, Choi BS, Lim MH, Jin M, et al: Characterization of the centromere and peri-centromere retrotransposons in Brassica rapa and their distribution in related Brassica species. Plant J. 2007, 49: 173-183. 10.1111/j.1365-313X.2006.02952.x.

    Article  PubMed  Google Scholar 

  40. Koo DH, Plaha P, Lim YP, Hur Y, Bang JW: A high-resolution karyotype of Brassica rapa ssp. pekinensis revealed by pachytene analysis and multicolor fluorescence in situ hybridization. Theor Appl Genet. 2004, 109: 1346-1352. 10.1007/s00122-004-1771-0.

    Article  PubMed  Google Scholar 

  41. Jackson SA, Wang ML, Goodman HM, Jiang J: Application of fiber-FISH in physical mapping of Arabidopsis thaliana. Genome. 1998, 41: 566-572. 10.1139/gen-41-4-566.

    Article  PubMed  Google Scholar 

  42. Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8: 186-194.

    Article  PubMed  Google Scholar 

  43. Schmidt-Puchta W, Gunther I, Sanger HL: Nucleotide sequence of the intergenic spacer (IGS) of the tomato ribosomal DNA. Plant Mol Biol. 1989, 13: 251-253. 10.1007/BF00016143.

    Article  PubMed  Google Scholar 

  44. Benson G: Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999, 27: 573-580. 10.1093/nar/27.2.573.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16: 111-120. 10.1007/BF01731581.

    Article  PubMed  Google Scholar 

  47. Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by grants from Crop Functional Genomics Center (CG1221) of the 21st Century Frontier Research Program funded by MOEST of the Korean government. We thank to Solanaceae Genome Network (Drs. Jim Giovanoni and Steve Tanksley) at Cornell University for providing all the tomato BAC resources.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Doil Choi.

Additional information

Authors' contributions

SJ, SL, JFK, DC complied results and wrote the manuscript. DK carried out all FISH analyses. SJ, CH, TY, SK carried out the bioinformatics analyses. DC directed the project. All authors read and approved the final manuscript.

Sung-Hwan Jo, Dal-Hoe Koo contributed equally to this work.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Jo, SH., Koo, DH., Kim, J.F. et al. Evolution of ribosomal DNA-derived satellite repeat in tomato genome. BMC Plant Biol 9, 42 (2009). https://doi.org/10.1186/1471-2229-9-42

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2229-9-42

Keywords