Dissection of artifactual and confounding glial signatures by single-cell sequencing of mouse and human brain

A key aspect of nearly all single-cell sequencing experiments is dissociation of intact tissues into single-cell suspensions. While many protocols have been optimized for optimal cell yield, they have often overlooked the effects that dissociation can have on ex vivo gene expression. Here, we demonstrate that use of enzymatic dissociation on brain tissue induces an aberrant ex vivo gene expression signature, most prominently in microglia, which is prevalent in published literature and can substantially confound downstream analyses. To address this issue, we present a rigorously validated protocol that preserves both in vivo transcriptional profiles and cell-type diversity and yield across tissue types and species. We also identify a similar signature in postmortem human brain single-nucleus RNA-sequencing datasets, and show that this signature is induced in freshly isolated human tissue by exposure to elevated temperatures ex vivo. Together, our results provide a methodological solution for preventing artifactual gene expression changes during fresh tissue digestion and a reference for future deeper analysis of the potential confounding states present in postmortem human samples

A Selection Index for Gene Expression Evolution and Its Application to the Divergence between Humans and Chimpanzees

The importance of gene regulation in animal evolution is a matter of long-standing interest, but measuring the impact of selection on gene expression has proven a challenge. Here, we propose a selection index of gene expression as a straightforward method for assessing the mode and strength of selection operating on gene expression levels. The index is based on the widely used McDonald-Kreitman test and requires the estimation of four quantities: the within-species and between-species expression variances as well as the sequence heterozygosity and divergence of neutrally evolving sequences. We apply the method to data from human and chimpanzee lymphoblastoid cell lines and show that gene expression is in general under strong stabilizing selection. We also demonstrate how the same framework can be used to estimate the proportion of adaptive gene expression evolution

A map of transcriptional heterogeneity and regulatory variation in human microglia.

Microglia, the tissue-resident macrophages of the central nervous system (CNS), play critical roles in immune defense, development and homeostasis. However, isolating microglia from humans in large numbers is challenging. Here, we profiled gene expression variation in primary human microglia isolated from 141 patients undergoing neurosurgery. Using single-cell and bulk RNA sequencing, we identify how age, sex and clinical pathology influence microglia gene expression and which genetic variants have microglia-specific functions using expression quantitative trait loci (eQTL) mapping. We follow up one of our findings using a human induced pluripotent stem cell-based macrophage model to fine-map a candidate causal variant for Alzheimer's disease at the BIN1 locus. Our study provides a population-scale transcriptional map of a critically important cell for human CNS development and disease

Common genetic variation drives molecular heterogeneity in human iPSCs.

Technology utilizing human induced pluripotent stem cells (iPS cells) has enormous potential to provide improved cellular models of human disease. However, variable genetic and phenotypic characterization of many existing iPS cell lines limits their potential use for research and therapy. Here we describe the systematic generation, genotyping and phenotyping of 711 iPS cell lines derived from 301 healthy individuals by the Human Induced Pluripotent Stem Cells Initiative. Our study outlines the major sources of genetic and phenotypic variation in iPS cells and establishes their suitability as models of complex human traits and cancer. Through genome-wide profiling we find that 5-46% of the variation in different iPS cell phenotypes, including differentiation capacity and cellular morphology, arises from differences between individuals. Additionally, we assess the phenotypic consequences of genomic copy-number alterations that are repeatedly observed in iPS cells. In addition, we present a comprehensive map of common regulatory variants affecting the transcriptome of human pluripotent cells

Association of genetic variants in complement factor H and factor H-related genes with systemic lupus erythematosus susceptibility

Systemic lupus erythematosus (SLE), a complex polygenic autoimmune disease, is associated with increased complement activation. Variants of genes encoding complement regulator factor H (CFH) and five CFH-related proteins (CFHR1-CFHR5) within the chromosome 1q32 locus linked to SLE, have been associated with multiple human diseases and may contribute to dysregulated complement activation predisposing to SLE. We assessed 60 SNPs covering the CFH-CFHRs region for association with SLE in 15,864 case-control subjects derived from four ethnic groups. Significant allelic associations with SLE were detected in European Americans (EA) and African Americans (AA), which could be attributed to an intronic CFH SNP (rs6677604, in intron 11, Pmeta = 6.6×10-8, OR = 1.18) and an intergenic SNP between CFHR1 and CFHR4 (rs16840639, Pmeta = 2.9×10-7, OR = 1.17) rather than to previously identified disease-associated CFH exonic SNPs, including I62V, Y402H, A474A, and D936E. In addition, allelic association of rs6677604 with SLE was subsequently confirmed in Asians (AS). Haplotype analysis revealed that the underlying causal variant, tagged by rs6677604 and rs16840639, was localized to a ~146 kb block extending from intron 9 of CFH to downstream of CFHR1. Within this block, the deletion of CFHR3 and CFHR1 (CFHR3-1Δ), a likely causal variant measured using multiplex ligation-dependent probe amplification, was tagged by rs6677604 in EA and AS and rs16840639 in AA, respectively. Deduced from genotypic associations of tag SNPs in EA, AA, and AS, homozygous deletion of CFHR3-1Δ (Pmeta = 3.2×10-7, OR = 1.47) conferred a higher risk of SLE than heterozygous deletion (Pmeta = 3.5×10-4, OR = 1.14). These results suggested that the CFHR3-1Δ deletion within the SLE-associated block, but not the previously described exonic SNPs of CFH, might contribute to the development of SLE in EA, AA, and AS, providing new insights into the role of complement regulators in the pathogenesis of SLE

Variation in the ICAM1-ICAM4-ICAM5 locus is associated with systemic lupus erythematosus susceptibility in multiple ancestries

Objective: Systemic lupus erythematosus (SLE; OMIM 152700) is a chronic autoimmune disease for which the aetiology includes genetic and environmental factors. ITGAM, integrin ?M(complement component 3 receptor 3 subunit) encoding a ligand for intracellular adhesion molecule (ICAM) proteins, is an established SLE susceptibility locus. This study aimed to evaluate the independent and joint effects of genetic variations in the genes that encode ITGAM and ICAM. Methods: The authors examined several markers in the ICAM1-ICAM4-ICAM5 locus on chromosome 19p13 and the single ITGAM polymorphism (rs1143679) using a large-scale case-control study of 17 481 unrelated participants from four ancestry populations. The singlemarker association and gene-gene interaction were analysed for each ancestry, and a meta-analysis across the four ancestries was performed. Results: The A-allele of ICAM1-ICAM4-ICAM5 rs3093030, associated with elevated plasma levels of soluble ICAM1, and the A-allele of ITGAM rs1143679 showed the strongest association with increased SLE susceptibility in each of the ancestry populations and the trans-ancestry meta-analysis (ORmeta=1.16, 95% CI 1.11 to 1.22; p=4.88 × 10-10 and ORmeta=1.67, 95% CI 1.55 to 1.79; p=3.32 × 10-46, respectively). The effect of the ICAM single-nucleotide polymorphisms (SNPs) was independent of the effect of the ITGAM SNP rs1143679, and carriers of both ICAM rs3093030-AA and ITGAM rs1143679-AA had an OR of 4.08 compared with those with no risk allele in either SNP (95% CI 2.09 to 7.98; p=3.91 × 10-5). Conclusion: These findings are the first to suggest that an ICAM-integrin-mediated pathway contributes to susceptibility to SLE

Genome-wide association study identifies Sjögren’s risk loci with functional implications in immune and glandular cells

Sjögren’s disease is a complex autoimmune disease with twelve established susceptibility loci. This genome-wide association study (GWAS) identifies ten novel genome-wide significant (GWS) regions in Sjögren’s cases of European ancestry: CD247, NAB1, PTTG1-MIR146A, PRDM1-ATG5, TNFAIP3, XKR6, MAPT-CRHR1, RPTOR-CHMP6-BAIAP6, TYK2, SYNGR1. Polygenic risk scores yield predictability (AUROC = 0.71) and relative risk of 12.08. Interrogation of bioinformatics databases refine the associations, define local regulatory networks of GWS SNPs from the 95% credible set, and expand the implicated gene list to >40. Many GWS SNPs are eQTLs for genes within topologically associated domains in immune cells and/or eQTLs in the main target tissue, salivary glands.Research reported in this publication was supported by the National Institutes of Health (NIH): R01AR073855 (C.J.L.), R01AR065953 (C.J.L.), R01AR074310 (A.D.F.), P50AR060804 (K.L.S.), R01AR050782 (K.L.S), R01DE018209 (K.L.S.), R33AR076803 (I.A.), R21AR079089 (I.A.); NIDCR Sjögren’s Syndrome Clinic and Salivary Disorders Unit were supported by NIDCR Division of Intramural Research at the National Institutes of Health funds - Z01-DE000704 (B.W.); Birmingham NIHR Biomedical Research Centre (S.J.B.); Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC 2155 – Projektnummer 390874280 (T.W.); Research Council of Norway (Oslo, Norway) – Grant 240421 (TR.R.), 316120 (M.W-H.); Western Norway Regional Health Authority (Helse Vest) – 911807, 912043 (R.O.); Swedish Research Council for Medicine and Health (L.R., G.N., M.W-H.); Swedish Rheumatism Association (L.R., G.N., M.W-H.); King Gustav V’s 80-year Foundation (G.N.); Swedish Society of Medicine (L.R., G.N., M.W-H.); Swedish Cancer Society (E.B.); Sjögren’s Syndrome Foundation (K.L.S.); Phileona Foundation (K.L.S.). The Stockholm County Council (M.W-H.); The Swedish Twin Registry is managed through the Swedish Research Council - Grant 2017-000641. The French ASSESS (Atteinte Systémique et Evolution des patients atteints de Syndrome de Sjögren primitive) was sponsored by Assistance Publique-Hôpitaux de Paris (Ministry of Health, PHRC 2006 P060228) and the French society of Rheumatology (X.M.).publishedVersio

Genome-wide association study identifies Sjögren's risk loci with functional implications in immune and glandular cells.

Sjögren’s disease is a complex autoimmune disease with twelve established susceptibility loci. This genome-wide association study (GWAS) identifies ten novel genome-wide significant (GWS) regions in Sjögren’s cases of European ancestry: CD247, NAB1, PTTG1-MIR146A, PRDM1-ATG5, TNFAIP3, XKR6, MAPT-CRHR1, RPTOR-CHMP6-BAIAP6, TYK2, SYNGR1. Polygenic risk scores yield predictability (AUROC = 0.71) and relative risk of 12.08. Interrogation of bioinformatics databases refine the associations, define local regulatory networks of GWS SNPs from the 95% credible set, and expand the implicated gene list to >40. Many GWS SNPs are eQTLs for genes within topologically associated domains in immune cells and/or eQTLs in the main target tissue, salivary glands.We thank all the research and clinical staff, consortium investigators, and study participants (detailed in Supplementary Information), and funding agencies who made this study possible. The content of this publication is solely the responsibility of the authors and does not represent the official views of the funding agencies listed below. Research reported in this publication was supported by the National Institutes of Health (NIH): R01AR073855 (C.J.L.), R01AR065953 (C.J.L.), R01AR074310 (A.D.F.), P50AR060804 (K.L.S.), R01AR050782 (K.L.S), R01DE018209 (K.L.S.), R33AR076803 (I.A.), R21AR079089 (I.A.); NIDCR Sjögren’s Syndrome Clinic and Salivary Disorders Unit were supported by NIDCR Division of Intramural Research at the National Institutes of Health funds - Z01-DE000704 (B.W.); Birmingham NIHR Biomedical Research Centre (S.J.B.); Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC 2155 – Projektnummer 390874280 (T.W.); Research Council of Norway (Oslo, Norway) – Grant 240421 (TR.R.), 316120 (M.W-H.); Western Norway Regional Health Authority (Helse Vest) – 911807, 912043 (R.O.); Swedish Research Council for Medicine and Health (L.R., G.N., M.W-H.); Swedish Rheumatism Association (L.R., G.N., M.W-H.); King Gustav V’s 80-year Foundation (G.N.); Swedish Society of Medicine (L.R., G.N., M.W-H.); Swedish Cancer Society (E.B.); Sjögren’s Syndrome Foundation (K.L.S.); Phileona Foundation (K.L.S.). The Stockholm County Council (M.W-H.); FOREUM Foundation for Research in Rheumatology (R.J., M.W-H). The Swedish Twin Registry is managed through the Swedish Research Council - Grant 2017-000641. The French ASSESS (Atteinte Systémique et Evolution des patients atteints de Syndrome de Sjögren primitive) was sponsored by Assistance Publique-Hôpitaux de Paris (Ministry of Health, PHRC 2006 P060228) and the French society of Rheumatology (X.M.). We want to acknowledge the following invesigators who recruited patients: Jacques-Eric Gottenberg, Valerie Devauchelle-Pensec, Jean Jacques Dubost, Anne-Laure Fauchais, Vincent Goeb, Eric Hachulla, Claire Larroche, Véronique Le Guern, Jacques Morel, Aleth Perdriger, Emmanuelle Dernis, Stéphanie Rist, Damien Sene, Olivier Vittecoq. We also thank Sarah Tubiana and all staff members of the Bichat Hospital Biological Resource Center (Paris) for centralizing and managing biological collection. We also thank Rezvan Kiani Dehkordi, Karolina Tandre, Käth Nilsson, Marianne Eidsheim, Kjerstin Jacobsen, Ingeborg Kvivik and Kjetil Bårdsen for collecting patient blood samples. We acknowledge the SNP&SEQ Technology Platform, Uppsala, part of National Genomics Infrastructure (NGI) Sweden, for genotyping of Scandinavian samples, and the Swedish Twin Registry for access to data. The SNP&SEQ Technology Platform was supported by Science for Life Laboratory, Uppsala University, the Knut and Alice Wallenberg Foundation and the Swedish Research Council. Last, we thank the investigators for the following dbGaP studies: Phs000428.v2.p2: This study used control data from the Health and Retirement Study in dbGaP (phs000428.v2.p2) submitted by David Weir, PhD at the University of Michigan and funded by the National Institute of Aging RC2 AG036495 and RC4 AG039029. Phs000672.v1.p1: Genotype data from the Sjögren’s International Collaborative Clinical Alliance (SICCA) Registry was obtained through dbGAP accession number phs000672.v1.p1. This study was supported by the National Institute of Dental and Craniofacial Research (NIDCR), the National Eye Institute, and the Office of Research on Women’s Health through contract number N01-DE-32636. Genotyping services were provided by the Center for Inherited Disease Research (CIDR). CIDR is fully funded through a federal contract from the National Institutes of Health (NIH) to the Johns Hopkins University (contract numbers HHSN268200782096C, HHSN268201100011I, HHSN268201200008I). Funds for genotyping were provided by the NIDCR through CIDR’s NIH contract. Assistance with data cleaning and imputation was provided by the University of Washington. We thank investigators from the following studies that provided DNA samples for genotyping: the Genetic Architecture of Smoking and Smoking Cessation, Collaborative Genetic Study of Nicotine Dependence (phs000404.v1.p1); Age-Related Eye Disease Study (AREDS) - Genetic Variation in Refractive Error Substudy (phs000429.v1.p1); and National Institute of Mental Health’s Human Genetics Initiative (phs000021.v3.p2, phs000167.v1.p1). We thank the many clinical collaborators and research participants who contributed to this research. Phs000196.v3.p1: Investigators and Parkinson Disease patients that contributed to this Genome-wide Association Study of Parkinson Disease. phs000187.v1.p1: Research support to collect data and develop an application to support the High Density SNP Association Analysis of Melanoma project was provided by 3P50CA093459, 5P50CA097007, 5R01ES011740, and 5R01CA133996

Variants at multiple loci implicated in both innate and adaptive immune responses are associated with Sjögren’s syndrome

Sjögren’s syndrome is a common autoimmune disease (~0.7% of European Americans) typically presenting as keratoconjunctivitis sicca and xerostomia. In addition to strong association within the HLA region at 6p21 (Pmeta=7.65×10−114), we establish associations with IRF5-TNPO3 (Pmeta=2.73×10−19), STAT4 (Pmeta=6.80×10−15), IL12A (Pmeta =1.17×10−10), FAM167A-BLK (Pmeta=4.97×10−10), DDX6-CXCR5 (Pmeta=1.10×10−8), and TNIP1 (Pmeta=3.30×10−8). Suggestive associations with Pmeta<5×10−5 were observed with 29 regions including TNFAIP3, PTTG1, PRDM1, DGKQ, FCGR2A, IRAK1BP1, ITSN2, and PHIP amongst others. These results highlight the importance of genes involved in both innate and adaptive immunity in Sjögren’s syndrome

Petabytes to science

Paper published on ArXiv to raise awareness and start discussions about data access in the era of large astronomical surveys.A Kavli foundation sponsored workshop on the theme Petabytes to Science was held 12th to 14th of February 2019 in Las Vegas. The aim of the this workshop was to discuss important trends and technologies which may support astronomy. We also tackled how to better shape the workforce for the new trends and how we should approach education and public outreach. This document was coauthored during the workshop and edited in the weeks after. It comprises the discussions and highlights many recommendations which came out of the workshop. We shall distill parts of this document and formulate potential white papers for the decadal survey.Preprin