1. Introduction
Methicillin-resistant
staphylococcus aureus (MRSA) is a bacteria resistant to antibiotics, particularly methicillin. It affects the subcutaneous tissue of the skin and often starts as a painful red bump that looks like a pimple or a spider bite. MRSA is difficult to treat due to its resistance to traditional antibiotics used in the treatment of Staph infections [
1,
2]; it can equally be deadly. The occurrence of MRSA started in 1960 just after two years of using it as a viable and effective treatment for methicillin-susceptible
staphylococcus aureus.
MRSA is transmitted when a healthy individual comes in direct contact with an infected person with an open wound or an infected tissue or secretions. It can also be contracted indirectly from contaminated surfaces, objects, or materials contaminated with the bacteria by an infected person; the contaminated fomites can remain infectious for months [
3]. A healthy individual becomes colonized with the bacteria until the MRSA bacteria finds an opportunity to infect the soft tissue of the skin [
4].
Substance abuse, including drug abuse, is a cultural and public health problem affecting many millions of persons in the United States [
5]. “Drug addiction is a complex health disorder characterized by its chronic and relapsing nature” [
6]; it is preventable and treatable, and it is not due to any moral or societal failure or a criminal behavior [
6].
MRSA is endemic among injection drug users (IDUs), its spread in the community started in 1981 and has since persisted in the community [
4]. Some of the research studies on MRSA among IDUs are usually carried out by tracking the registered IDUs [
7,
8,
9,
10]. For example, Fleisch et al. [
9] in their study followed 31 MRSA infected IDUs from August of 1994 to December 1999; they found that 19 of these individuals developed secondary, life-threatening ailments such as septic arthritis, endocarditis, pneumonia, and osteomyelitis. In another study, Binswanger et al. [
8] found that 29 individuals displayed subcutaneous inflammation and infection from a study involving 169 IDUs.
Injection drug users present a unique set of behavioral factors that all accumulate to increase the risk of MRSA transmission and infection. For instance, Kretzschmar [
10] found that 12% of the 494 IDUs in their study were HIV-positive while only 1.5% of 207 non-IDUs were HIV-positive and they concluded that HIV among these IDUs was not likely sexually transmitted [
10]. In recent times, there has been an increase in the proportion of invasive MRSA cases among IDUs due to the use of non-sterile injection drug another risky behavior engaged by IDUs [
11]. The number of IDUs increased from 4.1% in 2011 to 9.2% in 2016 [
11]. Furthermore, IDUs are 16.3 times more likely to develop invasive MRSA infections than others [
11].
Demographic studies of illicit drug use found a higher occurrence of illicit drug uses in large metropolitan populations [
12]. The Department of Health and Human Services defined a large metropolitan area in their 2015 national survey of drug use and health as an area having a population greater than one million [
13]. According to the Department, 8.3% of the individuals living in large metropolitan areas used drugs in the past month of the survey, followed by 8.2% of the populations living in small metropolitan areas with only 4.2% in rural areas [
13]. In 2010 National Institute on Drug Abuse found very large cities (population over 500,000) have a higher lifetime prevalence of drug use at 64.9% of drug users followed by 63.7% in large cities (population of 100,000–500,000), 59.3% in medium cities [
12]. According to the National Survey on Drug Use and Health (NSDUH) survey, the lifetime use of illicit drugs in large and small metropolitan areas was 113,967 while non-metropolitan areas were 16,644. These statistics bring to the light the great need for action in drug abuse interventions in cities, particularly in the metropolitan areas.
The goals in this paper are to determine the impact of movement on the control of MRSA when the communities are isolated and connected, and to also understand the effect of movement on drug rehabilitation programs and intervention as a means of controlling risky behaviors among IDUs when in isolated and connected communities. To achieve our goals, we extend the deterministic model in [
14] to include multiple patches and movement between these patches.
This paper is organized as follows: In
Section 2, we present the MRSA transmission model and carry out a theoretical study of its basic properties in
Appendix A. In
Section 2.2 we estimate the movement parameter values; values of the other model parameters are given in [
14]. In
Section 3, sensitivity analysis was carried out to identify the parameters that have significant impact on the model’s response function: The total number of infectious individuals in the community. In
Section 4, we use the results obtained from the sensitivity analysis to investigate the impact of movement on both disease and behavioral controls in connected and isolated communities. The research discussions and conclusions are given in
Section 5 and
Section 6 respectively.
2. Model Formulation
To model the disease transmission, we consider MRSA transmission in three different areas namely large metropolitan, suburban, and rural areas. We developed a deterministic model for the transmission when these areas are connected. We use the deterministic model developed in [
14] and extend it by including the movement of individuals between these connected areas. The methicillin-resistant
Staphylococcus aureus transmission model among injection drug users divides individuals in the community into three subgroups, non-drug users, low-risk IDUs, and high-risk IDUs. The low-risk IDUs use drugs at a limited to moderate level, while high-risk IDUs are intense to severe drug users [
14]. Each of these population subgroups is further divided into three compartments according to their disease status; namely uncolonized susceptible individuals (
), colonized individuals (
), and infected individuals (
), where
for non-injection drug users, low-risk IDUs, and high-risk IDUs and
for large metropolitan, suburban, and rural areas, respectively. These lead to the following total population.
The MRSA transmission dynamics with IDUs and mobility within three different areas is given by the following system of ordinary differential equations
where
and
for the different areas, i.e., large metropolitan, suburban, and rural communities. The flow diagram for the model is depicted in
Figure 1 for a single patch
i. The associated model variables and parameters are described in
Table 1.
The parameter
is the recruitment rate into the uncolonized non-injection drug users (
) subgroup for a particular area or patch
i. The parameter
is the natural death rate in each subpopulation in patch
i. The force of infection of the non-injection drug users in patch
i is given by
where,
for the different subgroups according to their behavior (i.e., non IDUs, low-and high-risk IDUs) and the parameter
is the probability in patch
i of an uncolonized individual becoming colonized following exposure to the MRSA bacteria through contact with colonized and infected individuals in either the non-injection drug users class or the low-risk class or the high-risk class. Colonized individuals in patch
i are decontaminated at the rate
[
15], after which they move back to the uncolonized class. However, as the disease progresses the colonized non-drug users (
) in patch
i move to the infected class at rate
. Infected individuals (
) recover at the rate
and thus move to the uncolonize class (
). The model does not include recovered compartment since there is no immunity associated with the infection. The parameters and transitions for the low-risk and high-risk populations in patch
i are similarly defined (with the subscript
replaced by
and
respectively).
If non-drug users in patch
i engage in risky behavior resulting in injection drug use, we assume that these individuals will initially start small with less frequent and less risky behavior and drug use; therefore, these individuals will move into the low-risk IDU subgroup from the non-drug user subgroup at rate
. However, these individuals can decrease these risky behaviors either by timely family interventions or by enrolling in drug rehabilitation programs [
7,
16] or by other kinds of interventions; when this happens, they go back to the non-drug user subgroup at the rate
.
Individuals in the low-risk subgroup in patch
i might progress and become high-risk injection drug users due to their persistent engagement in these risky behaviors, thereby moving into the high-risk IDUs subgroup at the rate
. As with the low-risk individual, these individuals may stop injecting drugs when they enter a drug rehabilitation program or utilizing needle exchange programs [
7,
16]. We assume these disease and behavioral processes are not instantaneous but rather a gradual process; hence these individuals due to their behavior first enter the low-risk subgroup at rate
.
Note that patch i refer to either large metropolitan, suburban or rural areas. We will investigate the effect of movement on both behavior and disease controls when these patches are either connected or isolated.
2.1. Reproductive Number
The basic reproductive number,
, is given by the result of the expression
where
is the spectral radius.
The basic reproductive number,
, is the average number of secondary infections produced by one infectious individual introduced into a wholly susceptible population [
17,
18,
19,
20]. If
, MRSA will die out in the community, and if
, MRSA it will persist and become endemic. The epidemiological implication of this result is that it will be possible to eliminate MRSA from the community provided the reproduction number (
) is brought to (and maintained at) value less than one.
The theoretical study of the basic quantitative properties of the model are presented in
Appendix A.1, and the derivation of the expression of basic reproductive number,
, are given in
Appendix A.2.
2.2. Model Parameter Estimation
2.2.1. MRSA Demographic Data from 2008–2013
To parameterize the MRSA model in [
14], we used the 2008–2013 demographic data obtained from the Agency for Healthcare and Research and Quality (AHRQ) [
21]. AHRQ compiles discharge data from hospitals throughout the United States; this data is classified as the International Classification of Diseases (ICD-9). Any patients discharged with the MRSA are recorded with ICD-9 code. To obtain the MRSA data from AHRQ discharge data, we set the demographic targets to look at large metro, large suburb, and rural areas. The large metropolitan population produced the highest number of patients with MRSA listed on their medical records upon discharge.
Figure 2 displays the data in tens of thousands of hospital discharges.
2.2.2. Parameter Estimation
The parameters in the single patch model in [
14] of the MRSA model Equation (
1) was obtained using a reduced model involving only the non-IDUs. The ICD-9 MRSA data for large metro, rural, and suburb regional areas was used to estimate the parameters. The results are given in
Table 2; the model simulation profile and the fitted data are depicted in
Figure 3.
To obtain the parameter estimates for the IDUs subgroups in the full MRSA model Equation (
1) for a single patch, we follow the approach in [
14] and multiplied the non-IDUs subgroup parameter estimates
and
by modifying parameters
and
representing low-risk and high-risk subgroups. These risk factors (modification parameters) listed in
Table 3 enable us to determine an estimate away (either in an increasing or decreasing form) from the non-IDUs parameters. Their derivation is given in [
14]. It is possible to also estimate the risk factors using addiction severity score [
22]. We assume that the values of the parameters in each subgroup are the same for the three different areas we are considering.
2.2.3. Estimation of the Movement Rates and
To estimate the movement rates
and
, we focus on the number of people migrating between three counties in Kansas, namely Shawnee, Douglass, and Jefferson counties. These counties were chosen to represent the Large metropolitan, suburban and rural areas. Although the sizes of these counties are smaller than the size of the metropolitan area defined by the Department of Health and Human Services in
Section 1 which puts the size of a large metro as a million. Nevertheless, the result and conclusion are not expected to deviate from the result obtained if we have used a larger population size. The population in these counties (Shawnee, Douglass, and Jefferson counties) from 2009–2013 according to the Census Bureau 2009–2013 data are 175,136, 110,973, and or 18,947, respectively. This period was used to match with the time range of the MRSA discharged data used to estimate the disease parameters. We have not included 2008 in this time range because the Census Bureau estimate of the gross number migrating for 2008 overlaps with other time range and cannot be tease apart.
Nevertheless, we believe the 2009–2013 estimate is a good estimate for the entire 2008–2013 time range. The number of people migrating between these counties over this period was obtained from the county-to-county 2009–2013 ins-outs-nets-gross migration data from the Census Bureau [
23].
Table 4 shows the number of individual migrating between each county. The migration rate is obtained by dividing the number of individuals migrating by the total population in each county, see
Table 5 and
Figure 4.
Figure 4 shows more people moving to the Shawnee county from Douglas, and Jefferson counties; furthermore, more people move to Shawnee than Douglas from Jefferson. However, fewer people migrate in the opposite direction from the larger counties to Jefferson.
.
3. Sensitivity Analysis
Outputs from deterministic models are governed by their input parameters, which often exhibit some level of uncertainty in the process of their selection. Exploration of the MRSA model Equation (
1) above requires the identification of parameters which significantly affect the number of infected individuals (the key model outputs). To assess the impact of the input parameters uncertainty and the sensitivity of the key model outputs from the numerical simulations to variations in each parameter, we carried out a global sensitivity analysis using Latin Hypercube Sampling (LHS) and partial rank correlation coefficients (PRCC).
The LHS is a stratified sampling technique without replacement; it allows for an efficient analysis of parameter variations across simultaneous uncertainty ranges in each parameter [
24,
25,
26,
27]. The PRCC method measures the strength of the relationship between the parameters and the model outcome. It gives the degree to which each parameter affects the model outcome [
24,
25,
26,
27].
To generate the LHS matrices, we assume a uniform distribution for all the model parameters. We then carried out a total of 1000 simulations (runs) of the models for the LHS matrix, using the parameter values given in
Table 2 and
Table 3 (with ranges varying from
the stated baseline values) and as response functions, the sum of carriers and infected in each region for the isolated case; for the connected case, the response function is the sum of the carriers and infected across the regions. These simulation runs were followed by the parameter ranking using PRCC.
The most significant parameters with the largest impact on the model outcome examined under each analysis have the largest partial rank correlation coefficient (PRCC) values. The parameters with positive PRCC values have a positive impact on the number of infected and increase in their values will lead to an increase in the number of infected. While the parameters with negative PRCC values will reduce the number of infected due to their negative impact on the carriers and infected individuals in the population.
The key parameters with negative PRCC values when the regions are isolated are
and
. While those with positive PRCC values are
, and
. It is interesting to note that the PRCC values are the same in all the patches when the patches are in isolation, see
Figure 5.
Figure 6 depict the PRCC values when the regions are connected. The key parameters with negative PRCC values are
and
. While those with positive PRCC values are
, and
. Under this case, the PRCC values are higher in the large metro area followed by the suburban area, the rural area have the lease values. This is due to the high movement rates into the large metro areas compare to the other two areas.
4. Control Measures
In this section we considered two types of control strategies; the first strategy in
Section 4.1 investigates the impact of controlling an individual’s risky behaviors which promotes injection drug use via the vertical downward movement parameters
and
, and the enrollment in rehabilitation programs via the vertical upward movement parameters
and
. The second strategy in
Section 4.2 investigates the control of the disease in the community using results of the highly effective strategy identified in [
14]. For the two different control strategies either control of risky behavior or disease, we used the intermediate risk factor values
% and
% as was done in [
14].
4.1. Control of Risky Injection Drug Use Behaviors
In this section, we investigate the impact of behavioral change due to vertical downward and upward movements between the sub-groups from risky behaviors that promote (upward movement) and discourages (downward movement) injection drug use in the community. Note, our focus here is to control the risky behaviors and not MRSA transmission within the subgroups of the different areas.
We set the baseline behavioral downward and upward movements parameters as , , , . Since we are interest in investigating the effect of behavioral control, on movement between the patches, we decrease the risky behavior by 2 and double the positive behavior. Thus, we set the behavioral parameters as , , , .
4.1.1. Behavioral Control When Patches Are Isolated
We start by considering the case where the patches are in isolation. We observe in
Figure 7 a reduction in the number of cases in each of the patches. We have normalized the number infected for proper comparison between the patches. We observed further from the figure an equal number of infected individuals regardless of the patch size.
Table 6 gives the proportion of infected individuals at the end of the simulation period.
4.1.2. Behavioral Control in a Patch One at a Time When the Patches are Connected
Under this scenario, the patches are connected, and behavioral control is implemented in a patch one at a time. By behavioral control, we mean discouraging risky behavior and encouraging IDUs to either practice safe needle use or to enter rehabilitation centers. In
Figure 8a–c, the control is in the Large metro area only, and we observed a reduction in the proportion of infected individuals; however, this number is more than when the patch is isolated (compare
Table 6 and
Table 7). In addion, we observed that controlling behavior in large metropolitan area leads to a decrease to some extent in the number of infected individuals in the suburban and rural. However, the effect is stronger in the rural area than the suburban. We observed similar dynamics in
Figure 8d–f when we are controlling behavior in the suburban patch. This, however, is not the case when we control behavior in the rural area; we see a negligible effect in the patches with large population even though more people move into these patches from the rural area, see
Figure 8g–i,
Table 7 gives the proportion of infected at the end of the simulation period.
The negligible effect of the control might be due to the fact that the number of individuals coming in from the rural area to the larger populated areas (large metropolitan and suburban areas) is small in comparison to the overall population size in the larger areas, even though the rate from rural area is higher (see
Figure 4). As a result, their overall influence might be weak. This study have not looked at the effect of super-spreaders or influencers, we have assumed that every individual in the community have the same global effect. This will be an interesting future work to consider.
4.2. Control of MRSA Transmission
In this section, we investigate the impact of controlling MRSA transmission within the subgroups when the patches are isolated and when they are connected. To do this, we use results obtained from the sensitivity analysis (see also [
14]). For the sensitivity analysis we obtained control strategies that reduce
and increase
,
, will impact MRSA in the community, since these parameters have strong positive (
) and negative (
) impact on the number of infected; the same holds for the parameters
,
,
, and
. In [
14], three different strategies were implemented namely low-effectiveness strategy, moderate-effectiveness strategy, and high-effectiveness strategy and it was concluded in [
14] that the high-effectiveness strategy did the best in reducing infection in the community. In this study, we focus on the high-effectiveness strategy and set in each patch the following disease-related parameters as
4.2.1. MRSA Control When the Patches Are Isolated
Considering the scenario when the patches are in isolation, we observe a substantial reduction in the number of cases in each of the patches (see
Figure 9 and
Table 7). The level of reduction is the same in each patch regardless of the patch size, and it’s more than when controlling behavior in the isolated patches case, see
Table 6 and
Table 7 for the proportion infected at the end of the simulation period.
4.2.2. MRSA Control in a Patch One at a Time When the Patches Are Connected
Next, we consider the case where the disease is controlled one at a time when the patches are connected. We see in
Figure 10a substantial reduction in infection level in the individual patch where the control is being implemented. The reduction is not as much as what we see when they are isolated (compare
Table 8 and
Table 9). Furthermore, the effect of the control in either the metro or suburban area is negligible when the control is in either place and slightly more when in the rural area, see
Table 9. For instance, when the control is implemented in the large metro area, we see a smaller reduction in the number of infected in the suburban area compare to the number in the large metro and rural areas. This is because the rate at which people are moving to the metro from the suburban is higher than the rate at which they are moving from the metro to the suburban see
Figure 4, and their population (metro and suburban) is relatively comparable in size. In addition, the rate at which people move to the metro from the rural area is higher than the rate from the metro to the rural; when converted to numbers, we see that more people actually move from the metro to the rural area than those that move from the rural area to the metro; this is the reason why the effect of the control in large metropolitan area is felt more in the rural area. In addition, the control of the disease in the rural area have no effect in the larger patches.
4.3. MRSA and Behavioral Controls
In this section, we investigate the effect of controlling both behavior and disease transmission at the same time when the patches are isolated and when they are connected. We set the behavioral parameters as , , , , and the disease parameters as
4.3.1. MRSA and Behavioral Control in Isolated Patches
We start with the case when the patches are isolated, and we are controlling both behavior and disease transmission. We see in
Figure 11 that the disease is just about eliminated, with very low infection levels. This, of course, happens in all the three isolated patches.
4.3.2. MRSA and Behavioral Control in a Patch One at a Time in Connected Patches
Next, we consider the case when the patches are connected, and both behavior and disease are being controlled in the patches one at a time. We observed under this scenario a significant reduction in the number of infected in both large metropolitan and suburban when control is implemented in either place. The proportion of infected is more than what was obtained when the patches are in isolation. Their effect on each other (large metro and suburban areas) is not as strong as with isolated patches. This effect carries through to the rural area as well with the control in the large metropolitan area having a stronger effect than the control in the suburban area, see
Figure 12 and
Table 10.
On the other hand, when the control is implemented in the rural area, the effect on the large metro and suburban areas is negligible. Although we are able to reduce the number of infection in the rural area under this strategy; the reduction is not as low as what was obtained under the isolated case.
4.4. MRSA and Behavioral Control in Connected Patches
Lastly, under this scenario, we implement both behavior and disease control when the region is connected, and we see in
Figure 13 an elimination of the disease across the three patches.
5. Discussion
In this study, we extended the model developed in [
14] to include the movement of individuals between different patches (metropolitan, suburban, and rural areas). The model further includes the movement of individuals within different subgroups into different patches due to behavioral changes. The goal of the study is to determine how such mobility of individuals within the different patches will affect interventions either those targeting individuals behavior or those aimed at curtailing the disease transmission and spread. Note, that we have not aimed to study the impact of movement on disease spread.
The result shows that when we control risky behaviors either by the use of programs like the clean needle exchange or encourage individuals to enter drug rehabilitation programs, we are able to reduce the disease but not eliminate the disease regardless of the connections between the regions. We see a substantial reduction in isolated patches unlike when the patches are connected, and we were doing a targeted control in one area at a time. These results agree with existing research [
28,
29,
30] which note that only a 25% to 40% partial reduction of HIV transmission risk among IDUs can be obtained from interventions which focus solely on individual behavioral change. However, some combination of individually oriented interventions may further reduce the incidence of HIV among IDUs but not completely eliminate it [
30,
31].
Mobility of individuals plays a vital role in the spread of communicable disease. The movement of individuals to an area increases the likelihood of disease occurrence. Population movement and mixing is one of the critical factors identified in the social structural production of HIV risk associated with drug injecting behaviors [
30]. Furthermore, movement of IDUs within social network significantly increases the likelihood of their risky injection drug use behavior over time [
32]. For instance, transient individuals due to homelessness play a role in the spread of HIV among IDUs as these individuals were more likely to engage in risky behaviors and sharing of needles and going to a shooting gallery than non-transient individuals [
33]. These residential transients often have more than one sleeping locations and a higher chance of syringe sharing [
34].
Hence, it is not surprising to see the results from our sensitivity analysis which showed high sensitivity indices for parameters associated with the large metropolitan areas which have a high flow of individuals into it. Although, our model assumes that all individuals from one local to another have the same mobility rates regardless of their behavior. Recall that the large metropolitan area has a large number of people compared to the rural area. This result implicitly aligns with Hoffmann et al. [
32], even though we have used a deterministic model unlike the network model used by Hoffmann et al.
The mobility of injection drug users is a key mechanism for the introduction of both IDUs and infectious diseases into low prevalence areas [
35,
36,
37]. The disease flow is usually from urban centers into rural areas; these rural localities often lack adequate disease prevention and harm reduction services [
35,
38,
39].
Figure 4 shows the migration rate from the larger Shawnee county to the smaller Jefferson county, a county with fewer residents. Although the rural to urban migration has a higher flow of migrating individuals, we still see a flow from the larger county to the smaller county which will impact the distribution of the disease in the smaller county.
Accessing the effect of the intervention on MRSA transmission in our current study shows a considerable reduction in infection level when each patch implements disease control strategies, see
Figure 10. Since there is this backward flow of individuals into the rural area from the larger populated metropolitan and suburban areas, it is not surprising to see the effect of this control trickling into the rural area leading to a noticeable reduction in the number of infection in the rural areas as risky behavior is being controlled in a large metropolitan and suburban areas. Even though these rural communities do not have infrastructures to deal with the disease spread and rehabilitation [
35,
38,
39]. From our study, we see that the effect of the control in the rural area is hardly felt in the large metro and suburban areas. This might be due to the overall people migrating. However, we see a stronger effect in the rural area when both control strategies (i.e., control behavior and disease) are implemented in one area at a time; first in the metropolitan area then in the suburban area, see
Figure 12. The effect in the metro area is stronger than the suburban area.
Thus, top-down rather than the bottom-up approach in targeting MRSA in the community with IDUs is a better strategy to controlling the disease in rural communities that might be struggling with resources or are in hard to reach places. As the effect of the treatment either behavioral or disease control will trickle down to surrounding smaller communities from larger communities implementing such controls. The same cannot be said for the reverse bottom-up approach as the effect of the interventions is only local to the rural area.
An underlining assumption in this study is that individuals have stable housing structures once they migrate from one locale to the other. As earlier noted, homeless transients often have more than one sleeping locations which might increase and encourage needle sharing practices and behaviors. In future studies, we will investigate the effect of transient mobility due to homelessness and the impact of various interventions on disease transmission and spread.