Silvia Naranjo

Obesity-associated variants within FTO form long-range functional connections with IRX3

Genome-wide association studies (GWAS) have reproducibly associated variants within introns of FTO with increased risk for obesity and type 2 diabetes (T2D). Although the molecular mechanisms linking these noncoding variants with obesity are not immediately obvious, subsequent studies in mice demonstrated that FTO expression levels influence body mass and composition phenotypes. However, no direct connection between the obesity-associated variants and FTO expression or function has been made. Here we show that the obesity-associated noncoding sequences within FTO are functionally connected, at megabase distances, with the homeobox gene IRX3. The obesity-associated FTO region directly interacts with the promoters of IRX3 as well as FTO in the human, mouse and zebrafish genomes. Furthermore, long-range enhancers within this region recapitulate aspects of IRX3 expression, suggesting that the obesity-associated interval belongs to the regulatory landscape of IRX3. Consistent with this, obesity-associated single nucleotide polymorphisms are associated with expression of IRX3, but not FTO, in human brains. A direct link between IRX3 expression and regulation of body mass and composition is demonstrated by a reduction in body weight of 25 to 30% in Irx3-deficient mice, primarily through the loss of fat mass and increase in basal metabolic rate with browning of white adipose tissue. Finally, hypothalamic expression of a dominant-negative form of Irx3 reproduces the metabolic phenotypes of Irx3-deficient mice. Our data suggest that IRX3 is a functional long-range target of obesity-associated variants within FTO and represents a novel determinant of body mass and composition.

Zebrafish enhancer detection (ZED) vector: A new tool to facilitate transgenesis and the functional analysis of cis-regulatory regions in zebrafish

The identification and characterization of the regulatory activity of genomic sequences is crucial for understanding how the information contained in genomes is translated into cellular function. The cis‐regulatory sequences control when, where, and how much genes are transcribed and can activate (enhancers) or repress (silencers) gene expression. Here, we describe a novel Tol2 transposon‐based vector for assessing enhancer activity in the zebrafish (Danio rerio). This Zebrafish Enhancer Detector (ZED) vector harbors several key improvements, among them a sensitive and specific minimal promoter chosen for optimal enhancer activity detection, insulator sequences to shield the minimal promoter from position effects, and a positive control for transgenesis. Additionally, we demonstrate that highly conserved noncoding sequences homologous between humans and zebrafish largely with enhancer activity largely retain their tissue‐specific enhancer activity during vertebrate evolution. More strikingly, insulator sequences from mouse and chicken, but not conserved in zebrafish, maintain their insulator capacity when tested in this model. Developmental Dynamics 238:2409–2417, 2009.

The colorectal cancer risk at 18q21 is caused by a novel variant altering SMAD7 expression

Recent genome-wide scans for colorectal cancer (CRC) have revealed the SMAD7 (mothers against decapentaplegic homolog 7) gene as a locus associated with a modest, but highly significant increase in CRC risk. To identify the causal basis of the association between 18q21 variation and CRC, we resequenced the 17-kb region of linkage disequilibrium and evaluated all variants in 2532 CRC cases and 2607 controls. A novel C to G single nucleotide polymorphism (SNP) at 44,703,563 bp was maximally associated with CRC risk (P = 5.98 x 10(-7); > or =1.5-fold more likely to be causal than other variants). Using transgenic assays in Xenopus laevis as a functional model, we demonstrate that the G risk allele leads to reduced reporter gene expression in the colorectum (P = 5.4 x 10(-3)). Electrophoretic mobility shift assays provided evidence for the role of Novel 1 in transcription factor binding. We propose that the novel SNP we have identified is the functional change leading to CRC predisposition through differential SMAD7 expression and, hence, aberrant TGF-beta signaling.

Allelic Variation at the 8q23.3 Colorectal Cancer Risk Locus Functions as a Cis-Acting Regulator of EIF3H

Common genetic variation at human 8q23.3 is significantly associated with colorectal cancer (CRC) risk. To elucidate the basis of this association we compared the frequency of common variants at 8q23.3 in 1,964 CRC cases and 2,081 healthy controls. Reporter gene studies showed that the single nucleotide polymorphism rs16888589 acts as an allele-specific transcriptional repressor. Chromosome conformation capture (3C) analysis demonstrated that the genomic region harboring rs16888589 interacts with the promoter of gene for eukaryotic translation initiation factor 3, subunit H (EIF3H). We show that increased expression of EIF3H gene increases CRC growth and invasiveness thereby providing a biological mechanism for the 8q23.3 association. These data provide evidence for a functional basis for the non-coding risk variant rs16888589 at 8q23.3 and provides novel insight into the etiological basis of CRC.

Evolutionary comparison reveals that diverging CTCF sites are signatures of ancestral topological associating domains borders

Significance Mammalian chromatin is compartmentalized in topologically associating domains (TADs), genomic regions within which sequences preferentially contact each other. This organization has been proposed to be essential to organize the regulatory information contained in mammalian genomes. We show that Six homeobox genes, essential developmental regulators organized in gene clusters across different animal phyla, share a deeply conserved chromatin organization formed by two abutting TADs that predates the Cambrian explosion. This organization is required to generate separate regulatory landscapes for neighboring genes within the cluster, resulting in very different gene expression patterns. Finally, we show that this extremely conserved 3D architecture is associated with a characteristic arrangement of CCCTC-binding factor (CTCF) binding sites in diverging orientations, revealing a genome-wide conserved signature for TAD borders. Increasing evidence in the last years indicates that the vast amount of regulatory information contained in mammalian genomes is organized in precise 3D chromatin structures. However, the impact of this spatial chromatin organization on gene expression and its degree of evolutionary conservation is still poorly understood. The Six homeobox genes are essential developmental regulators organized in gene clusters conserved during evolution. Here, we reveal that the Six clusters share a deeply evolutionarily conserved 3D chromatin organization that predates the Cambrian explosion. This chromatin architecture generates two largely independent regulatory landscapes (RLs) contained in two adjacent topological associating domains (TADs). By disrupting the conserved TAD border in one of the zebrafish Six clusters, we demonstrate that this border is critical for preventing competition between promoters and enhancers located in separated RLs, thereby generating different expression patterns in genes located in close genomic proximity. Moreover, evolutionary comparison of Six-associated TAD borders reveals the presence of CCCTC-binding factor (CTCF) sites with diverging orientations in all studied deuterostomes. Genome-wide examination of mammalian HiC data reveals that this conserved CTCF configuration is a general signature of TAD borders, underscoring that common organizational principles underlie TAD compartmentalization in deuterostome evolution.

The 14q22.2 colorectal cancer variant rs4444235 shows cis-acting regulation of BMP4

Common genetic variation at human 14q22.2 tagged by rs4444235 is significantly associated with colorectal cancer (CRC) risk. Re-sequencing was used to comprehensively annotate the 17kb region of strong linkage disequilibrium encompassing rs4444235. Through bioinformatic analyses using H3K4Me1, H3K4Me3, and DNase-I hypersensitivity chromatin signatures and evolutionary conservation we identified seven candidate disease-causing single-nucleotide polymorphisms mapping to six regions within the 17-kb region predicted to have regulatory potential. Reporter gene studies of these regions demonstrated that the element to which rs4444235 maps acts as an allele-specific transcriptional enhancer. Allele-specific expression studies in CRC cell lines heterozygous for rs4444235 showed significantly increased expression of bone morphogenetic protein-4 (BMP4) associated with the risk allele (P<0.001). These data provide evidence for a functional basis for the non-coding risk variant rs4444235 at 14q22.2 and emphasizes the importance of genetic variation in the BMP pathway genes as determinants of CRC risk.

Characterization of New Otic Enhancers of the Pou3f4 Gene Reveal Distinct Signaling Pathway Regulation and Spatio-Temporal Patterns

POU3F4 is a member of the POU-homedomain transcription factor family with a prominent role in inner ear development. Mutations in the human POU3F4 coding unit leads to X-linked deafness type 3 (DFN3), characterized by conductive hearing loss and progressive sensorineural deafness. Microdeletions found 1 Mb 5′ upstream of the coding region also displayed the same phenotype, suggesting that cis-regulatory elements might be present in that region. Indeed, we and others have recently identified several enhancers at the 1 Mb 5′ upstream interval of the pou3f4 locus. Here we characterize the spatio-temporal patterns of these regulatory elements in zebrafish transgenic lines. We show that the most distal enhancer (HCNR 81675) is activated earlier and drives GFP reporter expression initially to a broad ear domain to progressively restrict to the sensory patches. The proximal enhancer (HCNR 82478) is switched later during development and promotes expression, among in other tissues, in sensory patches from its onset. The third enhancer (HCNR 81728) is also active at later stages in the otic mesenchyme and in the otic epithelium. We also characterize the signaling pathways regulating these enhancers. While HCNR 81675 is regulated by very early signals of retinoic acid, HCNR 82478 is regulated by Fgf activity at a later stage and the HCNR 81728 enhancer is under the control of Hh signaling. Finally, we show that Sox2 and Pax2 transcription factors are bound to HCNR 81675 genomic region during otic development and specific mutations to these transcription factor binding sites abrogates HCNR 81675 enhancer activity. Altogether, our results suggest that pou3f4 expression in inner ear might be under the control of distinct regulatory elements that fine-tune the spatio-temporal activity of this gene and provides novel data on the signaling mechanisms controlling pou3f4 function.

A conserved Shh cis-regulatory module highlights a common developmental origin of unpaired and paired fins

Despite their evolutionary, developmental and functional importance, the origin of vertebrate paired appendages remains uncertain. In mice, a single enhancer termed ZRS is solely responsible for Shh expression in limbs. Here, zebrafish and mouse transgenic assays trace the functional equivalence of ZRS across the gnathostome phylogeny. CRISPR/Cas9-mediated deletion of the medaka (Oryzias latipes) ZRS and enhancer assays identify the existence of ZRS shadow enhancers in both teleost and human genomes. Deletion of both ZRS and shadow ZRS abolishes shh expression and completely truncates pectoral fin formation. Strikingly, deletion of ZRS results in an almost complete ablation of the dorsal fin. This finding indicates that a ZRS-Shh regulatory module is shared by paired and median fins and that paired fins likely emerged by the co-option of developmental programs established in the median fins of stem gnathostomes. Shh function was later reinforced in pectoral fin development with the recruitment of shadow enhancers, conferring additional robustness.The authors study the cis-regulatory evolution of the Shh locus in vertebrates. Using genomic editing and chromatin profiling, they conclude that paired fins emerged through the co-option of developmental programs for the median fins of gnathostomes.

Analysis of opo cis-regulatory landscape uncovers Vsx2 requirement in early eye morphogenesis

The self-organized morphogenesis of the vertebrate optic cup entails coupling the activation of the retinal gene regulatory network to the constriction-driven infolding of the retinal epithelium. Yet the genetic mechanisms underlying this coordination remain largely unexplored. Through phylogenetic footprinting and transgenesis in zebrafish, here we examine the cis-regulatory landscape of opo, an endocytosis regulator essential for eye morphogenesis. Among the different conserved enhancers identified, we isolate a single retina-specific element (H6_10137) and show that its activity depends on binding sites for the retinal determinant Vsx2. Gain- and loss-of-function experiments and ChIP analyses reveal that Vsx2 regulates opo expression through direct binding to this retinal enhancer. Furthermore, we show that vsx2 knockdown impairs the primary optic cup folding. These data support a model by which vsx2, operating through the effector gene opo, acts as a central transcriptional node that coordinates neural retina patterning and optic cup invagination in zebrafish.

A mobile insulator system to detect and disrupt cis-regulatory landscapes in vertebrates

In multicellular organisms, cis-regulation controls gene expression in space and time. Despite the essential implication of cis-regulation in the development and evolution of organisms and in human diseases, our knowledge about regulatory sequences largely derives from analyzing their activity individually and outside their genomic context. Indeed, the contribution of these sequences to the expression of their target genes in their genomic context is still largely unknown. Here we present a novel genetic screen designed to visualize and interrupt gene regulatory landscapes in vertebrates. In this screen, based on the random insertion of an engineered Tol2 transposon carrying a strong insulator separating two fluorescent reporter genes, we isolated hundreds of zebrafish lines containing insertions that disrupt the cis-regulation of tissue-specific expressed genes. We therefore provide a new easy-to-handle tool that will help to disrupt and chart the regulatory activity spread through the vast noncoding regions of the vertebrate genome.