Integration of Molecular Interactome and Targeted Interaction Analysis to Identify a COPD Disease Network Module

Correlation-based network integration of lung RNA sequencing and DNA methylation data in chronic obstructive pulmonary disease

Chronic Obstructive Pulmonary Disease (COPD) is a heterogeneous, chronic inflammatory process of the lungs and, like other complex diseases, is caused by both genetic and environmental factors. Detailed understanding of the molecular mechanisms of complex diseases requires the study of the interplay among different biomolecular layers, and thus the integration of different omics data types. In this study, we investigated COPD-associated molecular mechanisms through a correlation-based network integration of lung tissue RNA-seq and DNA methylation data of COPD cases (n = 446) and controls (n = 346) derived from the Lung Tissue Research Consortium. First, we performed a SWIM-network based analysis to build separate correlation networks for RNA-seq and DNA methylation data for our case-control study population. Then, we developed a method to integrate the results into a coupled network of differentially expressed and differentially methylated genes to investigate their relationships across both molecular layers. The functional enrichment analysis of the nodes of the coupled network revealed a strikingly significant enrichment in Immune System components, both innate and adaptive, as well as immune-system component communication (interleukin and cytokine-cytokine signaling). Our analysis allowed us to reveal novel putative COPD-associated genes and to analyze their relationships, both at the transcriptomics and epigenomics levels, thus contributing to an improved understanding of COPD pathogenesis.

Breathing new life into the study of COPD with genes identified from genome-wide association studies

COPD is a major cause of morbidity and mortality globally. While the significance of environmental exposures in disease pathogenesis is well established, the functional contribution of genetic factors has only in recent years drawn attention. Notably, many genes associated with COPD risk are also linked with lung function. Because reduced lung function precedes COPD onset, this association is consistent with the possibility that derangements leading to COPD could arise during lung development. In this review, we summarise the role of leading genes (HHIP, FAM13A, DSP, AGER and TGFB2) identified by genome-wide association studies in lung development and COPD. Because many COPD genome-wide association study genes are enriched in lung epithelial cells, we focus on the role of these genes in the lung epithelium in development, homeostasis and injury.

Precision Approaches to Chronic Obstructive Pulmonary Disease Management

Chronic obstructive pulmonary disease (COPD) is a leading cause of morbidity and mortality worldwide. COPD heterogeneity has hampered progress in developing pharmacotherapies that affect disease progression. This issue can be addressed by precision medicine approaches, which focus on understanding an individual's disease risk, and tailoring management based on pathobiology, environmental exposures, and psychosocial issues. There is an urgent need to identify COPD patients at high risk for poor outcomes and to understand at a mechanistic level why certain individuals are at high risk. Genetics, omics, and network analytic techniques have started to dissect COPD heterogeneity and identify patients with specific pathobiology. Drug repurposing approaches based on biomarkers of specific inflammatory processes (i.e., type 2 inflammation) are promising. As larger data sets, additional omics, and new analytical approaches become available, there will be enormous opportunities to identify high-risk individuals and treat COPD patients based on their specific pathophysiological derangements. These approaches show great promise for risk stratification, early intervention, drug repurposing, and developing novel therapeutic approaches for COPD.

Identification of Potential Differentially-Methylated/Expressed Genes in Chronic Obstructive Pulmonary Disease

Chronic obstructive pulmonary disease (COPD) is a chronic inflammatory lung disease that causes obstructed airflow from the lungs. DNA methylation can regulate gene expression. Understanding the potential molecular mechanism of COPD is of great importance. The aim of this study was to find differentially methylated/expressed genes in COPD. DNA methylation and gene expression profiles in COPD were downloaded from the dataset, followed by functional analysis of differentially-methylated/expressed genes. The potential diagnostic value of these differentially-methylated/expressed genes was determined by receiver operating characteristic (ROC) analysis. Expression validation of differentially-methylated/expressed genes was performed by in vitro experiment and extra online datasets. Totally, 81 hypermethylated-low expression genes and 121 hypomethylated-high expression genes were found in COPD. Among which, 9 core hypermethylated-low expression genes (CD247, CCR7, CD5, IKZF1, SLAMF1, IL2RB, CD3E, CD7 and IL7R) and 8 core hypomethylated-high expression genes (TREM1, AQP9, CD300LF, CLEC12A, NOD2, IRAK3, NLRP3 and LYZ) were identified in the protein-protein interaction (PPI) network. Moreover, these genes had a potential diagnostic utility for COPD. Some signaling pathways were identified in COPD, including T cell receptor signaling pathway, cytokine-cytokine receptor interaction, hematopoietic cell lineage, HTLV-I infection, endocytosis and Jak-STAT signaling pathway. In conclusion, differentially-methylated/expressed genes and involved signaling pathways are likely to be associated with the process of COPD.

De novo identification of complex traits associated with asthma

Introduction

Asthma is a heterogeneous inflammatory disease often associated with other complex phenotypes. Identifying asthma-associated diseases and uncovering the molecular mechanisms mediating their interaction can help detangle the heterogeneity of asthma. Network analysis is a powerful approach for untangling such inter-disease relationships.

Methods

Here, we integrated information on physical contacts between common single nucleotide polymorphisms (SNPs) and gene expression with expression quantitative trait loci (eQTL) data from the lung and whole blood to construct two tissue-specific spatial gene regulatory networks (GRN). We then located the asthma GRN (level 0) within each tissue-specific GRN by identifying the genes that are functionally affected by asthma-associated spatial eQTLs. Curated protein interaction partners were subsequently identified up to four edges or levels away from the asthma GRN. The eQTLs spatially regulating genes on levels 0-4 were queried against the GWAS Catalog to identify the traits enriched (hypergeometric test; FDR ≤ 0.05) in each level.

Results

We identified 80 and 82 traits significantly enriched in the lung and blood GRNs, respectively. All identified traits were previously reported to be comorbid or associated (positively or negatively) with asthma (e.g., depressive symptoms and lung cancer), except 8 traits whose association with asthma is yet to be confirmed (e.g., reticulocyte count). Our analysis additionally pinpoints the variants and genes that link asthma to the identified asthma-associated traits, a subset of which was replicated in a comorbidity analysis using health records of 26,781 asthma patients in New Zealand.

Discussion

Our discovery approach identifies enriched traits in the regulatory space proximal to asthma, in the tissue of interest, without a priori selection of the interacting traits. The predictions it makes expand our understanding of possible shared molecular interactions and therapeutic targets for asthma, where no cure is currently available.

Community Detection in Protein-Protein Interaction Networks and Applications

The ability to identify and characterize not only the protein-protein interactions but also their internal modular organization through network analysis is fundamental for understanding the mechanisms of biological processes at the molecular level. Indeed, the detection of the network communities can enhance our understanding of the molecular basis of disease pathology, and promote drug discovery and disease treatment in personalized medicine. This work gives an overview of recent computational methods for the detection of protein complexes and functional modules in protein-protein interaction networks, also providing a focus on some of its applications. We propose a systematic reformulation of frequently adopted taxonomies for these methods, also proposing new categories to keep up with the most recent research. We review the literature of the last five years (2017-2021) and provide links to existing data and software resources. Finally, we survey recent works exploiting module identification and analysis, in the context of a variety of disease processes for biomarker identification and therapeutic target detection. Our review provides the interested reader with an up-to-date and self-contained view of the existing research, with links to state-of-the-art literature and resources, as well as hints on open issues and future research directions in complex detection and its applications.

Deeper insights into long-term survival heterogeneity of pancreatic ductal adenocarcinoma (PDAC) patients using integrative individual- and group-level transcriptome network analyses

Pancreatic ductal adenocarcinoma (PDAC) is categorized as the leading cause of cancer mortality worldwide. However, its predictive markers for long-term survival are not well known. It is interesting to delineate individual-specific perturbed genes when comparing long-term (LT) and short-term (ST) PDAC survivors and integrate individual- and group-based transcriptome profiling. Using a discovery cohort of 19 PDAC patients from CHU-Liège (Belgium), we first performed differential gene expression analysis comparing LT to ST survivor. Second, we adopted systems biology approaches to obtain clinically relevant gene modules. Third, we created individual-specific perturbation profiles. Furthermore, we used Degree-Aware disease gene prioritizing (DADA) method to develop PDAC disease modules; Network-based Integration of Multi-omics Data (NetICS) to integrate group-based and individual-specific perturbed genes in relation to PDAC LT survival. We identified 173 differentially expressed genes (DEGs) in ST and LT survivors and five modules (including 38 DEGs) showing associations to clinical traits. Validation of DEGs in the molecular lab suggested a role of REG4 and TSPAN8 in PDAC survival. Via NetICS and DADA, we identified various known oncogenes such as CUL1 and TGFB1. Our proposed analytic workflow shows the advantages of combining clinical and omics data as well as individual- and group-level transcriptome profiling.

Genetics of chronic obstructive pulmonary disease: understanding the pathobiology and heterogeneity of a complex disorder

Chronic obstructive pulmonary disease (COPD) is a deadly and highly morbid disease. Susceptibility to and heterogeneity of COPD are incompletely explained by environmental factors such as cigarette smoking. Family-based and population-based studies have shown that a substantial proportion of COPD risk is related to genetic variation. Genetic association studies have identified hundreds of genetic variants that affect risk for COPD, decreased lung function, and other COPD-related traits. These genetic variants are associated with other pulmonary and non-pulmonary traits, demonstrate a genetic basis for at least part of COPD heterogeneity, have a substantial effect on COPD risk in aggregate, implicate early-life events in COPD pathogenesis, and often involve genes not previously suspected to have a role in COPD. Additional progress will require larger genetic studies with more ancestral diversity, improved profiling of rare variants, and better statistical methods. Through integration of genetic data with other omics data and comprehensive COPD phenotypes, as well as functional description of causal mechanisms for genetic risk variants, COPD genetics will continue to inform novel approaches to understanding the pathobiology of COPD and developing new strategies for management and treatment.

Protein interaction networks provide insight into fetal origins of chronic obstructive pulmonary disease

BACKGROUND

Chronic obstructive pulmonary disease (COPD) is a leading cause of death in adults that may have origins in early lung development. It is a complex disease, influenced by multiple factors including genetic variants and environmental factors. Maternal smoking during pregnancy may influence the risk for diseases during adulthood, potentially through epigenetic modifications including methylation.

METHODS

In this work, we explore the fetal origins of COPD by utilizing lung DNA methylation marks associated with in utero smoke (IUS) exposure, and evaluate the network relationships between methylomic and transcriptomic signatures associated with adult lung tissue from former smokers with and without COPD. To identify potential pathobiological mechanisms that may link fetal lung, smoke exposure and adult lung disease, we study the interactions (physical and functional) of identified genes using protein-protein interaction networks.

RESULTS

We build IUS-exposure and COPD modules, which identify connected subnetworks linking fetal lung smoke exposure to adult COPD. Studying the relationships and connectivity among the different modules for fetal smoke exposure and adult COPD, we identify enriched pathways, including the AGE-RAGE and focal adhesion pathways.

CONCLUSIONS

The modules identified in our analysis add new and potentially important insights to understanding the early life molecular perturbations related to the pathogenesis of COPD. We identify AGE-RAGE and focal adhesion as two biologically plausible pathways that may reveal lung developmental contributions to COPD. We were not only able to identify meaningful modules but were also able to study interconnections between smoke exposure and lung disease, augmenting our knowledge about the fetal origins of COPD.

Identification of Dysregulated Mechanisms and Candidate Gene Markers in Chronic Obstructive Pulmonary Disease

Purpose

This study aimed to identify candidate gene markers that may facilitate chronic obstructive pulmonary disease (COPD) diagnosis and treatment.

Methods

The GSE47460 and GSE151052 datasets were analyzed to identify differentially expressed mRNAs (DEmRs) between COPD patients and controls. DEmRs that were differentially expressed in the same direction in both datasets were analyzed for functional enrichment and for coexpression. Genes from the largest three modules were tested for their ability to diagnose COPD based on the area under the receiver operating characteristic curve (AUC). Genes with AUC > 0.7 in both datasets were used to perform regression based on the "least absolute shrinkage and selection operator" in order to identify feature genes. We also identified differentially expressed miRNAs (DEmiRs) between COPD patients and controls using the GSE38974 dataset, then constructed a regulatory network. We also examined associations between feature genes and immune cell infiltration in COPD, and we identified methylation markers of COPD using the GSE63704 dataset.

Results

A total of 1350 genes differentially regulated in the same direction in the GSE47460 and GSE151052 datasets were found. The genes were significantly enriched in immune-related biological functions. Of 186 modules identified using MEGENA, the largest were C1_ 6, C1_ 3, and C1_ 2. Of the 22 candidate genes screened based on AUC, 11 feature genes emerged from analysis of a subset of GSE47460 data, which we validated using another subset of GSE47460 data as well as the independent GSE151052 dataset. Feature genes correlated significantly with infiltration by immune cells. The feature genes GPC4 and RS1 were predicted to be regulated by miR-374a-3p. We identified 117 candidate methylation markers of COPD, including PRRG4.

Conclusion

The feature genes we identified may be potential diagnostic markers and therapeutic targets in COPD. These findings provide new leads for exploring disease mechanisms and targeted treatments.

Integrating molecular interactions and gene expression to identify biomarkers and network modules of chronic obstructive pulmonary disease

BACKGROUND:

Chronic obstructive pulmonary disease (COPD) causes chronic obstructive conditions, chronic bronchitis, and emphysema, and is a major cause of death worldwide. Although several efforts for identifying biomarkers and pathways have been made, specific causal COPD mechanism remains unknown.

OBJECTIVE:

This study combined biological interaction data with gene expression data for a better understanding of the biological process and network module for COPD.

METHODS:

Using a sparse network-based method, we selected 49 genes from peripheral blood mononuclear cell expression data of 136 subjects, including 42 ex-smoking controls and 94 subjects with COPD.

RESULTS:

These 49 genes might influence biological processes and molecular functions related to COPD. For example, our result suggests that FoxO signaling may contribute to the atrophy of COPD peripheral muscle tissues via oxidative stress.

CONCLUSIONS:

Our approach enhances the existing understanding of COPD disease pathogenesis and predicts new genetic markers and pathways that may influence COPD pathogenesis.

Integrating molecular interactions and gene expression to identify biomarkers to predict response to tumor necrosis factor inhibitor therapies in rheumatoid arthritis patients

BACKGROUND

Targeted therapy using anti-TNF (tumor necrosis factor) is the first option for patients with rheumatoid arthritis (RA). Anti-TNF therapy, however, does not lead to meaningful clinical improvement in many RA patients. To predict which patients will not benefit from anti-TNF therapy, clinical tests should be performed prior to treatment beginning.

OBJECTIVE

Although various efforts have been made to identify biomarkers and pathways that may be helpful to predict the response to anti-TNF treatment, gaps remain in clinical use due to the low predictive power of the selected biomarkers.

METHODS

In this paper, we used a network-based computational method to identify the select the predictive biomarkers to guide the treatment of RA patients.

RESULTS

We select 69 genes from peripheral blood expression data from 46 subjects using a sparse network-based method. The result shows that the selected 69 genes might influence biological processes and molecular functions related to the treatment.

CONCLUSIONS

Our approach advances the predictive power of anti-TNF therapy response and provides new genetic markers and pathways that may influence the treatment.

Gentamicin-induced hearing loss: A retrospective study using the Food and Drug Administration Adverse Event Reporting System and a toxicological study using drug-gene network analysis

The objectives of the study were to evaluate the relationship between gentamicin (GEN) and hearing loss using the Food and Drug Administration Adverse Event Reporting system (FAERS) database and elucidate the potential toxicological mechanism of GEN-induced hearing loss through a drug–gene network analysis.

Using the preferred terms and standardized queries from the Medical Dictionary for Regulatory Activities, we calculated the reporting odds ratios (RORs). We extracted GEN-associated genes (seed genes) and analyzed drug−gene interactions using the ClueGO plug-in in the Cytoscape software and the DIseAse MOdule Detection (DIAMOnD) algorithm.

The lower limit of the 95% confidence interval (CI) of the ROR for aminoglycosides (AG) antibacterials was over 1, and the ROR was 5.5 (5.1–6.0). We retrieved 17 seed genes related to GEN from the PharmGKB and Drug Gene Interaction databases. In total, 1018 human genes interacting with GEN were investigated using ClueGO. Through Molecular Complex Detection (MCODE) analysis, we identified 17 local gene clusters. The nodes and edges of the highest-ranked local gene cluster named “Cluster 1” were 30 and 433, respectively. According to the ClueGO analysis using the Kyoto Encyclopedia of Genes and Genomes (KEGG), Cluster 1 genes were highly enriched in “oxidative phosphorylation.” According to the ClueGO analysis using ClinVar, Cluster 1 genes were highly enriched in “mitochondrial diseases,” “mitochondrial complex I deficiency,” “hereditary hearing loss and deafness,” and “Leigh syndrome.” We identified 60 GEN-associated genes using the DIAMOnD algorithm. Several GEN-associated genes in the DIAMOnD algorithm were highly enriched in “PI3K-Akt signaling pathway,” “Ras signaling pathway,” “focal adhesion,” “MAPK signaling pathway,” “regulation of actin cytoskeleton,” “oxidative phosphorylation,” and “ECM-receptor interaction.”

Our analysis demonstrated an association between several AGs and hearing loss using the FAERS database. Drug−gene network analysis demonstrated that GEN may be associated with oxidative phosphorylation-associated genes and integrin genes, which may be associated with hearing loss.

Connecting COPD GWAS genes: FAM13A controls TGFβ2 secretion by modulating AP-3 transport

Chronic Obstructive Pulmonary Disease (COPD) is a common, complex disease and a major cause of morbidity and mortality. Although multiple genetic determinants of COPD have been implicated by genome-wide association studies (GWAS), the pathophysiologic significance of these associations remains largely unknown. From a COPD protein-protein interaction network module, we selected a network path between two COPD GWAS genes for validation studies: FAM13A-AP3D1-CTGF-TGFB2. We find that TGFβ2, FAM13A, and AP3D1 (but not CTGF) form a cellular protein complex. Functional characterization suggests that this complex mediates the secretion of TGFβ2 through an AP-3-dependent pathway, with FAM13A acting as a negative regulator by targeting a late stage of this transport that involves the dissociation of coat-cargo interaction. Moreover, we find that TGFβ2 is a transmembrane protein that engages the AP-3 complex for delivery to the late endosomal compartments for subsequent secretion through exosomes. These results identify a pathophysiologic context that unifies the biological network role of two COPD GWAS proteins and reveal novel mechanisms of cargo transport through an intracellular pathway.

Epigenetics and pulmonary diseases in the horizon of precision medicine: a review

Epigenetic mechanisms represent potential molecular routes which could bridge the gap between genetic background and environmental risk factors contributing to the pathogenesis of pulmonary diseases. In patients with COPD, asthma and pulmonary arterial hypertension (PAH), there is emerging evidence of aberrant epigenetic marks, mainly including DNA methylation and histone modifications which directly mediate reversible modifications to the DNA without affecting the genomic sequence. Post-translational events and microRNAs can be also regulated epigenetically and potentially participate in disease pathogenesis. Thus, novel pathogenic mechanisms and putative biomarkers may be detectable in peripheral blood, sputum, nasal and buccal swabs or lung tissue. Besides, DNA methylation plays an important role during the early phases of fetal development and may be impacted by environmental exposures, ultimately influencing an individual's susceptibility to COPD, asthma and PAH later in life. With the advances in omics platforms and the application of computational biology tools, modelling the epigenetic variability in a network framework, rather than as single molecular defects, provides insights into the possible molecular pathways underlying the pathogenesis of COPD, asthma and PAH. Epigenetic modifications may have clinical applications as noninvasive biomarkers of pulmonary diseases. Moreover, combining molecular assays with network analysis of epigenomic data may aid in clarifying the multistage transition from a "pre-disease" to "disease" state, with the goal of improving primary prevention of lung diseases and its subsequent clinical management.We describe epigenetic mechanisms known to be associated with pulmonary diseases and discuss how network analysis could improve our understanding of lung diseases.

Exploration of the sputum methylome and omics deconvolution by quadratic programming in molecular profiling of asthma and COPD: the road to sputum omics 2.0

BACKGROUND

To date, most studies involving high-throughput analyses of sputum in asthma and COPD have focused on identifying transcriptomic signatures of disease. No whole-genome methylation analysis of sputum cells has been performed yet. In this context, the highly variable cellular composition of sputum has potential to confound the molecular analyses.

METHODS

Whole-genome transcription (Agilent Human 4 × 44 k array) and methylation (Illumina 450 k BeadChip) analyses were performed on sputum samples of 9 asthmatics, 10 healthy and 10 COPD subjects. RNA integrity was checked by capillary electrophoresis and used to correct in silico for bias conferred by RNA degradation during biobank sample storage. Estimates of cell type-specific molecular profiles were derived via regression by quadratic programming based on sputum differential cell counts. All analyses were conducted using the open-source R/Bioconductor software framework.

RESULTS

A linear regression step was found to perform well in removing RNA degradation-related bias among the main principal components of the gene expression data, increasing the number of genes detectable as differentially expressed in asthma and COPD sputa (compared to controls). We observed a strong influence of the cellular composition on the results of mixed-cell sputum analyses. Exemplarily, upregulated genes derived from mixed-cell data in asthma were dominated by genes predominantly expressed in eosinophils after deconvolution. The deconvolution, however, allowed to perform differential expression and methylation analyses on the level of individual cell types and, though we only analyzed a limited number of biological replicates, was found to provide good estimates compared to previously published data about gene expression in lung eosinophils in asthma. Analysis of the sputum methylome indicated presence of differential methylation in genomic regions of interest, e.g. mapping to a number of human leukocyte antigen (HLA) genes related to both major histocompatibility complex (MHC) class I and II molecules in asthma and COPD macrophages. Furthermore, we found the SMAD3 (SMAD family member 3) gene, among others, to lie within differentially methylated regions which has been previously reported in the context of asthma.

CONCLUSIONS

In this methodology-oriented study, we show that methylation profiling can be easily integrated into sputum analysis workflows and exhibits a strong potential to contribute to the profiling and understanding of pulmonary inflammation. Wherever RNA degradation is of concern, in silico correction can be effective in improving both sensitivity and specificity of downstream analyses. We suggest that deconvolution methods should be integrated in sputum omics analysis workflows whenever possible in order to facilitate the unbiased discovery and interpretation of molecular patterns of inflammation.

Long noncoding RNA MCM3AP antisense RNA 1 is downregulated in chronic obstructive pulmonary disease and regulates human bronchial smooth muscle cell proliferation

OBJECTIVES

This study aimed to investigate the involvement of MCM3AP antisense RNA 1 (MCM3AP-AS1) in chronic obstructive pulmonary disease (COPD).

METHODS

The expression levels of plasma MCM3AP-AS1 in COPD patients and healthy controls were measured by quantitative PCR before treatment and at 3 months after the initiation of treatment (post-treatment) from COPD patients. The role of MCM3AP-AS1 in regulating the proliferation of human bronchial smooth muscle cells (HBSMCs) was explored by a cell proliferation assay.

RESULTS

We found that MCM3AP-AS1 expression was downregulated in the plasma of COPD patients compared with controls. Among controls, MCM3AP-AS1 expression was lower in smokers than never-smokers. A 3-year follow-up study showed that, among smokers, patients with low MCM3AP-AS1 expression showed a higher incidence of COPD. After treatment for COPD, MCM3AP-AS1 expression significantly increased. The cell proliferation assay showed that MCM3AP-AS1 overexpression decreased the proliferation rate of HBSMCs. MCM3AP-AS1 silencing had the opposite effect.

CONCLUSIONS

MCM3AP-AS1 appears to be downregulated in COPD and to predict its occurrence. MCM3AP-AS1 regulates the proliferation of HBSMCs to participate in airway remodeling.

The application of big data to cardiovascular disease: paths to precision medicine

Advanced phenotyping of cardiovascular diseases has evolved with the application of high-resolution omics screening to populations enrolled in large-scale observational and clinical trials. This strategy has revealed that considerable heterogeneity exists at the genotype, endophenotype, and clinical phenotype levels in cardiovascular diseases, a feature of the most common diseases that has not been elucidated by conventional reductionism. In this discussion, we address genomic context and (endo)phenotypic heterogeneity, and examine commonly encountered cardiovascular diseases to illustrate the genotypic underpinnings of (endo)phenotypic diversity. We highlight the existing challenges in cardiovascular disease genotyping and phenotyping that can be addressed by the integration of big data and interpreted using novel analytical methodologies (network analysis). Precision cardiovascular medicine will only be broadly applied to cardiovascular patients once this comprehensive data set is subjected to unique, integrative analytical strategies that accommodate molecular and clinical heterogeneity rather than ignore or reduce it.

SNPeffect: identifying functional roles of SNPs using metabolic networks

Genetic sources of phenotypic variation have been a focus of plant studies aimed at improving agricultural yield and understanding adaptive processes. Genome-wide association studies identify the genetic background behind a trait by examining associations between phenotypes and single-nucleotide polymorphisms (SNPs). Although such studies are common, biological interpretation of the results remains a challenge; especially due to the confounding nature of population structure and the systematic biases thus introduced. Here, we propose a complementary analysis (SNPeffect) that offers putative genotype-to-phenotype mechanistic interpretations by integrating biochemical knowledge encoded in metabolic models. SNPeffect is used to explain differential growth rate and metabolite accumulation in A. thaliana and P. trichocarpa accessions as the outcome of SNPs in enzyme-coding genes. To this end, we also constructed a genome-scale metabolic model for Populus trichocarpa, the first for a perennial woody tree. As expected, our results indicate that growth is a complex polygenic trait governed by carbon and energy partitioning. The predicted set of functional SNPs in both species are associated with experimentally characterized growth-determining genes and also suggest putative ones. Functional SNPs were found in pathways such as amino acid metabolism, nucleotide biosynthesis, and cellulose and lignin biosynthesis, in line with breeding strategies that target pathways governing carbon and energy partition.

Exploring the cross-phenotype network region of disease modules reveals concordant and discordant pathways between chronic obstructive pulmonary disease and idiopathic pulmonary fibrosis

Chronic obstructive pulmonary disease (COPD) and idiopathic pulmonary fibrosis (IPF) are two pathologically distinct chronic lung diseases that are associated with cigarette smoking. Genetic studies have identified shared loci for COPD and IPF, including several loci with opposite directions of effect. The existence of additional shared genetic loci, as well as potential shared pathobiological mechanisms between the two diseases at the molecular level, remains to be explored. Taking a network-based approach, we built disease modules for COPD and IPF using genome-wide association studies-implicated genes. The two disease modules displayed strong disease signals in an independent gene expression data set of COPD and IPF lung tissue and showed statistically significant overlap and network proximity, sharing 19 genes, including ARHGAP12 and BCHE. To uncover pathways at the intersection of COPD and IPF, we developed a metric, NetPathScore, which prioritizes the pathways of a disease by their network overlap with another disease. Applying NetPathScore to the COPD and IPF disease modules enabled the determination of concordant and discordant pathways between these diseases. Concordant pathways between COPD and IPF included extracellular matrix remodeling, Mitogen-activated protein kinase (MAPK) signaling and ALK pathways, whereas discordant pathways included advanced glycosylation end product receptor signaling and telomere maintenance and extension pathways. Overall, our findings reveal shared molecular interaction regions between COPD and IPF and shed light on the congruent and incongruent biological processes lying at the intersection of these two complex diseases.

Integrated transcriptomic correlation network analysis identifies COPD molecular determinants

Chronic obstructive pulmonary disease (COPD) is a complex and heterogeneous syndrome. Network-based analysis implemented by SWIM software can be exploited to identify key molecular switches - called "switch genes" - for the disease. Genes contributing to common biological processes or defining given cell types are usually co-regulated and co-expressed, forming expression network modules. Consistently, we found that the COPD correlation network built by SWIM consists of three well-characterized modules: one populated by switch genes, all up-regulated in COPD cases and related to the regulation of immune response, inflammatory response, and hypoxia (like TIMP1, HIF1A, SYK, LY96, BLNK and PRDX4); one populated by well-recognized immune signature genes, all up-regulated in COPD cases; one where the GWAS genes AGER and CAVIN1 are the most representative module genes, both down-regulated in COPD cases. Interestingly, 70% of AGER negative interactors are switch genes including PRDX4, whose activation strongly correlates with the activation of known COPD GWAS interactors SERPINE2, CD79A, and POUF2AF1. These results suggest that SWIM analysis can identify key network modules related to complex diseases like COPD.

Discovering the genes mediating the interactions between chronic respiratory diseases in the human interactome

The molecular and clinical features of a complex disease can be influenced by other diseases affecting the same individual. Understanding disease-disease interactions is therefore crucial for revealing shared molecular mechanisms among diseases and designing effective treatments. Here we introduce Flow Centrality (FC), a network-based approach to identify the genes mediating the interaction between two diseases in a protein-protein interaction network. We focus on asthma and COPD, two chronic respiratory diseases that have been long hypothesized to share common genetic determinants and mechanisms. We show that FC highlights potential mediator genes between the two diseases, and observe similar outcomes when applying FC to 66 additional pairs of related diseases. Further, we perform in vitro perturbation experiments on a widely replicated asthma gene, GSDMB, showing that FC identifies candidate mediators of the interactions between GSDMB and COPD-associated genes. Our results indicate that FC predicts promising gene candidates for further study of disease-disease interactions.

Recent advances in network medicine: From disease mechanisms to new treatment strategies

Conventional reductionist approaches have guided most of our understanding in disease diagnostic and treatment. However, most diseases are not consequence of perturbations in a single protein or metabolite, but rather of the effect that these perturbations have in their cellular context. The emerging field of network medicine offers a set of tools to explore molecular networks and to retrieve insights about mechanisms of different diseases. The study of the protein interactome, the map of physical interactions among human proteins, revealed that disease proteins tend to interact with each other, linking diseases to well-defined interactome neighborhoods. These disease-associated neighborhoods have been defined as disease modules, and they can uncover the biological significance of genes identified by genetic studies, reveal molecular mechanisms that connect different phenotypes, and help identify new pharmacological strategies for disease treatment. Therefore, network medicine offers a framework in which the complexity of different aspects of multiple sclerosis can be explored in an integrative fashion, which can ultimately provide insights about disease mechanisms and treatment.

Genetics of COPD

Although chronic obstructive pulmonary disease (COPD) risk is strongly influenced by cigarette smoking, genetic factors are also important determinants of COPD. In addition to Mendelian syndromes such as alpha-1 antitrypsin deficiency, many genomic regions that influence COPD susceptibility have been identified in genome-wide association studies. Similarly, multiple genomic regions associated with COPD-related phenotypes, such as quantitative emphysema measures, have been found. Identifying the functional variants and key genes within these association regions remains a major challenge. However, newly identified COPD susceptibility genes are already providing novel insights into COPD pathogenesis. Network-based approaches that leverage these genetic discoveries have the potential to assist in decoding the complex genetic architecture of COPD.

Network Diffusion Approach to Predict LncRNA Disease Associations Using Multi-Type Biological Networks: LION

Recently, long-non-coding RNAs (lncRNAs) have attracted attention because of their emerging role in many important biological mechanisms. The accumulating evidence indicates that the dysregulation of lncRNAs is associated with complex diseases. However, only a few lncRNA-disease associations have been experimentally validated and therefore, predicting potential lncRNAs that are associated with diseases become an important task. Current computational approaches often use known lncRNA-disease associations to predict potential lncRNA-disease links. In this work, we exploited the topology of multi-level networks to propose the LncRNA rankIng by NetwOrk DiffusioN (LION) approach to identify lncRNA-disease associations. The multi-level complex network consisted of lncRNA-protein, protein–protein interactions, and protein-disease associations. We applied the network diffusion algorithm of LION to predict the lncRNA-disease associations within the multi-level network. LION achieved an AUC value of 96.8% for cardiovascular diseases, 91.9% for cancer, and 90.2% for neurological diseases by using experimentally verified lncRNAs associated with diseases. Furthermore, compared to a similar approach (TPGLDA), LION performed better for cardiovascular diseases and cancer. Given the versatile role played by lncRNAs in different biological mechanisms that are perturbed in diseases, LION’s accurate prediction of lncRNA-disease associations helps in ranking lncRNAs that could function as potential biomarkers and potential drug targets.

Network Medicine in the Age of Biomedical Big Data

Network medicine is an emerging area of research dealing with molecular and genetic interactions, network biomarkers of disease, and therapeutic target discovery. Large-scale biomedical data generation offers a unique opportunity to assess the effect and impact of cellular heterogeneity and environmental perturbations on the observed phenotype. Marrying the two, network medicine with biomedical data provides a framework to build meaningful models and extract impactful results at a network level. In this review, we survey existing network types and biomedical data sources. More importantly, we delve into ways in which the network medicine approach, aided by phenotype-specific biomedical data, can be gainfully applied. We provide three paradigms, mainly dealing with three major biological network archetypes: protein-protein interaction, expression-based, and gene regulatory networks. For each of these paradigms, we discuss a broad overview of philosophies under which various network methods work. We also provide a few examples in each paradigm as a test case of its successful application. Finally, we delineate several opportunities and challenges in the field of network medicine. We hope this review provides a lexicon for researchers from biological sciences and network theory to come on the same page to work on research areas that require interdisciplinary expertise. Taken together, the understanding gained from combining biomedical data with networks can be useful for characterizing disease etiologies and identifying therapeutic targets, which, in turn, will lead to better preventive medicine with translational impact on personalized healthcare.