Metabolomics for early pancreatic cancer detection in plasma samples from a Swedish prospective population-based biobank.

Background: Pancreatic ductal adenocarcinoma (pancreatic cancer) is often detected at late stages resulting in poor overall survival. To improve survival, more patients need to be diagnosed early when curative surgery is feasible. We aimed to identify circulating metabolites that could be used as early pancreatic cancer biomarkers. Methods: We performed metabolomics by liquid and gas chromatography-mass spectrometry in plasma samples from 82 future pancreatic cancer patients and 82 matched healthy controls within the Northern Sweden Health and Disease Study (NSHDS). Logistic regression was used to assess univariate associations between metabolites and pancreatic cancer risk. Least absolute shrinkage and selection operator (LASSO) logistic regression was used to design a metabolite-based risk score. We used receiver operating characteristic (ROC) analyses to assess the discriminative performance of the metabolite-based risk score. Results: Among twelve risk-associated metabolites with a nominal P value <0.05, we defined a risk score of three metabolites [indoleacetate, 3-hydroxydecanoate (10:0-OH), and retention index (RI): 2,745.4] using LASSO. A logistic regression model containing these three metabolites, age, sex, body mass index (BMI), smoking status, sample date, fasting status, and carbohydrate antigen 19-9 (CA 19-9) yielded an internal area under curve (AUC) of 0.784 [95% confidence interval (CI): 0.714-0.854] compared to 0.681 (95% CI: 0.597-0.764) for a model without these metabolites (P value =0.007). Seventeen metabolites were significantly associated with pancreatic cancer survival [false discovery rate (FDR) <0.1]. Conclusions: Indoleacetate, 3-hydroxydecanoate (10:0-OH), and RI: 2,745.4 were identified as the top candidate biomarkers for early detection. However, continued efforts are warranted to determine the usefulness of these metabolites as early pancreatic cancer biomarkers.

Proteomics reveal biomarkers for diagnosis, disease activity and long-term disability outcomes in multiple sclerosis.

Sensitive and reliable protein biomarkers are needed to predict disease trajectory and personalize treatment strategies for multiple sclerosis (MS). Here, we use the highly sensitive proximity-extension assay combined with next-generation sequencing (Olink Explore) to quantify 1463 proteins in cerebrospinal fluid (CSF) and plasma from 143 people with early-stage MS and 43 healthy controls. With longitudinally followed discovery and replication cohorts, we identify CSF proteins that consistently predicted both short- and long-term disease progression. Lower levels of neurofilament light chain (NfL) in CSF is superior in predicting the absence of disease activity two years after sampling (replication AUC = 0.77) compared to all other tested proteins. Importantly, we also identify a combination of 11 CSF proteins (CXCL13, LTA, FCN2, ICAM3, LY9, SLAMF7, TYMP, CHI3L1, FYB1, TNFRSF1B and NfL) that predict the severity of disability worsening according to the normalized age-related MS severity score (replication AUC = 0.90). The identification of these proteins may help elucidate pathogenetic processes and might aid decisions on treatment strategies for persons with MS.

Bioinformatics analysis of miRNAs in the neuroblastoma 11q-deleted region reveals a role of miR-548l in both 11q-deleted and MYCN amplified tumour cells.

Neuroblastoma is a childhood tumour that is responsible for approximately 15% of all childhood cancer deaths. Neuroblastoma tumours with amplification of the oncogene MYCN are aggressive, however, another aggressive subgroup without MYCN amplification also exists; rather, they have a deleted region at chromosome arm 11q. Twenty-six miRNAs are located within the breakpoint region of chromosome 11q and have been checked for a possible involvement in development of neuroblastoma due to the genomic alteration. Target genes of these miRNAs are involved in pathways associated with cancer, including proliferation, apoptosis and DNA repair. We could show that miR-548l found within the 11q region is downregulated in neuroblastoma cell lines with 11q deletion or MYCN amplification. In addition, we showed that the restoration of miR-548l level in a neuroblastoma cell line led to a decreased proliferation of these cells as well as a decrease in the percentage of cells in the S phase. We also found that miR-548l overexpression suppressed cell viability and promoted apoptosis, while miR-548l knockdown promoted cell viability and inhibited apoptosis in neuroblastoma cells. Our results indicate that 11q-deleted neuroblastoma and MYCN amplified neuroblastoma coalesce by downregulating miR-548l.

MODalyseR-a novel software for inference of disease module hub regulators identified a putative multiple sclerosis regulator supported by independent eQTL data.

Motivation: Network-based disease modules have proven to be a powerful concept for extracting knowledge about disease mechanisms, predicting for example disease risk factors and side effects of treatments. Plenty of tools exist for the purpose of module inference, but less effort has been put on simultaneously utilizing knowledge about regulatory mechanisms for predicting disease module hub regulators. Results: We developed MODalyseR, a novel software for identifying disease module regulators and reducing modules to the most disease-associated genes. This pipeline integrates and extends previously published software packages MODifieR and ComHub and hereby provides a user-friendly network medicine framework combining the concepts of disease modules and hub regulators for precise disease gene identification from transcriptomics data. To demonstrate the usability of the tool, we designed a case study for multiple sclerosis that revealed IKZF1 as a promising hub regulator, which was supported by independent ChIP-seq data. Availability and implementation: MODalyseR is available as a Docker image at https://hub.docker.com/r/ddeweerd/modalyser with user guide and installation instructions found at https://gustafsson-lab.gitlab.io/MODalyseR/. Supplementary information: Supplementary data are available at Bioinformatics Advances online.

TFTenricher: a python toolbox for annotation enrichment analysis of transcription factor target genes.

BACKGROUND: Transcription factors (TFs) are the upstream regulators that orchestrate gene expression, and therefore a centrepiece in bioinformatics studies. While a core strategy to understand the biological context of genes and proteins includes annotation enrichment analysis, such as Gene Ontology term enrichment, these methods are not well suited for analysing groups of TFs. This is particularly true since such methods do not aim to include downstream processes, and given a set of TFs, the expected top ontologies would revolve around transcription processes. RESULTS: We present the TFTenricher, a Python toolbox that focuses specifically at identifying gene ontology terms, cellular pathways, and diseases that are over-represented among genes downstream of user-defined sets of human TFs. We evaluated the inference of downstream gene targets with respect to false positive annotations, and found an inference based on co-expression to best predict downstream processes. Based on these downstream genes, the TFTenricher uses some of the most common databases for gene functionalities, including GO, KEGG and Reactome, to calculate functional enrichments. By applying the TFTenricher to differential expression of TFs in 21 diseases, we found significant terms associated with disease mechanism, while the gene set enrichment analysis on the same dataset predominantly identified processes related to transcription. CONCLUSIONS AND AVAILABILITY: The TFTenricher package enables users to search for biological context in any set of TFs and their downstream genes. The TFTenricher is available as a Python 3 toolbox at https://github.com/rasma774/Tftenricher , under a GNU GPL license and with minimal dependencies.

A validated generally applicable approach using the systematic assessment of disease modules by GWAS reveals a multi-omic module strongly associated with risk factors in multiple sclerosis.

BACKGROUND: There exist few, if any, practical guidelines for predictive and falsifiable multi-omic data integration that systematically integrate existing knowledge. Disease modules are popular concepts for interpreting genome-wide studies in medicine but have so far not been systematically evaluated and may lead to corroborating multi-omic modules. RESULT: We assessed eight module identification methods in 57 previously published expression and methylation studies of 19 diseases using GWAS enrichment analysis. Next, we applied the same strategy for multi-omic integration of 20 datasets of multiple sclerosis (MS), and further validated the resulting module using both GWAS and risk-factor-associated genes from several independent cohorts. Our benchmark of modules showed that in immune-associated diseases modules inferred from clique-based methods were the most enriched for GWAS genes. The multi-omic case study using MS data revealed the robust identification of a module of 220 genes. Strikingly, most genes of the module were differentially methylated upon the action of one or several environmental risk factors in MS (n = 217, P = 10- 47) and were also independently validated for association with five different risk factors of MS, which further stressed the high genetic and epigenetic relevance of the module for MS. CONCLUSIONS: We believe our analysis provides a workflow for selecting modules and our benchmark study may help further improvement of disease module methods. Moreover, we also stress that our methodology is generally applicable for combining and assessing the performance of multi-omic approaches for complex diseases.

ComHub: Community predictions of hubs in gene regulatory networks.

BACKGROUND: Hub transcription factors, regulating many target genes in gene regulatory networks (GRNs), play important roles as disease regulators and potential drug targets. However, while numerous methods have been developed to predict individual regulator-gene interactions from gene expression data, few methods focus on inferring these hubs. RESULTS: We have developed ComHub, a tool to predict hubs in GRNs. ComHub makes a community prediction of hubs by averaging over predictions by a compendium of network inference methods. Benchmarking ComHub against the DREAM5 challenge data and two independent gene expression datasets showed a robust performance of ComHub over all datasets. CONCLUSIONS: In contrast to other evaluated methods, ComHub consistently scored among the top performing methods on data from different sources. Lastly, we implemented ComHub to work with both predefined networks and to perform stand-alone network inference, which will make the method generally applicable.

Whole-genome sequencing and gene network modules predict gemcitabine/carboplatin-induced myelosuppression in non-small cell lung cancer patients.

Gemcitabine/carboplatin chemotherapy commonly induces myelosuppression, including neutropenia, leukopenia, and thrombocytopenia. Predicting patients at risk of these adverse drug reactions (ADRs) and adjusting treatments accordingly is a long-term goal of personalized medicine. This study used whole-genome sequencing (WGS) of blood samples from 96 gemcitabine/carboplatin-treated non-small cell lung cancer (NSCLC) patients and gene network modules for predicting myelosuppression. Association of genetic variants in PLINK found 4594, 5019, and 5066 autosomal SNVs/INDELs with p ≤ 1 × 10-3 for neutropenia, leukopenia, and thrombocytopenia, respectively. Based on the SNVs/INDELs we identified the toxicity module, consisting of 215 unique overlapping genes inferred from MCODE-generated gene network modules of 350, 345, and 313 genes, respectively. These module genes showed enrichment for differentially expressed genes in rat bone marrow, human bone marrow, and human cell lines exposed to carboplatin and gemcitabine (p < 0.05). Then using 80% of the patients as training data, random LASSO reduced the number of SNVs/INDELs in the toxicity module into a feasible prediction model consisting of 62 SNVs/INDELs that accurately predict both the training and the test (remaining 20%) data with high (CTCAE 3-4) and low (CTCAE 0-1) maximal myelosuppressive toxicity completely, with the receiver-operating characteristic (ROC) area under the curve (AUC) of 100%. The present study shows how WGS, gene network modules, and random LASSO can be used to develop a feasible and tested model for predicting myelosuppressive toxicity. Although the proposed model predicts myelosuppression in this study, further evaluation in other studies is required to determine its reproducibility, usability, and clinical effect.

MODifieR: an Ensemble R Package for Inference of Disease Modules from Transcriptomics Networks.

MOTIVATION: Complex diseases are due to the dense interactions of many disease-associated factors that dysregulate genes that in turn form the so-called disease modules, which have shown to be a powerful concept for understanding pathological mechanisms. There exist many disease module inference methods that rely on somewhat different assumptions, but there is still no gold standard or best-performing method. Hence, there is a need for combining these methods to generate robust disease modules. RESULTS: We developed MODule IdentiFIER (MODifieR), an ensemble R package of nine disease module inference methods from transcriptomics networks. MODifieR uses standardized input and output allowing the possibility to combine individual modules generated from these methods into more robust disease-specific modules, contributing to a better understanding of complex diseases. AVAILABILITY AND IMPLEMENTATION: MODifieR is available under the GNU GPL license and can be freely downloaded from https://gitlab.com/Gustafsson-lab/MODifieR and as a Docker image from https://hub.docker.com/r/ddeweerd/modifier. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

miRFA: an automated pipeline for microRNA functional analysis with correlation support from TCGA and TCPA expression data in pancreatic cancer.

BACKGROUND: MicroRNAs (miRNAs) are small RNAs that regulate gene expression at a post-transcriptional level and are emerging as potentially important biomarkers for various disease states, including pancreatic cancer. In silico-based functional analysis of miRNAs usually consists of miRNA target prediction and functional enrichment analysis of miRNA targets. Since miRNA target prediction methods generate a large number of false positive target genes, further validation to narrow down interesting candidate miRNA targets is needed. One commonly used method correlates miRNA and mRNA expression to assess the regulatory effect of a particular miRNA. The aim of this study was to build a bioinformatics pipeline in R for miRNA functional analysis including correlation analyses between miRNA expression levels and its targets on mRNA and protein expression levels available from the cancer genome atlas (TCGA) and the cancer proteome atlas (TCPA). TCGA-derived expression data of specific mature miRNA isoforms from pancreatic cancer tissue was used. RESULTS: Fifteen circulating miRNAs with significantly altered expression levels detected in pancreatic cancer patients were queried separately in the pipeline. The pipeline generated predicted miRNA target genes, enriched gene ontology (GO) terms and Kyoto encyclopedia of genes and genomes (KEGG) pathways. Predicted miRNA targets were evaluated by correlation analyses between each miRNA and its predicted targets. MiRNA functional analysis in combination with Kaplan-Meier survival analysis suggest that hsa-miR-885-5p could act as a tumor suppressor and should be validated as a potential prognostic biomarker in pancreatic cancer. CONCLUSIONS: Our miRNA functional analysis (miRFA) pipeline can serve as a valuable tool in biomarker discovery involving mature miRNAs associated with pancreatic cancer and could be developed to cover additional cancer types. Results for all mature miRNAs in TCGA pancreatic adenocarcinoma dataset can be studied and downloaded through a shiny web application at https://emmbor.shinyapps.io/mirfa/ .

Batch-normalization of cerebellar and medulloblastoma gene expression datasets utilizing empirically defined negative control genes.

MOTIVATION: Medulloblastoma (MB) is a brain cancer predominantly arising in children. Roughly 70% of patients are cured today, but survivors often suffer from severe sequelae. MB has been extensively studied by molecular profiling, but often in small and scattered cohorts. To improve cure rates and reduce treatment side effects, accurate integration of such data to increase analytical power will be important, if not essential. RESULTS: We have integrated 23 transcription datasets, spanning 1350 MB and 291 normal brain samples. To remove batch effects, we combined the Removal of Unwanted Variation (RUV) method with a novel pipeline for determining empirical negative control genes and a panel of metrics to evaluate normalization performance. The documented approach enabled the removal of a majority of batch effects, producing a large-scale, integrative dataset of MB and cerebellar expression data. The proposed strategy will be broadly applicable for accurate integration of data and incorporation of normal reference samples for studies of various diseases. We hope that the integrated dataset will improve current research in the field of MB by allowing more large-scale gene expression analyses. AVAILABILITY AND IMPLEMENTATION: The RUV-normalized expression data is available through the Gene Expression Omnibus (GEO; https://www.ncbi.nlm.nih.gov/geo/) and can be accessed via the GSE series number GSE124814. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Gene Co-Expression Network Analysis for Identifying Modules and Functionally Enriched Pathways in Type 1 Diabetes.

Type 1 diabetes (T1D) is a complex disease, caused by the autoimmune destruction of the insulin producing pancreatic beta cells, resulting in the body's inability to produce insulin. While great efforts have been put into understanding the genetic and environmental factors that contribute to the etiology of the disease, the exact molecular mechanisms are still largely unknown. T1D is a heterogeneous disease, and previous research in this field is mainly focused on the analysis of single genes, or using traditional gene expression profiling, which generally does not reveal the functional context of a gene associated with a complex disorder. However, network-based analysis does take into account the interactions between the diabetes specific genes or proteins and contributes to new knowledge about disease modules, which in turn can be used for identification of potential new biomarkers for T1D. In this study, we analyzed public microarray data of T1D patients and healthy controls by applying a systems biology approach that combines network-based Weighted Gene Co-Expression Network Analysis (WGCNA) with functional enrichment analysis. Novel co-expression gene network modules associated with T1D were elucidated, which in turn provided a basis for the identification of potential pathways and biomarker genes that may be involved in development of T1D.

Using expression profiling to understand the effects of chronic cadmium exposure on MCF-7 breast cancer cells.

Cadmium is a metalloestrogen known to activate the estrogen receptor and promote breast cancer cell growth. Previous studies have implicated cadmium in the development of more malignant tumors; however the molecular mechanisms behind this cadmium-induced malignancy remain elusive. Using clonal cell lines derived from exposing breast cancer cells to cadmium for over 6 months (MCF-7-Cd4, -Cd6, -Cd7, -Cd8 and -Cd12), this study aims to identify gene expression signatures associated with chronic cadmium exposure. Our results demonstrate that prolonged cadmium exposure does not merely result in the deregulation of genes but actually leads to a distinctive expression profile. The genes deregulated in cadmium-exposed cells are involved in multiple biological processes (i.e. cell growth, apoptosis, etc.) and molecular functions (i.e. cadmium/metal ion binding, transcription factor activity, etc.). Hierarchical clustering demonstrates that the five clonal cadmium cell lines share a common gene expression signature of breast cancer associated genes, clearly differentiating control cells from cadmium exposed cells. The results presented in this study offer insights into the cellular and molecular impacts of cadmium on breast cancer and emphasize the importance of studying chronic cadmium exposure as one possible mechanism of promoting breast cancer progression.