Comparative genomics of the Hedgehog loci in chordates and the origins of Shh regulatory novelties

The origin and evolution of the complex regulatory landscapes of some vertebrate developmental genes, often spanning hundreds of Kbp and including neighboring genes, remain poorly understood. The Sonic Hedgehog (Shh) genomic regulatory block (GRB) is one of the best functionally characterized examples, with several discrete enhancers reported within its introns, vast upstream gene-free region and neighboring genes (Lmbr1 and Rnf32). To investigate the origin and evolution of this GRB, we sequenced and characterized the Hedgehog (Hh) loci from three invertebrate chordate amphioxus species, which share several early expression domains with Shh. Using phylogenetic footprinting within and between chordate lineages, and reporter assays in zebrafish probing >30 Kbp of amphioxus Hh, we report large sequence and functional divergence between both groups. In addition, we show that the linkage of Shh to Lmbr1 and Rnf32, necessary for the unique gnatostomate-specific Shh limb expression, is a vertebrate novelty occurred between the two whole-genome duplications

Extensive conservation of ancient microsynteny across metazoans due to cis-regulatory constraints

This article is distributed exclusively by Cold Spring Harbor Laboratory Press for the first six months after the full-issue publication date; after six months, it is available under a Creative Commons License.-- et al.The order of genes in eukaryotic genomes has generally been assumed to be neutral, since gene order is largely scrambled over evolutionary time. Only a handful of exceptional examples are known, typically involving deeply conserved clusters of tandemly duplicated genes (e.g., Hox genes and histones). Here we report the first systematic survey of microsynteny conservation across metazoans, utilizing 17 genome sequences. We identified nearly 600 pairs of unrelated genes that have remained tightly physically linked in diverse lineages across over 600 million years of evolution. Integrating sequence conservation, gene expression data, gene function, epigenetic marks, and other genomic features, we provide extensive evidence that many conserved ancient linkages involve (1) the coordinated transcription of neighboring genes, or (2) genomic regulatory blocks (GRBs) in which transcriptional enhancers controlling developmental genes are contained within nearby bystander genes. In addition, we generated ChIP-seq data for key histone modifications in zebrafish embryos, which provided further evidence of putative GRBs in embryonic development. Finally, using chromosome conformation capture (3C) assays and stable transgenic experiments, we demonstrate that enhancers within bystander genes drive the expression of genes such as Otx and Islet, critical regulators of central nervous system development across bilaterians. These results suggest that ancient genomic functional associations are far more common than previously thought—involving ∼12% of the ancestral bilaterian genome—and that cis-regulatory constraints are crucial in determining metazoan genome architecture.M.I., M.S.A., S.W.R., and H.B.F. were funded by NIH grant 1R21HG005240-01A1. H.B.F. is an Alfred P. Sloan Fellow and Pew Scholar in the Biomedical Sciences. J.J.T., A.F-M., O.B., E.C-M., and J.L.G-S. were funded by grants BFU2010-14839, CSD2007-00008, and Proyecto de Excelencia CVI-3488.Peer reviewe

Regulatory architecture of the RCA gene cluster captures an intragenic TAD boundary, CTCF-Mediated chromatin looping and a long-range intergenic enhancer

The Regulators of Complement Activation (RCA) gene cluster comprises several tandemly arranged genes with shared functions within the immune system. RCA members, such as complement receptor 2 (CR2), are well-established susceptibility genes in complex autoimmune diseases. Altered expression of RCA genes has been demonstrated at both the functional and genetic level, but the mechanisms underlying their regulation are not fully characterised. We aimed to investigate the structural organisation of the RCA gene cluster to identify key regulatory elements that influence the expression of CR2 and other genes in this immunomodulatory region. Using 4C, we captured extensive CTCF-mediated chromatin looping across the RCA gene cluster in B cells and showed these were organised into two topologically associated domains (TADs). Interestingly, an inter-TAD boundary was located within the CR1 gene at a well-characterised segmental duplication. Additionally, we mapped numerous gene-gene and gene-enhancer interactions across the region, revealing extensive co-regulation. Importantly, we identified an intergenic enhancer and functionally demonstrated this element upregulates two RCA members (CR2 and CD55) in B cells. We have uncovered novel, long-range mechanisms whereby autoimmune disease susceptibility may be influenced by genetic variants, thus highlighting the important contribution of chromatin topology to gene regulation and complex genetic disease.This work was supported by the National Institutes of Health [R01 AI24717 to JH], the Australian Government Research Training Program Scholarship at the University of Western Australia [to JC and JSC], the Spanish Government [BFU2016-74961-P to JG-S] and an institutional grant Unidad de Excelencia María de Maeztu [MDM-206-0687 to the Department of Gene Regulation and Morphogenesis, Centro Andaluz de Biología del Desarrol]

The role of Cis-Regulatory elements in morphological adaptation to cave environment in Astyanax mexicanus

Trabajo presentado en EMBO Workshop Enhanceropathies: Understanding enhancer function to understand human disease, celebrado en Santander (España) del 06 al 09 de octubre de 2021

Genome evolution in morphological adaptation to cave environment in Astyanax mexicanus

Trabajo presentado en el EMBO Workshop The evolution of animal genomes, celebrado en modalidad virtual del 13 al 17 de septiembre de 2021

Asthma-associated genetic variants induce IL33 differential expression through an enhancer-blocking regulatory region

Genome-wide association studies (GWAS) have implicated the IL33 locus in asthma, but the underlying mechanisms remain unclear. Here, we identify a 5 kb region within the GWAS-defined segment that acts as an enhancer-blocking element in vivo and in vitro. Chromatin conformation capture showed that this 5 kb region loops to the IL33 promoter, potentially regulating its expression. We show that the asthma-associated single nucleotide polymorphism (SNP) rs1888909, located within the 5 kb region, is associated with IL33 gene expression in human airway epithelial cells and IL-33 protein expression in human plasma, potentially through differential binding of OCT-1 (POU2F1) to the asthma-risk allele. Our data demonstrate that asthma-associated variants at the IL33 locus mediate allele-specific regulatory activity and IL33 expression, providing a mechanism through which a regulatory SNP contributes to genetic risk of asthma.This work was supported by NIH grants R01 HL118758, R01 HL128075, R01 HL119577, R01 HL085197, U19 AI095230, UG3 OD023282 and UM1 AI114271

Dissecting the Transcriptional Regulatory Properties of Human Chromosome 16 Highly Conserved Non-Coding Regions

Non-coding DNA conservation across species has been often used as a predictor for transcriptional enhancer activity. However, only a few systematic analyses of the function of these highly conserved non-coding regions (HCNRs) have been performed. Here we use zebrafish transgenic assays to perform a systematic study of 113 HCNRs from human chromosome 16. By comparing transient and stable transgenesis, we show that the first method is highly inefficient, leading to 40% of false positives and 20% of false negatives. When analyzed in stable transgenic lines, a great majority of HCNRs were active in the central nervous system, although some of them drove expression in other organs such as the eye and the excretory system. Finally, by testing a fraction of the HCNRs lacking enhancer activity for in vivo insulator activity, we find that 20% of them may contain enhancer-blocking function. Altogether our data indicate that HCNRs may contain different types of cis-regulatory activity, including enhancer, insulators as well as other not yet discovered functions

3D genomics across the tree of life reveals condensin II as a determinant of architecture type

We investigated genome folding across the eukaryotic tree of life. We find two types of three-dimensional(3D) genome architectures at the chromosome scale. Each type appears and disappears repeatedlyduring eukaryotic evolution. The type of genome architecture that an organism exhibits correlates with theabsence of condensin II subunits. Moreover, condensin II depletion converts the architecture of thehuman genome to a state resembling that seen in organisms such as fungi or mosquitoes. In this state,centromeres cluster together at nucleoli, and heterochromatin domains merge. We propose a physicalmodel in which lengthwise compaction of chromosomes by condensin II during mitosis determineschromosome-scale genome architecture, with effects that are retained during the subsequent interphase.This mechanism likely has been conserved since the last common ancestor of all eukaryotes.C.H. is supported by the Boehringer Ingelheim Fonds; C.H., Á.S.C., and B.D.R. are supported by an ERC CoG (772471, “CohesinLooping”); A.M.O.E. and B.D.R. are supported by the Dutch Research Council (NWO-Echo); and J.A.R. and R.H.M. are supported by the Dutch Cancer Society (KWF). T.v.S. and B.v.S. are supported by NIH Common Fund “4D Nucleome” Program grant U54DK107965. H.T. and E.d.W. are supported by an ERC StG (637597, “HAP-PHEN”). J.A.R., T.v.S., H.T., R.H.M., B.v.S., and E.d.W. are part of the Oncode Institute, which is partly financed by the Dutch Cancer Society. Work at the Center for Theoretical Biological Physics is sponsored by the NSF (grants PHY-2019745 and CHE-1614101) and by the Welch Foundation (grant C-1792). V.G.C. is funded by FAPESP (São Paulo State Research Foundation and Higher Education Personnel) grants 2016/13998-8 and 2017/09662-7. J.N.O. is a CPRIT Scholar in Cancer Research. E.L.A. was supported by an NSF Physics Frontiers Center Award (PHY-2019745), the Welch Foundation (Q-1866), a USDA Agriculture and Food Research Initiative grant (2017-05741), the Behavioral Plasticity Research Institute (NSF DBI-2021795), and an NIH Encyclopedia of DNA Elements Mapping Center Award (UM1HG009375). Hi-C data for the 24 species were created by the DNA Zoo Consortium (www.dnazoo.org). DNA Zoo is supported by Illumina, Inc.; IBM; and the Pawsey Supercomputing Center. P.K. is supported by the University of Western Australia. L.L.M. was supported by NIH (1R01NS114491) and NSF awards (1557923, 1548121, and 1645219) and the Human Frontiers Science Program (RGP0060/2017). The draft A. californica project was supported by NHGRI. J.L.G.-S. received funding from the ERC (grant agreement no. 740041), the Spanish Ministerio de Economía y Competitividad (grant no. BFU2016-74961-P), and the institutional grant Unidad de Excelencia María de Maeztu (MDM-2016-0687). R.D.K. is supported by NIH grant RO1DK121366. V.H. is supported by NIH grant NIH1P41HD071837. K.M. is supported by a MEXT grant (20H05936). M.C.W. is supported by the NIH grants R01AG045183, R01AT009050, R01AG062257, and DP1DK113644 and by the Welch Foundation. E.F. was supported by NHGR

SIRPB1 copy-number polymorphism as candidate quantitative trait locus for impulsive-disinhibited personality

Impulsive-disinhibited personality (IDP) is a behavioral trait mainly characterized by seeking immediate gratification at the expense of more enduring or long-term gains. This trait has a major role in the development of several disinhibitory behaviors and syndromes, including psychopathy, attention-deficit and hyperactivity disorder, cluster-B personality disorders, criminality and alcoholism. Available data consistently support a strong heritable component, accounting for 30-60% of the observed variance in personality traits. A genome-wide analysis of copy-number variants was designed to identify novel genetic pathways associated with the IDP trait, using a series of 261 male participants with maximized opposite IDP scores. Quantitative trait locus analysis of candidate copy-number variants (CNVs) was conducted across the entire IDP continuum. Functional effects of associated variants were evaluated in zebrafish embryos. A common CNV mapping to the immune-related gene SIRPB1 was significantly associated with IDP scores in a dose-dependent manner (β=-0.172, P<0.017). Expression quantitative trait locus analysis of the critical region revealed higher SIRPB1 mRNA levels associated with the haplotype containing the deleted allele (P<0.0007). Epigenetic marks highlighted the presence of two potential insulators within the deleted region, confirmed by functional assays in zebrafish embryos, which suggests that SIRPB1 expression rates are affected by the presence/absence of the insulator regions. Upregulation of SIRPB1 has been described in prefrontal cortex of patients with schizophrenia, providing a link between SIRPB1 and diseases involving disinhibition and failure to control impulsivity. We propose SIRPB1 as a novel candidate gene to account for phenotypic differences observed in the IDP trait.This study was founded by Fundació Marató TV3 (ref 200925-10).Peer Reviewe