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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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)

The impact of self-identified race on epidemiologic studies of gene expression.

Although population differences in gene expression have been established, the impact on differential gene expression studies in large populations is not well understood. We describe the effect of self-reported race on a gene expression study of lung function in asthma. We generated gene expression profiles for 254 young adults (205 non-Hispanic whites and 49 African Americans) with asthma on whom concurrent total RNA derived from peripheral blood CD4(+) lymphocytes and lung function measurements were obtained. We identified four principal components that explained 62% of the variance in gene expression. The dominant principal component, which explained 29% of the total variance in gene expression, was strongly associated with self-identified race (P<10(-16)). The impact of these racial differences was observed when we performed differential gene expression analysis of lung function. Using multivariate linear models, we tested whether gene expression was associated with a quantitative measure of lung function: pre-bronchodilator forced expiratory volume in one second (FEV(1)). Though unadjusted linear models of FEV(1) identified several genes strongly correlated with lung function, these correlations were due to racial differences in the distribution of both FEV(1) and gene expression, and were no longer statistically significant following adjustment for self-identified race. These results suggest that self-identified race is a critical confounding covariate in epidemiologic studies of gene expression and that, similar to genetic studies, careful consideration of self-identified race in gene expression profiling studies is needed to avoid spurious association.

Mapping of numerous disease-associated expression polymorphisms in primary peripheral blood CD4+ lymphocytes.

Genome-wide association studies of human gene expression promise to identify functional regulatory genetic variation that contributes to phenotypic diversity. However, it is unclear how useful this approach will be for the identification of disease-susceptibility variants. We generated gene expression profiles for 22 184 mRNA transcripts using RNA derived from peripheral blood CD4+ lymphocytes, and genome-wide genotype data for 516 512 autosomal markers in 200 subjects. We screened for cis-acting variants by testing variants mapping within 50 kb of expressed transcripts for association with transcript abundance using generalized linear models. Significant associations were identified for 1585 genes at a false discovery rate of 0.05 (corresponding to P-values ranging from 1 × 10(-91) to 7 × 10(-4)). Importantly, we identified evidence of regulatory variation for 119 previously mapped disease genes, including 24 examples where the variant with the strongest evidence of disease-association demonstrates strong association with specific transcript abundance. The prevalence of cis-acting variants among disease-associated genes was 63% higher than the genome-wide rate in our data set (P = 6.41 × 10(-6)), and although many of the implicated loci were associated with immune-related diseases (including asthma, connective tissue disorders and inflammatory bowel disease), associations with genes implicated in non-immune-related diseases including lipid profiles, anthropomorphic measurements, cancer and neurologic disease were also observed. Genetic variants that confer inter-individual differences in gene expression represent an important subset of variants that contribute to disease susceptibility. Population-based integrative genetic approaches can help identify such variation and enhance our understanding of the genetic basis of complex traits.

PRKCA: a positional candidate gene for body mass index and asthma.

Asthma incidence and prevalence are higher in obese individuals. A potential mechanistic basis for this relationship is pleiotropy. We hypothesized that significant linkage and candidate-gene association would be found for body mass index (BMI) in a population ascertained on asthma affection status. Linkage analysis for BMI was performed on 657 subjects in eight Costa Rican families enrolled in a study of asthma. Family-based association studies were conducted for BMI with SNPs within a positional candidate gene, PRKCA. SNPs within PRKCA were also tested for association with asthma. Association studies were conducted in 415 Costa Rican parent-child trios and 493 trios participating in the Childhood Asthma Management Program (CAMP). Although only modest evidence of linkage for BMI was obtained for the whole cohort, significant linkage was noted for BMI in females on chromosome 17q (peak LOD = 3.39). Four SNPs in a candidate gene in this region (PRKCA) had unadjusted association p values < 0.05 for BMI in both cohorts, with the joint p value for two SNPs remaining significant after adjustment for multiple comparisons (rs228883 and rs1005651, joint p values = 9.5 x 10(-)(5) and 5.6 x 10(-)(5)). Similarly, eight SNPs had unadjusted association p values < 0.05 for asthma in both populations, with one SNP remaining significant after adjustment for multiple comparisons (rs11079657, joint p value = 2.6 x 10(-)(5)). PRKCA is a pleiotropic locus that is associated with both BMI and asthma and that has been identified via linkage analysis of BMI in a population ascertained on asthma.

Genome-wide association analysis identifies PDE4D as an asthma-susceptibility gene.

Asthma, a chronic airway disease with known heritability, affects more than 300 million people around the world. A genome-wide association (GWA) study of asthma with 359 cases from the Childhood Asthma Management Program (CAMP) and 846 genetically matched controls from the Illumina ICONdb public resource was performed. The strongest region of association seen was on chromosome 5q12 in PDE4D. The phosphodiesterase 4D, cAMP-specific (phosphodiesterase E3 dunce homolog, Drosophila) gene (PDE4D) is a regulator of airway smooth-muscle contractility, and PDE4 inhibitors have been developed as medications for asthma. Allelic p values for top SNPs in this region were 4.3 x 10(-07) for rs1588265 and 9.7 x 10(-07) for rs1544791. Replications were investigated in ten independent populations with different ethnicities, study designs, and definitions of asthma. In seven white and Hispanic replication populations, two PDE4D SNPs had significant results with p values less than 0.05, and five had results in the same direction as the original population but had p values greater than 0.05. Combined p values for 18,891 white and Hispanic individuals (4,342 cases) in our replication populations were 4.1 x 10(-04) for rs1588265 and 9.2 x 10(-04) for rs1544791. In three black replication populations, which had different linkage disequilibrium patterns than the other populations, original findings were not replicated. Further study of PDE4D variants might lead to improved understanding of the role of PDE4D in asthma pathophysiology and the efficacy of PDE4 inhibitor medications.

MMP12, lung function, and COPD in high-risk populations.

BACKGROUND: Genetic variants influencing lung function in children and adults may ultimately lead to the development of chronic obstructive pulmonary disease (COPD), particularly in high-risk groups. METHODS: We tested for an association between single-nucleotide polymorphisms (SNPs) in the gene encoding matrix metalloproteinase 12 (MMP12) and a measure of lung function (prebronchodilator forced expiratory volume in 1 second [FEV(1)]) in more than 8300 subjects in seven cohorts that included children and adults. Within the Normative Aging Study (NAS), a cohort of initially healthy adult men, we tested for an association between SNPs that were associated with FEV(1) and the time to the onset of COPD. We then examined the relationship between MMP12 SNPs and COPD in two cohorts of adults with COPD or at risk for COPD. RESULTS: The minor allele (G) of a functional variant in the promoter region of MMP12 (rs2276109 [-82A-->G]) was positively associated with FEV(1) in a combined analysis of children with asthma and adult former and current smokers in all cohorts (P=2x10(-6)). This allele was also associated with a reduced risk of the onset of COPD in the NAS cohort (hazard ratio, 0.65; 95% confidence interval [CI], 0.46 to 0.92; P=0.02) and with a reduced risk of COPD in a cohort of smokers (odds ratio, 0.63; 95% CI, 0.45 to 0.88; P=0.005) and among participants in a family-based study of early-onset COPD (P=0.006). CONCLUSIONS: The minor allele of a SNP in MMP12 (rs2276109) is associated with a positive effect on lung function in children with asthma and in adults who smoke. This allele is also associated with a reduced risk of COPD in adult smokers.

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.

Regulatory haplotypes in ARG1 are associated with altered bronchodilator response.

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).

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.

Genome-wide association analysis of body mass in chronic obstructive pulmonary disease.

Cachexia, whether assessed by body mass index (BMI) or fat-free mass index (FFMI), affects a significant proportion of patients with chronic obstructive pulmonary disease (COPD), and is an independent risk factor for increased mortality, increased emphysema, and more severe airflow obstruction. The variable development of cachexia among patients with COPD suggests a role for genetic susceptibility. The objective of the present study was to determine genetic susceptibility loci involved in the development of low BMI and FFMI in subjects with COPD. A genome-wide association study (GWAS) of BMI was conducted in three independent cohorts of European descent with Global Initiative for Chronic Obstructive Lung Disease stage II or higher COPD: Evaluation of COPD Longitudinally to Identify Predictive Surrogate End-Points (ECLIPSE; n = 1,734); Norway-Bergen cohort (n = 851); and a subset of subjects from the National Emphysema Treatment Trial (NETT; n = 365). A genome-wide association of FFMI was conducted in two of the cohorts (ECLIPSE and Norway). In the combined analyses, a significant association was found between rs8050136, located in the first intron of the fat mass and obesity-associated (FTO) gene, and BMI (P = 4.97 × 10(-7)) and FFMI (P = 1.19 × 10(-7)). We replicated the association in a fourth, independent cohort consisting of 502 subjects with COPD from COPDGene (P = 6 × 10(-3)). Within the largest contributing cohort of our analysis, lung function, as assessed by forced expiratory volume at 1 second, varied significantly by FTO genotype. Our analysis suggests a potential role for the FTO locus in the determination of anthropomorphic measures associated with COPD.