Computational methods to identify bimodal gene expression and facilitate personalized treatment in cancer patients

mtDNA Single-Nucleotide Variants Associated with Type 2 Diabetes

Type 2 diabetes (T2D) is a chronic systemic disease with a complex etiology, characterized by insulin resistance and mitochondrial dysfunction in various cell tissues. To explore this relationship, we conducted a secondary analysis of complete mtDNA sequences from 1261 T2D patients and 1105 control individuals. Our findings revealed significant associations between certain single-nucleotide polymorphisms (SNPs) and T2D. Notably, the variants m.1438A>G (rs2001030) (controls: 32 [27.6%], T2D: 84 [72.4%]; OR: 2.46; 95%CI: 1.64-3.78; p < 0.001), m.14766C>T (rs193302980) (controls: 498 [36.9%], T2D: 853 [63.1%]; OR: 2.57, 95%CI: 2.18-3.04, p < 0.001), and m.16519T>C (rs3937033) (controls: 363 [43.4%], T2D: 474 [56.6%]; OR: 1.24, 95%CI: 1.05-1.47, p = 0.012) were significantly associated with the likelihood of developing diabetes. The variant m.16189T>C (rs28693675), which has been previously documented in several studies across diverse populations, showed no association with T2D in our analysis (controls: 148 [13.39] T2D: 171 [13.56%]; OR: 1.03; 95%CI: 0.815-1.31; p = 0.83). These results provide evidence suggesting a link between specific mtDNA polymorphisms and T2D, possibly related to association rules, topological patterns, and three-dimensional conformations associated with regions where changes occur, rather than specific point mutations in the sequence.

Cell Decision Making through the Lens of Bayesian Learning

Cell decision making refers to the process by which cells gather information from their local microenvironment and regulate their internal states to create appropriate responses. Microenvironmental cell sensing plays a key role in this process. Our hypothesis is that cell decision-making regulation is dictated by Bayesian learning. In this article, we explore the implications of this hypothesis for internal state temporal evolution. By using a timescale separation between internal and external variables on the mesoscopic scale, we derive a hierarchical Fokker-Planck equation for cell-microenvironment dynamics. By combining this with the Bayesian learning hypothesis, we find that changes in microenvironmental entropy dominate the cell state probability distribution. Finally, we use these ideas to understand how cell sensing impacts cell decision making. Notably, our formalism allows us to understand cell state dynamics even without exact biochemical information about cell sensing processes by considering a few key parameters.

Neurocan expression associates with better survival and viral positivity in Merkel cell carcinoma

Merkel cell carcinoma (MCC) is a rare cutaneous neuroendocrine carcinoma that is frequently divided into Merkel cell polyomavirus negative and positive tumors due their distinct genomic and transcriptomic profiles, and disease outcomes. Although some prognostic factors in MCC are known, tumorigenic pathways, which that explain outcome differences in MCC are not fully understood. We investigated transcriptomes of 110 tissue samples of a formalin-fixed, paraffin-embedded MCC series by RNA sequencing to identify genes showing a bimodal expression pattern and predicting outcome in cancer and that potentially could play a role in tumorigenesis. We discovered 19 genes among which IGHM, IGKC, NCAN, OTOF, and USH2A were associated also with overall survival (all p-values < 0.05). From these genes, NCAN (neurocan) expression was detected in all 144 MCC samples by immunohistochemistry. Increased NCAN expression was associated with presence of Merkel cell polyomavirus DNA (p = 0.001) and viral large T antigen expression in tumor tissue (p = 0.004) and with improved MCC-specific survival (p = 0.027) and overall survival (p = 0.034). We conclude that NCAN expression is common in MCC, and further studies are warranted to investigate its role in MCC tumorigenesis.

An integrated approach to identify bimodal genes associated with prognosis in câncer

Bimodal gene expression (where a gene expression distribution has two maxima) is associated with phenotypic diversity in different biological systems. A critical issue, thus, is the integration of expression and phenotype data to identify genuine associations. Here, we developed tools that allow both: i) the identification of genes with bimodal gene expression and ii) their association with prognosis in cancer patients from The Cancer Genome Atlas (TCGA). Bimodality was observed for 554 genes in expression data from 25 tumor types. Furthermore, 96 of these genes presented different prognosis when patients belonging to the two expression peaks were compared. The software to execute the method and the corresponding documentation are available at the Data access section.

Functional comparison of exome capture-based methods for transcriptomic profiling of formalin-fixed paraffin-embedded tumors

The availability of fresh frozen (FF) tissue is a barrier for implementing RNA sequencing (RNA-seq) in the clinic. The majority of clinical samples are stored as formalin-fixed, paraffin-embedded (FFPE) tissues. Exome capture platforms have been developed for RNA-seq from FFPE samples. However, these methods have not been systematically compared. We performed transcriptomic analysis of 32 FFPE tumor samples from 11 patients using three exome capture-based methods: Agilent SureSelect V6, TWIST NGS Exome, and IDT XGen Exome Research Panel. We compared these methods to the TruSeq RNA-seq of fresh frozen (FF-TruSeq) tumor samples from the same patients. We assessed the recovery of clinically relevant biological features. The Spearman's correlation coefficients between the global expression profiles of the three capture-based methods from FFPE and matched FF-TruSeq were high (rho = 0.72-0.9, p < 0.05). A significant correlation between the expression of key immune genes between individual capture-based methods and FF-TruSeq (rho = 0.76-0.88, p < 0.05) was observed. All exome capture-based methods reliably detected outlier expression of actionable gene transcripts, including ERBB2, MET, NTRK1, and PPARG. In urothelial cancer samples, the Agilent assay was associated with the highest molecular subtype concordance with FF-TruSeq (Cohen's k = 0.7, p < 0.01). The Agilent and IDT assays detected all the clinically relevant fusions that were initially identified in FF-TruSeq. All FFPE exome capture-based methods had comparable performance and concordance with FF-TruSeq. Our findings will enable the implementation of RNA-seq in the clinic to guide precision oncology approaches.

A comprehensive characterization of the cell-free transcriptome reveals tissue- and subtype-specific biomarkers for cancer detection

Cell-free RNA (cfRNA) is a promising analyte for cancer detection. However, a comprehensive assessment of cfRNA in individuals with and without cancer has not been conducted. We perform the first transcriptome-wide characterization of cfRNA in cancer (stage III breast [n = 46], lung [n = 30]) and non-cancer (n = 89) participants from the Circulating Cell-free Genome Atlas (NCT02889978). Of 57,820 annotated genes, 39,564 (68%) are not detected in cfRNA from non-cancer individuals. Within these low-noise regions, we identify tissue- and cancer-specific genes, defined as “dark channel biomarker” (DCB) genes, that are recurrently detected in individuals with cancer. DCB levels in plasma correlate with tumor shedding rate and RNA expression in matched tissue, suggesting that DCBs with high expression in tumor tissue could enhance cancer detection in patients with low levels of circulating tumor DNA. Overall, cfRNA provides a unique opportunity to detect cancer, predict the tumor tissue of origin, and determine the cancer subtype.

Cell-free RNA (cfRNA) is a promising analyte for cancer diagnosis. Here, the authors determine the baseline cell-free transcriptome in the absence of cancer and identify tissue- and subtype-specific cfRNA biomarkers in breast and lung cancer patients.

SLFN11 informs on standard of care and novel treatments in a wide range of cancer models

BACKGROUND

Schlafen 11 (SLFN11) has been linked with response to DNA-damaging agents (DDA) and PARP inhibitors. An in-depth understanding of several aspects of its role as a biomarker in cancer is missing, as is a comprehensive analysis of the clinical significance of SLFN11 as a predictive biomarker to DDA and/or DNA damage-response inhibitor (DDRi) therapies.

METHODS

We used a multidisciplinary effort combining specific immunohistochemistry, pharmacology tests, anticancer combination therapies and mechanistic studies to assess SLFN11 as a potential biomarker for stratification of patients treated with several DDA and/or DDRi in the preclinical and clinical setting.

RESULTS

SLFN11 protein associated with both preclinical and patient treatment response to DDA, but not to non-DDA or DDRi therapies, such as WEE1 inhibitor or olaparib in breast cancer. SLFN11-low/absent cancers were identified across different tumour types tested. Combinations of DDA with DDRi targeting the replication-stress response (ATR, CHK1 and WEE1) could re-sensitise SLFN11-absent/low cancer models to the DDA treatment and were effective in upper gastrointestinal and genitourinary malignancies.

CONCLUSION

SLFN11 informs on the standard of care chemotherapy based on DDA and the effect of selected combinations with ATR, WEE1 or CHK1 inhibitor in a wide range of cancer types and models.

Comparison of Methods for Feature Selection in Clustering of High-Dimensional RNA-Sequencing Data to Identify Cancer Subtypes

Cancer subtype identification is important to facilitate cancer diagnosis and select effective treatments. Clustering of cancer patients based on high-dimensional RNA-sequencing data can be used to detect novel subtypes, but only a subset of the features (e.g., genes) contains information related to the cancer subtype. Therefore, it is reasonable to assume that the clustering should be based on a set of carefully selected features rather than all features. Several feature selection methods have been proposed, but how and when to use these methods are still poorly understood. Thirteen feature selection methods were evaluated on four human cancer data sets, all with known subtypes (gold standards), which were only used for evaluation. The methods were characterized by considering mean expression and standard deviation (SD) of the selected genes, the overlap with other methods and their clustering performance, obtained comparing the clustering result with the gold standard using the adjusted Rand index (ARI). The results were compared to a supervised approach as a positive control and two negative controls in which either a random selection of genes or all genes were included. For all data sets, the best feature selection approach outperformed the negative control and for two data sets the gain was substantial with ARI increasing from (-0.01, 0.39) to (0.66, 0.72), respectively. No feature selection method completely outperformed the others but using the dip-rest statistic to select 1000 genes was overall a good choice. The commonly used approach, where genes with the highest SDs are selected, did not perform well in our study.

Disease-related cellular protein networks differentially affected under different EGFR mutations in lung adenocarcinoma

It is unclear how epidermal growth factor receptor EGFR major driver mutations (L858R or Ex19del) affect downstream molecular networks and pathways. This study aimed to provide information on the influences of these mutations. The study assessed 36 protein expression profiles of lung adenocarcinoma (Ex19del, nine; L858R, nine; no Ex19del/L858R, 18). Weighted gene co-expression network analysis together with analysis of variance-based screening identified 13 co-expressed modules and their eigen proteins. Pathway enrichment analysis for the Ex19del mutation demonstrated involvement of SUMOylation, epithelial and mesenchymal transition, ERK/mitogen-activated protein kinase signalling via phosphorylation and Hippo signalling. Additionally, analysis for the L858R mutation identified various pathways related to cancer cell survival and death. With regard to the Ex19del mutation, ROCK, RPS6KA1, ARF1, IL2RA and several ErbB pathways were upregulated, whereas AURK and GSKIP were downregulated. With regard to the L858R mutation, RB1, TSC22D3 and DOCK1 were downregulated, whereas various networks, including VEGFA, were moderately upregulated. In all mutation types, CD80/CD86 (B7), MHC, CIITA and IFGN were activated, whereas CD37 and SAFB were inhibited. Costimulatory immune-checkpoint pathways by B7/CD28 were mainly activated, whereas those by PD-1/PD-L1 were inhibited. Our findings may help identify potential therapeutic targets and develop therapeutic strategies to improve patient outcomes.