# Generation of full-length circular RNA libraries for Oxford Nanopore long-read sequencing

Steffen Fuchs, Loélia Babin, Elissa Andraos, Chloé Bessiere, Semjon Willier, Johannes Schulte, Christine Gaspin, Fabienne Meggetto

► To cite this version:

Steffen Fuchs, Loélia Babin, Elissa Andraos, Chloé Bessiere, Semjon Willier, et al.. Generation of full-length circular RNA libraries for Oxford Nanopore long-read sequencing. PLoS ONE, 2022, 17 (9), pp. e0273253. 10.1371/journal.pone.0273253. hal-03773269

HAL Id: hal-03773269

https://hal.science/hal-03773269v1

Submitted on 20 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
LAB PROTOCOL

# Generation of full-length circular RNA libraries for Oxford Nanopore long-read sequencing

Steffen Fuchs<sup>1,2,3,4,5,6*</sup>, Loélia Babin<sup>5,6</sup>, Elissa Andraos<sup>5,6</sup>, Chloé Bessiere<sup>5,6</sup>, Semjon Willier<sup>7</sup>, Johannes H. Schulte<sup>1,2,3,4</sup>, Christine Gaspin<sup>8,9</sup>, Fabienne Meggetto<sup>5,6*</sup>

1. Department of Pediatric Oncology and Hematology, Charité - Universitätsmedizin Berlin, Berlin, Germany,

2. German Cancer Consortium (DKTK), Partner Site Berlin, Berlin, Germany,

3. German Cancer Research Center (DKFZ), Heidelberg, Germany,

4. Berlin Institute of Health at Charité –Universitätsmedizin Berlin, BIH Biomedical Innovation Academy, BIH Charité Clinician Scientist Program, Berlin, Germany,

5. CRCT, Inserm, CNRS, Université Toulouse III-Paul Sabatier, Centre de Recherches en Cancérologie de Toulouse,

Université de Toulouse, Toulouse, France,

6. Laboratoire d'Excellence Toulouse Cancer-TOUCAN, Toulouse, France,

* steffen.fuchs@charite.de, steffen.fuchs@inserm.fr (SF); fabienne.meggetto@inserm.fr (FM)

<figure><img src="image_2.png" /></figure>

## Abstract

Circular RNA (circRNA) is a noncoding RNA class with important implications for gene expression regulation, mostly by interaction with other RNA species or RNA-binding proteins. While the commonly applied short-read Illumina RNA-sequencing techniques can be used to detect circRNAs, their full sequence is not revealed. However, the complete sequence information is needed to analyze potential interactions and thus the mechanism of action of circRNAs. Here, we present an improved protocol to enrich and sequence full-length circRNAs by using the Oxford Nanopore long-read sequencing platform. The protocol involves an enrichment of lowly abundant circRNAs by exonuclease treatment and negative selection of linear RNAs. Then, a cDNA library is created and amplified by PCR. This protocol provides enough material for several sequencing runs. The library is used as input for ligation-based sequencing together with native barcoding. Stringent quality control of the libraries is ensured by a combination of Qubit, Fragment Analyzer and qRT-PCR. Multiplexing of up to 4 libraries yields in total more than 1–2 Million reads per library, of which 1–2% are circRNA-specific reads with >99% of them full-length. The protocol works well with human cancer cell lines. We further provide suggestions for the bioinformatic analysis of the created data, as well as the limitations of our approach together with recommendations for troubleshooting and interpretation. Taken together, this protocol enables reliable full-length analysis of circRNAs, a noncoding RNA type involved in a growing number of physiologic and pathologic conditions.

### Metadata

**Associated content.** https://dx.doi.org/10.17504/protocols.io.rm7vzy8r4lx1/v2.

<figure><img src="image_3.png" /></figure>

### OPEN ACCESS

**Citation:** Fuchs S, Babin L, Andraos E, Bessiere C, Willier S, Schulte JH, et al. (2022) Generation of full-length circular RNA libraries for Oxford Nanopore long-read sequencing. PLoS ONE 17(9): e0273253. https://doi.org/10.1371/journal.pone.0273253

**Editor:** Eduardo Andrés-León, Institute of Parasitology and Biomedicine, SPAIN

**Received:** March 8, 2022

**Accepted:** July 29, 2022

**Published:** September 7, 2022

**Peer Review History:** PLOS recognizes the benefits of transparency in the peer review process; therefore, we enable the publication of all of the content of peer review and author responses alongside final, published articles. The editorial history of this article is available here: https://doi.org/10.1371/journal.pone.0273253

**Copyright:** © 2022 Fuchs et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** Sequencing data generated in this study are available at the NCBI Gene Expression Omnibus (GSE197872; https://www.ncbi.nlm.nih.gov/geo). All other data are
available from the corresponding authors on
reasonable request. Additionally, the DOI for the
protocol mentioned in the lab protocol paper. It is
the following: dx.doi.org/10.17504/protocols.io.
cbs9snh6 https://protocols.io/view/generation-of-
full-length-circrna-libraries-for-ox-cbs9snh6.html.

**Funding:** S.F. is a participant in the BIH-Charité Clinician Scientist Program funded by the Charité – Universitätsmedizin Berlin and the Berlin Institute of Health. S.F. was supported during work by a grant from the Berliner Krebsgesellschaft e.V. (grant no. FUFF201721KK), the Stiftung Tumorforschung Kopf-Hals (Wiesbaden, Germany) and the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation, grant no. 439441203). F.M. was supported by grants from Inserm, l'association Eva pour la vie, the Federation Grandir Sans Cancer and La ligue contre le cancer (Equipes Labelisée 2017-2021). E.A. is supported by a grant from Labex TOUCAN/Laboratoire d'excellence Toulouse Cancer. L.B. and C.B. are supported by fellowships from the Fondation de France. The funders had and will not have a role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared
that no competing interests exist.

**Introduction**

Circular RNAs (circRNAs) are a class of noncoding RNA, which is generated by a form of alternative splicing termed back-splicing. Their ring-like structure and lack of free 5' and 3' ends render them exonuclease resistant and more stable than linear RNAs. This is the reason why they escape detection by the highly used poly(A)-selected mRNA sequencing [1, 2]. circRNAs regulate gene expression by e.g. binding to microRNAs or RNA-binding proteins, or interact directly with the transcriptional machinery [3–5]. The commonly applied Illumina-based short-read sequencing techniques can be readily used (total RNA-seq) to identify circRNAs by their characteristic back-splice junction [6]. However, since circRNAs, which have an average size ranging from 200 to 800 nt [7, 8], share their remaining sequence with their cognate linear RNAs from the same gene, the full-length information cannot be confidently retrieved by these methods. The recently developed long-read RNA-sequencing techniques such as the Oxford-Nanopore sequencing platform bear great potential to obtain this missing information due to their ability to sequence transcripts in full-length.

Sequencing of circRNAs by a direct RNA-sequencing protocol with Oxford Nanopore is an attractive option to analyze circRNAs including their epigenetic modifications without introduction of bias by PCR. However, such an approach requires a linearization of circRNA molecules and has a reduced sensitivity, since the detection of lowly abundant circRNAs is limited and the sequencing coverage is not as high as for Illumina-based methods. Wang et al. followed this approach to analyze plant RNA [9], which they fragmented to be able to sequence it. While several circRNAs were detected, the high amount of input RNA remains a critical limitation in particular when analyzing human and especially patient samples that are often degraded.

Recently, a workflow was published to sequence circRNAs by creating a cDNA library that is amplified by PCR, without the need of fragmentation [8]. The approach of Zhang *et al.* uses a ribodepletion followed by an enrichment of circRNAs by exonuclease treatment and a size selection of transcripts longer than 1 kb. The created library is then used for ligation-based sequencing with Oxford Nanopore. With this approach the team obtained between 0.8 and 4 million reads per library, of which 1–6% were circRNA-specific reads that mapped to the back-splice junction. Here, we present a modified version of this workflow that we adapted to retrieve full-length sequencing information also of shorter circRNAs to cover the whole spectrum of circRNAs. The workflow produces an increased library output to sequence several times, if needed, to potentially detect also lowly abundant circRNAs. In more detail, we changed the ribodepletion method from a commercial kit to the published method of Baldwin *et al.* [10], which worked more efficiently in our hands. This ribodepletion method is based on a pool of DNA oligonucleotides that hybridize with ribosomal RNA and a subsequent digest of DNA:RNA hybrids by RNaseH. For further circRNA enrichment a negative selection of poly (A) transcripts was added by using oligo(dT)-conjugated magnetic beads. The final size selection of the amplified circRNA library was adapted to include shorter circRNAs. Finally, we introduced a quality control by qRT-PCR to detect the enrichment of circRNAs and the depletion of unwanted transcripts, such as ribosomal RNA, mitochondrial RNA and small-nucleolar RNAs. Further, we provide recommendations for the Nanopore library creation and sequencing together with recommendations for the bioinformatics analysis.

**Materials and methods**

The protocol described in this peer-reviewed article is published on protocols.io, https://dx.doi.org/10.17504/protocols.io.rm7vzy8r4lx1/v2 (version 2) and is included for printing as S1 File with this article.
<figure><img src="image_4.png" /><figcaption>**Fig 1.** Generated libraries have the size of the average circRNA length. Shown is the library created from RNA of the anaplastic large-cell lymphoma cell line SU-DHL-1. The library size was analyzed by Fragment Analyzer with the kit hs NGS. The average library size was 654 nt. RFU, relative fluorescence units.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0273253.g001

## circRNA validation

1 µg RNA isolated from COST anaplastic large-cell lymphoma (ALCL) cells was used to create cDNA with the ProtoScript II First Strand cDNA Synthesis kit (New England Biolabs, #E6560S) according to the manufacturer using random hexamer primers. An RT-PCR was performed using the FastStart Taq DNA Polymerase (Roche, #4738381001) according to the manufacturer. The used primer sequences are as follows: circARID1A (F: CTCCAGTAAGG-GAGGGCAAG, R: TGTTGCTGCGAGTATGGGTT), circNFATC3 (F: TGAAACTGAAGGTAG CCGAGG, R: ATGTGGTAAGCAAAGTGGTGTG), circASAP1 (F: GGATACAGCATGCATCAG CTC, R: TCAGCTCCTGTTATCTCTGTGC), circSETD3 (F: TCCTTTGGTGACACAGTTGCT, R: ACTCCTTCTGGCAGCCCTAT) and circZBTB46 (F: CCGGTAGTGGGACGTGATTT, R: ACTCGCTGTCCCAGTCTGTA). PCR products were analyzed by agarose gel electrophoresis (2% agarose, 100 V, 30 min.). The band of the expected size was cut and the DNA cleaned by NucleoSpin Gel and PCR Clean-up kit (Macherey Nagel, #740609.50) following the manufacturer's recommendations. The samples were sent for Sanger sequencing to Eurofins Genomics (Ebersberg, Germany). Obtained sequences were analyzed by NCBI BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi) and aligned to the human genome (GRCh38) and the respective circRNA sequence.

## Expected results

Using the described protocol, we prepared sequencing libraries of 4 different anaplastic large-cell lymphoma (ALCL) cell lines that served as a model to test the workflow (SU-DHL-1, Karpas-299, COST [11], SUP-M2). The obtained libraries had an average length of 606.8 nt and a concentration of 4.8 ng/µl (117–133 ng of library in total), respectively (Fig 1, Table 1). The library length corresponds with the average published size of circRNAs, which is between 200–800 nt [7, 8], thus showing that our workflow maintains the size of circRNAs and does not degrade them.

<table><caption>**Table 1.** Results of the library preparation. Libraries for Nanopore sequencing were prepared of 4 anaplastic large-cell lymphoma cell lines (SU-DHL1, Karpas-299, COST, SUP-M2). The concentration was measured with the Qubit BR dsDNA kit and library size by Fragment Analyzer with the hs NGS kit.</caption><thead><tr><th></th><th>SU-DHL-1</th><th>Karpas-299</th><th>COST</th><th>SUP-M2</th><th>Average</th></tr></thead><tbody><tr><td>Concentration [ng/µl]</td><td>5.12</td><td>4.68</td><td>4.92</td><td>4.5</td><td>4.81</td></tr><tr><td>Size [nt]</td><td>654</td><td>629</td><td>573</td><td>571</td><td>606.8</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0273253.t001
<figure><img src="image_5.png" /><figcaption>**Fig 2.** circRNAs get enriched by the library workflow. 4 different RNA samples from anaplastic large-cell lymphoma cell lines (SU-DHL-1, Karpas-299, COST and SUP-M2) were treated enzymatically to enrich for circRNAs as described in the protocol. The expression of circRNAs and unwanted transcripts (ribosomal RNA, 18S rRNA; mitochondrial RNA, mtRNR1; small-nucleolar RNA, RNU6B; signal recognition particle RNA, RN7SL2 and linear RNAs/mRNAs, linkZKSCAN1, linHIPK3) was analyzed by qRT-PCR and compared with an untreated Mock control.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0273253.g002

As part of our introduced quality control workflow, a qRT-PCR was performed to detect the enrichment of circRNAs and the depletion of unwanted RNA transcripts (Fig 2). We used 3 different circRNAs as indicator for an enrichment, of which we know from previous Illumina-based RNA-sequencing experiments that they are well expressed in our cell line models. The enrichment was on average between 4 and 17-fold in comparison to the non-enriched control. Ribosomal RNA was depleted more than 30,000-fold and also other unwanted transcripts were degraded (mitochondrial RNA, small nucleolar RNA, signal recognition particle RNA, linear RNAs and mRNAs).

The created library pool was used as input for Oxford Nanopore sequencing with the ligation-based sequencing kit (SQK-LSK109) together with the native barcoding kit (EXP-NBD104) according to the manufacturer and sequenced on one flow cell on a MinION MK1C. The sequencing output was on average 1,536,229 reads per library and the reads were of high quality (Table 2, mean Q-score 15).

<table><caption>**Table 2.** Sequencing results obtained with one MinION flow cell. circRNA-enriched libraries from 4 anaplastic large-cell lymphoma cell lines were sequenced by Oxford Nanopore. Calculations are based on the passed reads and the circRNA analysis was performed with CIRI-Long. BSJ, back-splice junction.</caption><thead><tr><th></th><th>SU-DHL1</th><th>Karpas-299</th><th>COST</th><th>SUP-M2</th><th>Average</th></tr></thead><tbody><tr><td>Raw reads</td><td>1,473,419</td><td>1,734,196</td><td>899,725</td><td>2,037,577</td><td>1,536,229</td></tr><tr><td>Mean read length [nt]</td><td>459.6</td><td>368.2</td><td>386.3</td><td>403.3</td><td>404.4</td></tr><tr><td>Maximum read length [nt]</td><td>4,006</td><td>3,889</td><td>3,538</td><td>3,455</td><td>3,722</td></tr><tr><td>BSJ-reads (% of reads)</td><td>1.05</td><td>0.95</td><td>0.95</td><td>1.06</td><td>1.00</td></tr><tr><td>Full-length circRNAs</td><td>15,673</td><td>16,725</td><td>8,750</td><td>21,918</td><td>15,767</td></tr><tr><td>Different circRNAs</td><td>3,143</td><td>3,195</td><td>1,426</td><td>4,017</td><td>2,945</td></tr><tr><td>Mean circRNA length [nt]</td><td>435.1</td><td>354.9</td><td>366.4</td><td>370.4</td><td>381.7</td></tr><tr><td>Maximum circRNA length [nt]</td><td>1,798</td><td>1,634</td><td>1,596</td><td>2,228</td><td>1,814</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0273253.t002
<figure><img src="image_6.png" /><figcaption>**Fig 3.** Distribution of the length of identified circRNAs. Shown are the results of circRNA sequencing from 4 anaplastic large-cell lymphoma (ALCL) cell lines with Oxford Nanopore. A) Distribution of read length and B) length of the identified circRNAs in the 4 different ALCL cell line samples. Gaussian kernel density estimation was used to calculate the distribution of circRNA length.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0273253.g003

Importantly, the pores were not completely saturated, so probably a longer sequencing run with more material would have been possible. Base calling of the raw sequencing data was performed with Guppy from the MinKNOW software (v22.05.8, part of the operation system of the MinION sequencing device) and fastq files were generated. Bioinformatics analysis of the fastq files involved cleaning the reads from adapter sequences with cutadapt (v.3.4, https://doi.org/10.14806/ej.17.1.200). Then the CIRI-Long software (v1.0.3, [8]) was used to detect circRNAs using default settings (detailed recommendations are included in our protocol as S1 File). Following the analysis workflow described for CIRI-Long in our protocol we could identify on average 15,767 circRNA-specific reads, thus 1.0% of the total reads, of which 99% covered the full length of the circRNA (Figs 3 and 4), similar to the study from Zhang *et al.* [8]. On average 2,945 different circRNAs were identified. Of note, it is visible that when more reads are generated, more different full-length circRNA isoforms are detected, which could be another argument for deeper sequencing. The results were comparable among the samples from the 4 different human cancer cell lines, which demonstrates the robustness of the workflow.

We then selected randomly 5 different circRNAs detected by Nanopore for validation (Fig 5 and S1 Fig). COST ALCL cells were used to isolate RNA and transcribe cDNA. Primers were designed to specifically amplify the circRNAs, which were submitted for Sanger sequencing. The characteristic back-splice junction essential for the formation of the RNA circle was detected for all of the circRNAs, which validates our sequencing workflow.

The limitations of this protocol involve the relatively low abundance of circRNAs in general. Since Nanopore sequencing creates less reads in comparison to Illumina-based techniques, especially lowly abundant circRNAs might not be detected by our protocol. Further, the recommended RNA input of 7 µg to have enough material for several rounds of
<figure><img src="image_7.png" /><figcaption>**Fig 4.** One long read captures the entire length of one circRNA several times. A) Shown is exemplified the alignment of sequencing reads obtained by Oxford Nanopore sequencing in the anaplastic large-cell lymphoma cell line SU-DHL-1. The alignment to exon 2–4 of the *ARID1A* gene supports a circRNA with a back-splice junction between exon 4 and 2 that was captured by one long-read (arrows). Track 1: genomic location, track 2: sequencing coverage, track 3: predicted splice junctions, track 4: selected sequencing reads and their alignment, track 5: genes and their exons and introns. The alignment was visualized with IGV Genomics Viewer. B) Scheme showing the long read from panel A aligning several times to circARID1A derived from *ARID1A*. Marked is the back-splice junction (BSJ).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0273253.g004

sequencing and also detect lower abundant circRNAs might be too high for situations, where the material is limited. However, we successfully tested our protocol with as little as 3 µg, without impairing the sequencing output. If the obtained amount of DNA by the circRNA enrichment workflow is not high enough, then the PCR amplification of the cDNA libraries can be adapted by increasing the volume of the PCR reaction and the number of cycles (step 12 of the protocol). Further, if the expected size of circRNAs in the cell model of interest is lower or higher than 200–800 nt, the size selection (step 13 of the protocol) can be easily adapted by changing the ratio of beads to DNA (lower ratio selects for larger fragments, higher ratio selects for shorter fragments). Limitations concerning the used tool CIRI-Long to detect circRNAs are that the tool does not detect circRNAs derived from fusion genes. Further, the alignment parameters are not modifiable and a bam file containing the aligned reads is not conserved. Therefore, we routinely perform a separate alignment using minimap2 (v.2.19, [12]) and a visualization by IGV Genomics Viewer (v2.9.4, [13]). By this means, chimeric alignments containing segments of the same read aligning to distant genes can be visualized and help to identify fusion gene derived circRNAs.

In summary, this protocol facilitates consistent full-length sequencing of circRNAs, which will help to study this noncoding RNA type in a variety of physiologic and pathologic contexts.
<figure><img src="image_8.png" /><figcaption>**Fig 5.** circRNAs detected by Nanopore were validated by Sanger sequencing. 3 different circRNAs (A: circZBTB46, B: circSETD3, C: circASAP1) detected by Nanopore-seq were amplified by RT-PCR using cDNA from COST ALCL cells and analyzed by Sanger sequencing. For each panel, the alignment of long reads against the exons that form the circRNA is shown in the upper part as visualized by IGV Genomics Viewer (tracks as in Fig 3). In the lower part a scheme of the circRNA is shown, together with the localization of the used primers and the BSJ-sequence obtained by Sanger sequencing. The circbase.org ID is mentioned [14]. BSJ, back-splice junction.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0273253.g005
**Supporting information**

**S1 Fig. Validation of circRNAs by Sanger sequencing.** 2 circRNAs detected by Nanopore-seq were validated by Sanger sequencing similar as in Fig 5. Representative alignments and the BSJ-sequence obtained by Sanger sequencing are shown for A) circNFATC3 and B) circAR-ID1A (alignments see Fig 4). The circbase.org ID is mentioned [14]. BSJ, back-splice junction. (TIF)

**S1 File. Step-by-step protocol, also available on protocols.io.**
(PDF)

**Acknowledgments**

The authors are grateful to Emeline Sarot and Nathalie Saint-Laurent (Genomic and Transcriptomic facility, Technology Cluster of the Cancer Research Center of Toulouse, INSERM-UMR1037) for their technical assistance. The authors thank further Falk Hertwig and Filippos Klironomos (Charité, Berlin) for help during the establishment of this protocol.

**Author Contributions**

**Conceptualization:** Steffen Fuchs, Fabienne Meggetto.

**Data curation:** Steffen Fuchs, Loélia Babin, Christine Gaspin.

**Formal analysis:** Steffen Fuchs, Loélia Babin, Christine Gaspin.

**Funding acquisition:** Steffen Fuchs, Fabienne Meggetto.

**Investigation:** Steffen Fuchs, Loélia Babin, Elissa Andraos, Chloé Bessiere, Semjon Willier, Johannes H. Schulte, Christine Gaspin.

**Methodology:** Steffen Fuchs.

**Project administration:** Steffen Fuchs, Fabienne Meggetto.

**Resources:** Fabienne Meggetto.

**Software:** Loélia Babin, Chloé Bessiere, Christine Gaspin.

**Supervision:** Steffen Fuchs, Fabienne Meggetto.

**Validation:** Steffen Fuchs, Loélia Babin, Christine Gaspin.

**Visualization:** Steffen Fuchs, Loélia Babin, Christine Gaspin.

**Writing – original draft:** Steffen Fuchs, Fabienne Meggetto.

**Writing – review & editing:** Steffen Fuchs, Fabienne Meggetto.

**References**

1. Franz A, Rabien A, Stephan C, Ralla B, Fuchs S, Jung K, et al. Circular RNAs: a new class of biomarkers as a rising interest in laboratory medicine. Clin Chem Lab Med. 2018; 56(12):1992–2003. https://doi.org/10.1515/cclm-2018-0231 PMID: 29804099.

2. Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013; 495(7441):333–8. Epub 2013/03/01. https://doi.org/10.1038/nature11928 PMID: 23446348.

3. Zheng Q, Bao C, Guo W, Li S, Chen J, Chen B, et al. Circular RNA profiling reveals an abundant circHIPK3 that regulates cell growth by sponging multiple miRNAs. Nature communications. 2016; 7:11215. https://doi.org/10.1038/ncomms11215 PMID: 27050392.
4. Abdelmohsen K, Panda AC, Munk R, Grammatikakis I, Dudekula DB, De S, et al. Identification of HuR target circular RNAs uncovers suppression of PABPN1 translation by CircPABPN1. RNA biology. 2017; 14(3):361-9. Epub 2017/01/13. https://doi.org/10.1080/15476286.2017.1279788 PMID: 28080204.

5. Li Z, Huang C, Bao C, Chen L, Lin M, Wang X, et al. Exon-intron circular RNAs regulate transcription in the nucleus. Nature structural & molecular biology. 2015; 22(3):256-64. https://doi.org/10.1038/nsmb. 2959 PMID: 25664725.

6. Gao Y, Zhang J, Zhao F. Circular RNA identification based on multiple seed matching. Brief Bioinform. 2017. https://doi.org/10.1093/bib/bbx014 PMID: 28334140.

7. Guo JU, Agarwal V, Guo H, Bartel DP. Expanded identification and characterization of mammalian circular RNAs. Genome Biol. 2014; 15(7):409. https://doi.org/10.1186/s13059-014-0409-z PMID: 25070500.

8. Zhang J, Hou L, Zuo Z, Ji P, Zhang X, Xue Y, et al. Comprehensive profiling of circular RNAs with nanopore sequencing and CIRI-long. Nature biotechnology. 2021. Epub 2021/03/13. https://doi.org/10.1038/s41587-021-00842-6 PMID: 33707777.

9. Wang Y, Wang H, Xi F, Wang H, Han X, Wei W, et al. Profiling of circular RNA N(6) -methyladenosine in moso bamboo (Phyllostachys edulis) using nanopore-based direct RNA sequencing. J Integr Plant Biol. 2020; 62(12):1823-38. Epub 2020/08/01. https://doi.org/10.1111/jipb.13002 PMID: 32735361.

10. Baldwin A, Morris AR, Mukherjee N. An Easy, Cost-Effective, and Scalable Method to Deplete Human Ribosomal RNA for RNA-seq. Curr Protoc. 2021; 1(6):e176. Epub 2021/06/25. https://doi.org/10.1002/ cpz1.176 PMID: 34165268.

11. Lamant L, Espinos E, Duplantier M, Dastugue N, Robert A, Allouche M, et al. Establishment of a novel anaplastic large-cell lymphoma-cell line (COST) from a 'small-cell variant' of ALCL. Leukemia. 2004; 18 (10):1693-8. Epub 2004/09/10. https://doi.org/10.1038/sj.leu.2403464 PMID: 15356659.

12. Li H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018; 34(18):3094-100. Epub 2018/05/12. https://doi.org/10.1093/bioinformatics/bty191 PMID: 29750242.

13. Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G, et al. Integrative genomics viewer. Nature biotechnology. 2011; 29(1):24-6. Epub 2011/01/12. https://doi.org/10.1038/nbt. 1754 PMID: 21221095.

14. Glazar P, Papavasileiou P, Rajewsky N. circBase: a database for circular RNAs. Rna. 2014; 20 (11):1666-70. https://doi.org/10.1261/rna.043687.113 PMID: 25234927.