Depression and coronary heart disease (CHD) have common risk mechanisms. Common single nucleotide polymorphisms (SNPs) may be associated with the risk of depression combined with coronary heart disease. This study was designed according to the PRISMA-P guidelines. We will include case-control studies and cohort studies investigating the relationship between gene SNPs and depression and coronary heart disease comorbidities. The Newcastle-Ottawa Scale (NOS) will be used to assess the risk of bias. When measuring dichotomous outcomes, we will use the odds ratio (OR) and 95% confidence interval (95%CIs) in a case-control study. Five genetic models (allele model, homozygous model, co-dominant model, dominant model, and recessive model) will be evaluated for each included study. Subgroup analysis by ethnicity will be performed. If necessary, post hoc analysis will be made according to different types. A total of 13 studies were included in this study, and the types of genes included are FKBP5 and SGK1 genes that act on glucocorticoid; miR-146a, IL-4-589, IL-6-174, TNF-α-308, CRP-717 genes that act on inflammatory mechanisms; eNOS genes from endothelial cells; HSP70 genes that act on the autoimmune response; ACE2 and MAS1 genes that act to mediate Ang(1-7) in the RAS system; 5-HTTLPR gene responsible for the transport of serotonin 5-HT and neurotrophic factor BDNF gene. There were three studies on 5-HTTLPR and BDNF genes, respectively, while there was only one study targeting FKBP5, SGK1, miR-146a, IL-4-589, IL-6-174, TNF-alpha-308, CRP-717, eNOS, HSP70, ACE2, and MAS1 genes. We did not perform a meta-analysis for genes reported in a single study, and meta-analysis was performed separately for studies exploring the 5-HTTLPR and BDNF genes. The results showed that for the 5-HTTLPR gene, there was a statistically significant association between 5-HTTLPR gene polymorphisms and depression in combination with coronary diseases (CHD-D) under the co-dominant model (LS vs LL: OR 1.76, 95%CI 1.20-2.59; SS vs LL: OR 2.80, 95%CI 1.45 to 5.41), the dominant model (LS+SS vs LL: OR 2.06, 95%CI 1.44 to 2.96), and the homozygous model (SS vs LL: OR 2.80 95%CI 1.45 to 5.5.41) were statistically significant for CHD-D, demonstrating that polymorphisms in the 5-HTTLPR gene are associated with the development of CHD-D and that the S allele in the 5-HTTLPR gene is likely to be a risk factor for CHD-D. For the BDNF gene, there were no significant differences between one of the co-dominant gene models (AA vs GG: OR 6.63, 95%CI 1.44 to 30.64), the homozygous gene model (AA vs GG: OR 6.63,95% CI 1.44 to 30.64), the dominant gene model (GA+AA vs GG: OR4.29, 95%CI 1.05 to 17.45), recessive gene model (AA vs GG+GA: OR 2.71, 95%CI 1.16 to 6.31), and allele model (A vs G: OR 2.59, 95%CI 1.18 to 5.67) were statistically significant for CHD-D, demonstrating that BDNFrs6265 gene polymorphisms are associated with the CHD-D development and that the A allele in the BDNFrs6265 gene is likely to be a risk factor for CHD-D. We analyzed the allele frequencies of SNPs reported in a single study and found that the SNPs in the microRNA146a gene rs2910164, the SNPs in the ACE2 gene rs2285666 and the SNPs in the SGK1 gene rs1743963 and rs1763509 were risk factors for the development of CHD-D. We performed a subgroup analysis of three studies involving the BDNFrs6265 gene. The results showed that European populations were more at risk of developing CHD-D than Asian populations in both dominant model (GA+AA vs GG: OR 10.47, 95%CI 3.53 to 31.08) and co-dominant model (GA vs GG: OR 6.40, 95%CI 1.98 to 20.73), with statistically significant differences. In contrast, the studies involving the 5-HTTLPR gene were all Asian populations, so subgroup analyses were not performed. We performed sensitivity analyses of studies exploring the 5-HTTLPR and BDNF rs6265 genes. The results showed that the results of the allele model, the dominant model, the recessive model, the homozygous model and the co-dominant model for both 5-HTTLPR and BDNF rs6265 genes were stable. Due to the limited number of studies of the 5-HTTLPR and BDNF genes, it was not possible to determine the symmetry of the funnel plot using Begg's funnel plot and Egger's test. Therefore, we did not assess publication bias. SNPs of the microRNA146a gene at rs2910164, the ACE2 gene at the rs2285666 and the SGK1 gene at rs1743963 and rs1763509, and the SNPs at the 5-HTTLPR and BDNF gene loci are associated with the onset of comorbid depression in coronary heart disease. We recommend that future research focus on studying SNPs' impact on comorbid depression in coronary heart disease, specifically targeting the 5-HTTLPR and BDNF gene at rs6265. https://www.crd.york.ac.uk/prospero/, identifier CRD42021229371.