1
|
Ricardo-Caldera D, Espitia-Pérez L, Avilés-Vergara PA, Benítez ÁJ, Chacón-Pacheco J, Ballesteros-Correa J, Negrette-Oquendo A, Soto-De León S, Tovar-Acero C. Trypanosoma cruzi in domestic and wild mammals in the northeast region of Colombia. Int J Parasitol Parasites Wildl 2024; 24:100940. [PMID: 38708188 PMCID: PMC11067360 DOI: 10.1016/j.ijppaw.2024.100940] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2024] [Revised: 04/23/2024] [Accepted: 04/24/2024] [Indexed: 05/07/2024]
Abstract
Trypanosoma cruzi hosts can serve as a source of infection for animals, vectors, and humans, contributing to the establishment of Chagas disease (CD) in a given area. Traditionally, the Department of Córdoba has not been considered a transmission area for CD; however, the report of several acute cases of Chagas disease highlights the importance of studying the dynamics of disease transmission in this region. This study aimed to detect T. cruzi in domestic and wild mammals in the department of Córdoba. In 2017, a cross-sectional descriptive study was conducted in six villages in two municipalities in the department of Córdoba. Blood samples from dogs living in the zones were collected in EDTA vacutainer tubes for domestic mammals. Wild mammals were collected using Sherman and Tomahawk traps and mist nets in crops and peridomiciles. T. cruzi DNA was detected using the kinetoplast DNA (kDNA) variable region and the tandem repeat satellite region of T. cruzi as molecular targets. We sampled 168 dogs and 146 wild mammals. The detected prevalence of T. cruzi was 6.37%; the TcI lineage was found in D. marsupialis, H. anomalus, and one canine. A specimen of D. marsupialis with TcI and TcII lineages was also identified. T. cruzi DNA was detected in domestic and wild animals in the study area, indicating the circulation of the parasite in peridomestic environments. D. marsupialis may represent an important host in maintaining this region's wild and domestic cycle.
Collapse
Affiliation(s)
- Dina Ricardo-Caldera
- Grupo de investigación en Enfermedades Tropicales y Resistencia Bacteriana, Facultad de Ciencias de la Salud, Universidad del Sinú, Montería, Colombia
| | - Lyda Espitia-Pérez
- Grupo de Investigación Biomédica y Biología Molecular, Facultad de Ciencias de la Salud, Universidad del Sinú, Montería, Colombia
| | - Paula A. Avilés-Vergara
- Grupo de investigación en Enfermedades Tropicales y Resistencia Bacteriana, Facultad de Ciencias de la Salud, Universidad del Sinú, Montería, Colombia
| | - Álvaro J. Benítez
- Grupo de investigación en Enfermedades Tropicales y Resistencia Bacteriana, Facultad de Ciencias de la Salud, Universidad del Sinú, Montería, Colombia
- Núcleo de Pesquisa em Ciências Biológicas & Departamento de Ciências Biológicas, Instituto de Ciências Exatas e Biológicas, Universidade Federal de Ouro Preto, Ouro Preto, MG, Brazil
| | - Julio Chacón-Pacheco
- Grupo de Investigación Biodiversidad Unicórdoba, Universidad de Córdoba, Montería, Colombia
- Laboratorio de Ecología Funcional Unidad de Ecología y Sistemática (UNESIS), Laboratorio de Ecología Funcional, Facultad de Ciencias, Pontificia Universidad Javeriana, Bogotá, Colombia
| | | | - Ana Negrette-Oquendo
- Grupo de investigación en Enfermedades Tropicales y Resistencia Bacteriana, Facultad de Ciencias de la Salud, Universidad del Sinú, Montería, Colombia
| | - Sara Soto-De León
- Grupo de investigación en Enfermedades Tropicales y Resistencia Bacteriana, Facultad de Ciencias de la Salud, Universidad del Sinú, Montería, Colombia
| | - Catalina Tovar-Acero
- Grupo de investigación en Enfermedades Tropicales y Resistencia Bacteriana, Facultad de Ciencias de la Salud, Universidad del Sinú, Montería, Colombia
| |
Collapse
|
2
|
Fiatsonu E, Busselman RE, Ndeffo-Mbah ML. A Scoping Review of Mathematical Models Used to Investigate the Role of Dogs in Chagas Disease Transmission. Animals (Basel) 2023; 13:ani13040555. [PMID: 36830342 PMCID: PMC9951694 DOI: 10.3390/ani13040555] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2023] [Revised: 01/24/2023] [Accepted: 02/02/2023] [Indexed: 02/08/2023] Open
Abstract
Chagas disease is a zoonotic vector-borne disease caused by the parasite Trypanosoma cruzi, which affects a variety of mammalian species across the Americas, including humans and dogs. Mathematical modeling has been widely used to investigate the transmission dynamics and control of vector-borne diseases. We performed a scoping review of mathematical models that investigated the role of dogs in T. cruzi transmission. We identified ten peer-reviewed papers that have explicitly modeled the role of dogs in Chagas transmission dynamics. We discuss the different methods employed in these studies, the different transmission metrics, disease transmission routes, and disease control strategies that have been considered and evaluated. In general, mathematical modeling studies have shown that dogs are not only at high risk of T. cruzi infection but are also major contributors to T. cruzi transmission to humans. Moreover, eliminating infected dogs from households or frequent use of insecticide was shown to be effective for curtailing T. cruzi transmission in both humans and dogs. However, when insecticide spraying is discontinued, T. cruzi infections in dogs were shown to return to their pre-spraying levels. We discuss the challenges and opportunities for future modeling studies to improve our understanding of Chagas disease transmission dynamics and control.
Collapse
|
3
|
Masud M, Islam MH, Kim BN. Understanding the Role of Environmental Transmission on COVID-19 Herd Immunity and Invasion Potential. Bull Math Biol 2022; 84:116. [PMID: 36088430 PMCID: PMC9464060 DOI: 10.1007/s11538-022-01070-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2022] [Accepted: 08/18/2022] [Indexed: 11/28/2022]
Abstract
AbstractCOVID-19 is caused by the SARS-CoV-2 virus, which is mainly transmitted directly between humans. However, it is observed that this disease can also be transmitted through an indirect route via environmental fomites. The development of appropriate and effective vaccines has allowed us to target and anticipate herd immunity. Understanding of the transmission dynamics and the persistence of the virus on environmental fomites and their resistive role on indirect transmission of the virus is an important scientific and public health challenge because it is essential to consider all possible transmission routes and route specific transmission strength to accurately quantify the herd immunity threshold. In this paper, we present a mathematical model that considers both direct and indirect transmission modes. Our analysis focuses on establishing the disease invasion threshold, investigating its sensitivity to both transmission routes and isolate route-specific transmission rate. Using the tau-leap algorithm, we perform a stochastic model simulation to address the invasion potential of both transmission routes. Our analysis shows that direct transmission has a higher invasion potential than that of the indirect transmission. As a proof of this concept, we fitted our model with early epidemic data from several countries to uniquely estimate the reproduction numbers associated with direct and indirect transmission upon confirming the identifiability of the parameters. As the indirect transmission possess lower invasion potential than direct transmission, proper estimation and necessary steps toward mitigating it would help reduce vaccination requirement.
Collapse
Affiliation(s)
- M.A Masud
- Natural Product Informatics Research Center, Korea Institute of Science and Technology, Gangneung, 25451 South Korea
| | - Md. Hamidul Islam
- Department of Applied Mathematics, University of Rajshahi, Rajshahi, 6205 Bangladesh
| | - Byul Nim Kim
- Institute for Mathematical Convergence, Kyungpook National University, Daegu, 41566 South Korea
| |
Collapse
|
4
|
Oizumi R, Inaba H, Takada T, Enatsu Y, Kinjo K. Sensitivity analysis on the declining population in Japan: Effects of prefecture-specific fertility and interregional migration. PLoS One 2022; 17:e0273817. [PMID: 36103457 PMCID: PMC9473415 DOI: 10.1371/journal.pone.0273817] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2022] [Accepted: 08/10/2022] [Indexed: 11/23/2022] Open
Abstract
Japan has been facing a population decline since 2010 due to low birth rates, interregional migration, and regional traits. In this study, we modeled the demographic dynamics of Japan using a transition matrix model. Then, from the mathematical structure of the model, we quantitatively evaluated the domestic factors of population decline. To achieve this, we constructed a multi-regional Leslie matrix model and developed a method for representing the reproductive value and stable age distribution using matrix entries. Our method enabled us to interpret the mathematical indices using the genealogies of the migration history of individuals and their ancestors. Furthermore, by combining our method with sensitivity analysis, we analyzed the effect of region-specific fertility rates and interregional migration rates on the population decline in Japan. We found that the sensitivity of the population growth rate to the migration rate from urban areas with large populations to prefectures with high fertility rates was greatest for people aged under 30. In addition, compared to other areas, the fertility rates of urban areas exhibited higher sensitivity for people aged over 30. Because this feature is robust in comparison with those in 2010 and 2015, it can be said to be a unique structure in Japan in recent years. We also established a method to represent the reproductive value and stable age distribution in an irreducible non-negative matrix population model by using the matrix entries. Furthermore, we show the effects of fertility and migration rates numerically in urban and non-urban areas on the population growth rates for each age group in a society with a declining population.
Collapse
Affiliation(s)
- Ryo Oizumi
- National Institute of Population and Social Security Research, Tokyo, Japan
- * E-mail:
| | - Hisashi Inaba
- Graduate School of Mathematical Science, The University of Tokyo, Tokyo, Japan
| | - Takenori Takada
- Graduate School of Environmental Earth Science, Hokkaido University, Sapporo, Hokkaido, Japan
| | - Youichi Enatsu
- Oshamambe Division, Institute of Arts and Sciences, Tokyo University of Science, Tokyo, Japan
| | - Kensaku Kinjo
- Academic Support Center, Kogakuin University, Hachioji-shi, Tokyo, Japan
| |
Collapse
|
5
|
Hanley BJ, Carstensen M, Walsh DP, Christensen SA, Storm DJ, Booth JG, Guinness J, Them CE, Ahmed MS, Schuler KL. Informing Surveillance through the Characterization of Outbreak Potential of Chronic Wasting Disease in White-Tailed Deer. Ecol Modell 2022. [DOI: 10.1016/j.ecolmodel.2022.110054] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
6
|
Brouwer AF. Why the Spectral Radius? An intuition-building introduction to the basic reproduction number. Bull Math Biol 2022; 84:96. [PMID: 35930076 PMCID: PMC9355935 DOI: 10.1007/s11538-022-01057-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2022] [Accepted: 07/21/2022] [Indexed: 11/26/2022]
Abstract
The basic reproduction number [Formula: see text] is a fundamental concept in mathematical epidemiology and infectious disease modeling. Loosely speaking, it describes the number of people that an infectious person is expected to infect. The basic reproduction number has profound implications for epidemic trajectories and disease control strategies. It is well known that the basic reproduction number can be calculated as the spectral radius of the next generation matrix, but why this is the case may not be intuitively obvious. Here, we walk through how the discrete, next generation process connects to the ordinary differential equation disease system of interest, linearized at the disease-free equilibrium. Then, we use linear algebra to develop a geometric explanation of why the spectral radius of the next generation matrix is an epidemic threshold. Finally, we work through a series of examples that help to build familiarity with the kinds of patterns that arise in parameter combinations produced by the next generation method. This article is intended to help new infectious disease modelers develop intuition for the form and interpretation of the basic reproduction number in their disease systems of interest.
Collapse
Affiliation(s)
- Andrew F Brouwer
- Department of Epidemiology, University of Michigan, 1415 Washington Heights, Ann Arbor, MI, 48109, USA.
| |
Collapse
|
7
|
Kajiwara T, Sasaki T, Otani Y. Global stability of an age-structured infection model in vivo with two compartments and two routes. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2022; 19:11047-11070. [PMID: 36124580 DOI: 10.3934/mbe.2022515] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
In this paper, for an infection age model with two routes, virus-to-cell and cell-to-cell, and with two compartments, we show that the basic reproduction ratio $ R_0 $ gives the threshold of the stability. If $ R_0 > 1 $, the interior equilibrium is unique and globally stable, and if $ R_0 \le 1 $, the disease free equilibrium is globally stable. Some stability results are obtained in previous research, but, for example, a complete proof of the global stability of the disease equilibrium was not shown. We give the proof for all the cases, and show that we can use a type reproduction number for this model.
Collapse
Affiliation(s)
- Tsuyoshi Kajiwara
- Graduate School of Environmental and Life Sciences, Okayama University, Japan
| | - Toru Sasaki
- Faculty of Environmental and Life Science, Okayama University, Japan
| | - Yoji Otani
- School of Engineering, Okayama University, Japan
| |
Collapse
|
8
|
Maity S, Mandal PS. A Comparison of Deterministic and Stochastic Plant-Vector-Virus Models Based on Probability of Disease Extinction and Outbreak. Bull Math Biol 2022; 84:41. [PMID: 35150332 DOI: 10.1007/s11538-022-01001-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2021] [Accepted: 01/24/2022] [Indexed: 11/02/2022]
Abstract
In this investigation, we formulate and analyse a stochastic epidemic model using the continuous-time Markov chain model for the propagation of a vector-borne cassava mosaic disease in a single population. The stochastic model is based upon a pre-existing deterministic plant-vector-virus model. To see how demographic stochasticity affects the vector-borne cassava mosaic disease dynamics, we compare the disease dynamics of both deterministic and stochastic models through disease extinction process. The probability of disease extinction and therefore the major outbreak are estimated analytically using the multitype Galton-Watson branching process (GWbp) approximation. Also, we have found the approximate probabilities of disease extinction numerically based on 30000 sample paths, and it is shown to be good estimate with the calculated probabilities from GWbp approximation. In particular, it is observed that there is a very high probability of disease extinction when the disease is introduced via the infected vectors rather than through infected plants.
Collapse
Affiliation(s)
- Sunil Maity
- Department of Mathematics, NIT Patna, Patna, Bihar, India
| | | |
Collapse
|
9
|
Pandey A, Mideo N, Platt TG. Virulence Evolution of Pathogens That Can Grow in Reservoir Environments. Am Nat 2022; 199:141-158. [DOI: 10.1086/717177] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2022]
|
10
|
Hayes S, Lushasi K, Sambo M, Changalucha J, Ferguson EA, Sikana L, Hampson K, Nouvellet P, Donnelly CA. Understanding the incidence and timing of rabies cases in domestic animals and wildlife in south-east Tanzania in the presence of widespread domestic dog vaccination campaigns. Vet Res 2022; 53:106. [PMID: 36510331 PMCID: PMC9743725 DOI: 10.1186/s13567-022-01121-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2021] [Accepted: 10/17/2022] [Indexed: 12/14/2022] Open
Abstract
The "Zero by 30" strategic plan aims to eliminate human deaths from dog-mediated rabies by 2030 and domestic dog vaccination is a vital component of this strategic plan. In areas where domestic dog vaccination has been implemented, it is important to assess the impact of this intervention. Additionally, understanding temporal and seasonal trends in the incidence of animal rabies cases may assist in optimizing such interventions. Data on the incidence of probable rabies cases in domestic and wild animals were collected between January 2011 and December 2018 in thirteen districts of south-east Tanzania where jackals comprise over 40% of reported rabies cases. Vaccination coverage was estimated over this period, as five domestic dog vaccination campaigns took place in all thirteen districts between 2011 and 2016. Negative binomial generalized linear models were used to explore the impact of domestic dog vaccination on the annual incidence of animal rabies cases, whilst generalized additive models were used to investigate the presence of temporal and/or seasonal trends. Increases in domestic dog vaccination coverage were significantly associated with a decreased incidence of rabies cases in both domestic dogs and jackals. A 35% increase in vaccination coverage was associated with a reduction in the incidence of probable dog rabies cases of between 78.0 and 85.5% (95% confidence intervals ranged from 61.2 to 92.2%) and a reduction in the incidence of probable jackal rabies cases of between 75.3 and 91.2% (95% confidence intervals ranged from 53.0 to 96.1%). A statistically significant common seasonality was identified in the monthly incidence of probable rabies cases in both domestic dogs and jackals with the highest incidence from February to August and lowest incidence from September to January. These results align with evidence supporting the use of domestic dog vaccination as part of control strategies aimed at reducing animal rabies cases in both domestic dogs and jackals in this region. The presence of a common seasonal trend requires further investigation but may have implications for the timing of future vaccination campaigns.
Collapse
Affiliation(s)
- Sarah Hayes
- grid.7445.20000 0001 2113 8111Department of Infectious Disease Epidemiology, Faculty of Medicine, School of Public Health, Imperial College London, London, UK ,grid.4991.50000 0004 1936 8948Department of Statistics, University of Oxford, Oxford, UK
| | - Kennedy Lushasi
- grid.414543.30000 0000 9144 642XIfakara Health Institute, Ifakara, Tanzania ,grid.8756.c0000 0001 2193 314XInstitute of Biodiversity, Animal Health and Comparative Medicine, University of Glasgow, Glasgow, UK ,grid.451346.10000 0004 0468 1595Nelson Mandela African Institution of Science and Technology, Arusha, Tanzania
| | - Maganga Sambo
- grid.414543.30000 0000 9144 642XIfakara Health Institute, Ifakara, Tanzania
| | - Joel Changalucha
- grid.414543.30000 0000 9144 642XIfakara Health Institute, Ifakara, Tanzania ,grid.8756.c0000 0001 2193 314XInstitute of Biodiversity, Animal Health and Comparative Medicine, University of Glasgow, Glasgow, UK
| | - Elaine A. Ferguson
- grid.8756.c0000 0001 2193 314XInstitute of Biodiversity, Animal Health and Comparative Medicine, University of Glasgow, Glasgow, UK
| | - Lwitiko Sikana
- grid.414543.30000 0000 9144 642XIfakara Health Institute, Ifakara, Tanzania ,grid.8756.c0000 0001 2193 314XInstitute of Biodiversity, Animal Health and Comparative Medicine, University of Glasgow, Glasgow, UK
| | - Katie Hampson
- grid.414543.30000 0000 9144 642XIfakara Health Institute, Ifakara, Tanzania ,grid.8756.c0000 0001 2193 314XInstitute of Biodiversity, Animal Health and Comparative Medicine, University of Glasgow, Glasgow, UK
| | - Pierre Nouvellet
- grid.12082.390000 0004 1936 7590School of Life Sciences, University of Sussex, Sussex, UK
| | - Christl A. Donnelly
- grid.7445.20000 0001 2113 8111Department of Infectious Disease Epidemiology, Faculty of Medicine, School of Public Health, Imperial College London, London, UK ,grid.4991.50000 0004 1936 8948Department of Statistics, University of Oxford, Oxford, UK
| |
Collapse
|
11
|
Lushasi K, Hayes S, Ferguson EA, Changalucha J, Cleaveland S, Govella NJ, Haydon DT, Sambo M, Mchau GJ, Mpolya EA, Mtema Z, Nonga HE, Steenson R, Nouvellet P, Donnelly CA, Hampson K. Reservoir dynamics of rabies in south-east Tanzania and the roles of cross-species transmission and domestic dog vaccination. J Appl Ecol 2021; 58:2673-2685. [PMID: 35221371 PMCID: PMC7612421 DOI: 10.1111/1365-2664.13983] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2020] [Accepted: 06/24/2021] [Indexed: 12/27/2022]
Abstract
Understanding the role of different species in the transmission of multi-host pathogens, such as rabies virus, is vital for effective control strategies. Across most of sub-Saharan Africa domestic dogs Canis familiaris are considered the reservoir for rabies, but the role of wildlife has been long debated. Here we explore the multi-host transmission dynamics of rabies across south-east Tanzania.Between January 2011 and July 2019, data on probable rabies cases were collected in the regions of Lindi and Mtwara. Hospital records of animal-bite patients presenting to healthcare facilities were used as sentinels for animal contact tracing. The timing, location and species of probable rabid animals were used to reconstruct transmission trees to infer who infected whom and the relative frequencies of within- and between-species transmission.During the study, 688 probable human rabies exposures were identified, resulting in 47 deaths. Of these exposures, 389 were from domestic dogs (56.5%) and 262 from jackals (38.1%). Over the same period, 549 probable animal rabies cases were traced: 303 in domestic dogs (55.2%) and 221 in jackals (40.3%), with the remainder in domestic cats and other wildlife species.Although dog-to-dog transmission was most commonly inferred (40.5% of transmission events), a third of inferred events involved wildlife-to-wildlife transmission (32.6%), and evidence suggested some sustained transmission chains within jackal populations.A steady decline in probable rabies cases in both humans and animals coincided with the implementation of widespread domestic dog vaccination during the first 6 years of the study. Following the lapse of this program, dog rabies cases began to increase in one of the northernmost districts. Synthesis and applications. In south-east Tanzania, despite a relatively high incidence of rabies in wildlife and evidence of wildlife-to-wildlife transmission, domestic dogs remain essential to the reservoir of infection. Continued dog vaccination alongside improved surveillance would allow a fuller understanding of the role of wildlife in maintaining transmission in this area. Nonetheless, dog vaccination clearly suppressed rabies in both domestic dog and wildlife populations, reducing both public health and conservation risks and, if sustained, has potential to eliminate rabies from this region.
Collapse
Affiliation(s)
- Kennedy Lushasi
- Ifakara Health Institute, Ifakara, Tanzania
- Institute of Biodiversity, Animal Health and Comparative Medicine, University of Glasgow, Glasgow, UK
- Nelson Mandela African Institution of Science and Technology, Arusha, Tanzania
| | - Sarah Hayes
- Department of Infectious Disease Epidemiology, Faculty of Medicine, School of Public Health, Imperial College London
| | - Elaine A. Ferguson
- Institute of Biodiversity, Animal Health and Comparative Medicine, University of Glasgow, Glasgow, UK
| | | | - Sarah Cleaveland
- Institute of Biodiversity, Animal Health and Comparative Medicine, University of Glasgow, Glasgow, UK
| | - Nicodem J. Govella
- Ifakara Health Institute, Ifakara, Tanzania
- Nelson Mandela African Institution of Science and Technology, Arusha, Tanzania
| | - Daniel T. Haydon
- Institute of Biodiversity, Animal Health and Comparative Medicine, University of Glasgow, Glasgow, UK
| | | | - Geofrey J. Mchau
- Ministry of Health, Community Development, Gender, Elderly and Children, Dodoma, Tanzania
| | - Emmanuel A. Mpolya
- Institute of Biodiversity, Animal Health and Comparative Medicine, University of Glasgow, Glasgow, UK
- Nelson Mandela African Institution of Science and Technology, Arusha, Tanzania
| | | | - Hezron E. Nonga
- Ministry of Livestock Development and Fisheries, Dodoma, Tanzania
| | - Rachel Steenson
- Institute of Biodiversity, Animal Health and Comparative Medicine, University of Glasgow, Glasgow, UK
| | | | - Christl A. Donnelly
- Department of Infectious Disease Epidemiology, Faculty of Medicine, School of Public Health, Imperial College London
- Department of Statistics, University of Oxford, Oxford, UK
| | - Katie Hampson
- Ifakara Health Institute, Ifakara, Tanzania
- Institute of Biodiversity, Animal Health and Comparative Medicine, University of Glasgow, Glasgow, UK
| |
Collapse
|
12
|
Spillover, hybridization, and persistence in schistosome transmission dynamics at the human-animal interface. Proc Natl Acad Sci U S A 2021; 118:2110711118. [PMID: 34615712 PMCID: PMC8521685 DOI: 10.1073/pnas.2110711118] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/06/2021] [Indexed: 12/24/2022] Open
Abstract
Zoonotic spillover and hybridization of parasites are major emerging public and veterinary health concerns at the interface of infectious disease biology, evolution, and control. Schistosomiasis is a neglected tropical disease of global importance caused by parasites of the Schistosoma genus, and the Schistosoma spp. system within Africa represents a key example of a system where spillover of animal parasites into human populations has enabled formation of hybrids. Combining model-based approaches and analyses of parasitological, molecular, and epidemiological data from northern Senegal, a region with a high prevalence of schistosome hybrids, we aimed to unravel the transmission dynamics of this complex multihost, multiparasite system. Using Bayesian methods and by estimating the basic reproduction number (R0 ), we evaluate the frequency of zoonotic spillover of Schistosoma bovis from livestock and the potential for onward transmission of hybrid S. bovis × S. haematobium offspring within human populations. We estimate R0 of hybrid schistosomes to be greater than the critical threshold of one (1.76; 95% CI 1.59 to 1.99), demonstrating the potential for hybridization to facilitate spread and establishment of schistosomiasis beyond its original geographical boundaries. We estimate R0 for S. bovis to be greater than one in cattle (1.43; 95% CI 1.24 to 1.85) but not in other ruminants, confirming cattle as the primary zoonotic reservoir. Through longitudinal simulations, we also show that where S. bovis and S. haematobium are coendemic (in livestock and humans respectively), the relative importance of zoonotic transmission is predicted to increase as the disease in humans nears elimination.
Collapse
|
13
|
Kuniya T. Structure of epidemic models: toward further applications in economics. JAPANESE ECONOMIC REVIEW (OXFORD, ENGLAND) 2021; 72:581-607. [PMID: 34483700 PMCID: PMC8405350 DOI: 10.1007/s42973-021-00094-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/30/2021] [Revised: 08/05/2021] [Accepted: 08/12/2021] [Indexed: 06/13/2023]
Abstract
In this paper, we review the structure of various epidemic models in mathematical epidemiology for the future applications in economics. The heterogeneity of population and the generalization of nonlinear terms play important roles in making more elaborate and realistic models. The basic, effective, control and type reproduction numbers have been used to estimate the intensity of epidemic, to evaluate the effectiveness of interventions and to design appropriate interventions. The advanced epidemic models includes the age structure, seasonality, spatial diffusion, mutation and reinfection, and the theory of reproduction numbers has been generalized to them. In particular, the existence of sustained periodic solutions has attracted much interest because they can explain the recurrent waves of epidemic. Although the theory of epidemic models has been developed in decades and the development has been accelerated through COVID-19, it is still difficult to completely answer the uncertainty problem of epidemic models. We would have to mind that there is no single model that can solve all questions and build a scientific attitude to comprehensively understand the results obtained by various researchers from different backgrounds.
Collapse
|
14
|
YANG HYUNMO. ARE THE BEGINNING AND ENDING PHASES OF EPIDEMICS CHARACTERIZED BY THE NEXT GENERATION MATRICES? – A CASE STUDY OF DRUG-SENSITIVE AND RESISTANT TUBERCULOSIS MODEL. J BIOL SYST 2021. [DOI: 10.1142/s0218339021500157] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
In epidemiological modeling, the basic reproduction number is usually defined as being the spectral radius of the next-generation matrix evaluated at the trivial equilibrium. The global stability of the trivial equilibrium point was determined by the left eigenvector associated with that next-generation matrix. More recently, the fraction of susceptible individuals was also obtained from the next generation matrix. The gross reproduction number and the fraction of susceptible individuals were calculated by revisiting the drug-sensitive and resistant tuberculosis model. Hence, the next-generation matrices shed light on the evolution of the dynamics: the beginning of the epidemic via the basic reproduction number and approaching the epidemic’s steady-state via the susceptible individuals’ asymptotic fraction.
Collapse
Affiliation(s)
- HYUN MO YANG
- UNICAMP – IMECC – DMA, Praça Sérgio Buarque de Holanda, 651, CEP 13083-859, Campinas, SP, Brazil
| |
Collapse
|
15
|
Understanding the Spread of COVID-19 in China: Spatial–Temporal Characteristics, Risk Analysis and the Impact of the Quarantine of Hubei Province on the Railway Transportation Network. SUSTAINABILITY 2021. [DOI: 10.3390/su13095163] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
The rapid spread of COVID-19 and its global growth constitutes an international public-health emergency, posing a serious threat to global health, safety, and social economy. In this paper, we systematically studied the temporal and spatial characteristics of COVID-19, infectivity, and the impact of Hubei province’s quarantine on the national railway system on the basis of epidemic and national train data. This study found the following: (1) The overall growth of the epidemic was exponential, and the outbreak of Hubei had a strong spread in the eastern and southern directions. The epidemic was generally more serious in the capital or developed cities in each province, and the epidemic outside Hubei was under control after the imported growth ended. (2) On the basis of analyzing the disturbance of the spread of the epidemic by traffic control, the average incubation period of COVID-19 was approximately 4 days. The ratio of the number of cured people to the number of deaths gradually increased, indicating that, given sufficient medical service, the cure rate can be greatly improved. (3) The quarantine of Hubei had greater impact on cities with higher centrality, especially in the Yangtze River Delta region, and smaller impact on the overall connectivity of the national railway network. For local people, quarantine had great impact on the outflow of local people to neighboring provinces.
Collapse
|
16
|
Hosack GR, Ickowicz A, Hayes KR. Quantifying the risk of vector-borne disease transmission attributable to genetically modified vectors. ROYAL SOCIETY OPEN SCIENCE 2021; 8:201525. [PMID: 33959322 PMCID: PMC8074930 DOI: 10.1098/rsos.201525] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/26/2020] [Accepted: 02/11/2021] [Indexed: 06/12/2023]
Abstract
The relative risk of disease transmission caused by the potential release of transgenic vectors, such as through sterile insect technique or gene drive systems, is assessed with comparison with wild-type vectors. The probabilistic risk framework is demonstrated with an assessment of the relative risk of lymphatic filariasis, malaria and o'nyong'nyong arbovirus transmission by mosquito vectors to human hosts given a released transgenic strain of Anopheles coluzzii carrying a dominant sterile male gene construct. Harm is quantified by a logarithmic loss function that depends on the causal risk ratio, which is a quotient of basic reproduction numbers derived from mathematical models of disease transmission. The basic reproduction numbers are predicted to depend on the number of generations in an insectary colony and the number of backcrosses between the transgenic and wild-type lineages. Analogous causal risk ratios for short-term exposure to a single cohort release are also derived. These causal risk ratios were parametrized by probabilistic elicitations, and updated with experimental data for adult vector mortality. For the wild-type, high numbers of insectary generations were predicted to reduce the number of infectious human cases compared with uncolonized wild-type. Transgenic strains were predicted to produce fewer infectious cases compared with the uncolonized wild-type.
Collapse
Affiliation(s)
- Geoffrey R. Hosack
- Commonwealth Scientific and Industrial Research Organisation, Data61, Hobart, Tasmania, Australia
| | - Adrien Ickowicz
- Commonwealth Scientific and Industrial Research Organisation, Data61, Hobart, Tasmania, Australia
| | - Keith R. Hayes
- Commonwealth Scientific and Industrial Research Organisation, Data61, Hobart, Tasmania, Australia
| |
Collapse
|
17
|
Mahmud AS, Kabir MI, Engø-Monsen K, Tahmina S, Riaz BK, Hossain MA, Khanom F, Rahman MM, Rahman MK, Sharmin M, Hossain DM, Yasmin S, Ahmed MM, Lusha MAF, Buckee CO. Megacities as drivers of national outbreaks: The 2017 chikungunya outbreak in Dhaka, Bangladesh. PLoS Negl Trop Dis 2021; 15:e0009106. [PMID: 33529229 PMCID: PMC7880496 DOI: 10.1371/journal.pntd.0009106] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2020] [Revised: 02/12/2021] [Accepted: 01/04/2021] [Indexed: 01/11/2023] Open
Abstract
BACKGROUND Several large outbreaks of chikungunya have been reported in the Indian Ocean region in the last decade. In 2017, an outbreak occurred in Dhaka, Bangladesh, one of the largest and densest megacities in the world. Population mobility and fluctuations in population density are important drivers of epidemics. Measuring population mobility during outbreaks is challenging but is a particularly important goal in the context of rapidly growing and highly connected cities in low- and middle-income countries, which can act to amplify and spread local epidemics nationally and internationally. METHODS We first describe the epidemiology of the 2017 chikungunya outbreak in Dhaka and estimate incidence using a mechanistic model of chikungunya transmission parametrized with epidemiological data from a household survey. We combine the modeled dynamics of chikungunya in Dhaka, with mobility estimates derived from mobile phone data for over 4 million subscribers, to understand the role of population mobility on the spatial spread of chikungunya within and outside Dhaka during the 2017 outbreak. RESULTS We estimate a much higher incidence of chikungunya in Dhaka than suggested by official case counts. Vector abundance, local demographics, and population mobility were associated with spatial heterogeneities in incidence in Dhaka. The peak of the outbreak in Dhaka coincided with the annual Eid holidays, during which large numbers of people traveled from Dhaka to other parts of the country. We show that travel during Eid likely resulted in the spread of the infection to the rest of the country. CONCLUSIONS Our results highlight the impact of large-scale population movements, for example during holidays, on the spread of infectious diseases. These dynamics are difficult to capture using traditional approaches, and we compare our results to a standard diffusion model, to highlight the value of real-time data from mobile phones for outbreak analysis, forecasting, and surveillance.
Collapse
Affiliation(s)
- Ayesha S. Mahmud
- Department of Demography, University of California, Berkeley, Berkeley, California, United States of America
- Harvard T. H. Chan School of Public Health, Boston, Massachusetts, United States of America
| | - Md. Iqbal Kabir
- National Institute of Preventive and Social Medicine, Dhaka, Bangladesh
- Directorate General of Health Services, Dhaka, Bangladesh
| | | | - Sania Tahmina
- Directorate General of Health Services, Dhaka, Bangladesh
| | | | - Md. Akram Hossain
- National Institute of Preventive and Social Medicine, Dhaka, Bangladesh
| | - Fahmida Khanom
- National Institute of Preventive and Social Medicine, Dhaka, Bangladesh
| | | | | | | | | | | | | | | | - Caroline O. Buckee
- Harvard T. H. Chan School of Public Health, Boston, Massachusetts, United States of America
| |
Collapse
|
18
|
Ives AR, Bozzuto C. Estimating and explaining the spread of COVID-19 at the county level in the USA. Commun Biol 2021; 4:60. [PMID: 33402722 PMCID: PMC7785728 DOI: 10.1038/s42003-020-01609-6] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2020] [Accepted: 12/15/2020] [Indexed: 11/25/2022] Open
Abstract
The basic reproduction number, R0, determines the rate of spread of a communicable disease and therefore gives fundamental information needed to plan public health interventions. Using mortality records, we estimated the rate of spread of COVID-19 among 160 counties and county-aggregates in the USA at the start of the epidemic. We show that most of the high among-county variance is explained by four factors (R2 = 0.70): the timing of outbreak, population size, population density, and spatial location. For predictions of future spread, population density and spatial location are important, and for the latter we show that SARS-CoV-2 strains containing the G614 mutation to the spike gene are associated with higher rates of spread. Finally, the high predictability of R0 allows extending estimates to all 3109 counties in the conterminous 48 states. The high variation of R0 argues for public health policies enacted at the county level for controlling COVID-19.
Collapse
Affiliation(s)
- Anthony R Ives
- Department of Integrative Biology, University of Wisconsin-Madison, Madison, WI, 53706, USA.
| | - Claudio Bozzuto
- Wildlife Analysis GmbH, Oetlisbergstrasse 38, 8053, Zurich, Switzerland
| |
Collapse
|
19
|
The Effect of Demographic Variability and Periodic Fluctuations on Disease Outbreaks in a Vector–Host Epidemic Model. INFECTIOUS DISEASES AND OUR PLANET 2021. [DOI: 10.1007/978-3-030-50826-5_2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
|
20
|
Masud MAB, Ahmed M, Rahman MH. Optimal control for COVID-19 pandemic with quarantine and antiviral therapy. SENSORS INTERNATIONAL 2021; 2:100131. [PMID: 34766063 PMCID: PMC8532375 DOI: 10.1016/j.sintl.2021.100131] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/09/2021] [Revised: 10/09/2021] [Accepted: 10/10/2021] [Indexed: 12/31/2022] Open
Abstract
In the absence of a proper cure for the disease, the recent pandemic caused by COVID-19 has been focused on isolation strategies and government measures to control the disease, such as lockdown, media coverage, and improve public hygiene. Mathematical models can help when these intervention mechanisms find some optimal strategies for controlling the spread of such diseases. We propose a set of nonlinear dynamic systems with optimal strategy including practical measures to limit the spread of the virus and to diagnose and isolate infected people while maintaining consciousness for citizens. We have used Pontryagin's maximum principle and solved our system by the finite difference method. In the end, several numerical simulations have been executed to verify the proposed model using Matlab. Also, we pursued the resilience of the parameters of control of the nonlinear dynamic systems, so that we can easily handle the pandemic situation.
Collapse
|
21
|
Wang X, Zhao XQ. Target reproduction numbers for reaction-diffusion population models. J Math Biol 2020; 81:625-647. [PMID: 32737586 DOI: 10.1007/s00285-020-01523-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2019] [Revised: 02/10/2020] [Indexed: 10/23/2022]
Abstract
A very important population threshold quantity is the target reproduction number, which is a measure of control effort required for a target prevention, intervention or control. This concept, as a generalization of type reproduction number, was first introduced in Shuai et al. (J Math Biol 67:1067-1082, 2013) for nonnegative matrices with immediate applications to compartmental population models of ordinary differential equations. The current paper is devoted to the study of all target reproduction numbers for reaction-diffusion population models with compartmental structure. It turns out that the target reproduction number can be regarded as the basic reproduction number of a modified system, where the state of newborn individuals is limited to the target control set and the offspring from the non-target set is regarded as a part of the transition. In other words, the target reproduction number can be interpreted as the expected number of offspring in a specific target set that a primary newborn individual of the same set would produce during its lifetime. We also characterize the target reproduction number so that it can be easily computed numerically for reaction-diffusion models. At the end, we demonstrate our theoretical observations using two examples.
Collapse
Affiliation(s)
- Xueying Wang
- Department of Mathematics and Statistics, Washington State University, Pullman, WA, 99164, USA
| | - Xiao-Qiang Zhao
- Department of Mathematics and Statistics, Memorial University of Newfoundland, St. John's, NL, A1C 5S7, Canada.
| |
Collapse
|
22
|
Zhao H, Wang L, Oliva SM, Zhu H. Modeling and Dynamics Analysis of Zika Transmission with Limited Medical Resources. Bull Math Biol 2020; 82:99. [PMID: 32705573 DOI: 10.1007/s11538-020-00776-1] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2019] [Accepted: 07/08/2020] [Indexed: 01/12/2023]
Abstract
Zika virus, a reemerging mosquito-borne flavivirus, posed a global public health emergency in 2016. Brazil is the most seriously affected country. Some measures have been implemented to control the Zika transmission, such as spraying mosquitoes, developing vaccines and drugs. However, because of the limited medical resources (LMRs) in the country, not every infected patient can be treated in time when infected with Zika virus. We aim to build a deterministic Zika model by introducing a piecewise smooth treatment recovery rate to research the effect of LMRs on the transmission and control of Zika. For the model without treatment, we analyze the global stability of equilibria. For the model with treatment, the model exhibits complex dynamics. We prove that the model with treatment undergoes backward bifurcation, Hopf bifurcation and Bogdanov-Takens bifurcation of codimension 2. It means that the model with LMRs is sensitive to parameters and initial conditions, which has important significance for control of Zika. We also apply the model to estimate the basic and control reproduction numbers for the Zika transmission by using the data on weekly reported accumulated Zika cases from March 25, 2016, to April 14, 2018, in Brazil.
Collapse
Affiliation(s)
- Hongyong Zhao
- Department of Mathematics, Nanjing University of Aeronautics and Astronautics, Nanjing, 210016, People's Republic of China.
| | - Liping Wang
- Department of Mathematics, Nanjing University of Aeronautics and Astronautics, Nanjing, 210016, People's Republic of China
| | - Sergio Muniz Oliva
- Departamento de Matemática Aplicada, Instituto de Matemática e Estatística, Universidade de São Paulo, Rua do Matão, 1010, Cidade Universitária, São Paulo, SP, CEP 05508-090, Brazil
| | - Huaiping Zhu
- Lamps and Department of Mathematics and Statistics, York University, Toronto, ON, M3J 1P3, Canada
| |
Collapse
|
23
|
Danbaba UA, Garba SM. Stability Analysis and Optimal Control for Yellow Fever Model with Vertical Transmission. INTERNATIONAL JOURNAL OF APPLIED AND COMPUTATIONAL MATHEMATICS 2020; 6:105. [PMID: 32835032 PMCID: PMC7336115 DOI: 10.1007/s40819-020-00860-z] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
In this study, a deterministic model for the transmission dynamics of yellow fever (YF) in a human-mosquito setting in the presence of control measures is constructed and rigorously analyzed. In addition to horizontal transmissions, vertical transmission within mosquito population is incorporated. Analysis of the mosquito-only component of the model shows that the reduced model has a mosquito-extinction equilibrium, which is globally-asymptotically stable whenever the basic offspring number ( N 0 ) is less than unity. The vaccinated and type reproduction numbers of the full-model are computed. Condition for global-asymptotic stability of the disease-free equilibrium of the model whenN 0 > 1 is presented. It is shown that, fractional dosing of YF vaccine does not meet YF vaccination requirements. Optimal control theory is applied to the model to characterize the controls parameters. Using Pontryagin's maximum principle and modified forward-backward sweep technique, the necessary conditions for existence of solutions to the optimal control problem is determined. Numerical simulations of the models to assess the effect of fractional vaccine dosing on the disease dynamics and global sensitivity analysis are presented.
Collapse
Affiliation(s)
- UA Danbaba
- Department of Mathematics and Applied Mathematics, University of Pretoria, Pretoria, 0002 South Africa
| | - SM Garba
- Department of Mathematics and Applied Mathematics, University of Pretoria, Pretoria, 0002 South Africa
| |
Collapse
|
24
|
Collins OC, Duffy KJ. Mathematical Analyses on the Effects of Control Measures for a Waterborne Disease Model with Socioeconomic Conditions. J Comput Biol 2020; 28:19-32. [PMID: 32471315 DOI: 10.1089/cmb.2019.0352] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Waterborne diseases are present major health problems to humanity especially in rural communities where many individuals belong to the lower socioeconomic classes (SECs). The impacts of introducing waterborne disease control measures for such communities are investigated by considering a waterborne disease model. The model is extended by introducing treatment of infected individuals and water purification as control measures. The possible benefits of considering these control measures for the various SECs are investigated. Further analyses show how different degrees of control impact the rate at which waterborne diseases are spread across SECs. The disease control model is validated by using it to study the cholera outbreak in Haiti.
Collapse
Affiliation(s)
| | - Kevin Jan Duffy
- Institute of Systems Science, Durban University of Technology, Durban, South Africa
| |
Collapse
|
25
|
LUPICA ANTONELLA, GUMEL ABBAB, PALUMBO ANNUNZIATA. THE COMPUTATION OF REPRODUCTION NUMBERS FOR THE ENVIRONMENT-HOST-ENVIRONMENT CHOLERA TRANSMISSION DYNAMICS. J BIOL SYST 2020. [DOI: 10.1142/s021833902040001x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
This study presents a new model for the environment-host-environment transmission dynamics of V. cholerae in a community with an interconnected aquatic pond–river water network. For the case when the human host is the sole target of anti-cholera control and the volume of water in the pond is maximum, the disease-free equilibrium of the model is shown to be globally asymptotically stable whenever a certain epidemiological threshold, known as the basic reproduction number [Formula: see text], is less than unity. The epidemiological implication of this result is that cholera can be eliminated from the community if the control strategies implemented can bring (and maintain) [Formula: see text] to a value less than unity. Four scenarios, that represent different interpretations of the role of the V. cholerea pathogen within the environment, were studied. The corresponding basic reproduction numbers were shown to exhibit the same threshold property with respect to the value unity (i.e., if one is less (equal, greater) than unity, then the three others are also less (equal, greater) than unity. Further, it was shown that for the case where anti-cholera control is focused on the human host population, the associated type reproduction number of the model (corresponding to each of the four transmission scenarios considered) is unique. The implication of this result is that the estimate of the effort needed for disease elimination (i.e., the required herd immunity threshold) is unique, regardless of which of the four transmission scenarios is considered. However, when any of the other two bacterial population types in the aquatic environment (i.e., bacterial in the pond or river) is the focus of the control efforts, this study shows that the associated type reproduction number is not unique. Extensive numerical simulations of the model, using a realistic set of parameters from the published literature, show that the community-wide implementation of a strategy that focus on improved water quality, sanitation, and hygiene (known as WASH-only strategy), using the current estimated coverage of 50% and efficacy of 60%, is unable to lead to the elimination of the disease. Such elimination is attainable if the coverage and efficacy are increased (e.g., to 80% and 90%, respectively). Further, elimination can be achieved using a strategy that focuses on oral rehydration therapy and the use of antibiotics to treat the infected humans (i.e., treatment-only strategy) for moderate effectiveness and coverage levels. The combined hybrid WASH-treatment strategy provides far better population-level impact vis a vis disease elimination. This study ranks the three interventions in the following order of population-level effectiveness: combined WASH-treatment, followed by treatment-only and then WASH-only strategy.
Collapse
Affiliation(s)
- ANTONELLA LUPICA
- Department of Mathematics and Computer Sciences, University of Catania, V.le A. Doria 6, 95125 Catania, Italy
- Department of Mathematical and Computer Sciences, Physical Sciences and Earth Sciences, University of Messina, V.le F. D’Alcontres 31, 98166 Messina, Italy
| | - ABBA B. GUMEL
- School of Mathematical and Statistical Sciences, Arizona State University, Tempe, Arizona, USA
- Department of Mathematics and Applied Mathematics, University of Pretoria, Pretoria 0002, South Africa
| | - ANNUNZIATA PALUMBO
- Department of Mathematical and Computer Sciences, Physical Sciences and Earth Sciences, University of Messina, V.le F. D’Alcontres 31, 98166 Messina, Italy
| |
Collapse
|
26
|
Dodero-Rojas E, Ferreira LG, Leite VBP, Onuchic JN, Contessoto VG. Modeling Chikungunya control strategies and Mayaro potential outbreak in the city of Rio de Janeiro. PLoS One 2020; 15:e0222900. [PMID: 31990920 PMCID: PMC6986714 DOI: 10.1371/journal.pone.0222900] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2019] [Accepted: 01/02/2020] [Indexed: 12/15/2022] Open
Abstract
Mosquito-borne diseases have become a significant health issue in many regions around the world. For tropical countries, diseases such as Dengue, Zika, and Chikungunya, became epidemic in the last decades. Health surveillance reports during this period were crucial in providing scientific-based information to guide decision making and resources allocation to control outbreaks. In this work, we perform data analysis of the last Chikungunya epidemics in the city of Rio de Janeiro by applying a compartmental mathematical model. Sensitivity analyses were performed in order to describe the contribution of each parameter to the outbreak incidence. We estimate the "basic reproduction number" for those outbreaks and predict the potential epidemic outbreak of the Mayaro virus. We also simulated several scenarios with different public interventions to decrease the number of infected people. Such scenarios should provide insights about possible strategies to control future outbreaks.
Collapse
Affiliation(s)
- Esteban Dodero-Rojas
- Center for Theoretical Biological Physics, Rice University, Houston, TX, United States of America
- Theoretical and Computational Physics Laboratory, University of Costa Rica, San José, Costa Rica
| | - Luiza G. Ferreira
- Department of Chemistry, Rice University, Houston, TX, United States of America
| | - Vitor B. P. Leite
- Department of Physics, Institute of Biosciences, Letters and Exact Sciences, São Paulo State University - UNESP, São José do Rio Preto, SP, Brazil
| | - José N. Onuchic
- Center for Theoretical Biological Physics, Rice University, Houston, TX, United States of America
- Department of Chemistry, Rice University, Houston, TX, United States of America
- Department of Physics & Astronomy, Rice University, Houston, TX, United States of America
- Department of Biosciences, Rice University, Houston, TX, United States of America
| | - Vinícius G. Contessoto
- Center for Theoretical Biological Physics, Rice University, Houston, TX, United States of America
- Brazilian Biorenewables National Laboratory - LNBR, Brazilian Center for Research in Energy and Materials - CNPEM, Campinas, SP, Brazil
- * E-mail:
| |
Collapse
|
27
|
Harrington PD, Lewis MA. A Next-Generation Approach to Calculate Source-Sink Dynamics in Marine Metapopulations. Bull Math Biol 2020; 82:9. [PMID: 31932972 DOI: 10.1007/s11538-019-00674-1] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2019] [Accepted: 12/02/2019] [Indexed: 11/26/2022]
Abstract
In marine systems, adult populations confined to isolated habitat patches can be connected by larval dispersal. Source-sink theory provides effective tools to quantify the effect of specific habitat patches on the dynamics of connected populations. In this paper, we construct the next-generation matrix for a marine metapopulation and demonstrate how it can be used to calculate the source-sink dynamics of habitat patches. We investigate the effect of environmental variables on the source-sink dynamics and demonstrate how the next-generation matrix can provide useful biological insight into transient as well as asymptotic dynamics of the model.
Collapse
Affiliation(s)
- Peter D Harrington
- Department of Mathematical and Statistical Sciences, University of Alberta, Edmonton, AB, T6G 2G1, Canada.
| | - Mark A Lewis
- Department of Mathematical and Statistical Sciences, University of Alberta, Edmonton, AB, T6G 2G1, Canada
- Department of Biological Sciences, University of Alberta, Edmonton, AB, T6G 2G1, Canada
| |
Collapse
|
28
|
Modelling Vector Transmission and Epidemiology of Co-Infecting Plant Viruses. Viruses 2019; 11:v11121153. [PMID: 31847125 PMCID: PMC6950130 DOI: 10.3390/v11121153] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2019] [Revised: 12/03/2019] [Accepted: 12/06/2019] [Indexed: 12/24/2022] Open
Abstract
Co-infection of plant hosts by two or more viruses is common in agricultural crops and natural plant communities. A variety of models have been used to investigate the dynamics of co-infection which track only the disease status of infected and co-infected plants, and which do not explicitly track the density of inoculative vectors. Much less attention has been paid to the role of vector transmission in co-infection, that is, acquisition and inoculation and their synergistic and antagonistic interactions. In this investigation, a general epidemiological model is formulated for one vector species and one plant species with potential co-infection in the host plant by two viruses. The basic reproduction number provides conditions for successful invasion of a single virus. We derive a new invasion threshold which provides conditions for successful invasion of a second virus. These two thresholds highlight some key epidemiological parameters important in vector transmission. To illustrate the flexibility of our model, we examine numerically two special cases of viral invasion. In the first case, one virus species depends on an autonomous virus for its successful transmission and in the second case, both viruses are unable to invade alone but can co-infect the host plant when prevalence is high.
Collapse
|
29
|
Ito H, Yamamoto T, Morita S. The type-reproduction number of sexually transmitted infections through heterosexual and vertical transmission. Sci Rep 2019; 9:17408. [PMID: 31758036 PMCID: PMC6874657 DOI: 10.1038/s41598-019-53841-8] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2019] [Accepted: 11/06/2019] [Indexed: 01/20/2023] Open
Abstract
Multiple sexually transmitted infections (STIs) have threatened human health for centuries. Most STIs spread not only through sexual (horizontal) transmission but also through mother-to-child (vertical) transmission. In a previous work (Ito et al. 2019), we studied a simple model including heterosexual and mother-to-child transmission and proposed a formulation of the basic reproduction number over generations. In the present study, we improved the model to take into account some factors neglected in the previous work: adult mortality from infection, infant mortality caused by mother-to-child transmission, infertility or stillbirth caused by infection, and recovery with treatment. We showed that the addition of these factors has no essential effect on the theoretical formulation. To study the characteristics of the epidemic threshold, we derived analytical formulas for three type-reproduction numbers for adult men, adult women and juveniles. Our result indicates that if an efficient vaccine exists for a prevalent STI, vaccination of females is more effective for containment of the STI than vaccination of males, because the type-reproduction number for adult men is larger than that for adult women when they are larger than one.
Collapse
Affiliation(s)
- Hiromu Ito
- Department of International Health and Medical Anthropology, Institute of Tropical Medicine, Nagasaki University, Nagasaki, 852-8523, Japan.,Department of Environmental Sciences, Zoology, University of Basel, 4051, Basel, Switzerland
| | - Taro Yamamoto
- Department of International Health and Medical Anthropology, Institute of Tropical Medicine, Nagasaki University, Nagasaki, 852-8523, Japan
| | - Satoru Morita
- Department of Mathematical and Systems Engineering, Shizuoka University, Hamamatsu, Shizuoka, 432-8561, Japan. .,Department of Environment and Energy Systems, Graduate School of Science and Technology, Shizuoka University, Hamamatsu, Shizuoka, 432-8561, Japan.
| |
Collapse
|
30
|
Johnstone-Robertson SP, Diuk-Wasser MA, Davis SA. Incorporating tick feeding behaviour into R 0 for tick-borne pathogens. Theor Popul Biol 2019; 131:25-37. [PMID: 31730874 DOI: 10.1016/j.tpb.2019.10.004] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2018] [Revised: 10/17/2019] [Accepted: 10/22/2019] [Indexed: 11/29/2022]
Abstract
Tick-borne pathogens pose a considerable disease burden in Europe and North America, where increasing numbers of human cases and the emergence of new tick-borne pathogens has renewed interest in resolving the mechanisms underpinning their geographical distribution and abundance. For Borrelia burgdorferi and tick-borne encephalitis (TBE) virus, transmission of infection from one generation of ticks to another occurs when older nymphal ticks infect younger larval ticks feeding on the same host, either indirectly via systemic infection of the vertebrate host or directly when feeding in close proximity. Here, expressions for the basic reproduction number, R0, and the related tick type-reproduction number, T, are derived that account for the observation that larval and nymphal ticks tend to aggregate on the same minority of hosts, a tick feeding behaviour known as co-aggregation. The pattern of tick blood meals is represented as a directed, acyclic, bipartite contact network, with individual vertebrate hosts having in-degree, kin, and out-degree, kout, that respectively represent cumulative counts of nymphal and larval ticks fed over the lifetime of the host. The in- and out-degree are not independent when co-aggregation occurs such that [Formula: see text] where 〈.〉 indicates expected value. When systemic infection in the vertebrate host is the dominant transmission route R02=T, whereas when direct transmission between ticks co-feeding on the same host is dominant then R0=T and the effect of co-aggregation on R0 is more pronounced. Simulations of B. burgdorferi and TBE virus transmission on theoretical tick-mouse contact networks revealed that aggregation and co-aggregation have a synergistic effect on R0 and T, that co-aggregation always increases R0 and T, and that aggregation only increases R0 and T when larvae and nymphs also co-aggregate. Co-aggregation has the greatest absolute effect on R0 and T when the mean larval burden of hosts is high, and the largest relative effect on R0 for pathogens sustained by co-feeding transmission, e.g. TBE virus in Europe, compared with those predominantly spread by systemic infection, e.g. B. burgdorferi. For both pathogens, though, co-aggregation increases the mean number of ticks infected per infectious tick, T, and so too the likelihood of pathogen persistence.
Collapse
Affiliation(s)
| | - Maria A Diuk-Wasser
- Department of Ecology, Evolution, and Environmental Biology, Columbia University, New York, NY, USA
| | - Stephen A Davis
- School of Science, RMIT University, Melbourne, Victoria, Australia
| |
Collapse
|
31
|
A Population Dynamics Model of Mosquito-Borne Disease Transmission, Focusing on Mosquitoes' Biased Distribution and Mosquito Repellent Use. Bull Math Biol 2019; 81:4977-5008. [PMID: 31595380 DOI: 10.1007/s11538-019-00666-1] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2019] [Accepted: 09/24/2019] [Indexed: 01/03/2023]
Abstract
We present an improved mathematical model of population dynamics of mosquito-borne disease transmission. Our model considers the effect of mosquito repellent use and the mosquito's behavior or attraction to the infected human, which cause mosquitoes' biased distribution around the human population. Our analysis of the model clearly shows the existence of thresholds for mosquito repellent efficacy and its utilization rate in the human population with respect to the elimination of mosquito-borne diseases. Further, the results imply that the suppression of mosquito-borne diseases becomes more difficult when the mosquitoes' distribution is biased to a greater extent around the human population.
Collapse
|
32
|
Canessa S, Bozzuto C, Pasmans F, Martel A. Quantifying the burden of managing wildlife diseases in multiple host species. CONSERVATION BIOLOGY : THE JOURNAL OF THE SOCIETY FOR CONSERVATION BIOLOGY 2019; 33:1131-1140. [PMID: 30868671 DOI: 10.1111/cobi.13313] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/26/2018] [Revised: 12/20/2018] [Accepted: 01/24/2019] [Indexed: 06/09/2023]
Abstract
Mitigation of infectious wildlife diseases is especially challenging where pathogens affect communities of multiple host species. Although most ecological studies recognize the challenge posed by multiple-species pathogens, the implications for management are typically assessed only qualitatively. Translating the intuitive understanding that multiple host species are important into practice requires a quantitative assessment of whether and how secondary host species should also be targeted by management and the effort this will require. Using a multiple-species compartmental model, we determined analytically whether and how intensively secondary host species should be managed to prevent outbreaks in focal hosts based on the reproduction number of individual host species and between-species transmission rates. We applied the model to the invasive pathogenic fungus Batrachochytrium salamandrivorans in a 2-host system in northern Europe. Avoiding a disease outbreak in the focal host (fire salamanders [Salamandra salamandra]) was impossible unless management also heavily targeted the secondary host (alpine newts [Ichthyosaura alpestris]). Preventing an outbreak in the community required targeted removal of at least 80% of each species. This proportion increased to 90% in the presence of an environmental reservoir of B. salamandrivorans and when the proportion of individuals removed could not be adjusted for different host species (e.g., when using traps that are not species specific). We recommend the focus of disease-mitigation plans should shift from focal species to the community level and calculate explicitly the management efforts required on secondary host species to move beyond the simple intuitive understanding that multiple host species may all influence the system. Failure to do so may lead to underestimating the magnitude of the effort required and ultimately to suboptimal or futile management attempts.
Collapse
Affiliation(s)
- Stefano Canessa
- Department of Pathology, Bacteriology and Avian Diseases, Faculty of Veterinary Medicine, Ghent University, Salisburylaan 133, 9820, Merelbeke, Belgium
| | - Claudio Bozzuto
- Wildlife Analysis GmbH, Oetlisbergstrasse 38, 8053, Zurich, Switzerland
| | - Frank Pasmans
- Department of Pathology, Bacteriology and Avian Diseases, Faculty of Veterinary Medicine, Ghent University, Salisburylaan 133, 9820, Merelbeke, Belgium
| | - An Martel
- Department of Pathology, Bacteriology and Avian Diseases, Faculty of Veterinary Medicine, Ghent University, Salisburylaan 133, 9820, Merelbeke, Belgium
| |
Collapse
|
33
|
Modelling the Effect of a Novel Autodissemination Trap on the Spread of Dengue in Shah Alam and Malaysia. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2019; 2019:1923479. [PMID: 31481976 PMCID: PMC6701433 DOI: 10.1155/2019/1923479] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/07/2018] [Revised: 04/03/2019] [Accepted: 04/18/2019] [Indexed: 11/17/2022]
Abstract
In this paper, we will start off by introducing the classical Ross-Macdonald model for vector-borne diseases which we use to describe the transmission of dengue between humans and Aedes mosquitoes in Shah Alam, which is a city and the state capital of Selangor, Malaysia. We will focus on analysing the effect of using the Mosquito Home System (MHS), which is an example of an autodissemination trap, in reducing the number of dengue cases by changing the Ross-Macdonald model. By using the national dengue data from Malaysia, we are able to estimate λ, which represents the initial growth rate of the dengue epidemic, and this allows us to estimate the number of mosquitoes in Malaysia. A mathematical expression is also constructed which allows us to estimate the potential number of breeding sites of Aedes mosquitoes. By using the data available from the MHS trial carried out in Section 15 of Shah Alam, we included the potential effect of the MHS into the dengue model and thus modelled the impact MHS has on the spread of dengue within the trial area. We then extended our results to analyse the effect of the MHSs on reducing the number of dengue cases in the whole of Malaysia. A new model was constructed with a basic reproduction number, R 0,Mala MHS, which allows us to identify the required MHSs coverage needed to achieve extinction in Malaysia. Numerical simulations and tables of results were also produced to illustrate our results.
Collapse
|
34
|
Dansu EJ, Seno H. A model for epidemic dynamics in a community with visitor subpopulation. J Theor Biol 2019; 478:115-127. [PMID: 31228488 PMCID: PMC7094103 DOI: 10.1016/j.jtbi.2019.06.020] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2019] [Revised: 06/16/2019] [Accepted: 06/19/2019] [Indexed: 11/04/2022]
Abstract
With a model consisting of SIR and SIS models, we affirm claims in previous works. We derive different basic reproduction numbers looking at varying perspectives. We discuss the biological meanings of these basic reproduction numbers. All the basic reproduction numbers coincide with respect to the critical condition. Relevant public health policies are proposed based on our findings.
With a five dimensional system of ordinary differential equations based on the SIR and SIS models, we consider the dynamics of epidemics in a community which consists of residents and short-stay visitors. Taking different viewpoints to consider public health policies to control the disease, we derive different basic reproduction numbers and clarify their common/different mathematical natures so as to understand their meanings in the dynamics of the epidemic. From our analyses, the short-stay visitor subpopulation could become significant in determining the fate of diseases in the community. Furthermore, our arguments demonstrate that it is necessary to choose one variant of basic reproduction number in order to formulate appropriate public health policies.
Collapse
Affiliation(s)
- Emmanuel J Dansu
- Research Center for Pure and Applied Mathematics, Graduate School of Information Sciences, Tohoku University, Aramaki-Aza-Aoba 6-3-09, Aoba-ku, Sendai 980-8579, Japan.
| | - Hiromi Seno
- Research Center for Pure and Applied Mathematics, Graduate School of Information Sciences, Tohoku University, Aramaki-Aza-Aoba 6-3-09, Aoba-ku, Sendai 980-8579, Japan
| |
Collapse
|
35
|
Bokil VA, Allen LJS, Jeger MJ, Lenhart S. Optimal control of a vectored plant disease model for a crop with continuous replanting. JOURNAL OF BIOLOGICAL DYNAMICS 2019; 13:325-353. [PMID: 31149889 DOI: 10.1080/17513758.2019.1622808] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/02/2018] [Accepted: 04/28/2019] [Indexed: 06/09/2023]
Abstract
Vector-transmitted diseases of plants have had devastating effects on agricultural production worldwide, resulting in drastic reductions in yield for crops such as cotton, soybean, tomato, and cassava. Plant-vector-virus models with continuous replanting are investigated in terms of the effects of selection of cuttings, roguing, and insecticide use on disease prevalence in plants. Previous models are extended to include two replanting strategies: frequencyreplanting and abundance-replanting. In frequency-replanting, replanting of infected cuttings depends on the selection frequency parameter ε, whereas in abundance-replanting, replanting depends on plant abundance via a selection rate parameter also denoted as ε. The two models are analysed and new thresholds for disease elimination are defined for each model. Parameter values for cassava, whiteflies, and African cassava mosaic virus serve as a case study. A numerical sensitivity analysis illustrates how the equilibrium densities of healthy and infected plants vary with parameter values. Optimal control theory is used to investigate the effects of roguing and insecticide use with a goal of maximizing the healthy plants that are harvested. Differences in the control strategies in the two models are seen for large values of ε. Also, the combined strategy of roguing and insecticide use performs better than a single control.
Collapse
Affiliation(s)
- V A Bokil
- a Department of Mathematics, Oregon State University , Corvallis , OR , USA
| | - L J S Allen
- b Department of Mathematics & Statistics, Texas Tech University , Lubbock , TX , USA
| | - M J Jeger
- c Centre for Environmental Policy, Imperial College , London , UK
| | - S Lenhart
- d Department of Mathematics, University of Tennessee , Knoxville , TN , USA
| |
Collapse
|
36
|
A general theory for target reproduction numbers with applications to ecology and epidemiology. J Math Biol 2019; 78:2317-2339. [PMID: 30854577 DOI: 10.1007/s00285-019-01345-4] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2018] [Revised: 02/16/2019] [Indexed: 10/27/2022]
Abstract
A general framework for threshold parameters in population dynamics is developed using the concept of target reproduction numbers. This framework identifies reproduction numbers and other threshold parameters in the literature in terms of their roles in population control. The framework is applied to the analysis of single and multiple control strategies in ecology and epidemiology, and this provides new biological insights.
Collapse
|
37
|
Champagne C, Cazelles B. Comparison of stochastic and deterministic frameworks in dengue modelling. Math Biosci 2019; 310:1-12. [PMID: 30735695 DOI: 10.1016/j.mbs.2019.01.010] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2018] [Revised: 01/28/2019] [Accepted: 01/30/2019] [Indexed: 11/16/2022]
Abstract
We perform estimations of compartment models for dengue transmission in rural Cambodia with increasing complexity regarding both model structure and the account for stochasticity. On the one hand, we successively account for three embedded sources of stochasticity: observation noise, demographic variability and environmental hazard. On the other hand, complexity in the model structure is increased by introducing vector-borne transmission, explicit asymptomatic infections and interacting virus serotypes. Using two sources of case data from dengue epidemics in Kampong Cham (Cambodia), models are estimated in the bayesian framework, with Markov Chain Monte Carlo and Particle Markov Chain Monte Carlo. We highlight the advantages and drawbacks of the different formulations in a practical setting. Although in this case the deterministic models provide a good approximation of the mean trajectory for a low computational cost, the stochastic frameworks better reflect and account for parameter and simulation uncertainty.
Collapse
Affiliation(s)
- Clara Champagne
- Institut de Biologie de l'Ecole Normale Supérieure (IBENS), Ecole Normale Supérieure, CNRS UMR 8197,46 rue d'Ulm, Paris 75005, France; CREST, ENSAE, Université Paris Saclay, 5, avenue Henry Le Chatelier, Palaiseau cedex 91764, France.
| | - Bernard Cazelles
- Institut de Biologie de l'Ecole Normale Supérieure (IBENS), Ecole Normale Supérieure, CNRS UMR 8197,46 rue d'Ulm, Paris 75005, France; International Center for Mathematical and Computational Modeling of Complex Systems (UMMISCO), UMI 209 Sorbonne Université - IRD, Bondy cedex, France
| |
Collapse
|
38
|
Uchiumi Y, Ohtsuki H, Sasaki A. Evolution of self-limited cell division of symbionts. Proc Biol Sci 2019; 286:20182238. [PMID: 30963950 DOI: 10.1098/rspb.2018.2238] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
In mutualism between unicellular hosts and their endosymbionts, symbiont's cell division is often synchronized with its host's, ensuring the permanent relationship between endosymbionts and their hosts. The evolution of synchronized cell division thus has been considered to be an essential step in the evolutionary transition from symbionts to organelles. However, if symbionts would accelerate their cell division without regard for the synchronization with the host, they would proliferate more efficiently. Thus, it is paradoxical that symbionts evolve to limit their own division for synchronized cell division. Here, we theoretically explore the condition for the evolution of self-limited cell division of symbionts, by assuming that symbionts control their division rate and that hosts control symbionts' death rate by intracellular digestion and nutrient supply. Our analysis shows that symbionts can evolve to limit their own cell division. Such evolution occurs if not only symbiont's but also host's benefit through symbiosis is large. Moreover, the coevolution of hosts and symbionts leads to either permanent symbiosis where symbionts proliferate to keep pace with their host, or the arms race between symbionts that behave as lytic parasites and hosts that resist them by rapid digestion.
Collapse
Affiliation(s)
- Yu Uchiumi
- 1 Department of Evolutionary Studies of Biosystems, SOKENDAI (The Graduate University for Advanced Studies) , Hayama, Kanagawa 240-0193 , Japan.,2 Meiji Institute for Advanced Study of Mathematical Sciences, Meiji University , 4-21-1 Nakano, Nakano-ku, Tokyo 164-8525 , Japan
| | - Hisashi Ohtsuki
- 1 Department of Evolutionary Studies of Biosystems, SOKENDAI (The Graduate University for Advanced Studies) , Hayama, Kanagawa 240-0193 , Japan
| | - Akira Sasaki
- 1 Department of Evolutionary Studies of Biosystems, SOKENDAI (The Graduate University for Advanced Studies) , Hayama, Kanagawa 240-0193 , Japan.,3 Evolution and Ecology Program, International Institute for Applied Systems Analysis , Laxenburg 2361 , Austria
| |
Collapse
|
39
|
Minter A, Diggle PJ, Costa F, Childs J, Ko AI, Begon M. A model for leptospire dynamics and control in the Norway rat (Rattus norvegicus) the reservoir host in urban slum environments. Epidemics 2018; 25:26-34. [PMID: 29773482 DOI: 10.1016/j.epidem.2018.05.002] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2017] [Revised: 05/04/2018] [Accepted: 05/04/2018] [Indexed: 12/24/2022] Open
Abstract
Leptospirosis is a zoonosis that humans can contract via contact with animal reservoirs directly or with water contaminated with their urine. The primary reservoir of pathogenic leptospires within urban slum environments is the Norway rat (Rattus norvegicus). Motivated by the annual outbreaks of human leptospirosis in slum urban settings, the within population infection dynamics of the Norway rat were investigated in Pau da Lima, an community in Salvador, Brazil. A mechanistic model of the dynamics of leptospire infection was informed by extensive field and laboratory data was developed and explored analytically. To identify the intraspecific transmission route of most importance, a global sensitivity analysis of the basic reproduction number to its components was performed. In addition, different methods of rodent control were investigated by calculating target reproduction numbers. Our results suggest environmental transmission plays an important role in the maintenance of infection in the rodent population. To control numbers of wild Norway rats, combinations of controls are recommended but environmental control should also be investigated to reduce prevalence of infection in rats.
Collapse
Affiliation(s)
- Amanda Minter
- Institute of Integrative Biology, The University of Liverpool, Liverpool, UK
| | - Peter J Diggle
- CHICAS, Lancaster Medical School, Lancaster University, Lancaster, UK
| | - Federico Costa
- Institute of Integrative Biology, The University of Liverpool, Liverpool, UK; Instituto de Saúde Coletiva, Universidade Federal da Bahia, Salvador, Bahia, Brazil; Centro de Pesquisas Gonçalo Moniz, Fundação Oswaldo Cruz, Ministério da Sáude, Salvador, Bahia, Brazil; Department of Epidemiology of Microbial Diseases, Yale School of Public Health, New Haven, CT, United States
| | - James Childs
- Department of Epidemiology of Microbial Diseases, Yale School of Public Health, New Haven, CT, United States
| | - Albert I Ko
- Centro de Pesquisas Gonçalo Moniz, Fundação Oswaldo Cruz, Ministério da Sáude, Salvador, Bahia, Brazil; Department of Epidemiology of Microbial Diseases, Yale School of Public Health, New Haven, CT, United States
| | - Mike Begon
- Institute of Integrative Biology, The University of Liverpool, Liverpool, UK.
| |
Collapse
|
40
|
Chong KC, Hu P, Lau S, Jia KM, Liang W, Wang MH, Zee BCY, Sun R, Zheng H. Monitoring the age-specificity of measles transmissions during 2009-2016 in Southern China. PLoS One 2018; 13:e0205339. [PMID: 30296273 PMCID: PMC6175510 DOI: 10.1371/journal.pone.0205339] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2018] [Accepted: 09/24/2018] [Indexed: 12/19/2022] Open
Abstract
Background Despite several immunization efforts, China saw a resurgence of measles in 2012. Monitoring of transmissions of individuals from different age groups could offer information that would be valuable for planning adequate disease control strategies. We compared the age-specific effective reproductive numbers (R) of measles during 2009–2016 in Guangdong, China. Methods We estimated the age-specific R values for 7 age groups: 0–8 months, 9–18 months, 19 months to 6 years, 7–15 years, 16–25 years, 26–45 years, and ≥46 years adapting the contact matrix of China. The daily numbers of laboratory and clinically confirmed cases reported to the Center for Disease Control and Prevention of Guangdong were used. Results The peak R values of the entire population were above unity from 2012 to 2016, indicating the persistence of measles in the population. In general, children aged 0–6 years and adults aged 26–45 years had larger values of R when comparing with other age groups after 2012. While the peaks of R values for children aged 0–6 years dropped steadily after 2013, the peaks of R values for adults aged 26–45 years kept at a high range every year. Conclusions Although the provincial supplementary immunization activities (SIAs) conducted in 2009 and 2010 were able to reduce the transmissions from 2009 to 2011, larger values of R for children aged 0–6 years were observed after 2012, indicating that the benefits of the SIAs were short-lived. In addition, the transmissions from adults aged between 26 and 45 years increased over time. Disease control strategies should target children and adult groups that carry high potential for measles transmission.
Collapse
Affiliation(s)
- Ka Chun Chong
- JC School of Public Health and Primary Care, The Chinese University of Hong Kong, Hong Kong, China
- Clinical Trials and Biostatistics Laboratory, Shenzhen Research Institute, The Chinese University of Hong Kong, Hong Kong, China
| | - Pei Hu
- Center for Disease Control and Prevention of Guangdong Province, Guangzhou, China
| | - Steven Lau
- JC School of Public Health and Primary Care, The Chinese University of Hong Kong, Hong Kong, China
| | - Katherine Min Jia
- JC School of Public Health and Primary Care, The Chinese University of Hong Kong, Hong Kong, China
| | - Wenjia Liang
- Center for Disease Control and Prevention of Guangdong Province, Guangzhou, China
| | - Maggie Haitian Wang
- JC School of Public Health and Primary Care, The Chinese University of Hong Kong, Hong Kong, China
- Clinical Trials and Biostatistics Laboratory, Shenzhen Research Institute, The Chinese University of Hong Kong, Hong Kong, China
| | - Benny Chung Ying Zee
- JC School of Public Health and Primary Care, The Chinese University of Hong Kong, Hong Kong, China
- Clinical Trials and Biostatistics Laboratory, Shenzhen Research Institute, The Chinese University of Hong Kong, Hong Kong, China
| | - Riyang Sun
- JC School of Public Health and Primary Care, The Chinese University of Hong Kong, Hong Kong, China
- * E-mail: (HZ); (RS)
| | - Huizhen Zheng
- Center for Disease Control and Prevention of Guangdong Province, Guangzhou, China
- * E-mail: (HZ); (RS)
| |
Collapse
|
41
|
Saldaña F, Barradas I. Control Strategies in Multigroup Models: The Case of the Star Network Topology. Bull Math Biol 2018; 80:2978-3001. [PMID: 30242634 DOI: 10.1007/s11538-018-0503-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2018] [Accepted: 09/05/2018] [Indexed: 11/24/2022]
Abstract
In this paper, we propose control strategies for multigroup epidemic models. We use compartmental [Formula: see text] models to study the dynamics of n host groups sharing the same source of infection in addition to the transmission among members of the same group. In particular, we consider a model for infectious diseases with free-living pathogens in the environment and a metapopulation model with a central patch. We give the detailed derivation of the target reproduction number under three public health interventions and provide the corresponding biological insights. Moreover, using the next-generation approach, we calculate the basic reproduction numbers associated with subsystems of our models and determine algebraic connections to the target reproduction number of the complete model. The analysis presented here illustrates that understanding the topological structure of the infection process and partitioning it into simple cycles is useful to design and evaluate the control strategies.
Collapse
Affiliation(s)
- Fernando Saldaña
- Centro de Investigación en Matemáticas, Apartado Postal 402, C.P. 36240, Guanajuato, Mexico.
| | - Ignacio Barradas
- Centro de Investigación en Matemáticas, Apartado Postal 402, C.P. 36240, Guanajuato, Mexico
| |
Collapse
|
42
|
Bai F, Huff KES, Allen LJS. The effect of delay in viral production in within-host models during early infection. JOURNAL OF BIOLOGICAL DYNAMICS 2018; 13:47-73. [PMID: 30021482 DOI: 10.1080/17513758.2018.1498984] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/19/2017] [Accepted: 06/29/2018] [Indexed: 06/08/2023]
Abstract
Delay in viral production may have a significant impact on the early stages of infection. During the eclipse phase, the time from viral entry until active production of viral particles, no viruses are produced. This delay affects the probability that a viral infection becomes established and timing of the peak viral load. Deterministic and stochastic models are formulated with either multiple latent stages or a fixed delay for the eclipse phase. The deterministic model with multiple latent stages approaches in the limit the model with a fixed delay as the number of stages approaches infinity. The deterministic model framework is used to formulate continuous-time Markov chain and stochastic differential equation models. The probability of a minor infection with rapid viral clearance as opposed to a major full-blown infection with a high viral load is estimated from a branching process approximation of the Markov chain model and the results are confirmed through numerical simulations. In addition, parameter values for influenza A are used to numerically estimate the time to peak viral infection and peak viral load for the deterministic and stochastic models. Although the average length of the eclipse phase is the same in each of the models, as the number of latent stages increases, the numerical results show that the time to viral peak and the peak viral load increase.
Collapse
Affiliation(s)
- Fan Bai
- a Department of Mathematics and Statistics, Texas Tech University , Lubbock , TX , USA
| | - Krystin E S Huff
- a Department of Mathematics and Statistics, Texas Tech University , Lubbock , TX , USA
| | - Linda J S Allen
- a Department of Mathematics and Statistics, Texas Tech University , Lubbock , TX , USA
| |
Collapse
|
43
|
Saad-Roy CM, Ma J, van den Driessche P. The effect of sexual transmission on Zika virus dynamics. J Math Biol 2018; 77:1917-1941. [PMID: 29696377 DOI: 10.1007/s00285-018-1230-1] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2017] [Revised: 03/27/2018] [Indexed: 01/22/2023]
Abstract
Zika virus is a human disease that may lead to neurological disorders in affected individuals, and may be transmitted vectorially (by mosquitoes) or sexually. A mathematical model of Zika virus transmission is formulated, taking into account mosquitoes, sexually active males and females, inactive individuals, and considering both vector transmission and sexual transmission from infectious males to susceptible females. Basic reproduction numbers are computed, and disease control strategies are evaluated. The effect of the incidence function used to model sexual transmission from infectious males to susceptible females is investigated. It is proved that for such functions that are sublinear, if the basic reproduction [Formula: see text], then the disease dies out and [Formula: see text] is a sharp threshold. Moreover, under certain conditions on model parameters and assuming mass action incidence for sexual transmission, it is proved that if [Formula: see text], there exists a unique endemic equilibrium that is globally asymptotically stable. However, under nonlinear incidence, it is shown that for certain functions backward bifurcation and Hopf bifurcation may occur, giving rise to subthreshold equilibria and periodic solutions, respectively. Numerical simulations for various parameter values are displayed to illustrate these behaviours.
Collapse
Affiliation(s)
- C M Saad-Roy
- Department of Mathematics and Statistics, University of Victoria, Victoria, BC, V8W 2Y2, Canada.
| | - Junling Ma
- Department of Mathematics and Statistics, University of Victoria, Victoria, BC, V8W 2Y2, Canada
| | - P van den Driessche
- Department of Mathematics and Statistics, University of Victoria, Victoria, BC, V8W 2Y2, Canada
| |
Collapse
|
44
|
Variation in Wolbachia effects on Aedes mosquitoes as a determinant of invasiveness and vectorial capacity. Nat Commun 2018; 9:1483. [PMID: 29662096 PMCID: PMC5902584 DOI: 10.1038/s41467-018-03981-8] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2017] [Accepted: 03/27/2018] [Indexed: 12/26/2022] Open
Abstract
Wolbachia has been introduced into Aedes aegypti mosquitoes to control the spread of arboviruses, such as dengue, chikungunya and Zika. Studies showed that certain Wolbachia strains (such as wMel) reduce replication of dengue viruses in the laboratory, prompting the release of mosquitoes carrying the bacterium into the field, where vectorial capacity can be realistically assessed in relation to native non-carriers. Here we apply a new analysis to two published datasets, and show that wMel increases the mean and the variance in Ae. aegypti susceptibility to dengue infection when introgressed into Brazil and Vietnam genetic backgrounds. In the absence of other processes, higher mean susceptibility should lead to enhanced viral transmission. The increase in variance, however, widens the basis for selection imposed by unexplored natural forces, retaining the potential for reducing transmission overall. Wolbachia infection in mosquitoes reduces dengue virus spread under specific lab conditions, prompting its use in disease control. Here, King et al. show that Wolbachia increases mean and variance in mosquito susceptibility and explain how this affects Wolbachia invasion and dengue transmission.
Collapse
|
45
|
Exploring contacts facilitating transmission of influenza A(H5N1) virus between poultry farms in West Java, Indonesia: A major role for backyard farms? Prev Vet Med 2018; 156:8-15. [PMID: 29891149 DOI: 10.1016/j.prevetmed.2018.04.008] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2017] [Revised: 03/11/2018] [Accepted: 04/11/2018] [Indexed: 11/21/2022]
Abstract
Highly pathogenic avian influenza virus (HPAIV) H5N1 has been reported in Asia, including Indonesia since 2003. Although several risk factors related to the HPAIV outbreaks in poultry in Indonesia have been identified, little is known of the contact structure of farms of different poultry production types (backyard chickens, broilers, layers, and ducks). This study aims to quantify the contact rates associated with the movement of people, and movements of live birds and products and equipment that affect the risk of HPAIV H5N1 transmission between poultry farms in Indonesia. On 124 poultry farms in 6 districts in West Java, logbooks were distributed to record the movements of farmers/staff and visitors and their poultry contacts. Most movements in backyard chicken, commercial native chicken, broiler and duck farms were visits to and from other poultry farms, whilst in layer farms visits to and from poultry companies, visits to egg collection houses and visit from other poultry farms were most frequent. Over 75% of persons visiting backyard chicken and duck farms had previously visited other poultry farms on the same day. Visitors of backyard chicken farms had the highest average contact rate, either direct contact with poultry on other farms before the visits (1.35 contact/day) or contact during their visits in the farms (10.03 contact/day). These results suggest that backyard chicken farms are most at risk for transmission of HPAIV compared to farms of the other poultry production types. Since visits of farm-to-farm were high, backyard farms could also a potential source for HPAIV transmission to commercial poultry farms.
Collapse
|
46
|
Bürli C, Harbrecht H, Odermatt P, Sayasone S, Chitnis N. Mathematical analysis of the transmission dynamics of the liver fluke, Opisthorchis viverrini. J Theor Biol 2018; 439:181-194. [DOI: 10.1016/j.jtbi.2017.11.020] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2016] [Revised: 11/10/2017] [Accepted: 11/28/2017] [Indexed: 12/31/2022]
|
47
|
Dynamics and control of infections on social networks of population types. Epidemics 2017; 23:11-18. [PMID: 29137859 DOI: 10.1016/j.epidem.2017.10.002] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2017] [Accepted: 10/18/2017] [Indexed: 11/24/2022] Open
Abstract
Random mixing in host populations has been a convenient simplifying assumption in the study of epidemics, but neglects important differences in contact rates within and between population groups. For HIV/AIDS, the assumption of random mixing is inappropriate for epidemics that are concentrated in groups of people at high risk, including female sex workers (FSW) and their male clients (MCF), injecting drug users (IDU) and men who have sex with men (MSM). To find out who transmits infection to whom and how that affects the spread and containment of infection remains a major empirical challenge in the epidemiology of HIV/AIDS. Here we develop a technique, based on the routine sampling of infection in linked population groups (a social network of population types), which shows how an HIV/AIDS epidemic in Can Tho Province of Vietnam began in FSW, was propagated mainly by IDU, and ultimately generated most cases among the female partners of MCF (FPM). Calculation of the case reproduction numbers within and between groups, and for the whole network, provides insights into control that cannot be deduced simply from observations on the prevalence of infection. Specifically, the per capita rate of HIV transmission was highest from FSW to MCF, and most HIV infections occurred in FPM, but the number of infections in the whole network is best reduced by interrupting transmission to and from IDU. This analysis can be used to guide HIV/AIDS interventions using needle and syringe exchange, condom distribution and antiretroviral therapy. The method requires only routine data and could be applied to infections in other populations.
Collapse
|
48
|
van den Driessche P. Reproduction numbers of infectious disease models. Infect Dis Model 2017; 2:288-303. [PMID: 29928743 PMCID: PMC6002118 DOI: 10.1016/j.idm.2017.06.002] [Citation(s) in RCA: 188] [Impact Index Per Article: 26.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2017] [Revised: 06/23/2017] [Accepted: 06/26/2017] [Indexed: 12/29/2022] Open
Abstract
This primer article focuses on the basic reproduction number, ℛ 0 , for infectious diseases, and other reproduction numbers related to ℛ 0 that are useful in guiding control strategies. Beginning with a simple population model, the concept is developed for a threshold value of ℛ 0 determining whether or not the disease dies out. The next generation matrix method of calculating ℛ 0 in a compartmental model is described and illustrated. To address control strategies, type and target reproduction numbers are defined, as well as sensitivity and elasticity indices. These theoretical ideas are then applied to models that are formulated for West Nile virus in birds (a vector-borne disease), cholera in humans (a disease with two transmission pathways), anthrax in animals (a disease that can be spread by dead carcasses and spores), and Zika in humans (spread by mosquitoes and sexual contacts). Some parameter values from literature data are used to illustrate the results. Finally, references for other ways to calculate ℛ 0 are given. These are useful for more complicated models that, for example, take account of variations in environmental fluctuation or stochasticity.
Collapse
|
49
|
Yang HM. Comparison between chikungunya and dengue viruses transmission based on a mathematical model. INT J BIOMATH 2017. [DOI: 10.1142/s1793524517500875] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Chikungunya and dengue viruses are transmitted by mosquitoes of genus Aedes. Based on a mathematical model dealing with arboviruses transmission that encompasses human and mosquito populations, the risks of dengue and chikungunya infections are compared. By the fact that chikungunya virus attains high viral load earlier than dengue virus in both humans and mosquitoes, the potential risk of chikungunya could be higher than the dengue infection. The risk of arboviruses infections is assessed by the reproduction number [Formula: see text], which is obtained by the next generation matrix method and Routh–Hurwitz criteria.
Collapse
Affiliation(s)
- Hyun Mo Yang
- Departamento de Matemática Aplicada, Universidade Estadual de Campinas, Praça Sérgio Buarque de Holanda, 651, CEP 13083-859, Campinas, Sao Paulo, Brazil
| |
Collapse
|
50
|
Abstract
The basic reproduction number R0 is a key parameter in plant disease epidemiology, which largely determines whether or not an epidemic will occur in a plant population. The next generation matrix approach to deriving and calculating the basic reproduction number of a plant virus epidemic is described. The approach is illustrated through a series of examples of increasing complexity, ranging from the simplest case of one vector transmitting one virus to a single host, to the case of multiple vectors, to combined horizontal (vector) and vertical (seed) transmission, and where vector control using insecticides is practised. The importance of parameters representing host and vector population dynamics and their interaction in the absence of disease is stressed, and the constraints these place on the calculation of the basic reproduction number. Finally, mention is made of further elaborations to the approach that could prove useful in plant virus epidemiology.
Collapse
Affiliation(s)
- Frank Van den Bosch
- Computational and Systems Biology, Rothamsted Research, Harpenden AL5 2JQ, United Kingdom
| | - Michael J Jeger
- Centre for Environmental Policy, Imperial College London, Silwood Park, Ascot SL5 7PY, United Kingdom.
| |
Collapse
|