All microarray data was subjected to QC and ERCC spike-in assessments, and any failing samples were omitted from the analysis. Biological outliers were identified by comparing samples from related structures using hierarchical clustering and inter-array correlation measures. Data for samples passing QC were normalized in three steps: 1) “within-batch” normalization to the 75th percentile expression values; 2) “cross-batch” bias reduction using ComBat57 ; and 3) "cross-brain" normalization as in step 1. Differential expression assessments were done using template vector correlation, where 1="in group" and 0="not in group", or by measuring the fold change, defined as mean expression in category divided by mean expression elsewhere. False discovery rates were estimated using permutation tests (Suppl. Methods). WGCNA was performed on all neocortical samples using the standard method36 ,58 , and on germinal layers by defining a consensus module in the 15 and 16pcw brains59 , only including genes differentially expressed across these layers (5494 genes; ANOVA p<0.01, Benjamini-Hochberg adjusted). Gene list characterizations were made using a combination of module eigengene / representative gene expression, gene ontology enrichment using DAVID60 , and enrichment for known brain-related categories (i.e.,61 ,62 ) using userListEnrichment63 . Module C31 is depicted using VisANT64 : the top 250 gene-gene connections based on topological overlap are shown, with histone genes removed for clarity. Rostral-caudal areal gradient genes were identified as follows: first, the center of each neocortical region was identified at 21pcw in Euclidean coordinates; second, the rostral/caudal region position was estimated as an angle along the lateral face of the brain centered at the temporal/frontal lobe juncture (ordering lobes roughly as frontal, parietal, occipital, temporal; Fig. 5c); third, for each brain gene expression in each layer was (Pearson) correlated with this region position; and finally, genes with R>0.5 in all four brains were identified. A similar strategy was used to identify unbiased areal gradient genes (Suppl. Methods). Enrichment of haCNSs was determined using hypergeometric tests. Samples in all plots are ordered in an anatomically relevant manner. Unless otherwise noted, all p-values are Bonferroni corrected for multiple comparisons.
Miller J.A., Ding S.L., Sunkin S.M., Smith K.A., Ng L., Szafer A., Ebbert A., Riley Z.L., Aiona K., Arnold J.M., Bennet C., Bertagnolli D., Brouner K., Butler S., Caldejon S., Carey A., Cuhaciyan C., Dalley R.A., Dee N., Dolbeare T.A., Facer B.A., Feng D., Fliss T.P., Gee G., Goldy J., Gourley L., Gregor B.W., Gu G., Howard R.E., Jochim J.M., Kuan C.L., Lau C., Lee C.K., Lee F., Lemon T.A., Lesnar P., McMurray B., Mastan N., Mosqueda N.F., Naluai-Cecchini T., Ngo N.K., Nyhus J., Oldre A., Olson E., Parente J., Parker P.D., Parry S.E., Player A.S., Pletikos M., Reding M., Royall J.J., Roll K., Sandman D., Sarreal M., Shapouri S., Shapovalova N.V., Shen E.H., Sjoquist N., Slaughterbeck C.R., Smith M., Sodt A.J., Williams D., Zöllei L., Fischl B., Gerstein M.B., Geschwind D.H., Glass I.A., Hawrylycz M.J., Hevner R.F., Huang H., Jones A.R., Knowles J.A., Levitt P., Phillips J.W., Sestan N., Wohnoutka P., Dang C., Bernard A., Hohmann J.G, & Lein E.S. (2014). Transcriptional Landscape of the Prenatal Human Brain. Nature, 508(7495), 199-206.