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).
Background
The genetic risk factors for susceptibility to chronic obstructive
pulmonary disease (COPD) are still largely unknown. Additional genetic
variants are likely to be identified by genome-wide association studies in
larger cohorts or specific subgroups.
Methods
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.
Findings
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
controls.
Interpretation
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.
Funding
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
Affairs.
The genetic risk factors for chronic obstructive pulmonary disease (COPD) are still largely unknown. To date, genome-wide association studies (GWASs) of limited size have identified several novel risk loci for COPD at CHRNA3/CHRNA5/IREB2, HHIP and FAM13A; additional loci may be identified through larger studies. We performed a GWAS using a total of 3499 cases and 1922 control subjects from four cohorts: the Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE); the Normative Aging Study (NAS) and National Emphysema Treatment Trial (NETT); Bergen, Norway (GenKOLS); and the COPDGene study. Genotyping was performed on Illumina platforms with additional markers imputed using 1000 Genomes data; results were summarized using fixed-effect meta-analysis. We identified a new genome-wide significant locus on chromosome 19q13 (rs7937, OR = 0.74, P = 2.9 × 10(-9)). Genotyping this single nucleotide polymorphism (SNP) and another nearby SNP in linkage disequilibrium (rs2604894) in 2859 subjects from the family-based International COPD Genetics Network study (ICGN) demonstrated supportive evidence for association for COPD (P = 0.28 and 0.11 for rs7937 and rs2604894), pre-bronchodilator FEV(1) (P = 0.08 and 0.04) and severe (GOLD 3&4) COPD (P = 0.09 and 0.017). This region includes RAB4B, EGLN2, MIA and CYP2A6, and has previously been identified in association with cigarette smoking behavior.
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.
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.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.