Cyclosporine a directly affects human and mouse b cell migration in vitro by disrupting a hIF-1 αdependent, o2 sensing, molecular switch

Background Hypoxia is a potent molecular signal for cellular metabolism, mitochondrial function, and migration. Conditions of low oxygen tension trigger regulatory cascades mediated via the highly conserved HIF-1 α post-translational modification system. In the adaptive immune response, B cells (Bc) are activated and differentiate under hypoxic conditions within lymph node germinal centers, and subsequently migrate to other compartments. During migration, they traverse through changing oxygen levels, ranging from 1-5% in the lymph node to 5-13% in the peripheral blood. Interestingly, the calcineurin inhibitor cyclosporine A is known to stimulate prolyl hydroxylase activity, resulting in HIF-1 α destabilization and may alter Bc responses directly. Over 60% of patients taking calcineurin immunosuppressant medications have hypo-gammaglobulinemia and poor vaccine responses, putting them at high risk of infection with significantly increased morbidity and mortality. Results We demonstrate that O 2 tension is a previously unrecognized Bc regulatory switch, altering CXCR4 and CXCR5 chemokine receptor signaling in activated Bc through HIF-1 α expression, and controlling critical aspects of Bc migration. Our data demonstrate that calcineurin inhibition hinders this O 2 regulatory switch in primary human Bc. Conclusion This previously unrecognized effect of calcineurin inhibition directly on human Bc has significant and direct clinical implications.

particular cellular migration, are sparsely described in the literature [5].
Several pieces of evidence suggest that hypoxia may modulate post-germinal center (GC) B cell migration, and that CNI could interfere with this process. Low tissue oxygen (O 2 ) tension (≤1% O 2 ; hypoxia) has emerged as an important immune modulating signal [6][7][8][9][10][11]. Intriguingly, low O 2 tension occurs in secondary lymphoid organs, particularly in draining lymph nodes (LN) and bone marrow (BM) [12,13]. Lymphocytes within these organs express hypoxia induced (transcription) factor-1α (HIF-1α) at low O 2 tensions (1% -4%) [14][15][16]. It has recently been shown that GC hypoxia itself is important for efficient B cell class switching and antibody production in mouse models [17,18], directly linking hypoxia to B cell function. Despite these results, the mechanism by which B cell function is directly affected by hypoxia during migration to regions of secondary lymphoid tissues with low oxygen O 2 tension, remains undefined [12,13].
Although current data support the premise that O 2 tension modulates transcriptome and molecular signaling events in human T [14,15,[19][20][21] and lymphoma cells [22][23][24], the effects of O 2 tension has on B cell responses are largely unknown. We have previously demonstrated that HIF-1α transcripts are upregulated in both human differentiating B cells in vitro and plasma cells migrating in vivo through peripheral blood to bone marrow post-vaccination [25,26]. Coordinated migration of B cells between GC, peripheral blood (PB), spleen and BM is critical for the B cell response [27][28][29][30], and is modulated in part by CXCR4 [31] and its ligand CXCL12 [27][28][29][30], which are known to be regulated by HIF-1α in other cells [14][15][16]. CXCR4 signaling is regulated by transcriptional control, protein expression, and receptor internalization [32]. Interestingly, GC B cells have been shown to express surface CXCR4, however, they are unresponsive to CXCL12 signaling [33,34]. As GC B cells encounter O 2 levels, at times <1%, it is likely that CXCR4 responsiveness is in part controlled by an O 2 dependent post-translational mechanism, independent of CXCR4 transcription, translation or surface expression.
Based on the above data, we hypothesize that changes in O 2 tension as B cells migrate within the GC may directly control the localization and functional activation and differentiation of B cells. This hypothesis is strongly supported by the O 2 dependent regulation of several CXCR4 signaling components, including RGS1, which mediates HIF-1α induced CXCR4 uncoupling, along with p44/p42 MAPK and MKP-1 [34]. Focal adhesion kinase (FAK) is also critical for CXCR4 dependent migration of B cells [16], and is modulated by O 2 tension in smooth muscle cells [35]. In addition, CNI are known to interact directly and indirectly with the HIF-1α signaling cascade, and may have a significant role in disrupting the normal hypoxia-induced regulation of B cell migration. For example, CNI destabilize HIF-1α in glioma cells by stimulating prolyl hydroxylase activity [36], suggesting CNI have the capacity to disrupt hypoxic responses. Thus, there is also strong support for the additional hypothesis that hypoxia induced pathways are involved in modulation of CXCR4 signaling in B cells and CNI may disrupt these pathways.
In the following study, we demonstrate that migration of human and mouse B cells is regulated by chemokine receptor (CXCR4 and CXCR5) responsiveness via an O 2 sensing molecular switch, controlled by HIF-1α, which is independent of CXCR4/CXCR5 receptor expression or surface localization. Our data support the hypothesis that B cell migration is modulated by chemokine receptor sensitivity, controlled by stabilization of the master regulator HIF-1α at low O 2 levels (<4%), and indeed, we show genetically that HIF-1 α is necessary for this effect. Significantly, CyA destabilizes HIF-1α in both human and mouse B cells, partially restoring chemokine receptor responsiveness at low O 2 levels. These identical findings in both human and mouse cells may allow for a highly correlated assessment of in vivo immunological responses developing in lymph node and spleen using mouse models, as direct assessments are not possible in humans for anatomical and ethical reasons. Additional unbiased proteomics data suggests a switch in several metabolic processes potentially facilitating migration. This is consistent with the regulation of extracellular matrix and intrinsic apoptosis observed in the proteomic analysis. Transient restabilization of HIF-1α in CyA treated B cells temporarily restores the O 2 dependent molecular switch modulating B cell migration. These novel findings identify a direct, and potentially therapeutically targetable effect of CNI on B cell function, independent of indirect helper T cell effects.

Human and mouse b cell chemokine receptor (CXCR4 and CXCR5) hypo-responsiveness is induced by low O 2 levels and this correlates with HIF-1α stabilization
In order to examine B cell CXCR4 and CXCR5 responsiveness at different O 2 levels, we developed a novel, high throughput, in vitro transmigration assay system that combines a 96 well transwell plate format with a rapid luminescent readout of migratory cell numbers. Precise O 2 level control was achieved using two separate C-Chamber O 2 controlled incubator chambers (Biospherix, Parish, NY) to measure cell migration under varying levels of hypoxia. B cells in a developing immune response typically experience 1-4% O 2 within a local draining lymph node. To replicate these conditions, multiple C-Chambers were used to incubate separate chambers in parallel under three different O 2 levels in a single experiment. Thus, O 2 levels from the traditional 37 • C, humidified, CO 2 incubator provided 19% O 2 , and two C-chambers within the same incubator provided 1% and 4% O 2 .
We first examined the affect of 1 -19% O 2 levels on freshly isolated human peripheral blood and murine splenic B cells. Primary human peripheral blood B cells (n=3 donors) incubated at 19% O 2 for 24 hours migrate in responses to CXCL12 or CXCL13 ( Fig. 1a and Supplementary Figure 1a). However, CXCL12 and CXCL13 induced human B cell migration was inhibited at O 2 tensions of 1% or 4% O 2 . Mouse B cells exhibited an identical inhibition of chemotaxis at lower O 2 levels ( Fig. 1b and Supplementary Figure 1b), with robust chemotaxis to CXCL12 or CXCL13 at 19% O 2 , which decreases with 24 hours of incubation at 1% or 4% O 2 . Importantly, we observed no effect on human or mouse B cell surface CXCR4/CXCR5 expression (  Figure 2ab) irrespective of O 2 levels. These data strongly suggest that the decrease in chemotaxis is not due to decreased cellular viability or to loss of CXCR4/CXCR5 surface expression. The decreased chemotaxis was highly reproducible and the decreased chemotactic activity observed at lower O 2 levels is highly significant (p<0.001 for all pairwise comparisons). This finding of virtually identical CXCR4/CXCL12 or CXCR5/CXCL13 chemotaxis in mouse and human B cells suggests a highly conserved regulatory system. CXCR4 chemokine receptor sensitivity has been shown to be modulated by stabilization of the transcription factor HIF-1α, which occurs at low O 2 levels in most cells and induces nuclear translocation [6]. Thus, we next performed experiments to confirm HIF-1α stabilization correlates with CXCR4 or CXCR5 hypo-responsiveness. Figure 1 shows a western blot analysis of HIF-1α protein levels in nuclear lysates prepared from human peripheral blood B cells Fig. 1c) and mouse splenic B cells (Fig. 1d), incubated at the same O 2 levels used in the chemotaxis assay. Neither human nor mouse B cells incubated at 19% O 2 exhibit detectable levels of HIF-1α protein in nuclear lysates (Fig. 1cd). Thus there is a negative correlation between nuclear HIF-1α protein levels and human or mouse B cell chemotactic activity.

Low O 2 levels induce CXCR4/CXCR5 hypo-responsiveness and HIF-1α stabilization in RAMOS human B cells
The quantities of nuclear extracts that can be obtained from primary human or mouse splenic B cells significantly limit comparisons of multiple experimental conditions. Thus, we choose to perform more detailed western and proteomics analysis using the RAMOS human B cell line. First we confirmed that relationship between chemotaxis and HIF-1α stabilization seen in primary mouse and human cells was recapitulated in RAMOS cells. Identical to what we observed using primary B cells, the RAMOS human B cell line migrates in response to a CXCL12 or CXCL13 chemokine gradient at 19% O 2 in a dose dependent fashion. In contrast, incubation for 24 hours at 1% or 4% O 2 results in significantly decreased migration for both CXCL12 and CXCL13 gradients ( Fig. 2a and Supplementary Figure 1c). As with the primary cells, we also observed no effect on RAMOS cell viability, proliferation or surface CXCR4/CXCR5 expression ( Fig. 2d and Supplementary  Figure 2ab and c). This was irrespective of O 2 levels, again suggesting that the decrease in chemotaxis is not due to decreased cellular viability and/or proliferation, nor due to loss of chemokine receptor surface expression. In addition to examining viability and proliferation, we also wanted to determine if the low O 2 effects permanently alter the ability of RAMOS cells to respond to a CXCL12 gradient. We incubated RAMOS cells at 1% O 2 levels for 24 hours and then re-exposed the cells to 19% O 2 for 2, 4 or 6 hours and then subjected the cells to our chemotaxis assay. As is evident from Figure 2e RAMOS cells progressively regain the capacity to migrate in response to a CXCL12 gradient through 6 hours of re-exposure to 19% O 2 .
Similar to primary human and mouse cells, western blot analysis of nuclear lysates showed that HIF-1α protein was barely detectable after incubation of human RAMOS cells at 19% O 2 for 24 hours (Fig. 2b and c). In contrast, levels of HIF-1α protein increase over time in nuclear lysates when cells are incubated at lower O 2 levels (4% and 1%, a 5.3 and 7.6 fold increase as compared to 19%, respectively), demonstrating O 2 dependent HIF-1α protein stabilization and nuclear localization.
In addition to HIF-1α stabilization, HIF-2α has also been shown to be stabilized upon hypoxic exposure of several different cell types [37,38], but there is little data in human B cells. We performed a western blot analysis of HIF-2α protein stabilization in the RAMOS human B cells under varying O 2 conditions. We did not observe any significant HIF-2α stabilization in the nuclear lysates at low O 2 levels (Supplementary Figure 3). These results suggest that HIF-2α does not play a significant role in O 2 sensitive human B cell migration within the experimental time frame.

HIF-1α is both necessary and sufficient for CXCR4 hypo-responsiveness
We hypothesized that HIF-1α expression is necessary for B cell chemotactic hypo-responsiveness at low O 2 levels, and evaluated this hypothesis using shRNA silencing of HIF-1α. Stable RAMOS B cell lentiviral transfectants were generated expressing either HIF-1α shRNA, or a control non-specific shRNA vector and evaluated at the peak chemotaxis dose of CXCL12 (166 μg/mL; Fig. 2a). We focused on CXCR4/CXCL12 chemotaxis for the remainder of the manuscript, as we have shown that the O 2 /chemotactic effect is similar when comparing CXCR4 and CXCR5. For simplicity, we analyzed chemotaxis at the peak response (166 μg/mL). RAMOS cells expressing HIF-1α shRNA had increased chemotactic activity at 1% O 2 levels as compared to either the control shRNA transfectants or non-transfected cells (Fig. 3a). As predicted, HIF-1α protein levels were significantly decreased in HIF-1α shRNA transfected cells ( Fig. 3b and c). These genetic silencing results clearly demonstrate that HIF-1α expression is necessary for the observed decrease in CXCR4 responsiveness at low O 2 levels.
We also tested the hypothesis that HIF-1α is sufficient to induce chemotactic hypo-responsiveness at high O 2 levels using a HIF-1α over-expression system. RAMOS cells were stably transfected with a lentiviral vector expressing either HIF-1α or a control firefly luciferase. Post-transfection chemotaxis was evaluated at the peak dose of CXCL12 (166 μg/mL). As expected, RAMOS cells over-expressing HIF-1α demonstrated a significant decrease in chemotactic capacity at 19% O 2 as compared to either the untransfected or luciferase expressing RAMOS cells (Fig. 3d). In addition, HIF-1α protein levels were significantly increased in transfected cells over-expressing the HIF-1α construct ( Fig. 3e and f). These genetic over-expressing results clearly demonstrate that HIF-1α expression is sufficient to account for the observed decrease in CXCR4 responsiveness at low O 2 levels.

Cyclosporine a inhibits hypoxia dependent CXCR4 hypo-responsiveness by destabilizing HIF-1α
Cyclosporine A (CyA) destabilizes HIF-1α in glioma cells [36]. We hypothesized that CyA could similarly disrupt B cell HIF-1α stabilization, altering chemokine-induced B cell migration. Consistent with this hypothesis, and the above data, exposing RAMOS B cells to 1% O 2 inhibits CXCL12 mediated chemotaxis relative to that observed at 19% O 2 . However, treatment with increasing concentrations of CyA progressively increases CXCL12 mediated chemotaxis of RAMOS cells at 1% O 2 (Fig. 4a). These results are preserved across both primary human and mouse B cells, which also exhibit a similar increase in chemotactic activity after treatment with CyA ( Fig. 4d and e). The levels of CyA used in these experiments are physiologically relevant, as serum levels of 10-2000 ng/mL are achieved in vivo in immunosuppressed renal transplant recipients [39].
To assess the correlation between CyA's effect on migration and HIF-1α stabilization, we performed western blot analysis on RAMOS nuclear lysates over a range of CyA levels. Figure 4b and c demonstrate that treatment of RAMOS cells, with increasing amounts of CyA at 1% O 2 leads to a progressive decrease in HIF-1α levels. These results again directly link HIF-1α with chemotactic activity, demonstrating that CyA treatment results in decreased HIF-1α protein levels with increased CXCR4 responsiveness and chemotaxis, despite the fact that the B cells are incubated at 1% O 2 .
We hypothesize that the observed effect CyA has on B cell chemotaxis is mediated through HIF-1α. As such, we would anticipate that the HIF shRNA expressing cells would be insensitive to the CyA effect at low O 2 levels. To confirm this hypothesis, we incubated untransfected, control shRNA, or HIF shRNA expressing RAMOS cells at 1% O 2 for 24 hours with or without CyA and then subjected the cells to our chemotaxis assay. As is evident from Fig. 4f both the untransfected and control shRNA expressing cells demonstrated increased chemotactic activity upon treatment with CyA. In contrast, no increase in chemotactic activity was observed with the HIF shRNA cells after treatment with CyA. These results strongly support the hypothesis that CyA alters B cell chemotactic activity directly through modulation of HIF-1α.

Pharmacological stabilization of HIF-1α at 19% O 2 mimics the effects of low O 2 levels
The above results strongly suggest that an O 2 sensing molecular switch modulates B cell migration, correlates with and requires HIF-1α stabilization, and is sensitive to CyA treatment. We thus tested the hypothesis that pharmacologically stabilization of HIF at 19% O 2 would mimic the effect of incubating the cells at low O 2 . We performed experiments stabilizing HIF-1α at 19% O 2 using the prolyl-4-hydroxylase (PHD) inhibitor dimethyloxallyl glycine (DMOG). DMOG has been reported to selectively stabilize HIF-1α, with minimal effects on HIF-2α [38]. We found that incubation of RAMOS cells with DMOG at 19% O 2 resulted in decreased CXCL12 responsiveness (Fig. 5a), similar to the effect observed at low O 2 levels. DMOG had a similar effect on primary human and mouse B cells ( Fig. 5c and d). Western blot analysis ( Fig. 5b) of nuclear lysates from RAMOS cells confirms that the CXCL12 chemotactic effects highly correlates with HIF-1α protein stabilization, even when the cells were incubated at 19% O 2 .   RAMOS B cells with CyA for 24 hours, and during the last 3 or 6 hours of the incubation, DMOG was added. These experiments demonstrated that DMOG partially restored CXCR4 hypo-responsiveness of RAMOS cells, despite the presence of CyA (Fig. 5e).
Moreover, HIF-1α coordinately represses and activates genes such as ALDOC and PFKB1 to effect metabolic switching to a glycolytic state, providing biosynthetic substrates that facilitate proliferation [43].
We were particularly interested in proteins that were restored after treatment with CyA at 1% oxygen. The treatment with CyA restored the abundance of 3 proteins at 1% oxygen to a level similar to that in cells grown at 19% oxygen, selenocysteine-specific elongation factor (EEFSEC), Bardet-Biedl syndrome 12 protein (BBS12) and breast cancer anti-estrogen resistance protein 3 (BCAR3). EEFSEC, BBS12 and BCAR3 have been shown to be responsive to cellular oxygen metabolism [44][45][46]. Pathway analysis of proteins differentially expressed with respect to 19% oxygen level were mapped to 21 pathways including HIF-1α and glucose metabolism (Supplementary File 1). The two genes BNIP3 and EGLN1 involved in the HIF-1α transcription factor pathway were only differentially expressed in the absence of CyA. Moreover, the CyA-modulated gene ALDH3A2 which is involved in many metabolic pathways including fatty acid and amino acid metabolism. Thus differential analysis revealed several signaling components that could be dysregulated in RAMOS cells due to hypoxia and further controlled by CyA.

Protein association networks identify signaling components perturbed by CyA
Differential analysis identified only 0.4% proteins to be responsive to the oxygen levels and/or CyA. To expand the search for relevant proteins we sought to investigate functionally related proteins both from (a) our dataset and (b) in the literature. We constructed a co-expression network in which high degree nodes are involved in several metabolic processes ( Fig. 7a and Supplementary File 2 contains an edge list representation of the coexpression network). For example N6-adenosine-methyltransferase non-catalytic subunit (METTL14) is required in hypoxic stabilization of mRNAs [47] This is consistent with the known role of HIF-1α in regulating metabolic switches.
We next compared the protein co-expression network with a knowledge-based network constructed using the STRING database [48] (Supplementary File 3). 72 nodes were present in both, the knowledge-based and coexpression networks. The top 20 proteins, sorted by node degree, and overlaid on a histogram of the degree distribution of the co-expression network highlights highdegree nodes from the two networks (Fig. 7a). Several overlapping high-degree nodes such as core-binding factor subunit beta (CBFB), clathrin light chain B (CLTB), phosphatidylinositol 4,5-bisphosphate 3-kinase catalytic subunit gamma isoform (PIK3CG)) suggest their key role in the regulation of HIF-1α and CyA signaling. Corebinding factor subunit beta is a master-regulator of genes involved in hematopoiesis (e.g., RUNX1) and osteogenesis (e.g., RUNX2). Clathrin light chain B is involved in receptor-mediated endocytosis [49].
To identify modules of highly associated proteins, and by inference their genes, we used the edge-betweenness community detection algorithm [50]. We identified 24 communities with ≥10 proteins. These modules were characterized by the significantly enriched pathways (Fig. 7b). One community was of particular interest since it was restored after treatment with CyA at 1% O 2 (Fig. 7c) and was enriched for proteins involved in regulation of the actin cytoskeleton (see Supplementary File 4 for a list of the top 3 enriched pathways, with associated pvalues from the hypergeometric test). Specifically, this module contains the proteins PFN2 (profilin 2), IQGAP2 (Ras GTPase-like activating protein), and ITGAL (integrin alpha L). Both PFN2 [51] and IQGAP2 [52] regulate cell motility in response to intracellular signals. PFN2 is an epigenetic regulator of SMAD2/3 (Mothers against decapentaplegic homolog 2/3) [53], which in turn has been shown to promote metastasis in cancer cells [54]. SMAD2 and SMAD3 are also regulated by HIF-1α [55]. In the knowledge-based STRING network, SMAD3 is directly linked to HIF-1α. ITGAL is a target of the immunosuppressant drug efalizumab, which acts by reducing lymphocyte migration to sites of inflammation [56]. Thus, these proteins may be involved in the restoration of B cell chemotactic activity under hypoxic conditions by CyA, as measured by our chemotaxis assay, and suggest targets for further investigation.

Discussion
In this report we provide strong evidence for an O 2 sensing-HIF-1α dependent molecular switch that alters human and mouse B cell chemokine receptor (CXCR4 and CXCR5) responsiveness, and that this switch is disrupted by the calcineurin inhibitor cyclosporine A (CyA). These results have broad implications for our understanding of normal B cell immunological processes. Critically, these effects occur naturally at differing O 2 levels in vivo in B cells, an aspect of immunological function that is often overlooked, and the particular importance of O 2 levels in B cell specific responses has only recently begun to come to light. In support of this view, recent work has shown that GC hypoxia itself is important for efficient B cell class switching, plasma cell differentiation and abundant antibody production [17]. In addition, hypoxia within the GC light zone has been reported to be critical for generation of high quality antibody responses, in part by creating a more stringent threshold for vital survival signals, leading to the selection and survival of B cell clones producing higher affinity immunoglobulin [57]. CXCR4 signaling has also been shown to be critical for normal B cell function and development [30,32], especially proper B cell migration to, and localization within, the GC light or dark zones [58], and is dependent on CXCR4 signaling [28]. Increasingly hypoxia has been connected to CXCR4 regulation and function [33,34]. Importantly, GC B cells reside within a highly hypoxic environment, express surface CXCR4, but are unresponsive to CXCL12 [34]. Thus, these reports support our experimental findings.
Combined, our results and previous reports link proper B cell GC localization through O 2 dependent differential CXCR4 responsiveness to optimum B cell responses. In direct support of this hypothesis, we have found a HIF-1α dependent O 2 sensing molecular switch that controls the responsiveness of CXCR4 to CXCL12, independent of CXCR4 surface expression, directly linking GC hypoxic conditions to CXCR4 signaling. Importantly, we show that CXCR4 responsiveness is modulated through O 2 levels that are physiological and observed within the LN in vivo. In particular, LN O 2 levels vary from 4-5% O 2 near afferent vessels to <1% in portions of the GC [17,18,34,59]. Thus, it is very likely that this process is naturally coordinated and necessary for appropriate nodal and GC B cell dark/light zone localization, and that any disruption to this process would likely result in substandard B cell responses.
In addition to our functional findings, proteomics analysis identified signaling components involved in CXCR4 and HIF-1α signaling that are modulated by CyA. The differential expression analysis also identified signaling pathways responding to hypoxia and their key regulators. Several of these pathways are associated with HIF-1α, in agreement with findings of other groups [60][61][62][63]. However, several of the proteins reversing the effect of hypoxia with drug had lower fold-changes upon treatment with CyA. Association network analysis allowed us to probe CXCR4 and HIF-1α signaling specifically even when CXCR4 expression itself was not changed. This analysis revealed key regulators involved in the response to CyA and hypoxia. Moreover, we identified a module that reversed the effect of hypoxia upon CyA treatment, as well as a novel, highly intra-correlated protein module potentially involved in the responsiveness to CyA. Further work will be needed to assess the dynamics of this protein network.
Our findings suggest a non-T cell dependent mechanism for calcineurin inhibitor (CNI) induced hypogammoglobulinemia and B cell hypo-responsiveness. CNI are commonly utilized to prevent allograft rejection in solid organ transplant recipients. [37][38][39] However, the overall immunosuppressive effects of CNI results in decreased antibody production and hypogammoglobulinemia. This can lead to greater morbidity and mortality associated with common infectious agents, such as influenza viruses [64][65][66]. Vaccination is generally recommended for transplant recipients, however, CNI-based immunosuppression also suppresses vaccine responses [2,3]. This effect is largely thought to be due to CNI dependent inhibition of NFAT in helper T cells, leading to decreased T cell proliferation and IL-2 production [4]. Direct CNI effects on B cells have remained largely unidentified. However, we have now demonstrated that CyA has a novel and direct effect on B cell migration by disrupting the natural O 2 sensing molecular switch, destabilizing HIF-1α, and allowing B cells to preserve CXCR4 responsiveness at low (<1%) O 2 levels. Preservation of CXCR4 responsiveness disrupts coordinated GC B cell dark zone vs. light zone localization, which would suppress B cell responses in vivo. This is the first time, to our knowledge, that CyA have been shown to directly affect B cell migration, identifying a novel, and targetable, pathway by which CyA may directly affect B cell responses.
Improvement of immune function in CyA treated transplant patients during vaccination has the potential to increase vaccine efficacy, resulting in decreased morbidity and mortality associated with infection. However, as immune suppression is required to avoid allograft rejection, such immunological improvement would have to be transient and applied only during vaccine immune responses. We have identified a potential pathway to target, namely, reversal of CyA effects by transiently restabilizing HIF-1α. Indeed our in vitro data demonstrates that transient re-stabilization of HIF-1α through treatment of CyA affected B cells with the prolyl-4-hydroxylase (PHD) inhibitor dimethyloxallyl glycine (DMOG) results in retention of O 2 dependent changes in CXCR4 responsiveness. It is important to note that PHD inhibitors are clinically available for study, including FG-4592 and GSK1278863. Further pre-clinical study in mouse models would be required prior to clinical trials. Our data clearly demonstrate that both changes on O 2 levels and CyA levels identically affect the migratory capacity of mouse B and human B cells, strongly suggesting that direct in vivo data would be highly applicable to designing pre-clinical studies.

Conclusions
We have characterized an O 2 sensing-HIF-1α dependent molecular switch that alters human and mouse B cell CXCR4 responsiveness. This switch likely plays a significant role in GC B cell development and function. We have also identified a novel and direct CyA B cell affect whereby CyA directly interferes with this molecular switch by destabilizing HIF-1α at low (<1%) O 2 levels, preserving both human and mouse B cell CXCR4 responsiveness. The clinical implications of these results are potentially profound, as they identify a readily targetable pathway, through the transient use of clinically available PHD inhibitors (e.g. FG-4592 and GSK1278863), which may improve vaccine responses in vulnerable immune suppressed transplant patients.

Human subjects protection
This study was approved by the Research Subjects Review Board at the University of Rochester Medical Center (RSRB #71460). Informed consent was obtained from all participants. Research data were coded such that subjects could not be identified, directly or through linked identifiers, in compliance with the Department of Health and Human Services Regulations for the Protection of Human Subjects (45 CFR 46.101(b)(4)).

Human cell lines and primary b cell isolation
Human cell line: RAMOS (ATCC, Manassas, VA; CRL-1596) Burkitt's lymphoma B cell line was maintained in complete RPMI 1640 media supplemented with 10% FBS (cR10; Invitrogen, Carlsbad, CA). Primary human B cells: Human peripheral blood was collected by routine phlebotomy in heparinized vacutainers (BD Biosciences, San Diego, CA) and diluted 1:2 with PBS containing 2% FBS. Diluted blood was then layered over Ficoll-Paque PLUS (GE Healthcare, Chicago, IL) for density centrifugation using Sepmate tubes (Stemcell Technologies, Vancouver, BC, Canada). The lymphocyte layer was then washed twice with PBS (Invitrogen, Carlsbad, CA). PBMCs were then enriched for B cells by using the EasySep Human Pan-B cell Enrichment Kit (Stemcell Technologies, Vancouver, BC, Canada). PBMCs were incubated with an antibody cocktail containing CD2, CD3, CD14, CD16, CD36, Cd42b, CD56, CD66b, CD123, glycophorin A and dextran for 10 minutes. Streptavidin magnetic beads were added to the mix and the sample subjected to magnetic separation. The untouched enriched B cells were transferred to a clean tube, washed twice with PBS. B cells were then analyzed by flow cytometry showing >95% purity of the isolates.

Mouse splenic B cell isolation
Spleens from euthanized (carbon dioxide inhalation followed by exanguination after severing the axillary artery and thoracotomy to ensure cessation of respiration and cardiac function) female 9-10 week old C57/Bl6 mice (Taconic Labs, Rensselaer, NY) were removed, mechanically dissociated and cell suspensions passed over a 100μm nylon cell strainer (BD Biosciences, San Diego, CA). Tissue remaining in the strainer was gently disrupted with the plunger of a 1 mL syringe, and the strainer rinsed with 5 mL of IMDM media (Invitrogen, Carlsbad, CA) containing 200 mM EDTA. The single cell suspension (SCS) was then passed over a 70-μm nylon cell strainer (BD Biosciences, San Diego, CA) to further remove clumps. Excess RBCs were removed by magnetic positive selection (Imag; BD Biosciences, San Diego, CA) after incubation of the cell suspension with 50 μL Ter119 magnetic beads (BD Biosciences, San Diego, CA) per 1 x 10 7 total cells at 4 • C for 30 min. Unbound cells were transferred to new tubes and washed once with PBS with 2% FBS and 0.4% EDTA. B cells were enriched from the SCS using Mouse Pan-B Cell Isolation kits (Stemcell Technologies, Vancouver, BC, Canada). Briefly, cells were suspended at 1 x 10 8 cells/mL in PBS with 2% FBS and 0.4% EDTA. Normal Mouse serum was added at 50 μL/mL followed by the Mouse Pan-B cell Isolation Cocktail at 50 μL/mL. Cells were incubated for 10 minutes at RT. RapidSpheres, at 50 μL/mL, were added and incubated for 5 min. at RT. Samples were resuspended to a volume 2.5 mL with PBS with 2% FBS and 0.4% EDTA and placed on an Imag (BD Biosciences, San Diego, CA) for 5 min. Unbound cells were transferred to new tubes; flow cytometry indicated >95% purity of the isolates.
All research involving live, vertebrate animals was conducted in accordance with the Public Health Service Policy on Human Care and Use of Laboratory Animals. Female C57/Bl6 mice were obtained from Taconic Labs (Taconic Labs, Rensselaer, NY) and mice were maintained at the University of Rochester Medical Center Vivarium, an AAALAC certified Vivarium (Animal Welfare Assurance Number is A-3292-01), under their established guidelines, including isolation, feeding, recovery procedures, and euthanasia in accordance with Federal regulations. All experimental procedures for animals were approved by the Institutional Animal Care and Use Committee (IUCAC; protocol number UCAR-2011-055E), and all personnel working with the animals were trained and certified by the IUCAC and Vivarium staff.

Generation of stable HIF-1α shRNA and HIF-1α over-expression RAMOS cell lines
Lentiviral particles containing with HIF-1α shRNA or a control shRNA encoding a scrambled shRNA sequence that will not lead to the specific degradation of any cellular message; or HIF-1α activation particles or control particles expressing firefly luciferase, were obtained commercially (Santa Cruz Biotechnology, Dallas, TX). Stable transfectants were generated and selected using either puromycin (shRNA constructs) or blasticidin S ( over-expression constructs) according to the manufacture's recommendations (Santa Cruz Biotechnology, Dallas, TX).

RAMOS B cells or human and mouse primary B cells
were resuspended in cR10 media at 5 × 10 4 cell/mL and added to appropriate sized tissue culture plates or flasks depending on total volume. To select flasks, CyA and/or DMOG (Sigma-Aldrich, St. Louis, MO) was added at the concentration indicated in the text and/or figure legends. Plates or flasks were then placed within controlled environment C-chamber 37 • C, 5% CO 2 incubator inserts (Biospherix, Parish, NY) that were equilibrated O/N at the indicated O 2 levels or placed within a traditional 5% CO 2 , 37 • C incubator (19% O 2 ). Cells were incubated for 24 hours, removed and subjected to either the chemotaxis, western blot or proteomics assays.

Chemotaxis assay
The human B cell line, RAMOS cells, or magnetic bead enriched primary human peripheral blood or mouse splenic B cells were resuspended in freshly prepared migration media (PBS + 1% BSA), and 5 x 10 5 cells (in 75 μL) loaded into the upper chambers of 5 μm polycarbonate 96 well transwell plates (Corning, Corning, NY). In triplicate wells, migration media containing varying concentrations of CXCL12 or CXCL13 was added to the bottom chambers of the transwells (225 μL per well). Plates were incubated for 1 hour at appropriate O 2 levels in the C-chamber incubator inserts (Biospherix, Parish, NY), at 37 • C and 5% CO 2 . Upper transwells were removed and 100 μL of Cell Titer Glo (Promega, Madison, WI) added to each bottom well and incubated at RT for 10 min. Luminescence was readout on a SynergyTM HT microplate reader (BioTek Instruments, Winooski, VT) and relative luminescent units (RLU) reported. Supplementary Figure 4 demonstrates the direct correlation between RLU and cell numbers, through a 4-log 10 range, independent of O 2 levels.

Western blot
Nuclear lysates were prepared from B cell pellets using NE-PER nuclear and cytoplasmic extraction kits and protein concentration determined using micro BCA protein assay kits, both according to the manufacture's recommendations (Thermo Fisher Scientific, Waltham, MA). Nuclear lysates (10μg per lane) were resolved on precast NuPAGE 4-12% Bis-Tris protein gels (Invitrogen, Carlsbad, CA) and proteins transferred to PVDF membranes (Bio-Rad, Hercules, CA). The membrane was then blocked with 1X TBS-T (Tris-Buffered Saline (Bio-Rad, Hercules, CA)) with 0.05% Tween-20 (Sigma-Aldrich, St. Louis, MO) + 5% blotting grade non-fat dry milk (Bio-Rad, Hercules, CA) for 1 hour at RT. Blots were washed 3 x with TBS-T and then primary antibody diluted as indicated in TBS-T + 5% milk. For human B cell blots a 1:250 dilution mouse anti-human HIF-1α (BD Biosciences, San Diego, CA) or for mouse blots a 1:100 dilution goat antimouse HIF-1α was added and the blots incubated O/N at 4 • C. As a loading control, blots were also probed with a 1:200 dilution mouse anti-actin primary antibody and the blots incubated O/N at 4 • C. For detection of the human HIF-1α or actin blots, a 1:1,000 dilution, in PBS + 5% milk, of horseradish peroxidase (HRP) conjugated goat anti-mouse Ig secondary was used. For mouse HIF-1α, a 1:2,000 dilution HRP conjugated donkey antigoat IgG was used. Secondaries were added and the blots incubated for 1 hour at RT. Blots were washed 3 X with TBS-T, rinsed with distilled water and freshly prepared ECL substrate (Thermo Fisher Scientific, Waltham, MA) added and the blots were imaged on a Chemi-Doc MP imaging system. Densitometry was performed using Image Lab software version 6.0.1 (Bio-Rad, Hercules, CA).

Sample preparation
RAMOS cells were re-suspended in cR10 media at 5 x 10 5 cell/mL at a final volume of 10 mL in T25 tissue culture flasks (BD Biosciences, San Diego, CA) and either left untreated or cyclosporine A (Sigma-Aldrich, St. Louis, MO) added to a final concentration of 1 μg/mL. Pairs of untreated or CyA treated cells were then incubated at either 19% or 1% O 2 levels for 24 hours. Cells were then rapidly harvested into 50 mL conical tubes (BD Biosciences, San Diego, CA) containing 30 mL of ice cold PBS. Cells were pelleted and washed 3 x 10 mL ice cold PBS. After the final wash, cells were re-suspended in 1 mL ice cold PBS and transferred to a 1.5 mL microcentrifuge tubes and pelleted. Supernatants were removed and the cell pellets flash frozen in liquid nitrogen. Cell pellets were stored at −80 • C until samples were lysed for analysis.
Lysis buffer consisted of 5% SDS, 50 mM triethylammonium bicarbonate (TEAB). Cell lysis or each pellet was done by adding 100 μL of lysis buffer per 10 6 cells, followed by vortexing and sonication using a QSonica sonicator. Sonication cycles consisted of 5 X 10 second sonications with one minute incubations on ice between each cycle. Samples were centrifuged for 5 minutes at 16,000 x g, and the supernatant was collected. Protein concentration was determined by BCA (Thermo Fisher Scientific, Waltham, MA). 25 μg of protein from each sample was removed, and brought up to 25 μL in lysis buffer. Disulfide bonds were reduced by addition of dithiothreitol (DTT) to 2 mM, followed by incubation at 55 • C for 1 hour. Alkylation was performed by adding iodoacetamide to 10 mM and incubating at room temperature in the dark for 30 minutes. 12% phosphoric acid was added to a final concentration of 1.2%, followed by the addition of 6 volumes of 90% methanol, 100 mM TEAB. The resulting solution was added to S-Trap micros (Protifi, Huntington NY), and centrifuged at 4,000 x g for 1 minute. The S-Traps were washed twice with 90% methanol, 100 mM TEAB. 20 μL of 100 mM TEAB containing 1 μg of trypsin was added to the S-Trap, followed by an additional 20 μl of TEAB.
Samples were placed in a humidity chamber at 37 • C and were allowed to digest overnight. The S-Trap was centrifuged at 4,000 x g for 1 minute to collect the digested peptides. Twenty μL of 0.1% trifluoroacetic acid (TFA) in acetonitrile was added to the S-Trap, and one more centrifugation step was done, and the solutions were pooled, frozen, and dried down in a Speed Vac (Labconco, Kansas City, MO). Tandem mass tag (TMT) ten-plex reagents (0.2 mg) (Thermo Fisher Scientific, Waltham, MA) were removed from 20 • C and allowed to reach room temperature prior to dissolving each tag in 20 μL of acetonitrile. Samples were re-constituted in 25 μL of TEAB, the TMT tags were added to the samples, and incubated at room temperature for one hour. The reaction was quenched by the addition of 3 μL of 5% hydroxylamine. 20% (5 μg) of each sample was combined, frozen, and dried in the Speed Vac.
To increase coverage, samples were fractionated using C18 spin columns. Columns were conditioned with acetonitrile, followed by equilibration with 100 mM ammonium formate (AF), pH 10. The samples were resuspended in 50 μL of AF and added to the spin columns. After washing the columns with water and then the AF buffer, samples were eluted with 10%, 12.5%, 15%, 17.5%, 20%, 22.5%, 25%, and then 50% acetonitrile in the AF buffer. These fractions were then frozen, dried down, and re-suspended in 0.1% TFA in water, and placed into autosampler vials.

Data acquisition
Peptides were injected onto a 30 cm C18 column packed with 1.8 μm beads (Sepax), with an Easy nLC-1000 HPLC (Thermo Fisher Scientific, Waltham, MA), connected to a Q Exactive Plus mass spectrometer (Thermo Fisher Scientific, Waltham, MA). Solvent A was 0.1% formic acid in water, while solvent B was 0.1% formic acid in acetonitrile. Ions were introduced to the mass spectrometer using a Nanospray Flex source operating at 2 kV. The gradient began at 6% B and held for 2 minutes, increased to 30% B over 85 minutes, increased to 50% B over 10 minutes, then ramped up to 70% B in 4 minutes and was held for 5 minutes, before returning to starting conditions in 4 minutes and re-equilibrating for 10 minutes, for a total run time of 120 minutes at a 300 nL/minute flow rate. The Q Exactive Plus was operated in data-dependent mode, with a full scan followed by 10 MS/MS scans. The full scan was performed over a range of 400-1700 m/z, with a resolution of 70,000 at m/z of 200, an automatic gain control (AGC) target of 10 6 , and a maximum injection time of 50 ms. Peptides with a charge state between 2-5 were picked for fragmentation. Precursor ions were fragmented by higher-energy collisional dissociation (HCD) using a collision energy of 35 and an isolation width of 1.0 m/z. MS2 scans were collected with a resolution of 35,000, a maximum injection time of 120 ms, and an AGC setting of 1e 5 . The fixed first mass for the MS2 scans was set to 110 m/z to ensure TMT reporter ions were always collected. Dynamic exclusion was set to 25 seconds.

Data analysis
Raw data was searched using the SEQUEST search engine within the Proteome Discoverer software platform, version 2.2 (Thermo Fisher Scientific, Waltham, MA), using the SwissProt human database [67]. Trypsin was selected as the enzyme allowing up to 2 missed cleavages, with an MS1 mass tolerance of 10 ppm, and an MS2 mass tolerance of 0.025 Da. Carbamidomethyl on cysteine, and TMT on lysine and peptide N-terminus were set as a fixed modifications, while oxidation of methionine was set as a variable modification. Percolator was used as the false discovery rate (FDR) calculator, filtering out peptides which had a q-value greater than 0.01. Reporter ions were quantified using the Reporter Ions Quantifier node, with an integration tolerance of 20 ppm, and the integration method being set to "most confident centroid". Protein abundances were calculated by summing the intensities of the reporter ions from each identified peptide, while excluding any peptides with an isolation interference of 30% or more. Low abundance proteins with less than one count per experiment/replicate were removed resulting in 5048 proteins. Proteins' abundance was log 2 -transformed, and normality was confirmed by diagnostic plots such as quantile-quantile plot. The differential expression analysis was performed using the implementation of the empirical Bayes statistic in the R package limma [68] across all conditions. Differentially expressed (DE) proteins are defined as p-value of ≤0.01 and a log 2 -fold change of ≥0.5. Pathway analysis of DE proteins was performed using enricher from R package clusterProfiler [69] and 1392 canonical gene sets curated in MSigDB [70].

Network analysis
To evaluate the associations among proteins and neighborhood of DE proteins, we assembled a co-expression network using 3919 proteins with coefficient of variation ≥0.01 (Supplementary File 2). The edges were weighted by the absolute Pearson correlation coefficient between the connected nodes. Each edge in the co-expression network was assigned a z-score by permutation of sample labels 50 times. We retained only edges with a weight (Pearson correlation coefficient) that had an absolute z-score of ≥1.5, i.e., the real correlation was at least 1.5 standard deviations away from the mean of correlations calculated by perturbing sample labels. Edges with a weight ≤0.5 were removed. This resulted in a network of 1408 nodes and 1211 edges. The implementation of the edgebetweenness community detection algorithm in the R package igraph [50,71] was used to identify communities with ≥10 genes. We performed a gene set enrichment analysis using the hypergeometric test. To further map genes/proteins known to be functionally related to HIF1A and CXCR4, 33 pathways from KEGG [72], Biocarta [73] and Reactome [74] were identified. The compiled list of 1504 unique genes from these pathways was used in the STRING database [48] to build a network where the edges are weighted based on strength of experimental evidence of interaction, co-expression measured using data in public domain, and database annotations. This network is referred to as the knowledge based network.