Phenotype Driven Analysis of Whole Genome Sequencing Identifies Deep Intronic Variants that Cause Retinal Dystrophies by Aberrant Exonization.

Purpose: To demonstrate the effectiveness of combining retinal phenotyping and focused variant filtering from genome sequencing (GS) in identifying deep intronic disease causing variants in inherited retinal dystrophies. Methods: Affected members from three pedigrees with classical enhanced S-cone syndrome (ESCS; Pedigree 1), congenital stationary night blindness (CSNB; Pedigree 2), and achromatopsia (ACHM; Pedigree 3), respectively, underwent detailed ophthalmologic evaluation, optical coherence tomography, and electroretinography. The probands underwent panel-based genetic testing followed by GS analysis. Minigene constructs (NR2E3, GPR179 and CNGB3) and patient-derived cDNA experiments (NR2E3 and GPR179) were performed to assess the functional effect of the deep intronic variants. Results: The electrophysiological findings confirmed the clinical diagnosis of ESCS, CSNB, and ACHM in the respective pedigrees. Panel-based testing revealed heterozygous pathogenic variants in NR2E3 (NM_014249.3; c.119-2A>C; Pedigree 1) and CNGB3 (NM_019098.4; c.1148delC/p.Thr383Ilefs*13; Pedigree 3). The GS revealed heterozygous deep intronic variants in Pedigrees 1 (NR2E3; c.1100+1124G>A) and 3 (CNGB3; c.852+4751A>T), and a homozygous GPR179 variant in Pedigree 2 (NM_001004334.3; c.903+343G>A). The identified variants segregated with the phenotype in all pedigrees. All deep intronic variants were predicted to generate a splice acceptor gain causing aberrant exonization in NR2E3 [89 base pairs (bp)], GPR179 (197 bp), and CNGB3 (73 bp); splicing defects were validated through patient-derived cDNA experiments and/or minigene constructs and rescued by antisense oligonucleotide treatment. Conclusions: Deep intronic mutations contribute to missing heritability in retinal dystrophies. Combining results from phenotype-directed gene panel testing, GS, and in silico splice prediction tools can help identify these difficult-to-detect pathogenic deep intronic variants.

Retrotransposon insertion as a novel mutational event in Bardet-Biedl syndrome.

BACKGROUND: Bardet-Biedl syndrome (BBS) is an autosomal recessive pleiotropic disorder of the primary cilia that leads to severe visual loss in the teenage years. Approximately 80% of BBS cases are explained by mutations in one of the 21 identified genes. Documented causative mutation types include missense, nonsense, copy number variation (CNV), frameshift deletions or insertions, and splicing variants. METHODS: Whole genome sequencing was performed on a patient affected with BBS for whom no mutations were identified using clinically approved genetic testing of the known genes. Analysis of the WGS was done using internal protocols and publicly available algorithms. The phenotype was defined by retrospective chart review. RESULTS: We document a female affected with BBS carrying the most common BBS1 mutation (BBS1: Met390Arg) on the maternal allele and an insertion of a ~1.7-kb retrotransposon in exon 13 on the paternal allele. This retrotransposon insertion was not automatically annotated by the standard variant calling protocols used. This novel variant was identified by visual inspection of the alignment file followed by specific genome analysis with an available algorithm for transposable elements. CONCLUSION: This report documents a novel mutation type associated with BBS and highlights the importance of systematically performing transposon detection analysis on WGS data of unsolved cases.

OTX2 mutations cause autosomal dominant pattern dystrophy of the retinal pigment epithelium.

PURPOSE: To identify the genetic cause of autosomal-dominant pattern dystrophy (PD) of the retinal pigment epithelium (RPE) in two families. METHODS AND RESULTS: Two families with autosomal-dominant PD were identified. Eight members of family 1 (five affected) were subjected to whole-genome SNP genotyping; multipoint genome-wide linkage analysis identified 7 regions of potential linkage, and genotyping four additional individuals from family 1 resulted in a maximum logarithm of odds score of 2.09 observed across four chromosomal regions. Exome sequencing of two affected family 1 members identified 15 shared non-synonymous rare coding sequence variants within the linked regions; candidate genes were prioritised and further analysed. Sanger sequencing confirmed a novel heterozygous missense variant (E79K) in orthodenticle homeobox 2 (OTX2) that segregated with the disease phenotype. Family 2 with PD (two affected) harboured the same missense variant in OTX2. A shared haplotype of 19.68 cM encompassing OTX2 was identified between affected individuals in the two families. Within the two families, all except one affected demonstrated distinct 'patterns' at the macula. In vivo structural retinal imaging showed discrete areas of RPE-photoreceptor separation at the macula in all cases. Electroretinogram testing showed generalised photoreceptor degeneration in three cases. Mild developmental anomalies were observed, including optic nerve head dysplasia (four cases), microcornea (one case) and Rathke's cleft cyst (one case); pituitary hormone levels were normal. CONCLUSIONS: This is the first report implicating OTX2 to underlie PD. The retinal disease resembles conditional mice models that show slow photoreceptor degeneration secondary to loss of Otx2 function in the adult RPE.

Mutational analysis of SDCCAG8 in Bardet-Biedl syndrome patients with renal involvement and absent polydactyly.

PURPOSE: To assess for SDCCAG8 mutations in Bardet-Biedl syndrome (BBS) subjects with renal involvement and no polydactyly, and to describe phenotypic characteristics of SDCCAG8-related disease. MATERIAL AND METHODS: Five patients (from 4 pedigrees) with clinical diagnosis of BBS, who had retinal and renal involvement and no polydactyly, were assessed. Sequence analysis of SDCCAG8 was undertaken and a detailed clinical review of an affected sibship was performed. RESULTS: A sibship of East Indian origin who carried a putative clinical diagnosis of BBS had compound heterozygous mutations in SDCCAG8 (p.Thr482LysfsX12/p.Asp543AlafsX24). The renal involvement was early and required transplant in both cases. Both were short statured and had asthma since childhood. The younger sister also had non-alcoholic fatty liver disease. Visual acuity and central fields were preserved in the teenage years in both patients. The optical coherence tomography showed preservation of the retinal lamination at the fovea; fundus autofluorescence demonstrated a perifoveal ring of hyperfluorescence as commonly observed in other forms of retinitis pigmentosa. Full-field electroretinogram revealed rod function to be more severely affected than cone function in both cases. CONCLUSION: Our results and prior literature suggest that SDCCAG8 could play an important role in presumed BBS patients affected with severe kidney disease and absent polydactyly. This report enhances the phenotypic description of SDCCAG8-related disease.

BBS mutational analysis: a strategic approach.

BACKGROUND: Bardet-Biedl syndrome (BBS, OMIM 209900) is a rare autosomal recessive, clinically and genetically heterogeneous disorder with 15 genes identified. The large amount of coding sequence challenges the cost effectiveness of mutational analysis of BBS. MATERIAL AND METHODS: We present our mutational analysis experience (83 BBS families) in the context of the literature published up to September 2010, to provide a comprehensive tabulation of all BBS1-BBS12 mutant alleles and optimize a screening approach. RESULTS: We identified two BBS disease alleles in 76% of probands. Together BBS1, BBS2, BBS10 and BBS12 account for 82.4% of published unrelated alleles. On average 82% of published alleles are private. The 267 published principal mutations were positioned and analysis of their distribution allowed the design of a mutation screening strategy. Starting by screening for recurrent mutations, for example BBS1 M390R (10% of our cohort) and BBS10 C91LfsX5 (6% of our cohort), allowed a capture of 23.5% of the principal mutated alleles. Following a hierarchy of frequently involved exons, subsequent sequencing of the 4 most commonly involved genes, BBS1, BBS10, BBS2 and BBS12 could bring this mutation detection to at least 62%. The 16 most frequently recurring alleles could be identified with the use of simple screening methods such as restriction enzyme digest and ARMS assay and require sequencing in only a few instances. CONCLUSION: Our results suggest that mutational analysis of such a "rare" genetically heterogeneous condition benefits from pooling of data. This allows the development of efficient and cost-conscious screening mutational analysis strategies.

BBS genotype-phenotype assessment of a multiethnic patient cohort calls for a revision of the disease definition.

Bardet-Biedl syndrome (BBS) is a ciliopathy characterized by retinal degeneration, obesity, polydactyly, renal abnormalities, and cognitive impairment for which 15 causative genes have been identified. Here we present the results of a mutational analysis of our multiethnic cohort of 83 families (105 cases); 75.9% of them have their mutations identified including 26 novel changes. Comprehensive phenotyping of these patients demonstrate that the spectrum of clinical features is greater than expected and overlapped with the features of other ciliopathies; specifically Alström and McKusick-Kauffman syndromes.

Oligocone trichromacy is part of the spectrum of CNGA3-related cone system disorders.

PURPOSE: To report the rare observation of CNGA3 mutation as a cause of oligocone trichromacy (OT) and present phenotypic characteristics. METHODS: A 20 year old male patient underwent ophthalmological evaluation including detailed color vision assessment using Ishihara pseudoisochromatic plates, American Optical Hardy Rand Rittler plates (HRR) and Mollon-Reffin Minimalist test (MRM). Optical coherence tomography (OCT), fundus autofluorescence (FAF), visual field assessment and electrophysiological testing was also performed. The patient's DNA was sequenced for mutations in the coding sequence of CNGA3 and CNGB3 genes. RESULTS: Best corrected visual acuity (BCVA) was 20/50 and 20/30 in the right and left eyes respectively. His color vision was normal to Ishihara, HRR and MRM tests. Fundus appearance, FAF, OCT and Goldmann visual fields (GVF) were all normal. Humphrey visual field analysis (HVF) demonstrated reduced sensitivity and paracentral scotomas (5-20°). The full-field electroretinogram (ERG) showed normal rod responses and severely reduced cone responses. The multifocal electroretinogram (mfERG) was non-recordable above noise. Compound heterozygous mutations in exon 8 of the CNGA3 coding sequence were identified; c.1070 A > G (Tyr357Cys; novel) and c.1694 C > T (Thr565Met). Allele-specific polymerase chain reaction confirmed that the mutations were located on separate alleles. No mutations were identified in CNGB3. CONCLUSION: This is the second reported case of CNGA3 associated OT. Mutations in CNGA3 have previously been associated with incomplete and complete achromatopsia. This report confirms that OT forms the mildest end of the spectrum of CNGA3 related diseases.

Mutations in chaperonin-like BBS genes are a major contributor to disease development in a multiethnic Bardet-Biedl syndrome patient population.

BACKGROUND: Bardet-Biedl syndrome is a pleiotropic disorder with 14 BBS genes identified. BBS1, BBS2, BBS4, BBS5, BBS7, BBS8, and BBS9 form a complex called the BBSome, which is believed to recruit Rab8(GTP) to the primary cilium and promote ciliogenesis. The second group, the chaperonin-like proteins BBS6, BBS10, and BBS12, have been defined as a vertebrate-specific branch of the type II chaperonin superfamily. These may play a role in the regulation of BBSome assembly. METHODS AND RESULTS: Using sequence analysis, the role of BBS6, 10 and 12 was assessed in the patient population comprising 93 cases from 74 families. Systemic and ocular phenotypes were defined. In the study, chaperonin-like BBS gene mutations accounted for the disease in approximately 36.5% of BBS families. A total of 38 different non-polymorphic exonic sequence variants were identified in 40.5% of BBS families (41.9% cases), of which 26 were novel (68%). Six cases had mutations present in more than one chaperonin-like BBS gene. One case with four mutations in BBS10 had a phenotype of overall greater severity. The phenotypes observed were beyond the classic BBS phenotype as they overlapped with characteristics of MKKS (congenital heart defect, vaginal atresia, hydrometrocolpos, cryptorchidism), as well as Alström syndrome (diabetes, hearing loss, liver abnormalities, endocrine anomalies, cardiomyopathy). CONCLUSIONS: While overlap between the MKKS and BBS phenotypes has previously been reported for cases with BBS6 mutations, we also observed MKKS phenotypes involving BBS10 and BBS12 and Alström-like phenotypes associated with mutations in BBS1, BBS2, BBS6, BBS7, BBS9, BBS10 and BBS12 for the first time.

Clinical phenotype of posterior polymorphous corneal dystrophy in a family with a novel ZEB1 mutation.

PURPOSE: To describe the clinical phenotype in a family with posterior polymorphous corneal dystrophy (PPCD) and a novel mutation in the ZEB1 gene. METHODS: Clinical examination, anterior segment photography, specular microscopy and electrophysiological investigations were performed and quantified. Genomic DNA extracted from peripheral blood was sequenced for ZEB1 exons. Cosegregation of identified mutation with the disease status in the family was confirmed using polymerase chain reaction and restriction fragment length polymorphism. RESULTS: Ocular examination was performed on five family members from two generations. Three had anomalies of the corneal endothelium that were consistent with PPCD. Endothelial cell counts ranged from 2306 to 2987 mm(2) (ref. 2000-4000 cells/mm(2) ). No evidence of glaucoma or retinal abnormalities was observed. Extraocular abnormalities such as inguinal herniation, hydrocoele and possible bony or connective tissue anomalies were part of the disease spectrum in this family. Mutation analysis revealed a novel change in exon 5 of ZEB1 (c.672delA) that cosegregated with the affected disease status. CONCLUSION: The detailed clinical features of PPCD associated with a novel ZEB1 mutation are supportive of the previously proposed range of phenotype parameters. Further phenotype-genotype correlations may provide insights into the clinical variability and pathological processes affecting the corneal endothelium, Descemet's membrane, retinal photoreceptor function and extraocular tissues of some patients.

BBS7 and TTC8 (BBS8) mutations play a minor role in the mutational load of Bardet-Biedl syndrome in a multiethnic population.

Bardet Biedl syndrome is a genetically heterogeneous ciliopathy with fourteen genes currently identified. To date, mutations in BBS7 and TTC8 (BBS8) were reported in 4.2% and 2.8% of BBS families respectively. We sequenced the coding regions of BBS7 and TTC8 in 35 BBS families of diverse ancestral backgrounds. In addition, the role of putative modifier genes on phenotype severity; NXNL1 and MGC1203 c.430C>T, was assessed. Genotype-phenotype correlation was explored in patients with identified mutations. Four novel pathogenic BBS7 changes were identified in 2/35 families (5.7%). In one family with two affected individuals with BBS7 mutations, a more severe phenotype was observed in association with a third mutation in BBS4. The overall retinal phenotype appeared more severe than that seen in patients with BBS1 mutations. This study confirms the small role of BBS7 and TTC8 in the overall mutational load of BBS patients. The variability of the ocular phenotype observed, could not be explained by the putative modifier genes; NXNL1 and MGC1203 c.430C>T.

CRYBA4, a novel human cataract gene, is also involved in microphthalmia.

Genetic analysis of a large Indian family with an autosomal dominant cataract phenotype allowed us to identify a novel cataract gene, CRYBA4. After a genomewide screen, linkage analysis identified a maximum LOD score of 3.20 (recombination fraction [theta] 0.001) with marker D22S1167 of the beta -crystallin gene cluster on chromosome 22. To date, CRYBA4 was the only gene in this cluster not associated with either human or murine cataracts. A pathogenic mutation was identified in exon 4 that segregated with the disease status. The c.317T-->C sequence change is predicted to replace the highly conserved hydrophobic amino acid phenylalanine94 with the hydrophilic amino acid serine. Modeling suggests that this substitution would significantly reduce the intrinsic stability of the crystalline monomer, which would impair its ability to form the association modes critical for lens transparency. Considering that CRYBA4 associates with CRYBB2 and that the latter protein has been implicated in microphthalmia, mutational analysis of CRYBA4 was performed in 32 patients affected with microphthalmia (small eye). We identified a c.242T-->C (Leu69Pro) sequence change in exon 4 in one patient, which is predicted here to disrupt the beta -sheet structure in CRYBA4. Protein folding would consequently be impaired, most probably leading to a structure with reduced stability in the mutant. This is the first report linking mutations in CRYBA4 to cataractogenesis and microphthalmia.

Further support of the role of CYP1B1 in patients with Peters anomaly.

PURPOSE: Peters anomaly is a developmental anomaly of the eye frequently associated with glaucoma. The aim of this study was to further define the molecular basis of this condition. METHODS: The role of four candidate genes implicated in ocular development or glaucoma, PAX6, PITX2, MYOC, and CYP1B1, was studied in 15 patients with Peters anomaly. Mutational analysis used a combination of single strand conformation polymorphism (SSCP) and direct cycle sequencing. RESULTS: Four mutations in CYP1B1 were found in 3/15 (20%) affected individuals compared with 1/140 (0.7%) control individuals. CONCLUSIONS: This study supports the role of CYP1B1 as a causative gene in Peters anomaly. Furthermore, this emphasizes the broad range of phenotypic expression for CYP1B1 mutations, and its role in eye development.

CRYBB1 mutation associated with congenital cataract and microcornea.

PURPOSE: The molecular characterization of a UK family with an autosomal dominant congenital cataract associated with microcornea is reported. METHODS: Family history and clinical data were recorded. This phenotype was linked to a 7.6 cM region of chromosome 22q11.2-q12.2, spanning the beta-crystallin gene cluster (ZMax of 3.91 for marker D22S1114 at theta=0). Candidate genes were PCR amplified and screened for mutations on both strands using direct sequencing. RESULTS: Sequencing of the coding regions and flanking intronic sequences of CRYBB2 and CRYBB1 showed the presence of a novel, heterozygous X253R change in exon 6 of CRYBB1. SSCP analysis confirmed that this sequence change segregated with the disease phenotype in all available family members and was not found in 109 ethnically matched controls. CONCLUSIONS: X253R is predicted to elongate the COOH-terminal extension of the protein and would be expected to disrupt beta-crystallin interactions. This is the first documented involvement of CRYBB1 in ocular development beyond cataractogenesis.

Defining the pathogenicity of optineurin in juvenile open-angle glaucoma.

PURPOSE: Juvenile open-angle glaucoma (JOAG) differs from primary open-angle glaucoma in that it is usually a more severe phenotype and has an earlier age of onset. Optineurin was recently associated with a variant of POAG that is characterized by intraocular pressure within normal limits: normal-tension glaucoma. The present study tested whether OPTN sequence changes play a role in early-onset glaucoma characterized by elevated intraocular pressure. METHODS: Sixty-six patients with JOAG characterized by high intraocular pressure were screened for mutations. Mutational analysis was performed with a combination of restriction enzyme digestion, single-strand conformation polymorphism, and direct sequencing. The effects of select changes on exon splicing were assessed using bioinformatic modeling approaches and RT-PCR. RESULTS: Ten sequence changes were identified, of which H486R was strongly suggestive of pathogenicity. H486R represents the first reported OPTN mutation associated with JOAG. Also, L41L is proposed to confer an increased susceptibility to the development of JOAG. Most of the other sequence changes observed were not thought to be biologically significant. The frequency of the previously reported M98K allele was not increased in the JOAG population studied but showed the previously reported skewed distribution in the POAG study population. The changes identified were not shown to affect the splicing machinery. CONCLUSIONS: The results of this work support the hypothesis that mutations in OPTN are not specifically associated with low-pressure glaucoma, but can play a role in JOAG.

Further genetic heterogeneity for autosomal dominant human sutural cataracts.

A unique sutural cataract was observed in a 4-generation German family to be transmitted as an isolated autosomal, dominant trait. Since mutations in the gamma-crystallin encoding CRYG genes have previously been demonstrated to be the most frequent reason for isolated congenital cataracts, all 4 active CRYG genes have been sequenced. A single base-pair change in the CRYGA gene has been shown, leading to a premature stop codon. This was not observed in 170 control individuals. However, it did not segregate with the disease phenotype. This is the first truncating mutation in an active CRYG gene without a dominant phenotype. As the CRYGA mutation did not explain the cataract, several other candidate loci (CCV, GJA8, CRYBB2, BFSP2, MIP, GJA8, CENTRAL POUCH-LIKE, CRYBA1) were investigated by microsatellite markers and linkage analysis, but they were excluded based on the combination of haplotype analysis and two-point linkage analysis. The phenotype in this family is due to a mutation in another sutural cataract gene yet to be identified.

VSX1: a gene for posterior polymorphous dystrophy and keratoconus.

We identified mutations in the VSX1 homeobox gene for two distinct inherited corneal dystrophies; posterior polymorphous dystrophy (PPD) and keratoconus. One of the mutation (R166W) responsible for keratoconus altered the homeodomain and impaired DNA binding. Two other sequence changes (L159M and G160D) were associated with keratoconus and PPD, respectively, and involved a region adjacent to the homeodomain. The G160D substitution, and a fourth defect affecting the highly conserved CVC domain (P247R), occurred in a child with very severe PPD who required a corneal transplant at 3 months of age. In this family, relatives with the G160D change alone had mild to moderate PPD, while P247R alone caused no corneal abnormalities. However, with either the G160D or P247R mutation, electroretinography detected abnormal function of the inner retina, where VSX1 is expressed. These data define the molecular basis of two important corneal dystrophies and reveal the importance of the CVC domain in the human retina.

Digenic inheritance of early-onset glaucoma: CYP1B1, a potential modifier gene.

"Early-onset glaucoma" refers to genetically heterogeneous conditions for which glaucoma manifests at age 5-40 years and for which only a small subset is molecularly characterized. We studied the role of MYOC, CYP1B1, and PITX2 in a population (n=60) affected with juvenile or early-onset glaucoma from the greater Toronto area. By a combination of single-strand conformation polymorphism and direct cycle sequencing, MYOC mutations were detected in 8 (13.3%) of the 60 individuals, CYP1B1 mutations were detected in 3 (5%) of the 60 individuals, and no PITX2 mutations were detected. The range of phenotypic expression associated with MYOC and CYP1B1 mutations was greater than expected. MYOC mutations included cases of juvenile glaucoma with or without pigmentary glaucoma and mixed-mechanism glaucoma. CYP1B1 mutations involved cases of juvenile open-angle glaucoma, as well as cases of congenital glaucoma. The study of a family with autosomal dominant glaucoma showed the segregation of both MYOC and CYP1B1 mutations with disease; however, in this family, the mean age at onset of carriers of the MYOC mutation alone was 51 years (range 48-64 years), whereas carriers of both the MYOC and CYP1B1 mutations had an average age at onset of 27 years (range 23-38 years) (P=.001). This work emphasizes the genetic heterogeneity of juvenile glaucoma and suggests, for the first time, that (1) congenital glaucoma and juvenile glaucoma are allelic variants and (2) the spectrum of expression of MYOC and CYP1B1 mutations is greater than expected. We also propose that CYP1B1 may act as a modifier of MYOC expression and that these two genes may interact through a common pathway.