Robustness of differential gene expression analysis of RNA-seq.

Loading...
Thumbnail Image

Embargo End Date

Authors

Stupnikov, A
McInerney, CE
Savage, KI
McIntosh, SA
Emmert-Streib, F
Kennedy, R
Salto-Tellez, M
Prise, KM
McArt, DG

Document Type

Journal Article

Date

2021-01-01

Date Accepted

2021-05-25

Abstract

RNA-sequencing (RNA-seq) is a relatively new technology that lacks standardisation. RNA-seq can be used for Differential Gene Expression (DGE) analysis, however, no consensus exists as to which methodology ensures robust and reproducible results. Indeed, it is broadly acknowledged that DGE methods provide disparate results. Despite obstacles, RNA-seq assays are in advanced development for clinical use but further optimisation will be needed. Herein, five DGE models (DESeq2, voom + limma, edgeR, EBSeq, NOISeq) for gene-level detection were investigated for robustness to sequencing alterations using a controlled analysis of fixed count matrices. Two breast cancer datasets were analysed with full and reduced sample sizes. DGE model robustness was compared between filtering regimes and for different expression levels (high, low) using unbiased metrics. Test sensitivity estimated as relative False Discovery Rate (FDR), concordance between model outputs and comparisons of a 'population' of slopes of relative FDRs across different library sizes, generated using linear regressions, were examined. Patterns of relative DGE model robustness proved dataset-agnostic and reliable for drawing conclusions when sample sizes were sufficiently large. Overall, the non-parametric method NOISeq was the most robust followed by edgeR, voom, EBSeq and DESeq2. Our rigorous appraisal provides information for method selection for molecular diagnostics. Metrics may prove useful towards improving the standardisation of RNA-seq for precision medicine.

Citation

Computational and structural biotechnology journal, 2021, 19 pp. 3470 - 3481

Source Title

Publisher

ELSEVIER

ISSN

2001-0370

eISSN

2001-0370

Collections

Research Team

Integrated Pathology
Integrated Pathology

Notes