{"data":{"id":"7bab332b-5c5b-40ab-a2c6-ce7da66e3417","title":"From reads to genes to pathways: differential expression analysis of RNA-Seq experiments using Rsubread and the edgeR quasi-likelihood pipeline","article_type":"Original research","article_type_stored":"original_research","record_class":"imported_paper","record_class_reading":"A paper the journal holds from the corpus, with its canonical record and every version of it. The journal has not analysed it.","record_kind":"paper","record_status":"imported","abstract":"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.","authors":["Yunshun Chen","Aaron T. L. Lun","Gordon K. Smyth"],"doi":"10.12688\/f1000research.8987.2","claim_family_id":"family-sci-104bd5caa2b578c9","subject_category":"F1000Research","intake":{"method":"Imported as a canonical record","method_stored":"sci_import","source_name":"F1000Research","source_identifier":"pmid:27508061","source_url":null,"recorded_by":{"name":"Journal Editor"},"ingested_at":"2026-11-08T18:43:31+00:00"},"lineage":{"relationship":"Original","relationship_stored":"original","parent":null,"descendants":[{"id":"292d8109-0f49-46b6-aadb-938d572f8dc2","title":"From reads to genes to pathways: differential expression analysis of RNA-Seq experiments using Rsubread and the edgeR quasi-likelihood pipeline","relationship":"A fork of the published master it was opened from","relationship_stored":"fork","url":"https:\/\/www.openjournal.science\/journal\/submissions\/292d8109-0f49-46b6-aadb-938d572f8dc2"}]},"record":{"id":"pmid:27508061","versions":2,"head":{"revision":7,"content_hash":"sha256:72abc1c558c8d6a3461ec0365234046ee2dc5f723b2f8b3b3ae0904ac6e4ab85","document_sha256":"c85b4752faf9753a139e48667e4b37b73ed5929c8e54c1318286a40c4f33ad84","sci_version":"1.5","authority":"corpus_engine","authority_reading":"offered by the Corpus Engine","held_at":"2026-11-08T19:04:22+00:00"},"url":"https:\/\/www.openjournal.science\/journal\/versions\/pmid:27508061"},"merge_request":null,"assays":[{"id":"assay:1","role":"data_generation","depends_on":[],"definition":{"id":"rna-seq-read-alignment","version":"0.1.0","digest":"sha256:781215e3b5e462988c0f0a214b7462f32b5e18c58a59ff0e7366fdf7059c67a9","operations":["rna_seq_read_alignment"],"operations_malformed":false,"operations_report":null,"support_state":"unsupported"},"request":null,"binding":{"source_record_id":null,"source_content_hash":null,"source_version_held":false,"names_this_record":false,"quoted_source_hash":null,"note":"This entry names no record, so the journal cannot say what it was run against."},"outcome":null,"words":{"status":null,"summary":null,"interpretation":null,"limitations":[]},"stated":[],"counts":[],"verdicts":[]},{"id":"assay:2","role":"data_generation","depends_on":["assay:1"],"definition":{"id":"rna-seq-gene-quantification","version":"0.1.0","digest":"sha256:a7731347e0ac054f8756d31ba6f586cb68a0e381b1946dfec53ef72bcddb1407","operations":["gene_count_quantification"],"operations_malformed":false,"operations_report":null,"support_state":"unsupported"},"request":null,"binding":{"source_record_id":null,"source_content_hash":null,"source_version_held":false,"names_this_record":false,"quoted_source_hash":null,"note":"This entry names no record, so the journal cannot say what it was run against."},"outcome":null,"words":{"status":null,"summary":null,"interpretation":null,"limitations":[]},"stated":[],"counts":[],"verdicts":[]},{"id":"assay:3","role":"normalisation","depends_on":["assay:2"],"definition":{"id":"rna-seq-count-filtering-and-normalisation","version":"0.1.0","digest":"sha256:029f48150fdbe3973df58cc1d9a63ea34bd5ba35c319469ee1cfd1249101c0fe","operations":["normalise_counts","normalise_counts_tmm"],"operations_malformed":false,"operations_report":null,"support_state":"partial"},"request":null,"binding":{"source_record_id":null,"source_content_hash":null,"source_version_held":false,"names_this_record":false,"quoted_source_hash":null,"note":"This entry names no record, so the journal cannot say what it was run against."},"outcome":null,"words":{"status":null,"summary":null,"interpretation":null,"limitations":[]},"stated":[],"counts":[],"verdicts":[]},{"id":"assay:4","role":"analysis","depends_on":["assay:3"],"definition":{"id":"rna-seq-differential-expression","version":"0.1.0","digest":"sha256:dfda866165c72d6cd1517a23e85867cdaa0adc1df88f9977b642792bdcdc7278","operations":["differential_expression_fit","fdr_bh","differential_expression_quasi_likelihood"],"operations_malformed":false,"operations_report":null,"support_state":"partial"},"request":null,"binding":{"source_record_id":null,"source_content_hash":null,"source_version_held":false,"names_this_record":false,"quoted_source_hash":null,"note":"This entry names no record, so the journal cannot say what it was run against."},"outcome":null,"words":{"status":null,"summary":null,"interpretation":null,"limitations":[]},"stated":[],"counts":[],"verdicts":[]},{"id":"assay:5","role":"analysis","depends_on":["assay:4"],"definition":{"id":"rna-seq-pathway-enrichment","version":"0.1.0","digest":"sha256:13560e312926f7eb0829d4fcf69d3ebb6bc2af081c63ea341def7d3ca69a1e27","operations":["pathway_enrich"],"operations_malformed":false,"operations_report":null,"support_state":"supported"},"request":null,"binding":{"source_record_id":null,"source_content_hash":null,"source_version_held":false,"names_this_record":false,"quoted_source_hash":null,"note":"This entry names no record, so the journal cannot say what it was run against."},"outcome":null,"words":{"status":null,"summary":null,"interpretation":null,"limitations":[]},"stated":[],"counts":[],"verdicts":[]}],"data_availability_statement":null,"data_sources":[{"kind":"sci_record","name":"Corpus Engine","uri":"https:\/\/corpus.bntrg.co.uk\/sci\/1.0\/sci.schema.json","identifier":"pmid:27508061","version":"sci 1.5 revision 7","checksum_sha256":"sha256:72abc1c558c8d6a3461ec0365234046ee2dc5f723b2f8b3b3ae0904ac6e4ab85","access":"open","metadata":{"gaps":["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","execution.not_run","method.rna_seq_read_alignment","review.human_critique","review.independent_reader","review.machine_critique","step:run_analysis blocked on method:rna_seq_read_alignment","target.alignment_figures_not_comparable_across_runs","target.count_matrix_cells_not_published","target.fry_and_camera_results_not_published","target.gene_ontology_results_not_published","target.kegg_p_values_not_comparable_across_runs"],"layers":{"data":"partial","paper":"partial","analysis":"partial","readiness":"blocked","capability":"complete","experiment":"partial"},"revision":7,"authority":"corpus_engine","readiness":{"tier":"R4","state":"blocked","blocked_on":"method:rna_seq_read_alignment","blocking_axis":"method"},"schema_uri":"https:\/\/corpus.bntrg.co.uk\/sci\/1.0\/sci.schema.json","sci_version":"1.5","content_hash":"sha256:72abc1c558c8d6a3461ec0365234046ee2dc5f723b2f8b3b3ae0904ac6e4ab85","document_sha256":"c85b4752faf9753a139e48667e4b37b73ed5929c8e54c1318286a40c4f33ad84","document_public_id":"a57bcded-356c-450a-90fe-498e4f5f2d42"}},{"kind":"source_text","name":"Source text","uri":"https:\/\/doi.org\/10.12688\/f1000research.8987.2","identifier":"https:\/\/doi.org\/10.12688\/f1000research.8987.2","version":"","checksum_sha256":"sha256:3bf8a3f4c0135bde16e66f7f232937e6a89868783968461f69da3b9c8795dceb","access":"open","metadata":{"media_type":"application\/xml","content_hash":"sha256:72abc1c558c8d6a3461ec0365234046ee2dc5f723b2f8b3b3ae0904ac6e4ab85","sci_record_id":"pmid:27508061"}},{"kind":"source_text","name":"Source text","uri":"https:\/\/www.ncbi.nlm.nih.gov\/geo\/download\/?acc=GSE60450&format=file&file=GSE60450_Lactation-GenewiseCounts.txt.gz","identifier":"https:\/\/www.ncbi.nlm.nih.gov\/geo\/download\/?acc=GSE60450&format=file&file=GSE60450_Lactation-GenewiseCounts.txt.gz","version":"","checksum_sha256":"sha256:906b6ec4c33e95ca0ba454a3e0d07ad74e9cb7528ee38b89a9b122f5674a5f0d","access":"open","metadata":{"media_type":"text\/tab-values","content_hash":"sha256:72abc1c558c8d6a3461ec0365234046ee2dc5f723b2f8b3b3ae0904ac6e4ab85","sci_record_id":"pmid:27508061"}},{"kind":"source_text","name":"Source text","uri":"https:\/\/bioconductor.org\/packages\/3.3\/data\/annotation\/src\/contrib\/org.Mm.eg.db_3.3.0.tar.gz","identifier":"https:\/\/bioconductor.org\/packages\/3.3\/data\/annotation\/src\/contrib\/org.Mm.eg.db_3.3.0.tar.gz","version":"3.3.0","checksum_sha256":"sha256:5dc93707b7ca319b295b5b19963c573d5a3817d836f458a3e30d635ddd4e2365","access":"open","metadata":{"media_type":"application\/x-tar+gzip","content_hash":"sha256:72abc1c558c8d6a3461ec0365234046ee2dc5f723b2f8b3b3ae0904ac6e4ab85","sci_record_id":"pmid:27508061"}},{"kind":"source_text","name":"Source text","uri":"http:\/\/bioinf.wehi.edu.au\/edgeR\/F1000Research2016\/targets.txt","identifier":"http:\/\/bioinf.wehi.edu.au\/edgeR\/F1000Research2016\/targets.txt","version":"","checksum_sha256":"sha256:88ae42c0539403073f0bbbd53287dd9c9cf61ddd5b1a42856243d7d99113cb8b","access":"open","metadata":{"media_type":"text\/plain","content_hash":"sha256:72abc1c558c8d6a3461ec0365234046ee2dc5f723b2f8b3b3ae0904ac6e4ab85","sci_record_id":"pmid:27508061"}},{"kind":"source_text","name":"Source text","uri":"http:\/\/bioinf.wehi.edu.au\/edgeR\/F1000Research2016\/edgeRQL.R","identifier":"http:\/\/bioinf.wehi.edu.au\/edgeR\/F1000Research2016\/edgeRQL.R","version":"","checksum_sha256":"sha256:3b11478ac314c4456af8b4ba7af7992e938488ffc7bec88d3f73f42aed4c0bba","access":"open","metadata":{"media_type":"text\/x-r-source","content_hash":"sha256:72abc1c558c8d6a3461ec0365234046ee2dc5f723b2f8b3b3ae0904ac6e4ab85","sci_record_id":"pmid:27508061"}},{"kind":"source_text","name":"Source text","uri":"http:\/\/bioinf.wehi.edu.au\/edgeR\/F1000Research2016\/edgeRQL.Rnw","identifier":"http:\/\/bioinf.wehi.edu.au\/edgeR\/F1000Research2016\/edgeRQL.Rnw","version":"","checksum_sha256":"sha256:defd3acb73dcbfe885ca02209ac9437b785b4ccf12c7ba09a8a8bb311173902d","access":"open","metadata":{"media_type":"text\/x-tex","content_hash":"sha256:72abc1c558c8d6a3461ec0365234046ee2dc5f723b2f8b3b3ae0904ac6e4ab85","sci_record_id":"pmid:27508061"}},{"kind":"source_text","name":"Source text","uri":"http:\/\/bioinf.wehi.edu.au\/edgeR\/F1000Research2016\/edgeRQL.pdf","identifier":"http:\/\/bioinf.wehi.edu.au\/edgeR\/F1000Research2016\/edgeRQL.pdf","version":"","checksum_sha256":"sha256:76ea4da65b151b1bb78ba42b0d80c0e864632df0e13219fe968525a43efd2276","access":"open","metadata":{"media_type":"application\/pdf","content_hash":"sha256:72abc1c558c8d6a3461ec0365234046ee2dc5f723b2f8b3b3ae0904ac6e4ab85","sci_record_id":"pmid:27508061"}}],"reproducibility_stages":[],"peer_review_forks":[{"id":"3721b743-9bff-4d9e-869d-e409a3f28d2d","type":"replication","version":1,"target_result":"The paper's own analysis finds more than 5000 differentially expressed genes at a false discovery rate of 5% between the basal pregnant and lactating groups, and its printout shows the top gene Csn1s2b with a log2 fold change of 6.08.","control_outcome":"not_established","url":"https:\/\/www.openjournal.science\/journal\/forks\/3721b743-9bff-4d9e-869d-e409a3f28d2d","api_url":"https:\/\/www.openjournal.science\/api\/v1\/forks\/3721b743-9bff-4d9e-869d-e409a3f28d2d"},{"id":"572d1745-73b3-4990-81ed-415f67c216e6","type":"replication","version":1,"target_result":"Between basal lactating and basal pregnant mice, the edgeR quasi-likelihood F-test calls 5245 differentially expressed genes at a false discovery rate of 5%, and the top gene is Csn1s2b with a log2 fold change of 6.08.","control_outcome":"not_established","url":"https:\/\/www.openjournal.science\/journal\/forks\/572d1745-73b3-4990-81ed-415f67c216e6","api_url":"https:\/\/www.openjournal.science\/api\/v1\/forks\/572d1745-73b3-4990-81ed-415f67c216e6"}],"state":"submitted","revision":0,"version":1,"first_receipt_at":"2026-11-08T18:43:31+00:00","days_open":1,"admitted_at":null,"decided_at":null,"package_hash":"509bdc13bb77e0d3d0d3df7c3de51f5b6ec998ea74675aa7d20236577adbeb36","plan_hash":null,"plan_version":null,"claims":[],"result":null,"url":"https:\/\/www.openjournal.science\/journal\/submissions\/7bab332b-5c5b-40ab-a2c6-ce7da66e3417","api_url":"https:\/\/www.openjournal.science\/api\/v1\/submissions\/7bab332b-5c5b-40ab-a2c6-ce7da66e3417","verdict_api_url":"https:\/\/www.openjournal.science\/api\/v1\/verdicts\/7bab332b-5c5b-40ab-a2c6-ce7da66e3417"},"contract_version":"1.1.0","api_version":"v1"}