Skip to main content

Interaction study of the effects of environmental exposure and gene polymorphisms of inflammatory and immune-active factors on chronic obstructive pulmonary disease

Abstract

Background

Chronic obstructive pulmonary disease (COPD) is a heterogeneous disease, influenced by both environmental and genetic factors. Single nucleotide polymorphism (SNP) in the human genome may influence the risk of developing COPD and the response to treatment. We assessed the effects of gene polymorphism of inflammatory and immune-active factors and gene-environment interaction on risk of COPD in middle-aged and older Chinese individuals.

Methods

In this community-based case–control study, 471 patients with COPD and 485 controls aged 40–76 years in Heilongjiang Province, China were included. Face-to-face interviews, lung function tests, and multiplex polymerase chain reaction were used to obtain data. Logistic regression model, generalized multifactor dimensionality reduction and crossover analysis were used to analyse the effects of SNPs, gene–gene interactions, and gene-environment interactions on COPD.

Results

CRP gene[rs1130864-A allele (OR, 1.77; 95% CI 1.11–2.81); G/A + A/A genotype (OR, 1.75; 95% CI 1.07–2.84)], FCAR gene[rs4806606-G (OR, 0.72; 95% CI 0.53–0.98); rs8112766-G (OR, 0.79; 95% CI 0.64–0.98)] and FCGR2A gene[rs4656308-C (OR, 0.74; 95% CI 0.55–1.00); rs4656309-T (OR, 0.81; 95% CI 0.66–0.99)] are independent influential factors for COPD. Rs1205 [RERI: 0.15 (0.07–1.00)] and rs1130864 [RERI: 2.45 (0.73–4.18)] of CRP gene, rs11084376 [OR: 0.54 (0.29–0.97)] of FCAR gene, rs844 of FCGR2B [SI: 0.30 (0.11–0.77); OR: 0.46 (0.24–0.90)] gene, rs4656308-rs4656309-rs2165088 haplotype [SI: 0.48 (0.26–0.89)] of FCGR2A gene and exposure to smoking index > 200, indoor coal/wood/straw use, and outdoor straw burning play synergistic or antagonistic roles in the development of COPD.

Conclusions

Alleles and genotypes of the CRP/FCAR/FCGR2A gene can increase the susceptibility to COPD in the northern Chinese population. For the first time, environmental exposure to the CRP/FCAR/FCGR2A/FCGR2B genes has been shown to have synergistic or antagonistic effects on COPD susceptibility on genotypes or haplotypes.

Introduction

Chronic obstructive pulmonary disease (COPD) is defined as a pulmonary or systemic inflammatory disease caused by harmful gases or particles and is characterized by progressive and irreversible airflow restriction, which is preventable and treatable [1]. Environmental exposure, represented by smoking, is considered to be the most common risk factor for COPD, and biofuels are an important cause of COPD in developing countries [2]. Our research group previously identified the smoking index, outdoor straw burning, and indoor use of coal/wood/straw/mixed fuel as risk factors for COPD in Heilongjiang Province, China [3]. However, not all smokers and individuals with a history of air pollution exposure develop COPD. Studies have shown that individuals with a family history of COPD, particularly those with parents who have the disease, are at an increased risk of developing COPD. Alpha one-antitrypsin deficiency has been identified as a genetic factor that contributes to COPD [4]. With the development of molecular biology and GWAS technology, genetic factors have gained increasing attention, although research on other genetic factors related to COPD still yields inconsistent results.

At present, studies have revealed that single nucleotide polymorphisms (SNPs) in the human genome play crucial roles in the occurrence of various diseases, clinical manifestations, and response to drug treatment. COPD is a heterogeneous disease influenced by both environmental and genetic factors. Extensive research has been conducted on the pathogenic mechanisms of COPD-related genes. In recent years, domestic and international scholars have discovered new susceptibility genes related to COPD through genetic polymorphism association analysis, whole-genome scanning, and other methods. Several studies have revealed significant associations between SNP sites such as serine peptidase inhibitor 2 (SERPINE2), hedgehog interacting protein (HHIP), β2-adrenergic receptor (ADRB2), vascular endothelial growth factor (VEGFA), vitamin D binding protein (VDBP), heat shock protein 70 (Hsp70), serine and arginine-rich with an additional glutamic acid-lysine-rich 1 (SREK1), Toll-like receptor (TLR), and human leukocyte antigen (HLA) [5,6,7,8,9,10,11], and COPD in different regions and ethnic groups. However, there are inconsistencies and even contradictions in the association studies of these candidate genes with COPD.

The main feature of COPD is chronic, progressive inflammation of the lung tissue. However, systemic immune and inflammatory responses are also important features of COPD. Our research group revealed that elevated C-reactive protein (CRP) levels are a risk factor for COPD; thus, it is necessary to conduct in-depth research on the associations between CRP expression-related genes and COPD. Fadi G HAge et al. reported that serum CRP levels are influenced by genetic factors. In family and twin studies, the heritability of elevated CRP ranges from 39 to 52% [12]. Wener MH et al. reported that the CRP levels of African descendants in the United States were higher than those of European descent. Aroon D Hingorani et al. reported that certain subtle changes in the CRP gene sequence, mainly SNPs, can affect blood CRP levels. By exploring the function of SNPs in the CRP gene in vitro, Carlson CS et al. reported that a variety of SNPs and haploids constructed in the CRP promoter region are closely related to the CRP level [13].

Our research group also reported that immunoglobulin G (IgG) and immunoglobulin A (IgA) are associated with the development of COPD. Both the Fc fragment of the IgG receptor (FCGR) and the Fc fragment of the IgA receptor (FCAR) genes are members of the immunoglobulin gene superfamily and are associated with the Fc receptors for IgA and IgE, which are also the Fc receptors for CRP [14]. The Fc fragment of IgG receptor IIa (FCGR2A) and the Fc fragment of IgG receptor IIb (FCGR2B) are expressed on the surface of monocytes, neutrophils, eosinophils, dendritic cells and macrophages and are the two receptors that cover most white blood cells. FCGR2B transmits inhibitory signals to play an anti-inflammatory role, while all other receptors trigger cell activation to play a proinflammatory role, and the combination of IgG and FCGR can play both proinflammatory and anti-inflammatory roles. Previous studies have shown that the genetic polymorphism of FCGR in IgG is associated with the risk and disease progression of systemic lupus erythematosus [15]. Through competitive binding and mutation analysis, the FCAR, the main receptor of IgA, was determined to be the receptor of pentraxin. The binding site of the FCAR identified by CRP is different from that identified by IgA. Therefore, CRP and IgA do not compete with each other for binding sites and receptors. Moreover, CRP promotes the ability of IgA to recognize the FCAR [16]. Simultaneous recognition of different FCAR regions by CRP and IgA can increase FCAR activation, improving the effect of the FCAR on the innate immune response.

Previous studies have shown that SNPs in the CRP gene are closely related to the occurrence of a variety of diseases, including cardiovascular diseases, autoimmune diseases and tumors. Immunoglobulin-related gene SNPs are also closely related to the occurrence of many diseases, including autoimmune diseases, infectious diseases and tumors. However, it is not clear whether CRP, IgG and IgA gene-associated SNPs are related to the development of COPD. Moreover, there is a lack of relevant studies on the impact of environmental factors and interactions between gene SNPs or gene haplotypes on COPD. Therefore, the effects of alleles, genotypes, haplotypes, gene‒gene interactions and environment‒gene interactions of CRP/FCAR/FCGR2A/FCGR2B polymorphic sites on COPD development need to be further studied.

Material and methods

Selection of cases and controls

The selection of study subjects, questionnaire surveys, physical examinations, serological testing and diagnostic criteria for COPD are consistent with the findings of a previously published study by our group on the effects of environmental, immune, and inflammatory risk factors on COPD [3]. The baseline survey for this prospective cohort study was conducted from November 2018 to September 2019, and the study subjects were selected from the community population in Mingshui County, Suihua, Heilongjiang Province, China. The study included 471 cases and 485 cases in the control group, with a total of 956 participants ranging in age from 40 to 76 years. Patients with chronic obstructive pulmonary disease were defined as individuals with a forced expiratory volume in one second (FEV1) to forced vital capacity (FVC) ratio of less than 70% (FEV1/FVC < 70%) to represent the definition of an epidemiological case; the remaining patients were included as eligible controls according to the inclusion and exclusion criteria. The investigation was conducted with the approval of the Ethics Committee of the Center for Endemic Disease Control of the Chinese Center for Disease Control and Prevention. All participants signed written informed consent.

Questionnaires and genotyping

The data collected through the questionnaire survey regarding demographics, disease history, smoking status, indoor fuel use, and outdoor straw burning information and the classification of the study subjects were in line with previously published research [3]. By referencing the 1000 Genomes Project data (https://www.internationalgenome.org) and using Haploview software (version 4.2), we employed a criterion of minor allele frequency (MAF) ≥ 5% for screening target gene SNPs and a linkage disequilibrium (LD) threshold of r2 ≥ 0.8 for the haplotypes in the Chinese population. Additionally, we referred to previous relevant studies to select SNP loci associated with inflammation and immune-related diseases. Genomic DNA (Deoxyribonucleic acid) was extracted from the blood via a DNA extraction kit (MagPure Blood DNA KF Kit), and the DNA concentration was quantified via a NanoDrop One microspectrophotometer. Genotype sequencing of DNA samples whose concentrations were above 20 µg/ml was completed in the laboratory of Genesky Biotechnologies (Shanghai, China), and SNP genotyping was performed via multiplex polymerase chain reaction (PCR).

Statistical analyses

All the data were gathered via dedicated software and analysed via the SAS statistical analysis system (ver 9.1.4), R software (ver 3.6.3) and Excel software. First, univariate analysis, including the Wilcoxon rank sum test or t test for continuous data and the chi-square test for classification variables, was performed to identify factors significantly associated with COPD. The chi-square test was used to assess the Hardy‒Weinberg equilibrium (HWE) of SNPs in control participants. The chi-square test and Fisher's exact test were applied to analyse the differences in allele frequency and genotype frequency in SNPs of each target gene. SNPStats software was used to adjust for the influence of confounding factors such as age, body mass index (BMI) and education level and analyse the genetic model of the gene SNPs associated with the risk of COPD. A logistic regression model was used to analyse the haplotypes associated with the risk of COPD after adjusting for confounding factors such as age, BMI and education level. Haploview software (version 4.2) was used to construct LD and haplotypes, and the differences in haplotype frequency distributions between groups were compared. Generalized multifactor dimensionality reducation (GMDR) software (ver 0.9) was used to analyse the higher-order interactions between the SNPs and the haploids.

We evaluated the impact of gene‒environment interactions on COPD by using odds ratios (ORs). In addition, the relative excess risk due to interaction (RERI) and synergy index (SI) were applied to evaluate the additive interaction effect between candidate gene SNPs and environmental exposure on COPD. The evaluation criteria were as follows: (1) if the RERI values and the upper and lower limits of the 95% CI are both greater than 0, a positive additive interaction exists; if both the RERI value and the upper and lower limits of the 95% CI are less than 0, then a negative additive interaction exists; otherwise, no additive interaction exists. (2) If the SI value and the upper and lower limits of the 95% CI are both greater than 1, then there is a positive additive interaction; if the SI value and the upper and lower limits of the 95% CI are both less than 1, then there is a negative additive interaction; otherwise, there is no additive interaction. A two-sided probability test with P < 0.05 was considered statistically significant.

Results

Study population characteristics and major risk factors

In our previous study [3], we compared the characteristics of 471 cases and 485 controls with respect to the main demographic characteristics, disease history, environmental exposure, lifestyle, dietary habits and disease history. Statistical differences were found in age, income, education, BMI, smoking, indoor fuel use (no fuel/natural gas/electricity vs. coal/wood/straw), outdoor straw burning, history of respiratory disease (RD), history of coronary heart disease and congestive heart failure (CHD/CHF) and history of cerebrovascular disease (CVD) between the two groups (P < 0.05) [3]. The correlations between genetics or the genetic or genetic-environmental factors and COPD were analysed after adjustment for the above covariates in the present study.

A HWE test of the SNP loci of the candidate genes

A chi-square test (Table S1 in the Supplementary Appendix) revealed that five SNPs (rs1341665, rs1865096, rs438228, rs56159502 and rs1771576) failed the HWE test (P < 0.05) in the controls and were excluded from the analysis. rs1800947 and rs16986050 were also excluded because the MAF was < 5% in the study population. Eventually, 18 SNPs, including rs1130864, rs1205, rs2808630, rs4806606, rs4806608, rs62123424, rs7259090, rs8112766, rs11084376, rs1865097, rs2304225, rs4656308, rs4656309, rs4657040, rs1801274, rs2165088, rs844 and rs16827592, were included in the subsequent analysis.

Association analysis of SNP genotypes

There were significant differences in the compositions of alleles and genotypes between the case and control groups (Table 1 and Figure S1 in the Supplementary Appendix). The rs1130864-A allele in the CRP gene was positively associated with susceptibility to COPD (OR, 1.77; 95% CI 1.11–2.81); the rs4806606-G and rs8112766-G alleles in the FCAR gene and the rs4656308-C and rs4656309-T alleles in the FCGR2A gene were negatively associated with susceptibility to COPD [OR, 0.72; 95% CI 0.53–0.98; OR, 0.79; 95% CI 0.64–0.98; OR, 0.74; 95% CI 0.55–1.00; OR, 0.81; 95% CI 0.66–0.99].

Table 1 Allele frequencies in COPD and control subjects for SNPs

Remarkable differences were found in the compositions of the rs1130864 genotypes in the CRP gene between the case and control groups (Table S2 in the Supplementary Appendix). In the dominant model, individuals with at least one mutation allele (G/A + A/A genotypes) had a greater risk of COPD than did wild-type homozygotes (G/G genotype) (OR, 1.75; 95% CI 1.07–2.84).

Owing to the limited frequency observed in the case and control groups, specifically with a count of three in each for the A/A genotype of the CRP gene in the codominant, recessive, and log-additive models, the reliability of the results may be compromised. Therefore, it is advisable to exclude these three models from consideration (Tables S3-S6 in the Supplementary Appendix).

Association analysis of haplotypes

Linkage disequilibrium (LD) analysis and haplotype construction were performed for each SNP site, and D’ > 0.8 was used as the criterion for meeting the strong linkage relationship. The following haplotypes were constructed: strong linkage disequilibrium was detected between rs2808630 and rs1205 (r2 = 0.264, D' > 0.8), and the haplotype block rs2808630-rs1205 in the CRP gene was constructed. Three haplotypes were identified. A second strong linkage disequilibrium was detected among rs2304225, rs11084376 and rs8112766 in the FCAR gene, and four haplotypes were identified. In addition, strong LD was found for rs1865097 and rs62123424 in the FCAR gene, and three haplotypes were identified. Additionally, a fourth strong linkage disequilibrium was constructed among rs4656308, rs4656309 and rs2165088 in the FCGR2A gene, and five haplotypes were identified. Nonetheless, no strong linkage disequilibrium was found between rs1801274 and rs16827592 in the FCGR2B gene (Fig. 1).

Fig. 1
figure 1

Linkage disequilibrium plot for SNPs of genes. (Shades in bright red indicate LOD ≥ 2 and D’ = 1. Shades in white indicate LOD < 2 and D’ < 1. The number in the block represents the value of D’)

The results revealed that none of the haplotypes were associated with COPD susceptibility (P > 0.05) (Table S7 in the Supplementary Appendix).

Association analysis of GMDR

The GMDR results for the gene SNP‒gene SNP interaction models revealed that the three-locus model (rs1205, rs8112766, and rs844) had the highest test accuracy (Table S8 in the Supplementary Appendix and Fig. 2). Therefore, the three-locus model composed of rs1205, rs8112766 and rs844 was chosen as the best model for predicting COPD risk, with a test accuracy of 0.5147 and a perfect CVC = 6/10.

Fig. 2
figure 2

Optimum model of the GMDR analysis results

The GMDR results for haplotype‒haplotype interaction models for COPD indicated that the three loci model (rs2808630‒rs1205, rs2304225‒rs11084376‒rs8112766, and rs1865097‒rs62123424) has the highest test accuracy (Table S9 in the Supplementary Appendix and Fig. 2). Thus, the three loci model composed of rs2808630-rs1205, rs2304225-rs11084376-rs8112766 and rs1865097-rs62123424 was identified as the best model for predicting COPD risk, with a test accuracy of 0.5315 and a perfect CVC = 5/10.

However, high-order interactions with significant differences between gene SNPs and haplotypes associated with COPD were not detected (P > 0.05).

Association analysis of the interaction between the environment and the genes

We identified four SNPs that interact with the investigated environmental exposures in COPD.

Two SNPs in the CRP gene significantly interact with the effects of environmental exposure on COPD. Logistic regression analysis and crossover analysis revealed that the RERI (95% CI) value of the interaction between exposure to a smoking index > 200 and the C/T + C/C genotype at the rs1205 locus in the CRP gene was 0.15 (0.07–1.00), and coexposure had positive and additive effects on COPD. The prevalence of COPD was significantly increased in people with a smoking index > 200 and C/T + C/C genotypes due to synergistic effects [OR: 1.81, 95% CI (1.183–2.79)] (Tables 2 and 3 and Figure S2 in the Supplementary Appendix). In addition, the RERI (95% CI) value of the interaction between exposure to indoor use of coal/wood/straw and the G/A + A/A genotypes at the rs1130864 locus in the CRP gene was 2.45 (0.73–4.18), and coexposure presented positive and additive interactions with COPD. The prevalence of COPD was significantly increased in people exposed to indoor use of the coal/wood/straw and G/A + A/A genotypes because of synergistic effects [OR: 3.01, 95% CI (1.39–6.75)] (Tables 3 and 4 and Figure S3 in the Supplementary Appendix).

Table 2 Interaction between smoking exposure and SNPs on COPD
Table 3 Interaction between smoking and SNPs on COPD
Table 4 Interaction between indoor fuel use and SNPs on COPD

The other two SNPs were located in the FCAR gene and FCGR2B gene. The OR (95% CI) value of the interaction between exposure to a smoking index > 200 and the C/G + G/G genotypes at the rs11084376 locus in the FCAR gene was 0.54 (0.29–0.97), and coexposure had negative and multiplicative effects on COPD. It can be inferred that when individuals are exposed to a smoking index > 200 alone, their susceptibility to COPD is significantly increased [OR: 1.89, 95% CI (1.25–2.88)], whereas when the C/G + G/G genotype is carried at the same time, the antagonistic effect of copresence leads to no difference in susceptibility to COPD. Additionally, the SI (95% CI) value of the interaction between exposure to straw burning outdoors and the A/G + A/A genotype at the rs844 locus in the FCGR2B gene was 0.295 (0.11–0.77), the OR (95% CI) value of the interaction was 0.463 (0.24–0.90), and coexposure had a negative multiplicative or additive interaction on COPD. It was suggested that when exposed to straw burning outdoors alone, the susceptibility to COPD is highly increased [OR: 2.18, 95% CI (1.39–3.46)], whereas when the A/G + A/A genotype is carried simultaneously, the effect of coexposure is not significant for susceptibility to COPD (Tables 3, 5 and Figure S4 in the Supplementary Appendix).

Table 5 Interaction between Straw burning outdoor exposure and SNPs on COPD

Association analysis of the interactions between the environment and the haplotypes

We identified only one haplotype block in the FCGR2A gene that interacted with the investigated environmental exposures to COPD (Tables 6 and 7). The SI (95% CI) value of the interaction between exposure to straw burning outdoors and the CTG + TCA + TTA + TTG haplotype in the FCGR2A gene was 0.48 (0.26–0.89), and coexposure had negative and additive effects on COPD. When exposed to straw burning outdoors alone, susceptibility to COPD is significantly increased [OR: 1.93, 95% CI (1.36–2.76)]; however, when the CTG + TCA + TTA + TTG haplotype is carried simultaneously, the antagonistic effect of coexposure causes a lower risk of COPD in the population [OR: 1.55, 95% CI (1.11–2.19).

Table 6 Interaction between environmental exposure and haplotypes on COPD
Table 7 Interaction between straw burning outdoors and the FCGR2A rs4656308-rs4656309-rs2165088 haplotypes on COPD

Discussion

In our previous study, several environmental exposure factors associated with COPD were identified. Considering the important role of genetic factors, the present study focused on the effects of the genotypes, alleles and haplotypes of the CRP/FCAR/FCGR2A/FCGR2B genes on COPD, and the environment‒gene interactions were explored.

The CRP gene is located at the 1q23.2 position. Rs1130864 (G > A) is a mutation located in the 3 ′UTR (untranslated region) of the CRP gene. CRP is an important cytokine involved in the inflammatory process of COPD. This study revealed that the rs1130864 (G > A) allele and dominant model in the CRP gene are associated with COPD; that is, rs1130864-A allele, G/A and A/A genotype carriers have increased COPD susceptibility. The concentration of CRP in the airway and blood is often elevated [3]. An increasing number of studies have confirmed that the rs1130864-A allele in the CRP gene is associated with an increase in CRP; these studies were conducted in Shanghai and Europe [17, 18]. Navarro P et al. reported that the rs1130864 genotype in the CRP gene was associated with high CRP levels through a cohort study in Spain. However, some of the findings are inconsistent with this study. Bolton CE et al. conducted a study on British men and reported that serum CRP was associated with decreased lung function, whereas the rs1130864 genotype was not associated with decreased lung function [19]. A clinical study conducted by Guo et al. in Shanghai in a Han population revealed that the rs1130864 locus was not correlated with COPD [20]. The reason for these inconsistent results may be that the populations selected for this study differed. This study chose people in cold rural areas in northern China. Therefore, this study highlights the influence of genetic differences between southern and northern Chinese people on COPD. Guo J et al. reported that rs1130864 was an independent predictor of the prognosis of ischemic stroke patients in the Chinese Han population [21]. A European study by Nimptsch K et al. revealed that the A allele at rs1130864 was associated with elevated CRP and colorectal cancer risk [22]. Our group previously demonstrated that serum CRP is a risk factor for COPD, and we speculated that the variable sequence of the rs1130864 locus in the untranslated region of the CRP gene may regulate the inflammatory response of the body and exacerbate airflow restriction by affecting posttranscriptional regulation and translation.

Rs1205 (T > C) is a mutation located in the untranslated 3 'UTR of the CRP gene. Rs2808630 (T > C) is a mutation located downstream of the CRP gene. In this study, no significant associations were found between the rs1205 or rs2808630 genotype or haplotype and COPD. However, a clinical study conducted by Yi Guo et al. in Shanghai in a Han population revealed that the CRP gene rs1205-T allele was a protective factor for COPD, and the TT haplotype frequency of rs1205 and rs2808630 was significantly greater than that of the control group [22]. The reason for these different results may be genetic differences between the northern and southern populations in China.

The FCAR gene is located at the 19q13.42 position. The eight loci included in this study are all mutation sites located in the introns of the FCAR gene. FCAR genes mediate a variety of immune system functions, including degranulation, endocytosis, phagocytosis, cytokine synthesis, and cytokine release. There are no studies on the associations of SNPs in the FCAR gene with COPD, and some evidence has revealed that these SNPs contribute to the pathogenesis of other inflammation-related diseases. The FCAR gene can increase the proinflammatory activity of IgA in serum [23] and thus affect many diseases, such as dermatitis herpetiformis, celiac disease, periodontal disease and IgA nephropathy [24]. Wu J et al. reported that the rs16986050 locus in the FCAR gene is involved in the production of proinflammatory cytokines and the release of IL-6 cytokines from neutrophils in vitro, which plays an important role in the regulation of the immune response in systemic lupus erythematosus [24]. In the present study, the rs4806606-G and rs8112766-G alleles were found to reduce susceptibility to COPD. Since these genes are located in introns of the FCAR gene, we hypothesized that they regulate gene expression through initiation and selective splicing, enhancing the anti-inflammatory effect mediated by the FCAR gene.

Both the FCGR2A gene and FCGR2B gene are located at the 1q23.3 position. Among the five FCGR2A gene loci included in this study, only rs1801274 was the mutation site of the FCGR2A gene exon, and the other four were located in introns. Among the two FCGR2B gene loci included in the study, rs844 is a mutation site located in the 3 'UTR (nontranslated region) of the FCGR2B gene, and rs16827592 is in an intergenic mutation site. At present, no studies have investigated the associations between SNPs in the FCGR2A and FCGR2B genes and COPD. Some studies have reported that FCGR2A and FCGR2B play roles in the pathogenesis of other immune-related diseases. FCGR2A is a susceptibility gene for several autoimmune diseases, such as systemic lupus erythematosus (SLE) [25], rheumatoid arthritis [26], multiple sclerosis [27], type 1 diabetes [28], and ulcerative colitis [28, 29]. The homozygous C/C at rs1801274 in the FCGR2A gene reduces susceptibility to ulcerative colitis in the Chinese population, whereas the TC haplotype consisting of the rs1801274 and rs511278 loci in the FCGR2A gene increases the risk of ulcerative colitis in the Chinese population [30]. The rs2099684 locus in the FCGR2A gene is a genetic risk factor for arteritis in the Chinese Han population [31]. The FCGR2A gene was found to be strongly associated with giant cell arteritis in Spain [32]. SNPs of the FCGR2B gene inhibit binding to its ligands IgG1, IgG2, and IgG3. The alleles rs4656308-C and rs4656309-T in the FCGR2A gene are located in introns, and our results revealed that both can reduce the susceptibility of the population to COPD. Therefore, we speculated that the regulation of gene expression through initiation and selective splicing involves FCGR2A receptors on the surface of monocytes, macrophages, neutrophils, natural killer cells, T lymphocytes and B lymphocytes and subsequently enhances the phagocytosis of FCGR2A-mediated immune complexes and the anti-inflammatory effect of B cells on antibody production.

This study revealed that a smoking index > 200 had a positive additive interaction with the rs1205 locus in the CRP gene and a multiplicative interaction with the rs11084376 locus in the FCAR gene in COPD patients. However, no independent role of the rs1205 locus or the rs11084376 locus in COPD pathogenesis was found, suggesting that susceptibility to COPD is significantly increased when individuals carry the C/T + C/C genotype at the rs1205 locus or the C/G + G/G genotype at the rs11084376 locus accompanied by a smoking index > 200. No previous study has explored the association of the rs1205 locus with COPD, and it was found to be involved in the pathogenesis of other inflammation-related diseases. A European population study by Nimptsch K et al. revealed that the rs1205-T allele was associated with reduced CRP and a lower risk of colorectal cancer [33]. Todendi PF et al. reported that in overweight and obese people, the CRP level was lower than that in those who carried the rs1205-T allele in the CRP gene [34]. A Chinese population study conducted by Tian Yaping revealed that the C/C genotype frequency at the rs1205 locus in the CRP gene in a healthy Han population was lower than that in T2DM patients [35]. Previous studies have suggested that the rs1205-C allele may promote the development of inflammatory diseases. Our results suggested that the C/T + C/C genotype at the rs1205 locus has no significant effect on COPD compared with the T/T genotype, with a synergistic effect between the C/T + C/C genotype and a smoking index > 200. The effects of the rs1205 locus and smoking alone can be amplified by additive interactions, increasing the risk of COPD.

Both exposure to coal/wood/straw and the rs1130864 locus in the CRP gene are independent risk factors for COPD, and the interaction is also a risk factor, suggesting that susceptibility to COPD increases through synergistic effects when coexposure to straw burning outdoors and that the rs844 locus in the FCGR2B gene has negative additive and multiplicative interactions with COPD. No significant effect of the rs844 locus in the FCGR2B gene was found on COPD, whereas straw burning outdoors was found to be an independent risk factor, and the interaction between the two factors was a protective factor against COPD. Previous studies on the rs844 locus and disease are lacking, but studies on the mechanism of the FCGR2B gene are lacking. The FCGR2B gene functions to produce inhibitory signals in cells and is highly expressed on B cells to inhibit the activation signals induced by immunoglobulin on the surface of the B-cell membrane. It is also expressed in monocytes, macrophages and dendritic cells to inhibit prophagocytosis signals [36,37,38]. Rs844 is a mutation site in the 3 'UTR untranslated region of the FCGR2B gene. Our results suggested that the rs844-G allele in the FCGR2B gene may promote the occurrence of inflammation by influencing posttranscriptional regulation and translation, inhibiting the activation signal induced by immunoglobulin on the surface of the B-cell membrane or inhibiting the phagocytosis signal from monocytes, macrophages and dendritic cells. However, the mutated rs844-A allele can change the function of inhibitory signals, influence the immune response, enhance immunity, inhibit inflammation, and have a weak effect on COPD alone. The smoke produced by outdoor straw burning may contain some chemical components, which can multiply the effect of the rs844 locus after it enters the body. When the two factors coexist, they exhibit antagonistic effects, which can reduce susceptibility to COPD.

Exposure to straw burning outdoors had a negative additive interaction with haplotypes at rs4656308-rs4656309-rs2165088 in the FCGR2A gene on COPD. Previous studies on the rs4656308, rs4656309, and rs2165088 loci and diseases have not been reported. The FCGR2A gene is widely expressed in monocytes, macrophages, dendritic cells, neutrophils and platelets and can induce cellular defense mechanisms [39]. This study revealed that straw burning outdoors was a risk factor for COPD. Both the rs4656308 and rs4656309 loci are protective factors for COPD. However, no effect of this haplotype on the occurrence of COPD was found. Rs4656308, rs4656309 and rs2165088 are located in introns of the FCGR2A gene. Compared with the TCG haplotype at rs4656308-rs4656309-rs2165088, the susceptibility of the population with the CTG + TCA + TTA + TTG haplotypes to COPD was significantly lower in those exposed to straw burning outdoors. This haplotype gene mutation may regulate changes in gene expression through initiation and selective splicing and enhance the phagocytosis of the FCGR2A-mediated immune complex and the immune and anti-inflammatory effects of B cells on antibody production, thus increasing the susceptibility of this population to COPD.

Admittedly, there are still limitations in this study. First, the inclusion of environmental exposure factors and genetic polymorphic loci in this study is limited, and the possibility of interactions between other environmental risk factors and genes cannot be ruled out. Second, in the genetic model construction and stratified analysis of this study, the mutation rates of some genes in the investigated population may have been lower, and the sample size after classification may have been affected to some extent. Although the sample size is sufficient, the results of this study need to be further verified in different areas and ethnicities.

Conclusions

Overall, the present study revealed that alleles and genotypes of the CRP, FCAR and FCGR2A genes can increase susceptibility to COPD in the northern Chinese population. In addition, exposure to a smoking index > 200, indoor coal/wood/straw use, and outdoor straw burning were found to have synergistic or antagonistic effects on the genotypes or haplotypes of the CRP/FCAR/FCGR2A/FCGR2B genes. The above findings can provide new ideas and references for personalized and COPD research.

Availability of data and materials

The data are available upon reasonable request.

Abbreviations

BMI:

Body mass index

CHD:

Coronary heart disease

CHF:

Congestive heart failure

COPD:

Chronic obstructive pulmonary diseases

CRP:

C-reactive protein

CVD:

Cerebrovascular disease

DNA:

Deoxyribonucleic acid

FCAR:

Fc fragment of IgA receptor

FCGR:

Fc fragment of IgG receptor

FCGR2A:

Fc fragment of IgG receptor IIa

FCGR2B:

Fc fragment of IgG receptor IIb

FEV1 :

Forced expiratory volume in one second

FVC:

Forced vital capacity

GMDR:

Generalized multifactor dimensionality reducation

HWE:

Hardy–weinberg equilibrium

IgA:

Immunoglobulin A

IgG:

Immunoglobulin G

LD:

Linkage disequilibrium

PCR:

Polymerase chain reaction

RD:

Respiratory disease

RERI:

Relative excess risk of interaction

SI:

The synergy index

SNP:

Single nucleotide polyorphism

References

  1. Calle Rubio M, Rodriguez Hermosa JL, Miravitlles M, Lopez-Campos JL. Knowledge of Chronic obstructive pulmonary disease, presence of chronic respiratory symptoms and use of spirometry among the Spanish population: CONOCEPOC 2019 study. Arch Bronconeumol. 2020. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.arbr.2021.10.003.

    Article  PubMed  Google Scholar 

  2. Hopkinson NS. COPD, smoking, and social justice. Lancet Respir Med. 2022;10(5):428–30.

    Article  PubMed  Google Scholar 

  3. Wang R, Zhang W, Li Y, et al. Evaluation of risk factors for chronic obstructive pulmonary disease in the middle-aged and elderly rural Population of Northeast China using logistic regression and principal component analysis. Risk Manag Healthc Policy. 2022;15:1717–26.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Gooptu BEU, Lomas DA. Mechanisms of emphysema in α1-antitrypsin deficiency: molecular and cellular insights. Eur Respir J. 2009;34:475–88.

    Article  PubMed  CAS  Google Scholar 

  5. Ortega-Martínez A, Pérez-Rubio G, Ambrocio-Ortiz E, Nava-Quiroz KJ, Hernández-Zenteno RJ, Abarca-Rojano E, Rodríguez-Llamazares S, Hernández-Pérez A, García-Gómez L, Ramírez-Venegas A, Falfán-Valencia R. The SNP rs13147758 in the HHIP gene is associated with COPD susceptibility, serum, and sputum protein levels in smokers. Front Genet. 2020;11:882.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Li JX, Fu WP, Zhang J, Zhang XH, Sun C, Dai LM, Zhong L, Yu L, Zhang YP. A functional SNP upstream of the ADRB2 gene is associated with COPD. Int J Chron Obstruct Pulmon Dis. 2018;13:917–25.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  7. Gao X, Wang X, Jiao N, Chen J, Sun D. Association of VEGFA polymorphisms with chronic obstructive pulmonary disease in Chinese Han and Mongolian populations. Exp Physiol. 2021;106(8):1839–48.

    Article  PubMed  CAS  Google Scholar 

  8. Rozmus D, Ciesielska A, Płomiński J, Grzybowski R, Fiedorowicz E, Kordulewska N, Savelkoul H, Kostyra E, Cieślińska A. Vitamin D binding protein (VDBP) and its gene polymorphisms-the risk of malignant tumors and other diseases. Int J Mol Sci. 2020;21(21):7822.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  9. Hlapčić I, Grdić Rajković M, Čeri A, Dabelić S, Popović-Grle S, Radić Antolic M, Knežević J, Försti A, Rumora L. Increased HSP70 and TLR2 gene expression and association of HSP70 rs6457452 single nucleotide polymorphism with the risk of chronic obstructive pulmonary disease in the Croatian population. Diagnostics. 2021;11(8):1412.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Abudureheman Z, Li L, Zhong X, Xu J, Gong H, Yilamujiang S, Ren J, Xie C, Zheng A, Tuerxun D, Abudukadeer A, Aini P, Xu A, Zou X. The rs74794265 SNP of the SREK1 Gene is associated with COPD in Kashi, China. Int J Chron Obstruct Pulmon Dis. 2021;16:2631–6.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Díaz-Peña R, Silva RS, Hosgood HD 3rd, Jaime S, Miravitlles M, Olloquequi J. HLA-DRB1 alleles are associated with COPD in a Latin American admixed population. Arch Bronconeumol. 2021;57(4):291–7.

    Article  PubMed  Google Scholar 

  12. Hage FG, Szalai AJ. C-reactive protein gene polymorphisms, C-reactive protein blood levels, and cardiovascular disease risk. J Am Coll Cardiol. 2007;50(12):1115–22.

    Article  PubMed  CAS  Google Scholar 

  13. Carlson CS, Aldred SF, Lee PK, et al. Polymorphisms within the C-reactive protein (CRP) promoter region are associated with plasma CRP levels. Am J Hum Genet. 2005;77(1):64–77.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Maliszewski CR, March CJ, Schoenborn MA, Gimpel S, Shen L. Expression cloning of a human Fc receptor for IgA. J Exp Med. 1990;172:1665–72. https://doiorg.publicaciones.saludcastillayleon.es/10.1084/jem.172.6.1665.

    Article  PubMed  CAS  Google Scholar 

  15. Li X, Ptacek TS, Brown EE, Edberg JC. Fcgamma receptors: structure, function and role as genetic risk factors in SLE. Genes Immun. 2009;10(5):380–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Lu J, Marjon KD, Marnell LL, et al. Recognition and functional activation of the human IgA receptor (FcalphaRI) by C-reactive protein. Proc Natl Acad Sci U S A. 2011;108(12):4974–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  17. Kong H, Qian YS, Tang XF, et al. C-reactive protein (CRP) gene polymorphisms, CRP levels and risk of incident essential hypertension: findings from an observational cohort of Han Chinese. Hypertens Res. 2012;35(10):1019–23.

    Article  PubMed  CAS  Google Scholar 

  18. Arouca AB, Meirhaeghe A, Dallongeville J, et al. Interplay between the Mediterranean diet and C-reactive protein genetic polymorphisms towards inflammation in adolescents. Clin Nutr. 2020;39(6):1919–26.

    Article  PubMed  CAS  Google Scholar 

  19. Bolton CE, Schumacher W, Cockcroft JR, et al. The CRP genotype, serum levels and lung function in men: the Caerphilly Prospective Study. Clin Sci. 2011;120(8):347–55.

    Article  Google Scholar 

  20. Guo Y, Gong Y, Pan C, et al. Association of genetic polymorphisms with chronic obstructive pulmonary disease in the Chinese Han population: a case-control study. BMC Med Genom. 2012;5:64.

    Article  CAS  Google Scholar 

  21. Guo J, Yu L, Zhang J, Chen N, et al. CRP gene polymorphism predicts post-stroke functional outcome in Han Chinese. Acta Neurol Scand. 2014;129(4):263–8.

    Article  PubMed  CAS  Google Scholar 

  22. Nimptsch K, Aleksandrova K, Boeing H, et al. Association of CRP genetic variants with blood concentrations of C-reactive protein and colorectal cancer risk. Int J Cancer. 2015;136(5):1181–92.

    Article  PubMed  CAS  Google Scholar 

  23. Wu J, Ji C, Xie F, et al. FcalphaRI (CD89) alleles determine the proinflammatory potential of serum IgA. J Immunol. 2007;178(6):3973–82.

    Article  PubMed  CAS  Google Scholar 

  24. Huang YT, Wright A, et al. Intraepithelial cell neutralization of HIV-1 replication by IgA. J Immunol. 2005;174:4828–35.

    Article  PubMed  CAS  Google Scholar 

  25. Harley JB, Alarcon-Riquelme ME, Criswell LA, et al. Genome-wide association scan in women with systemic lupus erythematosus identifies susceptibility variants in ITGAM, PXK, KIAA1542 and other loci. Nat Genet. 2008;40:204–10.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Alizadeh BZ, Valdigem G, Coenen MJ, et al. Association analysis of functional variants of the FcgRIIa and FcgRIIIa genes with type 1 diabetes, celiac disease and rheumatoid arthritis. Hum Mol Genet. 2007;16:2552–9.

    Article  PubMed  CAS  Google Scholar 

  27. Takai T. Roles of Fc receptors in autoimmunity. Nat Rev Immunol. 2002;2:580–92.

    Article  PubMed  CAS  Google Scholar 

  28. McGovern DP, Gardet A, Törkvist L, et al. Genome-wide association identifies multiple ulcerative colitis susceptibility loci. Nat Genet. 2010;42:332–7.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Asano K, Matsushita T, Umeno J, et al. A genome-wide association study identifies three new susceptibility loci for ulcerative colitis in the Japanese population. Nat Genet. 2009;41:1325–9.

    Article  PubMed  CAS  Google Scholar 

  30. Xia SL, Lin DP, Lin QR, et al. A case-control study on association of ulcerative colitis with FCGR2A gene polymorphisms in Chinese patients. Genet Test Mol Biomarkers. 2018;22(10):607–14.

    Article  PubMed  CAS  Google Scholar 

  31. Chen S, Wen X, Li J, et al. Association of FCGR2A/FCGR3A variant rs2099684 with Takayasu arteritis in the Han Chinese population. Oncotarget. 2017;8(10):17239–45.

    Article  PubMed  Google Scholar 

  32. Morgan AW, Robinson JI, Barrett JH, et al. Association of FCGR2A and FCGR2A-FCGR3A haplotypes with susceptibility to giant cell arteritis. Arthritis Res Ther. 2006;8:R109.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Nimptsch KAK, Boeing H, Janke J, et al. Association of CRP genetic variants with blood concentrations of C-reactive protein and colorectal cancer risk. Int J Cancer. 2015;136(5):1181–92.

    Article  PubMed  CAS  Google Scholar 

  34. Todendi PFPL, Klinger EI, Reuter CP, et al. Low-grade inflammation markers in children and adolescents: influence of anthropometric characteristics and CRP and IL6 polymorphisms. Cytokine. 2016;88:177–83.

    Article  PubMed  CAS  Google Scholar 

  35. Zhang Y, Tian Y, Dong J, Li LQ, Deng XX. Studies of CRP levels and CRP gene polymorphism rs1205 in type 2 diabetes mellitus patients. Label Immunoass Clin Med. 2009;16(01):19–22.

    Google Scholar 

  36. Tsang-A-Sjoe MWNS, Bultink IE, Geissler J, et al. Fc-gamma receptor polymorphisms differentially influence susceptibility to systemic lupus erythematosus and lupus nephritis. Rheumatology. 2016;55(5):939–48.

    Article  PubMed  CAS  Google Scholar 

  37. Su KYH, Li X, Li X, et al. Expression profile of FcgammaRIIb on leukocytes and its dysregulation in systemic lupus erythematosus. J Immunol. 2007;178(5):3272–80.

    Article  PubMed  CAS  Google Scholar 

  38. Van der Heijden JBW, Geissler J, de Boer M, et al. Phenotypic variation in IgG receptors by nonclassical FCGR2C alleles. J Immunol. 2012;188(3):1318–24.

    Article  PubMed  Google Scholar 

  39. Munde EOOW, Raballah E, Anyona SB, et al. Association between Fcγ receptor IIA, IIIA and IIIB genetic polymorphisms and susceptibility to severe malaria anemia in children in western Kenya. BMC Infect Dis. 2017;17(1):289.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank all the participants, staff, and investigators for their contributions to the study.

Funding

This work was supported by the National Health Commission of the People’s Republic of China and funded by the Central Transfer Payments Endemic Disease Project (N/A). This work was also funded by special development funds of local colleges from the central government (Study on prevention and control of major diseases in Heilongjiang Province (N/A). This work was also funded by the Heilongjiang Academy of Medical Sciences (Chronic Disease Cohort study in rural areas of Northeast China, CR201801).

Author information

Authors and Affiliations

Authors

Contributions

Conception, design, data curation, methodology, supervision, reviewing & editing the manuscript: WZ, DS. Conception, methodology, data acquisition, analysis, draft and interpretation, revising of the manuscript: RW, YL, YJ. Data acquisition & interpretation, reviewing the manuscript: XL, HF, ZJ, BL, CL, YS, FC, CZ. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Dianjun Sun or Wei Zhang.

Ethics declarations

Ethics approval and consent to participate

The study was conducted according to the guidelines of the Declaration of Helsinki and was approved in June 2018 by Harbin Medical University (research project: hrbmuecdc20180601). All participants provided informed consent prior to starting the study.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, R., Li, Y., Jiang, Y. et al. Interaction study of the effects of environmental exposure and gene polymorphisms of inflammatory and immune-active factors on chronic obstructive pulmonary disease. Respir Res 26, 5 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12931-024-03079-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12931-024-03079-4

Keywords