Research article Special Issues

Dynamics and optimal control of tuberculosis model with the combined effects of vaccination, treatment and contaminated environments


  • Received: 11 January 2024 Revised: 23 February 2024 Accepted: 28 February 2024 Published: 07 March 2024
  • Tuberculosis has affected human beings for thousands of years, and until today, tuberculosis still ranks third among 29 infectious diseases in China. However, most of the existing mathematical models consider a single factor, which is not conducive to the study of tuberculosis transmission dynamics. Therefore, this study considers the combined effects of vaccination, treatment, and contaminated environments on tuberculosis, and builds a new model with seven compartments of SVEITRW based on China's tuberculosis data. The study shows that when the basic reproduction number R0 is less than 1, the disease will eventually disappear, but when R0 is greater than 1, the disease may persist. In the numerical analysis part, we use Markov-chain Monte-Carlo method to obtain the optimal parameters of the model. Through the next generation matrix theory, we calculate that the R0 value of tuberculosis in China is 2.1102, that is, if not controlled, tuberculosis in China will not disappear over time. At the same time, through partial rank correlation coefficients, we find the most sensitive parameter to the basic reproduction number R0. On this basis, we combine the actual prevalence of tuberculosis in China, apply Pontryagin's maximum principle, and perform cost-effectiveness analysis to obtain the conditions required for optimal control. The analysis shows that four control strategies could effectively reduce the prevalence of TB, and simultaneously controlling u2,u3,u4 is the most cost-effective control strategy.

    Citation: Tao-Li Kang, Hai-Feng Huo, Hong Xiang. Dynamics and optimal control of tuberculosis model with the combined effects of vaccination, treatment and contaminated environments[J]. Mathematical Biosciences and Engineering, 2024, 21(4): 5308-5334. doi: 10.3934/mbe.2024234

    Related Papers:

    [1] Anastasiia D. Kurenkova, Peter S. Timashev . Mast cells: A dark horse in osteoarthritis treatment. AIMS Allergy and Immunology, 2022, 6(4): 228-247. doi: 10.3934/Allergy.2022017
    [2] Andrey Mamontov, Alexander Polevshchikov, Yulia Desheva . Mast cells in severe respiratory virus infections: insights for treatment and vaccine administration. AIMS Allergy and Immunology, 2023, 7(1): 1-23. doi: 10.3934/Allergy.2023001
    [3] Peter Jandus . Arachidonic acid metabolism and its use in the diagnosis of mastocytosis. AIMS Allergy and Immunology, 2023, 7(1): 82-91. doi: 10.3934/Allergy.2023006
    [4] Sammi Wong, Sara Hasan, Christina Parducci, Bernadette Ann Riley . The gastrointestinal effects amongst Ehlers-Danlos syndrome, mast cell activation syndrome and postural orthostatic tachycardia syndrome. AIMS Allergy and Immunology, 2022, 6(2): 19-24. doi: 10.3934/Allergy.2022004
    [5] Amolak S Bansal, Alex Nicholas, Nazira Sumar, Veronica Varney . Mast cells, mediators, and symptomatic activation. AIMS Allergy and Immunology, 2024, 8(1): 34-55. doi: 10.3934/Allergy.2024004
    [6] Haiyan Guo, Mingsheng Lei, Jinan Ma, Hongchun Du, Youming Zhang . Allergy and ferroptosis. AIMS Allergy and Immunology, 2025, 9(1): 8-26. doi: 10.3934/Allergy.2025002
    [7] Federica Serra, Lisa Aielli, Erica Costantini . The role of miRNAs in the inflammatory phase of skin wound healing. AIMS Allergy and Immunology, 2021, 5(4): 264-278. doi: 10.3934/Allergy.2021020
    [8] Masahiro Seike . Histamine suppresses T helper 17 responses mediated by transforming growth factor-β1 in murine chronic allergic contact dermatitis. AIMS Allergy and Immunology, 2018, 2(4): 180-189. doi: 10.3934/Allergy.2018.4.180
    [9] Omer Qibi, Severine Audusseau, Andrea Mogas, Zoulfia Allakhverdi, Abdelilah Soussi Gounni, Saba Al Heialy, Qutayba Hamid . No evidence for IgE receptor FcεRI expression on bronchial epithelial cells of asthmatic patients. AIMS Allergy and Immunology, 2018, 2(4): 165-179. doi: 10.3934/Allergy.2018.4.165
    [10] Ken S. Rosenthal, Daniel H. Zimmerman . J-LEAPS vaccines elicit antigen specific Th1 responses by promoting maturation of type 1 dendritic cells (DC1). AIMS Allergy and Immunology, 2017, 1(2): 89-100. doi: 10.3934/Allergy.2017.2.89
  • Tuberculosis has affected human beings for thousands of years, and until today, tuberculosis still ranks third among 29 infectious diseases in China. However, most of the existing mathematical models consider a single factor, which is not conducive to the study of tuberculosis transmission dynamics. Therefore, this study considers the combined effects of vaccination, treatment, and contaminated environments on tuberculosis, and builds a new model with seven compartments of SVEITRW based on China's tuberculosis data. The study shows that when the basic reproduction number R0 is less than 1, the disease will eventually disappear, but when R0 is greater than 1, the disease may persist. In the numerical analysis part, we use Markov-chain Monte-Carlo method to obtain the optimal parameters of the model. Through the next generation matrix theory, we calculate that the R0 value of tuberculosis in China is 2.1102, that is, if not controlled, tuberculosis in China will not disappear over time. At the same time, through partial rank correlation coefficients, we find the most sensitive parameter to the basic reproduction number R0. On this basis, we combine the actual prevalence of tuberculosis in China, apply Pontryagin's maximum principle, and perform cost-effectiveness analysis to obtain the conditions required for optimal control. The analysis shows that four control strategies could effectively reduce the prevalence of TB, and simultaneously controlling u2,u3,u4 is the most cost-effective control strategy.



    Atopic dermatitis (AD) affects nearly 31.6 million (10.1%) people in the United States [1] and of these, 18 million adults (7.2%) over the age of 18 have AD [2]. Gender disparities have been noted in the prevalence of AD with females being more prone during adulthood [3][12]. Data is lacking in males that supports mechanistic differences between the two sexes. Thus, many preclinical AD studies utilize female mice [13][14], with rare exceptions [14],[15].

    The pathophysiology of AD involves epidermal barrier dysfunction, a Th2/Th22 helper cell immune response and genetic polymorphisms which have a more uncertain role [16],[17]. Current therapies for AD are developed for lesional skin, although non-lesional skin is perturbed [17]. We previously reported pre-lesional skin remodeling and increased activation of mast cells (MC), skin-resident innate immune cells, via the signaling of a potently bioactive sphingolipid metabolite, sphingosine-1-phosphate (S1P), observed in female mice after 1 week of OVA treatment [18],[19]. This novel finding uncovered histologic and biochemical perturbations occurring in the skin before the development of any overt lesions. Here, we report pre-lesional skin responses in male mice upon a similar short 7-day epicutaneous exposure to OVA. We found that, conversely to female mice [18],[19], males did not exhibit any pre-lesional disease-inducing skin remodeling or MC activation. Importantly, samples from our previous report on female mice [18] and from the current report on male mice were collected and analyzed at the same time, following a gender-matched experimental design. While we did not ascertain the mechanism leading to the observed gender disparity, we established histological and biochemical differences in OVA-exposed male mouse skin compared to females that may partly elucidate the sex differences substantiated in adult AD prevalence.

    AD was induced in 8 to 12 weeks old male C57Bl/6J (WT) mice exactly as previously described and at the same time as the female mouse study [18],[20] (Charles River NCI, Frederick, MD, USA and The Jackson Laboratory Bar Harbor, ME, USA). Mouse randomization was performed using shuffled pieces of paper indicating treatment to assign mice to either experimental group. One by one centimeter square-sized gauze pad (“patch”) containing one hundred microliters of solution of either 100 µg OVA (Sigma-Aldrich, St Louis, MO, USA) or 0.9% saline vehicle only (saline control) was positioned on the shaved dorsal skin area after tape stripping. Patches were secured in position for seven days using a transparent dressing (Tegaderm, 3M HealthCare, St Paul, MN, USA) and bandages. Treated skin areas were harvested from euthanized mice. The AD model was validated as previously (our results [18] and [20]) through serum IgE level quantification in male mice that received three epicutaneous exposures, performed at the same time as in female mice. The full AD protocol consists in three rounds of one seven-day skin exposure via patch followed by fourteen days without patch. Circulating IgE levels were measured in sera from blood samples collected by cardiac puncture from euthanized mice and kept at −80 °C until use.

    All protocols were approved by the University of South Carolina Institutional Animal Care and Use Committee (Protocol 2403-101303-010818, “Mast cells and sphingosine-1-phosphate in allergic inflammation”). All experiments were performed in accordance with guidelines and regulations.

    Paraffin-embedded skin tissues were fixed with fresh paraformaldehyde (4%) and cut into 4 µm sections. Some sections were hematoxylin and eosin (H&E)-stained for morphometric analysis as previously described [18]. Others were stained with 0.1% methylene blue (MB) for 5 seconds, rinsed with water, dehydrated, and coverslip-mounted to locate MC and MC activation status [18],[21],[22]. All imaging and measurements were performed by blinded investigators. Imaging was performed at 10× and 40× magnifications using a Nikon E-600 microscope (Nikon Inc., Melville, NY) and a Micropublisher digital camera model 5.0. The QImaging software version 2.0.13 (QImaging Corp., Surrey, BC, Canada) was utilized for image analysis of 24-bit color or 8-bit monochrome images. The MetaMorph® 6.1 software (Molecular Devices, Sunnyvale, CA, USA) and the HarFa imaging software (www.fch.vutbr.cz/lectures/imagesci/) were employed to calculate morphometric and fractal parameters of nuclei allowing for infiltration quantification and measure MC status. Infiltration was measured in 25 to 30 ROIs per image using at least twelve 40× images per skin section and 4 skin sections per 2 mice from each treatment group. MC activation status (percent degranulation) was calculated in situ after analyzing fifty 40× images and 310 to 485 total MC present in 4 skin sections from 2 mice in each treatment group [18],[22].

    Epidermis, dermis, and hypodermis thickness was measured in saline- and OVA-treated samples with an ocular micrometer (Klarmann Rulings, Inc., Litchfield, NH, USA), exactly as described [18],[22]. The average thickness was calculated from 3 to 5 measurements (n = 3–5 mice per treatment group for each skin layer).

    We adapted our innovative imaging methods [22] to measure hypodermal cellular infiltration [18], using at least twelve high-magnification H&E images per individual mouse. First, the set-color-threshold subroutine (MetaMorph software) was used to isolate nuclei by color thresholding [18]. Next, the hypodermis in each section was surveyed using a fixed 75 µm diameter-sized region of interest (ROI) [18]. Thus, 25 to 30 ROIs were evaluated per skin section for a total of 2 sections per slide, using 2 slides for each individual mouse and 2 mice for each treatment group. Our previous work validated the morphometric parameters of area, perimeter, and integrated optical density (MetaMorph software's Integrated Morphometry Analysis subroutine) for accurate quantification of nuclei with elimination of nonspecific staining [18].

    After thresholding (Section 2.5), fractal dimensions (D) were measured after 8-bit grayscale image conversion exactly as previously described [18],[22].

    We established an imaging method to quantitatively described the cytoplasm of MC using fractal geometry [22], demonstrating that a fractal dimension (D) value of 1.378 ± 0.062 specifically designated intact/resting MC, whereas a D value of 1.484 ± 0.048 characterized degranulated/activated MC [22]. Here, we used these D values for nonsubjective measure of MC activation in male skin tissue samples (n = 2 to 5 mice per treatment group).

    Harvested skin tissue samples were immediately snap-frozen, and kept at −80 °C until use. Total RNA was extracted from 4 to 10 mice per treatment group using the miRNeasy kit (Qiagen, Valencia, CA, USA) and cDNA synthesized by reverse transcription (iScript cDNA synthesis kit, Bio-Rad, Hercules, CA, USA). A CFX Connect (Bio-Rad) was used to perform qPCR (duplicate determinations), using SensiFAST™ SYBR No-ROX Kit (Bioline, Taunton, MA, USA). We used the following primers (all purchased from Thermo Fisher Scientific, Inc., Waltham, MA, USA). for qPCR amplification: GAPDH forward primer CAGAAGGGGCGGAGATGAT and reverse primer AGGCCGGTGCTGCTGAGTATGTC, FcϵRIα forward primer ATTGTGAGTGCCACCGTTCA and reverse primer GCAGCCAATCTTGCGTTACA, CCL2 forward primer CACTCACCTGCTGCTACTCA and reverse primer GCTTGGTGACAAAAACTACAGC, CCL3 forward primer GCCATATGGAGCTGACACCC and reverse primer TAGTCAGGAAAATGACACCTGGC, and CCL5 forward primer TGCCCTCACCATCATCCTCACT and reverse primer GGCGGTTCCTTCGAGTGACA. We used the following qPCR conditions: initial step at 95 °C for 5 minutes and 40 cycles, as follows: 10 seconds at 95 °C, followed by annealing for 1 minute at 55 °C and extension for 1 minute at 72 °C. Data were analyzed with CFX Manager™ Software, after normalization to control (saline-treated) samples and to the amount of reference gene, GAPDH mRNA levels.

    CCL2, CCL3 and CCL5 chemokine quantifications were conducted in mouse serum and whole skin protein extracts prepared exactly as previously described [18], using a Bio-Plex Array Reader (LUMINEX 100; Bio-Rad Laboratories, Hercules, CA, USA) and Milliplex panels (EMD Millipore, Billerica, MA, USA). Whole skin protein extracts were weighed and immediately snap-frozen using liquid nitrogen. After homogenization with a mortar and pestle, skin samples were digested for 2 hours at 4 °C in reporter lysis buffer (Promega, Madison, WI, USA) supplemented with complete mini-protease inhibitors (Roche Diagnostics, Indianapolis, IN, USA). Supernatants were kept at −80 °C until use.

    Total IgE was quantified in mouse serum by ELISA (R&D Systems, Minneapolis, MN, USA). N = 6 and 9 mice per experimental group for 1 and 3 exposures, respectively.

    Lipids were extracted from snap-frozen skin tissues and S1P or sphingosine was quantified by liquid chromatography-electrospray ionization-tandem mass spectrometry (4000 QTRAP; AB Sciex, Foster City, CA, USA), as before [18], and normalized to skin sample weights. Saline and OVA solutions contained no detectable levels of S1P or sphingosine (data not shown). N = 3 to 4 samples per experimental group.

    Expressed as means ± SEM (unless otherwise indicated), data was analyzed using the unpaired 2-tailed Student's t-test, with Welch's correction for samples of unequal variance (Prism 6; GraphPad Software, La Jolla, CA, USA). Significance for all statistical tests is provided in figures. Each experiment was repeated at least 3 times in triplicate determinations with consistent results.

    We previously demonstrated that the epidermis and dermis were thickened in female mice after a single dose of OVA applied to the skin [18]. Moreover, there was an increase in cellular infiltration particularly around blood vessels (BV) [18]. Representative images of H&E-stained saline (Figure 1a) and OVA (Figure 1b) exposed male mice skins illustrate that no significant change was evidenced in epidermal, dermal or hypodermal thickness after 1 week of OVA treatment compared to controls (Figure 1c), in contrast to our reported findings in female mice collected at the same time [18]. We then measured cellular infiltration in the hypodermis by color thresholding the nuclei in the entire hypodermis of H&E-stained skin sections, exactly as previously described [18]. We report increased cellular infiltration in the hypodermis and around BV (Figure 1d–e) after OVA exposure of male dorsal skins, compared to saline controls. In order to independently confirm these results, we applied yet another of our previously published computerized imaging modalities to quantify cellular mass, a proxy for cell numbers/infiltration, by means of fractals, or space filling capacity [18],[22]. Cellular infiltration of the hypodermis (Figure 1d–e) was confirmed by an increase in cellular mass, expressed as augmented fractal dimension (D) values in the hypodermis and around BV (Figure 1f) after OVA treatment and compared to saline controls.

    Figure 1.  Skin remodeling quantitation in differentially treated male skin samples. Hypodermal cellular infiltration does not cause skin thickening after OVA exposure. Representative images of H&E-stained skin sections of (a) control (b) OVA treated mice. Original magnification 10× (scale bar = 50 µm). (c) Epidermal, dermal and hypodermal thickness. (d) Representative nuclei quantification in each region of interest across the hypodermis of one saline- (empty circles) and one OVA-treated (filled circles) mice, also distinguishing blood vessel (BV)-containing ROI (red symbols). (e) Quantification of cellular infiltration (nuclei per ROI) after saline (black symbols) or OVA (red symbols) skin treatment. (f) Validation of cell quantification by fractal dimension (D) of each ROI. N = 10–20 individual measurements collected from 3 to 5 mice per treatment group and per skin layer. ROI determinations: 25 to 30 ROI per image, at least 12 40×-images per skin section, 4 skin sections, 2 mice per treatment group.

    We previously reported that mRNA levels coding for FcϵRIα, the IgE-binding sub-unit of the high-affinity receptor for IgE reflected local MC numbers in this animal model of AD [18], yet were not significantly increased after OVA treatment, neither was MC degranulation (Figure 2b,d) compared to saline controls (Figure 2a,c) and in contrast to female mice [18].

    Figure 2.  Mast cell activation status in differentially treated male skin samples. OVA exposure does not cause skin MC activation after one OVA exposure in male mice. FcϵRIα (a) and percentage of degranulated MC (b) after OVA treatment (grey columns) compared to saline controls (black columns). Representative images of methylene blue-stained skin sections collected after saline (c) or OVA (d) exposure. Arrows indicate intact (red) and degranulated (yellow) MC, as detected by a validated computer-aided image analysis method. The insets represent magnifications of resting (red frame) or a degranulated (yellow frame) MC at 40×. Bar = 50 µm. Unpaired 2-tailed Student's t-test with Welch's correction. NS, not significant.

    As expected, the serum levels of IgE were not increased after one single OVA exposure but were increased after 3 OVA exposures (Figure 3a). Thus, we previously established that skin MC degranulation at this early stage of AD was not due to the classical IgE/Ag activation pathway but rather depended on S1P signaling in MC [18],[19]. Unaltered MC degranulation in male skin was concurrent to unchanged local levels of sphingolipid metabolites sphingosine (Figure 3b) and S1P (Figure 3c) and Sphk1 mRNA after OVA treatment compared to saline controls (Figure 3d), confirming the newly identified mechanistic insights we reported that pertained to AD onset [18],[19].

    Figure 3.  Circulating IgE levels, skin-associated sphingolipid sphingosine and S1P and SphK1 mRNA levels in differentially treated male skin samples. One OVA exposure does not elevate skin sphingosine or S1P levels, SphK1 mRNA expression or circulating IgE levels in male mice. (a) Serum IgE levels after 3 OVA or saline treatments. Skin levels of sphingosine (b), S1P (c) and SphK1 mRNA expression (d). Unpaired 2-tailed Student's t-test with Welch's correction. NS, not significant.

    The elevation of local cytokines/chemokines drives the development of AD, whereas their systemic augmentation correlates with disease progression and/or severity [17]. We previously demonstrated that airway [23] and skin [18] inflammatory responses were linked to local rather than systemic elevation of chemokines, thus demonstrating their pathogenic function in atopy. Whereas local levels of CCL2 (Figure 4a,b), CCL3 (Figure 4a,c) and CCL5 (Figure 4a,d) chemokines did not augment upon OVA skin treatment, circulating CCL2, CCL3 and CCL5 were significantly elevated in OVA-exposed male mice skin, compared to saline controls (Figure 4e).

    Figure 4.  Skin-associated versus circulating levels of chemokines in differentially treated male mice. Systemic but not local chemokines are elevated after one OVA exposure in male mice. (a) Skin CCL2, CCL3 and CCL5 mRNA expression. Skin protein levels of CCL2 (b), CCL3 (c) and CCL5 (d) and their circulating levels (e) in OVA- or saline-treated mice. Unpaired 2-tailed Student's t-test with Welch's correction. NS, not significant.

    Female mice are predominantly used in established preclinical models of AD that encompass epicutaneous applications of OVA [20], partly because they feature stronger antigen-dependent cutaneous inflammatory responses than males, similar to the female preponderance of AD observed in adolescent and adult humans [12],[24][26]. Although no single murine model fully recapitulates all dysregulated immune and barrier functions of human AD skin, OVA-challenged mice replicate activation of inflammatory pathways observed in patients with AD. Moreover, OVA as an AD-inducing agent can be used on different mouse backgrounds to account for genetic variability allowing to identify key players of cutaneous inflammation [14]. In gender-matched experiments, we noticed that male mice do not feature the same histopathological patterns seen in female mice. Pre-lesional skin changes occurring in males after one OVA exposure include hypodermal cell infiltration without epidermal or dermal thickening, suggesting a difference in the nature or the distribution of infiltrating cells between the two sexes. One limitation of our study is that we did not identify the cell types that were increased in the hypodermis. However, reports in mice identified eosinophils, MC and CD3+CD4+ cells [27] and human skin samples featured CD3+ T-cells and CD11c+ dendritic cells in non-lesional skin from patients with AD [28]. Skin changes occurring over time featured gradual chemokine increase, including CCL3 in female Balb/c mice after OVA treatment [27], in agreement with our data [18]. The paucity of reports analyzing the sequence of skin alterations at the onset of AD mouse models in males is noteworthy. Local CCL2, CCL3 and CCL5 mRNA and CCL2 and CCL3 protein levels were elevated in female mice [18] but not in male mice (our current study). These chemokines have been shown to induce recruitment of macrophages, granulocytes, CD3+ and CD4+ cells [23],[29]. Importantly, MC secrete these chemokines [30] and our group previously reported that S1P-mediated MC activation is essential to their local production in early allergic responses of the lung [23] and the skin [18]. In humans, data is mainly collected using skin samples from patients affected with moderate-to-severe AD, with overt skin lesions. Interestingly, many inflammatory chemokines, including CCL2, CCL3 and CCL5 are associated with human chronic, supporting leukocyte recruitment [31]. More recently, seminal studies focusing on the analysis of nonlesional skin specimens collected from patients with lesional AD, also demonstrated an increased inflammatory T cell infiltration, with elevated CCL5, CCL11, CCL17, CCL18 and CCL22 [28]. However, the authors acknowledged that their dataset does not reveal the key chemokines involved in primary pathogenesis of AD, which our studies are focusing on. Using immunohistochemistry, Su et al. established the chemokine profiles to study the evolution of T cell subsets in skin biopsies of 48 patients and their association with the progression of human AD from acute to chronic disease [32]. CCL5 was associated with sub-acute AD. Our current study established that male skin MC are not activated after a single OVA exposure, further supporting that MC activation does promote pre-lesional skin remodeling through local chemokine elevation. It has been postulated that increased local chemokine levels may result in systemic dissemination and elevated serum levels [33]. However, male mice did not display any elevation of local chemokines despite augmented serum levels. Brunner et al. suggested the occurrence of a mutual regulation with correlation of blood and skin markers in lesional and nonlesional human AD skin [34]. In this study, chemokines shared concomitant regulation between chronic lesional skin and serum in patients with moderate-to-severe AD. The authors found a strong systemic inflammation with elevated skin chemokine levels in this cohort of overweight patients, a plausible explanation for a pre-existing systemic inflammatory state. Because AD leads to systemic inflammation, Lauffer et al. proposed to use serum levels of chemokines, including CCL2, CCL3 and CCL5, as possible predictors of AD persistence in children [35]. However, no data was provided pertaining to skin-associated chemokine levels. Our data suggests that systemic detection of chemokines may perhaps mark resistance to AD development whereas local elevation may denote AD onset, susceptibility or a sign of active disease. Little is known about the dynamics of chemokine production in serum and in skin tissue or about the relative contributions of skin resident cells to inflammation development during the pre-lesional phase of AD, which both warrant further investigation.

    It is noteworthy that serum IgE levels are elevated in ~80% of patients with AD however, IgE levels do not necessarily correlate with disease severity [16]. In male mice, IgE levels were significantly elevated after 3 OVA patches, thus validating this AD model for both genders, although three skin exposures to OVA induced 2.5 times more serum IgE in female [18] compared to male mice (Figure 3a). Similar to our findings in female mice and as expected, a single OVA skin exposure did not suffice to elevate serum IgE levels in male mice (Figure 3a).

    Previous reports have demonstrated that epidermal thickening is observed in AD lesions with or without the presence of MC in filaggrin-deficient mice, rendering MC dispensable in this genetic model [36]. Notably, the disruption of the epidermal barrier due to filaggrin deficiency preceded any local MC activation in this model. It is tempting to speculate that MC may further perturb the skin barrier to enable AD. In contrast, other reports concluded that MC were protective to the epidermal barrier, also suggesting that MC mediators regulated skin barrier integrity [37]. However, the mouse gender was not revealed in their study where transgenic rather than allergen-induced AD models were employed, perhaps contributing to MC function discrepancies. Indeed, AD lesions feature neo-angiogenesis with MC-associated vascular endothelial growth factor-A (VEGF) [38]. In human AD lesions, epidermal thickening is observed, independently of MC localization to the epidermis [39]. In our current study, increased cellular infiltrate in the hypodermis and around BV was seen, as in female mice [18], but the absence of male skin remodeling was linked to lack of MC activation and of local CCL2, CCL3 and CCL5 mRNA and proteins.

    We previously identified a novel signaling pathway in MC by which SphK1-produced S1P from substrate sphingosine binds S1P receptor 2 (S1PR2) [23],[39][43], signaling through signal transducer and activator of transcription 3 to release CCL2, CCL3, CCL5, resulting in T-cell recruitment [23]. S1P/S1PR2 signaling in MC also triggered VEGF and matrix metalloproteinase 2 release, further pointing to an important role for MC in pre-lesional AD skin remodeling [44]. However, the MC/S1P axis was not disease inducing in male mice using the OVA-induced AD model, as skin levels of sphingosine, S1P or mRNA coding for SphK1 were not elevated. Moreover, MC degranulation was not augmented after OVA compared to control male mice, further supporting a critical role for S1P-mediated MC activation in pre-lesional AD skin remodeling, as observed in female mice. Nonetheless, OVA exposure did increase hypodermal cellular infiltration and systemic chemokine levels in male mice and three OVA exposures resulted in elevated serum IgE. Together, our results suggest that male mice may undergo an attenuated and/or delayed skin inflammation due to the absence of MC activation. Indeed, previous studies suggest that female sex hormones may exacerbate Th2 inflammation. Subcutaneous estradiol injection initiates allergic sensitization [45] and estradiol promotes in vitro MC degranulation [46],[47]. Anaphylactic responses are also more severe in female than male mice [48].

    There are potential limitations to this study. First, we are well aware that the mechanisms that underlie sex-related differences in skin disease and inflammation, while, mostly unknown, have not been addressed here. In addition to sex hormones as discussed above, genetic background and differences in skin barrier properties may also contribute to these differences. A better understanding of sex-related differences in skin inflammation may help better prevent, diagnose, and treat inflammatory skin diseases. Second, while we previously demonstrated an essential role for SphK1 in the upregulation of skin S1P levels after stimulation [18], SphK1 was similarly yet only assessed transcriptionally. Nonetheless, the unaltered levels of skin S1P after stimulation strongly supported an essential role for SphK1 in S1P elevation over any other isoform of SphK. Collectively, the current study demonstrates contrasted pathogenic mechanisms enabling skin inflammation in preclinical AD in male mice, compared to reports using female mice [18],[19].

    Clinical evidence had demonstrated that anti-IgE biologics do not affect the course of AD [16],[49], but may have therapeutic potential at high doses [50]. Thus, identifying S1P as a key MC stimulus of pre-lesional skin remodeling may be an important milestone in our understanding of factors that govern the transition from healthy to AD-prone skin. Estrogen exerts immunomodulatory functions on IgE-independent MC degranulation and chemokine production and higher susceptibility of human and mouse females to multiple sclerosis is due to higher production of S1PR2 compared to males [10], further supporting important inflammatory disease enabling functions of the MC/S1P axis.

    Our data suggests that S1P-mediated activation of MC is a critical event in early-stage AD, absent in male mice, perhaps relevant to higher prevalence of AD in females, supporting MC intervention to regulate the magnitude and/or perhaps time course [51] of inflammatory skin responses in a gender-dependent manner.



    [1] Wikipedia (2023) Tuberculosis, 2023. Available from: https://en.wikipedia.org/wiki/Tuberculosis.
    [2] J. Grange, M. Gandy, P. Farmer, A. Zumla, Historical declines in tuberculosis: nature, nurture and the biosocial model [Counterpoint], Int. J. Tuberc. Lung D, 5 (2001), 208–212. https://doi.org/10.1067/mhl.2001.118364 doi: 10.1067/mhl.2001.118364
    [3] World Health Organization, 2023. Available from: https://www.who.int/en/news-room/fact-sheets/detail/tuberculosis.
    [4] H. Waaler, A. Geser, S. Andersen, The use of mathematical models in the study of the epidemiology of tuberculosis, Am. J. Public Health, 52 (1962), 1002–1013. https://doi.org/10.2105/AJPH.52.6.1002 doi: 10.2105/AJPH.52.6.1002
    [5] C. Dye, G. P. Garnett, K. Sleeman, B. G. Williams, Prospects for worldwide tuberculosis control under the WHO DOTS strategy, Lancet, 352 (1998), 1886–1891. https://doi.org/10.1016/s0140-6736(98)03199-7 doi: 10.1016/s0140-6736(98)03199-7
    [6] B. Song, C. Castillo-Chavez, J. P. Aparicio, Tuberculosis models with fast and slow dynamics: the role of close and casual contacts, Math. Biosci., 180 (2002), 187–205. https://doi.org/10.1016/s0025-5564(02)00112-8 doi: 10.1016/s0025-5564(02)00112-8
    [7] T. C. Porco, S. M. Blower, Quantifying the intrinsic transmission dynamics of tuberculosis, Theor. Popul. Biol., 54 (1998), 117–132. http://dx.doi.org/10.1006%2Ftpbi.1998.1366
    [8] C. Bhunu, W. Garira, Z. Mukandavire, G. Magombedze, Modelling the effects of pre-exposure and post-exposure vaccines in tuberculosis control, J. Theor. Biol., 254 (2008), 633–649. https://doi.org/10.1016/j.jtbi.2008.06.023 doi: 10.1016/j.jtbi.2008.06.023
    [9] N. Mantilla, M. Gomes, Mycobacterial ecology as a modulator of tuberculosis vaccine success, Theor. Popul. Biol., 75 (2009), 142–152. http://dx.doi.org/10.1016/j.tpb.2009.01.006 doi: 10.1016/j.tpb.2009.01.006
    [10] Y. Yang, S. Tang, X. Ren, H. Zhao, C. Guo, Global stability and optimal control for a tuberculosis model with vaccination and treatment, Discrete Cnotin. Dyn. B, 21 (2016), 1009–1022. https://doi.org/10.3934/dcdsb.2016.21.1009 doi: 10.3934/dcdsb.2016.21.1009
    [11] E. Ziv, C. L. Daley, S. Blower, Potential public health impact of new tuberculosis vaccines, Emerg. Infect. Dis., 10 (2004), 1529–1535. https://doi.org/10.3201/eid1009.030921 doi: 10.3201/eid1009.030921
    [12] S. Ren, Global stability in a tuberculosis model of imperfect treatment with age-dependent latency and relapse, Math. Biosci. Eng., 14 (2017), 1337–1360. http://dx.doi.org/10.3934/mbe.2017069 doi: 10.3934/mbe.2017069
    [13] J. M. Trauer, J. T. Denholm, E. S. McBryde, Construction of a mathematical model for tuberculosis transmission in highly endemic regions of the Asia-Pacific, J. Theor. Biol., 358 (2014), 74–84. https://doi.org/10.1016/j.jtbi.2014.05.023 doi: 10.1016/j.jtbi.2014.05.023
    [14] C. Bhunu, W. Garira, Z. Mukandavire, Modeling HIV/AIDS and tuberculosis coinfection, Bull. Math. Biol., 71 (2009), 1745–1780. https://doi.org/10.1007/s11538-009-9423-9 doi: 10.1007/s11538-009-9423-9
    [15] Z. Guo, H. Xiang, H. Huo, Analysis of an age-structured tuberculosis model with treatment and relapse, J. Math. Biol., 82 (2021), 1–37. https://doi.org/10.1007/s00285-021-01595-1 doi: 10.1007/s00285-021-01595-1
    [16] L. Xue, S. Jing, H. Wang, Evaluating strategies for tuberculosis to achieve the goals of WHO in China: a seasonal age-structured model study, Bull. Math. Biol., 84 (2022), 61. http://dx.doi.org/10.1007/s11538-022-01019-1 doi: 10.1007/s11538-022-01019-1
    [17] S. Bowong, J. Kurths, Modeling and analysis of the transmission dynamics of tuberculosis without and with seasonality, Nonlinear Dyn., 67 (2012), 2027–2051. http://dx.doi.org/10.1007%2Fs11071-011-0127-y
    [18] T. Li, Y. Guo, Modeling and optimal control of mutated COVID-19 (Delta strain) with imperfect vaccination, Chaos Soliton Fractals, 156 (2022), 111825. https://doi.org/10.1016/j.chaos.2022.111825 doi: 10.1016/j.chaos.2022.111825
    [19] Y. Guo, T. Li, Modeling the competitive transmission of the Omicron strain and Delta strain of COVID-19, J. Math. Anal. Appl., 526 (2023), 127283. https://doi.org/10.1016/j.jmaa.2023.127283 doi: 10.1016/j.jmaa.2023.127283
    [20] T. R. Hawn, T. A. Day, T. J. Scriba, M. Hatherill, W. A. Hanekom, T. G. Evans, et al., Tuberculosis vaccines and prevention of infection, Microbiol. Mol. Biol. R., 78 (2014), 650–671. https://doi.org/10.1128/mmbr.00021-14
    [21] G. Disease, I. Incidence, L. Monasta, L. Ronfani, E. Beghi, B. Giussani, et al., Global, regional, and national incidence, prevalence, andyears lived with disability for 354 diseases and injuries for 195 countries and territories, 1990–2017: a systematic analysis for the Global Burden of Disease Study 2017, Lancet, 392 (2018), 1789–1858. https://doi.org/10.1016/S0140-6736(18)32279-7 doi: 10.1016/S0140-6736(18)32279-7
    [22] P. J. White, G. P. Garnett, Mathematical modelling of the epidemiology of tuberculosis, Adv. Exp. Med. Biol., 2010 (2010), 127–140. https://doi.org/10.1007/978-1-4419-6064-1_9 doi: 10.1007/978-1-4419-6064-1_9
    [23] Y. Li, X. Liu, Y. Yuan, J. Li, L. Wang, Global analysis of tuberculosis dynamical model and optimal control strategies based on case data in the United States, Appl. Math. Comput., 422 (2022), 126983. https://doi.org/10.1016/j.amc.2022.126983 doi: 10.1016/j.amc.2022.126983
    [24] Q. Jiang, Z. Liu, L. Wang, R. Tan, A tuberculosis model with early and late latency, imperfect vaccination, and relapse: An application to China, Math. Method. Appl. Sci., 46 (2023), 10929–10946. https://doi.org/10.1002/mma.9160 doi: 10.1002/mma.9160
    [25] A. Xu, Z. Wen, Y. Wang, W. Wang, Prediction of different interventions on the burden of drug-resistant tuberculosis in China: a dynamic modelling study, J. Glob. Antimicrob. Res., 29 (2022), 323–330. https://doi.org/10.21203/rs.3.rs-637762/v1 doi: 10.21203/rs.3.rs-637762/v1
    [26] Y. Cai, S. Zhao, Y. Niu, Z. Peng, K. Wang, D. He, et al., Modelling the effects of the contaminated environments on tuberculosis in Jiangsu, China, J. Math. Biol., 508 (2021) 110453. https://doi.org/10.1016/j.jtbi.2020.110453
    [27] P. Van den Driessche, J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci., 180 (2002), 29–48. https://doi.org/10.1016/s0025-5564(02)00108-6 doi: 10.1016/s0025-5564(02)00108-6
    [28] J. LaSalle, The Stability of Dynamical Systems, Society for Industrial and Applied Mathematics, Philadelphia Pennsylvania, 1976. https://doi.org/10.1137/1.9781611970432
    [29] X. Zhao, Dissipative Dynamical Systems, Dynamical Systems in Population Biology, Spring Link, New York, 2003. https://doi.org/10.1007/978-3-319-56433-3
    [30] X. Zhao, Uniform persistence and periodic coexistence states in infinite-dimensional periodic semiflows with applications, Canad. Appl. Math. Quart., 3 (1995), 473–495.
    [31] W. Wang, X. Zhao, An epidemic model in a patchy environment, Math. Biosci., 190 (2004), 97–112. https://doi.org/10.1016/j.mbs.2002.11.001 doi: 10.1016/j.mbs.2002.11.001
    [32] Y. Guo, T. Li, Fractional-order modeling and optimal control of a new online game addiction model based on real data, Commun. Nonlinear Sci., 121 (2023), 107221. https://doi.org/10.1016/j.cnsns.2023.107221 doi: 10.1016/j.cnsns.2023.107221
    [33] L. Pontryagin, V. Boltyanskii, R. Gamkrelidze, E. Mischenko, The Mathematical Theory of Optimal Processes, Gordon and Breach Science Publishers, New York, 1986. https://doi.org/10.1201/9780203749319
    [34] China Center for Disease Control and Prevention, 2021. Available from: https://www.ndcpa.gov.cn/jbkzzx/c100016/common/list.html.
    [35] China Population Statistic Yearbook, 2020. Available from: https://www.stats.gov.cn/sj/ndsj/2022/indexch.htm.
    [36] S. Liu, Y. Bi, Y. Liu, Modeling and dynamic analysis of tuberculosis in mainland China from 1998 to 2017: the effect of DOTS strategy and further control, Theor. Biol. Med. Model., 17 (2020), 6. https://doi.org/10.1186/s12976-020-00124-9 doi: 10.1186/s12976-020-00124-9
    [37] L. Xue, X. Ren, W. Sun, X. Zheng, Z. Peng, B. Singh, Seasonal transmission dynamics and optimal control strategies for tuberculosis in Jiangsu Province, China, Math. Method. Appl. Sci., 46 (2023), 2072–2092. https://doi.org/10.1002/mma.8629 doi: 10.1002/mma.8629
    [38] S. Marino, I. B. Hogue, C. J. Ray, D. E. Kirschner, A methodology for performing global uncertainty and sensitivity analysis in systems biology, J. Theor. Biol., 254 (2008), 178–196. https://doi.org/10.1016/j.jtbi.2008.04.011 doi: 10.1016/j.jtbi.2008.04.011
    [39] H. Pan, S. Bele, Y. Feng, S. Qiu, J. Lü, S. Tang, et al., Analysis of the economic burden of diagnosis and treatment of tuberculosis patients in rural China, Int. J. Tuberc. Lung D., 17 (2013), 1575–1580. https://doi.org/10.5588/ijtld.13.0144 doi: 10.5588/ijtld.13.0144
    [40] X. Li, Y. Ruan, X. Liu, C. Xu, W. Chen, H. D. Xin, Analysis on funding input and benefit output for tuberculosis control and prevention in China from 2011 to 2019, Chin. J. Antituberc., 43 (2021), 702–707. https://doi.org/10.3969/j.issn.1000-6621.2021.07.011 doi: 10.3969/j.issn.1000-6621.2021.07.011
  • This article has been cited by:

    1. Christopher D. Carlucci, Yvonne Hui, Alena P. Chumanevich, Piper A. Robida, John W. Fuseler, Mathew Sajish, Prakash Nagarkatti, Mitzi Nagarkatti, Carole A. Oskeritzian, Resveratrol Protects against Skin Inflammation through Inhibition of Mast Cell, Sphingosine Kinase-1, Stat3 and NF-κB p65 Signaling Activation in Mice, 2023, 24, 1422-0067, 6707, 10.3390/ijms24076707
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1670) PDF downloads(238) Cited by(4)

Figures and Tables

Figures(8)  /  Tables(4)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog