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

Using computational approaches to enhance the interpretation of missense variants in the PAX6 gene

The PAX6 gene encodes a highly-conserved transcription factor involved in eye development. Heterozygous loss-of-function variants in PAX6 can cause a range of ophthalmic disorders including aniridia. A key molecular diagnostic challenge is that many PAX6 missense changes are presently classified as variants of uncertain significance. While computational tools can be used to assess the effect of genetic alterations, the accuracy of their predictions varies. Here, we evaluated and optimised the performance of computational prediction tools in relation to PAX6 missense variants. Through inspection of publicly available resources (including HGMD, ClinVar, LOVD and gnomAD), we identified 241 PAX6 missense variants that were used for model training and evaluation. The performance of ten commonly used computational tools was assessed and a threshold optimization approach was utilized to determine optimal cut-off values. Validation studies were subsequently undertaken using PAX6 variants from a local database. AlphaMissense, SIFT4G and REVEL emerged as the best-performing predictors; the optimized thresholds of these tools were 0.967, 0.025, and 0.772, respectively. Combining the prediction from these top-three tools resulted in lower performance compared to using AlphaMissense alone. Tailoring the use of computational tools by employing optimized thresholds specific to PAX6 can enhance algorithmic performance. Our findings have implications for PAX6 variant interpretation in clinical settings.

The Systemic Genotype-Phenotype Characterization of PAX6-Related Eye Disease in 164 Chinese Families

Purpose

This study aims to evaluate the genetic and phenotypic characteristics and elucidate the genotype-phenotype correlations of a large Chinese cohort with PAX6-related disorders.

Methods

Variants detected with exome sequencing were filtered through multistep bioinformatic and co-segregation analyses, and validated by Sanger sequencing. The related clinical data were collected, and cluster analysis and statistical analysis of the PAX6-related phenotypes across different variant groups were carried out. Parental mosaicism was investigated using cloning analysis and Droplet digital PCR.

Results

A total of 119 pathogenic or likely pathogenic PAX6 variants, including 74 truncation, 31 missense, and 14 others, were identified in 228 patients from 164 unrelated families. The most common phenotypes were foveal hypoplasia (97.8%), nystagmus (92.6%), aniridia (76.7%), cataract (36.8%), and iris hypoplasia (22.4%). Mosaicism ranging from 13.9% to 18.8% was identified in 3 unrelated patients' parents with relatively mild phenotypes. Missense variants in the linker region of the paired domain were associated with high myopia, whereas truncation variants in the homeodomain and proline-serine-threonine-rich domain were associated with hyperopia. Similarly, the degree of iris defects, visual acuity, and associated ocular comorbidity varied among the different types and locations of PAX6 variants.

Conclusions

Our data indicate that foveal hypoplasia but not aniridia is the most common sign of PAX6-related disorders, contributing to subtle iris changes that might easily be overlooked in clinical practice. Recognition of mosaicism in atypical cases or parents with very mild phenotypes is important in genetic counseling as their offspring are at increased risk of typical aniridia. Recognition of the genotype-phenotype relationship emphasizes involvement of PAX6 regulation in shaping complex ocular phenotypes.

Identification and Characterization of ATOH7-Regulated Target Genes and Pathways in Human Neuroretinal Development

The proneural transcription factor atonal basic helix-loop-helix transcription factor 7 (ATOH7) is expressed in early progenitors in the developing neuroretina. In vertebrates, this is crucial for the development of retinal ganglion cells (RGCs), as mutant animals show an almost complete absence of RGCs, underdeveloped optic nerves, and aberrations in retinal vessel development. Human mutations are rare and result in autosomal recessive optic nerve hypoplasia (ONH) or severe vascular changes, diagnosed as autosomal recessive persistent hyperplasia of the primary vitreous (PHPVAR). To better understand the role of ATOH7 in neuroretinal development, we created ATOH7 knockout and eGFP-expressing ATOH7 reporter human induced pluripotent stem cells (hiPSCs), which were differentiated into early-stage retinal organoids. Target loci regulated by ATOH7 were identified by Cleavage Under Targets and Release Using Nuclease with sequencing (CUT&RUN-seq) and differential expression by RNA sequencing (RNA-seq) of wildtype and mutant organoid-derived reporter cells. Additionally, single-cell RNA sequencing (scRNA-seq) was performed on whole organoids to identify cell type-specific genes. Mutant organoids displayed substantial deficiency in axon sprouting, reduction in RGCs, and an increase in other cell types. We identified 469 differentially expressed target genes, with an overrepresentation of genes belonging to axon development/guidance and Notch signaling. Taken together, we consolidate the function of human ATOH7 in guiding progenitor competence by inducing RGC-specific genes while inhibiting other cell fates. Furthermore, we highlight candidate genes responsible for ATOH7-associated optic nerve and retinovascular anomalies, which sheds light to potential future therapy targets for related disorders.

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.

Disruption of common ocular developmental pathways in patient-derived optic vesicle models of microphthalmia

Genetic perturbations influencing early eye development can result in microphthalmia, anophthalmia, and coloboma (MAC). Over 100 genes are associated with MAC, but little is known about common disease mechanisms. In this study, we generated induced pluripotent stem cell (iPSC)-derived optic vesicles (OVs) from two unrelated microphthalmia patients and healthy controls. At day 20, 35, and 50, microphthalmia patient OV diameters were significantly smaller, recapitulating the "small eye" phenotype. RNA sequencing (RNA-seq) analysis revealed upregulation of apoptosis-initiating and extracellular matrix (ECM) genes at day 20 and 35. Western blot and immunohistochemistry revealed increased expression of lumican, nidogen, and collagen type IV, suggesting ECM overproduction. Increased apoptosis was observed in microphthalmia OVs with reduced phospho-histone 3 (pH3+) cells confirming decreased cell proliferation at day 35. Pharmacological inhibition of caspase-8 activity with Z-IETD-FMK decreased apoptosis in one patient model, highlighting a potential therapeutic approach. These data reveal shared pathophysiological mechanisms contributing to a microphthalmia phenotype.

Managing the apparently blind child presenting in the first year of life: A review

Severe vision impairment and blindness in childhood have a significant health burden on the child, family and society. This review article seeks to provide a structured framework for managing the apparently blind child presenting in the first year of life, starting from a comprehensive history and examination. Different investigation modalities and the increasingly important role of genetics will also be described, in addition to common causes of severe vision impairment. Crucially, a systematic approach to the blind infant is key to correct diagnoses and timely management. Incorrect diagnoses can be costly to all involved, however it is important to note that diagnoses can change with ongoing follow-up and investigations. Furthermore, the modern age of ophthalmology requires a multi-disciplinary approach and close collaboration with specialists including paediatricians, neurologists and geneticists, in addition to rehabilitation and low vision services, to ensure the best care for these vulnerable infants.

An Unusual Presentation of Novel Missense Variant in PAX6 Gene: NM_000280.4:c.341A>G, p.(Asn114Ser)

This study investigates a unique and complex eye phenotype characterized by minimal iris defects, foveal hypoplasia, optic nerve coloboma, and severe posterior segment damage. Through genetic analysis and bioinformatic tools, a specific nonsynonymous substitution, p.(Asn114Ser), within the PAX6 gene's paired domain is identified. Although this substitution is not in direct contact with DNA, its predicted stabilizing effect on the protein structure challenges the traditional understanding of PAX6 mutations, suggesting a gain-of-function mechanism. Contrary to classical loss-of-function effects, this gain-of-function hypothesis aligns with research demonstrating PAX6's dosage sensitivity. Gain-of-function mutations, though less common, can lead to diverse phenotypes distinct from aniridia. Our findings emphasize PAX6's multifaceted influence on ocular phenotypes and the importance of genetic variations. We contribute a new perspective on PAX6 mutations by suggesting a potential gain-of-function mechanism and showcasing the complexities of ocular development. This study sheds light on the intricate interplay of the genetic alterations and regulatory mechanisms underlying complex eye phenotypes. Further research, validation, and collaboration are crucial to unravel the nuanced interactions shaping ocular health and development.

Real-world clinical and molecular management of 50 prospective patients with microphthalmia, anophthalmia and/or ocular coloboma

BACKGROUND/AIMS

Microphthalmia, anophthalmia and coloboma (MAC) are clinically and genetically heterogenous rare developmental eye conditions, which contribute to a significant proportion of childhood blindness worldwide. Clear understanding of MAC aetiology and comorbidities is essential to providing patients with appropriate care. However, current management is unstandardised and molecular diagnostic rates remain low, particularly in those with unilateral presentation. To further understanding of clinical and genetic management of patients with MAC, we charted their real-world experience to ascertain optimal management pathways and yield from molecular analysis.

METHODS

A prospective cohort study of consecutive patients with MAC referred to the ocular genetics service at Moorfields Eye Hospital between 2017-2020.

RESULTS

Clinical analysis of 50 MAC patients (15 microphthalmia; 2 anophthalmia; 11 coloboma; and 22 mixed) from 44 unrelated families found 44% had additional ocular features (complex) and 34% had systemic involvement, most frequently intellectual/developmental delay (8/17). Molecular analysis of 39 families using targeted gene panels, whole genome sequencing and microarray comparative genomic hybridisation identified genetic causes in, 28% including novel variants in six known MAC genes (SOX2, KMT2D, MAB21L2, ALDH1A3, BCOR and FOXE3), and a molecular diagnostic rate of 33% for both bilateral and unilateral cohorts. New phenotypic associations were found for FOXE3 (bilateral sensorineural hearing loss) and MAB21L2 (unilateral microphthalmia).

CONCLUSION

This study highlights the importance of thorough clinical and molecular phenotyping of MAC patients to provide appropriate multidisciplinary care. Routine genetic testing for both unilateral and bilateral cases in the clinic may increase diagnostic rates in the future, helping elucidate genotype-phenotype correlations and informing genetic counselling.

Translation of paired box 6 (PAX6) mRNA is IRES-mediated and inhibited by cymarin in breast cancer cells

Paired box 6 (PAX6) is a member of the PAX family and plays an essential role in cancer cell cycle progression, colony formation, proliferation and invasion. Its expression is upregulated in many cancers including breast cancer, but the process of PAX6 mRNA translation has rarely been studied. We found that PAX6 translation level increased in MCF-7 breast cancer cells treated with the chemotherapeutic drug adriamycin (ADM), which might be attributable to internal ribosome entry site (IRES)-mediated translation. By modifying a bicistronic luciferase plasmid that is widely used to examine IRES activity, we found that the 469-base 5'-UTR of PAX6 mRNA contains an IRES element and that core IRES activity is located between nucleotides 159 and 333. Moreover, PAX6 IRES activity was induced during ADM treatment, which may be the main reason for the elevated level of PAX6 protein. We also found that cymarin, a cardiac glycoside, acts as an inhibitor of PAX6 protein expression by impairing its IRES-mediated translation. Furthermore, MCF-7 cell proliferation was suppressed during treatment with cymarin. These results provide novel insights into the translation mechanism of PAX6 in breast cancer cells and suggest that cymarin is a promising candidate for the treatment of breast cancer via targeting the expression of PAX6.

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 missense variant expands the phenotype and genotype of PAX6-associated foveal hypoplasia accompanied by various manifestations of anterior segment dysgenesis

BACKGROUND

According to previous reports, PAX6-associated foveal hypoplasia (FH) could usually be accompanied by various anterior segment anomalies including variable iris changes. This study aims to exhibit unusual phenotypes of a novel missense variant of PAX6 from a Chinese pedigree.

METHODS

Ophthalmic examinations including slit-lamp biomicroscopy, gonioscopy, ophthalmic ultrasound, ultrasonic biomicroscopy, optical coherence tomography, wide-field fundus imaging, and visual field test were performed to evaluate the clinical manifestations. Whole-exome sequencing (WES) and bioinformatics analysis were conducted in eight members from this pedigree to identify the causative mutation.

RESULTS

WES revealed a novel heterozygous substitution of PAX6 (NM_000280.5:c.157G > A, p.(Val53Met) (chr11:31823309 C > T, hg19)), which cosegregated with the phenotype of this pedigree. All the three patients (a pair of fraternal twins and their mother) exhibited bilateral FH and anterior segment dysgenesis (ASD) including microcornea, sclerocornea, obvious symmetrical corectopia, iris stromal dysplasia, goniodysgenesis, and abnormal distribution of fundus blood vessels. The girl of the fraternal twins also demonstrated bilateral temporal deviation of lenses and abnormal tissue membrane connecting anterior chamber angle and lens anterior capsule in the right eye. The mother additionally showed apparent cataract bilaterally and cupping of the optic disc in her left eye.

CONCLUSION

A novel missense variant in PAX6 gene was detected in a Chinese pedigree demonstrating bilateral FH and ASD. It is really distinctive that the ASD involves almost all parts of the anterior segment, and bilateral symmetrical corectopia is the most perceptible sign. This study expands the phenotypic and genotypic spectrum of PAX6-associated ocular diseases, and facilitates the understanding of the crucial role that PAX6 plays in the development of the eye. Meanwhile, PAX6 could be considered as a candidate pathogenic gene of bilateral symmetrical corectopia.

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.

Report on the 2021 Aniridia North America symposium on PAX6, aniridia, and beyond

The inaugural Aniridia North America (ANA) Symposium was held on the first weekend in November 2021 in Charlottesville, VA, at the University of Virginia. The purpose of this meeting was to bring together an international group of scientists, physicians, patient advocacy groups, and individuals with aniridia to discuss recent advances in knowledge about aniridia and other congenital eye diseases and the development of potential treatments for congenital eye disorders using personalized medicine. Leaders in several areas of eye research and clinical treatment provided a broad perspective on new research advances that impact an understanding of the causes of the damage to the eye associated with aniridia and the development of novel treatments for this and related disorders. Here we summarize the research discussed at the symposium.

Monoallelic variants resulting in substitutions of MAB21L1 Arg51 Cause Aniridia and microphthalmia

Classical aniridia is a congenital and progressive panocular disorder almost exclusively caused by heterozygous loss-of-function variants at the PAX6 locus. We report nine individuals from five families with severe aniridia and/or microphthalmia (with no detectable PAX6 mutation) with ultrarare monoallelic missense variants altering the Arg51 codon of MAB21L1. These mutations occurred de novo in 3/5 families, with the remaining families being compatible with autosomal dominant inheritance. Mice engineered to carry the p.Arg51Leu change showed a highly-penetrant optic disc anomaly in heterozygous animals with severe microphthalmia in homozygotes. Substitutions of the same codon (Arg51) in MAB21L2, a close homolog of MAB21L1, cause severe ocular and skeletal malformations in humans and mice. The predicted nucleotidyltransferase function of MAB21L1 could not be demonstrated using purified protein with a variety of nucleotide substrates and oligonucleotide activators. Induced expression of GFP-tagged wildtype and mutant MAB21L1 in human cells caused only modest transcriptional changes. Mass spectrometry of immunoprecipitated protein revealed that both mutant and wildtype MAB21L1 associate with transcription factors that are known regulators of PAX6 (MEIS1, MEIS2 and PBX1) and with poly(A) RNA binding proteins. Arg51 substitutions reduce the association of wild-type MAB21L1 with TBL1XR1, a component of the NCoR complex. We found limited evidence for mutation-specific interactions with MSI2/Musashi-2, an RNA-binding proteins with effects on many different developmental pathways. Given that biallelic loss-of-function variants in MAB21L1 result in a milder eye phenotype we suggest that Arg51-altering monoallelic variants most plausibly perturb eye development via a gain-of-function mechanism.

Loss-of-function, gain-of-function and dominant-negative mutations have profoundly different effects on protein structure

Most known pathogenic mutations occur in protein-coding regions of DNA and change the way proteins are made. Taking protein structure into account has therefore provided great insight into the molecular mechanisms underlying human genetic disease. While there has been much focus on how mutations can disrupt protein structure and thus cause a loss of function (LOF), alternative mechanisms, specifically dominant-negative (DN) and gain-of-function (GOF) effects, are less understood. Here, we investigate the protein-level effects of pathogenic missense mutations associated with different molecular mechanisms. We observe striking differences between recessive vs dominant, and LOF vs non-LOF mutations, with dominant, non-LOF disease mutations having much milder effects on protein structure, and DN mutations being highly enriched at protein interfaces. We also find that nearly all computational variant effect predictors, even those based solely on sequence conservation, underperform on non-LOF mutations. However, we do show that non-LOF mutations could potentially be identified by their tendency to cluster in three-dimensional space. Overall, our work suggests that many pathogenic mutations that act via DN and GOF mechanisms are likely being missed by current variant prioritisation strategies, but that there is considerable scope to improve computational predictions through consideration of molecular disease mechanisms.

Most known pathogenic mutations occur in protein-coding regions of DNA and change the way proteins are made. Here the authors analyse the locations of thousands of human disease mutations and their predicted effects on protein structure and show that,while loss-of-function mutations tend to be highly disruptive, non-loss-of-function mutations are in general much milder at a protein structural level.

Exome sequencing identifies genetic variants in anophthalmia and microphthalmia

Anophthalmia and microphthalmia (A/M) are rare birth defects affecting up to 2 per 10,000 live births. These conditions are manifested by the absence of an eye or reduced eye volumes within the orbit leading to vision loss. Although clinical case series suggest a strong genetic component in A/M, few systematic investigations have been conducted on potential genetic contributions owing to low population prevalence. To overcome this challenge, we utilized DNA samples and data collected as part of the National Birth Defects Prevention Study (NBDPS). The NBDPS employed multi-center ascertainment of infants affected by A/M. We performed exome sequencing on 67 family trios and identified numerous genes affected by rare deleterious nonsense and missense variants in this cohort, including de novo variants. We identified 9 nonsense changes and 86 missense variants that are absent from the reference human population (Genome Aggregation Database), and we suggest that these are high priority candidate genes for A/M. We also performed literature curation, single cell transcriptome comparisons, and molecular pathway analysis on the candidate genes and performed protein structure modeling to determine the potential pathogenic variant consequences on PAX6 in this disease.

Relative Frequencies of PAX6 Mutational Events in a Russian Cohort of Aniridia Patients in Comparison with the World's Population and the Human Genome

Genome-wide sequencing metadata allows researchers to infer bias in the relative frequencies of mutational events and to predict putative mutagenic models. In addition, much less data could be useful in the evaluation of the mutational frequency spectrum and the prevalent local mutagenic process. Here we analyzed the PAX6 gene locus for mutational spectra obtained in our own and previous studies and compared them with data on other genes as well as the whole human genome. MLPA and Sanger sequencing were used for mutation searching in a cohort of 199 index patients from Russia with aniridia and aniridia-related phenotypes. The relative frequencies of different categories of PAX6 mutations were consistent with those previously reported by other researchers. The ratio between substitutions, small indels, and chromosome deletions in the 11p13 locus was within the interval previously published for 20 disease associated genomic loci, but corresponded to a higher end due to very high frequencies of small indels and chromosome deletions. The ratio between substitutions, small indels, and chromosome deletions for disease associated genes, including the PAX6 gene as well as the share of PAX6 missense mutations, differed considerably from those typical for the whole genome.

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.

Diverse Molecular Mechanisms Underlying Pathogenic Protein Mutations: Beyond the Loss-of-Function Paradigm

Most known disease-causing mutations occur in protein-coding regions of DNA. While some of these involve a loss of protein function (e.g., through premature stop codons or missense changes that destabilize protein folding), many act via alternative molecular mechanisms and have dominant-negative or gain-of-function effects. In nearly all cases, these non-loss-of-function mutations can be understood by considering interactions of the wild-type and mutant protein with other molecules, such as proteins, nucleic acids, or small ligands and substrates. Here, we review the diverse molecular mechanisms by which pathogenic mutations can have non-loss-of-function effects, including by disrupting interactions, increasing binding affinity, changing binding specificity, causing assembly-mediated dominant-negative and dominant-positive effects, creating novel interactions, and promoting aggregation and phase separation. We believe that increased awareness of these diverse molecular disease mechanisms will lead to improved diagnosis (and ultimately treatment) of human genetic disorders. 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.

The properties of human disease mutations at protein interfaces

The assembly of proteins into complexes and their interactions with other biomolecules are often vital for their biological function. While it is known that mutations at protein interfaces have a high potential to be damaging and cause human genetic disease, there has been relatively little consideration for how this varies between different types of interfaces. Here we investigate the properties of human pathogenic and putatively benign missense variants at homomeric (isologous and heterologous), heteromeric, DNA, RNA and other ligand interfaces, and at different regions in proteins with respect to those interfaces. We find that different types of interfaces vary greatly in their propensity to be associated with pathogenic mutations, with homomeric heterologous and DNA interfaces being particularly enriched in disease. We also find that residues that do not directly participate in an interface, but are close in three-dimensional space, show a significant disease enrichment. Finally, we observe that mutations at different types of interfaces tend to have distinct property changes when undergoing amino acid substitutions associated with disease, and that this is linked to substantial variability in their identification by computational variant effect predictors.

Nearly all proteins interact with other molecules as part of their biological function. For example, proteins can interact with other copies of the same type of protein, with different proteins, with DNA, or with small ligand molecules. Many mutations at protein interfaces, the regions of proteins that interact with other molecules, are known to cause human genetic disease. In this study, we first investigate how different types of protein interfaces have different tendencies to be associated with disease. We also show that the closer a mutation is to an interface, the more likely it is to cause disease. Finally, we study how mutations at different types of interfaces tend to be associated with different changes in amino acid properties, which appears to influence our ability to computationally predict the effects of mutations. Ultimately, we hope that consideration of protein interface properties will eventually improve our ability to identify new disease-causing mutations.

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.

SAMPDI-3D: predicting the effects of protein and DNA mutations on protein-DNA interactions

MOTIVATION

Mutations that alter protein-DNA interactions may be pathogenic and cause diseases. Therefore, it is extremely important to quantify the effect of mutations on protein-DNA binding free energy to reveal the molecular origin of diseases and to assist the development of treatments. Although several methods that predict the change of protein-DNA binding affinity upon mutations in the binding protein were developed, the effect of DNA mutations was not considered yet.

RESULTS

Here, we report a new version of SAMPDI, the SAMPDI-3D, which is a gradient boosting decision tree machine learning method to predict the change of the protein-DNA binding free energy caused by mutations in both the binding protein and the bases of the corresponding DNA. The method is shown to achieve Pearson correlation coefficient of 0.76 and 0.80 in a benchmarking test against experimentally determined change of the binding free energy caused by mutations in the binding protein or DNA, respectively. Furthermore, three datasets collected from literature were used to do blind benchmark for SAMPDI-3D and it is shown that it outperforms all existing state-of-the-art methods. The method is very fast allowing for genome-scale investigations.

AVAILABILITY

It is available as a web server and a stand-code at http://compbio.clemson.edu/SAMPDI-3D/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genotype-Phenotype of Isolated Foveal Hypoplasia in a Large Cohort: Minor Iris Changes as an Indicator of PAX6 Involvement

Purpose

Foveal hypoplasia (FVH) is defined as the lack of fovea with a relatively preserved neuroretina, occurring either as an isolated FVH (IFVH) condition or associated with other diseases. This study aimed to systemically molecularly characterize IFVH.

Methods

Genetic defects in 33 families with IFVH were analyzed by exome sequencing. Variants in three genes (PAX6, SLC38A8, and AHR) were selected and evaluated with multistep bioinformatic tools.

Results

Mutations in the three genes were identified in 69.7% (23/33) of families with IFVH and infantile nystagmus, including 18 families with PAX6 mutations, 5 with SLC38A8 mutations, but none with AHR mutations. Clinical data from 32 patients in the 23 families showed FVH, infantile nystagmus, and full iris. Careful follow-up visits revealed subtle changes in iris in 9 of 14 patients with PAX6 variants. The PAX6 variants of the 18 families (15 missense and one stop-loss) were mostly located in the C-terminal region of the paired box domain. Variants in AHR, SLC38A8, and PAX6 contributed to IFVH in one (2%), 25 (45%), and 30 (53%) families with identified genetic defects (23 families in this study and 33 reported previously), respectively.

Conclusions

PAX6 and SLC38A8 mutations are the main cause of IFVH based on our data and a systematic review. IFVH-associated PAX6 variants are mostly missense with a specific location, indicating a specific correlation of these variants with IFVH but not with typical aniridia. Full iris with subtle structural abnormalities is more common in patients with PAX6-associated IFVH, suggesting a potential diagnostic indicator.

Longitudinal genotype-phenotype analysis in 86 patients with PAX6-related aniridia

Aniridia is most commonly caused by haploinsufficiency of the PAX6 gene, characterized by variable iris and foveal hypoplasia, nystagmus, cataracts, glaucoma, and aniridia-related keratopathy (ARK). Genotype-phenotype correlations have previously been described; however, detailed longitudinal studies of aniridia are less commonly reported. We identified 86 patients from 62 unrelated families with molecularly confirmed heterozygous PAX6 variants from a UK-based single-center ocular genetics service. They were categorized into mutation groups, and a retrospective review of clinical characteristics (ocular and systemic) from baseline to most recent was recorded. One hundred and seventy-two eyes were evaluated, with a mean follow-up period of 16.3 ± 12.7 years. Nystagmus was recorded in 87.2% of the eyes, and foveal hypoplasia was found in 75%. Cataracts were diagnosed in 70.3%, glaucoma in 20.6%, and ARK in 68.6% of eyes. Prevalence, age of diagnosis and surgical intervention, and need for surgical intervention varied among mutation groups. Overall, the missense mutation subgroup had the mildest phenotype, and surgically naive eyes maintained better visual acuity. Systemic evaluation identified type 2 diabetes in 12.8% of the study group, which is twice the UK prevalence. This is the largest longitudinal study of aniridia in the UK, and as such, it can provide insights into prognostic indicators for patients and guiding clinical management of both ocular and systemic features.

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.

Genetic Analysis in a Swiss Cohort of Bilateral Congenital Cataract

Importance

Identification of geographic population-based differences in genotype and phenotype heterogeneity are important for targeted and patient-specific diagnosis and treatment, counseling, and screening strategies.

Objective

To report disease-causing variants and their detailed phenotype in patients with bilateral congenital cataract from a single center in Switzerland and thereby draw a genetic map and perform a genotype-phenotype comparison of this cohort.

Design, Setting, and Participants

This clinical and molecular-genetic cohort study took place through the collaboration of the Department of Ophthalmology at the University Hospital Zurich and the Institute of Medical Molecular Genetics, University of Zurich, Schlieren, Switzerland. Thirty-seven patients from 25 families with different types of bilateral congenital cataract were included. All participating family members received a comprehensive eye examination. Whole exome sequencing was performed in the index patients, followed by a filtering process to detect possible disease-associated variants in genes previously described in association with congenital cataract. Probable disease-causing variants were confirmed by Sanger sequencing in available family members. All data were collected from January 2018 to June 2020, and the molecular-genetic analyses were performed from January 2019 to July 2020.

Main Outcomes and Measures

Identification of the underlying genetic causes of bilateral congenital cataract, including novel disease-causing variants and phenotype correlation.

Results

Among the 37 patients (18 [49%] male and 19 [51%] female; mean [SD] age, 17.3 [15.9] years) from 25 families, pathogenic variants were detected in 20 families (80% detection rate), which included 13 novel variants in the following genes: BCOR, COL4A1, CRYBA2, CRYBB2, CRYGC, CRYGS, GJA3, MAF, NHS, and WFS1. Putative disease-causing variants were identified in 14 of 20 families (70%) as isolated cases and in 6 of 20 families (30%) with syndromic cases. A recessive variant in the CRYBB2 gene in a consanguineous family with 2 affected siblings showing a nuclear and sutural cataract was reported in contrast to previously published reports. In addition, the effect on splicing in a minigene assay of a novel splice site variant in the NHS gene (c.[719-2A>G]) supported the pathogenicity of this variant.

Conclusions and Relevance

This study emphasizes the importance of genetic testing of congenital cataracts. Known dominant genes need to be considered for recessive inheritance patterns. Syndromic types of cataract may be underdiagnosed in patients with mild systemic features.

Expanding the Phenotypic Spectrum of PAX6 Mutations: From Congenital Cataracts to Nystagmus

Background: Congenital aniridia is a complex ocular disorder, usually associated with severe visual impairment, generally caused by mutations on the PAX6 gene. The clinical phenotype of PAX6 mutations is highly variable, making the genotype–phenotype correlations difficult to establish. Methods: we describe the phenotype of eight patients from seven unrelated families with confirmed mutations in PAX6, and very different clinical manifestations. Results: Only two patients had the classical aniridia phenotype while the other two presented with aniridia-related manifestations, such as aniridia-related keratopathy or partial aniridia. Congenital cataracts were the main manifestation in three of the patients in this series. All the patients had nystagmus and low visual acuity. Conclusions: The diagnosis of mild forms of aniridia is challenging, but these patients have a potentially blinding hereditary disease that might present with a more severe phenotype in future generations. Clinicians should be aware of the mild aniridia phenotype and request genetic testing to perform an accurate diagnosis.

Activation of cryptic donor splice sites by non-coding and coding PAX6 variants contributes to congenital aniridia

BACKGROUND

The paired-domain transcription factor paired box gene 6 (PAX6) causes a wide spectrum of ocular developmental anomalies, including congenital aniridia, Peters anomaly and microphthalmia. Here, we aimed to functionally assess the involvement of seven potentially non-canonical splicing variants on missplicing of exon 6, which represents the main hotspot region for loss-of-function PAX6 variants.

METHODS

By locus-specific analysis of PAX6 using Sanger and/or targeted next-generation sequencing, we screened a Spanish cohort of 106 patients with PAX6-related diseases. Functional splicing assays were performed by in vitro minigene approaches or directly in RNA from patient-derived lymphocytes cell line, when available.

RESULTS

Five out seven variants, including three synonymous changes, one small exonic deletion and one non-canonical splice variant, showed anomalous splicing patterns yielding partial exon skipping and/or elongation.

CONCLUSION

We describe new spliceogenic mechanisms for PAX6 variants mediated by creating or strengthening five different cryptic donor sites at exon 6. Our work revealed that the activation of cryptic PAX6 splicing sites seems to be a recurrent and underestimated cause of aniridia. Our findings pointed out the importance of functional assessment of apparently silent PAX6 variants to uncover hidden genetic alterations and to improve variant interpretation for genetic counselling in aniridia.

Structural and functional consequences of PAX6 mutations in the brain: Implications for aniridia

The paired-box 6 (PAX6) gene encodes a highly conserved transcription factor essential for the proper development of the eye and brain. Heterozygous loss-of-function mutations in PAX6 are causal for a condition known as aniridia in humans and the Small eye phenotype in mice. Aniridia is characterized by iris hypoplasia and other ocular abnormalities, but recent evidence of neuroanatomical, sensory, and cognitive impairments in this population has emerged, indicating brain-related phenotypes as a prevalent feature of the disorder. Determining the neurophysiological origins of brain-related phenotypes in this disorder presents a substantial challenge, as the majority of extra-ocular traits in aniridia demonstrate a high degree of heterogeneity. Here, we summarize and integrate findings from human and rodent model studies, which have focused on neuroanatomical and functional consequences of PAX6 mutations. We highlight novel findings from PAX6 central nervous system studies in adult mammals, and integrate these findings into what we know about PAX6's role in development of the central nervous system. This review presents the current literature in the field in order to inform clinical application, discusses what is needed in future studies, and highlights PAX6 as a lens through which to understand genetic disorders affecting the human nervous system.

Generation of human iPSC line (UCLi013-A) from a patient with microphthalmia and aniridia, carrying a heterozygous missense mutation c.372C>A p.(Asn124Lys) in PAX6

A human induced pluripotent stem cell (hiPSC) line (UCLi013-A) was generated from fibroblast cells of a 34-year-old donor with multiple ocular conditions including severe microphthalmia and aniridia. The patient had a heterozygous missense mutation in PAX6 c.372C>A, p.(Asn124Lys), validated in the fibroblasts through Sanger sequencing. Fibroblasts derived from a skin biopsy were reprogrammed using integration free episomal reprogramming. The established iPSC line was found to express pluripotency markers, exhibit differentiation potential in vitro and display a normal karyotype. This cell line will act as a tool for disease modelling of microphthalmia and aniridia, identification of therapeutic targets and drug screening.

Animal and cellular models of microphthalmia

Microphthalmia is a rare developmental eye disorder affecting 1 in 7000 births. It is defined as a small (axial length ⩾2 standard deviations below the age-adjusted mean) underdeveloped eye, caused by disruption of ocular development through genetic or environmental factors in the first trimester of pregnancy. Clinical phenotypic heterogeneity exists amongst patients with varying levels of severity, and associated ocular and systemic features. Up to 11% of blind children are reported to have microphthalmia, yet currently no treatments are available. By identifying the aetiology of microphthalmia and understanding how the mechanisms of eye development are disrupted, we can gain a better understanding of the pathogenesis. Animal models, mainly mouse, zebrafish and Xenopus, have provided extensive information on the genetic regulation of oculogenesis, and how perturbation of these pathways leads to microphthalmia. However, differences exist between species, hence cellular models, such as patient-derived induced pluripotent stem cell (iPSC) optic vesicles, are now being used to provide greater insights into the human disease process. Progress in 3D cellular modelling techniques has enhanced the ability of researchers to study interactions of different cell types during eye development. Through improved molecular knowledge of microphthalmia, preventative or postnatal therapies may be developed, together with establishing genotype-phenotype correlations in order to provide patients with the appropriate prognosis, multidisciplinary care and informed genetic counselling. This review summarises some key discoveries from animal and cellular models of microphthalmia and discusses how innovative new models can be used to further our understanding in the future.

Plain language summary

Animal and Cellular Models of the Eye Disorder, Microphthalmia (Small Eye) Microphthalmia, meaning a small, underdeveloped eye, is a rare disorder that children are born with. Genetic changes or variations in the environment during the first 3 months of pregnancy can disrupt early development of the eye, resulting in microphthalmia. Up to 11% of blind children have microphthalmia, yet currently no treatments are available. By understanding the genes necessary for eye development, we can determine how disruption by genetic changes or environmental factors can cause this condition. This helps us understand why microphthalmia occurs, and ensure patients are provided with the appropriate clinical care and genetic counselling advice. Additionally, by understanding the causes of microphthalmia, researchers can develop treatments to prevent or reduce the severity of this condition. Animal models, particularly mice, zebrafish and frogs, which can also develop small eyes due to the same genetic/environmental changes, have helped us understand the genes which are important for eye development and can cause birth eye defects when disrupted. Studying a patient's own cells grown in the laboratory can further help researchers understand how changes in genes affect their function. Both animal and cellular models can be used to develop and test new drugs, which could provide treatment options for patients living with microphthalmia. This review summarises the key discoveries from animal and cellular models of microphthalmia and discusses how innovative new models can be used to further our understanding in the future.

Morphometric analysis of the lens in human aniridia and mouse Small eye

Congenital aniridia is caused by heterozygous mutations in the PAX6 gene. In this disease, congenital iris and foveal hypoplasia is associated with juvenile onset cataract, glaucoma, and corneal keratopathy. In rodents, Pax6 mutations result in a congenital reduction in ocular size that is not typically described in human aniridia. Here, the ocular morphometry of aniridia patients is compared with the lens phenotype of Pax6^+/tm1/Pgr mice to reveal whether there are species differences in Pax6 regulation of lens development and homeostasis. Ultrasound biometry (UBM) revealed that eleven percent of aniridia patients exhibited mild microphthalmia while the anterior chamber depth of aniridic eyes was significantly reduced from 6 months of age onward. Although aniridic lens thickness was normal from birth, it was significantly decreased in aniridic lenses older than 30. Notably, 86% of aniridic lenses exhibited cataractous changes in this cohort. In addition, a significant proportion of aniridia patients develop lens subluxation as they age associated with reduced lens diameter as measured by anterior segment optical coherence tomography (AS-OCT). Analysis of young adult Pax6^+/tm1/Pgr mouse lenses by micro-computed tomography (microCT), bright field and dark field imaging revealed that they are reduced in size but did not exhibit overt cataracts at this age. Overall, this study reveals that congenital microphthalmia as assessed by axial length, or microphakia, as assessed by lens thickness, are not typical in human aniridia, although these are primary manifestations of Pax6 mutations in mice, suggesting that PAX6 regulates some aspects of lens development differently between these species.

PAX6 missense variants in two families with isolated foveal hypoplasia and nystagmus: evidence of paternal postzygotic mosaicism

PAX6 is considered the master regulator of eye development, the majority of variants affecting this gene cause the pan-ocular developmental eye disorder aniridia. Although no genotype-phenotype correlations are clearly established, missense variants affecting the DNA-binding paired domain of PAX6 are usually associated with non-aniridia phenotypes like microphthalmia, coloboma or isolated foveal hypoplasia. In this study, we report two missense heterozygous variants in the paired domain of PAX6 resulting in isolated foveal hypoplasia with nystagmus in two independent families: c.112 C > G; p.(Arg38Gly) and c.214 G > C; p.(Gly72Arg) in exons 5 and 6, respectively. Furthermore, we provide evidence that paternal postzygotic mosaicism is the cause of inheritance, with clinically unaffected fathers and reduced affected allele fraction. This work contributes to increase the phenotypic spectrum caused by PAX6 variants, and to our knowledge, is the first report to describe the presence of postzygotic parental mosaicism causing isolated foveal hypoplasia with nystagmus. These results support the growing evidence that suggest an overestimation of sporadic cases with PAX6 variants, which has strong implications for both genetic counselling and family planning.

Identification of pathogenic missense mutations using protein stability predictors

Attempts at using protein structures to identify disease-causing mutations have been dominated by the idea that most pathogenic mutations are disruptive at a structural level. Therefore, computational stability predictors, which assess whether a mutation is likely to be stabilising or destabilising to protein structure, have been commonly used when evaluating new candidate disease variants, despite not having been developed specifically for this purpose. We therefore tested 13 different stability predictors for their ability to discriminate between pathogenic and putatively benign missense variants. We find that one method, FoldX, significantly outperforms all other predictors in the identification of disease variants. Moreover, we demonstrate that employing predicted absolute energy change scores improves performance of nearly all predictors in distinguishing pathogenic from benign variants. Importantly, however, we observe that the utility of computational stability predictors is highly heterogeneous across different proteins, and that they are all inferior to the best performing variant effect predictors for identifying pathogenic mutations. We suggest that this is largely due to alternate molecular mechanisms other than protein destabilisation underlying many pathogenic mutations. Thus, better ways of incorporating protein structural information and molecular mechanisms into computational variant effect predictors will be required for improved disease variant prioritisation.

Systemic and ocular manifestations of a patient with mosaic ARID1A-associated Coffin-Siris syndrome and review of select mosaic conditions with ophthalmic manifestations

Mosaic genetic mutations may be somatic, germline, or "gonosomal" and have the potential to cause genetic syndromes, disorders, or malformations. Mutations can occur at any point in embryonic development and the timing determines the extent of distribution of the mutation throughout the body and different tissue types. The eye and visual pathway offer a unique opportunity to study somatic and gonosomal mosaic mutations as the eye consists of tissues derived from all three germ layers allowing disease pathology to be assessed with noninvasive imaging. In this review, we describe systemic and ocular manifestations in a child with mosaic Coffin-Siris syndrome. The patient presented with a significant medical history of accommodative esotropia and hyperopia, macrocephaly, polydactyly, global developmental delay, hypotonia, ureteropelvic junction (UPJ) obstruction, and brain MRI abnormalities. The ophthalmic findings in this patient were nonspecific, however, they are consistent with ocular manifestations reported in other patients with Coffin-Siris syndrome. We also review ophthalmic findings of select mosaic chromosomal and single-gene disorders. Ophthalmic assessment alongside clinical genetic testing may play an important role in diagnosis of genetic syndromes as well as understanding disease pathology, particularly when mosaicism plays a role.

Using deep mutational scanning to benchmark variant effect predictors and identify disease mutations

To deal with the huge number of novel protein-coding variants identified by genome and exome sequencing studies, many computational variant effect predictors (VEPs) have been developed. Such predictors are often trained and evaluated using different variant data sets, making a direct comparison between VEPs difficult. In this study, we use 31 previously published deep mutational scanning (DMS) experiments, which provide quantitative, independent phenotypic measurements for large numbers of single amino acid substitutions, in order to benchmark and compare 46 different VEPs. We also evaluate the ability of DMS measurements and VEPs to discriminate between pathogenic and benign missense variants. We find that DMS experiments tend to be superior to the top-ranking predictors, demonstrating the tremendous potential of DMS for identifying novel human disease mutations. Among the VEPs, DeepSequence clearly stood out, showing both the strongest correlations with DMS data and having the best ability to predict pathogenic mutations, which is especially remarkable given that it is an unsupervised method. We further recommend SNAP2, DEOGEN2, SNPs&GO, SuSPect and REVEL based upon their performance in these analyses.

Inherited Eye Diseases with Retinal Manifestations through the Eyes of Homeobox Genes

Retinal development is under the coordinated control of overlapping networks of signaling pathways and transcription factors. The paper was conceived as a review of the data and ideas that have been formed to date on homeobox genes mutations that lead to the disruption of eye organogenesis and result in inherited eye/retinal diseases. Many of these diseases are part of the same clinical spectrum and have high genetic heterogeneity with already identified associated genes. We summarize the known key regulators of eye development, with a focus on the homeobox genes associated with monogenic eye diseases showing retinal manifestations. Recent advances in the field of genetics and high-throughput next-generation sequencing technologies, including single-cell transcriptome analysis have allowed for deepening of knowledge of the genetic basis of inherited retinal diseases (IRDs), as well as improve their diagnostics. We highlight some promising avenues of research involving molecular-genetic and cell-technology approaches that can be effective for IRDs therapy. The most promising neuroprotective strategies are aimed at mobilizing the endogenous cellular reserve of the retina.

The Spectrum of PAX6 Mutations and Genotype-Phenotype Correlations in the Eye

The transcription factor PAX6 is essential in ocular development in vertebrates, being considered the master regulator of the eye. During eye development, it is essential for the correct patterning and formation of the multi-layered optic cup and it is involved in the developing lens and corneal epithelium. In adulthood, it is mostly expressed in cornea, iris, and lens. PAX6 is a dosage-sensitive gene and it is highly regulated by several elements located upstream, downstream, and within the gene. There are more than 500 different mutations described to affect PAX6 and its regulatory regions, the majority of which lead to PAX6 haploinsufficiency, causing several ocular and systemic abnormalities. Aniridia is an autosomal dominant disorder that is marked by the complete or partial absence of the iris, foveal hypoplasia, and nystagmus, and is caused by heterozygous PAX6 mutations. Other ocular abnormalities have also been associated with PAX6 changes, and genotype-phenotype correlations are emerging. This review will cover recent advancements in PAX6 regulation, particularly the role of several enhancers that are known to regulate PAX6 during eye development and disease. We will also present an updated overview of the mutation spectrum, where an increasing number of mutations in the non-coding regions have been reported. Novel genotype-phenotype correlations will also be discussed.