An Integrated TCGA Pan-Cancer Clinical Data Resource to Drive High-Quality Survival Outcome Analytics.

For a decade, The Cancer Genome Atlas (TCGA) program collected clinicopathologic annotation data along with multi-platform molecular profiles of more than 11,000 human tumors across 33 different cancer types. TCGA clinical data contain key features representing the democratized nature of the data collection process. To ensure proper use of this large clinical dataset associated with genomic features, we developed a standardized dataset named the TCGA Pan-Cancer Clinical Data Resource (TCGA-CDR), which includes four major clinical outcome endpoints. In addition to detailing major challenges and statistical limitations encountered during the effort of integrating the acquired clinical data, we present a summary that includes endpoint usage recommendations for each cancer type. These TCGA-CDR findings appear to be consistent with cancer genomics studies independent of the TCGA effort and provide opportunities for investigating cancer biology using clinical correlates at an unprecedented scale.

Cell-of-Origin Patterns Dominate the Molecular Classification of 10,000 Tumors from 33 Types of Cancer.

We conducted comprehensive integrative molecular analyses of the complete set of tumors in The Cancer Genome Atlas (TCGA), consisting of approximately 10,000 specimens and representing 33 types of cancer. We performed molecular clustering using data on chromosome-arm-level aneuploidy, DNA hypermethylation, mRNA, and miRNA expression levels and reverse-phase protein arrays, of which all, except for aneuploidy, revealed clustering primarily organized by histology, tissue type, or anatomic origin. The influence of cell type was evident in DNA-methylation-based clustering, even after excluding sites with known preexisting tissue-type-specific methylation. Integrative clustering further emphasized the dominant role of cell-of-origin patterns. Molecular similarities among histologically or anatomically related cancer types provide a basis for focused pan-cancer analyses, such as pan-gastrointestinal, pan-gynecological, pan-kidney, and pan-squamous cancers, and those related by stemness features, which in turn may inform strategies for future therapeutic development.

Perspective on Oncogenic Processes at the End of the Beginning of Cancer Genomics.

The Cancer Genome Atlas (TCGA) has catalyzed systematic characterization of diverse genomic alterations underlying human cancers. At this historic junction marking the completion of genomic characterization of over 11,000 tumors from 33 cancer types, we present our current understanding of the molecular processes governing oncogenesis. We illustrate our insights into cancer through synthesis of the findings of the TCGA PanCancer Atlas project on three facets of oncogenesis: (1) somatic driver mutations, germline pathogenic variants, and their interactions in the tumor; (2) the influence of the tumor genome and epigenome on transcriptome and proteome; and (3) the relationship between tumor and the microenvironment, including implications for drugs targeting driver events and immunotherapies. These results will anchor future characterization of rare and common tumor types, primary and relapsed tumors, and cancers across ancestry groups and will guide the deployment of clinical genomic sequencing.

Germline genetic contribution to the immune landscape of cancer.

Understanding the contribution of the host's genetic background to cancer immunity may lead to improved stratification for immunotherapy and to the identification of novel therapeutic targets. We investigated the effect of common and rare germline variants on 139 well-defined immune traits in ∼9000 cancer patients enrolled in TCGA. High heritability was observed for estimates of NK cell and T cell subset infiltration and for interferon signaling. Common variants of IFIH1, TMEM173 (STING1), and TMEM108 were associated with differential interferon signaling and variants mapping to RBL1 correlated with T cell subset abundance. Pathogenic or likely pathogenic variants in BRCA1 and in genes involved in telomere stabilization and Wnt-ß-catenin also acted as immune modulators. Our findings provide evidence for the impact of germline genetics on the composition and functional orientation of the tumor immune microenvironment. The curated datasets, variants, and genes identified provide a resource toward further understanding of tumor-immune interactions.

The Immune Landscape of Cancer.

We performed an extensive immunogenomic analysis of more than 10,000 tumors comprising 33 diverse cancer types by utilizing data compiled by TCGA. Across cancer types, we identified six immune subtypes-wound healing, IFN-γ dominant, inflammatory, lymphocyte depleted, immunologically quiet, and TGF-ß dominant-characterized by differences in macrophage or lymphocyte signatures, Th1:Th2 cell ratio, extent of intratumoral heterogeneity, aneuploidy, extent of neoantigen load, overall cell proliferation, expression of immunomodulatory genes, and prognosis. Specific driver mutations correlated with lower (CTNNB1, NRAS, or IDH1) or higher (BRAF, TP53, or CASP8) leukocyte levels across all cancers. Multiple control modalities of the intracellular and extracellular networks (transcription, microRNAs, copy number, and epigenetic processes) were involved in tumor-immune cell interactions, both across and within immune subtypes. Our immunogenomics pipeline to characterize these heterogeneous tumors and the resulting data are intended to serve as a resource for future targeted studies to further advance the field.

Validation and calibration of next-generation sequencing to identify Epstein-Barr virus-positive gastric cancer in The Cancer Genome Atlas.

The Epstein-Barr virus (EBV)-positive subtype of gastric adenocarcinoma is conventionally identified by in situ hybridization (ISH) for viral nucleic acids, but next-generation sequencing represents a potential alternative. We therefore determined normalized EBV read counts by whole-genome, whole-exome, mRNA and miRNA sequencing for 295 fresh-frozen gastric tumor samples. Formalin-fixed, paraffin-embedded tissue sections were retrieved for ISH confirmation of 13 high-EBV and 11 low-EBV cases. In pairwise comparisons, individual samples were either concordantly high or concordantly low by all genomic methods for which data were available. Empiric cutoffs of sequencing counts identified 26 (9 %) tumors as EBV positive. EBV positivity or negativity by molecular testing was confirmed by EBER-ISH in all but one tumor evaluated by both approaches (kappa = 0.91). EBV-positive gastric tumors can be accurately identified by quantifying viral sequences in genomic data. Simultaneous analyses of human and viral DNA, mRNA and miRNA could streamline tumor profiling for clinical care and research.

Systems biology approaches identify ATF3 as a negative regulator of Toll-like receptor 4.

The innate immune system is absolutely required for host defence, but, uncontrolled, it leads to inflammatory disease. This control is mediated, in part, by cytokines that are secreted by macrophages. Immune regulation is extraordinarily complex, and can be best investigated with systems approaches (that is, using computational tools to predict regulatory networks arising from global, high-throughput data sets). Here we use cluster analysis of a comprehensive set of transcriptomic data derived from Toll-like receptor (TLR)-activated macrophages to identify a prominent group of genes that appear to be regulated by activating transcription factor 3 (ATF3), a member of the CREB/ATF family of transcription factors. Network analysis predicted that ATF3 is part of a transcriptional complex that also contains members of the nuclear factor (NF)-kappaB family of transcription factors. Promoter analysis of the putative ATF3-regulated gene cluster demonstrated an over-representation of closely apposed ATF3 and NF-kappaB binding sites, which was verified by chromatin immunoprecipitation and hybridization to a DNA microarray. This cluster included important cytokines such as interleukin (IL)-6 and IL-12b. ATF3 and Rel (a component of NF-kappaB) were shown to bind to the regulatory regions of these genes upon macrophage activation. A kinetic model of Il6 and Il12b messenger RNA expression as a function of ATF3 and NF-kappaB promoter binding predicted that ATF3 is a negative regulator of Il6 and Il12b transcription, and this hypothesis was validated using Atf3-null mice. ATF3 seems to inhibit Il6 and Il12b transcription by altering chromatin structure, thereby restricting access to transcription factors. Because ATF3 is itself induced by lipopolysaccharide, it seems to regulate TLR-stimulated inflammatory responses as part of a negative-feedback loop.

The immune landscape of solid pediatric tumors.

BACKGROUND: Large immunogenomic analyses have demonstrated the prognostic role of the functional orientation of the tumor microenvironment in adult solid tumors, this variable has been poorly explored in the pediatric counterpart. METHODS: We performed a systematic analysis of public RNAseq data (TARGET) for five pediatric tumor types (408 patients): Wilms tumor (WLM), neuroblastoma (NBL), osteosarcoma (OS), clear cell sarcoma of the kidney (CCSK) and rhabdoid tumor of the kidney (RT). We assessed the performance of the Immunologic Constant of Rejection (ICR), which captures an active Th1/cytotoxic response. We also performed gene set enrichment analysis (ssGSEA) and clustered more than 100 well characterized immune traits to define immune subtypes and compared their outcome. RESULTS: A higher ICR score was associated with better survival in OS and high risk NBL without MYCN amplification but with poorer survival in WLM. Clustering of immune traits revealed the same five principal modules previously described in adult tumors (TCGA). These modules divided pediatric patients into six immune subtypes (S1-S6) with distinct survival outcomes. The S2 cluster showed the best overall survival, characterized by low enrichment of the wound healing signature, high Th1, and low Th2 infiltration, while the reverse was observed in S4. Upregulation of the WNT/Beta-catenin pathway was associated with unfavorable outcomes and decreased T-cell infiltration in OS. CONCLUSIONS: We demonstrated that extracranial pediatric tumors could be classified according to their immune disposition, unveiling similarities with adults' tumors. Immunological parameters might be explored to refine diagnostic and prognostic biomarkers and to identify potential immune-responsive tumors.

Analytic pipelines to assess the relationship between immune response and germline genetics in human tumors.

Germline genetic variants modulate human immune response. We present analytical pipelines for assessing the contribution of hosts' genetic background to the immune landscape of solid tumors using harmonized data from more than 9,000 patients in The Cancer Genome Atlas (TCGA). These include protocols for heritability, genome-wide association studies (GWAS), colocalization, and rare variant analyses. These workflows are developed around the structure of TCGA but can be adapted to explore other repositories or in the context of cancer immunotherapy. For complete details on the use and execution of this protocol, please refer to Sayaman et al. (2021).

Patient-Specific Cell Communication Networks Associate With Disease Progression in Cancer.

The maintenance and function of tissues in health and disease depends on cell-cell communication. This work shows how high-level features, representing cell-cell communication, can be defined and used to associate certain signaling "axes" with clinical outcomes. We generated a scaffold of cell-cell interactions and defined a probabilistic method for creating per-patient weighted graphs based on gene expression and cell deconvolution results. With this method, we generated over 9,000 graphs for The Cancer Genome Atlas (TCGA) patient samples, each representing likely channels of intercellular communication in the tumor microenvironment (TME). It was shown that cell-cell edges were strongly associated with disease severity and progression, in terms of survival time and tumor stage. Within individual tumor types, there are predominant cell types, and the collection of associated edges were found to be predictive of clinical phenotypes. Additionally, genes associated with differentially weighted edges were enriched in Gene Ontology terms associated with tissue structure and immune response. Code, data, and notebooks are provided to enable the application of this method to any expression dataset (https://github.com/IlyaLab/Pan-Cancer-Cell-Cell-Comm-Net).

Identification of tuberculosis susceptibility genes with human macrophage gene expression profiles.

Although host genetics influences susceptibility to tuberculosis (TB), few genes determining disease outcome have been identified. We hypothesized that macrophages from individuals with different clinical manifestations of Mycobacterium tuberculosis (Mtb) infection would have distinct gene expression profiles and that polymorphisms in these genes may also be associated with susceptibility to TB. We measured gene expression levels of >38,500 genes from ex vivo Mtb-stimulated macrophages in 12 subjects with 3 clinical phenotypes: latent, pulmonary, and meningeal TB (n = 4 per group). After identifying differentially expressed genes, we confirmed these results in 34 additional subjects by real-time PCR. We also used a case-control study design to examine whether polymorphisms in differentially regulated genes were associated with susceptibility to these different clinical forms of TB. We compared gene expression profiles in Mtb-stimulated and unstimulated macrophages and identified 1,608 and 199 genes that were differentially expressed by >2- and >5-fold, respectively. In an independent sample set of 34 individuals and a subset of highly regulated genes, 90% of the microarray results were confirmed by RT-PCR, including expression levels of CCL1, which distinguished the 3 clinical groups. Furthermore, 6 single nucleotide polymorphisms (SNPs) in CCL1 were found to be associated with TB in a case-control genetic association study with 273 TB cases and 188 controls. To our knowledge, this is the first identification of CCL1 as a gene involved in host susceptibility to TB and the first study to combine microarray and DNA polymorphism studies to identify genes associated with TB susceptibility. These results suggest that genome-wide studies can provide an unbiased method to identify critical macrophage response genes that are associated with different clinical outcomes and that variation in innate immune response genes regulate susceptibility to TB.

Multiplatform Integrative Analysis of Immunogenomic Data for Biomarker Discovery.

Genomic platforms are increasingly used to gain an understanding of how the immune system responds to cancer. In this chapter we describe steps applied in immunogenomic processing and data processing from multiple genomics platforms to enable study of immune response and the evaluation of candidate biomarkers. We also describe how publicly available web resources can be used to discover and evaluate candidate cancer immune biomarkers.

CRI iAtlas: an interactive portal for immuno-oncology research.

The Cancer Research Institute (CRI) iAtlas is an interactive web platform for data exploration and discovery in the context of tumors and their interactions with the immune microenvironment. iAtlas allows researchers to study immune response characterizations and patterns for individual tumor types, tumor subtypes, and immune subtypes. iAtlas supports computation and visualization of correlations and statistics among features related to the tumor microenvironment, cell composition, immune expression signatures, tumor mutation burden, cancer driver mutations, adaptive cell clonality, patient survival, expression of key immunomodulators, and tumor infiltrating lymphocyte (TIL) spatial maps. iAtlas was launched to accompany the release of the TCGA PanCancer Atlas and has since been expanded to include new capabilities such as (1) user-defined loading of sample cohorts, (2) a tool for classifying expression data into immune subtypes, and (3) integration of TIL mapping from digital pathology images. We expect that the CRI iAtlas will accelerate discovery and improve patient outcomes by providing researchers access to standardized immunogenomics data to better understand the tumor immune microenvironment and its impact on patient responses to immunotherapy.

Society for Immunotherapy of Cancer clinical and biomarkers data sharing resource document: Volume II-practical challenges.

The development of strongly predictive validated biomarkers is essential for the field of immuno-oncology (IO) to advance. The highly complex, multifactorial data sets required to develop these biomarkers necessitate effective, responsible data-sharing efforts in order to maximize the scientific knowledge and utility gained from their collection. While the sharing of clinical- and safety-related trial data has already been streamlined to a large extent, the sharing of biomarker-aimed clinical trial derived data and data sets has been met with a number of hurdles that have impaired the progression of biomarkers from hypothesis to clinical use. These hurdles include technical challenges associated with the infrastructure, technology, workforce, and sustainability required for clinical biomarker data sharing. To provide guidance and assist in the navigation of these challenges, the Society for Immunotherapy of Cancer (SITC) Biomarkers Committee convened to outline the challenges that researchers currently face, both at the conceptual level (Volume I) and at the technical level (Volume II). The committee also suggests possible solutions to these problems in the form of professional standards and harmonized requirements for data sharing, assisting in continued progress toward effective, clinically relevant biomarkers in the IO setting.

Society for Immunotherapy of Cancer clinical and biomarkers data sharing resource document: Volume I-conceptual challenges.

The sharing of clinical trial data and biomarker data sets among the scientific community, whether the data originates from pharmaceutical companies or academic institutions, is of critical importance to enable the development of new and improved cancer immunotherapy modalities. Through data sharing, a better understanding of current therapies in terms of their efficacy, safety and biomarker data profiles can be achieved. However, the sharing of these data sets involves a number of stakeholder groups including patients, researchers, private industry, scientific journals and professional societies. Each of these stakeholder groups has differing interests in the use and sharing of clinical trial and biomarker data, and the conflicts caused by these differing interests represent significant obstacles to effective, widespread sharing of data. Thus, the Society for Immunotherapy of Cancer (SITC) Biomarkers Committee convened to identify the current barriers to biomarker data sharing in immuno-oncology (IO) and to help in establishing professional standards for the responsible sharing of clinical trial data. The conclusions of the committee are described in two position papers: Volume I-conceptual challenges and Volume II-practical challenges, the first of which is presented in this manuscript. Additionally, the committee suggests actions by key stakeholders in the field (including organizations and professional societies) as the best path forward, encouraging the cultural shift needed to ensure responsible data sharing in the IO research setting.

Uncovering a macrophage transcriptional program by integrating evidence from motif scanning and expression dynamics.

Macrophages are versatile immune cells that can detect a variety of pathogen-associated molecular patterns through their Toll-like receptors (TLRs). In response to microbial challenge, the TLR-stimulated macrophage undergoes an activation program controlled by a dynamically inducible transcriptional regulatory network. Mapping a complex mammalian transcriptional network poses significant challenges and requires the integration of multiple experimental data types. In this work, we inferred a transcriptional network underlying TLR-stimulated murine macrophage activation. Microarray-based expression profiling and transcription factor binding site motif scanning were used to infer a network of associations between transcription factor genes and clusters of co-expressed target genes. The time-lagged correlation was used to analyze temporal expression data in order to identify potential causal influences in the network. A novel statistical test was developed to assess the significance of the time-lagged correlation. Several associations in the resulting inferred network were validated using targeted ChIP-on-chip experiments. The network incorporates known regulators and gives insight into the transcriptional control of macrophage activation. Our analysis identified a novel regulator (TGIF1) that may have a role in macrophage activation.

MITF has a central role in regulating starvation-induced autophagy in melanoma.

The MITF transcription factor is a master regulator of melanocyte development and a critical factor in melanomagenesis. The related transcription factors TFEB and TFE3 regulate lysosomal activity and autophagy processes known to be important in melanoma. Here we show that MITF binds the CLEAR-box element in the promoters of lysosomal and autophagosomal genes in melanocytes and melanoma cells. The crystal structure of MITF bound to the CLEAR-box reveals how the palindromic nature of this motif induces symmetric MITF homodimer binding. In metastatic melanoma tumors and cell lines, MITF positively correlates with the expression of lysosomal and autophagosomal genes, which, interestingly, are different from the lysosomal and autophagosomal genes correlated with TFEB and TFE3. Depletion of MITF in melanoma cells and melanocytes attenuates the response to starvation-induced autophagy, whereas the overexpression of MITF in melanoma cells increases the number of autophagosomes but is not sufficient to induce autophagic flux. Our results suggest that MITF and the related factors TFEB and TFE3 have separate roles in regulating a starvation-induced autophagy response in melanoma. Understanding the normal and pathophysiological roles of MITF and related transcription factors may provide important clinical insights into melanoma therapy.

Toward a comprehensive view of cancer immune responsiveness: a synopsis from the SITC workshop.

Tumor immunology has changed the landscape of cancer treatment. Yet, not all patients benefit as cancer immune responsiveness (CIR) remains a limitation in a considerable proportion of cases. The multifactorial determinants of CIR include the genetic makeup of the patient, the genomic instability central to cancer development, the evolutionary emergence of cancer phenotypes under the influence of immune editing, and external modifiers such as demographics, environment, treatment potency, co-morbidities and cancer-independent alterations including immune homeostasis and polymorphisms in the major and minor histocompatibility molecules, cytokines, and chemokines. Based on the premise that cancer is fundamentally a disorder of the genes arising within a cell biologic process, whose deviations from normality determine the rules of engagement with the host's response, the Society for Immunotherapy of Cancer (SITC) convened a task force of experts from various disciplines including, immunology, oncology, biophysics, structural biology, molecular and cellular biology, genetics, and bioinformatics to address the complexity of CIR from a holistic view. The task force was launched by a workshop held in San Francisco on May 14-15, 2018 aimed at two preeminent goals: 1) to identify the fundamental questions related to CIR and 2) to create an interactive community of experts that could guide scientific and research priorities by forming a logical progression supported by multiple perspectives to uncover mechanisms of CIR. This workshop was a first step toward a second meeting where the focus would be to address the actionability of some of the questions identified by working groups. In this event, five working groups aimed at defining a path to test hypotheses according to their relevance to human cancer and identifying experimental models closest to human biology, which include: 1) Germline-Genetic, 2) Somatic-Genetic and 3) Genomic-Transcriptional contributions to CIR, 4) Determinant(s) of Immunogenic Cell Death that modulate CIR, and 5) Experimental Models that best represent CIR and its conversion to an immune responsive state. This manuscript summarizes the contributions from each group and should be considered as a first milestone in the path toward a more contemporary understanding of CIR. We appreciate that this effort is far from comprehensive and that other relevant aspects related to CIR such as the microbiome, the individual's recombined T cell and B cell receptors, and the metabolic status of cancer and immune cells were not fully included. These and other important factors will be included in future activities of the taskforce. The taskforce will focus on prioritization and specific actionable approach to answer the identified questions and implementing the collaborations in the follow-up workshop, which will be held in Houston on September 4-5, 2019.