Public data and open source tools for multi-assay genomic investigation of disease

The mechanism underlying the oncogenic potential of AAA+ ATPase PSMC4 in cancer is revealed by mutations and copy number amplifications

Recent research has discovered a connection between the AAA+ ATPase PSMC4 (Proteasome 26S Subunit, ATPase 4) and several forms of cancer. However, a detailed analysis of the oncogenic potential of PSMC4 was elusive. In this study, we anticipate PSMC4's potential as a cancer biomarker. We aimed to comprehensively assess the expression profiles, prognostic significance, and relevant cellular pathways associated with it. Through our examination of various types of cancers, PSMC4 is found to be overexpressed. Interestingly, our result finds a positive correlation between PSMC4 overexpression and unfavourable overall survival rates in cancer. Further, we looked into the mutations and copy number amplifications of PSMC4 across various cancers. Our study reveals that missense mutations plays a great role behind the oncogenic potential of PSMC4. Several possible mutation sites are predicted. Interestingly, we found fifteen hotspot mutations in the ATPase domain of PSMC4. Additionally, PSMC4 has shown a high amplification percentage in various cancers. We are additionally attentive to the functional characteristics of the protein PSMC4 across various types of cancer. In the protein-protein interaction analyses, it was found that multiple oncoproteins were directly interacting with PSMC4. The top signaling pathways of PSMC4 also indicate that it plays a crucial role in cancer development. Overall, this study reveals that PSMC4 could be a potential diagnostic and prognostic marker for cancer, making it a promising biomarker and target.

Large-Scale Databases and Portals on Cancer Genome to Analyze Chaperone Genes Correlated to Patient Prognosis

Molecular chaperones, such as heat shock proteins (HSPs), have attracted attention as molecules involved in malignant events in cancers and are potential therapeutic targets and biomarkers for tumor therapy. Furthermore, mutations in chaperones can significantly impact cancer risk and prognosis. Bioinformatics is a particularly useful method for developing biomarkers as a practical consideration for the immediate clinical application of data. Many large-scale databases and portals on cancer genome are nowadays publicly available, including the International Cancer Genome Consortium (ICGC); The Cancer Genome Atlas (TCGA), renamed as Genomic Data Commons (GDC); Catalogue of Somatic Mutations in Cancer (COSMIC); and Cancer Cell Line Encyclopedia (CCLE). Referring to these databases, advanced web portals are publicized, including cBioPortal, Human Protein Atlas (HPA), Kaplan-Meier (KM) plotter, Gene Expression Profiling Interactive Analysis 2 (GEPIA2), Genomics of Drug Sensitivity in Cancer (GDSC), and Dependency Map (DepMap). Here, we assemble these databases and portals to clarify what is available and useful for current cancer research and provide protocols to utilize the HPA, KM plotter, and GEPIA2 for studies on chaperone genes in cancer patients. Utilizing these portals will reveal the correlation between tumor subtype-specific high expression of chaperone genes and patient prognosis. Our protocols are useful to increase systematic awareness of chaperones and find new biomarkers for diagnosis and prognosis and new targets for anticancer drugs.

Overview of Bioinformatics Software and Databases for Metabolic Engineering

The explosion of the "omics" era has introduced a growing number of sets and tools that facilitate molecular interrogation of the metabolome. These include various bioinformatics and pharmacogenomics resources that can be utilized independently or collectively to facilitate metabolic engineering across disease, clinical oncology, and understanding of molecular changes across larger systems. This review provides starting points for accessing publicly available data and computational tools that support assessment of metabolic profiles and metabolic regulation, providing both a depth-and-breadth approach toward understanding the metabolome. We focus in particular on pathway databases and tools, which provide in-depth analysis of metabolic pathways, which is at the heart of metabolic engineering.

STATegra: Multi-Omics Data Integration - A Conceptual Scheme With a Bioinformatics Pipeline

Technologies for profiling samples using different omics platforms have been at the forefront since the human genome project. Large-scale multi-omics data hold the promise of deciphering different regulatory layers. Yet, while there is a myriad of bioinformatics tools, each multi-omics analysis appears to start from scratch with an arbitrary decision over which tools to use and how to combine them. Therefore, it is an unmet need to conceptualize how to integrate such data and implement and validate pipelines in different cases. We have designed a conceptual framework (STATegra), aiming it to be as generic as possible for multi-omics analysis, combining available multi-omic anlaysis tools (machine learning component analysis, non-parametric data combination, and a multi-omics exploratory analysis) in a step-wise manner. While in several studies, we have previously combined those integrative tools, here, we provide a systematic description of the STATegra framework and its validation using two The Cancer Genome Atlas (TCGA) case studies. For both, the Glioblastoma and the Skin Cutaneous Melanoma (SKCM) cases, we demonstrate an enhanced capacity of the framework (and beyond the individual tools) to identify features and pathways compared to single-omics analysis. Such an integrative multi-omics analysis framework for identifying features and components facilitates the discovery of new biology. Finally, we provide several options for applying the STATegra framework when parametric assumptions are fulfilled and for the case when not all the samples are profiled for all omics. The STATegra framework is built using several tools, which are being integrated step-by-step as OpenSource in the STATegRa Bioconductor package.

A semi-automatic methodology for analysing distributed and private biobanks

Privacy issues limit the analysis and cross-exploration of most distributed and private biobanks, often raised by the multiple dimensionality and sensitivity of the data associated with access restrictions and policies. These characteristics prevent collaboration between entities, constituting a barrier to emergent personalized and public health challenges, namely the discovery of new druggable targets, identification of disease-causing genetic variants, or the study of rare diseases. In this paper, we propose a semi-automatic methodology for the analysis of distributed and private biobanks. The strategies involved in the proposed methodology efficiently enable the creation and execution of unified genomic studies using distributed repositories, without compromising the information present in the datasets. We apply the methodology to a case study in the current Covid-19, ensuring the combination of the diagnostics from multiple entities while maintaining privacy through a completely identical procedure. Moreover, we show that the methodology follows a simple, intuitive, and practical scheme.

Omics Integration Analysis Unravel the Landscape of Driving Mechanisms of Colorectal Cancer

Colorectal cancer (CRC) is one of the most malignant cancers and results in a substantial rate of morbidity and mortality. Diagnosis of this malignancy in early stages increases the chance of effective treatment. High-throughput data analyses reveal omics signatures and also provide the possibility of developing computational models for early detection of this disease. Such models would be able to use as complementary tools for early detection of different types of cancers including CRC. In this study, using gene expression data, the Flux balance analysis (FBA) applied to decode metabolic fluxes in cancer and normal cells. Moreover, transcriptome and genome analyses revealed driver agents of CRC in a biological network scheme. By applying comprehensive publicly available data from TCGA, different aspect of CRC regulome including the regulatory effect of gene expression, methylation, microRNA, copy number aberration and point mutation profile over protein levels investigated and the results provide a regulatory picture underlying CRC. Compiling omics profiles indicated snapshots of changes in different omics levels and flux rate of CRC. In conclusion, considering obtained CRC signatures and their role in biological operating systems of cells, the results suggest reliable driver regulatory modules that could potentially serve as biomarkers and therapeutic targets and furthermore expand our understanding of driving mechanisms of this disease.
.

Multiomic Integration of Public Oncology Databases in Bioconductor

PURPOSE

Investigations of the molecular basis for the development, progression, and treatment of cancer increasingly use complementary genomic assays to gather multiomic data, but management and analysis of such data remain complex. The cBioPortal for cancer genomics currently provides multiomic data from > 260 public studies, including The Cancer Genome Atlas (TCGA) data sets, but integration of different data types remains challenging and error prone for computational methods and tools using these resources. Recent advances in data infrastructure within the Bioconductor project enable a novel and powerful approach to creating fully integrated representations of these multiomic, pan-cancer databases.

METHODS

We provide a set of R/Bioconductor packages for working with TCGA legacy data and cBioPortal data, with special considerations for loading time; efficient representations in and out of memory; analysis platform; and an integrative framework, such as MultiAssayExperiment. Large methylation data sets are provided through out-of-memory data representation to provide responsive loading times and analysis capabilities on machines with limited memory.

RESULTS

We developed the curatedTCGAData and cBioPortalData R/Bioconductor packages to provide integrated multiomic data sets from the TCGA legacy database and the cBioPortal web application programming interface using the MultiAssayExperiment data structure. This suite of tools provides coordination of diverse experimental assays with clinicopathological data with minimal data management burden, as demonstrated through several greatly simplified multiomic and pan-cancer analyses.

CONCLUSION

These integrated representations enable analysts and tool developers to apply general statistical and plotting methods to extensive multiomic data through user-friendly commands and documented examples.

Paired like homeodomain 1 and SAM and SH3 domain-containing 1 in the progression and prognosis of head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) is an aggressive malignancy with high morbidity and mortality rates. In spite of numerous advancements have been made in therapeutic methods, the prognosis of HNSCC patients remains poor. Therefore, investigation of crucial genes during HNSCC tumorigenesis which could be exploited as biomarkers and therapeutic targets is greatly needed. In this study, original data of four independent datasets was downloaded from the Gene Expression Omnibus database and analyzed through R language to screen out differentially expressed genes. Paired like homeodomain 1 and SAM and SH3 domain-containing 1 were selected to be further explored through multiple online databases. Quantitative real-time polymerase chain reaction analysis and immunohistochemistry assay were adopted to validate the downregulation of paired like homeodomain 1 and SAM and SH3 domain-containing 1 in HNSCC and statistical analysis indicated their close associations with patient prognosis. In vitro experiments demonstrated the inhibitory effect of paired like homeodomain 1 and SAM and SH3 domain-containing 1 on HNSCC progression. Overall, we identified the aberrant downregulation of paired like homeodomain 1 and SAM and SH3 domain-containing 1 in HNSCC and suggested the potential of utilizing them as therapeutic targets or efficient biomarkers for diagnosis and prognosis evaluation. Our findings may provide novel evidences for the development of new strategies for HNSCC treatment.

RFCM3: Computational Method for Identification of miRNA-mRNA Regulatory Modules in Cervical Cancer

Cervical cancer is a leading severe malignancy throughout the world. Molecular processes and biomarkers leading to tumor progression in cervical cancer are either unknown or only partially understood. An increasing number of studies have shown that microRNAs play an important role in tumorigenesis so understanding the regulatory mechanism of miRNAs in gene-regulatory network will help elucidate the complex biological processes that occur during malignancy. Functional genomics data provides opportunities to study the aberrant microRNA-messenger RNA (miRNA-mRNA) interaction. Identification of miRNA-mRNA regulatory modules will aid deciphering aberrant transcriptional regulatory network in cervical cancer but is computationally challenging. In this regard, an algorithm, termed as relevant and functionally consistent miRNA-mRNA modules (RFCM³), is proposed. It integrates miRNA and mRNA expression data of cervical cancer for identification of potential miRNA-mRNA modules. It selects set of miRNA-mRNA modules by maximizing relation of mRNAs with miRNA and functional similarity between selected mRNAs. Later, using the knowledge of the miRNA-miRNA synergistic network different modules are fused and finally a set of modules are generated containing several miRNAs as well as mRNAs. This type of module explains the underlying biological pathways containing multiple miRNAs and mRNAs. The effectiveness of the proposed approach over other existing methods has been demonstrated on a miRNA and mRNA expression data of cervical cancer with respect to enrichment analyses and other standard metrices. The prognostic value of the genes in a module with respect to cervical cancer is also demonstrated. The approach was found to generate more robust, integrated, and functionally enriched miRNA-mRNA modules in cervical cancer.

Advancing from protein interactomes and gene co-expression networks towards multi-omics-based composite networks: approaches for predicting and extracting biological knowledge

Prediction of biological interaction networks from single-omics data has been extensively implemented to understand various aspects of biological systems. However, more recently, there is a growing interest in integrating multi-omics datasets for the prediction of interactomes that provide a global view of biological systems with higher descriptive capability, as compared to single omics. In this review, we have discussed various computational approaches implemented to infer and analyze two of the most important and well studied interactomes: protein-protein interaction networks and gene co-expression networks. We have explicitly focused on recent methods and pipelines implemented to infer and extract biologically important information from these interactomes, starting from utilizing single-omics data and then progressing towards multi-omics data. Accordingly, recent examples and case studies are also briefly discussed. Overall, this review will provide a proper understanding of the latest developments in protein and gene network modelling and will also help in extracting practical knowledge from them.

Vertical and horizontal integration of multi-omics data with miodin

BACKGROUND

Studies on multiple modalities of omics data such as transcriptomics, genomics and proteomics are growing in popularity, since they allow us to investigate complex mechanisms across molecular layers. It is widely recognized that integrative omics analysis holds the promise to unlock novel and actionable biological insights into health and disease. Integration of multi-omics data remains challenging, however, and requires combination of several software tools and extensive technical expertise to account for the properties of heterogeneous data.

RESULTS

This paper presents the miodin R package, which provides a streamlined workflow-based syntax for multi-omics data analysis. The package allows users to perform analysis of omics data either across experiments on the same samples (vertical integration), or across studies on the same variables (horizontal integration). Workflows have been designed to promote transparent data analysis and reduce the technical expertise required to perform low-level data import and processing.

CONCLUSIONS

The miodin package is implemented in R and is freely available for use and extension under the GPL-3 license. Package source, reference documentation and user manual are available at https://gitlab.com/algoromics/miodin.

Computational resources associating diseases with genotypes, phenotypes and exposures

The causes of a disease and its therapies are not only related to genotypes, but also associated with other factors, including phenotypes, environmental exposures, drugs and chemical molecules. Distinguishing disease-related factors from many neutral factors is critical as well as difficult. Over the past two decades, bioinformaticians have developed many computational resources to integrate the omics data and discover associations among these factors. However, researchers and clinicians are experiencing difficulties in choosing appropriate resources from hundreds of relevant databases and software tools. Here, in order to assist the researchers and clinicians, we systematically review the public computational resources of human diseases related to genotypes, phenotypes, environment factors, drugs and chemical exposures. We briefly describe the development history of these computational resources, followed by the details of the relevant databases and software tools. We finally conclude with a discussion of current challenges and future opportunities as well as prospects on this topic.

Bringing radiomics into a multi-omics framework for a comprehensive genotype-phenotype characterization of oncological diseases

Genomic and radiomic data integration, namely radiogenomics, can provide meaningful knowledge in cancer diagnosis, prognosis and treatment. Despite several data structures based on multi-layer architecture proposed to combine multi-omic biological information, none of these has been designed and assessed to include radiomic data as well. To meet this need, we propose to use the MultiAssayExperiment (MAE), an R package that provides data structures and methods for manipulating and integrating multi-assay experiments, as a suitable tool to manage radiogenomic experiment data. To this aim, we first examine the role of radiogenomics in cancer phenotype definition, then the current state of radiogenomics data integration in public repository and, finally, challenges and limitations of including radiomics in MAE, designing an extended framework and showing its application on a case study from the TCGA-TCIA archives. Radiomic and genomic data from 91 patients have been successfully integrated in a single MAE object, demonstrating the suitability of the MAE data structure as container of radiogenomic data.

The metaRbolomics Toolbox in Bioconductor and beyond

Metabolomics aims to measure and characterise the complex composition of metabolites in a biological system. Metabolomics studies involve sophisticated analytical techniques such as mass spectrometry and nuclear magnetic resonance spectroscopy, and generate large amounts of high-dimensional and complex experimental data. Open source processing and analysis tools are of major interest in light of innovative, open and reproducible science. The scientific community has developed a wide range of open source software, providing freely available advanced processing and analysis approaches. The programming and statistics environment R has emerged as one of the most popular environments to process and analyse Metabolomics datasets. A major benefit of such an environment is the possibility of connecting different tools into more complex workflows. Combining reusable data processing R scripts with the experimental data thus allows for open, reproducible research. This review provides an extensive overview of existing packages in R for different steps in a typical computational metabolomics workflow, including data processing, biostatistics, metabolite annotation and identification, and biochemical network and pathway analysis. Multifunctional workflows, possible user interfaces and integration into workflow management systems are also reviewed. In total, this review summarises more than two hundred metabolomics specific packages primarily available on CRAN, Bioconductor and GitHub.

In silico drug repositioning: from large-scale transcriptome data to therapeutics

Drug repositioning is an attractive alternative to conventional drug development when new beneficial effects of old drugs are clinically validated because pharmacokinetic and safety profiles are generally already available. Since ~ 30% of drugs newly approved by the US food and drug administration (FDA) are developed through drug repositioning, identifying novel usage for existing drugs is an emerging strategy for developing disease treatments. With advances in next-generation sequencing technologies, available transcriptome data related to diseases have expanded rapidly. Harnessing these resources enables a better understanding of disease mechanisms and drug mode of action (MOA), and moves toward personalized pharmacotherapy. In this review, we briefly outline publicly available large-scale transcriptome databases and tools for drug repositioning. We also highlight recent approaches leading to the discovery of novel drug targets, drug response biomarkers, drug indications, and drug MOA.

Bioinformatics-based discovery of PYGM and TNNC2 as potential biomarkers of head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) is an aggressive malignancy with high morbidity and mortality rates and ranks as the sixth most common cancer all over the world. Despite numerous advancements in therapeutic methods, the prognosis of HNSCC patients still remains poor. Therefore, there is an urgent need to have a better understanding of the molecular mechanisms underlying HNSCC progression and to identify essential genes that could serve as effective biomarkers and potential treatment targets. In the present study, original data of three independent datasets were downloaded from the Gene Expression Omnibus database (GEO) and R language was applied to screen out the differentially expressed genes (DEGs). PYGM and TNNC2 were finally selected from the overlapping DEGs of three datasets for further analyses. Transcriptional and survival data related to PYGM and TNNC2 was detected through multiple online databases such as Oncomine, Gene Expression Profiling Interactive Analysis (GEPIA), cBioportal, and UALCAN. Quantitative real-time polymerase chain reaction (qPCR) analysis was adopted for the validation of PYGM and TNNC2 mRNA level in HNSCC tissues and cell lines. Survival curves were plotted to evaluate the association of these two genes with HNSCC prognosis. It was demonstrated that PYGM and TNNC2 were significantly down-regulated in HNSCC and the aberrant expression of PYGM and TNNC2 were correlated with HNSCC prognosis, implying the potential of exploiting them as therapeutic targets for HNSCC treatment or potential biomarkers for diagnosis and prognosis.

MetaGxData: Clinically Annotated Breast, Ovarian and Pancreatic Cancer Datasets and their Use in Generating a Multi-Cancer Gene Signature

A wealth of transcriptomic and clinical data on solid tumours are under-utilized due to unharmonized data storage and format. We have developed the MetaGxData package compendium, which includes manually-curated and standardized clinical, pathological, survival, and treatment metadata across breast, ovarian, and pancreatic cancer data. MetaGxData is the largest compendium of curated transcriptomic data for these cancer types to date, spanning 86 datasets and encompassing 15,249 samples. Open access to standardized metadata across cancer types promotes use of their transcriptomic and clinical data in a variety of cross-tumour analyses, including identification of common biomarkers, and assessing the validity of prognostic signatures. Here, we demonstrate that MetaGxData is a flexible framework that facilitates meta-analyses by using it to identify common prognostic genes in ovarian and breast cancer. Furthermore, we use the data compendium to create the first gene signature that is prognostic in a meta-analysis across 3 cancer types. These findings demonstrate the potential of MetaGxData to serve as an important resource in oncology research, and provide a foundation for future development of cancer-specific compendia.

Metabolic gene NR4A1 as a potential therapeutic target for non-smoking female non-small cell lung cancer patients

BACKGROUND

Although cigarette smoking is considered one of the key risk factors for lung cancer, 15% of male patients and 53% of female patients with lung cancer are non-smokers. Metabolic changes are critical features of cancer. Therapeutic target identification from a metabolic perspective in non-small cell lung cancer (NSCLC) tissue of female non-smokers has long been ignored.

RESULTS

Based on microarray data retrieved from Affymetrix expression arrays E-GEOD-19804, we found that the downregulated genes in non-smoking female NSCLC patients tended to participate in protein/amino acid and lipid metabolism, while upregulated genes were more involved in protein/amino acid and carbohydrate metabolism. Combining nutrient metabolic co-expression, protein-protein interaction network construction and overall survival assessment, we identified NR4A1 and TIE1 as potential therapeutic targets for NSCLC in female non-smokers. To accelerate the drug development for non-smoking female NSCLC patients, we identified nilotinib as a potential agonist targeting NR4A1 encoded protein by molecular docking and molecular dynamic stimulation. We also show that nilotinib inhibited proliferation and induced senescence of cells in non-smoking female NSCLC patients in vitro.

CONCLUSIONS

These results not only uncover nutrient metabolic characteristics in non-smoking female NSCLC patients, but also provide a new paradigm for identifying new targets and drugs for novel therapy for such patients.

Epigenetic regulation of gene expression in cancer: techniques, resources and analysis

Cancer is a complex disease, driven by aberrant activity in numerous signaling pathways in even individual malignant cells. Epigenetic changes are critical mediators of these functional changes that drive and maintain the malignant phenotype. Changes in DNA methylation, histone acetylation and methylation, noncoding RNAs, posttranslational modifications are all epigenetic drivers in cancer, independent of changes in the DNA sequence. These epigenetic alterations were once thought to be crucial only for the malignant phenotype maintenance. Now, epigenetic alterations are also recognized as critical for disrupting essential pathways that protect the cells from uncontrolled growth, longer survival and establishment in distant sites from the original tissue. In this review, we focus on DNA methylation and chromatin structure in cancer. The precise functional role of these alterations is an area of active research using emerging high-throughput approaches and bioinformatics analysis tools. Therefore, this review also describes these high-throughput measurement technologies, public domain databases for high-throughput epigenetic data in tumors and model systems and bioinformatics algorithms for their analysis. Advances in bioinformatics data that combine these epigenetic data with genomics data are essential to infer the function of specific epigenetic alterations in cancer. These integrative algorithms are also a focus of this review. Future studies using these emerging technologies will elucidate how alterations in the cancer epigenome cooperate with genetic aberrations during tumor initiation and progression. This deeper understanding is essential to future studies with epigenetics biomarkers and precision medicine using emerging epigenetic therapies.

Genome, transcriptome and proteome: the rise of omics data and their integration in biomedical sciences

Advances in the technologies and informatics used to generate and process large biological data sets (omics data) are promoting a critical shift in the study of biomedical sciences. While genomics, transcriptomics and proteinomics, coupled with bioinformatics and biostatistics, are gaining momentum, they are still, for the most part, assessed individually with distinct approaches generating monothematic rather than integrated knowledge. As other areas of biomedical sciences, including metabolomics, epigenomics and pharmacogenomics, are moving towards the omics scale, we are witnessing the rise of inter-disciplinary data integration strategies to support a better understanding of biological systems and eventually the development of successful precision medicine. This review cuts across the boundaries between genomics, transcriptomics and proteomics, summarizing how omics data are generated, analysed and shared, and provides an overview of the current strengths and weaknesses of this global approach. This work intends to target students and researchers seeking knowledge outside of their field of expertise and fosters a leap from the reductionist to the global-integrative analytical approach in research.

Large-scale analysis reveals a novel risk score to predict overall survival in hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) is a major cause of cancer mortality and an increasing incidence worldwide; however, there are very few effective diagnostic approaches and prognostic biomarkers.

Materials and methods

One hundred forty-nine pairs of HCC samples from Gene Expression Omnibus (GEO) were obtained to screen differentially expressed genes (DEGs) between HCC and normal samples. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, Gene ontology enrichment analyses, and protein–protein interaction network were used. Cox proportional hazards regression analysis was used to identify significant prognostic DEGs, with which a gene expression signature prognostic prediction model was identified in The Cancer Genome Atlas (TCGA) project discovery cohort. The robustness of this panel was assessed in the GSE14520 cohort. We verified details of the gene expression level of the key molecules through TCGA, GEO, and qPCR and used immunohistochemistry for substantiation in HCC tissues. The methylation states of these genes were also explored.

Results

Ninety-eight genes, consisting of 13 upregulated and 85 downregulated genes, were screened out in three datasets. KEGG and Gene ontology analysis for the DEGs revealed important biological features of each subtype. Protein–protein interaction network analysis was constructed, consisting of 64 nodes and 115 edges. A subset of four genes (SPINK1, TXNRD1, LCAT, and PZP) that formed a prognostic gene expression signature was established from TCGA and validated in GSE14520. Next, the expression details of the four genes were validated with TCGA, GEO, and clinical samples. The expression panels of the four genes were closely related to methylation states.

Conclusion

This study identified a novel four-gene signature biomarker for predicting the prognosis of HCC. The biomarkers may also reveal molecular mechanisms underlying development of the disease and provide new insights into interventional strategies.

Systems Medicine Disease: Disease Classification and Scalability Beyond Networks and Boundary Conditions

In order to accommodate the forthcoming wealth of health and disease related information, from genome to body sensors to population and the environment, the approach to disease description and definition demands re-examination. Traditional classification methods remain trapped by history; to provide the descriptive features that are required for a comprehensive description of disease, systems science, which realizes dynamic processes, adaptive response, and asynchronous communication channels, must be applied (Wolkenhauer et al., 2013). When Disease is viewed beyond the thresholds of lines and threshold boundaries, disease definition is not only the result of reductionist, mechanistic categories which reluctantly face re-composition. Disease is process and synergy as the characteristics of Systems Biology and Systems Medicine are included. To capture the wealth of information and contribute meaningfully to medical practice and biology research, Disease classification goes beyond a single spatial biologic level or static time assignment to include the interface of Disease process and organism response (Bechtel, 2017a; Green et al., 2017).

Low expression of KCNN3 may affect drug resistance in ovarian cancer

Drug resistance is a principal contributor to the poor prognosis of ovarian cancer (OC). Therefore, identifying factors that affect drug resistance in OC is critical. In the present study, 51 OC specimens from lab collections were immunohistochemically tested, public data for 489 samples from The Cancer Genome Atlas cohort and 1,656 samples from the Kaplan‑Meier Plotter were downloaded, and data were retrieved from Oncomine. It was identified that the mRNA and protein expression of the potassium calcium‑activated channel subfamily N member 3 (KCNN3) was markedly lower in OC tissues compared with normal tissues, and in drug‑resistant OC tissues compared with sensitive OC tissues. Low KCNN3 expression consistently predicted shorter disease‑free and overall survival (OS). Specifically, low KCNN3 expression predicted shorter OS in 395 patients with low expression levels of mucin‑16. There was additional evidence that KCNN3 expression is mediated by microRNA‑892b. Furthermore, text mining and analyses of protein and gene interactions indicated that KCNN3 affects drug resistance. To the best of the authors' knowledge, this is the first report to associate KCNN3 with poor prognosis and drug resistance in OC. The present findings indicated that KCNN3 is a potential prognostic marker and therapeutic target for OC.

A review of connectivity map and computational approaches in pharmacogenomics

Large-scale perturbation databases, such as Connectivity Map (CMap) or Library of Integrated Network-based Cellular Signatures (LINCS), provide enormous opportunities for computational pharmacogenomics and drug design. A reason for this is that in contrast to classical pharmacology focusing at one target at a time, the transcriptomics profiles provided by CMap and LINCS open the door for systems biology approaches on the pathway and network level. In this article, we provide a review of recent developments in computational pharmacogenomics with respect to CMap and LINCS and related applications.

A review of connectivity map and computational approaches in pharmacogenomics

Large-scale perturbation databases, such as Connectivity Map (CMap) or Library of Integrated Network-based Cellular Signatures (LINCS), provide enormous opportunities for computational pharmacogenomics and drug design. A reason for this is that in contrast to classical pharmacology focusing at one target at a time, the transcriptomics profiles provided by CMap and LINCS open the door for systems biology approaches on the pathway and network level. In this article, we provide a review of recent developments in computational pharmacogenomics with respect to CMap and LINCS and related applications.

Training replicable predictors in multiple studies

This article considers replicability of the performance of predictors across studies. We suggest a general approach to investigating this issue, based on ensembles of prediction models trained on different studies. We quantify how the common practice of training on a single study accounts in part for the observed challenges in replicability of prediction performance. We also investigate whether ensembles of predictors trained on multiple studies can be combined, using unique criteria, to design robust ensemble learners trained upfront to incorporate replicability into different contexts and populations.

Review of Statistical Learning Methods in Integrated Omics Studies (An Integrated Information Science)

Integrated omics is becoming a new channel for investigating the complex molecular system in modern biological science and sets a foundation for systematic learning for precision medicine. The statistical/machine learning methods that have emerged in the past decade for integrated omics are not only innovative but also multidisciplinary with integrated knowledge in biology, medicine, statistics, machine learning, and artificial intelligence. Here, we review the nontrivial classes of learning methods from the statistical aspects and streamline these learning methods within the statistical learning framework. The intriguing findings from the review are that the methods used are generalizable to other disciplines with complex systematic structure, and the integrated omics is part of an integrated information science which has collated and integrated different types of information for inferences and decision making. We review the statistical learning methods of exploratory and supervised learning from 42 publications. We also discuss the strengths and limitations of the extended principal component analysis, cluster analysis, network analysis, and regression methods. Statistical techniques such as penalization for sparsity induction when there are fewer observations than the number of features and using Bayesian approach when there are prior knowledge to be integrated are also included in the commentary. For the completeness of the review, a table of currently available software and packages from 23 publications for omics are summarized in the appendix.

Epigenetic regulation of gene expression in cancer: techniques, resources and analysis

Cancer is a complex disease, driven by aberrant activity in numerous signaling pathways in even individual malignant cells. Epigenetic changes are critical mediators of these functional changes that drive and maintain the malignant phenotype. Changes in DNA methylation, histone acetylation and methylation, noncoding RNAs, posttranslational modifications are all epigenetic drivers in cancer, independent of changes in the DNA sequence. These epigenetic alterations were once thought to be crucial only for the malignant phenotype maintenance. Now, epigenetic alterations are also recognized as critical for disrupting essential pathways that protect the cells from uncontrolled growth, longer survival and establishment in distant sites from the original tissue. In this review, we focus on DNA methylation and chromatin structure in cancer. The precise functional role of these alterations is an area of active research using emerging high-throughput approaches and bioinformatics analysis tools. Therefore, this review also describes these high-throughput measurement technologies, public domain databases for high-throughput epigenetic data in tumors and model systems and bioinformatics algorithms for their analysis. Advances in bioinformatics data that combine these epigenetic data with genomics data are essential to infer the function of specific epigenetic alterations in cancer. These integrative algorithms are also a focus of this review. Future studies using these emerging technologies will elucidate how alterations in the cancer epigenome cooperate with genetic aberrations during tumor initiation and progression. This deeper understanding is essential to future studies with epigenetics biomarkers and precision medicine using emerging epigenetic therapies.

Software for the Integration of Multiomics Experiments in Bioconductor

Multiomics experiments are increasingly commonplace in biomedical research and add layers of complexity to experimental design, data integration, and analysis. R and Bioconductor provide a generic framework for statistical analysis and visualization, as well as specialized data classes for a variety of high-throughput data types, but methods are lacking for integrative analysis of multiomics experiments. The MultiAssayExperiment software package, implemented in R and leveraging Bioconductor software and design principles, provides for the coordinated representation of, storage of, and operation on multiple diverse genomics data. We provide the unrestricted multiple 'omics data for each cancer tissue in The Cancer Genome Atlas as ready-to-analyze MultiAssayExperiment objects and demonstrate in these and other datasets how the software simplifies data representation, statistical analysis, and visualization. The MultiAssayExperiment Bioconductor package reduces major obstacles to efficient, scalable, and reproducible statistical analysis of multiomics data and enhances data science applications of multiple omics datasets. Cancer Res; 77(21); e39-42. ©2017 AACR.

Guidelines for Developing Successful Short Advanced Courses in Systems Medicine and Systems Biology

Systems medicine and systems biology have inherent educational challenges. These have largely been addressed either by providing new masters programs or by redesigning undergraduate programs. In contrast, short courses can respond to a different need: they can provide condensed updates for professionals across academia, the clinic, and industry. These courses have received less attention. Here, we share our experiences in developing and providing such courses to current and future leaders in systems biology and systems medicine. We present guidelines for how to reproduce our courses, and we offer suggestions for how to select students who will nurture an interdisciplinary learning environment and thrive there.

Dissecting the genomic activity of a transcriptional regulator by the integrative analysis of omics data

In the study of genomic regulation, strategies to integrate the data produced by Next Generation Sequencing (NGS)-based technologies in a meaningful ensemble are eagerly awaited and must continuously evolve. Here, we describe an integrative strategy for the analysis of data generated by chromatin immunoprecipitation followed by NGS which combines algorithms for data overlap, normalization and epigenetic state analysis. The performance of our strategy is illustrated by presenting the analysis of data relative to the transcriptional regulator Estrogen Receptor alpha (ERα) in MCF-7 breast cancer cells and of Glucocorticoid Receptor (GR) in A549 lung cancer cells. We went through the definition of reference cistromes for different experimental contexts, the integration of data relative to co-regulators and the overlay of chromatin states as defined by epigenetic marks in MCF-7 cells. With our strategy, we identified novel features of estrogen-independent ERα activity, including FoxM1 interaction, eRNAs transcription and a peculiar ontology of connected genes.

Synthesizer: Expediting synthesis studies from context-free data with information retrieval techniques

Scientists have unprecedented access to a wide variety of high-quality datasets. These datasets, which are often independently curated, commonly use unstructured spreadsheets to store their data. Standardized annotations are essential to perform synthesis studies across investigators, but are often not used in practice. Therefore, accurately combining records in spreadsheets from differing studies requires tedious and error-prone human curation. These efforts result in a significant time and cost barrier to synthesis research. We propose an information retrieval inspired algorithm, Synthesize, that merges unstructured data automatically based on both column labels and values. Application of the Synthesize algorithm to cancer and ecological datasets had high accuracy (on the order of 85-100%). We further implement Synthesize in an open source web application, Synthesizer (https://github.com/lisagandy/synthesizer). The software accepts input as spreadsheets in comma separated value (CSV) format, visualizes the merged data, and outputs the results as a new spreadsheet. Synthesizer includes an easy to use graphical user interface, which enables the user to finish combining data and obtain perfect accuracy. Future work will allow detection of units to automatically merge continuous data and application of the algorithm to other data formats, including databases.

MultiDataSet: an R package for encapsulating multiple data sets with application to omic data integration

BACKGROUND

Reduction in the cost of genomic assays has generated large amounts of biomedical-related data. As a result, current studies perform multiple experiments in the same subjects. While Bioconductor's methods and classes implemented in different packages manage individual experiments, there is not a standard class to properly manage different omic datasets from the same subjects. In addition, most R/Bioconductor packages that have been designed to integrate and visualize biological data often use basic data structures with no clear general methods, such as subsetting or selecting samples.

RESULTS

To cover this need, we have developed MultiDataSet, a new R class based on Bioconductor standards, designed to encapsulate multiple data sets. MultiDataSet deals with the usual difficulties of managing multiple and non-complete data sets while offering a simple and general way of subsetting features and selecting samples. We illustrate the use of MultiDataSet in three common situations: 1) performing integration analysis with third party packages; 2) creating new methods and functions for omic data integration; 3) encapsulating new unimplemented data from any biological experiment.

CONCLUSIONS

MultiDataSet is a suitable class for data integration under R and Bioconductor framework.

Discovering MicroRNA-Regulatory Modules in Multi-Dimensional Cancer Genomic Data: A Survey of Computational Methods

MicroRNAs (miRs) are small single-stranded noncoding RNA that function in RNA silencing and post-transcriptional regulation of gene expression. An increasing number of studies have shown that miRs play an important role in tumorigenesis, and understanding the regulatory mechanism of miRs in this gene regulatory network will help elucidate the complex biological processes at play during malignancy. Despite advances, determination of miR–target interactions (MTIs) and identification of functional modules composed of miRs and their specific targets remain a challenge. A large amount of data generated by high-throughput methods from various sources are available to investigate MTIs. The development of data-driven tools to harness these multi-dimensional data has resulted in significant progress over the past decade. In parallel, large-scale cancer genomic projects are allowing new insights into the commonalities and disparities of miR–target regulation across cancers. In the first half of this review, we explore methods for identification of pairwise MTIs, and in the second half, we explore computational tools for discovery of miR-regulatory modules in a cancer-specific and pan-cancer context. We highlight strengths and limitations of each of these tools as a practical guide for the computational biologists.

Microarray-based identification of genes associated with cancer progression and prognosis in hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) is the third leading cause of cancer-related deaths. The average survival and 5-year survival rates of HCC patients still remains poor. Thus, there is an urgent need to better understand the mechanisms of cancer progression in HCC and to identify useful biomarkers to predict prognosis.

METHODS

Public data portals including Oncomine, The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) profiles were used to retrieve the HCC-related microarrays and to identify potential genes contributed to cancer progression. Bioinformatics analyses including pathway enrichment, protein/gene interaction and text mining were used to explain the potential roles of the identified genes in HCC. Quantitative real-time polymerase chain reaction analysis and Western blotting were used to measure the expression of the targets. The data were analysed by SPSS 20.0 software.

RESULTS

We identified 80 genes that were significantly dysregulated in HCC according to four independent microarrays covering 386 cases of HCC and 327 normal liver tissues. Twenty genes were consistently and stably dysregulated in the four microarrays by at least 2-fold and detection of gene expression by RT-qPCR and western blotting showed consistent expression profiles in 11 HCC tissues compared with corresponding paracancerous tissues. Eleven of these 20 genes were associated with disease-free survival (DFS) or overall survival (OS) in a cohort of 157 HCC patients, and eight genes were associated with tumour pathologic PT, tumour stage or vital status. Potential roles of those 20 genes in regulation of HCC progression were predicted, primarily in association with metastasis. INTS8 was specifically correlated with most clinical characteristics including DFS, OS, stage, metastasis, invasiveness, diagnosis, and age.

CONCLUSION

The significantly dysregulated genes identified in this study were associated with cancer progression and prognosis in HCC, and might be potential therapeutic targets for HCC treatment or potential biomarkers for diagnosis and prognosis.

TCGA Workflow: Analyze cancer genomics and epigenomics data using Bioconductor packages

Biotechnological advances in sequencing have led to an explosion of publicly available data via large international consortia such as The Cancer Genome Atlas (TCGA), The Encyclopedia of DNA Elements (ENCODE), and The NIH Roadmap Epigenomics Mapping Consortium (Roadmap). These projects have provided unprecedented opportunities to interrogate the epigenome of cultured cancer cell lines as well as normal and tumor tissues with high genomic resolution. The Bioconductor project offers more than 1,000 open-source software and statistical packages to analyze high-throughput genomic data. However, most packages are designed for specific data types (e.g. expression, epigenetics, genomics) and there is no one comprehensive tool that provides a complete integrative analysis of the resources and data provided by all three public projects. A need to create an integration of these different analyses was recently proposed. In this workflow, we provide a series of biologically focused integrative analyses of different molecular data. We describe how to download, process and prepare TCGA data and by harnessing several key Bioconductor packages, we describe how to extract biologically meaningful genomic and epigenomic data. Using Roadmap and ENCODE data, we provide a work plan to identify biologically relevant functional epigenomic elements associated with cancer. To illustrate our workflow, we analyzed two types of brain tumors: low-grade glioma (LGG) versus high-grade glioma (glioblastoma multiform or GBM). This workflow introduces the following Bioconductor packages: AnnotationHub, ChIPSeeker, ComplexHeatmap, pathview, ELMER, GAIA, MINET, RTCGAToolbox,  TCGAbiolinks.

TCGA Workflow: Analyze cancer genomics and epigenomics data using Bioconductor packages

Biotechnological advances in sequencing have led to an explosion of publicly available data via large international consortia such as The Cancer Genome Atlas (TCGA), The Encyclopedia of DNA Elements (ENCODE), and The NIH Roadmap Epigenomics Mapping Consortium (Roadmap). These projects have provided unprecedented opportunities to interrogate the epigenome of cultured cancer cell lines as well as normal and tumor tissues with high genomic resolution. The Bioconductor project offers more than 1,000 open-source software and statistical packages to analyze high-throughput genomic data. However, most packages are designed for specific data types (e.g. expression, epigenetics, genomics) and there is no one comprehensive tool that provides a complete integrative analysis of the resources and data provided by all three public projects. A need to create an integration of these different analyses was recently proposed. In this workflow, we provide a series of biologically focused integrative analyses of different molecular data. We describe how to download, process and prepare TCGA data and by harnessing several key Bioconductor packages, we describe how to extract biologically meaningful genomic and epigenomic data. Using Roadmap and ENCODE data, we provide a work plan to identify biologically relevant functional epigenomic elements associated with cancer. To illustrate our workflow, we analyzed two types of brain tumors: low-grade glioma (LGG) versus high-grade glioma (glioblastoma multiform or GBM). This workflow introduces the following Bioconductor packages: AnnotationHub, ChIPSeeker, ComplexHeatmap, pathview, ELMER, GAIA, MINET, RTCGAToolbox,  TCGAbiolinks.

Leveraging big data to transform target selection and drug discovery

The advances of genomics, sequencing, and high throughput technologies have led to the creation of large volumes of diverse datasets for drug discovery. Analyzing these datasets to better understand disease and discover new drugs is becoming more common. Recent open data initiatives in basic and clinical research have dramatically increased the types of data available to the public. The past few years have witnessed successful use of big data in many sectors across the whole drug discovery pipeline. In this review, we will highlight the state of the art in leveraging big data to identify new targets, drug indications, and drug response biomarkers in this era of precision medicine.