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.

Related papers:

title: Cell populations of the retina: the Proctor lecture.
abstract: It is now possible to describe—imperfectly, but in more than outline—the types of neurons by which the retina carries out its computations. This is the rough draft of a “neurome,” a conscious effort to list all the tissue's neurons. Identifying the cell populations has led directly to a revised conception of the retina's input–output functions and of the fundamental elements from which visual perception is built. 
 
The retina is the first substantially complex structure of the central nervous system for which such a more or less complete parts list is available . The resulting picture, the work of a community of investigators worldwide, resulted in a paradigm shift. A retina once imagined as based around a 6- to 12-cell module was replaced by a retina of multiple parallel arrays, containing at least 60 functionally different cell types and decomposing the visual input into more than a dozen different representations. The multiplex structure raises fundamental questions about the events that must occur later, in the central visual system. 
 
 
 
Figure 1. 
 
The major cell populations of a generic mammalian retina, shown semi-schematically in section view. Mammalian retinas contain in excess of 60 distinct cell types, each carrying out a distinct, dedicated function. The populations of photoreceptor, horizontal, ... 
 
 
 
Here, I will give a global view of the retina's cell populations, aimed, like the Proctor Lecture itself, at an audience of nonspecialist vision scientists. I will review how the cell populations were revealed and how they participate in the retina's functions. Along the way, I will point out some of the places where our concepts of the retina were overturned and comment on the next generation of questions.

title: Development and neurogenic potential of Müller glial cells in the vertebrate retina
abstract: Considerable research on normal and diseased states within the retina has focused on neurons. Recent research on glia throughout the central nervous system, including within the retina where Müller glia are the main type of glia, has provided a more in depth view of glial functions in health and disease. Glial cells have been recognized as being vital for the maintenance of a healthy tissue environment, where they actively participate in neuronal activity. More recently, Müller glia have been recognized as being very similar to retinal progenitor cells, particularly when compared at the molecular level using comprehensive expression profiling techniques. The molecular similarities, as well as the developmental events that occur at the end of the genesis period of retinal cells, have led us to propose that Müller glia are a form of late stage retinal progenitor cells. These late stage progenitor cells acquire some specialized glial functions, but do not irreversibly leave the progenitor state. Indeed, Müller glia appear to be able to behave as a progenitor in that they have been shown to proliferate and produce neurons in several instances when an acute injury has been applied to the retina. Enhancement of this response is thus an exciting strategy for retinal repair.

title: The genomic, biochemical, and cellular responses of the retina in inherited photoreceptor degenerations and prospects for the treatment of these disorders.
abstract: The association of more than 140 genes with human photoreceptor degenerations, together with studies of animal models of these monogenic diseases, has provided great insight into their pathogenesis. Here we review the responses of the retina to photoreceptor mutations, including mechanisms of photoreceptor death. We discuss the roles of oxidative metabolism, mitochondrial reactive oxygen species, metabolic stress, protein misfolding, and defects in ciliary proteins, as well as the responses of Müller glia, microglia, and the retinal vasculature. Finally, we report on potential pharmacologic and biologic therapies, the critical role of histopathology as a prerequisite to treatment, and the exciting promise of gene therapy in animal models and in phase 1 trials in humans.

title: AAV2 Gene Therapy Readministration in Three Adults with Congenital Blindness
abstract: Repeat administration of gene therapy to the contralateral retina of three congenitally blind patients was safe and resulted in improved vision. Shining a Light with Gene Therapy Gene therapy has great potential for treating certain diseases by providing therapeutic genes to target cells. Administration of a gene therapy vector carrying the RPE65 gene in 12 patients with congenital blindness due to RPE65 mutations led to improvements in retinal and visual function and proved to be a safe and stable procedure. In a follow-up study, the same group of researchers led by Jean Bennett set out to discover whether it would be possible to safely administer the vector and the therapeutic transgene to the contralateral eye of the patients. A big concern was whether the first gene therapy injection might have primed the patients’ immune system to respond to the adeno-associated virus (AAV) vector or the product of the therapeutic transgene that it had delivered. To test the safety and efficacy of a second administration of gene therapy to the second eye, the authors demonstrated that readministration was both safe and effective in animal models. Then, they selected 3 of the original 12 patients and readministered the AAV vector and its RPE65 transgene to the contralateral eye. They assessed safety by evaluating inflammatory responses, immune reactions, and extraocular exposure to the AAV vector. Efficacy was assessed through qualitative and quantitative measures of retinal and visual function including the ability to read letters, the extent of side vision, light sensitivity, the pupillary light reflex, the ability to navigate in dim light, and evidence from neuroimaging studies of cortical activation (which demonstrated that signals from the retina were recognized by the brain). The researchers did not discover any safety concerns and did not identify harmful immune responses to the vector or the transgene product. Before and after comparisons of psychophysical data and cortical responses provided the authors with evidence that gene therapy readministration was effective and mediated improvements in retinal and visual function in the three patients. The researchers report that the lack of immune response and the robust safety profile in this readministration gene therapy study may be due in part to the immune-privileged nature of the eye, and the low dose and very pure preparation of the AAV vector. Demonstration of safe and stable reversal of blindness after a single unilateral subretinal injection of a recombinant adeno-associated virus (AAV) carrying the RPE65 gene (AAV2-hRPE65v2) prompted us to determine whether it was possible to obtain additional benefit through a second administration of the AAV vector to the contralateral eye. Readministration of vector to the second eye was carried out in three adults with Leber congenital amaurosis due to mutations in the RPE65 gene 1.7 to 3.3 years after they had received their initial subretinal injection of AAV2-hRPE65v2.  (through 6 months) including evaluations of immune response, retinal and visual function testing, and functional magnetic resonance imaging indicate that readministration is both safe and efficacious after previous exposure to AAV2-hRPE65v2.

title: Safety and efficacy of gene transfer for Leber's congenital amaurosis.
abstract: Leber's congenital amaurosis (LCA) is a group of inherited blinding diseases with onset during childhood. One form of the disease, LCA2, is caused by mutations in the retinal pigment epithelium-specific 65-kDa protein gene (RPE65). We investigated the safety of subretinal delivery of a recombinant adeno-associated virus (AAV) carrying RPE65 complementary DNA (cDNA) (ClinicalTrials.gov number, NCT00516477 . Three patients with LCA2 had an acceptable local and systemic adverse-event profile after delivery of AAV2.hRPE65v2. Each patient had a modest improvement in measures of retinal function on subjective tests of visual acuity. In one patient, an asymptomatic macular hole developed, and although the occurrence was considered to be an adverse event, the patient had some return of retinal function. Although the follow-up was very short and normal vision was not achieved, this study provides the basis for further gene therapy studies in patients with LCA.

title: Effect of gene therapy on visual function in Leber's congenital amaurosis.
abstract: Early-onset, severe retinal dystrophy caused by mutations in the gene encoding retinal pigment epithelium-specific 65-kD protein (RPE65) is associated with poor vision at birth and complete loss of vision in early adulthood. We administered to three young adult patients subretinal injections of recombinant adeno-associated virus vector 2/2 expressing RPE65 complementary DNA (cDNA) under the control of a human RPE65 promoter. There were no serious adverse events. There was no clinically significant change in visual acuity or in peripheral visual fields on Goldmann perimetry in any of the three patients. We detected no change in retinal responses on electroretinography. One patient had significant improvement in visual function on microperimetry and on dark-adapted perimetry. This patient also showed improvement in a subjective test of visual mobility. These findings provide support for further clinical studies of this experimental approach in other patients with mutant RPE65.

title: Complete exon sequencing of all known Usher syndrome genes greatly improves molecular diagnosis
abstract: BackgroundUsher syndrome (USH) combines sensorineural deafness with blindness. It is inherited in an autosomal recessive mode. Early diagnosis is critical for adapted educational and patient management choices, and for genetic counseling. To date, nine causative genes have been identified for the three clinical subtypes (USH1, USH2 and USH3). Current diagnostic strategies make use of a genotyping microarray that is based on the previously reported mutations. The purpose of this study was to design a more accurate molecular diagnosis tool.MethodsWe sequenced the 366 coding exons and flanking regions of the nine known USH genes, in 54 USH patients (27 USH1, 21 USH2 and 6 USH3).ResultsBiallelic mutations were detected in 39 patients (72%) and monoallelic mutations in an additional 10 patients (18.5%). In addition to biallelic mutations in one of the USH genes, presumably pathogenic mutations in another USH gene were detected in seven patients (13%), and another patient carried monoallelic mutations in three different USH genes. Notably, none of the USH3 patients carried detectable mutations in the only known USH3 gene, whereas they all carried mutations in USH2 genes. Most importantly, the currently used microarray would have detected only 30 of the 81 different mutations that we found, of which 39 (48%) were novel.ConclusionsBased on these results, complete exon sequencing of the currently known USH genes stands as a definite improvement for molecular diagnosis of this disease, which is of utmost importance in the perspective of gene therapy.

title: A splice-site mutation in a retina-specific exon of BBS8 causes nonsyndromic retinitis pigmentosa.
abstract: Tissue-specific alternative splicing is an important mechanism for providing spatiotemporal protein diversity. Here we show that an in-frame splice mutation in BBS8, one of the genes involved in pleiotropic Bardet-Biedl syndrome (BBS), is sufficient to cause nonsyndromic retinitis pigmentosa (RP). A genome-wide scan of a consanguineous RP pedigree mapped the trait to a 5.6 Mb region; subsequent systematic sequencing of candidate transcripts identified a homozygous splice-site mutation in a previously unknown BBS8 exon. The allele segregated with the disorder, was absent from controls, was completely invariant across evolution, and was predicted to lead to the elimination of a 10 amino acid sequence from the protein. Subsequent studies showed the exon to be expressed exclusively in the retina and enriched significantly in the photoreceptor layer. Importantly, we found this exon to represent the major BBS8 mRNA species in the mammalian photoreceptor, suggesting that the encoded 10 amino acids play a pivotal role in the function of BBS8 in this organ. Understanding the role of this additional sequence might therefore inform the mechanism of retinal degeneration in patients with syndromic BBS or other related ciliopathies.

title: RPGR transcription studies in mouse and human tissues reveal a retina-specific isoform that is disrupted in a patient with X-linked retinitis pigmentosa.
abstract: X-linked retinitis pigmentosa (XLRP) is a genetically heterogeneous group of progressive retinal degenerations. The disease process is initiated by premature apoptosis of rod photoreceptor cells in the retina, which leads to reduced visual acuity and, eventually, complete blindness. Mutations in the retinitis pigmentosa GTPase regulator ( RPGR ), a ubiquitously expressed gene at the RP3 locus in Xp21.1, account for approximately 20% of all X-linked cases. We have analysed the expression of this gene by northern blot hybridization, cDNA library screening and RT-PCR in various organs from mouse and man. These studies revealed at least 12 alternatively spliced isoforms. Some of the transcripts are tissue specific and contain novel exons, which elongate or truncate the previously reported open reading frame of the mouse and human RPGR gene. One of the newly identified exons is expressed exclusively in the human retina and mouse eye and contains a premature stop codon. The deduced polypeptide lacks 169 amino acids from the C-terminus of the ubiquitously expressed variant, including an isoprenylation site. Moreover, this exon was found to be deleted in a family with XLRP. Our results indicate tissue-dependent regulation of alternative splicing of RPGR in mouse and man. The discovery of a retina-specific transcript may explain why phenotypic abberations in RP3 are confined to the eye.

title: Pleiotropic effects of CEP290 (NPHP6) mutations extend to Meckel syndrome.
abstract: Meckel syndrome (MKS) is a rare autosomal recessive lethal condition characterized by central nervous system malformations, polydactyly, multicystic kidney dysplasia, and ductal changes of the liver. Three loci have been mapped (MKS1-MKS3), and two genes have been identified (MKS1/FLJ20345 and MKS3/TMEM67), whereas the gene at the MKS2 locus remains unknown. To identify new MKS loci, a genomewide linkage scan was performed using 10-cM-resolution microsatellite markers in eight families. The highest heterogeneity LOD score was obtained for chromosome 12, in an interval containing CEP290, a gene recently identified as causative of Joubert syndrome (JS) and isolated Leber congenital amaurosis. In view of our recent findings of allelism, at the MKS3 locus, between these two disorders, CEP290 was considered a candidate, and homozygous or compound heterozygous truncating mutations were identified in four families. Sequencing of additional cases identified CEP290 mutations in two fetuses with MKS and in four families presenting a cerebro-reno-digital syndrome, with a phenotype overlapping MKS and JS, further demonstrating that MKS and JS can be variable expressions of the same ciliopathy. These data identify a fourth locus for MKS (MKS4) and the CEP290 gene as responsible for MKS.

title: Mechanistic insights into Bardet-Biedl syndrome, a model ciliopathy.
abstract: Bardet-Biedl syndrome (BBS) is a multisystemic disorder typified by developmental and progressive degenerative defects. A combination of genetic, in vitro, and in vivo studies have highlighted ciliary dysfunction as a primary cause of BBS pathology, which has in turn contributed to the improved understanding of the functions of the primary cilium in humans and other vertebrates. Here we discuss the evidence linking the clinical BBS phenotype to ciliary defects, highlight how the genetic and cellular characteristics of BBS overlap with and inform other ciliary disorders, and explore the possible mechanistic underpinnings of ciliary dysfunction.

title: Retinal dystrophy in Bardet–Biedl syndrome and related syndromic ciliopathies
abstract: Primary cilia are almost ubiquitously expressed in eukaryotic cells where they function as sensors relaying information either from the extracellular environment or between two compartments of the same cell, such as in the photoreceptor cell. In ciliopathies, a continuously growing class of genetic disorders related to ciliary defects, the modified primary cilium of the photoreceptor, also known as the connecting cilium, is frequently defective. Ciliary dysfunction involves disturbances in the trafficking and docking of specific proteins involved in its biogenesis or maintenance. The main well-conserved ciliary process, intraflagellar transport (IFT), is a complex process carried out by multimeric ciliary particles and molecular motors of major importance in the photoreceptor cell. It is defective in a growing number of ciliopathies leading to retinal degeneration. Retinitis pigmentosa related to ciliary dysfunction can be an isolated feature or a part of a syndrome such as Bardet-Biedl syndrome (BBS). Research on ciliopathies and BBS has led to the discovery of several major cellular processes carried out by the primary cilium structure and has highlighted their genetic heterogeneity.

title: Novel RPGR-ORF15 mutations in X-linked retinitis pigmentosa patients
abstract: X-linked retinitis pigmentosa (XLRP) is the most severe type of retinitis pigmentosa (RP), with patients consistently showing early onset and rapid deterioration. Obtaining a genetic diagnosis for a family with XLRP is important for counseling purposes. In this study, we aimed to identify disease-causing mutations in two unrelated XLRP families. Genetic analysis was performed on two unrelated XLRP families. Genomic DNA was extracted from peripheral blood or amniotic fluid samples. The coding regions and intron/exon boundaries of the Retinitis Pigmentosa GTPase Regulator (RPGR) and RP2 genes were amplified by PCR and then sequenced directly. A clinically unaffected pregnant female and the four month old fetus were found to have a hemizygous 2 base pair deletion (g.ORF15+484_485delAA) in the exon ORF15 of RPGR gene. In another XLRP family, a nonsense mutation (g.ORF15+810G>T) was identified. Neither mutation has been reported previously. Both are predicted to cause premature termination of the protein. In conclusion, we identified a micro-deletion through prenatal genetic diagnosis and another novel nonsense mutation in RPGR-ORF15. Identifying a disease-causing mutation facilitated early diagnosis and genetic counseling for the patients. Discovery of novel mutations also broadens knowledge of XLRP and the spectrum of its pathogenic genotypes.

title: RNA-sequence analysis of human B-cells.
abstract: RNA-sequencing (RNA-seq) allows quantitative measurement of expression levels of genes and their transcripts. In this study, we sequenced complementary DNA fragments of cultured human B-cells and obtained 879 million 50-bp reads comprising 44 Gb of sequence. The results allowed us to study the gene expression profile of B-cells and to determine experimental parameters for sequencing-based expression studies. We identified 20,766 genes and 67,453 of their alternatively spliced transcripts. More than 90% of the genes with multiple exons are alternatively spliced; for most genes, one isoform is predominantly expressed. We found that while chromosomes differ in gene density, the percentage of transcribed genes in each chromosome is less variable. In addition, genes involved in related biological processes are expressed at more similar levels than genes with different functions. Besides characterizing gene expression, we also used the data to investigate the effect of sequencing depth on gene expression measurements. While 100 million reads are sufficient to detect most expressed genes and transcripts, about 500 million reads are needed to measure accurately their expression levels. We provide examples in which deep sequencing is needed to determine the relative abundance of genes and their isoforms. With data from 20 individuals and about 40 million sequence reads per sample, we uncovered only 21 alternatively spliced, multi-exon genes that are not in databases; this result suggests that at this sequence coverage, we can detect most of the known genes. Results from this project are available on the UCSC Genome Browser to allow readers to study the expression and structure of genes in human B-cells.

title: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation.
abstract: High-throughput mRNA sequencing (RNA-Seq) promises simultaneous transcript discovery and abundance estimation. However, this would require algorithms that are not restricted by prior gene annotations and that account for alternative transcription and splicing. Here we introduce such algorithms in an open-source software program called Cufflinks. To test Cufflinks, we sequenced and analyzed >430 million paired 75-bp RNA-Seq reads from a mouse myoblast cell line over a differentiation time series. We detected 13,692 known transcripts and 3,724 previously unannotated ones, 62% of which are supported by independent expression data or by homologous genes in other species. Over the time series, 330 genes showed complete switches in the dominant transcription start site (TSS) or splice isoform, and we observed more subtle shifts in 1,304 other genes. These results suggest that Cufflinks can illuminate the substantial regulatory flexibility and complexity in even this well-studied model of muscle development and that it can improve transcriptome-based genome annotation.

title: A Global View of Gene Activity and Alternative Splicing by Deep Sequencing of the Human Transcriptome
abstract: The functional complexity of the human transcriptome is not yet fully elucidated. We report a high-throughput sequence of the human transcriptome from a human embryonic kidney and a B cell line. We used shotgun sequencing of transcripts to generate randomly distributed reads. Of these, 50% mapped to unique genomic locations, of which 80% corresponded to known exons. We found that 66% of the polyadenylated transcriptome mapped to known genes and 34% to nonannotated genomic regions. On the basis of known transcripts, RNA-Seq can detect 25% more genes than can microarrays. A global survey of messenger RNA splicing events identified 94,241 splice junctions (4096 of which were previously unidentified) and showed that exon skipping is the most prevalent form of alternative splicing.

title: Understanding mechanisms underlying human gene expression variation with RNA sequencing
abstract: Understanding the genetic mechanisms underlying natural variation in gene expression is a central goal of both medical and evolutionary genetics, and studies of expression quantitative trait loci (eQTLs) have become an important tool for achieving this goal. Although all eQTL studies so far have assayed messenger RNA levels using expression microarrays, recent advances in RNA sequencing enable the analysis of transcript variation at unprecedented resolution. We sequenced RNA from 69 lymphoblastoid cell lines derived from unrelated Nigerian individuals that have been extensively genotyped by the International HapMap Project. By pooling data from all individuals, we generated a map of the transcriptional landscape of these cells, identifying extensive use of unannotated untranslated regions and more than 100 new putative protein-coding exons. Using the genotypes from the HapMap project, we identified more than a thousand genes at which genetic variation influences overall expression levels or splicing. We demonstrate that eQTLs near genes generally act by a mechanism involving allele-specific expression, and that variation that influences the inclusion of an exon is enriched within and near the consensus splice sites. Our results illustrate the power of high-throughput sequencing for the joint analysis of variation in transcription, splicing and allele-specific expression across individuals.

title: Genome-wide RNA-seq analysis of human and mouse platelet transcriptomes.
abstract: Inbred mice are a useful tool for studying the in vivo functions of platelets. Nonetheless, the mRNA signature of mouse platelets is not known. Here, we use paired-end next-generation RNA sequencing (RNA-seq) to characterize the polyadenylated transcriptomes of human and mouse platelets. We report that RNA-seq provides unprecedented resolution of mRNAs that are expressed across the entire human and mouse genomes. Transcript expression and abundance are often conserved between the 2 species. Several mRNAs, however, are differentially expressed in human and mouse platelets. Moreover, previously described functional disparities between mouse and human platelets are reflected in differences at the transcript level, including protease activated receptor-1, protease activated receptor-3, platelet activating factor receptor, and factor V. This suggests that RNA-seq is a useful tool for predicting differences in platelet function between mice and humans. Our next-generation sequencing analysis provides new insights into the human and murine platelet transcriptomes. The sequencing dataset will be useful in the design of mouse models of hemostasis and a catalyst for discovery of new functions of platelets. Access to the dataset is found in the "

title: Transcriptome landscape of the human placenta
abstract: BackgroundThe placenta is a key component in understanding the physiological processes involved in pregnancy. Characterizing genes critical for placental function can serve as a basis for identifying mechanisms underlying both normal and pathologic pregnancies. Detailing the placental tissue transcriptome could provide a valuable resource for genomic studies related to placental disease.ResultsWe have conducted a deep RNA sequencing (RNA-Seq) study on three tissue components (amnion, chorion, and decidua) of 5 human placentas from normal term pregnancies. We compared the placental RNA-Seq data to that of 16 other human tissues and observed a wide spectrum of transcriptome differences both between placenta and other human tissues and between distinct compartments of the placenta. Exon-level analysis of the RNA-Seq data revealed a large number of exons with differential splicing activities between placenta and other tissues, and 79% (27 out of 34) of the events selected for RT-PCR test were validated. The master splicing regulator ESRP1 is expressed at a proportionately higher level in amnion compared to all other analyzed human tissues, and there is a significant enrichment of ESRP1-regulated exons with tissue-specific splicing activities in amnion. This suggests an important role of alternative splicing in regulating gene function and activity in specific placental compartments. Importantly, genes with differential expression or splicing in the placenta are significantly enriched for genes implicated in placental abnormalities and preterm birth. In addition, we identified 604-1007 novel transcripts and 494-585 novel exons expressed in each of the three placental compartments.ConclusionsOur data demonstrate unique aspects of gene expression and splicing in placental tissues that provide a basis for disease investigation related to disruption of these mechanisms. These data are publicly available providing the community with a rich resource for placental physiology and disease-related studies.

title: Integrative annotation of human large intergenic noncoding RNAs reveals global properties and specific subclasses.
abstract: Large intergenic noncoding RNAs (lincRNAs) are emerging as key regulators of diverse cellular processes. Determining the function of individual lincRNAs remains a challenge. Recent advances in RNA sequencing (RNA-seq) and computational methods allow for an unprecedented analysis of such transcripts. Here, we present an integrative approach to define a reference catalog of >8000 human lincRNAs. Our catalog unifies previously existing annotation sources with transcripts we assembled from RNA-seq data collected from ∼4 billion RNA-seq reads across 24 tissues and cell types. We characterize each lincRNA by a panorama of >30 properties, including sequence, structural, transcriptional, and orthology features. We found that lincRNA expression is strikingly tissue-specific compared with coding genes, and that lincRNAs are typically coexpressed with their neighboring genes, albeit to an extent similar to that of pairs of neighboring protein-coding genes. We distinguish an additional subset of transcripts that have high evolutionary conservation but may include short ORFs and may serve as either lincRNAs or small peptides. Our integrated, comprehensive, yet conservative reference catalog of human lincRNAs reveals the global properties of lincRNAs and will facilitate experimental studies and further functional classification of these genes.

title: Conserved Function of lincRNAs in Vertebrate Embryonic Development despite Rapid Sequence Evolution
abstract: Thousands of long intervening noncoding RNAs (lincRNAs) have been identified in mammals. To better understand the evolution and functions of these enigmatic RNAs, we used chromatin marks, poly(A)-site mapping and RNA-Seq data to identify more than 550 distinct lincRNAs in zebrafish. Although these shared many characteristics with mammalian lincRNAs, only 29 had detectable sequence similarity with putative mammalian orthologs, typically restricted to a single short region of high conservation. Other lincRNAs had conserved genomic locations without detectable sequence conservation. Antisense reagents targeting conserved regions of two zebrafish lincRNAs caused developmental defects. Reagents targeting splice sites caused the same defects and were rescued by adding either the mature lincRNA or its human or mouse ortholog. Our study provides a roadmap for identification and analysis of lincRNAs in model organisms and shows that lincRNAs play crucial biological roles during embryonic development with functionality conserved despite limited sequence conservation.

title: The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression.
abstract: The human genome contains many thousands of long noncoding RNAs (lncRNAs). While several studies have demonstrated compelling biological and disease roles for individual examples, analytical and experimental approaches to investigate these genes have been hampered by the lack of comprehensive lncRNA annotation. Here, we present and analyze the most complete human lncRNA annotation to date, produced by the GENCODE consortium within the framework of the ENCODE project and comprising 9277 manually annotated genes producing 14,880 transcripts. Our analyses indicate that lncRNAs are generated through pathways similar to that of protein-coding genes, with similar histone-modification profiles, splicing signals, and exon/intron lengths. In contrast to protein-coding genes, however, lncRNAs display a striking bias toward two-exon transcripts, they are predominantly localized in the chromatin and nucleus, and a fraction appear to be preferentially processed into small RNAs. They are under stronger selective pressure than neutrally evolving sequences-particularly in their promoter regions, which display levels of selection comparable to protein-coding genes. Importantly, about one-third seem to have arisen within the primate lineage. Comprehensive analysis of their expression in multiple human organs and brain regions shows that lncRNAs are generally lower expressed than protein-coding genes, and display more tissue-specific expression patterns, with a large fraction of tissue-specific lncRNAs expressed in the brain. Expression correlation analysis indicates that lncRNAs show particularly striking positive correlation with the expression of antisense coding genes. This GENCODE annotation represents a valuable resource for future studies of lncRNAs.

title: A novel resource for the study of genes expressed in the adult human retina.
abstract: To reconstruct the transcriptional profile of the human adult retina and the genomic map of the genes expressed in this tissue.



Original software was used for the retrieval and analysis of records from UniGene (http://www.ncbi.nlm.nih. gov/UniGene/) pertaining to selected cDNA libraries from adult human retina.



The 4974 genes reported so far to be expressed in retina were included in a catalog available on the Internet. For each entry, an estimation of the level of expression of the corresponding gene in the retina was provided. A high-resolution genomic map of the human retina was built up by inclusion of 3152 genes showing a precise and unique map assignment. The correspondence was established between 53 gene-orphan retinal diseases and clusters of genes expressed in the retina.



The in silico reconstruction of the transcriptional profile of the adult human retina provides preliminary information on the pattern of genomic expression in this tissue. The chromosomal location of many retinal genes, combined with their expression data, should speed up the identification of genes involved in retinal diseases.

title: Transcriptome analysis of the retina
abstract: The retina offers unique opportunities to define the molecular and cellular pathways mediating neuronal function and disease because of its morphological complexity, well-defined role in visual transduction and the availability of mutants. These investigations are being greatly facilitated by the ongoing identification of genes expressed in the retina using high-throughput methods.

title: The Retinome – Defining a reference transcriptome of the adult mammalian retina/retinal pigment epithelium
abstract: BackgroundThe mammalian retina is a valuable model system to study neuronal biology in health and disease. To obtain insight into intrinsic processes of the retina, great efforts are directed towards the identification and characterization of transcripts with functional relevance to this tissue.ResultsWith the goal to assemble a first genome-wide reference transcriptome of the adult mammalian retina, referred to as the retinome, we have extracted 13,037 non-redundant annotated genes from nearly 500,000 published datasets on redundant retina/retinal pigment epithelium (RPE) transcripts. The data were generated from 27 independent studies employing a wide range of molecular and biocomputational approaches. Comparison to known retina-/RPE-specific pathways and established retinal gene networks suggest that the reference retinome may represent up to 90% of the retinal transcripts. We show that the distribution of retinal genes along the chromosomes is not random but exhibits a higher order organization closely following the previously observed clustering of genes with increased expression.ConclusionThe genome wide retinome map offers a rational basis for selecting suggestive candidate genes for hereditary as well as complex retinal diseases facilitating elaborate studies into normal and pathological pathways. To make this unique resource freely available we have built a database providing a query interface to the reference retinome .

title: The NEIBank project for ocular genomics: Data-mining gene expression in human and rodent eye tissues
abstract: NEIBank is a project to gather and organize genomic resources for eye research. The first phase of this project covers the construction and sequence analysis of cDNA libraries from human and animal model eye tissues to develop an overview of the repertoire of genes expressed in the eye and a resource of cDNA clones for further studies. The sequence data are grouped and identified using the tools of bioinformatics and the results are displayed through a web site where they can be interrogated by keyword search, chromosome location, by Blast (sequence comparison) or by alignment on completed genomes. Many novel proteins and novel splice forms of known genes have already emerged from analysis of the accumulating data. This review provides an overview of the current state of the database for human eye tissues, with specific comparisons to some parallel data from mouse and rat, and with illustrative examples of the kinds of insights and discoveries these data can produce. One of the major themes that emerges is that at the molecular level human eye tissues have significant differences from those of rodents, encompassing species specific genes, alternative splice forms and great variation in levels of gene expression. These point to specific adaptations and mechanisms in the human eye and emphasize that care needs to be taken in the application of appropriate animal model systems.

title: Comprehensive Analysis of Photoreceptor Gene Expression and the Identification of Candidate Retinal Disease Genes
abstract: To identify the full set of genes expressed by mammalian rods, we conducted serial analysis of gene expression (SAGE) by using libraries generated from mature and developing mouse retina. We identified 264 uncharacterized genes that were specific to or highly enriched in rods. Nearly half of all cloned human retinal disease genes are selectively expressed in rod photoreceptors. In silico mapping of the human orthologs of genes identified in our screen revealed that 86 map within intervals containing uncloned retinal disease genes, representing 37 different loci. We expect these data will allow identification of many disease genes, and that this approach may be useful for cloning genes involved in classes of disease where cell type-specific expression of disease genes is observed.

title: Profile of the genes expressed in the human peripheral retina, macula, and retinal pigment epithelium determined through serial analysis of gene expression (SAGE)
abstract: We used the serial analysis of gene expression (SAGE) technique to catalogue and measure the relative levels of expression of the genes expressed in the human peripheral retina, macula, and retinal pigment epithelium (RPE) from one or both of two humans, aged 88 and 44 years. The cone photoreceptor contribution to all transcription in the retina was found to be similar in the macula versus the retinal periphery, whereas the rod contribution was greater in the periphery versus the macula. Genes encoding structural proteins for axons were found to be expressed at higher levels in the macula versus the retinal periphery, probably reflecting the large proportion of ganglion cells in the central retina. In comparison with the younger eye, the peripheral retina of the older eye had a substantially higher proportion of mRNAs from genes encoding proteins involved in iron metabolism or protection against oxidative damage and a substantially lower proportion of mRNAs from genes encoding proteins involved in rod phototransduction. These differences may reflect the difference in age between the two donors or merely interindividual variation. The RPE library had numerous previously unencountered tags, suggesting that this cell type has a large, idiosyncratic repertoire of expressed genes. Comparison of these libraries with 100 reported nonocular SAGE libraries revealed 89 retina-specific or enriched genes expressed at substantial levels, of which 14 are known to cause a retinal disease and 53 are RPE-specific genes. We expect that these libraries will serve as a resource for understanding the relative expression levels of genes in the retina and the RPE and for identifying additional disease genes.

title: Defining the human macula transcriptome and candidate retinal disease genes using EyeSAGE.
abstract: To develop large-scale, high-throughput annotation of the human macula transcriptome and to identify and prioritize candidate genes for inherited retinal dystrophies, based on ocular-expression profiles using serial analysis of gene expression (SAGE).



Two human retina and two retinal pigment epithelium (RPE)/choroid SAGE libraries made from matched macula or midperipheral retina and adjacent RPE/choroid of morphologically normal 28- to 66-year-old donors and a human central retina longSAGE library made from 41- to 66-year-old donors were generated. Their transcription profiles were entered into a relational database, EyeSAGE, including microarray expression profiles of retina and publicly available normal human tissue SAGE libraries. EyeSAGE was used to identify retina- and RPE-specific and -associated genes, and candidate genes for retina and RPE disease loci. Differential and/or cell-type specific expression was validated by quantitative and single-cell RT-PCR.



Cone photoreceptor-associated gene expression was elevated in the macula transcription profiles. Analysis of the longSAGE retina tags enhanced tag-to-gene mapping and revealed alternatively spliced genes. Analysis of candidate gene expression tables for the identified Bardet-Biedl syndrome disease gene (BBS5) in the BBS5 disease region table yielded BBS5 as the top candidate. Compelling candidates for inherited retina diseases were identified.



The EyeSAGE database, combining three different gene-profiling platforms including the authors' multidonor-derived retina/RPE SAGE libraries and existing single-donor retina/RPE libraries, is a powerful resource for definition of the retina and RPE transcriptomes. It can be used to identify retina-specific genes, including alternatively spliced transcripts and to prioritize candidate genes within mapped retinal disease regions.

title: Evaluation of statistical methods for normalization and differential expression in mRNA-Seq experiments
abstract: BackgroundHigh-throughput sequencing technologies, such as the Illumina Genome Analyzer, are powerful new tools for investigating a wide range of biological and medical questions. Statistical and computational methods are key for drawing meaningful and accurate conclusions from the massive and complex datasets generated by the sequencers. We provide a detailed evaluation of statistical methods for normalization and differential expression (DE) analysis of Illumina transcriptome sequencing (mRNA-Seq) data.ResultsWe compare statistical methods for detecting genes that are significantly DE between two types of biological samples and find that there are substantial differences in how the test statistics handle low-count genes. We evaluate how DE results are affected by features of the sequencing platform, such as, varying gene lengths, base-calling calibration method (with and without phi X control lane), and flow-cell/library preparation effects. We investigate the impact of the read count normalization method on DE results and show that the standard approach of scaling by total lane counts (e.g., RPKM) can bias estimates of DE. We propose more general quantile-based normalization procedures and demonstrate an improvement in DE detection.ConclusionsOur results have significant practical and methodological implications for the design and analysis of mRNA-Seq experiments. They highlight the importance of appropriate statistical methods for normalization and DE inference, to account for features of the sequencing platform that could impact the accuracy of results. They also reveal the need for further research in the development of statistical and computational methods for mRNA-Seq.

title: Comparative analysis of RNA-Seq alignment algorithms and the RNA-Seq unified mapper (RUM)
abstract: 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.



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).



ggrant@pcbi.upenn.edu; epierce@mail.med.upenn.edu



The RNA-Seq sequence reads described in the article are deposited at GEO, accession GSE26248.

title: Gene and noncoding RNA regulation underlying photoreceptor protection: microarray study of dietary antioxidant saffron and photobiomodulation in rat retina
abstract: To identify the genes and noncoding RNAs (ncRNAs) involved in the neuroprotective actions of a dietary antioxidant (saffron) and of photobiomodulation (PBM).  We used a previously published assay of photoreceptor damage, in which albino Sprague Dawley rats raised in dim cyclic illumination (12 h 5 lux, 12 h darkness) were challenged by 24 h exposure to bright (1,000 lux) light. Experimental groups were protected against light damage by pretreatment with dietary saffron (1 mg/kg/day for 21 days) or PBM (9 J/cm2 at the eye, daily for 5 days). RNA from one eye of four animals in each of the six experimental groups (control, light damage [LD], saffron, PBM, saffronLD, and PBMLD) was hybridized to Affymetrix rat genome ST arrays. Quantitative real-time PCR analysis of 14 selected genes was used to validate the microarray results.  LD caused the regulation of 175 entities (genes and ncRNAs) beyond criterion levels (p<0.05 in comparison with controls, fold-change >2). PBM pretreatment reduced the expression of 126 of these 175 LD-regulated entities below criterion; saffron pretreatment reduced the expression of 53 entities (50 in common with PBM). In addition, PBM pretreatment regulated the expression of 67 entities not regulated by LD, while saffron pretreatment regulated 122 entities not regulated by LD (48 in common with PBM). PBM and saffron, given without LD, regulated genes and ncRNAs beyond criterion levels, but in lesser numbers than during their protective action. A high proportion of the entities regulated by LD (>90%) were known genes. By contrast, ncRNAs were prominent among the entities regulated by PBM and saffron in their neuroprotective roles (73% and 62%, respectively).  Given alone, saffron and (more prominently) PBM both regulated significant numbers of genes and ncRNAs. Given before retinal exposure to damaging light, thus while exerting their neuroprotective action, they regulated much larger numbers of entities, among which ncRNAs were prominent. Further, the downregulation of known genes and of ncRNAs was prominent in the protective actions of both neuroprotectants. These comparisons provide an overview of gene expression induced by two neuroprotectants and provide a basis for the more focused study of their mechanisms.

title: The Ensembl automatic gene annotation system.
abstract: As more genomes are sequenced, there is an increasing need for automated first-pass annotation which allows timely access to important genomic information. The Ensembl gene-building system enables fast automated annotation of eukaryotic genomes. It annotates genes based on evidence derived from known protein, cDNA, and EST sequences. The gene-building system rests on top of the core Ensembl (MySQL) database schema and Perl Application Programming Interface (API), and the data generated are accessible through the Ensembl genome browser (http://www.ensembl.org). To date, the Ensembl predicted gene sets are available for the A. gambiae, C. briggsae, zebrafish, mouse, rat, and human genomes and have been heavily relied upon in the publication of the human, mouse, rat, and A. gambiae genome sequence analysis. Here we describe in detail the gene-building system and the algorithms involved. All code and data are freely available from http://www.ensembl.org.

title: AceView: a comprehensive cDNA-supported gene and transcripts annotation
abstract: BackgroundRegions covering one percent of the genome, selected by ENCODE for extensive analysis, were annotated by the HAVANA/Gencode group with high quality transcripts, thus defining a benchmark. The ENCODE Genome Annotation Assessment Project (EGASP) competition aimed at reproducing Gencode and finding new genes. The organizers evaluated the protein predictions in depth. We present a complementary analysis of the mRNAs, including alternative transcript variants.ResultsWe evaluate 25 gene tracks from the University of California Santa Cruz (UCSC) genome browser. We either distinguish or collapse the alternative splice variants, and compare the genomic coordinates of exons, introns and nucleotides. Whole mRNA models, seen as chains of introns, are sorted to find the best matching pairs, and compared so that each mRNA is used only once. At the mRNA level, AceView is by far the closest to Gencode: the vast majority of transcripts of the two methods, including alternative variants, are identical. At the protein level, however, due to a lack of experimental data, our predictions differ: Gencode annotates proteins in only 41% of the mRNAs whereas AceView does so in virtually all. We describe the driving principles of AceView, and how, by performing hand-supervised automatic annotation, we solve the combinatorial splicing problem and summarize all of GenBank, dbEST and RefSeq into a genome-wide non-redundant but comprehensive cDNA-supported transcriptome. AceView accuracy is now validated by Gencode.ConclusionRelative to a consensus mRNA catalog constructed from all evidence-based annotations, Gencode and AceView have 81% and 84% sensitivity, and 74% and 73% specificity, respectively. This close agreement validates a richer view of the human transcriptome, with three to five times more transcripts than in UCSC Known Genes (sensitivity 28%), RefSeq (sensitivity 21%) or Ensembl (sensitivity 19%).

title: The vertebrate genome annotation (Vega) database
abstract: The Vertebrate Genome Annotation (Vega) database (http://vega.sanger.ac.uk) was first made public in 2004 and has been designed to view manual annotation of human, mouse and zebrafish genomic sequences produced at the Wellcome Trust Sanger Institute. Since its initial release, the number of human annotated loci has more than doubled to close to 33 000 and now contains comprehensive annotation on 20 of the 24 human chromosomes, four whole mouse chromosomes and around 40% of the zebrafish Danio rerio genome. In addition, we offer manual annotation of a number of haplotype regions in mouse and human and regions of comparative interest in pig and dog that are unique to Vega.

title: The consensus coding sequence (CCDS) project: Identifying a common protein-coding gene set for the human and mouse genomes.
abstract: Effective use of the human and mouse genomes requires reliable identification of genes and their products. Although multiple public resources provide annotation, different methods are used that can result in similar but not identical representation of genes, transcripts, and proteins. The collaborative consensus coding sequence (CCDS) project tracks identical protein annotations on the reference mouse and human genomes with a stable identifier (CCDS ID), and ensures that they are consistently represented on the NCBI, Ensembl, and UCSC Genome Browsers. Importantly, the project coordinates on manually reviewing inconsistent protein annotations between sites, as well as annotations for which new evidence suggests a revision is needed, to progressively converge on a complete protein-coding set for the human and mouse reference genomes, while maintaining a high standard of reliability and biological accuracy. To date, the project has identified 20,159 human and 17,707 mouse consensus coding regions from 17,052 human and 16,893 mouse genes. Three evaluation methods indicate that the entries in the CCDS set are highly likely to represent real proteins, more so than annotations from contributing groups not included in CCDS. The CCDS database thus centralizes the function of identifying well-supported, identically-annotated, protein-coding regions.

title: The GENCODE exome: sequencing the complete human exome
abstract: Sequencing the coding regions, the exome, of the human genome is one of the major current strategies to identify low frequency and rare variants associated with human disease traits. So far, the most widely used commercial exome capture reagents have mainly targeted the consensus coding sequence (CCDS) database. We report the design of an extended set of targets for capturing the complete human exome, based on annotation from the GENCODE consortium. The extended set covers an additional 5594 genes and 10.3 Mb compared with the current CCDS-based sets. The additional regions include potential disease genes previously inaccessible to exome resequencing studies, such as 43 genes linked to ion channel activity and 70 genes linked to protein kinase activity. In total, the new GENCODE exome set developed here covers 47.9 Mb and performed well in sequence capture experiments. In the sample set used in this study, we identified over 5000 SNP variants more in the GENCODE exome target (24%) than in the CCDS-based exome sequencing.

title: NCBI Reference Sequences (RefSeq): current status, new features and genome annotation policy
abstract: The National Center for Biotechnology Information (NCBI) Reference Sequence (RefSeq) database is a collection of genomic, transcript and protein sequence records. These records are selected and curated from public sequence archives and represent a significant reduction in redundancy compared to the volume of data archived by the International Nucleotide Sequence Database Collaboration. The database includes over 16 000 organisms, 2.4 × 106 genomic records, 13 × 106 proteins and 2 × 106 RNA records spanning prokaryotes, eukaryotes and viruses (RefSeq release 49, September 2011). The RefSeq database is maintained by a combined approach of automated analyses, collaboration and manual curation to generate an up-to-date representation of the sequence, its features, names and cross-links to related sources of information. We report here on recent growth, the status of curating the human RefSeq data set, more extensive feature annotation and current policy for eukaryotic genome annotation via the NCBI annotation pipeline. More information about the resource is available online (see http://www.ncbi.nlm.nih.gov/RefSeq/).

title: Integrating genomic homology into gene structure prediction
abstract: TWINSCAN is a new gene-structure prediction system that directly extends the probability model of GENSCAN, allowing it to exploit homology between two related genomes. Separate probability models are used for conservation in exons, introns, splice sites, and UTRs, reflecting the differences among their patterns of evolutionary conservation. TWINSCAN is specifically designed for the analysis of high-throughput genomic sequences containing an unknown number of genes. In experiments on high-throughput mouse sequences, using homologous sequences from the human genome, TWINSCAN shows notable improvement over GENSCAN in exon sensitivity and specificity and dramatic improvement in exact gene sensitivity and specificity. This improvement can be attributed entirely to modeling the patterns of evolutionary conservation in genomic sequence.

title: Iterative gene prediction and pseudogene removal improves genome annotation.
abstract: Correct gene prediction is impaired by the presence of processed pseudogenes: nonfunctional, intronless copies of real genes found elsewhere in the genome. Gene prediction programs frequently mistake processed pseudogenes for real genes or exons, leading to biologically irrelevant gene predictions. While methods exist to identify processed pseudogenes in genomes, no attempt has been made to integrate pseudogene removal with gene prediction, or even to provide a freestanding tool that identifies such erroneous gene predictions. We have created PPFINDER (for Processed Pseudogene finder), a program that integrates several methods of processed pseudogene finding in mammalian gene annotations. We used PPFINDER to remove pseudogenes from N-SCAN gene predictions, and show that gene prediction improves substantially when gene prediction and pseudogene masking are interleaved. In addition, we used PPFINDER with gene predictions as a parent database, eliminating the need for libraries of known genes. This allows us to run the gene prediction/PPFINDER procedure on newly sequenced genomes for which few genes are known.

title: Most “Dark Matter” Transcripts Are Associated With Known Genes
abstract: Short-read RNA sequencing in mouse and human tissues shows that most transcripts are encoded within or nearby known genes and that most of the genome is not transcribed.

title: Noisy Splicing Drives mRNA Isoform Diversity in Human Cells
abstract: While the majority of multiexonic human genes show some evidence of alternative splicing, it is unclear what fraction of observed splice forms is functionally relevant. In this study, we examine the extent of alternative splicing in human cells using deep RNA sequencing and de novo identification of splice junctions. We demonstrate the existence of a large class of low abundance isoforms, encompassing approximately 150,000 previously unannotated splice junctions in our data. Newly-identified splice sites show little evidence of evolutionary conservation, suggesting that the majority are due to erroneous splice site choice. We show that sequence motifs involved in the recognition of exons are enriched in the vicinity of unconserved splice sites. We estimate that the average intron has a splicing error rate of approximately 0.7% and show that introns in highly expressed genes are spliced more accurately, likely due to their shorter length. These results implicate noisy splicing as an important property of genome evolution.

title: Transcriptome Sequencing to Detect Gene Fusions in Cancer
abstract: Recurrent gene fusions, typically associated with haematological malignancies and rare bone and soft-tissue tumours, have recently been described in common solid tumours. Here we use an integrative analysis of high-throughput long- and short-read transcriptome sequencing of cancer cells to discover novel gene fusions. As a proof of concept, we successfully used integrative transcriptome sequencing to ‘re-discover’ the BCR–ABL1 (ref. 10) gene fusion in a chronic myelogenous leukaemia cell line and the TMPRSS2–ERG gene fusion in a prostate cancer cell line and tissues. Additionally, we nominated, and experimentally validated, novel gene fusions resulting in chimaeric transcripts in cancer cell lines and tumours. Taken together, this study establishes a robust pipeline for the discovery of novel gene chimaeras using high-throughput sequencing, opening up an important class of cancer-related mutations for comprehensive characterization.

title: Landscape of transcription in human cells
abstract: Eukaryotic cells make many types of primary and processed RNAs that are found either in specific subcellular compartments or throughout the cells. A complete catalogue of these RNAs is not yet available and their characteristic subcellular localizations are also poorly understood. Because RNA represents the direct output of the genetic information encoded by genomes and a significant proportion of a cell’s regulatory capabilities are focused on its synthesis, processing, transport, modification and translation, the generation of such a catalogue is crucial for understanding genome function. Here we report evidence that three-quarters of the human genome is capable of being transcribed, as well as observations about the range and levels of expression, localization, processing fates, regulatory regions and modifications of almost all currently annotated and thousands of previously unannotated RNAs. These observations, taken together, prompt a redefinition of the concept of a gene.

title: Alternatively Spliced Human Genes by Exon Skipping - a database
abstract: Alternative splicing of mRNA allows many gene products with different functions to be produced from a single coding sequence. Exon skipping is the most commonly known alternative splicing mechanism. A comprehensive database of alternative splicing by exon skipping is made available for the human genome data. 1,229 human genes are identified to exhibit alternative splicing by exon skipping.

title: Alternative Isoform Regulation in Human Tissue Transcriptomes
abstract: Through alternative processing of pre-messenger RNAs, individual mammalian genes often produce multiple mRNA and protein isoforms that may have related, distinct or even opposing functions. Here we report an in-depth analysis of 15 diverse human tissue and cell line transcriptomes on the basis of deep sequencing of complementary DNA fragments, yielding a digital inventory of gene and mRNA isoform expression. Analyses in which sequence reads are mapped to exon–exon junctions indicated that 92–94% of human genes undergo alternative splicing, ∼86% with a minor isoform frequency of 15% or more. Differences in isoform-specific read densities indicated that most alternative splicing and alternative cleavage and polyadenylation events vary between tissues, whereas variation between individuals was approximately twofold to threefold less common. Extreme or ‘switch-like’ regulation of splicing between tissues was associated with increased sequence conservation in regulatory regions and with generation of full-length open reading frames. Patterns of alternative splicing and alternative cleavage and polyadenylation were strongly correlated across tissues, suggesting coordinated regulation of these processes, and sequence conservation of a subset of known regulatory motifs in both alternative introns and 3′ untranslated regions suggested common involvement of specific factors in tissue-level regulation of both splicing and polyadenylation.

title: Expansion of the eukaryotic proteome by alternative splicing
abstract: The collection of components required to carry out the intricate processes involved in generating and maintaining a living, breathing and, sometimes, thinking organism is staggeringly complex. Where do all of the parts come from? Early estimates stated that about 100,000 genes would be required to make up a mammal; however, the actual number is less than one-quarter of that, barely four times the number of genes in budding yeast. It is now clear that the 'missing' information is in large part provided by alternative splicing, the process by which multiple different functional messenger RNAs, and therefore proteins, can be synthesized from a single gene.

title: Targeted RNA sequencing reveals the deep complexity of the human transcriptome
abstract: Transcriptomic analyses have revealed an unexpected complexity to the human transcriptome, whose breadth and depth exceeds current RNA sequencing capability. Using tiling arrays to target and sequence select portions of the transcriptome, we identify and characterize unannotated transcripts whose rare or transient expression is below the detection limits of conventional sequencing approaches. We use the unprecedented depth of coverage afforded by this technique to reach the deepest limits of the human transcriptome, exposing widespread, regulated and remarkably complex noncoding transcription in intergenic regions, as well as unannotated exons and splicing patterns in even intensively studied protein-coding loci such as p53 and HOX. The data also show that intermittent sequenced reads observed in conventional RNA sequencing data sets, previously dismissed as noise, are in fact indicative of unassembled rare transcripts. Collectively, these results reveal the range, depth and complexity of a human transcriptome that is far from fully characterized.

title: Genome-wide detection of tissue-specific alternative splicing in the human transcriptome.
abstract: We have developed an automated method for discovering tissue-specific regulation of alternative splicing through a genome-wide analysis of expressed sequence tags (ESTs). Using this approach, we have identified 667 tissue-specific alternative splice forms of human genes. We validated our muscle-specific and brain-specific splice forms for known genes. A high fraction (8/10) were reported to have a matching tissue specificity by independent studies in the published literature. The number of tissue-specific alternative splice forms is highest in brain, while eye-retina, muscle, skin, testis and lymph have the greatest enrichment of tissue-specific splicing. Overall, 10-30% of human alternatively spliced genes in our data show evidence of tissue-specific splice forms. Seventy-eight percent of our tissue-specific alternative splices appear to be novel discoveries. We present bioinformatics analysis of several tissue-specific splice forms, including automated protein isoform sequence and domain prediction, showing how our data can provide valuable insights into gene function in different tissues. For example, we have discovered a novel kidney-specific alternative splice form of the WNK1 gene, which appears to specifically disrupt its N-terminal kinase domain and may play a role in PHAII hypertension. Our database greatly expands knowledge of tissue-specific alternative splicing and provides a comprehensive dataset for investigating its functional roles and regulation in different human tissues.

title: Variation in alternative splicing across human tissues
abstract: BackgroundAlternative pre-mRNA splicing (AS) is widely used by higher eukaryotes to generate different protein isoforms in specific cell or tissue types. To compare AS events across human tissues, we analyzed the splicing patterns of genomically aligned expressed sequence tags (ESTs) derived from libraries of cDNAs from different tissues.ResultsControlling for differences in EST coverage among tissues, we found that the brain and testis had the highest levels of exon skipping. The most pronounced differences between tissues were seen for the frequencies of alternative 3' splice site and alternative 5' splice site usage, which were about 50 to 100% higher in the liver than in any other human tissue studied. Quantifying differences in splice junction usage, the brain, pancreas, liver and the peripheral nervous system had the most distinctive patterns of AS. Analysis of available microarray expression data showed that the liver had the most divergent pattern of expression of serine-arginine protein and heterogeneous ribonucleoprotein genes compared to the other human tissues studied, possibly contributing to the unusually high frequency of alternative splice site usage seen in liver. Sequence motifs enriched in alternative exons in genes expressed in the brain, testis and liver suggest specific splicing factors that may be important in AS regulation in these tissues.ConclusionsThis study distinguishes the human brain, testis and liver as having unusually high levels of AS, highlights differences in the types of AS occurring commonly in different tissues, and identifies candidate cis-regulatory elements and trans-acting factors likely to have important roles in tissue-specific AS in human cells.

title: Estimation of alternative splicing variability in human populations.
abstract: DNA arrays have been widely used to perform transcriptome-wide analysis of gene expression, and many methods have been developed to measure gene expression variability and to compare gene expression between conditions. Because RNA-seq is also becoming increasingly popular for transcriptome characterization, the possibility exists for further quantification of individual alternative transcript isoforms, and therefore for estimating the relative ratios of alternative splice forms within a given gene. Changes in splicing ratios, even without changes in overall gene expression, may have important phenotypic effects. Here we have developed statistical methodology to measure variability in splicing ratios within conditions, to compare it between conditions, and to identify genes with condition-specific splicing ratios. Furthermore, we have developed methodology to deconvolute the relative contribution of variability in gene expression versus variability in splicing ratios to the overall variability of transcript abundances. As a proof of concept, we have applied this methodology to estimates of transcript abundances obtained from RNA-seq experiments in lymphoblastoid cells from Caucasian and Yoruban individuals. We have found that protein-coding genes exhibit low splicing variability within populations, with many genes exhibiting constant ratios across individuals. When comparing these two populations, we have found that up to 10% of the studied protein-coding genes exhibit population-specific splicing ratios. We estimate that ~60% of the total variability observed in the abundance of transcript isoforms can be explained by variability in transcription. A large fraction of the remaining variability can likely result from variability in splicing. Finally, we also detected that variability in splicing is uncommon without variability in transcription.

title: RNA integrity in post-mortem samples: influencing parameters and implications on RT-qPCR assays
abstract: Messenger RNA (mRNA) profiling in post-mortem human tissue might reveal information about gene expression at the time point of death or close to it. When working with post-mortem human tissue, one is confronted with a natural RNA degradation caused by several parameters which are not yet fully understood. The aims of the present study were to analyse the influence of impaired RNA integrity on the reliability of quantitative gene expression data and to identify ante- and post-mortem parameters that might lead to reduced RNA integrities in post-mortem human brain, cardiac muscle and skeletal muscle tissues. Furthermore, this study determined the impact of several parameters like type of tissue, age at death, gender and body mass index (BMI), as well as duration of agony, cause of death and post-mortem interval on the RNA integrity. The influence of RNA integrity on the reliability of quantitative gene expression data was analysed by generating degradation profiles for three gene transcripts. Based on the deduced cycle of quantification data, this study shows that reverse transcription quantitative polymerase chain reaction (RT-qPCR) performance is affected by impaired RNA integrity. Depending on the transcript and tissue type, a shift in cycle threshold values of up to two cycles was observed. Determining RNA integrity number of 136 post-mortem samples revealed significantly different RNA qualities among the three tissue types with brain revealing significantly lower integrities compared to skeletal and cardiac muscle. The body mass index was found to influence RNA integrity in skeletal muscle tissue (M. iliopsoas). Samples originating from deceased with a BMI > 25 were of significantly lower integrity compared to samples from normal weight donors. Correct data normalisation was found to partly diminish the effects caused by impaired RNA quality. Nevertheless, it can be concluded that in post-mortem tissue with low RNA integrity numbers, the detection of large differences in gene expression activities might still be possible, whereas small expression differences are prone to misinterpretation due to degradation. Thus, when working with post-mortem samples, we recommend generating degradation profiles for all transcripts of interest in order to reveal detection limits of RT-qPCR assays.

title: Race-and gender-related variation in NKp 46 expression associated with differential anti-HCV immunity
abstract: Major racial and gender differences have been documented in the natural history and treatment responses of chronic hepatitis C virus (HCV) infection, however, potential mechanisms have remained enigmatic. We hypothesized that racialand gender-related differences in NK cell populations may explain altered natural history and treatment responses. Our study cohort consisted of 29 African-American (AA, 55% male) and 29 Caucasian-American (CA, 48% male) normal uninfected control subjects. Multi-parameter flowcytometric analysis was used to characterize levels, phenotype with respect to 14 NK receptors, and lymphokine-activated killing (LAK) function. Gene expression was assessed by real-time RTPCR after 6 hour in vitro stimulation with TLR ligands. The ability to control HCV infection was assessed in the Huh-7.5/JFH-1 co-culture system. NK expression of natural cytotoxicity receptor NKp46 was strongly associated with CA race and female gender and correlated positively with LAK activity (p=0.0054). NKp46high NKs were more efficient at controlling HCV than their NKp46low counterparts (p<0.001). Similarly, ligation of NKp46 on isolated NK cells resulted in a significant reduction in the HCV copy number detected in Huh-7.5/JFH-1 co-culture (MOI 0.01) at an effector:target ratio of 5:1 (p<0.005). After TLR stimulation, genes involved in cytotoxicity but not cytokine genes were significantly upregulated in NKp46high NKs. Cytokine stimulation (IL-12 and IL-15) demonstrated that NKp46high NK cells have significantly higher IFN-γ production than NKp46low cells. TLR stimulation significantly induced degranulation as well as TRAIL, Fas and TNF-α protein expression in NKp46high NKs. NKp46 ligand was induced on HCV-infected hepatocytes . NKp46 expression may contribute to differential HCV responses. NKp46 expression correlates with anti-HCV activity in vitro and thus may prove to be a useful therapeutic target.

title: Genome-wide computational identification and manual annotation of human long noncoding RNA genes.
abstract: Experimental evidence suggests that half or more of the mammalian transcriptome consists of noncoding RNA. Noncoding RNAs are divided into short noncoding RNAs (including microRNAs) and long noncoding RNAs (lncRNAs). We defined complementary DNAs (cDNAs) lacking any positive-strand open reading frames (ORFs) longer than 30 amino acids, as well as cDNAs lacking any evidence of interspecies conservation of their longer-than-30-amino acid ORFs, as noncoding. We have identified 5446 lncRNA genes in the human genome from approximately 24,000 full-length cDNAs, using our new ORF-prediction pipeline. We combined them nonredundantly with lncRNAs from four published sources to derive 6736 lncRNA genes. In an effort to distinguish standalone and antisense lncRNA genes from database artifacts, we stratified our catalog of lncRNAs according to the distance between each lncRNA gene candidate and its nearest known protein-coding gene. We concurrently examined the protein-coding capacity of known genes overlapping with lncRNAs. Remarkably, 62% of known genes with "hypothetical protein" names actually lacked protein-coding capacity. This study has greatly expanded the known human lncRNA catalog, increased its accuracy through manual annotation of cDNA-to-genome alignments, and revealed that a large set of hypothetical-protein genes in GenBank lacks protein-coding capacity. In addition, we have developed, independently of existing NCBI tools, command-line programs with high-throughput ORF-finding and BLASTP-parsing functionality, suitable for future automated assessments of protein-coding capacity of novel transcripts.

title: Mapping and quantifying mammalian transcriptomes by RNA-Seq
abstract: We have mapped and quantified mouse transcriptomes by deeply sequencing them and recording how frequently each gene is represented in the sequence sample (RNA-Seq). This provides a digital measure of the presence and prevalence of transcripts from known and previously unknown genes. We report reference measurements composed of 41–52 million mapped 25-base-pair reads for poly(A)-selected RNA from adult mouse brain, liver and skeletal muscle tissues. We used RNA standards to quantify transcript prevalence and to test the linear range of transcript detection, which spanned five orders of magnitude. Although >90% of uniquely mapped reads fell within known exons, the remaining data suggest new and revised gene models, including changed or additional promoters, exons and 3′ untranscribed regions, as well as new candidate microRNA precursors. RNA splice events, which are not readily measured by standard gene expression microarray or serial analysis of gene expression methods, were detected directly by mapping splice-crossing sequence reads. We observed 1.45 × 105 distinct splices, and alternative splices were prominent, with 3,500 different genes expressing one or more alternate internal splices.

title: Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing
abstract: We carried out the first analysis of alternative splicing complexity in human tissues using mRNA-Seq data. New splice junctions were detected in ∼20% of multiexon genes, many of which are tissue specific. By combining mRNA-Seq and EST-cDNA sequence data, we estimate that transcripts from ∼95% of multiexon genes undergo alternative splicing and that there are ∼100,000 intermediate- to high-abundance alternative splicing events in major human tissues. From a comparison with quantitative alternative splicing microarray profiling data, we also show that mRNA-Seq data provide reliable measurements for exon inclusion levels.

title: Current-generation high-throughput sequencing: deepening insights into mammalian transcriptomes.
abstract: Recent papers have described the first application of high-throughput sequencing (HTS) technologies to the characterization of transcriptomes. These studies emphasize the tremendous power of this new technology, in terms of both profiling coverage and quantitative accuracy. Initial discoveries include the detection of substantial new transcript complexity, the elucidation of binding maps and regulatory properties of RNA-binding proteins, and new insights into the links between different steps in pre-mRNA processing. We review these findings, focusing on results from profiling mammalian transcriptomes. The strengths and limitations of HTS relative to microarray profiling are discussed. We also consider how future advances in HTS technology are likely to transform our understanding of integrated cellular networks operating at the RNA level.

title: Revealing the missing expressed genes beyond the human reference genome by RNA-Seq
abstract: BackgroundThe complete and accurate human reference genome is important for functional genomics researches. Therefore, the incomplete reference genome and individual specific sequences have significant effects on various studies.Resultswe used two RNA-Seq datasets from human brain tissues and 10 mixed cell lines to investigate the completeness of human reference genome. First, we demonstrated that in previously identified ~5 Mb Asian and ~5 Mb African novel sequences that are absent from the human reference genome of NCBI build 36, ~211 kb and ~201 kb of them could be transcribed, respectively. Our results suggest that many of those transcribed regions are not specific to Asian and African, but also present in Caucasian. Then, we found that the expressions of 104 RefSeq genes that are unalignable to NCBI build 37 in brain and cell lines are higher than 0.1 RPKM. 55 of them are conserved across human, chimpanzee and macaque, suggesting that there are still a significant number of functional human genes absent from the human reference genome. Moreover, we identified hundreds of novel transcript contigs that cannot be aligned to NCBI build 37, RefSeq genes and EST sequences. Some of those novel transcript contigs are also conserved among human, chimpanzee and macaque. By positioning those contigs onto the human genome, we identified several large deletions in the reference genome. Several conserved novel transcript contigs were further validated by RT-PCR.ConclusionOur findings demonstrate that a significant number of genes are still absent from the incomplete human reference genome, highlighting the importance of further refining the human reference genome and curating those missing genes. Our study also shows the importance of de novo transcriptome assembly. The comparative approach between reference genome and other related human genomes based on the transcriptome provides an alternative way to refine the human reference genome.

title: Exploring the Transcriptome of Ciliated Cells Using In Silico Dissection of Human Tissues
abstract: Cilia are cell organelles that play important roles in cell motility, sensory and developmental functions and are involved in a range of human diseases, known as ciliopathies. Here, we search for novel human genes related to cilia using a strategy that exploits the previously reported tendency of cell type-specific genes to be coexpressed in the transcriptome of complex tissues. Gene coexpression networks were constructed using the noise-resistant WGCNA algorithm in 12 publicly available microarray datasets from human tissues rich in motile cilia: airways, fallopian tubes and brain. A cilia-related coexpression module was detected in 10 out of the 12 datasets. A consensus analysis of this module's gene composition recapitulated 297 known and predicted 74 novel cilia-related genes. 82% of the novel candidates were supported by tissue-specificity expression data from GEO and/or proteomic data from the Human Protein Atlas. The novel findings included a set of genes (DCDC2, DYX1C1, KIAA0319) related to a neurological disease dyslexia suggesting their potential involvement in ciliary functions. Furthermore, we searched for differences in gene composition of the ciliary module between the tissues. A multidrug-and-toxin extrusion transporter MATE2 (SLC47A2) was found as a brain-specific central gene in the ciliary module. We confirm the localization of MATE2 in cilia by immunofluorescence staining using MDCK cells as a model. While MATE2 has previously gained attention as a pharmacologically relevant transporter, its potential relation to cilia is suggested for the first time. Taken together, our large-scale analysis of gene coexpression networks identifies novel genes related to human cell cilia.

title: Next-generation sequencing facilitates quantitative analysis of wild-type and Nrl−/− retinal transcriptomes
abstract: Next-generation sequencing (NGS) has revolutionized systems-based analysis of cellular pathways. The goals of this study are to compare NGS-derived retinal transcriptome profiling (RNA-seq) to microarray and quantitative reverse transcription polymerase chain reaction (qRT–PCR) methods and to evaluate protocols for optimal high-throughput data analysis.  Retinal mRNA profiles of 21-day-old wild-type (WT) and neural retina leucine zipper knockout (Nrl−/−) mice were generated by deep sequencing, in triplicate, using Illumina GAIIx. The sequence reads that passed quality filters were analyzed at the transcript isoform level with two methods: Burrows–Wheeler Aligner (BWA) followed by ANOVA (ANOVA) and TopHat followed by Cufflinks. qRT–PCR validation was performed using TaqMan and SYBR Green assays.  Using an optimized data analysis workflow, we mapped about 30 million sequence reads per sample to the mouse genome (build mm9) and identified 16,014 transcripts in the retinas of WT and Nrl−/− mice with BWA workflow and 34,115 transcripts with TopHat workflow. RNA-seq data confirmed stable expression of 25 known housekeeping genes, and 12 of these were validated with qRT–PCR. RNA-seq data had a linear relationship with qRT–PCR for more than four orders of magnitude and a goodness of fit (R2) of 0.8798. Approximately 10% of the transcripts showed differential expression between the WT and Nrl−/− retina, with a fold change ≥1.5 and p value <0.05. Altered expression of 25 genes was confirmed with qRT–PCR, demonstrating the high degree of sensitivity of the RNA-seq method. Hierarchical clustering of differentially expressed genes uncovered several as yet uncharacterized genes that may contribute to retinal function. Data analysis with BWA and TopHat workflows revealed a significant overlap yet provided complementary insights in transcriptome profiling.  Our study represents the first detailed analysis of retinal transcriptomes, with biologic replicates, generated by RNA-seq technology. The optimized data analysis workflows reported here should provide a framework for comparative investigations of expression profiles. Our results show that NGS offers a comprehensive and more accurate quantitative and qualitative evaluation of mRNA content within a cell or tissue. We conclude that RNA-seq based transcriptome characterization would expedite genetic network analyses and permit the dissection of complex biologic functions.

title: Genome-wide transcriptome analysis in murine neural retina using high-throughput RNA sequencing.
abstract: Genome-wide characterization of the retinal transcriptome is central to understanding development, physiology and disorders of the visual system. Massively parallel, short-read sequencing of mRNA libraries was used to generate an extensive map of the transcriptome of the adult, murine neural retina. RNA-seq data strongly corroborates prior transcriptome studies by microarray and SAGE. However, several novel features of the retinal transcriptome were discovered. For example, retinal disease genes were discovered to be among the most highly expressed in the transcriptome. We also demonstrate other interesting features of the retinal transcriptome, for example, that the retina appears to employ a very specific and restricted set of synaptic vesicle genes, and also that there is persistence of expression of a majority of "neurodevelopmental" genes into adulthood. Retina transcriptome studies utilizing novel sequencing methods have been highly informative and these data may also serve as a resource for the community of researchers.

title: Understanding alternative splicing: towards a cellular code
abstract: In violation of the 'one gene, one polypeptide' rule, alternative splicing allows individual genes to produce multiple protein isoforms — thereby playing a central part in generating complex proteomes. Alternative splicing also has a largely hidden function in quantitative gene control, by targeting RNAs for nonsense-mediated decay. Traditional gene-by-gene investigations of alternative splicing mechanisms are now being complemented by global approaches. These promise to reveal details of the nature and operation of cellular codes that are constituted by combinations of regulatory elements in pre-mRNA substrates and by cellular complements of splicing regulators, which together determine regulated splicing pathways.

title: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation.
abstract: High-throughput mRNA sequencing (RNA-Seq) promises simultaneous transcript discovery and abundance estimation. However, this would require algorithms that are not restricted by prior gene annotations and that account for alternative transcription and splicing. Here we introduce such algorithms in an open-source software program called Cufflinks. To test Cufflinks, we sequenced and analyzed >430 million paired 75-bp RNA-Seq reads from a mouse myoblast cell line over a differentiation time series. We detected 13,692 known transcripts and 3,724 previously unannotated ones, 62% of which are supported by independent expression data or by homologous genes in other species. Over the time series, 330 genes showed complete switches in the dominant transcription start site (TSS) or splice isoform, and we observed more subtle shifts in 1,304 other genes. These results suggest that Cufflinks can illuminate the substantial regulatory flexibility and complexity in even this well-studied model of muscle development and that it can improve transcriptome-based genome annotation.

title: Most “Dark Matter” Transcripts Are Associated With Known Genes
abstract: Short-read RNA sequencing in mouse and human tissues shows that most transcripts are encoded within or nearby known genes and that most of the genome is not transcribed.

title: The nonsense-mediated decay RNA surveillance pathway.
abstract: Nonsense-mediated mRNA decay (NMD) is a quality-control mechanism that selectively degrades mRNAs harboring premature termination (nonsense) codons. If translated, these mRNAs can produce truncated proteins with dominant-negative or deleterious gain-of-function activities. In this review, we describe the molecular mechanism of NMD. We first cover conserved factors known to be involved in NMD in all eukaryotes. We then describe a unique protein complex that is deposited on mammalian mRNAs during splicing, which defines a stop codon as premature. Interaction between this exon-junction complex (EJC) and NMD factors assembled at the upstream stop codon triggers a series of steps that ultimately lead to mRNA decay. We discuss whether these proofreading events preferentially occur during a "pioneer" round of translation in higher and lower eukaryotes, their cellular location, and whether they can use alternative EJC factors or act independent of the EJC.

title: Alternative Splicing and Nonsense-Mediated RNA Decay Contribute to the Regulation of SHOX Expression
abstract: The human SHOX gene is composed of seven exons and encodes a paired-related homeodomain transcription factor. SHOX mutations or deletions have been associated with different short stature syndromes implying a role in growth and bone formation. During development, SHOX is expressed in a highly specific spatiotemporal expression pattern, the underlying regulatory mechanisms of which remain largely unknown. We have analysed SHOX expression in diverse embryonic, fetal and adult human tissues and detected expression in many tissues that were not known to express SHOX before, e.g. distinct brain regions. By using RT-PCR and comparing the results with RNA-Seq data, we have identified four novel exons (exon 2a, 7-1, 7-2 and 7-3) contributing to different SHOX isoforms, and also established an expression profile for the emerging new SHOX isoforms. Interestingly, we found the exon 7 variants to be exclusively expressed in fetal neural tissues, which could argue for a specific role of these variants during brain development. A bioinformatical analysis of the three novel 3′UTR exons yielded insights into the putative role of the different 3′UTRs as targets for miRNA binding. Functional analysis revealed that inclusion of exon 2a leads to nonsense-mediated RNA decay altering SHOX expression in a tissue and time specific manner. In conclusion, SHOX expression is regulated by different mechanisms and alternative splicing coupled with nonsense-mediated RNA decay constitutes a further component that can be used to fine tune the SHOX expression level.

title: Specific regulation of low-abundance transcript variants encoding human Acyl-CoA binding protein (ACBP) isoforms
abstract: Despite intensive efforts on annotation of eukaryotic transcriptoms, little is known about the regulation of low‐abundance transcripts. To address this question, we analysed the regulation of novel low‐abundance transcript variants of human acyl‐CoA binding protein (ACBP), an important multifunctional housekeeping protein, which we have identified by screening of human expressed sequence tags in combination with ab initio gene prediction. By using RT‐, real‐time RT‐ and rapid amplification of cDNA ends‐PCR in five human tissues, we find these transcripts, which are generated by a consequent use of alternative promoters and alternate first or first two exons, to be authentic ones. They show a tissue‐specific distribution and intrinsic responsiveness to glucose and insulin. Promoter analyses of the corresponding transcripts revealed a differential regulation mediated by sterol regulatory element‐binding protein‐2, hepatocyte nuclear factor‐4α and nuclear factor κB (NF‐κB), central transcription factors of fat and glucose metabolism and inflammation. Subcellular localization studies of deduced isoforms in liver HepG2 cells showed that they are distributed in different compartments. By demonstrating that ACBP is a target of NF‐κB, our findings link fatty acid metabolism with inflammation. Furthermore, our findings show that low‐abundance transcripts are regulated in a similar mode than their high‐abundance counterparts.

title: Investigating the Integration of Alternative Splicing and Transcriptional Regulation in Mammalian Gene Expression
abstract: Alternative splicing functions to generate proteomic diversity and to regulate gene expression in higher eukaryotes. Genome-wide analyses suggest that alternative splicing and transcription typically regulate different gene sets to achieve cell- and tissue-type specificity. However, within individual cell-types, most alternative splicing events occur co-transcriptionally and are impacted by the transcriptional machinery. Despite many focused studies on co-transcriptional regulation of alternative splicing, its mechanisms and functions in regulation of gene expression are still poorly understood.To investigate relationships between transcription and alternative splicing, I performed microarray profiling of alternative splicing and transcript levels during activation of a T cell line. This experiment revealed that different sets of genes and associated functional categories are regulated by alternative splicing and transcription during T cell activation. I next employed inhibitors of RNA polymerase II (Pol II) elongation and microarray profiling to identify genes with coupled changes in splicing and transcript levels when transcription is impeded in activated T cell. Genes that were affected at both levels were significantly enriched in RNA binding and processing functions, and generally displayed increased alternative exon inclusion and decreased transcript levels when transcription elongation was disrupted. Similar effects were observed when transcription was driven by mutant polymerases with reduced elongation activity, and when cells were subjected to stress treatments. Many of the elongation inhibition-sensitive exons from the affected genes introduce premature termination codons into the mRNA, resulting in spliced mRNAs that are substrates of the nonsense-mediated decay pathway and further reduction in mRNA levels. ChIP-Seq experiment demonstrated that Pol II occupancy specifically increased in introns flanking the affected exons. These results provide evidence that a physiological function of transcription elongation-coupled alternative splicing regulation is to regulate the levels of RNA processing factors under conditions that reduce elongation activity, including cell stress. In summary, my thesis research has provided new insights into the integration of transcription and splicing control. While these two regulatory levels can control different gene sets during the activation of T cells, within a given cell type, they are closely coupled to control specific alternative splicing events that appear to coordinate mRNA and RNA processing factors levels.

title: Alternative Splicing of RNA Triplets Is Often Regulated and Accelerates Proteome Evolution
abstract: Inclusion or exclusion of single codons at the splice acceptor site of mammalian genes is regulated in a tissue-specific manner, is strongly conserved, and is associated with local accelerated protein evolution.

title: Exome sequencing identifies MLL2 mutations as a cause of Kabuki syndrome
abstract: We demonstrate the successful application of exome sequencing to discover a gene for an autosomal dominant disorder, Kabuki syndrome (OMIM%147920). We subjected the exomes of ten unrelated probands to massively parallel sequencing. After filtering against existing SNP databases, there was no compelling candidate gene containing previously unknown variants in all affected individuals. Less stringent filtering criteria allowed for the presence of modest genetic heterogeneity or missing data but also identified multiple candidate genes. However, genotypic and phenotypic stratification highlighted MLL2, which encodes a Trithorax-group histone methyltransferase: seven probands had newly identified nonsense or frameshift mutations in this gene. Follow-up Sanger sequencing detected MLL2 mutations in two of the three remaining individuals with Kabuki syndrome (cases) and in 26 of 43 additional cases. In families where parental DNA was available, the mutation was confirmed to be de novo (n = 12) or transmitted (n = 2) in concordance with phenotype. Our results strongly suggest that mutations in MLL2 are a major cause of Kabuki syndrome.

title: Spectrum of MLL2 (ALR) mutations in 110 cases of Kabuki syndrome
abstract: Kabuki syndrome is a rare, multiple malformation disorder characterized by a distinctive facial appearance, cardiac anomalies, skeletal abnormalities, and mild to moderate intellectual disability. Simplex cases make up the vast majority of the reported cases with Kabuki syndrome, but parent‐to‐child transmission in more than a half‐dozen instances indicates that it is an autosomal dominant disorder. We recently reported that Kabuki syndrome is caused by mutations in MLL2, a gene that encodes a Trithorax‐group histone methyltransferase, a protein important in the epigenetic control of active chromatin states. Here, we report on the screening of 110 families with Kabuki syndrome. MLL2 mutations were found in 81/110 (74%) of families. In simplex cases for which DNA was available from both parents, 25 mutations were confirmed to be de novo, while a transmitted MLL2 mutation was found in two of three familial cases. The majority of variants found to cause Kabuki syndrome were novel nonsense or frameshift mutations that are predicted to result in haploinsufficiency. The clinical characteristics of MLL2 mutation‐positive cases did not differ significantly from MLL2 mutation‐negative cases with the exception that renal anomalies were more common in MLL2 mutation‐positive cases. These results are important for understanding the phenotypic consequences of MLL2 mutations for individuals and their families as well as for providing a basis for the identification of additional genes for Kabuki syndrome.

title: Transcriptome sequencing across a prostate cancer cohort identifies PCAT-1, an unannotated lincRNA implicated in disease progression
abstract: Noncoding RNAs (ncRNAs) are emerging as key molecules in human cancer, with the potential to serve as novel markers of disease and to reveal uncharacterized aspects of tumor biology. Here we discover 121 unannotated prostate cancer–associated ncRNA transcripts (PCATs) by ab initio assembly of high-throughput sequencing of polyA+ RNA (RNA-Seq) from a cohort of 102 prostate tissues and cells lines. We characterized one ncRNA, PCAT-1, as a prostate-specific regulator of cell proliferation and show that it is a target of the Polycomb Repressive Complex 2 (PRC2). We further found that patterns of PCAT-1 and PRC2 expression stratified patient tissues into molecular subtypes distinguished by expression signatures of PCAT-1–repressed target genes. Taken together, our findings suggest that PCAT-1 is a transcriptional repressor implicated in a subset of prostate cancer patients. These findings establish the utility of RNA-Seq to identify disease-associated ncRNAs that may improve the stratification of cancer subtypes.

title: Transcriptome-wide sequencing reveals numerous APOBEC1 mRNA editing targets in transcript 3′ UTRs
abstract: Apolipoprotein B–editing enzyme, catalytic polypeptide-1 (APOBEC1) is a cytidine deaminase initially identified by its activity in converting a specific cytidine (C) to uridine (U) in apolipoprotein B (apoB) mRNA transcripts in the small intestine. Editing results in the translation of a truncated apoB isoform with distinct functions in lipid transport. To address the possibility that APOBEC1 edits additional mRNAs, we developed a transcriptome-wide comparative RNA sequencing (RNA-Seq) screen. We identified and validated 32 previously undescribed mRNA targets of APOBEC1 editing, all of which are located in AU-rich segments of transcript 3′ untranslated regions (3′ UTRs). Further analysis established several characteristic sequence features of editing targets, which were predictive for the identification of additional APOBEC1 substrates. The transcriptomics approach to RNA editing presented here dramatically expands the list of APOBEC1 mRNA editing targets and reveals a novel cellular mechanism for the modification of transcript 3′ UTRs.

title: Transcriptome profiling in neurodegenerative disease
abstract: Changes in gene expression and splicing patterns (that occur prior to the onset and during the progression of complex diseases) have become a major focus of neurodegenerative disease research. These signature patterns of gene expression provide clues about the mechanisms involved in the molecular pathogenesis of neurodegenerative disease and may facilitate the discovery of novel therapeutic drugs. With the development of array technologies and the very recent RNA-seq technique, our understanding of the pathogenesis of neurodegenerative disease is expanding exponentially. Here, we review the technologies involved in gene expression and splicing analysis and the related literature on three common neurodegenerative diseases: Alzheimer's disease, Parkinson's disease and Huntington's disease.

title: Analysis of next-generation genomic data in cancer: accomplishments and challenges.
abstract: The application of next-generation sequencing technology has produced a transformation in cancer genomics, generating large data sets that can be analyzed in different ways to answer a multitude of questions about the genomic alterations associated with the disease. Analytical approaches can discover focused mutations such as substitutions and small insertion/deletions, large structural alterations and copy number events. As our capacity to produce such data for multiple cancers of the same type is improving, so are the demands to analyze multiple tumor genomes simultaneously growing. For example, pathway-based analyses that provide the full mutational impact on cellular protein networks and correlation analyses aimed at revealing causal relationships between genomic alterations and clinical presentations are both enabled. As the repertoire of data grows to include mRNA-seq, non-coding RNA-seq and methylation for multiple genomes, our challenge will be to intelligently integrate data types and genomes to produce a coherent picture of the genetic basis of cancer.

title: Understanding the pathogenesis of Alzheimer’s disease: will RNA‐Seq realize the promise of transcriptomics?
abstract: (2011) 116, 937–946.

title: Exome RNA sequencing reveals rare and novel alternative transcripts
abstract: RNA sequencing has become an important method to perform hypothesis-free characterization of global gene expression. One of the limitations of RNA sequencing is that most sequence reads represent highly expressed transcripts, whereas low level transcripts are challenging to detect. To combine the benefits of traditional expression arrays with the advantages of RNA sequencing, we have used whole exome enrichment prior to sequencing of total RNA. We show that whole exome capture can be successfully applied to cDNA to study the transcriptional landscape in human tissues. By introducing the exome enrichment step, we are able to identify transcripts present at very low levels, which are below the level of detection in conventional RNA sequencing. Although the enrichment increases the ability to detect presence of transcripts, it also lowers the accuracy of quantification of expression levels. Our results yield a large number of novel exons and splice isoforms, suggesting that conventional RNA sequencing methods only detect a small fraction of the full transcript diversity. We propose that whole exome enrichment of RNA is a suitable strategy for genome-wide discovery of novel transcripts, alternative splice variants and fusion genes.

title: Missing lincs in the transcriptome
abstract: Are long, intervening noncoding (linc) RNAs a new class of functional transcripts?

title: LINCing chromatin remodeling to metastasis
abstract: A long intergenic noncoding RNA may promote metastatic progression by coordinating the activity of histone-modifying enzymes.

title: Linc-ing Long noncoding RNAs and enhancer function.
abstract: Enhancers are distal regulatory sequences that control gene expression in development. Ørom et al. now report in Cell that some long noncoding RNAs have functional properties of enhancers. Known enhancers are also transcribed in cells in which they are active, suggesting that noncoding RNAs are integral to enhancer action.

title: NCBI BLAST: a better web interface
abstract: Basic Local Alignment Search Tool (BLAST) is a sequence similarity search program. The public interface of BLAST, http://www.ncbi.nlm.nih.gov/blast, at the NCBI website has recently been reengineered to improve usability and performance. Key new features include simplified search forms, improved navigation, a list of recent BLAST results, saved search strategies and a documentation directory. Here, we describe the BLAST web application's new features, explain design decisions and outline plans for future improvement.