Refining susceptibility loci of chronic obstructive pulmonary disease with lung eqtls

Chronic obstructive pulmonary disease (COPD) is the fourth leading cause of mortality worldwide. Recent genome-wide association studies (GWAS) have identified robust susceptibility loci associated with COPD. However, the mechanisms mediating the risk conferred by these loci remain to be found. The goal of this study was to identify causal genes/variants within susceptibility loci associated with COPD. In the discovery cohort, genome-wide gene expression profiles of 500 non-tumor lung specimens were obtained from patients undergoing lung surgery. Blood-DNA from the same patients were genotyped for 1,2 million SNPs. Following genotyping and gene expression quality control filters, 409 samples were analyzed. Lung expression quantitative trait loci (eQTLs) were identified and overlaid onto three COPD susceptibility loci derived from GWAS; 4q31 (HHIP), 4q22 (FAM13A), and 19q13 (RAB4B, EGLN2, MIA, CYP2A6). Significant eQTLs were replicated in two independent datasets (n = 363 and 339). SNPs previously associated with COPD and lung function on 4q31 (rs1828591, rs13118928) were associated with the mRNA expression of HHIP. An association between mRNA expression level of FAM13A and SNP rs2045517 was detected at 4q22, but did not reach statistical significance. At 19q13, significant eQTLs were detected with EGLN2. In summary, this study supports HHIP, FAM13A, and EGLN2 as the most likely causal COPD genes on 4q31, 4q22, and 19q13, respectively. Strong lung eQTL SNPs identified in this study will need to be tested for association with COPD in case-control studies. Further functional studies will also be needed to understand the role of genes regulated by disease-related variants in COPD

Impact of Statins on Gene Expression in Human Lung Tissues

Statins are 3-hydroxy-3-methylglutaryl-coenzyme A reductase inhibitors that alter the synthesis of cholesterol. Some studies have shown a significant association of statins with improved respiratory health outcomes of patients with asthma, chronic obstructive pulmonary disease and lung cancer. Here we hypothesize that statins impact gene expression in human lungs and may reveal the pleiotropic effects of statins that are taking place directly in lung tissues. Human lung tissues were obtained from patients who underwent lung resection or transplantation. Gene expression was measured on a custom Affymetrix array in a discovery cohort (n = 408) and two replication sets (n = 341 and 282). Gene expression was evaluated by linear regression between statin users and non-users, adjusting for age, gender, smoking status, and other covariables. The results of each cohort were combined in a meta-analysis and biological pathways were studied using Gene Set Enrichment Analysis. The discovery set included 141 statin users. The lung mRNA expression levels of eighteen and three genes were up-regulated and down-regulated in statin users (FDR < 0.05), respectively. Twelve of the up-regulated genes were replicated in the first replication set, but none in the second (p-value < 0.05). Combining the discovery and replication sets into a meta-analysis improved the significance of the 12 up-regulated genes, which includes genes encoding enzymes and membrane proteins involved in cholesterol biosynthesis. Canonical biological pathways altered by statins in the lung include cholesterol, steroid, and terpenoid backbone biosynthesis. No genes encoding inflammatory, proteases, pro-fibrotic or growth factors were altered by statins, suggesting that the direct effect of statin in the lung do not go beyond its antilipidemic action. Although more studies are needed with specific lung cell types and different classes and doses of statins, the improved health outcomes and survival observed in statin users with chronic lung diseases do not seem to be mediated through direct regulation of gene expression in the lung

A Genome-Wide Linkage Study for Chronic Obstructive Pulmonary Disease in a Dutch Genetic Isolate Identifies Novel Rare Candidate Variants

Chronic obstructive pulmonary disease (COPD) is a complex and heritable disease, associated with multiple genetic variants. Specific familial types of COPD may be explained by rare variants, which have not been widely studied. We aimed to discover rare genetic variants underlying COPD through a genome-wide linkage scan. Affected-only analysis was performed using the 6K Illumina Linkage IV Panel in 142 cases clustered in 27 families from a genetic isolate, the Erasmus Rucphen Family (ERF) study. Potential causal variants were identified by searching for shared rare variants in the exome-sequence data of the affected members of the families contributing most to the linkage peak. The identified rare variants were then tested for association with COPD in a large meta-analysis of several cohorts. Significant evidence for linkage was observed on chromosomes 15q14-15q25 [logarithm of the odds (LOD) score = 5.52], 11p15.4-11q14.1 (LOD = 3.71) and 5q14.3-5q33.2 (LOD = 3.49). In the chromosome 15 peak, that harbors the known COPD locus for nicotinic receptors, and in the chromosome 5 peak we could not identify shared variants. In the chromosome 11 locus, we identified four rare (minor allele frequency (MAF) <0.02), predicted pathogenic, missense variants. These were shared among the affected family members. The identified variants localize to genes including neuroblast differentiation-associated protein (AHNAK), previously associated with blood biomarkers in COPD, phospholipase C Beta 3 (PLCB3), shown to increase airway hyper-responsiveness, solute carrier family 22-A11 (SLC22A11), involved in amino acid metabolism and ion transport, and metallothionein-like protein 5 (MTL5), involved in nicotinate and nicotinamide metabolism. Association of SLC22A11 and MTL5 variants were confirmed in the meta-analysis of 9,888 cases and 27,060 controls. In conclusion, we have identified novel rare variants in plausible genes related to COPD. Further studies utilizing large sample whole-genome sequencing should further confirm the associations at chromosome 11 and investigate the chromosome 15 and 5 linked regions

ITGB5 and AGFG1 variants are associated with severity of airway responsiveness

Background: Airway hyperresponsiveness (AHR), a primary characteristic of asthma, involves increased airway smooth muscle contractility in response to certain exposures. We sought to determine whether common genetic variants were associated with AHR severity. Methods: A genome-wide association study (GWAS) of AHR, quantified as the natural log of the dosage of methacholine causing a 20% drop in FEV1, was performed with 994 non-Hispanic white asthmatic subjects from three drug clinical trials: CAMP, CARE, and ACRN. Genotyping was performed on Affymetrix 6.0 arrays, and imputed data based on HapMap Phase 2, was used to measure the association of SNPs with AHR using a linear regression model. Replication of primary findings was attempted in 650 white subjects from DAG, and 3,354 white subjects from LHS. Evidence that the top SNPs were eQTL of their respective genes was sought using expression data available for 419 white CAMP subjects. Results: The top primary GWAS associations were in rs848788 (P-value 7.2E-07) and rs6731443 (P-value 2.5E-06), located within the ITGB5 and AGFG1 genes, respectively. The AGFG1 result replicated at a nominally significant level in one independent population (LHS P-value 0.012), and the SNP had a nominally significant unadjusted P-value (0.0067) for being an eQTL of AGFG1. Conclusions: Based on current knowledge of ITGB5 and AGFG1, our results suggest that variants within these genes may be involved in modulating AHR. Future functional studies are required to confirm that our associations represent true biologically significant findings

Lung eQTLs to Help Reveal the Molecular Underpinnings of Asthma

Genome-wide association studies (GWAS) have identified loci reproducibly associated with pulmonary diseases; however, the molecular mechanism underlying these associations are largely unknown. The objectives of this study were to discover genetic variants affecting gene expression in human lung tissue, to refine susceptibility loci for asthma identified in GWAS studies, and to use the genetics of gene expression and network analyses to find key molecular drivers of asthma. We performed a genome-wide search for expression quantitative trait loci (eQTL) in 1,111 human lung samples. The lung eQTL dataset was then used to inform asthma genetic studies reported in the literature. The top ranked lung eQTLs were integrated with the GWAS on asthma reported by the GABRIEL consortium to generate a Bayesian gene expression network for discovery of novel molecular pathways underpinning asthma. We detected 17,178 cis- and 593 trans- lung eQTLs, which can be used to explore the functional consequences of loci associated with lung diseases and traits. Some strong eQTLs are also asthma susceptibility loci. For example, rs3859192 on chr17q21 is robustly associated with the mRNA levels of GSDMA (P = 3.55 × 10(-151)). The genetic-gene expression network identified the SOCS3 pathway as one of the key drivers of asthma. The eQTLs and gene networks identified in this study are powerful tools for elucidating the causal mechanisms underlying pulmonary disease. This data resource offers much-needed support to pinpoint the causal genes and characterize the molecular function of gene variants associated with lung diseases

Multidrug resistance-associated protein-1 (MRP1) genetic variants, MRP1 protein levels and severity of COPD

<p>Abstract</p> <p>Background</p> <p>Multidrug resistance-associated protein-1 (MRP1) protects against oxidative stress and toxic compounds generated by cigarette smoking, which is the main risk factor for chronic obstructive pulmonary disease (COPD). We have previously shown that single nucleotide polymorphisms (SNPs) in <it>MRP1 </it>significantly associate with level of FEV₁in two independent population based cohorts. The aim of our study was to assess the associations of <it>MRP1 </it>SNPs with FEV₁level, MRP1 protein levels and inflammatory markers in bronchial biopsies and sputum of COPD patients.</p> <p>Methods</p> <p>Five SNPs (rs212093, rs4148382, rs504348, rs4781699, rs35621) in <it>MRP1 </it>were genotyped in 110 COPD patients. The effects of <it>MRP1 </it>SNPs were analyzed using linear regression models.</p> <p>Results</p> <p>One SNP, rs212093 was significantly associated with a higher FEV₁level and less airway wall inflammation. Another SNP, rs4148382 was significantly associated with a lower FEV₁level, higher number of inflammatory cells in induced sputum and with a higher MRP1 protein level in bronchial biopsies.</p> <p>Conclusions</p> <p>This is the first study linking <it>MRP1 </it>SNPs with lung function and inflammatory markers in COPD patients, suggesting a role of <it>MRP1 </it>SNPs in the severity of COPD in addition to their association with MRP1 protein level in bronchial biopsies.</p

CT-Based Local Distribution Metric Improves Characterization of COPD

Parametric response mapping (PRM) of paired CT lung images has been shown to improve the phenotyping of COPD by allowing for the visualization and quantification of non-emphysematous air trapping component, referred to as functional small airways disease (fSAD). Although promising, large variability in the standard method for analyzing PRM(fSAD) has been observed. We postulate that representing the 3D PRM(fSAD) data as a single scalar quantity (relative volume of PRM(fSAD)) oversimplifies the original 3D data, limiting its potential to detect the subtle progression of COPD as well as varying subtypes. In this study, we propose a new approach to analyze PRM. Based on topological techniques, we generate 3D maps of local topological features from 3D PRM(fSAD) classification maps. We found that the surface area of fSAD (S(fSAD)) was the most robust and significant independent indicator of clinically meaningful measures of COPD. We also confirmed by micro-CT of human lung specimens that structural differences are associated with unique S(fSAD) patterns, and demonstrated longitudinal feature alterations occurred with worsening pulmonary function independent of an increase in disease extent. These findings suggest that our technique captures additional COPD characteristics, which may provide important opportunities for improved diagnosis of COPD patients

Different genes interact with particulate matter and tobacco smoke exposure in affecting lung function decline in the general population

BACKGROUND: Oxidative stress related genes modify the effects of ambient air pollution or tobacco smoking on lung function decline. The impact of interactions might be substantial, but previous studies mostly focused on main effects of single genes. OBJECTIVES: We studied the interaction of both exposures with a broad set of oxidative-stress related candidate genes and pathways on lung function decline and contrasted interactions between exposures. METHODS: For 12679 single nucleotide polymorphisms (SNPs), change in forced expiratory volume in one second (FEV(1)), FEV(1) over forced vital capacity (FEV(1)/FVC), and mean forced expiratory flow between 25 and 75% of the FVC (FEF(25-75)) was regressed on interval exposure to particulate matter >10 microm in diameter (PM10) or packyears smoked (a), additive SNP effects (b), and interaction terms between (a) and (b) in 669 adults with GWAS data. Interaction p-values for 152 genes and 14 pathways were calculated by the adaptive rank truncation product (ARTP) method, and compared between exposures. Interaction effect sizes were contrasted for the strongest SNPs of nominally significant genes (p(interaction)>0.05). Replication was attempted for SNPs with MAF<10% in 3320 SAPALDIA participants without GWAS. RESULTS: On the SNP-level, rs2035268 in gene SNCA accelerated FEV(1)/FVC decline by 3.8% (p(interaction) = 2.5x10(-6)), and rs12190800 in PARK2 attenuated FEV1 decline by 95.1 ml p(interaction) = 9.7x10(-8)) over 11 years, while interacting with PM10. Genes and pathways nominally interacting with PM10 and packyears exposure differed substantially. Gene CRISP2 presented a significant interaction with PM10 (p(interaction) = 3.0x10(-4)) on FEV(1)/FVC decline. Pathway interactions were weak. Replications for the strongest SNPs in PARK2 and CRISP2 were not successful. CONCLUSIONS: Consistent with a stratified response to increasing oxidative stress, different genes and pathways potentially mediate PM10 and tobac smoke effects on lung function decline. Ignoring environmental exposures would miss these patterns, but achieving sufficient sample size and comparability across study samples is challengin

Genome-Wide Association Analysis in Asthma Subjects Identifies SPATS2L as a Novel Bronchodilator Response Gene

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