Comparative evolution history of SINEs in Arabidopsis thaliana and Brassica oleracea: evidence for a high rate of SINE loss

Brassica oleracea and Arabidopsis thaliana belong to the Brassicaceae(Cruciferae) family and diverged 16 to 19 million years ago. Although the genome size of B. oleracea (approximately 600 million base pairs) is more than four times that of A. thaliana (approximately 130 million base pairs), their gene content is believed to be very similar with more than 85% sequence identity in the coding region. Therefore, this important difference in genome size is likely to reflect a different rate of non-coding DNA accumulation. Transposable elements (TEs) constitute a major fraction of non-coding DNA in plant species. A different rate in TE accumulation between two closely related species can result in significant genome size variations in a short evolutionary period. Short interspersed elements (SINEs) are non-autonomous retroposons that have invaded the genome of most eukaryote species. Several SINE families are present in B. oleracea and A. thaliana and we found that two of them (called RathE1 and RathE2) are present in both species. In this study, the tempo of evolution of RathE1 and RathE2 SINE families in both species was compared. We observed that most B. oleracea RathE2 SINEs are "young" (close to the consensus sequence) and abundant while elements from this family are more degenerated and much less abundant in A. thaliana. However, the situation is different for the RathE1 SINE family for which the youngest elements are found in A. thaliana. Surprisingly, no SINE was found to occupy the same (orthologous) genomic locus in both species suggesting that either these SINE families were not amplified at a significant rate in the common ancestor of the two species or that older elements were lost and only the recent (lineage-specific) insertions remain. To test this latter hypothesis, loci containing a recently inserted SINE in the A. thaliana col-0 ecotype were selected and characterized in several other A. thaliana ecotypes. In addition to the expected SINE containing allele and the pre-integrative allele (i.e. the "empty" allele), we observed in the different ecotypes, alleles with truncated portions of the SINE (up to the complete loss of the element) and of the immediate genomic flanking sequences. The absence of SINEs in orthologous positions between B. oleracea and A. thaliana and the presence in recently diverged A. thaliana ecotypes of alleles containing severely truncated SINEs suggest a very high rate of SINE loss in these species.

The evolutionary origin and genomic organization of SINEs in Arabidopsis thaliana

We have characterized the two families of SINE retroposons present in Arabidopsis thaliana. The origin, distribution, organization, and evolutionary history of RAthE1 and RAthE2 elements were studied and compared to the well-characterized SINE S1 element from Brassica. Our studies show that RAthE1, RAthE2, and S1 retroposons were generated independently from three different tRNAs. The RAthE1 and RAthE2 families are older than the S1 family and are present in all tested Cruciferae species. The evolutionary history of the RAthE1 family is unusual for SINEs. The 144 RAthE1 elements of the Arabidopsis genome cannot be classified in distinct subfamilies of different evolutionary ages as is the case for S1, RAthE2, and mammalian SINEs. Instead, most RAthE1 elements were probably derived steadily from a single source gene that was maintained intact and active for at least 12-20 Myr, a result suggesting that the RAthE1 source gene was under selection. The distribution of RAthE1 and RAthE2 elements on the Arabidopsis physical map was studied. We observed that, in contrast to other Arabidopsis transposable elements, SINEs are not concentrated in the heterochromatic regions. Instead, SINEs are grouped in the euchromatic chromosome territories several hundred kilobase pairs long. In these territories, SINE elements are closely associated with genes. A retroposition partnership between Arabidopsis SINEs and LINEs is proposed.

Analysis of the SINE S1 Pol III promoter from Brassica; impact of methylation and influence of external sequences

Transcription is an important control point in the transposable element mobilization process. To better understand the regulation of the plant SINE (Short Interspersed Elements) S1, its promoter sequence was studied using an in vitro pol III transcription system derived from tobacco cells. We show that the internal S1 promoter can be functional although upstream external sequences were found to enhance this basal level of transcription. For one putative 'master' locus (na7), three CAA triplets (in positions -12, -7 and -2) and two overlapping TATA motifs (in positions -54 to -43) were important to stimulate transcription. For this locus, two transcription initiation regions were characterized, one centered on position + 1 (first nucleotide of the S1 element) and one centered on position - 19 independently of the internal motifs. The CAA triplets only influence transcription in + 1 and work in association with the internal motifs. We show that methylation can inhibit transcription at the na7 locus. We also observe that S1 RNA is cleaved in a smaller Poly (A) minus product by a process analogous to the maturation of mammalian SINEs.

SINE retroposons can be used in vivo as nucleation centers for de novo methylation

SINEs (short interspersed elements) are an abundant class of transposable elements found in a wide variety of eukaryotes. Using the genomic sequencing technique, we observed that plant S1 SINE retroposons mainly integrate in hypomethylated DNA regions and are targeted by methylases. Methylation can then spread from the SINE into flanking genomic sequences, creating distal epigenetic modifications. This methylation spreading is vectorially directed upstream or downstream of the S1 element, suggesting that it could be facilitated when a potentially good methylatable sequence is single stranded during DNA replication, particularly when located on the lagging strand. Replication of a short methylated DNA region could thus lead to the de novo methylation of upstream or downstream adjacent sequences.

A DNA target of 30 bp is sufficient for RNA-directed DNA methylation

In higher plants, RNA-DNA interactions can trigger de novo methylation of genomic sequences via a process that is termed RNA-directed DNA methylation (RdDM). In potato spindle tuber viroid (PSTVd)-infected tobacco plants, this process can potentially lead to methylation of all C residues at symmetrical and nonsymmetrical sites within chromosomal inserts that consist of multimers of the 359-bp-long PSTVd cDNA. Using PSTVd cDNA subfragments, we found that genomic targets with as few as 30 nt of sequence complementarity to the viroid RNA are detected and methylated. Genomic sequencing analyses of genome-integrated 30- and 60-bp-long PSTVd subfragments demonstrated that de novo cytosine methylation is not limited to the canonical CpG, CpNpG sites. Sixty-base-pair-long PSTVd cDNA constructs appeared to be densely methylated in nearly all tobacco leaf cells. With the 30-bp-long PSTVd-specific construct, the proportion of cells displaying dense transgene methylation was significantly reduced, suggesting that a minimal target size of about 30 bp is necessary for RdDM. The methylation patterns observed for two different 60-bp constructs further suggested that the sequence identity of the target may influence the methylation mechanism. Finally, a link between viroid pathogenicity and PSTVd RNA-directed methylation of host sequences is proposed.

Heavy de novo methylation at symmetrical and non-symmetrical sites is a hallmark of RNA-directed DNA methylation

Previous analysis of potato spindle tuber viroid (PSTVd) RNA-infected tobacco plants has suggested that an RNA-DNA interaction could trigger de novo methylation of PSTVd transgene sequences. Using the genomic sequencing technique, the methylation pattern associated with the RNA-directed DNA methylation process has been characterized. Three different PSTVd transgene constructs all showed a similar pattern of methylation. Most of the cytosines at symmetrical as well as non-symmetrical positions appeared to be methylated in both DNA strands of the viroid sequences. Heavy methylation was mostly restricted to the viroid cDNA sequences. Flanking DNA regions immediately adjacent to the viroid cDNA displayed a lower but significant level of cytosine methylation. The observation that the heavy methylation was essentially co-extensive with the length of the PSTVd cDNA sequences provided evidence that a direct RNA-DNA interaction can act as a strong and highly specific signal for de novo DNA methylation. These data also confirmed that de novo methylation was not limited to canonical CpG and CpNpG sites, but can also involve all the cytosine residues located in the genomic region where the RNA-DNA interaction takes place.

Isolation of an RNA-directed RNA polymerase-specific cDNA clone from tomato

A 3600-bp RNA-directed RNA polymerase (RdRP)-specific cDNA comprising an open reading frame (ORF) of 1114 amino acids was isolated from tomato. The putative protein encoded by this ORF does not share homology with any characterized proteins. Antibodies that were raised against synthetic peptides whose sequences have been deduced from the ORF were shown to specifically detect the 127-kD tomato RdRP protein. The immunoresponse to the antibodies correlated with the enzymatic activity profile of the RdRP after chromatography on Q-, poly(A)-, and poly(U)-Sepharose, hydroxyapatite, and Sephadex G-200 columns. DNA gel blot analysis revealed a single copy of the RdRP gene in tomato. RdRP homologs from petunia, Arabidopsis, tobacco, and wheat were identified by using polymerase chain reaction. A sequence comparison indicated that sequences homologous to RdRP are also present in the yeast Schizosaccharomyces pombe and in the nematode Caenorhabditis elegans. The previously described induction of RdRP activity upon viroid infection is shown to be correlated with an increased steady state level of the corresponding mRNA. The possible involvement of this heretofore functionally elusive plant RNA polymerase in homology-dependent gene silencing is discussed.

Isolation of an RNA-directed RNA polymerase-specific cDNA clone from tomato

A 3600-bp RNA-directed RNA polymerase (RdRP)-specific cDNA comprising an open reading frame (ORF) of 1114 amino acids was isolated from tomato. The putative protein encoded by this ORF does not share homology with any characterized proteins. Antibodies that were raised against synthetic peptides whose sequences have been deduced from the ORF were shown to specifically detect the 127-kD tomato RdRP protein. The immunoresponse to the antibodies correlated with the enzymatic activity profile of the RdRP after chromatography on Q-, poly(A)-, and poly(U)-Sepharose, hydroxyapatite, and Sephadex G-200 columns. DNA gel blot analysis revealed a single copy of the RdRP gene in tomato. RdRP homologs from petunia, Arabidopsis, tobacco, and wheat were identified by using polymerase chain reaction. A sequence comparison indicated that sequences homologous to RdRP are also present in the yeast Schizosaccharomyces pombe and in the nematode Caenorhabditis elegans. The previously described induction of RdRP activity upon viroid infection is shown to be correlated with an increased steady state level of the corresponding mRNA. The possible involvement of this heretofore functionally elusive plant RNA polymerase in homology-dependent gene silencing is discussed.

A model for RNA-mediated gene silencing in higher plants

Homology-dependent gene silencing (HdGS) which is the generic term for transcriptional gene silencing (TGS), post-transcriptional gene silencing (PTGS) and RNA-mediated virus-resistance (RmVR) has been shown to frequently occur in transgenic plants. The role of RNA as a target and initiator of PTGS and RmVR is more and more manifested. Because TGS is assumed to be induced by a DNA-DNA interaction-mediated promoter methylation, a possible involvement of RNA in TGS was not really considered up to now. In this review we attempt to demonstrate that all three types of HdGS could be triggered by one RNA-based mechanism. A model proposing TGS as a consequence of RNA-directed DNA methylation (RdDM) and a refined mRNA threshold mechanism are presented. In contrast to the view that high amounts of mRNA are required we assume that the concentration of RNAs that can serve as efficient templates for a plant-encoded RNA-directed RNA polymerase (RdRP) plays a key role in HdGS and possibly also in natural gene regulation of non-transformed cells. According to this idea a particular information must be encoded to render mRNA turn-over products a suitable RdRP substrate. It will be discussed that such a mechanism could account for the silencing phenomena of poorly transcribed transgenes. Finally, an explanation for the coherency between PTGS and DNA methylation is documented.

DNA regions flanking the major Arabidopsis thaliana satellite are principally enriched in Athila retroelement sequences

An analysis of Arabidopsis thaliana heterochromatic regions revealed that genomic sequences immediately flanking the major 180 bp satellite are essentially made of middle repetitive sequences and that most of these sequences correspond to defective Athila retroelements. Using YAC and lambda clones, we evaluated the distribution of Athila elements in the Arabidopsis genome and showed that, despite the presence of numerous euchromatic copies, these elements are especially concentrated in or near heterochromatic regions. Sequencing of the various DNA transitions between satellite and Athila repeats provides strong evidence that most of the heterochromatic elements retrotransposed directly into 180 bp satellite clusters.

Athila, a new retroelement from Arabidopsis thaliana

An analysis of Arabidopsis thaliana heterochromatic regions allowed the identification of a new family of retroelements called Athila. These 10.5 kb elements, representing ca. 0.3% of the genome, present several features of retrotransposons and retroviruses. Athila elements are flanked by 1.5 kb long terminal repeats (LTR) that are themselves bounded by 5 bp perfect inverted repeats. These LTRs start and end with the retroviral consensus 5'TG...CA3' nucleotides. A putative tRNA-binding site and a polypurine tract are found adjacent to the 5' and 3' LTR respectively. The central domain is composed of two long open reading frames (ORFs) of 935 and 694 amino acids. Despite several indications of recent transposition activity, the translation of these ORFs failed to reveal significant homology with proteins associated to retrotransposition. We suggest that the Athila family could result from the transduction and dispersion of a cellular gene by a retrotransposon.

A method to perform direct transcutaneous intrathecal injection in rats

This work describes a method for performing direct intrathecal (i.t.) injections in rats without introducing a spinal catheter. Its ease of use in awake animals yields rapid and reproducible results with no sign of motor impairment. The quality of each injection was ensured by the observation of an injection-induced tail-flick. A 10-microL injection of methylene blue was well localized yielding a very limited diffusion along the spinal cord. The method was validated by demonstrating that morphine (i.t.) had a marked antinociceptive effect and that naloxone (i.t.) blocked the effect of systemic (s.c.) morphine in mononeuropathic rat.

An analysis of retroposition in plants based on a family of SINEs from Brassica napus

The identification of a family of SINE retroposons dispersed in the genome of oilseed rape Brassica napus has provided the basis for an evolutionary analysis of retroposition in plants. The repetitive elements (called S1Bn) are 170 bp long and occupy roughly 500 loci by haploid genome. They present characteristic features of SINE retroposons such as a 3' terminal A-rich region, two conserved polymerase III motifs (box A and B), flanking direct repeats of variable sizes, and a primary and secondary sequence homology to several tRNA species. A consensus sequence was made from the alignment of 34 members of the family. The retroposon population was divided into five subfamilies based on several correlated sets of mutations from the consensus. These precise separations in subfamilies based on "diagnostic" mutations and the random distribution of mutations observed inside each subfamily are consistent with the master sequence model proposed for the dispersion of mammalian retroposons. An independent analysis of each subfamily provides strong evidence for the coexpression of at least three subfamily master sequences (SMS). In contrast to mammalian retroposition, diagnostic positions are not shared between SMS. We therefore propose that SMS were all derived from a general master sequence (GMS) and independently activated for retroposition after a variable period of random drift. Possible models for plant retroposition are discussed.

Characterization of minisatellites in Arabidopsis thaliana with sequence similarity to the human minisatellite core sequence

A strategy based on random PCR amplification was used to isolate new repetitive elements of Arabidopsis thaliana. One of the random PCR product analyzed by this approach contained a tandem repetitive minisatellite sequence composed of 33 bp repeated units. The genomic locus corresponding to this PCR product was isolated by screening a lambda genomic library. New related loci were also isolated from the genomic library by screening with a 14 mer oligonucleotide representing a region conserved among the different repeated units. Alignment of the consensus sequence for each minisatellite locus allowed the definition of an Arabidopsis thaliana core sequence that shows strong sequence similarities with the human core sequence and with the generalized recombination signal Chi of Escherichia coli. The minisatellites were tested for their ability to detect polymorphism, and their chromosomal position was established.