Loading [MathJax]/jax/output/SVG/jax.js
Research article Special Issues

The effect of the surface properties of poly(methyl methacrylate) on the attachment, adhesion and retention of fungal conidia

  • Received: 11 May 2020 Accepted: 12 June 2020 Published: 23 June 2020
  • Poly(methyl methacrylate) (PMMA) surfaces, (commercial PMMA (PMMAc), spin coated PMMA (PMMAsc) and a 90% methylmethacrylate/10% 3-methacryloxypropyltrimethoxysilane random copolymer (P(MMA-co-gMPS)) were used to determine the effect of surface properties on conidia biofouling. The contact angles of the substrates demonstrated that the PMMAsc and the P(MMA-co-gMPS) polymer (62.8°) were more wettable than the PMMAc surface (71.0°). The PMMAsc had the greatest roughness value (32.0 nm) followed by the PMMAc (3.0 nm), then P(MMA-co-gMPS) (1 nm). Aspergillus niger 1957 conidia were spherical, smooth and hydrophobic (12.1%). Aspergillus niger 1988 conidia were spherical with spikes and hydrophobic (17.1%). Aureobasidium pullulans was elliptical with longitudinal ridges and hydrophilic (79.9%). Following attachment assays, cPMMA attached the greatest numbers of conidia. Following the adhesion and retention assays (washing step included in the protocol), A. niger 1957 and A. niger 1988 were least adhered to the P(MMA-co-gMPS) surface, whilst A. pulluans was least adhered to the PMMAsc surface. This work demonstrated that in the absence of a washing step, only the surface properties influenced the conidia attachment, whilst in the presence of a washing step, both the properties of the surfaces and the conidia affected conidia adhesion and retention. Hence, the methodology used (with or without a washing step) should reflect the environment in which the surface is to be applied.

    Citation: Kathryn A. Whitehead, Christopher M. Liauw, Joels S. T. Wilson-Nieuwenhuis, Anthony J. Slate, Ted Deisenroth, Andrea Preuss, Joanna Verran. The effect of the surface properties of poly(methyl methacrylate) on the attachment, adhesion and retention of fungal conidia[J]. AIMS Bioengineering, 2020, 7(3): 165-178. doi: 10.3934/bioeng.2020015

    Related Papers:

    [1] Hsin-Lun Li . Leader–follower dynamics: stability and consensus in a socially structured population. AIMS Mathematics, 2025, 10(2): 3652-3671. doi: 10.3934/math.2025169
    [2] Xingyun Shi, Xuqiong Luo . Leader-follower consensus of a fractional-order multi-agent system based on event-triggered control and extended state observer. AIMS Mathematics, 2025, 10(7): 15785-15810. doi: 10.3934/math.2025707
    [3] Hongjie Li . Event-triggered bipartite consensus of multi-agent systems in signed networks. AIMS Mathematics, 2022, 7(4): 5499-5526. doi: 10.3934/math.2022305
    [4] Qiangqiang Zhang, Yiyan Han, Chuandong Li, Le You . Constraint impulsive consensus of nonlinear multi-agent systems with impulsive time windows. AIMS Mathematics, 2020, 5(4): 3682-3701. doi: 10.3934/math.2020238
    [5] Kairui Chen, Yongping Du, Shuyan Xia . Adaptive state observer event-triggered consensus control for multi-agent systems with actuator failures. AIMS Mathematics, 2024, 9(9): 25752-25775. doi: 10.3934/math.20241258
    [6] Jiaqi Liang, Zhanheng Chen, Zhiyong Yu, Haijun Jiang . Fixed-time consensus of second-order multi-agent systems based on event-triggered mechanism under DoS attacks. AIMS Mathematics, 2025, 10(1): 1501-1528. doi: 10.3934/math.2025070
    [7] Fang Zhu, Pengtong Li . Adaptive fuzzy consensus tracking control of multi-agent systems with predefined time. AIMS Mathematics, 2025, 10(3): 5307-5331. doi: 10.3934/math.2025245
    [8] Lu Zhi, Jinxia Wu . Adaptive constraint control for nonlinear multi-agent systems with undirected graphs. AIMS Mathematics, 2021, 6(11): 12051-12064. doi: 10.3934/math.2021698
    [9] Jiawen Li, Yi Zhang, Heung-wing Joseph Lee, Yingying Wang . Fuzzy tracking control of singular multi-agent systems under switching topology. AIMS Mathematics, 2024, 9(11): 29718-29735. doi: 10.3934/math.20241440
    [10] Yuming Chen, Jie Gao, Luan Teng . Distributed adaptive event-triggered control for general linear singular multi-agent systems. AIMS Mathematics, 2023, 8(7): 15536-15552. doi: 10.3934/math.2023792
  • Poly(methyl methacrylate) (PMMA) surfaces, (commercial PMMA (PMMAc), spin coated PMMA (PMMAsc) and a 90% methylmethacrylate/10% 3-methacryloxypropyltrimethoxysilane random copolymer (P(MMA-co-gMPS)) were used to determine the effect of surface properties on conidia biofouling. The contact angles of the substrates demonstrated that the PMMAsc and the P(MMA-co-gMPS) polymer (62.8°) were more wettable than the PMMAc surface (71.0°). The PMMAsc had the greatest roughness value (32.0 nm) followed by the PMMAc (3.0 nm), then P(MMA-co-gMPS) (1 nm). Aspergillus niger 1957 conidia were spherical, smooth and hydrophobic (12.1%). Aspergillus niger 1988 conidia were spherical with spikes and hydrophobic (17.1%). Aureobasidium pullulans was elliptical with longitudinal ridges and hydrophilic (79.9%). Following attachment assays, cPMMA attached the greatest numbers of conidia. Following the adhesion and retention assays (washing step included in the protocol), A. niger 1957 and A. niger 1988 were least adhered to the P(MMA-co-gMPS) surface, whilst A. pulluans was least adhered to the PMMAsc surface. This work demonstrated that in the absence of a washing step, only the surface properties influenced the conidia attachment, whilst in the presence of a washing step, both the properties of the surfaces and the conidia affected conidia adhesion and retention. Hence, the methodology used (with or without a washing step) should reflect the environment in which the surface is to be applied.



    1. Introduction

    Contact tracing is a mitigation strategy that aims at immediately detecting, testing, and treating the next-generation cases during the spreading of an infectious disease. Such local targeted control measure is very effective when the number of cases is limited, for example at the early stage of an outbreak. In 2014, West Africa experienced the most widespread Ebola epidemic in the history with more than 28,000 reported cases. Secondary infections were reported in several European countries and the United States. The Ebola virus is transmitted via physical contact with the infected individuals or their body fluids; the infected ones can transmit the virus to their contacts after becoming symptomatic [4]. In the case of the Ebola epidemic, the objective of contact tracing is to identify and monitor the individuals who have been exposed to the infectious ones for 21 days [27]. This procedure allows for a prompt isolation of the contacts of an infectious individual as soon as he/she becomes symptomatic. Contact tracing has shown effectiveness in several cases. In 2014, there was an Ebola virus disease (EVD) outbreak in Nigeria due to a traveler who returned from Liberia. An extensive contact tracing effort took place starting from day 3 and a total of 894 contacts were traced, all linked to the single index case [21]. Compared to other regions of West Africa, the outbreak of Nigeria was better contained with only 19 confirmed and 1 probable EVD cases out of which 8 died. The improved outcomes can be attributed to the early detection of the index case and effective isolation of infectious individuals due to contact tracing [6]. In September 2014, a person in Dallas - TX, returning from Liberia, was diagnosed with EVD. All of his contacts were traced and monitored for 21 days. Two healthcare workers who provided care for the index case were also diagnosed positive and one of them took a round trip to Cleveland, OH before detection. This prompted the Centers for Disease Control and Prevention (CDC) to trace all the passengers of the two flights. Another case was detected when a person returned to New York from Guinea in October 2014 and contacts of that individual were traced as well. A total of 458 contacts were traced in Texas, Ohio, and New York [26]. Contact tracing was an appropriate approach to stop the transmission of Ebola in the USA as the number of cases reported was quite small [5]. In this paper, we attempt to quantify the effect of contact tracing such early stages of the epidemic.

    In general, contact tracing can be carried out using different protocols, depending on the characteristics of the pathogen transmission. Definition of contact, duration and frequency of monitoring, are some examples of variables to consider while implementing contact tracing. Additionally, delays in implementing the contact tracing process are possible in realistic scenarios. Predicting the effectiveness of contact tracing as a function of the disease characteristics, disease stage, and protocol characteristics is a challenging task. Eames and Keeling have proposed a formula to correlate the effectiveness of contact tracing and the basic reproductive ratio by using detailed pairwise equations for a susceptible-infected-removed (SIR) model [5]. The basic reproductive ratio, R0, is a key indicator in epidemiology and represents the expected number of secondary infections over all possible initial infections during their infectious period [8][3]. This ratio is a crucial tool for a quantitative measurement of the severity of a disease outbreak and helps the public health authorities to evaluate the risk of an outbreak in the emergence of an infectious disease.

    Klinkenberg et al., evaluated the impact of time-related characteristics of the infection and the tracing process on the success of contact tracing [11]. They showed why contact tracing is effective for control of smallpox and SARS, only partially effective for foot-and-mouth disease and likely not effective for influenza [11]. The impact of contact pattern on the efficacy of contact tracing has also been studied. Researchers have concluded that not only the disease properties but also the contact network properties are crucial for the success of contact tracing. For example, Eames and Keeling indicated that contact tracing effectiveness increases with network clustering [5]. Kiss et al., showed that contact tracing is typically ineffective for random contact networks with high average node degrees and small clustering coefficients [10].

    Researchers have also attempted to analyze the impact of intervention strategies on the recent Ebola epidemic using mean-field compartmental models, which can be either stochastic or deterministic in nature. Browne et al., used a deterministic version of the compartment models and separated the infected individuals into different compartments based on whether they are hospitalized or unreported. They evaluated the impact of relevant epidemiological properties of Ebola on contact tracing efficiency and presented a formula to determine the minimum number of contacts to be traced per identified infectious individual in order to bring down the effective reproductive ratio below one [1]. In [16], Rizzo et al., adopted a susceptible-exposed-infected-removed (SEIR) compartmental model with additional compartments for hospitalized and dead people who had traditional funerals. Then, they adopted activity driven networks (ADN), where each individual has a network of contacts which depends on an activity potential and vary with time. ADN describes contact processes that evolve over time-varying networks [14]. The analysis of Rizzo et al., showed that contact tracing and other intervention policies adopted later in the West Africa would have drastically mitigated the epidemic spreading if used promptly.

    Despite being a relatively new area of research, there have been several works involving ADNs. Starnini et al., used ADNs to study temporal percolation properties and showed how SIR models can be mapped to the percolation problem [22]. In another work, they explored the relation between network topology and activity potential distribution to obtain analytical expressions for several topological properties of the integrated social networks [23]. Mata et al., studied a power law distribution of activity potentials and found very slow relaxation dynamics and aging in random walks [13]. Sun et al., investigated ADNs with Markovian and non-Markovian dynamics and found that memory slows down the spreading process in SIR model and boosts the spreading process in SIS models [24]. Perra et al., experimented random walks in time-varying networks and found that results vary significantly. They concluded that the network dynamics should be considered to avoid misleading results in practical situations [15].

    For a realistic study of contact tracing effectiveness in the early-stage of an Ebola outbreak, stochastic and individual-level models are needed. When considering a small number of cases, the localized and highly-structured contacts of infected individuals prominently influence the numbers, the timings, and the locations of the future cases. In this scenario, the accuracy of meta-population models, characterized by high levels of aggregation, dramatically deteriorates. The challenge is that successful modeling approaches to evaluate the effectiveness of contact tracing need to take into account the highly structured network of contacts and data on the network of contacts is often not available or too large. The real world networks are not static. The set of people with whom a person remain in contact changes with time. Temporal networks can incorporate these changes and accurately represent real behaviors in human populations. Agent-based models are accurate but computationally expensive. However, ADNs provide a tractable way to produce accurate results [18]. Our work focuses on the microscopic processes that occur at the beginning of an Ebola outbreak. We use ADNs to capture the contact dynamics using a constant underlying stochastic process. We also estimate the basic reproductive number of the disease spread in different scenarios.

    In this paper, we evaluate the effectiveness of contact tracing using a novel modeling framework. Our modeling framework has several features to characterize early-stage Ebola transmission: 1) the network model is patient-centric because when the number of infected cases is small, only the myopic networks of infected individuals matter and rest of the possible social contacts are irrelevant, 2) the Ebola disease model is individual-based and stochastic because during the early stages of a spreading process, the random fluctuations are significant and should be captured appropriately, 3) the contact tracing model is parameterizable to incorporate different critical aspects of the contact tracing protocols.

    Our model is built on susceptible-exposed-infected-hospitalized-removed compartments (SEIHR), where susceptible, exposed and infected individuals may become monitored-susceptible, monitored-exposed, and monitored-infected respectively as a consequence of contact tracing. We propose an activity driven network approach to consider the inherent time-varying nature of the contagion process in a host population, in the form of variations in connectivity pattern of contacts. Activity driven networks are characterized by an activity firing rate assigned to each node at each time step, so that a fraction of the nodes create new links with others and the contagion process develops over those new paths. This is a general model for evaluating contact tracing, intertwined within the transmission model, which can be applied to a wide variety of diseases beyond the specific case of Ebola. Using this modeling framework, we perform extensive simulations varying important delays: identification delay (local), starting delay (global), and hospitalization delay. Through simulations, we quantify the effectiveness of contact tracing and compute R0.

    The remainder of this paper is organized as follows. In section 2, we propose a compartmental model for Ebola transmission incorporating contact tracing, discuss an overview of activity driven network (ADN) and explain it based on the proposed compartmental model. Section 3 presents a method to compute the basic reproductive number in a heterogeneous network. In section 4, we define true positive and false positive ratios based on the proposed model to plot the receiver operating characteristic (ROC) curve. Section 5 summarizes the main results of this article and section 6 contains the concluding remarks.


    2. Mathematical modeling of Ebola disease spreading incorporating contact tracing

    Spreading of an infectious disease is a complex event with many interacting variables. One of the primary tools to analyze and predict the disease diffusion as well as the severity of infectious disease is the compartmental model. Compartmental models are the mathematical frameworks that can capture some major features of epidemic spreading such as pathogen transmission probabilities and host transition rates from one state to another [2]. In this work, we employ a discrete-time expression of the susceptible, exposed, infected, hospitalized and removed (SEIHR) compartmental model. Such model is compatible with the epidemiology of Ebola.


    2.1. Compartmental model

    SEIHR model is an extension of susceptible, exposed, infected and removed (SEIR) model [9] with an additional compartment H, where H stands for hospitalized. Each compartment variable denotes a fraction of individuals who belong to one of the following states: susceptible (S), exposed (E), infected (I), hospitalized (H) or removed (R). An individual can undergo transitions from one state to another during the disease evolution. To evaluate the impact of contact tracing for the detection of new Ebola patients, we add three more compartments: traced-susceptible (ST), traced-exposed (ET) and traced-infected (IT) to the (SEIHR) model. We follow the contact tracing implementation guideline published by the World Health Organization (WHO) and the Centers for Disease Control and Prevention (CDC). In the guideline it is mentioned that when an infectious individual enters a hospital and his/her laboratory results come out positive, any person who had contact with him/her in the last 21 days should be traced [27]. Therefore, in our proposed model, whenever an infectious individual enters the hospitalized (H) state, all the individuals exposed to the infected one will be identified and followed up for 21 days.

    We classify transitions between different epidemiological compartments of the proposed model into two groups: node-based transition and edge-based transition. In a node-based transition, a node moves from one state to another individually and the transition does not depend on the states of the node's neighbors. Contrary to the nodal transition, the edge-based transition is dependent on the states of a node's neighbors. Based on these definitions, we describe the transition processes of our proposed model as follows:

    ● Edge-based transition: When a susceptible individual has contact with an infectious or hospitalized one, he/she moves to exposed state with probability β. The probability is multiplied if there are multiple infected or hospitalized individuals in contact.

    ● Node-based transition: An exposed individual undergoes an average incubation period of 1/λ before proceeding to the infected (I) state [19]. An infected individual moves to the hospitalized (H) state with average delay of probability 1/γ and their susceptible, exposed and infectious contacts move to ST, ET and IT compartments respectively, with an average identification delay period of 1/α. Since a portion of the the contacts of the hospitalized individuals might be inaccessible or a portion of them might be unwilling to report all their contacts immediately, we assume that there is a delay to identify those contacts. An exposed individual who is traced, (referred to as the ET compartment) undergoes an average incubation period 1/λ days before progressing to the infectious compartment where infections are traced (IT). An infectious individual who is traced enters the hospitalized compartment with an average delay of 1/γT where γT>γ. Finally, a hospitalized individual moves to the removed compartment with an average delay of 1/δ.

    A schematic of the epidemiological transition processes of the proposed model is depicted in figure 1.

    Figure 1. Schematic of the transition processes in the Ebola progression with contact tracing model.

    In our model, we ignore demography and since the population is much greater than the number of infected people, we assume that S1. In addition to that, we assume that no individual will die or recover without hospitalization, despite the fact that Ebola has shown high mortality rates in the West African outbreaks. This is a realistic scenario for small-scale outbreaks in the United States. We implement stochastic transitions in our simulation program, based on the above-mentioned state transition rules.

    Based on the proposed model, we develop a quantitative approach to measure the effectiveness of contact tracing implementations. To asses the impact of contact tracing protocols in Ebola disease spreading before the epidemic phase, we propose two measures: missed-detection probability and contact tracing cost.

    Definition 2.1. To assess the risk detection capabilities of contact tracing efforts for Ebola, we introduce the missed-detection probability. Missed-detection probability denotes the probability that a secondary infected individual is not detected before transmitting the virus to others. Based on our model, we propose the missed-detection probability as follows: NEINEET+NEI, where, N represents the total number of individuals who move from one compartment to another indicated by the arrow.

    Definition 2.2. The aim of contact tracing is to facilitate the detection of secondary infections from the contacts of an infected person. However, a large proportion of an infectious individual's contacts could remain uninfected (susceptible). We define contact tracing cost as the number of detected individuals who had contact with infections but were not infected. Based on the proposed model, the definition of contact tracing cost is: Cost=NSST, where, NSST represents the total number of susceptible individuals who move to the ST compartment.


    2.2. Activity driven network

    Disease contagious process and network structure are two important elements which can have significant impacts on disease spreading [5]. Many intervention strategies such as contact tracing, target strategy and egocentric strategy aim at controlling the contagious process based on interactions between individuals in a social network [12] [1]. In particular, the contact tracing strategy or the identification of individuals who have contact with infections is fundamentally linked to potential transmission paths in the network [10] [5]. The goal of contact tracing is to identify all the potential routes in the network and isolate all the new infected individuals, before they become infectious [10]. Here, we implement activity driven network (ADN) to capture interactions between nodes in a network over a specific period of time and assess effectiveness of contact tracing strategy for a temporal network based on Ebola contagious process. Activity driven network is a random and memoryless process which can capture structural features of a network such as the evolution of contact patterns over time [14].

    Activity driven network considers an activity firing rate ai for each node which is the probability of establishing links with other nodes per unit of time [12]. Activity firing rates are assigned according to a probability distribution F(a), which can describe network dynamics and the corresponding structure [12]. Typically F(a) is a heavy tail density function: F(a)ac, where, 2c3 and a[ϵ,1] with ϵ=103 [14]. At each time increment Δt, an active node generates m links with m other nodes that are selected randomly. The generative network process in an increment time Δt is listed as follows [12],

    ● At time t, the network Gp(N,m) has N disconnected vertices.

    ● Each node i, with probability pi=aiηiΔt, becomes active and generates m links with m other nodes. Here, ηi is a constant scaling factor for node i.

    ● At time t+Δt, all the edges in network Gp are removed and the process is repeated. The activity potentials (ai) of nodes remain the same, so the underlying generator does not change.


    2.3. ADN for Ebola contagion process

    We implement the activity driven network (ADN) to generate a random network at each time step Δt. Then, we simulate Ebola contagion process using the proposed compartmental model. In our network generation process, we take discrete time steps with Δt=1, and at each step only the nodes which are in susceptible, infected, exposed or hospitalized states can generate new links with probabilities denoted by pi's. Here, piP, where P is the set of activation probabilities of the nodes that can become active. The contact tracing mechanism is constructed based on "CDC emergency guidelines of implementation and management of contact tracing" for Ebola virus disease [27]. Based on the guideline, any person who has a potential exposure to an Ebola Virus Disease (EVD) case, should go under observation for 21 days [27]. To implement contact tracing strategy in a temporal network such as ADN, we capture all the nodes that become neighbors of an infectious node j from the beginning. We also keep track of the time since the latest contact. We implement contact tracing after a delay of TCT since the first case is identified. It runs until all the infected people have recovered and the contacts have been monitored for 21 days. Algorithm 1 sets the rules to produce and simulate Ebola virus spreading in a host population for discrete time 1tT, where T is the end time of the simulation. All the steps and transitions shown in algorithm 1 are done in a stochastic manner to simulate realistic disease dynamics.

    Algorithm 1: ADN for Ebola Contagion Process
    1: Set NA=ϕ, NI=ϕ, INeigh=ϕ and tinitial=0
    2: while tT do
    3:   if i|xi{S,E,I,H} then
    4:      NAi
    5:   end if
         where xi represents the state of node i. And NA is the set of active nodes.
    6:   1. Network generation:
        Generate a network Gp(N,m), w.r.t (P,NA)
    7:   2. Contact identification:
    8:   if i|xi{H} then
    9:     tinitial=tinitial+1
    10:     NIi
    11:   end if
          where NI is the set of identified nodes. tinitial keeps track of the time since the first identified case.
    12:   3. Contact tracing:
    13:   if tinitialTCT then
    14:     for all iNI do
    15:       Update the states of node i's current contacts and untraced past contacts in (INeigh) since it became infectious to their respective new states with probability α (SST, EET, IIT).
    16:       end for
    17:   end if
    18:   4. Edge-based transition: Find susceptible nodes in contact with infectious nodes and update their state based on edge-based transition rule. Add all nodes that are neighbors of infected nodes to the set INeigh.
    19:   5. Node-based transition: Update the states of nodes other than susceptible; i|xiS based on node-based transition rules.
    20:   6. Tracing removal: Remove nodes belonging to INeigh which were traced for 21 days but not detected. Return these nodes from ST and ET compartments to S and E compartments respectively.
    21:   t=t+1
    22: end while

    3. Reproductive number in heterogeneous network

    The basic reproductive ratio, R0, is a descriptor of epidemic potential in the mathematical modeling of infectious diseases. This quantity helps the public health authorities to assess the risk of an outbreak in the emergence of an infectious disease [19]. Furthermore, early estimation of the basic reproductive number helps the healthcare authorities to plan appropriate control measures. A general definition of basic reproductive number is the expected number of secondary infections over all possible initial infections during their infectious period [8] [3]. Based on the general definition of reproductive number and characteristic of the heterogeneous network, we propose a suitable definition for R0 in a heterogeneous network. In heterogeneous network, contact patterns tend to have a high variability in prevalence and so, besides high degree nodes, there are some low degree nodes which may have no contacts with others [20]. Based on this characteristic, to compute R0 in a heterogeneous network, we only consider those initial infected nodes which establish links with other nodes and transmit infections. Therefore, we define R0 as follows:

    Definition 3.1. The basic reproductive number (R0) is the expected number of secondary infectious cases over all initial infections that establish interaction with others during the infectious period. In our activity driven network which is a heterogeneous network, we compute R0 as:

    R0=Total number of hosts infected by I0Active subset of I0 (1)

    where, I0 is the set of initially infected nodes at time t=0. We take care of the fact that every infected node may not be able to transmit infection due to network heterogeneity. The infected nodes that do not have any links to others will not be able to transmit infection. We do not include them in our calculation and only count the active subset (infected nodes with links to others) of the initially infected hosts.


    4. Receiver operating characteristic

    A receiver operating characteristic or ROC curve is a fundamental method to illustrate the performance of a system such as separating true positive results from false positive results in a test or comparing two alternative tasks [7]. In a ROC curve, we plot Sensitivity or true positive ratio (TPR) as a function of (1Specificity) or false positive ratio (FPR). TPR is defined as the fraction of samples that are detected correctly and FPR is defined as the fraction of samples that are identified as positive, incorrectly.

    The major aim of contact tracing is to identify exposed individuals before they become infectious, in order to halt the chain of pathogen transmission. We define the positive samples as those secondary infectious individuals who had contact with infections during the disease evolution. The true positives in the context of contact tracing are defined as the secondary infectious individuals who were traced. Similarly, false positives are defined as those individuals who were not infected but traced as possible secondary infections. Therefore, based on the different epidemiological compartments in our proposed model, we define Sensitivity and 1Specificity as the following:

    Definition 4.1. We define the true positive ratio (TPR) or Sensitivity as,

    TPR=Number of exposed hosts tracedTotal exposed hosts (2)

    Definition 4.2. We define the false positive ratio or 1Specificity as,

    FPR=Number of usceptible hosts tracedTotal susceptible hosts with infectious neighbors (3)

    A point (p,q) in the ROC shows that with probability p, the susceptible individuals who had contact with infections could be identified as healthy and with probability q, the infected individuals who had contacts with infections could be identified as infectious.


    5. Results

    To generate realizations for Ebola disease spreading without any immunization strategy, parameters of our proposed model for contagious process are given in table 1. The parameter values used in our proposed model were inspired from the works of Rizzo et al.[16]. The time unit is day, and all the rate/probability parameters are given as day1 values. We assume that the number of initially infected individuals is, I0=2 and each active node can generate m=7 links with other nodes where the total number of nodes is, N=1000. In a usual case, sickness from infection causes a reduction in an individuals social interaction. Therefore, we assume the scaling factors of the activity firing rates for the hospitalized, infected, and susceptible individuals as following: ηH<<ηI<ηS. Activity driven network's parameters are shown in table 2.

    Table 1. Time-invariant parameters of Ebola contagion process.
    Parameter Value
    Transmission probability(β) 0.11
    Incubation rate (λ) 0.095
    Recovery/removal probability (δ) 0.1
    Hospitalization probability in existence of contact tracing (γT) 0.9
    Hospitalization probability (γ) 0.33
     | Show Table
    DownLoad: CSV
    Table 2. Parameters of activity driven network generator.
    Parameter Value
    Density function exponent (c) 2.2
    Links per active node (m) 7
    Scaling factor for susceptible (ηS) 2.2
    Scaling factor for infected (ηI) 1.1
    Scaling factor for hospitalized (ηH) 0.005
     | Show Table
    DownLoad: CSV

    To assess the effectiveness of contact tracing in the early stage of the epidemic, we assume three different implementation-time scenarios for contact tracing. The first one is when we implement contact tracing from the beginning (TCT=1), the second scenario is when contact tracing is started on day 9 (TCT=9), and the third one is when contact tracing is implemented on day 22 (TCT=22). To evaluate the effectiveness of contact tracing in a more realistic situation, we consider five identification delay times, α1{1,2,5,10,20} days for the three above mentioned scenarios of contact tracing implementation.


    5.1. Effectiveness of contact tracing

    To measure the effectiveness of contact tracing, we use the term epidemic attack rate (AR), which is the cumulative total of exposed/infected individuals during the disease evolution. In figure 2, we plot the attack ratio, ARα/AR0 as a function of the identification delay, α1. Here, AR0 is the mean attack rate when no contact tracing strategy is used and ARα is the mean attack rate for a contact tracing strategy with an identification delay of α1. The plots are the averaged results of 10,000 simulations. Figure 2 shows that contact tracing is more effective in the first and the second scenario when the identification delay, α1, is less than 10 days. In the third scenario, we still observe a reduction in the attack ratio, although not as influential as in the first two. Making the identification rate faster could bring the attack ratio down significantly if contact tracing is started early (<10 days), otherwise, it is not so effective in controlling the epidemic.

    Figure 2. The epidemic attack ratio as a function of α1. The results are the averages of 10,000 simulations.

    In figure 3, we study the impact of hospitalization delay, γ1 on ARγ/AR0. Once again we perform 10,000 simulations and compute the average. Here, no contact tracing strategy has been used. ARγ is the mean attack rate with a hospitalization probability of γ and AR0 is the mean attack rate for the hospitalization probability given in table table 1. The hospitalization delay is the average time it takes to hospitalize an infectious individual. Figure 3 shows the importance of immediate access to hospitals for infected individuals. From the figure, we can see that quick hospitalization can bring down the attack ratio significantly. Comparing figure 3 and figure 2 shows that the epidemic is more sensitive to the hospitalization delay, γ1 than the identification delay, α1. Immediate hospitalization can keep the Ebola epidemic progression under control while contact tracing could be ineffective if it is started once the disease has progressed.

    Figure 3. The epidemic attack ratio as a function of γ1. The results are the averages of 10,000 simulations.

    5.2. Contact tracing performance

    To evaluate the performance of contact tracing on the Ebola contagious process, we employ the ROC curve approach. To plot the ROC curve, we compute the Sensitivity or TPR and (1Specificity) or FPR in each iteration from equations 2 and 3. We compute the average of these two ratios in 10,000 simulations. Figure 4 shows TPR as a function of FPR for the 5 identification delays, α1{1,2,5,10,20} days, in the three contact tracing implementation scenarios. We project the 3-D ROC plot to get a 2-D ROC plot, as shown in figure 5. In all those three scenarios, when α1 increases, the probability to identify secondary infected individuals produced by an infectious individual decreases. Missed-detection probability is the number of exposed individuals who are not traced as defined previously. Therefore, missed-detection is equal to (1Sensitivity). Figure 5 clearly shows that an increase in identification delay can increase the missed-detection probability. Furthermore, early implementation of contact tracing is needed for improving efficiency and reducing cost. It is evident from the figure that the second and third scenarios produce a lot more false positives than the true positives. This implies that the costs of contact tracing could overshadow its benefits if contact tracing is started late.

    Figure 4. 3-D ROC curve of contact tracing with 5 identification delay implemented in three different scenarios. The results are the averages of 10,000 simulations.
    Figure 5. 2-D ROC curve of contact tracing implementation in three different scenarios. The area under the curve (AUC) values for contact tracing starting on day 1, day 9 and day 22 are respectively 0.6550, 0.4060 and 0.1207. The results are the averages of 10,000 simulations.

    5.3. Basic reproductive number

    Using equation 1, we compute the basic reproductive number in a set of 10,000 simulations. In figure 6, we plot the average of the computed R0. It shows that for the same identification rate, α, earlier contact tracing implementations allow better reductions in the value of R0. However, even with the earliest contact tracing strategy, we never obtain a value of R0 smaller than 1. In figure 7, the basic reproductive number is plotted as a function of the hospitalization delay, γ1. Figure 7 shows that a rapid hospitalization strategy can bring the value of R0 below 1. If we can hospitalize infected individuals within 2 days (γ12), the epidemic can be controlled effectively.

    Figure 6. R0 as a function of the identification delay, α1 in three scenarios. The results are the averages of 10,000 simulations.
    Figure 7. R0 as a function of the hospitalization delay, γ1. The results are the averages of 10,000 simulations.

    6. Discussion and conclusion

    In this paper, we have simulated contact tracing on a compartmental model of Ebola in an activity driven network (ADN). We have performed simulations to analyze the effects of contact tracing initiation delay, contact identification delay, and hospitalization delay. Our results suggest that it is critical to start contact tracing within a few days (<10 days), if not immediately after the disease emergence. In addition to that, quick identification of contacts can reduce the epidemic attack ratio up to 50% compared to delayed identification, as shown in figure 2. Contact tracing is very effective when paired with immediate hospitalization in bringing down the reproductive ratio below one.

    The contacts are usually traced up to the maximum incubation period for the disease [27]. The duration of incubation period has counteracting effects on contact tracing. The contacts need to be monitored longer for diseases with longer incubation periods, which is expensive in terms of resources required. When the incubation time is long, we do not gain much by allocating all resources on immediate identification and monitoring of contacts. On the other hand, longer incubation period provides the tracing agencies more time to identify the potential exposed contacts before they become infectious and start to spread pathogens to other people. On this regard, a long incubation period has a positive impact on the disease control. In our model, "identification" and "monitoring" of contacts happen simultaneously. However, these two processes can be separated. A long incubation time will help better identification. However, monitoring individuals for a longer time is expensive. Due to the long incubation period of Ebola, contact tracing could be inefficient. Hence, immediate hospitalization of the infected cases could be a crucial factor in disease control and infection containment.

    A good collaboration between public health authorities and people can lead to rapid identification of secondary infections. Public health authorities should also keep the host population alert and run awareness campaigns to educate people about the disease. Since contact tracing protocols need to monitor all the contacts of the infected and hospitalized people, it cannot separate the exposed contacts from the healthy ones. Therefore, contact tracing might increase the financial burden on the public health authorities when many people are needed to be traced, as happened in the West African countries in 2014. Contact tracing is therefore, economically efficient during the early stage of the epidemic.

    Our work can be extended in the future by modifying the activity potentials based on the states of nodes. Hence, the behavioral changes of individuals in response to the disease can be incorporated and properly modeled. It has been found that these changes have a controlling effect on the epidemic [17]. The model can be made even closer to the reality by using temporal networks with memory (non-Markovian link creation process). There is an alternative approach that can be used with ADNs. Instead of using continuous distribution of activity potential and discrete time-steps, we can use a discrete activity potential distribution with continuous time-steps. This can overcome some limitations of the continuous distribution discrete time ADNs and enhance our capabilities. For example, non-exponential inter-event times can be incorporated, and the partition of nodes in several classes based on their activity potentials can help studying the non mean-field dynamics [25].

    When it comes to controlling invasion of a disease within a population, the best option is to contain the propagation processes in the early stage and at its source. Contact tracing is one of the important strategies towards this goal; and this applies to many diseases besides Ebola. Behavior of the early-stage dynamic process is very different from large scale outbreak. When number of cases are high, the infection process possesses a stable momentum for spreading which in turn makes modeling easier because mean-field assumption applies. When the number of cases are low, however, the infection process is characterized by extreme randomness and dynamic fluctuations. Furthermore, contact network dynamism is very influential. As a result, mean-field models or models based on quenched/averaged contact networks are not viable candidates. Our modeling effort in this paper and the use of ADN framework calls for further work on efficient and accurate modeling of epidemic processes in heterogeneous populations during early stages after introduction of the infection.




    Author contributions



    Kathryn A. Whitehead and Christopher M. Liauw carried out the laboratory work. Kathryn A. Whitehead produced the first draft of the work and oversaw subsequent manuscript preparation. Joels S. T. Wilson-Nieuwenhuis and Anthony J. Slate were involved with manuscript preparation and writing. Ted Deisenroth and Andrea Preuss designed the conceptualisation of the project and funded the work. Joanna Verran oversaw the project.

    Conflicts of interest



    The authors report no conflicts of interest.

    [1] Pawar E (2016) A review article on acrylic PMMA. IOSR J Mech Civil Eng 13: 1-4.
    [2] Ali U, Karim KJBA, Buang NA (2015) A review of the properties and applications of poly (methyl methacrylate) (PMMA). Polymer Rev 55: 678-705. doi: 10.1080/15583724.2015.1031377
    [3] John MJ (2017) Environmental degradation in biocomposites. Biocomposites for High-Performance Applications. Current Barriers and Future Needs Towards Industrial Development Woodhead Publishing, 181-194. doi: 10.1016/B978-0-08-100793-8.00007-7
    [4] Bell GM, Chadwick J (1994) Regulatory controls on biocides in the United Kingdom and restrictions on the use of triorganotin-containing antifouling products. Int Biodeter Biodegr 34: 375-386. doi: 10.1016/0964-8305(94)90095-7
    [5] Rajab FH, Liauw CM, Benson PS, et al. (2018) Picosecond laser treatment production of hierarchical structured stainless steel to reduce bacterial fouling. Food Bioprod Process 109: 29-40. doi: 10.1016/j.fbp.2018.02.009
    [6] Chaky J, Anderson K, Moss M, et al. (2001) Surface hydrophobicity and surface rigidity induce spore germination in Colletotrichum graminicolaPhytopathology 91: 558-564. doi: 10.1094/PHYTO.2001.91.6.558
    [7] Cappitelli F, Sorlini C (2008) Microorganisms attack synthetic polymers in items representing our cultural heritage. Appl Environ Microbiol 74: 564-569. doi: 10.1128/AEM.01768-07
    [8] Whitehead KA, Colligon J, Verran J (2005) Retention of microbial cells in substratum surface features of micrometer and sub-micrometer dimensions. Coll Surf B: Biointerfaces 41: 129-138. doi: 10.1016/j.colsurfb.2004.11.010
    [9] Ma H, Winslow CJ, Logan BE (2008) Spectral force analysis using atomic force microscopy reveals the importance of surface heterogeneity in bacterial and colloid adhesion to engineered surfaces. Coll Surf B: Biointerfaces 62: 232-237. doi: 10.1016/j.colsurfb.2007.10.007
    [10] Whitehead KA, Deisenroth T, Preuss A, et al. (2011) The effect of surface properties on the strength of attachment of fungal spores using AFM perpendicular force measurements. Coll Surf B: Biointerfaces 82: 483-489. doi: 10.1016/j.colsurfb.2010.10.001
    [11] Kipanga PN, Luyten W (2017) Influence of serum and polystyrene plate type on stability of Candida albicans biofilms. J Microbiol Meth 139: 8-11. doi: 10.1016/j.mimet.2017.04.006
    [12] Bigerelle M, Anselme K, Dufresne E, et al. (2002) An unscaled parameter to measure the order of surfaces: a new surface elaboration to increase cells adhesion. Biomol Eng 19: 79-83. doi: 10.1016/S1389-0344(02)00048-5
    [13] Beauvais A, Schmidt C, Guadagnini S, et al. (2007) An extracellular matrix glues together the aerial-grown hyphae of Aspergillus fumigatusCell Microbiol 9: 1588-1600. doi: 10.1111/j.1462-5822.2007.00895.x
    [14] Ramage G, Rajendran R, Gutierrez-Correa M, et al. (2011) Aspergillus biofilms: clinical and industrial significance. FEMS Microbiol Lett 324: 89-97. doi: 10.1111/j.1574-6968.2011.02381.x
    [15] Marques-Calvo MS (2002) In vitro colonization of hydrophilic contact lenses by Aspergillus nigerJ Ind Microbiol Biotech 29: 6-9. doi: 10.1038/sj.jim.7000255
    [16] Filler SG, Sheppard DC (2006) Fungal invasion of normally non-phagocytic host cells. PLOS Pathog 2: e129. doi: 10.1371/journal.ppat.0020129
    [17] Gostinčar C, Ohm RA, Kogej T, et al. (2014) Genome sequencing of four Aureobasidium pullulans varieties: biotechnological potential, stress tolerance, and description of new species. BMC Genomics 15: 549. doi: 10.1186/1471-2164-15-549
    [18] Lotrakul P, Deenarn P, Prasongsuk S, et al. (2009) Isolation of Aureobasidium pullulans from bathroom surfaces and their antifungal activity against some Aspergilli. Afr J Microbiol Res 3: 253-257.
    [19] Zalar P, Novak M, De Hoog GS, et al. (2011) Dishwashers—a man-made ecological niche accommodating human opportunistic fungal pathogens. Fungal Biol 115: 997-1007. doi: 10.1016/j.funbio.2011.04.007
    [20] Arvanitidou M, Kanellou K, Constantinides TC, et al. (1999) The occurrence of fungi in hospital and community potable waters. Lett Appl Microbiol 29: 81-84. doi: 10.1046/j.1365-2672.1999.00583.x
    [21] Shah AA, Hasan F, Hameed A, et al. (2008) Biological degradation of plastics: a comprehensive review. Biotech Adv 26: 246-265. doi: 10.1016/j.biotechadv.2007.12.005
    [22] Stuart MAC, Huck WTS, Genzer J, et al. (2010) Emerging applications of stimuli-responsive polymer materials. Nat Mat 9: 101-113. doi: 10.1038/nmat2614
    [23] Rosenberg M, Kjellerberg S (1986) Hydrophobic interactions: role in bacterial adhesion. Advances in Microbial Ecology Boston: Springer, 353-393. doi: 10.1007/978-1-4757-0611-6_8
    [24] Rezaei F, Abbasi-Firouzjah M, Shokri B (2014) Investigation of antibacterial and wettability behaviours of plasma-modified PMMA films for application in ophthalmology. J Phys D: Appl Phys 47: 085401. doi: 10.1088/0022-3727/47/8/085401
    [25] Sabia R, Stevens HJ, Varner JR (1999) Pitting of a glass-ceramic during polishing with cerium oxide. J Non-Cryst Solids 249: 123-130. doi: 10.1016/S0022-3093(99)00316-6
    [26] Söz CK, Yilgör E, Yilgör I (2015) Influence of the average surface roughness on the formation of superhydrophobic polymer surfaces through spin-coating with hydrophobic fumed silica. Polymer 62: 118-128. doi: 10.1016/j.polymer.2015.02.032
    [27] Nomura T, Minamiura M, Fukamachi K, et al. (2018) Adhesion control of fungal spores on solid surfaces using hydrophilic nanoparticles. Adv Powder Technol 29: 909-914. doi: 10.1016/j.apt.2018.01.007
    [28] Kuo KC, Hoch HC (1996) Germination of Phyllosticta ampelicida Pycnidiospores: Prerequisite of adhesion to the substratum and the relationship of substratum wettability. Fungal Genet Biol 20: 18-29. doi: 10.1006/fgbi.1996.0005
    [29] Hamer JE, Howard RJ, Chumley FG, et al. (1988) A mechanism for surface attachment in spores of a plant pathogenic fungus. Science 239: 288-290. doi: 10.1126/science.239.4837.288
    [30] Amiri A, Cholodowski D, Bompeix G (2005) Adhesion and germination of waterborne and airborne conidia of Penicillium expansum to apple and inert surfaces. Phys Mol Plant Pathol 67: 40-48. doi: 10.1016/j.pmpp.2005.07.003
    [31] Sammonds J, Jaspers MV, Jones EE (2016) Pre-infection processes of Botryosphaeriaceae spp.: adhesion of conidia to different substrata. Plant Pathol 65: 1142-1152. doi: 10.1111/ppa.12485
    [32] Faille C, Jullien C, Fontaine F, et al. (2002) Adhesion of Bacillus spores and Escherichia coli cells to inert surfaces: role of surface hydrophobicity. Can J Microbiol 48: 728-738. doi: 10.1139/w02-063
    [33] Foschino R, Picozzi C, Civardi A, et al. (2003) Comparison of surface sampling methods and cleanability assessment of stainless steel surfaces subjected or not to shot peening. J Food Eng 60: 375-381. doi: 10.1016/S0260-8774(03)00060-8
    [34] Luke B, Faull J, Bateman R (2015) Using particle size analysis to determine the hydrophobicity and suspension of fungal conidia with particular relevance to formulation of biopesticide. Biocontrol Sci Technol 25: 383-398. doi: 10.1080/09583157.2014.979396
    [35] Klis FM, Ram AFJ, De Groot PWJ (2007) A molecular and genomic view of the fungal cell wall. Biology of the fungal cell Heidelberg: Springer, 97-120. doi: 10.1007/978-3-540-70618-2_4
    [36] Temple MJ, Cuskin F, Baslé A, et al. (2017) A Bacteroidetes locus dedicated to fungal 1, 6-β-glucan degradation: unique substrate conformation drives specificity of the key endo-1, 6-β-glucanase. J Biol Chem 292: 10639-10650. doi: 10.1074/jbc.M117.787606
    [37] Rekha MR, Sharma CP (2007) Pullulan as a promising biomaterial for biomedical applications: a perspective. Trends Biomater Artif Organs 20: 116-121.
    [38] Prashantha KVH, Lakshmanb K, Shamalab TR, et al. (2005) Biodegradation of chitosan-graft-polymethylmethacrylate films. Int Biodeter Biodegr 56: 115-120. doi: 10.1016/j.ibiod.2005.06.007
  • This article has been cited by:

    1. Antonia Dobreva, 2021, 2321, 0094-243X, 030006, 10.1063/5.0040098
    2. Zhiwei Zhao, Ateeq Rehman, The Characteristics and Physical Fitness Characteristics of Track and Field Sprint Events Based on Discrete Dynamic Modeling of Complex Systems, 2022, 2022, 1875-905X, 1, 10.1155/2022/1722622
    3. Mohammed Said Souid, Krasimira Ivanova, 2022, 2459, 0094-243X, 030034, 10.1063/5.0083544
  • Reader Comments
  • © 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

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

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

Metrics

Article views(5584) PDF downloads(284) Cited by(6)

Article outline

Figures and Tables

Figures(6)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog