Detection of parechovirus-A in hospitalized children with acute lower respiratory infection in Myanmar, 2017-2018.

Parechovirus-A (PeV-A) causes emerging infection in children, and clinical presentation depends on genotype. The virus has been investigated mainly in developed countries; however, data from developing countries, especially in Asia, are sparse. This study investigated whether PeV-A circulated in children in Myanmar. This retrospective study evaluated PeV-A in nasopharyngeal samples from children aged 1 month to 12 years who were hospitalized with acute lower respiratory infection at Yankin Children Hospital, Yangon, Myanmar, during the period from May 2017 to April 2019. Real-time polymerase chain reaction (PCR) was used to detect PeV-A, and PCR-positive samples were used for genotyping and phylogenetic analysis. In total, 11/570 (1.9%) of samples were positive for PeV-A; 7 were successfully genotyped by sequencing the VP3/VP1 region, as follows: PeV-A1 (n = 4), PeV-A5 (n = 1), PeV-A6 (n = 1), and PeV-A14 (n = 1). Median age was 10.0 months (interquartile range 4.0-12.0 months), and other respiratory viruses were detected in all cases. Phylogenetic analysis showed that all detected PeV-A1 strains were in clade 1 A, which was a minor clade worldwide. Four PeV-A genotypes were detected in Myanmar. The clinical impact of PeV-A in children should be evaluated in future studies.

Clinical manifestations and outcome of viral acute lower respiratory infection in hospitalised children in Myanmar.

BACKGROUND: Acute lower respiratory infection (ALRI) remains the leading cause of death in children worldwide, and viruses have been the major cause of ALRI. In Myanmar, ALRI is associated with high morbidity and mortality in children, and detailed information on ALRI is currently lacking. METHODS: This prospective study investigated the viral aetiologies, clinical manifestations, and outcomes of ALRI in hospitalised children aged 1 month to 12 years at the Yankin Children Hospital, Yangon, Myanmar from May 2017 to April 2019. The sample size was set to 300 patients for each year. Two nasopharyngeal swabs were obtained for the patients with suspected viral ALRI; one for rapid tests for influenza and respiratory syncytial virus (RSV), and the other for real-time PCR for the 16 ALRI-causing viruses. Pneumococcal colonization rates were also investigated using real-time PCR. Clinical information was extracted from the medical records, and enrolled patients were categorised by age and severity for comparison. RESULTS: Among the 5463 patients admitted with a diagnosis of ALRI, 570 (10.4%) were enrolled in this study. The median age of the patients was 8 months (interquartile range, 4-15 months). The most common symptoms were cough (93%) and difficulty in breathing (73%), while the most common signs of ALRI were tachypnoea (78%) and chest indrawing (67%). A total of 16 viruses were detected in 502 of 570 patients' samples (88%), with RSV B (36%) and rhinovirus (28%) being the most commonly detected. Multiple viruses were detected in 221 of 570 samples (37%) collected from 570 patients. Severe ALRI was diagnosed in 107 of 570 patients (19%), and RSV B and human rhinovirus were commonly detected. The mortality rate was 5%; influenza virus A (29%) and RSV B (21%) were commonly detected, and stunting and lack of immunization were frequently observed in such cases. Additionally, 45% (259/570) of the patients had pneumococcal colonization. CONCLUSIONS: Viral ALRI in hospitalised children with a median of 8 months has significant morbidity and mortality rates in Myanmar. RSV and rhinovirus were the most commonly detected from nasopharyngeal swabs, while influenza virus and RSV were the most frequently associated with fatal cases.

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.

Detection of enterovirus D68 among children with severe acute respiratory infection in Myanmar.

BACKGROUND: Enterovirus D68 (EV-D68) is an important reemerging pathogen that causes severe acute respiratory infection and acute flaccid paralysis, mainly in children. Since 2014, EV-D68 outbreaks have been reported in the United States, Europe, and east Asia; however, no outbreaks have been reported in southeast Asian countries, including Myanmar, during the previous 10 years. METHODS: EV-D68 was detected in nasopharyngeal swabs from children with acute lower respiratory infections in Myanmar. The samples were previously collected from children aged 1 month to 12 years who had been admitted to the Yankin Children Hospital in Yangon, Myanmar, between May 2017 and January 2019. EV-D68 was detected with a newly developed EV-D68-specific real-time PCR assay. The clade was identified by using a phylogenetic tree created with the Bayesian Markov chain Monte Carlo method. RESULTS: During the study period, nasopharyngeal samples were collected from 570 patients. EV-D68 was detected in 42 samples (7.4 %)-11 samples from 2017 to 31 samples from 2018. The phylogenetic tree revealed that all strains belonged to clade B3, which has been the dominant clade worldwide since 2014. We estimate that ancestors of currently circulating genotypes emerged during the period 1980-2004. CONCLUSIONS: To our knowledge, this is the first report of EV-D68 detection in children with acute lower respiratory infections in Yangon, Myanmar, in 2017-2018. Detection and detailed virologic analyses of EV-D68 in southeast Asia is an important aspect of worldwide surveillance and will likely be useful in better understanding the worldwide epidemiologic profile of EV-D68 infection.

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.

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

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-3 substitutions/site/year) and neuraminidase (2.89 × 10-3 substitutions/site/year). The lowest rate was for non-structural protein segments (4.19 × 10-5 substitutions/site/year). The substantial genetic diversity that was revealed improved phylogenetic classification. This information will be particularly relevant for improving vaccine strain selection.

Epidemiology and Genetic Analysis of SARS-CoV-2 in Myanmar during the Community Outbreaks in 2020.

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 (Rt) 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.

Evolutionary analysis of human respiratory syncytial virus collected in Myanmar between 2015 and 2018.

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.

Rare influenza A (H3N2) variants with reduced sensitivity to antiviral drugs.

In 2007 and 2008 in Myanmar, we detected influenza viruses A (H3N2) that exhibited reduced sensitivity to both zanamivir and amantadine. These rare and naturally occurring viruses harbored a novel Q136K mutation in neuraminidase and S31N mutation in M2.

Epidemic of influenza A(H1N1)pdm09 analyzed by full genome sequences and the first case of oseltamivir-resistant strain in Myanmar 2017.

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 IC50 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.

Epidemiology of human influenza A and B viruses in Myanmar from 2005 to 2007.

OBJECTIVES: To perform genetic analysis of influenza A and B viruses in Myanmar from 2005 to 2007 and to determine the prevalence of amantadine-resistant influenza A viruses. METHODS: Phylogenies of the HA and NA genes were analyzed and mutations in M2 that confer resistance to amantadine were screened. RESULTS: Influenza in Myanmar exhibited seasonality, which coincided during the rainy season from June to August. Out of 2,618 samples, 76 influenza A and 132 influenza B viruses were isolated. Phylogenetic analysis showed that in 2005, 11 A/H1N1 isolates formed one cluster with A/Solomon Islands/3/2006 and were amantadine-sensitive strains. One A/H3N2 isolate was amantadine-resistant harboring S31N mutation in M2 and possessing S193F and D225N substitutions in HA (clade N), similar to A/Wisconsin/67/2005. No viruses were isolated in 2006 due to sample storage failure. In 2007, all 64 A/H3N2 isolates were amantadine-resistant and similar to A/Brisbane/10/2007. For influenza B, 3 Yamagata-lineage and 17 Victoria-lineage isolates were detected in 2005 and 112 Victoria-lineage viruses were isolated in 2007. All Victoria-lineage isolates were reassortants possessing NA derived from the Yamagata lineage. CONCLUSION: Continuous surveillance in tropical countries is important for elucidating the seasonality of influenza and determining the molecular characteristics of circulating strains.

Phylogeographic analysis of human influenza A and B viruses in Myanmar, 2010-2015.

We investigated the circulation patterns of human influenza A and B viruses in Myanmar between 2010 and 2015 by analyzing full HA genes. Upper respiratory tract specimens were collected from patients with symptoms of influenza-like illness. A total of 2,860 respiratory samples were screened by influenza rapid diagnostic test, of which 1,577 (55.1%) and 810 (28.3%) were positive for influenza A and B, respectively. Of the 1,010 specimens that were positive for virus isolation, 370 (36.6%) were A(H1N1)pdm09, 327 (32.4%) were A(H3N2), 130 (12.9%) B(Victoria), and 183 (18.1%) were B(Yamagata) viruses. Our data showed that influenza epidemics mainly occurred during the rainy season in Myanmar. Our three study sites, Yangon, Pyinmana, and Pyin Oo Lwin had similar seasonality and circulating type and subtype of influenza in a given year. Moreover, viruses circulating in Myanmar during the study period were closely related genetically to those detected in Thailand, India, and China. Phylogeographic analysis showed that A(H1N1)pdm09 viruses in Myanmar originated from Europe and migrated to other countries via Japan. Similarly, A(H3N2) viruses in Myanmar originated from Europe, and disseminated to the various countries via Australia. In addition, Myanmar plays a key role in reseeding of influenza B viruses to Southeast Asia and East Asia as well as Europe and Africa. Thus, we concluded that influenza virus in Myanmar has a strong link to neighboring Asian countries, Europe and Oceania.

Neuraminidase inhibitor susceptibility and evolutionary analysis of human influenza B isolates from three Asian countries during 2012-2015.

Influenza B viruses of both the Yamagata and the Victoria lineages are implicated in a large proportion of the morbidity and mortality associated with influenza outbreaks. In this study, we characterized the full genomes of 53 influenza B viruses isolated during 2012-2015 in three Asian countries: Japan, Myanmar, and Vietnam. Analysis of the hemagglutinin (HA) genes revealed co-circulation of both the Yamagata and Victoria lineages within the same season in these countries. Our analysis revealed, that a large proportion of viruses circulating during 2013-2014 in Japan and Vietnam were mismatched to the vaccine supporting the rationale for using quadrivalent vaccines. Molecular analysis of the neuraminidase (NA) genes did not reveal any of the previously reported substitutions associated with reduced susceptibility to neuraminidase inhibitors (NAIs). However, one isolate from Nagasaki displayed reduced inhibition by NAIs, associated with an NA-M426I substitution (N2-numbering). Phylogenetic analysis of the eight genome segments identified a 6 + 2 reassortant strain belonging to the Victoria lineage that circulated in Japan during the 2013-2014 season. This strain appears to have evolved from a descendent of a B/Brisbane/60/2008-like strain in an intra-lineage reassortment event involving the nucleoprotein (NP) and nonstructural (NS) genes. Therefore, influenza B strains circulating worldwide continue to evolve via complex reassortment events, which contribute to their survival and the emergence of new strains. These findings highlight the need for ongoing genome-wide studies of circulating viruses and assessing the implications of these evolutionary events on the vaccines.

Full Genome Characterization of Human Influenza A/H3N2 Isolates from Asian Countries Reveals a Rare Amantadine Resistance-Conferring Mutation and Novel PB1-F2 Polymorphisms.

Influenza A viruses evolve at a high rate requiring continuous monitoring to maintain the efficacy of vaccines and antiviral drugs. We performed next generation sequencing analysis of 100 influenza A/H3N2 isolates collected in four Asian countries (Japan, Lebanon, Myanmar, and Vietnam) during 2012-2015. Phylogenetic analysis revealed several reassortment events leading to the circulation of multiple clades within the same season. This was particularly evident during the 2013 and 2013/2014 seasons. Importantly, our data showed that certain lineages appeared to be fitter and were able to persist into the following season. The majority of A/H3N2 viruses continued to harbor the M2-S31N mutation conferring amantadine-resistance. In addition, an S31D mutation in the M2-protein, conferring a similar level of resistance as the S31N mutation, was detected in three isolates obtained in Japan during the 2014/2015 season. None of the isolates possessed the NA-H274Y mutation conferring oseltamivir-resistance, though a few isolates were found to contain mutations at the catalytic residue 151 (D151A/G/N or V) of the NA protein. These variations did not alter the susceptibility to neuraminidase inhibitors and were not detected in the original clinical specimens, suggesting that they had been acquired during their passage in MDCK cells. Novel polymorphisms were detected in the PB1-F2 open-reading frame resulting in truncations in the protein of 24-34 aminoacids in length. Thus, this study has demonstrated the utility of monitoring the full genome of influenza viruses to allow the detection of the potentially fittest lineages. This enhances our ability to predict the strain(s) most likely to persist into the following seasons and predict the potential degree of vaccine match or mismatch with the seasonal influenza season for that year. This will enable the public health and clinical teams to prepare for any related healthcare burden, depending on whether the vaccine match is predicted to be good or poor for that season.

Delayed emergence of oseltamivir-resistant seasonal influenza A (H1N1) and pandemic influenza A(H1N1)pdm09 viruses in Myanmar.

The prevalence and timing of emergence of oseltamivir-resistant seasonal and pandemic influenza A (H1N1) viruses in Myanmar in 2008 and 2009 are described in this report. In 2008, the oseltamivir-resistant seasonal H1N1 virus was detected at a lower rate (6%) and emerged at least 2 months later when compared with neighboring countries. Similarly, the prevalence of pandemic H1N1 virus was low (3%) and the timing of emergence was late (August 2009) in Myanmar. Interestingly, we detected three isolates that were resistant to both amantadine and oseltamivir. Limited movement of people into the country is attributed to the delayed emergence of drug-resistant seasonal and pandemic A(H1N1) viruses.