ORIGINAL RESEARCH article

Front. Cardiovasc. Med., 06 October 2022

Sec. Atherosclerosis and Vascular Medicine

Volume 9 - 2022 | https://doi.org/10.3389/fcvm.2022.964978

Rs420137, rs386360 and rs7763726 polymorphisms in fibronectin type III domain containing 1 are associated with susceptibility to coronary heart disease: Analysis in the Han population

    XH

    Xiaodan He

    XL

    Xuemei Li

    XD

    Xiaoyan Du

    JH

    Jianlun Han

    HZ

    Hui Zhang

    YZ

    Yan Zhu

    HM

    Honghong Ma *

  • Department of Cardiovasology, The First Hospital of Yulin, Yulin, China

Article metrics

View details

2

Citations

1,3k

Views

599

Downloads

Abstract

Background:

Numerous genetic studies have shown that genes are related to the pathogenesis of coronary heart disease (CHD). The main aim of this study was to confirm whether fibronectin type III domain containing 1 (FNDC1) polymorphisms correlate with the risk of CHD.

Methods:

In this study, in order to assess the association between three FNDC1 single nucleotide polymorphisms (SNPs) and the risk of CHD, we conducted a case-control study involving 630 patients with CHD and 568 healthy controls using Agena MassARRAY (Agena Bioscience, San Diego, CA, USA). Genotype distribution in case and control groups was analyzed by Chi square test. Odds ratios (ORs) and 95% confidence intervals (CIs) were calculated by logistic regression models adjusted for age, sex, smoking, and alcohol consumption to assess the correlation between SNPs and CHD risk.

Results:

Our results indicated that FNDC1-rs420137, -rs386360, and -rs7763726 played important roles in enhancing the risk of CHD. Subgroup analysis revealed that rs420137 increased the susceptibility to CHD in males, smokers, and patients aged ≤62 years. Rs360 had an increased risk of CHD in males, patients at aged ≤62 years, smokers, and non-drinkers. Furthermore, the association of rs7763726 with increased CHD risk was also observed in males, patients aged ≤62 years, smokers, and drinkers. Last but not least, these three SNPs we selected were protective factors against hypertension in CHD individuals.

Conclusion:

Our research suggest that FNDC1-rs420137, -rs386360, and -rs7763726 variants may be regarded as novel biomarkers for predicting CHD risk and other specific mechanisms of action of CHD need to be further studied.

Introduction

Coronary heart disease (CHD) is a very common cardiovascular disease. It is estimated that the total number of deaths from cardiovascular disease is expected to increase to 23.3 million by 2030 (1). Since the 1970 s, CHD has been the main cause of death all over the world. In China, 700,000 people die of CHD annually (2). In addition, the morbidity of CHD is rising. The occurrence and development of CHD correlate with a variety of environmental factors, including social factors (aging, employment, income, and education), behavioral factors (dietary habits, drinking, smoking, and exercise), and metabolic factors obesity, diabetes mellitus (DM), triglyceride (TG), and high-density lipoprotein cholesterol (HDL-C) (3). However, an increasing number of studies have shown that genetic factors play a critical role in the progression of CHD. By far, more than 60 genetic variants have been identified to be related to an increased susceptibility to CHD in genome-wide association study (GWAS) analyses (4). Interestingly, studies have shown that single nucleotide polymorphisms (SNPs) can affect the level of gene production or the binding of targets, thus influencing the normal expression of genes and leading to the occurrence of diseases (5, 6).

Fibronectin type III domain containing 1 (FNDC1) is a protein-encoding and disease-related gene located in the 6q25.3 region of the human chromosome and contains 102,713 bases, with a length of ~1,894 and a molecular weight of 205.6 kDa (7, 8). Moreover, FNDC1, also known as activator of G-protein signaling 8 (AGS8), is a positive factor in modulating the signal transduction system. It is expressed in various types of tissues, including thyroid, heart, fetal cartilage, adipose, and kidney tissues (911). As an important extracellular matrix protein, FNDC1 plays a regulatory role in tumorigenesis and participates in cell migration, proliferation, and apoptosis. It is also involved in the pathogenesis of multiple blood and cardiovascular diseases (12). Previous studies have also shown that FNDC1 has a certain relationship with hypertension (13).

To our knowledge, the effects of FNDC1 gene polymorphisms on CHD have not been studied. Therefore, the present study explored the linkage between FNDC1 polymorphisms and CHD susceptibility, which has considerable implications for the diagnosis and treatment of CHD.

Materials and methods

Study population

A total of 630 CHD cases and 568 controls were recruited from the First Hospital of Yulin City. All patients were independently judged by at least two cardiologists and diagnosed as CHD according to blood examination, coronary angiography, standardized electrocardiogram (ECG), and echocardiography. Out-patients in the control group were from the same hospital during the same period as the cases. According to previous medical history and angiographic diagnosis, they had no CHD, no obvious ischemic changes on ECG and no symptoms of chest pain. Participants with other heart diseases (rheumatic heart disease, cardiomyopathy, concomitant heart valve disease, or congenital heart disease) were eliminated from the study. Furthermore, it should be noted that patients with a history of intravenous thrombolysis, coronary stenting, angioplasty bypass surgery, or coronary artery disease were also excluded. Moreover, clinical indicators and basic characteristics of subjects were collected from medical records.

DNA extraction

We collected 5 mL of fasting peripheral venous blood from each participant and put it in a test tube containing ethylene diamine tetraacetic acid (EDTA). According to the manufacturer's instructions, we extracted DNA from the collected blood samples using GoldMag-Mini whole blood genomic DNA purification kits (GoldMag Co., Ltd., Xi'an City, China).

Genotyping of SNPs

First of all, three variants of FNDC1 (rs420137, rs386360, and rs7763726) were selected from the HapMap database according to the minor allele frequency (MAF) >5% and the unbalanced r2 value >0.8. Then, we genotyped them by the Agena MassARRAY platform (Agena Bioscience, San Diego, CA, USA). Finally, the results of genotyping were processed and analyzed by Agena Bioscience Typer software (version 4.0).

Data analysis

The basic statistical analysis was conducted by Microsoft Excel (Microsoft Corp., Redmond, WA, USA) and SPSS 20.0 software (SPSS, Chicago, IL, USA). The t-test was used to assess continuous variables, while the Chi-square test was utilized to access categorical variables between samples. Meanwhile, the Chi-square test was used to evaluate whether the genotype distribution of these SNPs in the control group was in accordance with Hardy-Weinberg equilibrium (HWE), and the differences in the genotype and allele frequencies of these SNPs in two groups were also analyzed. PLINK 1.9 software (version 1.07) and Haploview 4.2 software were applied to analyze the association of genetic variants with CHD susceptibility under five genetic models, and odds ratios (ORs) and 95% confidence intervals (CIs) were calculated to confirm differences in this study. Besides, multifactor dimensionality reduction (MDR) was used to analyze the effects of SNP-SNP interactions on CHD susceptibility. Finally, a p-value < 0.05 indicated a statistical difference.

Results

Characteristics of subjects

The information on CHD patients is presented in Table 1. Patients with an average age of 62.16 ± 10.37 years included 429 males and 201 females, while controls with an average age of 61.46 ± 8.96 years included 360 males and 208 females. There was no statistically significant difference in the distribution of sex and age between the two groups (p > 0.05). In addition, the levels of HDL-C, TG, urea, alanine amino transferase (ALT), and blood glucose (BG) had no significant differences (p > 0.05) between controls and cases (Table 1). Nevertheless, other indicators [white blood cell count (WBC), red blood cell count (RBC), low-density lipoprotein cholesterol (LDL-C), platelet (PLT), hemoglobin (Hgb), total cholesterol (TC), absolute value of lymphocytes (LYMP), absolute value of neutrophils (NEUT), creatinine, and aspartate aminotransferase (AST)] were significantly different (p < 0.05) between the case and control groups (Table 1).

Table 1

VariablesCases
(N = 630)
Controls
(N = 568)
P
Age, years (mean ± SD)62.16 ± 10.3761.46 ± 8.960.252
>62 years305 (48.4%)289 (50.9%)
≤62 years325 (51.6%)279 (49.1%)
Gender0.948
Male429 (68.1%)360 (63.4%)
Female201 (31.9%)208 (36.6%)
Smoking
Yes330 (52.4%)260 (45.8%)
No300 (47.6%)308 (54.2%)
Drinking
Yes293 (46.5%)266 (46.8%)
No337 (53.5%)302 (53.2%)
Complicating hypertension
Yes397 (63.0%)152 (50.0%)
No233 (37.0%)152 (50.0%)
Complicated with diabetes
Yes180 (28.6%)
No450 (71.4%)
WBC6.82 ± 2.045.83 ± 1.47<0.001
RBC4.43 ± 1.9956.5 ± 10.5<0.001
HDL-C (mmol/L)1.10 ± 0.261.13 ± 0.230.100
LDL-C (mmol/L)2.40 ± 0.892.61 ± 0.74<0.001
PLT (109/L)194.00 ± 56.87213.26 ± 59.85<0.001
HGB131.10 ± 27.44146.69 ± 14.34<0.001
TG (mmol/L)1.62 ± 1.001.75 ± 1.470.075
TC (mmol/L)4.04 ± 1.084.75 ± 0.98<0.001
Urea5.26 ± 2.055.20 ± 1.230.513
LYMP (10E10/L)8.06 ± 11.765.83 ± 1.47<0.001
NEUT (10E9/L)22.4 ± 28.437.55 ± 2.94<0.001
Creatinine (μmol/L)82.40 ± 20.9768.26 ± 13.15<0.001
AST29.37 ± 33.5325.59 ± 18.770.009
ALT25.73 ± 22.0325.63 ± 18.470.957
BG (mmol/L)5.96 ± 2.616.06 ± 1.630.459

Characteristic of patients with CHD and controls.

SD, standard deviation; WBC, white blood cell count; RBC, red blood cell count; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; PLT, platelet; HGB, hemoglobin; TG, triglyceride; TC, total cholesterol; LYMP, lymphocyte absolute value; NEUT, neutrophil absolute value; AST, aspartate aminotransferase; ALT, alanine aminotransferase; BG, blood glucose. Bold values are statistically significant (p <0.05).

Association between FNDC1 polymorphisms and CHD risk

Basic information on three SNPs (rs420137, rs386360, and rs7763726) in FNDC1 is displayed in Table 2. All SNPs in the control group were in line with HWE (p > 0.05) and could be further analyzed. At the same time, we found that the frequency of the rs420137-G allele in the case group was significantly higher than that in the control group (p = 0.013, OR: 1.25, 95% CI: 1.05–1.48). The “C” allele of rs386360 (p = 0.010, OR: 1.26, 95% CI: 1.05–1.48) and the “A” allele of rs7763726 (p = 0.001, OR: 1.36, 95% CI: 1.13–1.64) were obviously related to an increased susceptibility to CHD.

Table 2

SNP IDGeneChr: positionAlleleMAFHWEORP
(minor/major)CaseControlp(95% CI)
rs420137FNDC16: 159231899G / C0.3300.2830.8361.25 (1.05–1.48)0.013
rs386360FNDC16: 159239847C / A0.3160.2690.6711.26 (1.05–1.50)0.010
rs7763726FNDC16:159249068G / A0.2810.2230.6291.36 (1.13–1.64)0.001

Basic information of SNPs in FNDC1 gene.

SNP, single nucleotide polymorphism; MAF, minor allele frequency; HWE, Hardy-Weinberg equilibrium; OR, odds ratio; 95% CI, 95% confidence interval.

P-values were calculated from χ2 test.

Bold values are statistically significant (p < 0.05).

As can be seen from Table 3, logistic regression analysis revealed that rs420137-GC was correlated with an increased risk of CHD under the co-dominant (p = 0.037, OR: 1.29, 95% CI: 1.02–1.65), dominant (p = 0.015, OR: 1.34, 95% CI: 1.06–1.68), and additive models (p = 0.011, OR: 1.27, 95% CI: 1.06–1.51). Rs386360-CA had a strong linkage with a higher susceptibility to CHD (co-dominant model: p = 0.026, OR: 1.31, 95% CI: 1.03–1.67; dominant model: p = 0.011, OR: 1.35, 95% CI: 1.07–1.70; additive model: p = 0.009, OR: 1.28, 95% CI: 1.07–1.53). It also showed that rs7763726-GA elevated CHD risk (co-dominant model: p = 0.001, OR: 1.50, 95% CI: 1.17–1.90; dominant model: p = 4.0E-04, OR: 1.52, 95% CI: 1.21–1.92; additive model: p = 7.0E-04, OR: 1.39, 95% CI: 1.15–1.69).

Table 3

SNP-IDModelGenotypeCase (%)Control (%)OR (95% CI)PAdjusted OR (95% CI)P
rs420137Co-dominantCC290 (51.1)279 (44.3)11
GC234 (41.2)286 (45.4)1.27 (1.00–1.61)0.0491.29 (1.02–1.65)0.037
GG44 (7.8)65 (10.3)1.54 (1.01–2.33)1.55 (1.02–2.36)
DominantCC290 (51.1)279 (44.3)11
GC-GG278 (48.9)351 (55.7)1.31 (1.05–1.65)0.0191.34 (1.06–1.68)0.015
RecessiveCC-GC524 (92.2)565 (89.7)11
GG44 (7.8)65 (10.3)1.37 (0.92–2.05)0.1241.38 (0.92–2.06)0.132
Additive1.04 (0.65–1.68)0.0131.27 (1.06–1.51)0.011
rs386360Co-dominantAA301 (53)289 (45.9)11
CA228 (40.1)283 (44.9)1.29 (1.02–1.64)0.0351.31 (1.03–1.67)0.026
CC39 (6.9)58 (9.2)1.55 (1.00–2.40)1.56 (1.01–2.43)
DominantAA301 (53)289 (45.9)11
CA-CC267 (47)341 (54.1)1.33 (1.06–1.67)0.0141.35 (1.07–1.70)0.011
RecessiveAA-CA529 (93.1)572 (90.8)11
CC39 (6.9)58 (9.2)1.38 (0.90–2.10)0.1401.38 (0.90–2.11)0.149
Additive1.27 (1.06–1.52)0.0101.28 (1.07–1.53)0.009
rs7763726C-odominantAA345 (60.7)321 (51)11
GA193 (34)264 (41.9)1.47 (1.16–1.87)0.0021.50 (1.17–1.90)0.001
GG30 (5.3)45 (7.1)1.61 (0.99–2.62)1.68 (1.03–2.74)
DominantAA345 (60.7)321 (51)11
GA-GG223 (39.3)309 (49)1.49 (1.18–1.87)7.0E-041.52 (1.21–1.92)4.0E-04
RecessiveAA-GA538 (94.7)585 (92.9)11
GG30 (5.3)45 (7.1)1.38 (0.86–2.22)0.1861.43 (0.88–2.31)0.176
Additive1.37 (1.13–1.65)0.0011.39 (1.15–1.69)7.0E-04

Genotypes frequencies distribution and association with CHD risk.

CHD, coronary heart disease; SNP, single nucleotide polymorphism; OR, odds ratio; CI, confidence interval. Adjusted OR (95% CI) were calculated by logistic regression analysis with adjustments for age and gender. Bold values are statistically significant (p < 0.05).

Subgroup analysis of the correlation between FNDC1 polymorphisms and CHD risk

We also performed stratified analyses of the association of three SNPs with CHD risk according to sex, age, smoking, drinking, hypertension, and diabetes to determine whether traditional risk factors influence the role of genetic variants in CHD. According to Table 4, we found that all variants increased the incidence of CHD under the allele, co-dominant, dominant, and additive models in males (rs420137: p = 0.006, 0.010, 0.004, and 0.006, respectively; rs386360: p = 0.007, 0.017, 0.006, and 0.006, respectively; rs7763726: p < 0.001, 0.001, 0.001, and 0.001, respectively). Rs420137 correlated with CHD risk in the allele (p = 0.047, OR: 1.29, 95% CI: 1.00–1.65), co-dominant (GG vs. CC p = 0.026, OR: 1.83, 95% CI: 1.01–3.31), recessive (p = 0.043, OR: 1.68, 95% CI: 0.95–2.99), and additive (p = 0.031, OR: 1.29, 95% CI: 1.00–1.65) models in individuals aged 62 years and younger.

Table 4

SNP-IDModelGenotypeMaleFemale>6262
OR (95% CI)POR (95% CI)POR (95% CI)POR (95% CI)P
rs420137AlleleC1111
G1.36 (0.11–1.09)0.0061.01 (0.81–1.46)0.5731.21 (0.95–1.55)0.1211.29 (1.00–1.65)0.047
Co-dominantCC1111
GC1.48 (1.10–1.99)0.0101.00 (0.66–1.50)0.9971.37 (0.98–1.92)0.1421.20 (0.85–1.69)0.300
GG1.65 (1.00–2.74)1.38 (0.66–2.90)1.32 (0.72–2.39)1.83 (1.01–3.31)0.026
DominantCC1111
GC-GG1.51 (1.14–2.00)0.0041.05 (0.71–1.55)0.8141.36 (0.98–1.88)0.1251.30 (0.94–1.80)0.100
RecessiveCC-GC1111
GG1.38 (0.85–2.25)0.1941.38 (0.68–2.82)0.3721.12 (0.63–1.99)0.6491.68 (0.95–2.99)0.043
Additive1.36 (1.09–1.69)0.0061.10 (0.80–1.50)0.5621.23 (0.96–1.59)0.1621.29 (1.00–1.65)0.031
rs386360AlleleA1111
C1.36 (1.09–1.69)0.0071.11 (0.83–1.50)0.4981.24 (0.96–1.58)1.28 (0.13–1.00)0.051
Co-dominantA/A1111
C/A1.44 (1.06–1.93)0.0171.12 (0.75–1.68)0.5021.39 (0.99–1.95)0.1241.22 (0.87–1.72)0.279
C/C1.73 (1.01–2.94)1.31 (0.60–2.84)1.35 (0.72–2.52)1.79 (0.96–3.33)0.027
DominantA/A1111
C/A-C/C1.48 (1.12–1.96)0.0061.14 (0.77–1.69)0.5081.38 (1.00–1.91)0.1051.30 (0.94–1.80)0.095
RecessiveA/A-C/A1111
C/C1.47 (0.88–2.47)0.1421.23 (0.58–2.61)0.5841.15 (0.63–2.11)0.6011.64 (0.89–3.00)0.044
Additive1.37 (1.09–1.70)0.0061.13 (0.82–1.55)0.4451.26 (0.97–1.63)0.1311.29 (1.00–1.66)0.030
rs7763726AlleleA1111
G1.58 (0.12–1.25)< 0.0011.09 (0.80–1.48)0.2871.33 (1.03–1.73)0.0311.40 (0.14–1.08)0.013
Co-dominantA/A1111
A/G1.87 (1.38–2.53)< 0.0010.98 (0.65–1.48)0.3371.45 (1.04–2.04)0.0881.52 (1.08–2.15)0.020
G/G1.78 (0.96–3.29)1.48 (0.66–3.30)1.61 (0.80–3.24)1.75 (0.88–3.49)
DominantA/A1111
G/A-G/G1.86 (1.39–2.48)< 0.0011.04 (0.71–1.54)0.8311.47 (1.06–2.04)0.0551.55 (1.12–2.16)0.009
RecessiveA/A-G/A1111
G/G1.39 (0.76–2.56)0.2831.49 (0.68–3.26)0.3171.37 (1.00–1.93)0.3161.49 (0.76–2.94)0.182
Additive1.59 (1.25–2.03)< 0.0011.10 (0.80–1.50)0.5621.36 (1.04–1.78)0.0501.42 (1.09–1.86)0.008

FNDC1 SNPs associated with CHD risk stratified by gender and age status.

CHD, coronary heart disease; SNP, single nucleotide polymorphism; OR, odds ratio; CI, confidence interval.

OR (95% CI) values were calculated by logistic regression analysis with adjustments for gender and age.

Bold values are statistically significant (p < 0.05).

Table 4 also shows that rs386360 had a higher risk of CHD under the co-dominant, recessive, and additive models (p = 0.027, 0.044, and 0.030, respectively) in patients aged ≤62 years. Rs7763726 had an increased risk of CHD in patients aged ≤62 years in the allele (p = 0.013, OR: 1.40, 95% CI: 0.14–1.08), co-dominant (p = 0.020, OR: 1.52, 95% CI: 1.08–2.15), dominant (p = 0.009, OR: 1.55, 95% CI: 1.12–2.16), and additive models (p = 0.008, OR: 1.42, 95% CI: 1.09–1.86).

Stratified analyses by smoking and drinking (Table 5) demonstrated that rs420137 was significantly related to a higher risk of CHD under the allele model in smokers (p = 0.006) and drinkers (p = 0.044), and under the dominant (p = 0.017), recessive (p = 0.036), and additive models (p = 0.006) in smokers. Rs386360 increased the risk of CHD in smokers (allele model: p = 0.004, OR:1.46, 95% CI: 0.13–1.13; co-dominant model with the “CA vs. AA” genotype: p = 0.045, OR: 1.43, 95% CI: 1.00–2.02; dominant model: p = 0.010, OR: 1.55, 95% CI: 1.11–2.16; recessive model: p = 0.032, OR: 2.08, 95% CI: 1.06–4.08; additive model: p = 0.003, OR: 1.49, 95% CI: 1.14–1.95) and non-drinkers (co-dominant model with the “CA vs. AA” genotype: p = 0.027, OR: 1.45, 95% CI: 1.04–2.01; dominant model: p = 0.028, OR: 1.43, 95% CI: 1.04–1.95). In addition, rs7763726 enhanced the incidence of CHD in smokers (allele model: p < 0.001, OR: 1.83, 95% CI: 0.14–1.38; co-dominant model with the “AG vs. AA” genotype: p < 0.001, OR: 2.13, 95% CI: 1.48–3.05; dominant model: p < 0.001, OR: 2.19, 95% CI: 1.55–3.10; additive model: p < 0.001, OR: 1.91, 95% CI: 1.42–2.57) and drinkers (allele model: p = 0.002, OR: 1.56, 95% CI: 1.18–2.06; co-dominant model with the “AG vs. AA” genotype: p = 0.002, OR: 1.75, 95% CI: 1.22–2.51; dominant model: p = 0.001, OR: 1.75, 95% CI: 1.24–2.46; additive model: p = 0.003, OR: 1.51, 95% CI: 1.15–1.98). Meanwhile, in Table 6, we found that rs420137, rs386360, and rs7763726 were protective factors against the occurrence of hypertension in CHD patients (rs420137: allele model: p = 0.004, OR: 0.70, 95% CI: 0.55–0.89; co-dominant model with the “GC” genotype: p < 0.001, OR: 0.53, 95% CI: 0.38–0.76; dominant model: p < 0.001, OR: 0.54, 95% CI: 0.38–0.76; additive model: p = 0.002, OR: 0.68, 95% CI: 0.53–0.87; rs386360: allele model: p = 0.005, OR: 0.71, 95% CI: 0.55–0.90; co-dominant model: p = 0.003, OR: 0.59, 95% CI: 0.41–0.83; dominant model: p = 0.001, OR: 0.58, 95% CI: 0.41–0.81; additive model: p = 0.002, OR: 0.67, 95% CI: 0.52–0.87; rs7763726: allele model: p = 0.004, OR: 0.69, 95% CI: 0.54–0.89; co-dominant model: p = 0.001, OR: 0.56, 95% CI: 0.40–0.79; dominant model: p = 0.001, OR: 0.57, 95% CI: 0.41–0.79; additive model: p = 0.002, OR: 0.67, 95% CI: 0.51–0.87). Nevertheless, after stratified by diabetes, no other SNP except rs7763726 was found in the recessive model (p = 0.011, OR: 0.29, 95% CI: 0.11–0.75) to reduce the risk of diabetes in CHD individuals (Table 6).

Table 5

SNP-IDModelGenotypeSmokersNon-smokersDrinkersNon-drinkers
OR (95% CI)POR (95% CI)POR (95% CI)POR (95% CI)P
rs420137AlleleC1111
G1.43 (1.10–1.84)0.0061.12 (0.12–0.88)0.3531.30 (1.01–1.69)0.0441.20 (0.95–1.52)0.133
Co-dominantCC1111
GC1.38 (0.97–1.95)0.0721.18 (0.84–1.65)0.3331.19 (0.84–1.69)0.3361.37 (0.99–1.91)0.060
GG2.24 (1.18–4.25)1.13 (0.64–2.02)1.79 (0.98–3.26)1.28 (0.70–2.28)
DominantCC1111
GC-GG1.50 (1.07–2.09)0.0171.17 (0.85–1.61)0.3331.29 (0.92–1.90)0.1391.35 (0.99–1.86)0.060
RecessiveCC-GC1111
GG1.94 (1.04–3.61)0.0361.04 (0.60–1.81)0.8831.66 (0.93–2.95)0.0861.08 (0.61–1.91)0.785
Additive1.44 (1.11–1.88)0.0061.11 (0.86–1.42)0.4161.28 (0.99–1.65)0.0581.23 (0.95–1.57)0.111
rs386360AlleleA1111
C1.46 (0.13–1.13)0.0041.11 (0.12–0.87)0.4171.27 (0.13–0.97)0.0841.25 (0.12–0.98)0.069
Co-dominantA/A1111
C/A1.43 (1.00–2.02)0.0451.19 (0.85–1.66)0.3201.15 (0.81–1.64)0.4281.45 (1.04–2.01)0.027
C/C2.44 (1.22–4.87)1.08 (0.60–1.97)1.76 (0.93–3.34)1.31 (0.71–2.40)
DominantA/A1111
C/A-C/C1.55 (1.11–2.16)0.0101.17 (0.85–1.61)0.3421.24 (0.89–1.73)0.2111.43 (1.04–1.95)0.028
RecessiveA/A-C/A1111
C/C2.08 (1.06–4.08)0.0320.99 (0.56–1.78)0.9971.65 (0.89–3.07)0.1121.09 (0.61–1.97)0.764
Additive1.49 (1.14–1.95)0.0031.10 (0.86–1.41)0.4591.25 (0.96–1.62)0.0971.27 (0.99–1.64)0.061
rs7763726AlleleA1111
G1.83 (0.14–1.38)< 0.0010.53 (0.13–0.85)0.4781.56 (1.18–2.06)0.0021.22 (0.95–1.57)0.126
Co-dominantA/A1111
A/G2.13 (1.48–3.05)< 0.0011.11 (0.80–1.56)0.5281.75 (1.22–2.51)0.0021.28 (0.93–1.78)0.133
G/G2.66 (1.19–5.98)1.12 (0.59–2.12)1.76 (0.90–3.43)1.43 (0.69–2.94)
DominantA/A1111
G/A-G/G2.19 (1.55–3.10)< 0.0011.12 (0.81–1.54)0.5071.75 (1.24–2.46)0.0011.30 (0.95–1.78)0.104
RecessiveA/A-G/A1111
G/G2.02 (0.91–4.47)0.0851.07 (0.57–1.99)0.8381.42 (0.74–2.74)0.2881.28 (0.63–2.60)0.493
Additive1.91 (1.42–2.57)< 0.0011.08 (0.84–1.40)0.5391.51 (1.15–1.98)0.0031.24 (0.95–1.62)0.106

FNDC1 SNPs associated with CHD risk stratified by smoking and drinking statues.

CHD, coronary heart disease; SNP, single nucleotide polymorphism; OR, odds ratio; CI, confidence interval. OR (95% CI) values were calculated by logistic regression analysis with adjustments for gender. Bold values are statistically significant (p < 0.05).

Table 6

SNP-IDModelGenotypeHypertension vs. non-hypertensionDiabetes vs. non-diabetes
OR (95% CI)POR (95% CI)P
rs420137AlleleC11
G0.70 (0.55–0.89)0.0040.97 (0.75–1.26)0.806
Co-dominantCC11
GC0.53 (0.37–0.75)<0.0011.17 (0.81–1.68)0.405
GG0.58 (0.33–1.03)0.72 (0.37–1.37)
DominantCC11
GC-GG0.54 (0.38–0.76)<0.0011.08 (0.76–1.53)0.687
RecessiveCC-GC11
GG0.82 (0.48–1.38)0.4500.66 (0.36–1.23)0.191
Additive0.68 (0.53–0.87)0.0020.96 (0.74–1.25)0.760
rs386360AlleleA11
C0.71 (0.55–0.90)0.0050.96 (0.74–1.26)0.789
Co-dominantA/A11
C/A0.59 (0.41–0.83)0.0031.22 (0.85–1.76)0.272
C/C0.53 (0.30–0.96)0.62 (0.30–1.25)
DominantA/A11
C/A-C/C0.58 (0.41–0.81)0.0010.11 (0.78–1.57)0.575
RecessiveA/A-C/A11
C/C0.70 (0.40–1.21)0.2040.56 (0.28–1.10)0.092
Additive0.67 (0.52–0.87)0.0020.96 (0.73–1.25)0.740
rs7763726AlleleA11
G0.69 (0.54–0.89)0.0040.89 (0.67–1.17)0.394
Co-dominantA/A11
A/G0.56 (0.40–0.79)0.0011.21 (1.85–1.73)0.293
G/G0.59 (0.31–1.12)0.32 (0.12–0.83)
DominantA/A11
G/A-G/G0.57 (0.41–0.79)0.0011.05 (0.74–1.48)0.796
RecessiveA/A-G/A11
G/G0.77 (0.41–1.43)0.4070.29 (0.11–0.75)0.011
Additive0.67 (0.51–0.87)0.0020.88 (0.66–1.16)0.364

FNDC1 SNPs associated with CHD risk stratified by hypertension vs. non-hypertension and diabetes vs. non-diabetes.

CHD, coronary heart disease; SNP, single nucleotide polymorphism; OR, odds ratio; CI, confidence interval. OR (95% CI) values were calculated by logistic regression analysis with adjustments for age, gender, smoking, and drinking. Bold values are statistically significant (p < 0.05).

Haplotype and MDR analyses

Haplotype analysis was performed to analyze the association of FNDC1 polymorphisms with CHD risk, and one block including rs420137, rs386360 and rs7763726 was noted (Figure 1). Meanwhile, the connection between FNDC1 haplotypes and the risk of CHD is showed in Table 7. We found that the haplotype “GCG” (p = 0.002, OR: 1.34, 95% CI: 1.11–1.63) and the haplotype “CAA” (p = 0.004, OR: 1.29, 95% CI: 1.09–1.55) had a higher risk of CHD after adjustment.

Figure 1

Figure 1

Haplotype analysis was performed to analyze the association of FNDC1.

Table 7

SNPsHaplotypeF-AF-UORPa
rs420137|rs386360|rs7763726GCG0.2710.2711.34 (1.11–1.63)0.002
GCA0.0420.0500.82 (0.56–1.21)0.318
GAA0.0170.0141.20 (0.62–2.33)0.588
CAA0.3400.2871.29 (1.09–1.55)0.004

Haplotype analysis.

SNP, single nucleotide polymorphism; F-A, frequency in cases; F-U, frequency in controls; OR, odds ratio; 95% CI, 95% confidence interval.

pa

values were calculated using Pearson's chi-squared test. Bold values are statistically significant (p < 0.05).

MDR analysis showed that the best three-locus model was rs420137, rs386360, and rs7763726 (Bal. Acc. CV Training: 0.552; Bal. Acc. CV Testing: 0.548; CVC: 10/10; Testing Sensitivity: 0.647; Testing Specificity: 0.397; p = 0.0004, OR: 1.52; 95% CI: 1.20–1.92) (Table 8).

Table 8

ModelTraining Bal. Acc.Testing Bal. Acc.CVCOR (95% CI)P
rs77637260.5510.55110/101.51 (1.20–1.92)0.0005
rs420137,rs77637260.5510.5509/101.51 (1.20–1.92)0.0005
rs420137,rs386360,rs77637260.5520.54810/101.52 (1.20–1.92)0.0004

Analysis of SNP-SNP interaction models using MDR method.

MDR, multifactor dimensionality reduction; Bal. Acc., balanced accuracy; CVC, cross-validation consistency; OR, odds ratio; 95% CI, 95% confidence interval. P-values were calculated using χ2 tests. Bold values indicate statistical significance (p < 0.05).

Discussion

CHD is a disease affected by many factors. Diverse environmental and genetic factors have been reported to be associated with the etiology of CHD (14). So far, scores of genetic studies have revealed that genes are related to the pathogenesis of CHD (15). It has been reported that FNDC1 can activate the cell surface receptor AGS8, regulate the intracellular distribution of vascular endothelial growth factor receptor 2 (VEGFR2), promote VEGF-A-induced vascular endothelial cell signal transduction and growth, promote tumor angiogenesis, and lead to increased tumor stage (16). Hence, we investigated the relationship between FNDC1 polymorphisms and CHD susceptibility by a case-control study. Our results indicated that FNDC1-rs420137, -rs386360, and -rs7763726 played an important role in enhancing the risk of CHD. Subgroup analysis showed that rs420137 increased the susceptibility to CHD in males, patients aged ≤62 years, and smokers. Rs386360 had a higher risk of CHD in males, patients aged ≤62 years, smokers, and non-drinkers. Moreover, the association of rs7763726 with increased CHD susceptibility was also revealed in males, patients aged ≤62 years, smokers, and drinkers. Last but not least, these three SNPs we selected were protective factors against hypertension in CHD individuals.

FNDC1 contains a major component of the fibronectin domain (17). Recent studies have shown that FNDC1 is closely linked to the occurrence of multiple illnesses (8, 1821). It was also found that FNDC1 took part in VEGF-mediated endothelial cell angiogenesis, especially the distribution and transportation of VEGFR2, which had a strong tyrosine kinase activity and was a firsthand factor influencing angiogenesis (19, 22). The existing studies have also displayed that FNDC1 is a pathogenic gene of acute otitis media in children, and has a certain correlation with hypertension (8). In one study, the FNDC1 polymorphism (rs3003174) was associated with coronary artery aneurysm complications in Kawasaki disease (23). Previous study found that the expression of FNDC1 was up-regulated in cardiomyocytes induced by ischemia or hypoxia (24). Motohiko Sato et al. found FNDC1 isolated from rat hearts with repeated transient ischemia of collateral circulation may induce cardiomyocyte apoptosis under inflammatory or hypoxic conditions. Hypoxia induced FNDC1 mRNA expression in rat aortic smooth muscle cells, endothelial cells, or myocardial fibroblasts, suggesting a cardiomyocyte-specific adaptive mechanism including G protein signaling pathway remodeling (25). Bouchareb et al. (26) found that the mutations in the FNDC1 gene were associated with cardiovascular disease, indicating that FNDC1 polymorphisms may be involved in the development of CHD.

Sex and age are strong determinants of CHD risk (27), so is smoking (28). According to sex-stratified analysis, we noticed that the three SNPs in this study affected CHD susceptibility only in males, suggesting that there were gender differences in the influence of SNPs on the risk of CHD. At the same time, these three SNPs increased the risk of CHD in patients aged 62 years and younger and smokers. These results are consistent with previous studies indicating that genetic polymorphisms influence the susceptibility to CHD in males or patients aged ≤61 years (29, 30). Zhao et al. (31) found that the morbidity and mortality of cardiovascular diseases were correlated with the degree of blood pressure elevation. In our study, we observed that rs420137, rs386360 and rs7763726 could protect individuals with CHD from hypertension. These results indicate that patients with CHD need individualized treatment.

We need to emphasize that our study still has some limitations. First, our subjects were all Han Chinese. Second, the sample size was too small to validate the results. Third, CHD is a disease affected by a great number of risk factors. Therefore, we cannot completely exclude the potential impact of the other factors on its development.

Conclusions

In a word, this study indicated that FNDC1 polymorphisms (rs420137, rs386360, and rs7763726) were related to CHD risk, suggesting that these three variants may be new biomarkers for predicting CHD risk. To the best of our knowledge, this is the first study on the relationship between FNDC1 gene polymorphisms and CHD, and the specific mechanism of CHD still needs to be further studied.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.

Ethics statement

This study was reviewed and approved by the Ethics Committee of the First Hospital of Yulin City, and conducted in accordance with the standards of the Declaration of Helsinki. Written informed consent was obtained from each participant included in the study.

Author contributions

HM contributed to the study design and performed statistical analysis. XH was responsible for manuscript preparation. XL and XD performed experiments. JH was responsible for checking the data. HZ and YZ contributed to data collection. All authors were responsible for drafting the manuscript, they have read, and approved the final version.

Acknowledgments

We thank all individuals for their participation. We are also grateful to the hospital staff for collecting blood samples and clinical information.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1.

    MathersCDLoncarD. Projections of global mortality, and burden of disease from 2002 to 2030. PLoS Med. (2006) 3:201130. 10.1371/journal.pmed.0030442

  • 2.

    ZhangXHLuZLLiuL. Coronary heart disease in China. Heart. (2008) 94:112631. 10.1136/hrt.2007.132423

  • 3.

    FiatalSAdanyR. Application of single-nucleotide polymorphism-related risk estimates in identification of increased genetic susceptibility to cardiovascular diseases: a literature review. Front Public Health. (2017) 5:358. 10.3389/fpubh.2017.00358

  • 4.

    McPhersonRTybjaerg-HansenA. Genetics of coronary artery disease. Circ Res. (2016) 118:56478. 10.1161/CIRCRESAHA.115.306566

  • 5.

    LiQDuQ. Associations between nine candidate genetic polymorphisms with coronary heart disease: a meta-analysis. Herz. (2020) 45:1528. 10.1007/s00059-019-4806-7

  • 6.

    VogelLKSaeboMHoyerHKoppTIVogelUGodiksenSet al. Intestinal PTGS2 mRNA levels, PTGS2 gene polymorphisms, and colorectal carcinogenesis. PLoS ONE. (2014) 9:e105254. 10.1371/journal.pone.0105254

  • 7.

    HolanVTrosanPCejkaCJavorkovaEZajicovaAHermankovaBet al. A comparative study of the therapeutic potential of mesenchymal stem cells and limbal epithelial stem cells for ocular surface reconstruction. Stem Cells Transl Med. (2015) 4:105263. 10.5966/sctm.2015-0039

  • 8.

    van IngenGLiJGoedegebureAPandeyRLiYRMarchMEet al. Genome-wide association study for acute otitis media in children identifies FNDC1 as disease contributing gene. Nat Commun. (2016) 7:12792. 10.1038/ncomms12792

  • 9.

    LapointeJLiCHigginsJPvan de RijnMBairEMontgomeryKet al. Gene expression profiling identifies clinically relevant subtypes of prostate cancer. Proc Natl Acad Sci USA. (2004) 101:8116. 10.1073/pnas.0304146101

  • 10.

    ToWSMidwoodKS. Plasma and cellular fibronectin distinct and independent functions during tissue repair. Fibrogenesis Tissue Repair. (2011) 4:117. 10.1186/1755-1536-4-21

  • 11.

    WilliamsLVVeliceasaDVinokourEVolpertOV. miR-200b inhibits prostate cancer EMT, growth and metastasis. PLoS ONE. (2013) 8:e83991. 10.1371/journal.pone.0083991

  • 12.

    HalperJKjaerM. Basic components of connective tissues and extracellular matrix: elastin, fibrillin, fibulins, fibrinogen, fibronectin, laminin, tenascins and thrombospondins. Adv Exp Med Biol. (2014) 802:3147. 10.1007/978-94-007-7893-1_3

  • 13.

    DengAYChauvetCMenardA. Alterations in fibronectin type III domain containing 1 protein gene are associated with hypertension. PLoS ONE. (2016) 11:e0151399. 10.1371/journal.pone.0151399

  • 14.

    MussattoKAHoffmannRHoffmanGTweddellJSBearLCaoYet al. Risk factors for abnormal developmental trajectories in young children with congenital heart disease. Circulation. (2015) 132:75561. 10.1161/CIRCULATIONAHA.114.014521

  • 15.

    HeB-XLeiShi QiuJTaoLLiRYangLZhaoS-Jet al. Functional polymorphism in the CYP3A4 gene is associated with increased risk of coronary heart disease in the chinese han population. Basic Clin Pharmacol Toxicol. (2010) 108:20813. 10.1111/j.1742-7843.2010.00657.x

  • 16.

    SakimaMHayashiHMamunAASatoM. VEGFR-3 signaling is regulated by a G-protein activator, activator of G-protein signaling 8, in lymphatic endothelial cells. Exp Cell Res. (2018) 368:1323. 10.1016/j.yexcr.2018.04.007

  • 17.

    GaoMCraigDLequinOCampbellIDVogelVSchultenK. Structure and functional significance of mechanically unfolded fibronectin type III1 intermediates. Proc Natl Acad Sci USA. (2003) 100:147849. 10.1073/pnas.2334390100

  • 18.

    AndereggUBreitschwerdtKKöhlerMJSticherlingMHausteinU-FSimonJCet al. MEL4B3, a novel mRNA is induced in skin tumors and regulated by TGF-beta and pro-inflammatory cytokines. Exp Dermatol. (2005) 14:70918. 10.1111/j.0906-6705.2005.00349.x

  • 19.

    HayashiHMamunAASakimaMSatoM. Activator of G-protein signaling 8 is involved in VEGF-mediated signal processing in angiogenesis. J Cell Sci. (2016) 129:121022. 10.1242/jcs.181883

  • 20.

    SatoMJiaoQHondaTKurotaniRToyotaEOkumuraSet al. Activator of G protein signaling 8 (AGS8) is required for hypoxia-induced apoptosis of cardiomyocytes: role of G betagamma and connexin 43 (CX43). J Biol Chem. (2009) 284:3143140. 10.1074/jbc.M109.014068

  • 21.

    WuenschTWizentyJQuintJSpitzWBosmaMBeckerOet al. Expression analysis of fibronectin type III domain-containing (FNDC) genes in inflammatory bowel disease and colorectal cancer. Gastroenterol Res Pract. (2019) 2019:3784172. 10.1155/2019/3784172

  • 22.

    ShibuyaM. Vascular endothelial growth factor (VEGF) and its receptor (VEGFR) signaling in angiogenesis: a crucial target for anti- and pro-angiogenic therapies. Genes Cancer. (2011) 2:1097105. 10.1177/1947601911423031

  • 23.

    LinKZhangLWangYLiJXuYCheDet al. FNDC1 polymorphism (rs3003174 C > T) increased the incidence of coronary artery aneurysm in patients with kawasaki disease in a southern chinese population. J Inflamm Res. (2021) 14:263340. 10.2147/jir.S311956

  • 24.

    SatoMCismowskiMJToyotaESmrckaAVLucchesiPAChilianWMet al. Identification of a receptor-independent activator of G protein signaling (AGS8) in ischemic heart and its interaction with Gbetagamma. Proc Natl Acad Sci USA. (2006) 103:797802. 10.1073/pnas.0507467103

  • 25.

    YuanCSatoMLanierSMSmrckaAV. Signaling by a non-dissociated complex of G protein βγ and α subunits stimulated by a receptor-independent activator of G protein signaling, AGS8. J Biol Chem. (2007) 282:1993847. 10.1074/jbc.M700396200

  • 26.

    BoucharebRGuauque-OlarteSSniderJZaminskiDAnyanwuAStelzerPet al. Proteomic architecture of valvular extracellular matrix: FNDC1 and MXRA5 are new biomarkers of aortic stenosis. JACC Basic Transl Sci. (2021) 6:2539. 10.1016/j.jacbts.2020.11.008

  • 27.

    Barrett-ConnorE. Gender differences and disparities in all-cause and coronary heart disease mortality: epidemiological aspects. Best Pract Res Clin Endocrinol Metab. (2013) 27:481500. 10.1016/j.beem.2013.05.013

  • 28.

    HajarR. Risk factors for coronary artery disease: historical perspectives. Heart Views Offic J Gulf Heart Assoc. (2017) 18:10914. 10.4103/HEARTVIEWS.HEARTVIEWS_106_17

  • 29.

    ChenXZhangZWangXChenYWangC. NT5C2 gene polymorphisms and the risk of coronary heart disease. Public Health Genom. (2020) 23:909. 10.1159/000507714

  • 30.

    MaRZhuXBoY. SCARB1 rs5888 gene polymorphisms in coronary heart disease: a systematic review and a meta-analysis. Gene. (2018) 678:S0378111918308862-. 10.1016/j.gene.2018.08.024

  • 31.

    ZhaoXCYangSHYanYQXinZLinZBoJet al. Identification of differential gene expression profile from peripheral blood cells of military pilots with hypertension by RNA sequencing analysis. BMC Med Genom. (2018) 11:59. 10.1186/s12920-018-0378-2

Summary

Keywords

FNDC1, CHD, susceptibility, genetic polymorphism, SNPs

Citation

He X, Li X, Du X, Han J, Zhang H, Zhu Y and Ma H (2022) Rs420137, rs386360 and rs7763726 polymorphisms in fibronectin type III domain containing 1 are associated with susceptibility to coronary heart disease: Analysis in the Han population. Front. Cardiovasc. Med. 9:964978. doi: 10.3389/fcvm.2022.964978

Received

09 June 2022

Accepted

20 September 2022

Published

06 October 2022

Volume

9 - 2022

Edited by

Olorunseun O. Ogunwobi, Hunter College (CUNY), United States

Reviewed by

Priyanka Ghosh, The City University of New York, United States; Amrita Bhat, University of Jammu, India

Updates

Copyright

*Correspondence: Honghong Ma

†These authors share first authorship

This article was submitted to Atherosclerosis and Vascular Medicine, a section of the journal Frontiers in Cardiovascular Medicine

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics