Association study of genetic variation in the autophagy lysosome pathway genes and risk of eight kinds of cancers

Association of ATG10 rs1864183, ATG16L1 rs2241880 and miR-126 with esophageal cancer

BACKGROUND

In India, esophageal cancer (EC) is among the major cause of cancer-related deaths in both sexes. In recent past, autophagy has emerged as one of the crucial process associated with cancer. In the development of EC, the role of autophagy and the precise molecular mechanism involved has yet to be fully understood. Recently, a small number of studies have proposed how variations in autophagy genes affect the growth and development of EC. Micro-RNA's are also known to play a critical role in the development of EC. Here, we examined the relationship between the risk of EC and two single-nucleotide polymorphisms (SNPs) in the key autophagy genes, ATG10 rs1864183 and ATG16L1 rs2241880. We also analyzed the association of miR-107 and miR-126 with EC as these miRNA's are associated with autophagy.

METHODS AND RESULTS

A total of 230 EC patients and 230 healthy controls from North-west Indian population were enrolled. ATG10 rs1864183 and ATG16L1 rs2241880 polymorphism were analyzed using TaqMan genotyping assay. Expression levels of miR-107 and miR-126 were analyzed through quantitative PCR using SYBR green chemistry. We found significant association of CT + CC genotype (OR 0.64, p = 0.022) in recessive model for ATG10 rs1864183 polymorphism with decreased EC risk. For ATG16L1 rs2241880 polymorphism significant association for AG genotype (OR 1.48, p = 0.05) and G allele (OR 1.43, p = 0.025) was observed for increased EC risk. Expression levels of miR-126 were also found to be significantly up regulated (p = 0.008).

CONCLUSION

Our results suggest that ATG10 rs1864183, ATG16L1 rs2241880 and miR-126 may be associated with esophageal carcinogenesis and warrant further investigation.

Establishment of a lysosome-related prognostic signature in breast cancer to predict immune infiltration and therapy response

Background

Lysosomes are instrumental in intracellular degradation and recycling, with their functional alterations holding significance in tumor growth. Nevertheless, the precise role of lysosome-related genes (LRGs) in breast cancer (BC) remains elucidated. This study aimed to establish a prognostic model for BC based on LRGs.

Methods

Employing The Cancer Genome Atlas (TCGA) BC cohort as a training dataset, this study identified differentially expressed lysosome-related genes (DLRGs) through intersecting LRGs with differential expression genes (DEGs) between tumor and normal samples. A prognostic model of BC was subsequently developed using Cox regression analysis and validated within two Gene Expression Omnibus (GEO) external validation sets. Further analyses explored functional pathways, the immune microenvironment, immunotherapeutic responses, and sensitivity to chemotherapeutic drugs in different risk groups. Additionally, the mRNA and protein expression levels of genes within the risk model were examined by utilizing the Gene Expression Profiling Interactive Analysis (GEPIA) and Human Protein Atlas (HPA) databases. Clinical tissue specimens obtained from patients were gathered to validate the expression of the model genes via Real-Time Polymerase Chain Reaction (RT-PCR).

Results

We developed a risk model of BC based on five specific genes (ATP6AP1, SLC7A5, EPDR1, SDC1, and PIGR). The model was validated for overall survival (OS) in two GEO validation sets (p=0.00034 for GSE20685 and p=0.0095 for GSE58812). In addition, the nomogram incorporating clinical factors showed better predictive performance. Compared to the low-risk group, the high-risk group had a higher level of certain immune cell infiltration, including regulatory T cells (Tregs) and type 2 T helper cells (Th2). The high-risk patients appeared to respond less well to general immunotherapy and chemotherapeutic drugs, according to the Tumor Immune Dysfunction and Exclusion (TIDE), Immunophenotype Score (IPS), and drug sensitivity scores. The RT-PCR results validated the expression trends of some prognostic-related genes in agreement with the previous differential expression analysis.

Conclusion

Our innovative lysosome-associated signature can predict the prognosis for BC patients, offering insights for guiding subsequent immunotherapeutic and chemotherapeutic interventions. Furthermore, it has the potential to provide a scientific foundation for identifying prospective therapeutic targets.

LIMP-2 enhances cancer stem-like cell properties by promoting autophagy-induced GSK3β degradation in head and neck squamous cell carcinoma

Cancer stem cell-like cells (CSCs) play an integral role in the heterogeneity, metastasis, and treatment resistance of head and neck squamous cell carcinoma (HNSCC) due to their high tumor initiation capacity and plasticity. Here, we identified a candidate gene named LIMP-2 as a novel therapeutic target regulating HNSCC progression and CSC properties. The high expression of LIMP-2 in HNSCC patients suggested a poor prognosis and potential immunotherapy resistance. Functionally, LIMP-2 can facilitate autolysosome formation to promote autophagic flux. LIMP-2 knockdown inhibits autophagic flux and reduces the tumorigenic ability of HNSCC. Further mechanistic studies suggest that enhanced autophagy helps HNSCC maintain stemness and promotes degradation of GSK3β, which in turn facilitates nuclear translocation of β-catenin and transcription of downstream target genes. In conclusion, this study reveals LIMP-2 as a novel prospective therapeutic target for HNSCC and provides evidence for a link between autophagy, CSC, and immunotherapy resistance.

Transmembrane Protein 175, a Lysosomal Ion Channel Related to Parkinson's Disease

Lysosomes are membrane-bound organelles with an acidic lumen and are traditionally characterized as a recycling center in cells. Lysosomal ion channels are integral membrane proteins that form pores in lysosomal membranes and allow the influx and efflux of essential ions. Transmembrane protein 175 (TMEM175) is a unique lysosomal potassium channel that shares little sequence similarity with other potassium channels. It is found in bacteria, archaea, and animals. The prokaryotic TMEM175 consists of one six-transmembrane domain that adopts a tetrameric architecture, while the mammalian TMEM175 is comprised of two six-transmembrane domains that function as a dimer in lysosomal membranes. Previous studies have demonstrated that the lysosomal K⁺ conductance mediated by TMEM175 is critical for setting membrane potential, maintaining pH stability, and regulating lysosome-autophagosome fusion. AKT and B-cell lymphoma 2 regulate TMEM175's channel activity through direct binding. Two recent studies reported that the human TMEM175 is also a proton-selective channel under normal lysosomal pH (4.5-5.5) as the K⁺ permeation dramatically decreased at low pH while the H⁺ current through TMEM175 greatly increased. Genome-wide association studies and functional studies in mouse models have established that TMEM175 is implicated in the pathogenesis of Parkinson's disease, which sparks more research interests in this lysosomal channel.

Establishment of a novel lysosomal signature for the diagnosis of gastric cancer with in-vitro and in-situ validation

Background

Gastric cancer (GC) represents a malignancy with a multi-factorial combination of genetic, environmental, and microbial factors. Targeting lysosomes presents significant potential in the treatment of numerous diseases, while lysosome-related genetic markers for early GC detection have not yet been established, despite implementing this process by assembling artificial intelligence algorithms would greatly break through its value in translational medicine, particularly for immunotherapy.

Methods

To this end, this study, by utilizing the transcriptomic as well as single cell data and integrating 20 mainstream machine-learning (ML) algorithms. We optimized an AI-based predictor for GC diagnosis. Then, the reliability of the model was initially confirmed by the results of enrichment analyses currently in use. And the immunological implications of the genes comprising the predictor was explored and response of GC patients were evaluated to immunotherapy and chemotherapy. Further, we performed systematic laboratory work to evaluate the build-up of the central genes, both at the expression stage and at the functional aspect, by which we could also demonstrate the reliability of the model to guide cancer immunotherapy.

Results

Eight lysosomal-related genes were selected for predictive model construction based on the inclusion of RMSE as a reference standard and RF algorithm for ranking, namely ADRB2, KCNE2, MYO7A, IFI30, LAMP3, TPP1, HPS4, and NEU4. Taking into account accuracy, precision, recall, and F1 measurements, a preliminary determination of our study was carried out by means of applying the extra tree and random forest algorithms, incorporating the ROC-AUC value as a consideration, the Extra Tree model seems to be the optimal option with the AUC value of 0.92. The superiority of diagnostic signature is also reflected in the analysis of immune features.

Conclusion

In summary, this study is the first to integrate around 20 mainstream ML algorithms to construct an AI-based diagnostic predictor for gastric cancer based on lysosomal-related genes. This model will facilitate the accurate prediction of early gastric cancer incidence and the subsequent risk assessment or precise individualized immunotherapy, thus improving the survival prognosis of GC patients.

Autophagy mediates ER stress and inflammation in Helicobacter pylori-related gastric cancer

Autophagy is a cellular degradation mechanism, which is triggered by the bacterium Helicobacter pylori. A single nucleotide polymorphism (SNP) in the autophagy gene ATG16L1 (rs2241880, G-allele) has been shown to dysregulate autophagy and increase intestinal endoplasmic reticulum (ER) stress. Here, we investigate the role of this SNP in H.pylori-mediated gastric carcinogenesis and its molecular pathways. ATG16L1 rs2241880 was genotyped in subjects from different ethnic cohorts (Dutch and Australian) presenting with gastric (pre)malignant lesions of various severity. Expression of GRP78 (a marker for ER stress) was assessed in gastric tissues. The effect of ATG16L1 rs2241880 on H.pylori-mediated ER stress and pro-inflammatory cytokine induction was investigated in organoids and CRISPR/Cas9 modified cell lines. Development of gastric cancer was associated with the ATG16L1 rs2241880 G-allele. Intestinal metaplastic cells in gastric tissue of patients showed increased levels of ER-stress. In vitro models showed that H.pylori increases autophagy while reducing ER stress, which appeared partly mediated by the ATG16L1 rs2241880 genotype. H.pylori-induced IL-8 production was increased while TNF-α production was decreased, in cells homozygous for the G-allele. The ATG16L1 rs2241880 G-allele is associated with progression of gastric premalignant lesions and cancer. Modulation of H.pylori-induced ER stress pathways and pro-inflammatory mediators by ATG16L1 rs2441880 may underlie this increased risk.

Association of genetic variants in autophagy-lysosome pathway genes with susceptibility and survival to prostate cancer

BACKGROUND

Previous studies have indicated the connections between autophagy-lysosome pathway genes dysfunction and prostate cancer, but few studies have investigated whether single nucleotide polymorphisms (SNPs) in autophagy-lysosome pathway genes are implicated in prostate cancer risk and survival.

MATERIALS AND METHODS

Logistic regression analysis and stepwise Cox regression analysis were conducted in 4,662 cases and 3,114 controls from the Prostate, Lung, Colorectal and Ovarian (PLCO) Cancer Screening Trial. The false positive rate probability (FPRP) method was applied to correct for multiple comparisons. Gene-based analysis was calculated by versatile gene-based association study approach.

RESULTS

We found that SLC11A1 rs7573065 significantly increased the risk of prostate cancer [adjusted odds ratio (OR) = 1.24, 95% confidence interval (CI) = 1.06-1.46, P = 7.02 × 10^-3, FPRP = 0.082]. Furthermore, rs7573065 was confirmed as the independent predicator of overall survival (OS) for prostate cancer patients [Hazard ratio (HR) = 1.30, 95% CI = 1.01-1.66, P = 0.041]. The significant association between SLC11A1 and prostate cancer risk was calculated by gene-based analysis (P = 0.030). We also observed that the mRNA of SLC11A1 in prostate tumor tissues was significantly over-expressed than that in normal tissues.

CONCLUSION

This study suggested that rs7573065 in SLC11A1 was associated with an increased risk and poor OS of prostate cancer. Our findings may provide evidence for genetic variants in autophagy-lysosome pathway as the risk and prognostic biomarkers for prostate cancer.

Prognostic Significance of Cytoplasmic SPNS2 Expression in Patients with Oral Squamous Cell Carcinoma

Background and Objectives: Oral squamous cell carcinoma (OSCC) is a malignant disease with a particularly high incidence in Taiwan. Our objective in this study was to elucidate the involvement of sphingolipid transporter 2 (SPNS2) expression and SPNS2 protein expression in the clinicopathological indexes and the clinical outcomes of OSCC patients. Materials and Methods: Immunohistochemistry analysis was performed for SPNS2 protein expression in samples from 264 cases of OSCC. Correlations of SPNS2 expression with clinicopathological variables and patient survival were analyzed. Results: Our results revealed that the cytoplasmic protein expression of SPNS2 in OSCC tissue specimens was lower than in normal tissue specimens. Negative cytoplasmic protein expression of SPNS2 was significantly correlated with T status and stage. Kaplan-Meier survival curve analysis revealed that negative cytoplasmic SPNS2 expression was predictive of poorer overall survival of OSCC patients in stage III/IV. We also determined that low SPNS2 expression was an independent prognostic factor related to overall survival among OSCC patients in stage III/IV from univariate Cox proportional hazard models. Multivariate Cox proportional hazard models revealed that cytoplasmic SPNS2 expression, T status, lymph node metastasis, and histological grade were independent prognostic factors for survival. Conclusions: Overall, this study determined that SPNS2 protein may be a useful prognostic marker for OSCC patients and potential therapeutic target for OSCC treatment.

Integration of GWAS and eQTL Analysis to Identify Risk Loci and Susceptibility Genes for Gastric Cancer

Genome-wide association studies (GWAS) have identified several susceptibility loci for gastric cancer (GC), but the majority of identified single-nucleotide polymorphisms (SNPs) fall within the non-coding region and are likely to exert their biological function by modulating gene expression. To systematically estimate expression-associated SNPs (eSNPs) that confer genetic predisposition to GC, we evaluated the associations of 314,203 stomach tissue-specific eSNPs with GC risk in three GWAS datasets (2,631 cases and 4,373 controls). Subsequently, we conducted a gene-based analysis to calculate the cumulative effect of eSNPs through sequence kernel association combined test and Sherlock integrative analysis. At the SNP-level, we identified two novel variants (rs836545 at 7p22.1 and rs1892252 at 6p22.2) associated with GC risk. The risk allele carriers of rs836545-T and rs1892252-G exhibited higher expression levels of DAGLB (P = 3.70 × 10^–18) and BTN3A2 (P = 3.20 × 10^–5), respectively. Gene-based analyses identified DAGLB and FBXO43 as novel susceptibility genes for GC. DAGLB and FBXO43 were significantly overexpressed in GC tissues than in their adjacent tissues (P = 5.59 × 10^–7 and P = 3.90 × 10^–6, respectively), and high expression level of these two genes was associated with an unfavorable prognosis of GC patients (P = 1.30 × 10^–7 and P = 7.60 × 10^–3, respectively). Co-expression genes with these two novel genes in normal stomach tissues were significantly enriched in several cancer-related pathways, including P53, MAPK and TGF-beta pathways. In summary, our findings confirm the importance of eSNPs in dissecting the genetic basis of GC, and the identified eSNPs and relevant genes will provide new insight into the genetic and biological basis for the mechanism of GC development.

Prognostic significance of Spinster homolog gene family in acute myeloid leukemia

Acute myeloid leukemia (AML) is a clonal and heterogeneous disease characterized by proliferation of immature myeloid cells, with impaired differentiation and maturation. Spinster homolog (SPNS) is a widely distributed transmembrane transporter, which assists sphingolipids in playing their roles through the cell membrane. However, the expression and clinical implication of the SPNS family has not been investigated in AML. From the Cancer Genome Atlas database, a total of 155 AML patients with complete clinical characteristics and SPNS1-3 expression data were contained in our study. In patients who received chemotherapy only, high expressions of SPNS2 and SPNS3 had adverse effects on event-free survival (EFS) and overall survival (OS) (all P<0.05). However, in the allogeneic hematopoietic stem cell transplantation (allo-HSCT) group, we only found a significant difference in OS between the high and low SPNS3 expression groups (P=0.001), while other SPNS members showed no effect on survival. Multivariate analysis indicated that high SPNS2 expression was an independent risk factor for both EFS and OS in chemotherapy patients. The results confirmed that high expression of SPNS2 and SPNS3 were poor prognostic factors, and the effect of SPNS2 can be neutralized by allo-HSCT.

Functional variants of autophagy-related genes are associated with the development of hepatocellular carcinoma

AIMS

Hepatocellular carcinoma (HCC) is the most common primary liver cancer, and accounts for substantial morbidity and mortality. Autophagy plays an essential role in the development and progression of HCC. This study aims to evaluate whether genetic variants in autophagy-related genes (ATGs) affect the development of HCC.

MATERIALS AND METHODS

We conducted a case-control study with 986 HCC cases and 1000 healthy controls to analyze 14 functional variants of five ATGs (ATG3, ATG5, ATG10, ATG12 and ATG16L1) among a Chinese population.

KEY FINDINGS

We found ATG5 rs17067724 (G vs A: OR = 0.80; 95% CI = 0.65-0.98; P = 0.031), ATG10 rs1864183 (G vs A: OR = 1.29; 95% CI = 1.07-1.57; P = 0.009), ATG10 rs10514231 (C vs T: OR = 1.41; 95% CI = 1.15-1.73; P = 0.001), ATG12 rs26537 (C vs T: OR = 1.16; 95% CI = 1.02-1.33; P = 0.030), and ATG16L1 rs4663402 (T vs A: OR = 1.28; 95% CI = 1.01-1.63; P = 0.044) were significantly associated with HCC risk. Specifically, ATG10 rs10514231 kept significant association even adjusted for Bonferroni correction (P = 0.001 × 14 = 0.014). Bioinformatics analyses showed that allele C of ATG10 rs10514231 was significantly correlated with higher expression of ATG10 gene in both HCC tissues and normal liver tissues. Dual-luciferase reporter assay presented that cell lines transfected with vectors containing the risk allele C of rs10514231 showed higher relative luciferase activity compared to that containing the allele T.

SIGNIFICANCE

These results suggested that ATG10 rs10514231 might contribute to an allele-specific effect on the expression of host gene ATG10 and explain a fraction of HCC genetic susceptibility. Our study would benefit the construction of early warning model, early prevention, screening, even therapeutic target of HCC.