Precision cancer medicine promises to tailor clinical decisions to patients using genomic information. Indeed, successes of drugs targeting genetic alterations in tumors, such as imatinib that targets BCR-ABL in chronic myelogenous leukemia, have demonstrated the power of this approach. However, biological systems are complex, and patients may differ not only by the specific genetic alterations in their tumor, but also by more subtle interactions among such alterations. Systems biology and more specifically, network analysis, provides a framework for advancing precision medicine beyond clinical actionability of individual mutations. Here we discuss applications of network analysis to study tumor biology, early methods for N-of-1 tumor genome analysis, and the path for such tools to the clinic.

High-quality metadata annotations for data hosted in large public repositories are essential for research reproducibility and for conducting fast, powerful and scalable meta-analyses. Currently, a majority of sequencing samples in the National Center for Biotechnology Information's Sequence Read Archive (SRA) are missing metadata across several categories. In an effort to improve the metadata coverage of these samples, we leveraged almost 44 million attribute-value pairs from SRA BioSample to train a scalable, recurrent neural network that predicts missing metadata via named entity recognition (NER). The network was first trained to classify short text phrases according to 11 metadata categories and achieved an overall accuracy and area under the receiver operating characteristic curve of 85.2% and 0.977, respectively. We then applied our classifier to predict 11 metadata categories from the longer TITLE attribute of samples, evaluating performance on a set of samples withheld from model training. Prediction accuracies were high when extracting sample Genus/Species (94.85%), Condition/Disease (95.65%) and Strain (82.03%) from TITLEs, with lower accuracies and lack of predictions for other categories highlighting multiple issues with the current metadata annotations in BioSample. These results indicate the utility of recurrent neural networks for NER-based metadata prediction and the potential for models such as the one presented here to increase metadata coverage in BioSample while minimizing the need for manual curation. Database URL: https://github.com/cartercompbio/PredictMEE.

Anorexia nervosa (AN) is a psychiatric disorder affected by psychological, environmental, and biological factors. Individuals with AN avoid high-fat, high-calorie diets and have shown abnormal metabolism of fatty acids (FAs), which are essential for brain and cognitive/neuropsychiatric health. To clarify the relationship between FAs and AN, fasting and postprandial plasma FAs in AN patients and age-matched control women were analyzed via mass-spectrometry. Clinical phenotypes were assessed using Becker Anxiety Inventory and Becker Depression Inventory. AN patients and controls exhibited different FA signatures at both fasting and postprandial timepoints. Lauric acid, eicosapentaenoic acid (EPA), docosapentaenoic acid (DPA), and alpha-linoleic acid (ALA) were higher in AN than in controls (lauric acid: 15,081.6 ± 14,970.2 vs. 8257.4 ± 4740.2 pmol/mL; ALA at fasting: 2217.7 ± 1587.6 vs. 1087.9 ± 821.2 pmol/mL; ALA at postprandial: 1830.9 ± 1115.6 vs. 1159.4 ± 664.7 pmol/mL. EPA: 33,788.3 ± 17,487.5 vs. 22,860.6 ± 12,642.4 pmol/mL; DPA: 32,664.8 ± 16,215.0 vs. 20,969.0 ± 12,350.0 pmol/mL. FDR-adjusted p-values < 0.05). Food intake and AN status modified the correlations of FAs with body mass index (BMI), depression, and anxiety. Desaturases SCD-18 and D6D showed lower activities in AN compared to controls. Altered FA signature, specifically correlations between elevated n-3 FAs and worsened symptoms, illustrate metabolic underpinnings in AN. Future studies should investigate the mechanisms by which FA dysregulation, specifically elevated n-3 FAs, affects AN risk and outcome.

Tumor infiltrating lymphocytes (TILs) have been associated with favorable prognosis in multiple tumor types. The Cancer Genome Atlas (TCGA) represents the largest collection of cancer molecular data, but lacks detailed information about the immune environment. Here, we show that exome reads mapping to the complementarity-determining-region 3 (CDR3) of mature T-cell receptor beta (TCRB) can be used as an immune DNA (iDNA) signature. Specifically, we propose a method to identify CDR3 reads in a breast tumor exome and validate it using deep TCRB sequencing. In 1,078 TCGA breast cancer exomes, the fraction of CDR3 reads was associated with TILs fraction, tumor purity, adaptive immunity gene expression signatures and improved survival in Her2+ patients. Only 2/839 TCRB clonotypes were shared between patients and none associated with a specific HLA allele or somatic driver mutations. The iDNA biomarker enriches the comprehensive dataset collected through TCGA, revealing associations with other molecular features and clinical outcomes.

In vitro evolution and whole genome analysis has proven to be a powerful method for studying the mechanism of action of small molecules in many haploid microbes but has generally not been applied to human cell lines in part because their diploid state complicates the identification of variants that confer drug resistance. To determine if haploid human cells could be used in MOA studies, we evolved resistance to five different anticancer drugs (doxorubicin, gemcitabine, etoposide, topotecan, and paclitaxel) using a near-haploid cell line (HAP1) and then analyzed the genomes of the drug resistant clones, developing a bioinformatic pipeline that involved filtering for high frequency alleles predicted to change protein sequence, or alleles which appeared in the same gene for multiple independent selections with the same compound. Applying the filter to sequences from 28 drug resistant clones identified a set of 21 genes which was strongly enriched for known resistance genes or known drug targets (TOP1, TOP2A, DCK, WDR33, SLCO3A1). In addition, some lines carried structural variants that encompassed additional known resistance genes (ABCB1, WWOX and RRM1). Gene expression knockdown and knockout experiments of 10 validation targets showed a high degree of specificity and accuracy in our calls and demonstrates that the same drug resistance mechanisms found in diverse clinical samples can be evolved, discovered and studied in an isogenic background.

Cancer genomics has enabled the exhaustive molecular characterization of tumors and exposed hepatocellular carcinoma (HCC) as among the most complex cancers. This complexity is paralleled by dozens of mouse models that generate histologically similar tumors but have not been systematically validated at the molecular level. Accurate models of the molecular pathogenesis of HCC are essential for biomedical progress; therefore we compared genomic and transcriptomic profiles of four separate mouse models [MUP transgenic, TAK1-knockout, carcinogen-driven diethylnitrosamine (DEN), and Stelic Animal Model (STAM)] with those of 987 HCC patients with distinct etiologies. These four models differed substantially in their mutational load, mutational signatures, affected genes and pathways, and transcriptomes. STAM tumors were most molecularly similar to human HCC, with frequent mutations in Ctnnb1, similar pathway alterations, and high transcriptomic similarity to high-grade, proliferative human tumors with poor prognosis. In contrast, TAK1 tumors better reflected the mutational signature of human HCC and were transcriptionally similar to low-grade human tumors. DEN tumors were least similar to human disease and almost universally carried the Braf V637E mutation, which is rarely found in human HCC. Immune analysis revealed that strain-specific MHC-I genotype can influence the molecular makeup of murine tumors. Thus, different mouse models of HCC recapitulate distinct aspects of HCC biology, and their use should be adapted to specific questions based on the molecular features provided here.

Large-scale chromosomal translocations are frequent oncogenic drivers in acute myeloid leukemia (AML). These translocations often occur in critical transcriptional/epigenetic regulators and contribute to malignant cell growth through alteration of normal gene expression. Despite this knowledge, the specific gene expression alterations that contribute to the development of leukemia remain incompletely understood. Here, through characterization of transcriptional regulation by the RUNX1-ETO fusion protein, we have identified Ras-association domain family member 2 (RASSF2) as a critical gene that is aberrantly transcriptionally repressed in t(8;21)-associated AML. Re-expression of RASSF2 specifically inhibits t(8;21) AML development in multiple models. Through biochemical and functional studies, we demonstrate RASSF2-mediated functions to be dependent on interaction with Hippo kinases, MST1 and MST2, but independent of canonical Hippo pathway signaling. Using proximity-based biotin labeling we define the RASSF2-proximal proteome in leukemia cells and reveal association with Rac GTPase-related proteins, including an interaction with the guanine nucleotide exchange factor, DOCK2. Importantly, RASSF2 knockdown impairs Rac GTPase activation, and RASSF2 expression is broadly correlated with Rac-mediated signal transduction in AML patients. Together, these data reveal a previously unappreciated mechanistic link between RASSF2, Hippo kinases, and Rac activity with potentially broad functional consequences in leukemia.

Many cancers evade immune rejection by suppressing major histocompatibility class I (MHC-I) antigen processing and presentation (AgPP). Such cancers do not respond to immune checkpoint inhibitor therapies (ICIT) such as PD-1/PD-L1 [PD-(L)1] blockade. Certain chemotherapeutic drugs augment tumor control by PD-(L)1 inhibitors through potentiation of T-cell priming but whether and how chemotherapy enhances MHC-I-dependent cancer cell recognition by cytotoxic T cells (CTLs) is not entirely clear. We now show that the lysine acetyl transferases p300/CREB binding protein (CBP) control MHC-I AgPPM expression and neoantigen amounts in human cancers. Moreover, we found that two distinct DNA damaging drugs, the platinoid oxaliplatin and the topoisomerase inhibitor mitoxantrone, strongly up-regulate MHC-I AgPP in a manner dependent on activation of nuclear factor kappa B (NF-κB), p300/CBP, and other transcription factors, but independently of autocrine IFNγ signaling. Accordingly, NF-κB and p300 ablations prevent chemotherapy-induced MHC-I AgPP and abrogate rejection of low MHC-I-expressing tumors by reinvigorated CD8⁺ CTLs. Drugs like oxaliplatin and mitoxantrone may be used to overcome resistance to PD-(L)1 inhibitors in tumors that had "epigenetically down-regulated," but had not permanently lost MHC-I AgPP activity.