Search

Supplemental figures and tables reference in the manuscript

Detailed methods including analysis

A large Excel sheet of data and the complete set of scripts used to analyze the data.

Purpose: Colorectal cancer is the third most common cancer worldwide, causing approximately 700,000 deaths each year. The majority of colorectal cancers begin as adenomas. Definitive screening for colorectal adenomas is currently accomplished through colonoscopy but, owing largely to costs and invasiveness, is typically limited to patient groups at higher risk by virtue of age or family history. We sought to determine if blood-based small RNA markers could detect colorectal adenoma.Experimental Design: We applied high-depth small RNA sequencing to plasma from a large (n = 189) cohort of patients, balanced for age, sex, and ancestry. Our analytical methodology allowed for the detection of both microRNAs and other small RNA species. We replicated sequencing results by qPCR on plasma samples from an independent cohort (n = 140).Results: We found several small RNA species with significant associations to colorectal adenoma, including both microRNAs and non-microRNA small RNAs. These associations were robust to correction for patient covariates, including age. Among the adenoma-associated small RNAs, two, a miR-335-5p isoform and an un-annotated small RNA, were validated by qPCR in an independent cohort. A classifier trained on measures of these two RNAs in the discovery cohort yields an AUC of 0.755 (0.775 with age) for adenoma detection in the independent cohort. This classifier accurately detects adenomas in patients under 50 and is robust to sex or ancestry.Conclusions: Circulating small RNAs (including but not limited to miRNAs) discovered by sequencing and validated by qPCR identify patients with colorectal adenomas effectively. Clin Cancer Res; 24(9); 2092–9. ©2018 AACR.

BackgroundTauopathies are a group of neurodegenerative diseases driven by abnormal aggregates of tau, a microtubule associated protein encoded by theMAPTgene.MAPTexpression is absent in neural progenitor cells (NPCs) and increases during differentiation. This temporally dynamic expression pattern suggests thatMAPTexpression is controlled by transcription factors and cis-regulatory elements specific to differentiated cell types. Given the relevance ofMAPTexpression to neurodegeneration pathogenesis, identification of such elements is relevant to understanding genetic risk factors.MethodsWe performed HiC, chromatin conformation capture (Capture-C), single-nucleus multiomics (RNA-seq+ATAC-seq), bulk ATAC-seq, and ChIP-seq for H3K27Ac and CTCF in NPCs and neurons differentiated from human iPSC cultures. We nominated candidate cis-regulatory elements (cCREs) forMAPTin human NPCs, differentiated neurons, and pure cultures of inhibitory and excitatory neurons. We then assayed these cCREs using luciferase assays and CRISPR interference (CRISPRi) experiments to measure their effects onMAPTexpression. Finally, we integrated cCRE annotations into an analysis of genetic variation in AD cases and controls.ResultsUsing orthogonal genomics approaches, we nominated 94 cCREs forMAPT, including the identification of cCREs specifically active in differentiated neurons. Eleven regions enhanced reporter gene transcription in luciferase assays. Using CRISPRi, 5 of the 94 regions tested were identified as necessary forMAPTexpression as measured by RT-qPCR and RNA-seq. Rare and predicted damaging genetic variation in both nominated and confirmed CREs was depleted in AD cases relative to controls (OR = 0.40, p = 0.004), consistent with the hypothesis that variants that disruptMAPTenhancer activity, and thereby reduceMAPTexpression, may be protective against neurodegenerative disease.ConclusionsWe identified both proximal and distal regulatory elements forMAPTand confirmed the regulatory function for several regions, including three regions centromeric toMAPTbeyond the well-described H1/H2 haplotype inversion breakpoint. This study provides compelling evidence for pursuing detailed knowledge of CREs for genes of interest to permit better understanding of disease risk.

Massively parallel reporter assays (MPRAs) are useful tools to characterize regulatory elements in human genomes. An aspect of MPRAs that is not typically the focus of analysis is their intrinsic ability to differentiate activity levels for a given sequence element when placed in both of its possible orientations relative to the reporter construct. Here, we describe pervasive strand asymmetry of MPRA signals in data sets from multiple reporter configurations in both published and newly reported data. These effects are reproducible across different cell types and in different treatments within a cell type and are observed both within and outside of annotated regulatory elements. From elements in gene bodies, MPRA strand asymmetry favors the sense strand, suggesting that function related to endogenous transcription is driving the phenomenon. Similarly, we find that within Alu mobile element insertions, strand asymmetry favors the transcribed strand of the ancestral retrotransposon. The effect is consistent across the multiplicity of Alu elements in human genomes and is more pronounced in less diverged Alu elements. We find sequence features driving MPRA strand asymmetry and show its prediction from sequence alone. We see some evidence for RNA stabilization and transcriptional activation mechanisms and hypothesize that the effect is driven by natural selection favoring efficient transcription. Our results indicate that strand asymmetry is a pervasive and reproducible feature in MPRA data. More importantly, the fact that MPRA asymmetry favors naturally transcribed strands suggests that it stems from preserved biological functions that have a substantial, global impact on gene and genome evolution.

// Joy M. McDaniel 1, 2 , Katherine E. Varley 3 , Jason Gertz 3 , Daniel S. Savic 1 , Brian S. Roberts 1 , Sarah K. Bailey 4 , Lalita A. Shevde 4, 6 , Ryne C. Ramaker 1, 7 , Brittany N. Lasseigne 1 , Marie K. Kirby 1 , Kimberly M. Newberry 1 , E. Christopher Partridge 1 , Angela L. Jones 1 , Braden Boone 1 , Shawn E. Levy 1 , Patsy G. Oliver 5 , Katherine C. Sexton 6 , William E. Grizzle 6 , Andres Forero 6 , Donald J. Buchsbaum 5 , Sara J. Cooper 1 , Richard M. Myers 1 1 HudsonAlpha Institute for Biotechnology, Huntsville, AL 35806, USA 2 The University of Alabama in Huntsville, Huntsville, AL 35899, USA 3 Department of Oncological Sciences, Huntsman Cancer Institute, University of Utah, Salt Lake City, UT, 84112, USA 4 Department of Pathology, University of Alabama at Birmingham, Birmingham, AL, 35294, USA 5 Department of Radiation Oncology, University of Alabama at Birmingham, Birmingham, AL, 35294, USA 6 University of Alabama at Birmingham Comprehensive Cancer Center, Birmingham, AL 35294, USA 7 Department of Genetics, University of Alabama at Birmingham, Birmingham, AL, 35294, USA Correspondence to: Richard M. Myers, email: rmyers@hudsonalpha.org Keywords: TNBC, STAT3, ChIP-seq, RNA-seq, invasion Received: August 15, 2016 Accepted: November 14, 2016 Published: December 24, 2016 ABSTRACT Breast cancer is a heterogeneous disease comprised of four molecular subtypes defined by whether the tumor-originating cells are luminal or basal epithelial cells. Breast cancers arising from the luminal mammary duct often express estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth receptor 2 (HER2). Tumors expressing ER and/or PR are treated with anti-hormonal therapies, while tumors overexpressing HER2 are targeted with monoclonal antibodies. Immunohistochemical detection of ER, PR, and HER2 receptors/proteins is a critical step in breast cancer diagnosis and guided treatment. Breast tumors that do not express these proteins are known as “triple negative breast cancer” (TNBC) and are typically basal-like. TNBCs are the most aggressive subtype, with the highest mortality rates and no targeted therapy, so there is a pressing need to identify important TNBC tumor regulators. The signal transducer and activator of transcription 3 (STAT3) transcription factor has been previously implicated as a constitutively active oncogene in TNBC. However, its direct regulatory gene targets and tumorigenic properties have not been well characterized. By integrating RNA-seq and ChIP-seq data from 2 TNBC tumors and 5 cell lines, we discovered novel gene signatures directly regulated by STAT3 that were enriched for processes involving inflammation, immunity, and invasion in TNBC. Functional analysis revealed that STAT3 has a key role regulating invasion and metastasis, a characteristic often associated with TNBC. Our findings suggest therapies targeting STAT3 may be important for preventing TNBC metastasis.

MicroRNAs (miRNAs) play a critical role as post-transcriptional regulators of gene expression. The ENCODE project profiled the expression of miRNAs in a comprehensive set of tissues during a time-course of mouse embryonic development and captured the expression dynamics of 785 miRNAs. We found distinct tissue and developmental stage specific miRNA expression clusters, with an overall pattern of increasing tissue specific expression as development proceeds. Comparative analysis of conserved miRNAs in mouse and human revealed stronger clustering of expression patterns by tissue types rather than by species. An analysis of messenger RNA gene expression clusters compared with miRNA expression clusters identifies the potential role of specific miRNA expression clusters in suppressing the expression of mRNAs specific to other developmental programs in the tissue where these microRNAs are expressed during embryonic development. Our results provide the most comprehensive timecourse of miRNA expression as an integrated part of the ENCODE reference dataset for mouse embryonic development.

In small RNA (smRNA) sequencing studies, highly abundant molecules such as adapter dimer products and tissue-specific microRNAs (miRNAs) inhibit accurate quantification of lowly expressed species. We previously developed a method to selectively deplete highly abundant miRNAs. However, this method does not deplete adapter dimer ligation products that, unless removed by gel-separation, comprise most of the library. Here, we have adapted and modified recently described methods for CRISPR/Cas9–based Depletion of Abundant Species by Hybridization (‘DASH’) to smRNA-seq, which we have termed miRNA and Adapter Dimer—DASH (MAD-DASH). In MAD-DASH, Cas9 is complexed with single guide RNAs (sgRNAs) targeting adapter dimer ligation products, alongside highly expressed tissue-specific smRNAs, for cleavage in vitro. This process dramatically reduces adapter dimer and targeted smRNA sequences, can be multiplexed, shows minimal off-target effects, improves the quantification of lowly expressed miRNAs from human plasma and tissue derived RNA, and obviates the need for gel-separation, greatly increasing sample throughput. Additionally, the method is fully customizable to other smRNA-seq preparation methods. Like depletion of ribosomal RNA for mRNA-seq and mitochondrial DNA for ATAC-seq, our method allows for greater proportional read-depth of non-targeted sequences.

Summary The wide range of RNA-seq applications and their high-computational needs require the development of pipelines orchestrating the entire workflow and optimizing usage of available computational resources. We present aRNApipe, a project-oriented pipeline for processing of RNA-seq data in high-performance cluster environments. aRNApipe is highly modular and can be easily migrated to any high-performance computing (HPC) environment. The current applications included in aRNApipe combine the essential RNA-seq primary analyses, including quality control metrics, transcript alignment, count generation, transcript fusion identification, alternative splicing and sequence variant calling. aRNApipe is project-oriented and dynamic so users can easily update analyses to include or exclude samples or enable additional processing modules. Workflow parameters are easily set using a single configuration file that provides centralized tracking of all analytical processes. Finally, aRNApipe incorporates interactive web reports for sample tracking and a tool for managing the genome assemblies available to perform an analysis. Availability and documentation https://github.com/HudsonAlpha/aRNAPipe; DOI: 10.5281/zenodo.202950 Supplementary information Supplementary data are available at Bioinformatics online.

SummaryThe wide range of RNA-seq applications and their high computational needs require the development of pipelines orchestrating the entire workflow and optimizing usage of available computational resources. We present aRNApipe, a project-oriented pipeline for processing of RNA-seq data in high performance cluster environments. aRNApipe is highly modular and can be easily migrated to any high performance computing (HPC) environment. The current applications included in aRNApipe combine the essential RNA-seq primary analyses, including quality control metrics, transcript alignment, count generation, transcript fusion identification, alternative splicing, and sequence variant calling. aRNApipe is project-oriented and dynamic so users can easily update analyses to include or exclude samples or enable additional processing modules. Workflow parameters are easily set using a single configuration file that provides centralized tracking of all analytical processes. Finally, aRNApipe incorporates interactive web reports for sample tracking and a tool for managing the genome assemblies available to perform an analysis.Availability and documentationhttps://github.com/HudsonAlpha/aRNAPipe; DOI:10.5281/zenodo.202950Contactrmyers@hudsonalpha.orgSupplementary informationSupplementary data are available at Bioinformatics online.

Background The liver X receptors (LXRs, NR1H2 and NR1H3) and peroxisome proliferator-activated receptor gamma (PPARG, NR1C3) nuclear receptor transcription factors (TFs) are master regulators of energy homeostasis. Intriguingly, recent studies suggest that these metabolic regulators also impact tumor cell proliferation. However, a comprehensive temporal molecular characterization of the LXR and PPARG gene regulatory responses in tumor cells is still lacking. Methods To better define the underlying molecular processes governing the genetic control of cellular growth in response to extracellular metabolic signals, we performed a comprehensive, genome-wide characterization of the temporal regulatory cascades mediated by LXR and PPARG signaling in HT29 colorectal cancer cells. For this analysis, we applied a multi-tiered approach that incorporated cellular phenotypic assays, gene expression profiles, chromatin state dynamics, and nuclear receptor binding patterns. Results Our results illustrate that the activation of both nuclear receptors inhibited cell proliferation and further decreased glutathione levels, consistent with increased cellular oxidative stress. Despite a common metabolic reprogramming, the gene regulatory network programs initiated by these nuclear receptors were widely distinct. PPARG generated a rapid and short-term response while maintaining a gene activator role. By contrast, LXR signaling was prolonged, with initial, predominantly activating functions that transitioned to repressive gene regulatory activities at late time points. Conclusions Through the use of a multi-tiered strategy that integrated various genomic datasets, our data illustrate that distinct gene regulatory programs elicit common phenotypic effects, highlighting the complexity of the genome. These results further provide a detailed molecular map of metabolic reprogramming in cancer cells through LXR and PPARG activation. As ligand-inducible TFs, these nuclear receptors can potentially serve as attractive therapeutic targets for the treatment of various cancers. Electronic supplementary material The online version of this article (doi:10.1186/s13073-016-0328-6) contains supplementary material, which is available to authorized users.

Combination drug therapy is a widely used paradigm for managing numerous human malignancies. In cancer treatment, additive and/or synergistic drug combinations can convert weakly efficacious monotherapies into regimens that produce robust antitumor activity. This can be explained in part through pathway interdependencies that are critical for cancer cell proliferation and survival. However, identification of the various interdependencies is difficult due to the complex molecular circuitry that underlies tumor development and progression. Here, we present a high-throughput platform that allows for an unbiased identification of synergistic and efficacious drug combinations. In a screen of 22,737 experiments of 583 doublet combinations in 39 diverse cancer cell lines using a 4 by 4 dosing regimen, both well-known and novel synergistic and efficacious combinations were identified. Here, we present an example of one such novel combination, a Wee1 inhibitor (AZD1775) and an mTOR inhibitor (ridaforolimus), and demonstrate that the combination potently and synergistically inhibits cancer cell growth in vitro and in vivo. This approach has identified novel combinations that would be difficult to reliably predict based purely on our current understanding of cancer cell biology. Mol Cancer Ther; 15(6); 1155–62. ©2016 AACR.

Chromatin immunoprecipitation followed by next-generation DNA sequencing (ChIP-seq) is a widely used technique for identifying transcription factor (TF) binding events throughout an entire genome. However, ChIP-seq is limited by the availability of suitable ChIP-seq grade antibodies, and the vast majority of commercially available antibodies fail to generate usable data sets. To ameliorate these technical obstacles, we present a robust methodological approach for performing ChIP-seq through epitope tagging of endogenous TFs. We used clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9-based genome editing technology to develop CRISPR epitope tagging ChIP-seq (CETCh-seq) of DNA-binding proteins. We assessed the feasibility of CETCh-seq by tagging several DNA-binding proteins spanning a wide range of endogenous expression levels in the hepatocellular carcinoma cell line HepG2. Our data exhibit strong correlations between both replicate types as well as with standard ChIP-seq approaches that use TF antibodies. Notably, we also observed minimal changes to the cellular transcriptome and to the expression of the tagged TF. To examine the robustness of our technique, we further performed CETCh-seq in the breast adenocarcinoma cell line MCF7 as well as mouse embryonic stem cells and observed similarly high correlations. Collectively, these data highlight the applicability of CETCh-seq to accurately define the genome-wide binding profiles of DNA-binding proteins, allowing for a straightforward methodology to potentially assay the complete repertoire of TFs, including the large fraction for which ChIP-quality antibodies are not available.

Highly abundant microRNAs (miRNAs) in small RNA sequencing libraries make it difficult to obtain efficient measurements of more lowly expressed species. We present a new method that allows for the selective blocking of specific, abundant miRNAs during preparation of sequencing libraries. This technique is specific with little off-target effects and has no impact on the reproducibility of the measurement of non-targeted species. In human plasma samples, we demonstrate that blocking of highly abundant hsa-miR-16-5p leads to improved detection of lowly expressed miRNAs and more precise measurement of differential expression overall. Furthermore, we establish the ability to target a second abundant miRNA and to multiplex the blocking of two miRNAs simultaneously. For small RNA sequencing, this technique could fill a similar role as do ribosomal or globin removal technologies in messenger RNA sequencing.

Transcription factors (TFs) bind to thousands of DNA sequences in mammalian genomes, but most of these binding events appear to have no direct effect on gene expression. It is unclear why only a subset of TF bound sites are actively involved in transcriptional regulation. Moreover, the key genomic features that accurately discriminate between active and inactive TF binding events remain ambiguous. Recent studies have identified promoter-distal RNA polymerase II (RNAP2) binding at enhancer elements, suggesting that these interactions may serve as a marker for active regulatory sequences. Despite these correlative analyses, a thorough functional validation of these genomic co-occupancies is still lacking. To characterize the gene regulatory activity of DNA sequences underlying promoter-distal TF binding events that co-occur with RNAP2 and TF sites devoid of RNAP2 occupancy using a functional reporter assay, we performed cis-regulatory element sequencing (CRE-seq). We tested more than 1000 promoter-distal CCAAT/enhancer-binding protein beta (CEBPB)-bound sites in HepG2 and K562 cells, and found that CEBPB-bound sites co-occurring with RNAP2 were more likely to exhibit enhancer activity. CEBPB-bound sites further maintained substantial cell-type specificity, indicating that local DNA sequence can accurately convey cell-type–specific regulatory information. By comparing our CRE-seq results to a comprehensive set of genome annotations, we identified a variety of genomic features that are strong predictors of regulatory element activity and cell-type–specific activity. Collectively, our functional assay results indicate that RNAP2 occupancy can be used as a key genomic marker that can distinguish active from inactive TF bound sites.

Integrative ChIP-seq/Microarray Analysis Identifies a CTNNB1 Target Signature Enriched in Intestinal Stem Cells and Colon Cancer Kazuhide Watanabe 1 , Jacob Biesinger 2,3 , Michael L. Salmans 1 , Brian S. Roberts 4 , William T. Arthur 4 , Michele Cleary 4 , Bogi Andersen 1 , Xiaohui Xie 2,3 , Xing Dai 1,2 * 1 Department of Biological Chemistry, School of Medicine, University of California Irvine, Irvine, California, United States of America, 2 Institute for Genomics and Bioinformatics, University of California Irvine, Irvine, California, United States of America, 3 Department of Computer Science, University of California Irvine, Irvine, California, United States of America, 4 Rosetta Inpharmatics, LLC, Merck & Co Inc., Seattle, Washington, United States of America Abstract Background: Deregulation of canonical Wnt/CTNNB1 (beta-catenin) pathway is one of the earliest events in the pathogenesis of colon cancer. Mutations in APC or CTNNB1 are highly frequent in colon cancer and cause aberrant stabilization of CTNNB1, which activates the transcription of Wnt target genes by binding to chromatin via the TCF/LEF transcription factors. Here we report an integrative analysis of genome-wide chromatin occupancy of CTNNB1 by chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-seq) and gene expression profiling by microarray analysis upon RNAi-mediated knockdown of CTNNB1 in colon cancer cells. Results: We observed 3629 CTNNB1 binding peaks across the genome and a significant correlation between CTNNB1 binding and knockdown-induced gene expression change. Our integrative analysis led to the discovery of a direct Wnt target signature composed of 162 genes. Gene ontology analysis of this signature revealed a significant enrichment of Wnt pathway genes, suggesting multiple feedback regulations of the pathway. We provide evidence that this gene signature partially overlaps with the Lgr5 + intestinal stem cell signature, and is significantly enriched in normal intestinal stem cells as well as in clinical colorectal cancer samples. Interestingly, while the expression of the CTNNB1 target gene set does not correlate with survival, elevated expression of negative feedback regulators within the signature predicts better prognosis. Conclusion: Our data provide a genome-wide view of chromatin occupancy and gene regulation of Wnt/CTNNB1 signaling in colon cancer cells. Citation: Watanabe K, Biesinger J, Salmans ML, Roberts BS, Arthur WT, et al. (2014) Integrative ChIP-seq/Microarray Analysis Identifies a CTNNB1 Target Signature Enriched in Intestinal Stem Cells and Colon Cancer. PLoS ONE 9(3): e92317. doi:10.1371/journal.pone.0092317 Editor: Ted S. Acott, Casey Eye Institute, United States of America Received October 30, 2013; Accepted February 20, 2014; Published March 20, 2014 Copyright: s 2014 Watanabe et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Funding: Funding from Susan G. Komen grant KG110897 (XD), a U.S. Department of Defense BCRP Postdoctoral Fellowship (K.W. W81XWH-10-1-0383), and National Institutes of Health grant R01HG006870 (to XX). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing Interests: Brian S. Roberts and William T. Arthur are former employees of Rosetta Inpharmatics, LLC, Merck & Co Inc. Michele Cleary is currently an employee of Merck & Co., Inc. There are no patents, products in development or marketed products to declare. This does not alter the authors’ adherence to all the PLOS ONE policies on sharing data and materials, as detailed online in the guide for authors. * E-mail: xdai@uci.edu leads to abnormal proliferation of mutated cells through activation of target genes such as MYC and CCND1 [7]. However, given the divergent cellular roles of Wnt/CTNNB1 signaling, other target genes may also contribute to the pathogenesis of CRCs. Thus, identification and characterization of CTNNB1 target genes genome-wide have been an important pursuit in CRC studies. Transcriptional profiling was used to identify Wnt/ CTNNB1 target genes in colon cancer cells [8]. However, analysis of gene expression alone is limited due to inability to distinguish between primary and secondary effects of Wnt pathway activation. More recently, chromatin immunoprecipitation (ChIP) followed by large-scale DNA analysis such as DNA-chip (ChIP-chip) or high-throughput sequencing (ChIP-seq) [9,10] was used to identify the genomic loci to which CTNNB1 or TCF factors directly bind [11–14]. However, ChIP-based studies of various transcription factors suggest that not all binding events identified correlate with Background Wnt/CTNNB1 signaling is a conserved pathway that plays fundamental roles in embryonic development, tissue homeostasis and maintenance of stem cells. In normal intestine, this pathway is essential for the development and maintenance of intestinal stem cells [1,2]. Activation of canonical Wnt signaling involves stabilization of cytoplasmic CTNNB1, which is otherwise degrad- ed by the proteasome through a degradation complex composed of tumor suppressor protein APC, serine/threonine kinase GSK-3, and Axin. Stabilized CTNNB1 translocates into the nucleus where it binds to TCF/LEF transcription factors and activates target gene expression [2]. A key initiation event of colorectal cancers (CRCs) is CTNNB1 stabilization through loss of the APC gene or activating mutations in the CTNNB1 gene [3]. This genetic event occurs primarily in the intestinal stem cell population [4–6] and PLOS ONE | www.plosone.org March 2014 | Volume 9 | Issue 3 | e92317

Background Inhibition of kinases involved in the DNA damage response sensitizes cells to genotoxic agents by abrogating checkpoint-induced cell cycle arrest. CHK1 and WEE1 act in a pathway upstream of CDK1 to inhibit cell cycle progression in response to damaged DNA. Therapeutic targeting of either CHK1 or WEE1, in combination with chemotherapy, is under clinical evaluation. These studies examine the overlap and potential for synergy when CHK1 and WEE1 are inhibited in cancer cell models. Methods Small molecules MK-8776 and MK-1775 were used to selectively and potently inhibit CHK1 and WEE1, respectively. Results In vitro, the combination of MK-8776 and MK-1775 induces up to 50-fold more DNA damage than either MK-8776 or MK-1775 alone at a fixed concentration. This requires aberrant cyclin-dependent kinase activity but does not appear to be dependent on p53 status alone. Furthermore, DNA damage takes place primarily in S-phase cells, implying disrupted DNA replication. When dosed together, the combination of MK-8776 and MK-1775 induced more intense and more durable DNA damage as well as anti-tumor efficacy than either MK-8776 or MK-1775 dosed alone. DNA damage induced by the combination was detected in up to 40% of cells in a treated xenograft tumor model. Conclusions These results highlight the roles of WEE1 and CHK1 in maintaining genomic integrity. Importantly, the strong synergy observed upon inhibition of both kinases suggests unique yet complimentary anti-tumor effects of WEE1 and CHK1 inhibition. This demonstration of DNA double strand breaks in the absence of a DNA damaging chemotherapeutic provides preclinical rationale for combining WEE1 and CHK1 inhibitors as a cancer treatment regimen.

Although we have made great progress in understanding the complex genetic alterations that underlie human cancer, it has proven difficult to identify which molecularly targeted therapeutics will benefit which patients. Drug-specific modulation of oncogenic signaling pathways in specific patient subpopulations can predict responsiveness to targeted therapy. Here, we report a pathway-based phosphoprofiling approach to identify and quantify clinically relevant, drug-specific biomarkers for phosphatidylinositol 3-kinase (PI3K) pathway inhibitors that target AKT, phosphoinositide-dependent kinase 1 (PDK1), and PI3K-mammalian target of rapamycin (mTOR). We quantified 375 nonredundant PI3K pathway-relevant phosphopeptides, all containing AKT, PDK1, or mitogen-activated protein kinase substrate recognition motifs. Of these phosphopeptides, 71 were drug-regulated, 11 of them by all three inhibitors. Drug-modulated phosphoproteins were enriched for involvement in cytoskeletal reorganization (filamin, stathmin, dynamin, PAK4, and PTPN14), vesicle transport (LARP1, VPS13D, and SLC20A1), and protein translation (S6RP and PRAS40). We then generated phosphospecific antibodies against selected, drug-regulated phosphorylation sites that would be suitable as biomarker tools for PI3K pathway inhibitors. As proof of concept, we show clinical translation feasibility for an antibody against phospho-PRAS40(Thr246). Evaluation of binding of this antibody in human cancer cell lines, a PTEN (phosphatase and tensin homolog deleted from chromosome 10)-deficient mouse prostate tumor model, and triple-negative breast tumor tissues showed that phospho-PRAS40(Thr246) positively correlates with PI3K pathway activation and predicts AKT inhibitor sensitivity. In contrast to phosphorylation of AKT(Thr308), the phospho-PRAS40(Thr246) epitope is highly stable in tissue samples and thus is ideal for immunohistochemistry. In summary, our study illustrates a rational approach for discovery of drug-specific biomarkers toward development of patient-tailored treatments.

Background Hyperactivation of the Ras signaling pathway is a driver of many cancers, and RAS pathway activation can predict response to targeted therapies. Therefore, optimal methods for measuring Ras pathway activation are critical. The main focus of our work was to develop a gene expression signature that is predictive of RAS pathway dependence. Methods We used the coherent expression of RAS pathway-related genes across multiple datasets to derive a RAS pathway gene expression signature and generate RAS pathway activation scores in pre-clinical cancer models and human tumors. We then related this signature to KRAS mutation status and drug response data in pre-clinical and clinical datasets. Results The RAS signature score is predictive of KRAS mutation status in lung tumors and cell lines with high (> 90%) sensitivity but relatively low (50%) specificity due to samples that have apparent RAS pathway activation in the absence of a KRAS mutation. In lung and breast cancer cell line panels, the RAS pathway signature score correlates with pMEK and pERK expression, and predicts resistance to AKT inhibition and sensitivity to MEK inhibition within both KRAS mutant and KRAS wild-type groups. The RAS pathway signature is upregulated in breast cancer cell lines that have acquired resistance to AKT inhibition, and is downregulated by inhibition of MEK. In lung cancer cell lines knockdown of KRAS using siRNA demonstrates that the RAS pathway signature is a better measure of dependence on RAS compared to KRAS mutation status. In human tumors, the RAS pathway signature is elevated in ER negative breast tumors and lung adenocarcinomas, and predicts resistance to cetuximab in metastatic colorectal cancer. Conclusions These data demonstrate that the RAS pathway signature is superior to KRAS mutation status for the prediction of dependence on RAS signaling, can predict response to PI3K and RAS pathway inhibitors, and is likely to have the most clinical utility in lung and breast tumors.

High-throughput siRNA screens are now widely used for identifying novel drug targets and mapping disease pathways. Despite their popularity, there remain challenges related to data variability, primarily due to measurement errors, biological variance, uneven transfection efficiency, the efficacy of siRNA sequences, or off-target effects, and consequent high false discovery rates. Data variability can be reduced if siRNA screens are performed in replicate. Running a large-scale siRNA screen in replicate is difficult, however, because of the technical challenges related to automating complicated steps of siRNA transfection, often with multiplexed assay readouts, and controlling environmental humidity during long incubation periods. Small-molecule screens have greatly benefited in the past decade from assay miniaturization to high-density plates such that 1,536-well nanoplate screenings are now a routine process, allowing fast, efficient, and affordable operations without compromising underlying biology or important assay characteristics. Here, we describe the development of a 1,536-well nanoplate siRNA transfection protocol that utilizes the instruments commonly found in small-molecule high throughput screening laboratories. This protocol was then successfully demonstrated in a triplicate large-scale siRNA screen for the identification of regulators of the Wnt/beta-catenin pathway.

The multi-protein beta-catenin destruction complex tightly regulates beta-catenin protein levels by shuttling beta-catenin to the proteasome. Glycogen synthase kinase 3beta (GSK3beta), a key serine/threonine kinase in the destruction complex, is responsible for several phosphorylation events that mark beta-catenin for ubiquitination and subsequent degradation. Because modulation of both beta-catenin and GSK3beta activity may have important implications for treating disease, a complete understanding of the mechanisms that regulate the beta-catenin/GSK3beta interaction is warranted. We screened an arrayed lentivirus library expressing small hairpin RNAs (shRNAs) targeting 5,201 human druggable genes for silencing events that activate a beta-catenin pathway reporter (BAR) in synergy with 6-bromoindirubin-3'oxime (BIO), a specific inhibitor of GSK3beta. Top screen hits included shRNAs targeting dihydrofolate reductase (DHFR), the target of the anti-inflammatory compound methotrexate. Exposure of cells to BIO plus methotrexate resulted in potent synergistic activation of BAR activity, reduction of beta-catenin phosphorylation at GSK3-specific sites, and accumulation of nuclear beta-catenin. Furthermore, the observed synergy correlated with inhibitory phosphorylation of GSK3beta and was neutralized upon inhibition of phosphatidyl inositol 3-kinase (PI3K). Linking these observations to inflammation, we also observed synergistic inhibition of lipopolysaccharide (LPS)-induced production of pro-inflammatory cytokines (TNFalpha, IL-6, and IL-12), and increased production of the anti-inflammatory cytokine IL-10 in peripheral blood mononuclear cells exposed to GSK3 inhibitors and methotrexate. Our data establish DHFR as a novel modulator of beta-catenin and GSK3 signaling and raise several implications for clinical use of combined methotrexate and GSK3 inhibitors as treatment for inflammatory disease.

Wnts are secreted ligands that activate several receptor-mediated signal transduction cascades. Homeostatic Wnt signaling through beta-catenin is required in adults, because either elevation or attenuation of beta-catenin function has been linked to diverse diseases. To contribute to the identification of both protein and pharmacological regulators of this pathway, we describe a combinatorial screen that merged data from a high-throughput screen of known bioactive compounds with an independent focused small interfering RNA screen. Each screen independently revealed Bruton's tyrosine kinase (BTK) as an inhibitor of Wnt-beta-catenin signaling. Loss of BTK function in human colorectal cancer cells, human B cells, zebrafish embryos, and cells derived from X-linked agammaglobulinemia patients with a mutant BTK gene resulted in elevated Wnt-beta-catenin signaling, confirming that BTK acts as a negative regulator of this pathway. From affinity purification-mass spectrometry and biochemical binding studies, we found that BTK directly interacts with a nuclear component of Wnt-beta-catenin signaling, CDC73. Further, we show that BTK increased the abundance of CDC73 in the absence of stimulation and that CDC73 acted as a repressor of beta-catenin-mediated transcription in human colorectal cancer cells and B cells.

The identification and characterization of previously unidentified signal transduction molecules has expanded our understanding of biological systems and facilitated the development of mechanism-based therapeutics. We present a highly validated small interfering RNA (siRNA) screen that functionally annotates the human genome for modulation of the Wnt/beta-catenin signal transduction pathway. Merging these functional data with an extensive Wnt/beta-catenin protein interaction network produces an integrated physical and functional map of the pathway. The power of this approach is illustrated by the positioning of siRNA screen hits into discrete physical complexes of proteins. Similarly, this approach allows one to filter discoveries made through protein-protein interaction screens for functional contribution to the phenotype of interest. Using this methodology, we characterized AGGF1 as a nuclear chromatin-associated protein that participates in beta-catenin-mediated transcription in human colon cancer cells.

MicroRNAs (miRNAs) are a class of small noncoding RNAs that have important regulatory roles in multicellular organisms. The public miRNA database contains 321 human miRNA sequences, 234 of which have been experimentally verified. To explore the possibility that additional miRNAs are present in the human genome, we have developed an experimental approach called miRNA serial analysis of gene expression (miRAGE) and used it to perform the largest experimental analysis of human miRNAs to date. Sequence analysis of 273,966 small RNA tags from human colorectal cells allowed us to identify 200 known mature miRNAs, 133 novel miRNA candidates, and 112 previously uncharacterized miRNA* forms. To aid in the evaluation of candidate miRNAs, we disrupted the Dicer locus in three human colorectal cancer cell lines and examined known and novel miRNAs in these cells. These studies suggest that the human genome contains many more miRNAs than currently identified and provide an approach for the large-scale experimental cloning of novel human miRNAs in human tissues.

There has been a surge of interest in the biology of microRNAs and the technology of RNA interference. We describe a simple, robust, inexpensive assay for quantitative analysis of microRNAs and short-interfering RNAs. The method relies on primer extension conversion of RNA to cDNA by reverse transcription followed by quantitative, real-time PCR. Technical parameters critical to the success of the assay are presented. Measurements of microRNA levels are sensitive, with most assays allowing measurements in the femtomolar range, which corresponds to tens of copies per cell or less. The assay has a high dynamic range and provides linear readout over differences in microRNA concentrations that span 6–7 orders of magnitude. The assay is capable of discriminating between related microRNA family members that differ by subtle sequence differences. We used the method for quantitative analysis of six microRNAs across 12 tissue samples. The data confirm striking variation in the patterns of expression of these noncoding regulatory RNAs.

Objectives Human airway epithelial cells are the principal target of human rhinovirus (HRV), a common cold pathogen that triggers the majority of asthma exacerbations. The objectives of this study were 1) to evaluate an in vitro air liquid interface cultured human airway epithelial cell model for HRV infection, and 2) to identify gene expression patterns associated with asthma intrinsically and/or after HRV infection using this model. Methods Air-liquid interface (ALI) human airway epithelial cell cultures were prepared from 6 asthmatic and 6 non-asthmatic donors. The effects of rhinovirus RV-A16 on ALI cultures were compared. Genome-wide gene expression changes in ALI cultures following HRV infection at 24 hours post exposure were further analyzed using RNA-seq technology. Cellular gene expression and cytokine/chemokine secretion were further evaluated by qPCR and a Luminex-based protein assay, respectively. Main Results ALI cultures were readily infected by HRV. RNA-seq analysis of HRV infected ALI cultures identified sets of genes associated with asthma specific viral responses. These genes are related to inflammatory pathways, epithelial structure and remodeling and cilium assembly and function, including those described previously (e.g. CCL5, CXCL10 and CX3CL1, MUC5AC, CDHR3), and novel ones that were identified for the first time in this study (e.g. CCRL1). Conclusions ALI-cultured human airway epithelial cells challenged with HRV are a useful translational model for the study of HRV-induced responses in airway epithelial cells, given that gene expression profile using this model largely recapitulates some important patterns of gene responses in patients during clinical HRV infection. Furthermore, our data emphasize that both abnormal airway epithelial structure and inflammatory signaling are two important asthma signatures, which can be further exacerbated by HRV infection.

The inability of targeted BRAF inhibitors to produce long-lasting improvement in the clinical outcome of melanoma highlights a need to identify additional approaches to inhibit melanoma growth. Recent studies have shown that activation of the Wnt/β-catenin pathway decreases tumor growth and cooperates with ERK/MAPK pathway inhibitors to promote apoptosis in melanoma. Therefore, the identification of Wnt/β-catenin regulators may advance the development of new approaches to treat this disease. In order to move towards this goal we performed a large scale small-interfering RNA (siRNA) screen for regulators of β-catenin activated reporter activity in human HT1080 fibrosarcoma cells. Integrating large scale siRNA screen data with phosphoproteomic data and bioinformatics enrichment identified a protein, FAM129B, as a potential regulator of Wnt/β-catenin signaling. Functionally, we demonstrated that siRNA-mediated knockdown of FAM129B in A375 and A2058 melanoma cell lines inhibits WNT3A-mediated activation of a β-catenin-responsive luciferase reporter and inhibits expression of the endogenous Wnt/β-catenin target gene, AXIN2. We also demonstrate that FAM129B knockdown inhibits apoptosis in melanoma cells treated with WNT3A. These experiments support a role for FAM129B in linking Wnt/β-catenin signaling to apoptosis in melanoma.

Inhibitors of the checkpoint kinases CHK1 and WEE1 are under development as sensitizers to chemotherapy or radiation therapy. Disruption of either CHK1 or WEE1 activity compromises the cell cycle arrest that normally accompanies the DNA damage response, forcing cells to prematurely divide before repairing damaged DNA. However, inhibition of either CHK1 or WEE1 alone also disrupts DNA replication, resulting in DNA double strand breaks (DSBs) primarily in S-phase cells. Here we describe the ability of inhibitors against CHK1 and WEE1 to synergize both in vitro and in vivo. MK-8776 and MK-1775 are advancing through early clinical trials as potent and selective inhibitors of CHK1 and WEE1, respectively. MK-8776 increases the cell-based potency of MK-1775 in multiple cancer cell lines by as much as a 10-fold EC50 shift in proliferation assays. Staining for the DSB marker γH2AX at low concentrations of each drug demonstrates that MK-8776 and MK-1775 in combination induce up to 50-fold more DNA damage than either MK-8776 or MK-1775 alone. The ability of the combination to drive premature mitosis, determined by pHH3 quantitation, was variable, suggesting that generation of DSBs rather than G2/M checkpoint abrogation is the mechanism underlying combination cytotoxicity. Induction of DSBs requires cyclin-dependent kinase activity and takes place largely in S-phase cells, implying deregulation of DNA replication. Furthermore, staining for DNA damage response markers γH2AX and pCHK1-ser345 was both more intense and more durable in a LoVo xenograft model of cancer when MK-8776 and MK-1775 were dosed together compared to either compound alone. At tolerated doses, the combination demonstrated significantly greater anti-tumor efficacy than either MK-8776 or MK-1775 alone in LoVo xenografts (91% tumor growth inhibition for the combination versus 28% and 41% for MK-8776 and MK-1775, respectively). Similar observations were made in additional xenograft models. These results demonstrate in vitro and in vivo synergy between MK-8776 and MK-1775, indicating unique and complimentary anti-tumor effects of CHK1 and WEE1 inhibition. Our studies provide preclinical rationale for combining WEE1 and CHK1 inhibitors as a novel cancer treatment regimen. Citation 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 2051. doi:1538-7445.AM2012-2051

There has been increasing interest in developing cancer therapies targeting PI3K pathway nodes. However, inhibition of a single node in PI3K pathway such as mTOR by rapamycin analogues results in compensatory activation of survival signaling pathway such as AKT and thereby limiting monotherapy activity. Therefore, in most tumor types, multipathway inhibition, guided by an in-depth understanding of feedback and crosstalk between pathways, may be required for tumor regression. It is well known that that mTOR inhibition upregulates pAKT via the S6K-IRS2 negative feedback loop. Inhibiting both mTOR and IGF1R may therefore ablate such feedback upregulation and lead to clinical response. However, it is unclear which tumor types will respond to the combination in the clinic. Here, we present preclinical and molecular profiling data that supports LumB breast cancer as a potential indication for MK-8669/MK-0646 (mTOR/IGF1 R inhibitor combination). When a panel of over 60 breast cancer cell lines was treated with MK-8669, ER+ cell lines and HER2+ cell lines were clearly more responsive than the triple negative ones. Meanwhile, the status of the MAPK pathway, as captured by a RAS gene expression signature, correlated with resistance to MK-8669 across the panel. Moreover, MK-8669 treatment unregulated multiple key nodes of the IGF1R pathway in the breast cell line panel, especially IRS2, whose upregulation correlated with the response to MK-8669. Furthermore, IGF1R, whose mRNA level correlates with the response to MK-0646 (IGF1 Ri) in multiple internal and external studies, is the highest in a subset of LumB tumors, suggesting a potential dependency on the IGF1R pathway in those tumors. Finally, a comparison between the gene expression profiles of breast cancer cell lines and tumors revealed that the ER+ breast cancer cell lines only represented the LumB, not LumA, subtype of human breast tumors. Taken together, these findings allow us to hypothesize that LumB breast cancer may be enriched in responders to the MK-0646/MK-8669 combination. In 2009, Merck initiated a phase I trial of MK-8669/MK-0646 in which clinical response was observed in LumB breast cancer patients. A phase II clinical trial has been initiated to test the hypothesis retrospectively. This talk is also presented as Poster A52. Citation Information: Clin Cancer Res 2010;16(14 Suppl):PR3.

We have tested PI3K and MAPK targeting compounds in a panel of 69 breast cancer cell lines. Following the molecular stratification of breast cancer tumors, we have used ER and HER2 biomarkers to stratify cell lines into HER2 positive (HER2+), ER positive (ER+), and triple negative (TN; ER-, PR- and HER2-) subgroups. Both gene expression and response to compounds showed significant differences among the subtypes. Response to mTOR and AKT inhibition was the greatest in ER+ and HER2+ groups (86%, 85% for mTORi and 64%, 77% for AKTi response rates, respectively), whereas TN cell lines showed preferential sensitivity to MEK inhibition (30%). The combination of AKT plus MEK and mTOR plus MEK inhibitors showed synergy in TN tumors (56% and 80% response rates, respectively), particularly in the epithelial subgroup. We believe that RAS pathway deregulation, as reported by a gene expression signature of RAS activation, is the key latent variable that explains the difference in responses. This conclusion is based on pathway activation differences between subtypes. We observe that the TN group has both elevated PI3K and MAPK signaling, evident by high levels of RAS, pERK, and pMEK signatures, whereas the ER+/GGI+ group only has high levels of PI3K signaling. Thus, mTOR, AKT, or MEK alone is insufficient to inhibit the TN lines, whereas the combos are very effective. Finally, we have assessed the prevalence of the EMT and RAS signatures across subtypes of breast tumors. Similar to cell lines, the RAS signature is markedly higher in TN compared to ER+/GGI+ and, to a lesser extent, compared to HER2+ tumors. The EMT signature on the other hand is not significantly different among subtypes. This would suggest that the response patterns observed in ER+/GGI+ and TN/E subtypes of this breast panel will be applicable to a substantial subpopulation of breast tumors (∼38%). Importantly, these breast subtypes represent the largest unmet need and are a high priority for the clinical oncology. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 1986.

Read-through fusion transcripts that result from the splicing of two adjacent genes in the same coding orientation are a recently discovered type of chimeric RNA. We sought to determine if read-through fusion transcripts exist in breast cancer. We performed paired-end RNA-seq of 168 breast samples, including 28 breast cancer cell lines, 42 triple negative breast cancer primary tumors, 42 estrogen receptor positive (ER+) breast cancer primary tumors, and 56 non-malignant breast tissue samples. We analyzed the sequencing data to identify breast cancer associated read-through fusion transcripts. We discovered two recurrent read-through fusion transcripts that were identified in breast cancer cell lines, confirmed across breast cancer primary tumors, and were not detected in normal tissues (SCNN1A-TNFRSF1A and CTSD-IFITM10). Both fusion transcripts use canonical splice sites to join the last splice donor of the 5′ gene to the first splice acceptor of the 3′ gene, creating an in-frame fusion transcript. Western blots indicated that the fusion transcripts are translated into fusion proteins in breast cancer cells. Custom small interfering RNAs targeting the CTSD-IFITM10 fusion junction reduced expression of the fusion transcript and reduced breast cancer cell proliferation. Read-through fusion transcripts between adjacent genes with different biochemical functions represent a new type of recurrent molecular defect in breast cancer that warrant further investigation as potential biomarkers and therapeutic targets. Both breast cancer associated fusion transcripts identified in this study involve membrane proteins (SCNN1A-TNFRSF1A and CTSD-IFITM10), which raises the possibility that they could be breast cancer-specific cell surface markers. Electronic supplementary material The online version of this article (doi:10.1007/s10549-014-3019-2) contains supplementary material, which is available to authorized users.