Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use.

Tobacco and alcohol use are leading causes of mortality that influence risk for many complex diseases and disorders1. They are heritable2,3 and etiologically related4,5 behaviors that have been resistant to gene discovery efforts6-11. In sample sizes up to 1.2 million individuals, we discovered 566 genetic variants in 406 loci associated with multiple stages of tobacco use (initiation, cessation, and heaviness) as well as alcohol use, with 150 loci evidencing pleiotropic association. Smoking phenotypes were positively genetically correlated with many health conditions, whereas alcohol use was negatively correlated with these conditions, such that increased genetic risk for alcohol use is associated with lower disease risk. We report evidence for the involvement of many systems in tobacco and alcohol use, including genes involved in nicotinic, dopaminergic, and glutamatergic neurotransmission. The results provide a solid starting point to evaluate the effects of these loci in model organisms and more precise substance use measures

Autoimmunity-Associated LYP-W620 Does Not Impair Thymic Negative Selection of Autoreactive T Cells.

A C1858T (R620W) variation in the PTPN22 gene encoding the tyrosine phosphatase LYP is a major risk factor for human autoimmunity. LYP is a known negative regulator of signaling through the T cell receptor (TCR), and murine Ptpn22 plays a role in thymic selection. However, the mechanism of action of the R620W variant in autoimmunity remains unclear. One model holds that LYP-W620 is a gain-of-function phosphatase that causes alterations in thymic negative selection and/or thymic output of regulatory T cells (Treg) through inhibition of thymic TCR signaling. To test this model, we generated mice in which the human LYP-W620 variant or its phosphatase-inactive mutant are expressed in developing thymocytes under control of the proximal Lck promoter. We found that LYP-W620 expression results in diminished thymocyte TCR signaling, thus modeling a "gain-of-function" of LYP at the signaling level. However, LYP-W620 transgenic mice display no alterations of thymic negative selection and no anomalies in thymic output of CD4(+)Foxp3(+) Treg were detected in these mice. Lck promoter-directed expression of the human transgene also causes no alteration in thymic repertoire or increase in disease severity in a model of rheumatoid arthritis, which depends on skewed thymic selection of CD4(+) T cells. Our data suggest that a gain-of-function of LYP is unlikely to increase risk of autoimmunity through alterations of thymic selection and that LYP likely acts in the periphery perhaps selectively in regulatory T cells or in another cell type to increase risk of autoimmunity

Genetic diversity fuels gene discovery for tobacco and alcohol use

Tobacco and alcohol use are heritable behaviours associated with 15% and 5.3% of worldwide deaths, respectively, due largely to broad increased risk for disease and injury(1-4). These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries(5). Here we leveraged global genetic diversity across 3.4 million individuals from four major clines of global ancestry (approximately 21% non-European) to power the discovery and fine-mapping of genomic loci associated with tobacco and alcohol use, to inform function of these loci via ancestry-aware transcriptome-wide association studies, and to evaluate the genetic architecture and predictive power of polygenic risk within and across populations. We found that increases in sample size and genetic diversity improved locus identification and fine-mapping resolution, and that a large majority of the 3,823 associated variants (from 2,143 loci) showed consistent effect sizes across ancestry dimensions. However, polygenic risk scores developed in one ancestry performed poorly in others, highlighting the continued need to increase sample sizes of diverse ancestries to realize any potential benefit of polygenic prediction.Peer reviewe

Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use

Tobacco and alcohol use are leading causes of mortality that influence risk for many complex diseases and disorders 1 . They are heritable 2,3 and etiologically related 4,5 behaviors that have been resistant to gene discovery efforts 6–11 . In sample sizes up to 1.2 million individuals, we discovered 566 genetic variants in 406 loci associated with multiple stages of tobacco use (initiation, cessation, and heaviness) as well as alcohol use, with 150 loci evidencing pleiotropic association. Smoking phenotypes were positively genetically correlated with many health conditions, whereas alcohol use was negatively correlated with these conditions, such that increased genetic risk for alcohol use is associated with lower disease risk. We report evidence for the involvement of many systems in tobacco and alcohol use, including genes involved in nicotinic, dopaminergic, and glutamatergic neurotransmission. The results provide a solid starting point to evaluate the effects of these loci in model organisms and more precise substance use measures

Genome-wide analyses identify a role for SLC17A4 and AADAT in thyroid hormone regulation.

Thyroid dysfunction is an important public health problem, which affects 10% of the general population and increases the risk of cardiovascular morbidity and mortality. Many aspects of thyroid hormone regulation have only partly been elucidated, including its transport, metabolism, and genetic determinants. Here we report a large meta-analysis of genome-wide association studies for thyroid function and dysfunction, testing 8 million genetic variants in up to 72,167 individuals. One-hundred-and-nine independent genetic variants are associated with these traits. A genetic risk score, calculated to assess their combined effects on clinical end points, shows significant associations with increased risk of both overt (Graves' disease) and subclinical thyroid disease, as well as clinical complications. By functional follow-up on selected signals, we identify a novel thyroid hormone transporter (SLC17A4) and a metabolizing enzyme (AADAT). Together, these results provide new knowledge about thyroid hormone physiology and disease, opening new possibilities for therapeutic targets

AI is a viable alternative to high throughput screening: a 318-target study

: High throughput screening (HTS) is routinely used to identify bioactive small molecules. This requires physical compounds, which limits coverage of accessible chemical space. Computational approaches combined with vast on-demand chemical libraries can access far greater chemical space, provided that the predictive accuracy is sufficient to identify useful molecules. Through the largest and most diverse virtual HTS campaign reported to date, comprising 318 individual projects, we demonstrate that our AtomNet® convolutional neural network successfully finds novel hits across every major therapeutic area and protein class. We address historical limitations of computational screening by demonstrating success for target proteins without known binders, high-quality X-ray crystal structures, or manual cherry-picking of compounds. We show that the molecules selected by the AtomNet® model are novel drug-like scaffolds rather than minor modifications to known bioactive compounds. Our empirical results suggest that computational methods can substantially replace HTS as the first step of small-molecule drug discovery

Overexpression of the Cytokine BAFF and Autoimmunity Risk

$\textbf{BACKGROUND}$: Genomewide association studies of autoimmune diseases have mapped hundreds of susceptibility regions in the genome. However, only for a few association signals has the causal gene been identified, and for even fewer have the causal variant and underlying mechanism been defined. Coincident associations of DNA variants affecting both the risk of autoimmune disease and quantitative immune variables provide an informative route to explore disease mechanisms and drug-targetable pathways. $\textbf{METHODS}$: Using case-control samples from Sardinia, Italy, we performed a genomewide association study in multiple sclerosis followed by TNFSF13B locus-specific association testing in systemic lupus erythematosus (SLE). Extensive phenotyping of quantitative immune variables, sequence-based fine mapping, cross-population and cross-phenotype analyses, and gene-expression studies were used to identify the causal variant and elucidate its mechanism of action. Signatures of positive selection were also investigated. $\textbf{RESULTS}$: A variant in TNFSF13B, encoding the cytokine and drug target B-cell activating factor (BAFF), was associated with multiple sclerosis as well as SLE. The disease-risk allele was also associated with up-regulated humoral immunity through increased levels of soluble BAFF, B lymphocytes, and immunoglobulins. The causal variant was identified: an insertion-deletion variant, GCTGT→A (in which A is the risk allele), yielded a shorter transcript that escaped microRNA inhibition and increased production of soluble BAFF, which in turn up-regulated humoral immunity. Population genetic signatures indicated that this autoimmunity variant has been evolutionarily advantageous, most likely by augmenting resistance to malaria. $\textbf{CONCLUSIONS}$: A TNFSF13B variant was associated with multiple sclerosis and SLE, and its effects were clarified at the population, cellular, and molecular levels. (Funded by the Italian Foundation for Multiple Sclerosis and others.).Supported by grants (2011/R/13 and 2015/R/09, to Dr. Cucca) from the Italian Foundation for Multiple Sclerosis; contracts (N01-AG-1-2109 and HHSN271201100005C, to Dr. Cucca) from the Intramural Research Program of the National Institute on Aging, National Institutes of Health (NIH); a grant (FaReBio2011 “Farmaci e Reti Biotecnologiche di Qualità,” to Dr. Cucca) from the Italian Ministry of Economy and Finance; a grant (633964, to Dr. Cucca) from the Horizon 2020 Research and Innovation Program of the European Union; a grant (U1301.2015/AI.1157.BE Prat. 2015-1651, to Dr. Cucca) from Fondazione di Sardegna; grants (“Centro per la ricerca di nuovi farmaci per malattie rare, trascurate e della povertà” and “Progetto collezione di composti chimici ed attività di screening,” to Dr. Cucca) from Ministero dell’Istruzione, dell’Università e della Ricerca; grants (HG005581, HG005552, HG006513, and HG007022, to Dr. Abecasis) from the National Human Genome Research Institute; a grant (9-2011-253, to Dr. Todd) from JDRF; a grant (091157, to Dr. Todd) from the Wellcome Trust; a grant (to Dr. Todd) from the National Institute for Health Research (NIHR); and the NIHR Cambridge Biomedical Research Centre. Dr. Idda was a recipient of a Master and Back fellowship from the Autonomous Region of Sardinia

Population- and individual-specific regulatory variation in Sardinia

Genetic studies of complex traits have mainly identified associations with noncoding variants. To further determine the contribution of regulatory variation, we combined whole-genome and transcriptome data for 624 individuals from Sardinia to identify common and rare variants that influence gene expression and splicing. We identified 21,183 expression quantitative trait loci (eQTLs) and 6,768 splicing quantitative trait loci (sQTLs), including 619 new QTLs. We identified high-frequency QTLs and found evidence of selection near genes involved in malarial resistance and increased multiple sclerosis risk, reflecting the epidemiological history of Sardinia. Using family relationships, we identified 809 segregating expression outliers (median z score of 2.97), averaging 13.3 genes per individual. Outlier genes were enriched for proximal rare variants, providing a new approach to study large-effect regulatory variants and their relevance to traits. Our results provide insight into the effects of regulatory variants and their relationship to population history and individual genetic risk.M.P. is supported by the European Union’s Horizon 2020 Research and Innovation Programme under grant agreement 633964 (ImmunoAgeing). Z.Z. is supported by the National Science Foundation (NSF) GRFP (DGE- 114747) and by the Stanford Center for Computational, Evolutionary, and Human Genomics (CEHG). Z.Z., J.R.D., and G.T.H. also acknowledge support from the Stanford Genome Training Program (SGTP; NIH/NHGRI T32HG000044). J.R.D. is supported by the Stanford Graduate Fellowship. K.R.K. is supported by Department of Defense, Air Force Office of Scientific Research, National Defense Science and Engineering Graduate (NDSEQ) Fellowship 32 CFR 168a. S.J.S. is supported by the NIHR Cambridge Biomedical Research Centre. The SardiNIA project is supported in part by the intramural program of the National Institute on Aging through contract HHSN271201100005C to the Consiglio Nazionale delle Ricerche of Italy. The RNA sequencing was supported by the PB05 InterOmics MIUR Flagship grant; by the FaReBio2011 “Farmaci e Reti Biotecnologiche di Qualità” grant; and by Sardinian Autonomous Region (L.R. no. 7/2009) grant cRP3-154 to F. Cucca, who is also supported by the Italian Foundation for Multiple Sclerosis (FISM 2015/R/09) and by the Fondazione di Sardegna (ex Fondazione Banco di Sardegna, Prot. U1301.2015/AI.1157.BE Prat. 2015-1651). S.B.M. is supported by the US National Institutes of Health through R01HG008150, R01MH101814, U01HG007436, and U01HG009080. All of the authors would like to thank the CRS4 and the SCGPM for the computational infrastructure supporting this project

Autoimmunity Is a Significant Feature of Idiopathic Pulmonary Arterial Hypertension.

Rationale: Autoimmunity is believed to play a role in idiopathic pulmonary arterial hypertension (IPAH). It is not clear whether this is causative or a bystander of disease and if it carries any prognostic or treatment significance. Objectives: To study autoimmunity in IPAH using a large cross-sectional cohort. Methods: Assessment of the circulating immune cell phenotype was undertaken using flow cytometry, and the profile of serum immunoglobulins was generated using a standardized multiplex array of 19 clinically validated autoantibodies in 473 cases and 946 control subjects. Additional glutathione S-transferase fusion array and ELISA data were used to identify a serum autoantibody to BMPR2 (bone morphogenetic protein receptor type 2). Clustering analyses and clinical correlations were used to determine associations between immunogenicity and clinical outcomes. Measurements and Main Results: Flow cytometric immune profiling demonstrates that IPAH is associated with an altered humoral immune response in addition to raised IgG3. Multiplexed autoantibodies were significantly raised in IPAH, and clustering demonstrated three distinct clusters: "high autoantibody," "low autoantibody," and a small "intermediate" cluster exhibiting high concentrations of ribonucleic protein complex. The high-autoantibody cluster had worse hemodynamics but improved survival. A small subset of patients demonstrated immunoglobulin reactivity to BMPR2. Conclusions: This study establishes aberrant immune regulation and presence of autoantibodies as key features in the profile of a significant proportion of patients with IPAH and is associated with clinical outcomes