The cellular landscape of the endochondral bone during the transition to extrauterine life.

The cellular complexity of the endochondral bone underlies its essential and pleiotropic roles during organismal life. While the adult bone has received significant attention, we still lack a deep understanding of the perinatal bone cellulome. Here, we have profiled the full composition of the murine endochondral bone at the single-cell level during the transition from fetal to newborn life and in comparison with the adult tissue, with particular emphasis on the mesenchymal compartment. The perinatal bone contains different fibroblastic clusters with blastema-like characteristics in organizing and supporting skeletogenesis, angiogenesis and hematopoiesis. Our data also suggest dynamic inter- and intra-compartment interactions, as well as a bone marrow milieu that seems prone to anti-inflammation, which we hypothesize is necessary to ensure the proper program of lymphopoiesis and the establishment of central and peripheral tolerance in early life. Our study provides an integrative roadmap for the future design of genetic and cellular functional assays to validate cellular interactions and lineage relationships within the perinatal bone.

Trap-TRAP, a Versatile Tool for Tissue-Specific Translatomics in Zebrafish.

Developmental and physiological processes depend on the transcriptional and translational activity of heterogeneous cell populations. A main challenge in gene expression studies is dealing with this intrinsic complexity while keeping sequencing efficiency. Translating ribosome affinity purification (TRAP) methods have allowed cell-specific recovery of polyribosome-associated RNAs by genetic tagging of ribosomes in selected cell populations. Here we combined the TRAP approach with adapted enhancer trap methods (trap-TRAP) to systematically generate zebrafish transgenic lines suitable for tissue-specific translatome interrogation. Through the random integration of a GFP-tagged version of the large subunit ribosomal protein L10a (EGFP-Rpl10a), we have generated stable lines driving expression in a variety of tissues, including the retina, skeletal muscle, lateral line primordia, rhombomeres, or jaws. To increase the range of applications, a UAS:TRAP transgenic line compatible with available Gal4 lines was also generated and tested. The resulting collection of lines and applications constitutes a resource for the zebrafish community in developmental genetics, organ physiology and disease modelling.

Boundary sequences flanking the mouse tyrosinase locus ensure faithful pattern of gene expression.

Control of gene expression is dictated by cell-type specific regulatory sequences that physically organize the structure of chromatin, including promoters, enhancers and insulators. While promoters and enhancers convey cell-type specific activating signals, insulators prevent the cross-talk of regulatory elements within adjacent loci and safeguard the specificity of action of promoters and enhancers towards their targets in a tissue specific manner. Using the mouse tyrosinase (Tyr) locus as an experimental model, a gene whose mutations are associated with albinism, we described the chromatin structure in cells at two distinct transcriptional states. Guided by chromatin structure, through the use of Chromosome Conformation Capture (3C), we identified sequences at the 5' and 3' boundaries of this mammalian gene that function as enhancers and insulators. By CRISPR/Cas9-mediated chromosomal deletion, we dissected the functions of these two regulatory elements in vivo in the mouse, at the endogenous chromosomal context, and proved their mechanistic role as genomic insulators, shielding the Tyr locus from the expression patterns of adjacent genes.

Profiling of conserved non-coding elements upstream of SHOX and functional characterisation of the SHOX cis-regulatory landscape.

Genetic defects such as copy number variations (CNVs) in non-coding regions containing conserved non-coding elements (CNEs) outside the transcription unit of their target gene, can underlie genetic disease. An example of this is the short stature homeobox (SHOX) gene, regulated by seven CNEs located downstream and upstream of SHOX, with proven enhancer capacity in chicken limbs. CNVs of the downstream CNEs have been reported in many idiopathic short stature (ISS) cases, however, only recently have a few CNVs of the upstream enhancers been identified. Here, we set out to provide insight into: (i) the cis-regulatory role of these upstream CNEs in human cells, (ii) the prevalence of upstream CNVs in ISS, and (iii) the chromatin architecture of the SHOX cis-regulatory landscape in chicken and human cells. Firstly, luciferase assays in human U2OS cells, and 4C-seq both in chicken limb buds and human U2OS cells, demonstrated cis-regulatory enhancer capacities of the upstream CNEs. Secondly, CNVs of these upstream CNEs were found in three of 501 ISS patients. Finally, our 4C-seq interaction map of the SHOX region reveals a cis-regulatory domain spanning more than 1 Mb and harbouring putative new cis-regulatory elements.

MIR retrotransposon sequences provide insulators to the human genome.

Insulators are regulatory elements that help to organize eukaryotic chromatin via enhancer-blocking and chromatin barrier activity. Although there are several examples of transposable element (TE)-derived insulators, the contribution of TEs to human insulators has not been systematically explored. Mammalian-wide interspersed repeats (MIRs) are a conserved family of TEs that have substantial regulatory capacity and share sequence characteristics with tRNA-related insulators. We sought to evaluate whether MIRs can serve as insulators in the human genome. We applied a bioinformatic screen using genome sequence and functional genomic data from CD4(+) T cells to identify a set of 1,178 predicted MIR insulators genome-wide. These predicted MIR insulators were computationally tested to serve as chromatin barriers and regulators of gene expression in CD4(+) T cells. The activity of predicted MIR insulators was experimentally validated using in vitro and in vivo enhancer-blocking assays. MIR insulators are enriched around genes of the T-cell receptor pathway and reside at T-cell-specific boundaries of repressive and active chromatin. A total of 58% of the MIR insulators predicted here show evidence of T-cell-specific chromatin barrier and gene regulatory activity. MIR insulators appear to be CCCTC-binding factor (CTCF) independent and show a distinct local chromatin environment with marked peaks for RNA Pol III and a number of histone modifications, suggesting that MIR insulators recruit transcriptional complexes and chromatin modifying enzymes in situ to help establish chromatin and regulatory domains in the human genome. The provisioning of insulators by MIRs across the human genome suggests a specific mechanism by which TE sequences can be used to modulate gene regulatory networks.

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.

Long-range regulatory interactions at the 4q25 atrial fibrillation risk locus involve PITX2c and ENPEP.

BACKGROUND: Recent genome-wide association studies have uncovered genomic loci that underlie an increased risk for atrial fibrillation, the major cardiac arrhythmia in humans. The most significant locus is located in a gene desert at 4q25, approximately 170 kilobases upstream of PITX2, which codes for a transcription factor involved in embryonic left-right asymmetry and cardiac development. However, how this genomic region functionally and structurally relates to PITX2 and atrial fibrillation is unknown. RESULTS: To characterise its function, we tested genomic fragments from 4q25 for transcriptional activity in a mouse atrial cardiomyocyte cell line and in transgenic mouse embryos, identifying a non-tissue-specific potentiator regulatory element. Chromosome conformation capture revealed that this region physically interacts with the promoter of the cardiac specific isoform of Pitx2. Surprisingly, this regulatory region also interacts with the promoter of the next neighbouring gene, Enpep, which we show to be expressed in regions of the developing mouse heart essential for cardiac electrical activity. CONCLUSIONS: Our data suggest that de-regulation of both PITX2 and ENPEP could contribute to an increased risk of atrial fibrillation in carriers of disease-associated variants, and show the challenges that we face in the functional analysis of genome-wide disease associations.

Comparative epigenomics in distantly related teleost species identifies conserved cis-regulatory nodes active during the vertebrate phylotypic period.

The complex relationship between ontogeny and phylogeny has been the subject of attention and controversy since von Baer's formulations in the 19th century. The classic concept that embryogenesis progresses from clade general features to species-specific characters has often been revisited. It has become accepted that embryos from a clade show maximum morphological similarity at the so-called phylotypic period (i.e., during mid-embryogenesis). According to the hourglass model, body plan conservation would depend on constrained molecular mechanisms operating at this period. More recently, comparative transcriptomic analyses have provided conclusive evidence that such molecular constraints exist. Examining cis-regulatory architecture during the phylotypic period is essential to understand the evolutionary source of body plan stability. Here we compare transcriptomes and key epigenetic marks (H3K4me3 and H3K27ac) from medaka (Oryzias latipes) and zebrafish (Danio rerio), two distantly related teleosts separated by an evolutionary distance of 115-200 Myr. We show that comparison of transcriptome profiles correlates with anatomical similarities and heterochronies observed at the phylotypic stage. Through comparative epigenomics, we uncover a pool of conserved regulatory regions (≈700), which are active during the vertebrate phylotypic period in both species. Moreover, we show that their neighboring genes encode mainly transcription factors with fundamental roles in tissue specification. We postulate that these regulatory regions, active in both teleost genomes, represent key constrained nodes of the gene networks that sustain the vertebrate body plan.

Directional tissue migration through a self-generated chemokine gradient.

The directed migration of cell collectives is a driving force of embryogenesis. The predominant view in the field is that cells in embryos navigate along pre-patterned chemoattractant gradients. One hypothetical way to free migrating collectives from the requirement of long-range gradients would be through the self-generation of local gradients that travel with them, a strategy that potentially allows self-determined directionality. However, a lack of tools for the visualization of endogenous guidance cues has prevented the demonstration of such self-generated gradients in vivo. Here we define the in vivo dynamics of one key guidance molecule, the chemokine Cxcl12a, by applying a fluorescent timer approach to measure ligand-triggered receptor turnover in living animals. Using the zebrafish lateral line primordium as a model, we show that migrating cell collectives can self-generate gradients of chemokine activity across their length via polarized receptor-mediated internalization. Finally, by engineering an external source of the atypical receptor Cxcr7 that moves with the primordium, we show that a self-generated gradient mechanism is sufficient to direct robust collective migration. This study thus provides, to our knowledge, the first in vivo proof for self-directed tissue migration through local shaping of an extracellular cue and provides a framework for investigating self-directed migration in many other contexts including cancer invasion.

The developmental epigenomics toolbox: ChIP-seq and MethylCap-seq profiling of early zebrafish embryos.

Genome-wide profiling of DNA methylation and histone modifications answered many questions as to how the genes are regulated on a global scale and what their epigenetic makeup is. Yet, little is known about the function of these marks during early vertebrate embryogenesis. Here we provide detailed protocols for ChIP-seq and MethylCap-seq procedures applied to zebrafish (Danio rerio) embryonic material at four developmental stages. As a proof of principle, we have profiled on a global scale a number of post-translational histone modifications including H3K4me1, H3K4me3 and H3K27ac. We demonstrate that these marks are dynamic during early development and that such developmental transitions can be detected by ChIP-seq. In addition, we applied MethylCap-seq to show that developmentally-regulated DNA methylation remodeling can be detected by such a procedure. Our MethylCap-seq data concur with previous DNA methylation studies of early zebrafish development rendering this method highly suitable for the global assessment of DNA methylation in early vertebrate embryos.

Numb/Numbl-Opo antagonism controls retinal epithelium morphogenesis by regulating integrin endocytosis.

Polarized trafficking of adhesion receptors plays a pivotal role in controlling cellular behavior during morphogenesis. Particularly, clathrin-dependent endocytosis of integrins has long been acknowledged as essential for cell migration. However, little is known about the contribution of integrin trafficking to epithelial tissue morphogenesis. Here we show how the transmembrane protein Opo, previously described for its essential role during optic cup folding, plays a fundamental role in this process. Through interaction with the PTB domain of the clathrin adaptors Numb and Numbl via an integrin-like NPxF motif, Opo antagonizes Numb/Numbl function and acts as a negative regulator of integrin endocytosis in vivo. Accordingly, numb/numbl gain-of-function experiments in teleost embryos mimic the retinal malformations observed in opo mutants. We propose that developmental regulator Opo enables polarized integrin localization by modulating Numb/Numbl, thus directing the basal constriction that shapes the vertebrate retina epithelium.

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

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.

Dynamics of enhancer chromatin signatures mark the transition from pluripotency to cell specification during embryogenesis.

The generation of distinctive cell types that form different tissues and organs requires precise, temporal and spatial control of gene expression. This depends on specific cis-regulatory elements distributed in the noncoding DNA surrounding their target genes. Studies performed on mammalian embryonic stem cells and Drosophila embryos suggest that active enhancers form part of a defined chromatin landscape marked by histone H3 lysine 4 mono-methylation (H3K4me1) and histone H3 lysine 27 acetylation (H3K27ac). Nevertheless, little is known about the dynamics and the potential roles of these marks during vertebrate embryogenesis. Here, we provide genomic maps of H3K4me1/me3 and H3K27ac at four developmental time-points of zebrafish embryogenesis and analyze embryonic enhancer activity. We find that (1) changes in H3K27ac enrichment at enhancers accompany the shift from pluripotency to tissue-specific gene expression, (2) in early embryos, the peaks of H3K27ac enrichment are bound by pluripotent factors such as Nanog, and (3) the degree of evolutionary conservation is higher for enhancers that become marked by H3K27ac at the end of gastrulation, suggesting their implication in the establishment of the most conserved (phylotypic) transcriptome that is known to occur later at the pharyngula stage.

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.

A role for insulator elements in the regulation of gene expression response to hypoxia.

Hypoxia inducible factor (HIF) up-regulates the transcription of a few hundred genes required for the adaptation to hypoxia. This restricted set of targets is in sharp contrast with the widespread distribution of the HIF binding motif throughout the genome. Here, we investigated the transcriptional response of GYS1 and RUVBL2 genes to hypoxia to understand the mechanisms that restrict HIF activity toward specific genes. GYS1 and RUVBL2 genes are encoded by opposite DNA strands and separated by a short intergenic region (~1 kb) that contains a functional hypoxia response element equidistant to both genes. However, hypoxia induced the expression of GYS1 gene only. Analysis of the transcriptional response of chimeric constructs derived from the intergenic region revealed an inhibitory sequence whose deletion allowed RUVBL2 induction by HIF. Enhancer blocking assays, performed in cell culture and transgenic zebrafish, confirmed the existence of an insulator element within this inhibitory region that could explain the differential regulation of GYS1 and RUVBL2 by hypoxia. Hence, in this model, the selective response to HIF is achieved with the aid of insulator elements. This is the first report suggesting a role for insulators in the regulation of differential gene expression in response to environmental signals.

Genome-wide CTCF distribution in vertebrates defines equivalent sites that aid the identification of disease-associated genes.

Many genomic alterations associated with human diseases localize in noncoding regulatory elements located far from the promoters they regulate, making it challenging to link noncoding mutations or risk-associated variants with target genes. The range of action of a given set of enhancers is thought to be defined by insulator elements bound by the 11 zinc-finger nuclear factor CCCTC-binding protein (CTCF). Here we analyzed the genomic distribution of CTCF in various human, mouse and chicken cell types, demonstrating the existence of evolutionarily conserved CTCF-bound sites beyond mammals. These sites preferentially flank transcription factor-encoding genes, often associated with human diseases, and function as enhancer blockers in vivo, suggesting that they act as evolutionarily invariant gene boundaries. We then applied this concept to predict and functionally demonstrate that the polymorphic variants associated with multiple sclerosis located within the EVI5 gene impinge on the adjacent gene GFI1.

The Ste20 kinase misshapen is essential for the invasive behaviour of ovarian epithelial cells in Drosophila.

Stationary-to-migratory transitions of epithelial cells have a key role in development and tumour progression. Border cell migration is a powerful system in which to investigate this transition in living organisms. Here, we identify the Ste20-like kinase misshapen (msn) as a novel regulator of border-cell migration in Drosophila. Expression of msn in border cells is independent of the transcription factor slow border cells and of inputs from all pathways that are known to control border-cell migration. The msn gene functions to modulate the levels and/or distribution of Drosophila E-cadherin to promote the invasive migratory behaviour of border cells.

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.

Integrins contribute to the establishment and maintenance of cell polarity in the follicular epithelium of the Drosophila ovary.

The generation of epithelial cell polarity is a key process during development. Although the induction and orientation of cell polarity by cell-cell and cell-extracellular matrix (ECM) interactions is well established, the molecular mechanisms by which signals from the ECM control cell polarity in developing epithelial tissues remain poorly understood. Here, we have used the follicular epithelium of the Drosophila ovary to investigate the role that integrins, the main cell-ECM receptors, play in the establishment of apicobasal polarity. Mature follicle cells have an apical side facing the germ line and a basal side in contact with a basement membrane. Our results show that integrins - presumably via interactions with the basement membrane - play a reinforcing role in follicle cell polarization, as they are required to establish and/or maintain follicle cell membrane asymmetry only when contact with the germ line is prevented. We suggest that the primary cue for polarization of the follicular epithelium is contact with the germline cells. In addition, while interfering with apical and lateral polarization cues leads to apoptosis, we show here that inhibition of contact with the basement membrane mediated by integrins does not affect cell survival. Finally, we provide evidence to suggest that integrins are required to orientate epithelial polarity in vivo.

A role for the chaperone Hsp70 in the regulation of border cell migration in the Drosophila ovary.

Unravelling the molecular mechanisms that govern cell migration is of great importance towards understanding both normal embryogenesis and physiological and pathological processes occurring in the adult. Migration of border cells (BCs) during Drosophila oogenesis provides a simple and attractive model in which to address this problem. Here, we show that the molecular chaperone Hsp70 is required for BC migration. Thus, BCs lacking all Hsp70 genes present in the fly genome fail to reorganize their actin cytoskeleton, resulting in migration defects. Similar defects are found when the Hsp70 co-chaperone DnaJ-1, the Drosophila homolog of the human Hsp40, is overexpressed specifically in BCs. In addition, we provide biochemical and genetic evidence for an interaction between DnaJ-1 and PDGF/VEGF receptor (PVR), which is also required for actin-mediated BC migration. Furthermore, our results showing that PVR also interacts genetically with Hsp70 suggest that a mechanism by which the DnaJ-1/Hsp70 chaperone complex regulates BC migration is by modulating PVR function.