The Effects of Differentially-Expressed Homeobox Family Genes on the Prognosis and HOXC6 on Immune Microenvironment Orchestration in Colorectal Cancer

Construction of Prognostic Prediction Models for Colorectal Cancer Based on Ferroptosis-Related Genes: A Multi-Dataset and Multi-Model Analysis

Background

Colorectal cancer (CRC) remains a significant health burden globally, necessitating a deeper understanding of its molecular landscape and prognostic markers. This study characterized ferroptosis-related genes (FRGs) to construct models for predicting overall survival (OS) across various CRC datasets.

Methods

In TCGA-COAD dataset, differentially expressed genes (DEGs) were identified between tumor and normal tissues using DESeq2 package. Prognostic genes were identified associated with OS, disease-specific survival, and progression-free interval using survival package. Additionally, FRGs were downloaded from FerrDb website, categorized into unclassified, marker, and driver genes. Finally, multiple models (Coxboost, Elastic Net, Gradient Boosting Machine, LASSO Regression, Partial Least Squares Regression for Cox Regression, Ridge Regression, Random Survival Forest [RSF], stepwise Cox Regression, Supervised Principal Components analysis, and Support Vector Machines) were employed to predict OS across multiple datasets (TCGA-COAD, GSE103479, GSE106584, GSE17536, GSE17537, GSE29621, GSE39084, GSE39582, and GSE72970) using intersection genes across DEGs, OS, disease-specific survival, and progression-free interval, and FRG categories.

Results

Six intersection genes (ASNS, TIMP1, H19, CDKN2A, HOTAIR, and ASMTL-AS1) were identified, upregulated in tumor tissues, and associated with poor survival outcomes. In the TCGA-COAD dataset, the RSF model demonstrated the highest concordance index. Kaplan-Meier analysis revealed significantly lower OS probabilities in high-risk groups identified by the RSF model. The RSF model exhibited high accuracy with AUC values of 0.978, 0.985, and 0.965 for 1-, 3-, and 5-year survival predictions, respectively. Calibration curves demonstrated excellent agreement between predicted and observed survival probabilities. Decision curve analysis confirmed the clinical utility of the RSF model. Additionally, the model's performances were validated in GSE29621 dataset.

Conclusions

The study underscores the prognostic relevance of 6 intersection genes in CRC, providing insights into potential therapeutic targets and biomarkers for patient stratification. The RSF model demonstrates robust predictive performance, suggesting its utility in clinical risk assessment and personalized treatment strategies.

Prognostic and therapeutic potential of gene profiles related to tertiary lymphoid structures in colorectal cancer

The role of tertiary lymphoid structures (TLS) in oncology is gaining interest, particularly in colorectal carcinoma, yet a thorough analysis remains elusive. This study pioneered a novel TLS quantification system for prognostic and therapeutic response prediction in colorectal carcinoma, alongside a comprehensive depiction of the TLS landscape. Utilizing single-cell sequencing, we established a TLS score within the Tumor Immune Microenvironment (TIME). Analysis of tertiary lymphoid structure-related genes (TLSRGs) in 1,184 patients with colon adenocarcinoma/rectum adenocarcinoma (COADREAD) from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases led to the identification of two distinct molecular subtypes. Differentially expressed genes (DEGs) further segregated these patients into gene subtypes. A TLS score was formulated using gene set variation analysis (GSVA) and its efficacy in predicting immunotherapy outcomes was validated in two independent cohorts. High-scoring patients exhibited a 'hot' immune phenotype, correlating with enhanced immunotherapy efficacy. Key genes in our model, including C5AR1, APOE, CYR1P1, and SPP1, were implicated in COADREAD cell proliferation, invasion, and PD-L1 expression. These insights offer a novel approach to colorectal carcinoma treatment, emphasizing TLS targeting as a potential anti-tumor strategy.

Role of homeobox genes in cancer: immune system interactions, long non-coding RNAs, and tumor progression

The intricate interplay between Homeobox genes, long non-coding RNAs (lncRNAs), and the development of malignancies represents a rapidly expanding area of research. Specific discernible lncRNAs have been discovered to adeptly regulate HOX gene expression in the context of cancer, providing fresh insights into the molecular mechanisms that govern cancer development and progression. An in-depth comprehension of these intricate associations may pave the way for innovative therapeutic strategies in cancer treatment. The HOX gene family is garnering increasing attention due to its involvement in immune system regulation, interaction with long non-coding RNAs, and tumor progression. Although initially recognized for its crucial role in embryonic development, this comprehensive exploration of the world of HOX genes contributes to our understanding of their diverse functions, potentially leading to immunology, developmental biology, and cancer research discoveries. Thus, the primary objective of this review is to delve into these aspects of HOX gene biology in greater detail, shedding light on their complex functions and potential therapeutic applications.

Effects of HOX family regulator-mediated modification patterns and immunity characteristics on tumor-associated cell type in endometrial cancer

Endometrial cancer (UCEC) is one of three major malignant tumors in women. The HOX gene regulates tumor development. However, the potential roles of HOX in the expression mechanism of multiple cell types and in the development and progression of tumor microenvironment (TME) cell infiltration in UCEC remain unknown. In this study, we utilized both the The Cancer Genome Atlas (TCGA) database and International Cancer Genome Consortium (ICGC) database to analyze transcriptome data of 529 patients with UCEC based on 39 HOX genes, combing clinical information, we discovered HOX gene were a pivotal factor in the development and progression of UCEC and in the formation of TME diversity and complexity. Here, a new scoring system was developed to quantify individual HOX patterns in UCEC. Our study found that patients in the low HOX score group had abundant anti-tumor immune cell infiltration, good tumor differentiation, and better prognoses. In contrast, a high HOX score was associated with blockade of immune checkpoints, which enhances the response to immunotherapy. The Real-Time quantitative PCR (RT-qPCR) and Immunohistochemistry (IHC) exhibited a higher expression of the HOX gene in the tumor patients. We revealed that the significant upregulation of the HOX gene in the epithelial cells can activate signaling pathway associated with tumour invasion and metastasis through single-cell RNA sequencing (scRNA-seq), such as nucleotide metabolic proce and so on. Finally, a risk prognostic model established by the positive relationship between HOX scores and cancer-associated fibroblasts (CAFs) can predict the prognosis of individual patients by scRNA-seq and transcriptome data sets. In sum, HOX gene may serve as a potential biomarker for the diagnosis and prediction of UCEC and to develop more effective therapeutic strategies.

Downregulation of DIP2B as a prognostic marker inhibited cancer proliferation and migration and was associated with immune infiltration in lung adenocarcinoma via CCND1 and MMP2

Background

DIP2B is related to cancer progression. This study investigated the roles and pathways of DIP2B in lung adenocarcinoma (LUAD).

Methods

DIP2B expression and the relationship between survival time of cancer patients and DIP2B expression were analyzed. The relationship between DIP2B expression and survival time in LUAD patients was evaluated by a meta-analysis. Cox and survival analyses were used to evaluate the prognostic factors and construct a prognostic nomogram. The mechanisms and effects of DIP2B and the relationship between DIP2B expression and the immune microenvironment were investigated using bioinformatics, CCK-8, western blotting, and transwell experiments.

Results

DIP2B was overexpressed in LUAD tissues. DIP2B overexpression was associated with shorter prognosis and was an unfavorable risk factor for prognosis in LUAD patients. DIP2B co-expressed genes were involved in cell division, DNA repair, cell cycle, and others. Inhibition of DIP2B expression could downregulate the proliferation, migration, and invasion of LUAD A549 and H1299 cells, which was related to the decrease in CCND1 and MMP2 protein expression. BRCA1 overexpression was associated with short prognosis, and the nomogram formed by DIP2B and BRCA1 was associated with a poor prognosis in LUAD patients. DIP2B expression correlated with immune cells (such as CD8 T cells, Tcm, and iDCs) and cell markers.

Conclusion

DIP2B is a potential biomarker of poor prognosis and the immune microenvironment in LUAD. Inhibition of DIP2B expression downregulated cancer cell proliferation, migration, and invasion, which might be related to the decrease in CCND1 and MMP2 protein expression. DIP2B-related nomograms might be useful tools for predicting the prognosis of LUAD patients.

Prognostic and therapeutic insights into colorectal carcinoma through immunogenic cell death gene profiling

While the significance of immunogenic cell death (ICD) in oncology is acknowledged, its specific impact on colorectal carcinoma remains underexplored. In this study, we delved into the role of ICD in colorectal carcinoma, a topic not yet comprehensively explored. A novel ICD quantification system was developed to forecast patient outcomes and the effectiveness of immunotherapy. Utilizing single-cell sequencing, we constructed an ICD score within the tumor immune microenvironment (TIME) and examined immunogenic cell death related genes (ICDRGs). Using data from TCGA and GEO, we discovered two separate molecular subcategories within 1,184 patients diagnosed with colon adenocarcinoma/rectum adenocarcinoma (COADREAD). The ICD score was established by principal component analysis (PCA), which classified patients into groups with low and high ICD scores. Further validation in three independent cohorts confirmed the model's accuracy in predicting immunotherapy success. Patients with higher ICD scores exhibited a "hot" immune phenotype and showed increased responsiveness to immunotherapy. Key genes in the model, such as AKAP12, CALB2, CYR61, and MEIS2, were found to enhance COADREAD cell proliferation, invasion, and PD-L1 expression. These insights offered a new avenue for anti-tumor strategies by targeting ICD, marking advances in colorectal carcinoma treatment.

Unveiling HOXB7 as a novel diagnostic and prognostic biomarker through pan-cancer computer screening

We attempted to investigate the role of HOXB7 in tumor progression and evolution by means of an extensive computer screening analysis of various cancer types. We performed univariate Cox regression and Kaplan-Meier survival analyses to assess the impact of HOXB7 on overall survival (OS), disease-specific survival (DSS), and progression-free interval (PFI) in different types of cancer. Furthermore, we examined the relationship between HOXB7 and several clinical features: tumor microenvironment, immune regulatory genes, immune checkpoints, tumor mutational burden (TMB), and microsatellite instability (MSI). We performed gene set enrichment analysis to gain deeper insights into the potential molecular mechanisms of HOXB7, and validated our findings through functional assays in cells, including methyl thiazolyl tetrazolium cytotoxicity and Transwell invasion assays. HOXB7 expression was associated with different clinical characteristics in numerous malignancies. Higher HOXB7 expression was associated with worse OS, DSS, and PFI in some cancer types. In particular, HOXB7 expression was favorably associated with immune cell infiltration, immune regulatory genes, immunological checkpoints, TMB, and MSI in malignancies. Furthermore, we identified a strong link between copper death-associated gene expression and HOXB7 expression. According to the findings of this study, HOXB7 might serve as an appealing focus for tumor diagnosis and immunotherapy and a prospective indicator of prognosis.

HOXC6: A promising biomarker linked to an immunoevasive microenvironment in colorectal cancer based on TCGA analysis and cohort validation

HOXC6 plays an essential part of the carcinogenesis of solid tumors, but its functional relevance within the immune contexture in patients with colorectal cancer (CRC) is still uncertain. We intended to investigate the predictive value of HOXC6 expression for survival outcomes and its correlation with immune contexture in CRC patients by utilizing the Cancer Genome Atlas database (n = 619). Validation was performed in cohorts from Zhongshan Hospital (n = 200) and Shanghai Cancer Center (n = 300). Immunohistochemical (IHC) staining was utilized to compare the levels of immunocytes infiltrating the tumor between the groups with high and low expression of HOXC6. Elevated levels of HOXC6 expression in CRC tissues were linked to malignant progression and poor prognosis. HOXC6 as a risk factor for survival of CRC patients was confirmed. Receiver operating characteristic analysis confirmed its diagnostic value, and a reliable prognostic nomogram was constructed. KEGG analysis and GSEA showed that HOXC6 participated in immune regulation, and its expression was tightly linked to the abundance of infiltrating immunocytes. HOXC6 was upregulated in patients diagnosed with CRC within the two cohorts, and high HOXC6 levels were correlated with a worse prognosis. The high-HOXC6 expression group showed increased infiltration of Treg cells, CD68+ macrophages, CD66b+ neutrophils, and CD8+ T-cells and elevated levels of PD-L1 and PD-1, but decreased levels of granzyme B and perforin. These findings suggest that HOXC6 abundance in patients with CRC determines a poor prognosis, promotes an immunoevasive environment, and directs CD8+ T-cell dysfunction. HOXC6 is expected to become a prospective biomarker for the outcome of CRC.

Comprehensive analysis of HOX family genes in endometrial cancer

Background

Endometrial cancer (EC) is one of the most prevalent malignancies in the female population. Homeoboxes (HOXs) are a large family of transcription factors that have a variety of functions in biological processes (BPs), including developmental differentiation, and their dysregulated expression has been implicated in tumorigenesis. However, the involvement of HOXs in EC has received little attention. Thus, we aimed to identify the potential role of HOXs in EC from a multi-omics perspective through bioinformatics analysis.

Methods

We obtained transcriptome, mutation, and methylation data and the corresponding clinical data for normal and tumor tissues from The Cancer Genome Atlas (TCGA) database. Abnormal expression of HOXs in EC was identified via differential analysis, and the diagnostic value of HOXs in EC was assessed with the receiver operating characteristic (ROC) method. Univariate and multivariate Cox regression models were employed to evaluate the predictive efficacy of HOXs in EC. Methylation and mutation analyses revealed epigenetic and genetic sequence alterations in HOXs. Single-sample gene set enrichment analysis (ssGSEA) was used to explore the altered immune microenvironment in EC. Moreover, the gene activity and pathway enrichment of downstream key HOX genes were revealed through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis in EC.

Results

HOXs were found to be linked to the growth of EC and potentially playing a role in establishing the tumor immune microenvironment in patients with EC. HOXB9 was found to be a vital prognostic molecule in patients with EC and is expected to contribute to a novel treatment approach.

Conclusions

We used bioinformatics techniques to clarify the potential role of HOXs from a multi-omics perspective, and our findings provide a foundation for future investigations into the molecular mechanisms of HOXs in EC.

CFHTF2 Is Needed for Vegetative Growth, Conidial Morphogenesis and the Osmotic Stress Response in the Tea Plant Anthracnose (Colletotrichum fructicola)

Tea is an important cash crop worldwide, and its nutritional value has led to its high economic benefits. Tea anthracnose is a common disease of tea plants that seriously affects food safety and yield and has a far-reaching impact on the sustainable development of the tea industry. In this study, phenotypic analysis and pathogenicity analysis were performed on knockout and complement strains of HTF2-the transcriptional regulator of tea anthracnose homeobox-and the pathogenic mechanism of these strains was explored via RNA-seq. The MoHox1 gene sequence of the rice blast fungus was indexed, and the anthracnose genome was searched for CfHTF2. Evolutionary analysis recently reported the affinity of HTF2 for C. fructicola and C. higginsianum. The loss of CfHTF2 slowed the vegetative growth and spore-producing capacity of C. fructicola and weakened its resistance and pathogenesis to adverse conditions. The transcriptome sequencing of wild-type N425 and CfHTF2 deletion mutants was performed, and a total of 3144 differentially expressed genes (DEGs) were obtained, 1594 of which were upregulated and 1550 of which were downregulated. GO and KEGG enrichment analyses of DEGs mainly focused on signaling pathways such as the biosynthesis of secondary metabolites. In conclusion, this study lays a foundation for further study of the pathogenic mechanism of tea anthracnose and provides a molecular basis for the analysis of the pathogenic molecular mechanism of CfHTF2.

Engrailed 2 triggers the activation of multiple phosphorylation-induced signaling pathways in both transcription-dependent and -independent manners

Homeodomain (HD)-containing proteins are typically recognized as transcription factors. Engrailed 2 (EN2) is an HD-containing protein that is highly expressed in various types of cancers, however, the mechanism underlying the biological function of EN2 is not fully understood. Here, we report a transcription-independent function of EN2 in addition to its role as a transcription factor. EN2 expression leads to the activation of multiple signaling pathways mediated by phosphorylation cascades. A phosphoproteomic analysis revealed that the phosphorylation status of numerous protein sites was altered after EN2 is expressed. Notably, EN2 was shown to interact with a myriad of proteins implicated in phosphorylation signaling cascades, as determined by immunoprecipitation-mass spectrometry (IP-MS). We validated the interaction between EN2 and B55α, the regulatory subunit of the PP2A-B55α complex, and confirmed that the phosphatase activity of the complex was suppressed by EN2 binding. To target EN2-induced malignancy, two kinds of small molecules were utilized to inhibit the EN2-activated NF-κB and AKT signaling pathways. A clear synergistic effect was observed when the activation of the two pathways was simultaneously blocked. Collectively, the data show that EN2 functions in a transcription-independent manner in addition to its role as a transcription factor. This finding may have therapeutic implications in treating esophageal squamous cell carcinoma (ESCC).

High expression of HOXC6 predicts a poor prognosis and induces proliferation and inflammation in multiple myeloma cells

BACKGROUND

Multiple myeloma (MM) is a common blood system malignance accompanied by monoclonal plasma cell hyperplasia. Homeobox C6 (HOXC6) acts as an oncogene in various cancers, but its function on MM is elusive.

OBJECTIVE

The role of HOXC6 on MM development was clarified in this study.

METHODS

HOXC6 expression and its clinical significance were determined in the peripheral blood samples collected from forty MM patients and thirty healthy adult volunteers. The overall survival was evaluated by Kaplan-Meier analysis with the log-rank test. Cell viability, proliferation and apoptosis were measured by CCK-8, EdU assay and Flow cytometry in U266 and MM.1R cells. Tumor growth was estimated by a xenograft assay. The apoptosis of tumor tissues was evaluated using TUNEL staining. The protein level in tissues was tested by immunohistochemistry.

RESULTS

The HOXC6 expression was elevated in MM and high HOXC6 level was associated with the poor overall survival of MM. Besides, the HOXC6 expression was associated with hemoglobin level and ISS stage. Furthermore, silencing HOXC6 suppressed cell proliferation, induced cell apoptosis, and restrained the secretion of inflammatory factors (TNF-α, IL-6, and IL-8) in MM cells through inactivating the NF-κB pathway. Moreover, silencing HOXC6 suppressed the tumor growth of MM, the inflammatory factors levels, and the activation of NF-κB pathway but enhanced apoptosis in vivo.

CONCLUSION

HOXC6 was elevated in MM and associated with poor survival. Knockdown of HOXC6 suppressed proliferation, inflammation and tumorigenicity of MM cells via inactivating the NF-κB pathway. HOXC6 may be a meaningful target for MM therapy.

Identification of cuproptosis-based molecular subtypes, construction of prognostic signature and characterization of immune landscape in colon cancer

Background

Cuproptosis is a newly discovered form of cell death induced by targeting lipoacylated proteins involved in the tricarboxylic acid cycle. However, the roles of cuproptosis-related genes (CRGs) in the clinical outcomes and immune landscape of colon cancer remain unknown.

Methods

We performed bioinformatics analysis of the expression data of 13 CRGs identified from a previous study and clinical information of patients with colon cancer obtained from The Cancer Genome Atlas and Gene Expression Omnibus databases. Colon cancer cases were divided into two CRG clusters and prognosis-related differentially expressed genes. Patient data were separated into three corresponding distinct gene clusters, and the relationships between the risk score, patient prognosis, and immune landscape were analyzed. The identified molecular subtypes correlated with patient survival, immune cells, and immune functions. A prognostic signature based on five genes was identified, and the patients were divided into high- and low-risk groups based on the calculated risk score. A nomogram model for predicting patient survival was developed based on the risk score and other clinical features.

Results

The high-risk group showed a worse prognosis, and the risk score was related to immune cell abundance, microsatellite instability, cancer stem cell index, checkpoint expression, immune escape, and response to chemotherapeutic drugs and immunotherapy. Findings related to the risk score were validated in the imvigor210 cohort of patients with metastatic urothelial cancer treated with anti-programmed cell death ligand 1.

Conclusion

We demonstrated the potential of cuproptosis-based molecular subtypes and prognostic signatures for predicting patient survival and the tumor microenvironment in colon cancer. Our findings may improve the understanding of the role of cuproptosis in colon cancer and lead to the development of more effective treatment strategies.

Identification of necroptosis-related subtypes, development of a novel signature, and characterization of immune infiltration in colorectal cancer

Introduction

Necroptosis, a type of programmed cell death, has recently been extensively studied as an important pathway regulating tumor development, metastasis, and immunity. However, the expression patterns of necroptosis-related genes (NRGs) in colorectal cancer (CRC) and their potential roles in the tumor microenvironment (TME) have not been elucidated.

Methods

We explored the expression patterns of NRGs in 1247 colorectal cancer samples from genetics and transcriptional perspective. Based on a consensus clustering algorithm, we identified NRG molecular subtypes and gene subtypes, respectively. Furthermore, we constructed a necroptosis-related signature for predicting overall survival time and verified the predictive ability of the model. Using the ESTIMATE, CIBERSORT, and ssGSEA algorithms, we assessed the association between the above subtypes, scores and immune infiltration.

Results

Most NRGs were differentially expressed between CRC tissues and normal tissues. We found that distinct subtypes exhibited different NRGs expression, patients' prognosis, immune checkpoint gene expression, and immune infiltration characteristics. The scores calculated from the necroptosis-related signature can be used to classify patients into high-risk and low-risk groups, with the high-risk group corresponding to reduced immune cell infiltration and immune function, and a greater risk of immune dysfunction and immune escape.

Discussion

Our comprehensive analysis of NRGs in CRC demonstrated their potential role in clinicopathological features, prognosis, and immune infiltration in the TME. These findings help us deepen our understanding of NRGs and the tumor microenvironment landscape, and lay a foundation for effectively assessing patient outcomes and promoting more effective immunotherapy.

Ferroptosis regulator FANCD2 is associated with immune infiltration and predicts worse prognosis in lung adenocarcinoma

Lung adenocarcinoma (LUAD) remains one of the leading causes of cancer-related death. Although immunotherapy has been shown to improve survival in LUAD patients, only a select group of LUAD patients could benefit from it. The correlation between ferroptosis and the tumor immune environment requires further investigation in the setting of LUAD. An analysis using The Cancer Genome Atlas (TCGA)-LUAD cohort systematically evaluated the expression levels of ferroptosis regulators between LUAD and normal tissues and demonstrated the correlation of ferroptosis regulators with the immune checkpoint B7-H3 expression. Based on consensus clustering analysis, we divided LUAD patients into two subtypes according to the expression pattern of ferroptosis regulators. Cluster 2 patients showed more favorable overall survival (OS) (p < 0.001) and disease-free survival (DFS) (p < 0.001) than Cluster 1 patients. CIBERSORT analysis indicated that Cluster 1 patients harbored higher infiltrated levels of uncharacterized cells, CD4⁺ T cells (nonregulatory), and myeloid dendritic cells, while Cluster 2 patients were more correlated with B cells, M1 macrophages, natural killer cells (NK cells) and regulatory T cells (Tregs). More importantly, we identified FANCD2 as a potentially unfavorable prognostic factor that was overexpressed in LUAD and positively associated with the checkpoint molecule B7-H3 expression. In addition, higher FANCD2 expression was related to a higher tumor immune dysfunction and exclusion (TIDE) score, indicating lower responder rates to cancer immunotherapeutics. In summary, our study suggested a relationship between immune infiltration and ferroptosis and that FANCD2 is a potential biomarker for clinical outcomes and a therapeutic target for LUAD therapy concerning ferroptotic regulation. Our findings may help to advance personalized treatment and improve the prognosis of LUAD.

Cancer Risk Score Prediction Based on a Single-Nucleotide Polymorphism Network

Objectives

Genome-wide association studies (GWAS) are performed to study the associations between genetic variants with respect to certain phenotypic traits such as cancer. However, the method that is commonly used in GWAS assumes that certain traits are solely affected by a single mutation. We propose a network analysis method, in which we generate association networks of single-nucleotide polymorphisms (SNPs) that can differentiate case and control groups. We hypothesize that certain phenotypic traits are attributable to mutations in groups of associated SNPs.

Methods

We propose a method based on a network analysis framework to study SNP-SNP interactions related to cancer incidence. We employed logistic regression to measure the significance of all SNP pairs from GWAS for the incidence of colorectal cancer and computed a cancer risk score based on the generated SNP networks.

Results

We demonstrated our method in a dataset from a case-control study of colorectal cancer in the South Sulawesi population. From the GWAS results, 20,094 pairs of 200 SNPs were created. We obtained one cluster containing four pairs of five SNPs that passed the filtering threshold based on their p-values. A locus on chromosome 12 (12:54410007) was found to be strongly connected to the four variants on chromosome 1. A polygenic risk score was computed from the five SNPs, and a significant difference in colorectal cancer risk was obtained between the case and control groups.

Conclusions

Our results demonstrate the applicability of our method to understand SNP-SNP interactions and compute risk scores for various types of cancer.