Comparative 3D genome analysis between neural retina and retinal pigment epithelium reveals differential cis-regulatory interactions at retinal disease loci

BACKGROUND

Vision depends on the interplay between photoreceptor cells of the neural retina and the underlying retinal pigment epithelium (RPE). Most genes involved in inherited retinal diseases display specific spatiotemporal expression within these interconnected retinal components through the local recruitment of cis-regulatory elements (CREs) in 3D nuclear space.

RESULTS

To understand the role of differential chromatin architecture in establishing tissue-specific expression at inherited retinal disease loci, we mapped genome-wide chromatin interactions using in situ Hi-C and H3K4me3 HiChIP on neural retina and RPE/choroid from human adult donor eyes. We observed chromatin looping between active promoters and 32,425 and 8060 candidate CREs in the neural retina and RPE/choroid, respectively. A comparative 3D genome analysis between these two retinal tissues revealed that 56% of 290 known inherited retinal disease genes were marked by differential chromatin interactions. One of these was ABCA4, which is implicated in the most common autosomal recessive inherited retinal disease. We zoomed in on retina- and RPE-specific cis-regulatory interactions at the ABCA4 locus using high-resolution UMI-4C. Integration with bulk and single-cell epigenomic datasets and in vivo enhancer assays in zebrafish revealed tissue-specific CREs interacting with ABCA4.

CONCLUSIONS

Through comparative 3D genome mapping, based on genome-wide, promoter-centric, and locus-specific assays of human neural retina and RPE, we have shown that gene regulation at key inherited retinal disease loci is likely mediated by tissue-specific chromatin interactions. These findings do not only provide insight into tissue-specific regulatory landscapes at retinal disease loci, but also delineate the search space for non-coding genomic variation underlying unsolved inherited retinal diseases.

The little skate genome and the evolutionary emergence of wing-like fins

Skates are cartilaginous fish whose body plan features enlarged wing-like pectoral fins, enabling them to thrive in benthic environments1,2. However, the molecular underpinnings of this unique trait remain unclear. Here we investigate the origin of this phenotypic innovation by developing the little skate Leucoraja erinacea as a genomically enabled model. Analysis of a high-quality chromosome-scale genome sequence for the little skate shows that it preserves many ancestral jawed vertebrate features compared with other sequenced genomes, including numerous ancient microchromosomes. Combining genome comparisons with extensive regulatory datasets in developing fins-including gene expression, chromatin occupancy and three-dimensional conformation-we find skate-specific genomic rearrangements that alter the three-dimensional regulatory landscape of genes that are involved in the planar cell polarity pathway. Functional inhibition of planar cell polarity signalling resulted in a reduction in anterior fin size, confirming that this pathway is a major contributor to batoid fin morphology. We also identified a fin-specific enhancer that interacts with several hoxa genes, consistent with the redeployment of hox gene expression in anterior pectoral fins, and confirmed its potential to activate transcription in the anterior fin using zebrafish reporter assays. Our findings underscore the central role of genome reorganization and regulatory variation in the evolution of phenotypes, shedding light on the molecular origin of an enigmatic trait.

Multiomic atlas with functional stratification and developmental dynamics of zebrafish cis-regulatory elements

Zebrafish, a popular organism for studying embryonic development and for modeling human diseases, has so far lacked a systematic functional annotation program akin to those in other animal models. To address this, we formed the international DANIO-CODE consortium and created a central repository to store and process zebrafish developmental functional genomic data. Our data coordination center ( https://danio-code.zfin.org ) combines a total of 1,802 sets of unpublished and re-analyzed published genomic data, which we used to improve existing annotations and show its utility in experimental design. We identified over 140,000 cis-regulatory elements throughout development, including classes with distinct features dependent on their activity in time and space. We delineated the distinct distance topology and chromatin features between regulatory elements active during zygotic genome activation and those active during organogenesis. Finally, we matched regulatory elements and epigenomic landscapes between zebrafish and mouse and predicted functional relationships between them beyond sequence similarity, thus extending the utility of zebrafish developmental genomics to mammals.

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.

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.

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.

Reprogramming Nuclear Architecture: Just a TAD

In this issue of Cell Stem Cell, Laugsch et al. (2019) use direct reprogramming, epigenetics, and chromatin architecture studies to demonstrate that developmental defects observed in a BOFS patient are caused by reduced expression of TFAP2A in neural crest cells due to the spatial separation of the promoter from its neural crest enhancers.

Assaying Chromatin Accessibility Using ATAC-Seq in Invertebrate Chordate Embryos

Cis-regulatory elements (CREs) are non-coding DNA regions involved in the spatio-temporal regulation of gene expression. Gene regulatory changes drive animal development and play major roles during evolution of animal body plans. Therefore, we believe that determining CREs at different developmental stages and across animal lineages is critical to understand how evolution operates through development. The Assay for Transposase-Accessible Chromatin followed by high-throughput sequencing (ATAC-seq) is a powerful technique for the study of CREs that takes advantage of Tn5 transposase activity. Starting from fewer than 10⁵ cells, in a 1-day procedure, it is possible to detect, at a genome-wide level, CREs located in open chromatin regions with high resolution. Here, we describe a detailed step-by-step ATAC-seq protocol for invertebrate chordate marine embryos. We have successfully applied this technique to amphioxus and two species of tunicate embryos. We also show an easy workflow to analyze data generated with this technique. Moreover, we point out that this method and our bioinformatic pipeline are efficient to detect CREs associated with Wnt signaling pathway by simply using embryos treated with a drug that perturbs this pathway. This approach can be extended to other signaling pathways and also to embryo mutants for critical genes. Our results therefore demonstrate the power of ATAC-seq for the identification of CREs that play essential functions during animal development in a wide range of invertebrate or vertebrate animals.

Characterization of the accessible genome in the human malaria parasite Plasmodium falciparum

Human malaria is a devastating disease and a major cause of poverty in resource-limited countries. To develop and adapt within hosts Plasmodium falciparum undergoes drastic switches in gene expression. To identify regulatory regions in the parasite genome, we performed genome-wide profiling of chromatin accessibility in two culture-adapted isogenic subclones at four developmental stages during the intraerythrocytic cycle by using the Assay for Transposase-Accessible Chromatin by sequencing (ATAC-seq). Tn5 transposase hypersensitivity sites (THSSs) localize preferentially at transcriptional start sites (TSSs). Chromatin accessibility by ATAC-seq is predictive of active transcription and of the levels of histone marks H3K9ac and H3K4me3. Our assay allows the identification of novel regulatory regions including TSS and enhancer-like elements. We show that the dynamics in the accessible chromatin profile matches temporal transcription during development. Motif analysis of stage-specific ATAC-seq sites predicts the in vivo binding sites and function of multiple ApiAP2 transcription factors. At last, the alternative expression states of some clonally variant genes (CVGs), including eba, phist, var and clag genes, associate with a differential ATAC-seq signal at their promoters. Altogether, this study identifies genome-wide regulatory regions likely to play an essential function in the developmental transitions and in CVG expression in P. falciparum.

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.

An evolutionary perspective of regulatory landscape dynamics in development and disease

The organization of animal genomes into topologically associating domains (TADs) provides a structural scaffold in which cis-regulatory elements (CREs) operate on their target genes. Determining the position of CREs and genes relative to TADs has become instrumental to trace gene expression changes during evolution and in diseases. Here we will review recent studies and discuss TADs as structural units with respect to their conservation and stability during genome reorganization. Furthermore, we describe how TAD restructuring contributed to morphological novelties during evolution but also their deleterious effects associated with disease. Despite considering TADs as structural units, the nested and dynamic scaffold within TADs contributes to tissue-specific gene expression, implying that such changes can also account for gene expression differences during evolution.

Cis-regulatory landscapes in development and evolution

The recent advances in our understanding of the 3D organization of the chromatin together with an almost unlimited ability to detect cis-regulatory elements genome-wide using different biochemical signatures has provided us with an unprecedented power to study gene regulation. It is now possible to profile the complete regulatory apparatus controlling the spatio-temporal expression of any given gene, the so-called gene Regulatory Landscapes (RLs). Here we review several studies over the last two years demonstrating the functional consequences of altering RL structure in development, disease and evolution. These works clearly show that a deep understanding of transcriptional regulation is no longer conceivable without considering the 3D modular organization of animal genomes.

The Dynamic Regulatory Genome of Capsaspora and the Origin of Animal Multicellularity

The unicellular ancestor of animals had a complex repertoire of genes linked to multicellular processes. This suggests that changes in the regulatory genome, rather than in gene innovation, were key to the origin of animals. Here, we carry out multiple functional genomic assays in Capsaspora owczarzaki, the unicellular relative of animals with the largest known gene repertoire for transcriptional regulation. We show that changing chromatin states, differential lincRNA expression, and dynamic cis-regulatory sites are associated with life cycle transitions in Capsaspora. Moreover, we demonstrate conservation of animal developmental transcription-factor networks and extensive network interconnection in this premetazoan organism. In contrast, however, Capsaspora lacks animal promoter types, and its regulatory sites are small, proximal, and lack signatures of animal enhancers. Overall, our results indicate that the emergence of animal multicellularity was linked to a major shift in genome cis-regulatory complexity, most notably the appearance of distal enhancer regulation.

A single three-dimensional chromatin compartment in amphioxus indicates a stepwise evolution of vertebrate Hox bimodal regulation

The HoxA and HoxD gene clusters of jawed vertebrates are organized into bipartite three-dimensional chromatin structures that separate long-range regulatory inputs coming from the anterior and posterior Hox-neighboring regions. This architecture is instrumental in allowing vertebrate Hox genes to pattern disparate parts of the body, including limbs. Almost nothing is known about how these three-dimensional topologies originated. Here we perform extensive 4C-seq profiling of the Hox cluster in embryos of amphioxus, an invertebrate chordate. We find that, in contrast to the architecture in vertebrates, the amphioxus Hox cluster is organized into a single chromatin interaction domain that includes long-range contacts mostly from the anterior side, bringing distant cis-regulatory elements into contact with Hox genes. We infer that the vertebrate Hox bipartite regulatory system is an evolutionary novelty generated by combining ancient long-range regulatory contacts from DNA in the anterior Hox neighborhood with new regulatory inputs from the posterior side.

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.

Deletion 1q43-44 in a patient with clinical diagnosis of Warburg-Micro syndrome

Warburg-Micro syndrome (WARBM) is an autosomal recessive syndrome characterized by microcephaly, microphthalmia, microcornea, congenital cataracts, optic atrophy and central nervous system malformations. This syndrome is caused by mutations in the RAB3GAP1/2 and RAB18 genes, part of the Rab family, and in the TBC1D20 gene, which contributes to lipid droplet formation/metabolism. Here we present a patient with clinical diagnosis of WARBM syndrome, who did not have mutations in either the RAB3GAP1/2 genes, in the main exons of RAB18, nor in the TBC1D20 gene. However, the analysis with CGH-array detected a 9.6 Mb deletion at 1q43-qter. We performed a genotype-phenotype correlation using 20 previously published patients in whom the coordinates of the deleted regions were defined. The comparative analysis revealed that the current patient and three of the other 20 patients share the loss of six genes, four of which are related with the family of G proteins, and are strongly expressed in the brain, retina, heart and kidney. Consequently, their haploinsufficiency may result in different combinations of clinical alterations, including some of those of WARBM syndrome. In addition, the haploinsufficiency of other genes may contribute to other defects and clinical variability. Additionally, for the genotype-phenotype correlation, one must also consider molecular pathways that can result in the observed alterations. To early confirm a genetic diagnosis is essential for the patient and family. The current patient was considered as having a recessive syndrome, but since he had a "de novo" deletion, there was not an increased recurrence risk.

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.

The conserved barH-like homeobox-2 gene barhl2 acts downstream of orthodentricle-2 and together with iroquois-3 in establishment of the caudal forebrain signaling center induced by Sonic Hedgehog

In this study, we investigated the gene regulatory network that governs formation of the Zona limitans intrathalamica (ZLI), a signaling center that secretes Sonic Hedgehog (Shh) to control the growth and regionalization of the caudal forebrain. Using loss- and gain-of-function, explants and grafting experiments in amphibians, we demonstrate that barhl2 acts downstream of otx2 and together with the iroquois (irx)-3 gene in establishment of the ZLI compartment initiated by Shh influence. We find that the presumptive (pre)-ZLI domain expresses barhl2, otx2 and irx3, whereas the thalamus territory caudally bordering the pre-ZLI expresses barhl2, otx2 and irx1/2 and early on irx3. We demonstrate that Barhl2 activity is required for determination of the ZLI and thalamus fates and that within the p2 alar plate the ratio of Irx3 to Irx1/2 contributes to ZLI specification and size determination. We show that when continuously exposed to Shh, neuroepithelial cells coexpressing barhl2, otx2 and irx3 acquire two characteristics of the ZLI compartment-the competence to express shh and the ability to segregate from anterior neural plate cells. In contrast, neuroepithelial cells expressing barhl2, otx2 and irx1/2, are not competent to express shh. Noteworthy in explants, under Shh influence, ZLI-like cells segregate from thalamic-like cells. Our study establishes that Barhl2 activity plays a key role in p2 alar plate patterning, specifically ZLI formation, and provides new insights on establishment of the signaling center of the caudal forebrain.

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 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.

Deep conservation of cis-regulatory elements in metazoans

Despite the vast morphological variation observed across phyla, animals share multiple basic developmental processes orchestrated by a common ancestral gene toolkit. These genes interact with each other building complex gene regulatory networks (GRNs), which are encoded in the genome by cis-regulatory elements (CREs) that serve as computational units of the network. Although GRN subcircuits involved in ancient developmental processes are expected to be at least partially conserved, identification of CREs that are conserved across phyla has remained elusive. Here, we review recent studies that revealed such deeply conserved CREs do exist, discuss the difficulties associated with their identification and describe new approaches that will facilitate this search.

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.

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.

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.

Differences in enhancer activity in mouse and zebrafish reporter assays are often associated with changes in gene expression

Background

Phenotypic evolution in animals is thought to be driven in large part by differences in gene expression patterns, which can result from sequence changes in cis-regulatory elements (cis-changes) or from changes in the expression pattern or function of transcription factors (trans-changes). While isolated examples of trans-changes have been identified, the scale of their overall contribution to regulatory and phenotypic evolution remains unclear.

Results

Here, we attempt to examine the prevalence of trans-effects and their potential impact on gene expression patterns in vertebrate evolution by comparing the function of identical human tissue-specific enhancer sequences in two highly divergent vertebrate model systems, mouse and zebrafish. Among 47 human conserved non-coding elements (CNEs) tested in transgenic mouse embryos and in stable zebrafish lines, at least one species-specific expression domain was observed in the majority (83%) of cases, and 36% presented dramatically different expression patterns between the two species. Although some of these discrepancies may be due to the use of different transgenesis systems in mouse and zebrafish, in some instances we found an association between differences in enhancer activity and changes in the endogenous gene expression patterns between mouse and zebrafish, suggesting a potential role for trans-changes in the evolution of gene expression.

Conclusions

In total, our results: (i) serve as a cautionary tale for studies investigating the role of human enhancers in different model organisms, and (ii) suggest that changes in the trans environment may play a significant role in the evolution of gene expression in vertebrates.

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.

Suppression of Bmp4 signaling by the zinc-finger repressors Osr1 and Osr2 is required for Wnt/β-catenin-mediated lung specification in Xenopus

Embryonic development of the respiratory system is regulated by a series of mesenchymal-epithelial interactions that are only partially understood. Mesenchymal FGF and Wnt2/Wnt2b signaling are implicated in specification of mammalian pulmonary progenitors from the ventral foregut endoderm, but their epistatic relationship and downstream targets are largely unknown. In addition, how wnt2 and wnt2b are regulated in the developing foregut mesenchyme is unknown. We show that the Odd-skipped-related (Osr) zinc-finger transcriptional repressors Osr1 and Osr2 are redundantly required for Xenopus lung specification in a molecular pathway linking foregut pattering by FGFs to Wnt-mediated lung specification and RA-regulated lung bud growth. FGF and RA signals are required for robust osr1 and osr2 expression in the foregut endoderm and surrounding lateral plate mesoderm (lpm) prior to respiratory specification. Depletion of both Osr1 and Osr2 (Osr1/Osr2) results in agenesis of the lungs, trachea and esophagus. The foregut lpm of Osr1/Osr2-depleted embryos fails to express wnt2, wnt2b and raldh2, and consequently Nkx2.1(+) progenitors are not specified. Our data suggest that Osr1/Osr2 normally repress bmp4 expression in the lpm, and that BMP signaling negatively regulates the wnt2b domain. These results significantly advance our understanding of early lung development and may impact strategies to differentiate respiratory tissue from stem cells.

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.

Making reporter gene constructs to analyze cis-regulatory elements

Cis-regulatory sequences control when, where, and how much genes are transcribed. A better understanding on these elements is a fundamental keystone to better understand development, cell differentiation, and morphogenesis. Several methods based on in silico analysis or ChIP-seq experiments have been developed to detect cis-acting sequences. Here, we describe a protocol to isolate such sequences from genomic DNA and to clone them into expression vectors for functional assays using the Gateway cloning technology.

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.

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.

An evolutionarily conserved three-dimensional structure in the vertebrate Irx clusters facilitates enhancer sharing and coregulation

Developmental gene clusters are paradigms for the study of gene regulation; however, the mechanisms that mediate phenomena such as coregulation and enhancer sharing remain largely elusive. Here we address this issue by analysing the vertebrate Irx clusters. We first present a deep enhancer screen of a 2-Mbp span covering the IrxA cluster. Using chromosome conformation capture, we show that enhancer sharing is widespread within the cluster, explaining its evolutionarily conserved organization. We also identify a three-dimensional architecture, probably formed through interactions with CCCTC-binding factor, which is present within both Irx clusters of mouse, Xenopus and zebrafish. This architecture brings the promoters of the first two genes together in the same chromatin landscape. We propose that this unique and evolutionarily conserved genomic architecture of the vertebrate Irx clusters is essential for the coregulation of the first two genes and simultaneously maintains the third gene in a partially independent regulatory landscape.

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.

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.

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.

Common inherited variation on human chromosome 8q23 influences the risk of developing colorectal cancer (CRC). To understand the basis of this association we have compared the frequency of common genetic variants at 8q23 in ∼2,000 CRC cases and ∼2,000 healthy controls. Functional analyses of variants strongly associated with CRC risk showed that the single nucleotide polymorphism rs16888589 underscores the 8q23.3 association. The region of the genome harboring rs16888589 increases the expression of the gene for eukaryotic translation initiation factor 3, subunit H. We show that increased expression of this gene increases CRC growth thereby providing a biological mechanism for the 8q23.3 association. This finding is of particular importance in elucidating the etiological basis of CRC.

Conserved developmental expression of Fezf in chordates and Drosophila and the origin of the Zona Limitans Intrathalamica (ZLI) brain organizer

Background

The zona limitans intrathalamica (ZLI) and the isthmus organizer (IsO) are two major secondary organizers of vertebrate brain development. These organizers are located at the interface of the expression domains of key patterning genes (Fezf-Irx and Otx-Gbx, respectively). To gain insights into the evolutionary origin of the ZLI, we studied Fezf in bilaterians.

Results

In this paper, we identified a conserved sequence motif (Fezf box) in all bilaterians. We report the expression pattern of Fezf in amphioxus and Drosophila and compare it with those of Gbx, Otx and Irx. We found that the relative expression patterns of these genes in vertebrates are fully conserved in amphioxus and flies, indicating that the genetic subdivisions defining the location of both secondary organizers in early vertebrate brain development were probably present in the last common ancestor of extant bilaterians. However, in contrast to vertebrates, we found that Irx-defective flies do not show an affected Fezf expression pattern.

Conclusions

The absence of expression of the corresponding morphogens from cells at these conserved genetic boundaries in invertebrates suggests that the organizing properties might have evolved specifically in the vertebrate lineage by the recruitment of key morphogens to these conserved genetic locations.

Genome-wide profiling of p63 DNA-binding sites identifies an element that regulates gene expression during limb development in the 7q21 SHFM1 locus

Heterozygous mutations in p63 are associated with split hand/foot malformations (SHFM), orofacial clefting, and ectodermal abnormalities. Elucidation of the p63 gene network that includes target genes and regulatory elements may reveal new genes for other malformation disorders. We performed genome-wide DNA–binding profiling by chromatin immunoprecipitation (ChIP), followed by deep sequencing (ChIP–seq) in primary human keratinocytes, and identified potential target genes and regulatory elements controlled by p63. We show that p63 binds to an enhancer element in the SHFM1 locus on chromosome 7q and that this element controls expression of DLX6 and possibly DLX5, both of which are important for limb development. A unique micro-deletion including this enhancer element, but not the DLX5/DLX6 genes, was identified in a patient with SHFM. Our study strongly indicates disruption of a non-coding cis-regulatory element located more than 250 kb from the DLX5/DLX6 genes as a novel disease mechanism in SHFM1. These data provide a proof-of-concept that the catalogue of p63 binding sites identified in this study may be of relevance to the studies of SHFM and other congenital malformations that resemble the p63-associated phenotypes.

Mammalian embryonic development requires precise control of gene expression in the right place at the right time. One level of control of gene expression is through cis-regulatory elements controlled by transcription factors. Deregulation of gene expression by mutations in such cis-regulatory elements has been described in developmental disorders. Heterozygous mutations in the transcription factor p63 are found in patients with limb malformations, cleft lip/palate, and defects in skin and other epidermal appendages, through disruption of normal ectodermal development during embryogenesis. We reasoned that the identification of target genes and cis-regulatory elements controlled by p63 would provide candidate genes for defects arising from abnormally regulated ectodermal development. To test our hypothesis, we carried out a genome-wide binding site analysis and identified a large number of target genes and regulatory elements regulated by p63. We further showed that one of these regulatory elements controls expression of DLX6 and possibly DLX5 in the apical ectodermal ridge in the developing limbs. Loss of this element through a micro-deletion was associated with split hand foot malformation (SHFM1). The list of p63 binding sites provides a resource for the identification of mutations that cause ectodermal dysplasias and malformations in humans.

Understanding the regulatory genome

The sequencing of the whole genome of multiple species provides us with the instruction book of how to build an organism and make it work, plus a detailed history of how diversity was generated during evolution. Unfortunately, we still understand only a small fraction, which is locating where genes are and deciphering the proteins they code for. The next step is to understand how the correct amount of gene products are produced in space and time to obtain a fully functioning organism, from the egg to the adult. This is what is known as the regulatory genome, a term coined by Eric H. Davidson. In this review, we examine what we know about gene regulation from a genomic point of view, revise the current in silico, in vitro and in vivo methodological approaches to study transcriptional regulation, and point to the power of phylogenetic footprinting as a guide to regulatory element discovery. The advantages and limitations of each approach are considered, with the emerging view that only large-scale studies and data-crunching will give us insight into the language of genomic regulatory systems, and allow the discovery of regulatory codes in the genome.

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.

The Xenopus Irx genes are essential for neural patterning and define the border between prethalamus and thalamus through mutual antagonism with the anterior repressors Fezf and Arx

The Iroquois (Irx) genes encode homeoproteins conserved during evolution. Vertebrate genomes contain six Irx genes organized in two clusters, IrxA (which harbors Irx1, Irx2 and Irx4) and IrxB (which harbors Irx3, Irx5 and Irx6). To determine the precise role of these genes during development and their putative redundancies, we conducted a comparative expression analysis and a comprehensive loss-of-function study of all the early expressed Irx genes (Irx1-5) using specific morpholinos in Xenopus. We found that the five Irx genes display largely overlapping expression patterns and contribute to neural patterning. All Irx genes are required for proper formation of posterior forebrain, midbrain, hindbrain and, to a lesser an extent, spinal cord. Nevertheless, Irx1 and Irx3 seem to have a predominant role during regionalization of the neural plate. In addition, we find that the common anterior limit of Irx gene expression, which will correspond to the future border between the prethalamus and thalamus, is defined by mutual repression between Fezf and Irx proteins. This mutual repression is likely direct. Finally, we show that Arx, another anteriorly expressed repressor, also contribute to delineate the anterior border of Irx expression.

A dual requirement for Iroquois genes during Xenopus kidney development

The Iroquois (Irx) genes encode evolutionary conserved homeoproteins. We report that Xenopus genes Irx1 and Irx3 are expressed and required during different stages of Xenopus pronephros development. They are initially expressed during mid-neurulation in domains extending over most of the prospective pronephric territory. Expression onset takes place after kidney anlage specification, but before pronephric organogenesis occurs. Later, during nephron segmentation, expression becomes restricted to the intermediate tubule region of the proximal-distal axis. Loss- and gain-of-function analyses, performed with specific morpholinos and inducible wild-type and dominant-negative constructs, reveal a dual requirement for Irx1 and Irx3 during pronephros development. During neurula stages, these genes maintain the specification of the pronephric territory and define its size. This seems to occur, at least in part, through positive regulation of Bmp signalling. Subsequently, Irx genes are required for proper formation of the intermediate tubule. Finally, we find that retinoic acid signalling activates both Irx1 and Irx3 genes in the pronephros.