Dynamic and prognostic proteomic associations with FEV1 decline in chronic obstructive pulmonary disease.

Rationale: Identification and validation of circulating biomarkers for lung function decline in COPD remains an unmet need. Objective: Identify prognostic and dynamic plasma protein biomarkers of COPD progression. Methods: We measured plasma proteins using SomaScan from two COPD-enriched cohorts, the Subpopulations and Intermediate Outcomes Measures in COPD Study (SPIROMICS) and Genetic Epidemiology of COPD (COPDGene), and one population-based cohort, Multi-Ethnic Study of Atherosclerosis (MESA) Lung. Using SPIROMICS as a discovery cohort, linear mixed models identified baseline proteins that predicted future change in FEV1 (prognostic model) and proteins whose expression changed with change in lung function (dynamic model). Findings were replicated in COPDGene and MESA-Lung. Using the COPD-enriched cohorts, Gene Set Enrichment Analysis (GSEA) identified proteins shared between COPDGene and SPIROMICS. Metascape identified significant associated pathways. Measurements and Main Results: The prognostic model found 7 significant proteins in common (p < 0.05) among all 3 cohorts. After applying false discovery rate (adjusted p < 0.2), leptin remained significant in all three cohorts and growth hormone receptor remained significant in the two COPD cohorts. Elevated baseline levels of leptin and growth hormone receptor were associated with slower rate of decline in FEV1. Twelve proteins were nominally but not FDR significant in the dynamic model and all were distinct from the prognostic model. Metascape identified several immune related pathways unique to prognostic and dynamic proteins. Conclusion: We identified leptin as the most reproducible COPD progression biomarker. The difference between prognostic and dynamic proteins suggests disease activity signatures may be different from prognosis signatures.

Design of the SPIROMICS Study of Early COPD Progression: SOURCE Study.

The biological mechanisms leading some tobacco-exposed individuals to develop early-stage chronic obstructive pulmonary disease (COPD) are poorly understood. This knowledge gap hampers development of disease-modifying agents for this prevalent condition. Accord-ingly, with National Heart, Lung and Blood Institute support, we initiated the SPIROMICS Study of Early COPD Progression (SOURCE), a multicenter observational cohort study of younger individuals with a history of cigarette smoking and thus at-risk for, or with, early-stage COPD. Our overall objectives are to identify those who will develop COPD earlier in life, characterize them thoroughly, and by contrasting them to those not developing COPD, define mechanisms of disease progression. SOURCE utilizes the established SPIROMICS clinical network. Its goal is to enroll n=649 participants, ages 30-55 years, all races/ethnicities, with ≥10 pack-years cigarette smoking, in either Global Initiative for Chronic Obstructive Lung Disease (GOLD) groups 0-2 or with Preserved Ratio Impaired Spirometry (PRISm); and an additional n=40 never-smoker controls. Participants undergo baseline and three-year follow-up visits, each including high-resolution computed tomography; respiratory oscillometry and spirometry (pre- and post-bronchodilator administration), exhaled breath condensate (baseline only); and extensive biospecimen collection, including sputum induction. Symptoms, interim healthcare utilization, and exacerbations are captured every six months via follow-up phone calls. An embedded bronchoscopy sub-study involving n=100 participants (including all never-smokers) will allow collection of lower airway samples for genetic, epigenetic, genomic, immunological, microbiome, mucin analyses, and basal cell culture. SOURCE should provide novel insights into the natural history of lung disease in younger individuals with a smoking history, and its biological basis.

Dysregulated Airway Host Defense in Hyper IgE Syndrome due to STAT3 Mutations.

Rationale: Hyper IgE syndrome (STAT3-HIES), also known as Job's syndrome, is a rare immunodeficiency disease typically caused by dominant-negative STAT3 mutations. STAT3-HIES syndrome is characterized by chronic pulmonary infection and inflammation, suggesting impairment of pulmonary innate host defense. Objectives: To identify airway epithelial host defense defects consequent to STAT3 mutations that, in addition to reported mutant STAT3 immunologic abnormalities, produce pulmonary infection. Methods: STAT3-HIES sputum was evaluated for biochemical/biophysical properties. STAT3-HIES excised lungs were harvested for histology; bronchial brush samples were collected for RNA sequencing and in vitro culture. A STAT3-HIES-specific mutation (R382W), expressed by lentiviruses, and a STAT3 knockout, generated by CRISPR/Cas9, were maintained in normal human bronchial epithelia under basal or inflammatory (IL1ß) conditions. Effects of STAT3 deficiency on transcriptomics, and epithelial ion channel, secretory, antimicrobial, and ciliary functions were assessed. Measurements and Main Results: Mucus concentrations and viscoelasticity were increased in STAT3-HIES sputum. STAT3-HIES excised lungs exhibited mucus obstruction and elevated IL1ß expression. STAT3 deficiency impaired CFTR-dependent fluid and mucin secretion, inhibited expression of antimicrobial peptides, cytokines, and chemokines, and acidified airway surface liquid at baseline and post-IL1ß exposure in vitro. Notably, mutant STAT3 suppressed IL1R1 expression. STAT3 mutations also inhibited ciliogenesis in vivo and impaired mucociliary transport in vitro, a process mediated via HES6 suppression. Administration of a γ-secretase inhibitor increased HES6 expression and improved ciliogenesis in STAT3 R382W mutant cells. Conclusions: STAT3 dysfunction leads to multi-component defects in airway epithelial innate defense, which, in conjunction with STAT3-HIES immune deficiency, contributes to chronic pulmonary infection.

Structural and genetic diversity in the secreted mucins MUC5AC and MUC5B.

The secreted mucins MUC5AC and MUC5B are large glycoproteins that play critical defensive roles in pathogen entrapment and mucociliary clearance. Their respective genes contain polymorphic and degenerate protein-coding variable number tandem repeats (VNTRs) that make the loci difficult to investigate with short reads. We characterize the structural diversity of MUC5AC and MUC5B by long-read sequencing and assembly of 206 human and 20 nonhuman primate (NHP) haplotypes. We find that human MUC5B is largely invariant (5,761-5,762 amino acids [aa]); however, seven haplotypes have expanded VNTRs (6,291-7,019 aa). In contrast, 30 allelic variants of MUC5AC encode 16 distinct proteins (5,249-6,325 aa) with cysteine-rich domain and VNTR copy-number variation. We group MUC5AC alleles into three phylogenetic clades: H1 (46%, ∼5,654 aa), H2 (33%, ∼5,742 aa), and H3 (7%, ∼6,325 aa). The two most common human MUC5AC variants are smaller than NHP gene models, suggesting a reduction in protein length during recent human evolution. Linkage disequilibrium and Tajima's D analyses reveal that East Asians carry exceptionally large blocks with an excess of rare variation (p < 0.05) at MUC5AC. To validate this result, we use Locityper for genotyping MUC5AC haplogroups in 2,600 unrelated samples from the 1000 Genomes Project. We observe a signature of positive selection in H1 among East Asians and a depletion of the likely ancestral haplogroup (H3). In Europeans, H3 alleles show an excess of common variation and deviate from Hardy-Weinberg equilibrium (p < 0.05), consistent with heterozygote advantage and balancing selection. This study provides a generalizable strategy to characterize complex protein-coding VNTRs for improved disease associations.

The clinical utility of sequencing the entirety of CFTR.

BACKGROUND: Cystic fibrosis (CF) is caused by deleterious variants in each CFTR gene. We investigated the utility of whole-gene CFTR sequencing when fewer than two pathogenic or likely pathogenic (P/LP) variants were detected by conventional testing (sequencing of exons and flanking introns) of CFTR. METHODS: Individuals with features of CF and a CF-diagnostic sweat chloride concentration with zero or one P/LP variants identified by conventional testing enrolled in the CF Mutation Analysis Program (MAP) underwent whole-gene CFTR sequencing. Replication was performed on individuals enrolled in the CF Genome Project (CFGP), followed by phenotype review and interrogation of other genes. RESULTS: Whole-gene sequencing identified a second P/LP variant in 20/43 MAP enrollees (47 %) and 10/22 CFGP enrollees (45 %) who had one P/LP variant after conventional testing. No P/LP variants were detected when conventional testing was negative (MAP: n = 43; CFGP: n = 13). Genome-wide analysis was unable to find an alternative etiology in CFGP participants with fewer than two P/LP CFTR variants and CF could not be confirmed in 91 % following phenotype re-review. CONCLUSIONS: Whole-gene CFTR analysis is beneficial in individuals with one previously-identified P/LP variant and a CF-diagnostic sweat chloride. Negative conventional CFTR testing indicates that the phenotype should be re-evaluated.

MagicalRsq-X: A cross-cohort transferable genotype imputation quality metric.

Since genotype imputation was introduced, researchers have been relying on the estimated imputation quality from imputation software to perform post-imputation quality control (QC). However, this quality estimate (denoted as Rsq) performs less well for lower-frequency variants. We recently published MagicalRsq, a machine-learning-based imputation quality calibration, which leverages additional typed markers from the same cohort and outperforms Rsq as a QC metric. In this work, we extended the original MagicalRsq to allow cross-cohort model training and named the new model MagicalRsq-X. We removed the cohort-specific estimated minor allele frequency and included linkage disequilibrium scores and recombination rates as additional features. Leveraging whole-genome sequencing data from TOPMed, specifically participants in the BioMe, JHS, WHI, and MESA studies, we performed comprehensive cross-cohort evaluations for predominantly European and African ancestral individuals based on their inferred global ancestry with the 1000 Genomes and Human Genome Diversity Project data as reference. Our results suggest MagicalRsq-X outperforms Rsq in almost every setting, with 7.3%-14.4% improvement in squared Pearson correlation with true R2, corresponding to 85-218 K variant gains. We further developed a metric to quantify the genetic distances of a target cohort relative to a reference cohort and showed that such metric largely explained the performance of MagicalRsq-X models. Finally, we found MagicalRsq-X saved up to 53 known genome-wide significant variants in one of the largest blood cell trait GWASs that would be missed using the original Rsq for QC. In conclusion, MagicalRsq-X shows superiority for post-imputation QC and benefits genetic studies by distinguishing well and poorly imputed lower-frequency variants.

Structural and genetic diversity in the secreted mucins, MUC5AC and MUC5B.

The secreted mucins MUC5AC and MUC5B play critical defensive roles in airway pathogen entrapment and mucociliary clearance by encoding large glycoproteins with variable number tandem repeats (VNTRs). These polymorphic and degenerate protein coding VNTRs make the loci difficult to investigate with short reads. We characterize the structural diversity of MUC5AC and MUC5B by long-read sequencing and assembly of 206 human and 20 nonhuman primate (NHP) haplotypes. We find that human MUC5B is largely invariant (5761-5762aa); however, seven haplotypes have expanded VNTRs (6291-7019aa). In contrast, 30 allelic variants of MUC5AC encode 16 distinct proteins (5249-6325aa) with cysteine-rich domain and VNTR copy number variation. We grouped MUC5AC alleles into three phylogenetic clades: H1 (46%, ~5654aa), H2 (33%, ~5742aa), and H3 (7%, ~6325aa). The two most common human MUC5AC variants are smaller than NHP gene models, suggesting a reduction in protein length during recent human evolution. Linkage disequilibrium (LD) and Tajima's D analyses reveal that East Asians carry exceptionally large MUC5AC LD blocks with an excess of rare variation (p<0.05). To validate this result, we used Locityper for genotyping MUC5AC haplogroups in 2,600 unrelated samples from the 1000 Genomes Project. We observed signatures of positive selection in H1 and H2 among East Asians and a depletion of the likely ancestral haplogroup (H3). In Africans and Europeans, H3 alleles show an excess of common variation and deviate from Hardy-Weinberg equilibrium, consistent with heterozygote advantage and balancing selection. This study provides a generalizable strategy to characterize complex protein coding VNTRs for improved disease associations.

Genetic variation in severe cystic fibrosis liver disease is associated with novel mechanisms for disease pathogenesis.

BACKGROUND AND AIMS: It is not known why severe cystic fibrosis (CF) liver disease (CFLD) with portal hypertension occurs in only ~7% of people with CF. We aimed to identify genetic modifiers for severe CFLD to improve understanding of disease mechanisms. APPROACH AND RESULTS: Whole-genome sequencing was available in 4082 people with CF with pancreatic insufficiency (n = 516 with severe CFLD; n = 3566 without CFLD). We tested ~15.9 million single nucleotide polymorphisms (SNPs) for association with severe CFLD versus no-CFLD, using pre-modulator clinical phenotypes including (1) genetic variant ( SERPINA1 ; Z allele) previously associated with severe CFLD; (2) candidate SNPs (n = 205) associated with non-CF liver diseases; (3) genome-wide association study of common/rare SNPs; (4) transcriptome-wide association; and (5) gene-level and pathway analyses. The Z allele was significantly associated with severe CFLD ( p = 1.1 × 10 -4 ). No significant candidate SNPs were identified. A genome-wide association study identified genome-wide significant SNPs in 2 loci and 2 suggestive loci. These 4 loci contained genes [significant, PKD1 ( p = 8.05 × 10 -10 ) and FNBP1 ( p = 4.74 × 10 -9 ); suggestive, DUSP6 ( p = 1.51 × 10 -7 ) and ANKUB1 ( p = 4.69 × 10 -7 )] relevant to severe CFLD pathophysiology. The transcriptome-wide association identified 3 genes [ CXCR1 ( p = 1.01 × 10 -6 ) , AAMP ( p = 1.07 × 10 -6 ), and TRBV24 ( p = 1.23 × 10 -5 )] involved in hepatic inflammation and innate immunity. Gene-ranked analyses identified pathways enriched in genes linked to multiple liver pathologies. CONCLUSION: These results identify loci/genes associated with severe CFLD that point to disease mechanisms involving hepatic fibrosis, inflammation, innate immune function, vascular pathology, intracellular signaling, actin cytoskeleton and tight junction integrity and mechanisms of hepatic steatosis and insulin resistance. These discoveries will facilitate mechanistic studies and the development of therapeutics for severe CFLD.

Prolonged airway explant culture enables study of health, disease, and viral pathogenesis.

In vitro models play a major role in studying airway physiology and disease. However, the native lung's complex tissue architecture and non-epithelial cell lineages are not preserved in these models. Ex vivo tissue models could overcome in vitro limitations, but methods for long-term maintenance of ex vivo tissue has not been established. We describe methods to culture human large airway explants, small airway explants, and precision-cut lung slices for at least 14 days. Human airway explants recapitulate genotype-specific electrophysiology, characteristic epithelial, endothelial, stromal and immune cell populations, and model viral infection after 14 days in culture. These methods also maintain mouse, rabbit, and pig tracheal explants. Notably, intact airway tissue can be cryopreserved, thawed, and used to generate explants with recovery of function 14 days post-thaw. These studies highlight the broad applications of airway tissue explants and their use as translational intermediates between in vitro and in vivo studies.

Proximal and Distal Bronchioles Contribute to the Pathogenesis of Non-Cystic Fibrosis Bronchiectasis.

Rationale: Non-cystic fibrosis bronchiectasis (NCFB) may originate in bronchiolar regions of the lung. Accordingly, there is a need to characterize the morphology and molecular characteristics of NCFB bronchioles. Objectives: Test the hypothesis that NCFB exhibits a major component of bronchiolar disease manifest by mucus plugging and ectasia. Methods: Morphologic criteria and region-specific epithelial gene expression, measured histologically and by RNA in situ hybridization and immunohistochemistry, identified proximal and distal bronchioles in excised NCFB lungs. RNA in situ hybridization and immunohistochemistry assessed bronchiolar mucus accumulation and mucin gene expression. CRISPR-Cas9-mediated IL-1R1 knockout in human bronchial epithelial cultures tested IL-1α and IL-1ß contributions to mucin production. Spatial transcriptional profiling characterized NCFB distal bronchiolar gene expression. Measurements and Main Results: Bronchiolar perimeters and lumen areas per section area were increased in proximal, but not distal, bronchioles in NCFB versus control lungs, suggesting proximal bronchiolectasis. In NCFB, mucus plugging was observed in ectatic proximal bronchioles and associated nonectatic distal bronchioles in sections with disease. MUC5AC and MUC5B mucins were upregulated in NCFB proximal bronchioles, whereas MUC5B was selectively upregulated in distal bronchioles. Bronchiolar mucus plugs were populated by IL-1ß-expressing macrophages. NCFB sterile sputum supernatants induced human bronchial epithelial MUC5B and MUC5AC expression that was >80% blocked by IL-1R1 ablation. Spatial transcriptional profiling identified upregulation of genes associated with secretory cells, hypoxia, interleukin pathways, and IL-1ß-producing macrophages in mucus plugs and downregulation of epithelial ciliogenesis genes. Conclusions: NCFB exhibits distinctive proximal and distal bronchiolar disease. Both bronchiolar regions exhibit bronchiolar secretory cell features and mucus plugging but differ in mucin gene regulation and ectasia.

Chronic airway epithelial hypoxia exacerbates injury in muco-obstructive lung disease through mucus hyperconcentration.

Unlike solid organs, human airway epithelia derive their oxygen from inspired air rather than the vasculature. Many pulmonary diseases are associated with intraluminal airway obstruction caused by aspirated foreign bodies, virus infection, tumors, or mucus plugs intrinsic to airway disease, including cystic fibrosis (CF). Consistent with requirements for luminal O2, airway epithelia surrounding mucus plugs in chronic obstructive pulmonary disease (COPD) lungs are hypoxic. Despite these observations, the effects of chronic hypoxia (CH) on airway epithelial host defense functions relevant to pulmonary disease have not been investigated. Molecular characterization of resected human lungs from individuals with a spectrum of muco-obstructive lung diseases (MOLDs) or COVID-19 identified molecular features of chronic hypoxia, including increased EGLN3 expression, in epithelia lining mucus-obstructed airways. In vitro experiments using cultured chronically hypoxic airway epithelia revealed conversion to a glycolytic metabolic state with maintenance of cellular architecture. Chronically hypoxic airway epithelia unexpectedly exhibited increased MUC5B mucin production and increased transepithelial Na+ and fluid absorption mediated by HIF1α/HIF2α-dependent up-regulation of ß and γENaC (epithelial Na+ channel) subunit expression. The combination of increased Na+ absorption and MUC5B production generated hyperconcentrated mucus predicted to perpetuate obstruction. Single-cell and bulk RNA sequencing analyses of chronically hypoxic cultured airway epithelia revealed transcriptional changes involved in airway wall remodeling, destruction, and angiogenesis. These results were confirmed by RNA-in situ hybridization studies of lungs from individuals with MOLD. Our data suggest that chronic airway epithelial hypoxia may be central to the pathogenesis of persistent mucus accumulation in MOLDs and associated airway wall damage.

Genetic regulators of sputum mucin concentration and their associations with COPD phenotypes.

Hyper-secretion and/or hyper-concentration of mucus is a defining feature of multiple obstructive lung diseases, including chronic obstructive pulmonary disease (COPD). Mucus itself is composed of a mixture of water, ions, salt and proteins, of which the gel-forming mucins, MUC5AC and MUC5B, are the most abundant. Recent studies have linked the concentrations of these proteins in sputum to COPD phenotypes, including chronic bronchitis (CB) and acute exacerbations (AE). We sought to determine whether common genetic variants influence sputum mucin concentrations and whether these variants are also associated with COPD phenotypes, specifically CB and AE. We performed a GWAS to identify quantitative trait loci for sputum mucin protein concentration (pQTL) in the Sub-Populations and InteRmediate Outcome Measures in COPD Study (SPIROMICS, n = 708 for total mucin, n = 215 for MUC5AC, MUC5B). Subsequently, we tested for associations of mucin pQTL with CB and AE using regression modeling (n = 822-1300). Replication analysis was conducted using data from COPDGene (n = 5740) and by examining results from the UK Biobank. We identified one genome-wide significant pQTL for MUC5AC (rs75401036) and two for MUC5B (rs140324259, rs10001928). The strongest association for MUC5B, with rs140324259 on chromosome 11, explained 14% of variation in sputum MUC5B. Despite being associated with lower MUC5B, the C allele of rs140324259 conferred increased risk of CB (odds ratio (OR) = 1.42; 95% confidence interval (CI): 1.10-1.80) as well as AE ascertained over three years of follow up (OR = 1.41; 95% CI: 1.02-1.94). Associations between rs140324259 and CB or AE did not replicate in COPDGene. However, in the UK Biobank, rs140324259 was associated with phenotypes that define CB, namely chronic mucus production and cough, again with the C allele conferring increased risk. We conclude that sputum MUC5AC and MUC5B concentrations are associated with common genetic variants, and the top locus for MUC5B may influence COPD phenotypes, in particular CB.

Large scale proteomic studies create novel privacy considerations.

Privacy protection is a core principle of genomic but not proteomic research. We identified independent single nucleotide polymorphism (SNP) quantitative trait loci (pQTL) from COPDGene and Jackson Heart Study (JHS), calculated continuous protein level genotype probabilities, and then applied a naïve Bayesian approach to link SomaScan 1.3K proteomes to genomes for 2812 independent subjects from COPDGene, JHS, SubPopulations and InteRmediate Outcome Measures In COPD Study (SPIROMICS) and Multi-Ethnic Study of Atherosclerosis (MESA). We correctly linked 90-95% of proteomes to their correct genome and for 95-99% we identify the 1% most likely links. The linking accuracy in subjects with African ancestry was lower (~ 60%) unless training included diverse subjects. With larger profiling (SomaScan 5K) in the Atherosclerosis Risk Communities (ARIC) correct identification was > 99% even in mixed ancestry populations. We also linked proteomes-to-proteomes and used the proteome only to determine features such as sex, ancestry, and first-degree relatives. When serial proteomes are available, the linking algorithm can be used to identify and correct mislabeled samples. This work also demonstrates the importance of including diverse populations in omics research and that large proteomic datasets (> 1000 proteins) can be accurately linked to a specific genome through pQTL knowledge and should not be considered unidentifiable.

Impact of Bronchiectasis on COPD Severity and Alpha-1 Antitrypsin Deficiency as a Risk Factor in Individuals with a Heavy Smoking History.

Rationale: Bronchiectasis is common among those with heavy smoking histories, but risk factors for bronchiectasis, including alpha-1 antitrypsin deficiency, and its implications for COPD severity are uncharacterized in such individuals. Objectives: To characterize the impact of bronchiectasis on COPD and explore alpha-1antitrypsin as a risk factor for bronchiectasis. Methods: SubPopulations and InteRmediate Outcome Measures In COPD Study (SPIROMICS) participants (N=914; ages 40-80 years; ≥20-pack-year smoking) had high-resolution computed tomography (CT) scans interpreted visually for bronchiectasis, based on airway dilation without fibrosis or cicatrization. We performed regression-based models of bronchiectasis with clinical outcomes and quantitative CT measures. We deeply sequenced the gene encoding -alpha-1 antitrypsin, SERPINA1, in 835 participants to test for rare variants, focusing on the PiZ genotype (Glu366Lys, rs28929474). Measurements and Main Results: We identified bronchiectasis in 365 (40%) participants, more frequently in women (45% versus 36%, p=0.0045), older participants (mean age=66[standard deviation (SD)=8.3] versus 64[SD=9.1] years, p=0.0083), and those with lower lung function (forced expiratory volume in 1 second [FEV1 ] percentage predicted=66%[SD=27] versus 77%[SD=25], p<0.0001; FEV1 to forced vital capacity [FVC] ratio=0.54[0.17] versus 0.63[SD=0.16], p<0.0001). Participants with bronchiectasis had greater emphysema (%voxels ≤-950 Hounsfield units, 11%[SD=12] versus 6.3%[SD=9], p<0.0001) and parametric response mapping functional small airways disease (26[SD=15] versus 19[SD=15], p<0.0001). Bronchiectasis was more frequent in the combined PiZZ and PiMZ genotype groups compared to those without PiZ, PiS, or other rare pathogenic variants (N=21 of 40 [52%] versus N=283 of 707[40%], odds ratio [OR]=1.97; 95% confidence interval [CI]=1.002, 3.90, p=0.049), an association attributed to White individuals (OR=1.98; 95%CI = 0.9956, 3.9; p=0.051). Conclusions: Bronchiectasis was common in those with heavy smoking histories and was associated with detrimental clinical and radiographic outcomes. Our findings support alpha-1antitrypsin guideline recommendations to screen for alpha-1 antitrypsin deficiency in an appropriate bronchiectasis subgroup with a significant smoking history.

A blood and bronchoalveolar lavage protein signature of rapid FEV1 decline in smoking-associated COPD.

Accelerated progression of chronic obstructive pulmonary disease (COPD) is associated with increased risks of hospitalization and death. Prognostic insights into mechanisms and markers of progression could facilitate development of disease-modifying therapies. Although individual biomarkers exhibit some predictive value, performance is modest and their univariate nature limits network-level insights. To overcome these limitations and gain insights into early pathways associated with rapid progression, we measured 1305 peripheral blood and 48 bronchoalveolar lavage proteins in individuals with COPD [n = 45, mean initial forced expiratory volume in one second (FEV1) 75.6 ± 17.4% predicted]. We applied a data-driven analysis pipeline, which enabled identification of protein signatures that predicted individuals at-risk for accelerated lung function decline (FEV1 decline ≥ 70 mL/year) ~ 6 years later, with high accuracy. Progression signatures suggested that early dysregulation in elements of the complement cascade is associated with accelerated decline. Our results propose potential biomarkers and early aberrant signaling mechanisms driving rapid progression in COPD.

Reduced quantity and function of pneumococcal antibodies are associated with exacerbations of COPD in SPIROMICS.

While hypogammaglobulinemia is associated with COPD exacerbations, it is unknown whether frequent exacerbators have specific defects in antibody production/function. We hypothesized that reduced quantity/function of serum pneumococcal antibodies correlate with exacerbation risk in the SPIROMICS cohort. We measured total pneumococcal IgG in n = 764 previously vaccinated participants with COPD. In a propensity-matched subset of n = 200 with vaccination within five years (n = 50 without exacerbations in the previous year; n = 75 with one, n = 75 with ≥2), we measured pneumococcal IgG for 23 individual serotypes, and pneumococcal antibody function for 4 serotypes. Higher total pneumococcal IgG, serotype-specific IgG (17/23 serotypes), and antibody function (3/4 serotypes) were independently associated with fewer prior exacerbations. Higher pneumococcal IgG (5/23 serotypes) predicted lower exacerbation risk in the following year. Pneumococcal antibodies are inversely associated with exacerbations, supporting the presence of immune defects in frequent exacerbators. With further study, pneumococcal antibodies may be useful biomarkers for immune dysfunction in COPD.

Genetic Modifiers of Cystic Fibrosis Lung Disease Severity: Whole-Genome Analysis of 7,840 Patients.

Rationale: Lung disease is the major cause of morbidity and mortality in persons with cystic fibrosis (pwCF). Variability in CF lung disease has substantial non-CFTR (CF transmembrane conductance regulator) genetic influence. Identification of genetic modifiers has prognostic and therapeutic importance. Objectives: Identify genetic modifier loci and genes/pathways associated with pulmonary disease severity. Methods: Whole-genome sequencing data on 4,248 unique pwCF with pancreatic insufficiency and lung function measures were combined with imputed genotypes from an additional 3,592 patients with pancreatic insufficiency from the United States, Canada, and France. This report describes association of approximately 15.9 million SNPs using the quantitative Kulich normal residual mortality-adjusted (KNoRMA) lung disease phenotype in 7,840 pwCF using premodulator lung function data. Measurements and Main Results: Testing included common and rare SNPs, transcriptome-wide association, gene-level, and pathway analyses. Pathway analyses identified novel associations with genes that have key roles in organ development, and we hypothesize that these genes may relate to dysanapsis and/or variability in lung repair. Results confirmed and extended previous genome-wide association study findings. These whole-genome sequencing data provide finely mapped genetic information to support mechanistic studies. No novel primary associations with common single variants or rare variants were found. Multilocus effects at chr5p13 (SLC9A3/CEP72) and chr11p13 (EHF/APIP) were identified. Variant effect size estimates at associated loci were consistently ordered across the cohorts, indicating possible age or birth cohort effects. Conclusions: This premodulator genomic, transcriptomic, and pathway association study of 7,840 pwCF will facilitate mechanistic and postmodulator genetic studies and the development of novel therapeutics for CF lung disease.

Prolonged, physiologically relevant nicotine concentrations in the airways of smokers.

Nicotine from cigarette smoke is a biologically active molecule that has pleiotropic effects in the airway, which could play a role in smoking-induced lung disease. However, whether nicotine and its metabolites reach sustained, physiologically relevant concentrations on airway surfaces of smokers is not well defined. To address these issues, concentrations of nicotine, cotinine, and hydroxycotinine were measured by mass spectrometry (MS) in supernatants of induced sputum obtained from participants in the subpopulations and intermediate outcome measures in COPD study (SPIROMICS), an ongoing observational study that included never smokers, former smokers, and current smokers with and without chronic obstructive pulmonary disease (COPD). A total of 980 sputum supernatants were analyzed from 77 healthy never smokers, 494 former smokers (233 with COPD), and 396 active smokers (151 with COPD). Sputum nicotine, cotinine, and hydroxycotinine concentrations corresponded to self-reported smoking status and were strongly correlated to urine measures. A cutoff of ∼8-10 ng/mL of sputum cotinine distinguished never smokers from active smokers. Accounting for sample dilution during processing, active smokers had airway nicotine concentrations in the 70-850 ng/mL (∼0.5-5 µM) range, and concentrations remained elevated even in current smokers who had not smoked within 24 h. This study demonstrates that airway nicotine and its metabolites are readily measured in sputum supernatants and can serve as biological markers of smoke exposure. In current smokers, nicotine is present at physiologically relevant concentrations for prolonged periods, supporting a contribution to cigarette-induced airway disease.

Clinical Implications of Low Absolute Blood Eosinophil Count in the SPIROMICS COPD Cohort.

BACKGROUND: The Global Initiative for Chronic Obstructive Lung Disease (GOLD) considers blood eosinophil counts < 100 cells/µL (BEC≤100) in people with COPD to predict poor inhaled corticosteroid (ICS) responsiveness. However, the BEC≤100 phenotype is inadequately characterized, especially in advanced COPD. RESEARCH QUESTION: Are there differences between GOLD group D patients with high BEC and those with low BEC regarding baseline characteristics and longitudinal outcomes? STUDY DESIGN AND METHODS: We used multivariable mixed models and logistic regression to contrast clinical characteristics and outcomes of BEC≤100 vs BEC > 100 (BEC100+) in all subjects with COPD (n = 1,414) and GOLD group D subjects (n = 185) not receiving ICS. RESULTS: We identified n = 485 with BEC≤100 (n = 61 GOLD group D) and n = 929 people with BEC100+ (n = 124 GOLD group D). BEC≤100 status was stable at 6 weeks and approximately 52 weeks (intraclass correlations of 0.78 and 0.71, respectively). Compared with BEC100+, BEC≤100 comprised more women, with greater current smoking, and less frequent childhood asthma. Among all analyzed participants, the two BEC-defined subsets showed similar rates of lung function decline (mean slope, BEC≤100 vs BEC100+, -50 vs -39 mL/y; P = .140), exacerbations (0.40 vs 0.36/y; P = .098), subsequent ICS initiation (2.5% vs 4.4%; P = .071), and mortality (7.8% vs 8.4%; P = .715). However, in GOLD group D, people with BEC≤100 showed higher exacerbation rates within 365 days of enrollment (0.62 vs 0.33/y; P = .002) and total follow-up (1.16 vs 0.83/y; P = .014). They also had greater lung function decline (mean slope of -68 mL/y vs -23 mL/y; P = .036) and had greater emphysema at baseline (voxels < 950 Hounsfield units at total lung capacity of 7.46% vs 4.61%; P = .029). INTERPRETATION: In non-ICS-treated GOLD group D COPD, people with BEC≤100 had more baseline emphysema, prospective exacerbations, and lung function decline. Our analysis has identified a particularly vulnerable subpopulation of people with COPD, suggesting the need for studies focused specifically on their therapeutic treatment. CLINICAL TRIAL REGISTRATION: ClinicalTrials.gov; No.: NCT01969344; URL: www. CLINICALTRIALS: gov.

Genome-wide bidirectional CRISPR screens identify mucins as host factors modulating SARS-CoV-2 infection.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes a range of symptoms in infected individuals, from mild respiratory illness to acute respiratory distress syndrome. A systematic understanding of host factors influencing viral infection is critical to elucidate SARS-CoV-2-host interactions and the progression of Coronavirus disease 2019 (COVID-19). Here, we conducted genome-wide CRISPR knockout and activation screens in human lung epithelial cells with endogenous expression of the SARS-CoV-2 entry factors ACE2 and TMPRSS2. We uncovered proviral and antiviral factors across highly interconnected host pathways, including clathrin transport, inflammatory signaling, cell-cycle regulation, and transcriptional and epigenetic regulation. We further identified mucins, a family of high molecular weight glycoproteins, as a prominent viral restriction network that inhibits SARS-CoV-2 infection in vitro and in murine models. These mucins also inhibit infection of diverse respiratory viruses. This functional landscape of SARS-CoV-2 host factors provides a physiologically relevant starting point for new host-directed therapeutics and highlights airway mucins as a host defense mechanism.