Richard Sherva

Multiple Independent Loci at Chromosome 15q25.1 Affect Smoking Quantity: a Meta-Analysis and Comparison with Lung Cancer and COPD

Recently, genetic association findings for nicotine dependence, smoking behavior, and smoking-related diseases converged to implicate the chromosome 15q25.1 region, which includes the CHRNA5-CHRNA3-CHRNB4 cholinergic nicotinic receptor subunit genes. In particular, association with the nonsynonymous CHRNA5 SNP rs16969968 and correlates has been replicated in several independent studies. Extensive genotyping of this region has suggested additional statistically distinct signals for nicotine dependence, tagged by rs578776 and rs588765. One goal of the Consortium for the Genetic Analysis of Smoking Phenotypes (CGASP) is to elucidate the associations among these markers and dichotomous smoking quantity (heavy versus light smoking), lung cancer, and chronic obstructive pulmonary disease (COPD). We performed a meta-analysis across 34 datasets of European-ancestry subjects, including 38,617 smokers who were assessed for cigarettes-per-day, 7,700 lung cancer cases and 5,914 lung-cancer-free controls (all smokers), and 2,614 COPD cases and 3,568 COPD-free controls (all smokers). We demonstrate statistically independent associations of rs16969968 and rs588765 with smoking (mutually adjusted p-values<10−35 and <10−8 respectively). Because the risk alleles at these loci are negatively correlated, their association with smoking is stronger in the joint model than when each SNP is analyzed alone. Rs578776 also demonstrates association with smoking after adjustment for rs16969968 (p<10−6). In models adjusting for cigarettes-per-day, we confirm the association between rs16969968 and lung cancer (p<10−20) and observe a nominally significant association with COPD (p = 0.01); the other loci are not significantly associated with either lung cancer or COPD after adjusting for rs16969968. This study provides strong evidence that multiple statistically distinct loci in this region affect smoking behavior. This study is also the first report of association between rs588765 (and correlates) and smoking that achieves genome-wide significance; these SNPs have previously been associated with mRNA levels of CHRNA5 in brain and lung tissue.

Genome-wide association study of alcohol dependence: significant findings in African- and European-Americans including novel risk loci

We report a GWAS of alcohol dependence (AD) in European-American (EA) and African-American (AA) populations, with replication in independent samples of EAs, AAs and Germans. Our sample for discovery and replication was 16 087 subjects, the largest sample for AD GWAS to date. Numerous genome-wide significant (GWS) associations were identified, many novel. Most associations were population specific, but in several cases were GWS in EAs and AAs for different SNPs at the same locus,showing biological convergence across populations. We confirmed well-known risk loci mapped to alcohol-metabolizing enzyme genes, notably ADH1B (EAs: Arg48His, P=1.17 × 10−31; AAs: Arg369Cys, P=6.33 × 10−17) and ADH1C in AAs (Thr151Thr, P=4.94 × 10−10), and identified novel risk loci mapping to the ADH gene cluster on chromosome 4 and extending centromerically beyond it to include GWS associations at LOC100507053 in AAs (P=2.63 × 10−11), PDLIM5 in EAs (P=2.01 × 10−8), and METAP in AAs (P=3.35 × 10−8). We also identified a novel GWS association (1.17 × 10−10) mapped to chromosome 2 at rs1437396, between MTIF2 and CCDC88A, across all of the EA and AA cohorts, with supportive gene expression evidence, and population-specific GWS for markers on chromosomes 5, 9 and 19. Several of the novel associations implicate direct involvement of, or interaction with, genes previously identified as schizophrenia risk loci. Confirmation of known AD risk loci supports the overall validity of the study; the novel loci are worthy of genetic and biological follow-up. The findings support a convergence of risk genes (but not necessarily risk alleles) between populations, and, to a lesser extent, between psychiatric traits.

Fetal hemoglobin in sickle cell anemia: genome-wide association studies suggest a regulatory region in the 5′ olfactory receptor gene cluster

In a genome-wide association study of 848 blacks with sickle cell anemia, we identified single nucleotide polymorphisms (SNPs) associated with fetal hemoglobin concentration. The most significant SNPs in a discovery sample were tested in a replication set of 305 blacks with sickle cell anemia and in subjects with hemoglobin E or beta thalassemia trait from Thailand and Hong Kong. A novel region on chromosome 11 containing olfactory receptor genes OR51B5 and OR51B6 was identified by 6 SNPs (lowest P = 4.7E-08) and validated in the replication set. An additional olfactory receptor gene, OR51B2, was identified by a novel SNP set enrichment analysis. Genome-wide association studies also validated a previously identified SNP (rs766432) in BCL11A, a gene known to affect fetal hemoglobin levels (P = 2.6E-21) and in Thailand and Hong Kong subjects. Elements within the olfactory receptor gene cluster might play a regulatory role in gamma-globin gene expression.

Common CD36 SNPs reduce protein expression and may contribute to a protective atherogenic profile

Membrane CD36 functions in the uptake of fatty acids (FAs), oxidized lipoproteins and in signal transduction after binding these ligands. In rodents, CD36 is implicated in abnormal lipid metabolism, inflammation and atherosclerosis. In humans, CD36 variants have been identified to influence free FA and high-density lipoprotein (HDL) levels and to associate with the risk of the metabolic syndrome, coronary artery disease and stroke. In this study, 15 common lipid-associated CD36 single nucleotide polymorphisms (SNPs) were evaluated for the impact on monocyte CD36 expression (protein and transcript) in 104 African Americans. In a subset of subjects, the SNPs were tested for association with monocyte surface CD36 (n=65) and platelet total CD36 (n=57). The relationship between CD36 expression and serum HDL and very low-density lipoproteins (VLDLs) levels was also examined. After a permutation-based correction for multiple tests, four SNPs (rs1761667, rs3211909, rs3211913, rs3211938) influenced monocyte CD36 protein and two (rs3211909, rs3211938) platelet CD36. The effect of the HDL-associated SNPs on CD36 expression inversely related to the impact on serum HDL and potential causality was supported by Mendelian randomization analysis. Consistent with this, monocyte CD36 protein negatively correlated with total HDL and HDL subfractions. In contrast, positive correlations were documented between monocyte CD36 and VLDL lipid, particle number and apolipoprotein B. In conclusion, CD36 variants that reduce protein expression appear to promote a protective metabolic profile. The SNPs in this study may have predictive potential on CD36 expression and disease susceptibility in African Americans. Further studies are warranted to validate and determine whether these findings are population specific.

Variation in Nicotinic Acetylcholine Receptor Genes is Associated with Multiple Substance Dependence Phenotypes

There is shared genetic risk for dependence on multiple substances, and the nicotinic receptor gene cluster on chromosome 15 harbors multiple polymorphisms that associate to this risk. Here, we report the results of an association study with 21 SNPs genotyped across the CHRNA5, CHRNA3, and CHRNB4 loci on chromosome 15q25.1. The sample consists of a discovery set (N=1858) of European-American and African-American (AA) families, ascertained on the basis of a sibling pair with cocaine and/or opioid dependence, and a case–control replication sample (N=3388) collected for association studies of alcohol, cocaine, and opioid dependence. We tested the SNPs for association with lifetime cocaine, opioid, nicotine, and alcohol dependence. We replicated several previous findings, including associations between rs16969968 and nicotine dependence (P=0.002) and cocaine dependence (P=0.02), with opposite risk alleles for each substance. We observed these associations in AAs, which is a novel finding. The strongest association signal in either sample was between rs684513 in CHRNA5 and cocaine dependence (OR=1.43, P=0.0004) in the AA replication set. We also observed two SNPs associated with alcohol dependence, that is, rs615470 in CHRNA5 (OR=0.77, P=0.0006) and rs578776 (OR=0.78, P=0.001). The associations between CD and rs684513, AD and rs615470, and AD and rs578776 remained significant after a permutation-based correction for multiple testing. These data reinforce the importance of variation in the chromosome 15 nicotinic receptor subunit gene cluster for risk of dependence on multiple substances, although the direction of the effects may vary across substances.

Genome-Wide Association Study of Opioid Dependence: Multiple Associations Mapped to Calcium and Potassium Pathways

BACKGROUND We report a genome-wide association study (GWAS) of two populations, African-American and European-American (AA, EA) for opioid dependence (OD) in three sets of subjects, to identify pathways, genes, and alleles important in OD risk. METHODS The design employed three phases (on the basis of separate sample collections). Phase 1 included our discovery GWAS dataset consisting of 5697 subjects (58% AA) diagnosed with opioid and/or other substance dependence and control subjects. Subjects were genotyped with the Illumina OmniQuad microarray, yielding 890,000 single nucleotide polymorphisms (SNPs) suitable for analysis. Additional genotypes were imputed with the 1000 Genomes reference panel. Top-ranked findings were further evaluated in Phase 2 by incorporating information from the publicly available Study of Addiction: Genetics and Environment dataset, with GWAS data from 4063 subjects (32% AA). In Phase 3, the most significant SNPs from Phase 2 were genotyped in 2549 independent subjects (32% AA). Analyses were performed with case-control and ordinal trait designs. RESULTS Most significant results emerged from the AA subgroup. Genome-wide-significant associations (p < 5.0 × 10(-8)) were observed with SNPs from multiple loci-KCNG2*rs62103177 was most significant after combining results from datasets in every phase of the study. The most compelling results were obtained with genes involved in potassium signaling pathways (e.g., KCNC1 and KCNG2). Pathway analysis also implicated genes involved in calcium signaling and long-term potentiation. CONCLUSIONS This is the first study to identify risk variants for OD with GWAS. Our results strongly implicate risk pathways and provide insights into novel therapeutic and prevention strategies and might biologically bridge OD and other non-substance dependence psychiatric traits where similar pathways have been implicated.

A 3-bp deletion in the HBS1L-MYB intergenic region on chromosome 6q23 is associated with HbF expression

Fetal hemoglobin (HbF) is regulated as a multigenic trait. By genome-wide association study, we confirmed that HBS1L-MYB intergenic polymorphisms (HMIP) and BCL11A polymorphisms are highly associated with HbF in Chinese β-thalassemia heterozygotes. In this population, the variance in HbF resulting from the HMIP is 13.5%; that resulting from the BCL11A polymorphism is 6.4%. To identify the functional variant in HMIP, we used 1000 Genomes Project data, single nucleotide polymorphism imputation, comparisons of association results across populations, potential transcription factor binding sites, and analysis of phylogenetic conservation. Based on these studies, a hitherto unreported association between HbF expression and a 3-bp deletion, between 135 460 326 and 135 460 328 bp on chromosome 6q23 was found. This 3-bp deletion is in complete linkage disequilibrium with rs9399137, which is the single nucleotide polymorphism in HMIP most significantly associated with HbF among Chinese, Europeans, and Africans. Chromatin immunoprecipitation assays confirmed erythropoiesis-related transcription factors binding to this region in K562 cells. Based on transient expression of a luciferase reporter plasmid, the DNA fragment encompassing the 3-bp deletion polymorphism has enhancer-like activity that is further augmented by the introduction of the 3-bp deletion. This 3-bp deletion polymorphism is probably the most significant functional motif accounting for HMIP modulation of HbF in all 3 populations.

Genome-wide association study of cocaine dependence and related traits: FAM53B identified as a risk gene.

We report a genome-wide association study (GWAS) for cocaine dependence (CD) in three sets of African- and European-American subjects (AAs and EAs, respectively) to identify pathways, genes and alleles important in CD risk. The discovery GWAS data set (n=5697 subjects) was genotyped using the Illumina OmniQuad microarray (8 90 000 analyzed single-nucleotide polymorphisms (SNPs)). Additional genotypes were imputed based on the 1000 Genomes reference panel. Top-ranked findings were evaluated by incorporating information from publicly available GWAS data from 4063 subjects. Then, the most significant GWAS SNPs were genotyped in 2549 independent subjects. We observed one genome-wide-significant (GWS) result: rs2629540 at the FAM53B (‘family with sequence similarity 53, member B’) locus. This was supported in both AAs and EAs; P-value (meta-analysis of all samples)=4.28 × 10−8. The gene maps to the same chromosomal region as the maximum peak we observed in a previous linkage study. NCOR2 (nuclear receptor corepressor 2) SNP rs150954431 was associated with P=1.19 × 10−9 in the EA discovery sample. SNP rs2456778, which maps to CDK1 (‘cyclin-dependent kinase 1’), was associated with cocaine-induced paranoia in AAs in the discovery sample only (P=4.68 × 10−8). This is the first study to identify risk variants for CD using GWAS. Our results implicate novel risk loci and provide insights into potential therapeutic and prevention strategies.

Association of granulomatosis with polyangiitis (Wegener's) with HLA-DPB1*04 and SEMA6A gene variants: Evidence grom genome-wide analysis

OBJECTIVE To identify genetic determinants of granulomatosis with polyangiitis (Wegeners) (GPA). METHODS We carried out a genome-wide association study (GWAS) of 492 GPA cases and 1,506 healthy controls (white subjects of European descent), followed by replication analysis of the most strongly associated signals in an independent cohort of 528 GPA cases and 1,228 controls. RESULTS Genome-wide significant associations were identified in 32 single-nucleotide polymorphic (SNP) markers across the HLA region, the majority of which were located in the HLA-DPB1 and HLA-DPA1 genes encoding the class II major histocompatibility complex (MHC) DPβ chain 1 and DPα chain 1 proteins, respectively. Peak association signals in these 2 genes, emanating from SNPs rs9277554 (for DPβ chain 1) and rs9277341 (DPα chain 1) were strongly replicated in an independent cohort (in the combined analysis of the initial cohort and the replication cohort, P = 1.92 × 10(-50) and 2.18 × 10(-39) , respectively). Imputation of classic HLA alleles and conditional analyses revealed that the SNP association signal was fully accounted for by the classic HLA-DPB1*04 allele. An independent single SNP, rs26595, near SEMA6A (the gene for semaphorin 6A) on chromosome 5, was also associated with GPA, reaching genome-wide significance in a combined analysis of the GWAS and replication cohorts (P = 2.09 × 10(-8) ). CONCLUSION We identified the SEMA6A and HLA-DP loci as significant contributors to risk for GPA, with the HLA-DPB1*04 allele almost completely accounting for the MHC association. These two associations confirm the critical role of immunogenetic factors in the development of GPA.

Genome-wide association study of the rate of cognitive decline in Alzheimer's disease

Substantial interindividual variability exists in the disease trajectories of Alzheimers disease (AD) patients. Some decline rapidly whereas others decline slowly, and there are no known explanations for this variability. We describe the first genome‐wide association study to examine rate of cognitive decline in a sample of AD patients with longitudinal measures of cognition.