51
|
Oliveros W, Delfosse K, Lato DF, Kiriakopulos K, Mokhtaridoost M, Said A, McMurray BJ, Browning JW, Mattioli K, Meng G, Ellis J, Mital S, Melé M, Maass PG. Systematic characterization of regulatory variants of blood pressure genes. CELL GENOMICS 2023; 3:100330. [PMID: 37492106 PMCID: PMC10363820 DOI: 10.1016/j.xgen.2023.100330] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/15/2022] [Revised: 03/29/2023] [Accepted: 04/28/2023] [Indexed: 07/27/2023]
Abstract
High blood pressure (BP) is the major risk factor for cardiovascular disease. Genome-wide association studies have identified genetic variants for BP, but functional insights into causality and related molecular mechanisms lag behind. We functionally characterize 4,608 genetic variants in linkage with 135 BP loci in vascular smooth muscle cells and cardiomyocytes by massively parallel reporter assays. High densities of regulatory variants at BP loci (i.e., ULK4, MAP4, CFDP1, PDE5A) indicate that multiple variants drive genetic association. Regulatory variants are enriched in repeats, alter cardiovascular-related transcription factor motifs, and spatially converge with genes controlling specific cardiovascular pathways. Using heuristic scoring, we define likely causal variants, and CRISPR prime editing finally determines causal variants for KCNK9, SFXN2, and PCGF6, which are candidates for developing high BP. Our systems-level approach provides a catalog of functionally relevant variants and their genomic architecture in two trait-relevant cell lines for a better understanding of BP gene regulation.
Collapse
Affiliation(s)
- Winona Oliveros
- Life Sciences Department, Barcelona Supercomputing Center, 08034 Barcelona, Catalonia, Spain
| | - Kate Delfosse
- Genetics & Genome Biology Program, The Hospital for Sick Children, Toronto, ON M5G 0A4, Canada
| | - Daniella F. Lato
- Genetics & Genome Biology Program, The Hospital for Sick Children, Toronto, ON M5G 0A4, Canada
| | - Katerina Kiriakopulos
- Genetics & Genome Biology Program, The Hospital for Sick Children, Toronto, ON M5G 0A4, Canada
- Department of Molecular Genetics, University of Toronto, Toronto, ON, Canada
| | - Milad Mokhtaridoost
- Genetics & Genome Biology Program, The Hospital for Sick Children, Toronto, ON M5G 0A4, Canada
| | - Abdelrahman Said
- Genetics & Genome Biology Program, The Hospital for Sick Children, Toronto, ON M5G 0A4, Canada
| | - Brandon J. McMurray
- Genetics & Genome Biology Program, The Hospital for Sick Children, Toronto, ON M5G 0A4, Canada
| | - Jared W.L. Browning
- Genetics & Genome Biology Program, The Hospital for Sick Children, Toronto, ON M5G 0A4, Canada
- Department of Molecular Genetics, University of Toronto, Toronto, ON, Canada
| | - Kaia Mattioli
- Division of Genetics, Department of Medicine, Brigham & Women’s Hospital and Harvard Medical School, Boston, MA 02115, USA
| | - Guoliang Meng
- Developmental and Stem Cell Biology Program, The Hospital for Sick Children, Toronto, ON M5G 0A4, Canada
| | - James Ellis
- Developmental and Stem Cell Biology Program, The Hospital for Sick Children, Toronto, ON M5G 0A4, Canada
- Department of Molecular Genetics, University of Toronto, Toronto, ON, Canada
| | - Seema Mital
- Genetics & Genome Biology Program, The Hospital for Sick Children, Toronto, ON M5G 0A4, Canada
- Ted Rogers Centre for Heart Research, Toronto, ON M5G 1X8, Canada
- Department of Pediatrics, The Hospital for Sick Children, University of Toronto, Toronto, ON M5G 0A4, Canada
| | - Marta Melé
- Life Sciences Department, Barcelona Supercomputing Center, 08034 Barcelona, Catalonia, Spain
| | - Philipp G. Maass
- Genetics & Genome Biology Program, The Hospital for Sick Children, Toronto, ON M5G 0A4, Canada
- Department of Molecular Genetics, University of Toronto, Toronto, ON, Canada
| |
Collapse
|
52
|
Viñas R, Joshi CK, Georgiev D, Lin P, Dumitrascu B, Gamazon ER, Liò P. Hypergraph factorization for multi-tissue gene expression imputation. NAT MACH INTELL 2023; 5:739-753. [PMID: 37771758 PMCID: PMC10538467 DOI: 10.1038/s42256-023-00684-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2022] [Accepted: 06/02/2023] [Indexed: 09/30/2023]
Abstract
Integrating gene expression across tissues and cell types is crucial for understanding the coordinated biological mechanisms that drive disease and characterise homeostasis. However, traditional multitissue integration methods cannot handle uncollected tissues or rely on genotype information, which is often unavailable and subject to privacy concerns. Here we present HYFA (Hypergraph Factorisation), a parameter-efficient graph representation learning approach for joint imputation of multi-tissue and cell-type gene expression. HYFA is genotype-agnostic, supports a variable number of collected tissues per individual, and imposes strong inductive biases to leverage the shared regulatory architecture of tissues and genes. In performance comparison on Genotype-Tissue Expression project data, HYFA achieves superior performance over existing methods, especially when multiple reference tissues are available. The HYFA-imputed dataset can be used to identify replicable regulatory genetic variations (eQTLs), with substantial gains over the original incomplete dataset. HYFA can accelerate the effective and scalable integration of tissue and cell-type transcriptome biorepositories.
Collapse
Affiliation(s)
- Ramon Viñas
- Department of Computer Science and Technology, University of Cambridge
| | | | - Dobrik Georgiev
- Department of Computer Science and Technology, University of Cambridge
| | - Phillip Lin
- Division of Genetic Medicine, Vanderbilt University Medical Center
| | - Bianca Dumitrascu
- Department of Statistics and Irving Institute for Cancer Dynamics, Columbia University
| | - Eric R. Gamazon
- Vanderbilt Genetics Institute and Data Science Institute, MRC Epidemiology Unit, University of Cambridge
| | - Pietro Liò
- Department of Computer Science and Technology, University of Cambridge
| |
Collapse
|
53
|
Rummel CK, Gagliardi M, Herholt A, Ahmad R, Murek V, Weigert L, Hausruckinger A, Maidl S, Jimenez-Barron L, Trastulla L, Eder M, Rossner M, Ziller MJ. Cell type and condition specific functional annotation of schizophrenia associated non-coding genetic variants. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.06.27.545266. [PMID: 37425902 PMCID: PMC10326990 DOI: 10.1101/2023.06.27.545266] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/11/2023]
Abstract
Schizophrenia (SCZ) is a highly polygenic disease and genome wide association studies have identified thousands of genetic variants that are statistically associated with this psychiatric disorder. However, our ability to translate these associations into insights on the disease mechanisms has been challenging since the causal genetic variants, their molecular function and their target genes remain largely unknown. In order to address these questions, we established a functional genomics pipeline in combination with induced pluripotent stem cell technology to functionally characterize ~35,000 non-coding genetic variants associated with schizophrenia along with their target genes. This analysis identified a set of 620 (1.7%) single nucleotide polymorphisms as functional on a molecular level in a highly cell type and condition specific fashion. These results provide a high-resolution map of functional variant-gene combinations and offer comprehensive biological insights into the developmental context and stimulation dependent molecular processes modulated by SCZ associated genetic variation.
Collapse
Affiliation(s)
- Christine K. Rummel
- Max Planck Institute of Psychiatry, Munich, Germany
- International Max Planck Research School for Translational Psychiatry (IMPRS-TP), Munich, Germany
| | - Miriam Gagliardi
- Department of Psychiatry, University of Münster, Münster, Germany
| | - Alexander Herholt
- Department of Psychiatry and Psychotherapy, University Hospital, LMU Munich, Munich, Germany
| | - Ruhel Ahmad
- Max Planck Institute of Psychiatry, Munich, Germany
| | | | | | | | | | - Laura Jimenez-Barron
- Max Planck Institute of Psychiatry, Munich, Germany
- International Max Planck Research School for Translational Psychiatry (IMPRS-TP), Munich, Germany
| | - Lucia Trastulla
- Department of Psychiatry, University of Münster, Münster, Germany
| | - Mathias Eder
- Max Planck Institute of Psychiatry, Munich, Germany
| | - Moritz Rossner
- Department of Psychiatry and Psychotherapy, University Hospital, LMU Munich, Munich, Germany
| | - Michael J. Ziller
- Max Planck Institute of Psychiatry, Munich, Germany
- Department of Psychiatry, University of Münster, Münster, Germany
- Center for Soft Nanoscience, University of Münster, Münster, Germany
| |
Collapse
|
54
|
Brown AC, Cohen CJ, Mielczarek O, Migliorini G, Costantino F, Allcock A, Davidson C, Elliott KS, Fang H, Lledó Lara A, Martin AC, Osgood JA, Sanniti A, Scozzafava G, Vecellio M, Zhang P, Black MH, Li S, Truong D, Molineros J, Howe T, Wordsworth BP, Bowness P, Knight JC. Comprehensive epigenomic profiling reveals the extent of disease-specific chromatin states and informs target discovery in ankylosing spondylitis. CELL GENOMICS 2023; 3:100306. [PMID: 37388915 PMCID: PMC10300554 DOI: 10.1016/j.xgen.2023.100306] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/29/2022] [Revised: 01/30/2023] [Accepted: 03/27/2023] [Indexed: 07/01/2023]
Abstract
Ankylosing spondylitis (AS) is a common, highly heritable inflammatory arthritis characterized by enthesitis of the spine and sacroiliac joints. Genome-wide association studies (GWASs) have revealed more than 100 genetic associations whose functional effects remain largely unresolved. Here, we present a comprehensive transcriptomic and epigenomic map of disease-relevant blood immune cell subsets from AS patients and healthy controls. We find that, while CD14+ monocytes and CD4+ and CD8+ T cells show disease-specific differences at the RNA level, epigenomic differences are only apparent upon multi-omics integration. The latter reveals enrichment at disease-associated loci in monocytes. We link putative functional SNPs to genes using high-resolution Capture-C at 10 loci, including PTGER4 and ETS1, and show how disease-specific functional genomic data can be integrated with GWASs to enhance therapeutic target discovery. This study combines epigenetic and transcriptional analysis with GWASs to identify disease-relevant cell types and gene regulation of likely pathogenic relevance and prioritize drug targets.
Collapse
Affiliation(s)
- Andrew C. Brown
- Wellcome Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK
| | - Carla J. Cohen
- Wellcome Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK
- MRC WIMM Centre for Computational Biology, MRC Weatherall Institute of Molecular Medicine, University of Oxford, John Radcliffe Hospital, Oxford OX3 9DS, UK
- Nuffield Department of Orthopaedics, Rheumatology, and Musculoskeletal Sciences, University of Oxford, Oxford OX3 7LD, UK
| | - Olga Mielczarek
- Wellcome Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK
- Horizon Discovery (PerkinElmer) Cambridge Research Park, 8100 Beach Dr., Waterbeach, Cambridge CB25 9TL, UK
| | - Gabriele Migliorini
- Wellcome Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, UK
| | - Félicie Costantino
- Wellcome Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK
- UVSQ, INSERM UMR1173, Infection et Inflammation, Laboratory of Excellence INFLAMEX, Université Paris-Saclay, Paris, France
- Rheumatology Department, AP-HP, Ambroise Paré Hospital, Paris, France
| | - Alice Allcock
- Wellcome Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK
| | - Connor Davidson
- Wellcome Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK
- Nuffield Department of Orthopaedics, Rheumatology, and Musculoskeletal Sciences, University of Oxford, Oxford OX3 7LD, UK
| | | | - Hai Fang
- Wellcome Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK
- Shanghai Institute of Hematology, State Key Laboratory of Medical Genomics, National Research Centre for Translational Medicine at Shanghai, Ruijin Hospital affiliated with Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China
| | - Alicia Lledó Lara
- Wellcome Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK
| | - Alice C. Martin
- Wellcome Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK
| | - Julie A. Osgood
- Wellcome Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK
| | - Anna Sanniti
- Wellcome Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK
| | - Giuseppe Scozzafava
- Wellcome Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK
| | - Matteo Vecellio
- Nuffield Department of Orthopaedics, Rheumatology, and Musculoskeletal Sciences, University of Oxford, Oxford OX3 7LD, UK
- Centro Ricerche Fondazione Italiana Ricerca sull’Artrite (FIRA), Fondazione Pisana per la Scienza ONLUS, Via Ferruccio Giovannini 13, 56017 San Giuliano Terme (Pisa), Italy
| | - Ping Zhang
- Wellcome Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK
- Chinese Academy of Medical Sciences Oxford Institute, Nuffield Department of Medicine, University of Oxford, Oxford OX3 7BN, UK
| | - Mary Helen Black
- Data Science, Population Analytics, Janssen R&D, Spring House, PA 19002, USA
| | - Shuwei Li
- Data Science, Population Analytics, Janssen R&D, Spring House, PA 19002, USA
| | - Dongnhu Truong
- Data Science, Population Analytics, Janssen R&D, Spring House, PA 19002, USA
| | - Julio Molineros
- Data Science, Population Analytics, Janssen R&D, Spring House, PA 19002, USA
| | - Trevor Howe
- Data Science, External Innovation, Janssen R&D, London W1G 0BG, UK
| | - B. Paul Wordsworth
- Nuffield Department of Orthopaedics, Rheumatology, and Musculoskeletal Sciences, University of Oxford, Oxford OX3 7LD, UK
- National Institute for Health Research, Comprehensive Biomedical Research Centre, Oxford OX4 2PG, UK
| | - Paul Bowness
- Nuffield Department of Orthopaedics, Rheumatology, and Musculoskeletal Sciences, University of Oxford, Oxford OX3 7LD, UK
- National Institute for Health Research, Comprehensive Biomedical Research Centre, Oxford OX4 2PG, UK
| | - Julian C. Knight
- Wellcome Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK
- Chinese Academy of Medical Sciences Oxford Institute, Nuffield Department of Medicine, University of Oxford, Oxford OX3 7BN, UK
- National Institute for Health Research, Comprehensive Biomedical Research Centre, Oxford OX4 2PG, UK
| |
Collapse
|
55
|
Nott A, Holtman IR. Genetic insights into immune mechanisms of Alzheimer's and Parkinson's disease. Front Immunol 2023; 14:1168539. [PMID: 37359515 PMCID: PMC10285485 DOI: 10.3389/fimmu.2023.1168539] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2023] [Accepted: 04/17/2023] [Indexed: 06/28/2023] Open
Abstract
Microglia, the macrophages of the brain, are vital for brain homeostasis and have been implicated in a broad range of brain disorders. Neuroinflammation has gained traction as a possible therapeutic target for neurodegeneration, however, the precise function of microglia in specific neurodegenerative disorders is an ongoing area of research. Genetic studies offer valuable insights into understanding causality, rather than merely observing a correlation. Genome-wide association studies (GWAS) have identified many genetic loci that are linked to susceptibility to neurodegenerative disorders. (Post)-GWAS studies have determined that microglia likely play an important role in the development of Alzheimer's disease (AD) and Parkinson's disease (PD). The process of understanding how individual GWAS risk loci affect microglia function and mediate susceptibility is complex. A rapidly growing number of publications with genomic datasets and computational tools have formulated new hypotheses that guide the biological interpretation of AD and PD genetic risk. In this review, we discuss the key concepts and challenges in the post-GWAS interpretation of AD and PD GWAS risk alleles. Post-GWAS challenges include the identification of target cell (sub)type(s), causal variants, and target genes. Crucially, the prediction of GWAS-identified disease-risk cell types, variants and genes require validation and functional testing to understand the biological consequences within the pathology of the disorders. Many AD and PD risk genes are highly pleiotropic and perform multiple important functions that might not be equally relevant for the mechanisms by which GWAS risk alleles exert their effect(s). Ultimately, many GWAS risk alleles exert their effect by changing microglia function, thereby altering the pathophysiology of these disorders, and hence, we believe that modelling this context is crucial for a deepened understanding of these disorders.
Collapse
Affiliation(s)
- Alexi Nott
- Department of Brain Sciences, Imperial College London, London, United Kingdom
- UK Dementia Research Institute, Imperial College London, London, United Kingdom
| | - Inge R. Holtman
- Department of Biomedical Sciences of Cells and Systems, Section Molecular Neurobiology, University of Groningen, University Medical Center Groningen, Groningen, Netherlands
| |
Collapse
|
56
|
Fabo T, Khavari P. Functional characterization of human genomic variation linked to polygenic diseases. Trends Genet 2023; 39:462-490. [PMID: 36997428 PMCID: PMC11025698 DOI: 10.1016/j.tig.2023.02.014] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2022] [Revised: 02/22/2023] [Accepted: 02/23/2023] [Indexed: 03/30/2023]
Abstract
The burden of human disease lies predominantly in polygenic diseases. Since the early 2000s, genome-wide association studies (GWAS) have identified genetic variants and loci associated with complex traits. These have ranged from variants in coding sequences to mutations in regulatory regions, such as promoters and enhancers, as well as mutations affecting mediators of mRNA stability and other downstream regulators, such as 5' and 3'-untranslated regions (UTRs), long noncoding RNA (lncRNA), and miRNA. Recent research advances in genetics have utilized a combination of computational techniques, high-throughput in vitro and in vivo screening modalities, and precise genome editing to impute the function of diverse classes of genetic variants identified through GWAS. In this review, we highlight the vastness of genomic variants associated with polygenic disease risk and address recent advances in how genetic tools can be used to functionally characterize them.
Collapse
Affiliation(s)
- Tania Fabo
- Program in Epithelial Biology, Stanford University, Stanford, CA, USA; Stanford Cancer Institute, Stanford University, Stanford, CA, USA; Graduate Program in Genetics, Stanford University, Stanford, CA, USA; Stanford University School of Medicine, Stanford University, Stanford, CA, USA
| | - Paul Khavari
- Program in Epithelial Biology, Stanford University, Stanford, CA, USA; Stanford Cancer Institute, Stanford University, Stanford, CA, USA; Graduate Program in Genetics, Stanford University, Stanford, CA, USA; Stanford University School of Medicine, Stanford University, Stanford, CA, USA; Veterans Affairs Palo Alto Healthcare System, Palo Alto, CA, USA.
| |
Collapse
|
57
|
Shi FY, Wang Y, Huang D, Liang Y, Liang N, Chen XW, Gao G. Computational Assessment of the Expression-modulating Potential for Non-coding Variants. GENOMICS, PROTEOMICS & BIOINFORMATICS 2023; 21:662-673. [PMID: 34890839 PMCID: PMC10787178 DOI: 10.1016/j.gpb.2021.10.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/27/2021] [Revised: 10/13/2021] [Accepted: 11/01/2021] [Indexed: 06/13/2023]
Abstract
Large-scale genome-wide association studies (GWAS) and expression quantitative trait locus (eQTL) studies have identified multiple non-coding variants associated with genetic diseases by affecting gene expression. However, pinpointing causal variants effectively and efficiently remains a serious challenge. Here, we developed CARMEN, a novel algorithm to identify functional non-coding expression-modulating variants. Multiple evaluations demonstrated CARMEN's superior performance over state-of-the-art tools. Applying CARMEN to GWAS and eQTL datasets further pinpointed several causal variants other than the reported lead single-nucleotide polymorphisms (SNPs). CARMEN scales well with the massive datasets, and is available online as a web server at http://carmen.gao-lab.org.
Collapse
Affiliation(s)
- Fang-Yuan Shi
- State Key Laboratory of Protein and Plant Gene Research, School of Life Sciences, Biomedical Pioneering Innovative Center (BIOPIC) & Beijing Advanced Innovation Center for Genomics (ICG), Center for Bioinformatics (CBI), Peking University, Beijing 100871, China
| | - Yu Wang
- State Key Laboratory of Protein and Plant Gene Research, School of Life Sciences, Biomedical Pioneering Innovative Center (BIOPIC) & Beijing Advanced Innovation Center for Genomics (ICG), Center for Bioinformatics (CBI), Peking University, Beijing 100871, China
| | - Dong Huang
- State Key Laboratory of Membrane Biology, Institute of Molecular Medicine, Peking University, Beijing 100871, China
| | - Yu Liang
- Human Aging Research Institute, School of Life Science, Nanchang University, Nanchang 330031, China
| | - Nan Liang
- State Key Laboratory of Protein and Plant Gene Research, School of Life Sciences, Biomedical Pioneering Innovative Center (BIOPIC) & Beijing Advanced Innovation Center for Genomics (ICG), Center for Bioinformatics (CBI), Peking University, Beijing 100871, China
| | - Xiao-Wei Chen
- State Key Laboratory of Membrane Biology, Institute of Molecular Medicine, Peking University, Beijing 100871, China; Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
| | - Ge Gao
- State Key Laboratory of Protein and Plant Gene Research, School of Life Sciences, Biomedical Pioneering Innovative Center (BIOPIC) & Beijing Advanced Innovation Center for Genomics (ICG), Center for Bioinformatics (CBI), Peking University, Beijing 100871, China.
| |
Collapse
|
58
|
Rong S, Neil CR, Welch A, Duan C, Maguire S, Meremikwu IC, Meyerson M, Evans BJ, Fairbrother WG. Large-scale functional screen identifies genetic variants with splicing effects in modern and archaic humans. Proc Natl Acad Sci U S A 2023; 120:e2218308120. [PMID: 37192163 PMCID: PMC10214146 DOI: 10.1073/pnas.2218308120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2022] [Accepted: 04/12/2023] [Indexed: 05/18/2023] Open
Abstract
Humans coexisted and interbred with other hominins which later became extinct. These archaic hominins are known to us only through fossil records and for two cases, genome sequences. Here, we engineer Neanderthal and Denisovan sequences into thousands of artificial genes to reconstruct the pre-mRNA processing patterns of these extinct populations. Of the 5,169 alleles tested in this massively parallel splicing reporter assay (MaPSy), we report 962 exonic splicing mutations that correspond to differences in exon recognition between extant and extinct hominins. Using MaPSy splicing variants, predicted splicing variants, and splicing quantitative trait loci, we show that splice-disrupting variants experienced greater purifying selection in anatomically modern humans than that in Neanderthals. Adaptively introgressed variants were enriched for moderate-effect splicing variants, consistent with positive selection for alternative spliced alleles following introgression. As particularly compelling examples, we characterized a unique tissue-specific alternative splicing variant at the adaptively introgressed innate immunity gene TLR1, as well as a unique Neanderthal introgressed alternative splicing variant in the gene HSPG2 that encodes perlecan. We further identified potentially pathogenic splicing variants found only in Neanderthals and Denisovans in genes related to sperm maturation and immunity. Finally, we found splicing variants that may contribute to variation among modern humans in total bilirubin, balding, hemoglobin levels, and lung capacity. Our findings provide unique insights into natural selection acting on splicing in human evolution and demonstrate how functional assays can be used to identify candidate causal variants underlying differences in gene regulation and phenotype.
Collapse
Affiliation(s)
- Stephen Rong
- Center for Computational Molecular Biology, Brown University, Providence, RI02912
- Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI02912
| | - Christopher R. Neil
- Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI02912
| | - Anastasia Welch
- Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI02912
| | - Chaorui Duan
- Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI02912
| | - Samantha Maguire
- Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI02912
| | - Ijeoma C. Meremikwu
- Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI02912
| | - Malcolm Meyerson
- Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI02912
| | - Ben J. Evans
- Department of Biology, McMaster University, Hamilton, ONL8S 4K1, Canada
| | - William G. Fairbrother
- Center for Computational Molecular Biology, Brown University, Providence, RI02912
- Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI02912
- Hassenfeld Child Health Innovation Institute of Brown University, Providence, RI02912
| |
Collapse
|
59
|
Das M, Hossain A, Banerjee D, Praul CA, Girirajan S. Challenges and considerations for reproducibility of STARR-seq assays. Genome Res 2023; 33:479-495. [PMID: 37130797 PMCID: PMC10234304 DOI: 10.1101/gr.277204.122] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2022] [Accepted: 03/15/2023] [Indexed: 05/04/2023]
Abstract
High-throughput methods such as RNA-seq, ChIP-seq, and ATAC-seq have well-established guidelines, commercial kits, and analysis pipelines that enable consistency and wider adoption for understanding genome function and regulation. STARR-seq, a popular assay for directly quantifying the activities of thousands of enhancer sequences simultaneously, has seen limited standardization across studies. The assay is long, with more than 250 steps, and frequent customization of the protocol and variations in bioinformatics methods raise concerns for reproducibility of STARR-seq studies. Here, we assess each step of the protocol and analysis pipelines from published sources and in-house assays, and identify critical steps and quality control (QC) checkpoints necessary for reproducibility of the assay. We also provide guidelines for experimental design, protocol scaling, customization, and analysis pipelines for better adoption of the assay. These resources will allow better optimization of STARR-seq for specific research needs, enable comparisons and integration across studies, and improve the reproducibility of results.
Collapse
Affiliation(s)
- Maitreya Das
- Department of Biochemistry and Molecular Biology, Pennsylvania State University, University Park, Pennsylvania 16802, USA;
- Molecular and Cellular Integrative Biosciences Graduate Program, Pennsylvania State University, University Park, Pennsylvania 16802, USA
- Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - Ayaan Hossain
- Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, Pennsylvania 16802, USA
- Bioinformatics and Genomics Graduate Program, Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - Deepro Banerjee
- Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, Pennsylvania 16802, USA
- Bioinformatics and Genomics Graduate Program, Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - Craig Alan Praul
- Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - Santhosh Girirajan
- Department of Biochemistry and Molecular Biology, Pennsylvania State University, University Park, Pennsylvania 16802, USA;
- Molecular and Cellular Integrative Biosciences Graduate Program, Pennsylvania State University, University Park, Pennsylvania 16802, USA
- Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, Pennsylvania 16802, USA
- Bioinformatics and Genomics Graduate Program, Pennsylvania State University, University Park, Pennsylvania 16802, USA
- Department of Anthropology, Pennsylvania State University, University Park, Pennsylvania 16802, USA
| |
Collapse
|
60
|
Lu H, Ma L, Quan C, Li L, Lu Y, Zhou G, Zhang C. RegVar: Tissue-specific Prioritization of Non-coding Regulatory Variants. GENOMICS, PROTEOMICS & BIOINFORMATICS 2023; 21:385-395. [PMID: 34973416 PMCID: PMC10626172 DOI: 10.1016/j.gpb.2021.08.011] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/11/2021] [Revised: 06/11/2021] [Accepted: 09/27/2021] [Indexed: 06/14/2023]
Abstract
Non-coding genomic variants constitute the majority of trait-associated genome variations; however, the identification of functional non-coding variants is still a challenge in human genetics, and a method for systematically assessing the impact of regulatory variants on gene expression and linking these regulatory variants to potential target genes is still lacking. Here, we introduce a deep neural network (DNN)-based computational framework, RegVar, which can accurately predict the tissue-specific impact of non-coding regulatory variants on target genes. We show that by robustly learning the genomic characteristics of massive variant-gene expression associations in a variety of human tissues, RegVar vastly surpasses all current non-coding variant prioritization methods in predicting regulatory variants under different circumstances. The unique features of RegVar make it an excellent framework for assessing the regulatory impact of any variant on its putative target genes in a variety of tissues. RegVar is available as a web server at https://regvar.omic.tech/.
Collapse
Affiliation(s)
- Hao Lu
- Beijing Institute of Radiation Medicine, State Key Laboratory of Proteomics, Beijing 100850, China
| | - Luyu Ma
- Beijing Institute of Radiation Medicine, State Key Laboratory of Proteomics, Beijing 100850, China
| | - Cheng Quan
- Beijing Institute of Radiation Medicine, State Key Laboratory of Proteomics, Beijing 100850, China
| | - Lei Li
- Beijing Institute of Radiation Medicine, State Key Laboratory of Proteomics, Beijing 100850, China
| | - Yiming Lu
- Beijing Institute of Radiation Medicine, State Key Laboratory of Proteomics, Beijing 100850, China.
| | - Gangqiao Zhou
- Beijing Institute of Radiation Medicine, State Key Laboratory of Proteomics, Beijing 100850, China.
| | - Chenggang Zhang
- Beijing Institute of Radiation Medicine, State Key Laboratory of Proteomics, Beijing 100850, China.
| |
Collapse
|
61
|
Seah C, Huckins LM, Brennand KJ. Stem Cell Models for Context-Specific Modeling in Psychiatric Disorders. Biol Psychiatry 2023; 93:642-650. [PMID: 36658083 DOI: 10.1016/j.biopsych.2022.09.033] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/02/2022] [Revised: 09/27/2022] [Accepted: 09/27/2022] [Indexed: 01/21/2023]
Abstract
Genome-wide association studies reveal the complex polygenic architecture underlying psychiatric disorder risk, but there is an unmet need to validate causal variants, resolve their target genes(s), and explore their functional impacts on disorder-related mechanisms. Disorder-associated loci regulate transcription of target genes in a cell type- and context-specific manner, which can be measured through expression quantitative trait loci. In this review, we discuss methods and insights from context-specific modeling of genetically and environmentally regulated expression. Human induced pluripotent stem cell-derived cell type and organoid models have uncovered context-specific psychiatric disorder associations by investigating tissue-, cell type-, sex-, age-, and stressor-specific genetic regulation of expression. Techniques such as massively parallel reporter assays and pooled CRISPR (clustered regularly interspaced short palindromic repeats) screens make it possible to functionally fine-map genome-wide association study loci and validate their target genes at scale. Integration of disorder-associated contexts with these patient-specific human induced pluripotent stem cell models makes it possible to uncover gene by environment interactions that mediate disorder risk, which will ultimately improve our ability to diagnose and treat psychiatric disorders.
Collapse
Affiliation(s)
- Carina Seah
- Pamela Sklar Division of Psychiatric Genomics, Icahn School of Medicine at Mount Sinai, New York; Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York; Nash Family Department of Neuroscience, Friedman Brain Institute, Icahn School of Medicine at Mount Sinai, New York
| | - Laura M Huckins
- Pamela Sklar Division of Psychiatric Genomics, Icahn School of Medicine at Mount Sinai, New York; Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York; Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York; Seaver Autism Center for Research and Treatment, Icahn School of Medicine at Mount Sinai, New York; Department of Psychiatry, Yale University School of Medicine, New Haven, Connecticut.
| | - Kristen J Brennand
- Pamela Sklar Division of Psychiatric Genomics, Icahn School of Medicine at Mount Sinai, New York; Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York; Department of Psychiatry, Icahn School of Medicine at Mount Sinai, New York; Nash Family Department of Neuroscience, Friedman Brain Institute, Icahn School of Medicine at Mount Sinai, New York; Department of Psychiatry, Yale University School of Medicine, New Haven, Connecticut.
| |
Collapse
|
62
|
Jeong R, Bulyk ML. Colocalization of blood cell traits GWAS associations and variation in PU.1 genomic occupancy prioritizes causal noncoding regulatory variants. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.03.29.534582. [PMID: 37034747 PMCID: PMC10081269 DOI: 10.1101/2023.03.29.534582] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
Genome-wide association studies (GWAS) have uncovered numerous trait-associated loci across the human genome, most of which are located in noncoding regions, making interpretations difficult. Moreover, causal variants are hard to statistically fine-map at many loci because of widespread linkage disequilibrium. To address this challenge, we present a strategy utilizing transcription factor (TF) binding quantitative trait loci (bQTLs) for colocalization analysis to identify trait associations likely mediated by TF occupancy variation and to pinpoint likely causal variants using motif scores. We applied this approach to PU.1 bQTLs in lymphoblastoid cell lines and blood cell traits GWAS data. Colocalization analysis revealed 69 blood cell trait GWAS loci putatively driven by PU.1 occupancy variation. We nominate PU.1 motif-altering variants as the likely shared causal variants at 51 loci. Such integration of TF bQTL data with other GWAS data may reveal transcriptional regulatory mechanisms and causal noncoding variants underlying additional complex traits.
Collapse
Affiliation(s)
- Raehoon Jeong
- Division of Genetics, Department of Medicine, Brigham and Women’s Hospital and Harvard Medical School, Boston, MA 02115, USA
- Bioinformatics and Integrative Genomics Graduate Program, Harvard University, Cambridge, MA 02138, USA
| | - Martha L. Bulyk
- Division of Genetics, Department of Medicine, Brigham and Women’s Hospital and Harvard Medical School, Boston, MA 02115, USA
- Bioinformatics and Integrative Genomics Graduate Program, Harvard University, Cambridge, MA 02138, USA
- Department of Pathology, Brigham and Women’s Hospital and Harvard Medical School, Boston, MA 02115, USA
| |
Collapse
|
63
|
Lopdell TJ. Using QTL to Identify Genes and Pathways Underlying the Regulation and Production of Milk Components in Cattle. Animals (Basel) 2023; 13:ani13050911. [PMID: 36899768 PMCID: PMC10000085 DOI: 10.3390/ani13050911] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2023] [Revised: 02/23/2023] [Accepted: 02/28/2023] [Indexed: 03/06/2023] Open
Abstract
Milk is a complex liquid, and the concentrations of many of its components are under genetic control. Many genes and pathways are known to regulate milk composition, and the purpose of this review is to highlight how the discoveries of quantitative trait loci (QTL) for milk phenotypes can elucidate these pathways. The main body of this review focuses primarily on QTL discovered in cattle (Bos taurus) as a model species for the biology of lactation, and there are occasional references to sheep genetics. The following section describes a range of techniques that can be used to help identify the causative genes underlying QTL when the underlying mechanism involves the regulation of gene expression. As genotype and phenotype databases continue to grow and diversify, new QTL will continue to be discovered, and although proving the causality of underlying genes and variants remains difficult, these new data sets will further enhance our understanding of the biology of lactation.
Collapse
|
64
|
Gundry M, Sankaran VG. Hacking hematopoiesis - emerging tools for examining variant effects. Dis Model Mech 2023; 16:288409. [PMID: 36826849 PMCID: PMC9983777 DOI: 10.1242/dmm.049857] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/25/2023] Open
Abstract
Hematopoiesis is a continuous process of blood and immune cell production. It is orchestrated by thousands of gene products that respond to extracellular signals by guiding cell fate decisions to meet the needs of the organism. Although much of our knowledge of this process comes from work in model systems, we have learned a great deal from studies on human genetic variation. Considerable insight has emerged from studies on presumed monogenic blood disorders, which continue to provide key insights into the mechanisms critical for hematopoiesis. Furthermore, the emergence of large-scale biobanks and cohorts has uncovered thousands of genomic loci associated with blood cell traits and diseases. Some of these blood cell trait-associated loci act as modifiers of what were once thought to be monogenic blood diseases. However, most of these loci await functional validation. Here, we discuss the validation bottleneck and emerging methods to more effectively connect variant to function. In particular, we highlight recent innovations in genome editing, which have paved the path forward for high-throughput functional assessment of loci. Finally, we discuss existing barriers to progress, including challenges in manipulating the genomes of primary hematopoietic cells.
Collapse
Affiliation(s)
- Michael Gundry
- Division of Hematology/Oncology, Boston Children's Hospital, Harvard Medical School, Boston, MA 02115, USA
- Department of Pediatric Oncology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA 02115, USA
- Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA
| | - Vijay G. Sankaran
- Division of Hematology/Oncology, Boston Children's Hospital, Harvard Medical School, Boston, MA 02115, USA
- Department of Pediatric Oncology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA 02115, USA
- Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA
- Author for correspondence ()
| |
Collapse
|
65
|
Stankey CT, Lee JC. Translating non-coding genetic associations into a better understanding of immune-mediated disease. Dis Model Mech 2023; 16:297044. [PMID: 36897113 PMCID: PMC10040244 DOI: 10.1242/dmm.049790] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/11/2023] Open
Abstract
Genome-wide association studies have identified hundreds of genetic loci that are associated with immune-mediated diseases. Most disease-associated variants are non-coding, and a large proportion of these variants lie within enhancers. As a result, there is a pressing need to understand how common genetic variation might affect enhancer function and thereby contribute to immune-mediated (and other) diseases. In this Review, we first describe statistical and experimental methods to identify causal genetic variants that modulate gene expression, including statistical fine-mapping and massively parallel reporter assays. We then discuss approaches to characterise the mechanisms by which these variants modulate immune function, such as clustered regularly interspaced short palindromic repeats (CRISPR)-based screens. We highlight examples of studies that, by elucidating the effects of disease variants within enhancers, have provided important insights into immune function and uncovered key pathways of disease.
Collapse
Affiliation(s)
- Christina T Stankey
- Genetic Mechanisms of Disease Laboratory, The Francis Crick Institute, London NW1 1AT, UK
- Department of Immunology and Inflammation, Imperial College London, London W12 0NN, UK
| | - James C Lee
- Genetic Mechanisms of Disease Laboratory, The Francis Crick Institute, London NW1 1AT, UK
- Institute of Liver and Digestive Health, Royal Free Hospital, University College London, London NW3 2PF, UK
| |
Collapse
|
66
|
Carrasco Pro S, Hook H, Bray D, Berenzy D, Moyer D, Yin M, Labadorf AT, Tewhey R, Siggers T, Fuxman Bass JI. Widespread perturbation of ETS factor binding sites in cancer. Nat Commun 2023; 14:913. [PMID: 36808133 PMCID: PMC9938127 DOI: 10.1038/s41467-023-36535-8] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2022] [Accepted: 02/03/2023] [Indexed: 02/19/2023] Open
Abstract
Although >90% of somatic mutations reside in non-coding regions, few have been reported as cancer drivers. To predict driver non-coding variants (NCVs), we present a transcription factor (TF)-aware burden test based on a model of coherent TF function in promoters. We apply this test to NCVs from the Pan-Cancer Analysis of Whole Genomes cohort and predict 2555 driver NCVs in the promoters of 813 genes across 20 cancer types. These genes are enriched in cancer-related gene ontologies, essential genes, and genes associated with cancer prognosis. We find that 765 candidate driver NCVs alter transcriptional activity, 510 lead to differential binding of TF-cofactor regulatory complexes, and that they primarily impact the binding of ETS factors. Finally, we show that different NCVs within a promoter often affect transcriptional activity through shared mechanisms. Our integrated computational and experimental approach shows that cancer NCVs are widespread and that ETS factors are commonly disrupted.
Collapse
Affiliation(s)
| | - Heather Hook
- Department of Biology, Boston University, Boston, MA, USA
| | - David Bray
- Bioinformatics Program, Boston University, Boston, MA, USA
| | | | - Devlin Moyer
- Bioinformatics Program, Boston University, Boston, MA, USA
| | - Meimei Yin
- Department of Biology, Boston University, Boston, MA, USA
| | - Adam Thomas Labadorf
- Bioinformatics Hub, Boston University, Boston, MA, USA
- Boston University School of Medicine, Department of Neurology, Boston, MA, USA
| | | | - Trevor Siggers
- Bioinformatics Program, Boston University, Boston, MA, USA.
- Department of Biology, Boston University, Boston, MA, USA.
- Biological Design Center, Boston University, Boston, MA, USA.
| | - Juan Ignacio Fuxman Bass
- Bioinformatics Program, Boston University, Boston, MA, USA.
- Department of Biology, Boston University, Boston, MA, USA.
| |
Collapse
|
67
|
Deng C, Whalen S, Steyert M, Ziffra R, Przytycki PF, Inoue F, Pereira DA, Capauto D, Norton S, Vaccarino FM, Pollen A, Nowakowski TJ, Ahituv N, Pollard KS. Massively parallel characterization of psychiatric disorder-associated and cell-type-specific regulatory elements in the developing human cortex. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.02.15.528663. [PMID: 36824845 PMCID: PMC9949039 DOI: 10.1101/2023.02.15.528663] [Citation(s) in RCA: 13] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/17/2023]
Abstract
Nucleotide changes in gene regulatory elements are important determinants of neuronal development and disease. Using massively parallel reporter assays in primary human cells from mid-gestation cortex and cerebral organoids, we interrogated the cis-regulatory activity of 102,767 sequences, including differentially accessible cell-type specific regions in the developing cortex and single-nucleotide variants associated with psychiatric disorders. In primary cells, we identified 46,802 active enhancer sequences and 164 disorder-associated variants that significantly alter enhancer activity. Activity was comparable in organoids and primary cells, suggesting that organoids provide an adequate model for the developing cortex. Using deep learning, we decoded the sequence basis and upstream regulators of enhancer activity. This work establishes a comprehensive catalog of functional gene regulatory elements and variants in human neuronal development.
Collapse
Affiliation(s)
- Chengyu Deng
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco; San Francisco, CA, USA
- Institute for Human Genetics, University of California, San Francisco; San Francisco, CA, USA
| | - Sean Whalen
- Gladstone Institutes; San Francisco, CA, USA
| | - Marilyn Steyert
- Department of Anatomy, University of California, San Francisco; San Francisco, CA, USA
- Department of Psychiatry, University of California, San Francisco; San Francisco, CA, USA
- Department of Neurological Surgery, University of California, San Francisco; San Francisco, CA, USA
| | - Ryan Ziffra
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco; San Francisco, CA, USA
| | | | - Fumitaka Inoue
- Institute for the Advanced Study of Human Biology (WPI-ASHBi), Kyoto University; Kyoto, Japan
| | - Daniela A. Pereira
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco; San Francisco, CA, USA
- Institute for Human Genetics, University of California, San Francisco; San Francisco, CA, USA
- Graduate Program of Genetics, Institute of Biological Sciences, Federal University of Minas Gerais; Belo Horizonte, Minas Gerais, Brazil
| | | | - Scott Norton
- Child Study Center, Yale University; New Haven, CT, USA
| | - Flora M. Vaccarino
- Child Study Center, Yale University; New Haven, CT, USA
- Department of Neuroscience, Yale University; New Haven, CT, USA
| | - Alex Pollen
- Department of Neurology, University of California, San Francisco; San Francisco, CA, USA
- Eli and Edythe Broad Center for Regeneration Medicine and Stem Cell Research, University of California, San Francisco; San Francisco, CA, USA
| | - Tomasz J. Nowakowski
- Department of Anatomy, University of California, San Francisco; San Francisco, CA, USA
- Department of Psychiatry, University of California, San Francisco; San Francisco, CA, USA
- Department of Neurological Surgery, University of California, San Francisco; San Francisco, CA, USA
- Eli and Edythe Broad Center for Regeneration Medicine and Stem Cell Research, University of California, San Francisco; San Francisco, CA, USA
- Chan Zuckerberg Biohub, San Francisco; San Francisco, CA, USA
| | - Nadav Ahituv
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco; San Francisco, CA, USA
- Institute for Human Genetics, University of California, San Francisco; San Francisco, CA, USA
| | - Katherine S. Pollard
- Institute for Human Genetics, University of California, San Francisco; San Francisco, CA, USA
- Gladstone Institutes; San Francisco, CA, USA
- Chan Zuckerberg Biohub, San Francisco; San Francisco, CA, USA
- Department of Epidemiology and Biostatistics, University of California, San Francisco; San Francisco, CA, USA
| |
Collapse
|
68
|
Gallego Romero I, Lea AJ. Leveraging massively parallel reporter assays for evolutionary questions. Genome Biol 2023; 24:26. [PMID: 36788564 PMCID: PMC9926830 DOI: 10.1186/s13059-023-02856-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2022] [Accepted: 01/17/2023] [Indexed: 02/16/2023] Open
Abstract
A long-standing goal of evolutionary biology is to decode how gene regulation contributes to organismal diversity. Doing so is challenging because it is hard to predict function from non-coding sequence and to perform molecular research with non-model taxa. Massively parallel reporter assays (MPRAs) enable the testing of thousands to millions of sequences for regulatory activity simultaneously. Here, we discuss the execution, advantages, and limitations of MPRAs, with a focus on evolutionary questions. We propose solutions for extending MPRAs to rare taxa and those with limited genomic resources, and we underscore MPRA's broad potential for driving genome-scale, functional studies across organisms.
Collapse
Affiliation(s)
- Irene Gallego Romero
- Melbourne Integrative Genomics, University of Melbourne, Royal Parade, Parkville, Victoria, 3010, Australia. .,School of BioSciences, The University of Melbourne, Royal Parade, Parkville, 3010, Australia. .,The Centre for Stem Cell Systems, Faculty of Medicine, Dentistry and Health Sciences, The University of Melbourne, 30 Royal Parade, Parkville, Victoria, 3010, Australia. .,Center for Genomics, Evolution and Medicine, Institute of Genomics, University of Tartu, Riia 23b, 51010, Tartu, Estonia.
| | - Amanda J. Lea
- grid.152326.10000 0001 2264 7217Department of Biological Sciences, Vanderbilt University, Nashville, TN 37240 USA ,grid.152326.10000 0001 2264 7217Vanderbilt Genetics Institute, Vanderbilt University, Nashville, TN 37240 USA ,grid.152326.10000 0001 2264 7217Evolutionary Studies Initiative, Vanderbilt University, Nashville, TN 37240 USA ,Child and Brain Development Program, Canadian Institute for Advanced Study, Toronto, Canada
| |
Collapse
|
69
|
Wang QS, Ui-Tei K. Utilizing Large Functional and Population Genomics Resources for CRISPR/Cas Perturbation Experiment Design. Methods Mol Biol 2023; 2637:63-73. [PMID: 36773138 DOI: 10.1007/978-1-0716-3016-7_5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/12/2023]
Abstract
Genome sequencing technologies have rapidly evolved in the past decades, enabling us to interpret the human genome through multiple perspectives, ranging from cross-species comparisons, naturally occurring variation in health and disease state to regulatory mechanisms.Although such perspectives are all informative to narrow down the list of genes or variants for perturbation experiments based on specific biological aims, utilizing multiple sources of information is often challenging in practice.In this chapter, we provide an overview of major large-scale functional and population genomics resources, followed by a practical example of selecting target variants for genetic perturbation experiments involving genome engineering techniques such as CRISPR/Cas.
Collapse
Affiliation(s)
- Qingbo S Wang
- Department of Statistical Genetics, Graduate School of Medicine, Osaka University, Suita, Osaka, Japan.
- Laboratory of Statistical Immunology, Immunology Frontier Research Center (WPI-IFReC), Osaka University, Suita, Osaka, Japan.
| | - Kumiko Ui-Tei
- Department of Computational Biology and Medical Sciences, Graduate School of Frontier Sciences, The University of Tokyo, Chiba, Japan.
- Department of Biological Sciences, Graduate School of Science, The University of Tokyo, Tokyo, Japan.
| |
Collapse
|
70
|
Jagoda E, Marnetto D, Senevirathne G, Gonzalez V, Baid K, Montinaro F, Richard D, Falzarano D, LeBlanc EV, Colpitts CC, Banerjee A, Pagani L, Capellini TD. Regulatory dissection of the severe COVID-19 risk locus introgressed by Neanderthals. eLife 2023; 12:e71235. [PMID: 36763080 PMCID: PMC9917435 DOI: 10.7554/elife.71235] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2021] [Accepted: 01/26/2023] [Indexed: 02/11/2023] Open
Abstract
Individuals infected with the SARS-CoV-2 virus present with a wide variety of symptoms ranging from asymptomatic to severe and even lethal outcomes. Past research has revealed a genetic haplotype on chromosome 3 that entered the human population via introgression from Neanderthals as the strongest genetic risk factor for the severe response to COVID-19. However, the specific variants along this introgressed haplotype that contribute to this risk and the biological mechanisms that are involved remain unclear. Here, we assess the variants present on the risk haplotype for their likelihood of driving the genetic predisposition to severe COVID-19 outcomes. We do this by first exploring their impact on the regulation of genes involved in COVID-19 infection using a variety of population genetics and functional genomics tools. We then perform a locus-specific massively parallel reporter assay to individually assess the regulatory potential of each allele on the haplotype in a multipotent immune-related cell line. We ultimately reduce the set of over 600 linked genetic variants to identify four introgressed alleles that are strong functional candidates for driving the association between this locus and severe COVID-19. Using reporter assays in the presence/absence of SARS-CoV-2, we find evidence that these variants respond to viral infection. These variants likely drive the locus' impact on severity by modulating the regulation of two critical chemokine receptor genes: CCR1 and CCR5. These alleles are ideal targets for future functional investigations into the interaction between host genomics and COVID-19 outcomes.
Collapse
Affiliation(s)
- Evelyn Jagoda
- Department of Human Evolutionary Biology, Harvard UniversityCambridgeUnited States
| | - Davide Marnetto
- Estonian Biocentre, Institute of Genomics, University of TartuTartuEstonia
| | - Gayani Senevirathne
- Department of Human Evolutionary Biology, Harvard UniversityCambridgeUnited States
| | - Victoria Gonzalez
- Department of Veterinary Microbiology, University of SaskatchewanSaskatoonCanada
- Vaccine and Infectious Disease Organization, University of SaskatchewanSaskatoonCanada
| | - Kaushal Baid
- Vaccine and Infectious Disease Organization, University of SaskatchewanSaskatoonCanada
| | - Francesco Montinaro
- Estonian Biocentre, Institute of Genomics, University of TartuTartuEstonia
- Department of Biology, University of BariBariItaly
| | - Daniel Richard
- Department of Human Evolutionary Biology, Harvard UniversityCambridgeUnited States
| | - Darryl Falzarano
- Department of Veterinary Microbiology, University of SaskatchewanSaskatoonCanada
- Vaccine and Infectious Disease Organization, University of SaskatchewanSaskatoonCanada
| | - Emmanuelle V LeBlanc
- Department of Biomedical and Molecular Sciences, Queen’s UniversityKingstonCanada
| | - Che C Colpitts
- Department of Biomedical and Molecular Sciences, Queen’s UniversityKingstonCanada
| | - Arinjay Banerjee
- Department of Veterinary Microbiology, University of SaskatchewanSaskatoonCanada
- Vaccine and Infectious Disease Organization, University of SaskatchewanSaskatoonCanada
- Department of Biology, University of WaterlooWaterlooCanada
- Department of Laboratory Medicine and Pathobiology, University of TorontoTorontoCanada
| | - Luca Pagani
- Estonian Biocentre, Institute of Genomics, University of TartuTartuEstonia
- Department of Biology, University of PadovaPadovaItaly
| | - Terence D Capellini
- Department of Human Evolutionary Biology, Harvard UniversityCambridgeUnited States
- Broad Institute of MIT and HarvardCambridgeUnited States
| |
Collapse
|
71
|
Stikker BS, Hendriks RW, Stadhouders R. Decoding the genetic and epigenetic basis of asthma. Allergy 2023; 78:940-956. [PMID: 36727912 DOI: 10.1111/all.15666] [Citation(s) in RCA: 20] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2022] [Revised: 01/17/2023] [Accepted: 01/30/2023] [Indexed: 02/03/2023]
Abstract
Asthma is a complex and heterogeneous chronic inflammatory disease of the airways. Alongside environmental factors, asthma susceptibility is strongly influenced by genetics. Given its high prevalence and our incomplete understanding of the mechanisms underlying disease susceptibility, asthma is frequently studied in genome-wide association studies (GWAS), which have identified thousands of genetic variants associated with asthma development. Virtually all these genetic variants reside in non-coding genomic regions, which has obscured the functional impact of asthma-associated variants and their translation into disease-relevant mechanisms. Recent advances in genomics technology and epigenetics now offer methods to link genetic variants to gene regulatory elements embedded within non-coding regions, which have started to unravel the molecular mechanisms underlying the complex (epi)genetics of asthma. Here, we provide an integrated overview of (epi)genetic variants associated with asthma, focusing on efforts to link these disease associations to biological insight into asthma pathophysiology using state-of-the-art genomics methodology. Finally, we provide a perspective as to how decoding the genetic and epigenetic basis of asthma has the potential to transform clinical management of asthma and to predict the risk of asthma development.
Collapse
Affiliation(s)
- Bernard S Stikker
- Department of Pulmonary Medicine, Erasmus MC, University Medical Center, Rotterdam, The Netherlands
| | - Rudi W Hendriks
- Department of Pulmonary Medicine, Erasmus MC, University Medical Center, Rotterdam, The Netherlands
| | - Ralph Stadhouders
- Department of Pulmonary Medicine, Erasmus MC, University Medical Center, Rotterdam, The Netherlands.,Department of Cell Biology, Erasmus MC, University Medical Center, Rotterdam, The Netherlands
| |
Collapse
|
72
|
A single-cell massively parallel reporter assay detects cell-type-specific gene regulation. Nat Genet 2023; 55:346-354. [PMID: 36635387 PMCID: PMC9931678 DOI: 10.1038/s41588-022-01278-7] [Citation(s) in RCA: 24] [Impact Index Per Article: 24.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2021] [Accepted: 12/05/2022] [Indexed: 01/14/2023]
Abstract
Massively parallel reporter gene assays are key tools in regulatory genomics but cannot be used to identify cell-type-specific regulatory elements without performing assays serially across different cell types. To address this problem, we developed a single-cell massively parallel reporter assay (scMPRA) to measure the activity of libraries of cis-regulatory sequences (CRSs) across multiple cell types simultaneously. We assayed a library of core promoters in a mixture of HEK293 and K562 cells and showed that scMPRA is a reproducible, highly parallel, single-cell reporter gene assay that detects cell-type-specific cis-regulatory activity. We then measured a library of promoter variants across multiple cell types in live mouse retinas and showed that subtle genetic variants can produce cell-type-specific effects on cis-regulatory activity. We anticipate that scMPRA will be widely applicable for studying the role of CRSs across diverse cell types.
Collapse
|
73
|
Mouri K, Dewey HB, Castro R, Berenzy D, Kales S, Tewhey R. Whole-genome functional characterization of RE1 silencers using a modified massively parallel reporter assay. CELL GENOMICS 2023; 3:100234. [PMID: 36777181 PMCID: PMC9903721 DOI: 10.1016/j.xgen.2022.100234] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/11/2022] [Revised: 09/12/2022] [Accepted: 11/23/2022] [Indexed: 12/23/2022]
Abstract
Both upregulation and downregulation by cis-regulatory elements help modulate precise gene expression. However, our understanding of repressive elements is far more limited than activating elements. To address this gap, we characterized RE1, a group of transcriptional silencers bound by REST, at genome-wide scale using a modified massively parallel reporter assay (MPRAduo). MPRAduo empirically defined a minimal binding strength of REST (REST motif-intrinsic value [m-value]), above which cofactors colocalize and silence transcription. We identified 1,500 human variants that alter RE1 silencing and found that their effect sizes are predictable when they overlap with REST-binding sites above the m-value. Additionally, we demonstrate that non-canonical REST-binding motifs exhibit silencer function only if they precisely align half sites with specific spacer lengths. Our results show mechanistic insights into RE1, which allow us to predict its activity and effect of variants on RE1, providing a paradigm for performing genome-wide functional characterization of transcription-factor-binding sites.
Collapse
Affiliation(s)
| | | | | | | | - Susan Kales
- The Jackson Laboratory, Bar Harbor, ME 04609, USA
| | - Ryan Tewhey
- The Jackson Laboratory, Bar Harbor, ME 04609, USA
- Graduate School of Biomedical Sciences and Engineering, University of Maine, Orono, ME, USA
- Graduate School of Biomedical Sciences, Tufts University School of Medicine, Boston, MA, USA
| |
Collapse
|
74
|
Kircher M, Ludwig KU. Systematic assays and resources for the functional annotation of non-coding variants. MED GENET-BERLIN 2022; 34:275-286. [PMID: 37034418 PMCID: PMC10081529 DOI: 10.1515/medgen-2022-2161] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Abstract
Identification of genetic variation in individual genomes is now a routine procedure in human genetic research and diagnostics. For many variants, however, insufficient evidence is available to establish a pathogenic effect, particularly for variants in non-coding regions. Furthermore, the sheer number of candidate variants renders testing in individual assays virtually impossible. While scalable approaches are being developed, the selection of methods and resources and the application of a given framework to a particular disease or trait remain major challenges. This limits the translation of results from both genome-wide association studies and genome sequencing. Here, we discuss computational and experimental approaches available for functional annotation of non-coding variation.
Collapse
Affiliation(s)
- Martin Kircher
- Institute of Human Genetics , University of Lübeck , Lübeck , Germany
- Berlin Institute of Health at Charité – Universitätsmedizin Berlin , Berlin , Germany
| | - Kerstin U. Ludwig
- Institute of Human Genetics, University Hospital Bonn , University of Bonn , Venusberg-Campus 1, Building 76 , Bonn , Germany
| |
Collapse
|
75
|
Morova T, Ding Y, Huang CCF, Sar F, Schwarz T, Giambartolomei C, Baca S, Grishin D, Hach F, Gusev A, Freedman M, Pasaniuc B, Lack N. Optimized high-throughput screening of non-coding variants identified from genome-wide association studies. Nucleic Acids Res 2022; 51:e18. [PMID: 36546757 PMCID: PMC9943666 DOI: 10.1093/nar/gkac1198] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2022] [Revised: 11/19/2022] [Accepted: 12/06/2022] [Indexed: 12/24/2022] Open
Abstract
The vast majority of disease-associated single nucleotide polymorphisms (SNP) identified from genome-wide association studies (GWAS) are localized in non-coding regions. A significant fraction of these variants impact transcription factors binding to enhancer elements and alter gene expression. To functionally interrogate the activity of such variants we developed snpSTARRseq, a high-throughput experimental method that can interrogate the functional impact of hundreds to thousands of non-coding variants on enhancer activity. snpSTARRseq dramatically improves signal-to-noise by utilizing a novel sequencing and bioinformatic approach that increases both insert size and the number of variants tested per loci. Using this strategy, we interrogated known prostate cancer (PCa) risk-associated loci and demonstrated that 35% of them harbor SNPs that significantly altered enhancer activity. Combining these results with chromosomal looping data we could identify interacting genes and provide a mechanism of action for 20 PCa GWAS risk regions. When benchmarked to orthogonal methods, snpSTARRseq showed a strong correlation with in vivo experimental allelic-imbalance studies whereas there was no correlation with predictive in silico approaches. Overall, snpSTARRseq provides an integrated experimental and computational framework to functionally test non-coding genetic variants.
Collapse
Affiliation(s)
- Tunc Morova
- Vancouver Prostate Centre, Vancouver, BC V6H 3Z6, Canada
| | - Yi Ding
- Bioinformatics Interdepartmental Program, University of California, Los Angeles, Los Angeles, CA 90095, USA
| | | | - Funda Sar
- Vancouver Prostate Centre, Vancouver, BC V6H 3Z6, Canada
| | - Tommer Schwarz
- Bioinformatics Interdepartmental Program, University of California, Los Angeles, Los Angeles, CA 90095, USA
| | - Claudia Giambartolomei
- Central RNA Lab, Istituto Italiano di Tecnologia, Genova 16163, Italy,Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, CA 90095, USA
| | - Sylvan C Baca
- Department of Medical Oncology, The Center for Functional Cancer Epigenetics, Dana Farber Cancer Institute, Boston, MA 02215, USA
| | - Dennis Grishin
- Department of Medical Oncology, The Center for Functional Cancer Epigenetics, Dana Farber Cancer Institute, Boston, MA 02215, USA
| | - Faraz Hach
- Vancouver Prostate Centre, Vancouver, BC V6H 3Z6, Canada,Department of Urologic Science, University of British Columbia, Vancouver, BC V5Z 1M9, Canada
| | - Alexander Gusev
- Department of Medical Oncology, The Center for Functional Cancer Epigenetics, Dana Farber Cancer Institute, Boston, MA 02215, USA,Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA
| | - Matthew L Freedman
- Department of Medical Oncology, The Center for Functional Cancer Epigenetics, Dana Farber Cancer Institute, Boston, MA 02215, USA,The Center for Cancer Genome Discovery, Dana Farber Cancer Institute, Boston, MA 02215, USA
| | - Bogdan Pasaniuc
- Bioinformatics Interdepartmental Program, University of California, Los Angeles, Los Angeles, CA 90095, USA,Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, CA 90095, USA,Department of Human Genetics, David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, CA 90095, USA,Department of Computational Medicine, University of California, Los Angeles, Los Angeles, CA 90095, USA
| | - Nathan A Lack
- To whom correspondence should be addressed. Tel: +1 604 875 4411;
| |
Collapse
|
76
|
Han SK, Muto Y, Wilson PC, Humphreys BD, Sampson MG, Chakravarti A, Lee D. Quality assessment and refinement of chromatin accessibility data using a sequence-based predictive model. Proc Natl Acad Sci U S A 2022; 119:e2212810119. [PMID: 36508674 PMCID: PMC9907136 DOI: 10.1073/pnas.2212810119] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2022] [Accepted: 10/28/2022] [Indexed: 12/15/2022] Open
Abstract
Chromatin accessibility assays are central to the genome-wide identification of gene regulatory elements associated with transcriptional regulation. However, the data have highly variable quality arising from several biological and technical factors. To surmount this problem, we developed a sequence-based machine learning method to evaluate and refine chromatin accessibility data. Our framework, gapped k-mer SVM quality check (gkmQC), provides the quality metrics for a sample based on the prediction accuracy of the trained models. We tested 886 DNase-seq samples from the ENCODE/Roadmap projects to demonstrate that gkmQC can effectively identify "high-quality" (HQ) samples with low conventional quality scores owing to marginal read depths. Peaks identified in HQ samples are more accurately aligned at functional regulatory elements, show greater enrichment of regulatory elements harboring functional variants, and explain greater heritability of phenotypes from their relevant tissues. Moreover, gkmQC can optimize the peak-calling threshold to identify additional peaks, especially for rare cell types in single-cell chromatin accessibility data.
Collapse
Affiliation(s)
- Seong Kyu Han
- Department of Pediatrics, Division of Nephrology, Boston Children’s Hospital, Boston & Harvard Medical School, Boston, MA02115
- Kidney Disease Initiative, Broad Institute of MIT and Harvard, Cambridge, MA02142
| | - Yoshiharu Muto
- Division of Nephrology, Department of Medicine, Washington University in St. Louis, St. Louis, MO63130
| | - Parker C. Wilson
- Department of Pathology and Immunology, Washington University in St. Louis, St. Louis, MO63130
| | - Benjamin D. Humphreys
- Division of Nephrology, Department of Medicine, Washington University in St. Louis, St. Louis, MO63130
- Department of Developmental Biology, Washington University in St. Louis, St. Louis, MO63130
| | - Matthew G. Sampson
- Department of Pediatrics, Division of Nephrology, Boston Children’s Hospital, Boston & Harvard Medical School, Boston, MA02115
- Kidney Disease Initiative, Broad Institute of MIT and Harvard, Cambridge, MA02142
| | - Aravinda Chakravarti
- Center for Human Genetics and Genomics, New York University Grossman School of Medicine, New York, NY10016
| | - Dongwon Lee
- Department of Pediatrics, Division of Nephrology, Boston Children’s Hospital, Boston & Harvard Medical School, Boston, MA02115
- Manton Center for Orphan Disease Research, Boston Children’s Hospital, Boston, MA02115
| |
Collapse
|
77
|
Vespasiani DM, Jacobs GS, Cook LE, Brucato N, Leavesley M, Kinipi C, Ricaut FX, Cox MP, Gallego Romero I. Denisovan introgression has shaped the immune system of present-day Papuans. PLoS Genet 2022; 18:e1010470. [PMID: 36480515 PMCID: PMC9731433 DOI: 10.1371/journal.pgen.1010470] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2022] [Accepted: 10/10/2022] [Indexed: 12/13/2022] Open
Abstract
Modern humans have admixed with multiple archaic hominins. Papuans, in particular, owe up to 5% of their genome to Denisovans, a sister group to Neanderthals whose remains have only been identified in Siberia and Tibet. Unfortunately, the biological and evolutionary significance of these introgression events remain poorly understood. Here we investigate the function of both Denisovan and Neanderthal alleles characterised within a set of 56 genomes from Papuan individuals. By comparing the distribution of archaic and non-archaic variants we assess the consequences of archaic admixture across a multitude of different cell types and functional elements. We observe an enrichment of archaic alleles within cis-regulatory elements and transcribed regions of the genome, with Denisovan variants strongly affecting elements active within immune-related cells. We identify 16,048 and 10,032 high-confidence Denisovan and Neanderthal variants that fall within annotated cis-regulatory elements and with the potential to alter the affinity of multiple transcription factors to their cognate DNA motifs, highlighting a likely mechanism by which introgressed DNA can impact phenotypes. Lastly, we experimentally validate these predictions by testing the regulatory potential of five Denisovan variants segregating within Papuan individuals, and find that two are associated with a significant reduction of transcriptional activity in plasmid reporter assays. Together, these data provide support for a widespread contribution of archaic DNA in shaping the present levels of modern human genetic diversity, with different archaic ancestries potentially affecting multiple phenotypic traits within non-Africans.
Collapse
Affiliation(s)
- Davide M. Vespasiani
- Melbourne Integrative Genomics, University of Melbourne, Parkville, Australia
- School of Biosciences, University of Melbourne, Parkville, Australia
| | - Guy S. Jacobs
- Department of Archaeology, University of Cambridge, Cambridge, Uniteed Kingdom
| | - Laura E. Cook
- Melbourne Integrative Genomics, University of Melbourne, Parkville, Australia
- School of Biosciences, University of Melbourne, Parkville, Australia
| | - Nicolas Brucato
- Laboratoire de Evolution et Diversite Biologique, Université de Toulouse Midi-Pyrénées, Toulouse, France
| | - Matthew Leavesley
- School of Humanities and Social Sciences, University of Papua New Guinea, Port Moresby, Papua New Guinea
- College of Arts, Society and Education, James Cook University, Cairns, Australia
- ARC Centre of Excellence for Australian Biodiversity and Heritage, University of Wollongong, Wollongong, Australia
| | - Christopher Kinipi
- School of Humanities and Social Sciences, University of Papua New Guinea, Port Moresby, Papua New Guinea
| | - François-Xavier Ricaut
- Laboratoire de Evolution et Diversite Biologique, Université de Toulouse Midi-Pyrénées, Toulouse, France
| | - Murray P. Cox
- School of Natural Sciences, Massey University, Palmerston North, New Zealand
| | - Irene Gallego Romero
- Melbourne Integrative Genomics, University of Melbourne, Parkville, Australia
- School of Biosciences, University of Melbourne, Parkville, Australia
- Center for Stem Cell Systems, University of Melbourne, Parkville, Australia
- Center for Genomics, Evolution and Medicine, University of Tartu, Tartu, Estonia
- * E-mail:
| |
Collapse
|
78
|
Slight Variations in the Sequence Downstream of the Polyadenylation Signal Significantly Increase Transgene Expression in HEK293T and CHO Cells. Int J Mol Sci 2022; 23:ijms232415485. [PMID: 36555130 PMCID: PMC9779314 DOI: 10.3390/ijms232415485] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2022] [Revised: 11/28/2022] [Accepted: 12/01/2022] [Indexed: 12/13/2022] Open
Abstract
Compared to transcription initiation, much less is known about transcription termination. In particular, large-scale mutagenesis studies have, so far, primarily concentrated on promoter and enhancer, but not terminator sequences. Here, we used a massively parallel reporter assay (MPRA) to systematically analyze the influence of short (8 bp) sequence variants (mutations) located downstream of the polyadenylation signal (PAS) on the steady-state mRNA level of the upstream gene, employing an eGFP reporter and human HEK293T cells as a model system. In total, we evaluated 227,755 mutations located at different overlapping positions within +17..+56 bp downstream of the PAS for their ability to regulate the reporter gene expression. We found that the positions +17..+44 bp downstream of the PAS are more essential for gene upregulation than those located more distal to the PAS, and that the mutation sequences ensuring high levels of eGFP mRNA expression are extremely T-rich. Next, we validated the positive effect of a couple of mutations identified in the MPRA screening on the eGFP and luciferase protein expression. The most promising mutation increased the expression of the reporter proteins 13-fold and sevenfold on average in HEK293T and CHO cells, respectively. Overall, these findings might be useful for further improving the efficiency of production of therapeutic products, e.g., recombinant antibodies.
Collapse
|
79
|
Current challenges in understanding the role of enhancers in disease. Nat Struct Mol Biol 2022; 29:1148-1158. [PMID: 36482255 DOI: 10.1038/s41594-022-00896-3] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2022] [Accepted: 11/04/2022] [Indexed: 12/13/2022]
Abstract
Enhancers play a central role in the spatiotemporal control of gene expression and tend to work in a cell-type-specific manner. In addition, they are suggested to be major contributors to phenotypic variation, evolution and disease. There is growing evidence that enhancer dysfunction due to genetic, structural or epigenetic mechanisms contributes to a broad range of human diseases referred to as enhanceropathies. Such mechanisms often underlie the susceptibility to common diseases, but can also play a direct causal role in cancer or Mendelian diseases. Despite the recent gain of insights into enhancer biology and function, we still have a limited ability to predict how enhancer dysfunction impacts gene expression. Here we discuss the major challenges that need to be overcome when studying the role of enhancers in disease etiology and highlight opportunities and directions for future studies, aiming to disentangle the molecular basis of enhanceropathies.
Collapse
|
80
|
Koido M, Hon CC, Koyama S, Kawaji H, Murakawa Y, Ishigaki K, Ito K, Sese J, Parrish NF, Kamatani Y, Carninci P, Terao C. Prediction of the cell-type-specific transcription of non-coding RNAs from genome sequences via machine learning. Nat Biomed Eng 2022:10.1038/s41551-022-00961-8. [PMID: 36411359 DOI: 10.1038/s41551-022-00961-8] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2021] [Accepted: 10/12/2022] [Indexed: 11/22/2022]
Abstract
Gene transcription is regulated through complex mechanisms involving non-coding RNAs (ncRNAs). As the transcription of ncRNAs, especially of enhancer RNAs, is often low and cell type specific, how the levels of RNA transcription depend on genotype remains largely unexplored. Here we report the development and utility of a machine-learning model (MENTR) that reliably links genome sequence and ncRNA expression at the cell type level. Effects on ncRNA transcription predicted by the model were concordant with estimates from published studies in a cell-type-dependent manner, regardless of allele frequency and genetic linkage. Among 41,223 variants from genome-wide association studies, the model identified 7,775 enhancer RNAs and 3,548 long ncRNAs causally associated with complex traits across 348 major human primary cells and tissues, such as rare variants plausibly altering the transcription of enhancer RNAs to influence the risks of Crohn's disease and asthma. The model may aid the discovery of causal variants and the generation of testable hypotheses for biological mechanisms driving complex traits.
Collapse
Affiliation(s)
- Masaru Koido
- Laboratory for Statistical and Translational Genetics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan.,Division of Molecular Pathology, Department of Cancer Biology, Institute of Medical Science, The University of Tokyo, Tokyo, Japan.,Laboratory of Complex Trait Genomics, Department of Computational Biology and Medical Sciences, Graduate School of Frontier Sciences, The University of Tokyo, Tokyo, Japan
| | - Chung-Chau Hon
- Laboratory for Genome Information Analysis, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Satoshi Koyama
- Laboratory for Cardiovascular Genomics and Informatics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Hideya Kawaji
- Preventive Medicine and Applied Genomics Unit, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan.,Research Center for Genome & Medical Sciences, Tokyo Metropolitan Institute of Medical Science, Tokyo, Japan
| | - Yasuhiro Murakawa
- RIKEN-IFOM Joint Laboratory for Cancer Genomics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan.,IFOM ETS - The AIRC Institute of Molecular Oncology, Milan, Italy.,Institute for the Advanced Study of Human Biology, Kyoto University, Kyoto, Japan
| | - Kazuyoshi Ishigaki
- Laboratory for Statistical and Translational Genetics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan.,Divisions of Genetics and Rheumatology, Department of Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, USA.,Center for Data Sciences, Harvard Medical School, Boston, MA, USA.,Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Kaoru Ito
- Laboratory for Cardiovascular Genomics and Informatics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Jun Sese
- Artificial Intelligence Research Center, National Institute of Advanced Industrial Science and Technology, Aomi, Koto-ku, Tokyo, Japan.,Humanome Lab Inc., Tokyo, Japan
| | - Nicholas F Parrish
- Genome Immunobiology RIKEN Hakubi Research Team, RIKEN Cluster for Pioneering Research and RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Yoichiro Kamatani
- Laboratory for Statistical and Translational Genetics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan.,Laboratory of Complex Trait Genomics, Department of Computational Biology and Medical Sciences, Graduate School of Frontier Sciences, The University of Tokyo, Tokyo, Japan
| | - Piero Carninci
- Laboratory for Transcriptome Technology, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan.,Laboratory for Single Cell Technologies, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan.,Human Technopole, Milan, Italy
| | - Chikashi Terao
- Laboratory for Statistical and Translational Genetics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan. .,Clinical Research Center, Shizuoka General Hospital, Shizuoka, Japan. .,The Department of Applied Genetics, The School of Pharmaceutical Sciences, University of Shizuoka, Shizuoka, Japan.
| |
Collapse
|
81
|
Lu F, Sossin A, Abell N, Montgomery SB, He Z. Deep learning-assisted genome-wide characterization of massively parallel reporter assays. Nucleic Acids Res 2022; 50:11442-11454. [PMID: 36350674 PMCID: PMC9723615 DOI: 10.1093/nar/gkac990] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2021] [Revised: 10/04/2022] [Accepted: 10/19/2022] [Indexed: 11/10/2022] Open
Abstract
Massively parallel reporter assay (MPRA) is a high-throughput method that enables the study of the regulatory activities of tens of thousands of DNA oligonucleotides in a single experiment. While MPRA experiments have grown in popularity, their small sample sizes compared to the scale of the human genome limits our understanding of the regulatory effects they detect. To address this, we develop a deep learning model, MpraNet, to distinguish potential MPRA targets from the background genome. This model achieves high discriminative performance (AUROC = 0.85) at differentiating MPRA positives from a set of control variants that mimic the background genome when applied to the lymphoblastoid cell line. We observe that existing functional scores represent very distinct functional effects, and most of them fail to characterize the regulatory effect that MPRA detects. Using MpraNet, we predict potential MPRA functional variants across the genome and identify the distributions of MPRA effect relative to other characteristics of genetic variation, including allele frequency, alternative functional annotations specified by FAVOR, and phenome-wide associations. We also observed that the predicted MPRA positives are not uniformly distributed across the genome; instead, they are clumped together in active regions comprising 9.95% of the genome and inactive regions comprising 89.07% of the genome. Furthermore, we propose our model as a screen to filter MPRA experiment candidates at genome-wide scale, enabling future experiments to be more cost-efficient by increasing precision relative to that observed from previous MPRAs.
Collapse
Affiliation(s)
| | | | - Nathan Abell
- Department of Genetics, Stanford University, Stanford, CA 94305, USA
| | - Stephen B Montgomery
- Department of Genetics, Stanford University, Stanford, CA 94305, USA,Department of Pathology, Stanford University, Stanford, CA 94305, USA
| | - Zihuai He
- To whom correspondence should be addressed. Tel: +1 718 869 4929;
| |
Collapse
|
82
|
Cheng F, Geschwind D. Human Molecular Genetics Review Issue 2022. Hum Mol Genet 2022; 31:R1-R3. [PMID: 36268970 PMCID: PMC9989729 DOI: 10.1093/hmg/ddac219] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2022] [Revised: 08/30/2022] [Accepted: 08/30/2022] [Indexed: 11/04/2022] Open
Affiliation(s)
- Feixiong Cheng
- To whom correspondence should be addressed at: Feixiong Cheng, Ph.D., Genomic Medicine Institute, Cleveland Clinic Lerner Research Institute, 9500 Euclid Ave, Cleveland, OH 44195, USA. ; Daniel Geschwind, M.D., Ph.D., Departments of Neurology, Psychiatry and Human Genetics, and Institute of Precision Health, David Geffen School of Medicine, University of California, Los Angeles, 695 Charles E Young Dr South, Los Angeles, CA 90095.
| | - Daniel Geschwind
- To whom correspondence should be addressed at: Feixiong Cheng, Ph.D., Genomic Medicine Institute, Cleveland Clinic Lerner Research Institute, 9500 Euclid Ave, Cleveland, OH 44195, USA. ; Daniel Geschwind, M.D., Ph.D., Departments of Neurology, Psychiatry and Human Genetics, and Institute of Precision Health, David Geffen School of Medicine, University of California, Los Angeles, 695 Charles E Young Dr South, Los Angeles, CA 90095.
| |
Collapse
|
83
|
Cooper YA, Guo Q, Geschwind DH. Multiplexed functional genomic assays to decipher the noncoding genome. Hum Mol Genet 2022; 31:R84-R96. [PMID: 36057282 PMCID: PMC9585676 DOI: 10.1093/hmg/ddac194] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2022] [Revised: 08/08/2022] [Accepted: 08/09/2022] [Indexed: 11/14/2022] Open
Abstract
Linkage disequilibrium and the incomplete regulatory annotation of the noncoding genome complicates the identification of functional noncoding genetic variants and their causal association with disease. Current computational methods for variant prioritization have limited predictive value, necessitating the application of highly parallelized experimental assays to efficiently identify functional noncoding variation. Here, we summarize two distinct approaches, massively parallel reporter assays and CRISPR-based pooled screens and describe their flexible implementation to characterize human noncoding genetic variation at unprecedented scale. Each approach provides unique advantages and limitations, highlighting the importance of multimodal methodological integration. These multiplexed assays of variant effects are undoubtedly poised to play a key role in the experimental characterization of noncoding genetic risk, informing our understanding of the underlying mechanisms of disease-associated loci and the development of more robust predictive classification algorithms.
Collapse
Affiliation(s)
- Yonatan A Cooper
- Department of Human Genetics, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, CA, USA
- Medical Scientist Training Program, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, CA, USA
- Center for Neurobehavioral Genetics, Jane and Terry Semel Institute for Neuroscience and Human Behavior, University of California Los Angeles, Los Angeles, CA, USA
| | - Qiuyu Guo
- Center for Neurobehavioral Genetics, Jane and Terry Semel Institute for Neuroscience and Human Behavior, University of California Los Angeles, Los Angeles, CA, USA
| | - Daniel H Geschwind
- Department of Human Genetics, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, CA, USA
- Program in Neurogenetics, Department of Neurology, University of California Los Angeles, Los Angeles, CA, USA
- Center for Autism Research and Treatment, Semel Institute, University of California Los Angeles, Los Angeles, CA, USA
- Institute of Precision Health, University of California Los Angeles, Los Angeles, CA, USA
| |
Collapse
|
84
|
Yang H, Chen R, Wang Q, Wei Q, Ji Y, Zhong X, Li B. TVAR: assessing tissue-specific functional effects of non-coding variants with deep learning. Bioinformatics 2022; 38:4697-4704. [PMID: 36063453 PMCID: PMC9563698 DOI: 10.1093/bioinformatics/btac608] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2021] [Revised: 07/29/2022] [Accepted: 09/02/2022] [Indexed: 11/13/2022] Open
Abstract
MOTIVATION Analysis of whole-genome sequencing (WGS) for genetics is still a challenge due to the lack of accurate functional annotation of non-coding variants, especially the rare ones. As eQTLs have been extensively implicated in the genetics of human diseases, we hypothesize that rare non-coding variants discovered in WGS play a regulatory role in predisposing disease risk. RESULTS With thousands of tissue- and cell-type-specific epigenomic features, we propose TVAR. This multi-label learning-based deep neural network predicts the functionality of non-coding variants in the genome based on eQTLs across 49 human tissues in the GTEx project. TVAR learns the relationships between high-dimensional epigenomics and eQTLs across tissues, taking the correlation among tissues into account to understand shared and tissue-specific eQTL effects. As a result, TVAR outputs tissue-specific annotations, with an average AUROC of 0.77 across these tissues. We evaluate TVAR's performance on four complex diseases (coronary artery disease, breast cancer, Type 2 diabetes and Schizophrenia), using TVAR's tissue-specific annotations, and observe its superior performance in predicting functional variants for both common and rare variants, compared with five existing state-of-the-art tools. We further evaluate TVAR's G-score, a scoring scheme across all tissues, on ClinVar, fine-mapped GWAS loci, Massive Parallel Reporter Assay (MPRA) validated variants and observe the consistently better performance of TVAR compared with other competing tools. AVAILABILITY AND IMPLEMENTATION The TVAR source code and its scores on the ClinVar catalog, fine mapped GWAS Loci, high confidence eQTLs from GTEx dataset, and MPRA validated functional variants are available at https://github.com/haiyang1986/TVAR. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Hai Yang
- Department of Computer Science and Engineering, East China University of Science and Technology, Shanghai 200237, China
- Department of Molecular Physiology & Biophysics, Vanderbilt University, Nashville, TN 37232, USA
| | - Rui Chen
- Department of Molecular Physiology & Biophysics, Vanderbilt University, Nashville, TN 37232, USA
- Vanderbilt Genetics Institute, Vanderbilt University, Nashville, TN 37232, USA
| | - Quan Wang
- Department of Molecular Physiology & Biophysics, Vanderbilt University, Nashville, TN 37232, USA
- Vanderbilt Genetics Institute, Vanderbilt University, Nashville, TN 37232, USA
| | - Qiang Wei
- Department of Molecular Physiology & Biophysics, Vanderbilt University, Nashville, TN 37232, USA
- Vanderbilt Genetics Institute, Vanderbilt University, Nashville, TN 37232, USA
| | - Ying Ji
- Department of Molecular Physiology & Biophysics, Vanderbilt University, Nashville, TN 37232, USA
- Vanderbilt Genetics Institute, Vanderbilt University, Nashville, TN 37232, USA
| | - Xue Zhong
- Vanderbilt Genetics Institute, Vanderbilt University, Nashville, TN 37232, USA
- Department of Medicine, Vanderbilt University Medical Center, Nashville, TN 37232, USA
| | - Bingshan Li
- Department of Molecular Physiology & Biophysics, Vanderbilt University, Nashville, TN 37232, USA
- Vanderbilt Genetics Institute, Vanderbilt University, Nashville, TN 37232, USA
| |
Collapse
|
85
|
Nuytemans K, Vasquez ML, Wang L, Van Booven D, Griswold AJ, Rajabli F, Celis K, Oron O, Hofmann N, Rolati S, Garcia-Serje C, Zhang S, Jin F, Argenziano M, Grant SF, Chesi A, Brown CD, Young JI, Dykxhoorn DM, Pericak-Vance MA, Vance JM. Identifying differential regulatory control of APOE ɛ4 on African versus European haplotypes as potential therapeutic targets. Alzheimers Dement 2022; 18:1930-1942. [PMID: 34978147 PMCID: PMC9250552 DOI: 10.1002/alz.12534] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2021] [Revised: 09/29/2021] [Accepted: 10/19/2021] [Indexed: 01/28/2023]
Abstract
We previously demonstrated that in Alzheimer's disease (AD) patients, European apolipoprotein E (APOE) ε4 carriers express significantly more APOE ε4 in their brains than African AD carriers. We examined single nucleotide polymorphisms near APOE with significant frequency differences between African and European/Japanese APOE ε4 haplotypes that could contribute to this difference in expression through regulation. Two enhancer massively parallel reporter assay (MPRA) approaches were performed, supplemented with single fragment reporter assays. We used Capture C analyses to support interactions with the APOE promoter. Introns within TOMM40 showed increased enhancer activity in the European/Japanese versus African haplotypes in astrocytes and microglia. This region overlaps with APOE promoter interactions as assessed by Capture C analysis. Single variant analyses pinpoints rs2075650/rs157581, and rs59007384 as functionally different on these haplotypes. Identification of the mechanisms for differential regulatory function for APOE expression between African and European/Japanese haplotypes could lead to therapeutic targets for APOE ε4 carriers.
Collapse
Affiliation(s)
- Karen Nuytemans
- John P. Hussman Institute for Human Genomics, University of
Miami Miller School of Medicine, Miami, Florida, USA
- John T. Macdonald Foundation Department of Human Genetics,
University of Miami Miller School of Medicine, Miami, Florida, USA
| | - Marina Lipkin Vasquez
- John P. Hussman Institute for Human Genomics, University of
Miami Miller School of Medicine, Miami, Florida, USA
| | - Liyong Wang
- John P. Hussman Institute for Human Genomics, University of
Miami Miller School of Medicine, Miami, Florida, USA
- John T. Macdonald Foundation Department of Human Genetics,
University of Miami Miller School of Medicine, Miami, Florida, USA
| | - Derek Van Booven
- John P. Hussman Institute for Human Genomics, University of
Miami Miller School of Medicine, Miami, Florida, USA
| | - Anthony J. Griswold
- John P. Hussman Institute for Human Genomics, University of
Miami Miller School of Medicine, Miami, Florida, USA
- John T. Macdonald Foundation Department of Human Genetics,
University of Miami Miller School of Medicine, Miami, Florida, USA
| | - Farid Rajabli
- John P. Hussman Institute for Human Genomics, University of
Miami Miller School of Medicine, Miami, Florida, USA
- John T. Macdonald Foundation Department of Human Genetics,
University of Miami Miller School of Medicine, Miami, Florida, USA
| | - Katrina Celis
- John P. Hussman Institute for Human Genomics, University of
Miami Miller School of Medicine, Miami, Florida, USA
| | - Oded Oron
- John P. Hussman Institute for Human Genomics, University of
Miami Miller School of Medicine, Miami, Florida, USA
| | - Natalia Hofmann
- John P. Hussman Institute for Human Genomics, University of
Miami Miller School of Medicine, Miami, Florida, USA
| | - Sophie Rolati
- John P. Hussman Institute for Human Genomics, University of
Miami Miller School of Medicine, Miami, Florida, USA
| | - Catherine Garcia-Serje
- John P. Hussman Institute for Human Genomics, University of
Miami Miller School of Medicine, Miami, Florida, USA
| | - Shanshan Zhang
- Department of Genetics and Genome Sciences, School of
Medicine, Case Western Reserve University, Cleveland, Ohio, USA
- The Biomedical Sciences Training Program (BSTP), School of
Medicine, Case Western Reserve University, Cleveland, Ohio, USA
| | - Fulai Jin
- Department of Genetics and Genome Sciences, School of
Medicine, Case Western Reserve University, Cleveland, Ohio, USA
- Department of Computer and Data Sciences, Case
Comprehensive Cancer Center Case Western Reserve University, Cleveland, Ohio,
USA
| | | | - Struan F.A. Grant
- Division of Human Genetics, Children’s Hospital of
Philadelphia, Philadelphia, Pennsylvania, USA
- Department of Genetics, University of Pennsylvania,
Philadelphia, Pennsylvania, USA
- Center for Spatial and Functional Genomics,
Children’s Hospital of Philadelphia, Philadelphia, Pennsylvania, USA
- Division of Endocrinology and Diabetes, Children’s
Hospital of Philadelphia, Philadelphia, Pennsylvania, USA
- Department of Pediatrics, Perelman School of Medicine,
University of Pennsylvania, Philadelphia, Pennsylvania, USA
| | - Alessandra Chesi
- Division of Human Genetics, Children’s Hospital of
Philadelphia, Philadelphia, Pennsylvania, USA
- Department of Genetics, University of Pennsylvania,
Philadelphia, Pennsylvania, USA
- Department of Pathology and Laboratory Medicine,
University of PennsylvaniaPerelman School of Medicine, Philadelphia, Pennsylvania,
USA
| | - Christopher D. Brown
- Department of Genetics, University of
PennsylvaniaPerelman School of Medicine, Philadelphia, Pennsylvania, USA
| | - Juan I. Young
- John P. Hussman Institute for Human Genomics, University of
Miami Miller School of Medicine, Miami, Florida, USA
- John T. Macdonald Foundation Department of Human Genetics,
University of Miami Miller School of Medicine, Miami, Florida, USA
| | - Derek M. Dykxhoorn
- John P. Hussman Institute for Human Genomics, University of
Miami Miller School of Medicine, Miami, Florida, USA
- John T. Macdonald Foundation Department of Human Genetics,
University of Miami Miller School of Medicine, Miami, Florida, USA
| | - Margaret A. Pericak-Vance
- John P. Hussman Institute for Human Genomics, University of
Miami Miller School of Medicine, Miami, Florida, USA
- John T. Macdonald Foundation Department of Human Genetics,
University of Miami Miller School of Medicine, Miami, Florida, USA
| | - Jeffery M. Vance
- John P. Hussman Institute for Human Genomics, University of
Miami Miller School of Medicine, Miami, Florida, USA
- John T. Macdonald Foundation Department of Human Genetics,
University of Miami Miller School of Medicine, Miami, Florida, USA
| |
Collapse
|
86
|
Crinzi EA, Haley EK, Poppenberg KE, Jiang K, Tutino VM, Jarvis JN. Analysis of chromatin data supports a role for CD14+ monocytes/macrophages in mediating genetic risk for juvenile idiopathic arthritis. Front Immunol 2022; 13:913555. [PMID: 36248892 PMCID: PMC9559786 DOI: 10.3389/fimmu.2022.913555] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2022] [Accepted: 09/06/2022] [Indexed: 11/13/2022] Open
Abstract
IntroductionGenome wide association studies (GWAS) have identified multiple regions that confer genetic risk for the polyarticular/oligoarticular forms of juvenile idiopathic arthritis (JIA). However, genome-wide scans do not identify the cells impacted by genetic polymorphisms on the risk haplotypes or the genes impacted by those variants. We have shown that genetic variants driving JIA risk are likely to affect both innate and adaptive immune functions. We provide additional evidence that JIA risk variants impact innate immunity.Materials and methodsWe queried publicly available H3K4me1/H3K27ac ChIP-seq data in CD14+ monocytes to determine whether the linkage disequilibrium (LD) blocks incorporating the SNPs that tag JIA risk loci showed enrichment for these epigenetic marks. We also queried monocyte/macrophage GROseq data, a functional readout of active enhancers. We defined the topologically associated domains (TADs) encompassing enhancers on the risk haplotypes and identified genes within those TADs expressed in monocytes. We performed ontology analyses of these genes to identify cellular processes that may be impacted by these variants. We also used whole blood RNAseq data from the Genotype-Tissue Expression (GTEx) data base to determine whether SNPs lying within monocyte GROseq peaks influence plausible target genes within the TADs encompassing the JIA risk haplotypes.ResultsThe LD blocks encompassing the JIA genetic risk regions were enriched for H3K4me1/H3K27ac ChIPseq peaks (p=0.00021 and p=0.022) when compared to genome background. Eleven and sixteen JIA were enriched for resting and activated macrophage GROseq peaks, respectively risk regions (p=0.04385 and p=0.00004). We identified 321 expressed genes within the TADs encompassing the JIA haplotypes in human monocytes. Ontological analysis of these genes showed enrichment for multiple immune functions. Finally, we found that SNPs lying within the GROseq peaks are strongly associated with expression levels of plausible target genes in human whole blood.ConclusionsThese findings support the idea that both innate and adaptive immunity are impacted by JIA genetic risk variants.
Collapse
Affiliation(s)
- Elizabeth A. Crinzi
- Department of Pediatrics, University at Buffalo Jacobs School of Medicine & Biomedical Sciences, Buffalo, NY, United States
| | - Emma K. Haley
- Department of Pediatrics, University at Buffalo Jacobs School of Medicine & Biomedical Sciences, Buffalo, NY, United States
| | - Kerry E. Poppenberg
- Canon Stroke and Vascular Center, University at Buffalo Jacobs School of Medicine & Biomedical Sciences, Buffalo, NY, United States
- Department of Neurosurgery, University at Buffalo Jacobs School of Medicine & Biomedical Sciences, Buffalo, NY, United States
| | - Kaiyu Jiang
- Department of Pediatrics, University at Buffalo Jacobs School of Medicine & Biomedical Sciences, Buffalo, NY, United States
| | - Vincent M. Tutino
- Canon Stroke and Vascular Center, University at Buffalo Jacobs School of Medicine & Biomedical Sciences, Buffalo, NY, United States
- Department of Neurosurgery, University at Buffalo Jacobs School of Medicine & Biomedical Sciences, Buffalo, NY, United States
| | - James N. Jarvis
- Department of Pediatrics, University at Buffalo Jacobs School of Medicine & Biomedical Sciences, Buffalo, NY, United States
- Genetics, Genomics, & Bioinformatics Program, University at Buffalo Jacobs School of Medicine & Biomedical Sciences, Buffalo, NY, United States
- *Correspondence: James N. Jarvis,
| |
Collapse
|
87
|
Cooper YA, Teyssier N, Dräger NM, Guo Q, Davis JE, Sattler SM, Yang Z, Patel A, Wu S, Kosuri S, Coppola G, Kampmann M, Geschwind DH. Functional regulatory variants implicate distinct transcriptional networks in dementia. Science 2022; 377:eabi8654. [PMID: 35981026 DOI: 10.1126/science.abi8654] [Citation(s) in RCA: 48] [Impact Index Per Article: 24.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Predicting the function of noncoding variation is a major challenge in modern genetics. In this study, we used massively parallel reporter assays to screen 5706 variants identified from genome-wide association studies for both Alzheimer's disease (AD) and progressive supranuclear palsy (PSP), identifying 320 functional regulatory variants (frVars) across 27 loci, including the complex 17q21.31 region. We identified and validated multiple risk loci using CRISPR interference or excision, including complement 4 (C4A) and APOC1 in AD and PLEKHM1 and KANSL1 in PSP. Functional variants disrupt transcription factor binding sites converging on enhancers with cell type-specific activity in PSP and AD, implicating a neuronal SP1-driven regulatory network in PSP pathogenesis. These analyses suggest that noncoding genetic risk is driven by common genetic variants through their aggregate activity on specific transcriptional programs.
Collapse
Affiliation(s)
- Yonatan A Cooper
- Department of Human Genetics, David Geffen School of Medicine, University of California, Los Angeles, CA 90095, USA
- Medical Scientist Training Program, David Geffen School of Medicine, University of California, Los Angeles, CA 90095, USA
- Center for Neurobehavioral Genetics, Jane and Terry Semel Institute for Neuroscience and Human Behavior, University of California, Los Angeles, CA 90095, USA
| | - Noam Teyssier
- Institute for Neurodegenerative Diseases, University of California, San Francisco, CA 94158, USA
| | - Nina M Dräger
- Institute for Neurodegenerative Diseases, University of California, San Francisco, CA 94158, USA
| | - Qiuyu Guo
- Center for Neurobehavioral Genetics, Jane and Terry Semel Institute for Neuroscience and Human Behavior, University of California, Los Angeles, CA 90095, USA
| | - Jessica E Davis
- Department of Chemistry and Biochemistry, University of California, Los Angeles, CA 90095, USA
| | - Sydney M Sattler
- Institute for Neurodegenerative Diseases, University of California, San Francisco, CA 94158, USA
| | - Zhongan Yang
- Center for Neurobehavioral Genetics, Jane and Terry Semel Institute for Neuroscience and Human Behavior, University of California, Los Angeles, CA 90095, USA
| | - Abdulsamie Patel
- Center for Neurobehavioral Genetics, Jane and Terry Semel Institute for Neuroscience and Human Behavior, University of California, Los Angeles, CA 90095, USA
| | - Sarah Wu
- Center for Neurobehavioral Genetics, Jane and Terry Semel Institute for Neuroscience and Human Behavior, University of California, Los Angeles, CA 90095, USA
| | - Sriram Kosuri
- Department of Chemistry and Biochemistry, University of California, Los Angeles, CA 90095, USA
| | - Giovanni Coppola
- Center for Neurobehavioral Genetics, Jane and Terry Semel Institute for Neuroscience and Human Behavior, University of California, Los Angeles, CA 90095, USA
- Department of Psychiatry and Biobehavioral Sciences, David Geffen School of Medicine, University of California, Los Angeles, CA 90095, USA
| | - Martin Kampmann
- Institute for Neurodegenerative Diseases, University of California, San Francisco, CA 94158, USA
- Department of Biochemistry and Biophysics, University of California, San Francisco, CA 94143, USA
| | - Daniel H Geschwind
- Department of Human Genetics, David Geffen School of Medicine, University of California, Los Angeles, CA 90095, USA
- Program in Neurogenetics, Department of Neurology, University of California, Los Angeles, CA 90095, USA
- Center for Autism Research and Treatment, Jane and Terry Semel Institute for Neuroscience and Human Behavior, University of California, Los Angeles, CA 90095, USA
- Institute of Precision Health, University of California, Los Angeles, CA 90095, USA
| |
Collapse
|
88
|
Schipper M, Posthuma D. "Demystifying non-coding GWAS variants: an overview of computational tools and methods.". Hum Mol Genet 2022; 31:R73-R83. [PMID: 35972862 PMCID: PMC9585674 DOI: 10.1093/hmg/ddac198] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2022] [Revised: 08/11/2022] [Accepted: 08/11/2022] [Indexed: 02/01/2023] Open
Abstract
Genome-wide association studies (GWAS) have found the majority of disease-associated variants to be non-coding. Major efforts into the charting of the non-coding regulatory landscapes have allowed for the development of tools and methods which aim to aid in the identification of causal variants and their mechanism of action. In this review, we give an overview of current tools and methods for the analysis of non-coding GWAS variants in disease. We provide a workflow that allows for the accumulation of in silico evidence to generate novel hypotheses on mechanisms underlying disease and prioritize targets for follow-up study using non-coding GWAS variants. Lastly, we discuss the need for comprehensive benchmarks and novel tools for the analysis of non-coding variants.
Collapse
Affiliation(s)
- Marijn Schipper
- Department of Complex Trait Genetics, Center for Neurogenomics and Cognitive Research, Amsterdam Neuroscience, VU University Amsterdam, De Boelelaan 1105 1081HV Amsterdam, The Netherlands
| | - Danielle Posthuma
- Department of Complex Trait Genetics, Center for Neurogenomics and Cognitive Research, Amsterdam Neuroscience, VU University Amsterdam, De Boelelaan 1105 1081HV Amsterdam, The Netherlands
| |
Collapse
|
89
|
Mondal S, Ramanathan M, Miao W, Meyers RM, Rao D, Lopez-Pajares V, Siprashvili Z, Reynolds DL, Porter DF, Ferguson I, Neela P, Zhao Y, Meservey LM, Guo M, Yang YY, Li L, Wang Y, Khavari PA. PROBER identifies proteins associated with programmable sequence-specific DNA in living cells. Nat Methods 2022; 19:959-968. [PMID: 35927480 PMCID: PMC10202087 DOI: 10.1038/s41592-022-01552-w] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2021] [Accepted: 06/20/2022] [Indexed: 11/08/2022]
Abstract
DNA-protein interactions mediate physiologic gene regulation and may be altered by DNA variants linked to polygenic disease. To enhance the speed and signal-to-noise ratio (SNR) in the identification and quantification of proteins associated with specific DNA sequences in living cells, we developed proximal biotinylation by episomal recruitment (PROBER). PROBER uses high-copy episomes to amplify SNR, and proximity proteomics (BioID) to identify the transcription factors and additional gene regulators associated with short DNA sequences of interest. PROBER quantified both constitutive and inducible association of transcription factors and corresponding chromatin regulators to target DNA sequences and binding quantitative trait loci due to single-nucleotide variants. PROBER identified alterations in regulator associations due to cancer hotspot mutations in the hTERT promoter, indicating that these mutations increase promoter association with specific gene activators. PROBER provides an approach to rapidly identify proteins associated with specific DNA sequences and their variants in living cells.
Collapse
Affiliation(s)
- Smarajit Mondal
- Program in Epithelial Biology, Stanford University, Stanford, CA, USA
| | | | - Weili Miao
- Program in Epithelial Biology, Stanford University, Stanford, CA, USA
| | - Robin M Meyers
- Program in Epithelial Biology, Stanford University, Stanford, CA, USA
| | - Deepti Rao
- Program in Epithelial Biology, Stanford University, Stanford, CA, USA
| | | | - Zurab Siprashvili
- Program in Epithelial Biology, Stanford University, Stanford, CA, USA
| | - David L Reynolds
- Program in Epithelial Biology, Stanford University, Stanford, CA, USA
| | - Douglas F Porter
- Program in Epithelial Biology, Stanford University, Stanford, CA, USA
| | - Ian Ferguson
- Program in Epithelial Biology, Stanford University, Stanford, CA, USA
| | - Poornima Neela
- Program in Epithelial Biology, Stanford University, Stanford, CA, USA
| | - Yang Zhao
- Program in Epithelial Biology, Stanford University, Stanford, CA, USA
| | | | - Margaret Guo
- Program in Epithelial Biology, Stanford University, Stanford, CA, USA
- Program in Biomedical Informatics, Stanford University, Stanford, CA, USA
| | - Yen-Yu Yang
- Department of Chemistry, University of California, Riverside, CA, USA
| | - Lin Li
- Department of Chemistry, University of California, Riverside, CA, USA
| | - Yinsheng Wang
- Department of Chemistry, University of California, Riverside, CA, USA
| | - Paul A Khavari
- Program in Epithelial Biology, Stanford University, Stanford, CA, USA.
- Stanford Cancer Institute, Stanford University, Stanford, CA, USA.
- Veterans Affairs, Palo Alto Healthcare System, Palo Alto, CA, USA.
| |
Collapse
|
90
|
McQuerry JA, Mclaird M, Hartin SN, Means JC, Johnston J, Pastinen T, Younger ST. Massively parallel identification of functionally consequential noncoding genetic variants in undiagnosed rare disease patients. Sci Rep 2022; 12:7576. [PMID: 35534523 PMCID: PMC9085742 DOI: 10.1038/s41598-022-11589-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2021] [Accepted: 04/26/2022] [Indexed: 11/16/2022] Open
Abstract
Clinical whole genome sequencing has enabled the discovery of potentially pathogenic noncoding variants in the genomes of rare disease patients with a prior history of negative genetic testing. However, interpreting the functional consequences of noncoding variants and distinguishing those that contribute to disease etiology remains a challenge. Here we address this challenge by experimentally profiling the functional consequences of rare noncoding variants detected in a cohort of undiagnosed rare disease patients at scale using a massively parallel reporter assay. We demonstrate that this approach successfully identifies rare noncoding variants that alter the regulatory capacity of genomic sequences. In addition, we describe an integrative analysis that utilizes genomic features alongside patient clinical data to further prioritize candidate variants with an increased likelihood of pathogenicity. This work represents an important step towards establishing a framework for the functional interpretation of clinically detected noncoding variants.
Collapse
Affiliation(s)
- Jasmine A McQuerry
- Genomic Medicine Center, Children's Mercy Kansas City, Kansas City, MO, 64108, USA
- Children's Mercy Research Institute, Children's Mercy Kansas City, Kansas City, MO, 64108, USA
| | - Merry Mclaird
- Genomic Medicine Center, Children's Mercy Kansas City, Kansas City, MO, 64108, USA
- Children's Mercy Research Institute, Children's Mercy Kansas City, Kansas City, MO, 64108, USA
| | - Samantha N Hartin
- Genomic Medicine Center, Children's Mercy Kansas City, Kansas City, MO, 64108, USA
- Children's Mercy Research Institute, Children's Mercy Kansas City, Kansas City, MO, 64108, USA
| | - John C Means
- Genomic Medicine Center, Children's Mercy Kansas City, Kansas City, MO, 64108, USA
- Children's Mercy Research Institute, Children's Mercy Kansas City, Kansas City, MO, 64108, USA
| | - Jeffrey Johnston
- Genomic Medicine Center, Children's Mercy Kansas City, Kansas City, MO, 64108, USA
- Children's Mercy Research Institute, Children's Mercy Kansas City, Kansas City, MO, 64108, USA
| | - Tomi Pastinen
- Genomic Medicine Center, Children's Mercy Kansas City, Kansas City, MO, 64108, USA
- Children's Mercy Research Institute, Children's Mercy Kansas City, Kansas City, MO, 64108, USA
- Department of Pediatrics, University of Missouri-Kansas City School of Medicine, Kansas City, MO, 64110, USA
- Department of Cancer Biology, University of Kansas Medical Center, Kansas City, KS, 66160, USA
| | - Scott T Younger
- Genomic Medicine Center, Children's Mercy Kansas City, Kansas City, MO, 64108, USA.
- Children's Mercy Research Institute, Children's Mercy Kansas City, Kansas City, MO, 64108, USA.
- Department of Pediatrics, University of Missouri-Kansas City School of Medicine, Kansas City, MO, 64110, USA.
- Department of Pediatrics, University of Kansas Medical Center, Kansas City, KS, 66160, USA.
| |
Collapse
|
91
|
Liu Z, Hun J, Wang Y, Yang M, Niu L, Shi B. Association of serum C1Q/TNF-related protein 4 levels with carotid atherosclerosis in subjects with type 2 diabetes mellitus: a cross-sectional study. Clin Chim Acta 2022; 531:337-341. [PMID: 35525266 DOI: 10.1016/j.cca.2022.04.1004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2022] [Revised: 04/12/2022] [Accepted: 04/28/2022] [Indexed: 11/03/2022]
Abstract
BACKGROUND Carotid atherosclerosis (CAS) is a common manifestation of macroangiopathy in type 2 diabetes mellitus (T2DM). C1Q/TNF-related protein 4 (CTRP4) was found to be involved in regulation of food intake behaviors and glucolipid metabolism, which were also key factors in the development of CAS. However, the relationship between serum CTRP4 and CAS in T2DM remains unclear. METHODS A total of 111 participants with T2DM were enrolled in the study and were divided into 2 groups (T2DM group and T2DM + CAS group) according to the result of carotid ultrasound examinations. Serum CTRP4 levels were measured by enzyme linked immunosorbent assay (ELISA). Trend χ2 test and binary stepwise logistic regression were conducted to assess the association between serum CTRP4 and the risk of CAS in T2DM. RESULTS Serum CTRP4 concentrations in T2DM + CAS group were significantly lower compared with those in T2DM group [7.98 (5.53) vs. 11.29 (7.36) ng/ml, P < 0.01]. The risk of CAS in T2DM decreased with the increasing of CTRP4 quartiles (P for trend < 0.01). Binary stepwise logistic regression suggested that serum CTRP4 might be an independent influence factor for CAS in patients with T2DM (P < 0.01) and high concentrations of serum CTRP4 were related to low risk of CAS in T2DM. CONCLUSIONS The concentrations of serum CTRP4 are lower in T2DM patients with CAS compared to those without CAS. Serum CTRP4 levels are negatively related to the risk of CAS in T2DM.
Collapse
Affiliation(s)
- Zheng Liu
- Department of Endocrinology, First Affiliated Hospital University, Suzhou, Jiangsu, China
| | - Junxia Hun
- Department of Endocrinology, Renmin Hospital of Wuhan University, Wuhan, Hubeil, China
| | - Ye Wang
- Department of Endocrinology, First Affiliated Hospital University, Suzhou, Jiangsu, China
| | - Minyi Yang
- Department of Endocrinology, First Affiliated Hospital University, Suzhou, Jiangsu, China
| | - Lijuan Niu
- Department of Endocrinology, First Affiliated Hospital University, Suzhou, Jiangsu, China.
| | - Bimin Shi
- Department of Endocrinology, First Affiliated Hospital University, Suzhou, Jiangsu, China
| |
Collapse
|
92
|
Mouri K, Guo MH, de Boer CG, Lissner MM, Harten IA, Newby GA, DeBerg HA, Platt WF, Gentili M, Liu DR, Campbell DJ, Hacohen N, Tewhey R, Ray JP. Prioritization of autoimmune disease-associated genetic variants that perturb regulatory element activity in T cells. Nat Genet 2022; 54:603-612. [PMID: 35513721 PMCID: PMC9793778 DOI: 10.1038/s41588-022-01056-5] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2021] [Accepted: 03/18/2022] [Indexed: 12/30/2022]
Abstract
Genome-wide association studies (GWASs) have uncovered hundreds of autoimmune disease-associated loci; however, the causal genetic variants within each locus are mostly unknown. Here, we perform high-throughput allele-specific reporter assays to prioritize disease-associated variants for five autoimmune diseases. By examining variants that both promote allele-specific reporter expression and are located in accessible chromatin, we identify 60 putatively causal variants that enrich for statistically fine-mapped variants by up to 57.8-fold. We introduced the risk allele of a prioritized variant (rs72928038) into a human T cell line and deleted the orthologous sequence in mice, both resulting in reduced BACH2 expression. Naive CD8 T cells from mice containing the deletion had reduced expression of genes that suppress activation and maintain stemness and, upon acute viral infection, displayed greater propensity to become effector T cells. Our results represent an example of an effective approach for prioritizing variants and studying their physiologically relevant effects.
Collapse
Affiliation(s)
| | - Michael H Guo
- Department of Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Carl G de Boer
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
- School of Biomedical Engineering, University of British Columbia, Vancouver, BC, Canada
| | - Michelle M Lissner
- Fundamental Immunology, Benaroya Research Institute, Seattle, WA, USA
- Division of Rheumatology, Department of Medicine, University of Washington, Seattle, WA, USA
| | - Ingrid A Harten
- Systems Immunology, Benaroya Research Institute, Seattle, WA, USA
| | - Gregory A Newby
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Merkin Institute of Transformative Technologies in Healthcare, Broad Institute of Harvard and MIT, Cambridge, MA, USA
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA, USA
- Howard Hughes Medical Institute, Harvard University, Cambridge, MA, USA
| | - Hannah A DeBerg
- Systems Immunology, Benaroya Research Institute, Seattle, WA, USA
| | - Winona F Platt
- Systems Immunology, Benaroya Research Institute, Seattle, WA, USA
| | | | - David R Liu
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Merkin Institute of Transformative Technologies in Healthcare, Broad Institute of Harvard and MIT, Cambridge, MA, USA
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA, USA
- Howard Hughes Medical Institute, Harvard University, Cambridge, MA, USA
| | - Daniel J Campbell
- Fundamental Immunology, Benaroya Research Institute, Seattle, WA, USA
- Department of Immunology, University of Washington, Seattle, WA, USA
| | - Nir Hacohen
- Broad Institute of MIT and Harvard, Cambridge, MA, USA.
- Center for Cancer Research, Massachusetts General Hospital, Boston, MA, USA.
| | - Ryan Tewhey
- The Jackson Laboratory, Bar Harbor, ME, USA.
| | - John P Ray
- Broad Institute of MIT and Harvard, Cambridge, MA, USA.
- Systems Immunology, Benaroya Research Institute, Seattle, WA, USA.
- Department of Immunology, University of Washington, Seattle, WA, USA.
| |
Collapse
|
93
|
Wu T, Li G. An Improved EMSA-based Method to Prioritize Candidate cis-REs for Further Functional Validation. Bio Protoc 2022; 12:e4397. [PMID: 35800102 PMCID: PMC9081477 DOI: 10.21769/bioprotoc.4397] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Revised: 01/19/2022] [Accepted: 03/14/2022] [Indexed: 12/29/2022] Open
Abstract
Cells are the complex product of gene expression programs that involve the coordinated transcription of thousands of genes controlled by cis-regulatory elements (cis-REs). Therefore, identification of cis-REs is the key to decipher the mechanisms underlying the regulation of gene expression. Here, we describe a simple and time-effective protocol of fine-mapping cis-REs by using an electrophoresis mobility shift assay (EMSA)-based, in vitro, high-throughput (HTP) technique called regulatory element-sequencing (Reel-seq). Reel-seq can be applied to identify cis-REs at a high resolution and sensitivity over large genome regions, in a systematic and continuous manner. It can be used to prioritize candidate cis-REs as a complement to the existing approaches, such as massive parallel reporter assay (MPRA), chromatin immunoprecipitation DNA-sequencing (ChIP-seq), and the assay for transposase-accessible chromatin sequencing (ATAC-seq). Graphical abstract: Generation of the Reel-seq Library 1 and 2 (A) and identification of cis-REs by an electrophoresis mobility shift assay (EMSA)-based Reel-seq screen (B). NE: nuclear extract; NGS: next generation sequencing.
Collapse
Affiliation(s)
- Ting Wu
- Aging Institute, University of Pittsburgh, Pittsburgh, PA 15219, USA
- Department of Medicine, Xiangya School of Medicine, Central South University, Changsha, 410083, China
| | - Gang Li
- Aging Institute, University of Pittsburgh, Pittsburgh, PA 15219, USA
- Department of Medicine, Division of Cardiology, University of Pittsburgh Medical Center, Pittsburgh, PA 15223, USA
| |
Collapse
|
94
|
Redefining the hypotheses driving Parkinson's diseases research. NPJ Parkinsons Dis 2022; 8:45. [PMID: 35440633 PMCID: PMC9018840 DOI: 10.1038/s41531-022-00307-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2021] [Accepted: 03/04/2022] [Indexed: 12/20/2022] Open
Abstract
Parkinson’s disease (PD) research has largely focused on the disease as a single entity centred on the development of neuronal pathology within the central nervous system. However, there is growing recognition that PD is not a single entity but instead reflects multiple diseases, in which different combinations of environmental, genetic and potential comorbid factors interact to direct individual disease trajectories. Moreover, an increasing body of recent research implicates peripheral tissues and non-neuronal cell types in the development of PD. These observations are consistent with the hypothesis that the initial causative changes for PD development need not occur in the central nervous system. Here, we discuss how the use of neuronal pathology as a shared, qualitative phenotype minimises insights into the possibility of multiple origins and aetiologies of PD. Furthermore, we discuss how considering PD as a single entity potentially impairs our understanding of the causative molecular mechanisms, approaches for patient stratification, identification of biomarkers, and the development of therapeutic approaches to PD. The clear consequence of there being distinct diseases that collectively form PD, is that there is no single biomarker or treatment for PD development or progression. We propose that diagnosis should shift away from the clinical definitions, towards biologically defined diseases that collectively form PD, to enable informative patient stratification. N-of-one type, clinical designs offer an unbiased, and agnostic approach to re-defining PD in terms of a group of many individual diseases.
Collapse
|
95
|
Abdalla M, Abdalla M. A general framework for predicting the transcriptomic consequences of non-coding variation and small molecules. PLoS Comput Biol 2022; 18:e1010028. [PMID: 35421087 PMCID: PMC9041867 DOI: 10.1371/journal.pcbi.1010028] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2021] [Revised: 04/26/2022] [Accepted: 03/16/2022] [Indexed: 11/18/2022] Open
Abstract
Genome wide association studies (GWASs) for complex traits have implicated thousands of genetic loci. Most GWAS-nominated variants lie in noncoding regions, complicating the systematic translation of these findings into functional understanding. Here, we leverage convolutional neural networks to assist in this challenge. Our computational framework, peaBrain, models the transcriptional machinery of a tissue as a two-stage process: first, predicting the mean tissue specific abundance of all genes and second, incorporating the transcriptomic consequences of genotype variation to predict individual abundance on a subject-by-subject basis. We demonstrate that peaBrain accounts for the majority (>50%) of variance observed in mean transcript abundance across most tissues and outperforms regularized linear models in predicting the consequences of individual genotype variation. We highlight the validity of the peaBrain model by calculating non-coding impact scores that correlate with nucleotide evolutionary constraint that are also predictive of disease-associated variation and allele-specific transcription factor binding. We further show how these tissue-specific peaBrain scores can be leveraged to pinpoint functional tissues underlying complex traits, outperforming methods that depend on colocalization of eQTL and GWAS signals. We subsequently: (a) derive continuous dense embeddings of genes for downstream applications; (b) highlight the utility of the model in predicting transcriptomic impact of small molecules and shRNA (on par with in vitro experimental replication of external test sets); (c) explore how peaBrain can be used to model difficult-to-study processes (such as neural induction); and (d) identify putatively functional eQTLs that are missed by high-throughput experimental approaches.
Collapse
Affiliation(s)
- Moustafa Abdalla
- Wellcome Trust Centre for Human Genetics, Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom
- Oxford Centre for Diabetes, Endocrinology and Metabolism, Radcliffe Department of Medicine, University of Oxford, Oxford, United Kingdom
- Computational Statistics and Machine Learning, Department of Statistics, University of Oxford, Oxford, United Kingdom
- Department of Surgery, Harvard Medical School, Boston, Massachusetts, United States of America
- * E-mail: (MA); (MA)
| | - Mohamed Abdalla
- Vector Institute for Artificial Intelligence, Toronto, Canada
- Department of Computer Science, University of Toronto, Toronto, Canada
- * E-mail: (MA); (MA)
| |
Collapse
|
96
|
Chen L, Wang Y, Zhao F. Exploiting deep transfer learning for the prediction of functional non-coding variants using genomic sequence. Bioinformatics 2022; 38:3164-3172. [PMID: 35389435 PMCID: PMC9890318 DOI: 10.1093/bioinformatics/btac214] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2021] [Revised: 03/04/2022] [Accepted: 04/06/2022] [Indexed: 02/04/2023] Open
Abstract
MOTIVATION Though genome-wide association studies have identified tens of thousands of variants associated with complex traits and most of them fall within the non-coding regions, they may not be the causal ones. The development of high-throughput functional assays leads to the discovery of experimental validated non-coding functional variants. However, these validated variants are rare due to technical difficulty and financial cost. The small sample size of validated variants makes it less reliable to develop a supervised machine learning model for achieving a whole genome-wide prediction of non-coding causal variants. RESULTS We will exploit a deep transfer learning model, which is based on convolutional neural network, to improve the prediction for functional non-coding variants (NCVs). To address the challenge of small sample size, the transfer learning model leverages both large-scale generic functional NCVs to improve the learning of low-level features and context-specific functional NCVs to learn high-level features toward the context-specific prediction task. By evaluating the deep transfer learning model on three MPRA datasets and 16 GWAS datasets, we demonstrate that the proposed model outperforms deep learning models without pretraining or retraining. In addition, the deep transfer learning model outperforms 18 existing computational methods in both MPRA and GWAS datasets. AVAILABILITY AND IMPLEMENTATION https://github.com/lichen-lab/TLVar. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Li Chen
- To whom correspondence should be addressed.
| | | | - Fengdi Zhao
- Department of Biostatistics and Health Data Science, Indiana University School of Medicine, Indianapolis, IN 46202, USA,Center for Computational Biology and Bioinformatics, Indiana University School of Medicine, Indianapolis, IN 46202, USA
| |
Collapse
|
97
|
Spielmann M, Kircher M. Computational and experimental methods for classifying variants of unknown clinical significance. Cold Spring Harb Mol Case Stud 2022; 8:mcs.a006196. [PMID: 35483875 PMCID: PMC9059783 DOI: 10.1101/mcs.a006196] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022] Open
Abstract
The increase in sequencing capacity, reduction in costs, and national and international coordinated efforts have led to the widespread introduction of next-generation sequencing (NGS) technologies in patient care. More generally, human genetics and genomic medicine are gaining importance for more and more patients. Some communities are already discussing the prospect of sequencing each individual's genome at time of birth. Together with digital health records, this shall enable individualized treatments and preventive measures, so-called precision medicine. A central step in this process is the identification of disease causal mutations or variant combinations that make us more susceptible for diseases. Although various technological advances have improved the identification of genetic alterations, the interpretation and ranking of the identified variants remains a major challenge. Based on our knowledge of molecular processes or previously identified disease variants, we can identify potentially functional genetic variants and, using different lines of evidence, we are sometimes able to demonstrate their pathogenicity directly. However, the vast majority of variants are classified as variants of uncertain clinical significance (VUSs) with not enough experimental evidence to determine their pathogenicity. In these cases, computational methods may be used to improve the prioritization and an increasing toolbox of experimental methods is emerging that can be used to assay the molecular effects of VUSs. Here, we discuss how computational and experimental methods can be used to create catalogs of variant effects for a variety of molecular and cellular phenotypes. We discuss the prospects of integrating large-scale functional data with machine learning and clinical knowledge for the development of accurate pathogenicity predictions for clinical applications.
Collapse
Affiliation(s)
- Malte Spielmann
- Institute of Human Genetics, University of Lübeck, 23562 Lübeck, Germany;,Institute of Human Genetics, Christian-Albrechts-Universität, 24105 Kiel, Germany;,Human Molecular Genomics Group, Max Planck Institute for Molecular Genetics, 14195 Berlin, Germany;,DZHK (German Centre for Cardiovascular Research), partner site Hamburg/Lübeck/Kiel, 23562 Lübeck, Germany
| | - Martin Kircher
- Institute of Human Genetics, University of Lübeck, 23562 Lübeck, Germany;,Berlin Institute of Health at Charité—Universitätsmedizin Berlin, 10117 Berlin, Germany;,DZHK (German Centre for Cardiovascular Research), partner site Berlin, 10115 Berlin, Germany
| |
Collapse
|
98
|
Abell NS, DeGorter MK, Gloudemans MJ, Greenwald E, Smith KS, He Z, Montgomery SB. Multiple causal variants underlie genetic associations in humans. Science 2022; 375:1247-1254. [PMID: 35298243 DOI: 10.1126/science.abj5117] [Citation(s) in RCA: 65] [Impact Index Per Article: 32.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Associations between genetic variation and traits are often in noncoding regions with strong linkage disequilibrium (LD), where a single causal variant is assumed to underlie the association. We applied a massively parallel reporter assay (MPRA) to functionally evaluate genetic variants in high, local LD for independent cis-expression quantitative trait loci (eQTL). We found that 17.7% of eQTLs exhibit more than one major allelic effect in tight LD. The detected regulatory variants were highly and specifically enriched for activating chromatin structures and allelic transcription factor binding. Integration of MPRA profiles with eQTL/complex trait colocalizations across 114 human traits and diseases identified causal variant sets demonstrating how genetic association signals can manifest through multiple, tightly linked causal variants.
Collapse
Affiliation(s)
- Nathan S Abell
- Department of Genetics, School of Medicine, Stanford University, Stanford, CA 94305, USA
| | - Marianne K DeGorter
- Department of Pathology, School of Medicine, Stanford University, Stanford, CA 94305, USA
| | | | - Emily Greenwald
- Department of Genetics, School of Medicine, Stanford University, Stanford, CA 94305, USA
| | - Kevin S Smith
- Department of Pathology, School of Medicine, Stanford University, Stanford, CA 94305, USA
| | - Zihuai He
- Department of Neurology and Neurological Sciences, Stanford University, Stanford, CA 94305, USA.,Quantitative Sciences Unit, Department of Medicine, Stanford University, Stanford, CA 94305, USA
| | - Stephen B Montgomery
- Department of Genetics, School of Medicine, Stanford University, Stanford, CA 94305, USA.,Department of Pathology, School of Medicine, Stanford University, Stanford, CA 94305, USA
| |
Collapse
|
99
|
Li X, Yung G, Zhou H, Sun R, Li Z, Hou K, Zhang MJ, Liu Y, Arapoglou T, Wang C, Ionita-Laza I, Lin X. A multi-dimensional integrative scoring framework for predicting functional variants in the human genome. Am J Hum Genet 2022; 109:446-456. [PMID: 35216679 PMCID: PMC8948160 DOI: 10.1016/j.ajhg.2022.01.017] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2021] [Accepted: 01/26/2022] [Indexed: 12/26/2022] Open
Abstract
Attempts to identify and prioritize functional DNA elements in coding and non-coding regions, particularly through use of in silico functional annotation data, continue to increase in popularity. However, specific functional roles can vary widely from one variant to another, making it challenging to summarize different aspects of variant function with a one-dimensional rating. Here we propose multi-dimensional annotation-class integrative estimation (MACIE), an unsupervised multivariate mixed-model framework capable of integrating annotations of diverse origin to assess multi-dimensional functional roles for both coding and non-coding variants. Unlike existing one-dimensional scoring methods, MACIE views variant functionality as a composite attribute encompassing multiple characteristics and estimates the joint posterior functional probabilities of each genomic position. This estimate offers more comprehensive and interpretable information in the presence of multiple aspects of functionality. Applied to a variety of independent coding and non-coding datasets, MACIE demonstrates powerful and robust performance in discriminating between functional and non-functional variants. We also show an application of MACIE to fine-mapping and heritability enrichment analysis by using the lipids GWAS summary statistics data from the European Network for Genetic and Genomic Epidemiology Consortium.
Collapse
Affiliation(s)
- Xihao Li
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA
| | - Godwin Yung
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA; Methods, Collaboration and Outreach Group, Genentech/Roche, South San Francisco, CA 94080, USA
| | - Hufeng Zhou
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA
| | - Ryan Sun
- Department of Biostatistics, University of Texas M.D. Anderson Cancer Center, Houston, TX 77030, USA
| | - Zilin Li
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA
| | - Kangcheng Hou
- Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, CA 90095, USA
| | - Martin Jinye Zhang
- Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA; Program in Medical and Population Genetics, Broad Institute of Harvard and MIT, Cambridge, MA 02142, USA
| | - Yaowu Liu
- School of Statistics, Southwestern University of Finance and Economics, Chengdu, Sichuan, China
| | - Theodore Arapoglou
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA
| | - Chen Wang
- Department of Biostatistics, Columbia University Mailman School of Public Health, New York, NY 10032, USA
| | - Iuliana Ionita-Laza
- Department of Biostatistics, Columbia University Mailman School of Public Health, New York, NY 10032, USA.
| | - Xihong Lin
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA; Program in Medical and Population Genetics, Broad Institute of Harvard and MIT, Cambridge, MA 02142, USA; Department of Statistics, Harvard University, Cambridge, MA, 02138, USA.
| |
Collapse
|
100
|
Jimeno-Martín A, Sousa E, Brocal-Ruiz R, Daroqui N, Maicas M, Flames N. Joint actions of diverse transcription factor families establish neuron-type identities and promote enhancer selectivity. Genome Res 2022; 32:459-473. [PMID: 35074859 PMCID: PMC8896470 DOI: 10.1101/gr.275623.121] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2021] [Accepted: 01/19/2022] [Indexed: 11/24/2022]
Abstract
To systematically investigate the complexity of neuron specification regulatory networks, we performed an RNA interference (RNAi) screen against all 875 transcription factors (TFs) encoded in Caenorhabditis elegans genome and searched for defects in nine different neuron types of the monoaminergic (MA) superclass and two cholinergic motoneurons. We identified 91 TF candidates to be required for correct generation of these neuron types, of which 28 were confirmed by mutant analysis. We found that correct reporter expression in each individual neuron type requires at least nine different TFs. Individual neuron types do not usually share TFs involved in their specification but share a common pattern of TFs belonging to the five most common TF families: homeodomain (HD), basic helix loop helix (bHLH), zinc finger (ZF), basic leucine zipper domain (bZIP), and nuclear hormone receptors (NHR). HD TF members are overrepresented, supporting a key role for this family in the establishment of neuronal identities. These five TF families are also prevalent when considering mutant alleles with previously reported neuronal phenotypes in C. elegans, Drosophila, and mouse. In addition, we studied terminal differentiation complexity focusing on the dopaminergic terminal regulatory program. We found two HD TFs (UNC-62 and VAB-3) that work together with known dopaminergic terminal selectors (AST-1, CEH-43, CEH-20). Combined TF binding sites for these five TFs constitute a cis-regulatory signature enriched in the regulatory regions of dopaminergic effector genes. Our results provide new insights on neuron-type regulatory programs in C. elegans that could help better understand neuron specification and evolution of neuron types.
Collapse
Affiliation(s)
- Angela Jimeno-Martín
- Developmental Neurobiology Unit, Instituto de Biomedicina de Valencia IBV-CSIC, Valencia, 46010, Spain
| | - Erick Sousa
- Developmental Neurobiology Unit, Instituto de Biomedicina de Valencia IBV-CSIC, Valencia, 46010, Spain
| | - Rebeca Brocal-Ruiz
- Developmental Neurobiology Unit, Instituto de Biomedicina de Valencia IBV-CSIC, Valencia, 46010, Spain
| | - Noemi Daroqui
- Developmental Neurobiology Unit, Instituto de Biomedicina de Valencia IBV-CSIC, Valencia, 46010, Spain
| | - Miren Maicas
- Developmental Neurobiology Unit, Instituto de Biomedicina de Valencia IBV-CSIC, Valencia, 46010, Spain
| | - Nuria Flames
- Developmental Neurobiology Unit, Instituto de Biomedicina de Valencia IBV-CSIC, Valencia, 46010, Spain
| |
Collapse
|