Importin-13 genetic variation is associated with improved airway responsiveness in childhood asthma

doi:10.1186/1465-9921-10-67

Respir Res. 2009; 10(1): 67.

Published online 2009 July 20. doi: 10.1186/1465-9921-10-67

PMCID: PMC2724419

Importin-13 genetic variation is associated with improved airway responsiveness in childhood asthma

Benjamin A Raby,^1,^2,^3,⁴ Kristel Van Steen,⁵ Jessica Lasky-Su,^1,⁴ Kelan Tantisira,^1,^2,^3,⁴ Feige Kaplan,⁶ and Scott T Weiss^1,^2,^3,⁴

¹Channing Laboratory, Department of Medicine, Brigham and Women's Hospital, Boston, Massachusetts, USA

²Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital, Boston, Massachusetts, USA

³Harvard Medical School, Boston, Massachusetts, USA

⁴Center for Genomic Medicine, Brigham and Women's Hospital, Boston Massachusetts, USA

⁵Department of Oto-rhino-laryngology & Department of Applied Mathematics and Computer Science, University of Ghent, Belgium

⁶Departments of Human Genetics and Pediatrics, McGill University, Montreal Quebec, Canada

Corresponding author.

Benjamin A Raby: rebar/at/channing.harvard.edu; Kristel Van Steen: kvansteen/at/gmail.com; Jessica Lasky-Su: jessica.a.su/at/gmail.com; Kelan Tantisira: rekgt/at/channing.harvard.edu; Feige Kaplan: feige.kaplan/at/mcgill.ca; Scott T Weiss: restw/at/channing.harvard.edu

Received March 31, 2009; Accepted July 20, 2009.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article has been cited by other articles in PMC.

Abstract

Background

Glucocorticoid function is dependent on efficient translocation of the glucocorticoid receptor (GR) from the cytoplasm to the nucleus of cells. Importin-13 (IPO13) is a nuclear transport receptor that mediates nuclear entry of GR. In airway epithelial cells, inhibition of IPO13 expression prevents nuclear entry of GR and abrogates anti-inflammatory effects of glucocorticoids. Impaired nuclear entry of GR has been documented in steroid-non-responsive asthmatics. We hypothesize that common IPO13 genetic variation influences the anti-inflammatory effects of inhaled corticosteroids for the treatment of asthma, as measured by change in methacholine airway hyperresponsiveness (AHR-PC₂₀).

Methods

10 polymorphisms were evaluated in 654 children with mild-to-moderate asthma participating in the Childhood Asthma Management Program (CAMP), a clinical trial of inhaled anti-inflammatory medications (budesonide and nedocromil). Population-based association tests with repeated measures of PC₂₀were performed using mixed models and confirmed using family-based tests of association.

Results

Among participants randomized to placebo or nedocromil, IPO13 polymorphisms were associated with improved PC₂₀(i.e. less AHR), with subjects harboring minor alleles demonstrating an average 1.51–2.17 fold increase in mean PC₂₀at 8-months post-randomization that persisted over four years of observation (p = 0.01–0.005). This improvement was similar to that among children treated with long-term inhaled corticosteroids. There was no additional improvement in PC₂₀by IPO13 variants among children treated with inhaled corticosteroids.

Conclusion

IPO13 variation is associated with improved AHR in asthmatic children. The degree of this improvement is similar to that observed with long-term inhaled corticosteroid treatment, suggesting that IPO13 variation may improve nuclear bioavailability of endogenous glucocorticoids.

Background

Endogenous glucocorticoids (GCs) serve a broad range of biological and physiological processes, including metabolic control, induction of anti-inflammatory cascades, and fetal organ maturation. In the context of the lung, GCs modulate airway branching morphogenesis[1,2], augment production of surfactant in late gestation[3], and promote alveolar formation[4]. As a result of their potent anti-inflammatory properties, exogenous GCs also serve as the most commonly used treatment for the long-term control of asthma[5] by effectively reducing airway hyperresponsiveness (AHR) and asthma symptoms[6], preventing exacerbations[7], and reducing asthma-associated mortality[8].

Though the precise molecular mechanisms that explain the diverse effects of GC have yet to be completely defined, nearly all GC effects result through cell-specific transcriptional regulation following binding of GC (in complex with the glucocorticoid receptor – GR) to positive and negative GC response elements in the promoter region of genes[9,10]. In its unbound state, GC is typically sequestered in the cytoplasm by molecular chaperones[9]. To access its genomic targets, GC-GR complexes must first pass from the cytoplasm to the nucleus through nuclear pore complexes. This process of GC- receptor shuttling across the nuclear-cytoplasmic membrane is tightly regulated by the interaction of cell-specific nuclear transport factors with cognate nuclear localization sequences[11]. Importin-13 (IPO13) was initially cloned in a search for GC-regulated genes important in lung development[12], and found to be differentially expressed during fetal lung growth with enrichment in lung epithelium relative to the mesenchyme. More recently, we (FK) demonstrated that IPO13 silencing prevents GC transport across the cytoplasmic-nuclear membrane in airway epithelium and abrogates GC-induced anti-inflammatory responses, suggesting that IPO13 is a critical nuclear transporter of GC receptor in the airway epithelium[13].

Although most asthmatics demonstrate an improvement in asthma control with long-term inhaled GC therapy, large inter-individual variation in response to inhaled GCs is well documented [14-16]. While some patients do respond to higher doses than normally prescribed, the administration of these doses for prolonged periods can have marked adverse effects. The identification and characterization of genetic determinants of GC responsiveness would provide insight into GC pharmacology and asthma pathogenesis. Given its demonstrated ability to mediate anti-inflammatory GC effects (particularly in airway epithelium), we considered that IPO13 represented a compelling biologic candidate gene for pharmacogenetic responsiveness to glucocorticoid therapy for asthma. To assess whether common IPO13 DNA sequence variants influence treatment response to inhaled corticosteroids, we genotyped 10 common IPO13 variants in a cohort of children with asthma participating in a clinical trial evaluating the long-term efficacy of inhaled anti-inflammatory medication (including budesonide, a commonly prescribed inhaled GC)[17]. Herein we report that IPO13 variants differentially influenced airway hyperresponsiveness by treatment group, with improvements in AHR noted among subjects who were randomized to either placebo or nedocromil to levels similar to subjects who were randomized to budesonide, suggesting that common IPO13 variants may increase the nuclear bioavailability of endogenous GCs.

Methods

Population

CAMP is a multicenter, randomized, double-blinded clinical trial testing the safety and efficacy of inhaled budesonide (200 ug twice daily) vs. nedocromil (8 mg twice daily) vs. placebo over a mean of 4.3 years. Trial design and methodology have been published[17,18]. CAMP enrolled 1,041 children ages 5 to 12 years with mild to moderate asthma. Entry criteria included asthma symptoms and/or medication use for ≥ 6 months in the previous year and airway responsiveness with PC₂₀≤ 12.5 mg/ml. Follow-up visits with spirometry occurred at two and four months and every four months thereafter. Methacholine studies were performed during the run-in period, at 8 months post-randomization, then yearly thereafter. 968 children and 1,518 parents contributed DNA samples, including those of self-reported white (654 children and 950 parents), African-American (131 and 128), and Hispanic (86 and 94) ancestry[19]. Given the relatively small sample sizes of the non-white ethnic groups in CAMP and to avoid spurious association due to population stratification, association analyses were restricted to white probands (Table (Table11).

Table 1

Baseline characteristics of probands in CAMP

Human Subjects

The Institutional Review Boards of the Brigham and Women's Hospital and of the other CAMP study centers approved this study. Informed assent and consent were obtained from the study participants and their parents to collect DNA for genetic studies.

SNP genotyping

SNPs were genotyped using SEQUENOM^®(Sequenom, San Diego, CA). Primers and reaction conditions are available upon request. One SNP (rs2301993) that failed was genotyped using a TaqMAN™ assay (PE Biosystems, Foster City, CA). All SNP passed quality control, including high genotype completion rates (>95%), less than 1% genotype discordance upon repeat genotyping of a random sample of ~5–10% of the cohort, and lack of parental-offspring genotype incompatibilities.

SNP discovery

Bidirectional dye-terminator sequencing was performed according to protocol (Applied Biosystems, Foster City, CA) targeting all exons, intron-exon boundaries and 1 kb of flanking genomic sequence at the IPO13 locus in 23 white CAMP subjects selected to ensure representation of all four IPO13 haplotypes. Primers were designed using Primer3, and sequence analysis was performed using the 3130 DNA Analyzer (ABI).

Statistical analysis

Linkage disequilibrium (LD) and haplotype block analysis was performed using Haploview[20]. We assessed methacholine PC₂₀as the primary outcome of interest as it was the lung-function related phenotype most impacted by inhaled glucocorticoid treatment in the CAMP trial. Methacholine-PC₂₀– the dose of methacholine at which a 20% drop in the Forced Expiratory Volume in one second (FEV₁) from baseline was observed – was log₁₀transformed to achieve a normal distribution. We adjusted all association tests for baseline (pre-randomization) PC₂₀, age, gender, height, study center and visit. The primary clinical trial demonstrated equivalence of nedocromil and placebo with respect changes in airways responsiveness, lung function, and other clinical outcomes[17]. We therefore grouped subjects in these two groups for comparison with those participants randomized to budesonide in order to maximize statistical power. The primary analysis was a longitudinal analysis of methacholine-PC₂₀using the PROC MIXED procedure in SAS assuming a power spatial variance-covariance structure, with random slopes and intercepts estimated using Maximum Likelihood[21,22]. We confirmed results using family-based methods (FBAT-PC[23] in PBAT[24]) to exclude spurious association due to occult population stratification and to perform haplotype association testing. Primary hypothesis testing was performed assuming additive genetic effects. However, because of the relatively small numbers of rare homozygotes within treatment strata, subsequent estimates of genetic effect, tests for gene-by-treatment group interaction, and haplotype analysis were performed assuming a dominant genetic model. Point estimates of the effect of IPO13 polymorphism carrier status on methacholine PC₂₀were determined using generalized linear models (PROC GLM).

To control type I error, the primary analysis of association of IPO13 SNPs with PC₂₀was adjusted for multiple comparisons based on the methods of Nyholt[25] as modified by Li and Ji[26], as implemented in SNPSpD http://gump.qimr.edu.au/general/daleN/SNPSpD/. We first calculated the effective number of independent marker loci tested (M_effLi), defined by (i) calculating the correlation matrix (i.e. pairwise LD) across all markers using all available genotype data; (ii) measuring the collective correlation across the set of markers as the variance of the eigenvalues from this LD matrix; and (iii) using this measure of collective correlation to calculate the proportional reduction in the number of independent markers. In the CAMP population, the calculated collective correlation was high (variance of the observed eigenvalues = 2.4144) resulting in an estimated M_effLiof 6 markers, and an adjusted alpha of 0.0085 (see Table Table2).2). We note that because our primary hypothesis testing evaluated two treatment states (budesonide group vs. other), we further adjusted for two sets of tests, resulting in an experiment-wise significance threshold of 0.00425.

Table 2

IPO13 SNP pairwise LD correlation matrix used for MeffLi estimation

Results

SNP genotyping

The ten IPO13 SNPs evaluated in this study were selected from build 119 of dbSNP in order to achieve an average spacing of approximately 1 SNP every 5 kb across the IPO13 locus and its 20 kb flanking sequence (Figure (Figure1a).1a). Among non-Hispanic whites, all 10 SNPs had a minor allele frequency (MAF) of at least 0.05, and pair-wise LD high (median pair-wise D' 0.99, interquartile range of 0.98–1.0; see Figure Figure1b),1b), suggesting very limited haplotype diversity. Indeed, the 10 SNP form only four common haplotypes (Figure (Figure1c),1c), which can be unambiguously classified by several possible combinations of three tagging-SNPs (one representative set denoted with arrowheads in Figure Figure1c).1c). The genotyped SNP capture nearly all reported common variation at the IPO13 locus, in that they efficiently tag (r²> 0.80) all but two of the 21 SNPs with available genotype information in the HapMap families of northern and western European ancestry. This high degree of LD implies substantial redundancy with regard to the number of independent observations made in the association studies to be described below. Using the spectral decomposition of matrices (see Methods and references [25] and [26]) we estimated that the 10 genotyped variants effectively represent only 6 independent markers for association testing.

Figure 1

Genetic structure of IPO13. Panel A: Relative position of 10 variants genotyped on physical map of chromosome 1p34 region including IPO13. DPH2 = diphthamide biosynthesis protein 2. ATP6V0B = B component of vacuolar ATPase. Panel B: Pair-wise linkage (more ...)

Associations of IPO13 with airways hyperresponsiveness in CAMP

As outlined in the introduction, we initially hypothesized that IPO13 polymorphisms would impact individual responsiveness to inhaled corticosteroid responsiveness for the treatment of asthma given the documented role of IPO13 as an active nuclear transporter for receptor-bound glucocorticoid. The CAMP clinical trial demonstrated that inhaled corticosteroid therapy (budesonide 200 ug administered twice daily), as compared to placebo or nedocromil, most significantly impacted methacholine PC₂₀measurements rather than other spirometric measures of lung function[17]. We therefore focused our analysis on whether IPO13 variants impacted PC₂₀over the course of the clinical trial, and whether any observed effects were limited to subjects randomized to inhaled budesonide. As demonstrated in Table Table3,3, all IPO13 variants tested (with the exception of rs2240447) demonstrated some evidence of association with methacholine PC₂₀during the clinical trial, with two variants (rs6671164 and rs2301993) demonstrating significant association after multiple comparison adjustment for 12 tests (6 effective independent markers by 2 treatment strata, corrected alpha = 0.00425). However, in contrast to our assumptions that significant differences would be due to an effect in the budesonide treated group, stratified analysis clearly demonstrated that all of the observed effects were due to differences in the subjects who were not randomized to budesonide. These associations were not likely due to occult population stratification, as follow-up family-based association testing (which is immune to the effects of population stratification) demonstrated similar patterns of association (Table (Table44).

Table 3

Impact of IPO13 polymorphisms on methacholine hyperresponsiveness (PC₂₀) in childhood asthma

Table 4

Family-based association analysis of IPO13 polymorphisms on methacholine hyperresponsiveness (PC₂₀) in childhood asthma.

Though this longitudinal analysis of quantitative measures collected over the course of the clinical trial provides greatest statistical power, clinical interpretation is often more easily appreciated from cross-sectional analyses at discrete time points. As such, we next quantified the impact of IPO13 variants on methacholine PC₂₀in each treatment strata by estimating the SNP-specific fold-change in geometric mean methacholine PC₂₀at 8-months post-randomization (the time point during the clinical trial when maximal treatment response was noted and with fewest missing measurements). As shown in Table Table55 and in keeping with the repeated measures analysis described above, carriers of IPO13 variants who were randomized to placebo or nedocromil demonstrated between a 1.6 and 2.3 fold increase in methacholine PC₂₀compared to non-carriers for most SNP, whereas no significant difference in methacholine PC₂₀was noted for any SNP among subjects randomized to budesonide. Subtle differences in statistical significance (but not genetic effect estimates) were observed between the repeated measures and 8-month analyses (for example rs2301992 was not statistically significant in the 8-month cross-section analysis, while rs2240447 and rs2486014 demonstrated statistical significance only in the repeated measures analysis), suggesting instability in the variances in effect for these SNP in comparison to the others. Nonetheless, the general patterns across the locus were similar in both analyses and suggest that common IPO13 variants are associated with reduced airways responsiveness.

Table 5

Fold-change in methacholine PC₂₀by treatment group at 8-months post-randomization among IPO13 variant carriers

Figure Figure22 illustrates the differential effects of one representative polymorphism (rs2428953) on airway responsiveness over the course of the clinical trial. As can be seen, T-allele carriers on placebo or nedocromil had improvements in PC₂₀values over the course of the trial that approached those observed among subjects who were treated with budesonide. In keeping with the longitudinal analysis, similar profiles were observed for the other IPO13 SNP, with carriers of the minor alleles at all loci (with the exception of rs2240447 and rs2301993) randomized to placebo or nedocromil demonstrating significantly less airways responsiveness compared to non-carriers (data not shown). From these data, we conclude that IPO13 variation influences the natural progression of methacholine PC₂₀among children with asthma, resulting in significantly less severe airway hyperresponsiveness, and approaching levels achieved with inhaled corticosteroid therapy. These data also suggest that the addition of inhaled corticosteroids confers no additional benefit among IPO13 carriers, though our study is underpowered to formally test that specific interaction. We also note that airway hyperresponsiveness was lowest among subjects homozygous for IPO13 variants (for example, the mean logPC₂₀was 1.31 among SNP rs2428953 TT homozygotes compared to 0.77 among heterozygotes), suggesting an additive genetic relationship, though the number of subjects with rare homozygous genotypes was generally too small (13 in the non-budesonide group, 2 in the budesonide group) to make formal statements regarding the significance of this observation.

Figure 2

Impact of IPO13 polymorphism rs2428953 on airway hyperresponsiveness. Mean log(methacholine PC₂₀) values and SEM. Common genotype (CC) denoted by solid black lines, heterozygotes and TT homozygotes denoted by dashed red lines. Open circles denote placebo/nedocromil (more ...)

Haplotype association analysis

As described above, four common (frequency > 1%) IPO13 haplotypes are observed in this cohort. We assessed whether the associations observed above could be attributed to specific haplotypes, we preformed a family-based haplotype association test using one representative set of haplotype-tagging SNP (rs1990150, rs2240447, and rs1636879 – arrow heads in Figure Figure1c).1c). Similar to the single SNP analysis, treatment-stratified haplotype analysis demonstrated association with methacholine PC₂₀among subjects randomized to placebo or nedocromil only (global p-value for test of association of IPO13 locus = 0.02), but not among subjects randomized to budesonide (p = 0.58). Inspection of the haplotype-specific p-values (Table (Table6)6) suggests that the primary association is with haplotype B (haplotype specific p-value = 0.02), though the relatively small number of informative families available for evaluating haplotypes C and D preclude definitive statements regarding these later haplotypes.

Table 6

Family-based haplotype association analysis of IPO13 polymorphisms on methacholine hyperresponsiveness (PC20) in childhood asthma.

SNP discovery

Resequencing of the IPO13 locus in 23 CAMP subjects selected to ensure representation of all four common IPO13 haplotypes identified 8 variants that had not been genotyped, including four not present in the dbSNP database (Table (Table7).7). Two were potentially functional (a non-synonymous Asp891Ser substitution and a highly conserved non-coding variant), though their low frequency (each was observed only once) and presence only on the common haplotype A background (not the PC₂₀associated background) suggest that they are not responsible for the observed associations with PC₂₀.

Table 7

IPO13 polymorphisms identified through SNP discovery effort

Discussion

Herein, we observed significant associations of common IPO13 polymorphisms with airway responsiveness (the most dynamic treatment response phenotype in the CAMP clinical trial) among children with mild-to-moderate asthma. These associations were observed exclusively among those children who were not randomized to inhaled GCs, and persisted over the ~4.5 years of clinical observation. The genetic effects conferred by these IPO13 variants were clinically significant, with an average 1.5–2.1 fold increase in mean PC₂₀among carriers of IPO13 variants.

IPO13 has been functionally characterized as a primary regulator of GC-bound GR transport across the nuclear membrane. Inhibition of lung epithelial cell IPO13 production inhibits nuclear translocation of GR from the cytoplasm and subsequent GC-mediated silencing of inflammatory cytokine production[13], suggesting that the normal anti-inflammatory response induced by GC is dependent on normal IPO13 function. In light of these observations, it is curious that the genetic effects observed in the current study were observed only among subjects who were not taking daily corticosteroids. We propose two possible mechanisms to explain these findings. The first is developmental. IPO13 was first identified in studies of lung development, where IPO13 (initially known as LGL2) was found to be differentially expressed in fetal rat lung cell culture[12]. It is conceivable that IPO13 variation impacts airway hyperresponsiveness by altering airway anatomy through changes during airway morphogenesis and development. Though plausible, the absence of association of IPO13 variants with baseline lung function in the current study (data not shown) suggests that this is not likely the case. However, given the diverse roles of glucocorticoid during lung morphogenesis (see Introduction) and the considerable impact on lung development by aberrant perinatal glucocorticoid exposure[27,28], we are hesitant to completely discount this possibility at this time.

A second possible mechanism of action is that IPO13 variants improve airway hyperresponsiveness by enhancing the local anti-inflammatory effects of circulating, endogenous GCs by facilitating increased GR transport into the nucleus and thus increasing the effective bioavailability of endogenous GC. Supporting this hypothesis is the observation that carriers of IPO13 variants who were not on inhaled steroid exhibited improvements in airway responsiveness over the course of the clinical trial that approached those for subjects who were taking inhaled steroids (see Figure Figure2),2), suggesting that that IPO13 variation enhances endogenous GC nuclear availability to therapeutic levels. It has been previously demonstrated that nucleocytoplasmic shuttling of IPO13 is developmentally regulated and highly variable in rat lung[12,13]. Sequence variation could potentially influence this regulation by increasing nuclear membrane availability through increased IPO13 expression or by altering the kinetic properties of GR transport by influencing GR binding affinities. We note that the paucity of coding variation at the IPO13 identified through our resequencing efforts suggests that the genetic effects observed are likely due to regulatory variation rather than structural changes. Though studies are currently ongoing to explore these possibilities, we note that because none of the airway responsiveness-associated non-coding variants map to highly conserved genetic sequence and are not predicted to harbor transcription factor binding sites, it is unlikely that we have as of yet identified a putative functional variant.

It is noteworthy that the associated haplotype block not only spans the IPO13 locus, but two neighboring genes as well: ATP6V0B and DHP2 (Figure (Figure1).1). Though it is not possible to completely exclude these other genes as functionally responsible for the observed associations, it is unlikely to be the case. The motivation for studying these polymorphisms was the recognition of IPO13 as the primary nuclear transporter for steroid-bound glucocorticoid. Had we identified these variants through a hypothesis-free approach (i.e. a genome-wide study), the pretest probability for each gene in the region would be similar. However, because these SNPs were chosen due to the biologic prior on IPO13, it would be improbable that the true functional effects would be mediated through a neighboring gene. Unlike IPO13, there is little biological evidence to support either ATP6V0B or DHP2 in either the pathogenesis of airways responsiveness or glucocorticoid pharmacogenetics. ATP6V0B is a subunit of the vacuolar-type H(+)-ATPase (V-ATPase) multisubunit enzyme, responsible for organelle acidification. DHP2 encodes a protein involved in diphthamide biosynthesis that confers resistance in yeast to the effects of diphtheria toxin. Though surveys of genomic databases (including UniGene and GEO) suggest that ATP6V06 is ubiquitously expressed and DHP2 is weakly expressed in the lung, there is little reason to suspect either as the responsible locus.

In genetic association studies of complex traits such as airway hyperresponsiveness and pharmacogenetic responses, it is important to consider a variety of methodological and statistical issues that can hamper proper interpretation of observed findings. Two features of the current study warrant particular attention: phenotype misclassification and statistical power. Most observational studies of the genetics of asthma attempt to avoid confounding of lung phenotype measurements by medication use by performing spirometric and airway hyperresponsiveness measurements following a short period (less than 24–48 hours) off anti-asthma controller medications. However, due to safety consideration, long-term avoidance of asthma controller medication is typically not permitted. The randomized, placebo-controlled trial is perhaps the only setting in which such confounding can be eliminated, and is a major strength of the study presented here. It is directly a result of the availability of ~4.5 years of repeated methacholine challenge measurements off anti-inflammatory agents among more than two-thirds of CAMP participants who were randomized to placebo or nedocromil that enabled detection of the observed associations in this study. In light of the differences in genetic effect observed across treatment arms, we stress that future attempts to replicate our findings should use approaches that address this issues, as replication may only be possible when proper adjustment for glucocorticoid use are made.

Subjects randomized to inhaled budesonide represent only 1/3 of the cohort. We recognize that this relatively smaller sample size is inadequately powered to detect strong genetic effects with alleles of modest frequency, and that it is possible that an effect similar to that observed among those not on steroids could potentially have been observed among steroid users had the number of subjects in this latter group approached those of the former. However, we note that sample size impacts only the ability to claim a statistically significant difference in genotype effect but does not in any way influence the absolute effect observed. In this study, the trends of association among those on budesonide were quite dissimilar to those observed among subjects not on budesonide, with all variants actually conferring greater airway responsiveness (though the confidence intervals in this group are quite broad and all span the null). Therefore, while we stress that our conclusions reported herein apply only to those subjects in the non-budesonide arm, it is likely that differential effects of these alleles are present between subjects who were and were not taking steroids.

Conclusion

We have demonstrated that genetic variation in IPO13 is associated with reduced airway hyperresponsiveness among children with mild-to-moderate asthma who are not using long-term inhaled corticosteroids. The degree of this reduction is similar to the improvements noted among those children on long-term corticosteroids, supporting the notion that IPO13 variation improves endogenous glucocorticoid bioavailability in the cell nucleus. We are currently performing functional studies using patient-derived cell lines representing the spectrum of IPO13 genetic variation to explore this possibility.

Abbreviations

AHR: Airway hyperresponsiveness; CAMP: Childhood Asthma Management Program; GC: Glucocorticoid; GR: glucocorticoid receptor; IPO13: Importin-13; Kb: kilobase; LD: linkage disequilibrium; MAF: minor allele frequency; PC₂₀: dose of methacholine bonchoprovocation test resulting in 20% in baseline forced expiratory volume in one second; SNP: single nucleotide polymorphism.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

BAR designed and supervised the genotyping and sequencing studies, performed the primary data analysis and drafted the manuscript. KVS and JLS performed the statistical analysis. BAR, KT, FK and STW conceived of the study, participated in its design and coordination. All authors read and approved the final manuscript.

Acknowledgements

We thank all subjects for their ongoing participation in this study. We acknowledge the CAMP investigators and research team, supported by the NHLBI, for collection of CAMP Genetic Ancillary Study data. All work on data collected from the CAMP Genetic Ancillary Study was conducted at the Channing Laboratory of the Brigham and Women's Hospital under appropriate CAMP policies and human subjects' protections. CAMP is supported by U01 HL075419, U01 HL65899, P01 HL083069, and T32 HL07427 from the National Heart, Lung and Blood Institute, National Institutes of Health (NHLBI/NIH). F.K. is supported by grant MOP74683 from the Canadian Institutes of Health Research (FK). B.A.R. is a recipient of a Mentored Clinical Scientist Development Award from NHLBI/NIH (K08 HL074193).

References

Oshika E, Liu S, Ung LP, Singh G, Shinozuka H, Michalopoulos GK, Katyal SL. Glucocorticoid-induced effects on pattern formation and epithelial cell differentiation in early embryonic rat lungs. Pediatr Res. 1998;43(3):305–314. doi: 10.1203/00006450-199803000-00001. [PubMed] [Cross Ref]
Oshika E, Liu S, Singh G, Michalopoulos GK, Shinozuka H, Katyal SL. Antagonistic effects of dexamethasone and retinoic acid on rat lung morphogenesis. Pediatr Res. 1998;43(3):315–324. doi: 10.1203/00006450-199803000-00002. [PubMed] [Cross Ref]
Guettari N, Marin L, Bourbon J, Dufour ME, Rieutort M, Tordet C. Effects of the antiglucocorticoid RU 486 on the maturation of fetal rat lung surfactant. Exp Lung Res. 1989;15(2):151–166. doi: 10.3109/01902148909087850. [PubMed] [Cross Ref]
Massaro GD, Massaro D. Retinoic acid treatment partially rescues failed septation in rats and in mice. Am J Physiol Lung Cell Mol Physiol. 2000;278(5):L955–960. [PubMed]
Gauvreau GM, Doctor J, Watson RM, Jordana M, O'Byrne PM. Effects of inhaled budesonide on allergen-induced airway responses and airway inflammation. Am J Respir Crit Care Med. 1996;154(5):1267–1271. [PubMed]
Haahtela T, Jarvinen M, Kava T, Kiviranta K, Koskinen S, Lehtonen K, Nikander K, Persson T, Reinikainen K, Selroos O, et al. Comparison of a beta 2-agonist, terbutaline, with an inhaled corticosteroid, budesonide, in newly detected asthma. N Engl J Med. 1991;325(6):388–392. [PubMed]
Pauwels RA, Lofdahl CG, Postma DS, Tattersfield AE, O'Byrne P, Barnes PJ, Ullman A. Effect of inhaled formoterol and budesonide on exacerbations of asthma. Formoterol and Corticosteroids Establishing Therapy (FACET) International Study Group. N Engl J Med. 1997;337(20):1405–1411. doi: 10.1056/NEJM199711133372001. [PubMed] [Cross Ref]
Suissa S, Ernst P, Benayoun S, Baltzan M, Cai B. Low-dose inhaled corticosteroids and the prevention of death from asthma. N Engl J Med. 2000;343(5):332–336. doi: 10.1056/NEJM200008033430504. [PubMed] [Cross Ref]
Barnes PJ. How corticosteroids control inflammation: Quintiles Prize Lecture 2005. Br J Pharmacol. 2006;148(3):245–254. doi: 10.1038/sj.bjp.0706736. [PubMed] [Cross Ref]
Schaaf MJ, Cidlowski JA. Molecular determinants of glucocorticoid receptor mobility in living cells: the importance of ligand affinity. Mol Cell Biol. 2003;23(6):1922–1934. doi: 10.1128/MCB.23.6.1922-1934.2003. [PMC free article] [PubMed] [Cross Ref]
Picard D, Yamamoto KR. Two signals mediate hormone-dependent nuclear localization of the glucocorticoid receptor. Embo J. 1987;6(11):3333–3340. [PubMed]
Zhang C, Sweezey NB, Gagnon S, Muskat B, Koehler D, Post M, Kaplan F. A novel karyopherin-beta homolog is developmentally and hormonally regulated in fetal lung. Am J Respir Cell Mol Biol. 2000;22(4):451–459. [PubMed]
Tao T, Lan J, Lukacs GL, Hache RJ, Kaplan F. Importin 13 regulates nuclear import of the glucocorticoid receptor in airway epithelial cells. Am J Respir Cell Mol Biol. 2006;35(6):668–680. doi: 10.1165/rcmb.2006-0073OC. [PubMed] [Cross Ref]
Kerrebijn KF, van Essen-Zandvliet EE, Neijens HJ. Effect of long-term treatment with inhaled corticosteroids and beta-agonists on the bronchial responsiveness in children with asthma. J Allergy Clin Immunol. 1987;79(4):653–659. doi: 10.1016/S0091-6749(87)80163-X. [PubMed] [Cross Ref]
Szefler SJ, Martin RJ, King TS, Boushey HA, Cherniack RM, Chinchilli VM, Craig TJ, Dolovich M, Drazen JM, Fagan JK, Fahy JV, Fish JE, Ford JG, Israel E, Kiley J, Kraft M, Lazarus SC, Lemanske RF, Jr, Mauger E, Peters SP, Sorkness CA. Significant variability in response to inhaled corticosteroids for persistent asthma. J Allergy Clin Immunol. 2002;109(3):410–418. doi: 10.1067/mai.2002.122635. [PubMed] [Cross Ref]
Malmstrom K, Rodriguez-Gomez G, Guerra J, Villaran C, Pineiro A, Wei LX, Seidenberg BC, Reiss TF. Oral montelukast, inhaled beclomethasone, and placebo for chronic asthma. A randomized, controlled trial. Montelukast/Beclomethasone Study Group. Ann Intern Med. 1999;130(6):487–495. [PubMed]
The Childhood Asthma Management Program Research Group Long-term effects of budesonide or nedocromil in children with asthma. N Engl J Med. 2000;343(15):1054–1063. doi: 10.1056/NEJM200010123431501. [PubMed] [Cross Ref]
Childhood Asthma Management Program Research Group The Childhood Asthma Management Program (CAMP): design, rationale, and methods. Control Clin Trials. 1999;20(1):91–120. doi: 10.1016/S0197-2456(98)00044-0. [PubMed] [Cross Ref]
Raby BA, Klanderman B, Murphy A, Mazza S, Camargo CA, Jr, Silverman EK, Weiss ST. A common mitochondrial haplogroup is associated with elevated total serum IgE levels. J Allergy Clin Immunol. 2007;120(2):351–358. doi: 10.1016/j.jaci.2007.05.029. [PubMed] [Cross Ref]
Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21(2):263–265. doi: 10.1093/bioinformatics/bth457. [PubMed] [Cross Ref]
Singer JD. Using SAS PROC MIXED to fit multilevel models, hierarchical models, and individual growth models. Journal of Educational and Behavioral Statistics. 1998;24(4):323–355.
Diggle PJ, Heagerty P, Linang KY, Zeger SL. Analysis of longitudinal data. 2. Oxford: Oxford University Press; 2002.
Lange C, van Steen K, Andrew T, Lyon H, DeMeo DL, Raby BA, Murphy A, Silverman EK, MacGregor A, Weiss ST, Laird N. A Family-Based Association Test for Repeatedly Measured Quantitative Traits Adjusting for Unknown Environmental and/or Polygenic Effects. Statistical Applications in Genetics and Molecular Biology. 2004;3(1):17. doi: 10.2202/1544-6115.1067. [PubMed] [Cross Ref]
Lange C, DeMeo D, Silverman EK, Weiss ST, Laird NM. PBAT: tools for family-based association studies. Am J Hum Genet. 2004;74(2):367–369. doi: 10.1086/381563. [PubMed] [Cross Ref]
Nyholt DR. A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other. Am J Hum Genet. 2004;74(4):765–769. doi: 10.1086/383251. [PubMed] [Cross Ref]
Li J, Ji L. Adjusting multiple testing in multilocus analyses using the eigenvalues of a correlation matrix. Heredity. 2005;95(3):221–227. doi: 10.1038/sj.hdy.6800717. [PubMed] [Cross Ref]
Frank L, Roberts RJ. Effects of low-dose prenatal corticosteroid administration on the premature rat. Biol Neonate. 1979;36(1–2):1–9. [PubMed]
Jobe AH, Newnham JP, Willet KE, Sly P, Ervin MG, Bachurski C, Possmayer F, Hallman M, Ikegami M. Effects of antenatal endotoxin and glucocorticoids on the lungs of preterm lambs. Am J Obstet Gynecol. 2000;182(2):401–408. doi: 10.1016/S0002-9378(00)70231-6. [PubMed] [Cross Ref]
den Dunnen JT, Antonarakis SE. Nomenclature for the description of human sequence variations. Hum Genet. 2001;109(1):121–124. doi: 10.1007/s004390100505. [PubMed] [Cross Ref]
Siepel A, Bejerano G, Pedersen JS, Hinrichs AS, Hou M, Rosenbloom K, Clawson H, Spieth J, Hillier LW, Richards S, Weinstock GM, Wilson RK, Gibbs RA, Kent WJ, Miller W, Haussler D. Evolutionarily conserved elements in vertebrate, insect, worm, and yeast genomes. Genome Res. 2005;15(8):1034–1050. doi: 10.1101/gr.3715005. [PubMed] [Cross Ref]

Articles from Respiratory Research are provided here courtesy of
BioMed Central