1
|
Stephano MA, Mayengo MM, Irunde JI, Kuznetsov D. Sensitivity analysis and parameters estimation for the transmission of lymphatic filariasis. Heliyon 2023; 9:e20066. [PMID: 37810166 PMCID: PMC10559806 DOI: 10.1016/j.heliyon.2023.e20066] [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: 05/16/2023] [Revised: 08/28/2023] [Accepted: 09/10/2023] [Indexed: 10/10/2023] Open
Abstract
Lymphatic filariasis is a neglected tropical disease which poses public health concern and socio-economic challenges in developing and low-income countries. In this paper, we formulate a deterministic mathematical model for transmission dynamics of lymphatic filariasis to generate data by white noise and use least square method to estimate parameter values. The validity of estimated parameter values is tested by Gaussian distribution method. The residuals of model outputs are normally distributed and hence can be used to study the dynamics of Lymphatic filariasis. After deriving the basic reproduction number, R 0 by the next generation matrix approach, the Partial Rank Correlation Coefficient is employed to explore which parameters significantly affect and most influential to the model outputs. The analysis for equilibrium states shows that the Lymphatic free equilibrium is globally asymptotically stable when the basic reproduction number is less a unity and endemic equilibrium is globally asymptotically stable when R 0 ≥ 1 . The findings reveal that rate of human infection, recruitment rate of mosquitoes increase the average new infections for Lymphatic filariasis. Moreover, asymptomatic individuals contribute significantly in the transmission of Lymphatic filariasis.
Collapse
Affiliation(s)
- Mussa A. Stephano
- School of Computation and Communication Science and Engineering, The Nelson Mandela African Institution of Science and Technology (NM-AIST), P.O.BOX 447 Arusha, Tanzania
- Mkwawa University College of Education, Department of Mathematics, Physics and Informatics, P.O.Box 2513, Iringa, Tanzania
| | - Maranya M. Mayengo
- School of Computation and Communication Science and Engineering, The Nelson Mandela African Institution of Science and Technology (NM-AIST), P.O.BOX 447 Arusha, Tanzania
| | - Jacob I. Irunde
- Mkwawa University College of Education, Department of Mathematics, Physics and Informatics, P.O.Box 2513, Iringa, Tanzania
| | - Dmitry Kuznetsov
- School of Computation and Communication Science and Engineering, The Nelson Mandela African Institution of Science and Technology (NM-AIST), P.O.BOX 447 Arusha, Tanzania
| |
Collapse
|
2
|
On Spatiotemporal Overdispersion and Macroparasite Accumulation in Hosts Leading to Aggregation: A Quantitative Framework. Diseases 2022; 11:diseases11010004. [PMID: 36648869 PMCID: PMC9844341 DOI: 10.3390/diseases11010004] [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: 10/23/2022] [Revised: 12/19/2022] [Accepted: 12/24/2022] [Indexed: 12/29/2022] Open
Abstract
In many host-parasite systems, overdispersion in the distribution of macroparasites leads to parasite aggregation in the host population. This overdispersed distribution is often characterized by the negative binomial or by the power law. The aggregation is shown by a clustering of parasites in certain hosts, while other hosts have few or none. Here, I present a theory behind the overdispersion in complex spatiotemporal systems as well as a computational analysis for tracking the behavior of transmissible diseases with this kind of dynamics. I present a framework where heterogeneity and complexity in host-parasite systems are related to aggregation. I discuss the problem of focusing only on the average parasite burden without observing the risk posed by the associated variance; the consequences of under- or overestimation of disease transmission in a heterogenous system and environment; the advantage of including the network of social interaction in epidemiological modeling; and the implication of overdispersion in the management of health systems during outbreaks.
Collapse
|
3
|
Catano-Lopez A, Rojas-Diaz D, Vélez CM. The Influence of Anthropogenic and Environmental Disturbances on Parameter Estimation of a Dengue Transmission Model. Trop Med Infect Dis 2022; 8:tropicalmed8010005. [PMID: 36668912 PMCID: PMC9861738 DOI: 10.3390/tropicalmed8010005] [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: 08/25/2022] [Revised: 10/29/2022] [Accepted: 12/14/2022] [Indexed: 12/24/2022] Open
Abstract
Some deterministic models deal with environmental conditions and use parameter estimations to obtain experimental parameters, but they do not consider anthropogenic or environmental disturbances, e.g., chemical control or climatic conditions. Even more, they usually use theoretical or measured in-lab parameters without worrying about uncertainties in initial conditions, parameters, or changes in control inputs. Thus, in this study, we estimate parameters (including chemical control parameters) and confidence contours under uncertainty conditions using data from the municipality of Bello (Colombia) during 2010-2014, which includes two epidemic outbreaks. Our study shows that introducing non-periodic pulse inputs into the mathematical model allows us to: (i) perform parameter estimation by fitting real data of consecutive dengue outbreaks, (ii) highlight the importance of chemical control as a method of vector control, and (iii) reproduce the endemic behavior of dengue. We described a methodology for parameter and sub-contour box estimation under uncertainties and performed reliable simulations showing the behavior of dengue spread in different scenarios.
Collapse
|
4
|
Omitting age-dependent mosquito mortality in malaria models underestimates the effectiveness of insecticide-treated nets. PLoS Comput Biol 2022; 18:e1009540. [PMID: 36121847 PMCID: PMC9522293 DOI: 10.1371/journal.pcbi.1009540] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2021] [Revised: 09/29/2022] [Accepted: 08/08/2022] [Indexed: 11/19/2022] Open
Abstract
Mathematical models of vector-borne infections, including malaria, often assume age-independent mortality rates of vectors, despite evidence that many insects senesce. In this study we present survival data on insecticide-resistant Anopheles gambiae s.l. from experiments in Côte d’Ivoire. We fit a constant mortality function and two age-dependent functions (logistic and Gompertz) to the data from mosquitoes exposed (treated) and not exposed (control) to insecticide-treated nets (ITNs), to establish biologically realistic survival functions. This enables us to explore the effects of insecticide exposure on mosquito mortality rates, and the extent to which insecticide resistance might impact the effectiveness of ITNs. We investigate this by calculating the expected number of infectious bites a mosquito will take in its lifetime, and by extension the vectorial capacity. Our results show that the predicted vectorial capacity is substantially lower in mosquitoes exposed to ITNs, despite the mosquitoes in the experiment being highly insecticide-resistant. The more realistic age-dependent functions provide a better fit to the experimental data compared to a constant mortality function and, hence, influence the predicted impact of ITNs on malaria transmission potential. In models with age-independent mortality, there is a great reduction for the vectorial capacity under exposure compared to no exposure. However, the two age-dependent functions predicted an even larger reduction due to exposure, highlighting the impact of incorporating age in the mortality rates. These results further show that multiple exposures to ITNs had a considerable effect on the vectorial capacity. Overall, the study highlights the importance of including age dependency in mathematical models of vector-borne disease transmission and in fully understanding the impact of interventions. Interventions against malaria are most commonly targeted on the adult mosquitoes, which transmit the infection from person to person. One of the most important interventions are bed-nets, treated with insecticides. Unfortunately, extensive exposure of mosquitoes to insecticide has led to widespread evolution of insecticide resistance, which might threaten control strategies. Piecing together the overall impact of resistance on the efficacy of insecticide-treated nets is complex, but can be informed by the use of mathematical models. However, there are some assumptions that the models frequently use which are not realistic in terms of the mosquito biology. In this paper, we formulate a model that includes age-dependent mortality rates, an important parameter in vector control since control strategies most commonly aim to reduce the lifespan of the mosquitoes. By using novel data collected using field-derived insecticide-resistant mosquitoes, we explore the effects that the presence of insecticides on nets have on the mortality rates, as well as the difference incorporating age dependency in the model has on the results. We find that including age-dependent mortality greatly alters the anticipated effects of insecticide-treated nets on mosquito transmission potential, and that ignoring this realism potentially overestimates the negative impact of insecticide resistance.
Collapse
|
5
|
Malhotra A, Wüster W, Owens JB, Hodges CW, Jesudasan A, Ch G, Kartik A, Christopher P, Louies J, Naik H, Santra V, Kuttalam SR, Attre S, Sasa M, Bravo-Vega C, Murray KA. Promoting co-existence between humans and venomous snakes through increasing the herpetological knowledge base. Toxicon X 2021; 12:100081. [PMID: 34522881 PMCID: PMC8426276 DOI: 10.1016/j.toxcx.2021.100081] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2021] [Revised: 07/28/2021] [Accepted: 08/04/2021] [Indexed: 11/23/2022] Open
Abstract
Snakebite incidence at least partly depends on the biology of the snakes involved. However, studies of snake biology have been largely neglected in favour of anthropic factors, with the exception of taxonomy, which has been recognised for some decades to affect the design of antivenoms. Despite this, within-species venom variation and the unpredictability of the correlation with antivenom cross-reactivity has continued to be problematic. Meanwhile, other aspects of snake biology, including behaviour, spatial ecology and activity patterns, distribution, and population demography, which can contribute to snakebite mitigation and prevention, remain underfunded and understudied. Here, we review the literature relevant to these aspects of snakebite and illustrate how demographic, spatial, and behavioural studies can improve our understanding of why snakebites occur and provide evidence for prevention strategies. We identify the large gaps that remain to be filled and urge that, in the future, data and relevant metadata be shared openly via public data repositories so that studies can be properly replicated and data used in future meta-analyses.
Collapse
Affiliation(s)
- Anita Malhotra
- Molecular Ecology and Evolution @ Bangor, School of Natural Sciences, Bangor University, 3rd floor ECW, Deiniol Road, Bangor, LL57 2UW, UK
| | - Wolfgang Wüster
- Molecular Ecology and Evolution @ Bangor, School of Natural Sciences, Bangor University, 3rd floor ECW, Deiniol Road, Bangor, LL57 2UW, UK
| | - John Benjamin Owens
- Molecular Ecology and Evolution @ Bangor, School of Natural Sciences, Bangor University, 3rd floor ECW, Deiniol Road, Bangor, LL57 2UW, UK
- Captive & Field Herpetology Ltd, Wales, 13 Hirfron, Holyhead, Llaingoch, Anglesey, LL65 1YU, UK
| | - Cameron Wesley Hodges
- School of Biology, Institute of Science, Suranaree University of Technology, Muang Nakhon Ratchasima, Thailand
| | - Allwin Jesudasan
- Madras Crocodile Bank Trust, Centre for Herpetology, Post bag No.4, Vadanamelli Village, East Coast Road, Mamallapuram, 603 104, Tamil Nadu, India
| | - Gnaneswar Ch
- Madras Crocodile Bank Trust, Centre for Herpetology, Post bag No.4, Vadanamelli Village, East Coast Road, Mamallapuram, 603 104, Tamil Nadu, India
| | - Ajay Kartik
- Madras Crocodile Bank Trust, Centre for Herpetology, Post bag No.4, Vadanamelli Village, East Coast Road, Mamallapuram, 603 104, Tamil Nadu, India
| | - Peter Christopher
- Madras Crocodile Bank Trust, Centre for Herpetology, Post bag No.4, Vadanamelli Village, East Coast Road, Mamallapuram, 603 104, Tamil Nadu, India
| | | | - Hiral Naik
- School of Animal, Plant and Environmental Sciences, University of the Witwatersrand, Johannesburg. P. O. Wits, 2050, Gauteng, South Africa
- Save the Snakes, R527, Blyderus, Hoedspruit, 1380, South Africa
| | - Vishal Santra
- Captive & Field Herpetology Ltd, Wales, 13 Hirfron, Holyhead, Llaingoch, Anglesey, LL65 1YU, UK
- Society for Nature Conservation, Research and Community Engagement (CONCERN), Nalikul, Hooghly, West Bengal 712407, India
| | - Sourish Rajagopalan Kuttalam
- Society for Nature Conservation, Research and Community Engagement (CONCERN), Nalikul, Hooghly, West Bengal 712407, India
| | - Shaleen Attre
- Durrell Institute of Conservation and Ecology, School of Anthropology and Conservation, Marlowe Building, University of Kent, Canterbury, Kent, CT2 7NR, UK
| | - Mahmood Sasa
- Instituto Clodomiro Picado, Universidad de Costa Rica, San José, Costa Rica
- Escuela de Biología, Universidad de Costa Rica, San José, Costa Rica
| | - Carlos Bravo-Vega
- Research Group in Mathematical and Computational Biology (BIOMAC), Department of Biomedical Engineering, University of the Andes, Bogotá, Colombia
| | - Kris A. Murray
- MRC Centre for Global Infectious Disease Analysis, Imperial College London, UK
- MRC Unit The Gambia at London School of Hygiene and Tropical Medicine, Atlantic Boulevard, Fajara, Gambia
| |
Collapse
|
6
|
Pintor AF, Ray N, Longbottom J, Bravo-Vega CA, Yousefi M, Murray KA, Ediriweera DS, Diggle PJ. Addressing the global snakebite crisis with geo-spatial analyses - Recent advances and future direction. Toxicon X 2021; 11:100076. [PMID: 34401744 PMCID: PMC8350508 DOI: 10.1016/j.toxcx.2021.100076] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2021] [Revised: 07/13/2021] [Accepted: 07/14/2021] [Indexed: 02/08/2023] Open
Abstract
Venomous snakebite is a neglected tropical disease that annually leads to hundreds of thousands of deaths or long-term physical and mental ailments across the developing world. Insufficient data on spatial variation in snakebite risk, incidence, human vulnerability, and accessibility of medical treatment contribute substantially to ineffective on-ground management. There is an urgent need to collect data, fill knowledge gaps and address on-ground management problems. The use of novel, and transdisciplinary approaches that take advantage of recent advances in spatio-temporal models, 'big data', high performance computing, and fine-scale spatial information can add value to snakebite management by strategically improving our understanding and mitigation capacity of snakebite. We review the background and recent advances on the topic of snakebite related geospatial analyses and suggest avenues for priority research that will have practical on-ground applications for snakebite management and mitigation. These include streamlined, targeted data collection on snake distributions, snakebites, envenomings, venom composition, health infrastructure, and antivenom accessibility along with fine-scale models of spatio-temporal variation in snakebite risk and incidence, intraspecific venom variation, and environmental change modifying human exposure. These measures could improve and 'future-proof' antivenom production methods, antivenom distribution and stockpiling systems, and human-wildlife conflict management practices, while simultaneously feeding into research on venom evolution, snake taxonomy, ecology, biogeography, and conservation.
Collapse
Affiliation(s)
- Anna F.V. Pintor
- Division of Data, Analytics and Delivery for Impact (DDI), World Health Organization, Geneva, Switzerland
- Australian Institute of Tropical Health and Medicine, Division of Tropical Health and Medicine, James Cook University, Cairns, Australia
| | - Nicolas Ray
- GeoHealth Group, Institute of Global Health, Faculty of Medicine, University of Geneva, Geneva, Switzerland
- Institute for Environmental Sciences, University of Geneva, Geneva, Switzerland
| | - Joshua Longbottom
- Department of Vector Biology, Liverpool School of Tropical Medicine, Liverpool, United Kingdom
- Centre for Health Informatics, Computing and Statistics, Lancaster Medical School, Lancaster University, Lancaster, United Kingdom
| | - Carlos A. Bravo-Vega
- Research Group in Mathematical and Computational Biology (BIOMAC), Department of Biomedical Engineering, University of Los Andes, Bogotá, Colombia
| | - Masoud Yousefi
- School of Biology, College of Science, University of Tehran, Iran
| | - Kris A. Murray
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, UK
- MRC Unit the Gambia at London School of Hygiene and Tropical Medicine, Atlantic Blvd, Fajara, Gambia
| | - Dileepa S. Ediriweera
- Health Data Science Unit, Faculty of Medicine, University of Kelaniya, Ragama, Sri Lanka
| | - Peter J. Diggle
- Department of Vector Biology, Liverpool School of Tropical Medicine, Liverpool, United Kingdom
| |
Collapse
|
7
|
Variable bites and dynamic populations; new insights in Leishmania transmission. PLoS Negl Trop Dis 2021; 15:e0009033. [PMID: 33493192 PMCID: PMC7861551 DOI: 10.1371/journal.pntd.0009033] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2020] [Revised: 02/04/2021] [Accepted: 12/02/2020] [Indexed: 02/06/2023] Open
Abstract
Leishmaniasis is a neglected tropical disease which kills an estimated 50,000 people each year, with its deadly impact confined mainly to lower to middle income countries. Leishmania parasites are transmitted to human hosts by sand fly vectors during blood feeding. Recent experimental work shows that transmission is modulated by the patchy landscape of infection in the host's skin, and the parasite population dynamics within the vector. Here we assimilate these new findings into a simple probabilistic model for disease transmission which replicates recent experimental results, and assesses their relative importance. The results of subsequent simulations, describing random parasite uptake and dynamics across multiple blood meals, show that skin heterogeneity is important for transmission by short-lived flies, but that for longer-lived flies with multiple bites the population dynamics within the vector dominate transmission probability. Our results indicate that efforts to reduce fly lifespan beneath a threshold of around two weeks may be especially helpful in reducing disease transmission.
Collapse
|
8
|
Núñez-López M, Alarcón Ramos L, Velasco-Hernández JX. Migration rate estimation in an epidemic network. APPLIED MATHEMATICAL MODELLING 2021; 89:1949-1964. [PMID: 32952269 PMCID: PMC7486824 DOI: 10.1016/j.apm.2020.08.025] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/25/2020] [Revised: 07/31/2020] [Accepted: 08/06/2020] [Indexed: 05/07/2023]
Abstract
Most of the recent epidemic outbreaks in the world have as a trigger, a strong migratory component as has been evident in the recent Covid-19 pandemic. In this work we address the problem of migration of human populations and its effect on pathogen reinfections in the case of Dengue, using a Markov-chain susceptible-infected-susceptible (SIS) metapopulation model over a network. Our model postulates a general contact rate that represents a local measure of several factors: the population size of infected hosts that arrive at a given location as a function of total population size, the current incidence at neighboring locations, and the connectivity of the network where the disease spreads. This parameter can be interpreted as an indicator of outbreak risk at a given location. This parameter is tied to the fraction of individuals that move across boundaries (migration). To illustrate our model capabilities, we estimate from epidemic Dengue data in Mexico the dynamics of migration at a regional scale incorporating climate variability represented by an index based on precipitation data.
Collapse
Affiliation(s)
- M Núñez-López
- Department of Mathematics, ITAM Río Hondo 1, Ciudad de México 01080, México
| | - L Alarcón Ramos
- Departamento de Matemáticas Aplicadas y Sistemas, Universidad Autónoma Metropolitana, Cuajimalpa, Av. Vasco de Quiroga 4871, Cuajimalpa de Morelos, 05300, México
| | - J X Velasco-Hernández
- Instituto de Matemáticas, Universidad Nacional Autónoma de México, Boulevard Juriquilla No. 3001, Juriquilla, 76230, México
| |
Collapse
|
9
|
Simelane SM, Mwamtobe PM, Abelman S, Tchuenche JM. A Mathematical Model for the Transmission Dynamics of Lymphatic Filariasis with Intervention Strategies. Acta Biotheor 2020; 68:297-320. [PMID: 31758278 DOI: 10.1007/s10441-019-09370-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2018] [Accepted: 11/11/2019] [Indexed: 10/25/2022]
Abstract
This manuscript considers the transmission dynamics of lymphatic filariasis with some intervention strategies in place. Unlike previously developed models, our model takes into account both the exposed and infected classes in both the human and mosquito populations, respectively. We also consider vaccinated, treated and recovered humans in the presented model. The global dynamics of the proposed model are completely determined by the basic ([Formula: see text]) and effective reproduction numbers ([Formula: see text]). We then use Lyapunov function theory to find the sufficient conditions for global stability of both the disease-free equilibrium and endemic equilibrium. The Lyapunov functions show that when the basic reproduction number is less than or equal to unity, the disease-free equilibrium is globally asymptotically stable, and when it is greater than unity then the endemic equilibrium is also globally asymptotically stable. Finally, numerical simulations are carried out to investigate the effects of the intervention strategies and key parameters to the spread of lymphatic filariasis. The numerical simulations support the analytical results and illustrate possible model behavioral scenarios.
Collapse
|
10
|
Bravo-Vega CA, Cordovez JM, Renjifo-Ibáñez C, Santos-Vega M, Sasa M. Estimating snakebite incidence from mathematical models: A test in Costa Rica. PLoS Negl Trop Dis 2019; 13:e0007914. [PMID: 31790407 PMCID: PMC6907855 DOI: 10.1371/journal.pntd.0007914] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2019] [Revised: 12/12/2019] [Accepted: 11/09/2019] [Indexed: 11/18/2022] Open
Abstract
BACKGROUND Snakebite envenoming is a neglected public health challenge that affects mostly economically deprived communities who inhabit tropical regions. In these regions, snakebite incidence data is not always reliable, and access to health care is scare and heterogeneous. Thus, addressing the problem of snakebite effectively requires an understanding of how spatial heterogeneity in snakebite is associated with human demographics and snakes' distribution. Here, we use a mathematical model to address the determinants of spatial heterogeneity in snakebite and we estimate snakebite incidence in a tropical country such as Costa Rica. METHODS AND FINDINGS We combined a mathematical model that follows the law of mass action, where the incidence is proportional to the exposed human population and the venomous snake population, with a spatiotemporal dataset of snakebite incidence (Data from year 1990 to 2007 for 193 districts) in Costa Rica. This country harbors one of the most dangerous venomous snakes, which is the Terciopelo (Bothrops asper, Garman, 1884). We estimated B. asper distribution using a maximum entropy algorithm, and its abundance was estimated based on field data. Then, the model was adjusted to the data using a lineal regression with the reported incidence. We found a significant positive correlation (R2 = 0.66, p-value < 0.01) between our estimation and the reported incidence, suggesting the model has a good performance in estimating snakebite incidence. CONCLUSIONS Our model underscores the importance of the synergistic effect of exposed population size and snake abundance on snakebite incidence. By combining information from venomous snakes' natural history with census data from rural populations, we were able to estimate snakebite incidence in Costa Rica. The model was able to fit the incidence data at fine administrative scale (district level), which is fundamental for the implementation and planning of management strategies oriented to reduce snakebite burden.
Collapse
Affiliation(s)
- Carlos A. Bravo-Vega
- Research Group in Mathematical and Computational Biology (BIOMAC), Department of biomedical engineering, University of los Andes, Bogotá, Colombia
- * E-mail:
| | - Juan M. Cordovez
- Research Group in Mathematical and Computational Biology (BIOMAC), Department of biomedical engineering, University of los Andes, Bogotá, Colombia
| | | | - Mauricio Santos-Vega
- Research Group in Mathematical and Computational Biology (BIOMAC), Department of biomedical engineering, University of los Andes, Bogotá, Colombia
| | - Mahmood Sasa
- Instituto Clodomiro Picado and Escuela de Biología, Universidad de Costa Rica, San José, Costa Rica
| |
Collapse
|
11
|
Benelli G, Duggan MF. Management of arthropod vector data - Social and ecological dynamics facing the One Health perspective. Acta Trop 2018; 182:80-91. [PMID: 29454734 DOI: 10.1016/j.actatropica.2018.02.015] [Citation(s) in RCA: 73] [Impact Index Per Article: 12.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2018] [Revised: 02/12/2018] [Accepted: 02/12/2018] [Indexed: 12/16/2022]
Abstract
Emerging infectious diseases (EIDs) are spread by direct and/or indirect contacts between a pathogen or parasite and their hosts. Arthropod vectors have evolved as excellent bloodsuckers, providing an elegant transportation mode for a wide number of infectious agents. The nature of pathogen and parasite transfer and the models used to predict how a disease might spread are magnified in complexity when an arthropod vector is part of the disease cycle. One Health is a worldwide strategy for expanding interdisciplinary collaborations and communications in all aspects of health care for humans, animals and the environment. It would benefit from a structured analysis to address vectoring of arthropod-borne diseases as a dynamic transactional process. This review focused on how arthropod vector data can be used to better model and predict zoonotic disease outbreaks. With enhanced knowledge to describe arthropod vector disease transfer, researchers will have a better understanding about how to model disease outbreaks. As public health research evolves to include more social-ecological systems, the roles of society, ecology, epidemiology, pathogen/parasite evolution and animal behavior can be better captured in the research design. Overall, because of more collaborative data collection processes on arthropod vectors, disease modeling can better predict conditions where EIDs will occur.
Collapse
|
12
|
Abstract
BACKGROUND The incidence rate of scrub typhus has been increasing in the Republic of Korea. Previous studies have suggested that this trend may have resulted from the effects of climate change on the transmission dynamics among vectors and hosts, but a clear explanation of the process is still lacking. In this study, we applied mathematical models to explore the potential factors that influence the epidemiology of tsutsugamushi disease. METHODS We developed mathematical models of ordinary differential equations including human, rodent and mite groups. Two models, including simple and complex models, were developed, and all parameters employed in the models were adopted from previous articles that represent epidemiological situations in the Republic of Korea. RESULTS The simulation results showed that the force of infection at the equilibrium state under the simple model was 0.236 (per 100,000 person-months), and that in the complex model was 26.796 (per 100,000 person-months). Sensitivity analyses indicated that the most influential parameters were rodent and mite populations and contact rate between them for the simple model, and trans-ovarian transmission for the complex model. In both models, contact rate between humans and mites is more influential than morality rate of rodent and mite group. CONCLUSION The results indicate that the effect of controlling either rodents or mites could be limited, and reducing the contact rate between humans and mites is more practical and effective strategy. However, the current level of control would be insufficient relative to the growing mite population.
Collapse
Affiliation(s)
- Kyung Duk Min
- Department of Public Health Science, Graduate School of Public Health, Seoul National University, Seoul, Korea
| | - Sung Il Cho
- Department of Public Health Science, Graduate School of Public Health, Seoul National University, Seoul, Korea
- Institute of Health and Environment, Seoul National University, Seoul, Korea.
| |
Collapse
|
13
|
Guzzetta G, Trentini F, Poletti P, Baldacchino FA, Montarsi F, Capelli G, Rizzoli A, Rosà R, Merler S, Melegaro A. Effectiveness and economic assessment of routine larviciding for prevention of chikungunya and dengue in temperate urban settings in Europe. PLoS Negl Trop Dis 2017; 11:e0005918. [PMID: 28892499 PMCID: PMC5608415 DOI: 10.1371/journal.pntd.0005918] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2017] [Revised: 09/21/2017] [Accepted: 08/29/2017] [Indexed: 12/13/2022] Open
Abstract
In the last decades, several European countries where arboviral infections are not endemic have faced outbreaks of diseases such as chikungunya and dengue, initially introduced by infectious travellers from tropical endemic areas and then spread locally via mosquito bites. To keep in check the epidemiological risk, interventions targeted to control vector abundance can be implemented by local authorities. We assessed the epidemiological effectiveness and economic costs and benefits of routine larviciding in European towns with temperate climate, using a mathematical model of Aedes albopictus populations and viral transmission, calibrated on entomological surveillance data collected from ten municipalities in Northern Italy during 2014 and 2015.We found that routine larviciding of public catch basins can limit both the risk of autochthonous transmission and the size of potential epidemics. Ideal larvicide interventions should be timed in such a way to cover the month of July. Optimally timed larviciding can reduce locally transmitted cases of chikungunya by 20% - 33% for a single application (dengue: 18-22%) and up to 43% - 65% if treatment is repeated four times throughout the season (dengue: 31-51%). In larger municipalities (>35,000 inhabitants), the cost of comprehensive larviciding over the whole urban area overcomes potential health benefits related to preventing cases of disease, suggesting the adoption of more localized interventions. Small/medium sized towns with high mosquito abundance will likely have a positive cost-benefit balance. Involvement of private citizens in routine larviciding activities further reduces transmission risks but with disproportionate costs of intervention. International travels and the incidence of mosquito-borne diseases are increasing worldwide, exposing a growing number of European citizens to higher risks of potential outbreaks. Results from this study may support the planning and timing of interventions aimed to reduce the probability of autochthonous transmission as well as the nuisance for local populations living in temperate areas of Europe.
Collapse
Affiliation(s)
| | - Filippo Trentini
- Dondena Centre for Research on Social Dynamics and Public Policy, Bocconi University, Milan, Italy
| | - Piero Poletti
- Fondazione Bruno Kessler, Trento, Italy
- Dondena Centre for Research on Social Dynamics and Public Policy, Bocconi University, Milan, Italy
| | | | - Fabrizio Montarsi
- Laboratory of Parasitology, Istituto Zooprofilattico Sperimentale delle Venezie, Padova, Italy
| | - Gioia Capelli
- Laboratory of Parasitology, Istituto Zooprofilattico Sperimentale delle Venezie, Padova, Italy
| | - Annapaola Rizzoli
- Department of Biodiversity and Molecular Ecology, Fondazione Edmund Mach, San Michele all'Adige (Trento), Italy
| | - Roberto Rosà
- Department of Biodiversity and Molecular Ecology, Fondazione Edmund Mach, San Michele all'Adige (Trento), Italy
| | | | - Alessia Melegaro
- Dondena Centre for Research on Social Dynamics and Public Policy, Bocconi University, Milan, Italy
- Department of Policy Analysis and Public Management, Bocconi University, Milan, Italy
| |
Collapse
|
14
|
Mwamtobe PM, Simelane SM, Abelman S, Tchuenche JM. Mathematical analysis of a lymphatic filariasis model with quarantine and treatment. BMC Public Health 2017; 17:265. [PMID: 28302096 PMCID: PMC5356380 DOI: 10.1186/s12889-017-4160-8] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2015] [Accepted: 03/02/2017] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Lymphatic filariasis is a globally neglected tropical parasitic disease which affects individuals of all ages and leads to an altered lymphatic system and abnormal enlargement of body parts. METHODS A mathematical model of lymphatic filariaris with intervention strategies is developed and analyzed. Control of infections is analyzed within the model through medical treatment of infected-acute individuals and quarantine of infected-chronic individuals. RESULTS We derive the effective reproduction number, [Formula: see text] and its interpretation/investigation suggests that treatment contributes to a reduction in lymphatic filariasis cases faster than quarantine. However, this reduction is greater when the two intervention approaches are applied concurrently. CONCLUSIONS Numerical simulations are carried out to monitor the dynamics of the filariasis model sub-populations for various parameter values of the associated reproduction threshold. Lastly, sensitivity analysis on key parameters that drive the disease dynamics is performed in order to identify their relative importance on the disease transmission.
Collapse
Affiliation(s)
- Peter M Mwamtobe
- School of Computer Science and Applied Mathematics, University of the Witwatersrand, Johannesburg, Private Bag 3, Wits, Johannesburg, 2050, South Africa. .,DST-NRF Centre of Excellence in Mathematical and Statistical Sciences (CoE-MaSS), University of the Witwatersrand, Johannesburg, Private Bag 3, WitsJohannesburg, 2050, South Africa. .,Department of Mathematics and Statistics, University of Malawi, Chichiri, Blantyre, Malawi.
| | - Simphiwe M Simelane
- School of Computer Science and Applied Mathematics, University of the Witwatersrand, Johannesburg, Private Bag 3, Wits, Johannesburg, 2050, South Africa.,DST-NRF Centre of Excellence in Mathematical and Statistical Sciences (CoE-MaSS), University of the Witwatersrand, Johannesburg, Private Bag 3, WitsJohannesburg, 2050, South Africa
| | - Shirley Abelman
- School of Computer Science and Applied Mathematics, University of the Witwatersrand, Johannesburg, Private Bag 3, Wits, Johannesburg, 2050, South Africa.,DST-NRF Centre of Excellence in Mathematical and Statistical Sciences (CoE-MaSS), University of the Witwatersrand, Johannesburg, Private Bag 3, WitsJohannesburg, 2050, South Africa
| | - Jean M Tchuenche
- School of Computer Science and Applied Mathematics, University of the Witwatersrand, Johannesburg, Private Bag 3, Wits, Johannesburg, 2050, South Africa
| |
Collapse
|
15
|
Sevá AP, Ovallos FG, Amaku M, Carrillo E, Moreno J, Galati EAB, Lopes EG, Soares RM, Ferreira F. Canine-Based Strategies for Prevention and Control of Visceral Leishmaniasis in Brazil. PLoS One 2016; 11:e0160058. [PMID: 27471852 PMCID: PMC4966914 DOI: 10.1371/journal.pone.0160058] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2016] [Accepted: 07/13/2016] [Indexed: 11/19/2022] Open
Abstract
Visceral leishmaniasis (VL) is a zoonosis found worldwide. Its incidence has increased in Brazil in recent years, representing a serious public and animal health problem. The strategies applied in Brazil are questionable and are not sufficient to control the disease. Thus, we have compared the efficacy of some of the currently available strategies focused on dogs to prevent and control zoonotic VL in endemic areas by optimizing a mathematical model. The simulations showed that the elimination of seropositive dogs, the use of insecticide-impregnated dog collars, and the vaccination of dogs significantly contribute to reducing the prevalence of infection in both canines and humans. The use of insecticide-impregnated collars presented the highest level of efficacy mainly because it directly affected the force of infection and vector-dog contact. In addition, when used at a coverage rate of 90%, insecticide-impregnated collar was able to decrease the prevalence of seropositive dogs and humans to zero; moreover, because of the easy application and acceptance by the targeted population, these collars may be considered the most feasible for inclusion in public policies among the three simulated measures. Vaccination and euthanasia were efficacious, but the latter method is strongly criticized on ethical grounds, and both methods present difficulties for inclusion in public policies. When we compared the use of euthanasia and vaccination at coverages of 70 and 90%, respectively, the proportion of infected populations were similar. However, on evaluating the implications of both of these methods, particularly the negative aspects of culling dogs and the proportion of animals protected by vaccination, the latter measure appears to be the better option if the total cost is not significantly higher. The comparison of complications and advantages of different control strategies allows us to analyze the optimal measure and offer strategies to veterinary and public health authorities for making decisions to prevent and control zoonotic VL. Hence, improvements in both public and animal health can be achieved in regions with scenarios similar to that considered in the present study; such scenarios are characteristically found in some areas of Brazil and other countries.
Collapse
Affiliation(s)
- Anaiá P. Sevá
- Department of Preventive Veterinary Medicine and Animal Health, School of Veterinary Medicine and Animal Science, University of São Paulo, São Paulo, Brazil
- * E-mail:
| | - Fredy G. Ovallos
- Department of Epidemiology, School of Public Health, University of São Paulo, São Paulo, Brazil
| | - Marcus Amaku
- Department of Preventive Veterinary Medicine and Animal Health, School of Veterinary Medicine and Animal Science, University of São Paulo, São Paulo, Brazil
- Department of Pathology, School of Medicine, University of São Paulo, São Paulo, Brazil
| | - Eugenia Carrillo
- WHO Collaborating Centre for Leishmaniasis, Centro Nacional de Microbiología, Instituto de Salud Carlos III, Majadahonda, Spain
| | - Javier Moreno
- WHO Collaborating Centre for Leishmaniasis, Centro Nacional de Microbiología, Instituto de Salud Carlos III, Majadahonda, Spain
| | - Eunice A. B. Galati
- Department of Epidemiology, School of Public Health, University of São Paulo, São Paulo, Brazil
| | - Estela G. Lopes
- Department of Preventive Veterinary Medicine and Animal Health, School of Veterinary Medicine and Animal Science, University of São Paulo, São Paulo, Brazil
| | - Rodrigo M. Soares
- Department of Preventive Veterinary Medicine and Animal Health, School of Veterinary Medicine and Animal Science, University of São Paulo, São Paulo, Brazil
| | - Fernando Ferreira
- Department of Preventive Veterinary Medicine and Animal Health, School of Veterinary Medicine and Animal Science, University of São Paulo, São Paulo, Brazil
| |
Collapse
|
16
|
Brugger K, Köfer J, Rubel F. Outdoor and indoor monitoring of livestock-associated Culicoides spp. to assess vector-free periods and disease risks. BMC Vet Res 2016; 12:88. [PMID: 27259473 PMCID: PMC4893216 DOI: 10.1186/s12917-016-0710-z] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2016] [Accepted: 05/24/2016] [Indexed: 11/11/2022] Open
Abstract
Background Within the last few decades Culicoides spp. (Diptera: Ceratopogonidae) emerged Europe-wide as a major vector for epizootic viral diseases e.g. caused by Bluetongue (BT) or Schmallenberg virus. In accordance with the EU regulation 1266/2007, veterinary authorities are requested to determine vector-free periods for loosing trade and movement restrictions of susceptible livestock. Additionally, the widely used basic reproduction number \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$\mathcal {R}_{0}$\end{document}R0 is optionally applied for risk assessment of vector-borne diseases. Values of \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$\mathcal {R}_{0}<1$\end{document}R0<1 indicate periods with no disease transmission risk. For the determination of vector-free period and \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$\mathcal {R}_{0}$\end{document}R0 a continuously operating daily Culicoides spp. monitoring in Vienna (Austria) was established. It covered the period 2009–2013 and depicts the seasonal vector abundance indoor and outdoor. Future BT and African horse sickness (AHS) outbreak risks were estimated by projecting \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$\mathcal {R}_{0}$\end{document}R0 to climate change scenarios. Therefore, temperature-dependent vector parameters were applied. Results The vector-free period lasted about 100 days inside stables, while less than five Culicoides were trapped outdoors on 150 days per season, i.e. winter half year. Additionally, the potential outbreak risk was assessed for BT and AHS. For BT, a basic reproduction number of \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$\mathcal {R}_{0}>1$\end{document}R0>1 was found each year between June and August. The periods without transmission risk, i.e. \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$\mathcal {R}_{0}<1$\end{document}R0<1, were notably higher (200 days). Contrary, values of \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$\mathcal {R}_{0}<1$\end{document}R0<1 were estimated for AHS during the whole period. Finally, the basic reproduction numbers were projected to the future by using temperature forecasts for the period 2014–2100. While the mean summer peak values for BT increase from of \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$\mathcal {R}_{0}=2.3$\end{document}R0=2.3 to \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$\mathcal {R}_{0}=3.4$\end{document}R0=3.4 until 2100 (1.1/100 years), no risk for AHS was estimated even under climate warming assumptions. Conclusions Restrictions to trade and movement are always associated with an economic impact during epidemic diseases. To minimize these impacts, risk assessments based on the vector-free period or the basic reproduction number \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$\mathcal {R}_{0}$\end{document}R0 can essentially support veterinary authorities to improve protection and control measurements.
Collapse
Affiliation(s)
- Katharina Brugger
- Institute for Veterinary Public Health, University of Veterinary Medicine Vienna, Veterinaerplatz 1, Vienna, 1210, Austria.
| | - Josef Köfer
- Institute for Veterinary Public Health, University of Veterinary Medicine Vienna, Veterinaerplatz 1, Vienna, 1210, Austria
| | - Franz Rubel
- Institute for Veterinary Public Health, University of Veterinary Medicine Vienna, Veterinaerplatz 1, Vienna, 1210, Austria
| |
Collapse
|
17
|
Maran N, Gomes PS, Freire-de-Lima L, Freitas EO, Freire-de-Lima CG, Morrot A. Host resistance to visceral leishmaniasis: prevalence and prevention. Expert Rev Anti Infect Ther 2016; 14:435-42. [DOI: 10.1586/14787210.2016.1160779] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
|
18
|
Gorahava KK, Rosenberger JM, Mubayi A. Optimizing insecticide allocation strategies based on houses and livestock shelters for visceral leishmaniasis control in Bihar, India. Am J Trop Med Hyg 2015; 93:114-22. [PMID: 25940194 DOI: 10.4269/ajtmh.14-0612] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2014] [Accepted: 03/19/2015] [Indexed: 11/07/2022] Open
Abstract
Visceral leishmaniasis (VL) is the most deadly form of the leishmaniasis family of diseases, which affects numerous developing countries. The Indian state of Bihar has the highest prevalence and mortality rate of VL in the world. Insecticide spraying is believed to be an effective vector control program for controlling the spread of VL in Bihar; however, it is expensive and less effective if not implemented systematically. This study develops and analyzes a novel optimization model for VL control in Bihar that identifies an optimal (best possible) allocation of chosen insecticide (dichlorodiphenyltrichloroethane [DDT] or deltamethrin) based on the sizes of human and cattle populations in the region. The model maximizes the insecticide-induced sandfly death rate in human and cattle dwellings while staying within the current state budget for VL vector control efforts. The model results suggest that deltamethrin might not be a good replacement for DDT because the insecticide-induced sandfly deaths are 3.72 times more in case of DDT even after 90 days post spray. Different insecticide allocation strategies between the two types of sites (houses and cattle sheds) are suggested based on the state VL-control budget and have a direct implication on VL elimination efforts in a resource-limited region.
Collapse
Affiliation(s)
- Kaushik K Gorahava
- Industrial and Manufacturing Systems Engineering Department, The University of Texas at Arlington, Arlington, Texas; Department of Mathematics, Northeastern Illinois University, Chicago, Illinois; Department of Mathematics, The University of Texas at Arlington, Arlington, Texas; SAL Mathematical Computational and Modeling Science Center, Tempe, Arizona; School of Mathematical and Natural Sciences, Arizona State University, Phoenix, Arizona
| | - Jay M Rosenberger
- Industrial and Manufacturing Systems Engineering Department, The University of Texas at Arlington, Arlington, Texas; Department of Mathematics, Northeastern Illinois University, Chicago, Illinois; Department of Mathematics, The University of Texas at Arlington, Arlington, Texas; SAL Mathematical Computational and Modeling Science Center, Tempe, Arizona; School of Mathematical and Natural Sciences, Arizona State University, Phoenix, Arizona
| | - Anuj Mubayi
- Industrial and Manufacturing Systems Engineering Department, The University of Texas at Arlington, Arlington, Texas; Department of Mathematics, Northeastern Illinois University, Chicago, Illinois; Department of Mathematics, The University of Texas at Arlington, Arlington, Texas; SAL Mathematical Computational and Modeling Science Center, Tempe, Arizona; School of Mathematical and Natural Sciences, Arizona State University, Phoenix, Arizona
| |
Collapse
|
19
|
Nakagawa J, Ehrenberg JP, Nealon J, Fürst T, Aratchige P, Gonzales G, Chanthavisouk C, Hernandez LM, Fengthong T, Utzinger J, Steinmann P. Towards effective prevention and control of helminth neglected tropical diseases in the Western Pacific Region through multi-disease and multi-sectoral interventions. Acta Trop 2015; 141:407-18. [PMID: 23792012 DOI: 10.1016/j.actatropica.2013.05.010] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2012] [Revised: 04/22/2013] [Accepted: 05/21/2013] [Indexed: 12/13/2022]
Abstract
Neglected tropical diseases (NTDs) cause serious health, social and economic burdens in the countries of the World Health Organization Western Pacific Region. Among the NTDs, helminth infections are particularly prominent with regard to the number of infected individuals and health impact. Co-endemicity is common among impoverished and marginalized populations. To achieve effective and sustainable control of helminth NTDs, a deeper understanding of the social-ecological systems governing their endemicity and strategies beyond preventive chemotherapy are required to tackle the multiple causes of infection and re-infection. We discuss the feasibility of implementing multi-disease, multi-sectoral intervention packages for helminth NTDs in the Western Pacific Region. After reviewing the main determinants for helminth NTD endemicity and current control strategies, key control activities that involve or concern other programmes within and beyond the health sector are discussed. A considerable number of activities that have an impact on more than one helminth NTD are identified in a variety of sectors, suggesting an untapped potential for synergies. We also highlight the challenges of multi-sectoral collaboration, particularly of involving non-health sectors. We conclude that multi-sectoral collaboration for helminth NTD control is feasible if the target diseases and sectors are carefully selected. To do so, an incentive analysis covering key stakeholders in the sectors is crucial, and the disease-control strategies need to be well understood. The benefits of multi-disease, multi-sectoral approaches could go beyond immediate health impacts by contributing to sustainable development, raising educational attainment, increasing productivity and reducing health inequities.
Collapse
Affiliation(s)
- Jun Nakagawa
- World Health Organization, Regional Office for the Western Pacific, Division of Combating Communicable Diseases, P.O. Box 2932, 1000 Manila, Philippines
| | - John P Ehrenberg
- World Health Organization, Regional Office for the Western Pacific, Division of Combating Communicable Diseases, P.O. Box 2932, 1000 Manila, Philippines
| | - Joshua Nealon
- World Health Organization, Regional Office for the Western Pacific, Division of Combating Communicable Diseases, P.O. Box 2932, 1000 Manila, Philippines
| | - Thomas Fürst
- Department of Epidemiology and Public Health, Swiss Tropical and Public Health Institute, P.O. Box, CH-4002 Basel, Switzerland; University of Basel, P.O. Box, CH-4003 Basel, Switzerland
| | - Padmasiri Aratchige
- World Health Organization, Regional Office for the Western Pacific, Division of Combating Communicable Diseases, P.O. Box 2932, 1000 Manila, Philippines
| | - Glenda Gonzales
- World Health Organization, Regional Office for the Western Pacific, Division of Combating Communicable Diseases, P.O. Box 2932, 1000 Manila, Philippines
| | - Chitsavang Chanthavisouk
- World Health Organization, Regional Office for the Western Pacific, Division of Combating Communicable Diseases, P.O. Box 2932, 1000 Manila, Philippines
| | - Leda M Hernandez
- Infectious Disease Office, National Centre for Disease Prevention and Control, Department of Health, Sta. Cruz, 1000 Manila, Philippines
| | - Tayphasavanh Fengthong
- Department of Hygiene and Health Promotion, Ministry of Health, P.O. Box 1232, Vientiane, Lao People's Democratic Republic
| | - Jürg Utzinger
- Department of Epidemiology and Public Health, Swiss Tropical and Public Health Institute, P.O. Box, CH-4002 Basel, Switzerland; University of Basel, P.O. Box, CH-4003 Basel, Switzerland
| | - Peter Steinmann
- Department of Epidemiology and Public Health, Swiss Tropical and Public Health Institute, P.O. Box, CH-4002 Basel, Switzerland; University of Basel, P.O. Box, CH-4003 Basel, Switzerland.
| |
Collapse
|
20
|
Rebaudo F, Costa J, Almeida CE, Silvain JF, Harry M, Dangles O. Simulating population genetics of pathogen vectors in changing landscapes: guidelines and application with Triatoma brasiliensis. PLoS Negl Trop Dis 2014; 8:e3068. [PMID: 25102068 PMCID: PMC4125301 DOI: 10.1371/journal.pntd.0003068] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2014] [Accepted: 06/21/2014] [Indexed: 11/28/2022] Open
Abstract
Background Understanding the mechanisms that influence the population dynamics and spatial genetic structure of the vectors of pathogens infecting humans is a central issue in tropical epidemiology. In view of the rapid changes in the features of landscape pathogen vectors live in, this issue requires new methods that consider both natural and human systems and their interactions. In this context, individual-based model (IBM) simulations represent powerful yet poorly developed approaches to explore the response of pathogen vectors in heterogeneous social-ecological systems, especially when field experiments cannot be performed. Methodology/Principal Findings We first present guidelines for the use of a spatially explicit IBM, to simulate population genetics of pathogen vectors in changing landscapes. We then applied our model with Triatoma brasiliensis, originally restricted to sylvatic habitats and now found in peridomestic and domestic habitats, posing as the most important Trypanosoma cruzi vector in Northeastern Brazil. We focused on the effects of vector migration rate, maximum dispersal distance and attraction by domestic habitat on T. brasiliensis population dynamics and spatial genetic structure. Optimized for T. brasiliensis using field data pairwise fixation index (FST) from microsatellite loci, our simulations confirmed the importance of these three variables to understand vector genetic structure at the landscape level. We then ran prospective scenarios accounting for land-use change (deforestation and urbanization), which revealed that human-induced land-use change favored higher genetic diversity among sampling points. Conclusions/Significance Our work shows that mechanistic models may be useful tools to link observed patterns with processes involved in the population genetics of tropical pathogen vectors in heterogeneous social-ecological landscapes. Our hope is that our study may provide a testable and applicable modeling framework to a broad community of epidemiologists for formulating scenarios of landscape change consequences on vector dynamics, with potential implications for their surveillance and control. Worldwide, humans are modifying landscapes at an unprecedented rate. These modifications have an influence on the ecology of pathogen vectors, yet this issue has received relatively little input from modeling research. The current study presents guidelines for the use of a modeling framework for the representation of the dynamics and spatial genetic structure of pathogen vectors. It allows considering spatiotemporal landscape modifications explicitly, to represent human-altered modifications and consequences. We applied this modeling framework to Triatoma brasiliensis, vector of the pathogen Trypanosoma cruzi responsible for the Chagas disease, in the semi-arid Northeastern Brazil. Using field data of pairwise fixation index (FST) from microsatellite loci, we found that migration rate, maximum dispersal distance and attraction by domestic habitat were all key parameters to understand vector spatial genetic structure at the landscape level. At the interface across disciplines, this study provides to the community of epidemiologists a testable and applicable framework to foresee landscape modification consequences on vector dynamics and genetic structure, with potential implications for their surveillance and control.
Collapse
Affiliation(s)
- Francois Rebaudo
- BEI-UR072, IRD, Gif-sur-Yvette, France
- LEGS-UPR9034, CNRS-UPSud11, Gif-sur-Yvette, France
- * E-mail:
| | - Jane Costa
- Laboratório de Biodiversidade Entomológica, Instituto Oswaldo Cruz - Fiocruz, Rio de Janeiro, Rio de Janeiro, Brasil
| | - Carlos E. Almeida
- Departamento de Ciências Biológicas, Faculdade de Ciências Farmacêuticas, UNESP, Araraquara, Sao Paolo, Brasil
| | - Jean-Francois Silvain
- BEI-UR072, IRD, Gif-sur-Yvette, France
- LEGS-UPR9034, CNRS-UPSud11, Gif-sur-Yvette, France
| | - Myriam Harry
- LEGS-UPR9034, CNRS-UPSud11, Gif-sur-Yvette, France
| | - Olivier Dangles
- BEI-UR072, IRD, Gif-sur-Yvette, France
- LEGS-UPR9034, CNRS-UPSud11, Gif-sur-Yvette, France
- Instituto de Ecología, Campus Cotacota, Universidad Mayor San Andrés, La Paz, Bolivia
| |
Collapse
|
21
|
Xue L, Cohnstaedt LW, Scott HM, Scoglio C. A hierarchical network approach for modeling Rift Valley fever epidemics with applications in North America. PLoS One 2013; 8:e62049. [PMID: 23667453 PMCID: PMC3646918 DOI: 10.1371/journal.pone.0062049] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2012] [Accepted: 03/18/2013] [Indexed: 11/29/2022] Open
Abstract
Rift Valley fever is a vector-borne zoonotic disease which causes high morbidity and mortality in livestock. In the event Rift Valley fever virus is introduced to the United States or other non-endemic areas, understanding the potential patterns of spread and the areas at risk based on disease vectors and hosts will be vital for developing mitigation strategies. Presented here is a general network-based mathematical model of Rift Valley fever. Given a lack of empirical data on disease vector species and their vector competence, this discrete time epidemic model uses stochastic parameters following several PERT distributions to model the dynamic interactions between hosts and likely North American mosquito vectors in dispersed geographic areas. Spatial effects and climate factors are also addressed in the model. The model is applied to a large directed asymmetric network of 3,621 nodes based on actual farms to examine a hypothetical introduction to some counties of Texas, an important ranching area in the United States of America. The nodes of the networks represent livestock farms, livestock markets, and feedlots, and the links represent cattle movements and mosquito diffusion between different nodes. Cattle and mosquito (Aedes and Culex) populations are treated with different contact networks to assess virus propagation. Rift Valley fever virus spread is assessed under various initial infection conditions (infected mosquito eggs, adults or cattle). A surprising trend is fewer initial infectious organisms result in a longer delay before a larger and more prolonged outbreak. The delay is likely caused by a lack of herd immunity while the infection expands geographically before becoming an epidemic involving many dispersed farms and animals almost simultaneously. Cattle movement between farms is a large driver of virus expansion, thus quarantines can be efficient mitigation strategy to prevent further geographic spread.
Collapse
Affiliation(s)
- Ling Xue
- Kansas State Epicenter, Department of Electrical and Computer Engineering, Kansas State University, Manhattan, Kansas, United States of America
| | - Lee W. Cohnstaedt
- Center for Grain and Animal Health Research, United States Department of Agriculture, Manhattan, Kansas, United States of America
- * E-mail:
| | - H. Morgan Scott
- Department of Diagnostic Medicine/Pathobiology, Kansas State University, Manhattan, Kansas, United States of America
| | - Caterina Scoglio
- Kansas State Epicenter, Department of Electrical and Computer Engineering, Kansas State University, Manhattan, Kansas, United States of America
| |
Collapse
|
22
|
Multiscale analysis for a vector-borne epidemic model. J Math Biol 2013; 68:1269-93. [DOI: 10.1007/s00285-013-0666-6] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2011] [Revised: 02/14/2013] [Indexed: 11/30/2022]
|
23
|
Montresor A, Gabrielli AF, Yajima A, Lethanh N, Biggs BA, Casey GJ, Tinh TT, Engels D, Savioli L. Markov model to forecast the change in prevalence of soil-transmitted helminths during a control programme: a case study in Vietnam. Trans R Soc Trop Med Hyg 2013; 107:313-8. [PMID: 23471919 DOI: 10.1093/trstmh/trt019] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
BACKGROUND A mathematical model based on the Markov methodology to predict the change in prevalence of soil-transmitted helminth (STH) infections during public health control activities is not available, but would be an extremely efficient planning tool. METHOD We used the parasitological data collected during a deworming and iron supplementation programme for women of child-bearing age conducted in Vietnam between 2006 and 2011 to develop a Markov transition probability model. The transition probabilities were calculated from the observed changes in prevalence in the different classes of intensity for each STH species during the first year of intervention. The model was then developed and used to estimate the prevalence in year 2, 3, 4 and 5 for each STH species and for 'any STH infection'. The prevalence predicted by the model was then compared with the prevalence observed at different times during programme implementation. RESULTS The comparison between the model-predicted prevalence and the observed prevalence proved a good fit of the model. CONCLUSIONS We consider the Markov transition probability model to be a promising method of predicting changes in STH prevalence during control efforts. Further research to validate the model with observed data in different geographical and epidemiological settings is suggested to refine the prediction model.
Collapse
Affiliation(s)
- Antonio Montresor
- Department of Control of Neglected Tropical Diseases, World Heath Organization, Avenue Appia 20, 1211 Geneva, Switzerland.
| | | | | | | | | | | | | | | | | |
Collapse
|
24
|
K. Misra A, Sharma A, Li J. A mathematical model for control of vector borne diseases
through media campaigns. ACTA ACUST UNITED AC 2013. [DOI: 10.3934/dcdsb.2013.18.1909] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
25
|
Reck J, Marks FS, Guimarães JA, Termignoni C, Martins JR. Epidemiology of Ornithodoros brasiliensis (mouro tick) in the southern Brazilian highlands and the description of human and animal retrospective cases of tick parasitism. Ticks Tick Borne Dis 2012; 4:101-9. [PMID: 23238249 DOI: 10.1016/j.ttbdis.2012.09.004] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2012] [Revised: 09/21/2012] [Accepted: 09/23/2012] [Indexed: 10/27/2022]
Abstract
Ornithodoros brasiliensis, also known as the "mouro" tick, is an argasid tick found exclusively in the southern Brazilian highlands. O. brasiliensis parasitism is frequently associated with severe symptoms directly induced by the tick bite, a condition compatible with the definition of tick toxicosis. The objectives of this work include (i) the determination of the distribution of O. brasiliensis in farms located in the tick-endemic region, (ii) the description of the characteristics of O. brasiliensis habitats, (iii) the analysis of risk factors associated with O. brasiliensis, and (iv) the retrospective description of cases of human and animal parasitism by O. brasiliensis. Of the 30 farms included in this study, O. brasiliensis was identified on 5 farms (frequency 16.7%), in which several ticks found in high density buried in soil were collected. Information regarding the tick habitats and the local population was recorded. The data indicated that O. brasiliensis feeds on humans, dogs, armadillos (Dasypus hybridus), and possibly skunks (Conepatus chinga). The analysis of risk factors indicated that the presence of house basements with an unpaved (natural soil) floor on farms and insufficient sanitary conditions significantly enhanced the probability of identifying O. brasiliensis. Additionally, we describe retrospectively cases of tick parasitism in 28 humans and 11 dogs including the most common symptoms associated with tick toxicosis. This is the first study concerning O. brasiliensis epidemiology, distribution, and habitat, and the report represents the most comprehensive characterization of Ornithodoros bite-associated toxicosis syndrome.
Collapse
Affiliation(s)
- José Reck
- Instituto de Pesquisas Veterinárias Desidério Finamor, Fundação Estadual de Pesquisa Agropecuária, Eldorado do Sul, RS, Brazil.
| | | | | | | | | |
Collapse
|
26
|
An age-structured extension to the vectorial capacity model. PLoS One 2012; 7:e39479. [PMID: 22724022 PMCID: PMC3378582 DOI: 10.1371/journal.pone.0039479] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2012] [Accepted: 05/23/2012] [Indexed: 11/20/2022] Open
Abstract
Background Vectorial capacity and the basic reproductive number (R0) have been instrumental in structuring thinking about vector-borne pathogen transmission and how best to prevent the diseases they cause. One of the more important simplifying assumptions of these models is age-independent vector mortality. A growing body of evidence indicates that insect vectors exhibit age-dependent mortality, which can have strong and varied affects on pathogen transmission dynamics and strategies for disease prevention. Methodology/Principal Findings Based on survival analysis we derived new equations for vectorial capacity and R0 that are valid for any pattern of age-dependent (or age–independent) vector mortality and explore the behavior of the models across various mortality patterns. The framework we present (1) lays the groundwork for an extension and refinement of the vectorial capacity paradigm by introducing an age-structured extension to the model, (2) encourages further research on the actuarial dynamics of vectors in particular and the relationship of vector mortality to pathogen transmission in general, and (3) provides a detailed quantitative basis for understanding the relative impact of reductions in vector longevity compared to other vector-borne disease prevention strategies. Conclusions/Significance Accounting for age-dependent vector mortality in estimates of vectorial capacity and R0 was most important when (1) vector densities are relatively low and the pattern of mortality can determine whether pathogen transmission will persist; i.e., determines whether R0 is above or below 1, (2) vector population growth rate is relatively low and there are complex interactions between birth and death that differ fundamentally from birth-death relationships with age-independent mortality, and (3) the vector exhibits complex patterns of age-dependent mortality and R0∼1. A limiting factor in the construction and evaluation of new age-dependent mortality models is the paucity of data characterizing vector mortality patterns, particularly for free ranging vectors in the field.
Collapse
|
27
|
Auty HK, Picozzi K, Malele I, Torr SJ, Cleaveland S, Welburn S. Using molecular data for epidemiological inference: assessing the prevalence of Trypanosoma brucei rhodesiense in tsetse in Serengeti, Tanzania. PLoS Negl Trop Dis 2012; 6:e1501. [PMID: 22303496 PMCID: PMC3269424 DOI: 10.1371/journal.pntd.0001501] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2011] [Accepted: 12/12/2011] [Indexed: 11/19/2022] Open
Abstract
Background Measuring the prevalence of transmissible Trypanosoma brucei rhodesiense in tsetse populations is essential for understanding transmission dynamics, assessing human disease risk and monitoring spatio-temporal trends and the impact of control interventions. Although an important epidemiological variable, identifying flies which carry transmissible infections is difficult, with challenges including low prevalence, presence of other trypanosome species in the same fly, and concurrent detection of immature non-transmissible infections. Diagnostic tests to measure the prevalence of T. b. rhodesiense in tsetse are applied and interpreted inconsistently, and discrepancies between studies suggest this value is not consistently estimated even to within an order of magnitude. Methodology/Principal Findings Three approaches were used to estimate the prevalence of transmissible Trypanosoma brucei s.l. and T. b. rhodesiense in Glossina swynnertoni and G. pallidipes in Serengeti National Park, Tanzania: (i) dissection/microscopy; (ii) PCR on infected tsetse midguts; and (iii) inference from a mathematical model. Using dissection/microscopy the prevalence of transmissible T. brucei s.l. was 0% (95% CI 0–0.085) for G. swynnertoni and 0% (0–0.18) G. pallidipes; using PCR the prevalence of transmissible T. b. rhodesiense was 0.010% (0–0.054) and 0.0089% (0–0.059) respectively, and by model inference 0.0064% and 0.00085% respectively. Conclusions/Significance The zero prevalence result by dissection/microscopy (likely really greater than zero given the results of other approaches) is not unusual by this technique, often ascribed to poor sensitivity. The application of additional techniques confirmed the very low prevalence of T. brucei suggesting the zero prevalence result was attributable to insufficient sample size (despite examination of 6000 tsetse). Given the prohibitively high sample sizes required to obtain meaningful results by dissection/microscopy, PCR-based approaches offer the current best option for assessing trypanosome prevalence in tsetse but inconsistencies in relating PCR results to transmissibility highlight the need for a consensus approach to generate meaningful and comparable data. Human African trypanosomiasis is a fatal disease that is carried by a tsetse vector. Assessing the proportion of tsetse which carries human-infective trypanosomes is important in assessing human disease risk and understanding disease transmission dynamics. However, identifying flies which carry transmissible infections is difficult, due to potential presence of other trypanosome species in the same fly, and concurrent detection of immature infections which are not transmissible. We used three methods to estimate the proportion of flies carrying human-infective trypanosomes: dissection and microscopic examination of flies to visualise trypanosomes directly in the fly; PCR of fly midguts in which trypanosomes were observed by microscopy; and theoretical analysis using a mathematical model of disease transmission. All three methods found the prevalence to be extremely low. Given the low prevalence, dissection/microscopy requires prohibitively large sample sizes and therefore PCR-based approaches are likely to be of most value. However, interpretation of PCR data is not straightforward; whilst PCR identifies flies carrying pathogen genetic material it does not directly identify flies with transmissible infections. This study highlights the need for a consensus approach on the analysis and interpretation of PCR data to generate reliable and comparable measures of the proportion of flies which carry transmissible human-infective trypanosomes.
Collapse
Affiliation(s)
- Harriet K. Auty
- Division of Pathway Medicine and Centre for Infectious Diseases, School of Biomedical Sciences, College of Medicine and Veterinary Medicine, The University of Edinburgh, Edinburgh, United Kingdom
- Institute for Biodiversity, Animal Health and Comparative Medicine, College of Medicine, Veterinary Medicine and Life Sciences, University of Glasgow, Glasgow, United Kingdom
| | - Kim Picozzi
- Division of Pathway Medicine and Centre for Infectious Diseases, School of Biomedical Sciences, College of Medicine and Veterinary Medicine, The University of Edinburgh, Edinburgh, United Kingdom
| | - Imna Malele
- Tsetse and Trypanosomiasis Research Institute, Tanga, Tanzania
| | - Steve J. Torr
- Natural Resources Institute, University of Greenwich, Chatham Maritime, United Kingdom
| | - Sarah Cleaveland
- Institute for Biodiversity, Animal Health and Comparative Medicine, College of Medicine, Veterinary Medicine and Life Sciences, University of Glasgow, Glasgow, United Kingdom
| | - Sue Welburn
- Division of Pathway Medicine and Centre for Infectious Diseases, School of Biomedical Sciences, College of Medicine and Veterinary Medicine, The University of Edinburgh, Edinburgh, United Kingdom
- * E-mail:
| |
Collapse
|
28
|
Zheng W, Chen H, Liu X, Guo X, Fu R. Severe tick infestation in a hare and potential risk for transmitting pathogens to humans. THE KOREAN JOURNAL OF PARASITOLOGY 2011; 49:419-22. [PMID: 22355211 PMCID: PMC3279682 DOI: 10.3347/kjp.2011.49.4.419] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/31/2011] [Revised: 08/10/2011] [Accepted: 08/10/2011] [Indexed: 11/23/2022]
Abstract
Severe tick infestation was found in a hare in a suburban area of Nanchang, Jiangxi Province, China. We sampled ticks and identified them based on their morphologic characteristics. Three species, Ixodes sinensis, which is commonly found in China and can experimentally transmit Borrelia burgdorferi, Rhipicephalus haemaphysaloides, and Haemaphysalis longicornis which can transmit Lyme disease were detected with an optical microscope and a stereomicroscope. Risk of spreading ticks from suburban to urban areas exists due to human transportation and travel between the infested and non-infested areas around Nanchang.
Collapse
Affiliation(s)
- Weiqing Zheng
- Nangchang Center for Disease Control and Prevention, Nanchang 330038, China.
| | | | | | | | | |
Collapse
|