BackgroundPancreatic adenocarcinoma (PAAD) is a leading cause of malignancy-related deaths worldwide, and the efficacy of immunotherapy on PAAD is limited. Studies report that long non-coding RNAs (lncRNAs) play an important role in modulating genomic instability and immunotherapy. However, the identification of genome instability-related lncRNAs and their clinical significance has not been investigated in PAAD.MethodsThe current study developed a computational framework for mutation hypothesis based on lncRNA expression profile and somatic mutation spectrum in pancreatic adenocarcinoma genome. We explored the potential of GInLncRNAs(genome instability-related lncRNAs) through co-expression analysis and function enrichment analysis. We further analyzed GInLncRNAs by Cox regression and used the results to construct a prognostic lncRNA signature. Finally, we analyzed the relationship between GILncSig (genomic instability derived 3-lncRNA signature) and immunotherapy.ResultsA GILncSig was developed using bioinformatics analyses. It could divide patients into high-risk and low-risk groups, and there was a significant difference in OS between the two groups. In addition, GILncSig was associated with genome mutation rate in pancreatic adenocarcinoma, indicating its potential value as a marker for genomic instability. The GILncSig accurately grouped wild type patients of KRAS into two risk groups. The prognosis of the low-risk group was significantly improved. GILncSig was significantly correlated with the level of immune cell infiltration and immune checkpoint.ConclusionsIn summary, the current study provides a basis for further studies on the role of lncRNA in genomic instability and immunotherapy. The study provides a novel method for identification of cancer biomarkers related to genomic instability and immunotherapy.Supplementary InformationThe online version contains supplementary material available at 10.1186/s12885-023-10831-4.