Journal of Open Source Science
Menu

fork:pmid:27508061:pmid-27508061-from-reads-to-genes-to-pathways-f804cb05 Public

https://scm.piranha.sh/git/ada56728063306072cb24f17/sci-fork-pmid-27508061-pmid-27508061-from-reads-to-genes-to-pathways-f804cb05 · at 65c2fc923fb2

From reads to genes to pathways: differential expression analysis of RNA-Seq experiments using Rsubread and the edgeR quasi-likelihood pipeline

fork:pmid:27508061:pmid-27508061-from-reads-to-genes-to-pathways-f804cb05 · revision 2 · F1000Research

Original
2016
Created
2026-11-04
Ingested
2026-11-07
Paper record The research question, publication details and source material recorded for this paper.

In recent years, RNA sequencing (RNA-seq) has become a very widely used technology for profiling gene expression. One of the most common aims of RNA-seq profiling is to identify genes or molecular pathways that are differentially expressed (DE) between two or more biological conditions. This article demonstrates a computational workflow for the detection of DE genes and pathways from RNA-seq data by providing a complete analysis of an RNA-seq experiment profiling epithelial cell subsets in the mouse mammary gland. The workflow uses R software packages from the open-source Bioconductor project and covers all steps of the analysis pipeline, including alignment of read sequences, data exploration, differential expression analysis, visualization and pathway analysis. Read alignment and count quantification is conducted using the Rsubread package and the statistical analyses are performed using the edgeR package. The differential expression analysis uses the quasi-likelihood functionality of edgeR.

Yunshun Chen, Aaron T. L. Lun, Gordon K. Smyth

DOI
10.12688/f1000research.8987.2
PMID
27508061
Journal
F1000Research
Published
2016
Revision
2
Format version
1.5
Content hash
sha256:89ea19103aaddafa68668bf0a0e136a6db62f9f7b5bf522b0c5b62d5b53e274a

Recomputed from the file: it matches.

Opened from
pmid:27508061
Fork of pmid:27508061
Written
2026-11-04T13:22:05.000000Z
Updated
2026-11-04T15:31:17.276166Z
Made by
Automated Assays Administrator, Automated Assays
Source run
assays-session-f804cb05-b3cd-4b58-a0a2-9608e8b5289e
Data sources The sources the file names, with the access and the state the file records for each one.

manuscript · open · resolved

10.12688/f1000research.8987.2

https://doi.org/10.12688/f1000research.8987.2

xml · 60,617 bytes · sha256:3bf8a3f4c0135bde16e66f7f232937e6a89868783968461f69da3b9c8795dceb

Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.

raw · open · resolved

GSE60450

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

geo-series ·

Except for the targets file targets.txt, all data analyzed in the workflow is read automatically from public websites as part of the code. All software used is publicly available as part of Bioconductor 3.3, except for the fastq-dump utility, which can be downloaded from NCBI website as described in the text. The article includes the complete code necessary to reproduce the analyses shown.

supplementary · open · verified

GSE60450_Lactation-GenewiseCounts.txt.gz

https://www.ncbi.nlm.nih.gov/geo/download/?acc=GSE60450&format=file&file=GSE60450_Lactation-GenewiseCounts.txt.gz

geo-supplementary · 513,592 bytes · sha256:906b6ec4c33e95ca0ba454a3e0d07ad74e9cb7528ee38b89a9b122f5674a5f0d

The genewise read-count matrix the paper analyses: 27179 genes over twelve libraries, named by the repository. Fetched and hashed on 4 November 2026, and checked against the numbers the paper printed.

reference · open · verified

org.Mm.eg.db_3.3.0

https://bioconductor.org/packages/3.3/data/annotation/src/contrib/org.Mm.eg.db_3.3.0.tar.gz

bioconductor-package · 63,495,324 bytes · sha256:5dc93707b7ca319b295b5b19963c573d5a3817d836f458a3e30d635ddd4e2365

The annotation release the paper's own sessionInfo names, fetched and hashed on 4 November 2026 from the Bioconductor 3.3 archive. It decides which genes keep a symbol, and the paper discards the ones that do not, so it is an input to the filtering stage.

code · open · verified

targets.txt

http://bioinf.wehi.edu.au/edgeR/F1000Research2016/targets.txt

text/plain · 525 bytes · sha256:88ae42c0539403073f0bbbd53287dd9c9cf61ddd5b1a42856243d7d99113cb8b

The authors' own file, fetched and hashed on 3 November 2026. The publisher states no licence for the entry.

code · open · verified

edgeRQL.R

http://bioinf.wehi.edu.au/edgeR/F1000Research2016/edgeRQL.R

R source · 10,298 bytes · sha256:3b11478ac314c4456af8b4ba7af7992e938488ffc7bec88d3f73f42aed4c0bba

The authors' own file, fetched and hashed on 3 November 2026. The publisher states no licence for the entry.

code · open · verified

edgeRQL.Rnw

http://bioinf.wehi.edu.au/edgeR/F1000Research2016/edgeRQL.Rnw

Rnw source · 59,373 bytes · sha256:defd3acb73dcbfe885ca02209ac9437b785b4ccf12c7ba09a8a8bb311173902d

The authors' own file, fetched and hashed on 3 November 2026. The publisher states no licence for the entry.

code · open · verified

edgeRQL.pdf

http://bioinf.wehi.edu.au/edgeR/F1000Research2016/edgeRQL.pdf

PDF · 904,232 bytes · sha256:76ea4da65b151b1bb78ba42b0d80c0e864632df0e13219fe968525a43efd2276

The authors' own file, fetched and hashed on 3 November 2026. The publisher states no licence for the entry.

reference · open · declared

GCA_000001635.6_GRCm38.p4

fasta.gz ·

The mouse genome build the paper indexes and aligns to, as the paper names it.

The samples analysed The group sizes the file holds, and the group sizes the paper itself states. A run is only comparable when it used the same cohort.

In this file

Basal 18.5 dP
2
Basal 2 dL
2
Basal virgin
2
Luminal 18.5 dP
2
Luminal 2 dL
2
Luminal virgin
2

As the paper states them

B.lactating
2
B.pregnant
2
B.virgin
2
L.lactating
2
L.pregnant
2
L.virgin
2

12 samples in total.

Reference group: B.pregnant.

There are therefore six groups of RNA samples, one for each combination of cell type and mouse status. Two biological replicates were collected for each group.

Reproducibility map Follow the research stages and see which inputs each stage needs and what state the file records for it.

A fork of the exact Corpus revision, containing the five paper-defined assays in order. No assay has run until a recorded execution is bound to this fork.

  1. 1

    Read the inputs and check them against the hashes the paper and its sources publish

    complete

    computational

  2. 2

    Run the calls the authors published, in the order the authors published them

    complete

    computational

  3. 3

    Compare every number with the number the paper printed, at the precision it printed it

    complete

    computational

  4. 4

    Write the tables and the figures, and hash each file the assay produced

    complete

    computational

  5. 5

    Package the result with the dataset hash, the parameters and the artefact hashes

    complete

    computational

  6. 6

    Read the inputs and check them against the hashes the paper and its sources publish

    complete

    computational

  7. 7

    Run the calls the authors published, in the order the authors published them

    complete

    computational

  8. 8

    Compare every number with the number the paper printed, at the precision it printed it

    complete

    computational

  9. 9

    Write the tables and the figures, and hash each file the assay produced

    complete

    computational

  10. 10

    Package the result with the dataset hash, the parameters and the artefact hashes

    complete

    computational

  11. 11

    Read the inputs and check them against the hashes the paper and its sources publish

    complete

    computational

  12. 12

    Run the calls the authors published, in the order the authors published them

    complete

    computational

  13. 13

    Compare every number with the number the paper printed, at the precision it printed it

    complete

    computational

  14. 14

    Write the tables and the figures, and hash each file the assay produced

    complete

    computational

  15. 15

    Package the result with the dataset hash, the parameters and the artefact hashes

    complete

    computational

  16. 16

    Read the inputs and check them against the hashes the paper and its sources publish

    complete

    computational

  17. 17

    Run the calls the authors published, in the order the authors published them

    complete

    computational

  18. 18

    Compare every number with the number the paper printed, at the precision it printed it

    complete

    computational

  19. 19

    Write the tables and the figures, and hash each file the assay produced

    complete

    computational

  20. 20

    Package the result with the dataset hash, the parameters and the artefact hashes

    complete

    computational

  21. 21

    Read the inputs and check them against the hashes the paper and its sources publish

    complete

    computational

  22. 22

    Run the calls the authors published, in the order the authors published them

    complete

    computational

  23. 23

    Compare every number with the number the paper printed, at the precision it printed it

    complete

    computational

  24. 24

    Write the tables and the figures, and hash each file the assay produced

    complete

    computational

  25. 25

    Package the result with the dataset hash, the parameters and the artefact hashes

    complete

    computational

The methods the paper publishes Each method as the paper states it, with the evidence the file holds for it.

Read alignment

expressible

sequence_align

libraries
{"type":"array","description":"The libraries this run aligns. The twelve the targets file names, when the request names none. A subset is a declared change to the run and is recorded as one.","required":false}
Read alignment as the workflow of the paper runs it.
The assays in this file Each assay the file names, with the operations it needs and the definition it was read from.

rna-seq-read-alignment

supported

assay:1 · data_generation · version 0.1.0

needs: rna_seq_read_alignment, sequence_align

sha256:781215e3b5e462988c0f0a214b7462f32b5e18c58a59ff0e7366fdf7059c67a9

rna-seq-gene-quantification

supported

assay:2 · data_generation · version 0.1.0

needs: rna_seq_gene_quantification

depends on: assay:1

sha256:a7731347e0ac054f8756d31ba6f586cb68a0e381b1946dfec53ef72bcddb1407

rna-seq-count-filtering-and-normalisation

supported

assay:3 · normalisation · version 0.1.0

needs: rna_seq_count_filtering_and_normalisation, normalise_counts

depends on: assay:2

sha256:029f48150fdbe3973df58cc1d9a63ea34bd5ba35c319469ee1cfd1249101c0fe

rna-seq-differential-expression

supported

assay:4 · analysis · version 0.1.0

needs: rna_seq_differential_expression, differential_expression_fit

depends on: assay:3

sha256:dfda866165c72d6cd1517a23e85867cdaa0adc1df88f9977b642792bdcdc7278

rna-seq-pathway-enrichment

supported

assay:5 · analysis · version 0.1.0

needs: rna_seq_pathway_enrichment, pathway_enrich

depends on: assay:4

sha256:13560e312926f7eb0829d4fcf69d3ebb6bc2af081c63ea341def7d3ca69a1e27

Can this run here What the file says about running this experiment, and what still stands in the way.
State
blocked · tier R4
Blocked on
method:rna_seq_read_alignment

can run: normalise_counts, differential_expression_fit, fdr_bh, pathway_enrich

Still to build

  • rna_seq_read_alignment
  • gene_count_quantification
  • normalise_counts_tmm
  • differential_expression_quasi_likelihood

method

The DSL catalogue has no operation for: rna_seq_read_alignment, gene_count_quantification, normalise_counts_tmm, differential_expression_quasi_likelihood. The catalogue's normalise_counts is counts per million and its differential_expression_fit is a Welch t-test, so neither is edgeR's TMM normalisation nor its quasi-likelihood F-test.

capability

The corpus capability table owns computational analysis.

data

The GEO series resolves, its twelve samples are listed, and the twelve SRA runs resolve. No dataset was downloaded and no licence was read.

The result run What has run against this file, and the review the file carries.

Runs

85f9fdf4-4c2e-4155-8daa-da11d28021ef

rna_seq_differential_expression · complete

ran against
pmid:27508061
version
0.1.0
provider
automated-assays-computational
result
sha256:0da96b2fe1f4ab61ad9ec91f1a9130268fb314dd354d6b6f623388f5e7ae9d4c
evidence
899e82b58fd83ec21bfd48461d0797f7867bbd5dedf5744ed19b8654b456e291

cb681050-44e5-43c9-801c-87e82a09ed55

rna_seq_count_filtering_and_normalisation · complete

ran against
pmid:27508061
version
0.1.0
provider
automated-assays-computational
result
sha256:b681c5238c00d87068957502cabb988094a393818da50840119649867a1fb4df
evidence
f7288a5fdb3282fea31889580c03db42473612241a6fed77527aa211540c6df3

8573fa99-85cd-4f56-bdd5-3651d1cdfd58

rna_seq_gene_quantification · refused

ran against
pmid:27508061
version
0.1.0
provider
automated-assays-computational
result
sha256:43cac33f2a7f97f767d87937964e805d3ae6fab5b695818b5cf7d37efdb790fb
evidence
0e116046745cbd0ed15ae537d918f4776dfa6563908c644dbdf2221600f760ed

7daff85e-0e95-4e29-b922-1bf66aa96252

rna_seq_read_alignment · refused

ran against
pmid:27508061
version
0.1.0
provider
automated-assays-computational
result
sha256:be55ead2068cf19cc9efd01d113aa3057e4b3dde2343ea0ecd4562aaf2d1bd0e
evidence
48eb942441653112489a7ce4ff0821a719728c1188aeab486161fb93c17ec872

213bfb92-8ff7-4d1e-925c-19d9f7047551

rna_seq_count_filtering_and_normalisation · refused

ran against
pmid:27508061
version
0.1.0
provider
automated-assays-computational

c0d08ec6-0c1f-49db-b139-adcc2916656d

rna_seq_count_filtering_and_normalisation · refused

ran against
pmid:27508061
version
0.1.0
provider
automated-assays-computational

91daa9a9-2b77-45e0-8727-13a7643efe9d

rna_seq_pathway_enrichment · complete

ran against
pmid:27508061
version
0.1.0
provider
automated-assays-computational
result
sha256:6daad652773fe321cd885609fcc97c1378b769fa9173daaf172c33582fb30efb
evidence
d1d756c2c34777bc09aa0d236928ddd951666073cec97e52acd028342d3ed85f

a92cdcdd-e95c-4579-a4dd-7648d2c24eff

rna_seq_differential_expression · complete

ran against
pmid:27508061
version
0.1.0
provider
automated-assays-computational
result
sha256:194ba83e9468a04462de791795cd0bdf4bbdc992bb5cea0b420fa659fe8ecc1f
evidence
ffd49280b7bcfa17e938f077569600a4d38d3dce873176d5f4f8f5aa230659da

fa44f761-7022-4c3f-975f-d8ce85abd2a8

rna_seq_count_filtering_and_normalisation · complete

ran against
pmid:27508061
version
0.1.0
provider
automated-assays-computational
result
sha256:a1f4df7f4ba128a759d8da6db62614641875603a7474e4c02572b23b18563976
evidence
25d0e73c4e732aef23de129367c22fb8c3889275e672ff9c8be9c63fdd553a59

Review carried by the file

absent

No person has reviewed this file. The reader's assessment is an opinion and not a decision.

Scope and limitations What this file does not cover, as the file itself states it.
  • execution.unbound: 9 runs against this record are not bound to a version of it, so are not published as evidence on an assay entry.

data

  • asset:data:geo:GSE60450 licence not read
  • asset:data:geo:GSE60450 not downloaded
  • asset:data:sra:SRP045534 licence not read
  • asset:data:sra:SRP045534 not downloaded
  • asset:reference:mm10 not downloaded
  • data.licence_not_stated_in_the_repository_record

readiness

  • method:rna_seq_read_alignment

review

  • review.absent: this document is a working version. Review and publication are separate actions, and they belong to the Journal.
Forks opened from this paper Each fork is a repository of its own, with the paper it was opened from.

No fork has been opened from this paper yet.

Files

fork:pmid:27508061:pmid-27508061-from-reads-to-genes-to-pathways-f804cb05.json 90,559 bytes