Reference-free multiplexed single-cell sequencing identifies genetic modifiers of the human immune response

Multiplexed single-cell sequencing (mux-seq) using single-nucleotide polymorphisms (SNPs) has emerged as an efficient approach to perform expression quantitative trait loci (eQTL) studies that map interactions between genetic variants and cell types, cell states, or experimental perturbations. Here we introduce the clue framework, a novel approach to encode mux-seq experiments that eliminates the need for reference genotypes and experimental barcoding. The clue framework is made possible by the development of freemuxlet, an algorithm that clusters cells based on SNPs called from single-cell RNA-seq or ATAC-seq data. To demonstrate the feasibility of clue, we profiled the surface protein and RNA abundances of peripheral blood mononuclear cells from 64 individuals, stimulated with 5 distinct extracellular stimuli — all within a single day. Our analysis of the demultiplexed data identified rare immune cell types and cell type-specific responses to interferon and toll-like receptor stimulation. Furthermore, by integrating genotyping data, we mapped response eQTLs specific to certain cell types. These findings showcase the potential and scalability of the clue framework for reference-free multiplexed single-cell sequencing studies.

[1]  Chun Jimmie Ye,et al.  Cell-subtype specific effects of genetic variation in the aging and Alzheimer cortex , 2022, bioRxiv.

[2]  M. Foroutan,et al.  A new highly-specific Natural Killer cell-specific gene signature predicting recurrence in colorectal cancer patients , 2022, bioRxiv.

[3]  Chun Jimmie Ye,et al.  Single-cell eQTL mapping identifies cell type–specific genetic control of autoimmune disease , 2022, Science.

[4]  Michael J. T. Stubbington,et al.  Integrating T cell receptor sequences and transcriptional profiles by clonotype neighbor graph analysis (CoNGA) , 2021, Nature Biotechnology.

[5]  R. V. van Lier,et al.  Hobit identifies tissue-resident memory T cell precursors that are regulated by Eomes , 2021, Science Immunology.

[6]  Xuegong Zhang,et al.  Discovering single-cell eQTLs from scRNA-seq data only , 2021, bioRxiv.

[7]  T. Chan,et al.  High FRMD3 expression is prognostic for worse survival in rectal cancer patients treated with CCRT , 2021, International Journal of Clinical Oncology.

[8]  S. Raychaudhuri,et al.  Multimodally profiling memory T cells from a tuberculosis cohort identifies cell state associations with demographics, environment, and disease , 2021, Nature Immunology.

[9]  M. Frye,et al.  ERICH3: vesicular association and antidepressant treatment response , 2020, Molecular Psychiatry.

[10]  Casey S. Greene,et al.  Genetic demultiplexing of pooled single-cell RNA-sequencing samples in cancer facilitates effective experimental design , 2020, bioRxiv.

[11]  Daniel J. Gaffney,et al.  Souporcell: robust clustering of single-cell RNA-seq data by genotype without reference genotypes , 2020, Nature Methods.

[12]  Lindsey A. Muir,et al.  Cholesterol 25-hydroxylase (CH25H) as a promoter of adipose tissue inflammation in obesity and diabetes , 2020, Molecular metabolism.

[13]  C. Stancu,et al.  Ninjurin-1 upregulated by TNFα receptor 1 stimulates monocyte adhesion to human TNFα-activated endothelial cells; benefic effects of amlodipine. , 2020, Life sciences.

[14]  O. Blacque,et al.  ERICH3 in Primary Cilia Regulates Cilium Formation and the Localisations of Ciliary Transport and Sonic Hedgehog Signaling Proteins , 2019, Scientific Reports.

[15]  Melanie Bahlo,et al.  Inhibition of Upf2-Dependent Nonsense-Mediated Decay Leads to Behavioral and Neurophysiological Abnormalities by Activating the Immune Response , 2019, Neuron.

[16]  F. Vély,et al.  Imbalance of Circulating Innate Lymphoid Cell Subpopulations in Patients With Septic Shock , 2019, Front. Immunol..

[17]  J. Schultze,et al.  Human Monocyte Subsets and Phenotypes in Major Chronic Inflammatory Diseases , 2019, Front. Immunol..

[18]  M. Betts,et al.  Human CD4+CD103+ cutaneous resident memory T cells are found in the circulation of healthy individuals , 2019, Science Immunology.

[19]  Richa Bharti,et al.  Contribution of LTi and TH17 cells to B cell aggregate formation in the central nervous system in a mouse model of multiple sclerosis , 2019, Journal of Neuroinflammation.

[20]  L. Coin,et al.  Genotype-free demultiplexing of pooled single-cell RNA-seq , 2019, Genome Biology.

[21]  J. Niu,et al.  Hepatitis B Virus–Induced Imbalance of Inflammatory and Antiviral Signaling by Differential Phosphorylation of STAT1 in Human Monocytes , 2019, The Journal of Immunology.

[22]  Mingyao Liu,et al.  IFN-stimulated P2Y13 protects mice from viral infection by suppressing the cAMP/EPAC1 signaling pathway , 2018, Journal of molecular cell biology.

[23]  R. Locksley,et al.  Innate Lymphoid Cells: 10 Years On , 2018, Cell.

[24]  Evan Z. Macosko,et al.  Heritability enrichment of specifically expressed genes identifies disease-relevant tissues and cell types , 2018, Nature Genetics.

[25]  M. Swertz,et al.  Single-cell RNA sequencing identifies celltype-specific cis-eQTLs and co-expression QTLs , 2018, Nature Genetics.

[26]  Chun Jimmie Ye,et al.  Multiplexed droplet single-cell RNA-sequencing using natural genetic variation , 2017, Nature Biotechnology.

[27]  J. Derisi,et al.  IFN‐stimulated Gene Expression, Type 2 Inflammation, and Endoplasmic Reticulum Stress in Asthma , 2017, American journal of respiratory and critical care medicine.

[28]  Laura J. Scott,et al.  Genetic regulatory signatures underlying islet gene expression and type 2 diabetes , 2017, Proceedings of the National Academy of Sciences.

[29]  S. Sumitomo,et al.  Emerging roles of Egr2 and Egr3 in the control of systemic autoimmunity. , 2016, Rheumatology.

[30]  L. Steinman,et al.  IFN-β Treatment Requires B Cells for Efficacy in Neuroautoimmunity , 2015, The Journal of Immunology.

[31]  A. Skol,et al.  Enrichment of inflammatory bowel disease and colorectal cancer risk variants in colon expression quantitative trait loci , 2015, BMC Genomics.

[32]  Lijing Yao,et al.  Global loss of DNA methylation uncovers intronic enhancers in genes showing expression changes , 2014, Genome Biology.

[33]  Chun Jimmie Ye,et al.  Intersection of population variation and autoimmunity genetics in human T cell activation , 2014, Science.

[34]  Jingyuan Fu,et al.  Calling genotypes from public RNA-sequencing data enables identification of genetic variants that affect gene-expression levels , 2014, Genome Medicine.

[35]  Daphne Koller,et al.  Polarization of the Effects of Autoimmune and Neurodegenerative Risk Alleles in Leukocytes , 2014, Science.

[36]  Shing-Jyh Chang,et al.  Differential expression of distinct surface markers in early endothelial progenitor cells and monocyte-derived macrophages. , 2013, Gene expression.

[37]  G. Abecasis,et al.  Detecting and estimating contamination of human DNA samples in sequencing and array-based genotype data. , 2012, American journal of human genetics.

[38]  K. Hao,et al.  Bayesian method to predict individual SNP genotypes from gene expression data , 2012, Nature Genetics.

[39]  C. Carlson,et al.  Principles for the post-GWAS functional characterization of cancer risk loci , 2011, Nature Genetics.

[40]  S. Sánchez-Ramón,et al.  The chemokine CXCL12 regulates monocyte-macrophage differentiation and RUNX3 expression. , 2011, Blood.

[41]  Snehit Prabhu,et al.  Overlapping Pools for High Throughput Targeted Resequencing , 2009, RECOMB.

[42]  L. Staudt,et al.  Characterization of early stages of human B cell development by gene expression profiling , 2009, The Journal of Immunology.

[43]  Yi Li,et al.  Type I interferon modulates monocyte recruitment and maturation in chronic inflammation. , 2009, The American journal of pathology.

[44]  S. Swerdlow,et al.  Benign CD10-Positive T Cells in Reactive Lymphoid Proliferations and B-Cell Lymphomas , 2003, Modern Pathology.

[45]  P Green,et al.  Base-calling of automated sequencer traces using phred. II. Error probabilities. , 1998, Genome research.

[46]  P. Green,et al.  Base-calling of automated sequencer traces using phred. I. Accuracy assessment. , 1998, Genome research.

[47]  Margaret Comerford Freda,et al.  10 Years. , 2008, MCN. The American journal of maternal child nursing.

[48]  W. Tourtellotte,et al.  Early Growth Response Transcriptional Regulators Are Dispensable for Macrophage Differentiation , 2007 .