Ruchir Shah

The Transcription Factors Snail and Slug Activate the Transforming Growth Factor-Beta Signaling Pathway in Breast Cancer

The transcriptional repressors Snail and Slug are situated at the core of several signaling pathways proposed to mediate epithelial to mesenchymal transition or EMT, which has been implicated in tumor metastasis. EMT involves an alteration from an organized, epithelial cell structure to a mesenchymal, invasive and migratory phenotype. In order to obtain a global view of the impact of Snail and Slug expression, we performed a microarray experiment using the MCF-7 breast cancer cell line, which does not express detectable levels of Snail or Slug. MCF-7 cells were infected with Snail, Slug or control adenovirus, and RNA samples isolated at various time points were analyzed across all transcripts. Our analyses indicated that Snail and Slug regulate many genes in common, but also have distinct sets of gene targets. Gene set enrichment analyses indicated that Snail and Slug directed the transcriptome of MCF-7 cells from a luminal towards a more complex pattern that includes many features of the claudin-low breast cancer signature. Of particular interest, genes involved in the TGF-beta signaling pathway are upregulated, while genes responsible for a differentiated morphology are downregulated following Snail or Slug expression. Further we noticed increased histone acetylation at the promoter region of the transforming growth factor beta-receptor II (TGFBR2) gene following Snail or Slug expression. Inhibition of the TGF-beta signaling pathway using selective small-molecule inhibitors following Snail or Slug addition resulted in decreased cell migration with no impact on the repression of cell junction molecules by Snail and Slug. We propose that there are two regulatory modules embedded within EMT: one that involves repression of cell junction molecules, and the other involving cell migration via TGF-beta and/or other pathways.

RNA-Seq Profiling Reveals Novel Hepatic Gene Expression Pattern in Aflatoxin B1 Treated Rats

Deep sequencing was used to investigate the subchronic effects of 1 ppm aflatoxin B1 (AFB1), a potent hepatocarcinogen, on the male rat liver transcriptome prior to onset of histopathological lesions or tumors. We hypothesized RNA-Seq would reveal more differentially expressed genes (DEG) than microarray analysis, including low copy and novel transcripts related to AFB1’s carcinogenic activity compared to feed controls (CTRL). Paired-end reads were mapped to the rat genome (Rn4) with TopHat and further analyzed by DESeq and Cufflinks-Cuffdiff pipelines to identify differentially expressed transcripts, new exons and unannotated transcripts. PCA and cluster analysis of DEGs showed clear separation between AFB1 and CTRL treatments and concordance among group replicates. qPCR of eight high and medium DEGs and three low DEGs showed good comparability among RNA-Seq and microarray transcripts. DESeq analysis identified 1,026 differentially expressed transcripts at greater than two-fold change (p<0.005) compared to 626 transcripts by microarray due to base pair resolution of transcripts by RNA-Seq, probe placement within transcripts or an absence of probes to detect novel transcripts, splice variants and exons. Pathway analysis among DEGs revealed signaling of Ahr, Nrf2, GSH, xenobiotic, cell cycle, extracellular matrix, and cell differentiation networks consistent with pathways leading to AFB1 carcinogenesis, including almost 200 upregulated transcripts controlled by E2f1-related pathways related to kinetochore structure, mitotic spindle assembly and tissue remodeling. We report 49 novel, differentially-expressed transcripts including confirmation by PCR-cloning of two unique, unannotated, hepatic AFB1-responsive transcripts (HAfT’s) on chromosomes 1.q55 and 15.q11, overexpressed by 10 to 25-fold. Several potentially novel exons were found and exon refinements were made including AFB1 exon-specific induction of homologous family members, Ugt1a6 and Ugt1a7c. We find the rat transcriptome contains many previously unidentified, AFB1-responsive exons and transcripts supporting RNA-Seq’s capabilities to provide new insights into AFB1-mediated gene expression leading to hepatocellular carcinoma.

Analysis of Chromatin Dynamics during Glucocorticoid Receptor Activation

ABSTRACT Steroid hormone receptors initiate a genetic program tightly regulated by the chromatin environment of the responsive regions. Using the glucocorticoid receptor (GR) as a model factor for transcriptional initiation, we classified chromatin structure through formaldehyde-assisted isolation of regulatory elements (FAIRE). We looked at dynamic changes in FAIRE signals during GR activation specifically at regions of receptor interaction. We found a distribution of GR-responsive regions with diverse responses to activation and chromatin modulation. The majority of GR binding regions demonstrate increases in FAIRE signal in response to ligand. However, the majority GR-responsive regions shared a similar FAIRE signal in the basal chromatin state, suggesting a common chromatin structure for GR recruitment. Supporting this notion, global FAIRE sequencing (seq) data indicated an enrichment of signal surrounding the GR binding site prior to activation. Brg-1 knockdown showed response element-specific effects of ATPase-dependent chromatin remodeling. FAIRE induction was universally decreased by Brg-1 depletion, but to varying degrees in a target specific manner. Taken together, these data suggest classes of nuclear receptor response regions that react to activation through different chromatin regulatory events and identify a chromatin structure that classifies the majority of response elements tested.

Abstract 3928: Inhibition of the ubiquitin proteasome system differentially regulates glucocorticoid receptor-mediated transcriptional processes

Proceedings: AACR 103rd Annual Meeting 2012‐‐ Mar 31‐Apr 4, 2012; Chicago, ILnnThe ubiquitin-proteasome system (UPS) plays an integral role in controlling protein turnover to modulate cellular processes such as growth, proliferation, differentiation, and apoptosis. This occurs through multiple mechanisms that rely on the regulation of numerous transcription factors such as steroid hormone nuclear receptors. The glucocorticoid receptor (GR) is a member of this receptor family and its activation promotes nuclear translocation and downstream transcriptional activity through the recruitment of co-regulatory complexes. Our goal is to understand the mechanisms by which the proteasome regulates glucocorticoid receptor-mediated gene transcription. To explore the effects of proteasome inhibition upon GR-mediated transcription, we used ChIP-Seq and microarray studies using breast cancer cells to examine GR binding and gene expression profiles after treatment with dexamethasone (dex) or dex in combination with a proteasome inhibitor, MG132 (dex + mg), for 1 hour and 18 hours. More than 14000 and 5000 GR-binding sites were apparent after both treatments in the 1 hour and 18 hour samples, respectively. GR binding sites seen during 1 hour of treatment can be mapped to more than 9000 genes while the binding sites in the 18 hour treatment samples mapped to over 5000 genes. In addition, over 55% of binding sites in all treatment groups are within 50 kilobases of annotated transcription start sites. To verify these putative interactions, we assayed the transcript levels of genes at both time periods and examined the levels of bound GR throughout the gene. The results reveal a differential regulation of genes that can be divided into three classes: dex-responsive, dex + mg-responsive, or dually responsive. Further analyses will reveal if there are clearly defined chromatin states that demarcate each class. These data provide insights into the intricate mechanisms of proteasome-mediated GR regulation and highlight proteasomal-regulated genes involved in numerous diseases and developmental disorders.nnCitation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 3928. doi:1538-7445.AM2012-3928

SWIFT-Review: a text-mining workbench for systematic review

BackgroundThere is growing interest in using machine learning approaches to priority rank studies and reduce human burden in screening literature when conducting systematic reviews. In addition, identifying addressable questions during the problem formulation phase of systematic review can be challenging, especially for topics having a large literature base. Here, we assess the performance of the SWIFT-Review priority ranking algorithm for identifying studies relevant to a given research question. We also explore the use of SWIFT-Review during problem formulation to identify, categorize, and visualize research areas that are data rich/data poor within a large literature corpus.MethodsTwenty case studies, including 15 public data sets, representing a range of complexity and size, were used to assess the priority ranking performance of SWIFT-Review. For each study, seed sets of manually annotated included and excluded titles and abstracts were used for machine training. The remaining references were then ranked for relevance using an algorithm that considers term frequency and latent Dirichlet allocation (LDA) topic modeling. This ranking was evaluated with respect to (1) the number of studies screened in order to identify 95xa0% of known relevant studies and (2) the “Work Saved over Sampling” (WSS) performance metric. To assess SWIFT-Review for use in problem formulation, PubMed literature search results for 171 chemicals implicated as EDCs were uploaded into SWIFT-Review (264,588 studies) and categorized based on evidence stream and health outcome. Patterns of search results were surveyed and visualized using a variety of interactive graphics.ResultsCompared with the reported performance of other tools using the same datasets, the SWIFT-Review ranking procedure obtained the highest scores on 11 out of 15 of the public datasets. Overall, these results suggest that using machine learning to triage documents for screening has the potential to save, on average, more than 50xa0% of the screening effort ordinarily required when using un-ordered document lists. In addition, the tagging and annotation capabilities of SWIFT-Review can be useful during the activities of scoping and problem formulation.ConclusionsText-mining and machine learning software such as SWIFT-Review can be valuable tools to reduce the human screening burden and assist in problem formulation.

Brg1 Enables Rapid Growth of the Early Embryo by Suppressing Genes That Regulate Apoptosis and Cell Growth Arrest.

ABSTRACT SWI/SNF (switching/sucrose nonfermenting)-dependent chromatin remodeling establishes coordinated gene expression programs during development, yet important functional details remain to be elucidated. We show that the Brg1 (Brahma-related gene 1; Smarca4) ATPase is globally expressed at high levels during postimplantation development and its conditional ablation, beginning at gastrulation, results in increased apoptosis, growth retardation, and, ultimately, embryonic death. Global gene expression analysis revealed that genes upregulated in Rosa26CreERT2; Brg1flox/flox embryos (here referred to as Brg1d/d embryos to describe embryos with deletion of the Brg1flox/flox alleles) negatively regulate cell cycle progression and cell growth. In addition, the p53 (Trp53) protein, which is virtually undetectable in early wild-type embryos, accumulated in the Brg1d/d embryos and activated the p53-dependent pathways. Using P19 cells, we show that Brg1 and CHD4 (chromodomain helicase DNA binding protein 4) coordinate to control target gene expression. Both proteins physically interact and show a substantial overlap of binding sites at chromatin-accessible regions adjacent to genes differentially expressed in the Brg1d/d embryos. Specifically, Brg1 deficiency results in reduced levels of the repressive histone H3 lysine K27 trimethylation (H3K27me3) histone mark and an increase in the amount of open chromatin at the regulatory region of the p53 and p21 (Cdkn1a) genes. These results provide insights into the mechanisms by which Brg1 functions, which is in part via the p53 program, to constrain gene expression and facilitate rapid embryonic growth.

Fine-tuning of epigenetic regulation with respect to promoter CpG content in a cell type-specific manner.

Epigenetic regulation of gene expression is fundamental for cell type-specific gene expression. However, integrated comparative transcriptomic and epigenomic analyses in various adult primary differentiated cells remain underrepresented. We generated promoter landscapes of DNA methylation and three important histone methylation marks (H3K4me3, H3K9me2, and H3K27me3) in two primary cell types (B lymphocytes and liver) from adult mice. In line with previous studies, we also observed distinct H3K4me3 patterns at promoters dictated by CpG content in differentiated primary cells. We further explored the distribution of initiating RNA polymerase II and elongating RNA polymerase II across genes within different promoter classes, suggesting different rate-limiting steps at CpG-rich vs. CpG-poor genes. Examination of differentially expressed genes revealed that regulation of tissue-specific genes is closely related to gene function regardless of promoter type. Although repressive chromatin marks displayed differential preference to promoters based on CpG content, we observed fine-tuning of the pattern of association of these marks with specific promoter types in a cell type-specific manner. The distribution of H3K9me2 and H3K27me3, relative to CpG content, differed substantially between the two cell types. Cell-type specific accumulation of repressive chromatin marks was also observed at silent genes in both cell types, suggesting that differentiated primary cells may exhibit cell-type specificity in the distribution of repressive chromatin marks. Epigenetic regulation of gene expression and the association of specific histone marks with promoter sequence classes are fine-tuned in a cell type-specific manner. This unexpected finding underscores the value of extensive study of epigenetic marks across cell and tissue types.

Differential responses to retinoic acid and endocrine disruptor compounds of subpopulations within human embryonic stem cell lines

The heterogeneous nature of stem cells is an important issue in both research and therapeutic use in terms of directing cell lineage differentiation pathways, as well as self-renewal properties. Using flow cytometry we have identified two distinct subpopulations by size, large and small, within cultures of human embryonic stem (hES) cell lines. These two cell populations respond differentially to retinoic acid (RA) differentiation and several endocrine disruptor compounds (EDC). The large cell population responds to retinoic acid differentiation with greater than a 50% reduction in cell number and loss of Oct-4 expression, whereas the number of the small cell population does not change and Oct-4 protein expression is maintained. In addition, four estrogenic compounds altered SSEA-3 expression differentially between the two cell subpopulations changing their ratios relative to each other. Both populations express stem cell markers Oct-4, Nanog, Tra-1-60, Tra-1-80 and SSEA-4, but express low levels of differentiation markers common to the three germ layers. Cloning studies indicate that both populations can revive the parental population. Furthermore, whole genome microarray identified approximately 400 genes with significantly different expression between the two populations (p<0.01). We propose the differential response to RA in these populations is due to differential gene expression of Notch signaling members, CoupTF1 and CoupTF2, chromatin remodeling and histone modifying genes that render the small population resistant to RA differentiation. The findings that hES cells exist as heterogeneous populations with distinct responses to differentiation signals and environmental stimuli will be relevant for their use for drug discovery and disease therapy.

Transcriptome and DNA Methylome Analysis in a Mouse Model of Diet-Induced Obesity Predicts Increased Risk of Colorectal Cancer

SUMMARY Colorectal cancer (CRC) tends to occur at older age; however, CRC incidence rates have been rising sharply among young age groups. The increasing prevalence of obesity is recognized as a major risk, yet the mechanistic underpinnings remain poorly understood. Using a diet-induced obesity mouse model, we identified obesity-associated molecular changes in the colonic epithelium of young and aged mice, and we further investigated whether the changes were reversed after weight loss. Transcriptome analysis indicated that obesity-related colonic cellular metabolic switch favoring long-chain fatty acid oxidation happened in young mice, while obesity-associated downregulation of negative feedback regulators of pro-proliferative signaling pathways occurred in older mice. Strikingly, colonic DNA methylome was pre-programmed by obesity at young age, priming for a tumor-prone gene signature after aging. Furthermore, obesity-related changes were substantially preserved after short-term weight loss, but they were largely reversed after long-term weight loss. We provided mechanistic insights into increased CRC risk in obesity.

Exome Sequencing of Fresh-frozen or Formalin-fixed Paraffin-embedded B6C3F1/N Mouse Hepatocellular Carcinomas Arising Either Spontaneously or due to Chronic Chemical Exposure:

Hepatocellular carcinoma (HCC) is the third leading cause of cancer-related death worldwide; however, the mutational properties of HCC-associated carcinogens remain largely uncharacterized. We hypothesized that mechanisms underlying chemical-induced HCC can be characterized by evaluating the mutational spectra of these tumors. To test this hypothesis, we performed exome sequencing of B6C3F1/N HCCs that arose either spontaneously in vehicle controls (n = 3) or due to chronic exposure to gingko biloba extract (GBE; n = 4) or methyleugenol (MEG; n = 3). Most archived tumor samples are available as formalin-fixed paraffin-embedded (FFPE) blocks, rather than fresh-frozen (FF) samples; hence, exome sequencing from paired FF and FFPE samples was compared. FF and FFPE samples showed 63% to 70% mutation concordance. Multiple known (e.g., Ctnnb1T41A, BrafV637E) and novel (e.g., Erbb4C559S, Card10A700V, and Klf11P358L) mutations in cancer-related genes were identified. The overall mutational burden was greater for MEG than for GBE or spontaneous HCC samples. To characterize the mutagenic mechanisms, we analyzed the mutational spectra in the HCCs according to their trinucleotide motifs. The MEG tumors clustered closest to Catalogue of Somatic Mutations in Cancer signatures 4 and 24, which are, respectively, associated with benzo(a)pyrene- and aflatoxin-induced HCCs in humans. These results establish a novel approach for classifying liver carcinogens and understanding the mechanisms of hepatocellular carcinogenesis.