Human endogenous retroviruses (hERVs) are remnants of exogenous retroviruses that have integrated into the genome throughout evolution. We developed a computational workflow, hervQuant, which identified more than 3,000 transcriptionally active hERVs within The Cancer Genome Atlas (TCGA) pan-cancer RNA-Seq database. hERV expression was associated with clinical prognosis in several tumor types, most significantly clear cell renal cell carcinoma (ccRCC). We explored two mechanisms by which hERV expression may influence the tumor immune microenvironment in ccRCC: (i) RIG-I-like signaling and (ii) retroviral antigen activation of adaptive immunity. We demonstrated the ability of hERV signatures associated with these immune mechanisms to predict patient survival in ccRCC, independent of clinical staging and molecular subtyping. We identified potential tumor-specific hERV epitopes with evidence of translational activity through the use of a ccRCC ribosome profiling (Ribo-Seq) dataset, validated their ability to bind HLA in vitro, and identified the presence of MHC tetramer-positive T cells against predicted epitopes. hERV sequences identified through this screening approach were significantly more highly expressed in ccRCC tumors responsive to treatment with programmed death receptor 1 (PD-1) inhibition. hervQuant provides insights into the role of hERVs within the tumor immune microenvironment, as well as evidence that hERV expression could serve as a biomarker for patient prognosis and response to immunotherapy.

No abstract

Although a subset of clear cell renal cell carcinoma (ccRCC) patients respond to immune checkpoint blockade (ICB), predictors of response remain uncertain. We investigated whether abnormal expression of endogenous retroviruses (ERVs) in tumors is associated with local immune checkpoint activation (ICA) and response to ICB. Twenty potentially immunogenic ERVs (πERVs) were identified in ccRCC in The Cancer Genome Atlas data set, and tumors were stratified into 3 groups based on their expression levels. πERV-high ccRCC tumors showed increased immune infiltration, checkpoint pathway upregulation, and higher CD8+ T cell fraction in infiltrating leukocytes compared with πERV-low ccRCC tumors. Similar results were observed in ER+/HER2- breast, colon, and head and neck squamous cell cancers. ERV expression correlated with expression of genes associated with histone methylation and chromatin regulation, and πERV-high ccRCC was enriched in BAP1 mutant tumors. ERV3-2 expression correlated with ICA in 11 solid cancers, including the 4 named above. In a small retrospective cohort of 24 metastatic ccRCC patients treated with single-agent PD-1/PD-L1 blockade, ERV3-2 expression in tumors was significantly higher in responders compared with nonresponders. Thus, abnormal expression of πERVs is associated with ICA in several solid cancers, including ccRCC, and ERV3-2 expression is associated with response to ICB in ccRCC.

Response to immune checkpoint therapy can be associated with a high mutation burden, but other mechanisms are also likely to be important. We identified a patient with metastatic gastric cancer with meaningful clinical benefit from treatment with the anti-programmed death-ligand 1 (PD-L1) antibody avelumab. This tumor showed no evidence of high mutation burden or mismatch repair defect but was strongly positive for presence of Epstein-Barr virus (EBV) encoded RNA. Analysis of The Cancer Genome Atlas gastric cancer data (25 EBV+, 80 microsatellite-instable [MSI], 310 microsatellite-stable [MSS]) showed that EBV-positive tumors were MSS. Two-sided Wilcoxon rank-sum tests showed that: 1) EBV-positive tumors had low mutation burden (median = 2.07 vs 3.13 in log10 scale, P < 10-12) but stronger evidence of immune infiltration (median ImmuneScore 2212 vs 1295, P < 10-4; log2 fold-change of CD8A = 1.85, P < 10-6) compared with MSI tumors, and 2) EBV-positive tumors had higher expression of immune checkpoint pathway (PD-1, CTLA-4 pathway) genes in RNA-seq data (log2 fold-changes: PD-1 = 1.85, PD-L1 = 1.93, PD-L2 = 1.50, CTLA-4 = 1.31, CD80 = 0.89, CD86 = 1.31, P < 10-4 each), and higher lymphocytic infiltration by histology (median tumor-infiltrating lymphocyte score = 3 vs 2, P < .001) compared with MSS tumors. These data suggest that EBV-positive low-mutation burden gastric cancers are a subset of MSS gastric cancers that may respond to immune checkpoint therapy.

Purpose An association between mutational burden and response to immune checkpoint therapy has been documented in several cancer types. The potential for such a mutational burden threshold to predict response to immune checkpoint therapy was evaluated in several clinical datasets, where mutational burden was measured either by whole-exome sequencing (WXS) or using commercially available sequencing panels. Methods WXS and RNA-seq data of 33 solid cancer types from TCGA were analyzed to determine whether a robust immune checkpoint activating mutation (iCAM) burden threshold associated with evidence of immune checkpoint activation exists in these cancers that may serve as a biomarker for response to immune checkpoint blockade therapy. Results We find that a robust iCAM threshold, associated with signatures of immune checkpoint activation, exists in 8 of 33 solid cancers: melanoma, lung adenocarcinoma, colon adenocarcinoma, endometrial cancer, stomach adenocarcinoma, cervical cancer, ER+HER2− breast cancer, and bladder-urothelial cancer. Tumors with mutational burden higher than the threshold (iCAM+) also had clear histologic evidence of lymphocytic infiltration. In published datasets of melanoma, lung adenocarcinoma and colon cancer, patients with iCAM+ tumors had significantly better response to immune checkpoint therapy compared to those with iCAM− tumors. ROC analysis using TCGA predictions as gold standard showed that iCAM+ tumors are accurately identifiable using clinical sequencing assays, such as FoundationOne or StrandAdvantage. Using the FoundationOne derived threshold, analysis of 113 melanoma tumors, showed that iCAM+ patients have significantly better response to immune checkpoint therapy. iCAM+ and iCAM− tumors have distinct mutation patterns and different immune microenvironments. Conclusion In 8 solid cancers, a mutational burden threshold exists that may predict response to immune checkpoint blockade. This threshold is identifiable using available clinical sequencing assays.