Gillespie syndrome (GS) is characterized by bilateral iris hypoplasia, congenital hypotonia, non-progressive ataxia, and progressive cerebellar atrophy. Trio-based exome sequencing identified de novo mutations in ITPR1 in three unrelated individuals with GS recruited to the Deciphering Developmental Disorders study. Whole-exome or targeted sequence analysis identified plausible disease-causing ITPR1 mutations in 10/10 additional GS-affected individuals. These ultra-rare protein-altering variants affected only three residues in ITPR1: Glu2094 missense (one de novo, one co-segregating), Gly2539 missense (five de novo, one inheritance uncertain), and Lys2596 in-frame deletion (four de novo). No clinical or radiological differences were evident between individuals with different mutations. ITPR1 encodes an inositol 1,4,5-triphosphate-responsive calcium channel. The homo-tetrameric structure has been solved by cryoelectron microscopy. Using estimations of the degree of structural change induced by known recessive- and dominant-negative mutations in other disease-associated multimeric channels, we developed a generalizable computational approach to indicate the likely mutational mechanism. This analysis supports a dominant-negative mechanism for GS variants in ITPR1. In GS-derived lymphoblastoid cell lines (LCLs), the proportion of ITPR1-positive cells using immunofluorescence was significantly higher in mutant than control LCLs, consistent with an abnormality of nuclear calcium signaling feedback control. Super-resolution imaging supports the existence of an ITPR1-lined nucleoplasmic reticulum. Mice with Itpr1 heterozygous null mutations showed no major iris defects. Purkinje cells of the cerebellum appear to be the most sensitive to impaired ITPR1 function in humans. Iris hypoplasia is likely to result from either complete loss of ITPR1 activity or structure-specific disruption of multimeric interactions.

Biological differences between cell types and developmental processes are characterised by differences in gene expression profiles. Gene-distal enhancers are key components of the regulatory networks that specify the tissue-specific expression patterns driving embryonic development and cell fate decisions, and variations in their sequences are a major contributor to genetic disease and disease susceptibility. Despite advances in the methods for discovery of putative cis-regulatory sequences, characterisation of their spatio-temporal enhancer activities in a mammalian model system remains a major bottle-neck. We employed a strategy that combines gnathostome sequence conservation with transgenic mouse and zebrafish reporter assays to survey the genomic locus of the developmental control gene PAX6 for the presence of novel cis-regulatory elements. Sequence comparison between human and the cartilaginous elephant shark (Callorhinchus milii) revealed several ancient gnathostome conserved non-coding elements (agCNEs) dispersed widely throughout the PAX6 locus, extending the range of the known PAX6 cis-regulatory landscape to contain the full upstream PAX6-RCN1 intergenic region. Our data indicates that ancient conserved regulatory sequences can be tested effectively in transgenic zebrafish even when not conserved in zebrafish themselves. The strategy also allows efficient dissection of compound regulatory regions previously assessed in transgenic mice. Remarkable overlap in expression patterns driven by sets of agCNEs indicates that PAX6 resides in a landscape of multiple tissue-specific regulatory archipelagos.

Pax6 is a key transcriptional regulator in eye, olfactory system, forebrain, pituitary cerebellum, spinal cord and pancreas development. Alternative splicing, promoter usage and multiple enhancers regulate the complex Pax6 spatio-temporal expression pattern. Chromosomal rearrangements which abolish PAX6 gene expression have been characterised downstream of the coding region. Through evolutionary sequence comparison and transgenic reporter studies, we have identified a new Pax6 3' cis-regulatory region. This region, C1170 Box 123, contains three distinct modules of human-mouse sequence conservation, while only Box 1 is conserved to Fugu. Both the human and the orthologous Fugu sequence direct similar reporter gene expression in the developing pretectum, neural retina and olfactory region, indicating evolutionary conservation of Pax6 regulatory mechanisms despite the low level of overall sequence conservation.

Heterozygous loss-of-function (LOF) mutations in the gene encoding the DNA-binding protein, SATB2, result in micrognathia and cleft palate in both humans and mice. In three unrelated individuals, we show that translocation breakpoints (BPs) up to 896 kb 3' of SATB2 polyadenylation site cause a phenotype which is indistinguishable from that caused by SATB2 LOF mutations. This syndrome comprises long nose, small mouth, micrognathia, cleft palate, arachnodactyly and intellectual disability. These BPs map to a gene desert between PLCL1 and SATB2. We identified three putative cis-regulatory elements (CRE1-3) using a comparative genomic approach each of which would be placed in trans relative to SATB2 by all three BPs. CRE1-3 each bind p300 and mono-methylated H3K4 consistent with enhancer function. In silico analysis suggested that CRE1-3 contain one or more conserved SOX9-binding sites, and this binding was confirmed using chromatin immunoprecipitation on cells derived from mouse embryonic pharyngeal arch. Interphase bacterial artificial chromosome fluorescence in situ hybridization measurements in embryonic craniofacial tissues showed that the orthologous region in mice exhibits Satb2 expression-dependent chromatin decondensation consistent with Satb2 being a target gene of CRE1-3. To assess their in vivo function, we made multiple stable reporter transgenic lines for each enhancer in zebrafish. CRE2 was shown to drive SATB2-like expression in the embryonic craniofacial region. This expression could be eliminated by mutating the SOX9-binding site of CRE2. These observations suggest that SATB2 and SOX9 may be acting together via complex cis-regulation to coordinate the growth of the developing jaw.

PAX6 is a pleiotropic transcription factor essential for the development of several tissues including the eyes, central nervous system, and some endocrine glands. Recently it has also been shown to be important for the maintenance and functioning of corneal and pancreatic tissues in adults. We hypothesized that PAX6 is important for the maintenance of brain integrity in humans, and that adult heterozygotes may have abnormalities of cortical patterning analogous to those found in mouse models.

Cornelia de Lange syndrome (CdLS) is a multisystem disorder with distinctive facial appearance, intellectual disability and growth failure as prominent features. Most individuals with typical CdLS have de novo heterozygous loss-of-function mutations in NIPBL with mosaic individuals representing a significant proportion. Mutations in other cohesin components, SMC1A, SMC3, HDAC8 and RAD21 cause less typical CdLS.

Long noncoding RNAs (lncRNAs) have been implicated in various biological functions including the regulation of gene expression, however, the functionality of lncRNAs is not clearly understood and conflicting conclusions have often been reached when comparing different methods to investigate them. Moreover, little is known about the upstream regulation of lncRNAs. Here we show that the short isoform (p52) of a transcriptional co-activator-PC4 and SF2 interacting protein (Psip1), which is known to be involved in linking transcription to RNA processing, specifically regulates the expression of the lncRNA Hottip-located at the 5' end of the Hoxa locus. Using both knockdown and knockout approaches we show that Hottip expression is required for activation of the 5' Hoxa genes (Hoxa13 and Hoxa10/11) and for retaining Mll1 at the 5' end of Hoxa. Moreover, we demonstrate that artificially inducing Hottip expression is sufficient to activate the 5' Hoxa genes and that Hottip RNA binds to the 5' end of Hoxa. By engineering premature transcription termination, we show that it is the Hottip lncRNA molecule itself, not just Hottip transcription that is required to maintains active expression of posterior Hox genes. Our data show a direct role for a lncRNA molecule in regulating the expression of developmentally-regulated mRNA genes in cis.

Ocular coloboma (OC) is a defect in optic fissure closure and is a common cause of severe congenital visual impairment. Bilateral OC is primarily genetically determined and shows marked locus heterogeneity. Whole-exome sequencing (WES) was used to analyze 12 trios (child affected with OC and both unaffected parents). This identified de novo mutations in 10 different genes in eight probands. Three of these genes encoded proteins associated with actin cytoskeleton dynamics: ACTG1, TWF1, and LCP1. Proband-only WES identified a second unrelated individual with isolated OC carrying the same ACTG1 allele, encoding p.(Pro70Leu). Both individuals have normal neurodevelopment with no extra-ocular signs of Baraitser-Winter syndrome. We found this mutant protein to be incapable of incorporation into F-actin. The LCP1 and TWF1 variants each resulted in only minor disturbance of actin interactions, and no further plausibly causative variants were identified in these genes on resequencing 380 unrelated individuals with OC.

Pax6 is a developmental control gene essential for eye development throughout the animal kingdom. In addition, Pax6 plays key roles in other parts of the CNS, olfactory system, and pancreas. In mammals a single Pax6 gene encoding multiple isoforms delivers these pleiotropic functions. Here we provide evidence that the genomes of many other vertebrate species contain multiple Pax6 loci. We sequenced Pax6-containing BACs from the cartilaginous elephant shark (Callorhinchus milii) and found two distinct Pax6 loci. Pax6.1 is highly similar to mammalian Pax6, while Pax6.2 encodes a paired-less Pax6. Using synteny relationships, we identify homologs of this novel paired-less Pax6.2 gene in lizard and in frog, as well as in zebrafish and in other teleosts. In zebrafish two full-length Pax6 duplicates were known previously, originating from the fish-specific genome duplication (FSGD) and expressed in divergent patterns due to paralog-specific loss of cis-elements. We show that teleosts other than zebrafish also maintain duplicate full-length Pax6 loci, but differences in gene and regulatory domain structure suggest that these Pax6 paralogs originate from a more ancient duplication event and are hence renamed as Pax6.3. Sequence comparisons between mammalian and elephant shark Pax6.1 loci highlight the presence of short- and long-range conserved noncoding elements (CNEs). Functional analysis demonstrates the ancient role of long-range enhancers for Pax6 transcription. We show that the paired-less Pax6.2 ortholog in zebrafish is expressed specifically in the developing retina. Transgenic analysis of elephant shark and zebrafish Pax6.2 CNEs with homology to the mouse NRE/Pα internal promoter revealed highly specific retinal expression. Finally, morpholino depletion of zebrafish Pax6.2 resulted in a "small eye" phenotype, supporting a role in retinal development. In summary, our study reveals that the pleiotropic functions of Pax6 in vertebrates are served by a divergent family of Pax6 genes, forged by ancient duplication events and by independent, lineage-specific gene losses.

Neural stem cell self-renewal, neurogenesis, and cell fate determination are processes that control the generation of specific classes of neurons at the correct place and time. The transcription factor Pax6 is essential for neural stem cell proliferation, multipotency, and neurogenesis in many regions of the central nervous system, including the cerebral cortex. We used Pax6 as an entry point to define the cellular networks controlling neural stem cell self-renewal and neurogenesis in stem cells of the developing mouse cerebral cortex. We identified the genomic binding locations of Pax6 in neocortical stem cells during normal development and ascertained the functional significance of genes that we found to be regulated by Pax6, finding that Pax6 positively and directly regulates cohorts of genes that promote neural stem cell self-renewal, basal progenitor cell genesis, and neurogenesis. Notably, we defined a core network regulating neocortical stem cell decision-making in which Pax6 interacts with three other regulators of neurogenesis, Neurog2, Ascl1, and Hes1. Analyses of the biological function of Pax6 in neural stem cells through phenotypic analyses of Pax6 gain- and loss-of-function mutant cortices demonstrated that the Pax6-regulated networks operating in neural stem cells are highly dosage sensitive. Increasing Pax6 levels drives the system towards neurogenesis and basal progenitor cell genesis by increasing expression of a cohort of basal progenitor cell determinants, including the key transcription factor Eomes/Tbr2, and thus towards neurogenesis at the expense of self-renewal. Removing Pax6 reduces cortical stem cell self-renewal by decreasing expression of key cell cycle regulators, resulting in excess early neurogenesis. We find that the relative levels of Pax6, Hes1, and Neurog2 are key determinants of a dynamic network that controls whether neural stem cells self-renew, generate cortical neurons, or generate basal progenitor cells, a mechanism that has marked parallels with the transcriptional control of embryonic stem cell self-renewal.

Compromised heat shock protein 90 (Hsp90) function reveals cryptic phenotypes in flies and plants. These observations were interpreted to suggest that this molecular stress-response chaperone has a capacity to buffer underlying genetic variation. Conversely, the protective role of Hsp90 could account for the variable penetrance or severity of some heritable developmental malformations in vertebrates. Using zebrafish as a model, we defined Hsp90 inhibitor levels that did not induce a heat shock response or perturb phenotype in wild-type strains. Under these conditions the severity of the recessive eye phenotype in sunrise, caused by a pax6b mutation, was increased, while in dreumes, caused by a sufu mutation, it was decreased. In another strain, a previously unobserved spectrum of severe structural eye malformations, reminiscent of anophthalmia, microphthalmia, and nanophthalmia complex in humans, was uncovered by this limited inhibition of Hsp90 function. Inbreeding of offspring from selected unaffected carrier parents led to significantly elevated malformation frequencies and revealed the oligogenic nature of this phenotype. Unlike in Drosophila, Hsp90 inhibition can decrease developmental stability in zebrafish, as indicated by increased asymmetric presentation of anophthalmia, microphthalmia, and nanophthalmia and sunrise phenotypes. Analysis of the sunrise pax6b mutation suggests a molecular mechanism for the buffering of mutations by Hsp90. The zebrafish studies imply that mild perturbation of Hsp90 function at critical developmental stages may underpin the variable penetrance and expressivity of many developmental anomalies where the interaction between genotype and environment plays a major role.

Colobomata represent visually impairing ocular closure defects that are associated with a diverse range of developmental anomalies. Characterization of a chromosome 8q21.2-q22.1 segmental deletion in a patient with chorioretinal coloboma revealed elements of nonallelic homologous recombination and nonhomologous end joining. This genomic architecture extends the range of chromosomal rearrangements associated with human disease and indicates that a broader spectrum of human chromosomal rearrangements may use coupled homologous and nonhomologous mechanisms. We also demonstrate that the segmental deletion encompasses GDF6, encoding a member of the bone-morphogenetic protein family, and that inhibition of gdf6a in a model organism accurately recapitulates the proband's phenotype. The spectrum of disorders generated by morpholino inhibition and the more severe defects (microphthalmia and anophthalmia) observed at higher doses illustrate the key role of GDF6 in ocular development. These results underscore the value of integrated clinical and molecular investigation of patients with chromosomal anomalies.

Cilia assembly and disassembly are coupled to actin dynamics, ensuring a coherent cellular response during environmental change. How these processes are integrated remains undefined. The histone lysine demethylase KDM3A plays important roles in organismal homeostasis. Loss-of-function mouse models of Kdm3a phenocopy features associated with human ciliopathies, whereas human somatic mutations correlate with poor cancer prognosis. We demonstrate that absence of KDM3A facilitates ciliogenesis, but these resulting cilia have an abnormally wide range of axonemal lengths, delaying disassembly and accumulating intraflagellar transport (IFT) proteins. KDM3A plays a dual role by regulating actin gene expression and binding to the actin cytoskeleton, creating a responsive "actin gate" that involves ARP2/3 activity and IFT. Promoting actin filament formation rescues KDM3A mutant ciliary defects. Conversely, the simultaneous depolymerization of actin networks and IFT overexpression mimics the abnormal ciliary traits of KDM3A mutants. KDM3A is thus a negative regulator of ciliogenesis required for the controlled recruitment of IFT proteins into cilia through the modulation of actin dynamics.

Seventy-seven patients with aniridia, referred for cytogenetic analysis predominantly to assess Wilms tumor risk, were studied by fluorescence in situ hybridization (FISH), through use of a panel of cosmids encompassing the aniridia-associated PAX6 gene, the Wilms tumor predisposition gene WT1, and flanking markers, in distal chromosome 11p13. Thirty patients were found to be chromosomally abnormal. Cytogenetically visible interstitial deletions involving 11p13 were found in 13 patients, 11 of which included WT1. A further 13 patients had cryptic deletions detectable only by FISH, 3 of which included WT1. Six of these, with deletions <500 kb, share a similar proximal breakpoint within a cosmid containing the last 10 exons of PAX6 and part of the neighboring gene, ELP4. Two of these six patients were mosaic for the deletion. The remaining four had chromosomal rearrangements: an unbalanced translocation, t(11;13), with a deletion including the WAGR (Wilms' tumor, aniridia, genitourinary abnormalities, and mental retardation) region, and three balanced rearrangements with what appear to be position effect breakpoints 3' of PAX6: (a) a t(7;11) with the 11p13 breakpoint approximately 30 kb downstream of PAX6, (b) a dir ins(12;11) with a breakpoint >50 kb from PAX6, and (c) an inv(11)(p13q13) with a breakpoint >75 kb downstream of PAX6. The proportion and spectrum of chromosome anomalies in familial (4/14, or 28.5%) and sporadic (26/63, or 41%) cases are not significantly different. An unexpectedly high frequency of chromosomal rearrangements is associated with both sporadic and familial aniridia in this cohort.

Mutations in the coding sequence of SOX9 cause campomelic dysplasia (CD), a disorder of skeletal development associated with 46,XY disorders of sex development (DSDs). Translocations, deletions, and duplications within a ∼2 Mb region upstream of SOX9 can recapitulate the CD-DSD phenotype fully or partially, suggesting the existence of an unusually large cis-regulatory control region. Pierre Robin sequence (PRS) is a craniofacial disorder that is frequently an endophenotype of CD and a locus for isolated PRS at ∼1.2-1.5 Mb upstream of SOX9 has been previously reported. The craniofacial regulatory potential within this locus, and within the greater genomic domain surrounding SOX9, remains poorly defined. We report two novel deletions upstream of SOX9 in families with PRS, allowing refinement of the regions harboring candidate craniofacial regulatory elements. In parallel, ChIP-Seq for p300 binding sites in mouse craniofacial tissue led to the identification of several novel craniofacial enhancers at the SOX9 locus, which were validated in transgenic reporter mice and zebrafish. Notably, some of the functionally validated elements fall within the PRS deletions. These studies suggest that multiple noncoding elements contribute to the craniofacial regulation of SOX9 expression, and that their disruption results in PRS.

Exome sequence analysis of affected individuals from two families with autosomal-dominant inheritance of coloboma identified two different cosegregating heterozygous nonsense mutations (c.370C>T [p.Arg124*] and c. 1066G>T [p.Glu356*]) in YAP1. The phenotypes of the affected families differed in that one included no extraocular features and the other manifested with highly variable multisystem involvement, including hearing loss, intellectual disability, hematuria, and orofacial clefting. A combined LOD score of 4.2 was obtained for the association between YAP1 loss-of-function mutations and the phenotype in these families. YAP1 encodes an effector of the HIPPO-pathway-induced growth response, and whole-mount in situ hybridization in mouse embryos has shown that Yap1 is strongly expressed in the eye, brain, and fusing facial processes. RT-PCR showed that an alternative transcription start site (TSS) in intron 1 of YAP1 and Yap1 is widely used in human and mouse development, respectively. Transcripts from the alternative TSS are predicted to initiate at codon Met179 relative to the canonical transcript (RefSeq NM_001130145). In these alternative transcripts, the c.370C>T mutation in family 1305 is within the 5' UTR and cannot result in nonsense-mediated decay (NMD). The c. 1066G>T mutation in family 132 should result in NMD in transcripts from either TSS. Amelioration of the phenotype by the alternative transcripts provides a plausible explanation for the phenotypic differences between the families.

The lysine demethylase Kdm3a (Jhdm2a, Jmjd1a) is required for male fertility, sex determination, and metabolic homeostasis through its nuclear role in chromatin remodeling. Many histone-modifying enzymes have additional nonhistone substrates, as well as nonenzymatic functions, contributing to the full spectrum of events underlying their biological roles. We present two Kdm3a mouse models that exhibit cytoplasmic defects that may account in part for the globozoospermia phenotype reported previously. Electron microscopy revealed abnormal acrosome and manchette and the absence of implantation fossa at the caudal end of the nucleus in mice without Kdm3a demethylase activity, which affected cytoplasmic structures required to elongate the sperm head. We describe an enzymatically active new Kdm3a isoform and show that subcellular distribution, protein levels, and lysine demethylation activity of Kdm3a depended on Hsp90. We show that Kdm3a localizes to cytoplasmic structures of maturing spermatids affected in Kdm3a mutant mice, which in turn display altered fractionation of β-actin and γ-tubulin. Kdm3a is therefore a multifunctional Hsp90 client protein that participates directly in the regulation of cytoskeletal components.

Epigenetic cellular memory mechanisms that involve polycomb and trithorax group of proteins are well conserved across metazoans. The cis-acting elements interacting with these proteins, however, are poorly understood in mammals. In a directed search we identified a potential polycomb responsive element with 25 repeats of YY1 binding motifthatwe designate PRE-PIK3C2B as it occurs in the first intron of human PIK3C2B gene. It down regulates reporter gene expression in HEK cells and the repression is dependent on polycomb group of proteins (PcG). We demonstrate that PRE-PIK3C2B interacts directly with YY1 in vitro and recruits PRC2 complex in vivo. The localization of PcG proteins including YY1 to PRE-PIK3C2B in HEK cells is decreased on knock-down of either YY1 or SUZ12. Endogenous PRE-PIK3C2B shows bivalent marking having H3K27me3 and H3K4me3 for repressed and active state respectively. In transgenic Drosophila, PRE-PIK3C2B down regulates mini-white expression, exhibits variegation and pairing sensitive silencing (PSS), which has not been previously demonstrated for mammalian PRE. Taken together, our results strongly suggest that PRE-PIK3C2B functions as a site of interaction for polycomb proteins.

Major malformations of the human eye, including microphthalmia and anophthalmia, are examples of phenotypes that recur in families yet often show no clear Mendelian inheritance pattern. Defining loci by mapping is therefore rarely feasible. Using a candidate-gene approach, we have identified heterozygous coding-region changes in the homeobox gene OTX2 in eight families with ocular malformations. The expression pattern of OTX2 in human embryos is consistent with the eye phenotypes observed in the patients, which range from bilateral anophthalmia to retinal defects resembling Leber congenital amaurosis and pigmentary retinopathy. Magnetic resonance imaging scans revealed defects of the optic nerve, optic chiasm, and, in some cases, brain. In two families, the mutations appear to have occurred de novo in severely affected offspring, and, in two other families, the mutations have been inherited from a gonosomal mosaic parent. Data from these four families support a simple model in which OTX2 heterozygous loss-of-function mutations cause ocular malformations. Four additional families display complex inheritance patterns, suggesting that OTX2 mutations alone may not lead to consistent phenotypes. The high incidence of mosaicism and the reduced penetrance have implications for genetic counseling.

Gene duplication is a major driver of evolutionary divergence. In most vertebrates a single PAX6 gene encodes a transcription factor required for eye, brain, olfactory system, and pancreas development. In zebrafish, following a postulated whole-genome duplication event in an ancestral teleost, duplicates pax6a and pax6b jointly fulfill these roles. Mapping of the homozygously viable eye mutant sunrise identified a homeodomain missense change in pax6b, leading to loss of target binding. The mild phenotype emphasizes role-sharing between the co-orthologues. Meticulous mapping of isolated BACs identified perturbed synteny relationships around the duplicates. This highlights the functional conservation of pax6 downstream (3') control sequences, which in most vertebrates reside within the introns of a ubiquitously expressed neighbour gene, ELP4, whose pax6a-linked exons have been lost in zebrafish. Reporter transgenic studies in both mouse and zebrafish, combined with analysis of vertebrate sequence conservation, reveal loss and retention of specific cis-regulatory elements, correlating strongly with the diverged expression of co-orthologues, and providing clear evidence for evolution by subfunctionalization.