ARMT is auto RNA-seq data minning tool.
The aim of ARMT is :
i) integrate and facilitate the downstream analysis of RNA-seq data,
ii) provide the way to analyze geneset according to GSVA,
iii) explore the correlation of gene expression and mutation.
Author:
Guanda Huang —— 202010108432@mail.scut.edu.cn
Hongli Du —— hldu@scut.edu.cn
Year: 2021
License: GPL(>=3)
maftools
GSVA
GSVAdata
limma
GSEAbase
org.Hs.eg.db
edgeR
survival
clusterProfiler
DOSE
Provide TCGA clinical data
Create .gmt file for arbitary geneset
Normalization of counts matrix (TPM)
Gene set variant analysis
Survival analysis
Differential analysis
Correlation analysis
Analysis for multiple sets of data
Plot the mutation information
Easy to use
Automation
Visualization
Comprehensiveness
Integrating GSVA score analysis
R >=4.0.2, Rstudio, R packages(devtools or remotes)
The other dependencies are automatically installed. By default, outdated dependencies are automatically upgraded. In interactive sessions you can select a subset of the dependencies to upgrade.
To install this package from Github, please, use the code below.
if (!requireNamespace("devtools", quietly = TRUE))
install.packages("devtools")
devtools::install_github("Dulab2020/ARMT")
The following commands should be used in order to start the graphical user interface(GUI).
ARMT::run_app()
The GUI is developed based on 'shiny' package, and has four page including: Data, Normalization&GSVA, Integration&Analysis, Mutant mapping.
Figure1. GUI of ARMT: The GUI consists of two parts, sidebar panel at left and main panel at right. The sidebar panel is used to input data and adjust parameters; the main panel is used to demonstrate and save the results.
This page of ARMT provides TCGA clinical data and builds .gmt file of arbitrary geneset.
Through the internal datasets of ARMT, you can get TCGA clinical data of 33 cancer types.(Figure1)
To build a .gmt file, a .csv file should be provided, and it must contain two column as shown in Figure2.
Figure2. The .csv file to create .gmt file: This file has two columns. The first column contains genesets. The second one contains genes in corresponding geneset.
In this page, ARMT automatically normalize the expression counts matrix to TPM matrix, and use log2(TPM) to process gene set variant analysis(GSVA). The input data must be gene expression matrix with Ensemble ID, and the ID can be transfered to Symbole through normalization.
Figure3. Normalization&GSVA: There are 9 geneset from MSigDb in ARMT internal datasets available for direct selection to process GSVA, and it is also supported to input a .gmt file of arbitrary geneset from 'Data' page. If the input data is from TCGA, you can chose normal, tumor or all samples to process normalize and GSVA(default: tumor).
ARMT can integrate clinical, GSVA, gene expression(TPM), gene mutation data together, and use these data to carry out analysis including: survival analysis&cox proportional hazards regression analysis, differential analysis, enrichment analysis and correlation analysis.
The integrated data should contain the common samples. You should input the clinical, GSVA, TPM data by .csv file which can be produced by 'Data' and 'Integration&Analysis' page, and the mutation data should be entered by .maf file. The next analysis is carried out according to these integrated data. You can choose a column in integrated data table as the 'Group Column' to separate the data into multiple group to analyze respectively. Also, you can input your own data in the same format.
Figure4. Data Integration: In this page, clinical data and GSVA data will be merged automatically, and you can enter your interested gene name in your TPM and mutation data to integrate their expression and mutation information by clicking 'Integrate TPM' and 'Integrate Maf' button. The integrated data is showed in main panel.
In this page, you should choose two columns in integrated data table as the survival time and status of samples. Each column of integrated data table can be seem as a factor of samples to carried out survival analysis and cox proportional hazards regression analysis.
Figure5. Survival Analysis: If the survival factor is a continuous variable, the samples will be grouped into two parts(high&low). You can choose multiple factors to obtain multiple analysis results(K-M curve).
Figure6. Cox Proportional Hazards Regression Analysis: The cox analysis can focus on single-variable or multiple-variable. If the factor in cox analysis is not numeric variable, one value of this variable should be seem as a reference.
ARMT can carry out differential analysis to GSVA score, gene mutation and expression.
The difference factor should be selected in the integrated data, and the samples will be grouped by it. If the difference factor is a numeric variable, the samples will be grouped into high and low, and the low part will be seem as control group. The group cut off value 't' is a adjustable parameter in ARMT(upper and lower t*100%). If the difference factor is not a numeric variable, you should select two values used to group the samples into experiment group and control group.
You can screen the result by p-value, logFC and FDR(adj.p) after analysis and visualize it in main panel.
Figure7. The differential analysis of GSVA score: To analyze difference of GSVA score, the GSVA result of samples in integrated data should be entered through .csv file, and you can get this file in 'Normalization&GSVA' page.
Figure8. The differential analysis of gene mutation: To analyze difference of gene mutation, the mutation information of samples in integrated data should be entered through .maf file. The result is demonstrated with forest plot in main panel.
Figure9. The differential analysis of gene expression: To analyze difference of gene expression(DEG), the counts matrix of samples in integraated data should be entered through .csv file. The result is demonstrated with volcano plot in main panel.
The enrichment analysis ways in ARMT include GO and KEGG. This enrichment is processed for differential expression gene(DEG) from above differential analysis.
Figure10. Enrichment Analysis: You should set up the p-value and q-value cut off before enrichment analysis, and you can select your interested ontology(MF, CC, BP) of GO enrichment. The enriched genes are up regulation and down regulation in DEG, and their enrichment results(up regulation, down regulation and all DEG) are showed in different pages of main panel. These results are demonstrated with bar plot and dot plot.
ARMT can calculate the correlation between any continuous variable factors of integration data, such as TPM and GSVA score. The result is demonstrated in main panel with correlation coefficient matrix and heat map, and it can be descreened by p-value and correlation coefficient(r, Spearman or Pearson).
Figure11. ARMT can calculate correlation coefficient of each pair-factors.
Figure12. ARMT can calculate correlation coefficient between two list of factors.
If the integrated data is grouped by 'Group Column', ARMT can calculate correlation in each group of data by using single factor or all factors in the list 'Factor1'.
Figure13. For each group of data, the correlation coefficient between single factor and 'Factor2' list is put together in one result, and the single factor name is replaced by group name in correlation matrix. It is a case of pan-cancer correlation analysis.
Figure14. When all factors in 'Factor1' is demonstrated, you can choose one result of multiple groups in sidebar panel to show in main panel.
In corporation with 'maftools', ARMT can visualize the gene mutation in .maf file. There are two plot modes: Maftools Summary and Self-defined. The number of top mutant genes plot out can be set in sidebar panel.
Figure15. Maftools Summary: ARMT can plot the maf summary and illustrate the enrichment of known Oncogenic Signaling Pathways in TCGA cohorts. It is also supported to draw the oncoplot of interested pathway completely.
In Self-defined mode, ARMT can plot the map for any genes or genesets, and the mutation types (deletion, insert, SNV) can be specified.Figure16
Figure16. Plot for genes: About the specific genes, ARMT can plot out the VAF, and...
Figure17. Plot for genesets: Input the gmt file, ARMT shows the enrichment of mutation, and can plot a oncoplot of specific genesets.
The input data of ARMT include geneset, gene expression matrix, gene mutation, GSVA score matrix and clinical information.
In 'Data' page, the geneset file(.csv, Figure3) is used to create .gmt file, and the .gmt file is the main format of geneset input.
For gene expression, the counts matrix can be normalized to TPM matrix in ARMT.
Figure18. Counts matrix must be in a .csv file. The row name is gene Ensembl ID and the column name is the sample ID. It is used to normalization and differential analysis.
Figure19. TPM matrix must be in a .csv file. The row names is gene symbol ID and the column names is the sample ID. It can be integrated with other data.
To get the mutation information, ARMT requires standard .maf file like the mutation annotation format file in TCGA.
The GSVA score matrix can be obtained in 'Normalization&GSVA' page.
Figure20. The GSVA score is saved in a matrix of .csv file. The row name is the sample ID, and the column name is geneset name.
The TCGA clinical information can be obtained in 'Data' page.
Figure21. The clinical data must be in a .csv file, and the row name is sample ID and column name is sample characteristic. Any information about samples in this format file can be input as clinical data.
In 'example' folder, we have uploaded some example file to illustrate the input format of above data.