Global Genomes Research Articles

Intricate interplay of CRISPR-Cas systems, anti-CRISPR proteins, and antimicrobial resistance genes in a globally successful multi-drug resistant Klebsiella pneumoniae clone

BackgroundKlebsiella pneumoniae is one of the most prevalent pathogens responsible for multiple infections in healthcare settings and the community. K. pneumoniae CG147, primarily including ST147 (the founder ST), ST273, and ST392, is one of the most globally successful MDR clone linked to various carbapenemases.MethodsOne hundred and one CG147 strains were sequenced and additional 911 publicly available CG147 genome sequences were included for analysis. The molecular epidemiology, population structure, and time phylogeny were investigated. The virulome, resistome, and mobilome were analyzed, and the recombination in the capsular region was studied. The CRISPR-Cas and anti-CRISPR were identified. The interplay between CRISPR-Cas, anti-CRISPR, and carbapenemase-encoding plasmids was analyzed and experimentally validated.ResultsWe analyzed 1012 global CG147 genomes, with 80.4% encoding at least one carbapenemase (NDM [529/1012, 52.3%], OXA-48-like [182/1012, 17.7%], and KPC [105/1012, 10.4%]). Surprisingly, almost all CG147 strains (99.7%, 1009/1,012) harbor a chromosomal type I-E CRISPR-Cas system, with 41.8% (423/1012) containing an additional plasmid-borne type IV-A3 CRISPR-Cas system, and both target IncF plasmids, e.g., the most prevalent KPC-encoding pKpQIL-like plasmids. We found the presence of IV-A3 CRISPR-Cas system showed a negative correlation with the presence of KPC. Interestingly, a prophage-encoding anti-CRISPR AcrIE8.1 and a plasmid-borne anti-CRISPR AcrIE9.2 were detected in 40.1% (406/1012) and 54.2% (548/1012) of strains, respectively, which displayed positive correlations with the presence of a carbapenemase. Plasmid transfer experiments confirmed that the I-E and IV-A3 CRISPR-Cas systems significantly decreased (p < 0.001) KPC-encoding pKpQIL plasmid conjugation frequencies, while the AcrIE8.1 and AcrIE9.2 significantly increased (p < 0.001) pKpQIL conjugation frequencies and protected plasmids from elimination by CRISPR-Cas I-E system.ConclusionsOur results indicated a complex interplay between CRISPR-Cas, anti-CRISPR, and mobile genetic elements that shape the evolution of CG147. Our findings advance the understanding of multi-drug resistance mechanisms and will aid in preventing the emergence of future MDR clones.

Genomic and virulent characterization of a duck-associated Salmonella serovar Potsdam from China

Salmonella, a common zoonotic pathogen, is a significant concern for public health, particularly when it contaminates animal-borne products. The potential for Salmonella to infect duck embryos and disrupt their normal development not only causes substantial economic losses for the industry but also poses a severe threat to public health. However, there is a lack of understanding about the prevalence of Salmonella in duck embryos and their potential public health implications. Our study aims to fill this gap by providing genomic features of the antimicrobial resistance and virulence potential of Salmonella isolates from dead duck embryos using whole-genome sequencing and in silico toolkits. We also sought to assess the virulent characterization of the major serovar isolates by experimental infection of chicken and duck embryos. Our investigation of 195 duck embryo eggs led to the isolation of 40 (20.51%) Salmonella strains, with Salmonella serovar Potsdam being the most prevalent serovar. Most isolates were resistant to streptomycin (57.3%) and nalidixic acid (50%). Notably, our findings demonstrated that S. Potsdam exhibited a preference for ducks over chickens, suggesting potential host specificity. Additionally, global phylogenomic analysis, incorporating 180 global genomes, revealed a predominant association of S. Potsdam with ducks, supporting an adaptive process specific to the waterfowl. This study determined Salmonella serovars and antimicrobial resistance profiles in dead duck embryos, revealing a rare Salmonella serovar Potsdam with a potential for duck adaption.

Ceftriaxone-resistant Neisseria gonorrhoeae detected in England, 2015-24: an observational analysis.

Since June 2022, there has been a rise in the number of ceftriaxone-resistant Neisseria gonorrhoeae cases detected in England (n = 15), of which a third were XDR. We describe the demographic and clinical details of the recent cases and investigate the phenotypic and molecular characteristics of the isolates. For a comprehensive overview, we also reviewed 16 ceftriaxone-resistant cases previously identified in England since December 2015 and performed a global genomic comparison of all publicly available ceftriaxone-resistant N. gonorrhoeae strains with mosaic penA alleles. All N. gonorrhoeae isolates resistant to ceftriaxone (MIC > 0.125 mg/L) were whole-genome sequenced and compared with 142 global sequences of ceftriaxone-resistant N. gonorrhoeae. Demographic, behavioural and clinical data were collected. All cases were heterosexual, and most infections were associated with travel from the Asia-Pacific region. However, some had not travelled outside England within the previous few months. There were no ceftriaxone genital treatment failures, but three of five pharyngeal infections and the only rectal infection failed treatment. The isolates represented 13 different MLST STs, and most had the mosaic penA-60.001 allele. The global genomes clustered into eight major phylogroups, with regional associations. All XDR isolates belonged to the same phylogroup, represented by MLST ST16406. Most cases of ceftriaxone-resistant N. gonorrhoeae detected in England were associated with travel from the Asia-Pacific region. All genital infections were successfully treated with ceftriaxone, but there were extragenital treatment failures. Ceftriaxone resistance continues to be associated with the penA-60.001 allele within multiple genetic backgrounds and with widespread dissemination in the Asia-Pacific region.

Dissimilatory Iodate-Reducing Microorganisms Contribute to the Enrichment of Iodine in Groundwater.

Iodate reduction by dissimilatory iodate-reducing microorganisms (DIRMs) plays a crucial role in the biogeochemical cycling of iodine on Earth. However, the occurrence and distribution of DIRMs in iodine-rich groundwater remain unclear. In this study, we isolated the dissimilatory iodate-reducing bacteriumAzonexus hydrophilusstrain NCP973 from a geogenic high-iodine groundwater of China for the first time. The analysis of genome, transcriptome, and heterologous expression revealed that strain NCP973 uses the dissimilatory iodate-reducing enzyme IdrABP1P2 to reduce dissolved or in situ sediment-bound iodate to iodide. The location of IdrABP1P2 in the conjugative plasmid of strain NCP973 implies that IdrABP1P2 could be spread by horizontal gene transfer and allow the recipient microorganisms to participate in the enrichment of iodide in aquifers. Based on the global iodine-rich groundwater metagenomes and genomes, the identification of idrA showed that phylogenetically diverse DIRMs are widely distributed not only in geogenic high-iodine groundwater of China but also in radionuclide-contaminated groundwater of USA as well as in subsurface cavern waters in Germany and Italy. Moreover, the abundance of idrA was found to be higher in groundwater with a relatively high iodine content. Collectively, these results suggest that terrestrial iodine-affected groundwater systems are another important habitat for DIRMs in addition to marine environments, and their activity in aquifers triggers the mobilization and enrichment of iodine in groundwater worldwide.

Phylogeography and reassortment patterns of human influenza A viruses in sub-Saharan Africa

The role of sub-Saharan Africa in the global spread of influenza viruses remains unclear due to insufficient spatiotemporal sequence data. Here, we analyzed 222 codon-complete sequences of influenza A viruses (IAVs) sampled between 2011 and 2013 from five countries across sub-Saharan Africa (Kenya, Zambia, Mali, Gambia, and South Africa); these genomes were compared with 1209 contemporaneous global genomes using phylogeographical approaches. The spread of influenza in sub-Saharan Africa was characterized by (i) multiple introductions of IAVs into the region over consecutive influenza seasons, with viral importations originating from multiple global geographical regions, some of which persisted in circulation as intra-subtype reassortants for multiple seasons, (ii) virus transfer between sub-Saharan African countries, and (iii) virus export from sub-Saharan Africa to other geographical regions. Despite sparse data from influenza surveillance in sub-Saharan Africa, our findings support the notion that influenza viruses persist as temporally structured migrating metapopulations in which new virus strains can emerge in any geographical region, including in sub-Saharan Africa; these lineages may have been capable of dissemination to other continents through a globally migrating virus population. Further knowledge of the viral lineages that circulate within understudied sub-Saharan Africa regions is required to inform vaccination strategies in those regions.

Epidemiological and evolutionary analysis of canine circovirus from 1996 to 2023

BackgroundCanine circovirus (CanineCV), a non-enveloped virus with a circular DNA genome, has been identified in various avian and mammalian species, including domestic and wild canids. This study aimed to comprehensively analyze the prevalence of CanineCV across diverse animal species in 11 provinces of China.ResultsA total of 1,666 serum samples were collected, revealing a 5.82% prevalence of CanineCV in dogs, with the highest rates being observed in southern and eastern China. Phylogenetic analysis of 266 global CanineCV genomes sourced from the NCBI identified six distinct genotypes, elucidating the complex dynamics of their evolution. Evidence suggested a potential bat origin for CanineCV, with positive selection and high rates of evolution being observed. Recombination analysis revealed dynamic genetic exchange, highlighting the intricate nature of CanineCV evolution. Mutational analysis identified key amino acid substitutions likely to influence the virus’s adaptation. Additionally, glycosylation, palmitoylation, and SUMOylation sites were predicted, shedding light on crucial functional properties of the virus.ConclusionsThis study provides a global perspective on the origin, genetic diversity, and evolutionary dynamics of CanineCV. Understanding these factors is crucial for elucidating its epidemiology and potential health risks.

In silico gepotidacin target mining among 33 213 global Neisseria gonorrhoeae genomes from 1928 to 2023 combined with gepotidacin MIC testing of 22 gonococcal isolates with different GyrA and ParC substitutions.

The novel dual-target triazaacenaphthylene, gepotidacin, recently showed promising results in its Phase III randomized controlled trial for the treatment of gonorrhoea. We investigated alterations in the gepotidacin GyrA and ParC targets in gonococci by in silico mining of publicly available global genomes (n = 33 213) and determined gepotidacin MICs in isolates with GyrA A92 alterations combined with other GyrA and/or ParC alterations. We examined gonococcal gyrA and parC alleles available at the European Nucleotide Archive. MICs were determined using the agar dilution method (gepotidacin) or Etest (four antimicrobials). Models of DNA gyrase and topoisomerase IV were obtained from AlphaFold and used to model gepotidacin in the binding site. GyrA A92 alterations were identified in 0.24% of genomes: GyrA A92P/S/V + S91F + D95Y/A/N (0.208%), A92P + S91F (0.024%) and A92P (0.003%), but no A92T (previously associated with gepotidacin resistance) was found. ParC D86 alterations were found in 10.6% of genomes: ParC D86N/G (10.5%), D86N + S87I (0.051%), D86N + S88P (0.012%) and D86G + E91G (0.003%). One isolate had GyrA A92P + ParC D86N alterations, but remained susceptible to gepotidacin (MIC = 0.125 mg/L). No GyrA plus ParC alterations resulted in a gepotidacin MIC > 4 mg/L. Modelling of gepotidacin binding to GyrA A92/A92T/A92P suggested that gepotidacin resistance due to GyrA A92T might be linked to the formation of a new polar contact with DNA. In silico mining of 33 213 global gonococcal genomes (isolates from 1928 to 2023) showed that A92 is highly conserved in GyrA, while alterations in D86 of ParC are common. No GyrA plus ParC alterations caused gepotidacin resistance. MIC determination and genomic surveillance of potential antimicrobial resistance determinants are imperative.

Rapid identification and subsequent contextualization of an outbreak of methicillin-resistant Staphylococcus aureus in a neonatal intensive care unit using nanopore sequencing.

Outbreaks of methicillin-resistant Staphylococcus aureus (MRSA) are well described in the neonatal intensive care unit (NICU) setting. Genomics has revolutionized the investigation of such outbreaks; however, to date, this has largely been completed retrospectively and has typically relied on short-read platforms. In 2022, our laboratory established a prospective genomic surveillance system using Oxford Nanopore Technologies sequencing for rapid outbreak detection. Herein, using this system, we describe the detection and control of an outbreak of sequence-type (ST)97 MRSA in our NICU. The outbreak was identified 13 days after the first MRSA-positive culture and at a point where there were only two known cases. Ward screening rapidly defined the extent of the outbreak, with six other infants found to be colonized. There was minimal transmission once the outbreak had been detected and appropriate infection control measures had been instituted; only two further ST97 cases were detected, along with three unrelated non-ST97 MRSA cases. To contextualize the outbreak, core-genome single-nucleotide variants were identified for phylogenetic analysis after de novo assembly of nanopore data. Comparisons with global (n=45) and national surveillance (n=35) ST97 genomes revealed the stepwise evolution of methicillin resistance within this ST97 subset. A distinct cluster comprising nine of the ten ST97-IVa genomes from the NICU was identified, with strains from 2020 to 2022 national surveillance serving as outgroups to this cluster. One ST97-IVa genome presumed to be part of the outbreak formed an outgroup and was retrospectively excluded. A second phylogeny was created using Illumina sequencing, which considerably reduced the branch lengths of the NICU isolates on the phylogenetic tree. However, the overall tree topology and conclusions were unchanged, with the exception of the NICU outbreak cluster, where differences in branch lengths were observed. This analysis demonstrated the ability of a nanopore-only prospective genomic surveillance system to rapidly identify and contextualize an outbreak of MRSA in a NICU.

Prediction by genetic MATS of 4CMenB vaccine strain coverage of invasive meningococcal serogroup B isolates circulating in Taiwan between 2003 and 2020.

Neisseria meningitidis serogroup B (NmB) strains have diverse antigens, necessitating methods for predicting meningococcal serogroup B (MenB) vaccine strain coverage. The genetic Meningococcal Antigen Typing System (gMATS), a correlate of MATS estimates, predicts strain coverage by the 4-component MenB (4CMenB) vaccine in cultivable and non-cultivable NmB isolates. In Taiwan, 134 invasive, disease-causing NmB isolates were collected in 2003-2020 (23.1%, 4.5%, 5.2%, 29.8%, and 37.3% from individuals aged ≤11 months, 12-23 months, 2-4 years, 5-29 years, and ≥30 years, respectively). NmB isolates were characterized by whole-genome sequencing and vaccine antigen genotyping, and 4CMenB strain coverage was predicted using gMATS. Analysis of phylogenetic relationships with 502 global NmB genomes showed that most isolates belonged to three global hyperinvasive clonal complexes: ST-4821 (27.6%), ST-32 (23.9%), and ST-41/44 (14.9%). Predicted strain coverage by gMATS was 62.7%, with 27.6% isolates covered, 2.2% not covered, and 66.4% unpredictable by gMATS. Age group coverage point estimates ranged from 42.9% (2-4 years) to 66.1% (≤11 months). Antigen coverage estimates and percentages predicted as covered/not covered were highly variable, with higher estimates for isolates with one or more gMATS-positive antigens than for isolates positive for one 4CMenB antigen. In conclusion, this first study on NmB strain coverage by 4CMenB in Taiwan shows 62.7% coverage by gMATS, with predictable coverage for 29.8% of isolates. These could be underestimated since the gMATS calculation does not consider synergistic mechanisms associated with simultaneous antibody binding to multiple targets elicited by multicomponent vaccines or the contributions of minor outer membrane vesicle vaccine components.IMPORTANCEMeningococcal diseases, caused by the bacterium Neisseria meningitidis (meningococcus), include meningitis and septicemia. Although rare, invasive meningococcal disease is often severe and can be fatal. Nearly all cases are caused by six meningococcal serogroups (types), including meningococcal serogroup B. Vaccines are available against meningococcal serogroup B, but the antigens targeted by these vaccines have highly variable genetic features and expression levels, so the effectiveness of vaccination may vary depending on the strains circulating in particular countries. It is therefore important to test meningococcal serogroup B strains isolated from specific populations to estimate the percentage of bacterial strains that a vaccine can protect against (vaccine strain coverage). Meningococcal isolates were collected in Taiwan between 2003 and 2020, of which 134 were identified as serogroup B. We did further investigations on these isolates, including using a method (called gMATS) to predict vaccine strain coverage by the 4-component meningococcal serogroup B vaccine (4CMenB).

Exploring the genetic makeup of Xanthomonas species causing bacterial spot in Taiwan: evidence of population shift and local adaptation

The introduction of plant pathogens can quickly reshape disease dynamics in island agro-ecologies, representing a continuous challenge for local crop management strategies. Xanthomonas pathogens causing tomato bacterial spot were probably introduced in Taiwan several decades ago, creating a unique opportunity to study the genetic makeup and adaptive response of this alien population. We examined the phenotypic and genotypic identity of 669 pathogen entries collected across different regions of Taiwan in the last three decades. The analysis detected a major population shift, where X. euvesicatoria and X. vesicatoria races T1 and T2 were replaced by new races of X. perforans. After its introduction, race T4 quickly became dominant in all tomato-growing areas of the island. The genomic analysis of 317 global genomes indicates that the Xanthomonas population in Taiwan has a narrow genetic background, most likely resulting from a small number of colonization events. However, despite the apparent genetic uniformity, X. perforans race T4 shows multiple phenotypic responses in tomato lines. Additionally, an in-depth analysis of effector composition suggests diversification in response to local adaptation. These include unique mutations on avrXv3 which might allow the pathogen to overcome Xv3/Rx4 resistance gene. The findings underscore the dynamic evolution of a pathogen when introduced in a semi-isolated environment and provide insights into the potential management strategies for this important disease of tomato.

Analysis of global Aeromonas caviae genomes revealed that strains carrying T6SS are more common in human gastroenteritis than in environmental sources and are often phylogenetically related.

Aeromonas caviae is an emerging human enteric pathogen. However, the genomic features and virulence genes of A. caviae strains from human gastroenteritis and other sources have not been fully elucidated. Here, we conducted a genomic analysis of 565 global A. caviae strains isolated from different sources, including 261 strains isolated from faecal samples of gastroenteritis patients, of which 18 genomes were sequenced in this study. The presence of bacterial virulence genes and secretion systems in A. caviae strains from different sources was compared, and the phylogenetic relationship of A. caviae strains was assessed based on the core genome. The complete genome of A. caviae strain A20-9 isolated from a gastroenteritis patient was obtained in this study, from which 300 putative virulence factors and a T4SS-encoding plasmid, pAC, were identified. Genes encoding T4SS were also identified in a novel genomic island, ACI-1, from other T4SS-positive strains. The prevalence of T4SS was significantly lower in A. caviae strains from gastroenteritis patients than in environmental strains (3 %, P<0.0001 vs 14 %, P<0.01). Conversely, the prevalence of T6SS was significantly higher in A. caviae strains isolated from gastroenteritis patients than in environmental strains (25 %, P<0.05 vs 13 %, P<0.01). Four phylogenetic clusters were formed based on the core genome of 565 A. caviae strains, and strains carrying T6SS often showed close phylogenetic relationships. T3SS, aerolysin and thermostable cytotonic enterotoxin were absent in all 565 A. caviae strains. Our findings provide novel information on the genomic features of A. caviae and suggest that T6SS may play a role in A. caviae-induced human gastroenteritis.

Effect of molnupiravir on SARS-CoV-2 evolution in immunocompromised patients: a retrospective observational study

Continued SARS-CoV-2 infection among immunocompromised individuals is likely to play a role in generating genomic diversity and the emergence of novel variants. Antiviral treatments such as molnupiravir are used to mitigate severe COVID-19 outcomes, but the extended effects of these drugs on viral evolution in patients with chronic infections remain uncertain. This study investigates how molnupiravir affects SARS-CoV-2 evolution in immunocompromised patients with prolonged infections. The study included five immunocompromised patients treated with molnupiravir and four patients not treated with molnupiravir (two immunocompromised and two non-immunocompromised). We selected patients who had been infected by similar SARS-CoV-2 variants and with high-quality genomes across timepoints to allow comparison between groups. Throat and nasopharyngeal samples were collected in patients up to 44 days post treatment and were sequenced using tiled amplicon sequencing followed by variant calling. The UShER pipeline and University of California Santa Cruz genome viewer provided insights into the global context of variants. Treated and untreated patients were compared, and mutation profiles were visualised to understand the impact of molnupiravir on viral evolution. Patients treated with molnupiravir showed a large increase in low-to-mid-frequency variants in as little as 10days after treatment, whereas no such change was observed in untreated patients. Some of these variants became fixed in the viral population, including non-synonymous mutations in the spike protein. The variants were distributed across the genome and included unique mutations not commonly found in global omicron genomes. Notably, G-to-A and C-to-T mutations dominated the mutational profile of treated patients, persisting up to 44 days post treatment. Molnupiravir treatment in immunocompromised patients led to the accumulation of a distinctive pattern of mutations beyond the recommended 5 days of treatment. Treated patients maintained persistent PCR positivity for the duration of monitoring, indicating clear potential for transmission and subsequent emergence of novel variants. Australian Research Council.

Genomic investigation and nationwide tracking of pediatric invasive nontyphoidal Salmonella in China.

Invasive nontyphoidal Salmonella (iNTS) causes significant concern with ~15% morbidity, affecting populations mainly in African countries. However, iNTS infections among the Chinese pediatric population remain largely unknown. Here, we conducted a genomic investigation to study pediatric iNTS infections in a Chinese hospital. iNTS isolates accounted for 15.2% (18/119) of all nontyphoidal Salmonella (NTS) strains. Compared to non-iNTS isolates, iNTS isolates harbored a lower prevalence of antimicrobial-resistant genes of fluoroquinolones and β-lactams, as well as disinfectant determinants and plasmids, but carried a significantly higher prevalence of cdtB, faeCDE, and tcpC genes. Importantly, we detected an emerging serovar Goldcoast as the predominant iNTS serovar locally. By integrating 320 global Goldcoast genomes based on the One Health samplings, we conducted nationwide phylogenomic tracking and detected repeated human-to-human transmission events among iNTS cases caused by an underestimated serovar Goldcoast. Together, our exploratory genomic approach highlights a new trend in pediatric iNTS infections.

The rare Salmonella enterica serovar Isangi: genomic characterization of the antimicrobial resistance, virulence potential and epidemiology of Brazilian strains in comparison to global isolates.

Introduction. Salmonella enterica serovar Isangi (S. Isangi) is a rare non-typhoidal serovar, related to invasive nosocomial infections in various countries and to increasing antimicrobial resistance rates.Gap statement. Despite existing reports on S. Isangi, there is a lack of information of specific traits regarding this serovar, which could be improved through genomic analyses.Aim. Our goals were to characterize the antimicrobial resistance, virulence potential and genomic relatedness of 11 S. Isangi strains from Brazil in comparison to 185 genomes of global isolates using whole-genome sequencing (WGS) data.Methodology. Phenotypic resistance was determined by disc-diffusion. The search for resistance genes, plasmids, prophages, Salmonella pathogenicity islands (SPIs) and virulence genes, plus multi-locus sequence typing (MLST) and core-genome MLST (cgMLST) were performed using WGS.Results. Brazilian S. Isangi strains showed phenotypic resistance to nalidixic acid, ciprofloxacin and streptomycin, and harboured antimicrobial resistance [qnrB19, aac(6')-Iaa, mdsAB] and heavy metal tolerance (arsD, golST) genes. Col(pHAD28) and IncFII(S) plasmids, virulence genes related to adherence, macrophage induction, magnesium uptake, regulation and type III secretion systems, 12 SPIs and eight prophages were detected. The 185 additional global genomes analysed harboured resistance genes against 11 classes of antimicrobial compounds, 22 types of plasmids, 32 prophages, 14 SPIs, and additional virulence genes related to serum resistance, stress adaptation and toxins. Sequence type (ST)216 was assigned to genomes from Brazil and other countries, while ST335 was the most frequent ST, especially among South African genomes. cgMLST showed that Brazilian genomes were more closely related to genomes from European and African countries, the USA and Taiwan, while the majority of South African genomes were more closely related among each other.Conclusion. The presence of S. Isangi strains from Brazil and different countries showing a close genomic correlation, antimicrobial resistance profiles to drugs used in human therapy and a large number of virulence determinants reinforced the need for stronger initiatives to monitor rare non-typhoidal Salmonella serovars such as S. Isangi in order to prevent its dissemination among human and non-human sources.

Genome characterization based on the Spike-614 and NS8-84 loci of SARS-CoV-2 reveals two major possible onsets of the COVID-19 pandemic.

The global COVID-19 pandemic has lasted for 3 years since its outbreak, however its origin is still unknown. Here, we analyzed the genotypes of 3.14 million SARS-CoV-2 genomes based on the amino acid 614 of the Spike (S) and the amino acid 84 of NS8 (nonstructural protein 8), and identified 16 linkage haplotypes. The GL haplotype (S_614G and NS8_84L) was the major haplotype driving the global pandemic and accounted for 99.2% of the sequenced genomes, while the DL haplotype (S_614D and NS8_84L) caused the pandemic in China in the spring of 2020 and accounted for approximately 60% of the genomes in China and 0.45% of the global genomes. The GS (S_614G and NS8_84S), DS (S_614D and NS8_84S), and NS (S_614N and NS8_84S) haplotypes accounted for 0.26%, 0.06%, and 0.0067% of the genomes, respectively. The main evolutionary trajectory of SARS-CoV-2 is DS→DL→GL, whereas the other haplotypes are minor byproducts in the evolution. Surprisingly, the newest haplotype GL had the oldest time of most recent common ancestor (tMRCA), which was May 1 2019 by mean, while the oldest haplotype DS had the newest tMRCA with a mean of October 17, indicating that the ancestral strains that gave birth to GL had been extinct and replaced by the more adapted newcomer at the place of its origin, just like the sequential rise and fall of the delta and omicron variants. However, the haplotype DL arrived and evolved into toxic strains and ignited a pandemic in China where the GL strains had not arrived in by the end of 2019. The GL strains had spread all over the world before they were discovered, and ignited the global pandemic, which had not been noticed until the virus was declared in China. However, the GL haplotype had little influence in China during the early phase of the pandemic due to its late arrival as well as the strict transmission controls in China. Therefore, we propose two major onsets of the COVID-19 pandemic, one was mainly driven by the haplotype DL in China, the other was driven by the haplotype GL globally.

Persistence of a Frameshifting Deletion in SARS-CoV-2 ORF7a for the Duration of a Major Outbreak.

Australia experienced widespread COVID-19 outbreaks from infection with the SARS-CoV-2 Delta variant between June 2021 and February 2022. A 17-nucleotide frameshift-inducing deletion in ORF7a rapidly became represented at the consensus level (Delta-ORF7aΔ17del) in most Australian outbreak cases. Studies from early in the COVID-19 pandemic suggest that frameshift-inducing deletions in ORF7a do not persist for long in the population; therefore, Delta-ORF7aΔ17del genomes should have disappeared early in the Australian outbreak. In this study, we conducted a retrospective analysis of global Delta genomes to characterise the dynamics of Delta-ORF7aΔ17del over time, determined the frequency of all ORF7a deletions worldwide, and compared global trends with those of the Australian Delta outbreak. We downloaded all GISAID clade GK Delta genomes and scanned them for deletions in ORF7a. For each deletion we identified, we characterised its frequency, the number of countries it was found in, and how long it persisted. Of the 4,018,216 Delta genomes identified globally, 134,751 (~3.35%) possessed an ORF7a deletion, and ORF7aΔ17del was the most common. ORF7aΔ17del was the sole deletion in 28,014 genomes, of which 27,912 (~99.6%) originated from the Australian outbreak. During the outbreak, ~87% of genomes were Delta-ORF7aΔ17del, and genomes with this deletion were sampled until the outbreak's end. These data demonstrate that, contrary to suggestions early in the COVID-19 pandemic, genomes with frameshifting deletions in ORF7a can persist over long time periods. We suggest that the proliferation of Delta-ORF7aΔ17del genomes was likely a chance founder effect. Nonetheless, the frequency of ORF7a deletions in SARS-CoV-2 genomes worldwide suggests they might have some benefit for virus transmission.

GyrB in silico mining in 27 151 global gonococcal genomes from 1928-2021 combined with zoliflodacin in vitro testing of 71 international gonococcal isolates with different GyrB, ParC and ParE substitutions confirms high susceptibility.

Antimicrobial resistance (AMR) in Neisseria gonorrhoeae is a global threat and novel treatment alternatives are imperative. Herein, susceptibility to the novel antimicrobial zoliflodacin, currently in a global Phase 3 randomized controlled clinical trial for gonorrhoea treatment, was investigated by screening for zoliflodacin GyrB target mutations in publicly available gonococcal genomes and, where feasible, determination of the associated zoliflodacin MIC. The European Nucleotide Archive was queried using the search term 'Taxon: 485'. DNA sequences from 27 151 gonococcal isolates were analysed and gyrB, gyrA, parC and parE alleles characterized. GyrB amino acid alterations were rare (97.0% of isolates had a wild-type GyrB sequence). GyrB V470L (2.7% of isolates) was the most prevalent alteration, followed by S467N (0.12%), N. meningitidis GyrB (0.092%), V470I (0.059%), Q468R/P (0.015%), A466T (0.0074%), L425I + L465I (0.0037%), L465I (0.0037%), G482S (0.0037%) and D429V (0.0037%). Only one isolate (0.0037%) carried a substitution in a resistance-associated GyrB codon (D429V), resulting in a zoliflodacin MIC of 8 mg/L. None of the other detected gyrB, gyrA, parC or parE mutations caused a zoliflodacin MIC outside the wild-type MIC distribution. The zoliflodacin target GyrB was highly conserved among 27 151 global gonococcal isolates cultured in 1928-2021. The single zoliflodacin-resistant clinical isolate (0.0037%) was cultured from a male patient in Japan in 2000. Evidently, this strain has not clonally expanded nor has the gyrB zoliflodacin-resistance mutation disseminated through horizontal gene transfer to other strains. Phenotypic and genomic surveillance, including gyrB mutations, of zoliflodacin susceptibility are imperative.

Emergence of sulphonamide resistance in azithromycin-resistant pediatric strains of Salmonella Typhi and Paratyphi A: A genomics insight

Whole genome sequences of Salmonella enterica subspecies enterica serovar Typhi (S. Typhi) and Salmonella enterica subspecies enterica serovar Paratyphi A (S. Paratyphi A) from pediatric settings were used to assess the emerging Antimicrobial Resistance (AMR). The high throughput sequences of twenty pediatric clinical isolates of S. Typhi and S. Paratyphi A were retrieved and were screened for prevalent Antimicrobial Resistance Genes (ARGs) and Virulent Factors (VF). The resistance data was compared with the reference strains of S. Typhi and S. Paratyphi A. AMR studies identified sul1, sul2, dfrA7, tem-1, AH(6)-Id and APH(3″)-Ib as common ARGs. VFs were identified to understand the level of pathogenicity. The most prevalent AMR genes in the sequenced genomes were detected in phenotypically azithromycin-resistant S. Typhi. Correlation with the global genomes projected a trend of concurrent resistance to macrolides, β lactams, fluoroquinolones (FQs), tetracyclines, ansamycins, and aminoglycosides. Traces of sulphonamide-resistance were observed indicating the emergence of a currently non-prevalent S. Typhi resistance that could be a future threat. Hence new antibiotic regimen to treat azithromycin-resistant S. Typhi should be formulated by avoiding the risks of aggravating sulphonamide resistance. The identified ARGs in genomes from paediatric isolates will aid future studies to design anti-bacterial compounds against S. Typhi and S. Paratyphi A.

Thousands of small, novel genes predicted in global phage genomes.

Small genes (<150 nucleotides) have been systematically overlooked in phage genomes. We employ a large-scale comparative genomics approach to predict >40,000 small-gene families in ∼2.3 million phage genome contigs. We find that small genes in phage genomes are approximately 3-fold more prevalent than in host prokaryotic genomes. Our approach enriches for small genes that are translated in microbiomes, suggesting the small genes identified are coding. More than 9,000 families encode potentially secreted or transmembrane proteins, more than 5,000 families encode predicted anti-CRISPR proteins, and more than 500 families encode predicted antimicrobial proteins. By combining homology and genomic-neighborhood analyses, we reveal substantial novelty and diversity within phage biology, including small phage genes found in multiple host phyla, small genes encoding proteins that play essential roles in host infection, and small genes that share genomic neighborhoods and whose encoded proteins may share related functions.

Comparative genome analysis of Pasteurella multocida serogroup B:2 strains causing haemorrhagic septicaemia (HS) in bovines

Pasteurella multocida, a Gram-negative bacterium with ubiquitous nature, is known to affect wide range of host species worldwide with varied clinical manifestations including haemorrhagic septicaemia (HS) in bovines. Although, HS causing P. multocida strains were identified and characterized by conventional tools and PCR assays, diverse strains are indistinguishable by these tools in the face of disease outbreaks. In this study, draft genomes of three virulent P. multocida serotype B:2 strains (NIVEDIPm32, NIVEDIPm34 and NIVEDIPm35) were analyzed following whole genome sequencing, assembly, annotation and compared them with existing global genomes (n = 43) of bovine origin in the database. Three draft genomes of NIVEDIPm strains consisted of 40–52 contigs with GC content of ∼40.4%. The genome size and predicted genes content was ∼2.3 Mb and 2181–2189, respectively. Besides, the presence of various mobile genetic elements, antimicrobial resistance genes and biofilm related genes suggested their vital roles in virulence; further, adaptation to the host immune system as well as host pathogen interaction. Multi locus sequence analysis based on RIRDC scheme showed the presence of ST122 in all the three strains. wgMLST based phylogenic analysis suggested that HS causing Indian virulent field strains differed geographically and showed diversity from existing HS vaccine strain P52. The phylogenetic tree revealed that North Indian strains share high similarity with strains of Pakistan than South Indian Strain. Notably, a high divergence of SNPs between the HS causing circulating virulent strains of India and current HS vaccine strain P52 suggested an imminent need for relook in to HS vaccination strategy for livestock in India.