Research article Special Issues

A nonstandard finite difference scheme for a time-fractional model of Zika virus transmission


  • Received: 05 November 2023 Revised: 01 December 2023 Accepted: 05 December 2023 Published: 21 December 2023
  • In this work, we investigate the transmission dynamics of the Zika virus, considering both a compartmental model involving humans and mosquitoes and an extended model that introduces a non-human primate (monkey) as a second reservoir host. The novelty of our approach lies in the later generalization of the model using a fractional time derivative. The significance of this study is underscored by its contribution to understanding the complex dynamics of Zika virus transmission. Unlike previous studies, we incorporate a non-human primate reservoir host into the model, providing a more comprehensive representation of the disease spread. Our results reveal the importance of utilizing a nonstandard finite difference (NSFD) scheme to simulate the disease's dynamics accurately. This NSFD scheme ensures the positivity of the solution and captures the correct asymptotic behavior, addressing a crucial limitation of standard solvers like the Runge-Kutta Fehlberg method (ode45). The numerical simulations vividly demonstrate the advantages of our approach, particularly in terms of positivity preservation, offering a more reliable depiction of Zika virus transmission dynamics. From these findings, we draw the conclusion that considering a non-human primate reservoir host and employing an NSFD scheme significantly enhances the accuracy and reliability of modeling Zika virus transmission. Researchers and policymakers can use these insights to develop more effective strategies for disease control and prevention.

    Citation: Maghnia Hamou Maamar, Matthias Ehrhardt, Louiza Tabharit. A nonstandard finite difference scheme for a time-fractional model of Zika virus transmission[J]. Mathematical Biosciences and Engineering, 2024, 21(1): 924-962. doi: 10.3934/mbe.2024039

    Related Papers:

    [1] António M. Cunha, Ricardo Loureiro . Housing price dynamics and elasticities: Portugal's conundrum. National Accounting Review, 2024, 6(1): 75-94. doi: 10.3934/NAR.2024004
    [2] David Gray . Medium-term cycles in affordability: what does the house price to income ratio indicate?. National Accounting Review, 2021, 3(2): 204-217. doi: 10.3934/NAR.2021010
    [3] Alejandro Almeida, Julia Feria, Antonio Golpe, José Carlos Vides . Financial assets against inflation: Capturing the hedging properties of gold, housing prices, and equities. National Accounting Review, 2024, 6(3): 314-332. doi: 10.3934/NAR.2024014
    [4] Nuraddeen Umar Sambo, Ibrahim Sambo Farouq, Mukhtar Tijjani Isma'il . Asymmetric effect of exchange rate volatility on trade balance in Nigeria. National Accounting Review, 2021, 3(3): 342-359. doi: 10.3934/NAR.2021018
    [5] Jawad Saleemi . In COVID-19 outbreak, correlating the cost-based market liquidity risk to microblogging sentiment indicators. National Accounting Review, 2020, 2(3): 249-262. doi: 10.3934/NAR.2020014
    [6] Luccas Assis Attílio . Impact of Chinese financial shocks: A GVAR approach. National Accounting Review, 2024, 6(1): 27-49. doi: 10.3934/NAR.2024002
    [7] Irena Honková, Renáta Myšková . Techniques, objectives and motivations of creative accounting: evidence from the Czech Republic. National Accounting Review, 2024, 6(3): 333-351. doi: 10.3934/NAR.2024015
    [8] Jinhui Zhu, Zhenghui Li . Can digital financial inclusion effectively stimulate technological Innovation of agricultural enterprises?—A case study on China. National Accounting Review, 2021, 3(4): 398-421. doi: 10.3934/NAR.2021021
    [9] Mohammad Reza Abbaszadeh, Mehdi Jabbari Nooghabi, Mohammad Mahdi Rounaghi . Using Lyapunov's method for analysing of chaotic behaviour on financial time series data: a case study on Tehran stock exchange. National Accounting Review, 2020, 2(3): 297-308. doi: 10.3934/NAR.2020017
    [10] Changqing Luo, Zijing Li, Lan Liu . Does investor sentiment affect stock pricing? Evidence from seasoned equity offerings in China. National Accounting Review, 2021, 3(1): 115-136. doi: 10.3934/NAR.2021006
  • In this work, we investigate the transmission dynamics of the Zika virus, considering both a compartmental model involving humans and mosquitoes and an extended model that introduces a non-human primate (monkey) as a second reservoir host. The novelty of our approach lies in the later generalization of the model using a fractional time derivative. The significance of this study is underscored by its contribution to understanding the complex dynamics of Zika virus transmission. Unlike previous studies, we incorporate a non-human primate reservoir host into the model, providing a more comprehensive representation of the disease spread. Our results reveal the importance of utilizing a nonstandard finite difference (NSFD) scheme to simulate the disease's dynamics accurately. This NSFD scheme ensures the positivity of the solution and captures the correct asymptotic behavior, addressing a crucial limitation of standard solvers like the Runge-Kutta Fehlberg method (ode45). The numerical simulations vividly demonstrate the advantages of our approach, particularly in terms of positivity preservation, offering a more reliable depiction of Zika virus transmission dynamics. From these findings, we draw the conclusion that considering a non-human primate reservoir host and employing an NSFD scheme significantly enhances the accuracy and reliability of modeling Zika virus transmission. Researchers and policymakers can use these insights to develop more effective strategies for disease control and prevention.



    The primary sensory neurons that transmit information from muscle spindles in skeletal muscles of the limbs and trunk to spinal cord are clustered together in the dorsal root ganglion (DRG) located adjacent to spinal cord. The DRGs contain sensory neurons with central and peripheral branches projecting towards either terminate within spinal gray matter or ascends through gracile or cuneate funiculus along dorsal column and ends at the junctions of spinal cord with the medulla. The peripheral branches connect to motor neurons located in ventral grey matter. The central branches from hind limbs run a long way along gracile fasciculus and terminate in gracile nucleus, while those from fore limbs and upper trunk run in cuneate fasciculus and terminate in cuneate nucleus. The gracile and cuneate sensory tract is bundled up travel side by side within dorsal column of spinal cord (see a schematic view in Figure 1a) [1]. This review describes a mutant mouse showing abnormal behaviors due to severe sensory and motor neurological disorders.

    In January 1984, abnormally suffering gaits were found in three 80-day-old mice, two females and one male, in the Laboratory of Genetics, Nagoya University [2]. The clinical and pathological examinations were made promptly in the National Institute of Neuroscience in NCNP, Tokyo. The homozygous mice were not identified from normal controls until 30-day-old. The first sign of abnormal behavior was confirmed when they were suspended by tail. Although normal mice extend their limbs out from the trunk, the hind limbs of homozygotes are initially stretched spasmodically and then become flexed and folded tightly towards the trunk (sensory ataxia stage), By about 80-day-old, the hindlimb weakness and tremor, muscle atrophy became apparent (motor paresis stage). The neurological signs were progressively severe, and most of death occur by about 5 to 6 months of age. The inheritance followed autosomal recessive mode. Linkage tests showed that the genetic mutation located on chromosome 5 linking to Pgm-1 and W [3].

    Figure 1.  (a) Schematic view of nervous pathways shown on a series of spinal cord slices. The DRG axons have two branches, one projecting to the PNS (muscle spindle), and one projecting to the CNS (gracile nucleus). Axonal degeneration (*) begins retrogradely from axonal terminals of both CNS and PNS sides (arrows) along proprioceptive sensory neurons to parent cell bodies in DRG [1]. (b) Eosinophilic spheroids (arrowheads) in gracile nucleus at 78-day old gad mouse. H-E staining. Bar = 100 µm [4]. (c) Annulo-spiral structures at the sensory endings (Ia and II) of a muscle spindle are dissolved and their distal tips show partial swelling in EDL muscle of gad mouse aged 30-days. AchE-silver staining. Bar = 40 µm [5].

    This mutant mouse was named originally “gracile axonal dystrophy (GAD)”, but, since etiology became clear, a gene symbol (gad) came to be used widely. The earliest pathological sign was detected as early as 40-day-old in gracile nucleus at the junction of spinal cord with medulla oblongata. The swollen eosinophilic axons called “spheroid” scatter over this nucleus (Figure 1b). The spheroids are also called “dystrophic axons” showing enormous enlargement most striking in the terminal segment of affected axons that appear just distal to the normal termination of the myelin sheath [6],[7]. They are often filled a number of neurofilaments, dense bodies, mitochondria with dark matrix, and multiple-vesicles which often surround a whirlpool towards the center. Abundant astrocytic processes are often seen around spheroids in lesioned sites [4],[8].

    The ultrastructural observations revealed the first signs of spheroids in gad mice at 10 to 20-day-old after birth. Spheroids varying in size are identified by slight presynaptic boutons which are facing to secondary neurons contact to and ascend to thalamus, or by thin myelin sheath surrounding them. They maintain a synapse-complexes containing electron dense bodies, mitochondria at their axonal ends. With progress of axonal degeneration, the cytoplasmic neurofilaments and thicker filament-like structures with 40 nm gathered intensely, moreover, few lysosomes contained amorphous substances seen in some places which make electron density in cytoplasm darker than usual, and finally loose presynaptic boutons. The filament-like structures gather keeping in line complicatedly and rarely branched off [4].

    The anterograde transports from DRGs are often very long-distance. Many intracellular proteins and organelles have to be transported down them after synthesis in neuron body. The retrograde transport from axonal ends to neuron body is also necessary to degrade proteins, clear useless aggregates, and provide a signaling mechanism from the synapse side to the cell body [9]. The kinesin and dynein families move along tracks of cytoskeletal filaments, termed microtubules [10]. Any disturbance of axonal transport may have severe consequences for neuronal function causing generally known as axonal traffic jams that are the consequence of an aberrant axonal transport, and usually lead to local swellings and accumulations of needless organelles [10]. However, it is not clear enough what mechanism(s) and substances do cause the deficits in axonal transportation in long axons derived far away from the neuron cell body.

    The incidence of spheroids shifts from gracile nucleus to fasciculus with the disease progress, suggesting that the degenerating axonal terminals of DRG neurons back from the distal presynaptic parts in the gracile nucleus, toward the cell bodies in the DRGs. These results indicate that axonal degeneration is age dependent and follows a “dying-back” process of the primary sensory axons [4],[5].

    Axonal dystrophy, which shows the unique pathological feature including axonal swelling (spheroids) localized throughout the central nervous system (CNS) and peripheral nervous systems (PNS), occurs under many pathological and physiological conditions. It has been observed in the infantile neuroaxonal dystrophy in mutant mouse [11] and Hallervorden-Spatz disease [12], and in the vitamin E deficiency in rats [13],[14]. In addition, physiological ageing is known to accompany the formation of spheroids in dorsal column nuclei in humans [15],[16]. The swelling axons occur at preterminal or presynaptic portions of sensory axons in a progress of axonal dystrophy. However, why these axons in ascending spinal tract are vulnerable still remain enigmatic [17].

    During the course of our study on gad mice, it was speculated that a similar lesion seen in CNS may also occur in the opposite PNS derived from DRG neurons which have sensory axon ending in muscle spindle, and motor endplates on every muscle fiber. The axonal pathology was studied by using a method of combining acetylcholine-esterase and silver impregnation method [5],[18]. The sensory axons in gad mice show axonal degeneration and swelling at their distal ends in muscle spindles as early as around 20-day-old (Figure 1c). Similar degeneration occurs in the distal motor ends little later than that in sensory endings [18]. The progressive degeneration of motor nerves was recently confirmed in UCH-L1 knockout mouse model [19],[20]. There are several differences in spheroid formation in gad mice between CNS and PNS. Axonal dystrophy in ascending gracile tract are accompanied by marked swollen axons, whereas sensory and motor axons in PNS swell only slightly (Figure 1c) [5],[18]. The DRG neurons are morphologically intact until later stages as far as examined [5],[21]. Although motor neurons in the spinal ventral horn are morphologically normal even at advanced stages, the dying-back degeneration is seen at nerve terminals in the distal endplate zone. In pathological concepts, the distal spheroids with myelin loss develops first in the most distal parts of the axons, and, if abnormality persisted, the formation of spheroids is mostly limited at the tip of axons retained other parts intact and continue dying-back degeneration [22],[23]. The dying-back axonopathy seen in CNS of gad mice resembles Wallerian degeneration which shows spheroids and myelin loss [24]. It is important to understand mechanisms underlying axonal dystrophy in gad mice which may be caused by impaired axonal transport resulting inability to keep up with the metabolic demands of axon.

    The enhancement of substance P (SP)-like immunoreactivity was observed in axonal profiles of gracile nucleus and fasciculus in 13-month-old normal mice [25] and ageing rat [16], but no SP-positive axons were detected in gracile ascending tract of young 3- and 5-month-old mice. Although partial arrest of axonal transport of SP and/or enhanced SP-synthesis in DRG neurons are ascribable, the true reason is still unclear [26]. It was reported that some of spheroids in gracile nucleus contain cholecystokinin (CCK) which is one of neuropeptides thought to be a somatosensory neurotransmitter or neuromodulator [27]. Authors suggested that CCK may be partly involved in the formation of dystrophic axons in this region. Since axonal dystrophy of gracile sensory neurons in gad mice are likely to those in aged mice mentioned earlier [4],[20], it was expected that SP reactivity is enhanced in lesioned sites of gad mice.

    The SP reaction in gad mice is more intense throughout gracile tract, and is already visible by earlier stages than in normal aged mice (Figure 2a and 2b). As mentioned in chronic pain, reactive astrocytes are involved in the regulation of chronic pain [28]. Interestingly, the astrocytes in lesioned sites of gad mice are stained with both GFAP and SP, and are characterized by astrocytic hypertrophy and proliferation [29],[30]. Dacil N, et al. (2015) reported recently that, although SP detected within cytoplasm of astrocytes is taken in through SP-receptors, the preprotachykinin mRNA (precursor of SP) and SP are shown to be localized in reactive astrocytes around senile plaques, moreover, SP expression increases angiogenesis of brain vessels associated with β-amyloid plaques [31]. Application of SP prevents β-amyloid-induced impairment of cognitive processes through inhibition of β-amyloid-induced overexpression of potassium channel subunits [32]. It seems likely that astrocytes are stimulated to proliferate and grow rapidly by unknown factor(s) released from axons and/or other cells degenerating in nearby tissues to prevent triggering a perivascular inflammation. Noguchi et al. (1995) reported that the peripheral nerve axotomy of primary afferent neurons results in de novo expression of preprotchykinin mRNA in rat DRG neurons, and, after 2 weeks, an increase of SP immunoreactivity is observed in axons and their terminals of large myelinated fibers in gracile nucleus [33].

    From these reports, we postulated that the loss of peripheral nerve endings of gad mice may lead to abnormal pain sensations and other sensory abnormalities as with in peripheral nerve axotomy. The SP has both neuroprotective and neurodegenerative properties which are related to proper development of CNS and changes in progressive neurodegenerative diseases [34]. In this context, we remark a report that SP induces α-secretase-mediated amyloid precursor protein (APP) cleavage and to protect cerebellar granule neurons from attenuating amyloid-β protein (AβP) production, without modifying the steady state level of its precursor APP [35],[36]. Kowall et al. (1991) added a fact that SP prevents β-amyloid induced neuronal loss when coadministered it into cerebral cortex [37]. Moreover, the cognitive impairment of rat by AβP brain infusion was prevented by SP administration [38].

    Figure 2.  SP-like immunoreactivity in cervical spinal cord of normal (a and b) and gad (c and d) mice aged 122-day old. The (b) and (d) enlarged parts of gracile fasciculus (arrow) in normal and gad mice, respectively. In normal spinal cord, SP activity is more prominent in dorsal horn than in other parts (a). SP positive cells are rarely found in normal gracile fasciculus (b). By contrast, the dense networks of SP positive cells are seen not only in gracile fasciculus but also along marginal zone of lateral and ventral funiculi in gad mouse (c). Fine and dense networks of SP positive cells in the gracile fasciculus are shown in (d). Bar = 100 µm (a and c); 50 µm (b and d). The source is from [26].

    From several pathological features in gad mice, we thought that abnormal expression of APP may be involved in axonal dystrophy along the gracile tract, following an accumulation of AβP which appears to be generated by aberrant proteolytic cleavage of APP, and deposits as an insoluble form mainly in the extracellular matrix surrounded by activated astrocytes [39],[40]. Yamaguchi et al. (1990) proposed a “diffuse plaque” as a pathological term to define the early and immature stage of the senile plaque formation, and demonstrated that the majority of senile plaques in AD brain are diffuse rather than the classic-core plaques [40]. A striking feature of APP-IR and AβP-IR in gad mouse is a diffuse and intracellular distribution of the immunoreacted products, suggesting that there is an intracellular deficit to accumulate these proteins within the cytoplasm (Figure 3a and 3b) [41],[42]. It was also reported that APP accumulates within the swollen axons in the human brain lesions [43] and in ageing rat brain [44]. The APP synthesized in DRG neurons moves along axon by the fast-anterograde transport to the distal ends [45] and distributes in almost all neuron cytoplasm and some glial cells in the nervous system [46]. The astrocytes are known to express APP and AβP following experimentally induced neuronal damage [47][49] and in cell cultures [50]. The distribution and expression of SP in transgenic mice overexpressing human APP indicated that SP-immunoreactive β-amyloid plaques dramatically increased after 12 months, thereafter SP was expressed in a number of reactive astrocytes [51]. Furthermore, the morphological abnormalities at synapses in the cerebral cortices of aged macaques preceded the accumulation of APP within neurites and reactive astrocytes and microglia [52]. Indeed, aberrant processing of APP into AβP by these cells is implicated as the initiating factor in the formation of senile plaques [53]. Another finding that astrocytic cells are immunoreacted with both APP and AβP, leads to an idea that similar mechanism may function along the gracile tract of the gad mouse.

    The Ubiquitin (Ub) is a 76 amino acid polypeptide involved in an ATP-dependent and non-lysosomal degradation of abnormal and short-lived normal proteins [54], and in the pathogenesis of cellular inclusions, many of which have a fibrillar structure and contain cytoskeletal proteins [55]. The Ub also plays an important role in the proteolysis of abnormal proteins in aberrant cells [56].

    The gracile axonal dystrophy in gad mice is mainly represented by swelling of axonal ends and “dying-back” degeneration of afferent neurons [57], and by an accumulation of APP and AβP along gracile ascending tract [41],[42]. These abnormalities are first observed at the most distal ends of afferent axons in gracile nucleus as early as 20-day-old [4]. Despite a modest filamentous accumulation in swollen axons, Ub-positive products are not found in the cytoplasm by 4 weeks of age. This raises the possibility that there may not be enough abnormal filaments for their selective degradation by cytoplasmic Ub and/or a lysosomal proteolysis may already be activated to function cooperatively with Ub-proteolysis. The axonal terminals are not always Ub-positive in the earliest stages, even if they reveal an initial sign of the structural alterations due to axonal degeneration. Thereafter, when axons are filled with an excess amount of aberrant neurofilaments, the Ub immunoreactivity increases and becomes to be detected in these axonal terminals (Figure 3c and 3d) [58],[59].

    The distribution of ubiquitinated deposits in gracile tract of gad mice is closely related to the severity of axonal degeneration progressing retrogradely from distal ends of primary sensory neurons, followed by axonal dystrophy in the nervous tracts of the secondary sensory neurons in Clarke's nucleus and gracile nucleus [30]. Similar Ub distribution extents further throughout the thalamus to the sensory areas in cerebral cortex [59]. In contrast, The Ub deposits in the descending motor neuron system are seen in the grey matter surrounding motor neurons in the ventral horn of spinal cord, followed by an increase in the order of lumbar, thoracic, and cervical spinal cord. The axonal degeneration in the motor neuron system in gad mice starts from the distal axonal ends toward the upper cerebral cortex, especially in cortex motor areas [59]. Together, the transneuronal degeneration beyond PNS to the lower CNS extends further to the upper CNS [58],[59].

    Figure 3.  (a) APP and (b) AβP immunoreactivity in gracile nucleus (gr) at cervical spinal cord of gad mice aged 32-weeks old. Astrocytes are positive in AβP immunoreactivity (arrowheads) [42]. Ubiquitin immunoreactivity in gracile (Gr) nucleus of normal (c) and gad (d) mice aged 18-weeks old. A large number of ubiquitinated dot-like structures (arrow) are seen throughout the gracile tract. High magnification of Ubiquitin dot-like structures is shown in an inset in (d) [59]. Bars = 25 µm in (c) and (d) for (a) and (b). Bar = 5 µm in an inset of (d).

    In eleven years after the first report on gad mice, Saigoh et al. (1999) revealed gad mutation which is caused by an in-frame deletion including exons 7 and 8 of Uchl1 gene, encoding ubiquitin carboxy-terminal hydrolase (UCH) isozyme L1 (UCH-L1) selectively expressed in nervous system and testis/ovary. UCH-L1 is a 24.8 kDa, consisting of 223 amino acids and encoded by 9 exons with a transcript of 1172 bps in length, and located on human chromosome 4(4p14) [60]. UCH-L1 (originally known as PGP 9.5 or PARK5) plays a key role in ubiquitin turnover through its C-terminal hydrolase activity. Several functions of UCH-L1 were proposed including ubiquitin hydrolyze activity, ubiquitin ligase activity and stabilization of the mono-ubiquitin. Mutations in human UCH-L1 gene are associated with PD and other neurodegenerative disorders [62]. The gad mutation with approximately 3.5 Kb deletion exists at genomic DNA level, encoding a truncated Uch-L1 gene lacking a segment of 42 amino acids containing a catalytic residue.

    The UCH-L1 is, as mentioned above, a member of deubiquitinating enzymes which functions as a Ub ligase and also a role in stabilizing Ub monomers in vivo [61]. As a deubiquitinating enzyme activity, UCH-L1 facilitates Ub recycling and regulates the cellular pool of available Ub, giving the capacity to modulate many ubiquitin-dependent cellular processes. Thus, the lack of UCH-L1 expression in mice results in gad phenotype. Down-regulation and extensive oxidative modification of this enzyme are observed in the brains of AD patients as well as PD patients [63]. It is unclear which human diseases are related to the pathogenesis of gad mice.

    An isoleucine 93 to methionine amino acid mutation (I93M) of UCH-L1 was identified as a cause of autosomal dominant PD [64]. To make sure whether this I93M mutation occurs due to a gain of toxic function in vivo, a transgenic mice (Tg) expressing human UCH-L1I93M were cleated [65]. The midbrains of Tg mice contain increased amounts of insoluble UCH-L1 as compared with those of non-Tg mice, perhaps due to a toxic gain of function. Authors also revealed one of typical pathological phenotypes of human PD, namely “dopaminergic neuron loss” in the substantia nigra, but did not detect α-synuclein in dystrophic spheroids [65],[66] and a sign of “Lewy body formation” [62]. These observations suggest that the molecular mechanism of β- and γ-synuclein accumulation in gad spheroids is different from that in Lewy bodies. Alzheimer's disease (AD) is characterized by amyloid plaques (Senile plaques) composed of β-amyloid (Aβ) peptide surrounded by swollen presynaptic dystrophic neurites consisting of dysfunctional axons and terminals that accumulate the amyloid precursor protein (APP) required for Aβ generation. The primary lesions observed in gracile nucleus of gad mice, are characterized by dystrophic swollen axons, and surrounded by the projections of astrocytes. As mentioned earlier, swollen axonal dystrophy with a dying back-type degeneration occurs from both central and peripheral distal ends of DRG neurons.

    Leaving secondary neurons in gracile nuclei behind, axonal dystrophy gradually retreats towards parent DRG neurons. The secondary neurons then continue to ascend the midline of medulla and cross to the other side and ascend to thalamus. Since Ubiquitin accumulation is often seen in thalamus neurons in advanced stages, the functional and metabolic deficits in secondary neurons occur already to some extent, resulting in partial synaptic loss of 3rd order neurons in thalamus. The 3rd neurons project their axons to somatosensory cerebral cortex, where somatosensory stimuli are processed.

    On the other hand, the motor nerve fibers descend through corticospinal tract to motor neurons in spinal grey matter. The UCH-L1 in maintaining corticospinal motor neurons keeps viability and cellular integrity. Such motor neurons with a long distance from cortex, show early, selective, progressive, and profound cell loss in the absence of UCH-L1 [67],[68]. From the results that maintenance of memory in a passive avoidance test and exploratory behavior in an open field test is reduced, and the long-term potentiation (LTP) of field synaptic responses from Schaffer collaterals to CA1 pyramidal cells in hippocampal slices is also impaired in gad mice [69],[70], a collapse of the synapse-axon circulation around CNS from PNS seems to take place. Thus, the true mechanism(s) and functions to cause axonal dystrophy in gad mice should be clarified more.



    [1] G. W. Dick, S. F. Kitchen, A. J. Haddow, Zika virus. I. isolations and serological specificity, Trans. R. Soc. Trop. Med. Hyg., 46 (1952), 509–520. https://doi.org/10.1016/0035-9203(52)90042-4 doi: 10.1016/0035-9203(52)90042-4
    [2] G. W. Dick, S. F. Kitchen, A. J. Haddow, Zika virus. II. pathogenicity and physical properties, Trans. R. Soc. Trop. Med. Hyg., 46 (1952), 521–534. https://doi.org/10.1016/0035-9203(52)90043-6 doi: 10.1016/0035-9203(52)90043-6
    [3] M. R. Duffy, T. H. Chen, W. T. Hancock, A. M. Powers, J. L. Kool, R. S. Lanciotti, Zika virus outbreak on Yap Island, federated states of Micronesia, N. Engl. J. Med., 360 (2009), 2536–2543. https://doi.org/10.1056/NEJMoa0805715 doi: 10.1056/NEJMoa0805715
    [4] S. Ioos, H. P. Mallet, I. L. Goffart, V. Gauthier, T. Cardoso, M. Herida, Current Zika virus epidemiology and recent epidemics, Med. Mal. Infect., 44 (2014), 302–307. https://doi.org/10.1016/j.medmal.2014.04.008 doi: 10.1016/j.medmal.2014.04.008
    [5] V. M. C. Lormeau, C. Roche, A. Teissier, E. Robin, A. Berry, H. Mallet, et al., Zika virus, French Polynesia, South Pacific, 2013, Emerg. Infect. Dis., 20 (2014), 1085–1086. https://doi.org/10.3201/eid2006.140138
    [6] C. Zanluca, V. C. A. Melo, A. L. P. Mosimann, G. I. V. Santos, C. N. D. Santos, K. Luz, First report of autochthonous transmission of Zika virus in Brazil, Inst. Oswaldo Cruz., 110 (2015), 569–572. https://doi.org/10.1590/0074-02760150192 doi: 10.1590/0074-02760150192
    [7] D. Gatherer, A. Kohl, Zika virus: a previously slow pandemic spreads rapidly through the Americas, J. Gen. Virol., 97 (2016), 269–273. https://doi.org/10.1099/jgv.0.000381 doi: 10.1099/jgv.0.000381
    [8] J. P. Messina, M. Kraemer, O. J. Brady, D. M. Pigott, F. M. Shearer, D. J. Weiss, et al., Mapping global environmental suitability for Zika virus, elife, 5 (2016), e15272. https://doi.org/10.7554/eLife.15272 doi: 10.7554/eLife.15272
    [9] V. L. P. Junior, K. Luz, R. Parreira, P. Ferrinho, Zika virus: a review to clinicians, Acta Med. Port., 28 (2015), 760–765.
    [10] R. Becker, Missing link: Animal models to study whether Zika causes birth defects, Nat. Med., 22 (2016), 225–227. https://doi.org/10.1038/nm0316-225 doi: 10.1038/nm0316-225
    [11] L. Bouzid, O. Belhamiti, Effect of seasonal changes on predictive model of bovine babesiosis transmission, Int. J. Model. Simul. Sci. Comput., 8 (2017), 1–17. https://doi.org/10.1142/S1793962317500301 doi: 10.1142/S1793962317500301
    [12] E. B. Hayes, Zika virus outside Africa, Emerg. Infect. Dis. 15 (2009), 1347–1350. https://doi.org/10.3201/eid1509.090442
    [13] M. H. Maamar, L. Bouzid, O. Belhamiti, F. B. M. Belgacem, Stability and numerical study of theoretical model of Zika virus transmission, Int. J. Math. Modell. Numer. Optim., 10 (2020), 141–166. https://doi.org/10.1504/IJMMNO.2020.106528 doi: 10.1504/IJMMNO.2020.106528
    [14] F. L. H. Wertheim, P. Horby, J. P. Woodall, Atlas of Human Infectious Diseases, Wiley-Blackwell, Oxford, 2012.
    [15] W. O. Kermack, A. G. McKendrick, A contribution to the mathematical theory of epidemics, Proc. R. Soc. London, Ser. A, 115 (1927), 700–721. https://doi.org/10.1098/rspa.1927.0118 doi: 10.1098/rspa.1927.0118
    [16] C. Manore, M. Hyman, Mathematical models for fighting Zika virus, SIAM News, 49 (2016).
    [17] D. Gao, Y. Lou, D. He, T. C. Porco, Y. Kuang, G. Chowell, et al., Prevention and control of Zika as a mosquito-borne and sexually transmitted disease: a mathematical modeling analysis, Sci. Rep., 6 (2016), 28070. https://doi.org/10.1038/srep28070 doi: 10.1038/srep28070
    [18] E. K. Lee, Y. Liu, F. H. Pietz, A compartmental model for Zika virus with dynamic human and vector populations, AMIA Annu. Symp. Proc., 2016 (2011), 743–752.
    [19] H. Nishiura, R. Kinoshita, K. Mizumoto, Y. Yasuda, K. Nah, Transmission potential of Zika virus infection in the south pacific, Int. J. Infect. Dis., 45 (2016), 95–97. https://doi.org/10.1016/j.ijid.2016.02.017 doi: 10.1016/j.ijid.2016.02.017
    [20] S. Rezapour, H. Mohammadi, M. E. Samei, SEIR epidemic model for COVID-19 transmission by Caputo derivative of fractional order, Adv. Differ. Equations, 2020 (2020), 1–19. https://doi.org/10.1186/s13662-020-02952-y doi: 10.1186/s13662-020-02952-y
    [21] S. Ahmad, M. Rahman, M. Arfan, On the analysis of semi-analytical solutions of Hepatitis B epidemic model under the Caputo-Fabrizio operator, Chaos, Solitons Fractals, 146 (2021), 110892. https://doi.org/10.1016/j.chaos.2021.110892 doi: 10.1016/j.chaos.2021.110892
    [22] M. A. Taneco-Hernández, C. Vargas-De-León, Stability and Lyapunov functions for systems with Atangana-Baleanu Caputo derivative: an HIV/AIDS epidemic model, Chaos, Solitons Fractals, 132 (2020), 109586. https://doi.org/10.1016/j.chaos.2019.109586 doi: 10.1016/j.chaos.2019.109586
    [23] S. Qureshi, R. Jan, Modeling of measles epidemic with optimized fractional order under Caputo differential operator, Chaos, Solitons Fractals, 145 (2021), 110766. https://doi.org/10.1016/j.chaos.2021.110766 doi: 10.1016/j.chaos.2021.110766
    [24] W. Wang, M. Zhou, T. Zhang, Z. Feng, Dynamics of a Zika virus transmission model with seasonality and periodic delays, Commun. Nonl. Sci. Numer. Simul., 116 (2023), 106830. https://doi.org/10.1016/j.cnsns.2022.106830 doi: 10.1016/j.cnsns.2022.106830
    [25] M. A. Ibrahim, A. Dénes, A mathematical model for Zika Virus infection and microcephaly risk considering sexual and vertical transmission, Axioms, 12 (2023), 263. https://doi.org/10.3390/axioms12030263 doi: 10.3390/axioms12030263
    [26] M. Murugappan, R. Grienggarai, V. Govindan, Mathematical modelling on the transmission dynamics of Zika Virus, Int. J. Robot., Autom. Sci., 5 (2023), 79–84. https://doi.org/10.33093/ijoras.2023.5.2.9 doi: 10.33093/ijoras.2023.5.2.9
    [27] F. A. Oguntolu, O. J. Peter, A. Yusuf, B. I. Omede, G. Bolarin, T. A. Ayoola, Mathematical model and analysis of the soil-transmitted helminth infections with optimal control, Model. Earth Syst. Environm., 2023 (2023), 1–15. https://doi.org/10.1007/s40808-023-01815-1 doi: 10.1007/s40808-023-01815-1
    [28] O. J. Peter, H. S. Panigoro, A. Abidemi, M. M. Ojo, F.A. Oguntolu, Mathematical model of COVID-19 pandemic with double dose vaccination, Acta Biotheor., 71 (2023), 9. https://doi.org/10.1007/s10441-023-09460-y doi: 10.1007/s10441-023-09460-y
    [29] A. J. Kucharski, S. Funk, R. M. Eggo, H. P. Mallet, W. J. Edmunds, E. J. Nilles, Transmission dynamics of Zika virus in island populations: a modelling analysis of the 2013/14 French Polynesia outbreak, PLoS Negl. Trop. Dis., 10 (2016), e0004726. https://doi.org/10.1371/journal.pntd.0004726 doi: 10.1371/journal.pntd.0004726
    [30] J. P. T. Boorman, J. S. Porterfield, A simple technique for infection of mosquitoes with viruses, transmission of Zika virus, Trans. Roy. Soc. Trop. Med. Hyg., 50 (1956), 238–242. https://doi.org/10.1016/0035-9203(56)90029-3 doi: 10.1016/0035-9203(56)90029-3
    [31] M. Darwish, A seroepidemiological survey for bunyaviridae and certain other arboviruses in Pakistan, Trans. R. Soc. Trop. Med. Hyg., 77 (1983), 446–450. https://doi.org/10.1016/0035-9203(83)90108-6 doi: 10.1016/0035-9203(83)90108-6
    [32] C. Castillo-Chávez, H. Thieme, Asymptotically autonomous epidemic models, in Mathematical Population Dynamics: Analysis of Heterogeneity, Wuerz Publishing, Winnipeg, (1995), 33–50.
    [33] G. M. R. Costa, M. Lobosco, M. Ehrhardt, R. F. Reis, Mathematical analysis and a nonstandard scheme for a model of the immune response against COVID-19, in Mathematical and Computational Modeling of Phenomena Arising in Population Biology and Nonlinear Oscillations: In honour of the 80th birthday of Ronald E. Mickens, AMS Contemporary Mathematics, 2023.
    [34] K. Diethelm, A fractional calculus based model for the simulation of an outbreak of Dengue fever, Nonl. Dyn., 71 (2013), 613–619. https://doi.org/10.1007/s11071-012-0475-2 doi: 10.1007/s11071-012-0475-2
    [35] H. K. Khalil, Nonlinear Systems, Prentice-Hall, London, UK, 1996.
    [36] R. Gorenflo, A. A. Kilbas, F. Mainardi, S. V. Rogosin, Mittag-Leffler Functions, Related Topics and Applications, Springer, Berlin-Heidelberg, 2014.
    [37] C. Vargas-De-León, Volterra-type Lyapunov functions for fractional-order epidemic systems, Commun. Nonl. Sci. Numer. Simul., 24 (2015), 75–85. https://doi.org/10.1016/j.cnsns.2014.12.013 doi: 10.1016/j.cnsns.2014.12.013
    [38] R. E. Mickens, Applications of Nonstandard Finite Difference Schemes, World Scientific, 2000.
    [39] M. M. Khalsaraei, Positivity of an explicit Runge–Kutta method, Ain Shams Eng. J., 6 (2015), 1217–1223. https://doi.org/10.1016/j.asej.2015.05.018 doi: 10.1016/j.asej.2015.05.018
    [40] A. Gerisch, R. Weiner, The positivity of low-order explicit Runge-Kutta schemes applied in splitting methods, Comput. Math. Appl., 45 (2003), 53–67. https://doi.org/10.1016/S0898-1221(03)80007-X doi: 10.1016/S0898-1221(03)80007-X
    [41] R. E. Mickens, Calculation of denominator functions for nonstandard finite difference schemes for differential equations satisfying a positivity condition, Numer. Methods Partial Differ. Equations, 23 (2007), 672–691. https://doi.org/10.1002/num.20198 doi: 10.1002/num.20198
    [42] S. Berkhahn, M. Ehrhardt, A physics-informed neural network to model COVID-19 infection and hospitalization scenarios, Adv. Contin. Discrete Models, 2022 (2022), 61. https://doi.org/10.1186/s13662-022-03733-5 doi: 10.1186/s13662-022-03733-5
    [43] S. Treibert, H. Brunner, M. Ehrhardt, A nonstandard finite difference scheme for the SVICDR model to predict COVID-19 dynamics, Math. Biosci. Eng., 19 (2022), 1213–1238. https://doi.org/10.3934/mbe.2022056 doi: 10.3934/mbe.2022056
    [44] M. Ehrhardt, R. E. Mickens, A Nonstandard Finite Difference Scheme for Solving a Zika Virus Model, unpublished manuscript, 2017.
    [45] D. Matignon, Stability results for fractional differential equations with applications to control processing, Comput. Eng. Syst. Appl., 2 (1996), 963–968.
    [46] A. Ali, S. Islam, M. R. Khan, S. Rasheed, F. M. Allehiany, J. Baili, et al., Dynamics of a fractional order Zika virus model with mutant, Alexandria Eng. J., 61 (2022), 4821–4836. https://doi.org/10.1016/j.aej.2021.10.031 doi: 10.1016/j.aej.2021.10.031
    [47] A. Ali, F. S. Alshammari, S. Islam, M. A. Khan, S. Ullah, Modeling and analysis of the dynamics of novel coronavirus (COVID-19) with Caputo fractional derivative, Res. Phys., 20 (2021), 103669. https://doi.org/10.1016/j.rinp.2020.103669 doi: 10.1016/j.rinp.2020.103669
  • This article has been cited by:

    1. Suyan Tan, Yilin Guo, A study of the impact of scientific collaboration on the application of Large Language Model, 2024, 9, 2473-6988, 19737, 10.3934/math.2024963
    2. Tekin Tuncer, Küresel Ticaret Yolları ve Çin Doğu Türkistan'a Etkileri, 2024, 9, 2149-9535, 2370, 10.24186/vakanuvis.1567429
    3. Mumtaz Ahmed, Liviu Marian Matac, Naureen Maqbool, Asma Salman, Muthanna G. Abdul Razzaq, Lara Al-Haddad, Codruta Daniela Pavel, Bubbles Identification in an Emerging Economy and Within Stock Markets of its Trading Partners: Evidence from a GSADF Approach, 2025, 2010-4952, 10.1142/S2010495224500179
  • 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(1887) PDF downloads(85) Cited by(10)

Figures and Tables

Figures(20)  /  Tables(4)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog