A survey of ancient conserved non-coding elements in the PAX6 locus reveals a landscape of interdigitated cis-regulatory archipelagos

A bird's eye view on the use of whole exome sequencing in rare congenital ophthalmic diseases

Phenotypic and genotypic heterogeneity in congenital ocular diseases, especially in anterior segment dysgenesis (ASD), have created challenges for proper diagnosis and classification of diseases. Over the last decade, genomic research has indeed boosted our understanding in the molecular basis of ASD and genes associated with both autosomal dominant and recessive patterns of inheritance have been described with a wide range of expressivity. Here we describe the molecular characterization of a cohort of 162 patients displaying isolated or syndromic congenital ocular dysgenesis. Samples were analyzed with diverse techniques, such as direct sequencing, multiplex ligation-dependent probe amplification, and whole exome sequencing (WES), over 20 years. Our data reiterate the notion that PAX6 alterations are primarily associated with ASD, mostly aniridia, since the majority of the cohort (66.7%) has a pathogenic or likely pathogenic variant in the PAX6 locus. Unexpectedly, a high fraction of positive samples (20.3%) displayed deletions involving the 11p13 locus, either partially/totally involving PAX6 coding region or abolishing its critical regulatory region, underlying its significance. Most importantly, the use of WES has allowed us to both assess variants in known ASD genes (i.e., CYP1B1, ITPR1, MAB21L1, PXDN, and PITX2) and to identify rarer phenotypes (i.e., MIDAS, oculogastrointestinal-neurodevelopmental syndrome and Jacobsen syndrome). Our data clearly suggest that WES allows expanding the analytical portfolio of ocular dysgenesis, both isolated and syndromic, and that is pivotal for the differential diagnosis of those conditions in which there may be phenotypic overlaps and in general in ASD.

Short-read whole genome sequencing identifies causative variants in most individuals with previously unexplained aniridia

BACKGROUND

Classic aniridia is a highly penetrant autosomal dominant disorder characterised by congenital absence of the iris, foveal hypoplasia, optic disc anomalies and progressive opacification of the cornea. >90% of cases of classic aniridia are caused by heterozygous, loss-of-function variants affecting the PAX6 locus.

METHODS

Short-read whole genome sequencing was performed on 51 (39 affected) individuals from 37 different families who had screened negative for mutations in the PAX6 coding region.

RESULTS

Likely causative mutations were identified in 22 out of 37 (59%) families. In 19 out of 22 families, the causative genomic changes have an interpretable deleterious impact on the PAX6 locus. Of these 19 families, 1 has a novel heterozygous PAX6 frameshift variant missed on previous screens, 4 have single nucleotide variants (SNVs) (one novel) affecting essential splice sites of PAX6 5' non-coding exons and 2 have deep intronic SNV (one novel) resulting in gain of a donor splice site. In 12 out of 19, the causative variants are large-scale structural variants; 5 have partial or whole gene deletions of PAX6, 3 have deletions encompassing critical PAX6 cis-regulatory elements, 2 have balanced inversions with disruptive breakpoints within the PAX6 locus and 2 have complex rearrangements disrupting PAX6. The remaining 3 of 22 families have deletions encompassing FOXC1 (a known cause of atypical aniridia). Seven of the causative variants occurred de novo and one cosegregated with familial aniridia. We were unable to establish inheritance status in the remaining probands. No plausibly causative SNVs were identified in PAX6 cis-regulatory elements.

CONCLUSION

Whole genome sequencing proves to be an effective diagnostic test in most individuals with previously unexplained aniridia.

Unique activities of two overlapping PAX6 retinal enhancers

Enhancers play a critical role in development by precisely modulating spatial, temporal, and cell type-specific gene expression. Sequence variants in enhancers have been implicated in diseases; however, establishing the functional consequences of these variants is challenging because of a lack of understanding of precise cell types and developmental stages where the enhancers are normally active. PAX6 is the master regulator of eye development, with a regulatory landscape containing multiple enhancers driving the expression in the eye. Whether these enhancers perform additive, redundant or distinct functions is unknown. Here, we describe the precise cell types and regulatory activity of two PAX6 retinal enhancers, HS5 and NRE. Using a unique combination of live imaging and single-cell RNA sequencing in dual enhancer-reporter zebrafish embryos, we uncover differences in the spatiotemporal activity of these enhancers. Our results show that although overlapping, these enhancers have distinct activities in different cell types and therefore likely nonredundant functions. This work demonstrates that unique cell type-specific activities can be uncovered for apparently similar enhancers when investigated at high resolution in vivo.

A novel microdeletion of 517 kb downstream of the PAX6 gene in a Chinese family with congenital aniridia

BACKGROUND

To identify the disease-causing gene in a Chinese family affected with congenital aniridia.

METHODS

Patients underwent systematic ophthalmic examinations such as anterior segment photography, fundus photography, optical coherence tomography, and fundus fluorescein angiography. The proband was screened for pathogenic variants by whole exome sequencing (WES) and copy number variant (CNV) analysis. Real-time quantitative PCR (RT-qPCR) was applied to confirm the CNV results. Breakpoints were identified by long-range PCR followed by Sanger sequencing.

RESULTS

All seven members of this Chinese family, including four patients and three normal individuals, were recruited for this study. All patients showed bilateral congenital aniridia with nystagmus, except the son of the proband, who presented with bilateral partial coloboma of the iris. A novel heterozygous deletion (chr11:31,139,019-31,655,997) containing the 3' regulatory enhancers of the PAX6 gene was detected in this family. We also reviewed the reported microdeletions downstream of PAX6 in patients with aniridia.

CONCLUSIONS

We identified a novel microdeletion, 517 kb in size located about 133 kb downstream of the PAX6 gene, responsible for congenital aniridia in this Chinese family, which expands the spectrum of aniridia-associated mutations in PAX6.

Long-read genome sequencing identifies cryptic structural variants in congenital aniridia cases

BACKGROUND

Haploinsufficiency of the transcription factor PAX6 is the main cause of congenital aniridia, a genetic disorder characterized by iris and foveal hypoplasia. 11p13 microdeletions altering PAX6 or its downstream regulatory region (DRR) are present in about 25% of patients; however, only a few complex rearrangements have been described to date. Here, we performed nanopore-based whole-genome sequencing to assess the presence of cryptic structural variants (SVs) on the only two unsolved "PAX6-negative" cases from a cohort of 110 patients with congenital aniridia after unsuccessfully short-read sequencing approaches.

RESULTS

Long-read sequencing (LRS) unveiled balanced chromosomal rearrangements affecting the PAX6 locus at 11p13 in these two patients and allowed nucleotide-level breakpoint analysis. First, we identified a cryptic 4.9 Mb de novo inversion disrupting intron 7 of PAX6, further verified by targeted polymerase chain reaction amplification and sequencing and FISH-based cytogenetic analysis. Furthermore, LRS was decisive in correctly mapping a t(6;11) balanced translocation cytogenetically detected in a second proband with congenital aniridia and considered non-causal 15 years ago. LRS resolved that the breakpoint on chromosome 11 was indeed located at 11p13, disrupting the DNase I hypersensitive site 2 enhancer within the DRR of PAX6, 161 Kb from the causal gene. Patient-derived RNA expression analysis demonstrated PAX6 haploinsufficiency, thus supporting that the 11p13 breakpoint led to a positional effect by cleaving crucial enhancers for PAX6 transactivation. LRS analysis was also critical for mapping the exact breakpoint on chromosome 6 to the highly repetitive centromeric region at 6p11.1.

CONCLUSIONS

In both cases, the LRS-based identified SVs have been deemed the hidden pathogenic cause of congenital aniridia. Our study underscores the limitations of traditional short-read sequencing in uncovering pathogenic SVs affecting low-complexity regions of the genome and the value of LRS in providing insight into hidden sources of variation in rare genetic diseases.

Identification of clade-wide putative cis-regulatory elements from conserved non-coding sequences in Cucurbitaceae genomes

Cis-regulatory elements regulate gene expression and play an essential role in the development and physiology of organisms. Many conserved non-coding sequences (CNSs) function as cis-regulatory elements. They control the development of various lineages. However, predicting clade-wide cis-regulatory elements across several closely related species remains challenging. Based on the relationship between CNSs and cis-regulatory elements, we present a computational approach that predicts the clade-wide putative cis-regulatory elements in 12 Cucurbitaceae genomes. Using 12-way whole-genome alignment, we first obtained 632 112 CNSs in Cucurbitaceae. Next, we identified 16 552 Cucurbitaceae-wide cis-regulatory elements based on collinearity among all 12 Cucurbitaceae plants. Furthermore, we predicted 3 271 potential regulatory pairs in the cucumber genome, of which 98 were verified using integrative RNA sequencing and ChIP sequencing datasets from samples collected during various fruit development stages. The CNSs, Cucurbitaceae-wide cis-regulatory elements, and their target genes are accessible at http://cmb.bnu.edu.cn/cisRCNEs_cucurbit/. These elements are valuable resources for functionally annotating CNSs and their regulatory roles in Cucurbitaceae genomes.

Cell fate decisions, transcription factors and signaling during early retinal development

The development of the vertebrate eyes is a complex process starting from anterior-posterior and dorso-ventral patterning of the anterior neural tube, resulting in the formation of the eye field. Symmetrical separation of the eye field at the anterior neural plate is followed by two symmetrical evaginations to generate a pair of optic vesicles. Next, reciprocal invagination of the optic vesicles with surface ectoderm-derived lens placodes generates double-layered optic cups. The inner and outer layers of the optic cups develop into the neural retina and retinal pigment epithelium (RPE), respectively. In vitro produced retinal tissues, called retinal organoids, are formed from human pluripotent stem cells, mimicking major steps of retinal differentiation in vivo. This review article summarizes recent progress in our understanding of early eye development, focusing on the formation the eye field, optic vesicles, and early optic cups. Recent single-cell transcriptomic studies are integrated with classical in vivo genetic and functional studies to uncover a range of cellular mechanisms underlying early eye development. The functions of signal transduction pathways and lineage-specific DNA-binding transcription factors are dissected to explain cell-specific regulatory mechanisms underlying cell fate determination during early eye development. The functions of homeodomain (HD) transcription factors Otx2, Pax6, Lhx2, Six3 and Six6, which are required for early eye development, are discussed in detail. Comprehensive understanding of the mechanisms of early eye development provides insight into the molecular and cellular basis of developmental ocular anomalies, such as optic cup coloboma. Lastly, modeling human development and inherited retinal diseases using stem cell-derived retinal organoids generates opportunities to discover novel therapies for retinal diseases.

A Journey Through Genetics to Biology

Although my engagement with human genetics emerged gradually, and sometimes serendipitously, it has held me spellbound for decades. Without my teachers, students, postdocs, colleagues, and collaborators, I would not be writing this review of my scientific adventures. Early gene and disease mapping was a satisfying puzzle-solving exercise, but building biological insight was my main goal. The project trajectory was hugely influenced by the evolutionarily conserved nature of the implicated genes and by the pace of progress in genetic technologies. The rich detail of clinical observations, particularly in eye disease, makes humans an excellent model, especially when complemented by the use of multiple other animal species for experimental validation. The contributions of collaborators and rivals also influenced our approach. We are very fortunate to work in this era of unprecedented progress in genetics and genomics. Expected final online publication date for the Annual Review of Genomics and Human Genetics, Volume 23 is October 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Genetics and epidemiology of aniridia: Updated guidelines for genetic study

Aniridia is a panocular disease characterized by iris hypoplasia, accompanied by other ocular manifestations, with a high clinical variability and overlapping with different abnormalities of the anterior and posterior segment. This review focuses on the genetic features of this autosomal dominant pathology, which is caused by the haploinsufficiency of the PAX6 gene. Mutations causing premature stop codons are the most frequent among the wider mutational spectrum of PAX6, with more than 600 different mutations identified so far. Recent advances in next-generation sequencing (NGS) have increased the diagnostic yield in aniridia and contributed to elucidate new etiopathogenic mechanisms leading to PAX6 haploinsufficiency. Here, we also update good practices and recommendations to improve genetic testing and clinical management of aniridia using more cost-effective NGS analysis. Those new approaches also allow studying simultaneously both structural variants and point-mutations in PAX6 as well as other genes for differential diagnosis, simultaneously. Some patients with atypical phenotypes might present mutations in FOXC1 and PITX2, both genes causing a wide spectrum of anterior segment dysgenesis, or in ITPR1, which is responsible for a distinctive form of circumpupillary iris aplasia present in Gillespie syndrome, or other mutations in minor genes. Since aniridia can also associate extraocular anomalies, as it occurs in carriers of PAX6 and WT1 microdeletions leading to WAGR syndrome, genetic studies are crucial to assure a correct diagnosis and clinical management, besides allowing prenatal and preimplantational genetic testing in families.

dbCNS: A New Database for Conserved Noncoding Sequences

We developed dbCNS (http://yamasati.nig.ac.jp/dbcns), a new database for conserved noncoding sequences (CNSs). CNSs exist in many eukaryotes and are assumed to be involved in protein expression control. Version 1 of dbCNS, introduced here, includes a powerful and precise CNS identification pipeline for multiple vertebrate genomes. Mutations in CNSs may induce morphological changes and cause genetic diseases. For this reason, many vertebrate CNSs have been identified, with special reference to primate genomes. We integrated ∼6.9 million CNSs from many vertebrate genomes into dbCNS, which allows users to extract CNSs near genes of interest using keyword searches. In addition to CNSs, dbCNS contains published genome sequences of 161 species. With purposeful taxonomic sampling of genomes, users can employ CNSs as queries to reconstruct CNS alignments and phylogenetic trees, to evaluate CNS modifications, acquisitions, and losses, and to roughly identify species with CNSs having accelerated substitution rates. dbCNS also produces links to dbSNP for searching pathogenic single-nucleotide polymorphisms in human CNSs. Thus, dbCNS connects morphological changes with genetic diseases. A test analysis using 38 gnathostome genomes was accomplished within 30 s. dbCNS results can evaluate CNSs identified by other stand-alone programs using genome-scale data.

Genetics and epidemiology of aniridia: Updated guidelines for genetic study

Aniridia is a panocular disease characterized by iris hypoplasia, accompanied by other ocular manifestations, with a high clinical variability and overlapping with different abnormalities of the anterior and posterior segment. This review focuses on the genetic features of this autosomal dominant pathology, which is caused by the haploinsufficiency of the PAX6 gene. Mutations causing premature stop codons are the most frequent among the wider mutational spectrum of PAX6, with more than 600 different mutations identified so far. Recent advances in next-generation sequencing (NGS) have increased the diagnostic yield in aniridia and contributed to elucidate new etiopathogenic mechanisms leading to PAX6 haploinsufficiency. Here, we also update good practices and recommendations to improve genetic testing and clinical management of aniridia using more cost-effective NGS analysis. Those new approaches also allow studying simultaneously both structural variants and point-mutations in PAX6 as well as other genes for differential diagnosis, simultaneously. Some patients with atypical phenotypes might present mutations in FOXC1 and PITX2, both genes causing a wide spectrum of anterior segment dysgenesis, or in ITPR1, which is responsible for a distinctive form of circumpupillary iris aplasia present in Gillespie syndrome, or other mutations in minor genes. Since aniridia can also associate extraocular anomalies, as it occurs in carriers of PAX6 and WT1 microdeletions leading to WAGR syndrome, genetic studies are crucial to assure a correct diagnosis and clinical management, besides allowing prenatal and preimplantational genetic testing in families.

A novel measure of non-coding genome conservation identifies genomic regulatory blocks within primates

MOTIVATION

Clusters of extremely conserved non-coding elements (CNEs) mark genomic regions devoted to cis-regulation of key developmental genes in Metazoa. We have recently shown that their span coincides with that of topologically associating domains (TADs), making them useful for estimating conserved TAD boundaries in the absence of Hi-C data. The standard approach-detecting CNEs in genome alignments and then establishing the boundaries of their clusters-requires tuning of several parameters and breaks down when comparing closely related genomes.

RESULTS

We present a novel, kurtosis-based measure of pairwise non-coding conservation that requires no pre-set thresholds for conservation level and length of CNEs. We show that it performs robustly across a large span of evolutionary distances, including across the closely related genomes of primates for which standard approaches fail. The method is straightforward to implement and enables detection and comparison of clusters of CNEs and estimation of underlying TADs across a vastly increased range of Metazoan genomes.

AVAILABILITY AND IMPLEMENTATION

The data generated for this study, and the scripts used to generate the data, can be found at https://github.com/alexander-nash/kurtosis_conservation.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Understanding the genetics of neuropsychiatric disorders: the potential role of genomic regulatory blocks

Recent genome-wide association studies have identified numerous loci associated with neuropsychiatric disorders. The majority of these are in non-coding regions, and are commonly assigned to the nearest gene along the genome. However, this approach neglects the three-dimensional organisation of the genome, and the fact that the genome contains arrays of extremely conserved non-coding elements termed genomic regulatory blocks (GRBs), which can be utilized to detect genes under long-range developmental regulation. Here we review a GRB-based approach to assign loci in non-coding regions to potential target genes, and apply it to reanalyse the results of one of the largest schizophrenia GWAS (SWG PGC, 2014). We further apply this approach to GWAS data from two related neuropsychiatric disorders-autism spectrum disorder and bipolar disorder-to show that it is applicable to developmental disorders in general. We find that disease-associated SNPs are overrepresented in GRBs and that the GRB model is a powerful tool for linking these SNPs to their correct target genes under long-range regulation. Our analysis identifies novel genes not previously implicated in schizophrenia and corroborates a number of predicted targets from the original study. The results are available as an online resource in which the genomic context and the strength of enhancer-promoter associations can be browsed for each schizophrenia-associated SNP.

Lampreys, the jawless vertebrates, contain three Pax6 genes with distinct expression in eye, brain and pancreas

The transcription factor Pax6 is crucial for the development of the central nervous system, eye, olfactory system and pancreas, and is implicated in human disease. While a single Pax6 gene exists in human and chicken, Pax6 occurs as a gene family in other vertebrates, with two members in elephant shark, Xenopus tropicalis and Anolis lizard and three members in teleost fish such as stickleback and medaka. However, the complement of Pax6 genes in jawless vertebrates (cyclostomes), the sister group of jawed vertebrates (gnathostomes), is unknown. Using a combination of BAC sequencing and genome analysis, we discovered three Pax6 genes in lampreys. Unlike the paired-less Pax6 present in some gnathostomes, all three lamprey Pax6 have a highly conserved full-length paired domain. All three Pax6 genes are expressed in the eye and brain, with variable expression in other tissues. Notably, lamprey Pax6α transcripts are found in the pancreas, a vertebrate-specific organ, indicating the involvement of Pax6 in development of the pancreas in the vertebrate ancestor. Multi-species sequence comparisons revealed only a single conserved non-coding element, in the lamprey Pax6β locus, with similarity to the PAX6 neuroretina enhancer. Using a transgenic zebrafish enhancer assay we demonstrate functional conservation of this element over 500 million years of vertebrate evolution.

The Spectrum of PAX6 Mutations and Genotype-Phenotype Correlations in the Eye

The transcription factor PAX6 is essential in ocular development in vertebrates, being considered the master regulator of the eye. During eye development, it is essential for the correct patterning and formation of the multi-layered optic cup and it is involved in the developing lens and corneal epithelium. In adulthood, it is mostly expressed in cornea, iris, and lens. PAX6 is a dosage-sensitive gene and it is highly regulated by several elements located upstream, downstream, and within the gene. There are more than 500 different mutations described to affect PAX6 and its regulatory regions, the majority of which lead to PAX6 haploinsufficiency, causing several ocular and systemic abnormalities. Aniridia is an autosomal dominant disorder that is marked by the complete or partial absence of the iris, foveal hypoplasia, and nystagmus, and is caused by heterozygous PAX6 mutations. Other ocular abnormalities have also been associated with PAX6 changes, and genotype-phenotype correlations are emerging. This review will cover recent advancements in PAX6 regulation, particularly the role of several enhancers that are known to regulate PAX6 during eye development and disease. We will also present an updated overview of the mutation spectrum, where an increasing number of mutations in the non-coding regions have been reported. Novel genotype-phenotype correlations will also be discussed.

Zebrafish: A Powerful Model for Understanding the Functional Relevance of Noncoding Region Mutations in Human Genetic Diseases

Determining aetiology of genetic disorders caused by damaging mutations in protein-coding genes is well established. However, understanding how mutations in the vast stretches of the noncoding genome contribute to genetic abnormalities remains a huge challenge. Cis-regulatory elements (CREs) or enhancers are an important class of noncoding elements. CREs function as the primary determinants of precise spatial and temporal regulation of their target genes during development by serving as docking sites for tissue-specific transcription factors. Although a large number of potential disease-associated CRE mutations are being identified in patients, lack of robust methods for mechanistically linking these mutations to disease phenotype is currently hampering the understanding of their roles in disease aetiology. Here, we have described the various systems available for testing the CRE potential of stretches of noncoding regions harbouring mutations implicated in human disease. We highlight advances in the field leading to the establishment of zebrafish as a powerful system for robust and cost-effective functional assays of CRE activity, enabling rapid identification of causal variants in regulatory regions and the validation of their role in disruption of appropriate gene expression.

Implication of non-coding PAX6 mutations in aniridia

There is an increasing implication of non-coding regions in pathological processes of genetic origin. This is partly due to the emergence of sophisticated techniques that have transformed research into gene expression by allowing a more global understanding of the genome, both at the genomic, epigenomic and chromatin levels. Here, we implemented the analysis of PAX6, whose coding loss-of-function variants are mainly implied in aniridia, by studying its non-coding regions (untranslated regions, introns and cis-regulatory sequences). In particular, we have taken advantage of the development of high-throughput approaches to screen the upstream and downstream regulatory regions of PAX6 in 47 aniridia patients without identified mutation in the coding sequence. This was made possible through the use of custom targeted resequencing and/or CGH array to analyze the entire PAX6 locus on 11p13. We found candidate variants in 30 of the 47 patients. 9/30 correspond to the well-known described 3' deletions encompassing SIMO and other enhancer elements. In addition, we identified numerous different variants in various non-coding regions, in particular untranslated regions. Among these latter, most of them demonstrated an in vitro functional effect using a minigene strategy, and 12/21 are thus considered as causative mutations or very likely to explain the phenotypes. This new analysis strategy brings molecular diagnosis to more than 90% of our aniridia patients. This study revealed an outstanding mutation pattern in non-coding PAX6 regions confirming that PAX6 remains the major gene for aniridia.

A long range distal enhancer controls temporal fine-tuning of PAX6 expression in neuronal precursors

Proper embryonic development relies on a tight control of spatial and temporal gene expression profiles in a highly regulated manner. One good example is the ON/OFF switching of the transcription factor PAX6 that governs important steps of neurogenesis. In the neural tube PAX6 expression is initiated in neural progenitors through the positive action of retinoic acid signaling and downregulated in neuronal precursors by the bHLH transcription factor NEUROG2. How these two regulatory inputs are integrated at the molecular level to properly fine tune temporal PAX6 expression is not known. In this study we identified and characterized a 940-bp long distal cis-regulatory module (CRM), located far away from the PAX6 transcription unit and which conveys positive input from RA signaling pathway and indirect repressive signal(s) from NEUROG2. These opposing regulatory signals are integrated through HOMZ, a 94 bp core region within E940 which is evolutionarily conserved in distant organisms such as the zebrafish. We show that within HOMZ, NEUROG2 and RA exert their opposite temporal activities through a short 60 bp region containing a functional RA-responsive element (RARE). We propose a model in which retinoic acid receptors (RARs) and NEUROG2 repressive target(s) compete on the same DNA motif to fine tune temporal PAX6 expression during the course of spinal neurogenesis.

Conserved noncoding sequences conserve biological networks and influence genome evolution

Comparative genomics approaches have identified numerous conserved cis-regulatory sequences near genes in plant genomes. Despite the identification of these conserved noncoding sequences (CNSs), our knowledge of their functional importance and selection remains limited. Here, we used a combination of DNA methylome analysis, microarray expression analyses, and functional annotation to study these sequences in the model tree Populus trichocarpa. Methylation in CG contexts and non-CG contexts was lower in CNSs, particularly CNSs in the 5'-upstream regions of genes, compared with other sites in the genome. We observed that CNSs are enriched in genes with transcription and binding functions, and this also associated with syntenic genes and those from whole-genome duplications, suggesting that cis-regulatory sequences play a key role in genome evolution. We detected a significant positive correlation between CNS number and protein interactions, suggesting that CNSs may have roles in the evolution and maintenance of biological networks. The divergence of CNSs indicates that duplication-degeneration-complementation drives the subfunctionalization of a proportion of duplicated genes from whole-genome duplication. Furthermore, population genomics confirmed that most CNSs are under strong purifying selection and only a small subset of CNSs shows evidence of adaptive evolution. These findings provide a foundation for future studies exploring these key genomic features in the maintenance of biological networks, local adaptation, and transcription.

Topologically associating domains are ancient features that coincide with Metazoan clusters of extreme noncoding conservation

Developmental genes in metazoan genomes are surrounded by dense clusters of conserved noncoding elements (CNEs). CNEs exhibit unexplained extreme levels of sequence conservation, with many acting as developmental long-range enhancers. Clusters of CNEs define the span of regulatory inputs for many important developmental regulators and have been described previously as genomic regulatory blocks (GRBs). Their function and distribution around important regulatory genes raises the question of how they relate to 3D conformation of these loci. Here, we show that clusters of CNEs strongly coincide with topological organisation, predicting the boundaries of hundreds of topologically associating domains (TADs) in human and Drosophila. The set of TADs that are associated with high levels of noncoding conservation exhibit distinct properties compared to TADs devoid of extreme noncoding conservation. The close correspondence between extreme noncoding conservation and TADs suggests that these TADs are ancient, revealing a regulatory architecture conserved over hundreds of millions of years.

Metazoan genomes contain many clusters of conserved noncoding elements. Here, the authors provide evidence that these clusters coincide with distinct topologically associating domains in humans and Drosophila, revealing a conserved regulatory genomic architecture.

Signaling and Gene Regulatory Networks in Mammalian Lens Development

Ocular lens development represents an advantageous system in which to study regulatory mechanisms governing cell fate decisions, extracellular signaling, cell and tissue organization, and the underlying gene regulatory networks. Spatiotemporally regulated domains of BMP, FGF, and other signaling molecules in late gastrula-early neurula stage embryos generate the border region between the neural plate and non-neural ectoderm from which multiple cell types, including lens progenitor cells, emerge and undergo initial tissue formation. Extracellular signaling and DNA-binding transcription factors govern lens and optic cup morphogenesis. Pax6, c-Maf, Hsf4, Prox1, Sox1, and a few additional factors regulate the expression of the lens structural proteins, the crystallins. Extensive crosstalk between a diverse array of signaling pathways controls the complexity and order of lens morphogenetic processes and lens transparency.

Lampreys, the jawless vertebrates, contain only two ParaHox gene clusters

ParaHox genes (Gsx, Pdx, and Cdx) are an ancient family of developmental genes closely related to the Hox genes. They play critical roles in the patterning of brain and gut. The basal chordate, amphioxus, contains a single ParaHox cluster comprising one member of each family, whereas nonteleost jawed vertebrates contain four ParaHox genomic loci with six or seven ParaHox genes. Teleosts, which have experienced an additional whole-genome duplication, contain six ParaHox genomic loci with six ParaHox genes. Jawless vertebrates, represented by lampreys and hagfish, are the most ancient group of vertebrates and are crucial for understanding the origin and evolution of vertebrate gene families. We have previously shown that lampreys contain six Hox gene loci. Here we report that lampreys contain only two ParaHox gene clusters (designated as α- and β-clusters) bearing five ParaHox genes (Gsxα, Pdxα, Cdxα, Gsxβ, and Cdxβ). The order and orientation of the three genes in the α-cluster are identical to that of the single cluster in amphioxus. However, the orientation of Gsxβ in the β-cluster is inverted. Interestingly, Gsxβ is expressed in the eye, unlike its homologs in jawed vertebrates, which are expressed mainly in the brain. The lamprey Pdxα is expressed in the pancreas similar to jawed vertebrate Pdx genes, indicating that the pancreatic expression of Pdx was acquired before the divergence of jawless and jawed vertebrate lineages. It is likely that the lamprey Pdxα plays a crucial role in pancreas specification and insulin production similar to the Pdx of jawed vertebrates.

The Gene Regulatory Network of Lens Induction Is Wired through Meis-Dependent Shadow Enhancers of Pax6

Lens induction is a classical developmental model allowing investigation of cell specification, spatiotemporal control of gene expression, as well as how transcription factors are integrated into highly complex gene regulatory networks (GRNs). Pax6 represents a key node in the gene regulatory network governing mammalian lens induction. Meis1 and Meis2 homeoproteins are considered as essential upstream regulators of Pax6 during lens morphogenesis based on their interaction with the ectoderm enhancer (EE) located upstream of Pax6 transcription start site. Despite this generally accepted regulatory pathway, Meis1-, Meis2- and EE-deficient mice have surprisingly mild eye phenotypes at placodal stage of lens development. Here, we show that simultaneous deletion of Meis1 and Meis2 in presumptive lens ectoderm results in arrested lens development in the pre-placodal stage, and neither lens placode nor lens is formed. We found that in the presumptive lens ectoderm of Meis1/Meis2 deficient embryos Pax6 expression is absent. We demonstrate using chromatin immunoprecipitation (ChIP) that in addition to EE, Meis homeoproteins bind to a remote, ultraconserved SIMO enhancer of Pax6. We further show, using in vivo gene reporter analyses, that the lens-specific activity of SIMO enhancer is dependent on the presence of three Meis binding sites, phylogenetically conserved from man to zebrafish. Genetic ablation of EE and SIMO enhancers demostrates their requirement for lens induction and uncovers an apparent redundancy at early stages of lens development. These findings identify a genetic requirement for Meis1 and Meis2 during the early steps of mammalian eye development. Moreover, they reveal an apparent robustness in the gene regulatory mechanism whereby two independent "shadow enhancers" maintain critical levels of a dosage-sensitive gene, Pax6, during lens induction.

PAX6 MiniPromoters drive restricted expression from rAAV in the adult mouse retina

Current gene therapies predominantly use small, strong, and readily available ubiquitous promoters. However, as the field matures, the availability of small, cell-specific promoters would be greatly beneficial. Here we design seven small promoters from the human paired box 6 (PAX6) gene and test them in the adult mouse retina using recombinant adeno-associated virus. We chose the retina due to previous successes in gene therapy for blindness, and the PAX6 gene since it is: well studied; known to be driven by discrete regulatory regions; expressed in therapeutically interesting retinal cell types; and mutated in the vision-loss disorder aniridia, which is in need of improved therapy. At the PAX6 locus, 31 regulatory regions were bioinformatically predicted, and nine regulatory regions were constructed into seven MiniPromoters. Driving Emerald GFP, these MiniPromoters were packaged into recombinant adeno-associated virus, and injected intravitreally into postnatal day 14 mice. Four MiniPromoters drove consistent retinal expression in the adult mouse, driving expression in combinations of cell-types that endogenously express Pax6: ganglion, amacrine, horizontal, and Müller glia. Two PAX6-MiniPromoters drive expression in three of the four cell types that express PAX6 in the adult mouse retina. Combined, they capture all four cell types, making them potential tools for research, and PAX6-gene therapy for aniridia.

PAX6: 25th anniversary and more to learn

The DNA-binding transcription factor PAX6 was cloned 25 years ago by multiple teams pursuing identification of human and mouse eye disease causing genes, cloning vertebrate homologues of pattern-forming regulatory genes identified in Drosophila, or abundant eye-specific transcripts. Since its discovery in 1991, genetic, cellular, molecular and evolutionary studies on Pax6 mushroomed in the mid 1990s leading to the transformative thinking regarding the genetic program orchestrating both early and late stages of eye morphogenesis as well as the origin and evolution of diverse visual systems. Since Pax6 is also expressed outside of the eye, namely in the central nervous system and pancreas, a number of important insights into the development and function of these organs have been amassed. In most recent years, genome-wide technologies utilizing massively parallel DNA sequencing have begun to provide unbiased insights into the regulatory hierarchies of specification, determination and differentiation of ocular cells and neurogenesis in general. This review is focused on major advancements in studies on mammalian eye development driven by studies of Pax6 genes in model organisms and future challenges to harness the technology-driven opportunities to reconstruct, step-by-step, the transition from naïve ectoderm, neuroepithelium and periocular mesenchyme/neural crest cells into the three-dimensional architecture of the eye.

Genomic Locations of Conserved Noncoding Sequences and Their Proximal Protein-Coding Genes in Mammalian Expression Dynamics

Experimental studies have found the involvement of certain conserved noncoding sequences (CNSs) in the regulation of the proximal protein-coding genes in mammals. However, reported cases of long range enhancer activities and inter-chromosomal regulation suggest that proximity of CNSs to protein-coding genes might not be important for regulation. To test the importance of the CNS genomic location, we extracted the CNSs conserved between chicken and four mammalian species (human, mouse, dog, and cattle). These CNSs were confirmed to be under purifying selection. The intergenic CNSs are often found in clusters in gene deserts, where protein-coding genes are in paucity. The distribution pattern, ChIP-Seq, and RNA-Seq data suggested that the CNSs are more likely to be regulatory elements and not corresponding to long intergenic noncoding RNAs. Physical distances between CNS and their nearest protein coding genes were well conserved between human and mouse genomes, and CNS-flanking genes were often found in evolutionarily conserved genomic neighborhoods. ChIP-Seq signal and gene expression patterns also suggested that CNSs regulate nearby genes. Interestingly, genes with more CNSs have more evolutionarily conserved expression than those with fewer CNSs. These computationally obtained results suggest that the genomic locations of CNSs are important for their regulatory functions. In fact, various kinds of evolutionary constraints may be acting to maintain the genomic locations of CNSs and protein-coding genes in mammals to ensure proper regulation.

ChromContact: A web tool for analyzing spatial contact of chromosomes from Hi-C data

Background

Hi-C analysis has revealed the three-dimensional architecture of chromosomes in the nucleus. Although Hi-C data contains valuable information on long-range interactions of chromosomes, the data is not yet widely utilized by molecular biologists because of the quantity of data.

Results

We developed a web tool, ChromContact, to utilize the information obtained by Hi-C. The web tool is designed to be simple and easy to use. By specifying a locus of interest, ChromContact calculates contact profiles and generates links to the UCSC Genome Browser, enabling users to visually examine the contact information with various annotations.

Conclusion

ChromContact provides wide-range of molecular biologists with a user-friendly means to access high-resolution Hi-C data. One of the possible applications of ChromContact is investigating novel long-range promoter-enhancer interactions. This facilitates the functional interpretation of statistically significant markers identified by GWAS or ChIP-seq peaks that are located far from any annotated genes. ChromContact is freely accessible at http://bioinfo.sls.kyushu-u.ac.jp/chromcontact/.

The identification of cis-regulatory elements: A review from a machine learning perspective

The majority of the human genome consists of non-coding regions that have been called junk DNA. However, recent studies have unveiled that these regions contain cis-regulatory elements, such as promoters, enhancers, silencers, insulators, etc. These regulatory elements can play crucial roles in controlling gene expressions in specific cell types, conditions, and developmental stages. Disruption to these regions could contribute to phenotype changes. Precisely identifying regulatory elements is key to deciphering the mechanisms underlying transcriptional regulation. Cis-regulatory events are complex processes that involve chromatin accessibility, transcription factor binding, DNA methylation, histone modifications, and the interactions between them. The development of next-generation sequencing techniques has allowed us to capture these genomic features in depth. Applied analysis of genome sequences for clinical genetics has increased the urgency for detecting these regions. However, the complexity of cis-regulatory events and the deluge of sequencing data require accurate and efficient computational approaches, in particular, machine learning techniques. In this review, we describe machine learning approaches for predicting transcription factor binding sites, enhancers, and promoters, primarily driven by next-generation sequencing data. Data sources are provided in order to facilitate testing of novel methods. The purpose of this review is to attract computational experts and data scientists to advance this field.

PAX3 and PAX7 as upstream regulators of myogenesis

Like other subclasses within the PAX transcription factor family, PAX3 and PAX7 play important roles in the emergence of a number of different tissues during development. PAX3 regulates neural crest and, together with its orthologue PAX7, is also expressed in parts of the central nervous system. In this chapter we will focus on their role in skeletal muscle. Both factors are key regulators of myogenesis where Pax3 plays a major role during early skeletal muscle formation in the embryo while Pax7 predominates during post-natal growth and muscle regeneration in the adult. We review the expression and functions of these factors in the myogenic context. We also discuss mechanistic aspects of PAX3/7 function and modulation of their activity by interaction with other proteins, as well as the post-transcriptional and transcriptional regulation of their expression.

Transcriptional and epigenetic mechanisms of early cortical development: An examination of how Pax6 coordinates cortical development

The development of the cortex is an elaborate process that integrates a plethora of finely tuned molecular processes ranging from carefully regulated gradients of transcription factors, dynamic changes in the chromatin landscape, or formation of protein complexes to elicit and regulate transcription. Combined with cellular processes such as cell type specification, proliferation, differentiation, and migration, all of these developmental processes result in the establishment of an adult mammalian cortex with its typical lamination and regional patterning. By examining in-depth the role of one transcription factor, Pax6, on the regulation of cortical development, its integration in the regulation of chromatin state, and its regulation by cis-regulatory elements, we aim to demonstrate the importance of integrating each level of regulation in our understanding of cortical development.

Chromatin Dynamics in Lineage Commitment and Cellular Reprogramming

Dynamic structural properties of chromatin play an essential role in defining cell identity and function. Transcription factors and chromatin modifiers establish and maintain cell states through alteration of DNA accessibility and histone modifications. This activity is focused at both gene-proximal promoter regions and distally located regulatory elements. In the three-dimensional space of the nucleus, distal elements are localized in close physical proximity to the gene-proximal regulatory sequences through the formation of chromatin loops. These looping features in the genome are highly dynamic as embryonic stem cells differentiate and commit to specific lineages, and throughout reprogramming as differentiated cells reacquire pluripotency. Identifying these functional distal regulatory regions in the genome provides insight into the regulatory processes governing early mammalian development and guidance for improving the protocols that generate induced pluripotent cells.

Pax2/5/8 and Pax6 alternative splicing events in basal chordates and vertebrates: a focus on paired box domain

Paired box transcription factors play important role in development and tissue morphogenesis. The number of Pax homologs varies among species studied so far, due to genome and gene duplications that have affected PAX family to a great extent. Based on sequence similarity and functional domains, four Pax classes have been identified in chordates, namely Pax1/9, Pax2/5/8, Pax3/7, and Pax4/6. Numerous splicing events have been reported mainly for Pax2/5/8 and Pax6 genes. Of significant interest are those events that lead to Pax proteins with presumed novel properties, such as altered DNA-binding or transcriptional activity. In the current study, a thorough analysis of Pax2/5/8 splicing events from cephalochordates and vertebrates was performed. We focused more on Pax2/5/8 and Pax6 splicing events in which the paired domain is involved. Three new splicing events were identified in Oryzias latipes, one of which seems to be conserved in Acanthomorphata. Using representatives from deuterostome and protostome phyla, a comparative analysis of the Pax6 exon-intron structure of the paired domain was performed, during an attempt to estimate the time of appearance of the Pax6(5a) mRNA isoform. As shown in our analysis, this splicing event is characteristic of Gnathostomata and is absent in the other chordate subphyla. Moreover, expression pattern of alternative spliced variants was compared between cephalochordates and fish species. In summary, our data indicate expansion of alternative mRNA variants in paired box region of Pax2/5/8 and Pax6 genes during the course of vertebrate evolution.