James W. Baurley

Meta-analysis of genome-wide association studies of asthma in ethnically diverse North American populations

Asthma is a common disease with a complex risk architecture including both genetic and environmental factors. We performed a meta-analysis of North American genome-wide association studies of asthma in 5,416 individuals with asthma (cases) including individuals of European American, African American or African Caribbean, and Latino ancestry, with replication in an additional 12,649 individuals from the same ethnic groups. We identified five susceptibility loci. Four were at previously reported loci on 17q21, near IL1RL1, TSLP and IL33, but we report for the first time, to our knowledge, that these loci are associated with asthma risk in three ethnic groups. In addition, we identified a new asthma susceptibility locus at PYHIN1, with the association being specific to individuals of African descent (P = 3.9 × 10−9). These results suggest that some asthma susceptibility loci are robust to differences in ancestry when sufficiently large samples sizes are investigated, and that ancestry-specific associations also contribute to the complex genetic architecture of asthma.

Genome-wide Association Analysis Identifies PDE4D as an Asthma-Susceptibility Gene

Asthma, a chronic airway disease with known heritability, affects more than 300 million people around the world. A genome-wide association (GWA) study of asthma with 359 cases from the Childhood Asthma Management Program (CAMP) and 846 genetically matched controls from the Illumina ICONdb public resource was performed. The strongest region of association seen was on chromosome 5q12 in PDE4D. The phosphodiesterase 4D, cAMP-specific (phosphodiesterase E3 dunce homolog, Drosophila) gene (PDE4D) is a regulator of airway smooth-muscle contractility, and PDE4 inhibitors have been developed as medications for asthma. Allelic p values for top SNPs in this region were 4.3 x 10(-07) for rs1588265 and 9.7 x 10(-07) for rs1544791. Replications were investigated in ten independent populations with different ethnicities, study designs, and definitions of asthma. In seven white and Hispanic replication populations, two PDE4D SNPs had significant results with p values less than 0.05, and five had results in the same direction as the original population but had p values greater than 0.05. Combined p values for 18,891 white and Hispanic individuals (4,342 cases) in our replication populations were 4.1 x 10(-04) for rs1588265 and 9.2 x 10(-04) for rs1544791. In three black replication populations, which had different linkage disequilibrium patterns than the other populations, original findings were not replicated. Further study of PDE4D variants might lead to improved understanding of the role of PDE4D in asthma pathophysiology and the efficacy of PDE4 inhibitor medications.

TSLP polymorphisms are associated with asthma in a sex-specific fashion

To cite this article: Hunninghake GM, Soto‐Quirós ME, Avila L, Kim HP, Lasky‐Su J, Rafaels N, Ruczinski I, Beaty TH, Mathias RA, Barnes KC, Wilk JB, O’Connor GT, James Gauderman W, Vora H, Baurley JW, Gilliland F, Liang C, Sylvia JS, Klanderman BJ, Sharma SS, Himes BE, Bossley CJ, Israel E, Raby BA, Bush A, Choi AM, Weiss ST, Celedón JC. TSLP polymorphisms are associated with asthma in a sex‐specific fashion. Allergy 2010; 65: 1566–1575.

PDE11A associations with asthma: Results of a genome-wide association scan

Andrew T. DeWan, PhDa, Elizabeth W. Triche, PhDc, Xuming Xu, PhDa, Ling-I Hsu, MPHa, Connie Zhao, PhDd, Kathleen Belanger, PhDb, Karen Hellenbrand, MPHb, Saffron A. G. Willis-Owen, PhDe, Miriam Moffatt, PhDe, William O. C. Cookson, MD, PhDe, Blanca E. Himes, PhDf, Scott T. Weiss, MDf, W. James Gauderman, PhDg, James W. Baurley, MSg, Frank Gilliland, MD, PhDg, Jemma B. Wilk, DSch, George T. O’Connor, MDh, David P. Strachan, PhDi, Josephine Hoh, PhDa, and Michael B. Bracken, PhDb

Discovery of Complex Pathways from Observational Data

Unraveling complex interactions has been a challenge in epidemiologic research. We introduce a pathway modeling framework that discovers plausible pathways from observational data, and allows estimation of both the net effect of the pathway and the types of interactions occurring among genetic or environmental risk factors. Each discovered pathway structure links combinations of observed variables through intermediate latent nodes to a final node, the outcome. Biologic knowledge can be readily applied in this framework as a prior on pathway structure to give preference to more biologically plausible models, thereby providing more precise estimation of Bayes factors for pathways of greatest interest by Markov Chain Monte Carlo (MCMC) methods.Data were simulated for binary inputs of which only a subset was involved in different pathway topologies. Our algorithm was then used to recover the pathway from the simulated data. The posterior distributions of inputs, pairwise and higher-order interactions, and topologies were obtained by MCMC methods. The evidence in favor of a particular pathway or interaction was summarized using Bayes factors. Our method can correctly identify the risk factors and interactions involved in the simulated pathway. We apply our framework to an asthma case-control data set with polymorphisms in 12 genes.

Detecting Gene-Environment Interactions in Genome-Wide Association Data

Despite the importance of gene‐environment (G×E) interactions in the etiology of common diseases, little work has been done to develop methods for detecting these types of interactions in genome‐wide association study data. This was the focus of Genetic Analysis Workshop 16 Group 10 contributions, which introduced a variety of new methods for the detection of G×E interactions in both case‐control and family‐based data using both cross‐sectional and longitudinal study designs. Many of these contributions detected significant G×E interactions. Although these interactions have not yet been confirmed, the results suggest the importance of testing for interactions. Issues of sample size, quantifying the environmental exposure, longitudinal data analysis, family‐based analysis, selection of the most powerful analysis method, population stratification, and computational expense with respect to testing G×E interactions are discussed. Genet. Epidemiol. 33 (Suppl. 1):S68–S73, 2009.

Genetic risk factors for orofacial clefts in Central Africans and Southeast Asians

Genome‐wide association studies (GWAS) for orofacial clefts have identified several susceptibility regions, but have largely focused on non‐Hispanic White populations in developed countries. We performed a targeted genome‐wide study of single nucleotide polymorphisms (SNPs) in exons using the Illumina HumanExome+ array with custom fine mapping of 16 cleft susceptibility regions in three underserved populations: Congolese (87 case‐mother, 210 control‐mother pairs), Vietnamese (131 case‐parent trios), and Filipinos (42 case‐mother, 99 control‐mother pairs). All cases were children with cleft lip with or without cleft palate. Families were recruited from local hospitals and parental exposures were collected using interviewer‐administered questionnaires. We used logistic regression models for case‐control analyses, family‐based association tests for trios, and fixed‐effect meta‐analyses to determine individual SNP effects corrected for multiple testing. Of the 16 known susceptibility regions tested, SNPs in four regions reached statistical significance in one or more of these populations: 1q32.2 (IRF6), 10q25.3 (VAX1), and 17q22 (NOG). Due to different linkage disequilibrium patterns, significant SNPs in these regions differed between the Vietnamese and Filipino populations from the index SNP selected from previous GWAS studies. Among Africans, there were no significant associations identified for any of the susceptibility regions. rs10787738 near VAX1 (P = 4.98E−3) and rs7987165 (P = 6.1E−6) were significant in the meta‐analysis of all three populations combined. These results confirm several known susceptibility regions and identify novel risk alleles in understudied populations.

Genome-Wide Association of the Laboratory-Based Nicotine Metabolite Ratio in Three Ancestries

Introduction: Metabolic enzyme variation and other patient and environmental characteristics influence smoking behaviors, treatment success, and risk of related disease. Population-specific variation in metabolic genes contributes to challenges in developing and optimizing pharmacogenetic interventions. We applied a custom genome-wide genotyping array for addiction research (Smokescreen), to three laboratory-based studies of nicotine metabolism with oral or venous administration of labeled nicotine and cotinine, to model nicotine metabolism in multiple populations. The trans-3′-hydroxycotinine/cotinine ratio, the nicotine metabolite ratio (NMR), was the nicotine metabolism measure analyzed. Methods: Three hundred twelve individuals of self-identified European, African, and Asian American ancestry were genotyped and included in ancestry-specific genome-wide association scans (GWAS) and a meta-GWAS analysis of the NMR. We modeled natural-log transformed NMR with covariates: principal components of genetic ancestry, age, sex, body mass index, and smoking status. Results: African and Asian American NMRs were statistically significantly (P values ≤ 5E-5) lower than European American NMRs. Meta-GWAS analysis identified 36 genome-wide significant variants over a 43 kilobase pair region at CYP2A6 with minimum P = 2.46E-18 at rs12459249, proximal to CYP2A6. Additional minima were located in intron 4 (rs56113850, P = 6.61E-18) and in the CYP2A6-CYP2A7 intergenic region (rs34226463, P = 1.45E-12). Most (34/36) genome-wide significant variants suggested reduced CYP2A6 activity; functional mechanisms were identified and tested in knowledge-bases. Conditional analysis resulted in intergenic variants of possible interest (P values < 5E-5). Conclusions: This meta-GWAS of the NMR identifies CYP2A6 variants, replicates the top-ranked single nucleotide polymorphism from a recent Finnish meta-GWAS of the NMR, identifies functional mechanisms, and provides pan-continental population biomarkers for nicotine metabolism. Implications: This multiple ancestry meta-GWAS of the laboratory study-based NMR provides novel evidence and replication for genome-wide association of CYP2A6 single nucleotide and insertion–deletion polymorphisms. We identify three regions of genome-wide significance: proximal, intronic, and distal to CYP2A6. We replicate the top-ranking single nucleotide polymorphism from a recent GWAS of the NMR in Finnish smokers, identify a functional mechanism for this intronic variant from in silico analyses of RNA-seq data that is consistent with CYP2A6 expression measured in postmortem lung and liver, and provide additional support for the intergenic region between CYP2A6 and CYP2A7.

Drug Metabolizing Enzyme and Transporter Gene Variation, Nicotine Metabolism, Prospective Abstinence, and Cigarette Consumption.

The Nicotine Metabolite Ratio (NMR, ratio of trans-3’-hydroxycotinine and cotinine), has previously been associated with CYP2A6 activity, response to smoking cessation treatments, and cigarette consumption. We searched for drug metabolizing enzyme and transporter (DMET) gene variation associated with the NMR and prospective abstinence in 2,946 participants of laboratory studies of nicotine metabolism and of clinical trials of smoking cessation therapies. Stage I was a meta-analysis of the association of 507 common single nucleotide polymorphisms (SNPs) at 173 DMET genes with the NMR in 449 participants of two laboratory studies. Nominally significant associations were identified in ten genes after adjustment for intragenic SNPs; CYP2A6 and two CYP2A6 SNPs attained experiment-wide significance adjusted for correlated SNPs (CYP2A6 P ACT=4.1E-7, rs4803381 P ACT=4.5E-5, rs1137115, P ACT=1.2E-3). Stage II was mega-regression analyses of 10 DMET SNPs with pretreatment NMR and prospective abstinence in up to 2,497 participants from eight trials. rs4803381 and rs1137115 SNPs were associated with pretreatment NMR at genome-wide significance. In post-hoc analyses of CYP2A6 SNPs, we observed nominally significant association with: abstinence in one pharmacotherapy arm; cigarette consumption among all trial participants; and lung cancer in four case:control studies. CYP2A6 minor alleles were associated with reduced NMR, CPD, and lung cancer risk. We confirmed the major role that CYP2A6 plays in nicotine metabolism, and made novel findings with respect to genome-wide significance and associations with CPD, abstinence and lung cancer risk. Additional multivariate analyses with patient variables and genetic modeling will improve prediction of nicotine metabolism, disease risk and smoking cessation treatment prognosis.

Parental risk factors for oral clefts among Central Africans, Southeast Asians, and Central Americans

Background Several lifestyle and environmental exposures have been suspected as risk factors for oral clefts, although few have been convincingly demonstrated. Studies across global diverse populations could offer additional insight given varying types and levels of exposures. Methods We performed an international case–control study in the Democratic Republic of the Congo (133 cases, 301 controls), Vietnam (75 cases, 158 controls), the Philippines (102 cases, 152 controls), and Honduras (120 cases, 143 controls). Mothers were recruited from hospitals and their exposures were collected from interviewer‐administered questionnaires. We used logistic regression modeling to estimate odds ratios (OR) and 95% confidence intervals (CI). Results Family history of clefts was strongly associated with increased risk (maternal: OR = 4.7; 95% CI, 3.0–7.2; paternal: OR = 10.5; 95% CI, 5.9–18.8; siblings: OR = 5.3; 95% CI, 1.4–19.9). Advanced maternal age (5 year OR = 1.2; 95% CI, 1.0–1.3), pregestational hypertension (OR = 2.6; 95% CI, 1.3–5.1), and gestational seizures (OR = 2.9; 95% CI, 1.1–7.4) were statistically significant risk factors. Lower maternal (secondary school OR = 1.6; 95% CI, 1.2–2.2; primary school OR = 2.4, 95% CI, 1.6–2.8) and paternal education (OR = 1.9; 95% CI, 1.4–2.5; and OR = 1.8; 95% CI, 1.1–2.9, respectively) and paternal tobacco smoking (OR = 1.5, 95% CI, 1.1–1.9) were associated with an increased risk. No other significant associations between maternal and paternal factors were found; some environmental factors including rural residency, indoor cooking with wood, chemicals and water source appeared to be associated with an increased risk in adjusted models. Conclusion Our study represents one of the first international studies investigating risk factors for clefts among multiethnic underserved populations. Our findings suggest a multifactorial etiology including both maternal and paternal factors. Birth Defects Research (Part A) 103:863–879, 2015.