# Neurovascular microcirculatory vasodilation mediated by C-fibers and Transient receptor potential vanilloid-type-1 channels (TRPV 1) is impaired in type 1 diabetes

Pauline Marche, Séverine Dubois, Pierre Abraham, Elsa Parot-Schinkel, Lydie Gascoin, Anne Humeau-Heurtier, Pierre-Henry Ducluzeau-Fieloux, Guillaume Mahé

► To cite this version:

Pauline Marche, Séverine Dubois, Pierre Abraham, Elsa Parot-Schinkel, Lydie Gascoin, et al. Neurovascular microcirculatory vasodilation mediated by C-fibers and Transient receptor potential vanilloid-type-1 channels (TRPV 1) is impaired in type 1 diabetes. Scientific Reports, 2017, 7 (1), pp.44322. 10.1038/srep44322. hal-01507629

HAL Id: hal-01507629

https://univ-rennes.hal.science/hal-01507629

Submitted on 27 Jun 2018

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>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

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

OPEN

# Neurovascular microcirculatory vasodilation mediated by C-fibers and Transient receptor potential vanilloid-type-1 channels (TRPV 1) is impaired in type 1 diabetes

P. Marche<fnref n="1" />, S. Dubois<fnref n="1" />, P. Abraham², E. Parot-Schinkel³, L. Gascoin², A. Humeau-Heurtier⁴, PH. Ducluzeau⁵ & G. Mahe⁶

Microvascular dysfunction may have an early onset in type 1 diabetes (T1D) and can precede major complications. Our objectives were to assess the endothelial-dependent (acetylcholine, ACh; and post-occlusive hyperemia, PORH), non-endothelial-dependent (sodium nitroprusside, SNP) and neurovascular-dependent (local heating, LH and current induced vasodilation, CIV) microcirculatory vasodilation in T1D patients compared with matched control subjects using a laser speckle contrast imager. Seventeen T1D patients - matched with 17 subjects according to age, gender, Body-Mass-Index, and smoking status - underwent macro- and microvascular investigations. The LH early peak assessed the transient receptor potential vanilloid type 1 channels (TRPV1) mediated vasodilation, whereas the plateau assessed the Nitirc-Oxyde (NO) and endothelium-derived hyperpolarizing factor (EDHF) pathways. PORH explored sensory nerves and (EDHF), while CIV assessed sensory nerves (C-fibers) and prostaglandin-mediated vasodilation. Using neurological investigations, we observed that C-fiber and A-delta fiber functions in T1D patients were similar to control subjects. PORH, CIV, LH peak and plateau vasodilations were significantly decreased in T1D patients compared to controls, whereas there was no difference between the two groups for ACh and SNP vasodilations. Neurovascular microcirculatory vasodilations (C-fibers and TRPV 1-mediated vasodilations) are impaired in TD1 patients whereas no abnormalities were found using clinical neurological investigations. Clinicaltrials: No. NCT02538120.

Over 200 million people have diabetes worldwide, 10% of whom suffer from Type 1 diabetes (T1D). This disease preferentially affects young adults, and can lead to several complications such as sensorimotor neuropathy and macroangiopathy<fnref n="1" />. Diabetic neuropathy is one of the most common chronic complications of diabetes, affecting more than 50% at the advanced stage of the disease (more than 20 years of evolution)². This sensorimotor polyneuropathy can lead to diabetic foot ulcer, and in some cases, to Charcot disease (neuroarthropathy of the foot). Unfortunately, in the natural history of the disease, the neuropathy remains asymptomatic over a long period, inducing delayed diagnosis³. The rate of lower limb amputation is however 15 times higher in diabetic patients than in nondiabetic patients, and more than 50% of diabetic amputees require subsequent amputation of the contralateral limb within 4 years of loss of the first leg⁴. Moreover, macroangiopathy related to atherosclerosis is a major cause of morbidity and mortality in diabetic patients⁵. Endothelial dysfunction, while discussed in terms of vascular beds (conductance artery vs resistance artery) and the stage of the disease⁶,⁷, may appear early

Received: 14 November 2016

Accepted: 06 February 2017

Published: 13 March 2017

<footnote n="1">¹Endocrinology Department, University Hospital of Angers, Angers, France. ²Laboratory of Vascular Investigations, University Hospital of Angers, UMR CNRS6214/INSERM1083, LUNAM Université, Angers, France. ³Methodology and Biostatistics Unit, University Hospital of Angers, Angers, France. ⁴University of Angers, LARIS - Laboratoire Angevin de Recherche en Ingénierie des Systèmes, Angers, France. ⁵Endocrinology Department, University Hospital of Tours, Tours, France. ⁶INSERM Clinical Investigation Center (CIC 14 14), Rennes, France, Université de Rennes 1 and LUNAM University, Inserm 1083/CNRS 6214, Faculty of Medicine, Angers, France. Correspondence and requests for materials should be addressed to G.M. (email: maheguillaume@yahoo.fr)</footnote>
<table><caption>**Table 1.** Baseline clinical characteristics of type 1 diabetic patients and control subjects. Data presented as mean ± SD, median (interquartile range), or n (%). BMI: body mass index; SBP: systolic blood pressure. Bpm means beat per minute.</caption><thead><tr><th>Patient characteristics</th><th>Control subjects (n = 17)</th><th>Type 1 diabetic patients (n = 17)</th><th>p</th></tr></thead><tbody><tr><td>Clinical variables</td><td></td><td></td><td></td></tr><tr><td>Male n (%)</td><td>8 (47)</td><td>7 (41)</td><td>0.73</td></tr><tr><td>Smoker n (%)</td><td>6 (35)</td><td>7 (41)</td><td>0.72</td></tr><tr><td>Age (years)</td><td>24.6 (22.9–24.6)</td><td>27.4 (20.0–36.8)</td><td>0.39</td></tr><tr><td>Duration of diabetes (years)</td><td></td><td>9.29 ± 7.67</td><td></td></tr><tr><td>BMI (kg/m²)</td><td>21.67 ± 2.64</td><td>21.6 ± 2.28</td><td>0.94</td></tr><tr><td>SBP (mmHg)</td><td>113.53 ± 10.38</td><td>118.24 ± 10.66</td><td>0.20</td></tr><tr><td>Heart rate (bpm)</td><td>68.59 ± 10.04</td><td>70.82 ± 9.73</td><td>0.51</td></tr><tr><td>Pulse wave velocity (m/s)</td><td>4.95 (4.8–5.2)</td><td>4.98 (4.7–5.7)</td><td>0.69</td></tr><tr><td>Right ankle brachial index</td><td>1.17</td><td>1.16</td><td>0.70</td></tr><tr><td>Left ankle brachial index</td><td>1.15</td><td>1.16</td><td>0.92</td></tr></tbody></table>

on in the disease process and contributes to the induction and development of atherosclerotic plaques. Indeed, the endothelium plays a major role in blood flow regulation by producing several substances, such as nitric oxide (NO), prostacyclin, and endothelium-derived hyperpolarizing factor (EDHF), and participates in the control of many vascular functions: modulation of vasomotor tone, coagulation, cell adhesion mechanisms and migration, inflammatory processes, remodeling mechanisms, and vascular permeability⁸. Detecting all these complications as early as possible to improve patient management thus represents a major challenge for clinicians.

It has long been thought that microcirculation abnormalities in diabetes precede major complications⁹. Cutaneous microcirculation is considered as a mirror of generalized vascular function and is a useful translational model for investigating mechanisms of cutaneous physiology and cutaneous pathophysiology due either to skin disease or to other diseases such as diabetes¹⁰. Various microvascular tests for some of these, though not validated, have been used to assess both neurovascular abnormalities and endothelial dysfunction in T1D patients⁴,⁶,⁷,¹¹⁻¹⁸. The heterogeneity of both tests used and groups of patients studied does not allow us to understand the changes to the microvascular mechanisms involved in T1D or the time course over which microvascular abnormalities develop. And yet several dedicated tests are available and have been well-described for assessing the various mechanisms of vasodilation such as the NO-mediated pathway, the EDHF-mediated pathway and neurovascular pathways, whether mediated or not by transient receptor potential vanilloid Type 1 channels (TRPV-1)¹⁹. This latter pathway is of particular interest since it has been suggested that TRPV-1 may play a major role in the occurrence of T1D²⁰,²¹ and in cardiovascular mortality²². We hypothesized that microvascular abnormalities, especially impaired TRPV-1 mediated vasodilation, can be found in type 1 diabetes.

The aims of the present study were therefore to use five dedicated microvascular tests to assess neurovascular functions and endothelial functions in T1D patients compared to matched control subjects, and to determine the time course of the microvascular abnormalities.

## Results

Seventeen young adults with T1D (age 27.4 (20.0–36.8), 41% male, duration of diabetes 9.3 ± 7.8 years) and 17 controls matched for age, BMI, smoking status, and sex (age 24.6 (22.9–24.6), 47% male) were studied. All the diabetic patients were treated with insulin (insulin pens or subcutaneous insulin pump). We did not note any significant difference between the diabetic and control groups for blood pressure, pulse wave, and lipid levels (see Tables 1 and 2). Furthermore, as expected, fasting glycemia (8.3 (6.6–11.7) mmol/l vs. 4.7 (4.5–4.9) mmol/l, *p* = 0.0006) and glycosylated hemoglobin (8.4 (7.7–9.5)% vs 5.3 (5.2–5.5)%, *p* < 0.0001) were higher in the diabetic group than in the control group. Moreover, skin temperature at baseline and during each experiment was not significantly different between control and T1D subjects.

**Neurological investigations.** Our results show no difference on the upper limb and the lower limb for any of the parameters studied between T1D and controls, except for the diapason at the hallux and the vibration perception at the foot and hand: for them, the thresholds were significantly higher for T1D (see Table 3).

**Microvascular investigations.** We found similar baseline vascular conductance in control subjects and T1D patients, whatever the microvascular investigations. Moreover, no statistical difference was observed between the two groups for the peak value observed after ACh and for the peak or plateau observed after SNP iontophoresis (see Table 4).

Furthermore, the peak obtained during PORH was significantly lower in the T1D group than in the control subjects (1.01 ± 0.21 LSPU/mmHg and 1.26 ± 0.23 LSPU/mmHg, respectively, *p* = 0.003). Similarly, the peak (within 7 minutes following the start of heating) and the late plateau after heating (22 minutes following the start of recording) were both significantly decreased in the diabetic group compared to controls (0.81 ± 0.21 LSPU/mmHg and 1.07 ± 0.18 LSPU/mmHg for the peak; 0.81 ± 0.29 LSPU/mmHg and 1.13 ± 0.23 LSPU/mmHg for the plateau, respectively, *p* = 0.0004 and *p* = 0.002) (Fig. 1).
<table><caption>**Table 2.** Baseline biological characteristics of type 1 diabetic patients and control subjects. Data presented as mean ± SD or median (interquartile range). APTT: activated partial thromboplastin time; PR: Prothrombin ratio; MPV: mean platelet volume; Hb: Hemoglobin; Ht: Hematocrit. *Chronic Kidney Disease Epidemiology Collaboration.</caption><thead><tr><th>Patient characteristics</th><th>Control subjects (n = 17)</th><th>Type 1 diabetic patients (n = 17)</th><th>p</th></tr></thead><tbody><tr><td>Biological variables</td><td></td><td></td><td></td></tr><tr><td>APTT (s)</td><td>32 (31–33)</td><td>30 (28–31)</td><td>0.07</td></tr><tr><td>PR (%)</td><td>94.12 ± 7.34</td><td>93.82 ± 10.16</td><td>0.92</td></tr><tr><td>Platelets (G/L)</td><td>250.4 ± 68.49</td><td>276.5 ± 40.97</td><td>0.18</td></tr><tr><td>MPV (fl)</td><td>10.61 ± 0.57</td><td>10.28 ± 0.78</td><td>0.16</td></tr><tr><td>Hb (g/dL)</td><td>14.30 ± 1.23</td><td>14.02 ± 1.26</td><td>0.52</td></tr><tr><td>Ht (%)</td><td>43.59 ± 2.90</td><td>41.61 ± 3.31</td><td>0.07</td></tr><tr><td>White cells (G/L)</td><td>5.42 (5.20–6.76)</td><td>5.85 (5.25–7.05)</td><td>0.73</td></tr><tr><td>Creatinin (µmol/L)</td><td>77.18 ± 11.99</td><td>67.68 ± 8.54</td><td>0.01</td></tr><tr><td>CKD-EPI* (mL/mn)</td><td>107.5 ± 11.1</td><td>114.0 ± 16.1</td><td>0.18</td></tr><tr><td>LDL (mmol/L)</td><td>2.74 ± 0.50</td><td>2.91 ± 0.69</td><td>0.41</td></tr><tr><td>HDL (mmol/L)</td><td>1.31 ± 0.32</td><td>1.38 ± 0.25</td><td>0.77</td></tr><tr><td>Triglycerids (mmol/L)</td><td>0.79 (0.65–0.87)</td><td>0.77 (0.56–0.96)</td><td>0.88</td></tr><tr><td>Glucose (mmol/L)</td><td>4.7 (4.5–4.9)</td><td>8.3 (6.6–11.7)</td><td>0.0006</td></tr><tr><td>HbA1C (%)</td><td>5.3 (5.2–5.5)</td><td>8.4 (7.7–9.5)</td><td>&lt;0.0001</td></tr><tr><td>HbA1C (mmol/mol)</td><td>34.3 ± 2.5</td><td>71.1 ± 21.7</td><td>&lt;0.0001</td></tr><tr><td>Microalbuminuria (mg/L)</td><td>11 (10–17)</td><td>10 (6–15)</td><td>0.25</td></tr></tbody></table>

<table><caption>**Table 3.** Neurological characteristics of the populations. NSS is neuropathy symptom score, NDS is neuropathy disability score. AU means Arbitrary Unit and (Δ °C) means variation from 32 °Celsius.</caption><thead><tr><th></th><th>Control subjects</th><th>Type 1 diabetic patients</th><th>p</th></tr></thead><tbody><tr><td>Neuropathy scores</td><td></td><td></td><td></td></tr><tr><td>NSS</td><td>0 (0–0)</td><td>0 (0–4)</td><td>0.30</td></tr><tr><td>NDS hand</td><td>0 (0–0)</td><td>0 (0–0)</td><td>1</td></tr><tr><td>NDS foot</td><td>0 (0–0)</td><td>0 (0–1)</td><td>0.08</td></tr><tr><td>Sensory thresholds</td><td></td><td></td><td></td></tr><tr><td>Vibration hand (AU)</td><td>0.92 (0.55–1.24)</td><td>1.29 (1.03–2.74)</td><td><strong>0.049</strong></td></tr><tr><td>Vibration foot (AU)</td><td>1.45 (1.04–1.67)</td><td>3.04 (1.79–3.93)</td><td><strong>0.003</strong></td></tr><tr><td>Cold detection hand (Δ °C)</td><td>1.2 (0.8–1.4)</td><td>1.1 (1.0–1.6)</td><td>0.94</td></tr><tr><td>Cold detection foot (Δ °C)</td><td>2.5 (2–3.2)</td><td>3.4 (1.8–4.8)</td><td>0.23</td></tr><tr><td>Cold pain hand (Δ °C)</td><td>15.3 (7.6–26.9)</td><td>17.3 (9.1–22.5)</td><td>0.91</td></tr><tr><td>Cold pain foot (Δ °C)</td><td>17.4 (8.0–23.8)</td><td>18.2 (8.2–22.0)</td><td>0.51</td></tr><tr><td>Heat detection hand (Δ °C)</td><td>2.0 (1.2–2.7)</td><td>2.2 (1.4–2.8)</td><td>0.78</td></tr><tr><td>Heat detection foot (Δ °C)</td><td>4.3 (3.2–4.9)</td><td>4.0 (2.9–6.2)</td><td>0.89</td></tr><tr><td>Heat pain hand (Δ °C)</td><td>11.3 (8.6–13.5)</td><td>11.7 (10.5–13.6)</td><td>0.72</td></tr><tr><td>Heat pain foot (Δ °C)</td><td>11.6 (10.1–14.6)</td><td>12.7 (8.3–14.4)</td><td>0.74</td></tr><tr><td>Diapason hand (AU)</td><td>8 (8–8)</td><td>8 (7–8)</td><td>0.06</td></tr><tr><td>Diapason toe (AU)</td><td>8 (8–8)</td><td>7.5 (7–8)</td><td><strong>0.01</strong></td></tr></tbody></table>

For the CIV test, we did not note any difference in vascular conductance at rest between the T1D patients and the control subjects (see Table 4). However, the mean plateau value observed after the second 5 s cathodal current application was significantly lower in the T1D patients than in the matched controls (0.73 ± 0.33 LSPU/mmHg and 0.96 ± 0.27 LSPU/mmHg, respectively, *p* = 0.03).

Discussion

Our study, which compared five specific microvascular responses at the forearm level between T1D patients and control subjects matched for age, sex, BMI and smoking status, shows that neurovascular responses are significantly impaired in T1D patients with normal microalbuminuria and without macrovascular abnormalities compared with control subjects, whereas ACh-response and non-endothelial-dependent (ie, SNP) response are preserved.
<table><caption>**Table 4.** Microvascular characteristics of the population. Data presented as mean ± SD. Ach is Acetylcholine, SNP is Sodium Nitroprusside, PORH is post occlusive reactive hyperemia, CIV is current induced vasodilatation and T° cut is cutaneous temperature. LSPU means Laser Speckle Perfusion Unit.</caption><thead><tr><th>Microvascular investigations</th><th>Control subjects</th><th>Type 1 diabetic patients</th><th>p</th></tr></thead><tbody><tr><td><strong>ACh</strong></td><td></td><td></td><td></td></tr><tr><td>Base (LSPU/mmHg)</td><td>0.34 ± 0.09</td><td>0.34 ± 0.09</td><td>0.88</td></tr><tr><td>Peak (LSPU/mmHg)</td><td>0.78 ± 0.23</td><td>0.72 ± 0.25</td><td>0.49</td></tr><tr><td>Plateau (LSPU/mmHg)</td><td>0.49 ± 0.23</td><td>0.45 ± 0.19</td><td>0.73</td></tr><tr><td><strong>SNP</strong></td><td></td><td></td><td></td></tr><tr><td>Base (LSPU/mmHg)</td><td>0.38 ± 0.07</td><td>0.34 ± 0.09</td><td>0.16</td></tr><tr><td>Peak (LSPU/mmHg)</td><td>0.88 ± 0.19</td><td>0.83 ± 0.27</td><td>0.51</td></tr><tr><td>Plateau (LSPU/mmHg)</td><td>0.64 ± 0.15</td><td>0.63 ± 0.20</td><td>0.75</td></tr><tr><td><strong>Heat</strong></td><td></td><td></td><td></td></tr><tr><td>Base (LSPU/mmHg)</td><td>0.36 ± 0.11</td><td>0.33 ± 0.12</td><td>0.48</td></tr><tr><td>Early peak (LSPU/mmHg)</td><td>1.07 ± 0.18</td><td>0.81 ± 0.21</td><td><strong>0.0004</strong></td></tr><tr><td>Late plateau (LSPU/mmHg)</td><td>1.13 ± 0.23</td><td>0.81 ± 0.29</td><td><strong>0.0018</strong></td></tr><tr><td><strong>PORH</strong></td><td></td><td></td><td></td></tr><tr><td>Base (LSPU/mmHg)</td><td>0.44 ± 0.11</td><td>0.44 ± 0.15</td><td>0.91</td></tr><tr><td>Peak (LSPU/mmHg)</td><td>1.26 ± 0.23</td><td>1.01 ± 0.21</td><td><strong>0.0026</strong></td></tr><tr><td><strong>CIV</strong></td><td></td><td></td><td></td></tr><tr><td>Base (LSPU/mmHg)</td><td>0.34 ± 0.09</td><td>0.32 ± 0.08</td><td>0.53</td></tr><tr><td>Plateau (LSPU/mmHg)</td><td>0.96 ± 0.27</td><td>0.73 ± 0.33</td><td><strong>0.03</strong></td></tr><tr><td><strong>T° cut (°C)</strong></td><td>33.84 ± 1.44</td><td>34.26 ± 1.47</td><td>0.4</td></tr></tbody></table>

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Mean cutaneous vascular conductance (CVC: LSPU/mmHg) in type 1 diabetic patients in response to each microvascular stimulation compared to control subjects. Bars represent the mean ± SD. LSPU means Laser Speckle Perfusion Unit. \*p < 0.05.</figcaption></figure>

Diabetic polyneuropathy is linked to impaired large-diameter, myelinated (A-alpha and A-beta) fibers and small-diameter, unmyelinated nerve (A-delta and C) fibers¹⁶. Our clinical neurological results show that impairment, which was only found with A-beta fiber tests, is higher at the lower extremity level than at the upper extremity, confirming that nerve injury begins at the lower extremities, moving up to the upper extremities¹⁷. It has been suggested by previous authors that small-diameter, unmyelinated damage, especially C-fiber damage, appears earlier than large-diameter, myelinated (A-alpha and A-beta) fibers¹⁶. In our study using neurological investigations, C-fiber and A-delta fiber functions in T1D patients were similar to control subjects, whereas microvascular neurovascular responses were impaired in T1D patients compared with control subjects. This point is in accordance with Krishnan *et al.* whose work showed that C-fiber dysfunction assessed by “LDI flare” can be detected earlier than with clinical neurological investigations in type-2 diabetic (T2D) patients²³.

Our study shows that the axon reflex response assessed by two different neurovascular tests (peak local heating and cathodal CIV) is significantly impaired in T1D. As described elegantly by Minson *et al.* and Wong *et al.*, during the local thermal hyperemia test, the axon reflex mediates the early peak through transient receptor potential vanilloid type 1 (TRPV-1) channels located on sensory nerves, whereas the late plateau is dependent on NO.
(50-60%) and EDHFs (40-50%)¹⁹,²⁴. In our study, the significantly diminished peak response to local thermal hyperemia could suggest an abnormality (less activity) of the TRPV-1 channels at the skin level in T1D patients compared to control subjects. An important role of TRPV-1, which is distributed in the peripheral nervous system and in the network of sensory neurons that innervate the pancreatic islet, has already been evoked in T1D²⁰,²⁵. In the nerves of the islet beta cells, TRPV-1 can modulate T cell function. Inhibition of TRPV-1 has been shown to protect mice from autoimmune diabetes²¹, whereas TRPV-1 hypofunction initiates local inflammation and an autoimmune process in the pancreas leading to insulin resistance and beta cell destruction in non-obese diabetic mice²⁵. Razavi *et al.* have explained this controversy by a local feedback interaction between sensory nerves and beta cells²⁵. In 2013, Sadeh *et al.* also demonstrated that TRPV-1 might be a susceptible gene for T1D patients in a Jewish population²⁶, although this result needs to be confirmed in other populations. Interestingly, using TRPV-1 immunochemistry, Wilder-Smith *et al.* found increased keratinocyte TRPV-1 expression in diabetics compared to controls²⁷. We hypothesize that this overexpression could be linked to hypofunction of TRPV-1 that would induce a diminished microvascular peak response in our T1D patients, but this remains to be studied. In our study we also used the cathodal CIV, a second axon reflex test that is partially mediated by prostaglandins²⁸. This test is based on the anodal CIV test that has been widely used and studied in the literature²⁹,³⁰. Using this test, we confirm that T1D patients have an impairment of this response, suggesting either an impairment of the axon reflex response (C-fibers) or prostaglandin release as in T2D patients²⁹. As the late plateau phase of the ACh test is also partially mediated by prostaglandins³⁰ and was not different from the control subjects, we assume that the impairment of the CIV response in T1D patients is due to an impairment of the axon reflex response (C-fibers). This corroborates the results obtained with the local thermal peak suggesting a lesser activity of TRPV-1 channels. This lesser activity of TRPV-1 could also be a risk factor for cardiovascular diseases²².

Our study also demonstrates that the PORH peak is significantly reduced in T1D patients compared to con-
trol subjects. Several other authors have found the same results in T1D and T2D³¹, whereas others have found no
difference between diabetics and control subjects¹⁵,³² and, in one study, a surprising increase in PORH response³³.
The lack of difference between diabetics and control subjects could be due to a higher age of the population
studied³² or a higher duration of diabetes¹⁵ that may have minimized the difference between the populations. The
test used in our study enables us to assess the sensory nerves, the large-conductance calcium activated potassium
(BK<sub>CA</sub>) channels and EDHF, especially the 5,6-epoxyeicosatrienoic acid pathway³⁴,³⁵. This microvascular response
is mediated neither by NO nor by prostaglandins³⁵. The decrease of the PORH peak in T1D patients could suggest
a sensory nerve dysfunction, but also an EDHF dysfunction since the local thermal hyperemia plateau, that was
decreased in our T1D patients is also partially mediated by EDHF²⁴. Several authors have shown that the EDHF
pathway is altered in diabetes³⁶⁻³⁸. This decrease of the local thermal plateau response can be widely found in the
literature but the mechanism involved in such a decrease remains to be clarified¹²,¹⁴,³²,³⁹.

In our study, the ACh response, which is considered to be a specific test for endothelial function¹⁰, is similar between T1D patients and control subjects and is in accordance with previous publications¹⁴,¹⁷,⁴⁰,⁴¹. Contrary to our results, Khan *et al*.¹² showed a significantly reduced response to ACh in T1D. Two explanations can be given: i) skin blood flow measurements were not standardized for blood pressure whereas there was a difference between the systolic pressure of the diabetics compared to the controls; and ii) the subjects’ skin temperature was lower (30 °C) than in our study (33.8 °C) and this parameter is known to modify the microcirculatory response⁴²,⁴³. Gomes *et al*.¹⁵ also found a decreased response to ACh, although the difference was not significant for an electric charge of 2 mC, which is the one utilized in our study. Most of the other studies that have found a difference in ACh responses between diabetics and controls, used a higher electric charge (between 6 mC to 16 mC) or studied older populations⁶,¹²,¹³.

Furthermore, we did not find a significant difference in response to SNP that explores the endothelial-
independent pathway between T1D patients and controls. This NO donor acts directly on smooth muscle cells to
induce relaxation through an increase in cGMP formation⁴⁴. These findings with this electric charge (2 mC) are
consistent with previous studies⁶,¹²,¹⁵,¹⁸,⁴⁵ and demonstrate the integrity of the smooth muscle cell function. Other
authors found a decrease in SNP response, though using a higher electric charge or in older populations¹³,¹⁴.

Lastly, several studies were previously performed to evaluate microvascular abnormalities in T1D patients.
Most of these studies assessed one, two and up to four microvascular vasodilatory pathways at the same time,
limiting the knowledge of the time course for the onset of microvascular abnormalities⁷,¹¹,¹²,¹⁴,⁴⁶. According to
our results, neurovascular response impairment seems to appear earlier than ACh or SNP responses and before
impairment of macrovascular abnormalities. Therefore, for a clinician, the choice would have to focus either on
the local thermal test, the cathodal CIV test or the PORH test in order to assess microvascular dysfunction.

Our study has some limitations: The menstrual cycle stage was not checked in this study. This might have influenced the results since the role of the female hormonal status in skin reactivity is discussed⁴⁷. The proportion of females was not, however, statistically different between the groups. Moreover, while the absence of microangiopathy in T1D was confirmed by negative microalbuminuria, we had no information on diabetic retinopathy.

In conclusion, this study suggests that: i) neurovascular responses at the forearm level are impaired in T1D patients without impairment of the ACh and SNP responses compared with matched control-subjects; ii) microvascular response mediated by TRPV-1 channels is diminished; and iii) neurovascular abnormalities appears earlier than other microvascular abnormalities.

**Methods**

**Subjects.** Seventeen T1D patients were recruited in the diabetic care unit of the Angers University Hospital and matched with 17 healthy volunteers according to age, gender, Body Mass Index (BMI) and smoking status. Before their participation, all subjects were thoroughly informed of the methods and procedures and gave their written consent to participate in this prospective study approved by the “Comité Protection des Personnes-OUEST II” ethics committee (Number: 2012-A010008-35 and acceptance date: 2012). Informed consent was
obtained from all subjects. Experiments were carried out in accordance with the Declaration of Helsinki. The protocol was registered in the American National Institute of Health database under reference No. NCT02538120 (August, 2015).

Each subject had to meet the following inclusion criteria: be at least 18 years old, have read and signed the informed consent document. Diabetic subjects had to display at least one diagnostic criterion according to World Health Organization (Fasting plasma glucose ≥126 mg/l (7 mmol/l) twice after 8 h of fasting, or casual plasma glucose ≥200 mg/l (11.1 mmol/l) in presence of hyperglycemia symptoms (polyuria, polydipsia and unexplained weight loss), or 2-hour glucose ≥200 mg/l after an oral glucose tolerance test (OGTT) with 75 g glucose), and at least one positive antibody from among: ICA (against islet cytoplasm), GAD (against glutamic acid decarboxylase), IA2 (against tyrosine phosphatase) and IAA (anti-insulin). Healthy volunteers had glycated hemoglobin (HbA1c) equal to or lower than 5.6%. All treatments were allowed, except for anti-inflammatory drugs during the study and in the 10 days preceding inclusion, as they are known to abolish current-induced vasodilation (please refer to the microvascular tests section)²⁹.

**Protocol.** Each diabetic subject or volunteer was required to complete 2 distinct visits. The first one constituted the inclusion visit during which a medical examination (height, weight, BMI, blood pressure, tendon reflexes, Ankle Brachial Index [ABI]) and an interrogation searching for cardiovascular risk factors were performed. If the subject accepted to be included in the study, he/she was asked to sign the consent document and blood tests were conducted at the end of the first visit. The second visit consisted in neurological and microvascular investigations.

**Clinical neurological investigations.** The presence or absence of neuropathy was tested in accordance with the recommendations of the San Antonio Conference (1988) using the Neuropathy Symptom Score (NSS) and the Neuropathy Disability Score (NDS) evaluating large myelinated fibers. Both scores range from zero (absence of neuropathy) to ten⁴⁸. A technician ipsilateral to the vascular investigations performed neurological investigations. The quantitative evaluation was carried out by assessing the vibration perception threshold (VPT) with a neurothesiometer (Case-IV, Medical Electronics Co., USA) exploring A beta fibers. The VPT was also tested with a Rydel-Seiffer fork (64 Hz) exploring A beta fibers. The heat and cold detection thresholds (HDT and CDT) as well as the heat and cold pain thresholds (HPT and CPT) were tested (TSA II, Medical Electronics Co., USA), on the dorsum of the foot and hand examining C fibers and A delta fibers⁴⁹. For each measurement, the starting temperature of the probe was 32 degrees Celsius (°C) and the participants were asked to push a button when the test threshold was reached. For analysis, the relative heat deflection temperature, which is obtained by subtracting the skin temperature, was used.

**Macrovascular investigations.** Pulse wave velocity was assessed using a Mobil-O-Graph (I.E.M. GmbH, Stolberg, Germany) and the ankle-brachial index (ABI) using a hand-held Doppler on both limbs. The ABI was calculated as the ratio of the highest systolic pressure in the dorsalis pedis or posterior tibial artery to the highest systolic pressure between both arms according to guidelines⁵⁰.

**Microvascular investigations.** Experiments were performed with the participants resting in the supine position. All tests were performed in a quiet temperature-controlled room after a 20-minute acclimatization period.

Recently, laser speckle contrast imaging (LSCI) has been developed. LSCI uses the random speckle pattern generated by the illumination of the tissues under study by a coherent laser light. LSCI recordings of the cutaneous blood flow (CBF) in the forearm were performed using a 70-mW system (PeriCam PSI System®, Perimed, Järfälla, Sweden) with a laser wavelength of 785 nm⁴². The sampling frequency was 18 Hz. The distance between the laser head and skin surface was fixed at 15 cm taking into account environmental conditions⁴². The software expresses recorded values in laser speckle perfusion units (LSPU).

**Microvascular tests.** Five vasodilator microvascular tests involving different physiological vasodilation pathways were simultaneously performed on the same volar aspect of the forearm.

*ACh iontophoresis (endothelial-dependent pathway).* Iontophoresis is a non-invasive method that drives a pharmacologically-charged drug by electrorepulsion through the interstitium surrounding the blood vessels¹⁰. Transdermal iontophoresis of ACh (Sigma-Aldrich Corporation, L'Isle d'Abeau, France) was performed with an iontophoresis chamber (LI 611, Perimed Jarfalla, Sweden) and corresponded to the anode. The iontophoresis chamber was filled with ACh (2%) dissolved in deionized water. The peak, which is dependent on the muscarinic receptor, and the late plateau, which is partially prostaglandin-dependent, were studied.

*SNP iontophoresis (endothelial-independent pathway).* Transdermal iontophoresis of SNP (Nitriate®, SERB, Paris, France) was performed with an iontophoresis chamber (LI 611, Perimed Jarfalla, Sweden) and corresponded to the cathode. The iontophoresis chamber was filled with SNP (1%) dissolved in deionized water. The PeriIont Micropharmacology System (Perimed, Jarfalla, Sweden) was plugged into the iontophoresis chamber containing ACh and to the iontophoresis chamber containing SNP and delivered the current. The current stimulation was set to 0.1 mA for 20 seconds.
**PORH (endothelial-dependent and endothelial-independent pathway).** The PORH for exploring neurovascular interactions and probably the endothelium-derived hyperpolarizing factor (EDHF)⁵¹ was measured by producing ischemia for 3 minutes with a blood pressure cuff inflated to 50 mmHg over the subject's systolic blood pressure⁵¹.

**Current-induced vasodilation (CIV) (neurovascular pathway and prostaglandin dependent pathway).** Transdermal iontophoresis of deionized water was performed with an iontophoresis chamber (LI 611, Perimed Jarfalla, Sweden) and corresponded to the cathode. The current application was delivered on the active probe in two consecutive 5-second applications separated by 4 minutes²⁸. For CIV, we studied the late-plateau response.

**Local thermal hyperemia (neurovascular and endothelial-dependent pathway).** The hyperemic response to heat was evaluated by using a skin-heating probe filled with deionized water and heated to 44°C to assess heat-induced vasodilation. Heat-induced vasodilation is due to nerve fibers for the early peak and to the NO-dependent pathway and the EDHF-dependent pathway for the late plateau²⁴,⁵².

**Number of subjects.** Searching for a difference of at least 0.35 in CIV between the two groups with a SD of 0.30 for each group based on our previous study²⁹, the minimum number of subjects to be included for alpha = 0.05 and a 90% power was 16 per group.

**Data and statistical analysis.** The descriptive results are presented as the mean +/– standard deviation (SD) for variables with normal distribution and as the median (25 centile–75 centile) when distribution was not normal. Skin blood flow was measured over a region of interest (ROI) and a time of interest (TOI) from images recorded by the LSCI. The size of the ROI was 20 mm². The TOI was set to 5 s for the peak and 10 s for the baseline and plateau measurements⁵³. The results for the microvascular recordings (peak or plateau) are expressed as cutaneous vascular conductance (CVC). The CVC was calculated as the ratio of CBF to the mean arterial pressure (mmHg) and expressed in LSPU/mmHg¹⁰. It has previously been reported that, when using LSCI, expressing the results in CVC is the best way to obtain the lowest coefficient of variation (18.7% with an Intra-Class Correlation coefficient = 0.87 [0.65–0.95]) for the 7-day interval of reproducibility⁵⁴.

The normality of the distribution of variables was tested with the Shapiro-Wilk test. The microvascular results obtained in T1D patients and the control subjects were compared with a T-test when there was normal distribution and with a Mann-Whitney test when the distribution was not normal. Statistical analyses were performed using MedCalc for Windows, version 12.5 (MedCalc Software, Ostend, Belgium). For each statistical analysis, a two-tailed p value ≤ 0.05 was considered significant.

**References**

1. Harjutsalo, V., Forsblom, C. & Groop, P.-H. Time trends in mortality in patients with type 1 diabetes: nationwide population based cohort study. *BMJ* **343**, d5364 (2011).

2. Deli, G., Bosnyak, E., Pusch, G., Komoly, S. & Feher, G. Diabetic neuropathies: diagnosis and management. *Neuroendocrinology* **98**, 267–280 (2013).

3. Perkins, B. A., Olaleye, D., Zinman, B. & Bril, V. Simple screening tests for peripheral neuropathy in the diabetes clinic. *Diabetes Care* **24**, 250–256 (2001).

4. Veves, A. et al. Endothelial dysfunction and the expression of endothelial nitric oxide synthetase in diabetic neuropathy, vascular disease, and foot ulceration. *Diabetes* **47**, 457–463 (1998).

5. Joannides, R., Bellien, J. & Thuillez, C. Clinical methods for the evaluation of endothelial function– a focus on resistance arteries. *Fundam. Clin. Pharmacol.* **20**, 311–320 (2006).

6. Quattrini, C., Harris, N. D., Malik, R. A. & Tesfaye, S. Impaired skin microvascular reactivity in painful diabetic neuropathy. *Diabetes Care* **30**, 655–659 (2007).

7. Bellien, J. et al. Early stage detection of conduit artery endothelial dysfunction in patients with type 1 diabetes. *Diab. Vasc. Dis. Res.* **7**, 158–166 (2010).

8. Puissant, C. et al. [Endothelial function: role, assessment and limits]. *J. Mal. Vasc.* **39**, 47–56 (2014).

9. McMillan, D. E. Deterioration of the microcirculation in diabetes. *Diabetes* **24**, 944–957 (1975).

10. Mahé, G., Humeau-Heurtier, A., Durand, S., Leftheriotis, G. & Abraham, P. Assessment of skin microvascular function and dysfunction with laser speckle contrast imaging. *Circ. Cardiovasc. Imaging* **5**, 155–163 (2012).

11. Nabavi Nouri, M. et al. Diabetic neuropathy and axon reflex-mediated neurogenic vasodilatation in type 1 diabetes. *PloS One* **7**, e34807 (2012).

12. Khan, F., Elhadd, T. A., Greene, S. A. & Belch, J. J. Impaired skin microvascular function in children, adolescents, and young adults with type 1 diabetes. *Diabetes Care* **23**, 215–220 (2000).

13. Katz, A., Ekberg, K., Johansson, B. L. & Wahren, J. Diminished skin blood flow in Type I diabetes: evidence for non-endothelium-dependent dysfunction. *Clin. Sci. Lond. Engl.* **1979** **101**, 59–64 (2001).

14. Heimhalt-El Hamriti, M. et al. Impaired skin microcirculation in paediatric patients with type 1 diabetes mellitus. *Cardiovasc. Diabetol.* **12**, 115 (2013).

15. Gomes, M. B., Matheus, A. S. M. & Tibiriçá, E. Evaluation of microvascular endothelial function in patients with type 1 diabetes using laser-Doppler perfusion monitoring: which method to choose? *Microvasc. Res.* **76**, 132–133 (2008).

16. Breiner, A., Lovblom, L. E., Perkins, B. A. & Bril, V. Does the prevailing hypothesis that small-fiber dysfunction precedes large-fiber dysfunction apply to type 1 diabetic patients? *Diabetes Care* **37**, 1418–1424 (2014).

17. Arora, S. et al. Differences in foot and forearm skin microcirculation in diabetic patients with and without neuropathy. *Diabetes Care* **21**, 1339–1344 (1998).

18. Johnstone, M. T. et al. Impaired endothelium-dependent vasodilation in patients with insulin-dependent diabetes mellitus. *Circulation* **88**, 2510–2516 (1993).

19. Wong, B. J. & Fieger, S. M. Transient receptor potential vanilloid type 1 channels contribute to reflex cutaneous vasodilation in humans. *J. Appl. Physiol. Bethesda Md* **1985** **112**, 2037–2042 (2012).

20. Brito, R., Sheth, S., Mukherjea, D., Rybak, L. P. & Ramkumar, V. TRPV1: A Potential Drug Target for Treating Various Diseases. *Cells* **3**, 517–545 (2014).

21. Jancsó, G., Kiraly, E. & Jancsó-Gábor, A. Pharmacologically induced selective degeneration of chemosensitive primary sensory neurones. *Nature* **270**, 741–743 (1977).
22. Robbins, N., Koch, S. E. & Rubinstein, J. Targeting TRPV1 and TRPV2 for potential therapeutic interventions in cardiovascular disease. *Transl. Res. J. Lab. Clin. Med.* **161**, 469–476 (2013).

23. Krishnan, S. T. M., Quattrini, C., Jeziorska, M., Malik, R. A. & Rayman, G. Abnormal LDIflare but normal quantitative sensory testing and dermal nerve fiber density in patients with painful diabetic neuropathy. *Diabetes Care* **32**, 451–455 (2009).

24. Minson, C. T., Berry, L. T. & Joyner, M. J. Nitric oxide and neurally mediated regulation of skin blood flow during local heating. *J. Appl. Physiol. Bethesda Md* **1985** **91**, 1619–1626 (2001).

25. Razavi, R. et al. TRPV1+ sensory neurons control beta cell stress and islet inflammation in autoimmune diabetes. *Cell* **127**, 1123–1135 (2006).

26. Sadeh, M. et al. Association of the M3151 variant in the transient receptor potential vanilloid receptor-1 (TRPV1) gene with type 1 diabetes in an Ashkenazi Jewish population. *Isr. Med. Assoc. J. IMAJ* **15**, 477–480 (2013).

27. Wilder-Smith, E. P., Ong, W.-Y., Guo, Y. & Chow, A. W.-L. Epidermal transient receptor potential vanilloid 1 in idiopathic small nerve fibre disease, diabetic neuropathy and healthy human subjects. *Histopathology* **51**, 674–680 (2007).

28. Mahé, G. et al. Evidence for a vasomotor cyclo-oxygenase dependent mechanism of sensitization at the cutaneous level. *Br. J. Clin. Pharmacol.* **80**, 185–192 (2015).

29. Rousseau, P. et al. Axon-reflex cutaneous vasodilatation is impaired in type 2 diabetic patients receiving chronic low-dose aspirin. *Microvasc. Res.* **78**, 218–223 (2009).

30. Durand, S., Fromy, B., Bouyé, P., Saumet, J. L. & Abraham, P. Current-induced vasodilation during water iontophoresis (5 min, 0.10 mA) is delayed from current onset and involves aspirin sensitive mechanisms. *J. Vasc. Res.* **39**, 59–71 (2002).

31. Yamamoto-Suganuma, R. & Aso, Y. Relationship between post-occlusive forearm skin reactive hyperaemia and vascular disease in patients with Type 2 diabetes-a novel index for detecting micro- and macrovascular dysfunction using laser Doppler flowmetry. *Diabet. Med. J. Br. Diabet. Assoc.* **26**, 83–88 (2009).

32. Araszkiewicz, A. et al. In diabetic Charcot neuroarthropathy impaired microvascular function is related to long lasting metabolic control and low grade inflammatory process. *Microvasc. Res.* **101**, 143–147 (2015).

33. Schlager, O. et al. Microvascular autoregulation in children and adolescents with type 1 diabetes mellitus. *Diabetologia* **55**, 1633–1640 (2012).

34. Lorenzo, S. & Minson, C. T. Human cutaneous reactive hyperaemia: role of BKCa channels and sensory nerves. *J. Physiol.* **585**, 295–303 (2007).

35. Cracowski, J.-L. et al. Involvement of cytochrome epoxynase metabolites in cutaneous postocclusive hyperemia in humans. *J. Appl. Physiol. Bethesda Md* **1985** **114**, 245–251 (2013).

36. Park, Y. et al. Role of EDHF in type 2 diabetes-induced endothelial dysfunction. *Am. J. Physiol. Heart Circ. Physiol.* **295**, H1982–1988 (2008).

37. Fitzgerald, S. M., Kemp-Harper, B. K., Tare, M. & Parkington, H. C. Role of endothelium-derived hyperpolarizing factor in endothelial dysfunction during diabetes. *Clin. Exp. Pharmacol. Physiol.* **32**, 482–487 (2005).

38. Coppey, L. J., Gellett, J. S. & Yorek, M. A. Mediation of vascular relaxation in epineurial arterioles of the sciatic nerve: effect of diabetes in type 1 and type 2 diabetic rat models. *Endothel. J. Endothel. Cell Res.* **10**, 89–94 (2003).

39. Wilson, S. B., Jennings, P. E. & Belch, J. J. Detection of microvascular impairment in type I diabetics by laser Doppler flowmetry. *Clin. Physiol. Oxf. Engl.* **12**, 195–208 (1992).

40. Kilo, S., Berghoff, M., Hilz, M. & Freeman, R. Neural and endothelial control of the microcirculation in diabetic peripheral neuropathy. *Neurology* **54**, 1246–1252 (2000).

41. Brooks, B. A., McLennan, S. V., Twigg, S. M. & Yue, D. K. Detection and characterisation of microcirculatory abnormalities in the skin of diabetic patients with microvascular complications. *Diab. Vasc. Dis. Res.* **5**, 30–35 (2008).

42. Mahé, G., Durand, S., Humeau-Heurtier, A., Leftheriotis, G. & Abraham, P. Impact of experimental conditions on noncontact laser recordings in microvascular studies. *Microcirc. N. Y. N* **1994** **19**, 669–675 (2012).

43. Abraham, P. et al. Effect of skin temperature on skin endothelial function assessment. *Microvasc. Res.* **88**, 56–60 (2013).

44. Morris, S. J. & Shore, A. C. Skin blood flow responses to the iontophoresis of acetylcholine and sodium nitroprusside in man: possible mechanisms. *J. Physiol.* **496** (Pt 2), 531–542 (1996).

45. Koitka, A. et al. Impaired pressure-induced vasodilation at the foot in young adults with type 1 diabetes. *Diabetes* **53**, 721–725 (2004).

46. Rossi, M. et al. Peripheral microvascular dysfunction as an independent predictor of atherosclerotic damage in type 1 diabetes patients: a preliminary study. *Clin. Hemorheol. Microcirc.* **54**, 381–391 (2013).

47. Ketel, I. J. G. et al. Microvascular function has no menstrual-cycle-dependent variation in healthy ovulatory women. *Microcirc. N. Y. N* **1994** **16**, 714–724 (2009).

48. Report and recommendations of the San Antonio conference on diabetic neuropathy. Consensus statement. *Diabetes* **37**, 1000–1004 (1988).

49. Greene, D. A., Sima, A. A., Stevens, M. J., Feldman, E. L. & Lattimer, S. A. Complications: neuropathy, pathogenetic considerations. *Diabetes Care* **15**, 1902–1925 (1992).

50. Aboyans, V. et al. Measurement and interpretation of the ankle-brachial index: a scientific statement from the American Heart Association. *Circulation* **126**, 2890–2909 (2012).

51. Yvonne-Tee, G. B., Rasool, A. H. G., Halim, A. S. & Rahman, A. R. Reproducibility of different laser Doppler fluximetry parameters of postocclusive reactive hyperemia in human forearm skin. *J. Pharmacol. Toxicol. Methods* **52**, 286–292 (2005).

52. Vinik, A. I. et al. Dermal neurovascular dysfunction in type 2 diabetes. *Diabetes Care* **24**, 1468–1475 (2001).

53. Rousseau, P. et al. Increasing the 'region of interest' and 'time of interest', both reduce the variability of blood flow measurements using laser speckle contrast imaging. *Microvasc. Res.* **82**, 88–91 (2011).

54. Puissant, C. et al. Reproducibility of non-invasive assessment of skin endothelial function using laser Doppler flowmetry and laser speckle contrast imaging. *PloS One* **8**, e61320 (2013).

## Acknowledgements

The authors would like to thank Professor P. Rodien from the University Hospital of Angers for his support and Y. Onillon from the University Hospital of Angers for her assistance in the inclusion of subjects. G. Mahé takes full responsibility for the work as a whole, including the study design, access to data, and the decision to submit and publish the manuscript. This study was sponsored by the "Centre Hospitalo-Universitaire d'Angers" as the NEURODIAB 1 study.

## Author Contributions

P.M. researched data, wrote the manuscript. S.D. wrote manuscript, researched data. P.A. reviewed/edited manuscript. E.P.S. reviewed/edited manuscript. L.G. researched data. A.H.H. reviewed/edited manuscript. P.H.D. reviewed/edited manuscript. GM. researched data, wrote manuscript.

## Additional Information

**Competing Interests:** The authors declare no competing financial interests.
**How to cite this article:** Marche, P. et al. Neurovascular microcirculatory vasodilation mediated by C-fibers and Transient receptor potential vanilloid-type-1 channels (TRPV 1) is impaired in type 1 diabetes. *Sci. Rep.* 7, 44322; doi: 10.1038/srep44322 (2017).

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

© The Author(s) 2017