- Research
- Open access
- Published:
Genetic assessment of pathogenic germline alterations in lysosomal genes among Asian patients with pancreatic ductal adenocarcinoma
Journal of Translational Medicine volume 21, Article number: 730 (2023)
Abstract
Background
Lysosomes are closely linked to autophagic activity, which plays a vital role in pancreatic ductal adenocarcinoma (PDAC) biology. The survival of PDAC patients is still poor, and the identification of novel genetic factors for prognosis and treatment is highly required to prevent PDAC-related deaths. This study investigated the germline variants related to lysosomal dysfunction in patients with PDAC and to analyze whether they contribute to the development of PDAC.
Methods
The germline putative pathogenic variants (PPV) in genes involved in lysosomal storage disease (LSD) was compared between patients with PDAC (n = 418) and healthy controls (n = 845) using targeted panel and whole-exome sequencing. Furthermore, pancreatic organoids from wild-type and KrasG12D mice were used to evaluate the effect of lysosomal dysfunction on PDAC development. RNA sequencing (RNA-seq) analysis was performed with established PDAC patient-derived organoids (PDOs) according to the PPV status.
Results
The PPV in LSD-related genes was higher in patients with PDAC than in healthy controls (8.13 vs. 4.26%, Log2 OR = 1.65, P = 3.08 × 10–3). The PPV carriers of LSD-related genes with PDAC were significantly younger than the non-carriers (mean age 61.5 vs. 65.3 years, P = 0.031). We further studied a variant of the lysosomal enzyme, galactosylceramidase (GALC), which was the most frequently detected LSD variant in our cohort. Autophagolysosomal activity was hampered when GALC was downregulated, which was accompanied by paradoxically elevated autophagic flux. Furthermore, the number of proliferating Ki-67+ cells increased significantly in pancreatic organoids derived from Galc knockout KrasG12D mice. Moreover, GALC PPV carriers tended to show drug resistance in both PDAC cell line and PDAC PDO, and RNA-seq analysis revealed that various metabolism and gene repair pathways were upregulated in PDAC PDOs harboring a GALC variant.
Conclusions
Genetically defined lysosomal dysfunction is frequently observed in patients with young-onset PDAC. This might contribute to PDAC development by altering metabolism and impairing autophagolysosomal activity, which could be potentially implicated in therapeutic applications for PDAC.
Background
Despite medical advances, pancreatic cancer remains one of most fatal malignancies worldwide. According to Global Cancer Statistics 2020, pancreatic cancer, the seventh leading cause of cancer deaths worldwide, accounts for nearly as many deaths (466,000) as there are cases (496,000) because of its poor prognosis [1]. In most cases, delayed diagnosis is responsible for this dismal outcome: only 10–20% patients with pancreatic cancer are diagnosed with resectable disease [2]. Accordingly, early identification of high-risk populations and improvement in survival percentage are necessary. Approximately 10% of all pancreatic cancers are attributable to inherited risk factors [3, 4]. The genetic basis of familial or hereditary pancreatic cancer can be explained in 21% families based on previously described hereditary cancer-related genes and in 35% families based on low-frequency variants in other DNA repair genes [5]. More specifically, germline mutations in genes related to DNA instability, such as CDKN2A, TP53, MLH1, BRCA2, ATM, and BRCA1, are well known to be associated with the development of pancreatic adenocarcinoma (PDAC) [6].
Deleterious germline mutations are also evident in roughly 3.9–7% of pancreatic cancer cases lacking definitive family history [7, 8]. Consequently, knowledge regarding germline variant-driven susceptibility will be essential for understanding PDAC, as it will help in defining a high-risk population in cases of both familial and sporadic pancreatic cancer and will act as a rational background for novel anti-cancer drug development. Extensive research is warranted to elucidate the genetic mechanisms underlying PDAC. In addition to the aforementioned cancer-predisposing genes (CPGs), lysosomal storage diseases (LSDs) comprise more than 50 disorders caused by mutations in genes involved in the functioning of endosome–lysosome proteins [9]. Lysosomes are the main digestive compartments within cells and are closely linked to autophagy, a primary intracellular degradation system that derives its degradative abilities from lysosomes [10]. In cancer cells, lysosomes affect growth factor signaling via the endocytic degradation of growth factors, their receptors, or signal transduction mediators to modulate signaling output [11]. In addition, defective autophagy has been suggested to contribute to carcinogenesis, possibly owing to the reduced removal of defective organelles or damaged cells [12]. Accordingly, dysfunction and cancer development may be closely related. Indeed, a previous study has shown a possible association between rare variants of LSD-related genes and cancer [13]; compared with the average population, pathogenic variants of LSD genes were significantly enriched in the cancer cohort. PDAC is strongly associated with germline mutations in several LSD genes, including SGSH, MAN2B1, and IDUA. These results are consistent with evidence from a mouse model suggesting that autophagy suppresses cancer initiation [14]. Recently, adult-onset chronic diseases have been shown to originate from the heterozygote background of LSD gene variants. A good example is the association between Parkinson's disease and mutant GBA (encoding beta-glucocerebrosidase) heterozygote carriers [15].
We hypothesized that some heterozygote carriers of LSD genes may develop PDAC owing to the suppression of lysosomal dysfunction. In addition, once established, cancer cells use autophagy to promote survival during nutrient stress and recycle cell components to support a transformed phenotype [16], which is highly dependent on enhanced lysosomal function to facilitate the degradation, clearance, and recycling of cellular material delivered by increased rates of vesicle trafficking via autophagy and micropinocytosis [17,18,19]. Hence, impaired autophagic activity may contribute to cancer initiation; however, its role in established cancer cells may differ, as autophagy plays a biphasic role in cancer initiation and progression [20].
Taken together, the present study aimed to investigate the oncogenic effect of LSD heterozygote carrier status on PDAC. In particular, we aimed to (1) evaluate the clinical significance of rare LSD gene variants on large-scale in patient with PDAC and in healthy control cohorts, (2) focus on an ethnically homogenous population such as Koreans, since rare variant analysis is primarily affected by ethnicity, (3) evaluate the functional consequences of LSD gene dysfunction using mouse pancreatic organoids and human PDAC cell lines, and (4) assess characteristics of PDAC harboring rare LSD variants to understand the features of LSD-related PDAC using patient-derived organoids (PDOs).
Methods
Study cohort
In total, 418 patients diagnosed with PDAC between November 2011 and August 2020 at the Samsung Medical Center (SMC; n = 222) and Seoul National University Hospital (SNUH; n = 196) were prospectively enrolled and followed up until the end of 2021. Clinical and laboratory data were collected from the electronic health records (EHR). A cancer-free normal control (CFNC) cohort of 845 healthy volunteers was prospectively constructed at the SNUH Healthcare Checkup Center. The age distribution of the PDAC patients ranged from 35 to 87 years old. The age of all healthy volunteers was > 50 years, without a history of cancer, as proven by EHR, thereby excluding young age onset PDAC patients. This study was conducted in accordance with the principles of the Declaration of Helsinki, and the study protocol was approved by the SMC Institutional Review Board (IRB) and SNUH IRB (Seoul, South Korea). Written informed consent was obtained from all patients and volunteers, and all specimens were collected according to the IRB regulations and approval (IRB No. 2018-12-065, 1705-031-852).
DNA sequencing analysis
Genomic DNA collected from peripheral blood was extracted using a QIAamp DNA tissue kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions. Targeted panel or whole exome sequencing (WES) was used to evaluate the status of LSD germline variants in the study cohorts. First, targeted panel sequencing data were generated to identify LSD germline variants in 493 CFNC individuals. Then, the data generation method was changed for additional samples to identify variants in all genes. DNA from 352 additional CFNC individuals and 418 PDAC patients were used to generate WES data using an Illumina NovaSeq 6000 sequencer (Illumina, San Diego, CA, USA). Library capturing was performed using a probe from Integrated DNA Technologies. For the LSD panel sequencing, the DNA of 493 subjects from the CFNC cohort was sequenced using an Illumina HiSeq 2500 platform. Libraries were constructed using the ACCEL-NGS 2S DNA library kit (Swift Biosciences, Ann Arbor, MI, USA), which included 42 LSD genes (Additional file 4: Table S1).
Our in-house variant calling pipeline followed the Genome Analysis Toolkit (GATK) best practices recommended by the Broad institute [21, 22]. The sequence reads were aligned to the human reference genome (hg19) using the Burrows-Wheeler Aligner-MEM v0.7.10 algorithm [23]. Conversion to Binary Alignment Map (BAM) was performed using Picard v1.130. and deposited in the Sequence Read Archive (SRA bioproject accession # PRJNA929903). For insertions and deletions (indel) realignment, duplicated fragment elimination and base quality score recalibration were performed using the Genome Analysis Tool Kit (GATK v3.8.1, Broad Institute) [22]. Using the CollectHsMetrics of Picard tool, single nucleotide variants (SNVs) and indels were detected using the HaplotypeCaller of GATK v3.8.1 which is preferred for assessing germline variant calling by using de novo assembly of reads. Variants in each sample's genomic variant calling format (VCF) were merged, and a joint calling approach was performed using GATK to empower the variant discovery. Finally, sequencing errors were filtered out by assessing variant quality score recalibration (VQSR) and GATK's statistical modeling approach for variant filtration. Low-quality variants were discarded if the total coverage was less than 10, and sequencing reads with suspected bias were considered if the variant allele frequencies (VAF) were less than 20%. Additionally, variants listed in the ENCODE/DUKE [24] and DAC blacklist [25] regions were discarded, and only variants in the ENCODE/CRG GEM mappability region (75 mers) [26] were extracted to filter out well-known variants from low mappability regions. To generate a consistent set of variant filtration and functional annotations, the compiled variants in the VCF file were annotated using ANNOVAR [27] to perform filter-based functional annotation and a Variant Effect Predictor (VEP) [28] was used to explicate the gene-based information of canonical transcripts. In addition, to ensure that only rare germline variants were identified, we performed variant filtration using allele frequency (AF) information from The Genome Aggregation Database (gnomAD) [29] (gnomad.exomes.r2.1.1). The protein truncation variants and clinically validated variants in the ClinVar database was extracted as putative pathogenic variants (PPVs): (1) Tier1 variants were defined as protein truncating variants (PTVs), which included splice donor and acceptor site, frameshift indel, stop gain, and lost variants, as well as non-benign or likely benign loci annotated with Clinvar [30]; (2) Genetic variants with well-known clinical risks (pathogenic, likely pathogenic, association, and risk factor) and related phenotypes are clearly defined in ClinVar as Tier 2.
Most LSDs are caused by variants in genes that regulate lysosomal enzymes, which tend to be inherited in an autosomal recessive pattern. Somatic alterations were checked following the ‘two-hit’ hypothesis, in which a pathogenic variant occurs in one allele and the ‘second hit’ occurs in the other allele, leading to increased cancer risk. Therefore, somatic copy number alterations (CNA) across 42 LSD genes were determined using the CNVkit (https://github.com/etal/cnvkit) by comparing the tumor organoid DNA BAM files to germline samples (matched normal pancreatic blood). The mean read depths for each target (interval) were computed and normalized against a single reference of pooled standard samples, and the B-allele frequency was calculated. The observed log2 copy number ratios of the region < − 0.4 were derived as copy number loss, and CAN segments were visualized using the copy number package.
RNA sequencing analysis
For RNA sequencing (RNA-seq), total RNA was isolated using a RNeasy mini kit (Qiagen). TruSeq stranded mRNA (Illumina) was used to prepare the RNA-seq libraries. The 150 bp paired-end sequencing of these libraries was performed using a NovaSeq 6000 sequencing system (Illumina). The quality of the cDNA libraries was evaluated using an Agilent 2100 bioanalyzer (Agilent, Santa Clara, CA, USA). Quality was confirmed using the NGS QC Toolkit (version 2.3.3). RNA-seq reads of pancreatic tumor organoids were aligned to the human reference genome (hg19) using spliced transcript alignment to a reference (STAR) version 2.5.3 a [31] expression count was estimated using RNA-seq by Expectation Maximization (RSEM-1.3.0) [32] and normalized using the EdgeR TMM method before performing differential expression analysis. To investigate how variants in the LSD-related genes affected gene expression, we analyzed the differentially expressed genes (DEGs) in LSD carrier and non-carrier groups using the transcript expression level from the RNA-seq data. Differential expression between LSD gene variant carriers and non-carriers was estimated using gene-specific read counts and the R package platform DESeq2 [33]. The log2 scale variance stabilizing transformation was completed using DEGs. Gene set enrichment analyses (GSEA) which is a functional enrichment analysis to identify association with group of genes was performed along with LSD carriers and non-carriers using Java GSEA application version 3.0. Functionally meaningful pathways were explored using Kyoto Encyclopedia of Genes and Genomes (KEGG) [34] and Gene Ontology (GO) databases of The Molecular Signatures Database (MSigDB). [35].
Cell lines
Rare pathogenic germline variants in 42 LSD genes were screened in PDAC cell lines from the Cancer Cell Line Encyclopedia database [36]. The PK59 cell line with a GALC mutation (rs138577661 and rs137854543) was maintained in Roswell Park Memorial Institute (RPMI)-1640 medium with 10% fetal bovine serum (FBS). HPAFII, PANC1, AsPC1, and Capan-1 cells without GALC variants were maintained in Eagle’s minimum essential medium supplemented with 10% FBS, Dulbecco′s modified Eagle′s medium (DMEM) supplemented with 10% FBS, RPMI-1640 medium supplemented with 10% FBS, and Iscove's modified Dulbecco's medium supplemented with 20% FBS, respectively. All media and FBS were purchased from Gibco (Billings, MT, USA).
Analysis of GALC enzyme activity
Endogenous GALC activity was detected using a lysosomal GALC analysis kit (Marker Gene Technologies, Inc., Eugene, OR, USA), according to the manufacturer’s instructions. Briefly, 4 × 106 cells were harvested using the provided buffer and lysed using 2 × 30 s sonication cycles in a Bioruptor™ Pico. First, the lipidic fluorogenic substrate was incubated with 50 µg protein for 2 h at 37 °C. The fluorescence signals at excitation/emission = 365/454 nm were detected using a fluorometer (SYNERGY/HTS; BioTek Instruments, Winooski, VT, USA). The fluorescence of PK59 cells was normalized to the relative GALC activity.
Analysis of autophagic flux
To detect the LC3B signal, 3–4 × 105 cells were seeded in 6-well plates. The next day, the cells were washed twice with pre-warmed phosphate buffered saline (PBS) before treatment with 100 μM chloroquine (CQ; Sigma-Aldrich, St. Louis, MO, USA) for 1 h. After protein extraction, LC3B (Sigma-Aldrich) and β-actin (Santa Cruz) were detected using immunoblotting. The chemiluminescent signals were detected and visualized using a LAS-3000 luminescent image analyzer (FUJIFILM, Tokyo, Japan). The autophagy flux unit (A.F.U. = [LC3B-II/LC3B-I]CO(+)/[LC3B-II/LC3B-I]CO(-)) was calculated by analyzing the band intensities of LC3B-I and LC3B-II from three independent experiments using the ImageJ software.
PK59 and HPAF-II cells stably expressing GFP-LC3 were transfected with retrovirus (pBABEpuroGFP-LC3; Addgene, Watertown, MA, USA). After treatment with CQ (100 μM) for 5 h, GFP signals from LC3 were analyzed using a fluorescence microscope (Inverted Microscope Eclipse Ti-S; Nikon, Tokyo, Japan). Five images were selected to manually count the number of total and GFP-LC3 positive cells.
Mouse pancreatic organoid culture
The protocols for mouse experiments were approved by the Institutional Animal Care and Use of Committee (IACUC) of Seoul National University (SNU-150724-4). KrasG12D/+ (Jackson Laboratory, Bar Harbor, Maine, USA) mouse was crossed with Pdx1-Cre mouse (Jackson Laboratory) to obtain KrasG12D/+; Pdx1-Cre mouse (KrasG12D). Since an activating point mutation of the KRAS oncogene in codon 12 (exon 2) occurs in the majority of PDAC cases (70–95%), the transgenic KrasG12D mouse model is used in all pathophysiological studies [37]. Pancreatic ducts were isolated from wild-type and KrasG12D mice and lysed using collagenase P (Roche, Basel, Switzerland) and DNase I (Worthington Biochemical Corp., Lakewood, NJ, USA). The lysed pancreatic ductal cells were seeded in matrigel matrix (Corning, Corning, NY, USA) and grown in culture medium containing Advanced DMEM/F12, B27™ Supplement (Gibco), GlutaMAX™ (Gibco), R-spondin1, mEGF (Peprotech, Cranbury, NY, USA), mNoggin (Peprotech), hFGF10 (Peprotech), N-acetylcysteine (Sigma-Aldrich), A83-01 (Tocris, Bristol, UK), and nicotinamide (Sigma-Aldrich). Lentiviral CRISPR/Cas9 was used to generate Galc knockout (KO) mouse pancreatic organoids. The sgRNA targeting Galc was cloned into the lentiCRISPRv2 vector. LentiCRISPRv2-sgGalc, pMD2G, and psPAX2 were transfected into 293FT. The lentiviruses were collected from the cells and concentrated using a Lenti-X™ concentrator (Takara, Shiga, Japan). Lentiviruses were transduced into dissociated single cells using TrypLE (Gibco). The transduced cells were recovered in a matrigel matrix and selected using 2 µg/μL puromycin. The generation of Galc KO mouse pancreatic organoids was confirmed using semi-quantitative reverse transcription polymerase chain reaction (RT-PCR) with the following primers: Galc_F:5ʹ-AGG TCT CCA GCG AGT GAG AAT CAT AG-3,ʹ Galc_R:5ʹ-TGT GTG AGC TGA TAC CCA GAT AGG AG-3.ʹ To analyze ubiquitinated proteins in mouse pancreatic organoids, the cells were treated with 10 μM MG132 (Enzo, USA) for 24 h. For immunoblotting, the concentration of the extracted organoid protein was measured using the bicinchoninic acid (BCA)™ assay (Pierce, Celbio, Milan, Italy) and bovine serum albumin was used as the standard. Equal amounts of the protein extracts were electrophoresed on 10% sodium dodecyl sulfate–polyacrylamide gels and electro-blotted onto polyvinylidene difluoride membranes (Millipore SPA, Milan, Italy). The membranes were then incubated for 1 h at room temperature with blocking solution consisting of 5% skim milk in Tris-buffered saline solution and Tween 20 [TBST; 100 mmol/L Tris (pH 7.5), 0.9% NaCl, and 0.1% Tween 20] and probed overnight at 4 °C using anti-mTOR (Cell Signaling), anti-phospho-mTOR(Cell signaling), anti-rictor (Cell Signaling), anti-LC3B (Cell Signaling), anti-ubiquitin (Cell Signaling), anti-SQTM1/p62 (Abcam), and anti-beta actin (Abcam) antibodies (1:1000 in blocking solution). Horseradish peroxidase-conjugated IgG (1:5000 in blocking solution; Cell Signaling) was used to detect specific proteins. Immunodetection was performed using chemiluminescent substrates (Pierce).
Immunostaining and in situ hybridization
For the immunofluorescence assay, organoids were isolated from the matrigel using cell recovery solution (Corning) and washed with cold PBS. The isolated organoids were fixed in 4% paraformaldehyde for 1 h, permeabilized in 1% PBS-T (Triton X-100) for 1 h and blocked for 1 h at room temperature in a blocking solution (3% bovine serum albumin in 0.2% PBS-T). The samples were incubated overnight with primary and secondary antibodies (Additional file 5: Table S2) at 4 °C. Then, they were mounted using the VECTASHIELD antifade mounting solution with 4′,6-diamidino-2-phenylindole (DAPI; Vector Laboratories, Burlingame, CA, USA) and imaged using a Zeiss LSM 700 confocal microscope (Zeiss, Oberkochen, Germany). The images were processed using ImageJ software.
For immunohistochemistry (IHC) and in situ hybridization, the fixed organoids were embedded in paraffin blocks. IHC staining was performed using the OptiView DAB IHC detection kit (Ventana, Oro Valley, AZ, USA), according to the manufacturer's instructions. The Ki67 antibody (Invitrogen, Waltham, MA, USA) was used for IHC, and a mouse Galc probe (NM_008079.4) was used for in situ hybridization. All the cells were analyzed using the QuPath image analyzer.
Measuring autophagosome and autolysosome
Organoids were isolated from the matrigel and dissociated into single cells. FUW-mCherry-GFP-LC3 plasmids (Addgene plasmid # 110060) were transduced in the cells using lentiviruses. The transduced cells were then embedded and recovered in matrigel for 48 h. The organoids were fixed in 4% paraformaldehyde for 1 h and mounted using VECTASHIELD Antifade mounting solution with DAPI. Imaging was performed using the DeltaVision system (Applied Precision/GE Healthcare, Issa-quah, WA, USA) and processed using the ImageJ software.
Patient-derived organoid culture
Human PDAC specimens were homogenized using a GentleMACS™ tissue dissociator and a human tumor dissociation kit (Miltenyi Biotec, Bergisch Gladbach, Germany) according to the manufacturer's instructions. After filtering using a 70 μm strainer, the suspended cells were plated on matrigel (Corning) and grown in complete medium: advanced DMEM/F12 supplemented with GlutaMAX™, containing 10 mM HEPES, antibiotic antimycotics, B27 supplement, N2 supplement (all from Thermo Fisher Scientific, Waltham, MA, USA), 1 mM n-acetylcysteine, 10 mM nicotinamide (all from Sigma-Aldrich), 60 ng/mL murine Wnt-3a, 500 ng/mL human R-spondin 1, 10 nM human gastrin, 50 ng/mL human Noggin, 50 ng/mL human epidermal growth factor (EGF), 100 ng/mL, human fibroblast growth factor 10 (FGF10), 0.5 μM A83-01 (all from Peprotech), 1 × Primocin (InvivoGen, San Diego, CA, USA), and 10 μM Y-27632. The organoids were cultured in a 37 °C humidified incubator, and the culture medium was partially changed twice a week. Resuspended PDAC cells were seeded into 384-well plates (500 cells/well) with technical duplicates and treated with therapeutic drugs including gemcitabine, nab-paclitaxel, irinotecan and Olaparib (all from Selleckchem) for seven days. Cell viability was accessed using an adenosine triphosphate monitoring system based on firefly luciferase (ATPlite 1step; PerkinElmer) and estimated using the EnVision multilabel reader (PerkinElmer). The relative cell viability for each dose was obtained after normalization with dimethyl sulfoxide (DMSO) per plate.
Statistical analysis
A linear regression model was used to assess the association between the PPV of each gene and phenotypic characteristics, with a significance cut-off of P < 0.05.
where n = case (1) or control (0), germline variants = number of samples that carry rare pathogenic germline variants, gender = male (0) and female (1), and patient age at diagnosis were used as inputs for the regression model. Chi-squared and Fisher's exact tests were used to evaluate the associations between categorical variables. They were applied for the most prevalent comparisons for testing the association of variant incidence according to the disease phenotype. Independent t-test and one-way analysis of variance were used to assess the association between continuous variables. Survival analysis was performed using the Kaplan–Meier method. All statistical analyses were performed using either SPSS (version 25.0; IBM SPSS Statistics, Armonk, NY: IBM Corp) or R (http://www.r-project.org/).
Results
Baseline characteristics of the patients
The baseline characteristics of the 418 patients with PDAC are presented in Table 1. Patients’ ages ranged from 35 to 87 years, with a median age of 65 years. Among these, 243 (58.1%) were male, and 175 (41.9%) were female. The median value of body mass index was 22.7 kg/m2, and 146 patients (34.9%) had diabetes mellitus at the time of diagnosis. The median values of carcinoembryonic antigen and carbohydrate antigen 19-9 were 2.4 ng/mL and 150.6 U/mL, respectively. According to the 8th edition of the AJCC on Cancer staging system, 78 (18.7%), 123 (29.4%), 112 (26.8%), and 105 (25.1%) patients had stage I, II, III, and IV disease, respectively. Univariate analysis revealed a significant association between the stage and overall survival (OS), with median OS for stages I, II, III, and IV being 37.0, 24.0, 18.7, and 10.6 months (P < 0.001), respectively, implying that our cohort was a general pancreatic cancer cohort rather than a biased cohort with an unusual disease course. In addition, we recruited 845 individuals without any history of cancer through a health examination center of SNUH over 3 years who were grouped as the control set (named CFNC), respectively. These samples were selected only from those aged 50 years or older (age ranged from 51 to 91 years).
Enrichment of LSD germline variants in the PDAC cohort
Thirty-seven PPVs, including SNVs and indels, were detected in 21 of the 42 LSD genes (ARSA, GAA, GALC, GLB1, GNPTAB, GNPTG, GUSB, HEXA, HEXB, HGSNAT, HYAL1, IDUA, MAN2B1, MCOLN1, NAGLU, NPC1, PSAP, SGSH, SMPD1, SUMF1, and TPP1; Additional file 4: Table S1). Most variants were PTVs (24 of 37) and 48.7% were clinically proven pathogenic variants of ClinVar [30] (pathogenic or likely pathogenic; 18 of 37). Specific information on individual variants found in PDAC are listed in Additional file 6: Table S3. Among those, two typical variants in GALC were identified most frequently, although not statistically significant; rs138577661 was found in 2.4% PDAC patients (10 out of 418 PDAC patients), which was twice as many as that in the CFNC group (1.2%, 10 out of 845 individuals; P = 1.05 × 10–1 in Chi-squared test), and rs200607029 was found in 1.0% PDAC patients (4 out of 418 PDAC patients) compared to 0.4% in the CFNC group (3 out 845 CFNC group; P = 1.75 × 10–1 in Chi-squared test).
Overall, 34 Korean patients with PDAC (34/418, 8.13%) harbored at least one genomic variant, whereas 36 germline carriers (36/845, 4.26%) were identified among CFNC (P = 6.92 × 10–3 by Chi-squared test). The most frequently mutated gene in GALC (3.6%) was the most frequently mutated gene in Korean patients with PDAC, followed by HEXB (1.2%), GAA (0.5%), and NAGLU (0.5%; Additional file 4: Table S1). We also estimated the effect of PPV using a regression model after adjusting for gender and age using samples from WES data (418 patients with PDAC and 352 controls to avoid variant selection bias due to technical differences). Results revealed PDAC enrichment of LSD variants with Log2 odds ratio (OR) of 1.65 (P = 3.08 × 10–3). Of the 418 Korean patients with PDAC, 15 (3.6%) had CPG variants, including BRCA1/2, ATM, and COL7A1, and two patients harbored rare pathogenic germline variants in both LSD genes and CPGs (Table 1).
Early onset of PDAC in patients with LSD gene carriers
We tested whether pathogenic germline variant carriers in LSD gene developed PDAC at a younger age than non-carriers, similar to that observed in the well-known CPG carriers who are more likely to develop cancer at a younger age (Additional file 1: Figure S1 and Table 2). Two patients with both CPG and LSD gene variants were classified as CPG carriers, and the characteristics of each study group were compared. The mean ages at diagnosis of PDAC in CPG carriers and LSD carriers were 60.4 and 61.7 years, respectively, which were significantly lower than that in patients without CPG or LSD variants (65.3 years; P = 0.025). In addition, for all LSD gene carriers (n = 34), including two who also have CPG mutations, the mean age at diagnosis of PDAC was 61.5 years, which was significantly lower than 65.3 years for non-carriers of either LSD or CPG (n = 371) (P = 0.031). This indicates that LSD carriers could be one of the high-risk populations who showed early development of PDAC. Other than age, clinical characteristics, including tumor stage, did not differ significantly according to CPG or LSD gene carrier status (Table 2).
Defective lysosomal function in Kras G12D /Galc knockout mouse-derived pancreatic organoids
To investigate the effects of autophagy regulation and lysosomal function on PDAC development, a lysosomal enzyme GALC was eliminated from KrasG12D mice-derived pancreatic organoids using CRISPR/CRISPR-associated protein 9 (Cas9). As mentioned, Galc variants were the most frequent LSD PPV in our cohort. Galc expression was reduced in KrasG12D/Galc knockout (KO) organoid, which was confirmed using semi-quantitative RT-PCR (Additional file 2: Fig. S2a), RNA sequencing (as the median expression value of KrasG12D and KrasG12D/Galc KO organoid: 5.13 vs. 1.04, P = 0.05 in Wilcoxon rank sum test (Additional file 2: Fig. S2b), and in situ hybridization (Additional file 2: Fig. S2c, d).
The autophagic flux in KrasG12D organoids after Galc depletion was observed via western blotting of LC3B-II, a standard autophagosome marker (LC3B-I is lipidated upon autophagy activation to form LC3B-II), and SQSTM1/p62, a receptor involved in autophagy [38, 39]. LC3B-II level increased slightly in KrasG12D organoids and markedly in KrasG12D/Galc KO organoids (Fig. 1a). In contrast, SQSTM1/p62 level did not differ significantly between KrasG12D organoids and KrasG12D/Galc KO organoids. Then, we assessed the levels of mammalian target of rapamycin (mTOR), as mTORC2 (mTOR-rictor component) negatively regulates autophagy [40]. The results showed that the level of phospho-mTOR (p-mTOR) decreased considerably in the KrasG12D/Galc KO organoids, while that of mTOR was intact. Moreover, rictor (rapamycin-insensitive companion of mTOR), which is a component of mTORC2, was downregulated in KrasG12D/Galc KO organoids, suggesting that inactivation of mTORC2 led the autophagic flux (Fig. 1a).
In addition, lysosomal function was examined via immunostaining of lysosome-associated membrane protein 1 (LAMP1). The results showed that the intensity of LAMP1 decreased considerably in KrasG12D/Galc KO organoids (Fig. 1b), as represented in the bar graph (Fig. 1c), indicating that Galc depletion in KrasG12D causes lysosomal dysfunction.
As autophagic activity increased with defective lysosome function, we investigated the formation of autophagosomes and autolysosomes in wild-type, KrasG12D, and KrasG12D/Galc KO organoids. We transfected the mCherry-GFP-LC3 plasmid in three different types of organoids to visualize the formation of authophagosomes and autolysosomes, and their numbers were analyzed. Owing to the acidity of GFP, autopahgosomes are marked by both mCherry and GFP (yellow) signals, while the autolysosomes were marked with only mCherry (red) [41]. The number of autophagosomes was higher in KrasG12D/Galc KO organoids than in the wild-type or KrasG12D organoids, whereas the number of autolysosomes did not differ among organoids of the three different genotypes (Fig. 1e, f). Furthermore, we found that the accumulation of ubiquitinated proteins was enhanced in KrasG12D/Galc KO organoids compared to that in the KrasG12D organoids both before and after treatment with MG132, a proteasome inhibitor, for 24 h (Fig. 1d).
These results suggested that defective GALC in KrasG12D organoids enhanced autophagic flux and promoted autophagosome formation. However, autolysosome formation and degradation did not follow because of lower lysosomes and lysosomal dysfunction, causing accumulation of ubiquitinated proteins in KrasG12D/Galc KO organoids.
Increased proliferation of Kras G12D /Galc knockout mouse-derived pancreatic organoids
To understand the regulatory mechanism underlying lysosomal dysfunction that affect pancreatic tumorigenesis, we analyzed the gene expression patterns of KrasG12D and KrasG12D/Galc KO pancreatic organoids using RNA-seq. In DEG analysis, the majority (81.5%) of the genes, including MIRG, PIGR, SYT14, PLEKHD1, and SLC2A1, were found to be downregulated in KrasG12D/Galc KO organoids. In addition, low expression of monoamine oxidase A (MAOA), a key factor that regulates apoptosis/autophagy by targeting the repressor element-1 silencing transcription factor (REST1), was identified in KrasG12D/Galc KO organoids. In contrast, the expression of genes, including STAG3, GATA5, and CACNB2, which are associated with cell proliferation or tumorigenesis, and PRAP1, which is primarily involved in the inhibition of tumor cell apoptosis, increased in mouse KrasG12D/Galc KO pancreatic organoids (Fig. 2a).
To assess cell proliferation, we examined Ki-67 expression in mouse pancreatic organoids. The number of Ki-67-positive proliferating cells was considerably higher in KrasG12D/Galc KO organoids and lower in wild-type and KrasG12D organoids (Fig. 2b, c, Additional file 3: Fig. S3). Furthermore, GSEA revealed that KrasG12D/Galc KO organoids were strongly associated with RAS signaling, cytokine production, and cell death in GO terms by MSigDB. However, regulation of translation, RNA processing, protein ubiquitination, and glycolysis were significantly downregulated in the KrasG12D/Galc KO pancreatic organoids (Fig. 2d). In summary, transcriptomic analysis of Galc KO mouse-derived pancreatic organoid indicated that GALC dysfunction resulted in impaired autophagic function and markedly enhanced cell proliferation.
Enhanced autophagic flux in PDAC cell lines with LSD gene variants
Next, we examined the mechanisms underlying autophagy regulation in PDAC cells with LSD PPVs. The endogenous activity of lysosomal GALC enzyme was lower in PK59 (Galc variant) cells than in HPAFII and Capan-1 (Galc wild-type) cells (Fig. 3a). Furthermore, LC3B-II was observed after treatment with CQ, an inhibitor of autophagy, to monitor autophagy flux. The increase in the amount of LC3B-II after CQ treatment was prominent in PK59 cells (Fig. 3b and c). Moreover, we generated PK59 and HPAFII cells stably expressing GFPto a-LC3 and monitored the number of LC3-positive cells after CQ treatment. The relative increase in GFP-expressing cells was more pronounced in PK59 cells than that in HPAFII cells (Figs. 3d and e). These results indicated enhancement of autophagic flux in PDAC cells with low GALC activity, which is associated with Galc variants.
Drug responses in PDAC cells with LSD germline variants
We performed a drug test using PDAC PDOs and PDAC cell lines to examine their responses to treatment according to LSD PPV carriers (Fig. 4). PDAC PDOs and PDAC cell lines with or without GALC variants were treated with representative therapeutic drugs for PDAC, including gemcitabine, nab-paclitaxel, irinotecan, and Olaparib for 7 days, and the viability was assessed for calculating the half maximal inhibitory concentration (IC50). We observed that Galc carriers ORG2383 and PK59 showed higher IC50 values than other PDAC PDOs (Fig. 4a and b) and PDAC cell lines (Fig. 4c and d). Thus, Galc carriers tended to be more resistant to drugs used for treating PDAC.
Evaluation of pathways associated with LSD germline variants in human PDAC PDOs
To investigate the functional consequences of lysosomal dysfunction in established PDAC, we conducted RNA-seq of LSD PPV carriers (n = 2, both GALC heterozygote carriers) and non-carriers (n = 27) using our PDAC PDOs. Differential expression analysis identified 240 DEGs, of which 68 were upregulated, including S100P, AFAP1-AS1, TFCP2L1, and RASA3, and 172 were downregulated, including APOM, APOC1, ASB4, TF, and PLG in LSD carriers [Log2FoldChange > 1, false discovery rate (FDR) < 0.05] (Fig. 5a) compared to that in LSD non-carriers. However, the small sample size might be a limitation of this study, as it might affect the reproducibility of the results of DEG analysis. Interestingly, GSEA analysis to test the carcinogenic processes depending on LSD gene variation showed that amino acid biosynthesis-related pathways, including primary bile acid biosynthesis, steroid hormone biosynthesis, linoleic acid metabolism, alanine, aspartate, and glutamate metabolism, and glycerolipid metabolism were strongly upregulated in LSD carriers in KEGG (Fig. 5b). Furthermore, we identified upregulation of the DNA repair system (nucleotide excision repair, mismatch repair, and homologous recombination) in LSD carriers. In addition, metabolic upregulation was confirmed using the GO analysis from a biological perspective that occurs within cells in LSD carriers (Fig. 5c). However, most cell-to-cell signaling and organ development pathways were downregulated in PDAC PDOs from LSD carriers (Fig. 5c).
Discussion
Germline mutations in genes are associated with risk of cancer development and may cause early onset cancer [42]. Most genes involved in carcinogenesis have functions related to cell division and proliferation. The present study investigated the oncogenic effects of LSD heterozygosity in patients with PDAC. We evaluated the clinical implications of the rare variants of LSD genes on a large scale PDAC patients and healthy controls, and functional analyses were performed using mouse pancreatic organoids, PDAC cell lines, and PDAC PDOs. First, we found that LSD PPV carriers were enriched in the PDAC population compared to that in the healthy controls and they developed PDAC at a younger age than the non-carriers. Second, we showed paradoxically increased autophagic flux due to impaired autophagolysosome activity, which may contribute to PDAC development, using mouse pancreatic organoids. Third, PDAC PDOs with LSD PPV showed upregulation of metabolic pathways, which supports the idea that genes involved in lysosome function may also contribute to PDAC development via alterations in autophagic activity.
The enrichment of germline LSD PPV was approximately three times (Log2OR 1.6–1.8) in the Korean PDAC population compared to that in the healthy controls. This enrichment is also accompanied by a younger age of onset of PDAC by 3.6 years compared to that of non-LSD PPV carriers, similar to that observed with other CPGs. A previous study has reported the enrichment of six representative germline PPVs, including BRCA1/2 and TP53, in the PDAC population compared to that in the general population, with ORs ranging from 2.58 to 12.3 (all ORs were more significant than five except for 2.58 for BRCA1) [6] implying that the oncogenic potential of germline LSD PPV was weaker than that of conventional CPGs. Dysfunction of lysosomes and autophagy, which unlocks phenotypic plasticity, is an emerging hallmark of cancer [43]. Thus, LSD PPV may be used as a biomarker for distinguishing the high-risk population for PDAC after validation.
Technically, due to the genetic diversity of rare disease variants [44, 45], ethnic factors should be addressed in this type of analysis. In fact, several studies have revealed distinct germline predisposition according to ethnicity in many diseases [46,47,48,49]. We have eliminated this issue by focusing on a single ethnic Korean population for the germline PPV analysis. As Koreans have relatively homogenous genetic background, comparison of PPV frequencies across cohorts is feasible [50]. In the PPV analysis, GALC, HEXB, GAA, and NAGLU were identified most frequently in the Korean PDAC population. Among these, the contributions of GALC PPVs (rs138577661 and rs200607029), which are well-known variants of the Krabbe disease in the East Asian populations, was found to be considerable in our study [51]. The functionality of the rs138577661 variant was confirmed by performing enzymatic testing of a PDAC cell line (PK59) harboring the exact variant. Overall, the variant-level data also suggested the importance of ethnicity in this study on rare PPVs. Hence, extrapolation of our data to other ethnic populations should be performed with caution, and consideration of whether other rare PPVs are implicated in another ethnic population might be required.
Among the various cancers, we selected PDAC for investigation based on the results of a previous study that revealed the possible link between LSD PPV and PDAC [13] and its well-known relation with lysosomes, autophagy [52], and abnormal metabolism [53]. From the perspective of tumor initiation and development, dysfunctional lysosomes impede the autophagolysosomal activity of abnormal cells, resulting in cancer cell survival and growth [54]. From the perspective of cell death, lysosomal or autophagic cell death is not adequate for elimination of cancerous cells [55]. Indeed, evidence suggests that dysfunctional lysosomes contribute to cancer initiation [16, 56]. Our results show conclusively that Galc knockdown promotes cell proliferation (elevated Ki-67) in cooperation with Kras mutations in mouse pancreatic organoids. In addition, we observed that autophagolysosome function decreased, which subsequently increased autophagic flux as a feedback mechanism, suggesting alteration of autophagic activity with GALC downregulation. Considering the cellular function of lysosomes, we hypothesized that cooperation between oncogenic mutations such as RAS and lysosomal dysfunction is necessary for carcinogenesis. However, as we did not observe all the features of cancer in the mouse organoids, future studies using mice with knockdown of lysosome-related genes are necessary to confirm spontaneous cancer development. Furthermore, manipulation of specific variants such as rs138577661 is worthwhile to understand the carcinogenic mechanism of specific variants related to lysosome dysfunction.
However, once a tumor is established, autophagy plays a vital role in maintaining the cancer in many ways. In particular, for PDAC, nutrient-poor and hypoxic conditions allow cancer cells to thrive in harsh environments, which render autophagy important for recycling nutrients [53]. In addition, in PDAC, macropinocytosis plays a role similar to that of autophagy. Both autophagy and macropinocytosis rely on lysosomes for the final degradation of products. Macropinocytosis is induced by the blockade of autophagic activity in PDAC via NRF2 induction [57, 58], which is indicative of their compensatory roles. We interpreted that the increased autophagic flux in PDAC cell lines originates from lysosomal dysfunction, which compensates for the hampered autophagolysosomes. This increased autophagic flux may favor the survival of established cancer cell, with paradoxically increased autophagy and/or macropinocytosis in cancer cells. Furthermore, upregulation of autophagy not only enhances the tumor survival, but also induces drug resistance due to renewal of by cytoplasmic materials, gene repair, alterations in drug concentration and metabolism, and changes in the expression or activities of key proteins. In established cancers, autophagy increases metabolism to inactivate drugs and support drug resistance [59,60,61]. Moreover, autophagy actively participates in maintaining genomic integrity, as well as in repair processes [62], and the deregulation of DNA repair pathways is associated with the initiation, progression and resistance of cancer cells by promoting genomic instability and mutation [63]. Indeed, RNA-seq analysis revealed that several metabolic pathways and gene repair pathways were upregulated in PDAC PDOs from patients who were LSD PPV carriers, suggesting enhanced metabolic utility and dysregulation of DNA repair system in established PDAC with lysosome dysfunction. Further studies are required to generate direct evidence supporting our interpretations.
Moreover, the LSD PPV carriers tended to be less sensitive to drugs in both PDAC cell lines and human PDAC PDOs. The expression of S100 calcium-binding protein P (S100P) has been associated with drug resistance, metastasis, and poor clinical outcomes [64]. S100P promotes pancreatic cancer growth, survival, and invasion, and its intracellular levels affect resistance to 5-fluorouracil treatment in vitro [65]. In addition, high expression of the long non-coding RNA actin filament-associated protein 1 antisense RNA1 (AFAP1-AS1) is associated with poor survival and short-term recurrence in PDAC. Knockdown of AFAP1-AS1 attenuated PDAC cell proliferation, migration, and invasion [66]. Apolipoprotein M (APOM) can suppress the proliferation and invasion of hepatocellular carcinoma [67], breast cancer [68], and larynx carcinoma [69]. These results are in agreement with upregulation of S100P and AFAP1-AS1 and downregulation of APOM in the RNA-seq data of PDOs carrying LSD PPVs.
From a genetic perspective, our results imply that the heterozygous status of LSD genes is related to cancer development. Accordingly, we checked whether LSD genes follow Knudson’s two-hit hypothesis [70] similar to most CPGs. However, we did not observe any significant loss of heterozygosity in PDAC patients with LSD PPV from CNA analysis using matched tumor sequencing data (data not shown). Our results suggested that lysosomes and autophagy play biphasic roles in the initiation and progression of cancer. Complete loss of autophagy is not suitable for cancer cell survival in the PDAC microenvironment. Notably, we provide evidence showing that pathogenic variants in LSD genes can lead to the development of adult-onset chronic diseases such as PDAC other than the classic LSD phenotype, similar to that observed for the development of Parkinson’s disease and heterozygote carriers of GBA rare variants [15].
Finally, many attempts to develop therapeutics targeting autophagy in PDAC are underway, including a clinical trial (NCT03825289) exploring concurrent inhibition of the RAS-RAF-MEK-ERK pathway and autophagy [71,72,73,74,75]. The use of nanodrugs in PDAC have also been actively investigated [76, 77]. As all these modalities are closely related to lysosomal function, our results can act as a cornerstone for developing novel therapeutics for PDAC.
Conclusion
In conclusion, we showed the fact that genetically defined lysosomal dysfunction is frequently observed in young-onset PDAC. Lysosomal dysfunction may contribute to PDAC development by impairing autophagolysosome activity. In established PDAC, lysosomal dysfunction is closely associated with an increased autophagic flux and upregulated metabolism. This is the first study to report a relationship between lysosome dysfunction and PDAC development. Consistency in our results obtained using multiple methods, including DNA sequencing, RNA-seq analysis, and experiments using knockouts in mouse organoids and patient-derived human organoid, indicate the validity of our findings. We believe that our observations will act as a cornerstone for research on the role of lysosome dysfunction in PDAC. However, as our study has focused on a single ethnic group (Korean), generalization of our results to PDAC in non-Asian populations should be avoided till further investigations are performed. In addition, future studies should aim to understand the complete mechanism underlying lysosomal dysfunction in PDAC carcinogenesis, as we observed increased proliferation of pancreatic cells, and not complete carcinogenesis, after GALC knockdown.
Availability of data and materials
Sequencing BAM files were deposited in the Sequence Read Archive (SRA Bioproject Accession # PRJNA929903).
Abbreviations
- CFNC:
-
Cancer-free normal control
- CPG:
-
Cancer-predisposing gene
- CQ:
-
Chloroquine
- DEG:
-
Differentially expressed gene
- GALC:
-
Galactosylceramidase
- GO:
-
Gene ontology
- GSEA:
-
Gene set enrichment analysis
- KEGG:
-
Kyoto Encyclopedia of Genes and Genomes
- KO:
-
Knockout
- IC50:
-
Half maximal inhibitory concentration
- LSD:
-
Lysosomal storage diseases
- OS:
-
Overall survival
- PDAC:
-
Pancreatic ductal adenocarcinoma
- PDO:
-
Patient-derived organoid
- PPV:
-
Putative pathogenic variant
- PTV:
-
Protein truncating variant
- RNA-seq:
-
RNA sequencing
- VCF:
-
Variant calling format
- WES:
-
Whole exome sequencing
References
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer Statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.
Strobel O, Neoptolemos J, Jäger D, Büchler MW. Optimizing the outcomes of pancreatic cancer surgery. Nat Rev Clin Oncol. 2019;16(1):11–26.
Solomon S, Das S, Brand R, Whitcomb DC. Inherited pancreatic cancer syndromes. Cancer J. 2012;18(6):485–91.
Grover S, Syngal S. Hereditary pancreatic cancer. Gastroenterology. 2010;139(4):1076–80, 80.e1-1.
Earl J, Galindo-Pumariño C, Encinas J, Barreto E, Castillo ME, Pachón V, et al. A comprehensive analysis of candidate genes in familial pancreatic cancer families reveals a high frequency of potentially pathogenic germline variants. EBioMedicine. 2020;53:102675.
Hu C, Hart SN, Polley EC, Gnanaolivu R, Shimelis H, Lee KY, et al. Association between inherited germline mutations in cancer predisposition genes and risk of pancreatic cancer. JAMA. 2018;319(23):2401–9.
Hu C, Hart SN, Bamlet WR, Moore RM, Nandakumar K, Eckloff BW, et al. Prevalence of pathogenic mutations in cancer predisposition genes among pancreatic cancer patients. Cancer Epidemiol Biomarkers Prev. 2016;25(1):207–11.
Shindo K, Yu J, Suenaga M, Fesharakizadeh S, Cho C, Macgregor-Das A, et al. Deleterious germline mutations in patients with apparently sporadic pancreatic adenocarcinoma. J Clin Oncol. 2017;35(30):3382–90.
Parenti G, Andria G, Ballabio A. Lysosomal storage diseases: from pathophysiology to therapy. Annu Rev Med. 2015;66:471–86.
Yang C, Wang X. Lysosome biogenesis: regulation and functions. J Cell Biol. 2021;220(6):e202102001.
Goh LK, Sorkin A. Endocytosis of receptor tyrosine kinases. Cold Spring Harb Perspect Biol. 2013;5(5):a017459.
Kroemer G, Jaattela M. Lysosomes and autophagy in cell death control. Nat Rev Cancer. 2005;5(11):886–97.
Shin J, Kim D, Kim HL, Choi M, Koh Y, Yoon SS. Oncogenic effects of germline variants in lysosomal storage disease genes. Genet Med. 2019;21(12):2695–705.
Amaravadi R, Kimmelman AC, White E. Recent insights into the function of autophagy in cancer. Genes Dev. 2016;30(17):1913–30.
Sidransky E, Lopez G. The link between the GBA gene and parkinsonism. Lancet Neurol. 2012;11(11):986–98.
White E. The role for autophagy in cancer. J Clin Invest. 2015;125(1):42–6.
Commisso C, Davidson SM, Soydaner-Azeloglu RG, Parker SJ, Kamphorst JJ, Hackett S, et al. Macropinocytosis of protein is an amino acid supply route in Ras-transformed cells. Nature. 2013;497(7451):633–7.
Perera RM, Stoykova S, Nicolay BN, Ross KN, Fitamant J, Boukhali M, et al. Transcriptional control of autophagy-lysosome function drives pancreatic cancer metabolism. Nature. 2015;524(7565):361–5.
Yang S, Wang X, Contino G, Liesa M, Sahin E, Ying H, et al. Pancreatic cancers require autophagy for tumor growth. Genes Dev. 2011;25(7):717–29.
Wu WK, Coffelt SB, Cho CH, Wang XJ, Lee CW, Chan FK, et al. The autophagic paradox in cancer therapy. Oncogene. 2012;31(8):939–53.
Poplin R, Ruano-Rubio V, DePristo MA, Fennell TJ, Carneiro MO, Auwera GAV, et al. Scaling accurate genetic variant discovery to tens of thousands of samples. BioRxiv. 2018. https://doi.org/10.1101/201178.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754–60.
Scheinin I, Sie D, Bengtsson H, van de Wiel MA, Olshen AB, van Thuijl HF, et al. DNA copy number analysis of fresh and formalin-fixed specimens by shallow whole-genome sequencing with identification and exclusion of problematic regions in the genome assembly. Genome Res. 2014;24(12):2022–32.
Amemiya HM, Kundaje A, Boyle AP. The ENCODE blacklist: identification of problematic regions of the genome. Sci Rep. 2019;9(1):9354.
Derrien T, Estellé J, Marco Sola S, Knowles DG, Raineri E, Guigó R, et al. Fast computation and applications of genome mappability. PLoS ONE. 2012;7(1):e30377.
Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38(16):e164.
McLaren W, Gil L, Hunt SE, Riat HS, Ritchie GRS, Thormann A, et al. The ensembl variant effect predictor. Genome Biol. 2016;17(1):122.
Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alföldi J, Wang Q, et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature. 2020;581(7809):434–43.
Landrum MJ, Lee JM, Riley GR, Jang W, Rubinstein WS, Church DM, et al. ClinVar: public archive of relationships among sequence variation and human phenotype. Nucleic Acids Res. 2014;42(Database):D980–5.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinfor. 2011;12:323.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Liberzon A, Subramanian A, Pinchback R, Thorvaldsdóttir H, Tamayo P, Mesirov JP. Molecular signatures database (MSigDB) 3.0. Bioinformatics. 2011;27(12):1739–40.
Barretina J, Caponigro G, Stransky N, Venkatesan K, Margolin AA, Kim S, et al. The cancer cell line encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature. 2012;483(7391):603–7.
Buscail L, Bournet B, Cordelier P. Role of oncogenic KRAS in the diagnosis, prognosis and treatment of pancreatic cancer. Nat Rev Gastroenterol Hepatol. 2020;17(3):153–68.
Loos B, du Toit A, Hofmeyr JH. Defining and measuring autophagosome flux—concept and reality. Autophagy. 2014;10(11):2087–96.
Liu WJ, Ye L, Huang WF, Guo LJ, Xu ZG, Wu HL, et al. p62 links the autophagy pathway and the ubiqutin-proteasome system upon ubiquitinated protein degradation. Cell Mol Biol Lett. 2016;21:29.
Jung CH, Ro SH, Cao J, Otto NM, Kim DH. mTOR regulation of autophagy. FEBS Lett. 2010;584(7):1287–95.
Leeman DS, Hebestreit K, Ruetz T, Webb AE, McKay A, Pollina EA, et al. Lysosome activation clears aggregates and enhances quiescent neural stem cell activation during aging. Science. 2018;359(6381):1277–83.
Huang KL, Mashl RJ, Wu Y, Ritter DI, Wang J, Oh C, et al. Pathogenic germline variants in 10,389 adult cancers. Cell. 2018;173(2):355-70.e14.
Hanahan D. Hallmarks of cancer: new dimensions. Cancer Discov. 2022;12(1):31–46.
Risch N, Burchard E, Ziv E, Tang H. Categorization of humans in biomedical research: genes, race and disease. Genome Biol. 2002;3(7):comment2007.
Kittles RA, Weiss KM. Race, ancestry, and genes: implications for defining disease risk. Annu Rev Genomics Hum Genet. 2003;4:33–67.
Chen W, Wang S, Tithi SS, Ellison DW, Schaid DJ, Wu G. A rare variant analysis framework using public genotype summary counts to prioritize disease-predisposition genes. Nat Commun. 2022;13(1):2592.
Miyashita A, Kikuchi M, Hara N, Ikeuchi T. Genetics of Alzheimer’s disease: an East Asian perspective. J Hum Genet. 2023;68(3):115–24.
Peljto AL, Selman M, Kim DS, Murphy E, Tucker L, Pardo A, et al. The MUC5B promoter polymorphism is associated with idiopathic pulmonary fibrosis in a Mexican cohort but is rare among Asian ancestries. Chest. 2015;147(2):460–4.
Jang MA, Chung JW, Yeon JY, Kim JS, Hong SC, Bang OY, et al. Frequency and significance of rare RNF213 variants in patients with adult moyamoya disease. PLoS ONE. 2017;12(6):e0179689.
Jeon S, Bhak Y, Choi Y, Jeon Y, Kim S, Jang J, et al. Korean Genome Project: 1094 Korean personal genomes with clinical information. Sci Adv. 2020;6(22):eaaz7835.
Xu C, Sakai N, Taniike M, Inui K, Ozono K. Six novel mutations detected in the GALC gene in 17 Japanese patients with Krabbe disease, and new genotype-phenotype correlation. J Hum Genet. 2006;51(6):548–54.
Li J, Chen X, Kang R, Zeh H, Klionsky DJ, Tang D. Regulation and function of autophagy in pancreatic cancer. Autophagy. 2021;17(11):3275–96.
Perera RM, Bardeesy N. Pancreatic cancer metabolism: breaking it down to build it back up. Cancer Discov. 2015;5(12):1247–61.
Guo JY, Xia B, White E. Autophagy-mediated tumor promotion. Cell. 2013;155(6):1216–9.
Noguchi M, Hirata N, Tanaka T, Suizu F, Nakajima H, Chiorini JA. Autophagy as a modulator of cell death machinery. Cell Death Dis. 2020;11(7):517.
Kimmelman AC. The dynamic nature of autophagy in cancer. Genes Dev. 2011;25(19):1999–2010.
Autophagy-deficient pancreatic cancer cells depend on macropinocytosis. Cancer Discov. 2021;11(5):OF28.
Su H, Yang F, Fu R, Li X, French R, Mose E, et al. Cancer cells escape autophagy inhibition via NRF2-induced macropinocytosis. Cancer Cell. 2021;39(5):678-93.e11.
Chang H, Zou Z. Targeting autophagy to overcome drug resistance: further developments. J Hematol Oncol. 2020;13(1):159.
Smith AG, Macleod KF. Autophagy, cancer stem cells and drug resistance. J Pathol. 2019;247(5):708–18.
Sui X, Chen R, Wang Z, Huang Z, Kong N, Zhang M, et al. Autophagy and chemotherapy resistance: a promising therapeutic target for cancer treatment. Cell Death Dis. 2013;4(10):e838.
Ambrosio S, Majello B. Autophagy roles in genome maintenance. Cancers. 2020;12(7):1793.
Li LY, Guan YD, Chen XS, Yang JM, Cheng Y. DNA repair pathways in cancer therapy and resistance. Front Pharmacol. 2020;11:629266.
Arumugam T, Logsdon CD. S100P: a novel therapeutic target for cancer. Amino Acids. 2011;41(4):893–9.
Arumugam T, Simeone DM, Van Golen K, Logsdon CD. S100P promotes pancreatic cancer growth, survival, and invasion. Clin Cancer Res. 2005;11(15):5356–64.
Ye Y, Chen J, Zhou Y, Fu Z, Zhou Q, Wang Y, et al. High expression of AFAP1-AS1 is associated with poor survival and short-term recurrence in pancreatic ductal adenocarcinoma. J Transl Med. 2015;13:137.
Yu M, Pan L, Sang C, Mu Q, Zheng L, Luo G, et al. Apolipoprotein M could inhibit growth and metastasis of SMMC7721 cells via vitamin D receptor signaling. Cancer Manag Res. 2019;11:3691–701.
Zhou Y, Yao S, Yu M, Wei J, Fang Q, Xu N, et al. The effects and possible mechanism of action of apolipoprotein M on the growth of breast cancer cells. Mol Biol Rep. 2022;49(2):1171–9.
Xue H, Yu M, Zhou Y, Zhang J, Mu Q, Chen T, et al. Apolipoprotein M inhibits proliferation and migration of larynx carcinoma cells. Sci Rep. 2020;10(1):19424.
Knudson AG Jr. Mutation and cancer: statistical study of retinoblastoma. Proc Natl Acad Sci USA. 1971;68(4):820–3.
Samaras P, Tusup M, Nguyen-Kim TDL, Seifert B, Bachmann H, von Moos R, et al. Phase I study of a chloroquine-gemcitabine combination in patients with metastatic or unresectable pancreatic cancer. Cancer Chemother Pharmacol. 2017;80(5):1005–12.
Boone BA, Bahary N, Zureikat AH, Moser AJ, Normolle DP, Wu WC, et al. Safety and biologic response of pre-operative autophagy inhibition in combination with gemcitabine in patients with pancreatic adenocarcinoma. Ann Surg Oncol. 2015;22(13):4402–10.
Karasic TB, O’Hara MH, Loaiza-Bonilla A, Reiss KA, Teitelbaum UR, Borazanci E, et al. Effect of gemcitabine and nab-paclitaxel with or without hydroxychloroquine on patients with advanced pancreatic cancer: a phase 2 randomized clinical trial. JAMA Oncol. 2019;5(7):993–8.
Kinsey CG, Camolotto SA, Boespflug AM, Guillen KP, Foth M, Truong A, et al. Protective autophagy elicited by RAF→MEK→ERK inhibition suggests a treatment strategy for RAS-driven cancers. Nat Med. 2019;25(4):620–7.
Bryant KL, Stalnecker CA, Zeitouni D, Klomp JE, Peng S, Tikunov AP, et al. Combination of ERK and autophagy inhibition as a treatment approach for pancreatic cancer. Nat Med. 2019;25(4):628–40.
López-Méndez TB, Sánchez-Álvarez M, Trionfetti F, Pedraz JL, Tripodi M, Cordani M, et al. Nanomedicine for autophagy modulation in cancer therapy: a clinical perspective. Cell Biosci. 2023;13(1):44.
Ma Z, Li J, Lin K, Ramachandran M, Zhang D, Showalter M, et al. Pharmacophore hybridisation and nanoscale assembly to discover self-delivering lysosomotropic new-chemical entities for cancer therapy. Nat Commun. 2020;11(1):4615.
Acknowledgements
None.
Funding
This work was supported by the National Research Foundation of Korea funded by the Korean Government (MSIT) (Grant No. NRF-2021R1A2C3005360) (YK) and the Ministry of Health & Welfare, Republic of Korea (Grant No. HI18C1876) (SSY). This study was supported by the Future Medicine 20 × 30 Project of the Samsung Medical Center (Grant No. SMX1230041, SMO1230021) and a Samsung Medical Center Research and Development Grant (Grant No. SMO1230661) (JKP).
Author information
Authors and Affiliations
Contributions
YK (conceptualization: lead; funding acquisition: lead; project administration: lead; supervision: lead; writing—review and editing: lead). HK (data curation: lead; formal analysis: lead; investigation: lead; methodology: lead; writing the original draft). SYJ (data curation: lead; investigation: lead; methodology: lead; writing the original draft). SS (data curation: lead; formal analysis: lead; investigation: lead; methodology: lead; writing the original draft: lead). YHC (data curation: lead; formal analysis: lead; writing the original draft: lead). HRK (data curation: supporting; formal analysis: supporting). BM (investigation: supporting). JB (formal analysis: supporting). JH (formal analysis: supporting). D-YS (formal analysis: supporting). SP (formal analysis: supporting). KHL (resources: supporting). KTL (resources: supporting). JKL (resources: supporting). DP (formal analysis: supporting). S-HL (investigation: supporting). J-YJ (project administration: lead; supervision: lead; writing—review and editing: lead). HL (project administration: lead; supervision: lead; writing—review and editing: lead). J-AK (project administration: lead; supervision: lead; writing—review and editing: lead). S-SY (funding acquisition: lead; project administration: lead; supervision: lead). JKP (conceptualization: lead; funding acquisition: lead; project administration: lead; supervision: lead; writing—review and editing: lead).
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
All patients and volunteers provided written informed consent, and all specimens were collected according to the institutional review board regulations and approval (IRB No. 2018-12-065, 1705-031-852).
Consent for publication
Not applicable.
Competing interests
The authors have no potential conflicts of interest to disclose.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1
: Figure S1. Onset age in PDAC patients who were CPG or LSD carriers.
Additional file 2
: Figure S2. Evaluation of Galc knockout in mouse pancreatic organoids.
Additional file 3
: Figure S3. Increased proliferation in mouse KrasG12D/Galc knockout pancreatic organoids.
Additional file 4
: Table S1. Forty-two Lysosomal storage disorder-associated genes
Additional file 5
: Table S2. List of antibodies and drugs used in the study
Additional file 6
: Table S3. Pathogenic variants identified in the Korean PDAC and normal control.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Koh, Y., Kim, H., Joo, S.Y. et al. Genetic assessment of pathogenic germline alterations in lysosomal genes among Asian patients with pancreatic ductal adenocarcinoma. J Transl Med 21, 730 (2023). https://doi.org/10.1186/s12967-023-04549-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12967-023-04549-x