Abstract
Introduction Mucus pathology plays a critical role in airway diseases like chronic bronchitis (CB) and chronic obstructive pulmonary disease (COPD). Up to 32% of community-living persons report clinical manifestations of mucus pathology (e.g., cough and sputum production). However, airway mucus pathology has not been systematically studied in community-living individuals. In this study, we will use an objective, reproducible assessment of mucus pathology on chest computed tomography (CT) scans from community-living individuals participating in the Coronary Artery Risk Development in Young Adults (CARDIA) and Framingham Heart Study (FHS) cohorts.
Methods and analysis We will determine the clinical relevance of CT-based mucus plugs and modifiable and genetic risk and protective factors associated with this process. We will evaluate the associations of mucus plugs with lung function, respiratory symptoms, and chronic bronchitis and examine whether 5-yr. persistent CT-based mucus plugs are associated with the decline in FEV1 and future COPD. Also, we will assess whether modifiable factors, including air pollution and marijuana smoking are associated with increased odds of CT-based mucus plugs and whether cardiorespiratory fitness is related in an opposing manner. Finally, we will determine genetic resilience/susceptibility to mucus pathology. We will use CT data from the FHS and CARDIA cohorts and genome-wide sequencing data from the TOPMed initiative to identify common and rare variants associated with CT-based mucus plugging.
Ethics and Dissemination The Mass General Brigham Institutional Review Board approved the study. Findings will be disseminated through peer-reviewed journals and at professional conferences.
Strengths and limitations of this study
Utilization of data from two well-characterized large community-based US cohorts.
Use of chest CT scans to identify and quantify mucus plugs, providing a more objective and reproducible measure of airway pathology.
Use of whole-genome sequencing to identify common and rare genetic variants associated with mucus pathology.
Only the inclusion of participants self-identified as non-Hispanic white and non-Hispanic black.
A limitation of retrospective study design using prospectively collected data.
Introduction and Rationale
Mucus pathology plays a critical role in airway diseases, such as chronic bronchitis (CB), asthma, and chronic obstructive pulmonary disease (COPD) [1–9]. COPD affects ∼29 million people and is the 4th leading cause of death in the US [10, 11]. Mucus pathology and plug formation are complex processes in response to multiple inhaled noxious particles, including cigarette smoke and air pollution [1, 12, 13]. The public health and clinical relevance of airway mucus pathology are better understood through its clinical manifestations, including chronic cough and phlegm and chronic mucus hypersecretion/CB [14–19]. The prevalence of chronic respiratory symptoms in primary care is up to ∼40% [20]. In a large study, 32% of ∼98,000 community-living individuals with normal lung function reported chronic respiratory symptoms. More importantly, those symptoms were indicators of an increased risk of future hospitalization and mortality [14].
We and others have also demonstrated that in smokers with COPD, mucus plugs occluding the lumen of medium-to-large-sized airways (i.e., 2-10-mm-lumen diameter), those visualizable on CT, are associated with impaired lung function, worse quality of life, and increased all-cause mortality [6, 7, 21]. We posit that CT assessment of airway mucus pathology applies to community-living individuals [22, 23].
Extensive research has been conducted on chronic respiratory disorders. While this has brought new insights into diseases’ clinical variability and biological basis, attention has been limited to patients with established and even advanced disorders. Suppose clinical care is to identify early risk factors and prevent disease. In that case, new work is needed to identify those at high risk of future diseases, those who may have poor lung health but do not yet meet existing diagnostic criteria for an established airway disease [17]. We have shown that smokers have persistent mucus plugs on CT over a 5-year follow-up [6], but whether this pattern represents a risk factor for airway disease in community-living adults is unknown.
In this study, we will also determine the role of modifiable factors, such as air pollution, marijuana, and fitness, on mucus pathology, which may inform specific preventive and treatment strategies for lung health (e.g., air purification). Exposure to air pollution and cigarette smoke induce pathologic mucus production in the airways; however, not all individuals exposed to the same set of factors develop comparable mucus dysfunction and its consequences, suggesting genetic susceptibility or resilience may have a role [23].
Our team has implicated TP53 gene variants in increased mucus synthesis and sustained the lifespan of metaplastic bronchial mucous cells. Further, smokers with the p53Arg are at increased risk of CB (OR=1.59; p=0.03) [24].
We propose that CT-based mucus plugs are an intermediate phenotype that will enable us to investigate the clinical implications of mucus pathology and assess its contribution to the risk of airway disease. To test our hypothesis, we will use data from two population-based cohort studies: The Coronary Artery Risk Development in Young Adults (CARDIA) study [25] and the Framingham Heart Study (FHS) [26, 27]. The study protocol presented here was funded by the National Heart Lung and Blood Institute (Grant R01-HL164824).
Specific Aims
Our overall hypothesis is that CT-based mucus plugs are an intermediate phenotype that will enable us to examine the clinical implications of mucus pathology and gauge its contribution to airway disease risk in population-based studies. The study has three specific aims:
Determine whether CT-based mucus plugs are associated with reduced FEV1, respiratory symptoms, and CB in community-living adults.
Determine whether modifiable exposures (long-term air pollution, marijuana use, and fitness) are associated with CT-based mucus plugs.
Assess the genetic resilience/susceptibility to mucus pathology by examining whether common and rare variants are associated with CT-based mucus plugging in community-living adults.
We will use existing imaging, clinical, environmental, and genetics data already collected from the CARDIA and FHS cohort, thus not exerting additional burden to the study participants and avoiding further radiation exposure.
Methods and Analysis
Study design
The present will use prospectively collected data from the CARDIA and FHS cohorts. CARDIA is a United States (US) National Institutes of Health (NIH)-funded population-based cohort of 5,115 young adults recruited from 1985 to 1986 to explore the evolution of cardiovascular disease risk factors [25]. Participants were examined at four sites across the US during enrollment. In-person examinations occurred at baseline (year [Y] 0), Y2, Y5, Y7, Y10, Y15, Y20, Y25, and Y30. Examinations collected demographic, anthropometric, clinical (e.g., symptoms and tobacco and marihuana smoking history), fitness, chest CT, and spirometry data in a standardized manner. CARDIA also took and kept blood samples. Participants’ blood samples were sent to the NIH Trans-Omics for Precision Medicine (TopMed) program to perform whole-genome sequencing. We will use chest CT images from the Y20 (n = 3,122) and Y25 (n = 3,498) examinations to score airway mucus plugs.
The FHS started in 1948 as a prospective investigation of cardiovascular disease, and we will use participants from the FHS Multidetector Computed Tomography 2 (MDCT2) sub-study (n=2750). Participants of the MDCT2 sub-study are from the FHS Offspring Cohort (1971) and the FHS Third Generation (2002) cohort [27]. The FHS Third Generation participants are the children of the Offsprings. FHS also collected demographic, anthropometric, clinical (e.g., symptoms and tobacco smoking history), air pollution, spirometry, and chest CT data in a standardized manner. FHS participants’ blood samples were also sent to the NIH TopMed program to perform whole-genome sequencing. We will use chest CT images from the MDCT2 sub-study to identify and score airway mucus plugs. Characteristics of the Y20 CARDIA and FHS participants are depicted in Table 1.
Patient and Public Involvement Statement
Patients and/or the public were not involved in the design, conduct, reporting or dissemination plans of this research.
Eligibility
In the CARDIA cohort, eligibility was limited to young adults aged 18-30 years in 1985 to 1986 with no evidence of long-term diseases or disabilities. In the FHS cohort, eligibility was limited to participants aged 20 years or older and had at least one parent who was a member of the Offspring Cohort [27]. All individuals known to be eligible for the inclusion in the Third Generation Cohort were sent invitation letters along with response letters, and high interest was indicated by return of response cards, which led to more participation than anticipated.
Participant Background
At Y20, the mean age was 45 years; 38% were smokers, and 9% - 19% reported cough, phlegm, or dyspnea. CARDIA included a balanced number of male and female and Black and White participants, whereas FHS had well-balanced sexes (M/F) but only included White participants. In the FHS-MDCT2 cohort, the mean age was 59 years; 51% were former and current smokers, and 7-11% reported respiratory symptoms. Participants of both population-based cohort studies include varying cardiorespiratory fitness, air pollution exposure, and tobacco and marijuana smoking histories.
Study Data Collection
The CARDIA and FHS studies used standardized questionnaires and procedures to collect imaging, clinical (e.g., marijuana history), environmental, fitness, and genetic data. Also, the methodology to assess mucus plugs on CT scans for this study is detailed below. The Institutional Review Board of Mass General Brigham approved this study. All participants provided written consent.
Chest CT scans
All participants underwent chest CT scanning. CARDIA patients underwent a noncontrast electrocardiography-gated chest CT scan reconstructed at 2.5 to 3 mm slice thickness [28]. FHS-MDCT2 sub-study patients underwent lung volumetric CT scans in the supine position, reconstructed by a sharp lung algorithm with a 0.625 mm thickness and a 50 cm field of view [29]. Although there are differences in reconstructions between cohorts, the quality of CT scans reconstruction parameters is adequate for this study to survey the airways.
Defining Mucus Plugs
A mucus plug is defined as an opacity that completely occludes the lumen of an airway on sequential slices, regardless of the generation and orientation of the bronchus (Figure 1) [3,4]. The opacity can be round or cylindrical in shape and found in one or many airway branches. Plugging is identified in medium-to-large-sized airways on CT.
Mucus Plug Score
The scoring method was developed based on the number of pulmonary segments with mucus plugs regardless of the orientation of the bronchus (Figure 2). For instance, a mucus plug identified in the lower right lobe’s posterior and lateral basal segments yields a score of two. The score will range from 0 to 18, where 0 means no plugging and higher numbers indicate a greater extent of plugs.
Reading Method
The readers who will score mucus plugs comprise thoracic radiologists, a pulmonologist, physicians, and graduate students. Participating readers undergo training, which consisted of identifying a variety of mucus plugs on a set of typical CT images, scoring a training set of CT scans to learn the scoring system, and scoring a sample set of 20 CT scans to verify the inter-reader agreement. The concordance correlation for inter-reader agreement of mucus plug score previously reported was 0.84 (95% Cl 0.81-0.88) [6]. The Applied Chest Imaging Laboratory also has developed an automated method to detect mucus plugs, providing an alternate way to assess airway mucus pathology on CT images.
Spirometry
Spirometry protocols were designed following the American Thoracic Society (ATS)/European Respiratory Society (ERS) guidelines. While CARDIA protocols did not include the use of a bronchodilator, the FHS-Offspring cohort (i.e., at examination 9) performed postbronchodilator measurements in those with airflow limitation. Given the epidemiological nature of this study, we will use prebronchodilator measures of lung function. Airflow limitation is defined as prebronchodilator forced expiratory volume in 1 second (FEV1)/forced vital capacity (FVC) ratio < 0.70 or less than the lower limit of normal [30].
Symptoms
Chronic respiratory symptoms include asthma, airway wall thickening, dyspnea, phlegm, emphysema, and bronchiectasis. CT-based mucus plugging is associated with dyspnea, phlegm, cough, and CB.
Environmental Exposures + Air pollution
Information on FHS participants’ residential address, PM2.5 levels, and proximity to a major road in relation to air quality was determined with a geospatial model. Models contained daily air pollutants, weather, and chemical exposures that integrates satellite-based aerosol measurements, Goddard Earth Observing System atmospheric, weather data, and local land data. Machine learning techniques were employed to refine the model based on ground-level measured exposure data to estimate participant daily exposures at high resolution with excellent out-of-sample cross-validation performance. Proximity to major roads was calculated via the U.S. Census Feature Class A1, A2, and A3 data [31–33].
Marijuana Exposure
Information on CARDIA participants’ marijuana smoking habits was collected from Y0 to Y25. Lifetime marijuana use was defined as the number of joints smoked from Y0 onwards, including Y20 and Y25 [34].
Cardiorespiratory Fitness
Participants had a symptom-limited, graded treadmill exercise test to measure cardiorespiratory fitness during the Y20 CARDIA exam. The modified Balke protocol consists of nine 2-minute stages of increasing difficulty in both speed and grade to a maximum of 5.6 mph at 25% grade (19.0 metabolic equivalent task). The test duration and maximal metabolic equivalent task were recorded [28, 35].
Genetics
NIH TOPMed whole-genome sequencing will be used on both cohorts to identify common and rare variants associated with CT mucus plugging [36]. The NIH TOPMed Study performed alignment and variant calling using a standardized pipeline detailed at nhlbiwgs.org, yielding high-quality whole genome sequencing variant calls.
Hypotheses
There are three primary hypotheses:
♦ Hypothesis 1: CT-based mucus plugging is associated cross-sectionally with lower FEV1, higher burden of respiratory symptoms, and CB; persistent mucus plugs over time are associated with a decline in FEV1 and future development of COPD.
♦ Hypothesis 2: Air pollution and lifetime marijuana use are associated with increased odds of CT-based mucus plugging, whereas cardiorespiratory fitness is protective against airway mucus plug formation.
♦ Hypothesis 3: Common and rare variants are associated with CT-based mucus plugging.
Analysis Plan
To evaluate Hypothesis 1, we will conduct separate analyses for FHS and CARDIA. In CARDIA, linear and logistic multivariable models will be used to assess the cross-sectional relationships between CT-based mucus plugging and FEV1 % predicted (continuous), cough, phlegm, dyspnea, and bronchitis episodes (all binary). We will consider the following covariates: age, sex, race, center, education level, BMI, ever-smoking status, pack-year, and history of asthma. In FHS, linear mixed-effect models will be used for FEV1 and generalized mixed-effect models for binary outcomes: cough, phlegm, dyspnea, and CB. In these models, fixed effects, except for race, will be like those used in CARDIA. Since FHS includes participants with familial relatedness (parents and children), we will include a random effect to account for familial aggregation. Finally, we will conduct sex-stratified analyses for these outcomes, as preliminary studies suggest women have higher odds of mucus plugging.
We will use the CARDIA cohort to assess the association between CT-based mucus plugs and the outcome change in FEV1 over time, using Y20 and Y25 CT data. For these longitudinal analyses, we will classify changes in mucus plugging into two groups: Persistent Positive (mucus plugs at both Y20 and Y25) + Newly Formed (mucus plugs at Y25 but not at Y20) vs. Persistent Negative (no mucus plugs at both Y20 and Y25) + Resolved (mucus plugs at Y20 but not at Y25). The two latter groups combined are the referent group. We will use linear mixed models and generalized mixed models to assess whether that binary variable is associated with the outcomes decline in FEV1 and COPD at the Y30 exam. Decline in FEV1 will be defined as Y30 minus Y20 divided by 10, yielding changes per year. Unfortunately, CARDIA has no spirometry data at Y25. COPD is defined as a pre-bronchodilator FEV1/FVC ratio <0.7, an acceptable approach for epidemiologic studies [37, 38]. The cutoff of 0.7 has been validated in population-based studies [39]. Models for FEV1 decline will include age, sex, race, center, education level, BMI, ever-smoking status, pack-year, baseline FEV1, and history of asthma. Time-variant variables include pack-year, smoking status, and BMI. Also, we will perform analyses to explore sex and race differences in the associations of interest, drawing from the sex and racial-balanced number of CARDIA participants.
Power calculation and sample size for Hypothesis 1
Detectable effect sizes are calculated for FEV1 decline and future COPD to reach at least 80% power, assuming a sample size of 2,600 participants with CARDIA Y20 and Y25 CT data and an alpha value of 0.05. If the prevalence of persistent mucus plugging is 10%, we detect a 53mL FEV1 decline and a 1.87 odds ratio of future COPD. With a prevalence of persistent mucus plugging of 20%, we detect a 40 mL FEV1 decline and an odds ratio of 1.61 for future COPD.
To evaluate Hypothesis 2, we will perform separate analyses in CARDIA and FHS. For CARDIA, we will use a logistic model to assess the associated between lifetime marijuana smoking and the outcome CT-based mucus plugging (i.e., binary variable defined as mucus plug score 0 vs. >0). Based on prior studies, marijuana will be categorized to indicate high and low consumption. We will consider adjustments for age, race, education attained, BMI, FEV1, presence of emphysema on CT, and asthma history. Also, we will build separate models to examine the relationships between tobacco smoking (ever vs. never, pack-years) and mucus plugging and compared them. Finally, we will explore whether tobacco and marijuana smoking have independent or additive effects on mucus plug formation. In separate logistic models, we will test the association between fitness (i.e., treadmill duration) and mucus plugging. Further, we think that fitness may counter the effect of tobacco smoking on mucus formation. Thus, in secondary analyses, we will test whether fitness modifies (i.e., mitigates) the impact of cigarette smoking status and pack-years on CT-based mucus plugging using the appropriate interactions terms. For FHS, a generalized mixed model will test the association between air pollution and distance to roadways (exposures) with CT-based mucus plugging as the outcome. Based on previous studies, we will consider using two exposures: 5-year PM2.5 average concentration (2004-2008) before CT scanning (2008-2011) and distance to major roads as categories. We will consider the following covariates: age, sex, education level, BMI, smoking (ever vs. never), pack-years, the value of owner-occupied housing, population density in the census track, FEV1, and CT measures of emphysema and airway-wall thickening. A random effect will be used to account for FHS familial aggregation.
Power calculation and sample size for Hypothesis 2
Detectable effect sizes are also computed on the risk of mucus plugging (binary outcome) for fitness and PM2.5 and reach at least 80% power, assuming an alpha level of 0.05 and sample sizes of 2,677 (fitness, CARDIA) and 2,340 (PM2.5, FHS). CARDIA and FHS sample sizes were smaller due to missing fitness and PM2.5 data. In a prevalence of CT-based mucus plugs of 5%, there is a 1.29 relative risk of mucus plugging per 1 SD decrease in fitness and a 1.31 relative risk per 1 SD increase in PM2.5 mcg/m3. In a prevalence of CT-based mucus plugs of 20%, there is a 1.15 relative risk of mucus plugging per 1 SD decrease in fitness and a 1.16 relative risk per 1 SD increase in PM2.5 mcg/m3.
To evaluate Hypothesis 3, we will use NIH TOPMed data for common and rare variants. For common variants, we will examine variants with minor allele counts of above 30 in FHS and CARDIA. We will perform a combined analysis, leveraging the single dataset available for all TOPMed [40]. The main analysis will include all participants and adjust for age, sex, current smoking, pack-years, sequencing center/study, race, familial aggregation, and principal ancestry components in a linear mixed model implemented in SAIGE [41].To balance maximum power with the desire for replication, we will apply a tiered strategy previously used and consider ‘Tier 1’ signals those that reach pre-determined levels (see below) in FHS and replicate in CARDIA, and ‘Tier 2’ signals those that reach significance in the analysis of FHS and CARDIA, and that have nominally significant P values in both studies.
Given the importance of cigarette smoking and potential air pollution in mucus plug formation, we will also perform secondary analysis using common variants including both the main effect and a gene by environment interaction using a two-degree of freedom joint test that improves power in the setting of interaction [42]. Due to the differences in prevalence among races, we will also examine results by ancestry and perform admixture analysis in Black participants in CARDIA using LAMP-LD, as we have previously performed [43].
For rare variants, aggregating variants can increase power. We will use the recently developed ACAT-O (aggregated Cauchy association test, omnibus), which combines set-based tests including burden, SKAT, and ACAT-V, improves power in the setting of sparsity, and is superior to earlier methods such as SKAT-O [44]. We will test aggregating by gene or by sliding window region and also examine the effects of limiting our analysis to coding variants (nonsynonymous, stop, or splice) or those with putative regulatory function (e.g., based on overlap with open chromatin in lung-related cells) [45].
We also anticipate examining a specific set of ∼100 genes known to be important for mucus dysfunction (e.g., TP53) and muco-obstructive diseases such as COPD and bronchiectasis (e.g., FAM13A, MUC5B, MUC5AC, CFTR, and DNAH5). This secondary analysis decreases the multiple comparison penalty from genome-wide approaches while enabling a more detailed examination of genes with a higher likelihood of association. As a secondary analysis, we will perform a genome-wide association.
We will analyze heritability, genetic correlation, functional enrichment, and potential causal genes to analyze and interpret our genetic association results. We will use GCTA-SC and GCTA-LDMS in our European sample to estimate heritability and genetic correlation between this CT phenotype and lung function [46]. To identify potential causal genes, we will rely on multiple approaches, including expression quantitative trait loci (eQTL) in lung and blood, examination of DNAase hypersensitivity sites, chromosome conformation (Hi-C), and gene set and pathway methods (DEPICT), as previously performed [47].
Power calculation and sample size for Hypothesis 3
When assessing genetic variance to CT-based airway mucus plugs, in a total sample size of 5,427 (FHS + Y20 CARDIA) and assuming a prevalence of mucus plugging of 10%, we have ∼ 80% power to identify a variant at a P < 1 × 10-4 with an allele frequency of 40% and an odds ratio of ∼1.3. In addition, we have > 80% power at alpha = 5 × 10-8 in our genome-wide analysis to identify a variant that increases risk by 1.45.
Ethics and Dissemination
The present study was approved by the Mass General Brigham Institutional Review Board (2022P000723). Participants were recruited voluntarily for these studies. Written informed consent was obtained for CT scanning, clinical data, blood sample collection, and genomic analyses. Participants were under no obligation to participate. A potential risk is the risk of participant de-identification or the release of personal health information. We will use genotype and phenotype files, which are only identified by study ID numbers; the investigators are not members of the Data Coordinating Center of these studies and will therefore not have access to any information that is not already de-identified. To protect against inadvertent disclosure of identity, we will make all study data available to study investigators using the de-identified study identifications assigned by the CARDIA and FHS studies. The links to identifiers are not available to the investigators. All study data is stored on a secured, password-protected computer system at Brigham and Women’s Hospital (BWH) and CARDIA and FHS Data Coordinating Centers.
We plan to disseminate the findings of this study through peer-reviewed journals and professional scientific conferences, such as the American Thoracic Society’s annual Conference.
Data Availability
All data produced in the present study are available upon reasonable request to the authors.
Footnotes
Author Contributions
Guarantors of the integrity of the study, A.A.D.; study concepts/study design or data analysis/interpretation, A.A.D, R.E., M.C., R.K., M.R.; obtained funding, A.A.D., M.C.; data acquisition, M.R., G.O.,N.A., A.Y., M.H., M.A., R.E., A.A.D.; manuscript drafting or manuscript revision for important intellectual content, all authors; drafted the article and critically revised it for intellectual content, A.A.D., M.A.; approval of the final version of the submitted manuscript, all authors; agree to ensure any questions related to the work are appropriately resolved, all authors.
Funding
This work was supported by the National Heart Lung and Blood Institute (Grant R01-HL164824). The Coronary Artery Risk Development in Young Adults Study (CARDIA) is supported by contracts 75N92023D00002, 75N92023D00003, 75N92023D00004, 75N92023D00005, and 75N92023D00006 from the National Heart, Lung, and Blood Institute (NHLBI). The funders of this project were not involved in the writing, editing, or decision to publish this study. This manuscript was reviewed and approved by the CARDIA Publications and Presentations Committee.
Conflict of Interests
Maya Abdalla and Rim Elalami have no conflict of interests to disclose. Dr Cho reported receiving grants from Bayer. Dr. O’Connor reports no conflict of interests. Dr Rice reported receiving research grant funding from the NIH and receiving expert testimony fees from the Conservation Law Foundation. Dr. Horowitz reported consulting fees from Bayer and is a member of GE medical imaging advisory board. Dr. Akhoundi had no conflicts of interest to disclose. Dr. Yen reported receiving salary support. Dr Kalhan reported grants from AstraZeneca, PneumRx/BTG and Spiration to Northwestern University; consulting fees from CVS Caremark, AstraZeneca, GlaxoSmithKline, and CSA Medical; payment from GlaxoSmithKline, AstraZeneca, and Boehringer Ingelheim and is a member of CSA Medical advisory board. Dr Diaz reported receiving personal fees from Boehringer Ingelheim and having a patent for Methods and Compositions Relating to Airway Dysfunction pending (701586-190200USPT).