Calming the Cytokine Storm - Splenic Ultrasound for Treating Inflammatory Disorders and Potentially COVID-19

Non-invasive Splenic Ultrasound Reduces TNF in Humans

To investigate how non-invasive ultrasound stimulation of the spleen affects cytokine levels directly in human subjects, healthy participants were recruited into a study designed to investigate effect of different ultrasound parameters and locations of stimulation in the spleen on various immune and metabolic responses during acute stimulation (feasibility study described at clinicaltrials.gov: NCT03548116). Due to the urgency of identifying potential treatment options for COVID-19 patients, data pertaining specifically to ultrasound stimulation of the hilum of the spleen with the most effective intensity parameter in the study (corresponding to one cohort from the study) are presented in this paper.

Participants received pulsed ultrasound stimulation of the spleen for three minutes (2.2 MHz, 1 Hz pulse repetition frequency, 290.4 mW/cm² ISPTA) using the GE LOGIQ E9 ultrasound system with the C1-6 probe. An ultrasonographer first identified and targeted the hilum of the spleen, then delivered the pulsed ultrasound stimulus using a modified elastography software setting in the system. Sham stimulation control subjects underwent the same procedure as the stimulated group, including transducer placement and movement of the ultrasound probe to the proper splenic location trajectory, but with no ultrasound energy being transmitted to the spleen from the transducer. Since healthy subjects have low serum TNF levels, we revealed changes in inflammatory status using whole blood challenged with lipopolysaccharide (LPS) ex vivo. This ex vivo cytokine production assay has been previously used to verify splenic neuroimmune activation using implanted vagus nerve stimulators and is a well-established method for non-invasively assessing activation of neuroimmune pathways (18, 41-43). A reduction in LPS-induced TNF production in blood sampled from ultrasound stimulated subjects was tested at multiple LPS concentrations (Fig. 1A) and significant reduction is observed when comparing the ultrasound stimulated group to the non-stimulated controls (Fig. 1B; p-value = 0.006 using Wilcoxon unpaired rank-sum test) at the optimal LPS concentration (1 ng/mL; Fig.1A). These data are consistent with earlier reports using implanted nerve stimulators or pharmaceutical activators of the splenic pathway (18, 42). They are also consistent with the magnitude of TNF reduction using similar ultrasound stimulation parameters in a previous study activating the neuroimmune pathway in a rodent model of endotoxemia (31). Further supporting the human results, we also demonstrated that splenic ultrasound stimulation in a rodent model achieves a reduction across a number of pro-inflammatory cytokines beyond TNF, including IL-6, IFNγ, IL-1β, IL-1α, and IL-12 in splenic lysates (Fig. S1), which is consistent with previous reports using implanted vagus nerve stimulators (18, 44).

• Download figure
• Open in new tab

Figure 1. Non-invasive ultrasound stimulation of the spleen reduces TNF production in human subjects.

A) Whole blood from all subjects and samples was incubated with LPS at multiple concentrations. TNF response after LPS incubation is shown for one ultrasound stimulated subject from samples collected before (0 hours baseline) and after (2 hours post stimulation) ultrasound application. The optimal dose of LPS for measuring the immunomodulatory response to ultrasound was found to be at 1 ng/mL, and this is used for comparing the TNF response between the ultrasound stimulated versus control groups. B) Ultrasound stimulation of the spleen decreases whole-blood LPS-induced TNF release, where blood was obtained from 18 healthy human subjects prior to and then two hours following splenic ultrasound application (stimulated: n=9, control: n=9). The data shown were taken from the 1 ng/ml LPS concentration for each dose response curve for all subjects and samples. p-value was computed using the unpaired Wilcoxon rank-sum test.

Pro-inflammatory Cytokines and Signaling Pathways are Downregulated in Circulating Monocytes with Ultrasound Treatment

Previous animal research demonstrated that ultrasound can significantly reduce inflammation and improve clinical measures in a chronic inflammatory arthritis model (32). Expanding upon the findings in healthy human subjects in the previous section, we initiated a controlled, randomized clinical trial of splenic ultrasound treatment in RA patients (clinical study described at clinicaltrials.gov: NCT03690466). Participants received ultrasound stimulation to the spleen with the goal of lowering the over-active inflammatory response that occurs during RA flare-ups. Up to 20 RA patients are to be recruited into the study (enrollment is ongoing). Due to the urgency of identifying treatment options for this COVID-19 pandemic, data from single-cell RNA sequencing (scRNA-seq) in peripheral blood mononuclear cells (PBMCs) that has been collected thus far in five ultrasound-treated patients have been analyzed ahead of completion of the study for inclusion into this paper. PBMCs were isolated from whole blood taken from each research participant before ultrasound treatment (day 0, pre-stimulation) and after two weeks of daily 30-minute sessions of ultrasound to the spleen (day 14, post-stimulation). From these 10 samples (subjects = 5, timepoints = 2), a total of 53,343 PBMCs were successfully sequenced from single cells (Table S1). In silico analysis identified 13 cell-types based on gene expression patterns of standardly used marker genes (Fig. S2). Participants received stimulation with a 1 MHz transducer (1.2 W/cm²; SoundCare Plus, Roscoe Medical) that was moved continuously across a 5-inch by 5-inch square, centered on the spleen, as identified by an ultrasonographer. Further details on the clinical study design is presented in the Materials and Methods section.

We aimed to identify molecular changes that occur with splenic ultrasound in humans to better understand how this ultrasound treatment can drive an anti-inflammatory response. We observed reduced monocyte expression of key pro-inflammatory cytokine genes encoding for IL-1β and IL-8 from all five subjects post-stimulation compared with pre-stimulation (Fig. 2A, p-values given for transcripts showing significance of p<0.05). Although this result was determined by analyzing all monocytes, reduction of these two markers appears primarily in CD14+ monocytes as this is the cell type showing predominant expression of those genes (Fig. 2B). The same analysis in CD8+ T cells show a decrease in transcripts that encode for IFNγ (Fig. 2A). The changes in transcript level of these cytokines can only be detected when comparing within the same cell type before and after treatment, within the cell that expresses the cytokine (Fig. S2). Consistent with known expression patterns of IL1B and CXCL8 transcripts (encoding for IL-1β and IL-8), we see that these transcripts are only suppressed in the monocytes where they are primarily detected. Likewise, IFNG transcripts (encoding for IFNγ) are detected and found to be reduced in CD8+ T cells (Fig. S2). In our dataset, other cytokines like TNF and IL-6 did not show a significant change in monocytes; however, IL-6 transcripts were detected in only 66 of the 53,343 cells. Therefore, we cannot definitively conclude if there is or is not a change in transcript levels encoding IL-6 with ultrasound treatment.

• Download figure
• Open in new tab

Figure 2. Splenic ultrasound treatment reduces cytokines and downstream pathways in-vivo in RA patients.

Single cell RNA sequencing of PBMCs from RA patient receiving 14 days of 30-minute daily ultrasound stimulation of the spleen. Samples were collected pre-(day 0) and post-stimulation (day 14) with ultrasound. Dot plot shows percent of CD14+ and CD16+ monocytes (top) or CD8+ T cells (bottom) that express the transcript (labeled on the left) represented by size, and scaled expression level represented by color (obtained by subtracting the mean and dividing by the standard deviation). Key pro-inflammatory cytokines are shown with unadjusted p-values labeled for those transcripts significantly downregulated (p-value <0.05). B) UMAP showing top two dimensions (principal components) contributing to variability within the dataset. UMAP diagrams showing CD14+ and CD16+ monocytes represented by orange and green dots, respectively. IL-1β and CXCL8 expression on day 0 and day 14 are also mapped on these UMAP diagrams. C) Heatmap displaying the average relative expression (for all cells in specified group) of downregulated transcripts (genes listed on the right, p-value <0.05). Color denotes relative expression for post-ultrasound levels compared to (divided by) pre-ultrasound levels for each subject (represented by columns, n = 5, participant A-E shown left to right) in CD14+ cells or CD16+ cells. Genes encoding for cytokine receptors, components of the JAK/STAT pathway, MAP kinases, and transcription factor (TF) are grouped and labeled on the left of the plot.

Since monocytes play a critical role in the innate immune system and the inflammatory response, we further focused our analyses on transcriptional changes that occur in circulating monocytes after ultrasound treatment. Differential expression analysis of all monocytes across the five subjects, pre- and post-ultrasound stimulation, shows 938 differentially expressed genes (p-value<0.05) with 841 of those genes being down-regulated. These results reveal an overall shift in reduction of gene transcripts in monocytes with ultrasound treatment. Additionally, this reduction in transcript levels did not appear as a result of sequencing differences, as we observe no obvious difference in the number of cells sequenced, the number of sequencing reads, or number of average genes detected per cell between monocytes taken from day 0 as compared to day 14 (Table S1, ‘CD14 Mono’ and ‘CD16 Mono’). Although these gene transcripts are defined as differentially expressed when comparing monocytes from all subjects (Fig. 2A), we also observe that for most genes, this down-regulation is consistent across subjects (Fig. 2C; columns represent cell-type averages for CD14+ and CD16+ cells for each subject that are similar across columns for many of the genes).

These results are quite consistent given the high variability or inhomogeneity in age, sex, physical characteristics, disease properties and past treatments across subjects (participants A-E in Table S2 correspond to columns 1-5 in Fig. 2C, respectively). Within these down-regulated genes in monocytes, we also observe that many cytokine receptors, such as those for TNF, IL-6, IL-17, IL-13, and IFNγ were reduced with ultrasound treatment (Fig. 2C). Furthermore, components of signaling pathways downstream of these cytokines are also down regulated, including components of the JAK/STAT pathway, MAP kinases, and pro-inflammatory transcription factor NFκB (Fig. 2C). Together, these findings demonstrate a consistent and robust reduction in transcripts encoding key pro-inflammatory cytokines and cellular signaling molecules downstream of their receptors after 14 days of daily splenic ultrasound treatment across all five subjects.

Suppression of Multiple Pathways that Contribute to the Pro-inflammatory Response

Cytokines trigger multiple cellular pathways, which then contribute to systemic inflammation. Our transcriptomic data reveal that many genes known to be regulated by IFNγ and NFκB are subsequently downregulated with ultrasound treatment, further promoting a cascade of suppression across multiple pro-inflammatory pathways with ultrasound treatment (Fig. 3A). We sought to determine if the 840 down-regulated genes in monocytes showed statistical enrichment for functional Gene Ontology (GO)terms associated with inflammation. The GO system of classification assigns genes (for differentially expressed gene list) into bins based on their functional characteristics. Gene sets can then be analyzed for GO terms that are statistically over-represented among that set, allowing for a better understanding of which biological processes (e.g., inflammation) may be changed with treatment. Out of the top statistically enriched GO terms for the down-regulated genes, we observe many that contribute to pro-inflammatory cellular processes. Genes suppressed by ultrasound show enrichment of GO terms including “inflammatory response” (p-value = 1.50 E-5, Table S3), “cytokine-mediated signaling pathway,” and “cellular response to cytokine stimulus” (p-values = 9.7 E-5 and 6.7E-6, Fig. 3B). We also observe enriched GO terms associated with positive regulation of TNF, IL-6, and IL-8 among the down-regulated genes, suggesting that inducing signals for TNF or IL-6 are also suppressed, despite the absences of significant change in TNF or IL-6 transcript levels (Fig. 3B). These GO term results show that many gene transcripts encoding for cellular inflammatory functions go down after ultrasound treatment.

• Download figure
• Open in new tab

Figure 3. Reduction of pro-inflammatory pathway components in monocytes with ultrasound treatment.

Volcano plot showing all detected genes in monocytes, highlighting differentially expressed genes (p-value <0.05) and labeled IFNγ-regulated genes that are down regulated, as well as NFκB-regulated genes that are down regulated as characterized by KEGG pathway gene lists. B) Gene ontology (GO) analysis using Bioconductor package ‘topGO’, determining functional GO terms significantly enriched within downregulated genes in monocytes. Enriched GO terms associated with cytokine production are shown, where bars represent number of significant genes in a list relative to number of genes that would be expected by chance, and with number of genes contributing to each GO term represented by the shade of blue. p-values for each term shown is from Fisher exact test.

Down-regulation of Components Involved in Monocyte Migration

Within circulating monocytes, we observe ultrasound-treatment associated GO term enrichment of cellular components associated with IL-8 functionality (Fig. 3B), as well as a robust reduction in IL-8 across subjects (Fig. 4A, columns represent subject cell type averages). IL-8 is a chemokine (chemotactic factor), which stimulates neutrophil movement towards inflammatory sites. Therefore, we further assessed if ultrasound influences other signaling molecules involved with cellular migration (e.g., to sites of inflammation). GO term analysis of down-regulated genes in CD14+ monocytes show enrichment of “positive regulation of cell migration” among the top 15 significantly enriched GO terms. Other GO terms, such as “lamellipodium assembly” (facilitates cell mobility), “ephrin receptor signaling” (promotes cell adhesion), and “actin nucleation” (involved with cell locomotion), all suggest that ultrasound treatment induces down regulation of genes involved with cell migration (Fig. 4B). Furthermore, genes that more specifically contribute to “monocyte chemotaxis” are down-regulated post-ultrasound stimulation in CD14+ monocytes and many of these genes are consistently reduced across all five subjects (Fig. 4A), even with quite variable or inhomogeneous patient characteristics (Table S2).

• Download figure
• Open in new tab

Figure 4. Ultrasound results in downregulation of cellular migration components in CD14+ monocytes.

A) Heatmap displaying the average relative expression (for all cells in specified group) of significantly downregulated transcripts (genes listed on the right) for post-ultrasound compared to (divided by) pre-ultrasound levels for each subject (represented by columns, n = 5) in CD14+ cells. Genes encoding components involved in monocyte chemotaxis and cell migration are grouped and labeled to the left. B) Enriched GO terms within downregulated genes in CD14+ cells are shown, where bars represent the number of significant genes in the list relative to the number of genes that would be expected by chance, and with the number of genes contributing to each term represented by the shade of blue.. p-value for each term shown is based on Fisher exact test.

Immune Cell Repertoire Unchanged with Moderate Degree of Suppression

We sought to determine how ultrasound affects the circulating immune cell repertoire. In comparing the transcriptional profiles (based on principal component analysis and UMAP dimensional reduction) of all circulating peripheral blood mononuclear cells pre- and post-ultrasound stimulation, we observe remarkable consistency (Fig. 5A). Also, the calculated relative percent of each cell-type for each subject is displayed, and based on a paired Wilcoxon rank sum test for pre- and post-stimulation for each cell type, we observe no significant difference in the relative proportions of percent contribution of any of the 13 cell types (Fig. 5B, statistics shown in Table S4). Additionally, when comparing the number of cells sequenced for each cell type per participant there is no significant difference between timepoints (i.e. no technical bias, Table S4).

• Download figure
• Open in new tab

Figure 5. Ultrasound provides a moderate effect on the immune system.

A) UMAP showing top two dimensions (principal components) contributing to variability within the dataset. All cells are shown, consisting of 53,343 cells from 10 samples (subjects = 5, timepoints = 2) with cell type and timepoint represented by color showing transcriptional signatures and relative number of cells is consistent between timepoints. B) Box plot showing median percent of each cell type (subjects, n =5) captured before ultrasound stimulation (day 0, orange, open squares) and after ultrasound stimulation (day 14, blue, open triangles) with filled in circles representing outliers. No statistical difference was detected using a paired Wilcoxon ranked-sum test (Table S4). C) Percent down regulation range is shown after 14 days of ultrasound stimulation for each gene (represented by dots) in each functional category described previously.

We also determined what degree of suppression occurred after ultrasound stimulation for genes involved in inflammation (genes shown in Fig. 2C, 3A and 4A). Most genes show a ∼5-15% reduction in transcript levels, suggesting a moderate degree of suppression, with IL-8 and genes involved in cell migration showing the highest degree of suppression (>15%, Fig. 5C). These findings suggest that spleen-targeted ultrasound can potentially provide a way to blunt the elevated release of circulating cytokines and chemokines in inflammatory disorders or viral infections, without disrupting the overall immune cell repertoire.

B Cells Retain Antibody Production Abilities

One concern in treating acute hyperinflammation is that the therapy may suppress the adaptive immune response, and the ability of the immune system to provide pathogen clearance and/or protection from co-infection with a second pathogen. Therefore, we assessed whether non-invasive splenic ultrasound treatment suppresses the adaptive immune response or more specifically antibody production. Our transcriptomic data do not suggest a suppression of the adaptive immune response. In fact, we observe that among genes upregulated in B cells, there is an enrichment for GO terms associated with “B cell activation” (Fig. 6A). Genes associated with this GO term include transcripts that encode for IgA, IgG, and IgM antibodies. Interestingly, these gene transcripts are upregulated in some of the ultrasound-treated subjects (Fig. 6B), though with greater variability in results across subjects than was observed for suppression of cytokines and chemokines in monocytes. For example, we observe that for average transcript levels across B cells, transcripts encoding IgG are only upregulated in 2 or 3 subjects (IGHG2 and IGHG4, respectively) out of the 5 subjects. Regardless of the variability across subjects, these data support that ultrasound treatment does not generally inhibit B cell activation and antibody production across patients, and can even potentially enhance antibody production; thus, ultrasound is not expected to compromise the adaptive immune response.

• Download figure
• Open in new tab

Figure 6. Ultrasound treatment does not reduce adaptive immune response.

A) Enriched GO terms within upregulated genes in all B cells are shown, where bars represent number of significant genes in list relative to number of genes that would be expected by chance, and with number of genes contributing to each term represented by red hue. p-value for each term shown is based on Fisher exact test. B) Heatmap displaying genes that contribute to the GO term “positive regulation of B cell activation” within B cells. Each column represents relative average expression per subject for pre-stimulation (left 5 columns) and post-stimulation (right 5 columns) for each gene (labeled to the right). C) A rodent model of endotoxemia was used to further characterize the effects of splenic ultrasound stimulation on antibody production (further details of experiments are provided in Materials and Methods section). Box plot showing median levels of IgM (left) or IgG (right) concentrations in serum (open circles represent individual animals) for naïve animals (no LPS + sham), animals exposed to splenic ultrasound stimulation with no endotoxin (no LPS + ultrasound), animals exposed to the endotoxin but receiving sham stimulation (LPS + sham), and animals receiving both LPS exposure and splenic ultrasound stimulation (LPS + ultrasound). Animals receiving LPS + ultrasound showed significant increases in antibody production when compared to each of the other groups using an unpaired Wilcoxon ranked-sum test for IgM (no LPS + sham, p = 0.00145; no LPS + ultrasound, p = 0.00155; LPS + sham, p =0.01757) and for IgG (no LPS + sham, p =0.02930; no LPS + ultrasound, p = 0.00404; LPS + sham, p =0.04298). No other comparisons using this test showed significance, p < 0.05.

To further investigate the effects of splenic ultrasound stimulation on adaptive immune response and antibody production, as well as the variability we observed across human subjects for transcript levels that encode for IgA, IgG, and IgM antibodies, we performed additional experiments using prophylactic ultrasound stimulation (i.e., ultrasound stimulation prior to antigen exposure) in a rodent model of endotoxemia (31). We observe that splenic ultrasound stimulation prior to endotoxin exposure (three minutes daily, once a day for 3 days followed by LPS injection on the fourth day) results in a significant increase in IgM and IgG antibody production 24 hours after LPS exposure compared to sham and/or non-LPS controls (p-values listed in Fig. 6C description). Furthermore, this enhanced antibody output is dependent on the presence of the toxin/antigen, as there is no observed ultrasound effect on antibody production in the absence of LPS exposure. The increase in antibody transcript levels we observed in some of the human subjects, and not in the others, may be attributed to those participants also having encountered an antigen, such as a cold or other infection that was not identified during the study. This explanation is further supported by a recent report demonstrating that brain-spleen neural pathways are capable of modulating splenic plasma cells differentiation upon presentation of novel T Cell dependent antigens (45).

Splenic ultrasound in healthy subjects: Design of clinical study

For the splenic ultrasound study involving healthy human subjects (clinicaltrials.gov: NCT03548116), we provide the data from the two cohorts most relevant to the planned DARPA sponsored COVID study: the sham control cohort (i.e., no ultrasound, with only mechanical pressure over the spleen) and a splenic ultrasound stimulation cohort (i.e., 290.4 mW/cm² ISPTA or 1.4 mechanical index (MI) using the shear wave elastography setting on GE LOGIQ E9) that insonified a single site (i.e., the splenic hilum; as identified by a trained ultrasonographer performing the stimulation). The trial was carried out in accordance with International Conference on Harmonisation Good Clinical Practice (ICH GCP) and the United States Code of Federal Regulations (CFR) applicable to clinical studies (45 CFR Part 46, 21 CFR Part 50, 21 CFR Part 56, and 21 CFR Part 812). The protocol, informed consent form, recruitment materials, and all participant materials were submitted and approved by the Institutional Review Board (IRB) at Northwell Health.

Inclusion and exclusion criteria

To be eligible to participate in the study, an individual must have met the following criteria: aged between 18 and 45 years, be without physical disability or conditions that may make them incapable of undergoing the study procedures or otherwise place them at greater harm, be without significant past medical or surgical histories that would render them at a greater risk of harm, be considered English proficient due to the study requirement to follow verbal commands during the ultrasound session, be considered active as assessed by type of activity (e.g., walking or running) and number of hours a week performing the various activities, be able to attend all study visits at approximately the same time of day, be able to comprehend the study goals and procedures, and able to provide informed consent for participation. See Table S5 for subject characteristics used in the cytokine analysis shown in Fig. 1.

Study timeline

Subjects underwent a screening visit 1 to 11 days prior to baseline visit to assess eligibility to participate in the study. Eligible individuals that agreed to participate and provided written informed consent then underwent a physical and neurological examination. Women of childbearing potential were asked to provide a urine sample for pregnancy testing, and approximately 21 mL of blood was drawn. On the first protocol visit (day 0), participants again underwent a physical and neurological examination, and medical history review. A baseline blood draw was then collected for cytokine and blood chemistry analysis (∼35 mL). The blood draw was performed under sterile conditions using standard venipuncture techniques. Individuals then underwent the ultrasound procedure based on random assignment to one of 7 groups, including the two groups reported herein (i.e., the sham control group, and a stimulated group that received 100% power stimulus applied at a single splenic target location). For ultrasound stimulus delivery, individuals were asked to lie in the right lateral recumbent position with their arms above their heads to expose the splenic region of the abdomen. After ultrasound gel was applied to the region, the ultrasound probe was placed on the participant’s abdomen. This procedure was performed for both the sham control group (n=9) and ultrasound stimulated group (n=9). Note that one subject from each group was excluded, in which 10 subjects were initially recruited to the study for each group (see Table S5).

Stimulation paradigm

In the ultrasound stimulated cohort, the subject had the dimensions of their spleen assessed (using the B mode imaging setting on the GE LOGIQ E9), and the probe (C1-6) was positioned over the splenic hilum. Stimulation was performed using a modified elastography setting on the LOGIQ E9 operated in two modes: an imaging/targeting mode with short ultrasound pulses (on the order of 1 microsecond) and relatively low acoustic power levels (on the order of 10 mW/cm²), and a stimulation mode (from the clinical shear wave elastography setting) using longer ultrasound pulses (on the order of 1 ms) and higher acoustic power levels (290.4 mW/cm²). Stimulation was performed in twelve 15-second long epochs separated by 15 seconds each, and a total stimulation duration of 3 minutes. The 15 second epochs were performed during breath holds, during which an ultrasonographer was holding the probe in position above the spleen. Before each stimulation period, the ultrasonographer checked the location of the probe using the imaging mode to ensure that all stimulation pulses were on target. The focal landmark used for targeting and stimulation was the hilum of the spleen in the group included in the presented analysis. For the sham control group, the movements and probe contact for the ultrasound procedure were performed by the ultrasonographer, and the total duration of the procedure remained the same as the stimulated group, except that no ultrasound power was applied (i.e., ultrasound output via the probe was turned off). A post-ultrasound blood draw was again performed, and the 2-hour blood draw time-point data is reported herein.

Blood analysis procedures

Prior to the ex vivo cytokine test, the LPS was sonicated for 30 minutes. A 10 mL blood tube (heparin; green cap) was brought into a tissue culture/sterile room and blood was transferred to a 50 mL conical tube. A blinded researcher then diluted the LPS to the test concentrations in 15 pre-labeled tubes containing replicates of the different concentrations of LPS used for blood exposure (dilutions were made to achieve 0, 0.1, 1, and 10 ng/mL of LPS when mixed with blood). Upon mixing the LPS with blood, the final tubes were capped and placed on a rocker within an incubator and the samples were kept at 37 degrees Celsius for 4 hours during incubation. After the incubation, the samples were removed from the incubator, centrifuged at 6000 RPM for 5 minutes, and the plasma supernatant was transferred (via pipette) for storage (frozen prior to analysis). TNF analysis was performed in triplicate from each sample using DuoSet ELISA kit (R&D Systems), as per manufacturer’s instructions.

Splenic ultrasound in RA patients: Design of clinical study

This clinical study is a controlled, randomized, double-blinded trial of short-term (14 days) treatment and is listed at clinicaltrials.gov (NCT03690466) for which a portion of the results are presented in this paper. The study is still ongoing, but the RNAseq data obtained thus far in five ultrasound stimulated patients is presented in this paper for its relevance in revealing the use of splenic ultrasound for potentially treating COVID-19 patients. The protocol, informed consent form, recruitment materials, and all participant materials were submitted and approved by the University of Minnesota’s Institutional Review Board (IRB) and monitored by CTSI in accordance with its institutionally approved monitoring plan. The overall objective of the clinical trial is to demonstrate safety and efficacy of spleen ultrasound stimulation in the treatment of RA in 20 participants. For this ongoing study, 17 participants have been recruited, consented, and enrolled with all 17 patients completing the full study. The 17th participant completed the study, however the last ‘in-clinic’ visit was performed using clinician video conferencing due to the COVID-19 pandemic, while in-person blood sample collection was still able to be performed. Further details for this study can be found on clincialtrials.gov and will be presented in a future publication with complete outcome analyses. The description below only includes relevant methods and RNAseq analyses presented specifically in this paper.

Inclusion and exclusion criteria

To be eligible to participate in the study, an individual must have met the following criteria: 18 years of age or older, carried a diagnosis of seropositive (rheumatoid factor- or cyclic citrullinated peptide antibody-positive) rheumatoid arthritis (RA), and exhibited symptoms or signs of inadequate disease control (either modified HAQ score >0.3 or DAS-28-CRP >3.2). Exclusion criteria included active bacterial or viral infection, pregnancy, malignancy, or inability to provide daily self-care. The patient medical history and characteristics for each of the five subjects presented in this paper is provided in Table S2.

Study timeline

Participants were assessed during in-clinic visits on days 0 (before treatment began), 3, 7, 10, 14, and 21. Minor flexibility was allowed for scheduling clinic visits outside of holidays for special circumstances. On the first visit, the spleen was imaged by an ultrasonographer and the center of the spleen trajectory was marked on the participant’s skin. Participants were instructed on how to self-administer ultrasound treatment to the spleen. The procedure involved sweeping the ultrasound wand continuously across a 5-inch by 5-inch square centered on the spleen mark placed by the ultrasonographer and ensuring full contact was made with the skin during the daily 30 minute stimulation period. Researchers were able to visually confirm that daily treatments were administered through daily online video sessions with the patients. These video sessions took place at about the same time each day for the 14-day treatment period. With each in-clinic visit, participants received a joint examination, peripheral joint ultrasound, and blood draw, as well as completing several physical evaluations and questionnaires. Full study details including analysis of all patients and assessments will be presented in a future publication.

Treatment

After an in-clinic training session on Day 0, transcutaneous ultrasound was administered to the spleen for 30 minutes daily for 2 weeks (14 days total) via a portable device (1.2 W/cm²; SoundCare Plus, Roscoe Medical) in the patient’s home. The patients were randomized (1:1) to a treatment group or a control (sham) group, in which the latter received the exact same evaluations and stimulation paradigm as the treatment group, except that the device did not deliver any energy from the ultrasound transducer. RA patients and clinical assessors were blinded as to which patients received ultrasound or sham treatment.

RNA-sequencing and statistical analyses

4mL of blood was collected in a green top LiHep tube from each ultrasound stimulated participant (n=5) on day 0 (before the first splenic ultrasound stimulation was administered) and on day 14 (after the final treatment session). The samples were processed the same day as collection. Peripheral blood mononuclear cells were isolated from whole blood using density gradient centrifugation via SepMate tubes (Stemcell Technologies) and following manufacturer’s instructions. PBMCs were then gently frozen and stored at −80 °C as per 10X Genomics demonstrated protocols GC000039. Following this protocol, once all samples had been collected and frozen, samples were thawed quickly in 20% FBS in PBS, resuspended in 10% FBS, and strained before proceeding to 10X genomics single cell protocol. Samples were sequenced using UMI-based approach from10X Genomics, as previously described (32) and mapped to the human genome.

Transcriptomic analysis was performed using R statistical software, and primary filtering and normalization of each sample was done using Seurat package as previously described (32, 55, 56). Samples included two timepoints from five subjects and these individual samples were merged for cell type assignment and differential expression analysis. Cell types were assigned using standard marker genes (see Fig. S2) as previously described (57). Differential expression was determined between timepoints for all cells in each cell type using a negative binomial generalized linear model. GO term analysis was performed using the org.Hs.eg.db and topGO packages. GO term enrichment was assessed using the ‘elim’ algorithm with a Fisher exact test cutoff of 0.01, in which the top 40 significant GO terms in monocytes are shown in Table S3 (58).

Pre-clinical tests in LPS/endotoxin rodent model

Pre-clinical tests in rodents were performed as previously described (31). Experiments were performed under protocols approved by the Institutional Animal Care and Use Committee of GE Research. Briefly, a VIVID E9 (GE) ultrasound system was used to perform an ultrasound scan and locate the spleen. A HIFU transducer and system was then positioned on the target area, and a separate ultrasound scan performed using a smaller probe (3S, GE) that was placed within the opening of the HIFU transducer was used to verify alignment of the ultrasound beam with the spleen. Adult Sprague-Dawley rats that were 8-12 weeks old (250-300 g; Charles River Labs) were housed at 25 degrees Celsius on a 12-h light/dark cycle and acclimatized for 1 week, with handling before experiments to minimize the potential confounding measures due to stress response. Water and regular rodent chow were available ad libitum. LPS (0111:B4, Sigma Aldrich) was used to produce a significant state of inflammation in the naïve adult rodents. LPS was administered to animals in the amount of 10 mg/kg, which corresponds to a LD75 dose via intraperitoneal injection.

Treatment

For ultrasound stimulation, animals were anesthetized with 2-4% isoflurane, and the animals were laid on a water circulating warming pad to prevent hyperthermia during the procedures. The region designated for ultrasound stimulation was shaved with a disposable razor and animal clippers prior to stimulation. Ultrasound was applied to the designated area above the spleen using the probe and system described and calibrated from a previous study (31). An ultrasound stimulus using the previously found optimal stimulation parameters (1.1 MHz, 0.5 ms pulse repetition frequency, 136.36 microsecond pulse length, 568 mW/cm² ISPTA) was then applied for 2 minutes, and LPS was immediately injected post-stimulation. Animals were allowed to incubate under anesthesia for one hour, and after incubation the animal was euthanized, and tissue and blood samples were collected. For experiments testing antibody production, an ultrasound stimulus was applied daily (in the morning) for three days prior to LPS injection. On the fourth day, following the three consecutive treatment days, LPS (10mg/kg) was injected intraperitoneally. The animals were then incubated for 60 minutes in accordance with previous studies, and tissue and blood were harvested for analysis of cytokines and antibody production.

Tissue collection

An incision was made starting at the base of the peritoneal cavity extending up and through to the pleural cavity. The spleen was rapidly removed and homogenized in a solution of phosphate-buffered saline, containing phosphatase (0.2-mM phenylmethylsulfonyl fluoride, 5-µg/mL aprotinin, 1-mM benzamidine, 1-mM sodium orthovanadate, and 2-µM cantharidin) and protease (1-µL to 20 mg of tissue as per Roche Diagnostics) inhibitors. A targeted final concentration of 0.2-g tissue per mL PBS solution was applied in all samples. After collection of the whole blood, we allowed the blood to clot undisturbed at room temperature for 15–30 minutes. Samples were centrifuged at 1,000–2,000x g for 10 minutes in a refrigerated centrifuge with the resulting supernatant being collected for antibody testing. Samples were then stored at −80 °C until analysis. Splenic tissue lysate was analyzed for cytokine concentrations using ELISA kits for quantifying TNF alpha (Abcam, Ab100785), IFN gamma (Abcam, Ab239425), IL-13 (Abcam, Ab100766), IL-12 (LSBio, LS-F34357), IL-10 (Abcam, Ab100765), IL-2 (Abcam, Ab221834), IL-1 beta (Abcam, Ab100768), IL-1 alpha (Abcam, Ab113350), IL-6 (Abcam, Ab234570) and IL-4 (Abcam, Ab100771) as per manufacturer’s instructions for tissue samples. Serum was tested using ELISA kits for quantifying antibody production including IgG (Abcam, Ab189578), and IgM (Abcam, Ab215085) as per manufacturer’s instructions.