Genomic epidemiology and evolutionary analysis during XBB.1.16-predominant periods of SARS-CoV-2 omicron variant in Bangkok, Thailand: December 2022-August 2023.

The growing occurrence of novel recombinants, such as XBB.1.16, has emerged and become predominant, raising concerns about the impact of genomic recombination on the evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This study investigated the molecular epidemiological trends and evolution of the Omicron XBB.1.16 epidemic in Bangkok between December 2022 and August 2023. Partial spike and complete genome sequencing of SARS-CoV-2 samples collected from collaborating hospitals were performed. The analysis of 491 partial spike sequences identified 15 distinct lineages, with XBB.1.16 dominating the lineages beginning in March 2023. Phylogenetic analysis revealed at least four distinct XBB.1.16 lineages, suggesting multiple independent introductions into Bangkok. The estimated emergence of XBB.1.16 occurred approximately in January 2022, with an evolutionary rate of 0.79 × 10-3 substitutions per site per year. Monitoring the genomic epidemiology and evolution of XBB.1.16 is vital for the early detection of new strains or emerging variants, which may guide vaccine design and the inclusion of new vaccine strains.

In vitro modeling of isoniazid resistance mechanisms in Mycobacterium tuberculosis H37Rv.

Introduction: Mycobacterium tuberculosis (MTB), the causative agent of tuberculosis, has been a global threat to human beings for several decades. Treating tuberculosis has become more difficult as the prevalence of drug-resistant tuberculosis has increased globally. Evidence suggests that the comprehensive landscape of resistance mechanisms in MTB is ambiguous. More importantly, little is known regarding the series of events connected to resistance mechanisms in MTB before exposure to anti-TB drugs, during exposure to the drugs, and finally, when the MTB becomes resistant after exposure, upon analyses of its genome. Methods: We used the wild-type strain of MTB (H37Rv) in an in vitro model for generating induced resistance using a sub-inhibitory concentration of isoniazid, and the generated resistance-associated variants (RAVs) were identified using the whole genome sequencing method. Results: The detection of an inhA promoter mutation (fabG1-15C>T), which results in increased production of InhA protein, was found to be a major mechanism for developing resistance to isoniazid in the first place. We observed adaptation of MTB resistance mechanisms in high isoniazid stress by alteration and abolishment of KatG due to the detection of katG S315N, the common region of mutation that confers isoniazid resistance, along with katG K414N, katG N138S, and katG A162E. Furthermore, we detected the ahpC-72C>T and ahpC 21C>A mutations, but further investigation is needed to determine their role in compensating for the loss of KatG activity. Discussion: This suggests that increased InhA production is the main mechanism where there are low levels of isoniazid, whereas the alteration of KatG was found to be utilized in mycobacterium with a high concentration of isoniazid. Our work demonstrates that this in vitro approach of generating induced resistance could provide clinically relevant information after the fabG1-15C>T mutation, which is the common mutation found in clinical isolates. Moreover, other mutations detected in this work can also be found in clinical isolates. These findings may shed light on the impact of isoniazid in generating RAV and the resistance mechanism scenario that mycobacterium used under various isoniazid-pressuring conditions. More research is needed to understand better the role of RAV and mechanical resistance events within the mycobacterium genome in promoting a promising drug prediction platform that could lead to the right treatment for patients with MDR-TB and XDR-TB.

Investigation of the Molecular Epidemiology and Evolution of Circulating Severe Acute Respiratory Syndrome Coronavirus 2 in Thailand from 2020 to 2022 via Next-Generation Sequencing.

Coronavirus disease 2019 (COVID-19) is an infectious condition caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), which surfaced in Thailand in early 2020. The current study investigated the SARS-CoV-2 lineages circulating in Thailand and their evolutionary history. Complete genome sequencing of 210 SARS-CoV-2 samples collected from collaborating hospitals and the Institute of Urban Disease Control and Prevention over two years, from December 2020 to July 2022, was performed using next-generation sequencing technology. Multiple lineage introductions were observed before the emergence of the B.1.1.529 omicron variant, including B.1.36.16, B.1.351, B.1.1, B.1.1.7, B.1.524, AY.30, and B.1.617.2. The B.1.1.529 omicron variant was subsequently detected between January 2022 and June 2022. The evolutionary rate for the spike gene of SARS-CoV-2 was estimated to be between 0.87 and 1.71 × 10-3 substitutions per site per year. There was a substantial prevalence of the predominant mutations C25672T (L94F), C25961T (T190I), and G26167T (V259L) in the ORF3a gene during the Thailand outbreaks. Complete genome sequencing can enhance the prediction of future variant changes in viral genomes, which is crucial to ensuring that vaccine strains are protective against worldwide outbreaks.

Molecular characterisation and tracking of severe acute respiratory syndrome coronavirus 2 in Thailand, 2020-2022.

The global COVID-19 pandemic, caused by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was first detected in China in December 2019. To date, there have been approximately 3.4 million reported cases of COVID-19 and over 24,000 deaths in Thailand. In this study, we investigated the molecular characteristics and evolution of SARS-CoV-2 in Thailand from 2020 to 2022. Two hundred sixty-eight SARS-CoV-2 isolates, collected mostly in Bangkok from COVID-19 patients, were characterised by partial genome sequencing. Moreover, the viruses in 5,627 positive SARS-CoV-2 samples were identified as viral variants - B.1.1.7 (Alpha), B.1.617.2 (Delta), B.1.1.529 (Omicron/BA.1), or B.1.1.529 (Omicron/BA.2) - by multiplex real-time reverse transcription polymerase chain reaction (RT-PCR) assays. The results revealed that B.1.36.16 caused the predominant outbreak in the second wave (December 2020-January 2021), B.1.1.7 (Alpha) in the third wave (April-June 2021), B.1.617.2 (Delta) in the fourth wave (July-December 2021), and B.1.1.529 (Omicron) in the fifth wave (January-March 2022). The evolutionary rate of the viral genome was 2.60 × 10-3 (95% highest posterior density [HPD], 1.72 × 10-3 to 3.62 × 10-3) nucleotide substitutions per site per year. Continued molecular surveillance of SARS-CoV-2 is crucial for monitoring emerging variants with the potential to cause new COVID-19 outbreaks.

SARS-CoV-2 delta variant infection in domestic dogs and cats, Thailand.

In June-September 2021, we investigated severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infections in domestic dogs and cats (n = 225) in Bangkok and the vicinities, Thailand. SARS-CoV-2 was detected in a dog and a cat from COVID-19 positive households. Whole genome sequence analysis identified SARS-CoV-2 delta variant of concern (B.1.617.2). Phylogenetic analysis showed that SARS-CoV-2 isolated from dog and cat were grouped into sublineage AY.30 and AY.85, respectively. Antibodies against SARS-CoV-2 could be detected in both dog (day 9) and cat (day 14) after viral RNA detection. This study raises awareness on spill-over of variant of concern in domestic animals due to human-animal interface. Thus, surveillance of SARS-CoV-2 in domestic pets should be routinely conducted.

First cases of SARS-CoV-2 infection in dogs and cats in Thailand.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused the coronavirus disease 2019 (COVID-19) pandemic in humans since late 2019. Here, we investigated SARS-CoV-2 infection in dogs and cats during COVID-19 quarantine at private veterinary hospitals in Thailand. From April to May 2021, we detected SARS-CoV-2 in three out of 35 dogs and one out of nine cats from four out of 17 households with confirmed COVID-19 patients. SARS-CoV-2 RNA was detected from one of the nasal, oral, rectal and environmental swabs of dog-A (15 years old, mixed breed, male dog), cat-B (1 year old, domestic shorthair, male cat), dog-C (2 years old, mixed breed, female dog) and dog-D (4 years old, Pomeranian, female dog). The animals tested positive for SARS-CoV-2 RNA from 4 to 30 days after pet owners were confirmed to be COVID-19 positive. The animals consecutively tested positive for SARS-CoV-2 RNA for 4 to 10 days. One dog (dog-A) showed mild clinical signs, while the other dogs and a cat remained asymptomatic during quarantine at the hospitals. SARS-CoV-2 specific neutralizing antibodies were detected in both the dogs and cat by surrogate virus neutralization tests. Phylogenetic and genomic mutation analyses of whole genome sequences of three SARS-CoV-2 strains from the dogs and cat revealed SARS-CoV-2 of the Alpha variant (B.1.1.7 lineage). Our findings are suggestive of human-to-animal transmission of SARS-CoV-2 in COVID-19-positive households and contamination of viral RNA in the environment. Public awareness of SARS-CoV-2 infection in pet dogs and cats in close contact with COVID-19 patients should be raised.

Genetic diversity and evolution of enterovirus A71 subgenogroup C1 from children with hand, foot, and mouth disease in Thailand.

Enterovirus A71 (EV-A71) can cause hand, foot, and mouth disease (HFMD) in children and may be associated with severe neurological complications. There have been numerous reports of increased incidence of EV-A71 subgenogroup C1 (EV-A71 C1) infections associated with neurological diseases since the first occurrence in Germany in 2015. Here, we describe 11 full-length genome sequences of 2019 EV-A71 C1 strains isolated from HFMD patients in Thailand from 2019 to early 2020. The genetic evolution of 2019 EV-A71 C1 was traced in the outbreaks, and the emergence of multiple lineages was detected. Our results demonstrated that 2019 EV-A71 C1 from Thailand emerged through recombination between its nonstructural protein gene and those of other EV-A genotypes. Bayesian-based phylogenetic analysis showed that the 2019 EV-A71 C1 Thai strains share a common ancestor with variants in Europe (Denmark and France). The substitution rate for the 2019 EV-A71 C1 genome was estimated to be 4.38 × 10-3 substitutions/(site∙year-1) (95% highest posterior density interval: 3.84-4.94 × 10-3 substitutions/[site∙year-1]), approximating that observed between previous EV-A71 C1 outbreaks. These data are essential for understanding the evolution of EV-A C1 during the ongoing HFMD outbreak and may be relevant to disease outcomes in children worldwide.

Molecular characterization identifies intra-host recombination and zoonotic potential of canine rotavirus among dogs from Thailand.

From September 2016 to January 2019, we collected 710 rectal swabs from both healthy and sick dogs from small animal hospitals in 5 provinces of Thailand. The samples were tested for canine rotavirus group A (CRV) by using one-step RT-PCR specific to the VP6 gene. Our results showed that 0.70% (5/710) were positive for CRV. The five CRVs were then characterized by whole-genome sequencing. Our results showed that the genotype of Thai CRVs is G3P[3], which is the predominant genotype reported in dogs. The Thai CRVs posed a novel genetic constellation 'G3-P[3]-I3-R3-C3-M3-A9-N2-T3-E3-H6', which has never been reported in CRVs from dogs but has been reported in rotaviruses from humans. Based on phylogenetic analysis, the Thai CRVs are the result of multiple reassortments in which gene segments might have originated from human and bat rotaviruses and suggests the zoonotic potential of the virus.

Molecular epidemiology of the first wave of severe acute respiratory syndrome coronavirus 2 infection in Thailand in 2020.

The coronavirus disease 2019 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a major global concern. Several SARS-CoV-2 gene mutations have been reported. In the current study associations between SARS-CoV-2 gene variation and exposure history during the first wave of the outbreak in Thailand between January and May 2020 were investigated. Forty samples were collected at different time points during the outbreak, and parts of the SARS-CoV-2 genome sequence were used to assess genomic variation patterns. The phylogenetics of the 40 samples were clustered into L, GH, GR, O and T types. T types were predominant in Bangkok during the first local outbreak centered at a boxing stadium and entertainment venues in March 2020. Imported cases were infected with various types, including L, GH, GR and O. In southern Thailand introductions of different genotypes were identified at different times. No clinical parameters were significantly associated with differences in genotype. The results indicated local transmission (type T, Spike protein (A829T)) and imported cases (types L, GH, GR and O) during the first wave in Thailand. Genetic and epidemiological data may contribute to national policy formulation, transmission tracking and the implementation of measures to control viral spread.

Molecular modelling investigation for drugs and nutraceuticals against protease of SARS-CoV-2.

The widespread problem of a 2019-novel coronavirus (SARS-CoV-2) strain outbreak in Wuhan, China has prompted a search for new drugs to protect against and treat this disease. It is necessary to immediately investigate this due to the mutation of the viral genome and there being no current protective vaccines or therapeutic drugs. Molecular modelling and molecular docking based on in silico screening strategies were employed to determine the potential activities of seven HIV protease (HIV-PR) inhibitors, two flu drugs, and eight natural compounds. The computational approach was carried out to discover the structural modes with a high binding affinity for these drugs on the homology structure of the Wuhan coronavirus protease (SARS-CoV-2 PR). From the theoretical calculations, all the drugs and natural compounds demonstrated various favorable binding affinities. An interesting finding was that the natural compounds tested had a higher potential binding activity with the pocket sites of SARS-CoV-2 PR compared to the groups of HIV-PR inhibitors. The binding modes of each complex illustrated between the drugs and compounds interacted with the functional group of amino acids in the binding pocket via hydrophilic, hydrophobic, and hydrogen bond interactions using the molecular dynamics simulation technique. This result supports the idea that existing protease inhibitors and natural compounds could be used to treat the new coronavirus. This report sought to provide fundamental knowledge as preliminary experimental data to propose an existing nutraceutical material against viral infection. Collectively, it is suggested that molecular modelling and molecular docking are suitable tools to search and screen for new drugs and natural compounds that can be used as future treatments for viral diseases.

Human Norovirus Infection in Dogs, Thailand.

In July 2018, recombinant norovirus GII.Pe-GII.4 Sydney was detected in dogs who had diarrhea in a kennel and in children living on the same premises in Thailand. Whole-genome sequencing and phylogenetic analysis of 4 noroviruses from Thailand showed that the canine norovirus was closely related to human norovirus GII.Pe-GII.4 Sydney, suggesting human-to-canine transmission.

A Predominance of Clade 17 Candida albicans Isolated From Hemocultures in a Tertiary Care Hospital in Thailand.

Candida albicans is one of the most common human fungal pathogens. Candidemia has significant mortality globally. No epidemiological study of C. albicans based on multilocus sequence typing (MLST) has been conducted in Thailand. Therefore, MLST was used to study the molecular epidemiology of C. albicans blood strains in a large Thai teaching hospital. In vitro virulence phenotypes and antifungal susceptibility testing by broth microdilution were also conducted. Forty-six C. albicans blood strains from 37 patients were collected from the Department of Microbiology, Siriraj Hospital, in 2016 and 2017. Most patients (71.8%) were more than 60 years old, and the case fatality rate was 54.8%. The male-to-female ratio was 5:3. Thirty-four diploid sequence types (DSTs), including six new DSTs, were identified, with DST2514 (8.7%) and DST2876 (8.7%) as the most common DSTs. Strains were clustered into nine clades. Unlike other studies of C. albicans blood strains in Asia, clade 17 was the most common (13 strains, 28.3%). Sequential allelic changes were evident in sequential strains from one patient. All strains produced phospholipase and hemolysin, while none produced proteinase. The ability to form biofilm was found in 82.6% of the strains. Clade 17 strains showed significantly stronger hemolytic activity than non-clade 17 strains (69.2% versus 27.3%; p = 0.022). However, no significant association existed between clades and patient mortalities. All were susceptible or wild type to anidulafungin (MIC range = 0.015-0.12 and GM = 0.030), micafungin (MIC range = ≤ 0.008-0.015 and GM = 0.008), caspofungin (MIC range = 0.008-0.12 and GM = 0.036), and amphotericin B (MIC range = 0.25-0.5 and GM = 0.381). Only one strain was resistant to voriconazole (MIC range = ≤ 0.008 to ≥ 8 and GM = 0.010) and fluconazole (MIC range = 0.12-16 and GM = 0.398). In conclusion, a high prevalence of clade 17 C. albicans blood strains was found in Thailand, in contrast to other Asian countries. This unique finding might be explained by the strong hemolytic activity that is required for bloodstream infection of C. albicans.

Association of ubiA mutations and high-level of ethambutol resistance among Mycobacterium tuberculosis Thai clinical isolates.

Ethambutol (EMB) is the first-line antituberculosis drug and a potential supplementary agent for a treatment regimen of multidrug-resistant (MDR) and extensively drug-resistant (XDR) tuberculosis (TB). It has long been known that mutations in embCAB operon, encoding EMB target, arabinosyltransferase, confer resistance to EMB. Recently, ubiA was additionally reported to be specifically associated with high-level EMB resistance in Mycobacterium tuberculosis. However, such information on ubiA is very limited. This study aimed to investigate correlations between mutations in ubiA and phenotypic EMB resistance among EMB-resistant (EMBR) M. tuberculosis Thai clinical isolates. Minimum inhibitory concentration (MIC) level of EMB and ubiA sequences were determined and analyzed. Of 68 EMBR-MDR isolates, 8.9% harbored mutations in ubiA. However, 10.0% and 46.6% of EMB-sensitive (EMBS)-MDR and pan-susceptible isolates also had ubiA mutations detected, respectively. Most nonsynonymous mutations, L31P, A35S, and V55M were only found in the EMBR-MDR isolates except E149D which was also found in EMBS-MDR and pan-susceptible isolates. A further phylogenetic analysis based on spoligotyping and IS6110-RFLP illustrated that E149D was in fact associated to EAI-families rather than EMB resistance. By excluding synonymous mutations and the E149D, we found a high correlation between ubiA mutations and high-level of EMB resistance with 100.0% specificity. In conclusion, despite its rare occurrence, mutations in ubiA can potentially be a marker for a detection of high level of EMB resistance at least in the MDR M. tuberculosis background.

Species distribution, virulence factors, and antifungal susceptibility among Candida parapsilosis complex isolates from clinical specimens at Siriraj Hospital, Thailand, from 2011 to 2015.

The Candida parapsilosis complex has been described as the second or third most common yeast species isolated from patients with bloodstream infections worldwide. This complex consists of three species: C. parapsilosis sensu stricto, C. orthopsilosis, and C. metapsilosis. The distribution of species in this complex has never been studied in Thailand. Here we investigated the molecular epidemiology, in vitro on virulence factors, and antifungal susceptibility profiles of isolates of these three species collected from patients in Siriraj Hospital, Thailand, from 2011 to 2015. Of the 96 C. parapsilosis complex isolates analyzed, 66 (68.75%) were identified as C. parapsilosis sensu stricto, 28 (29.17%) as C. orthopsilosis, and two (2.08%) as C. metapsilosis. Most strains were isolated from blood (81.25%). Proteinase activity was only detected in four (6.06%) and two (7.14%) isolates of C. parapsilosis sensu stricto and C. orthopsilosis, respectively. Sixty (90.91%) isolates of C. parapsilosis sensu stricto, 12 (42.86%) isolates of C. orthopsilosis, and all C. metapsilosis isolates showed phospholipase activity. Psuedohyphae formation was only detected in 33 (50%) and 15 (53.57%) isolates of C. parapsilosis sensu stricto and C. orthopsilosis, respectively. All isolates were susceptible to caspofungin. Most (85-100%) isolates were susceptible to antifungal drugs, but 3.13 - 6.25% were resistant to voriconazole and fluconazole. In conclusion, our findings revealed that C. parapsilosis sensu stricto was the most common species among clinical isolates of the C. parapsilosis complex, and the most commonly used antifungal agents generally exhibited good in vitro activity against these strains.

Development of single-step multiplex real-time RT-PCR assays for rapid diagnosis of enterovirus 71, coxsackievirus A6, and A16 in patients with hand, foot, and mouth disease.

Real-time reverse-transcription polymerase chain reaction (rRT-PCR) to detect enterovirus 71 (EV-A71) and coxsackievirus A16 (CV-A16) has facilitated the rapid and accurate identification of the two most common etiological agents underlying hand, foot, and mouth disease (HFMD). However, the worldwide emergence of CV-A6 infection in HFMD necessitates development of an improved multiplex rRT-PCR method. To rapidly determine the etiology of HFMD, two rRT-PCR assays using TaqMan probes were developed to differentiate among three selected common enteroviruses (EV-A71, CV-A16 and CV-A6) and to enable broad detection of enteroviruses (pan-enterovirus assay). No cross-reactions were observed with other RNA viruses examined. The detection limits of both assays were 10 copies per microliter for EV-A71, CV-A6 and CV-A16, and pan-enterovirus. The methods showed high accuracy (EV-A71, 90.6%; CV-A6, 92.0%; CV-A16, 100%), sensitivity (EV-A71, 96.5%; CV-A6, 95.8%; CV-A16, 99.0%), and specificity (EV-A71, 100%; CV-A6, 99.9%; CV-A16, 99.9%) in testing clinical specimens (n=1049) during 2014-2016, superior to those of conventional RT-PCR. Overall, the multiplex rRT-PCR assays enabled highly sensitive detection and rapid simultaneous typing of EV-A71, CV-A6 and CV-A16, and enteroviruses, rendering them feasible and attractive methods for large-scale surveillance of enteroviruses associated with HFMD outbreaks.

Molecular epidemiology of cryptococcal genotype VNIc/ST5 in Siriraj Hospital, Thailand.

Despite the strong association between Cryptococcus neoformans infection and the Human immunodeficiency virus (HIV) status of patients globally, most cryptococcosis cases in Far East Asia occur in non-HIV individuals. Molecular epidemiological studies, using multilocus sequence typing (MLST), have shown that more than 95% of cryptococcal strains belong to a specific subtype of VNI. However, this association has never been specifically examined in other parts of Asia. Therefore, in this study, we investigated the VNIc/ST5 genotype distribution among cryptococcosis patients in Thailand. Fifty-one C. neoformans isolates were collected from clinical samples in Siriraj Hospital, Bangkok, Thailand. The strains were predominantly isolated from HIV-positive patients (88.57%) and all were molecular type VNI MATα. An MLST analysis identified five sequence types (ST) in Siriraj Hospital, of which ST4 (45.10%) and ST6 (35.29%) were most common, and ST5 (15.69%), ST32 (1.96%), and ST93 (1.96) were less common. Contrary to reports from Far East Asia, ST5 was predominantly (83.33%) found in HIV patients (P = 0.657), and there was no significant change in the prevalence of ST5 over the past 10 years (P = 0.548). A further analysis of comorbidities showed higher morbidity and delays in the cryptococcal diagnosis in patients with tuberculosis coinfection or without HIV. Our study suggests that although the Thai population is genetically closely related to the Far East Asian population, ST5 is not associated with non-HIV status in Thailand. Therefore, this association may not be related to the host's genetic background. However, its mechanism remains unclear.

Live Animal Markets in Minnesota: A Potential Source for Emergence of Novel Influenza A Viruses and Interspecies Transmission.

BACKGROUND: Live animal markets have been implicated in transmission of influenza A viruses (IAVs) from animals to people. We sought to characterize IAVs at 2 live animal markets in Minnesota to assess potential routes of occupational exposure and risk for interspecies transmission. METHODS: We implemented surveillance for IAVs among employees, swine, and environment (air and surfaces) during a 12-week period (October 2012-January 2013) at 2 markets epidemiologically associated with persons with swine-origin IAV (variant) infections. Real-time reverse transcription polymerase chain reaction (rRT-PCR), viral culture, and whole-genome sequencing were performed on respiratory and environmental specimens, and serology on sera from employees at beginning and end of surveillance. RESULTS: Nasal swabs from 11 of 17 (65%) employees tested positive for IAVs by rRT-PCR; 7 employees tested positive on multiple occasions and 1 employee reported influenza-like illness. Eleven of 15 (73%) employees had baseline hemagglutination inhibition antibody titers ≥40 to swine-origin IAVs, but only 1 demonstrated a 4-fold titer increase to both swine-origin and pandemic A/Mexico/4108/2009 IAVs. IAVs were isolated from swine (72/84), air (30/45), and pen railings (5/21). Whole-genome sequencing of 122 IAVs isolated from swine and environmental specimens revealed multiple strains and subtype codetections. Multiple gene segment exchanges among and within subtypes were observed, resulting in new genetic constellations and reassortant viruses. Genetic sequence similarities of 99%-100% among IAVs of 1 market customer and swine indicated interspecies transmission. CONCLUSIONS: At markets where swine and persons are in close contact, swine-origin IAVs are prevalent and potentially provide conditions for novel IAV emergence.

Molecular epidemiology and phylogenetic analyses of influenza B virus in Thailand during 2010 to 2014.

Influenza B virus remains a major contributor to the seasonal influenza outbreak and its prevalence has increased worldwide. We investigated the epidemiology and analyzed the full genome sequences of influenza B virus strains in Thailand between 2010 and 2014. Samples from the upper respiratory tract were collected from patients diagnosed with influenza like-illness. All samples were screened for influenza A/B viruses by one-step multiplex real-time RT-PCR. The whole genome of 53 influenza B isolates were amplified, sequenced, and analyzed. From 14,418 respiratory samples collected during 2010 to 2014, a total of 3,050 tested positive for influenza virus. Approximately 3.27% (471/14,418) were influenza B virus samples. Fifty three isolates of influenza B virus were randomly chosen for detailed whole genome analysis. Phylogenetic analysis of the HA gene showed clusters in Victoria clades 1A, 1B, 3, 5 and Yamagata clades 2 and 3. Both B/Victoria and B/Yamagata lineages were found to co-circulate during this time. The NA sequences of all isolates belonged to lineage II and consisted of viruses from both HA Victoria and Yamagata lineages, reflecting possible reassortment of the HA and NA genes. No significant changes were seen in the NA protein. The phylogenetic trees generated through the analysis of the PB1 and PB2 genes closely resembled that of the HA gene, while trees generated from the analysis of the PA, NP, and M genes showed similar topology. The NS gene exhibited the pattern of genetic reassortment distinct from those of the PA, NP or M genes. Thus, antigenic drift and genetic reassortment among the influenza B virus strains were observed in the isolates examined. Our findings indicate that the co-circulation of two distinct lineages of influenza B viruses and the limitation of cross-protection of the current vaccine formulation provide support for quadrivalent influenza vaccine in this region.

Prevalence and genetic characterization of human coronaviruses in southern Thailand from July 2009 to January 2011.

Emergence of viruses belonging to the coronavirus family has been widespread in the past, causing respiratory infections in humans, such as severe acute respiratory syndrome (SARS). This study investigated the prevalence of human coronavirus (HCoV) and characterized the molecular viral genetics. We collected 1,254 samples from patients diagnosed with respiratory infection in southern Thailand from July 2009 to January 2011 and screened for HCoV by RT-PCR and genotyped by BLAST analysis of nsp12 gene. Phylogenetic analysis was performed based on S gene sequences. Thirty-five of 1,254 samples were positive for HCoV. Viral genotyping revealed 4 genotypes with HCoV-OC43 being the predominant genotype. Viral prevalence and genotype distribution were not in accordance with seasonal distribution. Phylogenetic analysis and deduced amino acid sequences of the S gene showed amino acid variations in each genotype. The S gene sequence of HCoV-OC43 genotype indicated that it resulted from recombination between subgenotypes B and C. Viral genetics analysis disclosed genetic variations of HCoV and additionally, it can provide information suitable for monitoring and prevention of the emergence and re-emergence of various types of coronavirus.