CoChat
Sign Up
Cited by 2,239

The DOI matches its Crossref record; 1 of the 1 references it cites were found in a registry or an index.

Circular RNAs Are the Predominant Transcript Isoform from Hundreds of Human Genes in Diverse Cell Types

Salzman, Julia, Gawad, Charles, Wang, Peter Lincoln, Lacayo, Norman, Brown, Patrick O.·2012·PLoS ONE
DOI 10.1371/journal.pone.0030733PMC PMC3270023
Open this paper in CoChat — free account

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1371/journal.pone.0030733

Cites

1 of 1 references verified0 citation errors · 0 not checkable

Statistics

Figures

Figure 1 Models to explain exon scrambling. The canonical linear reference transcript is depicted with exons as colored boxes with four exons 1, 2, 3, and 4. Two simple models of RNA structure that could explain scrambled transcripts are depicted at left and right. At left, model 1 depicts how a scrambled exon 3-exon 2 junction could arise from a tandem duplication of exons 3 and 2, positioning the first copy of exon 3 upstream of exon 2. At the RNA level, this event could arise from post-transcriptional exon rearrangement, or a genomic duplication of exons 2 and 3. Under the model of tandem duplication, when one side of a paired-end read maps to the junction between exon 3 and 2, the other may map to any of exons 1, 2, 3 or 4 with probabilities determined by the library's insert length distribution and the exon lengths. Our data supports paired-end mapping between a junction and exons 2 or 3, but not exons 1 and 4. We note that in principle, the scrambled exon 3 - exon 2 junction could arise from other splicing events and does not necessarily entail tandem duplication. At right, model 2 depicts how a scrambled exon 3 - exon 2 junction could arise from splicing of exons 2 and 3 into a circular RNA molecule, again positioning exon 3 upstream of exon 2. In this model, when one side of a paired-end read maps to the junction between exon 3 and 2, the other will map to exon 2 or exon 3.
Figure 1 Models to explain exon scrambling. The canonical linear reference transcript is depicted with exons as colored boxes with four exons 1, 2, 3, and 4. Two simple models of RNA structure that could explain scrambled transcripts are depicted at left and right. At left, model 1 depicts how a scrambled exon 3-exon 2 junction could arise from a tandem duplication of exons 3 and 2, positioning the first copy of exon 3 upstream of exon 2. At the RNA level, this event could arise from post-transcriptional exon rearrangement, or a genomic duplication of exons 2 and 3. Under the model of tandem duplication, when one side of a paired-end read maps to the junction between exon 3 and 2, the other may map to any of exons 1, 2, 3 or 4 with probabilities determined by the library's insert length distribution and the exon lengths. Our data supports paired-end mapping between a junction and exons 2 or 3, but not exons 1 and 4. We note that in principle, the scrambled exon 3 - exon 2 junction could arise from other splicing events and does not necessarily entail tandem duplication. At right, model 2 depicts how a scrambled exon 3 - exon 2 junction could arise from splicing of exons 2 and 3 into a circular RNA molecule, again positioning exon 3 upstream of exon 2. In this model, when one side of a paired-end read maps to the junction between exon 3 and 2, the other will map to exon 2 or exon 3.CC BY · PubMed Central
Figure 2 Expression levels of scrambled exons. Analysis of paired-end RNA-Seq data from random primed libraries reveals evidence that scrambled exons are present at high stoichiometries compared to the canonical linear transcript transcribed from a large number of human genes. This phenomenon persists across cell types and is illustrated by the expression patterns of 3 leukocyte cell types: CD19 (B cells), CD34 (stem cells) and neutrophils. The fraction of each scrambled transcript as a fraction of total gene expression is computed. The bar plot depicts the number of circular isoforms with estimated abundance relative to all transcripts of the gene in the following ranges: between 0–25%, 25–50%, 50–75% and 75+%. Hundreds of isoforms in each cell type are estimated to represent more than half of all transcripts from each gene.
Figure 2 Expression levels of scrambled exons. Analysis of paired-end RNA-Seq data from random primed libraries reveals evidence that scrambled exons are present at high stoichiometries compared to the canonical linear transcript transcribed from a large number of human genes. This phenomenon persists across cell types and is illustrated by the expression patterns of 3 leukocyte cell types: CD19 (B cells), CD34 (stem cells) and neutrophils. The fraction of each scrambled transcript as a fraction of total gene expression is computed. The bar plot depicts the number of circular isoforms with estimated abundance relative to all transcripts of the gene in the following ranges: between 0–25%, 25–50%, 50–75% and 75+%. Hundreds of isoforms in each cell type are estimated to represent more than half of all transcripts from each gene.CC BY · PubMed Central
Figure 3 RNaseR assay confirms scrambled exons arise from circular RNA. Panel A: Total RNA from HeLa cells was digested with RNaseR at varying enzyme concentrations (0, 3, 10, and 100 units) after the RNA was depleted of ribosomal RNA. Primers capable of amplifying the canonical linear transcript and the predicted circular transcript (by outward facing primers within a single exon predicted in the scramble) were used in a RT-PCR experiment for each of the digestion conditions. Canonical transcripts were consistently degraded by RNaseR, only detectable by PCR at 0 units of RNaseR, whereas predicted circular transcripts consistently resisted the RNaseR challenge, providing strong evidence of circularity. FBXW4 and MAN1A2 respectively show 2 and 4 circular isoforms, both of which were predicted by the sequencing data. The predicted lengths of circular isoforms are respectively a 3-2 junction of CAMSAP1 (predicted to produce a 435 bp circle), a 4-2 and 5-2 junction of FBXW4 (predicted to produce 415 and 510 bp circles), a 4-2, 5-2 and 6-2 junction of MAN1A2 (predicted to produce 471, 553, and 648 bp circles), a 3-3 junction in REXO4 (predicted to produce a 338 bp circle), a 2-2 junction of RNF220 (predicted to produce a 742 bp circle) and a 3-2 junction of ZKSCAN1 (predicted to produce a 667 bp circle). Panel B: A northern blot on total and cytoplasmic lysate from HeLa cells shows hybridization of a 481 bp probe complementary to the MAN1A2 5-2 exon scramble. 3.7 and 6.2 ug of total and cytoplasmic RNA were loaded onto a 1% agarose gel and 10 pM of probe was hybridized for 24–48 hours. Detection was performed using the BrightStar BioDetect Kit (Ambion, Austin, TX). The specific band at 553 bp corresponds to the predicted size of a circular RNA containing exons 2,3,4 and 5 of MAN1A2.
Figure 3 RNaseR assay confirms scrambled exons arise from circular RNA. Panel A: Total RNA from HeLa cells was digested with RNaseR at varying enzyme concentrations (0, 3, 10, and 100 units) after the RNA was depleted of ribosomal RNA. Primers capable of amplifying the canonical linear transcript and the predicted circular transcript (by outward facing primers within a single exon predicted in the scramble) were used in a RT-PCR experiment for each of the digestion conditions. Canonical transcripts were consistently degraded by RNaseR, only detectable by PCR at 0 units of RNaseR, whereas predicted circular transcripts consistently resisted the RNaseR challenge, providing strong evidence of circularity. FBXW4 and MAN1A2 respectively show 2 and 4 circular isoforms, both of which were predicted by the sequencing data. The predicted lengths of circular isoforms are respectively a 3-2 junction of CAMSAP1 (predicted to produce a 435 bp circle), a 4-2 and 5-2 junction of FBXW4 (predicted to produce 415 and 510 bp circles), a 4-2, 5-2 and 6-2 junction of MAN1A2 (predicted to produce 471, 553, and 648 bp circles), a 3-3 junction in REXO4 (predicted to produce a 338 bp circle), a 2-2 junction of RNF220 (predicted to produce a 742 bp circle) and a 3-2 junction of ZKSCAN1 (predicted to produce a 667 bp circle). Panel B: A northern blot on total and cytoplasmic lysate from HeLa cells shows hybridization of a 481 bp probe complementary to the MAN1A2 5-2 exon scramble. 3.7 and 6.2 ug of total and cytoplasmic RNA were loaded onto a 1% agarose gel and 10 pM of probe was hybridized for 24–48 hours. Detection was performed using the BrightStar BioDetect Kit (Ambion, Austin, TX). The specific band at 553 bp corresponds to the predicted size of a circular RNA containing exons 2,3,4 and 5 of MAN1A2.CC BY · PubMed Central
Figure 4 Scrambled exons are enriched in poly-A depleted samples. Single-end 76-bp RNA-Seq was performed on matched experiments on HeLa, and H9 Human embryonic stem cell lysates were polyA selected and polyA depleted (data from Yang et al [22]). The numbers of scrambled exons detected in each sample which appeared in our curated database of scrambled junctions from the leukocyte data are depicted as colored bars. Roughly equal numbers of sequencing reads were available from each of 4 samples. Left panels of bar plot: both H9 and HeLa cells show markedly more exon scrambles in polyA depleted fractions compared to polyA enriched fractions, consistent with scrambles arising from circular transcripts which lack polyA tails. Right panels of bar plot: conversely, in the much smaller subset of scrambled exon pairs where we have evidence of internal tandem duplication (i.e. evidence against circularity), we find the opposite enrichment: more exon scrambles in polyA enriched fractions compared to polyA depleted fractions, consistent with this small subset of scrambles arising from linear, polyA transcripts.
Figure 4 Scrambled exons are enriched in poly-A depleted samples. Single-end 76-bp RNA-Seq was performed on matched experiments on HeLa, and H9 Human embryonic stem cell lysates were polyA selected and polyA depleted (data from Yang et al [22]). The numbers of scrambled exons detected in each sample which appeared in our curated database of scrambled junctions from the leukocyte data are depicted as colored bars. Roughly equal numbers of sequencing reads were available from each of 4 samples. Left panels of bar plot: both H9 and HeLa cells show markedly more exon scrambles in polyA depleted fractions compared to polyA enriched fractions, consistent with scrambles arising from circular transcripts which lack polyA tails. Right panels of bar plot: conversely, in the much smaller subset of scrambled exon pairs where we have evidence of internal tandem duplication (i.e. evidence against circularity), we find the opposite enrichment: more exon scrambles in polyA enriched fractions compared to polyA depleted fractions, consistent with this small subset of scrambles arising from linear, polyA transcripts.CC BY · PubMed Central
Figure 5 qPCR shows scrambled exons are enriched in the cytoplasm. HeLa whole cells lysates were fractioned into cytoplasmic and nuclear. The nuclear localized noncoding RNA XIST served as a control for fractionation:, and as expected, was enriched in the nuclear fraction. In addition, precursor ribosomal RNA bands were present in the nucleus but not the cytoplasm (see Figure S4). Using probes specific to each canonical and circular isoform (corresponding to those examples depicted in Figure 3), we compared Ct values calculated from qPCR on cDNA from the cytoplasmic fraction to the Ct value from qPCR on cDNA from the nuclear fraction. Bar heights show this average Ct value difference across 2 biological replicates. Error bars represent 2.5 standard deviations computed from biological variation of the qPCR assay. These results show that most circular isoforms are more enriched in the cytoplasm compared to the canonical linear isoforms.
Figure 5 qPCR shows scrambled exons are enriched in the cytoplasm. HeLa whole cells lysates were fractioned into cytoplasmic and nuclear. The nuclear localized noncoding RNA XIST served as a control for fractionation:, and as expected, was enriched in the nuclear fraction. In addition, precursor ribosomal RNA bands were present in the nucleus but not the cytoplasm (see Figure S4). Using probes specific to each canonical and circular isoform (corresponding to those examples depicted in Figure 3), we compared Ct values calculated from qPCR on cDNA from the cytoplasmic fraction to the Ct value from qPCR on cDNA from the nuclear fraction. Bar heights show this average Ct value difference across 2 biological replicates. Error bars represent 2.5 standard deviations computed from biological variation of the qPCR assay. These results show that most circular isoforms are more enriched in the cytoplasm compared to the canonical linear isoforms.CC BY · PubMed Central
Figure 6 Models for generation of circular RNA. At left: a schematic diagram of the canonical splicing process splicing out the first intron of the a pre-mRNA of a 4 exon gene, and subsequent removal of introns 2 and 3. Canonical splicing of exon 1 to exon 2 occurs when the splicing machinery catalyzes the formation of the intron lariat and the attack of the free 3′ OH of exon 1 on the 3′ splice site upstream of exon 2. This produces a lariat containing intron 1 and a pre-mRNA with exons 1 and 2 spliced together. At right: a model for the production of circular transcripts. If there is a canonical transcriptional start, and if intron excision does not proceed sequentially in time from the 5′ to 3′ direction of the pre-mRNA, non-canonical pairing of 3′ and 5′ splice sites could be generated. Since the sequences of each 5′ splice site of the pre-mRNA contain the same splicing signals, it is possible that the 3′ splice site upstream of exon 2 is paired with the 5′ splice site downstream of exon 3 and splicing proceeds as if this 5′ splice site were paired with the 3′ splice site upstream of exon 4. In this case, exon 3 would be spliced upstream of exon 2, creating a pre-mRNA intermediate comprised of these two exons and intron 2. Canonical splicing would be predicted to excise this intron, leaving a circular RNA composed of exons 2 and 3. Non-canonical transcription start, as suggested in [25], could produce an orphan 3′ splice site corresponding to the first transcribed exon. This splice site could be paired with a downstream 5′ splice site, generating a circular RNA. In both models, the excised intron would be linear and branched, and expected to be quickly degraded.
Figure 6 Models for generation of circular RNA. At left: a schematic diagram of the canonical splicing process splicing out the first intron of the a pre-mRNA of a 4 exon gene, and subsequent removal of introns 2 and 3. Canonical splicing of exon 1 to exon 2 occurs when the splicing machinery catalyzes the formation of the intron lariat and the attack of the free 3′ OH of exon 1 on the 3′ splice site upstream of exon 2. This produces a lariat containing intron 1 and a pre-mRNA with exons 1 and 2 spliced together. At right: a model for the production of circular transcripts. If there is a canonical transcriptional start, and if intron excision does not proceed sequentially in time from the 5′ to 3′ direction of the pre-mRNA, non-canonical pairing of 3′ and 5′ splice sites could be generated. Since the sequences of each 5′ splice site of the pre-mRNA contain the same splicing signals, it is possible that the 3′ splice site upstream of exon 2 is paired with the 5′ splice site downstream of exon 3 and splicing proceeds as if this 5′ splice site were paired with the 3′ splice site upstream of exon 4. In this case, exon 3 would be spliced upstream of exon 2, creating a pre-mRNA intermediate comprised of these two exons and intron 2. Canonical splicing would be predicted to excise this intron, leaving a circular RNA composed of exons 2 and 3. Non-canonical transcription start, as suggested in [25], could produce an orphan 3′ splice site corresponding to the first transcribed exon. This splice site could be paired with a downstream 5′ splice site, generating a circular RNA. In both models, the excised intron would be linear and branched, and expected to be quickly degraded.CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Show all 1 references

Research works printed but not linked: 27

  1. JM Nigro, 1991, Scrambled exons., Cell, 64, 607, 10.1016/0092-8674(91)90244-S · DOI
  2. C Cocquerelle, 1992, Splicing with inverted order of exons occurs proximal to large introns., The EMBO journal, 11, 1095, 10.1002/j.1460-2075.1992.tb05148.x · DOI
  3. FA Saad, 1992, A 3′ consensus splice mutation in the human dystrophin gene detected by a screening for intra-exonic deletions., Human molecular genetics, 1, 345, 10.1093/hmg/1.5.345 · DOI
  4. B Capel, 1993, Circular transcripts of the testis-determining gene Sry in adult mouse testis., Cell, 73, 1019, 10.1016/0092-8674(93)90279-Y · DOI
  5. QP Cao, 1995, Expression of the mouse testis-determining gene Sry in male preimplantation embryos., Molecular reproduction and development, 40, 196, 10.1002/mrd.1080400208 · DOI
  6. A Hacker, 1995, Expression of Sry, the mouse sex determining gene., Development, 121, 1603, 10.1242/dev.121.6.1603 · DOI
  7. YW Jeske, 1995, Expression of a linear Sry transcript in the mouse genital ridge., Nature genetics, 10, 480, 10.1038/ng0895-480 · DOI
  8. X-FLaJ Lytton, 1999, A Circularized Sodium-Calcium Exchanger Exon 2 Transcript., The Journal of Biological Chemistry, 274, 8153, 10.1074/jbc.274.12.8153 · DOI
  9. PG Zaphiropoulos, 1996, Circular RNAs from transcripts of the rat cytochrome P450 2C24 gene: correlation with exon skipping., Proceedings of the National Academy of Sciences of the United States of America, 93, 6536, 10.1073/pnas.93.13.6536 · DOI
  10. PG Zaphiropoulos, 1997, Exon skipping and circular RNA formation in transcripts of the human cytochrome P-450 2C18 gene in epidermis and of the rat androgen binding protein gene in testis., Molecular and cellular biology, 17, 2985, 10.1128/MCB.17.6.2985 · DOI
  11. JM Houseley, 2006, Noncanonical RNAs from transcripts of the Drosophila muscleblind gene., The Journal of heredity, 97, 253, 10.1093/jhered/esj037 · DOI
  12. CWS Carlos Caldasa, 1998, Exon scrambling of MLL transcripts occur commonly and mimic partial genomic duplication of the gene., Gene, 208, 167, 10.1016/S0378-1119(97)00640-9 · DOI
  13. B Bailleul, 1996, During in vivo maturation of eukaryotic nuclear mRNA, splicing yields excised exon circles., Nucleic acids research, 24, 1015, 10.1093/nar/24.6.1015 · DOI
  14. A Surono, 1999, Circular dystrophin RNAs consisting of exons that were skipped by alternative splicing., Human molecular genetics, 8, 493, 10.1093/hmg/8.3.493 · DOI
  15. CE Burd, 2010, Expression of linear and novel circular forms of an INK4/ARF-associated non-coding RNA correlates with atherosclerosis risk., PLoS genetics, 6, e1001233, 10.1371/journal.pgen.1001233 · DOI
  16. C Cocquerelle, 1993, Mis-splicing yields circular RNA molecules., The FASEB journal : official publication of the Federation of American Societies for Experimental Biology, 7, 155, 10.1096/fasebj.7.1.7678559 · DOI
  17. RJ Dixon, 2005, A genome-wide survey demonstrates widespread non-linear mRNA in expressed sequences from multiple species., Nucleic acids research, 33, 5904, 10.1093/nar/gki893 · DOI
  18. X Shao, 2006, Bioinformatic analysis of exon repetition, exon scrambling and trans-splicing in humans., Bioinformatics, 22, 692, 10.1093/bioinformatics/bti795 · DOI
  19. HH Al-Balool, 2011, Post-transcriptional exon shuffling events in humans can be evolutionarily conserved and abundant., Genome research
  20. J Salzman, 2011, ESRRA-C11orf20 Is a Recurrent Gene Fusion in Serous Ovarian Carcinoma., PLoS biology, 9, e1001156, 10.1371/journal.pbio.1001156 · DOI
  21. L Yang, 2011, Genomewide characterization of non-polyadenylated RNAs., Genome biology, 12, R16, 10.1186/gb-2011-12-2-r16 · DOI
  22. H Suzuki, 2006, Characterization of RNase R-digested cellular RNA source that consists of lariat and circular RNAs from pre-mRNA splicing., Nucleic acids research, 34, e63, 10.1093/nar/gkl151 · DOI
  23. M de la Mata, 2010, First come, first served revisited: factors affecting the same alternative splicing event have different effects on the relative rates of intron removal., RNA, 16, 904, 10.1261/rna.1993510 · DOI
  24. C Schindewolf, 1996, In vitro generation of a circular exon from a linear pre-mRNA transcript., Nucleic acids research, 24, 1260, 10.1093/nar/24.7.1260 · DOI
  25. CY Chen, 1995, Initiation of protein synthesis by the eukaryotic translational apparatus on circular RNAs., Science, 268, 415, 10.1126/science.7536344 · DOI · Paper page
  26. TB Hansen, 2011, miRNA-dependent gene silencing involving Ago2-mediated cleavage of a circular antisense RNA., The EMBO journal, 30, 4414, 10.1038/emboj.2011.359 · DOI
  27. L Poliseno, 2010, A coding-independent function of gene and pseudogene mRNAs regulates tumour biology., Nature, 465, 1033, 10.1038/nature09144 · DOI

Cited by

Sorted by citations

Metadata

Authors
Salzman, J., Gawad, C., Wang, P. L., Lacayo, N., & Brown, P. O.
Venue
PLoS ONE
Publisher
Public Library of Science (PLoS)
Published
2012-02-01
Volume
7
Issue
2
Pages
e30733
ISSN
1932-6203
DOI
10.1371/journal.pone.0030733
PMC
PMC3270023
Cited by
2,239
References
1
Language
en
Licence
creativecommons.org/licenses/by/4.0/ (unspecified)

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review