Gyan Bhanot

Autophagy Suppresses Tumorigenesis through Elimination of p62

Allelic loss of the essential autophagy gene beclin1 occurs in human cancers and renders mice tumor-prone suggesting that autophagy is a tumor-suppression mechanism. While tumor cells utilize autophagy to survive metabolic stress, autophagy also mitigates the resulting cellular damage that may limit tumorigenesis. In response to stress, autophagy-defective tumor cells preferentially accumulated p62/SQSTM1 (p62), endoplasmic reticulum (ER) chaperones, damaged mitochondria, reactive oxygen species (ROS), and genome damage. Moreover, suppressing ROS or p62 accumulation prevented damage resulting from autophagy defects indicating that failure to regulate p62 caused oxidative stress. Importantly, sustained p62 expression resulting from autophagy defects was sufficient to alter NF-kappaB regulation and gene expression and to promote tumorigenesis. Thus, defective autophagy is a mechanism for p62 upregulation commonly observed in human tumors that contributes directly to tumorigenesis likely by perturbing the signal transduction adaptor function of p62-controlling pathways critical for oncogenesis.

Single Cell Profiling of Circulating Tumor Cells: Transcriptional Heterogeneity and Diversity from Breast Cancer Cell Lines

Background To improve cancer therapy, it is critical to target metastasizing cells. Circulating tumor cells (CTCs) are rare cells found in the blood of patients with solid tumors and may play a key role in cancer dissemination. Uncovering CTC phenotypes offers a potential avenue to inform treatment. However, CTC transcriptional profiling is limited by leukocyte contamination; an approach to surmount this problem is single cell analysis. Here we demonstrate feasibility of performing high dimensional single CTC profiling, providing early insight into CTC heterogeneity and allowing comparisons to breast cancer cell lines widely used for drug discovery. Methodology/Principal Findings We purified CTCs using the MagSweeper, an immunomagnetic enrichment device that isolates live tumor cells from unfractionated blood. CTCs that met stringent criteria for further analysis were obtained from 70% (14/20) of primary and 70% (21/30) of metastatic breast cancer patients; none were captured from patients with non-epithelial cancer (n = 20) or healthy subjects (n = 25). Microfluidic-based single cell transcriptional profiling of 87 cancer-associated and reference genes showed heterogeneity among individual CTCs, separating them into two major subgroups, based on 31 highly expressed genes. In contrast, single cells from seven breast cancer cell lines were tightly clustered together by sample ID and ER status. CTC profiles were distinct from those of cancer cell lines, questioning the suitability of such lines for drug discovery efforts for late stage cancer therapy. Conclusions/Significance For the first time, we directly measured high dimensional gene expression in individual CTCs without the common practice of pooling such cells. Elevated transcript levels of genes associated with metastasis NPTN, S100A4, S100A9, and with epithelial mesenchymal transition: VIM, TGFß1, ZEB2, FOXC1, CXCR4, were striking compared to cell lines. Our findings demonstrate that profiling CTCs on a cell-by-cell basis is possible and may facilitate the application of ‘liquid biopsies’ to better model drug discovery.

The Quenched Eguchi-Kawai Model

Abstract It is argued that spontaneous symmetry breaking occurs in the recently proposed Eguchi-Kawai model. The argument is based on an analytic investigation and on Monte Carlo simulations. A quenched version of the model is proposed which gives good behavior at weak couplings.

Identification of a MicroRNA Panel for Clear-cell Kidney Cancer

OBJECTIVES To identify a robust panel of microRNA signatures that can classify tumor from normal kidney using microRNA expression levels. Mounting evidence suggests that microRNAs are key players in essential cellular processes and that their expression pattern can serve as diagnostic biomarkers for cancerous tissues. METHODS We selected 28 clear-cell type human renal cell carcinoma (ccRCC), samples from patient-matched specimens to perform high-throughput, quantitative real-time polymerase chain reaction analysis of microRNA expression levels. The data were subjected to rigorous statistical analyses and hierarchical clustering to produce a discrete set of microRNAs that can robustly distinguish ccRCC from their patient-matched normal kidney tissue samples with high confidence. RESULTS Thirty-five microRNAs were found that can robustly distinguish ccRCC from their patient-matched normal kidney tissue samples with high confidence. Among this set of 35 signature microRNAs, 26 were found to be consistently downregulated and 9 consistently upregulated in ccRCC relative to normal kidney samples. Two microRNAs, namely, MiR-155 and miR-21, commonly found to be upregulated in other cancers, and miR-210, induced by hypoxia, were also identified as overexpressed in ccRCC in our study. MicroRNAs identified as downregulated in our study can be correlated to common chromosome deletions in ccRCC. CONCLUSIONS Our analysis is a comprehensive, statistically relevant study that identifies the microRNAs dysregulated in ccRCC, which can serve as the basis of molecular markers for diagnosis.

Molecular Stratification of Clear Cell Renal Cell Carcinoma by Consensus Clustering Reveals Distinct Subtypes and Survival Patterns

Clear cell renal cell carcinoma (ccRCC) is the predominant RCC subtype, but even within this classification, the natural history is heterogeneous and difficult to predict. A sophisticated understanding of the molecular features most discriminatory for the underlying tumor heterogeneity should be predicated on identifiable and biologically meaningful patterns of gene expression. Gene expression microarray data were analyzed using software that implements iterative unsupervised consensus clustering algorithms to identify the optimal molecular subclasses, without clinical or other classifying information. ConsensusCluster analysis identified two distinct subtypes of ccRCC within the training set, designated clear cell type A (ccA) and B (ccB). Based on the core tumors, or most well-defined arrays, in each subtype, logical analysis of data (LAD) defined a small, highly predictive gene set that could then be used to classify additional tumors individually. The subclasses were corroborated in a validation data set of 177 tumors and analyzed for clinical outcome. Based on individual tumor assignment, tumors designated ccA have markedly improved disease-specific survival compared to ccB (median survival of 8.6 vs 2.0 years, P = 0.002). Analyzed by both univariate and multivariate analysis, the classification schema was independently associated with survival. Using patterns of gene expression based on a defined gene set, ccRCC was classified into two robust subclasses based on inherent molecular features that ultimately correspond to marked differences in clinical outcome. This classification schema thus provides a molecular stratification applicable to individual tumors that has implications to influence treatment decisions, define biological mechanisms involved in ccRCC tumor progression, and direct future drug discovery.

High expression of lymphocyte-associated genes in node-negative HER2+ breast cancers correlates with lower recurrence rates.

Gene expression analysis has identified biologically relevant subclasses of breast cancer. However, most classification schemes do not robustly cluster all HER2+ breast cancers, in part due to limitations and bias of clustering techniques used. In this article, we propose an alternative approach that first separates the HER2+ tumors using a gene amplification signal for Her2/neu amplicon genes and then applies consensus ensemble clustering separately to the HER2+ and HER2- clusters to look for further substructure. We applied this procedure to a microarray data set of 286 early-stage breast cancers treated only with surgery and radiation and identified two basal and four luminal subtypes in the HER2- tumors, as well as two novel and robust HER2+ subtypes. HER2+ subtypes had median distant metastasis-free survival of 99 months [95% confidence interval (95% CI), 83-118 months] and 33 months (95% CI, 11-54 months), respectively, and recurrence rates of 11% and 58%, respectively. The low recurrence subtype had a strong relative overexpression of lymphocyte-associated genes and was also associated with a prominent lymphocytic infiltration on histologic analysis. These data suggest that early-stage HER2+ cancers associated with lymphocytic infiltration are a biologically distinct subtype with an improved natural history.

Computerized Image-Based Detection and Grading of Lymphocytic Infiltration in HER2+ Breast Cancer Histopathology

The identification of phenotypic changes in breast cancer (BC) histopathology on account of corresponding molecular changes is of significant clinical importance in predicting disease outcome. One such example is the presence of lymphocytic infiltration (LI) in histopathology, which has been correlated with nodal metastasis and distant recurrence in HER2+ BC patients. In this paper, we present a computer-aided diagnosis (CADx) scheme to automatically detect and grade the extent of LI in digitized HER2+ BC histopathology. Lymphocytes are first automatically detected by a combination of region growing and Markov random field algorithms. Using the centers of individual detected lymphocytes as vertices, three graphs (Voronoi diagram, Delaunay triangulation, and minimum spanning tree) are constructed and a total of 50 image-derived features describing the arrangement of the lymphocytes are extracted from each sample. A nonlinear dimensionality reduction scheme, graph embedding (GE), is then used to project the high-dimensional feature vector into a reduced 3-D embedding space. A support vector machine classifier is used to discriminate samples with high and low LI in the reduced dimensional embedding space. A total of 41 HER2+ hematoxylin-and-eosin-stained images obtained from 12 patients were considered in this study. For more than 100 three-fold cross-validation trials, the architectural feature set successfully distinguished samples of high and low LI levels with a classification accuracy greater than 90%. The popular unsupervised Varma-Zisserman texton-based classification scheme was used for comparison and yielded a classification accuracy of only 60%. Additionally, the projection of the 50 image-derived features for all 41 tissue samples into a reduced dimensional space via GE allowed for the visualization of a smooth manifold that revealed a continuum between low, intermediate, and high levels of LI. Since it is known that extent of LI in BC biopsy specimens is a prognostic indicator, our CADx scheme will potentially help clinicians determine disease outcome and allow them to make better therapy recommendations for patients with HER2+ BC.

Patterns of Evolution and Host Gene Mimicry in Influenza and Other RNA Viruses

It is well known that the dinucleotide CpG is under-represented in the genomic DNA of many vertebrates. This is commonly thought to be due to the methylation of cytosine residues in this dinucleotide and the corresponding high rate of deamination of 5-methycytosine, which lowers the frequency of this dinucleotide in DNA. Surprisingly, many single-stranded RNA viruses that replicate in these vertebrate hosts also have a very low presence of CpG dinucleotides in their genomes. Viruses are obligate intracellular parasites and the evolution of a virus is inexorably linked to the nature and fate of its host. One therefore expects that virus and host genomes should have common features. In this work, we compare evolutionary patterns in the genomes of ssRNA viruses and their hosts. In particular, we have analyzed dinucleotide patterns and found that the same patterns are pervasively over- or under-represented in many RNA viruses and their hosts suggesting that many RNA viruses evolve by mimicking some of the features of their hosts genes (DNA) and likely also their corresponding mRNAs. When a virus crosses a species barrier into a different host, the pressure to replicate, survive and adapt, leaves a footprint in dinucleotide frequencies. For instance, since human genes seem to be under higher pressure to eliminate CpG dinucleotide motifs than avian genes, this pressure might be reflected in the genomes of human viruses (DNA and RNA viruses) when compared to those of the same viruses replicating in avian hosts. To test this idea we have analyzed the evolution of the influenza virus since 1918. We find that the influenza A virus, which originated from an avian reservoir and has been replicating in humans over many generations, evolves in a direction strongly selected to reduce the frequency of CpG dinucleotides in its genome. Consistent with this observation, we find that the influenza B virus, which has spent much more time in the human population, has adapted to its human host and exhibits an extremely low CpG dinucleotide content. We believe that these observations directly show that the evolution of RNA viral genomes can be shaped by pressures observed in the host genome. As a possible explanation, we suggest that the strong selection pressures acting on these RNA viruses are most likely related to the innate immune response and to nucleotide motifs in the host DNA and RNAs.

Expectation–Maximization-Driven Geodesic Active Contour With Overlap Resolution (EMaGACOR): Application to Lymphocyte Segmentation on Breast Cancer Histopathology

The presence of lymphocytic infiltration (LI) has been correlated with nodal metastasis and tumor recurrence in HER2+ breast cancer (BC). The ability to automatically detect and quantify extent of LI on histopathology imagery could potentially result in the development of an image based prognostic tool for human epidermal growth factor receptor-2 (HER2+) BC patients. Lymphocyte segmentation in hematoxylin and eosin (H&E) stained BC histopathology images is complicated by the similarity in appearance between lymphocyte nuclei and other structures (e.g., cancer nuclei) in the image. Additional challenges include biological variability, histological artifacts, and high prevalence of overlapping objects. Although active contours are widely employed in image segmentation, they are limited in their ability to segment overlapping objects and are sensitive to initialization. In this paper, we present a new segmentation scheme, expectation-maximization (EM) driven geodesic active contour with overlap resolution (EMaGACOR), which we apply to automatically detecting and segmenting lymphocytes on HER2+ BC histopathology images. EMaGACOR utilizes the expectation-maximization algorithm for automatically initializing a geodesic active contour (GAC) and includes a novel scheme based on heuristic splitting of contours via identification of high concavity points for resolving overlapping structures. EMaGACOR was evaluated on a total of 100 HER2+ breast biopsy histology images and was found to have a detection sensitivity of over 86% and a positive predictive value of over 64%. By comparison, the EMaGAC model (without overlap resolution) and GAC model yielded corresponding detection sensitivities of 42% and 19%, respectively. Furthermore, EMaGACOR was able to correctly resolve over 90% of overlaps between intersecting lymphocytes. Hausdorff distance (HD) and mean absolute distance (MAD) for EMaGACOR were found to be 2.1 and 0.9 pixels, respectively, and significantly better compared to the corresponding performance of the EMaGAC and GAC models. EMaGACOR is an efficient, robust, reproducible, and accurate segmentation technique that could potentially be applied to other biomedical image analysis problems.

Neurite Sprouting and Synapse Deterioration in the Aging Caenorhabditis elegans Nervous System

Caenorhabditis elegans is a powerful model for analysis of the conserved mechanisms that modulate healthy aging. In the aging nematode nervous system, neuronal death and/or detectable loss of processes are not readily apparent, but because dendrite restructuring and loss of synaptic integrity are hypothesized to contribute to human brain decline and dysfunction, we combined fluorescence microscopy and electron microscopy (EM) to screen at high resolution for nervous system changes. We report two major components of morphological change in the aging C. elegans nervous system: (1) accumulation of novel outgrowths from specific neurons, and (2) physical decline in synaptic integrity. Novel outgrowth phenotypes, including branching from the main dendrite or new growth from somata, appear at a high frequency in some aging neurons, but not all. Mitochondria are often associated with age-associated branch sites. Lowered insulin signaling confers some maintenance of ALM and PLM neuron structural integrity into old age, and both DAF-16/FOXO and heat shock factor transcription factor HSF-1 exert neuroprotective functions. hsf-1 can act cell autonomously in this capacity. EM evaluation in synapse-rich regions reveals a striking decline in synaptic vesicle numbers and a diminution of presynaptic density size. Interestingly, old animals that maintain locomotory prowess exhibit less synaptic decline than same-age decrepit animals, suggesting that synaptic integrity correlates with locomotory healthspan. Our data reveal similarities between the aging C. elegans nervous system and mammalian brain, suggesting conserved neuronal responses to age. Dissection of neuronal aging mechanisms in C. elegans may thus influence the development of brain healthspan-extending therapies.