Pathophysiological Mechanisms of the Onset, Development, and Disappearance Phases of Skin Eruptions in Chronic Spontaneous Urticaria ==================================================================================================================================== * Sungrim Seirin-Lee * Shunsuke Takahagi * Michihiro Hide ## Abstract Chronic spontaneous urticaria (CSU) is a typical example of an intractable skin disease with no clear cause and significantly affects daily life of patients. Because CSU is a human-specific disease and lacks proper animal model, there are many questions regarding its pathophysiological dynamics. On the other hand, most clinical symptoms of urticaria are notable as dynamic appearance of skin eruptions (wheals). In this study, we explored dynamics of wheal by dividing it into three phases using a mathematical model: onset, development, and disappearance. Our results suggest that CSU onset is critically associated with endovascular dynamics triggered by basophils positive feedback. In contrast, the development phase is regulated by mast cell dynamics via vascular gap formation. We also suggest a disappearance mechanism of skin eruptions in CSU through an extension of the mathematical model using qualitative and quantitative comparisons of wheal expansion data of real patients with urticaria. Our results suggest that the wheal dynamics of the three phases and CSU development are hierarchically related to endovascular and extravascular pathophysiological networks. Keywords * Chronic Spontaneous Urticaria * Skin eruption phases * CSU mathematical model ## 1 Introduction The skin is the most important organ in the human body for maintaining its homeostasis. It has physical and immunological barrier functions to protect the internal organs from the external environment. The skin is involved in the immune system and may develop life-threatening neoplasm, such as cancers and infectious diseases as well as various allergic and/or inflammatory diseases, such as atopic dermatitis, chronic urticaria and psoriasis that severely affect daily life of patients (Leiter et al. 2020; Seidel et al. 2018; Thandi and Whittam 2021; Weidinger et al. 2018; Kolkhir et al. 2022). Although most skin diseases are readily diagnosed by visual observation of skin eruptions, many of them lack appropriate animal models (Kolkhir et al. 2022). Even if animal models are available, many of them cannot precisely replicate correspodning human diseases, making it still difficult to elucidate the whole picture of underlying mechanisms *in vivo* and to develop effective treatments. Thus, it is challenging to fully understand the pathophysiology of skin diseases and find a way to cure by limited clinical data. Urticaria is one of the most common human skin disorders, affecting at least one in five people during their lifetime (Lee et al. 2017; Maxim et al. 2018). Chronic spontaneous urticaria (CSU), a major subtype of urticaria, is a typical example of an intractable skin disease with unknown cause (Greaves 2014). It is characterized by the appearance of notable skin eruptions called wheals, and is typically accompanied by distressing itching. Wheals of notable shapes recur on a daily basis in patients with CSU, and this significantly affects the quality of life (QoL) of patients’ daily lives over a long timescale, from months to decades (Itakura et al. 2018; Zuberbier et al. 2022). Wheal formation in urticaria is predominantly mediated by histamine which is released from skin mast cells and its action on the microvascular endothelium (Kolkhir et al. 2022; Seirin-Lee et al. 2023). In a certain population of patients with CSU, mast cells can be activated by cross-linkage of IgE against autoallergens, such as interleukin (IL)-24, double stand DNA and thyroid peroxidase, and/or IgG autoantibodies against IgE or the high affinity IgE receptor (Fc*ε*RI) (Maurer et al. 2021). A positive feedback reaction to histamine release from mast cells is suggested to play a role in the expansion of wheals (Bazilai et al. 2017). Antihistamines, the mainstay of treatments for CSU, may be effective in up to 70-80% of CSU patients (Maurer et al. 2011), but their efficacy varies among patients. This indicated that the pathomechanism of urticaria cannot be explained solely by histamine dynamics. The pathomechanism of wheal formation in CSU has been recently suggested to be involved by basophils in the peripheral blood circulation, blood coagulation, and complement systems (Yanase et al. 2021b; Huang et al. 2020; Rijavec et al. 2021). In our previous studies (Seirin-Lee et al. 2020, 2023), we successfully developed a mathematical model that integrates pathophysiological networks *in vivo* and recapitulated the geometric features of CSU skin eruptions *in silico*. Based on *in silico* experiments with changing parameters of this mathematical model, we obtained five types of eruption patterns with the development of the clinical criteria of CSU eruption geometry (EGe Criteria) and extracted critical pathophysiological networks for each type. Of note, eruptions appeared in almost of 105 patients with CSU were successfully classified into one of these five types based on EGe Criteria by six independent dermatologists. These observations suggest that morphological characteristics of wheals on the skin are predominantly, if not directly, and their related to dynamics of biological molecules and/or cells *in vivo*. Thus, the elucidation of these relationship by mathematical approach is a critical step to proceed forward to understand the overall mechanism of CSU and develop wheal shape-specific treatments. A daily or mostly daily dynamic phenomenon of the spontaneous appearance and disappearance is a notable characteristic of wheal formation in CSU (Zuberbier et al. 2022). The eruption shapes change, ranging in size from a few millimeters to several centimeters, and in time from several minutes to hours. Clinically, wheal dynamics can be divided into three phases; onset, development, and disappearance. The onset phase is critical as it is most likely linked to the triggering mechanism of the disease. Various patterns of eruptions are appeared during the developmental phase, suggesting the involvement of pathophysiological network in various ways during this phase (Seirin-Lee et al. 2023). However, the detailed mechanism of developing diversity in the patterns of eruptions in CSU remains mostly unexplored. The disappearance phase has received far more less attentions. However, the regulation of this phase could also open a way for therapies to prevent the progress and/or even to prevent the development of wheals. The mathematical model of previous study, Seirin-Lee et al. (2023), was developed based on pathophysiological network assumed to be involved in the onset and development phases, and the disappearance phase of wheal formation in CSU has not been well-explained yet. In the present study, we focus on more detailed dynamics of wheal formation in CSU and explore the involvement of pathophysiological dynamics in each phase of wheal formation. For this purpose, we choose the mathematical model developed by Seirin-Lee et al. (2023) and analyze it with reconstructing two independent models in the endovascular and extravascular sites. We also suggest a possible mechanism of eruption disappearance in CSU by extending a mathematical model with qualitative and quantitative comparisons of eruption expansion data of urticaria patients. Our results suggest that the dynamics of the three phases of wheal formation in CSU development are hierarchically related to endovascular and extravascular pathophysiological networks and that a global inhibitory mechanism of mast cells is essential in wheal dynamics of the disappearance phase. This study provides a better understanding of the mechanism of CSU by morphological dynamics in wheal formation and its relationship to cellular and molecular dynamics *in vivo*, thereby contributing to accurate diagnostic for individualized treatments decision making in clinical settings. ## 2 Method ### 2.1 Basal CSU model The pathophysiological structure that develops CSU *in vivo* is assumed as shown in Fig. 1A (Hide and Kaplan 2022; Seirin-Lee et al. 2023). Based on unexplained small stimuli, basophils that randomly adhere to the blood vessel endothelial tissue within capillaries begin to release histamine. The histamine released from basophils enhances the expression of tissue factor on the endothelial cell surface (*N*1) which activates the extrinsic coagulation pathway and histamine release from basophils (*N*2) (Karasuyama et al. 2009; Stone et al. 2010; Yanase et al. 2021a). This comprises a positive-feedback loop at the endovascular site (*N*1 − *N*2). At the same time, mast cells stimulated by another resource also start releasing histamine which enhances the expression of tissue factor (*N*1) and causes coagulation factors to leak from the blood vessels by creating a gap formation (*N*3) (Kurashima and Kiyono 2014; Yanase et al. 2018a,b). This further activates histamine release from the mast cells (*N*4) (Yanase et al. 2021a) and induces another positive-feedback loop at the extravascular site (*N*1 − *N*3 − *N*4). Networks *I*1, *I*2, and *I*3 are the networks inhibited by adenosine consequentially produced from adenosine triphosphate (ATP) which is released from basophils and mast cells simultaneously with histamine, and inhibits histamine release and tissue factor expression on endothelial cells (Rudich et al. 2012; Matsuo et al. 2018). ![Figure 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/11/12/2023.11.11.23298424/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2023/11/12/2023.11.11.23298424/F1) Figure 1: Summary of pathophysiological dynamics of CSU and network. (A) Schematic summary of pathophysiological CSU dynamics and details of endovascular and extravascular sites. (B) Network summary of pathophysiological dynamics in CSU. (A-B) The network *N*1−*N*2 consists of positive feedback in the endovascular site, and the network *N*1−*N*3−*N*4 consists of positive feedback in the extravascular site. *I*1, *I*2 and *I*3 are inhibition networks by adenosine. (C) Example of skin eruption patterns in a patient with CSU. (D) Regeneration of eruption patterns of CSU using the basal model (1)-(4). The picture was adopted from Seirin-Lee et al. (2022) under a CC-BY-NC-ND 4.0 International license. The dynamics of CSU development were mathematically modelled in Seirin-Lee et al. (2023) based on *in vitro* experimental data. We recalled and reconstructed these to explore the precise effects of each phase on CSU development. Specifically, the basal model is expressed as: **[Endovascular site]** ![Formula][1] ![Formula][2] **[Extravascular site]** ![Formula][3] ![Formula][4] where [*H**B*](**x**, *t*), [*TF*](**x**, *t*), [*C*](**x**, *t*), and [*H**M*](**x**, *t*) are concentrations of histamine released from basophils, tissue factor expressed on vascular endothelial cells, activated coagulation factors leaked from blood vessel, and histamine released from mast cells, respectively. *J**gap*(**x**, *t*) is a function of gap formation estimated from the experiments (Seirin-Lee et al. 2023), and its detailed form is: ![Formula][5] Each *N**i*(*i* = 1, · · ·, 4) and *I**j*(*j* = 1, · · ·, 3) indicate the main physiological networks shown in Fig. 1A and B. Ω*D* ⊂ ℝ2 is the region of the dermis, Ω*E* is the vascular endothelial tissue, and Ω*B*(⊂ Ω*E*) is a random region of the vascular endothelial tissue in which basophils affect the endothelial cells. A certain number of blood vessels are sufficiently and uniformly distributed in the skin and the depth scale at which the vascular reactions involved in wheal formation is negligibly small compared to the scale of the wheal patterns in patients with CSU. Thus, the model is assumed on Ω*D* = Ω*E* = [0, *L*] × [0, *L*] ⊂ ℝ2. All parameter values in the model are assumed to be non-negative constants and *T**sw* *>* [*TF*](0). *χ**b*(**x**, *t*) and *χ**m*(**x**, *t*) are defined as: ![Formula][6] ![Graphic][7] and ![Graphic][8] are the total amounts of histamine contained in basophil and mast cells, respectively (Hattori and Seifert 2017). Because skin eruptions reflect a condition in which substrates leaking from the blood vessels are visible on the skin surface, we define the eruption state function, *S**w*, on the skin (Ω*P* ⊂ℝ2) as defined in Seirin-Lee et al. (2023), such that: ![Formula][9] This function directly reflects the qualitative dynamics of the eruption patterns, which is consistent with the dynamics of the main factors in the model. A previous study by Seirin-Lee et al. (2023) found that the eruption pattern of a CSU can be classified into five types and successfully regenerated *in silico* patterns using the models (1)-(4). Here, we recall the actual eruptions of CSU patients and the representative *in silico* patterns generated by the models of (1)-(4) in Fig. 1C and D. ### 2.2 Basophil/Mast cell positive-feedback loop model By neglecting the extravascular or endovascular site network loops of either *N*1 − *N*3 − *N*4 or *N*1 − *N*2, we reconstructed the basal model (1)-(4) as the basophil/mast cell positive-feedback loop model (Fig. 1B). Specifically, we consider the CSU model (1)-(4) with either [*H**M*] ≡ 0 or [*H**B*] ≡ 0. Using these models, we explored which network loops play a key role in each CSU development phase. The former case corresponds to: **[Basophil positive-feedback loop model]** ![Formula][10] ![Formula][11] ![Formula][12] This represents the positive-feedback loop of basophil via vascular endothelium and coagulation factors. The latter case corresponds to the model of **[Mast cell positive-feedback loop model]** ![Formula][13] ![Formula][14] ![Formula][15] This represents the positive-feedback loop of mast cells via the tissue factor of the vascular endothelium and coagulation factors that leak from the blood vessels. ### 2.3 Histamine-only model of mast cells To investigate the role of mast cell in developing the CSU more precisely, we reduced the mast cell positive-feedback loop model to a core structure under certain assumptions. Let us assume that the concentration of tissue factors quickly approaches a steady state, namely, *d*[*TF*]*/dt* ≈ 0, and assume that the diffusion coefficient of coagulation factors is sufficiently slower than that of histamine released from mast cells and that it quickly approaches a steady state, namely, *d*[*C*]*/dt* ≈ 0. Then, models (8)-(9) can be approximated as a system with a sufficiently small *ε* as follows: ![Formula][16] ![Formula][17] ![Formula][18] Then, with *ε* → 0, we obtain the following equations from equations (11)-(12): ![Formula][19] ![Formula][20] Now, by instituting (11) and (12) into the equation (13), we obtain the reduction system for histamine alone as follows:: ![Formula][21] where ![Formula][22] ### 2.4 Numerical conditions and parameter values A numerical simulation was performed by using the alternating-direction implicit (ADI) method (Morton and Mayers 1994) with the zero-flux boundary conditions for [0, *L*] × [0, *L*]. For the initial conditions, we selected a spatially homogeneous state with small perturbations. The detailed forms are as follows: ![Formula][23] where *ϵ* ≪ 1 and we typically choose 0.01. *η**b* and *η**m* are stimulus intensities, and *ϕ*1(**x**) is a random variable function of uniform distribution [0, 1] on [0, *L*] × [0, *L*]. Because extravascular coagulation factors cannot leak from blood vessels without gap formation, we assumed that: ![Formula][24] Then, we set Ω*D* = Ω*E* = Ω*P* = [0, *L*] × [0, *L*]. To select a random region of the vascular endothelial tissue in which basophils affect endothelial cells, we selected a specific position for endothelial cells by defining Ω*B* such that Ω*B* = {**x**|*ϕ*1(**x**) *> p*} where 0 *< p <* 1. If *p* is close to 1, the position where basophils affect endothelial cells is restricted by a high stimulus point. By changing the value of *p*, the density (denseness or sparseness) of the initial stimuli is represented. The detailed parameter values were chosen from the original model of Seirin-Lee et al. (2023), in which the representative parameter set was estimated using *in vitro* experimental data. To explore the model dynamics, we tested the parameters in several regions. However, we restricted the inhibition rate to less than 1; namely, *α**b*, *α**f*, *α**m* ≤ 1, based on the original model development using *in vitro* experiments. The experiments showed that inhibition by adenosine suppressed the histamine release rate but not to a negative value (shown in Figure 2 of Seirin-Lee et al. (2023)). ![Figure 2:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/11/12/2023.11.11.23298424/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2023/11/12/2023.11.11.23298424/F2) Figure 2: Onset phase of wheal formation in CSU. (A) Simulation results of the basophil-feedback loop model (5)-(7). Dots pattern appears. *μ**f* = 1, *δ**f* = 0.01, *γ**f* = 7, *γ**f* = 2.2, *α**f* = 8.925, *α**f* = 1, *α**f* = 1, *δ**b* = 0.1, *γ**b* = 5.2, *α**b* = 0.00625, and *α**b* = 0.335. The values for initial conditions are *η**b* = 1 and *p* = 0.999. (B) Simulation results of the mast cell positive-feedback loop model (8)-(10). A pattern does not appear. *μ**f* = 1, *δ**f* = 0.01, *γ**f* = 7, *γ**f* = 2.2, *α**f* = 8.925, *α**f* = 1, *α**f* = 1, *δ**m* = 0.01, *γ**m* = 5.0, *α**m* = 0.0035, and *α**b* = 0.865. The value for initial condition is *η**m* = 1000. (C) The simulation result for ODE system (5)-(6) (upper panel) and equilibrium graph of [*H**B*] (lower panel). The same parameter values with (A) were used. The red dot indicates a monostable equilibrium. (D) Possible onset region depending on the inhibition rate *α**f* of TF expression. *μ**f* = 1, *δ**f* = 0.01, *γ**f* = 7, *γ**f* = 2.2, *α**f* = 8.925, *α**f* = 1(red line), 2(blue line), *and*50(green line). The graph is given by the equation (18). (E) Same parameters as those of (D), *α**f* = 1, *δ**b* = 0.1, *γ**b* = 5.2, *andα**b* = 0.00625. The value of the monostable equilibrium ![Graphic][25]was calculated from the equation (19). ### 2.5 Experiment of time courses of annular wheals in urticaria patients Multiple wheals in the CSU (Fig. 5A) and cholinergic urticaria (Fig. S1) groups with an annular eruption pattern (wheal) were analyzed with the approval of the Ethical Committee for Epidemiology of Hiroshima University, Hiroshima, Japan (approval number: E-1008). The picture data for patient A in Fig. 5A was adopted from Seirin-Lee et al. (2020) under a CC-BY 4.0 International license to confirm the overall changes of wheals with CSU. Wheals that emerged from the patient were recorded using a digital camera over multiple time points. The wheal area was calculated from the number of pixels occupying the target wheal on the digital image, with reference to the number of pixels in the standard length/area. ## 3 Result ### 3.1 Onset Phase: Positive-feedback loop of basophils is crucial for the onset of wheal formation in CSU We investigated whether CSU onset is due to basophils at the endovascular site or mast cells at the extravascular site. To explore this hypothesis, we considered the following: either the basophil positive-feedback loop model (5)-(7) or the mast cell positive-feedback loop model (8)-(10), with respect to the initial stage (namely, *χ**b*(**x**, *t*) = *χ**m*(**x**, *t*) = 1). We first tested whether histamine release from basophil is sufficient to induce TF to cause the onset of skin eruptions using the basophil-feedback loop model. We found that the onset of CSU was possible even without the reaction loop of mast cells, although we did not observe any other eruption patterns, except for the dot pattern, in the numerical tests (Fig. 2A). On the other hand, unexpectedly, we were not able to find any development of eruption patterns in the mast cell-feedback loop model, even for higher levels of the initial stimulus (Fig. 2B). These results indicate that the positive-feedback loop of basophils is crucial for the onset of wheal formation in CSU, and that a stimulus of mast cells may not be a sufficient trigger for the onset. To confirm our numerical results, we mathematically analysed both models. In the basophilfeedback loop model, the dynamics of CSU onset are determined by the first two ordinary differential equations (ODE) (5)-(6), because the dynamics of coagulation factors are only affected in one way from TF expression. Thus, the possibility of onset can be understood based on the equilibrium conditions of the ODE system at a local point where basophils stimulate endothelial cells. If a positive minimal equilibrium exists which is stable and greater than the initial concentration, onset always occurs. Thus, we first calculate the equilibrium and then investigate the parameter constraints for the existence of a positive equilibrium. By defining the equilibrium concentrations of [*H**B*] and [*TF*] as ![Graphic][26] and ![Graphic][27], respectively, we obtain: ![Formula][28] and ![Formula][29] from equations (6)≡ 0 and (5)≡ 0, respectively. From the equation (19), we see that a unique positive equilibrium exists if *α**b* *<* 1 and *α**f* *<* 1 as shown in Fig. 2C (See Appendix A for a more detailed analysis). From a mathematical rather than biological perspective, it follows that if *α**f* ≫ 1, the equilibrium concentration of [*TF*] becomes lower than the threshold value, below which gap formation cannot occur, namely, *T**** *< T**sw* (Fig.2D (Green line)). That is, CSU did not develop in this case. On the one hand, the onset should only be possible if the equilibrium is within an appropriate range (*T**** *> T**sw* and![Graphic][30]). Thus, we next investigate the minimal positive equilibrium ![Graphic][31] value from equation (19) by changing *α**b* under the condition *α**f* ≤ 1. We find a stable minimal equilibrium which satisfies ![Graphic][32] and ![Graphic][33] and confirm that onset is possible when *α**b* ≤ 1 (Fig. 2E (blue shaded region)). Taken together, we conclude that, without parameter sensitivity, onset is possible through the basophil positive-feedback loop. Finally, we confirm the simulation results shown in Fig. 2B by analysing the mast-cellfeedback loop model in more detail. If the initial TF expression is not sufficient (i.e. [*TF*] *< T**sw*), gap formation does not occur and coagulation factors cannot be leaked from the blood vessels. Thus, the coagulation factors in the extravascular site decay exponentially, which means [*C*] ≈ 0; consequently, the histamine from mast cells will remain in the initial equilibrium state, namely, *δ**m**/μ**m* from equation (13). This happened in Fig. 2B. However, if we assume that [*TF*] *> T**sw* and the coagulation factors quickly approach the steady state, namely [*C*](**x**, *t*) ≈ *γ**c**/μ**c* from equation (9), then, we obtained a single equation for histamine production from mast cells as follows: ![Formula][34] This assumption can be interpreted such that we consider the dynamics of histamine from mast cells under conditions where the expression of tissue factors is sufficiently strong with respect to even a small positive feedback from mast cells. Now, we consider a spatially perturbed stimulus around a constant steady state(≡ *u**** *>* 0), namely, [*H**M*](**x**, *t*) = *u****(1 + *ϵϕ*(**x**)), where *ϵ* ≪ 1. We can assess that the reaction term in equation (20) has a unique positive equilibrium and the calculation of the linear analysis for the homogeneous steady state confirms a negative eigenvalue (*λ*), such that: ![Formula][35] (See Appendix A for more details). This result indicates that all constant steady states are stable without parameter dependence and that the histamine level released from mast cells always reaches equilibrium if TF expression is large enough to cause gap formation. Therefore, the TF expression state is a critical factor for bifurcation to induce onset in the mast cell positive-feedback loop model. Taken together, locally strong stimulation TF expression via basophil positive-feedback loop induces gap formation and consequently triggers the wheal onset of CSU. However, globally weakly stimulated TF expression via the mast cell-positive-feedback loop may not be sufficient to form a gap, indicating that mast cells may not be fundamental players in the onset of wheal formation in CSU. ### 3.2 Development Phase In the previous section, we found that gap formation is crucial in the onset phase and the basophil dynamics related to TF expression are critical. However, we could not find any other types of patterns except for dots in the basophil positive-feedback loop model. These observations indicate that endovascular site dynamics play a role in triggering the CSU, while additional mechanisms may be involved during the developmental phase of wheal formation or the developmental phase involves the dynamics of extravascular site. Therefore, we asked the following questions. First, what is the key mechanism that induces the spatial heterogeneity (namely, spatially non-uniform appearance) of skin eruptions during the developmental phase? and, what is the main factor creating the diversity of eruption patterns? #### 3.2.1 The switch mechanism of gap formation is crucial for generating the spatial heterogeneity of skin eruptions To investigate the mechanism by which skin eruptions develop spatially and heterogeneously after CSU onset, we focused on the TF-switch role. As we can see in the analysis in the previous section, the level of TF expression is crucial in the onset of CSU, and gap formation plays a crucial role in connecting the endovascular and extravascular sites. In the original CSU model (1)-(4), the gap formation function was estimated by on-off switch type function based on the experimental data of Seirin-Lee et al. (2023). Thus, to understand how the on-off switch type plays a critical role in the overall pathophysiological dynamics of CSU wheal development, we assumed two other types of gap formation functions in addition to the original TF-switch threshold type: a linear and saturation types (Fig. 3A-D, Table 1). A linear type implies that the gap formation occurs in a linearly dependent manner on TF concentration; thus, that there is no threshold-like effect (Fig. 3D, red line). The saturation type represents the effect of a rapid increase and saturation of the gap formation (Fig. 3D, green line). View this table: [Table 1:](http://medrxiv.org/content/early/2023/11/12/2023.11.11.23298424/T1) Table 1: Effect of function types for gap formation on the heterogeneity of eruptions. ![Figure 3:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/11/12/2023.11.11.23298424/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2023/11/12/2023.11.11.23298424/F3) Figure 3: Gap formation functions. (A) Simulation results for the switch-like type in the full model. (B) Simulation results for the linear type in the full model. (C) Simulation results for the saturation type in the full model. (D) The three types of gap formation function. Each *J**gap*([*TF*]) is described in Table 1. *ϵ**κ* = 0.01. (E) Simulation results for the switch-like type in the reduced histamine-only model. (F) Simulation results for the liner type case in the reduced histamine-only model. (G) Comparison of the reaction functions of the reduced histamine-only model with switch-like and linear cases of gap formation. The switch-like case shows the bistability structure in contrast to the linear case of monostable structure. (H) Gap formation (left panels) and reaction functions (right panels) for the weak (blue line) and strong switch-like cases (red line). (I) Simulation results for weak switch-like cases (blue line in (G)) in the reduced histamine-only model. Representative parameters (Table S1) have been used in (A-C). For (E-I), the parameters are the same as those in Table S1 except for *δ**m* = 0.1. We tested the linear and saturation types of gap formation functions (Fig. 3B, C). We found a localized high expression of histamine released from basophil in the endothelial cell tissue, but other factors showed almost spatially homogeneous high concentrations, such as those observed in the case of anaphylaxis. Thus, we were unable to determine the heterogeneity of skin eruptions despite the wheals have been developed. This result indicates that the gap formation dynamics of endothelium are critical for the origin of eruption heterogeneity. The heterogeneous origin of the eruptions was generated by a switch-like system of gap formation. #### 3.2.2 The development phase is primarily regulated by the positive feedback of mast cell dynamics We found that the switch-like dynamics of gap formation is critical for the heterogeneous development of wheal. Thus, we hypothesized that the developmental phase may be regulated mainly by the mast cell positive-feedback loop at extravascular site via the vascular endothelium reaction. To confirm this, we reduced the mast cell-feedback loop model (8)-(9) to model (16), which describes only the dynamics of histamine released from mast cells. Let us recall the model (16): ![Formula][36] where ![Formula][37] Using the reduction system (21), we first confirmed the state of the pattern development using two gap formation functions: ![Formula][38] Similar to the full system of (1)-(4), we found a spatial pattern for the switch-like function ![Graphic][39] but not a linear gap formation function ![Graphic][40], as shown in Fig. 3E and 3F. This result indicates that the spatial heterogeneity of eruption dynamics in the developmental phase can be understood by the dynamics of histamine released from mast cells. Thus, we explored why switch-like dynamics of gap formation were required for the heterogeneous development of eruptions. We investigated the stability of homogeneous steady states to understand the dynamics of the initial histamine concentration. The linear case exhibited a monostable state (Fig. 3G, green line), indicating that a linear case will present either of the eruptions emerging spatially uniformly because the equilibrium state of histamine ![Graphic][41] is greater than the initial concentration, namely, ![Graphic][42](See Appendix B for more detailed mathematical proofs). However, we found that the switch-like function can create a bistable state (Fig. 3G, red line), which is a critical property that induces heterogeneity in the eruption development. To confirm this, we decreased the switch intensity and found that the bistable property changed to monostable (Fig. 3H). We also confirmed that uniform development (non-patterned state) occurred in the monostable state (Fig. 3I). Therefore, we assume that the switch intensity is sufficiently large, namely, *β* → ∞. Then, we have:![Graphic][43], ![Formula][44] and ![Formula][45] as shown in Fig. 3G (left-hand panel). Then, the following inequality holds: ![Formula][46] Therefore, the minimal and maximal positive stable steady states (![Graphic][47] and![Graphic][48], respectively) in equation (28) will be ![Graphic][49] and ![Graphic][50]based on the comparison theorem (Conway and Smoller 1977). This implies that the system has a bi-stability structure and the initial stimulus size *η**m**ϵϕ*3(**x**) of equation (17) is greater than ![Graphic][51] where ![Graphic][52] is an unstable equilibrium. [*H**M*](·, *t*) will go to the other stable equilibrium state![Graphic][53]. Otherwise, it remained at![Graphic][54]. Taken together, a bistable structure with diffusion is likely to cause the development of patterns and the positive feedback dynamics of histamine released from mast cells via switch-like gap formation play a crucial role in wheal development phase of CSU. #### 3.2.3 Diversity of eruption patterns is largely attributed to the positive-feedback loop of mast cell dynamics Finally, we investigated the essential mechanisms for creating diverse eruption patterns. As the heterogeneous development of an eruption is regulated by the positive feedback dynamics of mast cells, we hypothesized that the diversity of eruption patterns could be explained by mast cell dynamics. Thus, we numerically tested the types of eruption patterns that appeared in the mast cell histamine-only model (28) by changing the parameter values. Interestingly, we found four types of eruption patterns, which were classified by EGe Criteria developed in Seirin-Lee et al. (2023). The EGe criteria defines five types of eruption patterns: annular, broken-annular, geographic, circular, and dotted. As seen in Fig. 4A, we found four types of patterns except for the dot patterns. We also confirmed that pattern diversity occurs owing to the balance of the three parameters of histamine dynamics from mast cells and TF dynamics in endovascular endothelial cells; the histamine release rate (*γ**m*), histamine decay rate (*μ**m*), and TF activation rate (*γ**f*)) (Fig. 4B). Taken together, the positive-feedback loop of mast cell dynamics is essential in the diversity of eruption patterns, and the balance between activation and inhibition networks in mast cells may be critically involved in eruption diversity. ![Figure 4:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/11/12/2023.11.11.23298424/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2023/11/12/2023.11.11.23298424/F4) Figure 4: Diverse pattern of eruptions generated by the histamine-only model. (A) Representative simulation examples of eruption patterns. Annular (*γ**m* = 0.6, *μ**m* = 1.0, *γ**f* = 7.0), Broken-annular (*γ**m* = 0.42, *μ**m* = 1.0, *γ**f* = 7.0), Circular (*γ**m* = 0.6, *μ**m* = 0.1, *γ**f* = 4.2), and Geographic patterns (*γ**m* = 1.0, *μ**m* = 0.1, *γ**f* = 4.2) were used and the other parameters are listed in Table S1. (B) Example pattern for each parameter set. The initial condition of this model was given to [*H**M*](**x**, 0) = *δ**m**/μ**m* + *η**m**ϕ*3(**x**). ### 3.3 Disappearance Phase: Latent suppression of mast cells may be involved in disappearance of eruptions One of the characteristics of wheal formation in CSU is the rapid development and disappearance of skin eruptions on a minute or hour scale. Although the pathophysiological mechanism of CSU development was partially or integrally understood in our previous studies, including experiments and mathematical modeling, the underlying mechanism by which urticaria disappears remains elusive. The original models (1)-(4) suggested by Seirin-Lee et al. (2023) successfully regenerated the eruption patterns in the development phase; however, this model failed to regenerate the disappearance phase of the eruptions, except for the dot pattern case. The eruption patterns in this mathematical model generally spread throughout the entire spatial domain if there were pre-simulated mast cells, suggesting that there is a disappearance mechanism separate from the development mechanism of wheals in urticaria. Thus, we chose a mathematical approach to understand the disappearance mechanism. We inferred the mechanism by comparing the quantitative and qualitative dynamics of eruption patterns in a mathematical model with the actual dynamics of wheal patterns in urticaria patients. We first extracted three features of wheals from the photo data of CSU patients, both quantitatively and qualitatively, as shown in Fig. 5A. ![Figure 5:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/11/12/2023.11.11.23298424/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2023/11/12/2023.11.11.23298424/F5) Figure 5: Disappearance phase of wheal formation in CSU. (A) Evolutionary dynamics of wheal (eruption) patterns in a patient *P*1 and *P*2 with CSU (left panels) and the weal area expansion (right panel). **W**1, **W**2, **W**3 and **U**1, **U**2 indicate each wheal from the top in each patient figure. The picture of patient A and the data of weal area expansion were adopted from Seirin-Lee et al. (2020) under a CC-BY 4.0 International license. (B-H) *in silico* experiments for the cases of suppression of TF expression, plasma leakage (namely, the intercellular gaps of vascular endothelial cell in blood vessels become small), and mast cells in releasing histamine. (E) Multiple wheals emerge and disappear simultaneously. (F-H) Weal area expansion data calculated from Figs. B, C, D, respectively. The parameter values are listed in Table S1. 1. Wheals that appeared simultaneously will *disappear simultaneously* (the wheals **W**1, **W**2, and **W**3 in left panels in Fig. 5A, Patient *P*1). 2. The boundary of annular wheal becomes *broken* in the disappearance phase (see the enlarged pictures in left panels). 3. *The wheal expansion almost stops or slightly decrease* in the disappearance phase. (Fig.5A, grey zone (DP) in right graph for patient *P*1 and the right graph for patient *P*2). First, based on hints from feature (i), we developed the simplest suppression model that acts uniformly over a wide area after urticaria has progressed to a certain stage: ![Formula][55] where *μ**mc* is the suppression rate and *β**c* is a positive constant that controls the transition state at![Graphic][56]. Here, *I**s*(*t*) is defined as: ![Formula][57] where [*A*] is either a coagulation factor ([*A*] = [*C*]) leaked from the blood vessel or histamine ([*A*] = [*H**M*]) released from the mast cells in the dermis. This modeling formulation implies that a suppression network starts after the development phase of the CSU, and its activation timing is determined not locally for individual wheals, but globally for wheals in the area, depending on the state of each wheal progression. Next, we extend the basal model (1)-(4) by combining the suppression terms as follows: **Case 1** Suppression of TF expression on vascular endothelial cells: (**S****1**(*t*), **S****2**(*t*), **S****3**(*t*)) = (*S**p*[*I**s*(*t*)], 0, 0). **Case 2** Suppression of plasma leakage from blood vessels: (**S****1**(*t*), **S****2**(*t*), **S****3**(*t*)) = (0, *S**p*[*I**s*(*t*)], 0). **Case 3** Suppression of mast cells in releasing histamine: (**S****1**(*t*), **S****2**(*t*), **S****3**(*t*)) = (0, 0, *S**p*[*I**s*(*t*)]). such that: ![Formula][58] ![Formula][59] ![Formula][60] ![Formula][61] With the disappearance of the phase-combined model (22)-(25), we explored the model case that captured the characteristics of the disappearance phase shown in patients with urticaria in Fig. 5A. Interestingly, we found that only Case 3 succeeded in capturing the wheal evolution characteristics of the CSU patient data, whereas both Cases 1 and 2 failed to explain them (Fig. 5B-H). In Case 1 (Fig. 5B, F), the weal expansion continued in the disappearance phase and a broken pattern was not seen. Features (ii) and (iii) were not found. In Case 2 (Fig. 5C, G), the weal thickness becomes thinner during the disappearance phase, but the wheal patterns did not completely disappear and continuously expended with keeping low concentration of histamine. In contrast, Case 3 satisfied all three features (i)-(iii). The wheal disappeared with a broken annual pattern, and the expansion of wheal almost stops or slightly decreased during the disappearance phase (Fig. 5D, H) (See Fig. S1 in Appendix C for additional clinical data). We also confirmed that Case 3 satisfies wheal patterns disappeared globally (Fig. 5E). Taken together, we concluded that TF expression or plasma leakage dynamics by gap formation are not be involved, or do not dominantly contribute in the disappearance phase, while the suppression of histamine release from mast cells may play a critical role in the disappearance of wheals in many patients with urticaria. ## 4 Discussion CSU is an intractable skin disease that involves skin mast cells, basophils, and the blood coagulation and complement systems. IgG autoantibodies to IgE or Fc*ε*RI and/or IgE autoallergy to IL-24, thyroid peroxidase or other endogenous molecules may be involved in the activation of mast cells and basophils (Huang et al. 2020; Maurer et al. 2021; Rijavec et al. 2021; Yanase et al. 2021b). However, due to a lack of proper experimental animal model, the precise mechanism *in vivo* has remained elusive. No individualized medications for each CSU patient has been established, especially to achieve a cure of the disease. A recent mathematical approach (Seirin-Lee et al. 2023) suggested a CSU model, through which the skin eruption dynamics in CSU can be captured. To understand the precise mechanism of CSU development, we focused on the wheal dynamics of CSU by structurally analyzing a mathematical model. We explored the pathophysiological dynamics of wheal formation in three phases: the onset, development, and disappearance phases, and found that the pathophysiological networks play a different role in each phase. For the onset phase, we found that the basophil positive-feedback loop is critical for evoking the eruptions, but the small stimuli for mast cells were not sufficient to induce wheals in CSU. This observation indicates that the onset of CSU is mainly caused by the reaction of basophils on vascular endothelial cells. The increase of basophil chemoattractant released from skin tissues and/or activation of basophils by autoimmune mechanisms may sufficiently increase the probability of basophils to attach to and activate vascular endothelial cells triggering wheal formation via tissue factor expression. It can also be speculated that the number of basophils in the blood may decrease by moving out of blood vessels to dermal tissues. In fact, the number of basophils in the blood of CSU, especially those who develop with severe symptoms and refractory to treatments, tends to be low (Huang et al. 2020; Rijavec et al. 2021). During the developmental phase, we found that the mast cell positive-feedback loop was crucial for spatial heterogeneity and various patterns. In particular, we confirmed that the switch-like dynamics of gap formation play a critical role in generating spatial heterogeneity and creating a bistable structure in the mathematical model (Fig. 3). This result suggests that the regulation of vascular permeability by mast cell activation via gap formation dynamics of the vascular endothelium is important. Furthermore, the balance of positive feedback network of mast cells and dynamics of gap formation play an important role in generating various patterns of CSU eruption (Fig. 4B). This also supports the medical understanding that wheal formation in urticaria is primarily driven by mast cells (Zuberbier et al. 2022). During the development phase, we also found that all four patterns types of eruption reported in the previous study (EGe Criteria, Seirin-Lee et al. (2023)), except for the dot pattern, could be generated in the reduced mast cell model of this study. This indicates that the major, if not all, mechanism of the dot pattern is different from that of the other four patterns. Thus, the efficacy of antihistamine may vary depending on the type of wheals, and that the additional or alternative treatments may be required for patients with dots pattern. The relationship of the efficacy of antihistamine with different patterns of wheals is warranted for future study. In Fig. 4, we often find a high frequency tendency for geographic patterns in the chosen parameter regions. Indeed, the statistical data of patients with CSU on eruption types shown in a previous study, Seirin-Lee et al. (2023), support that geographic patterns were the most dominant patterns found in patients with CSU (Fig. 4B in Seirin-Lee et al. (2023)). The frequency of pattern types should be analyzed more precisely in a future study. We further confirmed that the positive-feedback loop of mast cell dynamics is essential in the diversity of wheal patterns, and the balance between activation and inhibition networks in mast cells may be critically involved in wheal diversity. This is consistent with the results proposed in the conceptual urticaria model of our previous work, Seirin-Lee et al. (2020). Mathematical insights in this study also showed that bistability based on gap formation function is important in creating spatial heterogeneity in wheal development. A similar structure was found in the previous conceptual model of Seirin-Lee et al. (2020), although the detailed form of the model equation differed. This suggests that an intrinsic and universal mechanism for the diversity of wheals in urticaria is due to mast cell-related dynamics, and that additional networks diverging from the intrinsic structure may give rise to urticarial subtypes. Finally, our study hypothesized that the disappearance phase of wheals in CSU occurs via a globally independent suppression mechanism followed by the development of wheals and we theoretically validated this by using wheal expansion data of urticaria patients. The *in silico* experiments showed that there may be a mechanism to suppress mast cells from releasing histamine over a wide area where the wheals have developed during the disappearance phase. This suggests that the disappearance mechanism is not a passive dynamic by a local histamine deficiency from stimulated mast cells but an active immune reaction which occurs in a global area over mast cells. Although there is no experimental evidence of the disappearance phase, this hypothesis suggests a novel concept to understand the mechanism of wheal disappearance in that there may exist a global suppression of mast cells after a time lag of wheal development. Mast cells express inhibitory molecules, such as Siglec-8, Siglec-6 and CD200R. The activation of these molecules expressed on mast cells and diminishes histamine release (Metz et al. 2023). The low-affinity IgG receptor (Fc*γ*R2b), which contributes to supra-optimal negative regulation of Fc*ε*RI-induced mast cell activation is another candidate of a player in the disappearance phase (Gast et al. 2018). It would be valuable to prove the precise of pathophysiological network involved in the disappearance phase. In this study, we explored the overall dynamics of the CSU with respect to the three-phase mechanism of wheal formation by precisely analyzing the CSU mathematical model. Although all the results are based on a mathematical perspective, they suggest a concept of narrowing down the possibility for experimental verification to reveal the precise pathology of CSU. Our theoretical observation of the relationship between the phases of wheal dynamics in CSU and pathophysiological networks may contribute to finding an optimal therapy for patient-specific treatment based on the spatiotemporal dynamics of wheal formation. ## Data Availability All data produced in the present work are contained in the manuscript [https://github.com/seirin-lee/urticaria-code](https://github.com/seirin-lee/urticaria-code) ## Competing interest statement The authors declare no conflicts of interest. ## Data availability All relevant data are included in the manuscript and the supporting information files. ## Code availability Numerical code is made available online on GitHub ([https://github.com/seirin-lee/urticaria-code](https://github.com/seirin-lee/urticaria-code)) and all other data are available from the corresponding author (SSL) on reasonable request. ## Author Contributions SSL initiated, designed, and executed the research. SSL developed a mathematical model and numerical algorithms for simulations and analyzed data. ST executed the clinical experiment and provided the data. MH provided photographs of patients. SSL wrote the first draft. SSL, ST, and MH revised the manuscript. ## Acknowledgements This study was supported by the Japan Science and Technology Agency(JST) CREST (JP-MJCR2111), and grants to MH from MEXT, Grant-in-Aid for Scientific Research (Grant Number 21K08346). ## Appendix ### A Proofs of stability in onset phase We consider the two stability problems, independently. One is the equation (19) which has been reduced from the basophil positive feedback loop model (5)-(6). We recall it here; ![Formula][62] The other is the equation (20) ![Formula][63] where ![Formula][64] which has been reduced from the mast cell positive feedback loop. Note that the solutions satisfying ℱ*B*([*H**B*]) = 0 and 𝒢 ℳ ([*H**M*]) = 0 are constant steady states. We have the following lemmas. Lemma A.1 *Suppose that α**f* ≤ 1 *and α**b* ≤ 1. *Then the minimal positive equilibrium of the equation* (26) *is stable. And the solution started from* [*H**B*](0) *converges to the minimal positive equilibrium*. Proof. Assume that there exists an equilibrium of *d*[*H**B*]*/dt* = 0 in the equation (26) and let us define the minimal positive equilibrium by![Graphic][65]. Since ℱ*B*(0) *>* 0 and ℱ*B*(∞) = −∞, ![Formula][66] should be hold. Hence, ![Graphic][67] is stable. Furthermore, from (26), we obtain ![Formula][68] because *α**f* ≤ 1 and *α**b* ≤ 1. Since ℱ*B*([*H**B*](0)) *>* 0, the solution started from [*H**B*](0) always converges to a stable ![Graphic][69]. Lemma A.2 *There exists a unique positive constant steady state in the equation* (27) *that is stable for all positive parameters*. Proof. Since 𝒢 ℳ (0) *>* 0 and ![Formula][70] for all positive parameters, 𝒢 ℳ ([*H**M*]) is decreasing in monotone and there exist a unique positive constant steady state. Now, let us denote a constant steady state by [*H**M*](**x**, *t*) = *u****. Instituting *z*(**x**, *t*) = *u**** + *ϕ* exp((2*π/L*)**n·x***i* + *λt*), where **n** is a wave mode, into the equation (27) for a standard linear analysis, we have ![Formula][71] without any constraint of parameter conditions except for the positiveness. Hence, *u**** is a stable steady state. Remark A.1 *The reduced model equation* (27) *has been reduced under the condition of sufficient expression of TF. However, if the TF expression is not sufficient, the coagulation factors* [*C*] *will quickly converge to* 0. *Therefore, the reduced model* (27) *becomes* ![Formula][72] *Thus, u**** = *δ**m**/μ**m*(≡ [*H**M*](**x**, 0)) *is always stable. Thus, onset will not occurs*. ### B Some properties in development phase We recall the histamine alone model on development phase so that *χ**m*(**x**, *t*) = 1 as the following. ![Formula][73] where ![Formula][74] Lemma B.1 *Let us define a constant steady state of the equation* (28) *with* ![Graphic][75] *by* ![Graphic][76]. *Then all positive equilibria satisfy* ![Graphic][77] *for α**f* ≤ 1 *and α**m* ≤ 1. Proof. From the reaction terms of (28), we have ![Formula][78] Since we have ![Formula][79] the following inequality is hold; ![Formula][80] Lemma B.2 *Let us suppose that*![Graphic][81] *is sufficiently small such that* ![Formula][82] *then the equation* (28) *with* ![Graphic][83] *has a unique constant steady state which is stable*. Proof. From the equation (28), we can easily check ℛ(0) *>* 0 and ℛ(∞) *<* 0. Thus, the proof is completed if we show that ℛ ([*H**M*]) is a function of monotone decreasing. The reaction term of the equation (28) denoted by ℛ ([*H**M*]) satisfies ![Formula][84] Then, we have ![Formula][85] Remark B.1 *The two lemmas suggest that if TF expression is constantly induced and the gap formation is linearly dependent on TF expression, there does not exist a spatial eruption patten. Namely, CSU occurs but it is spatially uniform such as the case of Anaphylactic shock*. ### C Additional experiment data and simulation results The dynamics of wheal expansion was observed in cholinergic urticaria patient and its dynamics shows that the expansion stops in the disappearance phase (Fig. S1 A). We measured the expansion of boundaries with wheal radius for both inner and outer boundaries. The patient data (Fig. S1 A lowest panels) shows that the boundary thickness is narrowing down. On the other hand, *in silico* data (Fig. S1 B) shows that the boundary thickness is narrowing down for the cases of Case 1 and Case 3 in disappearance phase but not in Case 2. Once the disappearance phase enters, the outer boundary stops in contrast that the inner boundary continuously expands as shown in Case 3. Thus, Case 3 show notable shrink of wheal thickness, as seen in both CUS patient in Fig. 5A and cholinergic urticaria patient (Fig. S1 A). ### D Table of parameter values used in simulations The representative parameter values used in Fig. 3-5 are shown in Table S1. ![Figure S1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/11/12/2023.11.11.23298424/F6.medium.gif) [Figure S1:](http://medrxiv.org/content/early/2023/11/12/2023.11.11.23298424/F6) Figure S1: Weal expansion and boundary thickness. (A) Wheal expansion dynamics of patients with cholinergic urticaria. The label of numbers indicates each wheal. The low four panels are quantitative data of wheal area and radius expansion which are calculated from the patient in top photos. The error bar is a standard deviation. (B) *in silico* data for the radius expansion of the wheals calculated in the same simulation of Fig. 5B-D with respect to outer and inner boundaries. Since the wheal expands isotopically, the radius of the wheal was calculated as the distance between the center of the wheal and the point of each boundary when *y*-axis value is the center of the wheal. View this table: [Table S1:](http://medrxiv.org/content/early/2023/11/12/2023.11.11.23298424/T2) Table S1: Representative Parameter set used in simulations * Received November 11, 2023. * Revision received November 11, 2023. * Accepted November 12, 2023. * © 2023, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. Bazilai, A., Sagi, L., Baum, S., Trau, H., Schvimer, M., Barshack, I., Solomon, M., 2017. The histopathology of urticaria revisited-clinical pathological study. Am J Dermatopathology 39, 753–759. 2. Conway, E. D., Smoller, J. A., 1977. A comparison technique for systems of reaction-diffusion equations. Communications in Partial Differential Equations 2 (7), 679–697. 3. Gast, M., Preisinger, C., Nimmerjahn, F., Huber, M., 2018. IgG-independent co-aggregation of FcεRI and FcγRIIB results in LYN- and SHIP1-dependent Tyrosine Phosphorylation of FcγRIIB in murine bone marrow-derived mast cells. Frontiers in immunology 9 (1937), doi: 10.3389/fimmu.2018.01937. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fimmu.2018.01937&link_type=DOI) 4. Greaves, W. M., 2014. Pathology and classification of urticaria. Immunol Allergy Clin North Am. 34 (1), 1–19. 5. Hattori, Y., Seifert, R., 2017. Histamine and Histamine Receptors in Health and Disease (Part IV: Church MK, Allergy, histamine and antihistamines). Springer, Cham. 6. Hide, M., Kaplan, A. P., 2022. Concise update on the pathogenesis of chronic spontaneous urticaria (csu). Journal of Allergy Clinical Immunology 141 (S0091-6749), 01172–1. 7. Huang, A. H., Chichester, K. L., Saini, S. S., 2020. Association of basophil parameters with disease severity and duration in chronic spontaneous urticaria (csu). J Allergy Clin Immunol Pract. 8 (2), 793–795. 8. Itakura, A., Tani, Y., Kaneko, N., Hide, M., 2018. Impact of chronic urticaria on quality of life and work in japan: Results of a real-world study. Journal of Dermatology 45, 963–970. 9. Karasuyama, H., Mukai, K., Tsujimura, Y., Obata, K., 2009. Newly discovered roles for basophils: a neglected minority gains new respect. Nature Review Immunology 9, 9–13,doi:10.1038/nri2458. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nri2458&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19039320&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F12%2F2023.11.11.23298424.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000261949700010&link_type=ISI) 10. Kolkhir, P., Gimenez-Arnau, A. M., Kulthanan, K., Peter, J., Metz, M., Maurer, M., 2022. Urticaria. Nature Review Disease Primers 8 (61), doi:10.1038/s41572–022–00389–z. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41572–022–00389–z&link_type=DOI) 11. Kurashima, Y., Kiyono, H., 2014. New era for mucosal mast cells: their roles in inflammation, allergic immune responses and adjuvant development. Experimental & Molecular Medicine 46, e83. 12. Lee, S. J., Ha, E. K., Jee, H. M., et al., 2017. Prevalence and risk factors of urticaria with a focus on chronic urticaria in children. Allergy Asthma Immunology Res. 9, 212–219. 13. Leiter, U., Keim, U., Garbe, C., 2020. Epidemiology of skin cancer: Update 2019. Sunlight, Vitamin D and Skin Cancer, Advances in Experimental Medicine and Biology 1268, Springer, Cham. doi: 10.1007/978–3–030–46227–7–6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/978–3–030–46227–7–6&link_type=DOI) 14. Matsuo, Y., Yanase, Y., Irifuku, R., Ishii, K., Kawaguchi, T., Takahagi, S., et al., 2018. The role of adenosine for ige receptor-dependent degranulation of human peripheral basophils and skin mast cells. Allergology International 67, 524–526, doi:10.1016/j.alit.2018.03.007. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.alit.2018.03.007&link_type=DOI) 15. Maurer, M., Khan, D. A., Elieh Ali Komi D., Kaplan, A. P., 2021. Biologics for the use in chronic spontaneous urticaria: When and which. J Allergy Clin Immunol Pract. 9 (3), 1067–1078. 16. Maurer, M., Weller, K., Bindslev-Jensen, C., Gimenex-Arnau, A., Bousquet, P. J., Bousquet, J., et al., 2011. Unmet clinical needs in chronic spontaneous urticaria. Allergy 66, 317–330, doi:10.1111/j.1398–9995.2010.02496.x. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1398-9995.2010.02496.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21083565&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F12%2F2023.11.11.23298424.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000286837100002&link_type=ISI) 17. Maxim, E., Aksut, C., Tsoi, D., Dellavalle, R., 2018. Global burden of urticaria: Insights from the 2016 global burden of disease study. J. Am. Acad. Dermatol 79, 567–569. 18. Metz, M., Kolkhir, P., Altrichter, S., Siebenhaar, F., Levi-Schaffer, F., Youngblood, B. A., Church, M. K., Maurer, M., 2023. Mast cell silencing: A novel therapeutic approach for urticaria and other mast cell-mediated diseases. Allergy 00, DOI: 10.1111/all.15850. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/all.15850&link_type=DOI) 19. Morton, K. W., Mayers, D. F., 1994. Numerical Solution of Partial Differential Equations. Cambridge University Press. 20. Rijavec, M. M M. K., Koren, A., Kopac, P., Selb, J., Vantur, R., et al., 2021. A very low number of circulating basophils is predictive of a poor response to omalizumab in chronic spontaneous urticaria. Allergy 76 (4), 1254–1257. 21. Rudich, N., Ravid, K., Sagi-Eisenberg, R., 2012. Mast cell adenosine receptors function: A focus on the a3 adenosine receptor and inflammation. Front Immunology 3, 134. 22. Seidel, J. A., Otsuka, A., Kabashima, K., 2018. Anti-PD-1 and Anti-CTLA-4 therapies in cancer: Mechanisms of action, efficacy, and limitations. Frontiers in Ongology 8 (86), doi: 10.3389/fonc.2018.00086. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fonc.2018.00086&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F12%2F2023.11.11.23298424.atom) 23. Seirin-Lee, S., Matsubara, D., Yanase, Y., Kunieda, T., Takahagi, S., Hide, M., 2022. Mathematical-structure based morphological classification of skin eruptions and linking to the pathophysiological state of chronic spontaneous urticaria. MedRxive, doi:10.1101/2022.11.04.22281917. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMi4xMS4wNC4yMjI4MTkxN3YxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTEvMTIvMjAyMy4xMS4xMS4yMzI5ODQyNC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 24. Seirin-Lee, S., Matsubara, D., Yanase, Y., Kunieda, T., Takahagi, S., Hide, M., 2023. Mathematical-based morphological classification of skin eruptions corresponding to the pathophysiological state of chronic spontaneous urticaria. Communications Medicine, In press. 25. Seirin-Lee, S., Yanase, Y., Takahagi, S., Hide, M., 2020. A single reaction-diffusion equation for the multifarious eruptions of urticaria. PLOS Computational Biology 16 (1), e1007590. 26. Stone, K. D., Prussin, C., Metcalfe, D. D., 2010. IgE, mast cells, basophils, and eosinophils. Journal of Allergy Clinical Immunology 125 (S73-80), doi:10.1016/j.jaci.2009.11.017. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaci.2009.11.017&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20176269&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F12%2F2023.11.11.23298424.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000280170600007&link_type=ISI) 27. Thandi, C. S., Whittam, L., 2021. Diagnosis and management of common viral skin infections. Prescriber 32 (4), 10–14, doi:10.1002/psb.1907. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/psb.1907&link_type=DOI) 28. Weidinger, S., Beck, L. A., Bieber, T., Kabashima, K., Irvine, A. D., 2018. Atopic dermatitis. Nature Reviews Disease Primers 4 (1), doi:10.1038/s41572–018–0001–z. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41572–018–0001–z&link_type=DOI) 29. Yanase, Y., Matsuo, Y., Takahagi, S., Kawaguchi, T., Uchida, K., Ishii, K., et al., 2021a. Coagulation factors induce human skin mast cell and basophil degranulation via activation of complement 5 and the c5a receptor. Journal of Allergy Clinical Immunology 147, 1101–4. 30. Yanase, Y., Morioke, S., Iwamoto, K., Takahagi, S., Uchida, K., Kawaguchi, T., et al., 2018a. Histamine and toll-like receptor ligands synergistically induce endothelial cell gap formation by the extrinsic coagulating pathway. Journal of Allergy Clinical Immunology 141 (e7), 1115–1118,doi:10.1016/j.jaci.2017.07.026. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaci.2017.07.026&link_type=DOI) 31. Yanase, Y., Takahagi, S., Hide, M., 2018b. Chronic spontaneous urticaria and the extrinsic coagulation system. Allergology International 67, 191–194, doi:10.1101/2022.11.04.22281917. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1101/2022.11.04.22281917&link_type=DOI) 32. Yanase, Y., Takahagi, S., Ozawa, K., Hide, M., 2021b. Role of coagulation and complement factors for mast cell activation in the pathogenesis of chronic spontaneous urticaria. Cells 10 (7), doi:10.3390/cells10071759. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/cells10071759&link_type=DOI) 33. Zuberbier, T., Latiff, A., Abuzakouk, M., Aquilina, S., et al., 2022. The international EAACI/GA2LEN/EuroGuiDerm/APAAACI guideline for the definition, classification, diagnosis, and management of urticaria. Allergy 77 (734-766), doi:10.1111/all.15090. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/all.15090&link_type=DOI) [1]: /embed/graphic-2.gif [2]: /embed/graphic-3.gif [3]: /embed/graphic-4.gif [4]: /embed/graphic-5.gif [5]: /embed/graphic-6.gif [6]: /embed/graphic-7.gif [7]: /embed/inline-graphic-1.gif [8]: /embed/inline-graphic-2.gif [9]: /embed/graphic-8.gif [10]: /embed/graphic-9.gif [11]: /embed/graphic-10.gif [12]: /embed/graphic-11.gif [13]: /embed/graphic-12.gif [14]: /embed/graphic-13.gif [15]: /embed/graphic-14.gif [16]: /embed/graphic-15.gif [17]: /embed/graphic-16.gif [18]: /embed/graphic-17.gif [19]: /embed/graphic-18.gif [20]: /embed/graphic-19.gif [21]: /embed/graphic-20.gif [22]: /embed/graphic-21.gif [23]: /embed/graphic-22.gif [24]: /embed/graphic-23.gif [25]: F2/embed/inline-graphic-3.gif [26]: /embed/inline-graphic-4.gif [27]: /embed/inline-graphic-5.gif [28]: /embed/graphic-25.gif [29]: /embed/graphic-26.gif [30]: /embed/inline-graphic-6.gif [31]: /embed/inline-graphic-7.gif [32]: /embed/inline-graphic-8.gif [33]: /embed/inline-graphic-9.gif [34]: /embed/graphic-27.gif [35]: /embed/graphic-28.gif [36]: /embed/graphic-31.gif [37]: /embed/graphic-32.gif [38]: /embed/graphic-33.gif [39]: /embed/inline-graphic-10.gif [40]: /embed/inline-graphic-11.gif [41]: /embed/inline-graphic-12.gif [42]: /embed/inline-graphic-13.gif [43]: /embed/inline-graphic-14.gif [44]: /embed/graphic-34.gif [45]: /embed/graphic-35.gif [46]: /embed/graphic-36.gif [47]: /embed/inline-graphic-15.gif [48]: /embed/inline-graphic-16.gif [49]: /embed/inline-graphic-17.gif [50]: /embed/inline-graphic-18.gif [51]: /embed/inline-graphic-19.gif [52]: /embed/inline-graphic-20.gif [53]: /embed/inline-graphic-21.gif [54]: /embed/inline-graphic-22.gif [55]: /embed/graphic-39.gif [56]: /embed/inline-graphic-23.gif [57]: /embed/graphic-40.gif [58]: /embed/graphic-41.gif [59]: /embed/graphic-42.gif [60]: /embed/graphic-43.gif [61]: /embed/graphic-44.gif [62]: /embed/graphic-45.gif [63]: /embed/graphic-46.gif [64]: /embed/graphic-47.gif [65]: /embed/inline-graphic-24.gif [66]: /embed/graphic-48.gif [67]: /embed/inline-graphic-25.gif [68]: /embed/graphic-49.gif [69]: /embed/inline-graphic-26.gif [70]: /embed/graphic-50.gif [71]: /embed/graphic-51.gif [72]: /embed/graphic-52.gif [73]: /embed/graphic-53.gif [74]: /embed/graphic-54.gif [75]: /embed/inline-graphic-27.gif [76]: /embed/inline-graphic-28.gif [77]: /embed/inline-graphic-29.gif [78]: /embed/graphic-55.gif [79]: /embed/graphic-56.gif [80]: /embed/graphic-57.gif [81]: /embed/inline-graphic-30.gif [82]: /embed/graphic-58.gif [83]: /embed/inline-graphic-31.gif [84]: /embed/graphic-59.gif [85]: /embed/graphic-60.gif