drug response; genetic variants; pharmacogenomics; vitamin D receptor
Rationale: Low vitamin D levels are associated with asthma and decreased airway responsiveness. Treatment with inhaled corticosteroids improves airway responsiveness and asthma control.
Objectives: To assess the effect of vitamin D levels on prebronchodilator FEV1, bronchodilator response, and responsiveness to methacholine (PC20, provocative concentration of methacholine producing a 20% decline in FEV1) in patients with asthma treated with inhaled corticosteroids.
Methods: We measured 25-hydroxyvitamin D levels in the serum of children with persistent asthma at the time of enrollment in the Childhood Asthma Management Program. We divided subjects into the vitamin D sufficiency (>30 ng/ml), insufficiency (20–30 ng/ml), and deficiency (<20 ng/ml) groups. Covariates included age, treatment, sex, body mass index, race, history of emergency department visits, hospitalizations, and season that vitamin D specimen was drawn. Our main outcome measures were change in prebronchodilator FEV1, bronchodilator response, and PC20 from enrollment to 8–12 months.
Measurements and Main Results: Of the 1,024 subjects, 663 (65%) were vitamin D sufficient, 260 (25%) were insufficient, and 101 (10%) were deficient. Vitamin D–deficient subjects were more likely to be older, African American, and have a higher body mass index compared with the vitamin D–sufficient and insufficient subjects. In the inhaled corticosteroid treatment group, prebronchodilator FEV1 increased from randomization to 12 months by 140 ml in the vitamin D–deficient group and prebronchodilator FEV1 increased by 330 ml in the vitamin D insufficiency group and by 290 ml in the vitamin D sufficiency group (P = 0.0072), in adjusted models.
Conclusions: In children with asthma treated with inhaled corticosteroids, vitamin D deficiency is associated with poorer lung function than in children with vitamin D insufficiency or sufficiency.
asthma; vitamin D; lung function; forced expiratory volume; children
Vitamin D deficiency is becoming more apparent in many populations. Genetic factors may play a role in the maintenance of vitamin D levels. The objective of this study was to perform a genome-wide analysis (GWAS) of vitamin D levels, including replication of prior GWAS results. We measured 25-hydroxyvitamin D (25(OH)D) levels in serum collected at the time of enrollment and at year 4 in 572 Caucasian children with asthma, who were part of a multi-center clinical trial, the Childhood Asthma Management Program. Replication was performed in a second cohort of 592 asthmatics from Costa Rica and a third cohort of 516 Puerto Rican asthmatics. In addition, we attempted replication of three SNPs that were previously identified in a large GWAS of Caucasian individuals. The setting included data from a clinical trial of childhood asthmatics and two cohorts of asthmatics recruited for genetic studies of asthma. The main outcome measure was circulating 25(OH)D levels. The 25(OH)D levels at the two time-points were only modestly correlated with each other (intraclass correlation coefficient = 0.33) in the CAMP population. We identified SNPs that were nominally associated with 25(OH)D levels at two time-points in CAMP, and replicated four SNPs in the Costa Rican cohort: rs11002969, rs163221, rs1678849, and rs4864976. However, these SNPs were not significantly associated with 25(OH)D levels in a third population of Puerto Rican asthmatics. We were able to replicate the SNP with the strongest effect, previously reported in a large GWAS: rs2282679 (GC), and we were able to replicate another SNP, rs10741657 (CYP2R1), to a lesser degree. We were able to replicate two of three prior significant findings in a GWAS of 25(OH)D levels. Other SNPs may be additionally associated with 25(OH)D levels in certain populations.
Patients with severe or difficult-to-treat asthma are an understudied population but account for considerable asthma morbidity, mortality, and costs. The Epidemiology and Natural History of Asthma: Outcomes and Treatment Regimens (TENOR) study was a large, 3-year, multicenter, observational cohort study of 4756 patients (n = 3489 adults ≥18 years of age, n = 497 adolescents 13-17 years of age, and n = 770 children 6-12 years of age) with severe or difficult-to-treat asthma. TENOR's primary objective was to characterize the natural history of disease in this cohort. Data assessed semiannually and annually included demographics, medical history, comorbidities, asthma control, asthma-related health care use, medication use, lung function, IgE levels, self-reported asthma triggers, and asthma-related quality of life. We highlight the key findings and clinical implications from more than 25 peer-reviewed TENOR publications. Regardless of age, patients with severe or difficult-to-treat asthma demonstrated high rates of health care use and substantial asthma burden despite receiving multiple long-term controller medications. Recent exacerbation history was the strongest predictor of future asthma exacerbations. Uncontrolled asthma, as defined by the 2007 National Heart, Lung, and Blood Institute guidelines’ impairment domain, was highly prevalent and predictive of future asthma exacerbations; this assessment can be used to identify high-risk patients. IgE and allergen sensitization played a role in the majority of severe or difficult-to-treat asthmatic patients.
TENOR; severe or difficult-to-treat asthma; asthma control; asthma exacerbations; burden; medication; quality of life; allergy; IgE
Rationale: Variability in pulmonary disease severity is found in patients with cystic fibrosis (CF) who have identical mutations in the CF transmembrane conductance regulator (CFTR) gene. We hypothesized that one factor accounting for heterogeneity in pulmonary disease severity is variation in the family of genes affecting the biology of interleukin-1 (IL-1), which impacts acquisition and maintenance of Pseudomonas aeruginosa infection in animal models of chronic infection. Methods: We genotyped 58 single nucleotide polymorphisms (SNPs) in the IL-1 gene cluster in 808 CF subjects from the University of North Carolina and Case Western Reserve University (UNC/CWRU) joint cohort. All were homozygous for ΔF508, and categories of “severe” (cases) or “mild” (control subjects) lung disease were defined by the lowest or highest quartile of forced expired volume (FEV1) for age in the CF population. After adjustment for age and gender, genotypic data were tested for association with lung disease severity. Odds ratios (ORs) comparing severe versus mild CF were also calculated for each genotype (with the homozygote major allele as the reference group) for all 58 SNPs. From these analyses, nine SNPs with a moderate effect size, OR ≤ 0.5or > 1.5, were selected for further testing. To replicate the case-control study results, we genotyped the same nine SNPs in a second population of CF parent-offspring trios (recruited from Children’s Hospital Boston), in which the offspring had similar pulmonary phenotypes. For the trio analysis, both family-based and population-based associations were performed. Results: SNPs rs1143634 and rs1143639 in the IL1B gene demonstrated a consistent association with lung disease severity categories (P < 0.10) and longitudinal analysis of lung disease severity (P < 0.10) in CF in both the case-control and family-based studies. In females, there was a consistent association (false discovery rate adjusted joint P-value < 0.06 for both SNPs) in both the analysis of lung disease severity in the UNC/CWRU cohort and the family-based analysis of affection status. Conclusion: Our findings suggest that IL1β is a clinically relevant modulator of CF lung disease.
gene modifiers; cystic fibrosis; CFTR; IL-1 gene family
IL10 is an anti-inflammatory cytokine that has been found to have lower production in macrophages and mononuclear cells from asthmatics. Since reduced IL10 levels may influence the severity of asthma phenotypes, we examined IL10 single-nucleotide polymorphisms (SNPs) for association with asthma severity and allergy phenotypes as quantitative traits. Utilizing DNA samples from 518 Caucasian asthmatic children from the Childhood Asthma Management Program (CAMP) and their parents, we genotyped six IL10 SNPs: 3 in the promoter, 2 in introns, and one in the 3′ UTR. Using family-based association tests, each SNP was tested for association with asthma and allergy phenotypes individually. Population-based association analysis was performed with each SNP locus, the promoter haplotypes and the 6-loci haplotypes. The 3′ UTR SNP was significantly associated with FEV1 as a percent of predicted (FEV1PP) (P=0.0002) in both the family and population analyses. The promoter haplotype GCC was positively associated with IgE levels and FEV1PP (P=0.007 and 0.012, respectively). The promoter haplotype ATA was negatively associated with lnPC20 and FEV1PP (P=0.008 and 0.043, respectively). Polymorphisms in IL10 are associated with asthma phenotypes in this cohort. Further studies of variation in the IL10 gene may help elucidate the mechanism of asthma development in children.
interleukin 10 (IL10); single nucleotide polymorphism (SNP); genetic association; family-based association test (FBAT); haplotype; promoter; 3′; untranslated region (3′UTR)
The adverse effects of corticosteroids on bone mineral accretion (BMA) have been well documented. Vitamin D insufficiency, a prevalent condition in the pediatric population, has also been associated with decreased bone mineral density (BMD).
To determine whether children with asthma who have lower vitamin D levels are more susceptible to the negative effects of corticosteroids on BMD over time.
Children aged 5–12 years with mild-to-moderate asthma who participated in the Childhood Asthma Management Program were followed for a mean of 4.3 years. Total doses of inhaled and oral corticosteroids (OCS) were recorded, serum 25-hydroxyvitamin D3 levels were measured at the beginning of the trial and serial DEXA scans of the lumbar spine were performed. Annual BMA rates were defined as: [(BMD at 4 years follow-up − BMD at baseline)/4 years].
BMA was calculated for 780 subjects. In boys, baseline vitamin D levels significantly modified the relationship between OCS and BMA (vitamin D x OCS interaction, p=0.023). Stratification by vitamin D levels showed a decrease in BMA with increased use of OCS in vitamin D insufficient boys only (p<0.001). Compared to vitamin D sufficient boys, vitamin D insufficient boys exposed to more than 2 courses of oral corticosteroids per year had twice the decrease in BMA rate (relative to boys who were OCS-unexposed).
Vitamin D levels significantly modified the effect of oral corticosteroids on bone mineral accretion in boys. Further research is needed to examine whether vitamin D supplementation in children with poorly controlled asthma may confer benefits to bone health.
Asthma; vitamin D; bone mineral density; corticosteroids
Childhood asthma is a complex disease with known heritability and phenotypic diversity. Although an earlier onset has been associated with more severe disease, there has been no genome-wide association study of the age of onset of asthma in children.
To identify genetic variants associated with earlier onset of childhood asthma.
We conducted the first genome-wide association study (GWAS) of the age of onset of childhood asthma among participants in the Childhood Asthma Management Program (CAMP), and used three independent cohorts from North America, Costa Rica, and Sweden for replication.
Two SNPs were associated with earlier onset of asthma in the combined analysis of CAMP and the replication cohorts: : rs9815663 (Fisher’s P value=2.31 × 10−8) and rs7927044 (P=6.54 × 10−9). Of these two SNPs, rs9815663 was also significantly associated with earlier asthma onset in an analysis including only the replication cohorts. Ten SNPs in linkage disequilibrium with rs9815663 were also associated with earlier asthma onset (2.24 × 10−7 < P < 8.22 ×10−6). Having ≥1 risk allele of the two SNPs of interest (rs9815663 and rs7927044) was associated with lower lung function and higher asthma medication use during 4 years of follow-up in CAMP.
We have identified two SNPs associated with earlier onset of childhood asthma in four independent cohorts.
Asthma; pediatrics; age of onset; asthma genetics; C1orf100; genome-wide association study; pediatric asthma
Rationale: Chronic mucoid Pseudomonas aeruginosa within the airway in cystic fibrosis (CF) patients can determine prognosis. Understanding the risk factors of mucoid P. aeruginosa acquisition may change how we deliver care. This study aims to evaluate whether presence of risk factors reported to predict disease severity including gender, CFTR genotype, bacterial organisms in airway cultures, and serum levels of vitamins A and E, albumin, C-reactive protein, alpha 1-antitrypsin, and immunoglobulins increased the risk of mucoid P. aeruginosa acquisition. Methods: Primary endpoint was age at first transition from negative to positive culture for mucoid P. aeruginosa. Cox proportional hazards regression with time-dependent covariates examined development of mucoid P. aeruginosa infection and its association with longitudinally measured serum biomarkers, pulmonary function, and culture results for other organisms. Results: Median ages at CF diagnosis and at first culture were 0.55 and 5.7 years, respectively. Median number of cultures/patient was 17. Of the 323 subjects, 150 developed mucoid P. aeruginosa during a median 8.1 years’ follow-up. In multivariate analysis, gender (relative hazard [RH] 0.55 for male vs. female, P=0.001), number of DF508 alleles (RH 1.66 for1 or 2 vs. 0, P=0.04), FEV1 % (RH 1.16 for 10% decrease, P=0.008), and most recent Staphylococcus aureus status (RH 0.24 for positive vs. negative, P< 0.0001) remained statistically significant. Conclusion: Female gender, number of DF508 alleles, decreased lung function, and lack of S. aureus on recent sputum culture are important risk factors for early detection of mucoid P. aeruginosa.
cystic fibrosis; Pseudomonas colonization
Rationale: To date, most studies aimed at discovering genetic factors influencing treatment response in asthma have focused on biologic candidate genes. Genome-wide association studies (GWAS) can rapidly identify novel pharmacogenetic loci.
Objectives: To investigate if GWAS can identify novel pharmacogenetic loci in asthma.
Methods: Using phenotypic and GWAS genotype data available through the NHLBI-funded Single-nucleotide polymorphism Health association-Asthma Resource Project, we analyzed differences in FEV1 in response to inhaled corticosteroids in 418 white subjects with asthma. Of the 444,088 single nucleotide polymorphisms (SNPs) analyzed, the lowest 50 SNPs by P value were genotyped in an independent clinical trial population of 407 subjects with asthma.
Measurements and Main Results: The lowest P value for the GWAS analysis was 2.09 × 10−6. Of the 47 SNPs successfully genotyped in the replication population, three were associated under the same genetic model in the same direction, including two of the top four SNPs ranked by P value. Combined P values for these SNPs were 1.06 × 10−5 for rs3127412 and 6.13 × 10−6 for rs6456042. Although these two were not located within a gene, they were tightly correlated with three variants mapping to potentially functional regions within the T gene. After genotyping, each T gene variant was also associated with lung function response to inhaled corticosteroids in each of the trials associated with rs3127412 and rs6456042 in the initial GWAS analysis. On average, there was a twofold to threefold difference in FEV1 response for those subjects homozygous for the wild-type versus mutant alleles for each T gene SNP.
Conclusions: Genome-wide association has identified the T gene as a novel pharmacogenetic locus for inhaled corticosteroid response in asthma.
polymorphism; genome; pharmacogenomics; glucocorticoid
Increasing evidence links altered intestinal flora in infancy to eczema and asthma. No studies have investigated the influence of maternal intestinal flora on wheezing and eczema in early childhood.
To investigate the link between maternal intestinal flora during pregnancy and development of wheeze and eczema in infancy.
Sixty pregnant women from the Boston area gave stool samples during the third trimester of their pregnancy and answered questions during pregnancy about their own health, and about their children’s health when the child was 2 and 6 months of age. Quantitative culture was performed on stool samples and measured in log10colony-forming units(CFU)/gram stool. Primary outcomes included infant wheeze and eczema in the first 6 months of life. Atopic wheeze, defined as wheeze and eczema, was analyzed as a secondary outcome.
In multivariate models adjusted for breastfeeding, daycare attendance and maternal atopy, higher counts of maternal total aerobes (TA) and enterococci (E) were associated with increased risk of infant wheeze (TA: OR 2.32 for 1 log increase in CFU/g stool [95% CI 1.22, 4.42]; E: OR 1.57 [95% CI 1.06, 2.31]). No organisms were associated with either eczema or atopic wheeze.
Conclusions & Clinical Relevance
In our cohort, higher maternal total aerobes and enterococci were related to increased risk of infant wheeze. Maternal intestinal flora may be an important environmental exposure in early immune system development.
infant wheeze; eczema; asthma; microbiota; intestinal flora; maternal flora
The response to treatment for asthma is characterized by wide interindividual variability, with a significant number of patients who have no response. We hypothesized that a genomewide association study would reveal novel pharmacogenetic determinants of the response to inhaled glucocorticoids.
We analyzed a small number of statistically powerful variants selected on the basis of a family-based screening algorithm from among 534,290 single-nucleotide polymorphisms (SNPs) to determine changes in lung function in response to inhaled glucocorticoids. A significant, replicated association was found, and we characterized its functional effects.
We identified a significant pharmacogenetic association at SNP rs37972, replicated in four independent populations totaling 935 persons (P = 0.0007), which maps to the glucocorticoid-induced transcript 1 gene (GLCCI1) and is in complete linkage disequilibrium (i.e., perfectly correlated) with rs37973. Both rs37972 and rs37973 are associated with decrements in GLCCI1 expression. In isolated cell systems, the rs37973 variant is associated with significantly decreased luciferase reporter activity. Pooled data from treatment trials indicate reduced lung function in response to inhaled glucocorticoids in subjects with the variant allele (P = 0.0007 for pooled data). Overall, the mean (± SE) increase in forced expiratory volume in 1 second in the treated subjects who were homozygous for the mutant rs37973 allele was only about one third of that seen in similarly treated subjects who were homozygous for the wild-type allele (3.2 ± 1.6% vs. 9.4 ± 1.1%), and their risk of a poor response was significantly higher (odds ratio, 2.36; 95% confidence interval, 1.27 to 4.41), with genotype accounting for about 6.6% of overall inhaled glucocorticoid response variability.
A functional GLCCI1 variant is associated with substantial decrements in the response to inhaled glucocorticoids in patients with asthma. (Funded by the National Institutes of Health and others; ClinicalTrials.gov number, NCT00000575.)
To assess the feasibility of developing a Combined Clinical and Pharmacogenetic Predictive Test, comprised of multiple single nucleotide polymorphisms (SNPs) that is associated with poor bronchodilator response (BDR).
We genotyped SNPs that tagged the whole genome of the parents and children in the Childhood Asthma Management Program (CAMP) and implemented an algorithm using a family-based association test that ranked SNPs by statistical power. The top eight SNPs that were associated with BDR comprised the Pharmacogenetic Predictive Test. The Clinical Predictive Test was comprised of baseline forced expiratory volume in 1 s (FEV1). We evaluated these predictive tests and a Combined Clinical and Pharmacogenetic Predictive Test in three distinct populations: the children of the CAMP trial and two additional clinical trial populations of asthma. Our outcome measure was poor BDR, defined as BDR of less than 20th percentile in each population. BDR was calculated as the percent difference between the prebronchodilator and postbronchodilator (two puffs of albuterol at 180 μg/puff) FEV1 value. To assess the predictive ability of the test, the corresponding area under the receiver operating characteristic curves (AUROCs) were calculated for each population.
The AUROC values for the Clinical Predictive Test alone were not significantly different from 0.50, the AUROC of a random classifier. Our Combined Clinical and Pharmacogenetic Predictive Test comprised of genetic polymorphisms in addition to FEV1 predicted poor BDR with an AUROC of 0.65 in the CAMP children (n= 422) and 0.60 (n= 475) and 0.63 (n= 235) in the two independent populations. Both the Combined Clinical and Pharmacogenetic Predictive Test and the Pharmacogenetic Predictive Test were significantly more accurate than the Clinical Predictive Test (AUROC between 0.44 and 0.55) in each of the populations.
Our finding that genetic polymorphisms with a clinical trait are associated with BDR suggests that there is promise in using multiple genetic polymorphisms simultaneously to predict which asthmatics are likely to respond poorly to bronchodilators.
asthma; bronchodilator response; personalized medicine; pharmacogenetic test; predictive medicine
Pharmacogenetic testing may change clinical medicine by allowing clinicians to tailor medications based on a patient’s genetic makeup, however, these tests must first be validated in large, real-life populations of subjects that include children. A dearth of knowledge exists for whether pediatric populations are as willing as adult populations to provide samples for such studies.
(1) To assess whether pediatric and adult patients with persistent asthma are willing to provide specimens for DNA extraction and genetic studies. (2) To assess whether patients’ willingness to provide blood as compared to buccal smear specimens differ.
Of 644 patients ages 4–38 years who had three or more prescription fills for inhaled corticosteroids in one year, 60% (385) were randomized to the blood specimen group and 40% (259) were randomized to the buccal smear group in order to study acceptance of different biospecimen collection methods. Research assistants contacted subjects to obtain consent, perform a phone survey, and request a specimen.
There were no baseline differences between subjects randomized to the blood specimen group versus buccal smear group with respect to age, gender, or number of dispensings of inhaled corticosteroids. Of 259 subjects in the buccal smear group, 30% (78) provided samples, and of 385 subjects in the blood specimen group, 16% (60) provided samples. Subjects randomized to the buccal smear group were more likely to provide specimens for genetic study compared to subjects randomized to the blood specimen group (RR 1.21; 95% CI 1.10 – 1.32), even after adjusting for age. Pediatric subjects were more likely to provide specimens for genetic study than adult subjects with 23% (113) of pediatric subjects providing samples and 15% (25) of adult subjects providing samples (p=0.03).
Children with asthma are as likely to participate in genetic studies as adults. Both children and adult subjects are more likely to provide buccal smear specimens rather than blood specimens for genetic study.
pharmacogenetics; patient recruitment; children; buccal swab; blood specimen
Asthma is a complex disease characterized by striking ethnic disparities not explained entirely by environmental, social, cultural, or economic factors. Of the limited genetic studies performed on populations of African descent, notable differences in susceptibility allele frequencies have been observed.
To test the hypothesis that some genes may contribute to the profound disparities in asthma.
We performed a genome-wide association study in two independent populations of African ancestry (935 African American asthma cases and controls from the Baltimore-Washington, D.C. area, and 929 African Caribbean asthmatics and their family members from Barbados) to identify single-nucleotide polymorphisms (SNPs) associated with asthma.
Meta-analysis combining these two African-ancestry populations yielded three SNPs with a combined P-value <10-5 in genes of potential biological relevance to asthma and allergic disease: rs10515807, mapping to alpha-1B-adrenergic receptor (ADRA1B) gene on chromosome 5q33 (3.57×10-6); rs6052761, mapping to prion-related protein (PRNP) on chromosome 20pter-p12 (2.27×10-6); and rs1435879, mapping to dipeptidyl peptidase 10 (DPP10) on chromosome 2q12.3-q14.2. The generalizability of these findings was tested in family and case-control panels of UK and German origin, respectively, but none of the associations observed in the African groups were replicated in these European studies.
Evidence for association was also examined in four additional case-control studies of African Americans; however, none of the SNPs implicated in the discovery population were replicated. This study illustrates the complexity of identifying true associations for a complex and heterogeneous disease such as asthma in admixed populations, especially populations of African descent.
Asthma; GWAS; ADRA1B; PRNP; DPP10; African ancestry; ethnicity; polymorphism; genetic association
Multiple intergenic single-nucleotide polymorphisms (SNPs) near hedgehog interacting protein (HHIP) on chromosome 4q31 have been strongly associated with pulmonary function levels and moderate-to-severe chronic obstructive pulmonary disease (COPD). However, whether the effects of variants in this region are related to HHIP or another gene has not been proven. We confirmed genetic association of SNPs in the 4q31 COPD genome-wide association study (GWAS) region in a Polish cohort containing severe COPD cases and healthy smoking controls (P = 0.001 to 0.002). We found that HHIP expression at both mRNA and protein levels is reduced in COPD lung tissues. We identified a genomic region located ∼85 kb upstream of HHIP which contains a subset of associated SNPs, interacts with the HHIP promoter through a chromatin loop and functions as an HHIP enhancer. The COPD risk haplotype of two SNPs within this enhancer region (rs6537296A and rs1542725C) was associated with statistically significant reductions in HHIP promoter activity. Moreover, rs1542725 demonstrates differential binding to the transcription factor Sp3; the COPD-associated allele exhibits increased Sp3 binding, which is consistent with Sp3's usual function as a transcriptional repressor. Thus, increased Sp3 binding at a functional SNP within the chromosome 4q31 COPD GWAS locus leads to reduced HHIP expression and increased susceptibility to COPD through distal transcriptional regulation. Together, our findings reveal one mechanism through which SNPs upstream of the HHIP gene modulate the expression of HHIP and functionally implicate reduced HHIP gene expression in the pathogenesis of COPD.
Atopy and plasma IgE concentration are genetically complex traits, and the specific genetic risk factors that lead to IgE dysregulation and clinical atopy are an area of active investigation.
To ascertain the genetic risk factors which lead to IgE dysregulation.
A genome wide association study (GWAS) was performed in 6,819 participants from the Framingham Heart Study (FHS). Seventy of the top SNPs were selected based on p-values and linkage disequilibrium among neighboring SNPs and evaluated in a meta-analysis with five independent populations from the KORA, B58C, and CAMP cohorts.
Thirteen SNPs located in the region of three genes, FCER1A, STAT6, and IL-13, were found to have genome-wide significance in the FHS GWAS. The most significant SNPs from the three regions were rs2251746 (FCER1A, p-value 2.11×10-12), rs1059513 (STAT6, p-value 2.87×10-08), and rs1295686 (IL-13, p-value 3.55×10-08). Four additional gene regions - HLA-G, HLA-DQA2, HLA-A, and DARC - reached genome-wide statistical significance in meta-analysis combining FHS and replication cohorts, although the DARC association did not appear independent of SNPs in the nearby FCER1A gene.
This GWAS of the FHS has identified genetic loci in HLA genes that may have a role in the pathogenesis of IgE dysregulation and atopy. It also confirmed the association of known susceptibility loci, FCER1A, STAT6, and IL-13, for the dysregulation of total IgE.
total IgE; atopy; asthma; GWAS
Distinct receptors likely exist for leukotriene(LT)E4, a potent mediator of airway inflammation. Purinergic receptor P2Y12 is needed for LTE4-induced airways inflammation, and P2Y12 antagonism attenuates house dust mite-induced pulmonary eosinophilia in mice. Although experimental data support a role for P2Y12 in airway inflammation, its role in human asthma has never been studied.
To test for association between variants in the P2Y12 gene (P2RY12) and lung function in human subjects with asthma, and to examine for gene-by-environment interaction with house dust mite exposure.
19 single nucleotide polymorphisms (SNPs) in P2RY12 were genotyped in 422 children with asthma and their parents (n=1266). Using family-based methods, we tested for associations between these SNPs and five lung function measures. We performed haplotype association analyses and tested for gene-by-environment interactions using house dust mite exposure. We used the false discovery rate to account for multiple comparisons.
Five SNPs in P2RY12 were associated with multiple lung function measures (P values 0.006–0.025). Haplotypes in P2RY12 were also associated with lung function (P values 0.0055–0.046). House dust mite exposure modulated associations between P2RY12 and lung function, with minor allele homozygotes exposed to house dust mite demonstrating worse lung function than those unexposed (significant interaction P values 0.0028–0.040).
Conclusions and clinical relevance
P2RY12 variants were associated with lung function in a large family-based asthma cohort. House dust mite exposure caused significant gene-by-environment effects. Our findings add the first human evidence to experimental data supporting a role for P2Y12 in lung function. P2Y12 could represent a novel target for asthma treatment.
Purinergic receptor; leukotriene; asthma; house dust mite; lung function
It has recently been shown that vitamin D deficiency can increase asthma development and severity and that variations in vitamin D receptor genes are associated with asthma susceptibility.
We sought to find genetic factors that might interact with vitamin D levels to affect the risk of asthma exacerbation. Methods: We conducted a genome-wide study of gene–vitamin D interaction on asthma exacerbations using population-based and family-based approaches on 403 subjects and trios from the Childhood Asthma Management Program. Twenty-three polymorphisms with significant interactions were studied in a replication analysis in 584 children from a Costa Rican cohort. Results: We identified 3 common variants in the class I MHC–restricted T cell–associated molecule gene (CRTAM) that were associated with an increased rate of asthma exacerbations based on the presence of a low circulating vitamin D level. These results were replicated in a second independent population (unadjusted combined interaction, P =.00028–.00097; combined odds ratio, 3.28–5.38). One variant, rs2272094, is a nonsynonymous coding polymorphism of CRTAM. Functional studies on cell lines confirmed the interaction of vitamin D and rs2272094 on CRTAM expression. CRTAM is highly expressed in activated human CD8+ and natural killer T cells, both of which have been implicated in asthmatic patients.
The findings highlight an important gene-environment interaction that elucidates the role of vitamin D and CD8+ and natural killer T cells in asthma exacerbation in a genome-wide gene-environment interaction study that has been replicated in an independent population. The results suggest the potential importance of maintaining adequate vitamin D levels in subsets of high-risk asthmatic patients.
Gene-environment interaction; genome-wide association study; vitamin D; asthma exacerbation
The Commonwealth of Massachusetts increased the copayment for prescription drugs by $1.50 for Medicaid (MassHealth) beneficiaries in 2003. We sought to determine the likely health outcomes and cost shifts attributable to this copayment increase using the example of inhaled corticosteroids (ICS) use among adult asthmatic Medicaid beneficiaries.
We compared the predicted costs and health outcomes projected over a 1-year time horizon with and without the increase in copayment from the perspective of MassHealth, providers, pharmacies, and MassHealth beneficiaries by employing decision analysis simulation model.
In a target population of 17,500 adult asthmatics, increased copayments from 50¢ to $2.00 would result in an additional 646 acute events per year, caused by increased drug nonadherence. Annual combined net savings for the state and federal governments would be $2.10 million. Projected MassHealth savings are attributable to both decreased drug utilization and lower pharmacy reimbursement rates; these more than offset the additional costs of more frequent acute exacerbations. Pharmacies would lose $1.98 million in net revenues, MassHealth beneficiaries would pay an additional $0.28 million, and providers would receive additional $0.16 million.
Over its first year of implementation, increase in the prescription drug copayment is expected to produce more frequent acute exacerbations among asthmatic MassHealth beneficiaries who use ICS and to shift the financial burden from government to other stakeholders.
asthma; copayment; medicaid; prescription drug
Omalizumab (trade name Xolair) is approved by the US Food and Drug Administration for treatment of moderate-to-severe allergic asthma. Given the high acquisition cost of omalizumab, its role and cost-effectiveness in disease management require definition.
We sought to identify the clinical and economic circumstances under which omalizumab might or might not be a cost-effective option by using a mathematic model.
We merged published data on clinical and economic outcomes (including acute event incidence, frequency/severity of hospitalizations, and health-related quality of life) to project 10-year costs, quality-adjusted life years (QALYs), and cost-effectiveness of treatment with omalizumab in addition to inhaled corticosteroids. Sensitivity analyses were conducted by using input data ranges from a variety of sources (published clinical trials and observational databases).
For patients with baseline acute event rates, omalizumab conferred an additional 1.7 quality-adjusted months at an incremental cost of $131,000 over a 10-year planning horizon, implying a cost-effectiveness ratio of $821,000 per QALY gained. For patients with 5 times the baseline acute event rate, the cost-effectiveness ratio was $491,000 per QALY gained. The projected cost-effectiveness ratio could fall within a range of other programs that are widely considered to be cost-effective if the cost of omalizumab decreases to less than $200.
Omalizumab is not cost-effective for most patients with severe asthma. The projected cost-effectiveness ratios could fall within a favorable range if the cost of omalizumab decreases significantly.
Based on the high cost of omalizumab, it is especially important that clinicians explore alternative medications for asthma before initiating omalizumab.
Omalizumab; cost-effectiveness; asthma; anti-IgE
Patients with severe or difficult-to-treat asthma account for substantial asthma morbidity, mortality, and healthcare burden despite comprising only a small proportion of the total asthma population. TENOR, a multicenter, observational, prospective cohort study was initiated in 2001. It enrolled 4,756 adults, adolescents and children with severe or difficult-to-treat asthma who were followed semi-annually and annually for three years, enabling insight to be gained into this understudied population. A broad range of demographic, clinical, and patient self-reported assessments were completed during the follow-up period. Here, we present key findings from the TENOR registry in relation to asthma control and exacerbations, including the identification of specific subgroups found to be at particularly high-risk. Identification of the factors and subgroups associated with poor asthma control and increased risk of exacerbations can help physicians design individual asthma management, and improve asthma-related health outcomes for these patients.
Severe asthma; Difficult-to-treat asthma; Asthma control; Exacerbation
Bronchodilator response (BDR) is an important asthma phenotype that measures reversibility of airway obstruction by comparing lung function (i.e. FEV1) before and after the administration of a short-acting β2-agonist, the most common rescue medications used for the treatment of asthma. BDR also serves as a test of β2-agonist efficacy. BDR is a complex trait that is partly under genetic control. A genome-wide association study (GWAS) of BDR, quantified as percent change in baseline FEV1 after administration of a β2-agonist, was performed with 1,644 non-Hispanic white asthmatic subjects from six drug clinical trials: CAMP, LOCCS, LODO, a medication trial conducted by Sepracor, CARE, and ACRN. Data for 469,884 single-nucleotide polymorphisms (SNPs) were used to measure the association of SNPs with BDR using a linear regression model, while adjusting for age, sex, and height. Replication of primary P-values was attempted in 501 white subjects from SARP and 550 white subjects from DAG. Experimental evidence supporting the top gene was obtained via siRNA knockdown and Western blotting analyses. The lowest overall combined P-value was 9.7E-07 for SNP rs295137, near the SPATS2L gene. Among subjects in the primary analysis, those with rs295137 TT genotype had a median BDR of 16.0 (IQR = [6.2, 32.4]), while those with CC or TC genotypes had a median BDR of 10.9 (IQR = [5.0, 22.2]). SPATS2L mRNA knockdown resulted in increased β2-adrenergic receptor levels. Our results suggest that SPATS2L may be an important regulator of β2-adrenergic receptor down-regulation and that there is promise in gaining a better understanding of the biological mechanisms of differential response to β2-agonists through GWAS.
Bronchodilator response (BDR) is an important asthma phenotype that measures reversibility of airway obstruction by comparing lung function before and after the administration of short-acting β2-agonists, common medications used for asthma treatment. We performed a genome-wide association study of BDR with 1,644 white asthmatic subjects from six drug clinical trials and attempted to replicate these findings in 1,051 white subjects from two independent cohorts. The most significant associated variant was near the SPATS2L gene. We knocked down SPATS2L mRNA in human airway smooth muscle cells and found that β2-adrenergic receptor levels increased, suggesting that SPATS2L may be a regulator of BDR. Our results highlight the promise of pursuing GWAS results that do not necessarily reach genome-wide significance and are an example of how results from pharmacogenetic GWAS can be studied functionally.
Thymic stromal lymphopoietin (TSLP) triggers dendritic cell–mediated T helper (Th) 2 inflammatory responses. A single-nucleotide polymorphism (SNP), rs3806933, in the promoter region of the TSLP gene creates a binding site for the transcription factor activating protein (AP)–1. The variant enhances AP-1 binding to the regulatory element, and increases the promoter–reporter activity of TSLP in response to polyinosinic-polycytidylic acid (poly[I:C]) stimulation in normal human bronchial epithelium (NHBE). We investigated whether polymorphisms including the SNP rs3806933 could affect the susceptibility to and clinical phenotypes of bronchial asthma. We selected three representative (i.e., Tag) SNPs and conducted association studies of the TSLP gene, using two independent populations (639 patients with childhood atopic asthma and 838 control subjects, and 641 patients with adult asthma and 376 control subjects, respectively). We further examined the effects of corticosteroids and a long-acting β2-agonist (salmeterol) on the expression levels of the TSLP gene in response to poly(I:C) in NHBE. We found that the promoter polymorphisms rs3806933 and rs2289276 were significantly associated with disease susceptibility in both childhood atopic and adult asthma. The functional SNP rs3806933 was associated with asthma (meta-analysis, P = 0.000056; odds ratio, 1.29; 95% confidence interval, 1.14–1.47). A genotype of rs2289278 was correlated with pulmonary function. Moreover, the induction of TSLP mRNA and protein expression induced by poly(I:C) in NHBE was synergistically impaired by a corticosteroid and salmeterol. TSLP variants are significantly associated with bronchial asthma and pulmonary function. Thus, TSLP may serve as a therapeutic target molecule for combination therapy.
asthma; TSLP; bronchial epithelial cells; combination therapy; genetic polymorphisms
For genomewide association studies with family-based designs, we propose a Bayesian approach. We show that standard TDT/FBAT statistics can naturally be implemented in a Bayesian framework. We construct a Bayes factor conditional on the offspring phenotype and parental genotype data and then use the data we conditioned on to inform the prior odds for each marker. In the construction of the prior odds, the evidence for association for each single marker is obtained at the population-level by estimating the genetic effect size in the conditional mean model. Since such genetic effect size estimates are statistically independent of the effect size estimation within the families, the actual data set can inform the construction of the prior odds without any statistical penalty. In contrast to Bayesian approaches that have recently been proposed for genomewide association studies, our approach does not require assumptions about the genetic effect size; this makes the proposed method entirely data-driven. The power of the approach was assessed through simulation. We then applied the approach to a genomewide association scan to search for associations between single nucleotide polymorphisms and body mass index in the Childhood Asthma Management Program data.
family-based association tests; Bayes factors; complex traits