The Overlap of Lung Tissue Transcriptome of Smoke Exposed Mice with Human Smoking and COPD

Genome-wide mRNA profiling in lung tissue from human and animal models can provide novel insights into the pathogenesis of chronic obstructive pulmonary disease (COPD). While 6 months of smoke exposure are widely used, shorter durations were also reported. The overlap of short term and long-term smoke exposure in mice is currently not well understood, and their representation of the human condition is uncertain. Lung tissue gene expression profiles of six murine smoking experiments (n = 48) were obtained from the Gene Expression Omnibus (GEO) and analyzed to identify the murine smoking signature. The 'human smoking' gene signature containing 386 genes was previously published in the lung eQTL study (n = 1,111). A signature of mild COPD containing 7 genes was also identified in the same study. The lung tissue gene signature of 'severe COPD' (n = 70) contained 4,071 genes and was previously published. We detected 3,723 differentially expressed genes in the 6 month-exposure mice datasets (FDR <0.1). Of those, 184 genes (representing 48% of human smoking) and 1,003 (representing 27% of human COPD) were shared with the human smoking-related genes and the COPD severity-related genes, respectively. There was 4-fold over-representation of human and murine smoking-related genes (P = 6.7 × 10-26) and a 1.4 fold in the severe COPD -related genes (P = 2.3 × 10-12). There was no significant enrichment of the mice and human smoking-related genes in mild COPD signature. These data suggest that murine smoke models are strongly representative of molecular processes of human smoking but less of COPD

Genome-Wide Interaction Analysis of Air Pollution Exposure and Childhood Asthma with Functional Follow-up

Rationale: The evidence supporting an association between traffic-related air pollution exposure and incident childhood asthma is inconsistent and may depend on genetic factors. Objectives: To identify gene–environment interaction effects on childhood asthma using genome-wide single-nucleotide polymorphism (SNP) data and air pollution exposure. Identified loci were further analyzed at epigenetic and transcriptomic levels. Methods: We used land use regression models to estimate individual air pollution exposure (represented by outdoor NO2 levels) at the birth address and performed a genome-wide interaction study for doctors’ diagnoses of asthma up to 8 years in three European birth cohorts (n = 1,534) with look-up for interaction in two separate North American cohorts, CHS (Children’s Health Study) and CAPPS/SAGE (Canadian Asthma Primary Prevention Study/Study of Asthma, Genetics and Environment) (n = 1,602 and 186 subjects, respectively). We assessed expression quantitative trait locus effects in human lung specimens and blood, as well as associations among air pollution exposure, methylation, and transcriptomic patterns. Measurements and Main Results: In the European cohorts, 186 SNPs had an interaction P < 1 × 10−4 and a look-up evaluation of these disclosed 8 SNPs in 4 loci, with an interaction P < 0.05 in the large CHS study, but not in CAPPS/SAGE. Three SNPs within adenylate cyclase 2 (ADCY2) showed the same direction of the interaction effect and were found to influence ADCY2 gene expression in peripheral blood (P = 4.50 × 10−4). One other SNP with P < 0.05 for interaction in CHS, rs686237, strongly influenced UDP-Gal:betaGlcNAc β-1,4-galactosyltransferase, polypeptide 5 (B4GALT5) expression in lung tissue (P = 1.18 × 10−17). Air pollution exposure was associated with differential discs, large homolog 2 (DLG2) methylation and expression. Conclusions: Our results indicated that gene–environment interactions are important for asthma development and provided supportive evidence for interaction with air pollution for ADCY2, B4GALT5, and DLG2

Genetic regulation of gene expression of MIF family members in lung tissue

Macrophage migration inhibitory factor (MIF) is a cytokine found to be associated with chronic obstructive pulmonary disease (COPD). However, there is no consensus on how MIF levels differ in COPD compared to control conditions and there are no reports on MIF expression in lung tissue. Here we studied gene expression of members of the MIF family MIF, D-Dopachrome Tautomerase (DDT) and DDT-like (DDTL) in a lung tissue dataset with 1087 subjects and identified single nucleotide polymorphisms (SNPs) regulating their gene expression. We found higher MIF and DDT expression in COPD patients compared to non-COPD subjects and found 71 SNPs significantly influencing gene expression of MIF and DDTL. Furthermore, the platform used to measure MIF (microarray or RNAseq) was found to influence the splice variants detected and subsequently the direction of the SNP effects on MIF expression. Among the SNPs found to regulate MIF expression, the major LD block identified was linked to rs5844572, a SNP previously found to be associated with lower diffusion capacity in COPD. This suggests that MIF may be contributing to the pathogenesis of COPD, as SNPs that influence MIF expression are also associated with symptoms of COPD. Our study shows that MIF levels are affected not only by disease but also by genetic diversity (i.e. SNPs). Since none of our significant eSNPs for MIF or DDTL have been described in GWAS for COPD or lung function, MIF expression in COPD patients is more likely a consequence of disease-related factors rather than a cause of the disease

Whole exome re-sequencing implicates CCDC38 and cilia structure and function in resistance to smoking related airflow obstruction

Chronic obstructive pulmonary disease (COPD) is a leading cause of global morbidity and mortality and, whilst smoking remains the single most important risk factor, COPD risk is heritable. Of 26 independent genomic regions showing association with lung function in genome-wide association studies, eleven have been reported to show association with airflow obstruction. Although the main risk factor for COPD is smoking, some individuals are observed to have a high forced expired volume in 1 second (FEV1) despite many years of heavy smoking. We # hypothesised that these ‘‘resistant smokers’’ may harbour variants which protect against lung function decline caused by smoking and provide insight into the genetic determinants of lung health. We undertook whole exome re sequencing of 100 heavy smokers who had healthy lung function given their age, sex, height and smoking history and applied three complementary approaches to explore the genetic architecture of smoking resistance. Firstly, we identified novel functional variants in the ‘‘resistant smokers’’ and looked for enrichment of these novel variants within biological pathways. Secondly, we undertook association testing of all exonic variants individually with two independent control sets. Thirdly, we undertook gene-based association testing of all exonic variants. Our strongest signal of association with smoking resistance for a non-synonymous SNP was for rs10859974 (P = 2.3461024) in CCDC38, a gene which has previously been reported to show association with FEV1/FVC, and we demonstrate moderate expression of CCDC38 in bronchial epithelial cells. We identified an enrichment of novel putatively functional variants in genes related to cilia structure and function in resistant smokers. Ciliary function abnormalities are known to be associated with both smoking and reduced mucociliary clearance in patients with COPD. We suggest that genetic influences on the development or function of cilia in the bronchial epithelium may affect growth of cilia or the extent of damage caused by tobacco smoke

The Pharmacogenomics of Inhaled Corticosteroids and Lung Function Decline in COPD

Inhaled corticosteroids (ICS) are widely prescribed for patients with chronic obstructive pulmonary disease (COPD), yet with variable outcomes and adverse reactions which may be genetically determined. The primary aim of the study was to identify the genetic determinants for FEV1 changes related to ICS therapy. In the Lung Health Study 2 (LHS-2), 1116 COPD patients were randomised to the ICS, triamcinolone acetonide (n=559), or placebo (n=557) with spirometry performed every 6 months for 3 years. We performed a pharmacogenomic genome-wide association study (GWAS) for the genotype-by-ICS treatment effect on 3 years of forced expiratory volume in 1 s (FEV1) changes (estimated as slope) in 802 genotyped LHS-2 participants. Replication was performed in 199 COPD patients randomised to the ICS, fluticasone or placebo. A total of five loci showed genotype-by-ICS interaction at p&lt;5×10-6; of these, SNP rs111720447 on chromosome 7 was replicated (discovery p=4.8×10-6, replication p=5.9×10-5) with the same direction of interaction effect. ENCODE data revealed that in glucocorticoid treated (dexamethasone) A549 alveolar cell line, glucocorticoid receptor binding sites were located near SNP rs111720447. In stratified analyses of LHS-2, genotype at SNP rs111720447 was significantly associated with rate of FEV1 decline in patients taking ICS (C allele beta=56.35 mL·year-1, 95% confidence interval (CI)=29.96, 82.76 mL·yr-1) and also in patients who were assigned to placebo, though the relationship was weaker and in the opposite direction than that in the ICS group (C allele beta=-27.57 mL·year-1, 95% CI=-53.27, -1.87 mL·yr-1). The study uncovered genetic factors associated with FEV1 changes related to ICS in COPD patients, which may provide new insight on the potential biology of steroid responsiveness in COPD.</p

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

Susceptibility to chronic mucus hypersecretion, a genome wide association study

Background: Chronic mucus hypersecretion (CMH) is associated with an increased frequency of respiratory infections, excess lung function decline, and increased hospitalisation and mortality rates in the general population. It is associated with smoking, but it is unknown why only a minority of smokers develops CMH. A plausible explanation for this phenomenon is a predisposing genetic constitution. Therefore, we performed a genome wide association (GWA) study of CMH in Caucasian populations. Methods: GWA analysis was performed in the NELSON-study using the Illumina 610 array, followed by replication and meta-analysis in 11 additional cohorts. In total 2,704 subjects with, and 7,624 subjects without CMH were included, all current or former heavy smokers (≥20 pack-years). Additional studies were performed to test the functional relevance of the most significant single nucleotide polymorphism (SNP). Results: A strong association with CMH, consistent across all cohorts, was observed with rs6577641 (p = 4.25x10-6, OR = 1.17), located in intron 9 of the special AT-rich sequence-binding protein 1 locus (SATB1) on chromosome 3. The risk allele (G) was associated with higher mRNA expression of SATB1 (4.3x10 -9) in lung tissue. Presence of CMH was associated with increased SATB1 mRNA expression in bronchial biopsies from COPD patients. SATB1 expression was induced during differentiation of primary human bronchial epithelial cells in culture. Conclusions: Our findings, that SNP rs6577641 is associated with CMH in multiple cohorts and is a cis-eQTL for SATB1, together with our additional observation that SATB1 expression increases during epithelial differentiation provide suggestive evidence that SATB1 is a gene that affects CMH