Search

Anophthalmia (missing eye) describes a failure of early embryonic ocular development. Mutations in a relatively small set of genes account for 75% of bilateral anophthalmia cases, yet 25% of families currently are left without a molecular diagnosis. Here, we report our experimental work that aimed to uncover the developmental and genetic basis of the anophthalmia characterising the X-linked Ie (eye-ear reduction) X-ray-induced allele in mouse that was first identified in 1947. Histological analysis of the embryonic phenotype showed failure of normal eye development after the optic vesicle stage with particularly severe malformation of the ventral retina. Linkage analysis mapped this mutation to a ~6 Mb region on the X chromosome. Short- and long-read whole-genome sequencing (WGS) of affected and unaffected male littermates confirmed the Ie linkage but identified no plausible causative variants or structural rearrangements. These analyses did reduce the critical candidate interval and revealed evidence of multiple variants within the ancestral DNA, although none were found that altered coding sequences or that were unique to Ie. To investigate early embryonic events at a genetic level, we then generated mouse ES cells derived from male Ie embryos and wild type littermates. RNA-seq and accessible chromatin sequencing (ATAC-seq) data generated from cultured optic vesicle organoids did not reveal any large differences in gene expression or accessibility of putative cis-regulatory elements between Ie and wild type. However, an unbiased TF-footprinting analysis of accessible chromatin regions did provide evidence of a genome-wide reduction in binding of transcription factors associated with ventral eye development in Ie, and evidence of an increase in binding of the Zic-family of transcription factors, including Zic3, which is located within the Ie-refined critical interval. We conclude that the refined Ie critical region at chrX: 56,145,000–58,385,000 contains multiple genetic variants that may be linked to altered cis regulation but does not contain a convincing causative mutation. Changes in the binding of key transcription factors to chromatin causing altered gene expression during development, possibly through a subtle mis-regulation of Zic3, presents a plausible cause for the anophthalmia phenotype observed in Ie, but further work is required to determine the precise causative allele and its genetic mechanism.

Specification of the eye field (EF) within the neural plate marks the earliest detectable stage of eye development. Experimental evidence, primarily from non-mammalian model systems, indicates that the stable formation of this group of cells requires the activation of a set of key transcription factors (TFs). This critical event is challenging to probe in mammals and, quantitatively, little is known regarding the regulation of the transition of cells to this ocular fate. Using optic vesicle organoids to model the onset of the EF, we generate timecourse transcriptomic data allowing us to identify dynamic gene-expression programs that characterise this cellular-state transition. Integrating this with chromatin accessibility data suggests a direct role of canonical EFTFs in regulating these gene-expression changes, and high-lights candidate cis-regulatory elements through which these TFs act. Finally, we begin to test a subset of these candidate enhancer elements, within the organoid system, by perturbing the underlying DNA sequence and measuring transcriptomic changes during EF activation.

Anophthalmia (missing eye) describes a failure of early embryonic ocular development. Mutations in a relatively small set of genes account for 75% of bilateral anophthalmia cases, yet 25% of families currently are left without a molecular diagnosis. Here we report our experimental work that aimed to uncover the developmental and genetic basis of the anophthalmia characterising the X-linked Ie (eye-ear reduction) X-ray induced allele in mouse that was first identified in 1947. Histological analysis of the embryonic phenotype showed failure of normal eye development after the optic vesicle stage with particularly severe malformation of the ventral retina. Linkage analysis mapped this mutation to a ∼ 6Mb region on the X chromosome. Short and long read whole-genome sequencing (WGS) of affected and unaffected male littermates confirmed the Ie linkage but identified no plausible causative variants or structural rearrangements. These analyses did reduce the critical candidate interval and revealed evidence of multiple variants within the ancestral DNA, although none were found that altered coding sequences or that were unique to Ie. To investigate early embryonic events at a genetic level, we then generated mouse ES cells derived from male Ie embryos and wild type littermates. RNA-seq and accessible chromatin sequencing (ATAC-seq) data generated from cultured optic vesicle organoids did not reveal any large differences in gene expression or accessibility of putative cis-regulatory elements between Ie and wild type. However, an unbiased TF-footprinting analysis of accessible chromatin regions did provide evidence of a genome-wide reduction in binding of transcription factors associated with ventral eye development in Ie, and evidence of an increase in binding of the Zic-family of transcription factors, including Zic3, which is located within the Ie-refined critical interval. We conclude that the refined Ie critical region at chrX: 56,145,000-58,385,000 contains multiple genetic variants that may be linked to altered cis regulation but does not contain a convincing causative mutation. Changes in the binding of key transcription factors to chromatin causing altered gene expression during development, possibly through a subtle mis-regulation of Zic3, presents a plausible cause for the anophthalmia phenotype observed in Ie, but further work is required to determine the precise causative allele and its genetic mechanism.

Funder: BBSRC studentship, Identifying causative variants in cis-regulatory elements (CRE) in neurodevelopmental disorders has proven challenging. We have used in vivo functional analyses to categorize rigorously filtered CRE variants in a clinical cohort that is plausibly enriched for causative CRE mutations: 48 unrelated males with a family history consistent with X-linked intellectual disability (XLID) in whom no detectable cause could be identified in the coding regions of the X chromosome (chrX). Targeted sequencing of all chrX CRE identified six rare variants in five affected individuals that altered conserved bases in CRE targeting known XLID genes and segregated appropriately in families. Two of these variants, FMR1CRE and TENM1CRE, showed consistent site- and stage-specific differences of enhancer function in the developing zebrafish brain using dual-color fluorescent reporter assay. Mouse models were created for both variants. In male mice Fmr1CRE induced alterations in neurodevelopmental Fmr1 expression, olfactory behavior and neurophysiological indicators of FMRP function. The absence of another likely causative variant on whole genome sequencing further supported FMR1CRE as the likely basis of the XLID in this family. Tenm1CRE mice showed no phenotypic anomalies. Following the release of gnomAD 2.1, reanalysis showed that TENM1CRE exceeded the maximum plausible population frequency of a XLID causative allele. Assigning causative status to any ultra-rare CRE variant remains problematic and requires disease-relevant in vivo functional data from multiple sources. The sequential and bespoke nature of such analyses renders them time-consuming and challenging to scale for routine clinical use.

Undiagnosed neurodevelopmental disease is significantly associated with rare variants in cis-regulatory elements (CRE) but demonstrating causality is challenging as target gene consequences may differ from a causative variant affecting the coding region. Here, we address this challenge by applying a procedure to discriminate likely diagnostic regulatory variants from those of neutral/low-penetrant effect. We identified six rare CRE variants using targeted and whole genome sequencing in 48 unrelated males with apparent X-linked intellectual disability (XLID) but without detectable coding region variants. These variants segregated appropriately in families and altered conserved bases in predicted CRE targeting known XLID genes. Three were unique and three were rare but too common to be plausibly causative for XLID. We compared the cis-regulatory activity of wild-type and mutant alleles in zebrafish embryos using dual-color fluorescent reporters. Two variants showed striking changes: one plausibly causative (FMR1CRE) and the other likely neutral/low-penetrant (TENM1CRE). These variants were “knocked-in” to mice and both altered embryonic neural expression of their target gene. Only Fmr1CRE mice showed disease-relevant behavioral defects. FMR1CRE is plausibly disease-associated resulting in complex misregulation of Fmr1/FMRP rather than loss-of-function. This is consistent both with absence of Fragile X syndrome in the probands and the observed electrophysiological anomalies in the FMR1CRE mouse brain. Although disruption of in vivo patterns of endogenous gene expression in disease-relevant tissues by CRE variants cannot be used as strong evidence for Mendelian disease association, in conjunction with extreme rarity in human populations and with relevant knock-in mouse phenotypes, such variants can become likely pathogenic.

Arhinia, or absence of the nose, is a rare malformation of unknown etiology that is often accompanied by ocular and reproductive defects. Sequencing of 40 people with arhinia revealed that 84% of probands harbor a missense mutation localized to a constrained region of SMCHD1 encompassing the ATPase domain. SMCHD1 mutations cause facioscapulohumeral muscular dystrophy type 2 (FSHD2) via a trans-acting loss-of-function epigenetic mechanism. We discovered shared mutations and comparable DNA hypomethylation patterning between these distinct disorders. CRISPR/Cas9-mediated alteration of smchd1 in zebrafish yielded arhinia-relevant phenotypes. Transcriptome and protein analyses in arhinia probands and controls showed no differences in SMCHD1 mRNA or protein abundance but revealed regulatory changes in genes and pathways associated with craniofacial patterning. Mutations in SMCHD1 thus contribute to distinct phenotypic spectra, from craniofacial malformation and reproductive disorders to muscular dystrophy, which we speculate to be consistent with oligogenic mechanisms resulting in pleiotropic outcomes.

In the version of this article initially published, Wendy Bickmore and Madapura Pradeepa were incorrectly not indicated as corresponding authors. The error has been corrected in the HTML and PDF versions of the paper.

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.

The strictly regulated expression of most pleiotropic developmental control genes is critically dependent on the activity of long-range cis-regulatory elements. This was revealed by the identification of individuals with a genetic condition lacking coding-region mutations in the gene commonly associated with the disease but having a variety of nearby chromosomal abnormalities, collectively described as cis-ruption disease cases. The congenital eye malformation aniridia is caused by haploinsufficiency of the developmental regulator PAX6. We discovered a de novo point mutation in an ultraconserved cis-element located 150 kb downstream from PAX6 in an affected individual with intact coding region and chromosomal locus. The element SIMO acts as a strong enhancer in developing ocular structures. The mutation disrupts an autoregulatory PAX6 binding site, causing loss of enhancer activity, resulting in defective maintenance of PAX6 expression. These findings reveal a distinct regulatory mechanism for genetic disease by disruption of an autoregulatory feedback loop critical for maintenance of gene expression through development.

Purpose: To characterize features associated with de novo mutations affecting SATB2 function in individuals ascertained on the basis of intellectual disability. Methods: Twenty previously unreported individuals with 19 different SATB2 mutations (11 loss-of-function and 8 missense variants) were studied. Fibroblasts were used to measure mutant protein production. Subcellular localization and mobility of wild-type and mutant SATB2 were assessed using fluorescently tagged protein. Results: Recurrent clinical features included neurodevelopmental impairment (19/19), absent/near absent speech (16/19), normal somatic growth (17/19), cleft palate (9/19), drooling (12/19), and dental anomalies (8/19). Six of eight missense variants clustered in the first CUT domain. Sibling recurrence due to gonadal mosaicism was seen in one family. A nonsense mutation in the last exon resulted in production of a truncated protein retaining all three DNA-binding domains. SATB2 nuclear mobility was mutation-dependent; p.Arg389Cys in CUT1 increased mobility and both p.Gly515Ser in CUT2 and p.Gln566Lys between CUT2 and HOX reduced mobility. The clinical features in individuals with missense variants were indistinguishable from those with loss of function. Conclusion: SATB2 haploinsufficiency is a common cause of syndromic intellectual disability. When mutant SATB2 protein is produced, the protein appears functionally inactive with a disrupted pattern of chromatin or matrix association. Genet Med advance online publication 02 February 2017

Folate-sensitive fragile sites (FSFS) are a rare cytogenetically visible subset of dynamic mutations. Of the eight molecularly characterized FSFS, four are associated with intellectual disability (ID). Cytogenetic expression results from CGG tri-nucleotide-repeat expansion mutation associated with local CpG hypermethylation and transcriptional silencing. The best studied is the FRAXA site in the FMR1 gene, where large expansions cause fragile X syndrome, the most common inherited ID syndrome. Here we studied three families with FRA2A expression at 2q11 associated with a wide spectrum of neurodevelopmental phenotypes. We identified a polymorphic CGG repeat in a conserved, brain-active alternative promoter of the AFF3 gene, an autosomal homolog of the X-linked AFF2/FMR2 gene: Expansion of the AFF2 CGG repeat causes FRAXE ID. We found that FRA2A-expressing individuals have mosaic expansions of the AFF3 CGG repeat in the range of several hundred repeat units. Moreover, bisulfite sequencing and pyrosequencing both suggest AFF3 promoter hypermethylation. cSNP-analysis demonstrates monoallelic expression of the AFF3 gene in FRA2A carriers thus predicting that FRA2A expression results in functional haploinsufficiency for AFF3 at least in a subset of tissues. By whole-mount in situ hybridization the mouse AFF3 ortholog shows strong regional expression in the developing brain, somites and limb buds in 9.5–12.5dpc mouse embryos. Our data suggest that there may be an association between FRA2A and a delay in the acquisition of motor and language skills in the families studied here. However, additional cases are required to firmly establish a causal relationship., Author Summary Some human genetic diseases are caused by dynamic mutations, or expansions of a short repeated sequence in the genome that can be unstably passed on from generation to generation. A subset of these dynamic mutations known as fragile sites can be seen as a break or gap on the chromosome when cells are cultured under specific conditions. To date eight folate-sensitive fragile sites (FSFS) have been characterized, and all are due to CGG-repeat expansions within the 5′ UTR or promoter region of the respective gene. When the repeat expands in size, it becomes hypermethylated and the adjacent gene or genes are transcriptionally silenced. For at least four of the eight known fragile sites this silencing of the associated gene(s) lead to intellectual disability syndromes such as fragile X. In this work we describe molecular characterization of an autosomal FSFS called FRA2A on chromosome 2. As the molecular cause of FRA2A, we identify an expansion of a CGG repeat which subsequently results in silencing of the neighbouring gene AFF3. This gene is one of the four autosomal paralogss of the AFF2/FMR2 gene which, when mutated, is the cause of the FRAXE syndrome. We find that FRA2A expression is associated with highly variable developmental anomalies in the three FRA2A families studied.

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.

Biallelic mutations in the gene encoding DHOdehase [dihydroorotate dehydrogenase (DHODH)], an enzyme required for de novo pyrimidine biosynthesis, have been identified as the cause of Miller (Genee-Weidemann or postaxial acrofacial dysostosis) syndrome (MIM 263750). We report compound heterozygous DHODH mutations in four additional families with typical Miller syndrome. Complementation in auxotrophic yeast demonstrated reduced pyrimidine synthesis and in vitro enzymatic analysis confirmed reduced DHOdehase activity in 11 disease-associated missense mutations, with 7 alleles showing discrepant activity between the assays. These discrepancies are partly explained by the domain structure of DHODH and suggest both assays are useful for interpretation of individual alleles. However, in all affected individuals, the genotype predicts that there should be significant residual DHOdehase activity. Urine samples obtained from two mutation-positive cases showed elevated levels of orotic acid (OA) but not dihydroorotate (DHO), an unexpected finding since these represent the product and the substrate of DHODH enzymatic activity, respectively. Screening of four unrelated cases with overlapping but atypical clinical features showed no mutations in either DHODH or the other de novo pyrimidine biosynthesis genes (CAD, UMPS), with these cases also showing normal levels of urinary OA and DHO. In situ analysis of mouse embryos showed Dhodh, Cad and Umps to be strongly expressed in the pharyngeal arch and limb bud, supporting a site- and stage-specific requirement for de novo pyrimidine synthesis. The developmental sensitivity to reduced pyrimidine synthesis capacity may reflect the requirement for an exceptional mitogenic response to growth factor signalling in the affected tissues.