1. Introduction
Lung function is usually assessed using spirometry, which reflects the physiological state of the lungs and airways [1]. Chronic obstructive pulmonary disease (COPD), characterised by airflow limitation as assessed by spirometric measurement, imposes a large burden in terms of disability and is a leading cause of death worldwide [2,3,4]. Several studies have investigated the genetics of lung function by estimating heritability and performing a genome-wide association study (GWAS). Family-based studies have provided evidence of significant genetic effects on lung function, with estimates of heritability () ranging from 0.28 to 0.52 for forced expiratory volume in 1 s (FEV1), 0.4 to 0.54 for forced vital capacity (FVC), and 0.24 to 0.45 for FEV1/FVC [5,6,7,8]. In addition, Zhou et al. estimated the single nucleotide polymorphism (SNP) heritability of FEV1 and FEV1/FVC in a non-Hispanic white cohort of smokers with and without COPD and found that both were approximately 0.37 [9].
Genome-wide association studies have identified specific SNPs associated with FEV1, FVC, and FEV1/FVC [10]. Several studies have detected genetic signals associated with lung functions, including variants in KANSL1, HLA-DQ, NPNT, TET2, TSEN54, and HHIP, in non-smokers, and RBM19-TBX5 in heavy smokers [11,12,13,14,15,16]. However, these studies focused on cross-sectional lung function, whereas just a few studies have focused on the genetic effects underlying the rate of decline in lung function [17].
Epidemiologic and genetic investigations of pulmonary function require longitudinal observation of lung function trajectories. Although a few studies followed subjects for more than 10 years [17,18], to the best of our knowledge, none of these studies assessed heritability in unrelated subjects, or subjects with non-European ancestry. The advantage of using unrelated subjects to estimate heritability is that they do not suffer from confounding problems caused by epistatic interactions or shared environment, which are present in family and twin studies [19,20,21,22].
In the present study, we utilised data from the Korean Genome and Epidemiology Study (KoGES) cohort [23]. All independent subjects in the KoGES cohort were enrolled when they were approximately 40 years of age and followed for 14 years. Eight biennially measured lung functions were analysed in this study. Using KoGES cohort data, we evaluated the importance of genetic components on the subject-specific means and annual change rates of lung function. Specifically, we defined two different heritability estimation approaches for the subject-specific mean, and annual change rates of longitudinally observed lung function changes were based on the method of Yang et al. [24]. We also performed GWAS on these lung functions.
2. Results
2.1. Characteristics of Study Subjects
Study subjects were obtained from the Ansan (urban community) and Ansung (rural community) cohorts, which were followed up biennially. Each subject underwent a maximum of eight measurements. The numbers of observed measurements of the eight lung functions at each period are shown in Figure 1a and Supplementary Table S1. Notably, the sample size for the post-bronchodilator tests post-FVC, post-FEV1, and post-FEV1/FVC was small; these test results were available only from subjects in Ansung. The average profile plots of the eight lung functions at each period showed a decreasing trend for all lung functions (Figure 1b,c).
The statistical analyses required estimates of the mean and annual change rates for each subject, for which we included 6622 participants (3569 Ansan and 3053 Ansung) who were followed up ≥3 times (Table 1). Since some of the participants had a missing smoking status, 5104 participants were analysed in the smoking-related studies (3009 never-smokers and 2095 ever-smokers). In addition, 3,352,722 SNPs were considered in the analysis after quality control and imputation.
2.2. Heritability Estimates of Subject-Specific Means and Annual Change Rates
To estimate the importance of genetic determinants for the eight lung function traits, we calculated the subject-specific mean () and annual change rate () of each trait. Supplementary Table S2 shows the statistics of subject-specific mean and annual change rate. The estimated mean value of FVC across the eight periods was 3.47 and the mean value of FEV1, across the eight periods, was 2.7. Both of the estimated mean values of these two traits were in the range of average normal values, even though both were slightly lower. The mean of FEV1/FVC ratio was 77.97%, which is also within the normal range above the cutoff of 70% by the GOLD criteria [25].
Both and were transformed using the rank-based inverse normal transformation and used to estimate and . Figure 2a shows the estimates of heritability represented by . All lung function traits were significant at the false-discovery rate (FDR)-adjusted at 0.05 significance level. Post-bronchodilator FVC exhibited the largest (0.325, p = 1.16 × 10−5) followed by post-bronchodilator FEV1/FVC (0.314, p = 1.86 × 10−5). Supplementary Figure S2 compares the cross-sectional SNP heritability at each period with . The results showed that the mean cross-sectional SNP heritability at each period and did not differ markedly. For each lung function trait, was < (Figure 2b). Post-bronchodilator FEV1/FVC exhibited the highest (0.176, p = 0.0099) followed by FEV1/FVC (0.158 p = 4.91 × 10−5; Detailed estimation results are listed in Supplementary Tables S3 and S14). Supplementary Table S4 shows the results when the subject-specific mean was also included as a covariate; here, was used as a response variable. The results showed that remained significant for FEV1/FVC and post-bronchodilator FEV1/FVC even after adjusting for the baseline effect.
For lung function traits with significant (FEV1/FVC, and post-bronchodilator FEV1/FVC), we calculated the correlation () between genetic components for the subject-specific mean () and annual change rate (). Supplementary Figure S3 shows the phenotypic correlation between subject-specific means and annual change rates of FEV1/FVC and post-bronchodilator FEV1/FVC; their correlations without any adjustments were 0.24 and 0.22, respectively. Table 2 shows and . indicates a genetic correlation where the relative proportions are shared between subject-specific means and annual change rates. Since we wanted to check genetic components between and and their shared genetic effects, we did not adjust the baseline in the genetic correlation model. The results showed that ≥50% of genetic components were significantly shared between subject-specific means and annual change rates ( = 0.628, p = 4.59 × 10−5 for FEV1/FVC; = 0.466, p = 0.022 for post-bronchodilator FEV1/FVC). Table 2 also shows the residual phenotypic correlations () between subject-specific means and annual change rates. The indicates the relative proportions of environmental variances shared between subject-specific means and annual change rates. Residual phenotypic correlations were much smaller than ( = 0.117 for FEV1/FVC; = 0.155 for post-bronchodilator FEV1/FVC), which indicates that subject-specific means and annual change rates may be affected by different environmental factors. We also checked other traits but did not find any significant genetic correlations.
2.3. Effect of Smoking on Heritability Estimates of the Eight Lung Function Traits
Subjects were separated into ever- and never-smokers and subgroup analyses were performed. Supplementary Table S5 shows the subject-specific means () and annual change rates () of subgroups. Significant differences between ever- and never-smokers were observed for all subject-specific means and annual change rates. For both groups, SNP heritability of subject-specific means () and SNP heritability of annual change rates () were separately estimated (Figure 3 and Supplementary Tables S6 and S7). The estimated in the never-smoker group was higher than that of ever-smokers, except for FEV1/FVC and post-bronchodilator FEV1/FVC. However, except for the post-bronchodilator FEV1/FVC of the never-smoker group, of the other traits was not significant at the 0.05 level.
We also evaluated the heritability of the SNP–smoking interaction (), which estimates the variances explained by the SNP and smoking interaction effect, for lung function traits with significant and . All eight lung function traits exhibited significant values but none showed a significant (Supplementary Table S8). For the annual change rate, only post-bronchodilator FEV1/FVC showed significant estimation in and . Table 3 shows that of post-bronchodilator FEV1/FVC was significant at the 0.05 significance level ( = 0.402, p = 0.02). The other traits are listed in Supplementary Table S8. The significant results in the interaction analyses indicate the amount of genetic variance that would be affected by smoking.
2.4. Genome-Wide Association Studies of Subject-Specific Means and Annual Change Rates
To identify the disease susceptibility loci (DSL) of FEV1, FEV1/FVC, post-bronchodilator FEV1, and post-bronchodilator FEV1/FVC, we conducted GWAS with the subject-specific mean () and annual change rate () for each trait. Supplementary Tables S9 and S12 show the genome-wide significant SNPs at the significance level (for each significant trait, top 40 variants after clumping are listed in Supplementary Tables S10, S11 and S13). For subject-specific means, four DSL were identified. The genome-wide significant results of CASC17 [26,27,28], FAM13A [12,26,29], PID1 [30], and TNXB [31] have been previously reported. Genome-wide significance of rs62201158 for FEV1/FVC was reported for the first time. For annual change rates, rs2445936 was significantly associated with post-bronchodilator FEV1/FVC. Supplementary Figures S4 and S5 show the Manhattan and quantile–quantile plots of significant traits, which indicates that our GWAS preserves the nominal significance level. In addition, the inflation factors (lambda) were close to 1.
3. Discussion
In this study, we used two different heritability analyses of subject-specific means and annual change rates to estimate heritability for eight lung function traits. Using a 14-year follow-up of Korean population-based cohort dataset, we estimated the genetic correlations between the mean and annual change rate; subgroup and interaction analyses were also performed according to smoking status. To the best of our knowledge, this is the first study to report the heritability of the decline in lung function and other measurements in a population-based setting and Asian ancestry.
In our analysis, we found that the heritability of subject-specific means was significant for all eight lung function traits. The heritability of change in lung function, FEV1/FVC, and post-bronchodilator FEV1/FVC showed significant results. By comparing the estimated heritability via subgroup analyses between never- and ever-smokers, we found that the heritability of subject-specific means was higher in the never-smoker group. In the interaction analysis, we observed that the annual change in post-bronchodilator FEV1/FVC, FEV1/FVC, and MVV showed a significant or near-significant SNP–smoking interaction, thus allowing us to infer the amount of genetic variance that would be affected by smoking.
Previous cross-sectional estimates of pedigree-based heritability, compared with SNP-based heritability, have shown that lung function traits range from approximately 20 to 40% [9,17,32]. In the present study, the estimation of SNP heritability of subject-specific means of the lung function traits showed values similar to those in the previous studies. The overall estimated heritability for the eight traits ranged from approximately 9% (MVV) to 33% (post-bronchodilator FVC). The SNP heritability of annual change rates was lower than that of subject-specific means and ranged from approximately 1% (MVV) to 18% (post-bronchodilator FEV1/FVC), similar to previous studies on heritability decline [17]. This may indicate fewer genetic effects on the annual changes in the lung function than on the mean values. Interestingly, post-bronchodilator FEV1/FVC exhibited the most substantial estimated heritability for annual change rate among the eight traits; FEV1/FVC also displayed significant heritability. Moreover, the heritability of these traits was more than half of the subject-specific means, which suggested that genetic effects played an essential role in the annual change of these lung function traits.
Generally, the lung function of healthy subjects peaked in early adulthood, followed by a steady decline. Subjects who failed to reach the predicted level of peak lung function at an early age had a higher chance of being affected by COPD [33]. Recently, several studies showed that the decline in the FEV1 and FEV1/FVC patterns differ in subjects, with four different possible trajectories obtained using 29-year follow-up data [34]. These observations showed that the traditional notion of COPD being primarily caused by smoking was incomplete; besides, the peak level of lung function attained in youth constituted a major determinant of disease susceptibility [35,36]. It was also shown that genetic effects substantively contribute to these trajectories (up to 83%) [34]. Thus, longitudinal and perspective trajectory analysis is important to understand the mechanism of lung function.
Environmental factors could influence the estimation of lung function heritability. Traditionally, smoking is considered the most crucial factor to influence the decline in lung function [37]. However, up to 30% of patients with COPD worldwide never smoked [38]; therefore, exposure to other inhaled particles and gases can also lead to a decline in lung function [39]. In our study, we estimated SNP heritability for both subject-specific mean and annual change rates in never- and ever-smoker groups. We found that the two types of heritability analyses yielded different results, and that for some traits, different results were obtained between the two smoking groups. The heritability estimates of subject-specific means for FVC was higher (4%) for never- than ever-smokers. Moreover, the estimated heritability of FEV1/FVC and post-bronchodilator FEV1/FVC of ever-smokers was slightly, albeit significantly (at the 0.05 significance level) higher than that of never-smokers, with differences of 7% and 9%, respectively. Conversely, the heritability estimates of annual change rates did not reach significance, potentially owing to insufficient sample size. We also performed SNP–smoking interaction analysis, and annual change rates of post-bronchodilator FEV1/FVC, FEV1/FVC, and MVV showed significant and near-significant results, respectively, suggesting that the genetic variance for these three traits was affected by smoking in the middle-aged population.
Subject-specific means were strongly and positively genetically correlated (≥0.466), with annual change rates for FEV1/FVC and post-bronchodilator FEV1/FVC. As these metrics consist of genetic and environmental components, the positive correlations between their SNP effects indicate that subjects with higher genetic risk for subject-specific means of traits, such as FEV1/FVC and post-bronchodilator FEV1/FVC, tend to have a higher genetic risk for annual change rates as well. If the genetic effect for subject-specific means is fixed, the conditional variance of genetic effect on annual change rates () becomes 0.0932 and 0.1364 for FEV1/FVC and post-bronchodilator FEV1/FVC, respectively, which indicates that the heritability of annual change rates is partially independent of the subject-specific mean genetic components. Therefore, we can conclude that a large proportion of the genetic effect on annual change rates is due to effects on subject-specific means. Our participants were aged 40 or older, with decreasing lung functions. Individuals with a higher lung function peak experience a lesser decline in lung function and our results can explain why the individuals with low peak lung function at an early age are often at a higher risk of developing disordered lung function, which is one of the parameters of lung function trajectories [34,40]. However, in view of different observations been reported, further studies are necessary to clarify this further.
We also performed genome-wide association analyses using subject-specific means and annual change rates of the lung function as the response, where sex, average age, and height were adjusted as covariates. Because of the small sample size, GWAS with the smoking subgroup were not conducted, and it is discussed in our future study as well as a validation study. The GWAS results for all samples, with a genome-wide significance level of 5 × 10−8, are shown in Supplementary Files (Supplementary Tables S9 and S12; Supplementary Figures S4 and S5). In the subject-specific means analysis, we detected that the association of rs4793538 with FEV1 was the most significant among the eight traits. This variant is located near CASC17, and interestingly, located in the upstream of SOX9, which was reported to be associated with FEV1 in previous studies [26,41]. SOX9 was shown to be upregulated in adenocarcinoma of the lung, and its gene expression was associated with cell proliferation and lung development [42]. For the FEV1/FVC ratio, we discovered rs2704589 within FAM13A, which was found to be associated with FEV1/FVC ratio and COPD susceptibility in the prior GWAS [12,43,44,45]. A previous study with a different Korean population found that expression of FAM13A is much higher in the lung tissue of COPD cases compared to controls, and the increased gene expression levels are associated with the risk allele [45]. We also detected that rs62201158, located between SPHKAP and PID1, are associated with FEV1/FVC for the first time. In the prior investigations, PID1 was associated with FEV1/FVC ratio [12] and was also found to be overexpressed in COPD lung tissues [45]. Previous studies have indicated that PID1 has a role in tissue homeostasis and cell growth [46], and that overexpression of PID1 induces mitochondrial dysfunction in adipocytes [47]. Furthermore, PID1 overexpression may contribute to COPD development via mitochondrial malfunction and excess reactive oxygen species [45,48]. Another identified significant SNP, rs9391733, which mapped to TNXB, was also previously reported to be associated with lung function traits (FEV1 and FEV1/FVC) and COPD [41,49,50].
In the annual change rates analysis, we identified that SNP rs2445936 near CEP164 was associated with post-bronchodilator FEV1/FVC. CEP164 was previously reported to be critical for ciliogenesis during differentiation of airway multiciliated cells implicated in respiratory function [51]. In COPD, the ciliary structure and function are impaired, affecting the clearance of harmful inhaled debris and particles, such as tobacco smoke and fumes, leading to chronic inflammation and hyperinflation of the lung [52,53]. Therefore, the link between post-bronchodilator FEV1/FVC and the locus in CEP164 in our analysis provides the evidence that the cilia in the airway play an important role in lung function decline.
One of the limitations of this study was sample size. Although the sample size was sufficiently large to support analysis of the complete data, it was limited for subgroup analysis. This led to significant standard errors for heritability estimation in the analysis. Previous studies have also reported that SNP heritability estimates are affected by the sample size when using GCTA [21]. Thus, to yield more reliable results, a larger sample size will be needed in future studies. In addition, the study population was based on a heterogenous community cohort, and smoking exposure was confounded with sex and height variables.
In summary, we estimated SNP heritability for eight lung function traits using a two-stage method to estimate cross-sectional averages and annual change rates in a Korean population-based longitudinal dataset. Our findings will help to understand aetiology of lung function decline and may facilitate the discovery of genetic factors influencing lung function-related traits. A larger sample size and novel statistical approaches will be required in future studies to validate and extend our findings.
4. Methods
4.1. Study Subjects
The KoGES cohort [23] consists of participants residing in Ansan (an urban area) and Ansung (a rural area) in the Gyeonggi Province of South Korea. KoGES was designed to investigate genetic, environmental, and behavioural risk factors of common complex diseases and causes of death in Koreans with long-term follow-up [54]. The baseline survey was completed in 2001–2002; 10,030 participants aged 40–69 years were recruited and followed for 14 years. Measurements were obtained from each participant every 2 years. To minimise potential bias and loss of power, we considered the data only from participants who were followed more than three times. Consequently, 6622 participants (3441 men and 3181 women; 3569 Ansan and 3053 Ansung) were identified, for whom both genotypic and phenotypic information were available. For the genotype data, 3,352,722 SNPs were considered in the analysis after quality control and imputation procedures. The details of the procedures are included in Supplementary Text S1. All participants provided written informed consent, and the ethics committee of the Korean Centre for Disease Control and institutional review boards of the Korea University Ansan Hospital approved the study (IRB no. 2020AS0356)
4.2. Lung Function
We considered pre-bronchodilator absolute values, which included FVC, FEV1, FEV1/FVC ratio, the average forced expiratory flow during the mid (25–75%) portion of the FVC (FEF 25–75%), and maximal voluntary ventilation (MVV). To investigate the exact COPD prevalence as well as bronchodilator response in a sub-population, we also performed a post-bronchodilator test for FVC, FEV1, and FEV1/FVC ratio, after short-acting bronchodilator treatment. These traits were measured by technicians using a portable spirometer (Vmax-2130, Sensor Medics, Yorba Linda, CA, USA) according to standardised protocols of the American Thoracic Society [55]. All participants performed the pre-bronchodilator spirometry test until at least three repeated measurements were completed; an acceptable measure was determined when the differences between the largest and the next largest FVC and FEV1 values were within 0.15 L [55]. Additionally, a post-bronchodilator test was performed in participants from Ansan-cohort only. Fenoterol hydrobromide (400 μg, 2001–2008) or salbutamol (400 μg, 2009 onwards) was administered and each test was repeated after a 15 to 30 min delay. Calibration and quality control of spirometric examinations were performed regularly based on the American Thoracic Society guidelines [55,56,57].
4.3. Epidemiological Concept of and
We assumed that the trait of subject i at time point j was yij, and yij was a subject-specific function of the subject’s age, ageij. For our data, j was a number from 1 to 8. The subject-specific function was fi and was a measurement error. Then we assumed the following:
(1)
Based on this definition, we estimated the subject-specific means (which indicate the mean of lung function traits across all visits) and annual change rates for the trait; two different SNP-based heritability parameters, and , were also defined for both (see Supplementary Text S2).
If is high, then the overall genetic effect on the cross-sectional average tends to be high, and is equivalent to the SNP heritability estimated using genome-wide complex trait analysis (GCTA) [24]. The average annual change affected by genetic components is indicated by > 0. Supplementary Figure S1 illustrates the practical concept for and .
4.4. Estimation of and
A two-stage method was used to calculate and based on the method of Yang et al. [24]. First, we fitted a simple linear regression model for the subjects during the same visiting time, adjusting for age for each trait, as in Equation (2), which was transformed from Equation (1) (see Supplementary Text S2). Each participant was measured up to eight times and participants with at least three measurements were analysed. Since variances in lung function traits were heterogeneous at the eight different time points, the residual variances for each time point were calculated from linear regressions after adjusting for age, sex, and height for absolute values. The inverse of the residual variances for trait k and time point j was taken as wjk, and the weighted linear regression was performed for each subject i using the following equation:
(2)
Here, indicates the mean age at the observed time points. indicates the subject-specific mean of subject i for trait k at , and indicates the annual change rate. Since we subtracted the mean age from each time point, the estimate of and is orthogonal. Since we considered several time points in the analysis, the subject-specific mean was more important than the baseline value. We then applied a rank-based inverse normal transformation to and , and the bivariate linear mixed model was applied to estimate and after adjusting for sex, mean age, and mean height. The genetic relationship matrix (GRM) between pairs of individuals from the genome-wide variants were applied in estimation of variance covariance matrix in the linear mixed model (Supplementary Text S2), and the PC plot, based on GRM, was also generated (Supplementary Figure S6). This model was estimated using the “--reml-bivar” option of GCTA [58]. The genetic correlation between and was also generated using this model and option. If the subject-specific mean () is smaller, the positive genetic correlation indicates the change rate (), which is a negative value, getting smaller, and vice versa.
4.5. Subgroup Analysis by Smoking Status
To conduct subgroup analyses, we separated subjects into ever- and never-smoker groups. The ever-smoking group included past and current smokers. The heritability of subject-specific means and annual change rates were estimated using GCTA for the two groups.
4.6. Estimating Heritability Attributed to SNP–Smoking Interaction
To estimate the variance in the SNP–smoking interaction for both subject-specific means and annual change rates, we applied the genotype-environment (GxE) interaction model, in which the main effects of environmental factors were included as fixed effects, while GxE interaction effects were treated as random effects (using the “--gxe” option of GCTA).
4.7. Genome-Wide Association Studies
Genome-wide association studies were conducted for eight lung function traits. Quality control and imputation were conducted, and the detailed procedure is described in Supplementary Text S1. Linear regression for and was conducted after adjusting for age, sex, height, and principle components (PCs, PC1 to PC10). The analyses were conducted using ONETOOL [59] and PLINK [60].
All authors contributed to this research study. S.W. developed the study concept, D.L. and S.W. developed the study design, S.K., S.L. (Seungku Lee), D.Y. and C.S. were responsible for patient recruitment and data acquisition. D.L. and S.W. conducted statistical analyses. D.L., S.W. and W.K. (Woojin Kim) wrote the manuscript. S.L. (Sanghun Lee), J.A., A.D., W.K. (Wonji Kim), K.L., S.H., E.K.S. and M.C. contributed to the critical revision of the manuscript. All authors have read and agreed to the published version of the manuscript.
All participants provided written informed consent and the ethics committee of the Korean Center for Disease Control, and institutional review boards of the Korea University Ansan Hospital approved the study (IRB no. 2020AS0356).
Not applicable.
The information of the datasets analysed during the current study is available in the KoGES website [
The authors thank Jungsoo Gim, Sangchul Park and Kyungtaek Park for their support in different aspects of the study.
Edwin K. Silverman has received institutional grant support from GSK and Bayer. Michael H. Cho has received institutional grant support from GSK and Bayer and consulting or speaking fees from AstraZeneca and Illumina.
COPD | (Chronic obstructive pulmonary disease) |
GWAS | (genome-wide association study) |
FEV1 | (forced expiratory volume in 1 s) |
FVC | (forced vital capacity) |
FEF | 25–75% (average forced expiratory flow during the mid (25–75%) portion of the FVC) |
MVV | (SNP (single nucleotide polymorphism) |
KoGES | (Korean Genome and Epidemiology Study) |
GCTA | (genome-wide complex trait analysis) |
GxE | (genotype-environment interaction) |
PC | (principle components) |
FDR | (false-discovery rate) |
DSL | (disease susceptibility loci) |
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure 1. Number of samples used for analyzing the eight lung function traits and the mean observed values across eight periods. (a) Numbers of samples in each period (p1–p8). (b,c) Profile plot of the mean values of the eight traits in each period. Since the ranges of the traits differed, they are presented as two separate plots. FVC, forced vital capacity; FEV1, forced expiratory volume in 1 s; FEV1/FVC, forced expiratory volume in 1 s to forced vital capacity ratio; FEF 25–75%, the average forced expiratory flow during the mid (25–75%) portion of the FVC; MVV, maximal voluntary ventilation; post-FVC, post-bronchodilator test for FVC; post-FEV1, post-bronchodilator test for FEV1; post- FEV1/FVC, post-bronchodilator test for FEV1/ FVC ratio.
Figure 2. Single nucleotide polymorphism (SNP) heritability of eight lung function traits. SNP heritability of subject-specific means (a) and the annual change rates (b) in eight lung function traits. Error bars correspond to standard error. The black dots on the bars represent p-values. The black dashed line indicates the 0.05 significance level. Red dots indicate significant findings at a false-discovery rate of 0.05. FVC, forced vital capacity; FEV1, forced expiratory volume in 1 s; FEV1/FVC, forced expiratory volume in 1 s to forced vital capacity ratio; FEF 25–75%, the average forced expiratory flow during the mid (25–75%) portion of the FVC; MVV, maximal voluntary ventilation; post-FVC, post-bronchodilator test for FVC; post- FEV1, post-bronchodilator test for FEV1; post- FEV1/FVC, post-bronchodilator test for FEV1/ FVC ratio.
Figure 3. Single nucleotide polymorphism (SNP) heritability of eight lung function traits in never- and ever-smoker groups. SNP heritability of subject-specific means (a) and the annual change rates (b). Error bars correspond to standard error. The black dot on the bars represent p-values. The black dashed line indicates the 0.05 significance level. Red dots indicate significant findings at a false-discovery rate of 0.05. FVC, forced vital capacity; FEV1, forced expiratory volume in 1 s; FEV1/FVC, forced expiratory volume in 1 s to forced vital capacity ratio; FEF 25–75%, the average forced expiratory flow during the mid (25–75%) portion of the FVC; MVV, maximal voluntary ventilation; post-FVC, post-bronchodilator test for FVC; post-FEV1, post-bronchodilator test for FEV1; post- FEV1/FVC, post-bronchodilator test for FEV1/FVC ratio.
Descriptive characteristics of study participants with more than three follow-up visits.
Characteristic | ALL | Never |
Ever |
p-Value |
---|---|---|---|---|
Sample Size, n | 6622 | 3009 | 2095 | - |
Female, n (%) | 3181 (48.04%) | 2588 (86%) | 104 (5%) | <2.2 × 10−16 |
Age, year (means ± SD) | 51.51 ± 8.52 | 51.19 ± 8.3 | 50.54 ± 7.94 | 0.005 |
Height, cm (means ± SD) | 160.3 ± 8.57 | 155.9 ± 7.02 | 166.6 ± 6.26 | <2.2 × 10−16 |
The “All” status included 5104 participants with a smoking history, and these participants were considered in the smoking subgroup analysis and SNP–smoking interaction analysis. Some participants had no smoking status. p-values were generated from statistical tests comparing never-smoker and ever-smoker groups. The Chi-square test was used for females and a t-test was performed for age and height. The sample sizes reported here are the number of samples used in the heritability estimation.
Genetic correlation of subject-specific means (
Trait |
|
|
|
|
|
|
|
|
---|---|---|---|---|---|---|---|---|
FEV1/FVC (%) | 0.1807 | 0.1538 | 0.1047 | 0.6279 | 0.1466 | 0.0000459 | 0.117 | 0.0121 |
post-FEV1/FVC (%) | 0.2593 | 0.1742 | 0.099 | 0.466 | 0.2156 | 0.0219 | 0.155 | 0.0165 |
Heritability of SNP–smoking interaction for annual change rates with significant p-value (p < 0.05).
Traits | Subject-Specific Means | Annual Change Rates | ||||
---|---|---|---|---|---|---|
|
s.e. | p-Value |
|
s.e. | p-Value | |
post-FEV1/FVC (%) | 0.206 | 0.191 | 1.37 × 10−1 | 0.402 | 0.199 | 2.12 × 10−2 |
SNP, single nucleotide polymorphism;
Supplementary Materials
The following supporting information can be downloaded at:
References
1. Miller, M.R.; Hankinson, J.; Brusasco, V.; Burgos, F.; Casaburi, R.; Coates, A.; Crapo, R.; Enright, P.; Van Der Grinten, C.P.M.; Gustafsson, P. et al. Standardisation of spirometry. Eur. Respir. J.; 2005; 26, pp. 319-338. [DOI: https://dx.doi.org/10.1183/09031936.05.00034805] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16055882]
2. Calverley, P.M.; Walker, P. Chronic obstructive pulmonary disease. Lancet; 2003; 362, pp. 1053-1061. [DOI: https://dx.doi.org/10.1016/S0140-6736(03)14416-9]
3. Halbert, R.J.; Natoli, J.L.; Gano, A.; Badamgarav, E.; Buist, A.S.; Mannino, D.M. Global burden of COPD: Systematic review and meta-analysis. Eur. Respir. J.; 2006; 28, pp. 523-532. [DOI: https://dx.doi.org/10.1183/09031936.06.00124605] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16611654]
4. Ferrer, M.; Alonso, J.; Morera, J.; Marrades, R.M.; Khalaf, A.; Aguar, M.C.; Plaza, V.; Prieto, L.; Antó, J.M. Chronic obstructive pulmonary disease stage and health-related quality of life. The Quality of Life of Chronic Obstructive Pulmonary Disease Study Group. Ann. Intern. Med.; 1997; 127, 1072. [DOI: https://dx.doi.org/10.7326/0003-4819-127-12-199712150-00003] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9412309]
5. Astemborski, J.A.; Beaty, T.H.; Cohen, B.H.; Opitz, J.M.; Reynolds, J.F. Variance components analysis of forced expiration in families. Am. J. Med. Genet.; 1985; 21, pp. 741-753. [DOI: https://dx.doi.org/10.1002/ajmg.1320210417] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/4025399]
6. Wilk, J.B.; Djousse, L.; Arnett, D.K.; Rich, S.S.; Province, M.A.; Hunt, S.C.; Crapo, R.O.; Higgins, M.; Myers, R.H. Evidence for major genes influencing pulmonary function in the NHLBI Family Heart Study. Genet. Epidemiol.; 2000; 19, pp. 81-94. [DOI: https://dx.doi.org/10.1002/1098-2272(200007)19:1<81::AID-GEPI6>3.0.CO;2-8]
7. Palmer, L.; Knuiman, M.; Divitini, M.; Burton, P.; James, A.; Bartholomew, H.; Ryan, G.; Musk, A. Familial aggregation and heritability of adult lung function: Results from the Busselton Health Study. Eur. Respir. J.; 2001; 17, pp. 696-702. [DOI: https://dx.doi.org/10.1183/09031936.01.17406960]
8. Lewitter, F.I.; Tager, I.B.; McGue, M.; Tishler, P.V.; Speizer, F.E.; Lewiiter, F.I. Genetic and environmental determinants of level of pulmonary function. Am. J. Epidemiol.; 1984; 120, pp. 518-530. [DOI: https://dx.doi.org/10.1093/oxfordjournals.aje.a113912]
9. Zhou, J.J.; Cho, M.H.; Castaldi, P.J.; Hersh, C.P.; Silverman, E.K.; Laird, N.M. Heritability of Chronic Obstructive Pulmonary Disease and Related Phenotypes in Smokers. Am. J. Respir. Crit. Care Med.; 2013; 188, pp. 941-947. [DOI: https://dx.doi.org/10.1164/rccm.201302-0263OC]
10. Wyss, A.B.; Sofer, T.; Lee, M.K.; Terzikhan, N.; Nguyen, J.N.; LaHousse, L.; Latourelle, J.C.; Smith, A.V.; Bartz, T.M.; Feitosa, M.F. et al. Multiethnic meta-analysis identifies ancestry-specific and cross-ancestry loci for pulmonary function. Nat. Commun.; 2018; 9, 2976. [DOI: https://dx.doi.org/10.1038/s41467-018-05369-0]
11. Soler Artigas, M.; Loth, D.W.; Wain, L.V.; Gharib, S.A.; Obeidat, M.; Tang, W.; Zhai, G.; Zhao, J.H.; Smith, A.V.; Huffman, J.E. et al. Genome-wide association and large-scale follow up identifies 16 new loci influencing lung function. Nat. Genet.; 2011; 43, pp. 1082-1090. [DOI: https://dx.doi.org/10.1038/ng.941]
12. Hancock, D.; Eijgelsheim, M.; Wilk, J.B.; Gharib, S.A.; Loehr, L.; Marciante, K.D.; Franceschini, N.; A Van Durme, Y.M.T.; Chen, T.-H.; Barr, R.G. et al. Meta-analyses of genome-wide association studies identify multiple loci associated with pulmonary function. Nat. Genet.; 2010; 42, pp. 45-52. [DOI: https://dx.doi.org/10.1038/ng.500]
13. Repapi, E.; Sayers, I.; Wain, L.V.; Burton, P.R.; Johnson, T.; Obeidat, M.; Zhao, J.H.; Ramasamy, A.; Zhai, G.; Vitart, V. et al. Genome-wide association study identifies five loci associated with lung function. Nat. Genet.; 2010; 42, pp. 36-44. [DOI: https://dx.doi.org/10.1038/ng.501]
14. Wain, L.V.; Shrine, N.; Miller, S.; E Jackson, V.; Ntalla, I.; Artigas, M.S.; Billington, C.K.; Kheirallah, A.K.; Allen, R.; Cook, J.P. et al. Novel insights into the genetics of smoking behaviour, lung function, and chronic obstructive pulmonary disease (UK BiLEVE): A genetic association study in UK Biobank. Lancet Respir. Med.; 2015; 3, pp. 769-781. [DOI: https://dx.doi.org/10.1016/S2213-2600(15)00283-0]
15. van der Plaat, D.A.; de Jong, K.; Lahousse, L.; Faiz, A.; Vonk, J.M.; van Diemen, C.C.; Nedeljkovic, I.; Amin, N.; Brusselle, G.G.; Hofman, A. et al. Genome-wide association study on the FEV1/FVC ratio in never-smokers identifies HHIP and FAM13A. J. Allergy Clin. Immunol.; 2017; 139, pp. 533-540. [DOI: https://dx.doi.org/10.1016/j.jaci.2016.06.062]
16. Park, B.; Koo, S.-M.; An, J.; Lee, M.; Kang, H.Y.; Qiao, D.; Cho, M.H.; Sung, J.; Silverman, E.K.; Yang, H.-J. et al. Genome-wide assessment of gene-by-smoking interactions in COPD. Sci. Rep.; 2018; 8, 9319. [DOI: https://dx.doi.org/10.1038/s41598-018-27463-5]
17. Gottlieb, D.J.; Wilk, J.B.; Harmon, M.; Evans, J.C.; Joost, O.; Levy, D.; O’Connor, G.T.; Myers, R.H. Heritability of Longitudinal Change in Lung Function. The Framingham study. Am. J. Respir. Crit. Care Med.; 2001; 164, pp. 1655-1659. [DOI: https://dx.doi.org/10.1164/ajrccm.164.9.2010122]
18. John, C.; Artigas, M.S.; Hui, J.; Nielsen, S.F.; Rafaels, N.; Paré, P.D.; Hansel, N.N.; Shrine, N.; Kilty, I.; Malarstig, A. et al. Genetic variants affecting cross-sectional lung function in adults show little or no effect on longitudinal lung function decline. Thorax; 2017; 72, pp. 400-408. [DOI: https://dx.doi.org/10.1136/thoraxjnl-2016-208448]
19. Li, D.; Kang, H.; Lee, S.; Won, S. Progressive effects of single-nucleotide polymorphisms on 16 phenotypic traits based on longitudinal data. Genes Genom.; 2020; 42, pp. 393-403. [DOI: https://dx.doi.org/10.1007/s13258-019-00902-x]
20. Terzikhan, N.; Sun, F.; Verhamme, F.M.; Adams, H.H.; Loth, D.; Bracke, K.; Stricker, B.H.; Lahousse, L.; Dupuis, J.; Brusselle, G. et al. Heritability and genome-wide association study of diffusing capacity of the lung. Eur. Respir. J.; 2018; 52, 1800647. [DOI: https://dx.doi.org/10.1183/13993003.00647-2018]
21. Robinson, M.R.; English, G.; Moser, G.; Lloyd-Jones, L.R.; Triplett, M.A.; Zhu, Z.; Nolte, I.M.; van Vliet-Ostaptchouk, J.; Snieder, H. LifeLines Cohort Studyet al. Genotype-covariate interaction effects and the heritability of adult body mass index. Nat. Genet.; 2017; 49, pp. 1174-1181. [DOI: https://dx.doi.org/10.1038/ng.3912]
22. Yang, J.; Bakshi, A.; Zhu, Z.; Hemani, G.; Vinkhuyzen, A.; Lee, S.H.; Robinson, M.; Perry, J.R.B.; Nolte, I.M. The LifeLines Cohort Studyet al. Genetic variance estimation with imputed variants finds negligible missing heritability for human height and body mass index. Nat. Genet.; 2015; 47, pp. 1114-1120. [DOI: https://dx.doi.org/10.1038/ng.3390]
23. Cho, Y.S.; Go, M.J.; Kim, Y.J.; Heo, J.Y.; Oh, J.H.; Ban, H.-J.; Yoon, D.; Lee, M.H.; Kim, D.-J.; Park, M. et al. A large-scale genome-wide association study of Asian populations uncovers genetic factors influencing eight quantitative traits. Nat. Genet.; 2009; 41, pp. 527-534. [DOI: https://dx.doi.org/10.1038/ng.357]
24. Yang, J.; Lee, S.H.; Goddard, M.E.; Visscher, P.M. GCTA: A Tool for Genome-wide Complex Trait Analysis. Am. J. Hum. Genet.; 2011; 88, pp. 76-82. [DOI: https://dx.doi.org/10.1016/j.ajhg.2010.11.011]
25. Vestbo, J.; Hurd, S.S.; Agustí, A.G.; Jones, P.W.; Vogelmeier, C.; Anzueto, A.; Barnes, P.J.; Fabbri, L.M.; Martinez, F.J.; Nishimura, M. et al. Global strategy for the diagnosis, management, and prevention of chronic obstructive pulmonary disease: GOLD executive summary. Am. J. Respir. Crit. Care Med.; 2013; 187, pp. 347-365. [DOI: https://dx.doi.org/10.1164/rccm.201204-0596PP]
26. Hancock, D.B.; Artigas, M.S.; Gharib, S.A.; Henry, A.; Manichaikul, A.; Ramasamy, A.; Loth, D.; Imboden, M.; Koch, B.; McArdle, W.L. et al. Genome-Wide Joint Meta-Analysis of SNP and SNP-by-Smoking Interaction Identifies Novel Loci for Pulmonary Function. PLoS Genet.; 2012; 8, e1003098. [DOI: https://dx.doi.org/10.1371/journal.pgen.1003098]
27. Sakornsakolpat, P.; Prokopenko, D.; Lamontagne, M.; Reeve, N.F.; Guyatt, A.L.; Jackson, V.E.; Shrine, N.; Qiao, D.; Bartz, T.M.; Kim, D.K. et al. Genetic landscape of chronic obstructive pulmonary disease identifies heterogeneous cell-type and phenotype associations. Nat. Genet.; 2019; 51, pp. 494-505. [DOI: https://dx.doi.org/10.1038/s41588-018-0342-2]
28. Hardin, M.; Cho, M.H.; McDonald, M.-L.; Wan, E.; Lomas, D.A.; Coxson, H.O.; MacNee, W.; Vestbo, J.; Yates, J.C.; Agusti, A. et al. A genome-wide analysis of the response to inhaled β2-agonists in chronic obstructive pulmonary disease. Pharm. J.; 2016; 16, pp. 326-335. [DOI: https://dx.doi.org/10.1038/tpj.2015.65]
29. Lutz, S.M.; Cho, M.H.; Young, K.; Hersh, C.P.; Castaldi, P.J.; McDonald, M.-L.; A Regan, E.; Mattheisen, M.; DeMeo, D.L.; Parker, M. et al. A genome-wide association study identifies risk loci for spirometric measures among smokers of European and African ancestry. BMC Genet.; 2015; 16, 138. [DOI: https://dx.doi.org/10.1186/s12863-015-0299-4]
30. Kichaev, G.; Bhatia, G.; Loh, P.-R.; Gazal, S.; Burch, K.; Freund, M.K.; Schoech, A.; Pasaniuc, B.; Price, A.L. Leveraging Polygenic Functional Enrichment to Improve GWAS Power. Am. J. Hum. Genet.; 2019; 104, pp. 65-75. [DOI: https://dx.doi.org/10.1016/j.ajhg.2018.11.008]
31. McKay, J.D. SpiroMeta Consortium Hung, R.J.; Han, Y.; Zong, X.; Carreras-Torres, R.; Christiani, D.C.; Caporaso, N.E.; Johansson, M.; Xiao, X. et al. Large-scale association analysis identifies new lung cancer susceptibility loci and heterogeneity in genetic susceptibility across histological subtypes. Nat. Genet.; 2017; 49, pp. 1126-1132. [DOI: https://dx.doi.org/10.1038/ng.3892] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28604730]
32. DeMeo, D.L.; Carey, V.J.; Chapman, H.A.; Reilly, J.J.; Ginns, L.C.; Speizer, F.E.; Weiss, S.T.; Silverman, E.K. Familial aggregation of FEF(25-75) and FEF(25-75)/FVC in families with severe, early onset COPD. Thorax; 2004; 59, pp. 396-400. [DOI: https://dx.doi.org/10.1136/thx.2003.012856] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15115866]
33. Agustí, A.; Hogg, J.C. Update on the Pathogenesis of Chronic Obstructive Pulmonary Disease. N. Engl. J. Med.; 2019; 381, pp. 1248-1256. [DOI: https://dx.doi.org/10.1056/NEJMra1900475] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31553836]
34. Ross, J.C.; Castaldi, P.J.; Cho, M.H.; Hersh, C.P.; Rahaghi, F.N.; Sanchez-Ferrero, G.V.; Parker, M.M.; Litonjua, A.A.; Sparrow, D.; Dy, J.G. et al. Longitudinal Modeling of Lung Function Trajectories in Smokers with and without Chronic Obstructive Pulmonary Disease. Am. J. Respir. Crit. Care Med.; 2018; 198, pp. 1033-1042. [DOI: https://dx.doi.org/10.1164/rccm.201707-1405OC]
35. Lange, P.; Celli, B.R.; Agustí, A.; Jensen, G.B.; Divo, M.; Faner, R.; Guerra, S.; Marott, J.L.; Martinez, F.D.; Martinez-Camblor, P. et al. Lung-Function Trajectories Leading to Chronic Obstructive Pulmonary Disease. N. Engl. J. Med.; 2015; 373, pp. 111-122. [DOI: https://dx.doi.org/10.1056/NEJMoa1411532]
36. McGeachie, M.J.; Yates, K.P.; Zhou, X.; Guo, F.; Sternberg, A.L.; Van Natta, M.L.; Wise, R.A.; Szefler, S.J.; Sharma, S.; Kho, A.T. et al. Patterns of Growth and Decline in Lung Function in Persistent Childhood Asthma. N. Engl. J. Med.; 2016; 374, pp. 1842-1852. [DOI: https://dx.doi.org/10.1056/NEJMoa1513737]
37. Hogg, J.C.; Chu, F.; Utokaparch, S.; Woods, R.; Elliott, W.M.; Buzatu, L.; Cherniack, R.M.; Rogers, R.M.; Sciurba, F.C.; Coxson, H.O. et al. The Nature of Small-Airway Obstruction in Chronic Obstructive Pulmonary Disease. N. Engl. J. Med.; 2004; 350, pp. 2645-2653. [DOI: https://dx.doi.org/10.1056/NEJMoa032158]
38. Salvi, S.S.; Barnes, P.J. Chronic obstructive pulmonary disease in non-smokers. Lancet; 2009; 374, pp. 733-743. [DOI: https://dx.doi.org/10.1016/S0140-6736(09)61303-9]
39. Olloquequi, J.; Jaime, S.; Parra, V.; Cornejo-Córdova, E.; Valdivia, G.; Agusti, A.; Silva, R.O. Correction to: Comparative analysis of COPD associated with tobacco smoking, biomass smoke exposure or both. Respir. Res.; 2018; 19, 77. [DOI: https://dx.doi.org/10.1186/s12931-018-0765-4]
40. Agustí, A.; Noell, G.; Brugada, J.; Faner, R. Lung function in early adulthood and health in later life: A transgenerational cohort analysis. Lancet Respir. Med.; 2017; 5, pp. 935-945. [DOI: https://dx.doi.org/10.1016/S2213-2600(17)30434-4]
41. Kim, W.J.; Lee, M.K.; Shin, C.; Cho, N.H.; Lee, S.D.; Oh, Y.-M.; Sung, J. Genome-wide association studies identify locus on 6p21 influencing lung function in the Korean population. Respirology; 2014; 19, pp. 360-368. [DOI: https://dx.doi.org/10.1111/resp.12230]
42. Jiang, S.S.; Fang, W.-T.; Hou, Y.-H.; Huang, S.-F.; Yen, B.L.; Chang, J.-L.; Li, S.-M.; Liu, H.-P.; Liu, Y.-L.; Huang, C.-T. et al. Upregulation of SOX9 in Lung Adenocarcinoma and Its Involvement in the Regulation of Cell Growth and Tumorigenicity. Clin. Cancer Res.; 2010; 16, pp. 4363-4373. [DOI: https://dx.doi.org/10.1158/1078-0432.CCR-10-0138]
43. Cho, M.H.; Boutaoui, N.; Klanderman, B.J.; Sylvia, J.S.; Ziniti, J.P.; Hersh, C.P.; DeMeo, D.L.; Hunninghake, G.M.; Litonjua, A.; Sparrow, D. et al. Variants in FAM13A are associated with chronic obstructive pulmonary disease. Nat. Genet.; 2010; 42, pp. 200-202. [DOI: https://dx.doi.org/10.1038/ng.535]
44. Fingerlin, T.E.; Murphy, E.; Zhang, W.; Peljto, A.L.; Brown, K.K.; Steele, M.P.; Loyd, J.E.; Cosgrove, G.P.; Lynch, D.; Groshong, S. et al. Genome-wide association study identifies multiple susceptibility loci for pulmonary fibrosis. Nat Genet.; 2013; 45, pp. 613-620. [DOI: https://dx.doi.org/10.1038/ng.2609]
45. Kim, W.J.; Lim, M.N.; Hong, Y.; Silverman, E.K.; Lee, J.-H.; Jung, B.H.; Ra, S.W.; Choi, H.S.; Jung, Y.J.; Park, Y.B. et al. Association of Lung Function Genes with Chronic Obstructive Pulmonary Disease. Lung; 2014; 192, pp. 473-480. [DOI: https://dx.doi.org/10.1007/s00408-014-9579-4]
46. Xu, H.; Xu, G.; Wang, D.; Zheng, C.; Wan, L. Molecular cloning and tissue distribution of the phosphotyrosine interaction domain containing 1 (PID1) gene in Tianfu goat. Gene; 2013; 515, pp. 71-77. [DOI: https://dx.doi.org/10.1016/j.gene.2012.11.025]
47. Zhao, Y.; Zhang, C.; Chen, X.; Gao, C.; Ji, C.; Chen, F.; Zhu, C.; Zhu, J.; Wang, J.; Qian, L. et al. Overexpression of NYGGF4 (PID1) induces mitochondrial impairment in 3T3-L1 adipocytes. Mol. Cell. Biochem.; 2010; 340, pp. 41-48. [DOI: https://dx.doi.org/10.1007/s11010-010-0398-5]
48. Kirkham, P.A.; Barnes, P.J. Oxidative Stress in COPD. Chest; 2013; 144, pp. 266-273. [DOI: https://dx.doi.org/10.1378/chest.12-2664]
49. Artigas, M.S.; Wain, L.V.; Shrine, N.; McKeever, T.M.; Sayers, I.; Hall, I.; Tobin, M.D.; Bileve, U. Targeted Sequencing of Lung Function Loci in Chronic Obstructive Pulmonary Disease Cases and Controls. PLoS ONE; 2017; 12, e0170222. [DOI: https://dx.doi.org/10.1371/journal.pone.0170222]
50. Lee, K.-S.; Kim, K.H.; Oh, Y.-M.; Han, B.; Kim, W.J. A genome wide association study for lung function in the Korean population using an exome array. Korean J. Intern. Med.; 2021; 36, (Suppl. S1), pp. S142-S150. [DOI: https://dx.doi.org/10.3904/kjim.2019.204]
51. Siller, S.S.; Sharma, H.; Li, S.; Yang, J.; Zhang, Y.; Holtzman, M.J.; Winuthayanon, W.; Colognato, H.; Holdener, B.C.; Li, F.-Q. et al. Conditional knockout mice for the distal appendage protein CEP164 reveal its essential roles in airway multiciliated cell differentiation. PLoS Genet.; 2017; 13, e1007128. [DOI: https://dx.doi.org/10.1371/journal.pgen.1007128]
52. Yaghi, A.; Dolovich, M.B. Airway Epithelial Cell Cilia and Obstructive Lung Disease. Cells; 2016; 5, 40. [DOI: https://dx.doi.org/10.3390/cells5040040]
53. Tilley, A.E.; Walters, M.S.; Shaykhiev, R.; Crystal, R.G. Cilia dysfunction in lung disease. Annu. Rev. Physiol.; 2015; 77, pp. 379-406. [DOI: https://dx.doi.org/10.1146/annurev-physiol-021014-071931]
54. Kim, Y.; Han, B.G. KoGES Group. Cohort Profile: The Korean Genome and Epidemiology Study (KoGES) Consortium. Int. J. Epidemiol.; 2017; 46, 1350. [DOI: https://dx.doi.org/10.1093/ije/dyx105] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28938752]
55. Graham, B.L.; Steenbruggen, I.; Miller, M.R.; Barjaktarevic, I.Z.; Cooper, B.G.; Hall, G.L.; Thompson, B.R. Standardization of Spirometry, 1994 Update. American Thoracic Society. Am. J. Respir. Crit. Care Med.; 1995; 152, pp. 1107-1136.
56. Kim, S.; Kim, H.; Cho, N.; Lee, S.K.; Han, B.-G.; Sull, J.W.; Jee, S.H.; Shin, C. Identification of FAM13A gene associated with the ratio of FEV1 to FVC in Korean population by genome-wide association studies including gene—Environment interactions. J. Hum. Genet.; 2015; 60, pp. 139-145. [DOI: https://dx.doi.org/10.1038/jhg.2014.118] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25608829]
57. Shin, C.; Kim, J.; Kim, J.; Lee, S.; Shim, J.; In, K.; Kang, K.; Yoo, S.; Cho, N.; Kimm, K. et al. Association of Habitual Snoring with Glucose and Insulin Metabolism in Nonobese Korean Adult Men. Am. J. Respir. Crit. Care Med.; 2005; 171, pp. 287-291. [DOI: https://dx.doi.org/10.1164/rccm.200407-906OC] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15542791]
58. Yang, J.; A Manolio, T.; Pasquale, L.R.; Boerwinkle, E.; Caporaso, N.; Cunningham, J.M.; de Andrade, M.; Feenstra, B.; Feingold, E.; Hayes, M.G. et al. Genome partitioning of genetic variation for complex traits using common SNPs. Nat. Genet.; 2011; 43, pp. 519-525. [DOI: https://dx.doi.org/10.1038/ng.823]
59. Song, Y.E.; Lee, S.; Park, K.; Elston, R.C.; Yang, H.-J.; Won, S. ONETOOL for the analysis of family-based big data. Bioinformatics; 2018; 34, pp. 2851-2853. [DOI: https://dx.doi.org/10.1093/bioinformatics/bty180]
60. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.W.; Daly, M.J. et al. PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am. J. Hum. Genet.; 2007; 81, pp. 559-575. [DOI: https://dx.doi.org/10.1086/519795]
61. Delaneau, O.; Marchini, J.; Zagury, J.F. A linear complexity phasing method for thousands of genomes. Nat. Methods; 2011; 9, pp. 179-181. [DOI: https://dx.doi.org/10.1038/nmeth.1785] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22138821]
62. Howie, B.N.; Donnelly, P.; Marchini, J. A flexible and accurate genotype imputation method for the next generation of genome-wide association studies. PLoS Genet.; 2009; 5, e1000529. [DOI: https://dx.doi.org/10.1371/journal.pgen.1000529] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19543373]
63. Howie, B.; Marchini, J.; Stephens, M. Genotype imputation with thousands of genomes. G3 (Bethesda); 2011; 1, pp. 457-470. [DOI: https://dx.doi.org/10.1534/g3.111.001198] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22384356]
64. Lee, S.H.; Yang, J.; Goddard, M.E.; Visscher, P.M.; Wray, N.R. Estimation of pleiotropy between complex diseases using single-nucleotide polymorphism-derived genomic relationships and restricted maximum likelihood. Bioinformatics; 2012; 28, pp. 2540-2542. [DOI: https://dx.doi.org/10.1093/bioinformatics/bts474]
65. Yang, J.; Benyamin, B.; McEvoy, B.P.; Gordon, S.; Henders, A.K.; Nyholt, D.R.; Madden, P.A.; Heath, A.C.; Martin, N.G.; Montgomery, G.W. et al. Common SNPs explain a large proportion of the heritability for human height. Nat. Genet.; 2010; 42, pp. 565-569. [DOI: https://dx.doi.org/10.1038/ng.608]
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Genetic influence on lung functions has been identified in previous studies; however, the relative longitudinal effects of genetic factors and their interactions with smoking on lung function remain unclear. Here, we identified the longitudinal effects of genetic variants on lung function by determining single nucleotide polymorphism (SNP) heritability and genetic correlations, and by analyzing interactions with smoking. Subject-specific means and annual change rates were calculated for eight spirometric measures obtained from 6622 Korean adults aged 40–69 years every two years for 14 years, and their heritabilities were estimated separately. Statistically significant (p < 0.05) heritability for the subject-specific means of all spirometric measures (8~32%) and change rates of forced expiratory volume in 1 s to forced vital capacity ratio (FEV1/FVC; 16%) and post-bronchodilator FEV1/FVC (17%) were detected. Significant genetic correlations of the change rate with the subject-specific mean were observed for FEV1/FVC (
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details








1 Interdisciplinary Program in Bioinformatics, Seoul National University, Seoul 08826, Korea;
2 Department of Internal Medicine and Environmental Health Center, School of Medicine, Kangwon National University, Chuncheon 24341, Korea;
3 Department of Public Health Science, Graduate School of Public Health, Seoul National University, Seoul 08826, Korea;
4 Institute for Human Genomic Study, College of Medicine, Korea University, Seoul 136701, Korea;
5 Interdisciplinary Program in Bioinformatics, Seoul National University, Seoul 08826, Korea;
6 Channing Division of Network Medicine and Division of Pulmonary and Critical Care Medicine, Brigham and Women’s Hospital, Boston, MA 02115, USA;
7 Department of Medical Consilience, Graduate School, Dankook University, Yongin 16890, Korea;
8 Division of Allergy and Respiratory Disease Research, Department of Chronic Disease Convergence Research, National Institute of Health, Korea Disease Control and Prevention Agency, Cheongju 28159, Korea;
9 Korea Medical Institute, Seoul 03173, Korea;
10 Institute for Human Genomic Study, College of Medicine, Korea University, Seoul 136701, Korea;
11 Interdisciplinary Program in Bioinformatics, Seoul National University, Seoul 08826, Korea;