Identification of germinal center (GC) B cells is typically reliant on the use of surface activation markers that exhibit a wide range of expression. Here, we identify Ephrin-B1, a ligand for Eph-related receptor tyrosine kinases, as a specific marker of mature GC B cells. The number of Ephrin-B1+ GC B cells increases during the course of an immune response with Ephrin-B1+ GC B cells displaying elevated levels of Bcl6, S1pr2, and Aicda relative to their Ephrin-B1– counterparts. We further identified a small proportion of recently dividing, somatically mutated Ephrin-B1+ GC B cells that have begun to down-regulate Bcl6 and S1pr2 and express markers associated with memory B cells, such as CD38 and EBI2. Transcriptional analysis indicates that these cells are developmentally related to memory B cells, and likely represent a population of GC memory precursor (PreMem) B cells. GC PreMem cells display enhanced survival relative to bulk GC B cells, localize near the edge of the GC, and are predominantly found within the light zone. These findings offer insight into the significant heterogeneity that exists within the GC B cell population and provide tools to further dissect signals regulating the differentiation of GC B cells.
Germinal centers (GCs) are tightly confined clusters of cells within the follicle, in which GC B cells compete for signals necessary for their survival and continued maturation into memory B cells or plasma cells. GC B cells highly express the transcription factor Bcl6 and the G protein–coupled receptor sphingosine-1-phosphate receptor (S1PR2) that promotes their confinement within the GC (Green et al., 2011; Muppidi et al., 2014; Huang and Melnick, 2015). The GC is divided into a light zone (LZ), where GC B cells interact with antigen-bearing follicular DCs (FDCs) and follicular helper T cells, and a dark zone (DZ) in which GC B cells rapidly divide and undergo somatic hypermutation (SHM). Through regulated expression of the chemokine receptor CXCR4, GC B cells rapidly transit between these compartments, allowing for continued selection of high affinity GC B cells via competition for T cell help (Allen et al., 2007; Victora and Nussenzweig, 2012).
Memory B cells can arise from both GC-independent and -dependent pathways, with the majority of memory B cells against T cell–dependent antigens thought to originate within the GC (McHeyzer-Williams et al., 2011; Tarlinton and Good-Jacobson, 2013; Kurosaki et al., 2015). Memory B cells emerge early during the GC response and derive from lower affinity GC B cells that receive less T cell help and, accordingly, maintain higher expression of the transcription factor Bach2 (Shinnakasu et al., 2016; Weisel et al., 2016). Expression of Bach2 predisposes GC B cell to differentiate into memory B cells, whereas expression of Blimp1 promotes the development of plasma cells (Turner et al., 1994; Shinnakasu et al., 2016). Memory B cells are a heterogeneous population with distinctly functioning subsets arising within the GC at different times (Zuccarino-Catania et al., 2014; Adachi et al., 2015; Weisel et al., 2016). The exact signals regulating GC B cell differentiation into memory B cells are poorly understood.
GC B cells are typically defined through their low expression of IgD or CD38 and their positive staining for one or two surface markers. Most studies use the rat monoclonal antibody GL7, which recognizes α2,6-linked N-acetylneuraminic acid (Neu5Ac) on glycan chains, Fas (CD95), and/or peanut agglutinin (PNA), a lectin that binds exposed galactose-β(1–3)-N-acetylgalactosamine (Gal-β(1–3)-GalNac; László et al., 1993; Naito et al., 2007). However, these markers have a wide range of expression and are also expressed on activated B cells before entry into the GC, and thus do not faithfully mark only cells that have already entered the GC (Wang et al., 1996; Shinall et al., 2000; Naito et al., 2007). A marker with discrete expression that specifically labels B cells residing in the GC would be a useful tool for separating maturation occurring within the GC from events occurring before GC entry. This level of resolution is particularly important for identifying GC B cells in the process of differentiating into memory B cells and could inform efforts to more precisely decipher the signals regulating this process.
In this study, we identify Ephrin-B1, a ligand for Eph-related receptor tyrosine kinases, as a specific marker of B cells residing within the GC. Ephrin-B1 functions as a repulsive guidance cue in mice and humans and is highly expressed on GC B cells (Bush and Soriano, 2009), but we find that it is not essential for GC B cell development or positioning. Using Ephrin-B1, we found that GC B cells, defined as B220+ IgDloGL7+CD95+ cells, represent at least four cell subsets of different maturation states based on expression of markers such as Bcl6, CXCR4, S1pr2, EBI2, CD38, and CD73. Importantly, we were able to identify a recently divided, somatically mutated GC B cell population in the process of down-regulating Bcl6 and S1pr2 and up-regulating CD38 and Ebi2. Transcriptional analysis indicates that these cells are developmentally related to memory B cells and likely represent a population of GC memory precursor (PreMem) cells. Finally, we found that PreMem cells localize near the edge of the GC and are predominantly found within the LZ.
Results and discussion
Ephrin-B1 is a specific marker of mature GC B cells
Microarray gene expression analysis of CXCR4hi (DZ) and CXCR4lo (LZ) GC B cells, gated as IgDloGL7+CD95+ cells, and IgD+CD23+GL7–CD95− follicular B cells, identified Efnb1 transcripts as being highly expressed in GC B cells relative to their follicular counterparts (Fig. 1 A). Ephrin-B1 protein was highly expressed on IgDloGL7+CD95+ cells after protein antigen or sheep RBC (SRBC) immunization, but was minimally expressed by other B cell subsets in the spleen or BM, including memory B cells (Fig. 1 A, Fig. S1 A, and not depicted). Ephrin-B1 began to become up-regulated after ∼7 cell divisions in B cells responding to a T cell–dependent antigen in vivo, with its expression preceded by loss of CD38 and IgD expression and occurring well after the start of CD95 up-regulation (Fig. 1 B). Ephrin-B1 has a critical role as a repulsive guidance cue during tissue development, and mutations in the gene result in a wide spectrum of developmental abnormalities constituting craniofrontonasal syndrome in humans and related defects in mice (Bush and Soriano, 2009). Ephrin-B1 is also important in bone formation and in thymocyte development (Xing et al., 2010; Luo et al., 2011; Cejalvo et al., 2013). To test whether Ephrin-B1 may have a functional role in GC B cell development we generated mice in which Efnb1 was specifically deleted in B cells (Efnb1f/fMb1Cre/+). GC B cells from conditional KO (cKO) mice did not express Ephrin-B1 (Fig. 1 A). After SRBC immunization, cKO GCs developed that were of normal size and location and that exhibited correct organization as determined by the polarized distribution of CD35+ FDCs and CXCR4hi GC B cells (Fig. 1 C). Ephrin-B1 staining of tissue sections confirmed that Ephrin-B1 was expressed by GC B cells while also being expressed by endothelial cells of blood vessels (Fig. 1 C). To test whether Ephrin-B1 influenced the competitiveness of B cells in the GC, adoptive cell transfer experiments were performed. We generated Efnb1f/fMb1Cre/+ Hy10 and control Mb1Cre/+ Hy10 mice, which have B cells that are specific for hen egg lysozyme (HEL). Equal numbers of cKO or control HEL-specific Hy10 cells were transferred, along with WT HEL-specific Hy10 cells and OVA-specific OT-II T cells, into congenically mismatched mice 1 d before immunization with duck egg lysozyme (DEL) conjugated to OVA. We detected no defect in the development of IgDloGL7+CD95+ B cells, plasmablasts, or class switching in cKO Hy10 cells (Fig. 1 D). Furthermore, female Efnb1f/+Mb1Cre/+ mice, which have mosaic expression of Efnb1 caused by the gene’s location on the X chromosome, displayed no apparent difference in development or positioning of Efnb1-deficient or -sufficient GC B cells after SRBC immunization (Fig. 1 E). We also did not detect defects in the GC B cell response in mediastinal LNs or spleen after influenza HKx31 infection (two experiments, cKO versus WT and mixed BM chimeras, with at least three mice per group), in the spleen after acute lymphocytic choriomeningitis virus strain Armstrong (LCMV Arm) infection (two mixed BM chimera experiments with three mice per group), in skin draining LNs after subcutaneous immunization with protein antigen in complete Freund’s adjuvant, or in the chronic GCs occurring in Peyer’s patches in response to commensal-derived antigens (three experiments, cKO versus WT and mixed BM chimeras, at least four mice per group; unpublished data). These observations suggest that Efnb1 is not essential for the development or organization of GCs.
Investigating the kinetics of Ephrin-B1 expression after LCMV Arm infection revealed that there was a marked increase in the proportion and numbers of IgDloGL7+CD95+ B cells that were Efnb1+ during the course of the GC response (Fig. 2 A). Conversely, the number of Efnb1– GL7+CD95+IgDlo B cells decreased over time (Fig. 2 A). Similar findings were obtained when using PNA in place of GL7 (unpublished data) and when analyzing Ephrin-B1 expression kinetics after SRBC immunization (Fig. 2 B). To examine whether Eprhin-B1 may mark distinct GC B cell populations, we determined the expression of markers associated with GC B cell fate and function in these populations. Efnb1+ cells displayed reduced expression of CD38 and increased expression of CD73 relative to their Efnb1– counterparts (Fig. 2, C and D, red, white and blue symbols in D correspond to days 7, 11, and 15, respectively). CD73 expression increases within GC B cells over time, whereas CD38 expression decreases, suggesting that Ephrin-B1 may distinguish more temporally mature GC B cells (Conter et al., 2014). Consistent with this notion, Efnb1– cells expressed intermediate levels of Bcl6 and heterogeneous amounts of S1pr2 (as detected using S1pr2Venus/+ reporter mice) and activation-induced cytidine deaminase (AID; as revealed using AID-GFP reporter mice), indicating that they were not fully confined to the GC or uniformly undergoing SHM (Fig. 2, C and D; Crouch et al., 2007; Moriyama et al., 2014). In contrast, Efnb1+ cells expressed high levels of Bcl6, S1pr2, and Aicda (Fig. 2, C and D). Collectively, these data indicate that Ephrin-B1 is a specific marker of mature GC-resident B cells.
Transitional GC B cells can be identified using Ephrin-B1 and
In the course of these studies, we observed that whereas the great majority of Efnb1+ cells were S1pr2hi, a small fraction of the cells (1–2%) had low S1pr2 expression (Fig. 2, C and D). Low S1pr2 expression in both Efnb1– and Efnb1+ IgDloGL7+CD95+ cells positively correlated with high CD38 expression, suggesting that these cells might represent transitional populations in the process of entering or exiting the GC, as CD38 is highly expressed on both follicular and memory B cells (Fig. 3 A). Both Efnb1– and Efnb1+ S1pr2lo cells (populations 1 and 4, respectively) expressed low levels of Bcl6 and CXCR4, consistent with these cells being in a distinct state from the bulk GC B cell population (Fig. 3 A). Efnb1– and Efnb1+ S1pr2hi cells (populations 2 and 3, respectively) conversely expressed high levels of Bcl6 and CXCR4, and had a similar fraction of DZ and LZ B cells (Fig. 3 A; detailed further in Fig. 5 A). S1pr2lo Efnb1– (Pop 1) cells had low expression of CD73, suggesting that they may largely be precursor GC cells (Fig. 3 A). Importantly, S1pr2lo Efnb1+ (Pop 4) cells expressed similar levels of CD73 as their S1pr2hi (Pop 3) counterparts, raising the possibility that they may represent a cell population in the process of exiting the GC (Fig. 3 A). The percentage of IgDloGL7+CD95+ B cells comprising Pop 1 and Pop 4 peaked at day 7 and 11, respectively, before declining, as Pop 3 increasingly dominated the GC response (Fig. 3 B). Pop 2 remained a small but stable fraction of the IgDloGL7+CD95+ B cells, even at late GC time points (Fig. 3 B). Together, these results suggest a model in which newly activated IgDloGL7+CD95+ B cells up-regulate S1pr2 before Ephrin-B1 as they enter the GC response, with Efnb1+ cells losing S1pr2 expression as they differentiate from the bulk GC state.
As plasma cells and memory B cells can both derive from the GC, we investigated whether S1pr2lo Efnb1+ (Pop 4) cells included plasma cell precursors by examining expression of the plasma cell marker CD138. Although the vast majority of CD138-expressing cells were B220lo/int (not depicted), we identified a small population of B220+ IgDloGL7+CD95+ B cells expressing CD138, likely representing plasma cell precursors (prePCs; Fig. 3 C; Fooksman et al., 2010). PrePCs were largely Efnb1– and found in the S1pr2lo Efnb1– (Pop 1) population, with only ∼2% of S1pr2lo Efnb1+ (Pop 4) cells expressing CD138 (Fig. 3 C). Similar results were found using Blimp1GFP mice to identify prePCs (unpublished data; three experiments, four mice per experiment). This result suggested that, at this time point of LCMV Arm infection, the majority of prePCs differentiate from activated rather than GC-resident B cells, with GC-derived prePCs perhaps arising at later GC time points (Weisel et al., 2016). Alternatively, GC-derived prePCs could lose expression of Ephrin-B1 more rapidly than those of GL7 or CD95. These findings also indicated that by defining transitional GC populations based on high expression of CD38 and low expression of CXCR4 (Fig. 3 A), we could evaluate these populations without significant contamination from prePCs, which express low CD38 levels and maintain elevated CXCR4 expression (Bhattacharya et al., 2007). Transitional populations defined in this way strongly enriched for S1pr2lo cells, allowing these populations to be evaluated in nonreporter mice (Fig. 3 D and Fig. S1, B and C). All four populations identified using CD38 and CXCR4 displayed a similar dependence on CD40-mediated signals for their survival, as treatment of mice with αCD40L antibody (MR1) from day 9–12 after LCMV Arm infection resulted in roughly equivalent decreases in the number of these cells (Fig. 3 E). This treatment only resulted in about a twofold decrease in memory B cells (CD73+CD95+CD38+IgDloGL7– B cells; Fig. S1 B), likely caused by a significant fraction of these cells exiting the GC before treatment (Fig. 3 E; Anderson et al., 2007; Dogan et al., 2009).
Efnb1+ S1pr2lo cells represent a post-GC B cell population
We next evaluated the extent of cell division occurring in these populations 24 h before analysis. To do this, we used a tet-off-H2b-GFP mouse that allows proliferation to be assessed over multiple cell divisions (Foudi et al., 2009; Gitlin et al., 2014; Bannard et al., 2016). Treatment of transgenic mice with doxycycline results in the H2b-GFP gene being turned off and dividing cells diluting their expression of GFP as histone segregation occurs during cell division. Within the IgDloGL7+CD95+ compartment, Efnb1–CD38+CXCR4lo (Pop 1) cells displayed minimal GFP dilution further indicating that these cells are precursor GC cells that have not yet entered a rapidly dividing state (Fig. 3 F). Efnb1–CD38– (Pop 2) cells displayed some dilution of GFP, indicative of a small amount of division, but did not have a sizeable population that had fully lost GFP expression (Fig. 3 F). Efnb1+CD38– (Pop 3) cells, in contrast, displayed robust loss of GFP, consistent with this population representing mature GC B cells (Fig. 3 F). Efnb1+CD38+CXCR4lo (Pop 4) cells demonstrated a similar dilution of GFP (Fig. 3 F). These data indicate that Pop 4 is derived from recently dividing GC B cells and could represent a population in the process of differentiating away from the GC B cell state.
We further tested the proximity of these populations to the GC state by examining the extent of SHM. SHM was assessed by sequencing of the IgH JH4-intronic enhancer downstream of the rearranged VJ558DJH4 element in DNA from the populations (Park et al., 2009; Bannard et al., 2013). The mutation frequency in this intronic region provides a measurement of AID activity (Jolly et al., 1997). As expected, follicular B cells were largely germline, with some sequences having a small number of mutations, likely caused by errors introduced during the nested PCR reactions (Fig. 3 G; Bannard et al., 2013). Memory B cells were split between cells with low and high levels of mutation consistent with some of these cells emerging from the GC during early time points after infection, before significant mutations have occurred (Fig. 3 G). PrePCs had predominantly low to intermediate numbers of mutations (Fig. 3 G). Efnb1– S1pr2hi (Pop 2) cells displayed low levels of mutations, suggesting that these cells have either only recently adopted a GC state or are intrinsically impaired in their ability to fully participate in the GC reaction and undergo rapid division and SHM. Both Efnb1+ S1pr2hi (Pop 3) and Efnb1+ S1pr2lo (Pop 4) cells had a high fraction of mutated cells (Fig. 3 G). The mutational frequency in Efnb1+ S1pr2lo (Pop 4) cells was slightly less than in that of the Efnb+ S1pr2hi (Pop 3) cells, perhaps indicative of AID activity being lost in this population as it transitions from the GC state. This slightly reduced mutational frequency could also be caused by GC B cells of lower affinity being more prone to differentiate into memory B cells (Shinnakasu et al., 2016).
Efnb1+ S1pr2lo cells are developmentally related to memory B cells
To further probe the developmental relationship between Efnb1+ S1pr2hi (Pop 3), Efnb1+ S1pr2lo (Pop 4), and memory B cells, we performed RNA-seq analysis at day 11 after LCMV Arm infection. Principal component analysis (PCA) indicated that Efnb1+ S1pr2hi (Pop 3) and Efnb1+ S1pr2lo (Pop 4) cells had vastly different transcriptional profiles, despite both expressing surface markers associated with GC B cells (Fig. 4 A). Instead, Efnb1+ S1pr2lo (Pop 4) cells clustered tightly with memory B cells and displayed striking similarities in their gene expression profile based on an unbiased analysis of genes with the greatest variance between groups (Fig. 4, A and B; and Fig. S2 A). Published RNA-seq data from splenic follicular B cells was included in the PCA as a reference (Shi et al., 2015). In total, there were 497 differentially expressed genes (DEGs) between Efnb1+ S1pr2hi (Pop 3) and Efnb1+ S1pr2lo (Pop 4) cells, and only 38 DEGs between Efnb1+ S1pr2lo (Pop 4) cells and memory B cells, with DEGs being defined as those with a Padj < 0.1, base count > 100, and log2foldchange > 1.5 (Fig. S2, B and C). Among the DEGs that were down-regulated in Efnb1+ S1pr2lo (Pop 4) cells and memory B cells were those associated with the GC B cell state, including Bcl6, Aicda, Bach2, Cxcr4, and Mki67 (Fig. 4 B). Although Bach2 has been reported to have higher expression on GC B cells predisposed to adopt a memory fate, its low expression on both Efnb1+ S1pr2lo (Pop 4) cells and memory B cells suggests that Bach2hi precursor cells are still in a GC state and that further development must occur within the GC before they can complete their differentiation into memory B cells (Shinnakasu et al., 2016). Up-regulated DEGs include many commonly associated with cell positioning (Ccr7, S1pr1, Ccr6, Gpr183, Cxcr3, and Sell), activation (Itgam, Cd44, and Tlr7), cytokine signaling (Il9r, Il27ra, Il15ra, and Il6ra), transcriptional regulation (Zeb2), and cell survival (Bcl2; Fig. 4 B). In several cases, these DEGs were validated by protein expression analysis (Fig. S2 D). Efnb1+CD38+CXCR4lo (Pop 4) cells displayed similar levels of death and apoptosis to memory B cells after in vitro culture, whereas Efnb1–CD38– (Pop 2) and Efnb1+CD38– (Pop 3) cells demonstrated poor survival (Fig. 4 C). Collectively, these data indicate that Efnb1+ S1pr2lo (Pop 4) cells are transcriptionally and functionally similar to memory B cells and likely represent a population of GC memory precursor (PreMem) B cells.
Ingenuity Pathway Analysis (IPA; Ingenuity Systems) was also used to explore upstream regulators of the gene signature. IPA found that the Efnb1+ S1pr2hi (Pop 3) to Efnb1+ S1pr2lo (Pop 4) transition was driven primarily by a loss in Bcl6 activity and an increase in genes associated with inflammatory exposure, possibly resulting from the loss of Bcl6-mediated repression of inflammatory signaling (Fig. S3 A; Dent et al., 1997; Cui et al., 2011). The Efnb1+ S1pr2lo (Pop 4) to memory B cell transition was marked by a loss in genes associated with cell cycle progression and an increase in those associated with a more quiescent phenotype, consistent with the notion that Efnb1+ S1pr2lo cells are still in the process of completing their differentiation into long-lived memory B cells (Fig. S3 B).
Efnb1+ S1pr2lo cells localize near the edge of the GC LZ
Finally, we sought to determine where PreMem B cells localize within the GC. Phenotypic analysis indicated that Efnb1+ S1pr2lo (Pop 4) cells had an increased propensity to adopt a LZ GC B cell phenotype and therefore might preferentially localize to this zone (Fig. 5 A). To positively identify these transitioning cells, we used EBI2GFP/+ mice, as the RNA-seq data indicated that Gpr183 (the gene encoding EBI2) was highly expressed within Efnb1+ S1pr2lo (Pop 4) cells. Indeed, use of EBI2-GFP along with Ephrin-B1 allowed for clear detection of four populations of GL7+CD95+IgDlo B cells with equivalent expression of CD38, CD73, Bcl6, and CXCR4 to those identified using S1pr2 and Ephrin-B1 (Fig. 5 B and not depicted). Identification of transitional cells using CD38 and CXCR4 also enriched for cells expressing EBI2-GFP (Fig. S1, B and C). Analysis of splenic sections from EBI2GFP/+ mice after LCMV Arm infection identified a small population of Efnb1+GL7+IgDlo cells within the GC that displayed intracellular expression of GFP, and these cells localized near the edge of the GC (Fig. 5 C). They were found within the LZ in 22 of 25 cases where they were identified, and GC polarization could be determined on serial sections (Fig. 5 D). These data are consistent with the notion that memory precursor cells predominantly arise in the LZ and exit the GC from this compartment.
In summary, we find that Ephrin-B1 is a specific marker of mature GC-resident B cells and that it facilitates the identification of transitional populations of GC B cells. Despite the high expression, Ephrin-B1 was not required for B cell participation in the GC response. However, because Eph receptors are expressed on various stromal and hematopoietic cells, including T cells (Alfaro et al., 2007; Maddigan et al., 2011), it remains possible that Ephrin-B1 will influence GC B cell selection or differentiation events under some conditions. We demonstrate that Efnb1+ S1pr2lo cells are somatically mutated, derived from recently dividing GC B cells, display enhanced survival, and are transcriptionally similar to memory B cells and thus likely represent a population of GC memory precursor B cells. However, we do not exclude the possibility that the PreMem B cell population is heterogeneous in its origins and properties. Future studies with lineage reporter mice and single cell RNA sequencing will help further define the pathway(s) leading to memory cell commitment in the GC. Based on the striking overlap in gene signature between the PreMem B cells and memory B cells, we speculate that memory B cell differentiation occurs to a large extent within the GC. As these cells lose GC confinement signals and become responsive to migration cues found outside the GC, they likely will exit the GC and migrate to niches, where they can complete their development into long-lived memory B cells. Several cell surface receptors, including Il9r, Il6ra, Il15ra, and Il27ra, were up-regulated on PreMem B cells and could be involved in promoting memory cell differentiation. Our work suggests that the GL7, PNA, CD95, and IgD marker combination is insufficient to precisely identify GC B cells. Incorporation of Ephrin-B1 in the marker combination highly enriches for GC B cells that display elevated expression of S1pr2, Bcl6, and Aicda, and which are undergoing rapid cell division. Further addition of markers such as CD38 and CXCR4 allow for segregation of GC precursor and PreMem B cells, the latter being identified as Efnb1+CD38hiCXCR4lo cells. The use of Ephrin-B1 together with these established markers should prove beneficial in studies aimed at gaining a more precise understanding of memory B cell differentiation.
Materials and methods
Adult C57BL/6 CD45.1+ (stock number 564) mice at least 6 wk of age were purchased from the National Cancer Institute (NCI) or Charles River. Mb1Cre mice were provided by M. Reth (Max Planck Institute of Immunobiology and Epigenetics, Freiburg, Germany; Hobeika et al., 2006). Efnb1f/f mice were provided by A. Soriano (Mount Sinai School of Medicine, New York, NY; Davy et al., 2004). S1pr2Venus/+ were generated as previously described (Moriyama et al., 2014). Tet-off H2b-GFP mice were generated by crossing TetOp-H2b-GFP, ROSA:LNL:tTA, and Mb1Cre mice (Hobeika et al., 2006; Wang et al., 2008; Foudi et al., 2009). EBI2GFP/+ (containing a GFP reporter in place of the Gpr183 coding exon), HEL-specific Hy10, OVA-specific OTII TCR-transgenic, and MD4-Ig transgenic mice were from an internal colony and have been previously described (Pereira et al., 2009; Yi et al., 2012; Yi and Cyster, 2013). AID-GFP mice were purchased from Jackson ImmunoResearch Laboratories. Mice were housed in a specific pathogen–free environment in the Laboratory Animal Research Center at the University of California (San Francisco [UCSF], CA), and all animal procedures were approved by the UCSF Institutional Animal Care and Use Committee.
Adoptive transfer, immunization, infections, and treatments
For experiments involving transfer of Hy10 or OT-II cells, 105 cells of each population were adoptively transferred into mice 1 d before immunization. To visualize proliferation of antigen-specific B cells, 30–50 × 106 MD4 cells were CFSE-labeled according to the manufacturer’s instructions (Invitrogen) and transferred 1 d before immunization. Mice were immunized with 2 × 108 SRBC (Colorado serum) by i.p. injection or 50 µg of duck egg lysozyme conjugated to OVA (DEL-OVA) in Sigma or Ribi Adjuvant System (Sigma-Aldrich). Mice were infected with 2 × 105 plaque-forming units of LCMV Armstrong administered i.p. For anti-CD40L treatment, mice were treated with neutralizing antibody by i.p. injection of 250 µg anti–mouse CD40L (clone MR1) daily for 4 d before sacrifice. For tet-off experiments, mice were treated with 1.6 mg doxycycline (Sigma-Aldrich) administered i.p. in saline, and then maintained by including doxycycline (2 mg/ml) and sucrose (2%) in their drinking water until analysis.
Antibodies for flow cytometry and microscopy staining
Spleens were mashed through a 70-µm cell strainer, and RBCs were lysed with RBC lysing buffer. Lymphocytes were then washed and counted. The following antibodies were used for flow cytometry and microscopy staining: Phycoerythrin (PE) anti-CD86 (105008), phycoerythrin–indotricarbocyanine (PE-Cy7) anti-CD38 (102718), Brilliant Violet 605 (BV605) anti-CD45.1 (110738), allophycocyanin (APC) anti-GL7 (144606), Pacific Blue (PacBlue) anti-GL7 (144614), peridinin chlorophyll protein Cy5.5 (PerCpCy5.5) anti-CD73 (127214), PacBlue anti-IgD (405712), PerCpCy5.5 anti-IgD (405710), APC anti-CD80 (104718), PE anti-CD11b (101208), PE anti-CD44 (103008), PE anti-IgD (405705), PE anti-CD62L (104418), APC Cy7 anti-B220 (103224), FITC anti-GL7 (144604), PE anti-CD138 (142504; all from BioLegend); APC anti-CXCR4 (558644), PE anti-CD95 (554258), PE Cy7 anti-CD95 (557653), APC anti-Bcl6 (561525), APC anti-CCR6 (557976), Biotin anti-CD35 (553816), Biotin anti-CXCR4 (551968), APC anti-TCRβ (17–5961-82), BV605 Streptavidin (563260), APC Cy7 anti-CD19 (115530), PE Texas Red Streptavidin (551487; all from BD); PerCPCy5.5 anti-CD45.2 (65–0454-U100; Tonbo Biosciences); Biotin donkey anti–mouse polyclonal Ephrin-B1 (BAF473; R&D Systems); Alexa Fluor 488 Rabbit polyclonal anti-GFP (A-21311; Invitrogen/Life Technologies); Goat anti–mouse IgD (goat polyclonal GAM/IGD(FC)/7S), APC anti-CD23 (CL8910APC; Cedarlane Labs); aminomethylcoumarin-donkey anti–goat (705–156-147), Cy3 Streptavidin (016–160-084; Jackson ImmunoResearch Laboratories). Flow cytometry data were acquired on a LSRII with FACSDiva software (BD) and were analyzed with FlowJo software (Tree Star).
Immunohistochemistry (IHC) and immunofluorescence (IF) microscopy
For IHC, 7-µm cryosections were acetone fixed and stained as previously described (Allen et al., 2007) For IF microscopy, 7-µm cryosections were prepared as previously described (Reboldi et al., 2016). For staining of mouse Ephrin-B1 for IF images, a tyramide kit was used according to the manufacturer’s instructions (TSA Biotin System; Perkin Elmer). Images were captured with a Zeiss AxioObserver Z1 inverted microscope. Identification of the center of the GC and distance from the center to specific cells was performed using Imaris (Bitplane). If the GC is assumed to be a circle with GC B cells evenly distributed throughout, the mean distance from the edge of the GC to an individual cell divided by the total distance from the cell to the GC center, and the cell to the GC edge is predicted to be ∼0.375. This is a result of more cells being located in the outer part of the GC, which has greater area than the inner GC.
Analysis of somatic mutations in JH558 intron
20,000–40,000 cells for each population were FACS sorted, and DNA was isolated using a QIAamp DNA Micro kit (QIAGEN). DNA was eluted in 20 µl and used as a template for nested PCR (Park et al., 2009; Bannard et al., 2013). New primers and reagents were added directly to the first PCR product for the secondary reaction. Primers are specific for the JH558 family members, and only PCR products of ∼700 bp were excised and purified from 1.2% agarose gels using a QIAquick Gel Extraction kit (QIAGEN). The purified gel product was then cloned into the pCR4 Bunt-TOPO vector according to the manufacturer’s instructions (Invitrogen). Colonies were then submitted for preparation and sequencing (TACgen), and sequences aligned to the germline JH558 intronic sequence using standard nucleotide BLAST. Primers are as follows: Nested Forward 1, 5′-AGCCTGACATCTGAGGAC-3′; nested reverse 1, 5′-TCTGATCGGCCATCTTGACTC-3′; nested forward 2, 5′-CATCTGAGGACTCTGCGGTCT-3′; nested reverse 2, 5′-CTGTGTTCCTTTGAAAGCTGG-3′.
RNA-seq library preparation and data analysis
Total RNA was purified from FACS sorted cells using the RNeasy Mini kit (QIAGEN). RNA quality was assessed with an Agilent 2100 Bioanalyzer (RNA integrity number > 9 for all samples). Barcoded sequencing libraries were generated with 100 ng of RNA with an Ovation RNA-seq System V2 (Nugen), KAPA Hyper Prep kit for Illumina (KAPA Biosystems), and NEXTflex DNA barcodes (Bioo Scientific). Single-end sequencing was performed on an HiSeq 2500 (Illumina; UCSF Center for Advanced Technology), and sequences reported as FASTQ files, which were aligned to the mm10 genome with STAR (Spliced Transcript Alignment to a Reference). Mappable reads were counted with HTseq and imported into RStudio software for analysis of differential expression with DESeq2 software. For the generation of the heat map, select genes with a difference in expression (log2) of 1.5-fold, p-adjusted value < 0.1, and mean count across all groups >100 were chosen for visualization with values scaled by row.
The pathway analysis was performed using the Ingenuity Pathway Analysis software (QIAGEN). The Upstream analysis was performed to identify upstream regulators across the Efnb1+S1pr2hi (Pop3) and Efnb1+S1pr2lo (Pop 4) groups as well as the Efnb1+S1pr2lo (Pop 4) and memory B cell groups. Upstream regulator predictions were made by the z-score algorithm. All upstream regulators shown have P < 10−9.
The RNA-seq data reported in this paper are available at the Gene Expression Omnibus under accession no.GSE89897.
Results represent the mean ± SEM unless indicated otherwise. Statistical significance was determined by the unpaired Student’s t test. Statistical analyses were performed using Prism GraphPad software v5.0. (*, P < 0.05; **, P < 0.01; ***, P < 0.001).
Online supplemental material
Fig. S1 shows an example flow cytometric gating scheme to identify GC B cell subsets and memory B cells. Fig. S2 shows a heat map and volcano plots derived from the RNA-seq analysis, along with protein validation of select DEGs. Fig. S3 shows an IPA upstream regulator analysis of the RNA-seq data.
We thank Phil Soriano for Efnb1f/f mice, M. Matloubian for help in providing LCMV Arm virus, and Y. Xu and J. An for expert technical assistance.
This work was supported by grants from the National Institutes of Health (R01AI045073 [J.G. Cyster], R01AI040098 [J.G. Cyster], and T32AI07019 [B.J. Laidlaw]). B.J. Laidlaw is a recipient of a postdoctoral fellowship from the Damon Runyon Cancer Institute. J.G. Cyster is an investigator of the Howard Hughes Medical Institute.
The authors declare no competing financial interests.
Author contributions: B.J. Laidlaw and J.G. Cyster conceived and designed the experiment. B.J. Laidlaw, T.H. Schmidt, J.A. Green and C.D.C. Allen performed the experiments. B.J. Laidlaw and J.G. Cyster analyzed the data. T. Okada provided tools for the experiments. B.J. Laidlaw and J.G. Cyster wrote the manuscript.
activation-induced cytidine deaminase
differentially expressed gene
hen egg lysozyme
lymphocytic choriomeningitis virus