Bhakti Dwivedi

Cellular and Molecular Identity of Tumor-Associated Macrophages in Glioblastoma

In glioblastoma (GBM), tumor-associated macrophages (TAM) represent up to one half of the cells of the tumor mass, including both infiltrating macrophages and resident brain microglia. In an effort to delineate the temporal and spatial dynamics of TAM composition during gliomagenesis, we used genetically engineered and GL261-induced mouse models in combination with CX3CR1GFP/WT;CCR2RFP/WT double knock-in mice. Using this approach, we demonstrated that CX3CR1LoCCR2Hi monocytes were recruited to the GBM, where they transitioned to CX3CR1HiCCR2Lo macrophages and CX3CR1HiCCR2- microglia-like cells. Infiltrating macrophages/monocytes constituted approximately 85% of the total TAM population, with resident microglia accounting for the approximately 15% remaining. Bone marrow-derived infiltrating macrophages/monocytes were recruited to the tumor early during GBM initiation, where they localized preferentially to perivascular areas. In contrast, resident microglia were localized mainly to peritumoral regions. RNA-sequencing analyses revealed differential gene expression patterns unique to infiltrating and resident cells, suggesting unique functions for each TAM population. Notably, limiting monocyte infiltration via genetic Ccl2 reduction prolonged the survival of tumor-bearing mice. Our findings illuminate the unique composition and functions of infiltrating and resident myeloid cells in GBM, establishing a rationale to target infiltrating cells in this neoplasm. Cancer Res; 77(9); 2266-78. ©2017 AACR.

Fatigue is associated with inflammation in patients with head and neck cancer before and after intensity-modulated radiation therapy

Patients with head and neck cancer (HNC) receiving intensity-modulated radiation therapy (IMRT) have particularly high rates of fatigue, and pre- and post-radiotherapy fatigue are prognostic factors for pathologic tumor responses and poor survival. Although inflammation has been proposed as one of the potential mechanisms of fatigue in cancer patients, findings have not been consistent, and there is a dearth of longitudinal studies. Accordingly, we conducted a prospective study in 46 HNC patients pre- and one-month post-IMRT. Fatigue was measured by the Multidimensional Fatigue Inventory (MFI)-20 at both time points along with the assessment of peripheral blood inflammatory markers including interleukin (IL)-6, soluble tumor necrosis factor receptor 2, and C-reactive protein (CRP) and gene expression. Generalized estimating equations were used to examine the association between inflammatory markers and fatigue. Gene enrichment analysis using MetaCore software was performed using up-regulated genes that were significantly associated with IMRT and fatigue. Significant associations between fatigue and IL-6 as well as CRP, which were independent of time, were observed. In addition the change in fatigue from pre- to post-IMRT was positively associated with the change in IL-6 and CRP. Analysis of up-regulated gene transcripts as a function of IMRT and fatigue revealed overrepresentation of transcripts related to the defense response and nuclear factor kappa B. In conclusion, our findings support the hypotheses that inflammation is associated with fatigue over time in HNC patients. Future studies on how inflammation contributes to fatigue as well as strategies targeting inflammation to reduce fatigue are warranted.

Image-guided genomics of phenotypically heterogeneous populations reveals vascular signalling during symbiotic collective cancer invasion

Phenotypic heterogeneity is widely observed in cancer cell populations. Here, to probe this heterogeneity, we developed an image-guided genomics technique termed spatiotemporal genomic and cellular analysis (SaGA) that allows for precise selection and amplification of living and rare cells. SaGA was used on collectively invading 3D cancer cell packs to create purified leader and follower cell lines. The leader cell cultures are phenotypically stable and highly invasive in contrast to follower cultures, which show phenotypic plasticity over time and minimally invade in a sheet-like pattern. Genomic and molecular interrogation reveals an atypical VEGF-based vasculogenesis signalling that facilitates recruitment of follower cells but not for leader cell motility itself, which instead utilizes focal adhesion kinase-fibronectin signalling. While leader cells provide an escape mechanism for followers, follower cells in turn provide leaders with increased growth and survival. These data support a symbiotic model of collective invasion where phenotypically distinct cell types cooperate to promote their escape.

Gene integrated set profile analysis: a context-based approach for inferring biological endpoints

The identification of genes with specific patterns of change (e.g. down-regulated and methylated) as phenotype drivers or samples with similar profiles for a given gene set as drivers of clinical outcome, requires the integration of several genomic data types for which an ‘integrate by intersection’ (IBI) approach is often applied. In this approach, results from separate analyses of each data type are intersected, which has the limitation of a smaller intersection with more data types. We introduce a new method, GISPA (Gene Integrated Set Profile Analysis) for integrated genomic analysis and its variation, SISPA (Sample Integrated Set Profile Analysis) for defining respective genes and samples with the context of similar, a priori specified molecular profiles. With GISPA, the user defines a molecular profile that is compared among several classes and obtains ranked gene sets that satisfy the profile as drivers of each class. With SISPA, the user defines a gene set that satisfies a profile and obtains sample groups of profile activity. Our results from applying GISPA to human multiple myeloma (MM) cell lines contained genes of known profiles and importance, along with several novel targets, and their further SISPA application to MM coMMpass trial data showed clinical relevance.

Factors affecting the persistence of drug-induced reprogramming of the cancer methylome.

ABSTRACT Aberrant DNA methylation is a critical feature of cancer. Epigenetic therapy seeks to reverse these changes to restore normal gene expression. DNA demethylating agents, including 5-aza-2′-deoxycytidine (DAC), are currently used to treat certain leukemias, and can sensitize solid tumors to chemotherapy and immunotherapy. However, it has been difficult to pin the clinical efficacy of these agents to specific demethylation events, and the factors that contribute to the durability of response remain largely unknown. Here we examined the genome-wide kinetics of DAC-induced DNA demethylation and subsequent remethylation after drug withdrawal in breast cancer cells. We find that CpGs differ in both their susceptibility to demethylation and propensity for remethylation after drug removal. DAC-induced demethylation was most apparent at CpGs with higher initial methylation levels and further from CpG islands. Once demethylated, such sites exhibited varied remethylation potentials. The most rapidly remethylating CpGs regained >75% of their starting methylation within a month of drug withdrawal. These sites had higher pretreatment methylation levels, were enriched in gene bodies, marked by H3K36me3, and tended to be methylated in normal breast cells. In contrast, a more resistant class of CpG sites failed to regain even 20% of their initial methylation after 3 months. These sites had lower pretreatment methylation levels, were within or near CpG islands, marked by H3K79me2 or H3K4me2/3, and were overrepresented in sites that become aberrantly hypermethylated in breast cancers. Thus, whereas DAC-induced demethylation affects both endogenous and aberrantly methylated sites, tumor-specific hypermethylation is more slowly regained, even as normal methylation promptly recovers. Taken together, these data suggest that the durability of DAC response is linked to its selective ability to stably reset at least a portion of the cancer methylome.

MERTK inhibition induces polyploidy and promotes cell death and cellular senescence in glioblastoma multiforme

Background MER receptor tyrosine kinase (MERTK) is expressed in a variety of malignancies, including glioblastoma multiforme (GBM). Our previous work demonstrated that inhibition of MERTK using RNA interference induced cell death and chemosensitivity in GBM cells, implicating MERTK as a potential therapeutic target. Here we investigate whether a novel MERTK-selective small molecule tyrosine kinase inhibitor, UNC2025, has similar anti-tumor effects in GBM cell lines. Methods Correlations between expression of GAS6, a MERTK ligand, and prognosis were determined using data from the TCGA database. GBM cell lines (A172, SF188, U251) were treated in vitro with increasing doses of UNC2025 (50-400nM). Cell count and viability were determined by trypan blue exclusion. Cell cycle profiles and induction of apoptosis were assessed by flow cytometric analysis after BrdU or Po-Pro-1/propidium iodide staining, respectively. Polyploidy was detected by propidium iodide staining and metaphase spread. Cellular senescence was determined by β-galactosidase staining and senescence-associated secretory cytokine analysis. Results Decreased overall survival significantly correlated with high levels of GAS6 expression in GBM, highlighting the importance of TAM kinase signaling in GBM tumorigenesis and/or therapy resistance and providing strong rationale for targeting these pathways in the clinic. All three GBM cell lines exhibited dose dependent reductions in cell number and colony formation (>90% at 200nM) after treatment with UNC2025. Cell cycle analysis demonstrated accumulation of cells in the G2/M phase and development of polyploidy. After extended exposure, 60–80% of cells underwent apoptosis. The majority of surviving cells (65–95%) were senescent and did not recover after drug removal. Thus, UNC2025 mediates anti-tumor activity in GBM by multiple mechanisms. Conclusions The findings described here provide further evidence of oncogenic roles for MERTK in GBM, demonstrate the importance of kinase activity for MERTK tumorigenicity and validate UNC2025, a novel MERTK inhibitor, as a potential therapeutic agent for treatment of GBM.

Low expression of pro-apoptotic Bcl-2 family proteins sets the apoptotic threshold in Waldenstrom macroglobulinemia

Waldenström macroglobulinemia (WM) is a proliferative disorder of IgM-secreting, lymphoplasmacytoid cells that inhabit the lymph nodes and bone marrow. The disease carries a high prevalence of activating mutations in MyD88 (91%) and CXCR4 (28%). Because signaling through these pathways leads to Bcl-xL induction, we examined Bcl-2 family expression in WM patients and cell lines. Unlike other B-lymphocyte-derived malignancies, which become dependent on expression of anti-apoptotic proteins to counter expression of pro-apoptotic proteins, WM samples expressed both pro- and anti-apoptotic Bcl-2 proteins at low levels similar to their normal B-cell and plasma cell counterparts. Three WM cell lines expressed pro-apoptotic Bcl-2 family members Bim or Bax and Bak at low levels, which determined their sensitivity to inducers of intrinsic apoptosis. In two cell lines, miR-155 upregulation, which is common in WM, was responsible for the inhibition of FOXO3a and Bim expression. Both antagonizing miR-155 to induce Bim and proteasome inhibition increased the sensitivity to ABT-737 in these lines indicating a lowering of the apoptotic threshold. In this manner, treatments that increase pro-apoptotic protein expression increase the efficacy of agents treated in combination in addition to direct killing.

Galactic Cosmic Radiation Induces Persistent Epigenome Alterations Relevant to Human Lung Cancer

Human deep space and planetary travel is limited by uncertainties regarding the health risks associated with exposure to galactic cosmic radiation (GCR), and in particular the high linear energy transfer (LET), heavy ion component. Here we assessed the impact of two high-LET ions 56Fe and 28Si, and low-LET X rays on genome-wide methylation patterns in human bronchial epithelial cells. We found that all three radiation types induced rapid and stable changes in DNA methylation but at distinct subsets of CpG sites affecting different chromatin compartments. The 56Fe ions induced mostly hypermethylation, and primarily affected sites in open chromatin regions including enhancers, promoters and the edges (“shores”) of CpG islands. The 28Si ion-exposure had mixed effects, inducing both hyper and hypomethylation and affecting sites in more repressed heterochromatic environments, whereas X rays induced mostly hypomethylation, primarily at sites in gene bodies and intergenic regions. Significantly, the methylation status of 56Fe ion sensitive sites, but not those affected by X ray or 28Si ions, discriminated tumor from normal tissue for human lung adenocarcinomas and squamous cell carcinomas. Thus, high-LET radiation exposure leaves a lasting imprint on the epigenome, and affects sites relevant to human lung cancer. These methylation signatures may prove useful in monitoring the cumulative biological impact and associated cancer risks encountered by astronauts in deep space.

TNFRSF10C copy number variation is associated with metastatic colorectal cancer.

BACKGROUND Genetic markers for distant metastatic disease in patients with colorectal cancer (CRC) are not well defined. Identification of genetic alterations associated with metastatic CRC could help to guide systemic and local treatment strategies. We evaluated the association of tumor necrosis factor receptor superfamily member 10C (TNFRSF10C) copy number variation (CNV) with distant metastatic disease in patients with CRC using The Cancer Genome Atlas (TCGA). METHODS Genetic sequencing data and clinical characteristics were obtained from TCGA for all available patients with CRC. There were 515 CRC patient samples with CNV and clinical outcome data, including a subset of 144 rectal adenocarcinoma patient samples. Using the TCGA CRC dataset, CNV of TNFRSF10C was evaluated for association with distant metastatic disease (M1 vs. M0). Multivariate logistic regression analysis with odds ratio (OR) using a 95% confidence interval (CI) was performed adjusting for age, T stage, N stage, adjuvant chemotherapy, gender, microsatellite instability (MSI), location, and surgical margin status. RESULTS TNFRSF10C CNV in patients with CRC was associated with distant metastatic disease [OR 4.81 (95% CI, 2.13-10.85) P<0.001] and positive lymph nodes [OR 18.83 (95% CI, 8.42-42.09)]; P<0.001) but not MSI (OR P=0.799). On multivariate analysis, after adjusting for pathologic T stage, N stage, adjuvant chemotherapy, gender, and MSI, TNFRSF10C CNV remained significantly associated with distant metastatic disease (OR P=0.018). Subset analysis revealed that TNFRSF10C CNV was also significantly associated with distant metastatic disease in patients with rectal adenocarcinoma (OR P=0.016). CONCLUSIONS TNFRSF10C CNV in patients with CRC is associated with distant metastatic disease. With further validation, such genetic profiles could be used clinically to support optimal systemic treatment strategies versus more aggressive local therapies in patients with CRC, including radiation therapy for rectal adenocarcinoma.

NOJAH: Not Just Another Heatmap for Genome-Wide Cluster Analysis

Since their inception, several tools have been developed for cluster analysis and heatmap construction. The application of such tools to the number and types of genome-wide data available from next generation sequencing (NGS) technologies requires the adaptation of statistical concepts, such as in defining a most variable gene set, and more intricate cluster analyses method to address multiple omic data types. Additionally, the growing number of publicly available datasets has created the desire to estimate the statistical significance of a gene signature derived from one dataset to similarly group samples based on another dataset. The currently available number of tools and their combined use for generating heatmaps, along with the several adaptations of statistical concepts for addressing the higher dimensionality of genome-wide NGS-derived data, has created a further challenge in the ability to replicate heatmap results. We introduce NOJAH (NOt Just Another Heatmap), an interactive tool that defines and implements a workflow for genome-wide cluster analysis and heatmap construction by creating and combining several tools into a single user interface. NOJAH includes several newly developed scripts for techniques that though frequently applied are not sufficiently documented to allow for replicability of results. These techniques include: defining a most variable gene set (a.k.a., ‘core genes’), estimating the statistical significance of a gene signature to separate samples into clusters, and performing a result merging integrated cluster analysis. With only a user uploaded dataset, NOJAH provides as output, among other things, the minimum documentation required for replicating heatmap results. Additionally, NOJAH contains five different existing R packages that are connected in the interface by their functionality as part of a defined workflow for genome-wide cluster analysis. The NOJAH application tool is available at http://bbisr.shinyapps.winship.emory.edu/NOJAH/ with corresponding source code available at https://github.com/bbisr-shinyapps/NOJAH/.