Deciphering whether actionable driver mutations are found in all or a subset of tumor cells will likely be required to improve drug development and precision medicine strategies. We analyzed nine ...cancer types to determine the subclonal frequencies of driver events, to time mutational processes during cancer evolution, and to identify drivers of subclonal expansions. Although mutations in known driver genes typically occurred early in cancer evolution, we also identified later subclonal "actionable" mutations, including BRAF (V600E), IDH1 (R132H), PIK3CA (E545K), EGFR (L858R), and KRAS (G12D), which may compromise the efficacy of targeted therapy approaches. More than 20% of IDH1 mutations in glioblastomas, and 15% of mutations in genes in the PI3K (phosphatidylinositol 3-kinase)-AKT-mTOR (mammalian target of rapamycin) signaling axis across all tumor types were subclonal. Mutations in the RAS-MEK (mitogen-activated protein kinase kinase) signaling axis were less likely to be subclonal than mutations in genes associated with PI3K-AKT-mTOR signaling. Analysis of late mutations revealed a link between APOBEC-mediated mutagenesis and the acquisition of subclonal driver mutations and uncovered putative cancer genes involved in subclonal expansions, including CTNNA2 and ATXN1. Our results provide a pan-cancer census of driver events within the context of intratumor heterogeneity and reveal patterns of tumor evolution across cancers. The frequent presence of subclonal driver mutations suggests the need to stratify targeted therapy response according to the proportion of tumor cells in which the driver is identified.
Validating prognostic or predictive candidate genes in appropriately powered breast cancer cohorts are of utmost interest. Our aim was to develop an online tool to draw survival plots, which can be ...used to assess the relevance of the expression levels of various genes on the clinical outcome both in untreated and treated breast cancer patients. A background database was established using gene expression data and survival information of 1,809 patients downloaded from GEO (Affymetrix HGU133A and HGU133+2 microarrays). The median relapse free survival is 6.43 years, 968/1,231 patients are estrogen-receptor (ER) positive, and 190/1,369 are lymph-node positive. After quality control and normalization only probes present on both Affymetrix platforms were retained (
n
= 22,277). In order to analyze the prognostic value of a particular gene, the cohorts are divided into two groups according to the median (or upper/lower quartile) expression of the gene. The two groups can be compared in terms of relapse free survival, overall survival, and distant metastasis free survival. A survival curve is displayed, and the hazard ratio with 95% confidence intervals and logrank
P
value are calculated and displayed. Additionally, three subgroups of patients can be assessed: systematically untreated patients, endocrine-treated ER positive patients, and patients with a distribution of clinical characteristics representative of those seen in general clinical practice in the US. Web address:
www.kmplot.com
. We used this integrative data analysis tool to confirm the prognostic power of the proliferation-related genes
TOP2A
and
TOP2B
,
MKI67
,
CCND2
,
CCND3
,
CCNDE2
, as well as
CDKN1A
, and
TK2
. We also validated the capability of microarrays to determine estrogen receptor status in 1,231 patients. The tool is highly valuable for the preliminary assessment of biomarkers, especially for research groups with limited bioinformatic resources.
Interpretation of gene expression microarrays requires a mapping from probe set to gene. On many Affymetrix gene expression microarrays, a given gene may be detected by multiple probe sets, which may ...deliver inconsistent or even contradictory measurements. Therefore, obtaining an unambiguous expression estimate of a pre-specified gene can be a nontrivial but essential task.
We developed scoring methods to assess each probe set for specificity, splice isoform coverage, and robustness against transcript degradation. We used these scores to select a single representative probe set for each gene, thus creating a simple one-to-one mapping between gene and probe set. To test this method, we evaluated concordance between protein measurements and gene expression values, and between sets of genes whose expression is known to be correlated. For both test cases, we identified genes that were nominally detected by multiple probe sets, and we found that the probe set chosen by our method showed stronger concordance.
This method provides a simple, unambiguous mapping to allow assessment of the expression levels of specific genes of interest.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Personalization of immunotherapies such as cancer vaccines and adoptive T cell therapy depends on identification of patient-specific neo-epitopes that can be specifically targeted. MuPeXI, the mutant ...peptide extractor and informer, is a program to identify tumor-specific peptides and assess their potential to be neo-epitopes. The program input is a file with somatic mutation calls, a list of HLA types, and optionally a gene expression profile. The output is a table with all tumor-specific peptides derived from nucleotide substitutions, insertions, and deletions, along with comprehensive annotation, including HLA binding and similarity to normal peptides. The peptides are sorted according to a priority score which is intended to roughly predict immunogenicity. We applied MuPeXI to three tumors for which predicted MHC-binding peptides had been screened for T cell reactivity, and found that MuPeXI was able to prioritize immunogenic peptides with an area under the curve of 0.63. Compared to other available tools, MuPeXI provides more information and is easier to use. MuPeXI is available as stand-alone software and as a web server at
http://www.cbs.dtu.dk/services/MuPeXI
.
Phages and cancer immunity
Gut bacteria are involved in the education of T cell immune responses, and the intestinal ecosystem influences anticancer immunity. Fluckiger
et al.
report microbial ...antigens that might cross-react with antigens associated with tumor cells. They found that a type of intestinal bacteria called enterococci harbor a bacteriophage that modulates immune responses. In mouse models, administration of enterococci containing the bacteriophage boosted T cell responses after treatment with chemotherapy or programmed cell death protein 1 (PD-1) blockade. In humans, the presence of the bacteriophage was associated with improved survival after PD-1 immunotherapy. A fraction of human T cells specific for naturally processed melanoma epitopes appeared to be able to recognize microbial peptides. This “molecular mimicry” may represent cross-reactivity between tumors and microbial antigens.
Science
, this issue p.
936
A bacteriophage found in commensal bacteria can modulate the effects of cancer immunotherapy.
Intestinal microbiota have been proposed to induce commensal-specific memory T cells that cross-react with tumor-associated antigens. We identified major histocompatibility complex (MHC) class I–binding epitopes in the tail length tape measure protein (TMP) of a prophage found in the genome of the bacteriophage
Enterococcus hirae
. Mice bearing
E. hirae
harboring this prophage mounted a TMP-specific H-2K
b
–restricted CD8
+
T lymphocyte response upon immunotherapy with cyclophosphamide or anti–PD-1 antibodies. Administration of bacterial strains engineered to express the TMP epitope improved immunotherapy in mice. In renal and lung cancer patients, the presence of the enterococcal prophage in stools and expression of a TMP–cross-reactive antigen by tumors correlated with long-term benefit of PD-1 blockade therapy. In melanoma patients, T cell clones recognizing naturally processed cancer antigens that are cross-reactive with microbial peptides were detected.
Human endogenous retroviruses (HERV) form a substantial part of the human genome, but mostly remain transcriptionally silent under strict epigenetic regulation, yet can potentially be reactivated by ...malignant transformation or epigenetic therapies. Here, we evaluate the potential for T cell recognition of HERV elements in myeloid malignancies by mapping transcribed HERV genes and generating a library of 1169 potential antigenic HERV-derived peptides predicted for presentation by 4 HLA class I molecules. Using DNA barcode-labeled MHC-I multimers, we find CD8
T cell populations recognizing 29 HERV-derived peptides representing 18 different HERV loci, of which HERVH-5, HERVW-1, and HERVE-3 have more profound responses; such HERV-specific T cells are present in 17 of the 34 patients, but less frequently in healthy donors. Transcriptomic analyses reveal enhanced transcription of the HERVs in patients; meanwhile DNA-demethylating therapy causes a small and heterogeneous enhancement in HERV transcription without altering T cell recognition. Our study thus uncovers T cell recognition of HERVs in myeloid malignancies, thereby implicating HERVs as potential targets for immunotherapeutic therapies.
Clear cell renal cell carcinoma (ccRCC) is characterized by near-universal loss of the short arm of chromosome 3, deleting several tumor suppressor genes. We analyzed whole genomes from 95 biopsies ...across 33 patients with clear cell renal cell carcinoma. We find hotspots of point mutations in the 5′ UTR of TERT, targeting a MYC-MAX-MAD1 repressor associated with telomere lengthening. The most common structural abnormality generates simultaneous 3p loss and 5q gain (36% patients), typically through chromothripsis. This event occurs in childhood or adolescence, generally as the initiating event that precedes emergence of the tumor’s most recent common ancestor by years to decades. Similar genomic changes drive inherited ccRCC. Modeling differences in age incidence between inherited and sporadic cancers suggests that the number of cells with 3p loss capable of initiating sporadic tumors is no more than a few hundred. Early development of ccRCC follows well-defined evolutionary trajectories, offering opportunity for early intervention.
Display omitted
•Novel hotspot of driver mutations in 5′-UTR repressor of TERT, expanding telomeres•Most common cause of 3p loss is a chromothripsis event, generating concurrent 5q gain•t(3;5) event occurs in childhood or adolescence, decades before tumor diagnosed•Initial clonal expansion after 3p loss starts from only a few hundred cells
Combination of whole-genome sequencing analysis and a multi-region sampling approach provides insights into the nature and timing of key oncogenic events in clear cell renal cell carcinoma, depicts the evolutionary trajectories of tumors in patients and highlights the opportunity for early intervention.
Spatial and temporal dissection of the genomic changes occurring during the evolution of human non–small cell lung cancer (NSCLC) may help elucidate the basis for its dismal prognosis. We sequenced ...25 spatially distinct regions from seven operable NSCLCs and found evidence of branched evolution, with driver mutations arising before and after subclonal diversification. There was pronounced intratumor heterogeneity in copy number alterations, translocations, and mutations associated with APOBEC cytidine deaminase activity. Despite maintained carcinogen exposure, tumors from smokers showed a relative decrease in smoking-related mutations over time, accompanied by an increase in APOBEC-associated mutations. In tumors from former smokers, genome-doubling occurred within a smoking-signature context before subclonal diversification, which suggested that a long period of tumor latency had preceded clinical detection. The regionally separated driver mutations, coupled with the relentless and heterogeneous nature of the genome instability processes, are likely to confound treatment success in NSCLC.
Defects in BRCA1, BRCA2 and other genes of the homology-dependent DNA repair (HR) pathway cause an elevated rate of mutagenesis, eliciting specific mutation patterns including COSMIC signature SBS3. ...Using genome sequencing of knock-out cell lines we show that Y family translesion synthesis (TLS) polymerases contribute to the spontaneous generation of base substitution and short insertion/deletion mutations in BRCA1 deficient cells, and that TLS on DNA adducts is increased in BRCA1 and BRCA2 mutants. The inactivation of 53BP1 in BRCA1 mutant cells markedly reduces TLS-specific mutagenesis, and rescues the deficiency of template switch-mediated gene conversions in the immunoglobulin V locus of BRCA1 mutant chicken DT40 cells. 53BP1 also promotes TLS in human cellular extracts in vitro. Our results show that HR deficiency-specific mutagenesis is largely caused by TLS, and suggest a function for 53BP1 in regulating the choice between TLS and error-free template switching in replicative DNA damage bypass.