Circulating microRNAs and prediction of asthma exacerbation in childhood asthma.

BACKGROUND: Circulating microRNAs have shown promise as non-invasive biomarkers and predictors of disease activity. Prior asthma studies using clinical, biochemical and genomic data have not shown excellent prediction of exacerbation. We hypothesized that a panel of circulating microRNAs in a pediatric asthma cohort combined with an exacerbation clinical score might predict exacerbation better than the latter alone. METHODS: Serum samples from 153 children at randomization in the Childhood Asthma Management Program were profiled for 754 microRNAs. Data dichotomized for asthma exacerbation one year after randomization to inhaled corticosteroid treatment were used for binary logistic regression with miRNA expressions and exacerbation clinical score. RESULTS: 12 of 125 well-detected circulating microRNAs had significant odd ratios for exacerbation with miR-206 being most significant. Each doubling of expression of the 12 microRNA corresponded to a 25-67% increase in exacerbation risk. Stepwise logistic regression yielded a 3-microRNA model (miR-146b, miR-206 and miR-720) that, combined with the exacerbation clinical score, had excellent predictive power with a 0.81 AUROC. These 3 microRNAs were involved in NF-kß and GSK3/AKT pathways. CONCLUSIONS: This combined circulating microRNA-clinical score model predicted exacerbation in asthmatic subjects on inhaled corticosteroids better than each constituent feature alone. TRIAL REGISTRATION: ClinicalTrials.gov Identifier: NCT00000575 .

Circulating microRNAs and association with methacholine PC20 in the Childhood Asthma Management Program (CAMP) cohort.

INTRODUCTION: Circulating microRNAs (miRNA) are promising biomarkers for human diseases. Our study hypothesizes that circulating miRNA would reveal candidate biomarkers related to airway hyperresponsiveness (AHR) and provide biologic insights into asthma epigenetic influences. METHODS: Serum samples obtained at randomization for 160 children in the Childhood Asthma Management Program were profiled using a TaqMan miRNA array set. The association of the isolated miRNA with methacholine PC20 was assessed. Network and pathway analyses were performed. Functional validation of two significant miRNAs was performed in human airway smooth muscle cells (HASMs). RESULTS: Of 155 well-detected circulating miRNAs, eight were significantly associated with PC20 with the strongest association with miR-296-5p. Pathway analysis revealed miR-16-5p as a network hub, and involvement of multiple miRNAs interacting with genes in the FoxO and Hippo signaling pathways by KEGG analysis. Functional validation of two miRNA in HASM showed effects on cell growth and diameter. CONCLUSION: Reduced circulatory miRNA expression at baseline is associated with an increase in PC20. These miRNA provide biologic insights into, and may serve as biomarkers of, asthma severity. miR-16-5p and -30d-5p regulate airway smooth muscle phenotypes critically involved in asthma pathogenesis, supporting a mechanistic link to these findings. Functional ASM phenotypes may be directly relevant to AHR.

Sex-Based Genetic Association Study Identifies CELSR1 as a Possible Chronic Obstructive Pulmonary Disease Risk Locus among Women.

Chronic obstructive pulmonary disease (COPD) is a complex disease with strong environmental and genetic influences and sexually dimorphic features. Although genetic risk factors for COPD have been identified, much of the heritability remains unexplained. Sex-based genetic association studies may uncover additional COPD genetic risk factors. We studied current and former smokers from COPD case-control cohorts (COPDGene non-Hispanic whites and African Americans, Evaluation of COPD Longitudinally to Identify Predictive Surrogate End-Points, and Genetics of Chronic Obstructive Lung Disease). COPD was defined as post-bronchodilator forced expiratory volume in 1 second/forced vital capacity less than 0.70 and forced expiratory volume in 1 second percent predicted less than 80. Testing was performed across all cohorts and combined in a meta-analysis adjusted for age, pack-years, and genetic ancestry. We first performed genome-wide single-nucleotide polymorphism (SNP)-by-sex interaction testing on the outcome of COPD affection status. We performed sex-stratified association testing for SNPs with interaction P less than 10-6. We examined over 8 million SNPs in four populations, including 6,260 subjects with COPD (40.6% female) and 5,269 smoking control subjects (47.3% female). The SNP rs9615358 in the cadherin gene CELSR1 approached genome-wide significance for an interaction with sex (P = 1.24 × 10-7). In the sex-stratified meta-analysis, this SNP was associated with COPD among females (odds ratio, 1.37 [95% confidence interval, 1.25-1.49]; P = 3.32 × 10-7) but not males (odds ratio, 0.90 [95% confidence interval, 0.79-1.01]; P = 0.06). CELSR1 is involved in fetal lung development. In a human fetal lung tissue dataset, we observed greater CELSR1 expression in female compared with male samples. This SNP-by-sex genome-wide association analysis identified the fetal lung development gene, CELSR1, as a potential sex-specific risk factor for COPD. Identifying sex-specific genetic risk factors may reveal new insights into sexually dimorphic features of COPD.

Circulating MicroRNAs: Association with Lung Function in Asthma.

BACKGROUND: MicroRNAs are key transcriptional and network regulators previously associated with asthma susceptibility. However, their role in relation to asthma severity has not been delineated. OBJECTIVE: We hypothesized that circulating microRNAs could serve as biomarkers of changes in lung function in asthma patients. METHODS: We isolated microRNAs from serum samples obtained at randomization for 160 participants of the Childhood Asthma Management Program. Using a TaqMan microRNA array containing 754 microRNA primers, we tested for the presence of known asthma microRNAs, and assessed the association of the individual microRNAs with lung function as measured by FEV1/FVC, FEV1% and FVC%. We further tested the subset of FEV1/FVC microRNAs for sex-specific and lung developmental associations. RESULTS: Of the 108 well-detected circulating microRNAs, 74 (68.5%) had previously been linked to asthma susceptibility. We found 22 (20.3%), 4 (3.7%) and 8 (7.4%) microRNAs to be associated with FEV1/FVC, FEV1% and FVC%, respectively. 8 (of 22) FEV1/FVC, 3 (of 4) FEV1% and 1 (of 8) FVC% microRNAs had functionally validated target genes that have been linked via genome wide association studies to asthma and FEV1 change. Among the 22 FEV1/FVC microRNAs, 9 (40.9%) remain associated with FEV1/FVC in boys alone in a sex-stratified analysis (compared with 3 FEV1/FVC microRNAs in girls alone), 7 (31.8%) were associated with fetal lung development, and 3 (13.6%) in both. Ontology analyses revealed enrichment for pathways integral to asthma, including PPAR signaling, G-protein coupled signaling, actin and myosin binding, and respiratory system development. CONCLUSIONS: Circulating microRNAs reflect asthma biology and are associated with lung function differences in asthmatics. They may represent biomarkers of asthma severity.

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.

Risk loci for chronic obstructive pulmonary disease: a genome-wide association study and meta-analysis.

BACKGROUND: The genetic risk factors for susceptibility to chronic obstructive pulmonary disease (COPD) are still largely unknown. Additional genetic variants are likely to be identified by genome-wide association studies in larger cohorts or specific subgroups. We sought to identify risk loci for moderate to severe and severe COPD with data from several cohort studies. METHODS: We combined genome-wide association analysis data from participants in the COPDGene study (non-Hispanic white and African-American ethnic origin) and the ECLIPSE, NETT/NAS, and Norway GenKOLS studies (self-described white ethnic origin). We did analyses comparing control individuals with individuals with moderate to severe COPD and with a subset of individuals with severe COPD. Single nucleotide polymorphisms yielding a p value of less than 5 × 10(-7) in the meta-analysis at loci not previously described were genotyped in individuals from the family-based ICGN study. We combined results in a joint meta-analysis (threshold for significance p<5 × 10(-8)). FINDINGS: Analysis of 6633 individuals with moderate to severe COPD and 5704 control individuals confirmed association at three known loci: CHRNA3 (p=6·38 × 10(-14)), FAM13A (p=1·12 × 10(-14)), and HHIP (p=1·57 × 10(-12)). We also showed significant evidence of association at a novel locus near RIN3 (p=5·25 × 10(-9)). In the overall meta-analysis (ie, including data from 2859 ICGN participants), the association with RIN3 remained significant (p=5·4 × 10(-9)). 3497 individuals were included in our analysis of severe COPD. The effect estimates for the loci near HHIP and CHRNA3 were significantly stronger in severe disease than in moderate to severe disease (p<0·01). We also identified associations at two additional loci: MMP12 (overall joint meta-analysis p=2·6 × 10(-9)) and TGFB2 (overall joint meta-analysis p=8·3 × 10(-9)). INTERPRETATION: We have confirmed associations with COPD at three known loci and identified three new genome-wide significant associations. Genetic variants other than in α-1 antitrypsin increase the risk of COPD. FUNDING: US National Heart, Lung, and Blood Institute; the Alpha-1 Foundation; the COPD Foundation through contributions from AstraZeneca, Boehringer Ingelheim, Novartis, and Sepracor; GlaxoSmithKline; Centers for Medicare and Medicaid Services; Agency for Healthcare Research and Quality; and US Department of Veterans Affairs.

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.

Mouse allergen, lung function, and atopy in Puerto Rican children.

OBJECTIVE: To examine the relation between mouse allergen exposure and asthma in Puerto Rican children. METHODS: Mus m 1, Der p 1, Bla g 2, and Fel d 1 allergens were measured in dust samples from homes of Puerto Rican children with (cases) and without (controls) asthma in Hartford, CT (n = 449) and San Juan (SJ), Puerto Rico (n = 678). Linear or logistic regression was used for the multivariate analysis of mouse allergen (Mus m 1) and lung function (FEV(1) and FEV(1)/FVC) and allergy (total IgE and skin test reactivity (STR) to ≥1 allergen) measures. RESULTS: Homes in SJ had lower mouse allergen levels than those in Hartford. In multivariate analyses, mouse allergen was associated with higher FEV(1) in cases in Hartford (+70.6 ml, 95% confidence interval (CI) = 8.6-132.7 ml, P = 0.03) and SJ (+45.1 ml, 95% CI =  -0.5 to 90.6 ml, P = 0.05). In multivariate analyses of controls, mouse allergen was inversely associated with STR to ≥1 allergen in non-sensitized children (odds ratio [OR] for each log-unit increment in Mus m 1 = 0.7, 95% CI = 0.5-0.9, P<0.01). In a multivariate analysis including all children at both study sites, each log-increment in mouse allergen was positively associated with FEV(1) (+28.3 ml, 95% CI = 1.4-55.2 ml, P = 0.04) and inversely associated with STR to ≥1 allergen (OR for each log-unit increment in Mus m 1 = 0.8, 95% CI = 0.6-0.9, P<0.01). CONCLUSIONS: Mouse allergen is associated with a higher FEV(1) and lower odds of STR to ≥1 allergen in Puerto Rican children. This may be explained by the allergen itself or correlated microbial exposures.

Vitamin D insufficiency and severe asthma exacerbations in Puerto Rican children.

RATIONALE: Vitamin D insufficiency (a serum 25(OH)D <30 ng/ml) has been associated with severe asthma exacerbations, but this could be explained by underlying racial ancestry or disease severity. Little is known about vitamin D and asthma in Puerto Ricans. OBJECTIVES: To examine whether vitamin D insufficiency is associated with severe asthma exacerbations in Puerto Rican children, independently of racial ancestry, atopy, and time outdoors. METHODS: A cross-sectional study was conducted of 560 children ages 6-14 years with (n = 287) and without (n = 273) asthma in San Juan, Puerto Rico. We measured plasma vitamin D and estimated the percentage of African racial ancestry among participants using genome-wide genotypic data. We tested whether vitamin D insufficiency is associated with severe asthma exacerbations, lung function, or atopy (greater than or equal to one positive IgE to allergens) using logistic or linear regression. Multivariate models were adjusted for African ancestry, time outdoors, atopy, and other covariates. MEASUREMENTS AND MAIN RESULTS: Vitamin D insufficiency was common in children with (44%) and without (47%) asthma. In multivariate analyses, vitamin D insufficiency was associated with higher odds of greater than or equal to one severe asthma exacerbation in the prior year (odds ratio [OR], 2.6; 95% confidence interval [CI], 1.5-4.9; P = 0.001) and atopy, and a lower FEV(1)/FVC in cases. After stratification by atopy, the magnitude of the association between vitamin D insufficiency and severe exacerbations was greater in nonatopic (OR, 6.2; 95% CI, 2-21.6; P = 0.002) than in atopic (OR, 2; 95% CI, 1-4.1; P = 0.04) cases. CONCLUSIONS: Vitamin D insufficiency is associated with severe asthma exacerbations in Puerto Rican children, independently of racial ancestry, atopy, or markers of disease severity or control.

African ancestry and lung function in Puerto Rican children.

BACKGROUND: Puerto Rican and African American subjects share a significant proportion of African ancestry. Recent findings suggest that African ancestry influences lung function in African American adults. OBJECTIVE: We sought to examine whether a greater proportion of African ancestry is associated with lower FEV(1) and forced vital capacity (FVC) in Puerto Rican children independently of socioeconomic status, health care access, or key environmental/lifestyle factors. METHODS: We performed a cross-sectional case-control study of 943 Puerto Rican children aged 6 to 14 years with (n= 520) and without (n= 423) asthma (defined as physician-diagnosed asthma and wheeze in the prior year) living in Hartford, Connecticut (n= 383), and San Juan, Puerto Rico (n= 560). We estimated the percentage of African racial ancestry in study participants using genome-wide genotypic data. We tested whether African ancestry is associated with FEV(1) and FVC using linear regression. Multivariate models were adjusted for indicators of socioeconomic status and health care and selected environmental/lifestyle exposures. RESULTS: After adjustment for household income and other covariates, each 20% increment in African ancestry was significantly associated with lower prebronchodilator FEV(1) (-105 mL; 95% CI, -159 to -51 mL; P< .001) and FVC (-133 mL; 95% CI, -197 to -69 mL; P< .001) and postbronchodilator FEV(1) (-152 mL; 95% CI, -210 to -94 mL; P< .001) and FVC (-145 mL; 95% CI, -211 to -79 mL; P< .001) in children with asthma. Similar but weaker associations were found for prebronchodilator and postbronchodilator FEV(1) (change for each 20% increment in African ancestry, -78 mL; 95% CI, -131 to -25 mL; P= .004) and for postbronchodilator FVC among children without asthma. CONCLUSIONS: Genetic factors, environmental/lifestyle factors, or both correlated with African ancestry might influence childhood lung function in Puerto Rican subjects.

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.

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.

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.

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.

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.

Analyses of shared genetic factors between asthma and obesity in children.

BACKGROUND: Epidemiologic studies consistently show associations between asthma and obesity. Shared genetics might account for this association. OBJECTIVE: We sought to identify genetic variants associated with both asthma and obesity. METHODS: On the basis of a literature search, we identified genes from (1) genome-wide association studies (GWASs) of body mass index (BMI; n = 17 genes), (2) GWASs of asthma (n = 14), and (3) candidate gene studies of BMI and asthma (n = 7). We used GWAS data from the Childhood Asthma Management Program to analyze associations between single nucleotide polymorphisms (SNPs) in these genes and asthma (n = 359 subjects) and BMI (n = 537). RESULTS: One top BMI GWAS SNP from the literature, rs10938397 near glucosamine-6-phosphate deaminase 2 (GNPDA2), was associated with both BMI (P = 4 x 10(-4)) and asthma (P = .03). Of the top asthma GWAS SNPs and the candidate gene SNPs, none was found to be associated with both BMI and asthma. Gene-based analyses that included all available SNPs in each gene found associations (P < .05) with both phenotypes for several genes: neuronal growth regulator 1 (NEGR1); roundabout, axon guidance receptor, homolog 1 (ROBO1); diacylglycerol kinase, gamma (DGKG); Fas apoptotic inhibitory molecule 2 (FAIM2); fat mass and obesity associated (FTO); and carbohydrate (N-acetylgalactosamine 4-0) sulfotransferase 8 (CHST8) among the BMI GWAS genes; interleukin 1 receptor-like 1 / interleukin 18 receptor 1 (IL1RL1/IL18R1), dipeptidyl-peptidase 10 (DPP10), phosphodiesterase 4D (PDE4D), V-myb myeloblastosis viral oncogene homolog (MYB), PDE10A, IL33, and especially protein tyrosine phosphatase, receptor type D (PTPRD) among the asthma GWAS genes; and protein kinase C, alpha (PRKCA) among the BMI and asthma candidate genes. CONCLUSIONS: SNPs within several genes showed associations to BMI and asthma at a genetic level, but none of these associations were significant after correction for multiple testing. Our analysis of known candidate genes reveals some evidence for shared genetics between asthma and obesity, but other shared genetic determinants are likely to be identified in novel loci.