Collaborative nowcasting of COVID-19 hospitalization incidences in Germany.

Real-time surveillance is a crucial element in the response to infectious disease outbreaks. However, the interpretation of incidence data is often hampered by delays occurring at various stages of data gathering and reporting. As a result, recent values are biased downward, which obscures current trends. Statistical nowcasting techniques can be employed to correct these biases, allowing for accurate characterization of recent developments and thus enhancing situational awareness. In this paper, we present a preregistered real-time assessment of eight nowcasting approaches, applied by independent research teams to German 7-day hospitalization incidences during the COVID-19 pandemic. This indicator played an important role in the management of the outbreak in Germany and was linked to levels of non-pharmaceutical interventions via certain thresholds. Due to its definition, in which hospitalization counts are aggregated by the date of case report rather than admission, German hospitalization incidences are particularly affected by delays and can take several weeks or months to fully stabilize. For this study, all methods were applied from 22 November 2021 to 29 April 2022, with probabilistic nowcasts produced each day for the current and 28 preceding days. Nowcasts at the national, state, and age-group levels were collected in the form of quantiles in a public repository and displayed in a dashboard. Moreover, a mean and a median ensemble nowcast were generated. We find that overall, the compared methods were able to remove a large part of the biases introduced by delays. Most participating teams underestimated the importance of very long delays, though, resulting in nowcasts with a slight downward bias. The accompanying prediction intervals were also too narrow for almost all methods. Averaged over all nowcast horizons, the best performance was achieved by a model using case incidences as a covariate and taking into account longer delays than the other approaches. For the most recent days, which are often considered the most relevant in practice, a mean ensemble of the submitted nowcasts performed best. We conclude by providing some lessons learned on the definition of nowcasting targets and practical challenges.

Mobile vaccination units to increase COVID-19 vaccination uptake in areas with lower coverage: a within-neighbourhood analysis using national registration data, the Netherlands, September-December 2021.

BackgroundVaccine uptake differs between social groups. Mobile vaccination units (MV-units) were deployed in the Netherlands by municipal health services in neighbourhoods with low uptake of COVID-19 vaccines.AimWe aimed to evaluate the impact of MV-units on vaccine uptake in neighbourhoods with low vaccine uptake.MethodsWe used the Dutch national-level registry of COVID-19 vaccinations (CIMS) and MV-unit deployment registrations containing observations in 253 neighbourhoods where MV-units were deployed and 890 contiguous neighbourhoods (total observations: 88,543 neighbourhood-days). A negative binomial regression with neighbourhood-specific temporal effects using splines was used to study the effect.ResultsDuring deployment, the increase in daily vaccination rate in targeted neighbourhoods ranged from a factor 2.0 (95% confidence interval (CI): 1.8-2.2) in urbanised neighbourhoods to 14.5 (95% CI: 11.6-18.0) in rural neighbourhoods. The effects were larger in neighbourhoods with more voters for the Dutch conservative Reformed Christian party but smaller in neighbourhoods with a higher proportion of people with non-western migration backgrounds. The absolute increase in uptake over the complete intervention period ranged from 0.22 percentage points (95% CI: 0.18-0.26) in the most urbanised neighbourhoods to 0.33 percentage point (95% CI: 0.28-0.37) in rural neighbourhoods.ConclusionDeployment of MV-units increased daily vaccination rate, particularly in rural neighbourhoods, with longer travel distance to permanent vaccination locations. This public health intervention shows promise to reduce geographic and social health inequalities, but more proactive and long-term deployment is required to identify its potential to substantially contribute to overall vaccination rates at country level.

Seroprevalence of Toxoplasma gondii and associated risk factors for infection in the Netherlands: third cross-sectional national study.

A third nationally representative serosurvey was performed to study the changes in Toxoplasma gondii (T. gondii) seroprevalence in the Netherlands over a 20-year time span and to identify and confirm risk factors for acquired toxoplasmosis. This cross-sectional study (conducted in 2016/2017) was designed similarly to the previous two studies (1995/1996 and 2006/2007) and included a questionnaire and serum sampling among Dutch residents. Factors associated with seropositivity for T. gondii were determined using multivariable analysis of the questionnaire-derived data. The earlier observed decrease in T. gondii seroprevalence between 1995/1996 and 2006/2007 (from 40.5% to 26.0%) did not continue into 2016/2017 (29.9%). Similarly to the previous studies, the seroprevalence increased with age and varied among regions. In all studies, higher T. gondii seropositivity was associated with increasing age, lower educational level, not living in the Southeast, and eating raw or semi-cooked pork. The incidence of congenital toxoplasmosis was estimated at 1.3/1000 (95% CI 0.9-1.8) live-born children in 2017. As the seroprevalence of T. gondii in the Netherlands did not decrease over the last decade, an increase in public health awareness is needed and prevention measures may need to be taken to achieve a further reduction in T. gondii infections in the Netherlands.

Contact patterns of older adults with and without frailty in the Netherlands during the COVID-19 pandemic.

BACKGROUND: During the COVID-19 pandemic, social distancing measures were imposed to protect the population from exposure, especially older adults and people with frailty, who have the highest risk for severe outcomes. These restrictions greatly reduced contacts in the general population, but little was known about behaviour changes among older adults and people with frailty themselves. Our aim was to quantify how COVID-19 measures affected the contact behaviour of older adults and how this differed between older adults with and without frailty. METHODS: In 2021, a contact survey was carried out among people aged 70 years and older in the Netherlands. A random sample of persons per age group (70-74, 75-79, 80-84, 85-89, and 90 +) and gender was invited to participate, either during a period with stringent (April 2021) or moderate (October 2021) measures. Participants provided general information on themselves, including their frailty, and they reported characteristics of all persons with whom they had face-to-face contact on a given day over the course of a full week. RESULTS: In total, 720 community-dwelling older adults were included (overall response rate of 15%), who reported 16,505 contacts. During the survey period with moderate measures, participants without frailty had significantly more contacts outside their household than participants with frailty. Especially for females, frailty was a more informative predictor of the number of contacts than age. During the survey period with stringent measures, participants with and without frailty had significantly lower numbers of contacts compared to the survey period with moderate measures. The reduction of the number of contacts was largest for the eldest participants without frailty. As they interact mostly with adults of a similar high age who are likely frail, this reduction of the number of contacts indirectly protects older adults with frailty from SARS-CoV-2 exposure. CONCLUSIONS: The results of this study reveal that social distancing measures during the COVID-19 pandemic differentially affected the contact patterns of older adults with and without frailty. The reduction of contacts may have led to the direct protection of older adults in general but also to the indirect protection of older adults with frailty.

Neighbourhood sociodemographic factors and COVID-19 vaccine uptake in the Netherlands: an ecological analysis.

BACKGROUND: While overall COVID-19 vaccine uptake is high in the Netherlands, it lags behind in certain subpopulations. AIM: We aimed to explore the characteristics of groups with lower COVID-19 vaccine uptake at neighbourhood level to inform the strategy to improve uptake and guide research into barriers for vaccination. METHODS: We performed an ecological study using national vaccination register and socio-demographic data at neighbourhood level. Using univariate and multivariable generalized additive models we examined the (potentially non-linear) effect of each determinant on uptake. We focused on those aged 50 years and older, since they are at highest risk of severe disease. RESULTS: In those over 50 years of age, a higher proportion of individuals with a non-Western migration background and higher voting proportions for right-wing Christian and conservative political parties were at neighbourhood level univariately associated with lower COVID-19 vaccine uptake. In contrast, higher socioeconomic status and higher voting proportions for right-wing liberal, progressive liberal and Christian middle political parties were associated with higher uptake. Multivariable results differed from univariate results in that a higher voting proportion for progressive left-wing political parties was also associated with higher uptake. In addition, with regard to migration background only a Turkish background remained significant. CONCLUSION: We identified determinants associated with COVID-19 vaccine uptake at neighbourhood level and observed heterogeneity in uptake between different subpopulations. Since the goal of vaccination is not only to reduce suffering and death by improving the average uptake, but also to reduce health inequity, it is important to focus on subpopulations with lower uptake.

A machine learning approach to small area estimation: predicting the health, housing and well-being of the population of Netherlands.

BACKGROUND: Local policymakers require information about public health, housing and well-being at small geographical areas. A municipality can for example use this information to organize targeted activities with the aim of improving the well-being of their residents. Surveys are often used to gather data, but many neighborhoods can have only few or even zero respondents. In that case, estimating the status of the local population directly from survey responses is prone to be unreliable. METHODS: Small Area Estimation (SAE) is a technique to provide estimates at small geographical levels with only few or even zero respondents. In classical individual-level SAE, a complex statistical regression model is fitted to the survey responses by using auxiliary administrative data for the population as predictors, the missing responses are then predicted and aggregated to the desired geographical level. In this paper we compare gradient boosted trees (XGBoost), a well-known machine learning technique, to a structured additive regression model (STAR) designed for the specific problem of estimating public health and well-being in the whole population of the Netherlands. RESULTS: We compare the accuracy and performance of these models using out-of-sample predictions with five-fold Cross Validation (5CV). We do this for three data sets of different sample sizes and outcome types. Compared to the STAR model, gradient boosted trees are able to improve both the accuracy of the predictions and the total time taken to get these predictions. Even though the models appear quite similar in overall accuracy, the small area predictions at neighborhood level sometimes differ significantly. It may therefore make sense to pursue slightly more accurate models for better predictions into small areas. However, one of the biggest benefits is that XGBoost does not require prior knowledge or model specification. Data preparation and modelling is much easier, since the method automatically handles missing data, non-linear responses, interactions and accounts for spatial correlation structures. CONCLUSIONS: In this paper we provide new nationwide estimates of health, housing and well-being indicators at neighborhood level in the Netherlands, see 'Online materials'. We demonstrate that machine learning provides a good alternative to complex statistical regression modelling for small area estimation in terms of accuracy, robustness, speed and data preparation. These results can be used to make appropriate policy decisions at a local level and make recommendations about which estimation methods are beneficial in terms of accuracy, time and budget constraints.

Excess Deaths during Influenza and Coronavirus Disease and Infection-Fatality Rate for Severe Acute Respiratory Syndrome Coronavirus 2, the Netherlands.

Since the 2009 influenza pandemic, the Netherlands has used a weekly death monitoring system to estimate deaths in excess of expectations. We present estimates of excess deaths during the ongoing coronavirus disease (COVID-19) epidemic and 10 previous influenza epidemics. Excess deaths per influenza epidemic averaged 4,000. The estimated 9,554 excess deaths (41% in excess) during the COVID-19 epidemic weeks 12-19 of 2020 appeared comparable to the 9,373 excess deaths (18%) during the severe influenza epidemic of 2017-18. However, these deaths occurred in a shorter time, had a higher peak, and were mitigated by nonpharmaceutical control measures. Excess deaths were 1.8-fold higher than reported laboratory-confirmed COVID-19 deaths (5,449). Based on excess deaths and preliminary results from seroepidemiologic studies, we estimated the infection-fatality rate to be 1%. Monitoring of excess deaths is crucial for timely estimates of disease burden for influenza and COVID-19. Our data complement laboratory-confirmed COVID-19 death reports and enable comparisons between epidemics.

Trade-off between local transmission and long-range dispersal drives infectious disease outbreak size in spatially structured populations.

Transmission of infectious diseases between immobile hosts (e.g., plants, farms) is strongly dependent on the spatial distribution of hosts and the distance-dependent probability of transmission. As the interplay between these factors is poorly understood, we use spatial process and transmission modelling to investigate how epidemic size is shaped by host clustering and spatial range of transmission. We find that for a given degree of clustering and individual-level infectivity, the probability that an epidemic occurs after an introduction is generally higher if transmission is predominantly local. However, local transmission also impedes transfer of the infection to new clusters. A consequence is that the total number of infections is maximal if the range of transmission is intermediate. In highly clustered populations, the infection dynamics is strongly determined by the probability of transmission between clusters of hosts, whereby local clusters act as multiplier of infection. We show that in such populations, a metapopulation model sometimes provides a good approximation of the total epidemic size, using probabilities of local extinction, the final size of infections in local clusters, and probabilities of cluster-to-cluster transmission. As a real-world example we analyse the case of avian influenza transmission between poultry farms in the Netherlands.

Vaccine effectiveness against SARS-CoV-2 transmission and infections among household and other close contacts of confirmed cases, the Netherlands, February to May 2021.

Several studies report high effectiveness of COVID-19 vaccines against SARS-CoV-2 infection and severe disease, however an important knowledge gap is the vaccine effectiveness against transmission (VET). We present estimates of the VET to household and other close contacts in the Netherlands, from February to May 2021, using contact monitoring data. The secondary attack rate among household contacts was lower for fully vaccinated than unvaccinated index cases (11% vs 31%), with an adjusted VET of 71% (95% confidence interval: 63-77).

Use of Ambulance Dispatch Calls for Surveillance of Severe Acute Respiratory Infections.

Ambulance dispatches for respiratory syndromes reflect incidence of influenza-like illness in primary care. Associations are highest in children (15%-34% of respiratory calls attributable to influenza), out-of-office hours (9%), and highest urgency-level calls (9%-11%). Ambulance dispatches might be an additional source of data for severe influenza surveillance.

Tracing the Origin of Food-borne Disease Outbreaks: A Network Model Approach.

BACKGROUND: Food-borne disease outbreaks constitute a large health burden on society. One of the challenges when investigating such outbreaks is to trace the origin of the outbreak. In this study, we consider a network model to determine the spatial origin of the contaminated food product that caused the outbreak. METHODS: The network model we use replaces the classic geographic distance of a network by an effective distance so that two nodes connected by a long-range link may be more strongly connected than their geographic distance would suggest. Furthermore, the effective distance transforms complex spatial patterns into regular topological patterns, creating a means for easier identification of the origin of the spreading phenomenon. Because detailed information on food distribution is generally not available, the model uses the gravity model from economics: the flow of goods from one node to another increases with population size and decreases with the geographical distance between them. RESULTS: This effective distance network approach has been shown to perform well in a large Escherichia coli O104:H4 outbreak in Germany in 2011. In this article, we apply the same method to various food-borne disease outbreaks in the Netherlands. We found the effective distance network approach to fail in certain scenarios. CONCLUSIONS: Great care should be taken as to whether the underlying network model correctly captures the spreading mechanism of the outbreak in terms of spatial scale and single or multiple source outbreak.

Nowcasting the Number of New Symptomatic Cases During Infectious Disease Outbreaks Using Constrained P-spline Smoothing.

During an infectious disease outbreak, timely information on the number of new symptomatic cases is crucial. However, the reporting of new cases is usually subject to delay due to the incubation period, time to seek care, and diagnosis. This results in a downward bias in the numbers of new cases by the times of symptoms onset towards the current day. The real-time assessment of the current situation while correcting for underreporting is called nowcasting. We present a nowcasting method based on bivariate P-spline smoothing of the number of reported cases by time of symptoms onset and delay. Our objective is to predict the number of symptomatic-but-not-yet-reported cases and combine these with the already reported symptomatic cases into a nowcast. We assume the underlying two-dimensional reporting intensity surface to be smooth. We include prior information on the reporting process as additional constraints: the smooth surface is unimodal in the reporting delay dimension, is (almost) zero at a predefined maximum delay and has a prescribed shape at the beginning of the outbreak. Parameter estimation is done efficiently by penalized iterative weighted least squares. We illustrate our method on a large measles outbreak in the Netherlands. We show that even with very limited information the method is able to accurately predict the number of symptomatic-but-not-yet-reported cases. This results in substantially improved monitoring of new symptomatic cases in real time.

"Sickenin' in the rain" - increased risk of gastrointestinal and respiratory infections after urban pluvial flooding in a population-based cross-sectional study in the Netherlands.

BACKGROUND: Climate change is expected to increase the chance of extreme rainfall events in the Northern Hemisphere and herewith, there is an increased chance of urban pluvial flooding. Urban pluvial flooding often consists of street flooding and/or flooding of combined sewerage systems, leading to contamination of the floodwater with several gastrointestinal and/or respiratory pathogens. An increase in flooding events therefore pose a health risk to those exposed to urban floodwater. We studied the association between exposure to pluvial floodwater and acute gastroenteritis (AGE) and acute respiratory infection (ARI). METHODS: We performed a retrospective, cross-sectional survey during the summer of 2015 in 60 locations in the Netherlands with reported flooding. Two weeks after the flooding, questionnaires were sent to households in these locations, collecting data on self-reported AGE and ARI and information on floodwater exposure in the previous 2 weeks. Multivariable generalized estimating equations (GEE) regression models, accounting for the clustered data structure, were used to identify risk factors for AGE and ARI. RESULTS: In total, 699 households with 1,656 participants (response rate 21%) returned the questionnaire. Contact with floodwater was significantly associated with AGE (aOR 4.2, 95%CI 2.1-8.4) and ARI (aOR 3.3, 95%CI 2.0-5.4). Risk factors for AGE were skin contact with floodwater (aOR 4.0, 95%CI 1.8-9.0), performing post-flooding cleaning operations (aOR 8.6, 95%CI 3.5-20.9) and cycling through floodwater (aOR 2.3, 95%CI 1.0-5.0). Skin contact with floodwater (aOR 3.6, 95%CI 1.9-6.9) and performing post-flooding cleaning operations (aOR 5.5, 95%CI 3.0-10.3) were identified as risk factors for ARI. CONCLUSIONS: Results suggest an association between direct exposure to pluvial floodwater and AGE and ARI. As it is predicted that the frequency of pluvial flooding events will increase in the future, there is a need for flood-proof solutions in urban development and increased awareness among stakeholders and the public about the potential health risks. Future prospective studies are recommended to confirm our results.

Infectious reactivation of cytomegalovirus explaining age- and sex-specific patterns of seroprevalence.

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.

Determinants of Rotavirus Transmission: A Lag Nonlinear Time Series Analysis.

Rotavirus is a common viral infection among young children. As in many countries, the infection dynamics of rotavirus in the Netherlands are characterized by an annual winter peak, which was notably low in 2014. Previous study suggested an association between weather factors and both rotavirus transmission and incidence. From epidemic theory, we know that the proportion of susceptible individuals can affect disease transmission. We investigated how these factors are associated with rotavirus transmission in the Netherlands, and their impact on rotavirus transmission in 2014. We used available data on birth rates and rotavirus laboratory reports to estimate rotavirus transmission and the proportion of individuals susceptible to primary infection. Weather data were directly available from a central meteorological station. We developed an approach for detecting determinants of seasonal rotavirus transmission by assessing nonlinear, delayed associations between each factor and rotavirus transmission. We explored relationships by applying a distributed lag nonlinear regression model with seasonal terms. We corrected for residual serial correlation using autoregressive moving average errors. We inferred the relationship between different factors and the effective reproduction number from the most parsimonious model with low residual autocorrelation. Higher proportions of susceptible individuals and lower temperatures were associated with increases in rotavirus transmission. For 2014, our findings suggest that relatively mild temperatures combined with the low proportion of susceptible individuals contributed to lower rotavirus transmission in the Netherlands. However, our model, which overestimated the magnitude of the peak, suggested that other factors were likely instrumental in reducing the incidence that year.

A summary index for antimicrobial resistance in food animals in the Netherlands.

BACKGROUND: The Dutch government has set targets for reduction of antimicrobial usage in food animals, stipulating a 50% reduction in usage (on a weight basis) in 2013 as compared to 2009 and a 70% decrease in 2015. A monitoring program has been instituted to evaluate the impact on antimicrobial resistance (AMR). The Dutch Ministry of Public Health Welfare and Sports has expressed the need for a summary index to present the results of the monitoring data concisely to policy makers. METHODS: We use data on AMR in bacteria from randomly collected samples from broiler chickens, fattening pigs, veal calves and dairy cows. Escherichia coli was selected for resistance monitoring because they are intrinsically susceptible to the antibiotics included in the test panel (ciprofloxacin, cefotaxime, tetracycline and ampicillin) and they are present in all samples, which facilitates proper randomization and trend analysis. The AMR summary index was calculated for each animal species as a weighted average over the four antibiotics, taking into account their clinical relevance. Weights were obtained by conjoint analysis, a pairwise comparison study involving infectious diseases professionals with clinical and public health backgrounds, with data analysis by conditional logistic regression. The AMR summary index was then computed by Monte Carlo simulation, accounting for sampling and regression uncertainty. RESULTS: The highest weights (0.35) were given to ciprofloxacin and cefotaxime followed by ampicillin (0.23) and tetracycline (0.07). Throughout the years, the AMR index was highest in broiler chickens, followed by pigs and veal calves, while the lowest values were consistently recorded in dairy cows. In all animal species, the index in 2014 was significantly lower than in 2009. CONCLUSIONS: We demonstrate that high-dimensional data on surveillance of antimicrobial resistance can be summarized in an index for evaluating trends between and within food animal species by a process involving decision makers and scientists to select and weight the most relevant antibiotics.

Estimating the prevalence of 26 health-related indicators at neighbourhood level in the Netherlands using structured additive regression.

BACKGROUND: Local policy makers increasingly need information on health-related indicators at smaller geographic levels like districts or neighbourhoods. Although more large data sources have become available, direct estimates of the prevalence of a health-related indicator cannot be produced for neighbourhoods for which only small samples or no samples are available. Small area estimation provides a solution, but unit-level models for binary-valued outcomes that can handle both non-linear effects of the predictors and spatially correlated random effects in a unified framework are rarely encountered. METHODS: We used data on 26 binary-valued health-related indicators collected on 387,195 persons in the Netherlands. We associated the health-related indicators at the individual level with a set of 12 predictors obtained from national registry data. We formulated a structured additive regression model for small area estimation. The model captured potential non-linear relations between the predictors and the outcome through additive terms in a functional form using penalized splines and included a term that accounted for spatially correlated heterogeneity between neighbourhoods. The registry data were used to predict individual outcomes which in turn are aggregated into higher geographical levels, i.e. neighbourhoods. We validated our method by comparing the estimated prevalences with observed prevalences at the individual level and by comparing the estimated prevalences with direct estimates obtained by weighting methods at municipality level. RESULTS: We estimated the prevalence of the 26 health-related indicators for 415 municipalities, 2599 districts and 11,432 neighbourhoods in the Netherlands. We illustrate our method on overweight data and show that there are distinct geographic patterns in the overweight prevalence. Calibration plots show that the estimated prevalences agree very well with observed prevalences at the individual level. The estimated prevalences agree reasonably well with the direct estimates at the municipal level. CONCLUSIONS: Structured additive regression is a useful tool to provide small area estimates in a unified framework. We are able to produce valid nationwide small area estimates of 26 health-related indicators at neighbourhood level in the Netherlands. The results can be used for local policy makers to make appropriate health policy decisions.

Estimating seroprevalence of human papillomavirus type 16 using a mixture model with smoothed age-dependent mixing proportions.

BACKGROUND: The presence in serum of antibodies to viral antigens is generally considered a well-defined marker of past infection or vaccination. However, analyses of serological data that use a cut-off value to classify individuals as seropositive are prone to misclassification bias, in particular when studying infections with a weak serological response, such as the sexually transmitted human papillomavirus (HPV). METHODS: We analyzed the serological concentrations of HPV type 16 (HPV16) antibodies in the general Dutch population in 2006-2007, before the introduction of mass vaccination against HPV. We used a 2-component mixture model to represent persons who were seronegative or seropositive for HPV16. Component densities were assumed to be log-normally distributed, with parameters possibly dependent on sex. The age-dependent mixing proportions were smoothed using penalized splines to obtain a flexible seroprevalence profile. RESULTS: Our results suggest that HPV16 seropositivity is associated with higher antibody concentrations in women as compared with men. Seroprevalence shows an increase starting from adolescence in men and women alike, coinciding with the age of sexual debut. Seroprevalence stabilizes in men around age 40, whereas it has a decreasing trend from age 50 onwards in women. Analyses that rely on a cut-off value to classify persons as seropositive yield substantially different seroprevalence profiles, leading to a qualitatively different interpretation of HPV16 infection dynamics. CONCLUSIONS: Our results provide a benchmark for examining the effect of HPV16 vaccination in future serological surveys. Our method may prove useful for estimating seroprevalence of other infections with a weak serological response.

Improved correlation of human Q fever incidence to modelled C. burnetii concentrations by means of an atmospheric dispersion model.

BACKGROUND: Atmospheric dispersion models (ADMs) may help to assess human exposure to airborne pathogens. However, there is as yet limited quantified evidence that modelled concentrations are indeed associated to observed human incidence. METHODS: We correlated human Q fever (caused by the bacterium Coxiella burnetii) incidence data in the Netherlands to modelled concentrations from three spatial exposure models: 1) a NULL model with a uniform concentration distribution, 2) a DISTANCE model with concentrations proportional to the distance between the source and residential addresses of patients, and 3) concentrations modelled by an ADM using three simple emission profiles. We used a generalized linear model to correlate the observed incidences to modelled concentrations and validated it using cross-validation. RESULTS: ADM concentrations generally correlated the best to the incidence data. The DISTANCE model always performed significantly better than the NULL model. ADM concentrations based on wind speeds exceeding threshold values of 0 and 2 m/s performed better than those based on 4 or 6 m/s. This might indicate additional exposure to bacteria originating from a contaminated environment. CONCLUSIONS: By adding meteorological information the correlation between modelled concentration and observed incidence improved, despite using three simple emission profiles. Although additional information is needed - especially regarding emission data - these results provide a basis for the use of ADMs to predict and to visualize the spread of airborne pathogens during livestock, industry and even bio-terroristic related outbreaks or releases to a surrounding human population.

Vaccine uptake determinants in The Netherlands.

BACKGROUND: Combining existing data on background characteristics with data from immunization registers might give insight into determinants of vaccine uptake, which can help to improve communication strategies and invitation policy of National Immunisation Programmes. METHODS: The study population consisted of children born in 2005 as registered in the Dutch national immunization register Præventis. A hierarchical logistic regression model was used to quantify associations between individual vaccination status and proxy variables for ethnic background (individual level), socio-economic status (postcode level) and religious objection to vaccination (municipal level). RESULTS: Most children whose both parents were not born in The Netherlands had a somewhat lower full vaccine uptake, for example, children whose both parents were born in Turkey [odds ratio = 0.7 (0.6-0.8)] or in Morocco [odds ratio = 0.8 (0.7-0.9)]. The partial uptake was also relatively high (3.7-8.0%) compared with children whose both parents were born in The Netherlands (3.1%). Municipalities with higher religious objection to vaccination and postcode areas with lower socio-economic status were also associated with a lower full uptake. CONCLUSIONS: Despite the high vaccination coverage in The Netherlands, we were able to identify determinants of vaccine uptake by combining existing data sets. This might be an example for other countries. The impact of ethnic background and socio-economic status is not as well known in The Netherlands as the effect of religious objection to vaccination, and deserves more attention. Groups that have a relatively high partial uptake deserve special attention because they do not reject vaccination in general.