Abstract
Severe asthma exacerbations are a major cause of school absences and healthcare costs in children, particularly those in high-risk racial/ethnic groups.
To identify susceptibility genes for severe asthma exacerbations in Latino children and adolescents, we conducted a meta-analysis of genome-wide association studies (GWAS) in 4010 Latino youth with asthma in four independent cohorts, including 1693 Puerto Ricans, 1019 Costa Ricans, 640 Mexicans, 256 Brazilians and 402 members of other Latino subgroups. We then conducted methylation quantitative trait locus, expression quantitative trait locus and expression quantitative trait methylation analyses to assess whether the top single nucleotide polymorphism (SNP) in the meta-analysis is linked to DNA methylation and gene expression in nasal (airway) epithelium in separate cohorts of Puerto Rican and Dutch children and adolescents.
In the meta-analysis of GWAS, an SNP in FLJ22447 (rs2253681) was significantly associated with 1.55 increased odds of severe asthma exacerbation (95% CI 1.34–1.79, p=6.3×10−9). This SNP was significantly associated with DNA methylation of a CpG site (cg25024579) at the FLJ22447 locus, which was in turn associated with increased expression of KCNJ2-AS1 in nasal airway epithelium from Puerto Rican children and adolescents (β=0.10, p=2.18×10−7).
SNP rs2253681 was significantly associated with both DNA methylation of a cis-CpG in FLJ22447 and severe asthma exacerbations in Latino youth. This may be partly explained by changes in airway epithelial expression of a gene recently implicated in atopic asthma in Puerto Rican children and adolescents (KCNJ2-AS1).
Abstract
A novel SNP in FLJ22447 is significantly associated with both severe asthma exacerbations in Latino youth and DNA methylation of a cis-CpG in nasal epithelium. This CpG is linked to nasal epithelial expression of a gene implicated in atopic asthma. https://bit.ly/33D9Joc
Introduction
Asthma is the most common chronic respiratory disease among children [1]. In the USA, total costs related to asthma exceed $81 billion per year [2]. Severe asthma exacerbations (SAEs), defined as episodes of disease worsening that require a change in treatment to prevent a serious outcome [3], are a major cause of school or work absences and healthcare costs. Of the ∼6.1 million children with asthma in the USA, 2.1% report ≥1 asthma-related hospitalisation and 10.7% report ≥1 asthma-related visit to the emergency department (ED) in the previous year [4]. Despite recent advances, the best predictor of SAEs is having had an SAE in the previous year [5].
Although genome-wide association studies (GWASs) have identified susceptibility loci for asthma [6–8], little is known about genetic determinants of asthma exacerbations, which may be distinct from those for asthma per se. A GWAS in Danish children with asthma (ages 2–6 years) identified cadherin-related family member 3 (CDHR3), a gene not previously associated with asthma, as a susceptibility locus for recurrent SAEs [9]. In a combined GWAS of two cohorts of non-Hispanic white children with asthma, four intronic single nucleotide polymorphisms (SNPs) in cadherin-associated protein alpha 3 (CTNNA3) were significantly associated with SAEs [10], but the association was not replicated in an independent cohort including 786 children and adults with asthma. Moreover, a GWAS of SAEs among 806 non-Hispanic white children and adults with asthma who were being treated with inhaled corticosteroids found no genome-wide significant results [11].
The burden of asthma varies across racial and ethnic groups in the USA and Latin America. For example, Puerto Rican children have a greater prevalence of and morbidity and mortality from asthma than non-Hispanic white children in the USA [12], and Costa Rican adolescents have a greater burden of asthma than those in other Latin American countries [13]. Moreover, recent evidence suggests that some susceptibility variants for asthma-related outcomes are ethnicity-specific [14, 15]. We hypothesised that there would be susceptibility variants for SAEs that would be more common or exert a greater effect in Latino subgroups at risk for morbidity from asthma. To test this hypothesis, we conducted a meta-analysis of GWAS of SAEs among Latino youth with asthma in four independent studies.
Methods
Please refer to the supplementary material for more details.
Study populations included in the meta-analysis of GWAS of severe exacerbations
Hartford-Puerto Rico study
The Hartford-Puerto Rico (HPR) study is a case–control study of childhood asthma in Puerto Ricans [6]. SAEs were defined as ≥1 hospitalisation for asthma or ≥1 ED/urgent care visit for asthma requiring treatment with systemic corticosteroids in the previous year, or ≥1 course of systemic corticosteroids for asthma in the previous year. After quality control (QC) measures, 554 independent children with asthma (236 of whom had ≥1 SAE in the previous year) were included in the genome-wide association analysis, which was conducted using logistic regression under an additive genetic model, adjusting for age, sex, inhaled steroid use and the first two principal components (PCs), calculated using smartPCA [16].
The Genetics of Asthma in Latino Americans study
The Genetics of Asthma in Latino Americans (GALA II) study is a case–control study of asthma in Latino children and youth [17]. An SAE was defined as having had ≥1 of the following events in the previous year: hospitalisations for asthma, ED visits or unscheduled and urgent doctor's visits because of asthma, or treatment with systemic corticosteroids for asthma. The current analysis focused on 2181 children with asthma (1283 of whom had ≥1 SAE) and self-reported Latino ethnicity (1139 Puerto Rican, 640 Mexican and 402 from other groups). Association testing was conducted using logistic regression under an additive genetic model, adjusting for age, sex, inhaled steroid use, ethnicity and the first two PCs.
The Genetics of Asthma in Costa Rica Study
The Genetics of Asthma in Costa Rica Study (GACRS) is a genetic study of nuclear families of children with asthma in Costa Rica [18, 19]. An SAE was defined as ≥1 ED visit for asthma or ≥1 hospitalisation for asthma in the previous year. After QC, 1019 independent children with asthma (851 with ≥1 SAE) were included in the analysis, which was also conducted using logistic regression under an additive genetic model, adjusting for age, sex, inhaled steroid use and the first two PCs.
The Social Changes, Asthma and Allergy in Latin America study, Bahia, Brazil
The Social Changes, Asthma and Allergy in Latin America (SCAALA) study is a longitudinal study of asthma and allergic diseases of Brazilian children [20]. Children who reported asthma symptoms and had data on SAEs (n=256) were included in this analysis. An SAE was defined as ≥1 ED visit or ≥1 hospitalisation due to wheeze in the previous year. The analysis was conducted using logistic regression under an additive genetic model, with adjustment for age, sex and the first two PCs from genotypic data.
Study populations included in molecular quantitative trait analyses in nasal epithelium
The Epigenetic Variation and Childhood Asthma in Puerto Rico study
Nasal (airway) epithelium can serve as a surrogate marker for DNA methylation and gene expression in bronchial (airway) epithelium [21]. In the Epigenetic Variation and Childhood Asthma in Puerto Rico (EVA-PR) study, whole-genome methylation assays were performed using HumanMethylation450 BeadChips (Illumina, San Diego, CA, USA). After QC, 227 901 CpG probes remained for the analysis of nasal epithelium, and M-values were used in all downstream analyses. RNA sequencing was performed with the Illumina NextSeq 500 platform, with paired-end reads at 75 cycles and 80 million reads per sample. After QC, 16 737 genes were retained for the analysis. The R function sva was used to estimate latent factors that capture unknown data heterogeneity [22]. Of the 543 study participants, 457 had complete genome-wide data for genotypes, methylation and transcriptomics in nasal epithelium (supplementary figure S1).
Prevention and Incidence of Asthma and Mite Allergy study
The Prevention and Incidence of Asthma and Mite Allergy (PIAMA) study is a birth cohort study of children born in the Netherlands in 1996 and 1997 [23, 24]. A total of 479 nasal epithelial samples were hybridised to the Infinium HumanMethylation450 BeadChip arrays. After QC, 455 samples and 436 824 probes remained, and 432 samples had matched genotype data. RNA sequencing was performed with the Illumina HiSeq 2500 platform with paired-end sequencing. After stringent QC, 17 156 genes and 326 samples remained, with 233 samples having matched genotype data [25].
Meta-analysis of GWAS of SAEs
METAL [26] software was used to perform the meta-analysis of GWAS of SAEs, using data from HPR, GALA II, GACRS and SCAALA.
Molecular quantitative trait locus analyses
To estimate the effects of our top SNP on DNA methylation and gene expression, we conducted a methylation quantitative trait locus (mQTL) analysis to test for an association between our top SNP and genome-wide DNA methylation in nasal epithelium, and an expression quantitative trait locus (eQTL) analysis to test for an association between our top SNP and genome-wide gene expression in nasal epithelium. We then conducted an expression quantitative trait methylation (eQTM) analysis to test whether the top CpG site identified in the mQTL analysis was associated with genome-wide gene expression in nasal epithelium. All analyses were adjusted for age, sex, asthma status, atopy status, the top five PCs from genotypic data and latent factors estimated from sva [22]. In addition, RNA sequencing batch and RNA sample sorting protocol (i.e. whole cells or CD326-positive nasal epithelial cells) [21] were adjusted for in the eQTM analyses, and methylation batch was adjusted for in the mQTL and eQTM analyses; both cis and trans effects were considered.
Pathway analysis
A pathway analysis was performed with MAGMA [27], which conducts SNP-wise gene analysis of summary statistics with correction for linkage disequilibrium between variants and genes, to test whether sets of genes are jointly associated with a phenotype (i.e. asthma exacerbation) compared to other genes across the genome. Adaptive permutation was used to produce an empirical p-value and false discovery rate (FDR). Gene sets used in the analyses were from Gene Ontology (GO) [28, 29], Kyoto Encyclopedia of Genes and Genomes (KEGG) [30, 31], Reactome [32, 33] and BioCarta pathways.
Results
The characteristics of the 4010 participants in the four studies included in the meta-analysis of GWAS of SAEs are shown in table 1. In total, there were 2509 children with asthma and ≥1 SAE (cases) and 1501 children with asthma but no SAEs (controls). Compared to subjects who participated in the HPR or SCAALA studies, those in GALA II were older and those in GACRS were younger. Compared to subjects in the other studies, those in GACRS were more likely to be male and to have had ≥1 severe asthma exacerbation in the previous year.
Approximately six million genotyped and imputed SNPs with a minor allele frequency (MAF) ≥0.05 were included in the meta-analysis of GWAS of SAEs. In this meta-analysis, one SNP (rs2253681, in FLJ22447 on chromosome 14q23.2) was significantly associated with SAEs (p<5×10−8) (figure 1a). This SNP was genotyped in HPR, GALA II and SCAALA, and imputed in GACRS with imputation quality r2=0.96. Each copy of the minor allele (A) of SNP rs2253681 was associated with 1.55 times increased odds of SAEs (95% CI 1.34–1.79, p=6.3×10−9, figure 1b). The Q-Q plot (figure 1c) showed no inflation for the results for each of the four individual cohorts or for the pooled analysis (figure 2). There was no significant interaction between SNP rs2253681 and inhaled steroid use on SAEs in the GWAS conducted in the HPR, GALA II or GACRS cohorts (p for interaction≥0.15 in all instances).
We then examined whether SNPs previously associated with asthma in either a multi-ancestry meta-analysis [7] or a meta-analysis from the UK Biobank [34] were associated with SAEs in our analysis (supplementary table S1). None of the previously reported asthma-susceptibility SNPs was significantly associated with SAEs in our meta-analysis of Latino youth. Moreover, no SNPs associated with SAEs or asthma hospitalisations in previous candidate-gene studies (e.g. IL13, IL4RA) or GWAS (e.g. CDHR3, CTNNA3) [9, 10, 18, 35–38] were significantly associated with SAEs in our meta-analysis (meta-p≥0.05 in all instances, supplementary table S2). Although the previously reported SNP rs7216389 in ORMDL3 on chromosome 17q21 was nominally associated with SAEs in HPR (p=1.6×10−3), it did not reach statistical significance in our meta-analysis (meta-p=0.06, supplementary table S2).
To examine whether SNP rs2253681 affects methylation of the FLJ22447 locus in nasal epithelium, we first conducted an mQTL analysis. In this analysis, rs2253681 was a significant cis-acting mQTL in nasal epithelium for cg25024579 in FLJ22447 (β=0.55, p=3.6×10−16 and FDR-p=8.1×10−11 in EVA-PR; β=0.34, p=1.2×10−11 in PIAMA; and meta-p=9.5×10−25 for a combined analysis of EVA-PR and PIAMA; table 2). In addition, among the top 20 mQTL, cg05223396 in SH3PXD2B, cg13127574 in ADSSL1 and cg21115391 in CCDC67 were associated with rs2253681 at p<0.05 in PIAMA, in the same direction as in EVA-PR (table 2). Next, we conducted an eQTL analysis of SNP rs2253681 in the FLJ22447 locus and gene expression in nasal epithelium, and found no cis or trans genome-wide significant results (supplementary table S3). Among the top 20 eQTL, SRF was associated with rs2253681 at p<0.05 in PIAMA, in the same direction as in EVA-PR (supplementary table S3).
Given the association between SNP rs2253681 and methylation of cg25024579, we then examined the effect of that CpG on gene expression in nasal epithelium by conducting an eQTM analysis (table 3). In this analysis, methylation of cg25024579 was significantly associated in trans with expression of KCNJ2 antisense RNA 1 (KCNJ2-AS1) on chromosome 17q24.3 (β=0.10, p=2.2×10−7 and FDR-p=3.2×10−3) in EVA-PR. Although this was not replicated in PIAMA, the observed association remained significant in the combined analysis of the two cohorts (β=0.09, p=1.4×10−6). The top cis-eQTM gene in the current analysis in EVA-PR was protein kinase C-η (PKCη), encoded by PRKCH (β=0.07, p=4.1×10−5 and FDR-p=1.5×10−1). PRKCH expression was also associated with cg2504579 methylation in PIAMA, in the same direction of association as in EVA-PR (β=0.11, p=1.4×10−2; and p=2.4×10−6 in the combined analysis of EVA-PR and PIAMA). Among the top 20 eQTM in EVA-PR, SLC27A2 expression was associated with cg25024579 methylation in PIAMA, in the same direction of association as in EVA-PR (β=0.22, p=4.2×10−5 in PIAMA; and β=0.09, p=4.9×10−7 in the combined analysis of the two cohorts). In addition, expression of FOCAD, PIGF and ERLEC1 was associated with cg25024579 methylation in PIAMA at p<0.05, in the same direction as in EVA-PR.
To further assess the biological relevance of the FLJ22447 locus, we conducted an SNP-wise pathway analysis of SAEs using the meta-analysis results and evaluated public repositories and databases. Although there were no pathways associated with SAEs after adjusting for multiple testing, there were a number of nominally significant pathways that included genes on chromosome 14q23.2, near FLJ22447 (supplementary table S4).
In a sensitivity analysis, we repeated the GWAS of SAEs after imputing genotypes using the 1000 Genomes Ad Mixed American (AMR) reference panel instead of the Haplotype Reference Consortium (HRC) r1.1 2016 reference panel (supplementary material), obtaining very similar results (supplementary figure S2 and supplementary table S5).
Discussion
Our combined analysis including 4010 youths in four study cohorts showed that a novel SNP in FLJ22447 (rs2253681) is significantly associated with SAEs (p=6.3×10−9) among children and adolescents in Latino subgroups affected with asthma (Puerto Ricans, Costa Ricans, Mexicans and Brazilians) [12]. In an mQTL analysis in nasal epithelium, this SNP was significantly associated with DNA methylation of a CpG site at the FLJ22447 locus (cg25024579), which was in turn significantly associated with the expression of KCNJ2-AS1 in nasal epithelium.
FLJ22447 codes for a long non-coding RNA, a type of RNA that does not get translated to a protein but may serve pre- and post-transcriptional regulatory functions. FLJ22447 partially overlaps with AL355916.3, which encodes a protein kinase C paralog expressed in epithelial tissues. Of note, an SNP in AL355916.3 has been associated with forced expiratory volume in 1 s/forced vital capacity and bronchodilator response in adults [34]. FLJ22447 mediates interleukin (IL)-33 upregulation in activated fibroblasts; together, FLI22447 and IL-33 can increase fibroblast expression of α-smooth muscle actin (α-SMA), vimentin and N-cadherin [35], which are associated with myofibroblast differentiation and airway remodelling [36, 37].
KCNJ2-AS1 also encodes a long non-coding RNA. Expression of KCNJ2-AS1 (log2(fold-change)=0.34, p=3.9×10−6, FDR-p=3.1×10−5) was significantly associated with atopic asthma in our recent transcriptome-wide association study of atopic asthma in Puerto Rican children and adolescents [38]. SNPs rs8066985 and rs312750 near KCNJ2-AS1 are associated with body mass index and waist-to-hip ratio [39, 40], which are associated with worse asthma outcomes in children and adults. To our knowledge, this is the first study linking FLJ22447 or KCNJ2-AS1 to SAEs in children and adolescents.
While not significant genome wide, the top cis-eQTM gene in EVA-PR was PRKCH, which is adjacent to FLJ22447 and most highly expressed in the lungs [41]. Both methylation and transcription of this gene in nasal (airway) epithelium were recently associated with atopic asthma in Puerto Rican children and adolescents [21]. PKCη, encoded by PRKCH, plays a key role in the assembly and maintenance of epithelial tight junctions. PKCη phosphorylates occludin on threonine residues (T403 and T404), and such phosphorylation is required for the assembly and/or maintenance of occludin in epithelial tight junctions that are key to the integrity and function of the airway epithelial barrier [42].
Recent GWAS of exacerbations in non-Hispanic white children identified SNPs in two genes, CDHR3 and CTNNA3 [9, 10], neither of which replicated in our analysis, despite our having similar or larger sample sizes than those in the original studies. Similarly, a previously reported association between an SNP on chromosome 17q21 (rs7216389) and SAEs was not statistically significant in our meta-analysis (p=0.06). While an association between the 17q21 locus and asthma is widely recognised, whether this locus is linked to severe disease exacerbations among children with asthma is less clear. The initial report in paediatric asthma analysed a cohort of 376 children and reported 80 with recurrent wheezing and 66 with asthma [43]. The authors then analysed the rate of severe exacerbations within the full cohort (57 among the 376 participants), rather than comparing children with asthma and an exacerbation to those with asthma but no exacerbation. A recent meta-analysis reported a significant association between rs7216389 and asthma hospitalisations and ED visits in children, but, while the pooled estimates were significant, only four of the 13 cohorts showed significant results [44]. Although the lack of an association between CDHR3 and SAEs in the current study may be due to differences in the age of participants and outcome definitions across studies, our results for FLJ22447 are novel and further highlight the importance of studying asthma outcomes in children and adolescents of diverse races and ethnicities.
We recognise several study limitations. First, we lack data on viral infections or air pollution, which may interact with genetic or epigenetic mechanisms causing SAEs in children. Second, we cannot assess temporal relationships between DNA methylation or gene expression and SAEs in our cross-sectional analysis. Third, we had insufficient statistical power to detect either uncommon risk alleles or alleles with modest genetic effects on SAEs. Moreover, we had limited statistical power to detect an interaction between our top SNP and inhaled steroid use on SAEs or for our molecular quantitative trait analyses, and lacked a replication cohort for analyses of DNA methylation or gene expression in nasal epithelium and SAEs (because PIAMA lacks adequate data on asthma exacerbations). Fourth, the definitions of SAEs varied across the study cohorts. However, we obtained similar results in a sensitivity analysis for the HPR and GALA II cohorts, in which we re-ran the GWAS after excluding subjects who received systemic corticosteroids but did not report an unscheduled and acute visit for asthma (data not shown).
In summary, we have identified a novel SNP in FLJ22447 that is significantly associated with both SAEs in Latino children and adolescents and DNA methylation of a cis-CpG in FLJ22447 in nasal epithelium. This CpG is, in turn, linked to nasal epithelial expression of a gene recently implicated in atopic asthma in children and adolescents (KCNJ2-AS1). Future longitudinal studies of asthma-omics should assess whether and how genetic variants affect airway epithelial function and SAEs in childhood.
Supplementary material
Supplementary Material
Please note: supplementary material is not edited by the Editorial Office, and is uploaded as it has been supplied by the author.
Supplementary material ERJ-02693-2020.SUPPLEMENT
Shareable PDF
Supplementary Material
This one-page PDF can be shared freely online.
Shareable PDF ERJ-02693-2020.Shareable
Footnotes
This article has supplementary material available from erj.ersjournals.com
Author contributions: Q. Yan, E. Forno and J.C. Celedón conceived and designed the study. E.G. Burchard, G.H. Koppelman and J.C. Celedón obtained funding. Q. Yan, E. Herrera-Luis, M. Pino-Yanes and R. Rios conducted the primary analysis and interpreted data. C. Qi, Y-Y. Han, S. Kim, S. Oh, E. Acosta-Pérez, R. Zhang, D. Hu, C. Eng, S. Huntsman, L. Avila, N. Boutaoui, M.M. Cloutier, M.E. Soto-Quiros, C-J. Xu, S.T. Weiss, J. Lasky-Su, M.R. Kiedrowski, C. Figueiredo, J. Bomberger, M.L. Barreto, G. Canino, W. Chen and E.G. Burchard participated in data collection and data analysis. Q. Yan, E. Forno and J.C. Celedón prepared the first draft of the manuscript. All authors reviewed the draft for intellectual content, and approved submission of the final version of the manuscript.
Support statement: This work was supported by grants HL079966, HL117191 and MD011764 from the US National Institutes of Health (NIH) to J.C. Celedón. Q. Yan's contribution was supported by grant HL138098 from the US NIH. E. Forno's contribution was supported by grant HL149693 from the US NIH. The GALA II study was supported by US NIH grants to E.G. Burchard: HL088133, HL004464, HL117004, ES015794, ES24844, TRDRP 24RT 0025, MD006902 and GM007546. M. Pino-Yanes was supported by the Ramón y Cajal Program (RYC-2015-17205) and by grant SAF2017-83417R from the Spanish Ministry of Economy, Industry and Competitiveness. The PIAMA study was supported by The Netherlands Organisation for Health Research and Development; The Netherlands Organisation for Scientific Research; the Lung Foundation of the Netherlands (with methylation and gene expression studies supported by AF 4.1.14.001); The Netherlands Ministry of Spatial Planning, Housing and the Environment; and The Netherlands Ministry of Health, Welfare and Sport. C. Qi was supported by a grant from the China Scholarship Council. Funding information for this article has been deposited with the Crossref Funder Registry.
Conflict of interest: Q. Yan reports grants from National Institutes of Health (K01HL138098), during the conduct of the study.
Conflict of interest: E. Forno has nothing to disclose.
Conflict of interest: E. Herrera-Luis reports grants from the Spanish Ministry of Science, Innovation and Universities (PRE2018-083837), during the conduct of the study.
Conflict of interest: M. Pino-Yanes reports grants from Spanish Ministry of Economy, Industry and Competitiveness, the State Research Agency and the European Regional Development Funds from the European Union (MINECO/AEI/FEDER, UE, SAF2017-83417R), and grants from Spanish Ministry of Economy, Industry and Competitiveness (Ramón y Cajal grant RYC-2015-17205), during the conduct of the study.
Conflict of interest: C. Qi has nothing to disclose.
Conflict of interest: R. Rios has nothing to disclose.
Conflict of interest: Y-Y. Han has nothing to disclose.
Conflict of interest: S. Kim has nothing to disclose.
Conflict of interest: S. Oh has nothing to disclose.
Conflict of interest: E. Acosta-Pérez has nothing to disclose.
Conflict of interest: R. Zhang has nothing to disclose.
Conflict of interest: D. Hu has nothing to disclose.
Conflict of interest: C. Eng has nothing to disclose.
Conflict of interest: S. Huntsman has nothing to disclose.
Conflict of interest: L. Avila has nothing to disclose.
Conflict of interest: N. Boutaoui has nothing to disclose.
Conflict of interest: M.M. Cloutier has nothing to disclose.
Conflict of interest: M.E. Soto-Quiros has nothing to disclose.
Conflict of interest: C-J. Xu has nothing to disclose.
Conflict of interest: S.T. Weiss has nothing to disclose.
Conflict of interest: J. Lasky-Su has nothing to disclose.
Conflict of interest: M.R. Kiedrowski has nothing to disclose.
Conflict of interest: C. Figueiredo has nothing to disclose.
Conflict of interest: J. Bomberger has nothing to disclose.
Conflict of interest: M.L. Barreto has nothing to disclose.
Conflict of interest: G. Canino has nothing to disclose.
Conflict of interest: W. Chen has nothing to disclose.
Conflict of interest: G.H. Koppelman reports grants from the Lung Foundation of the Netherlands, during the conduct of the study; grants from Ubbo Emmius Foundation, TETRI Foundation, GSK, Vertex, TEVA, the Lung Foundation of the Netherlands and the European Union H2020 Program, outside the submitted work; and has served on advisory boards to GSK and PURE IMS (fees paid to institution).
Conflict of interest: E.G. Burchard reports grants from National Heart, Lung, and Blood Institute, National Institute of Health, General Medical Sciences, National Institute on Minority Health and Health Disparities, Tobacco-Related Disease Research Program, Food and Drug Administration, and the Sandler Family Foundation, during the conduct of the study.
Conflict of interest: J.C. Celedón has received research materials from GSK and Merck (inhaled steroids) and Pharmavite (vitamin D and placebo capsules), in order to provide medications free of cost to participants in NIH-funded studies, unrelated to the current work.
- Received July 15, 2020.
- Accepted September 30, 2020.
- Copyright ©ERS 2021