Citations

Abstract RATIONALE. Significant sex-differences exist in asthma across the lifespan. Asthma prevalence and severity are greater in pre-pubescent boys compared to girls but greater in women after puberty. In women, asthma may improve after menopause when sex hormones wane. To address potential mechanisms underlying these changes, we compared the airway transcriptomics signatures of post-menopausal vs. premenopausal women. METHODS. We used bulk RNA seq data from bronchial epithelial cells from three different cohorts (SARP III, GEO-GSE85567, GEO-GSE201955) to perform a comprehensive analysis of differential gene expression and regulation. Data were combined into a meta-analysis, to compare pre-menopausal (n=114 cases with asthma, 48 controls) to post-menopausal women (n=52 cases with asthma, 13 controls). Differential gene expression between pre- and postmenopausal women was analyzed separately for asthma (Figure 1A) and control groups (not shown). Significant differentially expressed genes (DEGs) were identified by contrasting effect sizes (log fold-changes) between asthma and control groups using a Z-score interaction test. Accordingly, a scatter plot (Figure 1B) was created to represent whether differential gene expression between pre- and post-menopausal conditions is modified by disease status (asthma vs. control). The scatter plot also highlights that these findings are related to asthma per. se. and not to normal aging and menopausal transition. Pathway enrichment analysis was conducted using the differential expression data, followed by differential regulation analysis of the identified pathways. Statistical significance was assessed using a false discovery rate (FDR) threshold of 5%. RESULTS. Out of 15,468 genes that passed quality control, 19 genes were differentially expressed in bronchial epithelial cells from pre- and post-menopausal women with asthma (as compared to controls) (Figure 1A). Of these, 5 genes were down-regulated (ITLN1, POSTN, FCRL3, CCL22, PCSK6), and 14 genes were up-regulated (SPP1, S100A9, A2ML1, TNFRSF10C, CSF3R, OLFM1, C1QL2, NLRP3, HAL, PI3, RNASE1, TMEM178B, TNNI3, RNF175) when comparing pre- and post-menopausal women with asthma (Figure 1A). Of these 19 DEGs, 14 were differentially expressed in pre- as compared to post-menopausal women with asthma, but not in controls (Figure 1B). Comparing pre- to post-menopausal women with asthma, biologic pathways related to the innate and adaptive immune system activation and mitochondrial function (Electron transport Chain oxidative phosphorylation and ATP production) were both differentially regulated and expressed. CONCLUSIONS. This analysis highlights several DEGs and pathways that characterize asthma pathogenesis when comparing pre- and post-menopausal women, and emphasize the role of female sex hormones on the immune system response and mitochondrial function in asthma pathogenesis.

Abstract Rationale: SARS-CoV-2 infection may be associated with new or worsening respiratory symptoms that, if persistent, may contribute to “Long COVID.” Research is needed on the prevalence and risk determinants of Long COVID and related respiratory symptoms in adults with and without pre-existing obstructive lung disease or smoking history. Methods: C4R is a prospective, observational study that collected information on SARS-CoV-2 infection using standardized methods across 14 NIH-funded cohorts. The third wave of C4R questionnaires, administered June 2023 to January 2025, assessed self-reported Long COVID via the question, “Do you think that you are experiencing, or have ever experienced, what has been called “long COVID,” or symptoms related to COVID at least a month after your infection?” Questionnaires also assessed the presence (vs absence) of 22 symptoms that have been linked to Long COVID, including cough, shortness of breath, and post-exertional malaise. Pre-pandemic obstructive lung disease (COPD or asthma) and smoking status (ever/never) were harmonized from self-reported data collected at the most recent pre-pandemic cohort exam. Poisson regression was used to assess prevalence ratios (PRs) adjusting for age, sex, BMI, education, and pre-pandemic hypertension, diabetes, cardiovascular disease, and depressive symptoms. PRs for respiratory symptoms were tested according to self-reported Long COVID and history of non-hospitalized or hospitalized infection; effect modification by pre-pandemic lung disease and smoking was assessed by multiplicative interaction terms and in fully stratified models. Results: Among 23,434 questionnaire respondents, 1,484 (6.3%) self-reported Long COVID, including 246/1,262 (19.5%) with pre-pandemic lung disease and 664/4,018 (16.5%) with smoking history. The adjusted PR for Long COVID was 1.17 (95%CI:0.99-1.39;P=0.06) comparing participants with/without lung disease and 1.17 (95%CI:1.04-1.32;P=0.01) comparing participants with/without smoking history. Respiratory symptoms were more prevalent in those with self-reported Long COVID, as well as in those with a history of hospitalized or non-hospitalized infection (Table). Associations were generally stronger in participants without pre-pandemic lung disease (p-interaction<0.05 for all comparisons). Respiratory symptoms were also more prevalent with a history of repeated, and pre-Omicron infections, and were less prevalent among participants vaccinated prior to their first infection; findings were similar even after excluding participants with self-reported Long COVID. Conclusions: Long COVID is common, particularly in adults with pre-existing COPD or asthma or smoking history. Long COVID and history of SARS-CoV-2 infection, including non-hospitalized infections, was associated with greater respiratory symptom burden. Whether SARS-CoV-2 infection or Long COVID are associated with development or progression of obstructive lung diseases warrants future studies.

Abstract Rationale: Adverse social determinants of health (SDOH) are common in COPD and associated with increased mortality. However, cause-specific mortality is understudied. We evaluated associations of neighborhood-level SDOH with all-cause and cause-specific mortality among individuals with or at-risk for COPD. Methods: SPIROMICS is an ongoing, multi-center, longitudinal observational study of individuals across the US with current or former personal smoking histories (≥20 pack-years), with or without spirometrically-defined COPD. We established geocoded neighborhood-level SDOH at baseline using 2010 census tract data, including specifically, area deprivation index (ADI) (range 0-100, 100 = most disadvantaged), poverty level (percentage of households below the poverty line), poor food access (>1/3 of residents with no access to fresh food stores within 1 mile in urban areas and 10 miles in rural areas), and food desert (low income plus poor food access). We ascertained all-cause and cause-specific mortality through 12/31/2021 using the National Death Index. Using Cox proportional hazards models, we evaluated the risk of all-cause mortality for each SDOH. We utilized competing-risk regression models to evaluate the risk of cardiovascular, respiratory, COPD (a subset of respiratory), and neoplastic causes of death individually with each SDOH. Models were adjusted for demographics, smoking status, pack-years, BMI categories, coronary heart disease, FEV1% predicted and study site. Effect modification by lung function was assessed using multiplicative interactions between SDOH and FEV1% predicted. Results: Among 2977 participants, mean ADI was 40.1±27.8. Within a participant's neighborhood, on average 10.4±10.4% of households were below the poverty line, 933 (31.3%) had poor food access, and 220 (7.4%) resided in a food desert. During a median follow-up of 7.8 years (IQR 5.6-9.2), 744 (25.0%) participants died, most commonly from respiratory causes including COPD. Neither ADI nor neighborhood poverty level were associated with all-cause mortality or any specific cause of death (Table). Despite lack of associations between food access and mortality (Table), individuals residing in a food desert had elevated risks of all-cause (HR: 1.36 [95%CI: 1.03-1.79]) and cardiovascular (HR: 2.62 [95%CI: 1.5-4.35]) deaths. There was no significant effect modification by lung function. Conclusions: In a large cohort of individuals with personal smoking histories, with or without COPD, living in a food desert was associated with increased all-cause and cardiovascular mortality. However, neighborhood poverty level and food access, the two components of food desert, were not independently associated with mortality, suggesting the association of multimodal SDOH variables with mortality in this population.

Abstract Rationale: MUC5B and MUC5AC are neighboring genes on chromosome 11p15 that encode the two predominant mucins secreted in the respiratory airways. Increased airway mucin concentrations and an imbalance between these two mucins drive COPD symptoms and the development of airflow obstruction or COPD, but the epigenomic regulation of this loci has not been characterized in the airways. Methods: We isolated DNA from bronchial airway brushings of 72 participants with >20 pack-years smoking history from the Subpopulations and Intermediate Outcomes Measures in COPD Study (SPIROMICS). 36 participants had GOLD-based moderate-to-severe COPD (FEV1/FVC<0.70, FEV1<80%) and 36 were without or with mild COPD. We used the Illumina Infinium Methylation EPIC BeadChip Array to assess proportion of DNA methylation (β) spanning 689,937 CpG sites. We evaluated proportional DNA methylation (β) at 96 CpG loci within a 103kB region using matrixQTL for cis-expression quantitative trait methylation (eQTM: 10kB flanks relative to MUC5AC/MUC5B) integrating bulk RNA-seq data and mQTL models integrating gene sequencing data to identify DNA variation associated with proportion eQTM methylation. Models included age, sex, BMI, clinical site, current smoking vs past smoking, total pack-years, three ancestral components of ancestry, chipset, and position. Results: We identified a top eQTM, cg00524486, marginally associated with MUC5B expression (coefficient estimate=8.75, punadjusted=0.0014, pFDR=0.065, Figure 1B). This eQTM was associated with CT scan parametric response mapping-based functional small airways disease (PRM-fSAD, beta coefficient[β]=3.64, p=0.011). We tested for DNA variation associated with methylation of this top eQTM which identified two methylation quantitative trait loci (mQTL: rs885454 [β=0.044, minor allele frequency[MAF]=0.09, pFDR=0.01], rs868902 [β=0.020, MAF=0.43, pFDR=0.08, Figure 1A) associated with cg00524486 methylation, of which the more frequent mQTL, rs868902, was marginally associated with PRM-fSAD (p=0.07) in the larger SPIROMICS cohort with MUC5B DNA sequencing data. Conclusions: This preliminary multi-omic study of the lower airway epithelium demonstrates that the expression of MUC5B is epigenetically regulated by neighboring CpG probe methylation (eQTM) represented by DNA variation (mQTL). These studies require further confirmation in a larger dataset of airway samples and support a role for epigenetic variation as a biologic contributor and biomarker of early COPD pathogenesis.

Abstract Chronic obstructive pulmonary disease (COPD) and pulmonary emphysema were the third-leading cause of death globally in 2019. Besides alpha-1 antitrypsin deficiency, emphysema lacks disease-modifying medications due to its heterogeneity and there remains a need to further investigate emphysema subtypes to advance its diagnosis and treatment. Recent unsupervised machine learning of the texture and anatomical location of emphysema on computed tomography (CT) images of the lung defined six quantitative CT emphysema subtypes uniquely associated with specific gene variants, environmental exposures, symptoms, physiology and clinical events. We aim to identify associations of plasma proteins with each CT emphysema subtype cross-sectionally in a multi-ethnic general population sample. MESA is a prospective cohort study that recruited 6,814 participants from six communities who self-identified as White, Black, Hispanic, or Asian race/ethnicity, ages 45-84 years old and free of clinical cardiovascular disease in 2000-02. The MESA Lung Study performed CT scans among 3,137 participants in 2010-12, which were labeled with previously described (PMID: 37268414) CT emphysema subtypes. Proteomics was assessed using SomaScan version 4.1 assay (7,596 plasma proteins). Linear regression models were adjusted for age, height, weight, sex, race, educational attainment, smoking status, pack years, urine cotinine level, body-mass index (BMI), estimated glomerular filtration rate (eGFR), CT scanner manufacturer and other CT emphysema subtypes. Results were adjusted to reflect a False Discovery Rate (FDR) p-value of less than 5% using the Benjamini-Hochberg method. Among 2,343 participants with both CT and proteomic data (mean age 69.5 years, 52.3% women, 23.7% with COPD), we identified statistically significant associations of proteins with the vanishing lung subtype (42 proteins; top hit Natural cytotoxicity triggering receptor 3 ligand 1, q = 2.9[asterisk]10-9), restrictive combined pulmonary fibrosis/emphysema (CPFE) subtype (11 proteins; top hit Integrin alpha-2/b1, q=1.0[asterisk]10-5), the combined bronchitis-apical emphysema subtype (9 proteins; top hit Desmoglein-2, q=8.8[asterisk]10-13), and the diffuse emphysema subtype (4 proteins; top hit Pulmonary surfactant-associated protein D, q=5.7[asterisk]10-3). No significant proteins were associated with the obstructive CPFE or senile subtypes after FDR correction. Most proteins were uniquely associated with one CT emphysema subtype, while a few proteins were associated with more than one subtype. This study identified specific proteins for different CT emphysema subtypes, many of which have been related less specifically to emphysema, including pulmonary surfactant production, anti-inflammatory pathways, and mitochondrial deregulation, amongst others. Validation of these proteomic associations can provide insight into further pathway assessment, clinical biomarker discovery, and targeted prevention of emphysema in the general population.

Abstract Rationale Chronic obstructive pulmonary disease (COPD) manifests later in life and is a leading cause of death worldwide. While smoking-related lung pathology is a long-recognized COPD etiotype, emerging evidence suggests that early-life growth conditions also play an important role. A recently-developed index of early-life growth conditions – the difference between attained height and genotype-predicted height (height-GaP) – has been validated against prospectevely ascertained adverse early-life growth conditions. This study used height-GaP to evaluate the association of early-life growth conditions with COPD phenotype and prognosis later in life. Methods Among current and former smokers with spirometry-defined COPD in the SubPopulations and Intermediate Outcomes Measures in COPD Study (SPIROMICS), we computed height-GaP as the difference between attained height and genotype-predicted height. Regression models were fit for respiratory mortality (defined by primary diagnostic code), dyspnea prevalence (Medical Research Council dyspnea rating >1), chronic bronchitis prevalence, airflow obstruction severity as quantified by FEV1Q, and functional capacity (six-minute walk distance), adjusting for age, age2, sex, race-ethnicity, genotype-predicted height, cigarette smoking status and pack-years. Sensitivity analyses accounted for competing risks of death (Fine-Gray subdistribution hazard model) and accounted for later-life height loss (restricting sample to participants under 50 years). Results Among 1,258 participants with COPD (mean±SD age: 65.4±8.0 years, 42.4% female, mean±SD pack-years: 53.6±28.2, mean±SD FEV1%: 63.3±23.9%), there were 202 respiratory deaths over a median 7.5 years of follow-up (Figure). Larger height-GaP deficit was associated with higher respiratory mortality (adjusted hazard ratio [aHR]: 1.30 [95%CI: 1.12 to 1.50] per 1-SD increment in height-GaP deficit). This association remained significant after accounting for competing events (p=0.001) and when restricted to participants under 50 years (p=0.032). Larger height-GaP deficit was associated with higher dyspnea prevalence (adjusted prevalence ratio: 1.13 [95%CI: 1.02 to 1.25] per 1-SD increment in height-GaP deficit), lower FEV1Q (adjusted mean difference: -0.37 [95%CI: -0.45 to -0.29] per 1-SD height-GaP deficit), and shorter 6-minute walk (adjusted mean difference: -16.6 m [95%CI: -23.1 m to -10.0 m] per 1-SD increment in height-GaP deficit), but not chronic bronchitis or CAT score ≥10 (p>0.05 for each). Conclusions Among current and former smokers with COPD, a quantitative index of early-life conditions was associated with higher respiratory mortality, in addition to dyspnea, airflow obstruction severity and reduced functional capacity. These results suggest that early-life growth conditions may influence COPD phenotype and prognosis later in life.

Abstract Background:Chronic obstructive pulmonary disease (COPD) is a major global health burden, with genome-wide association studies explaining only a fraction of its heritability. Recent advancements, including large-scale, multi-ancestry whole-genome sequencing and a focus on severe, early-onset cases, offer promising avenues to identify additional COPD-related variants. Furthermore, recent studies show that adjusting for polygenic risk scores (PRS) enhances statistical power to detect rare variants, advancing efforts to understand COPD's genetic basis. Methods:We performed a whole genome sequencing association study of severe, early-onset COPD (FEV1 < 50%, age ≤ 55), FEV1 and FEV1/FVC in eight population-based studies (ARIC, CARDIA, CFS, CHS, FHS, HCHS/SOL, JHS, MESA) and seven COPD-enriched studies (COPDGene, ECLIPSE, Boston Early-Onset COPD Study, LTRC, TCGS (Poland), ICGN, SPIROMICS). We performed single-variant analyses using regenie software, applying a genome-wide significance threshold of P < 5.0 x 10-8. For rare variant set (MAF <1%) analysis, we annotated variants with loss of function or deleterious missense using WGS Annotator and prediction tools LOFTEE, SIFT, and PolyPhen2, then analyzed these variants with SKATO-ACAT method in regenie, adjusting for covariates including sex, sequencing center, and principal components of ancestry. We additionally evaluated adjusting the analyses for PRSs for lung function traits. Results: We analyzed 53,028 participants including 34,751 European, 12,508 African, 1,088 Latino American, and 674 East Asian ancestry. Across cohorts, we analyzed 5,063 severe (FEV1 < 50%) COPD cases, among whom 1,103 were early-onset (age ≤ 55). Excluding variants within 1 Mbp of known variants, single-variant analysis identified 68 novel associations (MAC ≥ 30), with 8 common associations (MAF ≥ 1%) in 6 genomic regions. Our most significant novel association for severe early-onset COPD was rs74297390 (intergenic to POM121L12, MAF = 9.6%, Beta = 0.41, P = 3.5 x 10-8). For FEV1 and FEV1/FVC ratio, variants near DKK2 and AKR1E2 were significantly associated and were replicated in an independent dataset (Shrine et al., 2019). DKK2 is an antagonist in the Wnt pathway. In a preliminary variant set analysis, CHRM3-AS1 was associated with severe early-onset COPD, GPAT with FEV1, and COL8A1I and HMCN1 with FEV1/FVC. With PRS adjustment, ST6GALNAC5 was associated with FEV1/FVC. Conclusions: In a whole genome sequencing study of lung function and severe early-onset COPD, we identified new loci for COPD-related phenotypes, including two novel loci that replicated in an independent dataset. Ongoing work involves variant set-based analyses with additional covariate adjustments, as well as PRS adjustments for binary traits.

Abstract RATIONALE The best described endotype of COPD is alpha-1 antitrypsin (AAT) deficiency, which is due to a genetic abnormality in the SERPINA1 gene. Homozygotes for the Z allele (PI[asterisk]ZZ individuals) carry the genotype most commonly associated with severe AAT deficiency, but a highly prevalent endotype is the heterozygous state (PI[asterisk]MZ individuals). Here we test the hypothesis that individuals with PI[asterisk]MZ genotype are at higher risk for COPD exacerbation and healthcare utilization compared to those with no pathogenic variants in the SERPINA1 gene (PI[asterisk]MM individuals). METHODS Real-world data on 4,456 individuals with available AAT genotyping were extracted from the Cleveland Clinic Health System (CCHS) electronic health records (EHR) and analyzed to compare the risk of moderate COPD exacerbation (defined as a short-term steroid prescription), all-cause emergent care (defined as an express care visit, urgent care visit, or emergency department visit), and all-cause hospitalization between individuals with PI[asterisk]MZ vs. PI[asterisk]MM genotypes. Cox proportional hazards models were used to adjust for age, sex, race, BMI, smoking status, comorbidity count, liver disease diagnosis, and median income by zip code.To assess the effect of AAT levels on outcomes, we categorized participants into three groups based on AAT level and used 90 mg/dL as a threshold given it is the lower limit of the normal reference range and genotype: MM individuals with an AAT level >90 mg/dL, MZ >90 mg/dL and MZ <90 mg/dL. To assess the effect of PI[asterisk]MZ genotype among individuals with pre-existing lung disease, we performed a subgroup analysis of individuals with at least a six-month inhaler prescription. RESULTS In this cohort, 308 individuals had a PI[asterisk]MZ genotype. PI[asterisk]MZ was associated with increased risk for moderate COPD exacerbations (adjusted HR [95% CI]: 1.66 [1.27, 2.17]) and hospitalizations (adjusted HR [95% CI]: 1.44 [1.19, 1.75]) compared to PI[asterisk]MM. The risk of hospitalization was higher among PI[asterisk]MZ individuals with AAT levels <90 mg/dL (adjusted HR [95% CI]: 1.59 [1.14, 2.23]) but not in those with AAT levels >90 mg/dL, as compared to PI[asterisk]MM. Sensitivity analysis in those treated with inhaler therapy showed overall similar results. CONCLUSIONS Increased hospitalizations and moderate COPD exacerbations among PI[asterisk]MZ individuals compared to PI[asterisk]MM individuals merit future prospective studies to better characterize the longitudinal course and healthcare utilization among those with an PI[asterisk]MZ genotype. Given the high prevalence of these carriers globally, PI[asterisk]MZ represents a COPD phenotype that invites additional investigation to identify biomarkers of worse disease and treatable traits.

Abstract Rationale: Individuals with chronic obstructive pulmonary disease (COPD) have different clinical trajectories, some with minimal symptom burden, others experiencing recurrent exacerbations, the most severe of which require hospitalization. Existing diagnostics, including spirometry and symptom assessments, lack the precision to identify high-risk patients who may benefit from timely preventive therapies. Nuclear cell-free DNA (ncfDNA) is a well-established marker of systemic tissue injury, implicated in amplifying inflammation by triggering type 1 interferon release by activating DNA-sensing pathways. Because COPD is accompanied by inflammation and tissue injury systemically, we hypothesized that ncfDNA is associated with increased risk of severe exacerbations requiring hospitalization. Methods: Design: We utilized baseline plasma samples randomly selected from 672 individuals with COPD from the Subpopulations and Intermediate Outcomes Measures in COPD Study (SPIROMICS) cohort, which prospectively collected exacerbation numbers and type. The primary endpoint was exacerbation requiring hospitalization. Secondary endpoints included total number of exacerbations (defined as use of oral steroids, antibiotics, or healthcare utilization) and severe exacerbations (defined as requiring emergency room care or hospitalization). Measure: We measured plasma ncfDNA using digital droplet PCR. Analysis: We assessed the association of ncfDNA with exacerbation outcomes using negative binomial regression or Cox regression, adjusting for age, percent post-bronchodilator forced expiratory volume in one second (FEV1), supplemental oxygen use, and smoking status. Results: In total, 57% (n=383) of COPD participants experienced exacerbations, 29% (n=197) at least one severe, and 22% (n=145) requiring hospitalization, the latter at a median of 622 days from enrollment (IQR=287-997). Higher ncfDNA levels at enrollment correlated with measures of COPD severity, including baseline FEV1 (beta = -0.12, 95% CI -0.18 to -0.049, p <0.001). Higher ncfDNA associated with increased risk of COPD exacerbation requiring hospitalization (incidence rate ratio (IRR) 1.40, 95% CI 1.12-1.75, p =0.003) and severe exacerbations (IRR 1.36, 95% CI 1.11-1.67, p=0.002) and a trend towards total exacerbations (IRR 1.11, 95% CI 0.98-1.26, p= 0.09). The primary outcome, risk of time-to-first hospitalization, correlated significantly with higher ncfDNA with hazard ratio (HR) 1.23, 95% CI 1.01-1.49, p=0.037. Conclusion: These results suggest that ncfDNA may serve as a predictor for severe and hospitalization-associated COPD exacerbations. To better understand the relationship between ncfDNA and COPD pathophysiology, follow-up studies should explore longitudinal trends and cfDNA tissue sources, and their link to COPD exacerbations.

Abstract Rationale: We have previously shown reduced total IgG, IgG subclasses, and levels/function of pneumococcal-specific antibodies are associated with increased risk for exacerbations of COPD (ECOPD). The binding of IgG to its receptor, the Fc receptor for IgG (FcγR) on immune cells, has a critical role in determining antibody levels and function. In addition to our research demonstrating that IgG heavy chain gene (IGHG) variants are associated with ECOPD, single nucleotide variants (SNVs) in genes encoding FcγR (FCGR1-3, which encode receptor subtypes FcγRI, FcγRIIa/b/c, and FcγRIIIa/b) may also result in antibody deficiency and increased exacerbation risk. Due to technical challenges in accurately identifying SNVs in FCGR, this region has not been well-characterized. The limited data that exist are from people of European descent, and prior studies have not included African American participants. We therefore sought to investigate the hypothesis that distinct SNVs in genes encoding FcγR are associated with disease severity and history of exacerbations in African American people with COPD. Methods: We developed a novel pipeline to identify SNVs in the coding (CDS) regions of chromosome 1 using short-read whole genome sequence data. The pipeline was benchmarked using paired long-read samples from Genome in a Bottle. SNVs were called within CDS regions and analyzed in a meta-analysis of self-reported African American participants in SPIROMICS and COPDGene (n=3063). Genetic variants were tested for association with COPD severity (GOLD spirometry stage 1-2 vs 3-4) as well as COPD exacerbations, adjusting for age, sex and principal components controlling for population structure. Results: Eight SNVs were significantly associated with frequent ECOPD (≥2 in one year) in the meta-analysis of African American COPDGene and SPIROMICS participants (Table). Seven of these 8 SNVs were found in FCGR3; 4 in the activating FCGR3A subtype, and 3 in FCGR3B. Seven of 11 SNVs associated with COPD severity were in FCGR3 (4 in FCGR3A; 3 in FCGR3B). One SNV (rs371849) was associated with both frequent exacerbations (P=0.002) and COPD severity (P=0.044). Conclusions: We used a cutting-edge variant calling pipeline to identify multiple SNVs in FCGR, primarily involving FcγRIII receptor subtypes, associated with COPD severity and exacerbations in African American participants in SPIROMICS and COPDGene. FcγRIII variants are known to cause impaired antibody levels and function in people of European descent. The structural and functional importance of newly recognized variants in African American cohorts require further study, as they may represent immunogenetic risk factors for antibody deficiency and exacerbations.