Opportunities and challenges for modelling epidemiological and evolutionary dynamics in a multihost, multiparasite system: Zoonotic hybrid schistosomiasis in West Africa

Multihost multiparasite systems are evolutionarily and ecologically dynamic, which presents substantial trans‐disciplinary challenges for elucidating their epidemiology and designing appropriate control. Evidence for hybridizations and introgressions between parasite species is gathering, in part in line with improvements in molecular diagnostics and genome sequencing. One major system where this is becoming apparent is within the Genus Schistosoma, where schistosomiasis represents a disease of considerable medical and veterinary importance, the greatest burden of which occurs in sub‐Saharan Africa. Interspecific hybridizations and introgressions bring an increased level of complexity over and above that already inherent within multihost, multiparasite systems, also representing an additional source of genetic variation that can drive evolution. This has the potential for profound implications for the control of parasitic diseases, including, but not exclusive to, widening host range, increased transmission potential and altered responses to drug therapy. Here, we present the challenging case example of haematobium group Schistosoma spp. hybrids in West Africa, a system involving multiple interacting parasites and multiple definitive hosts, in a region where zoonotic reservoirs of schistosomiasis were not previously considered to be of importance. We consider how existing mathematical model frameworks for schistosome transmission could be expanded and adapted to zoonotic hybrid systems, exploring how such model frameworks can utilize molecular and epidemiological data, as well as the complexities and challenges this presents. We also highlight the opportunities and value such mathematical models could bring to this and a range of similar multihost, multi and cross‐hybridizing parasites systems in our changing world

Влияние фенилгидразина и алкилсульфатов на осмотическую чувствительность эритроцитов млекопитающих

Показано, что чувствительность эритроцитов млекопитающих к гипертоническому стрессу при 37 °С после обработки фенилгидразином зависит от видовой принадлежности клеток. Установлено, что модификация фенилгидразином эритроцитов млекопитающих приводит к снижению антигемолитической активности алкилкилсульфатов в условиях гипертонического стресса.Показано, що чутливість еритроцитів ссавців до гіпертонічного стресу при 37 °С після обробки фенілгідразином залежить від видової приналежності клітин. Встановлено, що модифікація фенілгідразином еритроцитів ссавців призводить до зниження антигемолітичної активності алкілсульфатів в умовах гіпертонічного стресу.It is shown that the sensitivity of mammalian erythrocytes to hypertonic stress at 37 °C after the treatment with phenylhydrazine depends on the species of a cell. We have established that the modification with phenylhydrazine of mammalian erythrocytes leads to a decrease of the alkyl sulfates antihemolytic activity under hypertonic stress

Genetic analysis of infectious bronchitis virus (IBV) in vaccinated poultry populations over a period of 10 years

Infectious bronchitis virus (IBV) is an avian pathogen from the Coronavirus family causing major health issues in poultry flocks worldwide. Because of its negative impact on health, performance, and bird welfare, commercial poultry are routinely vaccinated by administering live attenuated virus. However, field strains are capable of rapid adaptation and may evade vaccine-induced immunity. We set out to describe dynamics within and between lineages and assess potential escape from vaccine-induced immunity. We investigated a large nucleotide sequence database of over 1700 partial sequences of the S1 spike protein gene collected from clinical samples of Dutch chickens submitted to the laboratory of Royal GD between 2011 and 2020. Relative frequencies of the two major lineages GI-13 (793B) and GI-19 (QX) did not change in the investigated period, but we found a succession of distinct GI-19 sublineages. Analysis of dN/dS ratio over all sequences demonstrated episodic diversifying selection acting on multiple sites, some of which overlap predicted N-glycosylation motifs. We assessed several measures that would indicate divergence from vaccine strains, both in the overall database and in the two major lineages. However, the frequency of vaccine-homologous lineages did not decrease, no increase in genetic variation with time was detected, and the sequences did not grow more divergent from vaccine sequences in the examined time window. Concluding, our results show sublineage turnover within the GI-19 lineage and we demonstrate episodic diversifying selection acting on the partial sequence, but we cannot confirm nor rule out escape from vaccine-induced immunity. RESEARCH HIGHLIGHTS Succession of GI-19 IBV variants in broiler populations. IBV lineages overrepresented in either broiler, or layer production chickens. Ongoing episodic selection at the IBV S1 spike protein gene sequence. Several positively selected codons coincident with N-glycosylation motifs

Model-consistent estimation of the basic reproduction number from the incidence of an emerging infection

We investigate the merit of deriving an estimate of 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} early in an outbreak of an (emerging) infection from estimates of the incidence and generation interval only. We compare such estimates 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$$ \end{document} with estimates incorporating additional model assumptions, and determine the circumstances under which the different estimates are consistent. We show that one has to be careful when using observed exponential growth rates to derive an estimate 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$$ \end{document} , and we quantify the discrepancies that arise

Extinction times in the subcritical stochastic SIS logistic epidemic

Many real epidemics of an infectious disease are not straightforwardly super- or sub-critical, and the understanding of epidemic models that exhibit such complexity has been identified as a priority for theoretical work. We provide insights into the near-critical regime by considering the stochastic SIS logistic epidemic, a well-known birth-and-death chain used to model the spread of an epidemic within a population of a given size $N$. We study the behaviour of the process as the population size $N$ tends to infinity. Our results cover the entire subcritical regime, including the "barely subcritical" regime, where the recovery rate exceeds the infection rate by an amount that tends to 0 as $N \to \infty$ but more slowly than $N^{-1/2}$. We derive precise asymptotics for the distribution of the extinction time and the total number of cases throughout the subcritical regime, give a detailed description of the course of the epidemic, and compare to numerical results for a range of parameter values. We hypothesise that features of the course of the epidemic will be seen in a wide class of other epidemic models, and we use real data to provide some tentative and preliminary support for this theory.Comment: Revised; 34 pages; 6 figure

Prioritizing Emerging Zoonoses in The Netherlands

Background: To support the development of early warning and surveillance systems of emerging zoonoses, we present a general method to prioritize pathogens using a quantitative, stochastic multi-criteria model, parameterized for the Netherlands. Methodology/Principal Findings: A risk score was based on seven criteria, reflecting assessments of the epidemiology and impact of these pathogens on society. Criteria were weighed, based on the preferences of a panel of judges with a background in infectious disease control. Conclusions/Significance: Pathogens with the highest risk for the Netherlands included pathogens in the livestock reservoir with a high actual human disease burden (e.g. Campylobacter spp., Toxoplasma gondii, Coxiella burnetii) or a low current but higher historic burden (e.g. Mycobacterium bovis), rare zoonotic pathogens in domestic animals with severe disease manifestations in humans (e.g. BSE prion, Capnocytophaga canimorsus) as well as arthropod-borne and wildlife associated pathogens which may pose a severe risk in future (e.g. Japanese encephalitis virus and West-Nile virus). These agents are key targets for development of early warning and surveillance.Infrastructures, Systems and ServicesTechnology, Policy and Managemen

Early Epidemiological Assessment of the Virulence of Emerging Infectious Diseases: A Case Study of an Influenza Pandemic

Background: The case fatality ratio (CFR), the ratio of deaths from an infectious disease to the number of cases, provides an assessment of virulence. Calculation of the ratio of the cumulative number of deaths to cases during the course of an epidemic tends to result in a biased CFR. The present study develops a simple method to obtain an unbiased estimate of confirmed CFR (cCFR), using only the confirmed cases as the denominator, at an early stage of epidemic, even when there have been only a few deaths. Methodology/Principal Findings: Our method adjusts the biased cCFR by a factor of underestimation which is informed by the time from symptom onset to death. We first examine the approach by analyzing an outbreak of severe acute respiratory syndrome in Hong Kong (2003) with known unbiased cCFR estimate, and then investigate published epidemiological datasets of novel swine-origin influenza A (H1N1) virus infection in the USA and Canada (2009). Because observation of a few deaths alone does not permit estimating the distribution of the time from onset to death, the uncertainty is addressed by means of sensitivity analysis. The maximum likelihood estimate of the unbiased cCFR for influenza may lie in the range of 0.16-4.48% within the assumed parameter space for a factor of underestimation. The estimates for influenza suggest that the virulence is comparable to the early estimate in Mexico. Even when there have been no deaths, our model permits estimating a conservative upper bound of the cCFR. Conclusions: Although one has to keep in mind that the cCFR for an entire population is vulnerable to its variations among sub-populations and underdiagnosis, our method is useful for assessing virulence at the early stage of an epidemic and for informing policy makers and the public. © 2009 Nishiura et al.published_or_final_versio

Modeling infectious disease dynamics in the complex landscape of global health.

Despite some notable successes in the control of infectious diseases, transmissible pathogens still pose an enormous threat to human and animal health. The ecological and evolutionary dynamics of infections play out on a wide range of interconnected temporal, organizational, and spatial scales, which span hours to months, cells to ecosystems, and local to global spread. Moreover, some pathogens are directly transmitted between individuals of a single species, whereas others circulate among multiple hosts, need arthropod vectors, or can survive in environmental reservoirs. Many factors, including increasing antimicrobial resistance, increased human connectivity and changeable human behavior, elevate prevention and control from matters of national policy to international challenge. In the face of this complexity, mathematical models offer valuable tools for synthesizing information to understand epidemiological patterns, and for developing quantitative evidence for decision-making in global health

Effects of infection-induced migration delays on the epidemiology of avian influenza in wild mallard populations

Wild waterfowl populations form a natural reservoir of Avian Influenza (AI) virus, and fears exist that these birds may contribute to an AI pandemic by spreading the virus along their migratory flyways. Observational studies suggest that individuals infected with AI virus may delay departure from migratory staging sites. Here, we explore the epidemiological dynamics of avian influenza virus in a migrating mallard (Anas platyrhynchos) population with a specific view to understanding the role of infection-induced migration delays on the spread of virus strains of differing transmissibility. We develop a host-pathogen model that combines the transmission dynamics of influenza with the migration, reproduction and mortality of the host bird species. Our modeling predicts that delayed migration of individuals influences both the timing and size of outbreaks of AI virus. We find that (1) delayed migration leads to a lower total number of cases of infection each year than in the absence of migration delay, (2) when the transmission rate of a strain is high, the outbreak starts at the staging sites at which birds arrive in the early part of the fall migration, (3) when the transmission rate is low, infection predominantly occurs later in the season, which is further delayed when there is a migration delay. As such, the rise of more virulent AI strains in waterfowl could lead to a higher prevalence of infection later in the year, which could change the exposure risk for farmed poultry. A sensitivity analysis shows the importance of generation time and loss of immunity for the effect of migration delays. Thus, we demonstrate, in contrast to many current transmission risk models solely using empirical information on bird movements to assess the potential for transmission, that a consideration of infection-induced delays is critical to understanding the dynamics of AI infection along the entire flyway.<br /