The Retinome - defining a reference transcriptome of the adult mammalian retina/retinal pigment epithelium

Functions and mechanisms of RNA tailing by metazoan terminal nucleotidyltransferases

Termini often determine the fate of RNA molecules. In recent years, 3' ends of almost all classes of RNA species have been shown to acquire nontemplated nucleotides that are added by terminal nucleotidyltransferases (TENTs). The best-described role of 3' tailing is the bulk polyadenylation of messenger RNAs in the cell nucleus that is catalyzed by canonical poly(A) polymerases (PAPs). However, many other enzymes that add adenosines, uridines, or even more complex combinations of nucleotides have recently been described. This review focuses on metazoan TENTs, which are either noncanonical PAPs or terminal uridylyltransferases with varying processivity. These enzymes regulate RNA stability and RNA functions and are crucial in early development, gamete production, and somatic tissues. TENTs regulate gene expression at the posttranscriptional level, participate in the maturation of many transcripts, and protect cells against viral invasion and the transposition of repetitive sequences. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Processing > 3' End Processing RNA Turnover and Surveillance > Regulation of RNA Stability.

Terminal nucleotidyl transferases (TENTs) in mammalian RNA metabolism

In eukaryotes, almost all RNA species are processed at their 3′ ends and most mRNAs are polyadenylated in the nucleus by canonical poly(A) polymerases. In recent years, several terminal nucleotidyl transferases (TENTs) including non-canonical poly(A) polymerases (ncPAPs) and terminal uridyl transferases (TUTases) have been discovered. In contrast to canonical polymerases, TENTs' functions are more diverse; some, especially TUTases, induce RNA decay while others, such as cytoplasmic ncPAPs, activate translationally dormant deadenylated mRNAs. The mammalian genome encodes 11 different TENTs. This review summarizes the current knowledge about the functions and mechanisms of action of these enzymes.

This article is part of the theme issue ‘5′ and 3′ modifications controlling RNA degradation’.

Characterization of lincRNA expression in the human retinal pigment epithelium and differentiated induced pluripotent stem cells

Long intervening non-coding RNAs (lincRNAs) are increasingly being implicated as important factors in many aspects of cellular development, function, and disease, but remain poorly understood. In this study, we examine the human retinal pigment epithelium (RPE) lincRNA transcriptome using RNA-Seq data generated from human fetal RPE (fRPE), RPE derived from human induced pluripotent stem cells (iPS-RPE), and undifferentiated iPS (iPS). In addition, we determine the suitability of iPS-RPE, from a transcriptome standpoint, as a model for use in future studies of lincRNA structure and function. A comparison of gene and isoform expression across the whole transcriptome shows only minimal differences between all sample types, though fRPE and iPS-RPE show higher concordance than either shows with iPS. Notably, RPE signature genes show the highest degree of fRPE to iPS-RPE concordance, indicating that iPS-RPE cells provide a suitable model for use in future studies. An analysis of lincRNAs demonstrates high concordance between fRPE and iPS-RPE, but low concordance between either RPE and iPS. While most lincRNAs are expressed at low levels (RPKM < 10), there is a high degree of concordance among replicates within each sample type, suggesting the expression is consistent, even at levels subject to high variability. Finally, we identified and annotated 180 putative novel genes in the fRPE samples, a majority of which are also expressed in the iPS-RPE. Overall, this study represents the first characterization of lincRNA expression in the human RPE, and provides a model for studying the role lincRNAs play in RPE development, function, and disease.

Targeted Capture of Complete Coding Regions across Divergent Species

Despite continued advances in sequencing technologies, there is a need for methods that can efficiently sequence large numbers of genes from diverse species. One approach to accomplish this is targeted capture (hybrid enrichment). While these methods are well established for genome resequencing projects, cross-species capture strategies are still being developed and generally focus on the capture of conserved regions, rather than complete coding regions from specific genes of interest. The resulting data is thus useful for phylogenetic studies, but the wealth of comparative data that could be used for evolutionary and functional studies is lost. Here, we design and implement a targeted capture method that enables recovery of complete coding regions across broad taxonomic scales. Capture probes were designed from multiple reference species and extensively tiled in order to facilitate cross-species capture. Using novel bioinformatics pipelines we were able to recover nearly all of the targeted genes with high completeness from species that were up to 200 myr divergent. Increased probe diversity and tiling for a subset of genes had a large positive effect on both recovery and completeness. The resulting data produced an accurate species tree, but importantly this same data can also be applied to studies of molecular evolution and function that will allow researchers to ask larger questions in broader phylogenetic contexts. Our method demonstrates the utility of cross-species approaches for the capture of full length coding sequences, and will substantially improve the ability for researchers to conduct large-scale comparative studies of molecular evolution and function.

MicroRNA changes through Müller glia dedifferentiation and early/late rod photoreceptor differentiation

Cell-type determination is a complex process driven by the combinatorial effect of extrinsic signals and the expression of transcription factors and regulatory genes. MicroRNAs (miRNAs) are non-coding RNAs that, generally, inhibit the expression of target genes and have been involved, among other processes, in cell identity acquisition. To search for candidate miRNAs putatively involved in mice rod photoreceptor and Müller glia (MG) identity, we compared miRNA expression profiles between late-stage retinal progenitor cells (RPCs), CD73-immunopositive (CD73+) rods and postnatal MG. We found a close similarity between RPCs and CD73+ miRNA expression profiles but a divergence between CD73+ and MG miRNA signatures. We validated preferentially expressed miRNAs in the CD73+ subpopulation (miR-182, 183, 124a, 9(∗), 181c and 301b(∗)) or MG (miR-143, 145, 214, 199a-5p, 199b(∗), and 29a). Taking advantage of the unique capacity of MG to dedifferentiate into progenitor-like cells that can be differentiated to a rod phenotype in response to external cues, we evaluated changes of selected miRNAs in MG-derived progenitors (MGDP) during neuronal differentiation. We found decreased levels of miR-143 and 145, but increased levels of miR-29a in MGDP. In MGDPs committed to early neuronal lineages we found increased levels of miR-124a and upregulation of miR-124a, 9(∗) and 181c during MGDP acquisition of rod phenotypes. Furthermore, we demonstrated that ectopic miR-124 expression is sufficient to enhance early neuronal commitment of MGDP. Our data reveal a dynamic regulation of miRNAs in MGDP through early and late neuronal commitment and miRNAs that could be potential targets to exploit the silent neuronal differentiation capacity of MG in mammals.

Vision from next generation sequencing: multi-dimensional genome-wide analysis for producing gene regulatory networks underlying retinal development, aging and disease

Genomics and genetics have invaded all aspects of biology and medicine, opening uncharted territory for scientific exploration. The definition of "gene" itself has become ambiguous, and the central dogma is continuously being revised and expanded. Computational biology and computational medicine are no longer intellectual domains of the chosen few. Next generation sequencing (NGS) technology, together with novel methods of pattern recognition and network analyses, has revolutionized the way we think about fundamental biological mechanisms and cellular pathways. In this review, we discuss NGS-based genome-wide approaches that can provide deeper insights into retinal development, aging and disease pathogenesis. We first focus on gene regulatory networks (GRNs) that govern the differentiation of retinal photoreceptors and modulate adaptive response during aging. Then, we discuss NGS technology in the context of retinal disease and develop a vision for therapies based on network biology. We should emphasize that basic strategies for network construction and analyses can be transported to any tissue or cell type. We believe that specific and uniform guidelines are required for generation of genome, transcriptome and epigenome data to facilitate comparative analysis and integration of multi-dimensional data sets, and for constructing networks underlying complex biological processes. As cellular homeostasis and organismal survival are dependent on gene-gene and gene-environment interactions, we believe that network-based biology will provide the foundation for deciphering disease mechanisms and discovering novel drug targets for retinal neurodegenerative diseases.

A Candidate Gene Association Study Identifies DAPL1 as a Female-Specific Susceptibility Locus for Age-Related Macular Degeneration (AMD)

Age-related macular degeneration (AMD) is the leading cause of blindness among white caucasians over the age of 50 years with a prevalence rate expected to increase markedly with an anticipated increase in the life span of the world population. To further expand our knowledge of the genetic architecture of the disease, we pursued a candidate gene approach assessing 25 genes and a total of 109 variants. Of these, synonymous single nucleotide polymorphism (SNP) rs17810398 located in death-associated protein-like 1 (DAPL1) was found to be associated with AMD in a joint analysis of 3,229 cases and 2,835 controls from five studies [combined P_ADJ = 1.15 × 10⁻⁶, OR 1.332 (1.187–1.496)]. This association was characterized by a highly significant sex difference (P_diff = 0.0032) in that it was clearly confined to females with genome-wide significance [P_ADJ = 2.62 × 10⁻⁸, OR 1.541 (1.324–1.796); males: P_ADJ = 0.382, OR 1.084 (0.905–1.298)]. By targeted resequencing of risk and non-risk associated haplotypes in the DAPL1 locus, we identified additional potentially functional risk variants, namely a common 897-bp deletion and a SNP predicted to affect a putative binding site of an exonic splicing enhancer. We show that the risk haplotype correlates with a reduced retinal transcript level of two, less frequent, non-canonical DAPL1 isoforms. DAPL1 plays a role in epithelial differentiation and may be involved in apoptotic processes thereby suggesting a possible novel pathway in AMD pathogenesis.

Transcriptome of the human retina, retinal pigmented epithelium and choroid

The retina and its adjacent supporting tissues - retinal pigmented epithelium (RPE) and choroid - are critical structures in human eyes required for normal visual perception. Abnormal changes in these layers have been implicated in diseases such as age-related macular degeneration and glaucoma. With the advent of high-throughput methods, such as serial analysis of gene expression, cDNA microarray, and RNA sequencing, there is unprecedented opportunity to facilitate our understanding of the normal retina, RPE, and choroid. This information can be used to identify dysfunction in age-related macular degeneration and glaucoma. In this review, we describe the current status in our understanding of these transcriptomes through the use of high-throughput techniques.

Generation of a foveomacular transcriptome

PURPOSE

Organizing molecular biologic data is a growing challenge since the rate of data accumulation is steadily increasing. Information relevant to a particular biologic query can be difficult to extract from the comprehensive databases currently available. We present a data collection and organization model designed to ameliorate these problems and applied it to generate an expressed sequence tag (EST)-based foveomacular transcriptome.

METHODS

Using Perl, MySQL, EST libraries, screening, and human foveomacular gene expression as a model system, we generated a foveomacular transcriptome database enriched for molecularly relevant data.

RESULTS

Using foveomacula as a gene expression model tissue, we identified and organized 6,056 genes expressed in that tissue. Of those identified genes, 3,480 had not been previously described as expressed in the foveomacula. Internal experimental controls as well as comparison of our data set to published data sets suggest we do not yet have a complete description of the foveomacula transcriptome.

CONCLUSIONS

We present an organizational method designed to amplify the utility of data pertinent to a specific research interest. Our method is generic enough to be applicable to a variety of conditions yet focused enough to allow for specialized study.

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

BACKGROUND

The retina is a complex tissue comprised of multiple cell types that is affected by a diverse set of diseases that are important causes of vision loss. Characterizing the transcripts, both annotated and novel, that are expressed in a given tissue has become vital for understanding the mechanisms underlying the pathology of disease.

RESULTS

We sequenced RNA prepared from three normal human retinas and characterized the retinal transcriptome at an unprecedented level due to the increased depth of sampling provided by the RNA-seq approach. We used a non-redundant reference transcriptome from all of the empirically-determined human reference tracks to identify annotated and novel sequences expressed in the retina. We detected 79,915 novel alternative splicing events, including 29,887 novel exons, 21,757 3' and 5' alternate splice sites, and 28,271 exon skipping events. We also identified 116 potential novel genes. These data represent a significant addition to the annotated human transcriptome. For example, the novel exons detected increase the number of identified exons by 3%. Using a high-throughput RNA capture approach to validate 14,696 of these novel transcriptome features we found that 99% of the putative novel events can be reproducibly detected. Further, 15-36% of the novel splicing events maintain an open reading frame, suggesting they produce novel protein products.

CONCLUSIONS

To our knowledge, this is the first application of RNA capture to perform large-scale validation of novel transcriptome features. In total, these analyses provide extensive detail about a previously uncharacterized level of transcript diversity in the human retina.

Comparison of an expanded ataxia interactome with patient medical records reveals a relationship between macular degeneration and ataxia

Spinocerebellar ataxias 6 and 7 (SCA6 and SCA7) are neurodegenerative disorders caused by expansion of CAG repeats encoding polyglutamine (polyQ) tracts in CACNA1A, the alpha1A subunit of the P/Q-type calcium channel, and ataxin-7 (ATXN7), a component of a chromatin-remodeling complex, respectively. We hypothesized that finding new protein partners for ATXN7 and CACNA1A would provide insight into the biology of their respective diseases and their relationship to other ataxia-causing proteins. We identified 118 protein interactions for CACNA1A and ATXN7 linking them to other ataxia-causing proteins and the ataxia network. To begin to understand the biological relevance of these protein interactions within the ataxia network, we used OMIM to identify diseases associated with the expanded ataxia network. We then used Medicare patient records to determine if any of these diseases co-occur with hereditary ataxia. We found that patients with ataxia are at 3.03-fold greater risk of these diseases than Medicare patients overall. One of the diseases comorbid with ataxia is macular degeneration (MD). The ataxia network is significantly (P= 7.37 × 10⁻⁵) enriched for proteins that interact with known MD-causing proteins, forming a MD subnetwork. We found that at least two of the proteins in the MD subnetwork have altered expression in the retina of Ataxin-7^266Q/+ mice suggesting an in vivo functional relationship with ATXN7. Together these data reveal novel protein interactions and suggest potential pathways that can contribute to the pathophysiology of ataxia, MD, and diseases comorbid with ataxia.

Characterizing light-regulated retinal microRNAs reveals rapid turnover as a common property of neuronal microRNAs

Adaptation to different levels of illumination is central to the function of the retina. Here, we demonstrate that levels of the miR-183/96/182 cluster, miR-204, and miR-211 are regulated by different light levels in the mouse retina. Concentrations of these microRNAs were downregulated during dark adaptation and upregulated in light-adapted retinas, with rapid decay and increased transcription being responsible for the respective changes. We identified the voltage-dependent glutamate transporter Slc1a1 as one of the miR-183/96/182 targets in photoreceptor cells. We found that microRNAs in retinal neurons decay much faster than microRNAs in nonneuronal cells. The high turnover is also characteristic of microRNAs in hippocampal and cortical neurons, and neurons differentiated from ES cells in vitro. Blocking activity reduced turnover of microRNAs in neuronal cells while stimulation with glutamate accelerated it. Our results demonstrate that microRNA metabolism in neurons is higher than in most other cells types and linked to neuronal activity.

A new strategy to identify and annotate human RPE-specific gene expression

Background

To identify and functionally annotate cell type-specific gene expression in the human retinal pigment epithelium (RPE), a key tissue involved in age-related macular degeneration and retinitis pigmentosa.

Methodology

RPE, photoreceptor and choroidal cells were isolated from selected freshly frozen healthy human donor eyes using laser microdissection. RNA isolation, amplification and hybridization to 44 k microarrays was carried out according to Agilent specifications. Bioinformatics was carried out using Rosetta Resolver, David and Ingenuity software.

Principal Findings

Our previous 22 k analysis of the RPE transcriptome showed that the RPE has high levels of protein synthesis, strong energy demands, is exposed to high levels of oxidative stress and a variable degree of inflammation. We currently use a complementary new strategy aimed at the identification and functional annotation of RPE-specific expressed transcripts. This strategy takes advantage of the multilayered cellular structure of the retina and overcomes a number of limitations of previous studies. In triplicate, we compared the transcriptomes of RPE, photoreceptor and choroidal cells and we deduced RPE specific expression. We identified at least 114 entries with RPE-specific gene expression. Thirty-nine of these 114 genes also show high expression in the RPE, comparison with the literature showed that 85% of these 39 were previously identified to be expressed in the RPE. In the group of 114 RPE specific genes there was an overrepresentation of genes involved in (membrane) transport, vision and ophthalmic disease. More fundamentally, we found RPE-specific involvement in the RAR-activation, retinol metabolism and GABA receptor signaling pathways.

Conclusions

In this study we provide a further specification and understanding of the RPE transcriptome by identifying and analyzing genes that are specifically expressed in the RPE.

Comparative genomic mapping of uncharacterized canine retinal ESTs to identify novel candidate genes for hereditary retinal disorders

PURPOSE

To identify the genomic location of previously uncharacterized canine retina-expressed expressed sequence tags (ESTs), and thus identify potential candidate genes for heritable retinal disorders.

METHODS

A set of over 500 retinal canine ESTs were mapped onto the canine genome using the RHDF(5000-2) radiation hybrid (RH) panel, and the resulting map positions were compared to their respective localization in the CanFam2 assembly of the canine genome sequence.

RESULTS

Unique map positions could be assigned for 99% of the mapped clones, of which only 29% showed significant homology to known RefSeq sequences. A comparison between RH map and sequence assembly indicated some areas of discrepancy. Retinal expressed genes were not concentrated in particular areas of the canine genome, and also were located on the canine Y chromosome (CFAY). Several of the EST clones were located within areas of conserved synteny to human retinal disease loci.

CONCLUSIONS

RH mapping of canine retinal ESTs provides insight into the location of potential candidate genes for hereditary retinal disorders, and, by comparison with the assembled canine genome sequence, highlights inconsistencies with the current assembly. Regions of conserved synteny between the canine and the human genomes allow this information to be extrapolated to identify potential positional candidate genes for mapped human retinal disorders. Furthermore, these ESTs can help identify novel or uncharacterized genes of significance for better understanding of retinal morphology, physiology, and pathology.

Functional annotation of the human retinal pigment epithelium transcriptome

Background

To determine level, variability and functional annotation of gene expression of the human retinal pigment epithelium (RPE), the key tissue involved in retinal diseases like age-related macular degeneration and retinitis pigmentosa. Macular RPE cells from six selected healthy human donor eyes (aged 63–78 years) were laser dissected and used for 22k microarray studies (Agilent technologies). Data were analyzed with Rosetta Resolver, the web tool DAVID and Ingenuity software.

Results

In total, we identified 19,746 array entries with significant expression in the RPE. Gene expression was analyzed according to expression levels, interindividual variability and functionality. A group of highly (n = 2,194) expressed RPE genes showed an overrepresentation of genes of the oxidative phosphorylation, ATP synthesis and ribosome pathways. In the group of moderately expressed genes (n = 8,776) genes of the phosphatidylinositol signaling system and aminosugars metabolism were overrepresented. As expected, the top 10 percent (n = 2,194) of genes with the highest interindividual differences in expression showed functional overrepresentation of the complement cascade, essential in inflammation in age-related macular degeneration, and other signaling pathways. Surprisingly, this same category also includes the genes involved in Bruch's membrane (BM) composition. Among the top 10 percent of genes with low interindividual differences, there was an overrepresentation of genes involved in local glycosaminoglycan turnover.

Conclusion

Our study expands current knowledge of the RPE transcriptome by assigning new genes, and adding data about expression level and interindividual variation. Functional annotation suggests that the RPE has high levels of protein synthesis, strong energy demands, and is exposed to high levels of oxidative stress and a variable degree of inflammation. Our data sheds new light on the molecular composition of BM, adjacent to the RPE, and is useful for candidate retinal disease gene identification or gene dose-dependent therapeutic studies.

Identification of candidate genes for human retinal degeneration loci using differentially expressed genes from mouse photoreceptor dystrophy models

PURPOSE

Retinal degeneration (RD) is a complex mechanism that appears to involve many biologic processes including oxidative stress, apoptosis, and cellular remodeling. Currently there are 51 mapped, but not identified, RD human disease loci.

METHODS

To assign possible disease genes to RD loci, we have used a comparative genomics procedure that incorporates microarray gene expression data of three independent mouse models for photoreceptor dystrophy (rd1, rd2, and constant light-damage in BALB/c mice), human ortholog data, and databases of known chromosomal locations involved in human RD. Immunohistochemistry and enzyme activity assays were used to further characterize a candidate gene product.

RESULTS

Our analysis yielded candidate genes for four mapped, but unsolved, human chromosomal locations and confirmed two previously identified monogenic disease loci for human RD, thus validating our approach. PLA2G7 (phospholipase A2, group VII; PAF-AH, Lp-PLA2), a candidate for a dominant form macular dystrophy (Benign Concentric Annular Macular Dystrophy [BCMAD]), was selected for further study. The PLA2G7 enzyme is known to mediate breakdown of oxidatively damaged phospholipids, a contributor to oxidative stress in the retina. PLA2G7 protein was enriched in mouse photoreceptor inner and outer segments. In the rd1, rd2, and BALB/c mice exposed to constant light, retinal tissue activity levels, but not plasma levels, were significantly reduced at the onset of photoreceptor cell death.

CONCLUSIONS

We have shown that this comparative genomics approach verified existing RD genes as well as identified novel RD candidate genes. The results on the characterization of the PLA2G7 protein, one of the novel RD genes, suggests that retinal tissue PLA2G7 levels may constitute an important risk factor for BCMAD. In summary, this reverse mapping approach, using accepted mouse models of human disease and known human RD loci, may prove useful in identifying possible novel disease candidates for RD and may be applicable to other human diseases.

Identification of candidate cancer genes involved in human retinoblastoma by data mining

OBJECTIVE

The objective of this study was to discover potential cancer-related genes involved in retinoblastoma (RB) tumorigenesis.

MATERIALS AND METHODS

Using a data-mining tool called cDNA Digital Gene Expression Displayer (DGED) and serial analysis of gene expression DGED from the Cancer Genome Anatomy Project (CGAP) database, eight cDNA libraries and five serial analysis of gene expression libraries from retinoblastoma (RB) solid tumors and normal retina tissues were analyzed. The deregulated genes were classified into major families using information from Gene Ontology. Several candidate cancer-related genes were analyzed by real-time reverse transcription polymerase chain reaction (RT-PCR) and immunohistochemistry (IHC) on tissue microarrays (TMA) of RB and human normal retina samples.

RESULTS

A total of 260 genes with deregulated expression emerged when examined by DGED from the CGAP database. Functional classification of these genes not only provided an interesting insight into RB tumorigenesis but also facilitated target identification for RB therapeutics. Several candidate genes were confirmed by real-time RT-PCR and IHC analysis on TMA and were found to be associated with RB genesis through text-mining in Information Hyperlinked over Proteins. The results also implicated MCM7 and WIF1 as promising therapeutic targets for RB, but further validation is needed.

Genetic analysis of FAM46A in Spanish families with autosomal recessive retinitis pigmentosa: characterisation of novel VNTRs

Retinitis pigmentosa (RP) is a group of retinal dystrophies characterised primarily by rod photoreceptor cell degeneration. Exhibiting great clinical and genetic heterogeneity, RP be inherited as an autosomal dominant (ad) and recessive (ar), X-linked (xl) and digenic disorder. RP25, a locus for arRP, was mapped to chromosome 6p12.1-q14.1 where several retinal dystrophy loci are located. A gene expressed in the retina, FAM46A, mapped within the RP25 locus, and computational data revealed its involvement in retinal signalling pathways. Therefore, we chose to perform molecular evaluation of this gene as a good candidate in arRP families linked to the RP25 interval. A comprehensive bioinformatic and retinal tissue expression characterisation of FAM46A was performed, together with mutation screening of seven RP25 families. Herein we present 4 novel sequence variants, of which one is a novel deletion within a low complexity region close to the initiation codon of FAM46A. Furthermore, we have characterised for the first time a coding tandem variation in the Caucasian population. This study reports on bioinformatic and moleculardata for the FAM46A gene that may give a wider insight into the putative function of this gene and its pathologic relevance to RP25 and other retinal diseases mapping within the 6q chromosomal interval.

Mapping of transcription start sites of human retina expressed genes

BACKGROUND

The proper assembly of the transcriptional initiation machinery is a key regulatory step in the execution of the correct program of mRNA synthesis. The use of alternative transcription start sites (TSSs) provides a mechanism for cell and tissue specific gene regulation. Our knowledge of transcriptional initiation sequences in the human genome is limited despite the availability of the complete genome sequence. While genome wide experimental and bioinformatic approaches are improving our knowledge of TSSs, they lack information concerning genes expressed in a restricted manner or at very low levels, such as tissue specific genes.

RESULTS

In this study we describe the mapping of TSSs of genes expressed in human retina. Genes have been selected on the basis of their physiological or developmental role in this tissue. Our work combines in silico analysis of ESTs and known algorithm predictions together with their experimental validation via Cap-finder RACE. We report here the TSSs mapping of 54 retina expressed genes: we retrieved new sequences for 41 genes, some of which contain un-annotated exons. Results can be grouped into five categories, compared to the RefSeq; (i) TSS located in new first exons, (ii) splicing variation of the second exon, (iii) extension of the annotated first exon, (iv) shortening of the annotated first exon, (v) confirmation of previously annotated TSS.

CONCLUSION

In silico and experimental analysis of the transcripts proved to be essential for the ultimate mapping of TSSs. Our results highlight the necessity of a tissue specific approach to complete the existing gene annotation. The new TSSs and transcribed sequences are essential for further exploration of the promoter and other cis-regulatory sequences at the 5'end of genes.

Genomic annotation of 15,809 ESTs identified from pooled early gestation human eyes

To complement cDNA libraries from the human eye at early gestation and to discover candidate genes associated with early ocular development, we used freshly dissected human eyeballs from week 9-14 of gestation to construct the early human fetal eye cDNA library. A total of 15,809 clones were isolated and sequenced from the unamplified and unnormalized library. We screened 11,246 good-quality ESTs, leading to the identification of 5,534 nonredundant clusters. Among them, 4,010 (72%) genes matched in the human protein database (Ensembl). The remaining 28% (1,524) corresponded to potentially novel or previously unidentified ESTs. We used BLASTX to compare our EST data with eight organisms and found common expression of a high portion of genes: Caenorhabditis briggsae (26%), Caenorhabditis elegans (27%), Anopheles gambiae (37%), Drosophila melanogaster (32%), Danio rerio (42%), Fugu rubripes (49%), Rattus norvegicusvalitus (52%), and Mus musculus (59%). Nevertheless, 48% (2,680 of 5,534) of the genes expressed in the early developing eye were not shared with current NEIBank human eye cDNA data. In addition, eight known retinal disease genes existed in our ESTs. Among them, six (COL11A1, BBS5, PDE6B, OAT, VMD2, and PGK1) were conserved among the genomes of other organisms, indicating that our annotated EST set provides not only a valuable resource for gene discovery and functional genomic analysis but also for phylogenetic analysis. Our foremost early gestation human eye cDNA library could provide detailed comparisons across species to identify physiological functions of genes and to elucidate evolutionary mechanisms.

Comprehensive in silico functional specification of mouse retina transcripts

BACKGROUND

The retina is a well-defined portion of the central nervous system (CNS) that has been used as a model for CNS development and function studies. The full specification of transcripts in an individual tissue or cell type, like retina, can greatly aid the understanding of the control of cell differentiation and cell function. In this study, we have integrated computational bioinformatics and microarray experimental approaches to classify the tissue specificity and developmental distribution of mouse retina transcripts.

RESULTS

We have classified a set of retina-specific genes using sequence-based screening integrated with computational and retina tissue-specific microarray approaches. 33,737 non-redundant sequences were identified as retina transcript clusters (RTCs) from more than 81,000 mouse retina ESTs. We estimate that about 19,000 to 20,000 genes might express in mouse retina from embryonic to adult stages. 39.1% of the RTCs are not covered by 60,770 RIKEN full-length cDNAs. Through comparison with 2 million mouse ESTs, spectra of neural, retinal, late-generated retinal, and photoreceptor -enriched RTCs have been generated. More than 70% of these RTCs have data from biological experiments confirming their tissue-specific expression pattern. The highest-grade retina-enriched pool covered almost all the known genes encoding proteins involved in photo-transduction.

CONCLUSION

This study provides a comprehensive mouse retina transcript profile for further gene discovery in retina and suggests that tissue-specific transcripts contribute substantially to the whole transcriptome.