Search

BACKGROUND: Patients with prolonged hospitalisation have a significant risk of carriage of and subsequent infection with extended spectrum β-lactamase (ESBL)-producing and carbapenemase-producing Klebsiella pneumoniae. However, the distinctive roles of the community and hospital environments in the transmission of ESBL-producing or carbapenemase-producing K pneumoniae remain elusive. We aimed to investigate the prevalence and transmission of K pneumoniae within and between the two tertiary hospitals in Hanoi, Viet Nam, using whole-genome sequencing. METHODS: We did a prospective cohort study of 69 patients in intensive care units (ICUs) from two hospitals in Hanoi, Viet Nam. Patients were included if they were aged 18 years or older, admitted for longer than the mean length of stay in their ICU, and cultured K pneumoniae from their clinical samples. Longitudinally collected samples from patients (collected weekly) and the ICU environment (collected monthly) were cultured on selective media, and whole-genome sequences from K pneumoniae colonies analysed. We did phylogenetic analyses and correlated phenotypic antimicrobial susceptibility testing with genotypic features of K pneumoniae isolates. We constructed transmission networks of patient samples, relating ICU admission times and locations with genetic similarity of infecting K pneumoniae. FINDINGS: Between June 1, 2017, and Jan 31, 2018, 69 patients were in the ICUs and eligible for inclusion, and a total of 357 K pneumoniae isolates were cultured and successfully sequenced. 228 (64%) of K pneumoniae isolates carried between two and four different ESBL-encoding and carbapenemase-encoding genes, with 164 (46%) isolates carrying genes encoding both, with high minimum inhibitory concentrations. We found a novel co-occurrence of blaKPC-2 and blaNDM-1 in 46·6% of samples from the globally successful ST15 lineage. Despite being physically and clinically separated, the two hospitals shared closely related strains carrying the same array of antimicrobial resistance genes. INTERPRETATION: These results highlight the high prevalence of ESBL-positive carbapenem-resistant K pneumoniae in ICUs in Viet Nam. Through studying K pneumoniae ST15 in detail, we showed how important resistance genes are contained within these strains that are carried broadly by patients entering the two hospitals directly or through referral. FUNDING: Medical Research Council Newton Fund, Ministry of Science and Technology, Wellcome Trust, Academy of Medical Sciences, Health Foundation, and National Institute for Health and Care Research Cambridge Biomedical Research Centre.

BackgroundSyphilis is a sexually transmitted bacterial infection caused byTreponema pallidumsubspeciespallidum, with approximately 6.3 million annual cases globally. Over the last decade, syphilis rates have risen dramatically in many high-income countries, including in England, which has seen a greater than 150% increase. Although this increase is known to be associated with high risk sexual activity in gay, bisexual and other men who have sex with men (GBMSM), cases are rising in heterosexual men and women, and congenital syphilis cases are now seen annually. The transmission dynamics within and between sexual networks of GBMSM and heterosexuals are not well understood.MethodsTo determine if whole genome sequencing could be used to identify discrete patterns of transmission, we linked national patient demographic, geospatial and behavioural metadata to wholeT. pallidumgenome sequences previously generated from 237 patient samples collected from across England between 2012 and 2018.FindingsPhylogenomic analysis and clustering revealed two of the eightT. pallidumsublineages detected in England dominated. These dominant sublineages exhibited different spatiotemporal trends linked to demography or behaviour, suggesting they represent different sexual networks: sublineage 1 was found throughout England and across all patient groups, whereas sublineage 14 occurred predominantly in older GBMSM and was absent from samples sequenced from the North of England. By focussing on different regions of England we were able to distinguish a local heterosexual transmission cluster from a background of transmission amongst GBMSM.InterpretationThese findings demonstrate that despite extremely close genetic relationships betweenT. pallidumgenomes globally, genomics can still be used to identify putative transmission clusters for epidemiological follow-up, and therefore has a role to play in informing public health interventions.FundingWellcome funding to the Sanger Institute (#206194 and 108413/A/15/D), UKRI and NIHR (COV0335; MR/V027956/1, NIHR200125), the EDCTP (RIA2018D-249), and UKHSA.Research in ContextEvidence before this studyDetailed phylogenomic analyses investigating the epidemiology and transmission dynamics ofTreponema pallidumare challenging due to low bacterial loads in clinical specimens, and difficulty in culturing the bacteria. We searched PubMed until August 9th2022 using the search terms “Syphilis” or “Treponema pallidum” and “genomic” or “genome(s)” or “sequencing”, finding 23 studies describing whole genome sequencing ofT. pallidumsubspeciespallidum, of which two used whole genome phylogenies to investigate sexual network epidemiology, with one large study of sexual networks conducted primarily in Victoria, Australia which characterised two major circulating sublineages in that setting, as well as putative sexual transmission networks with distinct sexual behavioural characteristics and potential bridging between networks.Added value of this studyIn this study, we linked national surveillance data toT. pallidumgenomes, and characterised the transmission dynamics of syphilis using samples from across a whole country, in a European setting (England). Integration of national-level sociodemographic, spatiotemporal and genomic data allowed the delineation of putative sexual networks at both the national and region levels, and revealed patterns not previously detected using epidemiological or genomic data alone.Implications of all the available evidenceOur findings are consistent with findings in Australia that demonstrate genomics can identify putative sociodemographic transmission clusters. However, in that study genomic clusters included samples separated by multiple single nucleotide polymorphisms, which could represent several years of evolution. Our study explored the value of linking identical genomes, and highlights that despite technical constraints, whole genome sequencing can be used to enable outbreak exclusion and identify putative local transmission clusters for epidemiological follow-up.

The humanitarian crisis in Yemen led in 2016 to the biggest cholera outbreak documented in modern history, with more than 2.5 million suspected cases to date. In late 2018, epidemiological surveillance showed that V. cholerae isolated from cholera patients had turned multi-drug resistant (MDR). We generated genomes from 260 isolates sampled in Yemen between 2018 and 2019 to identify a possible shift in circulating genotypes. 84% of V. cholerae isolates were serogroup O1 belonging to the seventh pandemic El Tor (7PET) lineage, sublineage T13 – same as in 2016 and 2017 – while the remaining 16% of strains were non-toxigenic and belonged to divergent V. cholerae lineages, likely reflecting sporadic gut colonisation by endemic strains. Phylogenomic analysis reveals a succession of T13 clones, with 2019 dominated by a clone that carried an IncC-type plasmid harbouring an MDR pseudo-compound transposon (PCT). Identical copies of these mobile elements were found independently in several unrelated lineages, suggesting exchange and recombination between endemic and epidemic strains. Treatment of severe cholera patients with macrolides in Yemen from 2016 to early 2019 coincides with the emergence of the plasmid-carrying T13 clone. The unprecedented success of this genotype where an SXT-family integrative and conjugative element (SXT/ICE) and an IncC plasmid coinhabit show the stability of this MDR plasmid in the 7PET background, which may durably reduce options for epidemic cholera case management. We advocate a heightened genomic epidemiology surveillance of cholera to help control the spread of this highly-transmissible, MDR clone.

Members of the Klebsiella pneumoniae species complex, particularly K. pneumoniae subsp. pneumoniae are antimicrobial resistance (AMR) associated pathogens of global importance, and polyvalent vaccines targeting Klebsiella O-antigens are in development. Whole-genome sequencing has provided insight into O-antigen distribution in the K. pneumoniae species complex, as well as population structure and virulence determinants, but genomes from sub-Saharan Africa are underrepresented in global sequencing efforts. We therefore carried out a genomic analysis of extended-spectrum beta-lactamase (ESBL)-producing K. pneumoniae species complex isolates colonizing adults in Blantyre, Malawi. We placed these isolates in a global genomic context, and compared colonizing to invasive isolates from the main public hospital in Blantyre. In total, 203 isolates from stool and rectal swabs from adults were whole-genome sequenced and compared to a publicly available multicounty collection and previously sequenced Malawian and Kenyan isolates from blood or sterile sites. We inferred phylogenetic relationships and analysed the diversity of genetic loci linked to AMR, virulence, capsule and LPS O-antigen (O-types). We find that the diversity of Malawian K. pneumoniae subsp. pneumoniae isolates represents the species’ population structure, but shows distinct local signatures concerning clonal expansions. Siderophore and hypermucoidy genes were more frequent in invasive versus colonizing isolates (present in 13 % vs 1 %) but still generally lacking in most invasive isolates. O-antigen population structure and distribution was similar in invasive and colonizing isolates, with O4 more common (14%) than in previously published studies (2–5 %). We conclude that host factors, pathogen opportunity or alternate virulence loci not linked to invasive disease elsewhere are likely to be the major determinants of invasive disease in Malawi. Distinct ST and O-type distributions in Malawi highlight the need to sample locations where the burden of invasive Klebsiella disease is greatest to robustly define secular trends in Klebsiella diversity to assist in the development of a useful vaccine. Colonizing and invasive isolates in Blantyre are similar, hence O-typing of colonizing Klebsiella isolates may be a rapid and cost-effective approach to describe global diversity and guide vaccine development.

BACKGROUND: Treponema pallidum subspecies pertenue causes yaws. Strategies to better control, eliminate, and eradicate yaws are needed. METHODS: In an open-label, cluster-randomized, community-based trial conducted in a yaws-endemic area of Papua New Guinea, we randomly assigned 38 wards (i.e., clusters) to receive one round of mass administration of azithromycin followed by two rounds of target treatment of active cases (control group) or three rounds of mass administration of azithromycin (experimental group); round 1 was administered at baseline, round 2 at 6 months, and round 3 at 12 months. The coprimary end points were the prevalence of active cases of yaws, confirmed by polymerase-chain-reaction assay, in the entire trial population and the prevalence of latent yaws, confirmed by serologic testing, in a subgroup of asymptomatic children 1 to 15 years of age; prevalences were measured at 18 months, and the between-group differences were calculated. RESULTS: Of the 38 wards, 19 were randomly assigned to the control group (30,438 persons) and 19 to the experimental group (26,238 persons). A total of 24,848 doses of azithromycin were administered in the control group (22,033 were given to the participants at round 1 and 207 and 2608 were given to the participants with yaws-like lesions and their contacts, respectively, at rounds 2 and 3 [combined]), and 59,852 doses were administered in the experimental group. At 18 months, the prevalence of active yaws had decreased from 0.46% (102 of 22,033 persons) at baseline to 0.16% (47 of 29,954 persons) in the control group and from 0.43% (87 of 20,331 persons) at baseline to 0.04% (10 of 25,987 persons) in the experimental group (relative risk adjusted for clustering, 4.08; 95% confidence interval [CI], 1.90 to 8.76). The prevalence of other infectious ulcers decreased to a similar extent in the two treatment groups. The prevalence of latent yaws at 18 months was 6.54% (95% CI, 5.00 to 8.08) among 994 children in the control group and 3.28% (95% CI, 2.14 to 4.42) among 945 children in the experimental group (relative risk adjusted for clustering and age, 2.03; 95% CI, 1.12 to 3.70). Three cases of yaws with resistance to macrolides were found in the experimental group. CONCLUSIONS: The reduction in the community prevalence of yaws was greater with three rounds of mass administration of azithromycin at 6-month intervals than with one round of mass administration of azithromycin followed by two rounds of targeted treatment. Monitoring for the emergence and spread of antimicrobial resistance is needed. (Funded by Fundació "la Caixa" and others; ClinicalTrials.gov number, NCT03490123.).

Background: \ud The SARS-CoV-2 delta (B.1.617.2) variant was first detected in England in March, 2021. It has since rapidly become the predominant lineage, owing to high transmissibility. It is suspected that the delta variant is associated with more severe disease than the previously dominant alpha (B.1.1.7) variant. We aimed to characterise the severity of the delta variant compared with the alpha variant by determining the relative risk of hospital attendance outcomes.\ud \ud Methods: \ud This cohort study was done among all patients with COVID-19 in England between March 29 and May 23, 2021, who were identified as being infected with either the alpha or delta SARS-CoV-2 variant through whole-genome sequencing. Individual-level data on these patients were linked to routine health-care datasets on vaccination, emergency care attendance, hospital admission, and mortality (data from Public Health England's Second Generation Surveillance System and COVID-19-associated deaths dataset; the National Immunisation Management System; and NHS Digital Secondary Uses Services and Emergency Care Data Set). The risk for hospital admission and emergency care attendance were compared between patients with sequencing-confirmed delta and alpha variants for the whole cohort and by vaccination status subgroups. Stratified Cox regression was used to adjust for age, sex, ethnicity, deprivation, recent international travel, area of residence, calendar week, and vaccination status.\ud \ud Findings: \ud Individual-level data on 43 338 COVID-19-positive patients (8682 with the delta variant, 34 656 with the alpha variant; median age 31 years [IQR 17–43]) were included in our analysis. 196 (2·3%) patients with the delta variant versus 764 (2·2%) patients with the alpha variant were admitted to hospital within 14 days after the specimen was taken (adjusted hazard ratio [HR] 2·26 [95% CI 1·32–3·89]). 498 (5·7%) patients with the delta variant versus 1448 (4·2%) patients with the alpha variant were admitted to hospital or attended emergency care within 14 days (adjusted HR 1·45 [1·08–1·95]). Most patients were unvaccinated (32 078 [74·0%] across both groups). The HRs for vaccinated patients with the delta variant versus the alpha variant (adjusted HR for hospital admission 1·94 [95% CI 0·47–8·05] and for hospital admission or emergency care attendance 1·58 [0·69–3·61]) were similar to the HRs for unvaccinated patients (2·32 [1·29–4·16] and 1·43 [1·04–1·97]; p=0·82 for both) but the precision for the vaccinated subgroup was low.\ud \ud Interpretation: \ud This large national study found a higher hospital admission or emergency care attendance risk for patients with COVID-19 infected with the delta variant compared with the alpha variant. Results suggest that outbreaks of the delta variant in unvaccinated populations might lead to a greater burden on health-care services than the alpha variant.\ud \ud Funding: \ud Medical Research Council; UK Research and Innovation; Department of Health and Social Care; and National Institute for Health Research.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections were rising during early summer 2021 in many countries as a result of the Delta variant. We assessed reverse transcription polymerase chain reaction swab positivity in the Real-time Assessment of Community Transmission–1 (REACT-1) study in England. During June and July 2021, we observed sustained exponential growth with an average doubling time of 25 days, driven by complete replacement of the Alpha variant by Delta and by high prevalence at younger, less-vaccinated ages. Prevalence among unvaccinated people [1.21% (95% credible interval 1.03%, 1.41%)] was three times that among double-vaccinated people [0.40% (95% credible interval 0.34%, 0.48%)]. However, after adjusting for age and other variables, vaccine effectiveness for double-vaccinated people was estimated at between ~50% and ~60% during this period in England. Increased social mixing in the presence of Delta had the potential to generate sustained growth in infections, even at high levels of vaccination.

Shortening courses of antimicrobials has been proposed to reduce risk of antimicrobial resistant (AMR) infections, but acquisition and selection dynamics under antimicrobial pressure at the individual level are poorly understood. We combine multi-state modelling and whole-genome sequencing to understand colonisation dynamics of extended-spectrum beta-lactamase producing Enterobacterales (ESBL-E) in Malawian adults. We demonstrate prolonged post-exposure antibiotic effect, meaning short courses exert similar colonisation pressure to longer ones. Genome data does not identify widespread hospital-associated ESBL-E transmission, hence apparent acquisitions may be selected from the patient microbiota by antimicrobial exposure. Understanding ESBL-E dynamics under antimicrobial pressure is crucial for evidence-based stewardship protocols.

Drug-resistant bacteria of the order Enterobacterales which produce extended-spectrum beta-lactamase enzymes (ESBL-Enterobacterales, ESBL-E) are global priority pathogens. Antimicrobial stewardship interventions proposed to curb their spread include shorter courses of antimicrobials to reduce selection pressure but individual-level acquisition and selection dynamics are poorly understood. We sampled stool of 425 adults (aged 16–76 years) in Blantyre, Malawi, over 6 months and used multistate modelling and whole-genome sequencing to understand colonization dynamics of ESBL-E. Models suggest a prolonged effect of antimicrobials such that truncating an antimicrobial course at 2 days has a limited effect in reducing colonization. Genomic analysis shows largely indistinguishable diversity of healthcare-associated and community-acquired isolates, hence some apparent acquisition of ESBL-E during hospitalization may instead represent selection from a patient’s microbiota by antimicrobial exposure. Our approach could help guide stewardship protocols; interventions that aim to review and truncate courses of unneeded antimicrobials may be of limited use in preventing ESBL-E colonization.

Klebsiella pneumoniae is an antimicrobial resistance (AMR) associated pathogen of global importance, and polyvalent vaccines targeting K. pneumoniae O-antigens are in development. Genomes from sub-Saharan Africa (sSA) are underrepresented in global sequencing efforts. We therefore carried out a genomic analysis of extended-spectrum beta-lactamase (ESBL)-producing K. pneumoniae complex isolates colonising adults in Blantyre, Malawi, placed these isolates in a global genomic context, and compared colonising to invasive isolates from the main public hospital in Blantyre. 203 isolates from stool and rectal swabs from adults were whole-genome sequenced and compared to a publicly available multicountry collection of 484 K. pneumoniae genomes sampled to cover maximum diversity of the species, 150 previously sequenced Malawian and 66 Kenyan isolates from blood or sterile sites. We inferred phylogenetic relationships and analysed the diversity of genetic loci linked to AMR, virulence, capsule (K-) and LPS O-antigen (O-types). We find that the diversity of Malawian Klebsiella isolates is representative of the species’ population structure, but with local success and expansion of sequence types (STs) ST14, ST15, ST340 and ST307. Siderophore and hypermucoidy genes were more frequent in invasive versus carriage isolates (present in 13% vs 1%, p < 0.001) but still generally lacking in most invasive isolates. The population structure and distribution of O-antigen types was similar in Malawian invasive and carriage isolates, with O4 being more common in Malawian isolates (14%) than in previously published studies (2-5%). We conclude that host factors, pathogen opportunity or alternate virulence loci not linked to invasive disease elsewhere are likely to be the major determinants of invasive disease in Malawi. Distinct ST and O-type distributions in Malawi highlights the need for geographically aware sampling to robustly define secular trends in Klebsiella diversity. Colonising and invasive isolates in Blantyre are similar and hence O-typing of colonising Klebsiella isolates may be a rapid and cost-effective approach to describe global diversity and guide vaccine development.Data SummaryAll data and code to replicate this analysis is available as the blantyreESBL v1.0.0 R package (https://doi.org/10.5281/zenodo.5554082) available at https://github.com/joelewis101/blantyreESBL. Reads from all isolates sequenced as part of this study have been deposited in the European Nucleotide Archive, and accession numbers (as well as accession numbers of publicly available genomes used in this analysis) are provided in the R package.

Escherichia coli is a ubiquitous bacterium and one of the most prevalent Gram-negative species associated with drug resistant infections. The large number of sequenced genomes available have provided us with a consistently growing knowledge base to further understand pathogenesis and epidemiology of this organism. However, data from sub-Saharan Africa (sSA) are underrepresented in global sequencing efforts and E. coli genetic diversity from this region is poorly described. To reduce this gap, we investigated extended-spectrum beta-lactamase (ESBL)-producing E. coli colonising adults in Blantyre, Malawi to assess the bacterial diversity and AMR determinants and to place these isolates in the context of the wider population structure. We performed short-read whole-genome sequencing of 473 colonising ESBL E. coli isolated from human stool and contextualised the genomes with a previously curated multi-country species wide collection of 10,146 genomes. The most frequently identified sequence types (STs) in our collection were the globally successful ST131, ST410 and ST167, and the dominant ESBL genes were blaCTX-M, mirroring global trends. However, 37% of Malawian isolates did not cluster with any isolates in the curated multicountry collection and a core gene phylogeny was consistent with locally spreading subclades within globally dominant clones, including in ST410 and ST167. We also found Carbapenemase genes in our collection at low frequency; we used long read sequencing to characterise selected ESBL and carbapenemase-associated plasmids, demonstrating the presence of globally distributed carbapenemase carrying plasmids. Increased genomic surveillance of E. coli from Malawi and sSA is necessary to understand local, regional and global transmission of both E. coli and the AMR genes they commonly carry.Impact StatementDrug-resistant Escherichia coli producing extended-spectrum beta lactamase (ESBL) or carbapenemase enzymes have been identified by the World Health Organisation as priority pathogens of global concern, and whole genome sequencing has provided insight into mechanisms of virulence, antimicrobial resistance, and the spread of high-risk clones. However, studies analysing large numbers of E. coli using whole-genome data often focus on opportunistic use of hospital diagnostic collections in high-income settings. Understanding how the genomic epidemiology of E. coli in low- and middle-income countries (including many of the nations of sub-Saharan Africa) differs is essential to provide insight into local, and global drivers of transmission. We therefore sequenced 473 ESBL-producing E. coli genomes colonising adults in Blantyre, Malawi. We analyse determinants of antimicrobial resistance and virulence and place the isolates in wider context using a previously published global E. coli collection that was generated to represent the whole species diversity of sequences publicly available at the time of generation. We find that there is diversity in Malawian isolates not reflected in the curated global collection: widely successful antimicrobial-resistance associated E. coli sequence types are represented in Blantyre, but locally circulating subclades are apparent. Furthermore, given the high number of ESBL producing pathogens causing infections there is an unmet need for carbapenem antimicrobials which are still active against ESBL-producers but are not yet widely available in our setting. We find that carbapenemases (enzymes that can render bacteria resistant to carbapenems) in our collection are unusual but present and carried on globally disseminated plasmids. So too are globally successful, stably carbapenemase-associated E. coli lineages. Although the Malawian isolates analysed typically lacked carbapenemases, carbapenem use is increasing in Malawi and their unstewarded use will accelerate selection for carbapememases in E. coli in the future. Our study highlights the need for robust stewardship protocols and ongoing genomic surveillance as these agents are introduced.Data SummaryAll data and code to replicate this analysis are available as the blantyreESBL v1.3 R package (https://doi.org/10.5281/zenodo.5554081) available at https://github.com/joelewis101/blantyreESBL. Reads from all isolates sequenced as part of this study have been deposited in the European Nucleotide Archive, under PRJEB26677, PRJEB28522 and PRJEB36486 (short reads) and PRJNA869071 (Nanopore reads and hybrid assemblies). Accession numbers (as well as accession numbers of publicly available genomes used in this analysis) linked to sample metadata are provided in the R package and as supplementary data to this manuscript.

The incidence of syphilis has increased markedly in the past decade in high-income countries, including Australia. To date, however, genomic studies of Treponema pallidum have focused mainly on the northern hemisphere. Here, we aimed to characterise the lineages of T pallidum driving the current syphilis epidemic in Australia.In this genomic epidemiological analysis, using phylogenomic and phylodynamic analyses, we analysed 456 high-quality T pallidum genomes collected from clinical samples in Australia between Oct 19, 2005, and Dec 31, 2020, and contextualised this information with publicly available sequence data. We also performed detailed genomic characterisation of putative antimicrobial resistance determinants, in addition to correlating single-locus typing of the TP0548 allele with the T pallidum phylogeny.Phylogenomic analyses identified four major sublineages circulating in Australia and globally, two belonging to the SS14 lineage, and two belonging to the Nichols lineage. Australian sublineages were further delineated into twelve subgroups, with five of the six largest subgroups associated with men who have sex with men, and the sixth lineage was predominantly associated with heterosexual people. Most Australian T pallidum genomes (398 [87%] of 456) were genotypically macrolide resistant, and TP0548 typing correlated significantly with T pallidum genomic subgroups.These findings show that the current syphilis epidemic in Australia is driven by multiple lineages of T pallidum, rather than one distinct outbreak. Major subgroups of T pallidum in Australia have emerged within the past 30 years, are closely related to global lineages, and circulate across different sexual networks. In conjunction with improved testing and treatment, these data could better inform the control of syphilis in Australia.National Health and Medical Research Council, Australian Research Council.

Syphilis is a sexually transmitted infection caused by Treponema pallidum subspecies pallidum and may lead to severe complications. Recent years have seen striking increases in syphilis in many countries. Previous analyses have suggested one lineage of syphilis, SS14, may have expanded recently, indicating emergence of a single pandemic azithromycin-resistant cluster. Here we use direct sequencing of T. pallidum combined with phylogenomic analyses to show that both SS14- and Nichols-lineages are simultaneously circulating in clinically relevant populations in multiple countries. We correlate the appearance of genotypic macrolide resistance with multiple independently evolved SS14 sub-lineages and show that genotypically resistant and sensitive sub-lineages are spreading contemporaneously. These findings inform our understanding of the current syphilis epidemic by demonstrating how macrolide resistance evolves in Treponema subspecies and provide a warning on broader issues of antimicrobial resistance., Syphilis is caused by the bacterium Treponema pallidum subspecies pallidum (TPA), and incidence has risen recently in many countries. Here, Beale et al. provide whole-genome TPA sequences from 73 clinical samples and show how antimicrobial resistance emerged independently in circulating lineages.

BackgroundExtended spectrum beta-lactamases-producing (ESBL-P) and/or carbapenem-resistant (CR) Klebsiella pneumoniae have severely restricted available treatment options in healthcare settings in Vietnam. Understanding the diversity and transmission mechanisms of ESBL- and carbapenemase-encoding K. pneumoniae is important in both hospital and community settings for patient management.MethodsWe conducted a 6-month prospective cohort study of 69 Intensive care unit (ICU) patients from two hospitals in Hanoi, Vietnam. Longitudinally collected samples from patients and the ICU environment were cultured on selective media, and 357 K. pneumoniae colonies were whole genome sequenced. We performed phylogenetic analyses, and correlated phenotypic antimicrobial susceptibility testing with genotypic features of K. pneumoniae isolates. We constructed transmission networks of patient samples, relating ICU admission times and locations with genetic similarity of infecting K. pneumoniae.FindingsDespite being geographically and clinically separated, the two hospitals shared closely related strains carrying the same array of antimicrobial resistance genes. Many patients carried the same resistant K. pneumoniae clone from admission to discharge. 45.9% of total isolates carried both ESBL- and carbapenemase-encoding genes, with high minimum inhibitory concentrations. We found a novel co-occurrence of blaKPC-2 and blaNDM-1 in 46. 6% of samples from the globally successful ST15 lineage.InterpretationThese results highlight the high prevalence of ESBL-positive carbapenem-resistant K. pneumoniae in Vietnamese ICUs. Through studying K. pneumoniae ST15 in detail, we illustrated how important resistance genes are coalescing in stains carried broadly by patients entering the two hospitals directly or through referral.FundingThis study was supported by the Medical Research Council Newton Fund, United Kingdom (grant MR/N029399/1); the Ministry of Science and Technology, Vietnam (grant HNQT/SPÐP/04.16); This research was funded in whole by the Wellcome Trust (grant 206194). For the purpose of Open Access, the author has applied a CC BY public copyright licence to any Author Accepted Manuscript version arising from this submission.

Syphilis is an important sexually transmitted infection caused by the bacterium Treponema pallidum subspecies pallidum. The last two decades have seen syphilis incidence rise in many high-income countries, yet the evolutionary and epidemiological relationships that underpin this are poorly understood, as is the global T. pallidum population structure. We assembled a geographically and temporally diverse collection of clinical and laboratory samples comprising 726 T. pallidum genomes. We used detailed phylogenetic analysis and clustering to show that syphilis globally can be described by only two deeply branching lineages, Nichols and SS14. We show that both of these lineages can be found circulating concurrently in 12 of the 23 countries sampled. To provide further phylodynamic resolution we subdivided Treponema pallidum subspecies pallidum into 17 distinct sublineages. Importantly, like SS14, we provide evidence that two Nichols sublineages have expanded clonally across 9 countries contemporaneously with SS14. Moreover, pairwise genome analysis showed that recent isolates circulating in 14 different countries were genetically identical in their core genome to those from other countries, suggesting frequent exchange through international transmission pathways. This contrasts with the majority of samples collected prior to 1983, which are phylogenetically distinct from these more recently isolated sublineages. Bayesian temporal analysis provided evidence of a population bottleneck and decline occurring during the late 1990s, followed by a rapid population expansion a decade later. This was driven by the dominant T. pallidum sublineages circulating today, many of which are resistant to macrolides. Combined we show that the population of contemporary syphilis in high-income countries has undergone a recent and rapid global expansion.

Background: The SARS-CoV-2 variant B.1.1.7 was first identified in December, 2020, in England. We aimed to investigate whether increases in the proportion of infections with this variant are associated with differences in symptoms or disease course, reinfection rates, or transmissibility. Methods: We did an ecological study to examine the association between the regional proportion of infections with the SARS-CoV-2 B.1.1.7 variant and reported symptoms, disease course, rates of reinfection, and transmissibility. Data on types and duration of symptoms were obtained from longitudinal reports from users of the COVID Symptom Study app who reported a positive test for COVID-19 between Sept 28 and Dec 27, 2020 (during which the prevalence of B.1.1.7 increased most notably in parts of the UK). From this dataset, we also estimated the frequency of possible reinfection, defined as the presence of two reported positive tests separated by more than 90 days with a period of reporting no symptoms for more than 7 days before the second positive test. The proportion of SARS-CoV-2 infections with the B.1.1.7 variant across the UK was estimated with use of genomic data from the COVID-19 Genomics UK Consortium and data from Public Health England on spike-gene target failure (a non-specific indicator of the B.1.1.7 variant) in community cases in England. We used linear regression to examine the association between reported symptoms and proportion of B.1.1.7. We assessed the Spearman correlation between the proportion of B.1.1.7 cases and number of reinfections over time, and between the number of positive tests and reinfections. We estimated incidence for B.1.1.7 and previous variants, and compared the effective reproduction number, R t, for the two incidence estimates. Findings: From Sept 28 to Dec 27, 2020, positive COVID-19 tests were reported by 36 920 COVID Symptom Study app users whose region was known and who reported as healthy on app sign-up. We found no changes in reported symptoms or disease duration associated with B.1.1.7. For the same period, possible reinfections were identified in 249 (0·7% [95% CI 0·6–0·8]) of 36 509 app users who reported a positive swab test before Oct 1, 2020, but there was no evidence that the frequency of reinfections was higher for the B.1.1.7 variant than for pre-existing variants. Reinfection occurrences were more positively correlated with the overall regional rise in cases (Spearman correlation 0·56–0·69 for South East, London, and East of England) than with the regional increase in the proportion of infections with the B.1.1.7 variant (Spearman correlation 0·38–0·56 in the same regions), suggesting B.1.1.7 does not substantially alter the risk of reinfection. We found a multiplicative increase in the R t of B.1.1.7 by a factor of 1·35 (95% CI 1·02–1·69) relative to pre-existing variants. However, R t fell below 1 during regional and national lockdowns, even in regions with high proportions of infections with the B.1.1.7 variant. Interpretation: The lack of change in symptoms identified in this study indicates that existing testing and surveillance infrastructure do not need to change specifically for the B.1.1.7 variant. In addition, given that there was no apparent increase in the reinfection rate, vaccines are likely to remain effective against the B.1.1.7 variant. Funding: Zoe Global, Department of Health (UK), Wellcome Trust, Engineering and Physical Sciences Research Council (UK), National Institute for Health Research (UK), Medical Research Council (UK), Alzheimer's Society.

We confirmed endemicity and autochthonous transmission of yaws in Liberia after a population-based, community-led burden estimation (56,825 participants). Serologically confirmed yaws was rare and focal at population level (24 cases; 2.6 [95% CI 1.4-3.9] cases/10,000 population) with similar clinical epidemiology to other endemic countries in West Africa. Unsupervised classification of spatially referenced case finding data indicated that yaws was more likely to occur in hard-to-reach communities; healthcare-seeking was low among communities, and clinical awareness of yaws was low among healthcare workers. We recovered whole bacterial genomes from 12 cases and describe a monophyletic clade of Treponema pallidum subspecies pertenue, phylogenetically distinct from known TPE lineages, including those affecting neighboring nonhuman primate populations (Tai Forest, Cote d'Ivoire). Yaws is endemic in Liberia but exhibits low focal population prevalence with evidence of a historical genetic bottleneck and subsequent local expansion. Reporting gaps appear attributable to challenging epidemiology and low disease awareness.

Syphilis is an important sexually transmitted infection caused by the bacterium Treponema pallidum subspecies pallidum. The last two decades have seen syphilis incidence rise in many high-income countries, yet the evolutionary and epidemiological relationships that underpin this are poorly understood, as is the global T. pallidum population structure. We assembled a geographically and temporally diverse collection of clinical and laboratory samples comprising 726 T. pallidum genomes. We used detailed phylogenetic analysis and clustering to show that syphilis globally can be described by only two deeply branching lineages, Nichols and SS14. We show that both of these lineages can be found circulating concurrently in 12 of the 23 countries sampled. To provide further phylodynamic resolution we subdivided Treponema pallidum subspecies pallidum into 17 distinct sublineages. Importantly, like SS14, we provide evidence that two Nichols sublineages have expanded clonally across 9 countries contemporaneously with SS14. Moreover, pairwise genome analysis showed that recent isolates circulating in 14 different countries were genetically identical in their core genome to those from other countries, suggesting frequent exchange through international transmission pathways. This contrasts with the majority of samples collected prior to 1983, which are phylogenetically distinct from these more recently isolated sublineages. Bayesian temporal analysis provided evidence of a population bottleneck and decline occurring during the late 1990s, followed by a rapid population expansion a decade later. This was driven by the dominant T. pallidum sublineages circulating today, many of which are resistant to macrolides. Combined we show that the population of contemporary syphilis in high-income countries has undergone a recent and rapid global expansion.

Global dispersal and increasing frequency of the SARS-CoV-2 spike protein variant D614G are suggestive of a selective advantage but may also be due to a random founder effect. We investigate the hypothesis for positive selection of spike D614G in the United Kingdom using more than 25,000 whole genome SARS-CoV-2 sequences. Despite the availability of a large dataset, well represented by both spike 614 variants, not all approaches showed a conclusive signal of positive selection. Population genetic analysis indicates that 614G increases in frequency relative to 614D in a manner consistent with a selective advantage. We do not find any indication that patients infected with the spike 614G variant have higher COVID-19 mortality or clinical severity, but 614G is associated with higher viral load and younger age of patients. Significant differences in growth and size of 614G phylogenetic clusters indicate a need for continued study of this variant., Graphical Abstract, Highlights • Increasing frequency of SARS-CoV-2 D614G is consistent with a selective advantage • Phylodynamic analyses do not show significantly different growth of D614G clusters • There is no association of D614G replacement with greater severity of infection • The D614G replacement is associated with higher viral loads and younger patient age, Analysis of the spread and frequency of SARS-CoV-2 D614G in the United Kingdom suggests a selective advantage for this strain that is associated with higher viral loads in younger patients but not higher COVID-19 clinical severity or mortality.

Histoplasmosis is a serious infectious disease in humans caused by Histoplasma spp. (Onygenales), whose natural reservoirs are thought to be soil enriched with bird and bat guano. The true global burden of histoplasmosis is underestimated and frequently the pulmonary manifestations are misdiagnosed as tuberculosis. Molecular data on epidemiology of Histoplasma are still scarce, even though there is increasing recognition of histoplasmosis in recent years in areas distant from the traditional endemic regions in the Americas. We used multi-locus sequence data from protein coding loci (ADP-ribosylation factor, H antigen precursor, and delta-9 fatty acid desaturase), DNA barcoding (ITS1/2+5.8s), AFLP markers and mating type analysis to determine the genetic diversity, population structure and recognise the existence of different phylogenetic species among 436 isolates of Histoplasma obtained globally. Our study describes new phylogenetic species and the molecular characteristics of Histoplasma lineages causing outbreaks with a high number of severe outcomes in Northeast Brazil between 2011 and 2015. Genetic diversity levels provide evidence for recombination, common ancestry and clustering of Brazilian isolates at different geographic scales with the emergence of LAm C, a new genotype assigned to a separate population cluster in Northeast Brazil that exhibited low diversity indicative of isolation. The global survey revealed that the high genetic variability among Brazilian isolates along with the presence of divergent cryptic species and/or genotypes may support the hypothesis of Brazil being the center of dispersion of Histoplasma in South America, possibly with the contribution of migratory hosts such as birds and bats. Outside Brazil, the predominant species depends on the region. We confirm that histoplasmosis has significantly broadened its area of occurrence, an important feature of emerging pathogens. From a practical point of view, our data point to the emergence of histoplasmosis caused by a plethora of genotypes, and will enable epidemiological analysis focused on understanding the processes that lead to histoplasmosis. Further, the description of this diversity opens avenues for comparative genomic studies, which will allow progress toward a consensus taxonomy, improve understanding of the presence of hybrids in natural populations of medically relevant fungi, test reproductive barriers and to explore the significance of this variation.

The origin of syphilis has been hotly debated for decades. Ancient pathogen DNA may provide new evidence to redefine our understanding of this mystery, but is the mystery itself flawed in its assumptions?

Summary Background In a longitudinal study assessing the WHO strategy for yaws eradication using mass azithromycin treatment, we observed resurgence of yaws cases with dominance of a single JG8 sequence type and emergence of azithromycin-resistant Treponema pallidum subspecies pertenue (T p pertenue). Here, we analyse genomic changes in the bacterial population using samples collected during the study. Methods We did whole bacterial genome sequencing directly on DNA extracted from 37 skin lesion swabs collected from patients on Lihir Island, Papua New Guinea, between April 1, 2013, and Nov 1, 2016. We produced phylogenies and correlated these with spatiotemporal information to investigate the source of new cases and the emergence of five macrolide-resistant cases. We used deep amplicon sequencing of surveillance samples to assess the presence of minority macrolide-resistant populations. Findings We recovered 20 whole T p pertenue genomes, and phylogenetic analysis showed that the re-emerging JG8 sequence type was composed of three bacterial sublineages characterised by distinct spatiotemporal patterns. Of five patients with resistant T p pertenue, all epidemiologically linked, we recovered genomes from three and found no variants. Deep sequencing showed that before treatment, the index patient had fixed macrolide-sensitive T p pertenue, whereas the post-treatment sample had a fixed resistant genotype, as did three of four contact cases. Interpretation In this study, re-emergence of yaws cases was polyphyletic, indicating multiple epidemiological sources. However, given the genomic and epidemiological linkage of resistant cases and the rarity of resistance alleles in the general population, azithromycin resistance is likely to have evolved only once in this study, followed by onward dissemination. Funding Wellcome and Provincial Deputation of Barcelona.

This Genome Watch highlights how the SARS-CoV-2 pandemic laid the groundwork for continued use of real-time genomic epidemiology for public health responses in low-and-middle-income countries.

Paracoccidioidomycosis (PCM) is a life-threatening systemic fungal infection acquired after inhalation ofParacoccidioidespropagules from the environment. The main agents include members of theP. brasiliensiscomplex (phylogenetically-defined species S1, PS2, PS3, and PS4) andP. lutzii. DNA-sequencing of protein-coding loci (e.g.,GP43,ARF, andTUB1) is the reference method for recognizingParacoccidioidesspecies due to a lack of robust phenotypic markers. Thus, developing new molecular markers that are informative and cost-effective is key to providing quality information to explore genetic diversity withinParacoccidioides. We report using new amplified fragment length polymorphism (AFLP) markers and mating-type analysis for genotypingParacoccidioidesspecies. The bioinformatic analysis generated 144in silicoAFLP profiles, highlighting two discriminatory primer pairs combinations (#1 EcoRI-AC/MseI-CT and #2 EcoRI-AT/MseI-CT). The combinations #1 and #2 were usedin vitroto genotype 165Paracoccidioidesisolates recovered from across a vast area of South America. Considering the overall scored AFLP markersin vitro(67-87 fragments), the values of polymorphism information content (PIC= 0.3345-0.3456), marker index (MI= 0.0018), effective multiplex ratio (E= 44.6788-60.3818), resolving power (Rp= 22.3152-34.3152), discriminating power (D= 0.5183-0.5553), expected heterozygosity (H= 0.4247-0.4443), and mean heterozygosity (Havp = 0.00002-0.00004), demonstrated the utility of AFLP markers to speciateParacoccidioidesand to dissect both deep and fine-scale genetic structures. Analysis of molecular variance (AMOVA) revealed that the total genetic variance (65-66 %) was due to variability amongP. brasiliensiscomplex andP. lutzii(PhiPT = 0.651-0.658,P MAT1-1orMAT1-2idiomorphs were not significantly skewed (1:1 ratio) forP. brasiliensis s. str.(χ2= 1.025;P= 0.3113),P. venezuelensis(χ2= 0.692;P= 0.4054), andP. lutzii(χ2= 0.027;P= 0.8694), supporting random mating within each species. In contrast, skewed distributions were found forP. americana(χ2= 8.909;P= 0.0028) andP. restrepiensis(χ2= 4.571;P= 0.0325) with a preponderance ofMAT1-1. Geographical distributions confirmed thatP. americana,P. restrepiensis, andP. lutziiare more widespread than previously thought.P. brasiliensis s. str.is by far the most widely occurring lineage in Latin America countries, occurring in all regions of Brazil. Our new DNA fingerprint assay proved to be rapid, reproducible, and highly discriminatory, to give insights into the taxonomy, ecology, and epidemiology ofParacoccidioidesspecies, guiding disease-control strategies to mitigate PCM.

Background The WHO recommends treatment of trachoma with community-wide mass drug administration (MDA) with a single dose of azithromycin as part of the SAFE strategy (surgery, antibiotics, facial cleanliness and environmental improvement). In the Solomon Islands, this programme had demonstrable collateral benefit by reducing the prevalence of urogenital Chlamydia trachomatis(Ct) infections. We evaluated the impact of this treatment on population genetics of urogenital Ct. Methods Two vaginal swabs were collected from consecutive women attending antenatal clinics during cross-sectional surveys before and after MDA. For every swab positive for Ctinfection, DNA was extracted from the second swab, enriched and sequenced using paired-end sequencing. Whole-genome sequences were aligned against selected references. Diversity was assessed using genome-wide pairwise diversity and a high-resolution multi-locus sequence typing (hr-MLST-6) scheme. ARIBA software was used to test for evidence of antimicrobial resistance to macrolides. Results Whole-genome sequence data was obtained from 23/49 (47%) pre-MDA and 32/41 (78%) post-MDA Ct-positive samples. Most strains were serotype E and F, and tissue tropism genes were consistent with their urogenital nature. Genetic diversity of Ctwas lower by both pairwise and hr-MLST-6 diversity metrics in the post-MDA sample than the pre-MDA sample. There was no evidence of mutations known to confer resistance to macrolides in any of the samples collected. Conclusions Reduced diversity after MDA may represent selection pressure from mass antibiotic delivery. The absence of antimicrobial resistance is encouraging. The collateral impact (both positive and negative) of large-scale preventative chemotherapy programmes should be considered when deciding whether should be implemented.

Herpes Simplex Virus type 1 (HSV-1) chronically infects over 70 per cent of the global population. Clinical manifestations are largely restricted to recurrent epidermal vesicles. However, HSV-1 also leads to encephalitis, the infection of the brain parenchyma, with high associated rates of mortality and morbidity. In this study, we performed target enrichment followed by direct sequencing of HSV-1 genomes, using target enrichment methods on the cerebrospinal fluid (CSF) of clinical encephalitis patients and from skin swabs of epidermal vesicles on non-encephalopathic patients. Phylogenetic analysis revealed high inter-host diversity and little population structure. In contrast, samples from different lesions in the same patient clustered with similar patterns of allelic variants. Comparison of consensus genome sequences shows HSV-1 has been freely recombining, except for distinct islands of linkage disequilibrium (LD). This suggests functional constraints prevent recombination between certain genes, notably those encoding pairs of interacting proteins. Distinct LD patterns characterised subsets of viruses recovered from CSF and skin lesions, which may reflect different evolutionary constraints in different body compartments. Functions of genes under differential constraint related to immunity or tropism and provide new hypotheses on tissue-specific mechanisms of viral infection and latency.

Background Gonorrhoea is the third most commonly notified sexually transmitted infection (STI) in Singapore. In 2012, there were 1781 notifications, with an incidence rate of 33.53 per 100,000 population with more males than females being diagnosed. While most of the cases occur in people aged between 20 and 39 years of age, gonorrhoea is the most common STI among teenagers and among men who have sex with men (MSM) in Singapore. The aim of this study was to use whole genome sequencing to gain insights into the patterns of transmission that exist within and between different subpopulations in Singapore and internationally. Methods We sequenced 676 samples from 544 patients infected with gonorrhoea between 2014–2016. Sequencing reads from N. gonorrhoeae genomes were mapped to a common reference and recombination masked, followed by phylogenomic, Bayesian clustering and pairwise network analyses. We correlated genetic relatedness with detailed clinical parameters. Results N. gonorrhoeae circulating in Singapore is polyphyletic, and we defined 31 circulating sub-lineages. We detected distinct patterns of sexual behaviour associated with different genetic lineages: some lineages are strongly associated with MSM groups, whilst other lineages have increased rates of reported contact of commercial sex workers. It is likely that these associations reflect the underlying population within the transmission networks. We further correlate these genomic and behavioural subpopulations according to genetically determined antimicrobial resistance patterns. Conclusion The analysis shows distinct transmission clustering within Singapore groups based on sexual preference and commercial sex worker use. Through the use of multiple isolates from single individuals, we established expected within patient diversity levels based on pairwise sequence differences and used this to infer both putative transmission events and also possible bridging between distinct transmission networks. Further work is required to increase the prediction accuracy of the transmission networks and relate this to predicted microbial resistance patterns. Disclosure No significant relationships.

Novel adenovirus genotypes are associated with outbreaks of disease, such as acute gastroenteritis, renal disease, upper respiratory tract infection and keratoconjunctivitis. Here, we identify novel and variant adenovirus genotypes in children coinfected with enterotoxigenic Escherichia coli, in Bangladesh. Metagenomic sequencing of stool was performed and whole adenovirus genomes were extracted. A novel species D virus, designated genotype 90 (P33H27F67) was identified, and the partial genome of a putative recombinant species B virus was recovered. Furthermore, the enteric types HAdV-A61 and HAdV-A40 were found in stool specimens. Knowledge of the diversity of adenovirus genomes circulating worldwide, especially in low-income countries where the burden of disease is high, will be required to ensure that future vaccination strategies cover the diversity of adenovirus strains associated with disease.

Genome sequence data are of great value in describing evolutionary processes in viral populations. However, in such studies, the extent to which data accurately describes the viral population is a matter of importance. Multiple factors may influence the accuracy of a dataset, including the quantity and nature of the sample collected, and the subsequent steps in viral processing. To investigate this phenomenon, we sequenced replica datasets spanning a range of viruses, and in which the point at which samples were split was different in each case, from a dataset in which independent samples were collected from a single patient to another in which all processing steps up to sequencing were applied to a single sample before splitting the sample and sequencing each replicate. We conclude that neither a high read depth nor a high template number in a sample guarantee the precision of a dataset. Measures of consistency calculated from within a single biological sample may also be insufficient; distortion of the composition of a population by the experimental procedure or genuine within-host diversity between samples may each affect the results. Where it is possible, data from replicate samples should be collected to validate the consistency of short-read sequence data.

Cryptococcus neoformans var. grubii is the causative agent of cryptococcal meningitis, a significant source of mortality in immunocompromised individuals, typically human immunodeficiency virus/AIDS patients from developing countries. Despite the worldwide emergence of this ubiquitous infection, little is known about the global molecular epidemiology of this fungal pathogen. Here we sequence the genomes of 188 diverse isolates and characterize the major subdivisions, their relative diversity, and the level of genetic exchange between them. While most isolates of C. neoformans var. grubii belong to one of three major lineages (VNI, VNII, and VNB), some haploid isolates show hybrid ancestry including some that appear to have recently interbred, based on the detection of large blocks of each ancestry across each chromosome. Many isolates display evidence of aneuploidy, which was detected for all chromosomes. In diploid isolates of C. neoformans var. grubii (serotype AA) and of hybrids with C. neoformans var. neoformans (serotype AD) such aneuploidies have resulted in loss of heterozygosity, where a chromosomal region is represented by the genotype of only one parental isolate. Phylogenetic and population genomic analyses of isolates from Brazil reveal that the previously “African” VNB lineage occurs naturally in the South American environment. This suggests migration of the VNB lineage between Africa and South America prior to its diversification, supported by finding ancestral recombination events between isolates from different lineages and regions. The results provide evidence of substantial population structure, with all lineages showing multi-continental distributions; demonstrating the highly dispersive nature of this pathogen.

Cryptococcus neoformansvar.grubiiis the causative agent of cryptococcal meningitis, a significant source of mortality in immunocompromised individuals, typically HIV/AIDS patients from developing countries. Despite the worldwide emergence of this ubiquitous infection, little is known about the global molecular epidemiology of this fungal pathogen. Here we sequence the genomes of 188 diverse isolates and characterized the major subdivisions, their relative diversity and the level of genetic exchange between them. While most isolates ofC. neoformansvar.grubiibelong to one of three major lineages (VNI, VNII, and VNB), some haploid isolates show hybrid ancestry including some that appear to have recently interbred, based on the detection of large blocks of each ancestry across each chromosome. Many isolates display evidence of aneuploidy, which was detected for all chromosomes. In diploid isolates ofC. neoformansvar.grubii (serotype A/A) and of hybrids withC. neoformansvar.neoformans(serotype A/D) such aneuploidies have resulted in loss of heterozygosity, where a chromosomal region is represented by the genotype of only one parental isolate. Phylogenetic and population genomic analyses of isolates from Brazil revealed that the previously ‘African’ VNB lineage occurs naturally in the South American environment. This suggests migration of the VNB lineage between Africa and South America prior to its diversification, supported by finding ancestral recombination events between isolates from different lineages and regions. The results provide evidence of substantial population structure, with all lineages showing multi-continental distributions demonstrating the highly dispersive nature of this pathogen.Author SummaryCryptococcus neoformansvar.grubiiis a human fungal pathogen of immunocompromised individuals that has global clinical impact, causing half a million deaths per year. Substantial genetic substructure exists for this pathogen, with two lineages found globally (VNI, VNII) whereas a third has appeared confined to sub-Saharan Africa (VNB). Here, we utilized genome sequencing of a large set of global isolates to examine the genetic diversity, hybridization, and biogeography of these lineages. We found that while the three major lineages are well separated, recombination between the lineages has occurred, notably resulting in hybrid isolates with segmented ancestry across the genome. In addition, we showed that isolates from South America are placed within the VNB lineage, formerly thought to be confined to Africa, and that there is phylogenetic separation between these geographies that substantially expands the diversity of these lineages. Our findings provide a new framework for further studies of the dynamics of natural populations ofC. neoformansvar.grubii.

Recurrence of meningitis due toCryptococcus neoformansafter treatment causes substantial mortality in HIV/AIDS patients across sub-Saharan Africa. In order to determine whether recurrence occurred due to relapse of the original infecting isolate or reinfection with a different isolate weeks or months after initial treatment, we used whole-genome sequencing to assess the genetic basis of infection in 17 HIV-infected individuals with recurrent cryptococcal meningitis. Comparisons revealed a clonal relationship for 15 pairs of isolates recovered before and after recurrence showing relapse of the original infection. The two remaining pairs showed high levels of genetic heterogeneity; in one pair we found this to be a result of infection by mixed genotypes, whilst the second was a result of nonsense mutations in the gene encoding the DNA mismatch repair proteinsMSH2, MSH5andRAD5. These nonsense mutations led to a hypermutator state, leading to dramatically elevated rates of synonymous and non-synonymous substitutions. Hypermutator phenotypes owing to nonsense mutations in these genes have not previously been reported inCryptococcus neoformansand represent a novel pathway for rapid within-host adaptation and evolution of resistance to firstline antifungal drugs.

Emerging infections caused by fungi have become a widely recognized global phenomenon and are causing an increasing burden of disease. Genomic techniques are providing new insights into the structure of fungal populations, revealing hitherto undescribed fine‐scale adaptations to environments and hosts that govern their emergence as infections. Cryptococcal meningitis is a neglected tropical disease that is responsible for a large proportion of AIDS‐related deaths across Africa; however, the ecological determinants that underlie a patient's risk of infection remain largely unexplored. Here, we use genome sequencing and ecological genomics to decipher the evolutionary ecology of the aetiological agents of cryptococcal meningitis, Cryptococcus neoformans and Cryptococcus gattii, across the central African country of Zambia. We show that the occurrence of these two pathogens is differentially associated with biotic (macroecological) and abiotic (physical) factors across two key African ecoregions, Central Miombo woodlands and Zambezi Mopane woodlands. We show that speciation of Cryptococcus has resulted in adaptation to occupy different ecological niches, with C. neoformans found to occupy Zambezi Mopane woodlands and C. gattii primarily recovered from Central Miombo woodlands. Genome sequencing shows that C. neoformans causes 95% of human infections in this region, of which over three‐quarters belonged to the globalized lineage VNI. We show that VNI infections are largely associated with urbanized populations in Zambia. Conversely, the majority of C. neoformans isolates recovered in the environment belong to the genetically diverse African‐endemic lineage VNB, and we show hitherto unmapped levels of genomic diversity within this lineage. Our results reveal the complex evolutionary ecology that underpins the reservoirs of infection for this, and likely other, deadly pathogenic fungi.

Cryptococcal meningitis is a major cause of mortality throughout the developing world, yet little is known about the genetic markers underlying Cryptococcal virulence and patient outcome. We studied a cohort of 230 Cryptococcus neoformans (Cn) isolates from HIV-positive South African clinical trial patients with detailed clinical follow-up using multi-locus sequence typing and in vitro phenotypic virulence assays, correlating these data with clinical and fungal markers of disease in the patient. South African Cn displayed high levels of genetic diversity and locus variability compared to globally distributed types, and we identified 50 sequence types grouped within the main molecular types VNI, VNII and VNB, with 72% of isolates typed into one of seven 'high frequency' sequence types. Spatial analysis of patients’ cryptococcal genotype was not shown to be clustered geographically, which might argue against recent local acquisition and in favour of reactivation of latent infection. Through comparison of MLST genotyping data with clinical parameters, we found a relationship between genetic lineage and clinical outcome, with patients infected with the VNB lineage having significantly worse survival (n=8, HR 3.35, CI 1.51-7.20, p=0.003), and this was maintained even after adjustment for known prognostic indicators and treatment regimen. Comparison of fungal genotype with in vitro phenotype (phagocytosis, laccase activity and CSF survival) performed on a subset of 89 isolates revealed evidence of lineage-associated virulence phenotype, with the VNII lineage displaying increased laccase activity (p=0.001) and ex vivo CSF survival (p=0.0001). These findings show that Cryptococcus neoformans is a phenotypically heterogeneous pathogen, and that lineage plays an important role in cryptococcal virulence during human infection. Furthermore, a detailed understanding of the genetic diversity in Southern Africa will support further investigation into how genetic diversity is structured across African environments, allowing assessment of the risks different ecotypes pose to infection., Author Summary Cryptococcus neoformans (Cn) is a yeast that commonly causes meningitis in HIV infected individuals in Africa, where it may account for up to 500,000 deaths every year. In this highly translational and multidisciplinary study, we used genetic analysis techniques to show that Cryptococcus found in Southern Africa represents a hotspot of genetic diversity. We combined this data with the results of microbiological techniques that assess the natural virulence traits that the yeast uses to survive and infect humans to further show that genetic diversity is associated with differences in cryptococcal phenotype. Finally, we analysed detailed clinical data on patients to investigate the clinical effects of infection with different lineages, and showed that one genetic lineage (VNB) is significantly associated with worse survival. Whilst much of our prior knowledge regarding the genetic basis of virulence is derived from studies on laboratory-adapted cryptococcal strains, our findings from this large and comprehensive MLST genotyping study of clinical isolates—linking genotype, phenotype, clinical presentation and outcome—provide direct insights into the contribution of pathogen lineage to virulence in human cryptococcal meningitis.

This work was made possible by funding from the Wellcome Trust (WT088148MF), the Medical Research Council (MR/J008176/1), the NIHR Surgical Reconstruction and Microbiology Research Centre and the Lister Institute for Preventive Medicine (to R.C. May), and a Wellcome Trust Intermediate fellowship (089966, to T. Bicanic). The C. neoformans isolates were collected within clinical trials funded by the British Infection Society (fellowship to T. Bicanic), the Wellcome Trust (research training fellowships WT069991, to A.E. Brouwer and WT081794, to J.N. Jarvis), and the Medical Research Council, United Kingdom (76201). Background. Cryptococcal meningitis (CM) is a leading cause of HIV-associated mortality globally. High fungal burden in cerebrospinal fluid (CSF) at diagnosis and poor fungal clearance during treatment are recognized adverse prognostic markers; however, the underlying pathogenic factors that drive these clinical manifestations are incompletely understood. We profiled a large set of clinical isolates for established cryptococcal virulence traits to evaluate the contribution of C. neoformans phenotypic diversity to clinical presentation and outcome in human cryptococcosis. Methods. Sixty-five C. neoformans isolates from clinical trial patients with matched clinical data were assayed in vitro to determine murine macrophage uptake, intracellular proliferation rate (IPR), capsule induction, and laccase activity. Analysis of the correlation between prognostic clinical and host immune parameters and fungal phenotypes was performed using Spearman's r, while the fungal-dependent impact on long-term survival was determined by Cox regression analysis. Results. High levels of fungal uptake by macrophages in vitro, but not the IPR, were associated with CSF fungal burden (r = 0.38, P = 0.002) and long-term patient survival (hazard ratio [HR] 2.6, 95% CI 1.2-5.5, P = 0.012). High-uptake strains were hypocapsular (r = -0.28, P = 0.05) and exhibited enhanced laccase activity (r = 0.36, P = 0.003). Fungal isolates with greater laccase activity exhibited heightened survival ex vivo in purified CSF (r = 0.49, P P = 0.05) and exhibited enhanced laccase activity (r = 0.36, P = 0.003). Fungal isolates with greater laccase activity exhibited heightened survival ex vivo in purified CSF (r = 0.49, P < 0.0001) and resistance to clearance following patient antifungal treatment (r = 0.39, P = 0.003). Conclusion. These findings underscore the contribution of cryptococcal-phagocyte interactions and laccase-dependent melanin pathways to human clinical presentation and outcome. Furthermore, characterization of fungal-specific pathways that drive clinical manifestation provide potential targets for the development of therapeutics and the management of CM. Publisher PDF

Background and Objective The risk of hepatitis E virus (HEV) to blood safety remains unknown in England. Reports of persistent HEV infection with serious disease sequelae indicate that transfusion transmitted HEV is not a trivial disease in immunosuppressed patients. Materials and Methods Samples from unselected blood donors and donors with a history of jaundice were tested for HEV antibody and RNA. Results Overall, 10% of the donor sera were anti-HEV IgG reactive. Four of the donor samples were anti-HEV IgM reactive but HEV RNA negative. Conclusion There is evidence of probable recent HEV infections in donors with a predicted attack rate of 2·8%.

A rapid and global emergence of azole resistance has been observed in the pathogenic fungus Aspergillus fumigatus over the past decade. The dominant resistance mechanism appears to be of environmental origin and involves mutations in the cyp51A gene, which encodes a protein targeted by triazole antifungal drugs. Whole-genome sequencing (WGS) was performed for high-resolution single-nucleotide polymorphism (SNP) analysis of 24 A. fumigatus isolates, including azole-resistant and susceptible clinical and environmental strains obtained from India, the Netherlands, and the United Kingdom, in order to assess the utility of WGS for characterizing the alleles causing resistance. WGS analysis confirmed that TR34/L98H (a mutation comprising a tandem repeat [TR] of 34 bases in the promoter of the cyp51A gene and a leucine-to-histidine change at codon 98) is the sole mechanism of azole resistance among the isolates tested in this panel of isolates. We used population genomic analysis and showed that A. fumigatus was panmictic, with as much genetic diversity found within a country as is found between continents. A striking exception to this was shown in India, where isolates are highly related despite being isolated from both clinical and environmental sources across >1,000 km; this broad occurrence suggests a recent selective sweep of a highly fit genotype that is associated with the TR34/L98H allele. We found that these sequenced isolates are all recombining, showing that azole-resistant alleles are segregating into diverse genetic backgrounds. Our analysis delineates the fundamental population genetic parameters that are needed to enable the use of genome-wide association studies to identify the contribution of SNP diversity to the generation and spread of azole resistance in this medically important fungus., IMPORTANCE Resistance to azoles in the ubiquitous ascomycete fungus A. fumigatus was first reported from clinical isolates collected in the United States during the late 1980s. Over the last decade, an increasing number of A. fumigatus isolates from the clinic and from nature have been found to show resistance to azoles, suggesting that resistance is emerging through selection by the widespread usage of agricultural azole antifungal compounds. Aspergillosis is an emerging clinical problem, with high rates of treatment failures necessitating the development of new techniques for surveillance and for determining the genome-wide basis of azole resistance in A. fumigatus.