A pro-asthmatic culture milieu and β2-agonist (isoproterenol) were previously shown to regulate the expression of select transcription factors (TFs) within human airway epithelial and smooth muscle cells. This study tests 1116 single nucleotide polymorphisms (SNPs) across 98 of these TF genes for association with bronchodilator response (BDR) in asthma patients. Genotyping was conducted using the Illumina HumanHap550v3 Beadchip in 403 non-Hispanic White asthmatic children and their parents. SNPs were evaluated for association with BDR using family and population-based analyses. Forty-two SNPs providing p values < 0.1 in both analyses were then genotyped in three adult asthma trials. One SNP 5’ of the thyroid hormone receptor beta gene was associated with BDR in the childhood population and two adult populations (p value = 0.0012). This investigation identified a novel locus for inter-individual variability in BDR and represents a translation of a cellular drug-response study to potential personalization of clinical asthma management.
Bronchodilator response; transcription factor; association; thyroid hormone receptor β; asthma; pharmacogenetics
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.
Multiple intergenic single-nucleotide polymorphisms (SNPs) near hedgehog interacting protein (HHIP) on chromosome 4q31 have been strongly associated with pulmonary function levels and moderate-to-severe chronic obstructive pulmonary disease (COPD). However, whether the effects of variants in this region are related to HHIP or another gene has not been proven. We confirmed genetic association of SNPs in the 4q31 COPD genome-wide association study (GWAS) region in a Polish cohort containing severe COPD cases and healthy smoking controls (P = 0.001 to 0.002). We found that HHIP expression at both mRNA and protein levels is reduced in COPD lung tissues. We identified a genomic region located ∼85 kb upstream of HHIP which contains a subset of associated SNPs, interacts with the HHIP promoter through a chromatin loop and functions as an HHIP enhancer. The COPD risk haplotype of two SNPs within this enhancer region (rs6537296A and rs1542725C) was associated with statistically significant reductions in HHIP promoter activity. Moreover, rs1542725 demonstrates differential binding to the transcription factor Sp3; the COPD-associated allele exhibits increased Sp3 binding, which is consistent with Sp3's usual function as a transcriptional repressor. Thus, increased Sp3 binding at a functional SNP within the chromosome 4q31 COPD GWAS locus leads to reduced HHIP expression and increased susceptibility to COPD through distal transcriptional regulation. Together, our findings reveal one mechanism through which SNPs upstream of the HHIP gene modulate the expression of HHIP and functionally implicate reduced HHIP gene expression in the pathogenesis of COPD.
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.
It has recently been shown that vitamin D deficiency can increase asthma development and severity and that variations in vitamin D receptor genes are associated with asthma susceptibility.
We sought to find genetic factors that might interact with vitamin D levels to affect the risk of asthma exacerbation. Methods: We conducted a genome-wide study of gene–vitamin D interaction on asthma exacerbations using population-based and family-based approaches on 403 subjects and trios from the Childhood Asthma Management Program. Twenty-three polymorphisms with significant interactions were studied in a replication analysis in 584 children from a Costa Rican cohort. Results: We identified 3 common variants in the class I MHC–restricted T cell–associated molecule gene (CRTAM) that were associated with an increased rate of asthma exacerbations based on the presence of a low circulating vitamin D level. These results were replicated in a second independent population (unadjusted combined interaction, P =.00028–.00097; combined odds ratio, 3.28–5.38). One variant, rs2272094, is a nonsynonymous coding polymorphism of CRTAM. Functional studies on cell lines confirmed the interaction of vitamin D and rs2272094 on CRTAM expression. CRTAM is highly expressed in activated human CD8+ and natural killer T cells, both of which have been implicated in asthmatic patients.
The findings highlight an important gene-environment interaction that elucidates the role of vitamin D and CD8+ and natural killer T cells in asthma exacerbation in a genome-wide gene-environment interaction study that has been replicated in an independent population. The results suggest the potential importance of maintaining adequate vitamin D levels in subsets of high-risk asthmatic patients.
Gene-environment interaction; genome-wide association study; vitamin D; asthma exacerbation
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
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.
Rationale: Chromosome 12p has been linked to chronic obstructive pulmonary disease (COPD) in the Boston Early-Onset COPD Study (BEOCOPD), but a susceptibility gene in that region has not been identified.
Objectives: We used high-density single-nucleotide polymorphism (SNP) mapping to implicate a COPD susceptibility gene and an animal model to determine the potential role of SOX5 in lung development and COPD.
Methods: On chromosome 12p, we genotyped 1,387 SNPs in 386 COPD cases from the National Emphysema Treatment Trial and 424 control smokers from the Normative Aging Study. SNPs with significant associations were then tested in the BEOCOPD study and the International COPD Genetics Network. Based on the human results, we assessed histology and gene expression in the lungs of Sox5−/− mice.
Measurements and Main Results: In the case-control analysis, 27 SNPs were significant at P ≤ 0.01. The most significant SNP in the BEOCOPD replication was rs11046966 (National Emphysema Treatment Trial–Normative Aging Study P = 6.0 × 10−4, BEOCOPD P = 1.5 × 10−5, combined P = 1.7 × 10−7), located 3′ to the gene SOX5. Association with rs11046966 was not replicated in the International COPD Genetics Network. Sox5−/− mice showed abnormal lung development, with a delay in maturation before the saccular stage, as early as E16.5. Lung pathology in Sox5−/− lungs was associated with a decrease in fibronectin expression, an extracellular matrix component critical for branching morphogenesis.
Conclusions: Genetic variation in the transcription factor SOX5 is associated with COPD susceptibility. A mouse model suggests that the effect may be due, in part, to its effects on lung development and/or repair processes.
chronic obstructive pulmonary disease; emphysema; knockout mice; lung development; single nucleotide polymorphism
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.
Rationale: Chitinases are enzymes that cleave chitin, which is present in fungal cells. Two types of human chitinases, chitotriosidase and acidic mammalian chitinase, and the chitinase-like protein, YKL-40, seem to play an important role in asthma. We hypothesized that exposure to environmental fungi may modulate the effect of chitinases in individuals with asthma.
Objectives: To explore whether interactions between high fungal exposure and common genetic variants in the two chitinases in humans, CHIT1 and CHIA, and the chitinase 3-like 1 gene, CHI3L1, are associated with severe asthma exacerbations and other asthma-related outcomes.
Methods: Forty-eight single nucleotide polymorphisms (SNPs) in CHIT1, CHIA, and CHI3L1 and one CHIT1 duplication were genotyped in 395 subjects and their parents as part of the Childhood Asthma Management Program. Household levels of mold (an index of fungal exposure) were determined on house dust samples. We conducted family-based association tests with gene–environment interactions. Our outcome was severe exacerbation, defined as emergency department visits and hospitalizations from asthma over a 4-year period, and our secondary outcomes included indices of lung function and allergy-related phenotypes.
Measurements and Main Results: Of the 395 subjects who had mold levels at randomization, 24% (95 subjects) had levels that were greater than 25,000 units per gram of house dust (high mold exposure). High mold exposure significantly modified the relation between three SNPs in CHIT1 (rs2486953, rs4950936, and rs1417149) and severe exacerbations (P for interaction 0.0010 for rs2486953, 0.0008 for rs4950936, and 0.0005 for rs1417149). High mold exposure did not significantly modify the relationship between any of the other variants and outcomes.
Conclusions: Environmental exposure to fungi, modifies the effect of CHIT1 SNPs on severe asthma exacerbations.
chitinase; asthma; CHIA; CHIT1; CHI13L1
Rationale: Several family-based studies have identified genetic linkage for lung function and airflow obstruction to chromosome 2q.
Objectives: We hypothesized that merging results of high-resolution single nucleotide polymorphism (SNP) mapping in four separate populations would lead to the identification of chronic obstructive pulmonary disease (COPD) susceptibility genes on chromosome 2q.
Methods: Within the chromosome 2q linkage region, 2,843 SNPs were genotyped in 806 COPD cases and 779 control subjects from Norway, and 2,484 SNPs were genotyped in 309 patients with severe COPD from the National Emphysema Treatment Trial and 330 community control subjects. Significant associations from the combined results across the two case-control studies were followed up in 1,839 individuals from 603 families from the International COPD Genetics Network (ICGN) and in 949 individuals from 127 families in the Boston Early-Onset COPD Study.
Measurements and Main Results: Merging the results of the two case-control analyses, 14 of the 790 overlapping SNPs had a combined P < 0.01. Two of these 14 SNPs were consistently associated with COPD in the ICGN families. The association with one SNP, located in the gene XRCC5, was replicated in the Boston Early-Onset COPD Study, with a combined P = 2.51 × 10−5 across the four studies, which remains significant when adjusted for multiple testing (P = 0.02). Genotype imputation confirmed the association with SNPs in XRCC5.
Conclusions: By combining data from COPD genetic association studies conducted in four independent patient samples, we have identified XRCC5, an ATP-dependent DNA helicase, as a potential COPD susceptibility gene.
emphysema; genetic linkage; metaanalysis; single nucleotide polymorphism
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
Despite the potential role of chitinases and chitinase-like proteins in the pathogenesis of asthma, variants in their respective genes are not associated with asthma, changes in lung physiology or allergy-related phenotypes in Caucasian children.
chitinase; chitinase-like protein; CHIA; CHIT1; CHI13L1; YKL-40; AMCase; SNPs; asthma
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
Chronic obstructive pulmonary disease (COPD) is a heterogeneous syndrome, including emphysema and airway disease. Phenotypes defined on the basis of chest computed tomography (CT) may decrease disease heterogeneity and aid in the identification of candidate genes for COPD subtypes. To identify these genes, we performed genome-wide linkage analysis in extended pedigrees from the Boston Early-Onset COPD Study, stratified by emphysema status (defined by chest CT scans) of the probands, followed by genetic association analysis of positional candidate genes. A region on chromosome 1p showed strong evidence of linkage to lung function traits in families of emphysema-predominant probands in the stratified analysis (LOD score = 2.99 in families of emphysema-predominant probands versus 1.98 in all families). Association analysis in 949 individuals from 127 early-onset COPD pedigrees revealed association for COPD-related traits with an intronic single-nucleotide polymorphism (SNP) in transforming growth factor-β receptor-3 (TGFBR3) (P = 0.005). This SNP was significantly associated with COPD affection status comparing 389 cases from the National Emphysema Treatment Trial to 472 control smokers (P = 0.04), and with FEV1 (P = 0.004) and CT emphysema (P = 0.05) in 3,117 subjects from the International COPD Genetics Network. Gene-level replication of association with lung function was seen in 427 patients with COPD from the Lung Health Study. In conclusion, stratified linkage analysis followed by association testing identified TGFBR3 (betaglycan) as a potential susceptibility gene for COPD. Published human microarray and murine linkage studies have also demonstrated the importance of TGFBR3 in emphysema and lung function, and our group and others have previously found association of COPD-related traits with TGFB1, a ligand for TGFBR3.
betaglycan; chronic obstructive pulmonary disease; computed tomography; linkage; single nucleotide polymorphism
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).
Genetic variants influencing lung function in children and adults may ultimately lead to the development of chronic obstructive pulmonary disease (COPD), particularly in high-risk groups.
We tested for an association between single-nucleotide polymorphisms (SNPs) in the gene encoding matrix metalloproteinase 12 (MMP12) and a measure of lung function (prebronchodilator forced expiratory volume in 1 second [FEV1]) in more than 8300 subjects in seven cohorts that included children and adults. Within the Normative Aging Study (NAS), a cohort of initially healthy adult men, we tested for an association between SNPs that were associated with FEV1 and the time to the onset of COPD. We then examined the relationship between MMP12 SNPs and COPD in two cohorts of adults with COPD or at risk for COPD.
The minor allele (G) of a functional variant in the promoter region of MMP12 (rs2276109 [−82A→G]) was positively associated with FEV1 in a combined analysis of children with asthma and adult former and current smokers in all cohorts (P=2×10−6). This allele was also associated with a reduced risk of the onset of COPD in the NAS cohort (hazard ratio, 0.65; 95% confidence interval [CI], 0.46 to 0.92; P = 0.02) and with a reduced risk of COPD in a cohort of smokers (odds ratio, 0.63; 95% CI, 0.45 to 0.88; P = 0.005) and among participants in a family-based study of early-onset COPD (P = 0.006).
The minor allele of a SNP in MMP12 (rs2276109) is associated with a positive effect on lung function in children with asthma and in adults who smoke. This allele is also associated with a reduced risk of COPD in adult smokers.
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
The destruction of elastic fibers has been implicated in the pathogenesis of chronic obstructive pulmonary disease (COPD). Emphysema has been described in autosomal dominant cutis laxa, which can be caused by mutations in the elastin gene. Previously, a rare functional mutation in the terminal exon of elastin was found in a case of severe, early-onset COPD. To test the hypothesis that other similar elastin mutations may predispose to COPD, we screened 90 probands from the Boston Early-Onset COPD Study and 90 smoking control subjects from the Normative Aging Study for mutations in elastin exons using high-resolution DNA melt analysis followed by resequencing. Rare nonsynonymous single-nucleotide polymorphisms (SNPs) seen only in cases were examined for segregation with airflow obstruction within pedigrees. Common nonsynonymous SNPs were tested for association with COPD in a family-based analysis of 949 subjects from the Boston Early-Onset COPD Study, and in a case–control analysis in 389 COPD cases from the National Emphysema Treatment Trial and 472 control subjects from the Normative Aging Study. Of 28 elastin variants found, 3 were nonsynonymous SNPs found only in cases. The previously described Gly773Asp mutation was found in another proband. The other two SNPs did not clearly segregate with COPD within families. Two common nonsynonymous SNPs did not demonstrate significant associations in either a family-based or case–control analysis. Exonic SNPs in the elastin gene do not appear to be common risk factors for severe COPD.
elastin; chronic obstructive pulmonary disease; emphysema; genetic polymorphism
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
High-resolution melting curve analysis is an accurate method for mutation detection in genomic DNA. However, performance in whole genome amplified (WGA) versus genomic DNA has not been directly compared.
23 amplicons from 9 genes were PCR amplified in 39 paired genomic and WGA samples and analyzed by high-resolution melting curve analysis using the 96-well LightScanner® (Idaho Technology). Genotyping and bidirectional resequencing were used to verify melting curve results.
Melting patterns were concordant between the genomic and WGA samples in 823/863 (95%) of analyzed sample pairs. Of the discordant patterns, there was an overrepresentation of alternate melting curve patterns in the WGA samples suggesting the presence of a mutation (false positives). Targeted resequencing in 135 genomic and 136 WGA samples revealed 43 single nucleotide polymorphisms (SNPs). All SNPs detected in genomic samples were also detected in WGA. Additional genotyping and sequencing allowed the classification of a total of 628 genomic and 614 WGA amplicon samples. Heterozygous variants were identified by non-wild type melting pattern in 98% of genomic and 97% of WGA samples (P = 0.11). Wild types were correctly classified in 99% of genomic and 91% of WGA samples (P < 0.001).
In whole genome amplified DNA, high-resolution DNA melting curve analysis is a sensitive tool for SNP discovery through detection of heterozygote variants; however, it may misclassify a greater number of wild type samples.
Mutation Detection; Single Nucleotide Polymorphism; Genome Sequencing; Candidate Gene
With the recent development of microarray technologies, the comparability of gene expression data obtained from different platforms poses an important problem. We evaluated two widely used platforms, Affymetrix U133 Plus 2.0 and the Illumina HumanRef-8 v2 Expression Bead Chips, for comparability in a biological system in which changes may be subtle, namely fetal lung tissue as a function of gestational age.
We performed the comparison via sequence-based probe matching between the two platforms. "Significance grouping" was defined as a measure of comparability. Using both expression correlation and significance grouping as measures of comparability, we demonstrated that despite overall cross-platform differences at the single gene level, increased correlation between the two platforms was found in genes with higher expression level, higher probe overlap, and lower p-value. We also demonstrated that biological function as determined via KEGG pathways or GO categories is more consistent across platforms than single gene analysis.
We conclude that while the comparability of the platforms at the single gene level may be increased by increasing sample size, they are highly comparable ontologically even for subtle differences in a relatively small sample size. Biologically relevant inference should therefore be reproducible across laboratories using different platforms.
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
Rare loss-of-function folliculin (FLCN) mutations are the genetic cause of Birt-Hogg-Dubé syndrome, a monogenic disorder characterized by spontaneous pneumothorax, fibrofolliculomas, and kidney tumors. Loss-of-function folliculin mutations have also been described in pedigrees with familial spontaneous pneumothorax. Because the majority of patients with folliculin mutations have radiographic evidence of pulmonary cysts, folliculin has been hypothesized to contribute to the development of emphysema.
To determine whether folliculin sequence variants are risk factors for severe COPD, we genotyped seven previously reported Birt-Hogg-Dubé or familial spontaneous pneumothorax associated folliculin mutations in 152 severe COPD probands participating in the Boston Early-Onset COPD Study. We performed bidirectional resequencing of all 14 folliculin exons in a subset of 41 probands and subsequently genotyped four identified variants in an independent sample of345 COPD subjects from the National Emphysema Treatment Trial (cases) and 420 male smokers with normal lung function from the Normative Aging Study (controls).
None of the seven previously reported Birt-Hogg-Dubé or familial spontaneous pneumothorax mutations were observed in the 152 severe, early-onset COPD probands. Exon resequencing identified 31 variants, including two non-synonymous polymorphisms and two common non-coding polymorphisms. No significant association was observed for any of these four variants with presence of COPD or emphysema-related phenotypes.
Genetic variation in folliculin does not appear to be a major risk factor for severe COPD. These data suggest that familial spontaneous pneumothorax and COPD have distinct genetic causes, despite some overlap in radiographic characteristics.
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