Next Article in Journal
Evaluating the Effects of Subnormothermic Perfusion with AP39 in a Novel Blood-Free Model of Ex Vivo Kidney Preservation and Reperfusion
Previous Article in Journal
Co-Adjuvant Therapy Efficacy of Catechin and Procyanidin B2 with Docetaxel on Hormone-Related Cancers In Vitro
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Upregulated Chemokine and Rho-GTPase Genes Define Immune Cell Emigration into Salivary Glands of Sjögren’s Syndrome-Susceptible C57BL/6.NOD-Aec1Aec2 Mice

1
Department of Infectious Diseases and Immunology, College of Veterinary Medicine, University of Florida, P.O. Box 100125, Gainesville, FL 32610, USA
2
Division of Allergy, Immunology and Rheumatology, SUNY Buffalo School of Medicine, 875 Ellicott Street, Buffalo, NY 14203, USA
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(13), 7176; https://doi.org/10.3390/ijms22137176
Submission received: 8 May 2021 / Revised: 16 June 2021 / Accepted: 17 June 2021 / Published: 2 July 2021
(This article belongs to the Section Molecular Pathology, Diagnostics, and Therapeutics)

Abstract

:
The C57BL/6.NOD-Aec1Aec2 mouse is considered a highly appropriate model of Sjögren’s Syndrome (SS), a human systemic autoimmune disease characterized primarily as the loss of lacrimal and salivary gland functions. This mouse model, as well as other mouse models of SS, have shown that B lymphocytes are essential for the development and onset of observed clinical manifestations. More recently, studies carried out in the C57BL/6.IL14α transgenic mouse have indicated that the marginal zone B (MZB) cell population is responsible for development of SS disease, reflecting recent observations that MZB cells are present in the salivary glands of SS patients and most likely initiate the subsequent loss of exocrine functions. Although MZB cells are difficult to study in vivo and in vitro, we have carried out an ex vivo investigation that uses temporal global RNA transcriptomic analyses to profile differentially expressed genes known to be associated with cell migration. Results indicate a temporal upregulation of specific chemokine, chemokine receptor, and Rho-GTPase genes in the salivary glands of C57BL/6.NOD-Aec1Aec2 mice that correlate with the early appearance of periductal lymphocyte infiltrations. Using the power of transcriptomic analyses to better define the genetic profile of lymphocytic emigration into the salivary glands of SS mice, new insights into the underlying mechanisms of SS disease development and onset begin to come into focus, thereby establishing a foundation for further in-depth and novel investigations of the covert and early overt phases of SS disease at the cellular level.

1. Introduction

Autoimmunity is generally recognized as a multi-step process initiated by environmental triggers that first activate an innate inflammatory reaction, then subsequently progress to an adaptive immune response, in genetically or physiologically predisposed hosts, that targets critical biological functions whose subsequent dysfunctions eventually result in an overt clinical pathology [1,2,3,4,5,6]. As such, autoimmune diseases are considered to have both an early covert disease phase and a late overt disease phase. Unfortunately, patients most often present in clinic only after the adaptive immune phase is active and irreversible pathology is occurring. For this reason, it remains necessary to identify molecular and cellular processes involved with the early covert inflammatory phase of autoimmunity as a basis to better understand the initiation of the adaptive immune response phase. This is especially true for Sjögren’s Syndrome (SS) where the time between disease onset and eventual diagnosis ranges from 4 to 10 or more years [7,8,9]. Furthermore, the pathological profile changes dramatically during the development and onset of disease in this time period.
SS is a highly debilitating, yet fascinating systemic autoimmune disease marked by leukocytic infiltrations of the salivary and lacrimal glands with a concomitant loss of exocrine secretion leading to clinical symptoms of severe dry mouth and dry eye diseases (reviewed in [9]). Occasionally, patients will develop kidney and lung infiltrations, central nervous system manifestations, and even lymphomas, the latter most commonly non-Hodgkin’s B cell lymphomas [6]. While considerable effort has long focused on defining the general pathophysiology and accompanying symptoms for disease diagnosis in patients, only recently has a more concerted effort surfaced to better define the molecular-based bioprocesses for the development and subsequent onset of SS autoimmunity and lymphomagenesis, research that has mostly been carried out utilizing a wide variety of mouse models [10,11,12].
Over the past couple decades, we have investigated multiple mouse models that exhibit progressive and spontaneous SS-like pathologies to define the early molecular events responsible for disease development and onset [13,14,15,16,17,18,19,20,21]. In general, these studies have defined this SS-like disease development as: (a) aberrant proteolytic enzyme activity most likely associated with acinar cell apoptosis, (b) a progressive loss of saliva and tear flow rates accompanied by increased protein content due in part to less fluid secretions, (c) decline in amylase and carbonic anhydrase activities, (d) appearance of autoantibodies, and (e) loss of acinar tissue, all manifestations occurring concomitantly with increasing glandular lymphocytic infiltrates. More recent studies have now shown a direct correlation between an upregulated expression of IL14α and late-stage B cell lymphomagenesis [13,14]. Importantly, global temporal transcriptome studies have verified these reported immuno-pathophysiological temporal changes, as well as the fact that full development of disease is clearly a multi-phase process involving an innate response, an adaptive response, and in some individuals, a late lymphomagenesis phase. This progression of autoimmune pathologies is quite significant when compared against sex- and aged-matched SS-non-susceptible (SSNS) control mice that do not develop clinical SS-like diseases.
While earlier studies revealed that inactivating the B cell receptor [14,15], interferon [16], IL-4 (interleukin-14) [17] and/or the alternate C’ system [18] prevents SS-like disease in our mouse models, our recent studies with B6.Il14α transgenic (TG) mice [19,20] have now indicated that elimination of marginal zone B (MZB) cells, or the blocking of lymphotoxin activity that is required for MZB cell ontogeny in marginal zones (MZ), prevents the development of SS-like disease, including lymphomagenesis [21,22]. MZB cells are an unique subpopulation of bone marrow derived B cells characterized by limited expression of immunoglobulin variable region genes that produce predominantly IgM antibodies, many of which are self-reactive [23]. In addition, they are strategically located within mucosal surfaces, function as innate and/or transitional cells capable of rapid responses to both T cell-independent and T cell-dependent antigens, and regulate the activation of subsequent adaptive immune responses by T and B2 lymphocytes in association with monocytic and neutrophilic antigen-presenting cells (APCs) (reviewed in [23]). Functionally, MZB cells differentiate from transitional type-1 (T1) B cells under the influence of low affinity B cell receptor (BCR) signals and transcription factors, especially Notch-2 [24]. They are enriched within splenic MZs, retained there by interactions between MZ integrin and MZB cell integrin receptors (e.g., MADCAM1 (mucosal cell adhesion molecule-1), LFA1 (lymphocyte function-association-1), ICAM (intercellular adhesion molecule-1), VLA-4 (very late antigen-4) and SIP1 (sphingosine-1-phosphate) [25]. MZB cells respond rapidly to Toll-like receptor (TLR) signaling, especially TLR2, 4 and 9 [25,26]. While MZB cells are important for responding to pathogens [27,28], they have been identified within the salivary glands of patients with salivary gland disease where they secrete cytokines cytotoxic to salivary gland cells [22,29,30].
Despite the noted difficulty in studying MZB cells either in vivo or in vitro, we have utilized ex vivo temporal global transcriptomic microarrays to identify MZB cell markers in order to better define MZB cell profiles and their temporal presence in SS-like disease of the SS-susceptible (SSS) C57BL/6.NOD-Aec1Aec2 mouse model. In the present report, we have analyzed global microarray data to identify the temporal genetic profile associated with the emigration of cells, predominantly MZB cells, into the salivary glands during early stages of SS-like disease development and onset in C57BL/6.NOD-Aec1Aec2 mice.

2. Results

2.1. Salivary Gland Histology Pre- and Post-Lymphocyte Infiltration in C57BL/6.Aec1Aec2 Mice

One of the main manifestations of SS/SS-like disease is the progressive and age-dependent leukocytic infiltrations of the salivary and lacrimal glands, resulting in the progressive formation of lymphocytic foci (comparison of Figure 1A to Figure 1B). While this emigration of B and T cells is considered an important bioprocess and is easily visualized histologically, little is known about its molecular basis. However, results from our earlier studies using several SS-susceptible mouse models revealed an absolute requirement for B lymphocytes in order for an overt clinical SS-like disease to develop, irrespective of whether T lymphocytes infiltrate the salivary or lacrimal glands [12]. Furthermore, multiple studies have now suggested that MZB cells are the B cell population responsible for disease development [22,31], a concept not only supported by the studies in B6.Il14α TG mice discussed above, but also by our recent findings in C57BL/6.NOD-Aec1Aec2 mice of a coordinated activation of the Notch2 and Type 1 interferon pathways in the absence of an upregulated Notch1 pathway [32]. This profile, in turn, predisposes mice to a prolonged upregulated Type1 interferon response characteristic of SS disease [33], our hypothesized anti-dsRNA reaction underlying SS-like disease [34], and induction of high expression of C´ (complement) receptors for the alternate C’ activity [35,36].

2.2. Gene Expression Profiles Defining, in Part, the Early Leukocyte Emigrations into the Salivary Glands of C57BL/6.NOD-Aec1Aec2 Mice

2.2.1. Coordinated Upregulation of Chemokine > Chemokine Receptor Genes with LFA1 and Rho-GTPase Upregulations

An interesting phenomenon that distinguishes MZB cells from follicular B cells, particularly in mice, is the relatively non-circulatory state of MZB cells compared to the high migratory propensity of follicular B cells [37]. This difference, as stated by Lu and Cyster [25], occurs because of the upregulated surface levels of integrin expression on MZB cells, especially LFA-1. LFA-1 is a heterodimeric integrin comprised of subunits ItgαL and Itgβ2 that binds to ICAM-1 and VCAM1. Development of splenic MZ areas during ontogeny requires the lymphotoxin-α1:β2 integrin-mediated induction of VCAM-1 and ICAM-1, the latter a primary receptor for LFA-1. At the same time, this high level of integrin expression on MZB cells requires higher threshold levels of Cxcl13 (chemokine (CXC motif) ligand) signaling for Cxcr5 (chemokine (CXC motif) receptor) expressing MZB cells to be dislodged from MZs and emigrate to other tissue sites. As presented in Figure 2, our transcriptomic data reveal that genes encoding Icam1, Vcam1, ItgαL, Itgβ2, Cxcl13, and Cxcr5 (the Cxcl13 receptor) are all upregulated in the salivary glands at the same time (i.e., 16 weeks of age), along with Ccr6, Ccr7, and Ptk2β, strongly suggesting an influx of immune cells actively recruited by glandular tissue signaling at the early stage of disease. In addition, the development and accumulation of MZB cells in MZs, followed by their subsequent migration, appear to rest in the requirement for the Dock2 > Rho-GEF > Ptk2 molecular pathway [28,38]. At the same time, it is important to note the published work of Gotoh et al. [39] showing that Dock2 is also indispensable for the co-migration of plasmacytoid dendritic cells (pDC), a highly relevant fact considering the extensive type 1 IFN-signature seen in exocrine glands of both human and mouse SS diseases.

2.2.2. Transcriptomic Profile of the Rho-GTPase Family of Proteins

Cell emigration is highly dependent on two events: chemokine signaling to direct leukocytes to sights of injury, and Rho-GTP (Ras homolog-G protein bound-guanine nucleotide) family protein (Rho-GTPase) activations. In a quiescent state, the Rho-GTP proteins mostly localize to cellular membranes but dissociate from these membranes during initiation of migration. In mice, the Rho-GTP family consists of 21 members subdivided into several subfamilies: 11 Rho proteins, 3 Rnd proteins, 3 Rac proteins, 3 Rhobtb proteins and Cdc42 [38]. Interestingly, as shown in Figure 3, a single molecular entity from each subfamily exhibits a temporal activation (i.e., Rhoc, Rhou, Cdc42, Rnd1, and Rhobtb1), except for Rsf3, where both Rac2 and Rac3 are upregulated. These activations temporally mimic the development of SS-like pathology within the salivary glands, with Rhoc, Cdc42, Rac3, and Rhou exhibiting early upregulation starting at 8 weeks of age, and with Rac2 and Rhobtb1 showing an upregulation at 16 weeks, an age when the adaptive response is observed in the salivary glands of these mice.

2.2.3. Transcriptomic Profile of the Rho-GAP and Rho-GEF Subfamilies of Rho-GTPase Proteins

The Rho-GTP proteins are further divided into two functional classes, those that are incapable of hydrolyzing GTP (i.e., Rhoh and Rnd subfamily members that constitutively bind GTP), and those that hydrolyze GTP through an opposing equilibrium of GTPase-activating proteins (GAPs) and Rho-specific guanine nucleotide exchange factors (GEFs) [38,39,40]. GEFs activate Rho-GTPases by exchanging a bound GDP with a GTP, while GAPs inactivate Rho GTPases vis a vis GTP by catalysis. Comparative temporal gene expression profiles for GAP and GEF molecules in the salivary glands of C57BL/6.NOD-Aec1Aec2 mice versus C57BL/6J mice, summarized in Figure 4, reveal a highly restricted activation with only 7 of the 17 Arhgef (Rho GTPase exchange family) gene members and 8 of the 23 Arhgap (Rho GTPase activating family) gene members upregulated on the arrays. These gene set profiles are markedly different from the profiles seen in the salivary glands of C57BL/6J mice, where there is a marked absence of upregulated gene expressions within the same gene sets.
In addition to the Arhgef proteins, a second subset of Rho-GEF proteins are the DOCK (dedicator of cytokinesis) family proteins (reviewed in [41]). DOCK molecules, also known as CZH proteins, possess a functional CZH2 domain that promotes the exchange of GDP to GTP in the formation of Rho-GTPases. Based on sequence homologies represented within the eleven identified DOCK molecules, this molecular family of proteins is currently divided into four major subgroups, (A) Dock1, 2 and 5, (B) Dock3 and 4, (C) Dock6, 7, and 8, and (D) Dock9, 10 and 11. Their temporal gene expression profiles in SS-susceptible C57BL/6.NOD-Aec1Aec2 mice relative to their gene expression profiles in SS-non-susceptible C57BL/6J control mice are shown in Figure 5. As might be expected in a family of proteins critical for multiple cellular functions, but especially cell migration and homeostasis, C57BL/6.NOD-Aec1Aec2 mice exhibit a distinct transcriptomic profile involving both in activation and in temporal expressions: Dock3, Dock4, Dock6, and Dock9 remain quiescent, Dock1, Dock5, Dock7 and Dock8 exhibit relatively long-term temporal upregulated expressions peaking between 8 and 12 weeks of age, while Dock2, Dock10 and Dock11 exhibit a unique short-term upregulated profile occurring around 16 weeks of age, a timepoint when the salivary glands are first being invaded by leukocytes. In contrast, SS-non-susceptible C57BL/6J mice show only four DOCK genes with upregulated expressions: Dock4, Dock7, Dock9 and Dock10, and in each case maximum expressions occurred between 8 and 12 weeks of age, consistent with the occasional leukocyte infiltrations observed in the salivary glands of C57BL/6J mice.

2.2.4. Transcriptomic Profiles of Factors Associated with Signal Transduction Pathways of Rho-GAP and Rho-GEF Subfamily Proteins of the Rho-GTPases

Arhgef molecules interact with several unique molecular families, including the Rac-GTPs, Rho-GTPs and cell division cycle-42 (Cdc42)-GTPs (reviewed in [42,43]). These molecules, in turn, interact with Wasf (WASP Family member 3) protein, Rock1 and/or Rock2 (Rho-associated coiled-coil containing protein kinase 1 and protein kinase 2) proteins, or Was (Wiskott–Aldrich molecule) proteins, respectively. Various functions regulated by these pathways include focal adhesions, myosin activations, lamellipodia formation, membrane ruffling, actin stabilization and polymerization, and filopodia formation, each an important factor in cytokinesis and cell migrations. Our transcriptome analysis, presented in Figure 6, reveals that the genes encoding for Wasf3, both Rock1 and Rock2, and Was are each upregulated in the salivary glands of C57BL/6.NOD-Aec1Aec2 mice. In contrast, only Wasf3 is highly upregulated and Rock1 minimally upregulated in the salivary glands of C57BL/6J control mice. These data point to an activation of Rho-GTP pathway(s) in C57BL/6.NOD-Aec1Aec2 mice at 8 weeks of age, followed by an activation of Cdc42 pathway(s) at 16 weeks of age in-line with the two early disease phases for SS-like disease in the salivary glands of our SSS C57BL/6.NOD-Aec1Aec2 mice.
Rho-GTPases are known to interact with multiple proteins that form signal transduction pathways that regulate activations of cellular biological functions. These include members of the Cdc42 protein family (e.g., Cdc42bpa, Cdc42bpb, Cdc42ep4, Cdc42se1 and Cdc42se2), the Vav family (e.g., Vav2, Vav3) and the Baiap family (e.g., Baiap2, Baiap2l1). In addition, genes encoding other interacting proteins include Setd6 (Prex1) (Ser domain containing 6), Scrib1 (Scribble homolog), Flii (Flii actin remodeling protein), Grb2, Creb3 (Camp response element binding protein), and Foxo3a, the majority of which exhibit upregulated gene profiles in the salivary glands of C57BL/6.NOD-Aec1Aec2 mice (Figure 7). In contrast, genes encoding Tiam1 (Rac1 associated GEF1) and Srgap1 (Synaptic Ras GTPase-activating protein 1), known to suppress Rho-GEFs, demonstrate no upregulated expressions. At the same time, the genes encoding Stat6 and Stat1 (Signal transduction and activator of transcription factors 6 and 1) show temporal upregulated expression profiles consistent with their known functional activities in SS, i.e., Stat6 being important during early-stage disease and Stat1 being important in the later stages of disease.

3. Discussion

The focus of the present transcriptome study has been to profile the underlying molecular bases for emigration of the autoimmune cellular infiltrates to the salivary glands during the early stages of disease development and onset of SS-like pathology observed in the C57BL/6.NOD-Aec1Aec2 mouse model. The data presented indicate an active upregulation of signal transduction pathways defining cell immigrations involving specific GTPases. This differentially expressed gene profile occurs from 8 weeks of age (Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7), a timepoint corresponding with the hypothesized appearance of MZB cells within the salivary glands. In contrast, an active chemokine expression, together with an environment for leukocyte adhesion in the salivary gland, appears to occur uniquely at 16 weeks of age (Figure 2), the time point corresponding to the infiltration of increasing numbers of leukocytes and the establishment of lymphocytic foci. Thus, we propose that this early disease stage can be divided into a ¨passive¨ emigration involving MZB cells and an ¨active¨ emigration involving recruitment of additional leukocyte populations. Nevertheless, the precise cellular functions attributed to MZB cells in SS still require better definition, since MZB cells are considered better responders to pathogens than to self-antigens [44,45], and responses toward microorganisms or other environmental triggers have yet to be ruled out as underlying initiators of SS. In this regard, our earlier transcriptome studies of the interferon-signature and activated TLR pathways point to the possibility of an anti-dsRNA virus response [11,13,14].
As an immune cell population, MZB cells show: (a) upregulated Notch2, Adam10, RBP-J and their associated molecular pathways for B cell ontogeny and differentiation [36]; (b) expressions of VLA-4 and SIP1 [45]; (c) rapid responses toward LPS while focusing T cell responses toward pathogens at restricted locations [46]; (d) strong responses to TLR2, TLR4 and TLR9 signaling [47,48]; (e) a capacity for delivering antigens to follicular dendritic cells in germinal centers [49], presentation of antigen efficiently to T cells [50], and regulation of both T cell-independent and T cell-dependent B and T cell responses [51]. Phenotypically, both mouse and human MZB cells are identified as IgMhi, IgDlow, CD21hi, and CD23−/−, but mouse MZB cells also express Cd1dhi, while human MZB cells express CD1c and CD27. Interestingly, although the microvasculature of mouse and human splenic marginal zones (MZs) are anatomically different, with human spleens apparently lacking a marginal sinus that can influence migratory pathways, both species provide an environment for MZB cells to interact with APCs in their respective stromal reticular cell networks [52], thereby permitting efficient immune surveillance of blood-borne antigens. Furthermore, both species separate MZB and follicular B cells via the periarteriolar lymphoid sheaths (PALS). Despite these important differences noted between mouse and human MZB cell ontogeny and behavior that include sites of potential maturation, resident locations, and antibody isotype repertoires (especially to T cell-dependent antigens), both mouse and human MZB cells tend to utilize similar activation signaling, e.g., BAFF, APRIL, CD40L, TLRs, IFN, IL6, IL10, and/or CXCL10 produced by APCs, to activate TACI and BCR for antibody production and immunoglobulin class switching [23].
For the present study, the overriding question is whether MZB cells function as the inducers of innate autoimmunity [53,54,55]. This would require a correct stromal cell environment within the MZ areas and, most likely, the presence of the immune-targeted autoantigen(s) to which extra-nodal lymphomagenesis can occur. Ontological development of MZs require the induced expression of ICAM1, VCAM1 and MADCAM that, in the presence of IL-7, recruits lymphoid tissue inducer cells (LTIs) that express RXRγt (retinoid X receptor gama) to sites surrounding blood vessels [56]. LTIs secrete lymphotoxin that, in turn, stimulates secretion of chemokines that attract dendritic and lymphoid cells, plus regulate endothelial cells delineating marginal sinuses [57]. Three additional stromal cell populations that dictate ontological effects on MZs are the marginal reticular cells (MRCs) which express RANKL (receptor activator of NF-kβ ligand), fibroblast reticular cells (FRCs) which along with MRCs regulate immune responses via production of reticular fibers and chemokines that guide lymphocyte trafficking, and follicular dendritic cells (FDCs) which secrete CXCL13 and BAFF to enhance B cell maturation, survival and emigration [58]. Importantly, these stromal cell populations secrete type 1 interferon (IFN1) in response to B cell secretion of lymphotoxin, giving rise to IFN-signatures that are common to the rheumatoid autoimmune diseases. As shown here, the salivary glands of C57BL/6.NOD-Aec1Aec2 mice exhibit upregulated Icam1, Vcam1 (the two LFA1 integrin units), Cxcl13 and chemokine receptors expressed by MZB cells, all at 16 weeks of age, the time frame when increasing numbers of lymphocytes are seen histologically entering the exocrine glands. As expected, Madcam1, not expressed in salivary glands, is not upregulated (data not shown).
The present analysis also reveals a distinct temporal profile for the Rho-GTPase families of proteins in the salivary glands of SSS C57BL/6.NOD-Aec1Aec2 mice as compared to SSNS C57BL/6J mice. Rho-GTPases play a critical role in signal transduction pathways involving multiple biological processes, especially cell migration and tissue homeostasis. Although we previously identified the Rac > Ras > Raf > Erk > AP gene pathway as a critically important component in development of glandular pathology [59], the current in-depth analyses of the Rho-GTPases points to a highly complex biological process indicated by the upregulation of genes encoding for multiple Rho-Gap, Rho-Gef, Dock and Cdc42 sub-family proteins of the Rho-GTPase family. Additional molecules include proteins involved in signal transduction pathways influencing cell functions and cell migrations, such as Was, Wasf and Rock [42,43].
Rho-GEF associated molecules are known to be critical for B cell functions. For example, Tedford et al. [60] reported that Vav1/Vav2-negative B cells are unresponsive to thymus-independent antigens in vivo and indicated a role for Vav-2 in BCR calcium signaling that is critical to B cell development and function. Similarly, Wang et al. [61] have reported that membrane-proximal BCR signaling molecules (including Vav3, Lyn, Syk, Btk, PLC-γ2, and Blnk), together with adaptor molecules Grb2, Cbl and Dok3, actually link BCR micro-clusters and motor proteins. Our transcriptome data indicate that Vav2 and Vav3, but not Vav1, exhibit upregulated activation profiles in the salivary glands of C57BL/6.NOD-Aec1Aec2 mice, while in C57BL/6J mice, Vav1 is strongly upregulated with Vav2 and Vav3 weakly upregulated. Sorting out the exact profiles for the various signal transduction pathways is the next step in better understanding these data.
In contrast, a direct correlation exists between upregulated gene expressions of Rho-GEF DOCK subfamily molecules encoding Dock2, 8, 10 and 11 with cellular homeostasis, development, and migration. Interactions of DOCK molecules with the Rho-family CDC42 proteins are also known to activate conformational changes in p21-activated PAK family molecules that regulate actin reorganization critical to cell adhesion and invasiveness, as well as to upregulate RAF and RAS pathways leading to activations of ERK, NF-kB, and AP1 (c-Jun, c-Fos and Atf2). As reported by Parrado and colleagues [62], expression of DOCK subgroup D molecules (i.e., Dock9, 10 and 11) occurs in peripheral blood lymphocytes; however, more recently, evidence has suggested that Dock10 expression is uniquely upregulated in B cells by IL-4. Furthermore, Dock10.1 isoform regulates T cell activities, while Dock10.2 isoform regulates CD23 expression on B cells, sustains B cell lymphopoiesis in secondary tissues, and up-regulates the IL4>Stat6 pathway [62], an essential pathway for SS development in C57BL/6.NOD-Aec1Aec2 mice [17]. This complexity has been expanded by Fukui et al., whose studies demonstrated that Dock2-Pi3k-delta activation regulates B cell migration and proliferation, while Dock2-Pi3k-gamma activation regulates T cell migration and proliferation [63,64]. Although the importance of Dock10 cannot be understated in B cell activation, Dock2 gene knockout mice are deficient in MZB cells [63], as are mice with deficiencies in Rac2, Ptk2 or Rho-GEFs.
Lastly, while Dock2, 10, and 11, as well as LFA1, Cxcr5, Cxcl13, and chemokine receptors Ccr6 and Ccr7 all show maximum expressions in the salivary glands of C57BL/6.NOD-Aec1Aec2 mice at around 16 weeks of age, Dock1, 5, 7 and 8 expressions occur earlier, at around 8 weeks of age. Based on histology, the expressions of Dock2, 10, and 11 are in line with the appearance of T and B cells involved with T cell-mediated cytotoxic responses, while expressions of Dock1, 5, 7 and 8 correspond with our predicted timeframe for physiological and structural changes within the exocrine glands and the apparent appearance of MZB cells. Furthermore, the involvement of Dock8 in B-cell development and function via activation of the BCR signaling molecules CD19 and BTK [65], together with LFA-1-dependent regulation of B and T cell positioning in germinal centers, suggest an environment that promotes IgG antibody responses to T-dependent antigens [66].

4. Materials and Methods

4.1. Animal Models

C57BL/6J and C57BL/6.NOD-Aec1Aec2 mice were bred and maintained under specific pathogen free (SPF) conditions within the Department of Pathology’s Mouse Facility with oversight by Animal Care Services at the University of Florida, Gainesville. The C57BL/6.NOD-Aec1Aec2 mouse is a well-characterized model of primary SS that spontaneously develops all major features of SS in human patients, except lymphomagenesis. The C57BL/6J mouse was used in generating the recombinant inbred C57BL/6.NOD-Aec1Aec2 mouse [67]. All animals were maintained on a 12-h light-dark schedule and provided food and acidified water ad libitum. Mice (n = 5 per experimental group) were euthanized at either 4, 8, 12, 16 or 20 weeks of age by cervical dislocation after deep anesthetization, the method approved by the Panel on Euthanasia of the American Veterinary Association. There were no indications that this procedure of euthanization affected the preparation of tissues or RNA. Both the breeding and use of these animals for the present studies were approved by the University of Florida Institutional Animal Care and Use Committee (IACUC).

4.2. Histology

Salivary glands surgically removed from each mouse at time of euthanasia were placed in 10% phosphate-buffered formalin for 24 h, embedded in paraffin and sectioned at 5 µm thickness. Following deparaffination and dehydration, each section was treated with blocking solution containing donkey serum. Sections were then stained with purified rat anti-mouse CD45R (Clone 30-F11, BD Pharmingen, San Jose, CA, USA) diluted 1:25 and goat polyclonal IgG anti-mouse CD3ε (Clone M-20, Santa Cruz Biotechnology, Santa Cruz, CA, USA) diluted 1:50 in an antibody diluent (Dako, Carpinteria, CA, USA) for 1 h at 25 °C. The slides were then washed with PBS followed by a 1 h incubation with Alexa Fluor 488 donkey anti-goat IgG (H+L) and Alexa Fluor 594 donkey anti-rat IgG (H+L) (Life Technologies, Grand Island, NY, USA). After a thorough wash with PBS, the slides were treated with a Vectashield DAPI-mounting medium (Vector Laboratory, Burlingame, CA, USA) and visualized microscopically at 200×.

4.3. RNA Preparations, Microarray Procedures and Microarray Data Analyses

Procedures for the isolation, preparation and quality testing of RNA samples, as well as analyses of microarray data, are described in full detail elsewhere [68,69]. In brief, salivary glands, free of lymph nodes, were excised from C57BL/6.NOD-Aec1Aec2 and C57BL/6J mice at 4-, 8-, 12-, 16- or 20 weeks of age (n = 5 per group), snap-frozen in liquid nitrogen and stored at −80 °C. Using one lobe from each explanted salivary gland, total RNA from each mouse was individually isolated using RNeasy Mini-Kits (Qiagen, Valencia, CA, USA), tested for purity and quantitated. Each RNA sample (n = 45 total) was hybridized on an Affymetrix 3′ Expression Array GeneChip Mouse Genome 430 2.0 array (Affymetrix, Santa Clara, CA, USA) and annotated (build 32; 06.09.2011). To verify array results, selected gene expressions were carried out using real-time polymerase chain reaction (rt-PCR) analyses. Microarray data were deposited with Expression Omnibus, accession number GSE15640.
Microarray data were normalized using the GCRMA (guanine-cytosine robust multi-array average) algorithm. Normalized data were analyzed using LIMMA (Linear Models for Microarray Analysis), available from the R Development Core Team [70], to identify differential gene expressions. The false discovery rate (FDR) method was used to adjust the p values to control for potential false values [71]. These temporal global data represent five equally spaced time points; thus, temporal patterns of gene expressions were used that included linear, quadratic, cubic, and quartic fit regression models. B-statistics were used to determine if a gene showed either positive or negative trends over time. To determine whether the observed number of gene counts exceeded the expected counts, one-tailed p values for enrichment of a particular biological process, molecular function, or pathway were calculated using the standard Fisher exact test. Pairwise comparisons between the five groups were made to assess each gene as differentially expressed at 8, 12, 16, and/or 20 weeks of age compared to their 4 weeks expressions. Furthermore, gene set profiling for specific biological processes and/or signal transduction pathways was carried out to match for identical temporal upregulated expressions as a means to support the concept that various individual genes comprise a functional gene set to help compensate for the fact that the encoded proteins may not be equally expressed and/or the RNA:aptamer binding on the arrays may be of different strengths for the represented genes.

5. Conclusions

In summary, the temporal transcriptome analyses presented in the current report provide new perspectives into several molecular mechanisms involved in lymphocyte emigrations, presumably of MZB cells, based on the unique time frame of their arrival in the salivary glands of our SSS C57BL/6.NOD-Aec1Aec2 mice (i.e., the earliest SS disease phase), plus the observed progressive pathophysiologic development of SS/SS-like disease thereafter. In addition, this study clearly shows the power of transcriptomic analyses to establish a foundation for further in-depth cellular investigations of SS-associated bioprocesses that may identify unique novel targets for possible intervention therapies. However, to better define and support the specific signal transduction pathways and their complementary interactive protein sets that culminate in the autoimmune attack against the salivary (and lacrimal) glands in Sjögren’s Syndrome, additional studies using single cells, single cell populations, and unique functional inhibitors must be carried out to confirm the data that are generated with transcriptome studies.

Author Contributions

Conceptualization, A.B.P., C.Q.N. and J.L.A.; methodology, A.B.P. and C.Q.N.; validation, A.B.P.; analyses, A.B.P.; investigation; A.B.P. and C.Q.N.; data curation, A.B.P.; writing, editing and visualization, A.B.P., C.Q.N. and J.L.A.; Funding and project management, A.B.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported in part by the National Institute of Dental and Craniofacial Research (NIDCR) grant R01 DE-014344 and the National Institute of Allergy and Infectious Diseases (NIAID) grant AI-0819529.

Institutional Review Board Statement

Both the breeding and use of these animals for the present studies were approved by the University of Florida’s Institutional Animal Care and Use Committee under protocols 2008-01-1756 and 2010-04-828.

Informed Consent Statement

Not applicable.

Data Availability Statement

Microarray data are publicly available and deposited with The Gene Expression Omnibus, Accession Numbers GSE15640 and GSE36378.

Conflicts of Interest

The authors declare no conflict of interest, and the NIH had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, nor in the decision to publish the results.

Abbreviations

APCAntigen-presenting cell
APRILTnfsf13 (Tumor necrosis factor superfamily, member 13
ARHGAPRho GTPase activating protein
ARHGEFRho GTPase exchange factor
BAFFTnfsf13b
BCRB cell receptor
CCRChemokine (C-C motif) receptor
CDCell differentiation antigen
CDC42Cell division cycle 42
CXCLChemokine (C-X-C motif) ligand
CXCRChemokine (C-X-C motif) receptor
DOCKDedicator of cytokinesis
ERKMitogen-activated protein kinase
FRCFibroblast reticular cell
GCRMAGuanine-cytosine robust multi-array average
GTP-GAPG protein bound-guanine nucleotide activating protein
GTP-GEFG protein bound-guanine nucleotide exchange factor
ICAM1Intercellular adhesion molecule 1
IFNInterferon
ILInterleukin
ITGIntegrin
LIMMALinear models for microarray analysis
LFLymphocytic foci
LFA1Lymphocyte function-association antigen 1
LTILymphoid tissue inducer cell
MADCAM1Mucosal cell adhesion molecule-1
MZMarginal zone
MZBMarginal zone B cell
NOTCHNotch gene homolog
PALSPeriarteriolar lymphoid sheaths
PTK2Protein tyrosine kinase
RACRas-related C3 botulinum toxin substrate
RANKLReceptor activator of NF-kB ligand
RBP-JRecombination signal binding protein for Igk J-region
RHO-GTPRas homolog member
RNARibonucleic acid
RNDRho-family GTPase
ROCKRho-associated coiled-coil containing protein
RSFRemodeling and spacing factor
RXRγRetinoid X receptor gamma
SSSjögren’s syndrome
S1PSphingosin-1-phosphate
TACITransmembrane activator and CAML interactor
TGTransgenic
TLRToll-like receptor
VCAM1Vascular cell adhesion molecule 1
VLA4Very late antigen 4
WASWiskott–Aldrich protein
WASFWiskott–Aldrich family member

References

  1. Vivino, F.B.; Bunya, V.Y.; Massaro-Giordano, G.; Johr, C.R.; Giattino, S.L.; Schorpion, A.; Shafer, B.; Peck, A.; Sivils, K.; Rasmussen, A.; et al. Sjogren’s Syndrome: An Update on Disease Pathogenesis, Clinical Manifestations and Treatment. Clin. Immunol. 2019, 203, 81–121. [Google Scholar] [CrossRef] [PubMed]
  2. Schmidt, R.E.; Grimbacher, B.; Witte, T. Autoimmunity and Primary Immunodeficiency: Two Sides of the Same Coin? Nat. Rev. Rheumatol. 2018, 14, 7–18. [Google Scholar] [CrossRef] [PubMed]
  3. Wu, H.; Chen, Y.; Zhu, H.; Zhao, M.; Lu, Q. The Pathogenic Role of Dysregulated Epigenetic Modifications in Autoimmune Diseases. Front. Immunol. 2019, 10, 2305. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Slight-Webb, S.; Bourn, R.L.; Holers, V.M.; James, J.A. Shared and Unique Immune Alterations in Pre-Clinical Autoimmunity. Curr. Opin. Immunol. 2019, 61, 60–68. [Google Scholar] [CrossRef]
  5. Teng, X.Y.; Cornaby, C.; Li, W.; Morel, L. Metabolic Regulation of Pathogenic Autoimmunity: Therapeutic Targeting. Curr. Opin. Immunol. 2019, 61, 10–16. [Google Scholar] [CrossRef]
  6. Marketos, N.; Cinoku, I.; Rapti, A.; Mavragani, C.P. Type I Interferon Signature in Sjogren’s Syndrome: Pathophysiological and Clinical Implications. Clin. Exp. Rheumatol. 2019, 37, S185–S191. [Google Scholar]
  7. Yin, H.; Wu, H.; Chen, Y.; Zhang, J.; Zheng, M.; Chen, G.; Li, L.; Lu, Q. The Therapeutic and Pathogenic Role of Autophagy in Autoimmune Diseases. Front. Immunol. 2018, 9, 1512. [Google Scholar] [CrossRef] [Green Version]
  8. Stojanov, S.; Kastner, D.L. Familial Autoinflammatory Diseases: Genetics, Pathogenesis and Treatment. Curr. Opin. Rheumatol. 2005, 17, 586–599. [Google Scholar] [CrossRef]
  9. Brito-Zeron, P.; Theander, E.; Baldini, C.; Seror, R.; Retamozo, S.; Quartuccio, L.; Bootsma, H.; Bowman, S.J.; Dorner, T.; Gottenberg, J.E.; et al. Early Diagnosis of Primary Sjogren’s Syndrome: EULAR-SS Task Force Clinical Recommendations. Expert Rev. Clin. Immunol. 2016, 12, 137–156. [Google Scholar] [CrossRef]
  10. Killedar, S.J.; Eckenrode, S.E.; McIndoe, R.A.; She, J.X.; Nguyen, C.Q.; Peck, A.B.; Cha, S. Early Pathogenic Events Associated with Sjogren’s Syndrome (Sjs)-Like Disease of the NOD Mouse Using Microarray Analysis. Lab. Investig. 2006, 86, 1243–1260. [Google Scholar] [CrossRef]
  11. Peck, A.B.; Nguyen, C.Q. Transcriptome Analysis of the Interferon-Signature Defining the Autoimmune Process of Sjogren’s Syndrome. Scand. J. Immunol. 2012, 76, 237–245. [Google Scholar] [CrossRef] [Green Version]
  12. Delaleu, N.; Nguyen, C.Q.; Tekle, K.M.; Jonsson, R.; Peck, A.B. Transcriptional Landscapes of Emerging Autoimmunity: Transient Aberrations in the Targeted Tissue’s Extracellular Milieu Precede Immune Responses in Sjogren’s Syndrome. Arthritis Res. Ther. 2013, 15, R174. [Google Scholar] [CrossRef] [Green Version]
  13. Nguyen, C.Q.; Cha, S.R.; Peck, A.B. Sjogren’s Syndrome (Sjs)-Like Disease of Mice: The Importance of B Lymphocytes and Autoantibodies. Front. Biosci. 2007, 12, 1767–1789. [Google Scholar] [CrossRef] [Green Version]
  14. Peck, A.B.; Nguyen, C.Q. What Can Sjogren’s Syndrome-Like Disease in Mice Contribute to Human Sjogren’s Syndrome? Clin. Immunol. 2017, 182, 14–23. [Google Scholar] [CrossRef] [Green Version]
  15. Robinson, C.P.; Yamachika, S.; Bounous, D.I.; Brayer, J.; Jonsson, R.; Holmdahl, R.; Peck, A.B.; Humphreys-Beher, M.G. A Novel NOD-Derived Murine Model of Primary Sjogren’s Syndrome. Arthritis Rheum. 1998, 41, 150–156. [Google Scholar] [CrossRef]
  16. Szczerba, B.M.; Rybakowska, P.D.; Dey, P.; Payerhin, K.M.; Peck, A.B.; Bagavant, H.; Deshmukh, U.S. Type I Interferon Receptor Deficiency Prevents Murine Sjogren’s Syndrome. J. Dent. Res. 2013, 92, 444–449. [Google Scholar] [CrossRef] [Green Version]
  17. Gao, J.H.; Killedar, S.; Cornelius, J.G.; Nguyen, C.; Cha, S.H.; Peck, A.B. Sjogren’s Syndrome in the NOD Mouse Model Is an Interleukin-4 Time-Dependent, Antibody Isotype-Specific Autoimmune Disease. J. Autoimmun. 2006, 26, 90–103. [Google Scholar] [CrossRef]
  18. Nguyen, C.; Cornelius, J.; Singson, E.; Killedar, S.; Cha, S.H.; Peck, A.B. Role of Complement and B Lymphocytes in Sjogren’s Syndrome-Like Autoimmune Exocrinopathy of NOD.B10-H2(B) Mice. Mol. Immunol. 2006, 43, 1332–1339. [Google Scholar] [CrossRef]
  19. Shen, L.; Zhang, C.; Wang, T.; Brooks, S.; Ford, R.J.; Lin-Lee, Y.C.; Kasianowicz, A.; Kumar, V.; Martin, L.; Liang, P. Development of Autoimmunity in IL-14 Alpha-Transgenic Mice. J. Immunol. 2006, 177, 5676–5686. [Google Scholar] [CrossRef] [Green Version]
  20. Shen, L.; Suresh, L.; Li, H.; Zhang, C.J.; Kumar, V.; Pankewycz, O.; Ambrus, J.L. IL-14 Alpha, the Nexus for Primary Sjogren’s Disease in Mice and Humans. Clin. Immunol. 2009, 130, 304–312. [Google Scholar] [CrossRef]
  21. Shen, L.; Suresh, L.; Wu, J.; Xuan, J.X.; Li, H.; Zhang, C.J.; Pankewycz, O.; Ambrus, J.L. A Role for Lymphotoxin in Primary Sjogren’s Disease. J. Immunol. 2010, 185, 6355–6363. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Shen, L.; Gao, C.; Suresh, L.; Xian, Z.; Song, N.; Chaves, L.D.; Yu, M.; Ambrus, J.L., Jr. Central Role for Marginal Zone B Cells in an Animal Model of Sjogren’s Syndrome. Clin. Immunol. 2016, 168, 30–36. [Google Scholar] [CrossRef] [Green Version]
  23. Cerutti, A.; Cols, M.; Puga, I. Marginal Zone B Cells: Virtues of Innate-Like Antibody-Producing Lymphocytes. Nat. Rev. Immunol. 2013, 13, 118–132. [Google Scholar] [CrossRef] [Green Version]
  24. Arruga, F.; Vaisitti, T.; Deaglio, S. The NOTCH Pathway and Its Mutations in Mature B Cell Malignancies. Front. Oncol. 2018, 8, 550. [Google Scholar] [CrossRef] [Green Version]
  25. Lu, T.T.; Cyster, J.G. Integrin-Mediated Long-Term B Cell Retention in the Splenic Marginal Zone. Science 2002, 297, 409–412. [Google Scholar] [CrossRef] [Green Version]
  26. Tan, J.B.; Xu, K.; Cretegny, K.; Visan, I.; Yuan, J.S.; Egan, S.E.; Guidos, C.J. Lunatic and Manic Fringe Cooperatively Marginal Zone B Cell Precursor Competition for Delta-like 1 in Splenic Endothelial Niches. Immunity 2009, 30, 254–263. [Google Scholar] [CrossRef] [Green Version]
  27. Baker, N.; Ehrenstein, M.R. Cutting Edge: Selection of B Lymphocyte Subsets Is Regulated by Natural IgM. J. Immunol. 2002, 169, 6686–6690. [Google Scholar] [CrossRef] [Green Version]
  28. Lawson, C.D.; Ridley, A.J. Rho GTPase Signaling Complexes in Cell Migration and Invasion. J. Cell Biol. 2018, 217, 447–457. [Google Scholar] [CrossRef]
  29. Kroese, F.G.M.; Abdulahad, W.H.; Haacke, E.; Bos, N.A.; Vissink, A.; Bootsma, H. B-cell Hyperactivity in Primary Sjogren’s Syndrome. Expert Rev. Clin. Immunol. 2014, 10, 483–499. [Google Scholar] [CrossRef]
  30. Karlsen, M.; Jakobsen, K.; Jonsson, R.; Hammenfors, D.; Hansen, T.; Appel, S. Expression of Toll-Like Receptors in Peripheral Blood Mononuclear Cells of Patients with Primary Sjogren’s Syndrome. Scand. J. Immunol. 2017, 85, 220–226. [Google Scholar] [CrossRef]
  31. Ambrus, J.L.; Suresh, L.; Peck, A. Multiple Roles for B-Lymphocytes in Sjogren’s Syndrome. J. Clin. Med. 2016, 5, 87. [Google Scholar] [CrossRef]
  32. Ambrus, J.L.; Nguyen, C.Q.; Peck, A.B. Early Covert Appearance of Marginal Zone B Cells in Sjögren’s Syndrome-Susceptible Mice: Initiators of Subsequent Overt Clinical Disease. Int. J. Mol. Sci. 2021, 22, 1919. [Google Scholar]
  33. Sozzani, S.; Bosisio, D.; Scarsi, M.; Tincani, A. Type I Interferons in Systemic Autoimmunity. Autoimmunity 2010, 43, 196–203. [Google Scholar] [CrossRef]
  34. Nguyen, C.Q.; Peck, A.B. The Interferon-Signature of Sjögren´S Signature: How Unique Biomarkers Can Identify Underlying Inflammatory and Immunopathological Mechanisms of Specific Disease. Front. Immunol. 2013, 4, 142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Maillard, I.; Weng, A.P.; Carpenter, A.C.; Rodriguez, C.G.; Sai, H.; Xu, L.; Allman, D.; Aster, J.C.; Pear, W.S. Mastermind Critically Regulates Notch-Mediated Lymphoid Cell Fate Decisions. Blood 2004, 104, 1696–1702. [Google Scholar] [CrossRef] [Green Version]
  36. Saito, T.; Chiba, S.; Ichikawa, M.; Kunisato, A.; Asai, T.; Shimizu, K.; Yamaguchi, T.; Yamamoto, G.; Seo, S.; Kumano, K.; et al. Notch2 Is Preferentially Expressed in Mature B Cells and Indispensable for Marginal Zone B Lineage Development. Immunity 2003, 18, 675–685. [Google Scholar] [CrossRef] [Green Version]
  37. LopesCarvalho, T.; Foote, J.; Kearney, J.F. Marginal Zone B Cells in Lymphocyte Activation and Regulation. Curr. Opin. Immunol. 2005, 17, 244–250. [Google Scholar] [CrossRef] [PubMed]
  38. Tybulewicz, V.L.; Henderson, R.B. Rho Family Gtpases and Their Regulators in Lymphocytes. Nat. Rev. Immunol. 2009, 9, 630–644. [Google Scholar] [CrossRef]
  39. Gotoh, K.; Tanaka, Y.; Nishikimi, A.; Inayoshi, A.; Enjoji, M.; Takayanagi, R.; Sasazuki, T.; Fukui, Y. Differential Requirement for DOCK2 in Migration of Plasmacytoid Dendritic Cells versus Myeloid Dendritic Cells. Blood 2008, 111, 2973–2976. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Namekata, K.; Kimura, A.; Kawamura, K.; Harada, C.; Harada, T. Dock GEFs and Their Therapeutic Potential: Neuroprotection and Axon Regeneration. Prog. Retin. Eye Res. 2014, 43, 1–16. [Google Scholar] [CrossRef] [PubMed]
  41. Gadea, G.; Blangy, A. Dock-Family Exchange Factors in Cell Migration and Disease. Eur. J. Cell Biol. 2014, 93, 466–477. [Google Scholar] [CrossRef]
  42. Heasman, S.J.; Ridley, A.J. Mammalian Rho GTPases: New Insights into Their Functions from In Vivo Studies. Nat. Rev. Mol. Cell Biol. 2008, 9, 690–701. [Google Scholar] [CrossRef] [PubMed]
  43. Spiering, D.; Hodgson, L. Dynamics of the Rho-Family Small Gtpases in Actin Regulation and Motility. Cell Adh. Migr. 2011, 5, 170–180. [Google Scholar] [CrossRef] [Green Version]
  44. Casola, S.; Otipoby, K.L.; Alimzhanov, M.; Humme, S.; Uyttersprot, N.; Kutok, J.L.; Carroll, M.C.; Rajewsky, K. B Cell Receptor Signal Strength Determines B Cell Fate. Nat. Immunol. 2004, 5, 317–327. [Google Scholar] [CrossRef] [PubMed]
  45. Pillai, S.; Cariappa, A.; Moran, S.T. Marginal Zone B Cells. Annu. Rev. Immunol. 2005, 23, 161–196. [Google Scholar] [CrossRef]
  46. Oliver, A.M.; Martin, F.; Kearney, J.F. IgM(High)CD21(High) lymphocytes Enriched in the Splenic Marginal Zone Generate Effector Cells More Rapidly than the Bulk of Follicular B Cells. J. Immunol. 1999, 162, 7198–7207. [Google Scholar] [PubMed]
  47. Lenert, P.; Brummel, R.; Field, E.H.; Ashman, R.F. TLR-9 Activation of Marginal Zone B Cells in Lupus Mice Regulates Immunity through Increased IL-10 Production. J. Clin. Immunol. 2005, 25, 29–40. [Google Scholar] [CrossRef] [PubMed]
  48. Treml, L.S.; Carlesso, G.; Hoek, K.L.; Stadanlick, J.E.; Kambayashi, T.; Bram, R.J.; Cancro, M.P.; Khan, W.N. TLR Stimulation Modifies Blys Receptor Expression in Follicular and Marginal Zone B Cells. J. Immunol. 2007, 178, 7531–7539. [Google Scholar] [CrossRef] [Green Version]
  49. Weill, J.C.; Weller, S.; Reynaud, C.A. Human Marginal Zone B Cells. Annu. Rev. Immunol. 2009, 27, 267–285. [Google Scholar] [CrossRef]
  50. Attanavanich, K.; Kearney, J.F. Marginal Zone, but Not Follicular B Cells, Are Potent Activators of Naive CD4 T Cells. J. Immunol. 2004, 172, 803–811. [Google Scholar] [CrossRef] [Green Version]
  51. Jamin, C.; Morva, A.; Lemoine, S.; Daridon, C.; de Mendoza, A.R.; Youinou, P. Regulatory B Lymphocytes in Humans—A Potential Role in Autoimmunity. Arthritis Rheum. 2008, 58, 1900–1906. [Google Scholar] [CrossRef] [PubMed]
  52. LopesCarvalho, T.; Kearney, J.F. Development and Selection of Marginal Zone B Cells. Immunol. Rev. 2004, 197, 192–205. [Google Scholar] [CrossRef]
  53. Allam, R.; Anders, H.J. The Role of Innate Immunity in Autoimmune Tissue Injury. Curr. Opin. Rheumatol. 2008, 20, 538–544. [Google Scholar] [CrossRef]
  54. Theofilopoulos, A.N.; Kono, D.H.; Baccala, R. The Multiple Pathways to Autoimmunity. Nat. Immunol. 2017, 18, 716–724. [Google Scholar] [CrossRef] [PubMed]
  55. Weidenbusch, M.; Kulkarni, O.P.; Anders, H.J. The Innate Immune System in Human Systemic Lupus Erythematosus. Clin. Sci. 2017, 131, 625–634. [Google Scholar] [CrossRef] [PubMed]
  56. Oliver, A.M.; Martin, F.; Gartland, G.L.; Carter, R.H.; Kearney, J.F. Marginal Zone B Cells Exhibit Unique Activation, Proliferative and Immunoglobulin Secretory Responses. Eur. J. Immunol. 1997, 27, 2366–2374. [Google Scholar] [CrossRef]
  57. Colonna, M. Innate Lymphoid Cells: Diversity, Plasticity, and Unique Functions in Immunity. Immunity 2018, 48, 1104–1117. [Google Scholar] [CrossRef] [Green Version]
  58. Magri, G.; Miyajima, M.; Bascones, S.; Mortha, A.; Puga, I.; Cassis, L.; Barra, C.M.; Comerma, L.; Chudnovskiy, A.; Gentile, M.; et al. Innate Lymphoid Cells Integrate Stromal and Immunological Signals to Enhance Antibody Production by Splenic Marginal Zone B Cells. Nat. Immunol. 2014, 15, 354–364. [Google Scholar] [CrossRef] [Green Version]
  59. Nguyen, C.Q.; Cornelius, J.G.; Cooper, L.; Neff, J.; Tao, J.; Lee, B.H.; Peck, A.B. Identification of Possible Candidate Genes Regulating Sjogren’s Syndrome-Associated Autoimmunity: A Potential Role for TNFSF4 in Autoimmune Exocrinopathy. Arthritis Res. Ther. 2008, 10, R137. [Google Scholar] [CrossRef] [Green Version]
  60. Tedford, K.; Nitschke, L.; Girkontaite, I.; Charlesworth, A.; Chan, G.; Sakk, V.; Barbacid, M.; Fischer, K.D. Compensation between Vav-1 and Vav-2 in B cell Development and Antigen Receptor Signaling. Nat. Immunol. 2001, 2, 548–555. [Google Scholar] [CrossRef]
  61. Wang, J.; Lin, F.; Wan, Z.; Sun, X.; Lu, Y.; Huang, J.; Wang, F.; Zeng, Y.; Chen, Y.H.; Shi, Y.; et al. Profiling the Origin, Dynamics, and Function of Traction Force in B Cell Activation. Sci. Signal. 2018, 11, eaai9192. [Google Scholar] [CrossRef] [Green Version]
  62. Garcia-Serna, A.M.; Alcaraz-Garcia, M.J.; Ruiz-Lafuente, N.; Sebastian-Ruiz, S.; Martinez, C.M.; Moya-Quiles, M.R.; Minguela, A.; Garcia-Alonso, A.M.; Martin-Orozco, E.; Parrado, A. Dock10 Regulates CD23 Expression and Sustains B-Cell Lymphopoiesis in Secondary Lymphoid Tissue. Immunobiology 2016, 221, 1343–1350. [Google Scholar] [CrossRef]
  63. NombelaArrieta, C.; Lacalle, R.A.; Montoya, M.C.; Kunisaki, Y.; Megias, D.; Marques, M.; Carrera, A.C.; Manes, S.; Fukui, Y.; MartinezA, C.; et al. Differential Requirements for DOCK2 and Phosphoinositide-3-Kinase Gamma during T and B Lymphocyte Homing. Immunity 2004, 21, 429–441. [Google Scholar] [CrossRef] [Green Version]
  64. Ushijima, M.; Uruno, T.; Nishikimi, A.; Sanematsu, F.; Kamikaseda, Y.; Kunimura, K.; Sakata, D.; Okada, T.; Fukui, Y. The Rac Activator DOCK2 Mediates Plasma Cell Differentiation and IgG Antibody Production. Front. Immunol. 2018, 9, 243. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Ma, C.S.; Tangye, S.G. Flow Cytometric-Based Analysis of Defects in Lymphocyte Differentiation and Function Due to Inborn Errors of Immunity. Front. Immunol. 2019, 10, 2108. [Google Scholar] [CrossRef]
  66. Randall, K.L.; Lambe, T.; Johnson, A.; Treanor, B.; Kucharska, E.; Domaschenz, H.; Whittle, B.; Tze, L.E.; Enders, A.; Crockford, T.L.; et al. Dock8 Mutations Cripple B Cell Immunological Synapses, Germinal Centers and Long-Lived Antibody Production. Nat. Immunol. 2009, 10, 1283–1291. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Cha, S.; Nagashima, H.; Brown, V.B.; Peck, A.B.; Humphreys-Beher, M.G. Two NOD Idd-Associated Intervals Contribute Synergistically to the Development of Autoimmune Exocrinopathy (Sjogren’s Syndrome) on a Healthy Murine Background. Arthritis Rheum. 2002, 46, 1390–1398. [Google Scholar] [CrossRef]
  68. Nguyen, C.Q.; Sharma, A.; She, J.X.; McIndoe, R.A.; Peck, A.B. Differential Gene Expressions in the Lacrimal Gland During Development and Onset of Keratoconjunctivitis Sicca in Sjogren’s Syndrome (SJS)-Like Disease of the C57BL/6.NOD-Aec1Aec2 Mouse. Exp. Eye Res. 2009, 88, 398–409. [Google Scholar] [CrossRef] [Green Version]
  69. Peck, A.B.; Saylor, B.T.; Nguyen, L.; Sharma, A.; She, J.X.; Nguyen, C.Q.; McIndoe, R.A. Gene Expression Profiling of Early-Phase Sjogren’s Syndrome in C57BL/6.NOD-Aec1Aec2 Mice Identifies Focal Adhesion Maturation Associated with Infiltrating Leukocytes. Investig. Ophthalmol. Vis. Sci. 2011, 52, 5647–5655. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. The R Project for Statistical Computing. Available online: www.r-project.org (accessed on 17 June 2021).
  71. Benjanini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R Statist. Soc. B 1995, 57, 289–300. [Google Scholar] [CrossRef]
Figure 1. Histological photomicrographs depicting lymphocytic infiltrations present within the salivary glands of C57BL/6.NOD-Aec1Aec2 mice at 16 (A) and 31 (B) weeks of age. Infiltrations are age-dependent and contain an ever-changing proportion of both T cells (green fluorescence) and B cells (red fluorescence). Lymphocytic infiltration, especially in SS-susceptible C57BL/6.NOD-Aec1Aec2 mice, is strongly periductal. Visualization is at 200X magnification.
Figure 1. Histological photomicrographs depicting lymphocytic infiltrations present within the salivary glands of C57BL/6.NOD-Aec1Aec2 mice at 16 (A) and 31 (B) weeks of age. Infiltrations are age-dependent and contain an ever-changing proportion of both T cells (green fluorescence) and B cells (red fluorescence). Lymphocytic infiltration, especially in SS-susceptible C57BL/6.NOD-Aec1Aec2 mice, is strongly periductal. Visualization is at 200X magnification.
Ijms 22 07176 g001
Figure 2. Coordinated temporal activation of factors important for MZB cell emigration from splenic marginal zone (MZ) areas to the salivary glands of C57BL/6.NOD-Aec1Aec2 mice at onset of SS-like disease. Transcriptomic data reveal that genes encoding ICAM1 and VCAM1 (Receptors for LFA1), LFA-1 (ItgaL and Itgβ2), Cxcl13, Cxcr5 (the Cxcl13 receptor), Dock2, Ccr6, Ccr7, and Ptk2β are all upregulated in the salivary glands of C57BL/6.NOD-Aec1Aec2 mice at 16 weeks of age (right panel), the approximate time point when the immune response is transitioning to the adaptive immune response phase. In contrast, no similar upregulated expressions of these genes are seen in the salivary glands of C57BL/6J mice (left panel).
Figure 2. Coordinated temporal activation of factors important for MZB cell emigration from splenic marginal zone (MZ) areas to the salivary glands of C57BL/6.NOD-Aec1Aec2 mice at onset of SS-like disease. Transcriptomic data reveal that genes encoding ICAM1 and VCAM1 (Receptors for LFA1), LFA-1 (ItgaL and Itgβ2), Cxcl13, Cxcr5 (the Cxcl13 receptor), Dock2, Ccr6, Ccr7, and Ptk2β are all upregulated in the salivary glands of C57BL/6.NOD-Aec1Aec2 mice at 16 weeks of age (right panel), the approximate time point when the immune response is transitioning to the adaptive immune response phase. In contrast, no similar upregulated expressions of these genes are seen in the salivary glands of C57BL/6J mice (left panel).
Ijms 22 07176 g002
Figure 3. Unique upregulation profile of genes encoding for members of the Rho-GTPase family of proteins in the salivary glands of C57BL/6.NOD-Aec1Aec2 mice. Rho-GTPase family proteins are important in regulating the cellular homeostasis of the GTP<>GDP system. The Rho-GTPase family consists of 21 members, of which 16 are divided among 5 subfamilies (Rsf1, Rsf2, Rsf3, Rsf4 and Rsf5), with 5 unassigned (designated Rsfx here). For cells in the non-activated state, Rho-GTP molecules are mostly associated with cellular membranes, but during cell activation the molecules dissociate to the cytoplasm by an as-yet unknown mechanism. In the salivary glands of C57BL/6.NOD-Aec1Aec2 mice, genes encoding for Rhoc, Rhou, Cdc42, Rac3, and Rnd1 are upregulated starting at 8 weeks of age, with Rac2 and Rhobtb1 upregulated starting at 16 weeks of age (lower panel). In contrast, the salivary glands of C57BL/6J mice showed a weak upregulated expression for Rnd1, Rnd2, Rhobtb3, Rhoh, Rhof and Rhov at 8 weeks of age, but this was not prolonged.
Figure 3. Unique upregulation profile of genes encoding for members of the Rho-GTPase family of proteins in the salivary glands of C57BL/6.NOD-Aec1Aec2 mice. Rho-GTPase family proteins are important in regulating the cellular homeostasis of the GTP<>GDP system. The Rho-GTPase family consists of 21 members, of which 16 are divided among 5 subfamilies (Rsf1, Rsf2, Rsf3, Rsf4 and Rsf5), with 5 unassigned (designated Rsfx here). For cells in the non-activated state, Rho-GTP molecules are mostly associated with cellular membranes, but during cell activation the molecules dissociate to the cytoplasm by an as-yet unknown mechanism. In the salivary glands of C57BL/6.NOD-Aec1Aec2 mice, genes encoding for Rhoc, Rhou, Cdc42, Rac3, and Rnd1 are upregulated starting at 8 weeks of age, with Rac2 and Rhobtb1 upregulated starting at 16 weeks of age (lower panel). In contrast, the salivary glands of C57BL/6J mice showed a weak upregulated expression for Rnd1, Rnd2, Rhobtb3, Rhoh, Rhof and Rhov at 8 weeks of age, but this was not prolonged.
Ijms 22 07176 g003
Figure 4. Unique upregulation profiles for genes encoding for members of the Rho-GTPase GAP and GEF families of proteins. Transcriptome data showing the temporal expressions of 8 of the 23 Arhgap genes and 7 of the 17 Arhgef genes upregulated in the salivary glands of C57BL/6.NOD-Aec1Aec2 mice (lower panels). In each case, upregulated expressions started around 8 weeks of age. In contrast, expressions of these Arhgap and Arhgef genes remain downregulated.
Figure 4. Unique upregulation profiles for genes encoding for members of the Rho-GTPase GAP and GEF families of proteins. Transcriptome data showing the temporal expressions of 8 of the 23 Arhgap genes and 7 of the 17 Arhgef genes upregulated in the salivary glands of C57BL/6.NOD-Aec1Aec2 mice (lower panels). In each case, upregulated expressions started around 8 weeks of age. In contrast, expressions of these Arhgap and Arhgef genes remain downregulated.
Ijms 22 07176 g004
Figure 5. Upregulated profiles of genes encoding for members of the Rho-GTPase Dock family of proteins. Transcriptome data showing the temporal expressions of the Dock genes upregulated in the salivary glands of C57BL/6.NOD-Aec1Aec2 mice (lower panel) versus C57BL/6J mice (upper panel). In the SS-susceptible C57BL/6.NOD-Aec1Aec2 mice three genes, Dock2, Dock10, and Dock11, reported to be uniquely associated with B and T lymphocyte functions, exhibit a concomitant upregulated short-term expression at 16 weeks of age (arrows), while Dock1, 5, 7 and 8 are each upregulated starting at 4 weeks of age and remain activated out to 16 or 20 weeks of age. In contrast, in the SS-non-susceptible C57BL/6J mice, Dock4, 7, 9, 10 and 11 exhibit an upregulated expression starting at 4 weeks of age, while Dock8 shows a unique expression pattern being upregulated at 16 weeks of age.
Figure 5. Upregulated profiles of genes encoding for members of the Rho-GTPase Dock family of proteins. Transcriptome data showing the temporal expressions of the Dock genes upregulated in the salivary glands of C57BL/6.NOD-Aec1Aec2 mice (lower panel) versus C57BL/6J mice (upper panel). In the SS-susceptible C57BL/6.NOD-Aec1Aec2 mice three genes, Dock2, Dock10, and Dock11, reported to be uniquely associated with B and T lymphocyte functions, exhibit a concomitant upregulated short-term expression at 16 weeks of age (arrows), while Dock1, 5, 7 and 8 are each upregulated starting at 4 weeks of age and remain activated out to 16 or 20 weeks of age. In contrast, in the SS-non-susceptible C57BL/6J mice, Dock4, 7, 9, 10 and 11 exhibit an upregulated expression starting at 4 weeks of age, while Dock8 shows a unique expression pattern being upregulated at 16 weeks of age.
Ijms 22 07176 g005
Figure 6. Rho-Arhgef pathway activations. Rho-Arhgef molecules are associated with several signal transduction pathways including Rac-GTP, Rho-GTP and CDC42-GP, each leading to downstream activation of a biological process important for cellular migrations. These three pathways interact, respectively, with Wasf, Rock, and Was protein family members. In the salivary glands of C57BL/6.NOD-Aec1Aec2 mouse these appear to be Wasf3, Rock1, Rock2 and Was (right panel). The genes encoding Wasf3, Rock1 and Rock2 show an activation starting at 4 weeks of age, while the gene encoding Was is upregulated uniquely at 16 weeks of age. In the salivary glands of C57BL/6J mice, a strong activation is observed with Wasf3, and possibly a weak response by Rock1.
Figure 6. Rho-Arhgef pathway activations. Rho-Arhgef molecules are associated with several signal transduction pathways including Rac-GTP, Rho-GTP and CDC42-GP, each leading to downstream activation of a biological process important for cellular migrations. These three pathways interact, respectively, with Wasf, Rock, and Was protein family members. In the salivary glands of C57BL/6.NOD-Aec1Aec2 mouse these appear to be Wasf3, Rock1, Rock2 and Was (right panel). The genes encoding Wasf3, Rock1 and Rock2 show an activation starting at 4 weeks of age, while the gene encoding Was is upregulated uniquely at 16 weeks of age. In the salivary glands of C57BL/6J mice, a strong activation is observed with Wasf3, and possibly a weak response by Rock1.
Ijms 22 07176 g006
Figure 7. Temporal and differential gene expressions of various factors involved with signal transduction pathways regulated by members of the Rho-GTPase families of proteins. Transcriptome data of multiple genes (i.e., Cdc42, Vav, and Baiap (brain-specific angiogenesis inhibitor -1 associated protein) family members, plus Scrib1, Setd6, Grb2, Flii, Creb3 and Foxo3a) exhibit upregulated expressions in the salivary glands of C57BL/6.NOD-Aec1Aec2 mice, virtually all starting at 8 weeks of age (lower panel). C57BL/6J mice exhibit a general lack of upregulated expression, except for Vav1 and possibly Cdc42bpg and Vav2 (upper panel).
Figure 7. Temporal and differential gene expressions of various factors involved with signal transduction pathways regulated by members of the Rho-GTPase families of proteins. Transcriptome data of multiple genes (i.e., Cdc42, Vav, and Baiap (brain-specific angiogenesis inhibitor -1 associated protein) family members, plus Scrib1, Setd6, Grb2, Flii, Creb3 and Foxo3a) exhibit upregulated expressions in the salivary glands of C57BL/6.NOD-Aec1Aec2 mice, virtually all starting at 8 weeks of age (lower panel). C57BL/6J mice exhibit a general lack of upregulated expression, except for Vav1 and possibly Cdc42bpg and Vav2 (upper panel).
Ijms 22 07176 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Peck, A.B.; Nguyen, C.Q.; Ambrus, J.L. Upregulated Chemokine and Rho-GTPase Genes Define Immune Cell Emigration into Salivary Glands of Sjögren’s Syndrome-Susceptible C57BL/6.NOD-Aec1Aec2 Mice. Int. J. Mol. Sci. 2021, 22, 7176. https://doi.org/10.3390/ijms22137176

AMA Style

Peck AB, Nguyen CQ, Ambrus JL. Upregulated Chemokine and Rho-GTPase Genes Define Immune Cell Emigration into Salivary Glands of Sjögren’s Syndrome-Susceptible C57BL/6.NOD-Aec1Aec2 Mice. International Journal of Molecular Sciences. 2021; 22(13):7176. https://doi.org/10.3390/ijms22137176

Chicago/Turabian Style

Peck, Ammon B., Cuong Q. Nguyen, and Julian L. Ambrus. 2021. "Upregulated Chemokine and Rho-GTPase Genes Define Immune Cell Emigration into Salivary Glands of Sjögren’s Syndrome-Susceptible C57BL/6.NOD-Aec1Aec2 Mice" International Journal of Molecular Sciences 22, no. 13: 7176. https://doi.org/10.3390/ijms22137176

APA Style

Peck, A. B., Nguyen, C. Q., & Ambrus, J. L. (2021). Upregulated Chemokine and Rho-GTPase Genes Define Immune Cell Emigration into Salivary Glands of Sjögren’s Syndrome-Susceptible C57BL/6.NOD-Aec1Aec2 Mice. International Journal of Molecular Sciences, 22(13), 7176. https://doi.org/10.3390/ijms22137176

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop