The impact of COVID-19 pandemic on influenza transmission: molecular and epidemiological evidence ================================================================================================= * Leon King Tran * Dai-Wei Huang * Nien-Kung Li * Lucy M. Li * Julia A. Palacios * Hsiao-Han Chang ## Abstract To quantify the impact of COVID-19-related control measures on the spread of human influenza virus, we analyzed case numbers, viral molecular sequences, personal behavior data, and policy stringency data from various countries, and found consistent evidence of decrease in influenza incidence after the emergence of COVID-19. **Article Summary** We quantify a noticeable decrease in H1N1 and H3N2 cases and genetic diversity in selected countries since the onset of the COVID-19 pandemic. Keywords * COVID-19 * influenza * case * effective population size * genetic diversity * epidemiology ## Introduction The emergence and spread of COVID-19 in 2020 led to a number of large-scale public health measures to limit international travel, reduce gatherings, and increase mask-wearing. While these preventative measures were implemented to curtail the COVID-19 pandemic, they seem to have also impacted the spread of other respiratory illnesses. There have been several reports on the decrease in case numbers during 2019-2020 influenza season in the northern hemisphere [6], and the lack of a 2020 influenza season in the southern hemisphere [8]. Here, we quantify the impact of the COVID-19 pandemic on the spread of influenza in terms of incidence and viral molecular diversity [1]. ## Method ### Case count data We analyzed weekly case count data of influenza available in FluNet [10] from various regions during the 2010-2020 influenza seasons. We defined *T**S* and *T**E* as the weeks during which the estimated number of cases reached 10% and 90% of the total case numbers in each influenza season, respectively. Since the influenza outbreaks for most regions started before the COVID-19 outbreak, we compared *T**E* and durations of influenza seasons pre- and post-COVID-19 pandemic. We defined the duration of an influenza season by the difference between *T**S* and *T**E*, and standardized the duration in the 2019-2020 season by the average and standard deviation of the duration from previous 9 seasons. ### Viral molecular sequences We analyzed the HA segment of human influenza A H1N1 and H3N2 sequences available in the GISAID EpiFlu database on November 1st 2020. The collection dates of the sequences ranged from January 2016 to December 2020. We used BEAST [3] to estimate the effective population size (Ne) from 2016 to 2020 for each location. The numbers of sequences analyzed are indicated in Appendix Table 1. ## Results and Discussion To examine the indirect impact of COVID-19 on influenza dynamics, we compared the 2019-2020 influenza season with previous 9 seasons in 21 locations across 5 continents (Table 1). We found that for all locations in Asia, the outbreak in 2019-2020 both ended earlier and lasted for a shorter duration than previous years. In the Americas, Europe, and Africa, 12 locations out of 15 showed an earlier end of the flu season in the 2019-2020 season than previous years; the rest remained similar. For locations where influenza seasons usually end later in the year, such as Brazil, Guatemala, and South Africa, the difference in duration between the 2019-2020 season and previous seasons was larger than other locations. The flu season in the Southern Hemisphere, which usually starts much later in the year, disappeared in several countries in 2020 [6,8]. In addition to case count data, we analyzed molecular data to evidence our findings. For each location and each influenza type, we calculated the within-location genetic diversity, Watterson’s θ [9], for the first half of 2019 and for the first of 2020 (Appendix Table 2). We found that for 10 out of 11 locations we analyzed, θ decreased from 2019 to 2020 for H1N1, and 9 out of 11 locations for H3N2. On the other hand, we calculated between-region genetic diversity for each pair of regions once in 2019 and again in 2020; 9 out of 11 locations for H1N1 and 6 out of 11 locations for H3N2 had their between-region genetic diversity increase from 2019 to 2020 in at least 50% of pairs for which the location was involved, reflecting reduced travel between regions in 2020 (Appendix Table 3). View this table: [Table 1.](http://medrxiv.org/content/early/2021/06/10/2021.06.08.21258434/T1) Table 1. Comparison of influenza outbreaks before and after the emergence of COVID-19 While Watterson’s θ measures overall viral diversity, the effective population size (Ne) quantifies genetic diversity over time [4]. We estimated Ne for H1N1 in 11 countries and found a decrease in Ne in 9 countries, including Italy and Taiwan (Figure 1; Appendix Figure 1; Appendix Table 1). For H3N2, we analyzed 5 countries, and found a decrease in Ne in South Africa and Taiwan in 2020 (Figure 1; Appendix Figure 2; Appendix Table 1). ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/10/2021.06.08.21258434/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2021/06/10/2021.06.08.21258434/F1) Figure 1. Effective Population Size Trajectory of H1N1 in Italy and H3N2 in Taiwan For personal measures and government policies against COVID-19, we also noticed that Asian countries tend to act earlier than countries in other continents, especially wearing masks (Appendix Table 4). Taken together, we observed earlier ends of flu seasons in Asia than in Europe and America, which could be explained by the earlier implementation of non-pharmaceutical interventions. The decrease in influenza incidence after the emergence of COVID-19 was confirmed by both case count and molecular data, providing stronger support for a genuine decrease in influenza incidence, even when count data might be incomplete or imprecise due to potentially less flu surveillance in 2020. ## Data Availability The flu sequences are found on GISAID. The accession numbers can be found at https://github.com/leonkt/accession-numbers. [https://github.com/leonkt/accession-numbers](https://github.com/leonkt/accession-numbers) [https://www.gisaid.org/](https://www.gisaid.org/) ## First Author Biography Leon King Tran is an undergraduate student in Mathematics and master’s student in Statistics at Stanford University and is interested in Biostatistics and Stochastic Processes. Dai-Wei Huang is an undergraduate student in the Interdisciplinary Program of Life Science, National Tsing Hua University, Hsinchu, Taiwan and is currently doing research related to epidemiology and biostatistics. ## Methods ### Personal measures and governmental policies against COVID-19 We analyzed two personal measures taken against the spread of COVID-19 from YouGov [7]– wearing masks when in public places and improving personal hygiene, from February 21, 2020 to December 17, 2020. For governmental responses to COVID-19, we used the data from OxCGRT, the Oxford COVID-19 Government Response Tracker [5], from January 1, 2020 to May 31, 2020. ### Molecular Epidemiology We estimated a mutation rate of H3N2 for New York State of 4 x 10−3 per base per year and 3.6 x 10−3 for H1N1. We then fixed the corresponding global mutation rates for the rest of the locations. We assumed the HKY prior mutation model with empirical frequencies and Gamma+Invariant site heterogeneity with 4 Gamma categories and 3 partitions into codon positions and assumed the Bayesian Skyride Coalescent prior on effective population size trajectories. We examined convergence by the likelihood trace plot and the effective sample size for each location and type. To estimate Ne, we divided the output (Ne * Gt) by the generation time (Gt) for H1N1 (2.3 days) and for H3N2 (3.1 days) [2]. We ran the MCMC for 10,000,000 iterations, thinning every 1000 iterations and with 10% of burnin. We removed any regions without convergence. The within-region diversities were measured by Watterson’s θ. The between-region diversities were measured as follows: Between-region diversity between region 1 and region ![Graphic][1], where *θ**pooled* was Watterson’s θ when combining two regions together. View this table: [Appendix Table 1.](http://medrxiv.org/content/early/2021/06/10/2021.06.08.21258434/T2) Appendix Table 1. Effective sample sizes of root age and number of sequences per country for influenza types H1N1 and H3N2. These effective sample size numbers were obtained with Tracer [3]. Effective sample sizes with N/A indicate lack of convergence and were omitted from the analysis. View this table: [Appendix Table 2.](http://medrxiv.org/content/early/2021/06/10/2021.06.08.21258434/T3) Appendix Table 2. Within-region genetic diversity during the first half of 2019 compared to the first half of 2020, for influenza types H1N1 and H3N2. Diversity was measured by Watterson’s θ. The decrease column indicates whether genetic diversity decreased from the first half of 2019 to the first half of 2020. Sequences with more than 5% gaps were removed from the analysis. N/A indicates a lack of data. View this table: [Appendix Table 3.](http://medrxiv.org/content/early/2021/06/10/2021.06.08.21258434/T4) Appendix Table 3. Between-region diversity during the first half of 2019 compared to the first half of 2020, for influenza types H1N1 and H3N2. For each pair of regions, we calculated the between-region diversity in the first half of 2019 and in the first half of 2020. The third column for each region denotes the number of pairs involving that region with an increase in between-region diversity from 2019 to 2020, divided by the total number of pairwise comparisons. Regions where this fraction is above 0.5 are marked in grey. N/A indicates that there was not enough data to calculate diversity. View this table: [Appendix Table 4.](http://medrxiv.org/content/early/2021/06/10/2021.06.08.21258434/T5) Appendix Table 4. Comparison of personal measures and governmental policies against COVID-19 ![Appendix Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/10/2021.06.08.21258434/F2.medium.gif) [Appendix Figure 1.](http://medrxiv.org/content/early/2021/06/10/2021.06.08.21258434/F2) Appendix Figure 1. Estimated Effective Population Size Trajectory for H1N1 in selected regions. ![Appendix Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/10/2021.06.08.21258434/F3.medium.gif) [Appendix Figure 2.](http://medrxiv.org/content/early/2021/06/10/2021.06.08.21258434/F3) Appendix Figure 2. Comparison of Effective Population Size Trajectory for H3N2 in selected regions. In addition, we would like to thank the following institutions for providing the sequences obtained from GISAID: Yokohama City Institute of Public Health., National Institute of Animal Health, China Center For Disease Control And Prevention, Akita Research Center for Public Health and Environment, Ningbo international travel healthcare center, University of Corsica, Tottori Prefectural Institute of Public Health and Environmental Science, Ishikawa Prefectural Institute of Public Health and Environmental Science, Zhejiang Provincial Center for Disease Control and Prevention, Research Institute for Environmental Sciences and Public Health of Iwate Prefecture, Nara Prefectural Institute for Hygiene and Environment, Royal Darwin Hospital, Kawasaki City Institute of Public Health, Wuhan Institute of Virology, Tochigi Prefectural Institute of Public Health and Environmental Science, National Institute for Biological Standards and Control (NIBSC), Kyoto City Institute of Health and Environmental Sciences, UOC Policlinico di Bari DIMO, California Department of Health Services, Laboratorio Aziendale di Microbiologia e Virologia, Cabrini Pathology, Laboratorio Central do Estado do Rio de Janeiro - LACEN, Okinawa Prefectural Institute of Health and Environment, Yamaguchi Prefectural Institute of Public Health and Environment, First Affiliated Hospital Medical School of Zhejiang University, New York State Department of Health, Hubei Provincial center for disease control and prevention, Shimane Prefectural Institute of Public Health and Environmental Science, Yamanashi Institute for Public Health, Beijing Center for Disease Prevention and Control, Kumamoto Prefectural Institute of Public Health and Environmental Science, National Institute of Biological Standards and Control, Dorevtich Pathology, Fondazione IRCCS Policlinico San Matteo, Shiga Prefectural Institute of Public Health, Austin Health, National Institute of Infectious Diseases (NIID), Instituto Adolfo Lutz, State Key Laboratory of Virology and Wuhan Institute of Virology, Chinese Academy of Sciences, Royal Melbourne Hospital, Fukui Prefectural Institute of Public Health, Kobe Institute of Health, Institute of Medical and Veterinary Science (IMVS), Prince of Wales Hospital, Laboratorio Nacional De Salud Guatemala, Melbourne Pathology, Hyogo Prefectural Institute of Public Health and Consumer Sciences, Osaka City Institute of Public Health and Environmental Sciences, Saitama Institute of Public Health, Universita degli Studi di Palermo, Dipartimento di Scienze per la Promozione della Salute, Ehime Prefecture Institute of Public Health and Environmental Science, Wakayama Prefectural Research Center of Environment and Public Health, Aichi Prefectural Institute of Public Health, WHO Collaborating Centre for Reference and Research on Influenza, Kagawa Prefectural Research Institute for Environmental Sciences and Public Health, Fukushima Prefectural Institute of Public Health, Sapporo City Institute of Public Health, Amedeo di Savoia Hospital, LACEN, AOU Meyer, Instituto Oswaldo Cruz FIOCRUZ - Laboratory of Respiratory Viruses and Measles (LVRS), Nagasaki Prefectural Institute for Environment Research and Public Health, University of Sassari, Tokushima Prefectural Centre for Public Health and Environmental Sciences, John Hunter Hospital, Friedrich-Loeffler-Institut, New York City Department of Health, Instituto de Salud Publica de Chile, WHO Chinese National Influenza Center, The Kyodoken Institute for Animal Science Research, Universita Cattolica del Sacro Cuore, Centers for Disease Control and Prevention, University of Florence, National Institute for Medical Research, Childrens Hospital Westmead, Saitama City Institute of Health Science and Research, Alfred Hospital, University of Trieste, Chiba Prefectural Institute of Public Health, Okayama Prefectural Institute for Environmental Science and Public Health, Laboratorio di Virologia - Azienda Ospedaliero-Universitaria Ospedali Riuniti - Ancona, Miyazaki Prefectural Institute for Public Health and Environment, Osaka Prefectural Institute of Public Health, Niigata City Institute of Public Health and Environment, Fukuoka Institute of Public Health and Environmental Sciences, Hiroshima City Institute of Public Health, University of Padova, Kochi Public Health and Sanitation Institute, Utsunomiya City Institute of Public Health and Environment Science, Beijing Mentougou District Center for Diseases Prevention and Control, FUNED MG, Hobart Pathology, University of Siena, Virus Research Center, Sendai Medical Center, Istituto Superiore di Sanita, Ningbo International Travel Healthcare Center, University of Milan, Kanagawa Prefectural Institute of Public Health, San Diego County Public Health Lab, Kyoto Prefectural Institute of Public Health and Environment, Ibaraki Prefectural Institute of Public Health, Nagano Environmental Conservation Research Institute, Gunma Prefectural Institute of Public Health and Environmental Sciences, Shanghai International Travel Healthcare Center, Saitama medical university, Chiba City Institute of Health and Environment, Istituto Superiore di Sanita (ISS), Niigata University (IDRC), CNR Virus des Infections Respiratoires - France SUD, Microbiology Services Colindale, Public Health England, The First Affiliated Hospital, College of Medicine, Zhejiang University, Niigata University (DPH), Sakai City Institute of Public Health, Tokyo Metropolitan Institute of Public Health, Institut Pasteur, Universita di Corsica Pasquale Paoli, Center for Public Health and Environment, Hiroshima Prefectural Technology Research Institute, Canberra Hospital, Clinical Virology Unit, CDIM, Seqirus Pty Ltd (CSL Group), National Institute for Communicable Diseases, WHO Centre for Reference, Beijing Institute of Microbiology and Epidemiology, Institut Pasteur of Shanghai, Victorian Infectious Diseases Reference Laboratory, University of Perugia, Mentougou District center for diseases prevention and control, Oita Prefectural Institute of Health and Environment, National Institute for Communicable Disease, The University of Hong Kong, Institut Pasteur of Shanghai, CAS, Laboratorio Central do Estado do Parana - LACEN, Queensland Health Forensic and Scientific Services, Shizuoka City Institute of Environmental Sciences and Public Health, Kumamoto City Environmental Research Center, Royal Chidrens Hospital, Amagasaki City Institute of Public Health, University of Genoa, Hubei Provincial Center for Disease Control and Prevention, CSL Ltd, icddr,b International Centre for Diarrhoeal Disease Research, Bangladesh, Laboratorio Central de Saude Publica Professor Goncalo Moniz, LACEN-BA, Evandro Chagas Institute, Laboratorio Central de Saude Publica (LACEN RJ), Kagoshima University, Fukuoka City Institute for Hygiene and the Environment, Guangzhou No.8 People, Lismore Base Hospital, Pathwest QE II Medical Centre, State Key Laboratory of Reipiratory Disease (Guangzhou Medical University), Sagamihara City Laboratory of Public Health, University of Parma, South China Agricultural University, Crick Worldwide Influenza Centre, Saga Prefectural Institute of Public Health and Pharmaceutical Research, Yamagata Prefectural Institute of Public Health, LACEN SC, Influenza Surveillance Centre for Disease Control, National Health Laboratory, Saitama Medical University, Kagoshima Prefectural Institute for Environmental Research and Public Health, Microbiology University Politecnica delle Marche, Mie Prefecture Health and Environment Research Institute, Institute of Epidemiology Disease Control and Research (IEDCR), Sendai City Institute of Public Health, Gifu Prefectural Institute of Health and Environmental Sciences, Guangdong Provincial Center for Disease Control and Prevention, Douglass Hanly Moir Pathology, Nagano City Health Center, Westmead Hospital, Shizuoka Institute of Environment and Hygiene, National Influenza Center, Yokosuka Institute of Public Health, Guangzhou Institute of Respiratory Diseases (GIRD), Royal Hobart Hospital, Miyagi Prefectural Institute of Public Health and Environment, Niigata Prefectural Institute of Public Health and Environmental Sciences, Hamamatsu City Health Environment Research Center, Aomori Prefectural Institute of Public Health and Environment, University of Study of Siena, South China Agricultural University Veterinary Medicine College, AO dei Colli Monaldi - Cotugno, U.O.C. Microbiologia e Virologia, Laboratorio Central de Saude Publica de Alagoas, LACEN-AL, Monash Medical Centre, Canterbury Health Services, Wakayama City Institute of Public Health, University of Bari, Universita degli Studi di Padova, Vanderbilt University Medical Center, U.S. Air Force School of Aerospace Medicine, Hokkaido Institute of Public Health, Kitakyusyu City Institute of Enviromental Sciences, Gifu Municipal Institute of Public Health, New York Medical College, and Toyama Institute of Health. GISAID accession numbers for all sequences used can be found at [https://github.com/leonkt/accession-numbers/](https://github.com/leonkt/accession-numbers/). ## Acknowledgements We thank GISAID and sequencing laboratories for making the influenza A viral molecular sequences available Specific laboratories and accession numbers are found at the end of the Appendix. HHC and DWH were supported by the Ministry of Science and Technology in Taiwan (MOST 110-2636-B-007-009); HHC was supported by the Yushan Scholar Program. The funders had no role in study design, data collection, data analysis, data interpretation, or writing of the report. ## Appendix * Received June 8, 2021. * Revision received June 8, 2021. * Accepted June 10, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. [1].Bedford, T., Riley, S., Barr, I. et al. Global circulation patterns of seasonal influenza viruses vary with antigenic drift. Nature 523, 217–220 (2015). [https://doi.org/10.1038/nature14460](https://doi.org/10.1038/nature14460) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature14460&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26053121&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F10%2F2021.06.08.21258434.atom) 2. [2].Carrat, F., Vergu, E., Ferguson, NM., Lemaitre, M., Cauchemez, S., Leach, S., Valleron, AJ. Timelines of infection and disease in human influenza: a review of volunteer challenge studies. Am J Epidemiol. 2008 Apr 1;167(7):775–85. [https://doi.org/10.1093/aje/kwm375](https://doi.org/10.1093/aje/kwm375) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwm375&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18230677&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F10%2F2021.06.08.21258434.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000254469500002&link_type=ISI) 3. [3].Drummond, A.J., Rambaut, A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol 7, 214 (2007). [https://doi.org/10.1186/1471-2148-7-214](https://doi.org/10.1186/1471-2148-7-214) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2148-7-214&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17996036&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F10%2F2021.06.08.21258434.atom) 4. [4].Frost, S. and Volz, EM. Viral phylodynamics and the search for an ‘effective number of infections’ Phil. Trans. R. Soc. [http://doi.org/10.1098/rstb.2010.0060](http://doi.org/10.1098/rstb.2010.0060) 5. [5].Hale, T., Angrist, N., Goldszmidt, R. et al. A global panel database of pandemic policies (Oxford COVID-19 Government Response Tracker). Nat Hum Behav 5, 529–538 (2021). [https://doi.org/10.1038/s41562-021-01079-8](https://doi.org/10.1038/s41562-021-01079-8) 6. [6].Kuo, S., Shih, S., Chien, L., & Hsiung, C. A. (2020). Collateral Benefit of COVID-19 Control Measures on Influenza Activity, Taiwan. Emerging Infectious Diseases, 26(8), 1928–1930. [https://doi.org/10.3201/eid2608.201192](https://doi.org/10.3201/eid2608.201192) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2608.201192&link_type=DOI) 7. [7].“Personal Measures Taken to Avoid COVID-19.” YouGov 17 Mar. 2020, [https://yougov.co.uk/topics/international/articles-reports/2020/03/17/personal-measures-taken-avoid-covid-19](https://yougov.co.uk/topics/international/articles-reports/2020/03/17/personal-measures-taken-avoid-covid-19). 8. [8].“The Southern Hemisphere Skipped Flu Season in 2020.” The Economist, The Economist Newspaper, [www.economist.com/graphic-detail/2020/09/12/the-southern-hemisphere-skipped-flu-season-in-2020](http://www.economist.com/graphic-detail/2020/09/12/the-southern-hemisphere-skipped-flu-season-in-2020) 9. [9].Watterson G. A. (1975). On the number of segregating sites in genetical models without recombination. Theoretical population biology, 7(2), 256–276. [https://doi.org/10.1016/0040-5809(75)90020-9](https://doi.org/10.1016/0040-5809(75)90020-9) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0040-5809(75)90020-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1145509&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F10%2F2021.06.08.21258434.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1975W295600011&link_type=ISI) 10. [10].“WHO FluNet.” World Health Organization, 2021, [https://www.who.int/influenza/en/](https://www.who.int/influenza/en/). [1]: /embed/inline-graphic-1.gif