Application of a computational model for complex fluvial ecosystems: The population dynamics of zebra mussel Dreissena polymorpha as a case study

The potential and adaptive flexibility of population dynamic P-systems (PDP) to study population dynamics suggests that they may be suitable for modelling complex fluvial ecosystems, characterized by a composition of dynamic habitats with many variables that interact simultaneously. Using as a model a reservoir occupied by the zebra mussel Dreissena polymorpha, we designed a computational model based on P systems to study the population dynamics of larvae, in order to evaluate management actions to control or eradicate this invasive species. The population dynamics of this species was simulated under different scenarios ranging from the absence of water flow change to a weekly variation with different flow rates, to the actual hydrodynamic situation of an intermediate flow rate. Our results show that PDP models can be very useful tools to model complex, partially desynchronized, processes that work in parallel. This allows the study of complex hydroecological processes such as the one presented, where reproductive cycles, temperature and water dynamics are involved in the desynchronization of the population dynamics both, within areas and among them. The results obtained may be useful in the management of other reservoirs with similar hydrodynamic situations in which the presence of this invasive species has been documented.

Spatial and temporal variation in population genetic structure of wild Nile tilapia (Oreochromis niloticus) across Africa

Background: Reconstructing the evolutionary history of a species is challenging. It often depends not only on the past biogeographic and climatic events but also the contemporary and ecological factors, such as current connectivity and habitat heterogeneity. In fact, these factors might interact with each other and shape the current species distribution. However, to what extent the current population genetic structure reflects the past and the contemporary factors is largely unknown. Here we investigated spatio-temporal genetic structures of Nile tilapia (Oreochromis niloticus) populations, across their natural distribution in Africa. While its large biogeographic distribution can cause genetic differentiation at the paleo-biogeographic scales, its restricted dispersal capacity might induce a strong genetic structure at micro-geographic scales. Results: Using nine microsatellite loci and 350 samples from ten natural populations, we found the highest genetic differentiation among the three ichthyofaunal provinces and regions (Ethiopian, Nilotic and Sudano-Sahelian) (R(ST) = 0.38 - 0.69). This result suggests the predominant effect of paleo-geographic events at macro-geographic scale. In addition, intermediate divergences were found between rivers and lakes within the regions, presumably reflecting relatively recent interruptions of gene flow between hydrographic basins (R(ST) = 0.24 - 0.32). The lowest differentiations were observed among connected populations within a basin (R(ST) = 0.015 in the Volta basin). Comparison of temporal sample series revealed subtle changes in the gene pools in a few generations (F = 0 - 0.053). The estimated effective population sizes were 23 - 143 and the estimated migration rate was moderate (m similar to 0.094 - 0.097) in the Volta populations. Conclusions: This study revealed clear hierarchical patterns of the population genetic structuring of O. niloticus in Africa. The effects of paleo-geographic and climatic events were predominant at macro-geographic scale, and the significant effect of geographic connectivity was detected at micro-geographic scale. The estimated effective population size, the moderate level of dispersal and the rapid temporal change in genetic composition might reflect a potential effect of life history strategy on population dynamics. This hypothesis deserves further investigation. The dynamic pattern revealed at micro-geographic and temporal scales appears important from a genetic resource management as well as from a biodiversity conservation point of view.

Individual and population level effects of partner notification for Chlamydia trachomatis

Partner notification (PN or contact tracing) is an important aspect of treating bacterial sexually transmitted infections (STIs), such as Chlamydia trachomatis. It facilitates the identification of new infected cases that can be treated through individual case management. PN also acts indirectly by limiting onward transmission in the general population. However, the impact of PN, both at the level of individuals and the population, remains unclear. Since it is difficult to study the effects of PN empirically, mathematical and computational models are useful tools for investigating its potential as a public health intervention. To this end, we developed an individual-based modeling framework called Rstisim. It allows the implementation of different models of STI transmission with various levels of complexity and the reconstruction of the complete dynamic sexual partnership network over any time period. A key feature of this framework is that we can trace an individual's partnership history in detail and investigate the outcome of different PN strategies for C. trachomatis. For individual case management, the results suggest that notifying three or more partners from the preceding 18 months yields substantial numbers of new cases. In contrast, the successful treatment of current partners is most important for preventing re-infection of index cases and reducing further transmission of C. trachomatis at the population level. The findings of this study demonstrate the difference between individual and population level outcomes of public health interventions for STIs.

Simulating the impact of four control strategies on the population dynamics of Neospora caninum infection in Swiss dairy cattle

A dynamic deterministic simulation model was developed to assess the impact of different putative control strategies on the seroprevalence of Neospora caninum in female Swiss dairy cattle. The model structure comprised compartments of "susceptible" and "infected" animals (SI-model) and the cattle population was divided into 12 age classes. A reference model (Model 1) was developed to simulate the current (status quo) situation (present seroprevalence in Switzerland 12%), taking into account available demographic and seroprevalence data of Switzerland. Model 1 was modified to represent four putative control strategies: testing and culling of seropositive animals (Model 2), discontinued breeding with offspring from seropositive cows (Model 3), chemotherapeutic treatment of calves from seropositive cows (Model 4), and vaccination of susceptible and infected animals (Model 5). Models 2-4 considered different sub-scenarios with regard to the frequency of diagnostic testing. Multivariable Monte Carlo sensitivity analysis was used to assess the impact of uncertainty in input parameters. A policy of annual testing and culling of all seropositive cattle in the population reduced the seroprevalence effectively and rapidly from 12% to <1% in the first year of simulation. The control strategies with discontinued breeding with offspring from all seropositive cows, chemotherapy of calves and vaccination of all cattle reduced the prevalence more slowly than culling but were still very effective (reduction of prevalence below 2% within 11, 23 and 3 years of simulation, respectively). However, sensitivity analyses revealed that the effectiveness of these strategies depended strongly on the quality of the input parameters used, such as the horizontal and vertical transmission factors, the sensitivity of the diagnostic test and the efficacy of medication and vaccination. Finally, all models confirmed that it was not possible to completely eradicate N. caninum as long as the horizontal transmission process was not interrupted.

Cost effectiveness of home based population screening for Chlamydia trachomatis in the UK: economic evaluation of chlamydia screening studies (ClaSS) project

OBJECTIVE: To investigate the cost effectiveness of screening for Chlamydia trachomatis compared with a policy of no organised screening in the United Kingdom. DESIGN: Economic evaluation using a transmission dynamic mathematical model. SETTING: Central and southwest England. PARTICIPANTS: Hypothetical population of 50,000 men and women, in which all those aged 16-24 years were invited to be screened each year. MAIN OUTCOME MEASURES: Cost effectiveness based on major outcomes averted, defined as pelvic inflammatory disease, ectopic pregnancy, infertility, or neonatal complications. RESULTS: The incremental cost per major outcome averted for a programme of screening women only (assuming eight years of screening) was 22,300 pounds (33,000 euros; $45,000) compared with no organised screening. For a programme screening both men and women, the incremental cost effectiveness ratio was approximately 28,900 pounds. Pelvic inflammatory disease leading to hospital admission was the most frequently averted major outcome. The model was highly sensitive to the incidence of major outcomes and to uptake of screening. When both were increased the cost effectiveness ratio fell to 6200 pound per major outcome averted for screening women only. CONCLUSIONS: Proactive register based screening for chlamydia is not cost effective if the uptake of screening and incidence of complications are based on contemporary empirical studies, which show lower rates than commonly assumed. These data are relevant to discussions about the cost effectiveness of the opportunistic model of chlamydia screening being introduced in England.

Epidemiological, social, diagnostic and economic evaluation of population screening for genital chlamydial infection

OBJECTIVES: To investigate epidemiological, social, diagnostic and economic aspects of chlamydia screening in non-genitourinary medicine settings. METHODS: Linked studies around a cross-sectional population-based survey of adult men and women invited to collect urine and (for women) vulvovaginal swab specimens at home and mail these to a laboratory for testing for Chlamydia trachomatis. Specimens were used in laboratory evaluations of an amplified enzyme immunoassay (PCE EIA) and two nucleic acid amplification tests [Cobas polymerase chain reaction (PCR), Becton Dickinson strand displacement amplification (SDA)]. Chlamydia-positive cases and two negative controls completed a risk factor questionnaire. Chlamydia-positive cases were invited into a randomised controlled trial of partner notification strategies. Samples of individuals testing negative completed psychological questionnaires before and after screening. In-depth interviews were conducted at all stages of screening. Chlamydia transmission and cost-effectiveness of screening were investigated in a transmission dynamic model. SETTING AND PARTICIPANTS: General population in the Bristol and Birmingham areas of England. In total, 19,773 women and men aged 16-39 years were randomly selected from 27 general practice lists. RESULTS: Screening invitations reached 73% (14,382/19,773). Uptake (4731 participants), weighted for sampling, was 39.5% (95% CI 37.7, 40.8%) in women and 29.5% (95% CI 28.0, 31.0%) in men aged 16-39 years. Chlamydia prevalence (219 positive results) in 16-24 year olds was 6.2% (95% CI 4.9, 7.8%) in women and 5.3% (95% CI 4.4, 6.3%) in men. The case-control study did not identify any additional factors that would help target screening. Screening did not adversely affect anxiety, depression or self-esteem. Participants welcomed the convenience and privacy of home-sampling. The relative sensitivity of PCR on male urine specimens was 100% (95% CI 89.1, 100%). The combined relative sensitivities of PCR and SDA using female urine and vulvovaginal swabs were 91.8% (86.1, 95.7, 134/146) and 97.3% (93.1, 99.2%, 142/146). A total of 140 people (74% of eligible) participated in the randomised trial. Compared with referral to a genitourinary medicine clinic, partner notification by practice nurses resulted in 12.4% (95% CI -3.7, 28.6%) more patients with at least one partner treated and 22.0% (95% CI 6.1, 37.8%) more patients with all partners treated. The health service and patients costs (2005 prices) of home-based postal chlamydia screening were 21.47 pounds (95% CI 19.91 pounds, 25.99) per screening invitation and 28.56 pounds (95% CI 22.10 pounds, 30.43) per accepted offer. Preliminary modelling found an incremental cost-effectiveness ratio (2003 prices) comparing screening men and women annually to no screening in the base case of 27,000 pounds/major outcome averted at 8 years. If estimated screening uptake and pelvic inflammatory disease incidence were increased, the cost-effectiveness ratio fell to 3700 pounds/major outcome averted. CONCLUSIONS: Proactive screening for chlamydia in women and men using home-collected specimens was feasible and acceptable. Chlamydia prevalence rates in men and women in the general population are similar. Nucleic acid amplification tests can be used on first-catch urine specimens and vulvovaginal swabs. The administrative costs of proactive screening were similar to those for opportunistic screening. Using empirical estimates of screening uptake and incidence of complications, screening was not cost-effective.

Predicting the population impact of chlamydia screening programmes: comparative mathematical modelling study

BACKGROUND: Published individual-based, dynamic sexual network modelling studies reach different conclusions about the population impact of screening for Chlamydia trachomatis. The objective of this study was to conduct a direct comparison of the effect of organised chlamydia screening in different models. METHODS: Three models simulating population-level sexual behaviour, chlamydia transmission, screening and partner notification were used. Parameters describing a hypothetical annual opportunistic screening program in 16-24 year olds were standardised, whereas other parameters from the three original studies were retained. Model predictions of the change in chlamydia prevalence were compared under a range of scenarios. RESULTS: Initial overall chlamydia prevalence rates were similar in women but not men and there were age and sex-specific differences between models. The number of screening tests carried out was comparable in all models but there were large differences in the predicted impact of screening. After 10 years of screening, the predicted reduction in chlamydia prevalence in women aged 16-44 years ranged from 4% to 85%. Screening men and women had a greater impact than screening women alone in all models. There were marked differences between models in assumptions about treatment seeking and sexual behaviour before the start of the screening intervention. CONCLUSIONS: Future models of chlamydia transmission should be fitted to both incidence and prevalence data. This meta-modelling study provides essential information for explaining differences between published studies and increasing the utility of individual-based chlamydia transmission models for policy making.

On the detection of dynamic responses in a drought-perturbed tropical rainforest in Borneo

The dynamics of aseasonal lowland dipterocarp forest in Borneo is influenced by perturbation from droughts. These events might be increasing in frequency and intensity in the future. This paper describes drought-affected dynamics between 1986 and 2001 in Sabah, Malaysia, and considers how it is possible, reliably and accurately, to measure both coarse- and fine-scale responses of the forest. Some fundamental concerns about methodology and data analysis emerge. In two plots forming 8 ha, mortality, recruitment, and stem growth rates of trees ≥10 cm gbh (girth at breast height) were measured in a ‘pre-drought’ period (1986–1996), and in a period (1996–2001) including the 1997–1998 ENSO-drought. For 2.56 ha of subplots, mortality and growth rates of small trees (10–<50 cm gbh) were found also for two sub-periods (1996–1999, 1999–2001). A total of c. 19 K trees were recorded. Mortality rate increased by 25% while both recruitment and relative growth rates increased by 12% for all trees at the coarse scale. For small trees, at the fine scale, mortality increased by 6% and 9% from pre-drought to drought and on to ‘post-drought’ sub-periods. Relative growth rates correspondingly decreased by 38% and increased by 98%. Tree size and topography interacted in a complex manner with between-plot differences. The forest appears to have been sustained by off-setting elevated tree mortality by highly resilient stem growth. This last is seen as the key integrating tree variable which links the external driver (drought causing water stress) and population dynamics recorded as mortality and recruitment. Suitably sound measurements of stem girth, leading to valid growth rates, are needed to understand and model tree dynamic responses to perturbations. The proportion of sound data, however, is in part determined by the drought itself.

Demographic model of the Swiss cattle population for the years 2009-2011 stratified by gender, age and production type

Demographic composition and dynamics of animal and human populations are important determinants for the transmission dynamics of infectious disease and for the effect of infectious disease or environmental disasters on productivity. In many circumstances, demographic data are not available or of poor quality. Since 1999 Switzerland has been recording cattle movements, births, deaths and slaughter in an animal movement database (AMD). The data present in the AMD offers the opportunity for analysing and understanding the dynamic of the Swiss cattle population. A dynamic population model can serve as a building block for future disease transmission models and help policy makers in developing strategies regarding animal health, animal welfare, livestock management and productivity. The Swiss cattle population was therefore modelled using a system of ordinary differential equations. The model was stratified by production type (dairy or beef), age and gender (male and female calves: 0-1 year, heifers and young bulls: 1-2 years, cows and bulls: older than 2 years). The simulation of the Swiss cattle population reflects the observed pattern accurately. Parameters were optimized on the basis of the goodness-of-fit (using the Powell algorithm). The fitted rates were compared with calculated rates from the AMD and differed only marginally. This gives confidence in the fitted rates of parameters that are not directly deductible from the AMD (e.g. the proportion of calves that are moved from the dairy system to fattening plants).

Growth in Virologically Suppressed HIV Positive Children on Antiretroviral Therapy: Individual and Population-Level References

BACKGROUND Combination antiretroviral therapy (ART) suppresses viral replication in HIV-infected children. The growth of virologically suppressed children on ART has not been well documented. We aimed to develop dynamic reference curves for weight-for-age z scores (WAZ) and height-for-age z scores (HAZ). RESULTS A total of 4,876 children were followed for 7,407 person-years. Analyses were stratified by baseline z-scores and age, which were the most important predictors of growth response. The youngest children showed the most pronounced increase in weight and height initially but catch-up growth stagnated after 1-2 years. Three years after starting ART, WAZ ranged from -2.2 (95% Prediction interval -5.6 to 0.8) in children with baseline age "5 years and z-score "-3 to 0.0 (-2.7 to 2.4) in children with baseline age "2 years and WAZ "-1. For HAZ the corresponding range was -2.3 (-4.9 to 0.3) in children with baseline age"5 years and z-score "-3 to 0.3 (-3.1 to 3.4) in children with baseline age 2-5 years and HAZ "-1. CONCLUSIONS We have developed an online tool to calculate reference trajectories in fully suppressed children. The web application could help to define 'optimal' growth response and identify children with treatment failure.