Abstract
Roughly 9% of African-Americans carry 1 copy of the sickle β-globin allele [sickle-cell trait (SCT)]. The sickle allele may affect the host responses in a way that protects against sequelae of malaria. Besides the association with sporadic episodes of erythrocyte sickling, acute chest syndrome, and sudden death, the SCT is generally regarded as asymptomatic. Separately, and for unclear reasons, the African-Americans have an elevated risk compared to whites for several immune-mediated diseases, including asthma, eczema, and other respiratory conditions. We propose to determine whether the SCT contributes to the elevated risk of these inflammatory conditions. We performed a retrospective cohort study of children with and without the sickle trait to assess the risk for common immune-mediated conditions. The cohorts were individually matched, and multiple logistic regression was used with variables selected using a backward stepwise approach. A second approach (case–control design) assessed the odds of the SCT among African-Americans with and without asthma. We found 2,481 children with and 4,962 matched patients without the SCT in the cohort design. The sickle trait was associated with greater odds for several immune-mediated conditions in a multivariable analysis, but not associated with asthma (odds ratio, OR=1.10, P=0.212). In an adjusted case–control analysis (n=20,000), the sickle trait was weakly associated with asthma [adjusted odds ratio (aOR)=1.46, 1.01–2.13]. In both the designs and in all statistical models used, we found that the SCT was associated with bronchitis (aOR=1.71, 1.09–2.67) and eczema (aOR=1.74, 1.23–2.46). The SCT may contribute to an increased risk for eczema and asthma-like symptoms in the African-American children.
Introduction
Since the African-Americans are at a greater risk for both the AS genotype and certain inflammatory disorders, and since the AS genotype is associated with an altered immune response to the Plasmodium falciparum parasite, it is rational to hypothesize that AS may promote the risk for other inflammatory conditions, as stated above. This question has been not been studied extensively. Palma-Carlos found the prevalence of asthma higher among 66 patients with milder hemoglobinopathies (predominantly the β-thalassemia carriers) compared to the controls. 32 Data from a second very small cohort suggested that AS was associated with asthma. 33 Utilizing a large multicenter electronic health record (EHR) database and 2 separate design approaches, we sought to determine if the SCT is associated with physician-diagnosed asthma, bronchitis, eczema, pneumonia, and sleep apnea.
Materials and Methods
This study was approved by the Nemours Florida Institutional Review Board [No. 183481-3]. We used a large EHR database involving Nemours hospitals and subspecialty clinics at and around Jacksonville, FL, Orlando, FL, Pensacola, FL, and Wilmington, DE, to construct a retrospective cohort study and a case–control study to determine the association between the SCT and select inflammatory conditions.
Study I: retrospective cohort design
Subject selection
We identified all patients seen between May, 2000, and December, 2010, at any of the 4 Nemours Children's Clinic outpatient sites with an International Classification of Disease (ICD)-9 diagnosis code 282.5 (SCT) in EHR diagnosis, problem list, medical history, or physician billing charge field (Epicare, Summer 2009, IU3, Epic Systems Corporation, Verona, WI). We excluded the patients who had an accompanying diagnosis code of 282.x (SCD, Sickle Cell-β-Thalassemia, and other Sickle Cell hemoglobinopathies) to target the patients likely to have the SCT. The SCT group was considered the exposure cohort (exposed to the S sickle allele). We collected data on race, ethnicity (Hispanic/non-Hispanic), age, Medicaid status, presence or absence of supplemental private insurance (as of November, 2010), and the clinic location. We categorized participants by age groupings: 0–4, 5–8, 9–12, 13–17, and 18 years and up.
We then created a matching no-exposure (non-SCT) group. These non-SCT patients were individually matched for the age group, sex, race, ethnicity, clinical location, Medicaid status, and supplemental insurance status. We randomly selected 2 controls for each SCT patient from a pool of non-SCT and non-SCD patients having the appropriate matching demographic profiles.
For both groups (SCT and non-SCT), we determined the cumulative incidence and odds of physician-diagnosed asthma, sinusitis, bronchiolitis, bronchitis, pneumonia, allergic rhinitis, snoring, sleep-disordered breathing, and obstructive sleep apnea (see the Supplementary data section for ICD-9 codes) using the EHR diagnosis and billing field.
Statistical analysis
We compared the baseline characteristics between the exposed and unexposed groups and used chi-square for any categorical outcome variables. We used Student's t-test, ANOVA, or Wilcoxon–Mann–Whitney U Rank-Sum test to compare the means or medians of continuous variables, as appropriate. SAS (SAS Institute, Cary, NC) and SPSS 11.0 (SPSS, Inc., Chicago, IL) were used for all computations.
We calculated the risk ratios and unadjusted odds ratios (ORs) for each diagnosis. We used multivariable logistic regression that included all variables associated with the SCT in a bivariate analysis. These included asthma, eczema, allergic rhinitis, sleep apnea, bronchitis, bronchiolitis, sinusitis, and snoring. A second model was used that involved variables selected by a stepwise backward selection in a simultaneous model. Multiple logistic regression models adjusted for multiple covariables to assess the influence of the SCT on individual odds of an inflammatory disease.
Study II: case–control design
Subject selection
We performed a second data query of African-American children, identifying all incident cases seen over the same period (May, 2000, through December, 2010) with an ICD-9 code of 493.x (Asthma) in the EHR diagnosis field, problem list field, medical history field, or physician billing field. We excluded the patients who also had a diagnosis code consistent with SCD, sickle-cell thalassemia, and other sickle-cell hemoglobinopathies/HbSC (282.41–282.42 and 282.6), so that we specifically identified the patients most likely to either have the SCT or be homozygous wild type (hemoglobin AA). Therefore, African-Americans with asthma constituted our case group. We collected information on ethnicity (Hispanic/non-Hispanic), age, Medicaid status, presence or absence of supplemental private insurance (as of November, 2010), and the most recently visited clinical location.
We then created a control group of African-American children without asthma based on individual matching. Controls could be either hemoglobin AS (SCT) or AA (non-SCT), but similar to cases, could not have any of the sickle-related hemoglobinopathies (SS, βS, or SC). We collected an individually matched control patient with the same characteristics for age group, gender, race, ethnicity, clinical location, Medicaid status, and supplemental insurance status. We randomly selected the controls from the pool of EHR patients without asthma and with the appropriate matching demographic profile.
Then for all cases and controls, we determined the prevalence and odds of having the SCT (282.5) documented among the EHR diagnosis field, problem list field, medical history field, or physician billing-charge field. In addition, we collected information on the presence or absence of sinusitis, bronchitis, pneumonia, eczema, allergic rhinitis, snoring, sleep-disordered breathing, and obstructive sleep apnea among the same EHR fields. See the Supplementary data section for ICD-9 codes.
Statistical analysis
We used stepwise logistic regression modeling to choose the predictor variables associated with increased odds of the SCT. Eligible predictor variables included the presence of sinusitis, bronchitis, pneumonia, allergic rhinitis, pneumonia, snoring, sleep-disordered breathing, obstructive sleep apnea, and eczema. We calculated the odds of having the SCT according to the asthma status. We created an adjusted multivariable logistic model that included bronchitis, eczema, body–mass index (BMI)-z score, and asthma status to assess the odds of the SCT.
Results
Table 1 lists the baseline characteristics for both the SCT and non-SCT cohorts from Study I. Because of individual matching, all characteristics between the 2 groups are similar and therefore are shown together. The exception was height and weight. Participants with the SCT were significantly shorter and weighed less despite being similar in age and sex. The mean heights [standard deviation (SD)] for the sickle trait group versus controls were 102.8 cm (32.8) and 107.2 cm (32.3), respectively, whereas the mean weights (SD) for the sickle trait group and controls were 22.7 kg (19.9) and 25.4 kg (22.1), respectively (P<0.001 for both comparisons).
P value comparing the sickle carrier cohort to noncarriers. P values were determined by chi-square or independent samples Student's t-test as appropriate.
SD, standard deviation; BMI, body–mass index.
Risk of respiratory disease in children by the sickle-allele status: Study I
Prevalences of any diagnosis of asthma, persistent asthma, and moderate or severe persistent asthma are shown in Table 2. Risk estimates for asthma, eczema, allergic rhinitis, obstructive sleep apnea, bronchiolitis, bronchitis, sleep-disordered breathing, snoring, and sinusitis were significantly elevated among children with the SCT (Table 2). In the adjusted models, we found significantly increased odds for eczema, allergic rhinitis, bronchitis, bronchiolitis, and sinusitis Table 3. The risk for persistent asthma was the same between the SCT and non-SCT patients (P=0.804). The relative risk and odds for all incident diagnoses were similar when we evaluated the African-Americans only and after adjusting for height, weight, and the BMI.
OR, odds ratio (unadjusted); CI, confidence interval; OSA, obstructive sleep apnea; SDB, sleep-disordered breathing; SE, standard error; SCT, sickle-cell trait.
All significant variables in the bivariable analysis (asthma, eczema, allergic rhinitis, OSA, bronchitis, bronchiolitis, SDB, snoring, and sinusitis) were used in the model. When asthma (presence/absence) was used in model one (instead of asthma severity), the OR (95% CI) was 1.10 (0.95, 1.26), P value=0.212.
Model includes all significant variables selected by stepwise backward logistic regression in a simultaneous model (and included eczema, allergic rhinitis, OSA, bronchiolitis, bronchitis, and sinusitis).
The number of clinical follow-up visits was similar among the SCT-versus-matched control patients (P>0.05, data not shown).
Odds of the SCT in African-Americans with asthma: Study II
We found more than 27,000 separate African-American children with asthma within our EHR database. Table 4 describes the background characteristics for the asthmatic cases and nonasthmatic controls combined. The baseline characteristics are shown combined, because cases and controls were matched for age group, gender, race, ethnicity, clinical location, Medicaid status, and supplemental insurance status. The groups were very similar due to the matching design; however, we found that the asthma group had a significantly higher BMI percentile and BMI-z score (data not shown, P<0.001 for both comparisons). Simple logistic regression showed several conditions, including physician-diagnosed asthma, bronchitis, eczema, sinusitis, and pneumonia, which were associated with increased odds of the SCT (Table 5). When we adjusted for these covariables, we found that only bronchitis and eczema were significantly associated with the SCT (Table 6). The point estimate of the OR for asthma was 1.20; however, this did not reach statistical significance in the multivariable model.
P value comparing the asthmatics versus nonasthmatics. P values were determined by chi-square or independent-sample Student's t-test as appropriate.
Associations with SCT among 10,000 asthmatic children and 10,000 nonasthmatics matched by race, age, sex, and insurance status.
Calculation for BMI-z score OR included 15,410 observations.
The variable BMI-z has 4,590 missing values.
A variable is removed from the model if the probability of the likelihood-ratio statistic using the conditional parameter estimates is>0.10.
Discussion
Our study suggests that children who have the SCT (β-globin genotype AS) may be at an increased risk for developing physician-diagnosed bronchitis and eczema, the 2 common inflammatory diseases of childhood. In both of the studies conducted, we saw a significantly increased risk estimate for bronchitis and eczema among children with the sickle trait compared to matched children without the SCT. In the cohort study (Study I), we assessed all physician visits involving children reporting to have the SCT. We performed a similar assessment on individually matched patients without the SCT. In both groups, we collected data on the physician diagnosis of asthma, bronchitis, eczema, and 6 other common inflammatory conditions. In 2 separate multivariable models, we found that eczema, bronchitis, bronchiolitis, allergic rhinitis, and sinusitis were more commonly diagnosed in the SCT group. In the second design (Study II), using a case–control approach, we found that only eczema and bronchitis were significantly associated with the SCT in the 2 separate multivariable models used. In total, we evaluated the data from more than 27,000 children from more than 10 sites in the Southeastern and mid-Atlantic regions of the United States. These epidemiologic data suggest the possibility of an association between the SCT and increased diagnosis of the inflammatory conditions, eczema, and bronchitis.
There are several limitations to the current study that must be considered carefully. It is important to note that this study merely suggests an epidemiologic link, and is far from establishing the causality or specific mechanisms linking the SCT with an immune or inflammatory dysfunction. Because these data were not collected prospectively and were not originally collected for the purpose of research, there exists the considerable possibility of bias and confounding. Bias commonly occurs in epidemiologic studies due to flawed measurement techniques, differential recall, and poor comparison group selection. An important limitation of the current study is the reliance on the EHR data that were created during the routine clinical care. For both the retrospective cohort analysis and the case–control analysis, the sickle trait status was not determined by genotyping or hemoglobin electrophoresis, but instead was determined by physician documentation of the SCT most likely coming from parental recall. There is some confusion among laypersons about the distinction between disease versus carrier status, and between different hemoglobinopathies that may have created a misclassification bias with regard to the risk group in the cohort study (and the exposure variable in the case–control study). Therefore, the risk-group (SCT versus normal) misclassification could have been present; however, the presence of a differential misclassification bias seems unlikely. The incident diseases in the cohort analysis were detected by a physician diagnosis documented in the EHR. These diagnoses also likely involved some degree of diagnostic misclassification, though differential diagnostic misclassification bias between the risk groups (SCT versus normal) again seems unlikely. Bias in this study was reduced by the use of an individual matching scheme that created near-identical control populations (for both study designs) in terms of race, age, sex, and access to healthcare. Potential factors of access to medical care, housing quality, and socioeconomic status should have been adequately controlled for and equal between the risk groups.
Replication of these findings is needed in other datasets and optimally in a prospective study that utilizes laboratory-based determination of the SCT status and supplemental testing to assess the markers of atopy and lung function. In addition, future studies using ancestral markers may be needed to reduce possible confounding by the African ancestry. Despite several methodologic limitations in the current study, there are several supporting lines of evidence that suggest that the links between the SCT and inflammatory conditions may be real. First, it is well established that SCD is a risk factor for small-airway obstruction,34,35 bronchial hyperresponsiveness, 11 and recurrent asthma-like symptoms.12,13,36,37 Since there are many reports of some children with the SCT displaying milder phenomena that are characteristic of true SCD (sickling and thrombosis), it would not be unexpected for the SCT to be associated with some element of airway dysfunction. Second, the few studies that have tried to evaluate the risk of asthma in the SCT have found results consistent with an epidemiologic link,32,33 although these studies have been small and could represent reporting bias and spurious findings.
Lastly, it is important to note that there is biologic plausibility to the SCT, promoting inflammatory disease. In vitro evidence suggests that the SCT may affect immunoglobulin, T-lymphocyte, and eosinophil functioning, all of which are components of asthma pathogenesis.38–46 In the malaria endemic regions, the SCT genotype appears to lessen the severity of malaria,38,39 possibly through accelerated antibody production, 40 enhanced antibody and complement binding, 41 and immune-mediated parasite clearance. 42 The SCT has been associated with upregulation of certain cell-mediated immune responses,43,44 eosinophilia, 45 and enhanced immune recognition of P. falciparum-infected RBCs. 46 It is possible that the heightened immune responses seen in the setting of endemic P. falciparum malaria may also contribute to altered inflammatory responses. One study looking at the SCT in a P. falciparum endemic region found no altered risk for upper- or lower-respiratory-tract infections. 40 However, within a nonendemic malarial region, Poehling found that children with the SCT may be at a greater risk for invasive pneumococcal disease, including blood, spinal, and parapneumonic infections. 47
The SCT has also been associated with reductions in somatic growth among children 48 and among newborns delivered by the SCT carriers. 49 Subsequent data have cast a doubt on whether the SCT truly affects somatic growth,50,51 and indeed, we also found that the SCT was not associated with somatic growth. Reduced somatic growth in early life (and subsequent reduced lung growth), therefore, is an unlikely explanation for our findings of greater bronchitis symptoms. Sensitive evaluations of the lung mechanics in children with the SCT may uncover growth alterations; however, research in this area is lacking. Similarly, it is rational to hypothesize that alterations in the immune function seen with the SCT in the context of parasitic infection may also be present during viral and bacterial infections that are common in early childhood. Immune modulation may be a mediating factor between the SCT and the risk for eczema and asthma-like symptoms. Since bronchitis in children lacks a clear set of diagnostic criteria, it is hard to know if the link with bronchitis represents some degree of misclassified asthma. Further prospective analysis with detailed phenotyping is needed to answer this question.
In summary, it remains unclear why African-American children (an ancestral group with a high S-allele frequency) have a higher prevalence for eczema, asthma symptoms, and respiratory infections. These common infectious and immune-related pediatric conditions represent a large burden to the primary care and acute care systems in the United States and worldwide. It is possible that a contributing factor is the S allele. The current study has significant limitations due to its retrospective design using EHR data, but it signals the need for more and higher-quality investigation, including a genotype-stratified cohort study prospectively after the relative incidence of asthma, bronchitis, and eczema.
Authors' Contributions
Dr. Lang has made substantial contributions to conception, design, and acquisition of data for this study. He contributed to the analysis and interpretation of data. He drafted the submitted article. Dr. Lang also vouches for the integrity of the data and accuracy of the data analysis. He is a guarantor of the article, taking responsibility for the integrity of the work as a whole, from inception to the final manuscript.
Dr. Hossain has made substantial contributions to the analysis and interpretation of data; He has revised the article critically for important intellectual content.
Dr. Blake made substantial contributions to the design and acquisition of data for this study. She helped critically to edit the manuscript for important intellectual content.
Arnel Mercado made substantial contributions to the design and acquisition of data for this study. He helped critically to edit the manuscript for important intellectual content.
Dr. Lima has made substantial contributions to the conception, design, and acquisition of data for this study. He contributed to the analysis and interpretation of data. He helped draft the submitted article and revised it critically for important intellectual content.
Footnotes
Author Disclosure Statement
The authors have no financial relationships or conflicts of interest relevant to this article to disclose.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
