PDF file - 242K, Table S2. Clinical characteristics of patients in the first replication set (UBC). Table S3. Clinical characteristics of patients in the second replication set (Groningen). Figure S1. Volcano plots showing the impact of smoking, COPD, and lung cancer on gene expression in the lung. Figure S2. Comparison of gene expression recovery following smoking cessation between the discovery set (Laval) and UBC. Figure S3. Comparison of gene expression recovery following smoking cessation between the discovery set (Laval) and Groningen. Figure S4. Expression of SERPIND1 in lung parenchyma of smoker and never-smoker by immunochemistry. Figure S5. Enrichment plot for slowly reversible gene in UBC. Figure S6. Enrichment plot for slowly reversible gene in Groningen

Rationale: Emerging data demonstrate that the smallest conducting airways, terminal bronchioles, are the early site of tissue destruction in chronic obstructive pulmonary disease (COPD) and are reduced by as much as 41% by the time someone is diagnosed with mild (Global Initiative for Chronic Obstructive Lung Disease [GOLD] stage 1) COPD. Objectives: To develop a single-cell atlas that describes the structural, cellular, and extracellular matrix alterations underlying terminal bronchiole loss in COPD. Methods: This cross-sectional study of 262 lung samples derived from 34 ex-smokers with normal lung function (n = 10) or GOLD stage 1 (n = 10), stage 2 (n = 8), or stage 4 (n = 6) COPD was performed to assess the morphology, extracellular matrix, single-cell atlas, and genes associated with terminal bronchiole reduction using stereology, micro-computed tomography, nonlinear optical microscopy, imaging mass spectrometry, and transcriptomics. Measurements and Main Results: The lumen area of terminal bronchioles progressively narrows with COPD severity as a result of the loss of elastin fibers within alveolar attachments, which was observed before microscopic emphysematous tissue destruction in GOLD stage 1 and 2 COPD. The single-cell atlas of terminal bronchioles in COPD demonstrated M1-like macrophages and neutrophils located within alveolar attachments and associated with the pathobiology of elastin fiber loss, whereas adaptive immune cells (naive, CD4, and CD8 T cells, and B cells) are associated with terminal bronchiole wall remodeling. Terminal bronchiole pathology was associated with the upregulation of genes involved in innate and adaptive immune responses, the interferon response, and the degranulation of neutrophils. Conclusions: This comprehensive single-cell atlas highlights terminal bronchiole alveolar attachments as the initial site of tissue destruction in centrilobular emphysema and an attractive target for disease modification.

<p>PDF file - 589K, Probe sets significantly associated with smoking in the discovery set</p>

Lung resistance ( R L ) and elastance ( E L ) can be measured during positive or negative pressure ventilation. Whether the different modes of ventilation produce different R L and E L is still being debated. Although negative pressure ventilation (NPV) is more physiological, positive pressure ventilation (PPV) is more commonly used for treating respiratory failure. In the present study, we measured lung volume, airway diameter, and airway volume, as well as R L and E L with PPV and NPV in explanted sheep lungs. We found that lung volume under a static pressure, either positive or negative, was not different. However, R L and E L were significantly higher in NPV at high inflation pressures. Interestingly, diameters of smaller airways (diameters <3.5 mm) and total airway volume were significantly greater at high negative inflation pressures compared with those at high positive inflation pressures. This suggests that NPV is more effective in distending the peripheral airways, likely due to the fact that negative pressure is applied through the pleural membrane and reaches the central airways via the peripheral airways, whereas positive pressure is applied in the opposite direction. More distension of lung periphery could explain why R L is higher in NPV (vs. PPV), because the peripheral parenchyma is a major source of tissue resistance, which is a part of the R L that increases with pressure. This explanation is consistent with the finding that during high frequency ventilation (>1 Hz, where R L reflects airway resistance more than tissue resistance), the difference in R L between NPV and PPV disappeared.

Lung resistance ( RL) is determined by airway and parenchymal tissue resistance, as well as the degree of heterogeneity in airway constriction. Deep inspirations (DIs) are known to reverse experimentally induced increase in RL, but the mechanism is not entirely clear. The first step toward understanding the effect of DI is to determine how each of the resistance components is affected by DI. In the present study, we measured RL and apparent airway resistance ( RAW, which combines the effects of airway resistance and airway heterogeneity) simultaneously before and after a DI in acetylcholine (ACh)-challenged ex vivo sheep lungs. We found that at normal breathing frequency (0.25 Hz) ACh-challenge led to a doubling of RL, 80.3% of that increase was caused by an increase in RAW; the increase in apparent tissue resistance ( RT) was insignificant. 57.7% of the increase in RAW was abolished by a single DI. After subtracting RAW from RL, the remaining RT was mostly independent of ACh-challenge and its reduction after a DI came mostly from the change in the mechanical properties of lung parenchyma. We conclude that at normal breathing frequency, RL in an unchallenged lung is mostly composed of RT, and the increase in RL due to ACh-challenge stems mostly from the increase in RAW and that both RAW and RT can be greatly reduced by a DI, likely due to a reduction in true airway resistance and heterogeneity, as well as parenchymal tissue hysteresis post DI.

Rationale: To improve disease outcomes in idiopathic pulmonary fibrosis ( IPF), it is essential to understand its early pathophysiology so that it can be targeted therapeutically. Objectives: Perform three-dimensional assessment of the IPF lung microstructure using stereology and multiresolution computed tomography (CT) imaging. Methods: Explanted lungs from patients with IPF (n= 8) and donor control subjects (n= 8) were inflated with air and frozen. CT scans were used to assess large airways. Unbiased, systematic uniform random samples (n = 8/lung) were scanned with microCT for stereological assessment of small airways (count number, and measure airway wall and lumen area) and parenchymal fibrosis (volume fraction of tissue, alveolar surface area, and septal wall thickness). Measurements and Main Results: The total number of airways on clinical CT was greater in IPF lungs than control lungs (P < 0.01), owing to an increase in the wall (P < 0.05) and lumen area (P < 0.05) resulting in more visible airways with a lumen larger than 2 mm. In IPF tissue samples without microscopic fibrosis, assessed by the volume fraction of tissue using microCT, there was a reduction in the number of the terminal (P < 0.01) and transitional (P < 0.001) bronchioles, and an increase in terminal bronchiole wall area (P < 0.001) compared with control lungs. In IPF tissue samples with microscopic parenchymal fibrosis, terminal bronchioles had increased airway wall thickness (P < 0.05) and dilated airway lumens (P < 0.001) leading to honeycomb cyst formations. Conclusions: This study has important implications for the current thinking on how the lung tissue is remodeled in IPF and highlights small airways as a potential target to modify IPF outcomes.

Deep inspiration (DI)-induced bronchodilation is the first line of defense against bronchoconstriction in healthy subjects. A hallmark of asthma is the lack of this beneficial effect of DI. The mechanism underlying the bronchodilatory effect of DI is not clear. Understanding the mechanism will help us unravel the mystery of asthma pathophysiology. It has been postulated that straining airway smooth muscle (ASM) during a DI could lead to bronchodilation and bronchoprotection. The hypothesis is currently under debate, and a central question is whether ASM is sufficiently stretched during a DI for its contractility to be compromised. Besides bronchoconstriction, another contributor to lung resistance is airway heterogeneity. The present study examines changes in airway diameter and heterogeneity at different lung volumes. Freshly explanted sheep lungs were used in plethysmographic measurements of lung resistance and elastance at different lung volumes, whereas the airway dimensions were measured by computed tomography (CT). The change in airway diameter informed by CT measurements was applied to isolated airway ring preparations to determine the strain-induced loss of ASM contractility. We found that changing the transpulmonary pressure from 5 to 30 cmH2O led to a 51% increase in lung volume, accompanied by a 46% increase in the airway diameter with no change in airway heterogeneity. When comparable airway strains measured in the whole lung were applied to isolated airway rings in either relaxed or contracted state, a significant loss of ASM contractility was observed, suggesting that DI-induced bronchodilation and bronchoprotection can result from strain-induced loss of ASM contractility.

Rationale: Although centrilobular emphysema (CLE) and paraseptal emphysema (PSE) are commonly identified on multidetector computed tomography (MDCT), little is known about the pathology associated with PSE compared with that of CLE.Objectives: To assess the pathological differences between PSE and CLE in chronic obstructive pulmonary disease (COPD).Methods: Air-inflated frozen lung specimens (n = 6) obtained from patients with severe COPD treated by lung transplantation were scanned with MDCT. Frozen tissue cores were taken from central (n = 8) and peripheral (n = 8) regions of each lung, scanned with micro-computed tomography (microCT), and processed for histology. The core locations were registered to the MDCT, and a percentage of PSE or CLE was assigned by radiologists to each of the regions. MicroCT scans were used to measure number and structural change of terminal bronchioles. Furthermore, microCT-based volume fractions of CLE and PSE allowed classifying cores into mild emphysema, CLE-dominant, and PSE-dominant.Measurements and Main Results: The percentages of PSE measured on MDCT and microCT were positively associated (P = 0.015). The number of terminal bronchioles per milliliter of lung and cross-sectional lumen area were significantly lower and wall area percentage was significantly higher in CLE-dominant regions compared with mild emphysema and PSE-dominant regions (all P < 0.05), whereas no difference was found between PSE-dominant and mild emphysema samples (all P > 0.5). Immunohistochemistry showed significantly higher infiltration of neutrophils (P = 0.002), but not of macrophages, CD4, CD8, or B cells, in PSE compared with CLE regions.Conclusions: The terminal bronchioles are relatively preserved, whereas neutrophilic inflammation is increased in PSE-dominant regions compared with CLE-dominant regions in patients with COPD.

Patients who have chronic obstructive pulmonary disease (COPD) and bronchial asthma (BA) share symptoms such as, dyspnoea, cough and wheeze. Differentiating these diseases in the ambulatory setting can be challenging especially in older adult smokers who are being treated with a variety of medications. The objective of this study was to test the value of adding a maximal inspiratory manoeuvre to basic spirometry to differentiate COPD and BA. One hundred forty-three COPD patients and 142 BA patients had measurements of maximal inspiratory and expiratory flow during routine spirometry. Parameters from these tests were used to assess diagnostic accuracy using receiver-operating characteristic (ROC) analyses followed by logistic regression. The association of two independent parameters were analyzed using linear regression analyses. Results show that forced expiratory volume in one second/forced vital capacity (FEV1/FVC%) <62.4 was the best independent predictor to diagnose COPD. The combination of FEV1/FVC% <62.4 and the ratio of peak inspiratory flow/maximal expiratory flow at 50% FVC (PIF/MEF50) >3.06 significantly predicted COPD. Post-test probability for prediction of COPD was 82.0% when patients had both parameters. When asthmatic patients with a smoking history were compared with COPD patients, FEV1/FVC% <63.4 and PIF/MEF50 >3.29 were both independent predictors of COPD. The post-test probability for COPD was 94.4% when patients had both parameters. The association between FEV1/FVC% and PIF/MEF50 was significantly different between COPD and BA. In conclusion, the addition of the maximal inspiratory effort to routine pulmonary function measurements provides a simple test to help differentiate COPD and BA.

The lack of bronchodilatory response to deep inspiration in asthmatics is thought to be partially due to reduced airway distensibility [1, 2], possibly caused by an increase in airway smooth muscle (ASM) tone and stiffness [3]. Rho-kinase (ROCK) is known to play a role in regulating ASM tone [4] and ASM cytoskeletal stiffness [5]. An important phenotypic change is associated with human asthma: the protein expression of both isoforms of Rho-kinase, ROCK1 and ROCK2, is increased in the smooth muscles of intra-parenchymal airways and in the pulmonary blood vessels http://bit.ly/2pPfyxH

Asthmatic airways are stiffer than normal. We have shown that the cytoskeletal passive stiffness of airway smooth muscle (ASM) can be regulated by intracellular signaling pathways, especially those associated with Rho kinase (ROCK). We have also shown that an oscillatory strain reduces the passive stiffness of ASM and its ability to generate force. Here, we investigated the combined effect of inhibiting the ASM contraction with β2 agonist and decreasing the ASM cytoskeletal stiffness with ROCK inhibitor and/or force oscillation (FO) on the relaxation of contracted ASM. We hypothesize that the ASM relaxation can be synergistically enhanced by the combination of these interventions, because drug-induced softening of the cytoskeleton enhances the FO-induced relaxation and vice versa. Sheep tracheal strips were isotonically contracted to acetylcholine (3 × 10−5 M). At the plateau of shortening, β2 agonist salbutamol (10−7 M), ROCK inhibitor H1152 (10−7 M), and FO (square wave, 1 Hz, amplitude 6% maximal active force) were applied either alone or in combination. After adjusting for nonspecific time-dependent variation, relengthening by individual interventions with low-dose salbutamol or H1152, or small amplitude FO was not significantly different from zero. However, significant relengthening was observed in all combination treatments. The relengthening was greater than the mathematical sum of relengthening caused by individual treatments thereby demonstrating synergistic relaxation. The ASM stiffness did not change with salbutamol or H1152 treatments, but was lower with FO in combination with H1152. The results suggest that the mechanopharmacological treatment can be an effective therapy for asthma.

Sex strongly affects the gene expression profile in lung tissue, and overall reflects the activity of different biological pathways that probably play a role in the susceptibility to develop COPD with smokinghttp://bit.ly/31drVkX

The cyclic interaction between myosin crossbridges and actin filaments underlies smooth muscle contraction. Phosphorylation of the 20-kDa myosin light chain (MLC20) is a crucial step in activating the crossbridge cycle. Our current understanding of smooth muscle contraction is based on observed correlations among MLC20 phosphorylation, maximal shortening velocity (Vmax), and isometric force over the time course of contraction. However, during contraction there are changes in the extent of phosphorylation of many additional proteins as well as changes in activation of enzymes associated with the signaling pathways. As a consequence, the mechanical manifestation of muscle contraction is likely to change with time. To simplify the study of these relationships, we measured the mechanical properties of airway smooth muscle at different levels of MLC20 phosphorylation at a fixed time during contraction. A simple correlation emerged when time-dependent variables were fixed. MLC20 phosphorylation was found to be directly and linearly correlated with the active stress, stiffness, and power of the muscle; the observed weak dependence of Vmax on MLC20 phosphorylation could be explained by the presence of an internal load in the muscle preparation. These results can be entirely explained by the Huxley crossbridge model. We conclude that when the influence of time-dependent events during contraction is held constant, the basic crossbridge mechanism in smooth muscle is the same as that in striated muscle.

Inhaled corticosteroids (ICS) are widely prescribed for patients with chronic obstructive pulmonary disease (COPD), yet have variable outcomes and adverse reactions, which may be genetically determined. The primary aim of the study was to identify the genetic determinants for forced expiratory volume in 1 s (FEV 1 ) changes related to ICS therapy. In the Lung Health Study (LHS)-2, 1116 COPD patients were randomised to the ICS triamcinolone acetonide (n=559) or placebo (n=557) with spirometry performed every 6 months for 3 years. We performed a pharmacogenomic genome-wide association study for the genotype-by-ICS treatment effect on 3 years of FEV 1 changes (estimated as slope) in 802 genotyped LHS-2 participants. Replication was performed in 199 COPD patients randomised to the ICS, fluticasone or placebo. A total of five loci showed genotype-by-ICS interaction at p<5×10 −6 ; of these, single nucleotide polymorphism (SNP) rs111720447 on chromosome 7 was replicated (discovery p=4.8×10 −6 , replication p=5.9×10 −5 ) with the same direction of interaction effect. ENCODE (Encyclopedia of DNA Elements) data revealed that in glucocorticoid-treated (dexamethasone) A549 alveolar cell line, glucocorticoid receptor binding sites were located near SNP rs111720447. In stratified analyses of LHS-2, genotype at SNP rs111720447 was significantly associated with rate of FEV 1 decline in patients taking ICS (C allele β 56.36 mL·year −1 , 95% CI 29.96–82.76 mL·year −1 ) and in patients who were assigned to placebo, although the relationship was weaker and in the opposite direction to that in the ICS group (C allele β −27.57 mL·year −1 , 95% CI −53.27– −1.87 mL·year −1 ). The study uncovered genetic factors associated with FEV 1 changes related to ICS in COPD patients, which may provide new insight on the potential biology of steroid responsiveness in COPD.

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. 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). 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. 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.

BACKGROUND:Smoking is the principal modifiable environmental risk factor for chronic obstructive pulmonary disease (COPD) which affects 300 million people and is the 3rd leading cause of death worldwide. Most of the genetic studies of smoking have relied on self-reported smoking status which is vulnerable to reporting and recall bias. Using data from the Lung Health Study (LHS), we sought to identify genetic variants associated with quantitative smoking and cessation in individuals with mild to moderate COPD. METHODS:The LHS is a longitudinal multicenter study of mild-to-moderate COPD subjects who were all smokers at recruitment. We performed genome-wide association studies (GWASs) for salivary cotinine (n = 4024), exhaled carbon monoxide (eCO) (n = 2854), cigarettes per day (CPD) (n = 2706) and smoking cessation at year 5 follow-up (n = 717 quitters and 2175 smokers). The GWAS analyses were adjusted for age, gender, and genetic principal components. RESULTS:For cotinine levels, SNPs near UGT2B10 gene achieved genome-wide significance (i.e. P < 5 × 10- 8) with top SNP rs10023464, P = 1.27 × 10- 11. For eCO levels, one significant SNP was identified which mapped to the CHRNA3 gene (rs12914385, P = 2.38 × 10- 8). A borderline region mapping to KCNMA1 gene was associated with smoking cessation (rs207675, P = 5.95 × 10- 8). Of the identified loci, only the CHRNA3/5 locus showed significant associations with lung function but only in heavy smokers. No regions met genome-wide significance for CPD. CONCLUSION:The study demonstrates that using objective measures of smoking such as eCO and/or salivary cotinine can more precisely capture the genetic contribution to multiple aspects of smoking behaviour. The KCNMA1 gene association with smoking cessation may represent a potential therapeutic target and warrants further studies. TRIAL REGISTRATION:The Lung Health Study ClinicalTrials.gov Identifier: NCT00000568 . Date of registration: October 28, 1999.