Multi-omics Networks Research Articles

Smccnet 2.0: a comprehensive tool for multi-omics network inference with shiny visualization

SummarySparse multiple canonical correlation network analysis (SmCCNet) is a machine learning technique for integrating omics data along with a variable of interest (e.g., phenotype of complex disease), and reconstructing multi-omics networks that are specific to this variable. We present the second-generation SmCCNet (SmCCNet 2.0) that adeptly integrates single or multiple omics data types along with a quantitative or binary phenotype of interest. In addition, this new package offers a streamlined setup process that can be configured manually or automatically, ensuring a flexible and user-friendly experience.AvailabilityThis package is available in both CRAN: https://cran.r-project.org/web/packages/SmCCNet/index.html and Github: https://github.com/KechrisLab/SmCCNet under the MIT license. The network visualization tool is available at https://smccnet.shinyapps.io/smccnetnetwork/.

Phenotype prediction using biologically interpretable neural networks on multi-cohort multi-omics data

Integrating multi-omics data into predictive models has the potential to enhance accuracy, which is essential for precision medicine. In this study, we developed interpretable predictive models for multi-omics data by employing neural networks informed by prior biological knowledge, referred to as visible networks. These neural networks offer insights into the decision-making process and can unveil novel perspectives on the underlying biological mechanisms associated with traits and complex diseases. We tested the performance, interpretability and generalizability for inferring smoking status, subject age and LDL levels using genome-wide RNA expression and CpG methylation data from the blood of the BIOS consortium (four population cohorts, Ntotal = 2940). In a cohort-wise cross-validation setting, the consistency of the diagnostic performance and interpretation was assessed. Performance was consistently high for predicting smoking status with an overall mean AUC of 0.95 (95% CI: 0.90–1.00) and interpretation revealed the involvement of well-replicated genes such as AHRR, GPR15 and LRRN3. LDL-level predictions were only generalized in a single cohort with an R2 of 0.07 (95% CI: 0.05–0.08). Age was inferred with a mean error of 5.16 (95% CI: 3.97–6.35) years with the genes COL11A2, AFAP1, OTUD7A, PTPRN2, ADARB2 and CD34 consistently predictive. For both regression tasks, we found that using multi-omics networks improved performance, stability and generalizability compared to interpretable single omic networks. We believe that visible neural networks have great potential for multi-omics analysis; they combine multi-omic data elegantly, are interpretable, and generalize well to data from different cohorts.

Multiomics Analysis Reveals Gut Virome-Bacteria-Metabolite Interactions and Their Associations with Symptoms in Patients with IBS-D.

The gut microbiota is involved in the pathogenesis of diarrhea-predominant irritable bowel syndrome (IBS-D), but few studies have focused on the role of the gut virome in IBS-D. We aimed to explore the characteristics of the gut virome in patients with IBS-D, its interactions with bacteria and metabolites, and the associations between gut multiomics profiles and symptoms. This study enrolled twelve patients with IBS-D and eight healthy controls (HCs). The stool samples were subjected to metavirome sequencing, 16S rRNA gene sequencing, and untargeted metabolomic analysis. The participants completed relevant scales to assess the severity of their gastrointestinal symptoms, depression, and anxiety. The results revealed unique DNA and RNA virome profiles in patients with IBS-D with significant alterations in the abundance of contigs from Siphoviridae, Podoviridae, Microviridae, Picobirnaviridae, and Tombusviridae. Single-omics co-occurrence network analyses demonstrated distinct differences in the gut virus, bacteria, and metabolite network patterns between patients with IBS-D and HCs. Multiomics networks revealed that short-chain fatty acid-producing bacteria occupied more core positions in IBS-D networks, but had fewer links to viruses. Amino acids and their derivatives exhibit unique connectivity patterns and centrality features within the IBS-D network. The gastrointestinal and psychological symptom factors of patients with IBS-D were highly clustered in the symptom-multiomics network compared with those of HCs. Machine learning models based on multiomics data can distinguish IBS-D patients from HCs and predict the scores of gastrointestinal and psychological symptoms. This study provides insights into the interactions among gut viruses, bacteria, metabolites, and clinical symptoms in patients with IBS-D, indicating further classification and personalized treatment for IBS-D.

Multiomic Signatures of Traffic-Related Air Pollution in London Reveal Potential Short-Term Perturbations in Gut Microbiome-Related Pathways.

This randomized crossover study investigated the metabolic and mRNA alterations associated with exposure to high and low traffic-related air pollution (TRAP) in 50 participants who were either healthy or were diagnosed with chronic pulmonary obstructive disease (COPD) or ischemic heart disease (IHD). For the first time, this study combined transcriptomics and serum metabolomics measured in the same participants over multiple time points (2 h before, and 2 and 24 h after exposure) and over two contrasted exposure regimes to identify potential multiomic modifications linked to TRAP exposure. With a multivariate normal model, we identified 78 metabolic features and 53 mRNA features associated with at least one TRAP exposure. Nitrogen dioxide (NO2) emerged as the dominant pollutant, with 67 unique associated metabolomic features. Pathway analysis and annotation of metabolic features consistently indicated perturbations in the tryptophan metabolism associated with NO2 exposure, particularly in the gut-microbiome-associated indole pathway. Conditional multiomics networks revealed complex and intricate mechanisms associated with TRAP exposure, with some effects persisting 24 h after exposure. Our findings indicate that exposure to TRAP can alter important physiological mechanisms even after a short-term exposure of a 2 h walk. We describe for the first time a potential link between NO2 exposure and perturbation of the microbiome-related pathways.

Progress on network modeling and analysis of gut microecology: a review.

The gut microecological network is a complex microbial community within the human body that plays a key role in linking dietary nutrition and host physiology. To understand the complex relationships among microbes and their functions within this community, network analysis has emerged as a powerful tool. By representing the interactions between microbes and their associated omics data as a network, we can gain a comprehensive understanding of the ecological mechanisms that drive the human gut microbiota. In addition, the network-based approach provides a more intuitive analysis of the gut microbiota, simplifying the study of its complex dynamics and interdependencies. This review provides a comprehensive overview of the methods used to construct and analyze networks in the context of gut microecological background. We discuss various types of network modeling approaches, including co-occurrence networks, causal networks, dynamic networks, and multi-omics networks, and describe the analytical techniques used to identify important network properties. We also highlight the challenges and limitations of network modeling in this area, such as data scarcity and heterogeneity, and provide future research directions to overcome these limitations. By exploring these network-based methods, researchers can gain valuable insights into the intricate relationships and functional roles of microbial communities within the gut, ultimately advancing our understanding of the gut microbiota's impact on human health.

Phospholipase A2 regulates autophagy in gouty arthritis: proteomic and metabolomic studies

BackgroundAcute gouty arthritis is inflammatory joint arthritis. Gouty arthritis (GA) involves multiple pathological processes. Deposition of joints by monosodium urate (MSU) crystals has been shown to play a critical role in the injury process. Due to the different effects of MSU stimulation on the joints, the exact changes in the synovial fluid are unknown. We want to explore the changes in proteins and metabolites in the joints of gouty arthritis. Regulating various functional substances in the joint can reduce inflammation and pain symptoms.Methods10 patients with gouty knee arthritis and 10 normal controls were selected from clinical, surgical cases. The biological function of the metabolome was assessed by co-expression network analysis. A molecular network based on metabolomic and proteomic data was constructed to study critical molecules. The fundamental molecular changes in the relevant pathways were then verified by western blot.ResultsProteomic analysis showed that the expressions of proteases Cathepsin B, Cathepsin D, Cathepsin G, and Cathepsin S in synovial fluid patients with gouty arthritis were significantly increased. Enrichment analysis showed a positive correlation between lysosomal and clinical inflammatory cell shape changes. Untargeted metabolomic analysis revealed that lipids and lipoids accumulate, inhibit autophagic flux, and modulate inflammation and immunity in gouty arthritis patients. It was determined that the accumulation of lipid substances such as phospholipase A2 led to the imbalanced state of the autophagy-lysosome complex, and the differentially expressed metabolites of Stearoylcarnitine, Tetradecanoylcarnitine, Palmitoylcarnitine were identified (|log2 fold change|> 1.5, adjusted P value < 0.05 and variable importance in prediction (VIP) > 1.5). The autophagy-lysosomal pathway was found to be associated with gouty knee arthritis. Essential molecular alterations of multi-omics networks in gouty knee arthritis patients compared with normal controls involve acute inflammatory response, exosomes, immune responses, lysosomes, linoleic acid metabolism, and synthesis.ConclusionsComprehensive analysis of proteomic and untargeted metabolomics revealed protein and characteristic metabolite alterations in gouty arthritis, it mainly involves lipids and lipid like molecules, phospholipase A2 and autophagic lysosomes. This study describes the pathological characteristics, pathways, potential predictors and treatment goals of gouty knee arthritis.

Multiomics profiling of human plasma and cerebrospinal fluid reveals ATN-derived networks and highlights causal links in Alzheimer's disease.

This study employed an integrative system and causal inference approach to explore molecular signatures in blood and CSF, the amyloid/tau/neurodegeneration [AT(N)] framework, mild cognitive impairment (MCI) conversion to Alzheimer's disease (AD), and genetic risk for AD. Using the European Medical Information Framework (EMIF)-AD cohort, we measured 696 proteins in cerebrospinal fluid (n=371), 4001 proteins in plasma (n=972), 611 metabolites in plasma (n=696), and genotyped whole-blood (7,778,465 autosomal single nucleotide epolymorphisms, n=936). We investigated associations: molecular modules to AT(N), module hubs with AD Polygenic Risk scores and APOE4 genotypes, molecular hubs to MCI conversion and probed for causality with AD using Mendelian randomization (MR). AT(N) framework associated with protein and lipid hubs. In plasma, Proprotein Convertase Subtilisin/Kexin Type 7 showed evidence for causal associations with AD. AD was causally associated with Reticulocalbin 2 and sphingomyelins, an association driven by the APOE isoform. This study reveals multi-omics networks associated with AT(N) and causal AD molecular candidates.

Integrative multi-omics networks identify PKCδ and DNA-PK as master kinases of glioblastoma subtypes and guide targeted cancer therapy

Despite producing a panoply of potential cancer-specific targets, the proteogenomic characterization of human tumors has yet to demonstrate value for precision cancer medicine. Integrative multi-omics using a machine-learning network identified master kinases responsible for effecting phenotypic hallmarks of functional glioblastoma subtypes. In subtype-matched patient-derived models, we validated PKCδ and DNA-PK as master kinases of glycolytic/plurimetabolic and proliferative/progenitor subtypes, respectively, and qualified the kinases as potent and actionable glioblastoma subtype-specific therapeutic targets. Glioblastoma subtypes were associated with clinical and radiomics features, orthogonally validated by proteomics, phospho-proteomics, metabolomics, lipidomics and acetylomics analyses, and recapitulated in pediatric glioma, breast and lung squamous cell carcinoma, including subtype specificity of PKCδ and DNA-PK activity. We developed a probabilistic classification tool that performs optimally with RNA from frozen and paraffin-embedded tissues, which can be used to evaluate the association of therapeutic response with glioblastoma subtypes and to inform patient selection in prospective clinical trials.

A Practical Guide to Inferring Multi-Omics Networks in Plant Systems.

The inference of gene regulatory networks can reveal molecular connections underlying biological processes and improve our understanding of complex biological phenomena in plants. Many previous network studies have inferred networks using only one type of omics data, such as transcriptomics. However, given more recent work applying multi-omics integration in plant biology, such as combining (phospho)proteomics with transcriptomics, it may be advantageous to integrate multiple omics data types into a comprehensive network prediction. Here, we describe a state-of-the-art approach for integrating multi-omics data with gene regulatory network inference to describe signaling pathways and uncover novel regulators. We detail how to download and process transcriptomics and (phospho)proteomics data for network inference, using an example dataset from the plant hormone signaling field. We provide a step-by-step protocol for inference, visualization, and analysis of an integrative multi-omics network using currently available methods. This chapter serves as an accessible guide for novice and intermediate bioinformaticians to analyze their own datasets and reanalyze published work.

Inferencing Bulk Tumor and Single-Cell Multi-Omics Regulatory Networks for Discovery of Biomarkers and Therapeutic Targets.

There are insufficient accurate biomarkers and effective therapeutic targets in current cancer treatment. Multi-omics regulatory networks in patient bulk tumors and single cells can shed light on molecular disease mechanisms. Integration of multi-omics data with large-scale patient electronic medical records (EMRs) can lead to the discovery of biomarkers and therapeutic targets. In this review, multi-omics data harmonization methods were introduced, and common approaches to molecular network inference were summarized. Our Prediction Logic Boolean Implication Networks (PLBINs) have advantages over other methods in constructing genome-scale multi-omics networks in bulk tumors and single cells in terms of computational efficiency, scalability, and accuracy. Based on the constructed multi-modal regulatory networks, graph theory network centrality metrics can be used in the prioritization of candidates for discovering biomarkers and therapeutic targets. Our approach to integrating multi-omics profiles in a patient cohort with large-scale patient EMRs such as the SEER-Medicare cancer registry combined with extensive external validation can identify potential biomarkers applicable in large patient populations. These methodologies form a conceptually innovative framework to analyze various available information from research laboratories and healthcare systems, accelerating the discovery of biomarkers and therapeutic targets to ultimately improve cancer patient survival outcomes.

Uncovering the gene regulatory network of type 2 diabetes through multi-omic data integration

BackgroundType 2 diabetes (T2D) onset is a complex, organized biological process with multilevel regulation, and its physiopathological mechanisms are yet to be elucidated. This study aims to find out the key drivers and pathways involved in the pathogenesis of T2D through multi-omics analysis.MethodsThe datasets used in the experiments comprise three groups: (1) genomic (2) transcriptomic, and (3) epigenomic categories. Then, a series of bioinformatics technologies including Marker set enrichment analysis (MSEA), weighted key driver analysis (wKDA) was performed to identify key drivers. The hub genes were further verified by the Receiver Operator Characteristic (ROC) Curve analysis, proteomic analysis, and Real-time quantitative polymerase chain reaction (RT-qPCR). The multi-omics network was applied to the Pharmomics pipeline in Mergeomics to identify drug candidates for T2D treatment. Then, we used the drug-gene interaction network to conduct network pharmacological analysis. Besides, molecular docking was performed using AutoDock/Vina, a computational docking program.ResultsModule-gene interaction network was constructed using MSEA, which revealed a significant enrichment of immune-related activities and glucose metabolism. Top 10 key drivers (PSMB9, COL1A1, COL4A1, HLA-DQB1, COL3A1, IRF7, COL5A1, CD74, HLA-DQA1, and HLA-DRB1) were selected by wKDA analysis. Among these, COL5A1, IRF7, CD74, and HLA-DRB1 were verified to have the capability to diagnose T2D, and expression levels of PSMB9 and CD74 had significantly higher in T2D patients. We further predict the co-expression network and transcription factor (TF) binding specificity of the key driver. Besides, based on module interaction networks and key driver networks, 17 compounds are considered to possess T2D-control potential, such as sunitinib.ConclusionsWe identified signature genes, biomolecular processes, and pathways using multi-omics networks. Moreover, our computational network analysis revealed potential novel strategies for pharmacologic interventions of T2D.

Effective Subject Representation based on Multi-omics Disease Networks using Graph Embedding.

The study of complex behavior of biological systems has become increasingly dependent on evolutionary network modeling. In particular, multi-omics networks capture interactions between biomolecules such as proteins and metabolites, providing a basis for predicting relationships between such biomolecules and various phenotypic traits of complex diseases. In this paper, we introduce an integrative framework that given a multi-omics network representing a cohort of subjects, learns expressive representations for network nodes, and combines the learned nodes representations with the biological profiles of individual subjects for enriched representation of the subjects. With extensive empirical evaluation using real-world multi-omics networks, we show that our proposed framework significantly outperforms existing and baseline methods in terms of subject representation accuracy, particularly when the multi-omics network representing the cohort is sparse and structured and therefore, more informative.

Hub Genes in Non-Small Cell Lung Cancer Regulatory Networks.

There are currently no accurate biomarkers for optimal treatment selection in early-stage non-small cell lung cancer (NSCLC). Novel therapeutic targets are needed to improve NSCLC survival outcomes. This study systematically evaluated the association between genome-scale regulatory network centralities and NSCLC tumorigenesis, proliferation, and survival in early-stage NSCLC patients. Boolean implication networks were used to construct multimodal networks using patient DNA copy number variation, mRNA, and protein expression profiles. T statistics of differential gene/protein expression in tumors versus non-cancerous adjacent tissues, dependency scores in in vitro CRISPR-Cas9/RNA interference (RNAi) screening of human NSCLC cell lines, and hazard ratios in univariate Cox modeling of the Cancer Genome Atlas (TCGA) NSCLC patients were correlated with graph theory centrality metrics. Hub genes in multi-omics networks involving gene/protein expression were associated with oncogenic, proliferative potentials and poor patient survival outcomes (p < 0.05, Pearson's correlation). Immunotherapy targets PD1, PDL1, CTLA4, and CD27 were ranked as top hub genes within the 10th percentile in most constructed multi-omics networks. BUB3, DNM1L, EIF2S1, KPNB1, NMT1, PGAM1, and STRAP were discovered as important hub genes in NSCLC proliferation with oncogenic potential. These results support the importance of hub genes in NSCLC tumorigenesis, proliferation, and prognosis, with implications in prioritizing therapeutic targets to improve patient survival outcomes.

Multi-Omics Immune Interaction Networks in Lung Cancer Tumorigenesis, Proliferation, and Survival.

There are currently no effective biomarkers for prognosis and optimal treatment selection to improve non-small cell lung cancer (NSCLC) survival outcomes. This study further validated a seven-gene panel for diagnosis and prognosis of NSCLC using RNA sequencing and proteomic profiles of patient tumors. Within the seven-gene panel, ZNF71 expression combined with dendritic cell activities defined NSCLC patient subgroups (n = 966) with distinct survival outcomes (p = 0.04, Kaplan-Meier analysis). ZNF71 expression was significantly associated with the activities of natural killer cells (p = 0.014) and natural killer T cells (p = 0.003) in NSCLC patient tumors (n = 1016) using Chi-squared tests. Overexpression of ZNF71 resulted in decreased expression of multiple components of the intracellular intrinsic and innate immune systems, including dsRNA and dsDNA sensors. Multi-omics networks of ZNF71 and the intracellular intrinsic and innate immune systems were computed as relevant to NSCLC tumorigenesis, proliferation, and survival using patient clinical information and in-vitro CRISPR-Cas9/RNAi screening data. From these networks, pan-sensitive and pan-resistant genes to 21 NCCN-recommended drugs for treating NSCLC were selected. Based on the gene associations with patient survival and in-vitro CRISPR-Cas9, RNAi, and drug screening data, MEK1/2 inhibitors PD-198306 and U-0126, VEGFR inhibitor ZM-306416, and IGF-1R inhibitor PQ-401 were discovered as potential targeted therapy that may also induce an immune response for treating NSCLC.

DrDimont: explainable drug response prediction from differential analysis of multi-omics networks.

MotivationWhile it has been well established that drugs affect and help patients differently, personalized drug response predictions remain challenging. Solutions based on single omics measurements have been proposed, and networks provide means to incorporate molecular interactions into reasoning. However, how to integrate the wealth of information contained in multiple omics layers still poses a complex problem.ResultsWe present DrDimont, Drug response prediction from Differential analysis of multi-omics networks. It allows for comparative conclusions between two conditions and translates them into differential drug response predictions. DrDimont focuses on molecular interactions. It establishes condition-specific networks from correlation within an omics layer that are then reduced and combined into heterogeneous, multi-omics molecular networks. A novel semi-local, path-based integration step ensures integrative conclusions. Differential predictions are derived from comparing the condition-specific integrated networks. DrDimont’s predictions are explainable, i.e. molecular differences that are the source of high differential drug scores can be retrieved. We predict differential drug response in breast cancer using transcriptomics, proteomics, phosphosite and metabolomics measurements and contrast estrogen receptor positive and receptor negative patients. DrDimont performs better than drug prediction based on differential protein expression or PageRank when evaluating it on ground truth data from cancer cell lines. We find proteomic and phosphosite layers to carry most information for distinguishing drug response.Availability and implementationDrDimont is available on CRAN: https://cran.r-project.org/package=DrDimont.Supplementary information Supplementary data are available at Bioinformatics online.

SUBATOMIC: a SUbgraph BAsed mulTi-OMIcs clustering framework to analyze integrated multi-edge networks

BackgroundRepresenting the complex interplay between different types of biomolecules across different omics layers in multi-omics networks bears great potential to gain a deep mechanistic understanding of gene regulation and disease. However, multi-omics networks easily grow into giant hairball structures that hamper biological interpretation. Module detection methods can decompose these networks into smaller interpretable modules. However, these methods are not adapted to deal with multi-omics data nor consider topological features. When deriving very large modules or ignoring the broader network context, interpretability remains limited. To address these issues, we developed a SUbgraph BAsed mulTi-OMIcs Clustering framework (SUBATOMIC), which infers small and interpretable modules with a specific topology while keeping track of connections to other modules and regulators.ResultsSUBATOMIC groups specific molecular interactions in composite network subgraphs of two and three nodes and clusters them into topological modules. These are functionally annotated, visualized and overlaid with expression profiles to go from static to dynamic modules. To preserve the larger network context, SUBATOMIC investigates statistically the connections in between modules as well as between modules and regulators such as miRNAs and transcription factors. We applied SUBATOMIC to analyze a composite Homo sapiens network containing transcription factor-target gene, miRNA-target gene, protein–protein, homologous and co-functional interactions from different databases. We derived and annotated 5586 modules with diverse topological, functional and regulatory properties. We created novel functional hypotheses for unannotated genes. Furthermore, we integrated modules with condition specific expression data to study the influence of hypoxia in three cancer cell lines. We developed two prioritization strategies to identify the most relevant modules in specific biological contexts: one considering GO term enrichments and one calculating an activity score reflecting the degree of differential expression. Both strategies yielded modules specifically reacting to low oxygen levels.ConclusionsWe developed the SUBATOMIC framework that generates interpretable modules from integrated multi-omics networks and applied it to hypoxia in cancer. SUBATOMIC can infer and contextualize modules, explore condition or disease specific modules, identify regulators and functionally related modules, and derive novel gene functions for uncharacterized genes. The software is available at https://github.com/CBIGR/SUBATOMIC.

Significant Subgraph Detection in Multi-omics Networks for Disease Pathway Identification.

Chronic obstructive pulmonary disease (COPD) is one of the leading causes of death in the United States. COPD represents one of many areas of research where identifying complex pathways and networks of interacting biomarkers is an important avenue toward studying disease progression and potentially discovering cures. Recently, sparse multiple canonical correlation network analysis (SmCCNet) was developed to identify complex relationships between omics associated with a disease phenotype, such as lung function. SmCCNet uses two sets of omics datasets and an associated output phenotypes to generate a multi-omics graph, which can then be used to explore relationships between omics in the context of a disease. Detecting significant subgraphs within this multi-omics network, i.e., subgraphs which exhibit high correlation to a disease phenotype and high inter-connectivity, can help clinicians identify complex biological relationships involved in disease progression. The current approach to identifying significant subgraphs relies on hierarchical clustering, which can be used to inform clinicians about important pathways involved in the disease or phenotype of interest. The reliance on a hierarchical clustering approach can hinder subgraph quality by biasing toward finding more compact subgraphs and removing larger significant subgraphs. This study aims to introduce new significant subgraph detection techniques. In particular, we introduce two subgraph detection methods, dubbed Correlated PageRank and Correlated Louvain, by extending the Personalized PageRank Clustering and Louvain algorithms, as well as a hybrid approach combining the two proposed methods, and compare them to the hierarchical method currently in use. The proposed methods show significant improvement in the quality of the subgraphs produced when compared to the current state of the art.

OmicsNet 2.0: a web-based platform for multi-omics integration and network visual analytics.

Researchers are increasingly seeking to interpret molecular data within a multi-omics context to gain a more comprehensive picture of their study system. OmicsNet (www.omicsnet.ca) is a web-based tool developed to allow users to easily build, visualize, and analyze multi-omics networks to study rich relationships among lists of ‘omics features of interest. Three major improvements have been introduced in OmicsNet 2.0, which include: (i) enhanced network visual analytics with eleven 2D graph layout options and a novel 3D module layout; (ii) support for three new ‘omics types: single nucleotide polymorphism (SNP) list from genetic variation studies; taxon list from microbiome profiling studies, as well as liquid chromatography–mass spectrometry (LC–MS) peaks from untargeted metabolomics; and (iii) measures to improve research reproducibility by coupling R command history with the release of the companion OmicsNetR package, and generation of persistent links to share interactive network views. We performed a case study using the multi-omics data obtained from a recent large-scale investigation on inflammatory bowel disease (IBD) and demonstrated that OmicsNet was able to quickly create meaningful multi-omics context to facilitate hypothesis generation and mechanistic insights.

Interpretation of network-based integration from multi-omics longitudinal data

Multi-omics integration is key to fully understand complex biological processes in an holistic manner. Furthermore, multi-omics combined with new longitudinal experimental design can unreveal dynamic relationships between omics layers and identify key players or interactions in system development or complex phenotypes. However, integration methods have to address various experimental designs and do not guarantee interpretable biological results. The new challenge of multi-omics integration is to solve interpretation and unlock the hidden knowledge within the multi-omics data. In this paper, we go beyond integration and propose a generic approach to face the interpretation problem. From multi-omics longitudinal data, this approach builds and explores hybrid multi-omics networks composed of both inferred and known relationships within and between omics layers. With smart node labelling and propagation analysis, this approach predicts regulation mechanisms and multi-omics functional modules. We applied the method on 3 case studies with various multi-omics designs and identified new multi-layer interactions involved in key biological functions that could not be revealed with single omics analysis. Moreover, we highlighted interplay in the kinetics that could help identify novel biological mechanisms. This method is available as an R package netOmics to readily suit any application.

Systems biomarkers for papillary thyroid cancer prognosis and treatment through multi-omics networks

The identification of biomolecules associated with papillary thyroid cancer (PTC) has upmost importance for the elucidation of the disease mechanism and the development of effective diagnostic and treatment strategies. Despite particular findings in this regard, a holistic analysis encompassing molecular data from different biological levels has been lacking. In the present study, a meta-analysis of four transcriptome datasets was performed to identify gene expression signatures in PTC, and reporter molecules were determined by mapping gene expression data onto three major cellular networks, i.e., transcriptional regulatory, protein-protein interaction, and metabolic networks. We identified 282 common genes that were differentially expressed in all PTC datasets. In addition, six proteins (FYN, JUN, LYN, PML, SIN3A, and RARA), two Erb-B2 receptors (ERBB2 and ERBB4), two cyclin-dependent receptors (CDK1 and CDK2), and three histone deacetylase receptors (HDAC1, HDAC2, and HDAC3) came into prominence as proteomic signatures in addition to several metabolites including lactaldehyde and proline at the metabolome level. Significant associations with calcium and MAPK signaling pathways and transcriptional and post-transcriptional activities of 12 TFs and 110 miRNAs were also observed at the regulatory level. Among them, six miRNAs (miR-30b-3p, miR-15b-5p, let-7a-5p, miR-130b-3p, miR-424-5p, and miR-193b-3p) were associated with PTC for the first time in the literature, and the expression levels of miR-30b-3p, miR-15b-5p, and let-7a-5p were found to be predictive of disease prognosis. Drug repositioning and molecular docking simulations revealed that 5 drugs (prochlorperazine, meclizine, rottlerin, cephaeline, and tretinoin) may be useful in the treatment of PTC. Consequently, we report here biomolecule candidates that may be considered as prognostic biomarkers or potential therapeutic targets for further experimental and clinical trials for PTC.