key: cord-0912435-vzpthw4j
authors: Radon, Katja; Bakuli, Abhishek; Pütz, Peter; Le Gleut, Ronan; Guggenbuehl Noller, Jessica Michelle; Olbrich, Laura; Saathoff, Elmar; Garí, Mercè; Schälte, Yannik; Frahnow, Turid; Wölfel, Roman; Pritsch, Michael; Rothe, Camilla; Pletschette, Michel; Rubio-Acero, Raquel; Beyerl, Jessica; Metaxa, Dafni; Forster, Felix; Thiel, Verena; Castelletti, Noemi; Rieß, Friedrich; Diefenbach, Maximilian N.; Fröschl, Günter; Bruger, Jan; Winter, Simon; Frese, Jonathan; Puchinger, Kerstin; Brand, Isabel; Kroidl, Inge; Wieser, Andreas; Hoelscher, Michael; Hasenauer, Jan; Fuchs, Christiane
title: From first to second wave: follow-up of the prospective COVID-19 cohort (KoCo19) in Munich (Germany)
date: 2021-09-08
journal: BMC Infect Dis
DOI: 10.1186/s12879-021-06589-4
sha: f69af80739e731764204be216a4dffa8bf4b79f7
doc_id: 912435
cord_uid: vzpthw4j

BACKGROUND: In the 2nd year of the COVID-19 pandemic, knowledge about the dynamics of the infection in the general population is still limited. Such information is essential for health planners, as many of those infected show no or only mild symptoms and thus, escape the surveillance system. We therefore aimed to describe the course of the pandemic in the Munich general population living in private households from April 2020 to January 2021. METHODS: The KoCo19 baseline study took place from April to June 2020 including 5313 participants (age 14 years and above). From November 2020 to January 2021, we could again measure SARS-CoV-2 antibody status in 4433 of the baseline participants (response 83%). Participants were offered a self-sampling kit to take a capillary blood sample (dry blood spot; DBS). Blood was analysed using the Elecsys® Anti-SARS-CoV-2 assay (Roche). Questionnaire information on socio-demographics and potential risk factors assessed at baseline was available for all participants. In addition, follow-up information on health-risk taking behaviour and number of personal contacts outside the household (N = 2768) as well as leisure time activities (N = 1263) were collected in summer 2020. RESULTS: Weighted and adjusted (for specificity and sensitivity) SARS-CoV-2 sero-prevalence at follow-up was 3.6% (95% CI 2.9–4.3%) as compared to 1.8% (95% CI 1.3–3.4%) at baseline. 91% of those tested positive at baseline were also antibody-positive at follow-up. While sero-prevalence increased from early November 2020 to January 2021, no indication of geospatial clustering across the city of Munich was found, although cases clustered within households. Taking baseline result and time to follow-up into account, men and participants in the age group 20–34 years were at the highest risk of sero-positivity. In the sensitivity analyses, differences in health-risk taking behaviour, number of personal contacts and leisure time activities partly explained these differences. CONCLUSION: The number of citizens in Munich with SARS-CoV-2 antibodies was still below 5% during the 2nd wave of the pandemic. Antibodies remained present in the majority of SARS-CoV-2 sero-positive baseline participants. Besides age and sex, potentially confounded by differences in behaviour, no major risk factors could be identified. Non-pharmaceutical public health measures are thus still important. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s12879-021-06589-4.

Background The SARS-CoV-2 virus affected almost all nations within a few weeks. Given the nature of the virus, a large proportion of infected individuals present only mild symptoms or no symptoms at all. Therefore, population-based sero-prevalence studies are necessary to estimate the true prevalence of the infection in the population. Starting in March 2020, such sero-prevalence studies have been conducted in many countries, mostly during or after the first wave of the pandemic [1] . Depending on the serological test used, the type of sample drawn, the timing of the study, and the region, general population sero-prevalence ranged from < 0.1% in Brazil to well over 20% in the USA [2] . For the German context, we reported a sero-prevalence of 1.8% in Munich, sampled towards the end of the first wave in Germany [3] .

Following the introduction of public health measures (lock-down including school closures) in March 2020 in Germany, the first wave of the pandemic was perceived as relatively mild with around 6000 cases registered in Munich during this period (Munich population ~ 1.5 Mio). Between June and October, public health measures were reduced, although physical distancing of 1.5 m between two persons, avoidance of mass events, and obligatory use of face masks, e.g. in restaurants and shops, were still required. Subsequently, officially registered monthly case numbers in Munich rose from 389 in June to 7181 in October 2020. A partial national lockdown was implemented on November 2nd, 2020. After a further rise in officially registered case numbers and COVID-19 related deaths, national lock-down measures were increased from December 16th, 2020 on, including closure of schools, shops (other than grocery and drug stores), restaurants, and hotels.

Given that asymptomatic and mildly symptomatic cases escape surveillance systems, prospective population-based cohort studies offer the chance to better understand the course of disease in the general population. They are independent of testing strategies and help to identify the population at risk over time. In addition, they provide an indication of population groups less well protected by public health measures. We therefore followed up the participants of the Munich COVID-19 cohort (KoCo19) to explore the SARS-CoV-2 antibody prevalence in the Munich general population at two time points: at the time the acute outbreaks happened and seven months later. In addition, we aimed at the identification of risk factors (demographic, social-economic, health status or individual risk behaviours factors) for acquiring SARS-CoV-2 infection defined by serology. The baseline study took place from April to June 2020, the questionnaire follow-up in summer 2020 and the 1st antibody follow-up was realised from early November 2020 to January 2021. On December 1st 2020 the KoCo19 cohort joined the ORCHESTRA (Connecting European Cohorts to Increase Common and Effective Response to SARS-CoV-2 Pandemic) project.

We described the baseline study in detail in [4] . In short, a random sample of the Munich population living in private households was drawn by random walk method. All household members older than 13 years were invited to provide a serum sample and to answer an online questionnaire. Serum samples were analysed for SARS-CoV-2 antibodies using the Elecsys ® Anti-SARS-CoV-2 (Roche) test [5] . Field work for the baseline study took place between April 5th and June 12th, 2020.

An online questionnaire covering risk behaviour, health related items, and psychosocial aspects (hereafter "behaviour questionnaire") was offered from June 4th to October 31st, 2020 to all 5240 participants who did not withdraw from the study. In parallel, an online-questionnaire on leisure time behaviour was available (hereafter "leisure time questionnaire"). We split the questionnaire into two, because long questionnaires are less likely to be completed [6] . Participants recruited in April (May to June) 2020 received an invitation via e-mail on June 4th (June 25th) with subsequent reminders and telephone follow-ups. In total, 3400 participants completed the behaviour questionnaire and 1390 participants the leisure time questionnaire.

On November 2nd 2020, we started the 1st antibody follow-up by sending out boxes with a self-sampling kit to take a capillary blood samples (dry blood spot; DBS) to the 5292 participants (2978 households) of the baseline study. Between baseline and follow-up, 77 participants withdrew from the study and were thus not contacted for the follow-up. Instructions for selfsampling were provided, including a video tutorial (https:// www. youtu be. com/ watch?v= vpZUz uQV10 E& featu re= emb_ title). Samples were collected using a barcode-labelled neonatal screening filter card (Euroimmun ZV 9701-0101) with circles indicating where the blood should be collected. Afterwards, participants should dry the filter card at least 12 h at room temperature, pack them in the sealable plastic pouch, place the plastic pouch into the prepaid envelope, and ship the envelope by mail to the laboratory. In case of handling difficulties, our telephone and e-mail hotline were available for any questions.

From November 2nd to January 31st, 2021, we received 4444 DBS samples from 2571 households (individual response 84%, household response 86%). Roughly half of the DBS samples (2372 of 4433; 54%) arrived at our laboratory within 1 week of mailing (November 2nd to November 8th). By week 2 (November 9th to November 15th), more than three quarters were received (3369 of 4433; 76%). Most of the remaining samples were turned in between week 3 (N = 372 from November 16th to November 22th) and week 4 (N = 343; November 23rd to November 29th). Few samples were received in December 2020 and January 2021 (N = 326; 7%). Participants not being able to collect a DBS on their own (N = 29) and those with intermediate results (N = 34, s. laboratory methods) were offered a full-blood test at our centre. For the latter group, this served to clarify the DBS result. However, 11 of the 34 participants with intermediate results in the DBS did not show up at our centre and thus had to be excluded from analyses, leaving 4433 subjects with baseline questionnaire, baseline serology and follow-up DBS data for the main analyses (Fig. 1 ).

The following items were considered for the analyses presented in this paper:

Baseline individual questionnaire: Follow-up questionnaire:

• Self-estimated health-related risk taking behaviour (10-level Likert scale from "not at all risk tolerant" to "very risk tolerant"): Dichotomised into not high (≤ 5, Quartile 3) and high self-estimated healthrelated risk taking behaviour (> 5) • Personal contacts: Five questions on places of personal contacts outside the own household during the two weeks before answering the questionnaire (meeting people, grocery shopping, shopping, use of public transport, work outside home), each assessed on a 5-level Likert scale: not at all (= 1); once per week (= 2); 2-4 times per week (= 3); 5 times per week (= 4); more often (= 5). Places of personal contacts were multiplied by frequency of contacts (0 contacts (= 0), 1 contact (= 1), 2-4 contacts (= 2) and 5 + contacts (= 3)) and summed up, resulting in a score ranging from 0 to 25. The score was dichotomised into lower number of personal contacts (≤ 8, Median) and higher number of personal contacts (> 8). The score was dichotomised into non-high leisure time activities (≤ 11, Quartile 3) and high leisure time activities (> 11).

• Number and intensity of leisure time activities before the pandemic (in February 2020): For that time, 16 activities assessed on a 5-level Likert scale from "never" (= 0) to "very often" (= 4): visit family member; visit friends; going out with friends; attend a party, festival, bar, pub or disco; go to the cinema; attend a theatre, opera or ballet performance; work out in a gym; visit a swimming pool; visit a sauna; skiing; train for a team sport or take part in sporting competitions; watch a sports game or event live outdoors; watch a sports game or event live indoors; worship attendance; play an instrument in an orchestra; sing in a choir. Activities were multi- Fig. 1 Flow chart of obtaining the study population. Light blue boxes indicate the total number of participants and households in the baseline study, the light green box refers to subjects included in the sero-prevalence analyses of the data, dark green boxes indicate analyses restricted to sero-incidence data for subjects with complete leisure time information (left) and risk behaviour data (right) Radon et al. BMC Infect Dis (2021) 21:925 plied by the Likert scores and summed up resulting in a score from 0 to 64. • Number and intensity of leisure time activities two weeks prior to the follow-up questionnaire: The score for leisure time activities at follow-up was built the same way as the score for leisure time activities before the pandemic. However, the number of leisure time activities was only seven at that time as many activities were not possible due to the restrictions related to the pandemic: visit family member; visit friends; going out with friends; visit a swimming pool; worship attendance; play an instrument in an orchestra; sing in a choir. Therefore, the resulting score only ranged from 0 to 28. The score was dichotomised into non-high leisure time activities at follow-up (≤ 5, Quartile 3) and high leisure time activities (> 5).

Filter paper cards were further processed if at least two of the five circles on the card were completely soaked with blood. Valid samples were stored at 4 °C until analysis. Before analysis, filter paper cards were equilibrated to room temperature and three blood-soaked smaller circles (diameter 3.2 mm) of each filter paper card were automatically punched into a 96-wells plate (Panthera-Puncher ™ 9, PerkinElmer). After elution, samples were transferred to a Cobas e801 module (Roche) compatible sample micro cup (Roche, 05085713001) for analysis using the Elecsys ® Anti-SARS-CoV-2 assay (Roche). Based on our validation study, DBS samples were considered positive if SARS-CoV-2 antibody levels were ≥ 0.12. Samples with SARS-CoV-2 antibody levels in the range between 0.09 and 0.12 were considered intermediate, and subsequently confirmed by plasma samples (s. Study population and field work). All other samples were considered negative. Compared to full blood samples, sensitivity of the DBS method was 99.2% and specificity 98.7%. Details of the laboratory methods are described in [7] .

All statistical analyses were performed using the statistical software R (version 4.0.3, R Development Core Team, 2020). The SARS-CoV-2 sero-prevalence was estimated primarily based on the DBS test results of the study participants applying the classification as described above (Laboratory methods). If the DBS test yielded an intermediate result, we considered the result of the full blood sampling. As described in [5] , an optimised cut-off of 0.4218 for the full blood sampling was used to predict SARS-CoV-2 sero-positivity with an estimated specificity and sensitivity of 99.7% and 88.6%, respectively (with regard to PCR test results considered as ground truth). We used these estimates to adjust the prevalence for the imperfect test performance [8] . The specificity and sensitivity of DBS with regard to full blood samples being very high, additional adjustment was omitted (Additional file 1: Appendix Text and Table S1 ).

The prevalence (adjusted or unadjusted for the specificity and the sensitivity of the test) was calculated in two different ways: including the information from the sampling design of the cohort [3] via the use of a weighting scheme, or without it. To account for the sampling design, the sampling weights computed at baseline (inverse of the probability of each individual to be included in the sample) were used for the follow-up analysis. These sampling weights were corrected for the attrition observed at follow-up by modelling the underlying non-response mechanism and estimating probabilities of response for each unit. Ten response homogeneity groups (where we assumed the non-response to be completely at random, [9] ) were created using the deciles of the estimated probabilities of response. These weights adjusted for the non-response were calibrated [10] on updated information from the Munich population (at 31.12.2020) in order to mirror the age, sex, country of birth, presence of children in the household and single member household structures. Moreover, to correct the sample for the loss of positive cases at follow-up, the sampling weights were calibrated on the estimated number of positive cases at baseline. Weighted prevalence estimates were calculated using these calibrated weights, and the associated 95% confidence intervals were computed based on variance estimators based on linearization [10] and residual [10, 11] techniques. These variance estimates were computed in order to account for every step in the selection process of the units, i.e., V = V1 + V2 with V1 the variance due to the sampling design and V2 the one due to the non-response [12] . For unweighted prevalence estimates, confidence intervals were determined by using a nonparametric cluster bootstrap procedure that accounts for household clustering [13] . To that end, 5000 bootstrap datasets were generated each by sampling n h households with replacement from the original sample of n h households. The sero-prevalence was estimated in each bootstrap sample and the 2.5 and 97.5 percentiles of the resulting 5000 estimates defined the 95% confidence intervals.

To analyse spatial clustering, we considered the mean within-cluster variance of the binary test results, with cluster variables being households, buildings, and geospatial clusters of different sizes. We performed a nonparametric approximate permutation test with 10,000 random permutations of cluster assignments. To account for household clustering, only full households were permuted when considering buildings and geospatial clusters [14] . In addition to this, we performed borough level sero-prevalence mapping using Conditional Auto Regressive Models which account for the spatial autocorrelation among neighbouring boroughs by using random effects. This allowed us to investigate if sero-prevalence was associated with the population density or not, as well as obtaining Borough/District level estimates within the city of Munich (Additional file 1: Appendix 2) [15] [16] [17] [18] .

We used generalised linear mixed models (GLMMs using the logit link function) to analyse the association between potential risk factors and SARS-CoV-2 sero-positivity at 1st follow-up, with a random effect for households to account for within household clustering of the data. Odds Ratio estimates and the corresponding confidence intervals were obtained applying a Bayesian framework with uniform priors on the regression estimates using the brms (Bayesian Regression Models using 'Stan') package in R [19, 20] . To account for missing data in covariates, we used the Joint Analysis and Imputation of Incomplete Data Framework (JointAI) in R for sensitivity analyses [21, 22] . In these sensitivity analyses, broad normal priors with mean zero and standard deviation 100 were used. The regression estimates were adjusted for the SARS-CoV-2 serology results at baseline, the time elapsed since baseline visit, age, and sex of the individual. Essentially, this adjustment for baseline positivity allowed us to obtain risk factors associated with newly incident cases within our cohort over and above the baseline positives.

To explore the importance of behavioural factors and leisure-time activities for the incidence of infection between baseline and follow-up, we used data of the 1st questionnaire follow-up combined with the DBS results. For these analyses, we included information of 2768 participants who responded to the behaviour questionnaire and had serology results; for the leisure time activities, we had questionnaire information for 1263 persons with serological results. Due to the large proportion of missing questionnaire data, we restricted these analyses to complete data and aggregated at the levels of the outcome variables. We analysed the incidence of new SARS-CoV-2 infections (as binomial outcome for proportions) between baseline and follow-up stratified for risk behaviour, leisure time activities, sex and age, using the count of new positives among the observed. Similar models were also applied to evaluate the association of the population densities at the constituency level and the trend in sero-prevalence estimates using aggregated data.

Follow-up participants compared to participants lost to follow-up were more likely to be between 35 and 79 years old, born in Germany (84% vs. 74%), and to have a higher socio-economic status ( Table 1 ). The latter was indicated, e.g., by level of education, household income, living area, and type of building. In addition, the sero-prevalence of SARS-CoV-2 at baseline was lower among follow-up participants (1.6%) compared to baseline only participants (2.6%; Table 1 ). These losses of positive cases at follow-up led to an underestimation of the total number of people tested positive at baseline in Munich (22,064 vs. 25,900 using all participants at baseline). To correct for this attrition bias, the weights at follow-up were calibrated on the estimate of positive cases at baseline, in addition to the other margins used to mirror the Munich structure.

The overall weighted and adjusted (for specificity and sensitivity) SARS-CoV-2 sero-prevalence at followup was 3.6% (95% CI 2.9-4.3%; Fig. 2 ). The overall unweighted and adjusted sero-prevalence was 3.1% (95% CI 2.5-3.8%), increasing from 2.5% (95% CI 1.7-3.3%) in the first week of November to 4.0% (95% CI 1.6-6.8%) in the last week of November (Fig. 3) . About half of the participants with intermediate result in the DBS test had a positive test result when considering the plasma sample. As plasma samples were collected in December and January, the prevalence estimates in the latest weeks were artificially high. Yet, the overall upward trend remained after excluding the participants with intermediate DBS result (Additional file 1: Figures S1, S2) .

Most participants who were SARS-CoV-2 sero-positive at baseline continued to be sero-positive at follow-up (64 out of 70 sero-positive subjects at baseline; Additional file 1: Table S3 ). The weighted and adjusted SARS-CoV-2 sero-incidence (negative at baseline, positive at followup) was estimated at 1.6% (95% CI 0.95-2.3%) (Fig. 2) .

Looking at the geospatial distribution of SARS-CoV-2 sero-prevalence by Munich city boroughs (Fig. 4) , an increase was visible from the South-East to the North-West of the city, although these differences were rather small. The estimates of Moran's I for spatial autocorrelation was 0.015 using the continuous distance based spatial neighbourhood matrix resulting in a p-value of 0.304, and hence was not statistically significant. Using a binary spatial neighbourhood matrix, the estimated Moran's I was 0.025 with a corresponding p-value of 0.269-thereby the conclusion remained unchanged. We also took population density as a potential risk factor into account and could not find any statistically significant association between population density in the constituency and SARS-CoV-2 sero-prevalence (Additional file 1: Figure S3 ).

The distribution of SARS-CoV-2 sero-positivity by covariates is shown in Additional file 1: Table S3 . Taking household clustering, time elapsed between baseline and follow-up, and baseline result into account, men had statistically significantly higher odds of sero-positivity at follow-up (OR adjusted for age: 2.4; 95% CI 1.0-6.0; Fig. 5 ). In addition, SARS-CoV-2 sero-prevalence decreased with increasing age group. It was lower in participants living in small apartment houses compared to participants living in single houses (OR adjusted for age and sex 0.0002; 95% CI 0.0-0.14).

The analysis of potential household and neighbourhood clustering indicated a highly significant within-household clustering of SARS-CoV-2 cases. In contrast, no indications for neighbourhood transmission of SARS-CoV-2 were observed (Fig. 6 , Additional file 1: Figure S5 ).

In order to identify which behavioural factors might be related to the different sero-prevalences in men and younger subjects, we compared these factors by sex and age group. The sum of contacts decreased significantly from the two younger age groups to the oldest age group in men and women (p < 0.001; Additional file 1: Table S4 ). In contrast, the self-estimated health-related risk taking behaviour and leisure time activity level were highest in Fig. 2 Weighted sero-prevalence and sero-incidence in % at follow-up adjusted (orange) and unadjusted (blue) for test specificity and sensitivity. The unadjusted weighted sero-prevalence was 3.4%, the relative number of new cases between baseline and follow-up 1.7%. Adjustment only slightly changed the unadjusted results

For sensitivity and specificity adjusted (left) and unadjusted (right) SARS-CoV-2 sero-prevalence over the follow-up period. The 95% confidence intervals for the weekly sero-prevalence were based on the 2.5 and 97.5 percentiles from 5000 repetitions of a cluster bootstrap that accounts for within household clustering. The estimates do not account for sample weights. The estimation without accounting for within-household clustering but considering sample weights produced similar trends (Additional file 1: Figure S1 ). A slight increase of sero-prevalence is indicated from the first to the fourth week of follow-up. The huge increase from week 4 to weeks 5-11 has to be taken with caution, as during these weeks, participants with intermediate results in the DBS (of whom 50% turned out to be positive in plasma sampling) were be retested by plasma-sampling during this time interval. The upward trend without these participants is shown in Additional file 1: Figures S1 and S2 the age group < 35 years. Comparing men and women, self-estimated health-related risk taking behaviour was statistically higher for men than for women in the age group 35-65 years only (p < 0.001, Additional file 1: Table S5 ). In contrast, sum of contacts and number of leisure time activities was similar for men and women by age strata (p > = 0.05).

In order to check for effect modification by sex, age and behavioural factors, we calculated the sex-and agestratified SARS-CoV-2 sero-incidence between baseline and follow-up by self-estimated health-related risk taking behaviour (Fig. 7A) , number of contacts outside own household (Fig. 7B) , number of leisure time activities in summer 2020 (Fig. 7C) . These data indicate a slightly higher risk of infection among men and women above the age of 34 years who indicated to have a high healthrelated risk taking behaviour compared to those with no high health-related risk taking behaviour. Men and women above the age of 64 years showed a higher SARS-CoV-2 sero-incidence if they had more personal contacts compared to participants having fewer personal contacts. Men with more leisure time activities in summer 2020 had a higher SARS-CoV-2 sero-incidence compared to less active men. However, all confidence intervals largely overlapped.

Our data indicate a low SARS-CoV-2 sero-prevalence for the Munich general population living in private households eight months after the start of the pandemic. The incidence between the end of the first wave and the middle of the second wave was about as high as the related number of infections acquired during the first wave. Almost all sero-positive participants at baseline remained sero-positive at follow-up, indicating a high validity of the antibody test. Additionally, this supports previous reports suggesting that the humoral SARS-CoV-2 immune response is stable at least over the first eight months after infection [23, 24] . We also showed a predominance in SARS-CoV-2 sero-positivity among male C Lower 95% confidence bounds of the weighted SARS-CoV-2 sero-prevalence; D Upper 95% confidence bounds of the weighted SARS-CoV-2 sero-prevalence. The sero-positivity varied slightly across the boroughs (as indicated by different colours), however, differences did not reach statistical significance compared to female participants, and a reduced antibody prevalence with increasing age group.

Based on our data, the sero-prevalence for the Munich population above the age of 13 years living in private households was 3.6% (95% CI 2.9-4.3%). Until the end of November 2020, a total of 30,180 SARS-CoV-2 cases were officially registered in Munich (https:// www. muenc hen. de/ ratha us/ Stadt infos/ Coron avirus-Fallz ahlen. html# Fallz ahlen; Access date: 19-April-2021) which results in a population prevalence of 1.9%. This prevalence increased to 44,377 registered cases by the end of December 2020 (population prevalence 3.0%). The data are not directly comparable, as the official data also include children and persons living in institutions. While the prevalence of infection in children was at that time considered to be smaller than in adults, it was unknown for people living in institutions (e.g., homes for the elderly). Nevertheless, the comparison gives an indication that the percentage of officially registered infections improved considerably compared to the beginning of the pandemic.

In a previous publication we estimated that solely one out of four infections was registered by the official infectious diseases surveillance system [3] . A few population-based SARS-CoV-2 sero-studies have been conducted since the beginning of the pandemic (for review see [2] ), most of them reporting sero-prevalences during or after the first wave. Up to now, only the Spanish national study reported the results of their follow-up data [25] with a sero-prevalence at follow-up (November 2020) of 5%, and thus comparable to our results.

In our study, one predictor of change in sero-prevalence from baseline to follow-up was male sex. While a higher risk of more severe COVID-19 among men was confirmed in several studies [26] , findings on sexdifferences in sero-prevalence are still inconsistent [2] . As younger age was also related to a larger increase in SARS-CoV-2 sero-prevalence at follow-up, one might assume that differences in behaviour may contribute to these findings. We could confirm differences in health-risk taking behaviour, frequency of leisure-time Association between potential risk factors and SARS-CoV-2 sero-positivity taking into account time between 1st and 2nd sampling, baseline result, age and sex. Age and sex was also adjusted for time between 1st and 2nd sampling, baseline result and each other (sex for age and age for sex). Unimputed (blue) and imputed (orange) GLM Models (Bayesian analysis). Main individual level risk factors were time and sex. Odds decreased by age and was lower for participants living in buildings with 3-4 apartments. Changes in the estimates by imputation were small activities, and number of contacts outside the own household especially by age. In the stratified analyses of the incidence of infection by age and sex, we observed a tendency that behaviour is related to higher seroincidence of infection; although the low incidence and the reduced number of respondents to the questionnaires limited the statistical power of these analyses and result interpretation has to be done cautiously. However, the hypothesis that specific behaviour, i.e., restriction of contacts, might reduce the risk of infection is also supported by our observation that patients with autoimmune disease were at reduced risk of SARS-CoV-2 sero-positivity. This finding is in line with studies among, e.g., patients with inflammatory bowel disease [27] . Overall, our findings support the notion that behavioural factors contribute to the spread of the pandemic, and therefore actions to increase adherence to public-health measures (such as information campaigns) are crucial especially in a time when acceptance of measures in the general population is faltering.

Our results also confirm the importance of household clustering while no indications for neighbourhood clustering were seen. The former finding is also supported by the observation that participants from higher income households were at non-significantly higher odds of SARS-CoV-2 sero-positivity. As we took into account total household income (not adjusted for number of persons in the household), single households were more likely to be in the lower income category and thus, at lower likelihood of household transmission.

We also saw a non-significant trend for lower odds of SARS-CoV-2 antibodies in smokers compared to nonsmokers (OR 0.2; 95% CI 0.02-1.1), confirming results of a meta-analysis [28] . Here, differences were mainly explained by differences in testing behaviour between smokers and non-smokers, which can be excluded in our study. One of the population-based studies published so far also indicated a lower SARS-CoV-2 sero-prevalence in smokers compared to non-smokers [29] . Whether this is a true effect of, e.g., nicotine [30] or vitamin D [31] or result of some form of bias needs to be evaluated in future studies. Of note is also the tendency for higher odds of SARS-CoV-2 in participants pursuing school, however, the wide confidence interval does not permit strong conclusions.

Among the strengths of our study are its populationbased, prospective nature in a large number of participants. Such population-based studies help authorities to plan public health measures based on the prevalence of exposure in the population, its spatial distribution and to further identify risk groups [32] . With increasing availability of vaccines, this study design with further followups will help public health authorities to understand the extent and duration of vaccine-induced immunity [33] . We previously showed a high sensitivity and specificity of the Elecsys ® Anti-SARS-CoV-2 assay (Roche) used in this sero-study [5] . For the follow-up, we developed and carefully validated a semi-automated protocol using self-sampled DBS for SARS-CoV-2 serology [7] . This approach facilitates field work to a very considerable extent and thus, makes studies with a higher frequency of follow-ups more feasible. Acceptance was high in our study population, and the percentage of participants lost to follow-up comparably low.

However, in the analyses we had to take selective participation into account by modelling the underlying non-response mechanism and calibrating the weights. This way, we could reduce attrition bias in our prevalence and incidence estimates. It is common in prospective cohort studies that baseline participants in younger age groups, with migration background and with lower socioeconomic status are less likely to participate at follow-up Fig. 7 Sero-incidence of SARS-CoV-2 between baseline and follow-up by A self-estimated health-related risk-taking behaviour, B sum of contacts and C leisure time activities in summer 2020 stratified for sex and age group. Results suggest some role of behavioural factors in the risk of SARS-CoV-2 sero-positivity but differences are not statistically significant [4] . While typically participants with positive outcome are also more likely to participate in follow-up studies (34) , our baseline participants who were SARS-CoV-2 antibody positive were less likely to take part at followup. This gives some indication that unknown sero-status motivated at least part of our baseline participants to take part in the study. Once positive sero-status was known to them, they might have lost interest. Further supported is this hypothesis by the fact that less participants were willing to complete the follow-up questionnaires than to take part in the SARS-CoV-2 antibody follow-up. As a consequence, statistical power to analyse the association between behavioural factor and SARS-CoV-2 sero-positivity was limited. Finally, the age group 14 to 19 years is of specific interest given that data on sero-prevalence in this age group is still limited. However, the number of participants in this age group was low (n = 212) and given the still low sero-prevalence at the time of the study, only 9 of them turned positive limiting the power of further analyses in this age group.

In conclusion, SARS-CoV-2 sero-prevalence in the Munich general population was still low by the end of 2020. Men and younger parts of the population were more likely to be affected. Risk-taking behaviour might be one reason for these differences. Therefore, non-pharmaceutical public health measures are still important.

95% CI: 95% Confidence interval; KoCo19: Munich COVID-19 cohort; LMU: Ludwig-Maximilians-University; OR: Odds ratio; ORCHESTRA : Connecting European Cohorts to Increase Common and Effective Response to SARS-CoV-2 Pandemic; SARS-CoV-2: Severe acute respiratory syndrome coronavirus type 2; Yrs: Years.

The online version contains supplementary material available at https:// doi. org/ 10. 1186/ s12879-021-06589-4.

Additional file 1: Figure S1 . For sensitivity and specificity adjusted (left) and unadjusted (right) SARS-CoV-2 sero-prevalence over the follow-up period excluding DBS intermediates. The 95% confidence intervals for the weekly sero-prevalence are based on the 2.5 and 97.5 percentiles from 5,000 repetitions of a cluster bootstrap that accounts for within household clustering. The estimates do not account for sample weights. Figure S2 . Comparison of the sero-positivity over the follow-up period accounting for sampling weights. Figure S3 . Sero-prevalence estimates across the different boroughs of Munich using CAR model priors with a single level of spatial autocorrelation as random effects. Figure S4 . Sero-prevalence estimates across the different boroughs of Munich with CAR model priors through the dissimilarity of the population density in neighbouring boroughs for spatial autocorrelation. Figure S5 . Distribution of mean within-cluster variance of test results under 10,000 random permutation of cluster assignments, with clusters being households, buildings, and geospatial clusters of different sizes. Household membership left invariant for building and geospatial clusters. Left: Value distribution. Right: 50%, 95%, 99% CIs. Black lines (left) and dots (right) indicate the observed values. Table S1 . Comparison of the follow-up plasma results, using DBS or Venous Blood Samples. Table S2 . Comparison of the estimates for the covariate of population density at the constituency level using the Poisson and the negative binomial model. Table S3 . Course of SARS-CoV-2 antibody status within the KoCo19 follow-up participants. Table S4 . Summary across age group and sex for the response to the behaviour questionnaires stratified by sex. p-values are from Pearson's Chi-squared test with simulated p-value (based on 10000 replicates). Table S5 . Summary across age group and sex for the response to the behaviour questionnaires stratified by age group. p-values are from Pearson's Chi-squared test with simulated p-value (based on 10000 replicates).

SARS-CoV-2 seroprevalence worldwide: a systematic review and meta-analysis

Population-based seroprevalence surveys of anti-SARS-CoV-2 antibody: an up-to-date review

Prevalence and risk factors of infection in the representative COVID-19 cohort Munich

Protocol of a population-based prospective COVID-19 cohort study Munich

A serology strategy for epidemiological studies based on the comparison of the performance of seven different test systems-the representative COVID-19 cohort Munich

Does length of questionnaire matter? A randomised trial of response rates to a mailed questionnaire

A dried blood spot protocol for high throughput analysis of SARS-CoV-2 Serology based on the Roche Elecsys anti-N assay

Adjusting coronavirus prevalence estimates for laboratory test kit error

Model assisted survey sampling

Calibration estimators in survey sampling

Variance estimation for complex statistics and estimators: linearization and residual techniques

Variance estimation under monotone nonresponse for a panel survey

Bootstrap-based improvements for inference with clustered errors

Nonparametric permutation tests for functional neuroimaging: a primer with examples

CARBayes: an R package for Bayesian spatial modeling with conditional autoregressive priors

Estimation of disease rates in small areas: a new mixed model for spatial dependence. In: Halloran ME, Berry D, editors. Models in epidemiology, the environment and clinical trials

Boundary detection in disease mapping studies

A comparison of conditional autoregressive models used in Bayesian disease mapping

Advances Bayesian multilevel modeling with the R package brms

convenient online submission • thorough peer review by experienced researchers in your field • rapid publication on acceptance • support for research data, including large and complex data types • gold Open Access which fosters wider collaboration and increased citations maximum visibility for your research: over 100M website views per year submit your research ?

brms: an R package for Bayesian multilevel models using stan

JointAI: joint analysis and imputation of incomplete data in R

Dealing with missing covariates in epidemiologic studies: a comparison between multiple imputation and a full Bayesian approach

Immunological memory to SARS-CoV-2 assessed for up to 8 months after infection. Science

Seroprevalence and humoral immune durability of anti-SARS-CoV-2 antibodies in Wuhan, China: a longitudinal, population-level, cross-sectional study

Evolution of antibodies against SARS-CoV-2 over seven months: experience of the Nationwide Seroprevalence ENE-COVID Study in Spain

Sex differences in immune responses that underlie COVID-19 disease outcomes

Addendum to S3-Guidelines Crohn's disease and ulcerative colitis: management of patients with inflammatory bowel disease in the COVID-19 pandemic-open questions and answers

The association of smoking status with SARS-CoV-2 infection, hospitalization and mortality from COVID-19: a living rapid evidence review with Bayesian meta-analyses (version 7)

Antibody status and cumulative incidence of SARS-CoV-2 infection among adults in three regions of France following the first lockdown and associated risk factors: a multicohort study

Editorial: Nicotine and SARS-CoV-2: COVID-19 may be a disease of the nicotinic cholinergic system

Vitamin D and COVID-19: is there a role?

The time to do serosurveys for COVID-19 is now

SARS-CoV-2 seroprevalence in COVID-19 hotspots

Third follow-up of the study on occupational allergy risks (SOLAR III) in Germany: design, methods, and initial data analysis

We gratefully thank all study participants for their trust, time, data, and specimens. This study would also not have been possible without the staff of the Division of Infectious Diseases and Tropical Medicine at the University Hospital of LMU Munich, Helmholtz Centre Munich, and Bundeswehr Institute of Microbiology, as well as all medical students involved. We thank Judith Eckstein for outstanding support regarding public relations. We thank the teams from the press offices of LMU, University Hospital of LMU Munich, and of Helmholtz Centre Munich. We thank the KoCo19 advisory board members 1 

Research (MoKoCo19, reference number 01KI20271). Euroimmun, Mikrogen, Roche, and Viramed provided kits and machines for analyses at discounted rates. The funders had no role in study design, data collection, data analyses, data interpretation, writing, or submission of this manuscript.

Our data are accessible to researchers upon reasonable request to the corresponding author taking data protection laws and privacy of study participants into account. To facilitate reproducibility and reuse, the analysis and figure generation code has been made available on GitHub (https:// github. com/ koco19/ epi2) and has been uploaded to ZENODO (https:// doi. org/ 10. 5281/ zenodo. 47070 37) for long-term storage.

The study was conducted in accordance with good clinical (GCP) and epidemiological practice (GEP) standards as well as the Declaration of Helsinki in its most recent form (as amended by the 64th WMA General Assembly, Fortaleza, Brazil, in October 2013). The study protocol was approved by the Institutional Review Board of the Medical Faculty at Ludwig Maximilian University Munich, Germany (opinion dated 31 March 2020; number 20-275; opinion date amendment: 10 October 2020), prior to study initiation. Written informed consent was obtained from all study participants prior to study inclusion.

Not applicable as no identifiable information was included in the manuscript.

In addition to the funding disclosed in the funding section, DM and VT are sub-investigators vaccines trials sponsored by Curevac AG. LO received non-financial support from Dr. Box Betrobox and grants from the Bavarian State Ministry of Science and the Arts during the conduct of the study. CR received royalities from Elsevier LTD for the textbook "Clinical Cases in Tropical Medicine". She also received honoraria for lectures done at non-private academic institutions. CR is also chair of the "Travel medicine" board of the of the German Society for the Tropical Medicine, Travel Medicine and global health (unpaid) and member of the permanent committee on vaccinations at the Robert-Koch-Institute work group on travel vaccinations (unpaid). AW is on Roche Diagnostics and Roche Pharma Advisory Boards related SARS-CoV-2. The Bundeswehr Institute of Microbiology (RW) received research funding from the Medical Biodefence program of the Bundeswehr. The funders had no role in study design, data collection, data analyses, data interpretation, writing, or submission of this manuscript.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.