Whole genome sequencing of pharmacogenetic drug response in racially and ethnically diverse children with asthma

Asthma is the most common chronic disease of children, with significant racial/ethnic differences in prevalence, morbidity, mortality and therapeutic response. Albuterol, a bronchodilator medication, is the first-line therapy for asthma treatment worldwide. We performed the largest whole genome sequencing (WGS) pharmacogenetics study to date using data from 1,441 minority children with asthma who had extremely high or low bronchodilator drug response (BDR). We identified population-specific and shared pharmacogenetic variants associated with BDR, including genome-wide significant (p < 3.53 x 10-7) and suggestive (p < 7.06 x 10-6) loci near genes previously associated with lung capacity (DNAH5), immunity (NFKB1 and PLCB1), and β-adrenergic signaling pathways (ADAMTS3 and COX18). Functional analyses centered on NFKB1 revealed potential regulatory function of our BDR-associated SNPs in bronchial smooth muscle cells. Specifically, these variants are in linkage disequilibrium with SNPs in a functionally active enhancer, and are also expression quantitative trait loci (eQTL) for a neighboring gene, SLC39A8. Given the lack of other asthma study populations with WGS data on minority children, replication of our rare variant associations is infeasible. We attempted to replicate our common variant findings in five independent studies with GWAS data. The age-specific associations previously found in asthma and asthma-related traits suggest that the over-representation of adults in our replication populations may have contributed to our lack of statistical replication, despite the functional relevance of the NFKB1 variants demonstrated by our functional assays. Our study expands the understanding of pharmacogenetic analyses in racially/ethnically diverse populations and advances the foundation for precision medicine in at-risk and understudied minority populations. AUTHOR SUMMARY Asthma is the most common chronic disease among children. Albuterol, a bronchodilator medication, is the first-line therapy for asthma treatment worldwide. In the U.S., asthma prevalence is the highest among Puerto Ricans, intermediate among African Americans and lowest in Whites and Mexicans. Asthma disparities extend to mortality, which is four- to five-fold higher in Puerto Ricans and African Americans compared to Mexicans [1]. Puerto Ricans and African Americans, the populations with the highest asthma prevalence and death rate, also have the lowest albuterol bronchodilator drug response (BDR). We conducted the largest pharmacogenetic study using whole genome sequencing data from 1,441 minority children with asthma who had extremely high or low albuterol bronchodilator drug response. We identified population-specific and shared pharmacogenetic variants associated with BDR. Our findings help inform the direction of future development of asthma medications and our study advances the foundation of precision medicine for at-risk, yet understudied, racially/ethnically diverse populations.

[1]  C. Dolea,et al.  World Health Organization , 1949, International Organization.

[2]  S. E. Avner Beta-adrenergic bronchodilators. , 1975, Pediatric clinics of North America.

[3]  D. McFadden Quantitative Methods for Analyzing Travel Behaviour of Individuals: Some Recent Developments , 1977 .

[4]  K. Saunders,et al.  Influence of age on response to ipratropium and salbutamol in asthma. , 1981, Thorax.

[5]  P. Holloway The missing patients , 1984, British Dental Journal.

[6]  J. Kaprio,et al.  A population-based study of bronchial asthma in adult twin pairs. , 1991, Chest.

[7]  G. Pugliese,et al.  Severe Streptococcus pyogenes Infections, United Kingdom, 2003–2004 , 2008, Emerging infectious diseases.

[8]  S. Leeder,et al.  A population based study , 1993, The Medical journal of Australia.

[9]  H. Nelson β-Adrenergic Bronchodilators , 1995 .

[10]  K. Little The saga continues. , 1996, Nursing standard (Royal College of Nursing (Great Britain) : 1987).

[11]  Malcolm W Johnson The β -Adrenoceptor , 1998 .

[12]  M. Johnson The beta-adrenoceptor. , 1998, American journal of respiratory and critical care medicine.

[13]  K. Kolodner,et al.  Medications Used by Children With Asthma Living in the Inner City , 1998, Pediatrics.

[14]  J L Hankinson,et al.  Spirometric reference values from a sample of the general U.S. population. , 1999, American journal of respiratory and critical care medicine.

[15]  J. Johnson,et al.  Effects of Gender and Race on Albuterol Pharmacokinetics , 1999, Pharmacotherapy.

[16]  E. Silverman,et al.  Heterogeneity of therapeutic responses in asthma. , 2000, British medical bulletin.

[17]  Elizabeth M. Smigielski,et al.  dbSNP: the NCBI database of genetic variation , 2001, Nucleic Acids Res..

[18]  Paula Lozano,et al.  Underuse of controller medications among Medicaid-insured children with asthma. , 2002, Archives of pediatrics & adolescent medicine.

[19]  L. Palmer,et al.  Pharmacogenetics of asthma. , 2002, American journal of respiratory and critical care medicine.

[20]  T. Nicolai,et al.  Longitudinal follow‐up of the changing gender ratio in asthma from childhood to adulthood: role of delayed manifestation in girls , 2003, Pediatric allergy and immunology : official publication of the European Society of Pediatric Allergy and Immunology.

[21]  Edwin K Silverman,et al.  Lower bronchodilator responsiveness in Puerto Rican than in Mexican subjects with asthma. , 2004, American journal of respiratory and critical care medicine.

[22]  Terrence S. Furey,et al.  The UCSC Table Browser data retrieval tool , 2004, Nucleic Acids Res..

[23]  Dmitri V Zaykin,et al.  Ranks of Genuine Associations in Whole-Genome Scans , 2005, Genetics.

[24]  F. Martinez Safety of long-acting beta-agonists--an urgent need to clear the air. , 2005, The New England journal of medicine.

[25]  J. Hankinson,et al.  Interpretative strategies for lung function tests , 2005, European Respiratory Journal.

[26]  J. Leong,et al.  Long acting beta agonists and asthma , 2005 .

[27]  C. Carroll,et al.  Childhood obesity increases duration of therapy during severe asthma exacerbations* , 2006, Pediatric critical care medicine : a journal of the Society of Critical Care Medicine and the World Federation of Pediatric Intensive and Critical Care Societies.

[28]  D. Reich,et al.  Principal components analysis corrects for stratification in genome-wide association studies , 2006, Nature Genetics.

[29]  D. Conrad,et al.  A worldwide survey of haplotype variation and linkage disequilibrium in the human genome , 2006, Nature Genetics.

[30]  I. Pe’er,et al.  Optimal two‐stage genotyping designs for genome‐wide association scans , 2006, Genetic epidemiology.

[31]  P. Dorinsky,et al.  The Salmeterol Multicenter Asthma Research Trial: a comparison of usual pharmacotherapy for asthma or usual pharmacotherapy plus salmeterol. , 2006, Chest.

[32]  David M. Shade,et al.  Effect of Obesity on Clinical Presentation and Response to Treatment in Asthma , 2006, The Journal of asthma : official journal of the Association for the Care of Asthma.

[33]  E. Burchard,et al.  Association between IgE levels and asthma severity among African American, Mexican, and Puerto Rican patients with asthma. , 2007, The Journal of allergy and clinical immunology.

[34]  E. Burchard,et al.  Ethnic-Specific Differences in Bronchodilator Responsiveness Among African Americans, Puerto Ricans, and Mexicans with Asthma , 2007, The Journal of asthma : official journal of the Association for the Care of Asthma.

[35]  G. Abecasis,et al.  Optimal designs for two‐stage genome‐wide association studies , 2007, Genetic epidemiology.

[36]  David J Erle,et al.  IL-13 and epidermal growth factor receptor have critical but distinct roles in epithelial cell mucin production. , 2007, American journal of respiratory cell and molecular biology.

[37]  J. Bousquet,et al.  Global surveillance, prevention and control of chronic respiratory diseases: a comprehensive approach. , 2007 .

[38]  Yee Hwa Yang,et al.  Genome-wide profiling identifies epithelial cell genes associated with asthma and with treatment response to corticosteroids , 2007, Proceedings of the National Academy of Sciences.

[39]  M. Daly,et al.  Estimation of the multiple testing burden for genomewide association studies of nearly all common variants , 2008, Genetic epidemiology.

[40]  Clifford A. Meyer,et al.  Model-based Analysis of ChIP-Seq (MACS) , 2008, Genome Biology.

[41]  Joan E Bailey-Wilson,et al.  Establishing an adjusted p-value threshold to control the family-wide type 1 error in genome wide association studies , 2008, BMC Genomics.

[42]  W. Sadee,et al.  The human zinc transporter SLC39A8 (Zip8) is critical in zinc-mediated cytoprotection in lung epithelia. , 2008, American journal of physiology. Lung cellular and molecular physiology.

[43]  Cole Trapnell,et al.  Ultrafast and memory-efficient alignment of short DNA sequences to the human genome , 2009, Genome Biology.

[44]  I. Annesi-Maesano,et al.  Heritability and Shared Genetic Effects of Asthma and Hay Fever: An Italian Study of Young Twins , 2008, Twin Research and Human Genetics.

[45]  M. Daly,et al.  Genetic Mapping in Human Disease , 2008, Science.

[46]  P. Donnelly,et al.  A Flexible and Accurate Genotype Imputation Method for the Next Generation of Genome-Wide Association Studies , 2009, PLoS genetics.

[47]  J. González,et al.  Positionally cloned genes and age-specific effects in asthma and atopy: an international population-based cohort study (ECRHS) , 2009, Thorax.

[48]  David H. Alexander,et al.  Fast model-based estimation of ancestry in unrelated individuals. , 2009, Genome research.

[49]  J. Kramer,et al.  Balancing the benefits and risks of inhaled long-acting beta-agonists--the influence of values. , 2009, The New England journal of medicine.

[50]  John P A Ioannidis,et al.  Meta-analysis in genome-wide association studies. , 2009, Pharmacogenomics.

[51]  E. Noguchi [Genetics of allergic diseases]. , 2009, Nihon rinsho. Japanese journal of clinical medicine.

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

[53]  Scott M. Williams,et al.  The Genetic Structure and History of Africans and African Americans , 2009, Science.

[54]  S. Johnston,et al.  Targeting the NF-κB pathway in asthma and chronic obstructive pulmonary disease , 2008, Pharmacology & Therapeutics.

[55]  M. Solon,et al.  Accumulation of intraepithelial mast cells with a unique protease phenotype in T(H)2-high asthma. , 2010, Journal of Allergy and Clinical Immunology.

[56]  Michael Boehnke,et al.  LocusZoom: regional visualization of genome-wide association scan results , 2010, Bioinform..

[57]  Josyf Mychaleckyj,et al.  Robust relationship inference in genome-wide association studies , 2010, Bioinform..

[58]  S. Weiss,et al.  Decreased Response To Inhaled Steroids In Overweight And Obese Asthmatic Children , 2010, ATS 2010.

[59]  J. Hankinson,et al.  Standardisation of lung function testing: the authors' replies to readers' comments , 2010, European Respiratory Journal.

[60]  M. DePristo,et al.  A framework for variation discovery and genotyping using next-generation DNA sequencing data , 2011, Nature Genetics.

[61]  D. Altshuler,et al.  Power in the phenotypic extremes: a simulation study of power in discovery and replication of rare variants , 2011, Genetic epidemiology.

[62]  Eleazar Eskin,et al.  Random-effects model aimed at discovering associations in meta-analysis of genome-wide association studies. , 2011, American journal of human genetics.

[63]  D. Conti,et al.  Using extreme phenotype sampling to identify the rare causal variants of quantitative traits in association studies , 2011, Genetic epidemiology.

[64]  C. Lamina Digging into the extremes: a useful approach for the analysis of rare variants with continuous traits? , 2011, BMC proceedings.

[65]  David J Balding,et al.  A genome-wide meta-analysis of genetic variants associated with allergic rhinitis and grass sensitization and their interaction with birth order. , 2011, The Journal of allergy and clinical immunology.

[66]  A. Dixon Long-acting β-agonists and asthma: the saga continues. , 2011, American journal of respiratory and critical care medicine.

[67]  L. Akinbami,et al.  Asthma prevalence, health care use, and mortality: United States, 2005-2009. , 2011, National health statistics reports.

[68]  E. R. Sutherland,et al.  Impact of race on asthma treatment failures in the asthma clinical research network. , 2011, American journal of respiratory and critical care medicine.

[69]  Francisco M. De La Vega,et al.  Genomics for the world , 2011, Nature.

[70]  Gonçalo R. Abecasis,et al.  The variant call format and VCFtools , 2011, Bioinform..

[71]  M. Stoltenberg,et al.  Apical Localization of Zinc Transporter ZnT4 in Human Airway Epithelial Cells and Its Loss in a Murine Model of Allergic Airway Inflammation , 2011, Nutrients.

[72]  Pui-Yan Kwok,et al.  Design and coverage of high throughput genotyping arrays optimized for individuals of East Asian, African American, and Latino race/ethnicity using imputation and a novel hybrid SNP selection algorithm. , 2011, Genomics.

[73]  T. Yao,et al.  Associations of Age, Gender, and BMI with Prevalence of Allergic Diseases in Children: PATCH Study , 2011, The Journal of asthma : official journal of the Association for the Care of Asthma.

[74]  J. Celedón,et al.  Decreased response to inhaled steroids in overweight and obese asthmatic children. , 2011, The Journal of allergy and clinical immunology.

[75]  Marc D. Perry,et al.  ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia , 2012, Genome research.

[76]  H. Kang,et al.  Mapping Genetic Variants Associated with Beta-Adrenergic Responses in Inbred Mice , 2012, PloS one.

[77]  M. Rieder,et al.  Optimal unified approach for rare-variant association testing with application to small-sample case-control whole-exome sequencing studies. , 2012, American journal of human genetics.

[78]  Yusuke Nakamura,et al.  Genome-Wide Association Analysis in Asthma Subjects Identifies SPATS2L as a Novel Bronchodilator Response Gene , 2012, PLoS genetics.

[79]  Ryan D. Hernandez,et al.  Case-control admixture mapping in Latino populations enriches for known asthma-associated genes. , 2012, The Journal of allergy and clinical immunology.

[80]  O. Delaneau,et al.  Supplementary Information for ‘ Improved whole chromosome phasing for disease and population genetic studies ’ , 2012 .

[81]  Yong Zhang,et al.  Identifying ChIP-seq enrichment using MACS , 2012, Nature Protocols.

[82]  Hua Tang,et al.  Estimating kinship in admixed populations. , 2012, American journal of human genetics.

[83]  H. Hakonarson,et al.  Gene Network Analysis In A Pediatric Cohort Identifies Novel Lung Function Genes , 2012, ATS 2012.

[84]  G. Abecasis,et al.  Detecting and estimating contamination of human DNA samples in sequencing and array-based genotype data. , 2012, American journal of human genetics.

[85]  O. Delaneau,et al.  A linear complexity phasing method for thousands of genomes , 2011, Nature Methods.

[86]  Xihong Lin,et al.  Optimal tests for rare variant effects in sequencing association studies. , 2012, Biostatistics.

[87]  H. Hakonarson,et al.  Gene Network Analysis in a Pediatric Cohort Identifies Novel Lung Function Genes , 2012, PloS one.

[88]  E. Stahl,et al.  Polygenic heritability estimates in pharmacogenetics: focus on asthma and related phenotypes , 2013, Pharmacogenetics and genomics.

[89]  S. S. Hall Genetics: A gene of rare effect , 2013, Nature.

[90]  Dan M Roden,et al.  Genome-wide association studies in pharmacogenomics: successes and lessons , 2013, Pharmacogenetics and genomics.

[91]  C. Bustamante,et al.  RFMix: a discriminative modeling approach for rapid and robust local-ancestry inference. , 2013, American journal of human genetics.

[92]  Esteban G Burchard,et al.  Early-life air pollution and asthma risk in minority children. The GALA II and SAGE II studies. , 2013, American journal of respiratory and critical care medicine.

[93]  Jiang Gui,et al.  Diverse convergent evidence in the genetic analysis of complex disease: coordinating omic, informatic, and experimental evidence to better identify and validate risk factors , 2014, BioData Mining.

[94]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[95]  Jason H. Moore,et al.  The limits of p-values for biological data mining , 2013, BioData Mining.

[96]  Christopher R. Gignoux,et al.  Socioeconomic status and childhood asthma in urban minority youths. The GALA II and SAGE II studies. , 2013, American journal of respiratory and critical care medicine.

[97]  E. Burchard,et al.  Childhood obesity and asthma control in the GALA II and SAGE II studies. , 2013, American journal of respiratory and critical care medicine.

[98]  Bernadette A. Thomas,et al.  Years lived with disability (YLDs) for 1160 sequelae of 289 diseases and injuries 1990–2010: a systematic analysis for the Global Burden of Disease Study 2010 , 2012, The Lancet.

[99]  David Haussler,et al.  The UCSC genome browser and associated tools , 2012, Briefings Bioinform..

[100]  L. Berthiaume,et al.  Wnt acylation: seeing is believing. , 2014, Nature chemical biology.

[101]  J. Celedón,et al.  A genome-wide association study of bronchodilator response in asthmatics , 2013, The Pharmacogenomics Journal.

[102]  E. Burchard Medical research: Missing patients , 2014, Nature.

[103]  Pui-Yan Kwok,et al.  Gene-based association identifies SPATA13-AS1 as a pharmacogenomic predictor of inhaled short-acting beta-agonist response in multiple population groups , 2014, The Pharmacogenomics Journal.

[104]  Ophthalmic statistics note: the perils of dichotomising continuous variables , 2014, British Journal of Ophthalmology.

[105]  D. Lynch,et al.  DNAH5 is associated with total lung capacity in chronic obstructive pulmonary disease , 2014, Respiratory Research.

[106]  Gregory M. Cooper,et al.  A general framework for estimating the relative pathogenicity of human genetic variants , 2014 .

[107]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[108]  Ryan D. Hernandez,et al.  A genome-wide association study of bronchodilator response in Latinos implicates rare variants. , 2014, The Journal of allergy and clinical immunology.

[109]  J. Shendure,et al.  A general framework for estimating the relative pathogenicity of human genetic variants , 2014, Nature Genetics.

[110]  A. Siepel,et al.  Probabilities of Fitness Consequences for Point Mutations Across the Human Genome , 2014, Nature Genetics.

[111]  Heng Li,et al.  Toward better understanding of artifacts in variant calling from high-coverage samples , 2014, Bioinform..

[112]  E. Burchard,et al.  Nocturnal asthma and the importance of race/ethnicity and genetic ancestry. , 2014, American journal of respiratory and critical care medicine.

[113]  Andrew Carroll,et al.  WGSA: an annotation pipeline for human genome sequencing studies , 2015, Journal of Medical Genetics.

[114]  Diversity in Clinical and Biomedical Research: A promise yet to be fulfilled , 2015, bioRxiv.

[115]  H. Zheng,et al.  Performance of Genotype Imputation for Low Frequency and Rare Variants from the 1000 Genomes , 2015, PloS one.

[116]  A. Siepel,et al.  Probabilities of Fitness Consequences for Point Mutations Across the Human Genome , 2014, Nature Genetics.

[117]  Gabor T. Marth,et al.  A global reference for human genetic variation , 2015, Nature.

[118]  J. Alcorn,et al.  Diet, interleukin-17, and childhood asthma in Puerto Ricans. , 2015, Annals of allergy, asthma & immunology : official publication of the American College of Allergy, Asthma, & Immunology.

[119]  R. Wu,et al.  Genome-wide association study of short-acting β2-agonists. A novel genome-wide significant locus on chromosome 2 near ASB3. , 2015, American journal of respiratory and critical care medicine.

[120]  E. Burchard,et al.  Moving toward true inclusion of racial/ethnic minorities in federally funded studies. A key step for achieving respiratory health equality in the United States. , 2015, American journal of respiratory and critical care medicine.

[121]  Jordan A. Ramilowski,et al.  Transcriptome analysis of controlled and therapy-resistant childhood asthma reveals distinct gene expression profiles. , 2015, The Journal of allergy and clinical immunology.

[122]  N. Powe,et al.  Diversity in Clinical and Biomedical Research: A Promise Yet to Be Fulfilled , 2015, bioRxiv.

[123]  Carson C Chow,et al.  Second-generation PLINK: rising to the challenge of larger and richer datasets , 2014, GigaScience.

[124]  Tom R. Gaunt,et al.  Improved imputation of low-frequency and rare variants using the UK10K haplotype reference panel , 2015, Nature Communications.

[125]  David Kulp,et al.  A guide to genome‐wide association analysis and post‐analytic interrogation , 2015, Statistics in medicine.

[126]  Wei Zheng,et al.  Addressing Population‐Specific Multiple Testing Burdens in Genetic Association Studies , 2015, Annals of human genetics.

[127]  Danny S. Park,et al.  Genetic and socioeconomic study of mate choice in Latinos reveals novel assortment patterns , 2015, Proceedings of the National Academy of Sciences.

[128]  Wen-Miin Liang,et al.  Genetic variants in PLCB4/PLCB1 as susceptibility loci for coronary artery aneurysm formation in Kawasaki disease in Han Chinese in Taiwan , 2015, Scientific Reports.

[129]  J. Lupski,et al.  Non-coding genetic variants in human disease. , 2015, Human molecular genetics.

[130]  G. Canino,et al.  Expression Quantitative Trait Loci (eQTL) Mapping in Puerto Rican Children , 2015, PloS one.

[131]  S. Gabriel,et al.  Phenotypic extremes in rare variant study designs , 2015, European Journal of Human Genetics.

[132]  Alan M. Kwong,et al.  Next-generation genotype imputation service and methods , 2016, Nature Genetics.

[133]  Christopher R. Gignoux,et al.  Making Precision Medicine Socially Precise. Take a Deep Breath. , 2016, American journal of respiratory and critical care medicine.

[134]  E. Burchard,et al.  Novel genetic risk factors for asthma in African American children: Precision Medicine and the SAGE II Study , 2016, bioRxiv.

[135]  E. Trevathan,et al.  Racial disparity in the association between body mass index and self-reported asthma in children: a population-based study , 2016, The Journal of asthma : official journal of the Association for the Care of Asthma.

[136]  Alan M. Kwong,et al.  A reference panel of 64,976 haplotypes for genotype imputation , 2015, Nature Genetics.

[137]  Kiang Liu,et al.  Pulmonary Disease and Age at Immigration among Hispanics. Results from the Hispanic Community Health Study/Study of Latinos. , 2016, American journal of respiratory and critical care medicine.

[138]  S. Fullerton,et al.  Genomics is failing on diversity , 2016, Nature.

[139]  John Novembre,et al.  Visualizing the geography of genetic variants , 2016, bioRxiv.

[140]  E. Burchard,et al.  Perceived Discrimination Associated With Asthma and Related Outcomes in Minority Youth: The GALA II and SAGE II Studies , 2017, Chest.