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.

This study aimed to analyze the genetic and evolutionary characteristics of the influenza A/H3N2 viruses circulating in Myanmar from 2015 to 2019. Whole genomes from 79 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. Based on the WHO clades classification, the A/H3N2 strains in Myanmar from 2015 to 2019 collectively belonged to clade 3c.2. These strains were further defined based on hemagglutinin substitutions as follows: clade 3C.2a (n = 39), 3C.2a1 (n = 2), and 3C.2a1b (n = 38). Genetic analysis revealed that the Myanmar strains differed from the Southern Hemisphere vaccine strains each year, indicating that the vaccine strains did not match the circulating strains. The highest rates of nucleotide substitution were estimated for hemagglutinin (3.37 × 10⁻³ substitutions/site/year) and neuraminidase (2.89 × 10⁻³ substitutions/site/year). The lowest rate was for non-structural protein segments (4.19 × 10⁻⁵ substitutions/site/year). The substantial genetic diversity that was revealed improved phylogenetic classification. This information will be particularly relevant for improving vaccine strain selection.

A community outbreak of human influenza A(H1N1)pdm09 virus strains was observed in Myanmar in 2017. We investigated the circulation patterns, antigenicity, and drug resistance of 2017 influenza A(H1N1)pdm09 viruses from Myanmar and characterized the full genome of influenza virus strains in Myanmar from in-patients and out-patients to assess the pathogenicity of the viruses. Nasopharyngeal swabs were collected from out-patients and in-patients with acute respiratory tract infections in Yangon and Pyinmana City in Myanmar during January-December 2017. A total of 215 out-patients and 18 in-patients infected with A(H1N1)pdm09 were detected by virus isolation and real-time RT-PCR. Among the positive patients, 90.6% were less than 14 years old. Hemagglutination inhibition (HI) antibody titers against A(H1N1)pdm09 viruses in Myanmar were similar to the recommended Japanese influenza vaccine strain for 2017–2018 seasons (A/Singapore/GP1908/2015) and WHO recommended 2017 southern hemisphere vaccine component (A/Michigan/45/2015). Phylogenetic analysis of the hemagglutinin sequence showed that the Myanmar strains belonged to the genetic subclade 6B.1, possessing mutations of S162N and S164T at potential antigenic sites. However, the amino acid mutation at position 222, which may enhance the severity of disease and mortality, was not found. One case with no prior history of oseltamivir treatment possessed H275Y mutated virus in neuraminidase (NA), which confers resistance to oseltamivir and peramivir with elevated IC₅₀ values. The full genome sequence of Myanmar strains showed no difference between samples from in-patients and out-patients, suggesting no additional viral mutations associated with patient severity. Several amino acid changes were observed in PB2, PB1, and M2 of Myanmar strains when compared to the vaccine strain and other Asian strains. However, no mutations associated with pathogenicity were found in the Myanmar strains, suggesting that viral factors cannot explain the underlying reasons of the massive outbreak in Myanmar. This study reported the first detection of an oseltamivir-resistant influenza virus in Myanmar, highlighting the importance of continuous antiviral monitoring and genetic characterization of the influenza virus in Myanmar.

We studied genetic variation in the second hypervariable region (HVR) of the G gene of human respiratory syncytial virus (HRSV) from 1701 nasal swab samples collected from outpatients with acute respiratory infections at two general hospitals in the cities Yangon and Pyinmana in Myanmar from 2015 to 2018. HRSV genotypes were characterized using phylogenetic trees constructed using the maximum likelihood method. Time-scale phylogenetic tree analyses were performed using the Bayesian Markov chain Monte Carlo method. In total, 244 (14.3%) samples were HRSV-positive and were classified as HRSV-A (n = 84, 34.4%), HRSV-B (n = 158, 64.8%), and co-detection of HRSV-A/HRSV-B (n = 2, 0.8%). HRSV epidemics occurred seasonally between July (1.9%, 15/785) and August (10.5%, 108/1028), with peak infections in September (35.8%, 149/416) and October (58.2%, 89/153). HRSV infection rate was higher in children ≥1 year of age than in those <1 year of age (70.5% vs. 29.5%). The most common HRSV symptoms in children were cough (80%-90%) and rhinorrhea (70%-100%). The predominant genotypes were ON1for HRSV-A (78%) and BA9 for HRSV-B (64%). Time to the most recent common ancestor was 2014 (95% highest posterior density [HPD], 2012-2015) for HRSV-A ON1 and 2009 (95% HPD, 2004-2012) for HRSV-B BA9. The mean evolutionary rate (substitutions/site/year) for HRSV-B (2.12 × 10^-2, 95% HPD, 8.53 × 10^-3-3.63 × 10^-2) was slightly higher than that for HRSV-A (1.39 × 10^-2, 95% HPD, 6.03 × 10^-3-2.12 × 10^-2). The estimated effective population size (diversity) for HRSV-A increased from 2015 to 2016 and declined in mid-2018, whereas HRSV-B diversity was constant in 2015 and 2016 and increased in mid-2017. In conclusion, the dominant HRSV-A and HRSV-B genotypes in Myanmar were ON1 and BA9, respectively, between 2015 and 2018. HRSV-B evolved slightly faster than HRSV-A and exhibited unique phylogenetic characteristics.

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.

We aimed to analyze the situation of the first two epidemic waves in Myanmar using the publicly available daily situation of COVID-19 and whole-genome sequencing data of SARS-CoV-2. From March 23 to December 31, 2020, there were 33,917 confirmed cases and 741 deaths in Myanmar (case fatality rate of 2.18%). The first wave in Myanmar from March to July was linked to overseas travel, and then a second wave started from Rakhine State, a western border state, leading to the second wave spreading countrywide in Myanmar from August to December 2020. The estimated effective reproductive number (R_t) nationwide reached 6–8 at the beginning of each wave and gradually decreased as the epidemic spread to the community. The whole-genome analysis of 10 Myanmar SARS-CoV-2 strains together with 31 previously registered strains showed that the first wave was caused by GISAID clade O or PANGOLIN lineage B.6 and the second wave was changed to clade GH or lineage B.1.36.16 with a close genetic relationship with other South Asian strains. Constant monitoring of epidemiological situations combined with SARS-CoV-2 genome analysis is important for adjusting public health measures to mitigate the community transmissions of COVID-19.