<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Immunol.</journal-id>
<journal-title>Frontiers in Immunology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Immunol.</abbrev-journal-title>
<issn pub-type="epub">1664-3224</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fimmu.2022.1039631</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Immunology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Epigenetic-related gene mutations serve as potential biomarkers for immune checkpoint inhibitors in microsatellite-stable colorectal cancer</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Liu</surname>
<given-names>Chao</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<xref ref-type="author-notes" rid="fn003">
<sup>&#x2020;</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2095276"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Xiao</surname>
<given-names>Huiting</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<xref ref-type="author-notes" rid="fn003">
<sup>&#x2020;</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1327811"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Cui</surname>
<given-names>Luying</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<xref ref-type="author-notes" rid="fn003">
<sup>&#x2020;</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2095114"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Fang</surname>
<given-names>Lin</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/788029"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Han</surname>
<given-names>Shuling</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2095146"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Ruan</surname>
<given-names>Yuli</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2095129"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Zhao</surname>
<given-names>Wenyuan</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/740065"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Zhang</surname>
<given-names>Yanqiao</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<xref ref-type="aff" rid="aff5">
<sup>5</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/789381"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Department of Gastrointestinal Medical Oncology, Harbin Medical University Cancer Hospital</institution>, <addr-line>Harbin</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Clinical Research Center for Colorectal Cancer in Heilongjiang, Harbin Medical University Cancer Hospital</institution>, <addr-line>Harbin</addr-line>, <country>China</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Key Laboratory of Tumor Immunology in Heilongjiang, Harbin Medical University Cancer Hospital</institution>, <addr-line>Harbin</addr-line>, <country>China</country>
</aff>
<aff id="aff4">
<sup>4</sup>
<institution>College of Bioinformatics Science and Technology, Harbin Medical University</institution>, <addr-line>Harbin</addr-line>, <country>China</country>
</aff>
<aff id="aff5">
<sup>5</sup>
<institution>Translational Medicine Research and Cooperation Center of Northern China, Heilongjiang Academy of Medical Sciences</institution>, <addr-line>Harbin</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Jia Wei, Nanjing Drum Tower Hospital, China</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Muzamil Yaqub Want, Roswell Park Comprehensive Cancer Center, United States; Koshy Nithin Thomas, Sanjay Gandhi Post Graduate Institute of Medical Sciences (SGPGI), India</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Yanqiao Zhang, <email xlink:href="mailto:yanqiaozhang@ems.hrbmu.edu.cn">yanqiaozhang@ems.hrbmu.edu.cn</email>; Wenyuan Zhao, <email xlink:href="mailto:zhaowenyuan@ems.hrbmu.edu.cn">zhaowenyuan@ems.hrbmu.edu.cn</email>
</p>
</fn>
<fn fn-type="equal" id="fn003">
<p>&#x2020;These authors have contributed equally to this work and share first authorship</p>
</fn>
<fn fn-type="other" id="fn002">
<p>This article was submitted to Cancer Immunity and Immunotherapy, a section of the journal Frontiers in Immunology</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>21</day>
<month>11</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>13</volume>
<elocation-id>1039631</elocation-id>
<history>
<date date-type="received">
<day>08</day>
<month>09</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>31</day>
<month>10</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2022 Liu, Xiao, Cui, Fang, Han, Ruan, Zhao and Zhang</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Liu, Xiao, Cui, Fang, Han, Ruan, Zhao and Zhang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<sec>
<title>Background</title>
<p>Combination therapy with immune checkpoint inhibitors (ICIs) may benefit approximately 10-20% of microsatellite-stable colorectal cancer (MSS-CRC) patients. However, there is a lack of optimal biomarkers. This study aims to understand the predictive value of epigenetic-related gene mutations in ICIs therapy in MSS-CRC patients.</p>
</sec>
<sec>
<title>Methods</title>
<p>We analyzed DNA sequences and gene expression profiles from The Cancer Genome Atlas (TCGA) to examine their immunological features. The Harbin Medical University Cancer Hospital (HMUCH) clinical cohort of MSS-CRC patients was used to validate the efficacy of ICIs in patients with epigenetic-related gene mutations (Epigenetic_Mut).</p>
</sec>
<sec>
<title>Results</title>
<p>In TCGA, 18.35% of MSS-CRC patients (78/425) had epigenetic-related gene mutations. The Epigenetic_Mut group had a higher tumor mutation burden (TMB) and frameshift mutation (FS_mut) rates. In all MSS-CRC samples, Epigenetic_Mut was elevated in the immune subtype (CMS1) and had a strong correlation with immunological features. Epigenetic_Mut was also associated with favorable clinical outcomes in MSS-CRC patients receiving anti-PD-1-based therapy from the HMUCH cohort. Using immunohistochemistry and flow cytometry, we demonstrated that Epigenetic_Mut samples were associated with increased anti-tumor immune cells both in tumor tissues and peripheral blood.</p>
</sec>
<sec>
<title>Conclusion</title>
<p>MSS-CRC patients with epigenetic regulation impairment exhibit an immunologically active environment and may be more susceptible to treatment strategies based on ICIs.</p>
</sec>
</abstract>
<kwd-group>
<kwd>immune checkpoint inhibitor therapy</kwd>
<kwd>biomarker</kwd>
<kwd>microsatellite-stable</kwd>
<kwd>colorectal cancer</kwd>
<kwd>epigenetic-related gene mutations</kwd>
</kwd-group>
<counts>
<fig-count count="6"/>
<table-count count="1"/>
<equation-count count="0"/>
<ref-count count="43"/>
<page-count count="12"/>
<word-count count="5068"/>
</counts>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<title>Introduction</title>
<p>Immune checkpoint inhibitor (ICI) therapy has achieved impressive success in deficient mismatch repair (dMMR)/microsatellite instability-high (MSI-H) colorectal cancer (CRC). ICI has been considered a standard therapy by the FDA, including the use of programmed death receptor 1 (PD-1) monoclonal antibodies and CTL-associated protein 4 (CTLA-4) monoclonal antibodies (<xref ref-type="bibr" rid="B1">1</xref>&#x2013;<xref ref-type="bibr" rid="B3">3</xref>). However, the vast majority of CRC cases (approximately 85%) are characterized by proficient mismatch repair (pMMR)/microsatellite stability (MSS) tumors that do not respond to ICIs (<xref ref-type="bibr" rid="B4">4</xref>). Recent studies suggest that a subgroup (approximately 10-20%) of MSS-CRC patients might benefit from combination regimens of ICIs (<xref ref-type="bibr" rid="B5">5</xref>&#x2013;<xref ref-type="bibr" rid="B7">7</xref>). Therefore, predictive biomarkers for screening these patients are urgently needed.</p>
<p>Current clinical and investigational studies of screening MSS-CRC patients who would benefit from ICIs treatment are limited. PD-L1 expression is a classic biomarker, but the Keynote-028 study demonstrated that PD-L1<sup>+</sup> MSS-CRC patients could not benefit from ICI therapy (<xref ref-type="bibr" rid="B8">8</xref>). <italic>POLD1</italic>/<italic>POLE</italic> mutations are predictive but occur in only 1% of MSS-CRC patients (<xref ref-type="bibr" rid="B9">9</xref>). Biomarkers such as tumor mutation burden (TMB), tumor-infiltrating lymphocytes (TIL), neo-antigen load (NAL), and immune-regulatory gene expression profiling (iGEP) may allow the selection of clinical patients for ICIs. However, the lack of uniform detection methods and validated cutoffs limit the use of these methods (<xref ref-type="bibr" rid="B5">5</xref>, <xref ref-type="bibr" rid="B10">10</xref>&#x2013;<xref ref-type="bibr" rid="B12">12</xref>). Several emerging biomarkers, such as gut microbiota and T-cell-receptor (TCR) sequencing, have also shown predictive value, although they are not yet clinically applicable (<xref ref-type="bibr" rid="B7">7</xref>, <xref ref-type="bibr" rid="B13">13</xref>). DNA damage response (DDR) gene mutations may induce a hypermutational phenotype (<xref ref-type="bibr" rid="B14">14</xref>), and recent studies have shown that patients with MSS-CRC and mutations in the DDR system have better immune responses and outcomes following ICI therapy (<xref ref-type="bibr" rid="B15">15</xref>, <xref ref-type="bibr" rid="B16">16</xref>). However, the pathogenicity of different DDR gene mutations in MSS-CRC remains unclear, and their incidence is significantly lower than in endometrial, ovarian, or biliary tract cancers (<xref ref-type="bibr" rid="B17">17</xref>).</p>
<p>Epigenomic alterations can affect tumor immunogenicity and anti-tumor responses by regulating genome stability and chromatin accessibility (<xref ref-type="bibr" rid="B18">18</xref>). Additionally, several epigenetic-related gene mutations have been shown to exhibit predictive functions in ICI therapy for multiple types of tumors. ARID1A, an AT-rich interactive domain-containing protein 1A, is a component of the switching defective/sucrose non-fermenting (SWI/SNF) complex that plays a role in chromatin remodeling (<xref ref-type="bibr" rid="B19">19</xref>), and increasing evidence suggests that ARID1A alterations are correlated with better outcomes after ICI therapy for bladder cancer, nonsmall-cell lung cancer (NSCLC), and gastric cancer (<xref ref-type="bibr" rid="B20">20</xref>, <xref ref-type="bibr" rid="B21">21</xref>). ARID1A mutation is defined as an immunologically active subgroup in MSS-CRC patients with abundant intra-tumoral T-cell infiltration (<xref ref-type="bibr" rid="B22">22</xref>). Lysine methyltransferase 2 (KMT2) family members facilitate transcription and gene accessibility by methylating lysine 4 on histone H3 (H3k4) (<xref ref-type="bibr" rid="B23">23</xref>), and <italic>KMT2</italic> family mutations have also been linked to a favorable response to ICIs in multiple cancers (<xref ref-type="bibr" rid="B24">24</xref>). Furthermore, as identified using clustered regularly interspaced short palindromic repeats (CRISPR), KMT2D mutant tumors exhibit an increased mutation burden, IFN-&#x3b3;-stimulated antigen presentation, and a higher sensitivity to ICIs. Moreover, disruption of DNA methylation signatures has been identified as a marker of anti-PD-1 therapy efficacy in NSCLC (<xref ref-type="bibr" rid="B25">25</xref>), and TET1, a DNA demethylase, enhances the immunotherapeutic effect (<xref ref-type="bibr" rid="B26">26</xref>). Although this evidence points to the role of epigenetic regulation in anti-tumor immune responses, there is no clinical data on the association between comprehensive epigenetic-related gene mutations (mutations in genes that are involved in epigenetic modifications) and the clinical benefit of ICIs in MSS-CRC.</p>
<p>Given the proposed role of epigenetic regulation impairment in predicting the response to ICIs, we hypothesize that epigenetic-related gene mutations in MSS-CRC cause hypermutation and improve the expression of immune response gene sets. As a result, we conducted this study to clarify the value of epigenetic-related gene mutations as an indicator of immunotherapy efficacy in patients with MSS-CRC. For this purpose, we analyzed whole-exome sequencing (WES) data from TCGA to study TMB, frameshift mutation (FS-mutation), and immune characteristics of Epigenetic_Mut and Epigenetic_Wt groups of MSS-CRC samples. Additionally, in a Chinese clinical MSS-CRC cohort of 89 patients who received PD-1-based treatment, we found that Epigenetic_Mut was associated with favorable clinical outcomes. Here, we report the relationships between epigenetic-related gene mutations and TMB, FS-mutation, immunomodulatory mRNA expression signature, and ICI therapy efficacy in patients with MSS-CRC.</p>
</sec>
<sec id="s2" sec-type="materials|methods">
<title>Materials and methods</title>
<sec id="s2_1">
<title>Patient information and sample collection</title>
<p>To determines the incidence of epigenetic-related gene mutations in MSS-CRC, we analyzed DNA sequencing and gene expression profiles of 514 MSS-CRC patients from two cohorts: (1) a TCGA cohort consisting of 425 MSS-CRC patients and (2) a HMUCH cohort comprising 89 Chinese patients with annotated response and mutational data from Harbin Medical University Cancer Hospital (the inclusion and exclusion criteria are shown in <xref ref-type="supplementary-material" rid="SF1">
<bold>Supplementary Figure&#xa0;1</bold>
</xref>). This study was approved by the Ethics Committee of the Harbin Medical University Cancer Hospital (No. KY2022-20).</p>
</sec>
<sec id="s2_2">
<title>Epigenetic-related gene status definition</title>
<p>Epigenetic-related gene status (Epigenetic_Wt or Epigenetic_Mut) was defined based on the presence of a loss-of-function (LOF) variant in 68 genes that have been proposed as core genes of epigenetic regulation (<xref ref-type="bibr" rid="B18">18</xref>). <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;1</bold>
</xref> presents a detailed description. Nonsense, frameshift, and splice site changes within consensus regions and start lost/gained variants were considered to be LOF variants. Missense and in-frame variants were excluded from the analysis.</p>
</sec>
<sec id="s2_3">
<title>DNA extraction and sequencing</title>
<p>For the TCGA cohort, gene mutation data were acquired using the GDC Data Portal. We assessed the mutational status of epigenetic-related genes in CRC using exome-sequencing data from HMUCH. For analysis, DNA was extracted using a DNA Kit (Applied Biosystems, Foster City, CA, USA), from whole blood samples or formalin-fixed paraffin-embedded (FFPE) tissues of each patient. The lymphocytes from the whole blood samples were isolated by centrifugation at 1,600 &#xd7; g for 10&#xa0;min in red cell lysis buffer (Tiangen, RT122, Beijing, China) at 25&#xb0;C, and DNA was extracted using a genomic DNA kit (Tiangen, DP304, Beijing, China). We sheared the DNA into fragments of 150-200 bp using an ultrasonicator and used a KAPA Kit (KAPA Biosystems, Wilmington, MA, USA) to prepare DNA fragment libraries for the Illumina platform (Illumina HiSeq X-Ten, Illumina, USA). Probe hybridization capture technology and Illumina high-throughput sequencing were used to detect the exonic regions and some intronic regions of 825 tumor-related genes (Genetron Health Co., Ltd. Beijing, China) (<xref ref-type="supplementary-material" rid="SM2">
<bold>Supplementary Table&#xa0;2</bold>
</xref>).</p>
</sec>
<sec id="s2_4">
<title>Analysis of MSI status, TMB, and FS-mutation in the TCGA and HMUCH cohorts</title>
<p>MSI status for the TCGA cohort was determined using the MSI sensor (version 0.5). In brief, for MSI sensor scores &lt; 3.5, samples were considered to be MSS; otherwise, they were considered MSI (<xref ref-type="bibr" rid="B27">27</xref>). Published studies using the TCGA cohort provided FS-mutation and TMB data (<xref ref-type="bibr" rid="B28">28</xref>&#x2013;<xref ref-type="bibr" rid="B30">30</xref>), and MSI status for the HMUCH cohort was determined using a 3730 sequencer (Life Technologies, Carlsbad, CA, USA). For this purpose, whole blood samples or prepared FFPE tissue were diluted to 2 ng/&#x3bc;L or 20 ng/&#x3bc;L, respectively, followed by the addition of 2.8 &#x3bc;L of ddH<sub>2</sub>O, 4 &#x3bc;L of 2.5&#xd7; Buffer A, 2 &#x3bc;L of 5&#xd7; MSI Primer Mix, and 0.2 &#x3bc;L of Taq DNA Polymerase I. PCR amplification was carried out as follows: pre-denaturation at 95&#xb0;C for 5&#xa0;min; followed by 30 cycles at 94&#xb0;C for 30 s, 60&#xb0;C for 1&#xa0;min, and 70&#xb0;C for 1&#xa0;min; and then a final extension at 60&#xb0;C for 30&#xa0;min. Finally, the temperature was reduced to 15&#xb0;C, and the samples were centrifuged at 3,000 &#xd7; g for 1&#xa0;min. NR-21 and BAT-26 were labeled with blue fluorescent dye, BAT-25 with green dye, and NR-24 and MONO-27 with yellow dye. Finally, tumors were classified as MSI-H if two or more markers showed instability; otherwise they were classified as MSS.</p>
</sec>
<sec id="s2_5">
<title>Analysis of the consensus molecular subtypes (CMSs) in the TCGA cohort</title>
<p>Consensus molecular subtypes (CMSs) are classification systems for CRC and include immune (CMS1), canonical (CMS2), metabolic (CMS3), and mesenchymal (CMS4) subtypes. These subtypes were identified through a large-scale analytical study and have unique molecular and metabolic characteristics (<xref ref-type="bibr" rid="B31">31</xref>).</p>
</sec>
<sec id="s2_6">
<title>Immune-related signature analysis</title>
<p>Our study compared the RNA expression of patients with Epigenetic_Mut and Epigenetic_Wt using gene signatures for the IFN-&#x3b3; pathway and other immunological responses (<xref ref-type="supplementary-material" rid="SM3">
<bold>Supplementary Table&#xa0;3</bold>
</xref>) (<xref ref-type="bibr" rid="B12">12</xref>, <xref ref-type="bibr" rid="B32">32</xref>). We obtained TCGA transcriptome profiles from the GDC data portal, and used transcripts per kilobase million (TPM) normalization to normalize gene expression. The geometric mean of gene expression levels in the log<sub>2</sub> (TPM + 1) format was used to evaluate immune signatures.</p>
</sec>
<sec id="s2_7">
<title>Clinical outcomes</title>
<p>The objective response rate (ORR), disease control rate (DCR), progression-free survival (PFS), and overall survival (OS) were the main clinical outcomes of interest. The Response Evaluation Criteria in Solid Tumors (RECIST) version 1.1 was used for the assessment of ORR and divided into complete response (CR) and partial response (PR). DCR was defined as CR, PR, or stable disease (SD) lasting more than six months. PFS was evaluated from when immunological therapy was initiated until progression or death, and patients who did not progress were examined at the last scan. OS was evaluated from the start of ICI therapy until patient death or the end of the trial, and the patients with whom we lost contact were classified based on the date of last contact.</p>
</sec>
<sec id="s2_8">
<title>Immunohistochemistry (IHC)</title>
<p>Primary tumor paraffin sections of 4 &#x3bc;m were processed for immunochemistry to evaluate CD8<sup>+</sup> and FOXP3<sup>+</sup> lymphocytes according to the following protocol: roast, deparaffination, and rehydration before performing heat-mediated antigen retrieval with EDTA buffer (pH 9.0), inactivation of endogenous peroxidase activity with 3% H<sub>2</sub>O<sub>2</sub>, incubation with antibody against CD8 (ab101500, 1:500; Abcam, Cambridge, UK) or against FOXP3 (ab200334, 1:500; Abcam) at 4&#xb0;C overnight, exposure to a DAB IHC Detection Kit after incubation with biotinylated secondary antibodies, and counterstaining with Mayer&#x2019;s hematoxylin solution. An open-source platform for biological-image analysis (Fiji/ImageJ) was used to estimate the densities of CD8<sup>+</sup> and FOXP3<sup>+</sup> lymphocytes.</p>
</sec>
<sec id="s2_9">
<title>Flow cytometry analysis</title>
<p>The peripheral blood mononuclear cells (PBMC) of CRC patients were isolated by centrifugation with erythrocytes lysate and were used to analyze PD1<sup>+</sup>CD8<sup>+</sup>T cells and CD3<sup>-</sup>CD56<sup>+</sup>CD16<sup>+</sup>NK cells by flow cytometry. The PBMC were stained for 30min on ice using the following antibodies: FITC anti-human CD8 (344704, Biolegend), PE anti-human PD1(367404, Biolegend), APC anti-human CD3 (300312, Biolegend), PE anti-human CD56 (985902, Biolegend), and PerCP anti-human CD16 (302030, Biolegend). Stained cell suspensions were analyzed using the BD flow cytometer (BD Accuri C6 Plus), and data analysis was performed using FlowJo_v10.8.1.</p>
</sec>
<sec id="s2_10">
<title>Statistical analysis</title>
<p>Fisher&#x2019;s exact test was used to analyze the relationship between epigenetic-related gene mutations and the ORR or DCR, and the Kaplan&#x2013;Meier method and log-rank test were employed to examine the PFS and OS probabilities of the Epigenetic_Mut and Epigenetic_Wt CRC groups. Based on the Mann&#x2013;Whitney U-test, TMB, FS-mutation, tumor-infiltrating lymphocytes, expression of immune-related genes, and immune signatures were compared between the Epigenetic_Mut and Epigenetic_Wt CRC groups. Statistical analysis was conducted using two-sided tests with a nominal significance level of 0.05 using R version 3.5.2.</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<title>Results</title>
<sec id="s3_1">
<title>The mutational landscape of epigenetic-related genes of MSS-CRC in the TCGA cohort</title>
<p>A total of 68 epigenetic-related genes involved in 13 different pathways were included in the current research, including genes involved in modifying DNA, histones, and protein complexes that reshape chromatin structure (<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;1</bold>
</xref>). In the TCGA cohort, MSI-H and MSS-CRCs had epigenetic-related gene mutation frequencies of 66.67% (50/75) and 18.35% (78/425), respectively. The three most frequently mutated pathways in the MSS-CRC cases from TCGA were SWI_SNF, Histone_methylase, and CHD (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1A</bold>
</xref>), and the epigenetic-related genes ARID1A, KMT2C, and RSF1 had the highest mutation rates in the TCGA cohort (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1B</bold>
</xref>).</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Mutational landscape of Epigenetic-related genes associated with MSS-CRC cases from the TCGA and HMUCH cohorts. <bold>(A)</bold> The frequency of epigenetic regulatory pathway alteration in MSS-CRC cases and MSI_H samples from the TCGA cohort. <bold>(B)</bold> The top 20 mutated epigenetic-related genes in MSS-CRC samples from the TCGA cohort.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-13-1039631-g001.tif"/>
</fig>
</sec>
<sec id="s3_2">
<title>Epigenetic-related gene mutations are linked with the TMB, FS-mutation, and molecular subtype of CRC</title>
<p>High levels of TMB and FS-mutations (FS_mut) reflect a high degree of genomic instability and potential immunogenicity of a tumor, and both of these are therefore potential biomarkers of immune checkpoint inhibitor responsiveness. Hence, we examined the relationships between TMB, FS_mut, and epigenetic-related gene mutation status. In TCGA cohort, epigenetic-related gene mutations were associated with an increased incidence of TMB in MSS-CRC (median mutation rate of 4.76/mb vs. 4.99/mb in Wt and Mut cases, respectively; <italic>p</italic> = 7.4e-05; <xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2A</bold>
</xref>). A higher rate of FS_mut was also linked with epigenetic-related gene mutations in MSS-CRC (median frameshift mutation rate of 1.39/mb vs. 1.79/mb in Wt and Mut cases, respectively, <italic>p</italic> = 3.5e-06; <xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2B</bold>
</xref>). Molecular subtypes of CRC (CMS) are currently a highly recognized classification method for CRC that can accurately guide patient treatment and prognosis. CMS1, also known as the immune subtype, has better immune activity and high reactivity to ICIs. Here, we analyzed the distribution of Epigenetic_Mut samples based on molecular subtype in all CRC and MSS-CRC cases (<xref ref-type="fig" rid="f2">
<bold>Figures&#xa0;2C, D</bold>
</xref>). Among the CMS1-CRC cases, 74.12% were Epigenetic_Mut samples (74.12%, 63/85), but in MSS CMS1-CRC cases, this rate was 40% (12/30). Both for all samples and MSS CRC specifically, Epigenetic_Mut samples were enriched in the CMS1 (immune subtype) group.</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Epigenetic-related gene mutations are linked with the TMB, FS-mutation, and molecular subtypes of CRC. <bold>(A)</bold> TMB violin plot of Epigenetic_Mut and Epigenetic_Wt from MSSCRC samples. <bold>(B)</bold> FS-mutation rate violin plot of Epigenetic_Mut and Epigenetic_Wt from MSS-CRC samples. <bold>(C)</bold> Molecular subtype-specific fold enrichment of epigenetic-related genes mutation in all CRC cases (MSI-H/MSS). <bold>(D)</bold> Molecular subtype-specific fold enrichment of ARID1A mutation in MSS-CRC.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-13-1039631-g002.tif"/>
</fig>
</sec>
<sec id="s3_3">
<title>Epigenetic_Mut is related to increased immune activity in MSS CRC</title>
<p>To identify the tumor immune microenvironment, we compared Epigenetic_Mut and Epigenetic_Wt for immune signatures, tumor-infiltrating lymphocytes, and expression of immune checkpoints and key genes. We demonstrated that epigenetic-related gene mutations increased the expression of immune response genes, including those involved in the IFN-&#x3b3; pathway, antigen presentation, and cytotoxic T-cell function (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3A</bold>
</xref>). In addition, the expression of NK cell-related genes was increased in the Epigenetic_Mut group. Other immune cells also showed an upward trend, but no statistical difference was observed due to the limited cohort size (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3B</bold>
</xref>). Finally, we compared the expression of immune checkpoints and key genes between the two groups. In line with the immune response pathway, several immune checkpoints and key genes were upregulated in the Epigenetic_Mut group. In particular, the expression of <italic>LAG3</italic> and HAVCR2 was significantly elevated, and elevated levels of <italic>TNFRSF4</italic>, <italic>PDCD1</italic>, and <italic>IL4l1</italic> were very nearly statistically significant (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3C</bold>
</xref>).</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Epigenetic_Mut is associated with increased immune activity in MSS-CRC. <bold>(A)</bold> The RNA expression of immune response gene sets in MSS-CRC based on the epigenetic-related genes&#x2019; mutational status. <bold>(B)</bold> The RNA expression of immune cells gene sets in MSS-CRC based on the epigenetic-related genes&#x2019; mutational status. <bold>(C)</bold> The RNA expression of a single immune response gene in MSS-CRC based on the epigenetic-related genes&#x2019; mutational status.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-13-1039631-g003.tif"/>
</fig>
</sec>
<sec id="s3_4">
<title>Epigenetic_Mut predicts favorable clinical outcomes following ICI therapy</title>
<p>Next, to validate the function of epigenetic-related gene mutations further in predicting responsiveness to ICI therapy in MSS-CRC, we collected a clinical cohort of 89 MSS-CRC patients who had received PD-1 mAb-based treatment. <xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref> shows the baseline patient characteristics based on epigenetic-related gene status. Of the 89 patients, 24 had Epigenetic_Mut, and 65 had Epigenetic_Wt. Using RECIST version 1.1, we evaluated the patients&#x2019; best overall responses. Compared to Epigenetic_Wt, Epigenetic_Mut had a significantly higher ORR (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4A</bold>
</xref>, 37.50% (9/24) vs. 15.38% (10/65), Fisher&#x2019;s exact test <italic>P</italic> = 0.039). As for DCR, the rate was 66.67% (16/24) in patients with epigenetic-related gene mutations from ICI treatment but only 36.92% (24/65) in patients without epigenetic-related gene mutations (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4B</bold>
</xref>, Fisher&#x2019;s exact test <italic>P</italic> = 0.017). As expected, PFS was greatly improved in patients with epigenetic-related gene mutations compared to those without epigenetic-related gene mutations in this cohort (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4C</bold>
</xref>, mPFS:6.00 vs. 3.17 months, Log_rank <italic>P</italic> = 0.002, HR = 0.4778), and ICI treatment also had a greater benefit on OS in the Epigenetic-Mut group than that in the Epigenetic-Wt group. (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4D</bold>
</xref>, mOS: 10.80 vs. 6.07 months, Log_rank <italic>P</italic> = 0.003, HR = 0.4279). In addition, we screened 9 genes with high mutation frequency from all epigenetic-related genes, whose predictive value has been demonstrated in other solid tumors, including ARID1A, ATRX, KMT2A/B/C/D, and TET1/2/3. The results showed that MSS-CRC with these gene mutations had more considerable ORR (<xref ref-type="supplementary-material" rid="SM4">
<bold>Supplementary Table&#xa0;4</bold>
</xref>, 8/16, 50%) and DCR (<xref ref-type="supplementary-material" rid="SM4">
<bold>Supplementary Table&#xa0;4</bold>
</xref>, 13/16, 81.25%).</p>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>Patient and disease characteristics of the validation set of MSS-CRC patients receiving ICI therapy.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="left">Characteristics</th>
<th valign="top" align="center">Epigenetic_Mut(n = 24)</th>
<th valign="top" align="center">Epigenetic_Wt (n = 65)</th>
<th valign="top" align="center">P-value*</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">
<bold>Age</bold>
</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.651</td>
</tr>
<tr>
<td valign="top" align="left">&lt;60</td>
<td valign="top" align="center">12</td>
<td valign="top" align="center">36</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">&#x2265;60</td>
<td valign="top" align="center">12</td>
<td valign="top" align="center">29</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>Sex</bold>
</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.423</td>
</tr>
<tr>
<td valign="top" align="left">Male</td>
<td valign="top" align="center">13</td>
<td valign="top" align="center">36</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">Female</td>
<td valign="top" align="center">11</td>
<td valign="top" align="center">29</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>ECOG PS</bold>
</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.857</td>
</tr>
<tr>
<td valign="top" align="left">0</td>
<td valign="top" align="center">16</td>
<td valign="top" align="center">42</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">&#x2265;1</td>
<td valign="top" align="center">6</td>
<td valign="top" align="center">23</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>ICI line</bold>
</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.967</td>
</tr>
<tr>
<td valign="top" align="left">1</td>
<td valign="top" align="center">0</td>
<td valign="top" align="center">2</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">2</td>
<td valign="top" align="center">6</td>
<td valign="top" align="center">16</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">&#x2265;3</td>
<td valign="top" align="center">18</td>
<td valign="top" align="center">47</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>Primary tumor sidedness</bold>
</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.334</td>
</tr>
<tr>
<td valign="top" align="left">Right</td>
<td valign="top" align="center">10</td>
<td valign="top" align="center">20</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">Left</td>
<td valign="top" align="center">14</td>
<td valign="top" align="center">45</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>Liver metastases</bold>
</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.683</td>
</tr>
<tr>
<td valign="top" align="left">With Liver</td>
<td valign="top" align="center">10</td>
<td valign="top" align="center">24</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">Without Liver</td>
<td valign="top" align="center">14</td>
<td valign="top" align="center">41</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>Regimen</bold>
</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.951</td>
</tr>
<tr>
<td valign="top" align="left">ICIs + TKIs</td>
<td valign="top" align="center">10</td>
<td valign="top" align="center">26</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">ICIs + Chemotherapy</td>
<td valign="top" align="center">11</td>
<td valign="top" align="center">32</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">ICIs + Chemoradiotherapy</td>
<td valign="top" align="center">3</td>
<td valign="top" align="center">7</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">
<bold>Best overall response</bold>**</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">0.026</td>
</tr>
<tr>
<td valign="top" align="left">CR/PR</td>
<td valign="top" align="center">9</td>
<td valign="top" align="center">10</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">SD</td>
<td valign="top" align="center">7</td>
<td valign="top" align="center">14</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">PD</td>
<td valign="top" align="center">8</td>
<td valign="top" align="center">41</td>
<td valign="top" align="center"/>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>* Fisher&#x2019;s exact test or Wilcoxon-Mann-Whitney test, as appropriate.</p>
</fn>
<fn>
<p>** CR, complete response; PR, partial response; SD, stable disease; PD, progressive disease.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Epigenetic_Mut predicts favorable clinical outcomes following ICI therapy. <bold>(A)</bold> Histogram presenting the proportion of patients that acquired ORR in the Epigenetic_Mut and Epigenetic_Wt groups. <bold>(B)</bold> Histogram presenting the proportion of patients that acquired DCR in the Epigenetic_Mut and Epigenetic_Wt groups. <bold>(C)</bold> Kaplan&#x2013;Meier estimates of PFS between Epigenetic_Mut or Epigenetic_Wt group patients in the discovery cohort. <bold>(D)</bold> Kaplan&#x2013;Meier estimates of OS between Epigenetic_Mut or Epigenetic_Wt group patients in the discovery cohort.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-13-1039631-g004.tif"/>
</fig>
</sec>
<sec id="s3_5">
<title>The abundance of immune cells in tumor tissue and peripheral blood of patients with or without epigenetic-related gene mutation</title>
<p>We explored the densities of CD8<sup>+</sup> and FOXP3<sup>+</sup> cells in MSS-CRC samples with different epigenetic-related gene statuses using IHC. Of the 34 MSS-CRC samples, 10 had epigenetic gene mutations. Further, we captured representative images of CD8<sup>+</sup> cells and FOXP3<sup>+</sup> cells from three samples. The first patient had an <italic>ARID1A</italic> mutation (<italic>ARID1A Frame_Shift_Del</italic>), and the second patient had a <italic>KMT2D</italic> mutation (<italic>KMT2D Nonsense_mutation</italic>). Both samples showed increased CD8<sup>+</sup> cell density and decreased FOXP3<sup>+</sup> cell density in tumor tissues (<xref ref-type="fig" rid="f5">
<bold>Figures&#xa0;5A, B</bold>
</xref>). However, in the third patient, who did not present any epigenetic-related gene mutations, the density of CD8<sup>+</sup> lymphocytes was lower, and the density of FOXP3<sup>+</sup> lymphocytes was higher than that of the other two (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5C</bold>
</xref>). CD8<sup>+</sup> and FOXP3<sup>+</sup> cell densities were counted in 38 patients (Epigenetic_Mut, N = 10; Epigenetic_Wt, N = 28), and from this we discovered that CD8<sup>+</sup> cell density increased in the Epigenetic_Mut group (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5D</bold>
</xref>) and that the FOXP3<sup>+</sup> cell density decreased in the Epigenetic_Mut group (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5E</bold>
</xref>). Furthermore, the ratio of CD8/FOXP3 cells in the Epigenetic_Mut group was significantly higher than that in the Epigenetic_Wt group (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5F</bold>
</xref>). Next, we collected peripheral blood from 12 patients with MSS-CRC (3 with epigenetic-related gene mutations) and measured the proportion of CD8<sup>+</sup>PD1<sup>+</sup>T cells and CD3<sup>-</sup>CD56<sup>+</sup>CD16<sup>+</sup>NK cells by flow cytometry. We found that both the proportion of CD8<sup>+</sup>PD1<sup>+</sup>T cells and CD3<sup>-</sup>CD56<sup>+</sup>CD16<sup>+</sup>NK cells was higher in the Epigenetic_Mut group (<xref ref-type="fig" rid="f6">
<bold>Figures&#xa0;6A, B</bold>
</xref>).</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Infiltration of CD8<sup>+</sup> and FOXP3<sup>+</sup> lymphocytes in the tumors of patients with or without epigenetic-related gene mutations. <bold>(A)</bold> A Representative image of CD8<sup>+</sup> and FOXP3<sup>+</sup> lymphocytes infiltrating the MSS-CRC with ARID1A Frame_Shift_Del. <bold>(B)</bold> A Representative image of CD8<sup>+</sup> and FOXP3<sup>+</sup> lymphocytes infiltrating the MSS-CRC with KMT2D Nonsense_Mutation. <bold>(C)</bold> A Representative image of CD8<sup>+</sup> and FOXP3<sup>+</sup> lymphocytes infiltrating the MSS-CRC without epigenetic-related gene mutations. <bold>(D)</bold> Tumors with epigenetic-related genes mutation had significantly higher levels of intra-tumoral CD8<sup>+</sup> lymphocytes than tumors with wild-type epigenetic-related genes. <bold>(E)</bold> Tumors with epigenetic-related genes mutation had significantly lower levels of intra-tumoral FOXP3<sup>+</sup> lymphocytes than tumors with wild-type epigenetic-related genes. <bold>(F)</bold> The Epigenetic_Mut group had a higher CD8/FOXP3 cell ratio than the Epigenetic_Wt group.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-13-1039631-g005.tif"/>
</fig>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Proportion of CD8<sup>+</sup>PD1+T cells and NK cells in the peripheral blood of patients with or without epigenetic-related gene mutations. <bold>(A)</bold> The Epigenetic_Mut group had a higher proportion of CD8<sup>+</sup>PD1<sup>+</sup>T cells compared to the Epigenetic_Wt group in peripheral blood. <bold>(B)</bold> The Epigenetic_Mut group had a higher proportion of NK cells compared to Epigenetic_Wt group in peripheral blood.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fimmu-13-1039631-g006.tif"/>
</fig>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<title>Discussion</title>
<p>Although ICI-based combination therapies have shown certain effectiveness in pMMR/MSS CRC, especially in combination with antiangiogenic agents (Lenvatinib or Regorafenib) that resulted in an ORR of 20-30% (<xref ref-type="bibr" rid="B6">6</xref>, <xref ref-type="bibr" rid="B33">33</xref>), most patients still cannot benefit from the combination therapy because of the high heterogeneity of pMMR/MSS CRC. Recently, the MAYA phase II trial (NCT03832621) showed that MSS-CRC patients with silenced MGMT could benefit from ICIs combined with temozolomide treatment (<xref ref-type="bibr" rid="B34">34</xref>). This trial showed 36% for 8-month PFS, 42% for ORR, and 18.4 months for the median OS. Therefore, screening the MSS CRC patients with active anti-tumor immune response may be the key to improving the efficacy of immunotherapy. However, the predictive biomarkers for ICI therapy in MSS-CRC patients are limited.</p>
<p>pMMR/MSS CRC is a cold tumor that contains few neoantigens and either no or inactive TILs (<xref ref-type="bibr" rid="B35">35</xref>). Meanwhile, CRC is a multilayered heterogeneous disease with specific treatment challenges and opportunities (<xref ref-type="bibr" rid="B36">36</xref>). Additionally, Previous studies have reported that epigenetic-related gene mutations affect both the tumor microenvironment and efficacy of ICIs (<xref ref-type="bibr" rid="B24">24</xref>&#x2013;<xref ref-type="bibr" rid="B26">26</xref>). Mechanistically, epigenetic modification can reshape the tumor microenvironment by affecting genomic instability and enhancing the immunogenicity of tumor cells. First, epigenetic modification can affect the DNA damage repair response by regulating the accessibility of chromatin. Studies have shown that epigenetic-related gene mutations can lead to increased TMB in tumor cells, such as ARID1A and KMT2D. ARID1A specifically has a 6.7% mutation rate in MSS-CRC (<xref ref-type="bibr" rid="B22">22</xref>) and may increase the instability of the genome by adjusting the MMR pathway (<xref ref-type="bibr" rid="B21">21</xref>, <xref ref-type="bibr" rid="B37">37</xref>). Mutations in the KMT2D gene are common in cancer patients, and their deficiency can increase the levels of genomic DNA damage and TMB, as well as increase transcription instability. Clinical studies have shown that individuals with mutations in genes from the KMT family are more likely to benefit from ICI therapy (<xref ref-type="bibr" rid="B24">24</xref>, <xref ref-type="bibr" rid="B25">25</xref>). Furthermore, epigenetic-related gene mutation enhances the immunogenicity of tumor cells. Accounting for 5%-10% of genomic DNA sequences, human endogenous retroviruses (ERVs) are remnants of the evolution of germline integrations of exogenous infectious retroviruses (<xref ref-type="bibr" rid="B38">38</xref>, <xref ref-type="bibr" rid="B39">39</xref>). These exogenous genes are not expressed in healthy tissues other than germ cells but are often abnormally expressed in tumors with epigenetic regulation defects. Here, neoantigen expression increases immunogenicity and triggers an innate immune response against tumors (<xref ref-type="bibr" rid="B40">40</xref>, <xref ref-type="bibr" rid="B41">41</xref>). Recently, genome-wide technologies have revealed frequent mutations in epigenetic modifier genes, particularly in cancers (<xref ref-type="bibr" rid="B42">42</xref>). It is therefore necessary to analyze systematically the immune activity and the effect of immunotherapy in MSS-CRC patients with epigenetic regulation impairment.</p>
<p>In our study, we systematically analyzed 68 epigenetic-related genes from 13 pathways involved in chromatin regulatory processes in MSS-CRC samples. The mutation rate of epigenetic-related genes in the TCGA cohort was 18.35%. This mutation frequency was higher than that of any previous marker in the population, such as <italic>POLE</italic> or DDR mutations, and was closer to the potential benefit ratio in MSS-CRC clinical trials. <italic>ARID1A</italic>, <italic>KMT2C</italic>, <italic>RSF1</italic>, <italic>CHD9, PBRM1</italic>, and <italic>ATRX</italic> were the most mutated genes in the TCGA cohort, accounting for approximately 75% of the epigenetic-related gene-mutated MSS-CRC patients. This is consistent with previous reports, and <italic>ARID1A</italic> is thus a marker gene that should be investigated in clinical practice.</p>
<p>Using bioinformatics algorithms, we also assessed whether the MSS-CRC samples with epigenetic-related gene mutations from TCGA had better immune activity, including immune signatures, tumor-infiltrating lymphocytes, and expression of immune checkpoints and key genes. Furthermore, we validated our bioinformatic findings using immunohistochemical analyses of CD8<sup>+</sup> and FOXP3<sup>+</sup> cells from a cohort of MSS-CRC patients, and similar results were obtained at the histopathological level. In the Epigenetic _Mut group, CD8<sup>+</sup> cells were higher and FOXP3<sup>+</sup> cells were lower. The Epigenetic _Mut group also had a higher proportion of CD8/FOXP3 cells than the Epigenetic_Wt group. The VOLTAGE trial demonstrated that among MSS-CRC patients receiving ICIs as neoadjuvant treatment, patients with an elevated CD8/FOXP3 cell ratio were more likely to achieve pathologic complete response (pCR), suggesting that the CD8/FOXP3 cell ratio may be a predictor for ICI therapy efficacy (<xref ref-type="bibr" rid="B43">43</xref>).</p>
<p>Finally, we validated the predictive power of epigenetic-related gene mutations in the HMUCH cohort of 89 MSS CRC patients who received immunotherapy and discovered that patients with epigenetic mutations were more likely to benefit from ICI-based combination therapy and had better clinical outcomes. These preliminary results demonstrate that epigenetic-related gene mutations can predict the response to ICIs in MSS-CRC patients.</p>
<p>This study has several limitations, including the validation cohort coming from a single-center, the small size of the cohort, and the lack of validation in other populations. This is because ICI-based regimens have not been recommended by any clinical guidelines for MSS-CRC. Numerous patients included in this study experienced the failure of standard treatment, and the treatment compliance and completeness of the clinical information in many of these patients, were not ideal. Additionally, since the genetic information in the HMUCH cohort was obtained from clinical testing, transcriptomic data were lacking. Thus, our TCGA cohort findings could not be validated. Instead, we performed immunohistochemical staining analysis of pathological sections to validate the immune activation status of the Epigenetic _Mut group, but a larger-scale validation remains necessary. Furthermore, the application of ICIs in MSS-CRC has not been standardized, and most patients enrolled in our study were patients who had experienced multiple failed lines of treatment, bringing considerable heterogeneity to the population of this study. Therefore, future prospective studies with larger cohort studies are needed.</p>
</sec>
<sec id="s5">
<title>Conclusion</title>
<p>In conclusion, our data suggest that identifying epigenetic-related gene mutations might help select the right immunotherapy for MSS-CRC patients and can be used as a biomarker to predict ICI therapy effectiveness. Importantly, the status of epigenetic-related gene mutations is highly accessible from clinical genetic testing, although it is often overlooked by clinicians. Further exploration of the molecular mechanisms underlying the increased effectiveness in specific MSS-CRC patients and prospective clinical trials are therefore warranted.</p>
</sec>
<sec id="s6" sec-type="data-availability">
<title>Data availability statement</title>
<p>The data presented in the study are deposited in the Genome Sequence Archive in National Genomics Data Center repository (<uri xlink:href="https://ngdc.cncb.ac.cn/gsa-human">https://ngdc.cncb.ac.cn/gsa-human</uri>), accession number GSA-Human: HRA003408.</p>
</sec>
<sec id="s7" sec-type="ethics-statement">
<title>Ethics statement</title>
<p>The study was approved by the Ethics Committee of Harbin Medical University Cancer Hospital (Harbin, China). Informed consents were obtained from all patients before surgery and during the experimental procedures.</p>
</sec>
<sec id="s8" sec-type="author-contributions">
<title>Author contributions</title>
<p>CL: Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Software, Validation, Visualization, Writing - original draft, Writing - review &amp; editing. HX: Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing - original draft, Writing - review &amp; editing. LC: Investigation, Methodology, Validation, Writing - original draft, Writing - review &amp; editing. LF: Investigation, Methodology, Writing - review &amp; editing. SH: Investigation, Writing - review &amp; editing. YR: Investigation, Writing - review &amp; editing. WZ: Conceptualization, Project administration, Resources, Supervision, Writing - review &amp; editing. YZ: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing - review &amp; editing. The work reported in the paper has been performed by the authors, unless clearly specified in the text. All authors contributed to the article and approved the submitted version.</p>
</sec>
<sec id="s9" sec-type="funding-information">
<title>Funding</title>
<p>This work was supported by grants from the National Natural Science Foundation of China (No. 82102858), the Natural Science Funding of Heilongjiang (No. YQ2022H017), nN10 Excellent Discipline Construction Program (No. Hepatobiliary and Pancreatic Tumor 2017); Basic scientific research projects of universities in Heilongjiang Province (No. 2020-KYYWF-1465); Top Young Talents Project of HMUCH (No. BJQN2021-01); Haiyan Research Fund of Harbin Medical University Cancer Hospital (No. JJMS2022-04), and Beijing Science and Technology Innovation Medical Development Foundation (No. KC2021-JX-0186-102, KC2021-JX-0186-69).</p>
</sec>
<sec id="s10" sec-type="acknowledgement">
<title>Acknowledgments</title>
<p>The results of the pan-cancer analysis were based on data generated by TCGA Research Network: <uri xlink:href="https://www.cancer.gov/tcga">https://www.cancer.gov/tcga</uri>. The authors thank AiMi Academic Services (<uri xlink:href="http://www.aimieditor.com">www.aimieditor.com</uri>) for English language editing services.</p>
</sec>
<sec id="s11" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="s12" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
</body>
<back>
<sec id="s13" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fimmu.2022.1039631/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fimmu.2022.1039631/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Image_1.pdf" id="SF1" mimetype="application/pdf"/>
<supplementary-material xlink:href="Table_1.docx" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document"/>
<supplementary-material xlink:href="Table_2.docx" id="SM2" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document"/>
<supplementary-material xlink:href="Table_3.docx" id="SM3" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document"/>
<supplementary-material xlink:href="Table_4.docx" id="SM4" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lenz</surname> <given-names>HJ</given-names>
</name>
<name>
<surname>Van Cutsem</surname> <given-names>E</given-names>
</name>
<name>
<surname>Luisa Limon</surname> <given-names>M</given-names>
</name>
<name>
<surname>Wong</surname> <given-names>KYM</given-names>
</name>
<name>
<surname>Hendlisz</surname> <given-names>A</given-names>
</name>
<name>
<surname>Aglietta</surname> <given-names>M</given-names>
</name>
<etal/>
</person-group>. <article-title>First-line nivolumab plus low-dose ipilimumab for microsatellite instability-high/mismatch repair-deficient metastatic colorectal cancer: the phase II CheckMate 142 study</article-title>. <source>J Clin Oncol</source> (<year>2022</year>) <volume>40</volume>:<page-range>161&#x2013;70</page-range>. doi: <pub-id pub-id-type="doi">10.1200/JCO.21.01015</pub-id>
</citation>
</ref>
<ref id="B2">
<label>2</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Le</surname> <given-names>DT</given-names>
</name>
<name>
<surname>Kim</surname> <given-names>TW</given-names>
</name>
<name>
<surname>Van Cutsem</surname> <given-names>E</given-names>
</name>
<name>
<surname>Geva</surname> <given-names>R</given-names>
</name>
<name>
<surname>J&#xe4;ger</surname> <given-names>D</given-names>
</name>
<name>
<surname>Hara</surname> <given-names>H</given-names>
</name>
<etal/>
</person-group>. <article-title>Phase II open-label study of pembrolizumab in treatment-refractory, microsatellite instability-high/mismatch repair-deficient metastatic colorectal cancer: KEYNOTE-164</article-title>. <source>J Clin Oncol</source> (<year>2020</year>) <volume>38</volume>:<page-range>11&#x2013;9</page-range>. doi: <pub-id pub-id-type="doi">10.1200/JCO.19.02107</pub-id>
</citation>
</ref>
<ref id="B3">
<label>3</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ganesh</surname> <given-names>K</given-names>
</name>
<name>
<surname>Stadler</surname> <given-names>ZK</given-names>
</name>
<name>
<surname>Cercek</surname> <given-names>A</given-names>
</name>
<name>
<surname>Mendelsohn</surname> <given-names>RB</given-names>
</name>
<name>
<surname>Shia</surname> <given-names>J</given-names>
</name>
<name>
<surname>Segal</surname> <given-names>NH</given-names>
</name>
<etal/>
</person-group>. <article-title>Immunotherapy in colorectal cancer: rationale, challenges and potential</article-title>. <source>Nat Rev Gastroenterol Hepatol</source> (<year>2019</year>) <volume>16</volume>:<page-range>361&#x2013;75</page-range>. doi: <pub-id pub-id-type="doi">10.1038/s41575-019-0126-x</pub-id>
</citation>
</ref>
<ref id="B4">
<label>4</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Eng</surname> <given-names>C</given-names>
</name>
<name>
<surname>Kim</surname> <given-names>TW</given-names>
</name>
<name>
<surname>Bendell</surname> <given-names>J</given-names>
</name>
<name>
<surname>Argil&#xe9;s</surname> <given-names>G</given-names>
</name>
<name>
<surname>Tebbutt</surname> <given-names>NC</given-names>
</name>
<name>
<surname>Di Bartolomeo</surname> <given-names>M</given-names>
</name>
<etal/>
</person-group>. <article-title>Atezolizumab with or without cobimetinib versus regorafenib in previously treated metastatic colorectal cancer (IMblaze370): a multicentre, open-label, phase 3, randomised, controlled trial</article-title>. <source>Lancet Oncol</source> (<year>2019</year>) <volume>20</volume>:<page-range>849&#x2013;61</page-range>. doi: <pub-id pub-id-type="doi">10.1016/S1470-2045(19)30027-0</pub-id>
</citation>
</ref>
<ref id="B5">
<label>5</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chalabi</surname> <given-names>M</given-names>
</name>
<name>
<surname>Fanchi</surname> <given-names>LF</given-names>
</name>
<name>
<surname>Dijkstra</surname> <given-names>KK</given-names>
</name>
<name>
<surname>Van den Berg</surname> <given-names>JG</given-names>
</name>
<name>
<surname>Aalbers</surname> <given-names>AG</given-names>
</name>
<name>
<surname>Sikorska</surname> <given-names>K</given-names>
</name>
<etal/>
</person-group>. <article-title>Neoadjuvant immunotherapy leads to pathological responses in MMR-proficient and MMR-deficient early-stage colon cancers</article-title>. <source>Nat Med</source> (<year>2020</year>) <volume>26</volume>:<page-range>566&#x2013;76</page-range>. doi: <pub-id pub-id-type="doi">10.1038/s41591-020-0805-8</pub-id>
</citation>
</ref>
<ref id="B6">
<label>6</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fukuoka</surname> <given-names>S</given-names>
</name>
<name>
<surname>Hara</surname> <given-names>H</given-names>
</name>
<name>
<surname>Takahashi</surname> <given-names>N</given-names>
</name>
<name>
<surname>Kojima</surname> <given-names>T</given-names>
</name>
<name>
<surname>Kawazoe</surname> <given-names>A</given-names>
</name>
<name>
<surname>Asayama</surname> <given-names>M</given-names>
</name>
<etal/>
</person-group>. <article-title>Regorafenib plus nivolumab in patients with advanced gastric or colorectal cancer: an open-label, dose-escalation, and dose-expansion phase ib trial (REGONIVO, EPOC1603)</article-title>. <source>J Clin Oncol</source> (<year>2020</year>) <volume>38</volume>:<page-range>2053&#x2013;61</page-range>. doi: <pub-id pub-id-type="doi">10.1200/JCO.19.03296</pub-id>
</citation>
</ref>
<ref id="B7">
<label>7</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>F</given-names>
</name>
<name>
<surname>He</surname> <given-names>MM</given-names>
</name>
<name>
<surname>Yao</surname> <given-names>YC</given-names>
</name>
<name>
<surname>Zhao</surname> <given-names>X</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>ZQ</given-names>
</name>
<name>
<surname>Jin</surname> <given-names>Y</given-names>
</name>
<etal/>
</person-group>. <article-title>Regorafenib plus toripalimab in patients with metastatic colorectal cancer: a phase Ib/II clinical trial and gut microbiome analysis</article-title>. <source>Cell Rep Med</source> (<year>2021</year>) <volume>2</volume>:<fpage>100383</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.xcrm.2021.100383</pub-id>
</citation>
</ref>
<ref id="B8">
<label>8</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>O'Neil</surname> <given-names>BH</given-names>
</name>
<name>
<surname>Wallmark</surname> <given-names>JM</given-names>
</name>
<name>
<surname>Lorente</surname> <given-names>D</given-names>
</name>
<name>
<surname>Elez</surname> <given-names>E</given-names>
</name>
<name>
<surname>Raimbourg</surname> <given-names>J</given-names>
</name>
<name>
<surname>Gomez-Roca</surname> <given-names>C</given-names>
</name>
<etal/>
</person-group>. <article-title>Safety and antitumor activity of the anti-PD-1 antibody pembrolizumab in patients with advanced colorectal carcinoma</article-title>. <source>PLos One</source> (<year>2017</year>) <volume>12</volume>:<fpage>e0189848</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.pone.0189848</pub-id>
</citation>
</ref>
<ref id="B9">
<label>9</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>F</given-names>
</name>
<name>
<surname>Zhao</surname> <given-names>Q</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>YN</given-names>
</name>
<name>
<surname>Jin</surname> <given-names>Y</given-names>
</name>
<name>
<surname>He</surname> <given-names>MM</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>ZX</given-names>
</name>
<etal/>
</person-group>. <article-title>Evaluation of POLE and POLD1 mutations as biomarkers for immunotherapy outcomes across multiple cancer types</article-title>. <source>JAMA Oncol</source> (<year>2019</year>) <volume>5</volume>:<page-range>1504&#x2013;6</page-range>. doi: <pub-id pub-id-type="doi">10.1001/jamaoncol.2019.2963</pub-id>
</citation>
</ref>
<ref id="B10">
<label>10</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ott</surname> <given-names>PA</given-names>
</name>
<name>
<surname>Bang</surname> <given-names>YJ</given-names>
</name>
<name>
<surname>Piha-Paul</surname> <given-names>SA</given-names>
</name>
<name>
<surname>Razak</surname> <given-names>ARA</given-names>
</name>
<name>
<surname>Bennouna</surname> <given-names>J</given-names>
</name>
<name>
<surname>Soria</surname> <given-names>JC</given-names>
</name>
<etal/>
</person-group>. <article-title>T-Cell-inflamed gene-expression profile, programmed death ligand 1 expression, and tumor mutational burden predict efficacy in patients treated with pembrolizumab across 20 cancers: KEYNOTE-028</article-title>. <source>J Clin Oncol</source> (<year>2019</year>) <volume>37</volume>:<page-range>318&#x2013;27</page-range>. doi: <pub-id pub-id-type="doi">10.1200/JCO.2018.78.2276</pub-id>
</citation>
</ref>
<ref id="B11">
<label>11</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Turajlic</surname> <given-names>S</given-names>
</name>
<name>
<surname>Litchfield</surname> <given-names>K</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>H</given-names>
</name>
<name>
<surname>Rosenthal</surname> <given-names>R</given-names>
</name>
<name>
<surname>McGranahan</surname> <given-names>N</given-names>
</name>
<name>
<surname>Reading</surname> <given-names>JL</given-names>
</name>
<etal/>
</person-group>. <article-title>Insertion-and-deletion-derived tumour-specific neoantigens and the immunogenic phenotype: a pan-cancer analysis</article-title>. <source>Lancet Oncol</source> (<year>2017</year>) <volume>18</volume>:<page-range>1009&#x2013;21</page-range>. doi: <pub-id pub-id-type="doi">10.1016/S1470-2045(17)30516-8</pub-id>
</citation>
</ref>
<ref id="B12">
<label>12</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ayers</surname> <given-names>M</given-names>
</name>
<name>
<surname>Lunceford</surname> <given-names>J</given-names>
</name>
<name>
<surname>Nebozhyn</surname> <given-names>M</given-names>
</name>
<name>
<surname>Murphy</surname> <given-names>E</given-names>
</name>
<name>
<surname>Loboda</surname> <given-names>A</given-names>
</name>
<name>
<surname>Kaufman</surname> <given-names>DR</given-names>
</name>
<etal/>
</person-group>. <article-title>IFN-gamma-related mRNA profile predicts clinical response to PD-1 blockade</article-title>. <source>J Clin Invest</source> (<year>2017</year>) <volume>127</volume>:<page-range>2930&#x2013;40</page-range>. doi: <pub-id pub-id-type="doi">10.1172/JCI91190</pub-id>
</citation>
</ref>
<ref id="B13">
<label>13</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bortolomeazzi</surname> <given-names>M</given-names>
</name>
<name>
<surname>Keddar</surname> <given-names>MR</given-names>
</name>
<name>
<surname>Montorsi</surname> <given-names>L</given-names>
</name>
<name>
<surname>Acha-Sagredo</surname> <given-names>A</given-names>
</name>
<name>
<surname>Benedetti</surname> <given-names>L</given-names>
</name>
<name>
<surname>Temelkovski</surname> <given-names>D</given-names>
</name>
<etal/>
</person-group>. <article-title>Immunogenomics of colorectal cancer response to checkpoint blockade: analysis of the KEYNOTE 177 trial and validation cohorts</article-title>. <source>Gastroenterology</source> (<year>2021</year>) <volume>161</volume>:<page-range>1179&#x2013;93</page-range>. doi: <pub-id pub-id-type="doi">10.1053/j.gastro.2021.06.064</pub-id>
</citation>
</ref>
<ref id="B14">
<label>14</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Song</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Huang</surname> <given-names>J</given-names>
</name>
<name>
<surname>Liang</surname> <given-names>D</given-names>
</name>
<name>
<surname>Hu</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Mao</surname> <given-names>B</given-names>
</name>
<name>
<surname>Li</surname> <given-names>Q</given-names>
</name>
<etal/>
</person-group>. <article-title>DNA Damage repair gene mutations are indicative of a favorable prognosis in colorectal cancer treated with immune checkpoint inhibitors</article-title>. <source>Front Oncol</source> (<year>2020</year>) <volume>10</volume>:<elocation-id>549777</elocation-id>. doi: <pub-id pub-id-type="doi">10.3389/fonc.2020.549777</pub-id>
</citation>
</ref>
<ref id="B15">
<label>15</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Arai</surname> <given-names>H</given-names>
</name>
<name>
<surname>Elliott</surname> <given-names>A</given-names>
</name>
<name>
<surname>Xiu</surname> <given-names>J</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>J</given-names>
</name>
<name>
<surname>Battaglin</surname> <given-names>F</given-names>
</name>
<name>
<surname>Kawanishi</surname> <given-names>N</given-names>
</name>
<etal/>
</person-group>. <article-title>The landscape of alterations in DNA damage response pathways in colorectal cancer</article-title>. <source>Clin Cancer Res</source> (<year>2021</year>) <volume>27</volume>:<page-range>3234&#x2013;42</page-range>. doi: <pub-id pub-id-type="doi">10.1158/1078-0432.CCR-20-3635</pub-id>
</citation>
</ref>
<ref id="B16">
<label>16</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname> <given-names>P</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>X</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>H</given-names>
</name>
<name>
<surname>Hu</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>H</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>L</given-names>
</name>
<etal/>
</person-group>. <article-title>The mutational pattern of homologous recombination-related (HRR) genes in Chinese colon cancer and its relevance to immunotherapy responses</article-title>. <source>Aging (Albany NY)</source> (<year>2020</year>) <volume>13</volume>:<page-range>2365&#x2013;78</page-range>. doi: <pub-id pub-id-type="doi">10.18632/aging.202267</pub-id>
</citation>
</ref>
<ref id="B17">
<label>17</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Heeke</surname> <given-names>AL</given-names>
</name>
<name>
<surname>Pishvaian</surname> <given-names>MJ</given-names>
</name>
<name>
<surname>Lynce</surname> <given-names>F</given-names>
</name>
<name>
<surname>Xiu</surname> <given-names>J</given-names>
</name>
<name>
<surname>Brody</surname> <given-names>JR</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>WJ</given-names>
</name>
<etal/>
</person-group>. <article-title>Prevalence of homologous recombination-related gene mutations across multiple cancer types</article-title>. <source>JCO Precis Oncol</source> (<year>2018</year>) <volume>2018</volume>:<page-range>1&#x2013;13</page-range>. doi: <pub-id pub-id-type="doi">10.1200/PO.17.00286</pub-id>
</citation>
</ref>
<ref id="B18">
<label>18</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Valencia</surname> <given-names>AM</given-names>
</name>
<name>
<surname>Kadoch</surname> <given-names>C</given-names>
</name>
</person-group>. <article-title>Chromatin regulatory mechanisms and therapeutic opportunities in cancer</article-title>. <source>Nat Cell Biol</source> (<year>2019</year>) <volume>21</volume>:<page-range>152&#x2013;61</page-range>. doi: <pub-id pub-id-type="doi">10.1038/s41556-018-0258-1</pub-id>
</citation>
</ref>
<ref id="B19">
<label>19</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Clapier</surname> <given-names>CR</given-names>
</name>
<name>
<surname>Cairns</surname> <given-names>BR</given-names>
</name>
</person-group>. <article-title>The biology of chromatin remodeling complexes</article-title>. <source>Annu Rev Biochem</source> (<year>2009</year>) <volume>78</volume>:<fpage>273</fpage>&#x2013;<lpage>304</lpage>. doi: <pub-id pub-id-type="doi">10.1146/annurev.biochem.77.062706.153223</pub-id>
</citation>
</ref>
<ref id="B20">
<label>20</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Goswami</surname> <given-names>S</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Anandhan</surname> <given-names>S</given-names>
</name>
<name>
<surname>Szabo</surname> <given-names>PM</given-names>
</name>
<name>
<surname>Basu</surname> <given-names>S</given-names>
</name>
<name>
<surname>Blando</surname> <given-names>JM</given-names>
</name>
<etal/>
</person-group>. <article-title>ARID1A mutation plus CXCL13 expression act as combinatorial biomarkers to predict responses to immune checkpoint therapy in mUCC</article-title>. <source>Sci Transl Med</source> (<year>2020</year>) <volume>12</volume>:<elocation-id>eabc4220</elocation-id>. doi: <pub-id pub-id-type="doi">10.1126/scitranslmed.abc4220</pub-id>
</citation>
</ref>
<ref id="B21">
<label>21</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>L</given-names>
</name>
<name>
<surname>Qu</surname> <given-names>J</given-names>
</name>
<name>
<surname>Zhou</surname> <given-names>N</given-names>
</name>
<name>
<surname>Hou</surname> <given-names>H</given-names>
</name>
<name>
<surname>Jiang</surname> <given-names>M</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>X</given-names>
</name>
</person-group>. <article-title>Effect and biomarker of immune checkpoint blockade therapy for ARID1A deficiency cancers</article-title>. <source>BioMed Pharmacother</source> (<year>2020</year>) <volume>130</volume>:<fpage>110626</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.biopha.2020.110626</pub-id>
</citation>
</ref>
<ref id="B22">
<label>22</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mehrvarz Sarshekeh</surname> <given-names>A</given-names>
</name>
<name>
<surname>Alshenaifi</surname> <given-names>J</given-names>
</name>
<name>
<surname>Roszik</surname> <given-names>J</given-names>
</name>
<name>
<surname>Manyam</surname> <given-names>GC</given-names>
</name>
<name>
<surname>Advani</surname> <given-names>SM</given-names>
</name>
<name>
<surname>Katkhuda</surname> <given-names>R</given-names>
</name>
<etal/>
</person-group>. <article-title>ARID1A mutation may define an immunologically active subgroup in patients with microsatellite stable colorectal cancer</article-title>. <source>Clin Cancer Res</source> (<year>2021</year>) <volume>27</volume>:<page-range>1663&#x2013;70</page-range>. doi: <pub-id pub-id-type="doi">10.1158/1078-0432.CCR-20-2404</pub-id>
</citation>
</ref>
<ref id="B23">
<label>23</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ruthenburg</surname> <given-names>AJ</given-names>
</name>
<name>
<surname>Allis</surname> <given-names>CD</given-names>
</name>
<name>
<surname>Wysocka</surname> <given-names>J</given-names>
</name>
</person-group>. <article-title>Methylation of lysine 4 on histone H3: intricacy of writing and reading a single epigenetic mark</article-title>. <source>Mol Cell</source> (<year>2007</year>) <volume>25</volume>:<fpage>15</fpage>&#x2013;<lpage>30</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.molcel.2006.12.014</pub-id>
</citation>
</ref>
<ref id="B24">
<label>24</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname> <given-names>P</given-names>
</name>
<name>
<surname>Huang</surname> <given-names>Y</given-names>
</name>
</person-group>. <article-title>Genomic alterations in KMT2 family predict outcome of immune checkpoint therapy in multiple cancers</article-title>. <source>J Hematol Oncol</source> (<year>2021</year>) <volume>14</volume>:<fpage>39</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s13045-021-01050-0</pub-id>
</citation>
</ref>
<ref id="B25">
<label>25</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>G</given-names>
</name>
<name>
<surname>Chow</surname> <given-names>RD</given-names>
</name>
<name>
<surname>Zhu</surname> <given-names>L</given-names>
</name>
<name>
<surname>Bai</surname> <given-names>Z</given-names>
</name>
<name>
<surname>Ye</surname> <given-names>L</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>F</given-names>
</name>
<etal/>
</person-group>. <article-title>CRISPR-GEMM pooled mutagenic screening identifies KMT2D as a major modulator of immune checkpoint blockade</article-title>. <source>Cancer Discovery</source> (<year>2020</year>) <volume>10</volume>:<page-range>1912&#x2013;33</page-range>. doi: <pub-id pub-id-type="doi">10.1158/2159-8290.CD-19-1448</pub-id>
</citation>
</ref>
<ref id="B26">
<label>26</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname> <given-names>HX</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>YX</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>ZX</given-names>
</name>
<name>
<surname>Zhao</surname> <given-names>Q</given-names>
</name>
<name>
<surname>He</surname> <given-names>MM</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>YN</given-names>
</name>
<etal/>
</person-group>. <article-title>Alteration in TET1 as potential biomarker for immune checkpoint blockade in multiple cancers</article-title>. <source>J Immunother Cancer</source> (<year>2019</year>) <volume>7</volume>:<fpage>264</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s40425-019-0737-3</pub-id>
</citation>
</ref>
<ref id="B27">
<label>27</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Niu</surname> <given-names>B</given-names>
</name>
<name>
<surname>Ye</surname> <given-names>K</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Q</given-names>
</name>
<name>
<surname>Lu</surname> <given-names>C</given-names>
</name>
<name>
<surname>Xie</surname> <given-names>M</given-names>
</name>
<name>
<surname>McLellan</surname> <given-names>MD</given-names>
</name>
<etal/>
</person-group>. <article-title>MSIsensor: microsatellite instability detection using paired tumor-normal sequence data</article-title>. <source>Bioinformatics</source> (<year>2014</year>) <volume>30</volume>:<page-range>1015&#x2013;6</page-range>. doi: <pub-id pub-id-type="doi">10.1093/bioinformatics/btt755</pub-id>
</citation>
</ref>
<ref id="B28">
<label>28</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Thorsson</surname> <given-names>V</given-names>
</name>
<name>
<surname>Gibbs</surname> <given-names>DL</given-names>
</name>
<name>
<surname>Brown</surname> <given-names>SD</given-names>
</name>
<name>
<surname>Wolf</surname> <given-names>D</given-names>
</name>
<name>
<surname>Bortone</surname> <given-names>DS</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>THO</given-names>
</name>
<etal/>
</person-group>. <article-title>The immune landscape of cancer</article-title>. <source>Immunity</source> (<year>2018</year>) <volume>48</volume>:<page-range>812&#x2013;30</page-range>. doi: <pub-id pub-id-type="doi">10.1016/j.immuni.2018.03.023</pub-id>
</citation>
</ref>
<ref id="B29">
<label>29</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hoadley</surname> <given-names>KA</given-names>
</name>
<name>
<surname>Yau</surname> <given-names>C</given-names>
</name>
<name>
<surname>Hinoue</surname> <given-names>T</given-names>
</name>
<name>
<surname>Wolf</surname> <given-names>DM</given-names>
</name>
<name>
<surname>Lazar</surname> <given-names>AJ</given-names>
</name>
<name>
<surname>Drill</surname> <given-names>E</given-names>
</name>
<etal/>
</person-group>. <article-title>Cell-of-origin patterns dominate the molecular classification of 10,000 tumors from 33 types of cancer</article-title>. <source>Cell</source> (<year>2018</year>) <volume>173</volume>:<fpage>291</fpage>&#x2013;<lpage>304</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.cell.2018.03.022</pub-id>
</citation>
</ref>
<ref id="B30">
<label>30</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Davoli</surname> <given-names>T</given-names>
</name>
<name>
<surname>Uno</surname> <given-names>H</given-names>
</name>
<name>
<surname>Wooten</surname> <given-names>EC</given-names>
</name>
<name>
<surname>Elledge</surname> <given-names>SJ</given-names>
</name>
</person-group>. <article-title>Tumor aneuploidy correlates with markers of immune evasion and with reduced response to immunotherapy</article-title>. <source>Science</source> (<year>2017</year>) <volume>355</volume>:<page-range>1&#x2013;14</page-range>. doi: <pub-id pub-id-type="doi">10.1126/science.aaf8399</pub-id>
</citation>
</ref>
<ref id="B31">
<label>31</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Guinney</surname> <given-names>J</given-names>
</name>
<name>
<surname>Dienstmann</surname> <given-names>R</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>X</given-names>
</name>
<name>
<surname>de Reyni&#xe8;s</surname> <given-names>A</given-names>
</name>
<name>
<surname>Schlicker</surname> <given-names>A</given-names>
</name>
<name>
<surname>Soneson</surname> <given-names>C</given-names>
</name>
<etal/>
</person-group>. <article-title>The consensus molecular subtypes of colorectal cancer</article-title>. <source>Nat Med</source> (<year>2015</year>) <volume>21</volume>:<page-range>1350&#x2013;6</page-range>. doi: <pub-id pub-id-type="doi">10.1038/nm.3967</pub-id>
</citation>
</ref>
<ref id="B32">
<label>32</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qin</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Ekmekcioglu</surname> <given-names>S</given-names>
</name>
<name>
<surname>Forget</surname> <given-names>MA</given-names>
</name>
<name>
<surname>Szekvolgyi</surname> <given-names>L</given-names>
</name>
<name>
<surname>Hwu</surname> <given-names>P</given-names>
</name>
<name>
<surname>Grimm</surname> <given-names>EA</given-names>
</name>
<etal/>
</person-group>. <article-title>Cervical cancer neoantigen landscape and immune activity is associated with human papillomavirus master regulators</article-title>. <source>Front Immunol</source> (<year>2017</year>) <volume>8</volume>:<elocation-id>689</elocation-id>. doi: <pub-id pub-id-type="doi">10.3389/fimmu.2017.00689</pub-id>
</citation>
</ref>
<ref id="B33">
<label>33</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gomez-Roca</surname> <given-names>CA</given-names>
</name>
<name>
<surname>Yanez</surname> <given-names>E</given-names>
</name>
<name>
<surname>Im</surname> <given-names>S-A</given-names>
</name>
<name>
<surname>Castanon Alvarez</surname> <given-names>EC</given-names>
</name>
<name>
<surname>Senellart</surname> <given-names>H</given-names>
</name>
<name>
<surname>Doherty</surname> <given-names>M</given-names>
</name>
<etal/>
</person-group>. <article-title>LEAP-005: A phase 2 multicohort study of lenvatinib plus pembrolizumab in patients with previously treated selected solid tumors&#x2013;results from the colorectal cancer cohort</article-title>. <source>J Clin Oncol</source> (<year>2021</year>) <volume>39</volume>:<fpage>3564</fpage>. doi: <pub-id pub-id-type="doi">10.1200/JCO.2021.39.15_suppl.3564</pub-id>
</citation>
</ref>
<ref id="B34">
<label>34</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pietrantonio</surname> <given-names>F</given-names>
</name>
<name>
<surname>Morano</surname> <given-names>F</given-names>
</name>
<name>
<surname>Lonardi</surname> <given-names>S</given-names>
</name>
<name>
<surname>Raimondi</surname> <given-names>A</given-names>
</name>
<name>
<surname>Salvatore</surname> <given-names>L</given-names>
</name>
<name>
<surname>Marmorino</surname> <given-names>F</given-names>
</name>
<etal/>
</person-group>. <article-title>383O MAYA trial: temozolomide (TMZ) priming followed by combination with low-dose ipilimumab and nivolumab in patients with microsatellite stable (MSS), MGMT silenced metastatic colorectal cancer (mCRC)</article-title>. <source>Ann Oncol</source> (<year>2021</year>) <volume>32</volume>:<page-range>S530&#x2013;1</page-range>. doi: <pub-id pub-id-type="doi">10.1016/j.annonc.2021.08.905</pub-id>
</citation>
</ref>
<ref id="B35">
<label>35</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Panarelli</surname> <given-names>NC</given-names>
</name>
<name>
<surname>Vaughn</surname> <given-names>CP</given-names>
</name>
<name>
<surname>Samowitz</surname> <given-names>WS</given-names>
</name>
<name>
<surname>Yantiss</surname> <given-names>RK</given-names>
</name>
</person-group>. <article-title>Sporadic microsatellite instability-high colon cancers rarely display immunohistochemical evidence of wnt signaling activation</article-title>. <source>Am J Surg Pathol</source> (<year>2015</year>) <volume>39</volume>:<page-range>313&#x2013;7</page-range>. doi: <pub-id pub-id-type="doi">10.1097/PAS.0000000000000380</pub-id>
</citation>
</ref>
<ref id="B36">
<label>36</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<collab>Cancer Genome Atlas Network</collab>
</person-group>. <article-title>Comprehensive molecular characterization of human colon and rectal cancer</article-title>. <source>Nature</source> (<year>2012</year>) <volume>487</volume>:<page-range>330&#x2013;7</page-range>. doi: <pub-id pub-id-type="doi">10.1038/nature11252</pub-id>
</citation>
</ref>
<ref id="B37">
<label>37</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hu</surname> <given-names>G</given-names>
</name>
<name>
<surname>Tu</surname> <given-names>W</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>L</given-names>
</name>
<name>
<surname>Peng</surname> <given-names>G</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>L</given-names>
</name>
</person-group>. <article-title>ARID1A deficiency and immune checkpoint blockade therapy: From mechanisms to clinical application</article-title>. <source>Cancer Lett</source> (<year>2020</year>) <volume>473</volume>:<page-range>148&#x2013;55</page-range>. doi: <pub-id pub-id-type="doi">10.1016/j.canlet.2020.01.001</pub-id>
</citation>
</ref>
<ref id="B38">
<label>38</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kassiotis</surname> <given-names>G</given-names>
</name>
</person-group>. <article-title>Endogenous retroviruses and the development of cancer</article-title>. <source>J Immunol</source> (<year>2014</year>) <volume>192</volume>:<page-range>1343&#x2013;9</page-range>. doi: <pub-id pub-id-type="doi">10.4049/jimmunol.1302972</pub-id>
</citation>
</ref>
<ref id="B39">
<label>39</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jansz</surname> <given-names>N</given-names>
</name>
<name>
<surname>Faulkner</surname> <given-names>GJ</given-names>
</name>
</person-group>. <article-title>Endogenous retroviruses in the origins and treatment of cancer</article-title>. <source>Genome Biol</source> (<year>2021</year>) <volume>22</volume>:<fpage>147</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s13059-021-02357-4</pub-id>
</citation>
</ref>
<ref id="B40">
<label>40</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chiappinelli</surname> <given-names>KB</given-names>
</name>
<name>
<surname>Strissel</surname> <given-names>PL</given-names>
</name>
<name>
<surname>Desrichard</surname> <given-names>A</given-names>
</name>
<name>
<surname>Li</surname> <given-names>H</given-names>
</name>
<name>
<surname>Henke</surname> <given-names>C</given-names>
</name>
<name>
<surname>Akman</surname> <given-names>B</given-names>
</name>
<etal/>
</person-group>. <article-title>Inhibiting DNA methylation causes an interferon response in cancer <italic>via</italic> dsrna including endogenous retroviruses</article-title>. <source>Cell</source> (<year>2015</year>) <volume>162</volume>:<page-range>974&#x2013;86</page-range>. doi: <pub-id pub-id-type="doi">10.1016/j.cell.2015.07.011</pub-id>
</citation>
</ref>
<ref id="B41">
<label>41</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Roulois</surname> <given-names>D</given-names>
</name>
<name>
<surname>Loo Yau</surname> <given-names>H</given-names>
</name>
<name>
<surname>Singhania</surname> <given-names>R</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Danesh</surname> <given-names>A</given-names>
</name>
<name>
<surname>Shen</surname> <given-names>SY</given-names>
</name>
<etal/>
</person-group>. <article-title>DNA-Demethylating agents target colorectal cancer cells by inducing viral mimicry by endogenous transcripts</article-title>. <source>Cell</source> (<year>2015</year>) <volume>162</volume>:<page-range>961&#x2013;73</page-range>. doi: <pub-id pub-id-type="doi">10.1016/j.cell.2015.07.056</pub-id>
</citation>
</ref>
<ref id="B42">
<label>42</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kaur</surname> <given-names>J</given-names>
</name>
<name>
<surname>Daoud</surname> <given-names>A</given-names>
</name>
<name>
<surname>Eblen</surname> <given-names>ST</given-names>
</name>
</person-group>. <article-title>Targeting chromatin remodeling for cancer therapy</article-title>. <source>Curr Mol Pharmacol</source> (<year>2019</year>) <volume>12</volume>:<page-range>215&#x2013;29</page-range>. doi: <pub-id pub-id-type="doi">10.2174/1874467212666190215112915</pub-id>
</citation>
</ref>
<ref id="B43">
<label>43</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yuki</surname> <given-names>S</given-names>
</name>
<name>
<surname>Bando</surname> <given-names>H</given-names>
</name>
<name>
<surname>Tsukada</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Inamori</surname> <given-names>K</given-names>
</name>
<name>
<surname>Komatsu</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Homma</surname> <given-names>S</given-names>
</name>
<etal/>
</person-group>. <article-title>Short-term results of VOLTAGE-a: nivolumab monotherapy and subsequent radical surgery following preoperative chemoradiotherapy in patients with microsatellite stable and microsatellite instability-high locally advanced rectal cancer</article-title>. <source>J Clin Oncol</source> (<year>2020</year>) <volume>38</volume>:<fpage>4100</fpage>. doi: <pub-id pub-id-type="doi">10.1200/JCO.2020.38.15_suppl.4100</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>