Missense Variant rs28362680 in BTNL2 Reduces Risk of Coronary Heart Disease

Introduction

Coronary heart disease (CHD), as a major cardiovascular disease, is the leading cause of death in the global population. The pathological basis of coronary heart disease is atherosclerosis (AS).¹ CHD is a complex disease that involves multiple mechanisms, multiple cell types, and is affected by multiple environmental factors, including smoking,² drinking,³ poor eating habits,⁴ lifestyle, etc.⁵ These life-related factors further promote the occurrence of traditional risk factors for coronary heart disease, including diabetes,⁶ dyslipidemia⁷ and hypertension.^8,9 In addition, with the development of molecular Biology and the improvement and application of genetic testing technology, the influence of genetic factors on disease risk has been highly valued. More and more studies have proved that genetic factors play an important role in the occurrence and development of CHD.¹⁰ CHD has an important genetic basis and is considered equivalent to environmental factors. At present, some genes single nucleotide polymorphism (SNP) have been found to be associated with CHD susceptibility, such as XKR6,¹¹ interleukin-6 receptor,¹² T-cell immunoglobulin and mucin domain 4 gene,¹³ etc. However, it should be noted that the genetic mechanism of coronary heart disease is not fully clearly. The discovery of SNPs associated with the occurrence and development of diseases plays an important role in exploring the susceptibility and pathogenesis of diseases, making the prevention, diagnosis and treatment of diseases in the future full of hope.

BTNL2 belongs to the butyrophilin-like (BTNL) family. Studies have found that BTNL2 acts as a co-inhibitory molecule for T cell activation.¹⁴ In addition, the genetic polymorphism of BTNL2 is associated with the susceptibility of Kawasaki disease¹⁵ and dilated cardiomyopathy.¹⁶ Kawasaki disease is a child vasculitis with a higher risk of coronary heart disease in adulthood.¹⁷ In summary, BTNL2 is very likely to participate in the occurrence and development of coronary heart disease. But so far, there is no relevant report on the association between BTNL2 genetic polymorphism and coronary heart disease.

Therefore, this study intends to conduct a “case-control” study with the Han population in southern China as the research subject. The purpose is to explore the association between BTNL2 genetic polymorphism and CHD susceptibility. At the same time, we conducted a stratified analysis based on other potential risk factors for CHD (age, smoking, drinking, hypertension, diabetes, etc.) to assess the association between BTNL2 genetic polymorphisms and these influencing factors. This study will lay a foundation for the study on mechanism of BTNL2 in the occurrence and development of coronary heart disease. And it will provide new ideas for the individualized prevention and treatment of CHD.

Materials and Methods

Sample Source

We recruited 710 patients with coronary heart disease and 709 healthy individuals for this case-control study. The 1419 participants were recruited in the same hospital (People’s Hospital of Wanning). We continuously recruited 710 CHD patients in the outpatient or inpatient department of cardiovascular. These patients were all newly diagnosed. Recruitment criteria: patients met the clear and diagnostic criteria of the World Health Organization (WHO) on coronary heart disease. Patients with coronary angiography indications were examined by coronary angiography, screening patients with coronary angiography at least one lumen diameter less than 50% or more. Exclusion criteria: patients with cardiomyopathy, heart valve disease, congenital heart disease or renal failure were excluded. During the same period, 709 healthy individuals were recruited in the health examination center as the control group in this study. Control group recruitment criteria: participants without family history of premature coronary heart disease, history of ischemic cardiomyopathy or history of chronic heart failure, and participants without hypertension, diabetes, hyperlipidemia, hyperuricemia and obesity will be enrolled in control group. We obtained the demographic information and basic epidemiological information of the participants by consulting medical records or professional questionnaire survey. After obtaining the informed consent of all participants, we collected their blood samples for further study. This study has been approved by the ethics committee of the Northwest University and People’s Hospital of Wanning before the beginning.

Selection of SNPs

First, we found that the physical position of the BTNL2 was on the Chromosome 6: 32393963–32407128 through the e!GRCh37 database (http://asia.ensembl.org/Homo_sapiens/Info/Index). Then, we use the online converter window (VCF to PED: http://grch37.ensembl.org/Homo_sapiens/Tools/VcftoPed) to download the related files of BTNL2 genetic variants. For the population, we choose CHB and CHS. Finally, after setting specific conditions, six candidate BTNL2 SNPs (rs35624343, rs117896888, rs41441651, rs41417449, rs28362680, rs2076523) were selected using Haploview software. These special conditions are as follows: Tagger r² > 0.8, Min Genotype > 75%, MAF> 0.05 and HWE> 0.01.

DNA Extraction and Genotyping

In this study, we chose the kit (GoldMag Co. Ltd. Xi’an, China) to complete the extraction and purification of genomic DNA, and the specific experimental steps were carried out according to the instructions. According to the DNA sequence of BTNL2, we designed specific amplification primers and extension primers for candidate genetic loci (All primers can be found in Supplemental Table 1) by MassARRAY Assay Design software. Genotyping was performed by the MassARRAY ^® -IPLEX SNP genotyping technology.

At the same time, we randomly selected 5% DNA samples for repeated experiments. The repetition rate of the results should reach >99% to ensure the reliability and repeatability of the experimental results.

Data Analysis

We used online software (HaploReg v4.1: https://pubs.broadinstitute.org/mammals/haploreg/haploreg.php) to predict the potential function of candidate SNPs, and the information related to SNPs was searched through dbSNP online database (https://www.ncbi.nlm.nih.gov/snp/). The statistical analysis of this study was mainly completed by SPSS 22.0 software (SPSS Inc., Chicago, IL, USA). In this study, the association between CHD risk and candidate genetic polymorphism under multiple genetic models was completed by SNPStats online software (https://www.snpstats.net/start.htm?q=snpstats/start.htm). We mainly evaluated the impact of candidate SNPs on the risk of coronary heart disease through the Odds ratios (OR) and 95% confidence intervals (CI). In order to avoid the influence of confounding factors on the results, all the results were adjusted by the confounding factors, such as age, gender, smoking or drinking. We also used false-positive report probability (FPRP) analysis to detect whether all positive results are noteworthy at a prior probability level of 0.25 and FPRP threshold of 0.2. At the same time, the statistic power analysis will be completed. Haploview 4.2 software and SNPStats online software were used to perform the haplotype analysis of candidate SNPs and evaluation of linkage disequilibrium (LD). Finally, the interaction of candidate SNPs in IS risk was evaluated by multi-factor dimensionality reduction (MDR). The p < 0.05 indicated statistically significant.

Results

There were 1419 participants in the case group (710) and the control group (709) enrolled in this study. The average ages of the case and control groups were 62.58 ± 10.49 years and 62.23 ± 7.44 years (Table 1). The number of males and females in the case group were 481 (67.7%) and 229 (32.3%), and in the control group were 479 (67.6%) and 230 (32.4%), respectively. There was no statistical difference in gender and age between the participants in the two groups, which indicated that the two groups were matched in age and gender. The basic information of the participants can be found in Table 1.

Table 1 Characteristics of Patients with CHD and Healthy Individuals

Genotyping and Information About Candidate SNPs

Six BTNL2 candidate genetic loci (rs35624343 A/G, rs117896888 G/C, rs41441651 T/C, rs41417449 C/T, rs28362680 A/G, rs2076523 C/T) have been successfully genotyped. These genetic variants are all missense variants in BTNL2, and dbSNP online database shows that these missense variants may lead to amino acid changes (Table 2). After testing, the distribution of all candidate genetic loci is in accordance with Hardy-Weinberg equilibrium (HWE p > 5%). HaploReg online software was used to predict the potential functions of candidate genetic loci and found that the 6 SNPs may be regulated by a variety of factors. These factors include Promoter histone marks; Enhancer histone marks; DNAse; Proteins bound; Motifs changed; Selected eQTL hits, etc. (Table 2).

Table 2 The Basic Information and HWE About the Selected SNPs of BTNL2

Association Analysis Between BTNL2 SNPS and CHD Risk (Overall Analysis)

The association analysis between candidate SNPs and the risk of CHD showed that BTNL2-rs35624343, -rs117896888, -rs41441651, -rs41417449 and -rs28362680 were significantly associated with the susceptibility of CHD among participants (Table 3). Specifically, compared with the GG genotype, the GA genotype of BTNL2-rs35624343 was a protective factor for CHD (GA vs GG, OR = 0.78, 95% CI = 0.62–0.99, p = 0.040). BTNL2-rs35624343 can significantly reduce the CHD risk under dominant (GA-AA vs GG, OR = 0.78, 95% CI = 0.62–0.98, p = 0.036), overdominant (GG-AA vs GA, OR = 0.79, 95% CI = 0.62–0.99, p = 0.045), log-additive genetic models (OR = 0.80, 95% CI = 0.65–0.99, p = 0.043). -rs117896888 can significantly reduce the CHD risk under allele (G vs C, OR = 0.74, 95% CI = 0.55–0.98, p = 0.036) and log-additive genetic models (OR = 0.73, 95% CI = 0.55–0.98, p = 0.035). -rs41441651 can significantly reduce the CHD risk under allele (T vs C, OR = 0.73, 95% CI = 0.54–0.97, p = 0.028), dominant (CT-TT vs CC, OR = 0.72, 95% CI = 0.53–0.98, p = 0.036) and log-additive genetic models (OR = 0.72, 95% CI = 0.54–0.97, p = 0.027). -rs41417449 can significantly reduce the CHD risk under allele (C vs T, OR = 0.73, 95% CI = 0.55–0.97, p = 0.031), dominant (CT-CC vs TT, OR = 0.73, 95% CI = 0.54–0.99, p = 0.040) and log-additive genetic models (OR = 0.73, 95% CI = 0.54–0.97, p = 0.030). Compared with wild allele G, allele A of -rs28362680 was a protective factor for the risk of CHD (A vs G, OR = 0.68, 95% CI = 0.57–0.81, p < 0.0001). Compared with the GG genotype, the AA (OR = 0.40, 95% CI = 0.23–0.67, p = 0.001) and GA (OR = 0.68, 95% CI = 0.54–0.84, p < 0.0001) genotypes were protective factors for CHD risk. In addition, -rs28362680 significantly reduced the risk of CHD among participants under all genetic models (dominant: OR = 0.64, 95% CI = 0.52–0.79, p < 0.0001; recessive: OR = 0.47, 95% CI = 0.28–0.78, p = 0.003; overdominant: OR = 0.73, 95% CI = 0.59–0.90, p = 0.004; log-additive: OR = 0.66, 95% CI = 0.55–0.79, p < 0.0001).

Table 3 The Association Analysis Between Susceptibility of CHD and Single Nucleotide Polymorphism of BTNL2

We found no evidence that BTNL2-rs28362680 was associated with the CHD risk among the participants.

Association Analysis Between BTNL2 SNPS and CHD Risk (Stratified Analysis)

In order to further avoid the influence of other factors on the reliability of the association analysis results, we divided the participants according to different influencing factors and carried out a stratified analysis (age, gender, smoking, drinking, or CHD complicated with hypertension/diabetes).

Age: the results showed that BTNL2-rs28362680 was significantly associated with CHD risk in people ≤ 62 years old or > 62 years old (Supplemental Table 2). Among participants ≤ 62 years old, -rs28362680 significantly reduced the risk of CHD under all genetic models (Allele: OR = 0.66, 95% CI = 0.52–0.84, p = 0.001. Genotype (AA): OR = 0.66, 95% CI = 0.52–0.84, p = 0.001; Genotype (GA): OR = 0.67, 95% CI = 0.49–0.92, p = 0.012. Dominant: OR = 0.62, 95% CI = 0.46–0.84, p = 0.002. Recessive: OR = 0.37, 95% CI = 0.17–0.80, p = 0.008. Overdominant: OR = 0.73, 95% CI = 0.54–0.98, p = 0.037. Log-additive: OR = 0.63, 95% CI = 0.48–0.81, p < 0.0001). Among participants > 62 years old, -rs28362680 significantly reduced the risk of CHD under multiple genetic models (Allele: OR = 0.70, 95% CI = 0.55–0.90, p = 0.005. Genotype (AA): OR = 0.44, 95% CI = 0.20–0.93, p = 0.032; Genotype (GA): OR = 0.67, 95% CI = 0.48–0.93, p = 0.017. Dominant: OR = 0.66, 95% CI = 0.49–0.89, p = 0.007. Overdominant: OR = 0.71, 95% CI = 0.52–0.97, p = 0.033. Log-additive: OR = 0.70, 95% CI = 0.54–0.90, p = 0.006). In the age stratified analysis, the remaining 5 candidate SNPs were not associated with the CHD risk.

Gender: the results showed that BTNL2-rs28362680 was significantly associated with CHD susceptibility in both male and female participants (Supplemental Table 2). Among female participants, -rs28362680 significantly reduced the risk of CHD under all genetic models (Allele: OR = 0.57, 95% CI = 0.42–0.79, p = 0.001. Genotype (AA): OR = 0.31, 95% CI = 0.12–0.84, p = 0.020; Genotype (GA): OR = 0.56, 95% CI = 0.38–0.83, p = 0.004. Dominant: OR = 0.47, 95% CI = 0.32–0.70, p < 0.0001. Recessive: OR = 0.38, 95% CI = 0.14–1.01, p = 0.042. Overdominant: OR = 0.55, 95% CI = 0.36–0.82, p = 0.004. Log-additive: OR = 0.51, 95% CI = 0.36–0.73, p < 0.0001). Among male participants, -rs28362680 significantly reduced the risk of CHD under multiple genetic models (Allele: OR = 0.73, 95% CI = 0.59–0.90, p = 0.003. Genotype (AA): OR = 0.45, 95% CI = 0.24–0.85, p = 0.013; Genotype (GA): OR = 0.73, 95% CI = 0.56–0.96, p = 0.022. Dominant: OR = 0.69, 95% CI = 0.53–0.89, p = 0.004. Recessive: OR = 0.51, 95% CI = 0.27–0.95, p = 0.030. Log-additive: OR = 0.70, 95% CI = 0.56–0.87, p = 0.001).

In addition, allele A (A vs G: OR = 0.76, 95% CI = 0.59–0.97, p = 0.031) and genotype GA of BTNL2-rs35624343 (GA vs GG: OR = 0.74, 95% CI = 0.55–0.99, p = 0.040) were protective factors for CHD in male participants. BTNL2-rs35624343 also significantly reduced the risk of CHD under dominant (OR = 0.74, 95% CI = 0.56–0.98, p = 0.039) and log-additive genetic models (OR = 0.76, 95% CI = 0.58–0.98, p = 0.037).

Allele G (G vs C: OR = 0.71, 95% CI = 0.50–1.00, p = 0.049) and genotype CG of BTNL2-rs117896888 (CG vs CC: OR = 0.68, 95% CI = 0.47–0.98, p = 0.038) were protective factors for CHD in male participants. BTNL2-rs117896888 also significantly reduced the risk of CHD under dominant (OR = 0.67, 95% CI = 0.47–0.97, p = 0.032), overdominant (OR = 0.66, 95% CI = 0.46–0.96, p = 0.030) and log-additive genetic models (OR = 0.69, 95% CI = 0.49–0.99, p = 0.040).

Allele T (T vs C: OR = 0.69, 95% CI = 0.49–0.98, p = 0.035) and genotype CT of BTNL2-rs41441651 (CT vs CC: OR = 0.66, 95% CI = 0.46–0.95, p = 0.026) were protective factors for CHD in male participants. BTNL2-rs41441651 also significantly reduced the risk of CHD under dominant (OR = 0.65, 95% CI = 0.45–0.94, p = 0.022), overdominant (OR = 0.65, 95% CI = 0.45–0.94, p = 0.020) and log-additive genetic models (OR = 0.68, 95% CI = 0.48–0.96, p = 0.028).

Allele C (C vs T: OR = 0.70, 95% CI = 0.49–0.99, p = 0.040) and genotype CT of BTNL2-rs41417449 (CT vs TT: OR = 0.67, 95% CI = 0.46–0.96, p = 0.030) were protective factors for CHD in male participants. BTNL2-rs41417449 also significantly reduced the risk of CHD under dominant (OR = 0.66, 95% CI = 0.46–0.95, p = 0.026), overdominant (OR = 0.65, 95% CI = 0.45–0.95, p = 0.024) and log-additive genetic models (OR = 0.68, 95% CI = 0.48–0.97, p = 0.032).

We found no evidence that BTNL2-rs2076523 was associated with the CHD risk in age or gender stratified analysis.

Smoking: the results showed that BTNL2-rs28362680 was significantly associated with CHD susceptibility among smokers or non-smokers (Supplemental Table 3). Among smoking participants, -rs28362680 significantly reduced the risk of CHD under all genetic models (Allele: OR =0.65, 95% CI = 0.51–0.82, p < 0.0001. Genotype (AA): OR = 0.37, 95% CI = 0.18–0.80, p = 0.011; Genotype (GA): OR = 0.59, 95% CI = 0.43–0.81, p = 0.001. Dominant: OR = 0.56, 95% CI = 0.41–0.77, p < 0.0001. Recessive: OR = 0.47, 95% CI = 0.22–0.99, p = 0.040. Overdominant: OR = 0.64, 95% CI = 0.47–0.87, p = 0.005. Log-additive: OR = 0.60, 95% CI = 0.46–0.78, p < 0.0001). Among non-smoking participants, -rs28362680 significantly reduced the risk of CHD under multiple genetic models (Allele: OR = 0.71, 95% CI = 0.56–0.91, p = 0.007. Genotype (AA): OR = 0.40, 95% CI = 0.19–0.85, p = 0.018. Dominant: OR = 0.70, 95% CI = 0.52–0.95, p = 0.020. Recessive: OR = 0.46, 95% CI = 0.22–0.96, p = 0.033. Log-additive: OR = 0.70, 95% CI = 0.54–0.91, p = 0.006).

Compared with wild allele G, allele A of BTNL2-rs35624343 was significantly associated with reducing CHD risk among smokers (A vs G, OR = 0.70, 95% CI = 0.52–0.94, p = 0.019). BTNL2-rs35624343 also significantly reduced the risk of CHD under log-additive genetic model (OR = 0.71, 95% CI = 0.52–0.97, p = 0.032) among smoking participants.

Compared with wild allele C, allele G of BTNL2-rs117896888 was significantly associated with reducing CHD risk among smokers (G vs C, OR = 0.64, 95% CI = 0.42–0.95, p = 0.027). BTNL2-rs117896888 also significantly reduced the risk of CHD under dominant (OR = 0.64, 95% CI = 0.42–0.99, p = 0.043) and log-additive genetic model (OR = 0.65, 95% CI = 0.43–0.97, p = 0.035) among smoking participants.

Compared with wild allele C, allele T of BTNL2-rs41441651 was significantly associated with reducing CHD risk among smokers (T vs C, OR = 0.63, 95% CI = 0.42–0.94, p = 0.023). Genotype (CT) of BTNL2-rs35624343 was significantly associated with reducing CHD risk among smokers (CT vs CC, OR = 0.64, 95% CI = 0.41–0.98, p = 0.041). BTNL2-rs35624343 also significantly reduced the risk of CHD under dominant (OR = 0.63, 95% CI = 0.41–0.97, p = 0.036) and log-additive genetic model (OR = 0.64, 95% CI = 0.42–0.96, p = 0.030) among smoking participants.

Compared with wild allele T, allele C of BTNL2-rs41417449 was significantly associated with reducing CHD risk among smokers (C vs T, OR = 0.63, 95% CI = 0.42–0.95, p = 0.026). Genotype (CT) of BTNL2-rs41417449 was significantly associated with reducing CHD risk among smokers (CT vs TT, OR = 0.64, 95% CI = 0.42–1.00, p = 0.048). BTNL2-rs41417449 also significantly reduced the risk of CHD under dominant (OR = 0.64, 95% CI = 0.41–0.99, p = 0.042) and log-additive genetic model (OR = 0.64, 95% CI = 0.43–0.97, p = 0.034) among smoking participants.

Compared with wild allele T, allele C of BTNL2-rs2076523 was significantly associated with reducing CHD risk among smokers (C vs T, OR = 0.78, 95% CI = 0.63–0.96, p = 0.020). Genotype (CC) of BTNL2-rs2076523 was significantly associated with reducing CHD risk among smokers (CC vs TT, OR = 0.58, 95% CI = 0.37–0.91, p = 0.017). BTNL2-rs2076523 also significantly reduced the risk of CHD under dominant (OR = 0.72, 95% CI = 0.52–0.99, p = 0.041), recessive (OR = 0.66, 95% CI = 0.44–0.99, p = 0.045) and log-additive genetic model (OR = 0.76, 95% CI = 0.61–0.95, p = 0.013) among smoking participants.

This study found that with the exception of BTNL2-rs28362680, the remaining five candidate SNPs were not associated with the risk of CHD in non-smoking participants.

Drinking: the results showed that BTNL2-rs28362680 was significantly associated with CHD susceptibility among drinking or non-drinking participants (Supplemental Table 3). Among drinking participants, -rs28362680 significantly reduced the risk of CHD under multiple genetic models (Allele: OR =0.65, 95% CI = 0.51–0.83, p = 0.001. Genotype (AA): OR = 0.29, 95% CI = 0.13–0.67, p = 0.004; Genotype (GA): OR = 0.68, 95% CI = 0.50–0.94, p = 0.018. Dominant: OR = 0.64, 95% CI = 0.47–0.87, p = 0.004. Recessive: OR = 0.33, 95% CI = 0.15–0.76, p = 0.006. Log-additive: OR = 0.63, 95% CI = 0.48–0.82, p = 0.001). Among non-drinking participants, -rs28362680 significantly reduced the risk of CHD under multiple genetic models (Allele: OR = 0.71, 95% CI = 0.56–0.90, p = 0.005. Genotype (GA): OR = 0.66, 95% CI = 0.48–0.89, p = 0.007. Dominant: OR = 0.64, 95% CI = 0.47–0.86, p = 0.003. Overdominant: OR = 0.69, 95% CI = 0.51–0.93, p = 0.015. Log-additive: OR = 0.68, 95% CI = 0.53–0.88, p = 0.003).

We found no evidence that the five remaining candidate SNPs were associated with the CHD risk in drinking stratified analysis.

CHD complicated with hypertension: the results showed that BTNL2-rs28362680 was significantly associated with the risk of CHD patients complicated with or without hypertension (Supplemental Table 4). Among CHD patients complicated with hypertension, -rs28362680 significantly associated with the CHD risk under all genetic models (Allele: OR =0.69, 95% CI = 0.56–0.84, p < 0.0001. Genotype (AA): OR = 0.35, 95% CI = 0.18–0.68, p = 0.002; Genotype (GA): OR = 0.72, 95% CI = 0.56–0.92, p = 0.010. Dominant: OR = 0.67, 95% CI = 0.53–0.86, p = 0.001. Recessive: OR = 0.40, 95% CI = 0.21–0.77, p = 0.003. Overdominant: OR = 0.78, 95% CI = 0.61–0.99, p = 0.044. Log-additive: OR = 0.67, 95% CI = 0.55–0.83, p < 0.0001). Among CHD patients complicated without hypertension, -rs28362680 significantly associated with the risk of CHD under multiple genetic models (Allele: OR = 0.66, 95% CI = 0.52–0.85, p = 0.001. Genotype (AA): OR = 0.47, 95% CI = 0.23–0.96, p = 0.038; Genotype (GA): OR = 0.62, 95% CI = 0.46–0.84, p = 0.002. Dominant: OR = 0.60, 95% CI = 0.45–0.80, p = 0.001. Overdominant: OR = 0.66, 95% CI = 0.49–0.89, p = 0.005. Log-additive: OR = 0.64, 95% CI = 0.50–0.83, p < 0.0001).

Genotype (CC) of BTNL2-rs2076523 was significantly associated with CHD risk among CHD patients complicated without hypertension (CT vs TT, OR = 0.71, 95% CI = 0.52–0.97, p = 0.032). BTNL2-rs2076523 were also significantly associated with CHD risk under overdominant genetic model (OR = 0.75, 95% CI = 0.56–1.00, p = 0.048) among CHD patients complicated without hypertension.

We found no evidence that the remaining candidate SNPs were associated with the CHD risk in stratified analysis (CHD complicated with hypertension).

CHD complicated with diabetes: the results showed that BTNL2-rs28362680 was significantly associated with the risk of CHD patients complicated with or without diabetes (Supplemental Table 4). Among CHD patients complicated with diabetes, -rs28362680 significantly associated with the CHD risk under multiple genetic models (Allele: OR =0.72, 95% CI = 0.55–0.94, p = 0.014. Genotype (AA): OR = 0.35, 95% CI = 0.13–0.89, p = 0.028. Dominant: OR = 0.72, 95% CI = 0.52–0.99, p = 0.042. Recessive: OR = 0.38, 95% CI = 0.15–0.98, p = 0.024. Log-additive: OR = 0.70, 95% CI = 0.53–0.93, p = 0.012). Among CHD patients complicated without diabetes, -rs28362680 significantly associated with the CHD risk under all genetic models (Allele: OR =0.66, 95% CI = 0.55–0.80, p < 0.0001. Genotype (AA): OR = 0.42, 95% CI = 0.23–0.75, p = 0.003; Genotype (GA): OR = 0.64, 95% CI = 0.51–0.82, p < 0.0001. Dominant: OR = 0.61, 95% CI = 0.48–0.77, p < 0.0001. Recessive: OR = 0.50, 95% CI = 0.28–0.89, p = 0.014. Overdominant: OR = 0.68, 95% CI = 0.54–0.87, p = 0.002. Log-additive: OR = 0.64, 95% CI = 0.53–0.78, p < 0.0001).

Compared with wild allele G, allele A of BTNL2-rs35624343 was significantly associated with the risk of CHD patients complicated without diabetes (A vs G, OR = 0.77, 95% CI = 0.61–0.97, p = 0.027). Genotype (GA) of BTNL2-rs35624343 was significantly associated with the risk of CHD patients complicated without diabetes (GA vs GG, OR = 0.74, 95% CI = 0.57–0.97, p = 0.026). BTNL2-rs35624343 was also significantly associated with the risk of CHD under dominant (OR = 0.74, 95% CI = 0.57–0.96, p = 0.020), Overdominant (OR = 0.75, 95% CI = 0.58–0.97, p = 0.030) and log-additive genetic models (OR = 0.76, 95% CI = 0.60–0.96, p = 0.022) among CHD patients complicated without diabetes.

Compared with wild allele C, allele G of BTNL2-rs117896888 was significantly associated with the risk of CHD patients complicated without diabetes (G vs C, OR = 0.65, 95% CI = 0.47–0.91, p = 0.010). Genotype (CG) of BTNL2-rs117896888 was significantly associated with the risk of CHD patients complicated without diabetes (CG vs CC, OR = 0.67, 95% CI = 0.47–0.95, p = 0.024). BTNL2-rs117896888 was also significantly associated with the risk of CHD under dominant (OR = 0.65, 95% CI = 0.46–0.92, p = 0.013), Overdominant (OR = 0.67, 95% CI = 0.47–0.95, p = 0.024) and log-additive genetic models (OR = 0.65, 95% CI = 0.47–0.90, p = 0.009) among CHD patients complicated without diabetes.

Compared with wild allele C, allele T of BTNL2-rs41441651 was significantly associated with the risk of CHD patients complicated without diabetes (T vs C, OR = 0.64, 95% CI = 0.46–0.89, p = 0.008). Genotype (CT) of BTNL2-rs41441651 was significantly associated with the risk of CHD patients complicated without diabetes (CT vs CC, OR = 0.66, 95% CI = 0.46–0.93, p = 0.018). BTNL2-rs41441651 was also significantly associated with the risk of CHD under dominant (OR = 0.64, 95% CI = 0.45–0.90, p = 0.010), Overdominant (OR = 0.66, 95% CI = 0.47–0.94, p = 0.018) and log-additive genetic models (OR = 0.64, 95% CI = 0.46–0.89, p = 0.007) among CHD patients complicated without diabetes.

Compared with wild allele T, allele C of BTNL2-rs41417449 was significantly associated with the risk of CHD patients complicated without diabetes (C vs T, OR = 0.65, 95% CI = 0.46–0.90, p = 0.009). Genotype (CT) of BTNL2-rs41417449 was significantly associated with the risk of CHD patients complicated without diabetes (CT vs TT, OR = 0.66, 95% CI = 0.47–0.94, p = 0.020). BTNL2-rs41417449 was also significantly associated with the risk of CHD under dominant (OR = 0.64, 95% CI = 0.45–0.91, p = 0.011), Overdominant (OR = 0.66, 95% CI = 0.47–0.94, p = 0.020) and log-additive genetic models (OR = 0.64, 95% CI = 0.46–0.90, p = 0.008) among CHD patients complicated without diabetes.

Compared with wild allele T, allele C of BTNL2-rs2076523 was significantly associated with the risk of CHD patients complicated without diabetes (C vs T, OR = 0.82, 95% CI = 0.69–0.96, p = 0.017). Genotype CT (CT vs TT, OR = 0.69, 95% CI = 0.48–0.99, p = 0.045) and CC (CC vs TT, OR = 0.74, 95% CI = 0.58–0.95, p = 0.020) of BTNL2-rs2076523 was significantly associated with the risk of CHD patients complicated without diabetes. BTNL2-rs2076523 was also significantly associated with the risk of CHD under dominant (OR = 0.73, 95% CI = 0.58–0.93, p = 0.010) and log-additive genetic models (OR = 0.81, 95% CI = 0.68–0.96, p = 0.014) among CHD patients complicated without diabetes.

This study found that with the exception of BTNL2-rs28362680, the remaining five candidate SNPs were not associated with the risk of CHD among CHD patients complicated with diabetes.

FPRP Analysis

The results of FPRP analysis showed that a very few positive results may be accidental, and most of them were noteworthy findings at the prior probability level of 0.25 and FPRP threshold of 0.2 (Supplemental Table 5). And the statistical power value of noteworthy positive results in this study reached more than 90%.

Specifically, under certain genetic models, BTNL2-rs28362680 was potentially associated with the CHD risk among women (genotype AA: prior probability = 0.269; recessive: prior probability = 0.351), > 62 years old (genotype AA: prior probability = 0.204), smoking (recessive: prior probability = 0.245), non-smoking participants (recessive: prior probability = 0.219), as well as CHD patients complicated with diabetes (genotype AA: prior probability = 0.266; recessive: prior probability = 0.323). The above positive results may not be worth noting. The remaining positive findings were all noteworthy in this study.

Haplotype Analysis

The results are shown in Figure 1, color represents the degree of linkage imbalance, and the redder the color, the stronger the linkage disequilibrium. We can see that the four candidate BTNL2 SNPs (rs117896888, rs41441651, rs41417449, rs28362680) can composed one linkage disequilibrium block. And the results of haplotype analysis showed that the haplotype “C_rs117896888C_rs41441651T_rs41417449A_rs28362680” (OR = 0.65, CI = 0.53–0.80, p < 0.0001) and “G_rs117896888T_rs41441651C_rs41417449A_rs28362680” (OR = 0.68, CI = 0.51–0.92, p = 0.013) all reduced the susceptibility of CHD among southern Chinese Han population (Table 4).

Table 4 Haplotype Analysis of Candidate BTNL2 Genetic Polymorphisms with CHD Risk

Figure 1 Haplotype block map for the BTNL2 SNPs (rs117896888, rs41441651, rs41417449, rs28362680). (A) The numbers inside the diamonds indicate The D’ for pairwise analyses. (B) The numbers inside the diamonds indicate the R2 for pairwise analyses. The colors represent the degree of linkage disequilibrium: the redder the color, the stronger the linkage disequilibrium.

MDR Analysis

The MDR results showed that the single locus model (rs28362680) has the highest test accuracy (Table 5). As is shown in Figure 2A, the dendrogram has described the interaction between the six candidate SNPs. The color of the lines in the dendrogram represents the level of redundancy or synergy. The closer to the red, the stronger the synergy, and the closer to the blue, the stronger the redundancy. In addition, information gain (IG) is used to evaluate attribute interactions. As shown in Figure 2B, the IG value of rs28362680 is the highest. Therefore, we chose the unit locus model composed of rs28362680 as the best model for predicting CHD risk, with a test accuracy of 0.558 and a perfect CVC = 10/10.

Table 5 BTNL2 SNP–SNP Interaction Models Analyzed by the MDR Method

Figure 2 Multifactor dimensionality reduction (MDR) analysis of interaction between the candidate SNPs of BTNL2 ((rs35624343, rs117896888, rs41441651, rs41417449, rs28362680 and rs2076523). (A) SNP-SNP Interaction Dendrogram: the color represents the degree of redundancy or synergy between SNP-SNP; the closer the color is to red, the more synergy, and the closer to blue, the more redundancy. (B) Fruchterman-Reingold: values in nodes represent the IGs of individual attribute (main effects). Values between nodes are IGs of each pair of attributes (interaction effects).

Discussion

In this study, we conducted an association study between BTNL2 genetic polymorphism and CHD susceptibility in southern Chinese Han population among 1419 participants. Participants were grouped according to the common CHD risk factors, and the association between BTNL2 genetic polymorphism and potential CHD risk factors was evaluated. Our results showed a strong association between BTNL2-rs28362680 and CHD susceptibility. Compared with the allele G and genotype GG, participants with allele A and genotype AA, GA of BTNL2-rs28362680 had lower risk of CHD (OR < 1). More importantly, the BTNL2-rs28362680 genetic polymorphism was significantly associated with a decreased risk of CHD in both overall and stratified analyses. Our study suggest that the allele A of BTNL2-rs28362680 can reduce the incidence of CHD in southern Chinese Han population.

It is worth noting that the results of all stratified analyses showed that the BTNL2-rs28362680 was significantly associated with the CHD susceptibility, which was not affected by other potential risk factors. That is to say, gender (male or female), age (≤62 years oldor >62 years old), smoking (smoking or non-smoking), drinking (drinking or non-drinking), hypertension (CHD complicated with hypertension or not), diabetes (CHD complicated with diabetes or not) cannot change the fact that BTNL2-rs28362680 was significantly associated with CHD susceptibility among participants. However, we all known that the occurrence and development of CHD will be affected by a variety of external factors.^18,19 Hajek et al have reported that the genetic risk score of CHD can predict the risk of cardiovascular disease in men but not women.²⁰ Many studies have shown that age and gender are the main determinants of CHD risk.^21,22 Male,^23,24 smokers,²⁵ and older people²⁶ are at higher risk of CHD. In addition, compared with non-diabetic patients, diabetes is the main risk factor for coronary artery disease (CAD).²⁷ Hypertension induces endothelial dysfunction and exacerbates the process of atherosclerosis.²⁸ And hypertension is also a risk factor for CHD.¹⁸ However, the results of our study showed that the A allele of BTNL2-rs28362680 was significantly associated with reduction of the CHD risk in all subgroup analyses. And this protection was not affected by other risk factors. So far, only one study has reported that BTNL2-rs28362680 was a new risk site for Crohn’s disease (CD) in Koreans.²⁹ This study is the first to find that BTNL2-rs28362680 is a protective factor of CHD in Southern Chinese Han population. Based on the FPRP analysis and large sample size of this study, our new findings are very noteworthy.

In addition to finding that BTNL2-rs28362680 was significantly associated with participants’ CHD susceptibility, we also found that BTNL2-rs35624343, -rs117896888, -rs41441651 and -rs41417449 were associated with CHD susceptibility in the overall analysis and multiple stratified analyses. In this study, the above four candidate SNPs were significantly associated CHD risk reduction in the overall analysis, and similar findings were found in stratified analysis (male, smoking, and CHD patients without diabetes). A large-scale genome-wide association study showed that CHD and diabetes share a common genetic background.³⁰ Studies have reported that BTNL2 genetic polymorphism is associated with the susceptibility of diabetes in Caucasians.³¹ Combined with the results of our study, we speculate that BTNL2-rs35624343, -rs117896888, -rs41441651 and -rs41417449 may reduce the risk of diabetes and thereby reduce the risk of CHD. Secondly, male and smoking have been reported as risk factors for CHD. Based on the results of this study, we speculate that these four candidate SNPs are protective genetic factors for CHD in southern Chinese Han population, and this protective effect will not be affected by potential risk factors for CHD. Our study was the first to find out the evidence that BTNL2-rs35624343, -rs117896888, -rs41441651 and -rs41417449 were significantly associated with CHD susceptibility. In this study, we identified some genetic variants that protect the risk of coronary heart disease in Southern Chinese Han population, and these SNPs are reported for the first time.

In previous studies on the association between BTNL2 and coronary artery, Hsueh et al found that BTNL2 gene polymorphism may be a genetic marker for the susceptibility of Taiwanese children to Kawasaki disease and the risk of coronary complications.³² Kawasaki disease is a childhood vasculitis that particularly affects the coronary arteries, and the risk of ischemic heart disease increases in the future.³² Previous studies have suggested that BTNL2 is potentially associated with the occurrence and development of coronary heart disease and this speculation has also been confirmed in the results of this study. Subsequently, Lu et al conducted a genome-wide association study on patients with coronary heart disease and found that the C6orf10-BTNL2 gene was a susceptibility gene for CHD in the Southern Chinese Han population, and the mutation of the BTNL2-rs9268402 locus can lead to the occurrence of coronary heart disease.³³ In previous studies, whether it is Kawasaki disease or coronary heart disease, BTNL2 genetic polymorphism is the susceptibility genetic locus of the disease (OR> 1). However, the candidate BTNL2 genetic locus was associated with a reduced risk of CHD in this study. We speculate that this may be due to the differences of genetic background and geographical of the study subjects.

The pathological basis of coronary heart disease is atherosclerosis, which is a chronic inflammatory disease. CHD is characterized by the formation of arterial plaques composed mainly of lipids, calcium, and inflammatory cells.³⁴ Atherosclerotic plaques are rich in immune cells, which can dominate and influence the inflammatory response.³⁵ T cells are crucial to the development of atherosclerotic plaques in the development of CHD.³⁶ Helper T cell 1 (Th1)/Th17-mediated pro-inflammatory response can aggravate atherosclerosis, while regulatory T cells (T reg) play a protective role against atherosclerosis by limiting inflammation and balancing plaque formation.³⁷ BTNL2 can inhibit the activation of T cells and reduce the secretion of pro-inflammatory cytokines, leading to immune dysfunction.³⁸ In chronic inflammation, CD4 T cells can secrete large amounts of Th1 cytokines, TNF-α and IFN-γ,³⁹ which are effective atherogenic cytokines.⁴⁰ Sinisalo et al found that the level of BTNL2 mRNA were higher in coronary artery samples from subjects with acute coronary artery disease. At the same time, the study also found that inhibition of BTNL2 significantly enhanced the proliferation of CD4(+)FOXP3(+) T reg cells.⁴¹ Considering the results of previous studies and our study, we analyzed the reasons why these missense variants in BTNL2 in this study are associated with the reduction of CHD risk among participants in this study. We speculate that these missense variants (rs35624343 A/G, rs117896888 G/C, rs41441651 T/C, rs41417449 C/T, rs28362680 A/G, rs2076523 C/T) in BTNL2 may play a protective role in atherosclerosis by regulating BTNL2 mRNA expression and participating in a certain inflammatory pathway, thereby limiting inflammatory response and balancing plaque formation. However, the specific molecular mechanism needs to be further studied.

Our study provides data supplement for the genetic loci associated with CHD susceptibility in Han population from southern China. Although the mechanism of these missense variants in coronary heart disease is still unclear, the study provides new ideas for the risk assessment and clinical individualized prevention and treatment of CHD in Han population from southern China. In addition, we found these missense variations in BTNL2 could cause changes in amino acid sequences through dbSNP database. Changes in amino acid sequence can lead to changes in protein structure, which is directly related to its function.⁴² Combined with the results of our study, we speculated that these missense variants caused changes in amino acid sequence, which resulted in structural changes of BTNL2 protein, thus affecting its biological functions. Accordingly, we proposed a hypothesis that these missense variants are protective factors for CHD in Chinese Han population, which may be achieved by affecting the amino acid sequence of BTNL2 and then changing its protein structure and function. But the above is just speculation, and further mechanistic studies are needed to confirm it. In any case, this study has laid a solid theoretical foundation for the mechanism of BTNL2 in the occurrence and development of CHD.

However, we must face the fact that there are certain limitations in this study. In order to ensure the reliability and repeatability of the results, a large sample size verification study is necessary. Secondly, studies on the mechanism of BTNL2 genetic polymorphism in CHD development need to be further determined. Nevertheless, to the best of our knowledge, this study is the first to investigate the association between BTNL2 genetic polymorphism and CHD susceptibility in China and to find positive results. This study will lay a reliable theoretical foundation for the study of the mechanism of BTNL2 in the occurrence and development of CHD, and provide new ideas for the risk assessment and clinical individualized prevention and treatment of CHD in Chinese Han population.

Conclusion

In summary, this study showed that missense variants (rs35624343, rs117896888, rs41441651, rs41417449, rs28362680, rs2076523) in BTNL2 were associated with coronary heart disease susceptibility. In particular, there is sufficient evidence that BTNL2-rs28362680 was associated with a reduced risk of coronary heart disease in the southern Chinese Han population. This study has laid a reliable theoretical foundation for the study of the mechanism of BTNL2 in the occurrence and development of coronary heart disease.