1
|
Sokolov V, Kyrchanova O, Klimenko N, Fedotova A, Ibragimov A, Maksimenko O, Georgiev P. New Drosophila promoter-associated architectural protein Mzfp1 interacts with CP190 and is required for housekeeping gene expression and insulator activity. Nucleic Acids Res 2024; 52:6886-6905. [PMID: 38769058 PMCID: PMC11229372 DOI: 10.1093/nar/gkae393] [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: 09/08/2023] [Revised: 04/20/2024] [Accepted: 05/10/2024] [Indexed: 05/22/2024] Open
Abstract
In Drosophila, a group of zinc finger architectural proteins recruits the CP190 protein to the chromatin, an interaction that is essential for the functional activity of promoters and insulators. In this study, we describe a new architectural C2H2 protein called Madf and Zinc-Finger Protein 1 (Mzfp1) that interacts with CP190. Mzfp1 has an unusual structure that includes six C2H2 domains organized in a C-terminal cluster and two tandem MADF domains. Mzfp1 predominantly binds to housekeeping gene promoters located in both euchromatin and heterochromatin genome regions. In vivo mutagenesis studies showed that Mzfp1 is an essential protein, and both MADF domains and the CP190 interaction region are required for its functional activity. The C2H2 cluster is sufficient for the specific binding of Mzfp1 to regulatory elements, while the second MADF domain is required for Mzfp1 recruitment to heterochromatin. Mzfp1 binds to the proximal part of the Fub boundary that separates regulatory domains of the Ubx and abd-A genes in the Bithorax complex. Mzfp1 participates in Fub functions in cooperation with the architectural proteins Pita and Su(Hw). Thus, Mzfp1 is a new architectural C2H2 protein involved in the organization of active promoters and insulators in Drosophila.
Collapse
Affiliation(s)
- Vladimir Sokolov
- Department of the Control of Genetic Processes, Institute of Gene Biology, Russian Academy of Sciences, Moscow 119334, Russia
| | - Olga Kyrchanova
- Department of the Control of Genetic Processes, Institute of Gene Biology, Russian Academy of Sciences, Moscow 119334, Russia
- Center for Precision Genome Editing and Genetic Technologies for Biomedicine, Institute of Gene Biology, Russian Academy of Sciences, Moscow 119334, Russia
| | - Natalia Klimenko
- Center for Precision Genome Editing and Genetic Technologies for Biomedicine, Institute of Gene Biology, Russian Academy of Sciences, Moscow 119334, Russia
| | - Anna Fedotova
- Department of the Control of Genetic Processes, Institute of Gene Biology, Russian Academy of Sciences, Moscow 119334, Russia
| | - Airat Ibragimov
- Center for Precision Genome Editing and Genetic Technologies for Biomedicine, Institute of Gene Biology, Russian Academy of Sciences, Moscow 119334, Russia
| | - Oksana Maksimenko
- Center for Precision Genome Editing and Genetic Technologies for Biomedicine, Institute of Gene Biology, Russian Academy of Sciences, Moscow 119334, Russia
| | - Pavel Georgiev
- Department of the Control of Genetic Processes, Institute of Gene Biology, Russian Academy of Sciences, Moscow 119334, Russia
| |
Collapse
|
2
|
Lodewijk GA, de Geus M, Guimarães RLFP, Jacobs FMJ. Emergence of the ZNF675 Gene During Primate Evolution-Influenced Human Neurodevelopment Through Changing HES1 Autoregulation. J Comp Neurol 2024; 532:e25648. [PMID: 38958676 DOI: 10.1002/cne.25648] [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: 12/31/2023] [Revised: 04/27/2024] [Accepted: 05/30/2024] [Indexed: 07/04/2024]
Abstract
In this study, we investigated recurrent copy number variations (CNVs) in the 19p12 locus, which are associated with neurodevelopmental disorders. The two genes in this locus, ZNF675 and ZNF681, arose via gene duplication in primates, and their presence in several pathological CNVs in the human population suggests that either or both of these genes are required for normal human brain development. ZNF675 and ZNF681 are members of the Krüppel-associated box zinc finger (KZNF) protein family, a class of transcriptional repressors important for epigenetic silencing of specific genomic regions. About 170 primate-specific KZNFs are present in the human genome. Although KZNFs are primarily associated with repressing retrotransposon-derived DNA, evidence is emerging that they can be co-opted for other gene regulatory processes. We show that genetic deletion of ZNF675 causes developmental defects in cortical organoids, and our data suggest that part of the observed neurodevelopmental phenotype is mediated by a gene regulatory role of ZNF675 on the promoter of the neurodevelopmental gene Hes family BHLH transcription factor 1 (HES1). We also find evidence for the recently evolved regulation of genes involved in neurological disorders, microcephalin 1 and sestrin 3. We show that ZNF675 interferes with HES1 auto-inhibition, a process essential for the maintenance of neural progenitors. As a striking example of how some KZNFs have integrated into preexisting gene expression networks, these findings suggest the emergence of ZNF675 has caused a change in the balance of HES1 autoregulation. The association of ZNF675 CNV with human developmental disorders and ZNF675-mediated regulation of neurodevelopmental genes suggests that it evolved into an important factor for human brain development.
Collapse
Affiliation(s)
- Gerrald A Lodewijk
- Swammerdam Institute for Life Sciences, Evolutionary Neurogenomics, University of Amsterdam, Amsterdam, The Netherlands
| | - Matthijs de Geus
- Swammerdam Institute for Life Sciences, Evolutionary Neurogenomics, University of Amsterdam, Amsterdam, The Netherlands
| | - Rita L F P Guimarães
- Swammerdam Institute for Life Sciences, Evolutionary Neurogenomics, University of Amsterdam, Amsterdam, The Netherlands
| | - Frank M J Jacobs
- Swammerdam Institute for Life Sciences, Evolutionary Neurogenomics, University of Amsterdam, Amsterdam, The Netherlands
| |
Collapse
|
3
|
Hong L, Kortemme T. An integrative approach to protein sequence design through multiobjective optimization. PLoS Comput Biol 2024; 20:e1011953. [PMID: 38991035 PMCID: PMC11265717 DOI: 10.1371/journal.pcbi.1011953] [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: 02/28/2024] [Revised: 07/23/2024] [Accepted: 06/25/2024] [Indexed: 07/13/2024] Open
Abstract
With recent methodological advances in the field of computational protein design, in particular those based on deep learning, there is an increasing need for frameworks that allow for coherent, direct integration of different models and objective functions into the generative design process. Here we demonstrate how evolutionary multiobjective optimization techniques can be adapted to provide such an approach. With the established Non-dominated Sorting Genetic Algorithm II (NSGA-II) as the optimization framework, we use AlphaFold2 and ProteinMPNN confidence metrics to define the objective space, and a mutation operator composed of ESM-1v and ProteinMPNN to rank and then redesign the least favorable positions. Using the two-state design problem of the foldswitching protein RfaH as an in-depth case study, and PapD and calmodulin as examples of higher-dimensional design problems, we show that the evolutionary multiobjective optimization approach leads to significant reduction in the bias and variance in RfaH native sequence recovery, compared to a direct application of ProteinMPNN. We suggest that this improvement is due to three factors: (i) the use of an informative mutation operator that accelerates the sequence space exploration, (ii) the parallel, iterative design process inherent to the genetic algorithm that improves upon the ProteinMPNN autoregressive sequence decoding scheme, and (iii) the explicit approximation of the Pareto front that leads to optimal design candidates representing diverse tradeoff conditions. We anticipate this approach to be readily adaptable to different models and broadly relevant for protein design tasks with complex specifications.
Collapse
Affiliation(s)
- Lu Hong
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, California, United States of America
| | - Tanja Kortemme
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, California, United States of America
- Quantitative Biosciences Institute, University of California, San Francisco, California, United States of America
- Chan Zuckerberg Biohub, San Francisco, California, United States of America
| |
Collapse
|
4
|
Hong L, Kortemme T. An integrative approach to protein sequence design through multiobjective optimization. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.03.01.582670. [PMID: 38496480 PMCID: PMC10942313 DOI: 10.1101/2024.03.01.582670] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/19/2024]
Abstract
With recent methodological advances in the field of computational protein design, in particular those based on deep learning, there is an increasing need for frameworks that allow for coherent, direct integration of different models and objective functions into the generative design process. Here we demonstrate how evolutionary multiobjective optimization techniques can be adapted to provide such an approach. With the established Non-dominated Sorting Genetic Algorithm II (NSGA-II) as the optimization framework, we use AlphaFold2 and ProteinMPNN confidence metrics to define the objective space, and a mutation operator composed of ESM-1v and ProteinMPNN to rank and then redesign the least favorable positions. Using the multistate design problem of the foldswitching protein RfaH as an in-depth case study, we show that the evolutionary multiobjective optimization approach leads to significant reduction in the bias and variance in RfaH native sequence recovery, compared to a direct application of ProteinMPNN. We suggest that this improvement is due to three factors: (i) the use of an informative mutation operator that accelerates the sequence space exploration, (ii) the parallel, iterative design process inherent to the genetic algorithm that improves upon the ProteinMPNN autoregressive sequence decoding scheme, and (iii) the explicit approximation of the Pareto front that leads to optimal design candidates representing diverse tradeoff conditions. We anticipate this approach to be readily adaptable to different models and broadly relevant for protein design tasks with complex specifications.
Collapse
Affiliation(s)
- Lu Hong
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Tanja Kortemme
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, San Francisco, CA 94158, USA
- Quantitative Biosciences Institute, University of California, San Francisco, San Francisco, CA 94158, USA
- Chan Zuckerberg Biohub, San Francisco, CA 94158, USA
| |
Collapse
|
5
|
Sukmak R, Suttinun C, Kovitvadhi U, Kovitvadhi A, Vongsangnak W. Uncovering nutrients and energy related gene functions of black soldier fly Hermetia illucens strain KUP. Gene 2024; 896:148045. [PMID: 38042219 DOI: 10.1016/j.gene.2023.148045] [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] [Received: 08/30/2023] [Revised: 11/18/2023] [Accepted: 11/28/2023] [Indexed: 12/04/2023]
Abstract
The black soldier fly (Hermetia illucens) has emerged as a significant insect species in the decomposition of organic waste for sustainable agricultural practices. Due to its remarkable characteristics and performance, H. illucens is increasingly utilised for insect farming, particularly for industrial-scale rearing throughout the world. In this study, we employed whole-genome sequencing to annotate the gene and protein functions of H. illucens and to explore the functional genomics related to nutrients and energy. As a result, a genome size of H. illucens strain KUP 1.68 Gb with a GC content of 42.13 % was achieved. Of the 14,036 coding sequences, we determined the function of 12,046 protein-coding genes. Based on metabolic functional assignment, we classified 4,218 protein-coding genes; the main category was metabolism (32.86 %). Comparative genomic analysis across the other H. illucens strain and insect species revealed that the major metabolic gene functions and pathways related to nutrient and energy sources of H. illucens KUP are involved in key amino acid metabolism (e.g., cysteine and methionine) as well as fatty acid biosynthesis and glycerolipid metabolism. These findings underscore the metabolic capability and versatility of H. illucens, which is regarded as a potential source of proteins and lipids. Our study contributes to the knowledge regarding the feed utilisation of H. illucens and offers insights into transforming waste into valuable products. H. illucens has the potential to create globally sustainable nutrients and environmentally friendly solutions, aligning with the goal of responsible resource utilisation.
Collapse
Affiliation(s)
- Rachrapee Sukmak
- Graduate Student in Animal Health and Biomedical Science Program, Faculty of Veterinary Medicine, Kasetsart University, Bangkok, Thailand
| | - Chanaporn Suttinun
- Department of Veterinary Biosciences and Veterinary Public Health, Faculty of Veterinary Medicine, Chiang Mai University, Chiang Mai, Thailand
| | - Uthaiwan Kovitvadhi
- Department of Zoology, Faculty of Science, Kasetsart University, Bangkok, Thailand
| | - Attawit Kovitvadhi
- Department of Physiology, Faculty of Veterinary Medicine, Kasetsart University, Bangkok, Thailand.
| | - Wanwipa Vongsangnak
- Department of Zoology, Faculty of Science, Kasetsart University, Bangkok, Thailand; Omics Center for Agriculture, Bioresources, Food, and Health, Kasetsart University (OmiKU), Bangkok, Thailand.
| |
Collapse
|
6
|
Wells JN, Chang NC, McCormick J, Coleman C, Ramos N, Jin B, Feschotte C. Transposable elements drive the evolution of metazoan zinc finger genes. Genome Res 2023; 33:1325-1339. [PMID: 37714714 PMCID: PMC10547256 DOI: 10.1101/gr.277966.123] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2023] [Accepted: 06/15/2023] [Indexed: 09/17/2023]
Abstract
Cys2-His2 zinc finger genes (ZNFs) form the largest family of transcription factors in metazoans. ZNF evolution is highly dynamic and characterized by the rapid expansion and contraction of numerous subfamilies across the animal phylogeny. The forces and mechanisms underlying rapid ZNF evolution remain poorly understood, but there is growing evidence that, in tetrapods, the targeting and repression of lineage-specific transposable elements (TEs) plays a critical role in the evolution of the Krüppel-associated box ZNF (KZNF) subfamily. Currently, it is unknown whether this function and coevolutionary relationship is unique to KZNFs or is a broader feature of metazoan ZNFs. Here, we present evidence that genomic conflict with TEs has been a central driver of the diversification of ZNFs in animals. Sampling from 3221 genome assemblies, we show that the copy number of retroelements correlates with that of ZNFs across at least 750 million years of metazoan evolution. Using computational predictions, we show that ZNFs preferentially bind TEs in diverse animal species. We further investigate the largest ZNF subfamily found in cyprinid fish, which is characterized by a conserved sequence we dubbed the fish N-terminal zinc finger-associated (FiNZ) domain. Zebrafish possess approximately 700 FiNZ-ZNFs, many of which are evolving adaptively under positive selection. Like mammalian KZNFs, most zebrafish FiNZ-ZNFs are expressed at the onset of zygotic genome activation, and blocking their translation using morpholinos during early embryogenesis results in derepression of transcriptionally active TEs. Together, these data suggest that ZNF diversification has been intimately connected to TE expansion throughout animal evolution.
Collapse
Affiliation(s)
- Jonathan N Wells
- Department of Molecular Biology and Genetics, Cornell University, Ithaca, New York 14850, USA;
| | - Ni-Chen Chang
- Department of Molecular Biology and Genetics, Cornell University, Ithaca, New York 14850, USA
| | - John McCormick
- Department of Molecular Biology and Genetics, Cornell University, Ithaca, New York 14850, USA
| | - Caitlyn Coleman
- Department of Cell Biology, Microbiology and Molecular Biology, University of South Florida, Tampa, Florida 33620, USA
| | - Nathalie Ramos
- Department of Molecular Biology and Genetics, Cornell University, Ithaca, New York 14850, USA
- Department of Genetics and Genomic Sciences, Center for Transformative Disease Modeling, Tisch Cancer Institute, Icahn Institute for Data Science and Genomic Technology, Icahn School of Medicine at Mount Sinai, New York, New York 10029, USA
| | - Bozhou Jin
- Department of Molecular Biology and Genetics, Cornell University, Ithaca, New York 14850, USA
| | - Cédric Feschotte
- Department of Molecular Biology and Genetics, Cornell University, Ithaca, New York 14850, USA;
| |
Collapse
|
7
|
de Tribolet-Hardy J, Thorball CW, Forey R, Planet E, Duc J, Coudray A, Khubieh B, Offner S, Pulver C, Fellay J, Imbeault M, Turelli P, Trono D. Genetic features and genomic targets of human KRAB-zinc finger proteins. Genome Res 2023; 33:1409-1423. [PMID: 37730438 PMCID: PMC10547255 DOI: 10.1101/gr.277722.123] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2023] [Accepted: 07/18/2023] [Indexed: 09/22/2023]
Abstract
Krüppel-associated box (KRAB) domain-containing zinc finger proteins (KZFPs) are one of the largest groups of transcription factors encoded by tetrapods, with 378 members in human alone. KZFP genes are often grouped in clusters reflecting amplification by gene and segment duplication since the gene family first emerged more than 400 million years ago. Previous work has revealed that many KZFPs recognize transposable element (TE)-embedded sequences as genomic targets, and that KZFPs facilitate the co-option of the regulatory potential of TEs for the benefit of the host. Here, we present a comprehensive survey of the genetic features and genomic targets of human KZFPs, notably completing past analyses by adding data on close to a hundred family members. General principles emerge from our study of the TE-KZFP regulatory system, which point to multipronged evolutionary mechanisms underlaid by highly complex and combinatorial modes of action with strong influences on human speciation.
Collapse
Affiliation(s)
- Jonas de Tribolet-Hardy
- School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Christian W Thorball
- School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Romain Forey
- School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Evarist Planet
- School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Julien Duc
- School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Alexandre Coudray
- School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Bara Khubieh
- School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Sandra Offner
- School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Cyril Pulver
- School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Jacques Fellay
- School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
- Precision Medicine Unit, Lausanne University Hospital (CHUV) and University of Lausanne, 1010 Lausanne, Switzerland
| | - Michael Imbeault
- School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Priscilla Turelli
- School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Didier Trono
- School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland;
| |
Collapse
|
8
|
Ray D, Laverty KU, Jolma A, Nie K, Samson R, Pour SE, Tam CL, von Krosigk N, Nabeel-Shah S, Albu M, Zheng H, Perron G, Lee H, Najafabadi H, Blencowe B, Greenblatt J, Morris Q, Hughes TR. RNA-binding proteins that lack canonical RNA-binding domains are rarely sequence-specific. Sci Rep 2023; 13:5238. [PMID: 37002329 PMCID: PMC10066285 DOI: 10.1038/s41598-023-32245-9] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2022] [Accepted: 03/23/2023] [Indexed: 04/03/2023] Open
Abstract
Thousands of RNA-binding proteins (RBPs) crosslink to cellular mRNA. Among these are numerous unconventional RBPs (ucRBPs)-proteins that associate with RNA but lack known RNA-binding domains (RBDs). The vast majority of ucRBPs have uncharacterized RNA-binding specificities. We analyzed 492 human ucRBPs for intrinsic RNA-binding in vitro and identified 23 that bind specific RNA sequences. Most (17/23), including 8 ribosomal proteins, were previously associated with RNA-related function. We identified the RBDs responsible for sequence-specific RNA-binding for several of these 23 ucRBPs and surveyed whether corresponding domains from homologous proteins also display RNA sequence specificity. CCHC-zf domains from seven human proteins recognized specific RNA motifs, indicating that this is a major class of RBD. For Nudix, HABP4, TPR, RanBP2-zf, and L7Ae domains, however, only isolated members or closely related homologs yielded motifs, consistent with RNA-binding as a derived function. The lack of sequence specificity for most ucRBPs is striking, and we suggest that many may function analogously to chromatin factors, which often crosslink efficiently to cellular DNA, presumably via indirect recruitment. Finally, we show that ucRBPs tend to be highly abundant proteins and suggest their identification in RNA interactome capture studies could also result from weak nonspecific interactions with RNA.
Collapse
Affiliation(s)
- Debashish Ray
- Donnelly Centre, University of Toronto, Toronto, ON, M5S 3E1, Canada
| | - Kaitlin U Laverty
- Donnelly Centre, University of Toronto, Toronto, ON, M5S 3E1, Canada
- Department of Molecular Genetics, University of Toronto, Toronto, ON, M5S 1A8, Canada
| | - Arttu Jolma
- Donnelly Centre, University of Toronto, Toronto, ON, M5S 3E1, Canada
| | - Kate Nie
- Donnelly Centre, University of Toronto, Toronto, ON, M5S 3E1, Canada
- Department of Molecular Genetics, University of Toronto, Toronto, ON, M5S 1A8, Canada
| | - Reuben Samson
- Donnelly Centre, University of Toronto, Toronto, ON, M5S 3E1, Canada
- Department of Molecular Genetics, University of Toronto, Toronto, ON, M5S 1A8, Canada
| | - Sara E Pour
- Donnelly Centre, University of Toronto, Toronto, ON, M5S 3E1, Canada
- Department of Molecular Genetics, University of Toronto, Toronto, ON, M5S 1A8, Canada
| | - Cyrus L Tam
- Computational and Systems Biology Program, Memorial Sloan Kettering Cancer Center, New York, NY, USA
- Tri-Institutional Training Program in Computational Biology and Medicine, Weill Cornell Medicine, New York, NY, USA
| | - Niklas von Krosigk
- Donnelly Centre, University of Toronto, Toronto, ON, M5S 3E1, Canada
- Department of Molecular Genetics, University of Toronto, Toronto, ON, M5S 1A8, Canada
| | - Syed Nabeel-Shah
- Donnelly Centre, University of Toronto, Toronto, ON, M5S 3E1, Canada
- Department of Molecular Genetics, University of Toronto, Toronto, ON, M5S 1A8, Canada
| | - Mihai Albu
- Donnelly Centre, University of Toronto, Toronto, ON, M5S 3E1, Canada
| | - Hong Zheng
- Donnelly Centre, University of Toronto, Toronto, ON, M5S 3E1, Canada
| | - Gabrielle Perron
- Department of Human Genetics, McGill University, Montréal, QC, H3A 0C7, Canada
- McGill Genome Centre, Montréal, QC, H3A 0G1, Canada
| | - Hyunmin Lee
- Donnelly Centre, University of Toronto, Toronto, ON, M5S 3E1, Canada
| | - Hamed Najafabadi
- Department of Human Genetics, McGill University, Montréal, QC, H3A 0C7, Canada
- McGill Genome Centre, Montréal, QC, H3A 0G1, Canada
| | - Benjamin Blencowe
- Donnelly Centre, University of Toronto, Toronto, ON, M5S 3E1, Canada
- Department of Molecular Genetics, University of Toronto, Toronto, ON, M5S 1A8, Canada
| | - Jack Greenblatt
- Donnelly Centre, University of Toronto, Toronto, ON, M5S 3E1, Canada
- Department of Molecular Genetics, University of Toronto, Toronto, ON, M5S 1A8, Canada
| | - Quaid Morris
- Donnelly Centre, University of Toronto, Toronto, ON, M5S 3E1, Canada.
- Department of Molecular Genetics, University of Toronto, Toronto, ON, M5S 1A8, Canada.
- Computational and Systems Biology Program, Memorial Sloan Kettering Cancer Center, New York, NY, USA.
- Tri-Institutional Training Program in Computational Biology and Medicine, Weill Cornell Medicine, New York, NY, USA.
| | - Timothy R Hughes
- Donnelly Centre, University of Toronto, Toronto, ON, M5S 3E1, Canada.
- Department of Molecular Genetics, University of Toronto, Toronto, ON, M5S 1A8, Canada.
| |
Collapse
|
9
|
Rozowsky J, Gao J, Borsari B, Yang YT, Galeev T, Gürsoy G, Epstein CB, Xiong K, Xu J, Li T, Liu J, Yu K, Berthel A, Chen Z, Navarro F, Sun MS, Wright J, Chang J, Cameron CJF, Shoresh N, Gaskell E, Drenkow J, Adrian J, Aganezov S, Aguet F, Balderrama-Gutierrez G, Banskota S, Corona GB, Chee S, Chhetri SB, Cortez Martins GC, Danyko C, Davis CA, Farid D, Farrell NP, Gabdank I, Gofin Y, Gorkin DU, Gu M, Hecht V, Hitz BC, Issner R, Jiang Y, Kirsche M, Kong X, Lam BR, Li S, Li B, Li X, Lin KZ, Luo R, Mackiewicz M, Meng R, Moore JE, Mudge J, Nelson N, Nusbaum C, Popov I, Pratt HE, Qiu Y, Ramakrishnan S, Raymond J, Salichos L, Scavelli A, Schreiber JM, Sedlazeck FJ, See LH, Sherman RM, Shi X, Shi M, Sloan CA, Strattan JS, Tan Z, Tanaka FY, Vlasova A, Wang J, Werner J, Williams B, Xu M, Yan C, Yu L, Zaleski C, Zhang J, Ardlie K, Cherry JM, Mendenhall EM, Noble WS, Weng Z, Levine ME, Dobin A, Wold B, Mortazavi A, Ren B, Gillis J, Myers RM, Snyder MP, Choudhary J, Milosavljevic A, Schatz MC, Bernstein BE, Guigó R, Gingeras TR, Gerstein M. The EN-TEx resource of multi-tissue personal epigenomes & variant-impact models. Cell 2023; 186:1493-1511.e40. [PMID: 37001506 PMCID: PMC10074325 DOI: 10.1016/j.cell.2023.02.018] [Citation(s) in RCA: 16] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2022] [Revised: 10/16/2022] [Accepted: 02/10/2023] [Indexed: 04/03/2023]
Abstract
Understanding how genetic variants impact molecular phenotypes is a key goal of functional genomics, currently hindered by reliance on a single haploid reference genome. Here, we present the EN-TEx resource of 1,635 open-access datasets from four donors (∼30 tissues × ∼15 assays). The datasets are mapped to matched, diploid genomes with long-read phasing and structural variants, instantiating a catalog of >1 million allele-specific loci. These loci exhibit coordinated activity along haplotypes and are less conserved than corresponding, non-allele-specific ones. Surprisingly, a deep-learning transformer model can predict the allele-specific activity based only on local nucleotide-sequence context, highlighting the importance of transcription-factor-binding motifs particularly sensitive to variants. Furthermore, combining EN-TEx with existing genome annotations reveals strong associations between allele-specific and GWAS loci. It also enables models for transferring known eQTLs to difficult-to-profile tissues (e.g., from skin to heart). Overall, EN-TEx provides rich data and generalizable models for more accurate personal functional genomics.
Collapse
Affiliation(s)
- Joel Rozowsky
- Section on Biomedical Informatics and Data Science, Yale University, New Haven, CT, USA; Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Jiahao Gao
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Beatrice Borsari
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA; Centre for Genomic Regulation, The Barcelona Institute of Science and Technology, Barcelona, Catalonia, Spain
| | - Yucheng T Yang
- Institute of Science and Technology for Brain-Inspired Intelligence; MOE Key Laboratory of Computational Neuroscience and Brain-Inspired Intelligence; MOE Frontiers Center for Brain Science, Fudan University, Shanghai 200433, China; Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Timur Galeev
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Gamze Gürsoy
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | | | - Kun Xiong
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Jinrui Xu
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Tianxiao Li
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Jason Liu
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Keyang Yu
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, USA
| | - Ana Berthel
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Zhanlin Chen
- Department of Statistics and Data Science, Yale University, New Haven, CT, USA
| | - Fabio Navarro
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Maxwell S Sun
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | | | - Justin Chang
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Christopher J F Cameron
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Noam Shoresh
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | | | - Jorg Drenkow
- Functional Genomics, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
| | - Jessika Adrian
- Department of Genetics, School of Medicine, Stanford University, Palo Alto, CA, USA
| | - Sergey Aganezov
- Departments of Computer Science and Biology, Johns Hopkins University, Baltimore, MD, USA
| | | | | | | | | | - Sora Chee
- Ludwig Institute for Cancer Research, University of California, San Diego, La Jolla, CA, USA
| | - Surya B Chhetri
- HudsonAlpha Institute for Biotechnology, Huntsville, AL, USA
| | - Gabriel Conte Cortez Martins
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Cassidy Danyko
- Functional Genomics, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
| | - Carrie A Davis
- Functional Genomics, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
| | - Daniel Farid
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | | | - Idan Gabdank
- Department of Genetics, School of Medicine, Stanford University, Palo Alto, CA, USA
| | - Yoel Gofin
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, USA
| | - David U Gorkin
- Ludwig Institute for Cancer Research, University of California, San Diego, La Jolla, CA, USA
| | - Mengting Gu
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Vivian Hecht
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Benjamin C Hitz
- Department of Genetics, School of Medicine, Stanford University, Palo Alto, CA, USA
| | - Robbyn Issner
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Yunzhe Jiang
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Melanie Kirsche
- Departments of Computer Science and Biology, Johns Hopkins University, Baltimore, MD, USA
| | - Xiangmeng Kong
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Bonita R Lam
- Department of Genetics, School of Medicine, Stanford University, Palo Alto, CA, USA
| | - Shantao Li
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Bian Li
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Xiqi Li
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, USA
| | - Khine Zin Lin
- Department of Genetics, School of Medicine, Stanford University, Palo Alto, CA, USA
| | - Ruibang Luo
- Department of Computer Science, The University of Hong Kong, Hong Kong, CHN
| | - Mark Mackiewicz
- HudsonAlpha Institute for Biotechnology, Huntsville, AL, USA
| | - Ran Meng
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Jill E Moore
- Program in Bioinformatics and Integrative Biology, University of Massachusetts Medical School, Worcester, MA, USA
| | - Jonathan Mudge
- European Bioinformatics Institute, Cambridge, Cambridgeshire, GB
| | | | - Chad Nusbaum
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Ioann Popov
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Henry E Pratt
- Program in Bioinformatics and Integrative Biology, University of Massachusetts Medical School, Worcester, MA, USA
| | - Yunjiang Qiu
- Ludwig Institute for Cancer Research, University of California, San Diego, La Jolla, CA, USA
| | - Srividya Ramakrishnan
- Departments of Computer Science and Biology, Johns Hopkins University, Baltimore, MD, USA
| | - Joe Raymond
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Leonidas Salichos
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA; Department of Biological and Chemical Sciences, New York Institute of Technology, Old Westbury, NY, USA
| | - Alexandra Scavelli
- Functional Genomics, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
| | - Jacob M Schreiber
- Department of Genome Sciences, University of Washington, Seattle, WA, USA
| | - Fritz J Sedlazeck
- Departments of Computer Science and Biology, Johns Hopkins University, Baltimore, MD, USA; Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA; Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX, USA
| | - Lei Hoon See
- Functional Genomics, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
| | - Rachel M Sherman
- Departments of Computer Science and Biology, Johns Hopkins University, Baltimore, MD, USA
| | - Xu Shi
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Minyi Shi
- Department of Genetics, School of Medicine, Stanford University, Palo Alto, CA, USA
| | - Cricket Alicia Sloan
- Department of Genetics, School of Medicine, Stanford University, Palo Alto, CA, USA
| | - J Seth Strattan
- Department of Genetics, School of Medicine, Stanford University, Palo Alto, CA, USA
| | - Zhen Tan
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Forrest Y Tanaka
- Department of Genetics, School of Medicine, Stanford University, Palo Alto, CA, USA
| | - Anna Vlasova
- Centre for Genomic Regulation, The Barcelona Institute of Science and Technology, Barcelona, Catalonia, Spain; Comparative Genomics Group, Life Science Programme, Barcelona Supercomputing Centre, Barcelona, Spain; Institute of Research in Biomedicine, Barcelona, Spain
| | - Jun Wang
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Jonathan Werner
- Functional Genomics, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
| | - Brian Williams
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, CA, USA
| | - Min Xu
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Chengfei Yan
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA
| | - Lu Yu
- Institute of Cancer Research, London, UK
| | - Christopher Zaleski
- Functional Genomics, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
| | - Jing Zhang
- Department of Computer Science, University of California, Irvine, Irvine, CA, USA
| | | | - J Michael Cherry
- Department of Genetics, School of Medicine, Stanford University, Palo Alto, CA, USA
| | | | - William S Noble
- Department of Genome Sciences, University of Washington, Seattle, WA, USA
| | - Zhiping Weng
- Program in Bioinformatics and Integrative Biology, University of Massachusetts Medical School, Worcester, MA, USA
| | - Morgan E Levine
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Pathology, Yale University School of Medicine, New Haven, CT, USA
| | - Alexander Dobin
- Functional Genomics, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
| | - Barbara Wold
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, CA, USA
| | - Ali Mortazavi
- Department of Developmental and Cell Biology, University of California, Irvine, Irvine, CA, USA
| | - Bing Ren
- Ludwig Institute for Cancer Research, University of California, San Diego, La Jolla, CA, USA
| | - Jesse Gillis
- Functional Genomics, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA; Department of Physiology, University of Toronto, Toronto, ON, Canada
| | - Richard M Myers
- HudsonAlpha Institute for Biotechnology, Huntsville, AL, USA
| | - Michael P Snyder
- Department of Genetics, School of Medicine, Stanford University, Palo Alto, CA, USA
| | | | | | - Michael C Schatz
- Departments of Computer Science and Biology, Johns Hopkins University, Baltimore, MD, USA; Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA.
| | - Bradley E Bernstein
- Broad Institute of MIT and Harvard, Cambridge, MA, USA; Department of Cancer Biology, Dana-Farber Cancer Institute, Boston, MA, USA.
| | - Roderic Guigó
- Centre for Genomic Regulation, The Barcelona Institute of Science and Technology, Barcelona, Catalonia, Spain; Universitat Pompeu Fabra, Barcelona, Catalonia, Spain.
| | - Thomas R Gingeras
- Functional Genomics, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA.
| | - Mark Gerstein
- Section on Biomedical Informatics and Data Science, Yale University, New Haven, CT, USA; Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, USA; Department of Statistics and Data Science, Yale University, New Haven, CT, USA; Department of Computer Science, Yale University, New Haven, CT, USA.
| |
Collapse
|
10
|
Wagner A. Adaptive evolvability through direct selection instead of indirect, second-order selection. JOURNAL OF EXPERIMENTAL ZOOLOGY. PART B, MOLECULAR AND DEVELOPMENTAL EVOLUTION 2022; 338:395-404. [PMID: 34254439 PMCID: PMC9786751 DOI: 10.1002/jez.b.23071] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/30/2020] [Revised: 05/11/2021] [Accepted: 06/04/2021] [Indexed: 12/30/2022]
Abstract
Can evolvability itself be the product of adaptive evolution? To answer this question is challenging, because any DNA mutation that alters only evolvability is subject to indirect, "second order" selection on the future effects of this mutation. Such indirect selection is weaker than "first-order" selection on mutations that alter fitness, in the sense that it can operate only under restrictive conditions. Here I discuss a route to adaptive evolvability that overcomes this challenge. Specifically, a recent evolution experiment showed that some mutations can enhance both fitness and evolvability through a combination of direct and indirect selection. Unrelated evidence from gene duplication and the evolution of gene regulation suggests that mutations with such dual effects may not be rare. Through such mutations, evolvability may increase at least in part because it provides an adaptive advantage. These observations suggest a research program on the adaptive evolution of evolvability, which aims to identify such mutations and to disentangle their direct fitness effects from their indirect effects on evolvability. If evolvability is itself adaptive, Darwinian evolution may have created more than life's diversity. It may also have helped create the very conditions that made the success of Darwinian evolution possible.
Collapse
Affiliation(s)
- Andreas Wagner
- Department of Evolutionary Biology and Environmental StudiesUniversity of ZurichZurichSwitzerland,Swiss Institute of BioinformaticsQuartier Sorge‐Batiment GenopodeLausanneSwitzerland,The Santa Fe InstituteSanta FeNew MexicoUSA,Stellenbosch Institute for Advanced Study, Wallenberg Research Centre at Stellenbosch UniversityStellenboschSouth Africa
| |
Collapse
|
11
|
Li P, Yu A, Sun R, Liu A. Function and Evolution of C1-2i Subclass of C2H2-Type Zinc Finger Transcription Factors in POPLAR. Genes (Basel) 2022; 13:genes13101843. [PMID: 36292728 PMCID: PMC9602059 DOI: 10.3390/genes13101843] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2022] [Revised: 09/24/2022] [Accepted: 10/11/2022] [Indexed: 11/16/2022] Open
Abstract
C2H2 zinc finger (C2H2-ZF) transcription factors participate in various aspects of normal plant growth regulation and stress responses. C1-2i C2H2-ZFs are a special subclass of conserved proteins that contain two ZnF-C2H2 domains. Some C1-2i C2H2-ZFs in Arabidopsis (ZAT) are involved in stress resistance and other functions. However, there is limited information on C1-2i C2H2-ZFs in Populus trichocarpa (PtriZATs). To analyze the function and evolution of C1-2i C2H2-ZFs, eleven PtriZATs were identified in P. trichocarpa, which can be classified into two subgroups. The protein structure, conserved ZnF-C2H2 domains and QALGGH motifs, showed high conservation during the evolution of PtriZATs in P. trichocarpa. The spacing between two ZnF-C2H2 domains, chromosomal locations and cis-elements implied the original proteins and function of PtriZATs. Furthermore, the gene expression of different tissues and stress treatment showed the functional differentiation of PtriZATs subgroups and their stress response function. The analysis of C1-2i C2H2-ZFs in different Populus species and plants implied their evolution and differentiation, especially in terms of stress resistance. Cis-elements and expression pattern analysis of interaction proteins implied the function of PtriZATs through binding with stress-related genes, which are involved in gene regulation by via epigenetic modification through histone regulation, DNA methylation, ubiquitination, etc. Our results for the origin and evolution of PtriZATs will contribute to understanding the functional differentiation of C1-2i C2H2-ZFs in P. trichocarpa. The interaction and expression results will lay a foundation for the further functional investigation of their roles and biological processes in Populus.
Collapse
|
12
|
Aizenshtein-Gazit S, Orenstein Y. DeepZF: improved DNA-binding prediction of C2H2-zinc-finger proteins by deep transfer learning. Bioinformatics 2022; 38:ii62-ii67. [PMID: 36124796 DOI: 10.1093/bioinformatics/btac469] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
MOTIVATION Cys2His2 zinc-finger (C2H2-ZF) proteins are the largest class of human transcription factors and hence play central roles in gene regulation and cell function. C2H2-ZF proteins are characterized by a DNA-binding domain containing multiple ZFs. A subset of the ZFs bind diverse DNA triplets. Despite their central roles, little is known about which of their ZFs are binding and how the DNA-binding preferences are encoded in the amino acid sequence of each ZF. RESULTS We present DeepZF, a deep-learning-based pipeline for predicting binding ZFs and their DNA-binding preferences given only the amino acid sequence of a C2H2-ZF protein. To the best of our knowledge, we compiled the first in vivo dataset of binding and non-binding ZFs for training the first ZF-binding classifier. Our classifier, which is based on a novel protein transformer, achieved an average AUROC of 0.71. Moreover, we took advantage of both in vivo and in vitro datasets to learn the recognition code of ZF-DNA binding through transfer learning. Our newly developed model, which is the first to utilize deep learning for the task, achieved an average Pearson correlation greater than 0.94 over each of the three DNA binding positions. Together, DeepZF outperformed extant methods in the task of C2H2-ZF protein DNA-binding preferences prediction: it achieved an average Pearson correlation of 0.42 in motif similarity compared with an average correlation smaller than 0.1 achieved by extant methods. By applying established interpretability techniques, we show that DeepZF inferred biologically relevant binding principles, such as the effect of amino acid residue positions on ZF DNA-binding potential. AVAILABILITY AND IMPLEMENTATION DeepZF code, model, and results are available via github.com/OrensteinLab/DeepZF. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Sofia Aizenshtein-Gazit
- School of Electrical and Computer Engineering, Ben-Gurion University of the Negev, Beer-Sheva 8410501, Israel
| | - Yaron Orenstein
- School of Electrical and Computer Engineering, Ben-Gurion University of the Negev, Beer-Sheva 8410501, Israel
| |
Collapse
|
13
|
Neiro J, Sridhar D, Dattani A, Aboobaker A. Identification of putative enhancer-like elements predicts regulatory networks active in planarian adult stem cells. eLife 2022; 11:79675. [PMID: 35997250 PMCID: PMC9522251 DOI: 10.7554/elife.79675] [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: 05/05/2022] [Accepted: 08/23/2022] [Indexed: 11/13/2022] Open
Abstract
Planarians have become an established model system to study regeneration and stem cells, but the regulatory elements in the genome remain almost entirely undescribed. Here, by integrating epigenetic and expression data we use multiple sources of evidence to predict enhancer elements active in the adult stem cell populations that drive regeneration. We have used ChIP-seq data to identify genomic regions with histone modifications consistent with enhancer activity, and ATAC-seq data to identify accessible chromatin. Overlapping these signals allowed for the identification of a set of high-confidence candidate enhancers predicted to be active in planarian adult stem cells. These enhancers are enriched for predicted transcription factor (TF) binding sites for TFs and TF families expressed in planarian adult stem cells. Footprinting analyses provided further evidence that these potential TF binding sites are likely to be occupied in adult stem cells. We integrated these analyses to build testable hypotheses for the regulatory function of TFs in stem cells, both with respect to how pluripotency might be regulated, and to how lineage differentiation programs are controlled. We found that our predicted GRNs were independently supported by existing TF RNAi/RNA-seq datasets, providing further evidence that our work predicts active enhancers that regulate adult stem cells and regenerative mechanisms.
Collapse
Affiliation(s)
- Jakke Neiro
- Department of Zoology, University of Oxford, Oxford, United Kingdom
| | - Divya Sridhar
- Department of Zoology, University of Oxford, Oxford, United Kingdom
| | - Anish Dattani
- Living Systems Institute, University of Exeter, Exeter, United Kingdom
| | - Aziz Aboobaker
- Department of Zoology, University of Oxford, Oxford, United Kingdom
| |
Collapse
|
14
|
Ray S, Tillo D, Assad N, Ufot A, Porollo A, Durell SR, Vinson C. Altering the Double-Stranded DNA Specificity of the bZIP Domain of Zta with Site-Directed Mutagenesis at N182. ACS OMEGA 2022; 7:129-139. [PMID: 35036684 PMCID: PMC8756438 DOI: 10.1021/acsomega.1c04148] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/03/2021] [Accepted: 11/23/2021] [Indexed: 06/14/2023]
Abstract
Zta, the Epstein-Barr virus bZIP transcription factor (TF), binds both unmethylated and methylated double-stranded DNA (dsDNA) in a sequence-specific manner. We studied the contribution of a conserved asparagine (N182) to sequence-specific dsDNA binding to four types of dsDNA: (i) dsDNA with cytosine in both strands ((DNA(C|C)), (ii, iii) dsDNA with 5-methylcytosine (5mC, M) or 5-hydroxymethylcytosine (5hmC, H) in one strand and cytosine in the second strand ((DNA(5mC|C) and DNA(5hmC|C)), and (iv) dsDNA with methylated cytosine in both strands in all CG dinucleotides ((DNA(5mCG)). We replaced asparagine with five similarly sized amino acids (glutamine (Q), serine (S), threonine (T), isoleucine (I), or valine (V)) and used protein binding microarrays to evaluate sequence-specific dsDNA binding. Zta preferentially binds the pseudo-palindrome TRE (AP1) motif (T-4G-3A-2G/C 0T2C3A4 ). Zta (N182Q) changes binding to A3 in only one half-site. Zta(N182S) changes binding to G3 in one or both halves of the motif. Zta(N182S) and Zta(N182Q) have 34- and 17-fold weaker median dsDNA binding, respectively. Zta(N182V) and Zta(N182I) have increased binding to dsDNA(5mC|C). Molecular dynamics simulations rationalize some of these results, identifying hydrogen bonds between glutamine and A3 , but do not reveal why serine preferentially binds G3 , suggesting that entropic interactions may mediate this new binding specificity.
Collapse
Affiliation(s)
- Sreejana Ray
- Laboratory
of Metabolism, National Cancer Institute,
National Institutes of Health, Room 5000, Building 37, Bethesda, Maryland 20892, United States
| | - Desiree Tillo
- Laboratory
of Metabolism, National Cancer Institute,
National Institutes of Health, Room 5000, Building 37, Bethesda, Maryland 20892, United States
- Cancer
Genetics Branch, National Cancer Institute,
National Institutes of Health, Building 37, Bethesda, Maryland 20892, United States
| | - Nima Assad
- Laboratory
of Metabolism, National Cancer Institute,
National Institutes of Health, Room 5000, Building 37, Bethesda, Maryland 20892, United States
| | - Aniekanabasi Ufot
- Laboratory
of Metabolism, National Cancer Institute,
National Institutes of Health, Room 5000, Building 37, Bethesda, Maryland 20892, United States
| | - Aleksey Porollo
- Center
for Autoimmune Genomics and Etiology, Division of Biomedical Informatics, Cincinnati Children’s Hospital Medical Center, Cincinnati, Ohio 45229, United States
- Department
of Pediatrics, University of Cincinnati
College of Medicine, Cincinnati, Ohio 45267, United States
| | - Stewart R. Durell
- Laboratory
of Cell Biology, National Cancer Institute,
National Institutes of Health, Building 37, Bethesda, Maryland 20892, United States
| | - Charles Vinson
- Laboratory
of Metabolism, National Cancer Institute,
National Institutes of Health, Room 5000, Building 37, Bethesda, Maryland 20892, United States
| |
Collapse
|
15
|
Schnepf M, von Reutern M, Ludwig C, Jung C, Gaul U. Transcription Factor Binding Affinities and DNA Shape Readout. iScience 2020; 23:101694. [PMID: 33163946 PMCID: PMC7607496 DOI: 10.1016/j.isci.2020.101694] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2020] [Revised: 09/30/2020] [Accepted: 10/13/2020] [Indexed: 12/16/2022] Open
Abstract
An essential event in gene regulation is the binding of a transcription factor (TF) to its target DNA. Models considering the interactions between the TF and the DNA geometry proved to be successful approaches to describe this binding event, while conserving data interpretability. However, a direct characterization of the DNA shape contribution to binding is still missing due to the lack of accurate and large-scale binding affinity data. Here, we use a binding assay we recently established to measure with high sensitivity the binding specificities of 13 Drosophila TFs, including dinucleotide dependencies to capture non-independent amino acid-base interactions. Correlating the binding affinities with all DNA shape features, we find that shape readout is widely used by these factors. A shape readout/TF-DNA complex structure analysis validates our approach while providing biological insights such as positively charged or highly polar amino acids often contact nucleotides that exhibit strong shape readout. The DNA shape contribution to Drosophila TFs-DNA binding is directly characterized Zeroth- and first-order TF-DNA binding specificities are measured with high accuracy DNA shape readout is widely used by these TFs A shape readout/structural correlation analysis provides biological insights
Collapse
Affiliation(s)
- Max Schnepf
- Gene Center and Department of Biochemistry, Center for Protein Science Munich (CIPSM), Ludwig-Maximilians-Universität München, Feodor-Lynen-Strasse 25, 81377 München, Germany
| | - Marc von Reutern
- Gene Center and Department of Biochemistry, Center for Protein Science Munich (CIPSM), Ludwig-Maximilians-Universität München, Feodor-Lynen-Strasse 25, 81377 München, Germany
| | - Claudia Ludwig
- Gene Center and Department of Biochemistry, Center for Protein Science Munich (CIPSM), Ludwig-Maximilians-Universität München, Feodor-Lynen-Strasse 25, 81377 München, Germany
| | - Christophe Jung
- Gene Center and Department of Biochemistry, Center for Protein Science Munich (CIPSM), Ludwig-Maximilians-Universität München, Feodor-Lynen-Strasse 25, 81377 München, Germany
| | - Ulrike Gaul
- Gene Center and Department of Biochemistry, Center for Protein Science Munich (CIPSM), Ludwig-Maximilians-Universität München, Feodor-Lynen-Strasse 25, 81377 München, Germany
| |
Collapse
|
16
|
Ferguson DCJ, Mokim JH, Meinders M, Moody ERR, Williams TA, Cooke S, Trakarnsanga K, Daniels DE, Ferrer-Vicens I, Shoemark D, Tipgomut C, Macinnes KA, Wilson MC, Singleton BK, Frayne J. Characterization and evolutionary origin of novel C 2H 2 zinc finger protein (ZNF648) required for both erythroid and megakaryocyte differentiation in humans. Haematologica 2020; 106:2859-2873. [PMID: 33054117 PMCID: PMC8561289 DOI: 10.3324/haematol.2020.256347] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2020] [Indexed: 01/01/2023] Open
Abstract
Human ZNF648 is a novel poly C-terminal C2H2 zinc finger protein identified amongst the most dysregulated proteins in erythroid cells differentiated from iPSC. Its nuclear localisation and structure indicate it is likely a DNA-binding protein. Using a combination of ZNF648 overexpression in an iPSC line and primary adult erythroid cells, ZNF648 knockdown in primary adult erythroid cells and megakaryocytes, comparative proteomics and transcriptomics we show that ZNF648 is required for both erythroid and megakaryocyte differentiation. Orthologues of ZNF648 were detected across Mammals, Reptilia, Actinopterygii, in some Aves, Amphibia and Coelacanthiformes suggesting the gene originated in the common ancestor of Osteichthyes (Euteleostomi or bony fish). Conservation of the C-terminal zinc finger domain is higher, with some variation in zinc finger number but a core of at least six zinc fingers conserved across all groups, with the N-terminus recognisably similar within but not between major lineages. This suggests the N-terminus of ZNF648 evolves faster than the C-terminus, however this is not due to exon-shuffling as the entire coding region of ZNF648 is within a single exon. As for other such transcription factors, the N-terminus likely carries out regulatory functions, but showed no sequence similarity to any known domains. The greater functional constraint on the zinc finger domain suggests ZNF648 binds at least some similar regions of DNA in the different organisms. However, divergence of the N-terminal region may enable differential expression, allowing adaptation of function in the different organisms.
Collapse
Affiliation(s)
- Daniel C. J. Ferguson
- School of Biochemistry, University of Bristol, Bristol, UK,*DCJF and JHM contributed equally as co-first authors
| | - Juraidah Haji Mokim
- School of Biochemistry, University of Bristol, Bristol, UK,*DCJF and JHM contributed equally as co-first authors
| | | | | | - Tom A. Williams
- School of Biological Sciences, University of Bristol, Bristol, UK
| | - Sarah Cooke
- School of Biochemistry, University of Bristol, Bristol, UK
| | - Kongtana Trakarnsanga
- Department of Biochemistry, Faculty of Medicine Siriraj Hospital, Mahidol University, Bangkok, Thailand
| | - Deborah E. Daniels
- School of Biochemistry, University of Bristol, Bristol, UK,NIHR Blood and Transplant Research Unit in Red Blood Cell Products, University of Bristol, Bristol, UK
| | | | | | - Chartsiam Tipgomut
- Department of Biochemistry, Faculty of Medicine Siriraj Hospital, Mahidol University, Bangkok, Thailand
| | - Katherine A. Macinnes
- School of Biochemistry, University of Bristol, Bristol, UK,NIHR Blood and Transplant Research Unit in Red Blood Cell Products, University of Bristol, Bristol, UK
| | | | - Belinda K. Singleton
- NIHR Blood and Transplant Research Unit in Red Blood Cell Products, University of Bristol, Bristol, UK,Bristol Institute for Transfusion Sciences, National Health Service Blood and Transplant (NHSBT), Bristol, UK
| | - Jan Frayne
- School of Biochemistry, University of Bristol, BS8 1TD, UK.; NIHR Blood and Transplant Research Unit in Red blood cell products, University of Bristol, Bristol BS8 1TD, UK.
| |
Collapse
|
17
|
Zhang Y, Zheng J. Bioinformatics of Metalloproteins and Metalloproteomes. Molecules 2020; 25:molecules25153366. [PMID: 32722260 PMCID: PMC7435645 DOI: 10.3390/molecules25153366] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2020] [Revised: 07/17/2020] [Accepted: 07/22/2020] [Indexed: 12/14/2022] Open
Abstract
Trace metals are inorganic elements that are required for all organisms in very low quantities. They serve as cofactors and activators of metalloproteins involved in a variety of key cellular processes. While substantial effort has been made in experimental characterization of metalloproteins and their functions, the application of bioinformatics in the research of metalloproteins and metalloproteomes is still limited. In the last few years, computational prediction and comparative genomics of metalloprotein genes have arisen, which provide significant insights into their distribution, function, and evolution in nature. This review aims to offer an overview of recent advances in bioinformatic analysis of metalloproteins, mainly focusing on metalloprotein prediction and the use of different metals across the tree of life. We describe current computational approaches for the identification of metalloprotein genes and metal-binding sites/patterns in proteins, and then introduce a set of related databases. Furthermore, we discuss the latest research progress in comparative genomics of several important metals in both prokaryotes and eukaryotes, which demonstrates divergent and dynamic evolutionary patterns of different metalloprotein families and metalloproteomes. Overall, bioinformatic studies of metalloproteins provide a foundation for systematic understanding of trace metal utilization in all three domains of life.
Collapse
Affiliation(s)
- Yan Zhang
- Shenzhen Key Laboratory of Marine Bioresources and Ecology, College of Life Sciences and Oceanography, Shenzhen University, Shenzhen 518055, China;
- Shenzhen-Hong Kong Institute of Brain Science-Shenzhen Fundamental Research Institutions, Shenzhen 518055, China
- Shenzhen Bay Laboratory, Shenzhen 518055, China
- Correspondence: ; Tel.: +86-755-2692-2024
| | - Junge Zheng
- Shenzhen Key Laboratory of Marine Bioresources and Ecology, College of Life Sciences and Oceanography, Shenzhen University, Shenzhen 518055, China;
- Shenzhen-Hong Kong Institute of Brain Science-Shenzhen Fundamental Research Institutions, Shenzhen 518055, China
- Shenzhen Bay Laboratory, Shenzhen 518055, China
| |
Collapse
|
18
|
Hemming IA, Blake S, Agostino M, Heng JI. General population ZBTB18 missense variants influence DNA binding and transcriptional regulation. Hum Mutat 2020; 41:1629-1644. [DOI: 10.1002/humu.24069] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2020] [Revised: 05/28/2020] [Accepted: 06/22/2020] [Indexed: 12/26/2022]
Affiliation(s)
- Isabel A. Hemming
- The Harry Perkins Institute of Medical Research QEII Medical Centre Nedlands Western Australia Australia
- The Centre for Medical Research The University of Western Australia Crawley Western Australia Australia
- The Faculty of Health and Medical Sciences, Medical School The University of Western Australia Crawley Western Australia Australia
- Ralph and Patricia Sarich Neuroscience Research Institute Nedlands Western Australia Australia
- Curtin Health Innovation Research Institute Curtin University Bentley Western Australia Australia
| | - Steven Blake
- Ralph and Patricia Sarich Neuroscience Research Institute Nedlands Western Australia Australia
- Curtin Health Innovation Research Institute Curtin University Bentley Western Australia Australia
- School of Pharmacy and Biomedical Science Curtin University Bentley Western Australia Australia
| | - Mark Agostino
- Curtin Health Innovation Research Institute Curtin University Bentley Western Australia Australia
- School of Pharmacy and Biomedical Science Curtin University Bentley Western Australia Australia
- Curtin Institute for Computation Curtin University Bentley Western Australia Australia
| | - Julian I‐T. Heng
- Ralph and Patricia Sarich Neuroscience Research Institute Nedlands Western Australia Australia
- Curtin Health Innovation Research Institute Curtin University Bentley Western Australia Australia
| |
Collapse
|
19
|
Mueller AL, Corbi-Verge C, Giganti DO, Ichikawa DM, Spencer JM, MacRae M, Garton M, Kim PM, Noyes MB. The geometric influence on the Cys2His2 zinc finger domain and functional plasticity. Nucleic Acids Res 2020; 48:6382-6402. [PMID: 32383734 PMCID: PMC7293014 DOI: 10.1093/nar/gkaa291] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2020] [Revised: 04/07/2020] [Accepted: 04/20/2020] [Indexed: 11/25/2022] Open
Abstract
The Cys2His2 zinc finger is the most common DNA-binding domain expanding in metazoans since the fungi human split. A proposed catalyst for this expansion is an arms race to silence transposable elements yet it remains poorly understood how this domain is able to evolve the required specificities. Likewise, models of its DNA binding specificity remain error prone due to a lack of understanding of how adjacent fingers influence each other's binding specificity. Here, we use a synthetic approach to exhaustively investigate binding geometry, one of the dominant influences on adjacent finger function. By screening over 28 billion protein–DNA interactions in various geometric contexts we find the plasticity of the most common natural geometry enables more functional amino acid combinations across all targets. Further, residues that define this geometry are enriched in genomes where zinc fingers are prevalent and specificity transitions would be limited in alternative geometries. Finally, these results demonstrate an exhaustive synthetic screen can produce an accurate model of domain function while providing mechanistic insight that may have assisted in the domains expansion.
Collapse
Affiliation(s)
- April L Mueller
- Institute for Systems Genetics, NYU Langone Health, New York, NY 10016, USA.,Department of Biochemistry and Molecular Pharmacology, NYU Langone Health, New York, NY 10016, USA
| | - Carles Corbi-Verge
- Donnelly Centre for Cellular and Biomolecular Research, University of Toronto, Toronto, Ontario M5S 3E1, Canada
| | - David O Giganti
- Institute for Systems Genetics, NYU Langone Health, New York, NY 10016, USA.,Department of Biochemistry and Molecular Pharmacology, NYU Langone Health, New York, NY 10016, USA
| | - David M Ichikawa
- Institute for Systems Genetics, NYU Langone Health, New York, NY 10016, USA.,Department of Biochemistry and Molecular Pharmacology, NYU Langone Health, New York, NY 10016, USA
| | - Jeffrey M Spencer
- Institute for Systems Genetics, NYU Langone Health, New York, NY 10016, USA.,Department of Biochemistry and Molecular Pharmacology, NYU Langone Health, New York, NY 10016, USA
| | - Mark MacRae
- Institute for Systems Genetics, NYU Langone Health, New York, NY 10016, USA.,Department of Biochemistry and Molecular Pharmacology, NYU Langone Health, New York, NY 10016, USA
| | - Michael Garton
- Institute of Biomaterials and Biomedical Engineering, University of Toronto, Toronto, Ontario M5S 3G9, Canada
| | - Philip M Kim
- Donnelly Centre for Cellular and Biomolecular Research, University of Toronto, Toronto, Ontario M5S 3E1, Canada.,Department of Molecular Genetics, University of Toronto, Toronto, Ontario M5S3E1, Canada.,Department of Computer Science, University of Toronto, Toronto, Ontario M5S3E1, Canada
| | - Marcus B Noyes
- Institute for Systems Genetics, NYU Langone Health, New York, NY 10016, USA.,Department of Biochemistry and Molecular Pharmacology, NYU Langone Health, New York, NY 10016, USA
| |
Collapse
|
20
|
Bruno M, Mahgoub M, Macfarlan TS. The Arms Race Between KRAB–Zinc Finger Proteins and Endogenous Retroelements and Its Impact on Mammals. Annu Rev Genet 2019; 53:393-416. [DOI: 10.1146/annurev-genet-112618-043717] [Citation(s) in RCA: 57] [Impact Index Per Article: 11.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Nearly half of the human genome consists of endogenous retroelements (EREs) and their genetic remnants, a small fraction of which carry the potential to propagate in the host genome, posing a threat to genome integrity and cell/organismal survival. The largest family of transcription factors in tetrapods, the Krüppel-associated box domain zinc finger proteins (KRAB-ZFPs), binds to specific EREs and represses their transcription. Since their first appearance over 400 million years ago, KRAB-ZFPs have undergone dramatic expansion and diversification in mammals, correlating with the invasions of new EREs. In this article we review our current understanding of the structure, function, and evolution of KRAB-ZFPs and discuss growing evidence that the arms race between KRAB-ZFPs and the EREs they target is a major driving force for the evolution of new traits in mammals, often accompanied by domestication of EREs themselves.
Collapse
Affiliation(s)
- Melania Bruno
- Eunice Kennedy Shriver National Institute of Child Health and Human Development, The National Institutes of Health, Bethesda, Maryland 20892, USA
| | - Mohamed Mahgoub
- Eunice Kennedy Shriver National Institute of Child Health and Human Development, The National Institutes of Health, Bethesda, Maryland 20892, USA
| | - Todd S. Macfarlan
- Eunice Kennedy Shriver National Institute of Child Health and Human Development, The National Institutes of Health, Bethesda, Maryland 20892, USA
| |
Collapse
|
21
|
Abstract
Evolvability is the ability of a biological system to produce phenotypic variation that is both heritable and adaptive. It has long been the subject of anecdotal observations and theoretical work. In recent years, however, the molecular causes of evolvability have been an increasing focus of experimental work. Here, we review recent experimental progress in areas as different as the evolution of drug resistance in cancer cells and the rewiring of transcriptional regulation circuits in vertebrates. This research reveals the importance of three major themes: multiple genetic and non-genetic mechanisms to generate phenotypic diversity, robustness in genetic systems, and adaptive landscape topography. We also discuss the mounting evidence that evolvability can evolve and the question of whether it evolves adaptively.
Collapse
|
22
|
Hemming IA, Clément O, Gladwyn-Ng IE, Cullen HD, Ng HL, See HB, Ngo L, Ulgiati D, Pfleger KDG, Agostino M, Heng JIT. Disease-associated missense variants in ZBTB18 disrupt DNA binding and impair the development of neurons within the embryonic cerebral cortex. Hum Mutat 2019; 40:1841-1855. [PMID: 31112317 DOI: 10.1002/humu.23803] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2018] [Revised: 05/08/2019] [Accepted: 05/19/2019] [Indexed: 12/12/2022]
Abstract
The activities of DNA-binding transcription factors, such as the multi-zinc-finger protein ZBTB18 (also known as RP58, or ZNF238), are essential to coordinate mammalian neurodevelopment, including the birth and radial migration of newborn neurons within the fetal brain. In humans, the majority of disease-associated missense mutations in ZBTB18 lie within the DNA-binding zinc-finger domain and are associated with brain developmental disorder, yet the molecular mechanisms explaining their role in disease remain unclear. To address this, we developed in silico models of ZBTB18, bound to DNA, and discovered that half of the missense variants map to residues (Asn461, Arg464, Glu486) predicted to be essential to sequence-specific DNA contact, whereas others map to residues (Leu434, Tyr447, Arg495) with limited contributions to DNA binding. We studied pathogenic variants to residues with close (N461S) and limited (R495G) DNA contact and found that each bound DNA promiscuously, displayed altered transcriptional regulatory activity in vitro, and influenced the radial migration of newborn neurons in vivo in different ways. Taken together, our results suggest that altered transcriptional regulation could represent an important pathological mechanism for ZBTB18 missense variants in brain developmental disease.
Collapse
Affiliation(s)
- Isabel A Hemming
- Molecular Medicine Division, QEII Medical Centre, The Harry Perkins Institute of Medical Research, Nedlands, Australia.,Centre for Medical Research, The University of Western Australia, Crawley, Australia.,Faculty of Health and Medical Sciences, Medical School, The University of Western Australia, Crawley, Australia
| | - Olivier Clément
- Molecular Medicine Division, QEII Medical Centre, The Harry Perkins Institute of Medical Research, Nedlands, Australia.,Centre for Medical Research, The University of Western Australia, Crawley, Australia.,Faculty of Health and Medical Sciences, Medical School, The University of Western Australia, Crawley, Australia.,School of Biomedical Sciences, Faculty of Health and Medical Sciences, The University of Western Australia, Crawley, Australia
| | - Ivan E Gladwyn-Ng
- Molecular Medicine Division, QEII Medical Centre, The Harry Perkins Institute of Medical Research, Nedlands, Australia.,Centre for Medical Research, The University of Western Australia, Crawley, Australia
| | - Hayley D Cullen
- Molecular Medicine Division, QEII Medical Centre, The Harry Perkins Institute of Medical Research, Nedlands, Australia.,Centre for Medical Research, The University of Western Australia, Crawley, Australia.,Faculty of Health and Medical Sciences, Medical School, The University of Western Australia, Crawley, Australia.,Curtin Health Innovation Research Institute, Curtin University, Bentley, Australia
| | - Han Leng Ng
- School of Biomedical Sciences, Faculty of Health and Medical Sciences, The University of Western Australia, Crawley, Australia
| | - Heng B See
- Molecular Medicine Division, QEII Medical Centre, The Harry Perkins Institute of Medical Research, Nedlands, Australia.,Centre for Medical Research, The University of Western Australia, Crawley, Australia
| | - Linh Ngo
- Molecular Medicine Division, QEII Medical Centre, The Harry Perkins Institute of Medical Research, Nedlands, Australia.,Centre for Medical Research, The University of Western Australia, Crawley, Australia
| | - Daniela Ulgiati
- School of Biomedical Sciences, Faculty of Health and Medical Sciences, The University of Western Australia, Crawley, Australia
| | - Kevin D G Pfleger
- Molecular Medicine Division, QEII Medical Centre, The Harry Perkins Institute of Medical Research, Nedlands, Australia.,Centre for Medical Research, The University of Western Australia, Crawley, Australia.,ARC Centre for Personalised Therapeutics Technologies, Melbourne, Australia
| | - Mark Agostino
- Curtin Health Innovation Research Institute, Curtin University, Bentley, Australia.,School of Pharmacy and Biomedical Sciences, Curtin University, Bentley, Australia.,Curtin Institute for Computation, Curtin University, Bentley, Australia
| | - Julian I-T Heng
- Curtin Health Innovation Research Institute, Curtin University, Bentley, Australia.,ARC Centre for Personalised Therapeutics Technologies, Melbourne, Australia.,School of Pharmacy and Biomedical Sciences, Curtin University, Bentley, Australia
| |
Collapse
|
23
|
Lambert SA, Yang AWH, Sasse A, Cowley G, Albu M, Caddick MX, Morris QD, Weirauch MT, Hughes TR. Similarity regression predicts evolution of transcription factor sequence specificity. Nat Genet 2019; 51:981-989. [PMID: 31133749 DOI: 10.1038/s41588-019-0411-1] [Citation(s) in RCA: 65] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2018] [Accepted: 04/04/2019] [Indexed: 11/09/2022]
Abstract
Transcription factor (TF) binding specificities (motifs) are essential for the analysis of gene regulation. Accurate prediction of TF motifs is critical, because it is infeasible to assay all TFs in all sequenced eukaryotic genomes. There is ongoing controversy regarding the degree of motif diversification among related species that is, in part, because of uncertainty in motif prediction methods. Here we describe similarity regression, a significantly improved method for predicting motifs, which we use to update and expand the Cis-BP database. Similarity regression inherently quantifies TF motif evolution, and shows that previous claims of near-complete conservation of motifs between human and Drosophila are inflated, with nearly half of the motifs in each species absent from the other, largely due to extensive divergence in C2H2 zinc finger proteins. We conclude that diversification in DNA-binding motifs is pervasive, and present a new tool and updated resource to study TF diversity and gene regulation across eukaryotes.
Collapse
Affiliation(s)
- Samuel A Lambert
- Department of Molecular Genetics, University of Toronto, Toronto, Ontario, Canada
| | - Ally W H Yang
- Donnelly Centre for Cellular and Biomolecular Research, University of Toronto, Toronto, Ontario, Canada
| | - Alexander Sasse
- Department of Molecular Genetics, University of Toronto, Toronto, Ontario, Canada
| | - Gwendolyn Cowley
- Institute of Integrative Biology, University of Liverpool, Liverpool, UK
| | - Mihai Albu
- Donnelly Centre for Cellular and Biomolecular Research, University of Toronto, Toronto, Ontario, Canada
| | - Mark X Caddick
- Institute of Integrative Biology, University of Liverpool, Liverpool, UK
| | - Quaid D Morris
- Department of Molecular Genetics, University of Toronto, Toronto, Ontario, Canada.,Donnelly Centre for Cellular and Biomolecular Research, University of Toronto, Toronto, Ontario, Canada.,Department of Computer Science, University of Toronto, Toronto, Ontario, Canada.,Canadian Institutes For Advanced Research (CIFAR) Artificial Intelligence Chair, Vector Institute, Toronto, Ontario, Canada.,Ontario Institute of Cancer Research, Toronto, Ontario, Canada
| | - Matthew T Weirauch
- Divisions of Biomedical Informatics and Developmental Biology, Center for Autoimmune Genomics and Etiology (CAGE), Cincinnati Children's Hospital Medical Center, Cincinnati, OH, USA.,Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati, OH, USA
| | - Timothy R Hughes
- Department of Molecular Genetics, University of Toronto, Toronto, Ontario, Canada. .,Donnelly Centre for Cellular and Biomolecular Research, University of Toronto, Toronto, Ontario, Canada. .,CIFAR, Toronto, Ontario, Canada.
| |
Collapse
|
24
|
Ray S, Ufot A, Assad N, Singh J, Durell SR, Porollo A, Tillo D, Vinson C. The bZIP mutant CEBPB (V285A) has sequence specific DNA binding propensities similar to CREB1. BIOCHIMICA ET BIOPHYSICA ACTA-GENE REGULATORY MECHANISMS 2019; 1862:486-492. [PMID: 30825655 DOI: 10.1016/j.bbagrm.2019.02.002] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/03/2018] [Revised: 01/09/2019] [Accepted: 02/05/2019] [Indexed: 12/25/2022]
Abstract
The bZIP homodimers CEBPB and CREB1 bind DNA containing methylated cytosines differently. CREB1 binds stronger to the C/EBP half-site GCAA when the cytosine is methylated. For CEBPB, methylation of the same cytosine does not affect DNA binding. The X-ray structure of CREB1 binding the half site GTCA identifies an alanine in the DNA binding region interacting with the methyl group of T, structurally analogous to the methyl group of methylated C. This alanine is replaced with a valine in CEBPB. To explore the contribution of this amino acid to binding with methylated cytosine of the GCAA half-site, we made the reciprocal mutants CEBPB(V285A) and CREB1(A297V) and used protein binding microarrays (PBM) to examine binding to four types of double-stranded DNA (dsDNA): 1) DNA with cytosine in both strands (DNA(C|C)), 2) DNA with 5-methylcytosine (M) in one strand and cytosine in the second strand (DNA(M|C)), 3) DNA with 5-hydroxymethylcytosine (H) in one strand and cytosine in the second strand (DNA(H|C)), and 4) DNA with both cytosines in all CG dinucleotides containing 5-methylcytosine (DNA(5mCG)). When binding to DNA(C|C), CEBPB (V285A) preferentially binds the CRE consensus motif (TGACGTCA), similar to CREB1. The reciprocal mutant, CREB1(A297V) binds DNA with some similarity to CEBPB, with strongest binding to the methylated PAR site 8-mer TTACGTAA. These data demonstrate that V285 residue inhibits CEBPB binding to methylated cytosine of the GCAA half-site.
Collapse
Affiliation(s)
- Sreejana Ray
- Laboratory of Metabolism, National Cancer Institute, National Institutes of Health, Bethesda, MD 20892, United States of America
| | - Aniekanabasi Ufot
- Laboratory of Metabolism, National Cancer Institute, National Institutes of Health, Bethesda, MD 20892, United States of America
| | - Nima Assad
- Laboratory of Metabolism, National Cancer Institute, National Institutes of Health, Bethesda, MD 20892, United States of America
| | - Jocelyn Singh
- Laboratory of Metabolism, National Cancer Institute, National Institutes of Health, Bethesda, MD 20892, United States of America
| | - Stewart R Durell
- Laboratory of Cell Biology, National Cancer Institute, National Institutes of Health, Bethesda, MD 20892, United States of America
| | - Aleksey Porollo
- Center for Autoimmune Genomics and Etiology, Division of Biomedical Informatics, Cincinnati Children's Hospital Medical Center, Cincinnati, OH 45229, United States of America; Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati, OH 45267, United States of America
| | - Desiree Tillo
- Laboratory of Metabolism, National Cancer Institute, National Institutes of Health, Bethesda, MD 20892, United States of America
| | - Charles Vinson
- Laboratory of Metabolism, National Cancer Institute, National Institutes of Health, Bethesda, MD 20892, United States of America.
| |
Collapse
|
25
|
Zhang Y, Ying H, Xu Y. Comparative genomics and metagenomics of the metallomes. Metallomics 2019; 11:1026-1043. [DOI: 10.1039/c9mt00023b] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Recent achievements and advances in comparative genomic and metagenomic analyses of trace metals were reviewed and discussed.
Collapse
Affiliation(s)
- Yan Zhang
- Shenzhen Key Laboratory of Marine Bioresources and Ecology
- College of Life Sciences and Oceanography
- Shenzhen University
- Shenzhen
- P. R. China
| | - Huimin Ying
- Department of Endocrinology
- Hangzhou Xixi Hospital
- Hangzhou
- P. R. China
| | - Yinzhen Xu
- Shenzhen Key Laboratory of Marine Bioresources and Ecology
- College of Life Sciences and Oceanography
- Shenzhen University
- Shenzhen
- P. R. China
| |
Collapse
|
26
|
Shi P, Zhang Y, Li Y, Bian L. Probing the interaction of l
-captopril with metallo-β-lactamase CcrA by fluorescence spectra and molecular dynamic simulation. LUMINESCENCE 2018; 33:954-961. [DOI: 10.1002/bio.3495] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2017] [Revised: 02/18/2018] [Accepted: 04/03/2018] [Indexed: 01/01/2023]
Affiliation(s)
- Penghui Shi
- College of Life Science; Northwest University; Xi'an Shaanxi People's Republic of China
| | - Yan Zhang
- College of Life Science; Northwest University; Xi'an Shaanxi People's Republic of China
| | - Yuhua Li
- College of Life Science; Northwest University; Xi'an Shaanxi People's Republic of China
| | - Liujiao Bian
- College of Life Science; Northwest University; Xi'an Shaanxi People's Republic of China
| |
Collapse
|
27
|
Doğan B, Najafabadi HS. Computational Methods for Analysis of the DNA-Binding Preferences of Cys2His2 Zinc-Finger Proteins. Methods Mol Biol 2018; 1867:15-28. [PMID: 30155812 DOI: 10.1007/978-1-4939-8799-3_2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Cys2His2 zinc-finger proteins (C2H2-ZFPs) constitute the largest class of human transcription factors (TFs) and also the least characterized one. Determining the DNA sequence preferences of C2H2-ZFPs is an important first step toward elucidating their roles in transcriptional regulation. Among the most promising approaches for obtaining the sequence preferences of C2H2-ZFPs are those that combine machine-learning predictions with in vivo binding maps of these proteins. Here, we provide a protocol and guidelines for predicting the DNA-binding preferences of C2H2-ZFPs from their amino acid sequences using a machine learning-based recognition code. This protocol also describes the tools and steps to combine these predictions with ChIP-seq data to remove inaccuracies, identify the zinc-finger domains within each C2H2-ZFP that engage with DNA in vivo, and pinpoint the genomic binding sites of the C2H2-ZFPs.
Collapse
Affiliation(s)
- Berat Doğan
- Department of Human Genetics, McGill University, Montreal, QC, Canada
- McGill University and Genome Quebec Innovation Centre, Montreal, QC, Canada
- Department of Biomedical Engineering, İnönü University, Malatya, Turkey
| | - Hamed S Najafabadi
- Department of Human Genetics, McGill University, Montreal, QC, Canada.
- McGill University and Genome Quebec Innovation Centre, Montreal, QC, Canada.
| |
Collapse
|
28
|
|