Samples were prepared using the SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Clontech), according to manufacturer’s protocols and using 10 ng of input. Final libraries were evaluated for size using a Fragment Analyzer (Advanced Analytical Technologies) and quantified using both Qubit (Thermo Fisher) and qPCR (Roche LightCycler 480, KAPA Illumina Library Quantification Kit) before sequencing. RNA-seq 75 bp paired-end reads from Illumina were checked with FastQC and fastq_screen. Reads were mapped to human genome hg19 using TopHat v 2.0.13 with Ensembl annotation (GRCh37.75) in gtf format. FeatureCounts was used to obtain gene counts with default option. The gene counts were normalized with DESeq, and PCA plot was created with plotPCA function. Bar charts represent the average normalized raw reads for different PSC lines and primary samples. Gene ontology differently selected top features for FMG, pMGLs and differentiated NPCs, ran on the online Gene Ontology Consortium tool (geneontology.org). Hierarchical clustering was performed on quantile normalized FPKMs from the different GEO datasets. For HMDM and iPSDM, we used GSE55536. For the comparison to adult human microglia, and other brain cells,
