RETRACTED ARTICLE: ITGA3 acts as a purity-independent biomarker of both immunotherapy and chemotherapy resistance in pancreatic cancer: bioinformatics and experimental analysis

Contribution of integrin superfamily genes to treatment resistance remains uncertain. Genome patterns of thirty integrin superfamily genes were analyzed of using bulk and single-cell RNA sequencing, mutation, copy number, methylation, clinical information, immune cell infiltration, and drug sensitivity data. To select the integrins that are most strongly associated with treatment resistance in pancreatic cancer, a purity-independent RNA regulation network including integrins were constructed using machine learning. The integrin superfamily genes exhibit extensive dysregulated expression, genome alterations, epigenetic modifications, immune cell infiltration, and drug sensitivity, as evidenced by multi-omics data. However, their heterogeneity varies among different cancers. After constructing a three-gene (TMEM80, EIF4EBP1, and ITGA3) purity-independent Cox regression model using machine learning, ITGA3 was identified as a critical integrin subunit gene in pancreatic cancer. ITGA3 is involved in the molecular transformation from the classical to the basal subtype in pancreatic cancer. Elevated ITGA3 expression correlated with a malignant phenotype characterized by higher PD-L1 expression and reduced CD8+ T cell infiltration, resulting in unfavorable outcomes in patients receiving either chemotherapy or immunotherapy. Our findings suggest that ITGA3 is an important integrin in pancreatic cancer, contributing to chemotherapy resistance and immune checkpoint blockade therapy resistance. Graphical abstract Supplementary Information The online version contains supplementary material available at 10.1007/s10142-023-01122-z.


Introduction
Pancreatic cancer is a highly aggressive and lethal tumor known for its resistance to treatment with a 5-year survival rate of approximately 10%.Currently, pancreatic cancer is the seventh leading cause of cancer-related deaths worldwide.However, by 2030, it has been projected to be the second leading cause of cancerrelated deaths in the USA (Badheeb et al. 2022;Cai et al. 2021).In China, both the incidence and mortality rates of pancreatic cancer have recently been on the rise (Chen et al. 2016).Radical resection is the most important approach for achieving long-term survival; however, patients in China are often diagnosed at an advanced stage and experience a short period of disease-free survival after surgery (Cai et al. 2021).
Resistance to chemotherapy, radiation therapy, and targeted therapies is a major challenge in pancreatic cancer treatment.Moreover, patients are naturally resistant to immunotherapy (Wu et al. 2022), with only 10% of patients classified as immune hot patients (Cao et al. 2021).Understanding the genomic factors contributing to treatment failure is crucial in pancreatic cancer research and can aid in determining personalized therapies based on individual genome alterations.Pancreatic cancer is characterized by a high frequency of genomic alterations, with common driver gene mutations including KRAS, TP53, SMAD4, and CDKN2A.However, DNA molecular subtyping could not reflect the complexity of pancreatic cancer (Waddell et al. 2015).Several subtyping schemas have been proposed, including Moffitt's two subtypes (Moffitt et al. 2015), Collisson's three subtypes (Collisson et al. 2011), Bailey's four subtypes (Bailey et al. 2016), and Puleo's five subtypes (Puleo et al. 2018).However, these classification approaches are complex and can be negatively affected by tumor purity or sampling location, possibly hindering accurate classification (Cancer Genome Atlas Research Network.Electronic address and Cancer Genome Atlas Research 2017; Collisson et al. 2019;Waddell et al. 2015).Therefore, these classifiers may not always be clinically implemented.Thus, it is important to identify a marker that is both purity-independent and easy to use for classifying molecular subtypes and guiding treatment decisions.
In the era of single-cell sequencing, the significance of the integrin superfamily genes in complex cellular communication within the tumor microenvironment has been rediscovered.Integrin superfamily genes play an important role in cancer (Miyamoto et al. 2022;

R E T R A C T E D A R T I C L E
196 Page 4 of 29 alterations, methylation, clinical information, and immune cell infiltration.For pan-cancer comparison of gene expression, processed transcripts per million data were used by combining TCGA and Genotype-Tissue Expression (GTEx) data (Consortium et al. 2015), which were downloaded from UCSC Xena (Goldman et al. 2017).
For the purity-independent network construction, patient cohorts with both pancreatic cancer and normal adjacent samples were included.Samples other than those identified as outliers during quality control checking were included for regulation prediction, obtained from TCGA datasets and four pathological pancreatic adenocarcinomas datasets: GSE32688, (Donahue et al. 2012), GSE41372, (Frampton et al. 2014), GSE119794, (Lin et al. 2019), and GSE43797 (Park et al. 2014).Tables 1 and 2 present the clinical and quality control information.All five existing studies had provided expression data for mRNAs, long noncoding RNAs (lncRNAs), and microRNAs (miRNAs) in paired pancreatic cancer and normal tissues.Samples with partial or whole neuroendocrine carcinoma were excluded from the analysis (Nagtegaal et al. 2020).After analyzing the prediction results, we only considered the RNA regulatory relationships observed in the cancer tissues obtained from TCGA datasets, GSE32688, and GSE119794 with sufficient cancer samples for constructing the purity-independent network.We excluded GSE41372 and GSE43797 because of insufficient cancer samples.
The patient cohorts of the five mRNA datasets, TCGA pancreatic ductal adenocarcinoma (PAAD), E-MTAB-6134 from ArrayExpress (Athar et al. 2019;Puleo et al. 2018), International Cancer Genome Consortium (ICGC)-PACA-AU (Bailey et al. 2016;Perez-Mancera et al. 2012) and ICGC-PACA-CA from ICGC (https:// daco.icgc.org/), and GSE21501 and GSE71729 from Gene Expression Omnibus were used to develop the prognostic model.E-MTAB-6134 was used for external validation, and the immunotherapy cohort (Parikh et al. 2021) was used to determine the prognostic values of the genes to guide the immunotherapy treatment.We used data that met the following criteria: availability of mRNA expression profiles of primary pancreatic cancer samples (excluding neuroendocrine carcinoma samples), patients with sufficient clinical data, and patients with at least 30 days of follow-up and survival status information.Table 3 presents the clinical information of the datasets used to construct and validate the purity-independent model.The median overall survival (OS) of the construction group was 540 days, with 1-, 2-, 3-, and 5-year survival rates of 68.7%, 33.7%, 22.9%, and 10.1%, respectively.The median OS of the validation group was 795 days, with 1-, 2-, 3-, and 5-year survival rates of 78.3%, 52.0%, 38.9%, and 27.5%, respectively.
Pancreatic cancer patient cohorts of single-cell RNA sequencing (RNA-seq) studies by Peng et al. (M Li et al

Cell cultures and RNA-seq
Primary cultures of human pancreatic cancer-associated fibroblasts (CAFs) were obtained from human pancreatic cancer tissues via pancreatectomy.The tissue was cut and incubated with dispase at 4 °C to separate the fibroblasts from other cells.Furthermore, we obtained other cell lines, including HPDE6C7, PANC1, BxPC3, AsPC1, and MIA-PACA-2, from the American Type Culture Collection (Manassas, VA, USA).RNA-seq was performed by the Annoroad Gene Company (Beijing, China) using the Illumina HiSeq X-ten system.For more detailed information, please refer to the Methods section of our previous work (Du et al. 2022).

Mouse model and single-cell sequencing
We purchased four C57BL/6J mice (12-week-old males) from GemPharmatech Co., Ltd.(Nanjing, China) and housed them under standard pathogen-free conditions.The mice were anesthetized by administering 5% isoflurane in the air, and each mouse received 2% isoflurane via a face mask.Cultured KPC cells (5 × 10 5 ; LSL-Kras G12D/+ ; LSL-Trp53 R172H/+ ; Pdx1-Cre) were injected into the tail of pancreas, and they were returned to their cages after recovery.The mice were euthanized after 4 weeks under anesthesia, and the tumors were removed for single-cell sequencing that was performed by the Annoroad Gene Company using NovaSeq 6000.

Preparation of single-cell samples
Mouse pancreatic tissue pieces were washed in phosphate-buffered saline (PBS) and digested in a mixture of collagenase (2 mg/mL), calcium chloride (8.5 mM), and DNA enzyme (50 μg/mL) at 4 °C for 30 min to 1 h.Subsequently, the suspension was filtered using a 70-μm strainer to remove undigested cell clumps, and the resulting liquid was centrifuged at 500 g for 7 min to obtain cell pellets.Dead cells were removed using the Dead Cell Removal Kit (Miltenyi Biotec), and the cells were counted using a LUNA counter.The supernatant was carefully discarded while retaining the precipitate after centrifugation.The cell pellets were resuspended in 1 mL of 1X PBS containing 0.04% bovine serum albumin and washed twice with PBS.
The cell concentration was determined using automatic cytometry, and the sample volume was calculated based on the optimal cell sampling concentration provided by 10× Genomics official website and the target capture number.Once we obtained the desired cell suspension, the samples were immediately cooled on ice for subsequent gel beads-in-emulsion (GEM) preparation and reverse transcription.

Library construction and single-cell sequencing
During the library preparation process, gel beads and cells containing a poly-T primer sequence linked with a cell barcode and a unique molecular identifier were encapsulated in oil drops, forming GEM.The cells were lysed within the GEM, and the released mRNA was subjected to reverse transcription, leading to the formation of full-length cDNA sequences.The oil drops were then crushed and purified, and the cDNA library was subjected to polymerase chain reaction amplification, with sequencing primers being connected to complete the library construction.Following library construction, the library concentration was determined using the Qubit® 3.0 Fluorometer (Life Technologies, Carlsbad, CA, USA), whereas the distribution of library product fragments was determined using the Agilent 2100 High Sensitivity DNA Assay Kit (Agilent Technologies, Santa Clara, CA, USA).The library product fragments were required to be > 10 nmol/L.The resulting libraries were sequenced on an Illumina NovaSeq 6000 platform using a 2 × 150cycle paired-end run.

Analysis and reanalysis of single-cell sequencing
We used the Cell Ranger software (10× Genomics) to perform unique molecular identifier and barcode counting for the mice data.The data were filtered and aligned to the GRCm38 (mm10) reference genome to generate the feature file, barcode file, and matrix.Subsequent analyses of the single-cell data were conducted using Seurat v4.3.0 (Hao et al. 2021) The mutation data were collected from TCGA.Gene set cancer analysis (GSCA) (C J Liu et al. 2023a) includes information regarding missense mutations, nonsense mutation, frame shift insertions, splice sites, frame shift deletions, in-frame deletions, and in-frame insertions as deleterious mutations for analysis.cBioPortal (Cerami et al. 2012;Gao et al. 2013) contains information regarding nonsilent mutations, heterozygous copy number variations, and structural variations as genomic alterations.Copy number data were processed using GISTICS 2.0 (Mermel et al. 2011) and classified into five statuses based on the standards of cBioPortal: no alteration, heterozygous amplification, heterozygous deletions, homozygous amplification, and homozygous deletions.DNA methylation data were compared using β values.

Quality control and batch correction before RNA regulatory network construction
Before analysis, the ComBat function in the sva package was used to adjust for batch effects (Leek et al. 2012).Samples outside the ellipses in the principal component analysis (PCA) plots were considered outliers, set at four standard deviations of the mean (99.9937% of values lie within one).Additional PCAs were conducted iteratively until no outliers were identified, and samples within the ellipses were selected for further network construction.PCA results were compared based on sample types (tumor vs. normal), studies (TCGA, GSE32688, GSE41372, GSE119794, and GSE43797), and sample types in individual studies.Based on these three factors, mRNA profiles before adjusting for batch effects are plotted in Supplementary Fig. 1A-C, those after adjusting for batch effects are in Supplementary Fig. 1D-F, and those after removing all outliers are in Supplementary Fig. 1G-I.Similarly, miRNA profiles before adjusting for batch effects are plotted in Supplementary Fig. 2A-C, those after adjusting for batch effects are in Supplementary Fig. 2D-F, and those after removing all outliers are in Supplementary Fig. 2G-I.We excluded two tumor samples and one normal sample from TCGA and one patient from GSE32688 after conducting PCA.Finally, 165 tumor samples from TCGA, 24 from GSE32688, 6 from GSE41372, 10 from GSE119794, and 6 from GSE43797 were retained for RNA regulation network prediction.Moreover, 3 control samples from TCGA, 7 from GSE32688, 6 from GSE41372, 10 from GSE119794, and 5 from GSE43797 were selected for analyses (Table 2).Although TCGA-PAAD dataset had the highest number of tumor samples, its purity limit was lower (≥ 30%) than that of the other datasets (Table 3).The expression profiles of the patients who did not meet the study standards were considered outliers in the pre-analysis quality control step and were hence excluded.

Purity-independent network construction
A purity-independent gene network was constructed using experimentally verified RNA regulatory pairs available in StarBase (Li et al. 2014), which implemented the competing endogenous RNA (ceRNA) mechanism among mRNAs, lncRNAs, and miRNAs.Theoretically, we expected the conservation of the actual RNA regulations at different tumor purities.Unverified relationships were eliminated to maintain network reliability.To ensure the relationship between the RNA regulatory pairs that have shared miRNAs, the expression profiles were subjected to correlation analyses and a hypergeometric test in the R package GDCRNAtools (Li et al. 2018).Only one RNA regulation pair with lncRNA-mRNA correlation P < 0.05, correlation coefficient > 0, and hypergeometric test P < 0.05 was considered to be preserved.The resulting coexpression networks were visualized using Cytoscape (Shannon et al. 2003).RNA relationships were calculated for tumor samples (Supplementary Fig. 3A), normal samples (Supplementary Fig. 3B), and all samples (Supplementary Fig. 3C).Significant differences were observed in the prediction results not only among the datasets but also among the sample types included.Only 10 common lncRNA-mRNA pairs were recognized in at least three tumor datasets and three control datasets (Supplementary Fig. 3D).No common predictive pairs were found if the criteria were set to at least four tumor datasets and four control datasets (Supplementary Fig. 3E).

Selection of prognosis-related genes and evaluation of the risk model
We observed differences in ceRNA regulation between cancerous and healthy tissues (Supplementary Fig. 3D and Supplementary Fig. 3E), which were in agreement with the results of prior studies (Li et al. 2018;Loganathan and Doss C 2023).Furthermore, we observed that the total number (red) of predicted ceRNA regulatory pairs significantly decreased with sample number (Supplementary Fig. 3F).Therefore, the prediction results from the tumor samples from GSE41372 and GSE43797 with only less than 10 tumor samples were removed.The prediction results from the tumor samples among the remaining three datasets (TCGA, GSE32688, and GSE119794) were observed.Finally, 201 lncRNA-mRNA pairs were considered to be preserved pairs for constructing the prognostic model (Supplementary Fig. 3F).The combined datasets (TCGA, ICGC-PACA-CA, ICGC-PACA-AU, GSE71729, and GSE21501) were used for constructing the model, and E-MTAB-6134 was used for external validation.To select the most useful prognostic mRNA models, least absolute shrinkage and selection operator (LASSO) (Gui and Li 2005;Tibshirani 1996;Tibshirani 1997) and Cox regression analyses were performed using the "glmnet" package in R (Friedman et al. 2010) for model construction, and the model with a lambda value of lambda.1sewas used.Parameters (λ) in the LASSO model were selected using three-fold cross-validation.Blasso, a variant of LASSO, was used to determine the frequency of the variables being selected among different seeds for steady selection of the prognostic genes (Zeng et al. 2020).

Survival analyses
We evaluated the survival data using Kaplan-Meier survival analyses (log-rank method).The optimal cut-off values were determined using the "survminer" package, which divided the patients into high-risk and low-risk groups based on the association between the survival time and risk mode in the cohorts (Kassambara et al. 2017).We conducted Cox regression model analyses as previously described, and P < 0.05 was considered statistically significant for all analyses.We used time-dependent area under the curve (AUC) and decision curve analyses to represent the advantages of the prognostic model.

Statistical analyses
Quantitative data are presented as the median ± interquartile range.Numeric data were analyzed using the Student's t test and Wilcoxon's rank-sum tests, whereas classified variables were analyzed using the chi-square and Fisher's exact tests.Significance of enrichment results was determined using a one-sided Fisher's exact test and corrected for multiple hypotheses using a Benjamini-Hochberg false discovery rate (FDR).Tests with a FDR or P value < 0.05 were considered significant.All statistical tests were performed in R v4.2.1 (http:// www.R-proje ct.org/).Correlation analyses were conducted as mentioned in the figures.

Pan-cancer view of integrin superfamily genes
Thirty integrin superfamily genes were comprehensively analyzed to investigate their pan-cancer expression patterns, revealing the heterogeneous expression of most integrins between tumor and normal tissue across human cancer types (Fig. 1A).As determined using the Wilcoxon's test, the expression of all integrins in pancreatic cancer, except that of ITGA1, increased in cancer tissues (Fig. 1A).Notably, differential expression analysis protocols may generate different results even when using the same datasets, as demonstrated by GEPIA2, wherein ITGA2B, ITGA7, ITGA8, ITGA10, ITGAD, and ITGB1BP2 did not exhibit significant difference (Supplementary Fig. 4A) (Du et al. 2022).Dysregulation of integrins is a characteristic feature of cancer, but the expression of specific genes varies with cancer and cell types.Single-nucleus RNA-seq analysis of untreated pancreatic cancer revealed heterogeneous distribution of integrins among cells.ITGA1, ITGA5, ITGA7, ITGA8, ITGA9, ITGA11, and ITGBL1 were found to be higher in stromal cells (fibroblasts, endothelial cells, nascent endothelial cells, and Schwann cells), while ITGA4, ITGAE, ITGAL, ITGAM, ITGAX, and ITGB2 were found to be higher in immune cells.ITGA2, ITGA3, ITGA10, ITGB4, and ITGB6 were more highly expressed in epithelial cells (tumor, atypical ductal, acinar, ductal, or endocrine cells), and ITGA2, ITGA3, and ITGB4 were highly expressed in malignant cells (tumor or atypical ductal cells) (Fig. 1B).
Genomic and epigenetic modifications such as DNA methylation can result in dysregulated gene expression, which can affect clinical outcomes in pancreatic cancer (X Li et al. 2022b).Figure 3A illustrates a pan-cancer overview of the deleterious mutations analyzed using GSCA.The patients with skin cutaneous melanoma (SKCM), uterine corpus endometrioid carcinoma (UCEC), colon adenocarcinoma (COAD), LUAD, STAD, lung squamous cell carcinoma (LUSC), BLCA, CESC, head and neck squamous cell carcinoma (HNSC), glioblastoma multiforme (GBM), and breast invasive carcinoma (BRCA) exhibited deleterious mutation rates of ≥ 10% among all cancers, whereas pancreatic cancer had a low rate.Copy number variants are believed to be an important source of transcriptomic dysregulation (Mielczarek et al. 2022).However, the cBioPortal oncoplot demonstrated that integrin superfamily genes have relatively low alteration rates (somatic mutations and homozygous alterations) ranging from 0 to 5% in pancreatic cancer (Supplementary Fig. 4B).Nevertheless, ITGA10 was identified as one of the most frequently altered integrin genes in patients with pancreatic cancer.When considering both homozygous and heterozygous alterations (homozygous deletions and homozygous amplifications), the number and types of alterations in integrin superfamily genes varied across cancers and most integrin superfamily genes were still not altered (Fig. 3B).Pan-cancer DNA methylation analysis revealed that ITGB3 was the most common hypomethylated integrin across cancer types, whereas ITGA8 was the most common hypermethylated integrin (Fig. 4A).We also observed a relationship between methylation and cell infiltration (Fig. 4B).ITGB6, ITGA10, ITGA3, ITGA6, ITGB3, and ITGA10 were positively correlated with infiltration scores, whereas ITGA5, ITGA9, ITGAM, ITGA7, ITGB1, ITGAE, ITGB2, ITGB7, ITGA4, ITGA8, ITGAL, ITGA11, and ITGAV were negatively correlated.Higher methylation of ITGA3, ITGA6, ITGB3, ITGB4, ITGB8, ITGBL1, and ITGB5 were associated with higher infiltration of CD8 + T cells into tumor sites, while ITGA9, ITGAM, ITGA7, ITGAE, ITGB2, ITGA4, ITGA8, and ITGAL had a negative impact in pancreatic cancer.TP53 is the second most commonly mutated gene in pancreatic cancer.The role of TP53 mutation in methylation regulation has Fig. 4 DNA methylation of integrins in human cancers.A Bubble plot comparing methylation levels (β values) of integrin superfamily genes in normal and tumor tissues across multiple cancer types from The Cancer Genome Atlas (TCGA).Cancer types with < 10 normal samples were excluded (n = 6560).Bubble color represents methylation level differences between tumor and normal samples, ranging from blue to red.P values were estimated using the two-tailed unpaired t test and adjusted using the false discovery rate (FDR).Bubbles with black borders indicate FDRs < 0.05.Bubble size indicates FDR significance, while color intensity increases with correlation.B Spearman's correlation between DNA methylation levels (β values) and cell infiltration degree for each cancer type from TCGA, evaluated using ImmuCellAI (n = 183).P values were adjusted using been extensively studied (Sadzeviciene et al. 2022;Xiao et al. 2022;Zhou et al. 2023).In this study, we compared methylation levels between patients with and without TP53 mutation.Our findings indicated that methylation levels not only differ between pancreatic cancer and normal tissues but also differ among patients with or without TP53 mutation.In particular, ITGAV, ITGB4, ITGB5, ITGA2, ITGA3, ITGAX, ITGAE, ITGAL, ITGB6, and ITGBL1 exhibited lower methylation levels, whereas ITGA4, ITGA8, ITGA11, ITGAL, ITGAM, and ITGB2 showed higher methylation levels (Fig. 4C).Methylation levels inversely correlated with the majority of integrin expression, suggesting the possibility of methylation-dependent regulation of integrin mRNA expression (Supplementary Fig. 4C).
From these analyses, we observed that integrins play a critical role in cancer omics, with pancreatic cancer exhibiting the highest number of prognostic integrins.Our focus lies in identifying an integrin with promising molecular and clinical applications in patients with pancreatic cancer.However, tumor purity poses a challenge to the analyses and selection process, given the high stromal content of pancreatic cancer.Therefore, we aimed to address this issue by developing a purity-independent prognostic model that incorporates integrins for clinical applications.

Purity-independent RNA regulation network uncovered the significance of ITGA3 in pancreatic cancer
Integrins are instrumental in the management of noncoding RNAs (Tan et al. 2018;Xu et al. 2022;Zheng et al. 2019).Therefore, we used experimentally confirmed RNA regulatory relationships to computationally predict these relationships under different levels of purity to establish a purity-independent prognostic model.The detailed protocol for constructing the purity-independent RNA regulatory network has been described in the "Materials and methods" section.Figure 5A presents a heatmap illustrating the total evidence and number of evidence for each lncRNA-mRNA or miRNA pair using prediction results from TCGA datasets, GSE32688 and GSE119794.We observed that the top 75% of the RNA regulatory relationships accounted for the majority of the regulatory pairs, suggesting that the shared prediction results are more credible.
To establish a purity-independent network, we utilized 201 mRNA-lncRNA pairs exclusively from calculation results obtained from three high-quality datasets using pancreatic cancer samples exclusively.After excluding genes absent in the survival cohorts, we obtained 194 qualified genes present in all prognostic datasets.These genes were used to establish the optimal prognosis model through LASSO regression analysis.The genes suitable for the Cox model were selected from the candidate genes using the LASSO algorithm.We performed data subsampling and determined tuning parameters based on the C-index, estimated via three-fold cross-validation.The lambda coefficients and C-indexes were plotted for different lambda values (Fig. 5 B and C).Finally, we selected lambda.1se(0.143) for further model selection and chose three genes-TMEM80, EIF4EBP1, and ITGA3-to construct the multiple Cox risk regression model for predicting OS (Fig. 5 B and C).
The optimal cut-off value for the risk score, which distinguishes the high-risk group from the low-risk group, was determined to be 0.02.In the construction group (Fig. 5 D  and F) and the validation groups (E-MTAB-6134, Fig. 5 E  and G), this value showed significance with P < 0.0001 and P = 0.0067, respectively.Higher scores in the model were associated with poor prognoses.The generalization of the model among each cohort was visualized at P < 0.05 (Supplementary Fig. 6A-E).The final mRNA-miRNA-lncRNA network genes, including three prognosis-related mRNAs, 15 miRNAs, and three lncRNAs, are shown in Fig. 6A.Time-dependent receiver operating characteristic (ROC) curves illustrated the 1-, 2-, 3-,4-and 5-year AUC values

R E T R A C T E D A R T I C L E
Page 17 of 29 196 (0.695, 0.699, 0.675, 0.701, and 0.614, respectively) for the construction cohort (Fig. 6B).Given the high annual death rates in pancreatic cancer, we compared the 1-year AUC values of the three-gene model to the clinical model for 1-year survival in the construction group.The three-gene model proved to be a superior model for determining patients' prognosis compared to commonly used clinical information markers, namely pathological tumor stage (pT), pathological nodal stage (pN), pathological stage, and resection margin (Fig. 6C). Figure 6D displays the 1-, 2-, 3-, 4-, and 5-year AUC values (0.703, 0.626, 0.612, 0.598, and 0.570) of the three-gene model for the validation group (E-MTAB-6134).The 1-year time-dependent AUC plots reveal that the threegene model also emerged as the best clinical marker for determining patients' prognosis in the validation groups (Fig. 6E).Decision curve analyses showed that the use of the three-gene predictive model was beneficial for patients in the construction group at 12 months after pancreatectomy (Fig. 6F).Furthermore, the three-gene model and the combined model (three-gene model high vs. low plus resection margin R1-2 vs. R0) showed a higher net benefit across all thresholds compared to other models for predicting OS of patients from the validation group at 12 months (Supplementary Fig. 6F).
To identify a reliable prognostic marker from the threegene purity-independent model for further investigations, we employed the Blasso algorithm (Zeng et al. 2020).This algorithm involves frequency-based screening of the stability of prognostic variables based on LASSO.We presented the top 20 variables selected in the TCGA-PAAD cohort with 1000 repetitions.We observed that ITGA3 (ENSG00000005884) was the only gene among the top 20 variables in the selection results of the three-gene purity-independent model (Fig. 6G).ITGA3 was chosen for further investigation because of its stability in the selection process (Fig. 6G) and its significant role in pancreatic cancer, as highlighted by various studies (Liu et al. 2021;Rashid et al. 2020;Wenric et al. 2022).Single-cell analyses of ITGA3 using the dataset provided by Peng et al. (2019) demonstrated its high expression in malignant cells (Fig. 6H), indicating its potential biological implications.

Implications of ITGA3 in genomics and association with treatment resistance
The mutation landscape differs between the high-ITGA3 and low-ITGA3 groups.The high-ITGA3 group exhibited higher rates of driver gene alterations, including KRAS, TP53, and CDKN2A.GNAS mutations, commonly found in intraductal papillary mucinous neoplasms, were more frequent in the low-ITGA3 group (Fig. 7A).We assessed the ability of ITGA3 participating in regulating the immune response.Although, there were no differences in Microsatellite Analysis for Normal-Tumor InStability (MANTIS) scores (regarding microsatellite instability, Fig. 7B) (Kautto et al. 2017).The high-ITGA3 group exhibited higher neoantigen loads (P < 0.05, Fig. 7C), cell components (CC) in T cell receptor complex and immunoglobulin complex (Fig. 7D), molecular function (MF) in antigen binding (Fig. 7E), and tumor mutational burden (P < 0.01, Fig. 7F) in high-ITGA3 samples.ITGA3 expression also differed among immune subtypes (P = 1.73 × 10 −4 , Fig. 7H).Most patients with the Moffit-basal subtype (Moffitt et al. 2015) exhibited high ITGA3 expression (P < 0.001).Patients with the Collisson quasi mesenchymal (QM-PDA) (Collisson et al. 2011) and Bailey squamous (Bailey et al. 2016) subtypes also exhibited high ITGA3 expression because of the overlapping classification of molecular subtypes (Fig. 7H).The correlation plots between ITGA3 and receptors, chemokines, tumorinfiltrating lymphocytes, major histocompatibility complex (MHC) molecules, immunostimulators, and immune inhibitors are shown in Fig. 7I.We observed that ITGA3 was positively correlated with classical and nonclassical MHC-I molecules and negatively correlated with MHC-II molecules and immune cell receptors.
Figure 8A demonstrates that ITGA3 is a risk factor for OS, PFS, DSS, and DFI.The expression of ITGA3 was strongly correlated with KRT19 (a marker for malignant ductal cells, correlation coefficient 0.781, P < 0.001), MKI67 (a marker for cell proliferation, correlation coefficient 0.433, P < 0.001), and TWIST (a marker for EMT, correlation coefficient 0.3, P < 0.001).Patients with proliferative cancer cells and EMT features exhibited ITGA3 expression independent of GATA6 (a marker for the basal-classical subtype, with a correlation P > 0.05).ITGA3 was a better prognostic marker than clinical factors (stage, gender, and age) in Fig. 6 Prognosis-related purity-independent three-gene network visualization.A Purity-independent regulatory network constructed using 3 prognosis-related mRNAs, 15 miRNAs, and 3 lncRNAs.B Time-dependent receiver operating characteristic (ROC) curves of the construction group (n = 540) for overall survival prediction using the three-gene model.C One-year time-dependent ROC curves for overall survival, showing higher area under the curve (AUC) values for the three-gene model compared to other clinic al factors in the validation groups (n = 288).D Time-dependent ROC curves of the validation group (n = 288) for overall survival prediction using the three-gene model.E One-year time-dependent ROC curves for overall survival, demonstrating higher AUC values for the three-gene model compared to other clinical factors in the validation groups (n = 288).F Decision curve analysis of the three-gene model, clinical models, and the combined model at 12 months post-pancreatectomy in the construction group (n = 540).G Frequency count of variables selected using the LASSO algorithm, repeated 1000 times in The Cancer Genome Atlas (TCGA) using BLASSO (n = 160) both univariate and multivariate Cox regression models (Fig. 8C) in TCGA dataset plotted using CAMOIP.Using our own processing protocols, we observed that ITGA3 (hazard ratio = 1.32, in multivariate analysis) was as an independent hazardous marker along with surgical margin (hazard ratio = 1.86, in multivariate analysis; Supplementary Fig. 7A) in TCGA-PAAD.GSEA results reveal that the keratinization and keratinocyte proliferation pathways differed between high-ITGA3 and low-ITGA3 groups, which explain the malignant phenotype in the high-ITGA3-high pancreatic cancer cohort with or without chemotherapy (Fig. 8D).No significant association between ITGA3 expression and tumor purity was found (Fig. 8E).Furthermore, the other two genes in the risk model did not exhibit a high correlation with tumor purity (Supplementary Fig. 7B).ITGA3 was not related to other compartment markers such as ACTA2 (a muscle and fibroblast marker), AMY1A (an acinar cell marker), or PTPRC (CD45, an immune cell marker), further supporting the idea that ITGA3 is not severely affected by environmental cells (Fig. 8E).Higher ITGA3 expression was correlated with poor prognosis in the Parikh et al. immunotherapy cohort, which includes pancreatic and colorectal cancer (Fig. 8F).ITGA3 was positively correlated with immune exhaustion markers CD274 (PD-L1, P < 0.05) and HAVCR2 (Tim-3, P < 0.05).Regarding the GSEA of immune pathways, downregulation of immune response, especially in adaptive immune response, was found to be more common in high-ITGA3 patients (Fig. 8H).The top 10 KEGG pathways using GSEA, GO-BP terms using GSEA, and Reactome terms using ORA reveal enrichment of various cancer and immune suppression pathways in high-ITGA3 pancreatic cancer samples (P < 0.05, Fig. 8I).
GSVA results show that high-ITGA3 patients exhibited a higher altered fraction, aneuploidy score, homologous recombination defects, proliferation, intratumor heterogeneity, wound healing, and number of segments.However, patients with high ITGA3 expression exhibited lower lymphocyte infiltration, TCR and BCR Shannon, TCR and BCR richness, macrophage regulation, and lymphocyte infiltration signature scores (Fig. 9A).We plotted current benchmark deconvolution methods in Fig. 9B-F and observed contradictory results because of different algorithms (MCPcounter, Fig. 9B; quanTIseq, Fig. 9C; CIBER-SORT, Fig. 9D; EPIC, Fig. 9E; IPS, Fig. 9F).Interestingly, we observed that CD8 + T cells were consistently downregulated across the methods in patients with high ITGA3 expression, suggesting that high ITGA3 expression is associated with an immune suppressive environment and a poorer response to immunotherapy.

Validation of ITGA3 expression in vivo and in vitro
In the present study, we analyzed ITGA3 in patient genomics data and validated its expression in cell lines and immunohistochemistry imaging, which are widely used by laboratory researchers and pathologists.ITGA3 expression was found to vary among pancreatic cancer cell lines from different molecular subtypes (Fig. 10A) but was relatively higher than in pan-organ normal fibroblasts based on data from the Cancer Dependency Map Project (DepMap) (Yu et al. 2019).Figure 10B depicted ITGA3 expression in normal pancreatic samples as well as low-, median-, and high-ITGA3 pancreatic cancer samples using the Human Protein Atlas (HPA, https:// www.prote inatl as.org/).While nonepithelial cells showed low but relatively constant ITGA3 expression, the greatest expression variances were observed in cancer cells (Fig. 10B).We also generated RNA-seq data from KRAS wild-type (BXPC3) and mutant (PANC1, MIA-PACA-2, and ASPC1) cancer cell lines, as well as the HPDE6C7 (an immortalized normal ductal cell line) and an immortalized primary CAF cell line (Fig. 10C, each cell line sequenced in triplicate).Comparing with single-cell (Fig. 1B) and single-nucleus (Fig. 6H) RNA-sequencing results, we observed notable ITGA3 expression in nonneoplastic cells in vitro (Fig. 10C).
The Moffit basal subtype cell line BXPC3 displayed the highest ITGA3 expression, supporting our inference (all P < 0.001).However, PANC1 exhibited a lower ITGA3 mRNA expression than the normal epithelial cell line HPDE6C7 (P = 4.3 × 10 −4 ) and CAFs (P = 7.8 × 10 −4 ).To address this seemingly contradictory result, we established an orthotopic transplantation model by injecting 5 × 10 5 cultured KPC cells into four C57BL/6J mice to investigate Itga3 expression in vivo.Subsequently, we performed single-cell sequencing of pancreatic cancer 28 days after orthotopic transplantation.We categorized the cell types in mice into epithelial, immune, and stromal compartments.Violin plots illustrating the expression of Itga3 (mouse gene version, Fig. 10D) using the Wilcoxon test revealed higher median expression in epithelial cells (including malignant cells and a few normal ductal, endocrine, and acinar cells) than in immune cells (mainly T, B, NK, and myeloid cells, P < 0.0001) and stromal cells (fibroblasts and endothelial cells, P < 0.0001).Our animal model results and human results are consistent with each other.We found that normal cells also express ITGA3; however, single-cell sequencing results in both humans and mice demonstrated that ITGA3 is primarily expressed in the epithelial compartment in vivo.Only a small number of epithelial cells exhibit weaker ITGA3 expression compared to immune and stromal compartments (Fig. 10D).

Discussion
Herein, the integrin superfamily genes were comprehensively analyzed.Using a purity-independent gene regulatory network, we demonstrated the importance of ITGA3 in pancreatic cancer as a prognostic marker for molecular transformation into the basal subtype and for patients who receive either chemotherapy or immunotherapy.
Studies have focused on the involvement of integrins in tumor stemness, proliferation, metastasis, and other malignant phenotypes (Cooper and Giancotti 2019;Dustin 2019;Weiss et al. 2022).Integrative analysis of a pan-cancer TCGA cohort, mostly comprising patients who received chemotherapy, revealed that the prognostic significance of distinct integrins varies across tumor types.The integrin/ FAK cascade controls YAP phosphorylation and nuclear localization, resulting in the activation of the YAP/TAZ signaling pathway, which controls stem cell expansion in normal and cancer tissues (Cooper and Giancotti 2019;Hu et al. 2017).The integrins α3β1 and α6β1 on the tumor surface mediate the transendothelial migration of circulating tumor cells.
The significance of integrins in immunotherapy has been reported (Cooper and Giancotti 2019;Ruan et al. 2020;Shaim et al. 2021).The integrins on immune cells play critical roles in adhesion, homing, and trafficking and can be used to specifically target a particular cell, thereby increasing or inhibiting the immune response (Dustin 2019;Ran et al. 2022).For example, in benign diseases, intestinal antigen-presenting cells can promote regulatory T cell activity via ITGB8, resulting in immune tolerance in the intestine to alleviate inflammatory bowel disease (Akagbosu et al. 2022).Furthermore, integrin αV/β8 on the surface of pancreatic cells induces DNA methylation and downregulates the expression of glucose metabolism-and oxidative phosphorylation-related genes in M1-like macrophages, but not in M2-like macrophages, resulting in immune suppression via direct cell-cell interaction.Blocking of integrin αV/β8 can suppress tumor-induced DNA methylation of Nqo-1 and reprogramming of macrophage metabolism (Zhang et al. 2021).In tumor therapy, T cell cytotoxicity is Fig. 8 Association of ITGA3 expression with chemotherapy and immunotherapy resistance.A ITGA3 is a risk factor for overall survival (OS, n = 178), progression-free survival (PFS, n = 178), disease-specific survival (DSS, n = 158), and disease-free interval (DFI, n = 58) (P values determined using log-rank tests).B Correlation of ITGA3 with KRT19 (malignant ductal cell marker, Spearman's correlation coefficient: 0.781, P < 0.001), MKI67 (cell proliferation marker, Spearman's correlation coefficient: 0.433, P < 0.001), and TWIST (EMT marker, Spearman's correlation coefficient: 0.3, P < 0.001).ITGA3 expression is independent of GATA6 expression (pancreatic cancer basal-classical subtype marker, Spearman correlation, P > 0.05).n = 179 C ITGA3 is a superior univariable and multivariable hazard factor compared to clinical markers (stage, age, and sex) based on CAMOIP.D Gene set enrichment analysis (GSEA) of Gene Ontology (GO) terms for biological processes in the ITGA3-high and ITGA3-low expression groups.Adjusted P values were calculated using the one-sided Fisher's exact test and corrected for multiple hypotheses using the Benjamini-Hochberg false discovery rate (FDR) of 5% (n = 179).E Spearman's correlation analysis of ITGA3 expression and cell component markers (ACTA2, AMY1A, and PTPRC), indicating the purity independence of ITGA3 expression in pancreatic cancer (P > 0.05, n = 179).F ITGA3 was associated with immunotherapy resistance for OS (Parikh et al., 2021, n = 23).P values determined using log-rank tests (P < 0.05).G Positive correlation between ITGA3 expression and the immunotherapy genes CD274 (PD-L1) and HAVCR2 (Spearman's correlation, P < 0.05).n = 179 H GSEA of GO terms for biological processes indicating the suppression of adaptive immunity in patients with high ITGA3 expression.n = 179 I GSEA of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and GO terms for biological processes, as well as Reactome analysis in patients with high ITGA3 expression (n = 179).Adjusted P values were calculated using the one-sided Fisher's exact test and corrected for multiple hypotheses using FDR enhanced by targeting the myeloid cell surface marker CD11b (ITGAM), resulting in the repolarization of macrophages, improvements in the responses of dendritic cells, and a decrease in myeloid-derived suppressor cells in pancreatic cancer (Bagati et al. 2021).We expanded the application of the integrin ITGA3, a basal subtype transforming gene associated with lesser degree of infiltration of CD8 + T cells and higher expression of PD-L1, to both chemotherapy and immunotherapy resistance.These lines of evidence confirm that using integrins to modulate immune responses may be a potential therapeutic target for treating pancreatic cancer.
Another aspect of the present study was the development of a purity-independent prognostic model, identification of the significance of ITGA3 as a single gene marker.Purity Independent Subtyping of Tumors (Pur-IST) and some studies have comprehensively explained risk stratification and treatment response from every aspect (Rashid et al. 2020;Shi et al. 2022;Tang et al. 2022).Compared to other models, our model possesses the advantage of being both simple and independent of tumor purity.Notably, the benchmark study PurIST also included ITGA3 in the 14-gene purity independent model to determine the basal subtype in pancreatic cancer (Rashid et al. 2020).Sample purity is an important issue hampering the identification of molecular subtypes in cancer.The application of spatial transcriptomics, bulk RNA sequencing, and bulk proteomics as well as other novel omics technologies for clinical sample analysis is currently affected by sample location, with different tumor heterogeneity and purity (Cao et al. 2021;Cui Zhou et al. 2022;Hwang et al. 2022).To construct the purity-independent regulatory model in pancreatic cancer, we adopted the ceRNA mechanism, which forms a link between the functions of protein-coding mRNAs and noncoding RNAs (such as lncRNAs, circRNAs, and miRNAs, including NamiRNAs) (Bailey et al. 2016;Liang et al. 2019;Yan et al. 2021) and has been considered a post-transcriptional regulation mechanism that is widely involved in the pathological processes of various cancers, including breast, colon, lung, gastric, and pancreatic cancers (Abdollahzadeh et al. 2019;X Liu et al. 2023b;Qi et al. 2015;Shuwen et al. 2018;Yan et al. 2021).By preserving the correlation between RNA interaction pairs under different tumor purities, we constructed the three-gene purity-independent network.For cancer research, it is important to select a feasible single-gene model (Ortega et al. 2022).Although the 14-gene model from the PurIST study and GATA6 as a single-gene model have been selected as benchmark classifiers for pancreatic cancer, we recommend ITGA3 as a promising supplementary marker for determining the immunotherapy or chemotherapy refractory status of patients with pancreatic cancer, regardless of the sampling location.We observed that HPDE6C7 cells and CAFs in pancreatic cancer also express ITGA3 to a certain extent in vitro.Similar results were observed in DepMap-based cell line analysis.Nonmalignant cells expressing ITGA3 contribute to steady ITGA3 expression, suggesting that ITGA3 serves as a marker for variations among molecular subtypes rather than for variations among tumor purities.However, we also observed a discrepancy in the expression of ITGA3 between laboratory cell lines and clinical samples, that the low expression of ITGA3 in some pancreatic cancer cells compared to non-malignant cells in vitro.However, using singlecell RNA-seq in humans and mice, and single-nucleus RNA-seq in humans, we observed that average ITGA3 expression was significantly higher in epithelial cells than that in immune and stroma cells.This situation dominates the majority in vivo.As cell lines represent only a fraction of the heterogeneous cancer cells found in vivo, and their gene expression in vitro may not align with in vivo observations.This highlights the issue of inadequate representativeness when relying solely on tumor cell lines for research purposes.Collaboration between clinicians and researchers are needed to bridge the gap between laboratory and clinical results (Sailer et al. 2023).The use of single-cell RNA sequencing has greatly expanded our understanding of precision therapy; the results of single-cell analysis can be combined to select the most efficient therapeutic targets (Yang et al. 2023;Yu et al. 2023).For example, a treatment regimen aimed at targeting the tumor stem cell marker ITGA4 may result in undesired effects owing to the high expression of ITGA4 in immune cells other than tumor cells, as illustrated in Fig. 1A (Ruan et al. 2020).In contrast, ITGA3 is highly localized to the tumor epithelium within tumor tissues, with weak expression in surrounding cells.
A limitation of the present study is that we analyzed public databases because we could not collect freely accessible data from single-cell and bulk sequencing data (Cui Zhou et al. 2022).To address this, we employed a combination of bulk, single-cell, and single-nucleus sequencing datasets, laboratory-collected cell lines, and mouse models to support the reliability of our conclusions.Another limitation is the inability to develop a highly generalized model from a single dataset with a small sample without combining the datasets, possibly introducing bias in the combined datasets used for model construction, which may introduce bias.However, we independently assessed the prognostic values of the three-gene model in each database.Finally, although we identified the advantages of ITGA3 in predicting chemotherapy and immunotherapy resistance, additional prospective clinical trials are warranted for further investigation (Dean et al. 2022;Principe et al. 2021;Slack et al. 2022;Zeltz et al. 2020).

Conclusions
Comprehensive multi-omics analyses showed dysregulation of integrin superfamily genes in pancreatic cancer and other cancers.Among all integrins, ITGA3 emerged as a promising single-gene biomarker with clinical relevance, specifically associated with resistance to chemotherapy and immunotherapy in pancreatic cancer owing to its purityindependent characteristics.
et al. 2019) and a single-nucleus RNA-seq study byHwang et al. (2022) were used together with annotated data from original articles.

Fig. 2
Fig.2Association between survival and prognosis and the enrichment pathways of integrin superfamily genes.A Cox proportional hazard regression modeling showing the association between overall survival (OS, n = 11505), progression-free survival (PFS, n = 11505), disease-specific survival (DSS, n = 10947), and disease-free interval (DFI, n = 5790) with integrin superfamily gene expression.B Bar plot displaying the top 10 Gene Ontology terms for biological processes associated with all integrins (n = 30).C Grid diagram illustrating the cancer pathways of the integrin superfamily genes.Red and blue shades represent the activated and inhibited pathways, respectively (n = 123).D Bubble plot illustrating the Spearman's correlation between integrin superfamily gene expression and immune infiltration in in pancreatic ductal adenocarcinoma (PAAD, n = 183), adjusted using false discovery rate (FDR).Higher color intensity reflects a stronger correlation, while bubble size represents FDR significance.Bubbles with black borders have FDRs < 0.05.Abbreviations for the cancer project are available at https:// gdc.cancer.gov/ resou rces-tcga-users/ tcga-code-tables/ tcga-study-abbre viati ons.Abbreviations and definitions for infiltrated immune cells were provided by the Immune Cell Abundance Identifier (ImmuCellAI) ◂

Fig. 3
Fig.3Mutation and alteration landscapes of integrin superfamily genes across multiple cancer types.A Deleterious mutations in integrin superfamily genes.The numbers within the boxes represent the number of deleterious gene mutations in different cancer types.An increase in color intensity indicates an increase in the percentage of deleterious mutations (n = 10234).B Pie charts showing the proportions of copy number alterations across cancer types (n = 11461).The red rectangle highlights information related to copy number variations in pancreatic ductal adenocarcinoma (PAAD, n = 184).CNV copy number variation, SNV single-nucleotide variation.The abbreviations for the cancer projects are available on the website: https:// gdc.cancer.gov/ resou rces-tcga-users/ tcga-code-tables/ tcga-study-abbre viati ons ◂ Fig.4DNA methylation of integrins in human cancers.A Bubble plot comparing methylation levels (β values) of integrin superfamily genes in normal and tumor tissues across multiple cancer types from The Cancer Genome Atlas (TCGA).Cancer types with < 10 normal samples were excluded (n = 6560).Bubble color represents methylation level differences between tumor and normal samples, ranging from blue to red.P values were estimated using the two-tailed unpaired t test and adjusted using the false discovery rate (FDR).Bubbles with black borders indicate FDRs < 0.05.Bubble size indicates FDR significance, while color intensity increases with correlation.B Spearman's correlation between DNA methylation levels (β values) and cell infiltration degree for each cancer type from TCGA, evaluated using ImmuCellAI (n = 183).P values were adjusted using FDR.Color intensity increases with correlation.Bubbles with black borders represent FDRs < 0.05.Bubble size indicates FDR significance.C Boxplots showing the distribution of DNA methylation levels (β values) of integrins among normal pancreatic tissues, TP53 wild-type, and TP53 mutant pancreatic cancer tissues (n = 192).P values were calculated using the two-tailed unpaired t-test.Boxes indicate the interquartile range, horizontal lines represent the median, and whiskers indicate values 1.5 times the interquartile range.Abbreviations for cancer projects are available at https:// gdc.cancer.gov/ resou rces-tcga-users/ tcga-code-tables/ tcga-study-abbre viati ons.Abbreviations and definitions for infiltrated cells were provided by Immune Cell Abundance Identifier (ImmuCellAI).Methy.diffmethylation difference, ns no significance ◂

Fig. 5
Fig.5Construction and validation of the purity-independent regulatory network.A Heatmap displaying eligible lncRNA-miRNA-mRNA pairs supported by GDCRNAtools.Only experiments confirming RNA relationships were included.RNA relationships with consistent evidence support across datasets with varying tumor purity levels were considered preserved in pancreatic cancer (n = 3).B Model coefficients (top axis) versus the logarithm of lambda (bottom axis) for LASSO analysis of the construction group (n = 540).The coefficients of three genes were non-zero at lambda.1se, which represents the lambda value generating the most regularized model with cross-validated errors within one standard error of the minimum.C Three-fold cross-validation C-indexes versus the logarithm of lambda (bottom axis) and the number of non-zero coefficients (top axis) in the construction group (n = 540).D Kaplan-Meier curve illustrating the high-and low-risk groups in the construction model using the three-gene model (n = 540, log-rank test, P < 0.0001).E Kaplan-Meier curve illustrating the high-and low-risk groups in the validation model using the three-gene model (n = 288, log-rank test, P = 0.0067).F Risk score (RS) panel of the construction model (n = 540).G Risk score (RS) panel of the validation model (n = 288) ◂ Fig.6Prognosis-related purity-independent three-gene network visualization.A Purity-independent regulatory network constructed using 3 prognosis-related mRNAs, 15 miRNAs, and 3 lncRNAs.B Time-dependent receiver operating characteristic (ROC) curves of the construction group (n = 540) for overall survival prediction using the three-gene model.C One-year time-dependent ROC curves for overall survival, showing higher area under the curve (AUC) values for the three-gene model compared to other clinic al factors in the validation groups (n = 288).D Time-dependent ROC curves of the validation group (n = 288) for overall survival prediction using the three-gene model.E One-year time-dependent ROC curves for overall survival, demonstrating higher AUC values for the three-gene model compared to other clinical factors in the validation groups (n = 288).F Decision curve analysis of the three-gene model, clinical models, and the combined model at 12 months post-pancreatectomy in the construction group (n = 540).G Frequency count of variables selected using the LASSO algorithm, repeated 1000 times in The Cancer Genome Atlas (TCGA) using BLASSO (n = 160).H Violin plot of ITGA3 expression constructed using the single-cell pancreatic cancer dataset of Peng et al. Cell types are arranged from left to right based on decreasing ITGA3 expression (n = 35) ◂

Fig. 7
Fig.7Association of ITGA3 with genomic and molecular characteristics of pancreatic cancer.A Mutation landscape of the ITGA3high and ITGA3-low expression groups.The Fisher's exact test was used to determine statistical differences between groups (n = 175).B Microsatellite Analysis for Normal-Tumor InStability (MANTIS) scores in the ITGA3-high and ITGA3-low groups (n = 175).C Neoantigen loads in the ITGA3-high and ITGA3-low groups (n = 160).D Gene Ontology (GO) analysis of cell component (CC) in the ITGA3high and ITGA3-low groups (n = 179).E GO analysis of molecular function (MF) in the ITGA3-high and ITGA3-low groups.Adjusted P values were calculated using the one-sided Fisher's exact test and corrected for multiple hypotheses using the Benjamini-Hochberg false discovery rate (FDR) of 5% (n = 179).F Tumor mutation burden (TMB) in the ITGA3-high and ITGA3-low groups (n = 167).G Differential expression of ITGA3 among immune subtypes (n = 151).H High expression of ITGA3 in common molecular subtypes, including Moffit basal subtypes, Collisson quasimesenchymal subtypes (QM-PDA), and Bailey squamous subtypes.Differences between groups were assessed using the chi-square test (n = 150).I Spearman's correlation analysis between ITGA3 and receptors, chemokines, tumorinfiltrating lymphocytes, MHC molecules, immune stimulators, and inhibitors in pan-cancer analysis.Blue squares indicate negative correlations, while red squares indicate positive correlations; darker colors indicate stronger correlations (n = 179).B, C, F, and G: Data represented as boxplots, and differences in immune-related scores or immune cells were compared using the Mann-Whitney U test.Boxes represent interquartile ranges, horizontal lines within each box represent the median, and whiskers represent values 1.5 times the interquartile range.Black dots indicate outliers.Statistical significance is indicated as follows: *P < 0.05, **P < 0.01, ***P < 0.001, and ****P < 0.0001 ◂

Fig. 9 001 ◂
Fig.9Effect of ITGA3 expression on immune response pathways and cell composition in pancreatic cancer tissue.A Gene set variation analysis (GSVA) comparing immune response pathways in ITGA3high and ITGA3-low groups (n = 176).B-E Deconvolution results of cell components using bulk RNA sequencing with The Cancer Genome Atlas-pancreatic ductal adenocarcinoma (TCGA-PAAD) dataset (n = 177).Cell composition analyzed using MCPcounter (B), quanTIseq (C), CIBERSORT (D), EPIC (E), and IPS (F) in samples with high and low ITGA3 expression.Boxplots illustrate the data, and the Mann-Whitney U test was used to compare immune-related scores or immune cells.Boxes represent the interquartile range, horizontal lines represent the median, and whiskers indicate values 1.5 times the interquartile range.Outliers are depicted as black dots.Note: ns no significance; *P < 0.05, **P < 0.01, and ***P < 0.001 ◂

Fig. 10
Fig. 10 Validation of ITGA3 expression in pancreatic cancer.A Expression of ITGA3 in pancreatic cancer cell lines and normal fibroblasts from the Cancer Dependency Map Project (DepMap, n = 88).B Immunohistochemistry of ITGA3 in pancreatic tissue and cancer samples from the Human Protein Atlas (HPA, n = 4).C ITGA3 expression in pancreatic cells (PANC-1, BxPC3, HPDE6C7, ASPC1, MIA-PACA-2, and cancer-associated fibroblasts[CAFs]) in our laboratory; each sequencing analysis was performed three times.Gene expression differences among cell lines were compared using the t-test, n = 3. D Validation of Itga3 expression using single-cell RNA sequencing in a mouse model developed in our laboratory.The violin plot illustrates significantly higher expression of Itga3 in epithelial cells than in immune cells (P < 0.0001) and stroma cells (P < 0.0001).Gene expression differences between epithelial cells and other cells were compared using the Mann-Whitney U test, n = 4

Table 1
Patient's baseline characteristics used to construct the purity-independent RNA network

Table 2
Number of samples and platforms used in purity-independent regulatory network prediction

Table 3
Clinical information of the datasets included in purity-independent three-gene model construction and validation TCGA The Cancer Genome Atlas, ICGC International Cancer Genome Consortium, NA missing information Dot plot illustrating integrin superfamily gene expression across multiple cell types from the single-cell RNA sequencing dataset of treatment naïve PAAD ofHwang et al (n  = 18).Each circle's size corresponds to the percentage of cells exhibiting gene expression.Gene expression refers to the average expression values of genes rescaled by all cell types.Statistical significance is denoted as follows: *P < 0.05, **P < 0.01, and ***P < 0.001.Cancer project abbreviations are available on the website: https:// gdc.cancer.gov/ resou rces-tcga-users/ tcga-code-tables/ tcga-study-abbre viati ons