1. Gene Expression Omnibus (GEO)
NOTE: TMOD3 expression in ovarian cancer, in ovarian cancer treated with platinum drugs, and in drug-resistant ovarian cancer were derived from the GEO datasets. The study type of all datasets was expression profiling by array, and the organisms were Homo sapiens.
- Go to the GEO database (see Table of Materials), and input keywords such as TMOD3, ovarian cancer, and drug-resistant or data accession in the search box (see Supplementary Figure 1A).
- Divide data into different groups in the Define group box. In the option menu, choose Benjamini and Hochberg (False discovery rate) in the Apply adjustment to the P-values box, and then click on analyzed in the GEO2R menu. Click on download full table to get the results.
- Investigate and plot the downloaded data with graphing and statistics software (see Supplementary Figure 1B).
- Use unpaired t-test for the comparison between two groups.
2. TNMplot
NOTE: TNMplot utilizes RNA-seq data from The Cancer Genome Atlas (TCGA), Therapeutic Research to Generate Effective Treatments (TARGET), and Genotype-Tissue Expression (GTEx) repositories20. The expression of TMOD3 in normal ovarian tissues and ovarian cancer was analyzed using TNMplot.
- Go to the TNMplot web tool (see Table of Materials), and click on Compare Tumor and Normal.
- Insert TMOD3 in Choose a gene box, and choose Ovarian Serous Cystadenocarcinoma in Choose tissue box.
- Click on start analysis to get the results (see Supplementary Figure 1C).
3. UALCAN
NOTE: The University of ALabama at Birmingham CANcer data analysis Portal (UALCAN) is a user-friendly online resource for analyzing publicly available cancer data21. Protein level expression of TMOD3 in normal tissue and ovarian cancer in CPTAC data was analyzed using UALCAN.
- Go to UALCAN (see Table of Materials) and click on the Proteomics menu.
- Insert TMOD3 in the Enter gene names box, and choose ovarian cancer in the CPTAC dataset box.
- Click on Explore to get the results (see Supplementary Figure 1D).
4. KM-plotter prognostic analysis
NOTE: The prognostic value of TMOD3 in ovarian cancer was analyzed using Kaplan Meier plotter (KM-plotter), including overall survival (OS) and progression-free survival (PFS)22.
- Go to KM-plotter (see Table of Materials), and click on Start KM-plotter for ovarian cancer.
- Insert TMOD3 in Affy id/ Gene symbol box.
- Choose Auto select best cutoff in Split patients by box.
- Choose contains platin in the chemotherapy option when performing prognosis analysis for platinum-based chemotherapy patients.
- Choose TCGA in Use following dataset(s) for the analysis box.
- Click on Draw Kaplan-Meier plot to get the results (see Supplementary Figure 2).
5. ROC plotter
NOTE: The Receiver Operating Characteristic Curve (ROC) Plotter was used to analyze the expression of TMOD3 in patients resistant or sensitive to platinum-based chemotherapy and allows validation of the interested gene as a predictive marker by ROC curves. Datasets of ROC plotter at the transcriptome level are mainly from the TCGA and GEO databases and contain treatment and response data from 1816 ovarian cancer patients23.
- Go to The ROC Plotter (see Table of Materials), and click on ROC plotter for ovarian cancer.
- Insert TMOD3 in the Gene symbol box.
- Choose Relapse-free survival at 6 months in the Responds box, and choose Platin in the Treatment box.
- Click on Calculate to get the results (see Supplementary Figure 3).
6. mRNA-miRNA analysis
NOTE: The miRNAs targeting TMOD3 were predicted by TargetScan24, and then the correlation of TMOD3 with these miRNAs in the TCGA ovarian cancer dataset was analyzed by cBioportal25. Then, the result above was visualized by Cytoscape26. MiRNA expression in cisplatin-sensitive and drug-resistant ovarian cancer patients was analyzed by LinkedOmics27.
- Go to TargetScan (see Table of Materials), and insert TMOD3 in the Enter a human gene symbol box (see Supplementary Figure 4A).
- Go to cBioportal (see Table of Materials).
- Choose Ovarian Serous Cystadenocarcinoma (TCGA, Nature 2011) dataset and insert TMOD3 and miRNA symbols in the Enter Genes box.
- Click on Submit Query to get the correlation data of TMOD3 with miRNAs in the TCGA ovarian cancer data set (see Supplementary Figure 4B), then visualize the result by Cytoscape (see Table of Materials) (see Supplementary Figure 4C).
- Go to LinkedOmics (see Table of Materials), and select TCGA_OV Sample cohort.
- Choose clinical in the Select Search Dataset box, and select platinum status in the Select Search Dataset Attribute box.
- Choose miRNASeq in the Select Target Dataset box, and choose t-test in Select Statistical Method.
- Click on Submit Query to get the results (see Supplementary Figure 4D).
7. Immuno-infiltration analysis
NOTE: The Human Protein Atlas (HPA) database was used to analyze the distribution of TMOD3 in various immune cells. TIMER is a convenient online database that analyzes immune infiltration associated with multiple cancer types28. This study used TIMER to analyze the relationship between TMOD3 mRNA expression and ovarian cancer purity, and immune cell infiltration. TISIDB is an online portal for tumor-immune system interactions29. This study used TISIDB to determine the correlation between TMOD3 and immunomodulators in ovarian cancer.
- Go to HPA (see Table of Materials), and Insert TMOD3 in the Search box, then get the result.
- Choose Immune to show the distribution of TMOD3 in various immune cells (see Supplementary Figure 5A).
- Go to TIMER (see Table of Materials) and insert TMOD3 in the Gene Symbol box.
- Choose OV in the cancer types box and click on Submit to get the result (see Supplementary Figure 5B).
- Go to TISIDB (see Table of Materials) and insert TMOD3 in the Gene Symbol box.
- Click on Submit to get the result (see Supplementary Figure 5C).
8. TMOD3 coexpression networks in ovarian cancer
NOTE: Genes co-expressed with TMOD3 were analyzed by LinkedOmics, and TOP50 genes were displayed by heat maps. TMOD3 interacting genes were predicted by STRING30. Then, the overlapping genes were displayed by the Venn diagram. The overlapping genes were functionally annotated by DAVID31for Gene Ontology Biological Process (GO-BP), Gene Ontology Cellular Component (GO-CC), Gene Ontology Molecular Function (GO-MF), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways.
- Go to LinkedOmics, and select TCGA_OV Sample cohort.
- Choose RNAseq in Select Search Dataset box and Select Target Dataset box.
- Click on Submit Query to get the results of genes co-expressed with TMOD3.
- Go to STRING (see Table of Materials), and insert TMOD3 in the Protein Name box.
- Choose homo sapiens in the Organisms box and click on search to get the result (see Supplementary Figure 6).
- Go to DAVID (see Table of Materials).
- Insert the overlapping genes in the Enter Gene List box and choose ENSEMBL _GENE_SYMBOL in the Select Identifier box.
- Click on Submit List to get the results (see Supplementary Figure 7).
9. CTD database
NOTE: The CTD database is a novel tool to analyze the relationships between chemistry, genes, phenotypes, disease, and the environment32. The CTD database predicts drugs that target TMOD3. The PubChem database is then used to determine the definitive molecular structure of the drug.
- Go to the CTD database (see Table of Materials) and choose Chemical - Gene Interaction Query in the Search menu.
- Insert TMOD3 in the GENE box, and click on search to get the results (see Supplementary Figure 8).
- Go to the PubChem database (see Table of Materials) and insert drugs in the search box to get the results (see Supplementary Figure 9).