An RNA-Sequencing Transcriptome and Splicing Database of Glia, Neurons, and Vascular Cells of the Cerebral Cortex

Comparison of RNA-Seq and microarray analyses showed that RNA-Seq identified more differentially expressed genes across cell types than microarray. A, Spearman's rank correlation between gene expression data obtained by the RNA-Seq and microarray methods across cell types. Gene expression profiles of the same cell type obtained by the RNA-Seq and microarray methods showed a high degree of correlation. B, Numbers of differentially expressed genes identified by the RNA-Seq and microarray methods. A snapshot of the data is summarized in the Venn diagram. Using a fourfold difference as the cutoff, RNA-Seq analysis identified 3129 genes as differentially expressed by astrocytes, whereas microarray identified only 1367 genes. The majority of genes identified by microarray as differentially expressed are similarly classified by RNA-Seq (1279). RNA-Seq identified an additional 1850 genes as differentially expressed that were not identified by microarray as differentially expressed. C–E, RNA-Seq versus microarray comparisons for neurons, oligodendrocytes, and microglia. F, The relationship between fold enrichment and expression level in astrocytes. There is a sharp cutoff line on the left (FPKM = 0.1 or −3.3 on the log2 scale) because we set any FPKM value <0.1 to 0.1 to avoid ratio inflation in fold enrichment calculations. G, The relationship between fold enrichment and expression level in astrocytes.

RNA-Seq analysis revealed cell type-specific gene expression profiles. A, Dendrogram and unsupervised hierarchical clustering heat map (using Euclidean distance) of purified cortical glia, neurons, and vascular cells. The vertical distances on each branch of the dendrogram represent the degree of similarity between gene expression profiles of various samples. Biological replicates showed the highest degree of correlation represented by short vertical distances. Cells in the oligodendrocyte lineage cluster closely together, and the order of the three oligodendrocyte-lineage cell types corresponds to their maturation stages (OPC–NFO–MO). Although astrocytes and oligodendrocytes are both glial cells, their gene expression profiles are as different between each other as they are different from neurons. Consistent with their embryonic origin, mesodermal-derived endothelial cells and microglia cluster farther away from ectoderm-derived neurons, astrocytes, and oligodendrocytes. B, The top 40 enriched genes per cell type are shown in a heat map. Only highly expressed genes with FPKM >20 are included in this analysis. Fold enrichment is calculated as FPKM of one cell type divided by the average FPKM of all other cell types. The majority of these genes showed specific expression by only one cell type, with the exception that some are expressed during more than one maturation stage in the oligodendrocyte lineage.

Cell-specific markers and transcription factors. The top 40 genes and top 10 transcription factors ranked by fold enrichment of each cell type are listed. The most highly expressed genes are highlighted. Green, FPKM >150; blue, FPKM >50 for transcription factors. Fold enrichment of astrocytes, neurons, microglia, and endothelial cells are calculated as FPKM of one cell type divided by the average FPKM of all other cell types. Fold enrichment of OPC, NFO, and MO are calculated as FPKM of one cell type divided by the average FPKM of all non-oligodendrocyte-lineage cells, to highlight top genes specifically expressed by a particular maturation stage during oligodendrocyte development. Only highly expressed genes with FPKM >20 are included in the ranking to highlight genes that are most likely to have significant cell type-specific functions.

Alternative splicing analysis of neurons, glia and vascular cells. A, Six types of alternative splicing events are detected by RNA-Seq. Boxes and black lines represent exons and introns, respectively. Blue lines and red lines represent alternative splicing events detected in the dataset. The 5′ end is to the left, and the 3′ end is to the right. Cassette, The inclusion or exclusion of an exon; Tandom Cassette, the inclusion or exclusion of two or more tandom exons; Mutually Exclusive, the inclusion of one exon in one transcript and inclusion of a different exon in another transcript; Intron Retention, the inclusion or exclusion of a segment previously annotated to be an intron; Alternative 5′ SS, the alternative usage of a splicing site on the 5′ end of an exon; Alternative 3′ SS, the alternative usage of a splicing site on the 3′ end of an exon. B, Frequencies of the six types of alternative splicing events detected in the entire dataset and in individual cell types. In all cell types, cassette exon events, i.e., the inclusion or exclusion of an exon, are the most frequently detected alternative splicing events. C, The numbers of genes that are alternatively spliced in each cell type and the union of these samples. The dotted line represents the total number of genes that are known to contain a potential splicing event in the mouse genome. The number of these genes that are expressed in a given cell type are represented by gray bars. The black bars indicate the numbers of genes that are alternatively spliced in a given cell type based on criteria outlined in Materials and Methods. D, The numbers of statistically significant cell type-specific alternative splicing events in each cell type. Neurons have the highest number of specific splicing events, whereas oligodendrocyte-lineage cells have the least amount of specific splicing events. E, Pkm2 is an example of a gene spliced uniquely in astrocytes and neurons. The traces represent raw data of the number of reads mapped to the Pkm2 gene from astrocytes and neurons. The height of the blue bars represents number of reads. The bottom schematic is the transcript model of Pkm2 gene from the UCSC Genome Browser. Boxes represent exons, and black lines represent introns. The exon shown in blue is predominantly included in neurons, whereas the exon shown in yellow is included only in astrocytes. This is an example of a mutually exclusive event. F, Validation of PKM1/2 splicing differences by PCR. We designed primers targeting exons unique to PKM1, unique to PKM2, and exons common to PKM1 and PKM2. PCR products were detected from neuron, astrocyte, and whole-brain samples in patterns predicted by the RNA-Seq data.

Validation of RNA-Seq results by qRT-PCR and in situ hybridization. A, qRT-PCR validation of cell type-enriched genes identified by RNA-Seq. We performed qRT-PCR with Fluidigm BioMark microfluidic technology. Expression of several genes identified by RNA-Seq as enriched in each cell type was examined by qRT-PCR. The housekeeping gene Gapdh was included for comparison. Twelve replicates of each purified cell type and three replicates of whole-brain samples were analyzed. Warmer colors represent lower Ct values (higher abundance of transcripts), and cooler colors represent higher Ct values (lower abundance of transcripts). Black indicates no amplification. Data of genes labeled in red were quantified in B. B, Ct differences of Atp13a4, Cpne7, Fam70b, Tmem88b, and Rcsd1 compared with Gapdh were plotted on a log2 scale. Error bar represents SD. RNA-Seq analysis showed that Atp13a4, Cpne7, Fam70b, Tmem88b, and Rcsd1 are enriched in astrocytes, neurons, OPCs, oligodendrocytes, and microglia, respectively. qRT-PCR validated these results. C–G, In situ hybridization validated novel cell type-specific genes. Left, Low-magnification image of the cortex and hippocampus. Scale bar, 200 μm. Right, High-magnification image of the cortex. Scale bar, 50 μm. Fluorescence in situ hybridization signals with probes against novel cell type-specific genes (red) and known cell type-specific markers (green) are shown. The regions in the yellow boxes are enlarged and shown as single channel and merged images on the right.