Genetic associations of T cell cancer immune response-related genes with T cell phenotypes and clinical outcomes of early-stage lung cancer

Germline variants of the genes involved in NF-kB activation are associated with the risk of COPD and lung cancer development

Chronic obstructive pulmonary disease (COPD) and lung cancer (LC) are closely related diseases associated with smoking history and dysregulated immune response. However, not all smokers develop the disease, indicating that genetic susceptibility could be important. Therefore, the aim of this study was to search for the potential overlapping genetic biomarkers, with a focus on single nucleotide polymorphisms (SNPs) located in the regulatory regions of immune-related genes. Additionally, the aim was to see if an identified SNP has potentially an effect on proinflamma-tory cytokine concentration in the serum of COPD patients. We extracted summary data of variants in 1511 immune-related genes from COPD and LC genome-wide association studies (GWAS) from the UK Biobank. The LC data had 203 cases, patients diagnosed with LC, and 360 938 controls, while COPD data had 1 897 cases and 359 297 controls. Assuming 1 association/gene, SNPs with a p-value < 3.3 × 10^-5 were considered statistically significantly associated with the disease. We identified seven SNPs located in different genes (BAG6, BTNL2, TNF, HCP5, MICB, NCR3, ABCF1, TCF7L1) to be associated with the COPD risk and two with the LC risk (HLA-C, HLA-B), with statistical significance. We also identified two SNPs located in the IL2RA gene associated with LC (rs2386841; p = 1.86 × 10^-4) and COPD (rs11256442; p = 9.79 × 10^-3) but with lower significance. Functional studies conducted on COPD patients showed that RNA expression of IL2RA, IFNγ and related proinflammatory cytokines in blood serum did not correlate with a specific genotype. Although results presented in this study do not fully support our hypothesis, it is worth to mention that the identified genes/SNPs that were associated with either COPD or LC risk, all were involved in the activation of the NF-κB transcription factor which is closely related to the regulation of the inflammatory response, a condition associated with both pathologies.

Polygenic Risk Scores Associated with Tumor Immune Infiltration in Common Cancers

It is largely unknown whether genetic susceptibility contributes to tumor immune infiltration in common cancers. We systematically investigated the association between polygenic risk scores (PRSs) and tumor immune infiltration in common cancers. First, we constructed a PRS for common cancers using the risk variants identified in previous genome-wide association studies. Then, we analyzed 139 immune traits predicted by previous studies by examining gene expression data in tumor tissues from The Cancer Genome Atlas (TCGA). We applied regression analyses to evaluate the associations between PRS and immune traits for each cancer overall and stratified by stage, including 2160 pathologically confirmed cases of breast, colorectal, lung, ovarian, pancreatic, and prostate cancers in the White population. At a nominal (p < 0.05) significance level, we identified 31 significant associations between PRS and immune traits. In the analyses stratified by stage for breast, colorectal, lung adenocarcinoma, and lung squamous cell carcinoma, we identified 65 significant associations, including 56 associations that were undetected by the overall analysis. This study provides evidence for genetic risk factors affecting immune infiltration and provides novel insights into the role of genetic susceptibility in immune responses, underlying cancer development, prognosis, and the potential role of an early diagnostic or therapeutic targeting strategy.

Evolutionary Origins of Metabolic Reprogramming in Cancer

Metabolic changes that facilitate tumor growth are one of the hallmarks of cancer. These changes are not specific to tumors but also take place during the physiological growth of tissues. Indeed, the cellular and tissue mechanisms present in the tumor have their physiological counterpart in the repair of tissue lesions and wound healing. These molecular mechanisms have been acquired during metazoan evolution, first to eliminate the infection of the tissue injury, then to enter an effective regenerative phase. Cancer itself could be considered a phenomenon of antagonistic pleiotropy of the genes involved in effective tissue repair. Cancer and tissue repair are complex traits that share many intermediate phenotypes at the molecular, cellular, and tissue levels, and all of these are integrated within a Systems Biology structure. Complex traits are influenced by a multitude of common genes, each with a weak effect. This polygenic component of complex traits is mainly unknown and so makes up part of the missing heritability. Here, we try to integrate these different perspectives from the point of view of the metabolic changes observed in cancer.

Identification of RHEX as a novel biomarker related to progression and immunity of non-small cell lung carcinoma

Background

The therapeutic response and prognosis of patients with non-small cell lung carcinoma (NSCLC) are widely related to immunity. To improve the prognosis of patients and provide reliable information to guide appropriate personalized treatment strategies, it is necessary to identify reliable prognostic or predictive indicators closely related to tumor phenotype and immune traits in NSCLC.

Methods

Based on The Cancer Genome Atlas (TCGA)-NSCLC mRNA expression profile data, a novel approach combining differential gene expression analysis, single-sample gene set enrichment analysis (ssGSEA), and weighted gene co-expression network analysis (WGCNA) was used to screen hub genes. Subsequently, the regulator of hemoglobinization and erythroid cell expansion (RHEX) was identified as a key gene using the log-rank test and confirmed in the ArrayExpress database. The relationship between RHEX and clinicopathological parameters was analyzed using the Wilcoxon rank-sum test. More importantly, through gene set enrichment analysis (GSEA) and cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT) algorithms, and with reference to the Tumor IMmune Estimation Resource (TIMER) database, we explored the relevant pathways of RHEX and its relationship with tumor-infiltrating immune cells (TICs). Finally, we depicted the association between RHEX and immunomodulators in the TCGA and a web portal TISIDB.

Results

The RHEX mRNA expression levels in tumor tissues were lower than those in normal tissues and declined with the progression of NSCLC. Meanwhile, RHEX overexpression was associated with high immune infiltration levels and a favorable clinical prognosis. RHEX may participate in tumor microenvironment (TME) regulation through multiple tumor-immune related pathways, especially the JAK-STAT signaling pathway. Furthermore, RHEX expression affected the infiltrating abundance of multiple TICs and positively correlated with most of the immunomodulators in NSCLC.

Conclusions

Our study is the first to propose that RHEX is an immune-related gene with prognostic value in NSCLC and reveals the underlying mechanism between RHEX and tumor-immune system interactions. These results ultimately provide guidance for prognosis and immunotherapy for NSCLC patients.

Programmed Death Ligand 2 Gene Polymorphisms Are Associated With Lung Adenocarcinoma Risk in Female Never-Smokers

Objectives

Lung cancer in never-smokers is a distinct disease associated with a different genomic landscape, pathogenesis, risk factors, and immune checkpoint inhibitor responses compared to those observed in smokers. This study aimed to identify novel single nucleotide polymorphisms (SNPs) of programmed death-1 (encoded by PDCD1) and its ligands, programmed death ligand 1 (CD274) and 2 (PDCD1LG2), associated with lung cancer risk in never-smoking women.

Materials and Methods

During September 2002 and July 2012, we enrolled never-smoking female patients with lung adenocarcinoma (LUAD) (n=1153) and healthy women (n=1022) from six tertiary hospitals in Taiwan. SNP data were obtained and analyzed from the genome-wide association study dataset and through an imputation method. The expression quantitative trait loci (eQTL) analysis was performed in both tumor and non-tumor tissues for the correlation between genetic expression and identified SNPs.

Results

A total of 12 PDCD1LG2 SNPs related to LUAD risk were identified in never-smoking women, including rs2381282, rs4742103, rs4237162, rs4742104, rs12237624, rs78096119, rs6476988, rs7857315, rs10975178, rs7854413, rs56001683, and rs7858319. Among them, six tagged PDCD1LG2 SNPs rs2381282, rs4742103, rs4237162, rs4742104, rs78096119, and rs56001683 were significantly associated with LUAD risk. Specifically, two PDCD1LG2 SNPs, rs12237624 and rs78096119, were associated with previous pulmonary tuberculosis infection in relation to LUAD susceptibility. Through an eQTL assay, we found that rs2381282 (p < 0.001), rs12237624 (p = 0.019), and rs78096119 (p = 0.019) were associated with the expression levels of programed death ligand 2.

Conclusions

Novel SNPs of programed death ligand 2 associated with lung adenocarcinoma risk were identified. Among them, two SNPs were associated with pulmonary tuberculosis infection in relation to lung adenocarcinoma susceptibility. These SNPs may help to stratify high-risk populations of never-smokers during lung cancer screening.

Clinical Relevance of Serum Kyn/Trp Ratio and Basal and IFNγ-Upregulated IDO1 Expression in Peripheral Monocytes in Early Stage Melanoma

Immune escape is an early phenomenon in cancer development/progression. Indoleamine 2,3-dioxygenase 1 (IDO1) is a normal endogenous mechanism of acquired peripheral immune tolerance and may therefore be tumor-promoting. This study investigated the clinical relevance of IDO1 expression by immune cells in the lymph nodes and blood and of the serum kynurenine/tryptophan (Kyn/Trp) ratio in 65 systemic treatment naïve stage I-III melanoma patients. Blood samples were collected within the first year of diagnosis. Patients had a median follow-up of 61 months. High basal IDO1 expression in peripheral monocytes and low IFNγ-induced IDO1 upregulation correlated with worse outcome independent from disease stage. Interestingly studied factors were not interrelated. During follow-up, the risk of relapse was 9% (2/22) in the subgroup with high IFNγ-induced IDO1 upregulation in monocytes. In contrast, if IDO1 upregulation was low, relapse occurred in 30% (3/10) of patients with low basal IDO1 expression in monocytes and in 61.5% (8/13) in the subgroup with high basal IDO1 expression in monocytes (Log-Rank test, p=0.008). This study reveals some immune features in the blood of early stage melanoma that may be of relevance for disease outcome. These may offer a target for sub-stratification and early intervention.

How to Better Understand the Influence of Host Genetics on Developing an Effective Immune Response to Thoracic Cancers

Thoracic cancers pose a significant global health burden. Immune checkpoint blockade therapies have improved treatment outcomes, but durable responses remain limited. Understanding how the host immune system interacts with a developing tumor is essential for the rational development of improved treatments for thoracic malignancies. Recent technical advances have improved our understanding of the mutational burden of cancer cells and changes in cancer-specific gene expression, providing a detailed understanding of the complex biology underpinning tumor-host interactions. While there has been much focus on the genetic alterations associated with cancer cells and how they may impact treatment outcomes, how host genetics affects cancer development is also critical and will greatly determine treatment response. Genome-wide association studies (GWAS) have identified genetic variants associated with cancer predisposition. This approach has successfully identified host genetic risk factors associated with common thoracic cancers like lung cancer, but is less effective for rare cancers like malignant mesothelioma. To assess how host genetics impacts rare thoracic cancers, we used the Collaborative Cross (CC); a powerful murine genetic resource designed to maximize genetic diversity and rapidly identify genes associated with any biological trait. We are using the CC in conjunction with our asbestos-induced MexTAg mouse model, to identify host genes associated with mesothelioma development. Once genes that moderate tumor development and progression are known, human homologues can be identified and human datasets interrogated to validate their association with disease outcome. Furthermore, our CC-MexTAg animal model enables in-depth study of the tumor microenvironment, allowing the correlation of immune cell infiltration and gene expression signatures with disease development. This strategy provides a detailed picture of the underlying biological pathways associated with mesothelioma susceptibility and progression; knowledge that is crucial for the rational development of new diagnostic and therapeutic strategies. Here we discuss the influence of host genetics on developing an effective immune response to thoracic cancers. We highlight current knowledge gaps, and with a focus on mesothelioma, describe the development and application of the CC-MexTAg to overcome limitations and illustrate how the knowledge gained from this unique study will inform the rational design of future treatments of mesothelioma.

CancerImmunityQTL: a database to systematically evaluate the impact of genetic variants on immune infiltration in human cancer

Tumor-infiltrating immune cells as integral component of the tumor microenvironment are associated with tumor progress, prognosis and responses to immunotherapy. Genetic variants have been demonstrated to impact tumor-infiltrating, underscoring the heritable character of immune landscape. Therefore, identification of immunity quantitative trait loci (immunQTLs), which evaluate the effect of genetic variants on immune cells infiltration, might present a critical step toward fully understanding the contribution of genetic variants in tumor development. Although emerging studies have demonstrated the determinants of germline variants on immune infiltration, no database has yet been developed to systematically analyze immunQTLs across multiple cancer types. Using genotype data from TCGA database and immune cell fractions estimated by CIBERSORT, we developed a computational pipeline to identify immunQTLs in 33 cancer types. A total of 913 immunQTLs across different cancer types were identified. Among them, 5 immunQTLs are associated with patient overall survival. Furthermore, by integrating immunQTLs with GWAS data, we identified 527 immunQTLs overlapping with known GWAS linkage disequilibrium regions. Finally, we constructed a user-friendly database, CancerImmunityQTL (http://www.cancerimmunityqtl-hust.com/) for users to browse, search and download data of interest. This database provides an informative resource to understand the germline determinants of immune infiltration in human cancer and benefit from personalized cancer immunotherapy.