Genome-wide association study across European and African American ancestries identifies a SNP in DNMT3B contributing to nicotine dependence.

Cigarette smoking is a leading cause of preventable mortality worldwide. Nicotine dependence, which reduces the likelihood of quitting smoking, is a heritable trait with firmly established associations with sequence variants in nicotine acetylcholine receptor genes and at other loci. To search for additional loci, we conducted a genome-wide association study (GWAS) meta-analysis of nicotine dependence, totaling 38,602 smokers (28,677 Europeans/European Americans and 9925 African Americans) across 15 studies. In this largest-ever GWAS meta-analysis for nicotine dependence and the largest-ever cross-ancestry GWAS meta-analysis for any smoking phenotype, we reconfirmed the well-known CHRNA5-CHRNA3-CHRNB4 genes and further yielded a novel association in the DNA methyltransferase gene DNMT3B. The intronic DNMT3B rs910083-C allele (frequency=44-77%) was associated with increased risk of nicotine dependence at P=3.7 × 10-8 (odds ratio (OR)=1.06 and 95% confidence interval (CI)=1.04-1.07 for severe vs mild dependence). The association was independently confirmed in the UK Biobank (N=48,931) using heavy vs never smoking as a proxy phenotype (P=3.6 × 10-4, OR=1.05, and 95% CI=1.02-1.08). Rs910083-C is also associated with increased risk of squamous cell lung carcinoma in the International Lung Cancer Consortium (N=60,586, meta-analysis P=0.0095, OR=1.05, and 95% CI=1.01-1.09). Moreover, rs910083-C was implicated as a cis-methylation quantitative trait locus (QTL) variant associated with higher DNMT3B methylation in fetal brain (N=166, P=2.3 × 10-26) and a cis-expression QTL variant associated with higher DNMT3B expression in adult cerebellum from the Genotype-Tissue Expression project (N=103, P=3.0 × 10-6) and the independent Brain eQTL Almanac (N=134, P=0.028). This novel DNMT3B cis-acting QTL variant highlights the importance of genetically influenced regulation in brain on the risks of nicotine dependence, heavy smoking and consequent lung cancer.

Collaborative meta-analysis finds no evidence of a strong interaction between stress and 5-HTTLPR genotype contributing to the development of depression.

The hypothesis that the S allele of the 5-HTTLPR serotonin transporter promoter region is associated with increased risk of depression, but only in individuals exposed to stressful situations, has generated much interest, research and controversy since first proposed in 2003. Multiple meta-analyses combining results from heterogeneous analyses have not settled the issue. To determine the magnitude of the interaction and the conditions under which it might be observed, we performed new analyses on 31 data sets containing 38 802 European ancestry subjects genotyped for 5-HTTLPR and assessed for depression and childhood maltreatment or other stressful life events, and meta-analysed the results. Analyses targeted two stressors (narrow, broad) and two depression outcomes (current, lifetime). All groups that published on this topic prior to the initiation of our study and met the assessment and sample size criteria were invited to participate. Additional groups, identified by consortium members or self-identified in response to our protocol (published prior to the start of analysis) with qualifying unpublished data, were also invited to participate. A uniform data analysis script implementing the protocol was executed by each of the consortium members. Our findings do not support the interaction hypothesis. We found no subgroups or variable definitions for which an interaction between stress and 5-HTTLPR genotype was statistically significant. In contrast, our findings for the main effects of life stressors (strong risk factor) and 5-HTTLPR genotype (no impact on risk) are strikingly consistent across our contributing studies, the original study reporting the interaction and subsequent meta-analyses. Our conclusion is that if an interaction exists in which the S allele of 5-HTTLPR increases risk of depression only in stressed individuals, then it is not broadly generalisable, but must be of modest effect size and only observable in limited situations.

Largest GWAS of PTSD (N=20 070) yields genetic overlap with schizophrenia and sex differences in heritability.

The Psychiatric Genomics Consortium-Posttraumatic Stress Disorder group (PGC-PTSD) combined genome-wide case-control molecular genetic data across 11 multiethnic studies to quantify PTSD heritability, to examine potential shared genetic risk with schizophrenia, bipolar disorder, and major depressive disorder and to identify risk loci for PTSD. Examining 20 730 individuals, we report a molecular genetics-based heritability estimate (h2SNP) for European-American females of 29% that is similar to h2SNP for schizophrenia and is substantially higher than h2SNP in European-American males (estimate not distinguishable from zero). We found strong evidence of overlapping genetic risk between PTSD and schizophrenia along with more modest evidence of overlap with bipolar and major depressive disorder. No single-nucleotide polymorphisms (SNPs) exceeded genome-wide significance in the transethnic (overall) meta-analysis and we do not replicate previously reported associations. Still, SNP-level summary statistics made available here afford the best-available molecular genetic index of PTSD-for both European- and African-American individuals-and can be used in polygenic risk prediction and genetic correlation studies of diverse phenotypes. Publication of summary statistics for ∼10 000 African Americans contributes to the broader goal of increased ancestral diversity in genomic data resources. In sum, the results demonstrate genetic influences on the development of PTSD, identify shared genetic risk between PTSD and other psychiatric disorders and highlight the importance of multiethnic/racial samples. As has been the case with schizophrenia and other complex genetic disorders, larger sample sizes are needed to identify specific risk loci.

Evidence of CNIH3 involvement in opioid dependence.

Opioid dependence, a severe addictive disorder and major societal problem, has been demonstrated to be moderately heritable. We conducted a genome-wide association study in Comorbidity and Trauma Study data comparing opioid-dependent daily injectors (N=1167) with opioid misusers who never progressed to daily injection (N=161). The strongest associations, observed for CNIH3 single-nucleotide polymorphisms (SNPs), were confirmed in two independent samples, the Yale-Penn genetic studies of opioid, cocaine and alcohol dependence and the Study of Addiction: Genetics and Environment, which both contain non-dependent opioid misusers and opioid-dependent individuals. Meta-analyses found five genome-wide significant CNIH3 SNPs. The A allele of rs10799590, the most highly associated SNP, was robustly protective (P=4.30E-9; odds ratio 0.64 (95% confidence interval 0.55-0.74)). Epigenetic annotation predicts that this SNP is functional in fetal brain. Neuroimaging data from the Duke Neurogenetics Study (N=312) provide evidence of this SNP's in vivo functionality; rs10799590 A allele carriers displayed significantly greater right amygdala habituation to threat-related facial expressions, a phenotype associated with resilience to psychopathology. Computational genetic analyses of physical dependence on morphine across 23 mouse strains yielded significant correlations for haplotypes in CNIH3 and functionally related genes. These convergent findings support CNIH3 involvement in the pathophysiology of opioid dependence, complementing prior studies implicating the α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) glutamate system.

Rare, low frequency and common coding variants in CHRNA5 and their contribution to nicotine dependence in European and African Americans.

The common nonsynonymous variant rs16969968 in the α5 nicotinic receptor subunit gene (CHRNA5) is the strongest genetic risk factor for nicotine dependence in European Americans and contributes to risk in African Americans. To comprehensively examine whether other CHRNA5 coding variation influences nicotine dependence risk, we performed targeted sequencing on 1582 nicotine-dependent cases (Fagerström Test for Nicotine Dependence score⩾4) and 1238 non-dependent controls, with independent replication of common and low frequency variants using 12 studies with exome chip data. Nicotine dependence was examined using logistic regression with individual common variants (minor allele frequency (MAF)⩾0.05), aggregate low frequency variants (0.05>MAF⩾0.005) and aggregate rare variants (MAF<0.005). Meta-analysis of primary results was performed with replication studies containing 12 174 heavy and 11 290 light smokers. Next-generation sequencing with 180 × coverage identified 24 nonsynonymous variants and 2 frameshift deletions in CHRNA5, including 9 novel variants in the 2820 subjects. Meta-analysis confirmed the risk effect of the only common variant (rs16969968, European ancestry: odds ratio (OR)=1.3, P=3.5 × 10(-11); African ancestry: OR=1.3, P=0.01) and demonstrated that three low frequency variants contributed an independent risk (aggregate term, European ancestry: OR=1.3, P=0.005; African ancestry: OR=1.4, P=0.0006). The remaining 22 rare coding variants were associated with increased risk of nicotine dependence in the European American primary sample (OR=12.9, P=0.01) and in the same risk direction in African Americans (OR=1.5, P=0.37). Our results indicate that common, low frequency and rare CHRNA5 coding variants are independently associated with nicotine dependence risk. These newly identified variants likely influence the risk for smoking-related diseases such as lung cancer.

ADH1B is associated with alcohol dependence and alcohol consumption in populations of European and African ancestry.

A coding variant in alcohol dehydrogenase 1B (ADH1B) (rs1229984) that leads to the replacement of Arg48 with His48 is common in Asian populations and reduces their risk for alcoholism, but because of very low allele frequencies the effects in European or African populations have been difficult to detect. We genotyped and analyzed this variant in three large European and African-American case-control studies in which alcohol dependence was defined by the Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition (DSM-IV) criteria, and demonstrated a strong protective effect of the His48 variant (odds ratio (OR) 0.34, 95% confidence interval (CI) 0.24, 0.48) on alcohol dependence, with genome-wide significance (6.6 × 10(-10)). The hypothesized mechanism of action involves an increased aversive reaction to alcohol; in keeping with this hypothesis, the same allele is strongly associated with a lower maximum number of drinks in a 24-hour period (lifetime), with P=3 × 10(-13). We also tested the effects of this allele on the development of alcoholism in adolescents and young adults, and demonstrated a significantly protective effect. This variant has the strongest effect on risk for alcohol dependence compared with any other tested variant in European populations.

Juxtaposed regions of extensive and minimal linkage disequilibrium in human Xq25 and Xq28.

Linkage disequilibrium (LD), or the non-random association of alleles, is poorly understood in the human genome. Population genetic theory suggests that LD is determined by the age of the markers, population history, recombination rate, selection and genetic drift. Despite the uncertainties in determining the relative contributions of these factors, some groups have argued that LD is a simple function of distance between markers. Disease-gene mapping studies and a simulation study gave differing predictions on the degree of LD in isolated and general populations. In view of the discrepancies between theory and experimental observations, we constructed a high-density SNP map of the Xq25-Xq28 region and analysed the male genotypes and haplotypes across this region for LD in three populations. The populations included an outbred European sample (CEPH males) and isolated population samples from Finland and Sardinia. We found two extended regions of strong LD bracketed by regions with no evidence for LD in all three samples. Haplotype analysis showed a paucity of haplotypes in regions of strong LD. Our results suggest that, in this region of the X chromosome, LD is not a monotonic function of the distance between markers, but is more a property of the particular location in the human genome.

A sex-adjusted and age-adjusted genome screen for nested alcohol dependence diagnoses.

Alcohol dependence is a complex disorder with a substantial genetic contribution to susceptibility. The Collaborative Study on the Genetics of Alcoholism is a multi-site study whose purpose is to detect, localize, and characterize genes contributing to this susceptibility. Previous linkage analyses of the trait of alcohol dependence in Collaborative Study on the Genetics of Alcoholism have used affected sib-pair methods with a dichotomous phenotype definition. In contrast, the analysis in this paper uses a sex-adjusted and age-adjusted multiple threshold liability model. The use of such a model, in that it includes unaffected as well as as affected subjects and in that it utilizes the differential severity of a diagnosis scale, should heuristically be more powerful than a straight affected sib-pair analysis. Three regions of interest are found on chromosome 1 (lod 5.17), chromosome 4 (lod 3.46), and chromosome 8 (lod 4.31). The region on chromosome 1 near the marker D1S532 is in the region previously reported as linked to alcohol dependence and correlated phenotypes in this dataset. The region on chromosome 4 near the alcohol dehydrogenase gene cluster has been reported to be linked to alcohol dependence in other studies, as well as to the alcohol consumption phenotype 'Maximum Number of Drinks in a 24-Hour Period' in this dataset. The region on chromosome 8 near the marker D8S1988 is homologous to a section of rat chromosome 5 to which an alcohol consumption phenotype has been linked.

Replication of ZNF804A gene variant associations with risk of heroin addiction.

Heroin addiction is heritable, but few specific genetic variants have been reproducibly associated with this disease. The zinc finger protein 804A (ZNF804A) gene is a biologically plausible susceptibility gene for heroin addiction, given its function as a transcription factor in human brain. Novel associations of two common ZNF804A single nucleotide polymorphisms (SNPs), rs7597593 and rs1344706, with heroin addiction have been reported in Han Chinese. Both SNPs have also been implicated for regulating ZNF804A expression in human brain, including the addiction-relevant dorsolateral prefrontal cortex. In this independent replication study, we tested the rs7597593 and rs1344706 SNP genotypes and their corresponding haplotypes for association with heroin addiction using cases drawn from the Urban Health Study and population controls: total N = 10 757 [7095 European Americans (EAs) and 3662 African Americans (AAs)]. We independently replicated both ZNF804A SNP associations in EAs: the rs7597593-T (P = 0.016) and rs1344706-A (P = 0.029) alleles both being associated with increased risk of heroin addiction, consistent with the prior report. Neither SNP was associated in AAs alone, but meta-analysis across both ancestry groups resulted in significant associations for rs1344706-A [P = 0.016, odds ratio (95% confidence interval) = 1.13 (1.02-1.25)] and its haplotype with rs7597593-T [P = 0.0067, odds ratio (95% confidence interval) = 1.16 (1.04-1.29)]. By showing consistent associations across independent studies and diverse ancestry groups, our study provides evidence that these two ZNF804A SNPs and their risk haplotype are among the few replicable genetic associations with heroin addiction.

Genome-wide meta-analysis reveals common splice site acceptor variant in CHRNA4 associated with nicotine dependence.

We conducted a 1000 Genomes-imputed genome-wide association study (GWAS) meta-analysis for nicotine dependence, defined by the Fagerström Test for Nicotine Dependence in 17 074 ever smokers from five European-ancestry samples. We followed up novel variants in 7469 ever smokers from five independent European-ancestry samples. We identified genome-wide significant association in the alpha-4 nicotinic receptor subunit (CHRNA4) gene on chromosome 20q13: lowest P=8.0 × 10(-9) across all the samples for rs2273500-C (frequency=0.15; odds ratio=1.12 and 95% confidence interval=1.08-1.17 for severe vs mild dependence). rs2273500-C, a splice site acceptor variant resulting in an alternate CHRNA4 transcript predicted to be targeted for nonsense-mediated decay, was associated with decreased CHRNA4 expression in physiologically normal human brains (lowest P=7.3 × 10(-4)). Importantly, rs2273500-C was associated with increased lung cancer risk (N=28 998, odds ratio=1.06 and 95% confidence interval=1.00-1.12), likely through its effect on smoking, as rs2273500-C was no longer associated with lung cancer after adjustment for smoking. Using criteria for smoking behavior that encompass more than the single 'cigarettes per day' item, we identified a common CHRNA4 variant with important regulatory properties that contributes to nicotine dependence and smoking-related consequences.

A multi-incident, Old-Order Amish family with PD.

BACKGROUND: PD is largely a sporadic condition of unknown etiology, but specific inherited mutations are a cause of PD. OBJECTIVE: To describe a large, multi-incident Amish pedigree with PD. METHODS: Case ascertainment, calculation of population prevalence, and calculation of kinship coefficients (a measure of relatedness between two individuals) for affecteds and subjects in a large kindred with PD were conducted. Sequencing of genes with known mutations sufficient to cause PD and marker-by-marker haplotype analysis in chromosomal regions flanking previously described genes with known mutations were performed. RESULTS: The authors have examined 113 members of this pedigree and classified 67 as normal (no evidence of PD), 17 as clinically definite PD, 6 as clinically probable PD, and 23 as clinically possible PD. The mean age at onset of the clinically definite subjects was 56.7 years. The phenotype in this family is typical of idiopathic PD, including rest tremor, rigidity, bradykinesia, postural instability, and response to levodopa. In addition, dementia occurred in six of the clinically definite subjects, and many subjects experienced levodopa-related motor complications including wearing off and dopa-induced dyskinesias. In the index Amish community, a minimum prevalence of PD in the population 40 years and older of 552/100,000 was calculated. The mean kinship coefficient in the subjects with PD and those with PD by history (0.036) was higher (p = 0.007) than in a group of age-matched normal Amish control subjects (0.016), providing evidence that PD is inherited in this family. Sequence analysis did not detect any mutations in known PD genes. No single haplotype cosegregates with the disease in any of the chromosomal regions previously found to be linked to PD, and no marker in these regions exhibits increased homozygosity among definite PD cases. CONCLUSIONS: PD in this community is more common than in the general population, and this increased prevalence may be due in part to a novel gene(s).

A genome screen of maximum number of drinks as an alcoholism phenotype.

The Collaborative Study on the Genetics of Alcoholism (COGA) is a multicenter research program to detect and map susceptibility genes for alcohol dependence and related phenotypes. The measure M of "maximum number of drinks consumed in a 24-hour period" is closely related to alcoholism diagnosis in this dataset and provides a quantitative measure to grade nonalcoholic individuals. Twin studies have shown log(M) to have a heritability of approximately 50%. Genome screens for this trait were performed in two distinct genotyped samples (wave 1 and wave 2), and in the combined sample. MAPMAKER/SIBS was used to carry out Haseman-Elston based regression analyses. On chromosome 4, an unweighted all-pairs multipoint LOD of 2.2 was obtained between D4S2407 and D4S1628 in wave 1; in wave 2, the region flanked by D4S2404 and D4S2407 gave a LOD of 1.5. In the combined sample, the maximal LOD was 3.5 very close to D4S2407. This evidence for linkage is in the region of the alcohol dehydrogenase gene cluster on chromosome 4. These findings on chromosome 4 are consistent with a prior report from COGA in which strictly defined nonalcoholic subjects in wave 1 were analyzed. The present analysis on log(M) allows more individuals to be included and thus is potentially more powerful.

Genome screen for platelet monoamine oxidase (MAO) activity.

To identify loci involved in the control of platelet monoamine oxidase B (MAO-B) activity, a genomewide linkage screen was performed using 291 markers in 148 nuclear families containing a total of 1,008 nonindependent sib-pairs. Participants were genotyped and their platelet MAO-B activity levels were measured as part of the Collaborative Study on the Genetics of Alcoholism (COGA). Sib-pair analysis using Haseman-Elston regression was carried out with two programs. Two-point analysis on all pairs with SIBPAL indicated three markers with p-values below 0.01: D6S1018 (p = 0.0004), D2S1328 (p = 0.008), and D2S408 (p = 0.003). MAPMAKER/SIBS multipoint analyses using independent pairs(N = 409) gave maximal lod scores of 2. 0 on chromosome 6 and 1.1 and 1.4 for the two regions on chromosome 2. These results are consistent with linkage, but do not provide definitive evidence. We are currently creating a denser map in these regions and have begun genotyping a second sample in COGA.

A bootstrapped commingling analysis of platelet monoamine oxidase activity levels corrected for cigarette smoking.

Monoamine oxidase (MAO) activity levels have been suggested as a possible biological marker for alcohol dependence and abuse, as well as for schizophrenia and other psychiatric conditions. Using platelet MAO activities in the Collaborative Study on the Genetics of Alcoholism data set, we applied bootstrapping methods as a novel way to test for admixture in families. This bootstrapping involved resampling in family units and hypothesis testing of the resampled datasets for commingling in the distribution of MAO activity levels. Prior to commingling analysis, we used linear models to find covariates of greatest effect on MAO activity levels. While an alcoholism diagnosis was significant in men (n = 1151, P < 0.0001) and women (n = 1254, P = 0.0003), the effect lost significance after controlling for cigarette smoking, indicating alcoholism and smoking behavior to be highly confounded. When smoking histories were compared, former smokers had levels (mean = 7.1) closer to those who never smoked (mean = 7.0) than to current smokers (mean = 5.4). Furthermore, current daily smoking and time since smoking cessation were significantly related to MAO levels, indicating smoking probably has a direct effect on MAO levels, rather than the reverse. These results suggest that studies using MAO levels as a biological marker should consider smoking as an important covariate. Finally, admixture was found in MAO levels controlled for smoking and sex, possibly indicating a major genetic locus; this confirms previous evidence for admixture.

Multiple cholinergic nicotinic receptor genes affect nicotine dependence risk in African and European Americans.

Several independent studies show that the chromosome 15q25.1 region, which contains the CHRNA5-CHRNA3-CHRNB4 gene cluster, harbors variants strongly associated with nicotine dependence, other smoking behaviors, lung cancer and chronic obstructive pulmonary disease. We investigated whether variants in other cholinergic nicotinic receptor subunit (CHRN) genes affect the risk of nicotine dependence in a new sample of African Americans (AAs) (N = 710). We also analyzed this AA sample together with a European American (EA) sample (N = 2062, 1608 of which have been previously studied), allowing for differing effects in the two populations. Cases are current nicotine-dependent smokers and controls are non-dependent smokers. Variants in or near CHRND-CHRNG, CHRNA7 and CHRNA10 show modest association with nicotine dependence risk in the AA sample. In addition, CHRNA4, CHRNB3-CHRNA6 and CHRNB1 show association in at least one population. CHRNG and CHRNA4 harbor single nucleotide polymorphisms (SNPs) that have opposite directions of effect in the two populations. In each of the population samples, these loci substantially increase the trait variation explained, although no loci meet Bonferroni-corrected significance in the AA sample alone. The trait variation explained by three key associated SNPs in CHRNA5-CHRNA3-CHRNB4 is 1.9% in EAs and also 1.9% in AAs; this increases to 4.5% in EAs and 7.3% in AAs when we add six variants representing associations at other CHRN genes. Multiple nicotinic receptor subunit genes outside chromosome 15q25 are likely to be important in the biological processes and development of nicotine dependence, and some of these risks may be shared across diverse populations.

Definition of the phenotype.

Definition of the phenotype is a key issue in designing any genetic study whose goal is to detect disease genes. This chapter describes strategies to increase the power to detect susceptibility loci for complex diseases. A narrowly defined disease phenotype can offer advantages over broad definitions. Studies of clinical disease can also benefit from judicious selection of endophenotypes and related quantitative traits for analysis. The effect of diagnostic and measurement error is also discussed; power is maximized when strategies to reduce error are incorporated into a study design.

The lod score method.

The lod score method originated in a seminal article by Newton Morton in 1955. The method is broadly concerned with issues of power and the posterior probability of linkage, ensuring that a reported linkage has a high probability of being a true linkage. In addition, the method is sequential, so that pedigrees or lod curves may be combined from published reports to pool data for analysis. This approach has been remarkably successful for 50 years in identifying disease genes for Mendelian disorders. After discussing these issues, we consider the situation for complex disorders, where the maximum lod score (MLS) statistic shares some of the advantages of the traditional lod score approach but is limited by unknown power and the lack of sharing of the primary data needed to optimally combine analytic results. We may still learn from the lod score method as we explore new methods in molecular biology and genetic analysis to utilize the complete human DNA sequence and the cataloging of all human genes.

Clustering methods applied to allele sharing data.

Here we focus on using clustering methods to disentangle the interacting factors that lead to the presentation of complex diseases. Relative pairs are placed in discrete subgroups, or classes, based upon their pattern of allele sharing at a sequence of markers and on concomitant risk factors. The relationship between the locus information and the affectation status of the relative pairs within each subgroup then can be assessed. Cluster analysis (CLA) and latent class analysis (LCA) were applied to sibling allele sharing data from GAW11 simulated data, and to an existing Alzheimer's disease (AD) dataset. Both methods were able to identify markers linked to all 3 disease loci in the GAW11 data. LCA and CLA also replicated regions of chromosomes identified in an analysis of the AD data using affected-sib-pair methods. These analyses indicate that classification tools may be useful for detecting susceptibility genes for complex traits.

Covariates in linkage analysis using sibling and cousin pairs.

Using simulated data from GAW 12, problem 2, we further develop a novel technique to detect and use significant covariates in linkage analysis. The method, first introduced by Rice et al. [Genet Epidemiol 17(Suppl. 1):S691-5, 1999], uses logistic regression to model perturbation in sharing as a function of covariate levels. The original method allows use of all sib pairs (concordant affected, concordant unaffected, and discordant). Here we extend this method to include cousin pairs in analysis.

Covariates in linkage analysis.

We apply a novel technique to detect significant covariates in linkage analysis using a logistic regression approach. An overall test of linkage is first performed to determine whether there is significant perturbation from the expected 50% sharing under the hypothesis of no linkage; if the overall test is significant, the importance of the individual covariate is assessed. In addition, association analyses were performed. These methods were applied to simulated data from multiple populations, and detected correct marker linkages and associations. No population heterogeneity was detected. These methods have the advantages of using all sib pairs and of providing a formal test for heterogeneity across populations.