Stochasticity in genetics and gene regulation.

Development from fertilized egg to functioning multi-cellular organism requires precision. There is no precision, and often no survival, without plasticity. Plasticity is conferred partly by stochastic variation, present inherently in all biological systems. Gene expression levels fluctuate ubiquitously through transcription, alternative splicing, translation and turnover. Small differences in gene expression are exploited to trigger early differentiation, conferring distinct function on selected individual cells and setting in motion regulatory interactions. Non-selected cells then acquire new functions along the spatio-temporal developmental trajectory. The differentiation process has many stochastic components. Meiotic segregation, mitochondrial partitioning, X-inactivation and the dynamic DNA binding of transcription factor assemblies-all exhibit randomness. Non-random X-inactivation generally signals deleterious X-linked mutations. Correct neural wiring, such as retina to brain, arises through repeated confirmatory activity of connections made randomly. In immune system development, both B-cell antibody generation and the emergence of balanced T-cell categories begin through stochastic trial and error followed by functional selection. Aberrant selection processes lead to immune dysfunction. DNA sequence variants also arise through stochastic events: some involving environmental fluctuation (radiation or presence of pollutants), or genetic repair system malfunction. The phenotypic outcome of mutations is also fluid. Mutations may be advantageous in some circumstances, deleterious in others. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

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.

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.

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.

Recurrent heterozygous PAX6 missense variants cause severe bilateral microphthalmia via predictable effects on DNA-protein interaction.

PURPOSE: Most classical aniridia is caused by PAX6 haploinsufficiency. PAX6 missense variants can be hypomorphic or mimic haploinsufficiency. We hypothesized that missense variants also cause previously undescribed disease by altering the affinity and/or specificity of PAX6 genomic interactions. METHODS: We screened PAX6 in 372 individuals with bilateral microphthalmia, anophthalmia, or coloboma (MAC) from the Medical Research Council Human Genetics Unit eye malformation cohort (HGUeye) and reviewed data from the Deciphering Developmental Disorders study. We performed cluster analysis on PAX6-associated ocular phenotypes by variant type and molecular modeling of the structural impact of 86 different PAX6 causative missense variants. RESULTS: Eight different PAX6 missense variants were identified in 17 individuals (15 families) with MAC, accounting for 4% (15/372) of our cohort. Seven altered the paired domain (p.[Arg26Gln]x1, p.[Gly36Val]x1, p.[Arg38Trp]x2, p.[Arg38Gln]x1, p.[Gly51Arg]x2, p.[Ser54Arg]x2, p.[Asn124Lys]x5) and one the homeodomain (p.[Asn260Tyr]x1). p.Ser54Arg and p.Asn124Lys were exclusively associated with severe bilateral microphthalmia. MAC-associated variants were predicted to alter but not ablate DNA interaction, consistent with the electrophoretic mobility shifts observed using mutant paired domains with well-characterized PAX6-binding sites. We found no strong evidence for novel PAX6-associated extraocular disease. CONCLUSION: Altering the affinity and specificity of PAX6-binding genome-wide provides a plausible mechanism for the worse-than-null effects of MAC-associated missense variants.

Genome sequencing-the dawn of a game-changing era.

The development of genome sequencing technologies has revolutionized the biological sciences in ways which could not have been imagined at the time. This article sets out to document the dawning of the age of genomics and to consider the impact of this revolution on biological investigation, our understanding of life, and the relationship between science and society.

Unique noncoding variants upstream of PRDM13 are associated with a spectrum of developmental retinal dystrophies including progressive bifocal chorioretinal atrophy.

The autosomal dominant progressive bifocal chorioretinal atrophy (PBCRA) disease locus has been mapped to chromosome 6q14-16.2 that overlaps the North Carolina macular dystrophy (NCMD) locus MCDR1. NCMD is a nonprogressive developmental macular dystrophy, in which variants upstream of PRDM13 have been implicated. Whole genome sequencing was performed to interrogate structural variants (SVs) and single nucleotide variants (SNVs) in eight individuals, six affected individuals from two families with PBCRA, and two individuals from an additional family with a related developmental macular dystrophy. A SNV (chr6:100,046,804T>C), located 7.8 kb upstream of the PRDM13 gene, was shared by all PBCRA-affected individuals in the disease locus. Haplotype analysis suggested that the variant arose independently in the two families. The two affected individuals from Family 3 were screened for rare variants in the PBCRA and NCMD loci. This revealed a de novo variant in the proband, 21 bp from the first SNV (chr6:100,046,783A>C). This study expands the noncoding variant spectrum upstream of PRDM13 and suggests altered spatio-temporal expression of PRDM13 as a candidate disease mechanism in the phenotypically distinct but related conditions, NCMD and PBCRA.

Duplication events downstream of IRX1 cause North Carolina macular dystrophy at the MCDR3 locus.

Autosomal dominant North Carolina macular dystrophy (NCMD) is believed to represent a failure of macular development. The disorder has been linked to two loci, MCDR1 (chromosome 6q16) and MCDR3 (chromosome 5p15-p13). Recently, non-coding variants upstream of PRDM13 (MCDR1) and a duplication including IRX1 (MCDR3) have been identified. However, the underlying disease-causing mechanism remains uncertain. Through a combination of sequencing studies on eighteen NCMD families, we report two novel overlapping duplications at the MCDR3 locus, in a gene desert downstream of IRX1 and upstream of ADAMTS16. One duplication of 43 kb was identified in nine families (with evidence for a shared ancestral haplotype), and another one of 45 kb was found in a single family. Three families carry the previously reported V2 variant (MCDR1), while five remain unsolved. The MCDR3 locus is thus refined to a shared region of 39 kb that contains DNAse hypersensitive sites active at a restricted time window during retinal development. Publicly available data confirmed expression of IRX1 and ADAMTS16 in human fetal retina, with IRX1 preferentially expressed in fetal macula. These findings represent a major advance in our understanding of the molecular genetics of NCMD and provide insights into the genetic pathways involved in human macular development.

Genome-wide linkage and haplotype sharing analysis implicates the MCDR3 locus as a candidate region for a developmental macular disorder in association with digit abnormalities.

BACKGROUND: Developmental macular disorders are a heterogeneous group of rare retinal conditions that can cause significant visual impairment from childhood. Among these disorders, autosomal dominant North Carolina macular dystrophy (NCMD) has been mapped to 6q16 (MCDR1) with recent support for a non-coding disease mechanism of PRDM13. A second locus on 5p15-5p13 (MCDR3) has been implicated in a similar phenotype, but the disease-causing mechanism still remains unknown. METHODS: Two families affected by a dominant developmental macular disorder that closely resembles NCMD in association with digit abnormalities were included in the study. Family members with available DNA were genotyped using the Affymetrix GeneChip Human Mapping 250K Sty array. A parametric multipoint linkage analysis assuming a fully penetrant dominant model was performed using MERLIN. Haplotype sharing analysis was carried out using the non-parametric Homozygosity Haplotype method. Whole-exome sequencing was conducted on selected affected individuals. RESULTS: Linkage analysis excluded MCDR1 from the candidate regions (LOD < -2). There was suggestive linkage (LOD = 2.7) at two loci, including 9p24.1 and 5p15.32 that overlapped with MCDR3. The haplotype sharing analysis in one of the families revealed a 5 cM shared IBD segment at 5p15.32 (p value = 0.004). Whole-exome sequencing did not provide conclusive evidence for disease-causing alleles. CONCLUSIONS: These findings do not exclude that this phenotype may be allelic with NCMD MCDR3 at 5p15 and leave the possibility of a non-coding disease mechanism, in keeping with recent findings on 6q16. Further studies, including whole-genome sequencing, may help elucidate the underlying genetic cause of this phenotype and shed light on macular development and function.

A recurrent de novo mutation in ACTG1 causes isolated ocular coloboma.

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.

KDM3A coordinates actin dynamics with intraflagellar transport to regulate cilia stability.

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.

PAX6, brain structure and function in human adults: advanced MRI in aniridia.

OBJECTIVE: 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. METHODS: We used advanced magnetic resonance imaging techniques, including surface-based morphometry and region-of-interest analysis in adult humans heterozygously mutated for PAX6 mutations (n = 19 subjects and n = 21 controls). Using immunohistochemistry, we also studied PAX6 expression in the adult brain tissue of healthy subjects (n = 4) and patients with epilepsy (n = 42), some of whom had focal injuries due to intracranial electrode track placement (n = 17). RESULTS: There were significant reductions in frontoparietal cortical area after correcting for age and intracranial volume. A greater decline in thickness of the frontoparietal cortex with age, in subjects with PAX6 mutations compared to controls, correlated with age-corrected, accelerated decline in working memory. These results also demonstrate genotypic effects: those subjects with the most severe genotypes have the most widespread differences compared with controls. We also demonstrated significant increases in PAX6-expressing cells in response to acute injury in the adult human brain. INTERPRETATION: These findings suggest a role for PAX6 in the maintenance and consequent functioning of the adult brain, homologous to that found in other tissues. This has significant implications for the understanding and treatment of neurodegenerative diseases.

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.

Functional assessment of disease-associated regulatory variants in vivo using a versatile dual colour transgenesis strategy in zebrafish.

Disruption of gene regulation by sequence variation in non-coding regions of the genome is now recognised as a significant cause of human disease and disease susceptibility. Sequence variants in cis-regulatory elements (CREs), the primary determinants of spatio-temporal gene regulation, can alter transcription factor binding sites. While technological advances have led to easy identification of disease-associated CRE variants, robust methods for discerning functional CRE variants from background variation are lacking. Here we describe an efficient dual-colour reporter transgenesis approach in zebrafish, simultaneously allowing detailed in vivo comparison of spatio-temporal differences in regulatory activity between putative CRE variants and assessment of altered transcription factor binding potential of the variant. We validate the method on known disease-associated elements regulating SHH, PAX6 and IRF6 and subsequently characterise novel, ultra-long-range SOX9 enhancers implicated in the craniofacial abnormality Pierre Robin Sequence. The method provides a highly cost-effective, fast and robust approach for simultaneously unravelling in a single assay whether, where and when in embryonic development a disease-associated CRE-variant is affecting its regulatory function.

Long-range evolutionary constraints reveal cis-regulatory interactions on the human X chromosome.

Enhancers can regulate the transcription of genes over long genomic distances. This is thought to lead to selection against genomic rearrangements within such regions that may disrupt this functional linkage. Here we test this concept experimentally using the human X chromosome. We describe a scoring method to identify evolutionary maintenance of linkage between conserved noncoding elements and neighbouring genes. Chromatin marks associated with enhancer function are strongly correlated with this linkage score. We test >1,000 putative enhancers by transgenesis assays in zebrafish to ascertain the identity of the target gene. The majority of active enhancers drive a transgenic expression in a pattern consistent with the known expression of a linked gene. These results show that evolutionary maintenance of linkage is a reliable predictor of an enhancer's function, and provide new information to discover the genetic basis of diseases caused by the mis-regulation of gene expression.

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.

Identification of novel craniofacial regulatory domains located far upstream of SOX9 and disrupted in Pierre Robin sequence.

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.

Kdm3a lysine demethylase is an Hsp90 client required for cytoskeletal rearrangements during spermatogenesis.

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.