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.34×10−4) 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.

[1]  J. Sunyer,et al.  The relation of circulating YKL-40 to levels and decline of lung function in adult life. , 2013, Respiratory medicine.

[2]  L. A. Daniels,et al.  Primary ciliary dyskinesia. Recent advances in diagnostics, genetics, and characterization of clinical disease. , 2013, American journal of respiratory and critical care medicine.

[3]  L. Wain,et al.  GSTCD and INTS12 Regulation and Expression in the Human Lung , 2013, PloS one.

[4]  Michael Boehnke,et al.  Recommended Joint and Meta‐Analysis Strategies for Case‐Control Association Testing of Single Low‐Count Variants , 2013, Genetic epidemiology.

[5]  Don D. Sin,et al.  Refining Susceptibility Loci of Chronic Obstructive Pulmonary Disease with Lung eqtls , 2013, PloS one.

[6]  D. MacArthur,et al.  Negligible impact of rare autoimmune-locus coding-region variants on missing heritability , 2013, Nature.

[7]  V. Kim,et al.  Chronic bronchitis and chronic obstructive pulmonary disease. , 2013, American journal of respiratory and critical care medicine.

[8]  T. Spector,et al.  Genes Contributing to Pain Sensitivity in the Normal Population: An Exome Sequencing Study , 2012, PLoS genetics.

[9]  David C. Nickle,et al.  Lung eQTLs to Help Reveal the Molecular Underpinnings of Asthma , 2012, PLoS genetics.

[10]  Edwin K Silverman,et al.  Genome-wide association studies identify CHRNA5/3 and HTR4 in the development of airflow obstruction. , 2012, American journal of respiratory and critical care medicine.

[11]  S. Jimenez,et al.  Chitinase 1 Is a Biomarker for and Therapeutic Target in Scleroderma-Associated Interstitial Lung Disease That Augments TGF-β1 Signaling , 2012, The Journal of Immunology.

[12]  C. Li,et al.  Genome‐wide search for replicable risk gene regions in alcohol and nicotine co‐dependence , 2012, American journal of medical genetics. Part B, Neuropsychiatric genetics : the official publication of the International Society of Psychiatric Genetics.

[13]  Eleftheria Zeggini,et al.  ARIEL and AMELIA: Testing for an Accumulation of Rare Variants Using Next-Generation Sequencing Data , 2012, Human Heredity.

[14]  E. Falconnet,et al.  Static respiratory cilia associated with mutations in Dnahc11/DNAH11: A mouse model of PCD , 2012, Human mutation.

[15]  T. Beaty,et al.  Genetic association between human chitinases and lung function in COPD , 2012, Human Genetics.

[16]  Blair H. Smith,et al.  Effect of five genetic variants associated with lung function on the risk of chronic obstructive lung disease, and their joint effects on lung function. , 2011, American journal of respiratory and critical care medicine.

[17]  Christian Gieger,et al.  Genome-wide association and large scale follow-up identifies 16 new loci influencing lung function , 2011, Nature Genetics.

[18]  S. Ware,et al.  Spectrum of clinical diseases caused by disorders of primary cilia. , 2011, Proceedings of the American Thoracic Society.

[19]  Xihong Lin,et al.  Rare-variant association testing for sequencing data with the sequence kernel association test. , 2011, American journal of human genetics.

[20]  N. Laird,et al.  The association of genome-wide significant spirometric loci with chronic obstructive pulmonary disease susceptibility. , 2011, American journal of respiratory cell and molecular biology.

[21]  Jeffrey A Lieberman,et al.  Genome-Wide Pharmacogenomic Study of Neurocognition As an Indicator of Antipsychotic Treatment Response in Schizophrenia , 2011, Neuropsychopharmacology.

[22]  K. Anderson,et al.  The coiled-coil domain containing protein CCDC40 is essential for motile cilia function and left-right axis formation , 2011, Nature Genetics.

[23]  J. Belmont,et al.  CCDC39 is required for assembly of inner dynein arms and the dynein regulatory complex and for normal ciliary motility in humans and dogs , 2011, Nature Genetics.

[24]  Inês Barroso,et al.  Genome-wide association study identifies five loci associated with lung function , 2010, Nature Genetics.

[25]  E. Lundberg,et al.  Towards a knowledge-based Human Protein Atlas , 2010, Nature Biotechnology.

[26]  C. Brightling,et al.  Ciliary dysfunction and ultrastructural abnormalities are features of severe asthma. , 2010, The Journal of allergy and clinical immunology.

[27]  M. DePristo,et al.  The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. , 2010, Genome research.

[28]  C. Klaver,et al.  Mutations in IMPG2, encoding interphotoreceptor matrix proteoglycan 2, cause autosomal-recessive retinitis pigmentosa. , 2010, American journal of human genetics.

[29]  Ayellet V. Segrè,et al.  Common Inherited Variation in Mitochondrial Genes Is Not Enriched for Associations with Type 2 Diabetes or Related Glycemic Traits , 2010, PLoS genetics.

[30]  H. Hakonarson,et al.  ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data , 2010, Nucleic acids research.

[31]  Tariq Ahmad,et al.  Meta-analysis and imputation refines the association of 15q25 with smoking quantity , 2010, Nature Genetics.

[32]  Ming D. Li,et al.  Genome-wide meta-analyses identify multiple loci associated with smoking behavior , 2010, Nature Genetics.

[33]  C. Gieger,et al.  Sequence variants at CHRNB3–CHRNA6 and CYP2A6 affect smoking behavior , 2010, Nature Genetics.

[34]  T. Ideker,et al.  Functional genomic screen for modulators of ciliogenesis and cilium length , 2010, Nature.

[35]  A. Hofman,et al.  Meta-analyses of genome-wide association studies identify multiple loci associated with pulmonary function , 2010, Nature Genetics.

[36]  C. Gieger,et al.  Sequence variants at CHRNB 3 – CHRNA 6 and CYP 2 A 6 affect smoking behavior , 2010 .

[37]  R. Crystal,et al.  Smoking Is Associated with Shortened Airway Cilia , 2009, PloS one.

[38]  Judy H. Cho,et al.  Finding the missing heritability of complex diseases , 2009, Nature.

[39]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[40]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[41]  Scott T. Weiss,et al.  A Genome-Wide Association Study of Pulmonary Function Measures in the Framingham Heart Study , 2009, PLoS genetics.

[42]  Brad T. Sherman,et al.  Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources , 2008, Nature Protocols.

[43]  Tim D Spector,et al.  The UK Adult Twin Registry (TwinsUK) , 2006, Twin Research and Human Genetics.

[44]  J Vestbo,et al.  Developing COPD: a 25 year follow up study of the general population , 2006, Thorax.

[45]  C. D. Mathers,et al.  Chronic obstructive pulmonary disease: current burden and future projections , 2006, European Respiratory Journal.

[46]  M. Sanderson,et al.  The Ciliary Rootlet Maintains Long-Term Stability of Sensory Cilia , 2005, Molecular and Cellular Biology.

[47]  E. Pierce,et al.  The Retinitis Pigmentosa 1 Protein Is a Photoreceptor Microtubule-Associated Protein , 2004, The Journal of Neuroscience.

[48]  P. Paré,et al.  The nature of small-airway obstruction in chronic obstructive pulmonary disease. , 2004, The New England journal of medicine.

[49]  B. Hogan,et al.  Novel Role for Netrins in Regulating Epithelial Behavior during Lung Branching Morphogenesis , 2004, Current Biology.

[50]  Rafael A Irizarry,et al.  Exploration, normalization, and summaries of high density oligonucleotide array probe level data. , 2003, Biostatistics.

[51]  Miguel Armengot,et al.  Mutations in the DNAH11 (axonemal heavy chain dynein type 11) gene cause one form of situs inversus totalis and most likely primary ciliary dyskinesia , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[52]  E. Broadfield,et al.  Prospective study of diet and decline in lung function in a general population. , 2002, American journal of respiratory and critical care medicine.

[53]  H. Lehrach,et al.  Mutations in DNAH5 cause primary ciliary dyskinesia and randomization of left–right asymmetry , 2002, Nature Genetics.

[54]  P. Burkhard,et al.  Coiled coils: a highly versatile protein folding motif. , 2001, Trends in cell biology.

[55]  B Rosner,et al.  Genetic epidemiology of severe, early-onset chronic obstructive pulmonary disease. Risk to relatives for airflow obstruction and chronic bronchitis. , 1998, American journal of respiratory and critical care medicine.

[56]  I. Pavord,et al.  Dietary antioxidant vitamin intake and lung function in the general population. , 1995, American journal of respiratory and critical care medicine.

[57]  J. Bernaudin,et al.  Ciliary abnormalities in bronchial epithelium of smokers, ex-smokers, and nonsmokers. , 1995, American journal of respiratory and critical care medicine.

[58]  C. Larsson.  Natural history and life expectancy in severe alpha1-antitrypsin deficiency, Pi Z. , 2009, Acta medica Scandinavica.