Transcriptome of embryonic and neonatal mouse cortex by high-throughput RNA sequencing

Expression of sorting nexin 12 is regulated in developing cerebral cortical neurons

The sorting nexins (SNXs) are a family of proteins functioning in diverse processes, including endocytosis, endosomal sorting, and endosomal signaling. Sorting nexin 12 (SNX12) is one of the SNXs family members; however, its function remains unknown. To clarify the function of SNX12, in this study, we first investigated the expression profiles in mice, particularly in the central nervous system (CNS), and then analyzed the functional role on neurite outgrowth. We found that SNX12 was widely expressed in the adult mouse CNS and that its expression level was higher in the cerebral cortex than in other examined regions. SNX12 expression was detected in the neurons but not the glial cells of the adult mouse cerebral cortex. In the fetal brain, SNX12 expression increased during the embryonic stage and gradually decreased after birth. Although the immunoreactivities of SNX12 were widespread in the cerebral cortical cells in the fetal brain, the immunopositive signals of SNX12 were more intense in the postmitotic neurons in the cortical plate than in the proliferating precursor cells in the ventricular zone, suggesting that SNX12 plays critical roles in the postmitotic neurons during cerebral cortical development. Furthermore, in mouse neuroblastoma and N1E-115 cells and rat primary cortical neurons, SNX12 expression was increased as neurite outgrowth progressed and the knockdown of SNX12 attenuated the outgrowth of neurites. These results suggest that SNX12 regulates neurite formation during cerebral cortical development.

Comparative analysis of neural transcriptomes and functional implication of unannotated intronic expression

Background

The transcriptome and its regulation bridge the genome and the phenome. Recent RNA-seq studies unveiled complex transcriptomes with previously unknown transcripts and functions. To investigate the characteristics of neural transcriptomes and possible functions of previously unknown transcripts, we analyzed and compared nine recent RNA-seq datasets corresponding to tissues/organs ranging from stem cell, embryonic brain cortex to adult whole brain.

Results

We found that the neural and stem cell transcriptomes share global similarity in both gene and chromosomal expression, but are quite different from those of liver or muscle. We also found an unusually high level of unannotated expression in mouse embryonic brains. The intronic unannotated expression was found to be strongly associated with genes annotated for neurogenesis, axon guidance, negative regulation of transcription, and neural transmission. These functions are the hallmarks of the late embryonic stage cortex, and crucial for synaptogenesis and neural circuit formation.

Conclusions

Our results revealed unique global and local landscapes of neural transcriptomes. It also suggested potential functional roles for previously unknown transcripts actively expressed in the developing brain cortex. Our findings provide new insights into potentially novel genes, gene functions and regulatory mechanisms in early brain development.

Sequence comparison of prefrontal cortical brain transcriptome from a tame and an aggressive silver fox (Vulpes vulpes)

Background

Two strains of the silver fox (Vulpes vulpes), with markedly different behavioral phenotypes, have been developed by long-term selection for behavior. Foxes from the tame strain exhibit friendly behavior towards humans, paralleling the sociability of canine puppies, whereas foxes from the aggressive strain are defensive and exhibit aggression to humans. To understand the genetic differences underlying these behavioral phenotypes fox-specific genomic resources are needed.

Results

cDNA from mRNA from pre-frontal cortex of a tame and an aggressive fox was sequenced using the Roche 454 FLX Titanium platform (> 2.5 million reads & 0.9 Gbase of tame fox sequence; >3.3 million reads & 1.2 Gbase of aggressive fox sequence). Over 80% of the fox reads were assembled into contigs. Mapping fox reads against the fox transcriptome assembly and the dog genome identified over 30,000 high confidence fox-specific SNPs. Fox transcripts for approximately 14,000 genes were identified using SwissProt and the dog RefSeq databases. An at least 2-fold expression difference between the two samples (p < 0.05) was observed for 335 genes, fewer than 3% of the total number of genes identified in the fox transcriptome.

Conclusions

Transcriptome sequencing significantly expanded genomic resources available for the fox, a species without a sequenced genome. In a very cost efficient manner this yielded a large number of fox-specific SNP markers for genetic studies and provided significant insights into the gene expression profile of the fox pre-frontal cortex; expression differences between the two fox samples; and a catalogue of potentially important gene-specific sequence variants. This result demonstrates the utility of this approach for developing genomic resources in species with limited genomic information.

Defective photoreceptor phagocytosis in a mouse model of enhanced S-cone syndrome causes progressive retinal degeneration

Enhanced S-cone syndrome (ESCS), featuring an excess number of S cones, manifests as a progressive retinal degeneration that leads to blindness. Here, through optical imaging, we identified an abnormal interface between photoreceptors and the retinal pigment epithelium (RPE) in 9 patients with ESCS. The neural retina leucine zipper transcription factor-knockout (Nrl(-/-)) mouse model demonstrates many phenotypic features of human ESCS, including unstable S-cone-positive photoreceptors. Using massively parallel RNA sequencing, we identified 6203 differentially expressed transcripts between wild-type (Wt) and Nrl(-/-) mouse retinas, with 6 highly significant differentially expressed genes of the Pax, Notch, and Wnt canonical pathways. Changes were also obvious in expression of 30 genes involved in the visual cycle and 3 key genes in photoreceptor phagocytosis. Novel high-resolution (100 nm) imaging and reconstruction of Nrl(-/-) retinas revealed an abnormal packing of photoreceptors that contributed to buildup of photoreceptor deposits. Furthermore, lack of phagosomes in the RPE layer of Nrl(-/-) retina revealed impairment in phagocytosis. Cultured RPE cells from Wt and Nrl(-/-) mice illustrated that the phagocytotic defect was attributable to the aberrant interface between ESCS photoreceptors and the RPE. Overcoming the retinal phagocytosis defect could arrest the progressive degenerative component of this disease.

The Role of FRMD7 in Idiopathic Infantile Nystagmus

Idiopathic infantile nystagmus (IIN) is an inherited disorder in which the nystagmus arises independently of any other symptoms, leading to the speculation that the disorder represents a primary defect in the area of the brain responsible for ocular motor control. The inheritance patterns are heterogeneous, however the most common form is X-linked. FRMD7 resides at Xq26-27 and approximately 50% of X-linked IIN families map to this region. Currently 45 mutations within FRMD7 have been associated with IIN, confirming the importance of FRMD7 in the pathogenesis of the disease. Although mutations in FRMD7 are known to cause IIN, very little is known about the function of the protein. FRMD7 contains a conserved N-terminal FERM domain suggesting that it may provide a link between the plasma membrane and actin cytoskeleton. Limited studies together with the knowledge of the function of other FERM domain containing proteins, suggest that FRMD7 may play a role in membrane extension during neuronal development through remodeling of the actin cytoskeleton.

Transcriptional programs in transient embryonic zones of the cerebral cortex defined by high-resolution mRNA sequencing

Characterizing the genetic programs that specify development and evolution of the cerebral cortex is a central challenge in neuroscience. Stem cells in the transient embryonic ventricular and subventricular zones generate neurons that migrate across the intermediate zone to the overlying cortical plate, where they differentiate and form the neocortex. It is clear that not one but a multitude of molecular pathways are necessary to progress through each cellular milestone, yet the underlying transcriptional programs remain unknown. Here, we apply differential transcriptome analysis on microscopically isolated cell populations, to define five transcriptional programs that represent each transient embryonic zone and the progression between these zones. The five transcriptional programs contain largely uncharacterized genes in addition to transcripts necessary for stem cell maintenance, neurogenesis, migration, and differentiation. Additionally, we found intergenic transcriptionally active regions that possibly encode unique zone-specific transcripts. Finally, we present a high-resolution transcriptome map of transient zones in the embryonic mouse forebrain.

A tissue-specific landscape of sense/antisense transcription in the mouse intestine

Background

The intestinal mucosa is characterized by complex metabolic and immunological processes driven highly dynamic gene expression programs. With the advent of next generation sequencing and its utilization for the analysis of the RNA sequence space, the level of detail on the global architecture of the transcriptome reached a new order of magnitude compared to microarrays.

Results

We report the ultra-deep characterization of the polyadenylated transcriptome in two closely related, yet distinct regions of the mouse intestinal tract (small intestine and colon). We assessed tissue-specific transcriptomal architecture and the presence of novel transcriptionally active regions (nTARs). In the first step, signatures of 20,541 NCBI RefSeq transcripts could be identified in the intestine (74.1% of annotated genes), thereof 16,742 are common in both tissues. Although the majority of reads could be linked to annotated genes, 27,543 nTARs not consistent with current gene annotations in RefSeq or ENSEMBL were identified. By use of a second independent strand-specific RNA-Seq protocol, 20,966 of these nTARs were confirmed, most of them in vicinity of known genes. We further categorized our findings by their relative adjacency to described exonic elements and investigated regional differences of novel transcribed elements in small intestine and colon.

Conclusions

The current study demonstrates the complexity of an archetypal mammalian intestinal mRNA transcriptome in high resolution and identifies novel transcriptionally active regions at strand-specific, single base resolution. Our analysis for the first time shows a strand-specific comparative picture of nTARs in two tissues and represents a resource for further investigating the transcriptional processes that contribute to tissue identity.

Functional genomics of the brain: uncovering networks in the CNS using a systems approach

The central nervous system (CNS) is undoubtedly the most complex human organ system in terms of its diverse functions, cellular composition, and connections. Attempts to capture this diversity experimentally were the foundation on which the field of neurobiology was built. Until now though, techniques were either painstakingly slow or insufficient in capturing this heterogeneity. In addition, the combination of multiple layers of information needed for a complete picture of neuronal diversity from the epigenome to the proteome requires an even more complex compilation of data. In this era of high-throughput genomics though, the ability to isolate and profile neurons and brain tissue has increased tremendously and now requires less effort. Both microarrays and next-generation sequencing have identified neuronal transcriptomes and signaling networks involved in normal brain development, as well as in disease. However, the expertise needed to organize and prioritize the resultant data remains substantial. A combination of supervised organization and unsupervised analyses are needed to fully appreciate the underlying structure in these datasets. When utilized effectively, these analyses have yielded striking insights into a number of fundamental questions in neuroscience on topics ranging from the evolution of the human brain to neuropsychiatric and neurodegenerative disorders. Future studies will incorporate these analyses with behavioral and physiological data from patients to more efficiently move toward personalized therapeutics.

SOBP is mutated in syndromic and nonsyndromic intellectual disability and is highly expressed in the brain limbic system

Intellectual disability (ID) affects 1%-3% of the general population. We recently reported on a family with autosomal-recessive mental retardation with anterior maxillary protrusion and strabismus (MRAMS) syndrome. One of the reported patients with ID did not have dysmorphic features but did have temporal lobe epilepsy and psychosis. We report on the identification of a truncating mutation in the SOBP that is responsible for causing both syndromic and nonsyndromic ID in the same family. The protein encoded by the SOBP, sine oculis binding protein ortholog, is a nuclear zinc finger protein. In mice, Sobp (also known as Jxc1) is critical for patterning of the organ of Corti; one of our patients has a subclinical cochlear hearing loss but no gross cochlear abnormalities. In situ RNA expression studies in postnatal mouse brain showed strong expression in the limbic system at the time interval of active synaptogenesis. The limbic system regulates learning, memory, and affective behavior, but limbic circuitry expression of other genes mutated in ID is unusual. By comparing the protein content of the +/jc to jc/jc mice brains with the use of proteomics, we detected 24 proteins with greater than 1.5-fold differences in expression, including two interacting proteins, dynamin and pacsin1. This study shows mutated SOBP involvement in syndromic and nonsyndromic ID with psychosis in humans.

Changes in molecular composition of rat medial prefrontal cortex synapses during adolescent development

Postnatal brain development continues throughout adolescence into young adulthood. In particular, synapse strengthening and elimination are prominent processes during adolescence. However, molecular data of this relatively late stage of synaptic development are sparse. In this study, we used iTRAQ (isobaric tag for relative and absolute quantification)-based proteomics and electron microscopy to investigate the molecular composition of a synaptic membrane fraction from adolescent postnatal day (P)34 and P44 and adult (P78) rat medial prefrontal cortex. Differential expression of proteins was most prominent between early adolescence and young adulthood (35%, P34-P78), with an over-representation of cell-membrane proteins during adolescent development (between P34 and P44), and synaptic vesicle proteins between late adolescence and young adulthood (P44-P78). Indicative of the critical period of development, we found that, between P34 and P44, a substantial number of proteins was differentially expressed (14%), much more than during the period after adolescence, i.e. between P44 and P78 (5%). A striking observation was the developmental non-stoichiometric regulation of distinct classes of proteins from the synaptic vesicle and the presynaptic release machinery. Electron microscopy demonstrated a small change in the number of docked vesicles between P34 and P44, but not in the total number of synaptic vesicles and in the size of the vesicle cluster. We conclude that the molecular composition of synapses, and more specifically the synaptic release machinery, of the medial prefrontal cortex changes drastically during adolescent development.

Transcriptome and expression profiling analysis revealed changes of multiple signaling pathways involved in immunity in the large yellow croaker during Aeromonas hydrophila infection

BACKGROUND

The large yellow croaker (Pseudosciaena crocea) is an economically important marine fish in China suffering from severe outbreaks of infectious disease caused by marine bacteria such as Aeromonas hydrophila (A. hydrophila), resulting in great economic losses. However, the mechanisms involved in the immune response of this fish to bacterial infection are not fully understood. To understand the molecular mechanisms underlying the immune response to such pathogenic bacteria, we used high-throughput deep sequencing technology to investigate the transcriptome and comparative expression profiles of the large yellow croaker infected with A. hydrophila.

RESULTS

A total of 13,611,340 reads were obtained and assembled into 26,313 scaffolds in transcriptional responses of the A. hydrophila-infected large yellow croaker. Via annotation to the NCBI database, we obtained 8216 identified unigenes. In total, 5590 (68%) unigenes were classified into Gene Ontology, and 3094 unigenes were found in 20 KEGG categories. These genes included representatives from almost all functional categories. By using Solexa/Illumina's DeepSAGE, 1996 differentially expressed genes (P value < 0.05) were detected in comparative analysis of the expression profiles between A. hydrophila-infected fish and control fish, including 727 remarkably upregulated genes and 489 remarkably downregulated genes. Dramatic differences were observed in genes involved in the inflammatory response. Bacterial infection affected the gene expression of many components of signaling cascades, including the Toll-like receptor, JAK-STAT, and MAPK pathways. Genes encoding factors involved in T cell receptor (TCR) signaling were also revealed to be regulated by infection in these fish.

CONCLUSION

Based on our results, we conclude that the inflammatory response may play an important role in the early stages of infection. The signaling cascades such as the Toll-like receptor, JAK-STAT, and MAPK pathways are regulated by A. hydrophila infection. Interestingly, genes encoding factors involved in TCR signaling were revealed to be downregulated by infection, indicating that TCR signaling was suppressed at this early period. These results revealed changes of multiple signaling pathways involved in immunity during A. hydrophila infection, which will facilitate our comprehensive understanding of the mechanisms involved in the immune response to bacterial infection in the large yellow croaker.

RNA-Seq Atlas of Glycine max: a guide to the soybean transcriptome

BACKGROUND

Next generation sequencing is transforming our understanding of transcriptomes. It can determine the expression level of transcripts with a dynamic range of over six orders of magnitude from multiple tissues, developmental stages or conditions. Patterns of gene expression provide insight into functions of genes with unknown annotation.

RESULTS

The RNA Seq-Atlas presented here provides a record of high-resolution gene expression in a set of fourteen diverse tissues. Hierarchical clustering of transcriptional profiles for these tissues suggests three clades with similar profiles: aerial, underground and seed tissues. We also investigate the relationship between gene structure and gene expression and find a correlation between gene length and expression. Additionally, we find dramatic tissue-specific gene expression of both the most highly-expressed genes and the genes specific to legumes in seed development and nodule tissues. Analysis of the gene expression profiles of over 2,000 genes with preferential gene expression in seed suggests there are more than 177 genes with functional roles that are involved in the economically important seed filling process. Finally, the Seq-atlas also provides a means of evaluating existing gene model annotations for the Glycine max genome.

CONCLUSIONS

This RNA-Seq atlas extends the analyses of previous gene expression atlases performed using Affymetrix GeneChip technology and provides an example of new methods to accommodate the increase in transcriptome data obtained from next generation sequencing. Data contained within this RNA-Seq atlas of Glycine max can be explored at http://www.soybase.org/soyseq.

Exploring plant transcriptomes using ultra high-throughput sequencing

Ultra high-throughput sequencing (UHTS) technologies offer the potential to interrogate transcriptomes in detail that has traditionally been restricted to single gene surveys. For instance, it is now possible to globally define transcription start sites, polyadenylation signals, alternative splice sites and generate quantitative data on gene transcript accumulation in single tissues or cell types. These technologies are thus paving the way for whole genome transcriptomics and will undoubtedly lead to novel insights into plant development and biotic and abiotic stress responses. However, several challenges exist to making this technology broadly accessible to the plant research community. These include the current need for a computationally intensive analysis of data sets, a lack of standardized alignment and formatting procedures and a relatively small number of analytical software packages to interpret UHTS outputs. In this review we summarize recent findings from UHTS and discuss potential opportunities and challenges for broad adoption of these technologies in the plant science community.