- Research article
- Open Access
Bead-probe complex capture a couple of SINE and LINE family from genomes of two closely related species of East Asian cyprinid directly using magnetic separation
BMC Genomics volume 10, Article number: 83 (2009)
Short and long interspersed elements (SINEs and LINEs, respectively), two types of retroposons, are active in shaping the architecture of genomes and powerful tools for studies of phylogeny and population biology. Here we developed special protocol to apply biotin-streptavidin bead system into isolation of interspersed repeated sequences rapidly and efficiently, in which SINEs and LINEs were captured directly from digested genomic DNA by hybridization to bead-probe complex in solution instead of traditional strategy including genomic library construction and screening.
A new couple of SINEs and LINEs that shared an almost identical 3'tail was isolated and characterized in silver carp and bighead carp of two closely related species. These SINEs (34 members), designated HAmo SINE family, were little divergent in sequence and flanked by obvious TSD indicated that HAmo SINE was very young family. The copy numbers of this family was estimated to 2 × 105 and 1.7 × 105 per haploid genome by Real-Time qPCR, respectively. The LINEs, identified as the homologs of LINE2 in other fishes, had a conserved primary sequence and secondary structures of the 3'tail region that was almost identical to that of HAmo SINE. These evidences suggest that HAmo SINEs are active and amplified recently utilizing the enzymatic machinery for retroposition of HAmoL2 through the recognition of higher-order structures of the conserved 42-tail region. We analyzed the possible structures of HAmo SINE that lead to successful amplification in genome and then deduced that HAmo SINE, SmaI SINE and FokI SINE that were similar in sequence each other, were probably generated independently and created by LINE family within the same lineage of a LINE phylogeny in the genomes of different hosts.
The presented results show the advantage of the novel method for retroposons isolation and a pair of young SINE family and its partner LINE family in two carp fishes, which strengthened the hypotheses containing the slippage model for initiation of reverse transcription, retropositional parasitism of SINEs on LINEs, the formation of the stem loop structure in 3'tail region of some SINEs and LINEs and the mechanism of template switching in generating new SINE family.
SINE and LINE are interspersed nucleotide repeats distributed widely in eukaryotic genomes and occupy a substantial fraction of genome. For example, Alu and LINE1 constitute more than 13% and 20% of human genome, respectively . They proliferate and replicate themselves through a "copy and paste" mechanism called retroposition involving transcription of their genomic copies followed by reverse transcription of an RNA intermediate and resulting cDNAs reintegration at a new location into the genome host [2–5]. Therefore, SINE and LINE, together with processed pseudogenes, are classified as retrotransposons .
Many biologists are interested and make great effort to isolate and characterize many retroposons (LINE and SINE) in various organisms, because it is an indispensable step for addressing the question that how they originate and evolve as well as their functioning and impact on the evolution of eukaryotic genomes. So far, over 100 LINE and nearly 100 SINE families have been described to date in various eukaryotic genomes . Moreover, retroposons insertions have proven to be nearly perfect tools for studies of phylogeny and population biology [8, 9] and have been successfully used to resolve phylogenetic relationships among various groups of different taxonomic rank [10–15]. Especially so far SINEs appear to have gained novel functions, acting for example as enhancers or silencers that regulate the expression of preexisting functional genes [16–18].
Currently methods allowing for isolation of SINE and LINE from an unknown genome mostly depend on construction of genomic library and subsequently screening by colony hybridization method using probe specific to particular region of repeated elements. Here we propose and describe a new strategy and method used to isolate SINEs and LINEs rapidly, in which library construction and screening is completely eliminated. This method is based on hybridization capture of repetitive elements from digested genomic DNA in solution using biotinylated oligonucleotide probes which have been pre-attached to streptavidin magnetic beads. Subsequently, the captured probe-target DNA fragment complex immobilized on the magnetic beads were selected and separated from other non-complementary fragments by magnetic separation, then released and amplified by adapter polymerase chain reaction (PCR), finally the PCR products enriched for SINE and LINE were cloned directly into T-vector for sequencing. The whole procedure was completed within only about a week.
LINEs are approximately 4–7 kilobase pair (kbp) in length and encode an endonuclease (EN) and a reverse transcriptase (RT), both of which are required for LINE retrotransposition . Luan et al. proposed the "target-primed reverse transcription" (TPRT) as the mechanism of LINE retrotransposition, in which the LINE EN creates a nick in the DNA of the host genome and the RT synthesizes cDNA in situ using a 3' OH of the DNA generated by the nick as a primer. In contrast, SINEs are relatively short (about 100–500 bp) and non-autonomous retroposons without ORFs and so lack the machinery to replicate themselves. It is suggested that SINE has recruited the enzymatic machinery for retroposition from the corresponding LINE through the common "tail"sequence [19, 20], based on the observation that the sequences of many couple of SINE and LINE pairs isolated from many organisms were similar in their 3'end regions [21–24]. This scenario is supported by recent experiments of retrotransposition assay of eel UnaL2 and human LINE1. UnaL2 can strictly recognize a specific sequence at their 3'tail and mobilize transcript that has the 3' tail of UnaSINE1 [25, 26], whereas Human LINE L1 can mobilize human SINE Alu via the poly A tail (no such 3'end-specific region) at the 3'end [27, 28].
The template switch during TPRT was proposed as possible mechanism to explain the formation of chimeric retrotranscripts from a full copy of U6 small nuclear RNA fused to the 3'terminus of L1 [29–31] and the observation that several SINE families have a common 5' half sequence but different 3'tails [32, 33]. So maybe the process of how SINE acquired the tail of partner LINE also resulted from the template switch between LINE and other RNA of SINE-to-be during TPRT [21, 29, 34].
The two closely related species of East Asian cyprinids, silver carp (Hypophthalmichthys molitrix) and bighead carp (Aristichthys nobilis), are believed to have origined recently . In the present paper, we successfully use magnetic-bead based system to isolate SINE and LINE in this two species by developing special protocol. The data show that the designated HAmoSINE family was successfully proliferated recently through borrowing the enzymatic machinery of partner LINEs for retrotransposition in the two genomes. After comparison and detail analysis, we deduced that the HAmo SINE, SmaI SINE and FokI SINE that are similar in sequence with each other, were probably generated independently through the switch template between the LINE2 and RNA of SINE-to-be in respective genome because of an existence of no-similar central region between them. At last, the advantages of this new method and possibility of combining with other protocols are provided.
All species DNA was isolated from ethanol-fixed tissues (fins or muscle) by incubation with proteinase K followed by phenol/chloroform extraction .
Preliminary PCR to detect and identify SINE
AB-PCR was performed in silver carp and bighead carp as described elsewhere . The reaction mixture (100 μl) contained 10 ng of genomic DNA and two 12-nucleotide primers ("A": 5'-TRGCTCAGTGGT-3', "B":5'-GGRATYGAACYC-3') specific to A and B boxes of RNA pol III promoter consensus, respectively. After 27 PCR cycles (95°C, 1 min; 34°C, 1 min; 72°C, 30s), the amplified ~55 bp DNA fragments were isolated by electrophoresis in 5% agarose gel.
Inverse PCR was carried out in silver carp as with modification of the method described elsewhere . A pair of inverse primers: primer IF, primer IR (Figure 1, Additional file 1) was designed according to consensus sequence of AB-PCR fragments in silver carp and bighead carp. In brief, HaeIII-digested genomic DNA fragments were self-circularized in a final concentration of 5 ng/μl in a 100 μl ligation reaction. Inverse PCR (94°C, 1 min; 52°C, 1 min; 72°C, 2 min) was carried out using the 100 ng of above circularized DNA as template. The resulting smear fragments were cloned and sequenced.
A pair of primers (primer ItF, primer ItR, Figure 1) corresponding to internal region of SINE was used to detect many individual SINE copies. The PCR was run in a total volume of 20 μl including 200 ng DNA template with 25 cycles of 95°C 40s, 62°C 40s, 72°C 40s. Nine clones were selected randomly and sequenced. The clone Hmo41_It was used as probe to conduct the next retroposons enrichment strategy.
Retroposons Enrichment Strategy
A. Preparation of genomic pool
1) Digestion of genomic DNA
Approximately 40 mg of genomic DNA was completely digested with HaeIII (20 U/μl, Promega, Madison, WI, USA) overnight in a total volume of 100 μl. The fragmented DNA were subsequently separated by electrophoresis in 1% agarose. Fragments ranging from 700–2000 bp were purified from the gel using gel extraction kit (Omega) and finally suspended in 40 μl of H2O.
2) Ligation of adapters
The adapter oligoA (5 P'-GGCAGGATCCACTGAATTCGC-3') and oligoB (5'-AGCGAATTCAGTGGATCCTGCC-3') were annealed by heating an equal volume of 10 μM oligonucleotides for 3 min at 95°C, 2 min at 65°C, 2 min at 45°C, 1 min at 25°C, conserved at 4°C. The annealed product is a double-stranded linker of which one end is blunt while the other has a 3' A overhang. Additionally, the oligoA was phosphorylated at the 5' base during manufacturing. Excess of annealed linkers were ligated to above 40 μl of prepared HaeIII-fragmented DNA (approximate 4 μg) in a 100 μl reaction containing 2 μM double-stranded linkers, l× ligase buffer, 20 units T4 DNA ligase (Fermentas, MBI), 10 μl 50% PEG4000. The reaction proceeded overnight at 22°C, and then purified through Takara column and resuspend in 150 μl H2O.
3) PCR enrichment
Set up 20 PCR reactions with 28 μL ddH2O, 4.0 μL 10 mM dNTP's, 5 μL 10× PCR Buffer, 10 μL of the 2 μM oligo B primer and 0.5 μL of Taq each. Then add 1 μL of the linker ligation product to each PCR tube. The PCR reaction profile began with a 5 min 72°C filling in the nick between each linker and size fraction fragment left by the ligation step, then followed by 12 cycles of 95°C for 45s, 55°C for 45s, 72°C for 1 min 50 s. A 10 min extension step concluded the reaction. Then the total 20 tube PCR products were purified using 3–4 columns column in order not to overload the columns and finally resuspended in 150 μL H2O. The reason for doing the 16–24 separate PCR reaction at only 12 cycles is to maintain the complexity of the linker ligation mixture. Otherwise, a lot of identical clones would produce in the end. Before being used, the pool must be heat-denatured at 95°C for 10 min to make target single-strand DNA accessible to probe.
B. Preparation of bead-probe complex
1) Probe biotinylation
Plasmid Hmo41_It corresponding to an internal region (18–144 bp) of an individual SINE was used as template to be biotinylated. In order to label one biotin at one terminus of the DNA fragment, we used the primers (primer ItF, ItR) with only primer ItF biotinylated to perform PCR. At last, the double strand PCR DNA with one strand biotinylated was purified in 100 μL H2O. Before bound to beads, the double strand biotinylated probe must be denatured at 95°C for 10 min.
2) Probe bound to beads
Following the manufacturer's recommendation, 200 μL of Streptavidin Magnetic Particles (10 mg/ml, Roche, Mannheim, Germany) was collected by removing the storage buffer then washed three times with 300 μL binding buffer TEN100 (10 mM Tris-HCl, 1 mM EDTA, 100 mM NaCl, pH 7.5) for 5 min. Each time remove the supernatant using magnetic particle concentrator.
Above 100 μL denatured probe is added to beads of 300 μL binding buffer TEN100, then incubate for 30 min at room temperature to specifically bind the biotinylated strand, then removed supernatant containing non-biotinylated DNA strands followed by washing two times for 5 min with TEN100, then the preparation of the single strand probe-beads complex is accomplished.
C. Capture of target sequences
The beads were washed once with 200 μL of hybridization buffer (5 × SSC, 0.1% SDS) for 5 min, then 150 μL of buffer (10 × SSC 0.2% SDS, preheated to 55°C) and 150 μL above denatured genomic pool were added to resuspend the beads followed by hybridization in 55°C for 2 hours, then non-complementary sequences were removed by washing successively with 400 μL TEN1000 (10 mM Tris-HCl, 1 mM EDTA, 1000 mM NaCl, pH 7.5) three times for 5 min; 400 μL buffer (0.2*SSC,0.1%SDS) three time for 5 min; 400 μL TEN1000 for 10 min, all these washing are conducted at room temperature, finally target DNA were release from the beads by elution at 95°C for 5 min in 50 μL H2O.
2) Adapter PCR
Set up separate 4 PCR reactions with 14 μL ddH2O, 2.0 μL 10 mM dNTPs, 2.5 μL 10× PCR Buffer, 5 μL of the 2 μM Er1Bh1Blunt primer and 0.5 μL of Taq polymerase each. Then add 1 μL of the linker ligation product to each PCR tube. The PCR reaction profile began with a 5 min 94°C then followed by 15 cycles of 95°C for 45s, 55°C for 45s, 72°C for 1 min 50 s. A 10 min extension step concluded the reaction. Then the total 4 tube PCR products were purified for cloning.
3) Cloning and Positive detection
The enriched PCR products were ligated directly into T- vector (Takara) using T4 DNA ligase, taking advantage of the 3'A overhangs often produced by Taq polymerase. Colony PCR amplification were performed directly on many single bacterial colony to determine the size of individual inserts, then select at random and sequence many clones in which inserted fragments were longer than 700 bp.
Quantification of HAmo SINE copy number in two genomes using Quantitative Real Time-PCR
Plasmid Hmo41_It corresponding to one individual copy of HAmo SINE and Genomic DNA of silver carp and bighead carp were prepared as standard and sample for Real-Time PCR, respectively. Then, concentration of them were measured using spectrophotometer and five-fold serial dilutions of them were prepared respectively as templates to perform Real-Time PCR in a PCR machine (Bio-Rad, Chromo4) one time. All Real-Time PCR reactions was performed with 40 cycles at 95°C 40s, 62°C 40s, 72°C 40s including Primer ItF and ItR (300 nM final concentration) and SYBR GREEN in a final volume of 25 μL. At last, a melting curve analysis was done after the amplification phase. The standard curve and data analysis were carried out in the software MJ Opticon Monitor 3.1.
Characterization of HAmo LINE family
To determine the 5' upstream sequence from a breakpoint at the HaeIII site of the HAmo LINE, we employed the method of genomic DNA walking in which TAIL-PCR (thermal asymmetric interlaced PCR) was conducted using one arbitrary degenerate prime provided by kit (Takara) and special primer designed according to the consensus of HaeIII-fragmented LINEs. The whole PCR processes are conducted according to manufacture's instruction and the last PCR products were cloned and sequenced.
Preliminary PCR to detect and identify SINE
AB-PCR was amplified with very small amount of genomic DNA as a template and two oligonucleotides specific to boxes A and B of the promoter of RNA polymerase III as primers as described elsewhere . Among many AB-PCR clones, we obtained some high similar sequences in two closely related species, silver carp and bighead carp, respectively (Figure 1A). The similar but not identical AB-PCR sequences and intact A box and B box, together with its reasonable sequence similarity to certain tRNAs (see below) indicated that they may have been amplified from different SINE copies of one same SINE family in the above two species.
To amplify the regions flanking AB-PCR fragment and test whether the above AB-PCR fragments belong to a part of a certain SINEs as we expected, a pair of primer: primer IvF and IvR which face in opposite orientations and correspond to consensus sequences of AB-PCR was designed to perform inverse PCR in silver carp (Figure 1A, Additional file 1). The whole procedure including recovery of origin sequence is shown schematically and described in Additional file 1. 7 of the recovered original sequences of inverse PCR shared a high similar region possessing characteristic features of typical SINEs, including LINE2-related region followed by short tandem repeats (TAAATG), but they differed in their flanking regions which indicate they may represent different retroposon locus (Figure 1B). These significant evidences implied that the shared region may be SINE and provided us a preliminary window to see the full structure of SINE.
In order to get the probe specific to the SINE for the next new non-library retroposons enrichment method, we designed a pair of primer (primer ItF, primer ItR) corresponding to the internal region of SINE to detect many individual SINE copies (Figure 1C). Nine clones selected by random show little sequence divergence indicate that this SINE family may be a young family. At last we selected plasmid Hmo41_It to be biotinylated as probe for the non-library enrichment strategy.
Non-library Retroposons Enrichment Strategy
The Non-library Retroposons Enrichment Strategy has been conducted using plasmid Hmo41_It corresponding to the internal region of an individual SINE copy as probe to directly capture the HaeIII-fragmented genomic DNA containing the SINE sequence in solution. The whole procedure was schemed in Figure 2 and described in the Materials and Methods section. The entire isolation procedure can be completed in about a week. Every important phase can be monitored by the electrophoresis (Figure 3). At last, captured specific DNA fragments were cloned and sequenced and most of them were 700–1100 bp in length (Figure 3). In this case, the efficiency of the non-library strategy reaches nearly 60% by calculating the ratio of unique positive clones (Additional file 2). Finally, 51 and 29 SINE loci are determined in silver carp and bighead carp, respectively. Simultaneously, 28 LINE elements (HaeIII-fragmented LINEs and 5' truncated LINEs) were isolated in the final products, because the probe contained an about 40 bp tail region shared by SINE and LINE.
Identification of young SINE family in silver carp and bighead carp
Using the above non-library enrichment method, we isolated and characterized a new SINE family containing 21 and 13 members from silver carp and bighead carp, respectively. We designated it HAmo SINE family for combining Hmo (Hypophthalmichthys molitrix) and Amo (Aristichthys mobilis). The consensus sequence of HAmo SINE is 150 bp in length and identical in the two species, which has typical structure of SINEs: a tRNALys-related promoter region at their 5'-end, a unique central family-specific region and an end with LINE2-derived 3'-terminus preceding the short tandem repeats TAAATG (Figure 4, 5).
We cannot find enough diagnostic nucleotides to divide them into subfamily, and only clone 65 in silver carp and clone 599 in bighead carp shared a common A9 insertion in the tRNA-unrelated region. Almost all members in silver carp (except clone1093 with uncompleted 5' sequence) and 7 out of 13 members in bighead carp are flanked by obvious TSD (target site duplications), which are thought to have been produced during retrotransposition. The small sequence divergence among the members of HAmo SINE and SD show that this SINE family seems to be very young and proliferated very recently
The HAmo SINE 5' End is derived from tRNA
BlastN homology search revealed that the tRNA-related region of the HAmo family was most similar to tRNALys in Rabbit (83%, not counting the acceptor stem) and showed equal similarity (80%) to tRNALys in Rat, chicken, mouse, Bombyx mori, Drosophila melanogaster, respectively. When compared the predicted secondary structure between the tRNA-derived region of HAmo family and Rabbit tRNALys (Figure 6), we found a feature that they have no homology in the acceptor stem region, which also happened to tRNALys-derived SmaI family . However, the significant homology in secondary structures and the numbers of nucleotides in the stem and loop structures suggested that the tRNALys species was the most likely candidate for the origin of HAmo SINEs. The obvious conservation of secondary structure in tRNA-related regions including conserved A and B boxes of the split promoter in HAmo SINE manifested their functional importance for transcription by polymerase III.
Characterization of partner LINE family of HAmo SINE in silver carp and bighead carp
Many isolated clones contain a region only matched to 3' tail of HAmo SINE probe (plasmid Hmo41_It) in the genomes of silver carp and bighead carp when isolating HAmo SINE. After Blast search and alignment of these sequences of clones, we characterize them as 5'HaeIII -fragmented and 5' truncated LINEs, which are homologs of CR1-2_DR in zebrafish. We designated these LINEs as HAmo LINE, since possibly they encode RTase responsible for retrotransposition of HAmo SINE through recognition of the common 3'tail conserved in nucleotide sequences and secondary structure (see next).
In order to determine the 5' upstream sequence from a breakpoint at the HaeIII site of these LINEs, we employed the genomic DNA walking method to determine many sequences of clones containing 5'-truncated partial LINE (see Additional file 3). At last, a consensus sequence of 1296 bp was deduced, corresponding to partial ORF that encode RTase and 3' UTR. The predicted partial amino acid sequences encoded by HAmo LINE are 72% identical to that of CR1-2_DR and homologous to that of other LINE2 (Figure 7).
Eickbush's group divided all identified LINEs into 11 distinct clades based on an extended sequence alignment of their RT domains . Recently a novel L2 clade is well separated from the CR1 clade which is widely distributed in eukaryotic genomes such as vertebrates, echinoderms and insects [24, 42]. When HAmo LINE sequence was added into analysis, the phylogenetic tree shows that HAmo LINE was most close to zebrafish CR1-2_DR and constituted a monophyletic group with zebrafish CR1-2_DR, salmon SalL2 and eel UnaL2 in the L2 clade of LINEs (Figure 8) [43–45]. So HAmo LINE are homolog of other L2 in various distantly related species that diverged over 300 million years ago .
The common tail conserved in primary and secondary structures between HAmo SINE and HAmo LINE
HAmo SINE have an approximately 42-bp-long conserved 3'-tail which are almost identical to HAmo LINE and also high similar to other SINE family and their paternal LINE families (see next, Figure 9, Figure 10). There is only one base different in this region between HAmo SINE and HAmo LINE, suggesting that the conserved 3'tail of HAmo SINE, which is important for the process of SINE retrotransposition, is derived from HAmo LINE in the same genomes of silver carp and bighead carp.
The predicted secondary structures for the 3' tail RNA of HAmoL2 and HAmo SINE forms a secondary structure consisting of a stem and a loop (Figure 9). HAmo SINE and HAmo LINE share a hairpin region with a GGAUA loop which is thought to be a recognition domain for the LINE RT in UnaL2 [26, 46]. So the clear and significant homolog in the 3'-tail between HAmo SINE and their partner HAmo LINE in primary and secondary structures suggested the HAmo SINE may borrow the enzymatic machinery of HAmo LINE to proliferate in the same genome through the conserved 3' -tail.
Interestingly, short tandem repeats (TAAATG) of variable numbers are observed in the 3'terminus of the tail in both HAmo SINE and HAmo LINE. Most of copies have more than one repetition of the repeat in the tail, which is probably required for the slippage reaction during reverse transcription initiation .
Estimation of copy numbers of HAmo SINE
The pair of internal primer: primer ItF and primer ItR also was used to amplify the HAmo SINE sequences in genomic DNA (sample) and plasmid Hmo41_It (standard) in Real-Time PCR. The results are summarized in Additional file 4. We used a serial of diluted genomic DNA as tested samples to perform PCR reaction, the final estimations of copy numbers using different concentrated genomic DNA template are very close, suggesting that the result of experiment are stable and efficient. At last, average copy numbers of HAmo SINE in haploid genome of silver carp and bighead carp were estimated to 2.22 × 105 and 1.37 × 105, respectively. Considering the possibility of mismatch between primers with more divergent HAmo SINE sequences, so the results of qRT-PCR were minimal estimates of HAmo SINE copy numbers in the two genomes of silver carp and bighead carp.
Possible structures of HAmo SINE leading to successful proliferation
The analysis of HAmo SINE shows that it was very young and proliferated recently to estimated about 2 × 105 and 1.7 × 105 copy numbers in the haploid genome of silver carp and bighead carp, respectively. In fact, most of the SINE loci isolated in this work are species-specific or even not fixed among fish populations when we detected the presence or absence of SINE insertions using flanking primers (our group, unpublished data). So HAmo SINE are highly efficient and successfully proliferated recently in the genome and it maybe owe to its overall structure and internal structure as described below.
Firstly, HAmo SINE keep the overall secondary structure and conserved A and B box in the tRNA-related region, which ensures the RNA III recognition and transcriptional activity of SINEs. More importantly, the irregularity of the acceptor stem, as same as SmaI family, seems to help to escape recognition by tRNA-processing or RNA-modifying enzymes and therefore prevent the RNA from being cleaved by the 3'-endonuclease.
Secondly, HAmo SINE share the almost identical 3'tail with HAmo LINE2 in primary sequence and secondary structures, which keep them to well utilize the LINE2 enzymatic machinery. Their shared same stem-loop region is thought to function as a recognition site for the UnaL2 protein (UnaL2p) when this region is transcribed in the RNA .
Moreover, more than one repetition of the short tandem repeat TAAATG appeared in most copies, which are revealed to be necessary for successful retrotransposition by mutational analyses in the experiments on other LINEs of the L2 clade and the initiation of reverse transcription of UnaL2 RNA in UnaL2 [25, 28].
Thirdly, RNA structure of the HAmo SINE is obviously composed of three parts: tRNA-related region, a family-specific region and LINE2-related region (Figure 10), that correspond to there parts of secondary structure of its RNA: the cloverleaf structure (the 5'domain), an unstructured region, the extended stem-loop (the 3'domain). This characteristic domain composition is experimentally probed in salmon SmaI SINE RNA and seems to have guaranteed successful and continuous amplification of SINEs in eukaryotic genomes during evolution. From this view, HAmo SINE may reveal some internal structures of SINE that may lead to its successful retrotransposition and proliferation.
Three similar but independently derived SINE families
When using HAmo SINE as query to tBlastN search, it shows unexpected high similarity to other two SINE families in salmon: SmaI family (77%) and FokI family (71%), both are young SINE families and have a limited distribution in several specific species belonging to the family Salmonidae . The two SINE families shared a common tail and parasitized SalL2 in salmon genome . After detailed comparison of the consensus sequences of them and their respective partner LINE families (Figure 10)[48, 49], we found that HAmo SINE is similar in tRNA-related region (1–76 bp) and LINE-related region (107–150 bp) with SmaI family and FokI family. But the existence of a central region (76–107 bp) which showed no similarity with each other and are specific to each family make us deduce that they are probably independently generated and evolved in respective evolutionary lineage other than horizontal transfer.
As noted in the Introduction section, the template switch during TPRT was proposed to explain how SINE acquired the tail from corresponding LINE. In this process, a short cDNA would first be generated by copying the 3'terminal LINE RNA sequence, and then RT landing pad will jump to another RNA parent of the SINE-to-be carrying an internal pol III promoter . So the above-mentioned tRNALys derived SINE may be born through template switch between respective LINE and ancestor RNA of SINE-to-be containing tRNA-derived region and family-specific region in respective genome of three fishes. Coincidently, the three young families are all derived from tRNALys or structurally related to tRNALys (Figure 6). Moreover, their parental LINE (HAmoL2 and Sal L2) of the above mentioned three SINE families are homologous and share a common tail.
In fact, tRNALys is the most common source of SINEs [7, 24, 51]. The possible reason is that maybe the ancestor tRNALys SINE RNA had special selective advantage in the above generation process or been preferentially transcripted and retrotransposed after generation among the population of RNA of SINE -to be.
So this finding suggested that the three similar but distantly related young SINE families were generated independently and created by LINE families within the same lineage of a LINE phylogeny in the genomes of different hosts.
Some aspects about the new retroposons enrichment strategy
Magnetic Bead-based isolation system has been widely used for the separation of several specific targets like cells, proteins, microsatellites and so on. However, our work is the first report about application of this system for isolation of SINEs and LINEs from fish genomes by developing new special protocol. The results demonstrate that this protocol is technically straightforward and permits the isolation of a large number of SINE and LINE from unknown genome in less time consumption and less cost and effort than is required to execute traditional protocol involving rounds of filter hybridization.
In general, if all steps work, the procedure takes only about a week from tissue to several hundred positive clones. Additionally, the purchase of the reagents needed for building and screening one library by traditional protocol will supply sufficient reagents for ten or more libraries applied by enrichment protocol. Moreover, the protocol can be easily controlled and handled since it requires little specialized equipment platform or technical expertise, May be the PCR and cloning be the most difficult step.
Our method, relying on solution hybridization, could greatly facilitate and speed up the interaction between probe and target DNA and result in better hybridization efficiencies in comparison with fixed solid supports [52, 53]. Moreover, this method can be useful in the case of low copy number SINEs and LINEs since at last only a population of sequences enriched for specific retroposons is cloned. Generally the frequency of positive clones can reach 50–90% if conditions were optimized . So it shows great advantage when usually a great number of retroposon insertions need be isolated as temporal landmarks of evolution for estimations of phylogeny.
Most steps in the protocol presented here can be readily modified to suit different experimental backgrounds and knowledge about SINE and LINE and can easily combine with other protocols. Okada's group successfully isolated many SINE families from many organisms by using the in vitro transcript of total genomic DNA as the probes utilizing the properties that SINEs are redundant in the genome and transcribed by RNA polymerase III [55, 56]. While Kramerov's group prefered to use AB-PCR product containing a 30–40 bp sequence located between boxes A and B of SINE as a probe [37, 57]. All these specific probes including known SINE sequence (this paper) can be biotinylated to join into this enrichment strategy.
But it is noted that there are many principles that should be kept in mind. Firstly, correct restriction enzyme should be selected to generate appropriate size fragments evenly and its recognized sites should not exist in the targeted repeat elements. In our work, although we isolated SINEs and LINEs simultaneously at one isolation reaction, we only obtained the HaeIII-fragmented partial LINEs because of the existed HaeIII site in full-size LINEs. Secondly, the amplification cycles of step PCR enrichment and adapter PCR (see Methods) should be optimized to generate a smear of the PCR products without specific bands. In this case, 15 and 12 of cycles were done in the two steps respectively to keep the complexity of DNA molecules for preventing the generation of a lot of identical clones at last. Moreover, the selectivity and specificity can be adjusted by altering specific probe and the stringency conditions (temperature and salinity of washing buffer). Thirdly, there are several methods for labeling one biotin at one terminus of the DNA fragment such as PCR method with one of the two primers biotinylated (this paper), end-labeling using terminal transferase, ligation reaction with a biotinylated adaptor  or direct generated by company service. No matter which method to be used, it is important to label only one biotin molecule at one terminus of the DNA fragment, otherwise magnetic beads will crosslink and clot through DNA bridges, which may result in poor reaction kinetics between beads and target molecules. In addition, isolation of large size of DNA fragments may be limited by beads binding ability and cloning efficiency of large fragments into T-vector. However, our procedures mainly base on PCR and hence could be use to track the progress of the entire process from step to step by gel electrophoresis.
The young HAmo SINE family and its partner HAmo LINE2 family shared a common 3'tail region in two carp fishes, indicated the retropositional parasitism of HAmo SINEs on HAmo LINEs and strengthened already proposed hypotheses including the formation of the stem loop structure in 3'tail region of some SINEs and LINEs and the mechanism of template switching in generating new SINE family. The finding of repeat sequences in the 3' tail of the HAmo SINEs strengthened the slippage model for initiation of reverse transcription. The obtained results show that the developed new protocol for isolation of SINE and LINE are advantages and technically straightforward. The characterization of new SINE and LINE pair is also beneficial to the future study about the molecular systematics of cyprinid fish.
short interspersed repetitive elements
long interspersed repetitive elements
Quantitative Real Time-PCR
target site duplication
Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, Devon K, Dewar K, Doyle M, FitzHugh W, et al: Initial sequencing and analysis of the human genome. Nature. 2001, 409 (6822): 860-921. 10.1038/35057062.
Singer MF: SINEs and LINEs: highly repeated short and long interspersed sequences in mammalian genomes. Cell. 1982, 28 (3): 433-434. 10.1016/0092-8674(82)90194-5.
Rogers J: Origins of repeated DNA. Nature. 1985, 317 (6040): 765-766. 10.1038/317765a0.
Weiner AM, Deininger PL, Efstratiadis A: Nonviral retroposons: genes, pseudogenes, and transposable elements generated by the reverse flow of genetic information. Annual review of biochemistry. 1986, 55: 631-661. 10.1146/annurev.bi.55.070186.003215.
Brosius J: Retroposons – seeds of evolution. Science. 1991, 251 (4995): 753-10.1126/science.1990437.
Nikaido M, Okada N: CetSINEs and AREs are not SINEs but are parts of cetartiodactyl L1. Mamm Genome. 2000, 11 (12): 1123-1126. 10.1007/s003350010221.
Kramerov DA, Vassetzky NS: Short retroposons in eukaryotic genomes. Int Rev Cytol. 2005, 247: 165-221. 10.1016/S0074-7696(05)47004-7.
Ray DA, Xing J, Salem AH, Batzer MA: SINEs of a nearly perfect character. Systematic biology. 2006, 55 (6): 928-935. 10.1080/10635150600865419.
Shedlock AM, Okada N: SINE insertions: powerful tools for molecular systematics. Bioessays. 2000, 22 (2): 148-160. 10.1002/(SICI)1521-1878(200002)22:2<148::AID-BIES6>3.0.CO;2-Z.
Sasaki T, Takahashi K, Nikaido M, Miura S, Yasukawa Y, Okada N: First application of the SINE (short interspersed repetitive element) method to infer phylogenetic relationships in reptiles: an example from the turtle superfamily Testudinoidea. Mol Biol Evol. 2004, 21 (4): 705-715. 10.1093/molbev/msh069.
Stoneking M, Fontius JJ, Clifford SL, Soodyall H, Arcot SS, Saha N, Jenkins T, Tahir MA, Deininger PL, Batzer MA: Alu insertion polymorphisms and human evolution: evidence for a larger population size in Africa. Genome research. 1997, 7 (11): 1061-1071.
Serdobova IM, Kramerov DA: Short retroposons of the B2 superfamily: evolution and application for the study of rodent phylogeny. Journal of molecular evolution. 1998, 46 (2): 202-214. 10.1007/PL00006295.
Nikaido M, Nishihara H, Hukumoto Y, Okada N: Ancient SINEs from African endemic mammals. Mol Biol Evol. 2003, 20 (4): 522-527. 10.1093/molbev/msg052.
Takahashi K, Terai Y, Nishida M, Okada N: A novel family of short interspersed repetitive elements (SINEs) from cichlids: the patterns of insertion of SINEs at orthologous loci support the proposed monophyly of four major groups of cichlid fishes in Lake Tanganyika. Mol Biol Evol. 1998, 15 (4): 391-407.
Schmitz J, Ohme M, Suryobroto B, Zischler H: The colugo (Cynocephalus variegatus, Dermoptera): the primates' gliding sister?. Mol Biol Evol. 2002, 19 (12): 2308-2312.
Sasaki T, Nishihara H, Hirakawa M, Fujimura K, Tanaka M, Kokubo N, Kimura-Yoshida C, Matsuo I, Sumiyama K, Saitou N, et al: Possible involvement of SINEs in mammalian-specific brain formation. Proc Natl Acad Sci USA. 2008, 105 (11): 4220-4225. 10.1073/pnas.0709398105.
Tomilin NV: Regulation of mammalian gene expression by retroelements and non-coding tandem repeats. Bioessays. 2008, 30 (4): 338-348. 10.1002/bies.20741.
Mariner PD, Walters RD, Espinoza CA, Drullinger LF, Wagner SD, Kugel JF, Goodrich JA: Human Alu RNA is a modular transacting repressor of mRNA transcription during heat shock. Molecular cell. 2008, 29 (4): 499-509. 10.1016/j.molcel.2007.12.013.
Eickbush TH: Transposing without ends: the non-LTR retrotransposable elements. The New biologist. 1992, 4 (5): 430-440.
Luan DD, Korman MH, Jakubczak JL, Eickbush TH: Reverse transcription of R2Bm RNA is primed by a nick at the chromosomal target site: a mechanism for non-LTR retrotransposition. Cell. 1993, 72 (4): 595-605. 10.1016/0092-8674(93)90078-5.
Ohshima K, Hamada M, Terai Y, Okada N: The 3' ends of tRNA-derived short interspersed repetitive elements are derived from the 3' ends of long interspersed repetitive elements. Mol Cell Biol. 1996, 16 (7): 3756-3764.
Ogiwara I, Miya M, Ohshima K, Okada N: Retropositional parasitism of SINEs on LINEs: identification of SINEs and LINEs in elasmobranchs. Mol Biol Evol. 1999, 16 (9): 1238-1250.
Okada N, Hamada M, Ogiwara I, Ohshima K: SINEs and LINEs share common 3' sequences: a review. Gene. 1997, 205 (1–2): 229-243. 10.1016/S0378-1119(97)00409-5.
Ohshima K, Okada N: SINEs and LINEs: symbionts of eukaryotic genomes with a common tail. Cytogenet Genome Res. 2005, 110 (1–4): 475-490. 10.1159/000084981.
Kajikawa M, Okada N: LINEs mobilize SINEs in the eel through a shared 3' sequence. Cell. 2002, 111 (3): 433-444. 10.1016/S0092-8674(02)01041-3.
Kajikawa M, Ichiyanagi K, Tanaka N, Okada N: Isolation and characterization of active LINE and SINEs from the eel. Mol Biol Evol. 2005, 22 (3): 673-682. 10.1093/molbev/msi054.
Moran JV, Holmes SE, Naas TP, DeBerardinis RJ, Boeke JD, Kazazian HH: High frequency retrotransposition in cultured mammalian cells. Cell. 1996, 87 (5): 917-927. 10.1016/S0092-8674(00)81998-4.
Dewannieux M, Esnault C, Heidmann T: LINE-mediated retrotransposition of marked Alu sequences. Nature genetics. 2003, 35 (1): 41-48. 10.1038/ng1223.
Bibillo A, Eickbush TH: The reverse transcriptase of the R2 non-LTR retrotransposon: continuous synthesis of cDNA on non-continuous RNA templates. Journal of molecular biology. 2002, 316 (3): 459-473. 10.1006/jmbi.2001.5369.
Buzdin A, Ustyugova S, Gogvadze E, Vinogradova T, Lebedev Y, Sverdlov E: A new family of chimeric retrotranscripts formed by a full copy of U6 small nuclear RNA fused to the 3' terminus of l1. Genomics. 2002, 80 (4): 402-406. 10.1006/geno.2002.6843.
Buzdin A, Gogvadze E, Kovalskaya E, Volchkov P, Ustyugova S, Illarionova A, Fushan A, Vinogradova T, Sverdlov E: The human genome contains many types of chimeric retrogenes generated through in vivo RNA recombination. Nucleic acids research. 2003, 31 (15): 4385-4390. 10.1093/nar/gkg496.
Kido Y, Himberg M, Takasaki N, Okada N: Amplification of distinct subfamilies of short interspersed elements during evolution of the Salmonidae. J Mol Biol. 1994, 241 (5): 633-644. 10.1006/jmbi.1994.1540.
Ogiwara I, Miya M, Ohshima K, Okada N: V-SINEs: a new superfamily of vertebrate SINEs that are widespread in vertebrate genomes and retain a strongly conserved segment within each repetitive unit. Genome Res. 2002, 12 (2): 316-324. 10.1101/gr.212302.
Weiner AM: SINEs and LINEs: the art of biting the hand that feeds you. Current opinion in cell biology. 2002, 14 (3): 343-350. 10.1016/S0955-0674(02)00338-1.
Wang X, Li J, He S: Molecular evidence for the monophyly of East Asian groups of Cyprinidae (Teleostei: Cypriniformes) derived from the nuclear recombination activating gene 2 sequences. Molecular phylogenetics and evolution. 2007, 42 (1): 157-170. 10.1016/j.ympev.2006.06.014.
Maniatis T, Fritsch EF, Sambrook J: Molecular Cloning: A Laboratory Manual. 1982, NY: Cold Spring Harbor Laboratory Press, Cold Spring Harbor
Borodulina OR, Kramerov DA: Wide distribution of short interspersed elements among eukaryotic genomes. FEBS Lett. 1999, 457 (3): 409-413. 10.1016/S0014-5793(99)01059-5.
Ochman H, Gerber AS, Hartl DL: Genetic applications of an inverse polymerase chain reaction. Genetics. 1988, 120 (3): 621-623.
Raba M, Limburg K, Burghagen M, Katze JR, Simsek M, Heckman JE, Rajbhandary UL, Gross HJ: Nucleotide sequence of three isoaccepting lysine tRNAs from rabbit liver and SV40-transformed mouse fibroblasts. European journal of biochemistry/FEBS. 1979, 97 (1): 305-318. 10.1111/j.1432-1033.1979.tb13115.x.
Matsumoto K, Murakami K, Okada N: Gene for lysine tRNA1 may be a progenitor of the highly repetitive and transcribable sequences present in the salmon genome. Proc Natl Acad Sci USA. 1986, 83 (10): 3156-3160. 10.1073/pnas.83.10.3156.
Malik HS, Burke WD, Eickbush TH: The age and evolution of non-LTR retrotransposable elements. Mol Biol Evol. 1999, 16 (6): 793-805.
Lovsin N, Gubensek F, Kordi D: Evolutionary dynamics in a novel L2 clade of non-LTR retrotransposons in Deuterostomia. Mol Biol Evol. 2001, 18 (12): 2213-2224.
Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J: Repbase Update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005, 110 (1–4): 462-467. 10.1159/000084979.
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599. 10.1093/molbev/msm092.
Terai Y, Takahashi K, Okada N: SINE cousins: the 3'-end tails of the two oldest and distantly related families of SINEs are descended from the 3' ends of LINEs with the same genealogical origin. Mol Biol Evol. 1998, 15 (11): 1460-1471.
Baba S, Kajikawa M, Okada N, Kawai G: Solution structure of an RNA stem-loop derived from the 3' conserved region of eel LINE UnaL2. Rna. 2004, 10 (9): 1380-1387. 10.1261/rna.7460104.
Kido Y, Aono M, Yamaki T, Matsumoto K, Murata S, Saneyoshi M, Okada N: Shaping and reshaping of salmonid genomes by amplification of tRNA-derived retroposons during evolution. Proc Natl Acad Sci USA. 1991, 88 (6): 2326-2330. 10.1073/pnas.88.6.2326.
Hamada M, Takasaki N, Reist JD, DeCicco AL, Goto A, Okada N: Detection of the ongoing sorting of ancestrally polymorphic SINEs toward fixation or loss in populations of two species of charr during speciation. Genetics. 1998, 150 (1): 301-311.
Takasaki N, Yamaki T, Hamada M, Park L, Okada N: The salmon SmaI family of short interspersed repetitive elements (SINEs): interspecific and intraspecific variation of the insertion of SINEs in the genomes of chum and pink salmon. Genetics. 1997, 146 (1): 369-380.
Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic acids research. 1997, 25 (5): 955-964. 10.1093/nar/25.5.955.
Okada N: SINEs. Current opinion in genetics & development. 1991, 1 (4): 498-504. 10.1016/S0959-437X(05)80198-4.
Lund V, Schmid R, Rickwood D, Hornes E: Assessment of methods for covalent binding of nucleic acids to magnetic beads, Dynabeads, and the characteristics of the bound nucleic acids in hybridization reactions. Nucleic acids research. 1988, 16 (22): 10861-10880. 10.1093/nar/16.22.10861.
Zammatteo N, Alexandre I, Ernest I, Le L, Brancart F, Remacle J: Comparison between microwell and bead supports for the detection of human cytomegalovirus amplicons by sandwich hybridization. Analytical biochemistry. 1997, 253 (2): 180-189. 10.1006/abio.1997.2352.
Zane L, Bargelloni L, Patarnello T: Strategies for microsatellite isolation: a review. Molecular ecology. 2002, 11 (1): 1-16. 10.1046/j.0962-1083.2001.01418.x.
Endoh H, Okada N: Total DNA transcription in vitro: a procedure to detect highly repetitive and transcribable sequences with tRNA-like structures. Proc Natl Acad Sci USA. 1986, 83 (2): 251-255. 10.1073/pnas.83.2.251.
Okada N, Shedlock AM, Nikaido M: Retroposon mapping in molecular systematics. Methods Mol Biol. 2004, 260: 189-226.
Borodulina OR, Kramerov DA: PCR-based approach to SINE isolation: simple and complex SINEs. Gene. 2005, 349: 197-205. 10.1016/j.gene.2004.12.035.
Xu H, Zhang S, Liu D, Liang CC: End-labeling of long DNA fragments with biotin and detection of DNA immobilized on magnetic beads. Molecular biotechnology. 2001, 17 (2): 183-185. 10.1385/MB:17:2:183.
This work was supported by National Natural Science Foundation of China [grant number 30530120, 30225008]. The authors thank W. Tao for help in experiments. We thank Dr. X. Wang, X. Ku of our Institute for critical reading of the manuscript. We also thank Prof. L. Li of Northwestern University, Chicago, for useful discussions.
SH and CT conceived and designed the experiments, CT and BG performed the experiments and analyzed the data, CT and SH wrote the paper. All authors have read and approved the final manuscript.
Electronic supplementary material
Additional file 3: A schematic representation of clones used to determine the consensus sequences of HAmo LINE. The clones obtained by the retroposons enrichment strategy and by genome walking method are separately located on the right and left. The GenBank accession numbers of HAmo LINE are as follows: FJ171663-FJ171689. (PPT 42 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.