Contents

1 Installation and Package Loading

The day-to-day development version from the Github repository can be installed.

# From Github 
if (!requireNamespace("devtools", quietly=TRUE))
  install.packages("devtools")
devtools::install_github("robogeno/SVMDO")
# From Biocodunctor
if(!requireNamespace("BiocManager", quietly = TRUE))
    install.packages("BiocManager")
BiocManager::install("SVMDO")
#library(SVMDO)
#Main screen of GUI
#runGUI()
# OR
#SVMDO::runGUI()

2 Introduction

Transcriptome-based supervised classification can be difficult to handle due to including vast amount of junk features reducing the efficiency of the process. To minimize the efficiency problem, feature selection methods have been applied [1]. Ensemble methodology in constructing a feature selection method has been mainly preferred for its better accuracy in classification processes [2]. Cancer-based studies have been reported to co-occur with several chronic diseases [3] and oncogenic viral infections [4]. Disease Ontology (DO) enrichment analysis can provide a filtration approach through reported diseases that interact with the gene sets used [5]. Wilk’s lambda criterion is one of the moderate techniques for feature selection. The method combines different features(genes) step by step based on their contributions to the discriminating power of overall model [6]. Thus, disease related genes with better discriminatory characteristics can be acquired. SVMDO is an easy-to-use GUI using disease information for detecting tumor/normal sample discriminating gene sets from differentially expressed genes. Our approach is based on an iterative algorithm filtering genes with disease ontology enrichment analysis and wilk’s lambda criterion connected to SVM classification model construction. Along with gene set extraction, SVMDO also provides individual prognostic marker detection. The algorithm is designed for FPKM and RPKM normalized RNA-Seq transcriptome datasets. During the development our algorithm, Bioconductor provided a robust approach for acquiring significant disease-gene interaction information from disease ontology database.

3 Implementation

SVMDO was developed by using Shiny package. It is available for Windows and Linux from Bioconductor website. SVMDO requires the following R packages: Shiny, ShinyFiles, golem, nortest, e1071, BSDA, data.table, sjmisc, klaR, caTools, dplyr, caret, survival, DOSE, AnnotationDbi, DOSE, org.Hs.eg.db

4 Dataset Preparation

RNA-Seq cancer transcriptome and clinical datasets must be prepared before applying them in SVMDO. The datasets must be used in txt format. The expected designs of input datasets are indicated in Table-1 and Table-2 for transcriptome and clinical datasets respectively.


Table 1: Structure of RNA-Seq Transcriptome Dataset
id tissue_type AB1G A2M
TCGA-AA-3662 Normal Exp_1 Exp_1
TCGA-AA-3514 Normal Exp_2 Exp_2
TCGA-D5-6541 Tumour Exp_3 Exp_3

Table 2: Structure of Clinical Dataset
id days_to_death vital_status
TCGA-AA-3662 49 Alive
TCGA-AA-3514 1331 Dead
TCGA-D5-6541 225 Dead

Using tissue_type and id as column names for representing tissue information and TCGA sample id is mandatory. If there is not any requirement for survival analysis, preparation of clinical dataset and involvement of id column in the transcriptome dataset are optional.

5 SVMDO GUI Description

The main dialog box is indicated in Figure-1. The GUI page consists of two main sections which are “Analysis” and “Result”. In the Analysis section, steps of acquiring discriminative gene set and preparations of individual survival plots are included. In the Result section, visualization and download of discriminative gene set and survival plots are included. At each step of the GUI sections, necessary variables are saved as objects in the workspace environment to be used in the following steps. To provide experience about the GUI usage, a test section involving dummy example using SummarizedExperiment objects of transcriptome (small form) and clinical datasets is also included.




Figure-1 SVMDO GUI Sections

Figure 1: Figure-1 SVMDO GUI Sections

6 How to open SVMDO Main Screen

To open GUI screen, you can directly write SVMDO::runGUI() in R console. If library is previously activated you can open GUI screen by writing runGUI() in R console.

7 How to use SVMDO Main Screen

(Important: Except file input processes, each step provides a message box indicating process success/failure. It disappears after clicking on any area in the GUI screen. This is necessary for continuing the steps.)

7.1 Steps of Analysis Screen

  1. To search your transcriptome dataset, use the file detection in Choose Your Expression Dataset section. The file will be automatically uploaded into the GUI.
  2. To prevent clashing with test datasets, select “None” option from the radio button section.
  3. By clicking on DEG Analysis button you further apply differential expression analysis. Labels of tissue_type column in dataset must contain “Nor” and “Tum” for determining normal/tumour (or tumor) samples.
  4. When the differential expression process is completed, a user-defined input size (n) is selected to filter the initial gene list (i.e., n number of upregulated and downregulated genes) by entering a number in Input Size section. It is predetermined as 50 in GUI which can be changed based on the user. A message window saying process completed will appear if there is not any problem. If there is problem with the value of input size, you will get a warning about inappropriate input size selection. If the input size remains, algorithm selects all of the differentiallly expressed genes to be used in the next process.
  5. To apply disease ontology-based gene filtration, click on DO Analysis button.
  6. To further apply the following feature selection and classification processes, click on the Classification button.
  7. Acquired discriminative gene set can be further used for survival analysis to detect individual prognostic genes. To apply this process, use the file detection in Choose Clinical Data section for searching clinical data about patient survival followed by clicking on Survival Analysis button.

7.2 Steps of Result Screen

  1. To visualize discriminative gene sets inside GUI screen, click on Show Gene Results button. When you click this button, a table of gene set will appear.
  2. To visualize survival plots of individual genes, two steps have to be applied. First of all, click on Prepare Plot Lists button to feed plot information to the visualization system. After that, click on Show Plots button to visualize survival plots.
  3. Before downloading files, you can adjust the output directory with Choose Directory button. It can be used for separating files by selecting a destination before clicking download buttons. If it is desired, files can be downloaded to the same folder by selecting an output directory just one time before the download steps. If you do not select any output directory, files will be downloaded to your working directory.
  4. To download the resulting discriminative gene set, it is obligatory to define a filename in the Enter Final Gene Set Filename section. After that, you can click on Download Gene List button to complete the process.
  5. To download survival plots, you have to click on Download Plot List button. Names of plot files are automatically done by assigning gene names.

8 Application of test datasets

SVMDO includes test datasets providing dummy examples for gaining experience on the GUI usage. Test datasets consists of summarized experiment objects including expression and clinical datasets. These objects are saved in RDA forms and called from extdata folder of the package. As expression datasets, test files includes simplified forms of TCGA-COAD (COAD) and TCGA-LUSC (LUSC) with 400 genes In test-based analysis, predetermined expression and clinical datasets are automatically uploaded into the GUI. Furthermore, predefined input size (n=50) is also automatically applied. Therefore, users have to continue with DO Analysis after DEG Analysis without adjusting input size.

9 Workspace Clearance

When the user task is completed, click on the Clear Environment button to remove the global variables created during the algorithm sections. To prevent error in the next usages of GUI, it is a necessary process. It can be applied at any moment without the necessity of completing all of the steps of algorithm.

10 Output Files of SVMDO

11 References

Appendix

  1. Cai,J. et al. (2018) Feature selection in machine learning: A new perspective. Neurocomputing, 300, 70–79.
  2. Gallagher,E.J. and LeRoith,D. (2015) Obesity and diabetes: The increased risk of cancer and cancer-related mortality. Physiological Reviews, 95, 727–748.
  3. Kori,M. and Arga,K.Y. (2020) Pathways involved in viral oncogenesis: New perspectives from virus-host protein interactomics. Biochimica et Biophysica Acta (BBA) - Molecular Basis of Disease, 1866, 165885.
  4. LePendu,P. et al. (2011) Enabling enrichment analysis with the human disease ontology. Journal of Biomedical Informatics, 44, S31–S38.
  5. Liu,B. et al. (2004) BMC Bioinformatics, 5, 136.
  6. Ouardighi,A. el et al. (2007) Feature selection on supervised classification using wilks lambda statistic. In, 2007 international symposium on computational intelligence and intelligent informatics. IEEE.

A How to install R and RStudio GUI

RStudio requires R 3.3.0+. Choose a version of R that matches your computer’s operating system

A.1 Windows Operating System

**R-base Install:

  1. Copy/Paste the address https://posit.co/download/rstudio-desktop/
  2. Click on “Download and install R”
  3. In the new page, click on “Download R for Windows”
  4. Click on “install R for the first time”
  5. Click on “Download R-X.X.X-win.exe for Windows”to install R-base

**RStudio Install:

  1. Copy/Paste the address https://posit.co/download/rstudio-desktop/
  2. Click on either download RStudio desktop for windows
  3. Alternatively, .exe /.zip versions can be downloaded from All Installers section

A.2 Ubuntu Operating System

**R-base Install:

  1. Copy/Paste the address https://posit.co/download/rstudio-desktop/
  2. Click on “Download and install R”.
  3. Click on “Download R for Linux”
  4. In the parent directory, click on ubuntu/
  5. Run the following lines in terminal (as root or by prefixing sudo) to provide access to recent version of R:
  • sudo apt install –no-install-recommends software-properties-common dirmngr # add the signing key (by Michael Rutter) for these repos # To verify key, run gpg –show-keys /etc/apt/trusted.gpg.d/cran_ubuntu_key.asc # Fingerprint: E298A3A825C0D65DFD57CBB651716619E084DAB9
  • wget -qO- https://cloud.r-project.org/bin/linux/ubuntu/marutter_pubkey.asc | sudo tee -a /etc/apt/trusted.gpg.d/cran_ubuntu_key.asc # add the R 4.0 repo from CRAN – adjust ‘focal’ to ‘groovy’ or ‘bionic’ as needed
  • sudo add-apt-repository “deb https://cloud.r-project.org/bin/linux/ubuntu $(lsb_release -cs)-cran40/”
  1. After the previous lines, run sudo apt install –no-install-recommends r-base

**RStudio Install:

  1. Copy/Paste the address https://posit.co/download/rstudio-desktop/
  2. In All Installers section, select .deb or .debian.tar.gz version for Ubuntu
  3. Right click on the download and/or unzipped .deb file to install RStudio (Alternatively, .deb file can be installed from terminal by typing sudo dpkg -i installed_Rstudio_file.deb)

B Session Info

#> R version 4.3.0 RC (2023-04-13 r84269)
#> Platform: x86_64-pc-linux-gnu (64-bit)
#> Running under: Ubuntu 22.04.2 LTS
#> 
#> Matrix products: default
#> BLAS:   /home/biocbuild/bbs-3.17-bioc/R/lib/libRblas.so 
#> LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.10.0
#> 
#> locale:
#>  [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
#>  [3] LC_TIME=en_GB              LC_COLLATE=C              
#>  [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
#>  [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
#>  [9] LC_ADDRESS=C               LC_TELEPHONE=C            
#> [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       
#> 
#> time zone: America/New_York
#> tzcode source: system (glibc)
#> 
#> attached base packages:
#> [1] stats     graphics  grDevices utils     datasets  methods   base     
#> 
#> other attached packages:
#> [1] BiocStyle_2.28.0
#> 
#> loaded via a namespace (and not attached):
#>  [1] digest_0.6.31       R6_2.5.1            codetools_0.2-19   
#>  [4] bookdown_0.33       fastmap_1.1.1       xfun_0.39          
#>  [7] cachem_1.0.7        knitr_1.42          htmltools_0.5.5    
#> [10] rmarkdown_2.21      cli_3.6.1           sass_0.4.5         
#> [13] jquerylib_0.1.4     compiler_4.3.0      highr_0.10         
#> [16] tools_4.3.0         evaluate_0.20       bslib_0.4.2        
#> [19] yaml_2.3.7          BiocManager_1.30.20 jsonlite_1.8.4     
#> [22] rlang_1.1.0