Whole genome sequencing of human Borrelia burgdorferi isolates reveals linked blocks of accessory genome elements located on plasmids and associated with human dissemination.

Lyme disease is the most common vector-borne disease in North America and Europe. The clinical manifestations of Lyme disease vary based on the genospecies of the infecting Borrelia burgdorferi spirochete, but the microbial genetic elements underlying these associations are not known. Here, we report the whole genome sequence (WGS) and analysis of 299 B. burgdorferi (Bb) isolates derived from patients in the Eastern and Midwestern US and Central Europe. We develop a WGS-based classification of Bb isolates, confirm and extend the findings of previous single- and multi-locus typing systems, define the plasmid profiles of human-infectious Bb isolates, annotate the core and strain-variable surface lipoproteome, and identify loci associated with disseminated infection. A core genome consisting of ~900 open reading frames and a core set of plasmids consisting of lp17, lp25, lp36, lp28-3, lp28-4, lp54, and cp26 are found in nearly all isolates. Strain-variable (accessory) plasmids and genes correlate strongly with phylogeny. Using genetic association study methods, we identify an accessory genome signature associated with dissemination in humans and define the individual plasmids and genes that make up this signature. Strains within the RST1/WGS A subgroup, particularly a subset marked by the OspC type A genotype, have increased rates of dissemination in humans. OspC type A strains possess a unique set of strongly linked genetic elements including the presence of lp56 and lp28-1 plasmids and a cluster of genes that may contribute to their enhanced virulence compared to other genotypes. These features of OspC type A strains reflect a broader paradigm across Bb isolates, in which near-clonal genotypes are defined by strain-specific clusters of linked genetic elements, particularly those encoding surface-exposed lipoproteins. These clusters of genes are maintained by strain-specific patterns of plasmid occupancy and are associated with the probability of invasive infection.

Whole genome sequencing of Borrelia burgdorferi isolates reveals linked clusters of plasmid-borne accessory genome elements associated with virulence.

Lyme disease is the most common vector-borne disease in North America and Europe. The clinical manifestations of Lyme disease vary based on the genospecies of the infecting Borrelia burgdorferi spirochete, but the microbial genetic elements underlying these associations are not known. Here, we report the whole genome sequence (WGS) and analysis of 299 patient-derived B. burgdorferi sensu stricto ( Bbss ) isolates from patients in the Eastern and Midwestern US and Central Europe. We develop a WGS-based classification of Bbss isolates, confirm and extend the findings of previous single- and multi-locus typing systems, define the plasmid profiles of human-infectious Bbss isolates, annotate the core and strain-variable surface lipoproteome, and identify loci associated with disseminated infection. A core genome consisting of ∻800 open reading frames and a core set of plasmids consisting of lp17, lp25, lp36, lp28-3, lp28-4, lp54, and cp26 are found in nearly all isolates. Strain-variable (accessory) plasmids and genes correlate strongly with phylogeny. Using genetic association study methods, we identify an accessory genome signature associated with dissemination and define the individual plasmids and genes that make up this signature. Strains within the RST1/WGS A subgroup, particularly a subset marked by the OspC type A genotype, are associated with increased rates of dissemination. OspC type A strains possess a unique constellation of strongly linked genetic changes including the presence of lp56 and lp28-1 plasmids and a cluster of genes that may contribute to their enhanced virulence compared to other genotypes. The patterns of OspC type A strains typify a broader paradigm across Bbss isolates, in which genetic structure is defined by correlated groups of strain-variable genes located predominantly on plasmids, particularly for expression of surface-exposed lipoproteins. These clusters of genes are inherited in blocks through strain-specific patterns of plasmid occupancy and are associated with the probability of invasive infection.

Delta-Omicron recombinant escapes therapeutic antibody neutralization.

The emergence of recombinant viruses is a threat to public health, as recombination may integrate variant-specific features that together result in escape from treatment or immunity. The selective advantages of recombinant SARS-CoV-2 isolates over their parental lineages remain unknown. We identified a Delta-Omicron (AY.45-BA.1) recombinant in an immunosuppressed transplant recipient treated with monoclonal antibody Sotrovimab. The single recombination breakpoint is located in the spike N-terminal domain adjacent to the Sotrovimab binding site. While Delta and BA.1 are sensitive to Sotrovimab neutralization, the Delta-Omicron recombinant is highly resistant. To our knowledge, this is the first described instance of recombination between circulating SARS-CoV-2 variants as a functional mechanism of resistance to treatment and immune escape.

Functional control and metabolic integration of stolen organelles in a photosynthetic ciliate.

Stealing prey plastids for metabolic gain is a common phenomenon among protists within aquatic ecosystems.1 Ciliates of the Mesodinium rubrum species complex are unique in that they also steal a transcriptionally active but non-dividing prey nucleus, the kleptokaryon, from certain cryptophytes.2 The kleptokaryon enables full control and replication of kleptoplastids but has a half-life of about 10 days.2 Once the kleptokaryon is lost, the ciliate experiences a slow loss of photosynthetic metabolism and eventually death.2,3,4 This transient ability to function phototrophically allows M. rubrum to form productive blooms in coastal waters.5,6,7,8 Here, we show, using multi-omics approaches, that an Antarctic strain of the ciliate not only depends on stolen Geminigera cryophila organelles for photosynthesis but also for anabolic synthesis of fatty acids, amino acids, and other essential macromolecules. Transcription of diverse pathways was higher in the kleptokaryon than that in G. cryophila, and many increased in higher light. Proteins of major biosynthetic pathways were found in greater numbers in the kleptokaryon relative to M. rubrum, implying anabolic dependency on foreign metabolism. We show that despite losing transcriptional control of the kleptokaryon, M. rubrum regulates kleptoplastid pigments with changing light, implying an important role for post-transcriptional control. These findings demonstrate that the integration of foreign organelles and their gene and protein expression, energy metabolism, and anabolism occur in the absence of a stable endosymbiotic association. Our results shed light on potential events early in the process of complex plastid acquisition and broaden our understanding of symbiogenesis.

Mycobacterial phosphodiesterase Rv0805 is a virulence determinant and its cyclic nucleotide hydrolytic activity is required for propionate detoxification.

Cyclic AMP (cAMP) signaling is essential to Mycobacterium tuberculosis (Mtb) pathogenesis. However, the roles of phosphodiesterases (PDEs) Rv0805, and the recently identified Rv1339, in cAMP homeostasis and Mtb biology are unclear. We found that Rv0805 modulates Mtb growth within mice, macrophages and on host-associated carbon sources. Mycobacterium bovis BCG grown on a combination of propionate and glycerol as carbon sources showed high levels of cAMP and had a strict requirement for Rv0805 cNMP hydrolytic activity. Supplementation with vitamin B12 or spontaneous genetic mutations in the pta-ackA operon restored the growth of BCGΔRv0805 and eliminated propionate-associated cAMP increases. Surprisingly, reduction of total cAMP levels by ectopic expression of Rv1339 restored only 20% of growth, while Rv0805 complementation fully restored growth despite a smaller effect on total cAMP levels. Deletion of an Rv0805 localization domain also reduced BCG growth in the presence of propionate and glycerol. We propose that localized Rv0805 cAMP hydrolysis modulates activity of a specialized pathway associated with propionate metabolism, while Rv1339 has a broader role in cAMP homeostasis. Future studies will address the biological roles of Rv0805 and Rv1339, including their impacts on metabolism, cAMP signaling and Mtb pathogenesis.

Cascading effects of prey identity on gene expression in a kleptoplastidic ciliate.

Kleptoplastidic, or chloroplast stealing, lineages transiently retain functional photosynthetic machinery from algal prey. This machinery, and its photosynthetic outputs, must be integrated into the host's metabolism, but the details of this integration are poorly understood. Here, we study this metabolic integration in the ciliate Mesodinium chamaeleon, a coastal marine species capable of retaining chloroplasts from at least six distinct genera of cryptophyte algae. To assess the effects of feeding history on ciliate physiology and gene expression, we acclimated M. chamaeleon to four different types of prey and contrasted well-fed and starved treatments. Consistent with previous physiological work on the ciliate, we found that starved ciliates had lower chlorophyll content, photosynthetic rates, and growth rates than their well-fed counterparts. However, ciliate gene expression mirrored prey phylogenetic relationships rather than physiological status, suggesting that, even as M. chamaeleon cells were starved of prey, their overarching regulatory systems remained tuned to the prey type to which they had been acclimated. Collectively, our results indicate a surprising degree of prey-specific host transcriptional adjustments, implying varied integration of prey metabolic potential into many aspects of ciliate physiology.

Combining genomic and epidemiological data to compare the transmissibility of SARS-CoV-2 variants Alpha and Iota.

SARS-CoV-2 variants shaped the second year of the COVID-19 pandemic and the discourse around effective control measures. Evaluating the threat posed by a new variant is essential for adapting response efforts when community transmission is detected. In this study, we compare the dynamics of two variants, Alpha and Iota, by integrating genomic surveillance data to estimate the effective reproduction number (Rt) of the variants. We use Connecticut, United States, in which Alpha and Iota co-circulated in 2021. We find that the Rt of these variants were up to 50% larger than that of other variants. We then use phylogeography to show that while both variants were introduced into Connecticut at comparable frequencies, clades that resulted from introductions of Alpha were larger than those resulting from Iota introductions. By monitoring the dynamics of individual variants throughout our study period, we demonstrate the importance of routine surveillance in the response to COVID-19.

Detection of SARS-CoV-2 Variants by Multiple Diagnostic Assays Second RESUBMISSION JCV-D-21-00675R2.

INTRODUCTION: The emergence of SARS-CoV-2 Variants of Concern (VOC) and Variants Being Monitored (VBM) have presented additional clinical and public health concerns regarding potential virus transmissibility, disease severity, and immune evasion. It is imperative that diagnostic assays can detect all such variants, and since commercial oligo sequences are commonly not available, empirical testing may be necessary to confirm this. To confirm the sensitivity of the SARS-CoV-2 assays used at the Wadsworth Center for the detection of VOC and VBM, relevant specimens were selected from the specimen archive and tested in the various platforms. MATERIALS AND METHODS: Patient respiratory specimens submitted from clinal laboratories across the state were selected; three samples per variant were chosen to account for inter assay and variant reproducibility. The four molecular diagnostic platforms for SARS-CoV-2 currently in use at our facility were examined. RESULTS: A total of 64 specimens were tested, representing 2 VOC, 8 VBM and 4 other variants circulating in New York State. For certain samples, original Ct values provided by sample submitters were much higher, or lower, than those obtained from this study. The investigation of submitter testing platforms, with consideration of the assay's viral targets, confirmed the differences in Ct were not variant specific. CONCLUSIONS: It was demonstrated that the diagnostic methods investigated in this study detected all the variants tested. Because of the continual evolution of the virus, it is vital to monitor new variants as they emerge for the ability of molecular diagnostic methods to detect them with acceptable sensitivity.

Delta-Omicron recombinant escapes therapeutic antibody neutralization.

Background: The emergence of recombinant viruses is a threat to public health. Recombination of viral variants may combine variant-specific features that together catalyze viral escape from treatment or immunity. The selective advantages of recombinant SARS-CoV-2 isolates over their parental lineages remain unknown. Methods: Multi-method amplicon and metagenomic sequencing of a clinical swab and the in vitro grown virus allowed for high-confidence detection of a novel recombinant variant. Mutational, phylogeographic, and structural analyses determined features of the recombinant genome and spike protein. Neutralization assays using infectious as well as pseudotyped viruses and point mutants thereof defined the recombinant's sensitivity to a panel of monoclonal antibodies and sera from vaccinated and/or convalescent individuals. Results: A novel Delta-Omicron SARS-CoV-2 recombinant was identified in an unvaccinated, immunosuppressed kidney transplant recipient treated with monoclonal antibody Sotrovimab. The recombination breakpoint is located in the spike N-terminal domain, adjacent to the Sotrovimab quaternary binding site, and results in a 5'-Delta AY.45 and a 3'-Omicron BA.1 mosaic spike protein. Delta and BA.1 are sensitive to Sotrovimab neutralization, whereas the Delta-Omicron recombinant is highly resistant to Sotrovimab, both with and without the RBD resistance mutation E340D. Conclusions: Recombination between circulating SARS-CoV-2 variants can functionally contribute to immune escape. It is critical to validate phenotypes of mosaic viruses and monitor immunosuppressed COVID-19 patients treated with monoclonal antibodies for the selection of recombinant and immune escape variants. (Funded by NYU, the National Institutes of Health, and others).

Depletion of the Microbiota Has a Modest but Important Impact on the Fungal Burden of the Heart and Lungs during Early Systemic Candida auris Infection in Neutropenic Mice.

The progression and systemic pathobiology of C. auris in the absence of a microbiota have not been described. Here, we describe the influence of the microbiota during the first 5 days of C. auris infection in germ-free or antibiotic-depleted mice. Depletion of the bacterial microbiota in both germ-free and antibiotic-depleted models results in a modest but important increase in the early stages of C. auris infection. Particularly the heart and lungs, followed by the cecum, uterus, and stomach, of intravenously (i.v.) infected neutropenic mice showed significant fungal organ burden. Understanding disease progression and pathobiology of C. auris in individuals with a depleted microbiota could potentially help in the development of care protocols that incorporate supplementation or restoration of the microbiota before invasive procedures, such as transplantation surgeries.

Spatiotemporal Analyses of 2 Co-Circulating SARS-CoV-2 Variants, New York State, USA.

The emergence of novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants in late 2020 and early 2021 raised alarm worldwide because of their potential for increased transmissibility and immune evasion. Elucidating the evolutionary and epidemiologic dynamics among novel SARS-CoV-2 variants is essential for understanding the trajectory of the coronavirus disease pandemic. We describe the interplay between B.1.1.7 (Alpha) and B.1.526 (Iota) variants in New York State, USA, during December 2020-April 2021 through phylogeographic analyses, space-time scan statistics, and cartographic visualization. Our results indicate that B.1.526 probably evolved in New York City, where it was displaced as the dominant lineage by B.1.1.7 months after its initial appearance. In contrast, B.1.1.7 became dominant earlier in regions with fewer B.1.526 infections. These results suggest that B.1.526 might have delayed the initial spread of B.1.1.7 in New York City. Our combined spatiotemporal methodologies can help disentangle the complexities of shifting SARS-CoV-2 variant landscapes.

The PqsE-RhlR Interaction Regulates RhlR DNA Binding to Control Virulence Factor Production in Pseudomonas aeruginosa.

Pseudomonas aeruginosa is an opportunistic pathogen that causes disease in immunocompromised individuals and individuals with underlying pulmonary disorders. P. aeruginosa virulence is controlled by quorum sensing (QS), a bacterial cell-cell communication mechanism that underpins transitions between individual and group behaviors. In P. aeruginosa, the PqsE enzyme and the QS receptor RhlR directly interact to control the expression of genes involved in virulence. Here, we show that three surface-exposed arginine residues on PqsE comprise the site required for interaction with RhlR. We show that a noninteracting PqsE variant [PqsE(NI)] possesses catalytic activity, but is incapable of promoting virulence phenotypes, indicating that interaction with RhlR, and not catalysis, drives these PqsE-dependent behaviors. Biochemical characterization of the PqsE-RhlR interaction coupled with RNA-seq analyses demonstrates that the PqsE-RhlR complex increases the affinity of RhlR for DNA, enabling enhanced expression of genes encoding key virulence factors. These findings provide the mechanism for PqsE-dependent regulation of RhlR and identify a unique regulatory feature of P. aeruginosa QS and its connection to virulence. IMPORTANCE Bacteria use a cell-cell communication process called quorum sensing (QS) to orchestrate collective behaviors. QS relies on the group-wide detection of molecules called autoinducers (AI). QS is required for virulence in the human pathogen Pseudomonas aeruginosa, which can cause fatal infections in patients with underlying pulmonary disorders. In this study, we determine the molecular basis for the physical interaction between two virulence-driving QS components, PqsE and RhlR. We find that the ability of PqsE to bind RhlR correlates with virulence factor production. Since current antimicrobial therapies exacerbate the growing antibiotic resistance problem because they target bacterial growth, we suggest that the PqsE-RhlR interface discovered here represents a new candidate for targeting with small molecule inhibition. Therapeutics that disrupt the PqsE-RhlR interaction should suppress virulence. Targeting bacterial behaviors such as QS, rather than bacterial growth, represents an attractive alternative for exploration because such therapies could potentially minimize the development of resistance.

Phenotypic and genetic changes associated with the seesaw effect in MRSA strain N315 in a bioreactor model.

OBJECTIVES: Over the past decade, daptomycin treatment of methicillin-resistant Staphylococcus aureus (MRSA) infections has led to the emergence of daptomycin nonsusceptible (DAP-NS) MRSA strains and a subsequent interest in combinatorial antibiotic therapies. We investigated the phenotypic and genetic changes associated with the seesaw effect, which describes the correlation between daptomycin resistance and increased ß-lactam susceptibility in DAP-NS MRSA and the reverse phenomenon of DAP-NS strains acquiring renewed susceptibility to daptomycin after ß-lactam exposure. METHODS: A continuous bioreactor model was used to study the effects of incremental doses of daptomycin followed by oxacillin on MRSA strain N315. Minimum inhibitory concentrations for daptomycin and oxacillin were determined for the bioreactor-derived samples. Transmission electron microscopy and cytochrome C binding assays were used to measure cell wall thickness and cell membrane charge, respectively, in the bioreactor-derived samples. Whole-genome sequencing was used to identify mutations associated with the seesaw effect. RESULTS: Although daptomycin resistance conferred enhanced susceptibility to oxacillin, oxacillin treatment of DAP-NS strains was accompanied by a lowered minimum inhibitory concentration for daptomycin. Additionally, there was a reduction in relative positive cell surface charge and cell wall thickness. However, the mutations acquired in our DAP-NS populations were not accompanied by additional genomic changes after treatment with oxacillin, implicating alternative mechanisms for the seesaw effect. CONCLUSION: In this study, we successfully produced and characterized the seesaw effect in MRSA strain N315 in a unique bioreactor model.

Whole-Genome Sequencing of SARS-CoV-2: Assessment of the Ion Torrent AmpliSeq Panel and Comparison with the Illumina MiSeq ARTIC Protocol.

Fast and effective methods are needed for sequencing of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome to track genetic mutations and to identify new and emerging variants during the ongoing pandemic. The objectives were to assess the performance of the SARS-CoV-2 AmpliSeq research panel and S5 plug-in analysis tools for whole-genome sequencing analysis of SARS-CoV-2 and to compare the results with those obtained with the MiSeq-based ARTIC analysis pipeline, using metrics such as depth, coverage, and concordance of single-nucleotide variant (SNV) calls. A total of 191 clinical specimens and a single cultured isolate were extracted and sequenced with AmpliSeq technology and analysis tools. Of the 191 clinical specimens, 83 (with threshold cycle [CT] values of 15.58 to 32.54) were also sequenced using an Illumina MiSeq-based method with the ARTIC analysis pipeline, for direct comparison. A total of 176 of the 191 clinical specimens sequenced on the S5XL system and prepared using the SARS-CoV-2 research panel had nearly complete coverage (>98%) of the viral genome, with an average depth of 5,031×. Similar coverage levels (>98%) were observed for 81/83 primary specimens that were sequenced with both methods tested. The sample with the lowest viral load (CT value of 32.54) achieved 89% coverage using the MiSeq method and failed to sequence with the AmpliSeq method. Consensus sequences produced by each method were identical for 81/82 samples in areas of equal coverage, with a single difference present in one sample. The AmpliSeq approach is as effective as the Illumina-based method using ARTIC v3 amplification for sequencing SARS-CoV-2 directly from patient specimens across a range of viral loads (CT values of 15.56 to 32.54 [median, 22.18]). The AmpliSeq workflow is very easily automated with the Ion Chef and S5 instruments and requires less training and experience with next-generation sequencing sample preparation than the Illumina workflow.

Combining genomic and epidemiological data to compare the transmissibility of SARS-CoV-2 lineages.

Emerging SARS-CoV-2 variants have shaped the second year of the COVID-19 pandemic and the public health discourse around effective control measures. Evaluating the public health threat posed by a new variant is essential for appropriately adapting response efforts when community transmission is detected. However, this assessment requires that a true comparison can be made between the new variant and its predecessors because factors other than the virus genotype may influence spread and transmission. In this study, we develop a framework that integrates genomic surveillance data to estimate the relative effective reproduction number (R t ) of co-circulating lineages. We use Connecticut, a state in the northeastern United States in which the SARS-CoV-2 variants B.1.1.7 and B.1.526 co-circulated in early 2021, as a case study for implementing this framework. We find that the R t of B.1.1.7 was 6-10% larger than that of B.1.526 in Connecticut in the midst of a COVID-19 vaccination campaign. To assess the generalizability of this framework, we apply it to genomic surveillance data from New York City and observe the same trend. Finally, we use discrete phylogeography to demonstrate that while both variants were introduced into Connecticut at comparable frequencies, clades that resulted from introductions of B.1.1.7 were larger than those resulting from B.1.526 introductions. Our framework, which uses open-source methods requiring minimal computational resources, may be used to monitor near real-time variant dynamics in a myriad of settings.

RNA-Seq analysis of blood meal induced gene-expression changes in Aedes aegypti ovaries.

BACKGROUND: Transmission of pathogens by vector mosquitoes is intrinsically linked with mosquito's reproductive strategy because anautogenous mosquitoes require vertebrate blood to develop a batch of eggs. Each cycle of egg maturation is tightly linked with the intake of a fresh blood meal for most species. Mosquitoes that acquire pathogens during the first blood feeding can transmit the pathogens to susceptible hosts during subsequent blood feeding and also vertically to the next generation via infected eggs. Large-scale gene-expression changes occur following each blood meal in various tissues, including ovaries. Here we analyzed mosquito ovary transcriptome following a blood meal at three different time points to investigate blood-meal induced changes in gene expression in mosquito ovaries. RESULTS: We collected ovaries from Aedes aegypti that received a sugar meal or a blood meal on days 3, 10 and 20 post blood meal for transcriptome analysis. Over 4000 genes responded differentially following ingestion of a blood meal on day 3, and 660 and 780 genes on days 10 and 20, respectively. Proteins encoded by differentially expressed genes (DEGs) on day 3 include odorant binding proteins (OBPs), defense-specific proteins, and cytochrome P450 detoxification enzymes. In addition, we identified 580 long non-coding RNAs that are differentially expressed at three time points. Gene ontology analysis indicated that genes involved in peptidase activity, oxidoreductase activity, extracellular space, and hydrolase activity, among others were enriched on day 3. Although most of the DEGs returned to the nonsignificant level compared to the sugar-fed mosquito ovaries following oviposition on days 10 and 20, there remained differences in the gene expression pattern in sugar-fed and blood-fed mosquitoes. CONCLUSIONS: Enrichment of OBPs following blood meal ingestion suggests that these genes may have other functions besides being part of the olfactory system. The enrichment of immune-specific genes and cytochrome P450 genes indicates that ovaries become well prepared to protect their germ line from any pathogens that may accompany the blood meal or from environmental contamination during oviposition, and to deal with the detrimental effects of toxic metabolites.

Early introductions and transmission of SARS-CoV-2 variant B.1.1.7 in the United States.

The emergence and spread of SARS-CoV-2 lineage B.1.1.7, first detected in the United Kingdom, has become a global public health concern because of its increased transmissibility. Over 2,500 COVID-19 cases associated with this variant have been detected in the United States (US) since December 2020, but the extent of establishment is relatively unknown. Using travel, genomic, and diagnostic data, we highlight that the primary ports of entry for B.1.1.7 in the US were in New York, California, and Florida. Furthermore, we found evidence for many independent B.1.1.7 establishments starting in early December 2020, followed by interstate spread by the end of the month. Finally, we project that B.1.1.7 will be the dominant lineage in many states by mid- to late March. Thus, genomic surveillance for B.1.1.7 and other variants urgently needs to be enhanced to better inform the public health response.

Early introductions and community transmission of SARS-CoV-2 variant B.1.1.7 in the United States.

The emergence and spread of SARS-CoV-2 lineage B.1.1.7, first detected in the United Kingdom, has become a global public health concern because of its increased transmissibility. Over 2500 COVID-19 cases associated with this variant have been detected in the US since December 2020, but the extent of establishment is relatively unknown. Using travel, genomic, and diagnostic data, we highlight the primary ports of entry for B.1.1.7 in the US and locations of possible underreporting of B.1.1.7 cases. Furthermore, we found evidence for many independent B.1.1.7 establishments starting in early December 2020, followed by interstate spread by the end of the month. Finally, we project that B.1.1.7 will be the dominant lineage in many states by mid to late March. Thus, genomic surveillance for B.1.1.7 and other variants urgently needs to be enhanced to better inform the public health response.

Divergent Mutational Landscapes of Consensus and Minority Genotypes of West Nile Virus Demonstrate Host and Gene-Specific Evolutionary Pressures.

Our current understanding of the natural evolution of RNA viruses comes largely from consensus level genetic analyses which ignore the diverse mutant swarms that comprise within-host viral populations. The breadth and composition of viral mutant swarms impact viral fitness and adaptation, and the capacity for swarm plasticity is likely to be particularly important for arthropod-borne viruses (arboviruses) that cycle between taxonomically divergent hosts. Despite this, characterization of the relationship between the selective pressures and genetic signatures of the mutant swarm and consensus sequences is lacking. To clarify this, we analyzed previously generated whole genome, deep-sequencing data from 548 West Nile virus samples isolated from avian tissues or mosquitoes in New York State from 1999-2018. Both consensus level (interhost) and minority level (intrahost) nucleotide and amino acid sequences were analyzed, and diversity at each position was calculated across the genome in order to assess the relationship between minority and consensus sequences for individual genes and hosts. Our results indicate that consensus sequences are an inept representation of the overall genetic diversity. Unique host and gene-specific signatures and selective pressures were identified. These data demonstrate that an accurate and comprehensive understanding of arbovirus evolution and adaptation within and between hosts requires consideration of minority genotypes.

Small-Molecule Antibiotics Inhibiting tRNA-Regulated Gene Expression Is a Viable Strategy for Targeting Gram-Positive Bacteria.

Bacterial infections and the rise of antibiotic resistance, especially multidrug resistance, have generated a clear need for discovery of novel therapeutics. We demonstrated that a small-molecule drug, PKZ18, targets the T-box mechanism and inhibits bacterial growth. The T-box is a structurally conserved riboswitch-like gene regulator in the 5' untranslated region (UTR) of numerous essential genes of Gram-positive bacteria. T-boxes are stabilized by cognate, unacylated tRNA ligands, allowing the formation of an antiterminator hairpin in the mRNA that enables transcription of the gene. In the absence of an unacylated cognate tRNA, transcription is halted due to the formation of a thermodynamically more stable terminator hairpin. PKZ18 targets the site of the codon-anticodon interaction of the conserved stem I and reduces T-box-controlled gene expression. Here, we show that novel analogs of PKZ18 have improved MICs, bactericidal effects against methicillin-resistant Staphylococcus aureus (MRSA), and increased efficacy in nutrient-limiting conditions. The analogs have reduced cytotoxicity against eukaryotic cells compared to PKZ18. The PKZ18 analogs acted synergistically with aminoglycosides to significantly enhance the efficacy of the analogs and aminoglycosides, further increasing their therapeutic windows. RNA sequencing showed that the analog PKZ18-22 affects expression of 8 of 12 T-box controlled genes in a statistically significant manner, but not other 5'-UTR regulated genes in MRSA. Very low levels of resistance further support the existence of multiple T-box targets for PKZ18 analogs in the cell. Together, the multiple targets, low resistance, and synergy make PKZ18 analogs promising drugs for development and future clinical applications.