# TRPM8: A therapeutic target for neuroinflammatory symptoms induced by severe dry eye disease

Darine Fakih, Christophe Baudouin, Annabelle Réaux Le Goazigo, Stéphane Mélik Parsadaniantz

► To cite this version:

Darine Fakih, Christophe Baudouin, Annabelle Réaux Le Goazigo, Stéphane Mélik Parsadaniantz.

TRPM8: A therapeutic target for neuroinflammatory symptoms induced by severe dry eye disease.

International Journal of Molecular Sciences, 2020, 21 (22), pp.1-21. 10.3390/ijms21228756. hal-03032264

HAL Id: hal-03032264

https://hal.science/hal-03032264

Submitted on 7 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

Article

# TRPM8: A Therapeutic Target for Neuroinflammatory Symptoms Induced by Severe Dry Eye Disease

Darine Fakih <fnref n="1" />,<fnref n="2" />, Christophe Baudouin <fnref n="1" />,<fnref n="3" />,<fnref n="4" /> ID, Annabelle Réaux-Le Goazigo <fnref n="1" />,† ID
and Stéphane Mélik Parsadaniantz <fnref n="1" />,*,† ID

Received: 14 October 2020; Accepted: 17 November 2020; Published: 19 November 2020

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** Dry eye disease (DED) is commonly associated with ocular surface inflammation and pain. In this study, we evaluated the effectiveness of repeated instillations of transient receptor potential melastatin 8 (TRPM8) ion channel antagonist M8-B on a mouse model of severe DED induced by the excision of extra-orbital lacrimal and Harderian glands. M8-B was topically administered twice a day from day 7 until day 21 after surgery. Cold and mechanical corneal sensitivities and spontaneous ocular pain were monitored at day 21. Ongoing and cold-evoked ciliary nerve activities were next evaluated by electrophysiological multi-unit extracellular recording. Corneal inflammation and expression of genes related to neuropathic pain and inflammation were assessed in the trigeminal ganglion. We found that DED mice developed a cold allodynia consistent with higher TRPM8 mRNA expression in the trigeminal ganglion (TG). Chronic M8-B instillations markedly reversed both the corneal mechanical allodynia and spontaneous ocular pain commonly associated with persistent DED. M8-B instillations also diminished the sustained spontaneous and cold-evoked ciliary nerve activities observed in DED mice as well as inflammation in the cornea and TG. Overall, our study provides new insight into the effectiveness of TRPM8 blockade for alleviating corneal pain syndrome associated with severe DED, opening a new avenue for ocular pain management.

**Keywords:** cornea; electrophysiology; inflammation; dry eye; TRPM8 antagonist

## 1. Introduction

Dry eye disease (DED) is a chronic ocular disease characterized by dryness, discomfort, burning sensation, and pain in the eye [1]. The prevalence of DED can reach up to 50% in people aged 65 years and over [2]. DED occurs due to deficient tear production or excessive evaporation, causing damage to the ocular surface. Persistent ocular surface pain, which can occur in moderate to severe DED, has gained recognition due to its increasing prevalence, morbidity, and the resulting social burden [3]. However, the underlying mechanisms are not fully understood and, to date, no therapeutic strategy is satisfactory for alleviating ocular pain associated with chronic DED.

The cornea is a transparent avascular tissue considered to have the densest nociceptive nerve terminals of all body tissues [3–5]. The nociceptive innervation of the cornea is provided by the ciliary nerve, where neurons are located in the ophthalmic branch of the trigeminal ganglion (TG) [3,6,7]. Different types of receptors (mechano-nociceptors, polymodal nociceptors, and thermoreceptors)

<footnote n="1">¹ Sorbonne Université, INSERM, CNRS, Institut de la Vision, 17 rue Moreau, F-75012 Paris, France; darine.fakih@inserm.fr (D.F.); cbaudouin@15-20.fr (C.B.); annabelle.reaux@inserm.fr (A.R.-L.G.)</footnote>

<footnote n="2">² R&D Department, Laboratoires Théa, 12 rue Louis Biérot, F-63000 Clermont-Ferrand, France</footnote>

<footnote n="3">³ CHNO des Quinze-Vings, INSERM-DGOS CIC 1423, 17 rue Moreau, F-75012 Paris, France</footnote>

<footnote n="4">⁴ Department of Ophthalmology, Ambroise Paré Hospital, AP-HP, University of Versailles Saint-Quentin-en-Yvelines, 9 avenue Charles de Gaulle, F-92100 Boulogne-Billancourt, France</footnote>

<footnote>* Correspondence: stephane.melik-parsadaniantz@inserm.fr; Tel.: +33-153-462-572</footnote>

<footnote>† These authors contributed equally to this work.</footnote>
coexist on a single corneal nerve ending [8]. Corneal cold thermoreceptors, which represent 50% of the total population of corneal sensory neurons in mice, detect ocular dryness and maintain tear homeostasis by adjusting the basal tearing rate and spontaneous blinking frequency [9–12]. The transient receptor potential melastatin 8 (TRPM8) receptor is a cold- and menthol-sensitive non-selective cation channel [13] that is expressed by a subpopulation of cold sensory trigeminal neurons [14–16]. Genetic deletion of TRPM8 revealed that TRPM8-/- mice are deficient in unpleasant cold sensitivity, confirming the implication of this channel in cold-pain sensation [17–19]. Cold allodynia is a common feature of neuropathic pain [20,21] and numerous preclinical studies demonstrated the essential function of TRPM8 in neuropathic pain perception [22–26]. In this regard, pharmacological studies have demonstrated that TRPM8 ligands or inhibitors/blockers attenuate pain sensation in numerous somatic pain models. TRPM8 antagonists were shown to efficiently alleviate acute and chronic pain [27–30], whereas the TRPM8 agonist may present significant antiallodynic activity through an excessive activation of TRPM8, leading to its downregulation [31].

In DED, noxious (thermal, mechanical, and chemical) stimuli and inflammation enhance the activation of corneal sensory neurons terminals, leading to the intensification of pain sensations [8]. Electrophysiological recording studies showed that DED increases ongoing corneal activities, cold nociceptors responsiveness, and TG cold neurons sensitivity in models of DED in guinea pigs and rats [32,33]. Acute ocular instillation of TRPM8 antagonist was reported to decrease the cold sensitivity of corneal neurons to dry stimulations in wild-type rat [34]. TRPM8 plays an essential role in the mechanisms of tearing and blinking [10,35,36]. Some studies supported the hypothesis that TRPM8 agonist may alleviate chronic pain and dryness in DED [37], whereas others reported that TRPM8 over-activation could result in the desensitization of trigeminal neurons, leading to reductions in tear production, which could have a more deleterious effect on DED patients [33]. Preclinical and clinical studies have reported that ocular menthol application at low doses produces a cooling sensation, whereas at higher doses, it evokes burning, irritation, and pain [35,38,39]. Topical application of menthol (0.2 mM) to human blunted cold-induced perceptions while the irritation sensation caused by low temperature stimuli endured [40].

The aim of this study was to provide new insight into the effectiveness of TRPM8 antagonist topical treatment on corneal neurosensory abnormality syndrome induced by severe DED. Herein, we report higher TRPM8 mRNA expression in the ophthalmic branch of TG and corneal cold hypersensitivity in a mouse model of severe DED [41]. We found in an ex vivo eye preparation of DED mice that single application of TRPM8 antagonist (M8-B) at 20 and 30 µM decreased cold-evoked ciliary nerve activity compared to a saline-treated DED group. We further found in vivo that chronic topical instillations of (20 µM) M8-B twice daily for 15 days in DED mice decreased (1) the spontaneous and cold-evoked activity of corneal nerve in ex vivo eyes preparations, (2) corneal inflammation, (3) expression of genes related to neuropathic and inflammatory pain in the TG, and (4) corneal mechanical allodynia and spontaneous ocular pain. The research design of this study is presented in Figure 1.
<figure><img src="image_4.png" /></figure>

**Figure 1.** Research design. (A) Evaluation of: corneal cold sensitivity using an ocular instillation of menthol in vivo, TRPM8 mRNA expression using RNAscope in situ hybridization in the ophthalmic branch of TG and the effectiveness of acute instillation of TRPM8 antagonist M8-B at 10, 20, and 30 µM using multi-unit extracellular recording of spontaneous (32 °C) and evoked (20 °C) ciliary nerve fiber activities in DED and sham mice 21 days post-surgery; (B) The effectiveness of the chronic ocular M8-B instillations (twice a day from day 7 until day 21 after surgery) in the mice on: spontaneous (32 °C) and evoked (20 °C) ciliary nerve fiber activities in ex-vivo eye preparations, mechanical corneal sensitivity using von Frey test in vivo, spontaneous eye pain measuring the eye closure ratio in vivo, corneal inflammation using in vivo confocal microscopy, and mRNA expression of genes implicated in pain conduction and inflammation using RT-qPCR in DED mice instilled with M8-B (20 µM) compared to DED mice instilled with PBS.

## 2. Results

### 2.1. Severe DED Induced Corneal Cold Hypersensitivity Correlated with an Increase in TRPM8 mRNA Expression in the TG

We first investigated whether a severe DED induced cold hypersensitivity at day 21. To do so, a drop of 50 µM menthol, a TRPM8 agonist, was directly instilled to the cornea and the eye closing ratio was then monitored. Figure 2A shows the eye closing ratio decreased by 23% in DED mice compared to sham mice (0.52 ± 0.03 vs. 0.39 ± 0.03, *p* < 0.05), suggesting that DED animals develop cold hypersensitivity. We next investigated whether this sensory abnormality could be due to changes in TRPM8 mRNA expression in the TG from DED mice using RNA scope in situ hybridization. Semi-quantitative analysis of the area covered by the TRPM8 mRNA fluorescent probe showed that TRPM8 mRNA expression was significantly higher (+254.84%) in the ophthalmic branch of the TG of DED mice compared to sham animals (sham: 0.31 ± 0.03 vs. DED: 1.10 ± 0.21, *p* < 0.05; Figure 2B).
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Evaluation of corneal cold sensitivity and TRPM8 expression in the ophthalmic branch of the TGs of sham and DED mice at day 21. (A) The cold corneal sensitivity was evaluated using a drop of menthol and recording the eye closing ratio after 5 min (min). (B) Staining (white arrows) and quantification of mRNA levels using in situ hybridization of TRPM8 (Opal 650) in the ophthalmic branch of TGs of sham and DED mice. All experiments were conducted at day 21. Sham and DED mice *n* = 4–5. Scale bar = 50 µM. \* *p* < 0.05 relative to sham group. Results are expressed as mean ± standard error of the mean (SEM). For statistical analysis, the Kolmogorov–Smirnov test was performed followed by a nonparametric Mann–Whitney using GraphPad Prism version 7.00 (GraphPad Software, La Jolla, CA, USA).</figcaption></figure>

## 2.2. Instillation of M8-B Reduced Ongoing and Cold-Induced Electrical Activities of Thermoreceptors in Ex Vivo Eye Preparations from DED Animals

Ongoing and evoked electrical activities of the ciliary nerve (Figure 3A) were recorded at day 21 post-surgery using a multi-unit extracellular recording in ex vivo eye preparations from sham and DED animals. First, we analyzed the spontaneous nerve activities while perfusing the cornea with superfusion saline solution at a temperature set to 32–34 °C. Figure 3B shows histograms
representing the number of impulses (imp) per second (sec) from sham and DED mice at 32–34 °C. DED mice exhibited a higher spontaneous activity (+72.29%) compared to sham animals (sham + vehicle: 55.25 ± 4.34 imp/sec, DED + vehicle: 95.19 ± 7.38 imp/sec, *p* < 0.00001; Figure 3B).

<figure><img src="image_6.png" /><figcaption>**Figure 3.** The acute instillation of TRPM8 antagonist M8-B in sham and DED ex vivo eye preparations. (A) Electrophysiological trace illustrating the extracellular spontaneous (32 °C) and evoked (20 °C) activity of the ciliary nerve fiber in control mice; time scale = 2 s. Histograms illustrating (B) the spontaneous activity of the ciliary nerve fiber after ex vivo instillation of 10, 20, and 30 µM M8-B in sham and DED at 32 °C; (C) the evoked activity of the ciliary nerve fiber after ex vivo instillation of 10, 20, and 30 µM M8-B in sham and DED at 20 °C at day 21. Sham and DED mice *n* = 5–8. \* *p* < 0.05 and \*\* *p* < 0.01 relative to DED vehicle group at 20 °C. \*\*\* *p* < 0.001 relative to sham vehicle group at 32 °C and \*\*\*\* *p* < 0.01 relative to sham vehicle group at 20 °C. Results are expressed as the mean ± SEM. For statistical analysis, the Kolmogorov–Smirnov test was performed followed by a parametric one-way ANOVA test using GraphPad Prism version 7.00 (GraphPad Software, La Jolla, CA, USA).</figcaption></figure>

Next, we recorded the spontaneous activity of the ciliary nerve while perfusing the cornea with saline solution at 32–34 °C containing 10, 20, or 30 µM of M8-B. The spontaneous activity of
the ciliary nerve remained unchanged for all M8-B concentrations tested at 32 °C (sham + vehicle: 55.25 ± 4.34 imp/sec vs. sham + M8-B (10 µM) 68.60 ± 10.09 imp/sec vs. sham + M8-B (20 µM) 68.40 ± 17.98 imp/sec vs. sham + M8-B (30 µM) 72.00 ± 14.45 imp/sec, *p* > 0.05; Figure 3A) and DED + vehicle (95.19 ± 7.38 imp/sec) vs. DED + M8-B (10 µM; 88.80 ± 12.18 imp/sec) vs. DED + M8-B (20 µM; 74.80 ± 19.00 imp/sec) vs. DED + M8-B (30 µM; 76.20 ± 20.27 imp/sec; *p* > 0.05; Figure 3B).

Then, we recorded the cold-evoked activities in response to a cold stimulation of the ciliary nerve while perfusing the cornea with superfusion saline solution with temperature set to 20 °C. Figure 3C shows histograms representing the number of imp per sec from sham and DED mice. At 20 °C, we observed a higher nerve electrical activity (+36.88%) in DED ex vivo eye preparations compared to sham animals (sham + vehicle 95.19 ± 7.38 imp/sec vs. DED + vehicle 130.30 ± 8.66 imp/sec, *p* < 0.001).

The acute application of M8-B at 10, 20, and 30 µM on ex vivo eye preparations from sham animals did not significantly alter the cold-evoked (at 20 °C) ciliary nerve activities (sham + vehicle 95.19 ± 7.38 imp/sec vs. sham + M8-B (10 µM) 88.80 ± 12.18 imp/sec vs. sham + M8-B (20 µM) 74.80 ± 19.00 imp/sec vs. sham + M8-B (30 µM) 76.20 ± 20.27 imp/sec, *p* > 0.05; Figure 3C). In contrast, the acute treatment of the ex vivo eye preparation from DED mice with M8-B at 20 or 30 µM in superfusion medium set to 20 °C significantly blocked the cold-induced activation of nociceptors by 44.05% and 40.42%, respectively (DED + vehicle 130.30 ± 8.66 imp/sec vs. DED + M8-B (10 µM) 88.12 ± 16.13 imp/sec, *p* > 0.05; vs. DED + M8-B (20 µM) 72.90 ± 16.34 imp/sec, *p* < 0.01; vs. DED + M8-B (30 µM) 77.63 ± 15.73 imp/sec, *p* < 0.05; Figure 3C). We found no statistical difference in the cold (20 °C)-evoked activities of the ciliary nerve between sham and DED groups while treating the cornea with M8-B at 10, 20, and 30 µM (Figure 3C).

## 2.3. Chronic Topical Treatment with M8-B Reduced Corneal Cold Nociceptor Sensitization in DED Animals

To decrease cold-sensitization of DED mice after an acute application to the ex vivo eye preparation, M8-B at 20 µM was the most efficient concentration. Thus, we decided to chronically treat DED mice twice per day from day 7 until day 21 with one drop of M8-B at 20 µM or phosphate buffer saline (PBS). At day 21 post-surgery, to test corneal cold nociceptors responsiveness, we assessed electrical activities of the ciliary nerve fiber in ex vivo eye preparations using a multi-unit extracellular recording.

When the temperature of superfusion medium was set to 32–34 °C, the ongoing nerve impulse activities decreased by 23% in DED mice instilled with 20 µM M8-B compared to DED mice instilled with PBS (DED + PBS 86.78 ± 5.57 imp/sec vs. DED + M8-B (20 µM) 66.00 ± 11.43 imp/sec, *p* < 0.05; Figure 4A). When the temperature of the superfusion medium was decreased to 20 °C, the cold-evoked nerve activities reduced by 35.33% in DED mice chronically instilled with M8-B (20 µM) compared to DED mice instilled with PBS. (DED + PBS 127.71 ± 13.53 imp/sec vs. DED + M8-B (20 µM) 82.33 ± 4.27 imp/sec, *p* < 0.05; Figure 4B).
<figure><img src="image_7.png" /></figure>

**Figure 4.** Ex vivo evaluation of spontaneous and evoked activity of ciliary nerve in DED mice treated with 20 µM M8-B or PBS. Histograms show mean value of (A) spontaneous activity of ciliary nerve at 32 °C and (B) evoked activity of ciliary nerve at 20 °C in DED mice treated with 20 µM M8-B vs. PBS. DED mice treated with 20 µM M8-B vs. PBS (n = 6). All experiments were conducted at day 21.

* *p* < 0.05 relative to DED mice treated with PBS. Results are expressed as mean ± SEM. For statistical analysis, the Kolmogorov-Smirnov test was performed followed by a parametric *t*-test using GraphPad Prism version 7.00 (GraphPad Software, La Jolla, CA, USA).

## 2.4. Chronic Topical Treatment with M8-B Reduced the Corneal Inflammation in the DED-Animals

At day 21 post-surgery, we assessed the effect of chronic M8-B (20 µM) topical treatment in DED mice on corneal layer alterations (i.e., superficial epithelium, sub-basal nerve plexus, and stroma) by in vivo confocal microscopy. Figure 5A shows that no difference was detected in the corneal superficial epithelium layer between PBS- and M8-B-treated DED mice. However, we noticed that M8-B topical chronic treatment diminished the number of long thin irregular dendritic cells in the sub-basal plexus
(Figure 5B, red arrows) and refractive cells in the stroma compared with DED mice topically instilled with PBS (Figure 5C, red stars).

<figure><img src="image_8.png" /><figcaption>**Figure 5.** In vivo confocal microscopy images showing the different corneal layers from DED mice treated with 20 µM M8-B vs. PBS at day 21. (A) Corneal superficial epithelium, (B) corneal sub-basal nerve plexus (dendritic cells shown by red arrows), and (C) corneal stroma (activated cells shown red stars). Scale bar =100 µm.</figcaption></figure>

## 2.5. Chronic Topical Treatment with M8-B Reduced Inflammatory Gene Expressions in the TG of DED Animals

To investigate the impact of chronic M8-B (20 µM) topical treatment on peripheral sensitization mechanisms affecting the TG of DED mice, we assessed the expression of genes related to inflammation and neuropathic pain. Thus, the levels of tachykinin precursor 1 (TAC1) mRNA (Figure 6A) decreased significantly by 98% in M8-B-treated DED mice compared with PBS-treated DED mice (DED + PBS 3.56 ± 0.85 vs. DED + M8-B 0.05 ± 0.10, *p* < 0.01). The levels of interleukins-1β (IL1β) and interleukins 18 (IL-18) mRNA decreased by 29% and 68%, respectively, in M8-B-treated compared to PBS-treated groups (IL1β mRNA: DED + PBS 1.89 ± 0.21 vs. DED + M8-B 1.35 ± 0.05, *p* < 0.05; IL18 mRNA: DED + PBS 3.69 ± 0.92 vs. DED + M8-B 1.18 ± 0.07, *p* < 0.01; Figure 6A). We noticed that levels of the chemokine (C-C motif) ligand 2 CCL2 and CX3C chemokine receptor 1 (CX3CR1) mRNA decreased by 65% and 44%, respectively, in M8-B-treated vs. PBS-treated DED mice (CCL2 mRNA: DED + PBS 4.28 ± 0.96 vs. DED + M8-B 1.46 ± 0.09, *p* < 0.01; CX3CR1: DED + PBS 1.74 ± 0.36 vs. DED + M8-B 0.96 ± 0.97, *p* < 0.05; Figure 6A). Levels of prostaglandin E synthase 3 (PTGES3) and prostaglandin E receptor 3 (PTGER3) mRNA decreased by 62.80% and 64.79%, respectively, in M8-B-treated vs. PBS-treated DED mice (PTGES3: DED + PBS 3.98 ± 0.85 vs. DED + M8-B 1.48 ± 0.11, *p* < 0.05; PTGER3: DED + PBS 2.84 ± 0.72 vs. DED + M8-B 1.00 ± 0.21, *p* < 0.01; Figure 6A).
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Evaluation of the expression of genes implicated in conduction and modulation of pain responses in TGs of DED mice treated with 20 µM M8-B vs. PBS. (A–C) RT-qPCR analysis of ipsilateral TGs of DED mice treated with 20 µM M-8B vs. PBS at day 21. DED mice treated with 20 µM M8-B vs. PBS (n = 4–5). \* *p* < 0.05 and \*\* *p* < 0.01 relative to DED mice treated with PBS. Results are expressed as mean ± SEM. For statistical analysis, the Kolmogorov–Smirnov test was performed followed by a nonparametric Mann–Whitney test using GraphPad Prism version 7.00 (GraphPad Software, La Jolla, CA, USA).</figcaption></figure>

## 2.6. Chronic Topical Treatment with M8-B Decreased the Expression of Genes Implicated in Pain Conduction

We investigated the impact of topical M8-B on genes implicated in pain conduction in the TG of DED mice treated with M8-B (20 µM) compared with DED mice treated with PBS.

### 2.6.1. Intracellular Signal Transduction

The levels of mitogen-activated protein kinase (MAPK8) mRNA decreased by 50% in M8-B-treated DED mice compared with PBS-treated DED mice (MAPK8: DED + PBS 2.26 ± 0.58 vs. DED + M8-B 1.12 ± 0.05, *p* < 0.05; Figure 6B).

### 2.6.2. Voltage-Gated Ion Channels

The levels of the sodium channel 1.8 (Nav1.8) and sodium channel 1.7 (Nav1.7) mRNA decreased by 39% and 23%, respectively, in M8-B-treated DED mice compared with PBS-treated DED mice (Nav1.8: DED + PBS 1.29 ± 0.27 vs. DED + M8-B 0.78 ± 0.03, *p* < 0.05; Nav1.7: DED + PBS 2.23 ± 0.70 vs. DED + M8-B 1.71 ± 0.01, *p* < 0.05; Figure 6). The levels of potassium channel 2 (KCNJ6) mRNA decreased by 70% (KCNJ6: DED + PBS 3.52 ± 1.91 vs. DED + M8-B 1.03 ± 0.06, *p* < 0.05; Figure 6C).

## 2.7. Pharmacological Blockage of TRPM8 Reduced Spontaneous and Mechanical-Evoked Ocular Pain Sensations in DED Mice

We previously reported that severe DED induced corneal mechanical allodynia and spontaneous pain behavior (decreased eye closing ratio) [41]. Here, we measured the mechanical sensitivity of the cornea in DED mice chronically instilled with M8-B (20 µM) or PBS. The corneal mechanical threshold
increased by 200% in DED mice treated with 20 µM M8-B compared with DED mice treated with PBS (DED + PBS 0.01 ± 0.01 g vs. DED + M8-B (20 µM) 0.03 ± 0.01 g, *p* < 0.01; Figure 7A). We found that the spontaneous eye pain increased by 15% in DED mice treated with 20 µM of M8-B compared with DED mice treated with PBS (DED + PBS 0.7 ± 0.02 vs. DED + M8-B (20 µM) 0.8 ± 0.03, *p* < 0.05; Figure 7B). Moreover, we evaluated corneal epithelial damage using slit lamp examinations, after instillation a drop of 0.5% fluorescein, and observed superficial punctate keratitis (red arrows, Figure 7C) in the cornea of the DED mice chronically treated with M8-B (20 µM) or PBS. However, the fluorescein staining score was not significantly different between DED groups (DED + PBS 2.61 ± 0.22 vs. DED + M8-B (20 µM) 2.62 ± 1.16, *p* > 0.05; Figure 7C).

<figure><img src="image_10.png" /><figcaption>**Figure 7.** Mechanical allodynia and corneal epithelial damage in DED mice treated with M8-B (20 µM) vs. PBS. (A) Corneal mechanical sensitivity measured with von Frey filaments (g: gram). (B) Images and quantification of spontaneous eye closing ratio calculated by measuring the width: height ratio. (C) Cornea alteration was examined using a slit-lamp microscope with a cobalt blue filter after placing a drop of fluorescein (0.5% Fluorescein Faure) in the right eye of the animals. The punctuate keratitis (red arrows) was quantified using standard grading scale, 21 days after surgery in DED mice treated with M-8B (20 µM) vs. PBS (n = 6–9). All tests were conducted at day 21. \* *p* < 0.05 and \*\* *p* < 0.01 relative to DED mice treated with PBS. Results are expressed as mean ± SEM. For statistical analysis, the Kolmogorov–Smirnov test was performed followed by a parametric *t*-test using GraphPad Prism version 7.00 (GraphPad Software, La Jolla, CA, USA).</figcaption></figure>
### 3. Discussion

DED patients suffer from abnormal hypersensitivity to normally harmless cold stimuli which is known as cold allodynia [42]. TRPM8 is a key channel for cold and pain perception [10,43–45]. Previous studies reported the effectiveness of M8-B, a selective and potent TRPM8 antagonist [46–48], on cold nociceptors via TRPM8 by decreasing deep body temperature in *Trpm8⁺⁺* mice and rats, but not in *Trpm8⁻⁻* mice, highlighting an on-target action [48]. Thus, in the present study, we explored whether M8-B could alleviate ocular pain associated with severe DED.

The relevant contribution of our study is the demonstration that chronic topical instillations of M8-B in mice suffering of severe DED blocks: (i) cold hyperresponsiveness of corneal nerves, (ii) corneal inflammation, and (iii) corneal mechanical allodynia and spontaneous ocular pain. So, our results confirmed the pivotal role of TRPM8 in acute pain and that the TRPM8 antagonist could be a promising treatment to alleviate chronic ocular pain caused by severe DED.

These results were obtained in adult male mice submitted to severe DED after the removal of two different functional (aqueous and lipid) glands which cause a drastic and permanent decrease of tear production by 97% [41]. Our DED model closely mimic both aqueous-deficient DED by excision of the extraorbital lacrimal gland (ELG) and evaporative DED by excision of the Harderian gland (HG), which produces an oily, lipid-enriched secretion [41]. The severe DED induced by this surgical procedure may be considered as a limitation since we cannot directly test the effectiveness of TRPM8 antagonist treatment on tear production. However, we developed a model of chronic DED associated with increased ocular nociception without the need to constantly inject chemical solutions with neurotropic activity (such as scopolamine) or to place mice in a desiccative stressful environment. Otherwise, a comparative study between adult female mice submitted to severe DED [41] with male DED and sham mice was previously performed on behavioral and electrophysiological responses related to the present study and no sex differences were noticed (data not shown).

Cool cells activity is increased by menthol, a specific TRPM8 agonist [18,19]. Here, we observed that DED mice developed a higher corneal menthol sensitivity (cold allodynia), which can be explained by the upregulation of TRPM8 mRNA expression observed in the ophthalmic branch of the TG. However, Hatta et al. showed that TRPM8 mRNA levels were not affected in corneal trigeminal neurons from a DED rat model obtained 10 days after extraorbital and infraorbital lacrimal glands removal [49]. In addition, TRPM8 protein levels in the cornea, conjunctiva, and TG were the same in sham animals and moderate DED rat model obtained 14 days after ELG excision [50]. This discrepancy of our results may be explained by (1) using a different DED model obtained by removing two functionally different glands (aqueous and lipid) vs. only aqueous glands. The excision of two functionally different glands at day 21 decreased tear production by 97% in our model [41] compared to 30% reported 10 days after ELG and infraorbital LG excision in rats [33] and 50% after 14 days ELG excision in rats [51], highlighting the severity of our DED model compared to moderate models of DED; (2) using post-surgery evaluation kinetics for the evaluation of TRPM8 expression at day 21 vs. day 10 or day 14 post ELG excision; and (3) monitoring TRPM8 mRNA expression in the ophthalmic branch of TG (responsible for both cornea and conjunctiva innervation) vs. only retrograde-labelled corneal trigeminal neurons. Altogether, our findings highlight an epigenetic regulation of the TRPM8 channel expression in the ophthalmic branch of TG during severe DED.

Ongoing ciliary nerve fiber electrical activity mostly results from cold nerve fiber activity, which represents half of the corneal sensory neurons [52,53]. Thus, the increase in ongoing ciliary nerve fiber activity in our severe mice model of DED may result from the triggering of corneal cold nociceptors, as was also reported in DED guinea pigs obtained after ELG removal [32]. Corneal cold-thermoreceptor fibers are classified into two groups: high-background low-threshold cold thermoreceptors (HB-LTCs) have a background firing rate at the normal temperature of the cornea (~34 °C), which increases with small temperature decreases (<1 °C); and low-background high-threshold cold thermoreceptors (LB-HTCs), which need large temperature reductions to increase their firing rate [15,54]. HB-LTCs control basal tearing and blinking rates, whereas LB-HTCs evoke sensations associated with corneal
dryness [55]. Whereas acute topical treatment with M8-B (20 and 30 µM) decreased only the spontaneous ciliary nerve fiber activity at 32 °C, we observed that 14 days' treatment by 20 µM M8-B significantly decreased the spontaneous and cold-evoked ciliary nerve fiber activity in DED mice. These observations confirm the effectiveness of M8-B on cold responsiveness by acting on both HB-LTCs and LB-HTCs, which may lead to modulating blinking rate and corneal dryness in DED mice. Our anatomical and electrophysiological results unveiled TRPM8-dependent changes in corneal neurons during severe DED, suggesting that increased ongoing and evoked nerve activities in DED animals are modulated via TRPM8 channels.

Ocular surface alterations and inflammation are the core mechanism of DED [56]. We previously found corneal superficial punctate keratitis and an increase in inflammatory cells in the sub-basal plexus and corneal stroma in DED animals [41]. In this study, we demonstrated that chronic pharmacological blockade of TRPM8 by chronic instillations of M8-B considerably decreased corneal inflammation in DED animals. However, this topical treatment was not effective in reducing severe ocular surface epitheliopathy observed in DED mice. This can be explained by the radical reduction in tear secretion (by 97%) and the lack of lubrication causing chronic mechanical alterations on the corneal epithelial integrity.

TRPM8 has been implicated in inflammation in a model of asthmatic mice where cold stimulus promoted inflammation via TRPM8 [57,58]. Thus, based on our observations and literature findings, we suggest a link between TRPM8 and corneal inflammation associated with DED. The anti-inflammatory effect of TRPM8 blockade was not limited to the cornea; it also led to changes in the ipsilateral TG. A drastic decrease in the immune receptors and inflammatory mediator's mRNA levels of IL-1β, IL-18, CCL2, prostaglandins, and CX3CR1 was observed in the TG of DED animals chronically instilled with M8-B. These receptors and inflammatory mediators are known to be increased in DED and have been implicated in the spread of inflammation from cornea to the TG [59–61]. Thus, our results clearly suggest that M8-B can decrease the spread of inflammation from the cornea to TG in DED animals via inhibition of TRPM8.

Sodium and potassium currents have been associated with cold responses [62–64]. In the context of DED, Kovacs et al. reported that ocular cold nociceptor activities are mediated by enhanced Na⁺ currents combined with diminished K⁺ currents in trigeminal corneal cold sensory neurons in a DED guinea pig model produced after ELG removal [32]. Our findings showed decreases in Nav1.8, Nav1.7, and KCNJ6 channel mRNA levels in the TG from DED mice treated with M8-B. Since Na⁺ and K⁺ channels play a role in ocular cold responsiveness, our result showing that TRPM8 antagonist decreases the gene expressions of Na⁺ and K⁺ channels in the TG highlights the effectiveness of the TRPM8 antagonist in blocking the exacerbated cold responses of the cornea observed in DED-mice.

Substance P (SP) was found to be required for corneal cold nociception [65]. SP is encoded by the *Tac1* gene and is known to modulate neurogenic inflammation [66,67]. Here, we report that M8-B decreased the *Tac1* mRNA level in the TG of DED animals, underlining a reduction in neurogenic inflammation and the excitation of peripheral free terminals in DED animals. This hypothesis is supported by the decrease in intracellular signaling pathway MAPK8 mRNA levels observed in the TG of DED animals treated with M8-B, emphasizing a lesser propagation of nociceptive messages.

Although much attention has been paid to the role of TRPM8 in cold allodynia and hyperalgesia, some studies have explored the role of TRPM8 in mechanical allodynia. Klein et al. reported that topical application of menthol induced mechanical allodynia in wild type (WT) rats [68], whereas TRPM8 inhibitors were shown to reduce cold and mechanical allodynia in acute and chronic pain models [29,30]. In line with these studies, we observed that chronic M8-B instillations decreased the mechanical corneal hypersensitivity and spontaneous eye pain developed in DED animals. One explanation for the anti-mechanical allodynic effect of TRMP8 antagonist is possible crosstalk between TRPM8 and a polymodal channel (TRPV1) that can be activated by mechanical stimulations. In agreement with our hypothesis, Hatta et al. reported increases in TRPM8 and TRPV1 co-expressions in the TG neurons of DED rats obtained after extraorbital and infraorbital lacrimal glands removal [49]. A
recent study demonstrated an overexpression of TRPV1 in TRPM8+ cold-sensing fibers, which caused cold allodynia in a mouse model of DED obtained by ELG excision [65]. For the first time, this finding showed that TRPV1 activity is required for corneal cold nociception. Altogether, we suggest that TRPM8 blockage influences TRPV1 activity and leads to a decrease in corneal mechanical hypersensitivity.

# 4. Methods

## 4.1. *Experimental Animals*

Seven- to eight-week-old adult male C57BL/6 mice (average weight 23.48 ± 0.04 g; Janvier Labs, Le Genest Saint Isle, France) were maintained under controlled conditions (22 ± 1 °C, 60% ± 10% relative humidity, 12/12 h (hour) light/dark cycle, food and water ad libitum). All animal procedures were performed in strict accordance with institutional guidelines for the care and use of experimental animals approved by the European Communities Council Directive 2010/63/UE (APAFIS #1501 2015081815454885 v2, 30 September 2015).

## 4.2. *Surgical Procedures*

Unilateral (right side) ELG and HG excision was performed under ketamine (80 mg/kg intraperitoneal (i.p.) and xylazine (8 mg/kg i.p.) anesthesia, as recently described [41]. Before surgery, a drop of lacrimal gel (Lubrithal™, Hamilton, Canada) was applied to both eyes. Under an operative microscope (Leica-Alcon II, Wetzlar, Germany), an 8 mm skin incision was made on the temporal side to expose and remove the ELG. After dissociating the conjunctival tissue above the orbital cavity near the internal canthus, the HG was carefully removed. Complete removal was verified by inspecting the surgical area for any remaining glandular tissue. The skin incision was then sutured using 6.0 braided silk sutures (6-0 Vicryl, Ethicon, Dulmen, Germany). A drop of iodine solution was applied to the incision to avoid bacterial infection. For sham mice, an incision was made in the same zone without touching the glands. The mice were placed in warm (30 °C) cages to recover from the surgery.

## 4.3. *Drugs*

M8-B hydrochlorid: N-(2-aminoethyl)-N-(4-(benzyloxy)-3-methoxybenzyl)thiophene-2-carboxamide hydrochloride, N-(2-Aminoethyl)-N-[[3-methoxy-4-phenylmethoxy) phenyl] methyl]-2-thiophenecarboxamide hydrochloride) (Scheme 1) (SML0893-5MG, Sigma-Aldrich, St. Louis, MO, USA) was dissolved in sterile distilled water at the concentration of 2 mg/mL and then diluted in phosphate-buffered saline (PBS) 1× to obtain 10, 20, and 30 µM solutions.

<figure><img src="image_11.png" /><figcaption>**Scheme 1.** M8-B hydrochlorid batch molecular structure.</figcaption></figure>

Batch Molecular Formula: C22H24N2O3S·HCl. Batch Molecular Weight: 432.96. HPLC: Shows ≥98% purity. M8-B is a selective and potent antagonist of the TRPM8 channel. M8-B blocked cold-induced and icilin or menthol-induced activation of rat, human, and murine TRPM8 channels with the half maximal IC50 values ranging from 7.8 to 64.3 nM, and it did not block other TRP channels (IC50 values > 20 mM).

Menthol (Sigma-Aldrich, St. Louis, Missouri, United States) was diluted in pure castor oil (Sigma-Aldrich, St. Louis, Missouri, United States) to obtain a 50 µM solution.
## 4.4. *Experimental Study*

Twice per day, DED mice received ocular instillations of M8-B (20 µM) or PBS from days 7 until 21. For topical ocular administration, the cage order was randomized daily, and the experimenter was blinded to the treatment group. The same experimenter performed all experiments, administrated the study treatment, and assessed the effect. The experimenter was blinded to the treatment group.

## 4.5. *Behavioral Tests*

For behavioral tests, mice were placed in the testing room (at least 30 min before the start of the experiments) and the order of testing was randomized. The immobilization (handling-restraint) lasted a maximum of 20 s to reduce any potential stress and anxiety. Behavioral experiments were conducted at day 21.

## 4.6. *Measurement of Corneal Sensitivity to Mechanical and Cold Stimulation*

Mechanical corneal sensitivity was monitored using von Frey filaments, as previously described [41]. Various forces of calibrated von Frey filaments (0.008 to 0.04 g) were applied to the center of the cornea of manually immobilized mice. The mechanical threshold corresponded to the eye-blink response. For corneal cold sensitivity, 10 µL of a 50 µM menthol solution was applied to the right eye. Animals were immediately placed in individual cages and the eye closing ratio was measured 5 min after menthol ocular application.

## 4.7. *Spontaneous Pain: Eye Closing Ratio Measurement*

Spontaneous eye closure is an accurate index used for monitoring spontaneous eye pain [41] and is one of the quantitative measures of the grimace scale, which is used to monitor spontaneous pain behavior [69,70]. The eye closing ratio was calculated based on photographs of awake and unconstrained mice and corresponds to the height/width ratio of palpebral closure. Width is the distance between the internal and external canthus and height is the distance between the edge of the upper and lower eyelids, going through the center of the cornea. Images were captured by a digital camera using EyeSuite™ software (Bern, Switzerland).

## 4.8. *Fluorescein Staining Score*

The corneas of lightly anesthetized mice were washed with sterile 0.9% NaCl. Then a drop of fluorescein (Fluorescein Faure 0.5%) was placed into the conjunctival sac of the right eye. Ocular surface staining was then evaluated using a slit lamp under a cobalt blue light (peak of approximately 400 nm). Images were captured by a digital camera using EyeSuite™ software (Koeniz, Switzerland). The extent of corneal damage (width and intensity of fluorescein take-up area) was scored according to the following scale: 0, absence of staining; 0.5, slight punctate staining; 1, diffuse punctate staining; 2, diffuse staining covering less than one third of the cornea; 3, diffuse staining covering more than one third of the cornea; and 4, staining covering more than two thirds of the cornea [41,71].

## 4.9. *In Vivo Confocal Microscopy*

An in vivo laser confocal microscope (IVCM, Heidelberg Retina Tomography (HRT) II/Rostock Cornea Module (RCM); Heidelberg Engineering GmbH, Heidelberg, Germany)) was used to examine the entire cornea of anesthetized mice [72]. The images covered a 400 × 400 µm area with a transversal optical resolution of 2 µm and an axial optical resolution of 4 µm (Heidelberg Engineering, Heidelberg, Germany). Approximately 100 images were acquired per animal. Image acquisitions always started in the center of the cornea at the level of the superficial epithelium. Peripheral acquisitions were then performed similarly.
## 4.10. Multi-Unit Extracellular Recording of Evoked Ciliary Nerve Fiber Activity in Ex Vivo Eye Preparations

Spontaneous ciliary nerve fiber activity was determined at day 21 as previously reported [41,72]. The spontaneous ciliary nerve activity's recordings were performed by superfusing the cornea with superfusion saline solution at 33 ± 1 °C and pH 7. Briefly, mice were euthanized and the eye was placed in a two-compartment chamber. The cornea was continuously superfused at a rate of 3 mL/min at 33 ± 1 °C with a superfused saline solution (133.4 mM NaCl, 4.7 mM KCl, 2 mM CaCl₂, 1.2 mM MgCl₂, 16.3 mM NaHCO₃, 1.3 mM NaH₂PO₄, and 7.8 mM glucose) saturated with O₂ and adjusted to pH 7.4 by bubbling with 95% O₂ and 5% CO₂. Multi-unit extracellular electrical activity of the ciliary nerve was recorded using a suction electrode (Ag/AgCl). The signal was filtered (300–5000 Hz), amplified (10,000×) (A-M Systems, Sequim, WA, USA), and digitalized by Spike 2 data analysis (CED Micro1401, Cambridge Electronic Design, Cambridge CB24 6AZ, United Kingdom) at a sampling frequency of 10,000 Hz. The cornea was superfused with superfusion saline solution for 30 min to stabilize the preparation before performing electrophysiological recordings. Extracellular spontaneous ciliary nerve fiber activity is defined as impulses per second (imp/s).

The ciliary nerve evoked activity's recordings were performed by superfusing the cornea with superfusion saline solution at 20 ± 1 °C and pH 7. A thermistor sensor (included in the CL-100 Bipolar Temperature Controller, Warner Instruments, CT, USA) monitored the temperature at the exit of the corneal superfusion. The acute application of M8-B was evaluated by superfusion of the cornea with 10, 20, and 30 µM M8-B diluted in the superfusion saline solution.

## 4.11. RT-qPCR Analysis

### 4.11.1. Tissue Preparation for RT-qPCR Analysis

Twenty-one days after surgery, the animals were deeply anesthetized with a 300 µL mixture of ketamine (80 mg/kg) and xylazine (8 mg/kg) and transcardially perfused with 10 mL 0.9% NaCl solution. Ipsilateral TGs were carefully removed and stored at -80 °C.

### 4.11.2. RT-qPCR Analysis

RNA extraction from the ipsilateral TG was performed with a NucleoSpin RNA Purification II Kit (NucleoSpin RNA S, Duren, Deutschland, Germany). RNA quality and concentration were then measured by the NanoDrop method (Thermo Scientific, Waltham, MA, USA). Then, reverse transcription was performed by iScript cDNA Synthesis Kit (Bio-Rad, Hercules, CA, USA) according to the manufacturer's instructions. PCR was performed with 300 ng cDNA for each simple. RT-qPCR was performed with SsoAdvanced Universal SYBR® Green Supermix (Bio-Rad, Hercules, CA, USA), and a pain, neuropathic, and inflammatory (SAB Target List) M96-well plate (Bio-Rad, Hercules, CA, USA; ref: 10034393). The *GAPDH* gene was used as the endogenous reference for each reaction; mRNA levels were calculated after normalizing the results for each sample with those for *GAPDH* mRNA. The 2⁻ΔΔCt method was used to analyze the relative differences in specific mRNA levels between groups.

## 4.12. *Tissue Preparation for Fluorescent In Situ Hybridization*

Twenty-one days after surgery, anesthetized mice were transcardially perfused with 10 mL 0.9% NaCl solution followed by 40 mL 4% (*w*/*v*) paraformaldehyde in PBS. Next, TGs were immersed in 10%, 20%, and 30% sucrose in PBS and then conserved in isopentane with liquid nitrogen and stored at -80 °C. TGs (12 µm) were cut on a cryostat (Leica CM 3050 S, Wetzlar, Germany) and mounted on Superfrost slides (Thermofisher scientific, Waltham, MA, USA).

## 4.13. Fluorescent RNAscope In Situ Hybridization

Fluorescent in situ hybridization studies was performed according to the protocol for fixed frozen tissue using the RNAscope Fluorescent Multiplex Reagent Kit v2 assay (Advanced Cell Diagnostics,
Newark, CA, USA). Tissues were washed with PBS and treated with hydrogen peroxide (RNAscope, Bio-techne, Minneapolis, Minnesota, United States, ref# 322335) for 10 min at RT and washed in autoclaved distilled water. Using a steamer, tissues were treated with distilled H₂O for 10 s at 99 °C and then moved to RNAscope 1x target Retrieval Reagent (RNAscope, Bio-techne, Minneapolis, Minnesota, USA, ref# 322000) for 5 min at 99 °C. Tissues were washed with autoclaved distilled water and transferred to 100% alcohol for 3 min. After, tissues were treated with RNAscope Protease III (RNAscope, Bio-techne, Minneapolis, MA, United States, ref#322337) for 30 min at 40 °C and washed with autoclaved distilled water. Species-specific target probe TRPM8-C3 (420451 C3) was used. Sections were treated with the probe; negative (ref# 320871) and positive (ref# 320881) controls and were hybridized for 2 h at 40 °C in a humidified oven (RNAscope, Bio-techne, Minneapolis, MA, USA, HybEZ oven with HybEZ humidity control tray). A series of incubations were then performed to amplify the hybridized probe signals and label target probes with the assigned fluorescence detection channel (target probe was labelled with the assigned fluorescence detection channel (PerkinElmer, Waltham, MA, USA. Opal 650 FP1496A). Nuclei were stained using a DAPI nuclear stain (RNAscope, Bio-techne, Minneapolis, MA, USA, Ref# 323108) for 30 s at RT. Slides were mounted with ProLong Gold Antifade Mountant (Thermofisher scientific, Waltham, MA, USA, Ref# P36934) onto glass sides and cover slipped.

### 4.14. Microscopic and Nanozoomer Analysis

TGs were examined using a nanozoomer. The area within the field of interest covered by the TRPM8 mRNA profiles relative to the total area of the measured field was analyzed to detect mRNA levels in the ophthalmic branch in the TG. The same gray threshold level was applied to all sections of the same series. TIFF images were analyzed using NIH ImageJ software.

### 4.15. Statistical Analyses

The data obtained in different groups were compared using the appropriate paired parametric or nonparametric statistical test, as indicated. For statistical analysis, the Kolmogorov-Smirnov test was performed, followed by the parametric t-test or nonparametric Mann-Whitney or parametric one-way ANOVA test (GraphPad Software, La Jolla, CA, USA). All *p*-values were considered statistically significant when values were <0.05. All results are presented as the mean ± standard error of the mean (SEM).

## 5. Conclusions

Overall, our study showed that pharmacological blockade of corneal-expressed TRPM8 (by chronic instillations of M8-B) significantly reduced corneal neurosensory abnormalities in a preclinical model of DED, thus reinforcing the pharmacological concept that TRPM8 antagonists may be of use as a local antalgic to relieve ocular pain in cases of severe DED.

**Author Contributions:** D.F. performed the experiments, analyzed and interpreted the data, generated the figures, discussed the results, and wrote the manuscript. C.B. interpreted the data, discussed the results and their significance, and commented on the manuscript. S.M.P. and A.R.-L.G. designed and supervised all experiments, interpreted the data, discussed the results and their significance, and wrote the manuscript. All authors approved the final article. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was funded by Laboratoire Théa, within the framework of research collaboration with Sorbonne Université/Institut de la Vision. Darine Fakih is supported by the Association Nationale de la Recherche et de la Technologie (Conventions Industrielles de Formation par la Recherche N° 2016/1641). This work was completed with the support of the Programme Investissements d'Avenir IHU FORESIGHT (ANR-18-IAHU-01).

**Acknowledgments:** The authors thank the imagery platform and the staff of the animal house facilities at the Institut de la Vision for their help. We also thank Sarah Gallet-Wimez from UMR-S1172 development and plasticity of the neuroendocrine brain, Lille, France, for RNAscope in situ hybridization training.
**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

**Abbreviations**

<table><tr><td>DED</td><td>Dry eye disease</td></tr><tr><td>TRPM8</td><td>Transient receptor potential melastatin 8</td></tr><tr><td>TG</td><td>Trigeminal ganglion</td></tr><tr><td>Min</td><td>Minutes</td></tr><tr><td>SEM</td><td>Standard error of the mean</td></tr><tr><td>Imp</td><td>Impulse</td></tr><tr><td>Sec</td><td>Seconds</td></tr><tr><td>TAC1</td><td>Tachykinin Precursor 1</td></tr><tr><td>IL-1β</td><td>Interleukin-1β</td></tr><tr><td>IL18</td><td>Interleukin 18</td></tr><tr><td>CCL2</td><td>Chemokine (C-C motif) ligand 2</td></tr><tr><td>CX3CR1</td><td>CX3C chemokine receptor 1</td></tr><tr><td>PTGES3</td><td>Prostaglandin E synthase 3</td></tr><tr><td>PTGER3</td><td>Prostaglandin E receptor 3</td></tr><tr><td>MAPK</td><td>Mitogen-activated protein kinase</td></tr><tr><td>Nav1.7</td><td>Sodium channel 1.7</td></tr><tr><td>Nav1.8</td><td>Sodium channel 1.8</td></tr><tr><td>KCNJ6</td><td>Potassium channel 2</td></tr><tr><td>g</td><td>Gram</td></tr><tr><td>ELG</td><td>Extraorbital lacrimal gland</td></tr><tr><td>HG</td><td>Harderian gland</td></tr><tr><td>HB-LTCs</td><td>high-background low-threshold cold thermoreceptors</td></tr><tr><td>LB-HTCs</td><td>Low-background high-threshold cold thermoreceptors</td></tr><tr><td>TRPV1</td><td>The transient receptor potential vanilloid 1</td></tr><tr><td>h</td><td>Hour</td></tr><tr><td>i.p.</td><td>Intraperitoneal</td></tr><tr><td>RT</td><td>Room temperature</td></tr></table>

**References**

1. Craig, J.P.; Nichols, K.K.; Akpek, E.K.; Caffery, B.; Dua, H.S.; Joo, C.K.; Liu, Z.; Nelson, J.D.; Nichols, J.J.; Tsubota, K.; et al. TFOS DEWS II Definition and Classification Report. *Ocul. Surf.* **2017**, *15*, 276–283. [CrossRef] [PubMed]

2. Stapleton, F.; Alves, M.; Bunya, V.Y.; Jalbert, I.; Lekhanont, K.; Malet, F.; Na, K.-S.; Schaumberg, D.; Uchino, M.; Vehof, J.; et al. TFOS DEWS II Epidemiology Report. *Ocul. Surf.* **2017**, *15*, 334–365. [CrossRef] [PubMed]

3. Belmonte, C.; Nichols, J.J.; Cox, S.M.; Brock, J.A.; Begley, C.G.; Bereiter, D.A.; Dartt, D.A.; Galor, A.; Hamrah, P.; Ivanusic, J.J.; et al. TFOS DEWS II pain and sensation report. *Ocul. Surf.* **2017**, *15*, 404–437. [CrossRef] [PubMed]

4. Goyal, S.; Hamrah, P. Understanding Neuropathic Corneal Pain-Gaps and Current Therapeutic Approaches. *Semin. Ophthalmol.* **2016**, *31*, 59–70. [CrossRef]

5. Cruzat, A.; Qazi, Y.; Hamrah, P. In Vivo Confocal Microscopy of Corneal Nerves in Health and Disease. *Ocul. Surf.* **2017**, *15*, 15–47. [CrossRef] [PubMed]

6. Marfurt, C.F.; Kingsley, R.E.; Echtenkamp, S.E. Sensory and sympathetic innervation of the mammalian cornea. A retrograde tracing study. *Investig. Ophthalmol. Vis. Sci.* **1989**, *30*, 461–472.

7. Marfurt, C.F.; Del Toro, D.R. Corneal sensory pathway in the rat: A horseradish peroxidase tracing study. *J. Comp. Neurol.* **1987**, *261*, 450–459. [CrossRef]

8. Belmonte, C.; Acosta, M.C.; Merayo-Lloves, J.; Gallar, J. What Causes Eye Pain? *Curr. Ophthalmol. Rep.* **2015**, *3*, 111–121. [CrossRef]
9. Hirata, H.; Meng, I.D. Cold-sensitive corneal afferents respond to a variety of ocular stimuli central to tear production: Implications for dry eye disease. *Investig. Ophthalmol. Vis. Sci.* **2010**, *51*, 3969–3976. [CrossRef]

10. Parra, A.; Madrid, R.; Echevarria, D.; del Olmo, S.; Morenilla-Palao, C.; Acosta, M.C.; Gallar, J.; Dhaka, A.; Viana, F.; Belmonte, C. Ocular surface wetness is regulated by TRPM8-dependent cold thermoreceptors of the cornea. *Nat. Med.* **2010**, *16*, 1396–1399. [CrossRef]

11. Quallo, T.; Vastani, N.; Horridge, E.; Gentry, C.; Parra, A.; Moss, S.; Viana, F.; Belmonte, C.; Andersson, D.A.; Bevan, S. TRPM8 is a neuronal osmosensor that regulates eye blinking in mice. *Nat. Commun.* **2015**, *6*, 7150. [CrossRef] [PubMed]

12. Belmonte, C.; Gallar, J. Cold thermoreceptors, unexpected players in tear production and ocular dryness sensations. *Investig. Ophthalmol. Vis. Sci.* **2011**, *52*, 3888–3892. [CrossRef]

13. González-Muñiz, R.; Bonache, M.A.; Martín-Escura, C.; Gómez-Monterrey, I. Recent Progress in TRPM8 Modulation: An Update. *Int. J. Mol. Sci.* **2019**, *20*, 2618. [CrossRef]

14. Alcalde, I.; Íñigo-Portugués, A.; González-González, O.; Almaraz, L.; Artime, E.; Morenilla-Palao, C.; Gallar, J.; Viana, F.; Merayo-Lloves, J.; Belmonte, C. Morphological and functional changes in TRPM8-expressing corneal cold thermoreceptor neurons during aging and their impact on tearing in mice. *J. Comp. Neurol.* **2018**, *526*, 1859–1874. [CrossRef]

15. Piña, R.; Ugarte, G.; Campos, M.; Íñigo-Portugués, A.; Olivares, E.; Orio, P.; Belmonte, C.; Bacigalupo, J.; Madrid, R. Role of TRPM8 Channels in Altered Cold Sensitivity of Corneal Primary Sensory Neurons Induced by Axonal Damage. *J. Neurosci. Off. J. Soc. Neurosci.* **2019**, *39*, 8177–8192. [CrossRef]

16. Huang, D.; Li, S.; Dhaka, A.; Story, G.M.; Cao, Y.-Q. Expression of the transient receptor potential channels TRPV1, TRPA1 and TRPM8 in mouse trigeminal primary afferent neurons innervating the dura. *Mol. Pain* **2012**, *8*, 66. [CrossRef]

17. Colburn, R.W.; Lubin, M.L.; Stone Jr, D.J.; Wang, Y.; Lawrence, D.; D'Andrea, M.R.; Brandt, M.R.; Liu, Y.; Flores, C.M.; Qin, N. Attenuated cold sensitivity in TRPM8 null mice. *Neuron* **2007**, *54*, 379–386. [CrossRef]

18. Dhaka, A.; Murray, A.N.; Mathur, J.; Earley, T.J.; Petrus, M.J.; Patapoutian, A. TRPM8 is required for cold sensation in mice. *Neuron* **2007**, *54*, 371–378. [CrossRef]

19. Bautista, D.M.; Siemens, J.; Glazer, J.M.; Tsuruda, P.R.; Basbaum, A.I.; Stucky, C.L.; Jordt, S.-E.; Julius, D. The menthol receptor TRPM8 is the principal detector of environmental cold. *Nature* **2007**, *448*, 204–208. [CrossRef]

20. Zimmermann, M. Pathobiology of neuropathic pain. *Eur. J. Pharmacol.* **2001**, *429*, 23–37. [CrossRef]

21. Jensen, T.S.; Gottrup, H.; Sindrup, S.H.; Bach, F.W. The clinical picture of neuropathic pain. *Eur. J. Pharmacol.* **2001**, *429*, 1–11. [CrossRef]

22. Marwaha, L.; Bansal, Y.; Singh, R.; Saroj, P.; Bhandari, R.; Kuhad, A. TRP channels: Potential drug target for neuropathic pain. *Inflammopharmacology* **2016**, *24*, 305–317. [CrossRef] [PubMed]

23. Caspani, O.; Zurborg, S.; Labuz, D.; Heppenstall, P.A. The Contribution of TRPM8 and TRPA1 Channels to Cold Allodynia and Neuropathic Pain. *PLOS ONE* **2009**, *4*, e7383. [CrossRef]

24. Frederick, J.; Buck, M.E.; Matson, D.J.; Cortright, D.N. Increased TRPA1, TRPM8, and TRPV2 Expression in Dorsal Root Ganglia by Nerve Injury. Available online: http://pubmed.ncbi.nlm.nih.gov/17517374/ (accessed on 31 August 2020).

25. Xing, H.; Chen, M.; Ling, J.; Tan, W.; Gu, J.G. TRPM8 mechanism of cold allodynia after chronic nerve injury. *J. Neurosci. Off. J. Soc. Neurosci.* **2007**, *27*, 13680–13690. [CrossRef]

26. Koh, W.U.; Choi, S.-S.; Kim, J.H.; Yoon, H.J.; Ahn, H.-S.; Lee, S.K.; Leem, J.G.; Song, J.G.; Shin, J.W. The preventive effect of resiniferatoxin on the development of cold hypersensitivity induced by spinal nerve ligation: Involvement of TRPM8. *BMC Neurosci.* **2016**, *17*, 38. [CrossRef]

27. Knowlton, W.M.; Daniels, R.L.; Palkar, R.; McCoy, D.D.; McKemy, D.D. Pharmacological Blockade of TRPM8 Ion Channels Alters Cold and Cold Pain Responses in Mice. *PLoS ONE* **2011**, *6*, [CrossRef] [PubMed]

28. Parks, D.J.; Parsons, W.H.; Colburn, R.W.; Meegalla, S.K.; Ballentine, S.K.; Illig, C.R.; Qin, N.; Liu, Y.; Hutchinson, T.L.; Lubin, M.L.; et al. Design and optimization of benzimidazole-containing transient receptor potential melastatin 8 (TRPM8) antagonists. *J. Med. Chem.* **2011**, *54*, 233–247. [CrossRef]

29. Sałat, K.; Filipek, B. Antinociceptive activity of transient receptor potential channel TRPV1, TRPA1, and TRPM8 antagonists in neurogenic and neuropathic pain models in mice. *J. Zhejiang Univ. Sci. B* **2015**, *16*, 167–178. [CrossRef]
30. De Caro, C.; Russo, R.; Avagliano, C.; Cristiano, C.; Calignano, A.; Aramini, A.; Bianchini, G.; Allegretti, M.; Brandolini, L. Antinociceptive effect of two novel transient receptor potential melastatin 8 antagonists in acute and chronic pain models in rat. *Br. J. Pharmacol.* **2018**, *175*, 1691–1706. [CrossRef]

31. De Caro, C.; Cristiano, C.; Avagliano, C.; Bertamino, A.; Ostacolo, C.; Campiglia, P.; Gomez-Monterrey, I.; La Rana, G.; Gualillo, O.; Calignano, A.; et al. Characterization of New TRPM8 Modulators in Pain Perception. *Int. J. Mol. Sci.* **2019**, *20*, 5544. [CrossRef]

32. Kovács, I.; Luna, C.; Quirce, S.; Mizerska, K.; Callejo, G.; Riestra, A.; Fernández-Sánchez, L.; Meseguer, V.M.; Cuenca, N.; Merayo-Lloves, J.; et al. Abnormal activity of corneal cold thermoreceptors underlies the unpleasant sensations in dry eye disease. *Pain* **2016**, *157*, 399–417. [CrossRef] [PubMed]

33. Kurose, M.; Meng, I.D. Dry eye modifies the thermal and menthol responses in rat corneal primary afferent cool cells. *J. Neurophysiol.* **2013**, *110*, 495–504. [CrossRef] [PubMed]

34. Hirata, H.; Oshinsky, M.L. Ocular dryness excites two classes of corneal afferent neurons implicated in basal tearing in rats: Involvement of transient receptor potential channels. *J. Neurophysiol.* **2012**, *107*, 1199–1209. [CrossRef] [PubMed]

35. Robbins, A.; Kurose, M.; Winterson, B.J.; Meng, I.D. Menthol activation of corneal cool cells induces TRPM8-mediated lacrimation but not nociceptive responses in rodents. *Investig. Ophthalmol. Vis. Sci.* **2012**, *53*, 7034–7042. [CrossRef] [PubMed]

36. Weyer, A.D.; Lehto, S.G. Development of TRPM8 Antagonists to Treat Chronic Pain and Migraine. *Pharmaceuticals* **2017**, *10*, 37. [CrossRef]

37. Yang, J.M.; Li, F.; Liu, Q.; Rüedi, M.; Wei, E.T.; Lentsman, M.; Lee, H.S.; Choi, W.; Kim, S.J.; Yoon, K.C. A novel TRPM8 agonist relieves dry eye discomfort. *BMC Ophthalmol.* **2017**, *17*, 101. [CrossRef]

38. Wasner, G.; Schattschneider, J.; Binder, A.; Baron, R. Topical menthol—A human model for cold pain by activation and sensitization of C nociceptors. *Brain* **2004**, *127*, 1159–1171. [CrossRef] [PubMed]

39. Hatem, S.; Attal, N.; Willer, J.-C.; Bouhassira, D. Psychophysical study of the effects of topical application of menthol in healthy volunteers. *Pain* **2006**, *122*, 190–196. [CrossRef] [PubMed]

40. Acosta, M.C.; Belmonte, C.; Gallar, J. Sensory experiences in humans and single-unit activity in cats evoked by polymodal stimulation of the cornea. *J. Physiol.* **2001**, *534*, 511–525. [CrossRef] [PubMed]

41. Fakih, D.; Zhao, Z.; Nicolle, P.; Reboussin, E.; Joubert, F.; Luzu, J.; Labbé, A.; Rostène, W.; Baudouin, C.; Mélik Parsadaniantz, S.; et al. Chronic dry eye induced corneal hypersensitivity, neuroinflammatory responses, and synaptic plasticity in the mouse trigeminal brainstem. *J. Neuroinflamm.* **2019**, *16*, 268. [CrossRef] [PubMed]

42. Situ, P.; Simpson, T.L.; Fonn, D.; Jones, L.W. Conjunctival and corneal pneumatic sensitivity is associated with signs and symptoms of ocular dryness. *Investig. Ophthalmol. Vis. Sci.* **2008**, *49*, 2971–2976. [CrossRef] [PubMed]

43. Knowlton, W.M.; Palkar, R.; Lippoldt, E.K.; McCoy, D.D.; Baluch, F.; Chen, J.; McKemy, D.D. A sensory-labeled line for cold: TRPM8-expressing sensory neurons define the cellular basis for cold, cold pain, and cooling-mediated analgesia. *J. Neurosci.* **2013**, *33*, 2837–2848. [CrossRef] [PubMed]

44. Namer, B.; Seifert, F.; Handwerker, H.O.; Maihöfner, C. TRPA1 and TRPM8 activation in humans: Effects of cinnamaldehyde and menthol. *Neuroreport* **2005**, *16*, 955–959. [CrossRef] [PubMed]

45. Proudfoot, C.J.; Garry, E.M.; Cottrell, D.F.; Rosie, R.; Anderson, H.; Robertson, D.C.; Fleetwood-Walker, S.M.; Mitchell, R. Analgesia mediated by the TRPM8 cold receptor in chronic neuropathic pain. *Curr. Biol.* **2006**, *16*, 1591–1605. [CrossRef] [PubMed]

46. Bouvier, V.; Roudaut, Y.; Osorio, N.; Aimonetti, J.-M.; Ribot-Ciscar, E.; Penalba, V.; Merrot, T.; Lebonvallet, N.; Le Gall-Ianotto, C.; Misery, L.; et al. Merkel Cells Sense Cooling with TRPM8 Channels. *J. Investig. Dermatol.* **2018**, *138*, 946–956. [CrossRef]

47. Urata, T.; Mori, N.; Fukuwatari, T. Vagus nerve is involved in the changes in body temperature induced by intragastric administration of 1,8-cineole via TRPM8 in mice. *Neurosci. Lett.* **2017**, *650*, 65–71. [CrossRef]

48. Almeida, M.C.; Hew-Butler, T.; Soriano, R.N.; Rao, S.; Wang, W.; Wang, J.; Tamayo, N.; Oliveira, D.L.; Nucci, T.B.; Aryal, P.; et al. Pharmacological Blockade of the Cold Receptor TRPM8 Attenuates Autonomic and Behavioral Cold Defenses and Decreases Deep Body Temperature. *J. Neurosci.* **2012**, *32*, 2086–2099. [CrossRef]

49. Hatta, A.; Kurose, M.; Sullivan, C.; Okamoto, K.; Fujii, N.; Yamamura, K.; Meng, I.D. Dry eye sensitizes cool cells to capsaicin-induced changes in activity via TRPV1. *J. Neurophysiol.* **2019**, *121*, 2191–2201. [CrossRef]
50. Bereiter, D.A.; Rahman, M.; Thompson, R.; Stephenson, P.; Saito, H. TRPV1 and TRPM8 Channels and Nocifensive Behavior in a Rat Model for Dry Eye. *Investig. Ophthalmol. Vis. Sci.* **2018**, *59*, 3739–3746. [CrossRef]

51. Rahman, M.; Okamoto, K.; Thompson, R.; Katagiri, A.; Bereiter, D.A. Sensitization of trigeminal brainstem pathways in a model for tear deficient dry eye. *Pain* **2015**, *156*, 942–950. [CrossRef]

52. Belmonte, C.; Aracil, A.; Acosta, M.C.; Luna, C.; Gallar, J. Nerves and sensations from the eye surface. *Ocul Surf.* **2004**, *2*, 248–253. [CrossRef]

53. González-González, O.; Bech, F.; Gallar, J.; Merayo-Lloves, J.; Belmonte, C. Functional Properties of Sensory Nerve Terminals of the Mouse Cornea. *Investig. Ophthalmol. Vis. Sci.* **2017**, *58*, 404–415. [CrossRef]

54. Brock, J.A.; Pianova, S.; Belmonte, C. Differences between nerve terminal impulses of polymodal nociceptors and cold sensory receptors of the guinea-pig cornea. *J. Physiol.* **2001**, *533*, 493–501. [CrossRef] [PubMed]

55. Bech, F.; González-González, O.; Artime, E.; Serrano, J.; Alcalde, I.; Gallar, J.; Merayo-Lloves, J.; Belmonte, C. Functional and Morphologic Alterations in Mechanical, Polymodal, and Cold Sensory Nerve Fibers of the Cornea Following Photorefractive Keratectomy. *Investig. Ophthalmol. Vis. Sci.* **2018**, *59*, 2281–2292. [CrossRef]

56. Baudouin, C.; Messmer, E.M.; Aragona, P.; Geerling, G.; Akova, Y.A.; Benítez-del-Castillo, J.; Boboridis, K.G.; Merayo-Lloves, J.; Rolando, M.; Labetoulle, M. Revisiting the vicious circle of dry eye disease: A focus on the pathophysiology of meibomian gland dysfunction. *Br. J. Ophthalmol.* **2016**, *100*, 300–306. [CrossRef]

57. Liu, H.; Hua, L.; Liu, Q.; Pan, J.; Bao, Y. Cold Stimuli Facilitate Inflammatory Responses Through Transient Receptor Potential Melastatin 8 (TRPM8) in Primary Airway Epithelial Cells of Asthmatic Mice. *Inflammation* **2018**, *41*, 1266–1275. [CrossRef]

58. Liu, H.; Liu, Q.; Hua, L.; Pan, J. Inhibition of transient receptor potential melastatin 8 alleviates airway inflammation and remodeling in a murine model of asthma with cold air stimulus. *Acta Biochim. Biophys. Sin.* **2018**, *50*, 499–506. [CrossRef]

59. Huang, J.-F.; Zhang, Y.; Rittenhouse, K.D.; Pickering, E.H.; McDowell, M.T. Evaluations of Tear Protein Markers in Dry Eye Disease: Repeatability of Measurement and Correlation with Disease. *Investig. Ophthalmol. Vis. Sci.* **2012**, *53*, 4556–4564. [CrossRef]

60. Willcox, M.D.P.; Argüeso, P.; Georgiev, G.A.; Holopainen, J.M.; Laurie, G.W.; Millar, T.J.; Papas, E.B.; Rolland, J.P.; Schmidt, T.A.; Stahl, U.; et al. TFOS DEWS II Tear Film Report. *Ocul. Surf.* **2017**, *15*, 366–403. [CrossRef]

61. Nicolle, P.; Liang, H.; Reboussin, E.; Rabut, G.; Warcoin, E.; Brignole-Baudouin, F.; Melik-Parsadaniantz, S.; Baudouin, C.; Labbe, A.; Reaux-Le Goazigo, A. Proinflammatory Markers, Chemokines, and Enkephalin in Patients Suffering from Dry Eye Disease. *Int. J. Mol. Sci.* **2018**, *19*, 1221. [CrossRef]

62. Noël, J.; Zimmermann, K.; Busserolles, J.; Deval, E.; Alloui, A.; Diochot, S.; Guy, N.; Borsotto, M.; Reeh, P.; Eschalier, A.; et al. The mechano-activated K+ channels TRAAK and TREK-1 control both warm and cold perception. *EMBO J.* **2009**, *28*, 1308–1318. [CrossRef] [PubMed]

63. Madrid, R.; de la Peña, E.; Donovan-Rodriguez, T.; Belmonte, C.; Viana, F. Variable threshold of trigeminal cold-thermosensitive neurons is determined by a balance between TRPM8 and Kv1 potassium channels. *J. Neurosci.* *Off. J. Soc. Neurosci.* **2009**, *29*, 3120–3131. [CrossRef] [PubMed]

64. Abrahamsen, B.; Zhao, J.; Asante, C.O.; Cendan, C.M.; Marsh, S.; Martinez-Barbera, J.P.; Nassar, M.A.; Dickenson, A.H.; Wood, J.N. The cell and molecular basis of mechanical, cold, and inflammatory pain. *Science* **2008**, *321*, 702–705. [CrossRef] [PubMed]

65. Li, F.; Yang, W.; Jiang, H.; Guo, C.; Huang, A.J.W.; Hu, H.; Liu, Q. TRPV1 activity and substance P release are required for corneal cold nociception. *Nat. Commun.* **2019**, *10*, 5678. [CrossRef]

66. Cao, Y.Q.; Mantyh, P.W.; Carlson, E.J.; Gillespie, A.M.; Epstein, C.J.; Basbaum, A.I. Primary afferent tachykinins are required to experience moderate to intense pain. *Nature* **1998**, *392*, 390–394. [CrossRef]

67. Amann, R.; Schuligoi, R.; Holzer, P.; Donnerer, J. The non-peptide NK1 receptor antagonist SR140333 produces long-lasting inhibition of neurogenic inflammation, but does not influence acute chemo- or thermonociception in rats. *Naunyn Schmiedebergs Arch. Pharmacol.* **1995**, *352*, 201–205. [CrossRef]

68. Klein, A.H.; Sawyer, C.M.; Carstens, M.I.; Tsagareli, M.G.; Tsiklauri, N.; Carstens, E. Topical application of L-menthol induces heat analgesia, mechanical allodynia, and a biphasic effect on cold sensitivity in rats. *Behav. Brain Res.* **2010**, *212*, 179–186. [CrossRef]
69. Langford, D.J.; Bailey, A.L.; Chanda, M.L.; Clarke, S.E.; Drummond, T.E.; Echols, S.; Glick, S.; Ingrao, J.; Klassen-Ross, T.; Lacroix-Fralish, M.L.; et al. Coding of facial expressions of pain in the laboratory mouse. *Nat. Methods* **2010**, *7*, 447–449. [CrossRef]

70. Matsumiya, L.C.; Sorge, R.E.; Sotocinal, S.G.; Tabaka, J.M.; Wieskopf, J.S.; Zaloum, A.; King, O.D.; Mogil, J.S. Using the Mouse Grimace Scale to reevaluate the efficacy of postoperative analgesics in laboratory mice. *J. Am. Assoc. Lab. Anim. Sci.* **2012**, *51*, 42–49.

71. Pauly, A.; Brignole-Baudouin, F.; Labbé, A.; Liang, H.; Warnet, J.-M.; Baudouin, C. New Tools for the Evaluation of Toxic Ocular Surface Changes in the Rat. *Investig. Ophthalmol. Vis. Sci.* **2007**, *48*, 5473–5483. [CrossRef]

72. Joubert, F.; Acosta, M.D.C.; Gallar, J.; Fakih, D.; Sahel, J.-A.; Baudouin, C.; Bodineau, L.; Mélik Parsadaniantz, S.; Réaux-Le Goazigo, A. Effects of corneal injury on ciliary nerve fibre activity and corneal nociception in mice: A behavioural and electrophysiological study. *Eur. J. Pain* **2019**, *23*, 589–602. [CrossRef] [PubMed]

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_12.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).