Chronic obstructive pulmonary disease (COPD) is a complex disease characterized by limited airflow and is influenced by genetic and environmental factors. The purpose of this study was to investigate the effects of gene polymorphisms in MIR5708 and MIR1208 on COPD risk.
MethodsFour single nucleotide polymorphisms (SNPs) in MIR5708 (rs6473227 and rs16907751) and MIR1208 (rs2608029 and rs13280095) were selected and genotyped among 315 COPD patients and 314 healthy controls using the Agena MassARRAY platform. SPSS 18.0 was used for statistical analysis and data processing. Odds ratios (ORs) and 95% confidence intervals (CIs) were calculated to assess the association between genetic variants of MIR1208 and MIR5708 and COPD risk.
ResultsThe results suggested that rs16907751 variants in MIR5708 contributed to an increased susceptibility to COPD in the allelic (P = 0.001), co-dominant (homozygous) (P = 0.001), dominant (P = 0.017), recessive (P = 0.002), and additive (P = 0.002) models. The effects of MIR5708 and MIR1208 gene polymorphisms on the risk of COPD were age-, sex-, smoking status-, and BMI-related. Furthermore, the C-A and G-A haplotypes of rs2608029 and rs13280095 in MIR1208 were identified as risk factors for COPD in the population over 70 years (P = 0.029) and in women (P = 0.049), respectively. Finally, significant associations between rs16907751genotypes with pulse rate and forced expiratory volume in 1 s were found among COPD patients.
ConclusionGenetic polymorphisms in MIR5708 and MIR1208 are associated with increased risk of COPD in China.
Chronic obstructive pulmonary disease (COPD) is a multifactorial chronic respiratory disease characterized by irreversible pulmonary airflow obstruction.1,2 Many studies have shown that the development of COPD is influenced by environmental factors and a complex set of genetic traits.2–4 Recently, genome-wide association and exome sequencing studies have been used to identify the hereditary factors of COPD among different populations, but the effect of genetic variants on COPD is still unclear.5
Micro RNAs (miRNAs) are a class of endogenous non-coding RNAs, 22–25 nucleotides in length that can bind to the 3′-UTR region of mRNA transcripts, thereby inhibiting their translation and expression.6,7 Variants in miRNAs are thought to be involved in post-transcriptional regulation and influence disease susceptibility.7,8 Zhou et al and Wang et al found that the miR-146a rs2910164 single nucleotide polymorphism (SNP) is associated with improved lung function in smokers with COPD.7,8 Related studies in the Korean population have suggested that gene polymorphisms in miR-196a2, miR-146a, and miR-499 could be associated with asthma phenotypes.6 Fawzy et al also found a significant association between the miR-196a2 rs11614913 polymorphism and the bronchodilator response in Egyptians with COPD.9 Therefore, it is important to explore the effect of miRNA variants on susceptibility to COPD.
MIR5708 and MIR1208 are miRNAs located on chromosome 8. MIR5708 has been identified as a risk locus for rheumatoid arthritis through gene-based association testing.10MIR1208 is thought to be involved in the regulation of tumor-related pathways by binding to the upstream gene circMTUS1.11 Relevant studies have also recognized MIR1208 as a tumor suppressor gene that can increase sensitivity of renal cancer cells to cisplatin.12 Moreover, genetic polymorphisms of MIR1208 rs2648841 have been found to be related to chemotherapy toxicity in children with acute lymphoblastic leukemia.13 However, to the best of our knowledge, the influence of MIR5708 and MIR1208 polymorphisms on the risk of COPD has not been studied.
Therefore, we performed this case-control study to investigate the relationship between MIR5708 (rs6473227 and rs16907751) and MIR1208 (rs2608029 and rs13280095) polymorphisms and susceptibility of Chinese to COPD, and to explain the roles of genetic factors in the pathogenesis of COPD through statistical analysis.
MethodsStudy populationThis study recruited 314 COPD patients (76% men) and 315 healthy controls (75% men). Patients with appropriate symptoms including dyspnea, cough, sputum production, wheezing, and chest tightness were considered as suspected cases of COPD, which were further confirmed by a post-bronchodilator fixed ratio of forced expiratory volume in 1 second/forced vital capacity (FEV1/FVC) < 0.70.14 Meanwhile, patients with other respiratory diseases, autoimmune diseases, cardiovascular diseases, and others were excluded. Healthy controls were recruited from the physical examination center of Hainan General Hospital, matched by age and sex with FEV1/FVC > 0.70, and no disease or family history. Demographic and clinical characteristics of all subjects were surveyed, and informed consent was obtained regarding the purpose and procedure of this study. This study was approved by the Hainan General Hospital and followed the Declaration of Helsinki guidelines.
SNP selection and genotypingThe selection of SNPs was based on haplotype or genotype data. Two SNPs in MIR1208 (rs2608029 and rs13280095) and two in MIR5708 (rs6473227 and rs16907751) were selected from the database of the 1000 Genomes Project (https://www.internationalgenome.org/) with the criteria of minor allele frequency (MAF) > 0.05, and r2 > 0.80. The online design of the corresponding primers was performed in AgenaCx Tools (https://agenacx.com/online-tools/), and the primer sequences are listed in Supplementary Table 1. Genomic DNA of all participants was extracted and purified from cryopreserved serum using a GoldMag-Beads Kit (GoldMag, Xi'an, Shaanxi, China), and spectrophotometry (Beckman, Fullerton, CA, USA) was used to determine DNA concentration. SNP genotyping and data collection were conducted using the Agena MessARRAY platform and Agena Bioscience TYPER software (Agena Bioscience, San Diego, CA, USA), respectively.
Statistical analysisThe statistical data of this case-control study were processed using SPSS software (version 22.0, SPSS Inc., Chicago, IL, USA), Microsoft Excel software, and PLINK software (version 1.07) (http://www.cog-genomics.org/plink2/). Independent-sample t-test and chi-squared test were used to assess differences in clinical indexes between COPD patients and healthy controls, where appropriate. The association between each SNP and the risk of COPD was evaluated by logistic regression analysis under different genetic models, and was represented by odds ratios (ORs) and 95 % confidence intervals (CIs) after adjusting for age and sex. Hardy-Weinberg equilibrium (HWE) for each SNP in the control group was determined using Fisher's exact test. Statistical significance was set at p < 0.05.
ResultsClinical indexes of participants and genotypic characteristics of selected SNPsDemographics and clinical indexes of all 629 participants comprising 315 COPD patients (mean age 71.93 ± 10.11 years) and 314 healthy controls (mean age 71.23 ± 6.83 years old) are displayed in Table 1. Statistically, cases and controls were matched by age (P = 0.306) and sex (P = 0.908).
Clinical characteristics of cases and controls.
Variables | Case (n = 315) | Control (n = 314) | P |
---|---|---|---|
Age (Mean ± SD) | 71.93 ± 10.11 | 71.23 ± 6.83 | 0.306a |
≤ 70 | 127 (40%) | 137 (44%) | |
> 70 | 188 (60%) | 177 (56%) | |
Gender | 0.908b | ||
Male | 239 (76%) | 237 (75%) | |
Female | 76 (24%) | 77 (25%) | |
BMI (kg/m2) | |||
≤ 24 | 251 (80%) | 67 (21%) | |
> 24 | 29 (9%) | 78 (25%) | |
Tobacco smoking status | |||
Yes | 147 (47%) | 52 (17%) | |
No | 166 (53%) | 118 (38%) | |
Comorbidity | |||
Yes | 93 (30%) | ||
No | 174 (55%) |
SD, standard deviation; BMI, body mass index.
Table 2 lists the detailed information of the four SNPs in MIR2708 and MIR1208; MAFs > 0.05 and SNPs in the control group conformed to HWE (P > 0.05). Moreover, a significant association between the rs16907751 polymorphism and COPD risk was found in the allelic model (P = 0.001).
Information regarding SNPs in MIR5708 and MIR1208.
SNP, single nucleotide polymorphism; MAF, minor allele frequency; HWE, Hardy-Weinberg equilibrium; OR, odds ratio; 95% CI, 95% confidence interval.
Bold P < 0.05 indicates statistical significance.
We carried out an association analysis on genetic variants and COPD susceptibility and the relevant data are presented in Table 3. The results suggested that the TT genotype at rs16907751 in MIR5708 was a risk factor for COPD in the co-dominant (OR = 4.62, 95% CI = 1.84–11.62, P = 0.001), dominant (OR = 1.49, 95% CI = 1.07–2.08, P = 0.017), recessive (OR = 4.27, 95% CI = 1.71–10.66, P = 0.002), and additive (OR = 1.57, 95% CI = 1.19–2.07, P = 0.002) models.
Association of MIR5708 and MIR1208 polymorphisms with COPD risk.
SNP, single nucleotide polymorphism; OR, odds ratio; 95% CI, 95% confidence interval; HOM, homozygous; HET, heterozygous
P value was calculated by logistic regression analysis with adjustment for age and sex.
Bold P < 0.05 indicates statistical significance.
To further explore the association between SNPs and COPD risk, we conducted a stratified analysis according to age, sex, BMI, and smoking status. According to the data listed in Table 4, rs6473227 in MIR5708 was associated with an increased risk of COPD in the population with a body mass index (BMI, kg/m2) less than or equal to 24 in the allelic (OR = 1.66, 95% CI = 1.13–2.44, P = 0.010), homozygous co-dominant (OR = 2.97, 95% CI = 1.34–6.62, P = 0.008), recessive (OR = 2.15, 95% CI = 1.06–4.35, P = 0.033), and additive (OR = 1.72, 95% CI = 1.16–2.55, P = 0.007) models. Moreover, rs16907751 in MIR5708 was found to be significantly associated with an increased risk of COPD in males, non-smokers, patients older than 70 years, and individuals with a BMI ≤ 24. In addition, females heterozygous at rs2608029 (OR = 2.24, 95% CI = 1.09–4.64, P = 0.029) were more susceptible to COPD risk in both dominant (OR = 2.19, 95% CI = 1.07–4.47, P = 0.031) and additive (OR = 1.99, 95% CI = 1.02–3.91, P = 0.045) models.
Stratification analysis of the associations between MIR1208 and MIR5708 polymorphisms and COPD risk.
SNP, single nucleotide polymorphism; OR, odds ratio; 95% CI, 95% confidence interval; HOM, homozygous; HET, heterozygous
“-” indicates no results.
P value was calculated by logistic regression analysis with adjustment for age and sex.
Bold P < 0.05 indicates statistical significance.
We continued with a haplotype analysis and found a block (rs2608029 and rs13280095) in MIR1208 (Table 5). The C-A haplotype of rs2608029 and rs13280095 was found to be related to an increased risk of COPD (OR = 2.59, 95% CI = 1.10–6.06, P = 0.029) in the population over 70 years, whereas Grs2608029Ars13280095 was found to be a risk factor for COPD in women.
Haplotype analysis of the association between MIR5708 and MIR1208 polymorphisms and COPD risk.
SNP, single nucleotide polymorphism; OR, odds ratio; 95% CI, 95% confidence interval; BMI, body mass index
P value was calculated by logistic regression analysis with adjustment for age and sex.
Bold P < 0.05 indicates statistical significance
Then, the effects of genetic variants in MIR5708 and MIR1208 on COPD were investigated in terms of clinical indices, including respiratory rate (RR), pulse rate (PR), FVC, FEV1, and FEV1/FVC (Table 6). The results suggested that rs16907751 in MIR5708 was significantly associated with PR (P = 0.009) and FEV1 (P = 0.049) in the genotypic model. Furthermore, no association was found between rs6473227, rs2608029, and rs13280095 and the clinical indices of COPD.
Association analysis on the clinical indexes of COPD and gene polymorphisms.
RR, respiratory rate; PR, pulse rate; FVC, forced vital capacity; FEV1, forced expiratory volume in 1s.
Bold P < 0.05 indicates statistical significance.
According to the latest Global Health Observatory data, COPD was ranked as the third leading cause of death worldwide, accounting for approximately 5% of all deaths globally in 2015 (3.17 million deaths). The prevalence and burden of COPD are projected to increase over the coming decades.15 It has been clearly demonstrated that the interaction between genotype and environment plays an important role in COPD phenotypes.16,17 Genetic background plays a critical role in the development of COPD, and related studies have proposed significant differences in genetic susceptibility to COPD among different races and ethnicities.18 Therefore, the effect of genetic polymorphisms of candidate genes on the risk of COPD among Chinese was examined.
As a result, statistical analysis of the genotyping results of 315 COPD patients and 314 healthy controls among Chinese patients identified MIR5708 rs16907751 as a risk factor of COPD, and significant associations between rs16907751 and clinical indexes including PR and FEV1 in COPD patients were found. The examined mutation loci have been rarely reported; only rs6473227 showed a close relationship with atopic dermatitis in a related case-control analysis.19 Our results suggest that rs6473227 may contribute to an increased risk of COPD in individuals with a BMI less than 24. A related study pointed out that neuroimmune interactions are implicated in chronic inflammations such as atopic dermatitis, COPD, and asthma20; therefore, we hypothesized that MIR5708 rs6473227 may increase susceptibility to inflammation-related disorders. In addition, we found that MIR5708 rs16907751 may increase the risk of COPD in patients with normal BMI. There have been studies supporting our findings and reported that low BMI is a risk factor for accelerated decline in lung function compared to normal BMI,21 thereby suggesting a protective role of high BMI in COPD patients.
Tobacco smoking is widely recognized as the most important risk factor for COPD,22 but it is not the only risk factor for COPD. Qian et al. concluded that both non-smokers and smokers are likely to suffer from COPD, which may be the result of differential miRNA expression between.23 At present, there has been no report on the association between miR5708 and miR1208 and the risk of disease COPD. However, our study further provided evidence that allele T of MIR5708 rs16907751 may increase susceptibility of nonsmokers to COPD. There have been reports confirming that up to 30% of patients with COPD do not have a smoking history.18 Therefore, multiple genes and exposure to environmental and occupational factors are thought to jointly affect COPD development.
Furthermore, we also found that the effects of MIR5708 and MIR1208 on the risk of COPD were age- and sex-dependent. Related epidemiological studies have shown that the prevalence of COPD increases with age,24 and Mercado et al. also showed that aging of the lung resulted in loss of lung elasticity and reduced ability to respond to environmental stress and damage.25 Our findings also indicated MIR5708 rs16907751 as a risk factor for COPD in patients over 70 years of age. Based on the stratified analysis, our results suggested that MIR5708 rs16907751 and MIR1208 rs2608029 are associated with increased COPD risk in men and women, respectively, and further confirmed gender differences in COPD susceptibility conferred by these two SNPs. Aryal et al. stated that the differences in sex were derived mainly by smoking status, hormone levels, and behavioral differences.26 Current data indicate that in China, the prevalence of COPD in men is higher than that in women, but the difference varies with exposure to risk factors and socioeconomic development in different geographic regions.27,28 However, in our current study, we did not focus on the effect of environment, occupational differences, and regional constraints. Moreover, the small sample size and incomplete sample information in some cases are also limitations of our study. Subsequent functional experiments are needed to verify the effects of these SNPs on COPD.
ConclusionIn summary, this study is the first to indicate that genetic variants in MIR5708 and MIR1208 are associated with an increased risk of COPD in the Chinese population, and the effects of these SNPs are related to age, sex, smoking status, and BMI.