key: cord-0816585-ioakytdl authors: Richter, E.; Liebl, D.; Schulte, B.; Lehmann, N.; Fuhrmann, C.; Joeckel, K.-H.; Ioannidis, J.; Streeck, H. title: Analysis of Fatality Impact and Seroprevalence Surveys in a Community Sustaining a SARS-CoV-2 Superspreading Event date: 2022-01-28 journal: nan DOI: 10.1101/2022.01.26.22269805 sha: f57e65e0d0720e0141889491672193c18d337562 doc_id: 816585 cord_uid: ioakytdl There is ongoing debate on the COVID-19 infection fatality rate (IFR) and the impact of COVID-19 on overall population mortality. Here, we addressed these issues in a community in Germany with a major superspreader event analyzing deaths over time as well as auditing death certificates in the community.18 deaths that occurred within the first 6 months of the pandemic in the community had a positive test for SARS-CoV-2. Six out of 18 SARS-CoV-2+ deaths had non-COVID-19 related causes of death (COD). Individuals with confirmed infection and COVID-19 COD typically died of respiratory failure (75%) and tended to have fewer reported comorbidities (p=0.029). Duration between first confirmed infection and death was negatively associated to COVID-19 being COD (p=0.04). Repeated seroprevalence essays on an original sample of 587 individuals in three visits showed modest increases in seroprevalence over time, and substantial seroreversion (30% [27/90] (95% CI: [20.5%; 39.5%])). IFR estimates accordingly varied depending on COVID-19 death attribution and seroprevalence caveats. Careful ascertainment and audit of COVID-19 deaths are important in understanding the impact of the pandemic. The COVID-19 pandemic has been associated with considerable mortality at a global level, but estimates about the total number of deaths have varied widely across populations, demographics and countries [1] [2] [3] [4] [5] . There is ongoing debate on the population specific COVID-19 fatality rates (IFRs) and on the impact of COVID-19 on overall population mortality in various settings 6 . The IFR would be a robust marker for the severity of a virus-induced disease if the number of infected individuals (IFR denominator) as well as the number of deaths (IFR numerator) could be correctly measured. The number of infected individuals can be inferred by seroprevalence studies, and many such studies have been done 7, 8 . However, biases often exist in such studies. Two main concerns are whether the surveyed sample is representative of the general population and whether people who get infected do not develop or lose their detectable antibodies over time. The validation of the number of COVID-19-related deaths has received less attention, but both over-and under-counting may occur 9 . COVID-19 deaths may have been underestimated in particular in some locations at the beginning of the pandemic when testing was limited, or because determining the cause of death (COD) in the elderly or those with multiple comorbidities was difficult. Conversely, COVID-19 deaths may have been over-attributed, e.g. if based on a positive SARS-CoV-2 test without a relevant clinical picture. Auditing of death certificates and/or medical records may help re-assess COVID-19 death attribution. For example, auditing efforts decreased the number of COVID-19 deaths in Alameda County and in Santa Clara county by about a quarter compared with those that had been originally reported 10, 11 . Here, we address these issues in a community that witnessed the first major SARS-CoV-2 outbreak in Germany, as carnival festivities around February 15, 2020 were followed by a massive outbreak. To investigate this outbreak, we previously reported a seroepidemiological observational study that found a modest IFR 12 . However, questions emerged whether early during the pandemic COVID-19-related deaths may have been missed, mislabeled or underestimated. Therefore, in the currently presented work we analyzed deaths in the studied community over time and audited death certificates during the early phase of the pandemic. Specifically, we requested official death certificates of individuals who died in the community between March and October of 2020 and analyzed in detail the underlying causes for their deaths. Furthermore, we also assessed seroprevalence with repeated surveys over a year and specifically assessed seroreversion. These analyses allow obtaining better insights on IFR estimates and on the impact of the pandemic in this intensively-studied community. Between March and October 2020 104 inhabitants of the studied community (n=12,597) deceased according to information we obtained from the state office for information and technology in North Rhine-Westphalia (Germany). With permission by the government of North Rhine-Westphalia (Germany), the state attorneys provided us with the death certificates of 77 subjects. Out of the remaining 27 deaths none was attributed to COVID-19. All death certificates were assessed by three of the authors (E.R., D.L., and H.S.) and consensus was achieved by applying the following procedures: We extracted information about gender, date of birth, date of death, COD, comorbidities and any SARS-CoV-2 infection on record with the corresponding date of the first positive test. We classified causes of deaths in the following categories: respiratory failure; sepsis; cardiovascular disease (CVD); liver failure; kidney failure; cancer; other as well as comorbidities in the categories: respiratory; sepsis; CVD; liver malfunction; kidney malfunction; cancer; neurological; diabetes; unknown; other. We also compared the number of deaths attributed to COVID-19 in the data submitted by the authorities to the national death counts versus the number of deaths attributed to COVID-19 after our audit. In attributing deaths to COVID-19 in our audit of death certificates, we followed the guidance of the WHO guidelines. Accordingly, a COVID-19-related death is defined as resulting from a clinically compatible illness in a probable or confirmed COVID-19 case, unless there is a clear alternative COD that cannot be related to COVID disease (e.g. trauma) and a death due to COVID-19 may not be attributed to another disease (e.g. cancer). In addition, there should be no period of complete recovery from COVID-19 between illness and death 13 . We could not obtain medical records for the deceased for further in-depth audit of these records. We also retrieved information on the number of individuals in the community who live in institutionalized settings (nursing homes). It was not possible to retrieve separate counts for deaths occurring among nursing home residents either for all deaths or for COVID-19 deaths specifically. 919 study participants, which we already enrolled in our previous study in April 2020 12 were contacted in October 2020 (Visit 1), January 2021 (Visit 2) and April 2021 (Visit 3) by letter and were invited to an study acquisition center. After having provided written and informed consent, study participants completed a questionnaire querying demographics, symptoms, underlying diseases and medication, as well as their SARS-CoV-2 vaccination status. For children under 18 years, written and informed consent was provided by the persons with care and custody of the children following aged-adapted participant information. The 411 people (n=297 of them from this specific community) who had participated in the main carnival session "Kappensitzung" 14 were enrolled in a separate study. Furthermore, study participants were asked to provide blood specimens and pharyngeal swabs. Over the course of all three follow-up visits, in total 587 participants from the originally study participated and provided blood specimens as well as pharyngeal swabs. Total participants for each visit are: Visit (1) 587; Visit (2) 488 and Visit (3) was assessed in each sample by real time RT-PCR (SuperScript™III One-Step RT-PCR System with Platinum™ TaqDNA Polymerase, Thermo Fisher). The following primers were used, for E gene: E_Sarbeco_F1 and R, and probe E_Sarbeco_P1, for RdRP gene: RdRP_SARSr_F, and R, and probe RdRP_SARSr-P2. In addition, an internal control for RNA extraction, reverse transcription, and amplification was applied to each sample (innuDETECT Internal Control RNA Assay, Analytik Jena #845-ID-0007100). If amplification occurred in both virus-specific reactions samples were considered positive. Plasma was used to determine Anti-SARS-CoV-2 IgG (against S1 domain) using enzyme-linked immunosorbent assays (ELISA) on the EUROIMMUN Analyzer I platform. According to the manufacturer's instructions a result was considered positive when a ratio (extinction of sample/extinction of calibrator) of 0.8 or higher was reached. In addition, the Roche Elecsys® N ELISA was used to determine antibody response against nucleocapsid (N) of SARS-CoV-2. The result is given in COI (cutoff index), which is positive between 1.0 and approximately 250. Further we estimate for each visit after the original one how many participants had undergone seroreversion, defined as ratio of IgG below 0.8. Finally, for each visit we report on the number and proportion of participants who had tested positive by RT-PCR and we compare this against the proportion of people in the overall community who had tested positive by RT-PCR and were thus considered to be documented cases. Analysis of Death certificate data: Parametric normality assumptions were not applied to avoid distributional misspecifications in moderate to small samples where a central limit theorem may not be applicable. Fisher's exact two-sided test was used to test for independence between two factor variables with two levels (2x2 designs). Pearson's Chisquared test was used to test for independence between factor variables and Pearson's' Chisquared goodness-of-fit test was used to test hypotheses about categorical distributions. For both chi-square tests, the p-values were simulated using Monte-Carlo simulations with 50,000 replications to consistently approximate the small sample distributions under the nullhypotheses 15 . The one-sided two-sample Kolmogorov-Smirnov test was applied to test for a location shift in the number of days from a positive COVID-19 test and death between All rights reserved. No reuse allowed without permission. perpetuity. preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in The copyright holder for this this version posted January 28, 2022. were adjusted for multiple testing following the Bonferroni-Holm procedure. Confidence Intervals (CIs) for the IFR were computed by dividing the number of deaths by the CI limits of the estimated number of infected. We estimated IFR using different acquisition periods (7, 20, 35 , and 60 days after the initial seroprevalence survey 12 ) as this has been a debated issue in All rights reserved. No reuse allowed without permission. perpetuity. preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in The copyright holder for this this version posted January 28, 2022. ; https://doi.org/10.1101/2022.01.26.22269805 doi: medRxiv preprint previous meta-analyses of IFR 18, 19 . The numbers of SARS-CoV-2-associated deaths in the studied community for given lengths of acquisition periods can be considered fixed as the data acquisition corresponded to a complete survey of all recorded SARS-CoV-2-associated deaths in the community during the study period. In addition, we generated also separate IFR estimates for the ages groups 0-54, 55-74 and > 75 years. To additionally account for possible uncertainty in the number of SARS-CoV-2-associated deaths we use a Bayesian credibility interval for the IFR which, however, can only be interpreted for populations for which the full survey in the community means a representative sample. This CI was computed as the perpetuity. preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in The copyright holder for this this version posted January 28, 2022. ; In total we analyzed the death certificates of 77 individuals who died between March 7 th and October 9 th 2020 in the studied community. The median age was 82 years [range 28 -98 ] and 55.8% were male. Among the 77 deaths, overall 18 individuals had been tested SARS-CoV-2 positive, of whom 61.1% were male. Among the remaining 59 individuals we screened for potential signs of COVID-19 infection but no case of acute respiratory disease (ARDS) or similar cause of deaths that may be related to COVID-19 was found, thus we excluded the possibility of an undiagnosed SARS-CoV-2 infection. For the 18 individuals who had tested positive for SARS-CoV-2, we analyzed the date of death as well as the main COD ( Table 1) . August 2020, and one individual in October 2020. Deaths per month and age stratum from the 77 deaths we were able to examine are shown in Supplemental Table 1 test. Overall, the likelihood to die of COVID-19 decreased markedly 100 days or more after confirmed SARS-CoV-2 infection (Supplemental Figure 1 ). Infection fatality rate Figure 3 shows the estimated IFR with different acquisition periods ranging from 7 days (7 deaths), 20 days (8 deaths), 35 days (9 deaths), and 60 days (12 deaths). Extending the observation period further did not lead to a change in the IFR depending on the acquisition period. The IFR estimates ranged accordingly from 0.36% to 0.62%, the 95% confidence intervals ranged from 0.28% to 0.84% and this confidence range would become 0.18% to 1.15% when introducing also uncertainty in the number of deaths. Table 2 shows the IFR estimates for age strata 0-50, 55-74, and >75 years. As no deaths occurred in the age stratum 0-50, the IFR estimates are all 0.0% For the 55-74 years the IFR estimates range from 0.18 to 0.74 and the 95% confidence intervals extend from 0.13 to 1.19 For people >75 years old, the IFR ranges from 3.21 to 4.28 and the 95% confidence intervals extend from 1.87 to 9.21. for sensitivity and specificity of IgG (sensitivity 90.9%; specificity 99.1%). Next, we wanted to understand how many participants had undergone seroreversion, defined as IgG values below the ratio of 0.8. Therefore, we analyzed participants who were IgG(+) at baseline and participated at least in one of the three follow-up visits. We could see differences in the IgG levels of the participants already 6 months after the initial visit. IgG was significantly decreased All rights reserved. No reuse allowed without permission. perpetuity. preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in The copyright holder for this this version posted January 28, 2022. ; https://doi.org/10.1101/2022.01.26.22269805 doi: medRxiv preprint over time and showed reduced antibody titers at 6 and 9 months (p<0.0001) (Figure 4) . In total 27 (30%) out of 90 participants (95% CI: [20.5%; 39.5%]) who were IgG(+) at time point 0 seroreversed and showed no detectable IgG levels after 12 months. Finally, we analyzed the pharyngeal swabs via RT-PCR to identify active SARS-CoV-2 infection. Overall, only 2 participants were tested positive in our cohort over the 12 months of follow-up. One participant at 6 months, who was IgG(+) at baseline but had no detectable IgG antibodies 6 months later. The second participant was RT-PCR(+) at 9 months and was IgG negative at 0 and 6 months. More importantly, the proportion with positive PCR in our seroprevalence-tested cohort was far smaller (2 out of 587 tested participants, i.e. 0.3%) than the respective proportion in the overall population of the community (n=374 PCR confirmed cases, i.e. 3.0%) in the same period of time, between October 2020 and April 2021. In this study, we audited and analyzed the deaths occurring in a German community between March and October of 2020 at the beginning of the SARS-CoV-2 pandemic. In total 18 dying individuals were reported positive for SARS-CoV-2, but COVID-19 was the underlying COD in only 12 out of the 18 deaths with confirmed SARS-CoV-2 infection. Therefore, one third of individuals with confirmed SARS-CoV-2 infection died with non-COVID-19 related COD. The duration between confirmed SARS-CoV-2 infection and death was significantly higher in individuals without COVID-19 related COD compared to those with COVID-19 related COD. Moreover, the number of deaths reported to national authorities as COVID-19 deaths was 25% higher than the truly "causal" COVID-19 deaths. This over-estimation of COVID-19 deaths is similar to the rates seen in audits in US counties 10, 11 . Over-estimation of COVID-19 deaths may apply also to other settings, especially in high-income countries where intensive testing is performed 9 . We also observed that more SARS-CoV-2 positive individuals with non-COVID-19 underlying COD had listed comorbidities on the death certificates compared to individuals with COVID-19 underlying COD. This finding needs to be interpreted with caution given the small All rights reserved. No reuse allowed without permission. perpetuity. preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in The copyright holder for this this version posted January 28, 2022. 19, 21 . Reported studies of COVID-19-associated deaths in Germany and the UK determined septic shock, multi organ failure and respiratory failure as the most common immediate COVID-19 underlying COD, often due to suppurative pulmonary infection and diffuse alveolar damage [22] [23] [24] . Our findings for our studied community are consistent with this, as 75% of the COVID-19-associated deaths died from respiratory failure in combination with septic shock and multi organ failure. Furthermore, it was described that smoking or some comorbidities like cancer and chronic liver disease had stronger associations with non-COVID than COVID-19 deaths in association with age 22 . Although we do not see significant associations with age or specific comorbidities (probably due to the low number of cases), we determined kidney failure, lactic acidosis or cancer as the main causes in non-COVID-19 associated deaths. These findings support the importance of looking more carefully into the deaths of individuals with a confirmed SARS-CoV-2 infection and taking into account the individuals' medical history along with their most recent medical data and symptoms. All rights reserved. No reuse allowed without permission. perpetuity. preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in The copyright holder for this this version posted January 28, 2022. ; https://doi.org/10.1101/2022.01.26.22269805 doi: medRxiv preprint Recent publications on COVID-19 disease progression analyzed the length of stay in the hospital and survival time [25] [26] [27] . It was reported that the median length between symptom onset and hospitalization ranged between 3 and 10.4 days, depending on the age of the patient 26 . However, length of stay in hospital for patients who died eventually were an additional six to seven days as reported in a Belgian study 26 , and for China 27 Evaluation of the deaths per month and age stratum showed a peak of deaths in March 2020 and it was concentrated in people over 70, consistent with the age of the early COVID-19 fatalities. We could not retrieve information on how many of these individuals might have been nursing home residents. However, this specific community has a high proportion of its population residing in nursing homes (about 2 times larger than the average for Germany [731,000 people, 0.9% of the German population]) 28 . It is possible that many of the COVID-19 fatalities in the community might have been in people with limited life expectancy, regardless of whether they were institutionalized or not. Obviously, one has to be cautious with these inferences, especially given the relatively small number of deaths. Of note, eventually the overall number of deaths in the community for calendar year 2020 (n=157) was similar to calendar years 2018 (n=156) and 2015 (n=158) but higher than other previous recent years. However, it is very likely that this was influenced by the hard lockdown in the studied community in March and April 2020. Furthermore, repeated seroprevalence surveys showed that seroprevalence changed only modestly in the fall of 2020 and beyond as compared with our previously published results on seroprevalence in April 2020 12 . However, we documented substantial rates of All rights reserved. No reuse allowed without permission. perpetuity. preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in The copyright holder for this this version posted January 28, 2022. seroreversion, which is in agreement with some other studies [29] [30] [31] [32] [33] [34] [35] [36] . Some studies have found even higher rates of seroreversion, e.g. one investigation found a median time to seroreversion for IgG being only 55 days 30 , but this may include false-positives at initial screening. Nevertheless, it is very likely that seroreversion causes the number of infected individuals to be substantially under-estimated and thus the IFR to be over-estimated 29, 30 . This may apply even to our early April 2020 survey that happened within two months of the superspreader event. Estimates of number of people infected in the population may also be affected by the representativeness of the surveyed sample. In our original seroprevalence survey in April 2020 12 , we had noticed that the tested sample had a lower proportion of documented RT-PCR infections than the overall community population (2.39% versus 3.08%). Correction for this factor would decrease IFR estimates by 29% (e.g. from 0.36-0.62% to 0.28-0.48%). The sampling deviation was seen also, even more prominently in the follow-up surveys, where only 2 participants tested positive by RT-PCR (10-fold less than in the general population). It may reflect the fact that participants (especially those who return also for follow-up visits) may be more health conscious; or alternatively, less likely to perform RT-PCR testing since they are tested for antibodies. We also observed that our original seroprevalence study sample was under-representative of individuals infected during the main carnival event, the Kappensitzung, where participants had very high infection rates 14 , three times higher than the general population of the community by April 2020. Due to the small number of cases, we were not able to exclude that this may have been entirely due to chance but the pattern is consistent with the possibility that seroprevalence is under-estimated (and, correspondingly, IFR overestimated). These observations highlight the difficulty of using seroprevalence samples and repeated surveys to assess the number of infected individuals in a population. Our observations show that the calculation of IFRs needs to be done very carefully and that their interpretations need to take into account all possible influencing factors. In our study, over-estimation of COVID-19 deaths and under-estimation of seroprevalence may have inflated the IFR estimate. Additional variability in the estimates may be introduced by the time All rights reserved. No reuse allowed without permission. perpetuity. preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in The copyright holder for this this version posted January 28, 2022. ; https://doi.org/10.1101/2022.01.26.22269805 doi: medRxiv preprint window used for capturing deaths, its relationship to the time window of the seroprevalence, and the assumptions made about the delay in developing antibodies and in dying after infection. In the case of this specific community, these uncertainties would probably still be captured by the original confidence interval that we reported for IFR in our original publication and that allowed for uncertainty in the number of deaths (95% CI, 0.17% to 0.77%). Moreover, the exact case mix of infected individuals can have a major impact on the IFR, given the extremely steep age gradient that has been documented before and which we also saw prominently in the community data 5, 37 . It should be noted that due to the study-design, the results of this study can only be representative for this specific community. Overall, in seroprevalence studies where only a tiny portion of the population is selected for serological testing (e.g. typically 0.01-1% in nation-wide surveys 38 ) one should remain cautious, as errors in counting deaths, seroprevalence, and other sources of uncertainty may have an larger impact. In conclusion, our in-depth assessment of the fatality impact of COVID-19 in the community shows a relatively low fatality rate in this community with deaths concentrated entirely in the elderly, but seroprevalence estimates that go into IFR calculations need careful considerations of their time window, potential seroreversion, and representativeness. Importantly, COVID-19 deaths may have been overcounted and the relative contribution of over-and under-counting of COVID-19 deaths needs careful auditing across multiple other locations. All rights reserved. No reuse allowed without permission. perpetuity. preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in The copyright holder for this this version posted January 28, 2022. ; Figure 1 : Age differences. All 77 deaths were differentiated in three groups: SARS-CoV-2 negative, SARS-CoV-2 positive and COVID-19 underlying COD as well as SARS-CoV-2 positive but no COVID-19 underlying COD. There were no significant differences in the age distribution among COVID-19 infected or non-infected individuals. All rights reserved. No reuse allowed without permission. perpetuity. preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in The copyright holder for this this version posted January 28, 2022. ; https://doi.org/10.1101/2022.01.26.22269805 doi: medRxiv preprint All rights reserved. No reuse allowed without permission. perpetuity. preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in The copyright holder for this this version posted January 28, 2022. ; https://doi.org/10.1101/2022.01.26.22269805 doi: medRxiv preprint Figure 3 : Newly calculated infection fatality rate (IFR) with 95% confidence intervals (blue) and 95% credibility intervals (gray) additionally accounting for uncertainty in the number of deaths. All rights reserved. No reuse allowed without permission. perpetuity. preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in The copyright holder for this this version posted January 28, 2022. ; Figure 4 : IgG levels over time of participants who were IgG(+) at baseline and participated at least in one of the three follow-up visits were analyzed. 27 participants who developed IgG values below the ratio of 0.8 (Limit of Quantification) were considered as seroreversed. Statistical significance was assessed using Wilcoxon sign-rank test and the resulting six pvalues were adjusted for multiple testing following the Bonferroni-Holm procedure. Visit 1 vs. 0: p<0.0001; Visit 2 vs. 0: p<0.0001; Visit 3 vs. 0: p=0.0016; Visit 1 vs. 2: p<0.0001; Visit 1 vs. 3: p=0.0053; Visit 2 vs. 3: p<0.0001. All comparisons show a decreasing trend with time. It is important to note, that although the graph only shows those with IgG ratio ≥ 0.8 at visit 0, all subjects who contribute data for visit 0 and at least one of the following visits are included in the test (without restriction at visit 0). All rights reserved. No reuse allowed without permission. perpetuity. preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in The copyright holder for this this version posted January 28, 2022. ; https://doi.org/10.1101/2022.01.26.22269805 doi: medRxiv preprint Excess mortality due to COVID-19 in Germany Estimation of Excess Deaths Associated With the COVID-19 Pandemic in the United States Excess all-cause mortality during the COVID-19 pandemic in Europe -preliminary pooled estimates from the EuroMOMO network Mortality impacts of the coronavirus disease (COVID-19) outbreak by sex and age: rapid mortality surveillance system Factors associated with COVID-19-related death using OpenSAFELY Reconciling estimates of global spread and infection fatality rates of COVID-19: An overview of systematic evaluations Global seroprevalence of SARS-CoV-2 antibodies: A systematic review and meta-analysis SeroTracker: a global SARS-CoV-2 seroprevalence dashboard Over-and under-estimation of COVID-19 deaths Alameda County revises COVID-19 death count by 25% County's COVID-19 Death Toll Drops by 505 After Data Change Infection fatality rate of SARS-CoV2 in a super-spreading event in Germany International guidelines for certifcation and classification (Coding) of COVID-19 as cause of death Analysis of the Dynamics, Outcome, and Prerequisites of the first German SARS-CoV-2 Superspreading Event. medRxiv Simplified Monte Carlo Significance Test Procedure Estimating prevalence from the results of a screening test Bootstrap Methods and their Application First provisional results of the 2020 causes of death statistics including data on COVID-19 and suicides Different approaches to the analysis of causes of death during the COVID-19 epidemic Coronavirus im Kreis Heinsberg COVID-19 as the underlying cause of death: disentangling facts and values Factors associated with deaths due to COVID-19 versus other causes: population-based cohort analysis of UK primary care data and linked national death registrations within the OpenSAFELY platform Causes of death and comorbidities in hospitalized patients with COVID-19 Effects of underlying morbidities on the occurrence of deaths in COVID-19 patients: A systematic review and meta-analysis Time-to-Death approach in revealing Chronicity and Severity of COVID-19 across the World Time between Symptom Onset, Hospitalisation and Recovery or Death: Statistical Analysis of Belgian COVID-19 Patients Estimates of the severity of coronavirus disease 2019: a model-based analysis Mehr als 730.000 Menschen leben in Deutschland in Pflegeeinrichtungen Levels of SARS-CoV-2 population exposure are considerably higher than suggested by seroprevalence surveys Duration of SARS-CoV-2 sero-positivity in a large longitudinal sero-surveillance cohort: the COVID-19 Community Research Partnership The Fall in Antibody Response to SARS-CoV-2: a Longitudinal Study of Asymptomatic to Critically Ill Patients Up to 10 Months after Recovery SARS-CoV-2 antibody magnitude and detectability are driven by disease severity, timing, and assay Estimating the Cumulative Incidence of SARS-CoV-2 Infection and the Infection Fatality Ratio in Light of Waning Antibodies Serological surveillance of SARS-CoV-2: Six-month trends and antibody response in a cohort of public health workers Time series analysis and mechanistic modelling of heterogeneity and sero-reversion in antibody responses to mild SARS-CoV-2 infection Decline in SARS-CoV-2 Antibodies After Mild Infection Among Frontline Health Care Personnel in a Multistate Hospital Network -12 States Population-level COVID-19 mortality risk for non-elderly individuals overall and for non-elderly individuals without underlying diseases in pandemic epicenters Germany's low SARS-CoV-2 seroprevalence confirms effective containment in 2020: Results of the nationwide RKI-SOEP study. medRxiv We would like to thank all the subjects who participated in our repeated seroprevalence surveys. Furthermore, we would like to thank the following people who helped with the study: The authors have declared no competing interest. However, the government of North Rhine-Westphalia (Germany) supported part of the study. The idea, the plan, the concept, protocol, the conduct, the data analysis, and the writing of the manuscript of this study were independent of any third parties, including the government of North Rhine-Westphalia, Germany.