Demonstration of the pathogenicity of a common non-exomic mutation in ABCA4 using iPSC-derived retinal organoids and retrospective clinical data.

Mutations in ABCA4 are the most common cause of Mendelian retinal disease. Clinical evaluation of this gene is challenging because of its extreme allelic diversity, the large fraction of non-exomic mutations, and the wide range of associated disease. We used patient-derived retinal organoids as well as DNA samples and clinical data from a large cohort of patients with ABCA4-associated retinal disease to investigate the pathogenicity of a variant in ABCA4 (IVS30 + 1321 A > G) that occurs heterozygously in 2% of Europeans. We found that this variant causes mis-splicing of the gene in photoreceptor cells such that the resulting protein contains 36 incorrect amino acids followed by a premature stop. We also investigated the phenotype of 10 patients with compound genotypes that included this mutation. Their median age of first vision loss was 39 years, which is in the mildest quintile of a large cohort of patients with ABCA4 disease. We conclude that the IVS30 + 1321 A > G variant can cause disease when paired with a sufficiently deleterious opposing allele in a sufficiently permissive genetic background.

A Retrospective Longitudinal Study of 460 Patients with ABCA4-Associated Retinal Disease.

PURPOSE: To investigate the distribution of genotypes and natural history of ABCA4-associated retinal disease in a large cohort of patients seen at a single institution. DESIGN: Retrospective, single-institution cohort review. PARTICIPANTS: Patients seen at the University of Iowa between November 1986 and August 2022 clinically suspected to have disease caused by sequence variations in ABCA4. METHODS: DNA samples from participants were subjected to a tiered testing strategy progressing from allele-specific screening to whole genome sequencing. Charts were reviewed, and clinical data were tabulated. The pathogenic severity of the most common alleles was estimated by studying groups of patients who shared 1 allele. Groups of patients with shared genotypes were reviewed for evidence of modifying factor effects. MAIN OUTCOME MEASURES: Age at first uncorrectable vision loss, best-corrected visual acuity, and the area of the I2e isopter of the Goldmann visual field. RESULTS: A total of 460 patients from 390 families demonstrated convincing clinical features of ABCA4-associated retinal disease. Complete genotypes were identified in 399 patients, and partial genotypes were identified in 61. The median age at first vision loss was 16 years (range, 4-76 years). Two hundred sixty-five families (68%) harbored a unique genotype, and no more than 10 patients shared any single genotype. Review of the patients with shared genotypes revealed evidence of modifying factors that in several cases resulted in a > 15-year difference in age at first vision loss. Two hundred forty-one different alleles were identified among the members of this cohort, and 161 of these (67%) were found in only a single individual. CONCLUSIONS: ABCA4-associated retinal disease ranges from a very severe photoreceptor disease with an onset before 5 years of age to a late-onset retinal pigment epithelium-based condition resembling pattern dystrophy. Modifying factors frequently impact the ABCA4 disease phenotype to a degree that is similar in magnitude to the detectable ABCA4 alleles themselves. It is likely that most patients in any cohort will harbor a unique genotype. The latter observations taken together suggest that patients' clinical findings in most cases will be more useful for predicting their clinical course than their genotype. FINANCIAL DISCLOSURE(S): Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Human photoreceptor cells from different macular subregions have distinct transcriptional profiles.

The human neural retina is a light sensitive tissue with remarkable spatial and cellular organization. Compared with the periphery, the central retina contains more densely packed cone photoreceptor cells with unique morphologies and synaptic wiring. Some regions of the central retina exhibit selective degeneration or preservation in response to retinal disease and the basis for this variation is unknown. In this study, we used both bulk and single-cell RNA sequencing to compare gene expression within concentric regions of the central retina. We identified unique gene expression patterns of foveal cone photoreceptor cells, including many foveal-enriched transcription factors. In addition, we found that the genes RORB1, PPFIA1 and KCNAB2 are differentially spliced in the foveal, parafoveal and macular regions. These results provide a highly detailed spatial characterization of the retinal transcriptome and highlight unique molecular features of different retinal regions.

Exome-based investigation of the genetic basis of human pigmentary glaucoma.

BACKGROUND: Glaucoma is a leading cause of visual disability and blindness. Release of iris pigment within the eye, pigment dispersion syndrome (PDS), can lead to one type of glaucoma known as pigmentary glaucoma. PDS has a genetic component, however, the genes involved with this condition are largely unknown. We sought to discover genes that cause PDS by testing cohorts of patients and controls for mutations using a tiered analysis of exome data. RESULTS: Our primary analysis evaluated melanosome-related genes that cause dispersion of iris pigment in mice (TYRP1, GPNMB, LYST, DCT, and MITF). We identified rare mutations, but they were not statistically enriched in PDS patients. Our secondary analyses examined PMEL (previously linked with PDS), MRAP, and 19 other genes. Four MRAP mutations were identified in PDS cases but not in controls (p = 0.016). Immunohistochemical analysis of human donor eyes revealed abundant MRAP protein in the iris, the source of pigment in PDS. However, analysis of MRAP in additional cohorts (415 cases and 1645 controls) did not support an association with PDS. We also did not confirm a link between PMEL and PDS in our cohorts due to lack of reported mutations and similar frequency of the variants in PDS patients as in control subjects. CONCLUSIONS: We did not detect a statistical enrichment of mutations in melanosome-related genes in human PDS patients and we found conflicting data about the likely pathogenicity of MRAP mutations. PDS may have a complex genetic basis that is not easily unraveled with exome analyses.

Spectacle: An interactive resource for ocular single-cell RNA sequencing data analysis.

Single-cell RNA sequencing has revolutionized ocular gene expression studies. This technology has enabled researchers to identify expression signatures for rare cell types and characterize how gene expression changes across biological conditions, such as topographic region or disease status. However, sharing single-cell RNA sequencing results remains a major obstacle, particular for individuals without a computational background. To address these limitations, we developed Spectacle, an interactive web-based resource for exploring previously published single-cell RNA sequencing data from ocular studies. Spectacle is powered by a locally developed R package, cellcuratoR, which utilizes the Shiny framework in R to generate interactive visualizations for single-cell expression data. Spectacle contains five pre-processed ocular single-cell RNA sequencing data sets and is accessible via the web at OcularGeneExpression.org/singlecell. With Spectacle, users can interactively identify which cell types express a gene of interest, detect transcriptomic subpopulations within a cell type, and perform highly flexible differential expression analyses. The freely-available Spectacle system reduces the bioinformatic barrier for interacting with rich single-cell RNA sequencing studies from ocular tissues, making it easy to quickly identify cell types that express a gene of interest.

Primary congenital and developmental glaucomas.

Glaucoma is the leading cause of irreversible blindness worldwide. Although most glaucoma patients are elderly, congenital glaucoma and glaucomas of childhood are also important causes of visual disability. Primary congenital glaucoma (PCG) is isolated, non-syndromic glaucoma that occurs in the first three years of life and is a major cause of childhood blindness. Other early-onset glaucomas may arise secondary to developmental abnormalities, such as glaucomas that occur with aniridia or as part of Axenfeld-Rieger syndrome. Congenital and childhood glaucomas have strong genetic bases and disease-causing mutations have been discovered in several genes. Mutations in three genes (CYP1B1, LTBP2, TEK) have been reported in PCG patients. Axenfeld-Rieger syndrome is caused by mutations in PITX2 or FOXC1 and aniridia is caused by PAX6 mutations. This review discusses the roles of these genes in primary congenital glaucoma and glaucomas of childhood.

Hypomorphic mutations in TRNT1 cause retinitis pigmentosa with erythrocytic microcytosis.

Retinitis pigmentosa (RP) is a highly heterogeneous group of disorders characterized by degeneration of the retinal photoreceptor cells and progressive loss of vision. While hundreds of mutations in more than 100 genes have been reported to cause RP, discovering the causative mutations in many patients remains a significant challenge. Exome sequencing in an individual affected with non-syndromic RP revealed two plausibly disease-causing variants in TRNT1, a gene encoding a nucleotidyltransferase critical for tRNA processing. A total of 727 additional unrelated individuals with molecularly uncharacterized RP were completely screened for TRNT1 coding sequence variants, and a second family was identified with two members who exhibited a phenotype that was remarkably similar to the index patient. Inactivating mutations in TRNT1 have been previously shown to cause a severe congenital syndrome of sideroblastic anemia, B-cell immunodeficiency, recurrent fevers and developmental delay (SIFD). Complete blood counts of all three of our patients revealed red blood cell microcytosis and anisocytosis with only mild anemia. Characterization of TRNT1 in patient-derived cell lines revealed reduced but detectable TRNT1 protein, consistent with partial function. Suppression of trnt1 expression in zebrafish recapitulated several features of the human SIFD syndrome, including anemia and sensory organ defects. When levels of trnt1 were titrated, visual dysfunction was found in the absence of other phenotypes. The visual defects in the trnt1-knockdown zebrafish were ameliorated by the addition of exogenous human TRNT1 RNA. Our findings indicate that hypomorphic TRNT1 mutations can cause a recessive disease that is almost entirely limited to the retina.

Clinically Focused Molecular Investigation of 1000 Consecutive Families with Inherited Retinal Disease.

PURPOSE: To devise a comprehensive multiplatform genetic testing strategy for inherited retinal disease and to describe its performance in 1000 consecutive families seen by a single clinician. DESIGN: Retrospective series. PARTICIPANTS: One thousand consecutive families seen by a single clinician. METHODS: The clinical records of all patients seen by a single retina specialist between January 2010 and June 2016 were reviewed, and all patients who met the clinical criteria for a diagnosis of inherited retinal disease were included in the study. Each patient was assigned to 1 of 62 diagnostic categories, and this clinical diagnosis was used to define the scope and order of the molecular investigations that were performed. The number of nucleotides evaluated in a given subject ranged from 2 to nearly 900 000. MAIN OUTCOME MEASURES: Sensitivity and false genotype rate. RESULTS: Disease-causing genotypes were identified in 760 families (76%). These genotypes were distributed across 104 different genes. More than 75% of these 104 genes have coding sequences small enough to be packaged efficiently into an adeno-associated virus. Mutations in ABCA4 were the most common cause of disease in this cohort (173 families), whereas mutations in 80 genes caused disease in 5 or fewer families (i.e., 0.5% or less). Disease-causing genotypes were identified in 576 of the families without next-generation sequencing (NGS). This included 23 families with mutations in the repetitive region of RPGR exon 15 that would have been missed by NGS. Whole-exome sequencing of the remaining 424 families revealed mutations in an additional 182 families, and whole-genome sequencing of 4 of the remaining 242 families revealed 2 additional genotypes that were invisible by the other methods. Performing the testing in a clinically focused tiered fashion would be 6.1% more sensitive and 17.7% less expensive and would have a significantly lower average false genotype rate than using whole-exome sequencing to assess more than 300 genes in all patients (7.1% vs. 128%; P < 0.001). CONCLUSIONS: Genetic testing for inherited retinal disease is now more than 75% sensitive. A clinically directed tiered testing strategy can increase sensitivity and improve statistical significance without increasing cost.

North Carolina Macular Dystrophy Is Caused by Dysregulation of the Retinal Transcription Factor PRDM13.

PURPOSE: To identify specific mutations causing North Carolina macular dystrophy (NCMD). DESIGN: Whole-genome sequencing coupled with reverse transcription polymerase chain reaction (RT-PCR) analysis of gene expression in human retinal cells. PARTICIPANTS: A total of 141 members of 12 families with NCMD and 261 unrelated control individuals. METHODS: Genome sequencing was performed on 8 affected individuals from 3 families affected with chromosome 6-linked NCMD (MCDR1) and 2 individuals affected with chromosome 5-linked NCMD (MCDR3). Variants observed in the MCDR1 locus with frequencies <1% in published databases were confirmed using Sanger sequencing. Confirmed variants absent from all published databases were sought in 8 additional MCDR1 families and 261 controls. The RT-PCR analysis of selected genes was performed in stem cell-derived human retinal cells. MAIN OUTCOME MEASURES: Co-segregation of rare genetic variants with disease phenotype. RESULTS: Five sequenced individuals with MCDR1-linked NCMD shared a haplotype of 14 rare variants spanning 1 Mb of the disease-causing allele. One of these variants (V1) was absent from all published databases and all 261 controls, but was found in 5 additional NCMD kindreds. This variant lies in a DNase 1 hypersensitivity site (DHS) upstream of both the PRDM13 and CCNC genes. Sanger sequencing of 1 kb centered on V1 was performed in the remaining 4 NCMD probands, and 2 additional novel single nucleotide variants (V2 in 3 families and V3 in 1 family) were identified in the DHS within 134 bp of the location of V1. A complete duplication of the PRDM13 gene was also discovered in a single family (V4). The RT-PCR analysis of PRDM13 expression in developing retinal cells revealed marked developmental regulation. Next-generation sequencing of 2 individuals with MCDR3-linked NCMD revealed a 900-kb duplication that included the entire IRX1 gene (V5). The 5 mutations V1 to V5 segregated perfectly in the 102 affected and 39 unaffected members of the 12 NCMD families. CONCLUSIONS: We identified 5 rare mutations, each capable of arresting human macular development. Four of these strongly implicate the involvement of PRDM13 in macular development, whereas the pathophysiologic mechanism of the fifth remains unknown but may involve the developmental dysregulation of IRX1.

A targeted sequencing study of glutamatergic candidate genes in suicide attempters with bipolar disorder.

Suicidal behavior has been shown to have a heritable component that is partly driven by psychiatric disorders [Brent and Mann, 2005]. However, there is also an independent factor contributing to the heritability of suicidal behavior. We previously conducted a genome-wide association study (GWAS) of bipolar suicide attempters and bipolar non-attempters to assess this independent factor [Willour et al., 2012]. This GWAS implicated glutamatergic neurotransmission in attempted suicide. In the current study, we have conducted a targeted next-generation sequencing study of the glutamatergic N-methyl-D-aspartate (NMDA) receptor, neurexin, and neuroligin gene families in 476 bipolar suicide attempters and 473 bipolar non-attempters. The goal of this study was to gather sequence information from coding and regulatory regions of these glutamatergic genes to identify variants associated with attempted suicide. We identified 186 coding variants and 4,298 regulatory variants predicted to be functional in these genes. No individual variants were overrepresented in cases or controls to a degree that was statistically significant after correction for multiple testing. Additionally, none of the gene-level results were statistically significant following correction. While this study provides no direct support for a role of the examined glutamatergic candidate genes, further sequencing in expanded gene sets and datasets will be required to ultimately determine whether genetic variation in glutamatergic signaling influences suicidal behavior. © 2016 Wiley Periodicals, Inc.

Non-exomic and synonymous variants in ABCA4 are an important cause of Stargardt disease.

Mutations in ABCA4 cause Stargardt disease and other blinding autosomal recessive retinal disorders. However, sequencing of the complete coding sequence in patients with clinical features of Stargardt disease sometimes fails to detect one or both mutations. For example, among 208 individuals with clear clinical evidence of ABCA4 disease ascertained at a single institution, 28 had only one disease-causing allele identified in the exons and splice junctions of the primary retinal transcript of the gene. Haplotype analysis of these 28 probands revealed 3 haplotypes shared among ten families, suggesting that 18 of the 28 missing alleles were rare enough to be present only once in the cohort. We hypothesized that mutations near rare alternate splice junctions in ABCA4 might cause disease by increasing the probability of mis-splicing at these sites. Next-generation sequencing of RNA extracted from human donor eyes revealed more than a dozen alternate exons that are occasionally incorporated into the ABCA4 transcript in normal human retina. We sequenced the genomic DNA containing 15 of these minor exons in the 28 one-allele subjects and observed five instances of two different variations in the splice signals of exon 36.1 that were not present in normal individuals (P < 10(-6)). Analysis of RNA obtained from the keratinocytes of patients with these mutations revealed the predicted alternate transcript. This study illustrates the utility of RNA sequence analysis of human donor tissue and patient-derived cell lines to identify mutations that would be undetectable by exome sequencing.

Cordova: web-based management of genetic variation data.

UNLABELLED: Cordova is an out-of-the-box solution for building and maintaining an online database of genetic variations integrated with pathogenicity prediction results from popular algorithms. Our primary motivation for developing this system is to aid researchers and clinician-scientists in determining the clinical significance of genetic variations. To achieve this goal, Cordova provides an interface to review and manually or computationally curate genetic variation data as well as share it for clinical diagnostics and the advancement of research. AVAILABILITY AND IMPLEMENTATION: Cordova is open source under the MIT license and is freely available for download at https://github.com/clcg/cordova.

Transcriptomic analysis across nasal, temporal, and macular regions of human neural retina and RPE/choroid by RNA-Seq.

Proper spatial differentiation of retinal cell types is necessary for normal human vision. Many retinal diseases, such as Best disease and male germ cell associated kinase (MAK)-associated retinitis pigmentosa, preferentially affect distinct topographic regions of the retina. While much is known about the distribution of cell types in the retina, the distribution of molecular components across the posterior pole of the eye has not been well-studied. To investigate regional difference in molecular composition of ocular tissues, we assessed differential gene expression across the temporal, macular, and nasal retina and retinal pigment epithelium (RPE)/choroid of human eyes using RNA-Seq. RNA from temporal, macular, and nasal retina and RPE/choroid from four human donor eyes was extracted, poly-A selected, fragmented, and sequenced as 100 bp read pairs. Digital read files were mapped to the human genome and analyzed for differential expression using the Tuxedo software suite. Retina and RPE/choroid samples were clearly distinguishable at the transcriptome level. Numerous transcription factors were differentially expressed between regions of the retina and RPE/choroid. Photoreceptor-specific genes were enriched in the peripheral samples, while ganglion cell and amacrine cell genes were enriched in the macula. Within the RPE/choroid, RPE-specific genes were upregulated at the periphery while endothelium associated genes were upregulated in the macula. Consistent with previous studies, BEST1 expression was lower in macular than extramacular regions. The MAK gene was expressed at lower levels in macula than in extramacular regions, but did not exhibit a significant difference between nasal and temporal retina. The regional molecular distinction is greatest between macula and periphery and decreases between different peripheral regions within a tissue. Datasets such as these can be used to prioritize candidate genes for possible involvement in retinal diseases with regional phenotypes.

Advancing genetic testing for deafness with genomic technology.

BACKGROUND: Non-syndromic hearing loss (NSHL) is the most common sensory impairment in humans. Until recently its extreme genetic heterogeneity precluded comprehensive genetic testing. Using a platform that couples targeted genomic enrichment (TGE) and massively parallel sequencing (MPS) to sequence all exons of all genes implicated in NSHL, we tested 100 persons with presumed genetic NSHL and in so doing established sequencing requirements for maximum sensitivity and defined MPS quality score metrics that obviate Sanger validation of variants. METHODS: We examined DNA from 100 sequentially collected probands with presumed genetic NSHL without exclusions due to inheritance, previous genetic testing, or type of hearing loss. We performed TGE using post-capture multiplexing in variable pool sizes followed by Illumina sequencing. We developed a local Galaxy installation on a high performance computing cluster for bioinformatics analysis. RESULTS: To obtain maximum variant sensitivity with this platform 3.2-6.3 million total mapped sequencing reads per sample were required. Quality score analysis showed that Sanger validation was not required for 95% of variants. Our overall diagnostic rate was 42%, but this varied by clinical features from 0% for persons with asymmetric hearing loss to 56% for persons with bilateral autosomal recessive NSHL. CONCLUSIONS: These findings will direct the use of TGE and MPS strategies for genetic diagnosis for NSHL. Our diagnostic rate highlights the need for further research on genetic deafness focused on novel gene identification and an improved understanding of the role of non-exonic mutations. The unsolved families we have identified provide a valuable resource to address these areas.

Exome sequencing and analysis of induced pluripotent stem cells identify the cilia-related gene male germ cell-associated kinase (MAK) as a cause of retinitis pigmentosa.

Retinitis pigmentosa (RP) is a genetically heterogeneous heritable disease characterized by apoptotic death of photoreceptor cells. We used exome sequencing to identify a homozygous Alu insertion in exon 9 of male germ cell-associated kinase (MAK) as the cause of disease in an isolated individual with RP. Screening of 1,798 unrelated RP patients identified 20 additional probands homozygous for this insertion (1.2%). All 21 affected probands are of Jewish ancestry. MAK encodes a kinase involved in the regulation of photoreceptor-connecting cilium length. Immunohistochemistry of human donor tissue revealed that MAK is expressed in the inner segments, cell bodies, and axons of rod and cone photoreceptors. Several isoforms of MAK that result from alternative splicing were identified. Induced pluripotent stem cells were derived from the skin of the proband and a patient with non-MAK-associated RP (RP control). In the RP control individual, we found that a transcript lacking exon 9 was predominant in undifferentiated cells, whereas a transcript bearing exon 9 and a previously unrecognized exon 12 predominated in cells that were differentiated into retinal precursors. However, in the proband with the Alu insertion, the developmental switch to the MAK transcript bearing exons 9 and 12 did not occur. In addition to showing the use of induced pluripotent stem cells to efficiently evaluate the pathogenicity of specific mutations in relatively inaccessible tissues like retina, this study reveals algorithmic and molecular obstacles to the discovery of pathogenic insertions and suggests specific changes in strategy that can be implemented to more fully harness the power of sequencing technologies.

Carcinoembryonic antigen-related cell adhesion molecule 16 interacts with alpha-tectorin and is mutated in autosomal dominant hearing loss (DFNA4).

We report on a secreted protein found in mammalian cochlear outer hair cells (OHC) that is a member of the carcinoembryonic antigen-related cell adhesion molecule (CEACAM) family of adhesion proteins. Ceacam16 mRNA is expressed in OHC, and its protein product localizes to the tips of the tallest stereocilia and the tectorial membrane (TM). This specific localization suggests a role in maintaining the integrity of the TM as well as in the connection between the OHC stereocilia and TM, a linkage essential for mechanical amplification. In agreement with this role, CEACAM16 colocalizes and coimmunoprecipitates with the TM protein α-tectorin. In addition, we show that mutation of CEACAM16 leads to autosomal dominant nonsyndromic deafness (ADNSHL) at the autosomal dominant hearing loss (DFNA4) locus. In aggregate, these data identify CEACAM16 as an α-tectorin-interacting protein that concentrates at the point of attachment of the TM to the stereocilia and, when mutated, results in ADNSHL at the DFNA4 locus.

Prioritization of retinal disease genes: an integrative approach.

The discovery of novel disease-associated variations in genes is often a daunting task in highly heterogeneous disease classes. We seek a generalizable algorithm that integrates multiple publicly available genomic data sources in a machine-learning model for the prioritization of candidates identified in patients with retinal disease. To approach this problem, we generate a set of feature vectors from publicly available microarray, RNA-seq, and ChIP-seq datasets of biological relevance to retinal disease, to observe patterns in gene expression specificity among tissues of the body and the eye, in addition to photoreceptor-specific signals by the CRX transcription factor. Using these features, we describe a novel algorithm, positive and unlabeled learning for prioritization (PULP). This article compares several popular supervised learning techniques as the regression function for PULP. The results demonstrate a highly significant enrichment for previously characterized disease genes using a logistic regression method. Finally, a comparison of PULP with the popular gene prioritization tool ENDEAVOUR shows superior prioritization of retinal disease genes from previous studies. The java source code, compiled binary, assembled feature vectors, and instructions are available online at https://github.com/ahwagner/PULP.

AudioGene: predicting hearing loss genotypes from phenotypes to guide genetic screening.

Autosomal dominant nonsyndromic hearing loss (ADNSHL) is a common and often progressive sensory deficit. ADNSHL displays a high degree of genetic heterogeneity and varying rates of progression. Accurate, comprehensive, and cost-effective genetic testing facilitates genetic counseling and provides valuable prognostic information to affected individuals. In this article, we describe the algorithm underlying AudioGene, a software system employing machine-learning techniques that utilizes phenotypic information derived from audiograms to predict the genetic cause of hearing loss in persons segregating ADNSHL. Our data show that AudioGene has an accuracy of 68% in predicting the causative gene within its top three predictions, as compared with 44% for a majority classifier. We also show that AudioGene remains effective for audiograms with high levels of clinical measurement noise. We identify audiometric outliers for each genetic locus and hypothesize that outliers may reflect modifying genetic effects. As personalized genomic medicine becomes more common, AudioGene will be increasingly useful as a phenotypic filter to assess pathogenicity of variants identified by massively parallel sequencing.

Exon-level expression profiling of ocular tissues.

The normal gene expression profiles of the tissues in the eye are a valuable resource for considering genes likely to be involved with disease processes. We profiled gene expression in ten ocular tissues from human donor eyes using Affymetrix Human Exon 1.0 ST arrays. Ten different tissues were obtained from six different individuals and RNA was pooled. The tissues included: retina, optic nerve head (ONH), optic nerve (ON), ciliary body (CB), trabecular meshwork (TM), sclera, lens, cornea, choroid/retinal pigment epithelium (RPE) and iris. Expression values were compared with publically available Expressed Sequence Tag (EST) and RNA-sequencing resources. Known tissue-specific genes were examined and they demonstrated correspondence of expression with the representative ocular tissues. The estimated gene and exon level abundances are available online at the Ocular Tissue Database.

Comprehensive genetic testing for hereditary hearing loss using massively parallel sequencing.

The extreme genetic heterogeneity of nonsyndromic hearing loss (NSHL) makes genetic diagnosis expensive and time consuming using available methods. To assess the feasibility of target-enrichment and massively parallel sequencing technologies to interrogate all exons of all genes implicated in NSHL, we tested nine patients diagnosed with hearing loss. Solid-phase (NimbleGen) or solution-based (SureSelect) sequence capture, followed by 454 or Illumina sequencing, respectively, were compared. Sequencing reads were mapped using GSMAPPER, BFAST, and BOWTIE, and pathogenic variants were identified using a custom-variant calling and annotation pipeline (ASAP) that incorporates publicly available in silico pathogenicity prediction tools (SIFT, BLOSUM, Polyphen2, and Align-GVGD). Samples included one negative control, three positive controls (one biological replicate), and six unknowns (10 samples total), in which we genotyped 605 single nucleotide polymorphisms (SNPs) by Sanger sequencing to measure sensitivity and specificity for SureSelect-Illumina and NimbleGen-454 methods at saturating sequence coverage. Causative mutations were identified in the positive controls but not in the negative control. In five of six idiopathic hearing loss patients we identified the pathogenic mutation. Massively parallel sequencing technologies provide sensitivity, specificity, and reproducibility at levels sufficient to perform genetic diagnosis of hearing loss.