Characterization of an eye field-like state during optic vesicle organoid development.

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. This crucial 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 time-course transcriptomic data allowing us to identify dynamic gene expression programmes that characterize this cellular-state transition. Integrating this with chromatin accessibility data suggests a direct role of canonical EF transcription factors in regulating these gene expression changes, and highlights candidate cis-regulatory elements through which these transcription factors 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.

A sensitive and affordable multiplex RT-qPCR assay for SARS-CoV-2 detection.

With the ongoing COVID-19 (Coronavirus Disease 2019) pandemic, caused by the novel coronavirus SARS-CoV-2 (Severe Acute Respiratory Syndrome Coronavirus 2), there is a need for sensitive, specific, and affordable diagnostic tests to identify infected individuals, not all of whom are symptomatic. The most sensitive test involves the detection of viral RNA using RT-qPCR (quantitative reverse transcription PCR), with many commercial kits now available for this purpose. However, these are expensive, and supply of such kits in sufficient numbers cannot always be guaranteed. We therefore developed a multiplex assay using well-established SARS-CoV-2 targets alongside a human cellular control (RPP30) and a viral spike-in control (Phocine Herpes Virus 1 [PhHV-1]), which monitor sample quality and nucleic acid extraction efficiency, respectively. Here, we establish that this test performs as well as widely used commercial assays, but at substantially reduced cost. Furthermore, we demonstrate >1,000-fold variability in material routinely collected by combined nose and throat swabbing and establish a statistically significant correlation between the detected level of human and SARS-CoV-2 nucleic acids. The inclusion of the human control probe in our assay therefore provides a quantitative measure of sample quality that could help reduce false-negative rates. We demonstrate the feasibility of establishing a robust RT-qPCR assay at approximately 10% of the cost of equivalent commercial assays, which could benefit low-resource environments and make high-volume testing affordable.

cGAS surveillance of micronuclei links genome instability to innate immunity.

DNA is strictly compartmentalized within the nucleus to prevent autoimmunity; despite this, cyclic GMP-AMP synthase (cGAS), a cytosolic sensor of double-stranded DNA, is activated in autoinflammatory disorders and by DNA damage. Precisely how cellular DNA gains access to the cytoplasm remains to be determined. Here, we report that cGAS localizes to micronuclei arising from genome instability in a mouse model of monogenic autoinflammation, after exogenous DNA damage and spontaneously in human cancer cells. Such micronuclei occur after mis-segregation of DNA during cell division and consist of chromatin surrounded by its own nuclear membrane. Breakdown of the micronuclear envelope, a process associated with chromothripsis, leads to rapid accumulation of cGAS, providing a mechanism by which self-DNA becomes exposed to the cytosol. cGAS is activated by chromatin, and consistent with a mitotic origin, micronuclei formation and the proinflammatory response following DNA damage are cell-cycle dependent. By combining live-cell laser microdissection with single cell transcriptomics, we establish that interferon-stimulated gene expression is induced in micronucleated cells. We therefore conclude that micronuclei represent an important source of immunostimulatory DNA. As micronuclei formed from lagging chromosomes also activate this pathway, recognition of micronuclei by cGAS may act as a cell-intrinsic immune surveillance mechanism that detects a range of neoplasia-inducing processes.

A Restricted Repertoire of De Novo Mutations in ITPR1 Cause Gillespie Syndrome with Evidence for Dominant-Negative Effect.

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.

Heterozygous loss-of-function mutations in YAP1 cause both isolated and syndromic optic fissure closure defects.

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.

Monoallelic and biallelic mutations in MAB21L2 cause a spectrum of major eye malformations.

We identified four different missense mutations in the single-exon gene MAB21L2 in eight individuals with bilateral eye malformations from five unrelated families via three independent exome sequencing projects. Three mutational events altered the same amino acid (Arg51), and two were identical de novo mutations (c.151C>T [p.Arg51Cys]) in unrelated children with bilateral anophthalmia, intellectual disability, and rhizomelic skeletal dysplasia. c.152G>A (p.Arg51His) segregated with autosomal-dominant bilateral colobomatous microphthalmia in a large multiplex family. The fourth heterozygous mutation (c.145G>A [p.Glu49Lys]) affected an amino acid within two residues of Arg51 in an adult male with bilateral colobomata. In a fifth family, a homozygous mutation (c.740G>A [p.Arg247Gln]) altering a different region of the protein was identified in two male siblings with bilateral retinal colobomata. In mouse embryos, Mab21l2 showed strong expression in the developing eye, pharyngeal arches, and limb bud. As predicted by structural homology, wild-type MAB21L2 bound single-stranded RNA, whereas this activity was lost in all altered forms of the protein. MAB21L2 had no detectable nucleotidyltransferase activity in vitro, and its function remains unknown. Induced expression of wild-type MAB21L2 in human embryonic kidney 293 cells increased phospho-ERK (pERK1/2) signaling. Compared to the wild-type and p.Arg247Gln proteins, the proteins with the Glu49 and Arg51 variants had increased stability. Abnormal persistence of pERK1/2 signaling in MAB21L2-expressing cells during development is a plausible pathogenic mechanism for the heterozygous mutations. The phenotype associated with the homozygous mutation might be a consequence of complete loss of MAB21L2 RNA binding, although the cellular function of this interaction remains unknown.

FRA2A is a CGG repeat expansion associated with silencing of AFF3.

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.

Disruption of SATB2 or its long-range cis-regulation by SOX9 causes a syndromic form of Pierre Robin sequence.

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.

A trans-acting protein effect causes severe eye malformation in the Mp mouse.

Mp is an irradiation-induced mouse mutation associated with microphthalmia, micropinna and hind limb syndactyly. We show that Mp is caused by a 660 kb balanced inversion on chromosome 18 producing reciprocal 3-prime gene fusion events involving Fbn2 and Isoc1. The Isoc1-Fbn2 fusion gene (Isoc1(Mp)) mRNA has a frameshift and early stop codon resulting in nonsense mediated decay. Homozygous deletions of Isoc1 do not support a significant developmental role for this gene. The Fbn2-Isoc1 fusion gene (Fbn2 (Mp)) predicted protein consists of the N-terminal Fibrillin-2 (amino acids 1-2646, exons 1-62) lacking the C-terminal furin-cleavage site with a short out-of-frame extension encoded by the final exon of Isoc1. The Mp limb phenotype is consistent with that reported in Fbn2 null embryos. However, severe eye malformations, a defining feature of Mp, are not seen in Fbn2 null animals. Fibrillin-2(Mp) forms large fibrillar structures within the rough endoplasmic reticulum (rER) associated with an unfolded protein response and quantitative mass spectrometry shows a generalised defect in protein secretion in conditioned media from mutant cells. In the embryonic eye Fbn2 is expressed within the peripheral ciliary margin (CM). Mp embryos show reduced canonical Wnt-signalling in the CM - known to be essential for ciliary body development - and show subsequent aplasia of CM-derived structures. We propose that the Mp "worse-than-null" eye phenotype plausibly results from a failure in normal trafficking of proteins that are co-expressed with Fbn2 within the CM. The prediction of similar trans-acting protein effects will be an important challenge in the medical interpretation of human mutations from whole exome sequencing.

A CGG-repeat expansion mutation in ZNF713 causes FRA7A: association with autistic spectrum disorder in two families.

We report de novo occurrence of the 7p11.2 folate-sensitive fragile site FRA7A in a male with an autistic spectrum disorder (ASD) due to a CGG-repeat expansion mutation (∼450 repeats) in a 5' intron of ZNF713. This expanded allele showed hypermethylation of the adjacent CpG island with reduced ZNF713 expression observed in a proband-derived lymphoblastoid cell line (LCL). His unaffected mother carried an unmethylated premutation (85 repeats). This CGG-repeat showed length polymorphism in control samples (five to 22 repeats). In a second unrelated family, three siblings with ASD and their unaffected father were found to carry FRA7A premutations, which were partially or mosaically methylated. In one of the affected siblings, mitotic instability of the premutation was observed. ZNF713 expression in LCLs in this family was increased in three of these four premutation carriers. A firm link cannot yet be established between ASD and the repeat expansion mutation but plausible pathogenic mechanisms are discussed.

Identification and functional modelling of plausibly causative cis-regulatory variants in a highly-selected cohort with X-linked intellectual disability.

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.

SMCHD1 mutations associated with a rare muscular dystrophy can also cause isolated arhinia and Bosma arhinia microphthalmia syndrome.

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.

Genetic Analysis of 'PAX6-Negative' Individuals with Aniridia or Gillespie Syndrome.

We report molecular genetic analysis of 42 affected individuals referred with a diagnosis of aniridia who previously screened as negative for intragenic PAX6 mutations. Of these 42, the diagnoses were 31 individuals with aniridia and 11 individuals referred with a diagnosis of Gillespie syndrome (iris hypoplasia, ataxia and mild to moderate developmental delay). Array-based comparative genomic hybridization identified six whole gene deletions: four encompassing PAX6 and two encompassing FOXC1. Six deletions with plausible cis-regulatory effects were identified: five that were 3' (telomeric) to PAX6 and one within a gene desert 5' (telomeric) to PITX2. Sequence analysis of the FOXC1 and PITX2 coding regions identified two plausibly pathogenic de novo FOXC1 missense mutations (p.Pro79Thr and p.Leu101Pro). No intragenic mutations were detected in PITX2. FISH mapping in an individual with Gillespie-like syndrome with an apparently balanced X;11 reciprocal translocation revealed disruption of a gene at each breakpoint: ARHGAP6 on the X chromosome and PHF21A on chromosome 11. In the other individuals with Gillespie syndrome no mutations were identified in either of these genes, or in HCCS which lies close to the Xp breakpoint. Disruption of PHF21A has previously been implicated in the causation of intellectual disability (but not aniridia). Plausibly causative mutations were identified in 15 out of 42 individuals (12/32 aniridia; 3/11 Gillespie syndrome). Fourteen of these mutations presented in the known aniridia genes; PAX6, FOXC1 and PITX2. The large number of individuals in the cohort with no mutation identified suggests greater locus heterogeneity may exist in both isolated and syndromic aniridia than was previously appreciated.

A syndromic form of Pierre Robin sequence is caused by 5q23 deletions encompassing FBN2 and PHAX.

Pierre Robin sequence (PRS) is an aetiologically distinct subgroup of cleft palate. We aimed to define the critical genomic interval from five different 5q22-5q31 deletions associated with PRS or PRS-associated features and assess each gene within the region as a candidate for the PRS component of the phenotype. Clinical array-based comparative genome hybridisation (aCGH) data were used to define a 2.08 Mb minimum region of overlap among four de novo deletions and one mother-son inherited deletion associated with at least one component of PRS. Commonly associated anomalies were talipes equinovarus (TEV), finger contractures and crumpled ear helices. Expression analysis of the orthologous genes within the PRS critical region in embryonic mice showed that the strongest candidate genes were FBN2 and PHAX. Targeted aCGH of the critical region and sequencing of these genes in a cohort of 25 PRS patients revealed no plausible disease-causing mutations. In conclusion, deletion of ∼2 Mb on 5q23 region causes a clinically recognisable subtype of PRS. Haploinsufficiency for FBN2 accounts for the digital and auricular features. A possible critical region for TEV is distinct and telomeric to the PRS region. The molecular basis of PRS in these cases remains undetermined but haploinsufficiency for PHAX is a plausible mechanism.

Clinical and mutation analysis of 51 probands with anophthalmia and/or severe microphthalmia from a single center.

Clinical evaluation and mutation analysis was performed in 51 consecutive probands with severe eye malformations - anophthalmia and/or severe microphthalmia - seen in a single specialist ophthalmology center. The mutation analysis consisted of bidirectional sequencing of the coding regions of SOX2, OTX2, PAX6 (paired domain), STRA6, BMP4, SMOC1, FOXE3, and RAX, and genome-wide array-based copy number assessment. Fifteen (29.4%) of the 51 probands had likely causative mutations affecting SOX2 (9/51), OTX2 (5/51), and STRA6 (1/51). Of the cases with bilateral anophthalmia, 9/12 (75%) were found to be mutation positive. Three of these mutations were large genomic deletions encompassing SOX2 (one case) or OTX2 (two cases). Familial inheritance of three intragenic, plausibly pathogenic, and heterozygous mutations was observed. An unaffected carrier parent of an affected child with an identified OTX2 mutation confirmed the previously reported nonpenetrance for this disorder. Two families with SOX2 mutations demonstrated a parent and child both with significant but highly variable eye malformations. Heterozygous loss-of-function mutations in SOX2 and OTX2 are the most common genetic pathology associated with severe eye malformations and bi-allelic loss-of-function in STRA6 is confirmed as an emerging cause of nonsyndromal eye malformations.