Search

Abstract Background Neoantigen-targeting therapies including personalized vaccines have shown promise in the treatment of cancers, particularly when used in combination with checkpoint blockade therapy. At least 100 clinical trials involving these therapies have been initiated globally. Accurate identification and prioritization of neoantigens is crucial for designing these trials, predicting treatment response, and understanding mechanisms of resistance. With the advent of massively parallel DNA and RNA sequencing technologies, it is now possible to computationally predict neoantigens based on patient-specific variant information. However, numerous factors must be considered when prioritizing neoantigens for use in personalized therapies. Complexities such as alternative transcript annotations, various binding, presentation and immunogenicity prediction algorithms, and variable peptide lengths/registers all potentially impact the neoantigen selection process. There has been a rapid development of computational tools that attempt to account for these complexities. While these tools generate numerous algorithmic predictions for neoantigen characterization, results from these pipelines are difficult to navigate and require extensive knowledge of the underlying tools for accurate interpretation. This often leads to over-simplification of pipeline outputs to make them tractable, for example, limiting prediction to a single RNA isoform or only summarizing the top ranked of many possible peptide candidates. In addition to variant detection, gene expression, and predicted peptide binding affinities, recent studies have also demonstrated the importance of mutation location, allele-specific anchor locations, and variation of T-cell response to long versus short peptides. Due to the intricate nature and number of salient neoantigen features, presenting all relevant information to facilitate candidate selection for downstream applications is a difficult challenge that current tools fail to address. Results We have created pVACview, the first interactive tool designed to aid in the prioritization and selection of neoantigen candidates for personalized neoantigen therapies including cancer vaccines. pVACview has a user-friendly and intuitive interface where users can upload, explore, select, and export their neoantigen candidates. The tool allows users to visualize candidates at multiple levels of detail including variant, transcript, peptide, and algorithm prediction information. Conclusions pVACview will allow researchers to analyze and prioritize neoantigen candidates with greater efficiency and accuracy in basic and translational settings. The application is available as part of the pVACtools software at pvactools.org and as an online server at pvacview.org.

Abstract Background Neoantigen vaccines can induce or enhance highly specific antitumor immune responses with minimal risk of autoimmunity. We have developed a neoantigen DNA vaccine platform capable of efficiently presenting both HLA class I and II epitopes and performed a phase 1 clinical trial in triple-negative breast cancer patients with persistent disease on surgical pathology following neoadjuvant chemotherapy, a patient population at high risk of disease recurrence. Methods Expressed somatic mutations were identified by tumor/normal exome sequencing and tumor RNA sequencing. The pVACtools software suite of neoantigen prediction algorithms was used to identify and prioritize cancer neoantigens and facilitate vaccine design for manufacture in an academic GMP facility. Neoantigen DNA vaccines were administered via electroporation in the adjuvant setting (i.e., following surgical removal of the primary tumor and completion of standard of care therapy). Vaccines were monitored for safety and immune responses via ELISpot, intracellular cytokine production via flow cytometry, and TCR sequencing. Results Eighteen subjects received three doses of a neoantigen DNA vaccine encoding on average 11 neoantigens per patient (range 4–20). The vaccinations were well tolerated with relatively few adverse events. Neoantigen-specific T cell responses were induced in 14/18 patients as measured by ELISpot and flow cytometry. At a median follow-up of 36 months, recurrence-free survival was 87.5% (95% CI: 72.7–100%) in the cohort of vaccinated patients. Conclusion Our study demonstrates neoantigen DNA vaccines are safe, feasible, and capable of inducing neoantigen-specific immune responses. Clinical trial registration number NCT02348320.

Abstract Spontaneous cancers in companion dogs are robust models of human disease. Tracking tumor-specific immune responses in these models requires reagents to perform species-specific single cell T cell receptor sequencing (scTCRseq). scTCRseq and integration with scRNA data have not been demonstrated on companion dogs with cancer. Here, five healthy dogs, two dogs with T cell lymphoma and four dogs with melanoma are selected to demonstrate applicability of scTCRseq in a cancer immunotherapy setting. Single-cell suspensions of PBMCs or lymph node aspirates are profiled using scRNA and dog-specific scTCRseq primers. In total, 77,809 V(D)J-expressing cells are detected, with an average of 3498 (348 - 5,971) unique clonotypes identified per sample. In total, 29/34, 40/40, 22/22 and 9/9 known functional TRAV, TRAJ, TRBV and TRBJ gene segments are observed respectively. Pseudogene or otherwise defective gene segments are also detected supporting re-annotation of several as functional. Healthy dogs exhibit highly diverse repertoires, T cell lymphomas exhibit clonal repertoires, and vaccine-treated melanoma dogs are dominated by a small number of highly abundant clonotypes. scRNA libraries define large clusters of V(D)J-expressing CD8+ and CD4 + T cells. Dominant clonotypes observed in melanoma PBMCs are predominantly CD8 + T cells, with activated phenotypes, suggesting possible anti-tumor T cell populations.

Abstract Background Although macrophages are now recognized as an essential part of the HIV latent reservoir, whether and how viral latency is established and reactivated in these cell types is poorly understood. To understand the fundamental mechanisms of viral latency in macrophages, there is an urgent need to develop latency models amenable to genetic manipulations and screening for appropriate latency-reversing agents (LRAs). Given that differentiated THP-1 cells resemble monocyte-derived macrophages in HIV replication mechanisms, we set out to establish a macrophage cell model for HIV latency using THP-1 cells. Methods We created single-cell clones of THP-1 cells infected with a single copy of the dual-labeled HIVGKO in which a codon switched eGFP (csGFP) is under the control of the HIV-1 5’ LTR promoter, and a monomeric Kusabira orange 2 (mKO2) under the control of cellular elongation factor one alpha promoter (EF1α). Latently infected cells are csGFP−, mKO2+, while cells with actively replicating HIV (or reactivated virus) are csGFP+,mKO2+. After sorting for latently infected cells, each of the THP-1 clones with unique integration sites for HIV was differentiated into macrophage-like cells with phorbol 12-myristate 13-acetate (PMA) and treated with established LRAs to stimulate HIV reactivation. Monocyte-derived macrophages (MDMs) harboring single copies of HIVGKO were used to confirm our findings. Results We obtained clones of THP-1 cells with latently infected HIV with unique integration sites. When the differentiated THP-1 or primary MDMs cells were treated with various LRAs, the bromodomain inhibitors JQ1 and I-BET151 were the most potent compounds. Knockdown of BRD4, the target of JQ1, resulted in increased reactivation, thus confirming the pharmacological effect. The DYRK1A inhibitor Harmine and lipopolysaccharide (LPS) also showed significant reactivation across all three MDM donors. Remarkably, LRAs like PMA/ionomycin, bryostatin-1, and histone deacetylase inhibitors known to potently reactivate latent HIV in CD4 + T cells showed little activity in macrophages. Conclusions Our results indicate that this model could be used to screen for appropriate LRAs for macrophages and show that HIV latency and reactivation mechanisms in macrophages may be distinct from those of CD4 + T cells.

Summary: While most dizygotic twins have a dichorionic placenta, rare cases of dizygotic twins with a monochorionic placenta have been reported. The monochorionic placenta in dizygotic twins allows in utero exchange of embryonic cells, resulting in chimerism in the twins. In practice, this chimerism is incidentally identified in mixed ABO blood types or in the presence of cells with a discordant sex chromosome. Here, we applied whole-genome sequencing to one triplet and one twin family to precisely understand their zygotic compositions, using millions of genomic variants as barcodes of zygotic origins. Peripheral blood showed asymmetrical contributions from two sister zygotes, where one of the zygotes was the major clone in both twins. Single-cell RNA sequencing of peripheral blood tissues further showed differential contributions from the two sister zygotes across blood cell types. In contrast, buccal tissues were pure in genetic composition, suggesting that in utero cellular exchanges were confined to the blood tissues. Our study illustrates the cellular history of twinning during human development, which is critical for managing the health of chimeric individuals in the era of genomic medicine.

Analysing the regulatory consequences of mutations and splice variants at large scale in cancer requires efficient computational tools. Here, the authors develop RegTools, a software package that can identify splice-associated variants from large-scale genomics and transcriptomics data with efficiency and flexibility.

Summary: To explore mechanisms of response to combined PD-1/CTLA-4 immune checkpoint blockade (ICB) treatment in individual cell types, we generated scRNA-seq using a mouse model of invasive urothelial carcinoma with three conditions: untreated tumor, treated tumor, and tumor treated after CD4+ T cell depletion. After classifying tumor cells based on detection of somatic variants and assigning non-tumor cell types using SingleR, we performed differential expression analysis, overrepresentation analysis, and gene set enrichment analysis (GSEA) within each cell type. GSEA revealed that endothelial cells were enriched for upregulated IFN-g response genes when comparing treated cells to both untreated cells and cells treated after CD4+ T cell depletion. Functional analysis showed that knocking out IFNgR1 in endothelial cells inhibited treatment response. Together, these results indicated that IFN-g signaling in endothelial cells is a key mediator of ICB induced anti-tumor activity.

Abstract Circulating tumor DNA (ctDNA) in peripheral blood has been used to predict prognosis and therapeutic response for triple-negative breast cancer (TNBC) patients. However, previous approaches typically use large comprehensive panels of genes commonly mutated across all breast cancers. Given the reduction in sequencing costs and decreased turnaround times associated with panel generation, the objective of this study was to assess the use of custom micro-panels for tracking disease and predicting clinical outcomes for patients with TNBC. Paired tumor-normal samples from patients with TNBC were obtained at diagnosis (T0) and whole exome sequencing (WES) was performed to identify somatic variants associated with individual tumors. Custom micro-panels of 4–6 variants were created for each individual enrolled in the study. Peripheral blood was obtained at baseline, during Cycle 1 Day 3, at time of surgery, and in 3–6 month intervals after surgery to assess variant allele fraction (VAF) at different timepoints during disease course. The VAF was compared to clinical outcomes to evaluate the ability of custom micro-panels to predict pathological response, disease-free intervals, and patient relapse. A cohort of 50 individuals were evaluated for up to 48 months post-diagnosis of TNBC. In total, there were 33 patients who did not achieve pathological complete response (pCR) and seven patients developed clinical relapse. For all patients who developed clinical relapse and had peripheral blood obtained ≤ 6 months prior to relapse (n = 4), the custom ctDNA micro-panels identified molecular relapse at an average of 4.3 months prior to clinical relapse. The custom ctDNA panel results were moderately associated with pCR such that during disease monitoring, only 11% of patients with pCR had a molecular relapse, whereas 47% of patients without pCR had a molecular relapse (Chi-Square; p-value = 0.10). In this study, we show that a custom micro-panel of 4–6 markers can be effectively used to predict outcomes and monitor remission for patients with TNBC. These custom micro-panels show high sensitivity for detecting molecular relapse in advance of clinical relapse. The use of these panels could improve patient outcomes through early detection of relapse with preemptive intervention prior to symptom onset.

Since the T-box transcription factors (TFs) T-BET and EOMES are necessary for initiation of NK cell development, their ongoing requirement for mature NK cell homeostasis, function, and molecular programming remains unclear. To address this, T-BET and EOMES were deleted in unexpanded primary human NK cells using CRISPR/Cas9. Deleting these TFs compromised in vivo antitumor response of human NK cells. Mechanistically, T-BET and EOMES were required for normal NK cell proliferation and persistence in vivo. NK cells lacking T-BET and EOMES also exhibited defective responses to cytokine stimulation. Single-cell RNA-Seq revealed a specific T-box transcriptional program in human NK cells, which was rapidly lost following T-BET and EOMES deletion. Further, T-BET– and EOMES-deleted CD56bright NK cells acquired an innate lymphoid cell precursor–like (ILCP-like) profile with increased expression of the ILC-3–associated TFs RORC and AHR, revealing a role for T-box TFs in maintaining mature NK cell phenotypes and an unexpected role of suppressing alternative ILC lineages. Our study reveals the critical importance of sustained EOMES and T-BET expression to orchestrate mature NK cell function and identity.

Background Interactions between immune and tumor cells are critical to determining cancer progression and response. In addition, preclinical prediction of immune-related drug efficacy is limited by interspecies differences between human and mouse, as well as inter-person germline and somatic variation. To address these gaps, we developed an autologous system that models the tumor microenvironment (TME) from individual patients with solid tumors.Method With patient-derived bone marrow hematopoietic stem and progenitor cells (HSPCs), we engrafted a patient’s hematopoietic system in MISTRG6 mice, followed by transfer of patient-derived xenograft (PDX) tissue, providing a fully genetically matched model to recapitulate the individual’s TME. We used this system to prospectively study tumor-immune interactions in patients with solid tumor.Results Autologous PDX mice generated innate and adaptive immune populations; these cells populated the TME; and tumors from autologously engrafted mice grew larger than tumors from non-engrafted littermate controls. Single-cell transcriptomics revealed a prominent vascular endothelial growth factor A (VEGFA) signature in TME myeloid cells, and inhibition of human VEGF-A abrogated enhanced growth.Conclusions Humanization of the interleukin 6 locus in MISTRG6 mice enhances HSPC engraftment, making it feasible to model tumor-immune interactions in an autologous manner from a bedside bone marrow aspirate. The TME from these autologous tumors display hallmarks of the human TME including innate and adaptive immune activation and provide a platform for preclinical drug testing.

Summary: Although tobacco use is an independent adverse prognostic feature in HPV(+) oropharyngeal squamous cell carcinoma (OPSCC), the biologic features associated with tobacco use have not been systematically investigated. We characterized genomic and immunologic features associated with tobacco use through whole exome sequencing, mRNA hybridization, and immunohistochemical staining in 47 HPV(+) OPSCC tumors. Low expression of transcripts in a T cell-inflamed gene expression profile (TGEP) was associated with tobacco use at diagnosis and lower overall and disease-free survival. Tobacco use was associated with an increased proportion of T > C substitutions and a lower proportion of expected mutational signatures, but not with increases in mutational burden or recurrent oncogenic mutations. Our findings suggest that rather than increased mutational burden, tobacco’s primary and clinically relevant association in HPV(+) OPSCC is immunosuppression of the tumor immune microenvironment. Quantitative assays of T cell infiltration merit further study as prognostic markers in HPV(+) OPSCC.

Abstract Background Precision oncology involves analysis of individual cancer samples to understand the genes and pathways involved in the development and progression of a cancer. To improve patient care, knowledge of diagnostic, prognostic, predisposing, and drug response markers is essential. Several knowledgebases have been created by different groups to collate evidence for these associations. These include the open-access Clinical Interpretation of Variants in Cancer (CIViC) knowledgebase. These databases rely on time-consuming manual curation from skilled experts who read and interpret the relevant biomedical literature. Methods To aid in this curation and provide the greatest coverage for these databases, particularly CIViC, we propose the use of text mining approaches to extract these clinically relevant biomarkers from all available published literature. To this end, a group of cancer genomics experts annotated sentences that discussed biomarkers with their clinical associations and achieved good inter-annotator agreement. We then used a supervised learning approach to construct the CIViCmine knowledgebase. Results We extracted 121,589 relevant sentences from PubMed abstracts and PubMed Central Open Access full-text papers. CIViCmine contains over 87,412 biomarkers associated with 8035 genes, 337 drugs, and 572 cancer types, representing 25,818 abstracts and 39,795 full-text publications. Conclusions Through integration with CIVIC, we provide a prioritized list of curatable clinically relevant cancer biomarkers as well as a resource that is valuable to other knowledgebases and precision cancer analysts in general. All data is publically available and distributed with a Creative Commons Zero license. The CIViCmine knowledgebase is available at http://bionlp.bcgsc.ca/civicmine/.

Summary: Maximizing the personal, public, research, and clinical value of genomic information will require the reliable exchange of genetic variation data. We report here the Variation Representation Specification (VRS, pronounced “verse”), an extensible framework for the computable representation of variation that complements contemporary human-readable and flat file standards for genomic variation representation. VRS provides semantically precise representations of variation and leverages this design to enable federated identification of biomolecular variation with globally consistent and unique computed identifiers. The VRS framework includes a terminology and information model, machine-readable schema, data sharing conventions, and a reference implementation, each of which is intended to be broadly useful and freely available for community use. VRS was developed by a partnership among national information resource providers, public initiatives, and diagnostic testing laboratories under the auspices of the Global Alliance for Genomics and Health (GA4GH).

Summary: The Global Alliance for Genomics and Health (GA4GH) aims to accelerate biomedical advances by enabling the responsible sharing of clinical and genomic data through both harmonized data aggregation and federated approaches. The decreasing cost of genomic sequencing (along with other genome-wide molecular assays) and increasing evidence of its clinical utility will soon drive the generation of sequence data from tens of millions of humans, with increasing levels of diversity. In this perspective, we present the GA4GH strategies for addressing the major challenges of this data revolution. We describe the GA4GH organization, which is fueled by the development efforts of eight Work Streams and informed by the needs of 24 Driver Projects and other key stakeholders. We present the GA4GH suite of secure, interoperable technical standards and policy frameworks and review the current status of standards, their relevance to key domains of research and clinical care, and future plans of GA4GH. Broad international participation in building, adopting, and deploying GA4GH standards and frameworks will catalyze an unprecedented effort in data sharing that will be critical to advancing genomic medicine and ensuring that all populations can access its benefits.

Abstract Manually curated variant knowledgebases and their associated knowledge models are serving an increasingly important role in distributing and interpreting variants in cancer. These knowledgebases vary in their level of public accessibility, and the complexity of the models used to capture clinical knowledge. CIViC (Clinical Interpretation of Variants in Cancer - www.civicdb.org) is a fully open, free-to-use cancer variant interpretation knowledgebase that incorporates highly detailed curation of evidence obtained from peer-reviewed publications and meeting abstracts, and currently holds over 6300 Evidence Items for over 2300 variants derived from over 400 genes. CIViC has seen increased adoption by, and also undertaken collaboration with, a wide range of users and organizations involved in research. To enhance CIViC’s clinical value, regular submission to the ClinVar database and pursuit of other regulatory approvals is necessary. For this reason, a formal peer reviewed curation guideline and discussion of the underlying principles of curation is needed. We present here the CIViC knowledge model, standard operating procedures (SOP) for variant curation, and detailed examples to support community-driven curation of cancer variants.

Abstract Neoantigens are newly formed peptides created from somatic mutations that are capable of inducing tumor-specific T cell recognition. Recently, researchers and clinicians have leveraged next generation sequencing technologies to identify neoantigens and to create personalized immunotherapies for cancer treatment. To create a personalized cancer vaccine, neoantigens must be computationally predicted from matched tumor–normal sequencing data, and then ranked according to their predicted capability in stimulating a T cell response. This candidate neoantigen prediction process involves multiple steps, including somatic mutation identification, HLA typing, peptide processing, and peptide-MHC binding prediction. The general workflow has been utilized for many preclinical and clinical trials, but there is no current consensus approach and few established best practices. In this article, we review recent discoveries, summarize the available computational tools, and provide analysis considerations for each step, including neoantigen prediction, prioritization, delivery, and validation methods. In addition to reviewing the current state of neoantigen analysis, we provide practical guidance, specific recommendations, and extensive discussion of critical concepts and points of confusion in the practice of neoantigen characterization for clinical use. Finally, we outline necessary areas of development, including the need to improve HLA class II typing accuracy, to expand software support for diverse neoantigen sources, and to incorporate clinical response data to improve neoantigen prediction algorithms. The ultimate goal of neoantigen characterization workflows is to create personalized vaccines that improve patient outcomes in diverse cancer types.

Human T-cell leukemia virus type 1 (HTLV-1) is a retrovirus that causes an aggressive T-cell malignancy and a variety of inflammatory conditions. The integrated provirus includes a single binding site for the epigenomic insulator, CCCTC-binding protein (CTCF), but its function remains unclear. In the current study, a mutant virus was examined that eliminates the CTCF-binding site. The mutation did not disrupt the kinetics and levels of virus gene expression, or establishment of or reactivation from latency. However, the mutation disrupted the epigenetic barrier function, resulting in enhanced DNA CpG methylation downstream of the CTCF binding site on both strands of the integrated provirus and H3K4Me3, H3K36Me3, and H3K27Me3 chromatin modifications both up- and downstream of the site. A majority of clonal cell lines infected with wild type HTLV-1 exhibited increased plus strand gene expression with CTCF knockdown, while expression in mutant HTLV-1 clonal lines was unaffected. These findings indicate that CTCF binding regulates HTLV-1 gene expression, DNA and histone methylation in an integration site dependent fashion.

Osteosarcoma is a rare disease in children but is one of the most common cancers in adult large breed dogs. The mutational landscape of both the primary and pulmonary metastatic tumor in two dogs with appendicular osteosarcoma (OSA) was comprehensively evaluated using an automated whole genome sequencing, exome and RNA-seq pipeline that was adapted for this study for use in dogs. Chromosomal lesions were the most common type of mutation. The mutational landscape varied substantially between dogs but the lesions within the same patient were similar. Copy number neutral loss of heterozygosity in mutant TP53 was the most significant driver mutation and involved a large region in the middle of chromosome 5. Canine and human OSA is characterized by loss of cell cycle checkpoint integrity and DNA damage response pathways. Mutational profiling of individual patients with canine OSA would be recommended prior to targeted therapy, given the heterogeneity seen in our study and previous studies.

Wikidata is a community-maintained knowledge base that has been assembled from repositories in the fields of genomics, proteomics, genetic variants, pathways, chemical compounds, and diseases, and that adheres to the FAIR principles of findability, accessibility, interoperability and reusability. Here we describe the breadth and depth of the biomedical knowledge contained within Wikidata, and discuss the open-source tools we have built to add information to Wikidata and to synchronize it with source databases. We also demonstrate several use cases for Wikidata, including the crowdsourced curation of biomedical ontologies, phenotype-based diagnosis of disease, and drug repurposing.

Unravelling the link between somatic mutation and prognosis in estrogen positive (ER+) breast cancer requires the use of long-term follow-up data. Here, combining archival formalin-fixed paraffin embedded tissue and targeted sequencing in three cohorts of ER+ breast cancer, the authors find associations with clinical outcome for NF1 frame-shift nonsense mutations, PIK3R1 mutation, and DDR1 mutations.

Small cell lung cancer (SCLC) patients frequently relapse and become resistant to chemotherapy. Here, the authors analyse the genomic and transcriptomic landscape of primary and relapsed SCLC patients as well as in vitro models, and discover that activation of WNT signalling can drive chemotherapy resistance.

Abstract Many immunotherapies rely on the presence of neoepitopes derived from somatic mutations that lead to altered peptide sequences. Several studies have now analyzed the neoepitope landscape of different cancer subtypes, predominantly for adult samples, which tend to feature significantly higher mutational burden. However, a new report publishing the first comprehensive analysis of the pediatric neoepitope landscape suggests that immunotherapies could also hold promise for pediatric cancers. See related research article 10.1186/s13073-017-0468-3

Estrogen receptor alpha-positive (ERα+) luminal tumors are the most frequent subtype of breast cancer. Stat1−/− mice develop mammary tumors that closely recapitulate the biological characteristics of this cancer subtype. To identify transforming events that contribute to tumorigenesis, we performed whole genome sequencing of Stat1−/− primary mammary tumors and matched normal tissues. This investigation identified somatic truncating mutations affecting the prolactin receptor (PRLR) in all tumor and no normal samples. Targeted sequencing confirmed the presence of these mutations in precancerous lesions, indicating that this is an early event in tumorigenesis. Functional evaluation of these heterozygous mutations in Stat1−/− mouse embryonic fibroblasts showed that co-expression of truncated and wild-type PRLR led to aberrant STAT3 and STAT5 activation downstream of the receptor, cellular transformation in vitro, and tumor formation in vivo. In conclusion, truncating mutations of PRLR promote tumor growth in a model of human ERα+ breast cancer and warrant further investigation.

Aromatase inhibitors are used to treat oestrogen-receptor-positive breast cancer. Here, the authors use genomic approaches to analyse tumours before and after neo-adjuvant treatment and find that treatment alters the clonal landscape of the tumours.

Solid cancers develop within a supportive microenvironment that promotes tumor formation and growth through the elaboration of mitogens and chemokines. Within these tumors, monocytes (macrophages and microglia) represent rich sources of these stromal factors. Leveraging a genetically engineered mouse model of neurofibromatosis type 1 (NF1) low-grade brain tumor (optic glioma), we have previously demonstrated that microglia are essential for glioma formation and maintenance. To identify potential tumor-associated microglial factors that support glioma growth (gliomagens), we initiated a comprehensive large-scale discovery effort using optimized RNA-sequencing methods focused specifically on glioma-associated microglia. Candidate microglial gliomagens were prioritized to identify potential secreted or membrane-bound proteins, which were next validated by quantitative real-time polymerase chain reaction as well as by RNA fluorescence in situ hybridization following minocycline-mediated microglial inactivation in vivo. Using these selection criteria, chemokine (C-C motif) ligand 5 (Ccl5) was identified as a chemokine highly expressed in genetically engineered Nf1 mouse optic gliomas relative to nonneoplastic optic nerves. As a candidate gliomagen, recombinant Ccl5 increased Nf1-deficient optic nerve astrocyte growth in vitro. Importantly, consistent with its critical role in maintaining tumor growth, treatment with Ccl5 neutralizing antibodies reduced Nf1 mouse optic glioma growth and improved retinal dysfunction in vivo. Collectively, these findings establish Ccl5 as an important microglial growth factor for low-grade glioma maintenance relevant to the development of future stroma-targeted brain tumor therapies.

Summary: Cell fate acquisition is heavily influenced by direct interactions between master regulators and tissue-specific enhancers. However, it remains unclear how lineage-specifying transcription factors, which are often expressed in both progenitor and mature cell populations, influence cell differentiation. Using in vivo mouse liver development as a model, we identified thousands of enhancers that are bound by the master regulators HNF4A and FOXA2 in a differentiation-dependent manner, subject to chromatin remodeling, and associated with differentially expressed target genes. Enhancers exclusively occupied in the embryo were found to be responsive to developmentally regulated TEAD2 and coactivator YAP1. Our data suggest that Hippo signaling may affect hepatocyte differentiation by influencing HNF4A and FOXA2 interactions with temporal enhancers. In summary, transcription factor-enhancer interactions are not only tissue specific but also differentiation dependent, which is an important consideration for researchers studying cancer biology or mammalian development and/or using transformed cell lines. : It is unclear how key transcription factors are critical for both lineage specification during embryonic development and maintenance of a differentiated, adult phenotype. By profiling the enhancer occupancy of the key transcription factors HNF4A and FOXA2 during mouse liver development, Alder et al. have found that YAP1 can influence enhancer interactions and target gene expression levels. Enhancer switching enables HNF4A and FOXA2 to fulfill distinct roles during organ development.

The original version of this Article contained errors in the depiction of confidence intervals in the NF1 BCSS data illustrated in Figure 3b. These have now been corrected in both the PDF and HTML versions of the Article. The incorrect version of Figure 3b is presented in the associated Author Correction.

Massively parallel RNA sequencing (RNA-seq) has rapidly become the assay of choice for interrogating RNA transcript abundance and diversity. This article provides a detailed introduction to fundamental RNA-seq molecular biology and informatics concepts. We make available open-access RNA-seq tutorials that cover cloud computing, tool installation, relevant file formats, reference genomes, transcriptome annotations, quality-control strategies, expression, differential expression, and alternative splicing analysis methods. These tutorials and additional training resources are accompanied by complete analysis pipelines and test datasets made available without encumbrance at www.rnaseq.wiki.

In this work, we present the Genome Modeling System (GMS), an analysis information management system capable of executing automated genome analysis pipelines at a massive scale. The GMS framework provides detailed tracking of samples and data coupled with reliable and repeatable analysis pipelines. The GMS also serves as a platform for bioinformatics development, allowing a large team to collaborate on data analysis, or an individual researcher to leverage the work of others effectively within its data management system. Rather than separating ad-hoc analysis from rigorous, reproducible pipelines, the GMS promotes systematic integration between the two. As a demonstration of the GMS, we performed an integrated analysis of whole genome, exome and transcriptome sequencing data from a breast cancer cell line (HCC1395) and matched lymphoblastoid line (HCC1395BL). These data are available for users to test the software, complete tutorials and develop novel GMS pipeline configurations. The GMS is available at https://github.com/genome/gms.

The sensitivity of massively-parallel sequencing has confirmed that most cancers are oligoclonal, with subpopulations of neoplastic cells harboring distinct mutations. A fine resolution view of this clonal architecture provides insight into tumor heterogeneity, evolution, and treatment response, all of which may have clinical implications. Single tumor analysis already contributes to understanding these phenomena. However, cryptic subclones are frequently revealed by additional patient samples (e.g., collected at relapse or following treatment), indicating that accurately characterizing a tumor requires analyzing multiple samples from the same patient. To address this need, we present SciClone, a computational method that identifies the number and genetic composition of subclones by analyzing the variant allele frequencies of somatic mutations. We use it to detect subclones in acute myeloid leukemia and breast cancer samples that, though present at disease onset, are not evident from a single primary tumor sample. By doing so, we can track tumor evolution and identify the spatial origins of cells resisting therapy.

To reveal the clonal architecture of melanoma and associated driver mutations, whole genome sequencing (WGS) and targeted extension sequencing were used to characterize 124 melanoma cases. Significantly mutated gene analysis using 13 WGS cases and 15 additional paired extension cases identified known melanoma genes such as BRAF, NRAS, and CDKN2A, as well as a novel gene EPHA3, previously implicated in other cancer types. Extension studies using tumors from another 96 patients discovered a large number of truncation mutations in tumor suppressors (TP53 and RB1), protein phosphatases (e.g., PTEN, PTPRB, PTPRD, and PTPRT), as well as chromatin remodeling genes (e.g., ASXL3, MLL2, and ARID2). Deep sequencing of mutations revealed subclones in the majority of metastatic tumors from 13 WGS cases. Validated mutations from 12 out of 13 WGS patients exhibited a predominant UV signature characterized by a high frequency of C->T transitions occurring at the 3' base of dipyrimidine sequences while one patient (MEL9) with a hypermutator phenotype lacked this signature. Strikingly, a subclonal mutation signature analysis revealed that the founding clone in MEL9 exhibited UV signature but the secondary clone did not, suggesting different mutational mechanisms for two clonal populations from the same tumor. Further analysis of four metastases from different geographic locations in 2 melanoma cases revealed phylogenetic relationships and highlighted the genetic alterations responsible for differential drug resistance among metastatic tumors. Our study suggests that clonal evaluation is crucial for understanding tumor etiology and drug resistance in melanoma.

Brain tumors (gliomas) contain large populations of infiltrating macrophages and recruited microglia, which in experimental murine glioma models promote tumor formation and progression. Among the barriers to understanding the contributions of these stromal elements to high-grade glioma (glioblastoma; GBM) biology is the relative paucity of tools to characterize infiltrating macrophages and resident microglia. In this study, we leveraged multiple RNA analysis platforms to identify new monocyte markers relevant to GBM patient outcome.High-confidence lists of mouse resident microglia- and bone marrow-derived macrophage-specific transcripts were generated using converging RNA-seq and microarray technologies and validated using qRT-PCR and flow cytometry. Expression of select cell surface markers was analyzed in brain-infiltrating macrophages and resident microglia in an induced GBM mouse model, while allogeneic bone marrow transplantation was performed to trace the origins of infiltrating and resident macrophages. Glioma tissue microarrays were examined by immunohistochemistry, and the Gene Expression Omnibus (GEO) database was queried to determine the prognostic value of identified microglia biomarkers in human GBM.We generated a unique catalog of differentially-expressed bone marrow-derived monocyte and resident microglia transcripts, and demonstrated that brain-infiltrating macrophages acquire F11R expression in GBM and following bone-marrow transplantation. Moreover, mononuclear cell F11R expression positively correlates with human high-grade glioma and additionally serves as a biomarker for GBM patient survival, regardless of GBM molecular subtype.These studies establish F11R as a novel monocyte prognostic marker for GBM critical for defining a subpopulation of stromal cells for future potential therapeutic intervention.

Gene fusions created by somatic genomic rearrangements are known to play an important role in the onset and development of some cancers, such as lymphomas and sarcomas. RNA-Seq (whole transcriptome shotgun sequencing) is proving to be a useful tool for the discovery of novel gene fusions in cancer transcriptomes. However, algorithmic methods for the discovery of gene fusions using RNA-Seq data remain underdeveloped. We have developed deFuse, a novel computational method for fusion discovery in tumor RNA-Seq data. Unlike existing methods that use only unique best-hit alignments and consider only fusion boundaries at the ends of known exons, deFuse considers all alignments and all possible locations for fusion boundaries. As a result, deFuse is able to identify fusion sequences with demonstrably better sensitivity than previous approaches. To increase the specificity of our approach, we curated a list of 60 true positive and 61 true negative fusion sequences (as confirmed by RT-PCR), and have trained an adaboost classifier on 11 novel features of the sequence data. The resulting classifier has an estimated value of 0.91 for the area under the ROC curve. We have used deFuse to discover gene fusions in 40 ovarian tumor samples, one ovarian cancer cell line, and three sarcoma samples. We report herein the first gene fusions discovered in ovarian cancer. We conclude that gene fusions are not infrequent events in ovarian cancer and that these events have the potential to substantially alter the expression patterns of the genes involved; gene fusions should therefore be considered in efforts to comprehensively characterize the mutational profiles of ovarian cancer transcriptomes.