The social contact hypothesis under the assumption of endemic equilibrium: Elucidating the transmission potential of VZV in Europe

Exposure-response relationship between temperature, relative humidity, and varicella: a multicity study in South China

Varicella is a rising public health issue. Several studies have tried to quantify the relationships between meteorological factors and varicella incidence but with inconsistent results. We aim to investigate the impact of temperature and relative humidity on varicella, and to further explore the effect modification of these relationships. In this study, the data of varicella and meteorological factors from 2011 to 2019 in 21 cities of Guangdong Province, China were collected. Distributed lag nonlinear models (DLNM) were constructed to explore the relationship between meteorological factors (temperature and relative humidity) and varicella in each city, controlling in school terms, holidays, seasonality, long-term trends, and day of week. Multivariate meta-analysis was applied to pool the city-specific estimations. And the meta-regression was used to explore the effect modification for the spatial heterogeneity of city-specific meteorological factors and social factors (such as disposable income per capita, vaccination coverage, and so on) on varicella. The results indicated that the relationship between temperature and varicella in 21 cities appeared nonlinear with an inverted S-shaped. The relative risk peaked at 20.8 ℃ (RR = 1.42, 95% CI: 1.22, 1.65). The relative humidity-varicella relationship was approximately L-shaped, with a peaking risk at 69.5% relative humidity (RR = 1.25, 95% CI: 1.04, 1.50). The spatial heterogeneity of temperature-varicella relationships may be caused by income or varicella vaccination coverage. And varicella vaccination coverage may contribute to the spatial heterogeneity of the relative humidity-varicella relationship. The findings can help us deepen the understanding of the meteorological factors-varicella association and provide evidence for developing prevention strategy for varicella epidemic.

Inferring age-specific differences in susceptibility to and infectiousness upon SARS-CoV-2 infection based on Belgian social contact data

Several important aspects related to SARS-CoV-2 transmission are not well known due to a lack of appropriate data. However, mathematical and computational tools can be used to extract part of this information from the available data, like some hidden age-related characteristics. In this paper, we present a method to investigate age-specific differences in transmission parameters related to susceptibility to and infectiousness upon contracting SARS-CoV-2 infection. More specifically, we use panel-based social contact data from diary-based surveys conducted in Belgium combined with the next generation principle to infer the relative incidence and we compare this to real-life incidence data. Comparing these two allows for the estimation of age-specific transmission parameters. Our analysis implies the susceptibility in children to be around half of the susceptibility in adults, and even lower for very young children (preschooler). However, the probability of adults and the elderly to contract the infection is decreasing throughout the vaccination campaign, thereby modifying the picture over time.

Estimates of the basic reproduction number for rubella using seroprevalence data and indicator-based approaches

The basic reproduction number (R0) of an infection determines the impact of its control. For many endemic infections, R0 is often estimated from appropriate country-specific seroprevalence data. Studies sometimes pool estimates from the same region for settings lacking seroprevalence data, but the reliability of this approach is unclear. Plausibly, indicator-based approaches could predict R0 for such settings. We calculated R0 for rubella for 98 settings and correlated its value against 66 demographic, economic, education, housing and health-related indicators. We also trained a random forest regression algorithm using these indicators as the input and R0 as the output. We used the mean-square error to compare the performances of the random forest, simple linear regression and a regional averaging method in predicting R0 using 4-fold cross validation. R0 was <5, 5–10 and >10 for 81, 14 and 3 settings respectively, with no apparent regional differences and in the limited available data, it was usually lower for rural than urban areas. R0 was most correlated with educational attainment, and household indicators for the Pearson and Spearman correlation coefficients respectively and with poverty-related indicators followed by the crude death rate considering the Maximum Information Coefficient, although the correlation for each was relatively weak (Pearson correlation coefficient: 0.4, 95%CI: (0.24,0.48) for educational attainment). A random forest did not perform better in predicting R0 than simple linear regression, depending on the subsets of training indicators and studies, and neither out-performed a regional averaging approach. R0 for rubella is typically low and using indicators to estimate its value is not straightforward. A regional averaging approach may provide as reliable an estimate of R0 for settings lacking seroprevalence data as one based on indicators. The findings may be relevant for other infections and studies estimating the disease burden and the impact of interventions for settings lacking seroprevalence data.

Modelling the early phase of the Belgian COVID-19 epidemic using a stochastic compartmental model and studying its implied future trajectories

Following the onset of the ongoing COVID-19 pandemic throughout the world, a large fraction of the global population is or has been under strict measures of physical distancing and quarantine, with many countries being in partial or full lockdown. These measures are imposed in order to reduce the spread of the disease and to lift the pressure on healthcare systems. Estimating the impact of such interventions as well as monitoring the gradual relaxing of these stringent measures is quintessential to understand how resurgence of the COVID-19 epidemic can be controlled for in the future. In this paper we use a stochastic age-structured discrete time compartmental model to describe the transmission of COVID-19 in Belgium. Our model explicitly accounts for age-structure by integrating data on social contacts to (i) assess the impact of the lockdown as implemented on March 13, 2020 on the number of new hospitalizations in Belgium; (ii) conduct a scenario analysis estimating the impact of possible exit strategies on potential future COVID-19 waves. More specifically, the aforementioned model is fitted to hospital admission data, data on the daily number of COVID-19 deaths and serial serological survey data informing the (sero)prevalence of the disease in the population while relying on a Bayesian MCMC approach. Our age-structured stochastic model describes the observed outbreak data well, both in terms of hospitalizations as well as COVID-19 related deaths in the Belgian population. Despite an extensive exploration of various projections for the future course of the epidemic, based on the impact of adherence to measures of physical distancing and a potential increase in contacts as a result of the relaxation of the stringent lockdown measures, a lot of uncertainty remains about the evolution of the epidemic in the next months.

Seasonal dynamics of varicella incidence in Bulgaria

Aim: Spatial and temporal distribution of varicella is heterogenic and insufficiently studied in Europe. The present study tries to fill the gap that exists about the seasonality of the infection in Bulgaria. Materials & methods: A 4-year retrospective study of the monthly and seasonal varicella epidemiology was performed at both national and district level. Results: In Bulgaria, varicella incidence peaked during winter (37% of the 2015–2018 cases), followed by spring (33%) and autumn (23%). Highly populated districts were more likely to follow this pattern, while less inhabited districts with smaller urbanized areas showed different periodicity of the infection. Conclusion: Winter peak in varicella incidence is positively associated with high accumulation of people in the large cities (>75,000 inhabitants).

Combining serological and contact data to derive target immunity levels for achieving and maintaining measles elimination

BACKGROUND

Vaccination has reduced the global incidence of measles to the lowest rates in history. However, local interruption of measles virus transmission requires sustained high levels of population immunity that can be challenging to achieve and maintain. The herd immunity threshold for measles is typically stipulated at 90-95%. This figure does not easily translate into age-specific immunity levels required to interrupt transmission. Previous estimates of such levels were based on speculative contact patterns based on historical data from high-income countries. The aim of this study was to determine age-specific immunity levels that would ensure elimination of measles when taking into account empirically observed contact patterns.

METHODS

We combined estimated immunity levels from serological data in 17 countries with studies of age-specific mixing patterns to derive contact-adjusted immunity levels. We then compared these to case data from the 10 years following the seroprevalence studies to establish a contact-adjusted immunity threshold for elimination. We lastly combined a range of hypothetical immunity profiles with contact data from a wide range of socioeconomic and demographic settings to determine whether they would be sufficient for elimination.

RESULTS

We found that contact-adjusted immunity levels were able to predict whether countries would experience outbreaks in the decade following the serological studies in about 70% of countries. The corresponding threshold level of contact-adjusted immunity was found to be 93%, corresponding to an average basic reproduction number of approximately 14. Testing different scenarios of immunity with this threshold level using contact studies from around the world, we found that 95% immunity would have to be achieved by the age of five and maintained across older age groups to guarantee elimination. This reflects a greater level of immunity required in 5-9-year-olds than established previously.

CONCLUSIONS

The immunity levels we found necessary for measles elimination are higher than previous guidance. The importance of achieving high immunity levels in 5-9-year-olds presents both a challenge and an opportunity. While such high levels can be difficult to achieve, school entry provides an opportunity to ensure sufficient vaccination coverage. Combined with observations of contact patterns, further national and sub-national serological studies could serve to highlight key gaps in immunity that need to be filled in order to achieve national and regional measles elimination.

Identifying human encounters that shape the transmission of Streptococcus pneumoniae and other acute respiratory infections

Although patterns of social contacts are believed to be an important determinant of infectious disease transmission, it remains unclear how the frequency and nature of human interactions shape an individual's risk of infection. We analysed data on daily social encounters individually matched to data on S. pneumoniae carriage and acute respiratory symptoms (ARS), from 566 individuals who took part in a survey in South-West Uganda. We found that the frequency of physical (i.e. skin-to-skin), long (≥1 h) and household contacts - which capture some measure of close (i.e. relatively intimate) contact - was higher among pneumococcal carriers than non-carriers, and among people with ARS compared to those without, irrespective of their age. With each additional physical encounter the age-adjusted risk of carriage and ARS increased by 6% (95%CI 2-9%) and 7% (2-13%) respectively. In contrast, the number of casual contacts (<5 min long) was not associated with either pneumococcal carriage or ARS. A detailed analysis by age of contacts showed that the number of close contacts with young children (<5 years) was particularly higher among older children and adult carriers than non-carriers, while the higher number of contacts among people suffering from ARS was more homogeneous across contacts of all ages. Our findings provide key evidence that the frequency of close interpersonal contact is important for transmission of respiratory infections, but not that of casual contacts. Those results are essential for both improving disease prevention and control efforts as well as informing research on infectious disease dynamics and transmission models, and more studies should be undertaken to further validate our results.

Modeling the impact of changes in day-care contact patterns on the dynamics of varicella transmission in France between 1991 and 2015

Annual incidence rates of varicella infection in the general population in France have been rather stable since 1991 when clinical surveillance started. Rates however show a statistically significant increase over time in children aged 0-3 years, and a decline in older individuals. A significant increase in day-care enrolment and structures' capacity in France was also observed in the last decade. In this work we investigate the potential interplay between an increase of contacts of young children possibly caused by earlier socialization in the community and varicella transmission dynamics. To this aim, we develop an age-structured mathematical model, informed with historical demographic data and contact matrix estimates in the country, accounting for longitudinal linear increase of early childhood contacts. While the reported overall varicella incidence is well reproduced independently of mixing variations, age-specific empirical trends are better captured by accounting for an increase in contacts among pre-school children in the last decades. We found that the varicella data are consistent with a 30% increase in the number of contacts at day-care facilities, which would imply a 50% growth in the contribution of 0-3y old children to overall yearly infections in 1991-2015. Our findings suggest that an earlier exposure to pathogens due to changes in day-care contact patterns, represents a plausible explanation for the epidemiological patterns observed in France. Obtained results suggest that considering temporal changes in social factors in addition to demographic ones is critical to correctly interpret varicella transmission dynamics.

The potential impact of the demographic transition in the Senegal-Gambia region of sub-Saharan Africa on the burden of infectious disease and its potential synergies with control programmes: the case of hepatitis B

BACKGROUND

Sub-Saharan Africa (SSA) continues to suffer high communicable disease burdens as its demographic transition (DT) proceeds. Although the consequent changes in population structures influence age-specific contact patterns relevant for transmission, the age distribution of immunity, and the disease burden, investigation of the potential of DT to affect infectious disease epidemiology in regions of SSA has hitherto been overlooked. With a substantial disease burden and complex epidemiology, hepatitis B virus (HBV) represents a prime example of an infection whose epidemiology may be significantly influenced by the DT.

METHODS

An age-structured mathematical model for HBV in the Senegal and Gambia (SG) region was set within a demographic framework with varying vital rates mirroring the entire course of the DT there over 1850-2100, to investigate the effects of the DT on HBV epidemiology, with and without the combined action of vaccination. The model was run from its reconstructed ancien régime (old order) demo-epidemiologic equilibrium and calibrated against SG 1950 age-distribution estimates and Gambian pre-vaccination HBV age-prevalence data.

RESULTS

The model, which reproduced well demographic and HBV age-prevalence data, predicted a complex transition of HBV epidemiology over the course of the DT. This included a prolonged epoch of expansion alongside population growth and rejuvenation until 1990-2000, followed by a dramatic retreat, mainly reflecting projected fertility decline during the twenty-first century. This transitional pattern was mostly explained by the underlying demographically driven changes in horizontal transmission resulting from the changes in the age structure of the population. During 2000-2150 the HBV burden is predicted to decline by more than 70% even in the absence of vaccination.

CONCLUSIONS

Demographic change alone may strongly affect HBV disease burden and shape HBV endemicity. The onset of the demographically driven decline in HBV prevalence, aligned with the expansion of HBV vaccination, forms a synergy potentially boosting effectiveness of control. Such a synergy currently appears to be presenting a "window of opportunity" facilitating HBV elimination which it would be important to exploit and which underlines the importance of taking demographic change into account when assessing the potential longer term impact of vaccination and other control measures.

Varicella zoster virus transmission dynamics in Vojvodina, Serbia

This study aimed at establishing baseline key epidemiological parameters for varicella zoster virus (VZV) infection in Vojvodina, Serbia, with the ultimate goal to quantify the VZV transmission potential in the population. Seroprevalence data generated during the first large cross-sectional VZV serosurvey were modelled, using a two-tiered modelling approach to calculate age-specific forces of infection (FOI), the basic reproduction number (R₀) and herd immunity threshold (H). Seroprevalence and modelling data were compared with corresponding pre-vaccination epidemiological parameters from 11 countries participating in the European Sero-Epidemiology Network 2 (ESEN2) project. Serbia fits into the general dynamic VZV transmission patterns in Europe in the pre-vaccine era, with estimated R₀ = 4.12, (95% CI: 2.69–7.07) and H = 0.76 (95% CI: 0.63–0.86). The highest VZV transmission occurs among preschool children, as evidenced by the estimation of the highest FOI (0.22, 95% CI: 0.11–0.34) in the 0.5–4 age group, with a peak FOI of 0.25 at 2.23 years. Seroprevalence was consistently lower in 5–14 year-olds, resulting in considerable shares of VZV-susceptible adolescents (7.3%), and young adults (6%), resembling the situation in a minority of European countries. The obtained key epidemiological parameters showed most intense VZV transmission in preschool children aged <4 years, justifying the consideration of universal childhood immunization in the future. National immunization strategy should consider programs for VZV serologic screening and immunization of susceptible groups, including adolescents and women of reproductive age. This work is an important milestone towards the evaluation of varicella immunization policy options in Serbia.

Varicella infection in a non-universally vaccinated population: Actual epidemiology in Bulgaria (2013-2015)

BACKGROUND

Varicella is a common and usually mild disease but it has great importance in regard to general infectious morbidity. The current study aimed to characterize possible risk factors of varicella epidemiology in Bulgaria, a country where infection follows its natural epidemiological pattern as no mandatory or recommended vaccine is currently applied.

METHODS

Administrative regions of Bulgaria were used as units of observation and a set of sociodemographic and economic determinants, as well as geographic location (south or north) were tested for associations with the mean 3-year varicella incidence rates (2013-2015).

RESULTS

The proportion of urban population, proportion of females, number of health care units and proportion of urban population aged <10 years were the four sociodemographic variables most strongly and significantly correlated (p<0.05) with varicella frequency (Spearman's rank correlation coefficients of 0.62, 0.47, 0.43, and 0.38, respectively). After reducing the number of intercorrelated factors to a few principal components and accounting for confounders, the demographic component and geographic location remained most robustly associated with varicella incidence in Bulgaria (adjusted R² of 0.51, p<0.001).

CONCLUSIONS

The results obtained identify important determinants in the local epidemiology of varicella and show that community characteristics should be considered, to improve our understanding of varicella distribution.

A systematic review of varicella seroprevalence in European countries before universal childhood immunization: deriving incidence from seroprevalence data

Surveillance systems for varicella in Europe are highly heterogeneous or completely absent. We estimated the varicella incidence based on seroprevalence data, as these data are largely available and not biased by under-reporting or underascertainment. We conducted a systematic literature search for varicella serological data in Europe prior to introduction of universal varicella immunization. Age-specific serological data were pooled by country and serological profiles estimated using the catalytic model with piecewise constant force of infection. From the estimated profiles, we derived the annual incidence of varicella infection (/100·000) for six age groups (<5, 5–9, 10–14, 15–19, 20–39 and 40–65 years). In total, 43 studies from 16 countries were identified. By the age of 15 years, over 90% of the population has been infected by varicella in all countries except for Greece (86·6%) and Italy (85·3%). Substantial variability across countries exists in the age-specific annual incidence of varicella primary infection among the <5 years old (from 7052 to 16 122 per 100 000) and 5–9 years old (from 3292 to 11 798 per 100 000). The apparent validity and robustness of our estimates highlight the importance of serological data for the characterization of varicella epidemiology, even in the absence of sampling or assay standardization.

Economic burden of varicella in children 1-12 Years of age in Hungary, 2011-2015

Background

Although live-attenuated varicella-zoster virus (VZV) vaccines have been proven to be safe and effective in preventing varicella and real-word evidence shows routine childhood immunization programs are effective in dramatically reducing varicella associated morbidity and mortality, varicella vaccine is not included in the National Immunization Program (NIP) in Hungary. The purpose of this study was to evaluate the clinical and economic burden associated with varicella in Hungary.

Methods

This was a multicenter, retrospective, chart review study of patients aged 1–12 years with a primary varicella diagnosis between 2011 and 2015. Healthcare resource utilization (HCRU) associated with varicella, unit costs, and work loss were used to estimate direct and indirect costs. All costs are presented in 2015 HUF / Euros (€).

Results

156 children with varicella were included (75 outpatients, 81 inpatients), with a mean age of 4.4 (SD: 2.0) and 3.7 (SD: 2.1) years, respectively. One or more complications were reported by 12.0% of outpatients and 92.6% of inpatients, the most common being dehydration, skin and soft tissue infections, pneumonia, keratoconjunctivitis, and cerebellitis. HCRU estimates included use of over-the-counter (OTC) medications (96.0% outpatients, 53.1% inpatients), prescription medications (9.3% outpatients, 70.4% inpatients), tests/procedures (4.0% outpatients, 97.5% inpatients), and consultation with allied health professionals (2.7% outpatients, 30.9% inpatients). The average duration of hospital stay (inpatients) was 3.6 (95% CI: 3.2, 4.1) days. The total combined direct and indirect cost per varicella case was 228,146.7 Hungarian Forint (HUF)/€ 736.0 for inpatients and 49,790.6 HUF/€ 106.6 for outpatients. The overall annual cost of varicella in Hungary for children aged <15 years in 2015 was estimated at 1,903,332,524.3 HUF/ € 6,139,980.4.

Conclusion

Varicella is associated with substantial clinical burden in Hungary, resulting in the utilization of a significant amount of healthcare resources. These results support the need for routine vaccination of all healthy children to reduce the varicella-associated disease burden.

Estimation of the burden of varicella in Europe before the introduction of universal childhood immunization

BACKGROUND

Varicella is generally considered a mild disease. Disease burden is not well known and country-level estimation is challenging. As varicella disease is not notifiable, notification criteria and rates vary between countries. In general, existing surveillance systems do not capture cases that do not seek medical care, and most are affected by underreporting and underascertainment. We aimed to estimate the overall varicella disease burden in Europe to provide critical information to support decision-making regarding varicella vaccination.

METHODS

We conducted a systematic literature review to identify all available epidemiological data on varicella IgG antibody seroprevalence, primary care and hospitalisation incidence, and mortality. We then developed methods to estimate age-specific varicella incidence and annual number of cases by different levels of severity (cases in the community, health care seekers in primary care and hospitals, and deaths) for all countries belonging to the European Medicines Agency (EMA) region and Switzerland.

RESULTS

In the absence of universal varicella immunization, the burden of varicella would be substantial with a total of 5.5 million (95% CI: 4.7-6.4) varicella cases occurring annually across Europe. Variation exists between countries but overall the majority of cases (3 million; 95% CI: 2.7-3.3) would occur in children <5 years. Annually, 3-3.9 million patients would consult a primary care physician, 18,200-23,500 patients would be hospitalised, and 80 varicella-related deaths would occur (95% CI: 19-822).

CONCLUSIONS

Varicella disease burden is substantial. Most cases occur in children <5 years old but adults require hospitalisation more often and are at higher risk of death. This information should be considered when planning and evaluating varicella control strategies. A better understanding of the driving factors of country-specific differences in varicella transmission and health care utilization is needed. Improving and standardizing varicella surveillance in Europe, as initiated by the European Centre for Disease Prevention and Control (ECDC), is important to improve data quality to facilitate inter-country comparison.

The natural history of varicella zoster virus infection in Norway: Further insights on exogenous boosting and progressive immunity to herpes zoster

We use age-structured models for VZV transmission and reactivation to reconstruct the natural history of VZV in Norway based on available pre-vaccination serological data, contact matrices, and herpes zoster incidence data. Depending on the hypotheses on contact and transmission patterns, the basic reproduction number of varicella in Norway ranges between 3.7 and 5.0, implying a vaccine coverage between 73 and 80% to effectively interrupt transmission with a 100% vaccine efficacy against infection. The varicella force of infection peaks during early childhood (3–5 yrs) and shows a prolonged phase of higher risk during the childbearing period, though quantitative variations can occur depending on contact patterns. By expressing the magnitude of exogenous boosting as a proportion of the force of infection, it is shown that reactivation is well described by a progressive immunity mechanism sustained by a large, though possibly below 100%, degree of exogenous boosting, in agreement with findings from other Nordic countries, implying large reproduction numbers of boosting. Moreover, magnitudes of exogenous boosting below 40% are robustly disconfirmed by data. These results bring further insight on the magnitude of immunity boosting and its relationship with reactivation.

Squamous epitheliotropism of Enterovirus A71 in human epidermis and oral mucosa

Hand-foot-and-mouth disease is a self-limiting paediatric infectious disease commonly caused by Enterovirus A71 (Genus: Enterovirus, Family: Picornaviridae). Typical lesions in and around the hands, feet, oral cavity and other places may rarely be complicated by acute flaccid paralysis and acute encephalomyelitis. Although virus is readily cultured from skin vesicles and oral secretions, the cellular target/s of Enterovirus A71 in human skin and oral mucosa are unknown. In Enterovirus A71-infected human skin and oral mucosa organotypic cultures derived from the prepuce and lip biopsies, focal viral antigens and viral RNA were localized to cytoplasm of epidermal and mucosal squamous cells as early as 2 days post-infection. Viral antigens/RNA were associated with cytoplasmic vacuolation and cellular necrosis. Infected primary prepuce epidermal keratinocyte cultures showed cytopathic effects with concomitant detection of viral antigens from 2 days post-infection. Supernatant and/or tissue homogenates from prepuce skin organotypic cultures and primary prepuce keratinocyte cultures showed viral titres consistent with active viral replication. Our data strongly support Enterovirus A71 squamous epitheliotropism in the human epidermis and oral mucosa, and suggest that these organs are important primary and/or secondary viral replication sites that contribute significantly to oral and cutaneous viral shedding resulting in person-to-person transmission, and viraemia, which could lead to neuroinvasion.

Estimating Age-Specific Immunity and Force of Infection of Varicella Zoster Virus in Norway Using Mixture Models

This study applies mixture modelling to examine age-specific immunity to varicella zoster virus (VZV) infection in Norway based on the first large-scale serological study in the general population. We estimated the seropositive proportions at different ages and calculated the underlying force of infection by using a sample of 2103 residual sera obtained from patients seeking primary and hospital care. A rapid increase in the VZV-associated immunity is observed in the first years of life with 63% of children being immune by age 5. The increase in the immunity levels slows down thereafter, with a large proportion of adults still susceptible by age 20 (around 14.5%), thus at risk of serious sequelae of varicella infection. The corresponding force of infection peaks during the preschool period, subsequently declines to a minimum between ages 10 and 20 years, and afterwards moderately increases to reach a plateau lasting throughout the childbearing period. In comparison with the traditional cut-off approach, mixture modelling used the whole data without producing any inconclusive cases, led to an unbiased classification of individuals between susceptible and immune, and provided a smoother immune profile by age. These findings represent an important step towards any decision about the introduction of varicella vaccination in Norway, as they are a primary input for mathematical transmission models aimed at evaluating potential vaccination scenarios.

Clustering of contacts relevant to the spread of infectious disease

OBJECTIVE

Infectious disease spread depends on contact rates between infectious and susceptible individuals. Transmission models are commonly informed using empirically collected contact data, but the relevance of different contact types to transmission is still not well understood. Some studies select contacts based on a single characteristic such as proximity (physical/non-physical), location, duration or frequency. This study aimed to explore whether clusters of contacts similar to each other across multiple characteristics could better explain disease transmission.

METHODS

Individual contact data from the POLYMOD survey in Poland, Great Britain, Belgium, Finland and Italy were grouped into clusters by the k medoids clustering algorithm with a Manhattan distance metric to stratify contacts using all four characteristics. Contact clusters were then used to fit a transmission model to sero-epidemiological data for varicella-zoster virus (VZV) in each country.

RESULTS AND DISCUSSION

Across the five countries, 9-15 clusters were found to optimise both quality of clustering (measured using average silhouette width) and quality of fit (measured using several information criteria). Of these, 2-3 clusters were most relevant to VZV transmission, characterised by (i) 1-2 clusters of age-assortative contacts in schools, (ii) a cluster of less age-assortative contacts in non-school settings. Quality of fit was similar to using contacts stratified by a single characteristic, providing validation that single stratifications are appropriate. However, using clustering to stratify contacts using multiple characteristics provided insight into the structures underlying infection transmission, particularly the role of age-assortative contacts, involving school age children, for VZV transmission between households.