Genomic diversity and antimicrobial resistance in clinical Klebsiella pneumoniae isolates from tertiary hospitals in Southern Ghana.

OBJECTIVES: Comprehensive data on the genomic epidemiology of hospital-associated Klebsiella pneumoniae in Ghana are scarce. This study investigated the genomic diversity, antimicrobial resistance patterns, and clonal relationships of 103 clinical K. pneumoniae isolates from five tertiary hospitals in Southern Ghana-predominantly from paediatric patients aged under 5 years (67/103; 65%), with the majority collected from urine (32/103; 31%) and blood (25/103; 24%) cultures. METHODS: We generated hybrid Nanopore-Illumina assemblies and employed Pathogenwatch for genotyping via Kaptive [capsular (K) locus and lipopolysaccharide (O) antigens] and Kleborate (antimicrobial resistance and hypervirulence) and determined clonal relationships using core-genome MLST (cgMLST). RESULTS: Of 44 distinct STs detected, ST133 was the most common, comprising 23% of isolates (n = 23/103). KL116 (28/103; 27%) and O1 (66/103; 64%) were the most prevalent K-locus and O-antigen types. Single-linkage clustering highlighted the global spread of MDR clones such as ST15, ST307, ST17, ST11, ST101 and ST48, with minimal allele differences (1-5) from publicly available genomes worldwide. Conversely, 17 isolates constituted novel clonal groups and lacked close relatives among publicly available genomes, displaying unique genetic diversity within our study population. A significant proportion of isolates (88/103; 85%) carried resistance genes for ≥3 antibiotic classes, with the blaCTX-M-15 gene present in 78% (n = 80/103). Carbapenem resistance, predominantly due to blaOXA-181 and blaNDM-1 genes, was found in 10% (n = 10/103) of the isolates. CONCLUSIONS: Our findings reveal a complex genomic landscape of K. pneumoniae in Southern Ghana, underscoring the critical need for ongoing genomic surveillance to manage the substantial burden of antimicrobial resistance.

AlbaTraDIS: Comparative analysis of large datasets from parallel transposon mutagenesis experiments.

Bacteria need to survive in a wide range of environments. Currently, there is an incomplete understanding of the genetic basis for mechanisms underpinning survival in stressful conditions, such as the presence of anti-microbials. Transposon directed insertion-site sequencing (TraDIS) is a powerful tool to identify genes and networks which are involved in survival and fitness under a given condition by simultaneously assaying the fitness of millions of mutants, thereby relating genotype to phenotype and contributing to an understanding of bacterial cell biology. A recent refinement of this approach allows the roles of essential genes in conditional stress survival to be inferred by altering their expression. These advancements combined with the rapidly falling costs of sequencing now allows comparisons between multiple experiments to identify commonalities in stress responses to different conditions. This capacity however poses a new challenge for analysis of multiple data sets in conjunction. To address this analysis need, we have developed 'AlbaTraDIS'; a software application for rapid large-scale comparative analysis of TraDIS experiments that predicts the impact of transposon insertions on nearby genes. AlbaTraDIS can identify genes which are up or down regulated, or inactivated, between multiple conditions, producing a filtered list of genes for further experimental validation as well as several accompanying data visualisations. We demonstrate the utility of our new approach by applying it to identify genes used by Escherichia coli to survive in a wide range of different concentrations of the biocide Triclosan. AlbaTraDIS identified all well characterised Triclosan resistance genes, including the primary target, fabI. A number of new loci were also implicated in Triclosan resistance and the predicted phenotypes for a selection of these were validated experimentally with results being consistent with predictions. AlbaTraDIS provides a simple and rapid method to analyse multiple transposon mutagenesis data sets allowing this technology to be used at large scale. To our knowledge this is the only tool currently available that can perform these tasks. AlbaTraDIS is written in Python 3 and is available under the open source licence GNU GPL 3 from https://github.com/quadram-institute-bioscience/albatradis.

A dengue outbreak on a floating village at Cat Ba Island in Vietnam.

BACKGROUND: A dengue outbreak in an ecotourism destination spot in Vietnam, from September to November 2013, impacted a floating village of fishermen on the coastal island of Cat Ba. The outbreak raises questions about how tourism may impact disease spread in rural areas. METHODS: Epidemiological data were obtained from the Hai Phong Preventive Medical Center (PMC), including case histories and residential location from all notified dengue cases from this outbreak. All household addresses were geo-located. Knox test, a spatio-temporal analysis that enables inference dengue clustering constrained by space and time, was performed on the geocoded locations. From the plasma available from two patients, positive for Dengue serotype 3 virus (DENV3), the Envelope (E) gene was sequenced, and their genetic relationships compared to other E sequences in the region. RESULTS: Of 192 dengue cases, the odds ratio of contracting dengue infections for people living in the floating villages compared to those living on the island was 4.9 (95 % CI: 3.6-6.7). The space-time analyses on 111 geocoded dengue residences found the risk of dengue infection to be the highest within 4 days and a radius of 20 m of a given case. Of the total of ten detected clusters with an excess risk greater than 2, the cluster with the highest number of cases was in the floating village area (24 patients for a total duration of 31 days). Phylogenetic analysis revealed a high homology of the two DENV3 strains (genotype III) from Cat Ba with DENV3 viruses circulating in Hanoi in the same year (99.1 %). CONCLUSIONS: Our study showed that dengue transmission is unlikely to be sustained on Cat Ba Island and that the 2013 epidemic likely originated through introduction of viruses from the mainland, potentially Hanoi. These findings suggest that prevention efforts should be focused on mainland rather than on the island.

Capturing clinically relevant Campylobacter attributes through direct whole genome sequencing of stool.

Campylobacter is the leading bacterial cause of infectious intestinal disease, but the pathogen typically accounts for a very small proportion of the overall stool microbiome in each patient. Diagnosis is even more difficult due to the fastidious nature of Campylobacter in the laboratory setting. This has, in part, driven a change in recent years, from culture-based to rapid PCR-based diagnostic assays which have improved diagnostic detection, whilst creating a knowledge gap in our clinical and epidemiological understanding of Campylobacter genotypes - no isolates to sequence. In this study, direct metagenomic sequencing approaches were used to assess the possibility of replacing genome sequences with metagenome sequences; metagenomic sequencing outputs were used to describe clinically relevant attributes of Campylobacter genotypes. A total of 37 diarrhoeal stool samples with Campylobacter and five samples with an unknown pathogen result were collected and processed with and without filtration, DNA was extracted, and metagenomes were sequenced by short-read sequencing. Culture-based methods were used to validate Campylobacter metagenome-derived genome (MDG) results. Sequence output metrics were assessed for Campylobacter genome quality and accuracy of characterization. Of the 42 samples passing quality checks for analysis, identification of Campylobacter to the genus and species level was dependent on Campylobacter genome read count, coverage and genome completeness. A total of 65% (24/37) of samples were reliably identified to the genus level through Campylobacter MDG, 73% (27/37) by culture and 97% (36/37) by qPCR. The Campylobacter genomes with a genome completeness of over 60% (n=21) were all accurately identified at the species level (100%). Of those, 72% (15/21) were identified to sequence types (STs), and 95% (20/21) accurately identified antimicrobial resistance (AMR) gene determinants. Filtration of stool samples enhanced Campylobacter MDG recovery and genome quality metrics compared to the corresponding unfiltered samples, which improved the identification of STs and AMR profiles. The phylogenetic analysis in this study demonstrated the clustering of the metagenome-derived with culture-derived genomes and revealed the reliability of genomes from direct stool sequencing. Furthermore, Campylobacter genome spiking percentages ranging from 0 to 2% total metagenome abundance in the ONT MinION sequencer, configured to adaptive sequencing, exhibited better assembly quality and accurate identification of STs, particularly in the analysis of metagenomes containing 2 and 1% of Campylobacter jejuni genomes. Direct sequencing of Campylobacter from stool samples provides clinically relevant and epidemiologically important genomic information without the reliance on cultured genomes.

Development and proof-of-concept demonstration of a clinical metagenomics method for the rapid detection of bloodstream infection.

BACKGROUND: The timely and accurate diagnosis of bloodstream infection (BSI) is critical for patient management. With longstanding challenges for routine blood culture, metagenomics is a promising approach to rapidly provide sequence-based detection and characterisation of bloodborne bacteria. Long-read sequencing technologies have successfully supported the use of clinical metagenomics for syndromes such as respiratory illness, and modified approaches may address two requisite factors for metagenomics to be used as a BSI diagnostic: depletion of the high level of host DNA to then detect the low abundance of microbes in blood. METHODS: Blood samples from healthy donors were spiked with different concentrations of four prevalent causative species of BSI. All samples were then subjected to a modified saponin-based host DNA depletion protocol and optimised DNA extraction, whole genome amplification and debranching steps in preparation for sequencing, followed by bioinformatical analyses. Two related variants of the protocol are presented: 1mL of blood processed without bacterial enrichment, and 5mL of blood processed following a rapid bacterial enrichment protocol-SepsiPURE. RESULTS: After first identifying that a large proportion of host mitochondrial DNA remained, the host depletion process was optimised by increasing saponin concentration to 3% and scaling the reaction to allow more sample volume. Compared to non-depleted controls, the 3% saponin-based depletion protocol reduced the presence of host chromosomal and mitochondrial DNA < 106 and < 103 fold respectively. When the modified depletion method was further combined with a rapid bacterial enrichment method (SepsiPURE; with 5mL blood samples) the depletion of mitochondrial DNA improved by a further > 10X while also increasing detectable bacteria by > 10X. Parameters during DNA extraction, whole genome amplification and long-read sequencing were also adjusted, and subsequently amplicons were detected for each input bacterial species at each of the spiked concentrations, ranging from 50-100 colony forming units (CFU)/mL to 1-5 CFU/mL. CONCLUSION: In this proof-of-concept study, four prevalent BSI causative species were detected in under 12 h to species level (with antimicrobial resistance determinants) at concentrations relevant to clinical blood samples. The use of a rapid and precise metagenomic protocols has the potential to advance the diagnosis of BSI.

A hybrid and poly-polish workflow for the complete and accurate assembly of phage genomes: a case study of ten przondoviruses.

Bacteriophages (phages) within the genus Przondovirus are T7-like podoviruses belonging to the subfamily Studiervirinae, within the family Autographiviridae, and have a highly conserved genome organisation. The genomes of these phages range from 37 to 42 kb in size, encode 50-60 genes and are characterised by the presence of direct terminal repeats (DTRs) flanking the linear chromosome. These DTRs are often deleted during short-read-only and hybrid assemblies. Moreover, long-read-only assemblies are often littered with sequencing and/or assembly errors and require additional curation. Here, we present the isolation and characterisation of ten novel przondoviruses targeting Klebsiella spp. We describe HYPPA, a HYbrid and Poly-polish Phage Assembly workflow, which utilises long-read assemblies in combination with short-read sequencing to resolve phage DTRs and correcting errors, negating the need for laborious primer walking and Sanger sequencing validation. Our assembly workflow utilised Oxford Nanopore Technologies for long-read sequencing for its accessibility, making it the more relevant long-read sequencing technology at this time, and Illumina DNA Prep for short-read sequencing, representing the most commonly used technologies globally. Our data demonstrate the importance of careful curation of phage assemblies before publication, and prior to using them for comparative genomics.

The use of representative community samples to assess SARS-CoV-2 lineage competition: Alpha outcompetes Beta and wild-type in England from January to March 2021.

Genomic surveillance for SARS-CoV-2 lineages informs our understanding of possible future changes in transmissibility and vaccine efficacy and will be a high priority for public health for the foreseeable future. However, small changes in the frequency of one lineage over another are often difficult to interpret because surveillance samples are obtained using a variety of methods all of which are known to contain biases. As a case study, using an approach which is largely free of biases, we here describe lineage dynamics and phylogenetic relationships of the Alpha and Beta variant in England during the first 3 months of 2021 using sequences obtained from a random community sample who provided a throat and nose swab for rt-PCR as part of the REal-time Assessment of Community Transmission-1 (REACT-1) study. Overall, diversity decreased during the first quarter of 2021, with the Alpha variant (first identified in Kent) becoming predominant, driven by a reproduction number 0.3 higher than for the prior wild-type. During January, positive samples were more likely to be Alpha in those aged 18 to 54 years old. Although individuals infected with the Alpha variant were no more likely to report one or more classic COVID-19 symptoms compared to those infected with wild-type, they were more likely to be antibody-positive 6 weeks after infection. Further, viral load was higher in those infected with the Alpha variant as measured by cycle threshold (Ct) values. The presence of infections with non-imported Beta variant (first identified in South Africa) during January, but not during February or March, suggests initial establishment in the community followed by fade-out. However, this occurred during a period of stringent social distancing. These results highlight how sequence data from representative community surveys such as REACT-1 can augment routine genomic surveillance during periods of lineage diversity.

Characterization of Genetic Elements Carrying mcr-1 Gene in Escherichia coli from the Community and Hospital Settings in Vietnam.

Colistin is widely used in agriculture and aquaculture as prophylaxis, particularly in Asia. Recently, mcr-1 and other mobilizable genes conferring colistin resistance have spread globally in community and hospital populations. Characterizing mcr-1 mobile genetic elements and host genetic background is important to understand the transmission of this resistance mechanism. We conducted whole-genome sequencing of 94 mcr-1-positive Escherichia coli isolates (Mcr1-Ec isolates) from human and animal feces, food, and water in a community cohort (N = 87) and from clinical specimens from a referral hospital (N = 7) in northern Vietnam. mcr-1 was plasmid-borne in 71 and chromosomally carried in 25 (2 isolates contain one copy on chromosome and one copy on a plasmid) of 94 E. coli isolates from the community and hospital settings. All seven clinical isolates carried mcr-1 on plasmids. Replicon types of mcr-1-carrying plasmids included IncI2, IncP, IncX4, and IncFIA single replicons and combinations of IncHI2, IncN, and IncX1 multireplicons. Alignment of a long-read sequence of an IncI2 plasmid from animal feces with short-read sequences of IncI2 plasmids from a healthy human, water, and hospitalized patients showed highly similar structures (query cover from 90% to 98%, overall identity of >81%). We detected the potential existence of multireplicon plasmids harboring mcr-1 regardless of sample setting, confirming 10/71 with long-read sequencing. An intact/conserved Tn6330 transposon sequence or its genetic context variants were found in 6/25 Mcr1-Ec isolates with chromosomally carried mcr-1. The dissemination of mcr-1 is facilitated by a high diversity of plasmid replicon types and a high prevalence of the chromosomal Tn6330 transposon. IMPORTANCE The article presented advances our understanding of genetic elements carrying mcr-1 in Escherichia coli in both community and hospital settings. We provide evidence to suggest that diverse plasmid types, including multireplicon plasmids, have facilitated the successful transmission of mcr-1 in different reservoirs. The widespread use of colistin in agriculture, where a high diversity of bacteria are exposed, has allowed the selection and evolution of various transmission mechanisms that will make it a challenge to get rid of. Colocalization of mcr-1 and other antibiotic resistance genes (ARGs) on multireplicon plasmids adds another layer of complexity to the rapid dissemination of mcr-1 genes among community and hospital bacterial populations and to the slow pandemic of antimicrobial resistance (AMR) in general.

Hybrid metagenome assemblies link carbohydrate structure with function in the human gut microbiome.

Complex carbohydrates that escape small intestinal digestion, are broken down in the large intestine by enzymes encoded by the gut microbiome. This is a symbiotic relationship between microbes and host, resulting in metabolic products that influence host health and are exploited by other microbes. However, the role of carbohydrate structure in directing microbiota community composition and the succession of carbohydrate-degrading microbes, is not fully understood. In this study we evaluate species-level compositional variation within a single microbiome in response to six structurally distinct carbohydrates in a controlled model gut using hybrid metagenome assemblies. We identified 509 high-quality metagenome-assembled genomes (MAGs) belonging to ten bacterial classes and 28 bacterial families. Bacterial species identified as carrying genes encoding starch binding modules increased in abundance in response to starches. The use of hybrid metagenomics has allowed identification of several uncultured species with the functional potential to degrade starch substrates for future study.

Genomic assessment of quarantine measures to prevent SARS-CoV-2 importation and transmission.

Mitigation of SARS-CoV-2 transmission from international travel is a priority. We evaluated the effectiveness of travellers being required to quarantine for 14-days on return to England in Summer 2020. We identified 4,207 travel-related SARS-CoV-2 cases and their contacts, and identified 827 associated SARS-CoV-2 genomes. Overall, quarantine was associated with a lower rate of contacts, and the impact of quarantine was greatest in the 16-20 age-group. 186 SARS-CoV-2 genomes were sufficiently unique to identify travel-related clusters. Fewer genomically-linked cases were observed for index cases who returned from countries with quarantine requirement compared to countries with no quarantine requirement. This difference was explained by fewer importation events per identified genome for these cases, as opposed to fewer onward contacts per case. Overall, our study demonstrates that a 14-day quarantine period reduces, but does not completely eliminate, the onward transmission of imported cases, mainly by dissuading travel to countries with a quarantine requirement.

Genomic epidemiological analysis of mcr-1-harboring Escherichia coli collected from livestock settings in Vietnam.

Livestock has been implicated as a reservoir for antimicrobial resistance (AMR) genes that can spread to humans when antimicrobials are used in animals for food production to treat clinical diseases and prevent and control common disease events. In Vietnam, mcr-1-harboring Escherichia coli (MCRPEC) strains have been isolated from humans, animals (chickens, pigs, and dogs) feces, flies, foods, and the environment (rainwater, well water, and irrigation water) in communities and from clinical specimens in hospitals. The relationship between levels of AMR in livestock and its occurrence in humans is complex and is driven by many factors. We conducted whole genome sequencing of MCRPEC to analyze the molecular epidemiological characteristics, history, and relatedness of 50 isolates obtained in 2019 from different reservoirs in farms and markets in Ha Nam province, Vietnam. 34 sequence types (STs) with 3 new STs were identified in multilocus sequence typing analysis: ST12945 and ST12946 from chicken feces, and ST12947 from flies. The AMR phenotypes of 50 MCRPEC isolates were as follows: ampicillin (100%, 50/50), cefotaxime (10%, 5/50), gentamicin (60%, 30/50), amikacin (8%, 4/50), meropenem (6%, 3/50), ceftazidime (18%, 9/50), colistin (24%, 12/50) and ciprofloxacin (80%, 40/50). All 50 MCRPEC isolates were identified as MDR. 100% (50/50) isolates carried AMR genes, ranging from 5 to 22 genes. The most prevalent plasmid replicon types carrying mcr-1 were IncP-1 (17/37, 45.9%), IncX4 (7/37, 18.9%), and IncHI2/IncHI2A (6/37, 16.2%). These data suggest that the epidemiology of the mcr-1 gene is mostly determined by plasmid spreading instead of clonal dissemination of MCRPE strains. The co-occurrence of several STs such as ST10, ST48, ST155, ST206, ST2705 in various sample types, joined to the higher prevalence of a few types of Inc plasmids, confirms the dissemination of the mcr-1 carrying plasmids in E. coli clones established in livestock. 5 over 8 STs identified in flies (ST206, ST2705, ST155, ST10, and ST48) suggested the fly contribution in the transmission of AMR bacteria in environments. These popular STs also occur in human samples and 100% of the human samples were positive for the mcr-1 gene.

Replacement of the Alpha variant of SARS-CoV-2 by the Delta variant in Lebanon between April and June 2021.

The COVID-19 pandemic continues to expand globally, with case numbers rising in many areas of the world, including the Eastern Mediterranean Region. Lebanon experienced its largest wave of COVID-19 infections from January to April 2021. Limited genomic surveillance was undertaken, with just 26 SARS-CoV-2 genomes available for this period, nine of which were from travellers from Lebanon detected by other countries. Additional genome sequencing is thus needed to allow surveillance of variants in circulation. In total, 905 SARS-CoV-2 genomes were sequenced using the ARTIC protocol. The genomes were derived from SARS-CoV-2-positive samples, selected retrospectively from the sentinel COVID-19 surveillance network, to capture diversity of location, sampling time, sex, nationality and age. Although 16 PANGO lineages were circulating in Lebanon in January 2021, by February there were just four, with the Alpha variant accounting for 97 % of samples. In the following 2 months, all samples contained the Alpha variant. However, this had changed dramatically by June and July 2021, when all samples belonged to the Delta variant. This study documents a ten-fold increase in the number of SARS-CoV-2 genomes available from Lebanon. The Alpha variant, first detected in the UK, rapidly swept through Lebanon, causing the country's largest wave to date, which peaked in January 2021. The Alpha variant was introduced to Lebanon multiple times despite travel restrictions, but the source of these introductions remains uncertain. The Delta variant was detected in Gambia in travellers from Lebanon in mid-May, suggesting community transmission in Lebanon several weeks before this variant was detected in the country. Prospective sequencing in June/July 2021 showed that the Delta variant had completely replaced the Alpha variant in under 6 weeks.

SARS-CoV-2 variants of concern dominate in Lahore, Pakistan in April 2021.

The SARS-CoV-2 pandemic continues to expand globally, with case numbers rising in many areas of the world, including the Indian sub-continent. Pakistan has one of the world's largest populations, of over 200 million people and is experiencing a severe third wave of infections caused by SARS-CoV-2 that began in March 2021. In Pakistan, during the third wave until now only 12 SARS-CoV-2 genomes have been collected and among these nine are from Islamabad. This highlights the need for more genome sequencing to allow surveillance of variants in circulation. In fact, more genomes are available among travellers with a travel history from Pakistan, than from within the country itself. We thus aimed to provide a snapshot assessment of circulating lineages in Lahore and surrounding areas with a combined population of 11.1 million. Within a week of April 2021, 102 samples were sequenced. The samples were randomly collected from two hospitals with a diagnostic PCR cutoff value of less than 25 cycles. Analysis of the lineages shows that the Alpha variant of concern (first identified in the UK) dominates, accounting for 97.9 % (97/99) of cases, with the Beta variant of concern (first identified in South Africa) accounting for 2.0 % (2/99) of cases. No other lineages were observed. In depth analysis of the Alpha lineages indicated multiple separate introductions and subsequent establishment within the region. Eight samples were identical to genomes observed in Europe (seven UK, one Switzerland), indicating recent transmission. Genomes of other samples show evidence that these have evolved, indicating sustained transmission over a period of time either within Pakistan or other countries with low-density genome sequencing. Vaccines remain effective against Alpha, however, the low level of Beta against which some vaccines are less effective demonstrates the requirement for continued prospective genomic surveillance.

Resistance mechanisms and genetic relatedness among carbapenem-resistant Pseudomonas aeruginosa isolates from three major hospitals in Hanoi, Vietnam (2011-15).

BACKGROUND: MDR bacteria including carbapenem-resistant Pseudomonas aeruginosa are recognized as an important cause of hospital-acquired infections worldwide. This investigation seeks to determine the molecular characterization and antibiotic resistance genes associated with carbapenem-resistant P. aeruginosa. METHODS: We conducted WGS and phylogenetic analysis of 72 carbapenem-resistant P. aeruginosa isolated from hospital-acquired infection patients from August 2011 to March 2015 in three major hospitals in Hanoi, Vietnam. RESULTS: We identified three variants of IMP gene, among which bla IMP-15 was the most frequent (n = 34) in comparison to bla IMP-26 (n = 2) and bla IMP-51 (n = 12). We observed two isolates with imipenem MIC >128 mg/L that co-harboured bla IMP-15 and bla DIM-1 genes and seven isolates (imipenem MIC > 128 mg/L) with a bla KPC-1 gene from the same hospital. MLST data shows that these 72 isolates belong to 18 STs and phylogenetic tree analysis has divided these isolates into nine groups. CONCLUSIONS: Our results provide evidence that not only bla IMP-26 but other IMP variants such as bla IMP-15 and bla IMP-51 genes and several STs (ST235, ST244, ST277, ST310, ST773 and ST3151) have been disseminating in healthcare settings in Vietnam. In addition, we report the emergence of two isolates belonging to ST1240 and ST3340 that harboured two important carbapenemase genes (bla IMP-15 and bla DIM-1) and seven isolates belonging to ST3151 of P. aeruginosa that carried the bla KPC-1 gene in Vietnam, which could potentially cause serious restricted availability of treatment options in healthcare settings.

Genomic epidemiology and the role of international and regional travel in the SARS-CoV-2 epidemic in Zimbabwe: a retrospective study of routinely collected surveillance data.

BACKGROUND: Advances in SARS-CoV-2 sequencing have enabled identification of new variants, tracking of its evolution, and monitoring of its spread. We aimed to use whole genome sequencing to describe the molecular epidemiology of the SARS-CoV-2 outbreak and to inform the implementation of effective public health interventions for control in Zimbabwe. METHODS: We performed a retrospective study of nasopharyngeal samples collected from nine laboratories in Zimbabwe between March 20 and Oct 16, 2020. Samples were taken as a result of quarantine procedures for international arrivals or to test for infection in people who were symptomatic or close contacts of positive cases. Samples that had a cycle threshold of less than 30 in the diagnostic PCR test were processed for sequencing. We began our analysis in July, 2020 (120 days since the first case), with a follow-up in October, 2020 (at 210 days since the first case). The phylogenetic relationship of the genome sequences within Zimbabwe and global samples was established using maximum likelihood and Bayesian methods. FINDINGS: Of 92 299 nasopharyngeal samples collected during the study period, 8099 were PCR-positive and 328 were available for sequencing, with 156 passing sequence quality control. 83 (53%) of 156 were from female participants. At least 26 independent introductions of SARS-CoV-2 into Zimbabwe in the first 210 days were associated with 12 global lineages. 151 (97%) of 156 had the Asp614Gly mutation in the spike protein. Most cases, 93 (60%), were imported from outside Zimbabwe. Community transmission was reported 6 days after the onset of the outbreak. INTERPRETATION: Initial public health interventions delayed onset of SARS-CoV-2 community transmission after the introduction of the virus from international and regional migration in Zimbabwe. Global whole genome sequence data are essential to reveal major routes of spread and guide intervention strategies. FUNDING: WHO, Africa CDC, Biotechnology and Biological Sciences Research Council, Medical Research Council, National Institute for Health Research, and Genome Research Limited.

CoronaHiT: high-throughput sequencing of SARS-CoV-2 genomes.

We present CoronaHiT, a platform and throughput flexible method for sequencing SARS-CoV-2 genomes (≤ 96 on MinION or > 96 on Illumina NextSeq) depending on changing requirements experienced during the pandemic. CoronaHiT uses transposase-based library preparation of ARTIC PCR products. Method performance was demonstrated by sequencing 2 plates containing 95 and 59 SARS-CoV-2 genomes on nanopore and Illumina platforms and comparing to the ARTIC LoCost nanopore method. Of the 154 samples sequenced using all 3 methods, ≥ 90% genome coverage was obtained for 64.3% using ARTIC LoCost, 71.4% using CoronaHiT-ONT and 76.6% using CoronaHiT-Illumina, with almost identical clustering on a maximum likelihood tree. This protocol will aid the rapid expansion of SARS-CoV-2 genome sequencing globally.

Large-scale sequencing of SARS-CoV-2 genomes from one region allows detailed epidemiology and enables local outbreak management.

The COVID-19 pandemic has spread rapidly throughout the world. In the UK, the initial peak was in April 2020; in the county of Norfolk (UK) and surrounding areas, which has a stable, low-density population, over 3200 cases were reported between March and August 2020. As part of the activities of the national COVID-19 Genomics Consortium (COG-UK) we undertook whole genome sequencing of the SARS-CoV-2 genomes present in positive clinical samples from the Norfolk region. These samples were collected by four major hospitals, multiple minor hospitals, care facilities and community organizations within Norfolk and surrounding areas. We combined clinical metadata with the sequencing data from regional SARS-CoV-2 genomes to understand the origins, genetic variation, transmission and expansion (spread) of the virus within the region and provide context nationally. Data were fed back into the national effort for pandemic management, whilst simultaneously being used to assist local outbreak analyses. Overall, 1565 positive samples (172 per 100 000 population) from 1376 cases were evaluated; for 140 cases between two and six samples were available providing longitudinal data. This represented 42.6 % of all positive samples identified by hospital testing in the region and encompassed those with clinical need, and health and care workers and their families. In total, 1035 cases had genome sequences of sufficient quality to provide phylogenetic lineages. These genomes belonged to 26 distinct global lineages, indicating that there were multiple separate introductions into the region. Furthermore, 100 genetically distinct UK lineages were detected demonstrating local evolution, at a rate of ~2 SNPs per month, and multiple co-occurring lineages as the pandemic progressed. Our analysis: identified a discrete sublineage associated with six care facilities; found no evidence of reinfection in longitudinal samples; ruled out a nosocomial outbreak; identified 16 lineages in key workers which were not in patients, indicating infection control measures were effective; and found the D614G spike protein mutation which is linked to increased transmissibility dominates the samples and rapidly confirmed relatedness of cases in an outbreak at a food processing facility. The large-scale genome sequencing of SARS-CoV-2-positive samples has provided valuable additional data for public health epidemiology in the Norfolk region, and will continue to help identify and untangle hidden transmission chains as the pandemic evolves.

Microbial composition and dynamics in environmental samples from a ready-to-eat food production facility with a long-term colonization of Listeria monocytogenes.

Listeria monocytogenes is a foodborne pathogen of significant concern for the food industry due to its remarkable ability to persist through safety control efforts, posing a subsequent health threat to consumers. Understanding the microbial communities coexisting with L. monocytogenes in food processing environments provides insights into its persistence mechanisms. We investigated the microbial communities on non-food contact surfaces in a facility producing ready-to-eat foods, known to harbour a ST121 L. monocytogenes strain over multiple years. A 10-week sampling period was coordinated with the company and public health authorities. Metagenomic analysis revealed a stable microbial composition dominated by Pseudomonas fluorescens. While highly related populations were present in high-care production zones, distinctive taxa characteristic of specific areas were observed (e.g., Sphingomonas aerolata). Although Listeria spp. were not detected in metagenomes, they were detected in cultured samples, suggesting low relative abundance in factory settings. The findings suggest that a stable resident microbiota, with distinct adaptations to different areas within the factory, was selected for by their collective ability to survive control efforts in this environment. Listeria spp. was a member of this microbial community, albeit at low abundance, and may likewise benefit from the mutualism of the overall microbial community.