Highly diverse chromoviruses of Beta vulgaris are classified by chromodomains and chromosomal integration
© Weber et al.; licensee BioMed Central Ltd. 2013
Received: 11 September 2012
Accepted: 22 January 2013
Published: 1 March 2013
Chromoviruses are one of the three genera of Ty3-gypsy long terminal repeat (LTR) retrotransposons, and are present in high copy numbers in plant genomes. They are widely distributed within the plant kingdom, with representatives even in lower plants such as green and red algae. Their hallmark is the presence of a chromodomain at the C-terminus of the integrase. The chromodomain exhibits structural characteristics similar to proteins of the heterochromatin protein 1 (HP1) family, which mediate the binding of each chromovirus type to specific histone variants. A specific integration via the chromodomain has been shown for only a few chromoviruses. However, a detailed study of different chromoviral clades populating a single plant genome has not yet been carried out.
We conducted a comprehensive survey of chromoviruses within the Beta vulgaris (sugar beet) genome, and found a highly diverse chromovirus population, with significant differences in element size, primarily caused by their flanking LTRs. In total, we identified and annotated full-length members of 16 families belonging to the four plant chromoviral clades: CRM, Tekay, Reina, and Galadriel. The families within each clade are structurally highly conserved; in particular, the position of the chromodomain coding region relative to the polypurine tract is clade-specific. Two distinct groups of chromodomains were identified. The group II chromodomain was present in three chromoviral clades, whereas families of the CRM clade contained a more divergent motif. Physical mapping using representatives of all four clades identified a clade-specific integration pattern. For some chromoviral families, we detected the presence of expressed sequence tags, indicating transcriptional activity.
We present a detailed study of chromoviruses, belonging to the four major clades, which populate a single plant genome. Our results illustrate the diversity and family structure of B. vulgaris chromoviruses, and emphasize the role of chromodomains in the targeted integration of these viruses. We suggest that the diverse sets of plant chromoviruses with their different localization patterns might help to facilitate plant-genome organization in a structural and functional manner.
KeywordsTy3-gypsy retrotransposon Chromovirus CR Chromodomain Heterochromatin Beta vulgaris
In plants, long terminal repeat (LTR) retrotransposons constitute the most abundant class of transposable elements (TEs), and represent more than 50% of the genome in some species [1, 2]. Their massive amplification and dynamic nature play a central role in the organization, function, and evolution of plant genomes [3–5]. The ‘copy and paste’ propagation by reverse transcription of RNA intermediates can rapidly increase the copy number of LTR retrotransposons and, hence, the genome size [6, 7]. However, the excision of LTR retrotransposons through diverse recombination processes counteracts genome expansion, and thus contributes to the dynamic balance of host genome size [8, 9].
The LTRs flanking the coding region contain promoter motifs that drive transcription. The proteins essential for the reverse transcription and integration of a new LTR retrotransposon copy are encoded as a gag-pol polyprotein organized in one or two open reading frames (ORFs). Based on gene arrangement, plant LTR retrotransposons can be further sub-classified into the Ty1-copia order (Pseudoviridae) and the Ty3-gypsy order (Metaviridae), with the latter having the same domain arrangement as retroviruses (Orthoviridae) [10, 11].
The physical mapping of LTR retrotransposons by fluorescence in situ hybridization (FISH) and the accessibility of complete plant-genome sequences have yielded considerable information on their chromosomal and genomic organization. Thereby, a non-random genomic distribution with accumulation in particular chromosomal regions is apparent, indicating different target specificities . Several studies have identified a role for the integrase in different patterns of retrotransposon integration [13, 14]. A specific lineage of Ty3-gypsy retrotransposons is characterized by the presence of a chromodomain (chr omatin o rganization mo difier domain) at the C-terminal region of the integrase, and members of this lineage are referred to as chromoviruses [15, 16]. Proteins containing chromodomains are strongly related to chromatin modification and gene regulation , and are able to interact with proteins, RNA, and DNA . Recent studies have shown that the chromodomain specifies the target site preference of LTR retrotransposons by the recognition of characteristic chromatin modifications [19, 20]. Chromoviruses represent the ancient and most widespread lineage of Ty3-gypsy retrotransposons, and the four chromoviral clades, Tekay, CRM, Galadriel, and Reina, are widely distributed throughout gymnosperms and angiosperms .
Recently, the chromosomal and genomic organization of centromere-specific chromoviruses in the wild beet genus Patellifolia, which is closely related to the beet genus Beta, were described . The genus Beta belongs to the family Amaranthaceae and is subdivided into the sections Beta, Corollinae, and Nanae. All cultivated species (sugar beet, fodder beet, garden beet, and chard) belong exclusively to the section Beta. The Beta vulgaris (sugar beet) genome is 758 Mb in size  and contains at least 63% repetitive sequences . The access to early assembly versions for a draft of the B. vulgaris genome sequence enabled bioinformatic identification of complete chromoviruses. These data have facilitated large-scale studies of the structure, variability, and evolution of LTR retrotransposons in the genus Beta to complement and extend the knowledge of the repetitive DNA fraction [21, 25, 26].
In this study, we characterize all four chromoviral clades populating the genome of B. vulgaris. Based on the amino acid composition of their chromodomain, all B. vulgaris chromoviruses can be attributed to one of two distinct groups. We investigated the evolution of the chromoviral clades within the beet genera Beta and Patellifolia, and found that the chromoviruses are highly divergent. Their chromosomal distribution patterns might result from a targeted integration mediated by different chromodomains.
Structure of the chromoviruses in Beta vulgaris
To analyze B. vulgaris chromoviral retrotransposons, we conducted tBLASTn searches against contigs of the B. vulgaris genome sequence draft, using domains of the gag-pol polyproteins of the chromoviral Tekay, CRM, Galadriel, and Reina clades.
Chromoviruses of Beta vulgaris
Total length, bp
LTR length 5′/3′, bp
Beetle 3 a
Beetle 4 a
Beetle 5 a
Beetle 6 a
Beetle 7 a
Bongo 1 a
Bongo 2 a
Bongo 3 a
Bingo 1 a
Bingo 2 a
Bingo 3 a
Bingo 4 a
Bingo 5 a
Bingo 6 a
Bingo 7 a
Beon 1 a
Members of the Galadriel clade were detectable only using a refined tBLASTn query based on the Galadriel integrase amino acid sequence . A complete internal chromoviral sequence, harboring no stop codons or frameshifts, but lacking the LTRs, was identified, and the reverse transcriptase (RT) gene sequence was used as a probe to detect the sugar-beet bacterial artificial chromosome (BAC) 54O24 by hybridization to high-density filters . The Galadriel sequence was obtained by an adaptor primed suppression PCR method  using the DNA of BAC 54O24 as template. The 21 retrotransposons belong to the 4 plant chromoviral clades, and were designated as Beetle (CRM clade), Bingo (Reina clade), Bongo (Tekay clade), and Beon (Galadriel clade).
Furthermore, the Ty3-gypsy retrotransposon SCHMIDT, previously detected in a B. vulgaris BAC sequence  and by sequence analysis of genomic clones from a c o t-1 library of B. vulgaris, represents an intact chromovirus of the Tekay clade. In accordance with the nomenclature of Beta chromoviruses, we assigned SCHMIDT to the Bongo chromoviruses, and designated it as Bongo 3.
In Bongo 3, the putative functional PPT is located immediately upstream of the 3′ LTR, and an additional PPT (PPTa) is located downstream of the chromodomain and overlaps with the TAA stop codon of the gag-pol ORF (Figure 1). Bongo 1 and Bongo 2, as well as members of other clades, do not contain duplicated PPTs.
Genome-wide analysis of chromoviral reverse transcriptases
Chromodomains distinguish the chromoviral clades
The position of the gag-pol ORF termini and the localization of the chromodomain in the integrase with reference to the PPT are the major characteristics facilitating the separation of chromoviruses into clades. Previous studies have shown that the gag-pol ORF of the CRM clade members extends into the 3′ LTR, with the chromodomain situated downstream of the PPT [21, 35]. For the remaining chromoviral clades, the chromodomain is located upstream of the PPT. In the chromoviruses Bingo 1 to 4 and Bingo 7, an extension of the ORF into the flanking LTR was found, whereas in Bingo 5 and Bingo 6, the gag-pol ORFs were seen to terminate upstream of the 3′ LTR.
A comprehensive investigation of the predicted integrase amino acid sequence harboring the chromodomain was carried out. The integrase sequence could be separated into two regions (see Additional file 2; Figure S2). The 5′ region spans the integrase region between the ZF and the GPF/GPY motif (positions 1 to 313) and is highly conserved throughout plant chromoviruses. The typical ZF HHCC is followed by the D, D35E motif. There is a marked enrichment of positively charged residues upstream of the GPF/GPY motif. In particular, the Beetle families are enriched in the amino acids arginine (K) and lysine (R). Interestingly, the proline (P) upstream of the RK-rich region implies similarity to nuclear localization signals (NLSs) with the motif PxR/KKxK.
Genomic organization of Beta chromoviruses
The strong signals indicated the high abundance of Beetle and Bongo retrotransposons (Figure 5A, B). By contrast, Bingo and Beon chromoviruses had either lower copy numbers or more diverged members (Figure 5C, D). Signals over a wide range of molecular weights were visible for Beetle, Bongo, and Bingo (Figure 5A–C), but not for Beon chromoviruses (Figure 5D). However, the presence of strongly hybridizing fragments in all clades indicates the existence of multiple copies with conserved restriction sites across most species. A similar banding pattern was detected in the cultivated sugar beet B. vulgaris ssp. vulgaris (lane 1) and the wild beet B. patula (lane 2), both belonging to the section Beta. Apart from Beetle 3 (Figure 5A), similar results were found for B. corolliflora (lane 3) and B. nana (lane 4), with conserved fragments in most of the autoradiograms, suggesting a close relationship of the various chromoviruses. Additionally, the wild beet Patellifolia procumbens (lane 5), formerly assigned to the genus Beta, and the more distant spinach (Spinacia oleracea) (lane 6) were analyzed as outgroups. The genomes of both species also contain chromoviruses of the clades CRM, Tekay, Galadriel and Reina.
The Southern blotting experiments indicate the widespread distribution of chromoviruses within the order Amaranthaceae. The differences in copy number are clade-specific rather than species-specific. The Bongo families of the Tekay clade showed the highest amplification, whereas Galadriel clade members were much less abundant throughout all analyzed species. This is consistent with the genome-wide analysis of RT sequences of B. vulgaris, for which more than 60% of the RTs were found to belong to the Tekay clade, followed by 23% and 14% for the CRM and Reina clades, respectively.
Estimation of the transposition time
To estimate the transposition time of individual chromoviral copies, we used the divergence of the LTR sequences, as these pairs are usually identical upon integration . Assessing the 22 full-length chromoviruses, we found that elements of the Reina clade had inserted most recently into the B. vulgaris genome, with an average transposition time of 0.35 Mya. Similarly, the Bingo families transposed between 26,000 and 600,000 years ago (Table 1). The average transposition times estimated for the CRM clade families and the Tekay chromoviruses were approximately equal (one million years ago). For the estimation of the transposition time of Beetle 6 and Bongo 2, which contain one partially deleted LTR, the matching LTR regions were used, resulting in respective timeframes of 0.97 and 1.29 million years ago.
Chromosomal localization of chromoviruses
The LTR probe of Beetle 7 hybridized to the centromeric and peri-centromeric heterochromatin of most chromosomes (Figure 6A).
The FISH image of Bongo 3 showed a pattern of widely dispersed signals along the length of all chromosome arms, including the intercalary and centromeric heterochromatin of both arms (Figure 6B); however, reduced hybridization to some centromeric regions was seen (Figure 6B, arrows). Strong signals were detected on both chromatids, indicating clusters most likely originating from nested organization. The results are consistent with the Southern hybridization results, which classified the Bongo3 family as a highly abundant component of the B. vulgaris genome. The Bingo 1-LTR probe produced signals that were predominantly located close to the pericentromeric regions of most chromosomes, although there were considerable differences in the strength of the signals between chromosomes (Figure 6C). The LTR probe of Beon 1 hybridizes to only one chromosome pair (Figure 6D). Multicolor FISH on interphase and pachytene chromosomes indicated the interspersion of Beon 1 copies (red fluorescence) in the 18S rRNA gene arrays (green fluorescence), which are located on chromosome 1 (Figure 6E, F). High-resolution fiber-FISH on extended genomic fibers revealed the genomic organization of rDNA arrays (green signals) within the genome of B. vulgaris. Two different types of rDNA subrepeats were identified: 1) those representing the canonical rDNA repeat without interruption by retrotransposon sequences and 2) those harboring Beon 1 (Figure 6G). Both repeat types were found to be organized in interspersed sequence stretches.
Transcriptional activity of Beta chromoviruses
More than 21,500 expressed sequence tags (ESTs; http://genomics.msu.edu/sugarbeet/blast.html) from sugar beet were screened for potential chromoviral transcripts. In general, chromoviruses seem to be barely expressed in the genome of B. vulgaris. Only 10 evaluable chromoviral transcripts were detected, which belonged to the CRM and Tekay clades. These ESTs had identity values greater than of 80% and were larger than 600 bp. Only four transcripts were associated with full-length chromoviral sequences; the others were generated from copies with internal deletions or they represent chimeric structures containing unrelated sequences. No ESTs were identified for Beon 1 or for one of the Bingo families.
A large amount of genomic data has been generated by high-throughput sequencing technologies, which enable comparative analyses of retrotransposons within or across genomes. Short-read assemblies offer new possibilities for the detection of TEs; however, they are often excluded from assemblies because of their repetitiveness, genome-wide distribution, truncation, and nested organization [1, 8, 41, 42].
We carried out a comparative analysis of chromoviruses in the beet genome, using data from a preliminary B. vulgaris genome draft sequence, which provides an insight into the genomic and chromosomal organization, distribution, and evolution of these viruses.
Chromoviruses are widely distributed within genera of the Amaranthaceae family
Although the coding region of chromoviruses is relatively conserved within a clade, the sequence identity between the LTRs has been used for the grouping of the 22 chromoviruses into separate families . In total, 16 different chromoviral families of the CRM, Reina, Tekay, and Galadriel clades have been identified, representing all known chromoviral clades in higher plants. Nevertheless, our analysis of genome-wide RT sequences clearly shows that more chromoviruses may be identified from the genome draft sequence, which is also evident from the signal strength in the blotting and in situ hybridization experiments. Only one similar, but less exhaustive analysis of chromoviruses has been conducted previously, which was for the Musa acuminata genome .
Our analysis shows that families from the same clade have conserved characteristics throughout the genera Beta and Patellifolia. In cultivated and wild beets, the Bongo families of the Tekay clade were the most abundant, followed by the Beetle families of the CRM clade. In comparison, Beon 1 of the Galadriel clade most probably comprises a single family harboring only a few elements. Similarly, all M. acuminata Galadriel chromoviruses are members of a single family, designated Monkey ; however, Monkey members constitute about 0.2 to 0.5% of the genome . In M. acuminata, the Reina clade constitutes more than half of all chromoviruses and makes up about 4% of the genome, followed by the Tekay clade, constituting about 2% .
Chromosomal localization of Beta chromoviruses
The localization of plant chromoviruses on chromosomes has been investigated in several plants, focusing on members of the CRM clade that show a specific accumulation in centromeric regions. In banana, the chromoviruses of the Reina and Tekay clade were physically mapped in centromeric and peri-centromeric associated heterochromatin, whereas Monkey elements belonging to the Galadriel clade preferentially inserted into the NOR and co-localized with the rRNA genes . This was also observed in this study for Beon 1 in B. vulgaris, where the physical mapping indicated exclusive localization of Beon 1 copies in the NOR. The results suggest that the localization might have been established by a single integration event of a Beon 1 copy into the 18S rRNA gene of a common ancestor, as supported by the presence of Beon 1 elements within the 18S rRNA genes of wild beets from the sections Beta and Corolliflora. Such integrations in the NOR region are not unusual for TEs, and have been reported for the long interspersed element (LINE) R2Bm of Bombyx mori. Furthermore, B. vulgaris ribosomal RNA genes seem to tolerate TE integrations, as an insertion of a single BNR1 LINE has previously been reported . Nevertheless, our FISH analysis clearly showed multiple Beon 1 copies interrupting the 18S rRNA genes. The integration of several members of the tomato rDNA-related retrotransposon (TRRT) family within 18S rRNA genes was recently shown, and the existence of segmental duplication events rather than targeted integration has been proposed to explain retrotransposon amplification . Our studies assigned TRRT elements to the Galadriel clade, forming a branch together with Beon 1. The sequence conservation of Beon 1 copies might result from strong purifying selection and homogenization of the coding sequences of the ribosomal DNA.
Integrations into genic regions were shown for families of the Reina and Tekay clade. A Bongo 3 copy was identified in the vicinity of a disease resistance-activating factor , and the recently described chromovirus Bert (in this paper assigned to Bingo elements) was found within an intron of the callose synthase gene . It is possible that these plant retrotransposons have the ability to alter the expression of nearby genes, as the importance of TEs for the epigenetic regulation of plant genomes has been stated previously [50–52].
The functional role of the integrase
Several studies have confirmed the functional role of integrase for targeted integrations [13, 14, 53]. During the integration process, the interaction of DNA and integrase is crucial , and the relevance of the chromodomain in modulating the interaction with diverse chromatin components has already been shown [19, 55]. An enrichment of positively charged amino acid residues was found within the integrase of all Beta chromoviruses. The differences in the extent of these residues possibly regulate the degree to which the integrase is capable of establishing electrostatic protein–protein interaction. This might enable the retrotransposon to sense target-specific chromatin states, as described by Roudier et al. . Furthermore, we detected a potential NLS in all Beetle families of the CRM clade . Such NLS signals are part of the integrase of several retroviruses and LTR retrotransposons, and are responsible for the transfer of the pre-integration complex to the nucleus [58–60].
The chromodomains of the plant chromoviral clades CRM, Reina, Tekay, and Galadriel in B. vulgaris are easily distinguishable by the presence or absence of conserved amino acid residues compared with HP1, and their position in relation to the gag-pol polyprotein. Although the chromodomain sequences have been classified into three groups [19, 38], only group II chromoviruses have been identified in plants, with all of them belonging to the clades Tekay, Reina, and Galadriel.
The chromodomain encoded by CRM clade members such as Beetle extends into the 3′ LTR [21, 35]. Because of this substantial difference from group II chromodomains, the Beetle chromodomains are referred to as CR motifs. A recent survey of CRM clade elements across diverse plant genomes assigned these chromoviruses to three different groups . CRM chromoviruses of group A carry a CR motif, and are genuine centromeric retrotransposons, which probably transpose actively into centromeric regions. By contrast, group B members are not localized at the centromere, whereas group C representatives, despite a lack of the CR motif, were also found in centromeric regions. Interestingly, the P. procumbens chromovirus Beetle 1  and the Beetle 7 chromovirus from B. vulgaris described in the current study share considerable amino acid identity (41%) within their chromodomain, and also have similar LTR lengths (1089 and 1086 bp, respectively), indicating a role for both Beetle 1 and Beetle 7 in the formation of functional centromeres. Amino acid conservations were also found within the chromodomains of CR elements from different grass species [40, 61], supporting the assumption that the CR elements of grasses were derived from a single ancient family , and that conservation of the chromodomain is also crucial for centromere stability and thus host genome integrity.
Transcriptional activity of Beta chromoviruses
Centromeres are thought to be determined epigenetically , including via a transcription-mediated mechanism [51, 64]. Several RT-PCR studies have identified the transcriptional activity of chromoviruses, in particular of CRM clade members [21, 40, 64–66]. The analysis of B. vulgaris EST datasets indicates the capability of Beetle chromoviruses for autonomous transposition. Thus, the chromodomain as a key component of genuine CRs facilitates the targeting process into centromeric regions, and might be therefore responsible for the generation of centromeric transcripts, which are involved in RNA interference-mediated centromere identity and function.
The rRNA genes have high transcriptional activity, thus it is possible that read-through transcripts of Beon 1 might be generated. Alternatively, as the Beon 1 copies harbor intact coding sequences, their reverse transcription and the integration of new copies into the genome is conceivable. However, corresponding transcripts were not detected in the EST database. Hence, epigenetic silencing mechanisms might prevent the reverse transcription and spreading of Beon 1 copies into other chromosomal regions. This could be caused by the insertion of Beon 1 in two orientations, as was shown in wild beet species, whereas transcription would result in double-stranded RNA, which would immediately initiate the RNA interference machinery. Subsequently, Dicer-generated small interfering RNAs would serve as substrates in RNA-induced transcriptional silencing (RITS) complexes or RNA-induced silencing complex (RISC). RITS would initiate the transcriptional silencing of Beon 1 copies by RNA-directed DNA methylation. Most likely, is the post-transcriptional silencing of Beon 1 copies ed by degradation of mRNA mediated by RISC. In plants, it has been shown that rDNA transcription is subject to dosage control , with only a subset of rDNA genes being transcribed. It might be possible therefore that rDNA genes containing Beon 1 copies are not transcribed.
Based on the accumulation of mutations within their LTRs , we calculated the age of the transposition events and concluded that members of the 16 families transposed less than two million years ago. Therefore, these transpositions are evolutionarily recent events. However, families of the four clades are likely to be much older, as deduced from their widespread distributions within the genera Beta and Patellifolia. Furthermore, nearly 70% of the analyzed contigs of the draft sequence of the beet genome contain incomplete or recombined copies, which over time have lost the typical retrotransposon hallmarks. Ma et al. found that LTR retrotransposons are subject to a genome-specific recombination rate that results in a half-life of less than 6 million years in rice or 3 million years in Arabidopsis.
Chromoviruses have been shown to make up a large fraction of LTR retrotransposons in the genome of B. vulgaris and related wild beets. In our study, we focused on complete members of all four chromoviral clades known in higher plants. Intact members could still be capable of autonomous transposition, which is in agreement with the relatively recent transposition for single elements such as Bingo 1, Bongo 3, and Beetle 7. Based on representative members of the four plant chromoviral clades, we were able to show their widespread presence within the family of Amaranthaceae, indicating their ancient origin. The analysis of the chromodomain-containing Ty3-gypsy retrotransposons provides valuable information for the annotation of the repeated DNA fraction of Beta genomes and is important for the understanding of the contribution of the chromodomain to retrotransposon guidance.
Plant material and preparation of genomic DNA
Plants of B. vulgaris ssp. vulgaris KWS2320, B. patula (BETA 548), B. corolliflora (BETA 846), B. nana (BETA 541), P. procumbens (BETA 951), and S. oleracea ‘Matador’ were grown under greenhouse conditions. Wild beet seeds were obtained from Dr L Frese (Julius Kühn Institute (JKI), Federal Centre for Breeding Research on Cultivated Plants, Quedlinburg, Germany, and from the Genbank of the Plant Genome Resources Center, Gatersleben, Germany), and genomic DNA was isolated from the young leaves using the cetyltrimethyl ammonium bromide (CTAB) protocol .
A local database of a preliminary B. vulgaris genome assembly was queried using clade-specific chromovirus gag-pol polyproteins. The unedited and non-public assembly RefBeet-0.1.1 draft comprises 628 Mb in 346,000 contigs with an N50 contig size of approximately 4000 bp (a current version of the sugar-beet genome draft (RefBeet-0.9) is available for download at http://bvseq.molgen.mpg.de). Additionally, for the in silico identification of chromoviral sequences, 20.4 Mbp of BAC end-sequence data from US H20 clones  were used.
For the identification of chromoviral RTs, an HMM-based strategy was used with hmmbuild and hmmsearch programs of the HMMER3 package (hmmer.janelia.org, ). Based on a Ty3-gypsy RT alignment containing 96 references of diverse eukaryotes (described by Llorens et al. , and downloadable from http://www.gydb.org), an HMM specific for Ty3-gypsy RTs was constructed and used to query a computationally translated B. vulgaris genome assembly. Hits with a bit score greater than 50 and a length of more than 200 amino acids were selected, aligned using the MUSCLE algorithm  and compared with the previously published chromoviral references . All hits that could be assigned to a chromoviral clade were kept for further analysis.
Phylogenetic analyses were conducted using MEGA version 4 . The constructed neighbor-joining tree of the ClustalW alignment of the complete gag-pol polyprotein sequences was calculated, using bootstrap values of 1000 cycles.
Integration times were calculated using the equation t = K/2r, where t is the age, K is the number of nucleotide substitutions per site between each LTR pair, and r is the nucleotide substitution rate. An average synonymous substitution rate of 1.5 × 10-8 mutations/site/ year as determined for the chalcone synthase and Adh loci in A. thaliana was used.
PCR analyses, cloning, and sequencing
Primers for amplification of specific chromoviral domains, LTR sequences and integration sites
Beetle 7 LTR
Beetle 3 RT
Bongo 3 LTR
Bongo 3 RT
Bingo 1 LTR
Bingo 1 RT
Beon 1 LTR
Beon 1 RT
For gel electrophoresis, genomic DNA was digested with the endonuclease Hin dIII and separated in 1.2% agarose gels, then transferred onto positively charged nylon membranes (GE Healthcare, Princeton, NJ, USA) using alkaline transfer. Southern hybridizations using 32P-labeled probes were performed in accordance with standard protocols . Filters were hybridized at 60°C overnight, then washed at 60°C in 2 × SSC/0.1% SDS and 1 × SSC/0.1% SDS for 10 minutes each. Signals were detected by autoradiography.
Fluorescence in situ hybridization
The meristem of young leaves was used for the preparation of mitotic chromosomes. Before fixation in methanol:acetic acid (3:1), leaves were incubated for 4 hours in 2 mmol/l 8-hydroxyquinoline. Fixed plant material was macerated in an enzyme mixture consisting of 0.3% (w/v) cytohelicase, 1.8% (w/v) cellulase from Aspergillus niger (Sigma-Aldrich Chemie GmbH, Taufkirchen, Germany), 0.2% (w/v) cellulase Onozuka-R10 (SERVA Electrophoresis GmbH, Heidelberg, Germany) and 20% (v/v) pectinase from A. niger, (Sigma-Aldrich Chemie GmbH) followed by dropping the nuclei suspension onto slides as described previously . FISH was performed on B. vulgaris chromosomes in accordance with the protocol of Heslop-Harrison et al. , modified for beet by Schmidt et al. . For high-resolution FISH, the chromatin fibers were prepared as described previously . The LTR fragments of Beetle 7, Bongo 3, Bingo 1, Beon 1 and the 18S rRNA-specific clone pZR18S [GenBank: HE578879] were used as probes after labeling by PCR in the presence of biotin-11-dUTP or digoxigenin-11-dUTP. Chromosome preparations were counterstained with 4′,6′-diamidino-2-phenylindole (DAPI) and mounted in antifade solution (CitiFluor Ltd, Leicester, UK). Slides were examined under a fluorescence microscope (Axioplan 2 imaging; Carl Zeiss, Jena, Germany equipped with filter numbers 02 (DAPI) ,09 (FITC), and15 (Cy3). Images were acquired directly with Applied Spectral Imaging software (v.3.3; Applied Spectral Imaging, Carlsbad, CA, USA) coupled with the high-resolution CCD camera ASI BV300-20A.
Accession numbers of analyzed chromoviruses and reference HP1 proteins
GenBank accession number
Bacterial artificial chromosome
Fluorescence in situ hybridization
Hidden Markov Model
Like heterochromatin protein (Arabidopsis homolog)
Long terminal repeat
Nuclear localization signal
Nucleolus organizer region
Open reading frame
Primer binding site
Sodium dodecyl sulphate
Saline sodium citrate
Standard tris acetate buffer
Target site duplication
We thank Ines Walter for excellent technical assistance in FISH analyses. This work was funded by the BMBF grant “Verbundprojekt GABI BeetSeq: Erstellungeiner Referenzsequenz für das Genom der Zuckerrübe (Beta vulgaris)”, FKZ 0315069A and 0315069B (to H.H. and B.W.) and by the BMBF grant “AnnoBeet: Annotation des Genoms der Zuckerrübe unter Berücksichtigung von Genfunktionen und struktureller Variabilität für Nutzung von Genomdaten in der Pflanzenbiotechnologie.", FKZ 0315962 A, 0315962 B and 0315962 C (To B.W, H.H., and T.S.).
- SanMiguel P, Tikhonov A, Jin YK, Motchoulskaia N, Zakharov D, Melake-Berhan A, Springer PS, Edwards KJ, Lee M, Avramova Z, Bennetzen JL: Nested retrotransposons in the intergenic regions of the maize genome. Science. 1996, 274: 765-768.View ArticlePubMed
- SanMiguel P, Gaut BS, Tikhonov A, Nakajima Y, Bennetzen JL: The paleontology of intergene retrotransposons of maize. Nat Genet. 1998, 20: 43-45.View ArticlePubMed
- Vitte C, Bennetzen JL: Eukaryotic Transposable Elements and Genome Evolution Special Feature: analysis of retrotransposon structural diversity uncovers properties and propensities in angiosperm genome evolution. Proc Natl Acad Sci U S A. 2006, 103: 17638-17643.PubMed CentralView ArticlePubMed
- Morgante M: Plant genome organisation and diversity: the year of the junk!. Curr Opin Biotechnol. 2006, 17: 168-173.View ArticlePubMed
- Lisch D, Bennetzen JL: Transposable element origins of epigenetic gene regulation. Curr Opin Plant Biol. 2011, 14: 156-161.View ArticlePubMed
- Bennetzen JL: Mechanisms and rates of genome expansion and contraction in flowering plants. Genetica. 2002, 115: 29-36.View ArticlePubMed
- Ma J, Bennetzen JL: Rapid recent growth and divergence of rice nuclear genomes. Proc Natl Acad Sci U S A. 2004, 101: 12404-12410.PubMed CentralView ArticlePubMed
- Devos KM, Brown JK, Bennetzen JL: Genome size reduction through illegitimate recombination counteracts genome expansion in Arabidopsis. Genome Res. 2002, 12: 1075-1079.PubMed CentralView ArticlePubMed
- Vitte C, Panaud O: Formation of solo-LTRs through unequal homologous recombination counterbalances amplifications of LTR retrotransposons in rice Oryza sativa L. Mol Biol Evol. 2003, 20: 528-540.View ArticlePubMed
- Hull R: Classifying reverse transcribing elements: a proposal and a challenge to the ICTV, International Committee on Taxonomy of Viruses. Arch Virol. 2001, 146: 2255-2261.View ArticlePubMed
- Wicker T, Sabot F, Hua-Van A, Bennetzen JL, Capy P, Chalhoub B, Flavell A, Leroy P, Morgante M, Panaud O: A unified classification system for eukaryotic transposable elements. Nat Rev Genet. 2007, 8: 973-982.View ArticlePubMed
- Pereira V: Insertion bias and purifying selection of retrotransposons in the Arabidopsis thaliana genome. Genome Biol. 2004, 5: R79-PubMed CentralView ArticlePubMed
- Zhu Y, Dai J, Fuerst PG, Voytas DF: Controlling integration specificity of a yeast retrotransposon. Proc Natl Acad Sci U S A. 2003, 100: 5891-5895.PubMed CentralView ArticlePubMed
- Lewinski MK, Yamashita M, Emerman M, Ciuffi A, Marshall H, Crawford G, Collins F, Shinn P, Leipzig J, Hannenhalli S: Retroviral DNA integration: viral and cellular determinants of target-site selection. PLoS Pathog. 2006, 2: e60-PubMed CentralView ArticlePubMed
- Marin I, Llorens C: Ty3/Gypsy retrotransposons: description of new Arabidopsis thaliana elements and evolutionary perspectives derived from comparative genomic data. Mol Biol Evol. 2000, 17: 1040-1049.View ArticlePubMed
- Gorinsek B, Gubensek F, Kordis D: Evolutionary genomics of chromoviruses in eukaryotes. Mol Biol Evol. 2004, 21: 781-798.View ArticlePubMed
- Koonin EV, Zhou S, Lucchesi JC: The chromo superfamily: new members, duplication of the chromo domain and possible role in delivering transcription regulators to chromatin. Nucleic Acids Res. 1995, 23: 4229-4233.PubMed CentralView ArticlePubMed
- Brehm A, Tufteland KR, Aasland R, Becker PB: The many colours of chromodomains. Bioessays. 2004, 26: 133-140.View ArticlePubMed
- Gao X, Hou Y, Ebina H, Levin HL, Voytas DF: Chromodomains direct integration of retrotransposons to heterochromatin. Genome Res. 2008, 18: 359-369.PubMed CentralView ArticlePubMed
- Chatterjee AG, Leem YE, Kelly FD, Levin HL: The chromodomain of Tf1 integrase promotes binding to cDNA and mediates target site selection. J Virol. 2009, 83: 2675-2685.PubMed CentralView ArticlePubMed
- Weber B, Schmidt T: Nested Ty3-gypsy retrotransposons of a single Beta procumbens centromere contain a putative chromodomain. Chromosome Res. 2009, 17: 379-396.View ArticlePubMed
- Tranzschel W: Die Arten der Gattung Beta. Bull Appl Bot and Pl Breed. 1927, 17: 203-223.
- Arumuganathan K, Earle E: Nuclear DNA content of some important plant species. Plant Mol Biol Rep. 1991, 9: 208-218.View Article
- Flavell RB, Bennett M, Smith J, Smith D: Genome size and the proportion of repeated nucleotide sequence DNA in plants. Biochem Genet. 1974, 12: 257-269.View ArticlePubMed
- Weber B, Wenke T, Frömmel U, Schmidt T, Heitkam T: The Ty1-copia families SALIRE and Cotzilla populating the Beta vulgaris genome show remarkable differences in abundance, chromosomal distribution, and age. Chromosome Res. 2010, 18: 247-263.View ArticlePubMed
- Wollrab C, Heitkam T, Holtgräwe D, Weisshaar B, Minoche AE, Dohm JC, Himmelbauer H, Schmidt T: Evolutionary reshuffling in the Errantivirus lineage Elbe within the Beta vulgaris genome. Plant J. 2012, 72: 636-651.View ArticlePubMed
- Parniske M, Wulff BBH, Bonnema G, Thomas CM, Jones DA, Jones JDG: Homologues of the Cf-9 disease resistance gene (Hcr9s) are present at multiple loci on the short arm of tomato chromosome 1. Mol Plant Microbe Interact. 1999, 12: 93-102.View ArticlePubMed
- Gindullis F, Dechyeva D, Schmidt T: Construction and characterization of a BAC library for the molecular dissection of a single wild beet centromere and sugar beet (Beta vulgaris) genome analysis. Genome. 2001, 44: 846-855.View ArticlePubMed
- Siebert PD, Chenchik A, Kellogg DE, Lukyanov KA, Lukyanov SA: An improved PCR method for walking in uncloned genomic DNA. Nucleic Acids Res. 1995, 23: 1087-1088.PubMed CentralView ArticlePubMed
- Kuykendall D, Shao J, Trimmer K: A nest of LTR retrotransposons adjacent the disease resistance-priming gene NPR1 in Beta vulgaris L. U.S. Hybrid H20. Int J Plant Genomics 2009. 2009, ID576742-
- Zakrzewski F, Wenke T, Holtgrawe D, Weisshaar B, Schmidt T: Analysis of a c0t-1 library enables the targeted identification of minisatellite and satellite families in Beta vulgaris. BMC Plant Biol. 2010, 10: 8-PubMed CentralView ArticlePubMed
- Poulter R, Butler M: A retrotransposon family from the pufferfish (fugu) Fugu rubripes. Gene. 1998, 215: 241-249.View ArticlePubMed
- Farman ML, Tosa Y, Nitta N, Leong SA: MAGGY, a retrotransposon in the genome of the rice blast fungus Magnaporthe grisea. Mol Gen Genet. 1996, 251: 665-674.PubMed
- Nakayashiki H, Matsuo H, Chuma I, Ikeda K, Betsuyaku S, Kusaba M, Tosa Y, Mayama S: Pyret, a Ty3/Gypsy retrotransposon in Magnaporthe grisea contains an extra domain between the nucleocapsid and protease domains. Nucleic Acids Res. 2001, 29: 4106-4113.PubMed CentralView ArticlePubMed
- 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: 845-855.View ArticlePubMed
- Jacobs SA, Khorasanizadeh S: Structure of HP1 chromodomain bound to a lysine 9-methylated histone H3 tail. Science. 2002, 295: 2080-2083.View ArticlePubMed
- Nielsen PR, Nietlispach D, Mott HR, Callaghan J, Bannister A, Kouzarides T, Murzin AG, Murzina NV, Laue ED: Structure of the HP1 chromodomain bound to histone H3 methylated at lysine 9. Nature. 2002, 416: 103-107.View ArticlePubMed
- Novikova O: Chromodomains and LTR retrotransposons in plants. Commun Integr Biol. 2009, 2: 158-PubMed CentralView ArticlePubMed
- Jiang J, Birchler JA, Parrott WA, Dawe RK: A molecular view of plant centromeres. Trends Plant Sci. 2003, 8: 570-575.View ArticlePubMed
- Neumann P, Navratilova A, Koblizkova A, Kejnovsky E, Hribova E, Hobza R, Widmer A, Dolezel J, Macas J: Plant centromeric retrotransposons: a structural and cytogenetic perspective. Mob DNA. 2011, 2: 4-PubMed CentralView ArticlePubMed
- 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.PubMed CentralView ArticlePubMed
- Alkan C, Sajjadian S, Eichler EE: Limitations of next-generation genome sequence assembly. Nat Methods. 2011, 8: 61-65.PubMed CentralView ArticlePubMed
- Wicker T, Keller B: Genome-wide comparative analysis of copia retrotransposons in Triticeae, rice, and Arabidopsis reveals conserved ancient evolutionary lineages and distinct dynamics of individual copia families. Genome Res. 2007, 17: 1072-1081.PubMed CentralView ArticlePubMed
- Balint-Kurti P: Identification and chromosomal localization of the monkey retrotransposon in Musa sp. Mol Gen Genet. 2000, 263: 908-915.View ArticlePubMed
- Hribova E, Neumann P, Matsumoto T, Roux N, Macas J, Dolezel J: Repetitive part of the banana (Musa acuminata) genome investigated by low-depth 454 sequencing. BMC Plant Biol. 2010, 10: 204-PubMed CentralView ArticlePubMed
- Burke WD, Calalang CC, Eickbush TH: The site-specific ribosomal insertion element type II of Bombyx mori (R2Bm) contains the coding sequence for a reverse transcriptase-like enzyme. Mol Cell Biol. 1987, 7: 2221-2230.PubMed CentralView ArticlePubMed
- Heitkam T, Schmidt T: BNR - a LINE family from Beta vulgaris - contains a RRM domain in open reading frame 1 and defines a L1 sub-clade present in diverse plant genomes. Plant J. 2009, 59: 872-882.View ArticlePubMed
- Jo SH, Koo DH, Kim JF, Hur CG, Lee S, Yang TJ, Kwon SY, Choi D: Evolution of ribosomal DNA-derived satellite repeat in tomato genome. BMC Plant Biol. 2009, 9: 42-PubMed CentralView ArticlePubMed
- Kuykendall LD, Shao JY: Genes encoding callose synthase and phytochrome a are adjacent to a MAP3Kalpha-like gene in Beta vulgaris US H20. Int J Plant Genomics. 2011, 2011: 370548-PubMed CentralView ArticlePubMed
- Kashkush K, Feldman M, Levy A: Transcriptional activation of retrotransposons alters the expression of adjacent genes in wheat. Nature Genet. 2002, 33: 102-106.View ArticlePubMed
- Slotkin RK, Martienssen R: Transposable elements and the epigenetic regulation of the genome. Nat Rev Genet. 2007, 8: 272-285.View ArticlePubMed
- Hollister JD, Smith LM, Guo YL, Ott F, Weigel D, Gaut BS: Transposable elements and small RNAs contribute to gene expression divergence between Arabidopsis thaliana and Arabidopsis lyrata. Proc Nat Acad Sci U S A. 2011, 108: 2322-2327.View Article
- Sandmeyer S: Integration by design. Proc Natl Acad Sci U S A. 2003, 100: 5586-5588.PubMed CentralView ArticlePubMed
- Lewinski MK, Bushman FD: Retroviral DNA integration – mechanism and consequences. Adv Genet. 2005, 55: 147-181.View ArticlePubMed
- Xie W, Gai X, Zhu Y, Zappulla DC, Sternglanz R, Voytas DF: Targeting of the yeast Ty5 retrotransposon to silent chromatin is mediated by interactions between integrase and Sir4p. Mol Cell Biol. 2001, 21: 6606-6614.PubMed CentralView ArticlePubMed
- Roudier F, Ahmed I, Berard C, Sarazin A, Mary-Huard T, Cortijo S, Bouyer D, Caillieux E, Duvernois-Berthet E, Al-Shikhley L: Integrative epigenomic mapping defines four main chromatin states in Arabidopsis. EMBO J. 2011, 30: 1928-1938.PubMed CentralView ArticlePubMed
- Kalderon D, Roberts BL, Richardson WD, Smith AE: A short amino acid sequence able to specify nuclear location. Cell. 1984, 39: 499-509.View ArticlePubMed
- Moore SP, Rinckel LA, Garfinkel DJ: A Ty1 integrase nuclear localization signal required for retrotransposition. Mol Cell Biol. 1998, 18: 1105-1114.PubMed CentralView ArticlePubMed
- Kenna MA, Brachmann CB, Devine SE, Boeke JD: Invading the yeast nucleus: a nuclear localization signal at the C terminus of Ty1 integrase is required for transposition in vivo. Mol Cell Biol. 1998, 18: 1115-1124.PubMed CentralView ArticlePubMed
- McLane LM, Pulliam KF, Devine SE, Corbett AH: The Ty1 integrase protein can exploit the classical nuclear protein import machinery for entry into the nucleus. Nucleic Acids Res. 2008, 36: 4317-4326.PubMed CentralView ArticlePubMed
- Nagaki K, Song J, Stupar RM, Parokonny AS, Yuan Q, Ouyang S, Liu J, Hsiao J, Jones KM, Dawe RK: Molecular and cytological analyses of large tracks of centromeric DNA reveal the structure and evolutionary dynamics of maize centromeres. Genetics. 2003, 163: 759-770.PubMed CentralPubMed
- Langdon T, Seago C, Mende M, Leggett M, Thomas H, Forster JW, Jones RN, Jenkins G: Retrotransposon evolution in diverse plant genomes. Genetics. 2000, 156: 313-325.PubMed CentralPubMed
- Dawe RK, Henikoff S: Centromeres put epigenetics in the driver’s seat. Trends Biochem Sci. 2006, 31: 662-669.View ArticlePubMed
- Topp CN, Zhong CX, Dawe RK: Centromere-encoded RNAs are integral components of the maize kinetochore. Proc Natl Acad Sci U S A. 2004, 101: 15986-15991.PubMed CentralView ArticlePubMed
- Neumann P, Yan H, Jiang J: The centromeric retrotransposons of rice are transcribed and differentially processed by RNA interference. Genetics. 2007, 176: 749-761.PubMed CentralView ArticlePubMed
- Fukai E, Umehara Y, Sato S, Endo M, Kouchi H, Hayashi M, Stougaard J, Hirochika H: Derepression of the plant chromovirus LORE1 induces germline transposition in regenerated plants. PLoS Genet. 2010, 6: e1000868-PubMed CentralView ArticlePubMed
- Tucker S, Vitins A, Pikaard C: Nucleolar dominance and ribosomal RNA gene silencing. Curr Opin Cell Biol. 2010, 22: 351-356.PubMed CentralView ArticlePubMed
- Saghai-Maroof MA, Soliman KM, Jorgensen RA, Allard RW: Ribosomal DNA spacer-length polymorphisms in barley: mendelian inheritance, chromosomal location, and population dynamics. Proc Natl Acad Sci U S A. 1984, 81: 8014-8018.PubMed CentralView ArticlePubMed
- McGrath JM, Shaw RS, de los Reyes BG, Weiland JJ: Construction of a sugar beet BAC library from a hybrid with diverse traits. Plant Mol Biol Rep. 2004, 22: 23-28.View Article
- Eddy SR: Profile hidden Markov models. Bioinformatics. 1998, 14: 755-763.View ArticlePubMed
- Llorens C, Munoz-Pomer A, Bernad L, Botella H, Moya A: Network dynamics of eukaryotic LTR retroelements beyond phylogenetic trees. Biol Direct. 2009, 4: 41-PubMed CentralView ArticlePubMed
- Llorens C, Futami R, Covelli L, Dominguez-Escriba L, Viu JM, Tamarit D, Aguilar-Rodriguez J, Vicente-Ripolles M, Fuster G, Bernet GP: The Gypsy Database (GyDB) of mobile genetic elements: release 2.0. Nucleic Acids Res. 2011, 39: 70-74.View Article
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797.PubMed CentralView ArticlePubMed
- 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.View ArticlePubMed
- Koch M, Haubold B, Mitchell-Olds T: Comparative evolutionary analysis of chalcone synthase and alcohol dehydrogenase loci in Arabidopsis, Arabis, and related genera (Brassicaceae). Mol Biol Evol. 2000, 17: 1483-1498.View ArticlePubMed
- Sambrook J, Fritsch EF, Maniatis T: Molecular cloning: A laboratory manual. 1989, Cold Spring Harbor Press: Cold Spring Harbor
- Desel C, Jung C, Cai D, Kleine M, Schmidt T: High-resolution mapping of YACs and the single-copy gene Hs1(pro-1) on Beta vulgaris chromosomes by multi-colour fluorescence in situ hybridization. Plant Mol Biol. 2001, 45: 113-122.View ArticlePubMed
- Heslop-Harrison JS: The molecular cytogenetics of plants. J Cell Sci. 1991, 100: 16-21.
- Schmidt T, Schwarzacher T, Heslop-Harrison JS: Physical mapping of rRNA genes by fluorescent in-situ hybridization and structural analysis of 5S rRNA genes and intergenic spacer sequences in sugar beet (Beta vulgaris). Theor Appl Genet. 1994, 88: 629-636.View ArticlePubMed
- Fransz PF, Alonso-Blanco C, Liharska TB, Peeters AJM, Zabel P, de Jong JH: High resolution physical mapping in Arabidopsis thaliana and tomato by fluorescence in situ hybridization to extended DNA fibres. Plant J. 1996, 9: 421-430.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.