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.
Genome-wide association analysis in COPDGene (non-Hispanic whites and
African-Americans) was combined with existing data from the ECLIPSE,
NETT/NAS, and GenKOLS (Norway) studies. Analyses were performed both using
all moderate-to-severe cases and the subset of severe cases. Top loci not
previously described as genome-wide significant were genotyped in the ICGN
study, and results combined in a joint meta-analysis.
Analysis of a total of 6,633 moderate-to-severe cases and 5,704
controls confirmed association at three known loci:
CHRNA3/CHRNA5/IREB2, FAM13A, and HHIP
(10−12 < P < 10−14),
and also showed significant evidence of association at a novel locus near
RIN3 (overall P, including ICGN =
5•4×10−9). In the severe COPD analysis
(n=3,497), the effects at two of three previously described loci were
significantly stronger; we also identified two additional loci previously
reported to affect gene expression of MMP12 and
TGFB2 (overall P = 2•6x10−9
and 8•3×10−9). RIN3 and
TGFB2 expression levels were reduced in a set of Lung
Tissue Research Consortium COPD lung tissue samples compared with
In a genome-wide study of COPD, we confirmed associations at three
known loci and found additional genome-wide significant associations with
moderate-to-severe COPD near RIN3 and with severe COPD near
MMP12 and TGFB2. Genetic variants,
apart from alpha-1 antitrypsin deficiency, increase the risk of COPD. Our
analysis of severe COPD suggests additional genetic variants may be
identified by focusing on this subgroup.
National Heart, Lung, and Blood Institute; 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; US Department of Veterans
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.
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.
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.
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.
Asthma; Airway hyperresponsiveness; Genome-wide association study; ITGB5; AGFG1
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 (FEV1). Though unadjusted linear models of FEV1 identified several genes strongly correlated with lung function, these correlations were due to racial differences in the distribution of both FEV1 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.
ancestry; gene expression; population stratification; self-identified race
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 FEV1/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.
vitamin D; asthma exacerbations; Puerto Ricans; childhood
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 FEV1 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 FEV1 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.
polymorphism; genome; pharmacogenomics; glucocorticoid
Puerto Ricans and African Americans share a significant proportion of African ancestry. Recent findings suggest that African ancestry influences lung function in African American adults.
To examine whether a greater proportion of African ancestry is associated with lower FEV1 and FVC in Puerto Rican children, independently of socioeconomic status (SES), healthcare access or key environmental/lifestyle (EL) factors.
Cross-sectional case-control study of 943 Puerto Rican children ages 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 (CT, n=383) and San Juan (PR, 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 FEV1 and FVC using linear regression. Multivariate models were adjusted for indicators of SES and healthcare, and selected EL exposures.
After adjustment for household income and other covariates, each 20% increment in African ancestry was significantly associated with lower pre-bronchodilator(BD) FEV1 (−105 ml, 95% confidence interval [CI] = −159 ml to −51 ml, P <0.001) and FVC (−133 ml, 95% CI −197 ml to −69 ml, P <0.001), and post-BD FEV1 (−152 ml, 95% CI=−210 ml to −94 ml, P <0.001) and FVC (−145 ml, 95% CI= −211 to −79 ml, P <0.001) in children with asthma. Similar but weaker associations were found for pre- and post-BD FEV1 (change for each 20% increment in African ancestry= −78 ml, 95% CI= −131 to −25 ml, P=0.004), and for post-BD FVC among children without asthma.
Genetic and/or EL factors correlated with African ancestry may influence childhood lung function in Puerto Ricans.
ancestry; FEV1; FVC; Puerto Ricans; childhood
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.
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.
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.
A functional GLCCI1 variant is associated with substantial decrements in the response to inhaled glucocorticoids in patients with asthma. (Funded by the National Institutes of Health and others; ClinicalTrials.gov number, NCT00000575.)
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 FEV1 (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.
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.
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.
chronic obstructive pulmonary disease genetics; chronic obstructive pulmonary disease epidemiology; chronic obstructive pulmonary disease metabolism; genome-wide association study
To examine the relation between mouse allergen exposure and asthma in Puerto Rican children.
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 (FEV1 and FEV1/FVC) and allergy (total IgE and skin test reactivity (STR) to ≥1 allergen) measures.
Homes in SJ had lower mouse allergen levels than those in Hartford. In multivariate analyses, mouse allergen was associated with higher FEV1 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 FEV1 (+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).
Mouse allergen is associated with a higher FEV1 and lower odds of STR to ≥1 allergen in Puerto Rican children. This may be explained by the allergen itself or correlated microbial exposures.
Bronchodilator response (BDR) is an important asthma phenotype that measures reversibility of airway obstruction by comparing lung function (i.e. FEV1) 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 FEV1 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.
Bronchodilator response (BDR) is an important asthma phenotype that measures reversibility of airway obstruction by comparing lung function before and after the administration of short-acting β2-agonists, common medications used for asthma treatment. We performed a genome-wide association study of BDR with 1,644 white asthmatic subjects from six drug clinical trials and attempted to replicate these findings in 1,051 white subjects from two independent cohorts. The most significant associated variant was near the SPATS2L gene. We knocked down SPATS2L mRNA in human airway smooth muscle cells and found that β2-adrenergic receptor levels increased, suggesting that SPATS2L may be a regulator of BDR. Our results highlight the promise of pursuing GWAS results that do not necessarily reach genome-wide significance and are an example of how results from pharmacogenetic GWAS can be studied functionally.
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.
Epidemiological studies consistently show associations between asthma and obesity. Shared genetics may account for this association.
To identify genetic variants associated with both asthma and obesity.
Based on a literature search, we identified genes from: 1) Genome-wide association studies (GWAS) of Body Mass Index (BMI) (n=17 genes), 2) GWAS 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 (CAMP) to analyze associations between single nucleotide polymorphisms (SNPs) in these genes and asthma (n=359 subjects) and BMI (n=537).
One top BMI GWAS SNP from the literature, rs10938397 near GNPDA2, was associated with both BMI (p=4 × 10−4) and asthma (p=0.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<0.05) with both phenotypes for several genes: NEGR1, ROBO1, DGKG, FAIM2, FTO and CHST8 among the BMI GWAS genes; ILRL1/IL18R1, DPP10, PDE4D, MYB, PDE10A, IL33 and especially PTPRD among the asthma GWAS genes; and PRKCA among the BMI and asthma candidate genes.
SNPs within several genes showed associations to BMI and asthma at a gene 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.
Association; Asthma; BMI; Children; Genetics; GWAS; Obesity; Polymorphism; SNP
Single nucleotide polymorphisms (SNPs) in thymic stromal lymphopoietin (TSLP) have been associated with IgE (in girls) and asthma (in general). We sought to determine whether TSLP SNPs are associated with asthma in a sex-specific fashion.
We conducted regular and sex-stratified analyses of association between SNPs in TSLP and asthma in families of asthmatic children in Costa Rica. Significant findings were replicated in white and African-American participants in the Childhood Asthma Management Program, in African Americans in the Genomic Research on Asthma in the African Diaspora study, in whites and Hispanics in the Children’s Health Study, and in whites in the Framingham Heart Study (FHS).
Two SNPs in TSLP (rs1837253 and rs2289276) were significantly associated with a reduced risk of asthma in combined analyses of all cohorts (p values of 2×10−5 and 1×10−5, respectively). In a sex-stratified analysis, the T allele of rs1837253 was significantly associated with a reduced risk of asthma in males only (p= 3×10−6). Alternately, the T allele of rs2289276 was significantly associated with a reduced risk of asthma in females only (p= 2×10−4). Findings for rs2289276 were consistent in all cohorts except the FHS.
TSLP variants are associated with asthma in a sex-specific fashion.
asthma; genetic association; sex-specific; thymic stromal lymphopoietin; TSLP
Substantial evidence suggests that there is genetic susceptibility to chronic obstructive pulmonary disease (COPD). To identify common genetic risk variants, we performed a genome-wide association study in 2940 cases and 1380 smoking controls with normal lung function. We demonstrate a novel susceptibility locus at 4q22.1 in FAM13A (rs7671167, OR=0.76, P=8.6×10−8) and provide evidence of replication in one case-control and two family-based cohorts (for all studies, combined P=1.2×10−11).
Rationale: Association studies have implicated many genes in asthma pathogenesis, with replicated associations between single-nucleotide polymorphisms (SNPs) and asthma reported for more than 30 genes. Genome-wide genotyping enables simultaneous evaluation of most of this variation, and facilitates more comprehensive analysis of other common genetic variation around these candidate genes for association with asthma.
Objectives: To use available genome-wide genotypic data to assess the reproducibility of previously reported associations with asthma and to evaluate the contribution of additional common genetic variation surrounding these loci to asthma susceptibility.
Methods: Illumina Human Hap 550Kv3 BeadChip (Illumina, San Diego, CA) SNP arrays were genotyped in 422 nuclear families participating in the Childhood Asthma Management Program. Genes with at least one SNP demonstrating prior association with asthma in two or more populations were tested for evidence of association with asthma, using family-based association testing.
Measurements and Main Results: We identified 39 candidate genes from the literature, using prespecified criteria. Of the 160 SNPs previously genotyped in these 39 genes, 10 SNPs in 6 genes were significantly associated with asthma (including the first independent replication for asthma-associated integrin β3 [ITGB3]). Evaluation of 619 additional common variants included in the Illumina 550K array revealed additional evidence of asthma association for 15 genes, although none were significant after adjustment for multiple comparisons.
Conclusions: We replicated asthma associations for a minority of candidate genes. Pooling genome-wide association study results from multiple studies will increase the power to appreciate marginal effects of genes and further clarify which candidates are true “asthma genes.”
asthma; replication; single-nucleotide polymorphism; integrin β3; association
Rationale: Maternal vitamin D intake during pregnancy has been inversely associated with asthma symptoms in early childhood. However, no study has examined the relationship between measured vitamin D levels and markers of asthma severity in childhood.
Objectives: To determine the relationship between measured vitamin D levels and both markers of asthma severity and allergy in childhood.
Methods: We examined the relation between 25-hydroxyvitamin D levels (the major circulating form of vitamin D) and markers of allergy and asthma severity in a cross-sectional study of 616 Costa Rican children between the ages of 6 and 14 years. Linear, logistic, and negative binomial regressions were used for the univariate and multivariate analyses.
Measurements and Main Results: Of the 616 children with asthma, 175 (28%) had insufficient levels of vitamin D (<30 ng/ml). In multivariate linear regression models, vitamin D levels were significantly and inversely associated with total IgE and eosinophil count. In multivariate logistic regression models, a log10 unit increase in vitamin D levels was associated with reduced odds of any hospitalization in the previous year (odds ratio [OR], 0.05; 95% confidence interval [CI], 0.004–0.71; P = 0.03), any use of antiinflammatory medications in the previous year (OR, 0.18; 95% CI, 0.05–0.67; P = 0.01), and increased airway responsiveness (a ≤8.58-μmol provocative dose of methacholine producing a 20% fall in baseline FEV1 [OR, 0.15; 95% CI, 0.024–0.97; P = 0.05]).
Conclusions: Our results suggest that vitamin D insufficiency is relatively frequent in an equatorial population of children with asthma. In these children, lower vitamin D levels are associated with increased markers of allergy and asthma severity.
Rationale: Polymorphisms in the gene for transforming growth factor-β1 (TGFB1) have been associated with asthma, but not with airway responsiveness or disease exacerbations in subjects with asthma.
Objectives: To test for association between single nucleotide polymorphisms (SNPs) in TGFB1 and markers of asthma severity in childhood.
Methods: We tested for the association between nine SNPs in TGFB1 and indicators of asthma severity (lung function, airway responsiveness, and disease exacerbations) in two cohorts: 416 Costa Rican parent-child trios and 465 families of non-Hispanic white children in the Childhood Asthma Management Program (CAMP). We also tested for the interaction between these polymorphisms and exposure to dust mite allergen on asthma severity.
Measurements and Main Results: The A allele of promoter SNP rs2241712 was associated with increased airway responsiveness in Costa Rica (P = 0.0006) and CAMP (P = 0.005), and the C allele of an SNP in the promoter region (rs1800469) was associated with increased airway responsiveness in both cohorts (P ≤ 0.01). Dust mite exposure modified the effect of the C allele of exonic SNP rs1800471 on airway responsiveness (P = 0.03 for interactions in both cohorts). The T allele of a coding SNP (rs1982073) was associated with a reduced risk of asthma exacerbations in Costa Rica (P = 0.009) and CAMP (P = 0.005). Dust mite exposure also significantly modified the effect of the A allele of the promoter SNP rs2241712 on asthma exacerbations in both cohorts.
Conclusions: SNPs in TGFB1 are associated with airway responsiveness and disease exacerbations in children with asthma. Moreover, dust mite exposure may modify the effect of TGFB1 SNPs on airway responsiveness and asthma exacerbations.
airway responsiveness; asthma; dust mite allergen; single nucleotide polymorphisms; transforming growth factor-β1
Although asthma is highly prevalent among certain Hispanic subgroups, genetic determinants of asthma and asthma‐related traits have not been conclusively identified in Hispanic populations. A study was undertaken to identify genomic regions containing susceptibility loci for pulmonary function and bronchodilator responsiveness (BDR) in Costa Ricans.
Eight extended pedigrees were ascertained through schoolchildren with asthma in the Central Valley of Costa Rica. Short tandem repeat (STR) markers were genotyped throughout the genome at an average spacing of 8.2 cM. Multipoint variance component linkage analyses of forced expiratory volume in 1 second (FEV1) and FEV1/ forced vital capacity (FVC; both pre‐bronchodilator and post‐bronchodilator) and BDR were performed in these eight families (pre‐bronchodilator spirometry, n = 640; post‐bronchodilator spirometry and BDR, n = 624). Nine additional STR markers were genotyped on chromosome 7. Secondary analyses were repeated after stratification by cigarette smoking.
Among all subjects, the highest logarithm of the odds of linkage (LOD) score for FEV1 (post‐bronchodilator) was found on chromosome 7q34–35 (LOD = 2.45, including the additional markers). The highest LOD scores for FEV1/FVC (pre‐bronchodilator) and BDR were found on chromosomes 2q (LOD = 1.53) and 9p (LOD = 1.53), respectively. Among former and current smokers there was near‐significant evidence of linkage to FEV1/FVC (post‐bronchodilator) on chromosome 5p (LOD = 3.27) and suggestive evidence of linkage to FEV1 on chromosomes 3q (pre‐bronchodilator, LOD = 2.74) and 4q (post‐bronchodilator, LOD = 2.66).
In eight families of children with asthma in Costa Rica, there is suggestive evidence of linkage to FEV1 on chromosome 7q34–35. In these families, FEV1/FVC may be influenced by an interaction between cigarette smoking and a locus (loci) on chromosome 5p.
Rationale: The basis for gender influences on allergen-specific IgEs is unclear.
Objectives: To perform regular and sex-stratified genomewide linkage analyses of IgE to each of three allergens (Ascaris lumbricoides, Blatella germanica [German cockroach]), and Dermatophagoides pteronyssinus [dust mite]) and to conduct an association study of a candidate gene in a linked genomic region.
Methods: Genomewide linkage analyses of allergen-specific IgEs were conducted in 653 members of eight large families of Costa Rican children with asthma. An analysis of the association between single-nucleotide polymorphisms in thymic stromal lymphopoietin (TSLP) and IgE measurements was conducted in 417 parent–child trios in Costa Rica. Significant results were replicated in 470 families of white children in the Childhood Asthma Management Program (CAMP).
Measurements and Main Results: Among all subjects, there was suggestive evidence of linkage (LOD ⩾ 2.72) to IgE to Ascaris (on chromosome 7q) and IgE to dust mite (on chromosomes 7p and 12q). In a sex-stratified analysis, there was significant evidence of linkage to IgE to cockroach on chromosome 5q23 (peak LOD, 4.14 at 127 cM) in female subjects. TSLP is located within the 1.5 LOD-unit support interval for this linkage peak and has female-specific effects on lung disease in mice. In a sex-stratified analysis, the T allele of single-nucleotide polymorphism rs2289276 in TSLP was associated with reductions in IgE to cockroach (in Costa Rican girls) and total IgE (in girls in Costa Rica and in CAMP; P value for sex-by-genotype interaction, <0.01 in both studies).
Conclusions: Consistent with findings in murine models, a variant in TSLP may have female-specific effects on allergic phenotypes.
immunoglobulin E; linkage; thymic stromal lymphopoietin; single-nucleotide polymorphisms
Rationale: Replication of gene-disease associations has become a requirement in complex trait genetics.
Objectives: In studies of childhood asthma from two different ethnic groups, we attempted to replicate associations with five potential asthma susceptibility genes previously identified by positional cloning.
Methods: We analyzed two family-based samples ascertained through an asthmatic proband: 497 European-American children from the Childhood Asthma Management Program and 439 Hispanic children from the Central Valley of Costa Rica. We genotyped 98 linkage disequilibrium–tagging single-nucleotide polymorphisms (SNPs) in five genes: ADAM33, DPP10, GPR154 (HUGO name: NPSR1), HLA-G, and the PHF11 locus (includes genes SETDB2 and RCBTB1). SNPs were tested for association with asthma and two intermediate phenotypes: airway hyperresponsiveness and total serum immunoglobulin E levels.
Measurements and Main Results: Despite differing ancestries, linkage disequilibrium patterns were similar in both cohorts. Of the five evaluated genes, SNP-level replication was found only for GPR154 (NPSR1). In this gene, three SNPs were associated with asthma in both cohorts, although the opposite alleles were associated in either study. Weak evidence for locus-level replication with asthma was found in the PHF11 locus, although there was no overlap in the associated SNP across the two cohorts. No consistent associations were observed for the three other genes.
Conclusions: These results provide some further support for the role of genetic variation in GPR154 (NPSR1) and PHF11 in asthma susceptibility and also highlight the challenges of replicating genetic associations in complex traits such as asthma, even for genes identified by linkage analysis.
bronchial hyperreactivity; immunoglobulin E; linkage disequilibrium; NPSR1; single-nucleotide polymorphism
Genetic epidemiology studies of end-stage lung disease are potentially hindered by low numbers of participants due to early death of patients from the underlying disease, or due to exclusion from studies after patients have had lung transplants, because of concern about bias of genotype data due to chimerism. The number of participants enrolled in genetic studies of end-stage lung disease could be increased by including those individuals who have undergone lung transplant. We hypothesized that individuals who have had lung transplants can be included in genetic epidemiology studies that use single nucleotide polymorphism and short tandem repeat marker data, without confounding due to chimerism. Ten probands with severe, early-onset chronic obstructive pulmonary disease were included in this analysis. Pre– and post–lung transplant DNA samples were used in the investigation of concordance of genotype results for 12 short tandem repeat markers and 23 single nucleotide polymorphisms. Concordance was observed for all genotypes before and after lung transplant. We conclude that the risk of biasing genetic epidemiology studies due to donor lung–related DNA microchimerism is low, and that the inclusion of post–lung transplantation participants will allow for larger genetic epidemiology studies of individuals with end-stage lung disease.
genetic epidemiology; lung; chimerism; transplantation
We describe a novel variant in the terminal exon of human elastin, c.2318 G>A, resulting in an amino acid substitution of glycine 773 to aspartate (G773D) in a pedigree with severe early-onset chronic obstructive pulmonary disease (COPD). Transfection studies with elastin cDNAs demonstrate that the glycine to aspartate change compromises the ability of the mutant protein to undergo normal elastin assembly. Other functional consequences of this amino acid substitution include altered proteolytic susceptibility of the C-terminal region of elastin and reduced interaction of the exon 36 sequence with matrix receptors on cells. These results suggest that the G773D variant confers structural and functional consequences relevant to the pathogenesis of COPD.
chronic obstructive pulmonary disease; elastin; extracellular matrix; genetics; mutation
Case-control studies have successfully identified many significant genetic associations for complex diseases, but lack of replication has been a criticism of case-control genetic association studies in general. We selected 12 candidate genes with reported associations to chronic obstructive pulmonary disease (COPD) and genotyped 29 polymorphisms in a family-based study and in a case-control study. In the Boston Early-Onset COPD Study families, significant associations with quantitative and/or qualitative COPD-related phenotypes were found for the tumor necrosis factor (TNF)-α −308G>A promoter polymorphism (P < 0.02), a coding variant in surfactant protein B (SFTPB Thr131Ile) (P = 0.03), and the (GT)31 allele of the heme oxygenase (HMOX1) promoter short tandem repeat (P = 0.02). In the case-control study, the SFTPB Thr131Ile polymorphism was associated with COPD, but only in the presence of a gene-by-environment interaction term (P = 0.01 for both main effect and interaction). The 30-repeat, but not the 31-repeat, allele of HMOX1 was associated (P = 0.04). The TNF −308G>A polymorphism was not significant. In addition, the microsomal epoxide hydrolase “fast” allele (EPHX1 His139Arg) was significantly associated in the case-control study (P = 0.03). Although some evidence for replication was found for SFTPB and HMOX1, none of the previously published COPD genetic associations was convincingly replicated across both study designs.
association studies; case-control studies; emphysema; genetics; single nucleotide polymorphism