A laboratory-based predictive pathway for the development of Neisseria gonorrhoeae high-level resistance to corallopyronin A, an inhibitor of bacterial RNA polymerase

The continued emergence of Neisseria gonorrhoeae strains that express resistance to multiple antibiotics, including the last drug for empiric monotherapy (ceftriaxone), necessitates the development of new treatment options to cure gonorrheal infections. Toward this goal, we recently reported that corallopyronin A (CorA), which targets the switch region of the β' subunit (RpoC) of bacterial DNA-dependent RNA polymerase (RNAP), has potent anti-gonococcal activity against a panel of multidrug-resistant clinical strains. Moreover, in that study, CorA could eliminate gonococcal infection of primary human epithelial cells and gonococci in a biofilm state. To determine if N. gonorrhoeae could develop high-level resistance to CorA in a single step, we sought to isolate spontaneous mutants expressing any CorA resistance phenotypes. However, no single-step mutants with high-level CorA resistance were isolated. High-level CorA resistance could only be achieved in this study through a multi-step pathway involving over-expression of the MtrCDE drug efflux pump and single amino acid changes in the β and β' subunits (RpoB and RpoC, respectively) of RNAP. Molecular modeling of RpoB and RpoC interacting with CorA was used to deduce how the amino acid changes in RpoB and RpoC could influence gonococcal resistance to CorA. Bioinformatic analyses of whole genome sequences of clinical gonococcal isolates indicated that the CorA resistance determining mutations in RpoB/C, identified herein, are very rare (≤ 0.0029%), suggesting that the proposed pathway for resistance is predictive of how this phenotype could potentially evolve if CorA is used therapeutically to treat gonorrhea in the future.

IMPORTANCE

The continued emergence of multi-antibiotic-resistant strains of Neisseria gonorrhoeae necessitates the development of new antibiotics that are effective against this human pathogen. We previously described that the RNA polymerase-targeting antibiotic corallopyronin A (CorA) has potent activity against a large collection of clinical strains that express different antibiotic resistance phenotypes including when such gonococci are in a biofilm state. Herein, we tested whether a CorA-sensitive gonococcal strain could develop spontaneous resistance. Our finding that CorA resistance could only be achieved by a multi-step process involving over-expression of the MtrCDE efflux pump and single amino acid changes in RpoB and RpoC suggests that such resistance may be difficult for gonococci to evolve if this antibiotic is used in the future to treat gonorrheal infections that are refractory to cure by other antibiotics.

Incomplete immunity in a natural animal-microbiota interaction selects for higher pathogen virulence

Incomplete immunity in recovered hosts is predicted to favor more virulent pathogens upon re-infection in the population.1 The microbiota colonizing animals can generate a similarly long-lasting, partial immune response, allowing for infection but dampened disease severity.2 We tracked the evolutionary trajectories of a widespread pathogen (Pseudomonas aeruginosa), experimentally passaged through populations of nematodes immune-primed by a natural microbiota member (P. berkeleyensis). This bacterium can induce genes regulated by a mitogen-activated protein kinase (MAPK) signaling pathway effective at conferring protection against pathogen-induced death despite infection.3 Across host populations, this incomplete immunity selected for pathogens more than twice as likely to kill as those evolved in non-primed (i.e., naive) or immune-compromised (mutants with a knockout of the MAPK ortholog) control populations. Despite the higher virulence, pathogen molecular evolution in immune-primed hosts was slow and constrained. In comparison, evolving pathogens in immune-compromised hosts were characterized by substantial genomic differentiation and attenuated virulence. These findings directly attribute the incomplete host immunity induced from microbiota as a significant force shaping the virulence and evolutionary dynamics of novel infectious diseases.

Experimentally evolved Staphylococcus aureus shows increased survival in the presence of Pseudomonas aeruginosa by acquiring mutations in the amino acid transporter, GltT

When cultured together under standard laboratory conditions Pseudomonas aeruginosa has been shown to be an effective inhibitor of Staphylococcus aureus. However, P. aeruginosa and S. aureus are commonly observed in coinfections of individuals with cystic fibrosis (CF) and in chronic wounds. Previous work from our group revealed that S. aureus isolates from CF infections are able to persist in the presence of P. aeruginosa strain PAO1 with a range of tolerances with some isolates being eliminated entirely and others maintaining large populations. In this study, we designed a serial transfer, evolution experiment to identify mutations that allow S. aureus to survive in the presence of P. aeruginosa. Using S. aureus USA300 JE2 as our ancestral strain, populations of S. aureus were repeatedly cocultured with fresh P. aeruginosa PAO1. After eight coculture periods, S. aureus populations that survived better in the presence of PAO1 were observed. We found two independent mutations in the highly conserved S. aureus aspartate transporter, gltT, that were unique to evolved P. aeruginosa-tolerant isolates. Subsequent phenotypic testing demonstrated that gltT mutants have reduced uptake of glutamate and outcompeted wild-type S. aureus when glutamate was absent from chemically defined media. These findings together demonstrate that the presence of P. aeruginosa exerts selective pressure on S. aureus to alter its uptake and metabolism of key amino acids when the two are cultured together.

The Burden and Impact of Early Post-transplant Multidrug-Resistant Organism Detection Among Renal Transplant Recipients, 2005-2021

Background

Reducing the burden of multidrug-resistant organism (MDRO) colonization and infection among renal transplant recipients (RTRs) may improve patient outcomes. We aimed to assess whether the detection of an MDRO or a comparable antibiotic-susceptible organism (CSO) during the early post-transplant (EPT) period was associated with graft loss and mortality among RTRs.

Methods

We conducted a retrospective cohort study of RTRs transplanted between 2005 and 2021. EPT positivity was defined as a positive bacterial culture within 30 days of transplant. The incidence and prevalence of EPT MDRO detection were calculated. The primary outcome was a composite of 1-year allograft loss or mortality following transplant. Multivariable Cox hazard regression, competing risk, propensity score-weighted sensitivity, and subgroup analyses were performed.

Results

Among 3507 RTRs, the prevalence of EPT MDRO detection was 1.3% (95% CI, 0.91%-1.69%) with an incidence rate per 1000 EPT-days at risk of 0.42 (95% CI, 0.31-0.57). Among RTRs who met survival analysis inclusion criteria (n = 3432), 91% (3138/3432) had no positive EPT cultures and were designated as negative controls, 8% (263/3432) had a CSO detected, and 1% (31/3432) had an MDRO detected in the EPT period. EPT MDRO detection was associated with the composite outcome (adjusted hazard ratio [aHR], 3.29; 95% CI, 1.21-8.92) and death-censored allograft loss (cause-specific aHR, 7.15; 95% CI, 0.92-55.5; subdistribution aHR, 7.15; 95% CI, 0.95-53.7). A similar trend was seen in the subgroup and sensitivity analyses.

Conclusions

MDRO detection during the EPT period was associated with allograft loss, suggesting the need for increased strategies to optimize prevention of MDRO colonization and infection.

Factors associated with foreign body infection in methicillin-resistant Staphylococcus aureus bacteremia

Background

We aimed to compare patient characteristics, MRSA sequence types, and biofilm production of MRSA strains that did and did not cause a foreign body infection in patients with MRSA bloodstream infections (BSI).

Methods

All adult patients with MRSA BSI hospitalized in two hospitals were identified by clinical microbiology laboratory surveillance. Only patients who had at least one implanted foreign body during the episode of BSI were included.

Results

In July 2018 - March 2022, of 423 patients identified with MRSA BSI, 118 (28%) had ≥1 foreign body. Among them, 51 (43%) had one or more foreign body infections. In multivariable analysis, factors associated with foreign body infection were history of MRSA infection in the last year (OR=4.7 [1.4-15.5], p=0.012) community-associated BSI (OR=68.1 [4.2-1114.3], p=0.003); surgical site infection as source of infection (OR=11.8 [2-70.4], p=0.007); presence of more than one foreign body (OR=3.4 [1.1-10.7], p=0.033); interval between foreign body implantation and infection <18 months (OR=3.3 [1.1-10], p=0.031); and positive blood culture ≥48h (OR=16.7 [4.3-65.7], p<0.001). The most prevalent sequence type was ST8 (39%), followed by ST5 (29%), and ST105 (20%) with no significant difference between patients with or without foreign body infection. Only 39% of MRSA isolates formed a moderate/strong biofilm. No significant difference was observed between patients with foreign body infection and those without foreign body infection. In multivariable analysis, subjects infected with a MRSA isolate producing moderate/strong in vitro biofilm were more likely to have a history of MRSA infection in the last year (OR=3.41 [1.23-9.43]), interval between foreign body implantation and MRSA BSI <18 months (OR=3.1 [1.05-9.2]) and ST8 (OR=10.64 [2-57.3]).

Conclusion

Most factors associated with foreign body infection in MRSA BSI were also characteristic of persistent infections. Biofilm-forming isolates were not associated with a higher risk of foreign-body infection but appeared to be associated with MRSA genetic lineage, especially ST8.

The role of hypermutation and collateral sensitivity in antimicrobial resistance diversity of Pseudomonas aeruginosa populations in cystic fibrosis lung infection

Pseudomonas aeruginosa is an opportunistic pathogen which causes chronic, drug-resistant lung infections in cystic fibrosis (CF) patients. In this study, we explore the role of genomic diversification and evolutionary trade-offs in antimicrobial resistance (AMR) diversity within P. aeruginosa populations sourced from CF lung infections. We analyzed 300 clinical isolates from four CF patients (75 per patient) and found that genomic diversity is not a consistent indicator of phenotypic AMR diversity. Remarkably, some genetically less diverse populations showed AMR diversity comparable to those with significantly more genetic variation. We also observed that hypermutator strains frequently exhibited increased sensitivity to antimicrobials, contradicting expectations from their treatment histories. Investigating potential evolutionary trade-offs, we found no substantial evidence of collateral sensitivity among aminoglycoside, beta-lactam, or fluoroquinolone antibiotics, nor did we observe trade-offs between AMR and growth in conditions mimicking CF sputum. Our findings suggest that (i) genomic diversity is not a prerequisite for phenotypic AMR diversity, (ii) hypermutator populations may develop increased antimicrobial sensitivity under selection pressure, (iii) collateral sensitivity is not a prominent feature in CF strains, and (iv) resistance to a single antibiotic does not necessarily lead to significant fitness costs. These insights challenge prevailing assumptions about AMR evolution in chronic infections, emphasizing the complexity of bacterial adaptation during infection.IMPORTANCEUpon infection in the cystic fibrosis (CF) lung, Pseudomonas aeruginosa rapidly acquires genetic mutations, especially in genes involved in antimicrobial resistance (AMR), often resulting in diverse, treatment-resistant populations. However, the role of bacterial population diversity within the context of chronic infection is still poorly understood. In this study, we found that hypermutator strains of P. aeruginosa in the CF lung undergoing treatment with tobramycin evolved increased sensitivity to tobramycin relative to non-hypermutators within the same population. This finding suggests that antimicrobial treatment may only exert weak selection pressure on P. aeruginosa populations in the CF lung. We further found no evidence for collateral sensitivity in these clinical populations, suggesting that collateral sensitivity may not be a robust, naturally occurring phenomenon for this microbe.

Average Nucleotide Identity based Staphylococcus aureus strain grouping allows identification of strain-specific genes in the pangenome

Staphylococcus aureus causes both hospital and community acquired infections in humans worldwide. Due to the high incidence of infection S. aureus is also one of the most sampled and sequenced pathogens today, providing an outstanding resource to understand variation at the bacterial subspecies level. We processed and downsampled 83,383 public S. aureus Illumina whole genome shotgun sequences and 1,263 complete genomes to produce 7,954 representative substrains. Pairwise comparison of core gene Average Nucleotide Identity (ANI) revealed a natural boundary of 99.5% that could be used to define 145 distinct strains within the species. We found that intermediate frequency genes in the pangenome (present in 10-95% of genomes) could be divided into those closely linked to strain background ("strain-concentrated") and those highly variable within strains ("strain-diffuse"). Non-core genes had different patterns of chromosome location; notably, strain-diffuse associated with prophages, strain-concentrated with the vSaβ genome island and rare genes (<10% frequency) concentrated near the origin of replication. Antibiotic genes were enriched in the strain-diffuse class, while virulence genes were distributed between strain-diffuse, strain-concentrated, core and rare classes. This study shows how different patterns of gene movement help create strains as distinct subspecies entities and provide insight into the diverse histories of important S. aureus functions.

Unique microbial diversity, community composition, and networks among Pacific Islander endocervical and vaginal microbiomes with and without Chlamydia trachomatis infection in Fiji

IMPORTANCE

Chlamydia trachomatis (Ct) is the most common sexually transmitted bacterium globally. Endocervical and vaginal microbiome interactions are rarely examined within the context of Ct or among vulnerable populations. We evaluated 258 vaginal and 92 paired endocervical samples from Fijian women using metagenomic shotgun sequencing. Over 37% of the microbiomes could not be classified into sub-community state types (subCSTs). We, therefore, developed subCSTs IV-D0, IV-D1, IV-D2, and IV-E-dominated primarily by Gardnerella vaginalis-to improve classification. Among paired microbiomes, the endocervix had a significantly higher alpha diversity and, independently, higher diversity for high-risk human papilloma virus (HPV) genotypes compared to low-risk and no HPV. Ct-infected endocervical networks had smaller clusters without interactions with potentially beneficial Lactobacillus spp. Overall, these data suggest that G. vaginalis may generate polymicrobial biofilms that predispose to and/or promote Ct and possibly HPV persistence and pathogenicity. Our findings expand on the existing repertoire of endocervical and vaginal microbiomes and fill in knowledge gaps regarding Pacific Islanders.

The Role of Hypermutation and Collateral Sensitivity in Antimicrobial Resistance Diversity of Pseudomonas aeruginosa Populations in Cystic Fibrosis Lung Infection

Pseudomonas aeruginosa is an opportunistic pathogen which causes chronic, drug-resistant lung infections in cystic fibrosis (CF) patients. In this study, we explore the role of genomic diversification and evolutionary trade-offs in antimicrobial resistance (AMR) diversity within P. aeruginosa populations sourced from CF lung infections. We analyzed 300 clinical isolates from four CF patients (75 per patient), and found that genomic diversity is not a consistent indicator of phenotypic AMR diversity. Remarkably, some genetically less diverse populations showed AMR diversity comparable to those with significantly more genetic variation. We also observed that hypermutator strains frequently exhibited increased sensitivity to antimicrobials, contradicting expectations from their treatment histories. Investigating potential evolutionary trade-offs, we found no substantial evidence of collateral sensitivity among aminoglycoside, beta-lactam, or fluoroquinolone antibiotics, nor did we observe trade-offs between AMR and growth in conditions mimicking CF sputum. Our findings suggest that (i) genomic diversity is not a prerequisite for phenotypic AMR diversity; (ii) hypermutator populations may develop increased antimicrobial sensitivity under selection pressure; (iii) collateral sensitivity is not a prominent feature in CF strains, and (iv) resistance to a single antibiotic does not necessarily lead to significant fitness costs. These insights challenge prevailing assumptions about AMR evolution in chronic infections, emphasizing the complexity of bacterial adaptation during infection.

Comparison of genomic diversity between single and pooled Staphylococcus aureus colonies isolated from human colonization cultures

The most common approach to sampling the bacterial populations within an infected or colonized host is to sequence genomes from a single colony obtained from a culture plate. However, it is recognized that this method does not capture the genetic diversity in the population. Sequencing a mixture of several colonies (pool-seq) is a better approach to detect population heterogeneity, but it is more complex to analyse due to different types of heterogeneity, such as within-clone polymorphisms, multi-strain mixtures, multi-species mixtures and contamination. Here, we compared 8 single-colony isolates (singles) and pool-seq on a set of 2286 Staphylococcus aureus culture samples to identify features that can distinguish pure samples, samples undergoing intraclonal variation and mixed strain samples. The samples were obtained by swabbing 3 body sites on 85 human participants quarterly for a year, who initially presented with a methicillin-resistant S. aureus skin and soft-tissue infection (SSTI). We compared parameters such as sequence quality, contamination, allele frequency, nucleotide diversity and pangenome diversity in each pool to those for the corresponding singles. Comparing singles from the same culture plate, we found that 18% of sample collections contained mixtures of multiple multilocus sequence types (MLSTs or STs). We showed that pool-seq data alone could predict the presence of multi-ST populations with 95% accuracy. We also showed that pool-seq could be used to estimate the number of intra-clonal polymorphic sites in the population. Additionally, we found that the pool may contain clinically relevant genes such as antimicrobial resistance markers that may be missed when only examining singles. These results highlight the potential advantage of analysing genome sequences of total populations obtained from clinical cultures rather than single colonies.

Integrative analysis of microbial 16S gene and shotgun metagenomic sequencing data improves statistical efficiency

Background

The most widely used technologies for profiling microbial communities are 16S marker-gene sequencing and shotgun metagenomic sequencing. Interestingly, many microbiome studies have performed both sequencing experiments on the same cohort of samples. The two sequencing datasets often reveal consistent patterns of microbial signatures, highlighting the potential for an integrative analysis to improve power of testing these signatures. However, differential experimental biases, partially overlapping samples, and differential library sizes pose tremendous challenges when combining the two datasets. Currently, researchers either discard one dataset entirely or use different datasets for different objectives.

Methods

In this article, we introduce the first method of this kind, named Com-2seq, that combines the two sequencing datasets for testing differential abundance at the genus and community levels while overcoming these difficulties. The new method is based on our LOCOM model (Hu et al., 2022), which employs logistic regression for testing taxon differential abundance while remaining robust to experimental bias. To benchmark the performance of Com-2seq, we introduce two ad hoc approaches: applying LOCOM to pooled taxa count data and combining LOCOM p-values from analyzing each dataset separately.

Results

Our simulation studies indicate that Com-2seq substantially improves statistical efficiency over analysis of either dataset alone and works better than the two ad hoc approaches. An application of Com-2seq to two real microbiome studies uncovered scientifically plausible findings that would have been missed by analyzing individual datasets.

Conclusions

Com-2seq performs integrative analysis of 16S and metagenomic sequencing data, which improves statistical efficiency and has the potential to accelerate the search of microbial communities and taxa that are involved in human health and diseases.

Transcriptional activation of ompA in Neisseria gonorrhoeae mediated by the XRE family member protein NceR

Increasing antibiotic resistance of Neisseria gonorrhoeae, the causative agent of gonorrhea, is a growing global concern that has renewed vaccine development efforts. The gonococcal OmpA protein was previously identified as a vaccine candidate due to its surface exposure, conservation, stable expression, and involvement in host-cell interactions. We previously demonstrated that the transcription of ompA can be activated by the MisR/MisS two-component system. Interestingly, earlier work suggested that the availability of free iron also influences ompA expression, which we confirmed in this study. In the present study, we found that iron regulation of ompA was independent of MisR and searched for additional regulators. A DNA pull-down assay with the ompA promoter from gonococcal lysates obtained from bacteria grown in the presence or absence of iron identified an XRE (Xenobiotic Response Element) family member protein encoded by NGO1982. We found that an NGO1982 null mutant of N. gonorrhoeae strain FA19 displayed a reduced level of ompA expression compared to the wild-type (WT) parent strain. Given this regulation, and the capacity of this XRE-like protein to regulate a gene involved in peptidoglycan biosynthesis (ltgA), along with its presence in other Neisseria sp., we termed the NGO1982-encoded protein as NceR (Neisseria cell envelope regulator). Critically, results from DNA-binding studies indicated that NceR regulates ompA through a direct mechanism. Thus, ompA expression is subject to both iron-dependent (NceR) and -independent (MisR/MisS) pathways. Hence, levels of the vaccine antigen candidate OmpA in circulating gonococcal strains could be influenced by transcriptional regulatory systems and the availability of iron. IMPORTANCE Herein, we report that the gene encoding a conserved gonococcal surface-exposed vaccine candidate (OmpA) is activated by a heretofore undescribed XRE family transcription factor, which we term NceR. We report that NceR regulation of ompA expression in N. gonorrhoeae is mediated by an iron-dependent mechanism, while the previously described MisR regulatory system is iron-independent. Our study highlights the importance of defining the complexity of coordinated genetic and physiologic systems that regulate genes encoding vaccine candidates to better understand their availability during infection.

Experimentally Evolved Staphylococcus aureus Survives in the Presence of Pseudomonas aeruginosa by Acquiring Mutations in the Amino Acid Transporter, GltT

Staphylococcus aureus and Pseudomonas aeruginosa are the most common bacterial pathogens isolated from cystic fibrosis (CF) related lung infections. When both of these opportunistic pathogens are found in a coinfection, CF patients tend to have higher rates of pulmonary exacerbations and experience a more rapid decrease in lung function. When cultured together under standard laboratory conditions, it is often observed that P. aeruginosa effectively inhibits S. aureus growth. Previous work from our group revealed that S. aureus from CF infections have isolate-specific survival capabilities when cocultured with P. aeruginosa. In this study, we designed a serial transfer evolution experiment to identify mutations that allow S. aureus to adapt to the presence of P. aeruginosa. Using S. aureus USA300 JE2 as our ancestral strain, populations of S. aureus were repeatedly cocultured with fresh P. aeruginosa strain, PAO1. After 8 coculture periods, S. aureus populations that survived better in the presence of PAO1 were observed. We found two independent mutations in the highly conserved S. aureus aspartate transporter, gltT, that were unique to evolved P. aeruginosa-tolerant isolates. Subsequent phenotypic testing demonstrated that gltT mutants have reduced uptake of glutamate and outcompete wild-type S. aureus when glutamate is absent from chemically-defined media. These findings together demonstrate that the presence of P. aeruginosa exerts selective pressure on S. aureus to alter its uptake and metabolism of key amino acids when the two bacteria are cultured together.

Integrative analysis of microbial 16S gene and shotgun metagenomic sequencing data improves statistical efficiency

The most widely used technologies for profiling microbial communities are 16S marker-gene sequencing and shotgun metagenomic sequencing. Interestingly, many microbiome studies have performed both sequencing experiments on the same cohort of samples. The two sequencing datasets often reveal consistent patterns of microbial signatures, highlighting the potential for an integrative analysis to improve power of testing these signatures. However, differential experimental biases, partially overlapping samples, and differential library sizes pose tremendous challenges when combining the two datasets. Currently, researchers either discard one dataset entirely or use different datasets for different objectives. In this article, we introduce the first method of this kind, named Com-2seq, that combines the two sequencing datasets for the objective of testing differential abundance at the genus and community levels while overcoming these difficulties. We demonstrate that Com-2seq substantially improves statistical efficiency over analysis of either dataset alone and works better than two ad hoc approaches.

Comparison of genomic diversity between single and pooled Staphylococcus aureus colonies isolated from human colonisation cultures

The most common approach to sampling the bacterial populations within an infected or colonised host is to sequence genomes from a single colony obtained from a culture plate. However, it is recognized that this method does not capture the genetic diversity in the population. An alternative is to sequence a mixture containing multiple colonies ("pool-seq"), but this has the disadvantage that it is a non-homogeneous sample, making it difficult to perform specific experiments. We compared differences in measures of genetic diversity between eight single-colony isolates (singles) and pool-seq on a set of 2286 S. aureus culture samples. The samples were obtained by swabbing three body sites on 85 human participants quarterly for a year, who initially presented with a methicillin-resistant S. aureus skin and soft-tissue infection (SSTI). We compared parameters such as sequence quality, contamination, allele frequency, nucleotide diversity and pangenome diversity in each pool to the corresponding singles. Comparing singles from the same culture plate, we found that 18% of sample collections contained mixtures of multiple Multilocus sequence types (MLSTs or STs). We showed that pool-seq data alone could predict the presence of multi-ST populations with 95% accuracy. We also showed that pool-seq could be used to estimate the number of polymorphic sites in the population. Additionally, we found that the pool may contain clinically relevant genes such as antimicrobial resistance markers that may be missed when only examining singles. These results highlight the potential advantage of analysing genome sequences of total populations obtained from clinical cultures rather than single colonies.

Detecting patterns of accessory genome coevolution in Staphylococcus aureus using data from thousands of genomes

Bacterial genomes exhibit widespread horizontal gene transfer, resulting in highly variable genome content that complicates the inference of genetic interactions. In this study, we develop a method for detecting coevolving genes from large datasets of bacterial genomes based on pairwise comparisons of closely related individuals, analogous to a pedigree study in eukaryotic populations. We apply our method to pairs of genes from the Staphylococcus aureus accessory genome of over 75,000 annotated gene families using a database of over 40,000 whole genomes. We find many pairs of genes that appear to be gained or lost in a coordinated manner, as well as pairs where the gain of one gene is associated with the loss of the other. These pairs form networks of rapidly coevolving genes, primarily consisting of genes involved in virulence, mechanisms of horizontal gene transfer, and antibiotic resistance, particularly the SCCmec complex. While we focus on gene gain and loss, our method can also detect genes that tend to acquire substitutions in tandem, or genotype-phenotype or phenotype-phenotype coevolution. Finally, we present the R package DeCoTUR that allows for the computation of our method.

Patterns of within-host spread of Chlamydia trachomatis between vagina, endocervix and rectum revealed by comparative genomic analysis

Introduction

Chlamydia trachomatis, a gram-negative obligate intracellular bacterium, commonly causes sexually transmitted infections (STIs). Little is known about C. trachomatis transmission within the host, which is important for understanding disease epidemiology and progression.

Methods

We used RNA-bait enrichment and whole-genome sequencing to compare rectal, vaginal and endocervical samples collected at the same time from 26 study participants who attended Fijian Ministry of Health and Medical Services clinics and tested positive for C. trachomatis at each anatomic site.

Results

The 78 C. trachomatis genomes from participants resolved into two major clades of the C. trachomatis phylogeny (the “prevalent urogenital and anorectal” clade and “non-prevalent urogenital and anorectal” clade). For 21 participants, genome sequences were almost identical in each anatomic site. For the other five participants, two distinct C. trachomatis strains were present in different sites; in two cases, the vaginal sample was a mixture of strains.

Discussion

The absence of large numbers of fixed SNPs between C. trachomatis genomes within many of the participants could indicate recent acquisition of infection prior to the clinic visit without sufficient time to accumulate significant genetic variation in different body sites. This model suggests that many C. trachomatis infections may be resolved relatively quickly in the Fijian population, possibly reflecting common prescription or over-the-counter antibiotics usage.

Development of an Amplicon Nanopore Sequencing Strategy for Detection of Mutations Conferring Intermediate Resistance to Vancomycin in Staphylococcus aureus Strains

Staphylococcus aureus is a major cause of bacteremia and other hospital-acquired infections. The cell-wall active antibiotic vancomycin is commonly used to treat both methicillin-resistant (MRSA) and sensitive (MSSA) infections. Vancomycin intermediate S. aureus (VISA) variants can arise through de novo mutations. Here, we performed pilot experiments to develop a combined PCR/long-read sequencing-based method for detection of previously known VISA-causing mutations. Primers were designed to generate 10 amplicons covering 16 genes associated with the VISA phenotype. We sequenced amplicon pools as long reads with Oxford Nanopore adapter ligation on Flongle flow cells. We then detected mutations by mapping reads against a parental consensus or known reference sequence and comparing called variants against a database of known VISA mutations from laboratory selection. Each amplicon in the pool was sequenced to high (>1,000×) coverage, and no relationship was found between amplicon length and coverage. We also were able to detect the causative mutation (walK 646C>G) in a VISA mutant derived from the USA300 strain (N384-3 from parental strain N384). Mixing mutant (N384-3) and parental (N384) DNA at various ratios from 0 to 1 mutant suggested a mutation detection threshold of the average minor allele frequency (6.5%) at 95% confidence (two standard errors above mean mutation frequency). The study lays the groundwork for direct S. aureus antibiotic resistance genotype inference using rapid nanopore sequencing from clinical samples. IMPORTANCE Bacteremia mortality is known to increase rapidly with time after infection, making rapid diagnostics and treatment necessary. Successful treatment depends on correct administration of antibiotics based on knowledge of strain antibiotic susceptibility. Staphylococcus aureus is a major causative agent of bacteremia that is also commonly antibiotic resistant. In this work, we develop a method to accelerate detection of a complex, polygenic antibiotic resistance phenotype in S. aureus, vancomycin-intermediate resistance (VISA), through long-read genomic sequencing of amplicons representing genes most commonly mutated in VISA selection. This method both rapidly identifies VISA genotypes and incorporates the most comprehensive database of VISA genetic determinants known to date.

Patterns of within-host spread of Chlamydia trachomatis between vagina, endocervix and rectum revealed by comparative genomic analysis

Chlamydia trachomatis , a gram-negative obligate intracellular bacterium, commonly causes sexually transmitted infections (STIs). Little is known about C. trachomatis transmission within the host, which is important for understanding disease epidemiology and progression. We used RNA-bait enrichment and whole-genome sequencing to compare rectal, vaginal and endocervical samples collected at the same time from 26 study participants who attended Fijian Ministry of Health and Medical Services clinics and tested positive for C. trachomatis at each anatomic site. The 78 C. trachomatis genomes from participants were from two major clades of the C. trachomatis phylogeny (the "prevalent urogenital and anorecta"l clade and "non-prevalent urogenital and anorectal" clade). For 21 participants, genome sequences were almost identical in each anatomic site. For the other five participants, two distinct C. trachomatis strains were present in different sites; in two cases, the vaginal sample was a mixture of strains. The absence of large numbers of fixed SNPs between C. trachomatis strains within many of the participants could indicate recent acquisition of infection prior to the clinic visit without sufficient time to accumulate significant variation in the different body sites. This model suggests that many C. trachomatis infections may be resolved relatively quickly in the Fijian population, possibly reflecting common prescription or over-the-counter antibiotics usage.

Importance

Chlamydia trachomatis is a bacterial pathogen that causes millions of sexually transmitted infections (STIs) annually across the globe. Because C. trachomatis lives inside human cells, it has historically been hard to study. We know little about how the bacterium spreads between body sites. Here, samples from 26 study participants who had simultaneous infections in their vagina, rectum and endocervix were genetically analyzed using an improved method to extract C. trachomatis DNA directly from clinical samples for genome sequencing. By analyzing patterns of mutations in the genomes, we found that 21 participants shared very similar C. trachomatis strains in all three anatomic sites, suggesting recent infection and spread. For five participants two C. trachomatis strains were evident, indicating multiple infections. This study is significant in that improved enrichment methods for genome sequencing provides robust data to genetically trace patterns of C. trachomatis infection and transmission within an individual for epidemiologic and pathogenesis interrogations.

Pre-epidemic evolution of the MRSA USA300 clade and a molecular key for classification

Introduction

USA300 has remained the dominant community and healthcare associated methicillin-resistant Staphylococcus aureus (MRSA) clone in the United States and in northern South America for at least the past 20 years. In this time, it has experienced epidemic spread in both of these locations. However, its pre-epidemic evolutionary history and origins are incompletely understood. Large sequencing databases, such as NCBI, PATRIC, and Staphopia, contain clues to the early evolution of USA300 in the form of sequenced genomes of USA300 isolates that are representative of lineages that diverged prior to the establishment of the South American epidemic (SAE) clade and North American epidemic (NAE) clade. In addition, historical isolates collected prior to the emergence of epidemics can help reconstruct early events in the history of this lineage.

Methods

Here, we take advantage of the accrued, publicly available data, as well as two newly sequenced pre-epidemic historical isolates from 1996, and a very early diverging ACME-negative NAE genome, to understand the pre-epidemic evolution of USA300. We use database mining techniques to emphasize genomes similar to pre-epidemic isolates, with the goal of reconstructing the early molecular evolution of the USA300 lineage.

Results

Phylogenetic analysis with these genomes confirms that the NAE and SAE USA300 lineages diverged from a most recent common ancestor around 1970 with high confidence, and it also pinpoints the independent acquisition events of the of the ACME and COMER loci with greater precision than in previous studies. We provide evidence for a North American origin of the USA300 lineage and identify multiple introductions of USA300 into South and North America. Notably, we describe a third major USA300 clade (the pre-epidemic branching clade; PEB1) consisting of both MSSA and MRSA isolates circulating around the world that diverged from the USA300 lineage prior to the establishment of the South and North American epidemics. We present a detailed analysis of specific sequence characteristics of each of the major clades, and present diagnostic positions that can be used to classify new genomes.

Unsuspected Clonal Spread of Methicillin-Resistant Staphylococcus aureus Causing Bloodstream Infections in Hospitalized Adults Detected Using Whole Genome Sequencing

BACKGROUND

Though detection of transmission clusters of methicillin-resistant Staphylococcus aureus (MRSA) infections is a priority for infection control personnel in hospitals, the transmission dynamics of MRSA among hospitalized patients with bloodstream infections (BSIs) has not been thoroughly studied. Whole genome sequencing (WGS) of MRSA isolates for surveillance is valuable for detecting outbreaks in hospitals, but the bioinformatic approaches used are diverse and difficult to compare.

METHODS

We combined short-read WGS with genotypic, phenotypic, and epidemiological characteristics of 106 MRSA BSI isolates collected for routine microbiological diagnosis from inpatients in 2 hospitals over 12 months. Clinical data and hospitalization history were abstracted from electronic medical records. We compared 3 genome sequence alignment strategies to assess similarity in cluster ascertainment. We conducted logistic regression to measure the probability of predicting prior hospital overlap between clustered patient isolates by the genetic distance of their isolates.

RESULTS

While the 3 alignment approaches detected similar results, they showed some variation. A gene family-based alignment pipeline was most consistent across MRSA clonal complexes. We identified 9 unique clusters of closely related BSI isolates. Most BSIs were healthcare associated and community onset. Our logistic model showed that with 13 single-nucleotide polymorphisms, the likelihood that any 2 patients in a cluster had overlapped in a hospital was 50%.

CONCLUSIONS

Multiple clusters of closely related MRSA isolates can be identified using WGS among strains cultured from BSI in 2 hospitals. Genomic clustering of these infections suggests that transmission resulted from a mix of community spread and healthcare exposures long before BSI diagnosis.

1377. Bloodstream and Skin Infection MRSA Isolates, 2019-2021: Strain Differences and Phylogenetic Clustering in a Single Health System

Background

Methicillin-resistant S. aureus (MRSA) is a common antibiotic-resistant human pathogen that spreads from person to person. Asymptomatic host colonization precedes both skin and soft tissue infections (SSTI) and bloodstream infections (BSI), but it is not known how closely related MRSA strains are that cause infections at different body sites.

Methods

Using Illumina whole genome sequencing, we compared the strain types and genomes of MRSA isolates from 132 SSTIs and 145 sequential BSIs at 3 Philadelphia hospitals in 7/2018-1/2021. We investigated the epidemiological links and genomic clusters among isolates causing BSI and SSTIs.

Results

Abscesses were the most common type of SSTIs (65/132, 49%). Clonal complex (CC) 8 was the most identified CC among both SSTI strains (102/132, 77%) and BSI strains (73/145, 50%). While CC5 was more commonly found among BSIs than SSTIs (50/145, 34% vs. 19/132, 13%). Outbreak clusters with 15 or fewer SNPs were identified. Three clusters that contained only strains from SSTIs were all CC8, and the five clusters that only contained strains from BSIs had three CC5 clusters and two CC8 clusters. Among eight clusters that included at least one SSTI and one BSI strain, only one was composed of CC5 and the others were all CC8 strains. Within these eight clusters the SSTIs were from cellulitis infections (4/8), an abscess (3/8), or an infected ulcer (1/8).

Conclusion

We found that CC8 strains were the most common among both SSTIs and BSIs. CC5 strains, however, were more commonly found in BSIs than SSTIs, and among outbreak clusters in the BSI strains. Eight outbreak clusters were identified that linked strains causing SSTIs or BSIs, seven of the clusters contained only CC8 strains. While abscesses were the most common infection type to cause an SSTI, within the outbreak clusters containing strains from both SSTIs and BSIs, both cellulitis and abscess were equally identified.

Disclosures

Michael Z. David, MD PhD, Contrafect: Grant/Research Support|GSK: Advisor/Consultant|Johnson and Johnson: Advisor/Consultant.

Metagenomic Shotgun Sequencing of Endocervical, Vaginal, and Rectal Samples among Fijian Women with and without Chlamydia trachomatis Reveals Disparate Microbial Populations and Function across Anatomic Sites: a Pilot Study

Chlamydia trachomatis is a sexually transmitted pathogen and a global public health concern. Little is known about the microbial composition and function across endocervical, vaginal, and rectal microbiomes in the context of C. trachomatis infection. We evaluated the microbiomes of 10 age-matched high-risk Fijian women with and without C. trachomatis using metagenomic shotgun sequencing (MSS). Lactobacillus iners and Lactobacillus crispatus dominated the vagina and endocervix of uninfected women. Species often found in higher relative abundance in bacterial vaginosis (BV)-Mageeibacillus indolicus, Prevotella spp., Sneathia spp., Gardnerella vaginalis, and Veillonellaceae spp.-were dominant in C. trachomatis-infected women. This combination of BV pathogens was unique to Pacific Islanders compared to previously studied groups. The C. trachomatis-infected endocervix had a higher diversity of microbiota and microbial profiles that were somewhat different from those of the vagina. However, community state type III (CST-III) and CST-IV predominated, reflecting pathogenic microbiota regardless of C. trachomatis infection status. Rectal microbiomes were dominated by Prevotella and Bacteroides, although four women had unique microbiomes with Gardnerella, Akkermansia, Bifidobacterium, and Brachyspira. A high level of microbial similarity across microbiomes in two C. trachomatis-infected women suggested intragenitorectal transmission. A number of metabolic pathways in the endocervix, driven by BV pathogens and C. trachomatis to meet nutritional requirements for survival/growth, 5-fold higher than that in the vagina indicated that endocervical microbial functions are likely more diverse and complex than those in the vagina. Our novel findings provide the impetus for larger prospective studies to interrogate microbial/microbiome interactions that promote C. trachomatis infection and better define the unique genitorectal microbiomes of Pacific Islanders. IMPORTANCE Chlamydia trachomatis is the primary cause of bacterial sexually transmitted infections worldwide, with a disturbing increase in annual rates. While there is a plethora of data on healthy and pathogenic vaginal microbiomes-defining microbial profiles and associations with sexually transmitted infections (STIs)-far fewer studies have similarly examined the endocervix or rectum. Further, vulnerable populations, such as Pacific Islanders, remain underrepresented in research. We investigated the microbial composition, structure, and function of these anatomic microbiomes using metagenomic shotgun sequencing among a Fijian cohort. We found, primarily among C. trachomatis-infected women, unique microbial profiles in endocervical, vaginal, and rectal microbiomes with an increased diversity and more complex microbial pathways in endocervical than vaginal microbiomes. Similarities in microbiome composition across sites for some women suggested intragenitorectal transmission. These novel insights into genitorectal microbiomes and their purported function require prospective studies to better define Pacific Islander microbiomes and microbial/microbiome interactions that promote C. trachomatis infection.

Association of Secretor Status and Recent Norovirus Infection With Gut Microbiome Diversity Metrics in a Veterans Affairs Population

Norovirus infection causing acute gastroenteritis could lead to adverse effects on the gut microbiome. We assessed the association of microbiome diversity with norovirus infection and secretor status in patients from Veterans Affairs medical centers. Alpha diversity metrics were lower among patients with acute gastroenteritis but were similar for other comparisons.

Species-Wide Phylogenomics of the Staphylococcus aureus Agr Operon Revealed Convergent Evolution of Frameshift Mutations

Staphylococcus aureus is a prominent nosocomial pathogen that causes several life-threatening diseases, such as pneumonia and bacteremia. S. aureus modulates the expression of its arsenal of virulence factors through sensing and integrating responses to environmental signals. The agr (accessory gene regulator) quorum sensing (QS) system is a major regulator of virulence phenotypes in S. aureus. There are four agr specificity groups each with a different autoinducer peptide sequence encoded by the agrD gene. Although agr is critical for the expression of many toxins, paradoxically, S. aureus strains often have nonfunctional agr activity due to loss-of-function mutations in the four-gene agr operon. To understand patterns in agr variability across S. aureus, we undertook a species-wide genomic investigation. We developed a software tool (AgrVATE; https://github.com/VishnuRaghuram94/AgrVATE) for typing and detecting frameshift mutations in the agr operon. In an analysis of over 40,000 S. aureus genomes, we showed a close association between agr type and S. aureus clonal complex. We also found a strong linkage between agrBDC alleles (encoding the peptidase, autoinducing peptide itself, and peptide sensor, respectively) but not agrA (encoding the response regulator). More than 5% of the genomes were found to have frameshift mutations in the agr operon. While 52% of these frameshifts occurred only once in the entire species, we observed cases where the recurring mutations evolved convergently across different clonal lineages with no evidence of long-term phylogenetic transmission, suggesting that strains with agr frameshifts were evolutionarily short-lived. Overall, genomic analysis of agr operon suggests evolution through multiple processes with functional consequences that are not fully understood. IMPORTANCE Staphylococcus aureus is a globally pervasive pathogen that produces a plethora of toxic molecules that can harm host immune cells. Production of these toxins is mainly controlled by an active agr quorum-sensing system, which senses and responds to bacterial cell density. However, there are many reports of S. aureus strains with genetic changes leading to impaired agr activity that are often found during chronic bloodstream infections and may be associated with increased disease severity. We developed an open-source software called AgrVATE to type agr systems and identify mutations. We used AgrVATE for a species-wide genomic survey of S. aureus, finding that more than 5% of strains in the public database had nonfunctional agr systems. We also provided new insights into the evolution of these genetic mutations in the agr system. Overall, this study contributes to our understanding of a common but relatively understudied means of virulence regulation in S. aureus.

Unrecognized introductions of SARS-CoV-2 into the US state of Georgia shaped the early epidemic

In early 2020, as diagnostic and surveillance responses for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) ramped up, attention focused primarily on returning international travelers. Here, we build on existing studies characterizing early patterns of SARS-CoV-2 spread within the USA by analyzing detailed clinical, molecular, and viral genomic data from the state of Georgia through March 2020. We find evidence for multiple early introductions into Georgia, despite relatively sparse sampling. Most sampled sequences likely stemmed from a single or small number of introductions from Asia three weeks prior to the state's first detected infection. Our analysis of sequences from domestic travelers demonstrates widespread circulation of closely related viruses in multiple US states by the end of March 2020. Our findings indicate that the exclusive focus on identifying SARS-CoV-2 in returning international travelers early in the pandemic may have led to a failure to recognize locally circulating infections for several weeks and point toward a critical need for implementing rapid, broadly targeted surveillance efforts for future pandemics.

Dynamic PET-facilitated modeling and high-dose rifampin regimens for Staphylococcus aureus orthopedic implant-associated infections

Staphylococcus aureus is a major human pathogen causing serious implant–associated infections. Combination treatment with rifampin (10 to 15 mg/kg per day), which has dose-dependent activity, is recommended to treat S. aureus orthopedic implant–associated infections. Rifampin, however, has limited bone penetration. Here, dynamic 11C-rifampin positron emission tomography (PET) performed in prospectively enrolled patients with confirmed S. aureus bone infection (n = 3) or without orthopedic infection (n = 12) demonstrated bone/plasma area under the concentration-time curve ratio of 0.14 (interquartile range, 0.09 to 0.19), exposures lower than previously thought. PET-based pharmacokinetic modeling predicted rifampin concentration-time profiles in bone and facilitated studies in a mouse model of S. aureus orthopedic implant infection. Administration of high-dose rifampin (human equipotent to 35 mg/kg per day) substantially increased bone concentrations (2 mg/liter versus <0.2 mg/liter with standard dosing) in mice and achieved higher bacterial killing and biofilm disruption. Treatment for 4 weeks with high-dose rifampin and vancomycin was noninferior to the recommended 6-week treatment of standard-dose rifampin with vancomycin in mice (risk difference, −6.7% favoring high-dose rifampin regimen). High-dose rifampin treatment ameliorated antimicrobial resistance (0% versus 38%; P = 0.04) and mitigated adverse bone remodeling (P < 0.01). Last, whole-genome sequencing demonstrated that administration of high-dose rifampin in mice reduced selection of bacterial mutations conferring rifampin resistance (rpoB) and mutations in genes potentially linked to persistence. These data suggest that administration of high-dose rifampin is necessary to achieve optimal bone concentrations, which could shorten and improve treatments for S. aureus orthopedic implant infections.

Comparing whole-genome shotgun sequencing and DNA metabarcoding approaches for species identification and quantification of pollen species mixtures

Molecular identification of mixed-species pollen samples has a range of applications in various fields of research. To date, such molecular identification has primarily been carried out via amplicon sequencing, but whole-genome shotgun (WGS) sequencing of pollen DNA has potential advantages, including (1) more genetic information per sample and (2) the potential for better quantitative matching. In this study, we tested the performance of WGS sequencing methodology and publicly available reference sequences in identifying species and quantifying their relative abundance in pollen mock communities. Using mock communities previously analyzed with DNA metabarcoding, we sequenced approximately 200Mbp for each sample using Illumina HiSeq and MiSeq. Taxonomic identifications were based on the Kraken k-mer identification method with reference libraries constructed from full-genome and short read archive data from the NCBI database. We found WGS to be a reliable method for taxonomic identification of pollen with near 100% identification of species in mixtures but generating higher rates of false positives (reads not identified to the correct taxon at the required taxonomic level) relative to rbcL and ITS2 amplicon sequencing. For quantification of relative species abundance, WGS data provided a stronger correlation between pollen grain proportion and sequence read proportion, but diverged more from a 1:1 relationship, likely due to the higher rate of false positives. Currently, a limitation of WGS-based pollen identification is the lack of representation of plant diversity in publicly available genome databases. As databases improve and costs drop, we expect that eventually genomics methods will become the methods of choice for species identification and quantification of mixed-species pollen samples.

Effect of genetic background on the evolution of Vancomycin-Intermediate Staphylococcus aureus (VISA)

Vancomycin-intermediate Staphylococcus aureus (VISA) typically arises through accumulation of chromosomal mutations that alter cell-wall thickness and global regulatory pathways. Genome-based prediction of VISA requires understanding whether strain background influences patterns of mutation that lead to resistance. We used an iterative method to experimentally evolve three important methicillin-resistant S. aureus (MRSA) strain backgrounds-(CC1, CC5 and CC8 (USA300)) to generate a library of 120 laboratory selected VISA isolates. At the endpoint, isolates had vancomycin MICs ranging from 4 to 10 μg/mL. We detected mutations in more than 150 genes, but only six genes (already known to be associated with VISA from prior studies) were mutated in all three background strains (walK, prs, rpoB, rpoC, vraS, yvqF). We found evidence of interactions between loci (e.g., vraS and yvqF mutants were significantly negatively correlated) and rpoB, rpoC, vraS and yvqF were more frequently mutated in one of the backgrounds. Increasing vancomycin resistance was correlated with lower maximal growth rates (a proxy for fitness) regardless of background. However, CC5 VISA isolates had higher MICs with fewer rounds of selection and had lower fitness costs than the CC8 VISA isolates. Using multivariable regression, we found that genes differed in their contribution to overall MIC depending on the background. Overall, these results demonstrated that VISA evolved through mutations in a similar set of loci in all backgrounds, but the effect of mutation in common genes differed with regard to fitness and contribution to resistance in different strains.

Generation of Tetracycline and Rifamycin Resistant Chlamydia Suis Recombinants

The Chlamydiaceae are a family of obligate intracellular, gram-negative bacteria known to readily exchange DNA by homologous recombination upon co-culture in vitro, allowing the transfer of antibiotic resistance residing on the chlamydial chromosome. Among all the obligate intracellular bacteria, only Chlamydia (C.) suis naturally integrated a tetracycline resistance gene into its chromosome. Therefore, in order to further investigate the readiness of Chlamydia to exchange DNA and especially antibiotic resistance, C. suis is an excellent model to advance existing co-culture protocols allowing the identification of factors crucial to promote homologous recombination in vitro. With this strategy, we co-cultured tetracycline-resistant with rifamycin group-resistant C. suis, which resulted in an allover recombination efficiency of 28%. We found that simultaneous selection is crucial to increase the number of recombinants, that sub-inhibitory concentrations of tetracycline inhibit rather than promote the selection of double-resistant recombinants, and identified a recombination-deficient C. suis field isolate, strain SWA-110 (1-28b). While tetracycline resistance was detected in field isolates, rifampicin/rifamycin resistance (RifR) had to be induced in vitro. Here, we describe the protocol with which RifR C. suis strains were generated and confirmed. Subsequent whole-genome sequencing then revealed that G530E and D461A mutations in rpoB, a gene encoding for the β-subunit of the bacterial RNA polymerase (RNAP), was likely responsible for rifampicin and rifamycin resistance, respectively. Finally, whole-genome sequencing of recombinants obtained by co-culture revealed that recombinants picked from the same plate may be sibling clones and confirmed C. suis genome plasticity by revealing variable, apparently non-specific areas of recombination.

Symbiont Genomic Features and Localization in the Bean Beetle Callosobruchus maculatus

A pervasive pest of stored leguminous products, the bean beetle Callosobruchus maculatus (Coleoptera: Chrysomelidae) associates with a simple bacterial community during adulthood. Despite its economic importance, little is known about the compositional stability, heritability, localization, and metabolic potential of the bacterial symbionts of C. maculatus. In this study, we applied community profiling using 16S rRNA gene sequencing to reveal a highly conserved bacterial assembly shared between larvae and adults. Dominated by Firmicutes and Proteobacteria, this community is localized extracellularly along the epithelial lining of the bean beetle's digestive tract. Our analysis revealed that only one species, Staphylococcus gallinarum (phylum Firmicutes), is shared across all developmental stages. Isolation and whole-genome sequencing of S. gallinarum from the beetle gut yielded a circular chromosome (2.8 Mb) and one plasmid (45 kb). The strain encodes complete biosynthetic pathways for the production of B vitamins and amino acids, including tyrosine, which is increasingly recognized as an important symbiont-supplemented precursor for cuticle biosynthesis in beetles. A carbohydrate-active enzyme search revealed that the genome codes for a number of digestive enzymes, reflecting the nutritional ecology of C. maculatus. The ontogenic conservation of the gut microbiota in the bean beetle, featuring a "core" community composed of S. gallinarum, may be indicative of an adaptive role for the host. In clarifying symbiont localization and metabolic potential, we further our understanding and study of a costly pest of stored products. IMPORTANCE From supplementing essential nutrients to detoxifying plant secondary metabolites and insecticides, bacterial symbionts are a key source of adaptations for herbivorous insect pests. Despite the pervasiveness and geographical range of the bean beetle Callosobruchus maculatus, the role of microbial symbioses in its natural history remains understudied. Here, we demonstrate that the bean beetle harbors a simple gut bacterial community that is stable throughout development. This community localizes along the insect's digestive tract and is largely dominated by Staphylococcus gallinarum. In elucidating symbiont metabolic potential, we highlight its possible adaptive significance for a widespread agricultural pest.

Shiftwork, functional bowel symptoms, and the microbiome

BACKGROUND

There are about 15 million Americans working full-time on evening, night, or rotating shifts. Between 48% and 81.9% of those working rotating or night shifts report abdominal pain, constipation, diarrhea and other symptoms of functional bowel disorders. The basis for this high prevalence of functional bowel disorders, including irritable bowel syndrome (IBS), among shift workers is unknown. Animal studies, however, suggest that circadian disruption, similar to that in shift workers, may contribute to the development of GI complaints among shift workers by altering the composition and normal diurnal rhythmicity of the resident intestinal microbes. Therefore, the present study was designed to determine if there were differences in (1) composition and diversity of the microbiome of night shift workers compared to day shift workers; and (2) the composition and diversity of the microbiome among shift workers experiencing functional bowel symptoms compared to shift workers who did not experience functional bowel symptoms.

METHODS

Fifty-one full time staff nurses who worked either 12-hour day or night shifts completed demographic information, and the Rome III IBS module. They also collected two samples of gut microbiota before the beginning and at the end of their last work shift on day 14, using validated field-tested methods consistent with the Human Microbiome Project. After DNA extraction, 16S rRNA sequencing and assignment to the genus level was completed, samples were then compared to determine if there were (1) differences in the diversity and profile of the microbiome by shift type; (2) if there were differences in the microbiome by time of day for collection; and (3) whether there were differences in the diversity and profile of the microbiome of nurses with IBS and those without IBS.

RESULTS

There were no differences in alpha or beta diversity of gut microbiota when specimens from day and night shift nurses were compared. There were however marginal differences in beta diversity when specimens collected at the beginning and end of the shifts were compared, with seven OTUs being differentially abundant when collected from day shift workers in the evening. There were also three OTUs to be differentially abundant in participants reporting IBS symptoms.

Vaginal Microbiome Composition in Early Pregnancy and Risk of Spontaneous Preterm and Early Term Birth Among African American Women

Objective

To evaluate the association between the early pregnancy vaginal microbiome and spontaneous preterm birth (sPTB) and early term birth (sETB) among African American women.

Methods

Vaginal samples collected in early pregnancy (8-14 weeks' gestation) from 436 women enrolled in the Emory University African American Vaginal, Oral, and Gut Microbiome in Pregnancy Study underwent 16S rRNA gene sequencing of the V3-V4 region, taxonomic classification, and community state type (CST) assignment. We compared vaginal CST and abundance of taxa for women whose pregnancy ended in sPTB (N = 44) or sETB (N= 84) to those who delivered full term (N = 231).

Results

Nearly half of the women had a vaginal microbiome classified as CST IV (Diverse CST), while one-third had CST III (L. iners dominated) and just 16% had CST I, II, or V (non-iners Lactobacillus dominated). Compared to vaginal CST I, II, or V (non-iners Lactobacillus dominated), both CST III (L. iners dominated) and CST IV (Diverse) were associated with sPTB with an adjusted odds ratio (95% confidence interval) of 4.1 (1.1, infinity) and 7.7 (2.2, infinity), respectively, in multivariate logistic regression. In contrast, no vaginal CST was associated with sETB. The linear decomposition model (LDM) based on amplicon sequence variant (ASV) relative abundance found a significant overall effect of the vaginal microbiome on sPTB (p=0.034) but not sETB (p=0.320), whereas the LDM based on presence/absence of ASV found no overall effect on sPTB (p=0.328) but a significant effect on sETB (p=0.030). In testing for ASV-specific effects, the LDM found that no ASV was significantly associated with sPTB considering either relative abundance or presence/absence data after controlling for multiple comparisons (FDR 10%), although in marginal analysis the relative abundance of Gardnerella vaginalis (p=0.011), non-iners Lactobacillus (p=0.016), and Mobiluncus curtisii (p=0.035) and the presence of Atopobium vaginae (p=0.049), BVAB2 (p=0.024), Dialister microaerophilis (p=0.011), and Prevotella amnii (p=0.044) were associated with sPTB. The LDM identified the higher abundance of 7 ASVs and the presence of 13 ASVs, all commonly residents of the gut, as associated with sETB at FDR < 10%.

Conclusions

In this cohort of African American women, an early pregnancy vaginal CST III or IV was associated with an increased risk of sPTB but not sETB. The relative abundance and presence of distinct taxa within the early pregnancy vaginal microbiome was associated with either sPTB or sETB.

Genes Influencing Phage Host Range in Staphylococcus aureus on a Species-Wide Scale

Staphylococcus aureus is a human pathogen that causes serious diseases, ranging from skin infections to septic shock. Bacteriophages (phages) are both natural killers of S. aureus, offering therapeutic possibilities, and important vectors of horizontal gene transfer (HGT) in the species. Here, we used high-throughput approaches to understand the genetic basis of strain-to-strain variation in sensitivity to phages, which defines the host range. We screened 259 diverse S. aureus strains covering more than 40 sequence types for sensitivity to eight phages, which were representatives of the three phage classes that infect the species. The phages were variable in host range, each infecting between 73 and 257 strains. Using genome-wide association approaches, we identified putative loci that affect host range and validated their function using USA300 transposon knockouts. In addition to rediscovering known host range determinants, we found several previously unreported genes affecting bacterial growth during phage infection, including trpA, phoR, isdB, sodM, fmtC, and relA We used the data from our host range matrix to develop predictive models that achieved between 40% and 95% accuracy. This work illustrates the complexity of the genetic basis for phage susceptibility in S. aureus but also shows that with more data, we may be able to understand much of the variation. With a knowledge of host range determination, we can rationally design phage therapy cocktails that target the broadest host range of S. aureus strains and address basic questions regarding phage-host interactions, such as the impact of phage on S. aureus evolution.IMPORTANCEStaphylococcus aureus is a widespread, hospital- and community-acquired pathogen, many strains of which are antibiotic resistant. It causes diverse diseases, ranging from local to systemic infection, and affects both the skin and many internal organs, including the heart, lungs, bones, and brain. Its ubiquity, antibiotic resistance, and disease burden make new therapies urgent. One alternative therapy to antibiotics is phage therapy, in which viruses specific to infecting bacteria clear infection. In this work, we identified and validated S. aureus genes that influence phage host range-the number of strains a phage can infect and kill-by testing strains representative of the diversity of the S. aureus species for phage host range and associating the genome sequences of strains with host range. These findings together improved our understanding of how phage therapy works in the bacterium and improve prediction of phage therapy efficacy based on the predicted host range of the infecting strain.

852. Genomic Clusters of Methicillin-Resistant Staphylococcus aureus (MRSA) Causing Bloodstream Infections (BSIs) in Hospitalized Adults, 2018-19

Background

MRSA BSIs have 15-50% mortality and are commonly diagnosed in US hospitals. However, the frequency of hospital transmission of MRSA causing BSI is unknown.

Methods

We performed Illumina shotgun whole genome sequencing (WGS) of 106 sequential MRSA isolates from different adults with a BSI at two Philadelphia academic hospitals in a single health system in July 2018-June 2019. We abstracted clinical data from the electronic medical record. Genomic data were analyzed preliminarily using the Staphopia Analysis Pipeline.

Results

Among 106 subjects, 51.9% were male, 47.2% were white, 46.2% were black, 23.6% were < 40 years of age, and mean age was 53.1 years (s.d. 17 years). One isolate had WGS data that were inadequate for analysis. Of 105 genomes, 52 were clonal cluster (CC) 8, 22 were sequence type (ST) 5, and 16 were ST105; the remaining 15 strains belonged to 8 other CCs. Of CC8 strains, 44 were USA300 and 6 were USA500. There were 6 clusters (i.e., < 35 SNP differences in the core genome) among the 105 isolates. Four clusters were CC5 and two were CC8 strains. One cluster of CC5 strains involved 3 subjects, and 5 clusters involved 2 subjects. One cluster of ST8/USA300 strains were separated by only 1 SNP (Fig a). This and two other clustered pairs were from subjects who had overlapping hospital stays. Two of these paired subjects had an overlap in the same unit while the third pair was in the hospital together on a number of occasions (total of 40 days overlap) but never simultaneously in the same unit. The other three clustered pairs did not have temporally overlapping hospital stays, suggesting transmission via a hospital reservoir. One of these three pairs had hospitalizations overlapping in time, one at each study hospital, before each of them had infections with the related MRSA strains. There was not a clear-cut clustering of SNP distances among the isolate genomes into transmission and non-transmission groups, with some pairs of patient isolates separated by 40-80 SNPs (Fig. b).

Figure 1.

Conclusion

We were able to discern from WGS data alone that some MRSA BSIs in 2 hospitals were likely due to strains transmitted between patients. Universal WGS of BSI strains may detect MRSA outbreaks in real time, even in the absence of overlapping hospitalizations, and is an emerging strategy to detect healthcare transmission of MRSA.

Disclosures

Michael Z. David, MD PhD, GSK (Consultant)

Bactopia: a Flexible Pipeline for Complete Analysis of Bacterial Genomes

Sequencing of bacterial genomes using Illumina technology has become such a standard procedure that often data are generated faster than can be conveniently analyzed. We created a new series of pipelines called Bactopia, built using Nextflow workflow software, to provide efficient comparative genomic analyses for bacterial species or genera. Bactopia consists of a data set setup step (Bactopia Data Sets [BaDs]), which creates a series of customizable data sets for the species of interest, the Bactopia Analysis Pipeline (BaAP), which performs quality control, genome assembly, and several other functions based on the available data sets and outputs the processed data to a structured directory format, and a series of Bactopia Tools (BaTs) that perform specific postprocessing on some or all of the processed data. BaTs include pan-genome analysis, computing average nucleotide identity between samples, extracting and profiling the 16S genes, and taxonomic classification using highly conserved genes. It is expected that the number of BaTs will increase to fill specific applications in the future. As a demonstration, we performed an analysis of 1,664 public Lactobacillus genomes, focusing on Lactobacillus crispatus, a species that is a common part of the human vaginal microbiome. Bactopia is an open source system that can scale from projects as small as one bacterial genome to ones including thousands of genomes and that allows for great flexibility in choosing comparison data sets and options for downstream analysis. Bactopia code can be accessed at https://www.github.com/bactopia/bactopiaIMPORTANCE It is now relatively easy to obtain a high-quality draft genome sequence of a bacterium, but bioinformatic analysis requires organization and optimization of multiple open source software tools. We present Bactopia, a pipeline for bacterial genome analysis, as an option for processing bacterial genome data. Bactopia also automates downloading of data from multiple public sources and species-specific customization. Because the pipeline is written in the Nextflow language, analyses can be scaled from individual genomes on a local computer to thousands of genomes using cloud resources. As a usage example, we processed 1,664 Lactobacillus genomes from public sources and used comparative analysis workflows (Bactopia Tools) to identify and analyze members of the L. crispatus species.

Glucocorticoid receptor sensitivity in early pregnancy in an African American cohort

PROBLEM

Disruption in homeostatic feedback loops between inflammatory mediators and the hypothalamic-pituitary-adrenal (HPA) axis is a key mechanism linking chronic stress to inflammation and adverse health outcomes, including those occurring during pregnancy. In particular, alterations in glucocorticoid sensitivity may occur as a result of chronic stress, including that due to racial discrimination, and may be implicated in the persistent adverse maternal and infant health outcomes experienced by African Americans. While there are a few large-scale studies in human pregnancy that measure both cytokines and HPA axis hormones, to our knowledge, none directly measure glucocorticoid sensitivity at the cellular level, especially in an African American population.

METHOD OF STUDY

We measured the full range of the dexamethasone (DEX) dose-response suppression of TNF-α in first-trimester blood samples from 408 African American women and estimated leukocyte cell type contribution to the production of TNF-α.

RESULTS

The mean (SD) DEX level needed to inhibit TNF-α production by 50% (ie, DEX IC50 ) was 9.8 (5.8) nmol/L. Monocytes appeared to be the main driver of Uninhibited TNF-α production, but monocyte counts explained only 14% of the variation. Monocyte counts were only weakly correlated with the DEX IC50 (r = -.11, P < .05). Moreover, there was no statistically significant correlation between the DEX IC50 and circulating pro-inflammatory (CRP, IL-6, IFN-γ) or anti-inflammatory (IL-10) mediators (P > .05).

CONCLUSION

These findings challenge some prior assumptions and position this comprehensive study of glucocorticoid sensitivity as an important anchor point in the growing recognition of interindividual variation in maternal HPA axis regulation and inflammatory responses.

Genotypic and Phenotypic Diversity of Staphylococcus aureus Isolates from Cystic Fibrosis Patient Lung Infections and Their Interactions with Pseudomonas aeruginosa

Staphylococcus aureus has recently overtaken Pseudomonas aeruginosa as the most commonly recognized bacterial pathogen that infects the respiratory tracts of individuals with the genetic disease cystic fibrosis (CF) in the United States. Most studies of S. aureus in CF patient lung infections have focused on a few isolates, often exclusively laboratory-adapted strains, and how they are killed by P. aeruginosa Less is known about the diversity of S. aureus CF patient lung isolates in terms of both their virulence and their interaction with P. aeruginosa To begin to address this gap, we recently sequenced 64 clinical S. aureus isolates and a reference isolate, JE2. Here, we analyzed the antibiotic resistance genotypes, sequence types, clonal complexes, spa types, agr types, and presence/absence of other known virulence factor genes of these isolates. We hypothesized that virulence phenotypes of S. aureus, namely, toxin production and the mucoid phenotype, would be lost in these isolates due to adaptation in the CF patient lung. In contrast to these expectations, we found that most isolates can lyse both rabbit and sheep blood (67.7%) and produce polysaccharide (69.2%), suggesting that these phenotypes were not lost during adaptation to the CF lung. We also identified three distinct phenotypic groups of S. aureus based on their survival in the presence of nonmucoid P. aeruginosa laboratory strain PAO1 and its mucoid derivative. Altogether, our work provides greater insight into the diversity of S. aureus isolates from CF patients, specifically the distribution of important virulence factors and their interaction with P. aeruginosa, all of which have implications in patient health.IMPORTANCEStaphylococcus aureus is now the most frequently detected recognized pathogen in the lungs of individuals who have cystic fibrosis (CF) in the United States, followed closely by Pseudomonas aeruginosa When these pathogens are found to coinfect the CF lung, patients have a significantly worse prognosis. While P. aeruginosa has been rigorously studied in the context of bacterial pathogenesis in CF, less is known about S. aureus Here, we present an in-depth study of 64 S. aureus clinical isolates from CF patients, for which we investigated genetic diversity utilizing whole-genome sequencing, virulence phenotypes, and interactions with P. aeruginosa We found that S. aureus isolated from CF lungs are phylogenetically diverse; most retain known virulence factors and vary in their interactions with P. aeruginosa (i.e., they range from being highly sensitive to P. aeruginosa to completely tolerant to it). Deepening our understanding of how S. aureus responds to its environment and other microbes in the CF lung will enable future development of effective treatments and preventative measures against these formidable infections.

Genotypic and Phenotypic Diversity of Staphylococcus aureus Isolates from Cystic Fibrosis Patient Lung Infections and Their Interactions with Pseudomonas aeruginosa

Staphylococcus aureus has recently overtaken Pseudomonas aeruginosa as the most commonly recognized bacterial pathogen that infects the respiratory tracts of individuals with the genetic disease cystic fibrosis (CF) in the United States. Most studies of S. aureus in CF patient lung infections have focused on a few isolates, often exclusively laboratory-adapted strains, and how they are killed by P. aeruginosa Less is known about the diversity of S. aureus CF patient lung isolates in terms of both their virulence and their interaction with P. aeruginosa To begin to address this gap, we recently sequenced 64 clinical S. aureus isolates and a reference isolate, JE2. Here, we analyzed the antibiotic resistance genotypes, sequence types, clonal complexes, spa types, agr types, and presence/absence of other known virulence factor genes of these isolates. We hypothesized that virulence phenotypes of S. aureus, namely, toxin production and the mucoid phenotype, would be lost in these isolates due to adaptation in the CF patient lung. In contrast to these expectations, we found that most isolates can lyse both rabbit and sheep blood (67.7%) and produce polysaccharide (69.2%), suggesting that these phenotypes were not lost during adaptation to the CF lung. We also identified three distinct phenotypic groups of S. aureus based on their survival in the presence of nonmucoid P. aeruginosa laboratory strain PAO1 and its mucoid derivative. Altogether, our work provides greater insight into the diversity of S. aureus isolates from CF patients, specifically the distribution of important virulence factors and their interaction with P. aeruginosa, all of which have implications in patient health.IMPORTANCEStaphylococcus aureus is now the most frequently detected recognized pathogen in the lungs of individuals who have cystic fibrosis (CF) in the United States, followed closely by Pseudomonas aeruginosa When these pathogens are found to coinfect the CF lung, patients have a significantly worse prognosis. While P. aeruginosa has been rigorously studied in the context of bacterial pathogenesis in CF, less is known about S. aureus Here, we present an in-depth study of 64 S. aureus clinical isolates from CF patients, for which we investigated genetic diversity utilizing whole-genome sequencing, virulence phenotypes, and interactions with P. aeruginosa We found that S. aureus isolated from CF lungs are phylogenetically diverse; most retain known virulence factors and vary in their interactions with P. aeruginosa (i.e., they range from being highly sensitive to P. aeruginosa to completely tolerant to it). Deepening our understanding of how S. aureus responds to its environment and other microbes in the CF lung will enable future development of effective treatments and preventative measures against these formidable infections.

Genomic analysis of variability in Delta-toxin levels between Staphylococcus aureus strains

BACKGROUND

The delta-toxin (δ-toxin) of Staphylococcus aureus is the only hemolysin shown to cause mast cell degranulation and is linked to atopic dermatitis, a chronic inflammatory skin disease. We sought to characterize variation in δ-toxin production across S. aureus strains and identify genetic loci potentially associated with differences between strains.

METHODS

A set of 124 S. aureus strains was genome-sequenced and δ-toxin levels in stationary phase supernatants determined by high performance liquid chromatography (HPLC). SNPs and kmers were associated with differences in toxin production using four genome-wide association study (GWAS) methods. Transposon mutations in candidate genes were tested for their δ-toxin levels. We constructed XGBoost models to predict toxin production based on genetic loci discovered to be potentially associated with the phenotype.

RESULTS

The S. aureus strain set encompassed 40 sequence types (STs) in 23 clonal complexes (CCs). δ-toxin production ranged from barely detectable levels to >90,000 units, with a median of >8,000 units. CC30 had significantly lower levels of toxin production than average while CC45 and CC121 were higher. MSSA (methicillin sensitive) strains had higher δ-toxin production than MRSA (methicillin resistant) strains. Through multiple GWAS approaches, 45 genes were found to be potentially associated with toxicity. Machine learning models using loci discovered through GWAS as features were able to predict δ-toxin production (as a high/low binary phenotype) with a precision of .875 and specificity of .990 but recall of .333. We discovered that mutants in the carA gene, encoding the small chain of carbamoyl phosphate synthase, completely abolished toxin production and toxicity in Caenorhabditis elegans.

CONCLUSIONS

The amount of stationary phase production of the toxin is a strain-specific phenotype likely affected by a complex interaction of number of genes with different levels of effect. We discovered new candidate genes that potentially play a role in modulating production. We report for the first time that the product of the carA gene is necessary for δ-toxin production in USA300. This work lays a foundation for future work on understanding toxin regulation in S. aureus and prediction of phenotypes from genomic sequences.

A Comparison of the Bacterial Nasal Microbiome in Allergic Rhinitis Patients Before and After Immunotherapy

OBJECTIVES/HYPOTHESIS

Recent research has examined the nasal microbiome in rhinosinusitis and nondiseased states. Given immunologic alterations in allergic rhinitis (AR) and after allergen immunotherapy (IT), we evaluated the nasal microbiome in these conditions.

STUDY DESIGN

Cross-sectional comparison.

METHODS

In this cross-sectional study, nasal swabs for microbiome analysis were collected from three patient groups: IT-naïve AR patients, AR patients undergoing IT for greater than 12 months, and a control group without sinonasal inflammatory disease.

RESULTS

Nasal swabs were successfully collected for 14 IT-naïve AR patients, 20 post-IT patients, and 17 controls. The α diversity showed a statistical difference in evenness but not in richness amongst samples, whereas the β-diversity was significantly different between groups. Corynebacterium and Staphylococcus were the most prevalent bacteria across all groups.

CONCLUSIONS

β-diversity was found to be significantly different across the three groups, but the AR groups were found to be more similar to each other than to the controls. Although there is symptomatic improvement in the AR group undergoing IT, the microbiome does not appear to transition to a healthy microbiome composition.

LEVEL OF EVIDENCE

4 Laryngoscope, 2020.

Human mAbs to Staphylococcus aureus IsdA Provide Protection Through Both Heme-Blocking and Fc-Mediated Mechanisms

The nutrient metal iron plays a key role in the survival of microorganisms. The iron-regulated surface determinant (Isd) system scavenges heme-iron from the human host, enabling acquisition of iron in iron-deplete conditions in Staphylococcus aureus during infection. The cell surface receptors IsdB and IsdH bind hemoproteins and transfer heme to IsdA, the final surface protein before heme-iron is transported through the peptidoglycan. To define the human B-cell response to IsdA, we isolated human monoclonal antibodies (mAbs) specific to the surface Isd proteins and determined their mechanism of action. We describe the first isolation of fully human IsdA and IsdH mAbs, as well as cross-reactive Isd mAbs. Two of the identified IsdA mAbs worked in a murine septic model of infection to reduce bacterial burden during staphylococcal infection. Their protection was a result of both heme-blocking and Fc-mediated effector functions, underscoring the importance of targeting S. aureus using diverse mechanisms.

Invasive Nontypeable Haemophilus influenzae Infection Among Adults With HIV in Metropolitan Atlanta, Georgia, 2008-2018

IMPORTANCE

Invasive nontypeable Haemophilus influenzae (NTHi) infection among adults is typically associated with bacteremic pneumonia. Nontypeable H influenzae is genetically diverse and clusters of infection are uncommon.

OBJECTIVE

To evaluate an increase in invasive NTHi infection from 2017-2018 among HIV-infected men who have sex with men in metropolitan Atlanta, Georgia.

DESIGN, SETTING, AND PARTICIPANTS

A population-based surveillance study with a cohort substudy and descriptive epidemiological analysis identified adults aged 18 years or older with invasive NTHi infection (isolation of NTHi from a normally sterile site) between January 1, 2008, and December 31, 2018 (final date of follow-up).

EXPOSURES

Time period, HIV status, and genetic relatedness (ie, cluster status) of available NTHi isolates.

MAIN OUTCOMES AND MEASURES

The primary outcome was incidence of invasive NTHi infection (from 2008-2016 and 2017-2018) among persons with HIV and compared with NTHi infection from 2008-2018 among those without HIV. The secondary outcomes were assessed among those aged 18 to 55 years with invasive NTHi infection and included epidemiological, clinical, and geographic comparisons by cluster status.

RESULTS

Among 553 adults with invasive NTHi infection (median age, 66 years [Q1-Q3, 48-78 years]; 52% male; and 38% black), 60 cases occurred among persons with HIV. Incidence of invasive NTHi infection from 2017-2018 among persons with HIV (41.7 cases per 100 000) was significantly greater than from 2008-2016 among those with HIV (9.6 per 100 000; P < .001) and from 2008-2018 among those without HIV (1.1 per 100 000; P < .001). Among adults aged 18 to 55 years with invasive NTHi infections from 2017-2018 (n = 179), persons with HIV (n = 31) were significantly more likely than those from 2008-2018 without HIV (n = 124) to be male (94% vs 49%, respectively; P < .001), black (100% vs 53%; P < .001), and have septic arthritis (35% vs 1%; P < .001). Persons with HIV who had invasive NTHi infection from 2017-2018 (n = 31) were more likely than persons with HIV who had invasive NTHi infection from 2008-2016 (n = 24) to have septic arthritis (35% vs 4%, respectively; P = .01). Pulsed-field gel electrophoresis of 174 of 179 NTHi isolates from 18- to 55-year-olds identified 2 genetically distinct clonal groups: cluster 1 (C1; n = 24) and cluster 2 (C2; n = 23). Whole-genome sequencing confirmed 2 clonal lineages of NTHi infection and revealed all C1 isolates (but none of the C2 isolates) carried IS1016 (an insertion sequence associated with H influenzae capsule genes). Persons with HIV were significantly more likely to have C1 or C2 invasive NTHi infection from 2017-2018 (28/31 [90%]) compared with from 2008-2016 among persons with HIV (10/24 [42%]; P < .001) and compared with from 2008-2018 among those without HIV (9/119 [8%]; P < .001). Among persons with C1 or C2 invasive NTHi infection who had HIV (n = 38) (median age, 34.5 years; 100% male; 100% black; 82% men who have sex with men), 32 (84%) lived in 2 urban counties and an area of significant spatial aggregation was identified compared with those without C1 or C2 invasive NTHi infection.

CONCLUSIONS AND RELEVANCE

Among persons with HIV in Atlanta, the incidence of invasive nontypeable H influenzae infection increased significantly from 2017-2018 compared with 2008-2016. Two unique but genetically related clonal strains were identified and were associated with septic arthritis among black men who have sex with men and who lived in geographic proximity.

Stability of the vaginal, oral, and gut microbiota across pregnancy among African American women: the effect of socioeconomic status and antibiotic exposure

OBJECTIVE

A growing body of research has investigated the human microbiota and pregnancy outcomes, especially preterm birth. Most studies of the prenatal microbiota have focused on the vagina, with fewer investigating other body sites during pregnancy. Although pregnancy involves profound hormonal, immunological and metabolic changes, few studies have investigated either shifts in microbiota composition across pregnancy at different body sites or variation in composition at any site that may be explained by maternal characteristics. The purpose of this study was to investigate: (1) the stability of the vaginal, oral, and gut microbiota from early (8-14 weeks) through later (24-30 weeks) pregnancy among African American women according to measures of socioeconomic status, accounting for prenatal antibiotic use; (2) whether measures of socioeconomic status are associated with changes in microbiota composition over pregnancy; and (3) whether exposure to prenatal antibiotics mediate any observed associations between measures of socioeconomic status and stability of the vaginal, oral, and gut microbiota across pregnancy.

METHODS

We used paired vaginal, oral, or gut samples available for 16S rRNA gene sequencing from two time points in pregnancy (8-14 and 24-30 weeks) to compare within-woman changes in measures of alpha diversity (Shannon and Chao1) and beta-diversity (Bray-Curtis dissimilarity) among pregnant African American women (n = 110). Multivariable linear regression was used to examine the effect of level of education and prenatal health insurance as explanatory variables for changes in diversity, considering antibiotic exposure as a mediator, adjusting for age, obstetrical history, and weeks between sampling.

RESULTS

For the oral and gut microbiota, there were no significant associations between measures of socioeconomic status or prenatal antibiotic use and change in Shannon or Chao1 diversity. For the vaginal microbiota, low level of education (high school or less) was associated with an increase in Shannon and Chao1 diversity over pregnancy, with minimal attenuation when controlling for prenatal antibiotic use. Conversely, for within-woman Bray-Curtis dissimilarity for early compared to later pregnancy, low level of education and prenatal antibiotics were associated with greater dissimilarity for the oral and gut sites, with minimal attenuation when controlling for prenatal antibiotics, and no difference in dissimilarity for the vaginal site.

CONCLUSIONS

Measures of maternal socioeconomic status are variably associated with changes in diversity across pregnancy for the vaginal, oral, and gut microbiota, with minimal attenuation by prenatal antibiotic exposure. Studies that evaluate stability of the microbiota across pregnancy in association with health outcomes themselves associated with socioeconomic status (such as preterm birth) should incorporate measures of socioeconomic status to avoid finding spurious relationships.

Determinants of Phage Host Range in Staphylococcus Species

Bacteria in the genus Staphylococcus are important targets for phage therapy due to their prevalence as pathogens and increasing antibiotic resistance. Here we review Staphylococcus outer surface features and specific phage resistance mechanisms that define the host range, the set of strains that an individual phage can potentially infect. Phage infection goes through five distinct phases: attachment, uptake, biosynthesis, assembly, and lysis. Adsorption inhibition, encompassing outer surface teichoic acid receptor alteration, elimination, or occlusion, limits successful phage attachment and entry. Restriction-modification systems (in particular, type I and IV systems), which target phage DNA inside the cell, serve as the major barriers to biosynthesis as well as transduction and horizontal gene transfer between clonal complexes and species. Resistance to late stages of infection occurs through mechanisms such as assembly interference, in which staphylococcal pathogenicity islands siphon away superinfecting phage proteins to package their own DNA. While genes responsible for teichoic acid biosynthesis, capsule, and restriction-modification are found in most Staphylococcus strains, a variety of other host range determinants (e.g., clustered regularly interspaced short palindromic repeats, abortive infection, and superinfection immunity) are sporadic. The fitness costs of phage resistance through teichoic acid structure alteration could make staphylococcal phage therapies promising, but host range prediction is complex because of the large number of genes involved, and the roles of many of these are unknown. In addition, little is known about the genetic determinants that contribute to host range expansion in the phages themselves. Future research must identify host range determinants, characterize resistance development during infection and treatment, and examine population-wide genetic background effects on resistance selection.

Whole-Genome Sequences of Staphylococcus aureus Isolates from Cystic Fibrosis Lung Infections

Staphylococcus aureus is an early colonizer in the lungs of individuals with cystic fibrosis (CF), but surprisingly, only a limited number of genomes from CF-associated S. aureus isolates have been sequenced. Here, we present the whole-genome sequences of 65 S. aureus isolates obtained from 50 individuals with CF.

USA300 Staphylococcus aureus persists on multiple body sites following an infection

BACKGROUND

USA300 methicillin-resistant Staphylococcus aureus (MRSA) is a community- and hospital-acquired pathogen that frequently causes infections but also can survive on the human body asymptomatically as a part of the normal microbiota. We devised a comparative genomic strategy to track colonizing USA300 at different body sites after an initial infection. We sampled ST8 S. aureus from subjects at the site of a first known MRSA infection. Within 60 days of this infection and again 12 months later, each subject was tested for asymptomatic colonization in the nose, throat and perirectal region. 93 S. aureus strains underwent whole genome shotgun sequencing.

RESULTS

Among 28 subjects at the initial sampling time, we isolated S. aureus from the nose, throat and perirectal sites from 15, 11 and 15 of them, respectively. Twelve months later we isolated S. aureus from 9 subjects, with 6, 3 and 3 strains from the nose, throat and perirectal area, respectively. Genome sequencing revealed that 23 patients (ages 0-66 years) carried USA300 intra-subject lineages (ISLs), defined as having an index infection isolate and closely related colonizing strains. Pairwise distance between strains in different ISLs was 48 to 162 single nucleotide polymorphisms (SNPs) across the core regions of the chromosome, whereas within the same ISL it was 0 to 26 SNPs. Strains in ISLs from the same subject differed in plasmid and prophage content, and contained deletions that removed the mecA-containing SCCmec and ACME regions. Five strains contained frameshift mutations in agr toxin-regulating genes. Persistence of an ISL was not associated with clinical or demographic subject characteristics. We inferred that colonization with the ISL occurred about 18 weeks before the first assessment of asymptomatic colonization.

CONCLUSIONS

Clonal lineages of USA300 may continue to colonize people at one or more anatomic sites up to a year after an initial infection and experience loss of the SCCmec, loss and gain of other mobile genetic elements, and mutations in the agr operon.

Whole genome sequencing for investigations of meningococcal outbreaks in the United States: a retrospective analysis

Although rare in the U.S., outbreaks due to Neisseria meningitidis do occur. Rapid, early outbreak detection is important for timely public health response. In this study, we characterized U.S. meningococcal isolates (N = 201) from 15 epidemiologically defined outbreaks (2009-2015) along with temporally and geographically matched sporadic isolates using multilocus sequence typing, pulsed-field gel electrophoresis (PFGE), and six whole genome sequencing (WGS) based methods. Recombination-corrected maximum likelihood (ML) and Bayesian phylogenies were reconstructed to identify genetically related outbreak isolates. All WGS analysis methods showed high degree of agreement and distinguished isolates with similar or indistinguishable PFGE patterns, or the same strain genotype. Ten outbreaks were caused by a single strain; 5 were due to multiple strains. Five sporadic isolates were phylogenetically related to 2 outbreaks. Analysis of 9 outbreaks using timed phylogenies identified the possible origin and estimated the approximate time that the most recent common ancestor emerged for outbreaks analyzed. U.S. meningococcal outbreaks were caused by single- or multiple-strain introduction, with organizational outbreaks mainly caused by a clonal strain and community outbreaks by divergent strains. WGS can infer linkage of meningococcal cases when epidemiological links are uncertain. Accurate identification of outbreak-associated cases requires both WGS typing and epidemiological data.

Fine-scale differentiation between Bacillus anthracis and Bacillus cereus group signatures in metagenome shotgun data

BACKGROUND

It is possible to detect bacterial species in shotgun metagenome datasets through the presence of only a few sequence reads. However, false positive results can arise, as was the case in the initial findings of a recent New York City subway metagenome project. False positives are especially likely when two closely related are present in the same sample. Bacillus anthracis, the etiologic agent of anthrax, is a high-consequence pathogen that shares >99% average nucleotide identity with Bacillus cereus group (BCerG) genomes. Our goal was to create an analysis tool that used k-mers to detect B. anthracis, incorporating information about the coverage of BCerG in the metagenome sample.

METHODS

Using public complete genome sequence datasets, we identified a set of 31-mer signatures that differentiated B. anthracis from other members of the B. cereus group (BCerG), and another set which differentiated BCerG genomes (including B. anthracis) from other Bacillus strains. We also created a set of 31-mers for detecting the lethal factor gene, the key genetic diagnostic of the presence of anthrax-causing bacteria. We created synthetic sequence datasets based on existing genomes to test the accuracy of a k-mer based detection model.

RESULTS

We found 239,503 B. anthracis-specific 31-mers (the Ba31 set), 10,183 BCerG 31-mers (the BCerG31 set), and 2,617 lethal factor k-mers (the lef31 set). We showed that false positive B. anthracis k-mers-which arise from random sequencing errors-are observable at high genome coverages of B. cereus. We also showed that there is a "gray zone" below 0.184× coverage of the B. anthracis genome sequence, in which we cannot expect with high probability to identify lethal factor k-mers. We created a linear regression model to differentiate the presence of B. anthracis-like chromosomes from sequencing errors given the BCerG background coverage. We showed that while shotgun datasets from the New York City subway metagenome project had no matches to lef31 k-mers and hence were negative for B. anthracis, some samples showed evidence of strains very closely related to the pathogen.

DISCUSSION

This work shows how extensive libraries of complete genomes can be used to create organism-specific signatures to help interpret metagenomes. We contrast "specialist" approaches to metagenome analysis such as this work to "generalist" software that seeks to classify all organisms present in the sample and note the more general utility of a k-mer filter approach when taxonomic boundaries lack clarity or high levels of precision are required.