Multilevel omic data integration in cancer cell lines: advanced annotation and emergent properties

Proteomics appending a complementary dimension to precision oncotherapy

Recent advances in high-throughput proteomic profiling technologies have facilitated the precise quantification of numerous proteins across multiple specimens concurrently. Researchers have the opportunity to comprehensively analyze the molecular signatures in plentiful medical specimens or disease pattern cell lines. Along with advances in data analysis and integration, proteomics data could be efficiently consolidated and employed to recognize precise elementary molecular mechanisms and decode individual biomarkers, guiding the precision treatment of tumors. Herein, we review a broad array of proteomics technologies and the progress and methods for the integration of proteomics data and further discuss how to better merge proteomics in precision medicine and clinical settings.

The downregulation of miR-509-3p expression by collagen type XI alpha 1-regulated hypermethylation facilitates cancer progression and chemoresistance via the DNA methyltransferase 1/Small ubiquitin-like modifier-3 axis in ovarian cancer cells

BACKGROUND

MicroRNAs are a group of small non-coding RNAs that are involved in development and diseases such as cancer. Previously, we demonstrated that miR-335 is crucial for preventing collagen type XI alpha 1 (COL11A1)-mediated epithelial ovarian cancer (EOC) progression and chemoresistance. Here, we examined the role of miR-509-3p in EOC.

METHODS

The patients with EOC who underwent primary cytoreductive surgery and postoperative platinum-based chemotherapy were recruited. Their clinic-pathologic characteristics were collected, and disease-related survivals were determined. The COL11A1 and miR-509-3p mRNA expression levels of 161 ovarian tumors were determined by real-time reverse transcription-polymerase chain reaction. Additionally, miR-509-3p hypermethylation was evaluated by sequencing in these tumors. The A2780CP70 and OVCAR-8 cells transfected with miR-509-3p mimic, while the A2780 and OVCAR-3 cells transfected with miR-509-3p inhibitor. The A2780CP70 cells transfected with a small interference RNA of COL11A1, and the A2780 cells transfected with a COL11A1 expression plasmid. Site-directed mutagenesis, luciferase, and chromatin immunoprecipitation assays were performed in this study.

RESULTS

Low miR-509-3p levels were correlated with disease progression, a poor survival, and high COL11A1 expression levels. In vivo studies reinforced these findings and indicated that the occurrence of invasive EOC cell phenotypes and resistance to cisplatin are decreased by miR-509-3p. The miR-509-3p promoter region (p278) is important for miR-509-3p transcription regulation via methylation. The miR-509-3p hypermethylation frequency was significantly higher in EOC tumors with a low miR-509-3p expression than in those with a high miR-509-3p expression. The patients with miR-509-3p hypermethylation had a significantly shorter overall survival (OS) than those without miR-509-3p hypermethylation. Mechanistic studies further indicated that miR-509-3p transcription was downregulated by COL11A1 through a DNA methyltransferase 1 (DNMT1) stability increase. Moreover, miR-509-3p targets small ubiquitin-like modifier (SUMO)-3 to regulate EOC cell growth, invasiveness, and chemosensitivity.

CONCLUSION

The miR-509-3p/DNMT1/SUMO-3 axis may be an ovarian cancer treatment target.

The Downregulation of miR-509-3p Expression by Collagen Type XI Alpha 1-Regulated Hypermethylation Facilitates Cancer Progression and Chemoresistance via the DNA Methyltransferase 1/Small Ubiquitin-like Modifier-3 Axis in Ovarian Cancer Cells

Background MicroRNAs are a group of small non-coding RNAs that are involved in development and diseases such as cancer. Previously, we demonstrated that miR-335 is crucial for preventing collagen type XI alpha 1 (COL11A1)-mediated epithelial ovarian cancer (EOC) progression and chemoresistance. Here, we examined the role of miR-509-3p in EOC. Methods The patients with EOC who underwent primary cytoreductive surgery and postoperative platinum-based chemotherapy were recruited. Their clinic-pathologic characteristics were collected, and disease-related survivals were determined. The COL11A1 and miR-509-3p mRNA expression levels of 161 ovarian tumors were determined by real-time reverse transcription-polymerase chain reaction. Additionally, miR-509-3p hypermethylation was evaluated by sequencing in these tumors. The A2780CP70 and OVCAR-8 cells transfected with miR-509-3p mimic, while the A2780 and OVCAR-3 cells transfected with miR-509-3p inhibitor. The A2780CP70 cells transfected with a small interference RNA of COL11A1, and the A2780 cells transfected with a COL11A1 expression plasmid. Site-directed mutagenesis, luciferase, and chromatin immunoprecipitation assays were performed in this study. Results Low miR-509-3p levels were correlated with disease progression, a poor survival, and high COL11A1 expression levels. In vivo studies reinforced these findings and indicated that the occurrence of invasive EOC cell phenotypes and resistance to cisplatin are decreased by miR-509-3p. The miR-509-3p promoter region (p278) is important for miR-509-3p transcription regulation via methylation. The miR-509-3p hypermethylation frequency was significantly higher in EOC tumors with a low miR-509-3p expression than in those with a high miR-509-3p expression. The patients with miR-509-3p hypermethylation had a significantly shorter overall survival (OS) than those without miR-509-3p hypermethylation. Mechanistic studies further indicated that miR-509-3p transcription was downregulated by COL11A1 through a DNA methyltransferase 1 (DNMT1) phosphorylation and stability increase. Moreover, miR-509-3p targets small ubiquitin-like modifier (SUMO)-3 to regulate EOC cell growth, invasiveness, and chemosensitivity. Conclusion The miR-509-3p/DNMT1/SUMO-3 axis may be an ovarian cancer treatment target.

Unsupervised Multi-Omics Data Integration Methods: A Comprehensive Review

Through the developments of Omics technologies and dissemination of large-scale datasets, such as those from The Cancer Genome Atlas, Alzheimer’s Disease Neuroimaging Initiative, and Genotype-Tissue Expression, it is becoming increasingly possible to study complex biological processes and disease mechanisms more holistically. However, to obtain a comprehensive view of these complex systems, it is crucial to integrate data across various Omics modalities, and also leverage external knowledge available in biological databases. This review aims to provide an overview of multi-Omics data integration methods with different statistical approaches, focusing on unsupervised learning tasks, including disease onset prediction, biomarker discovery, disease subtyping, module discovery, and network/pathway analysis. We also briefly review feature selection methods, multi-Omics data sets, and resources/tools that constitute critical components for carrying out the integration.

SMSPL: Robust Multimodal Approach to Integrative Analysis of Multiomics Data

With the recent advancement of technologies, it is progressively easier to collect diverse types of genome-wide data. It is commonly expected that by analyzing these data in an integrated way, one can improve the understanding of a complex biological system. Current methods, however, are prone to overfitting heavy noise such that their applications are limited. High noise is one of the major challenges for multiomics data integration. This may be the main cause of overfitting and poor performance in generalization. A sample reweighting strategy is typically used to cope with this problem. In this article, we propose a robust multimodal data integration method, called SMSPL, which can simultaneously predict subtypes of cancers and identify potentially significant multiomics signatures. Especially, the proposed method leverages the linkages between different types of data to interactively recommend high-confidence samples, adopts a new soft weighting scheme to assign weights to the training samples of each type, and then iterates between weights recalculating and classifiers updating. Simulation and five real experiments substantiate the capability of the proposed method for classification and identification of significant multiomics signatures with heavy noise. We expect SMSPL to take a small step in the multiomics data integration and help researchers comprehensively understand the biological process.

Evaluation and comparison of multi-omics data integration methods for cancer subtyping

Computational integrative analysis has become a significant approach in the data-driven exploration of biological problems. Many integration methods for cancer subtyping have been proposed, but evaluating these methods has become a complicated problem due to the lack of gold standards. Moreover, questions of practical importance remain to be addressed regarding the impact of selecting appropriate data types and combinations on the performance of integrative studies. Here, we constructed three classes of benchmarking datasets of nine cancers in TCGA by considering all the eleven combinations of four multi-omics data types. Using these datasets, we conducted a comprehensive evaluation of ten representative integration methods for cancer subtyping in terms of accuracy measured by combining both clustering accuracy and clinical significance, robustness, and computational efficiency. We subsequently investigated the influence of different omics data on cancer subtyping and the effectiveness of their combinations. Refuting the widely held intuition that incorporating more types of omics data always produces better results, our analyses showed that there are situations where integrating more omics data negatively impacts the performance of integration methods. Our analyses also suggested several effective combinations for most cancers under our studies, which may be of particular interest to researchers in omics data analysis.

Cancer is one of the most heterogeneous diseases, characterized by diverse morphological, phenotypic, and genomic profiles between tumors and their subtypes. Identifying cancer subtypes can help patients receive precise treatments. With the development of high-throughput technologies, genomics, epigenomics, and transcriptomics data have been generated for large cancer patient cohorts. It is believed that the more omics data we use, the more accurate identification of cancer subtypes. To examine this assumption, we first constructed three classes of benchmarking datasets to conduct a comprehensive evaluation and comparison of ten representative multi-omics data integration methods for cancer subtyping by considering their accuracy, robustness, and computational efficiency. Then, we investigated the influence of different omics data and their various combinations on the effectiveness of cancer subtyping. Our analyses showed that there are situations where integrating more omics data negatively impacts the performance of integration methods. We hope that our work may help researchers choose a proper method and an effective data combination when identifying cancer subtypes using data integration methods.

The Metallome as a Link Between the "Omes" in Autism Spectrum Disorders

Metal dyshomeostasis plays a significant role in various neurological diseases such as Alzheimer's disease, Parkinson's disease, Autism Spectrum Disorders (ASD), and many more. Like studies investigating the proteome, transcriptome, epigenome, microbiome, etc., for years, metallomics studies have focused on data from their domain, i.e., trace metal composition, only. Still, few have considered the links between other "omes," which may together result in an individual's specific pathologies. In particular, ASD have been reported to have multitudes of possible causal effects. Metallomics data focusing on metal deficiencies and dyshomeostasis can be linked to functions of metalloenzymes, metal transporters, and transcription factors, thus affecting the proteome and transcriptome. Furthermore, recent studies in ASD have emphasized the gut-brain axis, with alterations in the microbiome being linked to changes in the metabolome and inflammatory processes. However, the microbiome and other "omes" are heavily influenced by the metallome. Thus, here, we will summarize the known implications of a changed metallome for other "omes" in the body in the context of "omics" studies in ASD. We will highlight possible connections and propose a model that may explain the so far independently reported pathologies in ASD.

Candidate gene screening for lipid deposition using combined transcriptomic and proteomic data from Nanyang black pigs

Background

Lower selection intensities in indigenous breeds of Chinese pig have resulted in obvious genetic and phenotypic divergence. One such breed, the Nanyang black pig, is renowned for its high lipid deposition and high genetic divergence, making it an ideal model in which to investigate lipid position trait mechanisms in pigs. An understanding of lipid deposition in pigs might improve pig meat traits in future breeding and promote the selection progress of pigs through modern molecular breeding techniques. Here, transcriptome and tandem mass tag-based quantitative proteome (TMT)-based proteome analyses were carried out using longissimus dorsi (LD) tissues from individual Nanyang black pigs that showed high levels of genetic variation.

Results

A large population of Nanyang black pigs was phenotyped using multi-production trait indexes, and six pigs were selected and divided into relatively high and low lipid deposition groups. The combined transcriptomic and proteomic data identified 15 candidate genes that determine lipid deposition genetic divergence. Among them, FASN, CAT, and SLC25A20 were the main causal candidate genes. The other genes could be divided into lipid deposition-related genes (BDH2, FASN, CAT, DHCR24, ACACA, GK, SQLE, ACSL4, and SCD), PPARA-centered fat metabolism regulatory factors (PPARA, UCP3), transcription or translation regulators (SLC25A20, PDK4, CEBPA), as well as integrin, structural proteins, and signal transduction-related genes (EGFR).

Conclusions

This multi-omics data set has provided a valuable resource for future analysis of lipid deposition traits, which might improve pig meat traits in future breeding and promote the selection progress in pigs, especially in Nanyang black pigs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07764-2.

Radiomic and radiogenomic modeling for radiotherapy: strategies, pitfalls, and challenges

The power of predictive modeling for radiotherapy outcomes has historically been limited by an inability to adequately capture patient-specific variabilities; however, next-generation platforms together with imaging technologies and powerful bioinformatic tools have facilitated strategies and provided optimism. Integrating clinical, biological, imaging, and treatment-specific data for more accurate prediction of tumor control probabilities or risk of radiation-induced side effects are high-dimensional problems whose solutions could have widespread benefits to a diverse patient population-we discuss technical approaches toward this objective. Increasing interest in the above is specifically reflected by the emergence of two nascent fields, which are distinct but complementary: radiogenomics, which broadly seeks to integrate biological risk factors together with treatment and diagnostic information to generate individualized patient risk profiles, and radiomics, which further leverages large-scale imaging correlates and extracted features for the same purpose. We review classical analytical and data-driven approaches for outcomes prediction that serve as antecedents to both radiomic and radiogenomic strategies. Discussion then focuses on uses of conventional and deep machine learning in radiomics. We further consider promising strategies for the harmonization of high-dimensional, heterogeneous multiomics datasets (panomics) and techniques for nonparametric validation of best-fit models. Strategies to overcome common pitfalls that are unique to data-intensive radiomics are also discussed.

USP19 modulates cancer cell migration and invasion and acts as a novel prognostic marker in patients with early breast cancer

Tumor cell dissemination in cancer patients is associated with a significant reduction in their survival and quality of life. The ubiquitination pathway plays a fundamental role in the maintenance of protein homeostasis both in normal and stressed conditions and its dysregulation has been associated with malignant transformation and invasive potential of tumor cells, thus highlighting its value as a potential therapeutic target. In order to identify novel molecular targets of tumor cell migration and invasion we performed a genetic screen with an shRNA library against ubiquitination pathway-related genes. To this end, we set up a protocol to specifically enrich positive migration regulator candidates. We identified the deubiquitinase USP19 and demonstrated that its silencing reduces the migratory and invasive potential of highly invasive breast cancer cell lines. We extended our investigation in vivo and confirmed that mice injected with USP19 depleted cells display increased tumor-free survival, as well as a delay in the onset of the tumor formation and a significant reduction in the appearance of metastatic foci, indicating that tumor cell invasion and dissemination is impaired. In contrast, overexpression of USP19 increased cell invasiveness both in vitro and in vivo, further validating our findings. More importantly, we demonstrated that USP19 catalytic activity is important for the control of tumor cell migration and invasion, and that its molecular mechanism of action involves LRP6, a Wnt co-receptor. Finally, we showed that USP19 overexpression is a surrogate prognostic marker of distant relapse in patients with early breast cancer. Altogether, these findings demonstrate that USP19 might represent a novel therapeutic target in breast cancer.

Unsupervised subtyping and methylation landscape of pancreatic ductal adenocarcinoma

Pancreatic Ductal Adenocarcinoma (PDAC) is an aggressive form of pancreatic cancer that typically manifests itself at an advanced stage and does not respond to most treatment modalities. The survival rate of a PDAC patient is less than 5%, with a median survival of just a couple of months. A better understanding of the molecular pathology of PDAC is needed to guide research for the development of better clinical treatment modalities for PDAC patients. Gene expression studies performed to date have identified different subtypes of PDAC with prognostic and clinical relevance. Subtypes identified to date are highly heterogeneous since pancreatic cancer is heterogeneous cancer. Tumor microenvironment and stroma constitute a major chunk of PDAC and contribute to the heterogeneity. Better subtyping methods are need of the hour for better prognosis and classification of PDAC for future personalized treatment. In this work, we have performed an integrated analysis of DNA methylation and gene expression datasets to provide better mechanistic and molecular insights into Pancreatic cancers, especially PDAC. The use of varied and diverse datasets has provided valuable insights into different cancer types and can play an integral role in revealing the complex nature of underlying biological mechanisms. We performed subtyping of TCGA-PAAD gene expression and methylation datasets into different subtypes using state-of-the-art normalization methods and unsupervised clustering methods that reveal latent hidden factors, leading to additional insights for subtyping. Differential expression and differential methylation were performed for each of the subtypes obtained from clustering. Our analysis gave a consensus of five cluster solution with relevant pathways like MAPK, MET. The five subtypes corresponded to the tumor and stromal subtypes. This analysis helps in distinguishing and identifying different subtypes based on enriched putative genes. These results help propose novel experimentally-verifiable PDAC subtyping and demonstrate that using varied data sets and integrated methods can contribute to disease prognostication and precision medicine in PDAC treatment.

Integrative Analysis of Transcriptome-Wide Association Study and mRNA Expression Profiles Identifies Candidate Genes Associated With Idiopathic Pulmonary Fibrosis

Idiopathic pulmonary fibrosis (IPF) is a type of scarring lung disease characterized by a chronic, progressive, and irreversible decline in lung function. The genetic basis of IPF remains elusive. A transcriptome-wide association study (TWAS) of IPF was performed by FUSION using gene expression weights of three tissues combined with a large-scale genome-wide association study (GWAS) dataset, totally involving 2,668 IPF cases and 8,591 controls. Significant genes identified by TWAS were then subjected to gene ontology (GO) and pathway enrichment analysis. The overlapped GO terms and pathways between enrichment analysis of TWAS significant genes and differentially expressed genes (DEGs) from the genome-wide mRNA expression profiling of IPF were also identified. For TWAS significant genes, protein–protein interaction (PPI) network and clustering modules analyses were further conducted using STRING and Cytoscape. Overall, TWAS identified a group of candidate genes for IPF under the Bonferroni corrected P value threshold (0.05/14929 = 3.35 × 10^–6), such as DSP (P_TWAS = 1.35 × 10^–29 for lung tissue), MUC5B (P_TWAS = 1.09 × 10^–28 for lung tissue), and TOLLIP (P_TWAS = 1.41 × 10^–15 for whole blood). Pathway enrichment analysis identified multiple candidate pathways, such as herpes simplex infection (P value = 7.93 × 10^–5) and antigen processing and presentation (P value = 6.55 × 10^–5). 38 common GO terms and 8 KEGG pathways shared by enrichment analysis of TWAS significant genes and DEGs were identified. In the PPI network, 14 genes (DYNLL1, DYNC1LI1, DYNLL2, HLA-DRB5, HLA-DPB1, HLA-DQB2, HLA-DQA2, HLA-DQB1, HLA-DRB1, POLR2L, CENPP, CENPK, NUP133, and NUP107) were simultaneously detected by hub gene and module analysis. In conclusion, through integrative analysis of TWAS and mRNA expression profiles, we identified multiple novel candidate genes, GO terms and pathways for IPF, which contributes to the understanding of the genetic mechanism of IPF.

The Role of the Gut Microbiome in Cattle Production and Health: Driver or Passenger?

Ruminant production systems face significant challenges currently, driven by heightened awareness of their negative environmental impact and the rapidly rising global population. Recent findings have underscored how the composition and function of the rumen microbiome are associated with economically valuable traits, including feed efficiency and methane emission. Although omics-based technological advances in the last decade have revolutionized our understanding of host-associated microbial communities, there remains incongruence over the correct approach for analysis of large omic data sets. A global approach that examines host/microbiome interactions in both the rumen and the lower digestive tract is required to harness the full potential of the gastrointestinal microbiome for sustainable ruminant production. This review highlights how the ruminant animal production community may identify and exploit the causal relationships between the gut microbiome and host traits of interest for a practical application of omic data to animal health and production.

Multi-omics integration-a comparison of unsupervised clustering methodologies

With the recent developments in the field of multi-omics integration, the interest in factors such as data preprocessing, choice of the integration method and the number of different omics considered had increased. In this work, the impact of these factors is explored when solving the problem of sample classification, by comparing the performances of five unsupervised algorithms: Multiple Canonical Correlation Analysis, Multiple Co-Inertia Analysis, Multiple Factor Analysis, Joint and Individual Variation Explained and Similarity Network Fusion. These methods were applied to three real data sets taken from literature and several ad hoc simulated scenarios to discuss classification performance in different conditions of noise and signal strength across the data types. The impact of experimental design, feature selection and parameter training has been also evaluated to unravel important conditions that can affect the accuracy of the result.

Graph- and rule-based learning algorithms: a comprehensive review of their applications for cancer type classification and prognosis using genomic data

Cancer is well recognized as a complex disease with dysregulated molecular networks or modules. Graph- and rule-based analytics have been applied extensively for cancer classification as well as prognosis using large genomic and other data over the past decade. This article provides a comprehensive review of various graph- and rule-based machine learning algorithms that have been applied to numerous genomics data to determine the cancer-specific gene modules, identify gene signature-based classifiers and carry out other related objectives of potential therapeutic value. This review focuses mainly on the methodological design and features of these algorithms to facilitate the application of these graph- and rule-based analytical approaches for cancer classification and prognosis. Based on the type of data integration, we divided all the algorithms into three categories: model-based integration, pre-processing integration and post-processing integration. Each category is further divided into four sub-categories (supervised, unsupervised, semi-supervised and survival-driven learning analyses) based on learning style. Therefore, a total of 11 categories of methods are summarized with their inputs, objectives and description, advantages and potential limitations. Next, we briefly demonstrate well-known and most recently developed algorithms for each sub-category along with salient information, such as data profiles, statistical or feature selection methods and outputs. Finally, we summarize the appropriate use and efficiency of all categories of graph- and rule mining-based learning methods when input data and specific objective are given. This review aims to help readers to select and use the appropriate algorithms for cancer classification and prognosis study.

On fusion methods for knowledge discovery from multi-omics datasets

Recent years have witnessed the tendency of measuring a biological sample on multiple omics scales for a comprehensive understanding of how biological activities on varying levels are perturbed by genetic variants, environments, and their interactions. This new trend raises substantial challenges to data integration and fusion, of which the latter is a specific type of integration that applies a uniform method in a scalable manner, to solve biological problems which the multi-omics measurements target. Fusion-based analysis has advanced rapidly in the past decade, thanks to application drivers and theoretical breakthroughs in mathematics, statistics, and computer science. We will briefly address these methods from methodological and mathematical perspectives and categorize them into three types of approaches: data fusion (a narrowed definition as compared to the general data fusion concept), model fusion, and mixed fusion. We will demonstrate at least one typical example in each specific category to exemplify the characteristics, principles, and applications of the methods in general, as well as discuss the gaps and potential issues for future studies.

Multi-omics Data Integration for Identifying Osteoporosis Biomarkers and Their Biological Interaction and Causal Mechanisms

Osteoporosis is characterized by low bone mineral density (BMD). The advancement of high-throughput technologies and integrative approaches provided an opportunity for deciphering the mechanisms underlying osteoporosis. Here, we generated genomic, transcriptomic, methylomic, and metabolomic datasets from 119 subjects with high (n = 61) and low (n = 58) BMDs. By adopting sparse multiple discriminative canonical correlation analysis, we identified an optimal multi-omics biomarker panel with 74 differentially expressed genes (DEGs), 75 differentially methylated CpG sites (DMCs), and 23 differential metabolic products (DMPs). By linking genetic data, we identified 199 targeted BMD-associated expression/methylation/metabolite quantitative trait loci (eQTLs/meQTLs/metaQTLs). The reconstructed networks/pathways showed extensive biomarker interactions, and a substantial proportion of these biomarkers were enriched in RANK/RANKL, MAPK/TGF-β, and WNT/β-catenin pathways and G-protein-coupled receptor, GTP-binding/GTPase, telomere/mitochondrial activities that are essential for bone metabolism. Five biomarkers (FADS2, ADRA2A, FMN1, RABL2A, SPRY1) revealed causal effects on BMD variation. Our study provided an innovative framework and insights into the pathogenesis of osteoporosis.

Proteomics and phosphoproteomics in precision medicine: applications and challenges

Recent advances in proteomics allow the accurate measurement of abundances for thousands of proteins and phosphoproteins from multiple samples in parallel. Therefore, for the first time, we have the opportunity to measure the proteomic profiles of thousands of patient samples or disease model cell lines in a systematic way, to identify the precise underlying molecular mechanism and discover personalized biomarkers, networks and treatments. Here, we review examples of successful use of proteomics and phosphoproteomics data sets in as well as their integration other omics data sets with the aim of precision medicine. We will discuss the bioinformatics challenges posed by the generation, analysis and integration of such large data sets and present potential reasons why proteomics profiling and biomarkers are not currently widely used in the clinical setting. We will finally discuss ways to contribute to the better use of proteomics data in precision medicine and the clinical setting.

Integrated Omics: Tools, Advances, and Future Approaches

With the rapid adoption of high-throughput omic approaches to analyze biological samples such as genomics, transcriptomics, proteomics, and metabolomics, each analysis can generate tera- to peta-byte sized data files on a daily basis. These data file sizes, together with differences in nomenclature among these data types, make the integration of these multi-dimensional omics data into biologically meaningful context challenging. Variously named as integrated omics, multi-omics, poly-omics, trans-omics, pan-omics, or shortened to just 'omics', the challenges include differences in data cleaning, normalization, biomolecule identification, data dimensionality reduction, biological contextualization, statistical validation, data storage and handling, sharing, and data archiving. The ultimate goal is towards the holistic realization of a 'systems biology' understanding of the biological question in hand. Commonly used approaches in these efforts are currently limited by the 3 i's - integration, interpretation, and insights. Post integration, these very large datasets aim to yield unprecedented views of cellular systems at exquisite resolution for transformative insights into processes, events, and diseases through various computational and informatics frameworks. With the continued reduction in costs and processing time for sample analyses, and increasing types of omics datasets generated such as glycomics, lipidomics, microbiomics, and phenomics, an increasing number of scientists in this interdisciplinary domain of bioinformatics face these challenges. We discuss recent approaches, existing tools, and potential caveats in the integration of omics datasets for development of standardized analytical pipelines that could be adopted by the global omics research community.

MOBCdb: a comprehensive database integrating multi-omics data on breast cancer for precision medicine

BACKGROUND

Breast cancer is one of the most frequently diagnosed cancers among women worldwide, characterized by diverse biological heterogeneity. It is well known that complex and combined gene regulation of multi-omics is involved in the occurrence and development of breast cancer.

RESULTS

In this paper, we present the Multi-Omics Breast Cancer Database (MOBCdb), a simple and easily accessible repository that integrates genomic, transcriptomic, epigenomic, clinical, and drug response data of different subtypes of breast cancer. MOBCdb allows users to retrieve simple nucleotide variation (SNV), gene expression, microRNA expression, DNA methylation, and specific drug response data by various search fashions. The genome-wide browser /navigation facility in MOBCdb provides an interface for visualizing multi-omics data of multi-samples simultaneously. Furthermore, the survival module provides survival analysis for all or some of the samples by using data of three omics. The approved public drugs with genetic variations on breast cancer are also included in MOBCdb.

CONCLUSION

In summary, MOBCdb provides users a unique web interface to the integrated multi-omics data of different subtypes of breast cancer, which enables the users to identify potential novel biomarkers for precision medicine.

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.

Integrative Analysis of Omics Big Data

The diversity and huge omics data take biology and biomedicine research and application into a big data era, just like that popular in human society a decade ago. They are opening a new challenge from horizontal data ensemble (e.g., the similar types of data collected from different labs or companies) to vertical data ensemble (e.g., the different types of data collected for a group of person with match information), which requires the integrative analysis in biology and biomedicine and also asks for emergent development of data integration to address the great changes from previous population-guided to newly individual-guided investigations.Data integration is an effective concept to solve the complex problem or understand the complicate system. Several benchmark studies have revealed the heterogeneity and trade-off that existed in the analysis of omics data. Integrative analysis can combine and investigate many datasets in a cost-effective reproducible way. Current integration approaches on biological data have two modes: one is "bottom-up integration" mode with follow-up manual integration, and the other one is "top-down integration" mode with follow-up in silico integration.This paper will firstly summarize the combinatory analysis approaches to give candidate protocol on biological experiment design for effectively integrative study on genomics and then survey the data fusion approaches to give helpful instruction on computational model development for biological significance detection, which have also provided newly data resources and analysis tools to support the precision medicine dependent on the big biomedical data. Finally, the problems and future directions are highlighted for integrative analysis of omics big data.

Emergence of Intrahepatic Cholangiocarcinoma: How High-Throughput Technologies Expedite the Solutions for a Rare Cancer Type

Intrahepatic cholangiocarcinoma (ICC) is the cancer of the intrahepatic bile ducts, and together with hepatocellular carcinoma (HCC), constitute the majority of primary liver cancers. ICC is a rare disorder as its overall incidence is < 1/100,000 in the United States and Europe. However, it shows much higher incidence in particular geographical regions, such as northeastern Thailand, where liver fluke infection is the most common risk factor of ICC. Since the early stages of ICC are often asymptomatic, the patients are usually diagnosed at advanced stages with no effective treatments available, leading to the high mortality rate. In addition, unclear genetic mechanisms, heterogeneous nature, and various etiologies complicate the development of new efficient treatments. Recently, a number of studies have employed high-throughput approaches, including next-generation sequencing and mass spectrometry, in order to understand ICC in different biological aspects. In general, the majority of recurrent genetic alterations identified in ICC are enriched in known tumor suppressor genes and oncogenes, such as mutations in TP53, KRAS, BAP1, ARID1A, IDH1, IDH2, and novel FGFR2 fusion genes. Yet, there are no major driver genes with immediate clinical solutions characterized. Interestingly, recent studies utilized multi-omics data to classify ICC into two main subgroups, one with immune response genes as the main driving factor, while another is enriched with driver mutations in the genes associated with epigenetic regulations, such as IDH1 and IDH2. The two subgroups also show different hypermethylation patterns in the promoter regions. Additionally, the immune response induced by host-pathogen interactions, i.e., liver fluke infection, may further stimulate tumor growth through alterations of the tumor microenvironment. For in-depth functional studies, although many ICC cell lines have been globally established, these homogeneous cell lines may not fully explain the highly heterogeneous genetic contents of this disorder. Therefore, the advent of patient-derived xenograft and 3D patient-derived organoids as new disease models together with the understanding of evolution and genetic alterations of tumor cells at the single-cell resolution will likely become the main focus to fill the current translational research gaps of ICC in the future.

miR-509-3p is clinically significant and strongly attenuates cellular migration and multi-cellular spheroids in ovarian cancer

Ovarian cancer presents as an aggressive, advanced stage cancer with widespread metastases that depend primarily on multicellular spheroids in the peritoneal fluid. To identify new druggable pathways related to metastatic progression and spheroid formation, we integrated microRNA and mRNA sequencing data from 293 tumors from The Cancer Genome Atlas (TCGA) ovarian cancer cohort. We identified miR-509-3p as a clinically significant microRNA that is more abundant in patients with favorable survival in both the TCGA cohort (P = 2.3E–3), and, by in situ hybridization (ISH), in an independent cohort of 157 tumors (P < 1.0E–3). We found that miR-509-3p attenuated migration and disrupted multi-cellular spheroids in HEYA8, OVCAR8, SKOV3, OVCAR3, OVCAR4 and OVCAR5 cell lines. Consistent with disrupted spheroid formation, in TCGA data miR-509-3p's most strongly anti-correlated predicted targets were enriched in components of the extracellular matrix (ECM). We validated the Hippo pathway effector YAP1 as a direct miR-509-3p target. We showed that siRNA to YAP1 replicated 90% of miR-509-3p-mediated migration attenuation in OVCAR8, which contained high levels of YAP1 protein, but not in the other cell lines, in which levels of this protein were moderate to low. Our data suggest that the miR-509-3p/YAP1 axis may be a new druggable target in cancers with high YAP1, and we propose that therapeutically targeting the miR-509-3p/YAP1/ECM axis may disrupt early steps in multi-cellular spheroid formation, and so inhibit metastasis in epithelial ovarian cancer and potentially in other cancers.

mixOmics: An R package for 'omics feature selection and multiple data integration

The advent of high throughput technologies has led to a wealth of publicly available 'omics data coming from different sources, such as transcriptomics, proteomics, metabolomics. Combining such large-scale biological data sets can lead to the discovery of important biological insights, provided that relevant information can be extracted in a holistic manner. Current statistical approaches have been focusing on identifying small subsets of molecules (a 'molecular signature') to explain or predict biological conditions, but mainly for a single type of 'omics. In addition, commonly used methods are univariate and consider each biological feature independently. We introduce mixOmics, an R package dedicated to the multivariate analysis of biological data sets with a specific focus on data exploration, dimension reduction and visualisation. By adopting a systems biology approach, the toolkit provides a wide range of methods that statistically integrate several data sets at once to probe relationships between heterogeneous 'omics data sets. Our recent methods extend Projection to Latent Structure (PLS) models for discriminant analysis, for data integration across multiple 'omics data or across independent studies, and for the identification of molecular signatures. We illustrate our latest mixOmics integrative frameworks for the multivariate analyses of 'omics data available from the package.

Cancer as an ecomolecular disease and a neoplastic consortium

Current anticancer paradigms largely target driver mutations considered integral for cancer cell survival and tumor progression. Although initially successful, many of these strategies are unable to overcome the tremendous heterogeneity that characterizes advanced tumors, resulting in the emergence of resistant disease. Cancer is a rapidly evolving, multifactorial disease that accumulates numerous genetic and epigenetic alterations. This results in wide phenotypic and molecular heterogeneity within the tumor, the complexity of which is further amplified through specific interactions between cancer cells and the tumor microenvironment. In this context, cancer may be perceived as an "ecomolecular" disease that involves cooperation between several neoplastic clones and their interactions with immune cells, stromal fibroblasts, and other cell types present in the microenvironment. This collaboration is mediated by a variety of secreted factors. Cancer is therefore analogous to complex ecosystems such as microbial consortia. In the present article, we comment on the current paradigms and perspectives guiding the development of cancer diagnostics and therapeutics and the potential application of systems biology to untangle the complexity of neoplasia. In our opinion, conceptualization of neoplasia as an ecomolecular disease is warranted. Advances in knowledge pertinent to the complexity and dynamics of interactions within the cancer ecosystem are likely to improve understanding of tumor etiology, pathogenesis, and progression. This knowledge is anticipated to facilitate the design of new and more effective therapeutic approaches that target the tumor ecosystem in its entirety.

Multidimensional Integrative Genomics Approaches to Dissecting Cardiovascular Disease

Elucidating the mechanisms of complex diseases such as cardiovascular disease (CVD) remains a significant challenge due to multidimensional alterations at molecular, cellular, tissue, and organ levels. To better understand CVD and offer insights into the underlying mechanisms and potential therapeutic strategies, data from multiple omics types (genomics, epigenomics, transcriptomics, metabolomics, proteomics, microbiomics) from both humans and model organisms have become available. However, individual omics data types capture only a fraction of the molecular mechanisms. To address this challenge, there have been numerous efforts to develop integrative genomics methods that can leverage multidimensional information from diverse data types to derive comprehensive molecular insights. In this review, we summarize recent methodological advances in multidimensional omics integration, exemplify their applications in cardiovascular research, and pinpoint challenges and future directions in this incipient field.

Identifying the Metabolic Differences of a Fast-Growth Phenotype in Synechococcus UTEX 2973

The photosynthetic capabilities of cyanobacteria make them interesting candidates for industrial bioproduction. One obstacle to large-scale implementation of cyanobacteria is their limited growth rates as compared to industrial mainstays. Synechococcus UTEX 2973, a strain closely related to Synechococcus PCC 7942, was recently identified as having the fastest measured growth rate among cyanobacteria. To facilitate the development of 2973 as a model organism we developed in this study the genome-scale metabolic model iSyu683. Experimental data were used to define CO₂ uptake rates as well as the biomass compositions for each strain. The inclusion of constraints based on experimental measurements of CO₂ uptake resulted in a ratio of the growth rates of Synechococcus 2973 to Synechococcus 7942 of 2.03, which nearly recapitulates the in vivo growth rate ratio of 2.13. This identified the difference in carbon uptake rate as the main factor contributing to the divergent growth rates. Additionally four SNPs were identified as possible contributors to modified kinetic parameters of metabolic enzymes and candidates for further study. Comparisons against more established cyanobacterial strains identified a number of differences between the strains along with a correlation between the number of cytochrome c oxidase operons and heterotrophic or diazotrophic capabilities.

Informatics for Precision Medicine and Healthcare

The past decade has witnessed great advances in biomedical informatics. Biomedical informatics is an emerging field of healthcare that aims to translate the laboratory observation into clinical practice. Smart healthcare has also developed rapidly with ubiquitous sensor and communication technologies. It is able to capture the online patient-centric phenotypic variables, thus providing a rich information base for translational biomedical informatics. Biomedical informatics and smart healthcare represent two interrelated disciplines. On one hand, biomedical informatics translates the bench discoveries into bedside, and, on the other hand, it is reciprocally informed by clinical data generated from smart healthcare. In this chapter, we will introduce the major strategies and challenges in the application of biomedical informatics technology in precision medicine and healthcare. We highlight how the informatics technology will promote the precision medicine and therefore promise the improvement of healthcare.

How to make DNA methylome wide association studies more powerful

Genome-wide association studies had a troublesome adolescence, while researchers increased statistical power, in part by increasing subject numbers. Interrogating the interaction of genetic and environmental influences raised new challenges of statistical power, which were not easily bested by the addition of subjects. Screening the DNA methylome offers an attractive alternative as methylation can be thought of as a proxy for the combined influences of genetics and environment. There are statistical challenges unique to DNA methylome data and also multiple features, which can be exploited to increase power. We anticipate the development of DNA methylome association study designs and new analytical methods, together with integration of data from other molecular species and other studies, which will boost statistical power and tackle causality. In this way, the molecular trajectories that underlie disease development will be uncovered.

Methods for the integration of multi-omics data: mathematical aspects

BACKGROUND

Methods for the integrative analysis of multi-omics data are required to draw a more complete and accurate picture of the dynamics of molecular systems. The complexity of biological systems, the technological limits, the large number of biological variables and the relatively low number of biological samples make the analysis of multi-omics datasets a non-trivial problem.

RESULTS AND CONCLUSIONS

We review the most advanced strategies for integrating multi-omics datasets, focusing on mathematical and methodological aspects.

Current dichotomy between traditional molecular biological and omic research in cancer biology and pharmacology

There is currently a split within the cancer research community between traditional molecular biological hypothesis-driven and the more recent “omic” forms or research. While the molecular biological approach employs the tried and true single alteration-single response formulations of experimentation, the omic employs broad-based assay or sample collection approaches that generate large volumes of data. How to integrate the benefits of these two approaches in an efficient and productive fashion remains an outstanding issue. Ideally, one would merge the understandability, exactness, simplicity, and testability of the molecular biological approach, with the larger amounts of data, simultaneous consideration of multiple alterations, consideration of genes both of known interest along with the novel, cross-sample comparisons among cell lines and patient samples, and consideration of directed questions while simultaneously gaining exposure to the novel provided by the omic approach. While at the current time integration of the two disciplines remains problematic, attempts to do so are ongoing, and will be necessary for the understanding of the large cell line screens including the Developmental Therapeutics Program’s NCI-60, the Broad Institute’s Cancer Cell Line Encyclopedia, and the Wellcome Trust Sanger Institute’s Cancer Genome Project, as well as the the Cancer Genome Atlas clinical samples project. Going forward there is significant benefit to be had from the integration of the molecular biological and the omic forms or research, with the desired goal being improved translational understanding and application.

In vitro human cell line models to predict clinical response to anticancer drugs

In vitro human cell line models have been widely used for cancer pharmacogenomic studies to predict clinical response, to help generate pharmacogenomic hypothesis for further testing, and to help identify novel mechanisms associated with variation in drug response. Among cell line model systems, immortalized cell lines such as Epstein-Barr virus (EBV)-transformed lymphoblastoid cell lines (LCLs) have been used most often to test the effect of germline genetic variation on drug efficacy and toxicity. Another model, especially in cancer research, uses cancer cell lines such as the NCI-60 panel. These models have been used mainly to determine the effect of somatic alterations on response to anticancer therapy. Even though these cell line model systems are very useful for initial screening, results from integrated analyses of multiple omics data and drug response phenotypes using cell line model systems still need to be confirmed by functional validation and mechanistic studies, as well as validation studies using clinical samples. Future models might include the use of patient-specific inducible pluripotent stem cells and the incorporation of 3D culture which could further optimize in vitro cell line models to improve their predictive validity.

Investigating the utility of clinical outcome-guided mutual information network in network-based Cox regression

BACKGROUND

Network-based approaches have recently gained considerable popularity in high- dimensional regression settings. For example, the Cox regression model is widely used in expression analysis to predict the survival of patients. However, as the number of genes becomes substantially larger than the number of samples, the traditional Cox or L2-regularized Cox models are still prone to noise and produce unreliable estimations of regression coefficients. A recent approach called the network-based Cox (Net-Cox) model attempts to resolve this issue by incorporating prior gene network information into the Cox regression. The Net-Cox model has shown to outperform the models that do not use this network information.

RESULTS

In this study, we demonstrate an alternative network construction method for the outcome-guided gene interaction network, and we investigate its utility in survival analysis using Net-Cox regression as compared with conventional networks, such as co-expression or static networks obtained from the existing knowledgebase. Our network edges consist of gene pairs that are significantly associated with the clinical outcome. We measure the strength of this association using mutual information between the gene pair and the clinical outcome. We applied this approach to ovarian cancer patients' data in The Cancer Genome Atlas (TCGA) and compared the predictive performance of the proposed approach with those that use other types of networks.

CONCLUSIONS

We found that the alternative outcome-guided mutual information network further improved the prediction power of the network-based Cox regression. We expect that a modification of the network regularization term in the Net-Cox model could further improve its prediction power because the properties of our network edges are not optimally reflected in its current form.

Identifying and validating a combined mRNA and microRNA signature in response to imatinib treatment in a chronic myeloid leukemia cell line

Imatinib, a targeted tyrosine kinase inhibitor, is the gold standard for managing chronic myeloid leukemia (CML). Despite its wide application, imatinib resistance occurs in 20-30% of individuals with CML. Multiple potential biomarkers have been identified to predict imatinib response; however, the majority of them remain externally uncorroborated. In this study, we set out to systematically identify gene/microRNA (miRNA) whose expression changes are related to imatinib response. Through a Gene Expression Omnibus search, we identified two genome-wide expression datasets that contain expression changes in response to imatinib treatment in a CML cell line (K562): one for mRNA and the other for miRNA. Significantly differentially expressed transcripts/miRNAs post imatinib treatment were identified from both datasets. Three additional filtering criteria were applied 1) miRbase/miRanda predictive algorithm; 2) opposite direction of imatinib effect for genes and miRNAs; and 3) literature support. These criteria narrowed our candidate gene-miRNA to a single pair: IL8 and miR-493-5p. Using PCR we confirmed the significant up-regulation and down-regulation of miR-493-5p and IL8 by imatinib treatment, respectively in K562 cells. In addition, IL8 expression was significantly down-regulated in K562 cells 24 hours after miR-493-5p mimic transfection (p = 0.002). Furthermore, we demonstrated significant cellular growth inhibition after IL8 inhibition through either gene silencing or by over-expression of miR-493-5p (p = 0.0005 and p = 0.001 respectively). The IL8 inhibition also further sensitized K562 cells to imatinib cytotoxicity (p < 0.0001). Our study combined expression changes in transcriptome and miRNA after imatinib exposure to identify a potential gene-miRNA pair that is a critical target in imatinib response. Experimental validation supports the relationships between IL8 and miR-493-5p and between this gene-miRNA pair and imatinib sensitivity in a CML cell line. Our data suggests integrative analysis of multiple omic level data may provide new insight into biomarker discovery as well as mechanisms of imatinib resistance.

MicroRNA-205 targets tight junction-related proteins during urothelial cellular differentiation

The mammalian bladder urothelium classified as basal, intermediate, and terminally differentiated umbrella cells offers one of the most effective permeability barrier functions known to exist in nature because of the formation of apical uroplakin plaques and tight junctions. To improve our understanding of urothelial differentiation, we analyzed the microRNA (miRNA) expression profiles of mouse urinary tissues and by TaqMan miRNA analysis of microdissected urothelial layers and in situ miRNA-specific hybridization to determine the dependence of these miRNAs on the differentiation stage. Our in situ hybridization studies revealed that miR-205 was enriched in the undifferentiated basal and intermediate cell layers. We then used a quantitative proteomics approach to identify miR-205 target genes in primary cultured urothelial cells subjected to antagomir-mediated knockdown of specific miRNAs. Twenty-four genes were reproducibly regulated by miR-205; eleven of them were annotated as cell junction- and tight junction-related molecules. Western blot analysis demonstrated that antagomir-induced silencing of miR-205 in primary cultured urothelial cells elevated the expression levels of Tjp1, Cgnl1, and Cdc42. Ectopic expression of miR-205 in MDCK cells inhibited the expression of tight junction proteins and the formation of tight junctions. miR-205- knockdown urothelial cells showed alterations in keratin synthesis and increases of uroplakin Ia and Ib, which are the urothelial differentiation products. These results suggest that miR-205 may contribute a role in regulation of urothelial differentiation by modulating the expression of tight junction-related molecules.

A mathematical approach to emergent properties of metabolic networks: partial coupling relations, hyperarcs and flux ratios

Emergent properties in systems biology are those which arise only when the biological system passes a certain level of complexity. In this study, we introduce some of the emergent properties which appear in the constraint-based analysis of metabolic networks. These properties generally appear as a result of existence of hfdeyperarcs and irreversible reactions in networks. Here, we present examples of metabolic networks in which there exist at least two reactions whose fluxes cannot be written as products and/or ratios of the stoichiometric coefficients of the network. We show that any such network contains at least one hyperarc. Additionally, we prove that partial coupling cannot appear in consistent metabolic networks with less than four reactions, or with less than three irreversible reactions, or without hyperarc(s).

MicroRNAs: master regulators of drug resistance, stemness, and metastasis

MicroRNAs (miRNAs) are 20-22 nucleotides long small non-coding RNAs that regulate gene expression post-transcriptionally. Last decade has witnessed emerging evidences of active roles of miRNAs in tumor development, progression, metastasis, and drug resistance. Many factors contribute to their dysregulation in cancer, such as chromosomal aberrations, differential methylation of their own or host genes' promoters and alterations in miRNA biogenesis pathways. miRNAs have been shown to act as tumor suppressors or oncogenes depending on the targets they regulate and the tissue where they are expressed. Because miRNAs can regulate dozens of genes simultaneously and they can function as tumor suppressors or oncogenes, they have been proposed as promising targets for cancer therapy. In this review, we focus on the role of miRNAs in driving drug resistance and metastasis which are associated with stem cell properties of cancer cells. Furthermore, we discuss systems biology approaches to combine experimental and computational methods to study effects of miRNAs on gene or protein networks regulating these processes. Finally, we describe methods to target oncogenic or replace tumor suppressor miRNAs and current delivery strategies to sensitize refractory cells and to prevent metastasis. A holistic understanding of miRNAs' functions in drug resistance and metastasis, which are major causes of cancer-related deaths, and the development of novel strategies to target them efficiently will pave the way towards better translation of miRNAs into clinics and management of cancer therapy.

POMO--Plotting Omics analysis results for Multiple Organisms

Background

Systems biology experiments studying different topics and organisms produce thousands of data values across different types of genomic data. Further, data mining analyses are yielding ranked and heterogeneous results and association networks distributed over the entire genome. The visualization of these results is often difficult and standalone web tools allowing for custom inputs and dynamic filtering are limited.

Results

We have developed POMO (http://pomo.cs.tut.fi), an interactive web-based application to visually explore omics data analysis results and associations in circular, network and grid views. The circular graph represents the chromosome lengths as perimeter segments, as a reference outer ring, such as cytoband for human. The inner arcs between nodes represent the uploaded network. Further, multiple annotation rings, for example depiction of gene copy number changes, can be uploaded as text files and represented as bar, histogram or heatmap rings. POMO has built-in references for human, mouse, nematode, fly, yeast, zebrafish, rice, tomato, Arabidopsis, and Escherichia coli. In addition, POMO provides custom options that allow integrated plotting of unsupported strains or closely related species associations, such as human and mouse orthologs or two yeast wild types, studied together within a single analysis. The web application also supports interactive label and weight filtering. Every iterative filtered result in POMO can be exported as image file and text file for sharing or direct future input.

Conclusions

The POMO web application is a unique tool for omics data analysis, which can be used to visualize and filter the genome-wide networks in the context of chromosomal locations as well as multiple network layouts. With the several illustration and filtering options the tool supports the analysis and visualization of any heterogeneous omics data analysis association results for many organisms. POMO is freely available and does not require any installation or registration.

Systems biology in drug discovery and development

The complexity of human biology makes it challenging to develop safe and effective new medicines. Systems biology omics-based efforts have led to an explosion of high-throughput data and focus is now shifting to the integration of diverse data types to connect molecular and pathway information to predict disease outcomes. Better models of human disease biology, including more integrated network-based models that can accommodate multiple omics data types, as well as more relevant experimental systems, will help predict drug effects in patients, enabling personalized medicine, improvement of the success rate of new drugs in the clinic, and the finding of new uses for existing drugs.

Integrative analysis of longitudinal metabolomics data from a personal multi-omics profile

The integrative personal omics profile (iPOP) is a pioneering study that combines genomics, transcriptomics, proteomics, metabolomics and autoantibody profiles from a single individual over a 14-month period. The observation period includes two episodes of viral infection: a human rhinovirus and a respiratory syncytial virus. The profile studies give an informative snapshot into the biological functioning of an organism. We hypothesize that pathway expression levels are associated with disease status. To test this hypothesis, we use biological pathways to integrate metabolomics and proteomics iPOP data. The approach computes the pathways’ differential expression levels at each time point, while taking into account the pathway structure and the longitudinal design. The resulting pathway levels show strong association with the disease status. Further, we identify temporal patterns in metabolite expression levels. The changes in metabolite expression levels also appear to be consistent with the disease status. The results of the integrative analysis suggest that changes in biological pathways may be used to predict and monitor the disease. The iPOP experimental design, data acquisition and analysis issues are discussed within the broader context of personal profiling.