Michael H. Farkas

Comparative analysis of RNA-Seq alignment algorithms and the RNA-Seq unified mapper (RUM)

MOTIVATION A critical task in high-throughput sequencing is aligning millions of short reads to a reference genome. Alignment is especially complicated for RNA sequencing (RNA-Seq) because of RNA splicing. A number of RNA-Seq algorithms are available, and claim to align reads with high accuracy and efficiency while detecting splice junctions. RNA-Seq data are discrete in nature; therefore, with reasonable gene models and comparative metrics RNA-Seq data can be simulated to sufficient accuracy to enable meaningful benchmarking of alignment algorithms. The exercise to rigorously compare all viable published RNA-Seq algorithms has not been performed previously. RESULTS We developed an RNA-Seq simulator that models the main impediments to RNA alignment, including alternative splicing, insertions, deletions, substitutions, sequencing errors and intron signal. We used this simulator to measure the accuracy and robustness of available algorithms at the base and junction levels. Additionally, we used reverse transcription-polymerase chain reaction (RT-PCR) and Sanger sequencing to validate the ability of the algorithms to detect novel transcript features such as novel exons and alternative splicing in RNA-Seq data from mouse retina. A pipeline based on BLAT was developed to explore the performance of established tools for this problem, and to compare it to the recently developed methods. This pipeline, the RNA-Seq Unified Mapper (RUM), performs comparably to the best current aligners and provides an advantageous combination of accuracy, speed and usability. AVAILABILITY The RUM pipeline is distributed via the Amazon Cloud and for computing clusters using the Sun Grid Engine (http://cbil.upenn.edu/RUM). CONTACT ggrant@pcbi.upenn.edu; epierce@mail.med.upenn.edu SUPPLEMENTARY INFORMATION The RNA-Seq sequence reads described in the article are deposited at GEO, accession GSE26248.

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.

BackgroundThe 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.ResultsWe 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.ConclusionsTo 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.

Three Gene-Targeted Mouse Models of RNA Splicing Factor RP Show Late-Onset RPE and Retinal Degeneration

PURPOSE Mutations in genes that produce proteins involved in mRNA splicing, including pre-mRNA processing factors 3, 8, and 31 (PRPF3, 8, and 31), RP9, and SNRNP200 are common causes of the late-onset inherited blinding disorder retinitis pigmentosa (RP). It is not known how mutations in these ubiquitously expressed genes lead to retina-specific disease. To investigate the pathogenesis of the RNA splicing factor forms of RP, the authors generated and characterized the retinal phenotypes of Prpf3-T494M, Prpf8-H2309P knockin mice. The retinal ultrastructure of Prpf31-knockout mice was also investigated. METHODS The knockin mice have single codon alterations in their endogenous Prpf3 and Prpf8 genes that mimic the most common disease causing mutations in human PRPF3 and PRPF8. The Prpf31-knockout mice mimic the null alleles that result from the majority of mutations identified in PRPF31 patients. The retinal phenotypes of the gene targeted mice were evaluated by electroretinography (ERG), light, and electron microscopy. RESULTS The RPE cells of heterozygous Prpf3(+/T494M) and Prpf8(+/H2309P) knockin mice exhibited loss of the basal infoldings and vacuolization, with accumulation of amorphous deposits between the RPE and Bruch[b]s membrane at age two years. These changes were more severe in the homozygous mice, and were associated with decreased rod function in the Prpf3-T494M mice. Similar degenerative changes in the RPE were detected in Prpf31(±) mice at one year of age. CONCLUSIONS The finding of similar degenerative changes in RPE cells of all three mouse models suggests that the RPE may be the primary cell type affected in the RNA splicing factor forms of RP. The relatively late-onset phenotype observed in these mice is consistent with the typical adult onset of disease in patients with RP.

Development of a rapid biolistic assay to determine changes in relative levels of intracellular calcium in leaves following tetracycline uptake by pinto bean plants

Tetracycline antibiotics, such as chlortetracycline (CTC) and tetracycline (TC), are introduced into agricultural lands through the application of manure as fertilizer. These compounds are phytotoxic to certain crop plants, including pinto beans (Phaseolus vulgaris), the species used for this investigation. While the mechanism of this toxicity is not yet understood, CTC is known to be a calcium chelator. We describe here a novel method to show that CTC is taken up by pinto bean plants and chelates calcium in leaves. Cameleon fusion proteins can provide qualitative and quantitative imaging of intracellular calcium levels, but current methodology requires stable transformation. Many plant species, including pinto beans, are not yet transformable using standard Agrobacterium-based protocols. To determine the role of calcium chelation in this plant, a rapid, biolistic method was developed to transiently express the cameleon protein. This method can easily be adapted to other plant systems. Our findings provide evidence that chelation of intracellular calcium by CTC is related to phytotoxic effects caused by this antibiotic in pinto beans. Root uptake of CTC and TC by pinto beans and their translocation to leaves were further verified by fluorescence spectroscopy and liquid chromatography/mass spectrometry, confirming results of the biolistic method that showed calcium chelation by tetracyclines in leaves.

Mutations in Pre-mRNA Processing Factors 3, 8, and 31 Cause Dysfunction of the Retinal Pigment Epithelium

Mutations in the ubiquitously expressed pre-mRNA processing factors 3, 8, and 31 (PRPF3, PRPF8, and PRPF31) cause nonsyndromic dominant retinitis pigmentosa in humans, an inherited retinal degeneration. It is unclear what mechanisms, or which cell types of the retina, are affected. Transgenic mice with the human mutations in these genes display late-onset morphological changes in the retinal pigment epithelium (RPE). To determine whether the observed morphological changes are preceded by abnormal RPE function, we investigated its phagocytic function in Prpf3(T494M/T494M), Prpf8(H2309P/H2309P), and Prpf31(+/-) mice. We observe decreased phagocytosis in primary RPE cultures from mutant mice, and this is replicated by shRNA-mediated knockdown of PRPF31 in human ARPE-19 cells. The diurnal rhythmicity of phagocytosis is almost lost, indicated by the marked attenuation of the phagocytic burst 2 hours after light onset. The strength of adhesion between RPE apical microvilli and photoreceptor outer segments also declined during peak adhesion in all mutants. In all models, at least one of the receptors involved in binding and internalization of shed photoreceptor outer segments was subjected to changes in localization. Although the mechanism underlying these changes in RPE function is yet to be elucidated, these data are consistent with the mouse RPE being the primary cell affected by mutations in the RNA splicing factors, and these changes occur at an early age.

Targeted Exon Sequencing in Usher Syndrome Type I

PURPOSE Patients with Usher syndrome type I (USH1) have retinitis pigmentosa, profound congenital hearing loss, and vestibular ataxia. This syndrome is currently thought to be associated with at least six genes, which are encoded by over 180 exons. Here, we present the use of state-of-the-art techniques in the molecular diagnosis of a cohort of 47 USH1 probands. METHODS The cohort was studied with selective exon capture and next-generation sequencing of currently known inherited retinal degeneration genes, comparative genomic hybridization, and Sanger sequencing of new USH1 exons identified by human retinal transcriptome analysis. RESULTS With this approach, we were able to genetically solve 14 of the 47 probands by confirming the biallelic inheritance of mutations. We detected two likely pathogenic variants in an additional 19 patients, for whom family members were not available for cosegregation analysis to confirm biallelic inheritance. Ten patients, in addition to primary disease-causing mutations, carried rare likely pathogenic USH1 alleles or variants in other genes associated with deaf-blindness, which may influence disease phenotype. Twenty-one of the identified mutations were novel among the 33 definite or likely solved patients. Here, we also present a clinical description of the studied cohort at their initial visits. CONCLUSIONS We found a remarkable genetic heterogeneity in the studied USH1 cohort with multiplicity of mutations, of which many were novel. No obvious influence of genotype on phenotype was found, possibly due to small sample sizes of the genotypes under study.

Determination of enzyme kinetics and glutathione conjugates of chlortetracycline and chloroacetanilides using liquid chromatography–mass spectrometry

Glutathione S-transferases (GSTs) isolated from chlortetracycline (CTC)-treated maize catalyzed the conjugation of glutathione (GSH) with CTC, producing stable conjugates that were structurally characterized using liquid chromatography-ion-trap mass spectrometry (LC-IT-MS). Enzyme-mediated dechlorination of CTC resulted during GSH conjugation as revealed by the mass spectra of the CTC-GSH conjugate, which was characterized by the loss of the chlorine isotopic signature, and shorter chromatographic retention time relative to the chlorinated parent compound. Several fragmentation patterns in the mass spectrum of the CTC-GSH conjugate can be used to verify the identity of the enzyme reaction products. The expected molecular ion [M + H](+) of the CTC-GSH conjugate (m/z 751) with chlorine removal was not observed in the positive electrospray ionization. Instead, a base peak of m/z 677, corresponding to the loss of glycine (MW = 75 Da), was observed. When m/z 677 was subjected to further fragmentation, characteristic peaks corresponding to the loss of glutamic acid (m/z = 129) and water (m/z 18) were observed in the MS/MS spectrum. The catalytic activity of the CTC-induced GST towards dechlorination of chloroacetanilide herbicides (alachlor, metolachlor and propachlor), which are known to be detoxified in plants via the glutathione pathway, was also evaluated in vitro. Glutathione conjugates of chloroacetanilides also showed the losses of m/z 129 and m/z 18 that are characteristic of GSH conjugates when characterized by LC-IT-MS. Interestingly, the sensitivity of LC-IT-MS made it possible, for the first time, to detect chloroacetanilides that are conjugated with two GSH molecules, in addition to the known single GSH conjugates. This research demonstrates a more sensitive and specific method of measuring enzyme reaction products using LC-IT-MS.

Transcriptome Analyses to Investigate the Pathogenesis of RNA Splicing Factor Retinitis Pigmentosa

RNA-splicing factor retinitis pigmentosa (RP) is caused by mutations in components of the spliceosome. RP is an inherited blinding disorder characterized by late-onset retinal degeneration. Currently, mutations in five genes that encode components of the spliceosome have been identified to cause autosomal dominant RP. These are the pre-mRNA processing factors 3, 8, and 31 (PRPF3, 8, and 31), RP9, and SNRNP200. It is unknown how mutations in these ubiquitously expressed genes lead to retina-specific disease. It is hypothesized that mutations in these genes lead to aberrant splicing of pre-mRNA, which in turn causes retinal degeneration. To fully investigate this hypothesis requires the ability to accurately interrogate the transcriptomes of the affected tissue. The recent development of next-generation sequencing-based RNA sequencing (RNA-seq) makes these types of studies possible. This chapter will focus on the RNA splicing factor forms of RP and the application of RNA-seq to study the pathogenesis of these diseases.

RNA-Seq: Improving Our Understanding of Retinal Biology and Disease

Over the past several years, rapid technological advances have allowed for a dramatic increase in our knowledge and understanding of the transcriptional landscape, because of the ability to study gene expression in greater depth and with more detail than previously possible. To this end, RNA-Seq has quickly become one of the most widely used methods for studying transcriptomes of tissues and individual cells. Unlike previously favored analysis methods, RNA-Seq is extremely high-throughput, and is not dependent on an annotated transcriptome, laying the foundation for novel genetic discovery. Additionally, RNA-Seq derived transcriptomes provide a basis for widening the scope of research to identify potential targets in the treatment of retinal disease.

Disrupted alternative splicing for genes implicated in splicing and ciliogenesis causes PRPF31 retinitis pigmentosa

Mutations in pre-mRNA processing factors (PRPFs) cause autosomal-dominant retinitis pigmentosa (RP), but it is unclear why mutations in ubiquitously expressed genes cause non-syndromic retinal disease. Here, we generate transcriptome profiles from RP11 (PRPF31-mutated) patient-derived retinal organoids and retinal pigment epithelium (RPE), as well as Prpf31+/− mouse tissues, which revealed that disrupted alternative splicing occurred for specific splicing programmes. Mis-splicing of genes encoding pre-mRNA splicing proteins was limited to patient-specific retinal cells and Prpf31+/− mouse retinae and RPE. Mis-splicing of genes implicated in ciliogenesis and cellular adhesion was associated with severe RPE defects that include disrupted apical – basal polarity, reduced trans-epithelial resistance and phagocytic capacity, and decreased cilia length and incidence. Disrupted cilia morphology also occurred in patient-derived photoreceptors, associated with progressive degeneration and cellular stress. In situ gene editing of a pathogenic mutation rescued protein expression and key cellular phenotypes in RPE and photoreceptors, providing proof of concept for future therapeutic strategies.Mutations in pre-mRNA processing factors cause autosomal dominant retinitis pigmentosa. Here the authors provide insights into the pathophysiological mechanisms underlying non-syndromic retinal disease caused by heterozygous mutations in genes encoding ubiquitously expressed splicing factors.