Search

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.

Abstract Background Historically, the molecular classification of colorectal cancer (CRC) was based on the global genomic status, which identified microsatellite instability in mismatch repair (MMR) deficient CRC, and chromosomal instability in MMR proficient CRC. With the introduction of immune checkpoint inhibitors, the microsatellite and chromosomal instability classification regained momentum as the microsatellite instability condition predicted sensitivity to immune checkpoint inhibitors, possibly due to both high tumor mutation burden (TMB) and high levels of infiltrating lymphocytes. Conversely, proficient MMR CRC are mostly resistant to immunotherapy. To better understand the relationship between the microsatellite and chromosomal instability classification, and eventually discover additional CRC subgroups relevant for therapeutic decisions, we developed a computational pipeline that include molecular integrative analysis of genomic, epigenomic and transcriptomic data. Results The first step of the pipeline was based on unsupervised hierarchical clustering analysis of copy number variations (CNVs) versus hypermutation status that identified a first CRC cluster with few CNVs enriched in Hypermutated and microsatellite instability samples, a second CRC cluster with a high number of CNVs mostly including non-HM and microsatellite stable samples, and a third cluster (7.8% of the entire dataset) with low CNVs and low TMB, which shared clinical-pathological features with Hypermutated CRCs and thus defined Hypermutated-like CRCs. The mutational features, DNA methylation profile and base substitution fingerprints of these tumors revealed that Hypermutated-like patients are molecularly distinct from Hypermutated and non-Hypermutated tumors and are likely to develop and progress through different genetic events. Transcriptomic analysis highlighted further differences amongst the three groups and revealed an inflamed tumor microenvironment and modulation Immune Checkpoint Genes in Hypermutated-like CRCs. Conclusion Therefore, our work highlights Hypermutated-like tumors as a distinct and previously unidentified CRC subgroup possibly responsive to immune checkpoint inhibitors. If further validated, these findings can lead to expanding the fraction of patients eligible to immunotherapy.

Background Triple-negative breast cancers (TNBCs) display poor prognosis, have a high risk of tumour recurrence, and exhibit high resistance to drug treatments. Based on their gene expression profiles, the majority of TNBCs are classified as basal-like breast cancers. Currently, there are not available widely-accepted prognostic markers to predict outcomes in basal-like subtype, so the selection of new prognostic indicators for this BC phenotype represents an unmet clinical challenge. Results Here, we attempted to address this challenging issue by exploiting a bioinformatics pipeline able to integrate transcriptomic, genomic, epigenomic, and clinical data freely accessible from public repositories. This pipeline starts from the application of the well-established network-based SWIM methodology on the transcriptomic data to unveil important (switch) genes in relation with a complex disease of interest. Then, survival and linear regression analyses are performed to associate the gene expression profiles of the switch genes with both the patients’ clinical outcome and the disease aggressiveness. This allows us to identify a prognostic gene signature that in turn is fed to the last step of the pipeline consisting of an analysis at DNA level, to investigate whether variations in the expression of identified prognostic switch genes could be related to genetic (copy number variations) or epigenetic (DNA methylation differences) alterations in their gene loci, or to the activities of transcription factors binding to their promoter regions. Finally, changes in the protein expression levels corresponding to the so far identified prognostic switch genes are evaluated by immunohistochemical staining results taking advantage of the Human Protein Atlas. Conclusion The application of the proposed pipeline on the dataset of The Cancer Genome Atlas (TCGA)-Breast Invasive Carcinoma (BRCA) patients affected by basal-like subtype led to an in silico recognition of a basal-like specific gene signature composed of 11 potential prognostic biomarkers to be further investigated.

Breast cancer (BC) is a heterogeneous and complex disease characterized by different subtypes with distinct morphologies and clinical implications and for which new and effective treatment options are urgently demanded. The computational approaches recently developed for drug repurposing provide a very promising opportunity to offer tools that efficiently screen potential novel medical indications for various drugs that are already approved and used in clinical practice. Here, we started with disease-associated genes that were identified through a transcriptome-based analysis, which we used to predict potential repurposable drugs for various breast cancer subtypes by using an algorithm that we developed for drug repurposing called SAveRUNNER. Our findings were also in silico validated by performing a gene set enrichment analysis, which confirmed that most of the predicted repurposable drugs may have a potential treatment effect against breast cancer pathophenotypes.

Extensive molecular characterization of human colorectal cancer (CRC) via Next Generation Sequencing (NGS) indicated that genetic or epigenetic dysregulation of a relevant, but limited, number of molecular pathways typically occurs in this tumor. The molecular picture of the disease is significantly complicated by the frequent occurrence of individually rare genetic aberrations, which expand tumor heterogeneity. Inter- and intratumor molecular heterogeneity is very likely responsible for the remarkable individual variability in the response to conventional and target-driven first-line therapies, in metastatic CRC (mCRC) patients, whose median overall survival remains unsatisfactory. Implementation of an extensive molecular characterization of mCRC in the clinical routine does not yet appear feasible on a large scale, while multigene panel sequencing of most commonly mutated oncogene/oncosuppressor hotspots is more easily achievable. Here, we report that clinical multigene panel sequencing performed for anti-EGFR therapy predictive purposes in 639 formalin-fixed paraffin-embedded (FFPE) mCRC specimens revealed previously unknown pairwise mutation associations and a high proportion of cases carrying actionable gene mutations. Most importantly, a simple principal component analysis directed the delineation of a new molecular stratification of mCRC patients in eight groups characterized by non-random, specific mutational association patterns (MAPs), aggregating samples with similar biology. These data were validated on a The Cancer Genome Atlas (TCGA) CRC dataset. The proposed stratification may provide great opportunities to direct more informed therapeutic decisions in the majority of mCRC cases.