<?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. Oncol.</journal-id>
<journal-title>Frontiers in Oncology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Oncol.</abbrev-journal-title>
<issn pub-type="epub">2234-943X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fonc.2022.861392</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Oncology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Development and Validation of an Immune-Related Prognostic Signature in Cervical Cancer</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Su</surname><given-names>Rongjia</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1473780"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Jin</surname><given-names>Chengwen</given-names>
</name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1150973"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Bu</surname><given-names>Hualei</given-names>
</name>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1676237"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Xiang</surname><given-names>Jiangdong</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1205739"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Zhou</surname><given-names>Lina</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn001"><sup>*</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1631955"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Jin</surname><given-names>Chengjuan</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn001"><sup>*</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1152865"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Department of Obstetrics and Gynecology, Shanghai General Hospital, Shanghai Jiaotong University School of Medicine</institution>, <addr-line>Shanghai</addr-line>, <country>China</country></aff>
<aff id="aff2"><sup>2</sup><institution>Department of Gynecologic Oncology, International Peace Maternity and Child Health Hospital, Shanghai Jiaotong University School of Medicine</institution>, <addr-line>Shanghai</addr-line>, <country>China</country></aff>
<aff id="aff3"><sup>3</sup><institution>Department of Central Laboratory, Shanghai Chest Hospital, Shanghai Jiaotong University School of Medicine</institution>, <addr-line>Shanghai</addr-line>, <country>China</country></aff>
<aff id="aff4"><sup>4</sup><institution>Department of Obstetrics and Gynecology, Qilu Hospital, Shandong University</institution>, <addr-line>Jinan</addr-line>, <country>China</country></aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Florentia Fostira, National Centre of Scientific Research Demokritos, Greece</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Sarah Taylor, University of Pittsburgh, United States; Paraskevi Apostolou, National Centre of Scientific Research Demokritos, Greece</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Chengjuan Jin, <email xlink:href="mailto:chengjuanjin@126.com">chengjuanjin@126.com</email>; Lina Zhou, <email xlink:href="mailto:nanazhou1213@163.com">nanazhou1213@163.com</email>
</p>
</fn>
<fn fn-type="other" id="fn002">
<p>This article was submitted to Gynecological Oncology, a section of the journal Frontiers in Oncology</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>16</day>
<month>05</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>12</volume>
<elocation-id>861392</elocation-id>
<history>
<date date-type="received">
<day>24</day>
<month>01</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>08</day>
<month>04</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2022 Su, Jin, Bu, Xiang, Zhou and Jin</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Su, Jin, Bu, Xiang, Zhou and Jin</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>Cervical cancer is the fourth most frequent gynecological malignancy across the world. Immunotherapies have proved to improve prognosis of cervical cancer. However, few studies on immune-related prognostic signature had been reported in cervical cancer.</p>
</sec>
<sec>
<title>Methods</title>
<p>Raw data and clinical information of cervical cancer samples were downloaded from TCGA and UCSC Xena website. Immunophenoscore of immune infiltration cells in cervical cancer samples was calculated through the ssGSEA method using GSVA package. WGCNA, Cox regression analysis, LASSO analysis, and GSEA analysis were performed to classify cervical cancer prognosis and explore the biological signaling pathway.</p>
</sec>
<sec>
<title>Results</title>
<p>There were eight immune infiltration cells associated with prognosis of cervical cancer. Through WGCNA, 153 genes from 402 immune-related genes were significantly correlated with prognosis of cervical cancer. A 15-gene signature demonstrated powerful predictive ability in prognosis of cervical cancer. GSEA analysis showed multiple signaling pathways containing Programmed cell death ligand-1 (PD-L1) expression and PD-1 checkpoint pathway differences between high-risk and low-risk groups. Furthermore, the 15-gene signature was associated with multiple immune cells and immune infiltration in tumor microenvironment.</p>
</sec>
<sec>
<title>Conclusion</title>
<p>The 15-gene signature is an effective potential prognostic classifier in the immunotherapies and surveillance of cervical cancer.</p>
</sec>
</abstract>
<kwd-group>
<kwd>cervical cancer</kwd>
<kwd>weighted gene co-expression network analysis (WGCNA)</kwd>
<kwd>least absolute shrinkage and selection operator (LASSO)</kwd>
<kwd>prognostic signature</kwd>
<kwd>immune signature</kwd>
<kwd>gene set enrichment analysis (GSEA)</kwd>
</kwd-group>    <contract-num rid="cn001"> 82002736 </contract-num>    <contract-sponsor id="cn001">National Natural Science Foundation of China<named-content content-type="fundref-id">10.13039/501100001809</named-content>
</contract-sponsor>
<counts>
<fig-count count="9"/>
<table-count count="3"/>
<equation-count count="0"/>
<ref-count count="52"/>
<page-count count="14"/>
<word-count count="5112"/>
</counts>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<title>Introduction</title>
<p>According to estimates from GLOBOCAN 2018, cervical cancer was the fourth most common cancer among women, with approximately 570,000 new cases and 311,000 deaths (<xref ref-type="bibr" rid="B1">1</xref>). Cure rate for earliest stages is more than 90%, whereas locally advanced lesions treated with multimodality therapy can only achieve 65% cure rate in stage II lesions and 55% in stage III lesions (<xref ref-type="bibr" rid="B2">2</xref>). Therapies to improve survivorship are in desperate need for cervical cancer especially locally advanced disease.</p>
<p>With rapid development of precision medicines, novel therapeutic strategies, especially immunotherapies, have been proposed to significantly improve clinical outcomes of cervical cancer (<xref ref-type="bibr" rid="B3">3</xref>, <xref ref-type="bibr" rid="B4">4</xref>). Harnessing an antitumor immune response has been a fundamental strategy in cancer immunotherapy. A paradigm shift has appeared in cancer immunotherapy: from traditional immune enhancement with low-objective responses and frequent adverse events to more effective and less toxic reactions immune normalization (<xref ref-type="bibr" rid="B5">5</xref>, <xref ref-type="bibr" rid="B6">6</xref>).</p>
<p>Cancer immunotherapies, such as administration of the cytokine Interleukin-2 (IL-2), adoptive cell transfer, and the checkpoint modulators CTLA-4 and PD-1, have proved effective in clinical practice (<xref ref-type="bibr" rid="B6">6</xref>). Blockade of the checkpoint modulators cytotoxic T-lymphocyte-associated antigen 4 (CTLA-4) and programmed cell death protein-1 (PD-1) starts the field of immune normalization in immunotherapy. Upregulated PD-1 in tumor microenvironment inhibits an effector T cell antitumor immune response, and therapies blocking this pathway have proven effective against multiple tumor types (<xref ref-type="bibr" rid="B5">5</xref>). Anti-PD therapy perform antitumor immunity mainly through the following three principles:1) targeting a tumor-induced immune escape mechanism, 2) selectively modulating immunity in the tumor microenvironment, and 3) resetting immunity in the tumor microenvironment (<xref ref-type="bibr" rid="B5">5</xref>). It has been proved that patients with high mutation burden and burden of potential neoepitopes benefit more from immunological checkpoint blockade (<xref ref-type="bibr" rid="B7">7</xref>&#x2013;<xref ref-type="bibr" rid="B11">11</xref>). Immunotherapy plays a dispensable role in management of cervical cancer. In KEYNOTE-158, pembrolizumab has been approved by the US Food and Drug Administration for use in advanced cervical cancer with progressive disease either during or after chemotherapy (<xref ref-type="bibr" rid="B12">12</xref>). The objective response rate was 26.3% with a disease control rate of 68.4% for immunotherapy in cervical cancer (<xref ref-type="bibr" rid="B4">4</xref>). However, only small proportion patients benefited from immunotherapy, and this proportion will hopefully increase with better patient selection and combination therapy. Hence, it is urgent to find potential biomarkers for prediction of response to checkpoint immunotherapy and the rationale for the use of checkpoint immunotherapy.</p>
<p>In the present study, we qualify immune cells infiltration in cervical cancer and analyze the correlation between immune cells and cancer prognosis. Hub genes regulating prognosis through immune infiltration in cervical cancer were identified by weighted gene co-expression network analysis (WGCNA) and least absolute shrinkage and selection operator (LASSO). It was suggested by infiltrated immune cells and pathway enrichment analysis that our immune-related signature was closely related to tumor prognosis and could predict response of immunotherapy. A robust immune-related prognostic signature based on transcriptomics in cervical cancer was constructed and validated.</p>
</sec>
<sec id="s2" sec-type="materials|methods">
<title>Materials and Methods</title>
<sec id="s2_1">
<title>Data Source and Processing</title>
<p>The gene expression profiles and clinical information of cervical cancer were downloaded from The Cancer Genome Atlas (TCGA) Genomic Data Commons Data Portal (<uri xlink:href="https://portal.gdc.cancer.gov/">https://portal.gdc.cancer.gov/</uri>). Patients with pathologically confirmed cervical cancer and complete information about transcriptomics overall survivals (OSs) were included in this study. Finally, a total of 304 primary cervical cancer samples, two metastatic cervical cancer samples, and three normal cervix samples from TCGA were analyzed in our study. Details on clinical information of the included samples were summarized in <xref ref-type="table" rid="T1"><bold>Table&#xa0;1</bold></xref>. A workflow of this study was indicated in <xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1</bold></xref>.</p>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>Clinical characteristics of 309 samples from TCGA.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="left"/>
<th valign="top" align="center">Alive (N = 235)</th>
<th valign="top" align="center">Dead (N = 74)</th>
<th valign="top" align="center">Overall (N = 309)</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" colspan="4" align="left" style="background-color:#ffffff"><bold>Age&#xa0;(years)</bold>
</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color:#ffffff">&#x2003;Mean (SD)</td>
<td valign="top" align="center" style="background-color:#ffffff">47.4 (13.2)</td>
<td valign="top" align="center" style="background-color:#ffffff">50.8 (15.4)</td>
<td valign="top" align="center" style="background-color:#ffffff">48.2 (13.8)</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color:#ffffff">&#x2003;Median [Min, Max]</td>
<td valign="top" align="center" style="background-color:#ffffff">46.0 [20.0, 88.0]</td>
<td valign="top" align="center" style="background-color:#ffffff">48.0 [21.0, 79.0]</td>
<td valign="top" align="center" style="background-color:#ffffff">46.0 [20.0, 88.0]</td>
</tr>
<tr>
<td valign="top" colspan="4" align="left" style="background-color:#ffffff"><bold>Pathology</bold>
</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color:#ffffff">&#x2003;Adenocarcinoma</td>
<td valign="top" align="center" style="background-color:#ffffff">24 (10.2%)</td>
<td valign="top" align="center" style="background-color:#ffffff">5 (6.8%)</td>
<td valign="top" align="center" style="background-color:#ffffff">29 (9.4%)</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color:#ffffff">&#x2003;Adenosquamous carcinoma</td>
<td valign="top" align="center" style="background-color:#ffffff">4 (1.7%)</td>
<td valign="top" align="center" style="background-color:#ffffff">1 (1.4%)</td>
<td valign="top" align="center" style="background-color:#ffffff">5 (1.6%)</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color:#ffffff">&#x2003;Basaloid squamous cell carcinoma</td>
<td valign="top" align="center" style="background-color:#ffffff">1 (0.4%)</td>
<td valign="top" align="center" style="background-color:#ffffff">0 (0%)</td>
<td valign="top" align="center" style="background-color:#ffffff">1 (0.3%)</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color:#ffffff">&#x2003;Endometrioid adenocarcinoma</td>
<td valign="top" align="center" style="background-color:#ffffff">3 (1.3%)</td>
<td valign="top" align="center" style="background-color:#ffffff">0 (0%)</td>
<td valign="top" align="center" style="background-color:#ffffff">3 (1.0%)</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color:#ffffff">&#x2003;Mucinous adenocarcinoma</td>
<td valign="top" align="center" style="background-color:#ffffff">12 (5.1%)</td>
<td valign="top" align="center" style="background-color:#ffffff">5 (6.8%)</td>
<td valign="top" align="center" style="background-color:#ffffff">17 (5.5%)</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color:#ffffff">&#x2003;Papillary squamous cell carcinoma</td>
<td valign="top" align="center" style="background-color:#ffffff">1 (0.4%)</td>
<td valign="top" align="center" style="background-color:#ffffff">0 (0%)</td>
<td valign="top" align="center" style="background-color:#ffffff">1 (0.3%)</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color:#ffffff">&#x2003;Squamous cell carcinoma</td>
<td valign="top" align="center" style="background-color:#ffffff">190 (80.9%)</td>
<td valign="top" align="center" style="background-color:#ffffff">63 (85.1%)</td>
<td valign="top" align="center" style="background-color:#ffffff">253 (81.9%)</td>
</tr>
<tr>
<td valign="top" colspan="4" align="left" style="background-color:#ffffff"><bold>Stage</bold>
</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color:#ffffff">&#x2003;I</td>
<td valign="top" align="center" style="background-color:#ffffff">125 (53.2%)</td>
<td valign="top" align="center" style="background-color:#ffffff">38 (51.4%)</td>
<td valign="top" align="center" style="background-color:#ffffff">163 (52.8%)</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color:#ffffff">&#x2003;II</td>
<td valign="top" align="center" style="background-color:#ffffff">60 (25.5%)</td>
<td valign="top" align="center" style="background-color:#ffffff">11 (14.9%)</td>
<td valign="top" align="center" style="background-color:#ffffff">71 (23.0%)</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color:#ffffff">&#x2003;III</td>
<td valign="top" align="center" style="background-color:#ffffff">35 (14.9%)</td>
<td valign="top" align="center" style="background-color:#ffffff">11 (14.9%)</td>
<td valign="top" align="center" style="background-color:#ffffff">46 (14.9%)</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color:#ffffff">&#x2003;IV</td>
<td valign="top" align="center" style="background-color:#ffffff">8 (3.4%)</td>
<td valign="top" align="center" style="background-color:#ffffff">14 (18.9%)</td>
<td valign="top" align="center" style="background-color:#ffffff">22 (7.1%)</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color:#ffffff">&#x2003;Missing</td>
<td valign="top" align="center" style="background-color:#ffffff">7 (3.0%)</td>
<td valign="top" align="center" style="background-color:#ffffff">0 (0%)</td>
<td valign="top" align="center" style="background-color:#ffffff">7 (2.3%)</td>
</tr>
<tr>
<td valign="top" colspan="4" align="left" style="background-color:#ffffff"><bold>Sample</bold>
</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color:#ffffff">&#x2003;Metastatic</td>
<td valign="top" align="center" style="background-color:#ffffff">1 (0.4%)</td>
<td valign="top" align="center" style="background-color:#ffffff">1 (1.4%)</td>
<td valign="top" align="center" style="background-color:#ffffff">2 (0.6%)</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color:#ffffff">&#x2003;Primary tumor</td>
<td valign="top" align="center" style="background-color:#ffffff">232 (98.7%)</td>
<td valign="top" align="center" style="background-color:#ffffff">72 (97.3%)</td>
<td valign="top" align="center" style="background-color:#ffffff">304 (98.4%)</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color:#ffffff">&#x2003;Solid tissue normal</td>
<td valign="top" align="center" style="background-color:#ffffff">2 (0.9%)</td>
<td valign="top" align="center" style="background-color:#ffffff">1 (1.4%)</td>
<td valign="top" align="center" style="background-color:#ffffff">3 (1.0%)</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Flowchart of the 15-gene signature prognostic model construction process.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-12-861392-g001.tif"/>
</fig>
</sec>
<sec id="s2_2">
<title>Infiltration of Immune Cells</title>
<p>Tumor-infiltrating immune cells can be quantified from RNA sequencing (RNA-seq) data of human tumors using bioinformatics approaches. Single-sample gene set enrichment analysis (ssGSEA) calculated and qualified the infiltration of immune cells through RNA-seq data. The infiltration of 28 immune cells was obtained. Immunophenoscore (normalized enrichment score) of each TCGA cervical cancer sample was calculated through ssGSEA. The ssGSEA ranks the genes by their absolute expression in a sample and computes enrichment score by integrating the differences between the empirical cumulative distribution functions of the gene ranks (<xref ref-type="bibr" rid="B13">13</xref>).</p>
</sec>
<sec id="s2_3">
<title>Survival Analysis</title>
<p>Univariate Cox regression analysis was performed to identify the association between cervical cancer survival and immunophenoscore of infiltrated immune cells. Forest plot was drawn to demonstrate the influence of infiltrated immune cells on survival. The best separation statistic was performed by use of &#x201c;survminer&#x201d; package that divides gene expression into high groups and low groups based on best separation. Next, Kaplan&#x2013;Meier curve was made to further analyze relationship between survival and infiltration immune cells.</p>
</sec>
<sec id="s2_4">
<title>Construction of Weighted Gene Co-Expression Network Analysis</title>
<p>WGCNA is an algorithm for finding genetic interactions in a weighted manner. It is used to build a gene co-expression networks to mine network modules closely associated with clinical traits through systematic biological method (<xref ref-type="bibr" rid="B14">14</xref>). In this study, immunophenoscore of the infiltrated immune cells was regarded as target clinical traits. Genes expressing not available (NA) were removed. The top 25% genes with most median absolute deviation used as a robust measure of variability were selected for WGCNA analysis (<xref ref-type="bibr" rid="B15">15</xref>, <xref ref-type="bibr" rid="B16">16</xref>).</p>
<p>Using the WGCNA function adjacency, an adjacency matrix is constructed by computing the Pearson correlation between all pairs of genes in the selected sample. Genes were divided into different gene modules based on the dissimilarity measure. A hierarchical clustering tree was constructed with different branches of the tree representing different gene modules (<xref ref-type="bibr" rid="B15">15</xref>&#x2013;<xref ref-type="bibr" rid="B17">17</xref>).</p>
<p>The WGCNA network was built. There were eight immune cells associated with cervical cancer survival. The magenta module was the module that most related to infiltration of immune cells. A total of 402 genes contained in magenta module were analyzed for survival using the standardized expression data FPKM. Finally, we found that 153 genes were associated with survival (p &lt; 0.05) in magenta module.</p>
</sec>
<sec id="s2_5">
<title>Screening Hub Genes by LASSO and Survival Analysis</title>
<p>The 153 genes with the highest correlation associated with survival were picked out from WGCNA analysis. Then, hub genes were further screened from the 153 genes by the use of LASSO. Survival analysis of hub genes was performed using &#x201c;survival&#x201d; and &#x201c;survminer&#x201d; packages to verify whether differently expressed genes affected tumor prognosis. Data of cervical cancer from University of California, Santa Cruz (UCSC) Xena were downloaded and applied to compare the expression of hub genes in normal cervix and cervical cancer.</p>
</sec>
<sec id="s2_6">
<title>Construction of Prognostic Scoring Model</title>
<p>Patients with cervical cancer in TCGA were randomly divided into two groups by stage: the training group of 70% and the test group of 30% by the use of caret R package. Finally, a total of 199 patients were used as training group and 74 patients were regarded as test group. The selected key genes for support vector machine analysis were used to fit a LASSO Cox&#x2013;proportional hazards (Cox-PH) model for selecting an optimal panel of predictive genes with penalized package of R (<xref ref-type="bibr" rid="B18">18</xref>). Optimal lambda value was computed through a 10 cross-validations. Next, Cox-PH coefficients and infiltration of immune cells levels of these selected genes were used to calculate prognostic score as follows: Risk Score = &#x3a3; (coef<sub>gene</sub> &#xd7; immunophenoscore<sub>gene</sub>), where Coef<sub>gene</sub> and immunophenoscore<sub>gene</sub> suggest Cox-PH coefficient and immunophenoscore level of a gene, respectively.</p>
</sec>
<sec id="s2_7">
<title>Validation of Prognostic Scoring Model</title>
<p>We calculated risk score of every patient in the training and test group by the model. We separated the training group and test group into a high-risk group and a low-risk group with the median risk score as cutoff, respectively. Kaplan&#x2013;Meier curve was applied to obtain and compare the OS time of two risk groups. ROC (operating characteristic curve) was performed to evaluate the predictive accuracy of the model.</p>
</sec>
<sec id="s2_8">
<title>Comparison of Eight Immune Cell Subtypes Between High-Risk and Low-Risk Groups</title>
<p>To explore the differences of immune cell subtypes between high-risk and low-risk groups, the eight immune cell subtypes associated with OS in cervical cancer were assessed in test and train cohorts. Mann&#x2013;Whitney U-test was used to compare differences in immune cell subtypes in the high-risk and low-risk groups.</p>
</sec>
<sec id="s2_9">
<title>Gene Set Enrichment Analysis (GSEA)</title>
<p>All 304 patients with primary cervical cancer were evaluated by prognostic scoring model and then divided into high-risk group and low-risk group based on the standard cutoff. The global gene expression was analyzed and displayed with volcanic maps using the limma package in R. GSEA was conducted respectively to search &#x201c;all gene sets&#x201d; enriched in the samples with high-risk group and low-risk group. The differentially expressed genes were enriched in PD-L1 pathway functional sets (<xref ref-type="bibr" rid="B19">19</xref>).</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<title>Results</title>
<sec id="s3_1">
<title>Quantify Immune Cell Infiltration and Survival Analysis</title>
<p>We used ssGSEA to quantify mRNA data for immune cell infiltration. Finally, 28 infiltrating immune cells were included. Immunophenoscores of 304 primary cervical cancer samples, two metastatic cervical cancer samples, and three normal cervix samples for 28 immune cells were calculated and demonstrated in <xref ref-type="fig" rid="f2"><bold>Figure&#xa0;2A</bold></xref>. Univariate Cox regression analysis was performed to identify association between immune cells and OS of cervical cancer, and the results are shown in <xref ref-type="fig" rid="f2"><bold>Figure&#xa0;2B</bold></xref>. It indicated that eight immune cells were correlated with OS significantly. Moreover, these eight immune cells (activated B cell, activated CD8 T cell, eosinophil, monocyte, activated CD4 T cell, effector memory CD8 T cell, immature B cell, and plasmacytoid dendritic cell) were all protective factors for OS. Kaplan&#x2013;Meier curve validated that high expression of these eight immune cells were related to longer OS time suggested worse prognosis (<xref ref-type="fig" rid="f3"><bold>Figure&#xa0;3</bold></xref>).</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Immunophenoscores and hazard ratios of 28 immune cells in cervical cancer. <bold>(A)</bold> Boxplot of immunophenoscores from 28 immune cells TCGA cervical cancer and normal samples. <bold>(B)</bold> Forest plot of hazard ratios in each infiltrated immune cell for OS. Infiltrated immune cell subtypes significantly associated with overall survival were underlined with a red dashed line. *: P&lt;0.05, **: P&lt;0.01, ns, no significance.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-12-861392-g002.tif"/>
</fig>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Kaplan&#x2013;Meier curves for overall survival of eight infiltrated immune cells. Activated B cell, activated CD8 T cell, eosinophil, monocyte, activated CD4 T cell, effector memory CD8 T cell, immature B cell, and plasmacytoid dendritic cell were associated with overall survival of cervical cancer.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-12-861392-g003.tif"/>
</fig>
</sec>
<sec id="s3_2">
<title>The Weighted Gene Co-Expression Network Analysis Construction and Key Module Identification</title>
<p>The top 25% of the gene expression of variance was screened by the quartile of gene expression level, and 609 genes were screened out to construct co-expressed gene networks and the sample dendrogram and trait heatmap are constructed (<xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4A</bold></xref>). In this study, the power of &#x3b2; = 8 was selected to ensure a scale-free network (<xref ref-type="fig" rid="f4"><bold>Figures&#xa0;4B, C</bold></xref>) (scale-free R<sup>2</sup> = 0.9, slope = &#x2212;1.65).</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Sample tree and trait heatmap of 28 immune cells. <bold>(A)</bold> Sample dendrogram and trait heatmap of 28 immune cells. The clinical trait information is immunophenoscores of 28 infiltrated immune cells in each sample. <bold>(B)</bold> Analysis of the scale-free fit index for various soft-thresholding powers (&#x3b2;) and the mean connectivity for soft threshold powers. The soft-thresholding power in the WGCNA was determined on the basis of a scale-free R<sup>2</sup> (R<sup>2</sup> = 0.90). The left panel presents the relationship between the soft-threshold and scale-free R<sup>2</sup>. The right panel presents the relationship between the soft-threshold and mean connectivity. <bold>(C)</bold> Histogram of connectivity distribution when &#x3b2; = 8 and checking the scale-free topology when &#x3b2; = 8.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-12-861392-g004.tif"/>
</fig>
<p>After merging the modules with the high similarity of feature genes in the gene cluster dendrogram through a cutline (0.25) (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5A</bold></xref>), nine modules were identified by the dynamic tree cut method. The clustering dendrograms of genes were shown in <xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5B</bold></xref>. A heatmap illustrating the correlation between immunophenoscores of infiltrated immune cells and key genes in the module was created (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5C</bold></xref>).</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Weighted gene co-expression network in cervical cancer. <bold>(A)</bold> Clustering of module eigengenes. A cutline (0.25) was selected for module dendrogram and merged some modules according to dissimilarity of estimated module eigengenes. <bold>(B)</bold> A dendrogram of the differentially expressed genes clustered based on different metrics. Each branch in the figure represents one gene, and every color below represents one co-expression module. <bold>(C)</bold> A heatmap showing the correlation between the gene module and clinical traits. The correlation coefficient in each cell represented the correlation between gene module and the clinical traits. The magenta module showed the highest correlated with the infiltrated immune cells.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-12-861392-g005.tif"/>
</fig>
<p>It was obviously that magenta module was the most correlated with immune infiltration in the heatmap. Hence, magenta module was selected as the clinical significant module for further analysis. We found that the correlation coefficient between magenta module and 11 immune infiltration cells was &#x2265;&#xa0;0.7, which suggested strong correlation. This indicates that the genes in this module are most relevant to tumor OS. All 402 genes in magenta module were analyzed for survival of cervical cancer. Through univariate Cox regression survival analysis, 153 genes found to be significantly correlated with prognosis of cervical cancer were selected for further analysis.</p>
</sec>
<sec id="s3_3">
<title>Screening the Key Genes and Survival Analysis</title>
<p>To develop a prognostic scoring model based on WGCNA, the 153 key genes were used to fit the LASSO Cox-PH model. Using parameter lambda (0.038) obtained upon performing 1000 cross-validations, a combination of 15 genes was obtained. These 15 genes were as follows: LAG3, CD74, CCL22, CH25H, OLR1, MIAT, BATF, IKZF3, TRARG, ACSL6, C11orf21, GTSF1, APOL1, CD1C, and LINC00158 (<xref ref-type="table" rid="T2"><bold>Table&#xa0;2</bold></xref> and <xref ref-type="fig" rid="f6"><bold>Figure&#xa0;6</bold></xref>). Risk score of each sample in the training set was calculated. Then, the samples in the training cohort were divided into high- and low-risk groups according to median risk score. As demonstrated in <xref ref-type="fig" rid="f6"><bold>Figure&#xa0;6</bold></xref>, the high-risk group patients in the training cohort have shorter OS time than the low-risk group (p &lt; 0.001, HR: 5.1, 95% CI: 3.4&#x2013;7.5), with an AUC of 0.803 in 1 year, 0.809 in 3 years, and 0.800 in 5 years. Next, the 15-gene signature was then verified by the validation cohort from TCGA. Consistently, the high-risk group patients in the validation cohort had worse prognosis than the low-risk group (p &lt; 0.001, HR: 1.3, 95% CI: 1.1&#x2013;1.4), and the AUC was 0.845 in 1 year, 0.703 in 3 years, and 0.761 in 5 years (<xref ref-type="fig" rid="f7"><bold>Figure&#xa0;7</bold></xref>). All these suggested that the 15-gene signature associated with immune infiltration is able to predict prognosis in patients with cervical cancer.</p>
<table-wrap id="T2" position="float">
<label>Table&#xa0;2</label>
<caption>
<p>Fifteen genes from LASSO analysis.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="left">Gene </th>
<th valign="top" align="center">HR</th>
<th valign="top" align="center">95% CI</th>
<th valign="top" align="center">P-value</th>
<th valign="top" align="center">Lasso Coefficient</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">TRARG1</td>
<td valign="top" align="center">660</td>
<td valign="top" align="center">(0.00027&#x2013;1.6e+09)</td>
<td valign="top" align="center">0.39</td>
<td valign="top" align="center">4.7797999</td>
</tr>
<tr>
<td valign="top" align="left">ACSL6</td>
<td valign="top" align="center">3.4</td>
<td valign="top" align="center">(1.4&#x2013;8.6)</td>
<td valign="top" align="center">0.008</td>
<td valign="top" align="center">0.653170468</td>
</tr>
<tr>
<td valign="top" align="left">C11orf21</td>
<td valign="top" align="center">1.3</td>
<td valign="top" align="center">(1.1&#x2013;1.5)</td>
<td valign="top" align="center">0.00022</td>
<td valign="top" align="center">0.128185283</td>
</tr>
<tr>
<td valign="top" align="left">GTSF1</td>
<td valign="top" align="center">1.1</td>
<td valign="top" align="center">(1.0&#x2013;1.2)</td>
<td valign="top" align="center">0.002</td>
<td valign="top" align="center">0.081672124</td>
</tr>
<tr>
<td valign="top" align="left">CD74</td>
<td valign="top" align="center">1.0</td>
<td valign="top" align="center">(1.0&#x2013;1.0)</td>
<td valign="top" align="center">0.03</td>
<td valign="top" align="center">&#x2212;8.18E&#x2212;05</td>
</tr>
<tr>
<td valign="top" align="left">APOL1</td>
<td valign="top" align="center">1.0</td>
<td valign="top" align="center">(1.0&#x2013;1.0)</td>
<td valign="top" align="center">0.016</td>
<td valign="top" align="center">0.001359139</td>
</tr>
<tr>
<td valign="top" align="left">OLR1</td>
<td valign="top" align="center">1.0</td>
<td valign="top" align="center">(1.0&#x2013;1.0)</td>
<td valign="top" align="center">0.02</td>
<td valign="top" align="center">0.008057682</td>
</tr>
<tr>
<td valign="top" align="left">MIAT</td>
<td valign="top" align="center">1.0</td>
<td valign="top" align="center">(1.0&#x2013;1.0)</td>
<td valign="top" align="center">5.00E&#x2212;04</td>
<td valign="top" align="center">0.01207222</td>
</tr>
<tr>
<td valign="top" align="left">BATF</td>
<td valign="top" align="center">0.96</td>
<td valign="top" align="center">(0.93&#x2013;1.0)</td>
<td valign="top" align="center">0.025</td>
<td valign="top" align="center">&#x2212;0.006690799</td>
</tr>
<tr>
<td valign="top" align="left">LAG3</td>
<td valign="top" align="center">0.93</td>
<td valign="top" align="center">(0.86&#x2013;0.99)</td>
<td valign="top" align="center">0.029</td>
<td valign="top" align="center">&#x2212;0.012618</td>
</tr>
<tr>
<td valign="top" align="left">CCL22</td>
<td valign="top" align="center">0.90</td>
<td valign="top" align="center">(0.79&#x2013;1.0)</td>
<td valign="top" align="center">0.073</td>
<td valign="top" align="center">&#x2212;0.003652351</td>
</tr>
<tr>
<td valign="top" align="left">CH25H</td>
<td valign="top" align="center">0.88</td>
<td valign="top" align="center">(0.77&#x2013;1.0)</td>
<td valign="top" align="center">0.053</td>
<td valign="top" align="center">&#x2212;0.017915424</td>
</tr>
<tr>
<td valign="top" align="left">CD1C</td>
<td valign="top" align="center">0.78</td>
<td valign="top" align="center">(0.63&#x2013;0.97)</td>
<td valign="top" align="center">0.027</td>
<td valign="top" align="center">&#x2212;0.100839181</td>
</tr>
<tr>
<td valign="top" align="left">IKZF3</td>
<td valign="top" align="center">0.76</td>
<td valign="top" align="center">(0.63&#x2013;0.93)</td>
<td valign="top" align="center">0.0061</td>
<td valign="top" align="center">&#x2212;0.086985406</td>
</tr>
<tr>
<td valign="top" align="left">LINC00158</td>
<td valign="top" align="center">0.00062</td>
<td valign="top" align="center">(6.7e&#x2212;07&#x2013;0.58)</td>
<td valign="top" align="center">0.034</td>
<td valign="top" align="center">&#x2212;0.708766884</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Construction of the prognostic classifier. <bold>(A, B)</bold> Determination of the number of factors by the LASSO analysis. <bold>(C)</bold> The distribution of risk score in test cohort. <bold>(D)</bold> The survival duration and status of patients in test cohort. <bold>(E)</bold> A heatmap of immune related genes in the classifier in test cohort. <bold>(F)</bold> Kaplan&#x2013;Meier curve for patients with cervical cancer in test cohort. <bold>(G)</bold> ROC curve for cervical cancer in test cohort.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-12-861392-g006.tif"/>
</fig>
<fig id="f7" position="float">
<label>Figure&#xa0;7</label>
<caption>
<p>The distribution of time-dependent ROC curves and Kaplan&#x2013;Meier survival based on the integrated classifier in the train cohort. <bold>(A)</bold> The distribution of risk score in train cohort. <bold>(B)</bold> The survival duration and status of patients in train cohort. <bold>(C)</bold> A heatmap of immune related genes in the classifier in train cohort. <bold>(D)</bold> Kaplan&#x2013;Meier curve for patients with cervical cancer in train cohort. <bold>(E)</bold> ROC curve for cervical cancer in train cohort.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-12-861392-g007.tif"/>
</fig>
</sec>
<sec id="s3_4">
<title>Immune Cell Subtypes Between High-Risk and Low-Risk Groups</title>
<p>The expression levels of the 15 genes in test and validation cohorts are shown in <xref ref-type="fig" rid="f8"><bold>Figure&#xa0;8A</bold></xref>. Different immune scores had differential OS in patients with cervical cancer. In both test cohorts and train cohorts, these eight immune cell subtypes (activated B cell, activated CD8 T cell, eosinophil, monocyte, activated CD4 T cell, effector memory CD8 T cell, immature B cell, and plasmacytoid dendritic cell) expressed differentially in high-risk and low-risk groups. (<xref ref-type="fig" rid="f8"><bold>Figures&#xa0;8B, C</bold></xref>).</p>
<fig id="f8" position="float">
<label>Figure&#xa0;8</label>
<caption>
<p>Expression profile of 15 genes and immunophenoscores of eight immune cell subtypes in test and train sets. <bold>(A)</bold> Expression profile of 15 genes. The test cohort is shown in the left, and the train cohort is shown in the right. <bold>(B, C)</bold> Immunophenoscores of eight infiltrated immune cells in high-risk and low-risk groups based on the integrated classifier. The test cohort is shown in the left, and the train cohort is shown in the right. *: P&lt;0.05, **: P&lt;0.01, ***: P&lt;0.001****: P&lt;0.0001, ns, no significance.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-12-861392-g008.tif"/>
</fig>
</sec>
<sec id="s3_5">
<title>GSEA Analysis</title>
<p>All 304 primary cervical cancer samples were divided into high-risk group and low-risk group based on the 15-gene signature. To&#xa0;identify potential function of the 15 key genes, GSEA was conducted respectively to search &#x201c;all gene sets&#x201d; enriched in all 304 samples. GSEA showed 33 significant Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway associated with risk score, including PD-L1 expression and PD-1 checkpoint pathway in cancer (<xref ref-type="fig" rid="f9"><bold>Figure&#xa0;9</bold></xref> and <xref ref-type="table" rid="T3"><bold>Table&#xa0;3</bold></xref>). These results validated that the key genes in clinical significant module were mainly involved in the regulation of immune system. From GSEA, the 15-gene signature obviously participated in regulation of PD-1/PD-L1 pathway in cervical cancer.</p>
<fig id="f9" position="float">
<label>Figure&#xa0;9</label>
<caption>
<p>GSEA analysis and PD-L1 pathway. <bold>(A)</bold> GSEA analysis demonstrated 33 KEGG pathways associated with risk score in cervical cancer. <bold>(B)</bold> PD-L1 expression and PD-1 checkpoint pathway in cancer was one of the 33 pathways that significantly related with risk score.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-12-861392-g009.tif"/>
</fig>
<table-wrap id="T3" position="float">
<label>Table&#xa0;3</label>
<caption>
<p>Significant KEGG pathways from GSEA analysis.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="left">ID</th>
<th valign="top" align="center">Description</th>
<th valign="top" align="center">P-value</th>
<th valign="top" align="center">p-adjust</th>
<th valign="top" align="center">q-values</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">hsa04060</td>
<td valign="top" align="left">Cytokine-cytokine receptor interaction</td>
<td valign="top" align="center">3.29E&#x2212;05</td>
<td valign="top" align="center">0.000668</td>
<td valign="top" align="center">0.000554</td>
</tr>
<tr>
<td valign="top" align="left">hsa04062</td>
<td valign="top" align="left">Chemokine signaling pathway</td>
<td valign="top" align="center">3.40E&#x2212;05</td>
<td valign="top" align="center">0.000668</td>
<td valign="top" align="center">0.000554</td>
</tr>
<tr>
<td valign="top" align="left">hsa05152</td>
<td valign="top" align="left">Tuberculosis</td>
<td valign="top" align="center">3.42E&#x2212;05</td>
<td valign="top" align="center">0.000668</td>
<td valign="top" align="center">0.000554</td>
</tr>
<tr>
<td valign="top" align="left">hsa04514</td>
<td valign="top" align="left">Cell adhesion molecules</td>
<td valign="top" align="center">3.45E&#x2212;05</td>
<td valign="top" align="center">0.000668</td>
<td valign="top" align="center">0.000554</td>
</tr>
<tr>
<td valign="top" align="left">hsa05322</td>
<td valign="top" align="left">Systemic lupus erythematosus</td>
<td valign="top" align="center">3.49E&#x2212;05</td>
<td valign="top" align="center">0.000668</td>
<td valign="top" align="center">0.000554</td>
</tr>
<tr>
<td valign="top" align="left">hsa04659</td>
<td valign="top" align="left">Th17 cell differentiation</td>
<td valign="top" align="center">3.52E&#x2212;05</td>
<td valign="top" align="center">0.000668</td>
<td valign="top" align="center">0.000554</td>
</tr>
<tr>
<td valign="top" align="left">hsa04061</td>
<td valign="top" align="left">Viral protein interaction with cytokine and cytokine receptor</td>
<td valign="top" align="center">3.53E&#x2212;05</td>
<td valign="top" align="center">0.000668</td>
<td valign="top" align="center">0.000554</td>
</tr>
<tr>
<td valign="top" align="left">hsa04660</td>
<td valign="top" align="left">T cell receptor signaling pathway</td>
<td valign="top" align="center">3.53E&#x2212;05</td>
<td valign="top" align="center">0.000668</td>
<td valign="top" align="center">0.000554</td>
</tr>
<tr>
<td valign="top" align="left">hsa04640</td>
<td valign="top" align="left">Hematopoietic cell lineage</td>
<td valign="top" align="center">3.53E&#x2212;05</td>
<td valign="top" align="center">0.000668</td>
<td valign="top" align="center">0.000554</td>
</tr>
<tr>
<td valign="top" align="left">hsa04658</td>
<td valign="top" align="left">Th1 and Th2 cell differentiation</td>
<td valign="top" align="center">3.54E&#x2212;05</td>
<td valign="top" align="center">0.000668</td>
<td valign="top" align="center">0.000554</td>
</tr>
<tr>
<td valign="top" align="left">hsa04662</td>
<td valign="top" align="left">B cell receptor signaling pathway</td>
<td valign="top" align="center">3.56E&#x2212;05</td>
<td valign="top" align="center">0.000668</td>
<td valign="top" align="center">0.000554</td>
</tr>
<tr>
<td valign="top" align="left">hsa05416</td>
<td valign="top" align="left">Viral myocarditis</td>
<td valign="top" align="center">3.63E&#x2212;05</td>
<td valign="top" align="center">0.000668</td>
<td valign="top" align="center">0.000554</td>
</tr>
<tr>
<td valign="top" align="left">hsa05320</td>
<td valign="top" align="left">Autoimmune thyroid disease</td>
<td valign="top" align="center">3.64E&#x2212;05</td>
<td valign="top" align="center">0.000668</td>
<td valign="top" align="center">0.000554</td>
</tr>
<tr>
<td valign="top" align="left">hsa04672</td>
<td valign="top" align="left">Intestinal immune network for IgA production</td>
<td valign="top" align="center">3.64E&#x2212;05</td>
<td valign="top" align="center">0.000668</td>
<td valign="top" align="center">0.000554</td>
</tr>
<tr>
<td valign="top" align="left">hsa05340</td>
<td valign="top" align="left">Primary immunodeficiency</td>
<td valign="top" align="center">3.67E&#x2212;05</td>
<td valign="top" align="center">0.000668</td>
<td valign="top" align="center">0.000554</td>
</tr>
<tr>
<td valign="top" align="left">hsa05330</td>
<td valign="top" align="left">Allograft rejection</td>
<td valign="top" align="center">3.67E&#x2212;05</td>
<td valign="top" align="center">0.000668</td>
<td valign="top" align="center">0.000554</td>
</tr>
<tr>
<td valign="top" align="left">hsa05310</td>
<td valign="top" align="left">Asthma</td>
<td valign="top" align="center">3.70E&#x2212;05</td>
<td valign="top" align="center">0.000668</td>
<td valign="top" align="center">0.000554</td>
</tr>
<tr>
<td valign="top" align="left">hsa04630</td>
<td valign="top" align="left">JAK-STAT signaling pathway</td>
<td valign="top" align="center">6.88E&#x2212;05</td>
<td valign="top" align="center">0.001174</td>
<td valign="top" align="center">0.000974</td>
</tr>
<tr>
<td valign="top" align="left">hsa05145</td>
<td valign="top" align="left">Toxoplasmosis</td>
<td valign="top" align="center">0.000526</td>
<td valign="top" align="center">0.008506</td>
<td valign="top" align="center">0.007058</td>
</tr>
<tr>
<td valign="top" align="left">hsa04612</td>
<td valign="top" align="left">Antigen processing and presentation</td>
<td valign="top" align="center">0.000647</td>
<td valign="top" align="center">0.009924</td>
<td valign="top" align="center">0.008234</td>
</tr>
<tr>
<td valign="top" align="left">hsa05140</td>
<td valign="top" align="left">Leishmaniasis</td>
<td valign="top" align="center">0.000717</td>
<td valign="top" align="center">0.010486</td>
<td valign="top" align="center">0.008701</td>
</tr>
<tr>
<td valign="top" align="left">hsa05321</td>
<td valign="top" align="left">Inflammatory bowel disease</td>
<td valign="top" align="center">0.001157</td>
<td valign="top" align="center">0.01606</td>
<td valign="top" align="center">0.013326</td>
</tr>
<tr>
<td valign="top" align="left">hsa05235</td>
<td valign="top" align="left">PD-L1 expression and PD-1 checkpoint pathway in cancer</td>
<td valign="top" align="center">0.001203</td>
<td valign="top" align="center">0.01606</td>
<td valign="top" align="center">0.013326</td>
</tr>
<tr>
<td valign="top" align="left">hsa05169</td>
<td valign="top" align="left">Epstein-Barr virus infection</td>
<td valign="top" align="center">0.001524</td>
<td valign="top" align="center">0.01949</td>
<td valign="top" align="center">0.016172</td>
</tr>
<tr>
<td valign="top" align="left">hsa05170</td>
<td valign="top" align="left">Human immunodeficiency virus 1 infection</td>
<td valign="top" align="center">0.001867</td>
<td valign="top" align="center">0.022921</td>
<td valign="top" align="center">0.019019</td>
</tr>
<tr>
<td valign="top" align="left">hsa04950</td>
<td valign="top" align="left">Maturity onset diabetes of the young</td>
<td valign="top" align="center">0.002213</td>
<td valign="top" align="center">0.026135</td>
<td valign="top" align="center">0.021685</td>
</tr>
<tr>
<td valign="top" align="left">hsa03320</td>
<td valign="top" align="left">PPAR signaling pathway</td>
<td valign="top" align="center">0.002867</td>
<td valign="top" align="center">0.032598</td>
<td valign="top" align="center">0.027048</td>
</tr>
<tr>
<td valign="top" align="left">hsa04390</td>
<td valign="top" align="left">Hippo signaling pathway</td>
<td valign="top" align="center">0.003763</td>
<td valign="top" align="center">0.040506</td>
<td valign="top" align="center">0.03361</td>
</tr>
<tr>
<td valign="top" align="left">hsa03010</td>
<td valign="top" align="left">Ribosome</td>
<td valign="top" align="center">0.003912</td>
<td valign="top" align="center">0.040506</td>
<td valign="top" align="center">0.03361</td>
</tr>
<tr>
<td valign="top" align="left">hsa04080</td>
<td valign="top" align="left">Neuroactive ligand-receptor interaction</td>
<td valign="top" align="center">0.004057</td>
<td valign="top" align="center">0.040506</td>
<td valign="top" align="center">0.03361</td>
</tr>
<tr>
<td valign="top" align="left">hsa04020</td>
<td valign="top" align="left">Calcium signaling pathway</td>
<td valign="top" align="center">0.00416</td>
<td valign="top" align="center">0.040506</td>
<td valign="top" align="center">0.03361</td>
</tr>
<tr>
<td valign="top" align="left">hsa05166</td>
<td valign="top" align="left">Human T-cell leukemia virus 1 infection</td>
<td valign="top" align="center">0.004222</td>
<td valign="top" align="center">0.040506</td>
<td valign="top" align="center">0.03361</td>
</tr>
<tr>
<td valign="top" align="left">hsa04940</td>
<td valign="top" align="left">Type I diabetes mellitus</td>
<td valign="top" align="center">0.00475</td>
<td valign="top" align="center">0.044185</td>
<td valign="top" align="center">0.036663</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<title>Discussion</title>
<p>Cervical cancer represents a major public health problem even in developed countries such as United States, with 14,000 new cases and 4,285 deaths estimated in 2022 (<xref ref-type="bibr" rid="B20">20</xref>). Increasing incidence of cervical adenocarcinoma and attenuation of earlier declines for cervical squamous cell carcinoma emphasize the importance of the need for improved therapeutic options to reduce the burden of cervical cancer (<xref ref-type="bibr" rid="B21">21</xref>). The most important risk factors affecting prognosis of cervical cancer are stage, status of the lymph nodes, tumor volume, depth of tumor invasion into the cervical stroma, and Lymphovascular space invasion (LVSI) (<xref ref-type="bibr" rid="B22">22</xref>). The 5-year OS decreased significantly with rising stages (<xref ref-type="bibr" rid="B23">23</xref>). A substantial percentage of patients with advanced cervical cancer will undergo recurrence and poor prognosis (<xref ref-type="bibr" rid="B2">2</xref>, <xref ref-type="bibr" rid="B24">24</xref>), and the recurrence rate fluctuates from 10% to 74% (<xref ref-type="bibr" rid="B25">25</xref>). Tumor stage was one of the most pivotal factors related to recurrence: 10% for stage IB, 17% for stage IIA, 23% for stage IIB, 42% for stage III, and 74% for stage IVA (<xref ref-type="bibr" rid="B26">26</xref>). Prognosis of metastatic and recurrent cervical cancer was extremely poor with a median survival time of 12 months (<xref ref-type="bibr" rid="B27">27</xref>). Therefore, improved therapeutic options in management of cervical cancer, especially locally advanced cervical cancer is in urgent need.</p>
<p>High-risk subtypes of the human papilloma virus (HPV) are the cause of cervical cancer (<xref ref-type="bibr" rid="B28">28</xref>, <xref ref-type="bibr" rid="B29">29</xref>). Viral oncoproteins E6 and E7 leads to dysregulation of p53 and Hypoxia-inducible factor-1 (HIF-1) alpha, thus affecting cell cycle proteins and VEGF expression (<xref ref-type="bibr" rid="B28">28</xref>&#x2013;<xref ref-type="bibr" rid="B32">32</xref>). HPV is highly immunogenic and elicits immune responses in humans; thus, immune response might play important roles in carcinogenesis of cervical cancer. Incidence of cervical cancer was substantial declined in countries with high HPV vaccine coverage (<xref ref-type="bibr" rid="B33">33</xref>&#x2013;<xref ref-type="bibr" rid="B35">35</xref>). Immunotherapy fights against tumor cells through activating endogenous immune response, which seems to switch on the new frontier of the anticancer treatment (<xref ref-type="bibr" rid="B36">36</xref>). Immunotherapy includes different approaches, such as active immunotherapy (vaccine), passive immunotherapy (adoptive cellular transfer, antibodies, and cytokines), and immunomodulation (cyclooxygenase 2 inhibitor).</p>
<p>Discovery of immune checkpoint such as CTLA-4 and PD-1 plays an indispensable role in the development of cancer immunotherapy. It was surprising that immune checkpoint inhibitors anti&#x2013;CTLA-4 and anti&#x2013;PD-1 displayed enormous success in solid tumors (<xref ref-type="bibr" rid="B37">37</xref>). Similar to other solid tumors, novel immunotherapeutic approaches, such as immune checkpoint inhibitors, have shown encouraging results in cervical cancer (<xref ref-type="bibr" rid="B12">12</xref>). Implementing immunotherapeutic approaches earlier in advanced cervical cancer would seem to be most appropriate. However, the objective response rate with anti&#x2013;PD-1/PD-L1 monotherapy is hovering at 20%. Moreover, immune-related toxicities and severe adverse effects can occur during PD-1/PD-L1 blockade therapy (<xref ref-type="bibr" rid="B38">38</xref>). PD-1/PD-L1 blockade therapy does not demonstrate efficacy in almost 80% of patients with cervical cancer, suggesting that the potential mechanism of PD-1/PD-L1 in immunotherapy remains to be further clarified. Thus, new immune checkpoint inhibitors or comprehensive understanding of specific mechanism underlying PD-1/PD-L1 regulation in carcinogenesis is in urgent need.</p>
<p>Immune infiltration of cervical cancer determines the immune activation of tumor microenvironment and is related with clinical outcome of patients. In this study, immunophenoscore of 28 infiltrated immune cells in TCGA cervical cancer was calculated. Univariate Cox regression analysis showed that activated B cell, activated CD8 T cell, eosinophil, monocyte, activated CD4 T cell, effector memory CD8 T cell, immature B cell, and plasmacytoid dendritic cell were strongly associated with OS of cervical cancer. WGCNA revealed that magenta module was the most relevant module to immune infiltration. In total, we identified 10 kinds of infiltrated immune cells with strong correlation with magenta module. These were as follows: activated B cell (R = 0.88), activated CD8 T cell (R = 0.82), activated dendritic cell (R = 0.79), central memory CD4 T cell (R = 0.78), effector memory CD8 T cell (R = 0.85), immature B cell (R = 0.93), macrophage (R = 0.72), Myeloid-derived suppressor cells (MDSC) (R = 0.85), regulatory T cell (R = 0.71), T follicular helper cell (R = 0.72), and type 1 T helper cell (R = 0.83).</p>
<p>By the use of WGCNA, magenta module showed the highest correlation with the immune infiltration of cervical cancer. The magenta module contained 609 immune-related genes, and 153 genes were picked up for further analysis after taking intersection of prognosis and immune infiltration. Through LASSO analysis, 15 immune-related genes (LAG3, CD74, CCL22, CH25H, OLR1, MIAT, BATF, IKZF3, TRARG, ACSL6, C11orf21, GTSF1, APOL1, CD1C, and LINC00158) were included in prognostic classifier. In addition to CTLA-4&#x2013; and PD-1/PD-L1&#x2013;targeted cancer immunotherapy, LAG3 (lymphocyte activation gene-3, CD223) is the third clinically targeted inhibitory receptor (<xref ref-type="bibr" rid="B39">39</xref>). CD74 (invariant chain) plays a dispensable role in process of immune systems that it participates in antigen presentation, B-cell differentiation, and inflammatory signaling. CD74 has the potential to be a therapeutic target in cancer and autoimmune disease (<xref ref-type="bibr" rid="B40">40</xref>). The chemokine CCL22 promoted regulatory T cell communication with dendritic cells to control immunity and was associated with poor prognosis (<xref ref-type="bibr" rid="B41">41</xref>, <xref ref-type="bibr" rid="B42">42</xref>). CH25H produces 25-hydroxycholesterol, which inhibited tumor-derived extracellular vesicles uptake and correlated with prognosis in patients with melanoma (<xref ref-type="bibr" rid="B43">43</xref>). OLR1 (oxidized low-density lipoprotein (LDL) receptor 1) was a possible link between obesity, dyslipidemia, and cancer. OLR1 played carcinogenic role by activating Nuclear factor kappa B (NF-&#x3ba;B) pathway to promote proliferation and migration and to inhibit apoptosis and <italic>de novo</italic> lipogenesis (<xref ref-type="bibr" rid="B44">44</xref>). MIAT Long non-coding RNAs (lncRNA) was overexpressed in a number of malignancies and caused poor prognosis (<xref ref-type="bibr" rid="B45">45</xref>&#x2013;<xref ref-type="bibr" rid="B47">47</xref>). Basic leucine zipper transcription factor, ATF-like (BATF) was an important transcription factor regulating differentiation of early effector CD8<sup>+</sup> T cells (<xref ref-type="bibr" rid="B48">48</xref>) and was a prognostic indicator for patients with colon cancer (<xref ref-type="bibr" rid="B49">49</xref>). IKZF3 promoted growth of multiple tumors to cause poor prognosis (<xref ref-type="bibr" rid="B50">50</xref>).</p>
<p>Immune infiltration played important roles in the survival of cervical cancer. In this study, we identified 15 immune infiltration associated genes in cervical cancer and built a prognostic signature. Immune scores depended on the expression of these 15 genes and were associated with the survival of cervical cancer. High immune scores meant good prognosis. GSEA analysis showed that the 15-gene prognostic signature was obviously associated with PD-L1 expression and PD-1 checkpoint pathway in cancer. The prognostic signature could provide basis for potential immunotherapy in the future. Similarly, other researchers also constructed prognostic signatures for cervical cancer based on immune-related genes (<xref ref-type="bibr" rid="B51">51</xref>, <xref ref-type="bibr" rid="B52">52</xref>). Compared with our study, although these studies use different types and different quantity of genes, all the signatures can well predict the prognosis of cervical cancer. However, the study has several limitations. First, no <italic>in vitro</italic> or <italic>in vivo</italic> molecular experiment was performed to verify our analysis. Second, our study was a retrospective study. Thus, prospective study is in need to validate the findings of our study in the future.</p>
</sec>
<sec id="s5">
<title>Conclusion</title>
<p>In conclusion, we successfully constructed a 15-gene prognostic signature with powerful predictive function. Differences in the OS of high- and low-risk groups are implicated in immune infiltration, tumor microenvironment, PD-L1 expression, and PD-1 checkpoint pathway. These findings revealed the underlying mechanism of immunotherapy and provided basis for cervical cancer pathogenesis and clinical treatment.</p>
</sec>
<sec id="s6" sec-type="data-availability">
<title>Data Availability Statement</title>
<p>The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding authors.</p>
</sec>
<sec id="s7" sec-type="author-contributions">
<title>Author Contributions</title>
<p>RS is responsible for experimental design. CWJ and HB are responsible for instrument operation. JX is responsible for data analysis. LZ and CJJ are for providing overall ideas. All authors contributed to the article and approved the submitted version.</p>
</sec>
<sec id="s8" sec-type="funding-information">
<title>Funding</title>
<p>This work was supported by the National Natural Science Foundation of China (grant no. 82002736 to CJJ) and the clinical characteristic medical technology cultivation plan of Shanghai General Hospital (no. 02.DY12.06.22.07 to LZ and no. 02.DY12.06.22.12 to JX).</p>
</sec>
<sec id="s9" 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="s10" 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>
<ack>
<title>Acknowledgments</title>
<p>We thank Dr. Jianming Zeng (University of Macau) and all the members of his bioinformatics team, biotrainee, for generously sharing their experience and codes, especially Xiaojie Sun.</p>
</ack>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Arbyn</surname> <given-names>M</given-names>
</name>
<name>
<surname>Weiderpass</surname> <given-names>E</given-names>
</name>
<name>
<surname>Bruni</surname> <given-names>L</given-names>
</name>
<name>
<surname>de Sanjos&#xe9;</surname> <given-names>S</given-names>
</name>
<name>
<surname>Saraiya</surname> <given-names>M</given-names>
</name>
<name>
<surname>Ferlay</surname> <given-names>J</given-names>
</name>
<etal/>
</person-group>. <article-title>Estimates of Incidence and Mortality of Cervical Cancer in 2018: A Worldwide Analysis</article-title>. <source>Lancet Global Health</source> (<year>2020</year>) <volume>8</volume>(<issue>2</issue>):<page-range>e191&#x2013;203</page-range>. doi: <pub-id pub-id-type="doi">10.1016/S2214-109X(19)30482-6</pub-id>
</citation>
</ref>
<ref id="B2">
<label>2</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Monk</surname> <given-names>BJ</given-names>
</name>
<name>
<surname>Tewari</surname> <given-names>KS</given-names>
</name>
<name>
<surname>Koh</surname> <given-names>WJ</given-names>
</name>
</person-group>. <article-title>Multimodality Therapy for Locally Advanced Cervical Carcinoma: State of the Art and Future Directions</article-title>. <source>J Clin Oncol Off J Am Soc Clin Oncol</source> (<year>2007</year>) <volume>25</volume>(<issue>20</issue>):<page-range>2952&#x2013;65</page-range>. doi: <pub-id pub-id-type="doi">10.1200/JCO.2007.10.8324</pub-id>
</citation>
</ref>
<ref id="B3">
<label>3</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chung</surname> <given-names>HC</given-names>
</name>
<name>
<surname>Ros</surname> <given-names>W</given-names>
</name>
<name>
<surname>Delord</surname> <given-names>JP</given-names>
</name>
<name>
<surname>Perets</surname> <given-names>R</given-names>
</name>
<name>
<surname>Italiano</surname> <given-names>A</given-names>
</name>
<name>
<surname>Shapira-Frommer</surname> <given-names>R</given-names>
</name>
<etal/>
</person-group>. <article-title>Efficacy and Safety of Pembrolizumab in Previously Treated Advanced Cervical Cancer: Results From the Phase II KEYNOTE-158 Study</article-title>. <source>J Clin Oncol Off J Am Soc Clin Oncol</source> (<year>2019</year>) <volume>37</volume>(<issue>17</issue>):<page-range>1470&#x2013;8</page-range>. doi: <pub-id pub-id-type="doi">10.1200/JCO.18.01265</pub-id>
</citation>
</ref>
<ref id="B4">
<label>4</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cohen</surname> <given-names>PA</given-names>
</name>
<name>
<surname>Jhingran</surname> <given-names>A</given-names>
</name>
<name>
<surname>Oaknin</surname> <given-names>A</given-names>
</name>
<name>
<surname>Denny</surname> <given-names>L</given-names>
</name>
</person-group>. <article-title>Cervical Cancer</article-title>. <source>Lancet</source> (<year>2019</year>) <volume>393</volume>(<issue>10167</issue>):<page-range>169&#x2013;82</page-range>. doi: <pub-id pub-id-type="doi">10.1016/S0140-6736(18)32470-X</pub-id>
</citation>
</ref>
<ref id="B5">
<label>5</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sanmamed</surname> <given-names>MF</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>L</given-names>
</name>
</person-group>. <article-title>A Paradigm Shift in Cancer Immunotherapy: From Enhancement to Normalization</article-title>. <source>Cell</source> (<year>2018</year>) <volume>175</volume>(<issue>2</issue>):<page-range>313&#x2013;26</page-range>. doi: <pub-id pub-id-type="doi">10.1016/j.cell.2018.09.035</pub-id>
</citation>
</ref>
<ref id="B6">
<label>6</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tran</surname> <given-names>E</given-names>
</name>
<name>
<surname>Robbins</surname> <given-names>PF</given-names>
</name>
<name>
<surname>Rosenberg</surname> <given-names>SA</given-names>
</name>
</person-group>. <article-title>'Final Common Pathway' of Human Cancer Immunotherapy: Targeting Random Somatic Mutations</article-title>. <source>Nat Immunol</source> (<year>2017</year>) <volume>18</volume>(<issue>3</issue>):<page-range>255&#x2013;62</page-range>. doi: <pub-id pub-id-type="doi">10.1038/ni.3682</pub-id>
</citation>
</ref>
<ref id="B7">
<label>7</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Asaoka</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Ijichi</surname> <given-names>H</given-names>
</name>
<name>
<surname>Koike</surname> <given-names>K</given-names>
</name>
</person-group>. <article-title>PD-1 Blockade in Tumors With Mismatch-Repair Deficiency</article-title>. <source>N Engl J Med</source> (<year>2015</year>) <volume>373</volume>(<issue>20</issue>):<fpage>1979</fpage>. doi: <pub-id pub-id-type="doi">10.1056/NEJMc1510353</pub-id>
</citation>
</ref>
<ref id="B8">
<label>8</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hugo</surname> <given-names>W</given-names>
</name>
<name>
<surname>Zaretsky</surname> <given-names>JM</given-names>
</name>
<name>
<surname>Sun</surname> <given-names>L</given-names>
</name>
<name>
<surname>Song</surname> <given-names>C</given-names>
</name>
<name>
<surname>Moreno</surname> <given-names>BH</given-names>
</name>
<name>
<surname>Hu-Lieskovan</surname> <given-names>S</given-names>
</name>
<etal/>
</person-group>. <article-title>Genomic and Transcriptomic Features of Response to Anti-PD-1 Therapy in Metastatic Melanoma</article-title>. <source>Cell</source> (<year>2016</year>) <volume>165</volume>(<issue>1</issue>):<fpage>35</fpage>&#x2013;<lpage>44</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.cell.2016.02.065</pub-id>
</citation>
</ref>
<ref id="B9">
<label>9</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Powles</surname> <given-names>T</given-names>
</name>
<name>
<surname>Eder</surname> <given-names>JP</given-names>
</name>
<name>
<surname>Fine</surname> <given-names>GD</given-names>
</name>
<name>
<surname>Braiteh</surname> <given-names>FS</given-names>
</name>
<name>
<surname>Loriot</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Cruz</surname> <given-names>C</given-names>
</name>
<etal/>
</person-group>. <article-title>MPDL3280A (Anti-PD-L1) Treatment Leads to Clinical Activity in Metastatic Bladder Cancer</article-title>. <source>Nature</source> (<year>2014</year>) <volume>515</volume>(<issue>7528</issue>):<page-range>558&#x2013;62</page-range>. doi: <pub-id pub-id-type="doi">10.1038/nature13904</pub-id>
</citation>
</ref>
<ref id="B10">
<label>10</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rizvi</surname> <given-names>NA</given-names>
</name>
<name>
<surname>Hellmann</surname> <given-names>MD</given-names>
</name>
<name>
<surname>Snyder</surname> <given-names>A</given-names>
</name>
<name>
<surname>Kvistborg</surname> <given-names>P</given-names>
</name>
<name>
<surname>Makarov</surname> <given-names>V</given-names>
</name>
<name>
<surname>Havel</surname> <given-names>JJ</given-names>
</name>
<etal/>
</person-group>. <article-title>Cancer Immunology. Mutational Landscape Determines Sensitivity to PD-1 Blockade in Non-Small Cell Lung Cancer</article-title>. <source>Science</source> (<year>2015</year>) <volume>348</volume>(<issue>6230</issue>):<page-range>124&#x2013;8</page-range>. doi: <pub-id pub-id-type="doi">10.1126/science.aaa1348</pub-id>
</citation>
</ref>
<ref id="B11">
<label>11</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>van Rooij</surname> <given-names>N</given-names>
</name>
<name>
<surname>van Buuren</surname> <given-names>MM</given-names>
</name>
<name>
<surname>Philips</surname> <given-names>D</given-names>
</name>
<name>
<surname>Velds</surname> <given-names>A</given-names>
</name>
<name>
<surname>Toebes</surname> <given-names>M</given-names>
</name>
<name>
<surname>Heemskerk</surname> <given-names>B</given-names>
</name>
<etal/>
</person-group>. <article-title>Tumor Exome Analysis Reveals Neoantigen-Specific T-Cell Reactivity in an Ipilimumab-Responsive Melanoma</article-title>. <source>J Clin Oncol Off J Am Soc Clin Oncol</source> (<year>2013</year>) <volume>31</volume>(<issue>32</issue>):<page-range>e439&#x2013;442</page-range>. doi: <pub-id pub-id-type="doi">10.1200/JCO.2012.47.7521</pub-id>
</citation>
</ref>
<ref id="B12">
<label>12</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Marabelle</surname> <given-names>A</given-names>
</name>
<name>
<surname>Fakih</surname> <given-names>M</given-names>
</name>
<name>
<surname>Lopez</surname> <given-names>J</given-names>
</name>
<name>
<surname>Shah</surname> <given-names>M</given-names>
</name>
<name>
<surname>Shapira-Frommer</surname> <given-names>R</given-names>
</name>
<name>
<surname>Nakagawa</surname> <given-names>K</given-names>
</name>
<etal/>
</person-group>. <article-title>Association of Tumour Mutational Burden With Outcomes in Patients With Advanced Solid Tumours Treated With Pembrolizumab: Prospective Biomarker Analysis of the Multicohort, Open-Label, Phase 2 KEYNOTE-158 Study</article-title>. <source>Lancet Oncol</source> (<year>2020</year>) <volume>21</volume>(<issue>10</issue>):<page-range>1353&#x2013;65</page-range>. doi: <pub-id pub-id-type="doi">10.1016/S1470-2045(20)30445-9</pub-id>
</citation>
</ref>
<ref id="B13">
<label>13</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Finotello</surname> <given-names>F</given-names>
</name>
<name>
<surname>Trajanoski</surname> <given-names>Z</given-names>
</name>
</person-group>. <article-title>Quantifying Tumor-Infiltrating Immune Cells From Transcriptomics Data</article-title>. <source>Cancer Immunology Immunotherapy CII</source> (<year>2018</year>) <volume>67</volume>(<issue>7</issue>):<page-range>1031&#x2013;40</page-range>. doi: <pub-id pub-id-type="doi">10.1007/s00262-018-2150-z</pub-id>
</citation>
</ref>
<ref id="B14">
<label>14</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Langfelder</surname> <given-names>P</given-names>
</name>
<name>
<surname>Horvath</surname> <given-names>S</given-names>
</name>
</person-group>. <article-title>WGCNA: An R Package for Weighted Correlation Network Analysis</article-title>. <source>BMC Bioinf</source> (<year>2008</year>) <volume>9</volume>:<fpage>559</fpage>. doi: <pub-id pub-id-type="doi">10.1186/1471-2105-9-559</pub-id>
</citation>
</ref>
<ref id="B15">
<label>15</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dong</surname> <given-names>J</given-names>
</name>
<name>
<surname>Horvath</surname> <given-names>S</given-names>
</name>
</person-group>. <article-title>Understanding Network Concepts in Modules</article-title>. <source>BMC Syst Biol</source> (<year>2007</year>) <volume>1</volume>:<fpage>24</fpage>. doi: <pub-id pub-id-type="doi">10.1186/1752-0509-1-24</pub-id>
</citation>
</ref>
<ref id="B16">
<label>16</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Horvath</surname> <given-names>S</given-names>
</name>
<name>
<surname>Dong</surname> <given-names>J</given-names>
</name>
</person-group>. <article-title>Geometric Interpretation of Gene Coexpression Network Analysis</article-title>. <source>PloS Comput Biol</source> (<year>2008</year>) <volume>4</volume>(<issue>8</issue>):<elocation-id>e1000117</elocation-id>. doi: <pub-id pub-id-type="doi">10.1371/journal.pcbi.1000117</pub-id>
</citation>
</ref>
<ref id="B17">
<label>17</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Langfelder</surname> <given-names>P</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>B</given-names>
</name>
<name>
<surname>Horvath</surname> <given-names>S</given-names>
</name>
</person-group>. <article-title>Defining Clusters From a Hierarchical Cluster Tree: The Dynamic Tree Cut Package for R</article-title>. <source>Bioinf (Oxford England)</source> (<year>2008</year>) <volume>24</volume>(<issue>5</issue>):<page-range>719&#x2013;20</page-range>. doi: <pub-id pub-id-type="doi">10.1093/bioinformatics/btm563</pub-id>
</citation>
</ref>
<ref id="B18">
<label>18</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ferlay</surname> <given-names>J</given-names>
</name>
<name>
<surname>Soerjomataram</surname> <given-names>I</given-names>
</name>
<name>
<surname>Dikshit</surname> <given-names>R</given-names>
</name>
<name>
<surname>Eser</surname> <given-names>S</given-names>
</name>
<name>
<surname>Mathers</surname> <given-names>C</given-names>
</name>
<name>
<surname>Rebelo</surname> <given-names>M</given-names>
</name>
<etal/>
</person-group>. <article-title>Cancer Incidence and Mortality Worldwide: Sources, Methods and Major Patterns in GLOBOCAN 2012</article-title>. <source>Int J Cancer</source> (<year>2015</year>) <volume>136</volume>(<issue>5</issue>):<page-range>E359&#x2013;386</page-range>. doi: <pub-id pub-id-type="doi">10.1002/ijc.29210</pub-id>
</citation>
</ref>
<ref id="B19">
<label>19</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Subramanian</surname> <given-names>A</given-names>
</name>
<name>
<surname>Kuehn</surname> <given-names>H</given-names>
</name>
<name>
<surname>Gould</surname> <given-names>J</given-names>
</name>
<name>
<surname>Tamayo</surname> <given-names>P</given-names>
</name>
<name>
<surname>Mesirov</surname> <given-names>JP</given-names>
</name>
</person-group>. <article-title>GSEA-P: A Desktop Application for Gene Set Enrichment Analysis</article-title>. <source>Bioinformatics</source> (<year>2007</year>) <volume>23</volume>(<issue>23</issue>):<page-range>3251&#x2013;3</page-range>. doi: <pub-id pub-id-type="doi">10.1093/bioinformatics/btm369</pub-id>
</citation>
</ref>
<ref id="B20">
<label>20</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Siegel</surname> <given-names>RL</given-names>
</name>
<name>
<surname>Miller</surname> <given-names>KD</given-names>
</name>
<name>
<surname>Fuchs</surname> <given-names>HE</given-names>
</name>
<name>
<surname>Jemal</surname> <given-names>A</given-names>
</name>
</person-group>. <article-title>Cancer Statistics, 2022</article-title>. <source>CA Cancer J Clin</source> (<year>2022</year>) <volume>72</volume>(<issue>1</issue>):<fpage>7</fpage>&#x2013;<lpage>33</lpage>. doi: <pub-id pub-id-type="doi">10.3322/caac.21708</pub-id>
</citation>
</ref>
<ref id="B21">
<label>21</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Islami</surname> <given-names>F</given-names>
</name>
<name>
<surname>Fedewa</surname> <given-names>SA</given-names>
</name>
<name>
<surname>Jemal</surname> <given-names>A</given-names>
</name>
</person-group>. <article-title>Trends in Cervical Cancer Incidence Rates by Age, Race/Ethnicity, Histological Subtype, and Stage at Diagnosis in the United States</article-title>. <source>Prev Med</source> (<year>2019</year>) <volume>123</volume>:<page-range>316&#x2013;23</page-range>. doi: <pub-id pub-id-type="doi">10.1016/j.ypmed.2019.04.010</pub-id>
</citation>
</ref>
<ref id="B22">
<label>22</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Peters</surname> <given-names>WA</given-names>
<suffix>3rd</suffix>
</name>
<name>
<surname>Liu</surname> <given-names>PY</given-names>
</name>
<name>
<surname>Barrett</surname> <given-names>RJ</given-names>
<suffix>2nd</suffix>
</name>
<name>
<surname>Stock</surname> <given-names>RJ</given-names>
</name>
<name>
<surname>Monk</surname> <given-names>BJ</given-names>
</name>
<name>
<surname>Berek</surname> <given-names>JS</given-names>
</name>
<etal/>
</person-group>. <article-title>Concurrent Chemotherapy and Pelvic Radiation Therapy Compared With Pelvic Radiation Therapy Alone as Adjuvant Therapy After Radical Surgery in High-Risk Early-Stage Cancer of the Cervix</article-title>. <source>J Clin Oncol Off J Am Soc Clin Oncol</source> (<year>2000</year>) <volume>18</volume>(<issue>8</issue>):<page-range>1606&#x2013;13</page-range>. doi: <pub-id pub-id-type="doi">10.1200/JCO.2000.18.8.1606</pub-id>
</citation>
</ref>
<ref id="B23">
<label>23</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sturdza</surname> <given-names>A</given-names>
</name>
<name>
<surname>P&#xf6;tter</surname> <given-names>R</given-names>
</name>
<name>
<surname>Fokdal</surname> <given-names>LU</given-names>
</name>
<name>
<surname>Haie-Meder</surname> <given-names>C</given-names>
</name>
<name>
<surname>Tan</surname> <given-names>LT</given-names>
</name>
<name>
<surname>Mazeron</surname> <given-names>R</given-names>
</name>
<etal/>
</person-group>. <article-title>Image Guided Brachytherapy in Locally Advanced Cervical Cancer: Improved Pelvic Control and Survival in RetroEMBRACE, a Multicenter Cohort Study</article-title>. <source>Radiotherapy Oncol J Eur Soc Ther Radiol Oncol</source> (<year>2016</year>) <volume>120</volume>(<issue>3</issue>):<page-range>428&#x2013;33</page-range>. doi: <pub-id pub-id-type="doi">10.1016/j.radonc.2016.03.011</pub-id>
</citation>
</ref>
<ref id="B24">
<label>24</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liontos</surname> <given-names>M</given-names>
</name>
<name>
<surname>Kyriazoglou</surname> <given-names>A</given-names>
</name>
<name>
<surname>Dimitriadis</surname> <given-names>I</given-names>
</name>
<name>
<surname>Dimopoulos</surname> <given-names>MA</given-names>
</name>
<name>
<surname>Bamias</surname> <given-names>A</given-names>
</name>
</person-group>. <article-title>Systemic Therapy in Cervical Cancer: 30 Years in Review</article-title>. <source>Crit Rev Oncology/Hematology</source> (<year>2019</year>) <volume>137</volume>:<fpage>9</fpage>&#x2013;<lpage>17</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.critrevonc.2019.02.009</pub-id>
</citation>
</ref>
<ref id="B25">
<label>25</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Peiretti</surname> <given-names>M</given-names>
</name>
<name>
<surname>Zapardiel</surname> <given-names>I</given-names>
</name>
<name>
<surname>Zanagnolo</surname> <given-names>V</given-names>
</name>
<name>
<surname>Landoni</surname> <given-names>F</given-names>
</name>
<name>
<surname>Morrow</surname> <given-names>CP</given-names>
</name>
<name>
<surname>Maggioni</surname> <given-names>A</given-names>
</name>
</person-group>. <article-title>Management of Recurrent Cervical Cancer: A Review of the Literature</article-title>. <source>Surg Oncol</source> (<year>2012</year>) <volume>21</volume>(<issue>2</issue>):<page-range>e59&#x2013;66</page-range>. doi: <pub-id pub-id-type="doi">10.1016/j.suronc.2011.12.008</pub-id>
</citation>
</ref>
<ref id="B26">
<label>26</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Perez</surname> <given-names>CA</given-names>
</name>
<name>
<surname>Grigsby</surname> <given-names>PW</given-names>
</name>
<name>
<surname>Nene</surname> <given-names>SM</given-names>
</name>
<name>
<surname>Camel</surname> <given-names>HM</given-names>
</name>
<name>
<surname>Galakatos</surname> <given-names>A</given-names>
</name>
<name>
<surname>Kao</surname> <given-names>MS</given-names>
</name>
<etal/>
</person-group>. <article-title>Effect of Tumor Size on the Prognosis of Carcinoma of the Uterine Cervix Treated With Irradiation Alone</article-title>. <source>Cancer</source> (<year>1992</year>) <volume>69</volume>(<issue>11</issue>):<page-range>2796&#x2013;806</page-range>. doi: <pub-id pub-id-type="doi">10.1002/1097-0142(19920601)69:11&lt;2796::AID-CNCR2820691127&gt;3.0.CO;2-O</pub-id>
</citation>
</ref>
<ref id="B27">
<label>27</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Elit</surname> <given-names>LM</given-names>
</name>
<name>
<surname>Hirte</surname> <given-names>H</given-names>
</name>
</person-group>. <article-title>Management of Advanced or Recurrent Cervical Cancer: Chemotherapy and Beyond</article-title>. <source>Expert Rev Anticancer Ther</source> (<year>2014</year>) <volume>14</volume>(<issue>3</issue>):<page-range>319&#x2013;32</page-range>. doi: <pub-id pub-id-type="doi">10.1586/14737140.2014.866041</pub-id>
</citation>
</ref>
<ref id="B28">
<label>28</label>
<citation citation-type="journal">
<article-title>Wardak S: Human Papillomavirus (HPV) and Cervical Cancer</article-title>. <source>Medycyna doswiadczalna i mikrobiologia</source> (<year>2016</year>) <volume>68</volume>(<issue>1</issue>):<fpage>73</fpage>&#x2013;<lpage>84</lpage>.</citation>
</ref>
<ref id="B29">
<label>29</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Brianti</surname> <given-names>P</given-names>
</name>
<name>
<surname>De Flammineis</surname> <given-names>E</given-names>
</name>
<name>
<surname>Mercuri</surname> <given-names>SR</given-names>
</name>
</person-group>. <article-title>Review of HPV-Related Diseases and Cancers</article-title>. <source>New Microbiol</source> (<year>2017</year>) <volume>40</volume>(<issue>2</issue>):<page-range>80&#x2013;5</page-range>. doi: <pub-id pub-id-type="doi">10.1016/S0140-6736(13)60022-7</pub-id>
</citation>
</ref>
<ref id="B30">
<label>30</label>
<citation citation-type="journal">
<article-title>Burd EM: Human Papillomavirus and Cervical Cancer</article-title>. <source>Clin Microbiol Rev</source> (<year>2003</year>) <volume>16</volume>(<issue>1</issue>):<fpage>1</fpage>&#x2013;<lpage>17</lpage>. doi: <pub-id pub-id-type="doi">10.1128/CMR.16.1.1-17.2003</pub-id>
</citation>
</ref>
<ref id="B31">
<label>31</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Monk</surname> <given-names>BJ</given-names>
</name>
<name>
<surname>Willmott</surname> <given-names>LJ</given-names>
</name>
<name>
<surname>Sumner</surname> <given-names>DA</given-names>
</name>
</person-group>. <article-title>Anti-Angiogenesis Agents in Metastatic or Recurrent Cervical Cancer</article-title>. <source>Gynecologic Oncol</source> (<year>2010</year>) <volume>116</volume>(<issue>2</issue>):<page-range>181&#x2013;6</page-range>. doi: <pub-id pub-id-type="doi">10.1016/j.ygyno.2009.09.033</pub-id>
</citation>
</ref>
<ref id="B32">
<label>32</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nakamura</surname> <given-names>M</given-names>
</name>
<name>
<surname>Bodily</surname> <given-names>JM</given-names>
</name>
<name>
<surname>Beglin</surname> <given-names>M</given-names>
</name>
<name>
<surname>Kyo</surname> <given-names>S</given-names>
</name>
<name>
<surname>Inoue</surname> <given-names>M</given-names>
</name>
<name>
<surname>Laimins</surname> <given-names>LA</given-names>
</name>
</person-group>. <article-title>Hypoxia-Specific Stabilization of HIF-1alpha by Human Papillomaviruses</article-title>. <source>Virology</source> (<year>2009</year>) <volume>387</volume>(<issue>2</issue>):<page-range>442&#x2013;8</page-range>. doi: <pub-id pub-id-type="doi">10.1016/j.virol.2009.02.036</pub-id>
</citation>
</ref>
<ref id="B33">
<label>33</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Brotherton</surname> <given-names>JM</given-names>
</name>
<name>
<surname>Fridman</surname> <given-names>M</given-names>
</name>
<name>
<surname>May</surname> <given-names>CL</given-names>
</name>
<name>
<surname>Chappell</surname> <given-names>G</given-names>
</name>
<name>
<surname>Saville</surname> <given-names>AM</given-names>
</name>
<name>
<surname>Gertig</surname> <given-names>DM</given-names>
</name>
</person-group>. <article-title>Early Effect of the HPV Vaccination Programme on Cervical Abnormalities in Victoria, Australia: An Ecological Study</article-title>. <source>Lancet</source> (<year>2011</year>) <volume>377</volume>(<issue>9783</issue>):<page-range>2085&#x2013;92</page-range>. doi: <pub-id pub-id-type="doi">10.1016/S0140-6736(11)60551-5</pub-id>
</citation>
</ref>
<ref id="B34">
<label>34</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Drolet</surname> <given-names>M</given-names>
</name>
<name>
<surname>B&#xe9;nard</surname> <given-names>&#xc9;</given-names>
</name>
<name>
<surname>P&#xe9;rez</surname> <given-names>N</given-names>
</name>
<name>
<surname>Brisson</surname> <given-names>M</given-names>
</name>
</person-group>. <article-title>Population-Level Impact and Herd Effects Following the Introduction of Human Papillomavirus Vaccination Programmes: Updated Systematic Review and Meta-Analysis</article-title>. <source>Lancet</source> (<year>2019</year>) <volume>394</volume>(<issue>10197</issue>):<fpage>497</fpage>&#x2013;<lpage>509</lpage>. doi: <pub-id pub-id-type="doi">10.1016/S0140-6736(19)30298-3</pub-id>
</citation>
</ref>
<ref id="B35">
<label>35</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gargano</surname> <given-names>JW</given-names>
</name>
<name>
<surname>Zhou</surname> <given-names>F</given-names>
</name>
<name>
<surname>Stokley</surname> <given-names>S</given-names>
</name>
<name>
<surname>Markowitz</surname> <given-names>LE</given-names>
</name>
</person-group>. <article-title>Human Papillomavirus Vaccination in Commercially-Insured Vaccine-Eligible Males and Females, United States, 2007-2014</article-title>. <source>Vaccine</source> (<year>2018</year>) <volume>36</volume>(<issue>23</issue>):<page-range>3381&#x2013;6</page-range>. doi: <pub-id pub-id-type="doi">10.1016/j.vaccine.2018.03.045</pub-id>
</citation>
</ref>
<ref id="B36">
<label>36</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ventriglia</surname> <given-names>J</given-names>
</name>
<name>
<surname>Paciolla</surname> <given-names>I</given-names>
</name>
<name>
<surname>Pisano</surname> <given-names>C</given-names>
</name>
<name>
<surname>Cecere</surname> <given-names>SC</given-names>
</name>
<name>
<surname>Di Napoli</surname> <given-names>M</given-names>
</name>
<name>
<surname>Tambaro</surname> <given-names>R</given-names>
</name>
<etal/>
</person-group>. <article-title>Immunotherapy in Ovarian, Endometrial and Cervical Cancer: State of the Art and Future Perspectives</article-title>. <source>Cancer Treat Rev</source> (<year>2017</year>) <volume>59</volume>:<page-range>109&#x2013;16</page-range>. doi: <pub-id pub-id-type="doi">10.1016/j.ctrv.2017.07.008</pub-id>
</citation>
</ref>
<ref id="B37">
<label>37</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Balar</surname> <given-names>AV</given-names>
</name>
<name>
<surname>Weber</surname> <given-names>JS</given-names>
</name>
</person-group>. <article-title>PD-1 and PD-L1 Antibodies in Cancer: Current Status and Future Directions</article-title>. <source>Cancer Immunology Immunotherapy CII</source> (<year>2017</year>) <volume>66</volume>(<issue>5</issue>):<page-range>551&#x2013;64</page-range>. doi: <pub-id pub-id-type="doi">10.1007/s00262-017-1954-6</pub-id>
</citation>
</ref>
<ref id="B38">
<label>38</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Topalian</surname> <given-names>SL</given-names>
</name>
<name>
<surname>Hodi</surname> <given-names>FS</given-names>
</name>
<name>
<surname>Brahmer</surname> <given-names>JR</given-names>
</name>
<name>
<surname>Gettinger</surname> <given-names>SN</given-names>
</name>
<name>
<surname>Smith</surname> <given-names>DC</given-names>
</name>
<name>
<surname>McDermott</surname> <given-names>DF</given-names>
</name>
<etal/>
</person-group>. <article-title>Safety, Activity, and Immune Correlates of Anti-PD-1 Antibody in Cancer</article-title>. <source>N Engl J Med</source> (<year>2012</year>) <volume>366</volume>(<issue>26</issue>):<page-range>2443&#x2013;54</page-range>. doi: <pub-id pub-id-type="doi">10.1056/NEJMoa1200690</pub-id>
</citation>
</ref>
<ref id="B39">
<label>39</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Andrews</surname> <given-names>LP</given-names>
</name>
<name>
<surname>Marciscano</surname> <given-names>AE</given-names>
</name>
<name>
<surname>Drake</surname> <given-names>CG</given-names>
</name>
<name>
<surname>Vignali</surname> <given-names>DA</given-names>
</name>
</person-group>. <article-title>LAG3 (CD223) as a Cancer Immunotherapy Target</article-title>. <source>Immunol Rev</source> (<year>2017</year>) <volume>276</volume>(<issue>1</issue>):<fpage>80</fpage>&#x2013;<lpage>96</lpage>. doi: <pub-id pub-id-type="doi">10.1111/imr.12519</pub-id>
</citation>
</ref>
<ref id="B40">
<label>40</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Borghese</surname> <given-names>F</given-names>
</name>
<name>
<surname>Clanchy</surname> <given-names>FI</given-names>
</name>
</person-group>. <article-title>CD74: An Emerging Opportunity as a Therapeutic Target in Cancer and Autoimmune Disease</article-title>. <source>Expert Opin Ther Targets</source> (<year>2011</year>) <volume>15</volume>(<issue>3</issue>):<page-range>237&#x2013;51</page-range>. doi: <pub-id pub-id-type="doi">10.1517/14728222.2011.550879</pub-id>
</citation>
</ref>
<ref id="B41">
<label>41</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rapp</surname> <given-names>M</given-names>
</name>
<name>
<surname>Wintergerst</surname> <given-names>MWM</given-names>
</name>
<name>
<surname>Kunz</surname> <given-names>WG</given-names>
</name>
<name>
<surname>Vetter</surname> <given-names>VK</given-names>
</name>
<name>
<surname>Knott</surname> <given-names>MML</given-names>
</name>
<name>
<surname>Lisowski</surname> <given-names>D</given-names>
</name>
<etal/>
</person-group>. <article-title>CCL22 Controls Immunity by Promoting Regulatory T Cell Communication With Dendritic Cells in Lymph Nodes</article-title>. <source>J Exp Med</source> (<year>2019</year>) <volume>216</volume>(<issue>5</issue>):<page-range>1170&#x2013;81</page-range>. doi: <pub-id pub-id-type="doi">10.1084/jem.20170277</pub-id>
</citation>
</ref>
<ref id="B42">
<label>42</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>Q</given-names>
</name>
<name>
<surname>Schmoeckel</surname> <given-names>E</given-names>
</name>
<name>
<surname>Kost</surname> <given-names>BP</given-names>
</name>
<name>
<surname>Kuhn</surname> <given-names>C</given-names>
</name>
<name>
<surname>Vattai</surname> <given-names>A</given-names>
</name>
<name>
<surname>Vilsmaier</surname> <given-names>T</given-names>
</name>
<etal/>
</person-group>. <article-title>Higher CCL22+ Cell Infiltration is Associated With Poor Prognosis in Cervical Cancer Patients</article-title>. <source>Cancers</source> (<year>2019</year>) <volume>11</volume>(<issue>12</issue>). doi: <pub-id pub-id-type="doi">10.3390/cancers11122004</pub-id>
</citation>
</ref>
<ref id="B43">
<label>43</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ortiz</surname> <given-names>A</given-names>
</name>
<name>
<surname>Gui</surname> <given-names>J</given-names>
</name>
<name>
<surname>Zahedi</surname> <given-names>F</given-names>
</name>
<name>
<surname>Yu</surname> <given-names>P</given-names>
</name>
<name>
<surname>Cho</surname> <given-names>C</given-names>
</name>
<name>
<surname>Bhattacharya</surname> <given-names>S</given-names>
</name>
<etal/>
</person-group>. <article-title>An Interferon-Driven Oxysterol-Based Defense Against Tumor-Derived Extracellular Vesicles</article-title>. <source>Cancer Cell</source> (<year>2019</year>) <volume>35</volume>(<issue>1</issue>):<fpage>33</fpage>&#x2013;<lpage>45.e36</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.ccell.2018.12.001</pub-id>
</citation>
</ref>
<ref id="B44">
<label>44</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Khaidakov</surname> <given-names>M</given-names>
</name>
<name>
<surname>Mitra</surname> <given-names>S</given-names>
</name>
<name>
<surname>Kang</surname> <given-names>BY</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>X</given-names>
</name>
<name>
<surname>Kadlubar</surname> <given-names>S</given-names>
</name>
<name>
<surname>Novelli</surname> <given-names>G</given-names>
</name>
<etal/>
</person-group>. <article-title>Oxidized LDL Receptor 1 (OLR1) as a Possible Link Between Obesity, Dyslipidemia and Cancer</article-title>. <source>PloS One</source> (<year>2011</year>) <volume>6</volume>(<issue>5</issue>):<elocation-id>e20277</elocation-id>. doi: <pub-id pub-id-type="doi">10.1371/journal.pone.0020277</pub-id>
</citation>
</ref>
<ref id="B45">
<label>45</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname> <given-names>X</given-names>
</name>
<name>
<surname>Gao</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Qin</surname> <given-names>J</given-names>
</name>
<name>
<surname>Lu</surname> <given-names>S</given-names>
</name>
</person-group>. <article-title>lncRNA MIAT Promotes Proliferation and Invasion of HCC Cells via sponging miR-214</article-title>. <source>Am J Physiol Gastrointestinal liver Physiol</source> (<year>2018</year>) <volume>314</volume>(<issue>5</issue>):<fpage>G559</fpage>&#x2013;<lpage>65</lpage>. doi: <pub-id pub-id-type="doi">10.1152/ajpgi.00242.2017</pub-id>
</citation>
</ref>
<ref id="B46">
<label>46</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sha</surname> <given-names>M</given-names>
</name>
<name>
<surname>Lin</surname> <given-names>M</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>J</given-names>
</name>
<name>
<surname>Ye</surname> <given-names>J</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>J</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>N</given-names>
</name>
<etal/>
</person-group>. <article-title>Long Non-Coding RNA MIAT Promotes Gastric Cancer Growth and Metastasis Through Regulation of miR-141/DDX5 Pathway</article-title>. <source>J Exp Clin Cancer Res CR</source> (<year>2018</year>) <volume>37</volume>(<issue>1</issue>):<fpage>58</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s13046-018-0725-3</pub-id>
</citation>
</ref>
<ref id="B47">
<label>47</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname> <given-names>W</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Z</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>C</given-names>
</name>
<name>
<surname>Ai</surname> <given-names>Z</given-names>
</name>
</person-group>. <article-title>Long Non-Coding RNA MIAT Promotes Papillary Thyroid Cancer Progression Through Upregulating LASP1</article-title>. <source>Cancer Cell Int</source> (<year>2019</year>) <volume>19</volume>:<fpage>194</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s12935-019-0913-z</pub-id>
</citation>
</ref>
<ref id="B48">
<label>48</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kurachi</surname> <given-names>M</given-names>
</name>
<name>
<surname>Barnitz</surname> <given-names>RA</given-names>
</name>
<name>
<surname>Yosef</surname> <given-names>N</given-names>
</name>
<name>
<surname>Odorizzi</surname> <given-names>PM</given-names>
</name>
<name>
<surname>DiIorio</surname> <given-names>MA</given-names>
</name>
<name>
<surname>Lemieux</surname> <given-names>ME</given-names>
</name>
<etal/>
</person-group>. <article-title>The Transcription Factor BATF Operates as an Essential Differentiation Checkpoint in Early Effector CD8+ T Cells</article-title>. <source>Nat Immunol</source> (<year>2014</year>) <volume>15</volume>(<issue>4</issue>):<page-range>373&#x2013;83</page-range>. doi: <pub-id pub-id-type="doi">10.1038/ni.2834</pub-id>
</citation>
</ref>
<ref id="B49">
<label>49</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dai</surname> <given-names>L</given-names>
</name>
<name>
<surname>Cui</surname> <given-names>X</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>X</given-names>
</name>
<name>
<surname>Cheng</surname> <given-names>L</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>Y</given-names>
</name>
<etal/>
</person-group>. <article-title>SARI Inhibits Angiogenesis and Tumour Growth of Human Colon Cancer Through Directly Targeting Ceruloplasmin</article-title>. <source>Nat Commun</source> (<year>2016</year>) <volume>7</volume>:<fpage>11996</fpage>. doi: <pub-id pub-id-type="doi">10.1038/ncomms11996</pub-id>
</citation>
</ref>
<ref id="B50">
<label>50</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname> <given-names>XP</given-names>
</name>
<name>
<surname>He</surname> <given-names>L</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>QP</given-names>
</name>
<name>
<surname>Zeng</surname> <given-names>XT</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>SQ</given-names>
</name>
</person-group>. <article-title>Baicalein Inhibits Proliferation of Myeloma U266 Cells by Downregulating IKZF1 and IKZF3</article-title>. <source>Med Sci monitor Int Med J Exp Clin Res</source> (<year>2018</year>) <volume>24</volume>:<page-range>2809&#x2013;17</page-range>. doi: <pub-id pub-id-type="doi">10.12659/MSM.907058</pub-id>
</citation>
</ref>
<ref id="B51">
<label>51</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nie</surname> <given-names>H</given-names>
</name>
<name>
<surname>Bu</surname> <given-names>F</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>J</given-names>
</name>
<name>
<surname>Li</surname> <given-names>T</given-names>
</name>
<name>
<surname>Huang</surname> <given-names>J</given-names>
</name>
</person-group>. <article-title>29 Immune-Related Genes Pairs Signature Predict the Prognosis of Cervical Cancer Patients</article-title>. <source>Sci Rep</source> (<year>2020</year>) <volume>10</volume>(<issue>1</issue>):<fpage>14152</fpage>. doi: <pub-id pub-id-type="doi">10.1038/s41598-020-70500-5</pub-id>
</citation>
</ref>
<ref id="B52">
<label>52</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mei</surname> <given-names>J</given-names>
</name>
<name>
<surname>Xing</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Lv</surname> <given-names>J</given-names>
</name>
<name>
<surname>Gu</surname> <given-names>D</given-names>
</name>
<name>
<surname>Pan</surname> <given-names>J</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Y</given-names>
</name>
<etal/>
</person-group>. <article-title>Construction of an Immune-Related Gene Signature for Prediction of Prognosis in Patients With Cervical Cancer</article-title>. <source>Int Immunopharmacol</source> (<year>2020</year>) <volume>88</volume>:<fpage>106882</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.intimp.2020.106882</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>