Summary

Here we report the generation and analysis of genome-wide exon-level transcriptome data from 16 brain regions comprising the cerebellar cortex, mediodorsal nucleus of the thalamus, striatum, amygdala, hippocampus, and 11 areas of the neocortex. The dataset was generated from 1,340 tissue samples collected from one or both hemispheres of 57 postmortem human brains, spanning from embryonic development to late adulthood and representing males and females of multiple ethnicities. We also performed genotyping of 2.5 million SNPs and assessed copy number variations for all donors. Approximately 86% of protein-coding genes were found to be expressed using stringent criteria, and over 90% of these were differentially regulated at the whole transcript or exon level across regions and/or time. The majority of these spatiotemporal differences occurred before birth, followed by an increase in the similarity among regional transcriptomes during postnatal lifespan. Genes were organized into functionally distinct co-expression networks, and sex differences were present in gene expression and exon usage. Finally, we demonstrate how these results can be used to profile trajectories of genes associated with neurodevelopmental processes, cell types, neurotransmitter systems, autism, and schizophrenia, as well as to discover associations between SNPs and spatiotemporal gene expression. This study provides a comprehensive, publicly available dataset on the spatiotemporal human brain transcriptome and new insights into the transcriptional foundations of human neurodevelopment.

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.

Background

Age-related macular degeneration (AMD), a chronic neurodegenerative and neovascular retinal disease, is the leading cause of blindness in elderly people of western European origin. While structural and functional alterations in mitochondria (mt) and their metabolites have been implicated in the pathogenesis of chronic neurodegenerative and vascular diseases, the relationship of inherited variants in the mitochondrial genome and mt haplogroup subtypes with advanced AMD has not been reported in large prospective cohorts.

Methodology/Prinicipal Findings

We examined the relationship of inherited mtDNA variants with advanced AMD in 1168 people using a three-stage design on samples from 12-year and 10-year prospective studies on the natural history of age-related eye disease. In Stage I we resequenced the entire genome in 99 elderly AMD-free controls and 215 people with advanced AMD from the 12-year study. A consistent association with AMD in 14 of 17 SNPs characterizing the mtDNA T haplogroup emerged. Further analysis revealed these associations were driven entirely by the T2 haplogroup, and characterized by two variants in Complex I genes (A11812G of MT-ND4 and A14233G of MT-ND6). We genotyped T haplogroups in an independent sample of 490 cases and 61 controls from the same study (Stage II) and in 56 cases and 246 controls from the 10-year study (Stage III). People in the T2 haplogroup were approximately 2.5 times more likely to have advanced AMD than their peers (odds ratio [OR] = 2.54, 95%CI 1.36–4.80, P≤0.004) after considering the totality of evidence. Findings persisted after considering the impact of AMD-associated variants A69S and Y402H (OR = 5.19, 95%CI 1.19–22.69, P≤0.029).

Conclusion

Loci defining the mtDNA T2 haplogroup and Complex I are reasonable targets for novel functional analyses and therapeutic research in AMD.

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.

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.

Background

Pigment epithelium derived factor (PEDF), a member of the serpin family, regulates cell proliferation, promotes survival of neurons, and blocks growth of new blood vessels in mammals. Defining the molecular phylogeny of PEDF by bioinformatic analysis is one approach to understanding the link between its gene structure and its function in these biological processes.

Results

From a comprehensive search of available DNA databases we identified a single PEDF gene in all vertebrate species examined. These included four mammalian and six non-mammalian vertebrate species in which PEDF had not previously been described. A five gene cluster around PEDF was found in an approximate 100 kb region in mammals, birds, and amphibians. In ray-finned fish these genes are scattered over three chromosomes although only one PEDF gene was consistently found. The PEDF gene is absent in invertebrates including Drosophila melanogaster (D. melanogaster), Caenorhabditis elegans (C. elegans), and sea squirt (C. intestinalis). The PEDF gene is transcribed in all vertebrate phyla, suggesting it is biologically active throughout vertebrate evolution. The multiple actions of PEDF are likely conserved in evolution since it has the same gene structure across phyla, although the size of the gene ranges from 48.3 kb in X. tropicalis to 2.9 kb in fugu, with human PEDF at a size of 15.6 kb. A strong similarity in the proximal 200 bp of the PEDF promoter in mammals suggests the existence of a possible regulatory region across phyla. Using a non-synonymous/synonymous substitution rate ratio we show that mammalian and fish PEDFs have similar ratios of <0.13, reflecting a strong purifying selection of PEDF gene. A large number of repetitive transposable elements of the SINE and LINE class were found with random distribution in both the promoter and introns of mammalian PEDF.

Conclusion

The PEDF gene first appears in vertebrates and our studies suggest that the regulation and biological actions of this gene are preserved across vertebrates. This comprehensive analysis of the PEDF gene across phyla provides new information that will aid further characterization of common functional motifs of this serpin in biological processes.

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.