ABSTRACT
Fibroblastic reticular cells (FRCs) are specialized stromal cells that define tissue architecture and regulate lymphocyte compartmentalization, homeostasis, and innate and adaptive immunity in secondary lymphoid organs (SLOs). In the present study, we used single-cell RNA sequencing (scRNA-seq) of human and mouse lymph nodes (LNs) to identify a subset of T cell-zone FRCs defined by the expression of Gremlin1 (Grem1) in both species. Grem1-CreERT2 knock-in mice enabled localization, multi-omics characterization and genetic depletion of Grem1+ FRCs. Grem1+ FRCs primarily localize at T-B cell junctions of SLOs, neighboring pre-dendritic cells and conventional dendritic cells (cDCs). As such, their depletion resulted in preferential loss and decreased homeostatic proliferation and survival of resident cDCs and compromised T cell immunity. Trajectory analysis of human LN scRNA-seq data revealed expression similarities to murine FRCs, with GREM1+ cells marking the endpoint of both trajectories. These findings illuminate a new Grem1+ fibroblastic niche in LNs that functions to maintain the homeostasis of lymphoid tissue-resident cDCs.
Subject(s)
Dendritic Cells, Follicular/immunology , Fibroblasts/immunology , Lymph Nodes/immunology , Stromal Cells/immunology , Aged , Animals , Apoptosis/genetics , Apoptosis/immunology , Cell Proliferation/genetics , Cell Survival/genetics , Cell Survival/immunology , Dendritic Cells, Follicular/metabolism , Female , Fibroblasts/metabolism , Gene Expression Regulation/immunology , Gene Knock-In Techniques , Humans , Immunity, Cellular/genetics , Intercellular Signaling Peptides and Proteins/genetics , Intercellular Signaling Peptides and Proteins/metabolism , Lymph Nodes/cytology , Male , Mice , Mice, Transgenic , RNA-Seq , Single-Cell Analysis , Stromal Cells/metabolism , T-Lymphocytes/immunology , T-Lymphocytes/metabolismABSTRACT
Fibroblasts are non-haematopoietic structural cells that define the architecture of organs, support the homeostasis of tissue-resident cells and have key roles in fibrosis, cancer, autoimmunity and wound healing1. Recent studies have described fibroblast heterogeneity within individual tissues1. However, the field lacks a characterization of fibroblasts at single-cell resolution across tissues in healthy and diseased organs. Here we constructed fibroblast atlases by integrating single-cell transcriptomic data from about 230,000 fibroblasts across 17 tissues, 50 datasets, 11 disease states and 2 species. Mouse fibroblast atlases and a DptIRESCreERT2 knock-in mouse identified two universal fibroblast transcriptional subtypes across tissues. Our analysis suggests that these cells can serve as a reservoir thatĀ can yield specialized fibroblasts across a broad range of steady-state tissues and activated fibroblasts in disease. Comparison to an atlas of human fibroblasts from perturbed states showed that fibroblast transcriptional states are conserved between mice and humans, including universal fibroblasts and activated phenotypes associated with pathogenicity in human cancer, fibrosis, arthritis and inflammation. In summary, a cross-species and pan-tissue approach to transcriptomics at single-cell resolution has identified key organizing principles of the fibroblast lineage in health and disease.
Subject(s)
Fibroblasts/cytology , Transcriptome , Animals , Cells, Cultured , Disease , Female , Fibroblasts/classification , Gene Knock-In Techniques , Humans , Male , Mice , Mice, Inbred C57BL , Neoplasms , Organ Specificity , Phenotype , RNA-Seq , Single-Cell Analysis , Stromal CellsABSTRACT
Despite the resounding clinical success in cancer treatment of antibodies that block the interaction of PD1 with its ligand PDL11, the mechanisms involved remain unknown. A major limitation to understanding the origin and fate of T cells in tumour immunity is the lack of quantitative information on the distribution of individual clonotypes of T cells in patients with cancer. Here, by performing deep single-cell sequencing of RNA and T cell receptors in patients with different types of cancer, we survey the profiles of various populations of T cells and T cell receptors in tumours, normal adjacent tissue, and peripheral blood. We find clear evidence of clonotypic expansion of effector-like T cells not only within the tumour but also in normal adjacent tissue. Patients with gene signatures of such clonotypic expansion respond best to anti-PDL1 therapy. Notably, expanded clonotypes found in the tumour and normal adjacent tissue can also typically be detected in peripheral blood, which suggests a convenient approach to patient identification. Analyses of our data together with several external datasets suggest that intratumoural T cells, especially in responsive patients, are replenished with fresh, non-exhausted replacement cells from sites outside the tumour, suggesting continued activity of the cancer immunity cycle in these patients, the acceleration of which may be associated with clinical response.
Subject(s)
Lymphocytes, Tumor-Infiltrating/cytology , Lymphocytes, Tumor-Infiltrating/metabolism , Neoplasms/pathology , Pharmacogenomic Variants , Receptors, Antigen, T-Cell/genetics , T-Lymphocytes/cytology , Antibodies, Monoclonal, Humanized/therapeutic use , Antineoplastic Agents/therapeutic use , Clone Cells , Humans , Neoplasms/drug therapy , Neoplasms/immunology , T-Lymphocytes/metabolism , TranscriptomeABSTRACT
Therapeutic antibodies that block the programmed death-1 (PD-1)-programmed death-ligand 1 (PD-L1) pathway can induce robust and durable responses in patients with various cancers, including metastatic urothelial cancer. However, these responses only occur in a subset of patients. Elucidating the determinants of response and resistance is key to improving outcomes and developing new treatment strategies. Here we examined tumours from a large cohort of patients with metastatic urothelial cancer who were treated with an anti-PD-L1 agent (atezolizumab) and identified major determinants of clinical outcome. Response to treatment was associated with CD8+ T-effector cell phenotype and, to an even greater extent, high neoantigen or tumour mutation burden. Lack of response was associated with a signature of transforming growth factor Ć (TGFĆ) signalling in fibroblasts. This occurred particularly in patients with tumours, which showed exclusion of CD8+ T cells from the tumour parenchyma that were instead found in the fibroblast- and collagen-rich peritumoural stroma; a common phenotype among patients with metastatic urothelial cancer. Using a mouse model that recapitulates this immune-excluded phenotype, we found that therapeutic co-administration of TGFĆ-blocking and anti-PD-L1 antibodies reduced TGFĆ signalling in stromal cells, facilitated T-cell penetration into the centre of tumours, and provoked vigorous anti-tumour immunity and tumour regression. Integration of these three independent biological features provides the best basis for understanding patient outcome in this setting and suggests that TGFĆ shapes the tumour microenvironment to restrain anti-tumour immunity by restricting T-cell infiltration.
Subject(s)
Antibodies, Monoclonal/pharmacology , Antibodies, Monoclonal/therapeutic use , B7-H1 Antigen/antagonists & inhibitors , CD8-Positive T-Lymphocytes/drug effects , Transforming Growth Factor beta/metabolism , Urologic Neoplasms/drug therapy , Urologic Neoplasms/immunology , Urothelium/pathology , Animals , Antibodies/immunology , Antibodies/pharmacology , Antibodies/therapeutic use , Antibodies, Monoclonal, Humanized , Antigens, Neoplasm/analysis , Antigens, Neoplasm/immunology , Antigens, Neoplasm/metabolism , B7-H1 Antigen/immunology , CD8-Positive T-Lymphocytes/cytology , CD8-Positive T-Lymphocytes/immunology , Cell Cycle Checkpoints/drug effects , Cohort Studies , Collagen/metabolism , Disease Models, Animal , Drug Resistance, Neoplasm/drug effects , Fibroblasts/metabolism , Humans , Immunotherapy , Mice , Mutation , Neoplasm Metastasis , Phenotype , Signal Transduction/drug effects , Transforming Growth Factor beta/antagonists & inhibitors , Treatment Outcome , Tumor Microenvironment/immunology , Urologic Neoplasms/genetics , Urologic Neoplasms/pathology , Urothelium/drug effects , Urothelium/immunologyABSTRACT
Many cancers evade immune rejection by suppressing major histocompatibility class I (MHC-I) antigen processing and presentation (AgPP). Such cancers do not respond to immune checkpoint inhibitor therapies (ICIT) such as PD-1/PD-L1 [PD-(L)1] blockade. Certain chemotherapeutic drugs augment tumor control by PD-(L)1 inhibitors through potentiation of T-cell priming but whether and how chemotherapy enhances MHC-I-dependent cancer cell recognition by cytotoxic T cells (CTLs) is not entirely clear. We now show that the lysine acetyl transferases p300/CREB binding protein (CBP) control MHC-I AgPPM expression and neoantigen amounts in human cancers. Moreover, we found that two distinct DNA damaging drugs, the platinoid oxaliplatin and the topoisomerase inhibitor mitoxantrone, strongly up-regulate MHC-I AgPP in a manner dependent on activation of nuclear factor kappa B (NF-κB), p300/CBP, and other transcription factors, but independently of autocrine IFNĆĀ³ signaling. Accordingly, NF-κB and p300 ablations prevent chemotherapy-induced MHC-I AgPP and abrogate rejection of low MHC-I-expressing tumors by reinvigorated CD8+ CTLs. Drugs like oxaliplatin and mitoxantrone may be used to overcome resistance to PD-(L)1 inhibitors in tumors that had "epigenetically down-regulated," but had not permanently lost MHC-I AgPP activity.
Subject(s)
Antigen Presentation/immunology , Gene Expression Regulation, Neoplastic/drug effects , Histocompatibility Antigens Class I/immunology , Immune Checkpoint Inhibitors/pharmacology , NF-kappa B/metabolism , Neoplasms/drug therapy , p300-CBP Transcription Factors/metabolism , Animals , Antineoplastic Agents/pharmacology , Apoptosis , B7-H1 Antigen/genetics , B7-H1 Antigen/metabolism , Biomarkers, Tumor/genetics , Biomarkers, Tumor/metabolism , CD8-Positive T-Lymphocytes , Cell Proliferation , Drug Therapy, Combination , Humans , Immunotherapy/methods , Mice , NF-kappa B/genetics , Neoplasms/immunology , Neoplasms/metabolism , Neoplasms/pathology , Oxaliplatin/pharmacology , Prognosis , Survival Rate , Tumor Cells, Cultured , Xenograft Model Antitumor Assays , p300-CBP Transcription Factors/geneticsABSTRACT
The use of large-scale genomic and drug response screening of cancer cell lines depends crucially on the reproducibility of results. Here we consider two previously published screens, plus a later critique of these studies. Using independent data, we show that consistency is achievable, and provide a systematic description of the best laboratory and analysis practices for future studies.
Subject(s)
Antineoplastic Agents/pharmacology , Drug Screening Assays, Antitumor/methods , Drug Screening Assays, Antitumor/standards , Neoplasms/genetics , Neoplasms/pathology , Cell Line, Tumor , Cell Survival/drug effects , Genetic Markers/genetics , Genome, Human/genetics , Humans , Quality Control , Reproducibility of ResultsABSTRACT
Cell line misidentification, contamination and poor annotation affect scientific reproducibility. Here we outline simple measures to detect or avoid cross-contamination, present a framework for cell line annotation linked to short tandem repeat and single nucleotide polymorphism profiles, and provide a catalogue of synonymous cell lines. This resource will enable our community to eradicate the use of misidentified lines and generate credible cell-based data.
Subject(s)
Cell Line/classification , Cell Line/metabolism , Data Curation , Guidelines as Topic , Cell Separation , Genotype , Microsatellite Repeats/genetics , Polymorphism, Single Nucleotide/genetics , Quality Control , Reproducibility of Results , Species Specificity , Terminology as TopicABSTRACT
BACKGROUND: First-line chemotherapy for patients with cisplatin-ineligible locally advanced or metastatic urothelial carcinoma is associated with short response duration, poor survival, and high toxicity. This study assessed atezolizumab (anti-programmed death-ligand 1 [PD-L1]) as treatment for metastatic urothelial cancer in cisplatin-ineligible patients. METHODS: For this single-arm, multicentre, phase 2 study, in 47 academic medical centres and community oncology practices in seven countries in North America and Europe, we recruited previously untreated patients with locally advanced or metastatic urothelial cancer who were cisplatin ineligible. Patients were given 1200 mg intravenous atezolizumab every 21 days until progression. The primary endpoint was independently confirmed objective response rate per Response Evaluation Criteria in Solid Tumors version 1.1 (central review), assessed in prespecified subgroups based on PD-L1 expression and in all patients. All participants who received one or more doses of atezolizumab were included in the primary and safety analyses. This study was registered with ClinicalTrials.gov, number NCT02108652. FINDINGS: Between June 9, 2014, and March 30, 2015, we enrolled 123 patients, of whom 119 received one or more doses of atezolizumab. At 17Ā·2 months' median follow-up, the objective response rate was 23% (95% CI 16 to 31), the complete response rate was 9% (n=11), and 19 of 27 responses were ongoing. Median response duration was not reached. Responses occurred across all PD-L1 and poor prognostic factor subgroups. Median progression-free survival was 2Ā·7 months (2Ā·1 to 4Ā·2). Median overall survival was 15Ā·9 months (10Ā·4 to not estimable). Tumour mutation load was associated with response. Treatment-related adverse events that occurred in 10% or more of patients were fatigue (36 [30%] patients), diarrhoea (14 [12%] patients), and pruritus (13 [11%] patients). One treatment-related death (sepsis) occurred. Nine (8%) patients had an adverse event leading to treatment discontinuation. Immune-mediated events occurred in 14 (12%) patients. INTERPRETATION: Atezolizumab showed encouraging durable response rates, survival, and tolerability, supporting its therapeutic use in untreated metastatic urothelial cancer. FUNDING: F Hoffmann-La Roche, Genentech.
Subject(s)
Antibodies, Monoclonal/therapeutic use , Antineoplastic Agents/therapeutic use , Carcinoma, Transitional Cell/drug therapy , Carcinoma, Transitional Cell/secondary , Urologic Neoplasms/drug therapy , Aged , Aged, 80 and over , Antibodies, Monoclonal/administration & dosage , Antibodies, Monoclonal/adverse effects , Antibodies, Monoclonal, Humanized , Antineoplastic Agents/administration & dosage , Antineoplastic Agents/adverse effects , B7-H1 Antigen/blood , Biomarkers, Tumor/blood , Carcinoma, Transitional Cell/blood , Cisplatin , Contraindications , Female , Humans , Infusions, Intravenous , Kaplan-Meier Estimate , Lymphatic Metastasis , Male , Middle Aged , Response Evaluation Criteria in Solid Tumors , Urologic Neoplasms/bloodABSTRACT
Identifying and understanding changes in cancer genomes is essential for the development of targeted therapeutics. Here we analyse systematically more than 70 pairs of primary human colon tumours by applying next-generation sequencing to characterize their exomes, transcriptomes and copy-number alterations. We have identified 36,303 protein-altering somatic changes that include several new recurrent mutations in the Wnt pathway gene TCF7L2, chromatin-remodelling genes such as TET2 and TET3 and receptor tyrosine kinases including ERBB3. Our analysis for significantly mutated cancer genes identified 23 candidates, including the cell cycle checkpoint kinase ATM. Copy-number and RNA-seq data analysis identified amplifications and corresponding overexpression of IGF2 in a subset of colon tumours. Furthermore, using RNA-seq data we identified multiple fusion transcripts including recurrent gene fusions involving R-spondin family members RSPO2 and RSPO3 that together occur in 10% of colon tumours. The RSPO fusions were mutually exclusive with APC mutations, indicating that they probably have a role in the activation of Wnt signalling and tumorigenesis. Consistent with this we show that the RSPO fusion proteins were capable of potentiating Wnt signalling. The R-spondin gene fusions and several other gene mutations identified in this study provide new potential opportunities for therapeutic intervention in colon cancer.
Subject(s)
Colonic Neoplasms/genetics , Gene Fusion/genetics , Genes, Neoplasm/genetics , Intercellular Signaling Peptides and Proteins/genetics , Thrombospondins/genetics , Ataxia Telangiectasia Mutated Proteins , Base Sequence , Cell Cycle Proteins/genetics , Colonic Neoplasms/metabolism , Colonic Neoplasms/pathology , DNA Copy Number Variations/genetics , DNA-Binding Proteins/genetics , Dioxygenases/genetics , Exome/genetics , Gene Expression Profiling , Gene Expression Regulation, Neoplastic/genetics , Genes, APC , Humans , Insulin-Like Growth Factor II/genetics , Molecular Sequence Data , Mutation/genetics , Polymorphism, Single Nucleotide/genetics , Protein Serine-Threonine Kinases/genetics , Proto-Oncogene Proteins/genetics , Receptor, ErbB-3/genetics , Sequence Analysis, RNA , Signal Transduction/genetics , Transcription Factor 7-Like 2 Protein/genetics , Tumor Suppressor Proteins/genetics , Wnt Proteins/metabolismABSTRACT
TNF superfamily death ligands are expressed on the surface of immune cells and can trigger apoptosis in susceptible cancer cells by engaging cognate death receptors. A recombinant soluble protein comprising the ectodomain of Apo2 ligand/TNF-related apoptosis-inducing ligand (Apo2L/TRAIL) has shown remarkable preclinical anticancer activity but lacked broad efficacy in patients, possibly owing to insufficient exposure or potency. We observed that antibody cross-linking substantially enhanced cytotoxicity of soluble Apo2L/TRAIL against diverse cancer cell lines. Presentation of the ligand on glass-supported lipid bilayers enhanced its ability to drive receptor microclustering and apoptotic signaling. Furthermore, covalent surface attachment of Apo2L/TRAIL onto liposomes--synthetic lipid-bilayer nanospheres--similarly augmented activity. In vivo, liposome-displayed Apo2L/TRAIL achieved markedly better exposure and antitumor activity. Thus, covalent synthetic-membrane attachment of a cell-surface ligand enhances efficacy, increasing therapeutic potential. These findings have translational implications for liposomal approaches as well as for Apo2L/TRAIL and other clinically relevant TNF ligands.
Subject(s)
Antineoplastic Agents/chemistry , Cell Membrane/metabolism , TNF-Related Apoptosis-Inducing Ligand/metabolism , Animals , Apoptosis , Biotinylation , CD27 Ligand/metabolism , Caspase 8/metabolism , Caspases/metabolism , Cell Line, Tumor , Epitopes/chemistry , Fas Ligand Protein/metabolism , Humans , Immune System , Immunotherapy/methods , Inhibitory Concentration 50 , Ligands , Liposomes/chemistry , Mice , Mice, Nude , Microscopy, Fluorescence , Neoplasm Transplantation , Neoplasms/immunology , Neoplasms/metabolism , Recombinant Proteins/metabolismABSTRACT
BACKGROUND: Patients with metastatic urothelial carcinoma have few treatment options after failure of platinum-based chemotherapy. In this trial, we assessed treatment with atezolizumab, an engineered humanised immunoglobulin G1 monoclonal antibody that binds selectively to programmed death ligand 1 (PD-L1), in this patient population. METHODS: For this multicentre, single-arm, two-cohort, phase 2 trial, patients (aged ≥18 years) with inoperable locally advanced or metastatic urothelial carcinoma whose disease had progressed after previous platinum-based chemotherapy were enrolled from 70 major academic medical centres and community oncology practices in Europe and North America. Key inclusion criteria for enrolment were Eastern Cooperative Oncology Group performance status of 0 or 1, measurable disease defined by Response Evaluation Criteria In Solid Tumors version 1.1 (RECIST v1.1), adequate haematological and end-organ function, and no autoimmune disease or active infections. Formalin-fixed paraffin-embedded tumour specimens with sufficient viable tumour content were needed from all patients before enrolment. Patients received treatment with intravenous atezolizumab (1200 mg, given every 3 weeks). PD-L1 expression on tumour-infiltrating immune cells (ICs) was assessed prospectively by immunohistochemistry. The co-primary endpoints were the independent review facility-assessed objective response rate according to RECIST v1.1 and the investigator-assessed objective response rate according to immune-modified RECIST, analysed by intention to treat. A hierarchical testing procedure was used to assess whether the objective response rate was significantly higher than the historical control rate of 10% at an α level of 0Ā·05. This study is registered with ClinicalTrials.gov, number NCT02108652. FINDINGS: Between May 13, 2014, and Nov 19, 2014, 486 patients were screened and 315 patients were enrolled into the study. Of these patients, 310 received atezolizumab treatment (five enrolled patients later did not meet eligibility criteria and were not dosed with study drug). The PD-L1 expression status on infiltrating immune cells (ICs) in the tumour microenvironment was defined by the percentage of PD-L1-positive immune cells: IC0 (<1%), IC1 (≥1% but <5%), and IC2/3 (≥5%). The primary analysis (data cutoff May 5, 2015) showed that compared with a historical control overall response rate of 10%, treatment with atezolizumab resulted in a significantly improved RECIST v1.1 objective response rate for each prespecified immune cell group (IC2/3: 27% [95% CI 19-37], p<0Ā·0001; IC1/2/3: 18% [13-24], p=0Ā·0004) and in all patients (15% [11-20], p=0Ā·0058). With longer follow-up (data cutoff Sept 14, 2015), by independent review, objective response rates were 26% (95% CI 18-36) in the IC2/3 group, 18% (13-24) in the IC1/2/3 group, and 15% (11-19) overall in all 310 patients. With a median follow-up of 11Ā·7 months (95% CI 11Ā·4-12Ā·2), ongoing responses were recorded in 38 (84%) of 45 responders. Exploratory analyses showed The Cancer Genome Atlas (TCGA) subtypes and mutation load to be independently predictive for response to atezolizumab. Grade 3-4 treatment-related adverse events, of which fatigue was the most common (five patients [2%]), occurred in 50 (16%) of 310 treated patients. Grade 3-4 immune-mediated adverse events occurred in 15 (5%) of 310 treated patients, with pneumonitis, increased aspartate aminotransferase, increased alanine aminotransferase, rash, and dyspnoea being the most common. No treatment-related deaths occurred during the study. INTERPRETATION: Atezolizumab showed durable activity and good tolerability in this patient population. Increased levels of PD-L1 expression on immune cells were associated with increased response. This report is the first to show the association of TCGA subtypes with response to immune checkpoint inhibition and to show the importance of mutation load as a biomarker of response to this class of agents in advanced urothelial carcinoma. FUNDING: F Hoffmann-La Roche Ltd.
Subject(s)
Antibodies, Monoclonal/administration & dosage , Antineoplastic Agents/administration & dosage , B7-H1 Antigen/antagonists & inhibitors , Urologic Neoplasms/drug therapy , Adult , Aged , Aged, 80 and over , Antibodies, Monoclonal, Humanized/administration & dosage , Antineoplastic Combined Chemotherapy Protocols/therapeutic use , B7-H1 Antigen/immunology , Carboplatin/administration & dosage , Cisplatin/administration & dosage , Disease Progression , Drug Administration Schedule , Female , Humans , Immunoglobulins, Intravenous/administration & dosage , Infusions, Intravenous , Male , Middle Aged , Mutation/genetics , Neoplasm Metastasis , Response Evaluation Criteria in Solid Tumors , Treatment Outcome , Urologic Neoplasms/geneticsABSTRACT
The systematic characterization of somatic mutations in cancer genomes is essential for understanding the disease and for developing targeted therapeutics. Here we report the identification of 2,576 somatic mutations across approximately 1,800 megabases of DNA representing 1,507 coding genes from 441 tumours comprising breast, lung, ovarian and prostate cancer types and subtypes. We found that mutation rates and the sets of mutated genes varied substantially across tumour types and subtypes. Statistical analysis identified 77 significantly mutated genes including protein kinases, G-protein-coupled receptors such as GRM8, BAI3, AGTRL1 (also called APLNR) and LPHN3, and other druggable targets. Integrated analysis of somatic mutations and copy number alterations identified another 35 significantly altered genes including GNAS, indicating an expanded role for galpha subunits in multiple cancer types. Furthermore, our experimental analyses demonstrate the functional roles of mutant GNAO1 (a Galpha subunit) and mutant MAP2K4 (a member of the JNK signalling pathway) in oncogenesis. Our study provides an overview of the mutational spectra across major human cancers and identifies several potential therapeutic targets.
Subject(s)
Genes, Neoplasm/genetics , Mutation/genetics , Neoplasms/genetics , Neoplasms/metabolism , Signal Transduction/genetics , Breast Neoplasms/classification , Breast Neoplasms/genetics , DNA Copy Number Variations/genetics , DNA Mutational Analysis , Female , GTP-Binding Protein alpha Subunits/genetics , Humans , Lung Neoplasms/classification , Lung Neoplasms/genetics , MAP Kinase Kinase 4/genetics , Male , Neoplasms/enzymology , Neoplasms/pathology , Ovarian Neoplasms/classification , Ovarian Neoplasms/genetics , Prostatic Neoplasms/classification , Prostatic Neoplasms/genetics , Protein Kinases/genetics , Receptors, G-Protein-Coupled/geneticsABSTRACT
UNLABELLED: Connections between disease phenotypes and drug effects can be made by identifying commonalities in the associated patterns of differential gene expression. Searchable databases that record the impacts of chemical or genetic perturbations on the transcriptome--here referred to as 'connectivity maps'--permit discovery of such commonalities. We describe two R packages, gCMAP and gCMAPWeb, which provide a complete framework to construct and query connectivity maps assembled from user-defined collections of differential gene expression data. Microarray or RNAseq data are processed in a standardized way, and results can be interrogated using various well-established gene set enrichment methods. The packages also feature an easy-to-deploy web application that facilitates reproducible research through automatic generation of graphical and tabular reports. AVAILABILITY AND IMPLEMENTATION: The gCMAP and gCMAPWeb R packages are freely available for UNIX, Windows and Mac OS X operating systems at Bioconductor (http://www.bioconductor.org).
Subject(s)
Oligonucleotide Array Sequence Analysis/methods , User-Computer Interface , Animals , Cell Line , Gene Expression Profiling/methods , Humans , InternetABSTRACT
Polycomb group (PcG) complexes are multiprotein assemblages that bind to chromatin and establish chromatin states leading to epigenetic silencing. PcG proteins regulate homeotic genes in flies and vertebrates, but little is known about other PcG targets and the role of the PcG in development, differentiation and disease. Here, we determined the distribution of the PcG proteins PC, E(Z) and PSC and of trimethylation of histone H3 Lys27 (me3K27) in the D. melanogaster genome. At more than 200 PcG target genes, binding sites for the three PcG proteins colocalize to presumptive Polycomb response elements (PREs). In contrast, H3 me3K27 forms broad domains including the entire transcription unit and regulatory regions. PcG targets are highly enriched in genes encoding transcription factors, but they also include genes coding for receptors, signaling proteins, morphogens and regulators representing all major developmental pathways.
Subject(s)
Drosophila Proteins/genetics , Drosophila melanogaster/genetics , Genome , Animals , Cell Line , Chromatin Immunoprecipitation , DNA Methylation , Polycomb Repressive Complex 1 , Reverse Transcriptase Polymerase Chain ReactionABSTRACT
Meiotic recombination has a central role in the evolution of sexually reproducing organisms. The two recombination outcomes, crossover and non-crossover, increase genetic diversity, but have the potential to homogenize alleles by gene conversion. Whereas crossover rates vary considerably across the genome, non-crossovers and gene conversions have only been identified in a handful of loci. To examine recombination genome wide and at high spatial resolution, we generated maps of crossovers, crossover-associated gene conversion and non-crossover gene conversion using dense genetic marker data collected from all four products of fifty-six yeast (Saccharomyces cerevisiae) meioses. Our maps reveal differences in the distributions of crossovers and non-crossovers, showing more regions where either crossovers or non-crossovers are favoured than expected by chance. Furthermore, we detect evidence for interference between crossovers and non-crossovers, a phenomenon previously only known to occur between crossovers. Up to 1% of the genome of each meiotic product is subject to gene conversion in a single meiosis, with detectable bias towards GC nucleotides. To our knowledge the maps represent the first high-resolution, genome-wide characterization of the multiple outcomes of recombination in any organism. In addition, because non-crossover hotspots create holes of reduced linkage within haplotype blocks, our results stress the need to incorporate non-crossovers into genetic linkage analysis.
Subject(s)
Chromosome Mapping , Crossing Over, Genetic/genetics , Gene Conversion/genetics , Meiosis/genetics , Saccharomyces cerevisiae/genetics , Chromosomes, Fungal/genetics , DNA Breaks, Double-Stranded , DNA-Binding Proteins/genetics , Flap Endonucleases , Gene Expression Regulation, Fungal , Genetic Linkage/genetics , Genetic Markers/genetics , Genome, Fungal/genetics , Genotype , Haplotypes/genetics , Mutation/genetics , Saccharomyces cerevisiae/cytology , Saccharomyces cerevisiae Proteins/genetics , Trans-Activators/geneticsABSTRACT
The mutant forms of KRas, NRas and HRas drive the initiation and progression of a number of human cancers, but less is known about the role of WT (wild-type) Ras alleles and isoforms in cancer. We used zinc-finger nucleases targeting HRas and NRas to modify both alleles of these genes in the mutant KRas-driven Hec1A endometrial cancer cell line, which normally expresses WT copies of these genes. The disruption of either WT isoform of Ras compromised growth-factor-dependent signalling through the ERK (extracellular-signal-regulated kinase) pathway. In addition, the disruption of HRas hindered the activation of Akt and subsequent downstream signalling. This was associated with decreased proliferation, increased apoptosis and decreased anchorage-independent growth in the HRas-disrupted cells. However, xenograft tumour growth was not significantly affected by the disruption of either NRas or HRas. As expected, deleting the mutant allele of KRas abolished tumour growth, whereas deletion of the remaining WT copy of KRas increased the tumorigenic properties of these cells; deleting a single copy of either HRas or NRas did not mimic this effect. The present study demonstrates that the WT copies of HRas, NRas and KRas play unique roles in the context of mutant KRas-driven tumours.
Subject(s)
Cell Transformation, Neoplastic/genetics , Mutant Proteins/chemistry , Mutant Proteins/genetics , Proto-Oncogene Proteins/chemistry , Proto-Oncogene Proteins/genetics , Signal Transduction/genetics , ras Proteins/chemistry , ras Proteins/genetics , Amino Acid Sequence , Animals , Base Sequence , Cell Line, Tumor , Cell Transformation, Neoplastic/chemistry , Female , GTP Phosphohydrolases/chemistry , GTP Phosphohydrolases/genetics , Humans , Membrane Proteins/chemistry , Membrane Proteins/genetics , Mice , Mice, Inbred BALB C , Mice, Nude , Molecular Sequence Data , Protein Isoforms/chemistry , Protein Isoforms/genetics , Proto-Oncogene Proteins p21(ras)/chemistry , Proto-Oncogene Proteins p21(ras)/geneticsABSTRACT
The advent of high-throughput biotechnologies, which can efficiently measure gene expression on a global basis, has led to the creation and population of correspondingly rich databases and compendia. Such repositories have the potential to add enormous scientific value beyond that provided by individual studies which, due largely to cost considerations, are typified by small sample sizes. Accordingly, substantial effort has been invested in devising analysis schemes for utilizing gene-expression repositories. Here, we focus on one such scheme, the Connectivity Map (cmap), that was developed with the express purpose of identifying drugs with putative efficacy against a given disease, where the disease in question is characterized by a (differential) gene-expression signature. Initial claims surrounding cmap intimated that such tools might lead to new, previously unanticipated applications of existing drugs. However, further application suggests that its primary utility is in connecting a disease condition whose biology is largely unknown to a drug whose mechanisms of action are well understood, making cmap a tool for enhancing biological knowledge.The success of the Connectivity Map is belied by its simplicity. The aforementioned signature serves as an unordered query which is applied to a customized database of (differential) gene-expression experiments designed to elicit response to a wide range of drugs, across of spectrum of concentrations, durations, and cell lines. Such application is effected by computing a per experiment score that measures "closeness" between the signature and the experiment. Top-scoring experiments, and the attendant drug(s), are then deemed relevant to the disease underlying the query. Inference supporting such elicitations is pursued via re-sampling. In this paper, we revisit two key aspects of the Connectivity Map implementation. Firstly, we develop new approaches to measuring closeness for the common scenario wherein the query constitutes an ordered list. These involve using metrics proposed for analyzing partially ranked data, these being of interest in their own right and not widely used. Secondly, we advance an alternate inferential approach based on generating empirical null distributions that exploit the scope, and capture dependencies, embodied by the database. Using these refinements we undertake a comprehensive re-evaluation of Connectivity Map findings that, in general terms, reveal that accommodating ordered queries is less critical than the mode of inference.
Subject(s)
Data Mining/methods , Databases, Genetic , Gene Expression Profiling , Algorithms , Computational Biology/methods , Estrogens/pharmacology , Gene Expression/drug effects , Genetic Predisposition to Disease , Genomics/methods , Histone Deacetylase Inhibitors/pharmacology , Humans , Limonins/pharmacologyABSTRACT
With high-dimensional data, variable-by-variable statistical testing is often used to select variables whose behavior differs across conditions. Such an approach requires adjustment for multiple testing, which can result in low statistical power. A two-stage approach that first filters variables by a criterion independent of the test statistic, and then only tests variables which pass the filter, can provide higher power. We show that use of some filter/test statistics pairs presented in the literature may, however, lead to loss of type I error control. We describe other pairs which avoid this problem. In an application to microarray data, we found that gene-by-gene filtering by overall variance followed by a t-test increased the number of discoveries by 50%. We also show that this particular statistic pair induces a lower bound on fold-change among the set of discoveries. Independent filtering-using filter/test pairs that are independent under the null hypothesis but correlated under the alternative-is a general approach that can substantially increase the efficiency of experiments.
Subject(s)
Biometry/methods , Algorithms , Computational Biology , Models, GeneticABSTRACT
Polycomb (PcG) regulation has been thought to produce stable long-term gene silencing. Genomic analyses in Drosophila and mammals, however, have shown that it targets many genes, which can switch state during development. Genetic evidence indicates that critical for the active state of PcG target genes are the histone methyltransferases Trithorax (TRX) and ASH1. Here we analyze the repertoire of alternative states in which PcG target genes are found in different Drosophila cell lines and the role of PcG proteins TRX and ASH1 in controlling these states. Using extensive genome-wide chromatin immunoprecipitation analysis, RNAi knockdowns, and quantitative RT-PCR, we show that, in addition to the known repressed state, PcG targets can reside in a transcriptionally active state characterized by formation of an extended domain enriched in ASH1, the N-terminal, but not C-terminal moiety of TRX and H3K27ac. ASH1/TRX N-ter domains and transcription are not incompatible with repressive marks, sometimes resulting in a "balanced" state modulated by both repressors and activators. Often however, loss of PcG repression results instead in a "void" state, lacking transcription, H3K27ac, or binding of TRX or ASH1. We conclude that PcG repression is dynamic, not static, and that the propensity of a target gene to switch states depends on relative levels of PcG, TRX, and activators. N-ter TRX plays a remarkable role that antagonizes PcG repression and preempts H3K27 methylation by acetylation. This role is distinct from that usually attributed to TRX/MLL proteins at the promoter. These results have important implications for Polycomb gene regulation, the "bivalent" chromatin state of embryonic stem cells, and gene expression in development.
Subject(s)
Chromatin/metabolism , Drosophila Proteins/metabolism , Drosophila melanogaster/metabolism , Epigenesis, Genetic , Acetylation , Animals , Cell Line , Chromatin/genetics , DNA-Binding Proteins/genetics , DNA-Binding Proteins/metabolism , Drosophila Proteins/genetics , Drosophila melanogaster/genetics , Drosophila melanogaster/growth & development , Membrane Proteins/genetics , Membrane Proteins/metabolism , Polycomb Repressive Complex 1 , Protein Binding , Thioredoxins/genetics , Thioredoxins/metabolism , Transcription Factors/genetics , Transcription Factors/metabolismABSTRACT
Most colorectal (CRC) tumors are dependent on EGFR/KRAS/BRAF/MAPK signaling activation. ARID1A is an epigenetic regulator mutated in approximately 5% of non-hypermutated CRC tumors. Here we show that anti-EGFR but not anti-VEGF treatment enriches for emerging ARID1A mutations in CRC patients. In addition, we find that patients with ARID1A mutations, at baseline, are associated with worse outcome when treated with cetuximab- but not bevacizumab-containing therapies; thus, this suggests that ARID1A mutations may provide both an acquired and intrinsic mechanism of resistance to anti-EGFR therapies. We find that, ARID1A and EGFR-pathway genetic alterations are mutually exclusive across lung and colorectal cancers, further supporting a functional connection between these pathways. Our results not only suggest that ARID1A could be potentially used as a predictive biomarker for cetuximab treatment decisions but also provide a rationale for exploring therapeutic MAPK inhibition in an unexpected but genetically defined segment of CRC patients.