CoChat
Sign Up
Cited by 345MicroRNA in disease regulationPlant Molecular Biology ResearchRNA Research and Splicing

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

MicroRNAs and their isomiRs function cooperatively to target common biological pathways

Open access
Cloonan, Nicole, Wani, Shivangi, Xu, Qinying, Gu, Jian, Lea, Kristi, Heater, Sheila, Barbacioru, Catalin, Steptoe, Anita L, Martin, Hilary C, Nourbakhsh, Ehsan, Krishnan, Keerthana, Gardiner, Brooke, Wang, Xiaohui, Nones, Katia, Steen, Jason A, Matigian, Nicholas A, Wood, David L, Kassahn, Karin S, Waddell, Nic, Shepherd, Jill, Lee, Clarence, Ichikawa, Jeff, McKernan, Kevin, Bramlett, Kelli, Kuersten, Scott, Grimmond, Sean M·2011·Genome Biology
DOI 10.1186/gb-2011-12-12-r126PMC PMC3334621
Open this paper in CoChat — free accountOpen PDF

Abstract

Background: Variants of microRNAs (miRNAs), called isomiRs, are commonly reported in deep-sequencing studies; however, the functional significance of these variants remains controversial. Observational studies show that isomiR patterns are non-random, hinting that these molecules could be regulated and therefore functional, although no conclusive biological role has been demonstrated for these molecules. Results: To assess the biological relevance of isomiRs, we have performed ultra-deep miRNA-seq on ten adult human tissues, and created an analysis pipeline called miRNA-MATE to align, annotate, and analyze miRNAs and their isomiRs. We find that isomiRs share sequence and expression characteristics with canonical miRNAs, and are generally strongly correlated with canonical miRNA expression. A large proportion of isomiRs potentially derive from AGO2 cleavage independent of Dicer. We isolated polyribosome-associated mRNA, captured the mRNA-bound miRNAs, and found that isomiRs and canonical miRNAs are equally associated with translational machinery. Finally, we transfected cells with biotinylated RNA duplexes encoding isomiRs or their canonical counterparts and directly assayed their mRNA targets. These studies allow us to experimentally determine genome-wide mRNA targets, and these experiments showed substantial overlap in functional mRNA networks suppressed by both canonical miRNAs and their isomiRs. Conclusions: Together, these results find isomiRs to be biologically relevant and functionally cooperative partners of canonical miRNAs that act coordinately to target pathways of functionally related genes. This work exposes the complexity of the miRNA-transcriptome, and helps explain a major miRNA paradox: how specific regulation of biological processes can occur when the specificity of miRNA targeting is mediated by only 6 to 11 nucleotides.

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.1186/gb-2011-12-12-r126

Cites

78 of 78 references verified0 citation errors · 0 not checkable

Statistics

Figures

Figure 1 Correlation between miRNA-seq, qRT-PCR, and microarrays for placental miRNAs. Different alignment strategies will produce different correlations due to the characteristics of the validation technology. (a) The Agilent miRNA microarray typically will capture signal from isomiRs due to cross-hybridization, and therefore has a higher correlation to the recursive matching strategy, which identifies both miRNAs and isomiRs. (b) TaqMan LDA qRT-PCR is more sensitive to the exact sequence of the miRNA, and therefore has higher correlations to the adaptor trimming strategy, which was filtered for exact matches to canonical miRNA sequences. In all plots, points represent mean detection values, and bi-directional error bars represent standard error. ρ values indicated are Spearman correlations. These comparisons demonstrate that the miRNA-seq technology and our analysis approach are both sensitive and specific.
Figure 1 Correlation between miRNA-seq, qRT-PCR, and microarrays for placental miRNAs. Different alignment strategies will produce different correlations due to the characteristics of the validation technology. (a) The Agilent miRNA microarray typically will capture signal from isomiRs due to cross-hybridization, and therefore has a higher correlation to the recursive matching strategy, which identifies both miRNAs and isomiRs. (b) TaqMan LDA qRT-PCR is more sensitive to the exact sequence of the miRNA, and therefore has higher correlations to the adaptor trimming strategy, which was filtered for exact matches to canonical miRNA sequences. In all plots, points represent mean detection values, and bi-directional error bars represent standard error. ρ values indicated are Spearman correlations. These comparisons demonstrate that the miRNA-seq technology and our analysis approach are both sensitive and specific.CC BY · PubMed Central
Figure 2 Arm switching frequency in the miRNA-seq data. (a) Arm-specific expression of at least 21 miRNAs varied between tissues. For all 10 tissues (points), log ratios of 5p to 3p arm expression for 21 miRNAs (blue lines) are plotted according to the arm on which dominant expression occurred, clearly illustrating differences between the tissues for each miRNA. (b) Sequence logo plots for hsa-miR-126 showing the switch between expression dominance from the 5p arm (placenta) to the 3p arm (heart, liver, and kidney). The pre-miRNA hairpin is displayed along the bottom.
Figure 2 Arm switching frequency in the miRNA-seq data. (a) Arm-specific expression of at least 21 miRNAs varied between tissues. For all 10 tissues (points), log ratios of 5p to 3p arm expression for 21 miRNAs (blue lines) are plotted according to the arm on which dominant expression occurred, clearly illustrating differences between the tissues for each miRNA. (b) Sequence logo plots for hsa-miR-126 showing the switch between expression dominance from the 5p arm (placenta) to the 3p arm (heart, liver, and kidney). The pre-miRNA hairpin is displayed along the bottom.CC BY · PubMed Central
Figure 3 miRNA and isomiR characteristics. (a) Example of a canonical miRNA and each mutually exclusive isomiR category. IsomiRs that would fit more than one of these categories are classified as 'mixed' isomiRs. (b) Box-whisker plots showing expression of isomiR types within each tissue. (c) Ratio of each isomiR to its canonical miRNA grouped according to isomiR type within each tissue.
Figure 3 miRNA and isomiR characteristics. (a) Example of a canonical miRNA and each mutually exclusive isomiR category. IsomiRs that would fit more than one of these categories are classified as 'mixed' isomiRs. (b) Box-whisker plots showing expression of isomiR types within each tissue. (c) Ratio of each isomiR to its canonical miRNA grouped according to isomiR type within each tissue.CC BY · PubMed Central
Figure 4 Identification of ac-pre-miRNA candidates. (a) Sequence logo profile of hsa-miR-451 showing 5' arm dominance, almost no variation in the 3' end of the mature miRNA, and a trail of 3' isomiR extensions. (b) Schematic diagram illustrating the generation of end-site (5p-ac-pre-miRNA) and start-site (3p-ac-pre-miRNA) isomiRs via AGO2-mediated cleavage of a pre-miRNA. Dashed lines indicate variable regions. (c) Sequence logo profile for hsa-miR-31, which was identified by our analysis as a candidate for biogenesis through the ac-pre-miRNA pathway. Three representative tissues are shown. (d) Wild-type (WT) zebrafish whole embryo versus mutant AGO2 inactive (MZAgo2) whole embryo miRNA and isomiR expression. Each point represents a separate isomiR or miRNA. The thick blue line is the perfect correlation line, the thin blue lines are ± 2-fold visual guides. Red dots indicate isomiRs/miRNAs derived from candidate ac-pre-miRNA hairpins. A bias towards down-regulation in AGO2 mutants is evident in this plot. (e) Box and whiskers plots showing the expression of all isomiRs/miRNAs deriving from ac-pre-miRNA candidates in wild-type (WT) and mutant (MZago2) miRNA-seq samples; 75% of ac-pre-miRNA candidate molecules are not expressed in the mutant sample.
Figure 4 Identification of ac-pre-miRNA candidates. (a) Sequence logo profile of hsa-miR-451 showing 5' arm dominance, almost no variation in the 3' end of the mature miRNA, and a trail of 3' isomiR extensions. (b) Schematic diagram illustrating the generation of end-site (5p-ac-pre-miRNA) and start-site (3p-ac-pre-miRNA) isomiRs via AGO2-mediated cleavage of a pre-miRNA. Dashed lines indicate variable regions. (c) Sequence logo profile for hsa-miR-31, which was identified by our analysis as a candidate for biogenesis through the ac-pre-miRNA pathway. Three representative tissues are shown. (d) Wild-type (WT) zebrafish whole embryo versus mutant AGO2 inactive (MZAgo2) whole embryo miRNA and isomiR expression. Each point represents a separate isomiR or miRNA. The thick blue line is the perfect correlation line, the thin blue lines are ± 2-fold visual guides. Red dots indicate isomiRs/miRNAs derived from candidate ac-pre-miRNA hairpins. A bias towards down-regulation in AGO2 mutants is evident in this plot. (e) Box and whiskers plots showing the expression of all isomiRs/miRNAs deriving from ac-pre-miRNA candidates in wild-type (WT) and mutant (MZago2) miRNA-seq samples; 75% of ac-pre-miRNA candidate molecules are not expressed in the mutant sample.CC BY · PubMed Central
Figure 5 IsomiRs are associated with translational machinery. (a) Schematic diagram showing the laboratory workflow to isolate total miRNA, miRNA associated with mRNA, and miRNAs associated with polysomal mRNA. (b, c) Correlation scatter plots describing miRNA expression in the three preparations described in (a). The noise threshold was set at 25 tpm for these libraries. (b) Correlation scatter plot between miRNA associated with mRNA (x-axis), and miRNA associated with polysomal mRNA (y-axis). The dashed diagonal line represents perfect correlation, while the horizontal and vertical lines represent the likely noise thresholds for these data. (c) Correlation scatter plot between total cellular miRNAs, and polysomal associated miRNA. The dashed diagonal line represents perfect correlation, while the horizontal and vertical lines represent the likely noise thresholds for these data. (d) A sequence logo profile for hsa-miR-31 showing similar profiles for all three HeLa miRNA-seq samples. This example shows an ac-pre-miRNA candidate hairpin with no change in the proportion of intermediates in the translational machinery. (e) Proportional bar graph illustrating the similar distribution of isomiR categories in each of the three HeLa miRNA-seq samples. CHX, cycloheximide, nt, nucleotide.
Figure 5 IsomiRs are associated with translational machinery. (a) Schematic diagram showing the laboratory workflow to isolate total miRNA, miRNA associated with mRNA, and miRNAs associated with polysomal mRNA. (b, c) Correlation scatter plots describing miRNA expression in the three preparations described in (a). The noise threshold was set at 25 tpm for these libraries. (b) Correlation scatter plot between miRNA associated with mRNA (x-axis), and miRNA associated with polysomal mRNA (y-axis). The dashed diagonal line represents perfect correlation, while the horizontal and vertical lines represent the likely noise thresholds for these data. (c) Correlation scatter plot between total cellular miRNAs, and polysomal associated miRNA. The dashed diagonal line represents perfect correlation, while the horizontal and vertical lines represent the likely noise thresholds for these data. (d) A sequence logo profile for hsa-miR-31 showing similar profiles for all three HeLa miRNA-seq samples. This example shows an ac-pre-miRNA candidate hairpin with no change in the proportion of intermediates in the translational machinery. (e) Proportional bar graph illustrating the similar distribution of isomiR categories in each of the three HeLa miRNA-seq samples. CHX, cycloheximide, nt, nucleotide.CC BY · PubMed Central
Figure 6 IsomiRs target similar molecules to their canonical miRNAs. (a) Schematic diagram depicting the laboratory workflow for the pull-down of mRNA targets of a synthetic biotinylated miRNA. (b) The miRNA expression profiles of the two miRNAs and corresponding shifted isomiRs selected for pull-down experiments. (c) A four-way Venn diagram depicting the overlap between the target genes detected by microarray analysis of pulled-down fractions. (d) The distribution of the number of miRNAs/isomiRs that target genes from significantly enriched pathways. There are more genes targeted by three or more isomiRs than would be expected by chance, and less genes targeted by two or fewer isomiRs than would be expected by chance. (e) A four-way Venn diagram displaying the specific targeting of the genes used in (d).
Figure 6 IsomiRs target similar molecules to their canonical miRNAs. (a) Schematic diagram depicting the laboratory workflow for the pull-down of mRNA targets of a synthetic biotinylated miRNA. (b) The miRNA expression profiles of the two miRNAs and corresponding shifted isomiRs selected for pull-down experiments. (c) A four-way Venn diagram depicting the overlap between the target genes detected by microarray analysis of pulled-down fractions. (d) The distribution of the number of miRNAs/isomiRs that target genes from significantly enriched pathways. There are more genes targeted by three or more isomiRs than would be expected by chance, and less genes targeted by two or fewer isomiRs than would be expected by chance. (e) A four-way Venn diagram displaying the specific targeting of the genes used in (d).CC BY · PubMed Central
Figure 7 IsomiRs and miRNAs target similar biological pathways. The 'molecular mechanisms of cancer' pathway (Ingenuity Pathway Analysis) overlaid with rectangles corresponding to the miRNA or isomiR targeting that gene. hsa-miR-10a-5p|{hsa-miR-10a}|21_43| is colored red, hsa-miR-10b-5p|{hsa-miR-10b}|26_48|) is colored yellow, hsa-miR-10a-5p|{isomiR}|22_44| is colored blue, and hsa-miR-10b-5p|{isomiR}|27_49| is colored green.
Figure 7 IsomiRs and miRNAs target similar biological pathways. The 'molecular mechanisms of cancer' pathway (Ingenuity Pathway Analysis) overlaid with rectangles corresponding to the miRNA or isomiR targeting that gene. hsa-miR-10a-5p|{hsa-miR-10a}|21_43| is colored red, hsa-miR-10b-5p|{hsa-miR-10b}|26_48|) is colored yellow, hsa-miR-10a-5p|{isomiR}|22_44| is colored blue, and hsa-miR-10b-5p|{isomiR}|27_49| is colored green.CC BY · PubMed Central
Figure 8 A diagram illustrating how isomiRs may act cooperatively with the canonical miRNAs to increase the signal to noise ratio of mRNA targeting. (a, b) Where only a single miRNA molecule is used for targeting (a), there is no differentiation between potential 'off-target' effects, and targeting the 'desired' core molecular network when the dosage of that miRNA is increased. However, the dosage is increased using isomiRs (b), with slightly different targeting of both core network mRNAs and 'off-target' miRNAs, amplification of specific targeting of the 'core' network can be achieved relative to the 'off-target' effects. The more isomiRs present, the more differentiation could be achieved.
Figure 8 A diagram illustrating how isomiRs may act cooperatively with the canonical miRNAs to increase the signal to noise ratio of mRNA targeting. (a, b) Where only a single miRNA molecule is used for targeting (a), there is no differentiation between potential 'off-target' effects, and targeting the 'desired' core molecular network when the dosage of that miRNA is increased. However, the dosage is increased using isomiRs (b), with slightly different targeting of both core network mRNAs and 'off-target' miRNAs, amplification of specific targeting of the 'core' network can be achieved relative to the 'off-target' effects. The more isomiRs present, the more differentiation could be achieved.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 79 references

Research works printed but not linked: 1

  1. Maroney PA, Yu Y, Fisher J, Nilsen TW: Evidence that microRNAs are associated with translating messenger RNAs in human cells. Nat Struct Mol Biol. 2006, 13: 1102-1107. 10.1038/nsmb1174. · DOI

Cited by

Sorted by citations

Items 201–300 of 345

Related

See all

Metadata

Authors
Cloonan, N., Wani, S., Xu, Q., Gu, J., Lea, K., Heater, S., Barbacioru, C., Steptoe, A. L., Martin, H. C., Nourbakhsh, E., Krishnan, K., Gardiner, B., Wang, X., Nones, K., Steen, J. A., Matigian, N. A., Wood, D. L., Kassahn, K. S., Waddell, N., Shepherd, J., Lee, C., Ichikawa, J., McKernan, K., Bramlett, K., Kuersten, S., & Grimmond, S. M.
Affiliations
The University of Queensland (AU)
Corresponding author
  • Cloonan, Nicole
  • Kuersten, Scott
  • Grimmond, Sean M
ORCID
Venue
Genome Biology
Publisher
Springer Science and Business Media LLC
Published
2011-12-30
Volume
12
Issue
12
Pages
R126-R126
ISSN
  • 1465-6906 (ISSN-L)
  • 1474-760X
DOI
10.1186/gb-2011-12-12-r126
PMC
PMC3334621
OpenAlex
W2119942258
Cited by
345
References
79
Citations by year
2026: 6 · 2025: 13 · 2024: 12 · 2023: 13 · 2022: 24 · 2021: 21 · 2020: 29 · 2019: 27 · 2018: 21 · 2017: 35 · 2016: 23 · 2015: 29 · 2014: 34 · 2013: 39 · 2012: 19
FWCI
10.01
Type
journal article
Language
en
Open access
Licence
Locations
  • Genome biology (publishedVersion, cc-by)
  • Europe PMC (PubMed Central) (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • Kent Academic Repository (University of Kent) (acceptedVersion)
  • Minerva Access (University of Melbourne) (submittedVersion, cc-by)
  • Griffith Research Online (submittedVersion, cc-by)
Funders
  • Medical Research Council
  • Australian Research Council (DP0988754, DP1093164)
  • National Health and Medical Research Council
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Cancer Research › MicroRNA in disease regulation
  • Life Sciences › Agricultural and Biological Sciences › Plant Science › Plant Molecular Biology Research
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › RNA Research and Splicing
Keywords
microRNA isomiRs, miRNA target networks
MeSH terms
Base Sequence, HeLa Cells, Humans, Molecular Sequence Data, RNA, Messenger, Transfection, Sequence Alignment, Biotinylation, Oligonucleotide Array Sequence Analysis, Gene Expression Profiling, MicroRNAs, Ribonuclease III, Gene Regulatory Networks, DEAD-box RNA Helicases, HEK293 Cells, High-Throughput Nucleotide Sequencing, Transcriptome, Argonaute Proteins

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