A powerful framework for an integrative study with heterogeneous omics data: from univariate statistics to multi-block analysis

Systems-level computational modeling in ischemic stroke: from cells to patients

Ischemic stroke, a significant threat to human life and health, refers to a class of conditions where brain tissue damage is induced following decreased cerebral blood flow. The incidence of ischemic stroke has been steadily increasing globally, and its disease mechanisms are highly complex and involve a multitude of biological mechanisms at various scales from genes all the way to the human body system that can affect the stroke onset, progression, treatment, and prognosis. To complement conventional experimental research methods, computational systems biology modeling can integrate and describe the pathogenic mechanisms of ischemic stroke across multiple biological scales and help identify emergent modulatory principles that drive disease progression and recovery. In addition, by running virtual experiments and trials in computers, these models can efficiently predict and evaluate outcomes of different treatment methods and thereby assist clinical decision-making. In this review, we summarize the current research and application of systems-level computational modeling in the field of ischemic stroke from the multiscale mechanism-based, physics-based and omics-based perspectives and discuss how modeling-driven research frameworks can deliver insights for future stroke research and drug development.

Integrative multi-omics analysis of genomic, epigenomic, and metabolomics data leads to new insights for Attention-Deficit/Hyperactivity Disorder

The evolving field of multi-omics combines data and provides methods for simultaneous analysis across several omics levels. Here, we integrated genomics (transmitted and non-transmitted polygenic scores [PGSs]), epigenomics, and metabolomics data in a multi-omics framework to identify biomarkers for Attention-Deficit/Hyperactivity Disorder (ADHD) and investigated the connections among the three omics levels. We first trained single- and next multi-omics models to differentiate between cases and controls in 596 twins (cases = 14.8%) from the Netherlands Twin Register (NTR) demonstrating reasonable in-sample prediction through cross-validation. The multi-omics model selected 30 PGSs, 143 CpGs, and 90 metabolites. We confirmed previous associations of ADHD with glucocorticoid exposure and the transmembrane protein family TMEM, show that the DNA methylation of the MAD1L1 gene associated with ADHD has a relation with parental smoking behavior, and present novel findings including associations between indirect genetic effects and CpGs of the STAP2 gene. However, out-of-sample prediction in NTR participants (N = 258, cases = 14.3%) and in a clinical sample (N = 145, cases = 51%) did not perform well (range misclassification was [0.40, 0.57]). The results highlighted connections between omics levels, with the strongest connections between non-transmitted PGSs, CpGs, and amino acid levels and show that multi-omics designs considering interrelated omics levels can help unravel the complex biology underlying ADHD.

Genomic data integration tutorial, a plant case study

BACKGROUND

The ongoing evolution of the Next Generation Sequencing (NGS) technologies has led to the production of genomic data on a massive scale. While tools for genomic data integration and analysis are becoming increasingly available, the conceptual and analytical complexities still represent a great challenge in many biological contexts.

RESULTS

To address this issue, we describe a six-steps tutorial for the best practices in genomic data integration, consisting of (1) designing a data matrix; (2) formulating a specific biological question toward data description, selection and prediction; (3) selecting a tool adapted to the targeted questions; (4) preprocessing of the data; (5) conducting preliminary analysis, and finally (6) executing genomic data integration.

CONCLUSION

The tutorial has been tested and demonstrated on publicly available genomic data generated from poplar (Populus L.), a woody plant model. We also developed a new graphical output for the unsupervised multi-block analysis, cimDiablo_v2, available at https://forgemia.inra.fr/umr-gdec/omics-integration-on-poplar , and allowing the selection of master drivers in genomic data variation and interplay.

Maximizing the value of twin studies in health and behaviour

In the classical twin design, researchers compare trait resemblance in cohorts of identical and non-identical twins to understand how genetic and environmental factors correlate with resemblance in behaviour and other phenotypes. The twin design is also a valuable tool for studying causality, intergenerational transmission, and gene-environment correlation and interaction. Here we review recent developments in twin studies, recent results from twin studies of new phenotypes and recent insights into twinning. We ask whether the results of existing twin studies are representative of the general population and of global diversity, and we conclude that stronger efforts to increase representativeness are needed. We provide an updated overview of twin concordance and discordance for major diseases and mental disorders, which conveys a crucial message: genetic influences are not as deterministic as many believe. This has important implications for public understanding of genetic risk prediction tools, as the accuracy of genetic predictions can never exceed identical twin concordance rates.

Integrative Multi-omics Analysis of Childhood Aggressive Behavior

This study introduces and illustrates the potential of an integrated multi-omics approach in investigating the underlying biology of complex traits such as childhood aggressive behavior. In 645 twins (cases = 42%), we trained single- and integrative multi-omics models to identify biomarkers for subclinical aggression and investigated the connections among these biomarkers. Our data comprised transmitted and two non-transmitted polygenic scores (PGSs) for 15 traits, 78,772 CpGs, and 90 metabolites. The single-omics models selected 31 PGSs, 1614 CpGs, and 90 metabolites, and the multi-omics model comprised 44 PGSs, 746 CpGs, and 90 metabolites. The predictive accuracy for these models in the test (N = 277, cases = 42%) and independent clinical data (N = 142, cases = 45%) ranged from 43 to 57%. We observed strong connections between DNA methylation, amino acids, and parental non-transmitted PGSs for ADHD, Autism Spectrum Disorder, intelligence, smoking initiation, and self-reported health. Aggression-related omics traits link to known and novel risk factors, including inflammation, carcinogens, and smoking.

Computational approaches for network-based integrative multi-omics analysis

Advances in omics technologies allow for holistic studies into biological systems. These studies rely on integrative data analysis techniques to obtain a comprehensive view of the dynamics of cellular processes, and molecular mechanisms. Network-based integrative approaches have revolutionized multi-omics analysis by providing the framework to represent interactions between multiple different omics-layers in a graph, which may faithfully reflect the molecular wiring in a cell. Here we review network-based multi-omics/multi-modal integrative analytical approaches. We classify these approaches according to the type of omics data supported, the methods and/or algorithms implemented, their node and/or edge weighting components, and their ability to identify key nodes and subnetworks. We show how these approaches can be used to identify biomarkers, disease subtypes, crosstalk, causality, and molecular drivers of physiological and pathological mechanisms. We provide insight into the most appropriate methods and tools for research questions as showcased around the aetiology and treatment of COVID-19 that can be informed by multi-omics data integration. We conclude with an overview of challenges associated with multi-omics network-based analysis, such as reproducibility, heterogeneity, (biological) interpretability of the results, and we highlight some future directions for network-based integration.

Identification of Tumor Mutation Burden, Microsatellite Instability, and Somatic Copy Number Alteration Derived Nine Gene Signatures to Predict Clinical Outcomes in STAD

Genomic features, including tumor mutation burden (TMB), microsatellite instability (MSI), and somatic copy number alteration (SCNA), had been demonstrated to be involved with the tumor microenvironment (TME) and outcome of gastric cancer (GC). We obtained profiles of TMB, MSI, and SCNA by processing 405 GC data from The Cancer Genome Atlas (TCGA) and then conducted a comprehensive analysis though “iClusterPlus.” A total of two subgroups were generated, with distinguished prognosis, somatic mutation burden, copy number changes, and immune landscape. We revealed that Cluster1 was marked by a better prognosis, accompanied by higher TMB, MSIsensor score, TMEscore, and lower SCNA burden. Based on these clusters, we screened 196 differentially expressed genes (DEGs), which were subsequently projected into univariate Cox survival analysis. We constructed a 9-gene immune risk score (IRS) model using LASSO-penalized logistic regression. Moreover, the prognostic prediction of IRS was verified by receiver operating characteristic (ROC) curve analysis and nomogram plot. Another independent Gene Expression Omnibus (GEO) contained specimens from 109 GC patients was designed as an external validation. Our works suggested that the 9‐gene‐signature prediction model, which was derived from TMB, MSI, and SCNA, was a promising predictive tool for clinical outcomes in GC patients. This novel methodology may help clinicians uncover the underlying mechanisms and guide future treatment strategies.

An integrative Study Showing the Adaptation to Sub-Optimal Growth Conditions of Natural Populations of Arabidopsis thaliana: A Focus on Cell Wall Changes

In the global warming context, plant adaptation occurs, but the underlying molecular mechanisms are poorly described. Studying natural variation of the model plant Arabidopsisthaliana adapted to various environments along an altitudinal gradient should contribute to the identification of new traits related to adaptation to contrasted growth conditions. The study was focused on the cell wall (CW) which plays major roles in the response to environmental changes. Rosettes and floral stems of four newly-described populations collected at different altitudinal levels in the Pyrenees Mountains were studied in laboratory conditions at two growth temperatures (22 vs. 15 °C) and compared to the well-described Col ecotype. Multi-omic analyses combining phenomics, metabolomics, CW proteomics, and transcriptomics were carried out to perform an integrative study to understand the mechanisms of plant adaptation to contrasted growth temperature. Different developmental responses of rosettes and floral stems were observed, especially at the CW level. In addition, specific population responses are shown in relation with their environment and their genetics. Candidate genes or proteins playing roles in the CW dynamics were identified and will deserve functional validation. Using a powerful framework of data integration has led to conclusions that could not have been reached using standard statistical approaches.