Background Mathematical modelling of host-parasite systems has seen tremendous developments and broad applications in theoretical and applied ecology. The current study focuses on the infection dynamics of a gyrodactylid-fish system. Previous experimental studies have explored the infrapopulation dynamics of co-infecting ectoparasites, Gyrodactylus turnbulli and G. bullatarudis, on their fish host, Poecilia reticulata, but questions remain about parasite microhabitat preferences, host survival and parasite virulence over time. Here, we use more advanced statistics and a sophisticated mathematical model to investigate these questions based on empirical data to add to our understanding of this gyrodactylid-fish system. Methods A rank-based multivariate Kruskal-Wallis test coupled with its post-hoc tests and graphical summaries were used to investigate the spatial and temporal parasite distribution of different gyrodactylid strains across different host populations. By adapting a multi-state Markov model that extends the standard survival models, we improved previous estimates of survival probabilities. Finally, we quantified parasite virulence of three different strains as a function of host mortality and recovery across different fish stocks and sexes. Results We confirmed that the captive-bred G. turnbulli and wild G. bullatarudis strains preferred the caudal and rostral regions respectively across different fish stocks; however, the wild G. turnbulli strain changed microhabitat preference over time, indicating microhabitat preference of gyrodactylids is host and time dependent. The average time of host infection before recovery or death was between 6 and 14 days. For this gyrodactylid-fish system, a longer period of host infection led to a higher chance of host recovery. Parasite-related mortalities are host, sex and time dependent, whereas fish size is confirmed to be the key determinant of host recovery. Conclusion From existing empirical data, we provided new insights into the gyrodactylid-fish system. This study could inform the modelling of other host-parasite interactions where the entire infection history of the host is of interest by adapting multi-state Markov models. Such models are under-utilised in parasitological studies and could be expanded to estimate relevant epidemiological traits concerning parasite virulence and host survival. Graphical Abstract

Several mathematical and standard epidemiological models have been proposed in studying infectious disease dynamics. These models help to understand the spread of disease infections. However, most of these models are not able to estimate other relevant disease metrics such as probability of first infection and recovery as well as the expected time to infection and recovery for both susceptible and infected individuals. That is, most of the standard epidemiological models used in estimating transition probabilities (TPs) are not able to generalize the transition estimates of disease outcomes at discrete time steps for future predictions. This paper seeks to address the aforementioned problems through a discrete-time Markov chain model. Secondary datasets from cohort studies were collected on HIV, tuberculosis (TB), and hepatitis B (HB) cases from a regional hospital in Ghana. The Markov chain model revealed that hepatitis B was more infectious over time than tuberculosis and HIV even though the probability of first infection of these diseases was relatively low within the study population. However, individuals infected with HIV had comparatively lower life expectancies than those infected with tuberculosis and hepatitis B. Discrete-time Markov chain technique is recommended as viable for modeling disease dynamics in Ghana.

No abstract

Mathematical models can aid in elucidating the spread of infectious disease dynamics within a given population over time. In an attempt to model tuberculosis (TB) dynamics among high-burden districts in the Ashanti Region of Ghana, the SEIR epidemic model with demography was employed within both deterministic and stochastic settings for comparison purposes. The deterministic model showed success in modelling TB infection in the region to the transmission dynamics of the stochastic SEIR model over time. It predicted tuberculosis dying out in ten of twelve high-burden districts in the Ashanti Region, but an outbreak in Obuasi municipal and Amansie West district. The effect of introducing treatment at the incubation stage of TB transmission was also investigated, and it was discovered that treatment introduced at the exposed stage decreased the spread of TB. Branching process approximation was used to derive explicit forms of relevant epidemiological quantities of the deterministic SEIR model for stability analysis of equilibrium points. Numerical simulations were performed to validate the overall infection rate, basic reproductive number, herd immunity threshold, and Malthusian parameter based on bootstrapping, jackknife, and Latin Hypercube sampling schemes. It was recommended that the Ghana Health Service should find a good mechanism to detect TB in the early stages of infection in the region. Public health attention must also be given to districts with a potentially higher risk of experiencing endemic TB even though the estimates of the overall epidemic thresholds from our SEIR model suggested that the Ashanti Region as a whole had herd immunity against TB infection.