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.

A rare missense mutation in CHRNA4 associates with smoking behavior and its consequences.

Using Icelandic whole-genome sequence data and an imputation approach we searched for rare sequence variants in CHRNA4 and tested them for association with nicotine dependence. We show that carriers of a rare missense variant (allele frequency=0.24%) within CHRNA4, encoding an R336C substitution, have greater risk of nicotine addiction than non-carriers as assessed by the Fagerstrom Test for Nicotine Dependence (P=1.2 × 10(-4)). The variant also confers risk of several serious smoking-related diseases previously shown to be associated with the D398N substitution in CHRNA5. We observed odds ratios (ORs) of 1.7-2.3 for lung cancer (LC; P=4.0 × 10(-4)), chronic obstructive pulmonary disease (COPD; P=9.3 × 10(-4)), peripheral artery disease (PAD; P=0.090) and abdominal aortic aneurysms (AAAs; P=0.12), and the variant associates strongly with the early-onset forms of LC (OR=4.49, P=2.2 × 10(-4)), COPD (OR=3.22, P=2.9 × 10(-4)), PAD (OR=3.47, P=9.2 × 10(-3)) and AAA (OR=6.44, P=6.3 × 10(-3)). Joint analysis of the four smoking-related diseases reveals significant association (P=6.8 × 10(-5)), particularly for early-onset cases (P=2.1 × 10(-7)). Our results are in agreement with functional studies showing that the human α4ß2 isoform of the channel containing R336C has less sensitivity for its agonists than the wild-type form following nicotine incubation.

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 common biological basis of obesity and nicotine addiction.

Smoking influences body weight such that smokers weigh less than non-smokers and smoking cessation often leads to weight increase. The relationship between body weight and smoking is partly explained by the effect of nicotine on appetite and metabolism. However, the brain reward system is involved in the control of the intake of both food and tobacco. We evaluated the effect of single-nucleotide polymorphisms (SNPs) affecting body mass index (BMI) on smoking behavior, and tested the 32 SNPs identified in a meta-analysis for association with two smoking phenotypes, smoking initiation (SI) and the number of cigarettes smoked per day (CPD) in an Icelandic sample (N=34,216 smokers). Combined according to their effect on BMI, the SNPs correlate with both SI (r=0.019, P=0.00054) and CPD (r=0.032, P=8.0 × 10(-7)). These findings replicate in a second large data set (N=127,274, thereof 76,242 smokers) for both SI (P=1.2 × 10(-5)) and CPD (P=9.3 × 10(-5)). Notably, the variant most strongly associated with BMI (rs1558902-A in FTO) did not associate with smoking behavior. The association with smoking behavior is not due to the effect of the SNPs on BMI. Our results strongly point to a common biological basis of the regulation of our appetite for tobacco and food, and thus the vulnerability to nicotine addiction and obesity.

Genome-wide association analysis of coffee drinking suggests association with CYP1A1/CYP1A2 and NRCAM.

Coffee consumption is a model for addictive behavior. We performed a meta-analysis of genome-wide association studies (GWASs) on coffee intake from 8 Caucasian cohorts (N=18 176) and sought replication of our top findings in a further 7929 individuals. We also performed a gene expression analysis treating different cell lines with caffeine. Genome-wide significant association was observed for two single-nucleotide polymorphisms (SNPs) in the 15q24 region. The two SNPs rs2470893 and rs2472297 (P-values=1.6 × 10(-11) and 2.7 × 10(-11)), which were also in strong linkage disequilibrium (r(2)=0.7) with each other, lie in the 23-kb long commonly shared 5' flanking region between CYP1A1 and CYP1A2 genes. CYP1A1 was found to be downregulated in lymphoblastoid cell lines treated with caffeine. CYP1A1 is known to metabolize polycyclic aromatic hydrocarbons, which are important constituents of coffee, whereas CYP1A2 is involved in the primary metabolism of caffeine. Significant evidence of association was also detected at rs382140 (P-value=3.9 × 10(-09)) near NRCAM-a gene implicated in vulnerability to addiction, and at another independent hit rs6495122 (P-value=7.1 × 10(-09))-an SNP associated with blood pressure-in the 15q24 region near the gene ULK3, in the meta-analysis of discovery and replication cohorts. Our results from GWASs and expression analysis also strongly implicate CAB39L in coffee drinking. Pathway analysis of differentially expressed genes revealed significantly enriched ubiquitin proteasome (P-value=2.2 × 10(-05)) and Parkinson's disease pathways (P-value=3.6 × 10(-05)).

Screen for CACNA1A and ATP1A2 mutations in sporadic hemiplegic migraine patients.

The aim of this study was to investigate the involvement of the CACNA1A and ATP1A2 gene in a population-based sample of sporadic hemiplegic migraine (SHM). Patients with SHM (n = 105) were identified in a nationwide search in the Danish population. We sequenced all exons and promoter regions of the CACNA1A and ATP1A2 genes in 100 patients with SHM to search for possible SHM mutations. Novel DNA variants were discovered in eight SHM patients, four in exons of the CACNA1A gene and four in exons of the ATP1A2 gene. Six of the variants were considered non-pathogenic. The causal role of the two remaining DNA variants is unknown until functional studies have been made or independent genetic evidence is discovered. Only very few DNA variants were identified in 100 SHM patients, and regardless of whether the identified variants are causal the CACNA1A and ATP1A2 genes are not major genes in SHM.

CDKN2A mutations and melanoma risk in the Icelandic population.

BACKGROUND: Germline CDKN2A mutations have been observed in 20-40% of high risk, melanoma prone families; however, little is known about their prevalence in population based series of melanoma cases and controls. METHODS: We resequenced the CDKN2A gene, including the p14ARF variant and promoter regions, in approximately 703 registry ascertained melanoma cases and 691 population based controls from Iceland, a country in which the incidence of melanoma has increased rapidly. RESULTS: We identified a novel germline variant, G89D, that was strongly associated with increased melanoma risk and appeared to be an Icelandic founder mutation. The G89D variant was present in about 2% of Icelandic invasive cutaneous malignant melanoma cases. Relatives of affected G89D carriers were at significantly increased risk of melanoma, head and neck cancers, and pancreatic carcinoma compared to relatives of other melanoma patients. Nineteen other germline variants were identified, but none conferred an unequivocal risk of melanoma. CONCLUSIONS: This population based study of Icelandic melanoma cases and controls showed a frequency of disease related CDKN2A mutant alleles ranging from 0.7% to 1.0%, thus expanding our knowledge about the frequency of CDKN2A mutations in different populations. In contrast to North America and Australia where a broad spectrum of mutations was observed at a similar frequency, in Iceland, functional CDKN2A mutations consist of only one or two different variants. Additional genetic and/or environmental factors are likely critical for explaining the high incidence rates for melanoma in Iceland. This study adds to the geographic regions for which population based estimates of CDKN2A mutation frequencies are available.

The genetic spectrum of a population-based sample of familial hemiplegic migraine.

Familial hemiplegic migraine (FHM) is a rare subtype of migraine with aura and transient hemiplegia. FHM mutations are known in three genes, the CACNA1A (FHM1) gene, the ATP1A2 (FHM2) and the SCN1A (FHM3) gene and seem to have an autosomal-dominant mode of inheritance. The aim of this study was to search for FHM mutations in FHM families identified through a screen of the Danish population of 5.2 million people. FHM patients were diagnosed according to the International Classification of Headache Disorders and all FHM patients had a physical and neurological examination by a physician. A total of 147 FHM patients from 44 different families were identified; 43 FHM families participated in this study. Linkage analysis of these families shows clear linkage to the FHM locus (FHM1) on chromosome 19, supportive linkage to the FHM2 locus whereas no linkage was found to the FHM3 locus. Furthermore, we sequenced all exons and promoter regions of the CACNA1A and ATP1A2 genes and screened for the Q1489K mutation in the SCN1A gene. CACNA1A gene mutations were identified in three of the FHM families, two known FHM mutations, R583Q and T666M and one novel C1369Y mutation. Three FHM families were identified with novel mutations in the ATP1A2 gene; a family with a V138A mutation, a family with a R202Q mutation and a family with a R763C mutation. None of the Danish FHM families have the Q1489K mutation in the SCN1A gene. Our study shows that only 14% (6/42) of FHM families in the general Danish population have exonic FHM mutations in the CACNA1A or ATP1A2 gene. The families we identified with FHM mutations in the CACNA1A and ATP1A2 genes were extended, multiple affected families whereas the remaining FHM families were smaller. The existence of many small families in the Danish FHM cohort may reflect less bias in FHM family ascertainment and/or more locus heterogeneity than described previously.

Association between IL-1beta/TNF-alpha-induced glucocorticoid-sensitive changes in multiple gene expression and altered responsiveness in airway smooth muscle.

The pleiotropic cytokines interleukin (IL)-1beta and tumor necrosis factor (TNF)-alpha have been implicated in the pathophysiology of asthma. To elucidate the role of these cytokines in the pro-asthmatic state, the effects of IL-1beta and TNF-alpha on airway smooth muscle (ASM) responsiveness and ASM expression of multiple genes, assessed by high-density oligonucleotide array analysis, were examined in the absence and presence of the glucocorticoid dexamethasone (DEX). Administration of IL-1beta/TNF-alpha increased ASM contractility to acetylcholine and impaired ASM relaxation to isoproterenol. These pro-asthmatic- like changes in ASM responsiveness were associated with IL-1beta/ TNF-alpha-induced mRNA expression of a host of proinflammatory genes that regulate transcription, cytokines and chemokines, cellular adhesion molecules, and various signal transduction molecules that regulate ASM responsiveness. In the presence of DEX, the changes induced in ASM responsiveness were abrogated, and most of the IL-1beta/TNF-alpha-mediated changes in proinflammatory gene expression were repressed, although mRNA expression of a small number of genes was enhanced by DEX. Collectively, the observations support the concept that, together with its role as a regulator of airway tone, in response to IL-1beta/TNF-alpha, the ASM expresses a host of glucocorticoid-sensitive genes that contribute to the altered structure and function of the airways in the pro-asthmatic state. We speculate that glucocorticoid-sensitive, cytokine-induced pathways involved in ASM cell signaling represent important targets for new therapeutic interventions.

A large Icelandic family with early osteoarthritis of the hip associated with a susceptibility locus on chromosome 16p.

OBJECTIVE: To describe a large kinship with inherited hip osteoarthritis (OA) and its associated susceptibility locus. METHODS: Four generations of a kinship with familial hip OA were identified and characterized by family history and by clinical, radiographic, and histopathologic examination. In the genome-wide search for a susceptibility locus, OA cases were defined as those who had undergone total hip replacement associated with a clinical and radiographic diagnosis of hip OA. A genome-wide scan was performed using a framework set of microsatellite markers with an average spacing of 10 cM. RESULTS: The hip OA of this family was indistinguishable from that of idiopathic, nonfamilial hip OA. There was no apparent evidence of spondyloepiphyseal dysplasia or other dysplasias usually associated with mutations in collagen genes. The genome-wide scan revealed a locus on chromosome 16p between 28 cM and 47 cM from the telomere, and this locus met the criteria for suggestive linkage (multipoint allele-sharing logarithm of odds [LOD] score 2.58, P = 1.6 x 10(-4)). Two additional regions with LOD scores of >1.5 were obtained. CONCLUSION: We have identified and described the largest kinship with familial hip OA reported to date. Evidence for linkage in this family suggests that a gene for susceptibility to hip OA exists on chromosome 16p. This represents an independent identification of a susceptibility locus previously reported for hip OA in this geographic region.

The genealogic approach to human genetics of disease.

The goal of modern human genetics is to correlate genes with disease or, more specifically, relate genetic variation to phenotypic variation. Although this correlation is usually straightforward in the Mendelian disorders, it has proved to be much more difficult to find in the common diseases because they appear to be more complex, likely involving an interplay among multiple genes and between genes and the environment. Although the strategy of linkage mapping of families was very successful when it was applied to the rare monogenic diseases, few common diseases have been mapped to statistical significance. Many investigators are now abandoning linkage analysis altogether and are moving to a candidate gene case-control strategy. In this article, we describe a genealogic approach to mapping human disease genes and provide three examples of how we have used it to map common diseases to statistical significance. We focus on a simple population with little historic migration and use a computerized genealogy database to increase the number of patients who can be compared with other affected relatives through high-density microsatellite genotyping. The genealogy helps determine which phenotypic classification is inherited and therefore possible to map. It may represent a more efficient strategy than candidate gene case-control studies for determination of what alleles or haplotypes are shared by patients in a population. We suggest that the genetics community not give up on linkage analysis, nor should it assume that the common diseases are too complex to map.

The OMgp gene, a second growth suppressor within the NF1 gene.

The Oligodendrocyte-Myelin glycoprotein gene (OMgp) is placed within an intron of the NF1 gene. Neurofibromin, the product of NF1, acts as a RasGAP and suppresses growth; inactivating mutations in NF1 lead to neurofibromatosis type 1. We report that OMgp also has growth suppressive effects and downregulates mitogenic signaling pathways closely related to those influenced by neurofibromin. Overexpression of OMgp alters mitogenic signaling in NIH3T3 fibroblasts. Cells overexpressing OMgp grow more slowly in serum compared to controls and show a partial G1 block upon cell cycle analysis. PDGF is the primary mitogen for fibroblasts in serum. Overexpression of OMgp alters PDGF signaling in fibroblasts which results in a block of mitogenic signaling. PDGF induced activation of c-Src is blocked, as is the induction of c-Myc and c-Fos, while tyrosine phosphorylation of the PDGFbeta receptor, PLCgamma1 and induction of c-Jun are intact. Although a number of genes embedded within other genes have been described, the biological significance of this arrangement remains unknown. We demonstrate here that structurally unrelated products of two such genes may exercise closely related functions. Our data also raise the possibility of a role for OMgp in disorders of cell proliferation such as NF1.

Expression of the oligodendrocyte-myelin glycoprotein by neurons in the mouse central nervous system.

The oligodendrocyte-myelin glycoprotein (OMgp) is a 110-kDa glycosylphosphatidylinositol-linked protein that was initially identified as a myelin-specific protein but whose precise function remains unknown. In this study, immunohistochemistry, western blots, in situ hybridization, and northern blots were used to determine the distribution of OMgp in the mouse brain. OMgp is present in a concentration detectable on western blots in the brains of newborn mice, and its concentration gradually increases until day 24 of life. OMgp mRNA is also present in amounts detectable on northern blots in the brains of newborn mice, and its concentration gradually increases until day 21 of life, after which the concentration diminishes a little. Most of the OMgp in the mouse brain appears to be expressed in diverse groups of neurons, but it is particularly prominent in large projection neurons such as the pyramidal cells of the hippocampus, the Purkinje cells of the cerebellum, motoneurons in the brainstem, and anterior horn cells of the spinal cord. However, OMgp is not confined to these cells and is expressed in cells in the white matter as well. The OMgp gene is placed within an intron of the neurofibromatosis type I gene and on the opposite strand. This organization raises the possibility that there may be a relationship between the functions of the products of the two genes. In support of this possibility, we show that within the mouse CNS OMgp and neurofibromin are expressed in the same cell types.

The epidermal growth factor receptor associates with and recruits phosphatidylinositol 3-kinase to the platelet-derived growth factor beta receptor.

Receptor tyrosine kinases are classified into subfamilies, which are believed to function independently, with heterodimerization occurring only within the same subfamily. In this study, we present evidence suggesting a direct interaction between the epidermal growth factor (EGF) receptor (EGFR) and the platelet-derived growth factor beta (PDGFbeta) receptor (PDGFbetaR), members of different receptor tyrosine kinase subfamilies. We find that the addition of EGF to COS-7 cells and to human foreskin Hs27 fibroblasts results in a rapid tyrosine phosphorylation of the PDGFbetaR and results in the recruitment of phosphatidylinositol 3-kinase to the PDGFbetaR. In R1hER cells, which overexpress the EGFR, we find ligand-independent tyrosine phosphorylation of the PDGFbetaR and the constitutive binding of a substantial amount of PI-3 kinase activity to it, mimicking the effect of ligand in untransfected cells. In support of the possibility that this may be a direct interaction, we show that the two receptors can be coimmunoprecipitated from untransfected Hs27 fibroblasts and from COS-7 cells. This association can be reconstituted by introducing the two receptors into 293 EBNA cells. The EGFR/PDGFbetaR association is ligand-independent in all cell lines tested. We also demonstrate that the fraction of PDGFbetaR bound to the EGFR in R1hER cells undergoes an EGF-induced mobility shift on Western blots indicative of phosphorylation. Our findings indicate that direct interactions between receptor tyrosine kinases classified under different subfamilies may be more widespread than previously believed.

Beta-trace gene expression is regulated by a core promoter and a distal thyroid hormone response element.

We isolated and characterized the human beta-Trace protein (betaTP) gene promoter. betaTP, also known as prostaglandin D2 synthase, is a lipocalin secreted from the choroid plexus and meninges into cerebrospinal fluid. Basal transcription of the betaTP gene is directed from a core promoter found within the first 325 bases of the 5'-flanking sequence. The betaTP gene promoter is responsive to thyroid hormone (3,3',5-triiodothyronine, T3) and efficiently repressed by unliganded human thyroid hormone receptor beta (TRbeta). Functional analysis of the betaTP promoter in TE671 cells revealed that responsiveness to T3 occurs in sequences 2.5 kilobase pairs 5' of the start site. Within the hormone-responsive region we identified a thyroid hormone response element (TRE) located from -2576 to -2562 base pairs relative to the transcription start site. The betaTP TRE is composed of two directly repeated consensus half-sites separated by a 3-base pair space (DR3). The betaTP TRE forms specific complexes with TRbeta. We have shown that a gene active in the choroid plexus and meninges is responsive to T3. T3 may play a role in the regulated transport of substances into the cerebrospinal fluid and ultimately the brain.

Interferon-gamma-induced oligodendrocyte cell death: implications for the pathogenesis of multiple sclerosis.

BACKGROUND: The histopathology of multiple sclerosis (MS) is characterized by a loss of myelin and oligodendrocytes, relative preservation of axons, and a modest inflammatory response. The reasons for this selective oligodendrocyte death and demyelination are unknown. MATERIALS AND METHODS: In light of the T lymphocyte and macrophage infiltrates in MS lesions and the numerous cytokines these cells secrete, the direct influence of cytokines on survival of cultured oligodendrocytes and sensory neurons was investigated. Expression of cytokines in vivo was determined by immunolabeling cryostat sections of snap-frozen tissue containing chronic active lesions from four different patients. The samples were also analyzed for the presence of apoptotic nuclei by in situ labeling of 3'-OH ends of degraded nuclear DNA. RESULTS: The results showed: (i) interferon-gamma (IFN gamma) to be a potent inducer of apoptosis among oligodendrocytes in vitro and that this effect can be reversed by leukemia inhibitory factor (LIF); (ii) IFN gamma has a minimal effect on the survival of cultured neurons; (iii) IFN gamma at the margins of active MS plaques but not in unaffected white matter; (iv) evidence for apoptosis of oligodendrocytes at the advancing margins of chronic active MS plaques. CONCLUSIONS: Injury to a substantial number of oligodendrocytes in MS is the results of programmed cell death rather than necrotic cell death mechanisms. We postulate that IFN gamma plays a role in the pathogenesis of MS by activating apoptosis in oligodendrocytes.

Focal brain injury and upregulation of a developmentally regulated extracellular matrix protein.

Tenascin is an extracellular matrix glycoprotein expressed during both normal development and neoplastic growth in both neural and nonneural tissues. During development of the central nervous system (CNS), tenascin is synthesized by glial cells, in particular by immature astrocytes, and is concentrated in transient boundaries around emerging groups of functionally distinct neurons. In the mature CNS, only low levels of the glycoprotein can be detected. The present study demonstrates that following trauma to the adult human cerebral cortex, discrete populations of reactive astrocytes upregulate their expression of tenascin and dramatically increase their transcription of the tenascin gene. The enhanced expression of tenascin may be involved in CNS wound healing, and may also affect neurite growth within and around a brain lesion.

Inhibition of T cell activation by the extracellular matrix protein tenascin.

Tenascin (TN) is an extracellular matrix protein that is expressed widely in the fetus and sparingly in the adult, but reappears at high levels in certain areas of tissue insult such as tumor matrices and sites of wound healing. We show here that soluble TN inhibits proliferation of human T cells in response to alpha CD3 Ab co-immobilized with the extracellular matrix protein fibronectin (FN). TN also inhibits proliferation driven by alpha CD3/IL-2 or by phorbol ester/IL-2, and it prevents high level induction of IL-2R. The presence of TN in culture medium does not detectably alter the pattern of tyrosine phosphorylation resulting from T cell triggering with alpha CD3, but at later time points prevents the appearance of functional NF-AT1 transcription factor complexes in T cell nuclear extracts. These findings are consistent with the postulated role for TN as a natural antagonist to FN action, and suggest that T cell responses occurring at tissue sites in which TN is expressed could be influenced by its presence.