Air pollution, lung function and mortality: survival and mediation analyses in UK Biobank.

Background: Air pollution is associated with lower lung function, and both are associated with premature mortality and cardiovascular disease (CVD). Evidence remains scarce on the potential mediating effect of impaired lung function on the association between air pollution and mortality or CVD. Methods: We used data from UK Biobank (n∼200 000 individuals) with 8-year follow-up to mortality and incident CVD. Exposures to particulate matter <10 µm (PM10), particulate matter <2.5 µm (PM2.5) and nitrogen dioxide (NO2) were assessed by land-use regression modelling. Lung function (forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC) and the FEV1/FVC ratio) was measured between 2006 and 2010 and transformed to Global Lung Function Initiative (GLI) z-scores. Adjusted Cox proportional hazards and causal proportional hazards mediation analysis models were fitted, stratified by smoking status. Results: Lower FEV1 and FVC were associated with all-cause and CVD mortality, and incident CVD, with larger estimates in ever- than never-smokers (all-cause mortality hazard ratio per FEV1 GLI z-score decrease 1.29 (95% CI 1.24-1.34) for ever-smokers and 1.16 (95% CI 1.12-1.21) for never-smokers). Long-term exposure to PM2.5 or NO2 was associated with incident CVD, with similar effect sizes for ever- and never-smokers. Mediated proportions of the air pollution-all-cause mortality estimates driven by FEV1 were 18% (95% CI 2-33%) for PM2.5 and 27% (95% CI 3-51%) for NO2. Corresponding mediated proportions for incident CVD were 9% (95% CI 4-13%) for PM2.5 and 16% (95% CI 6-25%) for NO2. Conclusions: Lung function may mediate a modest proportion of associations between air pollution and mortality and CVD outcomes. Results likely reflect the extent of either shared mechanisms or direct effects relating to lower lung function caused by air pollution.

Genome-wide association study of thyroid-stimulating hormone highlights new genes, pathways and associations with thyroid disease.

Thyroid hormones play a critical role in regulation of multiple physiological functions and thyroid dysfunction is associated with substantial morbidity. Here, we use electronic health records to undertake a genome-wide association study of thyroid-stimulating hormone (TSH) levels, with a total sample size of 247,107. We identify 158 novel genetic associations, more than doubling the number of known associations with TSH, and implicate 112 putative causal genes, of which 76 are not previously implicated. A polygenic score for TSH is associated with TSH levels in African, South Asian, East Asian, Middle Eastern and admixed American ancestries, and associated with hypothyroidism and other thyroid disease in South Asians. In Europeans, the TSH polygenic score is associated with thyroid disease, including thyroid cancer and age-of-onset of hypothyroidism and hyperthyroidism. We develop pathway-specific genetic risk scores for TSH levels and use these in phenome-wide association studies to identify potential consequences of pathway perturbation. Together, these findings demonstrate the potential utility of genetic associations to inform future therapeutics and risk prediction for thyroid diseases.

Identification and characterisation of a rare MTTP variant underlying hereditary non-alcoholic fatty liver disease.

Background & Aims: Non-alcoholic fatty liver disease (NAFLD) is a complex trait with an estimated prevalence of 25% globally. We aimed to identify the genetic variant underlying a four-generation family with progressive NAFLD leading to cirrhosis, decompensation, and development of hepatocellular carcinoma in the absence of common risk factors such as obesity and type 2 diabetes. Methods: Exome sequencing and genome comparisons were used to identify the likely causal variant. We extensively characterised the clinical phenotype and post-prandial metabolic responses of family members with the identified novel variant in comparison with healthy non-carriers and wild-type patients with NAFLD. Variant-expressing hepatocyte-like cells (HLCs) were derived from human-induced pluripotent stem cells generated from homozygous donor skin fibroblasts and restored to wild-type using CRISPR-Cas9. The phenotype was assessed using imaging, targeted RNA analysis, and molecular expression arrays. Results: We identified a rare causal variant c.1691T>C p.I564T (rs745447480) in MTTP, encoding microsomal triglyceride transfer protein (MTP), associated with progressive NAFLD, unrelated to metabolic syndrome and without characteristic features of abetalipoproteinaemia. HLCs derived from a homozygote donor had significantly lower MTP activity and lower lipoprotein ApoB secretion than wild-type cells, while having similar levels of MTP mRNA and protein. Cytoplasmic triglyceride accumulation in HLCs triggered endoplasmic reticulum stress, secretion of pro-inflammatory mediators, and production of reactive oxygen species. Conclusions: We have identified and characterised a rare causal variant in MTTP, and homozygosity for MTTP p.I564T is associated with progressive NAFLD without any other manifestations of abetalipoproteinaemia. Our findings provide insights into mechanisms driving progressive NAFLD. Impact and Implications: A rare genetic variant in the gene MTTP has been identified as responsible for the development of severe non-alcoholic fatty liver disease in a four-generation family with no typical disease risk factors. A cell line culture created harbouring this variant gene was characterised to understand how this genetic variation leads to a defect in liver cells, which results in accumulation of fat and processes that promote disease. This is now a useful model for studying the disease pathways and to discover new ways to treat common types of fatty liver disease.

A Polygenic Risk Score for Idiopathic Pulmonary Fibrosis and Interstitial Lung Abnormalities.

Rationale: In addition to rare genetic variants and the MUC5B locus, common genetic variants contribute to idiopathic pulmonary fibrosis (IPF) risk. The predictive power of common variants outside the MUC5B locus for IPF and interstitial lung abnormalities (ILAs) is unknown. Objectives: We tested the predictive value of IPF polygenic risk scores (PRSs) with and without the MUC5B region on IPF, ILA, and ILA progression. Methods: We developed PRSs that included (PRS-M5B) and excluded (PRS-NO-M5B) the MUC5B region (500-kb window around rs35705950-T) using an IPF genome-wide association study. We assessed PRS associations with area under the receiver operating characteristic curve (AUC) metrics for IPF, ILA, and ILA progression. Measurements and Main Results: We included 14,650 participants (1,970 IPF; 1,068 ILA) from six multi-ancestry population-based and case-control cohorts. In cases excluded from genome-wide association study, the PRS-M5B (odds ratio [OR] per SD of the score, 3.1; P = 7.1 × 10-95) and PRS-NO-M5B (OR per SD, 2.8; P = 2.5 × 10-87) were associated with IPF. Participants in the top PRS-NO-M5B quintile had ∼sevenfold odds for IPF compared with those in the first quintile. A clinical model predicted IPF (AUC, 0.61); rs35705950-T and PRS-NO-M5B demonstrated higher AUCs (0.73 and 0.7, respectively), and adding both genetic predictors to a clinical model yielded the highest performance (AUC, 0.81). The PRS-NO-M5B was associated with ILA (OR, 1.25) and ILA progression (OR, 1.16) in European ancestry participants. Conclusions: A common genetic variant risk score complements the MUC5B variant to identify individuals at high risk of interstitial lung abnormalities and pulmonary fibrosis.

PCSK6 and Survival in Idiopathic Pulmonary Fibrosis.

Rationale: Idiopathic pulmonary fibrosis (IPF) is a devastating disease characterized by limited treatment options and high mortality. A better understanding of the molecular drivers of IPF progression is needed. Objectives: To identify and validate molecular determinants of IPF survival. Methods: A staged genome-wide association study was performed using paired genomic and survival data. Stage I cases were drawn from centers across the United States and Europe and stage II cases from Vanderbilt University. Cox proportional hazards regression was used to identify gene variants associated with differential transplantation-free survival (TFS). Stage I variants with nominal significance (P < 5 × 10-5) were advanced for stage II testing and meta-analyzed to identify those reaching genome-wide significance (P < 5 × 10-8). Downstream analyses were performed for genes and proteins associated with variants reaching genome-wide significance. Measurements and Main Results: After quality controls, 1,481 stage I cases and 397 stage II cases were included in the analysis. After filtering, 9,075,629 variants were tested in stage I, with 158 meeting advancement criteria. Four variants associated with TFS with consistent effect direction were identified in stage II, including one in an intron of PCSK6 (proprotein convertase subtilisin/kexin type 6) reaching genome-wide significance (hazard ratio, 4.11 [95% confidence interval, 2.54-6.67]; P = 9.45 × 10-9). PCSK6 protein was highly expressed in IPF lung parenchyma. PCSK6 lung staining intensity, peripheral blood gene expression, and plasma concentration were associated with reduced TFS. Conclusions: We identified four novel variants associated with IPF survival, including one in PCSK6 that reached genome-wide significance. Downstream analyses suggested that PCSK6 protein plays a potentially important role in IPF progression.

Longitudinal lung function and gas transfer in individuals with idiopathic pulmonary fibrosis: a genome-wide association study.

BACKGROUND: Idiopathic pulmonary fibrosis (IPF) is an incurable lung disease characterised by progressive scarring leading to alveolar stiffness, reduced lung capacity, and impeded gas transfer. We aimed to identify genetic variants associated with declining lung capacity or declining gas transfer after diagnosis of IPF. METHODS: We did a genome-wide meta-analysis of longitudinal measures of forced vital capacity (FVC) and diffusing capacity of the lung for carbon monoxide (DLCO) in individuals diagnosed with IPF. Individuals were recruited to three studies between June, 1996, and August, 2017, from across centres in the US, UK, and Spain. Suggestively significant variants were investigated further in an additional independent study (CleanUP-IPF). All four studies diagnosed cases following American Thoracic Society/European Respiratory Society guidelines. Variants were defined as significantly associated if they had a meta-analysis p<5 × 10-8 when meta-analysing across all discovery and follow-up studies, had consistent direction of effects across all four studies, and were nominally significant (p<0·05) in each study. FINDINGS: 1329 individuals with a total of 5216 measures were included in the FVC analysis. 975 individuals with a total of 3361 measures were included in the DLCO analysis. For the discovery genome-wide analyses, 7 611 174 genetic variants were included in the FVC analysis and 7 536 843 in the DLCO analysis. One variant (rs115982800) located in an antisense RNA gene for protein kinase N2 (PKN2) showed a genome-wide significant association with FVC decline (-140 mL/year per risk allele [95% CI -180 to -100]; p=9·14 × 10-12). INTERPRETATION: Our analysis identifies a genetic variant associated with disease progression, which might highlight a new biological mechanism for IPF. We found that PKN2, a Rho and Rac effector protein, is the most likely gene of interest from this analysis. PKN2 inhibitors are currently in development and signify a potential novel therapeutic approach for IPF. FUNDING: Action for Pulmonary Fibrosis, Medical Research Council, Wellcome Trust, and National Institutes of Health National Heart, Lung, and Blood Institute.

Genetic insights into smoking behaviours in 10,558 men of African ancestry from continental Africa and the UK.

Smoking is a leading risk factor for many of the top ten causes of death worldwide. Of the 1.3 billion smokers globally, 80% live in low- and middle-income countries, where the number of deaths due to tobacco use is expected to double in the next decade according to the World Health Organization. Genetic studies have helped to identify biological pathways for smoking behaviours, but have mostly focussed on individuals of European ancestry or living in either North America or Europe. We performed a genome-wide association study of two smoking behaviour traits in 10,558 men of African ancestry living in five African countries and the UK. Eight independent variants were associated with either smoking initiation or cessation at P-value < 5 × 10-6, four being monomorphic or rare in European populations. Gene prioritisation strategy highlighted five genes, including SEMA6D, previously described as associated with several smoking behaviour traits. These results confirm the importance of analysing underrepresented populations in genetic epidemiology, and the urgent need for larger genomic studies to boost discovery power to better understand smoking behaviours, as well as many other traits.

Risk factors associated with SARS-CoV-2 infection in a multiethnic cohort of United Kingdom healthcare workers (UK-REACH): A cross-sectional analysis.

BACKGROUND: Healthcare workers (HCWs), particularly those from ethnic minority groups, have been shown to be at disproportionately higher risk of infection with Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) compared to the general population. However, there is insufficient evidence on how demographic and occupational factors influence infection risk among ethnic minority HCWs. METHODS AND FINDINGS: We conducted a cross-sectional analysis using data from the baseline questionnaire of the United Kingdom Research study into Ethnicity and Coronavirus Disease 2019 (COVID-19) Outcomes in Healthcare workers (UK-REACH) cohort study, administered between December 2020 and March 2021. We used logistic regression to examine associations of demographic, household, and occupational risk factors with SARS-CoV-2 infection (defined by polymerase chain reaction (PCR), serology, or suspected COVID-19) in a diverse group of HCWs. The primary exposure of interest was self-reported ethnicity. Among 10,772 HCWs who worked during the first UK national lockdown in March 2020, the median age was 45 (interquartile range [IQR] 35 to 54), 75.1% were female and 29.6% were from ethnic minority groups. A total of 2,496 (23.2%) reported previous SARS-CoV-2 infection. The fully adjusted model contained the following dependent variables: demographic factors (age, sex, ethnicity, migration status, deprivation, religiosity), household factors (living with key workers, shared spaces in accommodation, number of people in household), health factors (presence/absence of diabetes or immunosuppression, smoking history, shielding status, SARS-CoV-2 vaccination status), the extent of social mixing outside of the household, and occupational factors (job role, the area in which a participant worked, use of public transport to work, exposure to confirmed suspected COVID-19 patients, personal protective equipment [PPE] access, aerosol generating procedure exposure, night shift pattern, and the UK region of workplace). After adjustment, demographic and household factors associated with increased odds of infection included younger age, living with other key workers, and higher religiosity. Important occupational risk factors associated with increased odds of infection included attending to a higher number of COVID-19 positive patients (aOR 2.59, 95% CI 2.11 to 3.18 for ≥21 patients per week versus none), working in a nursing or midwifery role (1.30, 1.11 to 1.53, compared to doctors), reporting a lack of access to PPE (1.29, 1.17 to 1.43), and working in an ambulance (2.00, 1.56 to 2.58) or hospital inpatient setting (1.55, 1.38 to 1.75). Those who worked in intensive care units were less likely to have been infected (0.76, 0.64 to 0.92) than those who did not. Black HCWs were more likely to have been infected than their White colleagues, an effect which attenuated after adjustment for other known risk factors. This study is limited by self-selection bias and the cross sectional nature of the study means we cannot infer the direction of causality. CONCLUSIONS: We identified key sociodemographic and occupational risk factors associated with SARS-CoV-2 infection among UK HCWs, and have determined factors that might contribute to a disproportionate odds of infection in HCWs from Black ethnic groups. These findings demonstrate the importance of social and occupational factors in driving ethnic disparities in COVID-19 outcomes, and should inform policies, including targeted vaccination strategies and risk assessments aimed at protecting HCWs in future waves of the COVID-19 pandemic. TRIAL REGISTRATION: The study was prospectively registered at ISRCTN (reference number: ISRCTN11811602).

Genome-wide gene-air pollution interaction analysis of lung function in 300,000 individuals.

BACKGROUND: Impaired lung function is predictive of mortality and is a key component of chronic obstructive pulmonary disease. Lung function has a strong genetic component but is also affected by environmental factors such as increased exposure to air pollution, but the effect of their interactions is not well understood. OBJECTIVES: To identify interactions between genetic variants and air pollution measures which affect COPD risk and lung function. Additionally, to determine whether previously identified lung function genetic association signals showed evidence of interaction with air pollution, considering both individual effects and combined effects using a genetic risk score (GRS). METHODS: We conducted a genome-wide gene-air pollution interaction analysis of spirometry measures with three measures of air pollution at home address: particulate matter (PM2.5 & PM10) and nitrogen dioxide (NO2), in approximately 300,000 unrelated European individuals from UK Biobank. We explored air pollution interactions with previously identified lung function signals and determined their combined interaction effect using a GRS. RESULTS: We identified seven new genome-wide interaction signals (P<5×10-8), and a further ten suggestive interaction signals (P<5×10-7). Additionally, we found statistical evidence of interaction for FEV1/FVC between PM2.5 and previously identified lung function signal, rs10841302, near AEBP2, suggesting increased susceptibility as copies of the G allele increased (but size of the impact was small - interaction beta: -0.363 percentage points, 95% CI: -0.523, -0.203 per 5 µg/m3). There was no observed interaction between air pollutants and the weighted GRS. DISCUSSION: We carried out the largest genome-wide gene-air pollution interaction study of lung function and identified potential effects of clinically relevant size and significance. We observed up to 440 ml lower lung function for certain genotypes when exposed to mean levels of outdoor air pollution, which is approximately equivalent to nine years of average normal loss of lung function in adults.

Interaction of Cigarette Smoking and Polygenic Risk Score on Reduced Lung Function.

Importance: The risk of airflow limitation and chronic obstructive pulmonary disease (COPD) is influenced by combinations of cigarette smoking and genetic susceptibility, yet it remains unclear whether gene-by-smoking interactions are associated with quantitative measures of lung function. Objective: To assess the interaction of cigarette smoking and polygenic risk score in association with reduced lung function. Design, Setting, and Participants: This UK Biobank cohort study included UK citizens of European ancestry aged 40 to 69 years with genetic and spirometry data passing quality control metrics. Data was analyzed from July 2020 to March 2021. Exposures: PRS of combined forced expiratory volume in 1 second (FEV1) and percent of forced vital capacity exhaled in the first second (FEV1/FVC), self-reported pack-years of smoking, ever- vs never-smoking status, and current- vs former- or never-smoking status. Main Outcomes and Measures: FEV1/FVC was the primary outcome. Models were used to test for interactions with models, including the main effects of PRS, different smoking variables, and their cross-product terms. The association between pack-years of smoking and FEV1/FVC were compared for those in the highest vs lowest decile of estimated genetic risk for low lung function. Results: We included 319 730 individuals, of whom 24 915 (8%) had moderate-to-severe COPD cases, and 44.4% were men. Participants had a mean (SD) age 56.5 of (8.02) years. The PRS and pack-years were significantly associated with lower FEV1/FVC (PRS: ß, -0.03; 95% CI, -0.031 to -0.03; pack-years: ß, -0.0064; 95% CI, -0.0064 to -0.0063) and the interaction term (ß, -0.0028; 95% CI, -0.0029 to -0.0026). A stepwise increment in estimated effect sizes for these interaction terms was observed per 10 pack-years of smoking exposure. The interaction of PRS with 11 to 20, 31 to 40, and more than 50 pack-years categories were ß (interaction) -0.0038 (95% CI, -0.0046 to -0.0031); -0.013 (95% CI, -0.014 to -0.012); and -0.017 (95% CI, -0.019 to -0.016), respectively. There was evidence of significant interaction between PRS with ever- or never- smoking status (ß, interaction; -0.0064; 95% CI, -0.0068 to -0.0060) and current or not-current smoking (ß, interaction; -0.0091; 95% CI, -0.0097 to -0.0084). For any given level of pack-years of smoking exposure, FEV1/FVC was significantly lower for individuals in the tenth decile (ie, highest risk) than the first decile (ie, lowest risk) of genetic risk. For every 20 pack-years of smoking, those in the tenth decile compared with the first decile of genetic risk showed nearly a 2-fold reduction in FEV1/FVC. Conclusions and Relevance: COPD is characterized by diminished lung function, and our analyses suggest there is substantial interaction between genome-wide PRS and smoking exposures. While smoking was associated with decreased lung function across all genetic risk categories, the associations were strongest in individuals with higher estimated genetic risk.

Pleiotropic associations of heterozygosity for the SERPINA1 Z allele in the UK Biobank.

Homozygosity for the SERPINA1 Z allele causes α1-antitrypsin deficiency, a rare condition that can cause lung and liver disease. However, the effects of Z allele heterozygosity on nonrespiratory phenotypes, and on lung function in the general population, remain unclear. We conducted a large, population-based study to determine Z allele effects on >2400 phenotypes in the UK Biobank (N=303 353). Z allele heterozygosity was strongly associated with increased height (ß=1.02 cm, p=3.91×10-68), and with other nonrespiratory phenotypes including increased risk of gall bladder disease, reduced risk of heart disease and lower blood pressure, reduced risk of osteoarthritis and reduced bone mineral density, increased risk of headache and enlarged prostate, as well as with blood biomarkers of liver function. Heterozygosity was associated with higher height-adjusted forced expiratory volume in 1 s (FEV1) (ß=19.36 mL, p=9.21×10-4) and FEV1/forced vital capacity (ß=0.0031, p=1.22×10-5) in nonsmokers, whereas in smokers, this protective effect was abolished. Furthermore, we show for the first time that sex modifies the association of the Z allele on lung function. We conclude that Z allele heterozygosity and homozygosity exhibit opposing effects on lung function in the UK population, and that these associations are modified by smoking and sex. In exploratory analyses, heterozygosity for the Z allele also showed pleiotropic associations with nonrespiratory health-related traits and disease risk.

Genome-Wide Gene-by-Smoking Interaction Study of Chronic Obstructive Pulmonary Disease.

Risk of chronic obstructive pulmonary disease (COPD) is determined by both cigarette smoking and genetic susceptibility, but little is known about gene-by-smoking interactions. We performed a genome-wide association analysis of 179,689 controls and 21,077 COPD cases from UK Biobank subjects of European ancestry recruited from 2006 to 2010, considering genetic main effects and gene-by-smoking interaction effects simultaneously (2-degrees-of-freedom (df) test) as well as interaction effects alone (1-df interaction test). We sought to replicate significant results in COPDGene (United States, 2008-2010) and SpiroMeta Consortium (multiple countries, 1947-2015) data. We considered 2 smoking variables: 1) ever/never and 2) current/noncurrent. In the 1-df test, we identified 1 genome-wide significant locus on 15q25.1 (cholinergic receptor nicotinic ß4 subunit, or CHRNB4) for ever- and current smoking and identified PI*Z allele (rs28929474) of serpin family A member 1 (SERPINA1) for ever-smoking and 3q26.2 (MDS1 and EVI1 complex locus, or MECOM) for current smoking in an analysis of previously reported COPD loci. In the 2-df test, most of the significant signals were also significant for genetic marginal effects, aside from 16q22.1 (sphingomyelin phosphodiesterase 3, or SMPD3) and 19q13.2 (Egl-9 family hypoxia inducible factor 2, or EGLN2). The significant effects at 15q25.1 and 19q13.2 loci, both previously described in prior genome-wide association studies of COPD or smoking, were replicated in COPDGene and SpiroMeta. We identified interaction effects at previously reported COPD loci; however, we failed to identify novel susceptibility loci.

Relative contributions of family history and a polygenic risk score on COPD and related outcomes: COPDGene and ECLIPSE studies.

INTRODUCTION: Family history is a risk factor for chronic obstructive pulmonary disease (COPD). We previously developed a COPD risk score from genome-wide genetic markers (Polygenic Risk Score, PRS). Whether the PRS and family history provide complementary or redundant information for predicting COPD and related outcomes is unknown. METHODS: We assessed the predictive capacity of family history and PRS on COPD and COPD-related outcomes in non-Hispanic white (NHW) and African American (AA) subjects from COPDGene and ECLIPSE studies. We also performed interaction and mediation analyses. RESULTS: In COPDGene, family history and PRS were significantly associated with COPD in a single model (PFamHx <0.0001; PPRS<0.0001). Similar trends were seen in ECLIPSE. The area under the receiver operator characteristic curve for a model containing family history and PRS was significantly higher than a model with PRS (p=0.00035) in NHWs and a model with family history (p<0.0001) alone in NHWs and AAs. Both family history and PRS were significantly associated with measures of quantitative emphysema and airway thickness. There was a weakly positive interaction between family history and the PRS under the additive, but not multiplicative scale in NHWs (relative excess risk due to interaction=0.48, p=0.04). Mediation analyses found that a significant proportion of the effect of family history on COPD was mediated through PRS in NHWs (16.5%, 95% CI 9.4% to 24.3%), but not AAs. CONCLUSION: Family history and the PRS provide complementary information for predicting COPD and related outcomes. Future studies can address the impact of obtaining both measures in clinical practice.

The sputum microbiome is distinct between COPD and health, independent of smoking history.

BACKGROUND: Airway bacterial dysbiosis is a feature of chronic obstructive pulmonary disease (COPD). However, there is limited comparative data of the lung microbiome between healthy smokers, non-smokers and COPD. METHODS: We compared the 16S rRNA gene-based sputum microbiome generated from pair-ended Illumina sequencing of 124 healthy subjects (28 smokers and 96 non-smokers with normal lung function), with single stable samples from 218 COPD subjects collected from three UK clinical centres as part of the COPDMAP consortium. RESULTS: In healthy subjects Firmicutes, Bacteroidetes and Actinobacteria were the major phyla constituting 88% of the total reads, and Streptococcus, Veillonella, Prevotella, Actinomyces and Rothia were the dominant genera. Haemophilus formed only 3% of the healthy microbiome. In contrast, Proteobacteria was the most dominant phylum accounting for 50% of the microbiome in COPD subjects, with Haemophilus and Moraxella at genus level contributing 25 and 3% respectively. There were no differences in the microbiome profile within healthy and COPD subgroups when stratified based on smoking history. Principal coordinate analysis on operational taxonomic units showed two distinct clusters, representative of healthy and COPD subjects (PERMANOVA, p = 0·001). CONCLUSION: The healthy and COPD sputum microbiomes are distinct and independent of smoking history. Our results underline the important role for Gammaproteobacteria in COPD.

Chronic obstructive pulmonary disease and related phenotypes: polygenic risk scores in population-based and case-control cohorts.

BACKGROUND: Genetic factors influence chronic obstructive pulmonary disease (COPD) risk, but the individual variants that have been identified have small effects. We hypothesised that a polygenic risk score using additional variants would predict COPD and associated phenotypes. METHODS: We constructed a polygenic risk score using a genome-wide association study of lung function (FEV1 and FEV1/forced vital capacity [FVC]) from the UK Biobank and SpiroMeta. We tested this polygenic risk score in nine cohorts of multiple ethnicities for an association with moderate-to-severe COPD (defined as FEV1/FVC <0·7 and FEV1 <80% of predicted). Associations were tested using logistic regression models, adjusting for age, sex, height, smoking pack-years, and principal components of genetic ancestry. We assessed predictive performance of models by area under the curve. In a subset of studies, we also studied quantitative and qualitative CT imaging phenotypes that reflect parenchymal and airway pathology, and patterns of reduced lung growth. FINDINGS: The polygenic risk score was associated with COPD in European (odds ratio [OR] per SD 1·81 [95% CI 1·74-1·88] and non-European (1·42 [1·34-1·51]) populations. Compared with the first decile, the tenth decile of the polygenic risk score was associated with COPD, with an OR of 7·99 (6·56-9·72) in European ancestry and 4·83 (3·45-6·77) in non-European ancestry cohorts. The polygenic risk score was superior to previously described genetic risk scores and, when combined with clinical risk factors (ie, age, sex, and smoking pack-years), showed improved prediction for COPD compared with a model comprising clinical risk factors alone (AUC 0·80 [0·79-0·81] vs 0·76 [0·75-0·76]). The polygenic risk score was associated with CT imaging phenotypes, including wall area percent, quantitative and qualitative measures of emphysema, local histogram emphysema patterns, and destructive emphysema subtypes. The polygenic risk score was associated with a reduced lung growth pattern. INTERPRETATION: A risk score comprised of genetic variants can identify a small subset of individuals at markedly increased risk for moderate-to-severe COPD, emphysema subtypes associated with cigarette smoking, and patterns of reduced lung growth. FUNDING: US National Institutes of Health, Wellcome Trust.

Smoking, DNA Methylation, and Lung Function: a Mendelian Randomization Analysis to Investigate Causal Pathways.

Whether smoking-associated DNA methylation has a causal effect on lung function has not been thoroughly evaluated. We first investigated the causal effects of 474 smoking-associated CpGs on forced expiratory volume in 1 s (FEV1) in UK Biobank (n = 321,047) by using two-sample Mendelian randomization (MR) and then replicated this investigation in the SpiroMeta Consortium (n = 79,055). Second, we used two-step MR to investigate whether DNA methylation mediates the effect of smoking on FEV1. Lastly, we evaluated the presence of horizontal pleiotropy and assessed whether there is any evidence for shared causal genetic variants between lung function, DNA methylation, and gene expression by using a multiple-trait colocalization ("moloc") framework. We found evidence of a possible causal effect for DNA methylation on FEV1 at 18 CpGs (p < 1.2 × 10-4). Replication analysis supported a causal effect at three CpGs (cg21201401 [LIME1 and ZGPAT], cg19758448 [PGAP3], and cg12616487 [EML3 and AHNAK] [p < 0.0028]). DNA methylation did not clearly mediate the effect of smoking on FEV1, although DNA methylation at some sites might influence lung function via effects on smoking. By using "moloc", we found evidence of shared causal variants between lung function, gene expression, and DNA methylation. These findings highlight potential therapeutic targets for improving lung function and possibly smoking cessation, although larger, tissue-specific datasets are required to confirm these results.

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.

High-Throughput Sequencing in Respiratory, Critical Care, and Sleep Medicine Research. An Official American Thoracic Society Workshop Report.

High-throughput, "next-generation" sequencing methods are now being broadly applied across all fields of biomedical research, including respiratory disease, critical care, and sleep medicine. Although there are numerous review articles and best practice guidelines related to sequencing methods and data analysis, there are fewer resources summarizing issues related to study design and interpretation, especially as applied to common, complex, nonmalignant diseases. To address these gaps, a single-day workshop was held at the American Thoracic Society meeting in May 2017, led by the American Thoracic Society Section on Genetics and Genomics. The aim of this workshop was to review the design, analysis, interpretation, and functional follow-up of high-throughput sequencing studies in respiratory, critical care, and sleep medicine research. This workshop brought together experts in multiple fields, including genetic epidemiology, biobanking, bioinformatics, and research ethics, along with physician-scientists with expertise in a range of relevant diseases. The workshop focused on application of DNA and RNA sequencing research in common chronic diseases and did not cover sequencing studies in lung cancer, monogenic diseases (e.g., cystic fibrosis), or microbiome sequencing. Participants reviewed and discussed study design, data analysis and presentation, interpretation, functional follow-up, and reporting of results. This report summarizes the main conclusions of the workshop, specifically addressing the application of these methods in respiratory, critical care, and sleep medicine research. This workshop report may serve as a resource for our research community as well as for journal editors and reviewers of sequencing-based manuscript submissions in our research field.