Disruption of autoregulatory feedback by a mutation in a remote, ultraconserved PAX6 enhancer causes aniridia

Deep mutational scanning quantifies DNA binding and predicts clinical outcomes of PAX6 variants

Nonsense and missense mutations in the transcription factor PAX6 cause a wide range of eye development defects, including aniridia, microphthalmia and coloboma. To understand how changes of PAX6:DNA binding cause these phenotypes, we combined saturation mutagenesis of the paired domain of PAX6 with a yeast one-hybrid (Y1H) assay in which expression of a PAX6-GAL4 fusion gene drives antibiotic resistance. We quantified binding of more than 2700 single amino-acid variants to two DNA sequence elements. Mutations in DNA-facing residues of the N-terminal subdomain and linker region were most detrimental, as were mutations to prolines and to negatively charged residues. Many variants caused sequence-specific molecular gain-of-function effects, including variants in position 71 that increased binding to the LE9 enhancer but decreased binding to a SELEX-derived binding site. In the absence of antibiotic selection, variants that retained DNA binding slowed yeast growth, likely because such variants perturbed the yeast transcriptome. Benchmarking against known patient variants and applying ACMG/AMP guidelines to variant classification, we obtained supporting-to-moderate evidence that 977 variants are likely pathogenic and 1306 are likely benign. Our analysis shows that most pathogenic mutations in the paired domain of PAX6 can be explained simply by the effects of these mutations on PAX6:DNA association, and establishes Y1H as a generalisable assay for the interpretation of variant effects in transcription factors.

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

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

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'.

Functional categorization of gene regulatory variants that cause Mendelian conditions

Much of our current understanding of rare human diseases is driven by coding genetic variants. However, non-coding genetic variants play a pivotal role in numerous rare human diseases, resulting in diverse functional impacts ranging from altered gene regulation, splicing, and/or transcript stability. With the increasing use of genome sequencing in clinical practice, it is paramount to have a clear framework for understanding how non-coding genetic variants cause disease. To this end, we have synthesized the literature on hundreds of non-coding genetic variants that cause rare Mendelian conditions via the disruption of gene regulatory patterns and propose a functional classification system. Specifically, we have adapted the functional classification framework used for coding variants (i.e., loss-of-function, gain-of-function, and dominant-negative) to account for features unique to non-coding gene regulatory variants. We identify that non-coding gene regulatory variants can be split into three distinct categories by functional impact: (1) non-modular loss-of-expression (LOE) variants; (2) modular loss-of-expression (mLOE) variants; and (3) gain-of-ectopic-expression (GOE) variants. Whereas LOE variants have a direct corollary with coding loss-of-function variants, mLOE and GOE variants represent disease mechanisms that are largely unique to non-coding variants. These functional classifications aim to provide a unified terminology for categorizing the functional impact of non-coding variants that disrupt gene regulatory patterns in Mendelian conditions.

Functional associations of evolutionarily recent human genes exhibit sensitivity to the 3D genome landscape and disease

Genome organization is intricately tied to regulating genes and associated cell fate decisions. In this study, we examine the positioning and functional significance of human genes, grouped by their evolutionary age, within the 3D organization of the genome. We reveal that genes of different evolutionary origin have distinct positioning relationships with both domains and loop anchors, and remarkably consistent relationships with boundaries across cell types. While the functional associations of each group of genes are primarily cell type-specific, such associations of conserved genes maintain greater stability across 3D genomic features and disease than recently evolved genes. Furthermore, the expression of these genes across various tissues follows an evolutionary progression, such that RNA levels increase from young genes to ancient genes. Thus, the distinct relationships of gene evolutionary age, function, and positioning within 3D genomic features contribute to tissue-specific gene regulation in development and disease.

Structural Variant Disrupting the Expression of the Remote FOXC1 Gene in a Patient with Syndromic Complex Microphthalmia

Ocular malformations (OMs) arise from early defects during embryonic eye development. Despite the identification of over 100 genes linked to this heterogeneous group of disorders, the genetic cause remains unknown for half of the individuals following Whole-Exome Sequencing. Diagnosis procedures are further hampered by the difficulty of studying samples from clinically relevant tissue, which is one of the main obstacles in OMs. Whole-Genome Sequencing (WGS) to screen for non-coding regions and structural variants may unveil new diagnoses for OM individuals. In this study, we report a patient exhibiting a syndromic OM with a de novo 3.15 Mb inversion in the 6p25 region identified by WGS. This balanced structural variant was located 100 kb away from the FOXC1 gene, previously associated with ocular defects in the literature. We hypothesized that the inversion disrupts the topologically associating domain of FOXC1 and impairs the expression of the gene. Using a new type of samples to study transcripts, we were able to show that the patient presented monoallelic expression of FOXC1 in conjunctival cells, consistent with the abolition of the expression of the inverted allele. This report underscores the importance of investigating structural variants, even in non-coding regions, in individuals affected by ocular malformations.

Multi-omics analysis in human retina uncovers ultraconserved cis-regulatory elements at rare eye disease loci

Cross-species genome comparisons have revealed a substantial number of ultraconserved non-coding elements (UCNEs). Several of these elements have proved to be essential tissue- and cell type-specific cis-regulators of developmental gene expression. Here, we characterize a set of UCNEs as candidate CREs (cCREs) during retinal development and evaluate the contribution of their genomic variation to rare eye diseases, for which pathogenic non-coding variants are emerging. Integration of bulk and single-cell retinal multi-omics data reveals 594 genes under potential cis-regulatory control of UCNEs, of which 45 are implicated in rare eye disease. Mining of candidate cis-regulatory UCNEs in WGS data derived from the rare eye disease cohort of Genomics England reveals 178 ultrarare variants within 84 UCNEs associated with 29 disease genes. Overall, we provide a comprehensive annotation of ultraconserved non-coding regions acting as cCREs during retinal development which can be targets of non-coding variation underlying rare eye diseases.

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.

Exploring noncoding variants in genetic diseases: from detection to functional insights

Previous studies on genetic diseases predominantly focused on protein-coding variations, overlooking the vast noncoding regions in the human genome. The development of high-throughput sequencing technologies and functional genomics tools has enabled the systematic identification of functional noncoding variants. These variants can impact gene expression, regulation, and chromatin conformation, thereby contributing to disease pathogenesis. Understanding the mechanisms that underlie the impact of noncoding variants on genetic diseases is indispensable for the development of precisely targeted therapies and the implementation of personalized medicine strategies. The intricacies of noncoding regions introduce a multitude of challenges and research opportunities. In this review, we introduce a spectrum of noncoding variants involved in genetic diseases, along with research strategies and advanced technologies for their precise identification and in-depth understanding of the complexity of the noncoding genome. We will delve into the research challenges and propose potential solutions for unraveling the genetic basis of rare and complex diseases.

A cell type-aware framework for nominating non-coding variants in Mendelian regulatory disorders

Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN) cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generated single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. Seventy-five percent of elements (44 of 59) validated in an in vivo transgenic reporter assay, demonstrating that single cell accessibility is a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieved significant reduction in our variant search space and nominated candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1, and EBF3 - as well as new candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work provides novel non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

Cis-Regulatory Elements in Mammals

In cis-regulatory elements, enhancers and promoters with complex molecular interactions are used to coordinate gene transcription through physical proximity and chemical modifications. These processes subsequently influence the phenotypic characteristics of an organism. An in-depth exploration of enhancers and promoters can substantially enhance our understanding of gene regulatory networks, shedding new light on mammalian development, evolution and disease pathways. In this review, we provide a comprehensive overview of the intrinsic structural attributes, detection methodologies as well as the operational mechanisms of enhancers and promoters, coupled with the relevant novel and innovative investigative techniques used to explore their actions. We further elucidated the state-of-the-art research on the roles of enhancers and promoters in the realms of mammalian development, evolution and disease, and we conclude with forward-looking insights into prospective research avenues.

EnhancerTracker: Comparing cell-type-specific enhancer activity of DNA sequence triplets via an ensemble of deep convolutional neural networks

Motivation

Transcriptional enhancers - unlike promoters - are unrestrained by distance or strand orientation with respect to their target genes, making their computational identification a challenge. Further, there are insufficient numbers of confirmed enhancers for many cell types, preventing robust training of machine-learning-based models for enhancer prediction for such cell types.

Results

We present EnhancerTracker , a novel tool that leverages an ensemble of deep separable convolutional neural networks to identify cell-type-specific enhancers with the need of only two confirmed enhancers. EnhancerTracker is trained, validated, and tested on 52,789 putative enhancers obtained from the FANTOM5 Project and control sequences derived from the human genome. Unlike available tools, which accept one sequence at a time, the input to our tool is three sequences; the first two are enhancers active in the same cell type. EnhancerTracker outputs 1 if the third sequence is an enhancer active in the same cell type(s) where the first two enhancers are active. It outputs 0 otherwise. On a held-out set (15%), EnhancerTracker achieved an accuracy of 64%, a specificity of 93%, a recall of 35%, a precision of 84%, and an F1 score of 49%.

Availability and implementation

https://github.com/BioinformaticsToolsmith/EnhancerTracker.

Contact

hani.girgis@tamuk.edu.

Unique activities of two overlapping PAX6 retinal enhancers

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

Disruption of regulatory domains and novel transcripts as disease-causing mechanisms

Deletions, duplications, insertions, inversions, and translocations, collectively called structural variations (SVs), affect more base pairs of the genome than any other sequence variant. The recent technological advancements in genome sequencing have enabled the discovery of tens of thousands of SVs per human genome. These SVs primarily affect non-coding DNA sequences, but the difficulties in interpreting their impact limit our understanding of human disease etiology. The functional annotation of non-coding DNA sequences and methodologies to characterize their three-dimensional (3D) organization in the nucleus have greatly expanded our understanding of the basic mechanisms underlying gene regulation, thereby improving the interpretation of SVs for their pathogenic impact. Here, we discuss the various mechanisms by which SVs can result in altered gene regulation and how these mechanisms can result in rare genetic disorders. Beyond changing gene expression, SVs can produce novel gene-intergenic fusion transcripts at the SV breakpoints.

The Genetic Determinants of Axial Length: From Microphthalmia to High Myopia in Childhood

The axial length of the eye is critical for normal visual function by enabling light to precisely focus on the retina. The mean axial length of the adult human eye is 23.5 mm, but the molecular mechanisms regulating ocular axial length remain poorly understood. Underdevelopment can lead to microphthalmia (defined as a small eye with an axial length of less than 19 mm at 1 year of age or less than 21 mm in adulthood) within the first trimester of pregnancy. However, continued overgrowth can lead to axial high myopia (an enlarged eye with an axial length of 26.5 mm or more) at any age. Both conditions show high genetic and phenotypic heterogeneity associated with significant visual morbidity worldwide. More than 90 genes can contribute to microphthalmia, and several hundred genes are associated with myopia, yet diagnostic yields are low. Crucially, the genetic pathways underpinning the specification of eye size are only now being discovered, with evidence suggesting that shared molecular pathways regulate under- or overgrowth of the eye. Improving our mechanistic understanding of axial length determination will help better inform us of genotype-phenotype correlations in both microphthalmia and myopia, dissect gene-environment interactions in myopia, and develop postnatal therapies that may influence overall eye growth.

A novel PAX6 variant as the cause of aniridia in a Chinese patient with SRRRD

BACKGROUND

The genotype characteristics and their associated clinical phenotypes in patients with aniridia were analyzed to explore pathogenic variants using whole-exome sequencing.

METHODS

One patient with aniridia was enrolled at the Beijing Tongren Hospital. Comprehensive ophthalmic and general examinations were performed on the patient. DNA was extracted from the patient, and whole-exome sequencing was performed to identify the causative variant. The pathogenicity of the variant was predicted using in silico analysis and evaluated according to American College of Medical Genetics and Genomics guidelines. Relationships between genetic variants and clinical features were analyzed.

RESULTS

In addition to the classical aniridia phenotype showing complete iris aplasia, foveal hypoplasia, and ectopic lentis, the patient also exhibited spontaneous reattachment rhegmatogenous retinal detachment (SRRRD). Whole-exome sequencing identified a novel heterozygous variant, exon8:c.640_646del:p.R214Pfs*28.

CONCLUSIONS

The present study broadens the range of genetic variants described in aniridia and presents an aniridia patient with SRRRD.

EyeG2P: an automated variant filtering approach improves efficiency of diagnostic genomic testing for inherited ophthalmic disorders

BACKGROUND

Genomic variant prioritisation is one of the most significant bottlenecks to mainstream genomic testing in healthcare. Tools to improve precision while ensuring high recall are critical to successful mainstream clinical genomic testing, in particular for whole genome sequencing where millions of variants must be considered for each patient.

METHODS

We developed EyeG2P, a publicly available database and web application using the Ensembl Variant Effect Predictor. EyeG2P is tailored for efficient variant prioritisation for individuals with inherited ophthalmic conditions. We assessed the sensitivity of EyeG2P in 1234 individuals with a broad range of eye conditions who had previously received a confirmed molecular diagnosis through routine genomic diagnostic approaches. For a prospective cohort of 83 individuals, we assessed the precision of EyeG2P in comparison with routine diagnostic approaches. For 10 additional individuals, we assessed the utility of EyeG2P for whole genome analysis.

RESULTS

EyeG2P had 99.5% sensitivity for genomic variants previously identified as clinically relevant through routine diagnostic analysis (n=1234 individuals). Prospectively, EyeG2P enabled a significant increase in precision (35% on average) in comparison with routine testing strategies (p<0.001). We demonstrate that incorporation of EyeG2P into whole genome sequencing analysis strategies can reduce the number of variants for analysis to six variants, on average, while maintaining high diagnostic yield.

CONCLUSION

Automated filtering of genomic variants through EyeG2P can increase the efficiency of diagnostic testing for individuals with a broad range of inherited ophthalmic disorders.

Congenital aniridia beyond black eyes: From phenotype and novel genetic mechanisms to innovative therapeutic approaches

Congenital PAX6-aniridia, initially characterized by the absence of the iris, has progressively been shown to be associated with other developmental ocular abnormalities and systemic features making congenital aniridia a complex syndromic disorder rather than a simple isolated disease of the iris. Moreover, foveal hypoplasia is now recognized as a more frequent feature than complete iris hypoplasia and a major visual prognosis determinant, reversing the classical clinical picture of this disease. Conversely, iris malformation is also a feature of various anterior segment dysgenesis disorders caused by PAX6-related developmental genes, adding a level of genetic complexity for accurate molecular diagnosis of aniridia. Therefore, the clinical recognition and differential genetic diagnosis of PAX6-related aniridia has been revealed to be much more challenging than initially thought, and still remains under-investigated. Here, we update specific clinical features of aniridia, with emphasis on their genotype correlations, as well as provide new knowledge regarding the PAX6 gene and its mutational spectrum, and highlight the beneficial utility of clinically implementing targeted Next-Generation Sequencing combined with Whole-Genome Sequencing to increase the genetic diagnostic yield of aniridia. We also present new molecular mechanisms underlying aniridia and aniridia-like phenotypes. Finally, we discuss the appropriate medical and surgical management of aniridic eyes, as well as innovative therapeutic options. Altogether, these combined clinical-genetic approaches will help to accelerate time to diagnosis, provide better determination of the disease prognosis and management, and confirm eligibility for future clinical trials or genetic-specific therapies.

A multi-omics atlas of the human retina at single-cell resolution

Cell classes in the human retina are highly heterogeneous with their abundance varying by several orders of magnitude. Here, we generated and integrated a multi-omics single-cell atlas of the adult human retina, including more than 250,000 nuclei for single-nuclei RNA-seq and 137,000 nuclei for single-nuclei ATAC-seq. Cross-species comparison of the retina atlas among human, monkey, mice, and chicken revealed relatively conserved and non-conserved types. Interestingly, the overall cell heterogeneity in primate retina decreases compared with that of rodent and chicken retina. Through integrative analysis, we identified 35,000 distal cis-element-gene pairs, constructed transcription factor (TF)-target regulons for more than 200 TFs, and partitioned the TFs into distinct co-active modules. We also revealed the heterogeneity of the cis-element-gene relationships in different cell types, even from the same class. Taken together, we present a comprehensive single-cell multi-omics atlas of the human retina as a resource that enables systematic molecular characterization at individual cell-type resolution.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Novel variants in the PAX6 gene related to isolated aniridia

Aniridia, which is a rare congenital defect of the eye, consists of iris hypoplasia or aplasia, and additional ocular abnormalities. It is most commonly caused by autosomal dominant PAX6 gene mutations. However, in about 30% of cases, it is associated with chromosomal rearrangements in the 11p13 region. The aim of this study was to identify the potential PAX6 gene variants, which could cause the isolated aniridia. Eight patients with isolated aniridia were included in this study. MLPA analysis allowed in the past to exclude large structural rearrangements of the PAX6 and adjacent genes like WT1. Blood samples were collected from the patients (and their families in familial cases) and genomic DNA was extracted from peripheral blood leukocytes and buccal cells. The amplification of the 11 exons of the PAX6 gene was performed. Bidirectional Sanger Sequencing was conducted for the identification of the potentially pathogenic variants, and for the segregation analysis of the identified variant in the family. The results were analyzed with the use of CodonCode Aligner software. In three patients, aniridia was sporadic, whereas in another five cases, the eye defect was familial. The potentially pathogenic variants in the PAX6 gene were found in 6 out of 8 patients with aniridia. We identified four known (c.781C > T, c.607C > T, and c.949C > T twice), and two novel variants (c.258_265del and c.495_496insG). Point mutations in the PAX6 gene are the most frequent cause of aniridia. The investigation of the genetic background of the disease is essential for patients to evaluate recurrence risk in the offspring.

Mendelian inheritance revisited: dominance and recessiveness in medical genetics

Understanding the consequences of genotype for phenotype (which ranges from molecule-level effects to whole-organism traits) is at the core of genetic diagnostics in medicine. Many measures of the deleteriousness of individual alleles exist, but these have limitations for predicting the clinical consequences. Various mechanisms can protect the organism from the adverse effects of functional variants, especially when the variant is paired with a wild type allele. Understanding why some alleles are harmful in the heterozygous state - representing dominant inheritance - but others only with the biallelic presence of pathogenic variants - representing recessive inheritance - is particularly important when faced with the deluge of rare genetic alterations identified by high throughput DNA sequencing. Both awareness of the specific quantitative and/or qualitative effects of individual variants and the elucidation of allelic and non-allelic interactions are essential to optimize genetic diagnosis and counselling.

Dynamic changes in whole genome DNA methylation, chromatin and gene expression during mouse lens differentiation

BACKGROUND

Cellular differentiation is marked by temporally and spatially coordinated gene expression regulated at multiple levels. DNA methylation represents a universal mechanism to control chromatin organization and its accessibility. Cytosine methylation of CpG dinucleotides regulates binding of methylation-sensitive DNA-binding transcription factors within regulatory regions of transcription, including promoters and distal enhancers. Ocular lens differentiation represents an advantageous model system to examine these processes as lens comprises only two cell types, the proliferating lens epithelium and postmitotic lens fiber cells all originating from the epithelium.

RESULTS

Using whole genome bisulfite sequencing (WGBS) and microdissected lenses, we investigated dynamics of DNA methylation and chromatin changes during mouse lens fiber and epithelium differentiation between embryos (E14.5) and newborns (P0.5). Histone H3.3 variant chromatin landscapes were also generated for both P0.5 lens epithelium and fibers by chromatin immunoprecipitation followed by next generation sequencing (ChIP-seq). Tissue-specific features of DNA methylation patterns are demonstrated via comparative studies with embryonic stem (ES) cells and neural progenitor cells (NPCs) at Nanog, Pou5f1, Sox2, Pax6 and Six3 loci. Comparisons with ATAC-seq and RNA-seq data demonstrate that reduced methylation is associated with increased expression of fiber cell abundant genes, including crystallins, intermediate filament (Bfsp1 and Bfsp2) and gap junction proteins (Gja3 and Gja8), marked by high levels of histone H3.3 within their transcribed regions. Interestingly, Pax6-binding sites exhibited predominantly DNA hypomethylation in lens chromatin. In vitro binding of Pax6 proteins showed Pax6's ability to interact with sites containing one or two methylated CpG dinucleotides.

CONCLUSIONS

Our study has generated the first data on methylation changes between two different stages of mammalian lens development and linked these data with chromatin accessibility maps, presence of histone H3.3 and gene expression. Reduced DNA methylation correlates with expression of important genes involved in lens morphogenesis and lens fiber cell differentiation.

The Power of Clinical Diagnosis for Deciphering Complex Genetic Mechanisms in Rare Diseases

Complex genetic disease mechanisms, such as structural or non-coding variants, currently pose a substantial difficulty in frontline diagnostic tests. They thus may account for most unsolved rare disease patients regardless of the clinical phenotype. However, the clinical diagnosis can narrow the genetic focus to just a couple of genes for patients with well-established syndromes defined by prominent physical and/or unique biochemical phenotypes, allowing deeper analyses to consider complex genetic origin. Then, clinical-diagnosis-driven genome sequencing strategies may expedite the development of testing and analytical methods to account for complex disease mechanisms as well as to advance functional assays for the confirmation of complex variants, clinical management, and the development of new therapies.

Multi-omics approach dissects cis-regulatory mechanisms underlying North Carolina macular dystrophy, a retinal enhanceropathy

North Carolina macular dystrophy (NCMD) is a rare autosomal-dominant disease affecting macular development. The disease is caused by non-coding single-nucleotide variants (SNVs) in two hotspot regions near PRDM13 and by duplications in two distinct chromosomal loci, overlapping DNase I hypersensitive sites near either PRDM13 or IRX1. To unravel the mechanisms by which these variants cause disease, we first established a genome-wide multi-omics retinal database, RegRet. Integration of UMI-4C profiles we generated on adult human retina then allowed fine-mapping of the interactions of the PRDM13 and IRX1 promoters and the identification of eighteen candidate cis-regulatory elements (cCREs), the activity of which was investigated by luciferase and Xenopus enhancer assays. Next, luciferase assays showed that the non-coding SNVs located in the two hotspot regions of PRDM13 affect cCRE activity, including two NCMD-associated non-coding SNVs that we identified herein. Interestingly, the cCRE containing one of these SNVs was shown to interact with the PRDM13 promoter, demonstrated in vivo activity in Xenopus, and is active at the developmental stage when progenitor cells of the central retina exit mitosis, suggesting that this region is a PRDM13 enhancer. Finally, mining of single-cell transcriptional data of embryonic and adult retina revealed the highest expression of PRDM13 and IRX1 when amacrine cells start to synapse with retinal ganglion cells, supporting the hypothesis that altered PRDM13 or IRX1 expression impairs interactions between these cells during retinogenesis. Overall, this study provides insight into the cis-regulatory mechanisms of NCMD and supports that this condition is a retinal enhanceropathy.

Natural and Experimental Rewiring of Gene Regulatory Regions

The successful development and ongoing functioning of complex organisms depend on the faithful execution of the genetic code. A critical step in this process is the correct spatial and temporal expression of genes. The highly orchestrated transcription of genes is controlled primarily by cis-regulatory elements: promoters, enhancers, and insulators. The medical importance of this key biological process can be seen by the frequency with which mutations and inherited variants that alter cis-regulatory elements lead to monogenic and complex diseases and cancer. Here, we provide an overview of the methods available to characterize and perturb gene regulatory circuits. We then highlight mechanisms through which regulatory rewiring contributes to disease, and conclude with a perspective on how our understanding of gene regulation can be used to improve human health.

Pancreas agenesis mutations disrupt a lead enhancer controlling a developmental enhancer cluster

Sequence variants in cis-acting enhancers are important for polygenic disease, but their role in Mendelian disease is poorly understood. Redundancy between enhancers that regulate the same gene is thought to mitigate the pathogenic impact of enhancer mutations. Recent findings, however, have shown that loss-of-function mutations in a single enhancer near PTF1A cause pancreas agenesis and neonatal diabetes. Using mouse and human genetic models, we show that this enhancer activates an entire PTF1A enhancer cluster in early pancreatic multipotent progenitors. This leading role, therefore, precludes functional redundancy. We further demonstrate that transient expression of PTF1A in multipotent progenitors sets in motion an epigenetic cascade that is required for duct and endocrine differentiation. These findings shed insights into the genome regulatory mechanisms that drive pancreas differentiation. Furthermore, they reveal an enhancer that acts as a regulatory master key and is thus vulnerable to pathogenic loss-of-function mutations.

•

The pancreas agenesis enhancer (EnhP) activates PTF1A in early pancreatic progenitors

•

EnhP also activates other progenitor PTF1A enhancers

•

This master key function explains why EnhP is vulnerable to loss-of-function mutations

•

Transient PTF1A expression in progenitors controls pancreas growth and endocrinogenesis

Unusual Presentation in WAGR Syndrome: Expanding the Phenotypic and Genotypic Spectrum of the Diseases

The deletion of chromosome 11p13 involving the WT1 and PAX6 genes has been shown to cause WAGR syndrome (OMIM #194072), a rare genetic disorder that features Wilms’ tumor, aniridia, genitourinary anomalies, as well as mental retardation. In this study, we expand the genotypic and phenotypic spectrum of WAGR syndrome by reporting on six patients from six unrelated families with different de novo deletions located on chromosome 11p13. Very rare phenotypes of lens automated absorption and lens thinning were detected in four of the six patients. We assessed the involvement of the ARL14EP gene in patients with and without severe lens abnormalities and found that its deletion may worsen the lens abnormalities in these patients.

Recommendations for clinical interpretation of variants found in non-coding regions of the genome

BACKGROUND

The majority of clinical genetic testing focuses almost exclusively on regions of the genome that directly encode proteins. The important role of variants in non-coding regions in penetrant disease is, however, increasingly being demonstrated, and the use of whole genome sequencing in clinical diagnostic settings is rising across a large range of genetic disorders. Despite this, there is no existing guidance on how current guidelines designed primarily for variants in protein-coding regions should be adapted for variants identified in other genomic contexts.

METHODS

We convened a panel of nine clinical and research scientists with wide-ranging expertise in clinical variant interpretation, with specific experience in variants within non-coding regions. This panel discussed and refined an initial draft of the guidelines which were then extensively tested and reviewed by external groups.

RESULTS

We discuss considerations specifically for variants in non-coding regions of the genome. We outline how to define candidate regulatory elements, highlight examples of mechanisms through which non-coding region variants can lead to penetrant monogenic disease, and outline how existing guidelines can be adapted for the interpretation of these variants.

CONCLUSIONS

These recommendations aim to increase the number and range of non-coding region variants that can be clinically interpreted, which, together with a compatible phenotype, can lead to new diagnoses and catalyse the discovery of novel disease mechanisms.

A Journey Through Genetics to Biology

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

Cell-specific cis-regulatory elements and mechanisms of non-coding genetic disease in human retina and retinal organoids

Cis-regulatory elements (CREs) play a critical role in the development and disease-states of all human cell types. In the retina, CREs have been implicated in several inherited disorders. To better characterize human retinal CREs, we performed single-nucleus assay for transposase-accessible chromatin sequencing (snATAC-seq) and single-nucleus RNA sequencing (snRNA-seq) on the developing and adult human retina and on induced pluripotent stem cell (iPSC)-derived retinal organoids. These analyses identified developmentally dynamic, cell-class-specific CREs, enriched transcription-factor-binding motifs, and putative target genes. CREs in the retina and organoids are highly correlated at the single-cell level, and this supports the use of organoids as a model for studying disease-associated CREs. As a proof of concept, we disrupted a disease-associated CRE at 5q14.3, confirming its principal target gene as the miR-9-2 primary transcript and demonstrating its role in neurogenesis and gene regulation in mature glia. This study provides a resource for characterizing human retinal CREs and showcases organoids as a model to study the function of CREs that influence development and disease.

Perfect and imperfect views of ultraconserved sequences

Across the human genome, there are nearly 500 'ultraconserved' elements: regions of at least 200 contiguous nucleotides that are perfectly conserved in both the mouse and rat genomes. Remarkably, the majority of these sequences are non-coding, and many can function as enhancers that activate tissue-specific gene expression during embryonic development. From their first description more than 15 years ago, their extreme conservation has both fascinated and perplexed researchers in genomics and evolutionary biology. The intrigue around ultraconserved elements only grew with the observation that they are dispensable for viability. Here, we review recent progress towards understanding the general importance and the specific functions of ultraconserved sequences in mammalian development and human disease and discuss possible explanations for their extreme conservation.

Deciphering the Retinal Epigenome during Development, Disease and Reprogramming: Advancements, Challenges and Perspectives

Retinal neurogenesis is driven by concerted actions of transcription factors, some of which are expressed in a continuum and across several cell subtypes throughout development. While seemingly redundant, many factors diversify their regulatory outcome on gene expression, by coordinating variations in chromatin landscapes to drive divergent retinal specification programs. Recent studies have furthered the understanding of the epigenetic contribution to the progression of age-related macular degeneration, a leading cause of blindness in the elderly. The knowledge of the epigenomic mechanisms that control the acquisition and stabilization of retinal cell fates and are evoked upon damage, holds the potential for the treatment of retinal degeneration. Herein, this review presents the state-of-the-art approaches to investigate the retinal epigenome during development, disease, and reprogramming. A pipeline is then reviewed to functionally interrogate the epigenetic and transcriptional networks underlying cell fate specification, relying on a truly unbiased screening of open chromatin states. The related work proposes an inferential model to identify gene regulatory networks, features the first footprinting analysis and the first tentative, systematic query of candidate pioneer factors in the retina ever conducted in any model organism, leading to the identification of previously uncharacterized master regulators of retinal cell identity, such as the nuclear factor I, NFI. This pipeline is virtually applicable to the study of genetic programs and candidate pioneer factors in any developmental context. Finally, challenges and limitations intrinsic to the current next-generation sequencing techniques are discussed, as well as recent advances in super-resolution imaging, enabling spatio-temporal resolution of the genome.

First evidence of SOX2 mutations in Peters' anomaly: lessons from molecular screening of 95 patients

Peters' anomaly (PA) is a rare anterior segment dysgenesis characterized by central corneal opacity and irido-lenticulo-corneal adhesions. Several genes are involved in syndromic or isolated PA (B3GLCT, PAX6, PITX3, FOXE3, CYP1B1). Some Copy Number Variations (CNVs) have also been occasionally reported. Despite this genetic heterogeneity, most of patients remain without genetic diagnosis. We retrieved a cohort of 95 individuals with PA and performed genotyping using a combination of Comparative genomic hybridization, whole genome, exome and targeted sequencing of 119 genes associated with ocular development anomalies. Causative genetic defects involving 12 genes and CNVs were identified for 1/3 of patients. Unsurprisingly, B3GLCT and PAX6 were the most frequently implicated genes, respectively in syndromic and isolated PA. Unexpectedly, the third gene involved in our cohort was SOX2, the major gene of micro-anophthalmia. Four unrelated patients with PA (isolated or with microphthalmia) were carrying pathogenic variants in this gene that was never associated with PA before. Here we described the largest cohort of PA patients ever reported. The genetic bases of PA are still to be explored as genetic diagnosis was unavailable for 2/3 of patients. Nevertheless, we showed here for the first time the involvement of SOX2 in PA, offering new evidence for its role in corneal transparency and anterior segment development. This article is protected by copyright. All rights reserved.

Developmental disorders caused by haploinsufficiency of transcriptional regulators: a perspective based on cell fate determination

Many human birth defects and neurodevelopmental disorders are caused by loss-of-function mutations in a single copy of transcription factor (TF) and chromatin regulator genes. Although this dosage sensitivity has long been known, how and why haploinsufficiency (HI) of transcriptional regulators leads to developmental disorders (DDs) is unclear. Here I propose the hypothesis that such DDs result from defects in cell fate determination that are based on disrupted bistability in the underlying gene regulatory network (GRN). Bistability, a crucial systems biology concept to model binary choices such as cell fate decisions, requires both positive feedback and ultrasensitivity, the latter often achieved through TF cooperativity. The hypothesis explains why dosage sensitivity of transcriptional regulators is an inherent property of fate decisions, and why disruption of either positive feedback or cooperativity in the underlying GRN is sufficient to cause disease. I present empirical and theoretical evidence in support of this hypothesis and discuss several issues for which it increases our understanding of disease, such as incomplete penetrance. The proposed framework provides a mechanistic, systems-level explanation of HI of transcriptional regulators, thus unifying existing theories, and offers new insights into outstanding issues of human disease. This article has an associated Future Leader to Watch interview with the author of the paper.

3D genome organization links non-coding disease-associated variants to genes

Genome sequencing has revealed over 300 million genetic variations in human populations. Over 90% of variants are single nucleotide polymorphisms (SNPs), the remainder include short deletions or insertions, and small numbers of structural variants. Hundreds of thousands of these variants have been associated with specific phenotypic traits and diseases through genome wide association studies which link significant differences in variant frequencies with specific phenotypes among large groups of individuals. Only 5% of disease-associated SNPs are located in gene coding sequences, with the potential to disrupt gene expression or alter of the function of encoded proteins. The remaining 95% of disease-associated SNPs are located in non-coding DNA sequences which make up 98% of the genome. The role of non-coding, disease-associated SNPs, many of which are located at considerable distances from any gene, was at first a mystery until the discovery that gene promoters regularly interact with distal regulatory elements to control gene expression. Disease-associated SNPs are enriched at the millions of gene regulatory elements that are dispersed throughout the non-coding sequences of the genome, suggesting they function as gene regulation variants. Assigning specific regulatory elements to the genes they control is not straightforward since they can be millions of base pairs apart. In this review we describe how understanding 3D genome organization can identify specific interactions between gene promoters and distal regulatory elements and how 3D genomics can link disease-associated SNPs to their target genes. Understanding which gene or genes contribute to a specific disease is the first step in designing rational therapeutic interventions.

Quantitative spatial and temporal assessment of regulatory element activity in zebrafish

Mutations or genetic variation in noncoding regions of the genome harbouring cis-regulatory elements (CREs), or enhancers, have been widely implicated in human disease and disease risk. However, our ability to assay the impact of these DNA sequence changes on enhancer activity is currently very limited because of the need to assay these elements in an appropriate biological context. Here, we describe a method for simultaneous quantitative assessment of the spatial and temporal activity of wild-type and disease-associated mutant human CRE alleles using live imaging in zebrafish embryonic development. We generated transgenic lines harbouring a dual-CRE dual-reporter cassette in a pre-defined neutral docking site in the zebrafish genome. The activity of each CRE allele is reported via expression of a specific fluorescent reporter, allowing simultaneous visualisation of where and when in development the wild-type allele is active and how this activity is altered by mutation.

Building a Mammalian Retina: An Eye on Chromatin Structure

Regulation of gene expression by chromatin structure has been under intensive investigation, establishing nuclear organization and genome architecture as a potent and effective means of regulating developmental processes. The substantial growth in our knowledge of the molecular mechanisms underlying retinogenesis has been powered by several genome-wide based tools that mapped chromatin organization at multiple cellular and biochemical levels. Studies profiling the retinal epigenome and transcriptome have allowed the systematic annotation of putative cis-regulatory elements associated with transcriptional programs that drive retinal neural differentiation, laying the groundwork to understand spatiotemporal retinal gene regulation at a mechanistic level. In this review, we outline recent advances in our understanding of the chromatin architecture in the mammalian retina during development and disease. We focus on the emerging roles of non-coding regulatory elements in controlling retinal cell-type specific transcriptional programs, and discuss potential implications in untangling the etiology of eye-related disorders.

Genetics and epidemiology of aniridia: Updated guidelines for genetic study

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

Role of Enhancers in Development and Diseases

Enhancers are cis-regulatory elements containing short DNA sequences that serve as binding sites for pioneer/regulatory transcription factors, thus orchestrating the regulation of genes critical for lineage determination. The activity of enhancer elements is believed to be determined by transcription factor binding, thus determining the cell state identity during development. Precise spatio-temporal control of the transcriptome during lineage specification requires the coordinated binding of lineage-specific transcription factors to enhancers. Thus, enhancers are the primary determinants of cell identity. Numerous studies have explored the role and mechanism of enhancers during development and disease, and various basic questions related to the functions and mechanisms of enhancers have not yet been fully answered. In this review, we discuss the recently published literature regarding the roles of enhancers, which are critical for various biological processes governing development. Furthermore, we also highlight that altered enhancer landscapes provide an essential context to understand the etiologies and mechanisms behind numerous complex human diseases, providing new avenues for effective enhancer-based therapeutic interventions.

Transcriptional enhancers and their communication with gene promoters

Transcriptional enhancers play a key role in the initiation and maintenance of gene expression programmes, particularly in metazoa. How these elements control their target genes in the right place and time is one of the most pertinent questions in functional genomics, with wide implications for most areas of biology. Here, we synthesise classic and recent evidence on the regulatory logic of enhancers, including the principles of enhancer organisation, factors that facilitate and delimit enhancer-promoter communication, and the joint effects of multiple enhancers. We show how modern approaches building on classic insights have begun to unravel the complexity of enhancer-promoter relationships, paving the way towards a quantitative understanding of gene control.

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.

Chromatin Conformation in Development and Disease

Chromatin domains and loops are important elements of chromatin structure and dynamics, but much remains to be learned about their exact biological role and nature. Topological associated domains and functional loops are key to gene expression and hold the answer to many questions regarding developmental decisions and diseases. Here, we discuss new findings, which have linked chromatin conformation with development, differentiation and diseases and hypothesized on various models while integrating all recent findings on how chromatin architecture affects gene expression during development, evolution and disease.

Molecular Characterization of Choroideremia-Associated Deletions Reveals an Unexpected Regulation of CHM Gene Transcription

Choroideremia (CHM) is a X-linked recessive chorioretinal dystrophy due to deficiency of the CHM gene product, i.e., Rab escort protein isoform 1 (REP1). To date, gene therapy for CHM has shown variable effectiveness, likely because the underlying pathogenic mechanisms as well as genotype-phenotype correlation are not yet fully known. Small nucleotide variants leading to premature termination codons (PTCs) are a major cause of CHM, but about 20% of patients has CHM gene deletions. To improve understanding of the disease mechanisms, we analyzed molecular features of seven deletions involving the CHM gene sequence. We mapped the deletion breakpoints by using polymerase chain reaction, sequencing and array comparative genomic hybridization; to identify rearrangement-promoting DNA sequences, we analyzed genomic architecture surrounding the breakpoint regions. Moreover, in some CHM patients with different mutation types, we measured transcript level of CHM and of CHML, encoding the REP2 isoform. Scattered along the whole CHM gene and in close proximity to the deletion breakpoints we found numerous repeat elements that generate a locus-specific rearrangement hot spot. Unexpectedly, patients with non-PTC variants had increased expression of the aberrant CHM mRNA; CHML expression was higher than normal in a patient lacking CHM and its putative regulatory sequences. This latest evidence suggests that mechanisms regulating CHM and CHML gene expression are worthy of further study, because their full knowledge could be also useful for developing effective therapies for this hitherto untreatable inherited retinal degeneration.

dbCNS: A New Database for Conserved Noncoding Sequences

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

Genetics and epidemiology of aniridia: Updated guidelines for genetic study

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

Genome sequencing in families with congenital limb malformations

The extensive clinical and genetic heterogeneity of congenital limb malformation calls for comprehensive genome-wide analysis of genetic variation. Genome sequencing (GS) has the potential to identify all genetic variants. Here we aim to determine the diagnostic potential of GS as a comprehensive one-test-for-all strategy in a cohort of undiagnosed patients with congenital limb malformations. We collected 69 cases (64 trios, 1 duo, 5 singletons) with congenital limb malformations with no molecular diagnosis after standard clinical genetic testing and performed genome sequencing. We also developed a framework to identify potential noncoding pathogenic variants. We identified likely pathogenic/disease-associated variants in 12 cases (17.4%) including four in known disease genes, and one repeat expansion in HOXD13. In three unrelated cases with ectrodactyly, we identified likely pathogenic variants in UBA2, establishing it as a novel disease gene. In addition, we found two complex structural variants (3%). We also identified likely causative variants in three novel high confidence candidate genes. We were not able to identify any noncoding variants. GS is a powerful strategy to identify all types of genomic variants associated with congenital limb malformation, including repeat expansions and complex structural variants missed by standard diagnostic approaches. In this cohort, no causative noncoding SNVs could be identified.

Diagnostic yield of targeted next-generation sequencing in infantile nystagmus syndrome

Background: Infantile nystagmus syndrome (INS) is a genetically heterogeneous disorder. Identifying genetic causes of INS would help clinicians to facilitate clinical diagnosis and provide appropriate treatment. The aim of this study was to determine the diagnostic utility of targeted next-generation sequencing (NGS) for INS.Materials and methods: We recruited 37 patients who were referred to the Neuro-ophthalmology clinics for evaluations of INS. NGS was performed using a targeted panel that included 98 candidate genes associated with INS. We identified pathogenic variants according to guidelines of the American College of Medical Genetics and Genomics. We also calculated the sensitivity and specificity of each clinical sign to assess the diagnostic yield of our gene panel.Results: After variant filtering, annotation, and interpretation, the potential pathogenic variants were detected in 13 of the 37 patients, achieving a molecular diagnostic rate of 35%. The identified genes were PAX6 (n = 4), FRMD7 (n = 4), GPR143 (n = 2), CACNA1F (n = 1), CNGA3 (n = 1) and GUCY2D (n = 1). In approximately 30% (n = 4) of the patients, the initial clinical diagnosis was revised after a molecular diagnosis was performed. The presence of a family history had the highest predictive power for a molecular diagnosis (sensitivity = 61.5%, specificity = 91.7%), and the sensitivity increased when the family history was considered together with one of two clinical signs such as pendular nystagmus waveforms or anterior segment dysgenesis.Conclusions: Our study shows that targeted NGS can be useful to determine a molecular diagnosis for patients with INS. Targeted NGS also helps to confirm a clinical diagnosis in atypical phenotypes or unresolved cases.

North Carolina Macular Dystrophy: Phenotypic Variability and Computational Analysis of Disease-Associated Noncoding Variants

Purpose

North Carolina macular dystrophy (NCMD) is an autosomal dominant, congenital disorder affecting the central retina. Here, we report clinical and genetic findings in three families segregating NCMD and use epigenomic datasets from human tissues to gain insights into the effect of NCMD-implicated variants.

Methods

Clinical assessment and genetic testing were performed. Publicly available transcriptomic and epigenomic datasets were analyzed and the activity-by-contact method for scoring enhancer elements and linking them to target genes was used.

Results

A previously described, heterozygous, noncoding variant upstream of the PRDM13 gene was detected in all six affected study participants (chr6:100,040,987G>C [GRCh37/hg19]). Interfamilial and intrafamilial variability were observed; the visual acuity ranged from 0.0 to 1.6 LogMAR and fundoscopic findings ranged from visually insignificant, confluent, drusen-like macular deposits to coloboma-like macular lesions. Variable degrees of peripheral retinal spots (which were easily detected on widefield retinal imaging) were observed in all study subjects. Notably, a 6-year-old patient developed choroidal neovascularization and required treatment with intravitreal bevacizumab injections. Computational analysis of the five single nucleotide variants that have been implicated in NCMD revealed that these noncoding changes lie within two putative enhancer elements; these elements are predicted to interact with PRDM13 in the developing human retina. PRDM13 was found to be expressed in the fetal retina, with greatest expression in the amacrine precursor cell population.

Conclusions

We provide further evidence supporting the role of PRDM13 dysregulation in the pathogenesis of NCMD and highlight the usefulness of widefield retinal imaging in individuals suspected to have this condition.

The rax homeobox gene is mutated in the eyeless axolotl, Ambystoma mexicanum

BACKGROUND

Vertebrate eye formation requires coordinated inductive interactions between different embryonic tissue layers, first described in amphibians. A network of transcription factors and signaling molecules controls these steps, with mutations causing severe ocular, neuronal, and craniofacial defects. In eyeless mutant axolotls, eye morphogenesis arrests at the optic vesicle stage, before lens induction, and development of ventral forebrain structures is disrupted.

RESULTS

We identified a 5-bp deletion in the rax (retina and anterior neural fold homeobox) gene, which was tightly linked to the recessive eyeless (e) axolotl locus in an F2 cross. This frameshift mutation, in exon 2, truncates RAX protein within the homeodomain (P154fs35X). Quantitative RNA analysis shows that mutant and wild-type rax transcripts are equally abundant in E/e embryos. Translation appears to initiate from dual start codons, via leaky ribosome scanning, a conserved feature among gnathostome RAX proteins. Previous data show rax is expressed in the optic vesicle and diencephalon, deeply conserved among metazoans, and required for eye formation in other species.

CONCLUSION

The eyeless axolotl mutation is a null allele in the rax homeobox gene, with primary defects in neural ectoderm, including the retinal and hypothalamic primordia.