Search

Purpose: The genetic intratumoral heterogeneity observed in human osteosarcomas poses challenges for drug development and the study of cell fate, plasticity, and differentiation, which are processes linked to tumor grade, cell metastasis, and survival., Experimental Design: To pinpoint errors in osteosarcoma differentiation, we transcriptionally profiled 31,527 cells from a tissue-engineered model that directs mesenchymal stem cells toward adipogenic and osteoblastic fates. Incorporating preexisting chondrocyte data, we applied trajectory analysis and non-negative matrix factorization to generate the first human mesenchymal differentiation atlas., Results: This "roadmap" served as a reference to delineate the cellular composition of morphologically complex osteosarcoma tumors and quantify each cell's lineage commitment. Projecting a bulk RNA-sequencing osteosarcoma dataset onto this roadmap unveiled a correlation between a stem-like transcriptomic phenotype and poorer survival outcomes., Conclusions: Our study quantifies osteosarcoma differentiation and lineage, a prerequisite to better understanding lineage-specific differentiation bottlenecks that might someday be targeted therapeutically., (©2024 American Association for Cancer Research.)

Melanoma response to immune-modulating therapy remains incompletely characterized at the molecular level. In this study, we assess melanoma immunotherapy response using a multi-scale network approach to identify gene modules with coordinated gene expression in response to treatment. Using gene expression data of melanoma before and after treatment with nivolumab, we modeled gene expression changes in a correlation network and measured a key network geometric property, dynamic Ollivier-Ricci curvature, to distinguish critical edges within the network and reveal multi-scale treatment-response gene communities. Analysis identified six distinct gene modules corresponding to sets of genes interacting in response to immunotherapy. One module alone, overlapping with the nuclear factor kappa-B pathway (NFkB), was associated with improved patient survival and a positive clinical response to immunotherapy. This analysis demonstrates the usefulness of dynamic Ollivier-Ricci curvature as a general method for identifying information-sharing gene modules in cancer., (© 2024. The Author(s).)

Directionally sensitive radiomic features including the histogram of oriented gradient (HOG) have been shown to provide objective and quantitative measures for predicting disease outcomes in multiple cancers. However, radiomic features are sensitive to imaging variabilities including acquisition differences, imaging artifacts and noise, making them impractical for using in the clinic to inform patient care. We treat the problem of extracting robust local directionality features by mapping via optimal transport a given local image patch to an iso-intense patch of its mean. We decompose the transport map into sub-work costs each transporting in different directions. To test our approach, we evaluated the ability of the proposed approach to quantify tumor heterogeneity from magnetic resonance imaging (MRI) scans of brain glioblastoma multiforme, computed tomography (CT) scans of head and neck squamous cell carcinoma as well as longitudinal CT scans in lung cancer patients treated with immunotherapy. By considering the entropy difference of the extracted local directionality within tumor regions, we found that patients with higher entropy in their images, had significantly worse overall survival for all three datasets, which indicates that tumors that have images exhibiting flows in many directions may be more malignant. This may seem to reflect high tumor histologic grade or disorganization. Furthermore, by comparing the changes in entropy longitudinally using two imaging time points, we found patients with reduction in entropy from baseline CT are associated with longer overall survival (hazard ratio = 1.95, 95% confidence interval of 1.4-2.8, p = 1.65e-5). The proposed method provides a robust, training free approach to quantify the local directionality contained in images.

As a generalization of the optimal mass transport (OMT) approach of Benamou and Brenier's, the regularized optimal mass transport (rOMT) formulates a transport problem from an initial mass configuration to another with the optimality defined by the total kinetic energy, but subject to an advection-diffusion constraint equation. Both rOMT and the Benamou and Brenier's formulation require the total initial and final masses to be equal; mass is preserved during the entire transport process. However, for many applications, e.g., in dynamic image tracking, this constraint is rarely if ever satisfied. Therefore, we propose to employ an unbalanced version of rOMT to remove this constraint together with a detailed numerical solution procedure and applications to analyzing fluid flows in the brain., (© 2024. The Author(s).)

Network properties account for the complex relationship between genes, making it easier to identify complex patterns in their interactions. In this work, we leveraged these network properties for dual purposes. First, we clustered pediatric sarcoma tumors using network information flow as a similarity metric, computed by the Wasserstein distance. We demonstrate that this approach yields the best concordance with histological subtypes, validated against three state-of-the-art methods. Second, to identify molecular targets that would be missed by more conventional methods of analysis, we applied a novel unsupervised method to cluster gene interactomes represented as networks in pediatric sarcoma. RNA-Seq data were mapped to protein-level interactomes to construct weighted networks that were then subjected to a non-Euclidean, multi-scale geometric approach centered on a discrete notion of curvature. This provides a measure of the functional association among genes in the context of their connectivity. In confirmation of the validity of this method, hierarchical clustering revealed the characteristic EWSR1-FLI1 fusion in Ewing sarcoma. Furthermore, assessing the effects of in silico edge perturbations and simulated gene knockouts as quantified by changes in curvature, we found non-trivial gene associations not previously identified., (© 2024. The Author(s).)

The plasma cell cancer multiple myeloma (MM) varies significantly in genomic characteristics, response to therapy, and long-term prognosis. To investigate global interactions in MM, we combined a known protein interaction network with a large clinically annotated MM dataset. We hypothesized that an unbiased network analysis method based on large-scale similarities in gene expression, copy number aberration, and protein interactions may provide novel biological insights. Applying a novel measure of network robustness, Ollivier-Ricci Curvature, we examined patterns in the RNA-Seq gene expression and CNA data and how they relate to clinical outcomes. Hierarchical clustering using ORC differentiated high-risk subtypes with low progression free survival. Differential gene expression analysis defined 118 genes with significantly aberrant expression. These genes, while not previously associated with MM, were associated with DNA repair, apoptosis, and the immune system. Univariate analysis identified 8/118 to be prognostic genes; all associated with the immune system. A network topology analysis identified both hub and bridge genes which connect known genes of biological significance of MM. Taken together, gene interaction network analysis in MM uses a novel method of global assessment to demonstrate complex immune dysregulation associated with shorter survival., (© 2023. The Author(s).)

The advance of sequencing technologies has enabled a thorough molecular characterization of the genome in human cancers. To improve patient prognosis predictions and subsequent treatment strategies, it is imperative to develop advanced computational methods to analyze large-scale, high-dimensional genomic data. However, traditional machine learning methods face a challenge in handling the high-dimensional, low-sample size problem that is shown in most genomic data sets. To address this, our group has developed geometric network analysis techniques on multi-omics data in connection with prior biological knowledge derived from protein-protein interactions (PPIs) or pathways. Geometric features obtained from the genomic network, such as Ollivier-Ricci curvature and the invariant measure of the associated Markov chain, have been shown to be predictive of survival outcomes in various cancers. In this study, we propose a novel supervised deep learning method called geometric graph neural network (GGNN) that incorporates such geometric features into deep learning for enhanced predictive power and interpretability. More specifically, we utilize a state-of-the-art graph neural network with sparse connections between the hidden layers based on known biology of the PPI network and pathway information. Geometric features along with multi-omics data are then incorporated into the corresponding layers. The proposed approach utilizes a local-global principle in such a manner that highly predictive features are selected at the front layers and fed directly to the last layer for multivariable Cox proportional-hazards regression modeling. The method was applied to multi-omics data from the CoMMpass study of multiple myeloma and ten major cancers in The Cancer Genome Atlas (TCGA). In most experiments, our method showed superior predictive performance compared to other alternative methods., Competing Interests: Declaration of Competing Interest Larry Norton discloses the following relationships and financial interests:, (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Background and Purpose: Late radiation-induced hematuria can develop in prostate cancer patients undergoing radiotherapy and can negatively impact the quality-of-life of survivors. If a genetic component of risk could be modeled, this could potentially be the basis for modifying treatment for high-risk patients. We therefore investigated whether a previously developed machine learning-based modeling method using genome-wide common single nucleotide polymorphisms (SNPs) can stratify patients in terms of the risk of radiation-induced hematuria., Materials and Methods: We applied a two-step machine learning algorithm that we previously developed for genome-wide association studies called pre-conditioned random forest regression (PRFR). PRFR includes a pre-conditioning step, producing adjusted outcomes, followed by random forest regression modeling. Data was from germline genome-wide SNPs for 668 prostate cancer patients treated with radiotherapy. The cohort was stratified only once, at the outset of the modeling process, into two groups: a training set (2/3 of samples) for modeling and a validation set (1/3 of samples). Post-modeling bioinformatics analysis was conducted to identify biological correlates plausibly associated with the risk of hematuria., Results: The PRFR method achieved significantly better predictive performance compared to other alternative methods (all p < 0.05). The odds ratio between the high and low risk groups, each of which consisted of 1/3 of samples in the validation set, was 2.87 (p = 0.029), implying a clinically useful level of discrimination. Bioinformatics analysis identified six key proteins encoded by CTNND2, GSK3B, KCNQ2, NEDD4L, PRKAA1, and TXNL1 genes as well as four statistically significant biological process networks previously shown to be associated with the bladder and urinary tract., Conclusion: The risk of hematuria is significantly dependent on common genetic variants. The PRFR algorithm resulted in a stratification of prostate cancer patients at differential risk levels of post-radiotherapy hematuria. Bioinformatics analysis identified important biological processes involved in radiation-induced hematuria., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Author(s). Published by Elsevier B.V. All rights reserved.)

Drug-induced liver injury (DILI) is an adverse hepatic drug reaction that can potentially lead to life-threatening liver failure. Previously published work in the scientific literature on DILI has provided valuable insights for the understanding of hepatotoxicity as well as drug development. However, the manual search of scientific literature in PubMed is laborious and time-consuming. Natural language processing (NLP) techniques along with artificial intelligence/machine learning approaches may allow for automatic processing in identifying DILI-related literature, but useful methods are yet to be demonstrated. To address this issue, we have developed an integrated NLP/machine learning classification model to identify DILI-related literature using only paper titles and abstracts. For prediction modeling, we used 14,203 publications provided by the Critical Assessment of Massive Data Analysis (CAMDA) challenge, employing word vectorization techniques in NLP in conjunction with machine learning methods. Classification modeling was performed using 2/3 of the data for training and the remainder for test in internal validation. The best performance was achieved using a linear support vector machine (SVM) model on the combined vectors derived from term frequency-inverse document frequency ( TF-IDF ) and Word2Vec , resulting in an accuracy of 95.0% and an F1-score of 95.0%. The final SVM model constructed from all 14,203 publications was tested on independent datasets, resulting in accuracies of 92.5%, 96.3%, and 98.3%, and F1-scores of 93.5%, 86.1%, and 75.6% for three test sets (T1-T3). Furthermore, the SVM model was tested on four external validation sets (V1-V4), resulting in accuracies of 92.0%, 96.2%, 98.3%, and 93.1%, and F1-scores of 92.4%, 82.9%, 75.0%, and 93.3%., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Oh, Tannenbaum and Deasy.)

Respiration can positively influence cerebrospinal fluid (CSF) flow in the brain, yet its effects on central nervous system (CNS) fluid homeostasis, including waste clearance function via glymphatic and meningeal lymphatic systems, remain unclear. Here, we investigated the effect of supporting respiratory function via continuous positive airway pressure (CPAP) on glymphatic-lymphatic function in spontaneously breathing anesthetized rodents. To do this, we used a systems approach combining engineering, MRI, computational fluid dynamics analysis, and physiological testing. We first designed a nasal CPAP device for use in the rat and demonstrated that it functioned similarly to clinical devices, as evidenced by its ability to open the upper airway, augment end-expiratory lung volume, and improve arterial oxygenation. We further showed that CPAP increased CSF flow speed at the skull base and augmented glymphatic transport regionally. The CPAP-induced augmented CSF flow speed was associated with an increase in intracranial pressure (ICP), including the ICP waveform pulse amplitude. We suggest that the augmented pulse amplitude with CPAP underlies the increase in CSF bulk flow and glymphatic transport. Our results provide insights into the functional crosstalk at the pulmonary-CSF interface and suggest that CPAP might have therapeutic benefit for sustaining glymphatic-lymphatic function.

We present a novel approach for volume exploration that is versatile yet effective in isolating semantic structures in both noisy and clean data. Specifically, we describe a hierarchical active contours approach based on Bhattacharyya gradient flow which is easier to control, robust to noise, and can incorporate various types of statistical information to drive an edge-agnostic exploration process. To facilitate a time-bound user-driven volume exploration process that is applicable to a wide variety of data sources, we present an efficient multi-GPU implementation that (1) is approximately 400 times faster than a single thread CPU implementation, (2) allows hierarchical exploration of 2D and 3D images, (3) supports customization through multidimensional attribute spaces, and (4) is applicable to a variety of data sources and semantic structures. The exploration system follows a 2-step process. It first applies active contours to isolate semantically meaningful subsets of the volume. It then applies transfer functions to the isolated regions locally to produce clear and clutter-free visualizations. We show the effectiveness of our approach in isolating and visualizing structures-of-interest without needing any specialized segmentation methods on a variety of data sources, including 3D optical microscopy, multi-channel optical volumes, abdominal and chest CT, micro-CT, MRI, simulation, and synthetic data. We also gathered feedback from a medical trainee regarding the usefulness of our approach and discussion on potential applications in clinical workflows.

Radiation exposure of healthy cells can halt cell cycle temporarily or permanently. In this work, we analyze the time evolution of p21 and p53 from two single cell datasets of retinal pigment epithelial cells exposed to several levels of radiation, and in particular, the effect of radiation on cell cycle arrest. Employing various quantification methods from signal processing, we show how p21 levels, and to a lesser extent p53 levels, dictate whether the cells are arrested in their cell cycle and how frequently these mitosis events are likely to occur. We observed that single cells exposed to the same dose of DNA damage exhibit heterogeneity in cellular outcomes and that the frequency of cell division is a more accurate monitor of cell damage rather than just radiation level. Finally, we show how heterogeneity in DNA damage signaling is manifested early in the response to radiation exposure level and has potential to predict long-term fate., (© 2022. The Author(s), under exclusive licence to ADMC Associazione Differenziamento e Morte Cellulare.)

The choroid plexus (ChP) is the blood-cerebrospinal fluid (CSF) barrier and the primary source of CSF. Acquired hydrocephalus, caused by brain infection or hemorrhage, lacks drug treatments due to obscure pathobiology. Our integrated, multi-omic investigation of post-infectious hydrocephalus (PIH) and post-hemorrhagic hydrocephalus (PHH) models revealed that lipopolysaccharide and blood breakdown products trigger highly similar TLR4-dependent immune responses at the ChP-CSF interface. The resulting CSF "cytokine storm", elicited from peripherally derived and border-associated ChP macrophages, causes increased CSF production from ChP epithelial cells via phospho-activation of the TNF-receptor-associated kinase SPAK, which serves as a regulatory scaffold of a multi-ion transporter protein complex. Genetic or pharmacological immunomodulation prevents PIH and PHH by antagonizing SPAK-dependent CSF hypersecretion. These results reveal the ChP as a dynamic, cellularly heterogeneous tissue with highly regulated immune-secretory capacity, expand our understanding of ChP immune-epithelial cell cross talk, and reframe PIH and PHH as related neuroimmune disorders vulnerable to small molecule pharmacotherapy., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 Elsevier Inc. All rights reserved.)

The regularized optimal mass transport (rOMT) problem adds a diffusion term to the continuity equation in the original dynamic formulation of the optimal mass transport (OMT) problem proposed by Benamou and Brenier. We show that the rOMT model serves as a powerful tool in computational fluid dynamics for visualizing fluid flows in the glymphatic system. In the present work, we describe how to modify the previous numerical method for efficient implementation, resulting in a significant reduction in computational runtime. Numerical results applied to synthetic and real-data are provided., Competing Interests: Declarations Conflict of interest The authors declare that they have no conflict of interest.

The emerging field of radiomics that transforms standard-of-care images to quantifiable scalar statistics endeavors to reveal the information hidden in these macroscopic images. The concept of texture is widely used and essential in many radiomic-based studies. Practice usually reduces spatial multidimensional texture matrices, e.g., gray-level co-occurrence matrices (GLCMs), to summary scalar features. These statistical features have been demonstrated to be strongly correlated and tend to contribute redundant information; and does not account for the spatial information hidden in the multivariate texture matrices. This study proposes a novel pipeline to deal with spatial texture features in radiomic studies. A new set of textural features that preserve the spatial information inherent in GLCMs is proposed and used for classification purposes. The set of the new features uses the Wasserstein metric from optimal mass transport theory (OMT) to quantify the spatial similarity between samples within a given label class. In particular, based on a selected subset of texture GLCMs from the training cohort, we propose new representative spatial texture features, which we incorporate into a supervised image classification pipeline. The pipeline relies on the support vector machine (SVM) algorithm along with Bayesian optimization and the Wasserstein metric. The selection of the best GLCM references is considered for each classification label and is performed during the training phase of the SVM classifier using a Bayesian optimizer. We assume that sample fitness is defined based on closeness (in the sense of the Wasserstein metric) and high correlation (Spearman's rank sense) with other samples in the same class. Moreover, the newly defined spatial texture features consist of the Wasserstein distance between the optimally selected references and the remaining samples. We assessed the performance of the proposed classification pipeline in diagnosing the coronavirus disease 2019 (COVID-19) from computed tomographic (CT) images. To evaluate the proposed spatial features' added value, we compared the performance of the proposed classification pipeline with other SVM-based classifiers that account for different texture features, namely: statistical features only, optimized spatial features using Euclidean metric, non-optimized spatial features with Wasserstein metric. The proposed technique, which accounts for the optimized spatial texture feature with Wasserstein metric, shows great potential in classifying new COVID CT images that the algorithm has not seen in the training step. The MATLAB code of the proposed classification pipeline is made available. It can be used to find the best reference samples in other data cohorts, which can then be employed to build different prediction models., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 The Authors. Published by Elsevier Ltd.. All rights reserved.)

This review focuses on the physiology of glymphatic solute transport and waste clearance, using evidence from experimental animal models as well as from human studies. Specific topics addressed include the biophysical characteristics of fluid and solute transport in the central nervous system, glymphatic-lymphatic coupling, as well as the role of cerebrospinal fluid movement for brain waste clearance. We also discuss the current understanding of mechanisms underlying increased waste clearance during sleep.

Background: Compositional systems, represented as parts of some whole, are ubiquitous. They encompass the abundances of proteins in a cell, the distribution of organisms in nature, and the stoichiometry of the most basic chemical reactions. Thus, a central goal is to understand how such processes emerge from the behaviors of their components and their pairwise interactions. Such a study, however, is challenging for two key reasons. Firstly, such systems are complex and depend, often stochastically, on their constituent parts. Secondly, the data lie on a simplex which influences their correlations., Results: To resolve both of these issues, we provide a general and data-driven modeling tool for compositional systems called Compositional Maximum Entropy (CME). By integrating the prior geometric structure of compositions with sample-specific information, CME infers the underlying multivariate relationships between the constituent components. We provide two proofs of principle. First, we measure the relative abundances of different bacteria and infer how they interact. Second, we show that our method outperforms a common alternative for the extraction of gene-gene interactions in triple-negative breast cancer., Conclusions: CME provides novel and biologically-intuitive insights and is promising as a comprehensive quantitative framework for compositional data., (© 2022. The Author(s).)

Diffusion tensor imaging (DTI) has been used as an outcome measure in clinical trials for several psychiatric disorders but has rarely been explored in autism clinical trials. This is despite a large body of research suggesting altered white matter structure in autistic individuals. The current study is a secondary analysis of changes in white matter connectivity from a double-blind placebo-control trial of a single intravenous cord blood infusion in 2-7-year-old autistic children (1). Both clinical assessments and DTI were collected at baseline and 6 months after infusion. This study used two measures of white matter connectivity: change in node-to-node connectivity as measured through DTI streamlines and a novel measure of feedback network connectivity, Ollivier-Ricci curvature (ORC). ORC is a network measure which considers both local and global connectivity to assess the robustness of any given pathway. Using both the streamline and ORC analyses, we found reorganization of white matter pathways in predominantly frontal and temporal brain networks in autistic children who received umbilical cord blood treatment versus those who received a placebo. By looking at changes in network robustness, this study examined not only the direct, physical changes in connectivity, but changes with respect to the whole brain network. Together, these results suggest the use of DTI and ORC should be further explored as a potential biomarker in future autism clinical trials. These results, however, should not be interpreted as evidence for the efficacy of cord blood for improving clinical outcomes in autism. This paper presents a secondary analysis using data from a clinical trial that was prospectively registered with ClinicalTrials.gov(NCT02847182)., Competing Interests: Authors KC, GS, and GD reported technology unrelated to the submitted work that has been licensed, have benefited financially from this license, and have a patent. Authors GD, AS, and JK had a patent and have developed technology, data, and/or products that have been licensed to Cryocell, Inc., from which they and Duke University have benefited financially. Author JK was the Director of the Carolinas Cord Blood Bank, Medical Director of Cryocell, Inc., and is a paid consultant for Neurogene. Author GS was affiliated with Apple Inc., the work here reported was initiated before such affiliation and it is independent of it. Allen Song has patents licensed unrelated to the submitted work and receives grants from GE Healthcare unrelated to the submitted work. Author GD was on the Scientific Advisory Boards of Akili Interactive, Inc., Zynerba, Nonverbal Learning Disability Project, and Tris Pharma, is a consultant to Apple, Gerson Lehrman Group, and Guidepoint Global, Inc., and receives book royalties from Guilford Press and Springer Nature. Author GD has stock interests in Neuvana, Inc. Author AT was a consultant for Polaris. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Simhal, Carpenter, Kurtzberg, Song, Tannenbaum, Zhang, Sapiro and Dawson.)

We review theoretical and numerical models of the glymphatic system, which circulates cerebrospinal fluid and interstitial fluid around the brain, facilitating solute transport. Models enable hypothesis development and predictions of transport, with clinical applications including drug delivery, stroke, cardiac arrest, and neurodegenerative disorders like Alzheimer's disease. We sort existing models into broad categories by anatomical function: Perivascular flow, transport in brain parenchyma, interfaces to perivascular spaces, efflux routes, and links to neuronal activity. Needs and opportunities for future work are highlighted wherever possible; new models, expanded models, and novel experiments to inform models could all have tremendous value for advancing the field., Competing Interests: We declare no conflicts of interest., (© 2022 The Author(s).)

The remarkable growth of multi-platform genomic profiles has led to the challenge of multiomics data integration. In this study, we present a novel network-based multiomics clustering founded on the Wasserstein distance from optimal mass transport. This distance has many important geometric properties making it a suitable choice for application in machine learning and clustering. Our proposed method of aggregating multiomics and Wasserstein distance clustering (aWCluster) is applied to breast carcinoma as well as bladder carcinoma, colorectal adenocarcinoma, renal carcinoma, lung non-small cell adenocarcinoma, and endometrial carcinoma from The Cancer Genome Atlas project. Subtypes were characterized by the concordant effect of mRNA expression, DNA copy number alteration, and DNA methylation of genes and their neighbors in the interaction network. aWCluster successfully clusters all cancer types into classes with significantly different survival rates. Also, a gene ontology enrichment analysis of significant genes in the low survival subgroup of breast cancer leads to the well-known phenomenon of tumor hypoxia and the transcription factor ETS1 whose expression is induced by hypoxia. We believe aWCluster has the potential to discover novel subtypes and biomarkers by accentuating the genes that have concordant multiomics measurements in their interaction network, which are challenging to find without the network inference or with single omics analysis.

In this paper, we present a network-based clustering method, called vector Wasserstein clustering (vWCluster), based on the vector-valued Wasserstein distance derived from optimal mass transport (OMT) theory. This approach allows for the natural integration of multi-layer representations of data in a given network from which one derives clusters via a hierarchical clustering approach. In this study, we applied the methodology to multi-omics data from the two largest breast cancer studies. The resultant clusters showed significantly different survival rates in Kaplan-Meier analysis in both datasets. CIBERSORT scores were compared among the identified clusters. Out of the 22 CIBERSORT immune cell types, 9 were commonly significantly different in both datasets, suggesting the difference of tumor immune microenvironment in the clusters. vWCluster can aggregate multi-omics data represented as a vectorial form in a network with multiple layers, taking into account the concordant effect of heterogeneous data, and further identify subgroups of tumors in terms of mortality., Competing Interests: The authors have declared that no competing interests exist.

The development of reliable predictive models for individual cancer cell lines to identify an optimal cancer drug is a crucial step to accelerate personalized medicine, but vast differences in cancer cell lines and drug characteristics make it quite challenging to develop predictive models that result in high predictive power and explain the similarity of cell lines or drugs. Our study proposes a novel network-based methodology that breaks the problem into smaller, more interpretable problems to improve the predictive power of anti-cancer drug responses in cell lines. For the drug-sensitivity study, we used the GDSC database for 915 cell lines and 200 drugs. The theory of optimal mass transport was first used to separately cluster cell lines and drugs, using gene-expression profiles and extensive cheminformatic drug features, represented in a form of data networks. To predict cell-line specific drug responses, random forest regression modeling was separately performed for each cell-line drug cluster pair. Post-modeling biological analysis was further performed to identify potential biological correlates associated with drug responses. The network-based clustering method resulted in 30 distinct cell-line drug cluster pairs. Predictive modeling on each cell-line-drug cluster outperformed alternative computational methods in predicting drug responses. We found that among the four drugs top-ranked with respect to prediction performance, three targeted the PI3K/mTOR signaling pathway. Predictive modeling on clustered subsets of cell lines and drugs improved the prediction accuracy of cell-line specific drug responses. Post-modeling analysis identified plausible biological processes associated with drug responses.

Network analysis methods can potentially quantify cancer aberrations in gene networks without introducing fitted parameters or variable selection. A new network curvature-based method is introduced to provide an integrated measure of variability within cancer gene networks. The method is applied to high-grade serous ovarian cancers (HGSOCs) to predict response to immune checkpoint inhibitors (ICIs) and to rank key genes associated with prognosis. Copy number alterations (CNAs) from targeted and whole-exome sequencing data were extracted for HGSOC patients (n = 45) treated with ICIs. CNAs at a gene level were represented on a protein-protein interaction network to define patient-specific networks with a fixed topology. A version of Ollivier-Ricci curvature was used to identify genes that play a potentially key role in response to immunotherapy and further to stratify patients at high risk of mortality. Overall survival (OS) was defined as the time from the start of ICI treatment to either death or last follow-up. Kaplan-Meier analysis with log-rank test was performed to assess OS between the high and low curvature classified groups. The network curvature analysis stratified patients at high risk of mortality with p = 0.00047 in Kaplan-Meier analysis in HGSOC patients receiving ICI. Genes with high curvature were in accordance with CNAs relevant to ovarian cancer. Network curvature using CNAs has the potential to be a novel predictor for OS in HGSOC patients treated with immunotherapy., (© 2021. The Author(s).)

Over the past decade there has been an enormous progress in our understanding of fluid and solute transport in the central nervous system (CNS). This is due to a number of factors, including important developments in whole brain imaging technology and computational fluid dynamics analysis employed for the elucidation of glymphatic transport function in the live animal and human brain. In this paper, we review the technical aspects of dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) in combination with administration of Gd-based tracers into the cerebrospinal fluid (CSF) for tracking glymphatic solute and fluid transport in the CNS as well as lymphatic drainage. Used in conjunction with advanced computational processing methods including optimal mass transport analysis, one gains new insights into the biophysical forces governing solute transport in the CNS which leads to intriguing new research directions. Considering drainage pathways, we review the novel T1 mapping technique for quantifying glymphatic transport and cervical lymph node drainage concurrently in the same subject. We provide an overview of knowledge gleaned from DCE-MRI studies of glymphatic transport and meningeal lymphatic drainage. Finally, we introduce positron emission tomography (PET) and CSF administration of radiotracers as an alternative method to explore other pharmacokinetic aspects of CSF transport into brain parenchyma as well as efflux pathways., (Copyright © 2020 IBRO. All rights reserved.)

Motivation: Convolutional neural networks (CNNs) have achieved great success in the areas of image processing and computer vision, handling grid-structured inputs and efficiently capturing local dependencies through multiple levels of abstraction. However, a lack of interpretability remains a key barrier to the adoption of deep neural networks, particularly in predictive modeling of disease outcomes. Moreover, because biological array data are generally represented in a non-grid structured format, CNNs cannot be applied directly., Results: To address these issues, we propose a novel method, called PathCNN, that constructs an interpretable CNN model on integrated multi-omics data using a newly defined pathway image. PathCNN showed promising predictive performance in differentiating between long-term survival (LTS) and non-LTS when applied to glioblastoma multiforme (GBM). The adoption of a visualization tool coupled with statistical analysis enabled the identification of plausible pathways associated with survival in GBM. In summary, PathCNN demonstrates that CNNs can be effectively applied to multi-omics data in an interpretable manner, resulting in promising predictive power while identifying key biological correlates of disease., Availability and Implementation: The source code is freely available at: https://github.com/mskspi/PathCNN., (© The Author(s) 2021. Published by Oxford University Press.)

Purpose: The goal of this study is to develop innovative methods for identifying radiomic features that are reproducible over varying image acquisition settings. Approach: We propose a regularized partial correlation network to identify reliable and reproducible radiomic features. This approach was tested on two radiomic feature sets generated using two different reconstruction methods on computed tomography (CT) scans from a cohort of 47 lung cancer patients. The largest common network component between the two networks was tested on phantom data consisting of five cancer samples. To further investigate whether radiomic features found can identify phenotypes, we propose a k -means clustering algorithm coupled with the optimal mass transport theory. This approach following the regularized partial correlation network analysis was tested on CT scans from 77 head and neck squamous cell carcinoma (HNSCC) patients in the Cancer Imaging Archive (TCIA) and validated using an independent dataset. Results: A set of common radiomic features was found in relatively large network components between the resultant two partial correlation networks resulting from a cohort of lung cancer patients. The reliability and reproducibility of those radiomic features were further validated on phantom data using the Wasserstein distance. Further analysis using the network-based Wasserstein k -means algorithm on the TCIA HNSCC data showed that the resulting clusters separate tumor subsites as well as HPV status, and this was validated on an independent dataset. Conclusion: We showed that a network-based analysis enables identifying reproducible radiomic features and use of the selected set of features can enhance clustering results., (© 2021 The Authors.)

Spiculations are important predictors of lung cancer malignancy, which are spikes on the surface of the pulmonary nodules. In this study, we proposed an interpretable and parameter-free technique to quantify the spiculation using area distortion metric obtained by the conformal (angle-preserving) spherical parameterization. We exploit the insight that for an angle-preserved spherical mapping of a given nodule, the corresponding negative area distortion precisely characterizes the spiculations on that nodule. We introduced novel spiculation scores based on the area distortion metric and spiculation measures. We also semi-automatically segment lung nodule (for reproducibility) as well as vessel and wall attachment to differentiate the real spiculations from lobulation and attachment. A simple pathological malignancy prediction model is also introduced. We used the publicly-available LIDC-IDRI dataset pathologists (strong-label) and radiologists (weak-label) ratings to train and test radiomics models containing this feature, and then externally validate the models. We achieved AUC = 0.80 and 0.76, respectively, with the models trained on the 811 weakly-labeled LIDC datasets and tested on the 72 strongly-labeled LIDC and 73 LUNGx datasets; the previous best model for LUNGx had AUC = 0.68. The number-of-spiculations feature was found to be highly correlated (Spearman's rank correlation coefficient ρ=0.44) with the radiologists' spiculation score. We developed a reproducible and interpretable, parameter-free technique for quantifying spiculations on nodules. The spiculation quantification measures was then applied to the radiomics framework for pathological malignancy prediction with reproducible semi-automatic segmentation of nodule. Using our interpretable features (size, attachment, spiculation, lobulation), we were able to achieve higher performance than previous models. In the future, we will exhaustively test our model for lung cancer screening in the clinic., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020 Elsevier B.V. All rights reserved.)

Background: Coronavirus disease 2019 (COVID-19) is a global public health concern. Recently, a genome-wide association study (GWAS) was performed with participants recruited from Italy and Spain by an international consortium group., Methods: Summary GWAS statistics for 1610 patients with COVID-19 respiratory failure and 2205 controls were downloaded. In the current study, we analyzed the summary statistics with the information of loci and p-values for 8,582,968 single-nucleotide polymorphisms (SNPs), using gene ontology analysis to determine the top biological processes implicated in respiratory failure in COVID-19 patients., Results: We considered the top 708 SNPs, using a p-value cutoff of 5 × 10 - 5 , which were mapped to the nearest genes, leading to 144 unique genes. The list of genes was input into a curated database to conduct gene ontology and protein-protein interaction (PPI) analyses. The top ranked biological processes were wound healing, epithelial structure maintenance, muscle system processes, and cardiac-relevant biological processes with a false discovery rate < 0.05. In the PPI analysis, the largest connected network consisted of 8 genes. Through a literature search, 7 out of the 8 gene products were found to be implicated in both pulmonary and cardiac diseases., Conclusion: Gene ontology and PPI analyses identified cardio-pulmonary processes that may partially explain the risk of respiratory failure in COVID-19 patients.

The brain's high bioenergetic state is paralleled by high metabolic waste production. Authentic lymphatic vasculature is lacking in brain parenchyma. Cerebrospinal fluid (CSF) flow has long been thought to facilitate central nervous system detoxification in place of lymphatics, but the exact processes involved in toxic waste clearance from the brain remain incompletely understood. Over the past 8 yr, novel data in animals and humans have begun to shed new light on these processes in the form of the "glymphatic system," a brain-wide perivascular transit passageway dedicated to CSF transport and interstitial fluid exchange that facilitates metabolic waste drainage from the brain. Here we will discuss glymphatic system anatomy and methods to visualize and quantify glymphatic system (GS) transport in the brain and also discuss physiological drivers of its function in normal brain and in neurodegeneration.

The purpose of this work is to make a case for epidemiological models with fractional exponent in the contribution of sub-populations to the incidence rate. More specifically, we question the standard assumption in the literature on epidemiological models, where the incidence rate dictating propagation of infections is taken to be proportional to the product between the infected and susceptible sub-populations; a model that relies on strong mixing between the two groups and widespread contact between members of the groups. We contend, that contact between infected and susceptible individuals, especially during the early phases of an epidemic, takes place over a (possibly diffused) boundary between the respective sub-populations. As a result, the rate of transmission depends on the product of fractional powers instead. The intuition relies on the fact that infection grows in geographically concentrated cells, in contrast to the standard product model that relies on complete mixing of the susceptible to infected sub-populations. We validate the hypothesis of fractional exponents (1) by numerical simulation for disease propagation in graphs imposing a local structure to allowed disease transmissions and (2) by fitting the model to the JHU CSSE COVID-19 Data for the period Jan-22-20 to April-30-20, for the countries of Italy, Germany, France, and Spain.

Variations in hematoxylin and eosin (H&E) stained images (due to clinical lab protocols, scanners, etc) directly impact the quality and accuracy of clinical diagnosis, and hence it is important to control for these variations for a reliable diagnosis. In this work, we present a new approach based on the multimarginal Wasserstein barycenter to normalize and augment H&E stained images given one or more references. Specifically, we provide a mathematically robust way of naturally incorporating additional images as intermediate references to drive stain normalization and augmentation simultaneously. The presented approach showed superior results quantitatively and qualitatively as compared to state-of-the-art methods for stain normalization. We further validated our stain normalization and augmentations in the nuclei segmentation task on a publicly available dataset, achieving state-of-the-art results against competing approaches.

In this work, a unified representation of all the time-varying dynamics is accomplished with a Lagrangian framework for analyzing Fisher-Rao regularized dynamical optimal mass transport (OMT) derived flows. While formally equivalent to the Eulerian based Schrödinger bridge OMT regularization scheme, the Fisher-Rao approach allows a simple and interpretable methodology for studying the flows of interest in the present work. The advantage of the proposed Lagrangian technique is that the time-varying particle trajectories and attributes are displayed in a single visualization. This provides a natural capability to identify and distinguish flows under different conditions. The Lagrangian analysis applied to the glymphatic system (brain waste removal pathway associated with Alzheimer's Disease) successfully captures known flows and distinguishes between flow patterns under two different anesthetics, providing deeper insights into altered states of waste drainage.

We present a technique to construct a simplification of a feature network which can be used for interactive data exploration, biological hypothesis generation, and the detection of communities or modules of cofunctional features. These are modules of features that are not necessarily correlated, but nevertheless exhibit common function in their network context as measured by similarity of relationships with neighboring features. In the case of genetic networks, traditional pathway analyses tend to assume that, ideally, all genes in a module exhibit very similar function, independent of relationships with other genes. The proposed technique explicitly relaxes this assumption by employing the comparison of relational profiles. For example, two genes which always activate a third gene are grouped together even if they never do so concurrently. They have common, but not identical, function. The comparison is driven by an average of a certain computationally efficient comparison metric between Gaussian mixture models. The method has its basis in the local connection structure of the network and the collection of joint distributions of the data associated with nodal neighborhoods. It is benchmarked on networks with known community structures. As the main application, we analyzed the gene regulatory network in lung adenocarcinoma, finding a cofunctional module of genes including the pregnancy-specific glycoproteins (PSGs). About 20% of patients with lung, breast, uterus, and colon cancer in The Cancer Genome Atlas (TCGA) have an elevated PSG+ signature, with associated poor group prognosis. In conjunction with previous results relating PSGs to tolerance in the immune system, these findings implicate the PSGs in a potential immune tolerance mechanism of cancers., Competing Interests: Competing interest statement: J.O.D. is a shareholder in PAIGE.AI. The authors have applied for patent protection related to this publication., (Copyright © 2020 the Author(s). Published by PNAS.)

Ollivier-Ricci curvature is a method for measuring the robustness of connections in a network. In this work, we use curvature to measure changes in robustness of brain networks in children with autism spectrum disorder (ASD). In an open label clinical trials, participants with ASD were administered a single infusion of autologous umbilical cord blood and, as part of their clinical outcome measures, were imaged with diffusion MRI before and after the infusion. By using Ricci curvature to measure changes in robustness, we quantified both local and global changes in the brain networks and their potential relationship with the infusion. Our results find changes in the curvature of the connections between regions associated with ASD that were not detected via traditional brain network analysis.

We consider damped stochastic systems in a controlled (time-varying) potential and study their transition between specified Gibbs-equilibria states in finite time. By the second law of thermodynamics, the minimum amount of work needed to transition from one equilibrium state to another is the difference between the Helmholtz free energy of the two states and can only be achieved by a reversible (infinitely slow) process. The minimal gap between the work needed in a finite-time transition and the work during a reversible one, turns out to equal the square of the optimal mass transport (Wasserstein-2) distance between the two end-point distributions times the inverse of the duration needed for the transition. This result, in fact, relates non-equilibrium optimal control strategies (protocols) to gradient flows of entropy functionals via the Jordan-Kinderlehrer-Otto scheme. The purpose of this paper is to introduce ideas and results from the emerging field of stochastic thermodynamics in the setting of classical regulator theory, and to draw connections and derive such fundamental relations from a control perspective in a multivariable setting.

We seek network routing towards a desired final distribution that can mediate possible random link failures. In other words, we seek a routing plan that utilizes alternative routes so as to be relatively robust to link failures. To this end, we provide a mathematical formulation of a relaxed transport problem where the final distribution only needs to be close to the desired one. The problem is cast as a maximum entropy problem for probability distributions on paths with an added terminal cost. The entropic regularizing penalty aims at distributing the choice of paths amongst possible alternatives. We prove that the unique solution may be obtained by solving a generalized Schrödinger system of equations. An iterative algorithm to compute the solution is provided. Each iteration of the algorithm contracts the distance (in the Hilbert metric) to the optimal solution by more than 1/2, leading to extremely fast convergence.

The Wasserstein distance is a powerful metric based on the theory of optimal mass transport. It gives a natural measure of the distance between two distributions with a wide range of applications. In contrast to a number of the common divergences on distributions such as Kullback-Leibler or Jensen-Shannon, it is (weakly) continuous, and thus ideal for analyzing corrupted and noisy data. Until recently, however, no kernel methods for dealing with nonlinear data have been proposed via the Wasserstein distance. In this work, we develop a novel method to compute the L 2 -Wasserstein distance in reproducing kernel Hilbert spaces (RKHS) called kernel L 2 -Wasserstein distance, which is implemented using the kernel trick. The latter is a general method in machine learning employed to handle data in a nonlinear manner. We evaluate the proposed approach in identifying computed tomography (CT) slices with dental artifacts in head and neck cancer, performing unsupervised hierarchical clustering on the resulting Wasserstein distance matrix that is computed on imaging texture features extracted from each CT slice. We further compare the performance of kernel Wasserstein distance with alternatives including kernel Kullback-Leibler divergence we previously developed. Our experiments show that the kernel approach outperforms classical non-kernel approaches in identifying CT slices with artifacts., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

Perivascular spaces include a variety of passageways around arterioles, capillaries and venules in the brain, along which a range of substances can move. Although perivascular spaces were first identified over 150 years ago, they have come to prominence recently owing to advances in knowledge of their roles in clearance of interstitial fluid and waste from the brain, particularly during sleep, and in the pathogenesis of small vessel disease, Alzheimer disease and other neurodegenerative and inflammatory disorders. Experimental advances have facilitated in vivo studies of perivascular space function in intact rodent models during wakefulness and sleep, and MRI in humans has enabled perivascular space morphology to be related to cognitive function, vascular risk factors, vascular and neurodegenerative brain lesions, sleep patterns and cerebral haemodynamics. Many questions about perivascular spaces remain, but what is now clear is that normal perivascular space function is important for maintaining brain health. Here, we review perivascular space anatomy, physiology and pathology, particularly as seen with MRI in humans, and consider translation from models to humans to highlight knowns, unknowns, controversies and clinical relevance.

The glymphatic system (GS) hypothesis states that advective driven cerebrospinal fluid (CSF) influx from the perivascular spaces into the interstitial fluid space rapidly transport solutes and clear waste from brain. However, the presence of advection in neuropil is contested and solutes are claimed to be transported by diffusion only. To address this controversy, we implemented a regularized version of the optimal mass transport (rOMT) problem, wherein the advection/diffusion equation is the only a priori assumption required. rOMT analysis with a Lagrangian perspective of GS transport revealed that solute speed was faster in CSF compared to grey and white matter. Further, rOMT analysis also demonstrated 2-fold differences in regional solute speed within the brain. Collectively, these results imply that advective transport dominates in CSF while diffusion and advection both contribute to GS transport in parenchyma. In a rat model of cerebral small vessel disease (cSVD), solute transport in the perivascular spaces (PVS) and PVS-to-tissue transfer was slower compared to normal rats. Thus, the analytical framework of rOMT provides novel insights in the local dynamics of GS transport that may have implications for neurodegenerative diseases. Future studies should apply the rOMT analysis approach to confirm GS transport reductions in humans with cSVD.

Purpose: Conventional radiotherapy for large lung tumors is given over several weeks, during which the tumor typically regresses in a highly nonuniform and variable manner. Adaptive radiotherapy would ideally follow these shape changes, but we need an accurate method to extrapolate tumor shape changes. We propose a computationally efficient algorithm to quantitate tumor surface shape changes that makes minimal assumptions, identifies fixed points, and can be used to predict future tumor geometrical response., Methods: A novel combination of nonrigid iterative closest point (ICP) and local shape-preserving map algorithms, LDeform, is developed to enable visualization, prediction, and categorization of both diffeomorphic and nondiffeomorphic tumor deformations during an extended course of radiotherapy., Results: We tested and validated our technique on 31 longitudinal CT/MRI subjects, with five to nine time points each. Based on this tumor deformation analysis, regions of local growth, shrinkage, and anchoring are identified and tracked across multiple time points. This categorization in turn represents a rational biomarker of local response. Results demonstrate useful predictive power, with an averaged Dice coefficient and surface mean-squared error of 0.85 and 2.8 mm, respectively, over all images., Conclusions: We conclude that the LDeform algorithm can facilitate the adaptive decision-making process during lung cancer radiotherapy., (© 2019 American Association of Physicists in Medicine.)

Unbalanced optimal mass transport (OMT) seeks to remove the conservation of mass constraint by adding a source term to the standard continuity equation in the Benamou-Brenier formulation of OMT. In this study, we show how the unbalanced case fits into the vector-valued OMT framework simply by adding an auxiliary source layer and taking the flow between the source layer and the original layer(s) as the source term. This allows for unbalanced models both in the scalar and vector-valued density settings. The results are demonstrated on a number of synthetic and real vector-valued data sets.