Transcriptome analyses of the human retina identify unprecedented transcript diversity and 3.5 Mb of novel transcribed sequence via significant alternative splicing and novel genes

Background: The retina is a complex tissue comprised of multiple cell types that is affected by a diverse set of diseases that are important causes of vision loss. Characterizing the transcripts, both annotated and novel, that are expressed in a given tissue has become vital for understanding the mechanisms underlying the pathology of disease. Results: We sequenced RNA prepared from three normal human retinas and characterized the retinal transcriptome at an unprecedented level due to the increased depth of sampling provided by the RNA-seq approach. We used a non-redundant reference transcriptome from all of the empirically-determined human reference tracks to identify annotated and novel sequences expressed in the retina. We detected 79,915 novel alternative splicing events, including 29,887 novel exons, 21,757 3′ and 5′ alternate splice sites, and 28,271 exon skipping events. We also identified 116 potential novel genes. These data represent a significant addition to the annotated human transcriptome. For example, the novel exons detected increase the number of identified exons by 3%. Using a high-throughput RNA capture approach to validate 14,696 of these novel transcriptome features we found that 99% of the putative novel events can be reproducibly detected. Further, 15-36% of the novel splicing events maintain an open reading frame, suggesting they produce novel protein products. Conclusions: To our knowledge, this is the first application of RNA capture to perform large-scale validation of novel transcriptome features. In total, these analyses provide extensive detail about a previously uncharacterized level of transcript diversity in the human retina

Transcriptome analyses of the human retina identify unprecedented transcript diversity and 3.5 Mb of novel transcribed sequence via significant alternative splicing and novel genes

Background: The retina is a complex tissue comprised of multiple cell types that is affected by a diverse set of diseases that are important causes of vision loss. Characterizing the transcripts, both annotated and novel, that are expressed in a given tissue has become vital for understanding the mechanisms underlying the pathology of disease. Results: We sequenced RNA prepared from three normal human retinas and characterized the retinal transcriptome at an unprecedented level due to the increased depth of sampling provided by the RNA-seq approach. We used a non-redundant reference transcriptome from all of the empirically-determined human reference tracks to identify annotated and novel sequences expressed in the retina. We detected 79,915 novel alternative splicing events, including 29,887 novel exons, 21,757 3′ and 5′ alternate splice sites, and 28,271 exon skipping events. We also identified 116 potential novel genes. These data represent a significant addition to the annotated human transcriptome. For example, the novel exons detected increase the number of identified exons by 3%. Using a high-throughput RNA capture approach to validate 14,696 of these novel transcriptome features we found that 99% of the putative novel events can be reproducibly detected. Further, 15-36% of the novel splicing events maintain an open reading frame, suggesting they produce novel protein products. Conclusions: To our knowledge, this is the first application of RNA capture to perform large-scale validation of novel transcriptome features. In total, these analyses provide extensive detail about a previously uncharacterized level of transcript diversity in the human retina

NMNAT1 Mutations Cause Leber Congenital Amaurosis

Leber congenital amaurosis (LCA) is an infantile-onset form of inherited retinal degeneration characterized by severe vision loss. Two-thirds of LCA cases are caused by mutations in 17 known disease genes (RetNet Retinal Information Network). Using exome sequencing, we identified a homozygous missense mutation (c.25G>A, p.Val9Met) in NMNAT1 as likely disease-causing in two siblings of a consanguineous Pakistani kindred affected by LCA. This mutation segregated with disease in their kindred, including in three other children with LCA. NMNAT1 resides in the previously identified LCA9 locus and encodes the nuclear isoform of nicotinamide mononucleotide adenylyltransferase, a rate-limiting enzyme in nicotinamide adenine dinucleotide $(NAD^+)$ biosynthesis. Functional studies showed the p.Val9Met mutation decreased NMNAT1 enzyme activity. Sequencing NMNAT1 in 284 unrelated LCA families identified 14 rare mutations in 13 additional affected individuals. These results are the first to link an NMNAT isoform to disease and indicate that NMNAT1 mutations cause LCA

Identification of the Photoreceptor Transcriptional Co-Repressor SAMD11 as Novel Cause of Autosomal Recessive Retinitis Pigmentosa

Retinitis pigmentosa (RP), the most frequent form of inherited retinal dystrophy is characterized by progressive photoreceptor degeneration. Many genes have been implicated in RP development, but several others remain to be identified. Using a combination of homozygosity mapping, whole-exome and targeted next-generation sequencing, we found a novel homozygous nonsense mutation in SAMD11 in five individuals diagnosed with adult-onset RP from two unrelated consanguineous Spanish families. SAMD11 is ortholog to the mouse major retinal SAM domain (mr-s) protein that is implicated in CRX-mediated transcriptional regulation in the retina. Accordingly, protein-protein network analysis revealed a significant interaction of SAMD11 with CRX. Immunoblotting analysis confirmed strong expression of SAMD11 in human retina. Immunolocalization studies revealed SAMD11 was detected in the three nuclear layers of the human retina and interestingly differential expression between cone and rod photoreceptors was observed. Our study strongly implicates SAMD11 as novel cause of RP playing an important role in the pathogenesis of human degeneration of photoreceptors.This work was supported by several grants from the Spanish Centre for Biomedical Network Research on Rare Diseases (CIBERER)(06/07/0036), Instituto de Salud Carlos III (ISCIII, Spanish Ministry of Health)/FEDER, including FIS (PI013/00226) and RETICS (RD09/0076/00101 and RD12/0034/0010), Ministry of Economy and Competitiveness (MINECO), including FEDER (BFU2012-36845), and BIO2011-27069, Conselleria de Educació of the Valencia Community (PROMETEOII/2014/025), Spanish National Organization of the Blind (ONCE) and the Spanish Fighting Blindness Foundation (FUNDALUCE). M.C. was sponsored by the Miguel Servet Program for Researchers in the Spanish National Health Service (CP12/03256) and RSA by Sara Borrel Postdoctoral Program (CD12/00676), both from the ISCIII/FEDER. A.A-F. was sponsored by CIBERER, RPC is supported by Fundación Conchita Rábago (FCR), L.C is sponsored by RETICS (RD12/0034/0010) from ISCIII and L.d.S. was supported by CAPES Foundation, Ministry of Education of Brazil

Efficient In Silico Identification of a Common Insertion in the MAK Gene which Causes Retinitis Pigmentosa

Background: Next generation sequencing (NGS) offers a rapid and comprehensive method of screening for mutations associated with retinitis pigmentosa and related disorders. However, certain sequence alterations such as large insertions or deletions may remain undetected using standard NGS pipelines. One such mutation is a recently-identified Alu insertion into the Male Germ Cell-Associated Kinase (MAK) gene, which is missed by standard NGS-based variant callers. Here, we developed an in silico method of searching NGS raw sequence reads to detect this mutation, without the need to recalculate sequence alignments or to screen every sample by PCR. Methods: The Linux program grep was used to search for a 23 bp “probe” sequence containing the known junction sequence of the insert. A corresponding search was performed with the wildtype sequence. The matching reads were counted and further compared to the known sequences of the full wildtype and mutant genomic loci. (See https://github.com/MEEIBioinformaticsCenter/grepsearch.) Results: In a test sample set consisting of eleven previously published homozygous mutants, detection of the MAK-Alu insertion was validated with 100% sensitivity and specificity. As a discovery cohort, raw NGS reads from 1,847 samples (including custom and whole exome selective capture) were searched in ~1 hour on a local computer cluster, yielding an additional five samples with MAK-Alu insertions and solving two previously unsolved pedigrees. Of these, one patient was homozygous for the insertion, one compound heterozygous with a missense change on the other allele (c. 46G>A; p.Gly16Arg), and three were heterozygous carriers. Conclusions: Using the MAK-Alu grep program proved to be a rapid and effective method of finding a known, disease-causing Alu insertion in a large cohort of patients with NGS data. This simple approach avoids wet-lab assays or computationally expensive algorithms, and could also be used for other known disease-causing insertions and deletions

Comprehensive Rare Variant Analysis via Whole-Genome Sequencing to Determine the Molecular Pathology of Inherited Retinal Disease

Inherited retinal disease is a common cause of visual impairment and represents a highly heterogeneous group of conditions. Here, we present findings from a cohort of 722 individuals with inherited retinal disease, who have had whole-genome sequencing (n = 605), whole-exome sequencing (n = 72), or both (n = 45) performed, as part of the NIHR-BioResource Rare Diseases research study. We identified pathogenic variants (single-nucleotide variants, indels, or structural variants) for 404/722 (56%) individuals. Whole-genome sequencing gives unprecedented power to detect three categories of pathogenic variants in particular: structural variants, variants in GC-rich regions, which have significantly improved coverage compared to whole-exome sequencing, and variants in non-coding regulatory regions. In addition to previously reported pathogenic regulatory variants, we have identified a previously unreported pathogenic intronic variant in $\textit{CHM}$ in two males with choroideremia. We have also identified 19 genes not previously known to be associated with inherited retinal disease, which harbor biallelic predicted protein-truncating variants in unsolved cases. Whole-genome sequencing is an increasingly important comprehensive method with which to investigate the genetic causes of inherited retinal disease.This work was supported by The National Institute for Health Research England (NIHR) for the NIHR BioResource – Rare Diseases project (grant number RG65966). The Moorfields Eye Hospital cohort of patients and clinical and imaging data were ascertained and collected with the support of grants from the National Institute for Health Research Biomedical Research Centre at Moorfields Eye Hospital, National Health Service Foundation Trust, and UCL Institute of Ophthalmology, Moorfields Eye Hospital Special Trustees, Moorfields Eye Charity, the Foundation Fighting Blindness (USA), and Retinitis Pigmentosa Fighting Blindness. M.M. is a recipient of an FFB Career Development Award. E.M. is supported by UCLH/UCL NIHR Biomedical Research Centre. F.L.R. and D.G. are supported by Cambridge NIHR Biomedical Research Centre

The Meckel-Gruber syndrome protein TMEM67 controls basal body positioning and epithelial branching morphogenesis in mice via the non-canonical Wnt pathway

Ciliopathies are a group of developmental disorders that manifest with multi-organ anomalies. Mutations in TMEM67 (MKS3) cause a range of human ciliopathies, including Meckel-Gruber and Joubert syndromes. In this study we describe multi-organ developmental abnormalities in the Tmem67(tm1Dgen/H1) knockout mouse that closely resemble those seen in Wnt5a and Ror2 knockout mice. These include pulmonary hypoplasia, ventricular septal defects, shortening of the body longitudinal axis, limb abnormalities, and cochlear hair cell stereociliary bundle orientation and basal body/kinocilium positioning defects. The basal body/kinocilium complex was often uncoupled from the hair bundle, suggesting aberrant basal body migration, although planar cell polarity and apical planar asymmetry in the organ of Corti were normal. TMEM67 (meckelin) is essential for phosphorylation of the non-canonical Wnt receptor ROR2 (receptor-tyrosine-kinase-like orphan receptor 2) upon stimulation with Wnt5a-conditioned medium. ROR2 also colocalises and interacts with TMEM67 at the ciliary transition zone. Additionally, the extracellular N-terminal domain of TMEM67 preferentially binds to Wnt5a in an in vitro binding assay. Cultured lungs of Tmem67 mutant mice failed to respond to stimulation of epithelial branching morphogenesis by Wnt5a. Wnt5a also inhibited both the Shh and canonical Wnt/β-catenin signalling pathways in wild-type embryonic lung. Pulmonary hypoplasia phenotypes, including loss of correct epithelial branching morphogenesis and cell polarity, were rescued by stimulating the non-canonical Wnt pathway downstream of the Wnt5a-TMEM67-ROR2 axis by activating RhoA. We propose that TMEM67 is a receptor that has a main role in non-canonical Wnt signalling, mediated by Wnt5a and ROR2, and normally represses Shh signalling. Downstream therapeutic targeting of the Wnt5a-TMEM67-ROR2 axis might, therefore, reduce or prevent pulmonary hypoplasia in ciliopathies and other congenital conditions

Clinical utility gene card for : inherited optic neuropathies including next-generation sequencing-based approaches

Non peer reviewe

Identification of novel mutations in Chinese Hans with autosomal dominant polycystic kidney disease

<p>Abstract</p> <p>Background</p> <p>Autosomal dominant polycystic kidney disease (ADPKD) is the most common inherited renal disease with an incidence of 1 in 400 to 1000. The disease is genetically heterogeneous, with two genes identified: <it>PKD1 </it>(16p13.3) and <it>PKD2 </it>(4q21). Molecular diagnosis of the disease in at-risk individuals is complicated due to the structural complexity of <it>PKD1 </it>gene and the high diversity of the mutations. This study is the first systematic ADPKD mutation analysis of both <it>PKD1 </it>and <it>PKD2 </it>genes in Chinese patients using denaturing high-performance liquid chromatography (DHPLC).</p> <p>Methods</p> <p>Both <it>PKD1 </it>and <it>PKD2 </it>genes were mutation screened in each proband from 65 families using DHPLC followed by DNA sequencing. Novel variations found in the probands were checked in their family members available and 100 unrelated normal controls. Then the pathogenic potential of the variations of unknown significance was examined by evolutionary comparison, effects of amino acid substitutions on protein structure, and effects of splice site alterations using online mutation prediction resources.</p> <p>Results</p> <p>A total of 92 variations were identified, including 27 reported previously. Definitely pathogenic mutations (ten frameshift, ten nonsense, two splicing defects and one duplication) were identified in 28 families, and probably pathogenic mutations were found in an additional six families, giving a total detection level of 52.3% (34/65). About 69% (20/29) of the mutations are first reported with a recurrent mutation rate of 31%.</p> <p>Conclusions</p> <p>Mutation study of <it>PKD1 </it>and <it>PKD2 </it>genes in Chinese Hans with ADPKD may contribute to a better understanding of the genetic diversity between different ethnic groups and enrich the mutation database. Besides, evaluating the pathogenic potential of novel variations should also facilitate the clinical diagnosis and genetic counseling of the disease.</p

High Resolution Melt analysis for mutation screening in PKD1 and PKD2

<p>Abstract</p> <p>Background</p> <p>Autosomal dominant polycystic kidney disease (ADPKD) is the most common hereditary kidney disorder. It is characterized by focal development and progressive enlargement of renal cysts leading to end-stage renal disease. <it>PKD1 </it>and <it>PKD2 </it>have been implicated in ADPKD pathogenesis but genetic features and the size of <it>PKD1 </it>make genetic diagnosis tedious.</p> <p>Methods</p> <p>We aim to prove that high resolution melt analysis (HRM), a recent technique in molecular biology, can facilitate molecular diagnosis of ADPKD. We screened for mutations in <it>PKD1 </it>and <it>PKD2 </it>with HRM in 37 unrelated patients with ADPKD.</p> <p>Results</p> <p>We identified 440 sequence variants in the 37 patients. One hundred and thirty eight were different. We found 28 pathogenic mutations (25 in <it>PKD1 </it>and 3 in <it>PKD2 </it>) within 28 different patients, which is a diagnosis rate of 75% consistent with literature mean direct sequencing diagnosis rate. We describe 52 new sequence variants in <it>PKD1 </it>and two in <it>PKD2</it>.</p> <p>Conclusion</p> <p>HRM analysis is a sensitive and specific method for molecular diagnosis of ADPKD. HRM analysis is also costless and time sparing. Thus, this method is efficient and might be used for mutation pre-screening in ADPKD genes.</p