Severe α1‐antitrypsin (AAT) deficiency is an autosomal recessive genetic condition associated with an increased but variable risk for chronic obstructive pulmonary disease (COPD). A study was undertaken to assess the impact of chronic bronchitis, pneumonia, asthma and sex on the development of COPD in individuals with severe AAT deficiency.
The AAT Genetic Modifier Study is a multicentre family‐based cohort study designed to study the genetic and epidemiological determinants of COPD in AAT deficiency. 378 individuals (age range 33–80 years), confirmed to be homozygous for the SERPINA1 Z mutation, were included in the analyses. The primary outcomes of interest were a quantitative outcome, forced expiratory volume in 1 s (FEV1) percentage predicted, and a qualitative outcome, severe airflow obstruction (FEV1 <50% predicted).
In multivariate analysis of the overall cohort, cigarette smoking, sex, asthma, chronic bronchitis and pneumonia were risk factors for reduced FEV1 percentage predicted and severe airflow obstruction (p<0.01). Index cases had lower FEV1 values, higher smoking histories and more reports of adult asthma, pneumonia and asthma before age 16 than non‐index cases (p<0.01). Men had lower pre‐ and post‐bronchodilator FEV1 percentage predicted than women (p<0.0001); the lowest FEV1 values were observed in men reporting a history of childhood asthma (26.9%). This trend for more severe obstruction in men remained when index and non‐index groups were examined separately, with men representing the majority of non‐index individuals with airflow obstruction (71%). Chronic bronchitis (OR 3.8, CI 1.8 to 12.0) and a physician's report of asthma (OR 4.2, CI 1.4 to 13.1) were predictors of severe airflow obstruction in multivariate analysis of non‐index men but not women.
In individuals with severe AAT deficiency, sex, asthma, chronic bronchitis and pneumonia are risk factors for severe COPD, in addition to cigarette smoking. These results suggest that, in subjects severely deficient in AAT, men, individuals with symptoms of chronic bronchitis and/or a past diagnosis of asthma or pneumonia may benefit from closer monitoring and potentially earlier treatment.
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.
68kDa (TGFBI); C-Reactive Protein (CRP); Chronic Obstructive Pulmonary Disease (COPD); Illumina GoldenGate Methylation Cancer Panel I; alpha-1 antitrypsin (AAT) deficiency; beta-induced; methylation; smoking behaviors; transforming growth factor
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; FEV1/FVC, 1.5 × 10−35; FEV1, 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.
chronic obstructive pulmonary disease; epigenetics; DNA methylation; smoking
Rationale: The characterization of young adults who develop late-onset diseases may augment the detection of novel genes and promote new pathogenic insights.
Methods: We analyzed data from 2,500 individuals of African and European ancestry in the COPDGene Study. Subjects with severe, early-onset chronic obstructive pulmonary disease (COPD) (n = 70, age < 55 yr, FEV1 < 50% predicted) were compared with older subjects with COPD (n = 306, age > 64 yr, FEV1 < 50% predicted).
Measurements and Main Results: Subjects with severe, early-onset COPD were predominantly females (66%), P = 0.0004. Proportionally, early-onset COPD was seen in 42% (25 of 59) of African Americans versus 14% (45 of 317) of non-Hispanic whites, P < 0.0001. Other risk factors included current smoking (56 vs. 17%, P < 0.0001) and self-report of asthma (39 vs. 25%, P = 0.008). Maternal smoking (70 vs. 44%, P = 0.0001) and maternal COPD (23 vs. 12%, P = 0.03) were reported more commonly in subjects with early-onset COPD. Multivariable regression analysis found association with African American race, odds ratio (OR), 7.5 (95% confidence interval [CI], 2.3–24; P = 0.0007); maternal COPD, OR, 4.7 (95% CI, 1.3–17; P = 0.02); female sex, OR, 3.1 (95% CI, 1.1–8.7; P = 0.03); and each pack-year of smoking, OR, 0.98 (95% CI, 0.96–1.0; P = 0.03).
Conclusions: These observations support the hypothesis that severe, early-onset COPD is prevalent in females and is influenced by maternal factors. Future genetic studies should evaluate (1) gene-by-sex interactions to address sex-specific genetic contributions and (2) gene-by-race interactions.
chronic obstructive pulmonary disease; female; African Americans
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
COPD; Genetics; Association analysis; Consortium
Severe α1-antitrypsin (AAT) deficiency is a proven genetic risk factor for chronic obstructive pulmonary disease (COPD), especially in individuals who smoke. There is marked variability in the development of lung disease in individuals homozygous (PI ZZ) for this autosomal recessive condition, suggesting that modifier genes could be important. We hypothesized that genetic determinants of obstructive lung disease may be modifiers of airflow obstruction in individuals with severe AAT deficiency. To identify modifier genes, we performed family-based association analyses for 10 genes previously associated with asthma and/or COPD, including IL10, TNF, GSTP1, NOS1, NOS3, SERPINA3, SERPINE2, SFTPB, TGFB1, and EPHX1. All analyses were performed in a cohort of 378 PI ZZ individuals from 167 families. Quantitative spirometric phenotypes included forced expiratory volume in one second (FEV1) and the ratio of FEV1/forced vital capacity (FVC). A qualitative phenotype of moderate-to-severe COPD was defined for individuals with FEV1 ⩽ 50 percent predicted. Six of 11 single-nucleotide polymorphisms (SNPs) in IL10 (P = 0.0005–0.05) and 3 of 5 SNPs in TNF (P = 0.01–0.05) were associated with FEV1 and/or FEV1/FVC. IL10 SNPs also demonstrated association with the qualitative COPD phenotype. When phenotypes of individuals with a physician's diagnosis of asthma were excluded, IL10 SNPs remained significantly associated, suggesting that the association with airflow obstruction was independent of an association with asthma. Haplotype analysis of IL10 SNPs suggested the strongest association with IL10 promoter SNPs. IL10 is likely an important modifier gene for the development of COPD in individuals with severe AAT deficiency.
chronic obstructive pulmonary disease; genetic modifiers; interleukin 10; family-based association analysis
Rationale: Computed tomography (CT) scanning of the lung may reduce phenotypic heterogeneity in defining subjects with chronic obstructive pulmonary disease (COPD), and allow identification of genetic determinants of emphysema severity and distribution.
Objectives: We sought to identify genes associated with CT scan distribution of emphysema in individuals without α1-antitrypsin deficiency but with severe COPD.
Methods: We evaluated baseline CT densitometry phenotypes in 282 individuals with emphysema enrolled in the Genetics Ancillary Study of the National Emphysema Treatment Trial, and used regression models to identify genetic variants associated with emphysema distribution.
Measurements and Main Results: Emphysema distribution was assessed by two methods—assessment by radiologists and by computerized density mask quantitation, using a threshold of −950 Hounsfield units. A total of 77 polymorphisms in 20 candidate genes were analyzed for association with distribution of emphysema. GSTP1, EPHX1, and MMP1 polymorphisms were associated with the densitometric, apical-predominant distribution of emphysema (p value range = 0.001–0.050). When an apical-predominant phenotype was defined by the radiologist scoring method, GSTP1 and EPHX1 single-nucleotide polymorphisms were found to be significantly associated. In a case–control analysis of COPD susceptibility limited to cases with densitometric upper-lobe–predominant cases, the EPHX1 His139Arg single-nucleotide polymorphism was associated with COPD (p = 0.005).
Conclusions: Apical and basal emphysematous destruction appears to be influenced by different genes. Polymorphisms in the xenobiotic enzymes, GSTP1 and EPHX1, are associated with apical-predominant emphysema. Altered detoxification of cigarette smoke metabolites may contribute to emphysema distribution, and these findings may lead to further insight into genetic determinants of emphysema.
COPD; genetics; association analysis; computed tomography; emphysema
Whether African Americans (AA) are more susceptible to COPD than non-Hispanic Whites (NHW) and whether racial differences in disease phenotype exist is controversial. The objective is to determine racial differences in the extent of emphysema and airway remodeling in COPD.
First, 2,500 subjects enrolled in the COPDGene study were used to evaluate racial differences in quantitative CT (QCT) parameters of % emphysema, air trapping and airway wall thickness. Independent variables studied included race, age, gender, education, BMI, pack-years, smoking status, age at smoking initiation, asthma, previous work in dusty job, CT scanner and center of recruitment.
Of the 1,063 subjects with GOLD Stage II-IV COPD, 200 self-reported as AA. AAs had a lower mean % emphysema (13.1 % vs. 16.1%, p = 0.005) than NHW and proportionately less emphysema in the lower lung zones. After adjustment for covariates, there was no statistical difference by race in air trapping or airway wall thickness. Measured QCT parameters were more predictive of poor functional status in NHWs compared to AAs.
AAs have less emphysema than NHWs but the same degree of airway disease. Additional factors not easily assessed by current QCT techniques may account for the poor functional status in AAs.
Airway wall thickness; Air trapping; Chronic obstructive pulmonary disease; Emphysema; Quantitative CT; Race
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
α1-Antitrypsin is primarily synthesised in the liver, circulates to the lung and protects pulmonary tissues from proteolytic damage. The Z mutant (Glu342Lys) undergoes inactivating conformational change and polymerises. Polymers are retained within the hepatocyte endoplasmic reticulum (ER) in homozygous (PiZZ) individuals, predisposing the individuals to hepatic cirrhosis and emphysema. Latency is an analogous process of inactivating, intra-molecular conformational change and may co-occur with polymerisation. However, the relationship between latency and polymerisation remained unexplored in the absence of a suitable probe. We have developed a novel monoclonal antibody specific for latent α1-antitrypsin and used it in combination with a polymer-specific antibody, to assess the association of both conformers in vitro, in disease and during augmentation therapy. In vitro kinetics analysis showed polymerisation dominated the pathway but latency could be promoted by stabilising monomeric α1-antitrypsin. Polymers were extensively produced in hepatocytes and a cell line expressing Z α1-antitrypsin but the latent protein was not detected despite manipulation of the secretory pathway. However, α1-antitrypsin augmentation therapy contains latent α1-antitrypsin, as did the plasma of 63/274 PiZZ individuals treated with augmentation therapy but 0/264 who were not receiving this medication (p < 10−14). We conclude that latent α1-antitrypsin is a by-product of the polymerisation pathway, that the intracellular folding environment is resistant to formation of the latent conformer but that augmentation therapy introduces latent α1-antitrypsin into the circulation. A suite of monoclonal antibodies and methodologies developed in this study can characterise α1-antitrypsin folding and conformational transitions, and screen methods to improve augmentation therapy.
α1-Antitrypsin; Latency; Polymerisation; Monoclonal antibodies; Augmentation therapy
Comorbidities are common in COPD, but quantifying their burden is difficult. Currently there is a COPD-specific comorbidity index to predict mortality and another to predict general quality of life. We sought to develop and validate a COPD-specific comorbidity score that reflects comorbidity burden on patient-centered outcomes.
Materials and Methods
Using the COPDGene study (GOLD II-IV COPD), we developed comorbidity scores to describe patient-centered outcomes employing three techniques: 1) simple count, 2) weighted score, and 3) weighted score based upon statistical selection procedure. We tested associations, area under the Curve (AUC) and calibration statistics to validate scores internally with outcomes of respiratory disease-specific quality of life (St. George's Respiratory Questionnaire, SGRQ), six minute walk distance (6MWD), modified Medical Research Council (mMRC) dyspnea score and exacerbation risk, ultimately choosing one score for external validation in SPIROMICS.
Associations between comorbidities and all outcomes were comparable across the three scores. All scores added predictive ability to models including age, gender, race, current smoking status, pack-years smoked and FEV1 (p<0.001 for all comparisons). Area under the curve (AUC) was similar between all three scores across outcomes: SGRQ (range 0·7624–0·7676), MMRC (0·7590–0·7644), 6MWD (0·7531–0·7560) and exacerbation risk (0·6831–0·6919). Because of similar performance, the comorbidity count was used for external validation. In the SPIROMICS cohort, the comorbidity count performed well to predict SGRQ (AUC 0·7891), MMRC (AUC 0·7611), 6MWD (AUC 0·7086), and exacerbation risk (AUC 0·7341).
Quantifying comorbidity provides a more thorough understanding of the risk for patient-centered outcomes in COPD. A comorbidity count performs well to quantify comorbidity in a diverse population with COPD.
There is growing evidence that many diseases develop, progress, and respond to therapy differently in men and women. This variability may manifest as a result of sex-specific structures in gene regulatory networks that influence how those networks operate. However, there are few methods to identify and characterize differences in network structure, slowing progress in understanding mechanisms driving sexual dimorphism.
Here we apply an integrative network inference method, PANDA (Passing Attributes between Networks for Data Assimilation), to model sex-specific networks in blood and sputum samples from subjects with Chronic Obstructive Pulmonary Disease (COPD). We used a jack-knifing approach to build an ensemble of likely networks for each sex. By adapting statistical methods to compare these network ensembles, we were able to identify strong differential-targeting patterns associated with functionally-related sets of genes, including those involved in mitochondrial function and energy metabolism. Network analysis also identified several potential sex- and disease-specific transcriptional regulators of these pathways.
Network analysis yielded insight into potential mechanisms driving sexual dimorphism in COPD that were not evident from gene expression analysis alone. We believe our ensemble approach to network analysis provides a principled way to capture sex-specific regulatory relationships and could be applied to identify differences in gene regulatory patterns in a wide variety of diseases and contexts.
Electronic supplementary material
The online version of this article (doi:10.1186/s12918-014-0118-y) contains supplementary material, which is available to authorized users.
Network modeling; Gene regulation; Regulatory networks; Sexual-dimorphism; Chronic Obstructive Lung Disease
Chronic obstructive pulmonary disease (COPD) has been classically divided into blue bloaters and pink puffers. The utility of these clinical subtypes is unclear. However, the broader distinction between airway-predominant and emphysema-predominant COPD may be clinically relevant. The objective was to define clinical features of emphysema-predominant and non-emphysematous COPD patients.
Current and former smokers from the Genetic Epidemiology of COPD Study (COPDGene) had chest computed tomography (CT) scans with quantitative image analysis. Emphysema-predominant COPD was defined by low attenuation area at -950 Hounsfield Units (LAA-950) ≥10%. Non-emphysematous COPD was defined by airflow obstruction with minimal to no emphysema (LAA-950 < 5%).
Out of 4197 COPD subjects, 1687 were classified as emphysema-predominant and 1817 as non-emphysematous; 693 had LAA-950 between 5–10% and were not categorized. Subjects with emphysema-predominant COPD were older (65.6 vs 60.6 years, p < 0.0001) with more severe COPD based on airflow obstruction (FEV1 44.5 vs 68.4%, p < 0.0001), greater exercise limitation (6-minute walk distance 1138 vs 1331 ft, p < 0.0001) and reduced quality of life (St. George’s Respiratory Questionnaire score 43 vs 31, p < 0.0001). Self-reported diabetes was more frequent in non-emphysematous COPD (OR 2.13, p < 0.001), which was also confirmed using a strict definition of diabetes based on medication use. The association between diabetes and non-emphysematous COPD was replicated in the ECLIPSE study.
Non-emphysematous COPD, defined by airflow obstruction with a paucity of emphysema on chest CT scan, is associated with an increased risk of diabetes. COPD patients without emphysema may warrant closer monitoring for diabetes, hypertension, and hyperlipidemia and vice versa.
Clinicaltrials.gov identifiers: COPDGene NCT00608764, ECLIPSE NCT00292552.
Electronic supplementary material
The online version of this article (doi:10.1186/1471-2466-14-164) contains supplementary material, which is available to authorized users.
Airway disease; CT scan; Diabetes mellitus; Emphysema; Spirometry
Preserved Ratio Impaired Spirometry (PRISm), defined as a reduced FEV1 in the setting of a preserved FEV1/FVC ratio, is highly prevalent and is associated with increased respiratory symptoms, systemic inflammation, and mortality. Studies investigating quantitative chest tomographic features, genetic associations, and subtypes in PRISm subjects have not been reported.
Data from current and former smokers enrolled in COPDGene (n = 10,192), an observational, cross-sectional study which recruited subjects aged 45–80 with ≥10 pack years of smoking, were analyzed. To identify epidemiological and radiographic predictors of PRISm, we performed univariate and multivariate analyses comparing PRISm subjects both to control subjects with normal spirometry and to subjects with COPD. To investigate common genetic predictors of PRISm, we performed a genome-wide association study (GWAS). To explore potential subgroups within PRISm, we performed unsupervised k-means clustering.
The prevalence of PRISm in COPDGene is 12.3%. Increased dyspnea, reduced 6-minute walk distance, increased percent emphysema and decreased total lung capacity, as well as increased segmental bronchial wall area percentage were significant predictors (p-value <0.05) of PRISm status when compared to control subjects in multivariate models. Although no common genetic variants were identified on GWAS testing, a significant association with Klinefelter’s syndrome (47XXY) was observed (p-value < 0.001). Subgroups identified through k-means clustering include a putative “COPD-subtype”, “Restrictive-subtype”, and a highly symptomatic “Metabolic-subtype”.
PRISm subjects are clinically and genetically heterogeneous. Future investigations into the pathophysiological mechanisms behind and potential treatment options for subgroups within PRISm are warranted.
Clinicaltrials.gov Identifier: NCT000608764.
Electronic supplementary material
The online version of this article (doi:10.1186/s12931-014-0089-y) contains supplementary material, which is available to authorized users.
Spirometry; Restriction; Lung diseases; Smoking
Rationale: Systemic glucocorticoids are used therapeutically to treat a variety of medical conditions. Epigenetic processes such as DNA methylation may reflect exposure to glucocorticoids and may be involved in mediating the responses and side effects associated with these medications.
Objectives: To test the hypothesis that differences in DNA methylation are associated with current systemic steroid use.
Methods: We obtained DNA methylation data at 27,578 CpG sites in 14,475 genes throughout the genome in two large, independent cohorts: the International COPD Genetics Network (ndiscovery = 1,085) and the Boston Early Onset COPD study (nreplication = 369). Sites were tested for association with current systemic steroid use using generalized linear mixed models.
Measurements and Main Results: A total of 511 sites demonstrated significant differential methylation by systemic corticosteroid use in all three of our primary models. Pyrosequencing validation confirmed robust differential methylation at CpG sites annotated to genes such as SLC22A18, LRP3, HIPK3, SCNN1A, FXYD1, IRF7, AZU1, SIT1, GPR97, ABHD16B, and RABGEF1. Functional annotation clustering demonstrated significant enrichment in intrinsic membrane components, hemostasis and coagulation, cellular ion homeostasis, leukocyte and lymphocyte activation and chemotaxis, protein transport, and responses to nutrients.
Conclusions: Our analyses suggest that systemic steroid use is associated with site-specific differential methylation throughout the genome. Differentially methylated CpG sites were found in biologically plausible and previously unsuspected pathways; these genes and pathways may be relevant in the development of novel targeted therapies.
DNA methylation; glucocorticoids; chronic obstructive pulmonary disease
We propose a method for testing gene–environment (G × E) interactions on a complex trait in family-based studies in which a phenotypic ascertainment criterion has been imposed. This novel approach employs G-estimation, a semiparametric estimation technique from the causal inference literature, to avoid modeling of the association between the environmental exposure and the phenotype, to gain robustness against unmeasured confounding due to population substructure, and to acknowledge the ascertainment conditions. The proposed test allows for incomplete parental genotypes. It is compared by simulation studies to an analogous conditional likelihood–based approach and to the QBAT-I test, which also invokes the G-estimation principle but ignores ascertainment. We apply our approach to a study of chronic obstructive pulmonary disorder.
Causal inference; COPD; Family-based association; G-estimation; Gene–environment interaction
The impact of cigarette smoking can persist for extended periods following smoking cessation and may involve epigenetic reprogramming. Changes in DNA methylation associated with smoking may help to identify molecular pathways that contribute to the latency between exposure and disease onset. Cross-sectional cohort data from subjects in the International COPD Genetics Network (n = 1085) and the Boston Early-Onset COPD study (n = 369) were analyzed as the discovery and replication cohorts, respectively. Genome-wide methylation data on 27 578 CpG sites in 14 475 genes were obtained on DNA from peripheral blood leukocytes using the Illumina HumanMethylation27K Beadchip in both cohorts. We identified 15 sites significantly associated with current smoking, 2 sites associated with cumulative smoke exposure, and, within the subset of former smokers, 3 sites associated with time since quitting cigarettes. Two loci, factor II receptor-like 3 (F2RL3) and G-protein-coupled receptor 15 (GPR15), were significantly associated in all three analyses and were validated by pyrosequencing. These findings (i) identify a novel locus (GPR15) associated with cigarette smoking and (ii) suggest the existence of dynamic, site-specific methylation changes in response to smoking which may contribute to the extended risks associated with cigarette smoking that persist after cessation.
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.
It is useful to have robust gene-environment interaction tests that can utilize a variety of family structures in an efficient way. This paper focuses on tests for gene-environment interaction in the presence of main genetic and environmental effects. The objective is to develop powerful tests that can combine trio data with parental genotypes and discordant sibships when parents genotypes are missing. We first make a modest improvement on a method for discordant sibs (discordant on phenotype), but the approach does not allow one to use families when all offspring are affected, e.g. trios. We then make a modest improvement on a Mendelian transmission-based approach that is inefficient when discordant sibs are available, but can be applied to any nuclear family. Finally, we propose a hybrid approach that utilizes the most efficient method for a specific family type, then combines over families. We utilize this hybrid approach to analyze a chronic obstructive pulmonary disorder dataset to test for gene-environment interaction in the Serpine2 gene with smoking. The methods are freely available in the R package fbati.
Gene-Environment Interaction; Family-Based Association Tests; Candidate Gene Analysis; Binary Trait; COPD; Serpine2
To investigate the association between methylation of transposable elements Alu and long-interspersed nuclear elements (LINE-1) and lung function.
Outpatient Veterans Administration facilities in greater Boston, Massachusetts, USA.
Individuals from the Veterans Administration Normative Aging Study, a longitudinal study of aging in men, evaluated between 1999 and 2007. The majority (97%) were white.
Primary and secondary outcome measures
Primary predictor was methylation, assessed using PCR-pyrosequencing after bisulphite treatment. Primary outcome was lung function as assessed by spirometry, performed according to American Thoracic Society/European Respiratory Society guidelines at the same visit as the blood draws.
In multivariable models adjusted for age, height, body mass index (BMI), pack-years of smoking, current smoking and race, Alu hypomethylation was associated with lower forced expiratory volume in 1 s (FEV1) (β=28 ml per 1% change in Alu methylation, p=0.017) and showed a trend towards association with a lower forced vital capacity (FVC) (β=27 ml, p=0.06) and lower FEV1/FVC (β=0.3%, p=0.058). In multivariable models adjusted for age, height, BMI, pack-years of smoking, current smoking, per cent lymphocytes, race and baseline lung function, LINE-1 hypomethylation was associated with more rapid decline of FEV1 (β=6.9 ml/year per 1% change in LINE-1 methylation, p=0.005) and of FVC (β=9.6 ml/year, p=0.002).
In multiple regression analysis, Alu hypomethylation was associated with lower lung function, and LINE-1 hypomethylation was associated with more rapid lung function decline in a cohort of older and primarily white men from North America. Future studies should aim to replicate these findings and determine if Alu or LINE-1 hypomethylation may be due to specific and modifiable environmental exposures.
Acute exacerbations are a significant source of morbidity and mortality associated with chronic obstructive pulmonary disease. Among patients with COPD, some patients suffer an inordinate number of exacerbations while others remain relatively protected. We undertook a study to determine the clinical factors associated with "frequent exacerbator" status within a population of subjects with severe COPD.
Case-control cohort recruited from two Boston-area practices. All subjects had GOLD stage 3 or 4 (FEV1 ≤50% predicted) COPD. "Frequent exacerbators" (n=192) had an average of ≥2 moderate-to-severe exacerbations per year while "non-exacerbators" (n=153) had no exacerbations in the preceding 12 months. Multivariate logistic regression was performed to determine the significant clinical predictors of "frequent exacerbator" status.
Physician-diagnosed asthma was a significant predictor of frequent exacerbations. Within a subset of our cohort, the modified Medical Research Council dyspnea score and FEF25–75 % predicted were also significant clinical predictors of frequent exacerbator status (p<0.05). Differences in exacerbation frequency were not found to be due to increased current tobacco use or decreased rates of maintenance medication use.
Within our severe COPD cohort, a history of physician-diagnosed asthma was found to be a significant clinical predictor of frequent exacerbations. Although traditional risk factors such as decreased FEV1% predicted were not significantly associated with frequent exacerbator status, lower mid-expiratory flow rates, as assessed by FEF 25–75 % predicted, were significantly associated with frequent exacerbations in a subset of our cohort.
Chronic obstructive pulmonary disease (COPD) is characterized by alveolar destruction and abnormal inflammatory responses to noxious stimuli. Surfactant protein–D (SFTPD) is immunomodulatory and essential to host defense. We hypothesized that polymorphisms in SFTPD could influence the susceptibility to COPD. We genotyped six single-nucleotide polymorphisms (SNPs) in surfactant protein D in 389 patients with COPD in the National Emphysema Treatment Trial (NETT) and 472 smoking control subjects from the Normative Aging Study (NAS). Case-control association analysis was performed using Cochran–Armitage trend tests and multivariate logistic regression. The replication of significant associations was attempted in the Boston Early-Onset COPD Study, the Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE) Study, and the Bergen Cohort. We also correlated SFTPD genotypes with serum concentrations of surfactant protein–D (SP-D) in the ECLIPSE Study. In the NETT–NAS case-control analysis, four SFTPD SNPs were associated with susceptibility to COPD: rs2245121 (P = 0.01), rs911887 (P = 0.006), rs6413520 (P = 0.004), and rs721917 (P = 0.006). In the family-based analysis of the Boston Early-Onset COPD Study, rs911887 was associated with prebronchodilator and postbronchodilator FEV1 (P = 0.003 and P = 0.02, respectively). An intronic SNP in SFTPD, rs7078012, was associated with COPD in the ECLIPSE Study and the Bergen Cohort. Multiple SFTPD SNPs were associated with serum SP-D concentrations in the ECLIPSE Study. We demonstrated an association of polymorphisms in SFTPD with COPD in multiple populations. We demonstrated a correlation between SFTPD SNPs and SP-D protein concentrations. The SNPs associated with COPD and SP-D concentrations differed, suggesting distinct genetic influences on susceptibility to COPD and SP-D concentrations.
COPD; surfactant protein–D; single-nucleotide polymorphisms; genetics
The development of COPD in subjects with alpha-1 antitrypsin (AAT) deficiency is likely to be influenced by modifier genes. Genome-wide association studies and integrative genomics approaches in COPD have demonstrated significant associations with SNPs in the chromosome 15q region that includes CHRNA3 (cholinergic nicotine receptor alpha3) and IREB2 (iron regulatory binding protein 2).
We investigated whether SNPs in the chromosome 15q region would be modifiers for lung function and COPD in AAT deficiency.
The current analysis included 378 PIZZ subjects in the AAT Genetic Modifiers Study and a replication cohort of 458 subjects from the UK AAT Deficiency National Registry. Nine SNPs in LOC123688, CHRNA3 and IREB2 were selected for genotyping. FEV1 percent of predicted and FEV1/FVC ratio were analyzed as quantitative phenotypes. Family-based association analysis was performed in the AAT Genetic Modifiers Study. In the replication set, general linear models were used for quantitative phenotypes and logistic regression models were used for the presence/absence of emphysema or COPD.
Three SNPs (rs2568494 in IREB2, rs8034191 in LOC123688, and rs1051730 in CHRNA3) were associated with pre-bronchodilator FEV1 percent of predicted in the AAT Genetic Modifiers Study. Two SNPs (rs2568494 and rs1051730) were associated with the post-bronchodilator FEV1 percent of predicted and pre-bronchodilator FEV1/FVC ratio; SNP-by-gender interactions were observed. In the UK National Registry dataset, rs2568494 was significantly associated with emphysema in the male subgroup; significant SNP-by-smoking interactions were observed.
IREB2 and CHRNA3 are potential genetic modifiers of COPD phenotypes in individuals with severe AAT deficiency and may be sex-specific in their impact.
CHRNA3; Chronic obstructive pulmonary disease; Genetic association analysis; Genetic modifiers; IREB2