Genome Reduction Is Associated with Bacterial Pathogenicity across Different Scales of Temporal and Ecological Divergence

Phenotypic and genetic heterogeneity of Acinetobacter baumannii in the course of an animal chronic infection

Acinetobacter baumannii is a nosocomial pathogen associated with various infections, including urinary tract infections (UTIs). In the course of an infection, A. baumannii is known to rapidly become resistant to antibiotic therapy, but much less is known about possible adaptation without antibiotic pressure. Through a retrospective study, we investigated within-host genetic diversity during a subclinical 5-year UTI in an animal-patient after withdrawal of colistin treatment. We conducted whole-genome sequencing and phenotypic assays on 17 clonally related isolates from the Sequence Type 25 lineage. Phylogenomic analysis revealed their proximity with animal and human strains from the same country suggesting zoonotic transmission (France). In this case study, the clonally related strains presented variations in genome sizes and nucleotide sequences. Over the course of the infection, A. baumannii underwent genome reduction through insertion sequence (IS) recombination, phage excision or plasmid curing. Alongside this global genome reduction, we observed an expansion of IS17, initially located on the endogenous large plasmid. Genetic variations were mainly located in biofilm formation and metabolism genes. We observed repeated variations affecting three biofilm genes and two adhesion operons associated with weak biofilm-forming capacity. Conversely, only two metabolic genes were recurrently affected, and phenotypic assays indicated a rather stable metabolism profile between the isolates suggesting minor adaptations to its host. Lastly, an overall decreased antibiotic resistance - expected in the absence of antibiotic treatment - contrasted with a conserved colistin resistance due to a pmrB mutation among the isolates.

A phylogenetic approach to comparative genomics

Comparative genomics, whereby the genomes of different species are compared, has the potential to address broad and fundamental questions at the intersection of genetics and evolution. However, species, genomes and genes cannot be considered as independent data points within statistical tests. Closely related species tend to be similar because they share genes by common descent, which must be accounted for in analyses. This problem of non-independence may be exacerbated when examining genomes or genes but can be addressed by applying phylogeny-based methods to comparative genomic analyses. Here, we review how controlling for phylogeny can change the conclusions of comparative genomics studies. We address common questions on how to apply these methods and illustrate how they can be used to test causal hypotheses. The combination of rapidly expanding genomic datasets and phylogenetic comparative methods is set to revolutionize the biological insights possible from comparative genomic studies.

Historical plague pandemics: perspectives from ancient DNA

Ancient DNA research has provided important insights into the evolutionary history of Yersinia pestis during the historical plague pandemics. Future work should prioritise a more diversified approach to sampling, to ensure a broader understanding of the factors underlying pandemic onset, spread, and impact across different regions and hosts.

Functional role of the biofilm regulator CsgD in Salmonella enterica sv. Typhi

Introduction

Typhoid fever is an infectious disease primarily caused by Salmonella enterica sv. Typhi (S. Typhi), a bacterium that causes as many as 20 million infections and 600,000 deaths annually. Asymptomatic chronic carriers of S. Typhi play a major role in the transmission of typhoid fever, as they intermittently shed the bacteria and can unknowingly infect humans in close proximity. An estimated 90% of chronic carriers have gallstones; biofilm formation on gallstones is a primary factor in the establishment and maintenance of gallbladder carriage. CsgD is a central biofilm regulator in Salmonella, but the S. Typhi csgD gene has a mutation that introduces an early stop codon, resulting in a protein truncated by 8 amino acids at the C-terminus. In this study, we investigate the role of role of CsgD in S. Typhi.

Methods

We introduced a fully functional copy of the csgD gene from S. Typhimurium into S. Typhi under both a native and a constitutive promoter and tested for red, dry, and rough (Rdar) colony morphology, curli fimbriae, cellulose, and biofilm formation.

Results and discussion

We demonstrate that although CsgD-regulated curli and cellulose production were partially restored, the introduction of the S. Typhimurium csgD did not induce the Rdar colony morphology. Interestingly, we show that CsgD does not have a significant role in S. Typhi biofilm formation, as biofilm-forming capacities depend more on the isolate than the CsgD regulator. This data suggests the presence of an alternative biofilm regulatory process in this human-restricted pathogen.

Comparative genomics and virulence potential of Campylobacter coli strains isolated from different sources over 25 years in Brazil

BACKGROUND

Campylobacter spp. have been reported as a common cause of gastroenteritis in humans in many countries. However, in Brazil there is insufficient data to estimate the impact of Campylobacter in public health. In light of the importance of this foodborne pathogen, the aim of this study was to perform comparative analyses on 80 Brazilian Campylobacter coli genomes isolated from human feces, animals, the environment, and food. Methods include Average Nucleotide Identity (ANI), Gegenees, genomic plasticity, presence of pathogenicity, resistance, and metabolic islands. In addition, virulence analysis in Galleria mellonella were also performed for 18 selected C. coli strains.

RESULTS

The ANI values confirmed that all strains belonged to the C. coli species. Phylogenetic analyses demonstrated the evolutionary relationships among the studied strains, highlighting the genetic diversity among them. The differences in shared and deleted regions of the studied genomes were demonstrated, with 16 genomic islands identified, including 4 metabolic islands, 4 resistance islands, and 8 pathogenicity islands. We detected genes associated with chemotaxis, exotoxin production, antimicrobial resistance, stress response, defense mechanisms, and intracellular survival among these islands, highlighting the pathogenic potential of these strains. Two strains isolated from human and one from animal showed high virulence, killing 100% of Galleria mellonella larvae. Two strains isolated from the environment and two isolated from food killed 70-90% of the larvae and were classified as virulent. Three strains isolated from animal, two from human, two from the environment and one from food killed 30% to 60% of the larvae and were considered of intermediate virulence. Campylobacter jejuni ATCC 33291, one strain isolated from human and one from food killed 10 to 20% of the larvae and were considered of low virulence. One strain isolated from food did not kill any larvae and was considered avirulent.

CONCLUSIONS

The results obtained highlighted the genetic diversity, pathogenic and virulence potential of many of the C. coli strains studied, contributing for a more complete characterization of this important pathogen recognized as a cause of human gastroenteritis.

Convergent reductive evolution in bee-associated lactic acid bacteria

Distantly related organisms may evolve similar traits when exposed to similar environments or engaging in certain lifestyles. Several members of the Lactobacillaceae [lactic acid bacteria (LAB)] family are frequently isolated from the floral niche, mostly from bees and flowers. In some floral LAB species (henceforth referred to as bee-associated LAB), distinctive genomic (e.g., genome reduction) and phenotypic (e.g., preference for fructose over glucose or fructophily) features were recently documented. These features are found across distantly related species, raising the hypothesis that specific genomic and phenotypic traits evolved convergently during adaptation to the floral environment. To test this hypothesis, we examined representative genomes of 369 species of bee-associated and non-bee-associated LAB. Phylogenomic analysis unveiled seven independent ecological shifts toward the bee environment in LAB. In these species, we observed significant reductions of genome size, gene repertoire, and GC content. Using machine leaning, we could distinguish bee-associated from non-bee-associated species with 94% accuracy, based on the absence of genes involved in metabolism, osmotic stress, or DNA repair. Moreover, we found that the most important genes for the machine learning classifier were seemingly lost, independently, in multiple bee-associated lineages. One of these genes, acetaldehyde-alcohol dehydrogenase (adhE), encodes a bifunctional aldehyde-alcohol dehydrogenase which has been associated with the evolution of fructophily, a rare phenotypic trait that is pervasive across bee-associated LAB species. These results suggest that the independent evolution of distinctive phenotypes in bee-associated LAB has been largely driven by independent losses of the same sets of genes.IMPORTANCESeveral LAB species are intimately associated with bees and exhibit unique biochemical properties with potential for food applications and honeybee health. Using a machine learning-based approach, our study shows that adaptation of LAB to the bee environment was accompanied by a distinctive genomic trajectory deeply shaped by gene loss. Several of these gene losses occurred independently in distantly related species and are linked to some of their unique biotechnologically relevant traits, such as the preference for fructose over glucose (fructophily). This study underscores the potential of machine learning in identifying fingerprints of adaptation and detecting instances of convergent evolution. Furthermore, it sheds light onto the genomic and phenotypic particularities of bee-associated bacteria, thereby deepening the understanding of their positive impact on honeybee health.

Genomic and functional divergence of Staphylococcus aureus strains from atopic dermatitis patients and healthy individuals: insights from global and local scales

Atopic dermatitis (AD) is the most common chronic inflammatory skin disease worldwide and is characterized by a complex interplay with skin microbiota, with Staphylococcus aureus often abnormally more abundant in AD patients than in healthy individuals (HE). S. aureus harbors diverse strains with varied genetic compositions and functionalities, which exhibit differential connections with the severity of AD. However, the differences in S. aureus strains between AD and HE remain unclear, with most variations seen at a specific geographic level, implying spontaneous adaptations rather than systematic distinctions. This study presents genomic and functional differences between these S. aureus strains from AD and HE on both global and local levels. We observed reduced gene content diversity but increased functional variation in the global AD-associated strains. Two additional AD-dominant clusters emerged, with Cluster 1 enriched in transposases and Cluster 2 showcasing genes linked to adaptability and antibiotic resistance. Particularly, robust evidence illustrates that the lantibiotic operon of S. aureus, involved in the biosynthesis of lantibiotics, was acquired via horizontal gene transfer from environmental bacteria. Comparisons of the gene abundance profiles in functional categories also indicate limited zoonotic potential between human and animal isolates. Local analysis mirrored global gene diversity but showed distinct functional variations between AD and HE strains. Overall, this research provides foundational insights into the genomic evolution, adaptability, and antibiotic resistance of S. aureus, with significant implications for clinical microbiology.IMPORTANCEOur study uncovers significant genomic variations in Staphylococcus aureus strains associated with atopic dermatitis. We observed adaptive evolution tailored to the disease microenvironment, characterized by a smaller pan-genome than strains from healthy skin both on the global and local levels. Key functional categories driving strain diversification include "replication and repair" and "transporters," with transposases being pivotal. Interestingly, the local strains predominantly featured metal-related genes, whereas global ones emphasized antimicrobial resistances, signifying scale-dependent diversification nuances. We also pinpointed horizontal gene transfer events, indicating interactions between human-associated and environmental bacteria. These insights expand our comprehension of S. aureus's genetic adaptation in atopic dermatitis, yielding valuable implications for clinical approaches.

Variation in bacterial pathotype is consistent with the sit-and-wait hypothesis

The sit-and-wait hypothesis predicts that bacteria can become more virulent when they survive and transmit outside of their hosts due to circumventing the costs of host mortality. While this hypothesis is largely supported theoretically and through comparative analysis, experimental validation is limited. Here we test this hypothesis in Streptococcus suis, an opportunistic zoonotic pig pathogen, where a pathogenic ecotype proliferated during the change to intensive pig farming that amplifies opportunities for fomite transmission. We show in an in vitro environmental survival experiment that pathogenic ecotypes survive for longer than commensal ecotypes, despite similar rates of decline. The presence of a polysaccharide capsule has no consistent effect on survival. Our findings suggest that extended survival in the food chain may augment the zoonotic capability of S. suis. Moreover, eliminating the long-term environmental survival of bacteria could be a strategy that will both enhance infection control and curtail the evolution of virulence.

Queuosine salvage in Bartonella henselae Houston 1: a unique evolutionary path

Queuosine (Q) stands out as the sole tRNA modification that can be synthesized via salvage pathways. Comparative genomic analyses identified specific bacteria that showed a discrepancy between the projected Q salvage route and the predicted substrate specificities of the two identified salvage proteins: (1) the distinctive enzyme tRNA guanine-34 transglycosylase (bacterial TGT, or bTGT), responsible for inserting precursor bases into target tRNAs; and (2) queuosine precursor transporter (QPTR), a transporter protein that imports Q precursors. Organisms such as the facultative intracellular pathogen Bartonella henselae, which possess only bTGT and QPTR but lack predicted enzymes for converting preQ1 to Q, would be expected to salvage the queuine (q) base, mirroring the scenario for the obligate intracellular pathogen Chlamydia trachomatis. However, sequence analyses indicate that the substrate-specificity residues of their bTGTs resemble those of enzymes inserting preQ1 rather than q. Intriguingly, MS analyses of tRNA modification profiles in B. henselae reveal trace amounts of preQ1, previously not observed in a natural context. Complementation analysis demonstrates that B. henselae bTGT and QPTR not only utilize preQ1, akin to their Escherichia coli counterparts, but can also process q when provided at elevated concentrations. The experimental and phylogenomic analyses suggest that the Q pathway in B. henselae could represent an evolutionary transition among intracellular pathogens - from ancestors that synthesized Q de novo to a state prioritizing the salvage of q. Another possibility that will require further investigations is that the insertion of preQ1 confers fitness advantages when B. henselae is growing outside a mammalian host.

Francisella sciaenopsi sp. nov. isolated from diseased red drum Sciaenops ocellatus in Florida, USA

Piscine francisellosis is one of the most important bacterial diseases affecting various fish species worldwide. Francisella orientalis, F. noatunensis, and F. salimarina (F. marina) have been reported as etiological agents of disease in fish. A Francisella sp. was isolated from several diseased red drum Sciaenops ocellatus experiencing morbidity in Florida, USA, in 2008. In this study, molecular and phenotypic characterization of the recovered isolate was conducted. Phenotypically, the isolate showed a biochemical reaction profile distinct from that of F. orientalis and F. salimarina. Although the 16S rRNA sequence of this isolate shared 99.61% identity to the type strain of F. philomiragia O#319LT, whole genome analysis (average nucleotide identity <95%; digital DNA-DNA hybridization <70%) and a multilocus sequence analysis of 8 concatenated housekeeping genes in comparison with other Francisella spp. indicated that this isolate was a novel Francisella species, more closely related to F. orientalis. Immersion, intracoelomic injection, and co-habitation challenges using a Nile tilapia Oreochromis niloticus fingerling model of infection were done to investigate virulence in a piscine model. Variably pigmented granulomas and pigmented macrophage aggregates were observed in the kidneys and spleens of the challenged fish, but no mortality was recorded during the 15 d challenge period, suggesting that this novel Francisella sp. might be an opportunistic pathogen of fish. Based on the phenotypic and genotypic differences from other Francisella spp. observed in this study, we propose the name Francisella sciaenopsi sp. nov. for this novel isolate.

Convergent reductive evolution in bee-associated lactic acid bacteria

Distantly related organisms may evolve similar traits when exposed to similar environments or engaging in certain lifestyles. Several members of the Lactobacillaceae (LAB) family are frequently isolated from the floral niche, mostly from bees and flowers. In some floral LAB species (henceforth referred to as bee-associated), distinctive genomic (e.g., genome reduction) and phenotypic (e.g., preference for fructose over glucose or fructophily) features were recently documented. These features are found across distantly related species, raising the hypothesis that specific genomic and phenotypic traits evolved convergently during adaptation to the floral environment. To test this hypothesis, we examined representative genomes of 369 species of bee-associated and non-bee-associated LAB. Phylogenomic analysis unveiled seven independent ecological shifts towards the floral niche in LAB. In these bee-associated LAB, we observed pervasive, significant reductions of genome size, gene repertoire, and GC content. Using machine leaning, we could distinguish bee-associated from non-bee-associated species with 94% accuracy, based on the absence of genes involved in metabolism, osmotic stress, or DNA repair. Moreover, we found that the most important genes for the machine learning classifier were seemingly lost, independently, in multiple bee-associated lineages. One of these genes, adhE, encodes a bifunctional aldehyde-alcohol dehydrogenase associated with the evolution of fructophily, a rare phenotypic trait that was recently identified in many floral LAB species. These results suggest that the independent evolution of distinctive phenotypes in bee-associated LAB has been largely driven by independent loss of the same set of genes.

Importance

Several lactic acid bacteria (LAB) species are intimately associated with bees and exhibit unique biochemical properties with potential for food applications and honeybee health. Using a machine-learning based approach, our study shows that adaptation of LAB to the bee environment was accompanied by a distinctive genomic trajectory deeply shaped by gene loss. Several of these gene losses occurred independently in distantly related species and are linked to some of their unique biotechnologically relevant traits, such as the preference of fructose over glucose (fructophily). This study underscores the potential of machine learning in identifying fingerprints of adaptation and detecting instances of convergent evolution. Furthermore, it sheds light onto the genomic and phenotypic particularities of bee-associated bacteria, thereby deepening the understanding of their positive impact on honeybee health.

On the evolution of chromosomal regions with high gene strand bias in bacteria

On circular bacterial chromosomes, the majority of genes are coded on the leading strand. This gene strand bias (GSB) ranges from up to 85% in some Bacillota to a little more than 50% in other phyla. The factors determining the extent of the strand bias remain to be found. Here, we report that species in the phylum Gemmatimonadota share a unique chromosome architecture, distinct from neighboring phyla: in a conserved 600-kb region around the terminus of replication, almost all genes were located on the leading strands, while on the remaining part of the chromosome, the strand preference was more balanced. The high strand bias (HSB) region harbors the rRNA clusters, core, and highly expressed genes. Selective pressure for reduction of collisions with DNA replication to minimize detrimental mutations can explain the conservation of essential genes in this region. Repetitive and mobile elements are underrepresented, suggesting reduced recombination frequency by structural isolation from other parts of the chromosome. We propose that the HSB region forms a distinct chromosomal domain. Gemmatimonadota chromosomes evolved mainly by expansion through horizontal gene transfer and duplications outside of the ancient high strand bias region. In support of our hypothesis, we could further identify two Spiroplasma strains on a similar evolutionary path.IMPORTANCEOn bacterial chromosomes, a preferred location of genes on the leading strand has evolved to reduce conflicts between replication and transcription. Despite a vast body of research, the question why bacteria show large differences in their gene strand bias is still not solved. The discovery of "hybrid" chromosomes in different phyla, including Gemmatimonadota, in which a conserved high strand bias is found exclusively in a region at ter, points toward a role of nucleoid structure, additional to replication, in the evolution of strand preferences. A fine-grained structural analysis of the ever-increasing number of available bacterial genomes could help to better understand the forces that shape the sequential and spatial organization of the cell's information content.

Bacterial lifestyle shapes pangenomes

Pangenomes vary across bacteria. Some species have fluid pangenomes, with a high proportion of genes varying between individual genomes. Other species have less fluid pangenomes, with different genomes tending to contain the same genes. Two main hypotheses have been suggested to explain this variation: differences in species' bacterial lifestyle and effective population size. However, previous studies have not been able to test between these hypotheses because the different features of lifestyle and effective population size are highly correlated with each other, and phylogenetically conserved, making it hard to disentangle their relative importance. We used phylogeny-based analyses, across 126 bacterial species, to tease apart the causal role of different factors. We found that pangenome fluidity was lower in i) host-associated compared with free-living species and ii) host-associated species that are obligately dependent on a host, live inside cells, and are more pathogenic and less motile. In contrast, we found no support for the competing hypothesis that larger effective population sizes lead to more fluid pangenomes. Effective population size appears to correlate with pangenome variation because it is also driven by bacterial lifestyle, rather than because of a causal relationship.

Queuosine Salvage in Bartonella henselae Houston 1: A Unique Evolutionary Path

Queuosine (Q) stands out as the sole tRNA modification that can be synthesized via salvage pathways. Comparative genomic analyses identified specific bacteria that showed a discrepancy between the projected Q salvage route and the predicted substrate specificities of the two identified salvage proteins: 1) the distinctive enzyme tRNA guanine-34 transglycosylase (bacterial TGT, or bTGT), responsible for inserting precursor bases into target tRNAs; and 2) Queuosine Precursor Transporter (QPTR), a transporter protein that imports Q precursors. Organisms like the facultative intracellular pathogen Bartonella henselae, which possess only bTGT and QPTR but lack predicted enzymes for converting preQ1 to Q, would be expected to salvage the queuine (q) base, mirroring the scenario for the obligate intracellular pathogen Chlamydia trachomatis. However, sequence analyses indicate that the substrate-specificity residues of their bTGTs resemble those of enzymes inserting preQ1 rather than q. Intriguingly, mass spectrometry analyses of tRNA modification profiles in B. henselae reveal trace amounts of preQ1, previously not observed in a natural context. Complementation analysis demonstrates that B. henselae bTGT and QPTR not only utilize preQ1, akin to their Escherichia coli counterparts, but can also process q when provided at elevated concentrations. The experimental and phylogenomic analyses suggest that the Q pathway in B. henselae could represent an evolutionary transition among intracellular pathogens-from ancestors that synthesized Q de novo to a state prioritizing the salvage of q. Another possibility that will require further investigations is that the insertion of preQ1 has fitness advantages when B. henselae is growing outside a mammalian host.

(Pan)genomic analysis of two Rhodococcus isolates and their role in phenolic compound degradation

The genus Rhodococcus is recognized for its potential to degrade a large range of aromatic substances, including plant-derived phenolic compounds. We used comparative genomics in the context of the broader Rhodococcus pan-genome to study genomic traits of two newly described Rhodococcus strains (type-strain Rhodococcus pseudokoreensis R79T and Rhodococcus koreensis R85) isolated from apple rhizosphere. Of particular interest was their ability to degrade phenolic compounds as part of an integrated approach to treat apple replant disease (ARD) syndrome. The pan-genome of the genus Rhodococcus based on 109 high-quality genomes was open with a small core (1.3%) consisting of genes assigned to basic cell functioning. The range of genome sizes in Rhodococcus was high, from 3.7 to 10.9 Mbp. Genomes from host-associated strains were generally smaller compared to environmental isolates which were characterized by exceptionally large genome sizes. Due to large genomic differences, we propose the reclassification of distinct groups of rhodococci like the Rhodococcus equi cluster to new genera. Taxonomic species affiliation was the most important factor in predicting genetic content and clustering of the genomes. Additionally, we found genes that discriminated between the strains based on habitat. All members of the genus Rhodococcus had at least one gene involved in the pathway for the degradation of benzoate, while biphenyl degradation was mainly restricted to strains in close phylogenetic relationships with our isolates. The ~40% of genes still unclassified in larger Rhodococcus genomes, particularly those of environmental isolates, need more research to explore the metabolic potential of this genus.IMPORTANCERhodococcus is a diverse, metabolically powerful genus, with high potential to adapt to different habitats due to the linear plasmids and large genome sizes. The analysis of its pan-genome allowed us to separate host-associated from environmental strains, supporting taxonomic reclassification. It was shown which genes contribute to the differentiation of the genomes based on habitat, which can possibly be used for targeted isolation and screening for desired traits. With respect to apple replant disease (ARD), our isolates showed genome traits that suggest potential for application in reducing plant-derived phenolic substances in soil, which makes them good candidates for further testing against ARD.

Differentiation of hypervirulent and classical Klebsiella pneumoniae with acquired drug resistance

Distinguishing hypervirulent (hvKp) from classical Klebsiella pneumoniae (cKp) strains is important for clinical care, surveillance, and research. Some combinations of iucA, iroB, peg-344, rmpA, and rmpA2 are most commonly used, but it is unclear what combination of genotypic or phenotypic markers (e.g., siderophore concentration, mucoviscosity) most accurately predicts the hypervirulent phenotype. Furthermore, acquisition of antimicrobial resistance may affect virulence and confound identification. Therefore, 49 K. pneumoniae strains that possessed some combinations of iucA, iroB, peg-344, rmpA, and rmpA2 and had acquired resistance were assembled and categorized as hypervirulent hvKp (hvKp) (N = 16) or cKp (N = 33) via a murine infection model. Biomarker number, siderophore production, mucoviscosity, virulence plasmid's Mash/Jaccard distances to the canonical pLVPK, and Kleborate virulence score were measured and evaluated to accurately differentiate these pathotypes. Both stepwise logistic regression and a CART model were used to determine which variable was most predictive of the strain cohorts. The biomarker count alone was the strongest predictor for both analyses. For logistic regression, the area under the curve for biomarker count was 0.962 (P = 0.004). The CART model generated the classification rule that a biomarker count = 5 would classify the strain as hvKP, resulting in a sensitivity for predicting hvKP of 94% (15/16), a specificity of 94% (31/33), and an overall accuracy of 94% (46/49). Although a count of ≥4 was 100% (16/16) sensitive for predicting hvKP, the specificity and accuracy decreased to 76% (25/33) and 84% (41/49), respectively. These findings can be used to inform the identification of hvKp.IMPORTANCEHypervirulent Klebsiella pneumoniae (hvKp) is a concerning pathogen that can cause life-threatening infections in otherwise healthy individuals. Importantly, although strains of hvKp have been acquiring antimicrobial resistance, the effect on virulence is unclear. Therefore, it is of critical importance to determine whether a given antimicrobial resistant K. pneumoniae isolate is hypervirulent. This report determined which combination of genotypic and phenotypic markers could most accurately identify hvKp strains with acquired resistance. Both logistic regression and a machine-learning prediction model demonstrated that biomarker count alone was the strongest predictor. The presence of all five of the biomarkers iucA, iroB, peg-344, rmpA, and rmpA2 was most accurate (94%); the presence of ≥4 of these biomarkers was most sensitive (100%). Accurately identifying hvKp is vital for surveillance and research, and the availability of biomarker data could alert the clinician that hvKp is a consideration, which, in turn, would assist in optimizing patient care.

PoMeLo: a systematic computational approach to predicting metabolic loss in pathogen genomes

BACKGROUND

Genome streamlining, the process by which genomes become smaller and encode fewer genes over time, is a common phenomenon among pathogenic bacteria. This reduction is driven by selection for minimized energy expenditure in a nutrient-rich environment. As pathogens evolve to become more reliant on the host, metabolic genes and resulting capabilities are lost in favor of siphoning metabolites from the host. Characterizing genome streamlining, gene loss, and metabolic pathway degradation can be useful in assessing pathogen dependency on host metabolism and identifying potential targets for host-directed therapeutics.

RESULTS

PoMeLo (Predictor of Metabolic Loss) is a novel evolutionary genomics-guided computational approach for identifying metabolic gaps in the genomes of pathogenic bacteria. PoMeLo leverages a centralized public database of high-quality genomes and annotations and allows the user to compare an unlimited number of genomes across individual genes and pathways. PoMeLo runs locally using user-friendly prompts in a matter of minutes and generates tabular and visual outputs for users to compare predicted metabolic capacity between groups of bacteria and individual species. Each pathway is assigned a Predicted Metabolic Loss (PML) score to assess the magnitude of genome streamlining. Optionally, PoMeLo places the results in an evolutionary context by including phylogenetic relationships in visual outputs. It can also initially compute phylogenetically-weighted mean genome sizes to identify genome streamlining events. Here, we describe PoMeLo and demonstrate its use in identifying metabolic gaps in genomes of pathogenic Treponema species.

CONCLUSIONS

PoMeLo represents an advance over existing methods for identifying metabolic gaps in genomic data, allowing comparison across large numbers of genomes and placing the resulting data in a phylogenetic context. PoMeLo is freely available for academic and non-academic use at https://github.com/czbiohub-sf/pomelo .

The streptococcal phase-variable type I restriction modification system SsuCC20p dictates the methylome of Streptococcus suis impacting the transcriptome and virulence in a zebrafish larvae infection model

IMPORTANCE

Phase variation allows a single strain to produce phenotypic diverse subpopulations. Phase-variable restriction modification (RM) systems are systems that allow for such phase variation via epigenetic regulation of gene expression levels. The phase-variable RM system SsuCC20p was found in multiple streptococcal species and was acquired by an emerging zoonotic lineage of Streptococcus suis. We show that the phase variability of SsuCC20p is dependent on a recombinase encoded within the SsuCC20p locus. We characterized the genome methylation profiles of the different phases of SsuCC20p and demonstrated the consequential impact on the transcriptome and virulence in a zebrafish infection model. Acquiring mobile genetic elements containing epigenetic regulatory systems, like phase-variable RM systems, enables bacterial pathogens to produce diverse phenotypic subpopulations that are better adapted to specific (host) environments encountered during infection.

Escherichia coli O157:H7 tir 255 T > A allele strains differ in chromosomal and plasmid composition

Shiga toxin-producing Escherichia coli (STEC) O157:H7 strains with the T allele in the translocated intimin receptor polymorphism (tir) 255 A > T gene associate with human disease more than strains with an A allele; however, the allele is not thought to be the direct cause of this difference. We sequenced a diverse set of STEC O157:H7 strains (26% A allele, 74% T allele) to identify linked differences that might underlie disease association. The average chromosome and pO157 plasmid size and gene content were significantly greater within the tir 255 A allele strains. Eighteen coding sequences were unique to tir 255 A allele chromosomes, and three were unique to tir 255 T allele chromosomes. There also were non-pO157 plasmids that were unique to each tir 255 allele variant. The overall average number of prophages did not differ between tir 255 allele strains; however, there were different types between the strains. Genomic and mobile element variation linked to the tir 255 polymorphism may account for the increased frequency of the T allele isolates in human disease.

The genomes of Scedosporium between environmental challenges and opportunism

Emerging fungal pathogens are a global challenge for humankind. Many efforts have been made to understand the mechanisms underlying pathogenicity in bacteria, and OMICs techniques are largely responsible for those advancements. By contrast, our limited understanding of opportunism and antifungal resistance is preventing us from identifying, limiting and interpreting the emergence of fungal pathogens. The genus Scedosporium (Microascaceae) includes fungi with high tolerance to environmental pollution, whilst some species can be considered major human pathogens, such as Scedosporium apiospermum and Scedosporium boydii. However, unlike other fungal pathogens, little is known about the genome evolution of these organisms. We sequenced two novel genomes of Scedosporium aurantiacum and Scedosporium minutisporum isolated from extreme, strongly anthropized environments. We compared all the available Scedosporium and Microascaceae genomes, that we systematically annotated and characterized ex novo in most cases. The genomes in this family were integrated in a Phylum-level comparison to infer the presence of putative, shared genomic traits in filamentous ascomycetes with pathogenic potential. The analysis included the genomes of 100 environmental and clinical fungi, revealing poor evolutionary convergence of putative pathogenicity traits. By contrast, several features in Microascaceae and Scedosporium were detected that might have a dual role in responding to environmental challenges and allowing colonization of the human body, including chitin, melanin and other cell wall related genes, proteases, glutaredoxins and magnesium transporters. We found these gene families to be impacted by expansions, orthologous transposon insertions, and point mutations. With RNA-seq, we demonstrated that most of these anciently impacted genomic features responded to the stress imposed by an antifungal compound (voriconazole) in the two environmental strains S. aurantiacum MUT6114 and S. minutisporum MUT6113. Therefore, the present genomics and transcriptomics investigation stands on the edge between stress resistance and pathogenic potential, to elucidate whether fungi were pre-adapted to infect humans. We highlight the strengths and limitations of genomics applied to opportunistic human pathogens, the multifactoriality of pathogenicity and resistance to drugs, and suggest a scenario where pressures other than anthropic contributed to forge filamentous human pathogens.

The emergence and diversification of a zoonotic pathogen from within the microbiota of intensively farmed pigs

The expansion and intensification of livestock production is predicted to promote the emergence of pathogens. As pathogens sometimes jump between species, this can affect the health of humans as well as livestock. Here, we investigate how livestock microbiota can act as a source of these emerging pathogens through analysis of Streptococcus suis, a ubiquitous component of the respiratory microbiota of pigs that is also a major cause of disease on pig farms and an important zoonotic pathogen. Combining molecular dating, phylogeography, and comparative genomic analyses of a large collection of isolates, we find that several pathogenic lineages of S. suis emerged in the 19th and 20th centuries, during an early period of growth in pig farming. These lineages have since spread between countries and continents, mirroring trade in live pigs. They are distinguished by the presence of three genomic islands with putative roles in metabolism and cell adhesion, and an ongoing reduction in genome size, which may reflect their recent shift to a more pathogenic ecology. Reconstructions of the evolutionary histories of these islands reveal constraints on pathogen emergence that could inform control strategies, with pathogenic lineages consistently emerging from one subpopulation of S. suis and acquiring genes through horizontal transfer from other pathogenic lineages. These results shed light on the capacity of the microbiota to rapidly evolve to exploit changes in their host population and suggest that the impact of changes in farming on the pathogenicity and zoonotic potential of S. suis is yet to be fully realized.

Differentiation of hypervirulent and classical Klebsiella pneumoniae with acquired drug resistance

Distinguishing hypervirulent (hvKp) from classical Klebsiella pneumoniae (cKp) strains is important for clinical care, surveillance, and research. Some combination of iucA, iroB, peg-344, rmpA, and rmpA2 are most commonly used, but it is unclear what combination of genotypic or phenotypic markers (e.g. siderophore concentration, mucoviscosity) most accurately predicts the hypervirulent phenotype. Further, acquisition of antimicrobial resistance may affect virulence and confound identification. Therefore, 49 K. pneumoniae strains that possessed some combination of iucA, iroB, peg-344, rmpA, and rmpA2 and had acquired resistance were assembled and categorized as hypervirulent hvKp (hvKp) (N=16) or cKp (N=33) via a murine infection model. Biomarker number, siderophore production, mucoviscosity, virulence plasmid's Mash/Jaccard distances to the canonical pLVPK, and Kleborate virulence score were measured and evaluated to accurately differentiate these pathotypes. Both stepwise logistic regression and a CART model were used to determine which variable was most predictive of the strain cohorts. The biomarker count alone was the strongest predictor for both analyses. For logistic regression the area under the curve for biomarker count was 0.962 (P = 0.004). The CART model generated the classification rule that a biomarker count = 5 would classify the strain as hvKP, resulting in a sensitivity for predicting hvKP of 94% (15/16), a specificity of 94% (31/33), and an overall accuracy of 94% (46/49). Although a count of ≥ 4 was 100% (16/16) sensitive for predicting hvKP, the specificity and accuracy decreased to 76% (25/33) and 84% (41/49) respectively. These findings can be used to inform the identification of hvKp.

Importance

Hypervirulent Klebsiella pneumoniae (hvKp) is a concerning pathogen that can cause life-threatening infections in otherwise healthy individuals. Importantly, although strains of hvKp have been acquiring antimicrobial resistance, the effect on virulence is unclear. Therefore, it is of critical importance to determine whether a given antimicrobial resistant K. pneumoniae isolate is hypervirulent. This report determined which combination of genotypic and phenotypic markers could most accurately identify hvKp strains with acquired resistance. Both logistic regression and a machine-learning prediction model demonstrated that biomarker count alone was the strongest predictor. The presence of all 5 of the biomarkers iucA, iroB, peg-344, rmpA, and rmpA2 was most accurate (94%); the presence of ≥ 4 of these biomarkers was most sensitive (100%). Accurately identifying hvKp is vital for surveillance and research, and the availability of biomarker data could alert the clinician that hvKp is a consideration, which in turn would assist in optimizing patient care.

Exploring the biosynthetic gene clusters in Brevibacterium: a comparative genomic analysis of diversity and distribution

Exploring Brevibacterium strains from various ecosystems may lead to the discovery of new antibiotic-producing strains. Brevibacterium sp. H-BE7, a strain isolated from marine sediments from Northern Patagonia, Chile, had its genome sequenced to study the biosynthetic potential to produce novel natural products within the Brevibacterium genus. The genome sequences of 98 Brevibacterium strains, including strain H-BE7, were selected for a genomic analysis. A phylogenomic cladogram was generated, which divided the Brevibacterium strains into four major clades. A total of 25 strains are potentially unique new species according to Average Nucleotide Identity (ANIb) values. These strains were isolated from various environments, emphasizing the importance of exploring diverse ecosystems to discover the full diversity of Brevibacterium. Pangenome analysis of Brevibacterium strains revealed that only 2.5% of gene clusters are included within the core genome, and most gene clusters occur either as singletons or as cloud genes present in less than ten strains. Brevibacterium strains from various phylogenomic clades exhibit diverse BGCs. Specific groups of BGCs show clade-specific distribution patterns, such as siderophore BGCs and carotenoid-related BGCs. A group of clade IV-A Brevibacterium strains possess a clade-specific Polyketide synthase (PKS) BGCs that connects with phenazine-related BGCs. Within the PKS BGC, five genes, including the biosynthetic PKS gene, participate in the mevalonate pathway and exhibit similarities with the phenazine A BGC. However, additional core biosynthetic phenazine genes were exclusively discovered in nine Brevibacterium strains, primarily isolated from cheese. Evaluating the antibacterial activity of strain H-BE7, it exhibited antimicrobial activity against Salmonella enterica and Listeria monocytogenes. Chemical dereplication identified bioactive compounds, such as 1-methoxyphenazine in the crude extracts of strain H-BE7, which could be responsible of the observed antibacterial activity. While strain H-BE7 lacks the core phenazine biosynthetic genes, it produces 1-methoxyphenazine, indicating the presence of an unknown biosynthetic pathway for this compound. This suggests the existence of alternative biosynthetic pathways or promiscuous enzymes within H-BE7's genome.

Genomic insights and comparative analysis of Flavobacterium bizetiae HJ-32-4 isolated from soil

In the late 1970s, Flavobacterium bizetiae was first isolated from diseased fish in Canada. After four decades of preservation, it was reported as a novel species in 2020. Here, we report the first complete genome sequence of HJ-32-4, a novel strain of F. bizetiae. Interestingly, HJ-32-4 was isolated from soil in Gangwon-do, Republic of Korea, unlike the other two previously reported F. bizetiae strains which were isolated from fish. We generated a single circular chromosome of HJ-32-4, comprising 5,745,280 bp with a GC content of 34.2%. The average nucleotide identity (ANI) value of 96.2% indicated that HJ-32-4 belongs to F. bizetiae CIP 105534T. The virulence factor was not detected in the genome. Comparative genomic analysis of F. bizetiae and major flavobacterial pathogens revealed that F. bizetiae had a larger genome size and the ratio of peptidases (PEP) and glycoside hydrolase (GH) genes of F. bizetiae were lower than those of the rest strains, implying that F. bizetiae exhibits similar characteristics with non-pathogenic strains from a genomic point of view. However, further experimental verification is required to ensure these in silico predictions. This study will provide insight into the overall characteristics of HJ-32-4 compared to other strains.

Genome characteristics related to the virulence of Streptococcus suis in Swedish pigs

The impact of S. suis on Swedish pig production has increased in recent years, and characterization of the strains present in the pig population is needed to aid in surveillance and prevention. Therefore, the aim of this study was to identify and characterize differences in the genomes between Swedish S. suis isolates associated with disease and isolates from healthy animals. Isolates categorized as being pathogenic (n = 100) or non-pathogenic (n = 117) were whole-genome sequenced, serotyped in silico, and sequence-typed using traditional MLST and core-genome MLST, and a genome-wide association study was performed to identify virulence-associated genes. In decreasing order, serotypes 2, 1, and 7 were the most common in the pathogenic group, and serotypes 15 and 12 were the most common in the non-pathogenic group. Among the commonly disease-associated sequence types, ST28 and ST25 were identified, whereas ST1 was scarcely found. The majority of isolates belonged to novel sequence types, revealing differences between Swedish isolates and those reported from other countries. The genomes of the pathogenic isolates were on average smaller and less heterogenic as compared to those of the non-pathogenic isolates. Although a majority of the previously published virulence-associated genes included in the study were found in the genomes of both pathogenic and non-pathogenic isolates, several new, significantly virulence-associated genes were identified.

Arcobacteraceae comparative genome analysis demonstrates genome heterogeneity and reduction in species isolated from animals and associated with human illness

The Arcobacteraceae family groups Gram-negative bacterial species previously included in the family Campylobacteraceae. These species of which some are considered foodborne pathogens, have been isolated from different environmental niches and hosts. They have been isolated from various types of foods, though predominantly from food of animal origin, as well as from stool of humans with enteritis. Their different abilities to survive in different hosts and environments suggest an evolutionary pressure with consequent variation in their genome content. Moreover, their different physiological and genomic characteristics led to the recent proposal to create new genera within this family, which is however criticized due to the lack of discriminatory features and biological and clinical relevance. Aims of the present study were to assess the Arcobacteraceae pangenome, and to characterize existing similarities and differences in 20 validly described species. For this, analysis has been conducted on the genomes of the corresponding type strains obtained by Illumina sequencing, applying several bioinformatic tools. Results of the present study do not support the proposed division into different genera and revealed the presence of pangenome partitions with numbers comparable to other Gram-negative bacteria genera, such as Campylobacter. Different gene class compositions in animal and human-associated species are present, including a higher percentage of virulence-related gene classes such as cell motility genes. The adaptation to environmental and/or host conditions of some species was identified by the presence of specific genes. Furthermore, a division into pathogenic and non-pathogenic species is suggested, which can support future research on food safety and public health.

Comparative genomics of the Liberibacter genus reveals widespread diversity in genomic content and positive selection history

'Candidatus Liberibacter' is a group of bacterial species that are obligate intracellular plant pathogens and cause Huanglongbing disease of citrus trees and Zebra Chip in potatoes. Here, we examined the extent of intra- and interspecific genetic diversity across the genus using comparative genomics. Our approach examined a wide set of Liberibacter genome sequences including five pathogenic species and one species not known to cause disease. By performing comparative genomics analyses, we sought to understand the evolutionary history of this genus and to identify genes or genome regions that may affect pathogenicity. With a set of 52 genomes, we performed comparative genomics, measured genome rearrangement, and completed statistical tests of positive selection. We explored markers of genetic diversity across the genus, such as average nucleotide identity across the whole genome. These analyses revealed the highest intraspecific diversity amongst the 'Ca. Liberibacter solanacearum' species, which also has the largest plant host range. We identified sets of core and accessory genes across the genus and within each species and measured the ratio of nonsynonymous to synonymous mutations (dN/dS) across genes. We identified ten genes with evidence of a history of positive selection in the Liberibacter genus, including genes in the Tad complex, which have been previously implicated as being highly divergent in the 'Ca. L. capsica' species based on high values of dN.

Characterization of Shigella flexneri Serotype 6 Strains Isolated from Bangladesh and Identification of a New Phylogenetic Cluster

A significant cause of shigellosis in Bangladesh and other developing countries is Shigella flexneri serotype 6. This serotype has been subtyped, on the basis of the absence or presence of a group-specific antigen, E1037, into S. flexneri 6a and 6b, respectively. Here, we provided rationales for the subclassification, using several phenotypic and molecular tools. A set of S. flexneri 6a and 6b strains isolated between 1997 and 2015 were characterized by analyzing their biochemical properties, plasmid profiles, virulence markers, pulsed-field gel electrophoresis (PFGE) results, and ribotype. Additionally, the genomic relatedness of these subserotypes was investigated with global isolates of serotype 6 using publicly available genomes. Both subserotypes of S. flexneri 6 agglutinated with monoclonal antiserum against S. flexneri (MASF) B and type VI-specific antiserum (MASF VI) and were PCR positive for O-antigen flippase-specific genes and virulence markers (ipaH, ial, sen, and sigA). Unlike S. flexneri 6a strains, S. flexneri 6b strains seroagglutinated with anti-E1037 antibodies, MASF IV-I. Notably, these two antigenically distinct subserotypes were clonally diverse, showing two distinct PFGE patterns following the digestion of chromosomal DNA with either XbaI or IceuI. In addition, hybridization of a 16S rRNA gene probe with HindIII-digested genomic DNA yielded two distinguishing ribotypes. Genomic comparison of S. flexneri subserotype 6a and 6b strains from Bangladesh indicated that, although these strains were in genomic synteny, the majority of them formed a unique phylogroup (PG-4) that was missing for the global isolates. This study supports the subserotyping and emphasizes the need for global monitoring of the S. flexneri subserotypes 6a and 6b. IMPORTANCE Shigella flexneri serotype 6 is one of the predominant serotypes among shigellosis cases in Bangladesh. Characterization of a novel subserotype of S. flexneri 6 (VI:E1037), agglutinated with type 6-specific antibody and anti-E1037, indicates a unique evolutionary ancestry. PFGE genotyping supports the finding that these two antigenically distinct subserotypes are clonally diverse. A phylogenetic study based on single-nucleotide polymorphism (SNP) data revealed that these two subserotypes were in genomic synteny, although their genomes were reduced. Interestingly, a majority of the S. flexneri 6 strains isolated from Bangladesh form a novel phylogenetic cluster. Therefore, this report underpins the global monitoring and tracking of the novel subserotype.

Transcriptomics in serum and culture medium reveal shared and differential gene regulation in pathogenic and commensal Streptococcus suis

Streptococcus suis colonizes the upper respiratory tract of healthy pigs at high abundance but can also cause opportunistic respiratory and systemic disease. Disease-associated S. suis reference strains are well studied, but less is known about commensal lineages. It is not known what mechanisms enable some S. suis lineages to cause disease while others persist as commensal colonizers, or to what extent gene expression in disease-associated and commensal lineages diverge. In this study we compared the transcriptomes of 21 S. suis strains grown in active porcine serum and Todd-Hewitt yeast broth. These strains included both commensal and pathogenic strains, including several strains of sequence type (ST) 1, which is responsible for most cases of human disease and is considered to be the most pathogenic S. suis lineage. We sampled the strains during their exponential growth phase and mapped RNA sequencing reads to the corresponding strain genomes. We found that the transcriptomes of pathogenic and commensal strains with large genomic divergence were unexpectedly conserved when grown in active porcine serum, but that regulation and expression of key pathways varied. Notably, we observed strong variation of expression across media of genes involved in capsule production in pathogens, and of the agmatine deiminase system in commensals. ST1 strains displayed large differences in gene expression between the two media compared to strains from other clades. Their capacity to regulate gene expression across different environmental conditions may be key to their success as zoonotic pathogens.

Role of Metabolic Adaptation of Streptococcus suis to Host Niches in Bacterial Fitness and Virulence

Streptococcus suis, both a common colonizer of the porcine upper respiratory tract and an invasive pig pathogen, successfully adapts to different host environments encountered during infection. Whereas the initial infection mainly occurs via the respiratory tract, in a second step, the pathogen can breach the epithelial barrier and disseminate within the whole body. Thereby, the pathogen reaches other organs such as the heart, the joints, or the brain. In this review, we focus on the role of S. suis metabolism for adaptation to these different in vivo host niches to encounter changes in nutrient availability, host defense mechanisms and competing microbiota. Furthermore, we highlight the close link between S. suis metabolism and virulence. Mutants deficient in metabolic regulators often show an attenuation in infection experiments possibly due to downregulation of virulence factors, reduced resistance to nutritive or oxidative stress and to phagocytic activity. Finally, metabolic pathways as potential targets for new therapeutic strategies are discussed. As antimicrobial resistance in S. suis isolates has increased over the last years, the development of new antibiotics is of utmost importance to successfully fight infections in the future.

Group B streptococcus virulence factors associated with different clinical syndromes: Asymptomatic carriage in pregnant women and early-onset disease in the newborn

Background

Group B streptococcus (GBS) harbors many virulence factors but there is limited data regarding their importance in colonization in pregnancy and early-onset disease (EOD) in the newborn. We hypothesized that colonization and EOD are associated with different distribution and expression of virulence factors.

Methods

We studied 36 GBS EOD and 234 GBS isolates collected during routine screening. Virulence genes (pilus-like structures-PI-1, PI-2a, PI-2b; rib and hvgA) presence and expression were identified by PCR and qRT-PCR. Whole genome sequencing (WGS) and comparative genomic analyses were used to compare coding sequences (CDSs) of colonizing and EOD isolates.

Results

Serotype III (ST17) was significantly associated with EOD and serotype VI (ST1) with colonization. hvgA and rib genes were more prevalent among EOD isolates (58.3 and 77.8%, respectively; p < 0.01). The pilus loci PI-2b and PI-2a were more prevalent among EOD isolates (61.1%, p < 0.01), while the pilus loci PI-2a and PI-1 among colonizing isolates (89.7 and 93.1% vs. 55.6 and 69.4%, p < 0.01). qRT PCR analysis revealed that hvgA was barely expressed in colonizing isolates, even though the gene was detected. Expression of the rib gene and PI-2b was two-fold higher in EOD isolates compared to colonizing isolates. Transcription of PI-2a was three-fold higher in colonizing isolates compared to EOD isolates. ST17 isolates (associated with EOD) had a smaller genome size compared ST1 and the genome was more conserved relative to the reference strain and ST17 isolates. In a multivariate logistic regression analysis virulence factors independently associated with EOD were serotype 3, and PI-1 and PI-2a was protective.

Conclusion

There was a significant difference in the distribution of hvg A, rib, and PI genes among EOD (serotype III/ST17) and colonizing (serotype VI/ST1) isolates suggesting an association between invasive disease and these virulence factors. Further study is needed to understand the contribution of these genes to GBS virulence.

Robust analysis of prokaryotic pangenome gene gain and loss rates with Panstripe

Horizontal gene transfer (HGT) plays a critical role in the evolution and diversification of many microbial species. The resulting dynamics of gene gain and loss can have important implications for the development of antibiotic resistance and the design of vaccine and drug interventions. Methods for the analysis of gene presence/absence patterns typically do not account for errors introduced in the automated annotation and clustering of gene sequences. In particular, methods adapted from ecological studies, including the pangenome gene accumulation curve, can be misleading as they may reflect the underlying diversity in the temporal sampling of genomes rather than a difference in the dynamics of HGT. Here, we introduce Panstripe, a method based on generalized linear regression that is robust to population structure, sampling bias, and errors in the predicted presence/absence of genes. We show using simulations that Panstripe can effectively identify differences in the rate and number of genes involved in HGT events, and illustrate its capability by analyzing several diverse bacterial genome data sets representing major human pathogens.

The compact genome of Caenorhabditis niphades n. sp., isolated from a wood-boring weevil, Niphades variegatus

BACKGROUND

The first metazoan genome sequenced, that of Caenorhabditis elegans, has motivated animal genome evolution studies. To date > 50 species from the genus Caenorhabditis have been sequenced, allowing research on genome variation.

RESULTS

In the present study, we describe a new gonochoristic species, Caenorhabditis niphades n. sp., previously referred as C. sp. 36, isolated from adult weevils (Niphades variegatus), with whom they appear to be tightly associated during its life cycle. Along with a species description, we sequenced the genome of C. niphades n. sp. and produced a chromosome-level assembly. A genome comparison highlighted that C. niphades n. sp. has the smallest genome (59 Mbp) so far sequenced in the Elegans supergroup, despite being closely related to a species with an exceptionally large genome, C. japonica.

CONCLUSIONS

The compact genome of C. niphades n. sp. can serve as a key resource for comparative evolutionary studies of genome and gene number expansions in Caenorhabditis species.

Biochemical Atypia in Russian Neisseria gonorrhoeae Clinical Isolates Belonging to the G807 NG-MAST Genogroup/ST1594 MLST

Many current gonococcal clinical isolates in Russia show atypical taxonomically significant biochemical activity, which leads to species misidentification. Molecular typing of such cultures according Neisseria gonorrhoeae multiantigen sequence typing (NG-MAST) and multilocus sequence typing (MLST) protocols assigned them to the G807 NG-MAST GENOGROUP/ST1594 MLST that has been predominant in Russia in recent years. The goal of the study was to analyze the molecular mechanisms of biochemical atypia in N. gonorrhoeae clinical isolates characterized as the members of G807 NG-MAST GENOGROUP/ST1594 MLST. Sixteen isolates of this genogroup were included in the study, eight showed defective amino acid metabolism or loss of D-glucose fermentation. Comparative bioinformatic analysis based on WGS data divided these isolates into two clusters strictly associated with typical or atypical biochemical activity. Cultures with defective amino acid metabolism had a 5-nucleotide insertion in the pip-gene that caused a stop codon and led to synthesis of the non-functional enzyme. Comparison of the sequenced genomes with publicly available N. gonorrhoeae genomes showed the rarity of this insertion. In the global N. gonorrhoeae phylogenetic tree the G807 NG-MAST GENOGROUP/ST1594 MLST forms a distinct branch characterized by 170 SNPs, most of which are non-synonymous. We hypothesized a unique strategy for G807 NG-MAST GENOGROUP/ST1594 MLST clone persistence in the global N. gonorrhoeae population via escape of antimicrobial therapy due to diagnostic misidentification.

Genetic diversity and variation in antimicrobial-resistance determinants of non-serotype 2 Streptococcus suis isolates from healthy pigs

Streptococcus suis is a leading cause of bacterial meningitis in South-East Asia, with frequent zoonotic transfer to humans associated with close contact with pigs. A small number of invasive lineages are responsible for endemic infection in the swine industry, causing considerable global economic losses. A lack of surveillance and a rising trend in clinical treatment failure has raised concerns of growing antimicrobial resistance (AMR) among invasive S. suis. Gene flow between healthy and disease isolates is poorly understood and, in this study, we sample and sequence a collection of isolates predominantly from healthy pigs in Chiang Mai province, Northern Thailand. Pangenome characterization identified extensive genetic diversity and frequent AMR carriage in isolates from healthy pigs. Multiple AMR genes were identified, conferring resistance to aminoglycosides, lincosamides, tetracycline and macrolides. All isolates were non-susceptible to three or more different antimicrobial classes, and 75 % of non-serotype 2 isolates were non-susceptible to six or more classes (compared to 37.5 % of serotype 2 isolates). AMR genes were found on integrative and conjugative elements previously observed in other species, suggesting a mobile gene pool that can be accessed by invasive disease isolates. This article contains data hosted by Microreact.

Structure and mechanism for streptococcal fatty acid kinase (Fak) system dedicated to host fatty acid scavenging

Staphylococcus and Streptococcus, two groups of major human pathogens, are equipped with a fatty acid kinase (Fak) machinery to scavenge host fatty acids. The Fak complex is contains an ATP-binding subunit FakA, which interacts with varied FakB isoforms, and synthesizes acyl-phosphate from extracellular fatty acids. However, how FakA recognizes its FakB partners and then activates different fatty acids is poorly understood. Here, we systematically describe the Fak system from the zoonotic pathogen, Streptococcus suis. The crystal structure of SsFakA complexed with SsFakB2 was determined at 2.6 Å resolution. An in vitro system of Fak-PlsX (phosphate: acyl-ACP transacylase) was developed to track acyl-phosphate intermediate and its final product acyl-ACP. Structure-guided mutagenesis enabled us to characterize a mechanism for streptococcal FakA working with FakB partners engaged in host fatty acid scavenging. These findings offer a comprehensive description of the Fak kinase machinery, thus advancing the discovery of attractive targets against deadly infections with Streptococcus.

Rickettsia-Host-Tick Interactions: Knowledge Advances and Gaps

Ticks are hematophagous ectoparasites capable of transmitting multiple human pathogens. Environmental changes have supported the expansion of ticks into new geographical areas that have become the epicenters of tick-borne diseases (TBDs). The spotted fever group (SFG) of Rickettsia frequently infects ticks and causes tick-transmitted rickettsioses in areas of endemicity where ixodid ticks support host transmission during blood feeding. Ticks also serve as a reservoir for SFG Rickettsia. Among the members of SFG Rickettsia, R. rickettsii causes Rocky Mountain spotted fever (RMSF), the most lethal TBD in the United States. Cases of RMSF have been reported for over a century in association with several species of ticks in the United States. However, the isolation of R. rickettsii from ticks has decreased, and recent serological and epidemiological studies suggest that novel species of SFG Rickettsia are responsible for the increased number of cases of RMSF-like rickettsioses in the United States. Recent analyses of rickettsial genomes and advances in genetic and molecular studies of Rickettsia provided insights into the biology of Rickettsia with the identification of conserved and unique putative virulence genes involved in the rickettsial life cycle. Thus, understanding Rickettsia-host-tick interactions mediating successful disease transmission and pathogenesis for SFG rickettsiae remains an active area of research. This review summarizes recent advances in understanding how SFG Rickettsia species coopt and manipulate ticks and mammalian hosts to cause rickettsioses, with a particular emphasis on newly described or emerging SFG Rickettsia species.

Evaluation of existing guidelines for their adequacy for the food and feed risk assessment of microorganisms obtained through synthetic biology

EFSA was asked by the European Commission to evaluate synthetic biology (SynBio) developments for agri-food use in the near future and to determine whether or not they are expected to constitute potential new hazards/risks. Moreover, EFSA was requested to evaluate the adequacy of existing guidelines for risk assessment of SynBio and if updated guidance is needed. The scope of this Opinion covers food and feed risk assessment, the variety of microorganisms that can be used in the food/feed chain and the whole spectrum of techniques used in SynBio. This Opinion complements a previously adopted Opinion with the evaluation of existing guidelines for the microbial characterisation and environmental risk assessment of microorganisms obtained through SynBio. The present Opinion confirms that microbial SynBio applications for food and feed use, with the exception of xenobionts, could be ready in the European Union in the next decade. New hazards were identified related to the use or production of unusual and/or new-to-nature components. Fifteen cases were selected for evaluating the adequacy of existing guidelines. These were generally adequate for assessing the product, the production process, nutritional and toxicological safety, allergenicity, exposure and post-market monitoring. The comparative approach and a safety assessment per se could be applied depending on the degree of familiarity of the SynBio organism/product with the non-genetically modified counterparts. Updated guidance is recommended for: (i) bacteriophages, protists/microalgae, (ii) exposure to plant protection products and biostimulants, (iii) xenobionts and (iv) feed additives for insects as target species. Development of risk assessment tools is recommended for assessing nutritional value of biomasses, influence of microorganisms on the gut microbiome and the gut function, allergenic potential of new-to-nature proteins, impact of horizontal gene transfer and potential risks of living cell intake. A further development towards a strain-driven risk assessment approach is recommended.

Three enigmatic BioH isoenzymes are programmed in the early stage of mycobacterial biotin synthesis, an attractive anti-TB drug target

Tuberculosis (TB) is one of the leading infectious diseases of global concern, and one quarter of the world’s population are TB carriers. Biotin metabolism appears to be an attractive anti-TB drug target. However, the first-stage of mycobacterial biotin synthesis is fragmentarily understood. Here we report that three evolutionarily-distinct BioH isoenzymes (BioH1 to BioH3) are programmed in biotin synthesis of Mycobacterium smegmatis. Expression of an individual bioH isoform is sufficient to allow the growth of an Escherichia coli ΔbioH mutant on the non-permissive condition lacking biotin. The enzymatic activity in vitro combined with biotin bioassay in vivo reveals that BioH2 and BioH3 are capable of removing methyl moiety from pimeloyl-ACP methyl ester to give pimeloyl-ACP, a cognate precursor for biotin synthesis. In particular, we determine the crystal structure of dimeric BioH3 at 2.27Å, featuring a unique lid domain. Apart from its catalytic triad, we also dissect the substrate recognition of BioH3 by pimeloyl-ACP methyl ester. The removal of triple bioH isoforms (ΔbioH1/2/3) renders M. smegmatis biotin auxotrophic. Along with the newly-identified Tam/BioC, the discovery of three unusual BioH isoforms defines an atypical ‘BioC-BioH(3)’ paradigm for the first-stage of mycobacterial biotin synthesis. This study solves a long-standing puzzle in mycobacterial nutritional immunity, providing an alternative anti-TB drug target.

Construction and Analysis of the Complete Genome Sequence of Leprosy Agent Mycobacterium lepromatosis

Leprosy is caused by Mycobacterium leprae and Mycobacterium lepromatosis. We report construction and analyses of the complete genome sequence of M. lepromatosis FJ924. The genome contained 3,271,694 nucleotides to encode 1,789 functional genes and 1,564 pseudogenes. It shared 1,420 genes and 885 pseudogenes (71.4%) with M. leprae but differed in 1,281 genes and pseudogenes (28.6%). In phylogeny, the leprosy bacilli started from a most recent common ancestor (MRCA) that diverged ~30 million years ago (Mya) from environmental organism Mycobacterium haemophilum. The MRCA then underwent reductive evolution with pseudogenization, gene loss, and chromosomal rearrangements. Analysis of the shared pseudogenes estimated the pseudogenization event ~14 Mya, shortly before species bifurcation. Afterwards, genomic changes occurred to lesser extent in each species. Like M. leprae, four major types of highly repetitive sequences were detected in M. lepromatosis, contributing to chromosomal rearrangements within and after MRCA. Variations in genes and copy numbers were noted, such as three copies of the gene encoding bifunctional diguanylate cyclase/phosphodiesterase in M. lepromatosis, but single copy in M. leprae; 6 genes encoding the TetR family transcriptional regulators in M. lepromatosis, but 11 such genes in M. leprae; presence of hemW gene in M. lepromatosis, but absence in M. leprae; and others. These variations likely aid unique pathogenesis, such as diffuse lepromatous leprosy associated with M. lepromatosis, while the shared genomic features should explain the common pathogenesis of dermatitis and neuritis in leprosy. Together, these findings and the genomic data of M. lepromatosis may facilitate future research and care for leprosy. IMPORTANCE Leprosy is a dreaded infection that still affects millions of people worldwide. Mycobacterium lepromatosis is a recently recognized cause in addition to the well-known Mycobacterium leprae. M. lepromatosis is likely specific for diffuse lepromatous leprosy, a severe form of the infection and endemic in Mexico. This study constructed and annotated the complete genome sequence of M. lepromatosis FJ924 and performed comparative genomic analyses with related mycobacteria. The results afford new and refined insights into the genome size, gene repertoire, pseudogenes, phylogenomic relationship, genome organization and plasticity, process and timing of reductive evolution, and genetic and proteomic basis for pathogenesis. The availability of the complete M. lepromatosis genome may prove to be useful for future research and care for the infection.

Genomic characterization of Streptococcus parasuis, a close relative of Streptococcus suis and also a potential opportunistic zoonotic pathogen

Streptococcus parasuis (S. parasuis) is a close relative of Streptococcus suis (S. suis), composed of former members of S. suis serotypes 20, 22 and 26. S. parasuis could infect pigs and cows, and recently, human infection cases have been reported, making S. parasuis a potential opportunistic zoonotic pathogen. In this study, we analysed the genomic characteristics of S. parasuis, using pan-genome analysis, and compare some phenotypic determinants such as capsular polysaccharide, integrative conjugative elements, CRISPR-Cas system and pili, and predicted the potential virulence genes by associated analysis of the clinical condition of isolated source animals and genotypes. Furthermore, to discuss the relationship with S. suis, we compared these characteristics of S. parasuis with those of S. suis. We found that the characteristics of S. parasuis are similar to those of S. suis, both of them have "open" pan-genome, their antimicrobial resistance gene profiles are similar and a srtF pilus cluster of S. suis was identified in S. parasuis genome. But S. parasuis still have its unique characteristics, two novel pilus clusters are and three different type CRISPR-Cas system were found. Therefore, this study provides novel insights into the interspecific and intraspecific genetic characteristics of S. parasuis, which can be useful for further study of this opportunistic pathogen, such as serotyping, diagnostics, vaccine development, and study of the pathogenesis mechanism.

Evolution of Tuberculosis Pathogenesis

Mycobacterium tuberculosis is a globally distributed, lethal pathogen of humans. The virulence armamentarium of M. tuberculosis appears to have been developed on a scaffold of antiphagocytic defenses found among diverse, mostly free-living species of Mycobacterium. Pathoadaptation was further aided by the modularity, flexibility, and interactivity characterizing mycobacterial effectors and their regulators. During emergence of M. tuberculosis, novel genetic material was acquired, created, and integrated with existing tools. The major mutational mechanisms underlying these adaptations are discussed in this review, with examples. During its evolution, M. tuberculosis lost the ability and/or opportunity to engage in lateral gene transfer, but despite this it has retained the adaptability that characterizes mycobacteria. M. tuberculosis exemplifies the evolutionary genomic mechanisms underlying adoption of the pathogenic niche, and studies of its evolution have uncovered a rich array of discoveries about how new pathogens are made. Expected final online publication date for the Annual Review of Microbiology, Volume 76 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

A dominant clonal lineage of Streptococcus uberis in cattle in Germany

Bovine mastitis causes enormous economic losses in the dairy industry with Streptococcus uberis as one of the most common bacterial pathogens causing clinical and subclinical variations. In most cases mastitis can be cured by intramammary administration of antimicrobial agents. However, the severity of the clinical manifestations can vary greatly from mild to severe symtoms. In this study, a comparative genomic analysis of 24 S. uberis isolates from three dairy farms in Germany, affected by different courses of infection was conducted. While there were sporadic mild infections in farm A and B, a large number of infections were observed within a very short period of time in farm C. The comparison of virulence genes, antimicrobial resistance genes and prophage regions revealed no features that might be responsible for this severe course. However, almost all isolates from farm C showed the same, novel MLST profile (ST1373), thus a clonal outbreak cannot be excluded, whereby the actual reason for the particular virulence remains unknown. This study demonstrates the importance of extensive metagenomic studies, including the host genomes and the environment, to gain further evidence on the pathogenicity of S. uberis.

A Large-Scale Genome-Based Survey of Acidophilic Bacteria Suggests That Genome Streamlining Is an Adaption for Life at Low pH

The genome streamlining theory suggests that reduction of microbial genome size optimizes energy utilization in stressful environments. Although this hypothesis has been explored in several cases of low-nutrient (oligotrophic) and high-temperature environments, little work has been carried out on microorganisms from low-pH environments, and what has been reported is inconclusive. In this study, we performed a large-scale comparative genomics investigation of more than 260 bacterial high-quality genome sequences of acidophiles, together with genomes of their closest phylogenetic relatives that live at circum-neutral pH. A statistically supported correlation is reported between reduction of genome size and decreasing pH that we demonstrate is due to gene loss and reduced gene sizes. This trend is independent from other genome size constraints such as temperature and G + C content. Genome streamlining in the evolution of acidophilic bacteria is thus supported by our results. The analyses of predicted Clusters of Orthologous Genes (COG) categories and subcellular location predictions indicate that acidophiles have a lower representation of genes encoding extracellular proteins, signal transduction mechanisms, and proteins with unknown function but are enriched in inner membrane proteins, chaperones, basic metabolism, and core cellular functions. Contrary to other reports for genome streamlining, there was no significant change in paralog frequencies across pH. However, a detailed analysis of COG categories revealed a higher proportion of genes in acidophiles in the following categories: "replication and repair," "amino acid transport," and "intracellular trafficking". This study brings increasing clarity regarding the genomic adaptations of acidophiles to life at low pH while putting elements, such as the reduction of average gene size, under the spotlight of streamlining theory.

Adaptation to genome decay in the structure of the smallest eukaryotic ribosome

The evolution of microbial parasites involves the counterplay between natural selection forcing parasites to improve and genetic drifts forcing parasites to lose genes and accumulate deleterious mutations. Here, to understand how this counterplay occurs at the scale of individual macromolecules, we describe cryo-EM structure of ribosomes from Encephalitozoon cuniculi, a eukaryote with one of the smallest genomes in nature. The extreme rRNA reduction in E. cuniculi ribosomes is accompanied with unparalleled structural changes, such as the evolution of previously unknown molten rRNA linkers and bulgeless rRNA. Furthermore, E. cuniculi ribosomes withstand the loss of rRNA and protein segments by evolving an ability to use small molecules as structural mimics of degenerated rRNA and protein segments. Overall, we show that the molecular structures long viewed as reduced, degenerated, and suffering from debilitating mutations possess an array of compensatory mechanisms that allow them to remain active despite the extreme molecular reduction.

Many parasitic organisms contain molecular structures that are drastically smaller than analogous structures in non-parasitic organisms. Here the authors describe a cryo-EM structure of the ribosome from E. cuniculi that reveals that it compensated rRNA truncations by evolving the ability to use small molecules as ribosomal building blocks.

Genome reorganization during emergence of host-associated Mycobacterium abscessus

Mycobacterium abscessus is a rapid growing, free-living species of bacterium that also causes lung infections in humans. Human infections are usually acquired from the environment; however, dominant circulating clones (DCCs) have emerged recently in both M. abscessus subsp. massiliense and subsp. abscessus that appear to be transmitted among humans and are now globally distributed. These recently emerged clones are potentially informative about the ecological and evolutionary mechanisms of pathogen emergence and host adaptation. The geographical distribution of DCCs has been reported, but the genomic processes underlying their transition from environmental bacterium to human pathogen are not well characterized. To address this knowledge gap, we delineated the structure of M. abscessus subspecies abscessus and massiliense using genomic data from 200 clinical isolates of M. abscessus from seven geographical regions. We identified differences in overall patterns of lateral gene transfer (LGT) and barriers to LGT between subspecies and between environmental and host-adapted bacteria. We further characterized genome reorganization that accompanied bacterial host adaptation, inferring selection pressures acting at both genic and intergenic loci. We found that both subspecies encode an expansive pangenome with many genes at rare frequencies. Recombination appears more frequent in M. abscessus subsp. massiliense than in subsp. abscessus, consistent with prior reports. We found evidence suggesting that phage are exchanged between subspecies, despite genetic barriers evident elsewhere throughout the genome. Patterns of LGT differed according to niche, with less LGT observed among host-adapted DCCs versus environmental bacteria. We also found evidence suggesting that DCCs are under distinct selection pressures at both genic and intergenic sites. Our results indicate that host adaptation of M. abscessus was accompanied by major changes in genome evolution, including shifts in the apparent frequency of LGT and impacts of selection. Differences were evident among the DCCs as well, which varied in the degree of gene content remodelling, suggesting they were placed differently along the evolutionary trajectory toward host adaptation. These results provide insight into the evolutionary forces that reshape bacterial genomes as they emerge into the pathogenic niche.

Mutation rate dynamics reflect ecological change in an emerging zoonotic pathogen

Mutation rates vary both within and between bacterial species, and understanding what drives this variation is essential for understanding the evolutionary dynamics of bacterial populations. In this study, we investigate two factors that are predicted to influence the mutation rate: ecology and genome size. We conducted mutation accumulation experiments on eight strains of the emerging zoonotic pathogen Streptococcus suis. Natural variation within this species allows us to compare tonsil carriage and invasive disease isolates, from both more and less pathogenic populations, with a wide range of genome sizes. We find that invasive disease isolates have repeatedly evolved mutation rates that are higher than those of closely related carriage isolates, regardless of variation in genome size. Independent of this variation in overall rate, we also observe a stronger bias towards G/C to A/T mutations in isolates from more pathogenic populations, whose genomes tend to be smaller and more AT-rich. Our results suggest that ecology is a stronger correlate of mutation rate than genome size over these timescales, and that transitions to invasive disease are consistently accompanied by rapid increases in mutation rate. These results shed light on the impact that ecology can have on the adaptive potential of bacterial pathogens.

Large-scale genomic analysis of antimicrobial resistance in the zoonotic pathogen Streptococcus suis

BACKGROUND

Antimicrobial resistance (AMR) is among the gravest threats to human health and food security worldwide. The use of antimicrobials in livestock production can lead to emergence of AMR, which can have direct effects on humans through spread of zoonotic disease. Pigs pose a particular risk as they are a source of zoonotic diseases and receive more antimicrobials than most other livestock. Here we use a large-scale genomic approach to characterise AMR in Streptococcus suis, a commensal found in most pigs, but which can also cause serious disease in both pigs and humans.

RESULTS

We obtained replicated measures of Minimum Inhibitory Concentration (MIC) for 16 antibiotics, across a panel of 678 isolates, from the major pig-producing regions of the world. For several drugs, there was no natural separation into 'resistant' and 'susceptible', highlighting the need to treat MIC as a quantitative trait. We found differences in MICs between countries, consistent with their patterns of antimicrobial usage. AMR levels were high even for drugs not used to treat S. suis, with many multidrug-resistant isolates. Similar levels of resistance were found in pigs and humans from regions associated with zoonotic transmission. We next used whole genome sequences for each isolate to identify 43 candidate resistance determinants, 22 of which were novel in S. suis. The presence of these determinants explained most of the variation in MIC. But there were also interesting complications, including epistatic interactions, where known resistance alleles had no effect in some genetic backgrounds. Beta-lactam resistance involved many core genome variants of small effect, appearing in a characteristic order.

CONCLUSIONS

We present a large dataset allowing the analysis of the multiple contributing factors to AMR in S. suis. The high levels of AMR in S. suis that we observe are reflected by antibiotic usage patterns but our results confirm the potential for genomic data to aid in the fight against AMR.