Progress is being made in the prevention and treatment of chronic obstructive pulmonary disease (COPD), but it is still unsatisfactory. With the development of genetic technology, validated genetic information can better explain COPD. The study utilized scRNA-seq and Mendelian randomization analysis of eQTLs to identify crucial genes and potential mechanistic pathways underlying COPD pathogenesis. Single-cell sequencing data were used to identify marker genes for immune cells in the COPD process. Data on eQTLs for immune cell marker genes were obtained from the eQTLGen consortium. To estimate the causal effect of marker genes on COPD, we selected an independent cohort (ukb-b-16751) derived from the UK Biobank database for two-sample Mendelian randomization analysis. Subsequently, we performed immune infiltration analysis, gene set enrichment analysis (GSEA), and co-expression network analysis on the key genes. The 154 immune cell-associated marker genes identified were mainly involved in pathways such as vacuolar cleavage, positive regulation of immune response and regulation of cell activation. Mendelian randomization analysis screened four pairs of marker genes (GZMH, COTL1, CSTA and CD14) were causally associated with COPD. These four key genes were significantly associated with immune cells. In addition, we have identified potential transcription factors associated with these key genes using the Cistrome database, thus contributing to a deeper understanding of the regulatory network of these gene expressions. This eQTLs Mendelian randomization study identified four key genes (GZMH, COTL1, CSTA, and CD14) causally associated with COPD, providing new insights for prevention and treatment of COPD.