Congruity of genomic and epidemiological data in modelling of local cholera outbreaks.

Cholera continues to be a global health threat. Understanding how cholera spreads between locations is fundamental to the rational, evidence-based design of intervention and control efforts. Traditionally, cholera transmission models have used cholera case-count data. More recently, whole-genome sequence data have qualitatively described cholera transmission. Integrating these data streams may provide much more accurate models of cholera spread; however, no systematic analyses have been performed so far to compare traditional case-count models to the phylodynamic models from genomic data for cholera transmission. Here, we use high-fidelity case-count and whole-genome sequencing data from the 1991 to 1998 cholera epidemic in Argentina to directly compare the epidemiological model parameters estimated from these two data sources. We find that phylodynamic methods applied to cholera genomics data provide comparable estimates that are in line with established methods. Our methodology represents a critical step in building a framework for integrating case-count and genomic data sources for cholera epidemiology and other bacterial pathogens.

Genomic Characterization of Two NDM-5-Producing Isolates of Klebsiella pneumoniae ST11 from a Single Patient.

Two metallo-ß-lactamase-producing Klebsiella pneumoniae (HA30 and HA31) were isolated in a hospital in Argentina during 2018. K. pneumoniae HA30 was isolated from a rectal swab during the epidemiological surveillance for carbapenemase-producing strains, while K. pneumoniae HA31 was collected from the same patient 4 days after hospitalization. The aim of the present study was to identify the clonal relationships and resistome of these two NDM-producing K. pneumoniae strains isolated from a patient with a fatal outcome. Whole-genome sequencing (WGS) was performed using Illumina MiSeq-I, and subsequent analysis involved genome assembly, annotation, antibiotic resistance gene identification, multilocus sequence typing (MLST), and plasmid characterization using bioinformatics tools. Conjugation assays to E. coli J53 was conducted as previously described. K. pneumoniae HA30 exhibited extensively drug-resistant phenotype, while HA31 was multidrug-resistant as defined by Magiorakos et al., including both resistance to carbapenems, aminoglycosides and ciprofloxacin with blaNDM-5, blaCTX-M-15 and rmtB genes found in both strains. MLST analysis showed that both strains belonged to ST11, differing by only 4 cgSNPs, indicating that K. pneumoniae HA30 and HA31 were the same strain. Conjugation assays revealed that K. pneumoniae HA31 strain possessed a transferable plasmid to E. coli J53. Bioinformatics studies identified that the same strain colonizing an inpatient during hospital admission subsequently caused the infection leading to a fatal outcome, being the first report of blaNDM-5, rmtB and blaCTX-M-15 genes in a K. pneumoniae ST11 strain from Latin America. Our results also highlighted the importance of focusing on epidemiological surveillance programs.

[Molecular characterization of Listeria monocytogenes isolates from human and food sources in Argentina, 2018-2023]. / Caracterización molecular de aislamientos de Listeria monocytogenes de origen humano y alimentario en Argentina, 2018-2023.

Human listeriosis is an infectious disease caused by Listeria monocytogenes. The invasive form of this disease leads to a high rate of hospitalizations and fatality. The main mode of transmission is through contaminated ready-to-eat foods such as dairy, vegetables and meat products. The knowledge of the diversity and population dynamics of isolates collected from human and food sources is essential for the detection of clusters and the identification of common sites of infection. The aim of this study was the molecular characterization of L. monocytogenes isolates in Argentina. We sequenced a total of 63 isolates, 35 from human and 28 from food sources, collected between 2018 and 2023. Our genomic study divided the isolates into two lineages, four serogroups, 17 sequence types and 15 clonal complexes (CCs). The hypervirulent clone CC1 (lineage I; serogroup IVb) predominated in human and food samples. The phylogenomic analysis showed a high and possible epidemiological relationship between isolates from human and/or food sources, suggesting the presence of transmission chains in our country. These findings highlight the need to strengthen genomic surveillance of L. monocytogenes in Argentina. The identification of geographic distribution and characteristics of predominant and emerging clones from human and food sources might help to focus action plans and public health policies better directed at the control and prevention of listeriosis.

Emergence of monkeypox virus in central Argentina: Epidemiological features and first complete genome sequences in the country.

Monkeypox (Mpox) is a zoonotic disease caused by the monkeypox virus (MPXV). MPXV can be transmitted by close contact with lesions, body fluids, respiratory droplets, and contaminated materials. A new pattern of spread among sexual networks has been recently described. The present work aimed to report the epidemiological and genomic characterization of the 2022 MPXV outbreak in central Argentina. A total of 113 scabs and/or lesion swab specimens were studied. MPXV infection was confirmed in 46.0% of the studied patients, all of whom were men. Varicella-zoster virus infection was the most frequent differential diagnosis. Eight complete viral genomes were obtained by next-generation sequencing. The Argentinian sequences were grouped intermingled with other sequences from the 2022 MPXV outbreak, related to samples from the USA, Europe, and Peru. Taken together, our study provided an initial assessment of the genetic and epidemiological characteristics of the 2022 MPXV outbreak in Córdoba, Argentina.

"Super-Spreaders" and Person-to-Person Transmission of Andes Virus in Argentina.

BACKGROUND: From November 2018 through February 2019, person-to-person transmission of Andes virus (ANDV) hantavirus pulmonary syndrome occurred in Chubut Province, Argentina, and resulted in 34 confirmed infections and 11 deaths. Understanding the genomic, epidemiologic, and clinical characteristics of person-to-person transmission of ANDV is crucial to designing effective interventions. METHODS: Clinical and epidemiologic information was obtained by means of patient report and from public health centers. Serologic testing, contact-tracing, and next-generation sequencing were used to identify ANDV infection as the cause of this outbreak of hantavirus pulmonary syndrome and to reconstruct person-to-person transmission events. RESULTS: After a single introduction of ANDV from a rodent reservoir into the human population, transmission was driven by 3 symptomatic persons who attended crowded social events. After 18 cases were confirmed, public health officials enforced isolation of persons with confirmed cases and self-quarantine of possible contacts; these measures most likely curtailed further spread. The median reproductive number (the number of secondary cases caused by an infected person during the infectious period) was 2.12 before the control measures were enforced and decreased to 0.96 after the measures were implemented. Full genome sequencing of the ANDV strain involved in this outbreak was performed with specimens from 27 patients and showed that the strain that was present (Epuyén/18-19) was similar to the causative strain (Epilink/96) in the first known person-to-person transmission of hantavirus pulmonary syndrome caused by ANDV, which occurred in El Bolsón, Argentina, in 1996. Clinical investigations involving patients with ANDV hantavirus pulmonary syndrome in this outbreak revealed that patients with a high viral load and liver injury were more likely than other patients to spread infection. Disease severity, genomic diversity, age, and time spent in the hospital had no clear association with secondary transmission. CONCLUSIONS: Among patients with ANDV hantavirus pulmonary syndrome, high viral titers in combination with attendance at massive social gatherings or extensive contact among persons were associated with a higher likelihood of transmission. (Funded by the Ministerio de Salud y Desarrollo Social de la Nación Argentina and others.).

First molecular characterization of Escherichia coli O157:H7 isolates from clinical samples in Paraguay using whole-genome sequencing.

Escherichia coli O157:H7 is a foodborne pathogen implicated in numerous outbreaks worldwide that has the ability to cause extra-intestinal complications in humans. The Enteropathogens Division of the Central Public Health Laboratory (CPHL) in Paraguay is working to improve the genomic characterization of Shiga toxin-producing E. coli (STEC) to enhance laboratory-based surveillance and investigation of foodborne disease outbreaks. Whole genome sequencing (WGS) is proposed worldwide to be used in the routine laboratory as a high-resolution tool that allows to have all the results in a single workflow. This study aimed to carry out for the first time, the genomic characterization by WGS of nine STEC O157:H7 strains isolated from human samples in Paraguay. We were able to identify virulence and resistance mechanisms, MLST subtype, and even establish the phylogenetic relationships between isolates. Furthermore, we detected the presence of strains belonging to hypervirulent clade 8 in most of the isolates studied.

Genomic Epidemiology of Azithromycin-Nonsusceptible Neisseria gonorrhoeae, Argentina, 2005-2019.

Azithromycin-nonsusceptible Neisseria gonorrhoeae strains are an emerging global public health threat. During 2015-2018, the prevalence of azithromycin-nonsusceptible gonococcal infection increased significantly in Argentina. To investigate the genomic epidemiology and resistance mechanisms of these strains, we sequenced 96 nonsusceptible isolates collected in Argentina during 2005-2019. Phylogenomic analysis revealed 2 main clades, which were characterized by a limited geographic distribution, circulating during January 2015-November 2019. These clades included the internationally spreading multilocus sequence types (STs) 1580 and 9363. The ST1580 isolates, which had MICs of 2-4 µg/mL, had mutations in the 23S rRNA. The ST9363 isolates, which had MICs of 2-4 or >256 µg/mL, had mutations in the 23S rRNA, a mosaic mtr locus, or both. Identifying the geographic dissemination and characteristics of these predominant clones will guide public health policies to control the spread of azithromycin-nonsusceptible N. gonorrhoeae in Argentina.

Five-year microevolution of a multidrug-resistant Mycobacterium tuberculosis strain within a patient with inadequate compliance to treatment.

BACKGROUND: Whole-genome sequencing has shown that the Mycobacterium tuberculosis infection process can be more heterogeneous than previously thought. Compartmentalized infections, exogenous reinfections, and microevolution are manifestations of this clonal complexity. The analysis of the mechanisms causing the microevolution -the genetic variability of M. tuberculosis at short time scales- of a parental strain into clonal variants with a patient is a relevant issue that has not been yet completely addressed. To our knowledge, a whole genome sequence microevolution analysis in a single patient with inadequate adherence to treatment has not been previously reported. CASE PRESENTATION: In this work, we applied whole genome sequencing analysis for a more in-depth analysis of the microevolution of a parental Mycobacterium tuberculosis strain into clonal variants within a patient with poor treatment compliance in Argentina. We analyzed the whole-genome sequence of 8 consecutive Mycobacterium tuberculosis isolates obtained from a patient within 57-months of intermittent therapy. Nineteen mutations (9 short-term, 10 fixed variants) emerged, most of them associated with drug resistance. The first isolate was already resistant to isoniazid, rifampicin, and streptomycin, thereafter the strain developed resistance to fluoroquinolones and pyrazinamide. Surprisingly, isolates remained susceptible to the pro-drug ethionamide after acquiring a frameshift mutation in ethA, a gene required for its activation. We also found a novel variant, (T-54G), in the 5' untranslated region of whiB7 (T-54G), a region allegedly related to kanamycin resistance. Notably, discrepancies between canonical and phage-based susceptibility testing to kanamycin were previously found for the isolate harboring this mutation. In our patient, microevolution was mainly driven by drug selective pressure. Rare short-term mutations fixed together with resistance-conferring mutations during therapy. CONCLUSIONS: This report highlights the relevance of whole-genome sequencing analysis in the clinic for characterization of pre-XDR and MDR resistance profile, particularly in patients with incomplete and/or intermittent treatment.

Economic evaluation of whole genome sequencing for pathogen identification and surveillance - results of case studies in Europe and the Americas 2016 to 2019.

BackgroundWhole genome sequencing (WGS) is increasingly used for pathogen identification and surveillance.AimWe evaluated costs and benefits of routine WGS through case studies at eight reference laboratories in Europe and the Americas which conduct pathogen surveillance for avian influenza (two laboratories), human influenza (one laboratory) and food-borne pathogens (five laboratories).MethodsThe evaluation focused on the institutional perspective, i.e. the 'investment case' for implementing WGS compared with conventional methods, based on costs and benefits during a defined reference period, mostly covering at least part of 2017. A break-even analysis estimated the number of cases of illness (for the example of Salmonella surveillance) that would need to be avoided through WGS in order to 'break even' on costs.ResultsOn a per-sample basis, WGS was between 1.2 and 4.3 times more expensive than routine conventional methods. However, WGS brought major benefits for pathogen identification and surveillance, substantially changing laboratory workflows, analytical processes and outbreaks detection and control. Between 0.2% and 1.1% (on average 0.7%) of reported salmonellosis cases would need to be prevented to break even with respect to the additional costs of WGS.ConclusionsEven at cost levels documented here, WGS provides a level of additional information that more than balances the additional costs if used effectively. The substantial cost differences for WGS between reference laboratories were due to economies of scale, degree of automation, sequencing technology used and institutional discounts for equipment and consumables, as well as the extent to which sequencers are used at full capacity.

Phylogeographic Analysis Reveals Multiple International transmission Events Have Driven the Global Emergence of Escherichia coli O157:H7.

BACKGROUND: Shiga toxin-producing Escherchia coli (STEC) O157:H7 is a zoonotic pathogen that causes numerous food and waterborne disease outbreaks. It is globally distributed, but its origin and the temporal sequence of its geographical spread are unknown. METHODS: We analyzed whole-genome sequencing data of 757 isolates from 4 continents, and performed a pan-genome analysis to identify the core genome and, from this, extracted single-nucleotide polymorphisms. A timed phylogeographic analysis was performed on a subset of the isolates to investigate its worldwide spread. RESULTS: The common ancestor of this set of isolates occurred around 1890 (1845-1925) and originated from the Netherlands. Phylogeographic analysis identified 34 major transmission events. The earliest were predominantly intercontinental, moving from Europe to Australia around 1937 (1909-1958), to the United States in 1941 (1921-1962), to Canada in 1960 (1943-1979), and from Australia to New Zealand in 1966 (1943-1982). This pre-dates the first reported human case of E. coli O157:H7, which was in 1975 from the United States. CONCLUSIONS: Inter- and intra-continental transmission events have resulted in the current international distribution of E. coli O157:H7, and it is likely that these events were facilitated by animal movements (eg, Holstein Friesian cattle). These findings will inform policy on action that is crucial to reduce the further spread of E. coli O157:H7 and other (emerging) STEC strains globally.

PulseNet Latin America and the Caribbean Network: Present and Future.

PulseNet Latin America and Caribbean (PNLAC) was established in 2003 and is one of seven Regional networks within PulseNet International. The main objectives of the network are to strengthen national and regional laboratory-based foodborne disease surveillance for early detection and investigation of outbreaks to setup control and prevention strategies in contribution to Public Health. Participants perform standardized pulsed-field gel electrophoresis (PFGE) protocols and analysis. For functioning, it is important for the network, the development of national and regional databases (RDBs) and the communication between countries, regionally and internationally. Metadata from over 8600 cases/outbreaks are profiled and isolated by PFGE and are incorporated into the RDB, hosted by the Pan American Health Organization. Currently PNLAC is moving toward whole-genome sequencing to use as a complementary strategy for surveillance. The aim of this article was to describe the experience of the construction of PNLAC, and its contribution to the surveillance of the foodborne diseases at the country and regional levels.

PulseNet International: Vision for the implementation of whole genome sequencing (WGS) for global food-borne disease surveillance.

PulseNet International is a global network dedicated to laboratory-based surveillance for food-borne diseases. The network comprises the national and regional laboratory networks of Africa, Asia Pacific, Canada, Europe, Latin America and the Caribbean, the Middle East, and the United States. The PulseNet International vision is the standardised use of whole genome sequencing (WGS) to identify and subtype food-borne bacterial pathogens worldwide, replacing traditional methods to strengthen preparedness and response, reduce global social and economic disease burden, and save lives. To meet the needs of real-time surveillance, the PulseNet International network will standardise subtyping via WGS using whole genome multilocus sequence typing (wgMLST), which delivers sufficiently high resolution and epidemiological concordance, plus unambiguous nomenclature for the purposes of surveillance. Standardised protocols, validation studies, quality control programmes, database and nomenclature development, and training should support the implementation and decentralisation of WGS. Ideally, WGS data collected for surveillance purposes should be publicly available, in real time where possible, respecting data protection policies. WGS data are suitable for surveillance and outbreak purposes and for answering scientific questions pertaining to source attribution, antimicrobial resistance, transmission patterns, and virulence, which will further enable the protection and improvement of public health with respect to food-borne disease.

Diagnosis of an outbreak of Salmonella typhimurium in chinchillas (Chinchilla lanigera) by pulsed-field gel electrophoresis.

Adult chinchillas (Chinchilla lanigera) that had suddenly died in a commercial farm located in La Plata City, Buenos Aires Province, Argentina, in July 2012 were macroscopically, histopathologically, and microbiologically examined. Salmonella enterica serovar Typhimurium (S. Typhimurium) was isolated from the liver, spleen, heart, lungs, kidneys and intestines from each of the five animals evaluated. The five strains were susceptible to ampicillin, cephalotin, cefotaxime, nalidixic acid, gentamicin, streptomycin, chloramphenicol, fosfomycin, nitrofurantoin and trimethoprim-sulfamethoxazole, and resistant to tetracycline. Each of the five S. Typhimurium isolates was analyzed by XbaI- pulsed-field gel electrophoresis (PFGE), showing an identical electrophoretic profile with 15 defined bands, which was found to be identical to pattern ARJPXX01.0220 of the PulseNet Argentine National database of Salmonella PFGE patterns. This is the first work describing the postmortem diagnosis of an outbreak of salmonellosis in chinchillas by using molecular methods such as PFGE.

Genetic and Phenotypic Changes Related to the Development of mec-Independent Oxacillin Non-Susceptibility in ST8 Staphylococcus aureus Recovered after Antibiotic Therapy in a Patient with Bacteremia.

The mec-independent oxacillin non-susceptible S. aureus (MIONSA) strains represent a great clinical challenge, as they are not easily detected and can lead to treatment failure. However, the responsible molecular mechanisms are still very little understood. Here, we studied four clinical ST8-MSSA-t024 isolates recovered during the course of antibiotic treatment from a patient suffering successive episodes of bacteremia. The first isolates (SAMS1, SAMS2, and SAMS3) were susceptible to cefoxitin and oxacillin. The last one (SA2) was susceptible to cefoxitin, resistant to oxacillin, lacked mec genes, and had reduced susceptibility to teicoplanin. SA2 showed higher ß-lactamase activity than SAMS1. However, ß-lactamase hyperproduction could not be linked to oxacillin resistance as it was not inhibited by clavulanic acid, and no genetic changes that could account for its hyperproduction were found. Importantly, we hereby report the in vivo acquisition and coexistence of different adaptive mutations in genes associated with peptidoglycan synthesis (pbp2, rodA, stp1, yjbH, and yvqF/vraT), which is possibly related with the development of oxacillin resistance and reduced susceptibility to teicoplanin in SA2. Using three-dimensional models and PBP binding assays, we demonstrated the high contribution of the SA2 PBP2 Ala450Asp mutation to the observed oxacillin resistance phenotype. Our results should be considered as a warning for physicians and microbiologists in the region, as MIONSA detection and treatment represent an important clinical challenge.

Revisiting typing systems for group B Streptococcus prophages: an application in prophage detection and classification in group B Streptococcus isolates from Argentina.

Group B Streptococcus (GBS) causes severe infections in neonates and adults with comorbidities. Prophages have been reported to contribute to GBS evolution and pathogenicity. However, no studies are available to date on the presence and diversity of prophages in GBS isolates from humans in South America. This study provides insights into the prophage content of 365 GBS isolates collected from clinical samples in the context of an Argentinean multicentric study. Using whole-genome sequence data, we implemented two previously proposed methods for prophage typing: a PCR approach (carried out in silico) coupled with a blastx-based method to classify prophages based on their prophage group and integrase type, respectively. We manually searched the genomes and identified 325 prophages. However, only 80% of prophages could be accurately categorized with the previous approaches. Integration of phylogenetic analysis, prophage group and integrase type allowed for all to be classified into 19 prophage types, which correlated with GBS clonal complex grouping. The revised prophage typing approach was additionally improved by using a blastn search after enriching the database with ten new genes for prophage group classification combined with the existing integrase typing method. This modified and integrated typing system was applied to the analysis of 615 GBS genomes (365 GBS from Argentina and 250 from public databases), which revealed 29 prophage types, including two novel integrase subtypes. Their characterization and comparative analysis revealed major differences in the lysogeny and replication modules. Genes related to bacterial fitness, virulence or adaptation to stressful environments were detected in all prophage types. Considering prophage prevalence, distribution and their association with bacterial virulence, it is important to study their role in GBS epidemiology. In this context, we propose the use of an improved and integrated prophage typing system suitable for rapid phage detection and classification with little computational processing.

Filling the gaps in the Argentinian distribution of orthohantavirus: First finding of Lechiguanas virus in rodents from Corrientes, Argentina.

Hantavirus Pulmonary Syndrome (HPS) is an emerging infectious disease caused by orthohantaviruses in the Americas. In Argentina, since 1995, several reservoirs and virus variants have been described, but the northeastern and central endemic zones in the country include an area without human or rodent infections, despite sharing rodent species with areas with that disease. The aim of this study was to search for orthohantavirus in rodent communities that inhabit this area, which borders two endemic areas of HPS. Small rodents were captured in June of 2022 through a total effort of 644 trap nights distributed in five grids located in the Iberá National Park, Corrientes, Northeastern Argentina. All rodents were sexed, weighed, and the species was recorded. Blood samples were extracted to detect ANDV-specific immunoglobulin G (IgG), and to extract the RNA virus. Trimmed sequences were mapped against reference sequences from GenBank. We captured a total of 36 Oligoryzomys flavescens and 15 Oxymycterus rufus. We detected the O. flavescens species infected with Lechiguanas orthohantavirus in the camping area of the National Park. A nucleotide comparison with previously published sequences shows a 98.34% similarity to the virus obtained from a human case of HPS reported in the adjacent Misiones province. This study demonstrated, for the first time, that O. flavescens is a host of the Lechiguanas orthohantavirus in this zone and contributes to closing information gaps on the distribution of orthohantavirus in Argentina. Additionally, the high similarity with the hantavirus found in the human case of Misiones suggests that the reservoir in that province would also be O. flavescens (not previously confirmed). This information permits us to focus on the preventive measurements to protect the human population.

Phytobacter spp: the emergence of a new genus of healthcare-associated Enterobacterales encoding carbapenemases in Argentina: a case series.

Members of the genus Phytobacter (order Enterobacterales) are isolated from the natural environment and clinical settings. Identification of Phytobacter strains based on biochemical characteristics is complicated due to taxonomic confusion, and they are often misidentified by automated identification systems in laboratories. In this study we describe the first three clinical cases associated with Phytobacter spp. reported in Argentina. We describe the identification, the molecular analysis using whole genome sequencing and the potential clinical relevance.

Genetic characterization and estimated 4CMenB vaccine strain coverage of 284 Neisseria meningitidis isolates causing invasive meningococcal disease in Argentina in 2010-2014.

Meningococcal (Neisseria meningitidis) serogroup B (MenB) strain antigens are diverse and a limited number of strains can be evaluated using the human serum bactericidal antibody (hSBA) assay. The genetic Meningococcal Antigen Typing System (gMATS) was developed to predict the likelihood of coverage for large numbers of isolates by the 4CMenB vaccine, which includes antigens Neisseria adhesin A (NadA), Neisserial Heparin-Binding Antigen (NHBA), factor H-binding protein (fHbp), and Porin A (PorA). In this study, we characterized by whole-genome analyses 284 invasive MenB isolates collected from 2010 to 2014 by the Argentinian National Laboratories Network (52-61 isolates per year). Strain coverage was estimated by gMATS on all isolates and by hSBA assay on 74 randomly selected isolates, representative of the whole panel. The four most common clonal complexes (CCs), accounting for 81.3% of isolates, were CC-865 (75 isolates, 26.4%), CC-32 (59, 20.8%), CC-35 (59, 20.8%), and CC-41/44 (38, 13.4%). Vaccine antigen genotyping showed diversity. The most prevalent variants/peptides were fHbp variant 2, NHBA peptides 24, 21, and 2, and PorA variable region 2 profiles 16-36 and 14. The nadA gene was present in 66 (23.2%) isolates. Estimated strain coverage by hSBA assay showed 78.4% of isolates were killed by pooled adolescent sera, and 51.4% and 64.9% (based on two different thresholds) were killed by pooled infant sera. Estimated coverage by gMATS (61.3%; prediction interval: 55.5%, 66.7%) was consistent with the infant hSBA assay results. Continued genomic surveillance is needed to evaluate the persistence of major MenB CCs in Argentina.

The most common clinical manifestations of invasive meningococcal disease include meningitis and septicemia, which can be deadly, and many survivors suffer long-term serious after-effects. Most cases of invasive meningococcal disease are caused by six meningococcal serogroups (types), including serogroup B. Although vaccines are available against meningococcal serogroup B infection, these vaccines target antigens that are highly diverse. Consequently, the effectiveness of vaccination may vary from country to country because the meningococcal serogroup B strains circulating in particular regions carry different forms of the target vaccine antigens. This means it is important to test serogroup B strains isolated from specific populations to estimate the percentage of strains that a vaccine is likely to be effective against (known as 'vaccine strain coverage'). The genetic Meningococcal Antigen Typing System (gMATS) was developed to predict strain coverage by the four-component meningococcal serogroup B vaccine, 4CMenB, against large numbers of serogroup B strains. In this study, we analyzed 284 invasive meningococcal serogroup B isolates collected between 2010 and 2014 in Argentina. Genetic analyses showed that the vaccine antigens of the isolates were diverse and some genetic characteristics had not been found in isolates from other countries. However, vaccine strain coverage estimated by gMATS was consistent with that reported in other parts of the world and with strain coverage results obtained for a subset via another method, the human serum bactericidal antibody (hSBA) assay. These results highlight the need for continued monitoring of circulating bacterial strains to assess the estimated strain coverage of meningococcal serogroup B vaccines.

Addressing pandemic-wide systematic errors in the SARS-CoV-2 phylogeny.

The SARS-CoV-2 genome occupies a unique place in infection biology - it is the most highly sequenced genome on earth (making up over 20% of public sequencing datasets) with fine scale information on sampling date and geography, and has been subject to unprecedented intense analysis. As a result, these phylogenetic data are an incredibly valuable resource for science and public health. However, the vast majority of the data was sequenced by tiling amplicons across the full genome, with amplicon schemes that changed over the pandemic as mutations in the viral genome interacted with primer binding sites. In combination with the disparate set of genome assembly workflows and lack of consistent quality control (QC) processes, the current genomes have many systematic errors that have evolved with the virus and amplicon schemes. These errors have significant impacts on the phylogeny, and therefore over the last few years, many thousands of hours of researchers time has been spent in "eyeballing" trees, looking for artefacts, and then patching the tree. Given the huge value of this dataset, we therefore set out to reprocess the complete set of public raw sequence data in a rigorous amplicon-aware manner, and build a cleaner phylogeny. Here we provide a global tree of 3,960,704 samples, built from a consistently assembled set of high quality consensus sequences from all available public data as of March 2023, viewable at https://viridian.taxonium.org. Each genome was constructed using a novel assembly tool called Viridian (https://github.com/iqbal-lab-org/viridian), developed specifically to process amplicon sequence data, eliminating artefactual errors and mask the genome at low quality positions. We provide simulation and empirical validation of the methodology, and quantify the improvement in the phylogeny. Phase 2 of our project will address the fact that the data in the public archives is heavily geographically biased towards the Global North. We therefore have contributed new raw data to ENA/SRA from many countries including Ghana, Thailand, Laos, Sri Lanka, India, Argentina and Singapore. We will incorporate these, along with all public raw data submitted between March 2023 and the current day, into an updated set of assemblies, and phylogeny. We hope the tree, consensus sequences and Viridian will be a valuable resource for researchers.