Lens development depends on a pair of highly conserved Sox21 regulatory elements

Functional associations of evolutionarily recent human genes exhibit sensitivity to the 3D genome landscape and disease

Genome organization is intricately tied to regulating genes and associated cell fate decisions. Here, we examine the positioning and functional significance of human genes, grouped by their lineage restriction level, within the 3D organization of the genome. We reveal that genes of different lineage restriction levels have distinct positioning relationships with both domains and loop anchors, and remarkably consistent relationships with boundaries across cell types. While the functional associations of each group of genes are primarily cell type-specific, associations of conserved genes maintain greater stability across 3D genomic features and disease than recently evolved genes. Furthermore, the expression of these genes across various tissues follows an evolutionary progression, such that RNA levels increase from young lineage restricted genes to ancient genes present in most species. Thus, the distinct relationships of gene evolutionary age, function, and positioning within 3D genomic features contribute to tissue-specific gene regulation in development and disease.

Congenital microcoria deletion in mouse links Sox21 dysregulation to disease and suggests a role for TGFB2 in glaucoma and myopia

Congenital microcoria (MCOR) is a rare hereditary developmental defect of the iris dilator muscle frequently associated with high axial myopia and high intraocular pressure (IOP) glaucoma. The condition is caused by submicroscopic rearrangements of chromosome 13q32.1. However, the mechanisms underlying the failure of iris development and the origin of associated features remain elusive. Here, we present a 3D architecture model of the 13q32.1 region, demonstrating that MCOR-related deletions consistently disrupt the boundary between two topologically associating domains (TADs). Deleting the critical MCOR-causing region in mice reveals ectopic Sox21 expression precisely aligning with Dct, each located in one of the two neighbor TADs. This observation is consistent with the TADs' boundary alteration and adoption of Dct regulatory elements by the Sox21 promoter. Additionally, we identify Tgfb2 as a target gene of SOX21 and show TGFΒ2 accumulation in the aqueous humor of an MCOR-affected subject. Accumulation of TGFB2 is recognized for its role in glaucoma and potential impact on axial myopia. Our results highlight the importance of SOX21-TGFB2 signaling in iris development and control of eye growth and IOP. Insights from MCOR studies may provide therapeutic avenues for this condition but also for glaucoma and high myopia conditions, affecting millions of people.

Sox21 homeologs autoregulate expression levels to control progression through neurogenesis

The SRY HMG box transcription factor Sox21 plays multiple critical roles in neurogenesis, with its function dependent on concentration and developmental stage. In the allotetraploid Xenopus laevis, there are two homeologs of sox21, namely sox21.S and sox21.L. Previous studies focused on Sox21.S, but its amino acid sequence is divergent, lacking conserved poly-A stretches and bearing more similarity with ancestral homologs. In contrast, Sox21.L shares higher sequence similarity with mouse and chick Sox21. To determine if Sox21.S and Sox21.L have distinct functions, we conducted gain and loss-of-function studies in Xenopus embryos. Our studies revealed that Sox21.S and Sox21.L are functionally redundant, but Sox21.L is more effective at driving changes than Sox21.S. These results also support our earlier findings in ectodermal explants, demonstrating that Sox21 function is dose-dependent. While Sox21 is necessary for primary neuron formation, high levels prevent their formation. Strikingly, these proteins autoregulate, with high levels of Sox21.L reducing sox21.S and sox21.L mRNA levels, and decreased Sox21.S promoting increased expression of sox21.L. Our findings shed light on the intricate concentration-dependent roles of Sox21 homeologs in Xenopus neurogenesis.

Zfp503/Nlz2 Is Required for RPE Differentiation and Optic Fissure Closure

Purpose

Uveal coloboma is a congenital eye malformation caused by failure of the optic fissure to close in early human development. Despite significant progress in identifying genes whose regulation is important for executing this closure, mutations are detected in a minority of cases using known gene panels, implying additional genetic complexity. We have previously shown knockdown of znf503 (the ortholog of mouse Zfp503) in zebrafish causes coloboma. Here we characterize Zfp503 knockout (KO) mice and evaluate transcriptomic profiling of mutant versus wild-type (WT) retinal pigment epithelium (RPE)/choroid.

Methods

Zfp503 KO mice were generated by gene targeting using homologous recombination. Embryos were characterized grossly and histologically. Patterns and level of developmentally relevant proteins/genes were examined with immunostaining/in situ hybridization. The transcriptomic profile of E11.5 KO RPE/choroid was compared to that of WT.

Results

Zfp503 is dynamically expressed in developing mouse eyes, and loss of its expression results in uveal coloboma. KO embryos exhibit altered mRNA levels and expression patterns of several key transcription factors involved in eye development, including Otx2, Mitf, Pax6, Pax2, Vax1, and Vax2, resulting in a failure to maintain the presumptive RPE, as evidenced by reduced melanin pigmentation and its differentiation into a neural retina-like lineage. Comparison of RNA sequencing data from WT and KO E11.5 embryos demonstrated reduced expression of melanin-related genes and significant overlap with genes known to be dynamically regulated at the optic fissure.

Conclusions

These results demonstrate a critical role of Zfp503 in maintaining RPE fate and optic fissure closure.

Genetic disruption of zebrafish mab21l1 reveals a conserved role in eye development and affected pathways

BACKGROUND

The male-abnormal 21 like (MAB21L) genes are important in human ocular development. Homozygous loss of MAB21L1 leads to corneal dystrophy in all affected individuals along with cataracts and buphthalmos in some. The molecular function and downstream pathways of MAB21L factors are largely undefined.

RESULTS

We generated the first mab21l1 zebrafish mutant carrying a putative loss-of-function allele, c.107delA p.(Lys36Argfs*7). At the final stages of embryonic development, homozygous mab21l1c.107delA fish displayed enlarged anterior chambers and corneal thinning which progressed with age. Additional studies revealed increased cell death in the mutant corneas, transformation of the cornea into a skin-like epithelium, and progressive lens degeneration with development of fibrous masses in the anterior chamber. RNA-seq of wild-type and mutant ocular transcriptomes revealed significant changes in expression of several genes, including irf1a and b, stat1, elf3, krt17, tlr9, and loxa associated with immunity and/or corneal function. Abnormal expression of lysyl oxidases have been previously linked with corneal thinning, fibrosis, and lens defects in mammals, suggesting a role for loxa misexpression in the progressive mab21l1c.107delA eye phenotype.

CONCLUSIONS

Zebrafish mab21l1 is essential for normal corneal development, similar to human MAB21L1. The identified molecular changes in mab21l1c.107delA mutants provide the first clues about possible affected pathways.

Chromosome-level assembly of southern catfish (silurus meridionalis) provides insights into visual adaptation to nocturnal and benthic lifestyles

The Southern catfish (Silurus meridionalis) is a nocturnal and benthic freshwater fish endemic to the Yangtze River and its tributaries. In this study, we constructed a chromosome-level draft genome of S. meridionalis using 69.7-Gb Nanopore long reads and 49.5-Gb Illumina short reads. The genome assembly was 741.2 Mb in size with a contig N50 of 13.19 Mb. An additional 116.4 Gb of Bionano and 77.4 Gb of Hi-C data were applied to assemble contigs into scaffolds and further into 29 chromosomes, resulting in a 738.9-Mb genome with a scaffold N50 of 28.04 Mb. A total of 22,965 protein-coding genes were predicted from the genome with 22,519 (98.06%) genes functionally annotated. Comparative genomic and transcriptomic analyses revealed a rod-dominated visual system which was responsible for scotopic vision. The absence of cone opsins SWS1 and SWS2 resulted in the lack of ultraviolet and blue violet sensitivity. Mutations at key amino acid sites of RH1.1, RH1.2 and RH2 resulted in spectral tuning good for dim light vision and narrow colour vision. A higher expression level of rod phototransduction genes than that of cone genes and higher rod-to-cone ratio led to higher optical sensitivity under dim light conditions. In addition, analysis of the genes involved in eye morphogenesis and development revealed the loss of some conserved noncoding elements, which might be associated with the small eyes in catfish. Together, our study provides important clues for the adaptation of the catfish visual system to the nocturnal and benthic lifestyles. The draft genome of S. meridionalis represents a valuable resource for studies of the molecular mechanisms of ecological adaptation.

Congenital Microcoria: Clinical Features and Molecular Genetics

Iris integrity is required to regulate both the amount of light reaching the retina and intraocular pressure (IOP), with elevated IOP being a major risk factor for glaucoma. Congenital microcoria (MCOR) is an extremely rare, autosomal dominant disease affecting iris development and hindering both of these functions. It is characterized by absent or underdeveloped dilator muscle fibers and immaturity of the iridocorneal angle—where the aqueous humor is drained—which play a central role in IOP regulation. The dilator muscle anomaly is manifested in pinhole pupils (<2 mm) and thin transilluminable irises, causing both hemeralopia and photoaversion. Axial myopia and juvenile open-angle glaucoma are very frequent (80% and 30% of all cases, respectively). It has been suggested that the immaturity of the chamber angle contributes to glaucoma, and myopia has been ascribed to photoaversion and elevated IOP. Though possible, these mechanisms are insufficient. The disease has been tied to chromosome 13q32.1 structural variations. In addition to compromising iris development, modification of the 13q32.1 architecture could alter signaling pathways for axial ocular length and IOP regulation. Here, we summarize the clinical, histological, and molecular features of this disease, and we discuss the possible etiology of associated anomalies.

Ultraconserved Elements Occupy Specific Arenas of Three-Dimensional Mammalian Genome Organization

This study explores the relationship between three-dimensional genome organization and ultraconserved elements (UCEs), an enigmatic set of DNA elements that are perfectly conserved between the reference genomes of distantly related species. Examining both human and mouse genomes, we interrogate the relationship of UCEs to three features of chromosome organization derived from Hi-C studies. We find that UCEs are enriched within contact domains and, further, that the subset of UCEs within domains shared across diverse cell types are linked to kidney-related and neuronal processes. In boundaries, UCEs are generally depleted, with those that do overlap boundaries being overrepresented in exonic UCEs. Regarding loop anchors, UCEs are neither overrepresented nor underrepresented, but those present in loop anchors are enriched for splice sites. Finally, as the relationships between UCEs and human Hi-C features are conserved in mouse, our findings suggest that UCEs contribute to interspecies conservation of genome organization and, thus, genome stability.

SoxB2 in sea urchin development: implications in neurogenesis, ciliogenesis and skeletal patterning

Background

Current studies in evolutionary developmental biology are focused on the reconstruction of gene regulatory networks in target animal species. From decades, the scientific interest on genetic mechanisms orchestrating embryos development has been increasing in consequence to the fact that common features shared by evolutionarily distant phyla are being clarified. In 2011, a study across eumetazoan species showed for the first time the existence of a highly conserved non-coding element controlling the SoxB2 gene, which is involved in the early specification of the nervous system. This discovery raised several questions about SoxB2 function and regulation in deuterostomes from an evolutionary point of view.

Results

Due to the relevant phylogenetic position within deuterostomes, the sea urchin Strongylocentrotus purpuratus represents an advantageous animal model in the field of evolutionary developmental biology. Herein, we show a comprehensive study of SoxB2 functions in sea urchins, in particular its expression pattern in a wide range of developmental stages, and its co-localization with other neurogenic markers, as SoxB1, SoxC and Elav. Moreover, this work provides a detailed description of the phenotype of sea urchin SoxB2 knocked-down embryos, confirming its key function in neurogenesis and revealing, for the first time, its additional roles in oral and aboral ectoderm cilia and skeletal rod morphology.

Conclusions

We concluded that SoxB2 in sea urchins has a neurogenic function; however, this gene could have multiple roles in sea urchin embryogenesis, expanding its expression in non-neurogenic cells. We showed that SoxB2 is functionally conserved among deuterostomes and suggested that in S. purpuratus this gene acquired additional functions, being involved in ciliogenesis and skeletal patterning.

Electronic supplementary material

The online version of this article (10.1186/s13227-018-0094-1) contains supplementary material, which is available to authorized users.

Genomic features of human limb specific enhancers

To elucidate important cellular and molecular interactions that regulate patterning and skeletal development, vertebrate limbs served as a model organ. A growing body of evidence from detailed studies on a subset of limb regulators like the HOXD cluster or SHH, reveals the importance of enhancers in limb related developmental and disease processes. Exploiting the recent genome-wide availability of functionally confirmed enhancer dataset, this study establishes regulatory interactions for dozens of human limb developmental genes. From these data, it appears that the long-range regulatory interactions are fairly common during limb development. This observation highlights the significance of chromosomal breaks/translocations in human limb deformities. Transcriptional factor (TF) analysis predicts that the differentiation of early nascent limb-bud into future territories entail distinct TF interaction networks. Conclusively, an important motivation for annotating the human limb specific regulatory networks is to pave way for the systematic exploration of their role in disease and evolution.

Evolution of lineage-specific functions in ancient cis-regulatory modules

Morphological evolution is driven both by coding sequence variation and by changes in regulatory sequences. However, how cis-regulatory modules (CRMs) evolve to generate entirely novel expression domains is largely unknown. Here, we reconstruct the evolutionary history of a lens enhancer located within a CRM that not only predates the lens, a vertebrate innovation, but bilaterian animals in general. Alignments of orthologous sequences from different deuterostomes sub-divide the CRM into a deeply conserved core and a more divergent flanking region. We demonstrate that all deuterostome flanking regions, including invertebrate sequences, activate gene expression in the zebrafish lens through the same ancient cluster of activator sites. However, levels of gene expression vary between species due to the presence of repressor motifs in flanking region and core. These repressor motifs are responsible for the relatively weak enhancer activity of tetrapod flanking regions. Ray-finned fish, however, have gained two additional lineage-specific activator motifs which in combination with the ancient cluster of activators and the core constitute a potent lens enhancer. The exploitation and modification of existing regulatory potential in flanking regions but not in the highly conserved core might represent a more general model for the emergence of novel regulatory functions in complex CRMs.

Identification and functional characterization of novel transcriptional enhancers involved in regulating human GLI3 expression during early development

The zinc-finger transcription factor GLI3 acts as a primary transducer of Sonic hedgehog (Shh) signaling in a context-dependent combinatorial fashion. GLI3 participates in the patterning and growth of many organs, including the central nervous system (CNS) and limbs. Previously, we reported a subset of human intronic cis-regulators controlling many known aspects of endogenous Gli3 expression in mouse and zebrafish. Here we demonstrate in a transgenic zebrafish assay the potential of two novel tetrapod-teleost conserved non-coding elements (CNEs) docking within GLI3 intronic intervals (intron 3 and 4) to induce reporter gene expression at known sites of endogenous Gli3 transcription in embryonic domains such as the central nervous system (CNS) and limbs. Interestingly, the cell culture based assays reveal harmony with the context dependent dual nature of intra-GLI3 conserved elements. Furthermore, a transgenic zebrafish assay of previously reported limb-specific GLI3 transcriptional enhancers (previously tested in mice and chicken limb buds) induced reporter gene expression in zebrafish blood precursor cells and notochord instead of fin. These results demonstrate that the appendage-specific activity of a subset of GLI3-associated enhancers might be a tetrapod innovation. Taken together with our recent data, these results suggest that during the course of vertebrate evolution Gli3 expression control acquired a complex cis-regulatory landscape for spatiotemporal patterning of CNS and limbs. Comparative data from fish and mice suggest that the functional aspects of a subset of these cis-regulators have diverged significantly between these two lineages.

A Simple Predictive Enhancer Syntax for Hindbrain Patterning Is Conserved in Vertebrate Genomes

Background

Determining the function of regulatory elements is fundamental for our understanding of development, disease and evolution. However, the sequence features that mediate these functions are often unclear and the prediction of tissue-specific expression patterns from sequence alone is non-trivial. Previous functional studies have demonstrated a link between PBX-HOX and MEIS/PREP binding interactions and hindbrain enhancer activity, but the defining grammar of these sites, if any exists, has remained elusive.

Results

Here, we identify a shared sequence signature (syntax) within a heterogeneous set of conserved vertebrate hindbrain enhancers composed of spatially co-occurring PBX-HOX and MEIS/PREP transcription factor binding motifs. We use this syntax to accurately predict hindbrain enhancers in 89% of cases (67/75 predicted elements) from a set of conserved non-coding elements (CNEs). Furthermore, mutagenesis of the sites abolishes activity or generates ectopic expression, demonstrating their requirement for segmentally restricted enhancer activity in the hindbrain. We refine and use our syntax to predict over 3,000 hindbrain enhancers across the human genome. These sequences tend to be located near developmental transcription factors and are enriched in known hindbrain activating elements, demonstrating the predictive power of this simple model.

Conclusion

Our findings support the theory that hundreds of CNEs, and perhaps thousands of regions across the human genome, function to coordinate gene expression in the developing hindbrain. We speculate that deeply conserved sequences of this kind contributed to the co-option of new genes into the hindbrain gene regulatory network during early vertebrate evolution by linking patterns of hox expression to downstream genes involved in segmentation and patterning, and evolutionarily newer instances may have continued to contribute to lineage-specific elaboration of the hindbrain.

Submicroscopic deletions at 13q32.1 cause congenital microcoria

Congenital microcoria (MCOR) is a rare autosomal-dominant disorder characterized by inability of the iris to dilate owing to absence of dilator pupillae muscle. So far, a dozen MCOR-affected families have been reported worldwide. By using whole-genome oligonucleotide array CGH, we have identified deletions at 13q32.1 segregating with MCOR in six families originating from France, Japan, and Mexico. Breakpoint sequence analyses showed nonrecurrent deletions in 5/6 families. The deletions varied from 35 kbp to 80 kbp in size, but invariably encompassed or interrupted only two genes: TGDS encoding the TDP-glucose 4,6-dehydratase and GPR180 encoding the G protein-coupled receptor 180, also known as intimal thickness-related receptor (ITR). Unlike TGDS which has no known function in muscle cells, GPR180 is involved in the regulation of smooth muscle cell growth. The identification of a null GPR180 mutation segregating over two generations with iridocorneal angle dysgenesis, which can be regarded as a MCOR endophenotype, is consistent with the view that deletions of this gene, with or without the loss of elements regulating the expression of neighboring genes, are the cause of MCOR.

Divergent evolutionary rates in vertebrate and mammalian specific conserved non-coding elements (CNEs) in echolocating mammals

BACKGROUND

The majority of DNA contained within vertebrate genomes is non-coding, with a certain proportion of this thought to play regulatory roles during development. Conserved Non-coding Elements (CNEs) are an abundant group of putative regulatory sequences that are highly conserved across divergent groups and thus assumed to be under strong selective constraint. Many CNEs may contain regulatory factor binding sites, and their frequent spatial association with key developmental genes - such as those regulating sensory system development - suggests crucial roles in regulating gene expression and cellular patterning. Yet surprisingly little is known about the molecular evolution of CNEs across diverse mammalian taxa or their role in specific phenotypic adaptations. We examined 3,110 vertebrate-specific and ~82,000 mammalian-specific CNEs across 19 and 9 mammalian orders respectively, and tested for changes in the rate of evolution of CNEs located in the proximity of genes underlying the development or functioning of auditory systems. As we focused on CNEs putatively associated with genes underlying the development/functioning of auditory systems, we incorporated echolocating taxa in our dataset because of their highly specialised and derived auditory systems.

RESULTS

Phylogenetic reconstructions of concatenated CNEs broadly recovered accepted mammal relationships despite high levels of sequence conservation. We found that CNE substitution rates were highest in rodents and lowest in primates, consistent with previous findings. Comparisons of CNE substitution rates from several genomic regions containing genes linked to auditory system development and hearing revealed differences between echolocating and non-echolocating taxa. Wider taxonomic sampling of four CNEs associated with the homeobox genes Hmx2 and Hmx3 - which are required for inner ear development - revealed family-wise variation across diverse bat species. Specifically within one family of echolocating bats that utilise frequency-modulated echolocation calls varying widely in frequency and intensity high levels of sequence divergence were found.

CONCLUSIONS

Levels of selective constraint acting on CNEs differed both across genomic locations and taxa, with observed variation in substitution rates of CNEs among bat species. More work is needed to determine whether this variation can be linked to echolocation, and wider taxonomic sampling is necessary to fully document levels of conservation in CNEs across diverse taxa.

Abnormal dosage of ultraconserved elements is highly disfavored in healthy cells but not cancer cells

Ultraconserved elements (UCEs) are strongly depleted from segmental duplications and copy number variations (CNVs) in the human genome, suggesting that deletion or duplication of a UCE can be deleterious to the mammalian cell. Here we address the process by which CNVs become depleted of UCEs. We begin by showing that depletion for UCEs characterizes the most recent large-scale human CNV datasets and then find that even newly formed de novo CNVs, which have passed through meiosis at most once, are significantly depleted for UCEs. In striking contrast, CNVs arising specifically in cancer cells are, as a rule, not depleted for UCEs and can even become significantly enriched. This observation raises the possibility that CNVs that arise somatically and are relatively newly formed are less likely to have established a CNV profile that is depleted for UCEs. Alternatively, lack of depletion for UCEs from cancer CNVs may reflect the diseased state. In support of this latter explanation, somatic CNVs that are not associated with disease are depleted for UCEs. Finally, we show that it is possible to observe the CNVs of induced pluripotent stem (iPS) cells become depleted of UCEs over time, suggesting that depletion may be established through selection against UCE-disrupting CNVs without the requirement for meiotic divisions.

Dissection of a Ciona regulatory element reveals complexity of cross-species enhancer activity

Vertebrate genomes share numerous conserved non-coding elements, many of which function as enhancer elements and are hypothesised to be under evolutionary constraint due to a need to be bound by combinations of sequence-specific transcription factors. In contrast, few such conserved elements can be detected between vertebrates and their closest invertebrate relatives. Despite this lack of sequence identity, cross-species transgenesis has identified some cases where non-coding DNA from invertebrates drives reporter gene expression in transgenic vertebrates in patterns reminiscent of the expression of vertebrate orthologues. Such instances are presumed to reflect the presence of conserved suites of binding sites in the regulatory regions of invertebrate and vertebrate orthologues, such that both regulatory elements can correctly interpret the trans-activating environment. Shuffling of binding sites has been suggested to lie behind loss of sequence conservation; however this has not been experimentally tested. Here we examine the underlying basis of enhancer activity for the Ciona intestinalis βγ-crystallin gene, which drives expression in the lens of transgenic vertebrates despite the Ciona lineage predating the evolution of the lens. We construct an interactive gene regulatory network (GRN) for vertebrate lens development, allowing network interactions to be robustly catalogued and conserved network components and features to be identified. We show that a small number of binding motifs are necessary for Ciona βγ-crystallin expression, and narrow down the likely factors that bind to these motifs. Several of these overlap with the conserved core of the vertebrate lens GRN, implicating these sites in cross species function. However when we test these motifs in a transgenic vertebrate they prove to be dispensable for reporter expression in the lens. These results show that current models depicting cross species enhancer function as dependent on conserved binding sites can be overly simplistic, with sound evolutionary inference requiring detailed dissection of underlying mechanisms.

•

Analysis of binding motifs in a Ciona enhancer that also works in vertebrate lens.

•

Establishment of candidate transcription factors that may regulate this enhancer.

•

Construction of a curated, interactive gene regulatory network of lens development.

•

Public accessibility of this via a dedicated web site.

•

Experimental test of binding motif function in cross species transgenesis.

Capturing the regulatory interactions of eukaryote genomes

A key finding from early genomics research is the remarkable consistency in the number of protein-coding regions across diverse species. This has led many researchers to look to the cis-regulatory elements of genes as the fundamental influence behind evolving gene function and subsequent species diversification. Historically, since these elements are often located in vast intergenic and intronic regions of the genome, their identification has been recalcitrant. Now, with the deluge of whole-genome data from representatives of numerous eukaryotic lineages, various approaches have enabled us to begin to recognize features that characterize regulatory regions of the genome. Here we endeavour to collate these approaches in order to give an overview of the complexities involved in extrapolating regulatory signatures. The resource provided by the escalating richness of whole-genome datasets enables more sophisticated modelling of these regulatory signatures yet at the same time introduces increasing potential for noise. While we are only at the advent of making these discoveries, the next decade promises to be a very exciting and rewarding time for genome researchers.