NCAE: data-driven representations using a deep network-coherent DNA methylation autoencoder identify robust disease and risk factor signatures.

Precision medicine relies on the identification of robust disease and risk factor signatures from omics data. However, current knowledge-driven approaches may overlook novel or unexpected phenomena due to the inherent biases in biological knowledge. In this study, we present a data-driven signature discovery workflow for DNA methylation analysis utilizing network-coherent autoencoders (NCAEs) with biologically relevant latent embeddings. First, we explored the architecture space of autoencoders trained on a large-scale pan-tissue compendium (n = 75 272) of human epigenome-wide association studies. We observed the emergence of co-localized patterns in the deep autoencoder latent space representations that corresponded to biological network modules. We determined the NCAE configuration with the strongest co-localization and centrality signals in the human protein interactome. Leveraging the NCAE embeddings, we then trained interpretable deep neural networks for risk factor (aging, smoking) and disease (systemic lupus erythematosus) prediction and classification tasks. Remarkably, our NCAE embedding-based models outperformed existing predictors, revealing novel DNA methylation signatures enriched in gene sets and pathways associated with the studied condition in each case. Our data-driven biomarker discovery workflow provides a generally applicable pipeline to capture relevant risk factor and disease information. By surpassing the limitations of knowledge-driven methods, our approach enhances the understanding of complex epigenetic processes, facilitating the development of more effective diagnostic and therapeutic strategies.

Generation and analysis of context-specific genome-scale metabolic models derived from single-cell RNA-Seq data.

Single-cell RNA sequencing combined with genome-scale metabolic models (GEMs) has the potential to unravel the differences in metabolism across both cell types and cell states but requires new computational methods. Here, we present a method for generating cell-type-specific genome-scale models from clusters of single-cell RNA-Seq profiles. Specifically, we developed a method to estimate the minimum number of cells required to pool to obtain stable models, a bootstrapping strategy for estimating statistical inference, and a faster version of the task-driven integrative network inference for tissues algorithm for generating context-specific GEMs. In addition, we evaluated the effect of different RNA-Seq normalization methods on model topology and differences in models generated from single-cell and bulk RNA-Seq data. We applied our methods on data from mouse cortex neurons and cells from the tumor microenvironment of lung cancer and in both cases found that almost every cell subtype had a unique metabolic profile. In addition, our approach was able to detect cancer-associated metabolic differences between cancer cells and healthy cells, showcasing its utility. We also contextualized models from 202 single-cell clusters across 19 human organs using data from Human Protein Atlas and made these available in the web portal Metabolic Atlas, thereby providing a valuable resource to the scientific community. With the ever-increasing availability of single-cell RNA-Seq datasets and continuously improved GEMs, their combination holds promise to become an important approach in the study of human metabolism.

Modeling glioblastoma heterogeneity as a dynamic network of cell states.

Tumor cell heterogeneity is a crucial characteristic of malignant brain tumors and underpins phenomena such as therapy resistance and tumor recurrence. Advances in single-cell analysis have enabled the delineation of distinct cellular states of brain tumor cells, but the time-dependent changes in such states remain poorly understood. Here, we construct quantitative models of the time-dependent transcriptional variation of patient-derived glioblastoma (GBM) cells. We build the models by sampling and profiling barcoded GBM cells and their progeny over the course of 3 weeks and by fitting a mathematical model to estimate changes in GBM cell states and their growth rates. Our model suggests a hierarchical yet plastic organization of GBM, where the rates and patterns of cell state switching are partly patient-specific. Therapeutic interventions produce complex dynamic effects, including inhibition of specific states and altered differentiation. Our method provides a general strategy to uncover time-dependent changes in cancer cells and offers a way to evaluate and predict how therapy affects cell state composition.

Dose-response relationships of intestinal organs and excessive mucus discharge after gynaecological radiotherapy.

BACKGROUND: The study aims to determine possible dose-volume response relationships between the rectum, sigmoid colon and small intestine and the 'excessive mucus discharge' syndrome after pelvic radiotherapy for gynaecological cancer. METHODS AND MATERIALS: From a larger cohort, 98 gynaecological cancer survivors were included in this study. These survivors, who were followed for 2 to 14 years, received external beam radiation therapy but not brachytherapy and not did not have stoma. Thirteen of the 98 developed excessive mucus discharge syndrome. Three self-assessed symptoms were weighted together to produce a score interpreted as 'excessive mucus discharge' syndrome based on the factor loadings from factor analysis. The dose-volume histograms (DVHs) for rectum, sigmoid colon, small intestine for each survivor were exported from the treatment planning systems. The dose-volume response relationships for excessive mucus discharge and each organ at risk were estimated by fitting the data to the Probit, RS, LKB and gEUD models. RESULTS: The small intestine was found to have steep dose-response curves, having estimated dose-response parameters: γ50: 1.28, 1.23, 1.32, D50: 61.6, 63.1, 60.2 for Probit, RS and LKB respectively. The sigmoid colon (AUC: 0.68) and the small intestine (AUC: 0.65) had the highest AUC values. For the small intestine, the DVHs for survivors with and without excessive mucus discharge were well separated for low to intermediate doses; this was not true for the sigmoid colon. Based on all results, we interpret the results for the small intestine to reflect a relevant link. CONCLUSION: An association was found between the mean dose to the small intestine and the occurrence of 'excessive mucus discharge'. When trying to reduce and even eliminate the incidence of 'excessive mucus discharge', it would be useful and important to separately delineate the small intestine and implement the dose-response estimations reported in the study.

DSAVE: Detection of misclassified cells in single-cell RNA-Seq data.

Single-cell RNA sequencing has become a valuable tool for investigating cell types in complex tissues, where clustering of cells enables the identification and comparison of cell populations. Although many studies have sought to develop and compare different clustering approaches, a deeper investigation into the properties of the resulting populations is lacking. Specifically, the presence of misclassified cells can influence downstream analyses, highlighting the need to assess subpopulation purity and to detect such cells. We developed DSAVE (Down-SAmpling based Variation Estimation), a method to evaluate the purity of single-cell transcriptome clusters and to identify misclassified cells. The method utilizes down-sampling to eliminate differences in sampling noise and uses a log-likelihood based metric to help identify misclassified cells. In addition, DSAVE estimates the number of cells needed in a population to achieve a stable average gene expression profile within a certain gene expression range. We show that DSAVE can be used to find potentially misclassified cells that are not detectable by similar tools and reveal the cause of their divergence from the other cells, such as differing cell state or cell type. With the growing use of single-cell RNA-seq, we foresee that DSAVE will be an increasingly useful tool for comparing and purifying subpopulations in single-cell RNA-Seq datasets.

Sources of variation in cell-type RNA-Seq profiles.

Cell-type specific gene expression profiles are needed for many computational methods operating on bulk RNA-Seq samples, such as deconvolution of cell-type fractions and digital cytometry. However, the gene expression profile of a cell type can vary substantially due to both technical factors and biological differences in cell state and surroundings, reducing the efficacy of such methods. Here, we investigated which factors contribute most to this variation. We evaluated different normalization methods, quantified the variance explained by different factors, evaluated the effect on deconvolution of cell type fractions, and examined the differences between UMI-based single-cell RNA-Seq and bulk RNA-Seq. We investigated a collection of publicly available bulk and single-cell RNA-Seq datasets containing B and T cells, and found that the technical variation across laboratories is substantial, even for genes specifically selected for deconvolution, and this variation has a confounding effect on deconvolution. Tissue of origin is also a substantial factor, highlighting the challenge of using cell type profiles derived from blood with mixtures from other tissues. We also show that much of the differences between UMI-based single-cell and bulk RNA-Seq methods can be explained by the number of read duplicates per mRNA molecule in the single-cell sample. Our work shows the importance of either matching or correcting for technical factors when creating cell-type specific gene expression profiles that are to be used together with bulk samples.

Molecular natural history of breast cancer: Leveraging transcriptomics to predict breast cancer progression and aggressiveness.

BACKGROUND: Characterizing breast cancer progression and aggressiveness relies on categorical descriptions of tumor stage and grade. Interpreting these categorical descriptions is challenging because stage convolutes the size and spread of the tumor and no consensus exists to define high/low grade tumors. METHODS: We address this challenge of heterogeneity in patient-specific cancer samples by adapting and applying several tools originally created for understanding heterogeneity and phenotype development in single cells (specifically, single-cell topological data analysis and Wanderlust) to create a continuous metric describing breast cancer progression using bulk RNA-seq samples from individual patient tumors. We also created a linear regression-based method to predict tumor aggressiveness in vivo from bulk RNA-seq data. RESULTS: We found that breast cancer proceeds along three convergent phenotype trajectories: luminal, HER2-enriched, and basal-like. Furthermore, 31 296 genes (for luminal cancers), 17 827 genes (for HER2-enriched), and 18 505 genes (for basal-like) are dynamically differentially expressed during breast cancer progression. Across progression trajectories, our results show that expression of genes related to ADP-ribosylation decreased as tumors progressed (while PARP1 and PARP2 increased or remained stable), suggesting the potential for a differential response to PARP inhibitors based on cancer progression. Additionally, we developed a 132-gene expression regression equation to predict mitotic index and a 23-gene expression regression equation to predict growth rate from a single breast cancer biopsy. CONCLUSION: Our results suggest that breast cancer dynamically changes during disease progression, and growth rate of the cancer cells is associated with distinct transcriptional profiles.

Integrative discovery of treatments for high-risk neuroblastoma.

Despite advances in the molecular exploration of paediatric cancers, approximately 50% of children with high-risk neuroblastoma lack effective treatment. To identify therapeutic options for this group of high-risk patients, we combine predictive data mining with experimental evaluation in patient-derived xenograft cells. Our proposed algorithm, TargetTranslator, integrates data from tumour biobanks, pharmacological databases, and cellular networks to predict how targeted interventions affect mRNA signatures associated with high patient risk or disease processes. We find more than 80 targets to be associated with neuroblastoma risk and differentiation signatures. Selected targets are evaluated in cell lines derived from high-risk patients to demonstrate reversal of risk signatures and malignant phenotypes. Using neuroblastoma xenograft models, we establish CNR2 and MAPK8 as promising candidates for the treatment of high-risk neuroblastoma. We expect that our method, available as a public tool (targettranslator.org), will enhance and expedite the discovery of risk-associated targets for paediatric and adult cancers.

Digital twins to personalize medicine.

Personalized medicine requires the integration and processing of vast amounts of data. Here, we propose a solution to this challenge that is based on constructing Digital Twins. These are high-resolution models of individual patients that are computationally treated with thousands of drugs to find the drug that is optimal for the patient.

Dose-response relationships of the sigmoid for urgency syndrome after gynecological radiotherapy.

PURPOSE: To find out what organs and doses are most relevant for 'radiation-induced urgency syndrome' in order to derive the corresponding dose-response relationships as an aid for avoiding the syndrome in the future. MATERIAL AND METHODS: From a larger group of gynecological cancer survivors followed-up 2-14 years, we identified 98 whom had undergone external beam radiation therapy but not brachytherapy and not having a stoma. Of those survivors, 24 developed urgency syndrome. Based on the loading factor from a factor analysis, and symptom frequency, 15 symptoms were weighted together to a score interpreted as the intensity of radiation-induced urgency symptom. On reactivated dose plans, we contoured the small intestine, sigmoid colon and the rectum (separate from the anal-sphincter region) and we exported the dose-volume histograms for each survivor. Dose-response relationships from respective risk organ and urgency syndrome were estimated by fitting the data to the Probit, RS, LKB and gEUD models. RESULTS: The rectum and sigmoid colon have steep dose-response relationships for urgency syndrome for Probit, RS and LKB. The dose-response parameters for the rectum were D50: 51.3, 51.4, and 51.3 Gy, γ50 = 1.19 for all models, s was 7.0e-09 for RS and n was 9.9 × 107 for LKB. For Sigmoid colon, D50 were 51.6, 51.6, and 51.5 Gy, γ50 were 1.20, 1.25, and 1.27, s was 2.8 for RS and n was 0.079 for LKB. CONCLUSIONS: Primarily the dose to sigmoid colon as well as the rectum is related to urgency syndrome among gynecological cancer survivors. Separate delineation of the rectum and sigmoid colon in order to incorporate the dose-response results may aid in reduction of the incidence of the urgency syndrome.

LASSIM-A network inference toolbox for genome-wide mechanistic modeling.

Recent technological advancements have made time-resolved, quantitative, multi-omics data available for many model systems, which could be integrated for systems pharmacokinetic use. Here, we present large-scale simulation modeling (LASSIM), which is a novel mathematical tool for performing large-scale inference using mechanistically defined ordinary differential equations (ODE) for gene regulatory networks (GRNs). LASSIM integrates structural knowledge about regulatory interactions and non-linear equations with multiple steady state and dynamic response expression datasets. The rationale behind LASSIM is that biological GRNs can be simplified using a limited subset of core genes that are assumed to regulate all other gene transcription events in the network. The LASSIM method is implemented as a general-purpose toolbox using the PyGMO Python package to make the most of multicore computers and high performance clusters, and is available at https://gitlab.com/Gustafsson-lab/lassim. As a method, LASSIM works in two steps, where it first infers a non-linear ODE system of the pre-specified core gene expression. Second, LASSIM in parallel optimizes the parameters that model the regulation of peripheral genes by core system genes. We showed the usefulness of this method by applying LASSIM to infer a large-scale non-linear model of naïve Th2 cell differentiation, made possible by integrating Th2 specific bindings, time-series together with six public and six novel siRNA-mediated knock-down experiments. ChIP-seq showed significant overlap for all tested transcription factors. Next, we performed novel time-series measurements of total T-cells during differentiation towards Th2 and verified that our LASSIM model could monitor those data significantly better than comparable models that used the same Th2 bindings. In summary, the LASSIM toolbox opens the door to a new type of model-based data analysis that combines the strengths of reliable mechanistic models with truly systems-level data. We demonstrate the power of this approach by inferring a mechanistically motivated, genome-wide model of the Th2 transcription regulatory system, which plays an important role in several immune related diseases.

Late radiation-induced bowel syndromes, tobacco smoking, age at treatment and time since treatment - gynecological cancer survivors.

BACKGROUND: It is unknown whether smoking; age at time of radiotherapy or time since radiotherapy influence the intensity of late radiation-induced bowel syndromes. MATERIAL AND METHODS: We have previously identified 28 symptoms decreasing bowel health among 623 gynecological-cancer survivors (three to twelve years after radiotherapy) and 344 matched population-based controls. The 28 symptoms were grouped into five separate late bowel syndromes through factor analysis. Here, we related possible predictors of bowel health to syndrome intensity, by combining factor analysis weights and symptom frequency on a person-incidence scale. RESULTS: A strong (p < .001) association between smoking and radiation-induced urgency syndrome was found with a syndrome intensity (normalized factor score) of 0.4 (never smoker), 1.2 (former smoker) and 2.5 (current smoker). Excessive gas discharge was also related to smoking (p = .001). Younger age at treatment resulted in a higher intensity, except for the leakage syndrome. For the urgency syndrome, intensity decreased with time since treatment. CONCLUSIONS: Smoking aggravates the radiation-induced urgency syndrome and excessive gas discharge syndrome. Smoking cessation may promote bowel health among gynecological-cancer survivors. Furthermore, by understanding the mechanism for the decline in urgency-syndrome intensity over time, we may identify new strategies for prevention and alleviation.

Identifying radiation-induced survivorship syndromes affecting bowel health in a cohort of gynecological cancer survivors.

BACKGROUND: During radiotherapy unwanted radiation to normal tissue surrounding the tumor triggers survivorship diseases; we lack a nosology for radiation-induced survivorship diseases that decrease bowel health and we do not know which symptoms are related to which diseases. METHODS: Gynecological-cancer survivors were followed-up two to 15 years after having undergone radiotherapy; they reported in a postal questionnaire the frequency of 28 different symptoms related to bowel health. Population-based controls gave the same information. With a modified factor analysis, we determined the optimal number of factors, factor loadings for each symptom, factor-specific factor-loading cutoffs and factor scores. RESULTS: Altogether data from 623 survivors and 344 population-based controls were analyzed. Six factors best explain the correlation structure of the symptoms; for five of these a statistically significant difference (P< 0.001, Mann-Whitney U test) was found between survivors and controls concerning factor score quantiles. Taken together these five factors explain 42 percent of the variance of the symptoms. We interpreted these five factors as radiation-induced syndromes that may reflect distinct survivorship diseases. We obtained the following frequencies, defined as survivors having a factor loading above the 95 percent percentile of the controls, urgency syndrome (190 of 623, 30 percent), leakage syndrome (164 of 623, 26 percent), excessive gas discharge (93 of 623, 15 percent), excessive mucus discharge (102 of 623, 16 percent) and blood discharge (63 of 623, 10 percent). CONCLUSION: Late effects of radiotherapy include five syndromes affecting bowel health; studying them and identifying the underlying survivorship diseases, instead of the approximately 30 long-term symptoms they produce, will simplify the search for prevention, alleviation and elimination.

Integrative Modeling Reveals Annexin A2-mediated Epigenetic Control of Mesenchymal Glioblastoma.

Glioblastomas are characterized by transcriptionally distinct subtypes, but despite possible clinical relevance, their regulation remains poorly understood. The commonly used molecular classification systems for GBM all identify a subtype with high expression of mesenchymal marker transcripts, strongly associated with invasive growth. We used a comprehensive data-driven network modeling technique (augmented sparse inverse covariance selection, aSICS) to define separate genomic, epigenetic, and transcriptional regulators of glioblastoma subtypes. Our model identified Annexin A2 (ANXA2) as a novel methylation-controlled positive regulator of the mesenchymal subtype. Subsequent evaluation in two independent cohorts established ANXA2 expression as a prognostic factor that is dependent on ANXA2 promoter methylation. ANXA2 knockdown in primary glioblastoma stem cell-like cultures suppressed known mesenchymal master regulators, and abrogated cell proliferation and invasion. Our results place ANXA2 at the apex of a regulatory cascade that determines glioblastoma mesenchymal transformation and validate aSICS as a general methodology to uncover regulators of cancer subtypes.

A Factor Analysis Approach for Clustering Patient Reported Outcomes.

BACKGROUND: In the field of radiation oncology, the use of extensive patient reported outcomes is increasingly common to measure adverse side effects after radiotherapy in cancer patients. Factor analysis has the potential to identify an optimal number of latent factors (i.e., symptom groups). However, the ultimate goal of treatment response modeling is to understand the relationship between treatment variables such as radiation dose and symptom groups resulting from FA. Hence, it is crucial to identify clinically more relevant symptom groups and improved response variables from those symptom groups for a quantitative analysis. OBJECTIVES: The goal of this study is to design a computational method for finding clinically relevant symptom groups from PROs and to test associations between symptom groups and radiation dose. METHODS: We propose a novel approach where exploratory factor analysis is followed by confirmatory factor analysis to determine the relevant number of symptom groups. We also propose to use a combination of symptoms in a symptom group identified as a new response variable in linear regression analysis to investigate the relationship between the symptom group and dose-volume variables. RESULTS: We analyzed patient-reported gastrointestinal symptom profiles from 3 datasets in prostate cancer patients treated with radiotherapy. The final structural model of each dataset was validated using the other two datasets and compared to four other existing FA methods. Our systematic EFA-CFA approach provided clinically more relevant solutions than other methods, resulting in new clinically relevant outcome variables that enabled a quantitative analysis. As a result, statistically significant correlations were found between some dose-volume variables to relevant anatomic structures and symptom groups identified by FA. CONCLUSIONS: Our proposed method can aid in the process of understanding PROs and provide a basis for improving our understanding of radiation-induced side effects.

Efficient exploration of pan-cancer networks by generalized covariance selection and interactive web content.

Statistical network modeling techniques are increasingly important tools to analyze cancer genomics data. However, current tools and resources are not designed to work across multiple diagnoses and technical platforms, thus limiting their applicability to comprehensive pan-cancer datasets such as The Cancer Genome Atlas (TCGA). To address this, we describe a new data driven modeling method, based on generalized Sparse Inverse Covariance Selection (SICS). The method integrates genetic, epigenetic and transcriptional data from multiple cancers, to define links that are present in multiple cancers, a subset of cancers, or a single cancer. It is shown to be statistically robust and effective at detecting direct pathway links in data from TCGA. To facilitate interpretation of the results, we introduce a publicly accessible tool (cancerlandscapes.org), in which the derived networks are explored as interactive web content, linked to several pathway and pharmacological databases. To evaluate the performance of the method, we constructed a model for eight TCGA cancers, using data from 3900 patients. The model rediscovered known mechanisms and contained interesting predictions. Possible applications include prediction of regulatory relationships, comparison of network modules across multiple forms of cancer and identification of drug targets.

DNA methylation changes separate allergic patients from healthy controls and may reflect altered CD4+ T-cell population structure.

Altered DNA methylation patterns in CD4(+) T-cells indicate the importance of epigenetic mechanisms in inflammatory diseases. However, the identification of these alterations is complicated by the heterogeneity of most inflammatory diseases. Seasonal allergic rhinitis (SAR) is an optimal disease model for the study of DNA methylation because of its well-defined phenotype and etiology. We generated genome-wide DNA methylation (N(patients) = 8, N(controls) = 8) and gene expression (N(patients) = 9, Ncontrols = 10) profiles of CD4(+) T-cells from SAR patients and healthy controls using Illumina's HumanMethylation450 and HT-12 microarrays, respectively. DNA methylation profiles clearly and robustly distinguished SAR patients from controls, during and outside the pollen season. In agreement with previously published studies, gene expression profiles of the same samples failed to separate patients and controls. Separation by methylation (N(patients) = 12, N(controls) = 12), but not by gene expression (N(patients) = 21, N(controls) = 21) was also observed in an in vitro model system in which purified PBMCs from patients and healthy controls were challenged with allergen. We observed changes in the proportions of memory T-cell populations between patients (N(patients) = 35) and controls (N(controls) = 12), which could explain the observed difference in DNA methylation. Our data highlight the potential of epigenomics in the stratification of immune disease and represents the first successful molecular classification of SAR using CD4(+) T cells.

Erratum: Music structure determines heart rate variability of singers.

[This corrects the article on p. 334 in vol. 4, PMID: 23847555.].

Searching for synergies: matrix algebraic approaches for efficient pair screening.

Functionally interacting perturbations, such as synergistic drugs pairs or synthetic lethal gene pairs, are of key interest in both pharmacology and functional genomics. However, to find such pairs by traditional screening methods is both time consuming and costly. We present a novel computational-experimental framework for efficient identification of synergistic target pairs, applicable for screening of systems with sizes on the order of current drug, small RNA or SGA (Synthetic Genetic Array) libraries (>1000 targets). This framework exploits the fact that the response of a drug pair in a given system, or a pair of genes' propensity to interact functionally, can be partly predicted by computational means from (i) a small set of experimentally determined target pairs, and (ii) pre-existing data (e.g. gene ontology, PPI) on the similarities between targets. Predictions are obtained by a novel matrix algebraic technique, based on cyclical projections onto convex sets. We demonstrate the efficiency of the proposed method using drug-drug interaction data from seven cancer cell lines and gene-gene interaction data from yeast SGA screens. Our protocol increases the rate of synergism discovery significantly over traditional screening, by up to 7-fold. Our method is easy to implement and could be applied to accelerate pair screening for both animal and microbial systems.

Music structure determines heart rate variability of singers.

Choir singing is known to promote wellbeing. One reason for this may be that singing demands a slower than normal respiration, which may in turn affect heart activity. Coupling of heart rate variability (HRV) to respiration is called Respiratory sinus arrhythmia (RSA). This coupling has a subjective as well as a biologically soothing effect, and it is beneficial for cardiovascular function. RSA is seen to be more marked during slow-paced breathing and at lower respiration rates (0.1 Hz and below). In this study, we investigate how singing, which is a form of guided breathing, affects HRV and RSA. The study comprises a group of healthy 18 year olds of mixed gender. The subjects are asked to; (1) hum a single tone and breathe whenever they need to; (2) sing a hymn with free, unguided breathing; and (3) sing a slow mantra and breathe solely between phrases. Heart rate (HR) is measured continuously during the study. The study design makes it possible to compare above three levels of song structure. In a separate case study, we examine five individuals performing singing tasks (1-3). We collect data with more advanced equipment, simultaneously recording HR, respiration, skin conductance and finger temperature. We show how song structure, respiration and HR are connected. Unison singing of regular song structures makes the hearts of the singers accelerate and decelerate simultaneously. Implications concerning the effect on wellbeing and health are discussed as well as the question how this inner entrainment may affect perception and behavior.