Effects of testosterone and sex hormone binding globulin on lung function in males and females: a multivariable Mendelian Randomisation study.

BACKGROUND: Observational studies suggest that total testosterone (TT) and sex hormone-binding globulin (SHBG) may have beneficial effects on lung function, but these findings might be spurious due to confounding and reverse causation. We addressed these limitations by using multivariable Mendelian randomisation (MVMR) to investigate the independent causal effects of TT and SHBG on lung function. METHODS: We first identified genetic instruments by performing genome-wide association analyses of TT and SHBG in the large UK Biobank, separately in males and females. We then assessed the independent effects of TT and SHBG on forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC) and FEV1/FVC using one-sample MVMR. We addressed pleiotropy, which could bias MVMR, using several methods that account for it. We performed subgroup MVMR analyses by obesity, physical activity and menopausal status, and assessed associations between TT and SHBG with lung function decline. Finally, we compared the MVMR results with those of observational analyses in the UK Biobank. FINDINGS: In the MVMR analyses, there was evidence of pleiotropy, but results were consistent when accounting for it. We found a strong beneficial effect of TT on FVC and FEV1 in both males and females, but a moderate detrimental effect of SHBG on FEV1 and FEV1/FVC in males only. Subgroup analyses suggested stronger effects of TT among obese and older males. The observational analyses, in line with previous studies, agreed with MRMV for TT, but not for SHBG. INTERPRETATION: These findings suggest that testosterone improves lung function in males and females, while SHBG has an opposite independent effect in males.

Pulmonary Function and Blood DNA Methylation: A Multiancestry Epigenome-Wide Association Meta-analysis.

Rationale: Methylation integrates factors present at birth and modifiable across the lifespan that can influence pulmonary function. Studies are limited in scope and replication. Objectives: To conduct large-scale epigenome-wide meta-analyses of blood DNA methylation and pulmonary function. Methods: Twelve cohorts analyzed associations of methylation at cytosine-phosphate-guanine probes (CpGs), using Illumina 450K or EPIC/850K arrays, with FEV1, FVC, and FEV1/FVC. We performed multiancestry epigenome-wide meta-analyses (total of 17,503 individuals; 14,761 European, 2,549 African, and 193 Hispanic/Latino ancestries) and interpreted results using integrative epigenomics. Measurements and Main Results: We identified 1,267 CpGs (1,042 genes) differentially methylated (false discovery rate, <0.025) in relation to FEV1, FVC, or FEV1/FVC, including 1,240 novel and 73 also related to chronic obstructive pulmonary disease (1,787 cases). We found 294 CpGs unique to European or African ancestry and 395 CpGs unique to never or ever smokers. The majority of significant CpGs correlated with nearby gene expression in blood. Findings were enriched in key regulatory elements for gene function, including accessible chromatin elements, in both blood and lung. Sixty-nine implicated genes are targets of investigational or approved drugs. One example novel gene highlighted by integrative epigenomic and druggable target analysis is TNFRSF4. Mendelian randomization and colocalization analyses suggest that epigenome-wide association study signals capture causal regulatory genomic loci. Conclusions: We identified numerous novel loci differentially methylated in relation to pulmonary function; few were detected in large genome-wide association studies. Integrative analyses highlight functional relevance and potential therapeutic targets. This comprehensive discovery of potentially modifiable, novel lung function loci expands knowledge gained from genetic studies, providing insights into lung pathogenesis.

Risks of COVID-19 by occupation in NHS workers in England.

OBJECTIVE: To quantify occupational risks of COVID-19 among healthcare staff during the first wave (9 March 2020-31 July 2020) of the pandemic in England. METHODS: We used pseudonymised data on 902 813 individuals employed by 191 National Health Service trusts to explore demographic and occupational risk factors for sickness absence ascribed to COVID-19 (n=92 880). We estimated ORs by multivariable logistic regression. RESULTS: With adjustment for employing trust, demographic characteristics and previous frequency of sickness absence, risk relative to administrative/clerical occupations was highest in 'additional clinical services' (care assistants and other occupations directly supporting those in clinical roles) (OR 2.31 (2.25 to 2.37)), registered nursing and midwifery professionals (OR 2.28 (2.23 to 2.34)) and allied health professionals (OR 1.94 (1.88 to 2.01)) and intermediate in doctors and dentists (OR 1.55 (1.50 to 1.61)). Differences in risk were higher after the employing trust had started to care for documented patients with COVID-19, and were reduced, but not eliminated, following additional adjustment for exposure to infected patients or materials, assessed by a job-exposure matrix. For prolonged COVID-19 sickness absence (episodes lasting >14 days), the variation in risk by staff group was somewhat greater. CONCLUSIONS: After allowance for possible bias and confounding by non-occupational exposures, we estimated that relative risks for COVID-19 among most patient-facing occupations were between 1.5 and 2.5. The highest risks were in those working in additional clinical services, nursing and midwifery and in allied health professions. Better protective measures for these staff groups should be a priority. COVID-19 may meet criteria for compensation as an occupational disease in some healthcare occupations. TRIAL REGISTRATION NUMBER: ISRCTN36352994.

Changing patterns of sickness absence among healthcare workers in England during the COVID-19 pandemic.

BACKGROUND: Patterns of sickness absence shed useful light on disease occurrence and illness-related behaviours in working populations. METHODS: We analysed prospectively collected, pseudonymized data on 959 356 employees who were continuously employed by National Health Service trusts in England from 1 January 2019 to 31 July 2020, comparing the frequency of new sickness absence in 2020 with that at corresponding times in 2019. RESULTS: After exclusion of episodes directly related to COVID-19, the overall incidence of sickness absence during the initial 10 weeks of the pandemic (March-May 2020) was more than 20% lower than in corresponding weeks of 2019. Trends for specific categories of illness varied substantially, with a fall by 24% for cancer, but an increase for mental illness. A doubling of new absences for pregnancy-related disorders during May-July of 2020 was limited to women with earlier COVID-19 sickness absence. CONCLUSIONS: Various factors will have contributed to the large and divergent changes that were observed. The findings reinforce concerns regarding delays in diagnosis and treatment of cancers and support a need to plan for a large backlog of treatment for many other diseases. Further research should explore the rise in absence for pregnancy-related disorders among women with earlier COVID-19 sickness absence.

Ethnic differences in risk of severe Covid-19: To what extent are they driven by exposure?

BACKGROUND: This study quantifies the risk of Covid-19 among ethnic groups of healthcare staff during the first pandemic wave in England. METHODS: We analysed data on 959 356 employees employed by 191 National Health Service trusts during 1 January 2019 to 31 July 2020, comparing rates of Covid-19 sickness absence in different ethnic groups. RESULTS: In comparison with White ethnic groups, the risk of short-duration Covid-19 sickness absence was modestly elevated in South Asian but not Black groups. However, all Black and ethnic minority groups were at higher risk of prolonged Covid-19 sickness absence. Odds ratios (ORs) relative to White ethnicity were more than doubled in South Asian groups (Indian OR 2.49, 95% confidence interval (CI) 2.36-2.63; Pakistani OR 2.38, 2.15-2.64; Bangladeshi OR 2.38, 1.98-2.86), while that for Black African ethnicity was 1.82 (1.71-1.93). In nursing/midwifery staff, the association of ethnicity with prolonged Covid-19 sickness absence was strong; the odds of South Asian nurses/midwives having a prolonged episode of Covid-19 sickness absence were increased 3-fold (OR 3.05, 2.82-3.30). CONCLUSIONS: Residual differences in risk of short term Covid-19 sickness absences among ethnic groups may reflect differences in non-occupational exposure to SARS-CoV-2. Our results indicate ethnic differences in vulnerability to Covid-19, which may be only partly explained by medical comorbidities.

Occupational exposure to gases/fumes and mineral dust affect DNA methylation levels of genes regulating expression.

Many workers are daily exposed to occupational agents like gases/fumes, mineral dust or biological dust, which could induce adverse health effects. Epigenetic mechanisms, such as DNA methylation, have been suggested to play a role. We therefore aimed to identify differentially methylated regions (DMRs) upon occupational exposures in never-smokers and investigated if these DMRs associated with gene expression levels. To determine the effects of occupational exposures independent of smoking, 903 never-smokers of the LifeLines cohort study were included. We performed three genome-wide methylation analyses (Illumina 450 K), one per occupational exposure being gases/fumes, mineral dust and biological dust, using robust linear regression adjusted for appropriate confounders. DMRs were identified using comb-p in Python. Results were validated in the Rotterdam Study (233 never-smokers) and methylation-expression associations were assessed using Biobank-based Integrative Omics Study data (n = 2802). Of the total 21 significant DMRs, 14 DMRs were associated with gases/fumes and 7 with mineral dust. Three of these DMRs were associated with both exposures (RPLP1 and LINC02169 (2×)) and 11 DMRs were located within transcript start sites of gene expression regulating genes. We replicated two DMRs with gases/fumes (VTRNA2-1 and GNAS) and one with mineral dust (CCDC144NL). In addition, nine gases/fumes DMRs and six mineral dust DMRs significantly associated with gene expression levels. Our data suggest that occupational exposures may induce differential methylation of gene expression regulating genes and thereby may induce adverse health effects. Given the millions of workers that are exposed daily to occupational exposures, further studies on this epigenetic mechanism and health outcomes are warranted.

COPD GWAS variant at 19q13.2 in relation with DNA methylation and gene expression.

Chronic obstructive pulmonary disease (COPD) is among the major health burdens in adults. While cigarette smoking is the leading risk factor, a growing number of genetic variations have been discovered to influence disease susceptibility. Epigenetic modifications may mediate the response of the genome to smoking and regulate gene expression. Chromosome 19q13.2 region is associated with both smoking and COPD, yet its functional role is unclear. Our study aimed to determine whether rs7937 (RAB4B, EGLN2), a top genetic variant in 19q13.2 region identified in genome-wide association studies of COPD, is associated with differential DNA methylation in blood (N = 1490) and gene expression in blood (N = 721) and lungs (N = 1087). We combined genetic and epigenetic data from the Rotterdam Study (RS) to perform the epigenome-wide association analysis of rs7937. Further, we used genetic and transcriptomic data from blood (RS) and from lung tissue (Lung expression quantitative trait loci mapping study), to perform the transcriptome-wide association study of rs7937. Rs7937 was significantly (FDR < 0.05) and consistently associated with differential DNA methylation in blood at 4 CpG sites in cis, independent of smoking. One methylation site (cg11298343-EGLN2) was also associated with COPD (P = 0.001). Additionally, rs7937 was associated with gene expression levels in blood in cis (EGLN2), 42% mediated through cg11298343, and in lung tissue, in cis and trans (NUMBL, EGLN2, DNMT3A, LOC101929709 and PAK2). Our results suggest that changes of DNA methylation and gene expression may be intermediate steps between genetic variants and COPD, but further causal studies in lung tissue should confirm this hypothesis.

A cross-omics integrative study of metabolic signatures of chronic obstructive pulmonary disease.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) is a common lung disorder characterized by persistent and progressive airflow limitation as well as systemic changes. Metabolic changes in blood may help detect COPD in an earlier stage and predict prognosis. METHODS: We conducted a comprehensive study of circulating metabolites, measured by proton Nuclear Magnetic Resonance Spectroscopy, in relation with COPD and lung function. The discovery sample consisted of 5557 individuals from two large population-based studies in the Netherlands, the Rotterdam Study and the Erasmus Rucphen Family study. Significant findings were replicated in 12,205 individuals from the Lifelines-DEEP study, FINRISK and the Prospective Investigation of the Vasculature in Uppsala Seniors (PIVUS) studies. For replicated metabolites further investigation of causality was performed, utilizing genetics in the Mendelian randomization approach. RESULTS: There were 602 cases of COPD and 4955 controls used in the discovery meta-analysis. Our logistic regression results showed that higher levels of plasma Glycoprotein acetyls (GlycA) are significantly associated with COPD (OR = 1.16, P = 5.6 × 10- 4 in the discovery and OR = 1.30, P = 1.8 × 10- 6 in the replication sample). A bi-directional two-sample Mendelian randomization analysis suggested that circulating blood GlycA is not causally related to COPD, but that COPD causally increases GlycA levels. Using the prospective data of the same sample of Rotterdam Study in Cox-regression, we show that the circulating GlycA level is a predictive biomarker of COPD incidence (HR = 1.99, 95%CI 1.52-2.60, comparing those in the highest and lowest quartile of GlycA) but is not significantly associated with mortality in COPD patients (HR = 1.07, 95%CI 0.94-1.20). CONCLUSIONS: Our study shows that circulating blood GlycA is a biomarker of early COPD pathology.

Age at menopause and lung function: a Mendelian randomisation study.

In observational studies, early menopause is associated with lower forced vital capacity (FVC) and a higher risk of spirometric restriction, but not airflow obstruction. It is, however, unclear if this association is causal. We therefore used a Mendelian randomisation (MR) approach, which is not affected by classical confounding, to assess the effect of age at natural menopause on lung function.We included 94 742 naturally post-menopausal women from the UK Biobank and performed MR analyses on the effect of age at menopause on forced expiratory volume in 1 s (FEV1), FVC, FEV1/FVC, spirometric restriction (FVC

Epigenome-wide association study of lung function level and its change.

Previous reports link differential DNA methylation (DNAme) to environmental exposures that are associated with lung function. Direct evidence on lung function DNAme is, however, limited. We undertook an agnostic epigenome-wide association study (EWAS) on pre-bronchodilation lung function and its change in adults.In a discovery-replication EWAS design, DNAme in blood and spirometry were measured twice, 6-15 years apart, in the same participants of three adult population-based discovery cohorts (n=2043). Associated DNAme markers (p<5×10-7) were tested in seven replication cohorts (adult: n=3327; childhood: n=420). Technical bias-adjusted residuals of a regression of the normalised absolute ß-values on control probe-derived principle components were regressed on level and change of forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC) and their ratio (FEV1/FVC) in the covariate-adjusted discovery EWAS. Inverse-variance-weighted meta-analyses were performed on results from discovery and replication samples in all participants and never-smokers.EWAS signals were enriched for smoking-related DNAme. We replicated 57 lung function DNAme markers in adult, but not childhood samples, all previously associated with smoking. Markers not previously associated with smoking failed replication. cg05575921 (AHRR (aryl hydrocarbon receptor repressor)) showed the statistically most significant association with cross-sectional lung function (FEV1/FVC: pdiscovery=3.96×10-21 and pcombined=7.22×10-50). A score combining 10 DNAme markers previously reported to mediate the effect of smoking on lung function was associated with lung function (FEV1/FVC: p=2.65×10-20).Our results reveal that lung function-associated methylation signals in adults are predominantly smoking related, and possibly of clinical utility in identifying poor lung function and accelerated decline. Larger studies with more repeat time-points are needed to identify lung function DNAme in never-smokers and in children.

DNA methylation is associated with lung function in never smokers.

BACKGROUND: Active smoking is the main risk factor for COPD. Here, epigenetic mechanisms may play a role, since cigarette smoking is associated with differential DNA methylation in whole blood. So far, it is unclear whether epigenetics also play a role in subjects with COPD who never smoked. Therefore, we aimed to identify differential DNA methylation associated with lung function in never smokers. METHODS: We determined epigenome-wide DNA methylation levels of 396,243 CpG-sites (Illumina 450 K) in blood of never smokers in four independent cohorts, LifeLines COPD&C (N = 903), LifeLines DEEP (N = 166), Rotterdam Study (RS)-III (N = 150) and RS-BIOS (N = 206). We meta-analyzed the cohort-specific methylation results to identify differentially methylated CpG-sites with FEV1/FVC. Expression Quantitative Trait Methylation (eQTM) analysis was performed in the Biobank-based Integrative Omics Studies (BIOS). RESULTS: A total of 36 CpG-sites were associated with FEV1/FVC in never smokers at p-value< 0.0001, but the meta-analysis did not reveal any epigenome-wide significant CpG-sites. Of interest, 35 of these 36 CpG-sites have not been associated with lung function before in studies including subjects irrespective of smoking history. Among the top hits were cg10012512, cg02885771, annotated to the gene LTV1 Ribosome Biogenesis factor (LTV1), and cg25105536, annotated to Kelch Like Family Member 32 (KLHL32). Moreover, a total of 11 eQTMS were identified. CONCLUSIONS: With the identification of 35 CpG-sites that are unique for never smokers, our study shows that DNA methylation is also associated with FEV1/FVC in subjects that never smoked and therefore not merely related to smoking.

Limited overlap in significant hits between genome-wide association studies on two airflow obstruction definitions in the same population.

BACKGROUND: Airflow obstruction is a hallmark of chronic obstructive pulmonary disease (COPD), and is defined as either the ratio between forced expiratory volume in one second and forced vital capacity (FEV1/FVC) < 70% or < lower limit of normal (LLN). This study aimed to assess the overlap between genome-wide association studies (GWAS) on airflow obstruction using these two definitions in the same population stratified by smoking. METHODS: GWASes were performed in the LifeLines Cohort Study for both airflow obstruction definitions in never-smokers (NS = 5071) and ever-smokers (ES = 4855). The FEV1/FVC < 70% models were adjusted for sex, age, and height; FEV1/FVC < LLN models were not adjusted. Ever-smokers models were additionally adjusted for pack-years and current-smoking. The overlap in significantly associated SNPs between the two definitions and never/ever-smokers was assessed using several p-value thresholds. To quantify the agreement, the Pearson correlation coefficient was calculated between the p-values and ORs. Replication was performed in the Vlagtwedde-Vlaardingen study (NS = 432, ES = 823). The overlapping SNPs with p < 10- 4 were validated in the Vlagtwedde-Vlaardingen and Rotterdam Study cohorts (NS = 1966, ES = 3134) and analysed for expression quantitative trait loci (eQTL) in lung tissue (n = 1087). RESULTS: In the LifeLines cohort, 96% and 93% of the never- and ever-smokers were classified concordantly based on the two definitions. 26 and 29% of the investigated SNPs were overlapping at p < 0.05 in never- and ever-smokers, respectively. At p < 10- 4 the overlap was 4% and 6% respectively, which could be change findings as shown by simulation studies. The effect estimates of the SNPs of the two definitions correlated strongly, but the p-values showed more variation and correlated only moderately. Similar observations were made in the Vlagtwedde-Vlaardingen study. Two overlapping SNPs in never-smokers (NFYC and FABP7) had the same direction of effect in the validation cohorts and the NFYC SNP was an eQTL for NFYC-AS1. NFYC is a transcription factor that binds to several known COPD genes, and FABP7 may be involved in abnormal pulmonary development. CONCLUSIONS: The definition of airflow obstruction and the population under study may be important determinants of which SNPs are associated with airflow obstruction. The genes FABP7 and NFYC(-AS1) could play a role in airflow obstruction in never-smokers specifically.

Age at puberty and risk of asthma: A Mendelian randomisation study.

BACKGROUND: Observational studies on pubertal timing and asthma, mainly performed in females, have provided conflicting results about a possible association of early puberty with higher risk of adult asthma, possibly due to residual confounding. To overcome issues of confounding, we used Mendelian randomisation (MR), i.e., genetic variants were used as instrumental variables to estimate causal effects of early puberty on post-pubertal asthma in both females and males. METHODS AND FINDINGS: MR analyses were performed in UK Biobank on 243,316 women using 254 genetic variants for age at menarche, and on 192,067 men using 46 variants for age at voice breaking. Age at menarche, recorded in years, was categorised as early (<12), normal (12-14), or late (>14); age at voice breaking was recorded and analysed as early (younger than average), normal (about average age), or late (older than average). In females, we found evidence for a causal effect of pubertal timing on asthma, with an 8% increase in asthma risk for early menarche (odds ratio [OR] 1.08; 95% CI 1.04 to 1.12; p = 8.7 × 10(-5)) and an 8% decrease for late menarche (OR 0.92; 95% CI 0.89 to 0.97; p = 3.4 × 10(-4)), suggesting a continuous protective effect of increasing age at puberty. In males, we found very similar estimates of causal effects, although with wider confidence intervals (early voice breaking: OR 1.07; 95% CI 1.00 to 1.16; p = 0.06; late voice breaking: OR 0.93; 95% CI 0.87 to 0.99; p = 0.03). We detected only modest pleiotropy, and our findings showed robustness when different methods to account for pleiotropy were applied. BMI may either introduce pleiotropy or lie on the causal pathway; secondary analyses excluding variants associated with BMI yielded similar results to those of the main analyses. Our study relies on self-reported exposures and outcomes, which may have particularly affected the power of the analyses on age at voice breaking. CONCLUSIONS: This large MR study provides evidence for a causal detrimental effect of early puberty on asthma, and does not support previous observational findings of a U-shaped relationship between pubertal timing and asthma. Common biological or psychological mechanisms associated with early puberty might explain the similarity of our results in females and males, but further research is needed to investigate this. Taken together with evidence for other detrimental effects of early puberty on health, our study emphasises the need to further investigate and address the causes of the secular shift towards earlier puberty observed worldwide.

From blood to lung tissue: effect of cigarette smoke on DNA methylation and lung function.

BACKGROUND: Genetic and environmental factors play a role in the development of COPD. The epigenome, and more specifically DNA methylation, is recognized as important link between these factors. We postulate that DNA methylation is one of the routes by which cigarette smoke influences the development of COPD. In this study, we aim to identify CpG-sites that are associated with cigarette smoke exposure and lung function levels in whole blood and validate these CpG-sites in lung tissue. METHODS: The association between pack years and DNA methylation was studied genome-wide in 658 current smokers with >5 pack years using robust linear regression analysis. Using mediation analysis, we subsequently selected the CpG-sites that were also associated with lung function levels. Significant CpG-sites were validated in lung tissue with pyrosequencing and expression quantitative trait methylation (eQTM) analysis was performed to investigate the association between DNA methylation and gene expression. RESULTS: 15 CpG-sites were significantly associated with pack years and 10 of these were additionally associated with lung function levels. We validated 5 CpG-sites in lung tissue and found several associations between DNA methylation and gene expression. CONCLUSION: This study is the first to validate a panel of CpG-sites that are associated with cigarette smoking and lung function levels in whole blood in the tissue of interest: lung tissue.

Occupational exposure to pesticides is associated with differential DNA methylation.

OBJECTIVES: Occupational pesticide exposure is associated with a wide range of diseases, including lung diseases, but it is largely unknown how pesticides influence airway disease pathogenesis. A potential mechanism might be through epigenetic mechanisms, like DNA methylation. Therefore, we assessed associations between occupational exposure to pesticides and genome-wide DNA methylation sites. METHODS: 1561 subjects of LifeLines were included with either no (n=1392), low (n=108) or high (n=61) exposure to any type of pesticides (estimated based on current or last held job). Blood DNA methylation levels were measured using Illumina 450K arrays. Associations between pesticide exposure and 420 938 methylation sites (CpGs) were assessed using robust linear regression adjusted for appropriate confounders. In addition, we performed genome-wide stratified and interaction analyses by gender, smoking and airway obstruction status, and assessed associations between gene expression and methylation for genome-wide significant CpGs (n=2802). RESULTS: In total for all analyses, high pesticide exposure was genome-wide significantly (false discovery rate P<0.05) associated with differential DNA methylation of 31 CpGs annotated to 29 genes. Twenty of these CpGs were found in subjects with airway obstruction. Several of the identified genes, for example, RYR1, ALLC, PTPRN2, LRRC3B, PAX2 and VTRNA2-1, are genes previously linked to either pesticide exposure or lung-related diseases. Seven out of 31 CpGs were associated with gene expression levels. CONCLUSIONS: We show for the first time that occupational exposure to pesticides is genome-wide associated with differential DNA methylation. Further research should reveal whether this differential methylation plays a role in the airway disease pathogenesis induced by pesticides.

Genome-wide association study on the FEV1/FVC ratio in never-smokers identifies HHIP and FAM13A.

BACKGROUND: Although a striking proportion (25% to 45%) of patients with chronic obstructive pulmonary disease are never-smokers, most genetic susceptibility studies have not focused on this group exclusively. OBJECTIVE: The aim of this study was to identify common genetic variants associated with FEV1 and its ratio to forced vital capacity (FVC) in never-smokers. METHODS: Genome-wide association studies were performed in 5070 never-smokers of the identification cohort LifeLines, and results (P < 10-5) were verified by using a meta-analysis of the Vlagtwedde-Vlaardingen study and the Rotterdam Study I-III (total n = 1966). Furthermore, we aimed to assess the effects of the replicated variants in more detail by performing genetic risk score, expression quantitative trait loci, and variant*ever-smoking interaction analyses. RESULTS: We identified associations between the FEV1/FVC ratio and 5 common genetic variants in the identification cohort, and 2 of these associations were replicated. The 2 variants annotated to the genes hedgehog interacting protein (HHIP) and family with sequence similarity 13 member A (FAM13A) were shown to have an additive effect on FEV1/FVC levels in the genetic risk score analysis; were associated with gene expression of HHIP and FAM13A in lung tissue, respectively; and were genome-wide significant in a meta-analysis including both identification and 4 verification cohorts (P < 2.19 × 10-7). Finally, we did not identify significant interactions between the variants and ever smoking. Results of the FEV1 identification analysis were not replicated. CONCLUSION: The genes HHIP and FAM13A confer a risk for airway obstruction in general that is not driven exclusively by cigarette smoking, which is the main risk factor for chronic obstructive pulmonary disease.