Next Article in Journal
Epidemic Waves in a Stochastic SIRVI Epidemic Model Incorporating the Ornstein–Uhlenbeck Process
Previous Article in Journal
Bifurcation Analysis Software and Chaotic Dynamics for Some Problems in Fluid Dynamics Laminar–Turbulent Transition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Universal Theory for Multiscale Modelling of Infectious Disease Dynamics

Modelling Health and Environmental Linkages Research Group (MHELRG), Department of Mathematical and Computational Sciences, University of Venda, Private Bag X5050, Thohoyandou 0950, South Africa
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(18), 3874; https://doi.org/10.3390/math11183874
Submission received: 27 July 2023 / Revised: 4 September 2023 / Accepted: 5 September 2023 / Published: 11 September 2023
(This article belongs to the Special Issue Multiscale Modeling of Infectious Disease Dynamics)

Abstract

:
The replication-transmission relativity theory, currently used to inform the development of multiscale models of infectious disease dynamics, needs a revision and extension to accommodate new basic science and clinical information about infectious disease dynamics. In this article, we revise and extend the replication-transmission relativity theory into a new scientific theory of infectious disease dynamics called the universal theory for the multiscale modelling of infectious disease dynamics. This new theory states that, for every host–pathogen interaction that results in an infectious disease system, there is no privileged or absolute scale of a disease system form that would determine the dynamics of the infectious disease system, only interactions between the scales of a level of organisation of the pathogen-centred disease system form and the scales of the corresponding levels of organisation of the host-centred disease system form. We further explain the utility of this theory, which is reflected in its flexibility and ability to incorporate new information and explain previous information that could not be accounted for by the replication-transmission relativity theory of infectious disease dynamics.

1. Introduction

Progress in the study of infectious disease systems dynamics as complex systems using multiscale modelling methods needs to be made over with a new scientific theory that represents a paradigm shift by revising and extending the replication-transmission relativity theory [1]. However, new scientific theories succeed when they re-order old knowledge into a new framework such that the new framework which they introduce makes it possible to break barriers in our understanding of many open questions and solve problems deemed important but barely addressed. The replication-transmission relativity theory [1], which states that at every level of organisation of an infectious disease system there is no privileged or absolute scale which would determine disease dynamics, only interactions between the microscale and macroscale, is presently used to inform the study of infectious disease dynamics using multiscale modelling methods. This theory is not well suited to capturing the dynamic complexity of infectious disease systems because it only considers their multiscale dynamics from a pathogen-centred perspective [2] while disregarding their multiscale dynamics from a host-centred perspective.
The aim of this article is to revise and extend the replication-transmission relativity theory [1] in a radical way into a new scientific theory of infectious disease dynamics called the universal theory for multiscale modelling of infectious dynamics, which incorporates both the pathogen-centred and the host-centred perspectives to the multiscale dynamics of infectious disease systems. This new theory states that, for every host-pathogen interaction that results in an infectious disease system, there is no privileged or absolute scale of a disease system form which would determine the dynamics of the infectious disease system, only interactions between the scales of a level of organisation of the pathogen-centred disease system form and the scales of the corresponding levels of organisation of the host-centred disease system form. This theory puts forth the view that the multiscale modelling of host-pathogen interaction in infectious disease dynamics is about considering the interaction of the scales of a level of organisation of the pathogen-centred disease system form on one hand and the scales of the corresponding levels of organisation of the host-centred disease system form on the other hand. Each of these two different disease system forms has a different multilevel and multiscale organisation and thus requires a different multilevel and multiscale modelling approach; that is, there are two different multiscale modelling approaches of the two different disease system forms of an infectious disease system, each of which is only partially true or provides partial information about the multiscale dynamics of an infectious disease system. Therefore, the pathogen-centred disease system form and the host-centred disease system form are considered to be mutually exclusive and complementary disease system forms of a single infectious disease system.
The weakness of the replication-transmission relativity theory arises from the fact that it privileges the scales of a level of organisation of the pathogen-centred disease system form while neglecting the scales of the corresponding levels of organisation of the host-centred disease system form in the multiscale dynamics of an infectious disease system. However, the universal theory for multiscale modelling of infectious disease dynamics requires that both the scales of a level of organisation of the pathogen-centred disease system form and the scales of the corresponding levels of organisation of the host-centred disease system form be integrated into a single multiscale model of infectious disease dynamics. The contents of the rest of this paper are organised as follows. In Section 2, we interpret the universal theory for the multiscale modelling of infectious disease dynamics. This is followed by Section 3, where we present the multilevel and multiscale organisation of the pathogen-centred disease system form. The multilevel and multiscale organisation of the host-centred disease system form is presented in Section 4. To illustrate ideas, we present a process-based multiscale model of the pathogen-centred disease system form of a malaria disease system as an example in Section 5. The dynamics of the pathogen-centred disease system form of the malaria disease system is analysed using the reproductive number in Section 6 and further analysed in Section 7 using numerical methods. The paper ends with concluding remarks in Section 8.

2. Interpretation of the Universal Theory for Multiscale Modelling of Disease Dynamics

The universal theory for the multiscale modelling of infectious disease dynamics, which states that, for every host–pathogen interaction that results in an infectious disease system, there is no privileged or absolute scale of a disease system form that would determine the dynamics of the infectious disease system, only interactions between the scales of a level of organisation of the pathogen-centred disease system form and the scales of the corresponding levels of organisation of the host-centred disease system form, is a radical theory for the multiscale modelling of disease dynamics based on four key assumptions, which are as follows:
[I.]
First, that infectious disease dynamics are a result of the interaction of three subsystems, which are the host subsystem, which we call the primary subsystem; the environment subsystem, which we call the secondary subsystem; and the pathogen subsystem, which we call the tertiary subsystem.
[II.]
Second, that the environment subsystem is considered to be an extended form of the host subsystem and thus an infectious disease system is considered to be a result of host–pathogen interaction only in such a way that every infectious disease system is considered to exist in two different forms: a host–centred disease system form, which is the form that the whole disease system takes when the host perspective of the interaction is considered, and a pathogen–centred disease system form, which is the form that the whole disease system takes when the pathogen perspective of the interaction is considered.
[III.]
Third, that the two different disease system forms; that is, the host-centred disease system form and the pathogen-centred disease system form, have different multilevel and multiscale organisations and are necessarily described by different types of multiscale models, which are process-based multiscale models for the pathogen-centred disease system form and mechanism-based multiscale models for the host-centred disease system form.
[IV.]
Fourth, that, together, the mechanism-based multiscale model of the host-centred disease system form and the process-based multiscale model of the pathogen-centred disease system form present a fuller description of the overall multiscale dynamics of an infectious disease system than either of the two when considered separately.
These four assumptions form a scaffold on which the universal theory, as a formal theory for the multiscale modelling of infectious disease dynamics, is built and tested on. Figure 1 shows a conceptual representation of the universal theory for the multiscale modelling of infectious disease dynamics.
From Figure 1, we note that, at the highest level, an infectious disease system is a result of the interaction of three subsystems: [a.] the host subsystem, which we call the primary subsystem, [b.] the environment subsystem, which we call the secondary subsystem, and [c.] the pathogen subsystem, which we call the tertiary subsystem. By invoking assumption II of the universal theory for the multiscale modelling of infectious disease dynamics, which considers the environment subsystem as an extended form of the host subsystem, we note that only two different perspectives are needed to describe the interaction of these three subsystems. First, there is the pathogen-centred perspective of the interaction of these three subsystems—which gives rise to the pathogen-centred disease system form. This disease system form is an outcome of the impact of the host subsystem on the pathogen subsystem in host–pathogen interaction. Second, there is the host-centred perspective of the interaction of these three subsystems, which gives rise to the host-centred disease system form. This disease system form is an outcome of the impact of the pathogen subsystem on the host subsystem in host–pathogen interaction. As illustrated in Figure 1, the pathogen-centred disease system form is organised into seven main levels of organisation, which are [1,2,3]: [a.] the cell level, [b.] the tissue level, [c.] the organ level or microcommunity level, [d.] the microecosystem level, [e.] the whole organism level, [f.] the macrocommunity level, and [g.] the macroecosystem level. However, the host-centred disease system form is organised into four main levels of organisation, which are: [a.] the host damage level, [b.] the host immune response level, [c.] the environmental change level, and [d.] the biodiversity change level. In addition, Figure 1 shows that each level of the pathogen-centred disease system form can be resolved into two hierarchically organised scales—a microscale and a macroscale—with both scales being hierarchically organised in both space and time. However, each level of the host-centred disease system form can be resolved into seven main scales—[a.] the molecular scale, [b.] the cell scale, [c.] the tissue scale, [d.] the organ scale, [e.] the whole organism scale, [f.] the community scale, and [g.] the ecosystem scale—that are hierarchically organised into physical size scale and time scale.
Further, we note from Figure 1 that the universal theory for the multiscale modelling of infectious disease dynamics is based on mathematical approaches to understanding the multiscale dynamics of infectious disease systems and expresses the mathematical point that the overall multiscale model of infectious disease dynamics M consists of two complementary multiscale models: M p , a process-based multiscale model that describes the scales of a level of organisation of the pathogen-centred disease system form, and M m , a mechanism-based multiscale model that describes the scales of the corresponding levels of organisation of the host-centred disease system form so that M = f ( M p , M m ) . As shown in Figure 1, the process-based multiscale model M p of the scales of a level of organisation of the pathogen-centred disease system form is described in terms of four main disease process, which are [a.] the infection or the super-infection by the pathogen process, which involves the movement of the pathogen from the macroscale to the microscale, and thus the term infection is defined as the acquisition of a pathogen at the microscale from the macroscale at any level of organisation of the pathogen-centred disease system form; [b.] the pathogen replication process at the microscale; [c.] the pathogen shedding or excretion process, which involves the movement of the pathogen from the microscale to the macroscale; and [d.] the pathogen transmission process at the macroscale. However, the mechanism-based multiscale model M m of the scales of the corresponding levels of organisation of the host-centred disease system form is described in terms of up to a maximum of seven main scale mechanisms, which are: [a.] the molecular-scale mechanism, [b.] the cellular-scale mechanism, [c.] the tissue-scale mechanism, [d.] the organ-scale mechanism, [e.] the organism-scale mechanism, [f.] the community-scale mechanism, and [g.] the ecosystem-scale mechanism.
Therefore, in the context of the universal theory for the multiscale modelling of infectious disease dynamics, conceptually represented by Figure 1, the two different disease system forms arise due to the need to account for the two reciprocal impacts of the pathogen subsystem and the host subsystem on each other in the multiscale dynamics of an infectious disease system as follows: [a.] the impact of the host subsystem on the pathogen subsystem, which gives rise to the pathogen-centred disease system form, and [b.] the impact of the pathogen subsystem on the host subsystem, which gives rise to the host-centred disease system form.
We then conclude from Figure 1 that the universal theory for the multiscale modelling of infectious disease dynamics can be interpreted to mean that, in host–pathogen interaction, the resulting infectious disease system from the mutual impacts of the host subsystem and pathogen subsystem on each other is a pair of interacting disease system forms, which are the pathogen-centred disease system form and the host-centred disease system form, with different multiscale and multilevel organisations and reciprocal influence on each other. This duality in interacting disease system forms with different multiscale and multilevel organisations for each infectious disease system can be interpreted to imply a law of multiscale dynamics of infectious disease systems that states that, for every set of disease processes, which are [a.] the infection/super-infection by pathogen process, [b.] the pathogen replication process, [c.] the pathogen shedding/excretion process, and [d.] the pathogen transmission process, at a level of organisation of the pathogen-centred disease system form, there is an associated set of disease mechanisms, which are [a.] the molecular-scale mechanism, [b.] the cell-scale mechanism, [c.] the tissue-scale mechanism, [d.] the organ-scale mechanism, [e.] the whole-organism-scale mechanism, [f.] the community-scale mechanism, and [g.] the ecosystem-scale mechanism, at corresponding levels of organisation of the host-centred disease system form with reciprocal influence between the processes and mechanisms. This law can be considered as an extension of Newton’s third law of motion in mechanics applied to infectious disease dynamics by requiring the disease processes at a level of organisation of the pathogen-centred disease system form and the disease mechanisms at the corresponding levels of organisation of the host-centred disease system form to be considered as characterising action and reaction force pairs in host–pathogen interaction. In the following two sections, we describe in detail the different multilevel and multiscale organisations and possible multiscale modelling frameworks of the pathogen-centred disease system form and the host-centred disease system form.

3. The Pathogen-Centred Disease System Form

In the context of the universal theory for the multiscale modelling of infectious disease dynamics, which is conceptually represented by Figure 1, the pathogen-centred disease system form arises due to the need to account for the impact of the host subsystem on the pathogen subsystem. In what follows, we discuss in detail the multilevel organisation of the pathogen-centred disease system form in Section 3.1 and the multiscale organisation of a level of organisation of the pathogen-centred disease system form in Section 3.2 together with the different categories of multiscale models that are used to describe the pathogen-centred disease system form at different levels of its organisation.

3.1. The Multilevel Organisation of the Pathogen-Centred Disease System Form

The pathogen-centred disease system form is the form that the whole infectious disease system takes when considering the pathogen-centred perspective of the interaction among the three subsystems of an infectious disease system; that is, the host subsystem, the environment subsystem, and the pathogen subsystem. In order to analyse the pathogen-centred disease system form using multiscale modelling methods, its complexity has to be brought down to manageable levels by discretising or decomposing it into different discrete levels of organisation. Based on the structural organisation of the host subsystem and the associated environment subsystem as an extended form of the host subsystem as a habitat for infectious agents, we establish that a pathogen-centred disease system form can be discretised into seven main hierarchical levels of organisation, which are [1,2,3]: [a.] the cell level, [b.] the tissue level, [c.] the organ or microcommunity level, [d.] the microecosystem level, [e.] the whole organism level, [f.] the macrocommunity level, and [g.] the macroecosystem level. Figure 2 shows a schematic representation of the seven main hierarchical levels of organisation of a pathogen-centred disease system form in a space-time diagram that portrays the hierarchical nature of the levels, showing a positive correlation in the spatial and temporal scales of varying disease processes.
In what follows, we briefly describe the seven main levels of organisation of the pathogen-centred disease system form and the associated scales; that is, the microscale and the macroscale, for each level of organisation.
[a.]
The cell level: The two hierarchically organised scales at this level of organisation of the pathogen-centred disease system form are the within-cell scale, which is the microscale, and the between-cell scale, which is the macroscale.
[b.]
The tissue level: The two hierarchically organised scales at this level of organisation of the pathogen-centred disease system form are the within-tissue scale, which is the microscale, and the between-tissue scale, which is the macroscale.
[c.]
The organ level or microcommunity level: For this level of organisation of the pathogen-centred disease system form, the two hierarchically organised scales are the within-organ or within-microcommunity scale, which is the microscale, and the between-organ or between-microcommunity scale, which is the macroscale.
[d.]
The microecosystem level: The two hierarchically organised scales at this level of organisation of the pathogen-centred disease system form are the within-microecosystem scale, which is the microscale, and the between-microecosystem scale, which is the macroscale.
[e.]
The whole organism level: At this level of organisation of the pathogen-centred disease system form, the two hierarchically organised scales; that is, the microscale and the macroscale, are the within-whole organism scale and the between-whole organism scale, respectively.
[f.]
The macrocommunity level: This level of organisation of the pathogen-centred disease system form consists of the within-macrocommunity scale as the microscale and the between-macrocommunity scale as the macroscale.
[g.]
The macroecosystem level: The two hierarchically organised scales at this level of organisation of the pathogen-centred disease system form are the within-macroecosystem scale, which is the microscale, and the between-macroecosystem scale, which is the macroscale.
As shown in Figure 2, the seven different levels of organisation of the pathogen-centred disease system form are hierarchically organised in both space and time. These seven levels of organisation can be demarcated into three main groups as follows. First, a group of primary levels of organisation consisting of the cell level, the tissue level, and the whole organism level, where the transmission process of the pathogen-centred disease system form is only through the local exchange of pathogen between the microscale and the macroscale in the context of single host species and single pathogen species/strain. Second, a group of secondary levels of organisation consisting of the organ/microcommunity level and the macrocommunity level, where the transmission process of the pathogen-centred disease system form is through both the local and global exchange of pathogen between the microscale and the macroscale, involving single host species and single pathogen species/strain in the context of multiple heterogeneous environments. Third, a group of tertiary levels of organisation consisting of the microecosystem level and the macroecosystem level, where the transmission process of the pathogen-centred disease system form is through a local or both local and global exchange of pathogen between the microscale and the macroscale in the context of multiple host species and or multiple strains/species of pathogens.

3.2. The Multiscale Organisation of the Pathogen-Centred Disease System Form

For the pathogen-centred disease system form, each of its seven levels of organisation illustrated in Figure 2 consists of two scales—a microscale and a macroscale. Figure 3 is a conceptual representation of a typical primary level of organisation of the pathogen-centred disease system form and the interaction of the microscale and macroscale through the four key disease processes [2]: [a.] the infection or super-infection by pathogen process, [b.] the pathogen replication process, [c.] the pathogen shedding or excretion process, and [d.] the pathogen transmission process. The interaction between the microscale and the macroscale of the pathogen-centred disease system form through these four disease processes is described by process-based multiscale models. These multiscale models belong to five different categories, which are [1,2,3]: [a.] individual-based multiscale models (IMSMs), which we also call phenomenological-based multiscale models, [b.] nested multiscale models (NMSMs), [c.] embedded multiscale models (EMSMs), [d.] hybrid multi-scale models (HMSMs), and [e.] coupled multiscale models (CMSMs). This taxonomic categorisation of multiscale models of infectious disease dynamics has been previously discussed in [1,2,3]. However, it is being further discussed in this article as part of the revision of the replication transmission relativity theory [1] to reflect the new knowledge based on the universal theory of disease dynamics, which dictates that these different categories of multiscale models are process-based multiscale models and thus describe the pathogen-centred disease system form only. In addition, and as illustrated in Figure 3, the development of these different categories of process-based multiscale models involves linking or integrating sub-models across scales by up-scaling and down-scaling variables associated with disease processes across the microscale and the macroscale. This is because, for the pathogen-centred disease system form, the hierarchical scales at any of its seven levels of organisation; that is, the microscale and the macroscale, indicate or represent shifts in the four disease processes; that is, the infection or super-infection by pathogen process, the pathogen replication process, the pathogen shedding or excretion process, and the pathogen transmission process.
In light of the new knowledge based on the universal theory of disease dynamics, we note that it is only the interaction of the microscale and the macroscale at a particular level of organisation of the pathogen-centred disease system form that necessarily assumes the replication-transmission relativity theory [1] and not the complete infectious disease system. Therefore, in the context of the universal theory for multiscale modelling of disease dynamics, process-based multiscale models refer to multiscale models of the pathogen-centred disease system form developed by explicitly incorporating four main disease processes that occur in the pathogen subsystem due to the impact of the host subsystem, which include [a.] the infection or super-infection by the pathogen process, [b.] the pathogen replication process, [c.] the pathogen excretion or shedding process, and [d.] the pathogen transmission process, to explain the multiscale temporal and spatial changes of the pathogen-centred disease system form at any of its seven main levels of organisation, which are [a.] the cell level, [b.] the tissue level, [c.] the organ/microcommunity level, [d.] the microecosystem level, [e.] the whole organism level, [f.] the macrocommunity level, and [g.] the macroecosystem level, without taking into consideration how the host subsystem causes these multiple scale processes. In what follows, we briefly describe the five different categories of the process-based multiscale models of the pathogen-centred disease system form and also give examples.
[a.]
Category I—Individual-Based Multiscale Models (IMSMs): A generic category of multiscale models that provides a simplified means of describing the multiscale dynamics of the pathogen-centred disease system form at a particular level of its organisation. We also sometimes refer to these individual-based multiscale models (IMSMs) as phenomenological-based multiscale models. This is because in this category of multiscale models, only the microscale is explicitly incorporated into the multiscale model and, as a result, the four main disease processes, which include the infection or super-infection by pathogen process, the pathogen replication process, the pathogen shedding or excretion process, and the pathogen transmission process, are only phenomenologically incorporated into the multiscale model with no explicit representation. In these phenomenological-based multiscale models, no macroscale sub-model is considered and so the macroscale is often observed as emergent behaviour of the microscale entities’ behaviour. In order to describe the macroscale, the microscale results are converted by summing, averaging, or performing some detailed statistical analysis of them and aggregating the information into macroscale variables for interpretation at that scale. By failing to explicitly incorporate details of disease processes, phenomenological-based multiscale models only make reference to the complexity of the multiscale dynamics of the pathogen-centred disease system form at a level of its organisation without incorporating the exact content of the complexity. Examples of phenomenological-based multiscale models are agent-based models [4] and graph-theoretic models [5]. Because they implicitly incorporate disease processes of the pathogen-centred disease system form, phenomenological-based multiscale models are typically simpler than the other four categories of multiscale models.
[b.]
Category II—Nested Multiscale Models (NMSMs): A generic category of process-based multiscale models developed to characterise the multiscale dynamics of the pathogen-centred disease system form at a particular level of its organisation by incorporating four main disease processes, which include the infection by pathogen process, the pathogen replication process, the pathogen shedding or excretion process, and the pathogen transmission process, across the microscale and the macroscale. In this category of process-based multiscale models, the microscale sub-model and the macroscale sub-model must be described by the same mathematical formalism or mathematical representation and are linked through the exchange of pathogen. For these multiscale models, the macroscale is linked to the microscale through the infection by pathogen process while the microscale is linked to the macroscale through the pathogen shedding or excretion process. A wide range of mathematical techniques can be used to integrate the microscale sub-model and the macroscale sub-model. This gives rise to different classes of process-based multiscale models of the pathogen-centred disease system form in this category of multiscale models. Typical examples of nested multiscale models include [6] for COVID-19 and [7] for paratuberculosis.
[c.]
Category III—Embedded Multiscale Models (EMSMs): A generic category of process-based multiscale models developed to characterise the multiscale dynamics of the pathogen-centred disease system form at a particular level of its organisation by incorporating four main disease processes, which include the super-infection by pathogen process, the pathogen replication process, the pathogen shedding or excretion process, and the pathogen transmission process, across the microscale and the macroscale. For this category of process-based multiscale models, the microscale sub-model and the macroscale sub-model must also be described by the same mathematical formalism or mathematical representation and are linked through the exchange of pathogen. Further, the macroscale is linked to the microscale through the super-infection by pathogen process while the microscale is linked to the macroscale through the pathogen shedding or excretion process. A wide range of mathematical techniques can also be used to integrate the microscale sub-model and the macroscale sub-model. This also gives rise to different classes of process-based multiscale models of the pathogen-centred disease system form in this category of multiscale models. Specific examples of embedded multiscale models are [7] for paratuberculosis and [1] for hookworm infection.
[d.]
Category IV: Hybrid multiscale models (HMSMs): A generic category of process-based multiscale models developed to characterise the multiscale dynamics of the pathogen-centred disease system form at a particular level of its organisation by incorporating four main disease processes, which include the infection or super-infection by pathogen process, the pathogen replication process, the pathogen shedding or excretion process, and the pathogen transmission process across the microscale and the macroscale. In this category of multiscale models, the microscale sub-model and macroscale sub-model are described by different mathematical formalisms or mathematical representations, depending on differences between scale variables of the microscale and the macroscale, which can be, for example, due to:
[i.]
The nature of variation in time scale variables, which can be discrete time at one scale and continuous time at the other scale;
[ii.]
The nature of variation in the state variables, which can be stochastic state variables at one scale and deterministic state variables at the other scale;
[iii.]
The nature of variation in the spatial scale, which can be homogeneous spatial variables at one scale described by ODEs, and heterogeneous spatial variables at the other scale described by PDEs, etc. In these process-based multiscale models, the macroscale is linked to the microscale through the infection by pathogen process or super-infection by pathogen process while the microscale is linked to the macroscale through the pathogen shedding or excretion process. A wide range of mathematical techniques can also be used to integrate the microscale sub-model and the macroscale sub-model. This gives rise to different classes of process-based multiscale models of the pathogen-centred disease system form in this category of multiscale models. Typical examples of hybrid multiscale models are [8,9,10].
[e.]
Category V—Coupled Multiscale Models (CMSMs): A generic category of multiscale models, which may be process-based multiscale models, phenomenological-based multiscale models, or a combination of these two types of multiscale models, that provides a simplified means of describing the multiscale dynamics of a pathogen-centred disease system form at a particular level of its organisation. For this category of multiscale models, the process-based multiscale models from any of categories I, II, III, and IV are used as sub-models in the development of multiscale models of the pathogen-centred disease system form. This is because the multiscale models developed in this category consider multiple pathogen strain infections, and/or multiple pathogen species infections, and/or multiple host group infections, and/or multiple host species infections, and/or multiple communities infections, and/or multiple organ/anatomical compartment infections. They are not like process-based multiscale models of categories I, II, III, and IV, which focus on a specific combination of (i) one-host and (ii) one-pathogen species/strain of a pathogen-centred disease system form relationships in multiscale modelling. Examples of coupled multiscale models are [11] for malaria and [12] for for river blindness, which are described in the context of multiple host species.
These different categories of multiscale models can be further demarcated into three groups. First, we have a group of primary process-based multiscale models, which are multiscale models developed at primary levels of organisation of the pathogen-centred disease system form, and are more suitable for an evaluation of the control of this disease system form. Second, we have a group of secondary process-based multiscale models, which are multiscale models developed at secondary levels of organisation of the pathogen-centred disease system form, and are more suitable for an evaluation of the elimination of this disease system form. Third and last, we have a group of tertiary process-based multiscale models, which are multiscale models developed at tertiary levels of organisation of the pathogen-centred disease system form, and are more suitable for an evaluation of the eradication of this disease system form. One of the pivotal concepts introduced in [11] is that process-based multiscale models can be used to translate existing knowledge about the efficacy of health interventions at the microscale into projected outcomes of the effectiveness of health interventions at the macroscale in the control, elimination, and even eradication of a pathogen-centred disease system form at any of its seven hierarchical levels of organisation. This has important applications in the evaluation of the effectiveness of drugs in different populations.

4. The Host-Centred Disease System Form

In the context of the universal theory for the multiscale modelling of infectious disease dynamics, which is conceptually represented by Figure 1, the host-centred disease system form arises due to the need to account for the impact of the pathogen subsystem on the host subsystem. We now discuss in detail the multilevel organisation of the host-centred disease system form in Section 4.1 and the multiscale organisation of a level of organisation of the host-centred disease system form in Section 4.2 together with the difficulties encountered in the development of mechanism-based multiscale models that are used to describe the host-centred disease system form.

4.1. The Multilevel Organisation of the Host-Centred Disease System Form

The host-centred disease system form is the form that the whole infectious disease system takes when considering the host-centred perspective of the interaction among the three subsystems of an infectious disease system; that is, the host subsystem, the environment subsystem, and the pathogen subsystem. In order to analyse a host-centred disease system form using multiscale modelling methods, its complexity also has to be brought down to manageable levels by discretising or decomposing it into different discrete levels of organisation so that, at each level of organisation, disease dynamics can be analysed in terms of their scales of organisation. Based on the functional organisation of living organisms and their associated environment in which pathogens generate effects, we establish that a host-centred disease system form can be discretised into four main hierarchical levels of organisation, which are: [a.] the host damage level, [b.] the host immune response level, [c.] the environmental change level, and [d.] the biodiversity change level. Figure 4 shows a conceptual representation of these four main levels of organisation of the host-centred disease system form.
Each of the four main levels of organisation of the host-centred disease system form consists of two main complementary level forms. In what follows, we describe the four main levels of organisation of the host-centred disease system form and the two main complementary level forms for each level of organisation.
[a.]
The host damage level: This functional level of organisation of the host-centred disease system form consists of two main complementary level forms: [i.] the pathogen-mediated damage level form and [ii.] the host-mediated damage level form. Usually, the host damage level begins with the pathogen-mediated damage level form, which then triggers the host-mediated damage level form. Currently, host damage can only be quantified in relative terms [13,14]. The term ‘relative’ is necessitated by the fact that, at present, host damage cannot be fully quantified because precise readouts of host damage remain limited and available mathematical modelling tools and computational platforms are insufficient for the quantification of host damage. When host damage surpasses a threshold that maintains host homeostasis, clinical disease occurs [14]. In the multiscale dynamics of the host-centred disease system form, the host damage level triggers the host immune response level [13,14].
[b.]
The host immune response level: This functional level of organisation of the host-centred disease system form also consists of two main complementary level forms, which are [i.] the innate immune response level form and [ii.] the adaptive immune response level form. Usually, the host immune response level begins with the innate immune response level form, which then triggers the adaptive host immune response level form. Each of these two complementary level forms of the host immune response level consists of three main scale mechanisms [15,16,17]: [i.] the molecular scale mechanism, which includes the proteome, lipodome, genome, metabolome, transcriptome, and complex molecular processes, such as gene expression, gene regulatory networks, signaling, and metabolic pathways involved in immunity and inflammation; [ii.] the cellular scale mechanism, which includes the activities and behaviour of the different immune cells, such as T-cells, B-cells, and different pathogen processes; and [iii.] the tissue scale mechanism, which includes inflammation mechanisms.
[c.]
The environmental change level: In the context of the universal theory for the multiscale modelling of infectious disease dynamics, this level of organisation of the host-centred disease system form arises because the environment subsystem is considered as an extended form of the host subsystem. This functional level of organisation of the host centred disease system form consists of two main complementary levels forms: [i.] the human-induced environmental-change-level form and [ii.] the naturally induced environmental-change-level form [18]. Usually, the environmental change level begins with the human-induced environmental-change-level form through mechanisms such as land use and population growth, which then triggers the naturally-induced environmental-change-level form such as extreme weather events, natural disasters, and climate change. Environmental change can significantly influence infectious disease dynamics through effects such as the survival and reproductive capacity of vectors and pathogens [19]. In the multiscale dynamics of the host-centred disease system form, the environmental change level triggers the biodiversity change level [20].
[d.]
The biodiversity change level: In the context of the universal theory for the multiscale modelling of infectious disease dynamics, the term biodiversity or biological diversity refers to the diversity of living organisms implicated in infectious disease dynamics, which are the pathogen subsystem and the host subsystem. This functional level of organisation of the host-centred disease system form also consists of two main complementary level forms: [i.] the pathogen-biodiversity-change-level form, which we also call the microbial-diversity-change-level form or microbiodiversity-change-level form [21], and [ii] the host-biodiversity-change-level form, which we also call the macrobiodiversity-change-level form [22,23]. In multiscale dynamics of the host-centred disease system form, biodiversity change occurs through various evolutionary mechanisms or adaptive mechanisms, which include [i.] mutation mechanisms, and/or [ii.] migration/dispersal mechanisms, and/or [iii.] genetic drift mechanisms, and/or [iv.] natural selection mechanisms. Usually, the biodiversity change level begins with the pathogen-biodiversity-change-level form, which then triggers the host-biodiversity-change-level form, resulting in the reciprocity of change occurring in both the pathogen subsystem and the host subsystem, whereby both the host subsystem and the pathogen subsystem impose evolutionary change on the other in a process involving co-evolutionary mechanisms. In this case, co-evolutionary mechanisms of the host subsystem and pathogen subsystem are a dynamic process of ongoing reciprocal change where a pathogen subsystem imposes an evolutionary influence on a host subsystem, which responds to the evolutionary pressure, in turn imposing an evolutionary influence on the pathogen subsystem, with this cycle potentially repeated over and over again. The outcome of this arms race may involve traits like parasite infectivity, host resistance, parasite host-finding ability, and parasite avoidance behaviour by the host [2]. The multiscale mechanisms of the biodiversity change level of the host-centred disease system form lie at the heart of the emergence and spread of pandemics, resulting in the emergence of new variants at the microbial-diversity-change-level form and new host species at the macrobiodiversity-change-level form.
These four main levels of organisation of the host-centred disease system form can be demarcated into two main groups. First, a group of primary levels of organisation of the host-centred disease system form consisting of the host damage level and the host immune response level. There is much stronger coupling between these two lower levels of organisation of the host-centred disease system form than with the higher levels of organisation of this disease system form since the host damage level triggers the host immune response level [13,14]. Second, a group of secondary levels of organisation of the host-centred disease system form consisting of the environmental change level and the biodiversity change level. There is also much stronger coupling between these two higher levels of organisation of the host-centred disease system form than with the lower levels of organisation of this disease system form; that is, the primary levels of organisation of the host-centred disease system form, since the environmental change level triggers the biodiversity change level [19,20,24]. The multiscale modelling of the host-centred disease system that incorporates the scales of both the primary levels of organisation and the secondary levels of organisation of the host-centred disease system form is at the core of the multiscale modelling of the ecology and evolution of infectious disease dynamics [22,23].

4.2. The Multiscale Organisation of the Host-Centred Disease System Form

Each level of organisation of the host-centred disease system form is dominated by multiple scale mechanisms that influence its dynamics at different scales. This is because, within each of the four main levels of organisation of the host-centred disease system form, which are [a.] the host damage level, [b.] the host immune response level, [c.] the environmental change level, and [d.] the biodiversity change level, a mechanism is considered compositional in the sense that the mechanism of each level of organisation as a whole can be broken down into organised interactions among the activities of different scale mechanisms at different scales of a level of organisation of this disease system form, with each scale associated with a particular scale mechanism. Based on this understanding, we establish that each of the four main levels of organisation of the host-centred disease system form is organised into seven main hierarchical scales of mechanism, which are: [a.] the molecular scale, [b.] the cell scale, [c.] the tissue scale, [d.] the organ scale, [e.] the whole organism scale, [f.] the community scale, and [g.] the ecosystem scale. Figure 5 is a conceptual representation of the seven main scales of a level of organisation of the host-centred disease system form.
Because of the different scale mechanisms at each of the four main levels of organisation, a host-centred disease system form is described by mechanism-based multiscale models. These mechanism-based multiscale models are described in terms of the following seven scale mechanisms at each of the four main levels of the host-centred disease system form:
[a.]
The molecular-scale mechanism: This involves molecular-scale mechanisms at each of the four main levels of organisation of the host-centred disease system form: the host damage level, the host immune response level, the environmental change level, and the biodiversity change level. Examples of molecular-scale mechanisms include mutation and genetic drift at the biodiversity change level [20], and cytokine, chemokine, and antibody mechanisms at the host immune response level [15,16,17]. Other examples of molecular mechanisms include nutrient cycle mechanisms such as nitrogen cycle and carbon cycle mechanisms at the environmental change level [25].
[b.]
The cellular-scale mechanism: This involves cellular=scale mechanisms at each of the four main levels of organisation of the host-centred disease system form: the host damage level, the host immune response level, the environmental change level, and the biodiversity change level. Examples of cellular-scale mechanisms include apoptosis and necrosis mechanisms at the host damage level [13], and T-cell, dendritic cell, and B-cell mechanisms, including macrophage mechanisms such as phagocytosis, at the host immune response level [15,16,17].
[c.]
The tissue-scale mechanism: This involves tissue-scale mechanisms at each of the four main levels of organisation of the host-centred disease system form: the host damage level, the host immune response level, the environmental change level, and the biodiversity change level. Examples of tissue-scale mechanisms include inflammation and fibrosis at the host immune response level [26].
[d.]
The organ-scale mechanism: This involves organ-scale mechanisms at each of the four main levels of organisation of the host-centred disease system form: the host damage level, the host immune response level, the environmental change level, and the biodiversity change level. Examples of organ-scale mechanisms include sepsis and gangrene formation [27] at the host damage level.
[e.]
The whole-organism-scale mechanism: This involves whole-organism-scale mechanisms at each of the four main levels of organisation of the host-centred disease system form: the host damage level, the host immune response level, the environmental change level, and the biodiversity change level. Examples of whole-organism-scale mechanisms include fever mechanisms and behavioural mechanisms such as hydrophobia at the host immune response level [13], and whole organism death at the host damage level.
[f.]
The community-scale mechanism: This involves community-scale mechanisms at each of the four main levels of organisation of the host-centred disease system form: the host damage level, the host immune response level, the environmental change level, and the biodiversity change level. Examples of community-scale mechanisms include adaptation and natural selection at the biodiversity change level and population extinction at the host damage level.
[g.]
The ecosystem-scale mechanism: This involves ecosystem-scale mechanisms at each of the four main levels of organisation of the host-centred disease system form: the host damage level, the host immune response level, the environmental change level, and the biodiversity change level. Examples of ecosystem-scale mechanisms include predation, competition, and the dilution effect—in which diverse host communities can reduce disease risk [28,29].
Therefore, mechanism-based multiscale models refer to multiscale models of the host-centred disease system form developed by incorporating multiple scale mechanisms that occur in the host subsystem due to the impact of the pathogen subsystem, which include [a.] molecular-scale mechanisms, [b.] cellular-scale mechanisms, [c.] tissue-scale mechanisms, [d.] organ-scale mechanisms, [e.] whole-organism-scale mechanisms, [f.] community-scale mechanisms, and [g.] ecosystem-scale mechanisms, in order to explain the multiscale temporal changes of the host-centred disease system form at each of its four main levels of organisation, which are [a.] the host damage level, [b.] the host immune response level, [c.] the environmental change level, and [d.] the biodiversity change level, taking into consideration how the pathogen subsystem causes these multiple scale mechanisms. The difficulty in the development of mechanism-based multiscale models of the host-centred disease system form when compared to the development of process-based multiscale models of the pathogen-centred disease system form is due to the following reasons: [a.] first, each level of the host-centred disease system form consists of two or more scales, with at most seven scales: [i.] the molecular scale, [ii.] cellular scale, [iii.] tissue scale, [iv.] organ scale, [v.] organism scale, [vi.] community scale, and [vii.] ecosystem scale]. This is unlike the pathogen-centred disease, where a level consists of only two scales, a macroscale and a microscale. [b.] Second, there is inter-dependence of levels, level forms, and scales with the simplest mechanism-based multiscale model required to incorporate at least the two primary levels of organisation of the host-centred disease system form: the host damage level and the host immune response level. However, in the pathogen-centred disease system form, only the two scales (i.e., the microscale and the macroscale) at a single level of this disease system form are assumed to be inter-dependent, while the levels are assumed to independent. [c.] Third, there are limitations imposed by observational and experimental tools to reveal the mechanisms at all the seven main scales of mechanism of each level of organisation of the host-centred disease system form, while the processes at each level of the pathogen-centred disease system form are well understood.

5. A Process-Based Multiscale of the Pathogen-Centred Disease System Form of Malaria Disease System as an Example

To illustrate ideas, we present a process-based multiscale model of the pathogen-centred disease system form of a malaria disease system. The universal theory of infectious disease dynamics demands both the process-based multiscale model of the pathogen-centred disease system form and the mechanism-based multiscale model of the host-centred disease system form to be integrated into a single multiscale model of infectious disease dynamics. However, it might not always be possible to account for both the pathogen-centred and the host-centred perspectives in the multiscale modelling of infectious disease dynamics because of current limitations in the availability of mathematical technology required to achieve that. But, the consideration of the possible contributions of both the pathogen-centred disease system form and the host-centred disease system form to the multiscale dynamics of an infectious disease system in line with the dictates of the universal theory for the multiscale modelling of infectious disease dynamics has the potential to provide a systems approach to its multiscale dynamics. The process-based multiscale model of the pathogen-centred disease system of a malaria disease system developed here integrates four sets of variables at the whole organism level, which are as follows:
[a.]
At the between-human scale, we have the following three variables: S H ( τ ) —susceptible human population size; I H ( τ ) —infected human population size; and P V ( τ ) —community sporozoite load.
[b.]
At the between-mosquito scale, we have the following three variables: S V ( τ ) —susceptible mosquito population size; I V ( τ ) —infected mosquito population size; and G H ( τ ) —community gametocyte load.
[c.]
At the within-human scale, we have the following four variables: R h ( t ) —susceptible erythrocytes (red blood cells); R m ( t ) —merozoite-infected erythrocytes; M h ( t ) —merozoites; and G h ( t ) —gametocyte-infected erythrocytes.
[d.]
At the within-mosquito scale, we have the following five variables: G v ( t ) —gametocyte-infected erythrocytes; G m ( t ) —sex cells called gametes; Z v ( t ) —zygotes; O v ( t ) —oocysts; and P v ( t ) —sporozoites.
Based on these fifteen variables, the coupled multiscale model of the pathogen-centred disease system form of a malaria disease system is given by the process-based multiscale model of scale order one (1). Details of the derivation of this multiscale model are given in Appendix A. Figure 6 shows a schematic representation of the coupled multiscale model of the pathogen-centred disease system form of a malaria disease system given by (1). The parameters of this coupled multiscale model are given and defined in Table A1 and Table A2 in Appendix A.
1 . d S H ( t ) d t = Λ H β V P V ( t ) S H ( t ) P 0 + P V ( t ) μ H S H ( t ) + γ H I H ( t ) , 2 . d I H ( t ) d t = β V P V ( t ) S H ( t ) P 0 + P V ( t ) μ H + γ H + δ H I H ( t ) , 3 . d P V ( t ) d t = P v ( t ) α v I V ( t ) + 1 α V P V ( t ) , 4 . d S V ( t ) d t = Λ V β H G H ( t ) S V ( t ) G 0 + G H ( t ) μ V S V ( t ) , 5 . d I V ( t ) d t = β H G H ( t ) S V ( t ) G 0 + G H ( t ) μ V + δ V I V ( t ) , 6 . d G H ( t ) d t = G h ( t ) α h I H ( t ) + 1 α H G H ( t ) , 7 . d G v ( t ) d t = β H G H ( t ) S V ( t ) 1 G 0 + G H ( t ) Φ V I V ( t ) + 1 α g + μ g G v ( t ) , 8 . d G m ( t ) d t = N g α g G v ( t ) α s + μ s G m ( t ) , 9 . d Z v ( t ) d t = ϕ s α s G m ( t ) α z + μ z Z v ( t ) , 10 . d O v ( t ) d t = α z Z v ( t ) α k + μ k O v ( t ) , 11 . d P v ( t ) d t = N k α k O v ( t ) α v + μ v P v ( t ) , 12 . d R h ( t ) d t = Λ h β h R h ( t ) M h ( t ) μ b R h ( t ) , 13 . d R m ( t ) d t = ( 1 ϕ h ) β h R h ( t ) M h ( t ) α m R m ( t ) , 14 . d M h ( t ) d t = β V P V ( t ) S H ( t ) 1 P 0 + P V ( t ) Φ H I H ( t ) + 1 + N m α m R m ( t ) μ m M h ( t ) , 15 . d G h ( t ) d t = ϕ h β h R h ( t ) M h ( t ) α h + μ h G h ( t ) .

6. Characterising the Multiscale Dynamics of the Pathogen-Centred Disease System Form of Malaria Disease System Using the Reproductive Number

In the context of the single-scale modelling of a malaria disease system at the whole organism level, the quantity R 0 quantifies the transmission process of this disease system at the macroscale only [30]. However, unlike the reproductive numbers derived from single-scale models of disease dynamics [30], which quantify the transmission process only, here, the reproductive number derived from the multiscale model (1) quantifies all four disease processes of the pathogen-centred disease system form of the malaria disease system—[a.] the super-infection by pathogen process, [b.] the pathogen replication process, [c.] the pathogen shedding or excretion process, and [d.] the pathogen transmission process—across the microscale and the macroscale. We use the next-generation operator approach to derive the basic reproduction number [31] from the process-based multiscale model (1). The process-based multiscale model (1) can be written in the form
1 . d X d t = f ( X , Y , Z ) , 2 . d Y d t = g ( X , Y , Z ) , 3 . d Z d t = h ( X , Y , Z ) ,
where
1 . X = S H , S V , R h , 2 . Y = I H , I V , G v , G m , Z v , O v , P v , R m , G h , 3 . Z = P V , G G , M h .
Components X denote the number of susceptibles, while components of Y represent the number of infected individuals that do not transmit the disease. Components of Z represent the number of individuals capable of transmitting the disease. Let
g ˜ ( X * , Z ) = [ g ˜ 1 ( X * , Z ) , g ˜ 2 ( X * , Z ) , g ˜ 3 ( X * , Z ) , g ˜ 4 ( X * , Z ) , g ˜ 5 ( X * , Z ) , g ˜ 6 ( X * , Z ) , g ˜ 7 ( X * , Z ) , g ˜ 8 ( X * , Z ) , g ˜ 9 ( X * , Z ) ] ,
with
1 . g ˜ 1 ( X * , Z ) = β V P V ( t ) Λ H μ H μ H + γ H + δ H P 0 + P V ( t ) , 2 . g ˜ 2 ( X * , Z ) = β H G H ( t ) Λ V μ V μ V + δ V G 0 + G H ( t ) , 3 . g ˜ 3 ( X * , Z ) = 1 α g + μ g . F V ( t ) μ V μ V + δ V , 4 . g ˜ 4 ( X * , Z ) = 1 α s + μ s . N g α g α g + μ g . F V ( t ) μ V μ V + δ V , 5 . g ˜ 5 ( X * , Z ) = 1 α z + μ z . ϕ s α s α s + μ s . N g α g α g + μ g . F V ( t ) μ V μ V + δ V , 6 . g ˜ 6 ( X * , Z ) = 1 α k + μ k . α z α z + μ n . ϕ s α s α s + μ s . N g α g α g + μ g . F V ( t ) μ V μ V + δ V , 7 . g ˜ 7 ( X * , Z ) = 1 α v + μ v . N k α k α k + μ k . α z α z + μ n . ϕ s α s α s + μ s . N g α g α g + μ g . F V ( t ) μ V μ V + δ V , 8 . g ˜ 8 ( X * , Z ) = ( 1 ϕ h ) β h M h Λ h α m μ b , 9 . g ˜ 9 ( X * , Z ) = ϕ h β h M h Λ h μ b α h + m u h ,
where
F V ( t ) = β H G H ( t ) ( Λ V μ V ) G 0 + G H ( t ) Φ V g ˜ 2 ( X * , Z ) + 1 ,
for Λ H > μ H and Λ V > μ V . We shall assume in all that follows, unless stated otherwise, that Λ H > μ H and Λ V > μ V .
Let
h ( X * , Z ) = h 1 = Ω v β H G H ( t ) ( Λ V μ V ) μ V μ V + δ V Φ V G 0 + G H ( t ) α V P V ( t ) , h 2 = ϕ h β h M h ( t ) Λ h α h [ g ˜ 1 ( X * , Z ) + 1 ] μ b α h + μ h α H G H ( t ) , h 3 = β V P V ( t ) ( Λ H μ H ) μ H P 0 + P V ( t ) Φ H [ g ˜ 1 ( X * , Z ) + 1 ] + N m α m ( 1 ϕ h ) β h M h ( t ) Λ h α m μ b μ m M h ( t )
where
Ω v = α v α v + μ v . N k α k α k + μ k . α z α z + μ z . ϕ s α s α s + μ s . N g α g α g + μ g .
Let A = D Z h ( X * , g ˜ ( X * , 0 ) , 0 ) and further assume that A can be written in the form A = M D , where M 0 and D > 0 , a diagonal matrix. Then, A becomes
A = α V Ω v β H ( Λ V μ V ) μ V μ V + δ V Φ V G 0 0 0 α H ϕ h β h Λ h α h μ b α h + μ h β V ( Λ H μ H ) Φ H μ H P 0 0 N m α m ( 1 ϕ h ) β h Λ h α m μ b μ m , ,
where Ω v is given by the expression (7). Since A = M D , we deduce matrices M and D to be
M = 0 Ω v β H ( Λ V μ V ) μ V μ V + δ V Φ V G 0 0 0 0 ϕ h β h Λ h α h μ b α h + μ h β V ( Λ H μ H ) Φ H μ H P 0 0 N m α m ( 1 ϕ h ) β h Λ h α m μ b , , D = α V 0 0 0 α H 0 0 0 μ m ,
where Ω v is given by the expression (7).
Then,
M D 1 = 0 Ω v β H ( Λ V μ V ) μ V μ V + δ V Φ V α H G 0 0 0 0 ϕ h β h Λ h α h μ b μ m α h + μ h β V ( Λ H μ H ) Φ H μ H α V P 0 0 N m α m ( 1 ϕ h ) β h Λ h α m μ b μ m , .
where Ω v is given by the expression (7).
The basic reproductive number is the spectral radius (dominant eigenvalue) of the matrix M D 1 ; that is,
R 0 = ρ ( M D 1 ) .
In this case, the reproductive number is obtained as the dominant eigenvalue from the cubic equation
λ 3 p λ 2 q = 0
Using the formula for solving cubic equations, the reproductive number can be shown to be
R 0 = p 3 + 27 q 2 p 2 54 + 27 q 2 p 2 54 2 p 6 729 2 3 + 27 q 2 p 2 54 27 q 2 p 2 54 2 p 6 729 2 3
where
p = N m α m ( 1 ϕ h ) β h Λ h α m μ b μ m q = Ω v β H ( Λ V μ V ) μ V μ V + δ V Φ V α H G 0 . β V ( Λ H μ H ) Φ H μ H α V P 0 . ϕ h β h Λ h α h μ b α h + μ h ,
and Ω v is given by the expression (7).
The reproductive number (12) derived from the coupled multiscale model (1) of the pathogen-centred disease system form of a malaria disease system indicates that it is a combination of all four main disease processes—[a.] the super-infection by pathogen process, [b.] the pathogen replication process, [c.] the pathogen shedding or excretion process, and [d.] the pathogen transmission process—across the microscale and the macroscale. This indicates that there is reciprocal influence between the microscale malaria disease dynamics and macroscale malaria disease dynamics for its pathogen-centred disease system form. Thus, the multiscale nature of the pathogen-centred disease system form of a malaria disease system described by the coupled multiscale model (1) is confirmed.

7. Characterising the Multiscale Dynamics of the Pathogen-Centred Disease System Form of Malaria Disease System Using Numerical Methods

In this section, we design a non-standard finite difference (NSFD) scheme that approximates the process-based multiscale model of scale order one (1). The non-standard finite difference (NSFD) scheme that we design for the process-based multiscale of a malaria disease system of scale order one (1) is based on the procedures and rules developed by Mickens [32,33,34,35]. Let X k = S H k , I H k , P V k , S V k , I V k , G H k , R h k , R m k , M h k , G h k T be the sequence of the solution, which must be non-negative and satisfy the underlying biological properties of the process-based multiscale model (1) for k N and the time step size h = Δ t . The NSFD scheme to solve the system (1) is designed so that it satisfies the conservation law property proposed by Mickens [32,35]. If we let λ V k = β V P H k P 0 + P H k and λ H k = β H G H k ( G 0 + G H k ) , then the NSFD scheme that approximates the process-based multiscale model (1) becomes:
1 . S H k + 1 = ϕ 1 ( h ) Λ H + S H k + ϕ 1 γ H I H k 1 + ϕ 1 ( h ) μ H + ϕ 1 ( h ) λ V k , 2 . I H k + 1 = ϕ 2 ( h ) λ V k S H k + 1 + I H k 1 + ϕ 2 ( h ) ( μ H + δ H + α H ) , 3 . S V k + 1 = ϕ 3 ( h ) Λ V + S V k 1 + ϕ 3 ( h ) μ V + ϕ 3 ( h ) * λ H k ) , 4 . I V k + 1 = ϕ 4 ( h ) λ H k S V k + 1 + I V k 1 + ϕ 4 ( h ) ( α V + δ V ) , 5 . P V k + 1 = ϕ 5 ( h ) I V k + 1 α h P v k + P V k 1 + ϕ 5 ( h ) α V , 6 . G H k + 1 = ϕ 6 ( h ) α h I V k + 1 + G H k 1 + ϕ 6 ( h ) α H , 7 . G v k + 1 = ϕ 7 ( h ) λ H k ( S V k + 1 1 ) Φ V ( I V k + 1 + 1 ) + G v k 1 + ϕ 7 ( h ) α g + μ g , 8 . G m k + 1 = ϕ 8 ( h ) N g α g G v k + 1 + G m k 1 + ϕ 8 ( h ) α s + μ s , 9 . Z v k + 1 = ϕ 9 ( h ) ϕ s α s G m k + 1 + G v k 1 + ϕ 9 ( h ) α z + μ z , 10 . O v k + 1 = ϕ 10 ( h ) α z Z v k + 1 + O v k 1 + ϕ 10 ( h ) α k + μ K , 11 . P v k + 1 = ϕ 11 ( h ) N k α k O v k + 1 + P v k 1 + ϕ 11 ( h ) α v + μ v , 12 . R h k + 1 = ϕ 12 ( h ) Λ h + G h k 1 + ϕ 12 ( h ) α m + | p h i 1 2 ( h ) M h k , 13 . R m k + 1 = ϕ 13 ( h ) ( 1 + ϕ h ) β h R h k + 1 M h k + R m k 1 + ϕ 13 ( h ) α m , 14 . G h k + 1 = ϕ 14 ( h ) ϕ h β h R h k + 1 M h k + G h k 1 + ϕ 14 ( h ) ( α h + μ h ) , 15 . M h k + 1 = ϕ 15 ( h ) λ V k ( S H k + 1 1 ) + M h k + N m α m R m k + 1 Φ H ( I H k + 1 + 1 ) ( 1 + ϕ 15 ( h ) μ m ) ,
where
  • ϕ 1 ( h ) = e h μ H 1 μ H ;
  • ϕ 2 ( h ) = e h ( δ H + δ H + μ H ) 1 δ H + δ H + μ H ;
  • ϕ 3 ( h ) = e h μ V 1 μ V ;
  • ϕ 4 ( h ) = e h ( α V + μ V ) 1 α V + μ V ;
  • ϕ 5 ( h ) = e h α V 1 α V ;
  • ϕ 6 ( h ) = e h α H 1 α H ;
  • ϕ 7 ( h ) = e h ( α g + μ g ) 1 α g + μ g ;
  • ϕ 8 ( h ) = e h ( α s + μ s ) 1 α s + μ s ;
  • ϕ 9 ( h ) = e h ( α z + μ z ) 1 α z + μ z ;
  • ϕ 10 ( h ) = e h ( α k + μ k ) 1 α k + μ k ;
  • ϕ 11 ( h ) = e h ( α v + μ v ) 1 α v + μ v ;
  • ϕ 12 ( h ) = e h μ h 1 μ h ;
  • ϕ 13 ( h ) = e h α m 1 α m ;
  • ϕ 14 ( h ) = e h ( α h + μ h ) 1 α h + μ h ;
  • ϕ 15 ( h ) = e h μ m 1 μ m .
The simulations for the process-based multiscale of a malaria disease system of scale order one (1) are performed using the NSFD scheme (14) and coded with MatLab. The parameters used in the simulations for the process-based multiscale of the pathogen-centred disease system form of a malaria disease system of scale order one (1) are summarised in Table A1 and Table A2 in Appendix A.
Figure 7 shows changes in (a) the community gametocyte load ( G H ), (b) population of infected mosquitoes ( I V ), (c) population of infected humans ( I H ), and (d) community sporozoite load ( P V ) for different values of the rate at which the gametocytes develop and become infectious to mosquitoes α h : α h = 0.02 , α h = 0.04 , α h = 0.08 . These results show that as that rate at which gametocytes develop and become infectious to mosquitoes at the within-human scale increases, so does the transmission rate of the pathogen-centred disease system form of the malaria disease system at both the between-human scale and the between-mosquito scale.
Figure 8 shows changes in (a) the community gametocyte load ( G H ), (b) population of infected mosquitoes ( I V ), (c) population of infected humans ( I H ), and (d) community sporozoite load ( P V ) for different values of the rate at which sporozoites become infectious to humans α v : α v = 0.025 , α v = 0.045 , α v = 0.085 . These results also show that as the rate at which sporozoites become infectious to humans α v increases at the within-human scale, so does the transmission rate of the pathogen-centred disease system form of the malaria disease system at both the between-human scale and the between-mosquito scale.
Figure 9 shows changes in (a) gametocyte-infected erythrocytes ( G h ), (b) the population of merozoites ( M h ), (c) merozoite-infected erythrocytes ( R m ), and (d) susceptible erythrocytes ( R h ) for different values of the contact rate of mosquitoes with the infectious reservoir of humans β H : β H = 0.00356 , β H = 0.0356 , β H = 0.356 . The results show that as the contact rate of mosquitoes with the infectious reservoir of humans β H increases at the macroscale, there is an insignificant increase in malaria parasite load at the within-human scale. This is expected because the malaria parasite has a replication cycle at the within-human scale that occurs at a very fast time scale. As a result, the malaria parasite load increase at the within-human scale is determined more by parasite replication than by the super-infection process.
Figure 10 shows changes in (a) gametocyte-infected erythtcytes ( G h ), (b) the population of merozoites ( M h ), (c) merozoite-infected erythrocytes ( R m ), and (d) susceptible erythrocytes ( R h ) for different values of contact rate of humans with the infectious reservoir of mosquitoes β V : β V = 0.003 , β V = 0.03 , β V = 0.3 . The results show that as the contact rate of humans with the infectious reservoir of mosquitoes β V increases at the macroscale, there is also an insignificant increase in malaria parasite load at the within-human scale. This is also expected because the malaria parasite has a replication cycle at the within-human scale that occurs at a very fast time scale. As a result, the malaria parasite load increase at the within-human scale is determined more by parasite replication than by the super-infection process.
Figure 11 shows changes in (a) sporozytes ( P v ), (b) oocysts ( O v ), (c) zygotes ( Z v ), and (d) gametocyte-infected erythrocytes ( G v ) for different values of the contact rate of humans with the infectious reservoir of mosquitoes β V : β V = 0.003 , β V = 0.03 , β V = 0.3 . The results show that as the contact rate of humans with the infectious reservoir of mosquitoes β V increases, there is also a significant increase in malaria parasite load at the within-mosquito scale. This also makes sense because the malaria parasite has no replication cycle at the within-mosquito scale. As a result, the malaria parasite load increase at the within-mosquito scale is determined more by the super-infection process of the within-mosquito scale.
Figure 12 shows changes in (a) sporozytes ( P v ), (b) oocysts ( O v ), (c) zygotes ( Z v ), and (d) gametocyte-infected erythrocytes ( G v ) for different values of the contact rate of mosquitoes with the infectious reservoir of humans β H : β H = 0.00356 , β H = 0.0356 , β H = 0.356 . The results show that as the contact rate of mosquitoes with the infectious reservoir of humans β H increases, there is a significant increase in malaria parasite load at the within-mosquito scale. This makes sense because the malaria parasite has no replication cycle at the within-mosquito scale. As a result, the malaria parasite load increase at the within-mosquito scale is determined more by the super-infection process of the within-mosquito scale.
Overall, all these numerical results show that unlike single-scale models of the pathogen-centred disease system form [30], which only consider the transmission process in malaria disease system dynamics, all four disease processes, which include [a.] the infection or super-infection by the pathogen process, [b.] the pathogen replication process, [c.] the pathogen excretion or shedding process, and [d.] the pathogen transmission process, that occur at the different scales; that is, the macroscale and microscale, of the pathogen-centred disease system form of malaria disease influence each in a reciprocal way for this disease system form. However, like all other modelling approaches, the process-based multiscale for malaria disease presented in this manuscript has its limitations. Further work carried out to improve the multiscale will require validation of the process-based multiscale model of the malaria disease system to account for uncertainties in data (clinical, experimental, demographic, and epidemiological) at both the microscale and microscale, as well as at the different temporal scales. The challenge is that there are no simple rules to follow when validating multiscale models of infectious disease dynamics. The parameters were not derived from a single epidemic situation. The demographic and epidemiology represent a general epidemic for developing countries. Therefore, the multiscale model presented models a general malaria epidemic situation for developing countries.

8. Discussion and Conclusions

The universal theory for multiscale modelling of infectious disease dynamics, which states that, for every host–pathogen interaction that results in an infectious disease system, there is no privileged or absolute scale of a disease system form that would determine the dynamics of the infectious disease system, only interactions between the scales of a level of organisation of the pathogen-centred disease system form and the scales of the corresponding levels of organisation of the host-centred disease system form, represents a historic opportunity to transform the way that we study the multiscale dynamics of infectious disease systems. It provides a wholly new paradigm for the multiscale modelling of infectious disease using mechanism-based multiscale models and process-based multiscale models. This new scientific theory revises and extends the replication-transmission relativity theory for the multiscale modelling of infectious disease dynamics [1] in a radical way and has great potential to transform mainstream thinking about the multiscale modelling of infectious systems. This is achieved by making possible the prospect that future progress in answering control, elimination, and even eradication questions of infectious disease systems is established on sound theoretical foundations using multiscale modelling methods. The theory explains the fact that an infectious disease system admits two different disease system form representations, which are the pathogen-centred disease system form and the host-centred disease system form. These two different disease system forms are considered to be mutually exclusive and complementary in their description of the multiscale dynamics of infectious disease systems. However, this scientific theory, as a conceptual breakthrough, is yet to reach its full potential in terms of scale integration and heuristic knowledge production once the mathematical and computational methods for the development of mechanism-based multiscale models of the host-centred disease system forms are fully established. Therefore, the research agenda for multiscale modelling as a complex systems science approach for the study of infectious disease dynamics in the coming years should place a greater emphasis on establishing the mathematical and computational methods for the development of mechanism-based multiscale models of the host-centred disease system form. We foresee the development of mechanism-based multiscale models of the host damage level and the biodiversity change level of the host-centred disease system form as grand challenges that will introduce gaps and bottlenecks along this future research path.
In order to illustrate the partial progress that has been achieved due to these grand challenges that introduce gaps and bottlenecks in the application of the universal theory for the multiscale modelling of infectious disease dynamics, we presented a coupled multiscale model of the pathogen-centred disease system form of a malaria disease system. For the pathogen-centred disease system form of the malaria disease system, the findings suggest that there are important medical planning consequences at the microscale to consider during the processes of epidemiological planning at the macroscale and vice versa in that a treatment that cures a patient of malaria at the microscale is equally good for both the patient and the general public because this patient no longer possesses a transmission risk for the malaria parasite towards both humans and mosquitoes at the macroscale.

Author Contributions

Conceptualization, W.G. methodology, W.G.; software, K.M.; validation, W.G. and K.M.; formal analysis, W.G.; investigation, K.M.; resources, W.G.; data curation, K.M.; writing—original draft preparation, W.G.; writing—review and editing, W.G. and K.M.; visualization, K.M.; supervision, W.G.; project administration, W.G.; funding acquisition, W.G. All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge with thanks financial support from South Africa National Research Foundation (NRF) Grant No. IPRR (UID 132 608).

Data Availability Statement

The authors declare that the data supporting the findings of this study are available within the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Derivation of the Process-Based Multiscale Model of the Pathogen-Centred Disease System Form of Malaria Disease System as an Example

The derivation of the process-based multiscale model of the pathogen-centred disease system form of a malaria disease system is an extension of a previously derived single-scale malaria model [30] and a previously derived coupled multiscale model of a malaria disease system [11]. Based on the taxonomic categorisation of process-based multiscale models into five main categories, which was first established in [2,36] and discussed further in this article in Section 3.2 as part of the revision of the replication-transmission relativity theory [1], coupled multiscale models are established by integrating multiscale models from the other four categories of multiscale models of the pathogen-centred disease system form. The coupled multiscale model in [11] of the malaria disease system was established by integrating a nested multiscale at the whole human level with an embedded multiscale model at the whole mosquito level. However, in this article, the coupled multiscale model is developed by integrating an embedded multiscale model at the whole human level with another embedded multiscale model at the whole mosquito level. This is because super-infection can occur in both the mosquito host and the human host. In an embedded multiscale model, the macroscale influences the microscale through super-infection; that is, repeated infection before the host recovers from the initial infectious episode [37]. From a mathematical point of view, the macroscale influences the microscale through macroscale variables and parameters in an embedded multiscale model [37]. However, in a nested multiscale, the macroscale influences the microscale through initial infective inoculum [1]. Mathematically, this means that the macroscale influences the microscale through initial conditions [7]. In the article [11], the derivation of the coupled multiscale model of the pathogen-centred disease system form of the malaria disease system was not fully presented, including its transformation from a multiscale model of scale order two to a multiscale of scale order one. In what follows, details of the derivation of the coupled multiscale model of the pathogen-centred disease system form of a malaria disease system are fully presented for the first time as a four-stage process as follows.

Appendix A.1. Stage I: Single-Scale Sub-models of the Process-Based Multiscale Model of the Pathogen-Centred Disease System Form of Malaria Disease System

At this stage, we present the derivation of the four single-scale sub-models of the pathogen-centred disease system form of a malaria disease system based on [11,30].
[a.]
The between-human scale sub-model: The development of a single-scale malaria sub-model at the macroscale of the whole human level of organisation of the pathogen-centred disease system form of a malaria disease system based on [30] involves establishing a new variable called community pathogen load P V ( τ ) . For the whole human level of organisation, the community pathogen load P V ( τ ) is the community sporozoite load. This community sporozoite load P V ( τ ) is a measure of the total infectious reservoir of mosquitoes in the community; that is, an aggregate population-level biomarker of a community’s sporozoite burden over a specific time period [11]. It is a useful metric for assessing the overall impact of malaria health interventions targeted at the mosquito vector or the uptake of malaria interventions targeted at the mosquito vector and quantifying their impact on the transmission of malaria from mosquitoes to humans as [30]: [i.] an indicator of a community’s level of infectiousness and transmission probability of malaria to humans, [ii.] a measure of the effectiveness of malaria interventions targeted at the mosquito vector, and [iii.] a proximal maker of malaria incidence among mosquitoes and their potential to propagate malaria to humans. Then, the single-scale malaria sub-model becomes an SIPS; that is, susceptible human hosts S H ( τ ) , infected human hosts I H ( τ ) , community sporozoite load P V ( τ ) , and susceptible human hosts S H ( τ ) , as follows:
1 . d S H ( τ ) d τ = Λ ˜ H β ˜ V P V ( τ ) S H ( τ ) P 0 + P V ( τ ) μ ˜ H S H ( τ ) + γ ˜ H I H ( τ ) , 2 . d I H ( τ ) d τ = β ˜ V P V ( τ ) S H ( τ ) P 0 + P V ( τ ) μ ˜ H + γ ˜ H + δ ˜ H I H ( τ ) , 3 . d P V ( τ ) d τ = N v α ˜ v I V ( τ ) α ˜ V P V ( τ ) .
The three variables of the between-human host scale sub-model (A1)
S H ( τ ) , I H ( τ ) , P V ( τ ) ,
vary at slow time scale τ . The units of the seven between-human host scale sub-model parameters
Λ ˜ H , β ˜ V , μ ˜ H , δ ˜ H , α ˜ v , γ ˜ H , α ˜ H ,
are per annum and have meanings as defined in Table A1. As indicated in [30], the expression
λ V ( P V ( τ ) ) = P V ( τ ) P 0 + P V ( τ ) ,
in the sub-model (A1) is the probability that a random bite by a mosquito vector in a particular community with a community sporozoite load P V ( τ ) will infect the human host with malaria in that community. Since the transmission rate of malaria (A4), which we also refer to as the infectivity response function of malaria, is a probability, it can be modelled by any function λ V ( P V ( τ ) ) with the specification that λ V : [ 0 , ) [ 0 , 1 ] . Further, since the function λ V ( P V ( τ ) ) is a probability, it must have the following properties [30]:
[i.]
Property I: The probability of infection vanishes in the absence of pathogen [i.e., λ V ( 0 ) = 0 ] and approaches 1 as the community sporozoite load becomes large [i.e., lim P V ( τ ) λ V ( P V ( τ ) ) = 1 ];
[ii.]
Property II: The probability of infection λ V ( P V ( τ ) ) increases with the community sporozoite load P V ( τ ) ; that is, λ V ( P V ( τ ) ) > 0 , where the prime denotes the derivative with respect to the argument.
Therefore, in the sub-model (A1), any function, λ V ( P V ( τ ) ) , that satisfies the above two properties can be used in place of the one derived from the Holling type I functional form (A4).
Table A1. Between-host (human and mosquito) parameter values and their description.
Table A1. Between-host (human and mosquito) parameter values and their description.
ParameterDescriptionParameter Value and RangeUnitsSource/Rational
α V Rate of elimination of community sporozoite load0.9 [0.09–0.99]day 1 Variable
Λ V Rate of supply of susceptible mosquitoes6000 [5000–7000]day 1 Variable
β V Contact rate of humans with the infectious reservoir of mosquitoes0.2 [0.1–0.5]day 1 [38]
G 0 Half-saturation constant for community gametocyte load 5 × 10 8 [ 1 × 10 8 10 × 10 8 ]day 1 Variable
μ V Natural death rate of mosquitoes0.12 [0.033–0.3]day 1 [39,40]
δ V Infection-induced death rate of mosquitoes0.00000426 [0.00000426–0.00000533]day 1 Assumed
Λ H Rate of supply of susceptible humans1000 [1000–2000]day 1 Assumed
β H Contact rate of mosquitoes with the infectious reservoir of humans0.3 [0.1–0.5]day 1 [41]
α H Rate of elimination of community gametocyte load0.0000913 [0.000467–0.000274]day 1 Variable
μ H Natural death rate of humans0.00004 [0.00001–0.00008]day 1 [39]
γ H Natural recovery rate of humans0.25 [0.1–0.5]day 1 Variable
P 0 Half-saturation constant for community sporozoite load 1 × 10 8 [ 1 × 10 7 5 × 10 8 ]day 1 Variable
δ H Disease-induced death rate of humans0.0027 [0.0001–0.5]day 1 Assumed
Φ H Proportion of new infected humans in the total infected human population0.001 [0.0001–0.5]NumberAssumed
Φ V Proportion of new infected mosquito vectors in the total infected mosquito vector population0.002 [0.0001–0.5]NumberAssumed
[b.]
The between-mosquito scale sub-model: Similarly, following [30], the community pathogen load G H ( τ ) for the whole mosquito level of organisation of the pathogen-centred disease system form of the malaria disease system is the community gametocyte load. This community gametocyte load G H ( τ ) is a measure of the total infectious reservoir of humans in the community; that is, an aggregate population-level biomarker of a community’s gametocyte burden over a specific time period [11]. It is a useful metric for assessing the overall impact of malaria health interventions targeted at the human host or the uptake of malaria interventions targeted at the human host and quantifying their impact on the transmission of malaria from humans to mosquitoes as [30]: [i.] an indicator of a community’s level of infectiousness and transmission probability of malaria to mosquitoes, [ii.] a measure of the effectiveness of malaria interventions targeted at the whole human host scale, and [iii.] a proximal maker of malaria incidence among humans and their potential to propagate malaria to mosquito vectors. Then, the single-scale malaria sub-model becomes an SIP; that is, susceptible mosquito hosts S V ( τ ) , infected mosquito hosts I V ( τ ) , and community gametocyte load G H ( τ ) , as follows:
1 . d S V ( t ) d τ = Λ ˜ V β ˜ H G H ( τ ) S V ( τ ) G 0 + G H ( τ ) μ ˜ V S V ( τ ) , 2 . d I V ( t ) d τ = β ˜ H G H ( τ ) S V ( τ ) G 0 + G H ( τ ) μ ˜ V + δ ˜ V I V ( τ ) , 3 . d G H ( τ ) d τ = N h α ˜ h I H ( τ ) α ˜ H G H ( τ ) .
The three variables of the between-mosquito host scale sub-model (A5)
S V ( τ ) , I V ( τ ) , G H ( τ ) ,
vary at slow time scale τ . The units of the six between-mosquito host scale sub-model parameters
Λ ˜ V , β ˜ H , μ ˜ V , δ ˜ V , α ˜ h , α ˜ V ,
are per annum and also have meanings as defined in Table A1.
Similarly, the expression
λ H ( G H ( τ ) ) = G H ( τ ) G 0 + G H ( τ ) ,
in the sub-model (A5) is the probability that a random bite of a human host by a mosquito vector in a particular community with a community gametocyte load G H ( τ ) will infect the mosquito host with malaria in that community. Since the transmission rate of malaria (A8), which we also refer to as the infectivity response function of malaria, is a probability, it must have the specification that λ H : [ 0 , ) [ 0 , 1 ] . Equally, since the function λ H ( G H ( τ ) ) is a probability, it must have the following properties [30]:
[i.]
Property I: The probability of infection vanishes in the absence of pathogen [i.e., λ H ( 0 ) = 0 ] and approaches 1 as the community gametocyte load becomes large [i.e., lim G H ( τ ) λ H ( G H ( τ ) ) = 1 ];
[ii.]
Property II: The probability of infection λ H ( G H ( τ ) ) increases with the community gametocyte load G H ( τ ) ; that is, λ H ( G H ( τ ) ) > 0 , where the prime denotes the derivative with respect to the argument.
Therefore, in the sub-model (A5), any function, λ H ( G H ( τ ) ) , that satisfies the above two properties can be used in place of those derived from the Holling type I functional form (A8).
[c.]
The within-human scale sub-model: Following [11], the within-human scale single-scale sub-model with variables R h ( t ) —susceptible erythrocytes (red blood cells); R m ( t ) —merozoite-infected erythrocytes; M h ( t ) —merozoites; and G h ( t ) —gametocyte-infected erythrocytes becomes:
1 . d R h ( t ) d t = Λ h β h R h ( t ) M h ( t ) μ b R h ( t ) , 2 . d R m ( t ) d t = ( 1 ϕ h ) β h R h ( t ) M h ( t ) α m R m ( t ) , 3 . d M h ( t ) d t = N m α m R m ( s ) μ m M h ( t ) , 4 . d G h ( t ) d t = ϕ h β h R h ( t ) M h ( t ) α h + μ h G h ( t ) .
The four variables of the within-human scale sub-model (A9)
R h ( t ) , R m ( t ) , M h ( t ) , G h ( t ) ,
vary at fast time scale t. The units of the eight within-human scale sub-model parameters
Λ h , β h , μ b , α m , N m , μ m , α h , μ h
are per day and have meanings as defined in Table A2.
Table A2. Within-mosquito and within-human parameter values and their description.
Table A2. Within-mosquito and within-human parameter values and their description.
ParameterDescriptionParameter Value and RangeUnitsSource/Rational
α g Rate at which gametocyte-infected erythrocytes burst96 [90–100]day 1 [42]
Λ v Rate of uptake of gametocytes through super-infection of mosquito300 [100–300]day 1 Variable
μ g Death rate of gametocytes0.0625 [0.0326–0.0725]day 1 [42]
N g Number of gametes produced per gametocyte-infected erythrocyte2 [1–3]day 1 Estimated
α s Fertilisation rate of gametes0.08 [0.01–0.2]no 1 day 1 [42]
μ s Natural decay rate of gametes58.0 [40–129]day 1 [42]
α z Rate at which zygotes develop into oocysts0.4240 [0.01–0.05]day 1 [43]
μ z Natural decay rate of zygotes1 [1–4]day 1 [42]
α k Bursting rate of oocysts to produce sporozoites0.2 [0.0–1.0]day 1 Variable
μ k Natural decay rate of oocysts0.01 [0.071–0.143]day 1 [44]
N k Number of sporozoites produced per bursting oocyst3000 [1000–10,000]day 1 [42]
α v Rate at which sporozoites become infectious to humans0.025 [0.167–1.00]day 1 [42]
μ v Natural decay rate of sporozoites0.0001 [0.0001–0.0]day 1 [42]
Λ h Rate of supply of susceptible red blood cells (erythrocytes)200 [100–300]day 1 [45]
β h Infection rate of erythrocytes by free merozoites0.1 [ 2 × 10 9 –0.2]day 1 [39,46]
μ b Natural decay rate of susceptible erythrocytes0.0083 [0.006–0.01]day 1 [46]
ϕ h Proportion of gametocyte-infected erythrocytes0.1 [0.1–0.5]NumberAssumed
μ m Natural decay rate of free merozoites0.001 [0.001–0.5]day 1 [39,46]
α m Rate at which erythrocytes burst to produce merozoites0.5 [0.1–7.0]day 1 [46,47]
N m Number of merozoites produced per bursting erythrocyte16 [10–30]day 1 [46]
α h Rate at which gametocytes develop and become infectious0.02 [0.01–0.9]day 1 [40]
μ h Natural decay rate of gametocyte-infected erythrocytes within infected humans0.0625 [0.0600–0.0625]day 1 [40]
  • The within-mosquito scale sub-model: Similarly, following [11], the within-mosquito scale sub-model with variables G v ( t ) —gametocyte infected erythrocytes; G m ( t ) —sex cells called gametes; Z v ( t ) —zygotes; O v ( t ) —oocysts; and P v ( t ) —sporozoites becomes:
    1 . d G v ( t ) d t = α g + μ g G v ( t ) , 2 . d G m ( t ) d t = N g α g G v ( t ) α s + μ s G m ( t ) , 3 . d Z v ( t ) d t = ϕ s α s G m ( t ) α z + μ z Z v ( t ) , 4 . d O v ( t ) d t = α z Z v ( t ) α k + μ k O v ( t ) , 5 . d P v ( t ) d t = N k α k O v ( t ) α v + μ v P v ( t ) .
    The main difference between the within-mosquito scale sub-model in [11] and the sub-model (A12) is that this sub-model does not incorporate super-infection in the mosquito vector. The five variables of the within-mosquito scale sub-model given by (A12)
    G v ( t ) , G m ( t ) , Z v ( t ) , O v ( t ) , P v ( t ) ,
    vary at fast time scale t. The units of the ten within-mosquito scale sub-model parameters
    α g , μ g , α p , μ s , α s , μ z , α z , μ k , α k , μ v , α v ,
    are per day and also have meanings as defined in Table A2.

Appendix A.2. Stage II: Integrating the Sub-models into a Single Process-Based Multiscale Model of the Pathogen-Centred Disease System Form of Malaria Disease System of Scale Order Two

The four sub-models of the pathogen-centred disease system form of the malaria disease system can be integrated into a coupled multiscale model by up-scaling individual infectiousness [ P v ( t ) , G h ( t ) ] to population infectiousness [ P V ( τ ) , G H ( τ ) ] as illustrated in Figure 3. This is achieved by replacing N v with P v ( t ) in the between-human scale submodel (A1) and N h with G h ( t ) in the between-mosquito scale submodel (A5). The integrated process-based multiscale model of the pathogen-centred disease system form of the malaria disease system is given by (A15) as follows:
1 . d S H ( τ ) d τ = Λ ˜ H β ˜ V P V ( τ ) S H ( τ ) P 0 + P V ( τ ) μ ˜ H S H ( τ ) + γ ˜ H I H ( τ ) , 2 . d I H ( τ ) d τ = β ˜ V P V ( τ ) S H ( τ ) P 0 + P V ( τ ) μ ˜ H + γ ˜ H + δ ˜ H I H ( τ ) , 3 . d P V ( τ ) d τ = P v ( t ) α ˜ v I V ( τ ) α ˜ V P V ( τ ) , 4 . d S V ( τ ) d τ = Λ ˜ V β ˜ H G H ( τ ) S V ( τ ) G 0 + G H ( τ ) μ ˜ V S V ( τ ) , 5 . d I V ( τ ) d τ = β ˜ H G H ( τ ) S V ( τ ) G 0 + G H ( τ ) μ ˜ V + δ ˜ V I V ( τ ) , 6 . d G H ( τ ) d τ = G h ( t ) α ˜ h I H ( τ ) α ˜ H G H ( τ ) , 7 . d G v ( t ) d t = α g + μ g G v ( t ) , 8 . d G m ( t ) d t = N g α g G v ( t ) α s + μ s G m ( t ) , 9 . d Z v ( t ) d t = ϕ s α s G m ( t ) α z + μ z Z v ( t ) , 10 . d O v ( t ) d t = α z Z v ( t ) α k + μ k O v ( t ) , 11 . d P v ( t ) d t = N k α k O v ( t ) α v + μ v P v ( t ) , 12 . d R h ( t ) d t = Λ h β h R h ( t ) M h ( t ) μ b R h ( t ) , 13 . d R m ( t ) d t = ( 1 ϕ h ) β h R h ( t ) M h ( t ) α m R m ( t ) , 14 . d M h ( t ) d t = N m α m R m ( t ) μ m M h ( t ) , 15 . d G h ( t ) d t = ϕ h β h R h ( t ) M h ( t ) α h + μ h G h ( t ) .
This is a coupled multiscale model of scale order two that integrates a nested multiscale model at the whole human level with another nested multiscale model at the whole mosquito level.

Appendix A.3. Stage III: Scale Order Reduction of the Coupled Multiscale Model of the Pathogen-Centred Disease System Form of Malaria Disease System from a Multiscale Model of Scale Order Two to a Multiscale Model of Scale Order One

In general, the scales; that is, the microscale and the macroscale, of a level of organisation of the pathogen-centred disease form are hierarchically organised in both space and time. Therefore, any process-based multiscale model of the pathogen-centred disease form that incorporates the hierarchy in spatial scales of the microscale and the macroscale and the hierarchy in the times scales of the processes that occur at the microscale and the macroscale, is said to be a process-based multiscale model of scale order two, which is the highest possible scale order for such process-based multiscale models. The scale order of two of a process-based multiscale model can be reduced to scale order one, which is the lowest possible scale order of such process-based multiscale models, by re-casting the multiscale model into a regular perturbation problem or singular perturbation problem [48] and then performing a slow and fast time scale analysis. This reduction in scale order of a process-based multiscale model relies on there being a dimensionless and relatively small parameter ϵ in the process-based multiscale model such that 0 < ϵ 1 . In such a reduced scale order process-based multiscale model, which does not incorporate the hierarchy of time scales of the processes that occur at the microscale and macroscale, disease processes at both the microscale and the macroscale occur at the same time scale. For an example of a scale order reduction of a process-based multiscale of scale order two to a process-based multiscale of scale order one, see [11] for the pathogen-centred disease system form of a malaria disease system by re-casting the process-based multiscale model into a singular perturbation problem. In this study, the scale order reduction is achieved by recasting the process-based multiscale model system (A15) into a regular perturbation problem by expressing the slow time scale of the processes at the macroscale and the fast time scale of the processes at the microscale in terms of each other using the relationship τ = ϵ t , where 0 < ϵ 1 and ϵ is a constant highlighting the fast time scale t of the within-host scale sub-model compared to the slow time scale τ of the between-host scale sub-model. By performing a fast–slow time scale analysis of the regular perturbation problem, the multiscale model of scale order two (A15) becomes a multiscale model of scale order one (A16) as follows.
1 . d S H ( t ) d t = ϵ Λ ˜ H ϵ β ˜ V P V ( t ) S H ( t ) P 0 + P V ( t ) ϵ μ ˜ H S H ( t ) + ϵ γ ˜ H I H ( t ) , 2 . d I H ( t ) d t = ϵ β ˜ V P V ( t ) S H ( t ) P 0 + P V ( t ) ϵ μ ˜ H + ϵ γ ˜ H + ϵ δ ˜ H I H ( t ) , 3 . d P V ( t ) d t = P v ( t ) ϵ α ˜ v I V ( t ) ϵ α ˜ V P V ( t ) , 4 . d S V ( t ) d t = ϵ Λ ˜ V ϵ β ˜ H G H ( t ) S V ( t ) G 0 + G H ( t ) S V ( t ) ϵ μ ˜ V , 5 . d I V ( t ) d t = ϵ β ˜ H G H ( t ) S V ( t ) G 0 + G H ( t ) ϵ μ ˜ V + ϵ δ ˜ V I V ( t ) , 6 . d G H ( τ ) d τ = G h ( t ) ϵ α ˜ h I H ( τ ) ϵ α ˜ H G H ( τ ) , 7 . d G v ( t ) d t = α g + μ g G v ( t ) , 8 . d G m ( t ) d t = N g α g G v ( t ) α s + μ s G m ( t ) , 9 . d Z v ( t ) d t = ϕ s α s G m ( t ) α z + μ z Z v ( t ) , 10 . d O v ( t ) d t = α z Z v ( t ) α k + μ k O v ( t ) , 11 . d P v ( t ) d t = N k α k O v ( t ) α v + μ v P v ( t ) , 12 . d R h ( t ) d t = Λ h β h R h ( t ) M h ( t ) μ b R h ( t ) , 13 . d R m ( t ) d t = ( 1 ϕ h ) β h R h ( t ) M h ( t ) α m R m ( t ) , 14 . d M h ( t ) d t = N m α m R m ( t ) μ m M h ( t ) , 15 . d G h ( t ) d t = ϕ h β h R h ( t ) M h ( t ) α h + μ h G h ( t ) .
In (A16), the process-based multiscale model is re-cast as a regular perturbation problem. For the process-based multiscale model (A16), choose ϵ = 1 / 365 so that
ϵ Λ ˜ H = Λ ˜ H 365 = Λ H , ϵ β ˜ H = β ˜ H 365 = β H , ϵ μ ˜ H = μ ˜ H 365 = μ H , ϵ δ ˜ H = δ ˜ H 365 = δ H , ϵ Λ ˜ V = Λ ˜ V 365 = Λ V , ϵ β ˜ V = β ˜ V 365 = β V , ϵ μ ˜ H = μ ˜ H 365 = μ V , ϵ δ ˜ V = δ ˜ H 365 = δ V , ϵ α ˜ h = α ˜ h 365 = α h , ϵ α ˜ v = α ˜ v 365 = α v , ϵ α ˜ V = α ˜ V 365 = α V , ϵ α ˜ H = α ˜ H 365 = α H ϵ γ ˜ H = γ ˜ H 365 = γ H .
The thirteen between-host scale parameters given by (A17) of the coupled multiscale model (A16) have units of per annum. The effect of multiplying by 1 / 365 is that the transformed new parameters
Λ H , β V , μ H , δ H , Λ V , β H , μ V , δ V , α h , γ H , α H , α v , α h .
now have units of per day. These units are the same as the units for the within-host scale parameters. The coupled multiscale model of the pathogen-centred disease system form of the malaria disease system of scale order two (A15) becomes a coupled multiscale model of scale order one (A19) as follows:
1 . d S H ( t ) d t = Λ H β V P V ( t ) S H ( t ) P 0 + P V ( t ) μ H S H ( t ) + γ H I H ( t ) , 2 . d I H ( t ) d t = β V P V ( t ) S H ( t ) P 0 + P V ( t ) μ H + γ H + δ H I H ( t ) , 3 . d P V ( t ) d t = P v ( t ) α v I V ( t ) α V P V ( t ) , 4 . d S V ( t ) d t = Λ V β H G H ( t ) S V ( t ) G 0 + G H ( t ) μ V S V ( t ) , 5 . d I V ( t ) d t = β H G H ( t ) S V ( t ) G 0 + G H ( t ) μ V + δ V I V ( t ) , 6 . d G H ( t ) d t = G h ( t ) α h I H ( t ) α H G H ( t ) , 7 . d G v ( t ) d t = α g + μ g G v ( t ) , 8 . d G m ( t ) d t = N g α g G v ( t ) α s + μ s G m ( t ) , 9 . d Z v ( t ) d t = ϕ s α s G m ( t ) α z + μ z Z v ( t ) , 10 . d O v ( t ) d t = α z Z v ( t ) α k + μ k O v ( t ) , 11 . d P v ( t ) d t = N k α k O v ( t ) α v + μ v P v ( t ) , 12 . d R h ( t ) d t = Λ h β h R h ( t ) M h ( t ) μ b R h ( t ) , 13 . d R m ( t ) d t = ( 1 ϕ h ) β h R h ( t ) M h ( t ) α m R m ( t ) , 14 . d M h ( t ) d t = N m α m R m ( t ) μ m M h ( t ) , 15 . d G h ( t ) d t = ϕ h β h R h ( t ) M h ( t ) α h + μ h G h ( t ) .
Unlike the reduced scale order multiscale model (A19), scale order reduction by re-casting the process-based multiscale into a singular perturbation problem as in [11] has the disadvantage of reducing the dimension of the process-based multiscale model and therefore reducing its accuracy. For the coupled multiscale model of scale order one (A19), the dynamics at the microscale and the macroscale vary at the same time scale, t, and the parameters are also expressed in the same units, which are per day. However, the coupled multiscale (A19) can be extended to incorporate the more realistic features of super-infection in both the human host and the mosquito host; that is, repeated infection before the host recovers from the initial infectious episode. In what follows, we convert the coupled multiscale model of scale order one (A19) into a corresponding coupled multiscale model of scale order one with super-infection in both the human host and the mosquito host.

Appendix A.4. Stage IV: Extending the Coupled Multiscale Model of the Pathogen-Centred Disease System Form of Malaria Disease System by Incorporating Variable Super-Infection in Both the Mosquito Host and the Human Host

The necessary mathematical basis for incorporating super-infection in a process-based multiscale model was established in [49]. The approach involves the down-scaling and up-scaling of variables and parameters at both the within-host scale and between-host scale as illustrated in Figure 3. In the context of the process-based multiscale model of the pathogen-centred disease system form of a malaria disease system, this can be achieved as follows:
[a.]
For the Super-infection of the Human Host: By down-scaling and up-scaling disease process variables and parameters across the microscale and the macroscale as illustrated in Figure 3, we have the following additional expressions:
[i.]
Influence of macroscale on microscale: The influence of the between-whole mosquito scale (the macroscale for the whole mosquito level) on the within-whole human scale (the microscale for the whole human level) is modeled by down-scaling the uptake of the pathogen in the community through the infection of humans at the between-human scale at a rate β V P V ( t ) S H ( t ) P 0 + P V ( t ) to the repeated infection of the within-human scale at a rate β V P V ( t ) S H ( t ) 1 P 0 + P V ( t ) Φ H I H ( t ) + 1 . However, this representation of super-infection is a refinement of the approach in [49] in the following way. The approach in [49] over-estimates the number of new infections, while, here, the number of new infections is assumed to be a proportion Φ H of the existing cumulative number of infections.
[ii.]
Influence of microscale on macroscale: The influence of the within-whole human scale (the microscale for the whole human level) on the between-whole human scale (the macroscale for the whole human level) is modeled by up-scaling the individual human excretion/shedding of the pathogen at the within-whole human scale at a rate α h G h to the between-whole mosquito scale at a rate I H ( t ) + 1 α h G h ( t ) .
[b.]
For the Super-infection of the Mosquito Host: Equally, by down-scaling and up-scaling disease process variables and parameters across the microscale and the macroscale as illustrated in Figure 3, we have the following additional expressions:
[i.]
Influence of macroscale on microscale: The influence of the between-whole human scale (the macroscale for the whole human level) on the within-whole mosquito scale (the microscale for the whole mosquito level) is modeled by down-scaling the uptake of the pathogen in the community through the infection of mosquitoes at the between-mosquito scale at a rate β H G H ( t ) S V ( t ) G 0 + G H ( t ) to the repeated infection of the within-mosquito scale at a rate β H G H ( t ) S V ( t ) 1 G 0 + G H ( t ) Φ V I V ( t ) + 1 . However, this representation of super-infection is also a refinement of the approach in [49] in the following way. The approach in [49] over-estimates the number of new infections, while, here, the number of new infections is assumed to be a proportion Φ V of the existing cumulative number of infections.
[ii.]
Influence of microscale on macroscale: The influence of the within-whole mosquito scale (the microscale for the whole mosquito level) on the between-whole human scale (the macroscale for the whole human level) is modeled by up-scaling the individual mosquito excretion/shedding of the pathogen at the within-whole mosquito scale at a rate α v P v to the between-whole human scale at a rate I V ( t ) + 1 α v P v ( t ) .
We obtain the coupled multiscale model of scale order one (A20) of the pathogen-centred disease system form of a malaria disease system with super-infection in both the human host and the mosquito host as follows.
1 . d S H ( t ) d t = Λ H β V P V ( t ) S H ( t ) P 0 + P V ( t ) μ H S H ( t ) + γ H I H ( t ) , 2 . d I H ( t ) d t = β V P V ( t ) S H ( t ) P 0 + P V ( t ) μ H + γ H + δ H I H ( t ) , 3 . d P V ( t ) d t = P v ( t ) α v I V ( t ) + 1 α V P V ( t ) , 4 . d S V ( t ) d t = Λ V β H G H ( t ) S V ( t ) G 0 + G H ( t ) μ V S V ( t ) , 5 . d I V ( t ) d t = β H G H ( t ) S V ( t ) G 0 + G H ( t ) μ V + δ V I V ( t ) , 6 . d G H ( t ) d t = G h ( t ) α h I H ( t ) + 1 α H G H ( t ) , 7 . d G v ( t ) d t = β H G H ( t ) S V ( t ) 1 G 0 + G H ( t ) Φ V I V ( t ) + 1 α g + μ g G v ( t ) , 8 . d G m ( t ) d t = N g α g G v ( t ) α s + μ s G m ( t ) , 9 . d Z v ( t ) d t = ϕ s α s G m ( t ) α z + μ z Z v ( t ) , 10 . d O v ( t ) d t = α z Z v ( t ) α k + μ k O v ( t ) , 11 . d P v ( t ) d t = N k α k O v ( t ) α v + μ v P v ( t ) , 12 . d R h ( t ) d t = Λ h β h R h ( t ) M h ( t ) μ b R h ( t ) , 13 . d R m ( t ) d t = ( 1 ϕ h ) β h R h ( t ) M h ( t ) α m R m ( t ) , 14 . d M h ( t ) d t = β V P V ( t ) S H ( t ) 1 P 0 + P V ( t ) Φ H I H ( t ) + 1 + N m α m R m ( t ) μ m M h ( t ) , 15 . d G h ( t ) d t = ϕ h β h R h ( t ) M h ( t ) α h + μ h G h ( t ) .
The study of the pathogen-centred disease system form of a malaria disease system using the coupled multiscale model (A20) is presented in Section 5 and Section 6.

References

  1. Garira, W. The Replication-Transmission Relativity theory for Multiscale Modelling of infectious Disease Systems. Sci. Rep. 2019, 9, 1–17. [Google Scholar] [CrossRef] [PubMed]
  2. Garira, W. A primer on multiscale modelling of infectious disease systems. Infect. Dis. Model. 2018, 3, 176–191. [Google Scholar] [CrossRef] [PubMed]
  3. Garira, W. The research and development process for multiscale models of infectious disease systems. PLoS Comput. Biol. 2020, 16, e1007734. [Google Scholar] [CrossRef] [PubMed]
  4. Hunter, E.; Mac Namee, B.; Kelleher, J.D. A Taxonomy for Agent-Based Models in Human Infectious Disease Epidemiology. J. Artif. Soc. Soc. Simul. 2017, 20, 1–2. [Google Scholar] [CrossRef]
  5. Keeling, M.J.; Eames, K.T. Networks and epidemic models. J. R. Soc. Interface 2005, 2, 295–307. [Google Scholar] [CrossRef]
  6. Chhetri, B.; Vamsi, D.K.; Sanjeevi, C. A Nested Multi-Scale Model for COVID-19 Viral Infection. arXiv 2021, arXiv:2108.12150. [Google Scholar]
  7. Netshikweta, R.; Garira, W. A nested multiscale model to study paratuberculosis in ruminants. Front. Appl. Math. Stat. 2022, 79, 8762. [Google Scholar] [CrossRef]
  8. Mideo, N.; Alizon, S.; Day, T. Linking within-and between-host dynamics in the evolutionary epidemiology of infectious diseases. Trends Ecol. Evol. 2008, 23, 511–517. [Google Scholar] [CrossRef]
  9. Quintela, B.M.; Conway, J.M.; Hyman, J.M.; Guedj, J.; Dos Santos, R.W.; Lobosco, M.; Perelson, A.S. A New Age-Structured Multiscale Model of the Hepatitis C Virus Life-Cycle During Infection and Therapy with Direct-Acting Antiviral Agents. Front. Microbiol. 2018, 9, 601. [Google Scholar] [CrossRef]
  10. Zitzmann, C.; Kaderali, L. Mathematical Analysis of Viral Replication Dynamics and Antiviral Treatment Strategies: From Basic Models to Age-Based Multi-Scale Modeling. Front. Microbiol. 2018, 9, 1546. [Google Scholar] [CrossRef]
  11. Garira, W.; Mathebula, D. Development and application of multiscale models of acute viral infections in intervention research. Math. Methods Appl. Sci. 2020, 43, 3280–3306. [Google Scholar] [CrossRef]
  12. Garira, W.; Chirove, F. A general method for multiscale modelling of vector-borne disease systems. Interface Focus 2020, 10, 20190047. [Google Scholar] [CrossRef]
  13. Casadevall, A.; Pirofski, L.A. The damage-response framework of microbial pathogenesis. Nat. Rev. Microbiol. 2003, 1, 17–24. [Google Scholar] [CrossRef]
  14. Pirofski, L.A.; Casadevall, A. The damage-response framework of microbial pathogenesis and infectious diseases. GI Microbiota Regul. Immune Syst. 2008, 4, 135–146. [Google Scholar]
  15. Waites, W.; Cavaliere, M.; Danos, V.; Datta, R.; Eggo, R.M.; Hallett, T.B.; Manheim, S.; Griffiths, J.P.; Russell, T.; Zarnitsyna, V.I. Compositional modelling of immune response and virus transmission dynamics. Philos. Trans. R. Soc. A 2022, 380, 20210307. [Google Scholar] [CrossRef] [PubMed]
  16. Holding, T.; Valletta, J.J.; Recker, M. Multiscale immune selection and the transmission-diversity feedback in antigenically diverse pathogen systems. Am. Nat. 2018, 192, E189–E201. [Google Scholar] [CrossRef] [PubMed]
  17. Cappuccio, A.; Tieri, P.; Castiglione, F. Multiscale modelling in immunology: A review. Brief. Bioinform. 2016, 17, 408–418. [Google Scholar] [CrossRef] [PubMed]
  18. Budria, A.; Candolin, U. How does human-induced environmental change influence host-parasite interactions? Parasitology 2014, 141, 462–474. [Google Scholar] [CrossRef]
  19. Johnson, P.T.; Wood, C.L.; Joseph, M.B.; Preston, D.L.; Haas, S.E.; Springer, Y.P. Habitat heterogeneity drives the host-diversity-begets-parasite-diversity relationship: Evidence from experimental and field studies. Ecol. Lett. 2016, 19, 752–761. [Google Scholar] [CrossRef]
  20. Brunner, F.S.; Eizaguirre, C. Can environmental change affect host/parasite-mediated speciation? Zoology 2016, 119, 384–394. [Google Scholar] [CrossRef]
  21. Gudelj, I.; Weitz, J.S.; Ferenci, T.; Claire Horner-Devine, M.; Marx, C.J.; Meyer, J.R.; Forde, S.E. An integrative approach to understanding microbial diversity: From intracellular mechanisms to community structure. Ecol. Lett. 2010, 13, 1073–1084. [Google Scholar] [CrossRef] [PubMed]
  22. Ostfeld, R.S.; Keesing, F. Effects of host diversity on infectious disease. Annu. Rev. Ecol. Evol. Syst. 2012, 43, 2012. [Google Scholar] [CrossRef]
  23. Johnson, P.T.; Preston, D.L.; Hoverman, J.T.; LaFonte, B.E. Host and parasite diversity jointly control disease risk in complex communities. Proc. Natl. Acad. Sci. USA 2013, 110, 16916–16921. [Google Scholar] [CrossRef] [PubMed]
  24. Halliday, F.W.; Rohr, J.R. Measuring the shape of the biodiversity-disease relationship across systems reveals new findings and key gaps. Nat. Commun. 2019, 10, 5032. [Google Scholar] [CrossRef] [PubMed]
  25. Johnson, P.T.; Townsend, A.R.; Cleveland, C.C.; Glibert, P.M.; Howarth, R.W.; McKenzie, V.J.; Rejmankova, P.; Ward, M.H. Linking environmental nutrient enrichment and disease emergence in humans and wildlife. Ecol. Appl. 2010, 20, 16–29. [Google Scholar] [CrossRef]
  26. Gangell, C.; Gard, S.; Douglas, T.; Park, J.; De Klerk, N.; Keil, T.; Brennan, A.; Ranganathan, S.; Browne, R.; Sly, P.D. Inflammatory responses to individual microorganisms in the lungs of children with cystic fibrosis. Clin. Infect. Dis. 2011, 53, 425–432. [Google Scholar] [CrossRef]
  27. Jiang, J.L.; Tseng, L.W.; Chang, H.R. Symmetrical peripheral gangrene in sepsis after treatment with inotropes. Tzu-Chi Med. J. 2017, 29, 121. [Google Scholar] [CrossRef]
  28. Luis, A.D.; Kuenzi, A.J.; Mills, J.N. Species diversity concurrently dilutes and amplifies transmission in a zoonotic host–pathogen system through competing mechanisms. Proc. Natl. Acad. Sci. USA 2018, 115, 7979–7984. [Google Scholar] [CrossRef]
  29. Halliday, F.W.; Rohr, J.R.; Laine, A.L. Biodiversity loss underlies the dilution effect of biodiversity. Ecol. Lett. 2020, 23, 1611–1622. [Google Scholar] [CrossRef]
  30. Garira, W.; Maregere, B. The transmission mechanism theory of disease dynamics: Its aims, assumptions and limitations. Infect. Dis. Model. 2023, 8, 122. [Google Scholar] [CrossRef]
  31. Castillo-Chavez, C.; Feng, Z.; Huang, W. On the computation of R0 and its role in global stability. In Mathematical Approaches for Emerging and Re-Emerging Infectious Diseases Part 1: An Introduction to Models, Methods and Theory; Castillo-Chavez, C., Blower, S., Driessche, P.V., Kirschner, D., Eds.; The IMA Volumes in Mathematics and Its Applications; Springer: Berlin/Heidelberg, Germany, 2002; Volume 125, pp. 229–250. [Google Scholar]
  32. Mickens, R.E. Nonstandard Finite Difference Models of Differential Equations; World Scientific: Singapore, 1994. [Google Scholar]
  33. Mickens, R.E. Advances in the Applications of Nonstandard Finite Difference Schemes; World Scientific: Singapore, 2005. [Google Scholar]
  34. Mickens, R.E. Nonstandard finite difference schemes for differential equations. J. Differ. Equ. Appl. 2002, 8, 823–847. [Google Scholar] [CrossRef]
  35. Mickens, R.E. Dynamic consistency: A fundamental principle for constructing nonstandard finite difference schemes for differential equations. J. Differ. Equ. Appl. 2005, 11, 645–653. [Google Scholar] [CrossRef]
  36. Garira, W. A complete categorization of multiscale models of infectious disease systems. J. Biol. Dyn. 2017, 11, 378–435. [Google Scholar] [CrossRef]
  37. Netshikweta, R.; Garira, W. An Embedded Multiscale Modelling to Guide Control and Elimination of Paratuberculosis in Ruminants. Comput. Math. Methods Med. 2021, 2021, 9919700. [Google Scholar] [CrossRef] [PubMed]
  38. Shapiro, L.L.M.; Whitehead, S.A.; Thomas, M.B. Quantifying the effects of temperature on mosquito and parasite traits that determine the transmission potential of human malaria. PLoS Biol. 2017, 15, e2003489. [Google Scholar] [CrossRef] [PubMed]
  39. Cai, L.; Tuncer, N.; Martcheva, M. How does within-host dynamics affect population level dynamics? Insights from an immuno-epidemiological model of malaria. Math. Methods Appl. Sci. 2017, 40, 6424–6450. [Google Scholar] [CrossRef]
  40. Legros, M.; Bonhoeffer, S. A combined within-host and between-hosts modelling framework for the evolution of resistance to antimalarial drugs. J. R. Soc. Interface 2016, 13, 148. [Google Scholar] [CrossRef]
  41. Smith, D.L.; McKenzie, F.E. Statics and dynamics of malaria infection in Anopheles mosquitoes. Malar. J. 2004, 3, 1–14. [Google Scholar] [CrossRef]
  42. Teboh-Ewungkem, M.I.; Yuster, T. A within-vector mathematical model of Plasmodium falciparum and implications of incomplete fertilization on optimal gametocyte sex ratio. J. Theor. Biol. 2010, 264, 273–286. [Google Scholar] [CrossRef]
  43. Gouagna, L.C.; Bonnet, S.; Gounoue, R.; Verhave, J.P.; Eling, W.; Sauerwein, R.; Boudin, C. Stage-specific effects of host plasma factors on the early sporogony of autologous Plasmodium falciparum isolates within Anopheles gambiae. Trop. Med. Int. Health 2004, 9, 937–948. [Google Scholar] [CrossRef]
  44. Howells, R.E.; Davies, E.E. Nuclear division in the oocyst of Plasmodium berghei. Ann. Trop. Med. Parasitol. 1971, 65, 451–459. [Google Scholar] [CrossRef] [PubMed]
  45. Tabo, Z.; Luboobi, L.S.; Ssebuliba, J. Mathematical modelling of the in-host dynamics of malaria and the effects of treatment. J. Math. Comput. Sci. 2017, 17, 1–21. [Google Scholar] [CrossRef]
  46. Hellriegel, B. Modelling the immune response to malaria with ecological concepts: Short-term behaviour against long-term equilibrium. Proc. R. Soc. Lond. B 1992, 250, 249–256. [Google Scholar]
  47. Chiyaka, C.; Garira, W.; Dube, S. Modelling immune response and drug therapy in human malaria infection. Comput. Math. Methods Med. 2008, 9, 143–163. [Google Scholar] [CrossRef]
  48. Holmes, M.H. Introduction to Perturbation Methods; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012; Volume 20. [Google Scholar]
  49. Garira, W.; Mathebula, D.; Netshikweta, R. A mathematical modelling framework for linked within-host and between-host dynamics for infections with free-living pathogens in the environment. Math. Biosci. 2014, 256, 58–78. [Google Scholar] [CrossRef]
Figure 1. A conceptual representation of the universal theory for multiscale modelling of infectious disease dynamics. In this case, an infectious disease system is considered to consist of two different disease system forms—a pathogen-centred disease system form and a host-centred disease system form—in such a way that the overall multiscale model of infectious disease dynamics M consists of two complementary multiscale models: M p , a process-based multiscale model that describes the scales of a level of organisation of the pathogen-centred disease system form, and M m , a mechanism-based multiscale model that describes the scales of the corresponding levels of organisation of the host-centred disease system form so that M = f ( M p , M m ) .
Figure 1. A conceptual representation of the universal theory for multiscale modelling of infectious disease dynamics. In this case, an infectious disease system is considered to consist of two different disease system forms—a pathogen-centred disease system form and a host-centred disease system form—in such a way that the overall multiscale model of infectious disease dynamics M consists of two complementary multiscale models: M p , a process-based multiscale model that describes the scales of a level of organisation of the pathogen-centred disease system form, and M m , a mechanism-based multiscale model that describes the scales of the corresponding levels of organisation of the host-centred disease system form so that M = f ( M p , M m ) .
Mathematics 11 03874 g001
Figure 2. A schematic representation of the seven main levels of organisation of the pathogen-centred disease system form based on the universal theory for multiscale dynamics of infectious disease dynamics, which are [a.] the cell level, [b.] the tissue level, [c.] the organ/microcommunity level, [d.] the microecosystem level, [e.] the whole organism level, [f.] the macrocommunity level, and [g.] the macroecosystem level, in a curved and discretised four-dimensional space-time that combines the three dimensions of space and one dimension of time into a single four-dimensional manifold. The seven levels of organisation of the pathogen-centred disease system form are hierarchically organised in space-time.
Figure 2. A schematic representation of the seven main levels of organisation of the pathogen-centred disease system form based on the universal theory for multiscale dynamics of infectious disease dynamics, which are [a.] the cell level, [b.] the tissue level, [c.] the organ/microcommunity level, [d.] the microecosystem level, [e.] the whole organism level, [f.] the macrocommunity level, and [g.] the macroecosystem level, in a curved and discretised four-dimensional space-time that combines the three dimensions of space and one dimension of time into a single four-dimensional manifold. The seven levels of organisation of the pathogen-centred disease system form are hierarchically organised in space-time.
Mathematics 11 03874 g002
Figure 3. A conceptual representation of a typical primary level of organisation of the pathogen-centred disease system form where the multiscale dynamics of this disease system form involve interaction between a microscale and a macroscale through four main disease processes, which are [a.] the infection or super-infection by pathogen process, [b.] the pathogen replication process, [c.] the pathogen shedding or excretion process, and [d.] the pathogen transmission process, in a curved four-dimensional space-time that combines the three dimensions of space and one dimension of time into a single four-dimensional manifold. The dotted line indicates a conceptual representation of the path of the pathogen along the multiscale cycle.
Figure 3. A conceptual representation of a typical primary level of organisation of the pathogen-centred disease system form where the multiscale dynamics of this disease system form involve interaction between a microscale and a macroscale through four main disease processes, which are [a.] the infection or super-infection by pathogen process, [b.] the pathogen replication process, [c.] the pathogen shedding or excretion process, and [d.] the pathogen transmission process, in a curved four-dimensional space-time that combines the three dimensions of space and one dimension of time into a single four-dimensional manifold. The dotted line indicates a conceptual representation of the path of the pathogen along the multiscale cycle.
Mathematics 11 03874 g003
Figure 4. A conceptual representation of the four main levels of organisation of the host-centred disease system form, which are: [a.] the host damage level, [b.] the host immune response level, [c.] the environmental change level, and [d.] the biodiversity change level. These levels of organisation of the host-centred disease system form are hierarchically organised into physical size scale and time scale in such a way that they can be placed at distinct and discrete positions in the hierarchy so that the larger the physical size scale of the biological functional structure that is central to a particular compositional physiological mechanism, the larger the time scale of the time at which the compositional physiological mechanism takes place in the hierarchy.
Figure 4. A conceptual representation of the four main levels of organisation of the host-centred disease system form, which are: [a.] the host damage level, [b.] the host immune response level, [c.] the environmental change level, and [d.] the biodiversity change level. These levels of organisation of the host-centred disease system form are hierarchically organised into physical size scale and time scale in such a way that they can be placed at distinct and discrete positions in the hierarchy so that the larger the physical size scale of the biological functional structure that is central to a particular compositional physiological mechanism, the larger the time scale of the time at which the compositional physiological mechanism takes place in the hierarchy.
Mathematics 11 03874 g004
Figure 5. A conceptual representation of the seven main scales that are used to describe each of the four main functional levels of organisation (i.e., the host damage level, the host immune response level, the environmental change level, and the biodiversity change level) of the host-centred disease system form, which are [a.] the molecular scale, [b.] the cell scale, [c.] the tissue scale, [d.] the organ scale, [e.] the whole organism scale, [f.] the community scale, and [g.] the ecosystem scale, using mechanism-based multiscale models. These scales of a level of organisation of the host-centred disease system form are hierarchically organised into physical size scale and time scale so that they can be placed at distinct and discrete positions in the hierarchy such that the larger the physical size scale of the biological functional structure that is central to a particular physiological mechanism, the larger the time scale at which the physiological mechanism takes place.
Figure 5. A conceptual representation of the seven main scales that are used to describe each of the four main functional levels of organisation (i.e., the host damage level, the host immune response level, the environmental change level, and the biodiversity change level) of the host-centred disease system form, which are [a.] the molecular scale, [b.] the cell scale, [c.] the tissue scale, [d.] the organ scale, [e.] the whole organism scale, [f.] the community scale, and [g.] the ecosystem scale, using mechanism-based multiscale models. These scales of a level of organisation of the host-centred disease system form are hierarchically organised into physical size scale and time scale so that they can be placed at distinct and discrete positions in the hierarchy such that the larger the physical size scale of the biological functional structure that is central to a particular physiological mechanism, the larger the time scale at which the physiological mechanism takes place.
Mathematics 11 03874 g005
Figure 6. A schematic representation of the process–based multiscale model of the pathogen-centred disease system form of malaria disease system (1). This coupled multiscale model incorporates super-infection in both the human host and the mosquito host in addition to the other three disease processes for both the human host level and the mosquito host level. The four disease processes for a pathogen-centred disease system form for each whole organism level are illustrated in Figure 3.
Figure 6. A schematic representation of the process–based multiscale model of the pathogen-centred disease system form of malaria disease system (1). This coupled multiscale model incorporates super-infection in both the human host and the mosquito host in addition to the other three disease processes for both the human host level and the mosquito host level. The four disease processes for a pathogen-centred disease system form for each whole organism level are illustrated in Figure 3.
Mathematics 11 03874 g006
Figure 7. Shows changes in (a) community gametocyte load ( G H ), (b) population of infected mosquitoes ( I V ), (c) population of infected humans ( I H ), and (d) community sporozoite load ( P V ) for different values of the rate at which the gametocytes develop and become infectious α h : α h = 0.02 , α h = 0.04 , α h = 0.08 .
Figure 7. Shows changes in (a) community gametocyte load ( G H ), (b) population of infected mosquitoes ( I V ), (c) population of infected humans ( I H ), and (d) community sporozoite load ( P V ) for different values of the rate at which the gametocytes develop and become infectious α h : α h = 0.02 , α h = 0.04 , α h = 0.08 .
Mathematics 11 03874 g007
Figure 8. Shows the changes in (a) community gametocyte load ( G H ), (b) population of infected mosquitoes ( I V ), (c) population of infected humans ( I H ), and (d) community sporozoite load ( P V ) for different values of the rate at which sporozoites become infectious to humans α v : α v = 0.025 , α v = 0.045 , α v = 0.085 .
Figure 8. Shows the changes in (a) community gametocyte load ( G H ), (b) population of infected mosquitoes ( I V ), (c) population of infected humans ( I H ), and (d) community sporozoite load ( P V ) for different values of the rate at which sporozoites become infectious to humans α v : α v = 0.025 , α v = 0.045 , α v = 0.085 .
Mathematics 11 03874 g008
Figure 9. Shows changes in (a) gametocyte–infected erythrocytes ( G h ), (b) population of merozoites ( M h ), (c) merozoite–infected erythrocytes ( R m ), and (d) susceptible erythrocytes ( R h ) for different values of contact rate of mosquitoes with the infectious reservoir of humans β H : β H = 0.00356 , β H = 0.0356 , β H = 0.356 .
Figure 9. Shows changes in (a) gametocyte–infected erythrocytes ( G h ), (b) population of merozoites ( M h ), (c) merozoite–infected erythrocytes ( R m ), and (d) susceptible erythrocytes ( R h ) for different values of contact rate of mosquitoes with the infectious reservoir of humans β H : β H = 0.00356 , β H = 0.0356 , β H = 0.356 .
Mathematics 11 03874 g009
Figure 10. Shows changes in (a) gametocyte-infected erythtcytes ( G h ), (b) population of merozoites ( M h ), (c) merozoite-infected erythrocytes ( R m ), and (d) susceptible erythrocytes ( R h ) for different values of contact rate of humans with the infectious reservoir of mosquitoes β V : β V = 0.003 , β V = 0.03 , β V = 0.3 .
Figure 10. Shows changes in (a) gametocyte-infected erythtcytes ( G h ), (b) population of merozoites ( M h ), (c) merozoite-infected erythrocytes ( R m ), and (d) susceptible erythrocytes ( R h ) for different values of contact rate of humans with the infectious reservoir of mosquitoes β V : β V = 0.003 , β V = 0.03 , β V = 0.3 .
Mathematics 11 03874 g010
Figure 11. Shows changes in (a) sporozytes ( P v ), (b) oocysts ( O v ), (c) zygotes ( Z v ), and (d) gametocyte-infected erythrocytes ( G v ) for different values of contact rate of mosquitoes with the infectious reservoir of humans β H : β H = 0.00356 , β H = 0.0356 , β H = 0.356 .
Figure 11. Shows changes in (a) sporozytes ( P v ), (b) oocysts ( O v ), (c) zygotes ( Z v ), and (d) gametocyte-infected erythrocytes ( G v ) for different values of contact rate of mosquitoes with the infectious reservoir of humans β H : β H = 0.00356 , β H = 0.0356 , β H = 0.356 .
Mathematics 11 03874 g011
Figure 12. Shows changes in (a) sporozytes ( P v ), (b) oocysts ( O v ), (c) zygotes ( Z v ), and (d) gametocyte-infected erythrocytes ( G v ) for different values of contact rate of humans with the infectious reservoir of mosquitoes β V : β V = 0.003 , β V = 0.03 , β V = 0.3 .
Figure 12. Shows changes in (a) sporozytes ( P v ), (b) oocysts ( O v ), (c) zygotes ( Z v ), and (d) gametocyte-infected erythrocytes ( G v ) for different values of contact rate of humans with the infectious reservoir of mosquitoes β V : β V = 0.003 , β V = 0.03 , β V = 0.3 .
Mathematics 11 03874 g012
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Garira, W.; Muzhinji, K. The Universal Theory for Multiscale Modelling of Infectious Disease Dynamics. Mathematics 2023, 11, 3874. https://doi.org/10.3390/math11183874

AMA Style

Garira W, Muzhinji K. The Universal Theory for Multiscale Modelling of Infectious Disease Dynamics. Mathematics. 2023; 11(18):3874. https://doi.org/10.3390/math11183874

Chicago/Turabian Style

Garira, Winston, and Kizito Muzhinji. 2023. "The Universal Theory for Multiscale Modelling of Infectious Disease Dynamics" Mathematics 11, no. 18: 3874. https://doi.org/10.3390/math11183874

APA Style

Garira, W., & Muzhinji, K. (2023). The Universal Theory for Multiscale Modelling of Infectious Disease Dynamics. Mathematics, 11(18), 3874. https://doi.org/10.3390/math11183874

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop