Journal of Immunology Research

Journal of Immunology Research / 2021 / Article
Special Issue

Cancer Immunotherapy: from Basic Research to Clinical Translations

View this Special Issue

Research Article | Open Access

Volume 2021 |Article ID 9960905 | https://doi.org/10.1155/2021/9960905

Qi Pan, Ying Cheng, Donghua Cheng, "Identification of CD8+ T Cell-Related Genes: Correlations with Immune Phenotypes and Outcomes of Liver Cancer", Journal of Immunology Research, vol. 2021, Article ID 9960905, 17 pages, 2021. https://doi.org/10.1155/2021/9960905

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

Academic Editor: Sainan Li
Received13 Mar 2021
Revised01 May 2021
Accepted18 May 2021
Published03 Jun 2021

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% H2O2 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.

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 (, ).

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.


IDModule colorGS.T.cells.CD8P.GS.T.cells.CD8

PSMB9Grey600.42727727
PSMB10Grey600.3960461
PSMB8Grey600.392423292
HLA-EGrey600.392387821
GBP4Grey600.38525836
RARRES3Grey600.379808931
CD74Grey600.379010979
HLA-DQB1Grey600.378859742
PSME2Grey600.378610599
GBP1Grey600.377173717
HLA-DPB1Grey600.376684402
HLA-DRB1Grey600.368162706
HLA-AGrey600.359477552
IRF1Grey600.345346667
HLA-DMAGrey600.337351935
HLA-FGrey600.334500299
UBE2L6Grey600.333949723
HLA-DRAGrey600.328624957
HLA-DPA1Grey600.318999761
HLA-DMBGrey600.315091809

GS: gene significance.

IDModule colorGS.T.cells.CD8P.GS.T.cells.CD8

NKG7Magenta0.626
CST7Magenta0.592
GZMAMagenta0.544
CCL5Magenta0.523
CXCL9Magenta0.470
HCSTMagenta0.463
CD3DMagenta0.421
CD2Magenta0.413
CALHM6Magenta0.406
GMFGMagenta0.405
MZB1Magenta0.397
C1QAMagenta0.395
PIM2Magenta0.389
CD3EMagenta0.384
C1QBMagenta0.375
IGLL5Magenta0.367
CORO1AMagenta0.362
AIF1Magenta0.355
C1QCMagenta0.341
JCHAINMagenta0.333

GS: gene significance.
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).

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.

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).

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)).

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.

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

  1. 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. View at: Publisher Site | Google Scholar
  2. 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. View at: Publisher Site | Google Scholar
  3. 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. View at: Publisher Site | Google Scholar
  4. 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. View at: Publisher Site | Google Scholar
  5. 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. View at: Publisher Site | Google Scholar
  6. 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. View at: Publisher Site | Google Scholar
  7. 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. View at: Publisher Site | Google Scholar
  8. 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. View at: Publisher Site | Google Scholar
  9. 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. View at: Publisher Site | Google Scholar
  10. 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. View at: Publisher Site | Google Scholar
  11. 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. View at: Publisher Site | Google Scholar
  12. 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. View at: Publisher Site | Google Scholar
  13. 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. View at: Publisher Site | Google Scholar
  14. 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. View at: Publisher Site | Google Scholar
  15. 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. View at: Publisher Site | Google Scholar
  16. P. Langfelder and S. Horvath, “WGCNA: an R package for weighted correlation network analysis,” BMC Bioinformatics, vol. 9, no. 1, p. 559, 2008. View at: Publisher Site | Google Scholar
  17. 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. View at: Publisher Site | Google Scholar
  18. 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. View at: Publisher Site | Google Scholar
  19. 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. View at: Publisher Site | Google Scholar
  20. 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. View at: Publisher Site | Google Scholar
  21. 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. View at: Publisher Site | Google Scholar
  22. 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. View at: Publisher Site | Google Scholar
  23. 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. View at: Publisher Site | Google Scholar
  24. 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. View at: Publisher Site | Google Scholar
  25. 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. View at: Publisher Site | Google Scholar
  26. 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. View at: Publisher Site | Google Scholar
  27. 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. View at: Publisher Site | Google Scholar
  28. 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. View at: Publisher Site | Google Scholar

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.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views379
Downloads406
Citations

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.