Phenotype-aware prioritisation of rare Mendelian disease variants.

A molecular diagnosis from the analysis of sequencing data in rare Mendelian diseases has a huge impact on the management of patients and their families. Numerous patient phenotype-aware variant prioritisation (VP) tools have been developed to help automate this process, and shorten the diagnostic odyssey, but performance statistics on real patient data are limited. Here we identify, assess, and compare the performance of all up-to-date, freely available, and programmatically accessible tools using a whole-exome, retinal disease dataset from 134 individuals with a molecular diagnosis. All tools were able to identify around two-thirds of the genetic diagnoses as the top-ranked candidate, with LIRICAL performing best overall. Finally, we discuss the challenges to overcome most cases remaining undiagnosed after current, state-of-the-art practices.

De novo and inherited monoallelic variants in TUBA4A cause ataxia and spasticity.

Alpha-tubulin 4A encoding gene (TUBA4A) has been associated with familial amyotrophic lateral sclerosis (fALS) and fronto-temporal dementia (FTD), based on identification of likely pathogenic variants in patients from distinct ALS and FTD cohorts. By screening a multicentric French cohort of 448 unrelated probands presenting with cerebellar ataxia, we identified ultra-rare TUBA4A missense variants, all being absent from public databases and predicted pathogenic by multiple in-silico tools. In addition, gene burden analyses in the 100,000 genomes project (100KGP) showed enrichment of TUBA4A rare variants in the inherited ataxia group compared to controls (OR: 57.0847 [10.2- 576.7]; p = 4.02 x10-07). Altogether, we report 12 patients presenting with spasticity and/or cerebellar ataxia and harboring a predicted pathogenic TUBA4A missense mutation, including 5 confirmed de novo cases and a mutation previously reported in a large family presenting with spastic ataxia. Cultured fibroblasts from 3 patients harboring distinct TUBA4A missense showed significant alterations in microtubule organisation and dynamics, providing insight of TUBA4A variants pathogenicity. Our data confirm the identification of a hereditary spastic ataxia disease gene with variable age of onset, expanding the clinical spectrum of TUBA4A associated phenotypes.

Beyond factor H: The impact of genetic-risk variants for age-related macular degeneration on circulating factor-H-like 1 and factor-H-related protein concentrations.

Age-related macular degeneration (AMD) is a leading cause of vision loss; there is strong genetic susceptibility at the complement factor H (CFH) locus. This locus encodes a series of complement regulators: factor H (FH), a splice variant factor-H-like 1 (FHL-1), and five factor-H-related proteins (FHR-1 to FHR-5), all involved in the regulation of complement factor C3b turnover. Little is known about how AMD-associated variants at this locus might influence FHL-1 and FHR protein concentrations. We have used a bespoke targeted mass-spectrometry assay to measure the circulating concentrations of all seven complement regulators and demonstrated elevated concentrations in 352 advanced AMD-affected individuals for all FHR proteins (FHR-1, p = 2.4 × 10-10; FHR-2, p = 6.0 × 10-10; FHR-3, p = 1.5 × 10-5; FHR-4, p = 1.3 × 10-3; FHR-5, p = 1.9 × 10-4) and FHL-1 (p = 4.9 × 10-4) when these individuals were compared to 252 controls, whereas no difference was seen for FH (p = 0.94). Genome-wide association analyses in controls revealed genome-wide-significant signals at the CFH locus for all five FHR proteins, and univariate Mendelian-randomization analyses strongly supported the association of FHR-1, FHR-2, FHR-4, and FHR-5 with AMD susceptibility. These findings provide a strong biochemical explanation for how genetically driven alterations in circulating FHR proteins could be major drivers of AMD and highlight the need for research into FHR protein modulation as a viable therapeutic avenue for AMD.

100,000 Genomes Pilot on Rare-Disease Diagnosis in Health Care - Preliminary Report.

BACKGROUND: The U.K. 100,000 Genomes Project is in the process of investigating the role of genome sequencing in patients with undiagnosed rare diseases after usual care and the alignment of this research with health care implementation in the U.K. National Health Service. Other parts of this project focus on patients with cancer and infection. METHODS: We conducted a pilot study involving 4660 participants from 2183 families, among whom 161 disorders covering a broad spectrum of rare diseases were present. We collected data on clinical features with the use of Human Phenotype Ontology terms, undertook genome sequencing, applied automated variant prioritization on the basis of applied virtual gene panels and phenotypes, and identified novel pathogenic variants through research analysis. RESULTS: Diagnostic yields varied among family structures and were highest in family trios (both parents and a proband) and families with larger pedigrees. Diagnostic yields were much higher for disorders likely to have a monogenic cause (35%) than for disorders likely to have a complex cause (11%). Diagnostic yields for intellectual disability, hearing disorders, and vision disorders ranged from 40 to 55%. We made genetic diagnoses in 25% of the probands. A total of 14% of the diagnoses were made by means of the combination of research and automated approaches, which was critical for cases in which we found etiologic noncoding, structural, and mitochondrial genome variants and coding variants poorly covered by exome sequencing. Cohortwide burden testing across 57,000 genomes enabled the discovery of three new disease genes and 19 new associations. Of the genetic diagnoses that we made, 25% had immediate ramifications for clinical decision making for the patients or their relatives. CONCLUSIONS: Our pilot study of genome sequencing in a national health care system showed an increase in diagnostic yield across a range of rare diseases. (Funded by the National Institute for Health Research and others.).

Evaluation of phenotype-driven gene prioritization methods for Mendelian diseases.

Yuan et al. recently described an independent evaluation of several phenotype-driven gene prioritization methods for Mendelian disease on two separate, clinical datasets. Although they attempted to use default settings for each tool, we describe three key differences from those we currently recommend for our Exomiser and PhenIX tools. These influence how variant frequency, quality and predicted pathogenicity are used for filtering and prioritization. We propose that these differences account for much of the discrepancy in performance between that reported by them (15-26% diagnoses ranked top by Exomiser) and previously published reports by us and others (72-77%). On a set of 161 singleton samples, we show using these settings increases performance from 34% to 72% and suggest a reassessment of Exomiser and PhenIX on their datasets using these would show a similar uplift.

Functional genomics provide key insights to improve the diagnostic yield of hereditary ataxia.

Improvements in functional genomic annotation have led to a critical mass of neurogenetic discoveries. This is exemplified in hereditary ataxia, a heterogeneous group of disorders characterised by incoordination from cerebellar dysfunction. Associated pathogenic variants in more than 300 genes have been described, leading to a detailed genetic classification partitioned by age-of-onset. Despite these advances, up to 75% of patients with ataxia remain molecularly undiagnosed even following whole genome sequencing, as exemplified in the 100 000 Genomes Project. This study aimed to understand whether we can improve our knowledge of the genetic architecture of hereditary ataxia by leveraging functional genomic annotations, and as a result, generate insights and strategies that raise the diagnostic yield. To achieve these aims, we used publicly-available multi-omics data to generate 294 genic features, capturing information relating to a gene's structure, genetic variation, tissue-specific, cell-type-specific and temporal expression, as well as protein products of a gene. We studied these features across genes typically causing childhood-onset, adult-onset or both types of disease first individually, then collectively. This led to the generation of testable hypotheses which we investigated using whole genome sequencing data from up to 2182 individuals presenting with ataxia and 6658 non-neurological probands recruited in the 100 000 Genomes Project. Using this approach, we demonstrated a high short tandem repeat (STR) density within childhood-onset genes suggesting that we may be missing pathogenic repeat expansions within this cohort. This was verified in both childhood- and adult-onset ataxia patients from the 100 000 Genomes Project who were unexpectedly found to have a trend for higher repeat sizes even at naturally-occurring STRs within known ataxia genes, implying a role for STRs in pathogenesis. Using unsupervised analysis, we found significant similarities in genomic annotation across the gene panels, which suggested adult- and childhood-onset patients should be screened using a common diagnostic gene set. We tested this within the 100 000 Genomes Project by assessing the burden of pathogenic variants among childhood-onset genes in adult-onset patients and vice versa. This demonstrated a significantly higher burden of rare, potentially pathogenic variants in conventional childhood-onset genes among individuals with adult-onset ataxia. Our analysis has implications for the current clinical practice in genetic testing for hereditary ataxia. We suggest that the diagnostic rate for hereditary ataxia could be increased by removing the age-of-onset partition, and through a modified screening for repeat expansions in naturally-occurring STRs within known ataxia-associated genes, in effect treating these regions as candidate pathogenic loci.

Missense variants in RPH3A cause defects in excitatory synaptic function and are associated with a clinically variable neurodevelopmental disorder.

PURPOSE: RPH3A encodes a protein involved in the stabilization of GluN2A subunit of N-methyl-D-aspartate (NMDA)-type glutamate receptors at the cell surface, forming a complex essential for synaptic plasticity and cognition. We investigated the effect of variants in RPH3A in patients with neurodevelopmental disorders. METHODS: By using trio-based exome sequencing, GeneMatcher, and screening of 100,000 Genomes Project data, we identified 6 heterozygous variants in RPH3A. In silico and in vitro models, including rat hippocampal neuronal cultures, have been used to characterize the effect of the variants. RESULTS: Four cases had a neurodevelopmental disorder with untreatable epileptic seizures [p.(Gln73His)dn; p.(Arg209Lys); p.(Thr450Ser)dn; p.(Gln508His)], and 2 cases [p.(Arg235Ser); p.(Asn618Ser)dn] showed high-functioning autism spectrum disorder. Using neuronal cultures, we demonstrated that p.(Thr450Ser) and p.(Asn618Ser) reduce the synaptic localization of GluN2A; p.(Thr450Ser) also increased the surface levels of GluN2A. Electrophysiological recordings showed increased GluN2A-dependent NMDA ionotropic glutamate receptor currents for both variants and alteration of postsynaptic calcium levels. Finally, expression of the Rph3AThr450Ser variant in neurons affected dendritic spine morphology. CONCLUSION: Overall, we provide evidence that missense gain-of-function variants in RPH3A increase GluN2A-containing NMDA ionotropic glutamate receptors at extrasynaptic sites, altering synaptic function and leading to a clinically variable neurodevelopmental presentation ranging from untreatable epilepsy to autism spectrum disorder.

Expanding SPTAN1 monoallelic variant associated disorders: From epileptic encephalopathy to pure spastic paraplegia and ataxia.

PURPOSE: Nonerythrocytic αII-spectrin (SPTAN1) variants have been previously associated with intellectual disability and epilepsy. We conducted this study to delineate the phenotypic spectrum of SPTAN1 variants. METHODS: We carried out SPTAN1 gene enrichment analysis in the rare disease component of the 100,000 Genomes Project and screened 100,000 Genomes Project, DECIPHER database, and GeneMatcher to identify individuals with SPTAN1 variants. Functional studies were performed on fibroblasts from 2 patients. RESULTS: Statistically significant enrichment of rare (minor allele frequency < 1 × 10-5) probably damaging SPTAN1 variants was identified in families with hereditary ataxia (HA) or hereditary spastic paraplegia (HSP) (12/1142 cases vs 52/23,847 controls, p = 2.8 × 10-5). We identified 31 individuals carrying SPTAN1 heterozygous variants or deletions. A total of 10 patients presented with pure or complex HSP/HA. The remaining 21 patients had developmental delay and seizures. Irregular αII-spectrin aggregation was noted in fibroblasts derived from 2 patients with p.(Arg19Trp) and p.(Glu2207del) variants. CONCLUSION: We found that SPTAN1 is a genetic cause of neurodevelopmental disorder, which we classified into 3 distinct subgroups. The first comprises developmental epileptic encephalopathy. The second group exhibits milder phenotypes of developmental delay with or without seizures. The final group accounts for patients with pure or complex HSP/HA.

Phenotype-driven approaches to enhance variant prioritization and diagnosis of rare disease.

Rare disease diagnostics and disease gene discovery have been revolutionized by whole-exome and genome sequencing but identifying the causative variant(s) from the millions in each individual remains challenging. The use of deep phenotyping of patients and reference genotype-phenotype knowledge, alongside variant data such as allele frequency, segregation, and predicted pathogenicity, has proved an effective strategy to tackle this issue. Here we review the numerous tools that have been developed to automate this approach and demonstrate the power of such an approach on several thousand diagnosed cases from the 100,000 Genomes Project. Finally, we discuss the challenges that need to be overcome if we are going to improve detection rates and help the majority of patients that still remain without a molecular diagnosis after state-of-the-art genomic interpretation.

Heterozygous UCHL1 loss-of-function variants cause a neurodegenerative disorder with spasticity, ataxia, neuropathy, and optic atrophy.

PURPOSE: Biallelic variants in UCHL1 have been associated with a progressive early-onset neurodegenerative disorder, autosomal recessive spastic paraplegia type 79. In this study, we investigated heterozygous UCHL1 variants on the basis of results from cohort-based burden analyses. METHODS: Gene-burden analyses were performed on exome and genome data of independent cohorts of patients with hereditary ataxia and spastic paraplegia from Germany and the United Kingdom in a total of 3169 patients and 33,141 controls. Clinical data of affected individuals and additional independent families were collected and evaluated. Patients' fibroblasts were used to perform mass spectrometry-based proteomics. RESULTS: UCHL1 was prioritized in both independent cohorts as a candidate gene for an autosomal dominant disorder. We identified a total of 34 cases from 18 unrelated families, carrying 13 heterozygous loss-of-function variants (15 families) and an inframe insertion (3 families). Affected individuals mainly presented with spasticity (24/31), ataxia (28/31), neuropathy (11/21), and optic atrophy (9/17). The mass spectrometry-based proteomics showed approximately 50% reduction of UCHL1 expression in patients' fibroblasts. CONCLUSION: Our bioinformatic analysis, in-depth clinical and genetic workup, and functional studies established haploinsufficiency of UCHL1 as a novel disease mechanism in spastic ataxia.

The Monarch Initiative in 2019: an integrative data and analytic platform connecting phenotypes to genotypes across species.

In biology and biomedicine, relating phenotypic outcomes with genetic variation and environmental factors remains a challenge: patient phenotypes may not match known diseases, candidate variants may be in genes that haven't been characterized, research organisms may not recapitulate human or veterinary diseases, environmental factors affecting disease outcomes are unknown or undocumented, and many resources must be queried to find potentially significant phenotypic associations. The Monarch Initiative (https://monarchinitiative.org) integrates information on genes, variants, genotypes, phenotypes and diseases in a variety of species, and allows powerful ontology-based search. We develop many widely adopted ontologies that together enable sophisticated computational analysis, mechanistic discovery and diagnostics of Mendelian diseases. Our algorithms and tools are widely used to identify animal models of human disease through phenotypic similarity, for differential diagnostics and to facilitate translational research. Launched in 2015, Monarch has grown with regards to data (new organisms, more sources, better modeling); new API and standards; ontologies (new Mondo unified disease ontology, improvements to ontologies such as HPO and uPheno); user interface (a redesigned website); and community development. Monarch data, algorithms and tools are being used and extended by resources such as GA4GH and NCATS Translator, among others, to aid mechanistic discovery and diagnostics.

Expansion of the Human Phenotype Ontology (HPO) knowledge base and resources.

The Human Phenotype Ontology (HPO)-a standardized vocabulary of phenotypic abnormalities associated with 7000+ diseases-is used by thousands of researchers, clinicians, informaticians and electronic health record systems around the world. Its detailed descriptions of clinical abnormalities and computable disease definitions have made HPO the de facto standard for deep phenotyping in the field of rare disease. The HPO's interoperability with other ontologies has enabled it to be used to improve diagnostic accuracy by incorporating model organism data. It also plays a key role in the popular Exomiser tool, which identifies potential disease-causing variants from whole-exome or whole-genome sequencing data. Since the HPO was first introduced in 2008, its users have become both more numerous and more diverse. To meet these emerging needs, the project has added new content, language translations, mappings and computational tooling, as well as integrations with external community data. The HPO continues to collaborate with clinical adopters to improve specific areas of the ontology and extend standardized disease descriptions. The newly redesigned HPO website (www.human-phenotype-ontology.org) simplifies browsing terms and exploring clinical features, diseases, and human genes.

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

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

The Human Phenotype Ontology in 2017.

Deep phenotyping has been defined as the precise and comprehensive analysis of phenotypic abnormalities in which the individual components of the phenotype are observed and described. The three components of the Human Phenotype Ontology (HPO; www.human-phenotype-ontology.org) project are the phenotype vocabulary, disease-phenotype annotations and the algorithms that operate on these. These components are being used for computational deep phenotyping and precision medicine as well as integration of clinical data into translational research. The HPO is being increasingly adopted as a standard for phenotypic abnormalities by diverse groups such as international rare disease organizations, registries, clinical labs, biomedical resources, and clinical software tools and will thereby contribute toward nascent efforts at global data exchange for identifying disease etiologies. This update article reviews the progress of the HPO project since the debut Nucleic Acids Research database article in 2014, including specific areas of expansion such as common (complex) disease, new algorithms for phenotype driven genomic discovery and diagnostics, integration of cross-species mapping efforts with the Mammalian Phenotype Ontology, an improved quality control pipeline, and the addition of patient-friendly terminology.

Phenopolis: an open platform for harmonization and analysis of genetic and phenotypic data.

SUMMARY: Phenopolis is an open-source web server providing an intuitive interface to genetic and phenotypic databases. It integrates analysis tools such as variant filtering and gene prioritization based on phenotype. The Phenopolis platform will accelerate clinical diagnosis, gene discovery and encourage wider adoption of the Human Phenotype Ontology in the study of rare genetic diseases. AVAILABILITY AND IMPLEMENTATION: A demo of the website is available at https://phenopolis.github.io . If you wish to install a local copy, source code and installation instruction are available at https://github.com/phenopolis . The software is implemented using Python, MongoDB, HTML/Javascript and various bash shell scripts. CONTACT: n.pontikos@ucl.ac.uk. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Benign Yellow Dot Maculopathy: A New Macular Phenotype.

PURPOSE: To describe a novel macular phenotype that is associated with normal visual function. DESIGN: Retrospective, observational case series. PARTICIPANTS: Thirty-six affected individuals from 23 unrelated families. METHODS: This was a retrospective study of patients who had a characteristic macular phenotype. Subjects underwent a full ocular examination, electrophysiologic studies, spectral-domain optical coherence tomography (OCT), and fundus autofluorescence imaging. Genomic analyses were performed using haplotype sharing analysis and whole-exome sequencing. MAIN OUTCOME MEASURES: Visual acuity, retinal features, electroretinography, and whole-exome sequencing. RESULTS: Twenty-six of 36 subjects were female. The median age of subjects at presentation was 15 years (range, 5-59 years). The majority of subjects were asymptomatic and presented after a routine eye examination (22/36 subjects) or after screening because of a positive family history (13/36 subjects) or by another ophthalmologist (1/36 subjects). Of the 3 symptomatic subjects, 2 had reduced visual acuity secondary to nonorganic visual loss and bilateral ametropic amblyopia with strabismus. Visual acuity was 0.18 logarithm of the minimum angle of resolution (logMAR) or better in 30 of 33 subjects. Color vision was normal in all subjects tested, except for the subject with nonorganic visual loss. All subjects had bilateral symmetric multiple yellow dots at the macula. In the majority of subjects, these were evenly distributed throughout the fovea, but in 9 subjects they were concentrated in the nasal parafoveal area. The dots were hyperautofluorescent on fundus autofluorescence imaging. The OCT imaging was generally normal, but in 6 subjects subtle irregularities at the inner segment ellipsoid band were seen. Electrophysiologic studies identified normal macular function in 17 of 19 subjects and normal full-field retinal function in all subjects. Whole-exome analysis across 3 unrelated families found no pathogenic variants in known macular dystrophy genes. Haplotype sharing analysis in 1 family excluded linkage with the North Carolina macular dystrophy (MCDR1) locus. CONCLUSIONS: A new retinal phenotype is described, which is characterized by bilateral multiple early-onset yellow dots at the macula. Visual function is normal, and the condition is nonprogressive. In familial cases, the phenotype seems to be inherited in an autosomal dominant manner, but a causative gene is yet to be ascertained.

Analysis of copy number variation at DMBT1 and age-related macular degeneration.

BACKGROUND: DMBT1 is a gene that shows extensive copy number variation (CNV) that alters the number of bacteria-binding domains in the protein and has been shown to activate the complement pathway. It lies next to the ARMS2/HTRA1 genes in a region of chromosome 10q26, where single nucleotide variants have been strongly associated with age-related macular degeneration (AMD), the commonest cause of blindness in Western populations. Complement activation is thought to be a key factor in the pathogenesis of this condition. We sought to investigate whether DMBT1 CNV plays any role in the susceptibility to AMD. METHODS: We analysed long-range linkage disequilibrium of DMBT1 CNV1 and CNV2 with flanking single nucleotide polymorphisms (SNPs) using our previously published CNV and HapMap Phase 3 SNP data in the CEPH Europeans from Utah (CEU). We then typed a large cohort of 860 AMD patients and 419 examined age-matched controls for copy number at DMBT1 CNV1 and CNV2 and combined these data with copy numbers from a further 480 unexamined controls. RESULTS: We found weak linkage disequilibrium between DMBT1 CNV1 and CNV2 with the SNPs rs1474526 and rs714816 in the HTRA1/ARMS2 region. By directly analysing copy number variation, we found no evidence of association of CNV1 or CNV2 with AMD. CONCLUSIONS: We have shown that copy number variation at DMBT1 does not affect risk of developing age-related macular degeneration and can therefore be ruled out from future studies investigating the association of structural variation at 10q26 with AMD.

Genome-wide association study of age-related macular degeneration identifies associated variants in the TNXB-FKBPL-NOTCH4 region of chromosome 6p21.3.

Age-related macular degeneration (AMD) is a leading cause of visual loss in Western populations. Susceptibility is influenced by age, environmental and genetic factors. Known genetic risk loci do not account for all the heritability. We therefore carried out a genome-wide association study of AMD in the UK population with 893 cases of advanced AMD and 2199 controls. This showed an association with the well-established AMD risk loci ARMS2 (age-related maculopathy susceptibility 2)-HTRA1 (HtrA serine peptidase 1) (P =2.7 × 10(-72)), CFH (complement factor H) (P =2.3 × 10(-47)), C2 (complement component 2)-CFB (complement factor B) (P =5.2 × 10(-9)), C3 (complement component 3) (P =2.2 × 10(-3)) and CFI (P =3.6 × 10(-3)) and with more recently reported risk loci at VEGFA (P =1.2 × 10(-3)) and LIPC (hepatic lipase) (P =0.04). Using a replication sample of 1411 advanced AMD cases and 1431 examined controls, we confirmed a novel association between AMD and single-nucleotide polymorphisms on chromosome 6p21.3 at TNXB (tenascin XB)-FKBPL (FK506 binding protein like) [rs12153855/rs9391734; discovery P =4.3 × 10(-7), replication P =3.0 × 10(-4), combined P =1.3 × 10(-9), odds ratio (OR) = 1.4, 95% confidence interval (CI) = 1.3-1.6] and the neighbouring gene NOTCH4 (Notch 4) (rs2071277; discovery P =3.2 × 10(-8), replication P =3.8 × 10(-5), combined P =2.0 × 10(-11), OR = 1.3, 95% CI = 1.2-1.4). These associations remained significant in conditional analyses which included the adjacent C2-CFB locus. TNXB, FKBPL and NOTCH4 are all plausible AMD susceptibility genes, but further research will be needed to identify the causal variants and determine whether any of these genes are involved in the pathogenesis of AMD.

Natural history and retinal structure in patients with Usher syndrome type 1 owing to MYO7A mutation.

PURPOSE: To evaluate the phenotypic variability and natural history of ocular disease in a cohort of 28 individuals with MYO7A-related disease. Mutations in the MYO7A gene are the most common cause of Usher syndrome type 1, characterized by profound congenital deafness, vestibular arreflexia, and progressive retinal degeneration. DESIGN: Retrospective case series. PARTICIPANTS: Twenty-eight patients from 26 families (age range, 3-65 years; median, 32) with 2 likely disease-causing variants in MYO7A. METHODS: Clinical investigations included fundus photography, optical coherence tomography, fundus autofluorescence (FAF) imaging, and audiologic and vestibular assessments. Longitudinal visual acuity and FAF data (over a 3-year period) were available for 20 and 10 study subjects, respectively. MAIN OUTCOME MEASURES: Clinical, structural, and functional characteristics. RESULTS: All patients with MYO7A mutations presented with features consistent with Usher type 1. The median visual acuity for the cohort was 0.39 logarithm of the minimum angle of resolution (logMAR; range, 0.0-2.7) and visual acuity in logMAR correlated with age (Spearman's rank correlation coefficient, r = 0.71; P<0.0001). Survival analysis revealed that acuity ≤ 0.22 logMAR was maintained in 50% of studied subjects until age 33.9; legal blindness based on loss of acuity (≥ 1.00 logMAR) or loss of field (≤ 20°) was reached at a median age of 40.6 years. Three distinct patterns were observed on FAF imaging: 13 of 22 patients tested had relatively preserved foveal autofluorescence surrounded by a ring of high density, 4 of 22 had increased signal in the fovea with no obvious hyperautofluorescent ring, and 5 of 22 had widespread hypoautofluorescence corresponding to retinal pigment epithelial atrophy. Despite a number of cases presenting with a milder phenotype, there seemed to be no obvious genotype-phenotype correlation. CONCLUSIONS: MYO7A-related ocular disease is variable. Central vision typically remains preserved at least until the third decade of life, with 50% of affected individuals reaching legal blindness by 40 years of age. Distinct phenotypic subsets were identified on FAF imaging. A specific allele, previously reported in nonsyndromic deafness, may be associated with a mild retinopathy.