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.

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.

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.

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.

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.

Polycysteine-encoding leaderless short ORFs function as cysteine-responsive attenuators of operonic gene expression in mycobacteria.

Genome-wide transcriptomic analyses have revealed abundant expressed short open reading frames (ORFs) in bacteria. Whether these short ORFs, or the small proteins they encode, are functional remains an open question. One quarter of mycobacterial mRNAs are leaderless, beginning with a 5'-AUG or GUG initiation codon. Leaderless mRNAs often encode unannotated short ORFs as the first gene of a polycistronic transcript. Here, we show that polycysteine-encoding leaderless short ORFs function as cysteine-responsive attenuators of operonic gene expression. Detailed mutational analysis shows that one polycysteine short ORF controls expression of the downstream genes. Our data indicate that ribosomes stalled in the polycysteine tract block mRNA structures that otherwise sequester the ribosome-binding site of the 3'gene. We assessed endogenous proteomic responses to cysteine limitation in Mycobacterium smegmatis using mass spectrometry. Six cysteine metabolic loci having unannotated polycysteine-encoding leaderless short ORF architectures responded to cysteine limitation, revealing widespread cysteine-responsive attenuation in mycobacteria. Individual leaderless short ORFs confer independent operon-level control, while their shared dependence on cysteine ensures a collective response mediated by ribosome pausing. We propose the term ribulon to classify ribosome-directed regulons. Regulon-level coordination by ribosomes on sensory short ORFs illustrates one utility of the many unannotated short ORFs expressed in bacterial genomes.

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.

Direct cell-cell contact activates SigM to express the ESX-4 secretion system in Mycobacterium smegmatis.

Conjugal cell-cell contact between strains of Mycobacterium smegmatis induces the esxUT transcript, which encodes the putative primary substrates of the ESAT-6 secretion system 4 (ESX-4) secretion system. This recipient response was required for conjugal transfer of chromosomal DNA from the donor strain. Here we show that the extracytoplasmic σ factor, SigM, is a cell contact-dependent activator of ESX-4 expression and is required for conjugal transfer of DNA in the recipient strain. The SigM regulon includes genes outside the seven-gene core esx4 locus that we show are also required for conjugation, and we show that some of these SigM-induced proteins likely function through ESX-4. A fluorescent reporter revealed that SigM is specifically activated in recipient cells in direct contact with donor cells. Coculture RNA-seq experiments indicated that SigM regulon induction occurred early and before transconjugants are detected. This work supports a model wherein donor contact with the recipient cell surface inactivates the transmembrane anti-SigM, thereby releasing SigM. Free SigM induces an extended ESX-4 secretion system, resulting in changes that facilitate chromosomal transfer. The contact-dependent inactivation of an extracytoplasmic σ-factor that tightly controls ESX-4 activity suggests a mechanism dedicated to detect, and appropriately respond to, external stimuli from mycobacteria.

Origin of 3 Rabid Terrestrial Animals in Raccoon Rabies Virus-Free Zone, Long Island, New York, USA, 2016-2017.

During 2016-2017, three rabid terrestrial animals were discovered in the raccoon rabies virus-free zone of Long Island, New York, USA. Whole-genome sequencing and phylogenetic analyses revealed the likely origins of the viruses, enabling the rabies outbreak response (often costly and time-consuming) to be done less expensively and more efficiently.

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.

Impact of Candida auris Infection in a Neutropenic Murine Model.

Candida auris has become a global public health threat due to its multidrug resistance and persistence. Currently, there are limited murine models to study C. auris infection. Those models use a combination of cyclophosphamide and cortisone acetate, suppressing both innate and adaptive immunity. Here, we compare C. auris infection in two neutrophil-depleted murine models in which innate immunity is targeted using the monoclonal antibodies 1A8 and RB6-8C5.

Whole-Genome Single-Nucleotide Polymorphism (SNP) Analysis Applied Directly to Stool for Genotyping Shiga Toxin-Producing Escherichia coli: an Advanced Molecular Detection Method for Foodborne Disease Surveillance and Outbreak Tracking.

Whole-genome sequencing (WGS) of pathogens from pure culture provides unparalleled accuracy and comprehensive results at a cost that is advantageous compared with traditional diagnostic methods. Sequencing pathogens directly from a primary clinical specimen would help circumvent the need for culture and, in the process, substantially shorten the time to diagnosis and public health reporting. Unfortunately, this approach poses significant challenges because of the mixture of multiple sequences from a complex fecal biomass. The aim of this project was to develop a proof of concept protocol for the sequencing and genotyping of Shiga toxin-producing Escherichia coli (STEC) directly from stool specimens. We have developed an enrichment protocol that reliably achieves a substantially higher DNA yield belonging to E. coli, which provides adequate next-generation sequencing (NGS) data for downstream bioinformatics analysis. A custom bioinformatics pipeline was created to optimize and remove non-E. coli reads, assess the STEC versus commensal E. coli population in the samples, and build consensus sequences based on population allele frequency distributions. Side-by-side analysis of WGS from paired STEC isolates and matched primary stool specimens reveal that this method can reliably be implemented for many clinical specimens to directly genotype STEC and accurately identify clusters of disease outbreak when no STEC isolate is available for testing.

The impact of HGT on phylogenomic reconstruction methods.

Supermatrix and supertree analyses are frequently used to more accurately recover vertical evolutionary history but debate still exists over which method provides greater reliability. Traditional methods that resolve relationships among organisms from single genes are often unreliable because of the frequent lack of strong phylogenetic signal and the presence of systematic artifacts. Methods developed to reconstruct organismal history from multiple genes can be divided into supermatrix and supertree approaches. A supermatrix analysis consists of the concatenation of multiple genes into a single, possibly partitioned alignment, from which phylogenies are reconstructed using a variety of approaches. Supertrees build consensus trees from the topological information contained within individual gene trees. Both methods are now widely used and have been demonstrated to solve previously ambiguous or unresolved phylogenies with high statistical support. However, the amount of misleading signal needed to induce erroneous phylogenies for both strategies is still unknown. Using genome simulations, we test the accuracy of supertree and supermatrix approaches in recovering the true organismal phylogeny under increased amounts of horizontally transferred genes and changes in substitution rates. Our results show that overall, supermatrix approaches are preferable when a low amount of gene transfer is suspected to be present in the dataset, while supertrees have greater reliability in the presence of a moderate amount of misleading gene transfers. In the face of very high or very low substitution rates without horizontal gene transfers, supermatrix approaches outperform supertrees as individual gene trees remain unresolved and additional sequences contribute to a congruent phylogenetic signal.

Insights into transcriptional changes that accompany organelle sequestration from the stolen nucleus of Mesodinium rubrum.

BACKGROUND: Organelle retention is a form of mixotrophy that allows organisms to reap metabolic benefits similar to those of photoautotrophs through capture of algal prey and sequestration of their plastids. Mesodinium rubrum is an abundant and broadly distributed photosynthetic marine ciliate that steals organelles from cryptophyte algae, such as Geminigera cryophila. M. rubrum is unique from most other acquired phototrophs because it also steals a functional nucleus that facilitates genetic control of sequestered plastids and other organelles. We analyzed changes in G. cryophila nuclear gene expression and transcript abundance after its incorporation into the cellular architecture of M. rubrum as an initial step towards understanding this complex system. METHODS: We compared Illumina-generated transcriptomes of the cryptophyte Geminigera cryophila as a free-living cell and as a sequestered nucleus in M. rubrum to identify changes in protein abundance and gene expression. After KEGG annotation, proteins were clustered by functional categories, which were evaluated for over- or under-representation in the sequestered nucleus. Similarly, coding sequences were grouped by KEGG categories/pathways, which were then evaluated for over- or under-expression via read count strategies. RESULTS: At the time of sampling, the global transcriptome of M. rubrum was dominated (~58-62 %) by transcription from its stolen nucleus. A comparison of transcriptomes from free-living G. cryophila cells to those of the sequestered nucleus revealed a decrease in gene expression and transcript abundance for most functional protein categories within the ciliate. However, genes coding for proteins involved in photosynthesis, oxidative stress reduction, and several other metabolic pathways revealed striking exceptions to this general decline. CONCLUSIONS: Major changes in G. cryophila transcript expression after sequestration by M. rubrum and the ciliate's success as a photoautotroph imply some level of control or gene regulation by the ciliate and at the very least reflect a degree of coordination between host and foreign organelles. Intriguingly, cryptophyte genes involved in protein transport are significantly under-expressed in M. rubrum, implicating a role for the ciliate's endomembrane system in targeting cryptophyte proteins to plastid complexes. Collectively, this initial portrait of an acquired transcriptome within a dynamic and ecologically successful ciliate highlights the remarkable cellular and metabolic chimerism of this system.

Distribution and evolution of the mobile vma-1b intein.

Inteins are self-splicing parasitic genetic elements found in all domains of life. These genetic elements are found in highly conserved positions in conserved proteins. One protein family that has been invaded by inteins is the vacuolar and archaeal catalytic ATPase subunits (vma-1). There are two intein insertion sites in this protein, "a" and "b." The b site was previously thought to be only invaded in archaeal lineages. Here we survey the distribution and evolutionary histories of the b site inteins and show that the intein is present in more lineages than previously annotated, including a bacterial lineage, Mahella australiensis 50-1 BON. We present evidence, through ancestral character state reconstruction and substitution ratios between host genes and inteins, for several transfers of this intein between divergent species, including an interdomain transfer between the archaea and bacteria. Although inteins may persist within a single population or species for long periods of time, transfer of the vma-1b intein between divergent species contributed to the distribution of this intein.

The effects of model choice and mitigating bias on the ribosomal tree of life.

Deep-level relationships within Bacteria, Archaea, and Eukarya as well as the relationships of these three domains to each other require resolution. The ribosomal machinery, universal to all cellular life, represents a protein repertoire resistant to horizontal gene transfer, which provides a largely congruent signal necessary for reconstructing a tree suitable as a backbone for life's reticulate history. Here, we generate a ribosomal tree of life from a robust taxonomic sampling of Bacteria, Archaea, and Eukarya to elucidate deep-level intra-domain and inter-domain relationships. Lack of phylogenetic information and systematic errors caused by inadequate models (that cannot account for substitution rate or compositional heterogeneities) or improper model selection compound conflicting phylogenetic signals from HGT and/or paralogy. Thus, we tested several models of varying sophistication on three different datasets, performed removal of fast-evolving or long-branched Archaea and Eukarya, and employed three different strategies to remove compositional heterogeneity to examine their effects on the topological outcome. Our results support a two-domain topology for the tree of life, where Eukarya emerges from within Archaea as sister to a Korarchaeota/Thaumarchaeota (KT) or Crenarchaeota/KT clade for all models under all or at least one of the strategies employed. Taxonomic manipulation allows single-matrix and certain mixture models to vacillate between two-domain and three-domain phylogenies. We find that models vary in their ability to resolve different areas of the tree of life, which does not necessarily correlate with model complexity. For example, both single-matrix and some mixture models recover monophyletic Crenarchaeota and Euryarchaeota archaeal phyla. In contrast, the most sophisticated model recovers a paraphyletic Euryarchaeota but detects two large clades that comprise the Bacteria, which were recovered separately but never together in the other models. Overall, models recovered consistent topologies despite dataset modifications due to the removal of compositional bias, which reflects either ineffective bias reduction or robust datasets that allow models to overcome reconstruction artifacts. We recommend a comparative approach for evolutionary models to identify model weaknesses as well as consensus relationships.

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.

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.