Emergence and clonal expansion of a qacA-harbouring sequence type 45 lineage of methicillin-resistant Staphylococcus aureus.

The past decade has seen an increase in the prevalence of sequence type (ST) 45 methicillin-resistant Staphylococcus aureus (MRSA), yet the underlying drivers for its emergence and spread remain unclear. To better understand the worldwide dissemination of ST45 S. aureus, we performed phylogenetic analyses of Australian isolates, supplemented with a global population of ST45 S. aureus genomes. Our analyses revealed a distinct lineage of multidrug-resistant ST45 MRSA harbouring qacA, predominantly found in Australia and Singapore. Bayesian inference predicted that the acquisition of qacA occurred in the late 1990s. qacA was integrated into a structurally variable region of the chromosome containing Tn552 (carrying blaZ) and Tn4001 (carrying aac(6')-aph(2")) transposable elements. Using mutagenesis and in vitro assays, we provide phenotypic evidence that qacA confers tolerance to chlorhexidine. These findings collectively suggest both antimicrobial resistance and the carriage of qacA may play a role in the successful establishment of ST45 MRSA.

Intratumoral presence of the genotoxic gut bacteria pks+ E. coli, Enterotoxigenic Bacteroides fragilis, and Fusobacterium nucleatum and their association with clinicopathological and molecular features of colorectal cancer.

BACKGROUND: This study aimed to investigate clinicopathological and molecular tumour features associated with intratumoral pks+ Escherichia coli (pks+E.coli+), pks+E.coli- (non-E.coli bacteria harbouring the pks island), Enterotoxigenic Bacteroides fragilis (ETBF) and Fusobacterium nucleatum (F. nucleatum). METHODS: We screened 1697 tumour-derived DNA samples from the Australasian Colorectal Cancer Family Registry, Melbourne Collaborative Cohort Study and the ANGELS study using targeted PCR. RESULTS: Pks+E.coli+ was associated with male sex (P < 0.01) and APC:c.835-8 A > G somatic mutation (P = 0.03). The association between pks+E.coli+ and APC:c.835-8 A > G was specific to early-onset CRCs (diagnosed<45years, P = 0.02). The APC:c.835-A > G was not associated with pks+E.coli- (P = 0.36). F. nucleatum was associated with DNA mismatch repair deficiency (MMRd), BRAF:c.1799T>A p.V600E mutation, CpG island methylator phenotype, proximal tumour location, and high levels of tumour infiltrating lymphocytes (Ps < 0.01). In the stratified analysis by MMRd subgroups, F. nucleatum was associated with Lynch syndrome, MLH1 methylated and double MMR somatic mutated MMRd subgroups (Ps < 0.01). CONCLUSION: Intratumoral pks+E.coli+ but not pks+E.coli- are associated with CRCs harbouring the APC:c.835-8 A > G somatic mutation, suggesting that this mutation is specifically related to DNA damage from colibactin-producing E.coli exposures. F. nucleatum was associated with both hereditary and sporadic MMRd subtypes, suggesting the MMRd tumour microenvironment is important for F. nucleatum colonisation irrespective of its cause.

Whole genome sequencing of drug-resistant Mycobacterium tuberculosis isolates in Victoria, Australia.

OBJECTIVES: Whole genome sequencing (WGS) can identify clusters, transmission patterns, and drug resistance mutations. This is important in low-burden settings such as Australia, as it can assist in efficient contact tracing and surveillance. METHODS: We conducted a retrospective cohort study using WGS from 155 genomically defined drug-resistant Mycobacterium tuberculosis (DR-TB) isolates collected between 2018-2021 in Victoria, Australia. Bioinformatic analysis was performed to identify resistance-conferring mutations, lineages, clusters and understand how local sequences compared with international context. RESULTS: Of the 155 sequences, 42% were identified as lineage 2 and 35% as lineage 1; 65.8% (102/155) were isoniazid mono-resistant, 8.4% were multi-drug resistant TB and 5.8% were pre-extensively drug-resistant / extensively drug-resistant TB. The most common mutations were observed in katG and fabG1 genes, especially at Ser315Thr and fabG1 -15 C>T for first-line drugs. Ser450Leu was the most frequent mutation in rpoB gene. Phylogenetic analysis confirmed that Victorian DR-TB were associated with importation events. There was little evidence of local transmission with only five isolate pairs. CONCLUSION: Isoniazid-resistant TB is the commonest DR-TB in Victoria, and the mutation profile is similar to global circulating DR-TB. Most cases are diagnosed among migrants with limited transmission. This study highlights the value of WGS in identification of clusters and resistance-conferring mutations. This information is crucial in supporting disease mitigation and treatment strategies.

Comparison of contemporary invasive and non-invasive Streptococcus pneumoniae isolates reveals new insights into circulating anti-microbial resistance determinants.

Streptococcus pneumoniae is a major human pathogen with a high burden of disease. Non-invasive isolates (those found in non-sterile sites) are thought to be a key source of invasive isolates (those found in sterile sites) and a reservoir of anti-microbial resistance (AMR) determinants. Despite this, pneumococcal surveillance has almost exclusively focused on invasive isolates. We aimed to compare contemporaneous invasive and non-invasive isolate populations to understand how they interact and identify differences in AMR gene distribution. We used a combination of whole-genome sequencing and phenotypic anti-microbial susceptibility testing and a data set of invasive (n = 1,288) and non-invasive (n = 186) pneumococcal isolates, collected in Victoria, Australia, between 2018 and 2022. The non-invasive population had increased levels of antibiotic resistance to multiple classes of antibiotics including beta-lactam antibiotics penicillin and ceftriaxone. We identified genomic intersections between the invasive and non-invasive populations and no distinct phylogenetic clustering of the two populations. However, this analysis revealed sub-populations overrepresented in each population. The sub-populations that had high levels of AMR were overrepresented in the non-invasive population. We determined that WamR-Pneumo was the most accurate in silico tool for predicting resistance to the antibiotics tested. This tool was then used to assess the allelic diversity of the penicillin-binding protein genes, which acquire mutations leading to beta-lactam antibiotic resistance, and found that they were highly conserved (≥80% shared) between the two populations. These findings show the potential of non-invasive isolates to serve as reservoirs of AMR determinants.

Public health implementation of pathogen genomics: the role for accreditation and application of ISO standards.

Pathogen genomics has transitioned rapidly from the research setting into a powerful tool now routinely used in public health microbiology, for surveillance, outbreak investigations and disease control. As these investigations can have significant public health, treatment and legal impacts, we must ensure the accuracy of these results through validation of testing processes. For laboratories working in this space, it is important to approach this work with a quality and accreditation framework in mind, working towards implementation of quality systems and test validation that meet international regulatory standards. Here we outline the key international standards and processes that lead toward accreditation for pathogen genomics.

Genomic epidemiology offers high resolution estimates of serial intervals for COVID-19.

Serial intervals - the time between symptom onset in infector and infectee - are a fundamental quantity in infectious disease control. However, their estimation requires knowledge of individuals' exposures, typically obtained through resource-intensive contact tracing efforts. We introduce an alternate framework using virus sequences to inform who infected whom and thereby estimate serial intervals. We apply our technique to SARS-CoV-2 sequences from case clusters in the first two COVID-19 waves in Victoria, Australia. We find that our approach offers high resolution, cluster-specific serial interval estimates that are comparable with those obtained from contact data, despite requiring no knowledge of who infected whom and relying on incompletely-sampled data. Compared to a published serial interval, cluster-specific serial intervals can vary estimates of the effective reproduction number by a factor of 2-3. We find that serial interval estimates in settings such as schools and meat processing/packing plants are shorter than those in healthcare facilities.

The importance of utilizing travel history metadata for informative phylogeographical inferences: a case study of early SARS-CoV-2 introductions into Australia.

Inferring the spatiotemporal spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) via Bayesian phylogeography has been complicated by the overwhelming sampling bias present in the global genomic dataset. Previous work has demonstrated the utility of metadata in addressing this bias. Specifically, the inclusion of recent travel history of SARS-CoV-2-positive individuals into extended phylogeographical models has demonstrated increased accuracy of estimates, along with proposing alternative hypotheses that were not apparent using only genomic and geographical data. However, as the availability of comprehensive epidemiological metadata is limited, many of the current estimates rely on sequence data and basic metadata (i.e. sample date and location). As the bias within the SARS-CoV-2 sequence dataset is extensive, the degree to which we can rely on results drawn from standard phylogeographical models (i.e. discrete trait analysis) that lack integrated metadata is of great concern. This is particularly important when estimates influence and inform public health policy. We compared results generated from the same dataset, using two discrete phylogeographical models: one including travel history metadata and one without. We utilized sequences from Victoria, Australia, in this case study for two unique properties. Firstly, the high proportion of cases sequenced throughout 2020 within Victoria and the rest of Australia. Secondly, individual travel history was collected from returning travellers in Victoria during the first wave (January to May) of the coronavirus disease 2019 (COVID-19) pandemic. We found that the implementation of individual travel history was essential for the estimation of SARS-CoV-2 movement via discrete phylogeography models. Without the additional information provided by the travel history metadata, the discrete trait analysis could not be fit to the data due to numerical instability. We also suggest that during the first wave of the COVID-19 pandemic in Australia, the primary driving force behind the spread of SARS-CoV-2 was viral importation from international locations. This case study demonstrates the necessity of robust genomic datasets supplemented with epidemiological metadata for generating accurate estimates from phylogeographical models in datasets that have significant sampling bias. For future work, we recommend the collection of metadata in conjunction with genomic data. Furthermore, we highlight the risk of applying phylogeographical models to biased datasets without incorporating appropriate metadata, especially when estimates influence public health policy decision making.

Population structure, serotype distribution and antibiotic resistance of Streptococcus pneumoniae causing invasive disease in Victoria, Australia.

Streptococcus pneumoniae is a major human pathogen and can cause a range of conditions from asymptomatic colonization to invasive pneumococcal disease (IPD). The epidemiology and distribution of IPD-causing serotypes in Australia has undergone large changes following the introduction of the 7-valent pneumococcal conjugate vaccine (PCV) in 2005 and the 13-valent PCV in 2011. In this study, to provide a contemporary understanding of the IPD causing population in Victoria, Australia, we aimed to examine the population structure and prevalence of antimicrobial resistance using whole-genome sequencing and comprehensive antimicrobial susceptibility data of 1288 isolates collected between 2018 and 2022. We observed high diversity among the isolates with 52 serotypes, 203 sequence types (STs) and 70 Global Pneumococcal Sequencing Project Clusters (GPSCs) identified. Serotypes contained in the 13v-PCV represented 35.3 % (n=405) of isolates. Antimicrobial resistance (AMR) to at least one antibiotic was identified in 23.8 % (n=358) of isolates with penicillin resistance the most prevalent (20.3 %, n=261 using meningitis breakpoints and 5.1 % n=65 using oral breakpoints). Of the AMR isolates, 28 % (n=101) were multidrug resistant (MDR) (resistant to three or more drug classes). Vaccination status of cases was determined for a subset of isolates with 34 cases classified as vaccine failure events (fully vaccinated IPD cases of vaccine serotype). However, no phylogenetic association with failure events was observed. Within the highly diverse IPD population, we identified six high-risk sub-populations of public health concern characterized by high prevalence, high rates of AMR and MDR, or serotype inclusion in vaccines. High-risk serotypes included serotypes 3, 19F, 19A, 14, 11A, 15A and serofamily 23. In addition, we present our data validating seroBA for in silico serotyping to facilitate ISO-accreditation of this test in routine use in a public health reference laboratory and have made this data set available. This study provides insights into the population dynamics, highlights non-vaccine serotypes of concern that are highly resistant, and provides a genomic framework for the ongoing surveillance of IPD in Australia which can inform next-generation IPD prevention strategies.

An ISO-certified genomics workflow for identification and surveillance of antimicrobial resistance.

Realising the promise of genomics to revolutionise identification and surveillance of antimicrobial resistance (AMR) has been a long-standing challenge in clinical and public health microbiology. Here, we report the creation and validation of abritAMR, an ISO-certified bioinformatics platform for genomics-based bacterial AMR gene detection. The abritAMR platform utilises NCBI's AMRFinderPlus, as well as additional features that classify AMR determinants into antibiotic classes and provide customised reports. We validate abritAMR by comparing with PCR or reference genomes, representing 1500 different bacteria and 415 resistance alleles. In these analyses, abritAMR displays 99.9% accuracy, 97.9% sensitivity and 100% specificity. We also compared genomic predictions of phenotype for 864 Salmonella spp. against agar dilution results, showing 98.9% accuracy. The implementation of abritAMR in our institution has resulted in streamlined bioinformatics and reporting pathways, and has been readily updated and re-verified. The abritAMR tool and validation datasets are publicly available to assist laboratories everywhere harness the power of AMR genomics in professional practice.

Hidden Resistances: How Routine Whole-Genome Sequencing Uncovered an Otherwise Undetected blaNDM-1 Gene in Vibrio alginolyticus from Imported Seafood.

Vibrio alginolyticus causes vibriosis of marine vertebrates, invertebrates, and humans, and while there have been several reports of multidrug resistance in V. alginolyticus, carbapenem resistance is rare. V. alginolyticus strain AUSMDU00064140 was isolated in Melbourne, Australia, from imported prawns. Routine genomic surveillance detected the presence of a full-length blaNDM-1 gene, subsequently shown to be collocated with additional acquired antimicrobial resistance genes on a resistance cassette on the largest chromosome, flanked by mobilization gene annotations. Comparisons to a previously described V. alginolyticus plasmid, pC1349, revealed differing gene content and arrangements between the resistance cassettes. Phylogenetic analysis was performed against a local and global data set (n = 109), demonstrating that AUSMDU00064140 was distinct and did not cluster with any other strains. Despite the presence of the complete blaNDM-1 gene and positive phenotypic assays for carbapenemase production, carbapenem MICs were low (meropenem MIC ≤0.5 mg/liter). However, it is still possible that this gene may be transferred to another species in the environment or a host, causing phenotypic carbapenem resistance and presenting a risk of great public health concern. IMPORTANCE Carbapenems are last-line antimicrobials, vital for use in human medicine. Antimicrobial resistance determinants such as blaNDM (New Delhi metallo-ß-lactamase producing) genes conferring resistance to the carbapenem class of antimicrobials, are typically found in Enterobacterales (first described in 2009 from a Klebsiella pneumoniae isolate). Our study shows that Vibrio alginolyticus isolated from cooked prawn is able to harbor antimicrobial resistance (AMR) genes of public health concern, specifically a chromosomally located blaNDM-1 gene, and there is the potential for transmission of resistance genes. This may be linked with antimicrobial use in low- and middle-income settings, which has typically been high, unregulated, or not reported. Many countries, including Thailand, have implemented national strategic plans to incorporate the World Health Organization (WHO)'s Global Action Plan (2015) recommendations of a global One Health approach, including increased resources for surveillance of antimicrobial usage and AMR; however, efficient antimicrobial surveillance systems incorporating genomic and phenotypic testing of isolates are still lacking in many jurisdictions.

State-wide genomic epidemiology investigations of COVID-19 in healthcare workers in 2020 Victoria, Australia: Qualitative thematic analysis to provide insights for future pandemic preparedness.

Background: COVID-19 has affected many healthcare workers (HCWs) globally. We performed state-wide SARS-CoV-2 genomic epidemiological investigations to identify HCW transmission dynamics and provide recommendations to optimise healthcare system preparedness for future outbreaks. Methods: Genome sequencing was attempted on all COVID-19 cases in Victoria, Australia. We combined genomic and epidemiologic data to investigate the source of HCW infections across multiple healthcare facilities (HCFs) in the state. Phylogenetic analysis and fine-scale hierarchical clustering were performed for the entire dataset including community and healthcare cases. Facilities provided standardised epidemiological data and putative transmission links. Findings: Between March-October 2020, approximately 1,240 HCW COVID-19 infection cases were identified; 765 are included here, requested for hospital investigations. Genomic sequencing was successful for 612 (80%) cases. Thirty-six investigations were undertaken across 12 HCFs. Genomic analysis revealed that multiple introductions of COVID-19 into facilities (31/36) were more common than single introductions (5/36). Major contributors to HCW acquisitions included mobility of staff and patients between wards and facilities, and characteristics and behaviours of patients that generated numerous secondary infections. Key limitations at the HCF level were identified. Interpretation: Genomic epidemiological analyses enhanced understanding of HCW infections, revealing unsuspected clusters and transmission networks. Combined analysis of all HCWs and patients in a HCF should be conducted, supported by high rates of sequencing coverage for all cases in the population. Established systems for integrated genomic epidemiological investigations in healthcare settings will improve HCW safety in future pandemics. Funding: The Victorian Government, the National Health and Medical Research Council Australia, and the Medical Research Future Fund.

SARS-CoV-2 infection results in immune responses in the respiratory tract and peripheral blood that suggest mechanisms of disease severity.

Respiratory tract infection with SARS-CoV-2 results in varying immunopathology underlying COVID-19. We examine cellular, humoral and cytokine responses covering 382 immune components in longitudinal blood and respiratory samples from hospitalized COVID-19 patients. SARS-CoV-2-specific IgM, IgG, IgA are detected in respiratory tract and blood, however, receptor-binding domain (RBD)-specific IgM and IgG seroconversion is enhanced in respiratory specimens. SARS-CoV-2 neutralization activity in respiratory samples correlates with RBD-specific IgM and IgG levels. Cytokines/chemokines vary between respiratory samples and plasma, indicating that inflammation should be assessed in respiratory specimens to understand immunopathology. IFN-α2 and IL-12p70 in endotracheal aspirate and neutralization in sputum negatively correlate with duration of hospital stay. Diverse immune subsets are detected in respiratory samples, dominated by neutrophils. Importantly, dexamethasone treatment does not affect humoral responses in blood of COVID-19 patients. Our study unveils differential immune responses between respiratory samples and blood, and shows how drug therapy affects immune responses during COVID-19.

Multi-site implementation of whole genome sequencing for hospital infection control: A prospective genomic epidemiological analysis.

Background: Current microbiological methods lack the resolution to accurately identify multidrug-resistant organism (MDRO) transmission, however, whole genome sequencing can identify highly-related patient isolates providing opportunities for precision infection control interventions. We investigated the feasibility and potential impact of a prospective multi-centre genomics workflow for hospital infection control. Methods: We conducted a prospective genomics implementation study across eight Australian hospitals over 15 months (2017,2018), collecting all clinical and screening isolates from inpatients with vanA VRE, MRSA, ESBL Escherichia coli (ESBL-Ec), or ESBL Klebsiella pneumoniae (ESBL-Kp). Genomic and epidemiologic data were integrated to assess MDRO transmission. Findings: In total, 2275 isolates were included from 1970 patients, predominantly ESBL-Ec (40·8%) followed by MRSA (35·6%), vanA VRE (15·2%), and ESBL-Kp (8·3%).Overall, hospital and genomic epidemiology showed 607 patients (30·8%) acquired their MDRO in hospital, including the majority of vanA VRE (266 patients, 86·4%), with lower proportions of ESBL-Ec (186 patients, 23·0%), ESBL-Kp (42 patients, 26·3%), and MRSA (113 patients, 16·3%). Complex patient movements meant the majority of MDRO transmissions would remain undetected without genomic data.The genomics implementation had major impacts, identifying unexpected MDRO transmissions prompting new infection control interventions, and contributing to vanA VRE becoming a notifiable condition. We identified barriers to implementation and recommend strategies for mitigation. Interpretation: Implementation of a multi-centre genomics-informed infection control workflow is feasible and identifies many unrecognised MDRO transmissions. This provides critical opportunities for interventions to improve patient safety in hospitals. Funding: Melbourne Genomics Health Alliance (supported by State Government of Victoria, Australia), and National Health and Medical Research Council (Australia).

Optimising genomic approaches for identifying vancomycin-resistant Enterococcus faecium transmission in healthcare settings.

Vancomycin-resistant Enterococcus faecium (VREfm) is a major nosocomial pathogen. Identifying VREfm transmission dynamics permits targeted interventions, and while genomics is increasingly being utilised, methods are not yet standardised or optimised for accuracy. We aimed to develop a standardized genomic method for identifying putative VREfm transmission links. Using comprehensive genomic and epidemiological data from a cohort of 308 VREfm infection or colonization cases, we compared multiple approaches for quantifying genetic relatedness. We showed that clustering by core genome multilocus sequence type (cgMLST) was more informative of population structure than traditional MLST. Pairwise genome comparisons using split k-mer analysis (SKA) provided the high-level resolution needed to infer patient-to-patient transmission. The more common mapping to a reference genome was not sufficiently discriminatory, defining more than three times more genomic transmission events than SKA (3729 compared to 1079 events). Here, we show a standardized genomic framework for inferring VREfm transmission that can be the basis for global deployment of VREfm genomics into routine outbreak detection and investigation.

Daptomycin Resistance Occurs Predominantly in vanA-Type Vancomycin-Resistant Enterococcus faecium in Australasia and Is Associated With Heterogeneous and Novel Mutations.

Healthcare associated infections caused by vancomycin-resistant Enterococcus faecium (VREfm) have a major impact on health outcomes. VREfm is difficult to treat because of intrinsic and acquired resistance to many clinically used antimicrobials, with daptomycin being one of the few last line therapeutic options for treating multidrug-resistant VREfm. The emergence of daptomycin-resistant VREfm is therefore of serious clinical concern. Despite this, the impact that daptomycin-resistant VREfm have on patient health outcomes is not clearly defined and knowledge on the mechanisms and genetic signatures linked with daptomycin resistance in VREfm remains incomplete. To address these knowledge gaps, phenotypic daptomycin susceptibility testing was undertaken on 324 E. faecium isolates from Australia and New Zealand. Approximately 15% of study isolates were phenotypically resistant to daptomycin. Whole genome sequencing revealed a strong association between vanA-VREfm and daptomycin resistance, with 95% of daptomycin-resistant study isolates harbouring vanA. Genomic analyses showed that daptomycin-resistant VREfm isolates were polyclonal and carried several previously characterised mutations in the liaR and liaS genes as well as several novel mutations within the rpoB, rpoC, and dltC genes. Overall, 70% of daptomycin-resistant study isolates were found to carry mutations within the liaR, rpoB, rpoC, or dltC genes. Finally, in a mouse model of VREfm bacteraemia, infection with the locally dominant daptomycin-resistant clone led to reduced daptomycin treatment efficacy in comparison to daptomycin-susceptible E. faecium. These findings have important implications for ongoing VREfm surveillance activities and the treatment of VREfm infections.

Defective Severe Acute Respiratory Syndrome Coronavirus 2 Immune Responses in an Immunocompromised Individual With Prolonged Viral Replication.

We describe severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-specific immune responses in a patient with lymphoma and recent programmed death 1 (PD-1) inhibitor therapy with late onset of severe coronavirus disease 2019 disease and prolonged SARS-CoV-2 replication, in comparison to age-matched and immunocompromised controls. High levels of HLA-DR+/CD38+ activation, interleukin 6, and interleukin 18 in the absence of B cells and PD-1 expression was observed. SARS-CoV-2-specific antibody responses were absent and SARS-CoV-2-specific T cells were minimally detected. This case highlights challenges in managing immunocompromised hosts who may fail to mount effective virus-specific immune responses.

Genomic Epidemiology and Antimicrobial Resistance Mechanisms of Imported Typhoid in Australia.

Typhoid fever is an invasive bacterial disease of humans that disproportionately affects low- and middle-income countries. Antimicrobial resistance (AMR) has been increasingly prevalent in recent decades in Salmonella enterica serovar Typhi (S. Typhi), the causative agent of typhoid fever, limiting treatment options. In Australia, most cases of typhoid fever are imported due to travel to regions where typhoid fever is endemic. Here, all 116 isolates of S. Typhi isolated in Victoria, Australia, between 1 July 2018 and 30 June 2020, underwent whole-genome sequencing and antimicrobial susceptibility testing. Genomic data were linked to international travel data collected from routine case interviews. Travel to South Asia accounted for most cases, with 92.2% imported from seven primary countries (the top two were India, n = 87, and Pakistan, n = 12). A total of 17 S. Typhi genotypes were detected in the 2-year cohort, with 48.2% genotyped as part of global AMR lineages. Ciprofloxacin resistance was detected in two lineages, 3.3 and 4.3.1.2, all from cases with reported travel to India. Nearly all multidrug and extensively drug resistant isolates (90%) were from cases with reported travel to Pakistan in genotypes 4.3.1.1 and 4.3.1.1.P1. Extended spectrum beta-lactamases, blaCTX-M-15 and blaSHV-12, were detected in cases with travel to Pakistan and India, respectively. Linking epidemiological data with genomic studies of S. Typhi provides an opportunity to improve understanding of the emergence, spread and risk of drug-resistant S. Typhi infections and to better inform empirical treatment guidelines in returned travelers.