Development and validation of multiplex real-time PCR for simultaneous detection of six bacterial pathogens causing lower respiratory tract infections and antimicrobial resistance genes.

BACKGROUND: Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, Escherichia coli, Streptococcus pneumoniae and Staphylococcus aureus are major bacterial causes of lower respiratory tract infections (LRTIs) globally, leading to substantial morbidity and mortality. The rapid increase of antimicrobial resistance (AMR) in these pathogens poses significant challenges for their effective antibiotic therapy. In low-resourced settings, patients with LRTIs are prescribed antibiotics empirically while awaiting several days for culture results. Rapid pathogen and AMR gene detection could prompt optimal antibiotic use and improve outcomes. METHODS: Here, we developed multiplex quantitative real-time PCR using EvaGreen dye and melting curve analysis to rapidly identify six major pathogens and fourteen AMR genes directly from respiratory samples. The reproducibility, linearity, limit of detection (LOD) of real-time PCR assays for pathogen detection were evaluated using DNA control mixes and spiked tracheal aspirate. The performance of RT-PCR assays was subsequently compared with the gold standard, conventional culture on 50 tracheal aspirate and sputum specimens of ICU patients. RESULTS: The sensitivity of RT-PCR assays was 100% for K. pneumoniae, A. baumannii, P. aeruginosa, E. coli and 63.6% for S. aureus and the specificity ranged from 87.5% to 97.6%. The kappa correlation values of all pathogens between the two methods varied from 0.63 to 0.95. The limit of detection of target bacteria was 1600 CFU/ml. The quantitative results from the PCR assays demonstrated 100% concordance with quantitative culture of tracheal aspirates. Compared to culture, PCR assays exhibited higher sensitivity in detecting mixed infections and S. pneumoniae. There was a high level of concordance between the detection of AMR gene and AMR phenotype in single infections. CONCLUSIONS: Our multiplex quantitative RT-PCR assays are fast and simple, but sensitive and specific in detecting six bacterial pathogens of LRTIs and their antimicrobial resistance genes and should be further evaluated for clinical utility.

Gallbladder carriage generates genetic variation and genome degradation in Salmonella Typhi.

Despite recent advances in typhoid fever control, asymptomatic carriage of Salmonella Typhi in the gallbladder remains poorly understood. Aiming to understand if S. Typhi becomes genetically adapted for long-term colonisation in the gallbladder, we performed whole genome sequencing on a collection of S. Typhi isolated from the gallbladders of typhoid carriers. These sequences were compared to contemporaneously sampled sequences from organisms isolated from the blood of acute patients within the same population. We found that S. Typhi carriage was not restricted to any particular genotype or conformation of antimicrobial resistance genes, but was largely reflective of S. Typhi circulating in the general population. However, gallbladder isolates showed a higher genetic variability than acute isolates, with median pairwise SNP distances of 21 and 13 SNPs (p = 2.8x10-9), respectively. Within gallbladder isolates of the predominant H58 genotype, variation was associated with a higher prevalence of nonsense mutations. Notably, gallbladder isolates displayed a higher frequency of non-synonymous mutations in genes encoding hypothetical proteins, membrane lipoproteins, transport/binding proteins, surface antigens, and carbohydrate degradation. Specifically, we identified several gallbladder-specific non-synonymous mutations involved in LPS synthesis and modification, with some isolates lacking the Vi capsular polysaccharide vaccine target due to the 134Kb deletion of SPI-7. S. Typhi is under strong selective pressure in the human gallbladder, which may be reflected phylogenetically by long terminal branches that may distinguish organisms from chronic and acute infections. Our work shows that selective pressures asserted by the hostile environment of the human gallbladder generate new antigenic variants and raises questions regarding the role of carriage in the epidemiology of typhoid fever.

The First 100 Days of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Control in Vietnam.

BACKGROUND: One hundred days after severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was first reported in Vietnam on 23 January, 270 cases were confirmed, with no deaths. We describe the control measures used by the government and their relationship with imported and domestically acquired case numbers, with the aim of identifying the measures associated with successful SARS-CoV-2 control. METHODS: Clinical and demographic data on the first 270 SARS-CoV-2 infected cases and the timing and nature of government control measures, including numbers of tests and quarantined individuals, were analyzed. Apple and Google mobility data provided proxies for population movement. Serial intervals were calculated from 33 infector-infectee pairs and used to estimate the proportion of presymptomatic transmission events and time-varying reproduction numbers. RESULTS: A national lockdown was implemented between 1 and 22 April. Around 200 000 people were quarantined and 266 122 reverse transcription polymerase chain reaction (RT-PCR) tests conducted. Population mobility decreased progressively before lockdown. In total, 60% (163/270) of cases were imported; 43% (89/208) of resolved infections remained asymptomatic for the duration of infection. The serial interval was 3.24 days, and 27.5% (95% confidence interval [CI], 15.7%-40.0%) of transmissions occurred presymptomatically. Limited transmission amounted to a maximum reproduction number of 1.15 (95% CI, .·37-2.·36). No community transmission has been detected since 15 April. CONCLUSIONS: Vietnam has controlled SARS-CoV-2 spread through the early introduction of mass communication, meticulous contact tracing with strict quarantine, and international travel restrictions. The value of these interventions is supported by the high proportion of asymptomatic and imported cases, and evidence for substantial presymptomatic transmission.

Novel Mutation of SARS-CoV-2, Vietnam, July 2020.

A cluster of severe acute respiratory syndrome coronavirus 2 infections in Danang, Vietnam, began July 25, 2020, and resulted in 551 confirmed cases and 35 deaths as of February 2021. We analyzed 26 sequences from this cluster and identified a novel shared mutation in nonstructural protein 9, suggesting a single introduction into Vietnam.

Emerging Coxsackievirus A6 Causing Hand, Foot and Mouth Disease, Vietnam.

Hand, foot and mouth disease (HFMD) is a major public health issue in Asia and has global pandemic potential. Coxsackievirus A6 (CV-A6) was detected in 514/2,230 (23%) of HFMD patients admitted to 3 major hospitals in southern Vietnam during 2011-2015. Of these patients, 93 (18%) had severe HFMD. Phylogenetic analysis of 98 genome sequences revealed they belonged to cluster A and had been circulating in Vietnam for 2 years before emergence. CV-A6 movement among localities within Vietnam occurred frequently, whereas viral movement across international borders appeared rare. Skyline plots identified fluctuations in the relative genetic diversity of CV-A6 corresponding to large CV-A6-associated HFMD outbreaks worldwide. These data show that CV-A6 is an emerging pathogen and emphasize the necessity of active surveillance and understanding the mechanisms that shape the pathogen evolution and emergence, which is essential for development and implementation of intervention strategies.

A universal genome sequencing method for rotavirus A from human fecal samples which identifies segment reassortment and multi-genotype mixed infection.

BACKGROUND: Genomic characterization of rotavirus (RoV) has not been adopted at large-scale due to the complexity of obtaining sequences for all 11 segments, particularly when feces are used as starting material. METHODS: To overcome these limitations, we developed a novel RoV capture and genome sequencing method combining commercial enzyme immunoassay plates and a set of routinely used reagents. RESULTS: Our approach had a 100% success rate, producing >90% genome coverage for diverse RoV present in fecal samples (Ct < 30). CONCLUSIONS: This method provides a novel, reproducible and comparatively simple approach for genomic RoV characterization and could be scaled-up for use in global RoV surveillance systems. TRIAL REGISTRATION (PROSPECTIVELY REGISTERED): Current Controlled Trials ISRCTN88101063 . Date of registration: 14/06/2012.

A Clinical and Epidemiological Investigation of the First Reported Human Infection With the Zoonotic Parasite Trypanosoma evansi in Southeast Asia.

BACKGROUND: Trypanosomais a genus of unicellular parasitic flagellate protozoa.Trypanosoma bruceispecies and Trypanosoma cruziare the major agents of human trypanosomiasis; other Trypanosomaspecies can cause human disease, but are rare. In March 2015, a 38-year-old woman presented to a healthcare facility in southern Vietnam with fever, headache, and arthralgia. Microscopic examination of blood revealed infection with Trypanosoma METHODS: Microscopic observation, polymerase chain reaction (PCR) amplification of blood samples, and serological testing were performed to identify the infecting species. The patient's blood was screened for the trypanocidal protein apolipoprotein L1 (APOL1), and a field investigation was performed to identify the zoonotic source. RESULTS: PCR amplification and serological testing identified the infecting species as Trypanosoma evansi.Despite relapsing 6 weeks after completing amphotericin B therapy, the patient made a complete recovery after 5 weeks of suramin. The patient was found to have 2 wild-type APOL1 alleles and a normal serum APOL1 concentration. After responsive animal sampling in the presumed location of exposure, cattle and/or buffalo were determined to be the most likely source of the infection, with 14 of 30 (47%) animal blood samples testing PCR positive forT. evansi. CONCLUSIONS: We report the first laboratory-confirmed case ofT. evansiin a previously healthy individual without APOL1 deficiency, potentially contracted via a wound while butchering raw beef, and successfully treated with suramin. A linked epidemiological investigation revealed widespread and previously unidentified burden ofT. evansiin local cattle, highlighting the need for surveillance of this infection in animals and the possibility of further human cases.

South Asia as a Reservoir for the Global Spread of Ciprofloxacin-Resistant Shigella sonnei: A Cross-Sectional Study.

BACKGROUND: Antimicrobial resistance is a major issue in the Shigellae, particularly as a specific multidrug-resistant (MDR) lineage of Shigella sonnei (lineage III) is becoming globally dominant. Ciprofloxacin is a recommended treatment for Shigella infections. However, ciprofloxacin-resistant S. sonnei are being increasingly isolated in Asia and sporadically reported on other continents. We hypothesized that Asia is a primary hub for the recent international spread of ciprofloxacin-resistant S. sonnei. METHODS AND FINDINGS: We performed whole-genome sequencing on a collection of 60 contemporaneous ciprofloxacin-resistant S. sonnei isolated in four countries within Asia (Vietnam, n = 11; Bhutan, n = 12; Thailand, n = 1; Cambodia, n = 1) and two outside of Asia (Australia, n = 19; Ireland, n = 16). We reconstructed the recent evolutionary history of these organisms and combined these data with their geographical location of isolation. Placing these sequences into a global phylogeny, we found that all ciprofloxacin-resistant S. sonnei formed a single clade within a Central Asian expansion of lineage III. Furthermore, our data show that resistance to ciprofloxacin within S. sonnei may be globally attributed to a single clonal emergence event, encompassing sequential gyrA-S83L, parC-S80I, and gyrA-D87G mutations. Geographical data predict that South Asia is the likely primary source of these organisms, which are being regularly exported across Asia and intercontinentally into Australia, the United States and Europe. Our analysis was limited by the number of S. sonnei sequences available from diverse geographical areas and time periods, and we cannot discount the potential existence of other unsampled reservoir populations of antimicrobial-resistant S. sonnei. CONCLUSIONS: This study suggests that a single clone, which is widespread in South Asia, is likely driving the current intercontinental surge of ciprofloxacin-resistant S. sonnei and is capable of establishing endemic transmission in new locations. Despite being limited in geographical scope, our work has major implications for understanding the international transfer of antimicrobial-resistant pathogens, with S. sonnei acting as a tractable model for studying how antimicrobial-resistant Gram-negative bacteria spread globally.

Evaluation of Luminex xTAG Gastrointestinal Pathogen Panel Assay for Detection of Multiple Diarrheal Pathogens in Fecal Samples in Vietnam.

Diarrheal disease is a complex syndrome that remains a leading cause of global childhood morbidity and mortality. The diagnosis of enteric pathogens in a timely and precise manner is important for making treatment decisions and informing public health policy, but accurate diagnosis is a major challenge in industrializing countries. Multiplex molecular diagnostic techniques may represent a significant improvement over classical approaches. We evaluated the Luminex xTAG gastrointestinal pathogen panel (GPP) assay for the detection of common enteric bacterial and viral pathogens in Vietnam. Microbiological culture and real-time PCR were used as gold standards. The tests were performed on 479 stool samples collected from people admitted to the hospital for diarrheal disease throughout Vietnam. Sensitivity and specificity were calculated for the xTAG GPP for the seven principal diarrheal etiologies. The sensitivity and specificity for the xTAG GPP were >88% for Shigellaspp.,Campylobacterspp., rotavirus, norovirus genotype 1/2 (GI/GII), and adenovirus compared to those of microbiological culture and/or real-time PCR. However, the specificity was low (∼60%) for Salmonella species. Additionally, a number of important pathogens that are not identified in routine hospital procedures in this setting, such as Cryptosporidiumspp. and Clostridium difficile, were detected with the GPP. The use of the Luminex xTAG GPP for the detection of enteric pathogens in settings, like Vietnam, would dramatically improve the diagnostic accuracy and capacity of hospital laboratories, allowing for timely and appropriate therapy decisions and a wider understanding of the epidemiology of pathogens associated with severe diarrheal disease in low-resource settings.

Host and viral features of human dengue cases shape the population of infected and infectious Aedes aegypti mosquitoes.

Dengue is the most prevalent arboviral disease of humans. The host and virus variables associated with dengue virus (DENV) transmission from symptomatic dengue cases (n = 208) to Aedes aegypti mosquitoes during 407 independent exposure events was defined. The 50% mosquito infectious dose for each of DENV-1-4 ranged from 6.29 to 7.52 log10 RNA copies/mL of plasma. Increasing day of illness, declining viremia, and rising antibody titers were independently associated with reduced risk of DENV transmission. High early DENV plasma viremia levels in patients were a marker of the duration of human infectiousness, and blood meals containing high concentrations of DENV were positively associated with the prevalence of infectious mosquitoes 14 d after blood feeding. Ambulatory dengue cases had lower viremia levels compared with hospitalized dengue cases but nonetheless at levels predicted to be infectious to mosquitoes. These data define serotype-specific viremia levels that vaccines or drugs must inhibit to prevent DENV transmission.

Deep sequencing of norovirus genomes defines evolutionary patterns in an urban tropical setting.

UNLABELLED: Norovirus is a highly transmissible infectious agent that causes epidemic gastroenteritis in susceptible children and adults. Norovirus infections can be severe and can be initiated from an exceptionally small number of viral particles. Detailed genome sequence data are useful for tracking norovirus transmission and evolution. To address this need, we have developed a whole-genome deep-sequencing method that generates entire genome sequences from small amounts of clinical specimens. This novel approach employs an algorithm for reverse transcription and PCR amplification primer design using all of the publically available norovirus sequence data. Deep sequencing and de novo assembly were used to generate norovirus genomes from a large set of diarrheal patients attending three hospitals in Ho Chi Minh City, Vietnam, over a 2.5-year period. Positive-selection analysis and direct examination of protein changes in the virus over time identified codons in the regions encoding proteins VP1, p48 (NS1-2), and p22 (NS4) under positive selection and expands the known targets of norovirus evolutionary pressure. IMPORTANCE: The high transmissibility and rapid evolutionary rate of norovirus, combined with a short-lived host immune responses, are thought to be the reasons why the virus causes the majority of pediatric viral diarrhea cases. The evolutionary patterns of this RNA virus have been described in detail for only a portion of the virus genome and never for a virus from a detailed urban tropical setting. We provide a detailed sequence description of the noroviruses circulating in three Ho Chi Minh City hospitals over a 2.5-year period. This study identified patterns of virus change in known sites of host immune response and identified three additional regions of the virus genome under selection that were not previously recognized. In addition, the method described here provides a robust full-genome sequencing platform for community-based virus surveillance.

The evolutionary consequences of blood-stage vaccination on the rodent malaria Plasmodium chabaudi.

Malaria vaccine developers are concerned that antigenic escape will erode vaccine efficacy. Evolutionary theorists have raised the possibility that some types of vaccine could also create conditions favoring the evolution of more virulent pathogens. Such evolution would put unvaccinated people at greater risk of severe disease. Here we test the impact of vaccination with a single highly purified antigen on the malaria parasite Plasmodium chabaudi evolving in laboratory mice. The antigen we used, AMA-1, is a component of several candidate malaria vaccines currently in various stages of trials in humans. We first found that a more virulent clone was less readily controlled by AMA-1-induced immunity than its less virulent progenitor. Replicated parasites were then serially passaged through control or AMA-1 vaccinated mice and evaluated after 10 and 21 rounds of selection. We found no evidence of evolution at the ama-1 locus. Instead, virulence evolved; AMA-1-selected parasites induced greater anemia in naïve mice than both control and ancestral parasites. Our data suggest that recombinant blood stage malaria vaccines can drive the evolution of more virulent malaria parasites.

Novel porcine-like human G26P[19] rotavirus identified in hospitalized paediatric diarrhoea patients in Ho Chi Minh City, Vietnam.

During a hospital-based diarrhoeal disease study conducted in Ho Chi Minh City, Vietnam from 2009 to 2010, we identified four symptomatic children infected with G26P[19] rotavirus (RV)--an atypical variant that has not previously been reported in human gastroenteritis. To determine the genetic structure and investigate the origin of this G26P[19] strain, the whole genome of a representative example was characterized, revealing a novel genome constellation: G26-P[19]-I5-R1-C1-M1-A8-N1-T1-E1-H1. The genome segments were most closely related to porcine (VP7, VP4, VP6 and NSP1) and Wa-like porcine RVs (VP1-3 and NSP2-5). We proposed that this G26P[19] strain was the product of zoonotic transmission coupled with one or more reassortment events occurring in human and/or animal reservoirs. The identification of such strains has potential implications for vaccine efficacy in south-east Asia, and outlines the utility of whole-genome sequencing for studying RV diversity and zoonotic potential during disease surveillance.

Genomic insights unveil the plasmid transfer mechanism and epidemiology of hypervirulent Klebsiella pneumoniae in Vietnam.

Hypervirulent Klebsiella pneumoniae (hvKp) is a significant cause of severe invasive infections in Vietnam, yet data on its epidemiology, population structure and dynamics are scarce. We screened hvKp isolates from patients with bloodstream infections (BSIs) at a tertiary infectious diseases hospital in Vietnam and healthy individuals, followed by whole genome sequencing and plasmid analysis. Among 700 BSI-causing Kp strains, 100 (14.3%) were hvKp. Thirteen hvKp isolates were identified from 350 rectal swabs of healthy adults; none from 500 rectal swabs of healthy children. The hvKp isolates were genetically diverse, encompassing 17 sequence types (STs), predominantly ST23, ST86 and ST65. Among the 113 hvKp isolates, 14 (12.6%) carried at least one antimicrobial resistance (AMR) gene, largely mediated by IncFII, IncR, and IncA/C plasmids. Notably, the acquisition of AMR conjugative plasmids facilitated horizontal transfer of the non-conjugative virulence plasmid between K. pneumoniae strains. Phylogenetic analysis demonstrated hvKp isolates from BSIs and human carriage clustered together, suggesting a significant role of intestinal carriage in hvKp transmission. Enhanced surveillance is crucial to understand the factors driving intestinal carriage and hvKp transmission dynamics for informing preventive measures. Furthermore, we advocate the clinical use of our molecular assay for diagnosing hvKp infections to guide effective management.

Spatiotemporal dynamics of dengue epidemics, southern Vietnam.

An improved understanding of heterogeneities in dengue virus transmission might provide insights into biological and ecologic drivers and facilitate predictions of the magnitude, timing, and location of future dengue epidemics. To investigate dengue dynamics in urban Ho Chi Minh City and neighboring rural provinces in Vietnam, we analyzed a 10-year monthly time series of dengue surveillance data from southern Vietnam. The per capita incidence of dengue was lower in Ho Chi Minh City than in most rural provinces; annual epidemics occurred 1-3 months later in Ho Chi Minh City than elsewhere. The timing and the magnitude of annual epidemics were significantly more correlated in nearby districts than in remote districts, suggesting that local biological and ecologic drivers operate at a scale of 50-100 km. Dengue incidence during the dry season accounted for 63% of variability in epidemic magnitude. These findings can aid the targeting of vector-control interventions and the planning for dengue vaccine implementation.

The clinical features and genomic epidemiology of carbapenem-resistant Acinetobacter baumannii infections at a tertiary hospital in Vietnam.

OBJECTIVES: To characterise the clinical features of Acinetobacter baumannii infections and investigate the phylogenetic structure and transmission dynamics of A. baumannii in Vietnam. METHODS: Between 2019 and 2020, a surveillance of A. baumannii (AB) infections was conducted at a tertiary hospital in Ho Chi Minh City, Vietnam. Risk factors for in-hospital mortality were analysed using logistic regressions. Whole-genome sequence data were used to characterise genomic species, sequence types (STs), antimicrobial resistance genes, surface antigens, and phylogenetic relatedness of AB isolates. RESULTS: Eighty-four patients with AB infections were enrolled in the study, 96% of whom were hospital-acquired. Half of the AB isolates were identified from ICU-admitted patients, while the remaining isolates were from non-ICU patients. The overall in-hospital mortality was 56%, with associated risk factors including advanced age, ICU stay, exposure to mechanical ventilation/central venous catheterization, pneumonia as source of AB infection, prior use of linezolid/aminoglycosides, and AB treatment with colistin-based therapy. Nearly 91% of isolates were carbapenem-resistant; 92% were multidrug-resistant; and 6% were colistin-resistant. ST2, ST571, and ST16 were the three dominant carbapenem-resistant A. baumannii (CRAB) genotypes, exhibiting distinct AMR gene profiles. Phylogenetic analysis of CRAB ST2 isolates together with previously published ST2 collection provided evidence of intra- and inter-hospital transmission of this clone. CONCLUSIONS: Our study highlights a high prevalence of carbapenem resistance and multidrug resistance in A. baumannii and elucidates the spread of CRAB within and between hospitals. Strengthening infection control measures and routine genomic surveillance are crucial to reducing the spread of CRAB and detecting novel pan-drug-resistant variants in a timely fashion.

Rhinovirus dynamics across different social structures.

Rhinoviruses (RV), common human respiratory viruses, exhibit significant antigenic diversity, yet their dynamics across distinct social structures remain poorly understood. Our study delves into RV dynamics within Kenya by analysing VP4/2 sequences across four different social structures: households, a public primary school, outpatient clinics in the Kilifi Health and Demographics Surveillance System (HDSS), and countrywide hospital admissions and outpatients. The study revealed the greatest diversity of RV infections at the countrywide level (114 types), followed by the Kilifi HDSS (78 types), the school (47 types), and households (40 types), cumulatively representing >90% of all known RV types. Notably, RV diversity correlated directly with the size of the population under observation, and several RV type variants occasionally fuelled RV infection waves. Our findings highlight the critical role of social structures in shaping RV dynamics, information that can be leveraged to enhance public health strategies. Future research should incorporate whole-genome analysis to understand fine-scale evolution across various social structures.

Genomic epidemiology, antimicrobial resistance and virulence factors of Enterobacter cloacae complex causing potential community-onset bloodstream infections in a tertiary care hospital of Nepal.

Objectives: Community-onset bloodstream infections (BSIs) caused by carbapenemase-producing Enterobacter cloacae complex (ECC) species are increasing internationally. This observation suggests that ECC are emerging pathogens, requiring for detailed understanding on their genomic epidemiology including transmission dynamics and antimicrobial resistance profiles. Patients and methods: We performed WGS on 79 Enterobacter spp. isolated from the patients with clinically significant BSIs and admitted to emergency department of a major tertiary hospital in Nepal between April 2016 and October 2017. Results: We identified 5 species and 13 STs of ECC. Enterobacter xiangfangensis ST171, one of the globally emerging carbapenem resistant ECC clones with epidemic potential, was the most prevalent (42%). Phylogenetic analysis showed a large (>19 400 SNPs) core genome SNP distance across major STs, which was minimal (<30 SNPs) among the isolates of each prevalent ST, suggesting the relatively recent importation of major STs followed by local clonal expansions. Genomic evidence for resistance to all major antimicrobial classes except for colistin and macrolides was detected. A limited number of isolates also carried bla NDM-1 (n = 2) and bla OXA-48 (n = 1) carbapenemase genes. Virulence factors encoding siderophores (24%), T6SSD (25%) and fimbriae (54%) were detected. Conclusions: Our study highlighted that MDR ECC clones are important pathogens of BSIs in community. Though of low prevalence, carbapenem resistance observed in our ECC isolates raised concern about further community dissemination, underscoring the need for community surveillance to identify MDR ECC clones with epidemic potential.

The evolutionary history of Shigella flexneri serotype 6 in Asia.

Shigella flexneri serotype 6 is an understudied cause of diarrhoeal diseases in developing countries, and has been proposed as one of the major targets for vaccine development against shigellosis. Despite being named as S. flexneri, Shigella flexneri serotype 6 is phylogenetically distinct from other S. flexneri serotypes and more closely related to S. boydii. This unique phylogenetic relationship and its low sampling frequency have hampered genomic research on this pathogen. Herein, by utilizing whole genome sequencing (WGS) and analyses of Shigella flexneri serotype 6 collected from epidemiological studies (1987-2013) in four Asian countries, we revealed its population structure and evolutionary history in the region. Phylogenetic analyses supported the delineation of Asian Shigella flexneri serotype 6 into two phylogenetic groups (PG-1 and -2). Notably, temporal phylogenetic approaches showed that extant Asian S. flexneri serotype 6 could be traced back to an inferred common ancestor arising in the 18th century. The dominant lineage PG-1 likely emerged in the 1970s, which coincided with the times to most recent common ancestors (tMRCAs) inferred from other major Southeast Asian S. flexneri serotypes. Similar to other S. flexneri serotypes in the same period in Asia, genomic analyses showed that resistance to first-generation antimicrobials was widespread, while resistance to more recent first-line antimicrobials was rare. These data also showed a number of gene inactivation and gene loss events, particularly on genes related to metabolism and synthesis of cellular appendages, emphasizing the continuing role of reductive evolution in the adaptation of the pathogen to an intracellular lifestyle. Together, our findings reveal insights into the genomic evolution of the understudied Shigella flexneri serotype 6, providing a new piece in the puzzle of Shigella epidemiology and evolution.

Evolutionary histories and antimicrobial resistance in Shigella flexneri and Shigella sonnei in Southeast Asia.

Conventional disease surveillance for shigellosis in developing country settings relies on serotyping and low-resolution molecular typing, which fails to contextualise the evolutionary history of the genus. Here, we interrogated a collection of 1,804 Shigella whole genome sequences from organisms isolated in four continental Southeast Asian countries (Thailand, Vietnam, Laos, and Cambodia) over three decades to characterise the evolution of both S. flexneri and S. sonnei. We show that S. sonnei and each major S. flexneri serotype are comprised of genetically diverse populations, the majority of which were likely introduced into Southeast Asia in the 1970s-1990s. Intranational and regional dissemination allowed widespread propagation of both species across the region. Our data indicate that the epidemiology of S. sonnei and the major S. flexneri serotypes were characterised by frequent clonal replacement events, coinciding with changing susceptibility patterns against contemporaneous antimicrobials. We conclude that adaptation to antimicrobial pressure was pivotal to the recent evolutionary trajectory of Shigella in Southeast Asia.