The porcine skin microbiome exhibits broad fungal antagonism.

The skin and its microbiome function to protect the host from pathogen colonization and environmental stressors. In this study, using the Wisconsin Miniature Swine™ model, we characterize the porcine skin fungal and bacterial microbiomes, identify bacterial isolates displaying antifungal activity, and use whole-genome sequencing to identify biosynthetic gene clusters encoding for secondary metabolites that may be responsible for the antagonistic effects on fungi. Through this comprehensive approach of paired microbiome sequencing with culturomics, we report the discovery of novel species of Corynebacterium and Rothia. Further, this study represents the first comprehensive evaluation of the porcine skin mycobiome and the evaluation of bacterial-fungal interactions on this surface. Several diverse bacterial isolates exhibit potent antifungal properties against opportunistic fungal pathogens in vitro. Genomic analysis of inhibitory species revealed a diverse repertoire of uncharacterized biosynthetic gene clusters suggesting a reservoir of novel chemical and biological diversity. Collectively, the porcine skin microbiome represents a potential unique source of novel antifungals.

Combined metagenomic- and culture-based approaches to investigate bacterial strain-level associations with medication-controlled mild-moderate atopic dermatitis.

Background: The skin microbiome is disrupted in atopic dermatitis (AD). Existing research focuses on moderate to severe, unmedicated disease. Objective: We sought to investigate metagenomic- and culture-based bacterial strain-level differences in mild, medicated AD and the effects these have on human keratinocytes (HKs). Methods: Skin swabs from anterior forearms were collected from 20 pediatric participants (11 participants with AD sampled at lesional and nonlesional sites and 9 age- and sex-matched controls). Participants had primarily mild to moderate AD and maintained medication use. Samples were processed for microbial metagenomic sequencing and bacterial isolation. Isolates identified as Staphylococcus aureus were tested for enterotoxin production. HK cultures were treated with cell-free conditioned media from representative Staphylococcus species to measure barrier effects. Results: Metagenomic sequencing identified significant differences in microbiome composition between AD and control groups. Differences were seen at the species and strain levels for Staphylococci, with S aureus found only in participants with AD and differences in Staphylococcus epidermidis strains between control and AD swabs. These strains showed differences in toxin gene presence, which was confirmed in vitro for S aureus enterotoxins. The strain from the participant with the most severe AD produced enterotoxin B levels more than 100-fold higher than the other strains (P < .001). Strains also displayed differential effects on HK metabolism and barrier function. Conclusions: Strain-level differences in toxin genes from Staphylococcus strains may explain varying effects on HK, with S aureus and non-aureus strains negatively affecting viability and barrier function. These differences are likely important in AD pathogenesis.

Global diversity of enterococci and description of 18 previously unknown species.

Enterococci are gut microbes of most land animals. Likely appearing first in the guts of arthropods as they moved onto land, they diversified over hundreds of millions of years adapting to evolving hosts and host diets. Over 60 enterococcal species are now known. Two species, Enterococcus faecalis and Enterococcus faecium, are common constituents of the human microbiome. They are also now leading causes of multidrug-resistant hospital-associated infection. The basis for host association of enterococcal species is unknown. To begin identifying traits that drive host association, we collected 886 enterococcal strains from widely diverse hosts, ecologies, and geographies. This identified 18 previously undescribed species expanding genus diversity by >25%. These species harbor diverse genes including toxins and systems for detoxification and resource acquisition. Enterococcus faecalis and E. faecium were isolated from diverse hosts highlighting their generalist properties. Most other species showed a more restricted distribution indicative of specialized host association. The expanded species diversity permitted the Enterococcus genus phylogeny to be viewed with unprecedented resolution, allowing features to be identified that distinguish its four deeply rooted clades, and the entry of genes associated with range expansion such as B-vitamin biosynthesis and flagellar motility to be mapped to the phylogeny. This work provides an unprecedentedly broad and deep view of the genus Enterococcus, including insights into its evolution, potential new threats to human health, and where substantial additional enterococcal diversity is likely to be found.

skDER: microbial genome dereplication approaches for comparative and metagenomic applications.

skDER (https://github.com/raufs/skDER) combines recent advances to efficiently estimate average nucleotide identity (ANI) between thousands of microbial genomes by skani1 with two low-memory methods for genomic dereplication. The first method implements a dynamic algorithm to determine a concise set of representative genomes. This approach is well-suited for selecting reference genomes to align metagenomic reads onto for tracking strain presence across related microbiome samples. This is because fewer representative genomes should alleviate the concern that reads belonging to the same strain get falsely partitioned across closely related genomes. The other method, which uses a greedy approach, is better suited for use in comparative genomics, where users might be overwhelmed with the high number of genomes available for certain taxa and aim to reduce redundancy and, therefore, computational requirements for downstream analytics. This method selects a larger number of representative genomes to comprehensively sample the pangenome space for the taxon of interest. To further aid usage for comparative genomics studies, skDER also features an option to automatically download genomes classified as a particular species or genus in the Genome Taxonomy Database2-4 and we provide precomputed representative genomes for commonly studied bacterial taxa5.

Genomic epidemiology and antibiotic susceptibility profiling of uropathogenic Escherichia coli among children in the United States.

Escherichia coli is the most common cause of urinary tract infections (UTIs) in children, and yet the underlying mechanisms of virulence and antibiotic resistance and the overall population structure of the species is poorly understood within this age group. To investigate whether uropathogenic E. coli (UPEC) from children who developed pyelonephritis carried specific genetic markers, we generated whole-genome sequence data for 96 isolates from children with UTIs. This included 57 isolates from children with either radiologically confirmed pyelonephritis or cystitis and 27 isolates belonging to the well-known multidrug-resistant sequence type ST131, selected to investigate their population structure and antibiotic resistance characteristics. We observed a UPEC population structure that is similar to those reported in adults. In comparison with prior investigations, we found that the full pap operon was more common among UPEC from pediatric cases of pyelonephritis. Further, in contrast with recent reports that the P-fimbriae adhesin-encoding papGII allele is substantially more prevalent in invasive UPEC from adults, we found papGII was common to both invasive and non-invasive UPEC from children. Among the set of ST131 isolates from children with UTIs, we found antibiotic resistance was correlated with known genetic markers of resistance, as in adults. Unexpectedly, we observed that fimH30, an allele of the fimbrial gene fimH often used as a proxy to type ST131 isolates into the most drug-resistant subclade C, was carried by some of the subclade A and subclade B isolates, suggesting that the fimH30 allele could confer a selective advantage for UPEC. IMPORTANCE Urinary tract infections (UTIs), which are most often caused by Escherichia coli, are not well studied in children. Here, we examine genetic characteristics that differentiate UTI-causing bacteria in children that either remain localized to the bladder or are involved in more serious kidney infections. We also examine patterns of antibiotic resistance among strains from children that are part of E. coli sequence type 131, a group of bacteria that commonly cause UTIs and are known to have high levels of drug resistance. This work provides new insight into the virulence and antibiotic resistance characteristics of the bacteria that cause UTIs in children.

zol & fai: large-scale targeted detection and evolutionary investigation of gene clusters.

Many universally and conditionally important genes are genomically aggregated within clusters. Here, we introduce fai and zol, which together enable large-scale comparative analysis of different types of gene clusters and mobile-genetic elements (MGEs), such as biosynthetic gene clusters (BGCs) or viruses. Fundamentally, they overcome a current bottleneck to reliably perform comprehensive orthology inference at large scale across broad taxonomic contexts and thousands of genomes. First, fai allows the identification of orthologous or homologous instances of a query gene cluster of interest amongst a database of target genomes. Subsequently, zol enables reliable, context-specific inference of protein-encoding ortholog groups for individual genes across gene cluster instances. In addition, zol performs functional annotation and computes a variety of statistics for each inferred ortholog group. These programs are showcased through application to: (i) longitudinal tracking of a virus in metagenomes, (ii) discovering novel population-genetic insights of two common BGCs in a fungal species, and (iii) uncovering large-scale evolutionary trends of a virulence-associated gene cluster across thousands of genomes from a diverse bacterial genus.

Global diversity of enterococci and description of 18 novel species.

Enterococci are commensal gut microbes of most land animals. They diversified over hundreds of millions of years adapting to evolving hosts and host diets. Of over 60 known enterococcal species, Enterococcus faecalis and E. faecium uniquely emerged in the antibiotic era among leading causes of multidrug resistant hospital-associated infection. The basis for the association of particular enterococcal species with a host is largely unknown. To begin deciphering enterococcal species traits that drive host association, and to assess the pool of Enterococcus-adapted genes from which known facile gene exchangers such as E. faecalis and E. faecium may draw, we collected 886 enterococcal strains from nearly 1,000 specimens representing widely diverse hosts, ecologies and geographies. This provided data on the global occurrence and host associations of known species, identifying 18 new species in the process expanding genus diversity by >25%. The novel species harbor diverse genes associated with toxins, detoxification, and resource acquisition. E. faecalis and E. faecium were isolated from a wide diversity of hosts highlighting their generalist properties, whereas most other species exhibited more restricted distributions indicative of specialized host associations. The expanded species diversity permitted the Enterococcus genus phylogeny to be viewed with unprecedented resolution, allowing features to be identified that distinguish its four deeply rooted clades as well as genes associated with range expansion, such as B-vitamin biosynthesis and flagellar motility. Collectively, this work provides an unprecedentedly broad and deep view of the genus Enterococcus, potential threats to human health, and new insights into its evolution.

Comparative Genomic and Metagenomic Investigations of the Corynebacterium tuberculostearicum Species Complex Reveals Potential Mechanisms Underlying Associations To Skin Health and Disease.

Corynebacterium are a diverse genus and dominant member of the human skin microbiome. Recently, we reported that the most prevalent Corynebacterium species found on skin, including Corynebacterium tuberculostearicum and Corynebacterium kefirresidentii, comprise a narrow species complex despite the diversity of the genus. Here, we apply high-resolution phylogenomics and comparative genomics to describe the structure of the C. tuberculostearicum species complex and highlight genetic traits which are enriched or depleted in it relative to other Corynebacterium. Through metagenomic investigations, we also find that individual species within the complex can associate with specific body sites. Finally, we discover that one species from the complex, C. kefirresidentii, increases in relative abundance during atopic dermatitis flares, and show that most genomes of this species encode a colocalized set of putative virulence genes. IMPORTANCE Corynebacterium are commonly found bacteria on the human skin. In this study, we perform comparative genomics to gain insight into genetic traits which differentiate a phylogenetically related group of Corynebacterium, the Corynebacterium tuberculostearicum species complex, that includes the most prevalent species from the genus in skin microbiomes. After resolving the presence of distinct species within the complex, we applied metagenomic analysis to uncover biogeographic associations of individual species within the complex with specific body sites and discovered that one species, commonly found in the nares of individuals, increases in abundance across multiple body sites during atopic dermatitis flares.

vRhyme enables binning of viral genomes from metagenomes.

Genome binning has been essential for characterization of bacteria, archaea, and even eukaryotes from metagenomes. Yet, few approaches exist for viruses. We developed vRhyme, a fast and precise software for construction of viral metagenome-assembled genomes (vMAGs). vRhyme utilizes single- or multi-sample coverage effect size comparisons between scaffolds and employs supervised machine learning to identify nucleotide feature similarities, which are compiled into iterations of weighted networks and refined bins. To refine bins, vRhyme utilizes unique features of viral genomes, namely a protein redundancy scoring mechanism based on the observation that viruses seldom encode redundant genes. Using simulated viromes, we displayed superior performance of vRhyme compared to available binning tools in constructing more complete and uncontaminated vMAGs. When applied to 10,601 viral scaffolds from human skin, vRhyme advanced our understanding of resident viruses, highlighted by identification of a Herelleviridae vMAG comprised of 22 scaffolds, and another vMAG encoding a nitrate reductase metabolic gene, representing near-complete genomes post-binning. vRhyme will enable a convention of binning uncultivated viral genomes and has the potential to transform metagenome-based viral ecology.

Inter-species geographic signatures for tracing horizontal gene transfer and long-term persistence of carbapenem resistance.

BACKGROUND: Carbapenem-resistant Enterobacterales (CRE) are an urgent global health threat. Inferring the dynamics of local CRE dissemination is currently limited by our inability to confidently trace the spread of resistance determinants to unrelated bacterial hosts. Whole-genome sequence comparison is useful for identifying CRE clonal transmission and outbreaks, but high-frequency horizontal gene transfer (HGT) of carbapenem resistance genes and subsequent genome rearrangement complicate tracing the local persistence and mobilization of these genes across organisms. METHODS: To overcome this limitation, we developed a new approach to identify recent HGT of large, near-identical plasmid segments across species boundaries, which also allowed us to overcome technical challenges with genome assembly. We applied this to complete and near-complete genome assemblies to examine the local spread of CRE in a systematic, prospective collection of all CRE, as well as time- and species-matched carbapenem-susceptible Enterobacterales, isolated from patients from four US hospitals over nearly 5 years. RESULTS: Our CRE collection comprised a diverse range of species, lineages, and carbapenem resistance mechanisms, many of which were encoded on a variety of promiscuous plasmid types. We found and quantified rearrangement, persistence, and repeated transfer of plasmid segments, including those harboring carbapenemases, between organisms over multiple years. Some plasmid segments were found to be strongly associated with specific locales, thus representing geographic signatures that make it possible to trace recent and localized HGT events. Functional analysis of these signatures revealed genes commonly found in plasmids of nosocomial pathogens, such as functions required for plasmid retention and spread, as well survival against a variety of antibiotic and antiseptics common to the hospital environment. CONCLUSIONS: Collectively, the framework we developed provides a clearer, high-resolution picture of the epidemiology of antibiotic resistance importation, spread, and persistence in patients and healthcare networks.

Genes Contributing to the Unique Biology and Intrinsic Antibiotic Resistance of Enterococcus faecalis.

The enterococci, which are among the leading causes of multidrug-resistant (MDR) hospital infection, are notable for their environmental ruggedness, which extends to intrinsic antibiotic resistance. To identify genes that confer this unique property, we used Tn-seq to comprehensively explore the genome of MDR Enterococcus faecalis strain MMH594 for genes important for growth in nutrient-containing medium and with low-level antibiotic challenge. As expected, a large core of genes for DNA replication, expression, and central metabolism, shared with other bacteria, are intolerant to transposon disruption. However, genes were identified that are important to E. faecalis that are either absent from or unimportant for Staphylococcus aureus and Streptococcus pneumoniae fitness when similarly tested. Further, 217 genes were identified that when challenged by sub-MIC antibiotic levels exhibited reduced tolerance to transposon disruption, including those previously shown to contribute to intrinsic resistance, and others not previously ascribed this role. E. faecalis is one of the few Gram-positive bacteria experimentally shown to possess a functional Entner-Doudoroff pathway for carbon metabolism, a pathway that contributes to stress tolerance in other microbes. Through functional genomics and network analysis we defined the unusual structure of this pathway in E. faecalis and assessed its importance. These approaches also identified toxin-antitoxin and related systems that are unique and active in E. faecalis Finally, we identified genes that are absent in the closest nonenterococcal relatives, the vagococci, and that contribute importantly to fitness with and without antibiotic selection, advancing an understanding of the unique biology of enterococci.IMPORTANCE Enterococci are leading causes of antibiotic-resistant infection transmitted in hospitals. The intrinsic hardiness of these organisms allows them to survive disinfection practices and then proliferate in the gastrointestinal tracts of antibiotic-treated patients. The objective of this study was to identify the underlying genetic basis for its unusual hardiness. Using a functional genomic approach, we identified traits and pathways of general importance for enterococcal survival and growth that distinguish them from closely related pathogens as well as ancestrally related species. We further identified unique traits that enable them to survive antibiotic challenge, revealing a large set of genes that contribute to intrinsic antibiotic resistance and a smaller set of uniquely important genes that are rare outside enterococci.

Key Transitions in the Evolution of Rapid and Slow Growing Mycobacteria Identified by Comparative Genomics.

Mycobacteria have been classified into rapid and slow growing phenotypes, but the genetic factors that underlie these growth rate differences are not well understood. We compared the genomes of 157 mycobacterial species, representing all major branches of the mycobacterial phylogenetic tree to identify genes and operons enriched among rapid and slow growing mycobacteria. Overlaying growth phenotype on a phylogenetic tree based on 304 core genes suggested that ancestral mycobacteria had a rapid growth phenotype with a single major evolutionary separation into rapid and slow growing sub-genera. We identified 293 genes enriched among rapid growing sub-genera, including genes encoding for amino acid transport/metabolism (e.g., livFGMH operon) and transcription, as well as novel ABC transporters. Loss of the livFGMH and ABC transporter operons among slow growing species suggests that reduced cellular amino acid transport may be growth limiting. Comparative genomic analysis suggests that horizontal gene transfer, from non-mycobacterial genera, may have contributed to niche adaptation and pathogenicity, especially among slow growing species. Interestingly, the mammalian cell entry (mce) operon was found to be ubiquitous, irrespective of growth phenotype or pathogenicity, although protein sequence homology between rapid and slow growing species was low (<50%). This suggests that the mce operon was present in ancestral rapid growing species, but later adapted by slow growing species for use as a mechanism to establish an intra-cellular lifestyle.

Strain-specific quantification of root colonization by plant growth promoting rhizobacteria Bacillus firmus I-1582 and Bacillus amyloliquefaciens QST713 in non-sterile soil and field conditions.

Bacillus amyloliquefaciens QST713 and B. firmus I-1582 are bacterial strains which are used as active ingredients of commercially-available soil application and seed treatment products Serenade® and VOTiVO®, respectively. These bacteria colonize plant roots promoting plant growth and offering protection against pathogens/pests. The objective of this study was to develop a qPCR protocol to quantitate the dynamics of root colonization by these two strains under field conditions. Primers and TaqMan® probes were designed based on genome comparisons of the two strains with publicly-available and unpublished bacterial genomes of the same species. An optimized qPCR protocol was developed to quantify bacterial colonization of corn roots after seed treatment. Treated corn seeds were planted in non-sterile soil in the greenhouse and grown for 28 days. Specific detection of bacteria was quantified weekly, and showed stable colonization between ~104-105 CFU/g during the experimental period for both bacteria, and the protocol detected as low as 103 CFU/g bacteria on roots. In a separate experiment, streptomycin-resistant QST713 and rifampicin-resistant I-1582 strains were used to compare dilution-plating on TSA with the newly developed qPCR method. Results also indicated that the presence of natural microflora and another inoculated strain does not affect root colonization of either one of these strains. The same qPCR protocol was used to quantitate root colonization by QST713 and I-1582 in two corn and two soybean varieties grown in the field. Both bacteria were quantitated up to two weeks after seeds were planted in the field and there were no significant differences in root colonization in either bacteria strain among varieties. Results presented here confirm that the developed qPCR protocol can be successfully used to understand dynamics of root colonization by these bacteria in plants growing in growth chamber, greenhouse and the field.

Genomic evidence of geographically widespread effect of gene flow from polar bears into brown bears.

Polar bears are an arctic, marine adapted species that is closely related to brown bears. Genome analyses have shown that polar bears are distinct and genetically homogeneous in comparison to brown bears. However, these analyses have also revealed a remarkable episode of polar bear gene flow into the population of brown bears that colonized the Admiralty, Baranof and Chichagof islands (ABC islands) of Alaska. Here, we present an analysis of data from a large panel of polar bear and brown bear genomes that includes brown bears from the ABC islands, the Alaskan mainland and Europe. Our results provide clear evidence that gene flow between the two species had a geographically wide impact, with polar bear DNA found within the genomes of brown bears living both on the ABC islands and in the Alaskan mainland. Intriguingly, while brown bear genomes contain up to 8.8% polar bear ancestry, polar bear genomes appear to be devoid of brown bear ancestry, suggesting the presence of a barrier to gene flow in that direction.

Genomic evidence for island population conversion resolves conflicting theories of polar bear evolution.

Despite extensive genetic analysis, the evolutionary relationship between polar bears (Ursus maritimus) and brown bears (U. arctos) remains unclear. The two most recent comprehensive reports indicate a recent divergence with little subsequent admixture or a much more ancient divergence followed by extensive admixture. At the center of this controversy are the Alaskan ABC Islands brown bears that show evidence of shared ancestry with polar bears. We present an analysis of genome-wide sequence data for seven polar bears, one ABC Islands brown bear, one mainland Alaskan brown bear, and a black bear (U. americanus), plus recently published datasets from other bears. Surprisingly, we find clear evidence for gene flow from polar bears into ABC Islands brown bears but no evidence of gene flow from brown bears into polar bears. Importantly, while polar bears contributed <1% of the autosomal genome of the ABC Islands brown bear, they contributed 6.5% of the X chromosome. The magnitude of sex-biased polar bear ancestry and the clear direction of gene flow suggest a model wherein the enigmatic ABC Island brown bears are the descendants of a polar bear population that was gradually converted into brown bears via male-dominated brown bear admixture. We present a model that reconciles heretofore conflicting genetic observations. We posit that the enigmatic ABC Islands brown bears derive from a population of polar bears likely stranded by the receding ice at the end of the last glacial period. Since then, male brown bear migration onto the island has gradually converted these bears into an admixed population whose phenotype and genotype are principally brown bear, except at mtDNA and X-linked loci. This process of genome erosion and conversion may be a common outcome when climate change or other forces cause a population to become isolated and then overrun by species with which it can hybridize.