Induced pluripotent stem cells (iPSCs), and cells derived from them, have become key tools to model biological processes, particularly in cell types that are difficult to access from living donors. We present the first map of regulatory variants in iPSC-derived neurons, based on 123 differentiations of iPSCs to a sensory neuronal fate. Gene expression was more variable across cultures than in primary dorsal root ganglion, particularly in genes related to nervous system development. Using single-cell RNA-sequencing, we found that the fraction of neuronal vs. contaminating cells was influenced by iPSC culture conditions prior to differentiation. Despite high differentiation-induced variability, using an allele-specific method we detected thousands of quantitative trait loci (QTLs) influencing gene expression, chromatin accessibility, and RNA splicing. Based on our QTLs, we estimate that recall-by-genotype studies using iPSC-derived cells will require at least 20-80 individuals to detect the effects of regulatory variants with moderately large effect sizes.

Highlights d Deconvolution of in vitro human dopamine neuron heterogeneity by scRNA-seq d Identification of dopamine subtypes with transcriptionally distinct profiles d Cell type-specific gene perturbations of cholesterol metabolism in stress response d PD GWAS associated gene expression analysis at the singlecell level

Induced pluripotent stem cells (iPSCs), and cells derived from them, have become key tools to model biological processes and disease mechanisms, particularly in cell types such as neurons that are difficult to access from living donors. Here, we present the first map of regulatory variants in an iPSC-derived cell type. To investigate genetic contributions to human sensory function, we performed 123 differentiations of iPSCs from 103 unique donors to a sensory neuronal fate, and measured gene expression, chromatin accessibility, and neuronal excitability. Compared with primary dorsal root ganglion, where sensory nerves collect near the spinal cord, gene expression was more variable across iPSC-derived neuronal cultures, particularly in genes related to differentiation and nervous system development. Single cell RNA-sequencing revealed that although the majority of cells are neuronal and express the expected marker genes, a substantial fraction have a fibroblast-like expression profile. By applying an allele-specific method we identify 3,778 quantitative trait loci influencing gene expression, 6,318 for chromatin accessibility, and 2,097 for RNA splicing at FDR 10%. A number of these overlap with common disease associations, and suggest candidate causal variants and target genes. These include known causal variants at SNCA for Parkinson's disease and TNFRSF1A for multiple sclerosis, as well as new candidates for migraine, Parkinson's disease, and schizophrenia.