MacroH2As regulate enhancer-promoter contacts affecting enhancer activity and sensitivity to inflammatory cytokines.

MacroH2A histone variants have a function in gene regulation that is poorly understood at the molecular level. We report that macroH2A1.2 and macroH2A2 modulate the transcriptional ground state of cancer cells and how they respond to inflammatory cytokines. Removal of macroH2A1.2 and macroH2A2 in hepatoblastoma cells affects the contact frequency of promoters and distal enhancers coinciding with changes in enhancer activity or preceding them in response to the cytokine tumor necrosis factor alpha. Although macroH2As regulate genes in both directions, they globally facilitate the nuclear factor κB (NF-κB)-mediated response. In contrast, macroH2As suppress the response to the pro-inflammatory cytokine interferon gamma. MacroH2A2 has a stronger contribution to gene repression than macroH2A1.2. Taken together, our results suggest that macroH2As have a role in regulating the response of cancer cells to inflammatory signals on the level of chromatin structure. This is likely relevant for the interaction of cancer cells with immune cells of their microenvironment.

The Shh/Gli3 gene regulatory network precedes the origin of paired fins and reveals the deep homology between distal fins and digits.

One of the central problems of vertebrate evolution is understanding the relationship among the distal portions of fins and limbs. Lacking comparable morphological markers of these regions in fish and tetrapods, these relationships have remained uncertain for the past century and a half. Here we show that Gli3 functions in controlling the proliferative expansion of distal progenitors are shared among dorsal and paired fins as well as tetrapod limbs. Mutant knockout gli3 fins in medaka (Oryzias latipes) form multiple radials and rays, in a pattern reminiscent of the polydactyly observed in Gli3-null mutant mice. In limbs, Gli3 controls both anterior-posterior patterning and cell proliferation, two processes that can be genetically uncoupled. In situ hybridization, quantification of proliferation markers, and analysis of regulatory regions reveal that in paired and dorsal fins, gli3 plays a main role in controlling proliferation but not in patterning. Moreover, gli3 down-regulation in shh mutant fins rescues fin loss in a manner similar to how Gli3 deficiency restores digits in the limbs of Shh mutant mouse embryos. We hypothesize that the Gli3/Shh gene pathway preceded the origin of paired appendages and was originally involved in modulating cell proliferation. Accordingly, the distal regions of dorsal fins, paired fins, and limbs retain a deep regulatory and functional homology that predates the origin of paired appendages.

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 repeatedly during eukaryotic evolution. The type of genome architecture that an organism exhibits correlates with the absence of condensin II subunits. Moreover, condensin II depletion converts the architecture of the human 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 physical model in which lengthwise compaction of chromosomes by condensin II during mitosis determines chromosome-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.

Ancient Genomic Regulatory Blocks Are a Source for Regulatory Gene Deserts in Vertebrates after Whole-Genome Duplications.

We investigated how the two rounds of whole-genome duplication that occurred at the base of the vertebrate lineage have impacted ancient microsyntenic associations involving developmental regulators (known as genomic regulatory blocks, GRBs). We showed that the majority of GRBs identified in the last common ancestor of chordates have been maintained as a single copy in humans. We found evidence that dismantling of the duplicated GRB copies occurred early in vertebrate evolution often through the differential retention of the regulatory gene but loss of the bystander gene's exonic sequences. Despite the large evolutionary scale, the presence of duplicated highly conserved noncoding regions provided unambiguous proof for this scenario for multiple ancient GRBs. Remarkably, the dismantling of ancient GRB duplicates has contributed to the creation of large gene deserts associated with regulatory genes in vertebrates, providing a potentially widespread mechanism for the origin of these enigmatic genomic traits.

yap1b, a divergent Yap/Taz family member, cooperates with yap1 in survival and morphogenesis via common transcriptional targets.

Yap1/Taz are well-known Hippo effectors triggering complex transcriptional programs controlling growth, survival and cancer progression. Here, we describe yap1b, a new Yap1/Taz family member with a unique transcriptional activation domain that cannot be phosphorylated by Src/Yes kinases. We show that yap1b evolved specifically in euteleosts (i.e. including medaka but not zebrafish) by duplication and adaptation of yap1. Using DamID-seq, we generated maps of chromatin occupancy for Yap1, Taz (Wwtr1) and Yap1b in gastrulating zebrafish and medaka embryos. Our comparative analyses uncover the genetic programs controlled by Yap family proteins during early embryogenesis, and show largely overlapping targets for Yap1 and Yap1b. CRISPR/Cas9-induced mutation of yap1b in medaka does not result in an overt phenotype during embryogenesis or adulthood. However, yap1b mutation strongly enhances the embryonic malformations observed in yap1 mutants. Thus yap1-/-; yap1b-/- double mutants display more severe body flattening, eye misshaping and increased apoptosis than yap1-/- single mutants, thus revealing overlapping gene functions. Our results indicate that, despite its divergent transactivation domain, Yap1b cooperates with Yap1 to regulate cell survival and tissue morphogenesis during early development.

4Cin: A computational pipeline for 3D genome modeling and virtual Hi-C analyses from 4C data.

The use of 3C-based methods has revealed the importance of the 3D organization of the chromatin for key aspects of genome biology. However, the different caveats of the variants of 3C techniques have limited their scope and the range of scientific fields that could benefit from these approaches. To address these limitations, we present 4Cin, a method to generate 3D models and derive virtual Hi-C (vHi-C) heat maps of genomic loci based on 4C-seq or any kind of 4C-seq-like data, such as those derived from NG Capture-C. 3D genome organization is determined by integrative consideration of the spatial distances derived from as few as four 4C-seq experiments. The 3D models obtained from 4C-seq data, together with their associated vHi-C maps, allow the inference of all chromosomal contacts within a given genomic region, facilitating the identification of Topological Associating Domains (TAD) boundaries. Thus, 4Cin offers a much cheaper, accessible and versatile alternative to other available techniques while providing a comprehensive 3D topological profiling. By studying TAD modifications in genomic structural variants associated to disease phenotypes and performing cross-species evolutionary comparisons of 3D chromatin structures in a quantitative manner, we demonstrate the broad potential and novel range of applications of our method.

Meis1 coordinates a network of genes implicated in eye development and microphthalmia.

Microphthalmos is a rare congenital anomaly characterized by reduced eye size and visual deficits of variable degree. Sporadic and hereditary microphthalmos have been associated with heterozygous mutations in genes fundamental for eye development. Yet, many cases are idiopathic or await the identification of molecular causes. Here we show that haploinsufficiency of Meis1, which encodes a transcription factor with evolutionarily conserved expression in the embryonic trunk, brain and sensory organs, including the eye, causes microphthalmic traits and visual impairment in adult mice. By combining analysis of Meis1 loss-of-function and conditional Meis1 functional rescue with ChIP-seq and RNA-seq approaches we show that, in contrast to its preferential association with Hox-Pbx BSs in the trunk, Meis1 binds to Hox/Pbx-independent sites during optic cup development. In the eye primordium, Meis1 coordinates, in a dose-dependent manner, retinal proliferation and differentiation by regulating genes responsible for human microphthalmia and components of the Notch signaling pathway. In addition, Meis1 is required for eye patterning by controlling a set of eye territory-specific transcription factors, so that in Meis1(-/-) embryos boundaries among the different eye territories are shifted or blurred. We propose that Meis1 is at the core of a genetic network implicated in eye patterning/microphthalmia, and represents an additional candidate for syndromic cases of these ocular malformations.

A novel chromatin insulator regulates the chicken folate receptor gene from the influence of nearby constitutive heterochromatin and the ß-globin locus.

The three-dimensional architecture of genomes provides new insights about genome organization and function, but many aspects remain unsolved at the local genomic scale. Here we investigate the regulation of two erythroid-specific loci, a folate receptor gene (FOLR1) and the ß-globin gene cluster, which are separated by 16kb of constitutive heterochromatin. We found that in early erythroid differentiation the FOLR1 gene presents a permissive chromatin configuration that allows its expression. Once the transition to the next differentiation state occurs, the heterochromatin spreads into the FOLR1 domain, concomitant with the dissociation of CTCF from a novel binding site, thereby resulting in irreversible silencing of the FOLR1 gene. We demonstrate that the sequences surrounding the CTCF-binding site possess classical insulator properties in vitro and in vivo. In contrast, the chicken cHS4 ß-globin insulator present on the other side of the heterochromatic segment is in a constitutive open chromatin configuration, with CTCF constantly bound from the early stages of erythroid differentiation. Therefore, this study demonstrates that the 16kb of constitutive heterochromatin contributes to silencing of the FOLR1 gene during erythroid differentiation.

TEAD and YAP regulate the enhancer network of human embryonic pancreatic progenitors.

The genomic regulatory programmes that underlie human organogenesis are poorly understood. Pancreas development, in particular, has pivotal implications for pancreatic regeneration, cancer and diabetes. We have now characterized the regulatory landscape of embryonic multipotent progenitor cells that give rise to all pancreatic epithelial lineages. Using human embryonic pancreas and embryonic-stem-cell-derived progenitors we identify stage-specific transcripts and associated enhancers, many of which are co-occupied by transcription factors that are essential for pancreas development. We further show that TEAD1, a Hippo signalling effector, is an integral component of the transcription factor combinatorial code of pancreatic progenitor enhancers. TEAD and its coactivator YAP activate key pancreatic signalling mediators and transcription factors, and regulate the expansion of pancreatic progenitors. This work therefore uncovers a central role for TEAD and YAP as signal-responsive regulators of multipotent pancreatic progenitors, and provides a resource for the study of embryonic development of the human pancreas.

Matricellular protein SPARC/osteonectin expression is regulated by DNA methylation in its core promoter region.

BACKGROUND: SPARC/osteonectin is an evolutionarily conserved matricellular protein that modulates cell-matrix interaction and cell function. In all vertebrates, SPARC is dynamically expressed during embryogenesis. However, the precise function of SPARC and the regulatory elements required for its expression in particular during early embryogenesis are largely unknown. RESULTS: The present study was undertaken to explore the molecular mechanisms that regulate sparc gene expression by in vivo functional characterization of the sparc promoter and identification of possible putative regulatory elements that govern basal promoter activity. We report here transient expression analyses of eGFP expression from transgenic zebrafish containing a Sparc-iTol2-eGFP-BAC and/or 7.25 kb-sparc-Tol2-eGFP constructs. eGFP expression was specifically found in the notochord, otic vesicle, fin fold, intermediate cell mass, and olfactory placode of BAC and Tol2 transposon vectors injected embryos. Deletion analysis revealed that promoter activity resides in the unique 5'-untranslated intronic region. Computer-based analysis revealed a putative CpG island immediately proximal to the translation start site within the intron sequence. Global inhibition of methylation with 5-Aza-2-deoxycytidine promoted sparc expression in association with decreasing CpG methylation. CONCLUSIONS: Taken together, these data identify a contributory role for DNA methylation in regulating sparc expression in zebrafish embryogenesis.

sox21a directs lateral line patterning by modulating FGF signaling.

The development of organs composed by repeated functional units is, in many cases, accomplished by the transition of cells from a progenitor to a differentiation domain, triggering a reiterated developmental program. Yet, how these discrete fields are formed during development is still a largely unresolved question. The posterior lateral line (pLL), a sensory organ present in fish and amphibians, develops from a primordium that migrates along the flanks of the animal periodically depositing neuromasts, the pLL functional units. In zebrafish (Danio rerio), the developmental program of the pLL is triggered by the transit of progenitor cells from a Wnt to a Fgf signaling domain. It has been proposed that these two fields are defined by the antagonistic activity of these two signaling pathways, but how they are formed and maintained is still an open question in the development of the pLL. In this work, we show that sox21a, an HMG -box transcription factor, is expressed within the Fgf domain. We demonstrate that, while the Fgf signaling pathway do not control sox21a, knockdown of sox21a causes impairment of Fgf signaling, expansion of the Wnt signaling domain and disruption of neuromast development. These results suggest that sox21a is a key player in the pLL primordium patterning, fine-tuning the border of the Fgf and Wnt signaling domains.

A polymorphic enhancer near GREM1 influences bowel cancer risk through differential CDX2 and TCF7L2 binding.

A rare germline duplication upstream of the bone morphogenetic protein antagonist GREM1 causes a Mendelian-dominant predisposition to colorectal cancer (CRC). The underlying disease mechanism is strong, ectopic GREM1 overexpression in the intestinal epithelium. Here, we confirm that a common GREM1 polymorphism, rs16969681, is also associated with CRC susceptibility, conferring ∼20% differential risk in the general population. We hypothesized the underlying cause to be moderate differences in GREM1 expression. We showed that rs16969681 lies in a region of active chromatin with allele- and tissue-specific enhancer activity. The CRC high-risk allele was associated with stronger gene expression, and higher Grem1 mRNA levels increased the intestinal tumor burden in Apc(Min) mice. The intestine-specific transcription factor CDX2 and Wnt effector TCF7L2 bound near rs16969681, with significantly higher affinity for the risk allele, and CDX2 overexpression in CDX2/GREM1-negative cells caused re-expression of GREM1. rs16969681 influences CRC risk through effects on Wnt-driven GREM1 expression in colorectal tumors.

Restless legs syndrome-associated intronic common variant in Meis1 alters enhancer function in the developing telencephalon.

Genome-wide association studies (GWAS) identified the MEIS1 locus for Restless Legs Syndrome (RLS), but causal single nucleotide polymorphisms (SNPs) and their functional relevance remain unknown. This locus contains a large number of highly conserved noncoding regions (HCNRs) potentially functioning as cis-regulatory modules. We analyzed these HCNRs for allele-dependent enhancer activity in zebrafish and mice and found that the risk allele of the lead SNP rs12469063 reduces enhancer activity in the Meis1 expression domain of the murine embryonic ganglionic eminences (GE). CREB1 binds this enhancer and rs12469063 affects its binding in vitro. In addition, MEIS1 target genes suggest a role in the specification of neuronal progenitors in the GE, and heterozygous Meis1-deficient mice exhibit hyperactivity, resembling the RLS phenotype. Thus, in vivo and in vitro analysis of a common SNP with small effect size showed allele-dependent function in the prospective basal ganglia representing the first neurodevelopmental region implicated in RLS.

Several cis-regulatory elements control mRNA stability, translation efficiency, and expression pattern of Prrxl1 (paired related homeobox protein-like 1).

The homeodomain transcription factor Prrxl1/DRG11 has emerged as a crucial molecule in the establishment of the pain circuitry, in particular spinal cord targeting of dorsal root ganglia (DRG) axons and differentiation of nociceptive glutamatergic spinal cord neurons. Despite Prrxl1 importance in the establishment of the DRG-spinal nociceptive circuit, the molecular mechanisms that regulate its expression along development remain largely unknown. Here, we show that Prrxl1 transcription is regulated by three alternative promoters (named P1, P2, and P3), which control the expression of three distinct Prrxl1 5'-UTR variants, named 5'-UTR-A, 5'-UTR-B, and 5'-UTR-C. These 5'-UTR sequences confer distinct mRNA stability and translation efficiency to the Prrxl1 transcript. The most conserved promoter (P3) contains a TATA-box and displays in vivo enhancer activity in a pattern that overlaps with the zebrafish Prrxl1 homologue, drgx. Regulatory modules present in this sequence were identified and characterized, including a binding site for Phox2b. Concomitantly, we demonstrate that zebrafish Phox2b is required for the expression of drgx in the facial, glossopharyngeal, and vagal cranial ganglia.

Identification of the first recurrent PAR1 deletion in Léri-Weill dyschondrosteosis and idiopathic short stature reveals the presence of a novel SHOX enhancer.

BACKGROUND: SHOX, located in the pseudoautosomal region 1 (PAR1) of the sexual chromosomes, encodes a transcription factor implicated in human growth. Defects in SHOX or its enhancers have been observed in ∼60% of Leri-Weill dyschondrosteosis (LWD) patients, a skeletal dysplasia characterised by short stature and/or the characteristic Madelung deformity, and in 2-5% of idiopathic short stature (ISS). To identify the molecular defect in the remaining genetically undiagnosed LWD and ISS patients, this study screened previously unanalysed PAR1 regions in 124 LWD and 576 ISS probands. METHODS: PAR1 screening was undertaken by multiplex ligation dependent probe amplification (MLPA). Copy number alterations were subsequently confirmed and delimited by locus-specific custom-designed MLPA, array comparative genomic hybridisation (CGH) and breakpoint junction PCR/sequencing. RESULTS: A recurrent PAR1 deletion downstream of SHOX spanning 47543 bp with identical breakpoints was identified in 19 LWD (15.3%) and 11 ISS (1.9%) probands, from 30 unrelated families. Eight evolutionarily conserved regions (ECRs 1-8) identified within the deleted sequence were evaluated for SHOX regulatory activity by means of chromosome conformation capture (3C) in chicken embryo limbs and luciferase reporter assays in human U2OS osteosarcoma cells. The 3C assay indicated potential SHOX regulatory activity by ECR1, which was subsequently confirmed to act as a SHOX enhancer, operating in an orientation and position independent manner, in human U2OS cells. CONCLUSIONS: This study has identified the first recurrent PAR1 deletion in LWD and ISS, which results in the loss of a previously uncharacterised SHOX enhancer. The loss of this enhancer may decrease SHOX transcription, resulting in LWD or ISS due to SHOX haploinsufficiency.

Identification and analysis of conserved cis-regulatory regions of the MEIS1 gene.

Meis1, a conserved transcription factor of the TALE-homeodomain class, is expressed in a wide variety of tissues during development. Its complex expression pattern is likely to be controlled by an equally complex regulatory landscape. Here we have scanned the Meis1 locus for regulatory elements and found 13 non-coding regions, highly conserved between humans and teleost fishes, that have enhancer activity in stable transgenic zebrafish lines. All these regions are syntenic in most vertebrates. The composite expression of all these enhancer elements recapitulate most of Meis1 expression during early embryogenesis, indicating they comprise a basic set of regulatory elements of the Meis1 gene. Using bioinformatic tools, we identify a number of potential binding sites for transcription factors that are compatible with the regulation of these enhancers. Specifically, HHc2:066650, which is expressed in the developing retina and optic tectum, harbors several predicted Pax6 sites. Biochemical, functional and transgenic assays indicate that pax6 genes directly regulate HHc2:066650 activity.

The Osr1 and Osr2 genes act in the pronephric anlage downstream of retinoic acid signaling and upstream of Wnt2b to maintain pectoral fin development.

Vertebrate odd-skipped related genes (Osr) have an essential function during the formation of the intermediate mesoderm (IM) and the kidney structures derived from it. Here, we show that these genes are also crucial for limb bud formation in the adjacent lateral plate mesoderm (LPM). Reduction of zebrafish Osr function impairs fin development by the failure of tbx5a maintenance in the developing pectoral fin bud. Osr morphant embryos show reduced wnt2b expression, and increasing Wnt signaling in Osr morphant embryos partially rescues tbx5a expression. Thus, Osr genes control limb bud development in a non-cell-autonomous manner, probably through the activation of Wnt2b. Finally, we demonstrate that Osr genes are downstream targets of retinoic acid (RA) signaling. Therefore, Osr genes act as a relay within the genetic cascade of fin bud formation: by controlling the expression of the signaling molecule Wnt2ba in the IM they play an essential function transmitting the RA signaling originated in the somites to the LPM.

Engineered Salmonella allows real-time heterologous gene expression monitoring within infected zebrafish embryos.

Microbial host-pathogen interactions have been traditionally well studied at genetic and physiological levels, but cell-resolution analyses have been particularly scarce. This has been especially remarkable for intracellular parasites for two major reasons: first, the inherent loss of bacteria traceability once infects its hosts; second and more important, the limited availability of genetic tools that allow a tight regulated expression of bacterial virulence genes once inside the host tissues. Here we present novel data supporting the use of zebrafish embryos to monitor Salmonella enterica serovar Thyphimurium infection. Intravenous infection of Salmonella can be easily monitored using in vivo fluorescence that allows the visualization of free-swimming bacteria through the circulatory system. Moreover, we have engineered Salmonella to voluntarily activate heterologous gene expression at any point during infection once inside the zebrafish macrophages using a salicylate-based expression system. This approach allows real-time cell-resolution in vivo monitoring of the infection. All together, this approach paves the road to cell-based resolution experiments that would be harder to mimic in other vertebrate infection models.

A conserved function of the chromatin ATPase Kismet in the regulation of hedgehog expression.

The development of the Drosophila melanogaster wing depends on its subdivision into anterior and posterior compartments, which constitute two independent cell lineages since their origin in the embryonic ectoderm. The anterior-posterior compartment boundary is the place where signaling by the Hedgehog pathway takes place, and this requires pathway activation in anterior cells by ligand expressed exclusively in posterior cells. Several mechanisms ensure the confinement of hedgehog expression to posterior cells, including repression by Cubitus interruptus, the co-repressor Groucho and Master of thick veins. In this work we identified Kismet, a chromodomain-containing protein of the SNF2-like family of ATPases, as a novel component of the hedgehog transcriptional repression mechanism in anterior compartment cells. In kismet mutants, hedgehog is ectopically expressed in a domain of anterior cells close to the anterior-posterior compartment boundary, causing inappropriate activation of the pathway and changes in the development of the central region of the wing. The contribution of Kismet to the silencing of hedgehog expression is limited to anterior cells with low levels of the repressor form of Cubitus interruptus. We also show that knockdown of CHD8, the kismet homolog in Xenopus tropicalis, is also associated with ectopic sonic hedgehog expression and up-regulation of one of its target genes in the eye, Pax2, indicating the evolutionary conservation of Kismet/CHD8 function in negatively controlling hedgehog expression.

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.