Genomic Surveillance and Evolutionary Dynamics of Influenza A Virus in Sri Lanka 

Genomic Surveillance and Evolutionary Dynamics of Influenza A Virus in Sri Lanka | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Genomic Surveillance and Evolutionary Dynamics of Influenza A Virus in Sri Lanka Tibutius Jayadas, Chandima Jeewandara, Bhagya Senadheera, Heshan Kuruppu, and 8 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4972640/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 26 Nov, 2024 Read the published version in Virology Journal → Version 1 posted 10 You are reading this latest preprint version Abstract Background Influenza A has been named as a priority pathogen by the WHO due to the potential to cause pandemics. Genomic sequencing of influenza strains is important to understand the evolution of the influenza strains and also to select the appropriate influenza vaccines to be used in the different influenza seasons in Sri Lanka. Therefore, we sought to understand the molecular epidemiology of the influenza viruses in the Western Province of Sri Lanka, including mutational analysis to investigate the evolutionary dynamics. Methods A total of 349 individuals presenting with fever and respiratory symptoms were enrolled in this study from November 2022 to May 2024. Nasopharyngeal and oropharyngeal specimens were collected and screened using quantitative PCR to detect Influenza A, Influenza B, and SARS-CoV-2. Subtyping and genomic sequencing was carried out on influenza A strains using Oxford Nanopore Technology. Results Influenza A was detected in 49 (14%) patients, influenza B in 20 (5.7%) and SARS-CoV-2 in 41 (11.7%). Co-infections were observed in five participants. The phylogenetic analysis assigned the H1N1 HA gene sequences within the 6B.1A.5a.2a clade. The HA gene of the H1N1 sequences in 2023 were assigned as belonging to the subclades C.1, C.1.2, and C.1.8, while the 2024 sequences were assigned to subclades C.1.8 and C.1.9. The H3N2 sequences from 2023 were assigned to the 3C.2a1b.2a.2a.1b clade and subclade G.1.1.2, while the 2024 sequences were assigned to the 3C.2a1b.2a.2a.3a.1 clade and subclade J.2. The K54Q, A186T, Q189E, E224A, R259K, K308R, I418V, and X215A amino acid substitutions were seen in the H1N1 in the 2023 and 2024 sequences. The 2024 H1N1 sequences additionally exhibited further substitutions, such as V47I, I96T, T120A, A139D, G339X, K156X, and T278S. Conclusion In this first study using genomic sequencing to characterize the influenza A strains in Sri Lanka, which showed different influenza A viruses circulating in an 18-month period. As the Sri Lankan strains also had certain mutations of unknown significance, it would be important to continue detailed surveillance of the influenza strains in Sri Lanka to choose the most suitable vaccines for the population and the timing of vaccine administration. Influenza A SARS-CoV-2 vaccines severe disease genomic sequencing clades subtypes molecular epidemiology Figures Figure 1 Figure 2 Figure 3 Figure 4 Introduction Seasonal influenza outbreaks are associated with significant morbidity and mortality, with estimated cases 3.2 million cases of severe disease each year, globally [ 1 ]. Due to the potential of influenza A strains causing pandemics, it has been included in the WHO pathogen prioritization list published in 2024 [ 2 ]. Despite the availability of effective vaccines and antivirals, the WHO estimates that 290,000 to 650,000 deaths occur annually due to this virus [ 3 ]. Those at extremes of age, pregnant women, individuals with comorbidities and immunocompromised individuals are at risk of developing severe disease and death [ 1 ]. Due to the rapid evolution of the virus and emergence of avian influenza in certain regions in the world, genomic surveillance of influenza strains in crucial to monitor the influenza strains that cause outbreaks in different countries. Influenza viruses belong to the Orthomyxoviridae family and are classified into four types: A, B, C, and D. Among these types, the influenza A virus has been responsible for several major pandemics in the past century, including the 1918 Spanish flu (H1N1), the 1957 Asian flu (H2N2), and the 1968 Hong Kong flu (H3N2), all of which caused a significant global health burden [ 4 ]. Influenza A viruses are categorized based on the properties of their surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA) [ 5 ]. The influenza A virus is classified into different subtypes based on the HA and NA glycoproteins. There are 18 known HA subtypes (H1 to H18) and 11 known NA subtypes (N1 to N11) [ 5 ]. The interplay of antigenic shift and drift among these subtypes results in generation of multiple strains due to varied combination of HA and NA subtypes [ 6 ]. Reassortment events, facilitated by natural reservoirs such as swine, birds, and horses, contribute to the emergence of novel strains with pandemic potential [ 6 ]. Globally, influenza A continues to exhibit seasonal patterns, with peaks typically occurring during the winter months in temperate regions and outbreaks often coinciding with the monsoon season in tropical and subtropical regions [ 7 ]. Although vaccination has been proven to be effective in providing some protection against influenza, they need to be given annually due to the changes in the circulating strains of the virus [ 8 ]. Therefore, the WHO Global Influenza Program recommends an evidence-based approach by grouping countries with similar seasonality patterns and virus antigenic characteristics into Influenza Vaccination Zones to address specific country needs [ 9 ]. The WHO encourages countries to conduct local surveillance to assess their seasonality patterns and circulating strains to facilitate the decision on selecting Northern Hemisphere (NH) and Southern Hemisphere (SH) vaccines and to determine the timing of vaccination campaigns [ 9 ]. Sri Lanka is a tropical country and influenza viruses circulate throughout the year, with two peaks typically occurring during the rainy seasons which are from May to July and from November to January[ 10 ]. As a part of the integrated SARS-CoV-2 and influenza surveillance platform, limited number of samples are subjected to testing for the presence of influenza A and SARS-CoV-2, which are then subjected to subtyping if influenza A is identified. However, genomic sequencing is not carried out, which is important to identify the origin and evolution of the influenza strains and also to select the appropriate influenza vaccines to be used in the different influenza seasons in Sri Lanka. In this study, we carried out proceed to understand the molecular epidemiology of the influenza viruses in the Western Province of Sri Lanka, including mutational analysis to investigate the evolutionary dynamics. Methods Recruitment of patients and collection of samples Nasopharyngeal and oropharyngeal specimens were collected from both 349 adult and paediatric patients presenting with an acute febrile illness with respiratory symptoms such as cough, sore throat, and rhinorrhea. Patients were recruited from two tertiary care hospitals, which were the National Institute of Infectious Disease and Colombo South Teaching Hospital, situated in the Western Province of Sri Lanka, between November 2022 to May 2024. Patients with a duration of illness of 7days were included in the study. Ethics approval Informed written consent was taken from all adult patients and in the case of paediatric patients, informed written consent was obtained from their parents/guardian. Ethics approval was obtained from the Ethics Review Committee, University of Sri Jayewardenepura. Screening for Respiratory Viruses Using Quantitative Polymerase Chain Reaction (qPCR) Viral RNA was extracted using Applied Biosystems™ MagMAX™ Viral/Pathogen Nucleic Acid Isolation Kit. All the collected samples were screened for Influenza A, Influenza B, using Respiratory Panel 1 qPCR Kit and when influenza A virus was detected it was also subtyped using the Viasure, Spain (VS-RPA112L v.03). Concurrently, each sample was tested for the presence of SARS-CoV-2 using TaqPath™ COVID-19 CE-IVD RT-PCR Kit (Thermo Fisher Scientific, USA). Library preparation and sequencing of the influenza A virus All 49 samples that were identified as being infected with influenza A by qPCR were chosen for ONT sequencing. The reaction mixture was prepared using 12.5 µL of Superscript III One-Step PCR reaction buffer, 0.5 µL of SuperScript III RT/Platinum Taq Mix (Thermo Fisher Scientific, USA), and primers (MBTuni-12 at 0.1 µM, MBTuni-12.4 at 0.1 µM, and MBTuni-13 at 0.2 µM). Additionally, 2.5 µL of RNA template was added, and PCR grade water was used to attain a final volume of 25 µL [ 11 ]. The PCR reactions were carried out with an initial incubation at 42°C for 60 minutes, followed by denaturation at 94°C for 2 minutes. This was succeeded by 5 cycles of denaturation at 94°C for 30 seconds, annealing at 45°C for 30 seconds, and extension at 68°C for 3 minutes. Subsequently, 20 cycles were performed with denaturation at 94°C for 30 seconds, annealing at 58°C for 30 seconds, and extension at 68°C for 3 minutes, concluding with a final extension at 68°C for 10 minutes. Libraries for sequencing were generated from the amplified samples using the ONT Rapid Barcoding Kit (SQK-RBK110.96), following the protocol version RBK_9126_v110_revO_24Mar2021. The pooled barcoded MinION library was subsequently loaded onto the MinION Mk1b sequencer from Oxford Nanopore Technologies, Oxford, United Kingdom, equipped with an R9.4 flow cell. Real-time base calling was performed using MinKNOW version 3.0.4 with the Guppy base calling software version 3.2.10. Generation of Consensus sequences (EPI2ME) Base calling was performed using the Guppy (version 6.5.7) with Fast model, 450 bps base calling model. The resulting reads were analyzed using the wf-flu workflow. Samples that were unclassified were excluded from further analysis. All samples that were successfully classified as Archetypes were subsequently submitted to the GISAID database. Construction of phylogenetic trees We used HA and NA genes sequences with > 90% coverage to create the phylogenetic trees. Accession numbers included in the analysis are included in Supplementary Table.1. From 2021 to 2024 sequences from WHO South-East Asian region (100 simple random samples), WHO global (100 random samples) and vaccine reference sequences in GISAID database was used to construct the phylogenetics trees for HA gene and NA gene. Phylogenetic analyzes for all IAV segments were performed. Sequence alignments were separately constructed for HA (H1 and H3 subtypes), NA (N1 and N2 subtypes). 222 sequences were included for analysis of H1, 209 for H3, 224 for N1 and 210 sequences for N2. Multiple sequence alignment was generated using MAFFT v.7.508 employing the FFT-NS-i algorithm. Subsequently, this multiple sequence alignment was used to infer a Randomized Axelerated Maximum Likelihood (RAxML) phylogenetic tree using RAxML (v.8.2.12) with GTRGAMMA substitution model and bootstrap of 1000 replicates. The best-fit model GTR + F + R5 was chosen using ModelFinder. Final visualizations of the phylogenetic tree were done using R\ggtree, R\ape and R\ggstar packages (R version 4.1.2). Mutational analysis Mutation analysis was carried out for the sequenced samples, prior to the variant calling, by removing the signal peptides in the H1 and H3 genes. To identify mutations in the H1N1 sequences, they were compared with the A/Wisconsin/588/2019 strain (EPI_ISL_19085699), which is the 2021–2022 Northern Hemisphere vaccine strain for H1NI [ 12 ]. To identify mutations in the Sri Lankan H3N2 sequences, they were compared with the A/Darwin/6/2021 (EPI_ISL_1563628) which was the 2022 Southern Hemisphere vaccine strain for H3N2[ 13 ]. The predicted position of the signal peptide in the sequences were identified with SignalP-5.0 tool. Based on the analysis using this predictive model, we identified that the predicted position of the signal peptide in the A/Wisconsin/588/2019 strain was in the positions in the amino acid positions, 1 to 17 (likelihood ratio 0.797) and for the A/Darwin/6/2021, amino acid positions 1 to 16 (likelihood ratio, 0.6971). Mutations were analyzed and visualized with R packages (R version 4.1.2) after removing the signal peptide region from the sequence of the protein. Results Of the 349 patients recruited in the study, 173 (49.5%) were males and 176 (50.4%) were females and 205 (58.7%) were adults. Influenza A was detected in 49 (14%) patients, influenza B in 20 (5.7%) and SARS-CoV-2 in 41 (11.7%). Co-infections were observed in five participants: four were co-infected with both Influenza A and B, and one individual was co-infected with Influenza A and SARS-CoV-2. The age distribution of these infections in different age groups is shown in Fig. 1 . Notably, the highest incidence of influenza A (42.8%) and influenza B (5.7%), was detected in children 60 years old, with 22.7% of the infections been detected in this age group, while 6/66 (9.1%) were infected with influenza A. Influenza A virus subtyping, and seasonal patterns of infection Infection due to influenza A, influenza B and SARS-CoV-2 was detected during the study period of November 2022 to May 2024. We paused the study during the months of July to October 2023, where very limited cases of respiratory infections were reported in both tertiary care hospitals. Of the individuals who tested positive for Influenza A, 23 identified as H1N1, 18 as H3N2, while 8 infections could not be classified. From December 2022 to February 2023, H1N1 was the predominant subtype of Influenza A (Supplementary Fig. 1). However, a significant shift occurred from early March 2023 to July 2023, with H3N2 becoming the dominant strain. By December 2023, a resurgence of the H1N1 subtype was observed. Phylogenetic Analysis of H1N1 viruses Out of the 49 influenza A samples, 21 were successfully sequenced, comprising 17 H1N1 samples and 4 H3N2 samples. Based on sequence quality, 14 HA genes and all 17 NA genes from the H1N1 viruses were included in the phylogenetic analysis, while 3 HA genes and 4 NA genes from the H3N2 viruses were analyzed. The phylogenetic analysis assigned the H1N1 HA gene sequences within the 6B.1A.5a.2a clade. The HA gene of the H1N1 sequences in 2023 were assigned as belonging to the subclades C.1, C.1.2, and C.1.8, while the 2024 sequences were assigned to subclades C.1.8 and C.1.9. Phylogenetic analysis of the H1N1 HA gene revealed that the 2023 sequences were most closely related to strains from Bangladesh and Bangkok, whereas the 2024 sequences were most similar to those from the Maldives (Fig. 2 ). The Sri Lankan H1N1 HA gene sequences and the A/Sydney/5/2021 Southern Hemisphere vaccine strain (used in the 2023 Southern influenza vaccine) belong to clade 6B.1A.5a.2a. Although A/Wisconsin/67/2022 and A/Victoria/4897/2022 from the Northern Hemisphere vaccine reference are in the 6B.1A.5a.2a.1 clade, the Sri Lankan HA gene sequences from 2023 and 2024 were more closely related to these Northern Hemisphere H1N1 strains than to the A/Sydney/5/2021 strain. The NA gene analysis showed that 2023 sequence showed close resemblance to the sequences from England and 2024 NA gene sequences closely resembled with sequences from Belgium and Bangladesh (Supplementary Fig. 2). Phylogenetic Analysis of H3N2 viruses The H3N2 sequences from 2023 were assigned to the 3C.2a1b.2a.2a.1b clade and subclade G.1.1.2, while the 2024 sequences were assigned to the 3C.2a1b.2a.2a.3a.1 clade and subclade J.2. HA gene analysis revealed that the 2023 sequences were closely related to those from Bangkok (Thailand), Cantabria (Spain), and England, whereas the 2024 sequences show similarity to those from Belgium and Nakhon Pathom (Thailand) (Fig. 3 ). The 2023 HA gene sequence is more closely related to the A/Darwin/6/2021 vaccine strain, while the 2024 HA gene sequence is more similar to A/Massachusetts/18/2022, both of which fall within the 3C.2a1b.2a.2a.3a.1 clade. NA gene analysis of the 2023 H3N2 samples indicated a close relationship with sequences from Catalonia (Spain), Rhode Island (USA) and England, while the 2024 samples showed the similarity to sequences from France (Supplementary Fig. 3). Mutation analysis The mutational analysis of the H1N1 hemagglutinin (HA) gene was carried out in reference to the A/Wisconsin/588/2019 (H1N1) strain. Accordingly, we identified amino acid substitutions, including K54Q, A186T, Q189E, E224A, R259K, K308R, I418V, and X215A across both the 2023 and 2024 sequences (Fig. 4 A). The 2024 H1N1 sequences additionally exhibited further substitutions, such as V47I, I96T, T120A, A139D, G339X, K156X, and T278S. The positions of these mutations and the function of these genes are shown in supplementary table 2. In the neuraminidase (NA) gene, H1N1 sequences identified in 2023 and 2024 shared the X136Q and V453M/V453T substitutions, with the 2024 sequences uniquely showing mutations at I264T, E433X, and E433K (Fig. 4 B). In comparison to the A/Wisconsin/67/2022 vaccine reference sequence, the HA gene of Sri Lankan H1N1 strains in 2023 and 2024, demonstrated substitutions including S137P, R142K, E260D, A277T, and D356T. The 2024 sequences also presented additional mutations, namely V47I, I96T, and T120A (Supplementary Fig. 4A). In the NA gene, the 2024 sequences revealed further substitutions at V13I, S200N, and L339S (Supplementary Fig. 4B). In the hemagglutinin (HA) gene of H3N2, using A/Darwin/6/2021 (H3N2) as the reference, the 2023 sequence revealed several amino acid substitutions, including I200V, M193I, N160D, E155G, R299K, D104G, and K276R. The I140K substitution was consistently observed across all analyzed sequences. In contrast, the 2024 sequences exhibited additional substitutions such as K189R, N49S, K276E, I260M, I223V, I192F, N122D, N96S, E50K, and G53N (Fig. 4 C). Due to the E50K and I223V substitutions, our H3N2 strains in 2024, are most similar to the A/Thailand/8/2022, subclade J). In comparison to the A/Massachusetts/18/2022 vaccine strain, both the 2023 and 2024 sequences shared the K276E/K276R substitutions, while the 2024 sequence uniquely exhibited the L86X substitution (Supplementary Fig. 4C). In the neuraminidase (NA) gene, the 2023 sequence displayed the D346G substitution, whereas the 2024 sequence showed additional substitutions, including M51I, I469T, R400K, S44X, and R150H (Fig. 4 D). Discussion In this study we have investigated the influenza strains circulating in the Western Province of Sri Lanka from 2023 to mid-2024, providing detailed analysis of the circulating clades of Influenza A. The frequency of both influenza A and influenza B was predominantly seen in children 60 years of age. Many studies have shown that individuals at extremes of age, including children, have shown to be vulnerable to be hospitalized due to influenza [ 14 ]. However, in our cohort SARS-CoV-2 accounted for most infections in those > 60 years of age (22.7%) compared to 9.1% of infections due to influenza A. Sri Lanka did not receive any COVID-19 vaccines as booster since 2022 [ 15 ] and therefore, elderly individuals and those with comorbidities are at increased risk of hospitalization due to COVID-19, possibly due to waning of immunity. In our cohort, 4 individuals had co-infection with influenza A and B, while one patient had co-infection with influenza and SARS-CoV-2. Co-infections with influenza A and B have been previously reported [ 16 – 18 ], and have shown to associate with a worse disease outcome [ 17 ]. We also reported one patient with co-infection with influenza and SAR-CoV-2, which has also previously been reported [ 19 ]. Due to the limited sample size, we could not determine if co-infections were associated with worse disease outcomes. Seasonal influenza outbreaks usually coinciding with the monsoon season in tropical and subtropical regions [ 7 ]. As reported in the Global Influenza Surveillance and Response System, of the WHO, a similar pattern is observed in the Western Province, Sri Lanka, where there are two influenza A seasons, which are from November to January and again from April to June [ 20 ]. During early 2023, the predominant influenza A subtype was H1N1, which was replaced by H3N2 as the predominant subtype by June 2023. In 2024, again H1N1 became the predominant subtype. These changes are consistent with the changes in the influenza A subtypes in India and Nepal, but different to the changes in subtypes seen in Bangladesh, Thailand and Bhutan [ 20 ]. In our study, all H1N1 sequences from 2023 and 2024 were classified within the 6B.1A.5a.2a clade. Our sequences, characterized by substitutions I418V and v47i, were placed within the C.1 subclade and its associated subclusters [ 21 ]. Similar H1N1 strains dominated in Southeast Asia, the Middle East, Africa, Central America, and parts of Europe [ 21 ]. The 5a.2a.1 clade, which has become more prevalent in the United States, Caribbean, Japan, and several European countries, marked by mutations like P137S and K142R, has significantly diverged from the 5a.2a clade in 2023 [ 21 ]. This antigenic drift resulted in reduced effectiveness of the 5a.2a-based vaccine, represented by the A/Sydney/5/2021 strain, leading the WHO to update the vaccine to target the 5a.2a.1 clade for the 2024 season, now represented by A/Wisconsin/67/2022 and A/Victoria/4897/2022 [ 13 ]. Our influenza A H1N1 strains in 2024 had the additional mutations I96T, T120A, A139D, G339X, K156X, and T278S. Although the positions and the function of these genes which carried these mutations are known, the significance of these mutations in relation to vaccine efficacy or virulence of the virus is not known. Therefore, it would be important to continue surveillance to understand if the influenza vaccine containing the strains of 5a.2a.1 provides protection against both 5a.2a and 5a.2a.1 viruses, currently circulating in Sri Lanka. Our phylogenetic analysis of the A(H3N2) HA gene sequences revealed the circulation of the 3C.2a1b.2a.2a.1b clade in 2023 and the 3C.2a1b.2a.2a.3a.1 clade in 2024. In 2023, the 3C.2a1b.2a.2 subclade, characterized by mutations such as I140K and K276R, was the most prevalence strain globally [ 20 ]. However, by 2024, the 3C.2a1b.2a.3 subclade, particularly the 2a.3a.1 lineage, emerged as the dominant strain [ 20 , 21 ]. Our 2024 sequences aligned with the .2a.3a.1 clade (clade J), marked by mutations such as K276E and V223I [ 21 ]. These changes led to significant antigenic drift, reducing the efficacy of the A/Darwin/9/2021-based vaccine, which was updated for 2024 including the A/Thailand/8/2022 and A/Massachusetts/18/ strains[ 21 ]. Our 2024 strains have additional mutations such as N122D and K276E. Although the effect of these mutations on vaccine efficacy is not clear, it would be important to continue surveillance to detect further emerging influenza strains. Conclusion In summary, in this study we have characterized the influenza A strains that circulated in Sri Lanka over a period of 18 months. We found that all H1N1 sequences from 2023 and 2024 were classified within the clade 6B.1A.5a.2a clade, while the H3N2 sequences in 2023 were assigned to clade 3C.2a1b.2a.2a.1b and the 2024 strains to clade 3C.2a1b.2a.2a.3a.1. As the Sri Lankan strains also had certain mutations of unknown significance, it would be important to continue detailed surveillance of the influenza strains in Sri Lanka to choose the most suitable vaccines for the population and the timing of vaccine administration. Declarations Ethics approval and consent to participate Ethics approval was obtained from the Ethics Review Committee of the University of Sri Jayewardenepura. All participants gave informed written consent. Consent for publication Not applicable. Availability of data and material All data is available within the manuscript and the supporting files. Competing interests The authors have no competing interests. Funding We are grateful to the NIH, USA (grant number 5U01AI151788-02) and the UK Medical Research Council for funding Authors' contributions Conceptualization: GNM, GSO, TTPJ Experiments and investigations: SD, MD, FB Data analysis: BS, SD, MD Project administration and supervision: GNM, CJ Funding acquisition: GNM, CJ, GSO Writing the manuscript: GNM, SD Reviewing the manuscript: GNM, GSO, CJ, BS Acknowledgements: also thankful to Asia Pathogen Genomic Initiative for providing Influenza sequencing primers. References Paget J, Staadegaard L, Wang X, Li Y, van Pomeren T, van Summeren J, et al. Global and national influenza-associated hospitalisation rates: Estimates for 40 countries and administrative regions. J Glob Health. 2023;13:04003. 10.7189/jogh.13.04003 . WHO: Pathogens Prioritization. A scientific framework for epidemic and pandemic research prepradeness. In. 2024: 60. WHO. Influenza (Seasonal). In.; 2023. Piret J, Boivin G. Pandemics Throughout History. Front Microbiol. 2020;11:631736. 10.3389/fmicb.2020.631736 . Uyeki TM, Hui DS, Zambon M, Wentworth DE, Monto AS, Influenza. Lancet. 2022;400(10353):693–706. 10.1016/S0140-6736(22)00982-5 . Wu NC, Wilson IA. Influenza Hemagglutinin Structures and Antibody Recognition. Cold Spring Harb Perspect Med. 2020;10(8). 10.1101/cshperspect.a038778 . Tamerius JD, Shaman J, Alonso WJ, Bloom-Feshbach K, Uejio CK, Comrie A, et al. Environmental predictors of seasonal influenza epidemics across temperate and tropical climates. PLoS Pathog. 2013;9(3):e1003194. 10.1371/journal.ppat.1003194 . Osterholm MT, Kelley NS, Sommer A, Belongia EA. Efficacy and effectiveness of influenza vaccines: a systematic review and meta-analysis. Lancet Infect Dis. 2012;12(1):36–44. 10.1016/S1473-3099(11)70295-X . WHO. Global Influenza Programme. https://www.who.int/teams/global-influenza-programme/surveillance-and-monitoring/influenza-surveillance-outputs (2024). Accessed 15th August 2024 2024. Rafeek RAM, Divarathna MVM, Noordeen F. History and current trends in influenza virus infections with special reference to Sri Lanka. Virusdisease. 2017;28(3):225–32. 10.1007/s13337-017-0390-8 . Zhou B, Donnelly ME, Scholes DT, St George K, Hatta M, Kawaoka Y, et al. Single-reaction genomic amplification accelerates sequencing and vaccine production for classical and Swine origin human influenza a viruses. J Virol. 2009;83(19):10309–13. 10.1128/JVI.01109-09 . WHO. Recommendations announced for influenza vaccine composition for the 2022–2023 northern hemisphere influenza season. In.; 2022. WHO. Recommended composition of influenza virus vaccines for use in the 2024 southern hemisphere influenza season. In.; 2022. Gilca R, Amini R, Carazo S, Doggui R, Frenette C, Boivin G, et al. The Changing Landscape of Respiratory Viruses Contributing to Hospitalizations in Quebec, Canada: Results From an Active Hospital-Based Surveillance Study. JMIR Public Health Surveill. 2024;10:e40792. 10.2196/40792 . Hannah Ritchie EO-O, Beltekian D, Mathieu E, Hasell J, Macdonald B. Charlie Giattino, Cameron Appel and Max Roser. Coronavirus (COVID-19) Vaccinations. https://ourworldindata.org/covid-vaccinations (2024). Accessed. He J, Hou S, Chen Y, Yu JL, Chen QQ, He L, et al. The Epidemiological Pattern and Co-infection of Influenza A and B by Surveillance Network From 2009 to 2014 in Anhui Province, China. Front Public Health. 2022;10:825645. 10.3389/fpubh.2022.825645 . Gregianini TS, Varella IRS, Fisch P, Martins LG, Veiga ABG. Dual and Triple Infections With Influenza A and B Viruses: A Case-Control Study in Southern Brazil. J Infect Dis. 2019;220(6):961–8. 10.1093/infdis/jiz221 . Miron VD, Banica L, Sandulescu O, Paraschiv S, Surleac M, Florea D, et al. Clinical and molecular epidemiology of influenza viruses from Romanian patients hospitalized during the 2019/20 season. PLoS ONE. 2021;16(11):e0258798. 10.1371/journal.pone.0258798 . Dadashi M, Khaleghnejad S, Abedi Elkhichi P, Goudarzi M, Goudarzi H, Taghavi A, et al. COVID-19 and Influenza Co-infection: A Systematic Review and Meta-Analysis. Front Med (Lausanne). 2021;8:681469. 10.3389/fmed.2021.681469 . WHO. Global Influenza Surveillance and Response System (GISRS). https://www.who.int/initiatives/global-influenza-surveillance-and-response-system (2024). Accessed 15th August 2024 2024. WHO: Influenza virus characterization: summary report, Europe, March 2024. In. Edited by Control CWROfEaSECfDPa. 2024: 29. Additional Declarations No competing interests reported. Supplementary Files Supplementarydata.docx Cite Share Download PDF Status: Published Journal Publication published 26 Nov, 2024 Read the published version in Virology Journal → Version 1 posted Editorial decision: Revision requested 04 Oct, 2024 Reviews received at journal 28 Sep, 2024 Reviews received at journal 19 Sep, 2024 Reviewers agreed at journal 16 Sep, 2024 Reviewers agreed at journal 14 Sep, 2024 Reviewers agreed at journal 07 Sep, 2024 Reviewers invited by journal 02 Sep, 2024 Editor assigned by journal 27 Aug, 2024 Submission checks completed at journal 27 Aug, 2024 First submitted to journal 25 Aug, 2024 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4972640","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":357627534,"identity":"a368c255-a616-4f6b-b1a9-8bb8e429ee52","order_by":0,"name":"Tibutius Jayadas","email":"","orcid":"","institution":"University of Sri Jayewardenepura","correspondingAuthor":false,"prefix":"","firstName":"Tibutius","middleName":"","lastName":"Jayadas","suffix":""},{"id":357627535,"identity":"9e3b3671-2f5c-421c-b571-a098833e0963","order_by":1,"name":"Chandima Jeewandara","email":"","orcid":"","institution":"University of Sri Jayewardenepura","correspondingAuthor":false,"prefix":"","firstName":"Chandima","middleName":"","lastName":"Jeewandara","suffix":""},{"id":357627536,"identity":"20aaf46d-e107-4570-a924-d468bbfe2368","order_by":2,"name":"Bhagya Senadheera","email":"","orcid":"","institution":"University of Sri Jayewardenepura","correspondingAuthor":false,"prefix":"","firstName":"Bhagya","middleName":"","lastName":"Senadheera","suffix":""},{"id":357627537,"identity":"f4b2f150-4fa0-4c89-996f-1081db20fc70","order_by":3,"name":"Heshan Kuruppu","email":"","orcid":"","institution":"University of Sri Jayewardenepura","correspondingAuthor":false,"prefix":"","firstName":"Heshan","middleName":"","lastName":"Kuruppu","suffix":""},{"id":357627538,"identity":"19388952-c107-4557-8ca6-428105628ebd","order_by":4,"name":"Rivindu Wickramanayaje","email":"","orcid":"","institution":"University of Sri Jayewardenepura","correspondingAuthor":false,"prefix":"","firstName":"Rivindu","middleName":"","lastName":"Wickramanayaje","suffix":""},{"id":357627539,"identity":"1e28aa74-17ba-431b-81fb-4ff9a09a1c2e","order_by":5,"name":"Farha Bary","email":"","orcid":"","institution":"University of Sri Jayewardenepura","correspondingAuthor":false,"prefix":"","firstName":"Farha","middleName":"","lastName":"Bary","suffix":""},{"id":357627540,"identity":"8dc9d30e-988d-43f7-bfa7-33bcfb429833","order_by":6,"name":"Ananda Wijewickrama","email":"","orcid":"","institution":"National Institute of Infectious Diseases, Sri Lanka","correspondingAuthor":false,"prefix":"","firstName":"Ananda","middleName":"","lastName":"Wijewickrama","suffix":""},{"id":357627541,"identity":"25e2c0f3-4ca2-4af6-97b2-0cb49ef4200e","order_by":7,"name":"Suranga Manilgama","email":"","orcid":"","institution":"National Institute of Infectious Diseases, Sri Lanka","correspondingAuthor":false,"prefix":"","firstName":"Suranga","middleName":"","lastName":"Manilgama","suffix":""},{"id":357627542,"identity":"2af7f777-2dc0-4815-860b-436945efcd15","order_by":8,"name":"Manouri Gamage","email":"","orcid":"","institution":"University of Sri Jayewardenepura","correspondingAuthor":false,"prefix":"","firstName":"Manouri","middleName":"","lastName":"Gamage","suffix":""},{"id":357627543,"identity":"4afecd2d-3c0d-4bfa-8d25-34729228e7a6","order_by":9,"name":"Nilanka Perera","email":"","orcid":"","institution":"University of Sri Jayewardenepura","correspondingAuthor":false,"prefix":"","firstName":"Nilanka","middleName":"","lastName":"Perera","suffix":""},{"id":357627544,"identity":"000a24ec-a788-408e-8899-56b00fbea843","order_by":10,"name":"Graham Ogg","email":"","orcid":"","institution":"University of Oxford","correspondingAuthor":false,"prefix":"","firstName":"Graham","middleName":"","lastName":"Ogg","suffix":""},{"id":357627545,"identity":"90d6c551-61b6-4d75-b73c-3ce0085c11b0","order_by":11,"name":"Gathsaurie MALAVIGE","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABQklEQVRIie2RMUvDQBTHXwi+LhezXoj0vsIrgUKpuPo1GgJ2tCCUQgUzxSU0q7tfwMk5cpAupV1bbmkpdOogCCWgiG0sUhIRR8H8hrt37/Hjf8cBlJT8TfBg7WSFNt8fM6hg6IcKZYVOv1IOx8h/UsT5JFmk1/LyuOIfPaR05g7M+6Tf6T6CyVvacwrSySm1kVdxWCIbAYtxFpLnBnx1oe7GCqww1q0QZD2v+B7agJKQt3DKSHeQj+rKCBTQ0AcbQJ7mlWiJVvq+VcQcZ290s1Umm6tMkaC/fqMI7iE3gl0KoGIkq2iGqO9TcJeSvxjxZd02Bm1C5gbqhIZV5OjYbKyYFT4FjZDa+eeLyF1Z6aZJ4lYms3Wvz0QkFy+sq6om9+Q07TVrfi4l/iq1/Yi3so19doofKfxCy4wLrZKSkpL/zQdf5mY6LLaAUQAAAABJRU5ErkJggg==","orcid":"","institution":"University of Sri Jayewardenepura","correspondingAuthor":true,"prefix":"","firstName":"Gathsaurie","middleName":"","lastName":"MALAVIGE","suffix":""}],"badges":[],"createdAt":"2024-08-25 12:22:38","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4972640/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4972640/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1186/s12985-024-02555-z","type":"published","date":"2024-11-26T15:57:27+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":66863462,"identity":"cdb37ddb-2158-4dca-b6b0-6a363d1f6acc","added_by":"auto","created_at":"2024-10-17 08:46:08","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":76709,"visible":true,"origin":"","legend":"\u003cp\u003eDistribution of influenza A, influenza B, and SARS-CoV-2 infections across different age groups among recruited participants\u003c/p\u003e","description":"","filename":"Figure1.png","url":"https://assets-eu.researchsquare.com/files/rs-4972640/v1/58bcfa301a825faa6bb5ce6f.png"},{"id":66863465,"identity":"cefdebc3-11a5-4177-81b2-fafa2acf8749","added_by":"auto","created_at":"2024-10-17 08:46:08","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":337539,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003ePhylogenetic tree of the H1N1 HA gene.\u003c/strong\u003e The phylogenetic tree was generated with the Sri Lankan H1N1 sequences (n=14) in comparison to the global H1N1 strains. All the Sri Lankan were assigned to clade 6B.1A.5a.2a. The H1N1 Sri Lankan sequence clusters are shaded in green, orange and grey shades, while the reference sequences are highlighted in green.\u003c/p\u003e","description":"","filename":"Figure2.png","url":"https://assets-eu.researchsquare.com/files/rs-4972640/v1/32183fa5bdcfb728b5c0f884.png"},{"id":66863463,"identity":"73964db9-a655-47e2-a97c-e02389e670f4","added_by":"auto","created_at":"2024-10-17 08:46:08","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":440450,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003ePhylogenetic tree of the H3N2 HA gene\u003c/strong\u003e. The phylogenetic tree was generated with the Sri Lankan H3N2 sequences (n=3) in comparison to the global H1N1 strains. The 2023 sequence was assigned to the 3C.2a1b.2a.2a.1b clade, while the 2024 sequences were assigned to the 3C.2a1b.2a.2a.3a.1 clade. The H1N1 Sri Lankan sequence clusters are shaded in green, orange and grey shades, while the reference sequences are highlighted in green.\u003c/p\u003e","description":"","filename":"Figure3.png","url":"https://assets-eu.researchsquare.com/files/rs-4972640/v1/c665fad472226fed43ad9916.png"},{"id":66864663,"identity":"3cf0211e-88d7-442f-9d97-c4ba1a2341dc","added_by":"auto","created_at":"2024-10-17 08:54:08","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":92841,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eHeatmap of amino acid substitutions in hemagglutinin (HA) and neuraminidase (NA) genes of influenza A H1N1 and H3N2 viruses\u003c/strong\u003e. (Panels A and B show mutations in the HA and NA genes of H1N1, respectively, while panels C and D display mutations in the HA and NA genes of H3N2. Each row represents an individual virus sequence, identified by its GISAID EPI_ISL accession number, and each column represents a specific amino acid position where mutations have occurred. The presence of a mutation is indicated by a blue square, and the absence by a white square).\u003c/p\u003e","description":"","filename":"Figure4.png","url":"https://assets-eu.researchsquare.com/files/rs-4972640/v1/f3cb176ccf94377957b97059.png"},{"id":70388579,"identity":"5409e4d0-2363-44fe-92e3-0be92454d559","added_by":"auto","created_at":"2024-12-02 17:26:32","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1298838,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4972640/v1/355fe9f2-90e7-43ee-babc-15d874cff98a.pdf"},{"id":66863466,"identity":"e64f9558-15a9-49ec-bc29-448f45a2d050","added_by":"auto","created_at":"2024-10-17 08:46:08","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":442272,"visible":true,"origin":"","legend":"","description":"","filename":"Supplementarydata.docx","url":"https://assets-eu.researchsquare.com/files/rs-4972640/v1/efa1328b9b57457de5a1dda9.docx"}],"financialInterests":"No competing interests reported.","formattedTitle":"Genomic Surveillance and Evolutionary Dynamics of Influenza A Virus in Sri Lanka ","fulltext":[{"header":"Introduction","content":"\u003cp\u003eSeasonal influenza outbreaks are associated with significant morbidity and mortality, with estimated cases 3.2\u0026nbsp;million cases of severe disease each year, globally [\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e]. Due to the potential of influenza A strains causing pandemics, it has been included in the WHO pathogen prioritization list published in 2024 [\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e]. Despite the availability of effective vaccines and antivirals, the WHO estimates that 290,000 to 650,000 deaths occur annually due to this virus [\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e]. Those at extremes of age, pregnant women, individuals with comorbidities and immunocompromised individuals are at risk of developing severe disease and death [\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e]. Due to the rapid evolution of the virus and emergence of avian influenza in certain regions in the world, genomic surveillance of influenza strains in crucial to monitor the influenza strains that cause outbreaks in different countries.\u003c/p\u003e \u003cp\u003eInfluenza viruses belong to the Orthomyxoviridae family and are classified into four types: A, B, C, and D. Among these types, the influenza A virus has been responsible for several major pandemics in the past century, including the 1918 Spanish flu (H1N1), the 1957 Asian flu (H2N2), and the 1968 Hong Kong flu (H3N2), all of which caused a significant global health burden [\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e]. Influenza A viruses are categorized based on the properties of their surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA) [\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e]. The influenza A virus is classified into different subtypes based on the HA and NA glycoproteins. There are 18 known HA subtypes (H1 to H18) and 11 known NA subtypes (N1 to N11) [\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e]. The interplay of antigenic shift and drift among these subtypes results in generation of multiple strains due to varied combination of HA and NA subtypes [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. Reassortment events, facilitated by natural reservoirs such as swine, birds, and horses, contribute to the emergence of novel strains with pandemic potential [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eGlobally, influenza A continues to exhibit seasonal patterns, with peaks typically occurring during the winter months in temperate regions and outbreaks often coinciding with the monsoon season in tropical and subtropical regions [\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e]. Although vaccination has been proven to be effective in providing some protection against influenza, they need to be given annually due to the changes in the circulating strains of the virus [\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e]. Therefore, the WHO Global Influenza Program recommends an evidence-based approach by grouping countries with similar seasonality patterns and virus antigenic characteristics into Influenza Vaccination Zones to address specific country needs [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e]. The WHO encourages countries to conduct local surveillance to assess their seasonality patterns and circulating strains to facilitate the decision on selecting Northern Hemisphere (NH) and Southern Hemisphere (SH) vaccines and to determine the timing of vaccination campaigns [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eSri Lanka is a tropical country and influenza viruses circulate throughout the year, with two peaks typically occurring during the rainy seasons which are from May to July and from November to January[\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e]. As a part of the integrated SARS-CoV-2 and influenza surveillance platform, limited number of samples are subjected to testing for the presence of influenza A and SARS-CoV-2, which are then subjected to subtyping if influenza A is identified. However, genomic sequencing is not carried out, which is important to identify the origin and evolution of the influenza strains and also to select the appropriate influenza vaccines to be used in the different influenza seasons in Sri Lanka. In this study, we carried out proceed to understand the molecular epidemiology of the influenza viruses in the Western Province of Sri Lanka, including mutational analysis to investigate the evolutionary dynamics.\u003c/p\u003e"},{"header":"Methods","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003eRecruitment of patients and collection of samples\u003c/h2\u003e \u003cp\u003eNasopharyngeal and oropharyngeal specimens were collected from both 349 adult and paediatric patients presenting with an acute febrile illness with respiratory symptoms such as cough, sore throat, and rhinorrhea. Patients were recruited from two tertiary care hospitals, which were the National Institute of Infectious Disease and Colombo South Teaching Hospital, situated in the Western Province of Sri Lanka, between November 2022 to May 2024. Patients with a duration of illness of 7days were included in the study.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003eEthics approval\u003c/h2\u003e \u003cp\u003eInformed written consent was taken from all adult patients and in the case of paediatric patients, informed written consent was obtained from their parents/guardian. Ethics approval was obtained from the Ethics Review Committee, University of Sri Jayewardenepura.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec5\" class=\"Section2\"\u003e \u003ch2\u003eScreening for Respiratory Viruses Using Quantitative Polymerase Chain Reaction (qPCR)\u003c/h2\u003e \u003cp\u003eViral RNA was extracted using Applied Biosystems\u0026trade; MagMAX\u0026trade; Viral/Pathogen Nucleic Acid Isolation Kit. All the collected samples were screened for Influenza A, Influenza B, using Respiratory Panel 1 qPCR Kit and when influenza A virus was detected it was also subtyped using the Viasure, Spain (VS-RPA112L v.03). Concurrently, each sample was tested for the presence of SARS-CoV-2 using TaqPath\u0026trade; COVID-19 CE-IVD RT-PCR Kit (Thermo Fisher Scientific, USA).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section2\"\u003e \u003ch2\u003eLibrary preparation and sequencing of the influenza A virus\u003c/h2\u003e \u003cp\u003eAll 49 samples that were identified as being infected with influenza A by qPCR were chosen for ONT sequencing. The reaction mixture was prepared using 12.5 \u0026micro;L of Superscript III One-Step PCR reaction buffer, 0.5 \u0026micro;L of SuperScript III RT/Platinum Taq Mix (Thermo Fisher Scientific, USA), and primers (MBTuni-12 at 0.1 \u0026micro;M, MBTuni-12.4 at 0.1 \u0026micro;M, and MBTuni-13 at 0.2 \u0026micro;M). Additionally, 2.5 \u0026micro;L of RNA template was added, and PCR grade water was used to attain a final volume of 25 \u0026micro;L [\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e]. The PCR reactions were carried out with an initial incubation at 42\u0026deg;C for 60 minutes, followed by denaturation at 94\u0026deg;C for 2 minutes. This was succeeded by 5 cycles of denaturation at 94\u0026deg;C for 30 seconds, annealing at 45\u0026deg;C for 30 seconds, and extension at 68\u0026deg;C for 3 minutes. Subsequently, 20 cycles were performed with denaturation at 94\u0026deg;C for 30 seconds, annealing at 58\u0026deg;C for 30 seconds, and extension at 68\u0026deg;C for 3 minutes, concluding with a final extension at 68\u0026deg;C for 10 minutes.\u003c/p\u003e \u003cp\u003eLibraries for sequencing were generated from the amplified samples using the ONT Rapid Barcoding Kit (SQK-RBK110.96), following the protocol version RBK_9126_v110_revO_24Mar2021. The pooled barcoded MinION library was subsequently loaded onto the MinION Mk1b sequencer from Oxford Nanopore Technologies, Oxford, United Kingdom, equipped with an R9.4 flow cell. Real-time base calling was performed using MinKNOW version 3.0.4 with the Guppy base calling software version 3.2.10.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section2\"\u003e \u003ch2\u003eGeneration of Consensus sequences (EPI2ME)\u003c/h2\u003e \u003cp\u003eBase calling was performed using the Guppy (version 6.5.7) with Fast model, 450 bps base calling model. The resulting reads were analyzed using the wf-flu workflow. Samples that were unclassified were excluded from further analysis. All samples that were successfully classified as Archetypes were subsequently submitted to the GISAID database.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003ch2\u003eConstruction of phylogenetic trees\u003c/h2\u003e \u003cp\u003eWe used HA and NA genes sequences with \u0026gt;\u0026thinsp;90% coverage to create the phylogenetic trees. Accession numbers included in the analysis are included in Supplementary Table.1. From 2021 to 2024 sequences from WHO South-East Asian region (100 simple random samples), WHO global (100 random samples) and vaccine reference sequences in GISAID database was used to construct the phylogenetics trees for HA gene and NA gene. Phylogenetic analyzes for all IAV segments were performed. Sequence alignments were separately constructed for HA (H1 and H3 subtypes), NA (N1 and N2 subtypes). 222 sequences were included for analysis of H1, 209 for H3, 224 for N1 and 210 sequences for N2.\u003c/p\u003e \u003cp\u003eMultiple sequence alignment was generated using MAFFT v.7.508 employing the FFT-NS-i algorithm. Subsequently, this multiple sequence alignment was used to infer a Randomized Axelerated Maximum Likelihood (RAxML) phylogenetic tree using RAxML (v.8.2.12) with GTRGAMMA substitution model and bootstrap of 1000 replicates. The best-fit model GTR\u0026thinsp;+\u0026thinsp;F\u0026thinsp;+\u0026thinsp;R5 was chosen using ModelFinder. Final visualizations of the phylogenetic tree were done using R\\ggtree, R\\ape and R\\ggstar packages (R version 4.1.2).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec9\" class=\"Section2\"\u003e \u003ch2\u003eMutational analysis\u003c/h2\u003e \u003cp\u003eMutation analysis was carried out for the sequenced samples, prior to the variant calling, by removing the signal peptides in the H1 and H3 genes. To identify mutations in the H1N1 sequences, they were compared with the A/Wisconsin/588/2019 strain (EPI_ISL_19085699), which is the 2021\u0026ndash;2022 Northern Hemisphere vaccine strain for H1NI [\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e]. To identify mutations in the Sri Lankan H3N2 sequences, they were compared with the A/Darwin/6/2021 (EPI_ISL_1563628) which was the 2022 Southern Hemisphere vaccine strain for H3N2[\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e]. The predicted position of the signal peptide in the sequences were identified with SignalP-5.0 tool. Based on the analysis using this predictive model, we identified that the predicted position of the signal peptide in the A/Wisconsin/588/2019 strain was in the positions in the amino acid positions, 1 to 17 (likelihood ratio 0.797) and for the A/Darwin/6/2021, amino acid positions 1 to 16 (likelihood ratio, 0.6971). Mutations were analyzed and visualized with R packages (R version 4.1.2) after removing the signal peptide region from the sequence of the protein.\u003c/p\u003e \u003c/div\u003e"},{"header":"Results","content":"\u003cp\u003eOf the 349 patients recruited in the study, 173 (49.5%) were males and 176 (50.4%) were females and 205 (58.7%) were adults. Influenza A was detected in 49 (14%) patients, influenza B in 20 (5.7%) and SARS-CoV-2 in 41 (11.7%). Co-infections were observed in five participants: four were co-infected with both Influenza A and B, and one individual was co-infected with Influenza A and SARS-CoV-2. The age distribution of these infections in different age groups is shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e. Notably, the highest incidence of influenza A (42.8%) and influenza B (5.7%), was detected in children\u0026thinsp;\u0026lt;\u0026thinsp;10 years of age. In contrast, the highest incidence of SARS-CoV-2 infection was seen in individuals\u0026thinsp;\u0026gt;\u0026thinsp;60 years old, with 22.7% of the infections been detected in this age group, while 6/66 (9.1%) were infected with influenza A.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cdiv id=\"Sec11\" class=\"Section2\"\u003e \u003ch2\u003eInfluenza A virus subtyping, and seasonal patterns of infection\u003c/h2\u003e \u003cp\u003eInfection due to influenza A, influenza B and SARS-CoV-2 was detected during the study period of November 2022 to May 2024. We paused the study during the months of July to October 2023, where very limited cases of respiratory infections were reported in both tertiary care hospitals. Of the individuals who tested positive for Influenza A, 23 identified as H1N1, 18 as H3N2, while 8 infections could not be classified. From December 2022 to February 2023, H1N1 was the predominant subtype of Influenza A (Supplementary Fig.\u0026nbsp;1). However, a significant shift occurred from early March 2023 to July 2023, with H3N2 becoming the dominant strain. By December 2023, a resurgence of the H1N1 subtype was observed.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec12\" class=\"Section2\"\u003e \u003ch2\u003ePhylogenetic Analysis of H1N1 viruses\u003c/h2\u003e \u003cp\u003eOut of the 49 influenza A samples, 21 were successfully sequenced, comprising 17 H1N1 samples and 4 H3N2 samples. Based on sequence quality, 14 HA genes and all 17 NA genes from the H1N1 viruses were included in the phylogenetic analysis, while 3 HA genes and 4 NA genes from the H3N2 viruses were analyzed. The phylogenetic analysis assigned the H1N1 HA gene sequences within the 6B.1A.5a.2a clade. The HA gene of the H1N1 sequences in 2023 were assigned as belonging to the subclades C.1, C.1.2, and C.1.8, while the 2024 sequences were assigned to subclades C.1.8 and C.1.9. Phylogenetic analysis of the H1N1 HA gene revealed that the 2023 sequences were most closely related to strains from Bangladesh and Bangkok, whereas the 2024 sequences were most similar to those from the Maldives (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThe Sri Lankan H1N1 HA gene sequences and the A/Sydney/5/2021 Southern Hemisphere vaccine strain (used in the 2023 Southern influenza vaccine) belong to clade 6B.1A.5a.2a. Although A/Wisconsin/67/2022 and A/Victoria/4897/2022 from the Northern Hemisphere vaccine reference are in the 6B.1A.5a.2a.1 clade, the Sri Lankan HA gene sequences from 2023 and 2024 were more closely related to these Northern Hemisphere H1N1 strains than to the A/Sydney/5/2021 strain. The NA gene analysis showed that 2023 sequence showed close resemblance to the sequences from England and 2024 NA gene sequences closely resembled with sequences from Belgium and Bangladesh (Supplementary Fig.\u0026nbsp;2).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec13\" class=\"Section2\"\u003e \u003ch2\u003ePhylogenetic Analysis of H3N2 viruses\u003c/h2\u003e \u003cp\u003eThe H3N2 sequences from 2023 were assigned to the 3C.2a1b.2a.2a.1b clade and subclade G.1.1.2, while the 2024 sequences were assigned to the 3C.2a1b.2a.2a.3a.1 clade and subclade J.2. HA gene analysis revealed that the 2023 sequences were closely related to those from Bangkok (Thailand), Cantabria (Spain), and England, whereas the 2024 sequences show similarity to those from Belgium and Nakhon Pathom (Thailand) (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). The 2023 HA gene sequence is more closely related to the A/Darwin/6/2021 vaccine strain, while the 2024 HA gene sequence is more similar to A/Massachusetts/18/2022, both of which fall within the 3C.2a1b.2a.2a.3a.1 clade. NA gene analysis of the 2023 H3N2 samples indicated a close relationship with sequences from Catalonia (Spain), Rhode Island (USA) and England, while the 2024 samples showed the similarity to sequences from France (Supplementary Fig.\u0026nbsp;3).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec14\" class=\"Section2\"\u003e \u003ch2\u003eMutation analysis\u003c/h2\u003e \u003cp\u003eThe mutational analysis of the H1N1 hemagglutinin (HA) gene was carried out in reference to the A/Wisconsin/588/2019 (H1N1) strain. Accordingly, we identified amino acid substitutions, including K54Q, A186T, Q189E, E224A, R259K, K308R, I418V, and X215A across both the 2023 and 2024 sequences (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eA). The 2024 H1N1 sequences additionally exhibited further substitutions, such as V47I, I96T, T120A, A139D, G339X, K156X, and T278S. The positions of these mutations and the function of these genes are shown in supplementary table 2. In the neuraminidase (NA) gene, H1N1 sequences identified in 2023 and 2024 shared the X136Q and V453M/V453T substitutions, with the 2024 sequences uniquely showing mutations at I264T, E433X, and E433K (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eB). In comparison to the A/Wisconsin/67/2022 vaccine reference sequence, the HA gene of Sri Lankan H1N1 strains in 2023 and 2024, demonstrated substitutions including S137P, R142K, E260D, A277T, and D356T. The 2024 sequences also presented additional mutations, namely V47I, I96T, and T120A (Supplementary Fig.\u0026nbsp;4A). In the NA gene, the 2024 sequences revealed further substitutions at V13I, S200N, and L339S (Supplementary Fig.\u0026nbsp;4B).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eIn the hemagglutinin (HA) gene of H3N2, using A/Darwin/6/2021 (H3N2) as the reference, the 2023 sequence revealed several amino acid substitutions, including I200V, M193I, N160D, E155G, R299K, D104G, and K276R. The I140K substitution was consistently observed across all analyzed sequences. In contrast, the 2024 sequences exhibited additional substitutions such as K189R, N49S, K276E, I260M, I223V, I192F, N122D, N96S, E50K, and G53N (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eC). Due to the E50K and I223V substitutions, our H3N2 strains in 2024, are most similar to the A/Thailand/8/2022, subclade J). In comparison to the A/Massachusetts/18/2022 vaccine strain, both the 2023 and 2024 sequences shared the K276E/K276R substitutions, while the 2024 sequence uniquely exhibited the L86X substitution (Supplementary Fig.\u0026nbsp;4C). In the neuraminidase (NA) gene, the 2023 sequence displayed the D346G substitution, whereas the 2024 sequence showed additional substitutions, including M51I, I469T, R400K, S44X, and R150H (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eD).\u003c/p\u003e \u003c/div\u003e"},{"header":"Discussion","content":"\u003cp\u003eIn this study we have investigated the influenza strains circulating in the Western Province of Sri Lanka from 2023 to mid-2024, providing detailed analysis of the circulating clades of Influenza A. The frequency of both influenza A and influenza B was predominantly seen in children\u0026thinsp;\u0026lt;\u0026thinsp;10 years of age while SARS-CoV-2 infection was seen in adults\u0026thinsp;\u0026gt;\u0026thinsp;60 years of age. Many studies have shown that individuals at extremes of age, including children, have shown to be vulnerable to be hospitalized due to influenza [\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e]. However, in our cohort SARS-CoV-2 accounted for most infections in those\u0026thinsp;\u0026gt;\u0026thinsp;60 years of age (22.7%) compared to 9.1% of infections due to influenza A. Sri Lanka did not receive any COVID-19 vaccines as booster since 2022 [\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e] and therefore, elderly individuals and those with comorbidities are at increased risk of hospitalization due to COVID-19, possibly due to waning of immunity. In our cohort, 4 individuals had co-infection with influenza A and B, while one patient had co-infection with influenza and SARS-CoV-2. Co-infections with influenza A and B have been previously reported [\u003cspan additionalcitationids=\"CR17\" citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e], and have shown to associate with a worse disease outcome [\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e]. We also reported one patient with co-infection with influenza and SAR-CoV-2, which has also previously been reported [\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e]. Due to the limited sample size, we could not determine if co-infections were associated with worse disease outcomes.\u003c/p\u003e \u003cp\u003eSeasonal influenza outbreaks usually coinciding with the monsoon season in tropical and subtropical regions [\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e]. As reported in the Global Influenza Surveillance and Response System, of the WHO, a similar pattern is observed in the Western Province, Sri Lanka, where there are two influenza A seasons, which are from November to January and again from April to June [\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e]. During early 2023, the predominant influenza A subtype was H1N1, which was replaced by H3N2 as the predominant subtype by June 2023. In 2024, again H1N1 became the predominant subtype. These changes are consistent with the changes in the influenza A subtypes in India and Nepal, but different to the changes in subtypes seen in Bangladesh, Thailand and Bhutan [\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eIn our study, all H1N1 sequences from 2023 and 2024 were classified within the 6B.1A.5a.2a clade. Our sequences, characterized by substitutions I418V and v47i, were placed within the C.1 subclade and its associated subclusters [\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e]. Similar H1N1 strains dominated in Southeast Asia, the Middle East, Africa, Central America, and parts of Europe [\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e]. The 5a.2a.1 clade, which has become more prevalent in the United States, Caribbean, Japan, and several European countries, marked by mutations like P137S and K142R, has significantly diverged from the 5a.2a clade in 2023 [\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e]. This antigenic drift resulted in reduced effectiveness of the 5a.2a-based vaccine, represented by the A/Sydney/5/2021 strain, leading the WHO to update the vaccine to target the 5a.2a.1 clade for the 2024 season, now represented by A/Wisconsin/67/2022 and A/Victoria/4897/2022 [\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e]. Our influenza A H1N1 strains in 2024 had the additional mutations I96T, T120A, A139D, G339X, K156X, and T278S. Although the positions and the function of these genes which carried these mutations are known, the significance of these mutations in relation to vaccine efficacy or virulence of the virus is not known. Therefore, it would be important to continue surveillance to understand if the influenza vaccine containing the strains of 5a.2a.1 provides protection against both 5a.2a and 5a.2a.1 viruses, currently circulating in Sri Lanka.\u003c/p\u003e \u003cp\u003eOur phylogenetic analysis of the A(H3N2) HA gene sequences revealed the circulation of the 3C.2a1b.2a.2a.1b clade in 2023 and the 3C.2a1b.2a.2a.3a.1 clade in 2024. In 2023, the 3C.2a1b.2a.2 subclade, characterized by mutations such as I140K and K276R, was the most prevalence strain globally [\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e]. However, by 2024, the 3C.2a1b.2a.3 subclade, particularly the 2a.3a.1 lineage, emerged as the dominant strain [\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e, \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e]. Our 2024 sequences aligned with the .2a.3a.1 clade (clade J), marked by mutations such as K276E and V223I [\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e]. These changes led to significant antigenic drift, reducing the efficacy of the A/Darwin/9/2021-based vaccine, which was updated for 2024 including the A/Thailand/8/2022 and A/Massachusetts/18/ strains[\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e]. Our 2024 strains have additional mutations such as N122D and K276E. Although the effect of these mutations on vaccine efficacy is not clear, it would be important to continue surveillance to detect further emerging influenza strains.\u003c/p\u003e"},{"header":"Conclusion","content":"\u003cp\u003eIn summary, in this study we have characterized the influenza A strains that circulated in Sri Lanka over a period of 18 months. We found that all H1N1 sequences from 2023 and 2024 were classified within the clade 6B.1A.5a.2a clade, while the H3N2 sequences in 2023 were assigned to clade 3C.2a1b.2a.2a.1b and the 2024 strains to clade 3C.2a1b.2a.2a.3a.1. As the Sri Lankan strains also had certain mutations of unknown significance, it would be important to continue detailed surveillance of the influenza strains in Sri Lanka to choose the most suitable vaccines for the population and the timing of vaccine administration.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eEthics approval was obtained from the Ethics Review Committee of the University of Sri Jayewardenepura. All participants gave informed written consent.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsent for publication\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable. \u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and material\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAll data is available within the manuscript and the supporting files.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors have no competing interests.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe are grateful to the NIH, USA (grant number 5U01AI151788-02) and the UK Medical Research Council for funding\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthors\u0026apos; contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eConceptualization: GNM, GSO, TTPJ\u003c/p\u003e\n\u003cp\u003eExperiments and investigations: SD, MD, FB\u003c/p\u003e\n\u003cp\u003eData analysis: BS, SD, MD\u003c/p\u003e\n\u003cp\u003eProject administration and supervision: GNM, CJ\u003c/p\u003e\n\u003cp\u003eFunding acquisition: GNM, CJ, GSO\u003c/p\u003e\n\u003cp\u003eWriting the manuscript: GNM, SD\u003c/p\u003e\n\u003cp\u003eReviewing the manuscript: GNM, GSO, CJ, BS\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgements:\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003ealso thankful to Asia Pathogen Genomic Initiative for providing Influenza sequencing primers.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003ePaget J, Staadegaard L, Wang X, Li Y, van Pomeren T, van Summeren J, et al. Global and national influenza-associated hospitalisation rates: Estimates for 40 countries and administrative regions. J Glob Health. 2023;13:04003. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.7189/jogh.13.04003\u003c/span\u003e\u003cspan address=\"10.7189/jogh.13.04003\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWHO: Pathogens Prioritization. A scientific framework for epidemic and pandemic research prepradeness. In. 2024: 60.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWHO. Influenza (Seasonal). In.; 2023.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePiret J, Boivin G. Pandemics Throughout History. Front Microbiol. 2020;11:631736. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.3389/fmicb.2020.631736\u003c/span\u003e\u003cspan address=\"10.3389/fmicb.2020.631736\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eUyeki TM, Hui DS, Zambon M, Wentworth DE, Monto AS, Influenza. Lancet. 2022;400(10353):693\u0026ndash;706. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/S0140-6736(22)00982-5\u003c/span\u003e\u003cspan address=\"10.1016/S0140-6736(22)00982-5\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWu NC, Wilson IA. Influenza Hemagglutinin Structures and Antibody Recognition. Cold Spring Harb Perspect Med. 2020;10(8). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1101/cshperspect.a038778\u003c/span\u003e\u003cspan address=\"10.1101/cshperspect.a038778\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTamerius JD, Shaman J, Alonso WJ, Bloom-Feshbach K, Uejio CK, Comrie A, et al. Environmental predictors of seasonal influenza epidemics across temperate and tropical climates. PLoS Pathog. 2013;9(3):e1003194. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1371/journal.ppat.1003194\u003c/span\u003e\u003cspan address=\"10.1371/journal.ppat.1003194\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eOsterholm MT, Kelley NS, Sommer A, Belongia EA. Efficacy and effectiveness of influenza vaccines: a systematic review and meta-analysis. Lancet Infect Dis. 2012;12(1):36\u0026ndash;44. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/S1473-3099(11)70295-X\u003c/span\u003e\u003cspan address=\"10.1016/S1473-3099(11)70295-X\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWHO. Global Influenza Programme. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.who.int/teams/global-influenza-programme/surveillance-and-monitoring/influenza-surveillance-outputs\u003c/span\u003e\u003cspan address=\"https://www.who.int/teams/global-influenza-programme/surveillance-and-monitoring/influenza-surveillance-outputs\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e (2024). Accessed 15th August 2024 2024.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRafeek RAM, Divarathna MVM, Noordeen F. History and current trends in influenza virus infections with special reference to Sri Lanka. Virusdisease. 2017;28(3):225\u0026ndash;32. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1007/s13337-017-0390-8\u003c/span\u003e\u003cspan address=\"10.1007/s13337-017-0390-8\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhou B, Donnelly ME, Scholes DT, St George K, Hatta M, Kawaoka Y, et al. Single-reaction genomic amplification accelerates sequencing and vaccine production for classical and Swine origin human influenza a viruses. J Virol. 2009;83(19):10309\u0026ndash;13. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1128/JVI.01109-09\u003c/span\u003e\u003cspan address=\"10.1128/JVI.01109-09\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWHO. Recommendations announced for influenza vaccine composition for the 2022\u0026ndash;2023 northern hemisphere influenza season. In.; 2022.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWHO. Recommended composition of influenza virus vaccines for use in the 2024 southern hemisphere influenza season. In.; 2022.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGilca R, Amini R, Carazo S, Doggui R, Frenette C, Boivin G, et al. The Changing Landscape of Respiratory Viruses Contributing to Hospitalizations in Quebec, Canada: Results From an Active Hospital-Based Surveillance Study. JMIR Public Health Surveill. 2024;10:e40792. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.2196/40792\u003c/span\u003e\u003cspan address=\"10.2196/40792\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHannah Ritchie EO-O, Beltekian D, Mathieu E, Hasell J, Macdonald B. Charlie Giattino, Cameron Appel and Max Roser. Coronavirus (COVID-19) Vaccinations. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://ourworldindata.org/covid-vaccinations\u003c/span\u003e\u003cspan address=\"https://ourworldindata.org/covid-vaccinations\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e (2024). Accessed.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHe J, Hou S, Chen Y, Yu JL, Chen QQ, He L, et al. The Epidemiological Pattern and Co-infection of Influenza A and B by Surveillance Network From 2009 to 2014 in Anhui Province, China. Front Public Health. 2022;10:825645. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.3389/fpubh.2022.825645\u003c/span\u003e\u003cspan address=\"10.3389/fpubh.2022.825645\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGregianini TS, Varella IRS, Fisch P, Martins LG, Veiga ABG. Dual and Triple Infections With Influenza A and B Viruses: A Case-Control Study in Southern Brazil. J Infect Dis. 2019;220(6):961\u0026ndash;8. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1093/infdis/jiz221\u003c/span\u003e\u003cspan address=\"10.1093/infdis/jiz221\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMiron VD, Banica L, Sandulescu O, Paraschiv S, Surleac M, Florea D, et al. Clinical and molecular epidemiology of influenza viruses from Romanian patients hospitalized during the 2019/20 season. PLoS ONE. 2021;16(11):e0258798. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1371/journal.pone.0258798\u003c/span\u003e\u003cspan address=\"10.1371/journal.pone.0258798\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDadashi M, Khaleghnejad S, Abedi Elkhichi P, Goudarzi M, Goudarzi H, Taghavi A, et al. COVID-19 and Influenza Co-infection: A Systematic Review and Meta-Analysis. Front Med (Lausanne). 2021;8:681469. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.3389/fmed.2021.681469\u003c/span\u003e\u003cspan address=\"10.3389/fmed.2021.681469\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWHO. Global Influenza Surveillance and Response System (GISRS). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.who.int/initiatives/global-influenza-surveillance-and-response-system\u003c/span\u003e\u003cspan address=\"https://www.who.int/initiatives/global-influenza-surveillance-and-response-system\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e (2024). Accessed 15th August 2024 2024.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWHO: Influenza virus characterization: summary report, Europe, March 2024. In. Edited by Control CWROfEaSECfDPa. 2024: 29.\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"virology-journal","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"virj","sideBox":"Learn more about [Virology Journal](http://virologyj.biomedcentral.com/)","snPcode":"12985","submissionUrl":"https://submission.nature.com/new-submission/12985/3","title":"Virology Journal","twitterHandle":"@VirologyJ","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"BMC/SO AJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"Influenza A, SARS-CoV-2, vaccines, severe disease, genomic sequencing, clades, subtypes, molecular epidemiology","lastPublishedDoi":"10.21203/rs.3.rs-4972640/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4972640/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003ch2\u003eBackground\u003c/h2\u003e \u003cp\u003eInfluenza A has been named as a priority pathogen by the WHO due to the potential to cause pandemics. Genomic sequencing of influenza strains is important to understand the evolution of the influenza strains and also to select the appropriate influenza vaccines to be used in the different influenza seasons in Sri Lanka. Therefore, we sought to understand the molecular epidemiology of the influenza viruses in the Western Province of Sri Lanka, including mutational analysis to investigate the evolutionary dynamics.\u003c/p\u003e\u003ch2\u003eMethods\u003c/h2\u003e \u003cp\u003eA total of 349 individuals presenting with fever and respiratory symptoms were enrolled in this study from November 2022 to May 2024. Nasopharyngeal and oropharyngeal specimens were collected and screened using quantitative PCR to detect Influenza A, Influenza B, and SARS-CoV-2. Subtyping and genomic sequencing was carried out on influenza A strains using Oxford Nanopore Technology.\u003c/p\u003e\u003ch2\u003eResults\u003c/h2\u003e \u003cp\u003eInfluenza A was detected in 49 (14%) patients, influenza B in 20 (5.7%) and SARS-CoV-2 in 41 (11.7%). Co-infections were observed in five participants. The phylogenetic analysis assigned the H1N1 HA gene sequences within the 6B.1A.5a.2a clade. The HA gene of the H1N1 sequences in 2023 were assigned as belonging to the subclades C.1, C.1.2, and C.1.8, while the 2024 sequences were assigned to subclades C.1.8 and C.1.9. The H3N2 sequences from 2023 were assigned to the 3C.2a1b.2a.2a.1b clade and subclade G.1.1.2, while the 2024 sequences were assigned to the 3C.2a1b.2a.2a.3a.1 clade and subclade J.2. The K54Q, A186T, Q189E, E224A, R259K, K308R, I418V, and X215A amino acid substitutions were seen in the H1N1 in the 2023 and 2024 sequences. The 2024 H1N1 sequences additionally exhibited further substitutions, such as V47I, I96T, T120A, A139D, G339X, K156X, and T278S.\u003c/p\u003e\u003ch2\u003eConclusion\u003c/h2\u003e \u003cp\u003eIn this first study using genomic sequencing to characterize the influenza A strains in Sri Lanka, which showed different influenza A viruses circulating in an 18-month period. As the Sri Lankan strains also had certain mutations of unknown significance, it would be important to continue detailed surveillance of the influenza strains in Sri Lanka to choose the most suitable vaccines for the population and the timing of vaccine administration.\u003c/p\u003e","manuscriptTitle":"Genomic Surveillance and Evolutionary Dynamics of Influenza A Virus in Sri Lanka ","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-10-17 08:46:04","doi":"10.21203/rs.3.rs-4972640/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2024-10-04T09:07:03+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2024-09-28T12:43:28+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2024-09-19T06:46:12+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"95167001048026703132943217686871789979","date":"2024-09-16T13:13:07+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"147746276093035669157109367281511640325","date":"2024-09-15T03:31:28+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"168715046027511354379861055463688935042","date":"2024-09-08T01:52:34+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2024-09-03T01:08:37+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2024-08-27T10:18:26+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2024-08-27T10:04:39+00:00","index":"","fulltext":""},{"type":"submitted","content":"Virology Journal","date":"2024-08-25T12:21:05+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"virology-journal","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"virj","sideBox":"Learn more about [Virology Journal](http://virologyj.biomedcentral.com/)","snPcode":"12985","submissionUrl":"https://submission.nature.com/new-submission/12985/3","title":"Virology Journal","twitterHandle":"@VirologyJ","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"BMC/SO AJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"2fbf2967-87be-478e-8da8-fdc7d3a8d702","owner":[],"postedDate":"October 17th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[],"tags":[],"updatedAt":"2024-12-02T17:21:26+00:00","versionOfRecord":{"articleIdentity":"rs-4972640","link":"https://doi.org/10.1186/s12985-024-02555-z","journal":{"identity":"virology-journal","isVorOnly":false,"title":"Virology Journal"},"publishedOn":"2024-11-26 15:57:27","publishedOnDateReadable":"November 26th, 2024"},"versionCreatedAt":"2024-10-17 08:46:04","video":"","vorDoi":"10.1186/s12985-024-02555-z","vorDoiUrl":"https://doi.org/10.1186/s12985-024-02555-z","workflowStages":[]},"version":"v1","identity":"rs-4972640","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4972640","identity":"rs-4972640","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}