Journal of Immunology Research
Volume 2021 (2021), Article ID 9960905, 17 pages
https://doi.org/10.1155/2021/9960905

Research Article

Identification of CD8+ T Cell-Related Genes: Correlations with Immune Phenotypes and Outcomes of Liver Cancer

Qi Pan, Ying Cheng, and Donghua Cheng

Department of Hepatobiliary Surgery and Organ Transplantation, Key Laboratory of Organ Transplantation of Liaoning Province, First Hospital of China Medical University, China

Correspondence should be addressed to Donghua Cheng; dhcheng@cmu.edu.cn

Received 13 March 2021; Revised 1 May 2021; Accepted 18 May 2021; Published 3 June 2021

Academic Editor: Sainan Li

Copyright © 2021 Qi Pan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Purpose. Treatment outcomes for advanced liver cancer are poor. Immunotherapy is a treatment strategy that has been widely used to treat other cancers. Studies have shown that CD8+ T lymphocytes are essential factors affecting the efficacy of immunotherapy. We used computational biology methods to determine the coexpressed gene network that promotes CD8+ T lymphocyte infiltration. Method. We obtained the liver cancer gene matrix and clinical follow-up information data from TCGA liver hepatocellular carcinoma FPKM. We obtained single nucleotide polymorphism (SNP) data to evaluate the tumor mutation burden. The “estimate” package and the CIBERSORT algorithm were used to evaluate tumor purity and the proportion of CD8+ T lymphocytes in the liver cancer cohort. We used the gene expression matrix of liver cancer and the relative proportion of CD8+ T lymphocytes as input files and performed WGCNA based on this analysis. The weighted coexpression network identified the most CD8+ T lymphocyte-related coexpression modules in liver cancer. Then, we analyzed the biological processes involved in the module. We determined the coexpression module with CD8+ T lymphocyte infiltration in terms of data and function. We then screened the factors in the coexpression module correlated with CD8+ T lymphocyte content greater than 0.4. Finally, the expression levels of these factors were verified at the protein level using immunohistochemistry and single-cell sequencing. Results. We determined the CD8+ T lymphocyte proportions that correlated with coexpression networks. Four coexpressed genes (C1QC, CD3D, GZMA, and PSMB9) were identified as CD8+ T cell coexpression genes that promoted infiltration of CD8+ T cells. Because the factors in the coexpression network often participate in similar biological processes, we found that these factors were most related to antigen processing and presentation of peptide antigen through functional enrichment. In the clinical phenotype analysis, we found that 18 factors can be used as independent prognostic protective factors. We found that these factors were significantly negatively correlated with tumor purity and negatively correlated with M2 macrophages in the immunophenotyping analysis. Using immunohistochemistry and single-cell sequencing analysis, we found that CD3D antibody staining was weaker in tumor tissues than normal tissues and was related to CD8+ T cells. Conclusion. These coexpressed genes were positively related to the high infiltration proportion of CD8+ T lymphocytes in an antigen presentation process. The biological process might provide new directions for patients who are insensitive to immune therapy.

1. Introduction

In recent years, breakthroughs have been made in using immune checkpoint inhibitors [1], ushering in a new era to treat advanced tumors. The emergence of immunotherapy provides options for the treatment of liver cancer; these include immune checkpoint inhibitors, adoptive cell transfer, tumor vaccines, and cytokine therapy. Immune checkpoint inhibitors enhance antitumor immune responses by reversing the exhaustion of T cell function and restoring immune recognition and immune attack. Immune checkpoint inhibitor targets include programmed death-ligand 1 (PD-L1) and its receptor PD-1 (programmed cell death protein 1) and cytotoxic T lymphocyte-related antigen 4. PD-1 is a member of the CD28 family and is expressed on the surface of most immune cells, mainly on CD8+ T cells [2]. It binds to PD-L1 and PD-L2 to cause inhibitory signals to be transmitted to T cells and induce tolerance [3]. PD-L1 is abnormally expressed in various tumors, including liver cancer; tumor cells achieve immune escape by abnormally expressing PD-L1 or PD-L2 [4]. Nivolumab (PD-1 monoclonal antibody) [5], pembrolizumab (PD-1 monoclonal antibody), and atezolizumab (PD-L1 monoclonal antibody) have been used to treat advanced liver cancer. Although immunotherapy has achieved gratifying results and has been approved as a treatment option for liver cancer, the objective response rate of the PD-1/PD-L1 antibody alone rarely exceeds 40%. The objective response rate of nivolumab and pembrolizumab in liver cancer was not significantly effective [6]. It is currently believed that the main reason for the low objective response rate of immunotherapy is the emergence of drug resistance [7]. The drug resistance of immunotherapy is a complex and multimechanism interdependent dynamic process, explained by impaired immune infiltration in tumors, depletion of T cells, or recruitment of immunosuppressive cells.

Studies have shown that low PD-1 expression in tumor tissues and low CD8+ T lymphocyte infiltration can cause this insensitivity [8]. According to the results of a meta-analysis of non-small-cell lung cancer, increased numbers of CD8+ tumor-infiltrating lymphocytes are associated with better overall survival [9]. In patients with advanced melanoma treated with pembrolizumab, the densities of CD8+ T cells in the invasion margin and tumor center of the tissue specimens of responders were higher than those of nonresponders [10].

Therefore, this article intends to identify the coexpression modules and functions related to the content of CD8+ T lymphocytes in the tumor microenvironment and define the coexpression network related to the content of CD8+ T cells, to provide a basis for improving the efficacy of immunotherapy.

2. Methods

2.1. Data Collection and Preprocessing

TCGA expression matrix data from LIHC samples were downloaded from The Cancer Genome Atlas (http://cancergenome.nih.gov/). Tumor transcriptomic profiles of 19530 mRNA were measured in 377 liver hepatocellular carcinoma patients and were brought into the subsequent analysis. The hepatocellular carcinoma single-cell mice sequencing cohort GSE129516 [11] was also downloaded from the GEO database, whose platform is GPL24247.

2.2. CD8+ T Cell Proportion Evaluation

To obtain the relative proportion of CD8+ T lymphocytes in each liver hepatocellular carcinoma sample, we used the CIBERSORT [12] method. CIBERSORT evaluated the proportion of 22 tumor-infiltrating immune cells in each sample. The samples with were brought into the WGCNA.

2.3. Tumor Microenvironment Score and Tumor Mutation Burden (TMB)

The estimation of stromal and immune cells in malignant tumor tissues [13] is a method that evaluates the proportion of stromal and immune cells by gene expression signatures [13]. To perform a correlation analysis with the tumor mutation burden [14, 15], we obtained TCGA-LIHC SNP data and evaluated the tumor mutation burden of each sample.

2.4. Weighted Gene Coexpression Network Analysis (WGCNA)

We used WGCNA [16] to explore the correlations between gene expression and CD8+ T lymphocytes. First, we matched the mRNA matrix of TCGA-LIHC and the corresponding CD8+ T lymphocyte content of the sample and included them in the subsequent WGCNA. To ensure the nonscale of the network, we built a scale-free topology network [17, 18] and set the soft threshold at 5 (, ) and the number of genes in the minimum module at 30.

2.5. Gene Ontology Functional Analysis

Gene ontology (GO) [19] analysis was performed to show the biological processes and molecular functions based on the different modules. In this study, we performed GO analysis based on the clusterProfiler [20] package in R3.6.2.

2.6. Cox Hazard Proportion Regression Model and Subgroup Evaluation

We applied the Cox risk proportional regression model for the factors in the grey module. The prognostic risk model related to CD8+ T lymphocytes in liver cancer was determined using the forward screening method. We used survival analysis and the area under the curve (AUC) to evaluate the accuracy of the prognostic risk model after rain. Next, we divided TCGA-LIHC cohort into various subgroups according to age, gender, clinical stage, and tumor purity and evaluated the CD8+ T lymphocyte-related risk models in various subgroups.

2.7. Gene Set Enrichment Analysis

Gene set enrichment analysis (GSEA) [21] was used to calculate the most involved pathway to these coexpression genes.

2.8. Immunohistochemistry

The extracted human tissues were fixed with a 4% formaldehyde buffer. Deparaffinized specimens were then sectioned into 4 μm slices. Tissue slices were incubated at 60°C for 2 h before dewaxing; the sections were autoclaved at 115°C for 3 min for antigen retrieval in a citric acid buffer (pH 6.0) and quenched for endogenous peroxidase activity with 0.3% H₂O₂ solution for 15 min. Then, the slices were blocked for nonspecific binding with normal goat serum for 45 min and incubated with the specific primary antibody against C1QC, GZMA, CD3D, and PSMB9 (dilution 1 : 200) overnight at 4°C. Subsequently, the sections were treated with the goat anti-mouse secondary antibody for 30 min at room temperature. Protein expression was visualized using 3,3-diaminobenzidine. Images were captured using a Nikon Eclipse 80i microscope (Nikon Corporation).

3. Results

A flow chart is displayed in Figure 1, which illustrates the logic of our analysis.

Figure 1: Flow chart of this paper.

3.1. CD8+ T Lymphocyte, Tumor Purity, and Tumor Mutation Burden Evaluation

We measured the CD8+ T lymphocyte proportion, tumor purity, matrix score, immune score, and tumor mutation burden of each LIHC sample. Using the screening principle of , we obtained 377 liver cancer samples accurately evaluated by CD8+ T lymphocytes. By integrating the immune-related file with TCGA-LIHC mRNA expression files, we determined WGCNA phenotype entry files.

3.2. CD8+ T Lymphocyte Coexpression Network Conduction in TCGA-LIHC

WGCNA was performed on TCGA liver hepatocellular carcinomas. To construct a CD8+ T lymphocyte coexpression network, we used a dynamic hybrid cutting method to construct a hierarchical clustering tree (Figure 2(a)). Each leaf on the tree represents a gene, and each branch represents a coexpression module. A total of 22 expression modules were obtained (Figure 2(b)). Then, we calculated the correlation coefficient between each module and CD8+ T lymphocytes and selected the grey60 and purple modules according to the correlation coefficient (Figure 2(c)). The grey60 module had the strongest correlation with CD8+ T lymphocyte proportion in TCGA-LIHC cohort (; ) (Figure 2(c)). The magenta module correlated with the CD8+ T lymphocyte proportion in TCGA-LIHC cohort (; ) (Figure 2(c)). Based on these findings, we supplemented the heat map of the correlation between the factors in the grey60 and magenta modules (Figures 2(d)–2(g)). The grey60 module showed a significant correlation with CD8+ T cells (, ) and tumor purity (, ). The magenta module showed a significant correlation with CD8+ T cells (, ) and tumor purity (, ).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 2: WGCNA. (a) TCGA-LIHC sample clustering results. (b) Twenty-two coexpression modules were conducted by WGCNA. (c) Correlation heat map among different coexpression modules and clinical phenotypes and CD8+ T levels. (d, e) The correlation of genes in grey60 modules with CD8+ T cells and tumor purity. (f, g) The correlation of genes in magenta modules with CD8+ T cells and tumor purity.

3.3. CD8+ T Lymphocyte Coexpression Module Functional Analysis

We determined the top 20 CD8+ T lymphocyte proportions positively coexpressing mRNA in TCGA-LIHC grey60 and magenta modules (Tables 1 and 2). The 20 CD8+ T lymphocyte proportions positively coexpressing mRNA in the magenta module were most significantly enriched in antigen processing and presentation. The 20 CD8+ T lymphocyte proportions positively coexpressing mRNA in the grey60 module were most significantly enriched in regulating lymphocyte activation, suggesting that these biological processes might promote CD8+ T lymphocyte infiltration in the liver hepatocellular cancer microenvironment (Figures 3(a) and 3(b)). The protein-protein interaction network of yellow and green modules is shown in Figure 3.

Table 1: The module and gene significance for CD8+ T cell-related genes in the grey60 module.

ID
Module color
GS.T.cells.CD8
P.GS.T.cells.CD8

PSMB9
Grey60
0.42727727

PSMB10
Grey60
0.3960461

PSMB8
Grey60
0.392423292

HLA-E
Grey60
0.392387821

GBP4
Grey60
0.38525836

RARRES3
Grey60
0.379808931

CD74
Grey60
0.379010979

HLA-DQB1
Grey60
0.378859742

PSME2
Grey60
0.378610599

GBP1
Grey60
0.377173717

HLA-DPB1
Grey60
0.376684402

HLA-DRB1
Grey60
0.368162706

HLA-A
Grey60
0.359477552

IRF1
Grey60
0.345346667

HLA-DMA
Grey60
0.337351935

HLA-F
Grey60
0.334500299

UBE2L6
Grey60
0.333949723

HLA-DRA
Grey60
0.328624957

HLA-DPA1
Grey60
0.318999761

HLA-DMB
Grey60
0.315091809

GS: gene significance.

Table 2: The module and gene significance for CD8+ T cell-related genes in the magenta module.

ID
Module color
GS.T.cells.CD8
P.GS.T.cells.CD8

NKG7
Magenta
0.626

CST7
Magenta
0.592

GZMA
Magenta
0.544

CCL5
Magenta
0.523

CXCL9
Magenta
0.470

HCST
Magenta
0.463

CD3D
Magenta
0.421

CD2
Magenta
0.413

CALHM6
Magenta
0.406

GMFG
Magenta
0.405

MZB1
Magenta
0.397

C1QA
Magenta
0.395

PIM2
Magenta
0.389

CD3E
Magenta
0.384

C1QB
Magenta
0.375

IGLL5
Magenta
0.367

CORO1A
Magenta
0.362

AIF1
Magenta
0.355

C1QC
Magenta
0.341

JCHAIN
Magenta
0.333

GS: gene significance.

(a)

(b)

Figure 3: Functional analysis and protein-protein network of grey60 and magenta modules. (a) Grey60 module. (b) Magenta module.

3.4. Clinical Phenotype and Immune Phenotype Correlation of Coexpression Genes

To determine the outcome status correlation of these coexpression genes, we performed survival analysis. The patients in low expression groups for CCL5 (TCGA: ), CST7 (TCGA: ), HLA-DPA1 (TCGA: ), HLA-E (TCGA: ), NKG7 (TCGA: ), CD2 (TCGA: ), GBP1 (TCGA: ), HLA-DPB1 (TCGA: ), IGLL5 (TCGA: ), HLA-DRB1 (TCGA: ), CD3E (TCGA: ), GZMA (TCGA: ), HLA-DRA (TCGA: ), JCHAIN (TCGA: ), MZB1 (TCGA: ), CORO1A (TCGA: ), and HLA-DMA (TCGA: ) showed survival risk against high expression groups (Figure 4). These results suggest that these coexpression genes in grey60 and magenta modules protect against liver hepatocellular cancer. Next, we found that these factors were positively correlated with CD8+ T lymphocytes (Figure 5) and negatively correlated with tumor purity (Figure 6).

Figure 4: The survival analysis of genes in magenta and grey60 modules. The grey60 module genes acted as prognostic protective genes.

Figure 5: The Pearson correlation test among genes in magenta and grey60 modules with CD8+ T cell proportions.

Figure 6: The Pearson correlation test among genes in magenta and grey60 modules with tumor purity.

3.5. Cox Regression Hazard Model of CD8+ T Lymphocyte Coexpression Genes

A CD8+ T lymphocyte coexpression gene Cox regression hazard model was conducted based on these liver hepatocellular outcome protective factors.

The samples in high-risk samples for liver hepatocellular cancer patients (TCGA: ; ) (Figure 7) showed survival risk against low-risk groups, with (Figure 7). The risk score was evaluated in various subgroups, including age, gender, stage, tumor purity, and tumor mutation burden. The results were significant in these subgroups.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

(m)

(n)

Figure 7: Cox regression analysis based on these CD8+ T cell-related genes. (a) The patients in the higher risk group showed worse outcome levels. (b) The receiver operating characteristic curve of the risk score in TCGA-LIHC cohort. (c–n) The subgroup survival analysis of CD8+ T cell-related gene risk scores.

3.6. GSEA

Antigen processing and presentation, the chemokine signaling pathway, the cytokine-cytokine receptor interaction pathway, and the T cell receptor signaling pathway were related to the high expression group in C1QC, CD3D, GZMA, and PSMB9 (Figure 8).

Figure 8: GSEA.

3.7. Tissue Verification and Single-Cell Markers

Immunohistochemical analysis of the protein expression levels of CD3D in and around carcinomas was performed in the liver cancer cohorts of the First Affiliated Hospital of China Medical University. The results showed that CD3D had lower staining intensity (Figure 9(a)). We also performed UMAP dimensional reduction clustering in the data of a single cell of liver cancer in GEO. After annotating the subsets with “SingleR,” we obtained the subsets containing T cell macrophages and other cells. We found that the expression content of CD3D was relatively high in the T cell subsets, which confirmed our previous conclusion in TCGA-LIHC (Figure 9(b)).

(a)

(b)

Figure 9: (a) Immunohistochemical analysis in tumor and paracarcinoma tissues from China Medical University clinical samples. (b) Single-cell cohort validation in GSE129516.

We analyzed the relationship between CD3D and currently known gene sets in the tumor microenvironment and verified the expression association between CD3D and CD8A in two other cohorts, GSE29721 [22] and GSE121248 [23]. We also marked the distribution of different T cell subtypes in single-cell cohorts, and the results showed that CD3D was more strongly associated with CD8A and less with CD4+ T lymphocytes (Supplementary Figure 1).

4. Discussion

Immunotherapy is a late-stage cancer treatment strategy widely used in clinical practice; nevertheless, there are many unsatisfactory results, including the low success rate of immunotherapy, complications associated with immunotherapy, and the super progression of immunotherapy. Studies have shown that the low success rate of immunotherapy is related to the low infiltration of CD8+ T lymphocytes, the weakening of the antigen presentation process, and the reduction of PD-1 expression. In this study, we explored the coexpression network that promotes the infiltration of CD8+ T lymphocytes in liver cancer by combining computational biology and experiments. The factors in the coexpression network are considered to have similar biological functions in organisms. Therefore, mining these coexpression factors helps us understand the biological processes most closely related to the infiltration of CD8+ T lymphocytes in liver cancer.

The grey60 module was the most relevant coexpression module for CD8+ T lymphocytes. The factors in this module are related to the regulation of T cell proliferation. To date, we have demonstrated that these factors are related to the infiltration of CD8+ T lymphocytes at the sequencing and functional levels. Based on our inferences, we believe that these factors may improve outcomes by increasing the infiltration of CD8+ T lymphocytes. Then, we conducted survival analysis on these factors and successfully constructed a prognostic risk score for liver cancer based on CD8+ T lymphocyte coexpression factors. The differences in protein levels of C1QC, CD3D, GZMA, and PSMB9 in different liver cancer stages were determined using immunohistochemistry.

C1QC encodes the complement C1q C chain, which associates with C1r and C1s to yield the first component of the serum complement system [24, 25]. C1q is composed of 18 polypeptide chains which include six A-chains, six B-chains, and six C-chains. Each chain contains an N-terminal collagen-like region and a C-terminal C1q globular domain. The protein encoded by CD3D is part of the T cell receptor/CD3 complex (TCR/CD3 complex) and is involved in T cell development and signal transduction [26, 27]. The encoded membrane protein represents the delta subunit of the CD3 complex. Along with four other CD3 subunits, the encoded membrane protein binds either TCR alpha/beta or TCR gamma/delta to form the TCR/CD3 complex on the surface of T cells [28]. In this study, transcriptome, histological, and single-cell cohorts were used to demonstrate the importance of C1QC and CD3D in CD8+ T lymphocytes.

This article has some limitations. We used TCGA and liver cancer tissues from China Medical University to conduct a joint analysis. More external cohorts still need to be cross-validated. Due to the limited computing power of computers, we only included factors with variances in the top 25% when performing WGCNA. This may cause some false-negative results. More advanced computers need to be used to repeat this screening method. Although we hypothesize that these factors can improve the therapeutic effect of liver cancer immunotherapy by promoting the infiltration of CD8+ T lymphocytes, due to the lack of immunotherapy efficacy evaluation in the follow-up data of TCGA, we only explored the passage of these coexpression factors in liver cancer. More cohorts with immune follow-up data need to be added to promote the infiltration of CD8+ T lymphocytes and improve the outcome of patients.

In summary, we constructed a prognostic risk scoring model for liver cancer based on a CD8+ T lymphocyte content coexpression molecular network, determined that the factors in the risk scoring model can be used as independent prognostic factors for liver cancer, and determined the levels of these factors at the mRNA and protein levels. These prognosis-related CD8+ T lymphocyte coexpression factors and their related biological functions may provide new directions for improving the efficacy of immunotherapy.

Data Availability

The datasets TCGA-LIHC for this study can be found at http://cancergenome.nih.gov/. The datasets GSE129516, GSE29721, and GSE121248 in this study can be found at http://www.ncbi.nlm.nih.gov/geo/.

Ethical Approval

The ethical number was China Medical University AF-0G-03-1.0-02.

Consent

All the patients signed the informed consent.

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

Qi Pan wrote the manuscript. Ying Cheng and Donghua Cheng revised the manuscript.

Acknowledgments

This work was supported by the Liaoning Provincial Department of Education’s scientific research project in 2020 (FWZR2020002). We thank the authors of those references for their excellent work. A preprint version has been uploaded in https://www.researchsquare.com/article/rs-146843/v1.

Supplementary Materials

Supplementary Figure 1: CD3D immune environment analysis and verification. (A) CD3D expression positively correlates to the intensity immune response. (B) CD3D expression positively correlates to the CD8A in GSE29721 and GSE121248. (C) The distribution of CD8A and CD4. (Supplementary Materials)

References

J. M. Kirkwood, L. H. Butterfield, A. A. Tarhini, H. Zarour, P. Kalinski, and S. Ferrone, “Immunotherapy of cancer in 2012,” CA: a Cancer Journal for Clinicians, vol. 62, no. 5, pp. 309–335, 2012.

F. Xu, T. Jin, Y. Zhu, and C. Dai, “Immune checkpoint therapy in liver cancer,” Journal of Experimental & Clinical Cancer Research : CR, vol. 37, no. 1, p. 110, 2018.

T. Hato, L. Goyal, T. F. Greten, D. G. Duda, and A. X. Zhu, “Immune checkpoint blockade in hepatocellular carcinoma: current progress and future directions,” Hepatology, vol. 60, no. 5, pp. 1776–1782, 2014.

T. Okazaki and T. Honjo, “PD-1 and PD-1 ligands: from discovery to clinical application,” International Immunology, vol. 19, no. 7, pp. 813–824, 2007.

A. B. el-Khoueiry, B. Sangro, T. Yau et al., “Nivolumab in patients with advanced hepatocellular carcinoma (CheckMate 040): an open-label, non-comparative, phase 1/2 dose escalation and expansion trial,” Lancet, vol. 389, no. 10088, pp. 2492–2502, 2017.

B. Sangro, P. Sarobe, S. Hervás-Stubbs, and I. Melero, “Advances in immunotherapy for hepatocellular carcinoma,” Nature Reviews. Gastroenterology & Hepatology, pp. 1–19, 2021.

P. Sharma, S. Hu-Lieskovan, J. A. Wargo, and A. Ribas, “Primary, adaptive, and acquired resistance to cancer immunotherapy,” Cell, vol. 168, no. 4, pp. 707–723, 2017.

D. J. Byun, J. D. Wolchok, L. M. Rosenberg, and M. Girotra, “Cancer immunotherapy - immune checkpoint blockade and associated endocrinopathies,” Nature Reviews. Endocrinology, vol. 13, no. 4, pp. 195–207, 2017.

Y. Yu, D. Zeng, Q. Ou et al., “Association of survival and immune-related biomarkers with immunotherapy in patients with non-small cell lung cancer: a meta-analysis and individual patient-level analysis,” JAMA Network Open, vol. 2, no. 7, article e196879, 2019.

P. C. Tumeh, C. L. Harview, J. H. Yearley et al., “PD-1 blockade induces responses by inhibiting adaptive immune resistance,” Nature, vol. 515, no. 7528, pp. 568–571, 2014.

X. Xiong, H. Kuang, S. Ansari et al., “Landscape of intercellular crosstalk in healthy and NASH liver revealed by single-cell secretome gene analysis,” Molecular Cell, vol. 75, no. 3, pp. 644–660.e5, 2019.

B. Chen, M. S. Khodadoust, C. L. Liu, A. M. Newman, and A. A. Alizadeh, “Profiling tumor infiltrating immune cells with CIBERSORT,” Methods in Molecular Biology, vol. 1711, pp. 243–259, 2018.

K. Yoshihara, M. Shahmoradgoli, E. Martínez et al., “Inferring tumour purity and stromal and immune cell admixture from expression data,” Nature Communications, vol. 4, no. 1, p. 2612, 2013.

Y. Li, Z. Chen, L. Wu, and W. Tao, “Novel tumor mutation score versus tumor mutation burden in predicting survival after immunotherapy in pan-cancer patients from the MSK-IMPACT cohort,” Annals of Translational Medicine, vol. 8, no. 7, p. 446, 2020.

Z. Yang, S. Wei, Y. Deng, Z. Wang, and L. Liu, “Clinical significance of tumour mutation burden in immunotherapy across multiple cancer types: an individual meta-analysis,” Japanese Journal of Clinical Oncology, vol. 50, no. 9, pp. 1023–1031, 2020.

P. Langfelder and S. Horvath, “WGCNA: an R package for weighted correlation network analysis,” BMC Bioinformatics, vol. 9, no. 1, p. 559, 2008.

J. Jiang, X. Sun, W. Wu et al., “Correction: Corrigendum: Construction and application of a co-expression network in Mycobacterium tuberculosis,” Scientific Reports, vol. 7, no. 1, article 40563, 2017.

J. Jiang, X. Sun, W. Wu et al., “Construction and application of a co-expression network in Mycobacterium tuberculosis,” Scientific Reports, vol. 6, no. 1, article 28422, 2016.

M. Ashburner, C. A. Ball, J. A. Blake et al., “Gene ontology: tool for the unification of biology. The Gene Ontology Consortium,” Nature Genetics, vol. 25, no. 1, pp. 25–29, 2000.

G. Yu, L. G. Wang, Y. Han, and Q. Y. He, “clusterProfiler: an R package for comparing biological themes among gene clusters,” Omics : a journal of integrative biology, vol. 16, no. 5, pp. 284–287, 2012.

A. Subramanian, P. Tamayo, V. K. Mootha et al., “Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 43, pp. 15545–15550, 2005.

B. Stefanska, J. Huang, B. Bhattacharyya et al., “Definition of the landscape of promoter DNA hypomethylation in liver cancer,” Cancer Research, vol. 71, no. 17, pp. 5891–5903, 2011.

S. M. Wang, L. L. P. J. Ooi, and K. M. Hui, “Identification and validation of a novel gene signature associated with the recurrence of human hepatocellular carcinoma,” Clinical Cancer Research, vol. 13, no. 21, pp. 6275–6283, 2007.

K. K. Hosszu, F. Santiago-Schwarz, E. I. B. Peerschke, and B. Ghebrehiwet, “Evidence that a C1q/C1qR system regulates monocyte-derived dendritic cell differentiation at the interface of innate and acquired immunity,” Innate Immunity, vol. 16, no. 2, pp. 115–127, 2010.

Y. Qiu, L. A. Marquez-Curtis, and A. Janowska-Wieczorek, “Mesenchymal stromal cells derived from umbilical cord blood migrate in response to complement C1q,” Cytotherapy, vol. 14, no. 3, pp. 285–295, 2012.

J. Gil, E. M. Busto, B. Garcillán et al., “A leaky mutation in CD3D differentially affects αβ and γδ T cells and leads to a Tαβ-Tγδ+B+NK+ human SCID,” The Journal of Clinical Investigation, vol. 121, no. 10, pp. 3872–3876, 2011.

J. Dagher, B. Delahunt, N. Rioux-Leclercq et al., “Clear cell renal cell carcinoma: validation of World Health Organization/International Society of Urological Pathology grading,” Histopathology, vol. 71, no. 6, pp. 918–925, 2017.

N. Manolios, O. Kemp, and Z. G. Li, “The T cell antigen receptor alpha and beta chains interact via distinct regions with CD3 chains,” European Journal of Immunology, vol. 24, no. 1, pp. 84–92, 1994.