Corticosteroids exert their anti-inflammatory action by binding and activating the intracellular the glucocorticoid receptor (GR) hetero-complex.
Evaluate the genes HSPCB, HSPCA, STIP1, HSPA8, DNAJB1, PTGES3, FKBP5, and FKBP4 on corticosteroid response.
Caucasian asthmatics (382) randomized to once daily flunisolide or conventional inhaled corticosteroid therapy were genotyped. Outcome measures were baseline FEV1, % predicted FEV1, and % change in FEV1 after corticosteroid treatment. Multivariable analyses adjusted for age, gender, and height, were performed fitting the most appropriate genetic model based on quantitative mean derived from ANOVA models to determine if there was an independent effect of polymorphisms on change in FEV1 independent of baseline level.
Positive recessive model correlations for STIP1 SNPs were observed for baseline FEV1 [rs4980524, p=0.009; rs6591838, p=0.0045; rs2236647, p=0.002; and rs2236648; p=0.013], baseline % predicted FEV1 [rs4980524, p=0.002; rs6591838, p=0.017; rs2236647, p=0.003; and rs2236648; p=0.008] ; % change in FEV1 at 4 weeks [rs4980524, p=0.044; rs6591838, p=0.016; rs2236647; p=0.01] and 8 weeks therapy [rs4980524, p=0.044; rs6591838, p=0.016; rs2236647; p=0.01]. Haplotypic associations were observed for baseline FEV1 and % change in FEV1 at 4 weeks therapy [p=0.05 and p=0.01, respectively]. Significant trends towards association were observed for baseline % predicted FEV1 and % change in FEV1 at 8 weeks therapy. Positive correlations between haplotypes and % change in FEV1 were also observed.
STIP1 genetic variations may play a role in regulating corticosteroid response in asthmatics with reduced lung function. Replication in a second asthma population is required to confirm these observations.
Identifying genes that regulate corticosteroid responses could allow a priori determination of individual responses to corticosteroid therapy, leading to more effective dosing and/or selection of drug therapies for treating asthma.
corticosteroid; pharmacogenetics; glucocorticoid receptor; SNP; heat shock protein; heat shock organizing protein; immunophilin
The bronchodilator response (BDR) reflects the reversibility of airflow obstruction and is recommended as an adjunctive test to diagnose asthma. The validity of the commonly used definition of BDR, a 12% or greater change in FEV1 from baseline, has been questioned in childhood.
We sought to examine the diagnostic accuracy of the BDR test by using 3 large pediatric cohorts.
Cases include 1041 children with mild-to-moderate asthma from the Childhood Asthma Management Program.
Control subjects (nonasthmatic and nonwheezing) were chosen from Project Viva and Home Allergens, 2 population-based pediatric cohorts. Receiver operating characteristic curves were constructed, and areas under the curve were calculated for different BDR cutoffs.
A total of 1041 cases (59.7% male; mean age, 8.9 ± 2.1 years) and 250 control subjects (46.8% male; mean age, 8.7 ± 1.7 years) were analyzed, with mean BDRs of 10.7% ± 10.2% and 2.7% ± 8.4%, respectively. The BDR test differentiated asthmatic patients from nonasthmatic patients with a moderate accuracy (area under the curve, 73.3%).
Despite good specificity, a cutoff of 12% was associated with poor sensitivity (35.6%). A cutoff of less than 8% performed significantly better than a cutoff of 12% (P = .03, 8% vs 12%).
Our findings highlight the poor sensitivity associated with the commonly used 12% cutoff for BDR. Although our data show that a threshold of less than 8% performs better than 12%, given the variability of this test in children, we conclude that it might be not be appropriate to choose a specific BDR cutoff as a criterion for the diagnosis of asthma.
Asthma; bronchodilator response; diagnosis
Reversibility of airway obstruction in response to β2-agonists is highly variable among asthmatics, which is partially attributed to genetic factors. In a genome-wide association study of acute bronchodilator response (BDR) to inhaled albuterol, 534,290 single nucleotide polymorphisms (SNPs) were tested in 403 white trios from the Childhood Asthma Management Program using five statistical models to determine the most robust genetic associations. The primary replication phase included 1397 polymorphisms in three asthma trials (pooled n=764). The second replication phase tested 13 SNPs in three additional asthma populations (n=241, n=215, and n=592). An intergenic SNP on chromosome 10, rs11252394, proximal to several excellent biological candidates, significantly replicated (p=1.98×10−7) in the primary replication trials. An intronic SNP (rs6988229) in the collagen (COL22A1) locus also provided strong replication signals (p=8.51×10−6). This study applied a robust approach for testing the genetic basis of BDR and identified novel loci associated with this drug response in asthmatics.
pharmacogenetics; asthma; bronchodilator response; genome-wide association study; albuterol
Asthma is a chronic inflammatory respiratory disease that affects over 300 million people worldwide. Glucocorticoids are a mainstay therapy for asthma because they exert anti-inflammatory effects in multiple lung tissues, including the airway smooth muscle (ASM). However, the mechanism by which glucocorticoids suppress inflammation in ASM remains poorly understood. Using RNA-Seq, a high-throughput sequencing method, we characterized transcriptomic changes in four primary human ASM cell lines that were treated with dexamethasone—a potent synthetic glucocorticoid (1 µM for 18 hours). Based on a Benjamini-Hochberg corrected p-value <0.05, we identified 316 differentially expressed genes, including both well known (DUSP1, KLF15, PER1, TSC22D3) and less investigated (C7, CCDC69, CRISPLD2) glucocorticoid-responsive genes. CRISPLD2, which encodes a secreted protein previously implicated in lung development and endotoxin regulation, was found to have SNPs that were moderately associated with inhaled corticosteroid resistance and bronchodilator response among asthma patients in two previously conducted genome-wide association studies. Quantitative RT-PCR and Western blotting showed that dexamethasone treatment significantly increased CRISPLD2 mRNA and protein expression in ASM cells. CRISPLD2 expression was also induced by the inflammatory cytokine IL1β, and small interfering RNA-mediated knockdown of CRISPLD2 further increased IL1β-induced expression of IL6 and IL8. Our findings offer a comprehensive view of the effect of a glucocorticoid on the ASM transcriptome and identify CRISPLD2 as an asthma pharmacogenetics candidate gene that regulates anti-inflammatory effects of glucocorticoids in the ASM.
African Americans have worse asthma outcomes compared to whites. Adrenoceptor beta 2, surface gene (ADRB2) Gly16Arg genotypes have been associated with β2-agonist bronchodilator response, asthma exacerbation rate, response to methacholine, and lung function decline but not specifically in African Americans.
We sought to compare the provocative concentration of methacholine that causes a 20% fall in FEV1 (PC20) in African Americans and whites with asthma who were ADRB2 homozygous at codon16 (Arg16Arg or Gly16Gly).
African Americans and whites whose parents and grandparents were of the same race, aged ≥ 10 years, with baseline FEV1 of ≥60% predicted, and no upper or lower respiratory tract infection within the previous 2 weeks meeting genotype criteria were enrolled. PC20 was measured after withholding short-acting and long-acting β2-agonists for 8 and 12 hours respectively, montelukast for 24 hours, ipratropium bromide and inhaled corticosteroids for 12 hours, and antihistamines for 72 hours.
423 participants were screened and 88 had a positive challenge. Participants were 32yrs ± 19yrs (mean ± SD), 70% female, 51% White (vs. African American), 6% Hispanic. Similar numbers of participants were using inhaled corticosteroids by race and genotype. There were significant differences in log PC20 between race/genotype groups (p=0.012). African American Arg16Arg participants had a lower log PC20 than White Gly16Gly (p=0.009) and African American Gly16Gly (p=0.041) participants. Both race and genotype contributed significantly to the model (p=0.037 and p=0.014, respectively) but there was no interaction between race and genotype on log PC20.
CONCLUSIONS AND CLINICAL RELEVANCE
Airway hyperresponsiveness is influenced by race and the ADRB2 codon 16 polymorphism. African Americans with the Arg16Arg genotype have increased airway reactivity and may be at risk for worse asthma outcomes. Inclusion of genetic information as an additional clinical tool may aid in the personalization of asthma management decisions.
[ClinicalTrials.gov Identifier: NCT00708227].
African American; White; methacholine; asthma; ADRB2; genotype
While accurate measures of heritability are needed to understand the pharmacogenetic basis of drug treatment response, these are generally not available, since it is unfeasible to give medications to individuals for which treatment is not indicated. Using a polygenic linear mixed modeling approach, we estimated lower-bounds on asthma heritability and the heritability of two related drug-response phenotypes, bronchodilator response and airway hyperreactivity, using genome-wide SNP data from existing asthma cohorts. Our estimate of the heritability for bronchodilator response is 28.5% (se 16%, p = 0.043) and airway hyperresponsiveness is 51.1% (se 34%, p = 0.064), while we estimate asthma genetic liability at 61.5% (se 16%, p < 0.001). Our results agree with previously published estimates of the heritability of these traits, suggesting that the LMM method is useful for computing the heritability of other pharmacogenetic traits. Furthermore, our results indicate that multiple SNP main-effects, including SNPs as yet unidentified by GWAS methods, together explain a sizable portion of the heritability of these traits.
Asthma; Pharmacogenetics; Heritability; Bronchodilator Response; Airway Hyperresponsiveness
Genome-wide association studies (GWAS) have emerged as a powerful tool to identify loci that affect drug response or susceptibility to adverse drug reactions. However, current GWAS based on a simple analysis of associations between genotype and phenotype ignores the biochemical reactions of drug response, thus limiting the scope of inference about its genetic architecture. To facilitate the inference of GWAS in pharmacogenomics, we sought to undertake the mathematical integration of the pharmacodynamic process of drug reactions through computational models. By estimating and testing the genetic control of pharmacodynamic and pharmacokinetic parameters, this mechanistic approach does not only enhance the biological and clinical relevance of significant genetic associations, but also improve the statistical power and robustness of gene detection. This report discusses the general principle and development of pharmacodynamics-based GWAS, highlights the practical use of this approach in addressing various pharmacogenomic problems, and suggests that this approach will be an important method to study the genetic architecture of drug responses or reactions.
Inhaled corticosteroids are the most commonly used controller medications prescribed for asthma. Two single-nucleotide polymorphisms (SNPs), rs1876828 in CRHR1 and rs37973 in GLCCI1, have previously been associated with corticosteroid efficacy. We studied data from four existing clinical trials of asthmatics who received inhaled corticosteroids and had lung function measured by forced expiratory volume in one second (FEV1) before and after the period of such treatment. We combined the two SNPs rs37973 and rs1876828 into a predictive test of FEV1 change using a Bayesian model, which identified patients with good or poor steroid response (highest or lowest quartile, respectively) with predictive performance of 65.7% (p = 0.039 vs. random) area under the receiver-operator characteristic curve in the training population and 65.9% (p = 0.025 vs. random) in the test population. These findings show that two genetic variants can be combined into a predictive test that achieves similar accuracy and superior replicability compared with single SNP predictors.
Pharmacogenetics; Asthma; Glucocorticoids; Predictive Modeling
Poor maternal vitamin D intake is a risk factor for subsequent childhood asthma, suggesting that in utero changes related to vitamin D responsive genes might play a crucial role in later disease susceptibility. We hypothesized that vitamin D pathway genes are developmentally active in the fetal lung and that these developmental genes would be associated with asthma susceptibility and regulation in asthma.
Vitamin D pathway genes were derived from PubMed and Gene Ontology surveys. Principal component analysis was used to identify characteristic lung development genes.
Vitamin D regulated genes were markedly over-represented in normal human (odds ratio OR 2.15, 95% confidence interval CI: 1.69-2.74) and mouse (OR 2.68, 95% CI: 2.12-3.39) developing lung transcriptomes. 38 vitamin D pathway genes were in both developing lung transcriptomes with >63% of genes more highly expressed in the later than earlier stages of development. In immortalized B-cells derived from 95 asthmatics and their unaffected siblings, 12 of the 38 (31.6%) vitamin D pathway lung development genes were significantly differentially expressed (OR 3.00, 95% CI: 1.43-6.21), whereas 11 (29%) genes were significantly differentially expressed in 43 control versus vitamin D treated immortalized B-cells from Childhood Asthma Management Program subjects (OR 2.62, 95% CI: 1.22-5.50). 4 genes, LAMP3, PIP5K1B, SCARB2 and TXNIP were identified in both groups; each displays significant biologic plausibility for a role in asthma.
Our findings demonstrate a significant association between early lung development and asthma–related phenotypes for vitamin D pathway genes, supporting a genomic mechanistic basis for the epidemiologic observations relating maternal vitamin D intake and childhood asthma susceptibility.
Vitamin D; Cholecalciferol; Lung development; Asthma; Fetal programming
A pro-asthmatic culture milieu and β2-agonist (isoproterenol) were previously shown to regulate the expression of select transcription factors (TFs) within human airway epithelial and smooth muscle cells. This study tests 1116 single nucleotide polymorphisms (SNPs) across 98 of these TF genes for association with bronchodilator response (BDR) in asthma patients. Genotyping was conducted using the Illumina HumanHap550v3 Beadchip in 403 non-Hispanic White asthmatic children and their parents. SNPs were evaluated for association with BDR using family and population-based analyses. Forty-two SNPs providing p values < 0.1 in both analyses were then genotyped in three adult asthma trials. One SNP 5’ of the thyroid hormone receptor beta gene was associated with BDR in the childhood population and two adult populations (p value = 0.0012). This investigation identified a novel locus for inter-individual variability in BDR and represents a translation of a cellular drug-response study to potential personalization of clinical asthma management.
Bronchodilator response; transcription factor; association; thyroid hormone receptor β; asthma; pharmacogenetics
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
Airway hyperresponsiveness (AHR), a primary characteristic of asthma, involves increased airway smooth muscle contractility in response to certain exposures. We sought to determine whether common genetic variants were associated with AHR severity.
A genome-wide association study (GWAS) of AHR, quantified as the natural log of the dosage of methacholine causing a 20% drop in FEV1, was performed with 994 non-Hispanic white asthmatic subjects from three drug clinical trials: CAMP, CARE, and ACRN. Genotyping was performed on Affymetrix 6.0 arrays, and imputed data based on HapMap Phase 2, was used to measure the association of SNPs with AHR using a linear regression model. Replication of primary findings was attempted in 650 white subjects from DAG, and 3,354 white subjects from LHS. Evidence that the top SNPs were eQTL of their respective genes was sought using expression data available for 419 white CAMP subjects.
The top primary GWAS associations were in rs848788 (P-value 7.2E-07) and rs6731443 (P-value 2.5E-06), located within the ITGB5 and AGFG1 genes, respectively. The AGFG1 result replicated at a nominally significant level in one independent population (LHS P-value 0.012), and the SNP had a nominally significant unadjusted P-value (0.0067) for being an eQTL of AGFG1.
Based on current knowledge of ITGB5 and AGFG1, our results suggest that variants within these genes may be involved in modulating AHR. Future functional studies are required to confirm that our associations represent true biologically significant findings.
Asthma; Airway hyperresponsiveness; Genome-wide association study; ITGB5; AGFG1
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
In order to improve therapeutic outcomes, there is a tremendous need to identify patients who are likely to respond to a given asthma treatment. Pharmacogenomic studies have explained a portion of the variability in drug response and provided an increasing list of candidate genes and SNPs. However, as phenotypic variation arises from a network of complex interactions among genetic and environmental factors, rather than individual genes or SNPs, a multidisciplinary, systems-level approach is required in order to understand the inter-relationships among these factors. Systems biology, which seeks to capture interactions between genetic factors and other variables, offers a promising approach to improved therapeutic outcomes in asthma. This aritcle will review and update progress in the pharmacogenomics of asthma and then discuss the application of systems biology approaches to asthma pharmacogenomics.
asthma; genes; GWAS; network medicine; pharmacogenomics; SNP; systems biology
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
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.)
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 common chronic respiratory disease characterized by airway hyperresponsiveness (AHR). The genetics of asthma have been widely studied in mouse and human, and homologous genomic regions have been associated with mouse AHR and human asthma-related phenotypes. Our goal was to identify asthma-related genes by integrating AHR associations in mouse with human genome-wide association study (GWAS) data. We used Efficient Mixed Model Association (EMMA) analysis to conduct a GWAS of baseline AHR measures from males and females of 31 mouse strains. Genes near or containing SNPs with EMMA p-values <0.001 were selected for further study in human GWAS. The results of the previously reported EVE consortium asthma GWAS meta-analysis consisting of 12,958 diverse North American subjects from 9 study centers were used to select a subset of homologous genes with evidence of association with asthma in humans. Following validation attempts in three human asthma GWAS (i.e., Sepracor/LOCCS/LODO/Illumina, GABRIEL, DAG) and two human AHR GWAS (i.e., SHARP, DAG), the Kv channel interacting protein 4 (KCNIP4) gene was identified as nominally associated with both asthma and AHR at a gene- and SNP-level. In EVE, the smallest KCNIP4 association was at rs6833065 (P-value 2.9e-04), while the strongest associations for Sepracor/LOCCS/LODO/Illumina, GABRIEL, DAG were 1.5e-03, 1.0e-03, 3.1e-03 at rs7664617, rs4697177, rs4696975, respectively. At a SNP level, the strongest association across all asthma GWAS was at rs4697177 (P-value 1.1e-04). The smallest P-values for association with AHR were 2.3e-03 at rs11947661 in SHARP and 2.1e-03 at rs402802 in DAG. Functional studies are required to validate the potential involvement of KCNIP4 in modulating asthma susceptibility and/or AHR. Our results suggest that a useful approach to identify genes associated with human asthma is to leverage mouse AHR association data.
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
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
Despite the availability of several classes of asthma medications and their overall effectiveness, a significant portion of patients fail to respond to these therapeutic agents. Evidence suggests that genetic factors may partly mediate the heterogeneity in asthma treatment response. This review discusses important findings in asthma pharmacogenetics and pharmacogenomics studies conducted to date, examines limitations of these studies and finally, proposes future research directions in this field. The focus will be on the three major classes of asthma medications: β-adrenergic receptor agonists, inhaled corticosteroids and leukotriene modifiers. Although many studies are limited by small sample sizes and replication of the findings is needed, several candidate genes have been identified. High-throughput technologies is also allowing for large-scale genetic investigations. Thus, the future is promising for a personalized treatment of asthma, which will improve therapeutic outcomes, minimize side effects and lead to a more cost-effective care.
asthma; pharmacogenetics; pharmacogenomics
Rationale: β2-agonists, the most common treatment for asthma, have a wide interindividual variability in response, which is partially attributed to genetic factors. We previously identified single nucleotide polymorphisms in the arginase 1 (ARG1) gene, which are associated with β2-agonist bronchodilator response (BDR).
Objectives: To identify cis-acting haplotypes in the ARG1 locus that are associated with BDR in patients with asthma and regulate gene expression in vitro.
Methods: We resequenced ARG1 in 96 individuals and identified three common, 5′ haplotypes (denoted 1, 2, and 3). A haplotype-based association analysis of BDR was performed in three independent, adult asthma drug trial populations. Next, each haplotype was cloned into vectors containing a luciferase reporter gene and transfected into human airway epithelial cells (BEAS-2B) to ascertain its effect on gene expression.
Measurements and Main Results: BDR varied by haplotype in each of the three populations with asthma. Individuals with haplotype 1 were more likely to have higher BDR, compared to those with haplotypes 2 and 3, which is supported by odds ratios of 1.25 (95% confidence interval, 1.03–1.71) and 2.18 (95% confidence interval, 1.34–2.52), respectively. Luciferase expression was 50% greater in cells transfected with haplotype 1 compared to haplotypes 2 and 3.
Conclusions: The identified ARG1 haplotypes seem to alter BDR and differentially regulate gene expression with a concordance of decreased BDR and reporter activity from haplotypes 2 and 3. These findings may facilitate pharmacogenetic tests to predict individuals who may benefit from other therapeutic agents in addition to β2-agonists for optimal asthma management.
Clinical trial registered with www.clinicaltrials.gov (NCT00156819, NCT00046644, and NCT00073840).
pharmacogenetics; asthma; β2-agonist