Harmonizing Clinical Sequencing and Interpretation for the eMERGE III Network

The advancement of precision medicine requires new methods to coordinate and deliver genetic data from heterogeneous sources to physicians and patients. The eMERGE III Network enrolled >25,000 participants from biobank and prospective cohorts of predominantly healthy individuals for clinical genetic testing to determine clinically actionable findings. The network developed protocols linking together the 11 participant collection sites and 2 clinical genetic testing laboratories. DNA capture panels targeting 109 genes were used for testing of DNA and sample collection, data generation, interpretation, reporting, delivery, and storage were each harmonized. A compliant and secure network enabled ongoing review and reconciliation of clinical interpretations, while maintaining communication and data sharing between clinicians and investigators. A total of 202 individuals had positive diagnostic findings relevant to the indication for testing and 1,294 had additional/secondary findings of medical significance deemed to be returnable, establishing data return rates for other testing endeavors. This study accomplished integration of structured genomic results into multiple electronic health record (EHR) systems, setting the stage for clinical decision support to enable genomic medicine. Further, the established processes enable different sequencing sites to harmonize technical and interpretive aspects of sequencing tests, a critical achievement toward global standardization of genomic testing. The eMERGE protocols and tools are available for widespread dissemination.

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.

A genome-wide association study of bronchodilator response in asthmatics

Reversibility of airway obstruction in response to β2-agonists is highly variable among asthmatics, which is partially attributed to genetic factors. In a genome-wide association study of acute bronchodilator response (BDR) to inhaled albuterol, 534 290 single-nucleotide polymorphisms (SNPs) were tested in 403 white trios from the Childhood Asthma Management Program using five statistical models to determine the most robust genetic associations. The primary replication phase included 1397 polymorphisms in three asthma trials (pooled n=764). The second replication phase tested 13 SNPs in three additional asthma populations (n=241, n=215 and n=592). An intergenic SNP on chromosome 10, rs11252394, proximal to several excellent biological candidates, significantly replicated (P=1.98 × 10(-7)) in the primary replication trials. An intronic SNP (rs6988229) in the collagen (COL22A1) locus also provided strong replication signals (P=8.51 × 10(-6)). This study applied a robust approach for testing the genetic basis of BDR and identified novel loci associated with this drug response in asthmatics.

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

CONCLUSIONS

We have identified 2 SNPs associated with earlier onset of childhood asthma in 4 independent cohorts.

Collapse

Key Words

• asthma
• pediatrics
• age of onset
• asthma genetics
• c1orf100
• genome-wide association study
• pediatric asthma

Collapse

MESH Headings

• Adolescent
• Age of Onset
• Asthma/epidemiology
• Asthma/genetics
• Child
• Cohort Studies
• Costa Rica/epidemiology
• Female
• Genetic Predisposition to Disease
• Genome-Wide Association Study/methods
• Genotype
• Humans
• Linkage Disequilibrium
• Male
• North America/epidemiology
• Polymorphism, Single Nucleotide/genetics
• Sweden/epidemiology

Collapse

Grants

• R00 HL096840 NHLBI NIH HHS
• N01-HR-16051 NHLBI NIH HHS
• M01 RR002719 NCRR NIH HHS
• N01-HR-16044 NHLBI NIH HHS
• N01-HR-16048 NHLBI NIH HHS
• M01 RR000036 NCRR NIH HHS
• N01-HR-16045 NHLBI NIH HHS
• M01RR02719-14 NCRR NIH HHS
• M01 RR000051 NCRR NIH HHS
• HL66289 NHLBI NIH HHS
• R01 HL066289 NHLBI NIH HHS
• N01-HR-16047 NHLBI NIH HHS
• N01-HR-16046 NHLBI NIH HHS
• N02-HL-6-4278 NHLBI NIH HHS
• K01 HL004370 NHLBI NIH HHS
• N01-HR-16049 NHLBI NIH HHS
• RR00036 NCRR NIH HHS
• N01-HR-16052 NHLBI NIH HHS
• R37 HL066289 NHLBI NIH HHS
• N01-HR-16050 NHLBI NIH HHS
• Wellcome Trust
• WT084703MA Wellcome Trust
• N01HR16044 NHLBI NIH HHS
• M01RR00051 NCRR NIH HHS
• M01RR0099718-24 NCRR NIH HHS

Collapse

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

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.

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.

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.

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.

Regulatory haplotypes in ARG1 are associated with altered bronchodilator response

RATIONALE

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

Asthma-susceptibility variants identified using probands in case-control and family-based analyses

BACKGROUND

Asthma is a chronic respiratory disease whose genetic basis has been explored for over two decades, most recently via genome-wide association studies. We sought to find asthma-susceptibility variants by using probands from a single population in both family-based and case-control association designs.

METHODS

We used probands from the Childhood Asthma Management Program (CAMP) in two primary genome-wide association study designs: (1) probands were combined with publicly available population controls in a case-control design, and (2) probands and their parents were used in a family-based design. We followed a two-stage replication process utilizing three independent populations to validate our primary findings.

RESULTS

We found that single nucleotide polymorphisms with similar case-control and family-based association results were more likely to replicate in the independent populations, than those with the smallest p-values in either the case-control or family-based design alone. The single nucleotide polymorphism that showed the strongest evidence for association to asthma was rs17572584, which replicated in 2/3 independent populations with an overall p-value among replication populations of 3.5E-05. This variant is near a gene that encodes an enzyme that has been implicated to act coordinately with modulators of Th2 cell differentiation and is expressed in human lung.

CONCLUSIONS

Our results suggest that using probands from family-based studies in case-control designs, and combining results of both family-based and case-control approaches, may be a way to augment our ability to find SNPs associated with asthma and other complex diseases.

Predicting response to short-acting bronchodilator medication using Bayesian networks

AIMS

Bronchodilator response tests measure the effect of beta(2)-agonists, the most commonly used short-acting reliever drugs for asthma. We sought to relate candidate gene SNP data with bronchodilator response and measure the predictive accuracy of a model constructed with genetic variants.

MATERIALS & METHODS

Bayesian networks, multivariate models that are able to account for simultaneous associations and interactions among variables, were used to create a predictive model of bronchodilator response using candidate gene SNP data from 308 Childhood Asthma Management Program Caucasian subjects.

RESULTS

The model found that 15 SNPs in 15 genes predict bronchodilator response with fair accuracy, as established by a fivefold cross-validation area under the receiver-operating characteristic curve of 0.75 (standard error: 0.03).

CONCLUSION

Bayesian networks are an attractive approach to analyze large-scale pharmacogenetic SNP data because of their ability to automatically learn complex models that can be used for the prediction and discovery of novel biological hypotheses.

A role for Wnt signaling genes in the pathogenesis of impaired lung function in asthma

RATIONALE

Animal models demonstrate that aberrant gene expression in utero can result in abnormal pulmonary phenotypes.

OBJECTIVES

We sought to identify genes that are differentially expressed during in utero airway development and test the hypothesis that variants in these genes influence lung function in patients with asthma.

METHODS

Stage 1 (Gene Expression): Differential gene expression analysis across the pseudoglandular (n = 27) and canalicular (n = 9) stages of human lung development was performed using regularized t tests with multiple comparison adjustments. Stage 2 (Genetic Association): Genetic association analyses of lung function (FEV(1), FVC, and FEV(1)/FVC) for variants in five differentially expressed genes were conducted in 403 parent-child trios from the Childhood Asthma Management Program (CAMP). Associations were replicated in 583 parent-child trios from the Genetics of Asthma in Costa Rica study.

MEASUREMENTS AND MAIN RESULTS

Of the 1,776 differentially expressed genes between the pseudoglandular (gestational age: 7-16 wk) and the canalicular (gestational age: 17-26 wk) stages, we selected 5 genes in the Wnt pathway for association testing. Thirteen single nucleotide polymorphisms in three genes demonstrated association with lung function in CAMP (P < 0.05), and associations for two of these genes were replicated in the Costa Ricans: Wnt1-inducible signaling pathway protein 1 with FEV(1) (combined P = 0.0005) and FVC (combined P = 0.0004), and Wnt inhibitory factor 1 with FVC (combined P = 0.003) and FEV(1)/FVC (combined P = 0.003).

CONCLUSIONS

Wnt signaling genes are associated with impaired lung function in two childhood asthma cohorts. Furthermore, gene expression profiling of human fetal lung development can be used to identify genes implicated in the pathogenesis of lung function impairment in individuals with asthma.

IL1B polymorphisms modulate cystic fibrosis lung disease

RATIONALE

Variability in pulmonary disease severity is found in patients with cystic fibrosis (CF) who have identical mutations in the CF transmembrane conductance regulator (CFTR) gene. We hypothesized that one factor accounting for heterogeneity in pulmonary disease severity is variation in the family of genes affecting the biology of interleukin-1 (IL-1), which impacts acquisition and maintenance of Pseudomonas aeruginosa infection in animal models of chronic infection.

METHODS

We genotyped 58 single nucleotide polymorphisms (SNPs) in the IL-1 gene cluster in 808 CF subjects from the University of North Carolina and Case Western Reserve University (UNC/CWRU) joint cohort. All were homozygous for DeltaF508, and categories of "severe" (cases) or "mild" (control subjects) lung disease were defined by the lowest or highest quartile of forced expired volume (FEV(1)) for age in the CF population. After adjustment for age and gender, genotypic data were tested for association with lung disease severity. Odds ratios (ORs) comparing severe versus mild CF were also calculated for each genotype (with the homozygote major allele as the reference group) for all 58 SNPs. From these analyses, nine SNPs with a moderate effect size, OR < or =0.5 or >1.5, were selected for further testing. To replicate the case-control study results, we genotyped the same nine SNPs in a second population of CF parent-offspring trios (recruited from Children's Hospital Boston), in which the offspring had similar pulmonary phenotypes. For the trio analysis, both family-based and population-based associations were performed.

RESULTS

SNPs rs1143634 and rs1143639 in the IL1B gene demonstrated a consistent association with lung disease severity categories (P < 0.10) and longitudinal analysis of lung disease severity (P < 0.10) in CF in both the case-control and family-based studies. In females, there was a consistent association (false discovery rate adjusted joint P-value <0.06 for both SNPs) in both the analysis of lung disease severity in the UNC/CWRU cohort and the family-based analysis of affection status.

CONCLUSION

Our findings suggest that IL1beta is a clinically relevant modulator of CF lung disease.

5-lipoxygenase pharmacogenetics in asthma: overlap with Cys-leukotriene receptor antagonist loci

The two classes of leukotriene modifiers work by inhibiting different portions of the same pathway. We hypothesized that single nucleotide polymorphisms (SNPs) in genes associated with response to montelukast (a cys-leukotriene receptor antagonist) would also be associated with response to zileuton (a 5-lipoxygenase inhibitor). We genotyped 26 SNPs that had previously been interrogated for association with montelukast response in five candidate genes (ABCC1, ALOX5, CYSLTR1, LTA4H, LTC4S) in a population of 577 asthmatics who participated in a clinical trial comparing intermittent and continuous-release zileuton to placebo. After adjusting for age and sex, six SNPs in three genes were associated with longitudinal forced expiratory volume at 1 s in response to zileuton (P values 0.005-0.05). After adjusting for age and sex, six SNPs in three genes were associated with longitudinal forced expiratory volume at 1 s in response to zileuton (P values 0.005-0.05), including two SNPs (ALOX5 rs2115819 and ABCC1 rs119774) that we had previously reported as associated with FEV1 response to montelukast. Thus, the lung function response to zileuton is modulated by several of the loci that also influence montelukast response.

ARG1 is a novel bronchodilator response gene: screening and replication in four asthma cohorts

RATIONALE

Inhaled beta-agonists are one of the most widely used classes of drugs for the treatment of asthma. However, a substantial proportion of patients with asthma do not have a favorable response to these drugs, and identifying genetic determinants of drug response may aid in tailoring treatment for individual patients.

OBJECTIVES

To screen variants in candidate genes in the steroid and beta-adrenergic pathways for association with response to inhaled beta-agonists.

METHODS

We genotyped 844 single nucleotide polymorphisms (SNPs) in 111 candidate genes in 209 children and their parents participating in the Childhood Asthma Management Program. We screened the association of these SNPs with acute response to inhaled beta-agonists (bronchodilator response [BDR]) using a novel algorithm implemented in a family-based association test that ranked SNPs in order of statistical power. Genes that had SNPs with median power in the highest quartile were then taken for replication analyses in three other asthma cohorts.

MEASUREMENTS AND MAIN RESULTS

We identified 17 genes from the screening algorithm and genotyped 99 SNPs from these genes in a second population of patients with asthma. We then genotyped 63 SNPs from four genes with significant associations with BDR, for replication in a third and fourth population of patients with asthma. Evidence for association from the four asthma cohorts was combined, and SNPs from ARG1 were significantly associated with BDR. SNP rs2781659 survived Bonferroni correction for multiple testing (combined P value = 0.00048, adjusted P value = 0.047).

CONCLUSIONS

These findings identify ARG1 as a novel gene for acute BDR in both children and adults with asthma.

On the analysis of copy-number variations in genome-wide association studies: a translation of the family-based association test

Though there is an increasing support for an important contribution of copy number variation (CNV) to the genetic architecture of complex disease, few methods have been developed for the analysis of such variation in the context of genetic association studies. In this paper, we propose a generalization of family-based association tests (FBATs) to allow for the analysis of CNVs at a genome-wide level. We translate the popular FBAT approach so that, instead of genotypes, raw intensity values that reflect copy number are used directly in the test statistic, thereby bypassing the need for a CNV genotyping algorithm. Moreover, both inherited and de novo CNVs can be analyzed without any prior knowledge about the type of CNV, making it easily applicable to large-scale association studies. All robustness properties of the genotype FBAT approach are maintained and all previously developed FBAT extensions, including FBATs for time-to-onset, multivariate FBATs, and FBAT-testing strategies, can be directly transferred to the analysis of CNVs. Using simulation studies, we evaluate the power and the robustness of the new approach. Furthermore, for those CNVs that can be genotyped, we compare FBATs based on genotype calls with FBATs that are directly based on the intensity data. An application to one of the first CNV genome-wide-association studies of asthma identifies a very plausible candidate gene. A software implementation of the approach is freely available at http://www.hsph.harvard.edu/research/iuliana-ionita/software. The approach has also been completely integrated in the PBAT software package.

On the replication of genetic associations: timing can be everything!

The failure of researchers to replicate genetic-association findings is most commonly attributed to insufficient statistical power, population stratification, or various forms of between-study heterogeneity or environmental influences.(1) Here, we illustrate another potential cause for nonreplications that has so far not received much attention in the literature. We illustrate that the strength of a genetic effect can vary by age, causing "age-varying associations." If not taken into account during the design and the analysis of a study, age-varying genetic associations can cause nonreplication. By using the 100K SNP scan of the Framingham Heart Study, we identified an age-varying association between a SNP in ROBO1 and obesity and hypothesized an age-gene interaction. This finding was followed up in eight independent samples comprising 13,584 individuals. The association was replicated in five of the eight studies, showing an age-dependent relationship (one-sided combined p = 3.92 x 10(-9), combined p value from pediatric cohorts = 2.21 x 10(-8), combined p value from adult cohorts = 0.00422). Furthermore, this study illustrates that it is difficult for cross-sectional study designs to detect age-varying associations. If the specifics of age- or time-varying genetic effects are not considered in the selection of both the follow-up samples and in the statistical analysis, important genetic associations may be missed.

A common mitochondrial haplogroup is associated with elevated total serum IgE levels

BACKGROUND

Maternal history of asthma or atopy is among the most consistently reported risk factors for asthma and atopy in children, yet the molecular basis for this observation remains unclear. Mitochondria are inherited exclusively through the maternal line, raising the possibility that sequence variation in the mitochondrial genome contributes to the pathogenesis of asthma and atopy.

OBJECTIVE

We set out to determine whether common European mitochondrial haplogroups are associated with asthma-related atopic phenotypes.

METHODS

We studied 654 self-reported white children age 5 to 12 years with mild to moderate asthma participating in the Childhood Asthma Management Program. Eight haplogroup-tagging polymorphisms were genotyped with TaqMan probe hybridization assays, and mitochondrial haplogroup tests of association with asthma-related and atopy-related phenotypes were performed with haplo.stats.

RESULTS

We found significant evidence of mitochondrial haplogroup association with total serum IgE levels (global significance, P = .04), with carriers of European haplogroup U (frequency 11%) having higher total serum IgE levels (median level, 684 IU/L) compared with noncarriers (389 IU/L; P = .001). Haplogroup U carriers also had trends of greater skin prick test reactivity (P = .03) and higher frequency of atopic dermatitis (P = .07), although global haplogroup tests for these later 2 phenotypes were not significant at an alpha level of 0.05.

CONCLUSION

These data are the first to suggest that common mitochondrial haplogroups influence the atopic diathesis.

CLINICAL IMPLICATIONS

These findings may provide a molecular explanation for the prominent influence of maternal history of atopy on the development of atopy in offspring.

The association of a SNP upstream of INSIG2 with body mass index is reproduced in several but not all cohorts

A SNP upstream of the INSIG2 gene, rs7566605, was recently found to be associated with obesity as measured by body mass index (BMI) by Herbert and colleagues. The association between increased BMI and homozygosity for the minor allele was first observed in data from a genome-wide association scan of 86,604 SNPs in 923 related individuals from the Framingham Heart Study offspring cohort. The association was reproduced in four additional cohorts, but was not seen in a fifth cohort. To further assess the general reproducibility of this association, we genotyped rs7566605 in nine large cohorts from eight populations across multiple ethnicities (total n = 16,969). We tested this variant for association with BMI in each sample under a recessive model using family-based, population-based, and case-control designs. We observed a significant (p < 0.05) association in five cohorts but saw no association in three other cohorts. There was variability in the strength of association evidence across examination cycles in longitudinal data from unrelated individuals in the Framingham Heart Study Offspring cohort. A combined analysis revealed significant independent validation of this association in both unrelated (p = 0.046) and family-based (p = 0.004) samples. The estimated risk conferred by this allele is small, and could easily be masked by small sample size, population stratification, or other confounders. These validation studies suggest that the original association is less likely to be spurious, but the failure to observe an association in every data set suggests that the effect of SNP rs7566605 on BMI may be heterogeneous across population samples.

Obesity is an epidemic in the United States of America and developing world, portending an epidemic of related diseases such as diabetes and heart disease. While diet and lifestyle contribute to obesity, half of the population variation in body mass index, a common measure of obesity, is determined by inherited factors. Many studies have reported that common sequence variants in genes are associated with an increased risk for obesity, yet most of these are not reproducible in other study cohorts, suggesting that some are false. Recently, Herbert et al. reported a slightly increased risk of obesity for people carrying two copies of the minor allele at a common variant near INSIG2. We present our attempts to further evaluate this potential association with obesity in additional populations. We find evidence of increased risk of obesity for people carrying two copies of the minor allele in five out of nine cohorts tested, using both family- and population-based testing. We indicate possible reasons for the varied results, with the hope of encouraging a combined analysis across study cohorts to more precisely define the effect of this INSIG2 gene variant.

Sequence, haplotype, and association analysis of ADRbeta2 in a multiethnic asthma case-control study

RATIONALE

The comprehensive evaluation of gene variation, haplotype structure, and linkage disequilibrium is important in understanding the function of beta2-adrenergic receptor gene (ADRbeta2) on disease susceptibility, pulmonary function, and therapeutic responses in different ethnic groups with asthma.

OBJECTIVES

To identify ADRbeta2 polymorphisms and haplotype structure in white and African American subjects and to test for genotype and haplotype association with asthma phenotypes.

METHODS

A 5.3-kb region of ADRbeta2 was resequenced in 669 individuals from 429 whites and 240 African Americans. A total of 12 polymorphisms, representing an optimal haplotype tagging set, were genotyped in whites (338 patients and 326 control subjects) and African Americans (222 patients and 299 control subjects).

RESULTS

A total of 49 polymorphisms were identified, 21 of which are novel; 31 polymorphisms (frequency > 0.03) were used to identify 24 haplotypes (frequency > 0.01) and assess linkage disequilibrium. Association with ratio (FEV1/FVC)2 for single-nucleotide polymorphism +79 (p < 0.05) was observed in African Americans. Significant haplotype association for (FEV1/FVC)2 was also observed in African Americans.

CONCLUSIONS

There are additional genetic variants besides +46 (Gly16Arg) that are important in determining asthma phenotypes. These data suggest that the length of a poly-C repeat (+1269) in the 3' untranslated region of ADRbeta2 may influence lung function, and may be important in delineating variation in beta-agonist responses, especially in African Americans.