The number of new HIV infections among men who have sex with men (MSM) in the Netherlands has been decreasing, but additional efforts are required to bring it further down. This study aims to assess the impact of increased diagnosis of early HIV infection combined with immediate antiretroviral treatment (ART) initiation on reducing HIV transmission among MSM. We developed an agent-based model calibrated to HIV surveillance and sexual behavior data for MSM in the Netherlands in 2017-2022. Starting in 2023, we simulated a 10-year intervention that accelerates HIV diagnosis during the first 3 or 6 months after HIV acquisition across five levels of increased diagnosis rates (2, 4, 8, 16, and 32-fold), followed by immediate ART initiation. The upper limit of the intervention's impact over 10 years is projected to result in the cumulative 298 (95-th QI: 162-451) HIV infections averted. A 32-fold increase in the diagnosis rate within 3 months after HIV acquisition (corresponding to 100% of all new HIV infections diagnosed within 3 months of acquisition) results in 269 (95-th QI: 147-400) infections averted, approaching closely maximum impact. By extending the scope of the intervention to individuals who acquired HIV infection within the previous 6 months, a smaller 8-fold increase in the diagnosis rate (corresponding to 97% of new HIV infections diagnosed within 6 months of acquisition) approaches closely the maximum impact of the intervention by averting 256 (95-th QI: 122-411) HIV infections. Our sensitivity analyses showed that, in an epidemiological context similar to the modern-day the Netherlands, immediate initiation of ART accompanying accelerated diagnosis of individuals with early HIV infection does not significantly affect HIV transmission dynamics. Accelerating early HIV diagnosis through increased awareness, screening, and testing can further reduce transmission among MSM. Meeting this goal necessitates a stakeholder needs assessment.

Household studies provide an efficient means to study transmission of infectious diseases, enabling estimation of susceptibility and infectivity by person-type. A main inclusion criterion in such studies is usually the presence of an infected person. This precludes estimation of the hazards of pathogen introduction into the household. Here we estimate age- and time-dependent household introduction hazards together with within household transmission rates using data from a prospective household-based study in the Netherlands. A total of 307 households containing 1,209 persons were included from August 2020 until March 2021. Follow-up of households took place between August 2020 and August 2021 with maximal follow-up per household mostly limited to 161 days. Almost 1 out of 5 households (59/307) had evidence of an introduction of SARS-CoV-2. We estimate introduction hazards and within-household transmission rates in our study population with penalized splines and stochastic epidemic models, respectively. The estimated hazard of introduction of SARS-CoV-2 in the households was lower for children (0-12 years) than for adults (relative hazard: 0.62; 95%CrI: 0.34-1.0). Estimated introduction hazards peaked in mid October 2020, mid December 2020, and mid April 2021, preceding peaks in hospital admissions by 1-2 weeks. Best fitting transmission models included increased infectivity of children relative to adults and adolescents, such that the estimated child-to-child transmission probability (0.62; 95%CrI: 0.40-0.81) was considerably higher than the adult-to-adult transmission probability (0.12; 95%CrI: 0.057-0.19). Scenario analyses indicate that vaccination of adults can strongly reduce household infection attack rates and that adding adolescent vaccination offers limited added benefit.

Knowledge of who infected whom during an outbreak of an infectious disease is important to determine risk factors for transmission and to design effective control measures. Both whole-genome sequencing of pathogens and epidemiological data provide useful information about the transmission events and underlying processes. Existing models to infer transmission trees usually assume that the pathogen is introduced only once from outside into the population of interest. However, this is not always true. For instance, SARS-CoV-2 is suggested to be introduced multiple times in mink farms in the Netherlands from the SARS-CoV-2 pandemic among humans. Here, we developed a Bayesian inference method combining whole-genome sequencing data and epidemiological data, allowing for multiple introductions of the pathogen in the population. Our method does not a priori split the outbreak into multiple phylogenetic clusters, nor does it break the dependency between the processes of mutation, within-host dynamics, transmission, and observation. We implemented our method as an additional feature in the R-package phybreak. On simulated data, our method correctly identifies the number of introductions, with an accuracy depending on the proportion of all observed cases that are introductions. Moreover, when a single introduction was simulated, our method produced similar estimates of parameters and transmission trees as the existing package. When applied to data from a SARS-CoV-2 outbreak in Dutch mink farms, the method provides strong evidence for independent introductions of the pathogen at 13 farms, infecting a total of 63 farms. Using the new feature of the phybreak package, transmission routes of a more complex class of infectious disease outbreaks can be inferred which will aid infection control in future outbreaks.

The aim of this study is to analyze patient movement patterns between hospital departments to derive the underlying intra-hospital movement network, and to assess if movement patterns differ between patients at high or low risk of colonization. For that purpose, we analyzed patient electronic medical record data from five hospitals to extract information on risk stratification and patient intra-hospital movements. Movement patterns were visualized as networks, and network centrality measures were calculated. Next, using an agent-based model where agents represent patients and intra-hospital patient movements were explicitly modeled, we simulated the spread of multidrug resistant enterobacteriacae (MDR-E) inside a hospital. Risk stratification of patients according to certain ICD-10 codes revealed that length of stay, patient age, and mean number of movements per admission were higher in the high-risk groups. Movement networks in all hospitals displayed a high variability among departments concerning their network centrality and connectedness with a few highly connected departments and many weakly connected peripheral departments. Simulating the spread of a pathogen in one hospital network showed positive correlation between department prevalence and network centrality measures. This study highlights the importance of intra-hospital patient movements and their possible impact on pathogen spread. Targeting interventions to departments of higher (weighted) degree may help to control the spread of MDR-E. Moreover, when the colonization status of patients coming from different departments is unknown, a ranking system based on department centralities may be used to design more effective interventions that mitigate pathogen spread.

In the year 2020, there were 105 different statutory insurance companies in Germany with heterogeneous regional coverage. Obtaining data from all insurance companies is challenging, so that it is likely that projects will have to rely on data not covering the whole population. Consequently, the study of epidemic spread in hospital referral networks using data-driven models may be biased. We studied this bias using data from three German regional insurance companies covering four federal states: AOK (historically "general local health insurance company", but currently only the abbreviation is used) Lower Saxony (in Federal State of Lower Saxony), AOK Bavaria (in Bavaria), and AOK PLUS (in Thuringia and Saxony). To understand how incomplete data influence network characteristics and related epidemic simulations, we created sampled datasets by randomly dropping a proportion of patients from the full datasets and replacing them with random copies of the remaining patients to obtain scale-up datasets to the original size. For the sampled and scale-up datasets, we calculated several commonly used network measures, and compared them to those derived from the original data. We found that the network measures (degree, strength and closeness) were rather sensitive to incompleteness. Infection prevalence as an outcome from the applied susceptible-infectious-susceptible (SIS) model was fairly robust against incompleteness. At incompleteness levels as high as 90% of the original datasets the prevalence estimation bias was below 5% in scale-up datasets. Consequently, a coverage as low as 10% of the local population of the federal state population was sufficient to maintain the relative bias in prevalence below 10% for a wide range of transmission parameters as encountered in clinical settings. Our findings are reassuring that despite incomplete coverage of the population, German health insurance data can be used to study effects of patient traffic between institutions on the spread of pathogens within healthcare networks.

Inter-hospital patient transfers (direct transfers) between healthcare facilities have been shown to contribute to the spread of pathogens in a healthcare network. However, the impact of indirect transfers (patients re-admitted from the community to the same or different hospital) is not well studied. This work aims to study the contribution of indirect transfers to the spread of pathogens in a healthcare network. To address this aim, a hybrid network-deterministic model to simulate the spread of multiresistant pathogens in a healthcare system was developed for the region of Lower Saxony (Germany). The model accounts for both, direct and indirect transfers of patients. Intra-hospital pathogen transmission is governed by a SIS model expressed by a system of ordinary differential equations. Our results show that the proposed model reproduces the basic properties of healthcare-associated pathogen spread. They also show the importance of indirect transfers: restricting the pathogen spread to direct transfers only leads to 4.2% system wide prevalence. However, adding indirect transfers leads to an increase in the overall prevalence by a factor of 4 (18%). In addition, we demonstrated that the final prevalence in the individual healthcare facilities depends on average length of stay in a way described by a non-linear concave function. Moreover, we demonstrate that the network parameters of the model may be derived from administrative admission/discharge records. In particular, they are sufficient to obtain inter-hospital transfer probabilities, and to express the patients' transfers as a Markov process. Using the proposed model, we show that indirect transfers of patients are equally or even more important as direct transfers for the spread of pathogens in a healthcare network.

Pseudomonas aeruginosa (P. aeruginosa) is an important cause of healthcare-associated infections, particularly in immunocompromised patients. Understanding how this multi-drug resistant pathogen is transmitted within intensive care units (ICUs) is crucial for devising and evaluating successful control strategies. While it is known that moist environments serve as natural reservoirs for P. aeruginosa, there is little quantitative evidence regarding the contribution of environmental contamination to its transmission within ICUs. Previous studies on other nosocomial pathogens rely on deploying specific values for environmental parameters derived from costly and laborious genotyping. Using solely longitudinal surveillance data, we estimated the relative importance of P. aeruginosa transmission routes by exploiting the fact that different routes cause different pattern of fluctuations in the prevalence. We developed a mathematical model including background transmission, cross-transmission and environmental contamination. Patients contribute to a pool of pathogens by shedding bacteria to the environment. Natural decay and cleaning of the environment lead to a reduction of that pool. By assigning the bacterial load shed during an ICU stay to cross-transmission, we were able to disentangle environmental contamination during and after a patient's stay. Based on a data-augmented Markov Chain Monte Carlo method the relative importance of the considered acquisition routes is determined for two ICUs of the University hospital in Besançon (France). We used information about the admission and discharge days, screening days and screening results of the ICU patients. Both background and cross-transmission play a significant role in the transmission process in both ICUs. In contrast, only about 1% of the total transmissions were due to environmental contamination after discharge. Based on longitudinal surveillance data, we conclude that cleaning improvement of the environment after discharge might have only a limited impact regarding the prevention of P.A. infections in the two considered ICUs of the University hospital in Besançon. Our model was developed for P. aeruginosa but can be easily applied to other pathogens as well.

Human cytomegalovirus (CMV) is a herpes virus with poorly understood transmission dynamics. Person-to-person transmission is thought to occur primarily through transfer of saliva or urine, but no quantitative estimates are available for the contribution of different infection routes. Using data from a large population-based serological study (n = 5,179), we provide quantitative estimates of key epidemiological parameters, including the transmissibility of primary infection, reactivation, and re-infection. Mixture models are fitted to age- and sex-specific antibody response data from the Netherlands, showing that the data can be described by a model with three distributions of antibody measurements, i.e. uninfected, infected, and infected with increased antibody concentration. Estimates of seroprevalence increase gradually with age, such that at 80 years 73% (95%CrI: 64%-78%) of females and 62% (95%CrI: 55%-68%) of males are infected, while 57% (95%CrI: 47%-67%) of females and 37% (95%CrI: 28%-46%) of males have increased antibody concentration. Merging the statistical analyses with transmission models, we find that models with infectious reactivation (i.e. reactivation that can lead to the virus being transmitted to a novel host) fit the data significantly better than models without infectious reactivation. Estimated reactivation rates increase from low values in children to 2%-4% per year in women older than 50 years. The results advance a hypothesis in which transmission from adults after infectious reactivation is a key driver of transmission. We discuss the implications for control strategies aimed at reducing CMV infection in vulnerable groups.

The WHO's early-release guideline for antiretroviral treatment (ART) of HIV infection based on a recent trial conducted in 34 countries recommends starting treatment immediately upon an HIV diagnosis. Therefore, the test-and-treat strategy may become more widely used in an effort to scale up HIV treatment and curb further transmission. Here we examine behavioural determinants of HIV transmission and how heterogeneity in sexual behaviour influences the outcomes of this strategy. Using a deterministic model, we perform a systematic investigation into the effects of various mixing patterns in a population of men who have sex with men (MSM), stratified by partner change rates, on the elimination threshold and endemic HIV prevalence. We find that both the level of overdispersion in the distribution of the number of sexual partners and mixing between population subgroups have a large influence on endemic prevalence before introduction of ART and on possible long term effectiveness of ART. Increasing heterogeneity in risk behavior may lead to lower endemic prevalence levels, but requires higher coverage levels of ART for elimination. Elimination is only feasible for populations with a rather low degree of assortativeness of mixing and requires treatment coverage of almost 80% if rates of testing and treatment uptake by all population subgroups are equal. In this case, for fully assortative mixing and 80% coverage endemic prevalence is reduced by 57%. In the presence of heterogeneity in ART uptake, elimination is easier to achieve when the subpopulation with highest risk behavior is tested and treated more often than the rest of the population, and vice versa when it is less. The developed framework can be used to extract information on behavioral heterogeneity from existing data which is otherwise hard to determine from population surveys.

There are four major quantities that are measured in sexual behavior surveys that are thought to be especially relevant for the performance of sexual network models in terms of disease transmission. These are (i) the cumulative distribution of lifetime number of partners, (ii) the distribution of partnership durations, (iii) the distribution of gap lengths between partnerships, and (iv) the number of recent partners. Fitting a network model to these quantities as measured in sexual behavior surveys is expected to result in a good description of Chlamydia trachomatis transmission in terms of the heterogeneity of the distribution of infection in the population. Here we present a simulation model of a sexual contact network, in which we explored the role of behavioral heterogeneity of simulated individuals on the ability of the model to reproduce population-level sexual survey data from the Netherlands and UK. We find that a high level of heterogeneity in the ability of individuals to acquire and maintain (additional) partners strongly facilitates the ability of the model to accurately simulate the powerlaw-like distribution of the lifetime number of partners, and the age at which these partnerships were accumulated, as surveyed in actual sexual contact networks. Other sexual network features, such as the gap length between partnerships and the partnership duration, could-at the current level of detail of sexual survey data against which they were compared-be accurately modeled by a constant value (for transitional concurrency) and by exponential distributions (for partnership duration). Furthermore, we observe that epidemiological measures on disease prevalence in survey data can be used as a powerful tool for building accurate sexual contact networks, as these measures provide information on the level of mixing between individuals of different levels of sexual activity in the population, a parameter that is hard to acquire through surveying individuals.