Lung Tissue Multi-Layer Network Analysis Uncovers the Molecular Heterogeneity of COPD.

BACKGROUND: Chronic Obstructive Pulmonary Disease (COPD) is a heterogeneous condition. We hypothesized that the unbiased integration of different COPD lung omics using a novel multi-layer approach may unravel mechanisms associated with clinical characteristics. METHODS: We profiled mRNA, miRNA and methylome in lung tissue samples from 135 former smokers with COPD. For each omic (layer) we built a patient network based on molecular similarity. The three networks were used to build a multi-layer network, and optimization of multiplex-modularity was employed to identify patient communities across the three distinct layers. Uncovered communities were related to clinical features. RESULTS: We identified five patient communities in the multi-layer network which were molecularly distinct and related to clinical characteristics, such as FEV1 and blood eosinophils. Two communities (C#3 and C#4) had both similarly low FEV1 values and emphysema, but were molecularly different: C#3, but not C#4, presented B and T cell signatures and a downregulation of secretory (SCGB1A1/SCGB3A1) and ciliated cells. A machine learning model was set up to discriminate C#3 and C#4 in our cohort, and to validate them in an independent cohort. Finally, using spatial transcriptomics we characterized the small airway differences between C#3 and C#4, identifying an upregulation of T/B cell homing chemokines, and bacterial response genes in C#3. CONCLUSIONS: A novel multi-layer network analysis is able to identify clinically relevant COPD patient communities. Patients with similarly low FEV1 and emphysema can have molecularly distinct small airways and immune response patterns, indicating that different endotypes can lead to similar clinical presentation.

Rare disease research workflow using multilayer networks elucidates the molecular determinants of severity in Congenital Myasthenic Syndromes.

Exploring the molecular basis of disease severity in rare disease scenarios is a challenging task provided the limitations on data availability. Causative genes have been described for Congenital Myasthenic Syndromes (CMS), a group of diverse minority neuromuscular junction (NMJ) disorders; yet a molecular explanation for the phenotypic severity differences remains unclear. Here, we present a workflow to explore the functional relationships between CMS causal genes and altered genes from each patient, based on multilayer network community detection analysis of complementary biomedical information provided by relevant data sources, namely protein-protein interactions, pathways and metabolomics. Our results show that CMS severity can be ascribed to the personalized impairment of extracellular matrix components and postsynaptic modulators of acetylcholine receptor (AChR) clustering. This work showcases how coupling multilayer network analysis with personalized -omics information provides molecular explanations to the varying severity of rare diseases; paving the way for sorting out similar cases in other rare diseases.

The PENGUIN approach to reconstruct protein interactions at enhancer-promoter regions and its application to prostate cancer.

We introduce Promoter-Enhancer-Guided Interaction Networks (PENGUIN), a method for studying protein-protein interaction (PPI) networks within enhancer-promoter interactions. PENGUIN integrates H3K27ac-HiChIP data with tissue-specific PPIs to define enhancer-promoter PPI networks (EPINs). We validated PENGUIN using cancer (LNCaP) and benign (LHSAR) prostate cell lines. Our analysis detected EPIN clusters enriched with the architectural protein CTCF, a regulator of enhancer-promoter interactions. CTCF presence was coupled with the prevalence of prostate cancer (PrCa) single nucleotide polymorphisms (SNPs) within the same EPIN clusters, suggesting functional implications in PrCa. Within the EPINs displaying enrichments in both CTCF and PrCa SNPs, we also show enrichment in oncogenes. We substantiated our identified SNPs through CRISPR/Cas9 knockout and RNAi screens experiments. Here we show that PENGUIN provides insights into the intricate interplay between enhancer-promoter interactions and PPI networks, which are crucial for identifying key genes and potential intervention targets. A dedicated server is available at https://penguin.life.bsc.es/ .

Artificial intelligence in cancer research: learning at different levels of data granularity.

From genome-scale experimental studies to imaging data, behavioral footprints, and longitudinal healthcare records, the convergence of big data in cancer research and the advances in Artificial Intelligence (AI) is paving the way to develop a systems view of cancer. Nevertheless, this biomedical area is largely characterized by the co-existence of big data and small data resources, highlighting the need for a deeper investigation about the crosstalk between different levels of data granularity, including varied sample sizes, labels, data types, and other data descriptors. This review introduces the current challenges, limitations, and solutions of AI in the heterogeneous landscape of data granularity in cancer research. Such a variety of cancer molecular and clinical data calls for advancing the interoperability among AI approaches, with particular emphasis on the synergy between discriminative and generative models that we discuss in this work with several examples of techniques and applications.

The multilayer community structure of medulloblastoma.

Multilayer networks allow interpreting the molecular basis of diseases, which is particularly challenging in rare diseases where the number of cases is small compared with the size of the associated multi-omics datasets. In this work, we develop a dimensionality reduction methodology to identify the minimal set of genes that characterize disease subgroups based on their persistent association in multilayer network communities. We use this approach to the study of medulloblastoma, a childhood brain tumor, using proteogenomic data. Our approach is able to recapitulate known medulloblastoma subgroups (accuracy >94%) and provide a clear characterization of gene associations, with the downstream implications for diagnosis and therapeutic interventions. We verified the general applicability of our method on an independent medulloblastoma dataset (accuracy >98%). This approach opens the door to a new generation of multilayer network-based methods able to overcome the specific dimensionality limitations of rare disease datasets.