Prerpints.org logo
Preprint
Article

Electroacupuncture Reduces Fibromyalgia Pain via Neuronal/Microglial Inactivation and Toll-Like Receptor 4 in the Mouse Brain: Precise Interpretation of Chemogenetics

Altmetrics

Downloads

84

Views

40

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

24 January 2024

Posted:

24 January 2024

You are already at the latest version

Alerts
Abstract
Fibromyalgia (FM) is a complex, chronic, widespread pain syndrome that can cause significant health and economic burden. Emerging evidence has shown that neuroinflammation is an underlying pathological mechanism in fibromyalgia. Toll-like receptors (TLRs) are key mediators of the immune system. TLR4 is expressed primarily in microglia and regulates downstream signaling pathways, such as MyD88/NF-κB and TRIF/IRF3. It remains unknown whether electroacupuncture (EA) has therapeutic benefit in attenuating fibromyalgia pain and what role the TLR4 pathway may play in this effect. We compared EA with sham EA to eliminate the placebo effect due to acupuncture. We demonstrated that intermittent cold stress significantly induced an increase in mechanical and thermal FM pain in mice (mechanical: 2.48 ± 0.53 g; thermal: 5.64 ± 0.32 s). EA but not sham EA has analgesic effect on FM mice. TLR4 and inflammatory mediators-related molecules were increased in FM mice’s thalamus, medial prefrontal cortex, somatosensory cortex, and amygdala, indicating neuroinflammation and microglial activation. These molecules were reduced by EA but not sham EA. Furthermore, a new chemogenetics method was used to precisely inhibit SSC activity displayed an anti-nociceptive effect through TLR4 pathway. Our results imply that the analgesic effect of EA is associated with TLR4 downregulation. We provide novel evidence that EA can modulate the TLR4 signaling pathway, revealing potential therapeutic targets for FM pain.
Keywords: 
Subject: Biology and Life Sciences  -   Neuroscience and Neurology

1. Introduction

Fibromyalgia (FM) is the third most common cause of musculoskeletal pain, next to low back pain and osteoarthritis in terms of prevalence [1]. FM prevalence is around 2.7% globally when using the 2010 modified classification criteria of the American College of Rheumatology (ACR) [2]. The total health care costs are estimated to be three times higher for fibromyalgia patients than in other populations [3]. Furthermore, fibromyalgia can cause inability to work in patients and productivity loss in society [4]. FM is a complex, chronic pain condition that primarily involves the musculoskeletal system [5]. FM pain control cannot be achieved by medication alone [6]. This unsatisfactory pharmacologic response may be related to the complex mechanism of FM pain. Current evidence suggests that immune cells, glia, and neurons form an integrated network, which coordinates immune responses and modulates the excitability of pain pathways and results in central sensitization [7]. A human functional imaging study demonstrated higher neuroinflammation in fibromyalgia patients compared to healthy controls across several brain areas (such as primary somatosensory cortex, medial frontal gyrus) [8]. In 2008, a Japanese team published a mouse model of fibromyalgia [9]. They used intermittent cold stress (ICS) to induce abnormal pain, mechanical allodynia, and hyperalgesia, which lasts for more than two weeks. This model exhibits sex hormone-independent female predominance of chronic pain, which is similar to clinical fibromyalgia patients. Our previous research showed aggravated neuroinflammation in this ICS FM model [10,11].
Toll-like receptors (TLRs) recognize microbial components and initiates signal transduction pathways to induce gene expression. These gene products control innate immune responses and the direct the development of antigen-specific acquired immunity [12]. TLRs can be classified based on their location in the cell. TLR1, TLR2, TLR4, TLR5, and TLR6 are located on the cell surface. TLR3, TLR7, TLR8, and TLR9 are found within endosomes [13]. TLR4 is a type I transmembrane protein that can be activated by bacterial lipopolysaccharide (LPS) and can induce the production of pro-inflammatory cytokines and interferons [14,15,16]. TLR4 pathway is expressed primarily in microglia and is crucial for regulating inflammatory responses [17,18,19,20]. During inflammation in the CNS, the innate immune system, including microglia, is activated by TLR4. Activated TLR4 then triggers intracellular signaling in two distinct pathways including MyD88/NF-kB and TRIF/IRF3 (15). TLR4 signaling pathway was reported to be involved in major depressive disorders (21) and peripheral neuropathic pain [22]. Loss-of-function mutation in the Tlr4 gene reliably reduced microglial activation and led to lower levels of inflammatory mediators [23]. Antagonism of TLR4 attenuated NF-κB pathways and thus decreased the inflammatory mediators secreted by astrocytes [24]. Furthermore, TLR4 plays a crucial role in central neuroinflammation including IL-1β upregulation [25].
Acupuncture is a traditional Chinese medicine technique, representing an ancient physiological system that believes health to be the result of harmony among bodily functions as well as between body and nature [26]. Acupuncture was first developed around 3000 years ago and has been shown to provide therapeutic effects in many diseases [27], especially pain management. Compared to conventional treatment, several potential advantages are associated with acupuncture, including its low cost, low number of complications [28], the ability to offer a personalized treatment program, and fewer adverse effects [29]. Clinically, two meta-analyses summarized the efficacy of acupuncture in relieving stiffness and pain in FM patients [30,31]. Regarding its mechanism, our previous article showed electroacupuncture can reduce cold stress pain through downregulation of interleukins, TNF-α, and IFN-γ in mouse plasma [11]. Electroacupuncture can reduce mechanical and thermal hyperalgesia in the same ICS mouse model of FM by attenuating brain TRPV1 signaling [10].
In this study, we aimed to determine the role of inflammatory mediators and TLR4 signaling pathway in ICS-induced FM pain mouse model. We hypothesized that inflammation underlies FM in mice through the activation of receptors such as TLR4. Furthermore, the therapeutic effect and detailed mechanisms of EA in this model are still unknown. We tested whether EA can relieve FM by regulating inflammatory mediators and related mechanisms. We also aimed to provide new evidence for the roles of EA and TLR4 in central sensitization. Chemogenetics precisely block somatosensory cortex (SSC) activity displayed an analgesic effect through TLR4 pathway. Our data provide novel evidence to support the clinical use of EA in treating FM.

2. Materials and Methods

2.1. Animals

There are totally 40 female C57BL/6 mice, aged 8-12 weeks, were used in this study. After arriving, the mice were kept in a 12h light-dark cycle with food and water ad libitum. A sample size of ten animals per group was calculated as the number required for an alpha of 0.05 and a power of 80%. In addition, the number of animals used here and their suffering were minimized. The laboratory workers were blind to treatment allocation during the experiments and analysis. The use of these animals was approved by the Institute of Animal Care and Use Committee of China Medical University (Permit no. CMUIACUC-2021-317), Taiwan, following the Guide for the use of Laboratory Animals (National Academy Press). Mice were subdivided into four groups: Normal group (Group 1: Normal); ICS-induced FM group (Group 2: FM); 2Hz Electroacupuncture group (Group 3: 2Hz EA), and Sham EA group (Group 4: Sham EA).

2.2. FM Model and Bio-Plex ELISA

All mice were host at room temperature, 24 ± 1°C, before experiments. In the ICS-induced FM model, not in normal group, 2 mice were caged in a plexiglass cage (13 × 18.8 × 29.5 cm) covered with a stainless-steel mesh. On the first day (Day 0), the mice were kept in a cold room at 4 °C overnight (from 4:00 pm-10am). The mice were next moved to 24 °C for 30 min at 10 am. After 30 min, mice were then moved back to the cold room at 4 °C for 30 min. This process was repeated for till 4:00 pm. The mice were then placed in the 4 °C cold room overnight. Normal mice were kept at room temperature from Day 0 to 7 of the experiment, with no interventions applied. Mice plasma was collected and analyzed on Bio-Plex cytokine assays (BIO-RAD, CA, USA).

2.3. EA Treatments

The mice were anaesthetized with 5% isoflurane for induction, and then maintained in 1% isoflurane. Under anesthesia, a pair of stainless-steel acupuncture needles (1” inch, 36G, YU KUANG, Taiwan) were bilaterally inserted at a depth of 3-4 mm into the murine equivalent of the human ZuSanLi (ST36) acupoints. The murine ST36 is located approximately 4 mm below and 1–2 mm lateral to the midpoint of the knee in mice. A bilateral subcutaneous transverse insertion of acupuncture needles into the scapular region was considered to be the sham acupoint 10. In both the EA and Sham EA groups, electrical stimuli were delivered by Trio 300 stimulator (Ito, Japan) at an intensity of 1mA for 20 minutes at 2Hz with a pulse width of 100 μs. The EA and Sham EA treatment caused slight visible muscle twitching around the area of insertion. The EA and Sham EA stimulation was applied thrice from day 5 to 7, following the FM protocol.

2.4. Pain Behavior Test

The mechanical and thermal pain behaviors were determined 3 times from day 5 to 7 throughout the experiment after the induction of the FM model. All mice were moved to the behavior analysis room, and were adapted to the environment for at least 30 min before behavior tests. All experiments were performed at room temperature and the stimuli were applied only when the animals were calm but not sleeping or grooming. First, the von Frey filament test was conducted. Mechanical sensitivity was measured by testing the force of responses to stimulation with 3 applications of the electronic, calibrated von Frey filament (IITC Life Science Inc., USA). Mice were placed onto a metal mesh (75x25x45 cm) and covered with a plexiglass cage (10x6x11 cm). Subjects were then mechanically stimulated by the tip of the filament at the plantar region of the right hind paw. The filament gram counts were recorded when the stimulation caused the subject to withdraw its hind paw. Second, the Hargreaves’ assessment was used to measure thermal pain behavior by testing the time of response to thermal stimulation with 3 applications using Hargreaves' test IITC analgesiometer (IITC Life Sciences, SERIES8, Model 390G). The mice were placed in a plexiglass cage on top of a glass sheet. The thermal stimulator was positioned under the glass sheet and the focus of the projection bulb was aimed exactly at the middle of the plantar surface of the right hind paw. A cut-off time of 20s was set to prevent tissue damage. In the thermal paw withdrawal test, the nociception threshold was assessed using the latency of paw withdrawal upon stimulus, and was recorded when the constant applied heat stimulation caused the subject to withdraw its hind paw.

2.5. Western Blot Analysis

The mice were anaesthetized with 1% isoflurane and cervical dislocation. The thalamus, medial frontal cortex (mPFC), SSC, and amygdala tissues were immediately excised to extract proteins. Tissues were initially placed on ice and later stored at -80°C, pending protein extraction. Total proteins were homogenized in cold radioimmunoprecipitation (RIPA) lysis buffer containing 50 mM Tris-HCl pH 7.4, 250 mM NaCl, 1% NP-40, 5 mM EDTA, 50 mM NaF, 1 mM Na3VO4, 0.02% NaN3, and 1× protease inhibitor cocktail (AMRESCO). The extracted proteins were subjected to 8% SDS-Tris glycine gel electrophoresis and transferred to a PVDF membrane. The membrane was blocked with 5% non-fat milk in TBS-T buffer (10 mM Tris pH 7.5, 100 mM NaCl, 0.1% Tween 20), incubated with a primary antibody in TBS-T with 1% bovine serum albumin (BSA) for 1 hour at room temperature antibody against TLR4 (∼35 kDa, 1 : 1000, Millipore, USA), MyD88 (∼35 kDa, 1 : 1000, Millipore, USA), TRAF6 (∼87 kDa, 1 : 1000, Millipore, USA), pERK1/2 (∼42–44 kDa, 1 : 1000, Millipore, USA), pp38 (∼41 kDa, 1 : 1000, Millipore, USA), pJNK (∼42 kDa, 1 : 1000, Millipore, USA), and pNF𝜅B (∼65 kDa, 1 : 1000, Millipore, USA), in TBS-T with 1% bovine serum albumin. Peroxidase-conjugated anti-rabbit antibody, anti-mouse antibody or anti-goat antibody (1: 5000) was used as the appropriate secondary antibody. The bands were visualized by an enhanced chemiluminescent substrate kit (PIERCE) with LAS-3000 Fujifilm (Fuji Photo Film Co., Ltd.). Where applicable, the image intensities of specific bands were quantified with NIH ImageJ software (Bethesda, MD, USA).

2.6. Immunofluorescence

Mice were euthanized with a 5% isoflurane via inhalation and intracardially perfused with normal saline followed by 4% paraformaldehyde. The brain was immediately dissected and post fixed with 4% paraformaldehyde at 4 ºC for 3 days. The tissues were placed in 30% sucrose for cryoprotection overnight at 4 ºC. The brain was embedded in an Optimal cutting temperature (OCT) compound and rapidly frozen using liquid nitrogen before storing the tissues at -80 ºC. Frozen segments were cut at 20-μm width on a cryostat then instantaneously placed on glass slides. The samples were fixed with 4% paraformaldehyde, then incubated with a blocking solution, consisting of 3% BSA, 0.1% Triton X-100, and 0.02% sodium azide, for 1 h at room temperature. After blocking, the samples were incubated with the primary antibody (1:200, Alomone), TLR4 and Iba1, prepared in 1% bovine serum albumin solution at 4ºC overnight. The samples were then incubated with the secondary antibody (1:500), 488-conjugated AffiniPure donkey anti-rabbit IgG (H + L), 594-conjugated AffiniPure donkey anti-goat IgG (H + L) and Peroxidase-conjugated AffiniPure donkey anti-mouse IgG (H + L) for 2 h at room temperature before being fixed with cover slips for immunofluorescence visualization. The samples were observed by an epi-fluorescent microscope (Olympus, BX-51, Japan) with 20 x numerical aperture (NA = 1.4) objective. The images were analyzed by NIH ImageJ software (Bethesda, MD, USA).

2.7. Chemogenetics Operation

All mice were anesthetized with isoflurane and fixed in a stereotaxic device, a 23-gauge, 2 mm stainless cannula was inserted into SSC, 0.5 mm posterior and 1.5 mm lateral of bregma at 175 μm below the cortical superficial and fixed to the skull with dental glue. The inoculation cannula was inserted and connected with the Hamilton needle through a PE duct to inject 0.3 μL of viral solution more than 3 min through the pump (KD Scientific). After injection, the cannula was well-maintained at SSC for additional 2 min to permit fluid to diffuse. Next, 0.3 μL of hM4iD DREADD (designer receptors exclusively activated by designer drugs: AAV8-hSyn-hM4D(Gi)-mCherry; Addgene Plasmid #50477) were injected into the SSC over two weeks. Clozapine N-oxide (CNO; Sigma C0832) was injected to active the DREADD. CNO was thawed in 5% dimethyl sulfoxide (DMSO; Sigma D2650) and diluted with normal saline before intraperitoneal injection of 1 mg/kg at day 0.

2.8. Statistical Analysis

Statistical analysis was performed using the SPSS statistic program. All statistic data are presented as the mean ± standard error (SEM). Statistical significance among all groups was tested using the one-way ANOVA test, followed by a post hoc Tukey’s test. Values of p < 0.05 were considered statistically significant.

3. Results

3.1. Electroacupuncture Inhibits Cold Stress-Induced FM Pain in Mice

First, we examined the effects of EA in a cold stress-induced FM pain model to address the function of EA in the subacute phase. Before FM induction, all mice had similar mechanical responses that showed normal distribution and no statistical significance between each group. Intermediate cold stress-induced typical mechanical hyperalgesia (Figure 1 A, D4: 2.48 ± 0.53, n = 10), determined by von Frey test, and was substantially attenuated by EA manipulation but not sham EA (Figure 1 A, D7: 4.89 ± 0.46 and 2.74 ± 0.51, respectively). Next, we examined whether EA or sham EA would also alter thermal hyperalgesia in FM mice. The Hargraves’ test revealed significant thermal hyperalgesia in paw withdrawal latency after cold stress induction (Figure 1 B, D4: 5.64 ± 0.32, n = 9). The latency decrease was reversed by EA but not sham EA (Figure 1 B, D7: 10.21 ± 0.74 and 6.18 ± 0.69, respectively).

3.2. Inflammatory Mediators Were Increased in FM Model and Reduced by EA but Not Sham EA

To test the role of inflammatory mediators in fibromyalgia mice, we used a Bio-Plex ELISA technique to detect these molecules in mouse plasma. Before FM induction, inflammatory mediator levels were low in normal mice. Furthermore, IL-1β, IL-2, IL-5, IL-6, IL-9, IL-12, IL-17A, TNF-α, IFN-γ, and MCP-1 were increased in FM mice (Figure 2, *P < 0.05, n = 6). EA but not sham EA dramatically reversed these phenomena (Figure 2, #P < 0.05, n = 6).

3.3. EA but Not Sham EA Reduced Fibromyalgia Pain through Microglia and TLR4 Signaling Pathways in the Mouse Thalamus

We first used Western blot to measure ionized calcium-binding adaptor molecule 1 (Iba1), a microglia/macrophage-specific calcium-binding protein [32]. We found increased amounts of Iba1 in FM mouse brain (Figure 2A, black column, *P ˂ 0.05, n = 6). EA but not sham EA significantly attenuated this phenomenon (Figure 2A, black column, #P ˂ 0.05, n = 6). Furthermore, we explored TLR4 and its downstream effectors: MyD88, TRAF6, pERK, pp38, pJNK, and the transcriptional factor pNFκB. We found increased expression of TLR4 and downstream proteins in FM mice brain (Figure 2, *P ˂ 0.05, n = 6). Similarly, EA but not sham EA significantly attenuated this phenomenon (Figure 2, *P ˂ 0.05, n = 6). We performed immunofluorescence staining of the mice thalamus. TLR4 and Iba1 were present at lower levels in normal mice but augmented in FM mice (Figure 2C, n = 4). Dealing with EA but not sham EA meaningfully abridged these phenomena (Figure 2C, n = 4). Significantly increased double-positive signals were also gotten in FM group, suggesting a colocolization of TLR4 and Iba1 (Figure 2C, n = 4). This pattern was abrogated by EA but not sham EA.

3.4. EA but Not Sham EA Reduced FM through Microglia and TLR4 Signaling Pathways in the Mouse Medial Prefrontal Cortex

We quantified the inflammatory mediators and TLR4 signaling pathway in the mouse medial prefrontal cortex. We found increased expression of Iba1 and TLR4-associated proteins in FM mouse brain: MyD88, TRAF6, pERK, pp38, pJNK, and pNFκB (Figure 3, *P ˂ 0.05, n = 6). EA but not sham EA significantly attenuated this phenomenon (Figure 3, * P ˂ 0.05, n = 6). TLR4 and Iba1 were existing at minor levels in normal mice but amplified in FM mice (Figure 3C, n = 4). Management with EA but not sham EA significantly alleviated these indications (Figure 3C, n = 4). Meaningfully greater than before double-positive signals were also achieved in FM group, telling a colocolization of TLR4 and Iba1 (Figure 3C, n = 4). The phenomenon was abolished by EA but not sham EA.

3.5. EA but Not Sham EA Reduced FM through Microglia and TLR4 Signaling Pathways in the Mouse Somatosensory Cortex

We also found increased expression of Iba1 and TLR4 downward pathway proteins in fibromyalgia mouse somatosensory cortex (SSC): MyD88, TRAF6, pERK, pp38, pJNK, and pNFκB (Figure 4, * P ˂ 0.05, n = 6). EA but not sham EA significantly attenuated this phenomenon (Figure 4, * P ˂ 0.05, n = 6). Immunofluorance signals of TLR4 and Iba1 were extant low in normal mice but augmented in FM mice (Figure 4C, n = 4). EA but not sham EA management attenuated these signals (Figure 4C, n = 4). Increased dual signals were observed in the FM group (Figure 4C, n = 4). The increased signals were alleviated by EA but not sham EA.

3.6. EA but Not Sham EA Reduced FM through Microglia and TLR4 Signaling Pathways in the Mice Amygdala

We found similar increases in expression of Iba1 and TLR4 downward pathway proteins in the amygdala of FM mice: MyD88, TRAF6, pERK, pp38, pJNK, and pNFκB (Figure 5, * P ˂ 0.05, n = 6). EA but not sham EA significantly attenuated this phenomenon (Figure 5, * P ˂ 0.05, n = 6). Immunofluorescence results in amygdala indicated that TLR4 and Iba1 were low in normal mice and increased in FM mice (Figure 5C, n = 4). EA but not sham EA relieve these increase (Figure 5C, n = 4). Increased double-positive indications were shown in FM group but abrogated by EA but not sham EA.

3.7. Chemogenetics Inhibition of SSC Attenuated Fibromyalgia Pain through TLR4 Pathway

Figure 6A shows a significant fibromyalgia-induced mechanical hyperalgesia after induction (Figure 6A, black column, n = 6). Chemogenetics inhibition of the SSC meaningfully diminished mechanical hyperalgesia (Figure 6A, red column, n = 6). In thermal hyperalgesia, mice under fibromyalgia induction presented noteworthy thermal hyperalgesia (Figure 6B, black column, n = 6). After CNO injection, thermal hyperalgesia was alleviated in mice under SSC inhibition (Figure 6B, red column, n = 6). Likewise, our results showed that TLR4 and associated molecules were reduced in SSC of fibromyalgia mice subjected to chemogenetics manipulation (Figure 6C, n = 6). Iba1 was not altered after CNO injection. TLR4, MyD88, TRAF6, pERK, pp38, pJNK, and pNFκB levels were reduced after CNO injection in the mice SSC. (Figure 6C, *P < 0.05, n = 6).

4. Discussion

We used the classical intermittent cold stress mice model to induce abnormal pain, mechanical allodynia and hyperalgesia, with gender hormone-independent female predominance of chronic pain, which was similar to clinical fibromyalgia patients. The behavior test confirmed the mechanical and thermal hyperalgesia in the FM mice. And the EA but not sham EA treatment reversed this hyperalgesia. The serum analysis showed increase pro-inflammatory cytokines: IL-1β, IL-2, IL-5, IL-6, IL-9, IL-12, IL-17A, TNF-α, IFN-γ, MCP-1 in FM mice. Similarly, EA but not sham EA treatment could reduce these pro-inflammatory cytokines. Furthermore, the molecular analysis found increased Iba1 and TLR4 related molecules in FM mice brain, indicating neuroinflammation. And EA treatment could reduce this neuroinflammation in FM mice.
Clinically, fibromyalgia is complex and difficult to treat. It was known as "nociplastic pain" [33], that chronic pain is maintained in part by central sensitization, a phenomenon of synaptic plasticity, and increased neuronal responsiveness in central pain pathways after a stimulus [34]. Fibromyalgia has features of central sensitization [35], hyperalgesia, allodynia [36], and hypersensitivity to various external stimuli such as sounds or lights [37,38]. A meta-analysis summarized 37 fibromyalgia syndrome papers, including 1264 subjects, and found consistent results of differences between FM patients and controls in pain-related brain areas [39]. So, in our study, we focused on four mice brain areas: thalamus, prefrontal cortex, somatosensory cortex, and amygdala. The thalamus and somatosensory cortex are well-known pain-related brain region. Where prefrontal cortex receives projections from the mediodorsal nucleus of the thalamus [40]. And the amygdala plays a crucial rule in emotional processing of pain. The pain modulatory role of amygdala is based on its projections to descending pain regions in the brain [41]. We found increase expression of microglia and TLR4 pathway in these four brain regions.
Central sensitization is an important phenomenon and maintains the chronicity of pain. Accumulating evidence suggests that central sensitization is driven by neuroinflammation in the peripheral and central nervous system [34]. A characteristic feature of neuroinflammation is the activation of glial cells, such as microglia and astrocytes, in the spinal cord and brain, leading to the release of proinflammatory cytokines and chemokines. Currently some researchers used proton magnetic resonance spectroscopy (MRS) to evaluate the brain neuroinflammation in patients with chronic pain [42]. Another study used 18kDa translocator protein (TSPO) as a marker of CNS inflammation, used radioligand [11C]PBR28 bonded to TPSO. Then they could visualize the CNS inflammation by PET and MRI imaging. Most importantly, they found the signal of neuroinflammation decreased after local steroid injection, consistent with clinical pain release [43]. Furthermore, a Japanese team created human-induced microglia-like (iMG) cells from human peripheral blood monocytes [44]. They found higher expression of tumor necrosis factor (TNF)-α at mRNA and protein levels in ATP-stimulated iMG cells from patients with fibromyalgia compared to cells from healthy individuals. Interestingly, there was a moderate correlation between ATP-induced upregulation of TNF-α expression and clinical parameters of subjective pain and other mental manifestations of patients with fibromyalgia. Their findings suggest that microglia in patients with fibromyalgia are hypersensitive to ATP. And TNF-α from microglia may be a key factor underlying the complex pathology of fibromyalgia. According to above research, we used the classical central sensitization mice model of fibromyalgia to investigate the neuroinflammation in mice brain, focused on the microglia.
Toll-like receptors (TLRs) are pivotal components of the innate immune response, and shaping the adaptive immune response. However, TLR signaling pathways must be tightly regulated because undue TLR stimulation may disrupt the fine balance between pro- and anti-inflammatory responses [45]. TLR4 had the role in chronic pain [46] and neuropathic pain [47], but not specific linking to fibromyalgia. Our previous research of fibromyalgia found alteration of TLR4 by electroacupuncture treatment [6].
Complementary and alternative medicine (CAM) is commonly used to treat patients with fibromyalgia [48]. About 91% of fibromyalgia patients have been identified as ever users of CAM [49]. One of the most commonly used forms of CAM is acupuncture [50]. For example, in Japan, the annual acupuncture utilization percentage ranges from 4.8 to 6.7%, with lifetime experiences of acupuncture ranging 19.4%–26.7% [51]. A recent meta-analysis summarized 15 neuroimaging studies of musculoskeletal pain, including fibromyalgia, and revealed activated clusters in multiple cortical and subcortical brain structures in response to acupuncture, including the thalamus, insula, caudate, claustrum, and lentiform nucleus [52]. However, controlling placebo effect is extremely important in acupuncture studies. In 1995, Beecher’s published an article “The powerful placebo,” which started the discussion of placebo effect [53]. In patients with fibromyalgia, an interesting review including 124 studies found that placebo treatment affected pain and other patient-centered outcomes such as sleep, fatigue, function, and overall quality of life [54]. Therefore, we compared true electroacupuncture and sham electroacupuncture even in our study. Moreover, we selected the acupoint ST36 (Zusanli) because of its well-recognized analgesic effect in mouse pain models [55].

5. Conclusions

In conclusion, we used intermittent cold stress to induce neuroinflammation and fibromyalgia-like pain in mice. In this model, electroacupuncture reduced fibromyalgia pain via neuronal/microglial inactivation and decreased TLR4 signaling in the mouse brain. This study opens a window of opportunity for the discovery and development of novel TLR4-targeted pain treatments in patients with fibromyalgia. In this study, we analyzed the specific TLR4 pathway in the mice model of fibromyalgia. Our data provide clinical implications of TLR4 in acupuncture analgesia (Figure 7).

Author Contributions

Conceptualization, methodology, software, validation, and investigation, ST Tsai and CC Yang; resources, data curation, writing—original draft preparation, writing—review and editing, and funding acquisition, HY Liao and YW Lin. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the following grants: NSTC 112-2314-B-039-021, CMU112-MF-99, and the "Chinese Medicine Research Center, China Medical University" from The Featured Areas Research Center Program within the framework of the Higher Education Sprout Project by the Ministry of Education (MOE) in Taiwan.

Institutional Review Board Statement

The use of these animals was approved by the Institute of Animal Care and Use Committee of China Medical University (Permit no. CMUIACUC-2021-317), Taiwan, following the Guide for the use of Laboratory Animals (National Academy Press).

Data Availability Statement

The datasets supporting the conclusions of this article are included within the article.

Acknowledgments

The authors would like to thank Enago (www.enago.tw) for the English language review.

Conflicts of Interest

There are no financial or other relationships that might lead to a conflict of interest for all authors.

References

  1. Spaeth, M. Epidemiology, costs, and the economic burden of fibromyalgia. Arthritis Res. Ther. 2009, 11, 117. [Google Scholar] [CrossRef] [PubMed]
  2. Queiroz, L.P. Worldwide epidemiology of fibromyalgia. Curr. Pain. Headache Rep. 2013, 17, 356. [Google Scholar] [CrossRef] [PubMed]
  3. Knight, T.; et al. Health-resource use and costs associated with fibromyalgia in France, Germany, and the United States. Clin. Outcomes Res. 2013, 5, 171–180. [Google Scholar] [CrossRef] [PubMed]
  4. Guymer, E.K.; Littlejohn, G.O.; Brand, C.K.; Kwiatek, R.A. Fibromyalgia onset has a high impact on work ability in Australians. Intern. Med. J. 2016, 46, 1069–1074. [Google Scholar] [CrossRef] [PubMed]
  5. Sarzi-Puttini, P.; Giorgi, V.; Marotto, D.; Atzeni, F. Fibromyalgia: An update on clinical characteristics, aetiopathogenesis and treatment. Nat. Rev. Rheumatol. 2020, 16, 645–660. [Google Scholar] [CrossRef] [PubMed]
  6. Liao, H.-Y.; Lin, Y.-W. Electroacupuncture reduces cold stress-induced pain through microglial inactivation and transient receptor potential V1 in mice. Chin. Med. 2021, 16, 43. [Google Scholar] [CrossRef] [PubMed]
  7. Ren, K.; Dubner, R. Interactions between the immune and nervous systems in pain. Nat. Med. 2010, 16, 1267–1276. [Google Scholar] [CrossRef] [PubMed]
  8. Seo, S.; et al. Abnormal neuroinflammation in fibromyalgia and CRPS using [11C]-(R)-PK11195 PET. PLoS ONE 2021, 16, e0246152. [Google Scholar] [CrossRef] [PubMed]
  9. Nishiyori, M.; Ueda, H. Prolonged gabapentin analgesia in an experimental mouse model of fibromyalgia. Mol. Pain. 2008, 4, 52. [Google Scholar] [CrossRef]
  10. Lottering, B.; Lin, Y.-W. Functional characterization of nociceptive mechanisms involved in fibromyalgia and electroacupuncture. Brain Res. 2021, 1755, 147260. [Google Scholar] [CrossRef]
  11. Lin, Y.-W.; Chou, A.I.W.; Su, H.; Su, K.-P. Transient receptor potential V1 (TRPV1) modulates the therapeutic effects for comorbidity of pain and depression: The common molecular implication for electroacupuncture and omega-3 polyunsaturated fatty acids. Brain Behav. Immun. 2020, 89, 604–614. [Google Scholar] [CrossRef] [PubMed]
  12. Takeda, K.; Akira, S. Toll-like receptors in innate immunity. Int. Immunol. 2005, 17, 1–14. [Google Scholar] [CrossRef] [PubMed]
  13. Kumar, H.; Kawai, T.; Akira, S. Toll-like receptors and innate immunity. Biochem. Biophys. Res. Commun. 2009, 388, 621–625. [Google Scholar] [CrossRef]
  14. Kawai, T.; Akira, S. The role of pattern-recognition receptors in innate immunity: Update on Toll-like receptors. Nat. Immunol. 2010, 11, 373–384. [Google Scholar] [CrossRef]
  15. Kuzmich, N.N.; et al. TLR4 Signaling Pathway Modulators as Potential Therapeutics in Inflammation and Sepsis. Vaccines 2017, 5. [Google Scholar] [CrossRef] [PubMed]
  16. Hsu, H.C.; Hsieh, C.L.; Wu, S.Y.; Lin, Y.W. Toll-like receptor 2 plays an essential role in electroacupuncture analgesia in a mouse model of inflammatory pain. Acupunct. Med. 2019, 37, 356–364. [Google Scholar] [CrossRef] [PubMed]
  17. Capiralla, H.; et al. Resveratrol mitigates lipopolysaccharide- and Aβ-mediated microglial inflammation by inhibiting the TLR4/NF-κB/STAT signaling cascade. J. Neurochem. 2012, 120, 461–472. [Google Scholar] [CrossRef] [PubMed]
  18. Facci, L.; et al. Toll-Like Receptors 2, -3 and -4 Prime Microglia but not Astrocytes Across Central Nervous System Regions for ATP-Dependent Interleukin-1β Release. Sci. Rep. 2014, 4, 6824. [Google Scholar] [CrossRef]
  19. Ghowsi, M.; Qalekhani, F.; Farzaei, M.H.; Mahmudii, F.; Yousofvand, N.; Joshi, T. Inflammation, oxidative stress, insulin resistance, and hypertension as mediators for adverse effects of obesity on the brain: A review. Biomedicine 2021, 11, 13–22. [Google Scholar] [CrossRef]
  20. Liao, E.-T.; Lin, Y.-W.; Huang, C.-P.; Tang, N.-Y.; Hsieh, C.-L. Electric Stimulation of Ear Reduces the Effect of Toll-Like Receptor 4 Signaling Pathway on Kainic Acid-Induced Epileptic Seizures in Rats. BioMed Res. Int. 2018, 2018, 5407256. [Google Scholar] [CrossRef]
  21. Wu, M.K.; Huang, T.L.; Huang, K.W.; Huang, Y.L.; Hung, Y.Y. Association between toll-like receptor 4 expression and symptoms of major depressive disorder. Neuropsychiatr. Dis. Treat. 2015, 11, 1853–1857. [Google Scholar]
  22. Thakur, K.K.; et al. Therapeutic implications of toll-like receptors in peripheral neuropathic pain. Pharmacol. Res. 2017, 115, 224–232. [Google Scholar] [CrossRef]
  23. Walter, S.; et al. Role of the toll-like receptor 4 in neuroinflammation in Alzheimer's disease. Cell Physiol. Biochem. 2007, 20, 947–956. [Google Scholar] [CrossRef] [PubMed]
  24. Alfonso-Loeches, S.; Pascual-Lucas, M.; Blanco, A.M.; Sanchez-Vera, I.; Guerri, C. Pivotal role of TLR4 receptors in alcohol-induced neuroinflammation and brain damage. J. Neurosci. 2010, 30, 8285–8295. [Google Scholar] [CrossRef]
  25. Chakravarty, S.; Herkenham, M. Toll-Like Receptor 4 on Nonhematopoietic Cells Sustains CNS Inflammation during Endotoxemia, Independent of Systemic Cytokines. J. Neurosci. 2005, 25, 1788–1796. [Google Scholar] [CrossRef]
  26. Su, K.-P.; Chou, L.-W.; Sun, M.-F.; Lin, J.-G. Acupuncture Treatment in Depression. in Experimental Acupuncturology (ed. Lin, J.-G.) 43-66 (Springer Singapore, Singapore, 2018).
  27. Tsai, S.T.; et al. Transient receptor potential V1 modulates neuroinflammation in Parkinson's disease dementia: Molecular implications for electroacupuncture and rivastigmine. Iran. J. Basic. Med. Sci. 2021, 24, 1336–1345. [Google Scholar] [PubMed]
  28. Lin, S.-K.; et al. Incidence of iatrogenic pneumothorax following acupuncture treatments in Taiwan. Acupunct. Med. 2019, 37, 332–339. [Google Scholar] [CrossRef] [PubMed]
  29. NIH Consensus Conference. Acupuncture. Jama 1998, 280, 1518–1524.
  30. Deare, J.C.; et al. Acupuncture for treating fibromyalgia. Cochrane Database Syst. Rev. 2013, 2013, Cd007070. [Google Scholar] [CrossRef]
  31. Yang, B.; et al. Efficacy of acupuncture on fibromyalgia syndrome: A meta-analysis. J. Tradit. Chin. Med. 2014, 34, 381–391. [Google Scholar]
  32. Ohsawa, K.; Imai, Y.; Sasaki, Y.; Kohsaka, S. Microglia/macrophage-specific protein Iba1 binds to fimbrin and enhances its actin-bundling activity. J. Neurochem. 2004, 88, 844–856. [Google Scholar] [CrossRef] [PubMed]
  33. Kosek, E.; et al. Do we need a third mechanistic descriptor for chronic pain states? Pain 2016, 157, 1382–1386. [Google Scholar] [CrossRef] [PubMed]
  34. Ji, R.-R.; Nackley, A.; Huh, Y.; Terrando, N.; Maixner, W. Neuroinflammation and Central Sensitization in Chronic and Widespread Pain. Anesthesiology 2018, 129, 343–366. [Google Scholar] [CrossRef] [PubMed]
  35. Borchers, A.T.; Gershwin, M.E. Fibromyalgia: A Critical and Comprehensive Review. Clin. Rev. Allergy Immunol. 2015, 49, 100–151. [Google Scholar] [CrossRef] [PubMed]
  36. Desmeules, J.A.; et al. Neurophysiologic evidence for a central sensitization in patients with fibromyalgia. Arthritis Rheum. 2003, 48, 1420–1429. [Google Scholar] [CrossRef] [PubMed]
  37. Geisser, M.E.; et al. A psychophysical study of auditory and pressure sensitivity in patients with fibromyalgia and healthy controls. J. Pain. 2008, 9, 417–422. [Google Scholar] [CrossRef] [PubMed]
  38. Martenson, M.E.; et al. A possible neural mechanism for photosensitivity in chronic pain. Pain. 2016, 157, 868–878. [Google Scholar] [CrossRef] [PubMed]
  39. Dehghan, M.; et al. Coordinate-based (ALE) meta-analysis of brain activation in patients with fibromyalgia. Hum. Brain Mapp. 2016, 37, 1749–1758. [Google Scholar] [CrossRef]
  40. Murillo-Garcia, A.; Leon-Llamas, J.L.; Villafaina, S.; Gusi, N. Fibromyalgia impact in the prefrontal cortex subfields: An assessment with MRI. Clin. Neurol. Neurosurg. 2022, 219, 107344. [Google Scholar] [CrossRef]
  41. Simons, L.E.; et al. The human amygdala and pain: Evidence from neuroimaging. Hum. Brain Mapp. 2014, 35, 527–538. [Google Scholar] [CrossRef]
  42. Jung, C.; et al. Magnetic resonance imaging of neuroinflammation in chronic pain: A role for astrogliosis? Pain 2020, 161, 1555–1564. [Google Scholar] [CrossRef]
  43. Albrecht, D.S.; et al. Neuroinflammation of the spinal cord and nerve roots in chronic radicular pain patients. Pain. 2018, 159, 968–977. [Google Scholar] [CrossRef]
  44. Ohgidani, M.; et al. Fibromyalgia and microglial TNF-α: Translational research using human blood induced microglia-like cells. Sci. Rep. 2017, 7, 11882. [Google Scholar] [CrossRef] [PubMed]
  45. Anwar, M.A.; Basith, S.; Choi, S. Negative regulatory approaches to the attenuation of Toll-like receptor signaling. Exp. Mol. Med. 2013, 45, e11. [Google Scholar] [CrossRef] [PubMed]
  46. Lacagnina, M.J.; Watkins, L.R.; Grace, P.M. Toll-like receptors and their role in persistent pain. Pharmacol. Ther. 2018, 184, 145–158. [Google Scholar] [CrossRef] [PubMed]
  47. Liu, X.; et al. Toll-like receptors and their role in neuropathic pain and migraine. Mol. Brain 2022, 15, 73. [Google Scholar] [CrossRef] [PubMed]
  48. Branco, J.; et al. Fibromyalgia Syndrom: A european epidemiological survey. in Annals of the Rheumatic Diseases, Vol. 64 547-+ (BMJ PUBLISHING GROUP BRITISH MED ASSOC HOUSE, TAVISTOCK SQUARE, LONDON WC1H …, 2005).
  49. Pioro-Boisset, M.; Esdaile, J.M.; Fitzcharles, M.A. Alternative medicine use in fibromyalgia syndrome. Arthritis Care Res. 1996, 9, 13–17. [Google Scholar] [CrossRef] [PubMed]
  50. Itoh, K.; Kitakoji, H. Effects of acupuncture to treat fibromyalgia: A preliminary randomised controlled trial. Chin. Med. 2010, 5, 11. [Google Scholar] [CrossRef] [PubMed]
  51. Ishizaki, N.; Yano, T.; Kawakita, K. Public status and prevalence of acupuncture in Japan. Evid. Based Complement. Altern. Med. 2010, 7, 493–500. [Google Scholar] [CrossRef]
  52. Ha, G.; et al. Coordinate-based (ALE) meta-analysis of acupuncture for musculoskeletal pain. Front. Neurosci. 2022, 16. [Google Scholar] [CrossRef]
  53. Beecher, H.K. THE POWERFUL PLACEBO. J. Am. Med. Assoc. 1955, 159, 1602–1606. [Google Scholar] [CrossRef] [PubMed]
  54. Chen, X.; et al. The placebo effect and its determinants in fibromyalgia: Meta-analysis of randomised controlled trials. Clin. Rheumatol. 2017, 36, 1623–1630. [Google Scholar] [CrossRef] [PubMed]
  55. Wu, S.-Y.; Chen, W.-H.; Hsieh, C.-L.; Lin, Y.-W. Abundant expression and functional participation of TRPV1 at Zusanli acupoint (ST36) in mice: Mechanosensitive TRPV1 as an “acupuncture-responding channel”. BMC Complement. Altern. Med. 2014, 14, 96. [Google Scholar] [CrossRef] [PubMed]
Preprints 97182 g001
Figure 1. Mechanical withdrawal, thermal latency, and the concentration of inflammatory mediators in normal, FM, EA, and sham EA mice. (A) Mechanical threshold from the von Frey tests. (B) Thermal latency from the Hargreaves’ test. (C) IL-1β, IL-2, IL-5, IL-6, and IL-9 (D) IL-12, IL-17A, TNF-α, IFN-γ, and MCP-1 in the mice plasma. *means statistical difference when compared to the normal mice. #indicates statistical significance when compared to the FM groups. IL= Interleukin; IFN= Interferon; TNF= Tumor necrosis factor; MCP= Monocyte Chemoattractant Protein. n=6 in all groups.
Figure 1. Mechanical withdrawal, thermal latency, and the concentration of inflammatory mediators in normal, FM, EA, and sham EA mice. (A) Mechanical threshold from the von Frey tests. (B) Thermal latency from the Hargreaves’ test. (C) IL-1β, IL-2, IL-5, IL-6, and IL-9 (D) IL-12, IL-17A, TNF-α, IFN-γ, and MCP-1 in the mice plasma. *means statistical difference when compared to the normal mice. #indicates statistical significance when compared to the FM groups. IL= Interleukin; IFN= Interferon; TNF= Tumor necrosis factor; MCP= Monocyte Chemoattractant Protein. n=6 in all groups.
Preprints 97182 g001
Preprints 97182 g002
Figure 2. The levels of Iba1 and TLR4 related molecules in the mice thalamus. The western blot bands contain four lanes of protein expression in order of Normal, FM, 2Hz EA, and sham EA. (A) Iba1, TLR4, MyD88, TRAF6. (B) pERK, pp38, pJNK, and pNFκB. *indicates statistical significance when compared with the normal group. #indicates statistical significance when compared to the FM group. n = 6 in all groups. (C) Immunofluorescence staining of TLR4, Iba1, and double staining protein expression in the mice thalamus. Scale bar means 100 μm. n = 4 in all groups.
Figure 2. The levels of Iba1 and TLR4 related molecules in the mice thalamus. The western blot bands contain four lanes of protein expression in order of Normal, FM, 2Hz EA, and sham EA. (A) Iba1, TLR4, MyD88, TRAF6. (B) pERK, pp38, pJNK, and pNFκB. *indicates statistical significance when compared with the normal group. #indicates statistical significance when compared to the FM group. n = 6 in all groups. (C) Immunofluorescence staining of TLR4, Iba1, and double staining protein expression in the mice thalamus. Scale bar means 100 μm. n = 4 in all groups.
Preprints 97182 g002
Preprints 97182 g003
Figure 3. The levels of Iba1 and TLR4 related molecules in the mice mPFC. The western blot bands contain four lanes of protein expression in order of Normal, FM, 2Hz EA, and sham EA. (A) Iba1, TLR4, MyD88, TRAF6. (B) pERK, pp38, pJNK, and pNFκB. *indicates statistical significance when compared with the normal group. #indicates statistical significance when compared to the FM group. n = 6 in all groups. (C) Immunofluorescence staining of TLR4, Iba1, and double staining protein expression in the mice mPFC. Scale bar means 100 μm. n = 4 in all groups.
Figure 3. The levels of Iba1 and TLR4 related molecules in the mice mPFC. The western blot bands contain four lanes of protein expression in order of Normal, FM, 2Hz EA, and sham EA. (A) Iba1, TLR4, MyD88, TRAF6. (B) pERK, pp38, pJNK, and pNFκB. *indicates statistical significance when compared with the normal group. #indicates statistical significance when compared to the FM group. n = 6 in all groups. (C) Immunofluorescence staining of TLR4, Iba1, and double staining protein expression in the mice mPFC. Scale bar means 100 μm. n = 4 in all groups.
Preprints 97182 g003
Preprints 97182 g004
Figure 4. The levels of Iba1 and TLR4 related molecules in the mice SSC. The western blot bands contain four lanes of protein expression in order of Normal, FM, 2Hz EA, and sham EA. (A) Iba1, TLR4, MyD88, TRAF6. (B) pERK, pp38, pJNK, and pNFκB. *indicates statistical significance when compared with the normal group. #indicates statistical significance when compared to the FM group. n = 6 in all groups. (C) Immunofluorescence staining of TLR4, Iba1, and double staining protein expression in the mice SSC. Scale bar means 100 μm. n = 4 in all groups.
Figure 4. The levels of Iba1 and TLR4 related molecules in the mice SSC. The western blot bands contain four lanes of protein expression in order of Normal, FM, 2Hz EA, and sham EA. (A) Iba1, TLR4, MyD88, TRAF6. (B) pERK, pp38, pJNK, and pNFκB. *indicates statistical significance when compared with the normal group. #indicates statistical significance when compared to the FM group. n = 6 in all groups. (C) Immunofluorescence staining of TLR4, Iba1, and double staining protein expression in the mice SSC. Scale bar means 100 μm. n = 4 in all groups.
Preprints 97182 g004
Preprints 97182 g005
Figure 5. The levels of Iba1 and TLR4 related molecules in the mice amygdala. The western blot bands contain four lanes of protein expression in order of Normal, FM, 2Hz EA, and sham EA. (A) Iba1, TLR4, MyD88, TRAF6. (B) pERK, pp38, pJNK, and pNFκB. *indicates statistical significance when compared with the normal group. #indicates statistical significance when compared to the FM group. n = 6 in all groups. (C) Immunofluorescence staining of TLR4, Iba1, and double staining protein expression in the mice amygdala. Scale bar means 100 μm. n = 4 in all groups.
Figure 5. The levels of Iba1 and TLR4 related molecules in the mice amygdala. The western blot bands contain four lanes of protein expression in order of Normal, FM, 2Hz EA, and sham EA. (A) Iba1, TLR4, MyD88, TRAF6. (B) pERK, pp38, pJNK, and pNFκB. *indicates statistical significance when compared with the normal group. #indicates statistical significance when compared to the FM group. n = 6 in all groups. (C) Immunofluorescence staining of TLR4, Iba1, and double staining protein expression in the mice amygdala. Scale bar means 100 μm. n = 4 in all groups.
Preprints 97182 g005
Preprints 97182 g006
Figure 6. Pain behavior of FM and FM mice treated with the chemogenetics technique. Black column: FM group, red column: FM treated with chemogenetics. * P < 0.05 compared with the FM group. (A) Mechanical hyperalgesia (von Frey test). (B) Thermal hyperalgesia (Hargreaves’ test). (C) Protein levels of Iba1, TLR4, MyD88, TRAF6, pERK, pJNK, pp38, and pNFkB were measured in mice SSC.
Figure 6. Pain behavior of FM and FM mice treated with the chemogenetics technique. Black column: FM group, red column: FM treated with chemogenetics. * P < 0.05 compared with the FM group. (A) Mechanical hyperalgesia (von Frey test). (B) Thermal hyperalgesia (Hargreaves’ test). (C) Protein levels of Iba1, TLR4, MyD88, TRAF6, pERK, pJNK, pp38, and pNFkB were measured in mice SSC.
Preprints 97182 g006
Preprints 97182 g007
Figure 7. Toll-like receptor (TLR) signaling pathway, focusing on TLR4. The red boxes are the proteins we analyzed in this study. Abbreviations: ERK= extracellular signal-regulated kinase; IFN= interferon; JNK= c-Jun N-terminal kinase; MyD88= myeloid differentiation primary-response protein 88; TLR, Toll-like receptors; TRAF= tumor necrosis factor receptor-associated factor.
Figure 7. Toll-like receptor (TLR) signaling pathway, focusing on TLR4. The red boxes are the proteins we analyzed in this study. Abbreviations: ERK= extracellular signal-regulated kinase; IFN= interferon; JNK= c-Jun N-terminal kinase; MyD88= myeloid differentiation primary-response protein 88; TLR, Toll-like receptors; TRAF= tumor necrosis factor receptor-associated factor.
Preprints 97182 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated