Evolutionary Dynamics of Whole-Genome Influenza A/H3N2 Viruses Isolated in Myanmar from 2015 to 2019

Molecular epidemiological surveillance of respiratory syncytial virus infection in Myanmar from 2019 to 2023

To evaluate genetic changes in respiratory syncytial virus (RSV) between 2019 and 2023, we analyzed RSV strains from Myanmar before and after the COVID- 19 pandemic. Real-time polymerase chain reaction (RT-PCR) positive samples from children presenting with acute respiratory infections at outpatient clinics in Yangon were sequenced to determine the genotype. Phylogenetic and molecular evolutionary analyses were conducted using the Bayesian Markov Chain Monte Carlo method to construct the time-scale Maximum Clade Credibility tree. Of 1127 samples, 104 (9.2%) RSV-A and 233 (20.7%) RSV-B were positive by RT-PCR. There was an absence of a notable epidemic in 2020, a temporal shift with a surge of RSV-A in the 2021 outbreak, a lack of expected cases in 2022 and a substantial resurgence of RSV-B in 2023. The genotype of RSV-A was mainly A.D.3 lineage through the study period, while RSV-B were B.D.4.1.1 and B.D.E.1. RSV-A showed that the same lineage persisted within Myanmar throughout the pandemic, leading to a large outbreak post-COVID. In contrast, RSV-B strains appear to have temporarily disappeared during the pandemic, but subsequently, globally circulating strains likely entered Myanmar, resulting in a major outbreak in 2023. The estimated evolutionary rate at the G-ectodomain for RSV-A was 7.76 × 10⁻³ and RSV-B was 5.67 × 10⁻³ substitutions/site/year. Strengthening genomic surveillance will likely support comparisons of circulating strains with those in other countries and facilitate the introduction of vaccines and other interventions.

The Ongoing Epidemics of Seasonal Influenza A(H3N2) in Hangzhou, China, and Its Viral Genetic Diversity

This study examined the genetic and evolutionary features of influenza A/H3N2 viruses in Hangzhou (2010-2022) by analyzing 28,651 influenza-like illness samples from two sentinel hospitals. Influenza A/H3N2 coexisted with other subtypes, dominating seasonal peaks (notably summer). Whole-genome sequencing of 367 strains was performed on GridION platforms. Phylogenetic analysis showed they fell into 16 genetic groups, with multiple clades circulating simultaneously. Shannon entropy indicated HA, NA, and NS gene segments exhibited significantly higher variability than other genomic segments, with HA glycoprotein mutations concentrated in antigenic epitopes A-E. Antiviral resistance showed no inhibitor resistance mutations in PA, PB1, or PB2, but NA mutations were detected in some strains, and most strains harbored M2 mutations. A Bayesian molecular clock showed the HA segment exhibited the highest nucleotide substitution rate (3.96 × 10-3 substitutions/site/year), followed by NA (3.77 × 10-3) and NS (3.65 × 10-3). Selective pressure showed A/H3N2 strains were predominantly under purifying selection, with only sporadic positive selection at specific sites. The Pepitope model demonstrated that antigenic epitope mismatches between circulating H3N2 variants and vaccine strains led to a significant decline in influenza vaccine effectiveness (VE), particularly in 2022. Overall, the study underscores the complex circulation patterns of influenza in Hangzhou and the global importance of timely vaccine strain updates.

Reducing amplification cycles to improve the coverage of influenza A virus genome sequencing in heterosubtypic co-infection

This study delineates the enhancement of a Reverse Transcription Polymerase Chain Reaction (RT-PCR) method for the amplification of the complete genome of the influenza A virus during heterosubtypic co-infection, relying on the amplification of intact gene segments. The precision of the method was assessed using all amplicons, which underwent both capillary electrophoresis and DNA sequencing. Five samples featuring co-infection of Influenza A viruses with H1N1 and H3N2 subtypes were evaluated. The improved strategy successfully amplified all eight segments of H3N2 strains in four samples, and the entire genome of H1N1 strains in three samples.

Genetic and antigenic characterization of influenza A(H3N2) virus after 13 consecutive years of influenza surveillance in Senegal, 2010-2022

Despite decades of influenza surveillance in many African countries, little is known about the evolutionary dynamics of seasonal influenza viruses. This study aimed to characterize the epidemiological, genetic and antigenic profiles of A/H3N2 viruses in Senegal from 2010 to 2022. A/H3N2 infection was confirmed using reverse transcription-polymerase chain reaction. Subsequently, a representative of A/H3N2 isolates was selected for genome sequencing. Predicted vaccine efficacy was measured using the Pepitope model. During the study period, 22638 samples were tested and influenza was detected in 31.8%, among which type A was confirmed in 78.1%. Of the Influenza A cases, the H3N2 subtype was detected in 29.8%, peaking at expected times during the rainy season. Genome sequencing of 123A/H3N2 isolates yielded 24 complete and 99 partial genomic sequences. Phylogenetic analysis revealed the circulation of multiple clades of A/H3N2 in Senegal, including 2a.3, 3C.2 and 3C.3a. A/H3N2 isolates were mainly susceptible to the influenza antiviral drugs oseltamivir and zanamivir, but the primary adamantine-resistance marker, S31N was encountered in all isolates. At least nine potential N-linked glycosylation sites were predicted among A/H3N2 strains, six of which (at positions 24, 38, 79, 181, 262 and 301) remains conserved among all isolates. Antigenic distances between circulating strains and vaccine viruses indicated varying vaccine efficacies, from suboptimal to moderate protection. The findings emphasize the need to enhance local genomic and antigenic surveillance and further research on influenza epidemiology and genetic evolution in sub-Saharan Africa.

Evolutionary Insights from Association Rule Mining of Co-Occurring Mutations in Influenza Hemagglutinin and Neuraminidase

Seasonal influenza viruses continuously evolve via antigenic drift. This leads to recurring epidemics, globally significant mortality rates, and the need for annually updated vaccines. Co-occurring mutations in hemagglutinin (HA) and neuraminidase (NA) are suggested to have synergistic interactions where mutations can increase the chances of immune escape and viral fitness. Association rule mining was used to identify temporal relationships of co-occurring HA-NA mutations of influenza virus A/H3N2 and its role in antigenic evolution. A total of 64 clusters were found. These included well-known mutations responsible for antigenic drift, as well as previously undiscovered groups. A majority (41/64) were associated with known antigenic sites, and 38/64 involved mutations across both HA and NA. The emergence and disappearance of N-glycosylation sites in the pattern of N-X-[S/T] were also identified, which are crucial post-translational processes to maintain protein stability and functional balance (e.g., emergence of NA:339ASP and disappearance of HA:187ASP). Our study offers an alternative approach to the existing mutual-information and phylogenetic methods used to identify co-occurring mutations, enabling faster processing of large amounts of data. Our approach can facilitate the prediction of critical mutations given their occurrence in a previous season, facilitating vaccine development for the next flu season and leading to better preparation for future pandemics.

Whole-Genome Analysis of the Influenza A(H1N1)pdm09 Viruses Isolated from Influenza-like Illness Outpatients in Myanmar and Community-Acquired Oseltamivir-Resistant Strains Present from 2015 to 2019

In this study, we describe the genetic characteristics of influenza A(H1N1)pdm09 strains detected in Myanmar from 2015 to 2019. Whole genomes from 60 A(H1N1)pdm09 virus isolates were amplified using real-time polymerase chain reaction and successfully sequenced using the Illumina iSeq100 platforms. Eight individual phylogenetic trees were retrieved for each segment along with those of the World Health Organization (WHO)-recommended Southern Hemisphere vaccine strains for the respective years. A(H1N1)pdm09 viruses from 2015 were found to belong to clade 6B, those from 2016 to 6B.1, 2017 to 6B.1A, and 2019 to 6B.1A.5a, and were genetically distinct from the Southern Hemisphere vaccine strains for the respective seasons, A/California/7/2009 and A/Michigan/45/2015. We observed one virus with intra-subtype reassortment, collected in the 2015 season. Importantly, three viruses possessed the H275Y substitution in the neuraminidase protein, appearing to be community-acquired without the prior administration of neuraminidase inhibitors. These viruses exhibited highly reduced susceptibility to oseltamivir and peramivir. This study demonstrates the importance of monitoring genetic variations in influenza viruses that will contribute to the selection of global influenza vaccines.

Characteristics and evolution of hemagglutinin and neuraminidase genes of Influenza A(H3N2) viruses in Thailand during 2015 to 2018

Background

Influenza A(H3N2) virus evolves continuously. Its hemagglutinin (HA) and neuraminidase (NA) genes have high genetic variation due to the antigenic drift. This study aimed to investigate the characteristics and evolution of HA and NA genes of the influenza A(H3N2) virus in Thailand.

Methods

Influenza A positive respiratory samples from 2015 to 2018 were subtyped by multiplex real-time RT-PCR. Full-length HA and NA genes from the positive samples of influenza A(H3N2) were amplified and sequenced. Phylogenetic analysis with the maximum likelihood method was used to investigate the evolution of the virus compared with the WHO-recommended influenza vaccine strain. Homology modeling and N-glycosylation site prediction were also performed.

Results

Out of 443 samples, 147 (33.18%) were A(H1N1)pdm09 and 296 (66.82%) were A(H3N2). The A(H3N2) viruses circulating in 2015 were clade 3C.2a whereas sub-clade 3C.2a1 and 3C.2a2 dominated in 2016-2017 and 2018, respectively. Amino acid substitutions were found in all antigenic sites A, B, C, D, and E of HA but the majority of the substitutions were located at antigenic sites A and B. The S245N and N329S substitutions in the NA gene affect the N-glycosylation. None of the mutations associated with resistance to NA inhibitors were observed. Mean evolutionary rates of the HA and NA genes were 3.47 × 10 -3 and 2.98 × 10-3 substitutions per site per year.

Conclusion

The influenza A(H3N2) virus is very genetically diverse and is always evolving to evade host defenses. The HA and NA gene features including the evolutionary rate of the influenza A(H3N2) viruses that were circulating in Thailand between 2015 and 2018 are described. This information is useful for monitoring the genetic characteristics and evolution in HA and NA genes of influenza A(H3N2) virus in Thailand which is crucial for predicting the influenza vaccine strains resulting in high vaccine effectiveness.

Genomic Analyses Uncover Evolutionary Features of Influenza A/H3N2 Viruses in Yunnan Province, China, from 2017 to 2022

Influenza A viruses evolve at a high rate of nucleotide substitution, thereby requiring continuous monitoring to determine the efficacy of vaccines and antiviral drugs. In the current study, we performed whole-genome sequencing analyses of 253 influenza A/H3N2 strains from Yunnan Province, China, during 2017-2022. The hemagglutinin (HA) segments of Yunnan A/H3N2 strains isolated during 2017-2018 harbored a high genetic diversity due to heterogeneous distribution across branches. The mutation regularity of the predominant antigenic epitopes of HA segments in Yunnan was inconsistent in different years. Some important functional mutations in gene segments associated with viral adaptation and drug tolerance were revealed. The rapid genomic evolution of Yunnan A/H3N2 strains from 2017 to 2022 mainly concentrated on segments, i.e., matrix protein 2 (M2), non-structural protein 1 (NS1), neuraminidase (NA), NS2, and HA, with a high overall non-synonymous/synonymous substitution ratio (dN/dS). Our results highlighted a decline in vaccine efficacy against the A/H3N2 circulating strains, particularly against the Yunnan 2021-2022 A/H3N2 strains. These findings aid our understanding of evolutionary characteristics and epidemiological monitoring of the A/H3N2 viruses and provide in-depth insights into the protective efficacy of influenza vaccines.

Evolution of Influenza A(H3N2) Viruses in 2 Consecutive Seasons of Genomic Surveillance, 2021-2023

Background

The circulation and the genomic evolution of influenza A(H3N2) viruses during the 2021/2022 and 2022/2023 seasons were studied and associated with infection outcomes.

Methods

Remnant influenza A-positive samples following standard-of-care testing from patients across the Johns Hopkins Health System (JHHS) were used for the study. Samples were randomly selected for whole viral genome sequencing. The sequence-based pEpitope model was used to estimate the predicted vaccine efficacy (pVE) for circulating H3N2 viruses. Clinical data were collected and associated with viral genomic data.

Results

A total of 121 683 respiratory specimens were tested for influenza at JHHS between 1 September 2021 and 31 December 2022. Among them, 6071 (4.99%) tested positive for influenza A. Of these, 805 samples were randomly selected for sequencing, with hemagglutinin (HA) segments characterized for 610 samples. Among the characterized samples, 581 were H3N2 (95.2%). Phylogenetic analysis of HA segments revealed the exclusive circulation of H3N2 viruses with HA segments of the 3C.2a1b.2a.2 clade. Analysis of a total of 445 complete H3N2 genomes revealed reassortments; 200 of 227 of the 2022/2023 season genomes (88.1%) were found to have reassorted with clade 3C.2a1b.1a. The pVE was estimated to be -42.53% for the 2021/2022 season and 30.27% for the 2022/2023 season. No differences in clinical presentations or admissions were observed between the 2 seasons.

Conclusions

The increased numbers of cases and genomic diversity of influenza A(H3N2) during the 2022/2023 season were not associated with a change in disease severity compared to the previous influenza season.

Investigating the Genetic Diversity of H5 Avian Influenza Viruses in the United Kingdom from 2020-2022

Since 2020, the United Kingdom and Europe have experienced annual epizootics of high-pathogenicity avian influenza virus (HPAIV). The first epizootic, during the autumn/winter of 2020-2021, involved six H5Nx subtypes, although H5N8 HPAIV dominated in the United Kingdom. While genetic assessments of the H5N8 HPAIVs within the United Kingdom demonstrated relative homogeneity, there was a background of other genotypes circulating at a lower degree with different neuraminidase and internal genes.  Following a small number of detections of H5N1 in wild birds over the summer of 2021, the autumn/winter of 2021-2022 saw another European H5 HPAIV epizootic that dwarfed the prior epizootic. This second epizootic was dominated almost exclusively by H5N1 HPAIV, although six distinct genotypes were defined. We have used genetic analysis to evaluate the emergence of different genotypes and proposed reassortment events that have been observed. The existing data suggest that the H5N1 viruses circulating in Europe during late 2020 continued to circulate in wild birds throughout 2021, with minimal adaptation, but then went on to reassort with AIVs in the wild bird population. We have undertaken an in-depth genetic assessment of H5 HPAIVs detected in the United Kingdom over two winter seasons and demonstrate the utility of in-depth genetic analyses in defining the diversity of H5 HPAIVs circulating in avian species, the potential for zoonotic risk, and whether incidents of lateral spread can be defined over independent incursions of infections from wild birds. This provides key supporting data for mitigation activities. IMPORTANCE High-pathogenicity avian influenza virus (HPAIV) outbreaks devastate avian species across all sectors, having both economic and ecological impacts through mortalities in poultry and wild birds, respectively. These viruses can also represent a significant zoonotic risk. Since 2020, the United Kingdom has experienced two successive outbreaks of H5 HPAIV. While H5N8 HPAIV was predominant during the 2020-2021 outbreak, other H5 subtypes were also detected. The following year, there was a shift in the subtype dominance to H5N1 HPAIV, but multiple H5N1 genotypes were detected. Through the thorough utilization of whole-genome sequencing, it was possible to track and characterize the genetic evolution of these H5 HPAIVs in United Kingdom poultry and wild birds. This enabled us to assess the risk posed by these viruses at the poultry-wild bird and the avian-human interfaces and to investigate the potential lateral spread between infected premises, a key factor in understanding the threat to the commercial sector.

Factors associated with viral RNA shedding and evaluation of potential viral infectivity at returning to school in influenza outpatients after treatment with baloxavir marboxil and neuraminidase inhibitors during 2013/2014-2019/2020 seasons in Japan: an observational study

BACKGROUND

This study assessed the differences in daily virus reduction and the residual infectivity after the recommended home stay period in Japan in patients infected with influenza and treated with baloxavir (BA), laninamivir (LA), oseltamivir (OS), and zanamivir (ZA).

METHODS

We conducted an observational study on children and adults at 13 outpatient clinics in 11 prefectures in Japan during seven influenza seasons from 2013/2014 to 2019/2020. Virus samples were collected twice from influenza rapid test-positive patients at the first and second visit 4-5 days after the start of treatment. The viral RNA shedding was quantified using quantitative RT-PCR. Neuraminidase (NA) and polymerase acidic (PA) variant viruses that reduce susceptibility to NA inhibitors and BA, respectively, were screened using RT-PCR and genetic sequencing. Daily estimated viral reduction was evaluated using univariate and multivariate analyses for the factors such as age, treatment, vaccination status, or the emergence of PA or NA variants. The potential infectivity of the viral RNA shedding at the second visit samples was determined using the Receiver Operator Curve based on the positivity of virus isolation.

RESULTS

Among 518 patients, 465 (80.0%) and 116 (20.0%) were infected with influenza A (189 with BA, 58 with LA, 181 with OS, 37 with ZA) and influenza B (39 with BA, 10 with LA, 52 with OS, 15 with ZA). The emergence of 21 PA variants in influenza A was detected after BA treatment, but NA variants were not detected after NAIs treatment. Multiple linear regression analysis showed that the daily viral RNA shedding reduction in patients was slower in the two NAIs (OS and LA) than in BA, influenza B infection, aged 0-5 years, or the emergence of PA variants. The residual viral RNA shedding potentially infectious was detected in approximately 10-30% of the patients aged 6-18 years after five days of onset.

CONCLUSIONS

Viral clearance differed by age, type of influenza, choice of treatment, and susceptibility to BA. Additionally, the recommended homestay period in Japan seemed insufficient, but reduced viral spread to some extent since most school-age patients became non-infectious after 5 days of onset.

Exploring the Correlation Between Influenza A Virus (H3N2) Infections and Neurological Manifestations: A Scoping Review

Influenza A virus (IAV), particularly the H3N2 variant, is known to cause respiratory manifestations, but it can also lead to neurological complications ranging from mild symptoms like headache and dizziness to severe conditions such as encephalitis and acute necrotizing encephalopathy (ANE). In this article, the correlation between the H3N2 variant of the IAV and neurological manifestations is discussed. Additionally, prompt recognition and treatment of influenza-associated neurological manifestations are highlighted to prevent infection-related long-term complications. This review briefly discusses various neurological complications linked to IAV infections, such as encephalitis, febrile convulsions, and acute disseminated encephalomyelitis, and the potential mechanisms involved in the development of neurological complications.

Whole-Genome Analysis of Influenza A(H3N2) and B/Victoria Viruses Detected in Myanmar during the COVID-19 Pandemic in 2021

An influenza circulation was observed in Myanmar between October and November in 2021. Patients with symptoms of influenza-like illness were screened using rapid diagnostic test (RDT) kits, and 147/414 (35.5%) upper respiratory tract specimens presented positive results. All RDT-positive samples were screened by a commercial multiplex real-time polymerase chain reaction (RT-PCR) assay, and 30 samples positive for influenza A(H3N2) or B underwent further typing/subtyping for cycle threshold (Ct) value determination based on cycling probe RT-PCR. The majority of subtyped samples (n = 13) were influenza A(H3N2), while only three were B/Victoria. Clinical samples with low Ct values obtained by RT-PCR were used for whole-genome sequencing via next-generation sequencing technology. All collected viruses were distinct from the Southern Hemisphere vaccine strains of the corresponding season but matched with vaccines of the following season. Influenza A(H3N2) strains from Myanmar belonged to clade 2a.3 and shared the highest genetic proximity with Bahraini strains. B/Victoria viruses belonged to clade V1A.3a.2 and were genetically similar to Bangladeshi strains. This study highlights the importance of performing influenza virus surveillance with genetic characterization of the influenza virus in Myanmar, to contribute to global influenza surveillance during the COVID-19 pandemic.

Successful Whole Genome Nanopore Sequencing of Swine Influenza A Virus (swIAV) Directly from Oral Fluids Collected in Polish Pig Herds

Influenza A virus (IAV) is a single-stranded, negative-sense RNA virus and a common cause of seasonal flu in humans. Its genome comprises eight RNA segments that facilitate reassortment, resulting in a great variety of IAV strains. To study these processes, the genetic code of each segment should be unraveled. Fortunately, new third-generation sequencing approaches allow for cost-efficient sequencing of IAV segments. Sequencing success depends on various factors, including proper sample storage and processing. Hence, this work focused on the effect of storage of oral fluids and swIAV sequencing. Oral fluids (n = 13) from 2017 were stored at -22 °C and later transferred to -80 °C. Other samples (n = 21) were immediately stored at -80 °C. A reverse transcription quantitative PCR (RT-qPCR) pre- and post-storage was conducted to assess IAV viral loads. Next, samples were subjected to two IAV long-read nanopore sequencing methods to evaluate success in this complex matrix. A significant storage-associated loss of swIAV loads was observed. Still, a total of 17 complete and 6 near-complete Polish swIAV genomes were obtained. Genotype T, (H1avN2, seven herds), P (H1N1pdm09, two herds), U (H1avN1, three herds), and A (H1avN1, 1 herd) were circulated on Polish farms. In conclusion, oral fluids can be used for long-read swIAV sequencing when considering appropriate storage and segment amplification protocols, which allows us to monitor swIAV in an animal-friendly and cost-efficient manner.