Association between the COVID-19 pandemic and pertussis in France using multiple nationwide data sources

Background: Interventions to mitigate coronavirus disease 19 (COVID-19) pandemic may impact other respiratory diseases such as pertussis. We aimed to study the course of pertussis in France over an 8-year period and its association with COVID-19 mitigation strategies, using multiple nationwide data sources. Methods: We analyzed the number of French pertussis cases between 2013 and 2020, using the PCR test results from nationwide outpatient laboratories (Source 1) and the pediatric network of 41 hospitals (Source 2), and using the reports of an office-based pediatric national network (Source 3). We conducted a quasi-experimental interrupted time-series analysis, relying on negative binomial regression models. The models accounted for seasonality, longterm cycles, and secular trend, and included a binary variable for the first national lockdown (ordered on March 16, 2021). Results: We identified 19,039 cases of pertussis from the three data sources during the study period. There was a significant decrease of pertussis cases following the implementation of mitigation measures, with adjusted incidence rate ratios of 0.102 (95% CI 0.040-0.256) and 0.216 (95% CI 0.071-0.656) for Source 1 and Source 2, respectively. The association was confirmed in Source 3 (median of 1 [IQR 0-2] vs. 0 [IQR 0-0] pertussis cases per month before and after lockdown, respectively, p=0.0048). Conclusion: The drastic reduction of outpatient and hospitalized cases of pertussis strongly suggests an impact of COVID-19 mitigation measures and their consequences on pertussis epidemiology. Pertussis vaccination recommendations should be carefully followed, and disease monitoring should be continued to detect any resurgence after relaxation of mitigation measures.


Methods:
We analyzed the number of French pertussis cases between 2013 and 2020, using the PCR test results from nationwide outpatient laboratories (Source 1) and the pediatric network of 41 hospitals (Source 2), and using the reports of an office-based pediatric national network (Source 3). We conducted a quasi-experimental interrupted time-series analysis, relying on negative binomial regression models. The models accounted for seasonality, longterm cycles, and secular trend, and included a binary variable for the first national lockdown (ordered on March 16, 2021).

Results:
We identified 19,039 cases of pertussis from the three data sources during the study period. There was a significant decrease of pertussis cases following the implementation of mitigation measures, with adjusted incidence rate ratios of 0.102 (95% CI 0.040-0.256) and 0.216 (95% CI 0.071-0.656) for Source 1 and Source 2, respectively. The association was confirmed in Source 3 (median of 1 [IQR 0-2] vs. 0 [IQR 0-0] pertussis cases per month before and after lockdown, respectively, p=0.0048).

Conclusion:
The drastic reduction of outpatient and hospitalized cases of pertussis strongly suggests an impact of COVID-19 mitigation measures and their consequences on pertussis epidemiology. Pertussis vaccination recommendations should be carefully followed, and disease monitoring should be continued to detect any resurgence after relaxation of mitigation measures.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Funding
There was no specific funding for the study.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

INTRODUCTION
Whooping cough or pertussis is a highly contagious respiratory disease transmitted from human to human via aerosolized respiratory droplets, and is mainly caused by Bordetella pertussis. A resurgence of pertussis has been observed in many countries during the last decade despite widespread vaccine implementation. 1 Most pertussis morbidity and mortality are due to severe clinical forms in young infants that usually require admission to intensive care units. In 2014, the World Health Organization (WHO) estimated pertussis as a cause of 160,700 deaths in children less than five years-old. 2 Since December 2019, the world has been facing another respiratory infectious disease, the Coronavirus disease 2019 (COVID-19) pandemic due to SARS-CoV-2. 3 In France, the lockdown was ordered at the beginning of the first wave of the COVID-19 pandemic on March 16, 2020. Measures included school closures, shortly followed by face covering in public places on May 11, 2020. 4 During the following months, the French government successively introduced several bundles of COVID-19 mitigation measures based upon the dynamics of SARS-CoV-2 infection. The French population initially faced difficulties of access to clinical and biological diagnosis, and delay of childhood vaccinations despite the mandatory infant vaccinations since 2018 and booster vaccination guidelines. 5 We hypothesized that such measures and their consequences might have impacted pertussis epidemiology, as suggested for other transmissible airborne infectious diseases. 6 Thanks to a well-established surveillance system of pertussis cases through networks of outpatient laboratories and voluntary hospitals, and the French ambulatory surveillance for outpatient pediatric cases, we assessed the course of pertussis epidemiology in France between January 1, 2013, and December 31, 2020, and its potential association with COVID-19 mitigation strategies.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted July 20, 2021.

Data sources and case definitions
We performed a retrospective analysis using three nationwide data sources from pertussis surveillance systems in France. Those 9 We also included data from the French ambulatory surveillance of whooping cough for outpatient pediatric cases. Implemented in 2002 by the Association Clinique et Thérapeutique Infantile du Val de Marne (ACTIV) and the French NRC for Bordetella infections, this surveillance (office-based pediatrician network pertussis surveillance, Source 3) does not aim to estimate the exhaustive number of pertussis cases but to estimate the duration of immunity conferred by pertussis vaccines in children in France. 10 From 2013 to 2020, 76 pediatricians from the ACTIV network throughout France reported patients with suspected pertussis (cough illness). All suspected patients were invited to perform biological confirmation of pertussis. Cases were defined as either 'confirmed' (i.e., positive on culture and/or PCR, or 'epidemiological' (i.e., acute cough lasting at least 14 days in a child or adult who was in contact with a confirmed case). 10 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted July 20, 2021. ; https://doi.org/10.1101/2021.07. 16.21260367 doi: medRxiv preprint Strengthening the Reporting of Observational studies in Epidemiology (STROBE) guidelines were followed to report the study (Supplementary Table1).

Statistical analysis
We conducted a quasi-experimental interrupted time-series analysis relying on negative binomial regression models, 11 with a time-unit of one month, to analyze the changes in the incidence of positive Bordetella PCRs over time. The models accounted for seasonality (using pairs of sine/cosine terms), secular linear trend, a binary variable to define periods before and after the 1 st lockdown (i.e., April 1 st , 2020, Supplementary Figure), and over-dispersion of data. The models also included a dummy variable (with 36-month periods) to adjust for longterm cycles commonly observed in pertussis epidemiology. 8 First, the analysis was carried out on our most comprehensive data source (Source 1), overall, and then stratified across three pre-specified age groups: pre-school (0-5 years), primary and secondary school children (6-17 years), and adults (≥18 years). We further performed the following sensitivity analyses to assess the robustness of our findings: (1) by fitting the data from Lab 1 and Lab 2 separately, (2) by including the negative samples when available (i.e., modeling the proportion of positive PCR results), (3) by using a different modeling strategy relying on segmented linear regression with autoregressive errors to account for autocorrelation in the data. 11 Second, the analysis was carried out on data from Source 2, relying on the same negative binomial regression model used for Source 1. Third, the data from Source 3 were not included in the time-series modeling because these data, collected to evaluate the duration of vaccine-induced immunity, did not fulfill model requirements; an analysis comparing the number of pertussis cases per month before and after the lockdown was carried out. Differences between groups were assessed by the Mann-. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 20, 2021.

Ethics approval
The data collections received approval by Institutional Review and Data protection Board

Outpatient laboratory pertussis surveillance (Source 1)
Over the study period, there were 17,912 pertussis cases reported from the two private laboratories (Source 1). The median age of patients with a positive PCR was 20.0 (interquartile range [IQR] 6.5-46.2) years, and distribution among age groups was 23% for the 0 -5 years, 24% for the 6 -17 years, and 53% for the is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 20, 2021. In sensitivity analysis, our findings did not differ when data from Lab 1 and Lab 2 were fitted separately ( Table 1) Table 2). Segmented linear regression models yielded results similar to negative binomial modeling, with a significant drop in the number of positive PCR results across all age groups ( Table 1).

Hospital-based pertussis surveillance (Source 2)
Over the study period, 992 positive PCR tests were collected through the Renacoq hospital network (Source 2). Distribution among age groups was 61% for the 0 -2 months, 25% for the 3-5 months, and 14% for the 6-11 months. We also observed a seasonal pattern ( Figure   1C). Before lockdown, the average number of pertussis cases per month was 11.3 (SD 7.3).
After the lockdown, a sharp decrease in the number of hospitalized pertussis cases was confirmed, with an average number of the case reaching 1.1 (SD 2.3), p<0.001. In time-series modeling, there was a significant association between COVID-19 mitigation measures and pertussis cases, with a reduction of 78.4% (95%CI 34.4; 92.9); aIRR 0.216 (95% CI 0.071; . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 20, 2021. ; https://doi.org/10.1101/2021.07.16.21260367 doi: medRxiv preprint 1 1 of pertussis cases close to null after the lockdown, a level never observed since the beginning of the French pertussis surveillance in 1996 7 .
All age groups were affected by the decrease of pertussis cases following the mitigation measures. However, a lower decrease was observed in the 0-5-year group tested by the two private laboratories. One explanation could be that the youngest children were more likely to visit their pediatrician despite lockdown measures. Another explanation could be that transmission might still have been active in this age group during the first months of 2020.
We did not find a significant association between the proportion of positive PCRs and the mitigation measures against SARS-CoV-2. Our data show that the overall number of PCR tests decreased right after the lockdown, while pertussis might still have circulated in households. However, despite a steady number of PCR tests performed during the whole postlockdown period, the number of positive cases reached almost zero within the last 4 months of 2020; of note access to clinical and biological diagnosis went back to normal during that period. All these findings suggest an important decrease of pertussis circulation in France.
Previous studies revealed a major impact of COVID-19 pandemic mitigation measures on transmissible diseases, such as, in the pediatric population, gastroenteritis, common cold, bronchiolitis, and acute otitis media. 6 However, a resurgence of some other airborne infectious diseases has been observed in the setting of relaxed physical distancing recommendations, such as the delayed inter-seasonal Respiratory Syncytial Virus (RSV) epidemic described by Foley et al. 13 For now, we did not observe any pertussis resurgence despite the end of the first lockdown and school fully reopening in France. Our data suggest that, unlike RSV, slight relaxation of public health measures (with maintained social distancing and mask-wearing) has no impact on pertussis dynamics, as for the influenza virus for which no outbreak has been identified worldwide this year. 14, 15 We identified five studies from Australia, South Korea, and Italy, describing the trends of several infectious diseases, . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 20, 2021. including pertussis, before and after implementing their respective COVID-19 mitigation measures. [16][17][18][19] Although the power of these studies was much lower compared to ours (lower duration of surveillance and description of a maximum of ~1500 cases, vs. 19,039 cases in our study), they all found a decrease in pertussis cases after the implementation of mitigation measures, with a reduction ranging from 28% to 87%. [16][17][18][19][20] Even if the incidence of pertussis has strongly decreased since the introduction of booster vaccines, B. pertussis is still circulating in France, especially in older children, adults and the elderly, as we have previously shown. 7,8,10 This may be due to non-optimal vaccination coverage in these populations, 21 and waning of immunity with the acellular pertussis vaccine. 7,22 The limited circulation of B. pertussis in the community due to the measures imposed by the COVID-19 pandemic is a clear positive collateral effect. However, this reduced bacterial circulation, raises concern about a potential decrease in herd immunity and the risk of a larger epidemic in the coming years. During the lockdown, pediatricians involved in data Source 3 declared following vaccinations guidelines. 23 However, as in many countries, delays in vaccination were observed in France in the general population since the beginning of the pandemic; they mostly concerned booster vaccinations for children and adults. 5,24 Such delays can be dangerous for infants under 6 months of age who are non-or partially vaccinated, notably if there is no catch-up before a potential resurgence of pertussis.
Indeed, as many other community-acquired infections, the dramatic decrease of pertussis incidence could suggest an epidemic rebound after relaxation of all mitigation measures. 25 Therefore, vaccination guidelines (infants, children, adolescents, adults and cocooning strategy) should be strictly followed to increase herd immunity, and maternal vaccination against pertussis during pregnancy, not yet implemented in France, should presently be considered by French national health institutions. 26 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 20, 2021. ; https://doi.org/10.1101/2021.07.16.21260367 doi: medRxiv preprint 1 3 Our study has limitations. First, measurement bias is present in Sources 1 and 2 as our pertussis case definition only considered positive PCRs. We therefore might have missed cases diagnosed by serology (sometimes prescribed even if not recommended), or based on clinical grounds, especially during the 1 st lockdown (from mid-March to May, 2020), when transportation was limited, office-based physicians unreachable, and private laboratories overwhelmed by the implementation of large-scale SARS-CoV-2 testing. Second, we cannot determine whether the decrease of pertussis cases observed here was due to decreased pertussis circulation or reduced testing. However, the similar decrease observed in hospitalized cases in the youngest population does not favor the latter hypothesis, and access retriction to outpatient testing was mostly limited to the first lockdown period. Third, we may have lacked statistical power to detect a significant decrease in the proportion of positive cases, as the post-lockdown period was relatively short compared to the pre-lockdown period (i.e., 9 months vs. 87 months, respectively), and negative PCR tests results were not available from data Source 2.
In conclusion, this national-level study shows a strong association between the COVID-19 pandemic and pertussis in France, with an unprecedented drop of pertussis cases.
Pertussis should be closely monitored to detect any resurgence in the community when social distancing restrictions will be relaxed, as well as compliance to vaccination of infants, children, and adults.

Data sharing
Anonymised database will be made available on reasonable request.

Authors license statement
This research was funded, in whole or in part, by Institut Pasteur and Santé publique France. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 20, 2021. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted July 20, 2021. ; https://doi.org/10.1101/2021.07.16.21260367 doi: medRxiv preprint 1 9 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. IRR: Incidence Rate Ratio; 1 the estimated IRR was adjusted for long-term cycles, seasonality (using pairs of sine/cosine terms) and secular trend; 2 the estimated change in level was adjusted for long-term cycles, seasonality (adjustment by calendar month), year, and secular trend. *Proportion calculated using Lab 1 data only for 2013-2015, and then data from Lab 1 and Lab 2 for 2016-2020, because of missing data on negative test results for Lab 2 for 2013-2015.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted July 20, 2021. ; https://doi.org/10.1101/2021.07.16.21260367 doi: medRxiv preprint