Search

Widespread antibiotic usage has resulted in the rapid evolution of drug-resistant bacterial pathogens and poses significant threats to public health. Resolving how pathogens respond to antibiotics under different contexts is critical for understanding disease emergence and evolution going forward. The impact of antibiotics has been demonstrated most directly through in vitro pathogen passaging experiments. Independent from antibiotic selection, interactions with hosts have also altered the evolutionary trajectories and fitness landscapes of pathogens, shaping infectious disease outcomes. However, it is unclear how interactions between hosts and antibiotics impact the evolution of pathogen virulence. Here, we evolved and re-sequenced Staphylococcus aureus , a major bacterial pathogen, varying exposure to host and antibiotics to tease apart the contributions of these selective pressures on pathogen adaptation. After 12 passages, S. aureus evolving in Caenorhabditis elegans nematodes exposed to a sub-minimum inhibitory concentration of antibiotic (oxacillin) became highly virulent, regardless of whether the ancestral pathogen was methicillin-resistant (MRSA) or methicillin-sensitive (MSSA). Host and antibiotic exposure selected for reduced drug susceptibility in MSSA lineages while increasing MRSA total growth outside hosts. We identified mutations in genes involved in complex regulatory networks linking virulence and metabolism, including codY , agr , and gdpP , suggesting that rapid adaptation to infect hosts may have pleiotropic effects. In particular, MSSA populations under selection from host and antibiotic accumulated mutations in the global regulator gene codY , which controls biofilm formation in S. aureus. These populations had indeed evolved more robust biofilms-a trait linked to both virulence and antibiotic resistance-suggesting evolution of one trait can confer multiple adaptive benefits. Despite evolving in similar environments, MRSA and MSSA populations proceeded on divergent evolutionary paths, with MSSA populations exhibiting more similarities across replicate populations. Our results underscore the importance of considering multiple and concurrent selective pressures as drivers of pervasive pathogen traits.

Background: Tracking the spread of antibiotic resistant bacteria is critical to reduce global morbidity and mortality associated with human and animal infections. There is a need to understand the role that wild animals in maintenance and transfer of antibiotic resistance genes (ARGs)., Methods: This study used metagenomics to identify and compare the abundance of bacterial species and ARGs detected in the gut microbiomes from sympatric humans and wild mouse lemurs in a forest-dominated, roadless region of Madagascar near Ranomafana National Park. We examined the contribution of human geographic location toward differences in ARG abundance and compared the genomic similarity of ARGs between host source microbiomes., Results: Alpha and beta diversity of species and ARGs between host sources were distinct but maintained a similar number of detectable ARG alleles. Humans were differentially more abundant for four distinct tetracycline resistance-associated genes compared to lemurs. There was no significant difference in human ARG diversity from different locations. Human and lemur microbiomes shared 14 distinct ARGs with highly conserved in nucleotide identity. Synteny of ARG-associated assemblies revealed a distinct multidrug-resistant gene cassette carrying dfrA1 and aadA1 present in human and lemur microbiomes without evidence of geographic overlap, suggesting that these resistance genes could be widespread in this ecosystem. Further investigation into intermediary processes that maintain drug-resistant bacteria in wildlife settings is needed., Competing Interests: Timothy D. Read and Patricia C. Wright are Academic Editors for PeerJ., (©2024 Talbot et al.)

Background: The obligate intracellular bacterial family Chlamydiaceae comprises a number of different species that cause disease in various vertebrate hosts including humans. Chlamydia suis, primarily found in the gastrointestinal tract of pigs, is the only species of the Chlamydiaceae family to have naturally gained tetracycline resistance (TetR), through a genomic island (Tet-island), integrated into the middle of chromosomal invasin-like gene inv. Previous studies have hypothesised that the uptake of the Tet-island from a host outside the Chlamydiaceae family was a unique event, followed by spread among C. suis through homologous recombination. In vitro recombination studies have shown that Tet-island exchange between C. suis strains is possible. Our aim in this study was to gain a deeper understanding of the interclade recombination of the Tet-island, among currently circulating C. suis field strains compared to in vitro-generated recombinants, using published whole genome sequences of C. suis field strains (n = 35) and in vitro-generated recombinants (n = 63)., Results: We found that the phylogeny of inv better reflected the phylogeny of the Tet-island than that of the whole genome, supporting recombination rather than site-specific insertion as the means of transfer. There were considerable differences between the distribution of recombinations within in vitro-generated strains compared to that within the field strains. These differences are likely because in vitro-generated recombinants were selected for a tetracycline and rifamycin/rifampicin resistant background, leading to the largest peak of recombination across the Tet-island. Finally, we found that interclade recombinations across the Tet-island were more variable in length downstream of the Tet-island than upstream., Conclusions: Our study supports the hypothesis that the occurrence of TetR strains in both clades of C. suis came about through interclade recombination after a single ancestral horizontal gene transfer event., (© 2024. The Author(s).)

The metabolic intimacy of symbiosis often demands the work of specialists. Natural products and defensive secondary metabolites can drive specificity by ensuring infection and propagation across host generations. But in contrast to bacteria, little is known about the diversity and distribution of natural product biosynthetic pathways among fungi and how they evolve to facilitate symbiosis and adaptation to their host environment. In this study, we define the secondary metabolism of Escovopsis and closely related genera, symbionts in the gardens of fungus-farming ants. We ask how the gain and loss of various biosynthetic pathways correspond to divergent lifestyles. Long-read sequencing allowed us to define the chromosomal features of representative Escovopsis strains, revealing highly reduced genomes composed of seven to eight chromosomes. The genomes are highly syntenic with macrosynteny decreasing with increasing phylogenetic distance, while maintaining a high degree of mesosynteny. An ancestral state reconstruction analysis of biosynthetic pathways revealed that, while many secondary metabolites are shared with non-ant-associated Sordariomycetes , 56 pathways are unique to the symbiotic genera. Reflecting adaptation to diverging ant agricultural systems, we observe that the stepwise acquisition of these pathways mirrors the ecological radiations of attine ants and the dynamic recruitment and replacement of their fungal cultivars. As different clades encode characteristic combinations of biosynthetic gene clusters, these delineating profiles provide important insights into the possible mechanisms underlying specificity between these symbionts and their fungal hosts. Collectively, our findings shed light on the evolutionary dynamic nature of secondary metabolism in Escovopsis and its allies, reflecting adaptation of the symbionts to an ancient agricultural system.IMPORTANCEMicrobial symbionts interact with their hosts and competitors through a remarkable array of secondary metabolites and natural products. Here, we highlight the highly streamlined genomic features of attine-associated fungal symbionts. The genomes of Escovopsis species, as well as species from other symbiont genera, many of which are common with the gardens of fungus-growing ants, are defined by seven chromosomes. Despite a high degree of metabolic conservation, we observe some variation in the symbionts' potential to produce secondary metabolites. As the phylogenetic distribution of the encoding biosynthetic gene clusters coincides with attine transitions in agricultural systems, we highlight the likely role of these metabolites in mediating adaptation by a group of highly specialized symbionts., Competing Interests: The authors declare no conflict of interest.

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 average nucleotide identity 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 genes were associated with prophages; strain-concentrated genes were associated with the vSaβ genome island and rare genes (<10% frequency) concentrated near the origin of replication. Antibiotic resistance 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., Importance: We analyzed the genomic diversity of Staphylococcus aureus , a globally prevalent bacterial species that causes serious infections in humans. Our goal was to build a genetic picture of the different strains of S. aureus and which genes may be associated with them. We reprocessed >84,000 genomes and subsampled to remove redundancy. We found that individual samples sharing >99.5% of their genome could be grouped into strains. We also showed that a portion of genes that are present in intermediate frequency in the species are strongly associated with some strains but completely absent from others, suggesting a role in strain specificity. This work lays the foundation for understanding individual gene histories of the S. aureus species and also outlines strategies for processing large bacterial genomic data sets., Competing Interests: The authors declare no conflict of interest.

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., Competing Interests: The authors declare no conflict of interest.

Staphylococcus aureus is a common cause of skin and soft-tissue infections (SSTIs) and has become the most common cause of bloodstream infections (BSIs) in recent years, but whether the strains causing these two clinical syndromes overlap has not been studied adequately. USA300/500 (clonal complex [CC] 8-sequence type [ST] 8) and USA100 (CC5-ST5) have dominated among methicillin-resistant S aureus (MRSA) strains in the United States since the early 2000s. We compared the genomes of unselected MRSA isolates from 131 SSTIs with those from 145 BSIs at a single US center in overlapping periods in 2018-2021. CC8 MRSA was more common among SSTIs, and CC5 was more common among BSIs, consistent with prior literature. Based on clustering genomes with a threshold of 15 single-nucleotide polymorphisms, we identified clusters limited to patients with SSTI and separate clusters exclusively comprising patients with BSIs. However, we also identified eight clusters that included at least one SSTI and one BSI isolate. This suggests that virulent MRSA strains are transmitted from person to person locally in the healthcare setting or the community and that single lineages are often capable of causing both SSTIs and BSIs., Competing Interests: Potential conflicts of interest. M. Z. D. reports support from GSK and Covance for clinical trials; service on a GSK advisor board; a contract from Johnson & Johnson to prepare a research report; and an honorarium for a lecture at the University of Iowa. All other authors report no potential conflicts., (© The Author(s) 2024. Published by Oxford University Press on behalf of Infectious Diseases Society of America.)

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., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

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.

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Bouiller, Jacko, Shumaker, Talbot, Read and David.)

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., Competing Interests: The authors declare no conflict of interest.

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., Competing Interests: Potential conflicts of interest. A.B. has served on a clinical advisory board for Beckman Coulter. All other authors: no reported conflicts., (© The Author(s) 2024. Published by Oxford University Press on behalf of Infectious Diseases Society of America.)

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., Competing Interests: The authors declare no conflict of interest.

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.

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 ( N eisseria c ell e nvelope r egulator). 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., Competing Interests: The authors declare no conflict of interest.

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., (© 2023. The Author(s).)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Joseph, Bommana, Ziklo, Kama, Dean and Read.)

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.

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Bianco, Moustafa, O’Brien, Martin, Read, Kreiswirth and Planet.)

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., Competing Interests: Potential conflicts of interest. D. A. P. is an associate editor for Clinical Infectious Diseases. M. Z. D. reports grants unrelated to this work from Johnson & Johnson and ContraFect (clinical trial); consulting fees from GSK; support from GSK for travel and lodging costs for a lecture given at GSK in Siena, Italy, December 2019; and participation on advisory board for GSK. T. D. R. reports payment or honoraria as an NIH peer reviewer; unpaid position on the editorial board for the Journal of Clinical Microbiology; and position as editor of PeerJ. All other authors report no potential conflicts. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed., (© The Author(s) 2022. Published by Oxford University Press on behalf of the Infectious Diseases Society of America. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Gonorrhea remains a major global public health problem because of the high incidence of infection (estimated 82 million cases in 2020) and the emergence and spread of Neisseria gonorrhoeae strains resistant to previous and current antibiotics used to treat infections. Given the dearth of new antibiotics that are likely to enter clinical practice in the near future, there is concern that cases of untreatable gonorrhea might emerge. In response to this crisis, the World Health Organization (WHO), in partnership with the Global Antibiotic Research and Development Partnership (GARDP), has made the search for and development of new antibiotics against N. gonorrhoeae a priority. Ideally, these antibiotics should also be active against other sexually transmitted organisms, such as Chlamydia trachomatis and/or Mycoplasma genitalium, which are often found with N. gonorrhoeae as co-infections. Corallopyronin A is a potent antimicrobial that exhibits activity against Chlamydia spp. and inhibits transcription by binding to the RpoB switch region. Accordingly, we tested the effectiveness of corallopyronin A against N. gonorrhoeae. We also examined the mutation frequency and modes of potential resistance against corallopyronin A. We report that corallopyronin A has potent antimicrobial action against antibiotic-susceptible and antibiotic-resistant N. gonorrhoeae strains and could eradicate gonococcal infection of cultured, primary human cervical epithelial cells. Critically, we found that spontaneous corallopyronin A-resistant mutants of N. gonorrhoeae are exceedingly rare (≤10 -10 ) when selected at 4× the MIC. Our results support pre-clinical studies aimed at developing corallopyronin A for gonorrheal treatment regimens. IMPORTANCE The high global incidence of gonorrhea, the lack of a protective vaccine, and the emergence of N. gonorrhoeae strains expressing resistance to currently used antibiotics demand that new treatment options be developed. Accordingly, we investigated whether corallopyronin A, an antibiotic which is effective against other pathogens, including C. trachomatis, which together with gonococci frequently cause co-infections in humans, could exert anti-gonococcal action in vitro and ex vivo , and potential resistance emergence. We propose that corallopyronin A be considered a potential future treatment option for gonorrhea because of its potent activity, low resistance development, and recent advances in scalable production.

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.

The continued emergence of Neisseria gonorrhoeae isolates which are resistant to first-line antibiotics has reinvigorated interest in alternative therapies such as expanded use of gentamicin (Gen). We hypothesized that expanded use of Gen promotes emergence of gonococci with clinical resistance to this aminoglycoside. To understand how decreased susceptibility of gonococci to Gen might develop, we selected spontaneous low-level Gen-resistant (Gen R ) mutants (Gen MIC = 32 μg/mL) of the Gen-susceptible strain FA19. Consequently, we identified a novel missense mutation in fusA , which encodes elongation factor G (EF-G), causing an alanine (A) to valine (V) substitution at amino acid position 563 in domain IV of EF-G; the mutant allele was termed fusA2 . Transformation analysis showed that fusA2 could increase the Gen MIC by 4-fold. While possession of fusA2 did not impair either in vitro gonococcal growth or protein synthesis, it did result in a fitness defect during experimental infection of the lower genital tract in female mice. Through bioinformatic analysis of whole-genome sequences of 10,634 international gonococcal clinical isolates, other fusA alleles were frequently detected, but genetic studies revealed that they could not decrease Gen susceptibility in a similar manner to fusA2 . In contrast to these diverse international fusA alleles, the fusA2- encoded A563V substitution was detected in only a single gonococcal clinical isolate. We hypothesize that the rare occurrence of fusA2 in N. gonorrhoeae clinical isolates is likely due to a fitness cost during infection, but compensatory mutations which alleviate this fitness cost could emerge and promote Gen R in global strains.

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., (© The Author(s) 2022. Published by Oxford University Press on behalf of Infectious Diseases Society of America.)

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.

Phage therapy has been proposed as a possible alternative treatment for infections caused by the ubiquitous bacterial pathogen Staphylococcus aureus. However, successful therapy requires understanding the genetic basis of host range-the subset of strains in a species that could be killed by a particular phage. We searched diverse sets of S. aureus public genome sequences against a database of genes suggested from prior studies to influence host range to look for patterns of variation across the species. We found that genes encoding biosynthesis of molecules that were targets of S. aureus phage adsorption to the outer surface of the cell were the most conserved in the pangenome. Putative phage resistance genes that were core components of the pangenome genes had similar nucleotide diversity, ratio of nonsynonymous to synonymous substitutions, and functionality (measured by delta-bitscore) to other core genes. However, phage resistance genes that were not part of the core genome were significantly less consistent with the core genome phylogeny than all noncore genes in this set, suggesting more frequent movement between strains by horizontal gene transfer. Only superinfection immunity genes encoded by temperate phages inserted in the genome correlated with experimentally determined temperate phage resistance. Taken together, these results suggested that, while phage adsorption genes are heavily conserved in the S. aureus species, HGT may play a significant role in strain-specific evolution of host range patterns. IMPORTANCE Staphylococcus aureus is a widespread, hospital- and community-acquired pathogen that is commonly antibiotic resistant. It causes diverse diseases affecting both the skin and internal organs. Its ubiquity, antibiotic resistance, and disease burden make new therapies urgent, such as phage therapy, in which viruses specific to infecting bacteria clear infection. S. aureus phage host range not only determines whether phage therapy will be successful by killing bacteria but also horizontal gene transfer through transduction of host genetic material by phages. In this work, we comprehensively reviewed existing literature to build a list of S. aureus phage resistance genes and searched our database of almost 43,000 S. aureus genomes for these genes to understand their patterns of evolution, finding that prophages' superinfection immunity correlates best with phage resistance and HGT. These findings improved our understanding of the relationship between known phage resistance genes and phage host range in the species.

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., (© The Author(s) 2022. Published by Oxford University Press.)

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 11 C-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.

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., Competing Interests: None declared.sss, (© 2021 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

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., Competing Interests: Timothy D. Read is an Academic Editor for PeerJ., (© 2021 Su et al.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Marti, Bommana, Read, Pesch, Prähauser, Dean and Borel.)

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.

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., Competing Interests: Timothy D. Read is an Academic Editor for PeerJ., (©2021 Rogers et al.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Dunlop, Satten, Hu, Knight, Hill, Wright, Smith, Read, Pearce and Corwin.)

Mutations within the mtrR gene are commonly found amongst multidrug resistant clinical isolates of Neisseria gonorrhoeae, which has been labelled a superbug by the Centers for Disease Control and Prevention. These mutations appear to contribute to antibiotic resistance by interfering with the ability of MtrR to bind to and repress expression of its target genes, which include the mtrCDE multidrug efflux transporter genes and the rpoH oxidative stress response sigma factor gene. However, the DNA-recognition mechanism of MtrR and the consensus sequence within these operators to which MtrR binds has remained unknown. In this work, we report the crystal structures of MtrR bound to the mtrCDE and rpoH operators, which reveal a conserved, but degenerate, DNA consensus binding site 5'-MCRTRCRN4YGYAYGK-3'. We complement our structural data with a comprehensive mutational analysis of key MtrR-DNA contacts to reveal their importance for MtrR-DNA binding both in vitro and in vivo. Furthermore, we model and generate common clinical mutations of MtrR to provide plausible biochemical explanations for the contribution of these mutations to multidrug resistance in N. gonorrhoeae. Collectively, our findings unveil key biological mechanisms underlying the global stress responses of N. gonorrhoeae., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Chlamydia trachomatis , an obligately intracellular bacterium, is the most prevalent cause of bacterial sexually transmitted infections (STIs) worldwide. Numbers of U.S. infections of the urogenital tract and rectum have increased annually. Because C. trachomatis is not easily cultured, comparative genomic studies are limited, restricting our understanding of strain diversity and emergence among populations globally. While Agilent SureSelect XT target enrichment RNA bait libraries have been developed for whole-genome enrichment and sequencing of C. trachomatis directly from clinical urine, vaginal, conjunctival, and rectal samples, public access to these libraries is not available. We therefore designed an RNA bait library (34,795 120-mer probes based on 85 genomes, versus 33,619 probes using 74 genomes in a previous one) to augment organism sequencing from clinical samples that can be shared with the scientific community, enabling comparison studies. We describe the library and limit of detection for genome copy input, and we present results of 100% efficiency and high-resolution determination of recombination and identical genomes within vaginal-rectal specimen pairs in women. This workflow provides a robust approach for discerning genomic diversity and advancing our understanding of the molecular epidemiology of contemporary C. trachomatis STIs across sample types, geographic populations, sexual networks, and outbreaks associated with proctitis/proctocolitis among women and men who have sex with men. IMPORTANCE Chlamydia trachomatis is an obligate intracellular bacterium that is not easily cultured, which limits our understanding of urogenital and rectal C. trachomatis transmission and impact on morbidity. To provide a publicly available workflow for whole-genome target enrichment and sequencing of C. trachomatis directly from clinical urine, vaginal, conjunctival, and rectal specimens, we developed and report on an RNA bait library to enrich the organism from clinical samples for sequencing. We demonstrate an increased efficiency in the percentage of reads mapping to C. trachomatis and identified recombinant and identical C. trachomatis genomes in paired vaginal-rectal samples from women. Our workflow provides a robust genomic epidemiologic approach to advance our understanding of C. trachomatis strains causing ocular, urogenital, and rectal infections and to explore geo-sexual networks, outbreaks of colorectal infections among women and men who have sex with men, and the role of these strains in morbidity.

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. IMPORTANCE Staphylococcus 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., (Copyright © 2021 Moller et al.)

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., (© 2020 The American Laryngological, Rhinological and Otological Society, Inc.)

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/bactopia IMPORTANCE 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., (Copyright © 2020 Petit and Read.)

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 IC 50 ) 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 IC 50 (r = -.11, P < .05). Moreover, there was no statistically significant correlation between the DEX IC 50 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., (© 2020 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

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. IMPORTANCE Staphylococcus 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., (Copyright © 2020 Bernardy et al.)

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., Competing Interests: Timothy D. Read is an Academic Editor for PeerJ., (©2020 Su et al.)

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.