Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2020 Jul 14;53(1):217-232.e5.
doi: 10.1016/j.immuni.2020.06.013.

An Integrated Multi-omic Single-Cell Atlas of Human B Cell Identity

Affiliations

An Integrated Multi-omic Single-Cell Atlas of Human B Cell Identity

David R Glass et al. Immunity. .

Abstract

B cells are capable of a wide range of effector functions including antibody secretion, antigen presentation, cytokine production, and generation of immunological memory. A consistent strategy for classifying human B cells by using surface molecules is essential to harness this functional diversity for clinical translation. We developed a highly multiplexed screen to quantify the co-expression of 351 surface molecules on millions of human B cells. We identified differentially expressed molecules and aligned their variance with isotype usage, VDJ sequence, metabolic profile, biosynthesis activity, and signaling response. Based on these analyses, we propose a classification scheme to segregate B cells from four lymphoid tissues into twelve unique subsets, including a CD45RB+CD27- early memory population, a class-switched CD39+ tonsil-resident population, and a CD19hiCD11c+ memory population that potently responds to immune activation. This classification framework and underlying datasets provide a resource for further investigations of human B cell identity and function.

Keywords: B cells; cell atlas; human; mass cytometry; multi-omic; single cell.

PubMed Disclaimer

Conflict of interest statement

Declaration of Interests The authors declare no competing interests.

Figures

None
Graphical abstract
Figure 1
Figure 1
A highly multiplexed single-cell surface screen reveals the human B cell surface proteome (A) Experimental overview (n = 2 donors). (B) Representative gating of canonical populations. (C) Representative thresholding of positivity for molecules on the screen. (D) Percent positive of total B cells (top row) and median expression by subset (bottom rows) of molecules expressed by B cells.
Figure 2
Figure 2
Differential expression analysis reveals the anergic profile of naive B cells (A) Difference in median expression for each pairwise comparison of subsets. All non-white tiles are significant (p < 0.005). (B) Volcano plots of the comparisons, colored by association with the GO term “transport.” Significantly different molecules listed in boxes are ordered by decreasing magnitude of difference of expression. (C) Median expression of transport molecules (colors). Mean of median expression of all transport molecules (black). (D) Mean of median expression of molecules associated with GO terms (color). Mean of median expression of all molecules (black). (E) Expression of six molecules more highly expressed in naive cells (p < 0.005).
Figure 3
Figure 3
CD45RB marks human memory B cells and identifies an early memory population (A) Computational workflow. (B) UMAP plots of molecules correlated with the conserved molecule in their row. Arrows indicate RB+CD27 cells. (C) UMAP plot colored by subset. (D) Biaxial plot colored by subset. Percent of CD38lo B cells in each quadrant is quantified (red text). (E) Percent of cells in each quadrant from (D). (F) Experimental workflow. (G) Mean IgH mutation frequencies. All pairwise comparisons were significantly different (p < 0.005), except where indicated. (H) Sequence diversity across a range of diversity orders. Shaded regions indicate 95% confidence interval. (I) Computational workflow. (J) Z-scores of frequencies of shared clonality.
Figure 4
Figure 4
Segregating B cells into phenotypically and isotypically distinct subsets (A) Experimental workflow (n = 3 donors). (B) Median expression by subset. (C) Percent of B cells for all subsets, colored as in (B). (D) UMAP plot generated from an equal subsampling by subset. (E) IgH isotype usage by subset. ND denotes “not determined”; IgMD denotes co-expression of IgM and IgD. (F) Euclidean distance between each B cell subset based on median expression profile. White boxes denote column minimum. (G) Subset composition by isotype as determined by canonical gating or meta-clustering. (H) Contour plots by IgH isotype. Dots and error bars indicate mean and SEM of individual donors. (I) Relative contribution of phenotype and isotype to the variance explained by linear models created to predict single-cell expression of CD79b or surface Ig.
Figure 5
Figure 5
Interrogation of B cell subset function reveals differential metabolic, biosynthesis, and immune signaling activity (A) Experimental workflow (n = 9 donors). (B) Expression of metabolic enzymes. Stars indicate significance (p < 0.005). Boxes represent interquartile range (IQR) and whiskers represent IQR +/− 1.5*IQR (C) Biaxial of biosynthesis activity. (D) Violin plots of biosynthesis activity. Bimodality indicated by arrows. (E) Violin plots of significantly differentially expressed molecules (p < 0.005). (F) Signaling diagram. (G) Median expression of Igλ B cells, grouped by phenotype or isotype. (H) Contour plots of signaling molecules. (I) Quantification of earth mover’s distance from baseline samples to stimulated samples (1 μg/mL) for pPLCγ2 and pSyk and pp38. (J) Relative contribution of phenotype and isotype to the variance explained by linear models created to predict single-cell expression of metabolic pathways, biosynthesis activity, or cell signaling.
Figure 6
Figure 6
Characterization of lymphoid tissue-specific B cell populations (A) Experimental workflow (n = 11 donors). (B) Violin plots of molecules significantly differentially expressed by at least two tissues (p < 0.005). (C) IgH isotype usage by tissue. ND denotes “not determined”; IgMD denotes co-expression of IgM and IgD. (D) Subset composition by tissue. (E) Manhattan distance between each tissue based on subset composition. (F) Biaxial of B cells colored as GC or other. (G) Biaxial of B cells colored as CD39+ tonsillar or other. (H) UMAP plot generated from an equal subsampling by tissue, then an equal subsampling by subset.

Similar articles

Cited by

  • Whole transcriptome landscape in HAPE under the stress of environment at high altitudes: new insights into the mechanisms of hypobaric hypoxia tolerance.
    Li Q, Fang F, Yang C, Yu D, Gong Q, Shen X. Li Q, et al. Front Immunol. 2024 Sep 12;15:1444666. doi: 10.3389/fimmu.2024.1444666. eCollection 2024. Front Immunol. 2024. PMID: 39328420 Free PMC article.
  • Longitudinal Multi-omic Immune Profiling Reveals Age-Related Immune Cell Dynamics in Healthy Adults.
    Gong Q, Sharma M, Kuan EL, Glass MC, Chander A, Singh M, Graybuck LT, Thomson ZJ, LaFrance CM, Zaim SR, Peng T, Okada LY, Genge PC, Henderson KE, Dornisch EM, Layton ED, Wittig PJ, Heubeck AT, Mukuka NM, Reading J, Roll CR, Hernandez V, Parthasarathy V, Stuckey TJ, Musgrove B, Swanson E, Lord C, Weiss MDA, Phalen CG, Mettey RR, Lee KJ, Johanneson JB, Kawelo EK, Garber J, Krishnan U, Smithmyer M, Wherry EJ, Vella L, Henrickson SE, Kopp MS, Savage AK, Becker LA, Meijer P, Coffey EM, Goronzy JJ, Speake C, Bumol TF, Goldrath AW, Torgerson TR, Li XJ, Skene PJ, Buckner JH, Gustafson CE. Gong Q, et al. bioRxiv [Preprint]. 2024 Sep 14:2024.09.10.612119. doi: 10.1101/2024.09.10.612119. bioRxiv. 2024. PMID: 39314416 Free PMC article. Preprint.
  • HLA-E and NKG2A Mediate Resistance to M. bovis BCG Immunotherapy in Non-Muscle-Invasive Bladder Cancer.
    Ranti D, Yu H, Wang YA, Bieber C, Strandgaard T, Salomé B, Houghton S, Kim J, Ravichandran H, Okulate I, Merritt E, Bang S, Demetriou A, Li Z, Lindskrog SV, Ruan DF, Daza J, Rai R, Hegewisch-Solloa E, Mace EM, Fernandez-Rodriguez R, Izadmehr S, Doherty G, Narasimhan A, Farkas AM, Cruz-Encarnacion P, Shroff S, Patel F, Tran M, Park SJ, Qi J, Patel M, Geanon D, Kelly G, de Real RM, Lee B, Nie K, Miake-Iye S, Angeliadis K, Radkevich E, Thin TH, Garcia-Barros M, Brown H, Martin B, Mateo A, Soto A, Sussman R, Shiwlani S, Francisco-Simon S, Beaumont KG, Hu Y, Wang YC, Wang L, Sebra RP, Smith S, Skobe M, Clancy-Thompson E, Palmer D, Hammond S, Hopkins BD, Wiklund P, Zhu J, Bravo-Cordero JJ, Brody R, Hopkins B, Chen Z, Kim-Schulze S, Dyrskjøt L, Elemento O, Tocheva A, Song WM, Bhardwaj N, Galsky MD, Sfakianos JP, Horowitz A. Ranti D, et al. bioRxiv [Preprint]. 2024 Sep 3:2024.09.02.610816. doi: 10.1101/2024.09.02.610816. bioRxiv. 2024. PMID: 39282294 Free PMC article. Preprint.
  • Deciphering bone marrow engraftment after allogeneic stem cell transplantation in humans using single-cell analyses.
    Bordenave J, Gajda D, Michonneau D, Vallet N, Chevalier M, Clappier E, Lemaire P, Mathis S, Robin M, Xhaard A, Sicre de Fontbrune F, Corneau A, Caillat-Zucman S, Peffault de Latour R, Curis E, Socié G. Bordenave J, et al. J Clin Invest. 2024 Aug 29;134(20):e180331. doi: 10.1172/JCI180331. J Clin Invest. 2024. PMID: 39207851 Free PMC article. Clinical Trial.
  • Human circulating CD24hi marginal zone B cells produce IgM targeting atherogenic antigens and confer protection from vascular disease.
    Pattarabanjird T, Nguyen AT, McSkimming C, Dinh HQ, Marshall MA, Ghosheh Y, Gulati R, Durant C, Vallejo J, Saigusa R, Drago F, Guy TV, Premo K, Taylor AM, Paul S, Kundu B, Berr S, Gonen A, Tsimikas S, Miller Y, Pillai S, Ley K, Hedrick CC, McNamara CA. Pattarabanjird T, et al. Nat Cardiovasc Res. 2023 Nov;2(11):1003-1014. doi: 10.1038/s44161-023-00356-1. Epub 2023 Oct 30. Nat Cardiovasc Res. 2023. PMID: 39196097

References

    1. Abdel-Wahab N., Shah M., Suarez-Almazor M.E. Adverse Events Associated with Immune Checkpoint Blockade in Patients with Cancer: A Systematic Review of Case Reports. PLoS ONE. 2016;11:e0160221. - PMC - PubMed
    1. Antonioli L., Pacher P., Vizi E.S., Haskó G. CD39 and CD73 in immunity and inflammation. Trends Mol. Med. 2013;19:355–367. - PMC - PubMed
    1. Becht E., McInnes L., Healy J., Dutertre C.A., Kwok I.W.H., Ng L.G., Ginhoux F., Newell E.W. Dimensionality reduction for visualizing single-cell data using UMAP. Nat. Biotechnol. 2018;37 doi: 10.1038/nbt.4314. - DOI - PubMed
    1. Bendall S.C., Davis K.L., Amir A.D., Tadmor M.D., Simonds E.F., Chen T.J., Shenfeld D.K., Nolan G.P., Pe’er D. Single-cell trajectory detection uncovers progression and regulatory coordination in human B cell development. Cell. 2014;157:714–725. - PMC - PubMed
    1. Boyd S.D., Crowe J.E. Deep sequencing and human antibody repertoire analysis. Current Opinion in Immunology. 2016;40:103–109. - PMC - PubMed

Publication types

MeSH terms