Automated Pharmacogenomic Reports for Clinical Genome Sequencing.

Clinical laboratories offering genome sequencing have the opportunity to return pharmacogenomic findings to patients, providing the added benefit of preemptive testing that could help inform medication selection or dosing throughout the lifespan. Implementation of pharmacogenomic reporting must address several challenges, including inherent limitations in short-read genome sequencing methods, gene and variant selection, standardization of genotype and phenotype nomenclature, and choice of guidelines and drugs to report. An automated pipeline, lmPGX, was developed as an end-to-end solution that produces two versions of a pharmacogenomic report, presenting either Clinical Pharmacogenetics Implementation Consortium or US Food and Drug Administration guidelines for 12 genes. The pipeline was validated for performance using reference samples and pharmacogenetic data from the Genetic Testing Reference Materials Coordination Program. To determine performance and limitations, lmPGX was compared with three additional publicly available pharmacogenomic pipelines. The lmPGX pipeline offers clinical laboratories an opportunity for seamless integration of pharmacogenomic results with genome reporting.

Vitamin D Modulates Expression of the Airway Smooth Muscle Transcriptome in Fatal Asthma.

Globally, asthma is a chronic inflammatory respiratory disease affecting over 300 million people. Some asthma patients remain poorly controlled by conventional therapies and experience more life-threatening exacerbations. Vitamin D, as an adjunct therapy, may improve disease control in severe asthma patients since vitamin D enhances glucocorticoid responsiveness and mitigates airway smooth muscle (ASM) hyperplasia. We sought to characterize differences in transcriptome responsiveness to vitamin D between fatal asthma- and non-asthma-derived ASM by using RNA-Seq to measure ASM transcript expression in five donors with fatal asthma and ten non-asthma-derived donors at baseline and with vitamin D treatment. Based on a Benjamini-Hochberg corrected p-value <0.05, 838 genes were differentially expressed in fatal asthma vs. non-asthma-derived ASM at baseline, and vitamin D treatment compared to baseline conditions induced differential expression of 711 and 867 genes in fatal asthma- and non-asthma-derived ASM, respectively. Functional gene categories that were highly represented in all groups included extracellular matrix, and responses to steroid hormone stimuli and wounding. Genes differentially expressed by vitamin D also included cytokine and chemokine activity categories. Follow-up qPCR and individual analyte ELISA experiments were conducted for four cytokines (i.e. CCL2, CCL13, CXCL12, IL8) to measure TNFα-induced changes by asthma status and vitamin D treatment. Vitamin D inhibited TNFα-induced IL8 protein secretion levels to a comparable degree in fatal asthma- and non-asthma-derived ASM even though IL8 had significantly higher baseline levels in fatal asthma-derived ASM. Our findings identify vitamin D-specific gene targets and provide transcriptomic data to explore differences in the ASM of fatal asthma- and non-asthma-derived donors.

Genome-wide association study of short-acting ß2-agonists. A novel genome-wide significant locus on chromosome 2 near ASB3.

RATIONALE: ß2-Agonists are the most common form of treatment of asthma, but there is significant variability in response to these medications. A significant proportion of this responsiveness may be heritable. OBJECTIVES: To investigate whether a genome-wide association study (GWAS) could identify novel pharmacogenetic loci in asthma. METHODS: We performed a GWAS of acute bronchodilator response (BDR) to inhaled ß2-agonists. A total of 444,088 single-nucleotide polymorphisms (SNPs) were examined in 724 individuals from the SNP Health Association Resource (SHARe) Asthma Resource Project (SHARP). The top 50 SNPs were carried forward to replication in a population of 444 individuals. MEASUREMENTS AND MAIN RESULTS: The combined P value for four SNPs reached statistical genome-wide significance aftercorrecting for multiple comparisons. Combined P values for rs350729, rs1840321, rs1384918, and rs1319797 were 2.21 × 10(-10), 5.75 × 10(-8), 9.3 × 10(-8), and 3.95 × 10(-8), respectively. The significant variants all map to a novel genetic region on chromosome 2 near the ASB3 gene, a region associated with smooth muscle proliferation. As compared with the wild type, the presence of the minor alleles reduced the degree of BDR by 20% in the original population and by a similar percentage in the confirmatory population. CONCLUSIONS: These GWAS findings for BDR in subjects with asthma suggest that a gene associated with smooth muscle proliferation may influence a proportion of the smooth muscle relaxation that occurs in asthma.

RNA-Seq transcriptome profiling identifies CRISPLD2 as a glucocorticoid responsive gene that modulates cytokine function in airway smooth muscle cells.

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.

Pharmacogenomics: novel loci identification via integrating gene differential analysis and eQTL analysis.

Nearly one-half of asthmatic patients do not respond to the most commonly prescribed controller therapy, inhaled corticosteroids (ICS). We conducted an expression quantitative trait loci (eQTL) analysis using >300 expression microarrays (from 117 lymphoblastoid cell lines) in corticosteroid (dexamethasone) treated and untreated cells derived from asthmatic subjects in the Childhood Asthma Management Program (CAMP) clinical trial. We then tested the associations of eQTL with longitudinal change in airway responsiveness to methacholine (LnPC20) on ICS. We identified 2484 cis-eQTL affecting 767 genes following dexamethasone treatment. A significant over-representation of lnPC20-associated cis-eQTL [190 single-nucleotide polymorphisms (SNPs)] among differentially expressed genes (odds ratio = 1.76, 95% confidence interval: 1.35-2.29) was noted in CAMP Caucasians. Forty-six of these 190 clinical associations were replicated in CAMP African Americans, including seven SNPs near six genes meeting criteria for genome-wide significance (P < 2 × 10(-7)). Notably, the majority of genome-wide findings would not have been uncovered via analysis of untreated samples. These results indicate that identifying eQTL after relevant environmental perturbation enables identification of true pharmacogenetic variants.

Vitamin D related genes in lung development and asthma pathogenesis.

BACKGROUND: 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. METHODS: Vitamin D pathway genes were derived from PubMed and Gene Ontology surveys. Principal component analysis was used to identify characteristic lung development genes. RESULTS: 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. CONCLUSIONS: 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.

ITGB5 and AGFG1 variants are associated with severity of airway responsiveness.

BACKGROUND: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: 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.

Integration of mouse and human genome-wide association data identifies KCNIP4 as an asthma gene.

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.

Genome-wide association analysis in asthma subjects identifies SPATS2L as a novel bronchodilator response gene.

Bronchodilator response (BDR) is an important asthma phenotype that measures reversibility of airway obstruction by comparing lung function (i.e. FEV(1)) 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 FEV(1) 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.

Association of cigarette smoking and CRP levels with DNA methylation in α-1 antitrypsin deficiency.

Alpha-1 antitrypsin (AAT) deficiency and tobacco smoking are confirmed risk factors for Chronic Obstructive Pulmonary Disease. We hypothesized that variable DNA methylation would be associated with smoking and inflammation, as reflected by the level of C-Reactive Protein (CRP) in AAT-deficient subjects. Methylation levels of 1,411 autosomal CpG sites from the Illumina GoldenGate Methylation Cancer Panel I were analyzed in 316 subjects. Associations of five smoking behaviors and CRP levels with individual CpG sites and average methylation levels were assessed using non-parametric testing, linear regression and linear mixed effect models, with and without adjustment for age and gender. Univariate linear regression analysis revealed that methylation levels of 16 CpG sites significantly associated with ever-smoking status. A CpG site in the TGFBI gene was the only site associated with ever-smoking after adjustment for age and gender. No highly significant associations existed between age at smoking initiation, pack-years smoked, duration of smoking, and time since quitting smoking as predictors of individual CpG site methylation levels. However, ever-smoking and younger age at smoking initiation associated with lower methylation level averaged across all sites. DNA methylation at CpG sites in the RUNX3, JAK3 and KRT1 genes associated with CRP levels. The most significantly associated CpG sites with gender and age mapped to the CASP6 and FZD9 genes, respectively. In summary, this study identified multiple potential candidate CpG sites associated with ever-smoking and CRP level in AAT-deficient subjects. Phenotypic variability in Mendelian diseases may be due to epigenetic factors.

Genome-wide association study of the age of onset of childhood asthma.

BACKGROUND: 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. OBJECTIVE: We sought to identify genetic variants associated with earlier onset of childhood asthma. METHODS: We conducted the first genome-wide association study of the age of onset of childhood asthma among participants in the Childhood Asthma Management Program (CAMP) and used 3 independent cohorts from North America, Costa Rica, and Sweden for replication. RESULTS: Two single nucleotide polymorphisms (SNPs) were associated with earlier onset of asthma in the combined analysis of CAMP and the replication cohorts: rs9815663 (Fisher P= 2.31 × 10(-8)) and rs7927044 (P= 6.54 × 10(-9)). Of these 2 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)

Genome-wide association identifies the T gene as a novel asthma pharmacogenetic locus.

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 FEV(1) 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 FEV(1) 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.

A genome-wide association study of COPD identifies a susceptibility locus on chromosome 19q13.

The genetic risk factors for chronic obstructive pulmonary disease (COPD) are still largely unknown. To date, genome-wide association studies (GWASs) of limited size have identified several novel risk loci for COPD at CHRNA3/CHRNA5/IREB2, HHIP and FAM13A; additional loci may be identified through larger studies. We performed a GWAS using a total of 3499 cases and 1922 control subjects from four cohorts: the Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE); the Normative Aging Study (NAS) and National Emphysema Treatment Trial (NETT); Bergen, Norway (GenKOLS); and the COPDGene study. Genotyping was performed on Illumina platforms with additional markers imputed using 1000 Genomes data; results were summarized using fixed-effect meta-analysis. We identified a new genome-wide significant locus on chromosome 19q13 (rs7937, OR = 0.74, P = 2.9 × 10(-9)). Genotyping this single nucleotide polymorphism (SNP) and another nearby SNP in linkage disequilibrium (rs2604894) in 2859 subjects from the family-based International COPD Genetics Network study (ICGN) demonstrated supportive evidence for association for COPD (P = 0.28 and 0.11 for rs7937 and rs2604894), pre-bronchodilator FEV(1) (P = 0.08 and 0.04) and severe (GOLD 3&4) COPD (P = 0.09 and 0.017). This region includes RAB4B, EGLN2, MIA and CYP2A6, and has previously been identified in association with cigarette smoking behavior.

Genome-wide association study reveals class I MHC-restricted T cell-associated molecule gene (CRTAM) variants interact with vitamin D levels to affect asthma exacerbations.

BACKGROUND: 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. OBJECTIVE: 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. CONCLUSION: 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.

Variable DNA methylation is associated with chronic obstructive pulmonary disease and lung function.

RATIONALE: Chronic obstructive pulmonary disease (COPD) is associated with local (lung) and systemic (blood) inflammation and manifestations. DNA methylation is an important regulator of gene transcription, and global and specific gene methylation marks may vary with cigarette smoke exposure. OBJECTIVES: To perform a comprehensive assessment of methylation marks in DNA from subjects well phenotyped for nonneoplastic lung disease. METHODS: We conducted array-based methylation screens, using a test-replication approach, in two family-based cohorts (n = 1,085 and 369 subjects). MEASUREMENTS AND MAIN RESULTS: We observed 349 CpG sites significantly associated with the presence and severity of COPD in both cohorts. Seventy percent of the associated CpG sites were outside of CpG islands, with the majority of CpG sites relatively hypomethylated. Gene ontology analysis based on these 349 CpGs (330 genes) suggested the involvement of a number of genes responsible for immune and inflammatory system pathways, responses to stress and external stimuli, as well as wound healing and coagulation cascades. Interestingly, our observations include significant, replicable associations between SERPINA1 hypomethylation and COPD and lower average lung function phenotypes (combined P values: COPD, 1.5 × 10(-23); FEV(1)/FVC, 1.5 × 10(-35); FEV(1), 2.2 × 10(-40)). CONCLUSIONS: Genetic and epigenetic pathways may both contribute to COPD. Many of the top associations between COPD and DNA methylation occur in biologically plausible pathways. This large-scale analysis suggests that DNA methylation may be a biomarker of COPD and may highlight new pathways of COPD pathogenesis.

Identification of a chronic obstructive pulmonary disease genetic determinant that regulates HHIP.

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.

Genomewide association between GLCCI1 and response to glucocorticoid therapy in asthma.

BACKGROUND: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: A functional GLCCI1 variant is associated with substantial decrements in the response to inhaled glucocorticoids in patients with asthma.

Association of SERPINE2 with asthma.

BACKGROUND: The "Dutch hypothesis" suggests that asthma and COPD have common genetic determinants. The serpin peptidase inhibitor, clade E (nexin, plasminogen activator inhibitor type 1), member 2 (SERPINE2) gene previously has been associated with COPD. We sought to determine whether SERPINE2 is associated with asthma and asthma-related phenotypes. METHODS: We measured the association of 39 SERPINE2 single-nucleotide polymorphisms (SNPs) with asthma-related phenotypes in 655 parent-child trios from the Childhood Asthma Management Program (CAMP), and we measured the association of 19 SERPINE2 SNPs with asthma in a case-control design of 359 CAMP probands and 846 population control subjects. We attempted to replicate primary asthma-related phenotype findings in one independent population and primary asthma affection status findings in two independent populations. We compared association results with CAMP proband expression quantitative trait loci. RESULTS: Nine of 39 SNPs had P < .05 for at least one phenotype in CAMP, and two of these replicated in an independent population of 426 people with childhood asthma. Six of 19 SNPs had P < .05 for association with asthma in CAMP/Illumina. None of these replicated in two independent populations. The expression quantitative trait loci revealed that five SNPs associated with asthma in CAMP/Illumina and one SNP associated with FEV(1) in CAMP are strongly correlated with SERPINE2 expression levels. Comparison of results to previous COPD studies identified five SNPs associated with both asthma- and COPD-related phenotypes. CONCLUSIONS: Our results weakly support SERPINE2 as a Dutch hypothesis candidate gene through nominally significant associations with asthma and related traits. Further study of SERPINE2 is necessary to verify its involvement in asthma and COPD.

SOX5 is a candidate gene for chronic obstructive pulmonary disease susceptibility and is necessary for lung development.

RATIONALE: Chromosome 12p has been linked to chronic obstructive pulmonary disease (COPD) in the Boston Early-Onset COPD Study (BEOCOPD), but a susceptibility gene in that region has not been identified. OBJECTIVES: We used high-density single-nucleotide polymorphism (SNP) mapping to implicate a COPD susceptibility gene and an animal model to determine the potential role of SOX5 in lung development and COPD. METHODS: On chromosome 12p, we genotyped 1,387 SNPs in 386 COPD cases from the National Emphysema Treatment Trial and 424 control smokers from the Normative Aging Study. SNPs with significant associations were then tested in the BEOCOPD study and the International COPD Genetics Network. Based on the human results, we assessed histology and gene expression in the lungs of Sox5(-/-) mice. MEASUREMENTS AND MAIN RESULTS: In the case-control analysis, 27 SNPs were significant at P ≤ 0.01. The most significant SNP in the BEOCOPD replication was rs11046966 (National Emphysema Treatment Trial-Normative Aging Study P = 6.0 × 10(-4), BEOCOPD P = 1.5 × 10(-5), combined P = 1.7 × 10(-7)), located 3' to the gene SOX5. Association with rs11046966 was not replicated in the International COPD Genetics Network. Sox5(-/-) mice showed abnormal lung development, with a delay in maturation before the saccular stage, as early as E16.5. Lung pathology in Sox5(-/-) lungs was associated with a decrease in fibronectin expression, an extracellular matrix component critical for branching morphogenesis. CONCLUSIONS: Genetic variation in the transcription factor SOX5 is associated with COPD susceptibility. A mouse model suggests that the effect may be due, in part, to its effects on lung development and/or repair processes.