Functional genomic screens with death rate analyses reveal mechanisms of drug action.

A common approach for understanding how drugs induce their therapeutic effects is to identify the genetic determinants of drug sensitivity. Because 'chemo-genetic profiles' are performed in a pooled format, inference of gene function is subject to several confounding influences related to variation in growth rates between clones. In this study, we developed Method for Evaluating Death Using a Simulation-assisted Approach (MEDUSA), which uses time-resolved measurements, along with model-driven constraints, to reveal the combination of growth and death rates that generated the observed drug response. MEDUSA is uniquely effective at identifying death regulatory genes. We apply MEDUSA to characterize DNA damage-induced lethality in the presence and absence of p53. Loss of p53 switches the mechanism of DNA damage-induced death from apoptosis to a non-apoptotic death that requires high respiration. These findings demonstrate the utility of MEDUSA both for determining the genetic dependencies of lethality and for revealing opportunities to potentiate chemo-efficacy in a cancer-specific manner.

p53 controls choice between apoptotic and non-apoptotic death following DNA damage.

DNA damage can activate apoptotic and non-apoptotic forms of cell death; however, it remains unclear what features dictate which type of cell death is activated. We report that p53 controls the choice between apoptotic and non-apoptotic death following exposure to DNA damage. In contrast to the conventional model, which suggests that p53-deficient cells should be resistant to DNA damage-induced cell death, we find that p53-deficient cells die at high rates following DNA damage, but exclusively using non-apoptotic mechanisms. Our experimental data and computational modeling reveal that non-apoptotic death in p53-deficient cells has not been observed due to use of assays that are either insensitive to cell death, or that specifically score apoptotic cells. Using functional genetic screening - with an analysis that enables computational inference of the drug-induced death rate - we find in p53-deficient cells that DNA damage activates a mitochondrial respiration-dependent form of cell death, called MPT-driven necrosis. Cells deficient for p53 have high basal respiration, which primes MPT-driven necrosis. Finally, using metabolite profiling, we identified mitochondrial activity-dependent metabolic vulnerabilities that can be targeted to potentiate the lethality of DNA damage specifically in p53-deficient cells. Our findings reveal how the dual functions of p53 in regulating mitochondrial activity and the DNA damage response combine to facilitate the choice between apoptotic and non-apoptotic death.

Spontaneous Propionibacterium Acnes abscess with intraventricular rupture in an immunocompetent adult without prior neurosurgical intervention.

Previously viewed as a culture contaminant, Propionibacterium Acnes can cause infection following neurosurgical intervention. Its role in brain abscess in the immunocompetent, surgically naïve population has been infrequently reported. Herein, we describe an immunocompetent 55-year-old man with no risk factors found to have a thalamic abscess with intraventricular rupture.

Utilization of anterior lumbar interbody fusion for severe kyphotic deformity secondary to Pott's disease: illustrative case.

BACKGROUND: Spinal tuberculosis may result in severe kyphotic deformity. Effective restoration of lordosis and correction of sagittal balance often requires invasive osteotomies associated with significant morbidity. The advantages of focusing on symptomatic management and staging in the initial treatment of these deformities have not been well reported to date. OBSERVATIONS: The authors reported the case of a 64-year-old Vietnamese woman with a history of spinal tuberculosis who underwent anterior lumbar interbody fusion (ALIF) for symptomatic treatment of L5-S1 radiculopathy resulting from fixed kyphotic deformity. Postoperatively, the patient experienced near immediate symptom improvement, and radiographic evidence at 1-year follow-up showed continued lordotic correction of 30° as well as stable sagittal balance. LESSONS: In this case, an L5-S1 ALIF was sufficient to treat the patient's acute symptoms and provided satisfactory correction of a tuberculosis-associated fixed kyphotic deformity while effectively delaying more invasive measures, such as a vertebral column resection. Patients with adult spinal deformity may benefit from less invasive staging procedures before treating these deformities with larger surgeries.

ELP-dependent expression of MCL1 promotes resistance to EGFR inhibition in triple-negative breast cancer cells.

Targeted therapeutics for cancer generally exploit "oncogene addiction," a phenomenon in which the growth and survival of tumor cells depend on the activity of a particular protein. However, the efficacy of oncogene-targeted therapies varies substantially. For instance, targeting epidermal growth factor receptor (EGFR) signaling is effective in some non-small cell lung cancer (NSCLC) but not in triple-negative breast cancer (TNBC), although these cancers show a similar degree of increase in EGFR activity. Using a genome-wide CRISPR-Cas9 genetic knockout screen, we found that the Elongator (ELP) complex mediates insensitivity to the EGFR inhibitor erlotinib in TNBC cells by promoting the synthesis of the antiapoptotic protein Mcl-1. Depleting ELP proteins promoted apoptotic cell death in an EGFR inhibition-dependent manner. Pharmacological inhibition of Mcl-1 synergized with EGFR inhibition in a panel of genetically diverse TNBC cells. The findings indicate that TNBC "addiction" to EGFR signaling is masked by the ELP complex and that resistance to EGFR inhibitors in TNBC might be overcome by cotargeting Mcl-1.

Drug antagonism and single-agent dominance result from differences in death kinetics.

Cancer treatment generally involves drugs used in combinations. Most previous work has focused on identifying and understanding synergistic drug-drug interactions; however, understanding antagonistic interactions remains an important and understudied issue. To enrich for antagonism and reveal common features of these combinations, we screened all pairwise combinations of drugs characterized as activators of regulated cell death. This network is strongly enriched for antagonism, particularly a form of antagonism that we call 'single-agent dominance'. Single-agent dominance refers to antagonisms in which a two-drug combination phenocopies one of the two agents. Dominance results from differences in cell death onset time, with dominant drugs acting earlier than their suppressed counterparts. We explored mechanisms by which parthanatotic agents dominate apoptotic agents, finding that dominance in this scenario is caused by mutually exclusive and conflicting use of Poly(ADP-ribose) polymerase 1 (PARP1). Taken together, our study reveals death kinetics as a predictive feature of antagonism, due to inhibitory crosstalk between cell death pathways.

Beclin 1 Promotes Endosome Recruitment of Hepatocyte Growth Factor Tyrosine Kinase Substrate to Suppress Tumor Proliferation.

Beclin 1 has nonautophagic functions that include its ability to regulate endocytic receptor trafficking. However, the contribution of this function to tumor suppression is poorly understood. Here, we provide in vivo evidence that Beclin 1 suppresses tumor proliferation by regulating the endocytic trafficking and degradation of the EGFR and transferrin (TFR1) receptors. Beclin 1 promoted endosomal recruitment of hepatocyte growth factor tyrosine kinase substrate (HRS), which was necessary for sorting surface receptors to intraluminal vesicles for signal silencing and lysosomal degradation. In tumors with low Beclin 1 expression, endosomal HRS recruitment was diminished and receptor function was sustained. Collectively, our results demonstrate a novel role for Beclin 1 in impeding tumor growth by coordinating the regulation of key growth factor and nutrient receptors. These data provide an explanation for how low levels of Beclin 1 facilitate tumor proliferation and contribute to poor cancer outcomes. SIGNIFICANCE: Beclin 1 controls the trafficking fate of growth regulatory receptors to suppress tumor proliferation.

Tumor-stroma interactions differentially alter drug sensitivity based on the origin of stromal cells.

Due to tumor heterogeneity, most believe that effective treatments should be tailored to the features of an individual tumor or tumor subclass. It is still unclear, however, what information should be considered for optimal disease stratification, and most prior work focuses on tumor genomics. Here, we focus on the tumor microenvironment. Using a large-scale coculture assay optimized to measure drug-induced cell death, we identify tumor-stroma interactions that modulate drug sensitivity. Our data show that the chemo-insensitivity typically associated with aggressive subtypes of breast cancer is not observed if these cells are grown in 2D or 3D monoculture, but is manifested when these cells are cocultured with stromal cells, such as fibroblasts. Furthermore, we find that fibroblasts influence drug responses in two distinct and divergent manners, associated with the tissue from which the fibroblasts were harvested. These divergent phenotypes occur regardless of the drug tested and result from modulation of apoptotic priming within tumor cells. Our study highlights unexpected diversity in tumor-stroma interactions, and we reveal new principles that dictate how fibroblasts alter tumor drug responses.

Whole-exome sequencing of circulating tumor cells provides a window into metastatic prostate cancer.

Comprehensive analyses of cancer genomes promise to inform prognoses and precise cancer treatments. A major barrier, however, is inaccessibility of metastatic tissue. A potential solution is to characterize circulating tumor cells (CTCs), but this requires overcoming the challenges of isolating rare cells and sequencing low-input material. Here we report an integrated process to isolate, qualify and sequence whole exomes of CTCs with high fidelity using a census-based sequencing strategy. Power calculations suggest that mapping of >99.995% of the standard exome is possible in CTCs. We validated our process in two patients with prostate cancer, including one for whom we sequenced CTCs, a lymph node metastasis and nine cores of the primary tumor. Fifty-one of 73 CTC mutations (70%) were present in matched tissue. Moreover, we identified 10 early trunk and 56 metastatic trunk mutations in the non-CTC tumor samples and found 90% and 73% of these mutations, respectively, in CTC exomes. This study establishes a foundation for CTC genomics in the clinic.

Widespread genetic heterogeneity in multiple myeloma: implications for targeted therapy.

We performed massively parallel sequencing of paired tumor/normal samples from 203 multiple myeloma (MM) patients and identified significantly mutated genes and copy number alterations and discovered putative tumor suppressor genes by determining homozygous deletions and loss of heterozygosity. We observed frequent mutations in KRAS (particularly in previously treated patients), NRAS, BRAF, FAM46C, TP53, and DIS3 (particularly in nonhyperdiploid MM). Mutations were often present in subclonal populations, and multiple mutations within the same pathway (e.g., KRAS, NRAS, and BRAF) were observed in the same patient. In vitro modeling predicts only partial treatment efficacy of targeting subclonal mutations, and even growth promotion of nonmutated subclones in some cases. These results emphasize the importance of heterogeneity analysis for treatment decisions.

Discovery and prioritization of somatic mutations in diffuse large B-cell lymphoma (DLBCL) by whole-exome sequencing.

To gain insight into the genomic basis of diffuse large B-cell lymphoma (DLBCL), we performed massively parallel whole-exome sequencing of 55 primary tumor samples from patients with DLBCL and matched normal tissue. We identified recurrent mutations in genes that are well known to be functionally relevant in DLBCL, including MYD88, CARD11, EZH2, and CREBBP. We also identified somatic mutations in genes for which a functional role in DLBCL has not been previously suspected. These genes include MEF2B, MLL2, BTG1, GNA13, ACTB, P2RY8, PCLO, and TNFRSF14. Further, we show that BCL2 mutations commonly occur in patients with BCL2/IgH rearrangements as a result of somatic hypermutation normally occurring at the IgH locus. The BCL2 point mutations are primarily synonymous, and likely caused by activation-induced cytidine deaminase-mediated somatic hypermutation, as shown by comprehensive analysis of enrichment of mutations in WRCY target motifs. Those nonsynonymous mutations that are observed tend to be found outside of the functionally important BH domains of the protein, suggesting that strong negative selection against BCL2 loss-of-function mutations is at play. Last, by using an algorithm designed to identify likely functionally relevant but infrequent mutations, we identify KRAS, BRAF, and NOTCH1 as likely drivers of DLBCL pathogenesis in some patients. Our data provide an unbiased view of the landscape of mutations in DLBCL, and this in turn may point toward new therapeutic strategies for the disease.

Extensive changes in the expression of the opioid genes between humans and chimpanzees.

The various means by which the body perceives, transmits, and resolves the experiences of pain and nociception are mediated by a host of molecules, including neuropeptides within the opioid gene signaling pathway. The peptide ligands and receptors encoded by this group of genes have been linked to behavioral disorders as well as a number of psychiatric affective disorders. Our aim was to explore the recent evolutionary history of these two gene families by taking a comparative genomics approach, specifically through a comparison between humans and chimpanzees. Our analyses indicate differential expression of these genes between the two species, more than expected based on genome-wide comparisons, indicating that differential expression is pervasive among the opioid genes. Of the 8 family members, three genes showed significant expression differences (PENK, PNOC, and OPRL1), with two others marginally significant (OPRM1 and OPRD1). Accelerated substitution rates along human and chimpanzee lineages within the putative regulatory regions of OPRM1, POMC, and PDYN between the human and chimpanzee branches are consistent with positive selection. Collectively, these results suggest that there may have been a selective advantage to modulating the expression of the opioid genes in humans compared with our closest living relatives. Information about the cognitive roles mediated by these genes in humans may help to elucidate the trait consequences of these putatively adaptive expression changes.

A pipeline to determine RT-QPCR control genes for evolutionary studies: application to primate gene expression across multiple tissues.

Because many species-specific phenotypic differences are assumed to be caused by differential regulation of gene expression, many recent investigations have focused on measuring transcript abundance. Despite the availability of high-throughput platforms, quantitative real-time polymerase chain reaction (RT-QPCR) is often the method of choice because of its low cost and wider dynamic range. However, the accuracy of this technique heavily relies on the use of multiple valid control genes for normalization. We created a pipeline for choosing genes potentially useful as RT-QPCR control genes for measuring expression between human and chimpanzee samples across multiple tissues, using published microarrays and a measure of tissue-specificity. We identified 13 genes from the pipeline and from commonly used control genes: ACTB, USP49, ARGHGEF2, GSK3A, TBP, SDHA, EIF2B2, GPDH, YWHAZ, HPTR1, RPL13A, HMBS, and EEF2. We then tested these candidate genes and validated their expression stability across species. We established the rank order of the most preferable set of genes for single and combined tissues. Our results suggest that for at least three tissues (cerebral cortex, liver, and skeletal muscle), EIF2B2, EEF2, HMBS, and SDHA are useful genes for normalizing human and chimpanzee expression using RT-QPCR. Interestingly, other commonly used control genes, including TBP, GAPDH, and, especially ACTB do not perform as well. This pipeline could be easily adapted to other species for which expression data exist, providing taxonomically appropriate control genes for comparisons of gene expression among species.