IntroductionCurrently, programmed cell death-1 (PD-1)-targeted treatment is ineffective for a sizable minority of patients, and drug resistance still cannot be overcome.MethodsTo explore the mechanisms of immunotherapy and identify new therapeutic opportunities in lung adenocarcinoma (LUAD), data from patients who did and did not respond to the anti-PD-1 treatment were evaluated using single-cell RNA sequencing, and bulk RNA sequencing were collected.ResultsWe investigated the gene expression that respond or not respond to immunotherapy in diverse cell types and revealed transcriptional characteristics at the single-cell level. To ultimately explore the molecular response or resistance to anti-PD-1 therapy, cell-cell interactions were carried out to identify the different LRIs (ligand-receptor interactions) between untreated patients vs. no-responders, untreated patients vs. responders, and responders vs. non-responders. Next, two molecular subgroups were proposed based on 73 LRI genes, and subtype 1 had a poor survival status and was likely to be the immunosuppressive tumor subtype. Furthermore, based on the LASSO Cox regression analysis results, we found that TNFSF13, AXL, KLRK1, FAS, PROS1, and CDH1 can be distinct prognostic biomarkers, immune infiltration levels, and responses to immunotherapy in LUAD.DiscussionAltogether, the effects of immunotherapy were connected to LRIs scores, indicating that potential medications targeting these LRIs could contribute to the clinical benefit of immunotherapy. Our integrative omics analysis revealed the mechanisms underlying the anti-PD-1 therapy response and offered abundant clues for potential strategies to improve precise diagnosis and immunotherapy.