Transposases exist in both prokaryotes and eukaryotes and catalyze the movement of defined DNA elements (transposon) to another part of the genome in a ‘cut and paste’ mechanism (Kleckner, 1981; Finnegan, 1989; Curcio and Derbyshire, 2003). Taking advantage of this catalytic activity, transposases are widely used in many biomedical applications: for instance, an engineered, hyperactive Tn5 transposase from E. coli can bind to synthetic 19 bp mosaic end-recognition sequences appended to Illumina sequencing adapters (termed ‘Tn5 transposome’) (Adey et al., 2010) and has been utilized in an in vitro double-stranded DNA (dsDNA) tagmentation reaction (namely simultaneously fragment and tag a target sequence with sequencing adaptors) to achieve rapid and low-input library construction for next-generation sequencing (Adey et al., 2010; Goryshin and Reznikoff, 1998; Picelli et al., 2014a; Caruccio, 2011; Ramsköld et al., 2012; Gertz et al., 2012). In addition, Tn5 was also used for in vivo transposition of native chromatin to profile open chromatin, DNA-binding proteins and nucleosome position (‘ATAC-seq’) (Buenrostro et al., 2013). While Tn5 has been broadly adopted in high-throughput sequencing, bioinformatic analysis and structural studies reveal that it belongs to the retroviral integrase superfamily that act on not only dsDNA but also RNA/DNA hybrids (for instance, RNase H). Despite the distinct substrates, these proteins all share a conserved catalytic RNase H-like domain (Figure 1a; Yang and Steitz, 1995; Savilahti et al., 1995; Nowotny, 2009; Rice and Baker, 2001). Given their structural and mechanistic similarity, we attempted to ask whether or not Tn5 is able to catalyze co-tagmentation reactions to both the RNA and DNA strands of RNA/DNA hybrids (Figure 1b), in addition to its canonical function of dsDNA transposition. In this study, we tested this hypothesis and found that indeed Tn5 possesses in vitro tagmentation activity towards both strands of RNA/DNA hybrids. As a proof of concept, we apply such Transposase-assisted RNA/DNA hybrids Co-tagmEntation (TRACE-seq) to achieve rapid and low-cost RNA sequencing starting from total RNA extracted from 10,000 to 100 cells. We find that TRACE-seq performs well when compared with conventional RNA-seq methods in terms of detected gene number, gene expression measurement, library complexity, GC content and differential expression analysis, although TRACE-seq shows bias in gene body coverage and is not strand-specific. At the same time, it avoids many laborious and time-consuming steps in traditional RNA-seq experiments. Such Tn5-assisted tagmentation of RNA/DNA hybrids could have broad applications in RNA biology and chromatin research.

Figure 1 with 1 supplement see all

Download asset Open asset

Figure 1—source data 1

qPCR Ct values of tagmentation products of samples under different conditions in Figure 1d and e.

https://cdn.elifesciences.org/articles/54919/elife-54919-fig1-data1-v2.xlsx

Download elife-54919-fig1-data1-v2.xlsx

To test whether Tn5 transposase has tagmentation activity on RNA/DNA hybrids, we prepared RNA/DNA duplexes by performing mRNA reverse transcription. We first validated the efficiency of reverse transcription and the presence of RNA/DNA duplexes using a model mRNA sequence (IRF9,~1000 nt) as template (Figure 1—figure supplement 1a). We then subjected the prepared RNA/DNA hybrids from 293T mRNA to Tn5 transposome, heat-inactivated Tn5 transposome and a blank control (without Tn5), respectively (see Methods). The hybrids were then recovered and their length distribution was analyzed by Fragment Analyzer (Figure 1c). Comparing with the heat-inactivated Tn5 sample or the blank control sample, the Tn5 transposome sample exhibited a modest but clear smear signal corresponding to small fragments ranging from ~30–650 base-pair (bp) (the blue patches in Figure 1c). Consistent with the fragmentation event, we also observed a down shift of large fragments ranging from ~700–4000 bp (the orange patches in Figure 1c). In addition, the fragmentation efficiency increased in a dose-dependent manner with the transposome, suggesting that fragmentation of RNA/DNA hybrids is dependent on Tn5 (Figure 1—figure supplement 1b).

We next asked whether RNA/DNA hybrids are tagged by Tn5 and performed quantitative polymerase chain reaction (qPCR) quantification for the three samples. We observed that cycle threshold (Ct) value of the Tn5 transposome sample is about eight cycles smaller than the heat inactivated Tn5 sample or the control sample, indicating approximately 256 times more amplifiable products (Figure 1d). We also tested different buffer conditions and found that the performance of Tn5 remained similar, indicating the robustness of the Tn5 tagmentation activity (Figure 1—figure supplement 1c). Using Sanger sequencing, we validated that the adaptor sequences are indeed ligated to the insert sequences (Figure 1—figure supplement 1d).

To compare Tn5 tagmentation efficiency between RNA/DNA hybrids and dsDNA, we performed tagmentation and qPCR on equal amount of mRNA RT products and genomic DNA (gDNA). Average Ct value of the hybrids samples was about four cycles more than gDNA samples, indicating the efficiency of Tn5 toward hybrids is about 1/16 of that of dsDNA (Figure 1e). It is known that natural RNA/DNA hybrids favor A-form conformation. Interestingly, in the presence of PEG200, hybrids were found to favor B-form conformation (Pramanik et al., 2011), which we expected to make the hybrids a better substrate of Tn5. Indeed, addition of PEG200 diminishes this difference by greatly improving the Tn5 tagmentation efficiency towards hybrids (Figure 1e). This result indicates that the conformation of substrates certainly affects the preference of Tn5. We also ruled out the possibility of gDNA contamination in RT products (Figure 1—figure supplement 1e). Taken together, under optimized conditions, Tn5 shows significantly improved efficiency towards RNA/DNA hybrids.

As reverse transcriptase could produce dsDNA from RNA/DNA hybrids, we next designed an experiment by eliminating the RT component and directly assess tagmentation activity using annealed RNA/DNA hybrids where no dsDNA is possible. We annealed in vitro transcribed and purified ssRNA (CLuc, 150 nt, GC% = 51%) with chemically synthesized complementary ssDNA. We confirmed the successful production and purity of the RNA/DNA hybrids by dot-blot assay and native-PAGE (Figure 1—figure supplement 1f and g). A 8-cycle difference between Tn5 transposome sample (Ct = 22.68) and the control sample (Ct = 30.40) was reproducibly observed (Figure 1—figure supplement 1h). As a positive control, we also annealed two complementary ssDNA strands of the same 150bp-CLuc sequence to produce dsDNA and observed that a Ct value of 18.08 for the dsDNA sample (Figure 1—figure supplement 1h). While it is unclear this difference in tagmentation efficiency obtained from short oligos can be applied to long oligos, this result clearly demonstrates that Tn5 has a direct tagmentation activity towards RNA/DNA hybrids.

Having demonstrated the tagmentation activity of Tn5 on RNA/DNA hybrids, we then thought about its potential applications. RNA/DNA duplexes can be found in many in vivo scenarios, including but not limited to R-loop and chromatin-bound lncRNAs (Santos-Pereira and Aguilera, 2015; Li and Fu, 2019). Under in vitro conditions, RNA/DNA hybrids are also key intermediates in various molecular biology and genomics experiments. For instance, RNA has to be first reverse transcribed into cDNA in a traditional RNA-seq experiment so as to construct a library for sequencing. Because traditional RNA-seq library construction involves many laborious and time-consuming steps, including mRNA purification, fragmentation, reverse transcription, second-strand synthesis, end-repair and adaptor ligation, we attempted to replace the process using the tagmentation activity towards RNA/DNA duplexes. With the help of TRACE-seq, these steps are replaced with a ‘one-tube’ protocol (Figure 2a), which uses total RNA as input material and involves just three seamless steps (reverse transcription, tagmentation and strand extension and PCR), without the need for a second strand synthesis step. We first conducted TRACE-seq with 200 ng total RNA as input and tested several enzymes and conditions (Supplementary file 1); we observed very high correlation in gene-expression levels among three replicates, indicating TRACE-seq is highly reproducible (Figure 2b). To test the robustness of TRACE-seq, we performed the experiments with 20 ng and 2 ng total RNA. TRACE-seq results are again highly reproducible among replicates (Figure 2—figure supplement 1a,b). More importantly, gene expression levels measured using different amount of starting materials remain consistent with each other (Figure 2c).

Figure 2 with 1 supplement see all

Download asset Open asset

Figure 2—source data 1

Distribution of reads across known genome features for NEBNext Ultra II RNA kit, Smart-seq2 and TRACE-seq with different amount of RNA as input.

https://cdn.elifesciences.org/articles/54919/elife-54919-fig2-data1-v2.xls

Download elife-54919-fig2-data1-v2.xls

We then compared the library quality between TRACE-seq and NEBNext Ultra II RNA library prep kit, a commonly used kit for RNA-seq library construction. In addition, we conducted a comparison to Smart-seq2, which is a similar method in its use of oligo(dT) primed cDNA synthesis and Tn5 tagmentation. The HEK293T RNA used in these libraries was all from the same batch of cells. We found that TRACE-seq libraries exhibited similar percentage of reads mapped to annotated transcripts, rRNA contamination and gene numbers to NEBNext data when mRNA was used as input, but a higher rRNA contamination when total RNA was used as input (~9%, Supplementary file 2). In addition, a similar percentage of rRNA contamination was also observed in Smart-seq2 libraries (Supplementary file 3). Most of the genes detected by TRACE-seq overlap with that of NEBNext and Smart-seq2 (Figure 2d). In addition, TRACE-seq showed comparable performance in terms of gene expression measurement, using either a set of housekeeping genes (Figure 2e) or all genes (Figure 2—figure supplement 1c). The insert size of TRACE-seq library was moderately shorter (Figure 2—figure supplement 1d); in the meanwhile, TRACE-seq shows a higher coefficient of variation of gene coverage (0.54–0.70 vs 0.42–0.44, Figure 2—figure supplement 1e). We further found that TRACE-seq showed a slight tendency to 3’ end of the gene body (Figure 2f, Figure 2—figure supplement 1f). When transcripts were grouped according to annotated lengths, we found comparable gene body coverage for transcripts shorter than 1 kb among TRACE-seq, NEBNext kit and Smart-seq2 libraries. For transcripts with length between 1 and 4 kb, a slight 3’ end bias was observed in TRACE-seq library, while for transcripts longer than 4 kb, the central regions of transcripts were less covered by both TRACE-seq and Smart-seq2. We also performed TRACE-seq by using rRNA depletion together with random-primed cDNA synthesis. While this solved the 3’ end bias, a 5’ end bias appeared, which is a common phenomenon when using random primers during reverse transcription (Figure 2—figure supplement 1f). In spite of the gene body coverage bias, GC content (Figure 2g) and library complexity (Figure 2—figure supplement 1g) are unnoticeably affected. In addition, the gene expression measurement (Figure 2e) is also unaffected here because of the use of RNA with high quality (RIN: 9.5, Figure 2—figure supplement 1h); yet, cautions might be taken when the quality of RNA is compromised. Further inspection of reads distribution of TRACE-seq over genome features revealed similar pattern for that of NEBNext and Smart-seq2 (Figure 2h). Coverages of some representative transcripts are shown in Figure 2i.

One of the most important goals of RNA-seq is to detect differentially expressed genes among different samples. Having assessed the library quality of TRACE-seq, we next compared the performance of TRACE-seq in detecting differentially expressed genes between undifferentiated and differentiated mESCs to NEBNext Ultra II RNA library prep kit. As shown in Figure 3a, TRACE-seq successfully detected 4577 differentially expressed genes (3264 up-regulated genes and 1313 down-regulated genes), while NEBNext detected 4452 differentially expressed genes (3157 up-regulated genes and 1295 down-regulated genes). The overlapping gene number is 4,071, showing very high consistency between methods (Figure 3b). Besides, the fold change of the 4071 overlapping genes is highly correlated between the two methods (R > 0.99, Figure 3c). Therefore, TRACE-seq shows excellent performance in differential gene expression analysis.

Figure 3

Download asset Open asset

Previous studies found that Tn5 exhibits a slight insertion bias on dsDNA substrates (Goryshin et al., 1998; Green et al., 2012; Lodge et al., 1988). To further investigate whether potential bias exists for TRACE-seq, we thus characterized sites of Tn5-catalyzed adaptor insertion by calculating nucleotide composition of the first 30 bp of each sequence read per library. Similar to dsDNA substrates, we also observed an apparent insertion signature on RNA/DNA hybrids (Figure 2—figure supplement 1i). Nevertheless, per-position information contents were extremely low, suggesting such insertion bias is less likely to affect gene body coverage (Figure 2—figure supplement 1j). Overall, in spite of gene body coverage bias, TRACE-seq allows construction of high complexity RNA libraries and demonstrates similar performance as traditional RNA library preparation methods in terms of detected genes (97% and 93% overlapped with NEBNext and Smart-seq2 library respectively), gene expression measurement (R > 0.90) and differential expression analysis (R > 0.99), but outcompetes the traditional methods in terms of speed, convenience and cost.

Based on substrate diversity and the conserved catalytic domain of the retroviral integrase superfamily including the Tn5 transposase, we envision in this study that Tn5 may be able to directly tagment RNA/DNA hybrid duplexes, in addition to its canonical dsDNA substrates. Having validated such in vitro tagmentation activity, we developed TRACE-seq, which enables one-tube, low-input and low-cost library construction for RNA-seq experiments and demonstrates excellent performance in DE analysis. Compared to conventional RNA-seq methods, TRACE-seq does not need to pre-extract mRNA and synthesize a second DNA chain after mRNA reverse transcription. Therefore, TRACE-seq bypasses laborious and time-consuming processes, is compatible with low input, and reduces reagent cost (Supplementary file 4). During the preparation of this paper, an independent study also reported the similar finding and developed a RNA-seq method named SHERRY (Di et al., 2020). The major conclusions are very consistent between the two studies.

Despite its unique advantages, there is room to further improve TRACE-seq. For instance, the libraries generated by TRACE-seq in its current form are not strand-specific, which is a significant drawback for RNA-Seq experiments. Yet, TRACE-seq should be able to be converted to 5’ RNA-seq or 3’ RNA-seq, which can directionally preserve the 5’ or 3’ end information of transcripts (Cole et al., 2018; Pallares et al., 2020). In addition, TRACE-seq could be used in multiplex profiling when utilizing Tn5 transposase containing barcoded adaptors (Cusanovich et al., 2015; Zhu et al., 2019). Besides, if home-made Tn5 can be used (as have done in Picelli et al., 2014a; Kia et al., 2017; Kaya-Okur et al., 2019), the costs will be further cut down. The in vitro tagmentation efficiency of Tn5 on RNA/DNA hybrids could also be further improved. We have shown that the addition of PEG200 substantially enhanced the tagmentation efficiency of hybrids. It is also tempting to speculate that hyperactive mutants towards RNA/DNA hybrids could also be obtained through screening and protein engineering, as wild-type Tn5 transposase has been engineered to obtain hyperactive forms (Goryshin and Reznikoff, 1998; Wiegand and Reznikoff, 1992; Weinreich et al., 1994; Zhou and Reznikoff, 1997). Such hyperactive mutants are expected to have immediate utility in single-cell RNA-seq experiments, for instance. Moreover, Tn5 transposition in vivo has been harnessed to profile chromatin accessibility in ATAC-seq (Buenrostro et al., 2013); it remains to be seen whether or not an equivalent version may exist to enable in vivo detection of R-loop, chromatin bound long non-coding RNA and epitranscriptome analysis (Santos-Pereira and Aguilera, 2015; Li and Fu, 2019; Li et al., 2016; Song et al., 2020). To summarize, TRACE-seq manifests a ‘cryptic’ activity of the Tn5 transposase as a powerful tool, which may have broad biomedical applications in the future.

High-throughput sequence data has been deposited in GEO under accession code GSE143422.

1. 1. Lu B
   2. Dong L
   3. Yi D
   4. Zhang M
   5. Yi C

   (2020) NCBI Gene Expression Omnibus

   ID GSE143422. Transposase assisted tagmentation of RNA/DNA hybrid duplexes.

   https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE143422

1. 1. Adey A
   2. Morrison HG
   3. Xun X
   4. Kitzman JO
   5. Turner EH
   6. Stackhouse B
   7. MacKenzie AP
   8. Caruccio NC
   9. Zhang X
   10. Shendure J

   (2010) Rapid, low-input, low-bias construction of shotgun fragment libraries by high-density in vitro transposition

   Genome Biology 11:R119.

   https://doi.org/10.1186/gb-2010-11-12-r119

   • PubMed
   • Google Scholar
2. 1. Anders S
   2. Pyl PT
   3. Huber W

   (2015) HTSeq--a Python framework to work with high-throughput sequencing data

   Bioinformatics 31:166–169.

   https://doi.org/10.1093/bioinformatics/btu638

   • PubMed
   • Google Scholar
3. 1. Buenrostro JD
   2. Giresi PG
   3. Zaba LC
   4. Chang HY
   5. Greenleaf WJ

   (2013) Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position

   Nature Methods 10:1213–1218.

   https://doi.org/10.1038/nmeth.2688

   • PubMed
   • Google Scholar
4. 1. Caruccio N

   (2011) Preparation of next-generation sequencing libraries using nextera technology: simultaneous DNA fragmentation and adaptor tagging by in vitro transposition

   Methods in Molecular Biology 733:241–255.

   https://doi.org/10.1007/978-1-61779-089-8_17

   • PubMed
   • Google Scholar
5. 1. Cole C
   2. Byrne A
   3. Beaudin AE
   4. Forsberg EC
   5. Vollmers C

   (2018) Tn5Prime, a Tn5 based 5' capture method for single cell RNA-seq

   Nucleic Acids Research 46:e62.

   https://doi.org/10.1093/nar/gky182

   • PubMed
   • Google Scholar
6. 1. Curcio MJ
   2. Derbyshire KM

   (2003) The outs and ins of transposition: from mu to kangaroo

   Nature Reviews Molecular Cell Biology 4:865–877.

   https://doi.org/10.1038/nrm1241

   • PubMed
   • Google Scholar
7. 1. Cusanovich DA
   2. Daza R
   3. Adey A
   4. Pliner HA
   5. Christiansen L
   6. Gunderson KL
   7. Steemers FJ
   8. Trapnell C
   9. Shendure J

   (2015) Multiplex single cell profiling of chromatin accessibility by combinatorial cellular indexing

   Science 348:910–914.

   https://doi.org/10.1126/science.aab1601

   • PubMed
   • Google Scholar
8. 1. Daley T
   2. Smith AD

   (2013) Predicting the molecular complexity of sequencing libraries

   Nature Methods 10:325–327.

   https://doi.org/10.1038/nmeth.2375

   • PubMed
   • Google Scholar
9. 1. Di L
   2. Fu Y
   3. Sun Y
   4. Li J
   5. Liu L
   6. Yao J
   7. Wang G
   8. Wu Y
   9. Lao K
   10. Lee RW
   11. Zheng G
   12. Xu J
   13. Oh J
   14. Wang D
   15. Xie XS
   16. Huang Y
   17. Wang J

   (2020) RNA sequencing by direct tagmentation of RNA/DNA hybrids

   PNAS 117:2886–2893.

   https://doi.org/10.1073/pnas.1919800117

   • PubMed
   • Google Scholar
10. 1. Dobin A
    2. Davis CA
    3. Schlesinger F
    4. Drenkow J
    5. Zaleski C
    6. Jha S
    7. Batut P
    8. Chaisson M
    9. Gingeras TR

    (2013) STAR: ultrafast universal RNA-seq aligner

    Bioinformatics 29:15–21.

    https://doi.org/10.1093/bioinformatics/bts635

    • PubMed
    • Google Scholar
11. 1. Eisenberg E
    2. Levanon EY

    (2013) Human housekeeping genes, revisited

    Trends in Genetics 29:569–574.

    https://doi.org/10.1016/j.tig.2013.05.010

    • PubMed
    • Google Scholar
12. 1. Finnegan DJ

    (1989) Eukaryotic transposable elements and genome evolution

    Trends in Genetics 5:103–107.

    https://doi.org/10.1016/0168-9525(89)90039-5

    • PubMed
    • Google Scholar
13. 1. Gertz J
    2. Varley KE
    3. Davis NS
    4. Baas BJ
    5. Goryshin IY
    6. Vaidyanathan R
    7. Kuersten S
    8. Myers RM

    (2012) Transposase mediated construction of RNA-seq libraries

    Genome Research 22:134–141.

    https://doi.org/10.1101/gr.127373.111

    • PubMed
    • Google Scholar
14. 1. Goryshin IY
    2. Miller JA
    3. Kil YV
    4. Lanzov VA
    5. Reznikoff WS

    (1998) Tn5/IS50 target recognition

    PNAS 95:10716–10721.

    https://doi.org/10.1073/pnas.95.18.10716

    • PubMed
    • Google Scholar
15. 1. Goryshin IY
    2. Reznikoff WS

    (1998) Tn5 in vitro transposition

    The Journal of Biological Chemistry 273:7367–7374.

    https://doi.org/10.1074/jbc.273.13.7367

    • PubMed
    • Google Scholar
16. 1. Green B
    2. Bouchier C
    3. Fairhead C
    4. Craig NL
    5. Cormack BP

    (2012) Insertion site preference of mu, Tn5, and Tn7 transposons

    Mobile DNA 3:3.

    https://doi.org/10.1186/1759-8753-3-3

    • Google Scholar
17. 1. Hartley SW
    2. Mullikin JC

    (2015) QoRTs: a comprehensive toolset for quality control and data processing of RNA-Seq experiments

    BMC Bioinformatics 16:224.

    https://doi.org/10.1186/s12859-015-0670-5

    • PubMed
    • Google Scholar
18. 1. Kaya-Okur HS
    2. Wu SJ
    3. Codomo CA
    4. Pledger ES
    5. Bryson TD
    6. Henikoff JG
    7. Ahmad K
    8. Henikoff S

    (2019) CUT&tag for efficient epigenomic profiling of small samples and single cells

    Nature Communications 10:1930.

    https://doi.org/10.1038/s41467-019-09982-5

    • PubMed
    • Google Scholar
19. 1. Kia A
    2. Gloeckner C
    3. Osothprarop T
    4. Gormley N
    5. Bomati E
    6. Stephenson M
    7. Goryshin I
    8. He MM

    (2017) Improved genome sequencing using an engineered transposase

    BMC Biotechnology 17:6.

    https://doi.org/10.1186/s12896-016-0326-1

    • PubMed
    • Google Scholar
20. 1. Kleckner N

    (1981) Transposable elements in prokaryotes

    Annual Review of Genetics 15:341–404.

    https://doi.org/10.1146/annurev.ge.15.120181.002013

    • PubMed
    • Google Scholar
21. 1. Langmead B
    2. Salzberg SL

    (2012) Fast gapped-read alignment with bowtie 2

    Nature Methods 9:357–359.

    https://doi.org/10.1038/nmeth.1923

    • PubMed
    • Google Scholar
22. 1. Li H
    2. Handsaker B
    3. Wysoker A
    4. Fennell T
    5. Ruan J
    6. Homer N
    7. Marth G
    8. Abecasis G
    9. Durbin R
    10. 1000 Genome Project Data Processing Subgroup

    (2009) The sequence alignment/Map format and SAMtools

    Bioinformatics 25:2078–2079.

    https://doi.org/10.1093/bioinformatics/btp352

    • PubMed
    • Google Scholar
23. 1. Li X
    2. Xiong X
    3. Yi C

    (2016) Epitranscriptome sequencing technologies: decoding RNA modifications

    Nature Methods 14:23–31.

    https://doi.org/10.1038/nmeth.4110

    • PubMed
    • Google Scholar
24. 1. Li X
    2. Fu XD

    (2019) Chromatin-associated RNAs as facilitators of functional genomic interactions

    Nature Reviews Genetics 20:503–519.

    https://doi.org/10.1038/s41576-019-0135-1

    • PubMed
    • Google Scholar
25. 1. Lodge JK
    2. Weston-Hafer K
    3. Berg DE

    (1988)

    Transposon Tn5 target specificity: preference for insertion at G/C pairs

    Genetics 120:645–650.

    • PubMed
    • Google Scholar
26. 1. Love MI
    2. Huber W
    3. Anders S

    (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2

    Genome Biology 15:550.

    https://doi.org/10.1186/s13059-014-0550-8

    • PubMed
    • Google Scholar
27. 1. Nowotny M

    (2009) Retroviral integrase superfamily: the structural perspective

    EMBO Reports 10:144–151.

    https://doi.org/10.1038/embor.2008.256

    • PubMed
    • Google Scholar
28. 1. Pallares LF
    2. Picard S
    3. Ayroles JF

    (2020) TM3’seq: a tagmentation-mediated 3’sequencing approach for improving scalability of RNAseq experiments

    G3: Genes, Genomes, Genetics 10:143–150.

    https://doi.org/10.1534/g3.119.400821

    • Google Scholar
29. 1. Picelli S
    2. Björklund AK
    3. Reinius B
    4. Sagasser S
    5. Winberg G
    6. Sandberg R

    (2014a) Tn5 transposase and tagmentation procedures for massively scaled sequencing projects

    Genome Research 24:2033–2040.

    https://doi.org/10.1101/gr.177881.114

    • PubMed
    • Google Scholar
30. 1. Picelli S
    2. Faridani OR
    3. Björklund AK
    4. Winberg G
    5. Sagasser S
    6. Sandberg R

    (2014b) Full-length RNA-seq from single cells using Smart-seq2

    Nature Protocols 9:171–181.

    https://doi.org/10.1038/nprot.2014.006

    • PubMed
    • Google Scholar
31. 1. Pramanik S
    2. Nagatoishi S
    3. Saxena S
    4. Bhattacharyya J
    5. Sugimoto N

    (2011) Conformational flexibility influences degree of hydration of nucleic acid hybrids

    The Journal of Physical Chemistry B 115:13862–13872.

    https://doi.org/10.1021/jp207856p

    • PubMed
    • Google Scholar
32. 1. Ramsköld D
    2. Luo S
    3. Wang YC
    4. Li R
    5. Deng Q
    6. Faridani OR
    7. Daniels GA
    8. Khrebtukova I
    9. Loring JF
    10. Laurent LC
    11. Schroth GP
    12. Sandberg R

    (2012) Full-length mRNA-Seq from single-cell levels of RNA and individual circulating tumor cells

    Nature Biotechnology 30:777–782.

    https://doi.org/10.1038/nbt.2282

    • PubMed
    • Google Scholar
33. 1. Rice PA
    2. Baker TA

    (2001) Comparative architecture of transposase and integrase complexes

    Nature Structural Biology 8:302–307.

    https://doi.org/10.1038/86166

    • PubMed
    • Google Scholar
34. 1. Robinson JT
    2. Thorvaldsdóttir H
    3. Winckler W
    4. Guttman M
    5. Lander ES
    6. Getz G
    7. Mesirov JP

    (2011) Integrative genomics viewer

    Nature Biotechnology 29:24–26.

    https://doi.org/10.1038/nbt.1754

    • PubMed
    • Google Scholar
35. 1. Santos-Pereira JM
    2. Aguilera A

    (2015) R loops: new modulators of genome dynamics and function

    Nature Reviews Genetics 16:583–597.

    https://doi.org/10.1038/nrg3961

    • PubMed
    • Google Scholar
36. 1. Savilahti H
    2. Rice PA
    3. Mizuuchi K

    (1995) The phage mu transpososome core: dna requirements for assembly and function

    The EMBO Journal 14:4893–4903.

    https://doi.org/10.1002/j.1460-2075.1995.tb00170.x

    • PubMed
    • Google Scholar
37. 1. Song J
    2. Zhuang Y
    3. Zhu C
    4. Meng H
    5. Lu B
    6. Xie B
    7. Peng J
    8. Li M
    9. Yi C

    (2020) Differential roles of human PUS10 in miRNA processing and tRNA pseudouridylation

    Nature Chemical Biology 16:160–169.

    https://doi.org/10.1038/s41589-019-0420-5

    • PubMed
    • Google Scholar
38. 1. Trapnell C
    2. Williams BA
    3. Pertea G
    4. Mortazavi A
    5. Kwan G
    6. van Baren MJ
    7. Salzberg SL
    8. Wold BJ
    9. Pachter L

    (2010) Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation

    Nature Biotechnology 28:511–515.

    https://doi.org/10.1038/nbt.1621

    • PubMed
    • Google Scholar
39. 1. Wang L
    2. Wang S
    3. Li W

    (2012) RSeQC: quality control of RNA-seq experiments

    Bioinformatics 28:2184–2185.

    https://doi.org/10.1093/bioinformatics/bts356

    • PubMed
    • Google Scholar
40. 1. Weinreich MD
    2. Gasch A
    3. Reznikoff WS

    (1994) Evidence that the Cis preference of the Tn5 transposase is caused by nonproductive multimerization

    Genes & Development 8:2363–2374.

    https://doi.org/10.1101/gad.8.19.2363

    • PubMed
    • Google Scholar
41. 1. Wiegand TW
    2. Reznikoff WS

    (1992) Characterization of two hypertransposing Tn5 mutants

    Journal of Bacteriology 174:1229–1239.

    https://doi.org/10.1128/JB.174.4.1229-1239.1992

    • PubMed
    • Google Scholar
42. 1. Yang W
    2. Steitz TA

    (1995) Recombining the structures of HIV integrase, RuvC and RNase H

    Structure 3:131–134.

    https://doi.org/10.1016/S0969-2126(01)00142-3

    • PubMed
    • Google Scholar
43. 1. Zhou M
    2. Reznikoff WS

    (1997) Tn5 transposase mutants that alter DNA binding specificity

    Journal of Molecular Biology 271:362–373.

    https://doi.org/10.1006/jmbi.1997.1188

    • PubMed
    • Google Scholar
44. 1. Zhu C
    2. Yu M
    3. Huang H
    4. Juric I
    5. Abnousi A
    6. Hu R
    7. Lucero J
    8. Behrens MM
    9. Hu M
    10. Ren B

    (2019) An ultra high-throughput method for single-cell joint analysis of open chromatin and transcriptome

    Nature Structural & Molecular Biology 26:1063–1070.

    https://doi.org/10.1038/s41594-019-0323-x

    • PubMed
    • Google Scholar

Author details

1. Bo Lu

   1. State Key Laboratory of Protein and Plant Gene Research, School of Life Sciences, Peking University, Beijing, China
   2. Peking-Tsinghua Center for Life Sciences, Peking University, Beijing, China

   Contribution

   Conceptualization, Validation, Investigation, Visualization, Methodology, Writing - original draft

   Contributed equally with

   Liting Dong, Danyang Yi and Meiling Zhang

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-5852-0477

2. Liting Dong

   1. State Key Laboratory of Protein and Plant Gene Research, School of Life Sciences, Peking University, Beijing, China
   2. Peking-Tsinghua Center for Life Sciences, Peking University, Beijing, China

   Contribution

   Validation, Investigation, Visualization, Methodology, Writing - original draft

   Contributed equally with

   Bo Lu, Danyang Yi and Meiling Zhang

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0001-8396-374X

3. Danyang Yi

   State Key Laboratory of Protein and Plant Gene Research, School of Life Sciences, Peking University, Beijing, China

   Contribution

   Validation, Investigation, Visualization, Methodology, Writing - original draft

   Contributed equally with

   Bo Lu, Liting Dong and Meiling Zhang

   Competing interests

   No competing interests declared

4. Meiling Zhang

   State Key Laboratory of Protein and Plant Gene Research, School of Life Sciences, Peking University, Beijing, China

   Contribution

   Validation, Investigation, Methodology, Writing - review and editing

   Contributed equally with

   Bo Lu, Liting Dong and Danyang Yi

   Competing interests

   No competing interests declared

5. Chenxu Zhu

   State Key Laboratory of Protein and Plant Gene Research, School of Life Sciences, Peking University, Beijing, China

   Contribution

   Conceptualization

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0003-4216-6562

6. Xiaoyu Li

   State Key Laboratory of Protein and Plant Gene Research, School of Life Sciences, Peking University, Beijing, China

   Contribution

   Investigation

   Competing interests

   No competing interests declared

7. Chengqi Yi

   1. State Key Laboratory of Protein and Plant Gene Research, School of Life Sciences, Peking University, Beijing, China
   2. Peking-Tsinghua Center for Life Sciences, Peking University, Beijing, China
   3. Department of Chemical Biology and Synthetic and Functional Biomolecules Center, College of Chemistry and Molecular Engineering, Peking University, Beijing, China

   Contribution

   Conceptualization, Supervision, Funding acquisition, Project administration, Writing - review and editing

   For correspondence

   chengqi.yi@pku.edu.cn

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0003-2540-9729

• 14,685

  views

• 1,438

  downloads

• 45

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.