The Cancer Genome Atlas (TCGA) research network has made public a large collection of clinical and molecular phenotypes of more than 10 000 tumor patients across 33 different tumor types. Using this cohort, TCGA has published over 20 marker papers detailing the genomic and epigenomic alterations associated with these tumor types. Although many important discoveries have been made by TCGA's research network, opportunities still exist to implement novel methods, thereby elucidating new biological pathways and diagnostic markers. However, mining the TCGA data presents several bioinformatics challenges, such as data retrieval and integration with clinical data and other molecular data types (e.g. RNA and DNA methylation). We developed an R/Bioconductor package called TCGAbiolinks to address these challenges and offer bioinformatics solutions by using a guided workflow to allow users to query, download and perform integrative analyses of TCGA data. We combined methods from computer science and statistics into the pipeline and incorporated methodologies developed in previous TCGA marker studies and in our own group. Using four different TCGA tumor types (Kidney, Brain, Breast and Colon) as examples, we provide case studies to illustrate examples of reproducibility, integrative analysis and utilization of different Bioconductor packages to advance and accelerate novel discoveries.

SUMMARY Therapy development for adult diffuse glioma is hindered by incomplete knowledge of somatic glioma driving alterations and suboptimal disease classification. We defined the complete set of genes associated with 1,122 diffuse grade II-III-IV gliomas from The Cancer Genome Atlas and used molecular profiles to improve disease classification, identify molecular correlations, and provide insights into the progression from low- to high-grade disease. Whole genome sequencing data analysis determined that ATRX but not TERT promoter mutations are associated with increased telomere length. Recent advances in glioma classification based on IDH mutation and 1p/19q co-deletion status were recapitulated through analysis of DNA methylation profiles, which identified clinically relevant molecular subsets. A subtype of IDH-mutant glioma was associated with DNA demethylation and poor outcome; a group of IDH-wildtype diffuse glioma showed molecular similarity to pilocytic astrocytoma and relatively favorable survival. Understanding of cohesive disease groups may aid improved clinical outcomes.

Oncogenic TACC-tics Human cancers exhibit many types of genomic rearrangements—including some that juxtapose sequences from two unrelated genes—thereby creating fusion proteins with oncogenic activity. Functional analysis of these fusion genes can provide mechanistic insights into tumorigenesis and potentially lead to effective drugs, as famously illustrated by the BCR-ABL gene in chronic myelogenous leukemia. Singh et al. (p. 1231 , published online 26 July) identify and characterize a fusion gene present in 3% of human glioblastomas, a deadly brain cancer. In the resultant fusion protein, the tyrosine kinase region of the fibroblast growth factor receptor (FGFR) is joined to a domain from a transforming acidic coiled-coil (TACC) protein. The TACC-FGFR protein is oncogenic, shows unregulated kinase activity, localizes to the mitotic spindle, and disrupts chromosome segregation. In mice, FGFR inhibitors slowed the growth of tumors driven by the TACC-FGFR gene, suggesting that a subset of glioblastoma patients may benefit from these types of drugs.

Summary The molecular characterization of immune subsets is important for designing effective strategies to understand and treat diseases. We characterized 29 immune cell types within the peripheral blood mononuclear cell (PBMC) fraction of healthy donors using RNA-seq (RNA sequencing) and flow cytometry. Our dataset was used, first, to identify sets of genes that are specific, are co-expressed, and have housekeeping roles across the 29 cell types. Then, we examined differences in mRNA heterogeneity and mRNA abundance revealing cell type specificity. Last, we performed absolute deconvolution on a suitable set of immune cell types using transcriptomics signatures normalized by mRNA abundance. Absolute deconvolution is ready to use for PBMC transcriptomic data using our Shiny app ( https://github.com/giannimonaco/ABIS ). We benchmarked different deconvolution and normalization methods and validated the resources in independent cohorts. Our work has research, clinical, and diagnostic value by making it possible to effectively associate observations in bulk transcriptomics data to specific immune subsets.

No abstract

We examined how the immune microenvironment molds tumor evolution at different metastatic organs in a longitudinal dataset of colorectal cancer. Through multiplexed analyses, we showed that clonal evolution patterns during metastatic progression depend on the immune contexture at the metastatic site. Genetic evidence of neoantigen depletion was observed in the sites with high Immunoscore and spatial proximity between Ki67 + tumor cells and CD3 + cells. The immunoedited tumor clones were eliminated and did not recur, while progressing clones were immune privileged, despite the presence of tumor-infiltrating lymphocytes. Characterization of immune-privileged metastases revealed tumor-intrinsic and tumorextrinsic mechanisms of escape. The lowest recurrence risk was associated with high Immunoscore, occurrence of immunoediting, and low tumor burden. We propose a parallel selection model of metastatic progression, where branched evolution could be traced back to immune-escaping clones. The findings could inform the understanding of cancer dissemination and the development of immunotherapeutics.

No abstract

BackgroundDespite extensive molecular characterization, we lack a comprehensive understanding of lineage identity, differentiation, and proliferation in high-grade gliomas (HGGs).MethodsWe sampled the cellular milieu of HGGs by profiling dissociated human surgical specimens with a high-density microwell system for massively parallel single-cell RNA-Seq. We analyzed the resulting profiles to identify subpopulations of both HGG and microenvironmental cells and applied graph-based methods to infer structural features of the malignantly transformed populations.ResultsWhile HGG cells can resemble glia or even immature neurons and form branched lineage structures, mesenchymal transformation results in unstructured populations. Glioma cells in a subset of mesenchymal tumors lose their neural lineage identity, express inflammatory genes, and co-exist with marked myeloid infiltration, reminiscent of molecular interactions between glioma and immune cells established in animal models. Additionally, we discovered a tight coupling between lineage resemblance and proliferation among malignantly transformed cells. Glioma cells that resemble oligodendrocyte progenitors, which proliferate in the brain, are often found in the cell cycle. Conversely, glioma cells that resemble astrocytes, neuroblasts, and oligodendrocytes, which are non-proliferative in the brain, are generally non-cycling in tumors.ConclusionsThese studies reveal a relationship between cellular identity and proliferation in HGG and distinct population structures that reflects the extent of neural and non-neural lineage resemblance among malignantly transformed cells.Electronic supplementary materialThe online version of this article (10.1186/s13073-018-0567-9) contains supplementary material, which is available to authorized users.