Comprehensive in silico functional specification of mouse retina transcripts

Identification and prediction of alternative transcription start sites that generate rod photoreceptor-specific transcripts from ubiquitously expressed genes

Transcriptome complexity is substantially increased by the use of multiple transcription start sites for a given gene. By utilizing a rod photoreceptor-specific chromatin signature, and the RefSeq database of established transcription start sites, we have identified essentially all known rod photoreceptor genes as well as a group of novel genes that have a high probability of being expressed in rod photoreceptors. Approximately half of these novel rod genes are transcribed into multiple mRNA and/or protein isoforms through alternative transcriptional start sites (ATSS), only one of which has a rod-specific epigenetic signature and gives rise to a rod transcript. This suggests that, during retina development, some genes use ATSS to regulate cell type and temporal specificity, effectively generating a rod transcript from otherwise ubiquitously expressed genes. Biological confirmation of the relationship between epigenetic signatures and gene expression, as well as comparison of our genome-wide chromatin signature maps with available data sets for retina, namely a ChIP-on-Chip study of Polymerase-II (Pol-II) binding sites, ChIP-Seq studies for NRL- and CRX- binding sites and DHS (University of Washington data, available on UCSC mouse Genome Browser as a part of ENCODE project) fully support our hypothesis and together accurately identify and predict an array of new rod transcripts. The same approach was used to identify a number of TSS that are not currently in RefSeq. Biological conformation of the use of some of these TSS suggests that this method will be valuable for exploring the range of transcriptional complexity in many tissues. Comparison of mouse and human genome-wide data indicates that most of these alternate TSS appear to be present in both species, indicating that our approach can be useful for identification of regulatory regions that might play a role in human retinal disease.

Unravelling the genetics of inherited retinal dystrophies: Past, present and future

The identification of the genes underlying monogenic diseases has been of interest to clinicians and scientists for many years. Using inherited retinal dystrophies as an example of monogenic disease we describe the history of molecular genetic techniques that have been pivotal in the discovery of disease causing genes. The methods that were developed in the 1970's and 80's are still in use today but have been refined and improved. These techniques enabled the concept of the Human Genome Project to be envisaged and ultimately realised. When the successful conclusion of the project was announced in 2003 many new tools and, as importantly, many collaborations had been developed that facilitated a rapid identification of disease genes. In the post-human genome project era advances in computing power and the clever use of the properties of DNA replication has allowed the development of next-generation sequencing technologies. These methods have revolutionised the identification of disease genes because for the first time there is no need to define the position of the gene in the genome. The use of next generation sequencing in a diagnostic setting has allowed many more patients with an inherited retinal dystrophy to obtain a molecular diagnosis for their disease. The identification of novel genes that have a role in the development or maintenance of retinal function is opening up avenues of research which will lead to the development of new pharmacological and gene therapy approaches. Neither of which can be used unless the defective gene and protein is known. The continued development of sequencing technologies also holds great promise for the advent of truly personalised medicine.

Microarray and morphological analysis of early postnatal CRB2 mutant retinas on a pure C57BL/6J genetic background

In humans, the Crumbs homologue-1 (CRB1) gene is mutated in progressive types of autosomal recessive retinitis pigmentosa and Leber congenital amaurosis. The severity of the phenotype due to human CRB1 or mouse Crb1 mutations is dependent on the genetic background. Mice on C57BL/6J background with Crb1 mutations show late onset of retinal spotting phenotype or no phenotype. Recently, we showed that conditional deletion of mouse Crb2 in the retina results in early retinal disorganization leading to severe and progressive retinal degeneration with concomitant visual loss that mimics retinitis pigmentosa due to mutations in the CRB1 gene. Recent studies in the fruit fly and zebrafish suggest roles of the Crumbs (CRB) complex members in the regulation of cellular signalling pathways including the Notch1, mechanistic target of rapamycin complex 1 (mTORC1) and the Hippo pathway. Here, we demonstrate that mice backcrossed to C57BL/6J background with loss of CRB2 in the retina show a progressive disorganization and degeneration phenotype during late retinal development. We used microarray gene profiling to study the transcriptome of retinas lacking CRB2 during late retinal development. Unexpectedly, the retinas of newborn mice lacking CRB2 showed no changes in the transcriptome during retinal development. These findings suggest that loss of CRB2 in the developing retina results in retinal disorganization and subsequent degeneration without major changes in the transcriptome of the retina. These mice might be an interesting model to study the onset of retinal degeneration upon loss of CRB proteins.

CD4 positive T helper cells contribute to retinal ganglion cell death in mouse model of ischemia reperfusion injury

Neuron degeneration is a common pathological process associated with many disease conditions in the central nervous system including retina. Although immune responses have been proposed as one potential element in triggering neural damage, the mechanism of action of specific immune components underlying the pathogenesis is unclear. In this study we focus on adaptive immune activities to evaluate CD4 positive helper cells in the retinal ganglion cell (RGC) degeneration in response to transient retinal ischemic/reperfusion (I/R) injury. Transient retinal ischemia was induced in four mouse strains with different immune backgrounds, including wild type mice from C57BL/6 and BABL/c strains, severe combined immunodeficient (SCID) mice lacking T and B lymphocytes, SCID mice with transferred wild type CD4+ T cells, and the STAT6 deficient mice without T helper 2 (TH2) cells. In SCID mice RGCs showed a strong resistance to cell death in response to I/R injury (89% ± 3% of the survival cells in contralateral eye) compared with C57BL/6 (p = 0.018) and BALB/C (p = 0.038) wild types. By transferring the mature CD4+ T cells from matched wild type into SCID mice, the resistance of RGCs to injury was significantly compromised (p < 0.05). Furthermore a significant resistance of RGCs to cell death (p < 0.05) accompanied with an overexpression of STAT1 and STAT3 was confirmed in STAT6 deficient mice in response to I/R injury compared with the wild type controls, indicating that TH2 cells maturation might be involved in RGC damage. Adaptive immunity carried by CD4 T cells plays an essential role in RGC degeneration.

Stage and gene specific signatures defined by histones H3K4me2 and H3K27me3 accompany mammalian retina maturation in vivo

The epigenetic contribution to neurogenesis is largely unknown. There is, however, growing evidence that posttranslational modification of histones is a dynamic process that shows many correlations with gene expression. Here we have followed the genome-wide distribution of two important histone H3 modifications, H3K4me2 and H3K27me3 during late mouse retina development. The retina provides an ideal model for these studies because of its well-characterized structure and development and also the extensive studies of the retinal transcriptome and its development. We found that a group of genes expressed only in mature rod photoreceptors have a unique signature consisting of de-novo accumulation of H3K4me2, both at the transcription start site (TSS) and over the whole gene, that correlates with the increase in transcription, but no accumulation of H3K27me3 at any stage. By in silico analysis of this unique signature we have identified a larger group of genes that may be selectively expressed in mature rod photoreceptors. We also found that the distribution of H3K4me2 and H3K27me3 on the genes widely expressed is not always associated with their transcriptional levels. Different histone signatures for retinal genes with the same gene expression pattern suggest the diversities of epigenetic regulation. Genes without H3K4me2 and H3K27me3 accumulation at any stage represent a large group of transcripts never expressed in retina. The epigenetic signatures defined by H3K4me2 and H3K27me3 can distinguish cell-type specific genes from widespread transcripts and may be reflective of cell specificity during retina maturation. In addition to the developmental patterns seen in wild type retina, the dramatic changes of histone modification in the retinas of mutant animals lacking rod photoreceptors provide a tool to study the epigenetic changes in other cell types and thus describe a broad range of epigenetic events in a solid tissue in vivo.

Transcription factories in the context of the nuclear and genome organization

In the eukaryotic nucleus, genes are transcribed in transcription factories. In the present review, we re-evaluate the models of transcription factories in the light of recent and older data. Based on this analysis, we propose that transcription factories result from the aggregation of RNA polymerase II-containing pre-initiation complexes assembled next to each other in the nuclear space. Such an aggregation can be triggered by the phosphorylation of the C-terminal domain of RNA polymerase II molecules and their interaction with various transcription factors. Individual transcription factories would thus incorporate tissue-specific, co-regulated as well as housekeeping genes based only on their initial proximity to each other in the nuclear space. Targeting genes to be transcribed to protein-dense factories that contain all factors necessary for transcription initiation and elongation through chromatin templates clearly favors a more economical utilization and better recycling of the transcription machinery.

Genetic networks in the mouse retina: growth associated protein 43 and phosphatase tensin homolog network

PURPOSE

The present study examines the structure and covariance of endogenous variation in gene expression across the recently expanded family of C57BL/6J (B) X DBA/2J (D) Recombinant Inbred (BXD RI) strains of mice. This work is accompanied by a highly interactive database that can be used to generate and test specific hypotheses. For example, we define the genetic network regulating growth associated protein 43 (Gap43) and phosphatase tensin homolog (Pten).

METHODS

The Hamilton Eye Institute (HEI) Retina Database within GeneNetwork features the data analysis of 346 Illumina Sentrix BeadChip Arrays (mouse whole genome-6 version 2). Eighty strains of mice are presented, including 75 BXD RI strains, the parental strains (C57BL/6J and DBA/2J), the reciprocal crosses, and the BALB/cByJ mice. Independent biologic samples for at least two animals from each gender were obtained with a narrow age range (48 to 118 days). Total RNA was prepared followed by the production of biotinylated cRNAs, which were pipetted into the Mouse WG-6V2 arrays. The data was globally normalized with rank invariant and stabilization (2z+8).

RESULTS

The HEI Retina Database is located on the GeneNetwork website. The database was used to extract unique transcriptome signatures for specific cell types in the retina (retinal pigment epithelial, amacrine, and retinal ganglion cells). Two genes associated with axonal outgrowth (Gap43 and Pten) were used to display the power of this new retina database. Bioinformatic tools located within GeneNetwork in conjunction with the HEI Retina Database were used to identify the unique signature Quantitative Trait Loci (QTLs) for Gap43 and Pten on chromosomes 1, 2, 12, 15, 16, and 19. Gap43 and Pten possess networks that are similar to ganglion cell networks that may be associated with axonal growth in the mouse retina. This network involves high correlations of transcription factors (SRY sex determining region Y-box 2 [Sox2], paired box gene 6 [Pax6], and neurogenic differentiation 1 [Neurod1]), and genes involved in DNA binding (proliferating cell nuclear antigen [Pcna] and zinc finger, BED-type containing 4 [Zbed4]), as well as an inhibitor of DNA binding (inhibitor of DNA binding 2, dominant negative helix-loop-helix protein [Id2]). Furthermore, we identified the potential upstream modifiers on chromosome 2 (teashirt zinc finger homeobox 2 [Tshz2], RNA export 1 homolog [Rae1] and basic helix-loop-helix domain contatining, class B4 [Bhlhb4]) on chromosome 15 (RAB, member of RAS oncogene family-like 2a [Rabl2a], phosphomannomutase 1 [Pmm1], copine VIII [Cpne8], and fibulin 1 [Fbln1]).

CONCLUSIONS

The endogenous variation in mRNA levels among BXD RI strains can be used to explore and test expression networks underlying variation in retina structure, function, and disease susceptibility. The Gap43 and Pten network highlights the covariance of gene expression and forms a molecular network associated with axonal outgrowth in the adult retina.

Comparison of gene expression during in vivo and in vitro postnatal retina development

Retina explants are widely used as a model of neural development. To define the molecular basis of differences between the development of retina in vivo and in vitro during the early postnatal period, we carried out a series of microarray comparisons using mouse retinas. About 75% of 8,880 expressed genes from retina explants kept the same expression volume and pattern as the retina in vivo. Fewer than 6% of the total gene population was changed at two consecutive time points, and only about 1% genes showed more than a threefold change at any time point studied. Functional Gene Ontology (GO) mapping for both changed and unchanged genes showed similar distribution patterns, except that more genes were changed in the GO clusters of response to stimuli and carbohydrate metabolism. Three distinct expression patterns of genes preferentially expressed in rod photoreceptors were observed in the retina explants. Some genes showed a lag in increased expression, some showed no change, and some continued to have a reduced level of expression. An early downregulation of cyclin D1 in the explanted retina might explain the reduction in numbers of precursors in explanted retina and suggests that external factors are required for maintenance of cyclin D1. The global view of gene profiles presented in this study will help define the molecular changes in retina explants over time and will provide criteria to define future changes that improve this model system.

ELECTRONIC SUPPLEMENTARY MATERIAL

The online version of this article (doi:10.1007/s12177-008-9009-z) contains supplementary material, which is available to authorized users.

Identification of molecular markers of bipolar cells in the murine retina

Retinal bipolar neurons serve as relay interneurons that connect rod and cone photoreceptor cells to amacrine and ganglion cells. They exhibit diverse morphologies essential for correct routing of photoreceptor cell signals to specific postsynaptic amacrine and ganglion cells. The development and physiology of these interneurons have not been completely defined molecularly. Despite previous identification of genes expressed in several bipolar cell subtypes, molecules that mark each bipolar cell type still await discovery. In this report, novel genetic markers of murine bipolar cells were found. Candidates were initially generated by using microarray analysis of single bipolar cells and mining of retinal serial analysis of gene expression (SAGE) data. These candidates were subsequently tested for expression in bipolar cells by RNA in situ hybridization. Ten new molecular markers were identified, five of which are highly enriched in their expression in bipolar cells within the adult retina. Double-labeling experiments using probes for previously characterized subsets of bipolar cells were performed to identify the subtypes of bipolar cells that express the novel markers. Additionally, the expression of bipolar cell genes was analyzed in Bhlhb4 knockout retinas, in which rod bipolar cells degenerate postnatally, to delineate further the identity of bipolar cells in which novel markers are found. From the analysis of Bhlhb4 mutant retinas, cone bipolar cell gene expression appears to be relatively unaffected by the degeneration of rod bipolar cells. Identification of molecular markers for the various subtypes of bipolar cells will lead to greater insights into the development and function of these diverse interneurons.

Lim1 is essential for the correct laminar positioning of retinal horizontal cells

Although much is known about the transcriptional regulation that coordinates retinal cell fate determination, very little is known about the developmental processes that establish the characteristic laminar architecture of the retina, in particular, the specification of neuronal positioning. The LIM class homeodomain transcription factor Lim1 (Lhx1) is expressed in postmitotic, differentiating, and mature retinal horizontal cells. We show that conditional ablation of Lim1 results in the ectopic localization of horizontal cells to inner aspects of the inner nuclear layer, among the retinal amacrine cells. The ectopic cells maintain a molecular phenotype consistent with horizontal cell identity; however, these neurons adopt a unique morphology more reminiscent of amacrine cells, including a dendritic arbor positioned within the inner plexiform layer. All other retinal cell populations appear unaltered. Our data suggest a model whereby Lim1 lies downstream of horizontal cell fate determination factors and functions cell autonomously to instruct differentiating horizontal cells to the appropriate laminar position in the developing retina. This study is the first to describe a cell type-specific genetic program that is essential for targeting a discrete retinal neuron population to the proper lamina.

A biphasic pattern of gene expression during mouse retina development

BACKGROUND

Between embryonic day 12 and postnatal day 21, six major neuronal and one glia cell type are generated from multipotential progenitors in a characteristic sequence during mouse retina development. We investigated expression patterns of retina transcripts during the major embryonic and postnatal developmental stages to provide a systematic view of normal mouse retina development,

RESULTS

A tissue-specific cDNA microarray was generated using a set of sequence non-redundant EST clones collected from mouse retina. Eleven stages of mouse retina, from embryonic day 12.5 (El2.5) to postnatal day 21 (PN21), were collected for RNA isolation. Non-amplified RNAs were labeled for microarray experiments and three sets of data were analyzed for significance, hierarchical relationships, and functional clustering. Six individual gene expression clusters were identified based on expression patterns of transcripts through retina development. Two developmental phases were clearly divided with postnatal day 5 (PN5) as a separate cluster. Among 4,180 transcripts that changed significantly during development, approximately 2/3 of the genes were expressed at high levels up until PN5 and then declined whereas the other 1/3 of the genes increased expression from PN5 and remained at the higher levels until at least PN21. Less than 1% of the genes observed showed a peak of expression between the two phases. Among the later increased population, only about 40% genes are correlated with rod photoreceptors, indicating that multiple cell types contributed to gene expression in this phase. Within the same functional classes, however, different gene populations were expressed in distinct developmental phases. A correlation coefficient analysis of gene expression during retina development between previous SAGE studies and this study was also carried out.

CONCLUSION

This study provides a complementary genome-wide view of common gene dynamics and a broad molecular classification of mouse retina development. Different genes in the same functional clusters are expressed in the different developmental stages, suggesting that cells might change gene expression profiles from differentiation to maturation stages. We propose that large-scale changes in gene regulation during development are necessary for the final maturation and function of the retina.

Proteomic trajectory mapping of biological transformation: Application to developmental mouse retina

In this report we introduce a new concept "proteomic trajectory mapping" for the investigation of a complex phenomenon underlying biological transformation and transition. We define proteomic trajectory to be the kinetic trace of protein expression and present a successful proteomic trajectory mapping of complex molecular events underlying postnatal development of mouse retina. Cluster analysis of the trajectory data using a two-state model identified four proteomic trajectory types: two distinct trajectory types accounting for the decline or the rise of protein molecules actively expressed in the juvenile stage (J-type) or in the adult stage (A-type), a class of transient trajectories that mediate the transformation from the juvenile to the adult stage (T-type), and the steady trajectories throughout the entire process of transformation (C-type). The dominance of particular protein categories expressed in each trajectory characterizes the stage of retinal development. Proteomic trajectory mapping will be a powerful tool to study the systematic changes of protein expression caused by physiological, genetic, or pathological agents and the reverse of such changes to the norm by a treatment. The proteomic trajectory mapping is applicable to any biological transformation and, therefore, will be a powerful tool in biomedical sciences.

Gene expression profiles of mouse retinas during the second and third postnatal weeks

Mouse retina undergoes crucial changes during early postnatal development. By using Affymetrix microarrays, we analyzed gene expression profiles of wild-type 129SvEv/C57BL/6 mouse retinas at postnatal days (P) 7, 10, 14, 18, and 21 and found significantly altered expression of 355 genes. Characterization of these 355 genes provided insight into physiologic and pathologic processes of mouse retinal development during the second and third postnatal weeks, a period that corresponds to human embryogenesis between weeks 12 and 28. These genes formed 6 groups with similar change patterns. Among the genes, sixteen cause retinal diseases when mutated; most of these 16 genes were upregulated in retina during this period. Using the PathArt program, we identified the biological processes in which many of the 355 gene products function. Among the most active processes in the P7-P21 retina are those involved in neurogenesis, obesity, diabetes type II, apoptosis, growth and differentiation, and protein kinase activity. We examined the expression patterns of 58 genes in P7 and adult retinas by searching the Brain Gene Expression Map database. Although most genes were present in various cell types in retinas, many displayed high levels of expression specifically in the outer nuclear, inner nuclear, and/or ganglion cell layers. By combining our 3 analyses, we demonstrated that during this period of mouse retinal development, many genes play important roles in various cell types, multiple pathways are involved, and some genes in a pathway are expressed in coordinated patterns. Our results thus provide foundation for future detailed studies of specific genes and pathways in various genetic and environmental conditions during retinal development.