- 1Cleveland Clinic Lerner College of Medicine, Department of Molecular Medicine, Case Western Reserve University, Cleveland, OH, United States
- 2Department of Inflammation and Immunity, Lerner Research Institute, NE40, Cleveland Clinic Foundation, Cleveland, OH, United States
- 3Cleveland Clinic Lerner College of Medicine at Case Western Reserve University, Cleveland, OH, United States
- 4Institute of Immunology, University of Muenster, Muenster, Germany
Systemic lupus erythematosus (SLE) is an autoimmune disorder disproportionally affecting women. A similar sex difference exists in the murine New Zealand Black/White hybrid model (NZBWF1) of SLE with all females, but only 30-40% of males, developing disease within the first year of life. Myeloid-derived suppressor cells (MDSCs) are prominent in NZBWF1 males and while depletion of these cells in males, but not females, promotes disease development, the mechanism of suppression remains unknown. S100a9, expressed by neutrophils and MDSCs, has previously been shown to exert immunosuppressive functions in cancer and inflammation. Here we investigated if S100a9 exerts immunosuppressive functions in NZBWF1 male and female mice. S100a9+/+, S100a9+/- and S100a9-/- NZBWF1 mice were followed for disease development for up to 8 months of age. Serum autoantibody levels, splenomegaly, lymphocyte activation, glomerulonephritis and proteinuria were measured longitudinally or at the time of harvest. In accordance with an immunosuppressive function of MDSCs in male mice, S100a9-deficient male NZBWF1 mice developed accelerated autoimmunity as indicated by increased numbers of differentiated effector B and T cells, elevated serum autoantibody levels, increased immune-complex deposition and renal inflammation, and accelerated development of proteinuria. In contrast, female mice showed either no response to S100a9-deficiency or even a slight reduction in disease symptoms. Furthermore, male, but not female, S100a9-/- NZBWF1 mice displayed an elevated type I interferon-induced gene signature, suggesting that S100a9 may dampen a pathogenic type I interferon signal in male mice. Taken together, S100a9 exerts an immunosuppressive function in male NZBWF1 mice effectively moderating lupus-like disease development via inhibition of type I interferon production, lymphocyte activation, autoantibody production and the development of renal disease.
Introduction
Systemic lupus erythematosus (SLE) is a chronic autoimmune inflammatory disease presenting with variable manifestations (1). Common indicators of SLE are elevated serum antinuclear autoantibodies (ANAs) and immune cell infiltration of affected organs such as the skin, heart, joints, or kidneys. The disease presents with a strong female predominance especially during child-bearing years, which has prompted much interest into the role of sex hormones (2–5). As such, numerous animal studies have shown that androgens exert a protective effect, while estrogens exacerbate disease development, the latter likely via direct effects on lymphocyte development and tolerance (6–12). Few studies however have looked into the mechanism driving immunoprotective effects of testosterone [reviewed in (13)]. Interestingly, immunosuppressive cells such as regulatory T and B cells, M2 macrophages and myeloid-derived suppressor cells (MDSCs) appear to be regulated by testosterone either directly during development or functionally during an immune response (14–20), suggesting that the presence of testosterone could lead to dysregulation of these cell populations and thus immune protection.
MDSCs have been ascribed a role in multiple immune disorders, but are best known for their role in cancer (21–25). The cells have a remarkable ability to suppress T cell responses by several methods including the production of reactive oxygen- and nitrogen-species (ROS, NOS), release of immunosuppressive cytokines (e.g. TGFβ, IL-10), production of indolamine 2’,3’-dioxygenase (IDO), and the production of S100A8/S100A9 heterodimers (also known as calprotectin) (26–29). While most of these effector mechanisms are shared with many other cell subsets, S100A8 and S100A9 proteins are exclusively expressed by myeloid cells constituting up to 45% of all cytosolic proteins in human neutrophils (30–32). Levels of heterodimers of S100A8 and S100A9 (S100A8/A9) are often elevated during inflammation, although whether the complex exerts immune stimulatory or immunoprotective properties depends on the underlying pathology (24, 25, 33–39). For example, S100a9-/- mice showed increased susceptibility to respiratory infection and increased renal fibrosis and damage supporting an immunoprotective or immune-resolution role (33, 34). Oppositely, S100a9-deficiency promoted reduced disease in animal models of arthritis and Alzheimer’s disease (24, 36). Importantly, none of these studies compared the response to S100a9-deficiency between males and females. In SLE patients (97% female cohort), serum levels of S100a9 have been found to be elevated and associated with SLEDAI scores and active bacterial infections (40). As for most studies, it remains unknown whether elevated S100a9 levels reflect proinflammatory or inflammation-resolution properties.
The (New Zealand black × New Zealand white)F1 (NZBWF1) mouse model develops a lupus-like disease characterized by hyperactive B cells, abnormal autoantibody production, glomerulonephritis, IgG immune complex (IC) deposition in the kidney glomeruli with complement fixation, and eventual renal failure resembling human SLE (41). Importantly, the model displays a similar female predominance with 100% of female NZBWF1 mice, but only 30-40% of male NZBWF1 mice, developing disease within one year. We and others previously reported an immunoprotective role for MDSC-like cells in male and young, prepubescent female NZBWF1 mice (14, 15, 42). While ROS/NOS has been identified as a mechanism of suppression in young females (14), the mechanism used by MDSCs in male NZBWF1 mice has not been identified. Based on studies suggesting a role for S100a9 in MDSC-driven immunosuppression in cancer (25), we hypothesized that MDSCs utilized S100a9 to inhibit immune activation and disease progression in NZBWF1 lupus-prone males. We report here that MDSC-like cells from S100a9-/- NZBWF1 mice failed to suppress B cell differentiation in vitro, suggesting that S100a9 is required for the immunosuppressive function of MDSC-like cells in these mice. Furthermore, male, but not female, S100a9-/- NZBWF1 mice showed increased disease development including elevated splenomegaly, hyperactive B and T cells, and accelerated renal disease. Interestingly, disease development in S100a9-deficient male mice was associated with elevated levels of type I interferon-stimulated gene transcripts and accumulation of IFNα-producing low density granulocytes (LDGs), suggesting that S100a9 limits type I interferon production in male lupus-prone mice hereby protecting the mice from disease development.
Materials and Methods
Animals and Cells
S100a9-deficient C57Bl/6 mice were obtained from Dr. Thomas Vogl (Westfälische Wilhelms-Universität Münster). S100a9-deficiency was backcrossed onto the NZB/BinJ and NZW/LacZ backgrounds for at least 8 generations. Crosses between NZB.S100a9+/- and NZW.S100a9+/- mice were performed to generate S100a9+/+, S100a9+/- and S100a9-/- NZBWF1 mice. For all analyses, littermates were used for comparison. Proteinuria was obtained from mice monthly starting at 3-4 months of age using dipsticks (Roche Biotech). Readings are presented on a scale from 0-4 (0, trace, 1, 2, 3, 4) corresponding with albumin levels of 0-400mg/dL urine as described by the manufacturer (Roche). When a mouse showed levels at or above 250 mg/dL (≥ 3), it was retested one week later. If consistently high, the mouse was euthanized for immediate analysis. Mice that did not develop proteinuria were euthanized at 8 months of age. For immunizations, S100a9+/+, S100a9+/- and S100a9-/- male and female NZBWF1 mice were immunized i.p. with 20µg (4-Hydroxy-3-nitrophenylacetyl)27 conjugated chicken γ-globulin (NP27-CGG) in complete Freund’s adjuvant (day 0). Immunized mice were bled on days -1, 7, 14, 21 and 28, and all animals were euthanized thereafter. Mice were maintained in the Biological Research Unit at the Lerner Research Institute, in accordance with Cleveland Clinic Foundation Animal Research Committee guidelines and all procedures were approved by the Institutional Animal Care and Use Committee of the Lerner Research Institute of the Cleveland Clinic Foundation and conducted in compliance with guidelines issued by the National Institutes of Health. For in vitro studies Gr1highCD11b+, Gr1lowCD11b+, B220highCD19+CD138- cells (in vitro differentiation assay) and Ly6ChighCD11b+SSClow, Ly6ClowCD11b+SSClow cells (real-time RT-PCR) were isolated by fluorescence activated cell sorting (FACS) on a FACSAria I (BD Biosciences, San Jose, CA). All sorted cells were confirmed negative for CD3 and CD11c expression.
Flow Cytometry
Spleen single cell suspensions were prepared by gently separating single cells between the frosted areas of two microscopy slides and red blood cells were lysed using 1x ACK buffer (0.15 M NH4Cl, 0.01 M KHCO3, 0.1 mM EDTA, pH 7.3). Single cells were incubated with unlabeled anti-CD16/32 antibodies in 1x phosphate buffered saline (PBS) for 20 min. after which fluorescently labeled antibodies specific for CD3, CD4, CD8, CD11b, CD11c, CD19, CD21/35, CD23, CD25, CD38, CD44, CD62L, CD69, B220, GL-7, Gr1, IgD, IgM, SiglecH (all from eBiosciences Inc., San Diego, CA), CD138 (BD Pharmingen, San Jose, CA), were added in combinations and samples were incubated for an additional 30 minutes. Flow cytometry was performed on a LSR Fortessa (BD Biosciences) and all analyses were performed using FlowJo version 9.7.5 or later (Tree Star Inc., San Carlos, CA).
In Vitro B Cell Differentiation
B220+CD19highCD138- B cells from S100a9+/+ NZBWF1 male mice were plated (2 x 105/well) with 5x104 flow sorted Gr1highCD11b+ cells from either S100a9+/+ or S100a9-/- NZBWF1 male mice in the presence or absence of recombinant IFNαA (500 units/ml, PBL InterferonSource, Piscataway, NJ) and recombinant CD40L (10µg/ml, EBiosciences) as previously described (14). After 72 hrs, 25,000 cells per well were transferred to ELISPOT plates (EMDMillipore, Burlington, MA) pre-coated with anti-Ig antibody (Southern Biotech, Birmingham, AL) for additional 20 hrs of culture. Plates were developed using HRP-conjugated anti-IgG and anti-IgM specific secondary antibodies (KPL, Milford, MA and Southern Biotech, respectively) and AEC Substrate set (BD Biosciences). The numbers of IgM and IgG-secreting cells were detected on an ELISPOT reader (CTL Immunospot, New York, NY) and number of spots per 105 B cells were calculated.
Immunohistochemistry
For histological analyses, kidneys were harvested at the end of the study (8 months of age) or earlier if mice presented with severe proteinuria. One half kidney was fixed in 10% formalin for at least 24 hours, transferred to 80% ethanol and embedded in paraffin. Five µm sections were cut and stained with hematoxylin/eosin (Newcomer Supply, Middleton, WI) or Masson’s Trichrome (Thermo Fisher Scientific, Waltham, MA) according to the manufacturer’s guidelines. Whole sections were scored in a blinded fashion by a renal pathologist (JN) at the Cleveland Clinic. Renal score: kidneys were evaluated on a scale of 0-5 for each of the following characteristics: mesangial hypercellularity, endocapillary hypercellularity, extracapillary proliferation (crescents), immune deposits, tubular atrophy, tubular casts, tubular dilation and interstitial fibrosis and inflammation for a maximum score of 40. 0 = absent, 1 = 1-5%, 2 = 6-10%, 3 = 11-20%, 4 = 21-50%, 5 > 50%. Glomerular area was calculated by measuring the area of 5-15 individual glomeruli per section for each mouse using the Keyence BZ-X analysis software (Keyence, Osaka, Osaka, Japan) and the Keyence BZ-X700 All-in-one microscope (Keyence, Osaka, Osaka, Japan). For identification of Mac2, S100a8 and S100a9, sections were blocked in 2% fetal bovine serum (FBS) in Hank’s Balanced Salt Solution (HBSS) for 30 min. Primary rat anti-mouse MAC2 (clone 125401)(BioLegend), rat anti-mouse S100a8 (clone MAB3059) or mouse S100A9 clone AF2065) (both from (R & D Systems, Minneapolis, MN) antibodies were added at 1:1500, 1:250 and 1:250 dilution, respectively, and sections were left to incubate for 1 hour in a humidified chamber at room temperature. After washing with 1x PBS, slides were applied rat-on-mouse horse radish peroxidase (HRP) polymer (#RT517) (Biocare medical, Pacheco, CA) for 20 minutes at room temperature, washed with 1x PBS, after which 3,3’-diaminobenzidine (DAB) substrate was added for 1 minute. Slides were counterstained with hematoxylin 7211 (Thermo Fisher Scientific) for 1 minute, treated with Clear Rite™ and mounted for microscopy. All light microscopy images were taken on an Eclipse 55i Nikon microscope equipped with a 12 megapixel DS-Ri1 Digital Camera.
Immunofluorescence Staining
Half kidneys and 2mm cross sections of spleens were isolated and immediately frozen in OCT™. Five µm sections were cut and sections were stained for the presence of IgG and complement C’3, or B220 and GL7, respectively. Briefly, sections were fixed with cold acetone and blocked with unlabeled anti-mouse CD16/CD32 (1:200, EBiosciences) in 10% non-immune goat serum (Invitrogen). Texas-red conjugated anti-mouse IgG (1:500, Southern Biotech), FITC-conjugated anti-mouse C’3 (1:500, ICL), FITC-conjugated anti-B220 antibodies (EBiosciences), biotinylated anti-GL7 antibodies (eBioscience), and Alexa Fluor 568–conjugated streptavidin (Invitrogen) were added as indicated for each staining combination and sections were incubated overnight at room temperature. The next day, sections were washed and mounted with 70% glycerol. Imaging was done on a Keyence BZ-X700 All-in-one microscope (Keyence, Osaka, Japan) and images were quantified using the Keyence BZ-X analysis software (Keyence, Osaka, Japan). Colocalization of IgG/C’3 in renal samples is displayed by the color yellow.
Hep2-Assay
The Hep2 assay was performed according to the manufacturer’s guidelines (Bio-Rad, Hercules, CA). Positive and negative controls provided with the manufacturers’ kit were added (1:64 dilution) to two wells on each slide. Serum from S100a9+/+, S100a9+/- and S100a9-/- mice were diluted 1:10 in 1x phosphate buffered saline and added to remaining wells. Images were obtained on a Leica Leitz fluorescence microscope and processed using Image-Pro Plus software.
ELISA
Serum anti-chromatin and anti-histone IgG levels were determined as previously described using serum diluted 1:300 in serum diluent (43). Briefly, microtiter plates (Immulon 2HD, Thermo Fisher Scientific) were coated with purified chromatin or total histones overnight at 4°C, blocked in 5% gelatin/PBS for >2 hours, and incubated with diluted serum samples for 2 hours. Secondary horseradish peroxidase–conjugated anti-mouse IgG antibodies (Invitrogen, Calsbad, CA) were added for 1.5 hours, and the plates were developed using 10 mg/ml 2,2’-Azinobis [3-ethylbenzothiazoline-6-sulfonic acid (ABTS) in McIlwain’s buffer (0.09M Na2HPO4, 0.06M citric acid, pH 4.6) or a 3,3′,5,5′-Tetramethylbenzidine (TMB) Substrate Kit (Thermo Fisher Scientific). Anti-dsDNA IgG levels were determined using serum diluted 1:100 according to the manufacturer’s guidelines (Alpha Diagnostics, Santa Monica, CA). Total serum IgG and total serum IgM levels were determined on serum diluted 1:100,000 as previously described (43). Serum S100a8/a9 heterodimer was detected using the mouse S100a8/S100a9 ELISA kit on serum diluted 2-4 times (MyBiosource, CA, USA). IgG1 specific to 4-Hydroxy-3-nitrophenylacetyl-chicken gamma globulin (NP-CGG) was determined as previously described using serum diluted 1:10,000 (42). Serum levels of soluble RAGE (sRAGE) and HMGB1 were performed on 1:10 diluted serum from 4 month old mice according to the manufacturer’s instructions (MyBiosource). Serum BAFF levels were determined on 1:3-1:12 diluted serum samples according to the manufacturer’s instructions (R&D Systems). All samples were run in duplicate. All ELISA plates were read on a Victor™ plate reader (Perkin Elmer) at 405nm or 450 nm. When possible, concentrations were calculated based on standard curves provided with the individual kits. All analyses were done using GraphPad Prism v. 5.02 (San Diego, CA).
Real Time RT-PCR
For detection of S100a4, S100a8 and S100a9 transcript levels Gr1highCD11b+ cells were flow sorted from male and female 4 week old NZBWF1 mice. RNA was isolated using the Micro RNeasy kit (Qiagen, Germantown, MD) and cDNA was prepared using qScript cDNA Supermix (QuantaBio, Beverly, MA). Transcript levels were determined in a SyBr-green based assay (PerfecCTa ® SYBR ® Green Fastmix ® ROX (QuantaBio)) using the following primers: S100a4 Forward 5’-ctactgaccagggagctg c-3’, S100a4 Reverse 5’-tgttgctgtccaagttgctca-3’, S100a8 Forward 5’-gagtgtcctcagtttgtgcagaa-3’, S100a8 Reverse 5’-tgagatgccacacccactttt-3’, S100a9 Forward 5’-gaagcacagttggcaacctt-3’, S100a9 Reverse 5’-caggtcctccatgatgtcat-3’, β2M Forward 5’-tcagtcgcggtcgcttc-3’, β2M Reverse 5’-caagcaccagaaagactagggtc-3’. For detection of Irf7, Isg15, Ifi202, and Ifna total splenocytes from 8 months old S100a9+/+ and S100a9-/- NZBWF1 male and female mice were isolated, RNA extracted and cDNA prepared as described above. βeta 2-microglubulin (β2M) was used as the internal control gene and the following gene specific primers were used: Irf7 Forward 5’-gcgtaccctggaagcatt tc-3’, Irf7 Reverse 5’-gcacagcggaagttggtct-3’, Isg15 Forward 5’-ggtgtccgtgactaactccat-3’, Isg15 Reverse 5’- tggaaagggtaagaccgt cct-3’, Ifi202 Forward 5’-caagcctctcctggacctaa-3’, Ifi202 Reverse 5’-ctaggatgccactgctgttg-3’, Pan-Ifna Forward 5’-cttccacaggatcactgtgta cct-3’, Pan-Ifna Reverse 5’-ttctgctctgaccacctccc-3’.
Statistical Analyses
GraphPad Prism was used for all statistical analyses. One-way ANOVA tests were done for all group comparisons. One-way ANOVA results are provided in figure legends. The Kruskal-Wallis test followed by Dunn’s multiple comparisons test was done for one-time point proteinuria scores (non-Gaussian distribution). Longitudinal proteinuria data were evaluated using a log rank test. Student’s t-test with Welch’s correction was done for pair-wise comparisons. p < 0.05 was considered statistically significant for all analyses.
Results
S100a9 mRNA Expression Is Elevated in Gr1highCD11b+ Cells From Male NZBWF1 Mice
Gr1highCD11b+ MDSC-like cells exert immunosuppressive functions and are increased in male NZBWF1 mice as compared to female littermates (14). We have previously shown that female Gr1highCD11b+ cells from 4 week old NZBWF1 mice utilized ROS/NOS production to inhibit cytokine-driven B cell differentiation, but failed to identify the mechanism used by male NZBWF1-derived cells (14). S100a8/a9 heterodimers are produced predominantly by neutrophils (30, 31) and have previously been associated with immunosuppression (25). We isolated Gr1highCD11b+ and Gr1lowCD11b+ cells from 9 week old male and female NZBWF1 lupus-prone mice and determined levels of S100a4, S100a8, and S100a9 mRNA. S100a9 mRNA levels were significantly elevated in male Gr1highCD11b+ cells as compared to female-derived cells (p < 0.05) (Figure 1A). A similar trend was seen in Gr1lowCD11b+ cells (p = 0.08). Levels of S100a8 mRNA were also elevated in male Gr1highCD11b+ cells as compared to female-derived cells (p = 0.08) (Figure 1B), while levels of S100a4 mRNA were unchanged between males and females (Figure 1C).
Figure 1 S100a9 expression is elevated within Gr1highCD11b+ cells of male lupus-prone mice and regulates B cell differentiation in vitro and antibody production in vivo. Gr1highCD11b+ and Gr1lowCD11b+ cells were isolated from 9 week old male and female NZBWF1 lupus-prone mice by flow cytometry and levels of S100a9 (A), S100a8 (B), and S100a4 (C) mRNA were determined by RT-PCR analysis (n = 4). Graphs show mean +/- SEM. (D) Supernatants from overnight cultured Gr1highCD11b+ cells from 4 wk old NZBWF1 male mice suppress IFNα/CD40L-driven B cell differentiation. Shown is the average (± SEM) of 4 independent experiments. (E) S100A8/A9 heterodimer levels in serum, spleen and bone marrow from 9 wk old NZBWF1 mice. Each symbol represents one mouse. Male vs. female: p = 0.07; two-way ANOVA. (F) Gr1highCD11b+ cells from 9 wk old S100a9-deficienct male NZBWF1 mice fail to suppress cytokine-driven B cell differentiation in vitro. IgM/IgG-secreting cells were enumerated by ELISPOT and are presented as number of cells per 105 plated B cells. Data shown represent the mean ± SEM of 4 independent assays. *P < 0.05; **P < 0.01, Student’s unpaired t test.
To investigate whether the previously reported mechanism of suppression by male Gr1highCD11b+ cells depended on a secreted factor, we cultured flow sorted NZBWF1 male CD19+CD138- B cells with supernatants from overnight cultured flow-sorted male Gr1highCD11b+ cells in the presence of differentiation-inducing cytokines (CD40L and IFNα). B cell differentiation was effectively blocked in a dose-dependent manner (Figure 1D), suggesting that the inhibitory factor was indeed secreted by Gr1highCD11b+ cells. As S100a8/S100a9 heterodimers can be secreted by neutrophils, we evaluated levels of S100a8/S100a9 heterodimers in serum, spleen, and bone marrow samples. S100a8/S100a9 heterodimers were slightly increased in male NZBWF1 mice at all sites (p = 0.07) (Figure 1E), prompting us to further investigate a possible role for S100a9 in Gr1highCD11b+ MDSC-like cell-mediated immunosuppression during lupus-like disease development via the generation of S100a9-deficient NZBWF1 (see Supplemental Figure 1). Gr1highCD11b+ cells were sort purified from S100a9+/+ and S100a9-/- male NZBWF1 mice and tested for their ability to suppress B cell differentiation in vitro as previously described (14). Supporting an immunosuppressive role, S100a9-sufficient, but not S100a9-deficient, male Gr1highCD11b+ cells inhibited B cell differentiation in vitro (Figure 1F).
Elevated Antibody Responses in Male S100a9-/- NZBWF1 Mice
To determine the in vivo effect of S100a9-deficiency, S100a9+/+, S100a9+/- and S100a9-/- NZBWF1 mice were immunized with NP-CGG in CFA, and serum evaluated for the levels of NP-specific IgG1 levels. Supporting an immunosuppressive role in vivo, Male S100a9-/- and S100a9+/- mice displayed significantly increased anti-NP IgG1 levels as compared to male S100a9+/+ mice (p < 0.01; Figure 2A). In contrast, immunization of S100a9+/+, S100a9+/- and S100a9-/- female NZBWF1 mice resulted in similar levels of NP-specific antibodies between the strains (Figure 2B), suggesting a sex-specific effect of S100a9 in NZBWF1 lupus-prone mice. Four weeks post immunization, animals were euthanized and levels of splenic germinal center (GC) B cells and T follicular helper (Tfh) cells were analyzed. Further supporting a male-specific effect, percentages of GC B cells were significantly elevated in male, but not female, S100a9-/- NZBWF1 mice (Figure 2C). In contrast, the frequencies of Tfh cells were unchanged between all the groups, likely reflecting the advanced time point of evaluation (data not shown).
Figure 2 Enhanced response to T-dependent Ag immunization in male S100a9-/- NZBWF1 mice. Male (A) S100a9+/+ (black squares), S100a9+/- (grey diamonds), S100a9-/- (symbol X), and female (B) S100a9+/+ (black triangles, up), S100a9+/- (grey triangles, down), S100a9-/- (open circles) NZBWF1 mice (n = 3-4) were immunized on day 0 with NP27-CGG in CFA. Levels of NP-specific IgG1 antibodies were determined on days -1, 7, 14, 21 and 28. *p < 0.05, two-way ANOVA. (C) GC B cells were identified by flow cytometry four weeks post immunization. *p < 0.05, Student’s unpaired t test with Welch’s correction. One-way ANOVA: p < 0.05.
Male S100a9-/- NZBWF1 Mice Develop Splenomegaly And Present With Elevated Anti-Nuclear Autoantibodies
Splenomegaly is a well-established indicator of mouse lupus-like disease. Male S100a9-/- and S100a9+/- NZBWF1 mice displayed significantly larger spleens than male S100a9+/+ mice (p < 0.05, Figure 3A). In contrast, no significant difference was observed between S100A9-sufficient and S100A9-deficient female mice. A similar trend was observed for total splenocyte count, although this measure only reached statistical significance when comparing S100a9-/- and S100a9+/+ male mice (p < 0.05, Figure 3B). As expected, S100a9-sufficient female mice exhibited significantly larger spleens than S100a9-sufficient male mice (p < 0.05-0.01, Figures 3A, B) and S100a9-/- females (p < 0.05), suggesting that S100a9 might play a different role in disease progression in older female NZBWF1 mice.
Figure 3 S100a9-/- male BWF1 mice develop splenomegaly and elevated anti-nuclear autoantibodies. Spleens were harvested from eight month old BWF1 lupus-prone mice. Spleen weight (A) and splenocyte count (B) were measured (n = 5-11). Each symbol represents one individual mouse, horizontal lines represent mean values. C-G) Serum was obtained from five months old BWF1 lupus-prone mice and serum autoantibody levels were determined by ELISA: serum IgM (C), serum IgG (D), anti-histone IgG (E), anti-chromatin IgG (F), anti-dsDNA (G). Each symbol represents one individual mouse (n = 5-11). (H) Hep2 Assay depicting nuclear and perinuclear antibody staining. Representative pictures shown from 5-11 mice/group. *p < 0.05; **p < 0.01; ***p < 0.001, Student’s unpaired t test with Welch’s correction. One-way ANOVA: p = 0.21 (A), p < 0.05 (B), p < 0.01 (C), p < 0.01 (D), p < 0.01 (E), p < 0.0001 (F), p = 0.87 (G).
We hypothesized that S100a9-deficient male, but not female, NZBWF1 mice would present with elevated levels of serum anti-nuclear autoantibodies; a hallmark of lupus-like disease in this animal model. Cohorts of male and female S100a9+/+, S100a9+/-, and S100a9-/- NZBWF1 mice were followed for up to 8 months of age or until severe renal disease was present (n = 5-11). Serum was obtained monthly and tested for total IgG, total IgM, and anti-nuclear autoantibodies by ELISA. At 4-5 months of age, serum IgM levels were significantly increased in all female mice as compared to male mice (p < 0.05), but there was no major difference between mice of similar sex (Figure 3C). In contrast, serum total IgG was significantly elevated in male S100a9-/- NZBWF1 mice as compared to S100a9+/- and S100a9+/+ male littermates (p < 0.05, Figure 3D), while there was no difference between female mice. We also observed a modest elevation in serum anti-histones IgG and serum anti-dsDNA IgG levels and a significant increase in anti-chromatin IgG levels in male S100a9-/- mice (p < 0.05), with concentrations approaching those measured in the female cohorts (Figures 3E–G). Interestingly, there was a trend towards less serum anti-chromatin IgG in female S100a9-/- mice as compared with female S100a9+/+ mice (Figure 3F, p = 0.07). Further analysis of the nature of serum autoantibodies from S100a9+/+, S100a9+/- and S100a9-/- male and female NZBWF1 mice supported stronger and more wide-spread recognition of both cytoplasmic and nuclear antigens by serum antibodies from male S100a9-/- mice as compared with antibodies from male S100a9+/+ mice, but no difference between female S100a9+/+, S100a9+/- and S100a9-/- mice (Figure 3H).
S100a9-/- Male NZBWF1 Mice Develop Spontaneous B Cell Hyper-Activation
Lupus-like disease in NZBWF1 mice is characterized by hyperactive B cells, abnormal germinal center (GC) formation, and an accumulation of memory B cells and plasma cells (44, 45). Despite the difference in spleen size, there was no difference in the percentages of total B220+ cells and follicular mature B cells (CD23highCD21lowIgMlow) in the spleens of male and female S100a9-deficient and -sufficient NZBWF1 mice (Figures 4A, B). In alignment with previous findings (46), female S100a9+/+ NZBWF1 mice displayed reduced percentages of MZ B cells (B220+CD21highCD23low) as compared with S100a9+/+ male mice, however a similar pattern was not observed in male S100a9-/- mice (Figure 4C). As expected, the percentage of GC B cells was increased in female S100a9+/+ mice as compared to male S100a9+/+ mice (p = 0.057), but there was no significant difference in the percentages of GC B cells between S100a9-sufficient and -deficient male mice (Figure 4D), despite a significant increase in GC area (p < 0.05) (Figures 4G, H). Female S100a9-/- mice showed a decreased frequency of GC B cells as compared with female S100a9+/+ mice (p < 0.05)(Figure 4D), a pattern that was mimicked by smaller size GCs in female S100a9-/- mice (Figures 4E, F). Finally, frequencies of memory B cells and plasma cells (PCs) were elevated in male S100a9-/- mice as compared with male S100a9+/+ mice (p = 0.08 and p = 0.05, respectively) and approached levels observed in S100a9+/+ female mice (Figures 4G, H). Taken together, several B cell subsets, plasma cells and GCs from male S100a9-/- NZBWF1 mice acquired an intra-splenic balance approaching that observed in female S100a9+/+ NZBWF1 mice, supporting a role for S100a9 in the regulation of B cell differentiation in male NZBWF1 mice.
Figure 4 S100a9-deficiency results in increased spontaneous germinal reactions in male BWF1 mice. Spleens were harvested from eight month old BWF1 lupus-prone mice and percentages of B cell subsets were determined by flow cytometry: (A) total B220+ cells, (B) Follicular mature (B220+CD23highCD21lowIgMlow) B cells, (C) Marginal Zone B cells (B220+CD21highCD23low IgMhigh), (D) Germinal center B cells (B220+GL7+CD38lowIgMlow), (E) memory B cells (B220+ CD38highGL7-IgM-) and (F) Plasma cells (B220-/lowCD138+IgM-IgD-). G) Frozen spleen sections were stained for B cells (B220-FITC) and germinal centers (GL-7-TexasRed). Representative pictures are shown. H) Germinal center areas from (G) were determined using the Keyence BZ-X analysis software. (A–F, H) Each symbol represents one individual mouse (n = 5-11). *p < 0.05; **p < 0.01, Student’s unpaired t test with Welch’s correction. One-way ANOVA: p = 0.13 (A), p = 0.14 (B), p = 0.33 (C), p < 0.01 (D), p = 0.36 (E), p < 0.01 (F), p < 0.01 (H).
S100A9-Deficiency Drives Spontaneous T Cell Activation and Differentiation in Male NZBWF1 Mice
It has previously been described that B cell hyper-activation and autoantibody production in female NZBWF1 mice depends on T cell help (47). We determined the percentages of T helper cell subsets in S100a9-sufficient and -deficient, male and female mice. There was no statistically significant difference in the percentages of total CD4+ T cells among the cohorts of male and female mice (Figure 5A). The expression of CD25 and CD69, was significantly elevated on female S100a9+/+ CD4+ cells as compared to male S100a9+/+ CD4+ cells (p<0.01), while male S100a9-/- NZBWF1 mice displayed elevated percentages of CD25+CD4+ cells, but not CD69+CD4+ cells, as compared with male S100a9+/+ NZBWF1 mice (p < 0.05) (Figures 5B, C). Interestingly, we again observed a partial normalization in the levels of CD69+CD4+ T cells in S100a9-/- female mice as these approached levels in S100a9+/+ males, but no significant change in the frequency of CD25+CD4+ T cells between female mice (Figures 5B, C).
Figure 5 S100a9-deficiency promotes the activation and differentiation of CD4+ T naïve cells to effector memory cells in S100a9-/- male BWF1 mice. Spleens were harvested from eight month old BWF1 lupus-prone mice and percentages of T lymphocytes and their subsets were determined by flow cytometry: total CD4+ cells (A), CD25+ CD4+ cells (B), CD69+ CD4+ cells (C), naïve (CD44lowCD62Lhigh) CD4+ cells (D), effector memory (CD44highCD62Llow) CD4+ cells (E), and the ratio between naïve and effector memory cells (F). Each symbol represents one individual mouse. *p < 0.05; **p < 0.01, Student’s unpaired t test with Welch’s correction. One-way ANOVA: p < 0.05 (A), p < 0.05 (B), p < 0.0001 (C), < 0.0001 (D), p < 0.05 (E), p < 0.01 (F).
Finally, male S100a9+/+ NZBWF1 mice expressed high levels of naïve CD4+ T cells (CD62LhighCD44low) and low levels of effector/memory CD4+ T cells (CD62LlowCD44high), while the balance was shifted in both male S100a9-/- and female S100a9+/+ NZBWF1 mice (p = 0.05 and p < 0.01, respectively) (Figures 5D–F). Thus, male S100a9-/- NZBWF1 mice presented with a hyperactive T cell phenotype similar to that of female S100a9+/+ NZBWF1 mice, further supporting an immunosuppressive role for S100a9 in male NZBWF1 mice.
S100a9-/- Male NZBWF1 Mice Develop Enhanced Glomerulonephritis and Renal Damage
Lupus-like renal disease in female NZBWF1 mice is characterized by the development of glomerulonephritis, the deposition of IgG-immune complexes (IgG-IC) in the glomeruli, complement fixation, and renal damage (44, 48). Kidneys were harvested from all mice at the end of the study and evaluated by immunostainings. Male S100a9+/+ mice appeared relatively healthy at eight months of age with small glomeruli and minimal inflammation (Figure 6A, panels a, g). In contrast, S100a9+/- and S100a9-/- male mice displayed a progressive worsening as characterized by increased mesangial proliferation (black open-head arrows) and interstitial fibrosis (Figure 6A, panels b, c), increased collagen deposits (closed-head blue arrows), tubular atrophy, and visible proteinaceous material (black asterisk) (Figure 6A, panels h, i). Irrespective of their genotype, all females displayed significantly more damaged kidneys with significantly increased mesangial cell proliferation, interstitial fibrosis, and tubular casts than their male counterparts (Figure 6A, panels d-f and j-m).
Figure 6 S100a9-deficiency promotes proteinuria and renal damage in male BWF1 mice. Kidneys were harvested from eight month old BWF1 lupus-prone mice (n = 5-11) and stained. (A) H&E (a-f) and Masson’s Trichrome (g-m) were performed to assess glomerulonephritis, mesangial proliferation (black open-head arrows), interstitial inflammation, tubular atrophy, collagen deposits (blue closed-head arrows), and cast formation (red asterisk). All images were taken at similar settings and magnifications. (B, C) Deposition of IgG-immune complexes (red) and fixation of complement factor 3 (green) were determined by immunofluorescent staining of kidneys (B, n-s). Colocalization of IgG and C’3 can be seen as yellow. IgG deposition was quantified (C). (D) Proteinuria readings at the time of harvest. (E) Renal histology scores. Each symbol represents one individual mouse. *p < 0.05; **p < 0.01, Student’s unpaired t test with Welch’s correction (C, D), Student’s paired t- test (E). One-way ANOVA: p < 0.01 (C), p < 0.01 (D) and p < 0.01 (E).
Kidneys were further evaluated for the presence of IgG-IC deposition and complement factor C’3 fixation by immunostaining. S100a9-/- males displayed increased levels of IgG-IC deposition (red color), increased C’3 fixation to the glomeruli (green color), and expanded areas of co-localization with C’3 (yellow color) as compared with S100a9+/+ males (Figure 6B, panels n-p). In contrast, all female mice showed comprehensive overlapping IgG-IC deposition and complement C’3 fixation (Figure 6B, panels q-s, yellow color). Quantification of IgG-IC deposition showed statistically more depositions in S100a9-/- males and S100a9+/+ females as compared with S100a9+/+ males (p < 0.05 and p < 0.01, respectively) (Figure 6C). As renal damage progresses to renal failure in female NZBWF1 mice, proteinuria levels were measured monthly in all animals starting at 3-4 months of age to monitor the onset of renal failure. At 8 months of age, levels of proteinuria were significantly increased in male S100a9-/- and S100a9+/- NZBWF1 mice as compared to S100a9+/+ male littermates (p < 0.05), while there was no statistical difference between the female mice (Figure 6D). The difference among male mice was not apparent at earlier time points (Supplemental Figure 2A). There was also no difference in the average time for disease to develop among female mice: S100a9+/+: t = 6.8 months (n = 5), S100a9+/-: t = 6.8 months (n = 5), and S100a9-/-: t = 6.875 months (n = 8) months, indicating that S100a9-deficiency did not accelerate disease development in the females (Supplemental Figure 2A). Finally, renal damage was supported by an elevated renal score in female S100a9+/+ mice (p < 0.01), and a trend towards a higher score in male S100a9-/- mice (p = 0.08). There was no statistically significant difference among the groups of female mice (Figure 6E). Given the presence of areas of fibrosis, we analyzed kidney sections for levels of pro-inflammatory Mac2+ macrophages. Numbers of Mac2+ macrophages within glomeruli increased slightly in both S100a9+/- and S100a9-/- male mice as compared with S100a9+/+ male mice, albeit significant differences were not observed (Supplemental Figures 2B, C). It should be noted that a significant increase in tubular Mac2 expression was similarly observed in S100a9-deficient mice (males and females), although the significance of this finding remains unknown.
S100a9-Deficiency Does Not Alter Levels of HMGB1 and RAGE
While S100a9-deficiency drives disease development in male NZBWF1 mice, several of the humoral readouts tested suggested an immunostimulatory function of S100a9 in female animals. S100a9 has been shown to bind to multiple inflammatory molecules including high-mobility group box 1 (HMGB1) and receptor for advanced glycation end products (RAGE). We therefore tested levels of serum HMGB1 and soluble RAGE in male and female S100a9+/+, S100a9+/- and S100a9-/- NZBWF1 mice at 4-5 months of age. Neither molecule was affected by expression of S100a9 and neither differed between male and female mice (Supplemental Figure 3), suggesting that these pathways were not accountable for the different outcome of S100a9-deficiency in males and females.
S100a9-Deficiency Drives Neutrophil Accumulation in Spleens of Male NZBWF1 Mice, but No Change in BAFF Levels
Further understanding the mechanism behind lupus-like disease development in S100a9-deficient male NZBWF1 mice is imperative for the identification of de novo molecular or cellular therapeutic targets. Previously, a population of B cell activating factor (BAFF)-producing B helper neutrophils (NBH) was described located in the marginal zone in secondary lymphoid organs in response to inflammation, either induced or chronic as in SLE (49). We therefore tested levels of mature neutrophils (Gr1highCD11b+SSChigh) in spleens of S100a9-deficient and –sufficient, male and female NZBWF1 mice along with systemic levels of BAFF (Figures 7A, B). While levels of Gr1highCD11b+SSChigh cells were ~ 3 fold higher in S100a9-/- male NZBWF1 mice than in S100a9+/+ males and all female mice, we found no differences in levels of serum BAFF, suggesting that another mechanism may be responsible for driving lupus-like disease in male S100a9-deficient mice.
Figure 7 Male S100a9-deficient BWF1 mice express elevated levels of mature neutrophils, but unaltered levels of BAFF. (A) Levels of Gr1highCC11b+SSChigh mature neutrophils was measured in spleens of male and female S100a9+/+, S100a9+/- and S100a9-/- mice by flow cytometry. (B) Serum BAFF levels were measured by ELISA on serum from 4 month old mice. Each dot represent one mouse. Shown is Mean ± SEM. *p < 0.05; **p < 0.01. One-way ANOVA: p < 0.05 (A), p = 0.14 (B).
S100a9-Deficiency Results in Elevated ISG Signature in Male NZBWF1 Mice
It has been suggested that S100a9, granulocytic MDSCs (Gr1highCD11b+) and the interferon-related factor 7 (IRF7) exist in an auto-regulatory loop (50). Since IRF7 is essential for type I interferon production and type I interferons are known to drive lupus-like disease in NZBWF1 mice (51, 52), we tested transcript levels of Irf7 and interferon-stimulated genes Ifi202 and Isg15 in a small subset of S100a9-sufficient and -deficient male and female NZBWF1 mice. Interestingly, all three genes were significantly upregulated in S100a9-/- male NZBWF1 mice as compared to S100a9+/+ male NZBWF1 mice (p < 0.05-0.01) (Figures 8A–C). In contrast, neither gene was upregulated in S100a9-/- female NZBWF1 mice as compared with S100a9+/+ females. Both Irf7 and Isg15 displayed slightly increased levels in WT females as compared to WT males. In lupus, type I interferons have been shown to be produced by either SiglecH+ plasmacytoid dendritic cells (pDCs) or low density granulocytes (LDGs) (53, 54). To evaluate if levels of these cell subsets were affected by S100a9-deficiency and could account for the suspected increased levels of type I interferons, we analyzed splenic frequencies of both cell subsets. We found no differences in the levels of SigH+ pDCs between S100a9+/+ and S100a9-/- mice of either sex, despite elevated levels of total pDCs in female S100a9+/+ mice as compared with male S100a9+/+ mice (Figures 8D, E). In contrast, levels of total Gr1+CD11b+SSClow LDGs were significantly elevated in male S100a9-deficient NZBWF1 mice (p < 0.05), but not in female S100a9-deficient mice (p = 0.12)(Figure 8F). Further analysis of LDGs in NZBWF1 mice identified elevated Ifna transcripts in both Ly6Chigh (Gr1low) and Ly6Clow (Gr1high) LDG subsets, although differences only reached statistical significance in the Ly6ChighCD11b+SSClow population (Supplemental Figure 4).
Figure 8 Male S100a9-deficient BWF1 mice express elevated expression of interferon-stimulated genes and increased levels of low-density granulocytes. (A–C) Real-time RT-PCR was performed on total splenocytes obtained from a subset of S100a9+/+ and S100a9-/- male and female BWF1 mice at the time of harvest (≥ 7 months of age) and levels of Ifi202 (A), Irf7 (B), and Isg15 (C) were determined. Levels of total pDCs (D), SigH+ pDCs (E), and Low-density granulocytes (LDGs) (F), were determined by flow cytometry. Each dot represent one mouse. Shown is Mean ± SEM. *p < 0.05; **p < 0.01, ns, not statistically significant. One-way ANOVA: p < 0.05 (A), p = 0.08 (B), p = 0.07 (C), p = 0.09 (D), p = 0.67 (E), p = 0.0001 (F).
Discussion
Lupus, along with many other autoimmune disorders, is more prevalent in females than in males. While this discrepancy has been extensively investigated from the point of view that the immune system of females is hyperactive and therefore more likely to drive autoimmunity, only a few studies have evaluated why otherwise genetically predisposed males fail to develop disease. We previously reported that lupus-prone male NZBWF1 mice exhibited elevated levels of immunosuppressive neutrophil-like cells (Gr1highCD11b+ granulocytic MDSCs) and that depletion of these cells stimulated ANA production and IgG-IC deposition within the kidney glomeruli of lupus-prone male NZBWF1 mice (14). While we showed that the corresponding female MDSCs were immunosuppressive via the production of ROS/NOS only before puberty (at 4 weeks of age), male MDSCs remained immunosuppressive until at least 16 weeks of age (14). Here we show that MDSCs from male NZBWF1 mice utilize in part S100a9 as a mechanism of suppression, as cells from S100a9-/- male NZBWF1 mice were unable to inhibit cytokine-driven B cell differentiation in vitro, control antibody responses to T-dependent NP-CGG antigen, and most importantly, control lupus-like disease development. As such, B and T cell hyperactivation, ANA production, IgG-IC deposition and complement C’3 fixation in kidney glomeruli, as well as renal damage, were all significantly enhanced in S100a9-deficient NZBWF1 male mice, but largely unaffected female S100a9-deficient NZBWF1 mice.
We previously suggested that the inflammatory milieu in lupus-prone female NZBWF1 mice may drive the maturation of MDSCs into pro-inflammatory neutrophils and/or macrophages, as previously shown (55). Subsequent studies by others showed that MDSC-like cells from female NZBWF1 mice were indeed effectively eliminated in a process involving IFNα, IFNγ, IL-6 and ROS production (15, 49, 56). Oppositely, depletion of neutrophils (including MDSCs) resulted in increased IFNγ production by NK cells (57). In lupus, both IFNα and IFNγ have been suggested as drivers of disease. For example, blocking or elimination of the IFNγ receptor was previously shown to significantly reduce disease development in female NZBWF1 mice (58), showing that IFNγ contributes to disease. Interestingly, NK cell driven IFNγ production is induced by pDC-derived IFNα, and inhibition of IFNα or the IFNα/β-receptor has similarly been shown to ameliorate lupus-like disease (52, 59). Finally, MDSCs were shown to be upregulated in the absence of Irf7 in a tumor model (50). Irf7 is induced downstream of both Toll-like receptor 7 (TLR7), TLR9, and type I interferon receptor ligation, all key molecules in mouse lupus pathogenesis. Thus, it is likely that in the presence of high levels of IFNα, as seen in some lupus patients and many mouse models of lupus, Irf7 is upregulated and MDSCs downregulated. Interestingly, Irf7 binding sites have been identified in the S100a9 promoter (50), further supporting an autoregulatory loop between S100a9, MDSCs and Irf7. In our study, we found that Irf7 mRNA was upregulated in male S100a9-/- NZBWF1 mice as compared to their S100a9-sufficient littermates. We suggest that Irf7 mRNA is upregulated due to elevated levels of IFNα, which is supported by the concomitant upregulation of two other ISGs: Isg15 and Ifi202, and the accumulation of Ifna-expressing LDGs in S100a9-/- NZBWF1 males. Thus, in this system, a lack of S100a9 interferes with the function of MDSCs, which in turn may lead to an accumulation of IFNα-producing LDGs, driving T and B cell activation, autoantibody production and end-organ inflammation and damage.
Surprisingly, female S100a9-deficient NZBWF1 mice displayed some evidence for reduced humoral disease as determined by reduced levels of activated T and B cells, reduced GC formation and diminished splenomegaly, although no differences in renal parameters were observed. In that regard, S100a9 has been found to possess both pro- and anti-inflammatory properties under different conditions. For example, some studies indicate that S100a8/a9 possess pro-inflammatory functions consistent with the complex being a TLR4 ligand (60, 61), while other evidence has suggested an immunoprotective role (33, 34, 39, 62). The two conclusions are not mutually exclusive and may depend on one of several variables such as the timing of investigation, whether inflammation is infectious or sterile, which cytokines are produced, which organ system is involved, and the sex of the animal. More studies are needed to determine if the latter is indeed involved, as the sex of the animals studied was not consistently reported in several of these studies, and no study compared the response in males and females. Both lupus patients and female NZBWF1 mice express elevated levels of several pro-inflammatory cytokines including IFNα, TNFα, IFNγ that may act to differentiate immunosuppressive MDSCs into proinflammatory MDSCs. In fact, we observed previously that Ly6ChighCD11b+ MDSCs became immunostimulatory in female NZBWF1 mice as the mice aged (14). Here we observed elevated levels of Gr1highCD11b+SSChigh neutrophilic cells in male S100a9-deficient mice, suggesting that these cells may be involved in the pathogenesis. Others have found that a subset of B helper neutrophils capable of producing BAFF was significantly increased in autoimmune mice and people with autoimmune diseases (49) however, we found no evidence for differential levels of BAFF in S100a9-sufficient and –deficient NBWF1 mice, and thus a potential pathogenic mechanism by mature neutrophils in S100a9-deficient male mice remains unknown. Finally, it should also be mentioned that the conformation of S100a9-contaning multimers can affect the function of the complex. As such, it has been shown that S100a8/S100a9 heterodimers can drive an inflammatory immune response via binding to TLR4/Mdm2 receptors, while (S100a8/a9)2 tetramers are unable to bind TLR4 and thus unable to drive inflammation (63). Whether S100a9-containing complexes are different in male and female NZBWF1 mice, and how such complex formation is regulated, remains to be fully determined.
Besides transcriptional regulation as discussed above, post-translational modifications have been implicated in the function of S100a8 and/or S100a9. Post-translational modifications, including S-nitrosylation, S-gluthathionylation and phosphorylation, have been proposed to be differentially associated with pro- and anti-inflammatory functions of S100a9 (64–66). Recent studies have identified dysregulated miR-146a-5p and miR-155-5p miRNAs as drivers of phosphorylated S100a8/a9 heterodimers and the production of proinflammatory cytokines in female rheumatoid arthritis patients (67). Both of these microRNAs have also been found to be dysregulated in SLE patients (68, 69). No study has yet evaluated if this observation correlates with the presence of phosphorylated, proinflammatory S100a9 in lupus, although phosphorylated S100a9 has been found upregulated in SLE patients (70) and to exert pro-inflammatory functions (64). Oppositely, S-nitrosylated S100a8 (and S100a9 to a lesser extend) has been proposed to inhibit leukocyte-endothelial cell interactions and leukocyte extravasation due to regulation of CD11b expression (66). Interestingly, S-nitrosylation is estrogen-dependent in endothelial cells (71), although whether estrogens have a similar effect on S100a8 and S100a9 in neutrophils remains unknown. Finally, the long-non-coding RNA Hotairm1 was found to drive the immunosuppressive function of MDSCs via its binding to S100a9 which functionally blocked secretion and led to an accumulation of the protein in the nucleus (72). Further studies are needed to identify the subcellular localization and post-translational modification status of S100a9 (and S100a8) in male and female lupus-prone NZBWF1 mice and SLE patients.
In summary, S100a9 is expressed and secreted by myeloid cells, including MDSCs, and exerts an immunosuppressive role in male, but not female, lupus-prone NZBWF1 mice. S100a9-expressing MDSCs functionally suppress B cell differentiation in vitro, B and T cell activation in vivo, and spontaneous autoantibody production. Interestingly, S100a9 may play a different role in female mice, as S100a9-deficiency results in slightly reduced disease patterns. Thus, as future studies evaluate the effect(s) of therapeutically targeting S100a9 or S100a8/S100a9 in various diseases, studying responses in both males and females will be of utmost importance to elucidate if such therapy will be equally effective in both sexes.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Ethics Statement
The animal study was reviewed and approved by Cleveland Clinic Foundation Animal Research Committee guidelines and the Institutional Animal Care and Use Committee of the Lerner Research Institute of the Cleveland Clinic Foundation.
Author Contributions
LD developed the mouse model, helped with mouse harvests, ran flow cytometry, did immunization assays, and ran the NP-ELISA. AA helped with mouse harvests, ran S100a8/a9 ELISA, and performed immunofluorescence staining on kidney tissue. AD performed Hep2G assay, dsDNA ELISA, ISH imaging, and helped write the manuscript. MP did anti-chromatin, anti-histone, IgG, and IgM ELISA, helped with mouse harvests and flow cytometry, and analyzed immunostainings of the kidneys. EK performed HMGB1, RAGE, and BAFF ELISA, ISG RT-PCR, performed and analyzed all GC immunofluorescence staining on spleen tissue, did select flow analyses, and helped revise the manuscript. ND performed trichrome staining, and S100a8 and S100a9 immunohistochemistry. AW helped with ISG RT-PCR. AK did all B cell differentiation assays. NS performed S100a4/a8/a9 RT-PCR on sorted cells. LL did flow analyses. TV generated S100a9-deficient mice and helped revise the manuscript. TJ designed the study, supervised all analyses, mentored students, and wrote the manuscript.
Funding
This study was supported in part by the Department of Defense Grant #W81XWH-11-1-0667 (TJ).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Dr. William M Baldwin 3rd for discussions and help with immunohistochemistry, and Jennifer Powers for flow assisted cell sorting.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2021.681503/full#supplementary-material
References
1. Kotzin BL, West SG. Systemic Lupus Erythrematosus. In: Rich R, Shearer WT, Kotzin BL, Schroeder JHW, editors. Clinical Immunology. Principles and Practice, 2nd edn. London, UK: Mosby International Limited (2001). pp. 60.61–60.24.
2. Gubbels Bupp MR, Jorgensen TN. Androgen-Induced Immunosuppression. Front Immunol (2018) 9:794. doi: 10.3389/fimmu.2018.00794
3. Gordon C, Wallace DJ, Shinada S, Kalunian KC, Forbess L, Braunstein GD, et al. Testosterone Patches in the Management of Patients with Mild/Moderate Systemic Lupus Erythematosus. Rheumatol (Oxford) (2008) 47:334–8. doi: 10.1093/rheumatology/kem342
4. Lahita RG, Cheng CY, Monder C, Bardin CW. Experience With 19-Nortestosterone in the Therapy of Systemic Lupus Erythematosus: Worsened Disease After Treatment With 19-Nortestosterone in Men and Lack of Improvement in Women. J Rheumatol (1992) 19:547–55.
5. van Vollenhoven RF. Dehydroepiandrosterone for the Treatment of Systemic Lupus Erythematosus. Expert Opin Pharmacother (2002) 3:23–31. doi: 10.1517/14656566.3.1.23
6. Roubinian JR, Papoian R, Talal N. Androgenic Hormones Modulate Autoantibody Responses and Improve Survival in Murine Lupus. J Clin Invest (1977) 59:1066–70. doi: 10.1172/JCI108729
7. Roubinian JR, Talal N, Greenspan JS, Goodman JR, Siiteri PK. Effect of Castration and Sex Hormone Treatment on Survival, Anti- Nucleic Acid Antibodies, and Glomerulonephritis in NZB/NZW F1 Mice. J Exp Med (1978) 147:1568–83. doi: 10.1084/jem.147.6.1568
8. Roubinian JR, Talal N, Greenspan JS, Goodman JR, Siiteri PK. Delayed Androgen Treatment Prolongs Survival in Murine Lupus. J Clin Invest (1979) 63:902–11. doi: 10.1172/JCI109390
9. Gubbels Bupp MR, Jorgensen TN, Kotzin BL. Identification of Candidate Genes That Influence Sex Hormone-Dependent Disease Phenotypes in Mouse Lupus. Genes Immun (2008) 9:47–56. doi: 10.1038/sj.gene.6364447
10. Bynote KK, Hackenberg JM, Korach KS, Lubahn DB, Lane PH, Gould KA. Estrogen Receptor-Alpha Deficiency Attenuates Autoimmune Disease in (NZB x NZW)F1 Mice. Genes Immun (2008) 9:137–52. doi: 10.1038/sj.gene.6364458
11. Dai R, McReynolds S, Leroith T, Heid B, Liang Z, Ahmed SA. Sex Differences in the Expression of Lupus-Associated miRNAs in Splenocytes From Lupus-Prone NZB/WF1 Mice. Biol Sex Differ (2013) 4:19. doi: 10.1186/2042-6410-4-19
12. Keisler LW, Vom Saal FS, Keisler DH, Walker SE. Hormonal Manipulation of the Prenatal Environment Alters Reproductive Morphology and Increases Longevity in Autoimmune NZB/W Mice. Biol Reprod (1991) 44:707–16. doi: 10.1095/biolreprod44.4.707
13. Jones JM, Jørgensen TN. Androgen-Mediated Anti-inflammatory Cellular Processes as Therapeutic Targets in Lupus. Front Immunol (2020) 11:1271:1271. doi: 10.3389/fimmu.2020.01271
14. Trigunaite A, Khan A, Der E, Song A, Varikuti S, Jorgensen TN. Gr1(high)CD11b(+) Cells Suppress B Cell Differentiation and Lupus-Like Disease in Lupus-Prone Male Mice. Arthritis Rheum (2013) 65:2392–402. doi: 10.1002/art.38048
15. Bird AK, Chang M, Barnard J, Goldman BI, Meednu N, Rangel-Moreno J, et al. Neutrophils Slow Disease Progression in Murine Lupus via Modulation of Autoreactive Germinal Centers. J Immunol (2017) 199:458–66. doi: 10.4049/jimmunol.1700354
16. Walecki M, Eisel F, Klug J, Baal N, Paradowska-Dogan A, Wahle E, et al. Androgen Receptor Modulates Foxp3 Expression in CD4+CD25+Foxp3+ Regulatory T-Cells. Mol Biol Cell (2015) 26:2845–57. doi: 10.1091/mbc.E14-08-1323
17. Frisancho-Kiss S, Coronado MJ, Frisancho JA, Lau VM, Rose NR, Klein SL, et al. Gonadectomy of Male BALB/c Mice Increases Tim-3(+) Alternatively Activated M2 Macrophages, Tim-3(+) T cells, Th2 Cells and Treg in the Heart During Acute Coxsackievirus-Induced Myocarditis. Brain Behav Immun (2009) 23:649–57. doi: 10.1016/j.bbi.2008.12.002
18. Lapierre P, Béland K, Martin C, Alvarez F Jr., Alvarez F. Forkhead Box p3+ Regulatory T Cell Underlies Male Resistance to Experimental Type 2 Autoimmune Hepatitis. Hepatology (2010) 51:1789–98. doi: 10.1002/hep.23536
19. Ma W, Zhang J, Guo L, Wang Y, Lu S, Wang Z, et al. Suppressed Androgen Receptor Expression Promotes M2 Macrophage Reprogramming Through the STAT3/SOCS3 pathway. Excli J (2019) 18:21–9.
20. Zhu W, Zhao Z, Chou F, Zuo L, Liu T, Yeh S, et al. Loss of the androgen Receptor Suppresses Intrarenal Calcium Oxalate Crystals Deposition via Altering Macrophage Recruitment/M2 Polarization with Change of the miR-185-5p/CSF-1 Signals. Cell Death Dis (2019) 10:275. doi: 10.1038/s41419-019-1358-y
21. Gabrilovich DI, Nagaraj S. Myeloid-Derived Suppressor Cells as Regulators of the Immune System. Nat Rev Immunol (2009) 9:162–74. doi: 10.1038/nri2506
22. Aochi S, Tsuji K, Sakaguchi M, Huh N, Tsuda T, Yamanishi K, et al. Markedly Elevated Serum Levels of Calcium-Binding S100A8/A9 Proteins in Psoriatic Arthritis are Due to Activated Monocytes/Macrophages. J Am Acad Dermatol (2011) 64:879–87. doi: 10.1016/j.jaad.2010.02.049
23. Boso M, Emanuele E, Minoretti P, Arra M, Politi P, Ucelli di NS, et al. Alterations of Circulating Endogenous Secretory RAGE and S100A9 Levels Indicating Dysfunction of the AGE-RAGE Axis in Autism. Neurosci Lett (2006) 410:169–73. doi: 10.1016/j.neulet.2006.08.092
24. Cesaro A, Anceriz N, Plante A, Page N, Tardif MR, Tessier PA. An inflammation Loop Orchestrated by S100A9 and Calprotectin is Critical for Development of Arthritis. PLoS. One (2012) 7:e45478. doi: 10.1371/journal.pone.0045478
25. Cheng P, Corzo CA, Luetteke N, Yu B, Nagaraj S, Bui MM, et al. Inhibition of Dendritic Cell Differentiation and Accumulation of Myeloid-Derived Suppressor Cells in Cancer is Regulated by S100A9 Protein. J Exp Med (2008) 205:2235–49. doi: 10.1084/jem.20080132
26. Jitschin R, Braun M, Buttner M, Dettmer-Wilde K, Bricks J, Berger J, et al. CLL-Cells Induce IDOhi CD14+HLA-DRlo Myeloid-Derived Suppressor Cells That Inhibit T-Cell Responses and Promote TRegs. Blood (2014) 124:750–60. doi: 10.1182/blood-2013-12-546416
27. Moline-Velazquez V, Cuervo H, Vila-Del Sol V, Ortega MC, Clemente D. & de, C.F. Myeloid-Derived Suppressor Cells Limit the Inflammation by Promoting T Lymphocyte Apoptosis in the Spinal Cord of a Murine Model of Multiple Sclerosis. Brain Pathol (2011) 21:678–91. doi: 10.1111/j.1750-3639.2011.00495.x
28. Ochoa AC, Zea AH, Hernandez C, Rodriguez PC. Arginase, Prostaglandins, and Myeloid-Derived Suppressor Cells in Renal Cell Carcinoma. Clin Cancer Res (2007) 13:721s–6s. doi: 10.1158/1078-0432.CCR-06-2197
29. Tacke RS, Lee HC, Goh C, Courtney J, Polyak SJ, Rosen HR, et al. Myeloid Suppressor Cells Induced by Hepatitis C Virus Suppress T-Cell Responses Through the Production of Reactive Oxygen Species. Hepatology (2012) 55:343–53. doi: 10.1002/hep.24700
30. Stroncek DF, Shankar RA, Skubitz KM. The Subcellular Distribution of Myeloid-Related Protein 8 (MRP8) and MRP14 in Human Neutrophils. J Transl Med (2005) 3:36. doi: 10.1186/1479-5876-3-36
31. Hessian PA, Edgeworth J, Hogg N. MRP-8 and MRP-14, Two Abundant Ca(2+)-Binding Proteins of Neutrophils and Monocytes. J Leukoc Biol (1993) 53:197–204.
32. Hobbs JA, May R, Tanousis K, McNeill E, Mathies M, Gebhardt C, et al. Myeloid Cell Function in MRP-14 (S100A9) Null Mice. Mol Cell Biol (2003) 23:2564–76. doi: 10.1128/mcb.23.7.2564-2576.2003
33. De FK, Neill DR, Mathies M, Bangert M, McNeill E, Kadioglu A, et al. A New Protective Role for S100A9 in Regulation of Neutrophil Recruitment During Invasive Pneumococcal Pneumonia. FASEB J (2014) 28:3600–8. doi: 10.1096/fj.13-247460
34. Dessing MC, Tammaro A, Pulskens WP, Teske GJ, Butter LM, Claessen N, et al. The Calcium-Binding Protein Complex S100A8/A9 Has a Crucial Role in Controlling Macrophage-Mediated Renal Repair Following Ischemia/Reperfusion. Kidney Int (2015) 87:85–94. doi: 10.1038/ki.2014.216
35. Hiratsuka S, Watanabe A, Aburatani H, Maru Y. Tumour-Mediated Upregulation of Chemoattractants and Recruitment of Myeloid Cells Predetermines Lung Metastasis. Nat Cell Biol (2006) 8:1369–75. doi: 10.1038/ncb1507
36. Kim HJ, Chang KA, Ha TY, Kim J, Ha S, Shin KY, et al. S100A9 Knockout Decreases the Memory Impairment and Neuropathology in Crossbreed Mice of Tg2576 and S100A9 Knockout Mice Model. PLoS. One (2014) 9:e88924. doi: 10.1371/journal.pone.0088924
37. Zhang X, Wei L, Wang J, Qin Z, Wang J, Lu Y, et al. Suppression Colitis and Colitis-Associated Colon Cancer by Anti-S100a9 Antibody in Mice. Front Immunol (2017) 8:1774. doi: 10.3389/fimmu.2017.01774
38. Dai J, Kumbhare A, Youssef D, McCall CE, El Gazzar M. Intracellular S100A9 Promotes Myeloid-Derived Suppressor Cells During Late Sepsis. Front Immunol (2017) 8:1565. doi: 10.3389/fimmu.2017.01565
39. Hiroshima Y, Hsu K, Tedla N, Wong SW, Chow S, Kawaguchi N, et al. S100A8/A9 and S100A9 Reduce Acute Lung Injury. Immunol Cell Biol (2017) 95:461–72. doi: 10.1038/icb.2017.2
40. Soyfoo MS, Roth J, Vogl T, Pochet R, Decaux G. Phagocyte-specific S100A8/A9 Protein Levels During Disease Exacerbations and Infections in Systemic Lupus Erythematosus. J Rheumatol (2009) 36:2190–4. doi: 10.3899/jrheum.081302
41. Burnet M, Holmes MC. Genetic Investigations of Autoimmune Disease in Mice. Nature (1965) 207:368–71. doi: 10.1038/207368a0
42. Der E, Dimo J, Trigunaite A, Jones J, Jorgensen TN. Gr1+ Cells Suppress T-Dependent Antibody Responses in (NZB x NZW)F1 Male Mice Through Inhibition of T Follicular Helper Cells and Germinal Center Formation. J Immunol (2014) 192:1570–6. doi: 10.4049/jimmunol.1302479
43. Jorgensen TN, Roper E, Thurman JM, Marrack P, Kotzin BL. Type I interferon Signaling is Involved in the Spontaneous Development of Lupus-Like Disease in B6.Nba2 and (B6.Nba2 x NZW)F(1) Mice. Genes Immun (2007) 8:653–62. doi: 10.1038/sj.gene.6364430
44. Dubois EL, Horowitz RE, Demopoulos HB, Teplitz R. NZB/NZW Mice as a Model of Systemic Lupus Erythematosus. JAMA (1966) 195:285–9.
45. Cohen PL, Ziff M. Abnormal Polyclonal B Cell Activation in NZB/NZW F1 Mice. J Immunol (1977) 119:1534–7.
46. Atencio S, Amano H, Izui S, Kotzin BL. Separation of the New Zealand Black Genetic Contribution to Lupus from New Zealand Black Determined Expansions of Marginal Zone B and B1a Cells. J Immunol (2004) 172:4159–66. doi: 10.4049/jimmunol.172.7.4159
47. Mathis KW, Taylor EB, Ryan MJ. Anti-CD3 Antibody Therapy Attenuates the Progression of Hypertension in Female Mice with Systemic Lupus Erythematosus. Pharmacol Res (2017) 120:252–7. doi: 10.1016/j.phrs.2017.04.005
48. Mellors RC. Autoimmune and Immunoproliferative Diseases of NZB/Bl mice and hybrids. Int Rev Exp Pathol (1966) 5:217–52.
49. Puga I, Cols M, Barra CM, He B, Cassis L, Gentile M, et al. B Cell-Helper Neutrophils Stimulate the Diversification and Production of Immunoglobulin in the Marginal Zone of the Spleen. Nat Immunol (2011) 13:170–80. doi: 10.1038/ni.2194
50. Yang Q, Li X, Chen H, Cao Y, Xiao Q, He Y, et al. IRF7 Regulates the Development of Granulocytic Myeloid-Derived Suppressor Cells Through S100A9 Transrepression in Cancer. Oncogene (2017) 36:2969–80. doi: 10.1038/onc.2016.448
51. Liu Z, Bethunaickan R, Huang W, Lodhi U, Solano I, Madaio MP, et al. Interferon-Alpha Accelerates Murine Systemic Lupus Erythematosus in a T Cell-Dependent Manner. Arthritis Rheum (2011) 63:219–29. doi: 10.1002/art.30087
52. Mathian A, Weinberg A, Gallegos M, Banchereau J, Koutouzov S. IFN-Alpha Induces Early Lethal Lupus in Preautoimmune (New Zealand Black x New Zealand White) F1 but Not in BALB/c mice. J Immunol (2005) 174:2499–506. doi: 10.4049/jimmunol.174.5.2499
53. Denny MF, Yalavarthi S, Zhao W, Thacker SG, Anderson M, Sandy AR, et al. A Distinct Subset of Proinflammatory Neutrophils Isolated From Patients With Systemic Lupus Erythematosus Induces Vascular Damage and Synthesizes Type I IFNs. J Immunol (2010) 184:3284–97. doi: 10.4049/jimmunol.0902199
54. Blasius A, Vermi W, Krug A, Facchetti F, Cella M, Colonna M. A Cell-Surface Molecule Selectively Expressed on Murine Natural Interferon-Producing Cells That Blocks Secretion of Interferon-Alpha. Blood (2004) 103:4201–6. doi: 10.1182/blood-2003-09-3108
55. Trigunaite A, Dimo J, Jorgensen TN. Suppressive Effects of Androgens on the Immune System. Cell Immunol (2015) 294:87–94. doi: 10.1016/j.cellimm.2015.02.004
56. Vlachou K, Fanouriakis A, Glymenaki M, Ioannou M, Mastorodemos V, Verginis P, et al. Defects Within the Myeloid Derived Supressor Cells (MDSCx) Compartment May Facilitate Aberrant Immune Responses in Systemic Lupus Erythematosus (SLE). Ann Rheum Dis (2012) 71(suppl. 1):A84–5. doi: 10.1136/annrheumdis-2011-201238.29
57. Huang X, Li J, Dorta-Estremera S, Di Domizio J, Anthony SM, Watowich SS, et al. Neutrophils Regulate Humoral Autoimmunity by Restricting Interferon-γ Production via the Generation of Reactive Oxygen Species. Cell Rep (2015) 12:1120–32. doi: 10.1016/j.celrep.2015.07.021
58. Haas C, Ryffel B, Le Hir M. IFN-Gamma Receptor Deletion Prevents Autoantibody Production and Glomerulonephritis in Lupus-Prone (NZB x NZW)F1 Mice. J Immunol (1998) 160:3713–8.
59. Santiago-Raber ML, Baccala R, Haraldsson KM, Choubey D, Stewart TA, Kono DH, et al. Type-I Interferon Receptor Deficiency Reduces Lupus-like Disease in NZB Mice. J Exp Med (2003) 197(6):777–88. doi: 10.1084/jem.20021996
60. Vogl T, Tenbrock K, Ludwig S, Leukert N, Ehrhardt C, van Zoelen MA, et al. Mrp8 and Mrp14 are Endogenous Activators of Toll-Like Receptor 4, Promoting Lethal, Endotoxin-Induced Shock. Nat Med (2007) 13:1042–9. doi: 10.1038/nm1638
61. Wache C, Klein M, Ostergaard C, Angele B, Häcker H, Pfister HW, et al. Myeloid-Related Protein 14 Promotes Inflammation and Injury in Meningitis. J Infect Dis (2015) 212:247–57. doi: 10.1093/infdis/jiv028
62. Huang M, Wu R, Chen L, Peng Q, Li S, Zhang Y, et al. S100A9 Regulates MDSCs-Mediated Immune Suppression via the RAGE and TLR4 Signaling Pathways in Colorectal Carcinoma. Front Immunol (2019) 10:2243. doi: 10.3389/fimmu.2019.02243
63. Vogl T, Stratis A, Wixler V, Völler T, Thurainayagam S, Jorch SK, et al. Autoinhibitory Regulation of S100A8/S100A9 Alarmin Activity Locally Restricts Sterile Inflammation. J Clin Invest (2018) 128:1852–66. doi: 10.1172/JCI89867
64. Schenten V, Plançon S, Jung N, Hann J, Bueb JL, Bréchard S, et al. Secretion of the Phosphorylated Form of S100A9 from Neutrophils Is Essential for the Proinflammatory Functions of Extracellular S100A8/A9. Front Immunol (2018) 9:447. doi: 10.3389/fimmu.2018.00447
65. Lim SY, Raftery MJ, Goyette J, Geczy CL. S-Glutathionylation Regulates Inflammatory Activities of S100A9. J Biol Chem (2010) 285:14377–88. doi: 10.1074/jbc.M109.075242
66. Lim SY, Raftery M, Cai H, Hsu K, Yan WX, Hseih HL, et al. S-Nitrosylated S100A8: Novel Anti-Inflammatory Properties. J Immunol (2008) 181:5627–36. doi: 10.4049/jimmunol.181.8.5627
67. Jung N, Schenten V, Bueb JL, Tolle F, Bréchard S. miRNAs Regulate Cytokine Secretion Induced by Phosphorylated S100A8/A9 in Neutrophils. Int J Mol Sci (2019) 20:5699. doi: 10.3390/ijms20225699
68. Chen JQ, Papp G, Póliska S, Szabó K, Tarr T, Bálint BL, et al. MicroRNA Expression profiles Identify Disease-Specific Alterations in Systemic Lupus Erythematosus and Primary Sjögren’s Syndrome. PloS One (2017) 12:e0174585. doi: 10.1371/journal.pone.0174585
69. Su YJ, Lin IC, Wang L, Lu CH, Huang YL, Kuo HC. Next Generation Sequencing Identifies miRNA-Based Biomarker Panel for Lupus Nephritis. Oncotarget (2018) 9:27911–9. doi: 10.18632/oncotarget.25575
70. Pavon EJ, Garcia-Rodriguez S, Zumaquero E, Perandres-Lopez R, Rosal-Vela A, Lario A, et al. Increased Expression and Phosphorylation of the Two S100A9 Isoforms in Mononuclear Cells From Patients With Systemic Lupus Erythematosus: A Proteomic Signature for Circulating Low-Density Granulocytes. J Proteomics (2012) 75:1778–91. doi: 10.1016/j.jprot.2011.12.020
71. Chakrabarti S, Lekontseva O, Peters A, Davidge ST. 17beta-Estradiol Induces Protein S-Nitrosylation in the Endothelium. Cardiovasc Res (2010) 85:796–805. doi: 10.1093/cvr/cvp368
Keywords: lupus, MDSC, S100A9, autoantibody, type I interferon, sex specific effects, mouse model, proteinuria
Citation: Davison LM, Alberto AA, Dand HA, Keller EJ, Patt M, Khan A, Dvorina N, White A, Sakurai N, Liegl LN, Vogl T and Jorgensen TN (2021) S100a9 Protects Male Lupus-Prone NZBWF1 Mice From Disease Development. Front. Immunol. 12:681503. doi: 10.3389/fimmu.2021.681503
Received: 16 March 2021; Accepted: 01 June 2021;
Published: 17 June 2021.
Edited by:
Carlo Perricone, Sapienza University of Rome, ItalyCopyright © 2021 Davison, Alberto, Dand, Keller, Patt, Khan, Dvorina, White, Sakurai, Liegl, Vogl and Jorgensen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Trine N. Jorgensen, am9yZ2VudEBjY2Yub3Jn
†These authors have contributed equally to this work