Cardinal is an R package for statistical analysis of mass spectrometry-based imaging (MSI) experiments of biological samples such as tissues. Cardinal supports both Matrix-Assisted Laser Desorption/Ionization (MALDI) and Desorption Electrospray Ionization-based MSI workflows, and experiments with multiple tissues and complex designs. The main analytical functionalities include (1) image segmentation, which partitions a tissue into regions of homogeneous chemical composition, selects the number of segments and the subset of informative ions, and characterizes the associated uncertainty and (2) image classification, which assigns locations on the tissue to pre-defined classes, selects the subset of informative ions, and estimates the resulting classification error by (cross-) validation. The statistical methods are based on mixture modeling and regularization.

Active neurons exert a mitogenic effect on normal neural precursor and oligodendroglial precursor cells, the putative cellular origins of high-grade glioma (HGG). By using optogenetic control of cortical neuronal activity in a patient-derived pediatric glioblastoma xenograft model, we demonstrate that active neurons similarly promote HGG proliferation and growth in vivo. Conditioned medium from optogenetically stimulated cortical slices promoted proliferation of pediatric and adult patient-derived HGG cultures, indicating secretion of activity-regulated mitogen(s). The synaptic protein neuroligin-3 (NLGN3) was identified as the leading candidate mitogen, and soluble NLGN3 was sufficient and necessary to promote robust HGG cell proliferation. NLGN3 induced PI3K-mTOR pathway activity and feedforward expression of NLGN3 in glioma cells. NLGN3 expression levels in human HGG negatively correlated with patient overall survival. These findings indicate the important role of active neurons in the brain tumor microenvironment and identify secreted NLGN3 as an unexpected mechanism promoting neuronal activity-regulated cancer growth.

The genetic origins of chemotherapy resistance are well established; however, the role of epigenetics in drug resistance is less well understood. To investigate mechanisms of drug resistance, we performed systematic genetic, epigenetic, and transcriptomic analyses of an alkylating agent-sensitive murine lymphoma cell line and a series of resistant lines derived by drug dose escalation.

The acquisition of oncogenic drivers is a critical feature of cancer progression. For some carcinomas, it is clear that certain genetic drivers occur early in neoplasia and others late. Why these drivers are selected and how these changes alter the neoplasia's fitness is less understood.

Nasopharyngeal carcinoma is a multifactorial disease mainly affecting the Asian and North African populations including Morocco. This study aimed to determine the epidemiological profile of nasopharyngeal carcinoma in Northern Morocco as well as its clinicopathological, therapeutic, and prognostic characteristics.

Assessing the efficacy of cancer therapeutics in mouse models is a critical step in treatment development. However, low-resolution measurement tools and small sample sizes make determining drug efficacy in vivo a difficult and time-intensive task. Here, we present a commercially scalable wearable electronic strain sensor that automates the in vivo testing of cancer therapeutics by continuously monitoring the micrometer-scale progression or regression of subcutaneously implanted tumors at the minute time scale. In two in vivo cancer mouse models, our sensor discerned differences in tumor volume dynamics between drug- and vehicle-treated tumors within 5 hours following therapy initiation. These short-term regression measurements were validated through histology, and caliper and bioluminescence measurements taken over weeklong treatment periods demonstrated the correlation with longer-term treatment response. We anticipate that real-time tumor regression datasets could help expedite and automate the process of screening cancer therapies in vivo.

Biomarkers remain the highest value proposition in cancer medicine today-especially protein biomarkers. Despite decades of evolving regulatory frameworks to facilitate the review of emerging technologies, biomarkers have been mostly about promise with very little to show for improvements in human health. Cancer is an emergent property of a complex system, and deconvoluting the integrative and dynamic nature of the overall system through biomarkers is a daunting proposition. The last 2 decades have seen an explosion of multiomics profiling and a range of advanced technologies for precision medicine, including the emergence of liquid biopsy, exciting advances in single-cell analysis, artificial intelligence (machine and deep learning) for data analysis, and many other advanced technologies that promise to transform biomarker discovery. Combining multiple omics modalities to acquire a more comprehensive landscape of the disease state, we are increasingly developing biomarkers to support therapy selection and patient monitoring. Furthering precision medicine, especially in oncology, necessitates moving away from the lens of reductionist thinking toward viewing and understanding that complex diseases are, in fact, complex adaptive systems. As such, we believe it is necessary to redefine biomarkers as representations of biological system states at different hierarchical levels of biological order. This definition could include traditional molecular, histologic, radiographic, or physiological characteristics, as well as emerging classes of digital markers and complex algorithms. To succeed in the future, we must move past purely observational individual studies and instead start building a mechanistic framework to enable integrative analysis of new studies within the context of prior studies. Identifying information in complex systems and applying theoretical constructs, such as information theory, to study cancer as a disease of dysregulated communication could prove to be "game changing" for the clinical outcome of cancer patients.

Glioblastoma multiforme (GBM) is an aggressive, malignant cancer Johnson and O'Neill (J Neurooncol 107: 359-364, 2012). An extract from the winter cherry plant (Withania somnifera ), AshwaMAX, is concentrated (4.3 %) for Withaferin A; a steroidal lactone that inhibits cancer cells Vanden Berghe et al. (Cancer Epidemiol Biomark Prev 23: 1985-1996, 2014). We hypothesized that AshwaMAX could treat GBM and that bioluminescence imaging (BLI) could track oral therapy in orthotopic murine models of glioblastoma. Human parietal-cortical glioblastoma cells (GBM2, GBM39) were isolated from primary tumors while U87-MG was obtained commercially. GBM2 was transduced with lentiviral vectors that express Green Fluorescent Protein (GFP)/firefly luciferase fusion proteins. Mutational, expression and proliferative status of GBMs were studied. Intracranial xenografts of glioblastomas were grown in the right frontal regions of female, nude mice (n = 3-5 per experiment). Tumor growth was followed through BLI. Neurosphere cultures (U87-MG, GBM2 and GBM39) were inhibited by AshwaMAX at IC50 of 1.4, 0.19 and 0.22 µM equivalent respectively and by Withaferin A with IC50 of 0.31, 0.28 and 0.25 µM respectively. Oral gavage, every other day, of AshwaMAX (40 mg/kg per day) significantly reduced bioluminescence signal (n = 3 mice, p < 0.02, four parameter non-linear regression analysis) in preclinical models. After 30 days of treatment, bioluminescent signal increased suggesting onset of resistance. BLI signal for control, vehicle-treated mice increased and then plateaued. Bioluminescent imaging revealed diffuse growth of GBM2 xenografts. With AshwaMAX, GBM neurospheres collapsed at nanomolar concentrations. Oral treatment studies on murine models confirmed that AshwaMAX is effective against orthotopic GBM. AshwaMAX is thus a promising candidate for future clinical translation in patients with GBM.

The genome of Mycobacterium tuberculosis was analyzed using recently developed computational approaches to infer protein function and protein linkages. We evaluated and employed a method to infer genes likely to belong to the same operon, as judged by the nucleotide distance between genes in the same genomic orientation, and combined this method with those of the Rosetta Stone, Phylogenetic Profile and conserved Gene Neighbor computational methods for the inference of protein function.

We present the Saccharomyces cerevisiae PeptideAtlas composed from 47 diverse experiments and 4.9 million tandem mass spectra. The observed peptides align to 61% of Saccharomyces Genome Database (SGD) open reading frames (ORFs), 49% of the uncharacterized SGD ORFs, 54% of S. cerevisiae ORFs with a Gene Ontology annotation of 'molecular function unknown', and 76% of ORFs with Gene names. We highlight the use of this resource for data mining, construction of high quality lists for targeted proteomics, validation of proteins, and software development.

Convolutional neural networks (CNNs) have gained steady popularity as a tool to perform automatic classification of whole slide histology images. While CNNs have proven to be powerful classifiers in this context, they fail to explain this classification, as the network engineered features used for modeling and classification are ONLY interpretable by the CNNs themselves. This work aims at enhancing a traditional neural network model to perform histology image modeling, patient classification, and interpretation of the distinctive features identified by the network within the histology whole slide images (WSIs). We synthesize a workflow which (a) intelligently samples the training data by automatically selecting only image areas that display visible disease-relevant tissue state and (b) isolates regions most pertinent to the trained CNN prediction and translates them to observable and qualitative features such as color, intensity, cell and tissue morphology and texture. We use the Cancer Genome Atlas's Breast Invasive Carcinoma (TCGA-BRCA) histology dataset to build a model predicting patient attributes (disease stage and node status) and the tumor proliferation challenge (TUPAC 2016) breast cancer histology image repository to help identify disease-relevant tissue state (mitotic activity). We find that our enhanced CNN based workflow both increased patient attribute predictive accuracy (~2% increase for disease stage and ~10% increase for node status) and experimentally proved that a data-driven CNN histology model predicting breast invasive carcinoma stages is highly sensitive to features such as color, cell size, and shape, granularity, and uniformity. This work summarizes the need for understanding the widely trusted models built using deep learning and adds a layer of biological context to a technique that functioned as a classification only approach till now.

While technological advances have made it possible to profile the immune system at high resolution, translating high-throughput data into knowledge of immune mechanisms has been challenged by the complexity of the interactions underlying immune processes. Tools to explore the immune network are critical for better understanding the multi-layered processes that underlie immune function and dysfunction, but require a standardized network map of immune interactions. To facilitate this we have developed ImmunoGlobe, a manually curated intercellular immune interaction network extracted from Janeway's Immunobiology textbook.

Many initially successful anticancer therapies lose effectiveness over time, and eventually, cancer cells acquire resistance to the therapy. Acquired resistance remains a major obstacle to improving remission rates and achieving prolonged disease-free survival. Consequently, novel approaches to overcome or prevent resistance are of significant clinical importance. There has been considerable interest in treating non-small cell lung cancer (NSCLC) with combinations of EGFR-targeted therapeutics (e.g., erlotinib) and cytotoxic therapeutics (e.g., paclitaxel); however, acquired resistance to erlotinib, driven by a variety of mechanisms, remains an obstacle to treatment success. In about 50% of cases, resistance is due to a T790M point mutation in EGFR, and T790M-containing cells ultimately dominate the tumor composition and lead to tumor regrowth. We employed a combined experimental and mathematical modeling-based approach to identify treatment strategies that impede the outgrowth of primary T790M-mediated resistance in NSCLC populations. Our mathematical model predicts the population dynamics of mixtures of sensitive and resistant cells, thereby describing how the tumor composition, initial fraction of resistant cells, and degree of selective pressure influence the time until progression of disease. Model development relied upon quantitative experimental measurements of cell proliferation and death using a novel microscopy approach. Using this approach, we systematically explored the space of combination treatment strategies and demonstrated that optimally timed sequential strategies yielded large improvements in survival outcome relative to monotherapies at the same concentrations. Our investigations revealed regions of the treatment space in which low-dose sequential combination strategies, after preclinical validation, may lead to a tumor reduction and improved survival outcome for patients with T790M-mediated resistance.

Tumor-derived, circulating proteins are potentially useful as biomarkers for detection of cancer, for monitoring of disease progression, regression and recurrence, and for assessment of therapeutic response. Here we interrogated how a protein's stability, cellular localization, and abundance affect its observability in blood by mass-spectrometry-based proteomics techniques. We performed proteomic profiling on tumors and plasma from two different xenograft mouse models. A statistical analysis of this data revealed protein properties indicative of the detection level in plasma. Though 20% of the proteins identified in plasma were tumor-derived, only 5% of the proteins observed in the tumor tissue were found in plasma. Both intracellular and extracellular tumor proteins were observed in plasma; however, after normalizing for tumor abundance, extracellular proteins were seven times more likely to be detected. Although proteins that were more abundant in the tumor were also more likely to be observed in plasma, the relationship was nonlinear: Doubling the spectral count increased detection rate by only 50%. Many secreted proteins, even those with relatively low spectral count, were observed in plasma, but few low abundance intracellular proteins were observed. Proteins predicted to be stable by dipeptide composition were significantly more likely to be identified in plasma than less stable proteins. The number of tryptic peptides in a protein was not significantly related to the chance of a protein being observed in plasma. Quantitative comparison of large versus small tumors revealed that the abundance of proteins in plasma as measured by spectral count was associated with the tumor size, but the relationship was not one-to-one; a 3-fold decrease in tumor size resulted in a 16-fold decrease in protein abundance in plasma. This study provides quantitative support for a tumor-derived marker prioritization strategy that favors secreted and stable proteins over all but the most abundant intracellular proteins.

Mass spectrometry imaging is a powerful tool for investigating the spatial distribution of chemical compounds in a biological sample such as tissue. Two common goals of these experiments are unsupervised segmentation of images into newly discovered homogeneous segments and supervised classification of images into predefined classes. In both cases, the important secondary goals are to characterize the uncertainty associated with the segmentation and with the classification and to characterize the spectral features that define each segment or class. Recent analysis methods have focused on the spatial structure of the data to improve results. However, they either do not address these secondary goals or do this with separate post hoc procedures.We introduce spatial shrunken centroids, a statistical model-based framework for both supervised classification and unsupervised segmentation. It takes as input sets of previously detected, aligned, quantified, and normalized spectral features and expresses both spatial and multivariate nature of the data using probabilistic modeling. It selects informative subsets of spectral features that define each unsupervised segment or supervised class and quantifies and visualizes the uncertainty in spatial segmentations and in tissue classification. In the unsupervised setting, it also guides the choice of an appropriate number of segments. We demonstrate the usefulness of this framework in a supervised human renal cell carcinoma experimental dataset and several unsupervised experimental datasets, including a pig fetus cross-section, three rodent brains, and a controlled image with known ground truth. This framework is available for use within the open-source R package Cardinal as part of a full pipeline for the processing, visualization, and statistical analysis of mass spectrometry imaging experiments.

A major goal of proteomics research is the accurate and sensitive identification and quantification of a broad range of proteins within a sample. Data-independent acquisition (DIA) approaches that acquire MS/MS spectra independently of precursor information have been developed to overcome the reproducibility challenges of data-dependent acquisition and the limited breadth of targeted proteomics strategies. Typical DIA implementations use wide MS/MS isolation windows to acquire comprehensive fragment ion data. However, wide isolation windows produce highly chimeric spectra, limiting the achievable sensitivity and accuracy of quantification and identification. Here, we present a DIA strategy in which spectra are collected with overlapping (rather than adjacent or random) windows and then computationally demultiplexed. This approach improves precursor selectivity by nearly a factor of 2, without incurring any loss in mass range, mass resolution, chromatographic resolution, scan speed, or other key acquisition parameters. We demonstrate a 64% improvement in sensitivity and a 17% improvement in peptides detected in a 6-protein bovine mix spiked into a yeast background. To confirm the method's applicability to a realistic biological experiment, we also analyze the regulation of the proteasome in yeast grown in rapamycin and show that DIA experiments with overlapping windows can help elucidate its adaptation toward the degradation of oxidatively damaged proteins. Our integrated computational and experimental DIA strategy is compatible with any DIA-capable instrument. The computational demultiplexing algorithm required to analyze the data has been made available as part of the open-source proteomics software tools Skyline and msconvert (Proteowizard), making it easy to apply as part of standard proteomics workflows. Graphical Abstract.

Live cell imaging has improved our ability to measure phenotypic heterogeneity. However, bottlenecks in imaging and image processing often make it difficult to differentiate interesting biological behavior from technical artifact. Thus there is a need for new methods that improve data quality without sacrificing throughput. Here we present a 3-step workflow to improve dynamic phenotype measurements of heterogeneous cell populations. We provide guidelines for image acquisition, phenotype tracking, and data filtering to remove erroneous cell tracks using the novel Tracking Aberration Measure (TrAM). Our workflow is broadly applicable across imaging platforms and analysis software. By applying this workflow to cancer cell assays, we reduced aberrant cell track prevalence from 17% to 2%. The cost of this improvement was removing 15% of the well-tracked cells. This enabled detection of significant motility differences between cell lines. Similarly, we avoided detecting a false change in translocation kinetics by eliminating the true cause: varied proportions of unresponsive cells. Finally, by systematically seeking heterogeneous behaviors, we detected subpopulations that otherwise could have been missed, including early apoptotic events and pre-mitotic cells. We provide optimized protocols for specific applications and step-by-step guidelines for adapting them to a variety of biological systems.

Currently prostate-specific antigen is used for prostate cancer (PCa) screening, however it lacks the necessary specificity for differentiating PCa from other diseases of the prostate such as benign prostatic hyperplasia (BPH), presenting a clinical need to distinguish these cases at the molecular level. Protein glycosylation plays an important role in a number of cellular processes involved in neoplastic progression and is aberrant in PCa. In this study, we systematically interrogate the alterations in the circulating levels of hundreds of serum proteins and their glycoforms in PCa and BPH samples using multi-lectin affinity chromatography and quantitative mass spectrometry-based proteomics. Specific lectins (AAL, PHA-L and PHA-E) were used to target and chromatographically separate core-fucosylated and highly-branched protein glycoforms for analysis, as differential expression of these glycan types have been previously associated with PCa. Global levels of CD5L, CFP, C8A, BST1, and C7 were significantly increased in the PCa samples. Notable glycoform-specific alterations between BPH and PCa were identified among proteins CD163, C4A, and ATRN in the PHA-L/E fraction and among C4BPB and AZGP1 glycoforms in the AAL fraction. Despite these modest differences, substantial similarities in glycoproteomic profiles were observed between PCa and BPH sera.

To investigate the transition from non-cancerous to metastatic from a physical sciences perspective, the Physical Sciences-Oncology Centers (PS-OC) Network performed molecular and biophysical comparative studies of the non-tumorigenic MCF-10A and metastatic MDA-MB-231 breast epithelial cell lines, commonly used as models of cancer metastasis. Experiments were performed in 20 laboratories from 12 PS-OCs. Each laboratory was supplied with identical aliquots and common reagents and culture protocols. Analyses of these measurements revealed dramatic differences in their mechanics, migration, adhesion, oxygen response, and proteomic profiles. Model-based multi-omics approaches identified key differences between these cells' regulatory networks involved in morphology and survival. These results provide a multifaceted description of cellular parameters of two widely used cell lines and demonstrate the value of the PS-OC Network approach for integration of diverse experimental observations to elucidate the phenotypes associated with cancer metastasis.

Despite tremendous advances in targeted therapies against lung adenocarcinoma, the majority of patients do not benefit from personalized treatments. A deeper understanding of potential therapeutic targets is crucial to increase the survival of patients. One promising target, ADAR, is amplified in 13% of lung adenocarcinomas and in-vitro studies have demonstrated the potential of its therapeutic inhibition to inhibit tumor growth. ADAR edits millions of adenosines to inosines within the transcriptome, and while previous studies of ADAR in cancer have solely focused on protein-coding edits, >99% of edits occur in non-protein coding regions. Here, we develop a pipeline to discover the regulatory potential of RNA editing sites across the entire transcriptome and apply it to lung adenocarcinoma tumors from The Cancer Genome Atlas. This method predicts that 1413 genes contain regulatory edits, predominantly in non-coding regions. Genes with the largest numbers of regulatory edits are enriched in both apoptotic and innate immune pathways, providing a link between these known functions of ADAR and its role in cancer. We further show that despite a positive association between ADAR RNA expression and apoptotic and immune pathways, ADAR copy number is negatively associated with apoptosis and several immune cell types' signatures.