Cited by 119MicroRNA in disease regulationCancer-related molecular mechanisms researchRNA modifications and cancer

The DOI matches its Crossref record; 140 of the 140 references it cites were found in a registry or an index; 53 of the 69 sources its claims cite were checked for support.

A mouse tissue atlas of small noncoding RNA

Open access
Isakova, Alina, Fehlmann, Tobias, Keller, Andreas, Quake, Stephen R.·2020·Proceedings of the National Academy of Sciences

Abstract

Significance: We report a systematic unbiased analysis of small RNA molecule expression in 11 different tissues of the model organism mouse. We discovered uncharacterized noncoding RNA molecules and identified that ∼30% of total noncoding small RNA transcriptome are distributed across the body in a tissue-specific manner with some also being sexually dimorphic. Distinct distribution patterns of small RNA across the body suggest the existence of tissue-specific mechanisms involved in noncoding RNA processing.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

This study presents a small noncoding RNA atlas across 11 adult mouse tissues, identifying tissue-specific and sex-dimorphic expression across miRNA, snoRNA, snRNA, scaRNA, tRF, and other ncRNA classes. It also reports novel miRNAs, tissue-dependent miRNA arm switching, and integration with scATAC-seq to infer cell-type origins.

Key findings

  • From 140 small RNA-seq libraries across 11 tissues, the authors identified 1,317 miRNAs, 733 snRNAs, 583 snoRNAs, 25 scaRNAs, 346 tRNAs, 22 mitochondrial tRNAs, and 193 other small ncRNAs; ∼30% of the total noncoding small RNA transcriptome was tissue-specific, with 897 of 3,219 detected genes (28%) contributing to tissue-specific signatures at FDR < 1%, and brain containing the highest number of unique transcripts (∼400).
  • Tissue specificity extended beyond miRNAs: >200 snoRNAs showed tissue-specific patterns, Snora35 and Snord116 were brain-exclusive (TSI = 1), Rny1 was high in brain, Terc was highest in bone marrow and spleen, and ∼400 miRNAs were differentially expressed across tissues, nearly one-quarter of them brain-specific; ∼5% of miRNAs switched arm preference between tissues, including mir-141 and mir-350.
  • Sex-biased miRNAs were identified at FDR < 0.01, with kidney and lung containing the most (27 and 18), followed by brain (7), muscle (5), and heart (2); four of five male-specific brain miRNAs came from a 5-kb region of the imprinted Dlk1-Dio3 locus, and all six sex-specific novel miRNAs were male-dominant.
  • Novel miRNA discovery yielded 473 candidate regions, 214 supported by AGO-CLIP, 87 previously reported, and 12 validated; tRFs made up about a quarter of libraries, with nuclear tRNA fragments mostly 33 nt and mitochondrial fragments spanning 18–54 nt; an SVM tissue classifier reached 0.96 accuracy on postnatal ENCODE data and 0.69 on embryonic data.

Methods

  • Small RNA-seq was performed on 11 tissues from adult C57BL/6J mice (4 males, 10 females, ∼3 mo old), generating 140 libraries; RNA was extracted, prepared with Illumina TruSeq Small RNA, and sequenced on NextSeq500.
  • Reads were mapped to GRCm38 with STAR and quantified using GENCODE M20, miRBase v22, and GtRNAdb annotations; tRFs were analyzed with unitas.
  • Analyses included DESeq2 normalization, t-SNE, tissue specificity index (TSI), likelihood-ratio test differential expression, sex comparisons with Wald tests and permutation, target correlation, miRDeep2 for novel miRNAs with AGO-CLIP/literature/RT-qPCR validation, SVM tissue classification trained on 1,973 miRNAs and tested on ENCODE data, and integration with scATAC-seq/Cicero and 10X mouse brain data.

Limitations

  • The paper states that resolution limitations of the mouse ATAC atlas prevented unambiguous cell-type assignment for the majority of brain-specific ncRNAs, and that lung open chromatin data did not provide sufficient resolution to infer cell types for most transcripts; direct validation of cell-specific transcripts remains desirable.
  • The paper notes technical limitations of small RNA sequencing, which has not yet been applied to single neurons and still relies on enrichment of certain cell types.
  • Evident from the text: sex-specific comparisons used an unbalanced sample (4 males vs 10 females), and tissue-level bulk profiling cannot fully resolve cell-type composition.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.1073/pnas.2002277117

Cites

140 of 140 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Figures

Fig. 1. Small ncRNA expression across mouse tissues. (A) Tissues and ncRNA classes profiled in the current study. Ten somatic tissues were collected from adult mice (n = 14). Testes were collected from male mice (n = 4). (B) ncRNAs identified in the current study. Numbers indicate detected and total annotated within GENCODE M20 miRNA, snoRNA, snRNA, scaRNA, Mt_tRNA, as well as high-confidence tRNA listed in GtRNAdb. (C) Coverage of ncRNA types within the profiled tissues. ncRNA was considered transcribed in a tissue if detected at >1 cpm. (D) Genomic map of small RNAs (sRNAs) expression across mouse genome. The bars show the log-transformed normalized expression count of ncRNAs. The red and gray bars around each circle represent the variance of each sRNA across 10 mouse tissues. Red denotes highly (the SD of expression above 25% of the mean value), and gray, low, variable ncRNAs (SD below 25% of the mean value).
Fig. 1. Small ncRNA expression across mouse tissues. (A) Tissues and ncRNA classes profiled in the current study. Ten somatic tissues were collected from adult mice (n = 14). Testes were collected from male mice (n = 4). (B) ncRNAs identified in the current study. Numbers indicate detected and total annotated within GENCODE M20 miRNA, snoRNA, snRNA, scaRNA, Mt_tRNA, as well as high-confidence tRNA listed in GtRNAdb. (C) Coverage of ncRNA types within the profiled tissues. ncRNA was considered transcribed in a tissue if detected at >1 cpm. (D) Genomic map of small RNAs (sRNAs) expression across mouse genome. The bars show the log-transformed normalized expression count of ncRNAs. The red and gray bars around each circle represent the variance of each sRNA across 10 mouse tissues. Red denotes highly (the SD of expression above 25% of the mean value), and gray, low, variable ncRNAs (SD below 25% of the mean value).CC BY-NC-ND · PubMed Central
Fig. 2. Tissue-specific patterns of small ncRNA expression. (A) t-SNE projection of ncRNA expression patterns performed on ∼4,000 ncRNA genes of various classes detected in 11 mouse tissues. (B) Dot plot of tissue-specific snoRNA, Rny1 and Terc, identified as the most tissue-specific (Benjamini–Hochberg [BH]-adjusted P value < 0.01 in LRT test). The size of the dot represents scaled log-transformed normalized counts. (C) Dot plot of tissue-specific miRNA. Only a subset of miRNA passing the specificity threshold (BH-adjusted P value < 0.01 in LRT test) is shown. “Arm” column denotes whether 5p-, 3p-, or both arms are passing the specificity threshold. (D) Levels of miRNA arms detected across tissues. Representative miRNAs, for which we consistently detect either one of the arms, both, or switched arm between tissues. The y axis represents normalized scaled counts. (E) Examples of ubiquitous ncRNA present in all tissues.
Fig. 2. Tissue-specific patterns of small ncRNA expression. (A) t-SNE projection of ncRNA expression patterns performed on ∼4,000 ncRNA genes of various classes detected in 11 mouse tissues. (B) Dot plot of tissue-specific snoRNA, Rny1 and Terc, identified as the most tissue-specific (Benjamini–Hochberg [BH]-adjusted P value < 0.01 in LRT test). The size of the dot represents scaled log-transformed normalized counts. (C) Dot plot of tissue-specific miRNA. Only a subset of miRNA passing the specificity threshold (BH-adjusted P value < 0.01 in LRT test) is shown. “Arm” column denotes whether 5p-, 3p-, or both arms are passing the specificity threshold. (D) Levels of miRNA arms detected across tissues. Representative miRNAs, for which we consistently detect either one of the arms, both, or switched arm between tissues. The y axis represents normalized scaled counts. (E) Examples of ubiquitous ncRNA present in all tissues.CC BY-NC-ND · PubMed Central
Fig. 3. miRNAs uniquely detected in the present study. (A) Pie chart of predicted and verified miRNA. (B) Examples of RT-qPCR verified miRNAs identified by this and previous studies. The levels of miRNA were determined from small RNA-seq data (DESeq2 normalized counts) and through RT-qPCR (Cq, adjusted for the sample-to-sample variability using cel-mir-39 spike in control). (C) Same as C but for miRNAs detected uniquely within the present study. (D) RT-qPCR verified miRNA, 14_6588, transcribed from the negative strand of mir-208a. (E) Tissue specificity of putative miRNA.
Fig. 3. miRNAs uniquely detected in the present study. (A) Pie chart of predicted and verified miRNA. (B) Examples of RT-qPCR verified miRNAs identified by this and previous studies. The levels of miRNA were determined from small RNA-seq data (DESeq2 normalized counts) and through RT-qPCR (Cq, adjusted for the sample-to-sample variability using cel-mir-39 spike in control). (C) Same as C but for miRNAs detected uniquely within the present study. (D) RT-qPCR verified miRNA, 14_6588, transcribed from the negative strand of mir-208a. (E) Tissue specificity of putative miRNA.CC BY-NC-ND · PubMed Central
Fig. 4. tRFs detected in mouse tissues. (A) Schematic depiction of tRNA cleavage and the resulting fragments. (B) Average tRF length identified across tissues for either ntRNA-derived (Top) or mtRNA-derived (Bottom) fragments. (C) Fragment type abundance across tissues. (D) Heatmap of tRF levels in tissues. For each of the 23 tRNA types, the sum among its tRFs is plotted. (E) Heatmap of relative abundance of ntRNA fragment types (5′-tR-halves, 3′-tR-halves, 5′-tRFs, 3′-tRFs, or 3′-CCA-tRFs) for each tRNA isoacceptor across 10 somatic tissues. Relative abundance is represented by row-wise scaled fractionated scores of tRFs computed by unitas. (F) Same as E but for mtRNA.
Fig. 4. tRFs detected in mouse tissues. (A) Schematic depiction of tRNA cleavage and the resulting fragments. (B) Average tRF length identified across tissues for either ntRNA-derived (Top) or mtRNA-derived (Bottom) fragments. (C) Fragment type abundance across tissues. (D) Heatmap of tRF levels in tissues. For each of the 23 tRNA types, the sum among its tRFs is plotted. (E) Heatmap of relative abundance of ntRNA fragment types (5′-tR-halves, 3′-tR-halves, 5′-tRFs, 3′-tRFs, or 3′-CCA-tRFs) for each tRNA isoacceptor across 10 somatic tissues. Relative abundance is represented by row-wise scaled fractionated scores of tRFs computed by unitas. (F) Same as E but for mtRNA.CC BY-NC-ND · PubMed Central
Fig. 5. Sex-specific miRNA expression. (A) Sex-dimorphic miRNAs identified in this study. The y axis represents the log10 of the difference between mean miRNA levels computed for males and females in each tissue (in counts per million [cpm]). Error bars denote SD of miRNA levels across a tissue within each sex. (B) Volcano plot showing miRNAs differentially expressed between the female and the male brain. (C) Expression and genomic location of male-biased miRNA in the brain.
Fig. 5. Sex-specific miRNA expression. (A) Sex-dimorphic miRNAs identified in this study. The y axis represents the log10 of the difference between mean miRNA levels computed for males and females in each tissue (in counts per million [cpm]). Error bars denote SD of miRNA levels across a tissue within each sex. (B) Volcano plot showing miRNAs differentially expressed between the female and the male brain. (C) Expression and genomic location of male-biased miRNA in the brain.CC BY-NC-ND · PubMed Central
Fig. 6. (A) Confusion matrices obtained from SVM tissue classifier for postnatal and embryonic datasets. (B) Correlation between small RNA-seq scores derived for top 400 tissue-specific ncRNA genes with the respective activity scores obtained from mouse ATAC atlas. (C–F) Average gene activity scores of tissue-specific ncRNAs within each cell type resident to the respective tissue. Gene activity scores for the ncRNAs of interest were retrieved from the mouse ATAC atlas. (G) Log10 gene activity scores of brain-specific miRNA across individual cells. (H) Gene activity scores computed for brain-specific ncRNA from the droplet-based scATAC-seq data (10XGenomics; Materials and Methods).
Fig. 6. (A) Confusion matrices obtained from SVM tissue classifier for postnatal and embryonic datasets. (B) Correlation between small RNA-seq scores derived for top 400 tissue-specific ncRNA genes with the respective activity scores obtained from mouse ATAC atlas. (C–F) Average gene activity scores of tissue-specific ncRNAs within each cell type resident to the respective tissue. Gene activity scores for the ncRNAs of interest were retrieved from the mouse ATAC atlas. (G) Log10 gene activity scores of brain-specific miRNA across individual cells. (H) Gene activity scores computed for brain-specific ncRNA from the droplet-based scATAC-seq data (10XGenomics; Materials and Methods).CC BY-NC-ND · PubMed Central

Verification

Claims the paper makes

107 Own result · 19 Partly supported · 14 Own method · 9 Own interpretation · 7 Supported · 5 No full text of the sources · 2 Supported so far · 2 Mixed support · 2 Couldn’t confirm · 1 Not supported · accuracy not yet measured

  • Probably not a substantial claim. 34
  • It may state two separate claims as one. 43
  • The quote only partly states this claim. 1
76 of 76
  • The paper relies on this without support.

    No prior study encompasses a spectrum of mammalian tissues from both female and male individuals, nor includes the other noncoding RNA types that have recently been identified to carry out tissue- and cell type-specific functions.

    “However, no prior study encompasses a spectrum of mammalian tissues from both female and male individuals, nor includes the other noncoding RNA types that have recently been identified to carry out tissue- and cell type-specific functions (26, 27).”
  • The paper relies on this without support.

    Surprisingly little is known about cell type-specific small ncRNA in the brain.

    “surprisingly little is known about cell type-specific small ncRNA in the brain.”
  • The paper relies on this without support.

    The brain remains one of the most underrepresented tissue within available tissue-level ncRNA resources.

    “the brain remains one of the most underrepresented tissue.”
  • The paper relies on this without support.

    Small ncRNAs have long been known to regulate the development and function of the brain.

    “Small ncRNAs have been long known to regulate the development and function of the brain.”
4 of 4

Self-check

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

Each contradiction between two passages listed here was confirmed by an independent check.

2 of 2

Citations

Cites

Show all 138 references

Cited by

Sorted by citations

See all 119

Metadata

Authors
Isakova, A., Fehlmann, T., Keller, A., & Quake, S. R.
Affiliations
  • Stanford University (US)
  • Saarland University (DE)
  • Stanford Medicine (US)
  • Chan Zuckerberg Biohub San Francisco (US)
Corresponding author
Quake, Stephen R.
Venue
Proceedings of the National Academy of Sciences
Published
2020-09-25
Volume
117
Issue
41
Pages
25634-25645
ISSN
  • 0027-8424 (ISSN-L)
  • 1091-6490
OpenAlex
W3088117294
Cited by
119
References
138
Citations by year
2026: 19 · 2025: 18 · 2024: 16 · 2023: 13 · 2022: 18 · 2021: 33 · 2020: 1
FWCI
5.19
Type
journal article
Language
en
Open access
Locations
  • Proceedings of the National Academy of Sciences (publishedVersion, cc-by-nc-nd)
  • bioRxiv (Cold Spring Harbor Laboratory) (acceptedVersion, cc-by-nc-nd)
  • PubMed Central (publishedVersion, cc-by-nc-nd)
  • PubMed (publishedVersion)
  • Publications of the UdS (Saarland University)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Cancer Research › MicroRNA in disease regulation
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Cancer Research › Cancer-related molecular mechanisms research
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › RNA modifications and cancer
Keywords
small RNA, tissue-specific gene expression, mouse tissues, microRNA, small nucleolar RNA, small nuclear RNA, tRNA-derived fragments, sexual dimorphism, single-cell chromatin accessibility, tissue atlas
MeSH terms
Animals, Brain, Female, Male, Mice, Inbred C57BL, Organ Specificity, Sex Characteristics, Mice, RNA, Small Untranslated, Single-Cell Analysis, Transcriptome

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.

A language model drew the claims from the paper on 2026-11-06, under version 3 of the claim definition (model deepseek-v4.1-flash, prompt extract-v3).

A second model labelled each claim on 2026-11-06 (model jev-1.13-20260917, definition revision gold-guideline-v0.2).

Cited sources are checked with a language model’s help (gpt-6-sol).

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