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.

Residual malaria among migrant workers in Myanmar: why still persistent and how to eliminate it?

BACKGROUND: Residual malaria is probably an important source for the re-emergence of malaria infection in the elimination era. Assessment to identify the factors influencing residual malaria in high-risk groups is needed to develop evidence-based decisions by stakeholders and policymakers. METHODS: This study was conducted to explore the factors influencing the residual malaria infection among migrant workers in two sentinel sites (endemic vs. pre-elimination areas) in Myanmar using the mixed-model method. RESULTS: A total of 102 migrant respondents (65 in Bamauk and 37 in Shwegyin) were included for the quantitative assessment using pretested questionnaires during household visits. Although 87.3% of them had insecticidal bed nets (ITNs/LLINs), only 68.3% of the migrants in Bamauk and 57.9% in Shwegyin used it regularly. The use of any bed net was high (79.9% in Bamauk vs. 91.0% in Shwegyin). The mean LLINs in their families were 1.64 (95%CI: 1.48-1.81) in Bamauk and 2.89 (95%CI: 2.67-3.11) in Shwegyin. Most of them received no health information for malaria prevention within the last year and their knowledge about malaria was low. Their working nature was a challenge for control measures against malaria in migrants. CONCLUSION: The strategy for distributing LLINs and health promotion activities for mobile/migrant populations should be reviewed, and an appropriate action plan should be developed for the specific migrant group. Moreover, health promotion activities for behavior change communication should be strengthened in the migrant population in Myanmar.

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.

Biochemical and biological characterization of the venoms of Naja kaouthia and Naja mandalayensis from Myanmar and neutralization effects of BPI cobra antivenom.

Snakebite is a neglected public health issue, with many scientific and medical issues to be solved. Cobras are among the most common venomous snakes in Myanmar and are responsible for a considerable number of severe snakebite envenoming. There are three species of cobra (Naja kaouthia, Naja mandalayensis and Ophiophagus hannah) in Myanmar. The study aims to characterize the N. kaouthia and N. mandalayensis venoms and to investigate the efficacy of anti-cobra antivenom (BPI) against the two venoms. Protein components and fibrinogenolytic activity were determined by SDS-PAGE. Enzymatic activities for PLA2, protease and acetylcholinesterase were determined by spectrophotometric method. Anticoagulant activity was determined by recalcification time of citrated human plasma. Myotoxicity, necrotizing activity, median lethal dose (LD50) and median effective dose (ED50) were determined by WHO recommended methods. The SDS-PAGE displayed the proteins and enzymes containing in two venoms were different. N. kaouthia venom exhibited more in PLA2, acetylcholinesterase, anticoagulant, fibrinogenolytic and necrotizing activities than N. mandalayensis venom. N. mandalayensis venom had more protease activity and myotoxicity than N. kaouthia venom. The median lethal dose (LD50) of N. kaouthia and N. mandalayensis venom was 4.33 µg/mouse and 5.04 µg/mouse respectively. Both venoms induced fibrinogen Aα chain degradation in 30 min (N. kaouthia) and in 6 h (N. mandalayensis). The same median effective dose (ED50) (19.56 µg/mouse) showed that anti-NK antivenom can neutralize against lethal effect of N. mandalayensis venom. It can also neutralize the protease activity, anticoagulant activity and fibrinogenolytic activity of both venoms. Immunodiffusion and immunoblotting studies showed that the antivenom recognized its homologous venom (N. kaouthia) and cross-reacted against the heterologous venom (N. mandalayensis). The anti-NK antivenom is suitable to use for N. mandalayensis bite if monospecific antivenom is not available.

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.

Surge of severe acute respiratory syndrome coronavirus 2 infections linked to single introduction of a virus strain in Myanmar, 2020.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is a major health concern globally. Genomic epidemiology is an important tool to assess the pandemic of coronavirus disease 2019 (COVID-19). Several mutations have been reported by genome analysis of the SARS-CoV-2. In the present study, we investigated the mutational and phylogenetic analysis of 30 whole-genome sequences for the virus's genomic characteristics in the specimens collected in the early phase of the pandemic (March-June, 2020) and the sudden surge of local transmission (August-September, 2020). The four samples in the early phase of infection were B.6 lineage and located within a clade of the samples collected at the same time in Singapore and Malaysia, while five returnees by rescue flights showed the lineage B. 1.36.1 (three from India), B.1.1 (one from India) and B.1.80 (one from China). However, there was no evidence of local spread from these returnees. Further, all 19 whole-genome sequences collected in the sudden surge of local transmission showed lineage B.1.36. The surge of the second wave on SARS-CoV-2 infection was linked to the single-introduction of a variant (B.1.36) that may result from the strict restriction of international travel and containment efforts. These genomic data provides the useful information to disease control and prevention strategy.

Pyrazinamide resistance and pncA mutations in drug resistant Mycobacterium tuberculosis clinical isolates from Myanmar.

Pyrazinamide (PZA) is an important anti-tuberculosis drug, which is active against semi-dormant bacilli and used as a component of first-line drugs and drug-resistant tuberculosis regimens. Mutations in pncA and its promoter region are main cause of PZA resistance. There are limited PZA susceptibility data as there is no routine drug susceptibility testing (DST) for PZA. This study was aimed to determine the proportion of PZA resistance among rifampicin-resistant tuberculosis patients and to identify mutations which are responsible for PZA resistance in pncA and its promoter region. Liquid-based DST was performed to detect PZA susceptibility on 192 culture positive rifampicin-resistant isolates collected from National Tuberculosis Reference Laboratory. Sequencing on pncA including its promoter region was performed and analysis was done on 157 isolates. Phenotypic PZA resistance was detected in 58.9% of isolates. Sixty-five different mutations were distributed in pncA or promoter region of 82 isolates. Sensitivity and specificity of pncA sequencing in detection of PZA resistance showed 89.8% and 95.6% respectively. High proportion of PZA resistance among rifampicin-resistant cases highlighted the need for effective treatment regimen development for PZA-resistant MDR-TB. It is also suggested that routine PZA susceptibility test should be incorporated to treatment monitoring regimen and National Drug Resistance surveys.

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.

Genotypes and genetic characters of Mycobacterium tuberculosis from Myanmar using three typing methods.

Knowledge on basic characteristics of Mycobacterium tuberculosis (MTB) is helpful to understand the disease epidemiology and support the prediction of clinical outcome of the disease. The aim of this study was to detect the genotypes and genotypic characters of clinical Mycobacterium tuberculosis (MTB) isolates from new and retreatment rifampicin-resistant patients using three different genotyping methods. Mycobacterial interspersed repetitive units-variable number tandem repeat (MIRU-VNTR) typing was used to determine the diversity of 222 clinical isolates. Spoligotyping and IS6110-restriction fragment length polymorphism (RFLP) typing were also used to investigate the genetic characters of 105 MTB strains. Among the 15 genotypes detected by MIRU-VNTR, Beijing strains were the most prevalent of all strains (54.8%); new cases (40.5%) and retreatment cases (69.4%), followed by EAI strain. Spoligotyping categorized the strains into 11 lineages and 13 orphans whereas 96 different IS6110 patterns were identified using RFLP method. The mode number of IS6110 was 18 and 20. Higher band numbers were found in Beijing genotype (p < 0.001). Clustering rates by spoligotyping, MIRU-VNTR and IS6110-RFLP typing were 0.714, 0.004 and 0.085, respectively. Discriminatory powers of spoligotyping, MIRU-VNTR typing and IS6110-RFLP typing were 0.637, 1.000 and 0.997, respectively. Dominant Beijing genotype in both new and retreatment cases denoting that prevailing tuberculosis in Myanmar changed from EAI to Beijing lineage.