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

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

Mycobacterium tuberculosis epidemiology in Oman: whole-genome sequencing uncovers transmission pathways.

Tuberculosis (TB) originating from expatriates that hail from high TB-burden countries is hypothesized to play a role in continued TB transmission in Oman. Here, we used whole-genome sequencing (WGS) to assess national TB transmission dynamics. The annual incidence per 100,000 population per year was calculated for nationals and expatriates. A convenience sample of Mycobacterium tuberculosis (MTB) isolates from 2018 to 2019 was sequenced and analyzed with publicly available TB sequences from Bangladesh, Tanzania, the Philippines, India, and Pakistan. Relatedness was assessed by generating core-genome single nucleotide polymorphism (SNP) distances. The incidence of TB was five cases per 100,000 persons in 2018 and seven cases per 100,000 persons in 2020 (R2 = 0.34, P = 0.60). Incidence among nationals was 3.9 per 100,000 persons in 2018 and 3.5 per 100,000 persons in 2020 (R2 = 0.20, P = 0.70), and incidence among expatriates was 7.2 per 100,000 persons in 2018 and 12.7 per 100,000 persons in 2020 (R2 = 0.74, P = 0.34). Sixty-eight local MTB isolates were sequenced and analyzed with 393 global isolates. Isolates belonged to nine distinct spoligotypes. Two isolates, originating from an expatriate and an Omani national, were grouped into a WGS-based cluster (SNP distance < 12), which was corroborated by an epidemiological investigation. Relatedness of local and global isolates (SNP distance < 100) was also seen. The relatedness between MTB strains in Oman and those in expatriate countries of origin can aid inform TB control policy. Our results provide evidence that WGS can complement epidemiological analysis to achieve the End TB strategy goal in Oman. IMPORTANCE Tuberculosis (TB) incidence in Oman remains above national program control targets. TB transmission originating from expatriates from high TB-burden countries has been hypothesized to play a role. We used whole-genome sequencing (WGS) to assess TB transmission dynamics between expatriates and Omani nationals to inform TB control efforts. Available Mycobacterium tuberculosis isolates from 2018 to 2019 underwent WGS and analysis with publicly available TB sequences from Bangladesh, the Philippines, India, and Pakistan to assess for genetic relatedness. Our analysis revealed evidence of previously unrecognized transmission between an expatriate and an Omani national, which was corroborated by epidemiological investigation. Analysis of local and global isolates revealed evidence of distant relatedness between local and global isolates. Our results provide evidence that WGS can complement classic public health surveillance to inform targeted interventions to achieve the End TB strategy goal in Oman.

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

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

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

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

Clinical and Genomic Epidemiology of mcr-9-Carrying Carbapenem-Resistant Enterobacterales Isolates in Metropolitan Atlanta, 2012 to 2017.

Colistin is a last-resort antibiotic for multidrug-resistant Gram-negative infections. Recently, the ninth allele of the mobile colistin resistance (mcr) gene family, designated mcr-9, was reported. However, its clinical and public health significance remains unclear. We queried genomes of carbapenem-resistant Enterobacterales (CRE) for mcr-9 from a convenience sample of clinical isolates collected between 2012 and 2017 through the Georgia Emerging Infections Program, a population- and laboratory-based surveillance program. Isolates underwent phenotypic characterization and whole-genome sequencing. Phenotypic characteristics, genomic features, and clinical outcomes of mcr-9-positive and -negative CRE cases were then compared. Among 235 sequenced CRE genomes, 13 (6%) were found to harbor mcr-9, all of which were Enterobacter cloacae complex. The median MIC and rates of heteroresistance and inducible resistance to colistin were similar between mcr-9-positive and -negative isolates. However, rates of resistance were higher among mcr-9-positive isolates across most antibiotic classes. All cases had significant health care exposures. The 90-day mortality was similarly high in both mcr-9-positive (31%) and -negative (7%) CRE cases. Nucleotide identity and phylogenetic analysis did not reveal geotemporal clustering. mcr-9-positive isolates had a significantly higher number of median [range] antimicrobial resistance (AMR) genes (16 [4 to 22] versus 6 [2 to 15]; P < 0.001) than did mcr-9-negative isolates. Pangenome tests confirmed a significant association of mcr-9 detection with mobile genetic element and heavy metal resistance genes. Overall, the presence of mcr-9 was not associated with significant changes in colistin resistance or clinical outcomes, but continued genomic surveillance to monitor for emergence of AMR genes is warranted. IMPORTANCE Colistin is a last-resort antibiotic for multidrug-resistant Gram-negative infections. A recently described allele of the mobile colistin resistance (mcr) gene family, designated mcr-9, has been widely reported among Enterobacterales species. However, its clinical and public health significance remains unclear. We compared characteristics and outcomes of mcr-9-positive and -negative CRE cases. All cases were acquired in the health care setting and associated with a high rate of mortality. The presence of mcr-9 was not associated with significant changes in colistin resistance, heteroresistance, or inducible resistance but was associated with resistance to other antimicrobials and antimicrobial resistance (AMR), virulence, and heavy metal resistance (HMR) genes. Overall, the presence of mcr-9 was not associated with significant phenotypic changes or clinical outcomes. However, given the increase in AMR and HMR gene content and potential clinical impact, continued genomic surveillance of multidrug-resistant organisms to monitor for emergence of AMR genes is warranted.

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

BACKGROUND: Though detection of transmission clusters of methicillin-resistant Staphylococcus aureus (MRSA) infections is a priority for infection control personnel in hospitals, the transmission dynamics of MRSA among hospitalized patients with bloodstream infections (BSIs) has not been thoroughly studied. Whole genome sequencing (WGS) of MRSA isolates for surveillance is valuable for detecting outbreaks in hospitals, but the bioinformatic approaches used are diverse and difficult to compare. METHODS: We combined short-read WGS with genotypic, phenotypic, and epidemiological characteristics of 106 MRSA BSI isolates collected for routine microbiological diagnosis from inpatients in 2 hospitals over 12 months. Clinical data and hospitalization history were abstracted from electronic medical records. We compared 3 genome sequence alignment strategies to assess similarity in cluster ascertainment. We conducted logistic regression to measure the probability of predicting prior hospital overlap between clustered patient isolates by the genetic distance of their isolates. RESULTS: While the 3 alignment approaches detected similar results, they showed some variation. A gene family-based alignment pipeline was most consistent across MRSA clonal complexes. We identified 9 unique clusters of closely related BSI isolates. Most BSIs were healthcare associated and community onset. Our logistic model showed that with 13 single-nucleotide polymorphisms, the likelihood that any 2 patients in a cluster had overlapped in a hospital was 50%. CONCLUSIONS: Multiple clusters of closely related MRSA isolates can be identified using WGS among strains cultured from BSI in 2 hospitals. Genomic clustering of these infections suggests that transmission resulted from a mix of community spread and healthcare exposures long before BSI diagnosis.

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

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

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

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

Species-Scale Genomic Analysis of Staphylococcus aureus Genes Influencing Phage Host Range and Their Relationships to Virulence and Antibiotic Resistance Genes.

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.

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

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

Bactopia: a Flexible Pipeline for Complete Analysis of Bacterial Genomes.

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

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

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

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

Importance: Invasive nontypeable Haemophilus influenzae (NTHi) infection among adults is typically associated with bacteremic pneumonia. Nontypeable H influenzae is genetically diverse and clusters of infection are uncommon. Objective: To evaluate an increase in invasive NTHi infection from 2017-2018 among HIV-infected men who have sex with men in metropolitan Atlanta, Georgia. Design, Setting, and Participants: A population-based surveillance study with a cohort substudy and descriptive epidemiological analysis identified adults aged 18 years or older with invasive NTHi infection (isolation of NTHi from a normally sterile site) between January 1, 2008, and December 31, 2018 (final date of follow-up). Exposures: Time period, HIV status, and genetic relatedness (ie, cluster status) of available NTHi isolates. Main Outcomes and Measures: The primary outcome was incidence of invasive NTHi infection (from 2008-2016 and 2017-2018) among persons with HIV and compared with NTHi infection from 2008-2018 among those without HIV. The secondary outcomes were assessed among those aged 18 to 55 years with invasive NTHi infection and included epidemiological, clinical, and geographic comparisons by cluster status. Results: Among 553 adults with invasive NTHi infection (median age, 66 years [Q1-Q3, 48-78 years]; 52% male; and 38% black), 60 cases occurred among persons with HIV. Incidence of invasive NTHi infection from 2017-2018 among persons with HIV (41.7 cases per 100 000) was significantly greater than from 2008-2016 among those with HIV (9.6 per 100 000; P < .001) and from 2008-2018 among those without HIV (1.1 per 100 000; P < .001). Among adults aged 18 to 55 years with invasive NTHi infections from 2017-2018 (n = 179), persons with HIV (n = 31) were significantly more likely than those from 2008-2018 without HIV (n = 124) to be male (94% vs 49%, respectively; P < .001), black (100% vs 53%; P < .001), and have septic arthritis (35% vs 1%; P < .001). Persons with HIV who had invasive NTHi infection from 2017-2018 (n = 31) were more likely than persons with HIV who had invasive NTHi infection from 2008-2016 (n = 24) to have septic arthritis (35% vs 4%, respectively; P = .01). Pulsed-field gel electrophoresis of 174 of 179 NTHi isolates from 18- to 55-year-olds identified 2 genetically distinct clonal groups: cluster 1 (C1; n = 24) and cluster 2 (C2; n = 23). Whole-genome sequencing confirmed 2 clonal lineages of NTHi infection and revealed all C1 isolates (but none of the C2 isolates) carried IS1016 (an insertion sequence associated with H influenzae capsule genes). Persons with HIV were significantly more likely to have C1 or C2 invasive NTHi infection from 2017-2018 (28/31 [90%]) compared with from 2008-2016 among persons with HIV (10/24 [42%]; P < .001) and compared with from 2008-2018 among those without HIV (9/119 [8%]; P < .001). Among persons with C1 or C2 invasive NTHi infection who had HIV (n = 38) (median age, 34.5 years; 100% male; 100% black; 82% men who have sex with men), 32 (84%) lived in 2 urban counties and an area of significant spatial aggregation was identified compared with those without C1 or C2 invasive NTHi infection. Conclusions and Relevance: Among persons with HIV in Atlanta, the incidence of invasive nontypeable H influenzae infection increased significantly from 2017-2018 compared with 2008-2016. Two unique but genetically related clonal strains were identified and were associated with septic arthritis among black men who have sex with men and who lived in geographic proximity.

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

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

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

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

USA300 Staphylococcus aureus persists on multiple body sites following an infection.

BACKGROUND: USA300 methicillin-resistant Staphylococcus aureus (MRSA) is a community- and hospital-acquired pathogen that frequently causes infections but also can survive on the human body asymptomatically as a part of the normal microbiota. We devised a comparative genomic strategy to track colonizing USA300 at different body sites after an initial infection. We sampled ST8 S. aureus from subjects at the site of a first known MRSA infection. Within 60 days of this infection and again 12 months later, each subject was tested for asymptomatic colonization in the nose, throat and perirectal region. 93 S. aureus strains underwent whole genome shotgun sequencing. RESULTS: Among 28 subjects at the initial sampling time, we isolated S. aureus from the nose, throat and perirectal sites from 15, 11 and 15 of them, respectively. Twelve months later we isolated S. aureus from 9 subjects, with 6, 3 and 3 strains from the nose, throat and perirectal area, respectively. Genome sequencing revealed that 23 patients (ages 0-66 years) carried USA300 intra-subject lineages (ISLs), defined as having an index infection isolate and closely related colonizing strains. Pairwise distance between strains in different ISLs was 48 to 162 single nucleotide polymorphisms (SNPs) across the core regions of the chromosome, whereas within the same ISL it was 0 to 26 SNPs. Strains in ISLs from the same subject differed in plasmid and prophage content, and contained deletions that removed the mecA-containing SCCmec and ACME regions. Five strains contained frameshift mutations in agr toxin-regulating genes. Persistence of an ISL was not associated with clinical or demographic subject characteristics. We inferred that colonization with the ISL occurred about 18 weeks before the first assessment of asymptomatic colonization. CONCLUSIONS: Clonal lineages of USA300 may continue to colonize people at one or more anatomic sites up to a year after an initial infection and experience loss of the SCCmec, loss and gain of other mobile genetic elements, and mutations in the agr operon.

Staphylococcus aureus viewed from the perspective of 40,000+ genomes.

Low-cost Illumina sequencing of clinically-important bacterial pathogens has generated thousands of publicly available genomic datasets. Analyzing these genomes and extracting relevant information for each pathogen and the associated clinical phenotypes requires not only resources and bioinformatic skills but organism-specific knowledge. In light of these issues, we created Staphopia, an analysis pipeline, database and application programming interface, focused on Staphylococcus aureus, a common colonizer of humans and a major antibiotic-resistant pathogen responsible for a wide spectrum of hospital and community-associated infections. Written in Python, Staphopia's analysis pipeline consists of submodules running open-source tools. It accepts raw FASTQ reads as an input, which undergo quality control filtration, error correction and reduction to a maximum of approximately 100× chromosome coverage. This reduction significantly reduces total runtime without detrimentally affecting the results. The pipeline performs de novo assembly-based and mapping-based analysis. Automated gene calling and annotation is performed on the assembled contigs. Read-mapping is used to call variants (single nucleotide polymorphisms and insertion/deletions) against a reference S. aureus chromosome (N315, ST5). We ran the analysis pipeline on more than 43,000 S. aureus shotgun Illumina genome projects in the public European Nucleotide Archive database in November 2017. We found that only a quarter of known multi-locus sequence types (STs) were represented but the top 10 STs made up 70% of all genomes. methicillin-resistant S. aureus (MRSA) were 64% of all genomes. Using the Staphopia database we selected 380 high quality genomes deposited with good metadata, each from a different multi-locus ST, as a non-redundant diversity set for studying S. aureus evolution. In addition to answering basic science questions, Staphopia could serve as a potential platform for rapid clinical diagnostics of S. aureus isolates in the future. The system could also be adapted as a template for other organism-specific databases.

Draft sequencing and assembly of the genome of the world's largest fish, the whale shark: Rhincodon typus Smith 1828.

BACKGROUND: The whale shark (Rhincodon typus) has by far the largest body size of any elasmobranch (shark or ray) species. Therefore, it is also the largest extant species of the paraphyletic assemblage commonly referred to as fishes. As both a phenotypic extreme and a member of the group Chondrichthyes - the sister group to the remaining gnathostomes, which includes all tetrapods and therefore also humans - its genome is of substantial comparative interest. Whale sharks are also listed as an endangered species on the International Union for Conservation of Nature's Red List of threatened species and are of growing popularity as both a target of ecotourism and as a charismatic conservation ambassador for the pelagic ecosystem. A genome map for this species would aid in defining effective conservation units and understanding global population structure. RESULTS: We characterised the nuclear genome of the whale shark using next generation sequencing (454, Illumina) and de novo assembly and annotation methods, based on material collected from the Georgia Aquarium. The data set consisted of 878,654,233 reads, which yielded a draft assembly of 1,213,200 contigs and 997,976 scaffolds. The estimated genome size was 3.44Gb. As expected, the proteome of the whale shark was most closely related to the only other complete genome of a cartilaginous fish, the holocephalan elephant shark. The whale shark contained a novel Toll-like-receptor (TLR) protein with sequence similarity to both the TLR4 and TLR13 proteins of mammals and TLR21 of teleosts. The data are publicly available on GenBank, FigShare, and from the NCBI Short Read Archive under accession number SRP044374. CONCLUSIONS: This represents the first shotgun elasmobranch genome and will aid studies of molecular systematics, biogeography, genetic differentiation, and conservation genetics in this and other shark species, as well as providing comparative data for studies of evolutionary biology and immunology across the jawed vertebrate lineages.

Transmission and microevolution of USA300 MRSA in U.S. households: evidence from whole-genome sequencing.

UNLABELLED: Methicillin-resistant Staphylococcus aureus (MRSA) USA300 is a successful S. aureus clone in the United States and a common cause of skin and soft tissue infections (SSTIs). We performed whole-genome sequencing (WGS) of 146 USA300 MRSA isolates from SSTIs and colonization cultures obtained from an investigation conducted from 2008 to 2010 in Chicago and Los Angeles households that included an index case with an S. aureus SSTI. Identifying unique single nucleotide polymorphisms (SNPs) and analyzing whole-genome phylogeny, we characterized isolates to understand transmission dynamics, genetic relatedness, and microevolution of USA300 MRSA within the households. We also compared the 146 USA300 MRSA isolates from our study with the previously published genome sequences of the USA300 MRSA isolates from San Diego (n = 35) and New York City (n = 277). We found little genetic variation within the USA300 MRSA household isolates from Los Angeles (mean number of SNPs ± standard deviation, 17.6 ± 35; π nucleotide diversity, 3.1 × 10(-5)) or from Chicago (mean number of SNPs ± standard deviation, 12 ± 19; π nucleotide diversity, 3.1 × 10(-5)). The isolates within a household clustered into closely related monophyletic groups, suggesting the introduction into and transmission within each household of a single common USA300 ancestral strain. From a Bayesian evolutionary reconstruction, we inferred that USA300 persisted within households for 2.33 to 8.35 years prior to sampling. We also noted that fluoroquinolone-resistant USA300 clones emerged around 1995 and were more widespread in Los Angeles and New York City than in Chicago. Our findings strongly suggest that unique USA300 MRSA isolates are transmitted within households that contain an individual with an SSTI. Decolonization of household members may be a critical component of prevention programs to control USA300 MRSA spread in the United States. IMPORTANCE: USA300, a virulent and easily transmissible strain of methicillin-resistant Staphylococcus aureus (MRSA), is the predominant community-associated MRSA clone in the United States. It most commonly causes skin infections but also causes necrotizing pneumonia and endocarditis. Strategies to limit the spread of MRSA in the community can only be effective if we understand the most common sources of transmission and the microevolutionary processes that provide a fitness advantage to MRSA. We performed a whole-genome sequence comparison of 146 USA300 MRSA isolates from Chicago and Los Angeles. We show that households represent a frequent site of transmission and a long-term reservoir of USA300 strains; individuals within households transmit the same USA300 strain among themselves. Our study also reveals that a large proportion of the USA300 isolates sequenced are resistant to fluoroquinolone antibiotics. The significance of this study is that if households serve as long-term reservoirs of USA300, household MRSA eradication programs may result in a uniquely effective control method.