ABSTRACT Epitranscriptomics has recently gained significant momentum due to technological advances and translational applications, however, studies on bacterial RNA modifications remain limited. Bacterial RNA remains notoriously prone to degradation and methodologies to investigate the epitranscriptome are challenging. Prior research has shown RNA modifications modulate antimicrobial resistance, virulence and pathogenicity. This research employed CRISPR interference to knock down five known Escherichia coli rRNA modification genes ( rlmF , rlmJ , rluD , rsmF and rsmG ) in three E. coli strains. These isolates underwent growth curves, proteome analysis and native RNA sequencing CRISPRi adequately silenced the majority of RNA modification genes in E. coli (>80% reduction). Significant growth delays were associated with rlmF , rsmF and rsmG repression. Unique protein pathways corresponding with RNA modification loss were found for rlmJ (TreB, XylF), rluD (CysH, HycB, PutP, TrpB), rsmF (EvgA) and rsmG (OppC). Known rRNA modification sites for rluD (Ψ) and rsmG (m7G) were detected from analysis of nanopore electrical signal, however, only a weak signal was apparent for m6A ( rlmF, rlmJ ) and m5C ( rsmF ) modifications. The inhibition of rRNA modifications resulted in mRNA modification changes including for genes ompC , cspC , dbhA , dbhB and secY . Our work provides an approach for unravelling the epitranscriptome of E. coli and gain insight into its functional role.

Burkholderia pseudomallei (Bp), an environmental bacterium and opportunistic pathogen endemic to tropical regions, is highly adaptive and thrives in diverse environments, from soil to human hosts. Bacterial adaptation is critical for survival, virulence modulation, and persistence during infection and can manifest as colony morphotype variation (CMV). Although Bp adaptation has been studied, CMV remains poorly understood. Here, we characterized five clinical Bp isolates exhibiting heterogeneous populations with rough and smooth colony morphologies. We used phenotypic assays, whole-genome sequencing, and proteomics to investigate the molecular pathways reflecting CMV, by comparing smooth and rough morphotypes. Although phenotypic differences in protease activity, hemolysis, mucoidy, iron uptake, and antibiotic sensitivity-including to antimicrobial agents commonly used to treat infections-were rare, these traits alone could not distinguish morphotypes or groups of isolates. Genomic comparisons revealed either no differences or limited isolate-specific mutations, which do not explain the overall difference in phenotypes. In contrast, proteomic analysis uncovered consistent shifts in protein abundance related to virulence, including quorum sensing, DNA methylation, and secretion systems. Rough variants showed higher abundance of EPS-associated proteins, the BpsI3/R3 quorum-sensing system, and the global regulator ScmR, whereas smooth variants displayed higher abundances of proteins belonging to type III/VI secretion and siderophore biosynthesis pathways. These findings suggest that CMV is driven by phase variation and regulatory mechanisms rather than punctual genomic modifications. Our study underscores the limitations of phenotype or genome-based classification alone in the context of CMV and highlights the value of integrated multi-omics approaches to uncover CMV-associated biomarkers, with potential applications in diagnostics and the development of targeted therapies against persistent and drug-resistant Bp infections.IMPORTANCEBurkholderia pseudomallei (Bp), the causative agent of melioidosis, is endemic to Australia, Asia, Africa, and the Americas. It predominantly affects Indigenous populations and individuals suffering from diabetes, chronic lung or kidney disease, or alcoholism. Bp is known for its exceptional genomic and phenotypic plasticity, enabling rapid adaptation to diverse environments. This adaptability is reflected by colony morphotype variation (CMV), including reversible phase variation between smooth and rough colonies. In this study, we report rough and smooth colonies from clinical samples and emphasize the importance of characterizing CMV through multi-omics approaches rather than relying solely on genomics and phenotypic traits. By integrating genomic, phenotypic, and proteomic data, we identified that a limited number of mutations, including one in a putative regulatory element, likely drive major molecular changes between morphotypes. These affect the expression of quorum-sensing systems, the transcriptional regulator ScmR, DNA methyltransferase, and virulence-associated genes.

Long-read RNA sequencing has been broadly utilized to examine the diversity of transcriptomes, understand differential expression and discover novel transcript isoforms. One of the major limitations of whole transcriptome sequencing is the difficulty in obtaining sufficient depth for low abundant transcripts. Methods which address this are either difficult to scale or customize: long- range PCR is customizable but difficult to scale beyond a few targets; probe hybridization panels are suited for scaling but require substantial investment to customize. In this study, we adopted RNA-guided CRISPR-Cas9 nuclease-based enrichment to target specific human and SARS-CoV-2 transcripts followed by long-read sequencing, utilizing minimal number of guide RNAs per target isoform. Our findings demonstrate that the CRISPR-Cas system is a highly effective method for customizable long-read sequencing of target transcripts while maintaining the accuracy of relative gene expression levels. The results highlight a valuable method for future research on transcript enrichment for isoform identification and low abundance transcript detection in infectious disease diagnosis. ### Competing Interest Statement LC has received travel funding and research funding from Oxford Nanopore Technologies unrelated to this manuscript.

BACKGROUND:Streptococcus dysgalactiae subspecies equisimilis (SDSE) is closely related to Streptococcus pyogenes, with overlapping disease manifestations. We compared the clinical and genomic epidemiology of invasive SDSE with invasive S pyogenes across different settings in Australia and phylogenetically contextualised the SDSE sequences within a global cohort of genomes. METHODS:In this retrospective clinical and genomic epidemiology study, cases of invasive SDSE isolated from normally sterile sites were identified and whole-genome sequenced across five hospital networks in temperate southeast Australia (Melbourne and Sydney) and the tropical Top End of the Northern Territory. SDSE disease incidence, case demographics, clinical outcomes, and longitudinal lineage dynamics were compared between southeast Australia and the Top End and to co-collected invasive S pyogenes cases in each region. SDSE genomes and lineages were also contextualised within 1166 global SDSE sequences. Genomic transmission clusters (not necessarily direct transmission) were inferred between isolates from different individuals by single-linkage clustering at a single nucleotide polymorphism threshold of less than or equal to seven for SDSE and less than or equal to five for S pyogenes based on previous transmission analyses. FINDINGS:Between Jan 1, 2011, and Feb 28, 2023, there were 693 invasive SDSE cases and 995 invasive S pyogenes cases. Invasive SDSE occurred almost exclusively in adults. The overall invasive SDSE incidence in southeast Australia was similar to invasive S pyogenes (incidence rate ratio [IRR] 1·15, 95% CI 0·91-1·46; p=0·26) and increased over the study period (IRR 1·06 per year, 95% CI 1·05-1·08; p<0·0001) from 1·30 cases per 10 000 admissions in 2011 to 3·72 cases per 10 000 admissions in the first 2 months of 2023 (95% CI 2·13-6·07). In southeast Australia, where stringent COVID-19 non-pharmaceutical interventions (NPIs) were implemented between 2020 and 2021, the SDSE incidence plateaued during 2020-21 but did not significantly decline (IRR 1·09 compared with 2017-19, 95% CI 0·88-1·35; p=0·47). By contrast, S pyogenes incidence substantially declined in 2020-21 in southeast Australia (IRR 0·35 compared to 2017-19, 95% CI 0·22-0·52; p=0·017). In the Top End, SDSE incidence was lower than S pyogenes (IRR 0·24, 95% CI 0·19-0·31; p<0·0001). However, crude incidence remained higher than southeast Australia (crude IRR 1·24, 95% CI 1·07-1·42; p=0·0037) and disproportionately affected First Nations Australians in the Top End compared with non-First Nations individuals (IRR 3·36, 95% CI 2·33-4·85; p<0·0001). Comparing 2020-21 with 2017-19, there was no decline in SDSE (IRR 1·27, 95% CI 0·73-2·24; p=0·45) or S pyogenes (IRR 0·97, 95% CI 0·80-1·18; p=0·81) incidence in the Top End, which did not implement prolonged stringent COVID-19 NPIs. Analysing the available genomes of invasive cases and in lineages for which more than or equal to five invasive cases occurred, only 24 (6%) of 384 SDSE cases were assigned to genomic transmission clusters, compared with 271 (52%) of 524 S pyogenes cases. An stG62647 lineage encompassed 113 (26%) of 436 sequenced SDSE genomes. Analysis of available SDSE sequences from Australia, western Europe, and North America inferred concurrent international expansion of the stG62647 lineage in all three regions between 1990 and 2005. INTERPRETATION:We identified a substantial burden of invasive SDSE, dominated by the emergent stG62647 lineage. The contrasting epidemiology between species in the different Australian regions, during COVID-19 NPIs, and genomic infection patterns indicates transmission dynamic, pathogen population, and host-pathogen interaction differences between SDSE and S pyogenes and indicates implications for disease control measures. FUNDING:Australian National Health and Medical Research Council.

Background and objectives The increasing occurrence of MRSA clinical isolates harbouring reduced susceptibility to mainstay antibiotics has escalated the use of second and last line antibiotics. Hence, it is critical to evaluate the likelihood of MRSA developing clinical resistance to these antibiotics. Our study sought to characterize the development of resistance to vancomycin (VAN), daptomycin (DAP) and linezolid (LZD) in MRSA ATCC 43300 in vitro and further determine the mechanisms underpinning resistance. Methods MRSA was exposed to increasing concentrations of VAN, DAP and LZD for 20 days, with eight replicates for each antibiotic conducted in parallel. The resulting day 20 (D20) isolates were subjected to antimicrobial susceptibility testing, whole genome sequencing, autolysis assays, and growth curves to determine bacterial fitness. Results Exposure to VAN or LZD for 20 days resulted in a subtle 2-fold increase in the MIC, whereas DAP exposure yielded DAP-non-susceptible isolates with up to 16-fold MIC increase. The MIC increase was accompanied by variable changes in relative fitness and reduced resistance to autolysis in some isolates. D20 isolates harboured mutations in genes commonly associated with resistance to the respective antibiotics (e.g. walK for VAN, mprF and rpoB for DAP, rplC for LZD), along with several previously unreported variants. Introduction of key mutations to these identified genes in the parental strain via allelic exchange confirmed their role in the development of resistance. Conclusions In vitro selection against VAN, DAP or LZD resulted in the acquisition of mutations similar to those correlated with clinical resistance, including the associated phenotypic alterations.

SYNOPSISBackgroundThe increasing occurrence of MRSA clinical isolates harbouring reduced susceptibility to mainstay antibiotics has escalated the use of second and last line antibiotics. Hence, it is critical to evaluate the likelihood of MRSA developing clinical resistance to these antibiotics.ObjectivesOur study sought to identify the rate in which MRSA develop resistance to vancomycin, daptomycin and linezolidin vitroand further determine the mechanisms underpinning resistance.MethodsMRSA was exposed to increasing concentrations of vancomycin, daptomycin, and linezolid for 20 days, with eight replicates for each antibiotic conducted in parallel. The resulting day 20 (D20) isolates were subjected to antimicrobial susceptibility testing, whole genome sequencing, autolysis assays, and growth curves to determine bacterial fitness.ResultsExposure to vancomycin or linezolid for 20 days resulted in a subtle two-fold increase in the MIC, whereas daptomycin exposure yielded daptomycin-nonsusceptible isolates with up to 16-fold MIC increase. The MIC increase was accompanied by variable changes in relative fitness and reduced resistance to autolysis in some isolates. D20 isolates harboured mutations in genes commonly associated with resistance to the respective antibiotics (e.g.walKfor vancomycin,mprFandrpoBfor daptomycin,rplCfor linezolid), along with several previously unreported variants. Introduction of key mutations to these identified genes in the parental strain via allelic exchange confirmed their role in the development of resistance.ConclusionsIn vitroselection against vancomycin, daptomycin, or linezolid resulted in the acquisition of mutations similar to those correlated with clinical resistance, including the associated phenotypic alterations.

Background: The increasing occurrence of MRSA clinical isolates harbouring reduced susceptibility to mainstay antibiotics has escalated the use of second and last line antibiotics. Hence, it is critical to evaluate the likelihood of MRSA developing clinical resistance to these antibiotics. Objectives: Our study sought to identify the rate in which MRSA develop resistance to vancomycin, daptomycin and linezolid in vitro and further determine the mechanisms underpinning resistance. Methods: MRSA was exposed to increasing concentrations of vancomycin, daptomycin, and linezolid for 20 days, with eight replicates for each antibiotic conducted in parallel. The resulting day 20 (D20) isolates were subjected to antimicrobial susceptibility testing, whole genome sequencing, autolysis assays, and growth curves to determine bacterial fitness. Results: Exposure to vancomycin or linezolid for 20 days resulted in a subtle two-fold increase in the MIC, whereas daptomycin exposure yielded daptomycin-nonsusceptible isolates with up to 16-fold MIC increase. The MIC increase was accompanied by variable changes in relative fitness and reduced resistance to autolysis in some isolates. D20 isolates harboured mutations in genes commonly associated with resistance to the respective antibiotics (e.g. walK for vancomycin, mprF and rpoB for daptomycin, rplC for linezolid), along with several previously unreported variants. Introduction of key mutations to these identified genes in the parental strain via allelic exchange confirmed their role in the development of resistance. Conclusions: In vitro selection against vancomycin, daptomycin, or linezolid resulted in the acquisition of mutations similar to those correlated with clinical resistance, including the associated phenotypic alterations. ### Competing Interest Statement M. A. C. currently holds a fractional Professorial Research Fellow appointment at the University of Queensland with his remaining time as CEO of Inflazome Ltd, a company with headquarters in Dublin, Ireland that is developing drugs to address clinical unmet needs in inflammatory disease by targeting the inflammasome.

Continuous assessment of the impact of SARS-CoV-2 on the host at the cell-type level is crucial for understanding key mechanisms involved in host defense responses to viral infection. We investigated host response to ancestral-strain and Alpha-variant SARS-CoV-2 infections within air-liquid-interface human nasal epithelial cells from younger adults (26-32 Y) and older children (12-14 Y) using single-cell RNA-sequencing. Ciliated and secretory-ciliated cells formed the majority of highly infected cell-types, with the latter derived from ciliated lineages. Strong innate immune responses were observed across lowly infected and uninfected bystander cells and heightened in Alpha-infection. Alpha highly infected cells showed increased expression of protein-refolding genes compared with ancestral-strain-infected cells in children. Furthermore, oxidative phosphorylation-related genes were down-regulated in bystander cells versus infected and mock-control cells, underscoring the importance of these biological functions for viral replication. Overall, this study highlights the complexity of cell-type-, age- and viral strain-dependent host epithelial responses to SARS-CoV-2.

A new variant of Streptococcus pyogenes serotype M1 (designated 'M1UK') has been reported in the United Kingdom, linked with seasonal scarlet fever surges, marked increase in invasive infections, and exhibiting enhanced expression of the superantigen SpeA. The progenitor S. pyogenes 'M1global' and M1UK clones can be differentiated by 27 SNPs and 4 indels, yet the mechanism for speA upregulation is unknown. Here we investigate the previously unappreciated expansion of M1UK in Australia, now isolated from the majority of serious infections caused by serotype M1 S. pyogenes. M1UK sub-lineages circulating in Australia also contain a novel toxin repertoire associated with epidemic scarlet fever causing S. pyogenes in Asia. A single SNP in the 5' transcriptional leader sequence of the transfer-messenger RNA gene ssrA drives enhanced SpeA superantigen expression as a result of ssrA terminator read-through in the M1UK lineage. This represents a previously unappreciated mechanism of toxin expression and urges enhanced international surveillance.

Motivation; PE/PPE proteins, highly abundant in the Mycobacterium genome, play a vital role in virulence and immune modulation. Understanding their functions is key to comprehending the internal mechanisms of Mycobacterium. However, a lack of dedicated resources has limited research into PE/PPE proteins. Results: Addressing this gap, we introduce MERITS, a comprehensive 3D structure database specifically designed for PE/PPE proteins. MERITS hosts 22,353 non-redundant PE/PPE proteins, encompassing details like physicochemical properties, subcellular localisation, post-translational modification sites, protein functions, and measures of antigenicity, toxicity, and allergenicity. MERITS also includes data on their secondary and tertiary structure, along with other relevant biological information. MERITS is designed to be user-friendly, offering interactive search and data browsing Features to aid researchers in exploring the potential functions of PE/PPE proteins. MERITS is expected to become a crucial resource in the field, aiding in developing new diagnostics and vaccines by elucidating the sequence-structure-functional relationships of PE/PPE proteins.

The diversity of COVID-19 disease in otherwise healthy people, from seemingly asymptomatic infection to severe life-threatening disease, is not clearly understood. We passaged a naturally occurring near-ancestral SARS-CoV-2 variant, capable of infecting wild-type mice, and identified viral genomic mutations coinciding with the acquisition of severe disease in young adult mice and lethality in aged animals. Transcriptomic analysis of lung tissues from mice with severe disease elucidated a host antiviral response dominated mainly by interferon and IL-6 pathway activation in young mice, while in aged animals, a fatal outcome was dominated by TNF and TGF-β signaling. Congruent with our pathway analysis, we showed that young TNF-deficient mice had mild disease compared to controls and aged TNF-deficient animals were more likely to survive infection. Emerging clinical correlates of disease are consistent with our preclinical studies, and our model may provide value in defining aberrant host responses that are causative of severe COVID-19.

Summary Assessing the impact of SARS-CoV-2 variants on the host is crucial with continuous emergence of new variants. We employed single-cell sequencing to investigate host transcriptomic response to ancestral and Alpha-strain SARS-CoV-2 infections within air-liquid-interface human nasal epithelial cells from adults and adolescents. Strong innate immune responses were observed across lowly-infected and bystander cell-types, and heightened in Alpha-infection. Contrastingly, the innate immune response of highly-infected cells was like mock-control cells. Alpha highly-infected cells showed increased expression of protein refolding genes compared with ancestral-strain-infected adolescent cells. Oxidative phosphorylation- and translation-related genes were down-regulated in bystander cells versus infected and mock-control cells, suggesting that the down-regulation is protective and up-regulation supports viral activity. Infected adult cells revealed up-regulation of these pathways compared with infected adolescents, implying enhanced pro-viral states in infected adults. Overall, this highlights the complexity of cell-type-, age- and viral-strain-dependent host epithelial responses to SARS-CoV-2 and the value of air-liquid-interface cultures.

Contents:nanopolish: direct RNA datasets analysed by nanopolish - raw_data: merged data between nanopolish and npTranscript results - analysis: testing changes in poly(A) length between control and infected using Wilcoxon's test of ranks with log-transformed data of Calu-3 48hpi direct RNA dataset- scripts: example scripts for testing changes in poly(A) length between control and infected using Wilcoxon's test of ranks, and also with log-transformed datatailfindr: direct cDNA datasets analysed by tailfindr - scripts: example scripts for poly(A) vs poly(T) correlation, mixed_model and generating raincloud plots for poly(A) length- raw_data: merged data between tailfindr and npTranscript results- analysis: tailfindr outputs and mixed model analysis outputs where files with names ending with 'int' shows the coefficient (changes in log tail length relative to control (intercept), and files with names ending with 'log_nofilter' shows the padj values for changes between control and infected poly(A) lengths.nanopolish_vs_tailfindr:- raw_data: median poly(A/T) lengths from nanopolish and tailfindr (before log transformation) - analysis: Spearman's correlation tests between nanopolish poly(A) and tailfindr poly(A/T) median lengths- scripts: example scripts for generating Spearman's correlation tests between nanopolish poly(A) and tailfindr poly(A/T) median lengthstailfindr_scripts: -tailfindr-related scripts onlytailfindr_calu_48hpi_analysis: - analysis results files from calu 48 hpi dataset including 'meancalc' - working sheet for backtransforming and calculating difference in log-tranformed poly(A) lengths between control and infected - control: intercepts for control - infected: coefficients for infected - mt: mitochondrial genes - nonmt: non-mitochondrial genes

The genome of Mycobacterium tuberculosis contains a relatively high percentage (10%) of genes that are poorly characterised because of their highly repetitive nature and high GC content. Some of these genes encode proteins of the PE/PPE family, which are thought to be involved in host-pathogen interactions, virulence, and disease pathogenicity. Members of this family are genetically divergent and challenging to both identify and classify using conventional computational tools. Thus, advanced in silico methods are needed to identify proteins of this family for subsequent functional annotation efficiently. In this study, we developed the first deep learning-based approach, termed Digerati, for the rapid and accurate identification of PE and PPE family proteins. Digerati was built upon a multipath parallel hybrid deep learning framework, which equips multi-layer convolutional neural networks with bidirectional, long short-term memory, equipped with a self-attention module to effectively learn the higher-order feature representations of PE/PPE proteins. Empirical studies demonstrated that Digerati achieved a significantly better performance (∼18-20%) than alignment-based approaches, including BLASTP, PHMMER, and HHsuite, in both prediction accuracy and speed. Digerati is anticipated to facilitate community-wide efforts to conduct high-throughput identification and analysis of PE/PPE family members. The webserver and source codes of Digerati are publicly available at http://web.unimelb-bioinfortools.cloud.edu.au/Digerati/.

Better methods to interrogate host-pathogen interactions during Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infections are imperative to help understand and prevent this disease. Here we implemented RNA-sequencing (RNA-seq) using Oxford Nanopore Technologies (ONT) long-reads to measure differential host gene expression, transcript polyadenylation and isoform usage within various epithelial cell lines permissive and non-permissive for SARS-CoV-2 infection. SARS-CoV-2-infected and mock-infected Vero (African green monkey kidney epithelial cells), Calu-3 (human lung adenocarcinoma epithelial cells), Caco-2 (human colorectal adenocarcinoma epithelial cells) and A549 (human lung carcinoma epithelial cells) were analyzed over time (0, 2, 24, 48 hours). Differential polyadenylation was found to occur in both infected Calu-3 and Vero cells during a late time point (48 hpi), with Gene Ontology (GO) terms such as viral transcription and translation shown to be significantly enriched in Calu-3 data. Poly(A) tails showed increased lengths in the majority of the differentially polyadenylated transcripts in Calu-3 and Vero cell lines (up to ~101 nt in mean poly(A) length, padj = 0.029). Of these genes, ribosomal protein genes such as RPS4X and RPS6 also showed downregulation in expression levels, suggesting the importance of ribosomal protein genes during infection. Furthermore, differential transcript usage was identified in Caco-2, Calu-3 and Vero cells, including transcripts of genes such as GSDMB and KPNA2 , which have previously been implicated in SARS-CoV-2 infections. Overall, these results highlight the potential role of differential polyadenylation and transcript usage in host immune response or viral manipulation of host mechanisms during infection, and therefore, showcase the value of long-read sequencing in identifying less-explored host responses to disease.

Mycobacterium tuberculosis genome comprises approximately 10% of two families of poorly characterised genes due to their high GC content and highly repetitive nature. The largest sub-group, the proline-glutamic acid polymorphic guanine-cytosine-rich sequence (PE_PGRS) family, is thought to be involved in host response and disease pathogenicity. Due to their high genetic variability and complexity of analysis, they are typically disregarded for further research in genomic studies. There are currently limited online resources and homology computational tools that can identify and analyse PE_PGRS proteins. In addition, they are computational-intensive and time-consuming, and lack sensitivity. Therefore, computational methods that can rapidly and accurately identify PE_PGRS proteins are valuable to facilitate the functional elucidation of the PE_PGRS family proteins. In this study, we developed the first machine learning-based bioinformatics approach, termed PEPPER, to allow users to identify PE_PGRS proteins rapidly and accurately. PEPPER was built upon a comprehensive evaluation of 13 popular machine learning algorithms with various sequence and physicochemical features. Empirical studies demonstrated that PEPPER achieved significantly better performance than alignment-based approaches, BLASTP and PHMMER, in both prediction accuracy and speed. PEPPER is anticipated to facilitate community-wide efforts to conduct high-throughput identification and analysis of PE_PGRS proteins.

Drug-resistant Gram-positive bacterial infections are still a substantial burden on the public health system, with two bacteria ( Staphylococcus aureus and Streptococcus pneumoniae ) accounting for over 1.5 million drug-resistant infections in the United States alone in 2017. In 2019, 250,000 deaths were attributed to these pathogens globally. We have developed a preclinical glycopeptide antibiotic, MCC5145, that has excellent potency (MIC 90 ≤ 0.06 μg/ml) against hundreds of isolates of methicillin-resistant S. aureus (MRSA) and other Gram-positive bacteria, with a greater than 1000-fold margin over mammalian cell cytotoxicity values. The antibiotic has therapeutic in vivo efficacy when dosed subcutaneously in multiple murine models of established bacterial infections, including thigh infection with MRSA and blood septicemia with S. pneumoniae , as well as when dosed orally in an antibiotic-induced Clostridioides difficile infection model. MCC5145 exhibited reduced nephrotoxicity at microbiologically active doses in mice compared to vancomycin. MCC5145 also showed improved activity against biofilms compared to vancomycin, both in vitro and in vivo, and a low propensity to select for drug resistance. Characterization of drug action using a transposon library bioinformatic platform showed a mechanistic distinction from other glycopeptide antibiotics.

Genomic neighbor typing enables heuristic inference of bacterial lineages and phenotypes from nanopore sequencing data. However, small reference databases may not be sufficiently representative of the diversity of lineages and genotypes present in a collection of isolates. In this study, we explore the use of genomic neighbor typing for surveillance of community-associated Staphylococcus aureus outbreaks in Papua New Guinea (PNG) and Far North Queensland, Australia (FNQ). We developed Sketchy , an implementation of genomic neighbor typing that queries exhaustive whole genome reference databases using MinHash. Evaluations were conducted using nanopore read simulations and six species-wide reference sketches (4832 - 47616 genomes), as well as two S. aureus outbreak data sets sequenced at low depth using a sequential multiplex library protocol on the MinION (n = 160, with matching Illumina data). Heuristic inference of lineages and antimicrobial resistance profiles allowed us to conduct multiplex genotyping in situ at the Papua New Guinea Institute of Medical Research in Goroka, on low-throughput Flongle adapters and using multiple successive libraries on the same MinION flow cell (n = 24 - 48). Comparison to phylogenetically informed genomic neighbor typing with RASE on the dominant outbreak sequence type suggests slightly better performance at predicting lineage-scale genotypes using large sketch sizes, but inferior performance in resolving clade-specific genotypes (methicillin resistance). Sketchy can be used for large-scale bacterial outbreak surveillance and in challenging sequencing scenarios, but improvements to clade-specific genotype inference are needed for diagnostic applications. Sketchy is available open-source at: https://github.com/esteinig/sketchy

Summary SARS-CoV-2 uses subgenomic (sg)RNA to produce viral proteins for replication and immune evasion. We applied long-read RNA and cDNA sequencing to in vitro human and primate infection models to study transcriptional dynamics. Transcription-regulating sequence (TRS)-dependent sgRNA was upregulated earlier in infection than TRS-independent sgRNA. An abundant class of TRS-independent sgRNA consisting of a portion of ORF1ab containing nsp1 joined to ORF10 and 3’UTR was upregulated at 48 hours post infection in human cell lines. We identified double-junction sgRNA containing both TRS-dependent and independent junctions. We found multiple sites at which the SARS-CoV-2 genome is consistently more modified than sgRNA, and that sgRNA modifications are stable across transcript clusters, host cells and time since infection. Our work highlights the dynamic nature of the SARS-CoV-2 transcriptome during its replication cycle. Our results are available via an interactive web-app at http://coinlab.mdhs.unimelb.edu.au/ .