A Clustering Procedure for Three-way RNA Sequencing Data Using Data Transformations and Matrix-variate Gaussian Mixture Models

Overview

Journal BMC Bioinformatics

Publisher Biomed Central

Specialty Biology

Date 2024 Mar 1

PMID 38429687

Authors

Theresa Scharl

Bettina Grun

Affiliations

Soon will be listed here.

Abstract

RNA sequencing of time-course experiments results in three-way count data where the dimensions are the genes, the time points and the biological units. Clustering RNA-seq data allows to extract groups of co-expressed genes over time. After standardisation, the normalised counts of individual genes across time points and biological units have similar properties as compositional data. We propose the following procedure to suitably cluster three-way RNA-seq data: (1) pre-process the RNA-seq data by calculating the normalised expression profiles, (2) transform the data using the additive log ratio transform to map the composition in the D-part Aitchison simplex to a -dimensional Euclidean vector, (3) cluster the transformed RNA-seq data using matrix-variate Gaussian mixture models and (4) assess the quality of the overall cluster solution and of individual clusters based on cluster separation in the transformed space using density-based silhouette information and on compactness of the cluster in the original space using cluster maps as a suitable visualisation. The proposed procedure is illustrated on RNA-seq data from fission yeast and results are also compared to an analogous two-way approach after flattening out the biological units.

References

Love M, Anders S, Kim V, Huber W . RNA-Seq workflow: gene-level exploratory analysis and differential expression. F1000Res. 2015; 4:1070. PMC: 4670015. DOI: 10.12688/f1000research.7035.1. View

Pontes B, Giraldez R, Aguilar-Ruiz J . Biclustering on expression data: A review. J Biomed Inform. 2015; 57:163-80. DOI: 10.1016/j.jbi.2015.06.028. View

Bourgon R, Gentleman R, Huber W . Independent filtering increases detection power for high-throughput experiments. Proc Natl Acad Sci U S A. 2010; 107(21):9546-51. PMC: 2906865. DOI: 10.1073/pnas.0914005107. View

Love M, Huber W, Anders S . Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15(12):550. PMC: 4302049. DOI: 10.1186/s13059-014-0550-8. View

Nueda M, Tarazona S, Conesa A . Next maSigPro: updating maSigPro bioconductor package for RNA-seq time series. Bioinformatics. 2014; 30(18):2598-602. PMC: 4155246. DOI: 10.1093/bioinformatics/btu333. View

Rau A, Maugis-Rabusseau C . Transformation and model choice for RNA-seq co-expression analysis. Brief Bioinform. 2017; 19(3):425-436. DOI: 10.1093/bib/bbw128. View

Scharl T, Voglhuber I, Leisch F . Exploratory and inferential analysis of gene cluster neighborhood graphs. BMC Bioinformatics. 2009; 10:288. PMC: 2761902. DOI: 10.1186/1471-2105-10-288. View

Androulakis I, Yang E, Almon R . Analysis of time-series gene expression data: methods, challenges, and opportunities. Annu Rev Biomed Eng. 2007; 9:205-28. PMC: 4181347. DOI: 10.1146/annurev.bioeng.9.060906.151904. View

Anders S, Huber W . Differential expression analysis for sequence count data. Genome Biol. 2010; 11(10):R106. PMC: 3218662. DOI: 10.1186/gb-2010-11-10-r106. View

10.

Robinson M, McCarthy D, Smyth G . edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2009; 26(1):139-40. PMC: 2796818. DOI: 10.1093/bioinformatics/btp616. View

11.

Si Y, Liu P, Li P, Brutnell T . Model-based clustering for RNA-seq data. Bioinformatics. 2013; 30(2):197-205. DOI: 10.1093/bioinformatics/btt632. View

12.

Huber W, Carey V, Gentleman R, Anders S, Carlson M, Carvalho B . Orchestrating high-throughput genomic analysis with Bioconductor. Nat Methods. 2015; 12(2):115-21. PMC: 4509590. DOI: 10.1038/nmeth.3252. View

13.

Srivastava H, Ferrell D, Popescu G . NetSeekR: a network analysis pipeline for RNA-Seq time series data. BMC Bioinformatics. 2022; 23(1):54. PMC: 8796424. DOI: 10.1186/s12859-021-04554-1. View

14.

Harris M, Rutherford K, Hayles J, Lock A, Bahler J, Oliver S . Fission stories: using PomBase to understand Schizosaccharomyces pombe biology. Genetics. 2022; 220(4). PMC: 9209812. DOI: 10.1093/genetics/iyab222. View

15.

Amar D, Yekutieli D, Maron-Katz A, Hendler T, Shamir R . A hierarchical Bayesian model for flexible module discovery in three-way time-series data. Bioinformatics. 2015; 31(12):i17-26. PMC: 4765869. DOI: 10.1093/bioinformatics/btv228. View

16.

Scrucca L, Fop M, Murphy T, Raftery A . mclust 5: Clustering, Classification and Density Estimation Using Gaussian Finite Mixture Models. R J. 2016; 8(1):289-317. PMC: 5096736. View

17.

Sun J, Nishiyama T, Shimizu K, Kadota K . TCC: an R package for comparing tag count data with robust normalization strategies. BMC Bioinformatics. 2013; 14:219. PMC: 3716788. DOI: 10.1186/1471-2105-14-219. View

18.

Leong H, Dawson K, Wirth C, Li Y, Connolly Y, Smith D . A global non-coding RNA system modulates fission yeast protein levels in response to stress. Nat Commun. 2014; 5:3947. PMC: 4050258. DOI: 10.1038/ncomms4947. View

19.

Russo M, Singer B, Dunson D . MULTIVARIATE MIXED MEMBERSHIP MODELING: INFERRING DOMAIN-SPECIFIC RISK PROFILES. Ann Appl Stat. 2022; 16(1):391-413. PMC: 9222983. DOI: 10.1214/21-aoas1496. View

20.

Jung I, Jo K, Kang H, Ahn H, Yu Y, Kim S . TimesVector: a vectorized clustering approach to the analysis of time series transcriptome data from multiple phenotypes. Bioinformatics. 2017; 33(23):3827-3835. DOI: 10.1093/bioinformatics/btw780. View