progressiveMauve: multiple genome alignment with gene gain, loss and rearrangement

Integrative and conjugative elements of Pasteurella multocida: Prevalence and signatures in population evolution

Pasteurella multocida (P. multocida) is a bacterial pathogen responsible for a range of infections in humans and various animal hosts, causing significant economic losses in farming. Integrative and conjugative elements (ICEs) are important horizontal gene transfer elements, potentially enabling host bacteria to enhance adaptability by acquiring multiple functional genes. However, the understanding of ICEs in P. multocida and their impact on the transmission of this pathogen remains limited. In this study, 42 poultry-sourced P. multocida genomes obtained by high-throughput sequencing together with 393 publicly available P. multocida genomes were used to analyse the horizontal transfer of ICEs. Eighty-two ICEs were identified in P. multocida, including SXT/R391 and Tn916 subtypes, as well as three subtypes of ICEHin1056 family, with the latter being widely prevalent in P. multocida and carrying multiple resistance genes. The correlations between insertion sequences and resistant genes in ICEs were also identified, and some ICEs introduced the carbapenem gene blaOXA-2 and the bleomycin gene bleO to P. multocida. Phylogenetic and collinearity analyses of these bioinformatics found that ICEs in P. multocida were transmitted vertically and horizontally and have evolved with host specialization. These findings provide insight into the transmission and evolution mode of ICEs in P. multocida and highlight the importance of understanding these elements for controlling the spread of antibiotic resistance.

Identification of highly conserved Trypanosoma cruzi antigens for the development of a universal serological diagnostic assay

Chagas Disease is an important neglected tropical disease caused by Trypanosoma cruzi. There is no gold standard for diagnosis and commercial serological tests perform poorly in certain locations. By aligning T. cruzi genomes covering parasite genetic and geographic diversity, we identified highly conserved proteins that could serve as universal antigens for improved diagnosis. Their antigenicity was tested in high-density peptide microarrays using well-characterized plasma samples, including samples presenting true infections but discordant serology. Individual and combination of epitopes were also evaluated in peptide-ELISAs. We identified >1400 highly conserved T. cruzi proteins evaluated in microarrays. Remarkably, T. cruzi positive controls had a different epitope recognition profile compared to serologically discordant samples. In particular, multiple T. cruzi antigens used in current tests and their strain-variants, and novel epitopes thought to be broadly antigenic failed to be recognized by discordant samples. Nonetheless, >2000 epitopes specifically recognized by IgGs from both positive controls and discordant samples were identified. Evaluation of selected peptides in ELISA further illustrated the extensive variation in antibody profiles among subjects and a peptide combination could outperform a commercial ELISA, increasing assay sensitivity from 52.3% to 72.7%. Individual variation in antibody profiles rather than T. cruzi diversity appears to be the main factor driving differences in serological diagnostic performance according to geography, which will be important to further elucidate. ELISA with a combination of peptides recognized by a greater number of individuals could better capture infections, and further development may lead to an optimal antigen mixture for a universal diagnostic assay.

Comparative analysis of mitochondrial genomes in Ceratocystis fimbriata complex across diverse hosts

The decline ofAcacia mangiumWilld. in Malaysia, especially in Sabah since 2010, is primarily due to Ceratocystiswilt and canker disease (CWCD) caused by theCeratocystis fimbriataEllis & Halst. complex. This study was aimed to investigate the mitochondrial genome architecture of two differentC. fimbriatacomplex isolates from Malaysia: one fromA. mangiumin Pahang (FRIM1162) and another fromEucalyptus pellitain Sarawak (FRIM1441). This research employed Next-Generation Sequencing (NGS) to contrast genomes from diverse hosts with nine additional mitochondrial sequences, identifying significant genetic diversity and mutational hotspots in the mitochondrial genome alignment. The mitochondrial genome-based phylogenetic analysis revealed a significant genetic relationship between the studied isolates and theC. fimbriatacomplex in the South American Subclade, indicating that theC. fimbriatacomplex discovered in Malaysia isC. manginecans. The comparative mitochondrial genome demonstrates the adaptability of the complex due to mobile genetic components and genomic rearrangements in the studiedfungal isolates. This research enhances our knowledge of the genetic diversity and evolutionary patterns within theC. fimbriatacomplex, aiding in a deeper understanding of fungal disease development and host adaption processes. The acquired insights are crucial for creating specific management strategies for CWCD, improving the overall understanding of fungal disease evolution and control.

Clonal background and routes of plasmid transmission underlie antimicrobial resistance features of bloodstream Klebsiella pneumoniae

Bloodstream infections caused by the opportunistic pathogen Klebsiella pneumoniae are associated with adverse health complications and high mortality rates. Antimicrobial resistance (AMR) limits available treatment options, thus exacerbating its public health and clinical burden. Here, we aim to elucidate the population structure of K. pneumoniae in bloodstream infections from a single medical center and the drivers that facilitate the dissemination of AMR. Analysis of 136 short-read genome sequences complemented with 12 long-read sequences shows the population consisting of 94 sequence types (STs) and 99 clonal groups, including globally distributed multidrug resistant and hypervirulent clones. In vitro antimicrobial susceptibility testing and in silico identification of AMR determinants reveal high concordance (90.44-100%) for aminoglycosides, beta-lactams, carbapenems, cephalosporins, quinolones, and sulfonamides. IncF plasmids mediate the clonal (within the same lineage) and horizontal (between lineages) transmission of the extended-spectrum beta-lactamase gene blaCTX-M-15. Nearly identical plasmids are recovered from isolates over a span of two years indicating long-term persistence. The genetic determinants for hypervirulence are carried on plasmids exhibiting genomic rearrangement, loss, and/or truncation. Our findings highlight the importance of considering both the genetic background of host strains and the routes of plasmid transmission in understanding the spread of AMR in bloodstream infections.

An enterococcal phage protein inhibits type IV restriction enzymes involved in antiphage defense

The prevalence of multidrug resistant (MDR) bacterial infections continues to rise as the development of antibiotics needed to combat these infections remains stagnant. MDR enterococci are a major contributor to this crisis. A potential therapeutic approach for combating MDR enterococci is bacteriophage (phage) therapy, which uses lytic viruses to infect and kill pathogenic bacteria. While phages that lyse some strains of MDR enterococci have been identified, other strains display high levels of resistance and the mechanisms underlying this resistance are poorly defined. Here, we use a CRISPR interference (CRISPRi) screen to identify a genetic locus found on a mobilizable plasmid from Enterococcus faecalis involved in phage resistance. This locus encodes a putative serine recombinase followed by a Type IV restriction enzyme (TIV-RE) that we show restricts the replication of phage phi47 in vancomycin-resistant E. faecalis. We further find that phi47 evolves to overcome restriction by acquiring a missense mutation in a TIV-RE inhibitor protein. We show that this inhibitor, termed type IV restriction inhibiting factor A (tifA), binds and inactivates diverse TIV-REs. Overall, our findings advance our understanding of phage defense in drug-resistant E. faecalis and provide mechanistic insight into how phages evolve to overcome antiphage defense systems.

Plastome Evolution, Phylogenomics, and DNA Barcoding Investigation of Gastrochilus (Aeridinae, Orchidaceae), with a Focus on the Systematic Position of Haraella retrocalla

Gastrochilus is an orchid genus containing about 70 species in tropical and subtropical Asia with high morphological diversity. The phylogenetic relationships among this genus have not been fully resolved, and the plastome evolution has not been investigated either. In this study, five plastomes of Gastrochilus were newly reported, and sixteen plastomes of Gastrochilus were used to conduct comparative and phylogenetic analyses. Our results showed that the Gastrochilus plastomes ranged from 146,183 to 148,666 bp, with a GC content of 36.7-36.9%. There were 120 genes annotated, consisting of 74 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. No contraction and expansion of IR borders, gene rearrangements, or inversions were detected. Additionally, the repeat sequences and codon usage bias of Gastrochilus plastomes were highly conserved. Twenty hypervariable regions were selected as potential DNA barcodes. The phylogenetic relationships within Gastrochilus were well resolved based on the whole plastome, especially among main clades. Furthermore, both molecular and morphological data strongly supported Haraella retrocalla as a member of Gastrochilus (G. retrocallus).

An enterococcal phage protein broadly inhibits type IV restriction enzymes involved in antiphage defense

The prevalence of multidrug resistant (MDR) bacterial infections continues to rise as the development of antibiotics needed to combat these infections remains stagnant. MDR enterococci are a major contributor to this crisis. A potential therapeutic approach for combating MDR enterococci is bacteriophage (phage) therapy, which uses lytic viruses to infect and kill pathogenic bacteria. While phages that lyse some strains of MDR enterococci have been identified, other strains display high levels of resistance and the mechanisms underlying this resistance are poorly defined. Here, we use a CRISPR interference (CRISPRi) screen to identify a genetic locus found on a mobilizable plasmid from Enterococcus faecalis involved in phage resistance. This locus encodes a putative serine recombinase followed by a Type IV restriction enzyme (TIV-RE) that we show restricts the replication of phage phi47 in E. faecalis. We further find that phi47 evolves to overcome restriction by acquiring a missense mutation in a TIV-RE inhibitor protein. We show that this inhibitor, termed type IV restriction inhibiting factor A (tifA), binds and inactivates diverse TIV-REs. Overall, our findings advance our understanding of phage defense in drug-resistant E. faecalis and provide mechanistic insight into how phages evolve to overcome antiphage defense systems.

Discovery of a novel spotted fever group Rickettsia, "Candidatus Rickettsia kedanie," in unfed larval chigger mites, Leptotrombidium scutellare

Spotted fever group (SFG) rickettsia, the causative agent of SFG rickettsiosis, is predominantly carried by ticks, whereas Orientia tsutusgamushi, the causative agent of scrub typhus, is primarily transmitted by chigger mites in Japan. In this study, we attempted to isolate intracellular eubacteria from Leptotrombidium scutellare, a major vector of O. tsutsugamushi; moreover, we isolated an SFG rickettsia using a mosquito-derived cell line. Draft genome sequences of this unique isolate, by applying criteria for species delimitation, classified this isolate as a novel strain, proposed as "Rickettsia kedanie." Further genetic analysis identified conserved virulence factors, and the isolate successfully propagated in mammalian cells, suggesting its ability to cause diseases in humans. The presence of SFG rickettsia in unfed larvae implies potential dual-pathogen carriage and reflects a symbiotic relationship similar to that between the mites and O. tsutsugamushi, indicating possibility of its transovarial transmission from female adults. Furthermore, conserved genomic similarity of the novel isolate to known SFG rickettsia suggests potential multiple hosts, including chiggers and ticks. In the natural environment, ticks, chigger mites, and wild animals may carry new isolates, complicating the infection cycle and increasing the transmission risks to humans. This discovery challenges the conventional association of SFG rickettsia with ticks, emphasizing its implications for research and disease control. However, this study was confined to a particular species of chigger mites and geographic area, underscoring the necessity for additional studies to comprehend the ecological dynamics, host interactions, and health implications linked to this newly identified SFG rickettsia.

Mitochondrial genome of the fluke pond snail, Austropeplea cf. brazieri (Gastropoda: Lymnaeidae)

BACKGROUND

Lymnaeid snails of the genus Austropeplea are an important vector of the liver fluke (Fasciola hepatica), contributing to livestock production losses in Australia and New Zealand. However, the species status within Austropeplea is ambiguous due to heavy reliance on morphological analysis and a relative lack of genetic data. This study aimed to characterise the mitochondrial genome of A. cf. brazieri, an intermediate host of liver fluke in eastern Victoria.

METHODS

The mitochondrial genome was assembled and annotated from a combination of second- and third-generation sequencing data. For comparative purposes, we performed phylogenetic analyses of the concatenated nucleotide sequences of the mitochondrial protein-coding genes, cytochrome c oxidase subunit 1 and 16S genes.

RESULTS

The assembled mt genome was 13,757 base pairs and comprised 37 genes, including 13 protein-coding genes, 22 transfer RNA genes and 2 ribosomal RNA genes. The mt genome length, gene order and nucleotide compositions were similar to related species of lymnaeids. Phylogenetic analyses of the mt nucleotide sequences placed A. cf. brazieri within the same clade as Orientogalba ollula with strong statistical supports. Phylogenies of the cox1 and 16S mt sequences were constructed due to the wide availability of these sequences representing the lymnaeid taxa. As expected in both these phylogenies, A. cf. brazieri clustered with other Austropeplea sequences, but the nodal supports were low.

CONCLUSIONS

The representative mt genome of A. cf. brazieri should provide a useful resource for future molecular, epidemiology and parasitological studies of this socio-economically important lymnaeid species.

De Novo Hybrid Assembly Unveils Multi-Chromosomal Mitochondrial Genomes in Ludwigia Species, Highlighting Genomic Recombination, Gene Transfer, and RNA Editing Events

Biological invasions have been identified as the fifth cause of biodiversity loss, and their subsequent dispersal represents a major ecological challenge. The aquatic invasive species Ludwigia grandiflora subsp. hexapetala (Lgh) and Ludwigia peploides subsp. montevidensis (Lpm) are largely distributed in aquatic environments in North America and in Europe. However, they also present worrying terrestrial forms that are able to colonize wet meadows. To comprehend the mechanisms of the terrestrial adaptation of Lgh and Lpm, it is necessary to develop their genomic resources, which are currently poorly documented. We performed de novo assembly of the mitogenomes of Lgh and Lpm through hybrid assemblies, combining short reads (SR) and/or long reads (LR) before annotating both mitogenomes. We successfully assembled the mitogenomes of Lgh and Lpm into two circular molecules each, resulting in a combined total length of 711,578 bp and 722,518 bp, respectively. Notably, both the Lgh and Lpm molecules contained plastome-origin sequences, comprising 7.8% of the mitochondrial genome length. Additionally, we identified recombinations that were mediated by large repeats, suggesting the presence of multiple alternative conformations. In conclusion, our study presents the first high-quality mitogenomes of Lpm and Lgh, which are the only ones in the Myrtales order found as two circular molecules.

Genomic analysis of a novel phage vB_SenS_ST1UNAM with lytic activity against Salmonella enterica serotypes

In this study, we present the complete annotated genome of a novel Salmonella phage, vB_SenS_ST1UNAM. This phage exhibits lytic activity against several Salmonella enterica serotypes, such as S. Typhi, S. Enteritidis, and S. Typhimurium strains, which are major causes of foodborne illness worldwide. Its genome consists of a linear, double-stranded DNA of 47,877 bp with an average G+C content of 46.6%. A total of 85 coding regions (CDS) were predicted, of which only 43 CDS were functionally assigned. Neither genes involved in the regulation of lysogeny, nor antibiotic resistance genes were identified. This phage harbors a lytic cassette that encodes a type II-holin and a Rz/Rz1-like spanin complex, along with a restriction-modification evasion system and a depolymerase that degrades Salmonella exopolysaccharide. Moreover, the comparative analysis with closely related phage genomes revealed that vB_SenS_ST1UNAM represents a novel genus, for which the genus "Gomezvirus" within the subfamily "ST1UNAM-like" is proposed.

Time-calibrated phylogenetic and chromosomal mobilome analyses of Staphylococcus aureus CC398 reveal geographical and host-related evolution

An international collection of Staphylococcus aureus of clonal complex (CC) 398 from diverse hosts spanning all continents and a 30 year-period is studied based on whole-genome sequencing (WGS) data. The collection consists of publicly available genomic data from 2994 strains and 134 recently sequenced Swiss methicillin-resistant S. aureus (MRSA) CC398 strains. A time-calibrated phylogeny reveals the presence of distinct phylogroups present in Asia, North and South America and Europe. European MRSA diverged from methicillin-susceptible S. aureus (MSSA) at the beginning of the 1950s. Two major European phylogroups (EP4 and EP5), which diverged approximately 1974, are the main drivers of MRSA CC398 spread in Europe. Within EP5, an emergent MRSA lineage spreading among the European horse population (EP5-Leq) diverged approximately 1996 from the pig lineage (EP5-Lpg), and also contains human-related strains. EP5-Leq is characterized by staphylococcal cassette chromosome mec (SCCmec) IVa and spa type t011 (CC398-IVa-t011), and EP5-Lpg by CC398-SCCmecVc-t011. The lineage-specific antibiotic resistance and virulence gene patterns are mostly mediated by the acquisition of mobile genetic elements like SCCmec, S. aureus Genomic Islands (SaGIs), prophages and transposons. Different combinations of virulence factors are present on S. aureus pathogenicity islands (SaPIs), and novel antimicrobial resistance gene containing elements are associated with certain lineages expanding in Europe. This WGS-based analysis reveals the actual evolutionary trajectory and epidemiological trend of the international MRSA CC398 population considering host, temporal, geographical and molecular factors. It provides a baseline for global WGS-based One-Health studies of adaptive evolution of MRSA CC398 as well as for local outbreak investigations.

Complete chloroplast genomes of the hemiparasitic genus Cymbaria: Insights into comparative analysis, development of molecular markers, and phylogenetic relationships

The hemiparasitic tribe Cymbarieae (Orobanchaceae) plays a crucial role in elucidating the initial stage of the transition from autotrophism to heterotrophism. However, the complete chloroplast genome of the type genus Cymbaria has yet to be reported. In addition, the traditional Mongolian medicine Cymbaria daurica is frequently subjected to adulteration or substitution because of the minor morphological differences with Cymbaria mongolica. In this study, the complete chloroplast genomes of the two Cymbaria species were assembled and annotated, and those of other published 52 Orobanchaceae species were retrieved for comparative analyses. We found that the Cymbaria chloroplast genomes are characterized by pseudogenization or loss of stress-relevant genes (ndh) and a unique rbcL-matK inversion. Unlike the high variability observed in holoparasites, Cymbaria and other hemiparasites exhibit high similarity to autotrophs in genome size, guanine-cytosine (GC) content, and intact genes. Notably, four pairs of specific DNA barcodes were developed and validated to distinguish the medicinal herb from its adulterants. Phylogenetic analyses revealed that the genus Cymbaria and the Schwalbea-Siphonostegia clade are grouped into the tribe Cymbarieae, which forms a sister clade to the remaining Orobanchaceae parasitic lineages. Moreover, the diversification of monophyletic Cymbaria occurred during the late Miocene (6.72 Mya) in the Mongol-Chinese steppe region. Our findings provide valuable genetic resources for studying the phylogeny of Orobanchaceae and plant parasitism, and genetic tools to validate the authenticity of the traditional Mongolian medicine "Xinba.".

Plastid genome evolution in leafless members of the orchid subfamily Orchidoideae, with a focus on Degranvillea dermaptera

PREMISE

Leafless, heterotrophic plants are prime examples of organismal modification, the genomic consequences of which have received considerable interest. In particular, plastid genomes (plastomes) are being sequenced at a high rate, allowing continual refinement of conceptual models of reductive evolution in heterotrophs. However, numerous sampling gaps exist, hindering the ability to conduct comprehensive phylogenomic analyses in these plants.

METHODS

Using floral tissue from an herbarium specimen, we sequenced and analyzed the plastome of Degranvillea dermaptera, a rarely collected, leafless orchid species from South America about which little is known, including its phylogenetic affinities.

RESULTS

The plastome is the most reduced of those sequenced among the orchid subfamily Orchidoideae. In Degranvillea, it has lost the majority of genes found in leafy autotrophic species, is structurally rearranged, and has similar gene content to the most reduced plastomes among the orchids. We found strong evidence for the placement of Degranvillea within the subtribe Spiranthinae using models that explicitly account for heterotachy, or lineage-specific evolutionary rate variation over time. We further found evidence of relaxed selection on several genes and of correlations among substitution rates and several other "traits" of the plastome among leafless members of orchid subfamily Orchidoideae.

CONCLUSIONS

Our findings advance knowledge on the phylogenetic relationships and paths of plastid genome evolution among the orchids, which have experienced more independent transitions to heterotrophy than any other plant family. This study demonstrates the importance of herbarium collections in comparative genomics of poorly known species of conservation concern.

Development and validation of a multi-target TaqMan qPCR method for detection of Borrelia burgdorferi sensu lato

Reliable detection of bacteria belonging to the Borrelia burgdorferi sensu lato species complex in vertebrate reservoirs, tick vectors, and patients is key to answer questions regarding Lyme borreliosis epidemiology. Nevertheless, the description of characteristics of qPCRs for the detection of B. burgdorferi s. l. are often limited. This study covers the development and validation of two duplex taqman qPCR assays used to target four markers on the chromosome of genospecies of B. burgdorferi s. l. Analytical specificity was determined with a panel of spirochete strains. qPCR characteristics were specified using water or tick DNA spiked with controlled quantities of the targeted DNA sequences of B. afzelii, B. burgdorferi sensu stricto or B. bavariensis. The effectiveness of detection results was finally evaluated using DNA extracted from ticks and biopsies from mammals whose infectious status had been determined by other detection assays. The developed qPCR assays allow exclusive detection of B. burgdorferi s. l. with the exception of the M16 marker which also detect relapsing fever Borreliae. The limit of detection is between 10 and 40 copies per qPCR reaction depending on the sample type, the B. burgdorferi genospecies and the targeted marker. Detection tests performed on various kind of samples illustrated the accuracy and robustness of our qPCR assays. Within the defined limits, this multi-target qPCR method allows a versatile detection of B. burgdorferi s. l., regardless of the genospecies and the sample material analyzed, with a sensitivity that would be compatible with most applications and a reproducibility of 100% under measurement conditions of limits of detection, thereby limiting result ambiguities.

A phosphorylation signal activates genome-wide transcriptional control by BfmR, the global regulator of Acinetobacter resistance and virulence

The nosocomial pathogen Acinetobacter baumannii is a major threat to human health. The sensor kinase-response regulator system, BfmS-BfmR, is essential to multidrug resistance and virulence in the bacterium and represents a potential antimicrobial target. Important questions remain about how the system controls resistance and pathogenesis. Although BfmR knockout alters expression of >1000 genes, its direct regulon is undefined. Moreover, how phosphorylation controls the regulator is unclear. Here, we address these problems by combining mutagenesis, ChIP-seq, and in vitro phosphorylation to study the functions of phospho-BfmR. We show that phosphorylation is required for BfmR-mediated gene regulation, antibiotic resistance, and sepsis development in vivo. Consistent with activating the protein, phosphorylation induces dimerization and target DNA affinity. Integrated analysis of genome-wide binding and transcriptional profiles of BfmR led to additional key findings: (1) Phosphorylation dramatically expands the number of genomic sites BfmR binds; (2) DNA recognition involves a direct repeat motif widespread across promoters; (3) BfmR directly regulates 303 genes as activator (eg, capsule, peptidoglycan, and outer membrane biogenesis) or repressor (pilus biogenesis); (4) BfmR controls several non-coding sRNAs. These studies reveal the centrality of a phosphorylation signal in driving A. baumannii disease and disentangle the extensive pathogenic gene-regulatory network under its control.

Evidence of Lactobacillus strains shared between the female urinary and vaginal microbiota

Lactobacillus species are common inhabitants of the 'healthy' female urinary and vaginal communities, often associated with a lack of symptoms in both anatomical sites. Given identification by prior studies of similar bacterial species in both communities, it has been hypothesized that the two microbiotas are in fact connected. Here, we carried out whole-genome sequencing of 49 Lactobacillus strains, including 16 paired urogenital samples from the same participant. These strains represent five different Lactobacillus species: L. crispatus, L. gasseri, L. iners, L. jensenii, and L. paragasseri. Average nucleotide identity (ANI), alignment, single-nucleotide polymorphism (SNP), and CRISPR comparisons between strains from the same participant were performed. We conducted simulations of genome assemblies and ANI comparisons and present a statistical method to distinguish between unrelated, related, and identical strains. We found that 50 % of the paired samples have identical strains, evidence that the urinary and vaginal communities are connected. Additionally, we found evidence of strains sharing a common ancestor. These results establish that microbial sharing between the urinary tract and vagina is not limited to uropathogens. Knowledge that these two anatomical sites can share lactobacilli in females can inform future clinical approaches.

Assembly and Genome Annotation of Different Strains of Apple Fruit Moth Virus (Cydia pomonella granulovirus)

Cydia pomonella granulovirus is a natural pathogen for Cydia pomonella that is used as a biocontrol agent of insect populations. The study of granulovirus virulence is of particular interest since the development of resistance in natural populations of C. pomonella has been observed during the long-term use of the Mexican isolate CpGV. In our study, we present the genomes of 18 CpGV strains endemic to southern Russia and from Kazakhstan, as well as a strain included in the commercial preparation "Madex Twin", which were sequenced and analyzed. We performed comparative genomic analysis using several tools. From comparisons at the level of genes and protein products that are involved in the infection process of virosis, synonymous and missense substitution variants have been identified. The average nucleotide identity has demonstrated a high similarity with other granulovirus genomes of different geographic origins. Whole-genome alignment of the 18 genomes relative to the reference revealed regions of low similarity. Analysis of gene repertoire variation has shown that BZR GV 4, BZR GV 6, and BZR GV L-7 strains have been the closest in gene content to the commercial "Madex Twin" strain. We have confirmed two deletions using read depth coverage data in regions lacking genes shown by homology analysis for granuloviruses BZR GV L-4 and BZR GV L-6; however, they are not related to the known genes causing viral pathogenicity. Thus, we have isolated novel CpGV strains and analyzed their potential as strains producing highly effective bioinsecticides against C. pomonella.

In Silico Analysis of the Phylogenetic and Physiological Characteristics of Sphingobium indicum B90A: A Hexachlorocyclohexane-Degrading Bacterium

The study focuses on the in silico genomic characterization of Sphingobium indicum B90A, revealing a wealth of genes involved in stress response, carbon monoxide oxidation, β-carotene biosynthesis, heavy metal resistance, and aromatic compound degradation, suggesting its potential as a bioremediation agent. Furthermore, genomic adaptations among nine Sphingomonad strains were explored, highlighting shared core genes via pangenome analysis, including those related to the shikimate pathway and heavy metal resistance. The majority of genes associated with aromatic compound degradation, heavy metal resistance, and stress response were found within genomic islands across all strains. Sphingobium indicum UT26S exhibited the highest number of genomic islands, while Sphingopyxis alaskensis RB2256 had the maximum fraction of its genome covered by genomic islands. The distribution of lin genes varied among the strains, indicating diverse genetic responses to environmental pressures. Additionally, in silico evidence of horizontal gene transfer (HGT) between plasmids pSRL3 and pISP3 of the Sphingobium and Sphingomonas genera, respectively, has been provided. The manuscript offers novel insights into strain B90A, highlighting its role in horizontal gene transfer and refining evolutionary relationships among Sphingomonad strains. The discovery of stress response genes and the czcABCD operon emphasizes the potential of Sphingomonads in consortia development, supported by genomic island analysis.

Comparative chloroplast genomes study of five officinal Ardisia Species: Unraveling interspecific diversity and evolutionary insights in Ardisia

Ardisia S.W. (Primulaceae), naturally distributed in tropical and subtropical regions, has edible and medicinal values and is prevalent in clinical and daily use in China. More genetic information for distinct species delineation is needed to support the development and utilization of the genus Ardisia. We sequenced, annotated, and compared the chloroplast genomes of five Ardisia species: A. brunnescens, A. pusilla, A. squamulosa, A. crenata, and A. brevicaulis in this study. We found a typical quadripartite structure in all five chloroplast genomes, with lengths ranging from 155,045 to 156,943 bp. Except for A. pusilla, which lacked the ycf15 gene, the other four Ardisia species contained 114 unique genes, including 79 protein-coding genes, 30 tRNAs, and four rRNAs. In addition, the rps19 pseudogene gene was present only in A. brunnescens. Five highly variable DNA barcodes were identified for five Ardisia species, including trnT-GGU-psbD, trnT-UGU-trnL-UAA, rps4-trnT-UGU, rpl32-trnL-UAG, and rpoB-trnC-GAA. The RNA editiing sites of protein-coding genes in the five Ardisia plastome were characterized and compared, and 274 (A. crenata)-288 (A. brevicaulis) were found. The results of the phylogenetic analysis were consistent with the morphological classification. Sequence alignment and phylogenetic analysis showed that ycf15 genes were highly divergent in Primulaceae. Reconstructions of ancestral character states indicated that leaf margin morphology is critical for classifying the genus Ardisia, with a rodent-like character being the most primitive. These results provide valuable information on the taxonomy and evolution of Ardisia plants.

Strain tracking in complex microbiomes using synteny analysis reveals per-species modes of evolution

Microbial species diversify into strains through single-nucleotide mutations and structural changes, such as recombination, insertions and deletions. Most strain-comparison methods quantify differences in single-nucleotide polymorphisms (SNPs) and are insensitive to structural changes. However, recombination is an important driver of phenotypic diversification in many species, including human pathogens. We introduce SynTracker, a tool that compares microbial strains using genome synteny-the order of sequence blocks in homologous genomic regions-in pairs of metagenomic assemblies or genomes. Genome synteny is a rich source of genomic information untapped by current strain-comparison tools. SynTracker has low sensitivity to SNPs, has no database requirement and is robust to sequencing errors. It outperforms existing tools when tracking strains in metagenomic data and is particularly suited for phages, plasmids and other low-data contexts. Applied to single-species datasets and human gut metagenomes, SynTracker, combined with an SNP-based tool, detects strains enriched in either point mutations or structural changes, providing insights into microbial evolution in situ.

Completed genome and emergence scenario of the multidrug-resistant nosocomial pathogen Staphylococcus epidermidis ST215

BACKGROUND

A multidrug-resistant lineage of Staphylococcus epidermidis named ST215 is a common cause of prosthetic joint infections and other deep surgical site infections in Northern Europe, but is not present elsewhere. The increasing resistance among S. epidermidis strains is a global concern. We used whole-genome sequencing to characterize ST215 from healthcare settings.

RESULTS

We completed the genome of a ST215 isolate from a Swedish hospital using short and long reads, resulting in a circular 2,676,787 bp chromosome and a 2,326 bp plasmid. The new ST215 genome was placed in phylogenetic context using 1,361 finished public S. epidermidis reference genomes. We generated 10 additional short-read ST215 genomes and 11 short-read genomes of ST2, which is another common multidrug-resistant lineage at the same hospital. We studied recombination's role in the evolution of ST2 and ST215, and found multiple recombination events averaging 30-50 kb. By comparing the results of antimicrobial susceptibility testing for 31 antimicrobial drugs with the genome content encoding antimicrobial resistance in the ST215 and ST2 isolates, we found highly similar resistance traits between the isolates, with 22 resistance genes being shared between all the ST215 and ST2 genomes. The ST215 genome contained 29 genes that were historically identified as virulence genes of S. epidermidis ST2. We established that in the nucleotide sequence stretches identified as recombination events, virulence genes were overrepresented in ST215, while antibiotic resistance genes were overrepresented in ST2.

CONCLUSIONS

This study features the extensive antibiotic resistance and virulence gene content in ST215 genomes. ST215 and ST2 lineages have similarly evolved, acquiring resistance and virulence through genomic recombination. The results highlight the threat of new multidrug-resistant S. epidermidis lineages emerging in healthcare settings.

Experimental evolution of Staphylococcus aureus in macrophages: dissection of a conditional adaptive trait promoting intracellular survival

Staphylococcus aureus is a major pathogen associated with important diseases in humans and animals. Macrophages are a key component of the innate immune response to S. aureus infection and play a major role in disease outcomes. To investigate the adaptive evolution of S. aureus in response to macrophages, we developed an experimental infection assay. S. aureus strains representing major human epidemic clones were passaged many times in a macrophage cell line, accumulating mutations in an array of genomic loci. Phenotypic analysis revealed the emergence of a lineage exhibiting increased survival in macrophages and human blood, and resistance to vancomycin. The evolved lineage exhibited a previously undescribed small colony variant (SCV) phenotype characterized by hyper-pigmentation, which resulted from a missense mutation in rsbW. Notably, the novel SCV was a conditional adaptive trait that was unstable in nutrient-replete conditions in vitro, rapidly converting from hyper-pigmented SCV to a non-pigmented large colony variant via spontaneous sigB deletion events. Importantly, we identified similar deletions in the genome sequences of a limited number of clinical S. aureus isolates from public databases, indicating that related events may occur during clinical infection. Experimental infection of zebrafish did not reveal a difference in virulence between parent and novel SCV but demonstrated an in vivo fitness cost for the compensatory sigB deletion events. Taken together, we report an experimental evolutionary approach for investigating bacterial innate immune cell interactions, revealing a conditional adaptation that promotes S. aureus survival in macrophages and resistance to vancomycin.

IMPORTANCE

Staphylococcus aureus is an important human bacterial pathogen. The host response to S. aureus involves the production of innate immune cells such as macrophages which are important for fighting infection. Here we report a new model of experimental evolution for studying how S. aureus can evade killing by macrophages. We identified a novel adaptive phenotype that promotes survival in macrophages and blood and resistance to antibiotics. The phenotype is lost rapidly upon growth in nutrient-rich conditions via disruption of the alternative sigma factor sigB, revealing a conditional niche-specific fitness advantage. Genomic analysis of clinical isolates suggests similar adaptations may occur during human infections. Our model may be used broadly to identify adaptations of S. aureus to the innate immune response.

Comparative genomics of Giardia duodenalis sub-assemblage AI beaver (Be-2) and human (WB-C6) strains show remarkable homozygosity, sequence similarity, and conservation of VSP genes

Giardia duodenalis, a major cause of waterborne infection, infects a wide range of mammalian hosts and is subdivided into eight genetically well-defined assemblages named A through H. However, fragmented genomes and a lack of comparative analysis within and between the assemblages render unclear the molecular mechanisms controlling host specificity and differential disease outcomes. To address this, we generated a near-complete de novo genome of AI assemblage using the Oxford Nanopore platform by sequencing the Be-2 genome. We generated 148,144 long-reads with quality scores of > 7. The final genome assembly consists of only nine contigs with an N50 of 3,045,186 bp. This assembly agrees closely with the assembly of another strain in the AI assemblage (WB-C6). However, a critical difference is that a region previously placed in the five-prime region of Chr5 belongs to Chr4 of Be-2. We find a high degree of conservation in the ploidy, homozygosity, and the presence of cysteine-rich variant-specific surface proteins (VSPs) within the AI assemblage. Our assembly provides a nearly complete genome of a member of the AI assemblage of G. duodenalis, aiding population genomic studies capable of elucidating Giardia transmission, host range, and pathogenicity.

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.

Characterization of the family-level Borreliaceae pan-genome and development of an episomal typing protocol

Background

The Borreliaceae family includes many obligate parasitic bacterial species which are etiologically associated with a myriad of zoonotic borrelioses including Lyme disease and vector-borne relapsing fevers. Infections by the Borreliaceae are difficult to detect by both direct and indirect methods, often leading to delayed and missed diagnoses. Efforts to improve diagnoses center around the development of molecular diagnostics (MDx), but due to deep tissue sequestration of the causative spirochaetes and the lack of persistent bacteremias, even MDx assays suffer from a lack of sensitivity. Additionally, the highly extensive genomic heterogeneity among isolates, even within the same species, contributes to the lack of assay sensitivity as single target assays cannot provide universal coverage. This within-species heterogeneity is partly due to differences in replicon repertoires and genomic structures that have likely arisen to support the complex Borreliaceae lifecycle in which these parasites have to survive in multiple hosts each with unique immune responses.

Results

We constructed a Borreliaceae family-level pangenome and characterized the phylogenetic relationships among the constituent taxa which supports the recent taxonomy of splitting the family into at least two genera. Gene content pro les were created for the majority of the Borreliaceae replicons, providing for the first time their unambiguous molecular typing.

Conclusion

Our characterization of the Borreliaceae pan-genome supports the splitting of the former Borrelia genus into two genera and provides for the phylogenetic placement of several non-species designated isolates. Mining this family-level pangenome will enable precision diagnostics corresponding to gene content-driven clinical outcomes while also providing targets for interventions.

GenoFig: a user-friendly application for the visualization and comparison of genomic regions

MOTIVATION

Understanding the molecular evolutionary history of organisms usually requires visual comparison of genomic regions from related species or strains. Although several applications already exist to achieve this task, they are either too old, too limited, or too complex for most user's needs.

RESULTS

GenoFig is a graphical application for the visualization of prokaryotic genomic regions, intended to be as easy to use as possible and flexible enough to adapt to a variety of needs. GenoFig allows the personalized representation of annotations extracted from GenBank files in a consistent way across sequences, using regular expressions. It also provides several unique options to optimize the display of homologous regions between sequences, as well as other more classical features such as sequence GC percent or GC-skew representations. In summary, GenoFig is a simple, free, and highly configurable tool to explore the evolution of specific genomic regions in prokaryotes and to produce publication-ready figures.

AVAILABILITY AND IMPLEMENTATION

Genofig is fully available at https://forgemia.inra.fr/public-pgba/genofig under a GPL 3.0 license.

A taxonomic proposal for cedratviruses, orpheoviruses, and pithoviruses

Orpheoviruses, cedratviruses, and pithoviruses are large DNA viruses that cluster together taxonomically within the order Pimascovirales of the phylum Nucleocytoviricota. However, they were not classified previously by the International Committee on Taxonomy of Viruses (ICTV). Here, we present a comprehensive analysis of the gene content, morphology, and phylogenomics of these viruses, providing data that underpinned the recent proposal to establish new taxa for their initial classification. The new taxonomy, which has now been ratified by the ICTV, includes the family Orpheoviridae and genus Alphaorpheovirus, the family Pithoviridae and genus Alphapithovirus, and the family Cedratviridae and genus Alphacedratvirus, aiming to formally catalogue the isolates covered in this study. Additionally, as per the newly adopted rules, we applied standardized binomial names for the virus species created to classify isolates with complete genome sequences available in public databases at the time of the proposal. The specific epithet of each virus species was chosen as a reference to the location where the exemplar virus was isolated.

Isolation of Saccharomycopsis species from plant material

Saccharomycopsis species are natural organic sulphur auxotrophs. Their genomes do not encode genes for the uptake and assimilation of sulphate and thus these species cannot grow on media lacking e.g. methionine. Due to the similarity between sulphate and selenate, uptake and assimilation of selenate occurs through the same pathway starting from sulphate transporters encoded by the homologs of the SUL1 and SUL2 genes in S. cerevisiae. Lack of these transporters renders Saccharomycopsis species resistant to selenate levels that are toxic to other microorganisms. We used this feature to enrich environmental samples for Saccharomycopsis species. This led to the isolation of S. schoenii, S. lassenensis and a hitherto undescribed Saccharomycopsis species with limited by-catch of other yeasts, mainly belonging to Metschnikowia and Hanseniaspora. We performed growth and predation assays to characterize the potential of these new isolates as predacious yeasts. Most Saccharomycopsis species are temperature sensitive and cannot grow at 37°C; with the exception of S. lassenensis strains. Predation assays with S. schoenii and S. cerevisiae as prey indicated that predation was enhanced at 20°C compared to 30°C. We crossed an American isolate of S. schoenii with our German isolate using marker directed breeding. Viable progeny indicated that both strains are interfertile and belong to the same biological species. S. lassenensis is heterothallic, while S. schoenii and the new Saccharomycopsis isolate, for which we suggest the name S. geisenheimensis sp. nov., are homothallic.

Phylogenetic Analysis and Comparative Genomics of Brucella abortus and Brucella melitensis Strains in Egypt

Brucellosis is a notifiable disease induced by a facultative intracellular Brucella pathogen. In this study, eight Brucella abortus and eighteen Brucella melitensis strains from Egypt were annotated and compared with RB51 and REV1 vaccines respectively. RAST toolkit in the BV-BRC server was used for annotation, revealing genome length of 3,250,377 bp and 3,285,803 bp, 3289 and 3323 CDS, 48 and 49 tRNA genes, the same number of rRNA (3) genes, 583 and 586 hypothetical proteins, 2697 and 2726 functional proteins for B. abortus and B. melitensis respectively. B. abortus strains exhibit a similar number of candidate genes, while B. melitensis strains showed some differences, especially in the SRR19520422 Faiyum strain. Also, B. melitensis clarified differences in antimicrobial resistance genes (KatG, FabL, MtrA, MtrB, OxyR, and VanO-type) in SRR19520319 Faiyum and (Erm C and Tet K) in SRR19520422 Faiyum strain. Additionally, the whole genome phylogeny analysis proved that all B. abortus strains were related to vaccinated animals and all B. melitensis strains of Menoufia clustered together and closely related to Gharbia, Dameitta, and Kafr Elshiek. The Bowtie2 tool identified 338 (eight B. abortus) and 4271 (eighteen B. melitensis) single nucleotide polymorphisms (SNPs) along the genomes. These variants had been annotated according to type and impact. Moreover, thirty candidate genes were predicted and submitted at GenBank (24 in B. abortus) and (6 in B. melitensis). This study contributes significant insights into genetic variation, virulence factors, and vaccine-related associations of Brucella pathogens, enhancing our knowledge of brucellosis epidemiology and evolution in Egypt.

Genomic characterization of a novel ureolytic bacteria, Lysinibacillus capsici TSBLM, and its application to the remediation of acidic heavy metal-contaminated soil

Soil heavy metal contamination is an essential challenge in ecological and environmental management, especially for acidic soils. Microbially induced carbonate precipitation (MICP) is an effective and environmentally friendly remediation technology for heavy metal contaminated sites, and one of the key factors for its realization lies in the microorganisms. In this study, Lysinibacillus capsici TSBLM was isolated from heavy metal contaminated soil around a gold mine, and inferred to be a novel ureolytic bacteria after phylogenomic inference and genome characterization. The urease of L. capsici TSBLM was analyzed by genetic analysis and molecular docking, and further applied this bacteria to the remediation of Cu and Pb in solution and acidic soils to investigate its biomineralization mechanism and practical application. The results revealed L. capsici TSBLM possessed a comprehensive urease gene cluster ureABCEFGD, and the encoded urease docked with urea at the lowest binding energy site (ΔG = -3.43 kcal/mol) connected to three amino acids threonine, aspartic, and alanine. The urease of L. capsici TSBLM is synthesized intracellularly but mainly functions extracellularly. L. capsici TSBLM removes Cu/Pb from the solution by generating heavy metal carbonates or co-precipitating with CaCO3 vaterite. For acidic heavy metal-contaminated soil, the carbonate-bound states of Cu and Pb increased significantly from 7 % to 16 % and from 23 % to 35 % after 30 days by L. capsici TSBLM. Soil pH improved additionally. L. capsici TSBLM maintained the dominant status in the remediated soil after 30 days, demonstrating good environmental adaptability and curing persistence. The results provided new strain resources and practical application references for the remediation of acidic heavy metal contaminated soil based on MICP.

Facultative symbiont virulence determines horizontal transmission rate without host specificity in Dictyostelium discoideum social amoebas

In facultative symbioses, only a fraction of hosts are associated with symbionts. Specific host and symbiont pairings may be the result of host-symbiont coevolution driven by reciprocal selection or priority effects pertaining to which potential symbiont is associated with a host first. Distinguishing between these possibilities is important for understanding the evolutionary forces that affect facultative symbioses. We used the social amoeba, Dictyostelium discoideum, and its symbiont, Paraburkholderia bonniea, to determine whether ongoing coevolution affects which host-symbiont strain pairs naturally cooccur within a facultative symbiosis. Relative to other Paraburkholderia, including another symbiont of D. discoideum, P. bonniea features a reduced genome size that indicates a significant history of coevolution with its host. We hypothesized that ongoing host-symbiont coevolution would lead to higher fitness for naturally cooccurring (native) host and symbiont pairings compared to novel pairings. We show for the first time that P. bonniea symbionts can horizontally transmit to new amoeba hosts when hosts aggregate together during the social stage of their life cycle. Here we find evidence for a virulence-transmission trade-off without host specificity. Although symbiont strains were significantly variable in virulence and horizontal transmission rate, hosts and symbionts responded similarly to associations in native and novel pairings. We go on to identify candidate virulence factors in the genomes of P. bonniea strains that may contribute to variation in virulence. We conclude that ongoing coevolution is unlikely for D. discoideum and P. bonniea. The system instead appears to represent a stable facultative symbiosis in which naturally cooccurring P. bonniea host and symbiont pairings are the result of priority effects.

Conservation of a Chromosome 8 Inversion and Exon Mutations Confirm Common Gulonolactone Oxidase Gene Evolution Among Primates, Including H. Neanderthalensis

Ascorbic acid functions as an antioxidant and facilitates other biochemical processes such as collagen triple helix formation, and iron uptake by cells. Animals which endogenously produce ascorbic acid have a functional gulonolactone oxidase gene (GULO); however, humans have a GULO pseudogene (GULOP) and depend on dietary ascorbic acid. In this study, the conservation of GULOP sequences in the primate haplorhini suborder were investigated and compared to the GULO sequences belonging to the primates strepsirrhini suborder. Phylogenetic analysis suggested that the conserved GULOP exons in the haplorhini primates experienced a high rate of mutations following the haplorhini/strepsirrhini divergence. This high mutation rate has decreased during the evolution of the haplorhini primates. Additionally, indels of the haplorhini GULOP sequences were conserved across the suborder. A separate analysis for GULO sequences and well-conserved GULOP sequences focusing on placental mammals identified an in-frame GULO sequence in the Brazilian guinea pig, and a potential GULOP sequence in the pika. Similar to haplorhini primates, the guinea pig and lagomorph species have experienced a high substitution rate when compared to the mammals used in this study. A shared synteny to examine the conservation of local genes near GULO/GULOP identified a conserved inversion around the GULO/GULOP locus between the haplorhini and strepsirrhini primates. Fischer's exact test did not support an association between GULOP and the chromosomal inversion. Mauve alignment showed that the inversion of the length of the syntenic block that the GULO/GULOP genes belonged to was variable. However, there were frequent rearrangements around ~ 2 million base pairs adjacent to GULOP involving the KIF13B and MSRA genes. These data may suggest that genes acquiring deleterious mutations in the coding sequence may respond to these deleterious mutations with rapid substitution rates.

Long-term evolution of Streptococcus mitis and Streptococcus pneumoniae leads to higher genetic diversity within rather than between human populations

Evaluation of the apportionment of genetic diversity of human bacterial commensals within and between human populations is an important step in the characterization of their evolutionary potential. Recent studies showed a correlation between the genomic diversity of human commensal strains and that of their host, but the strength of this correlation and of the geographic structure among human populations is a matter of debate. Here, we studied the genomic diversity and evolution of the phylogenetically related oro-nasopharyngeal healthy-carriage Streptococcus mitis and Streptococcus pneumoniae, whose lifestyles range from stricter commensalism to high pathogenic potential. A total of 119 S. mitis genomes showed higher within- and among-host variation than 810 S. pneumoniae genomes in European, East Asian and African populations. Summary statistics of the site-frequency spectrum for synonymous and non-synonymous variation and ABC modelling showed this difference to be due to higher ancestral bacterial population effective size (Ne) in S. mitis, whose genomic variation has been maintained close to mutation-drift equilibrium across (at least many) generations, whereas S. pneumoniae has been expanding from a smaller ancestral bacterial population. Strikingly, both species show limited differentiation among human populations. As genetic differentiation is inversely proportional to the product of effective population size and migration rate (Nem), we argue that large Ne have led to similar differentiation patterns, even if m is very low for S. mitis. We conclude that more diversity within than among human populations and limited population differentiation must be common features of the human microbiome due to large Ne.

Genomic investigation of multispecies and multivariant blaNDM outbreak reveals key role of horizontal plasmid transmission

OBJECTIVES

New Delhi metallo-β-lactamases (NDMs) are major contributors to the spread of carbapenem resistance globally. In Australia, NDMs were previously associated with international travel, but from 2019 we noted increasing incidence of NDM-positive clinical isolates. We investigated the clinical and genomic epidemiology of NDM carriage at a tertiary-care Australian hospital from 2016 to 2021.

METHODS

We identified 49 patients with 84 NDM-carrying isolates in an institutional database, and we collected clinical data from electronic medical record. Short- and long-read whole genome sequencing was performed on all isolates. Completed genome assemblies were used to assess the genetic setting of blaNDM genes and to compare NDM plasmids.

RESULTS

Of 49 patients, 38 (78%) were identified in 2019-2021 and only 11 (29%) of 38 reported prior travel, compared with 9 (82%) of 11 in 2016-2018 (P = .037). In patients with NDM infection, the crude 7-day mortality rate was 0% and the 30-day mortality rate was 14% (2 of 14 patients). NDMs were noted in 41 bacterial strains (ie, species and sequence type combinations). Across 13 plasmid groups, 4 NDM variants were detected: blaNDM-1, blaNDM-4, blaNDM-5, and blaNDM-7. We noted a change from a diverse NDM plasmid repertoire in 2016-2018 to the emergence of conserved blaNDM-1 IncN and blaNDM-7 IncX3 epidemic plasmids, with interstrain spread in 2019-2021. These plasmids were noted in 19 (50%) of 38 patients and 35 (51%) of 68 genomes in 2019-2021.

CONCLUSIONS

Increased NDM case numbers were due to local circulation of 2 epidemic plasmids with extensive interstrain transfer. Our findings underscore the challenges of outbreak detection when horizontal transmission of plasmids is the primary mode of spread.

Algal-bacterial consortium promotes carbon sink formation in saline environment

INTRODUCTION

The level of atmospheric CO2 has continuously been increasing and the resulting greenhouse effects are receiving attention globally. Carbon removal from the atmosphere occurs naturally in various ecosystems. Among them, saline environments contribute significantly to the global carbon cycle. Carbonate deposits in the sediments of salt lakes are omnipresent, and the biological effects, especially driven by halophilic microalgae and bacteria, on carbonate formation remain to be elucidated.

OBJECTIVES

The present study aims to characterize the carbonates formed in saline environments and demonstrate the mechanisms underlying biological-driven CO2 removal via microalgal-bacterial consortium.

METHODS

The carbonates naturally formed in saline environments were collected and analyzed. Two saline representative organisms, the photosynthetic microalga Dunaliella salina and its mutualistic halophilic bacteria Nesterenkonia sp. were isolated from the inhabiting saline environment and co-cultivated to study their biological effects on carbonates precipitation and isotopic composition. During this process, electrochemical parameters and Ca2+ flux, and expression of genes related to CaCO3 formation were analyzed. Genome sequencing and metagenomic analysis were conducted to provide molecular evidence.

RESULTS

The results showed that natural saline sediments are enriched with CaCO3 and enrichment of genes related to photosynthesis and ureolysis. The co-cultivation stimulated 54.54% increase in CaCO3 precipitation and significantly promoted the absorption of external CO2 by 49.63%. A pH gradient was formed between the bacteria and algae culture, creating 150.22 mV of electronic potential, which might promote Ca2+ movement toward D. salina cells. Based on the results of lab-scale induction and 13C analysis, a theoretical calculation indicates a non-negligible amount of 0.16 and 2.3 Tg C/year carbon sequestration in China and global saline lakes, respectively.

CONCLUSION

The combined effects of these two typical representative species have contributed to the carbon sequestration in saline environments, by promoting Ca2+ influx and increase of pH via microalgal and bacterial metabolic processes.

Genomic and taxonomic evaluation of 38 Treponema prophage sequences

BACKGROUND

Despite Spirochetales being a ubiquitous and medically important order of bacteria infecting both humans and animals, there is extremely limited information regarding their bacteriophages. Of the genus Treponema, there is just a single reported characterised prophage.

RESULTS

We applied a bioinformatic approach on 24 previously published Treponema genomes to identify and characterise putative treponemal prophages. Thirteen of the genomes did not contain any detectable prophage regions. The remaining eleven contained 38 prophage sequences, with between one and eight putative prophages in each bacterial genome. The prophage regions ranged from 12.4 to 75.1 kb, with between 27 and 171 protein coding sequences. Phylogenetic analysis revealed that 24 of the prophages formed three distinct sequence clusters, identifying putative myoviral and siphoviral morphology. ViPTree analysis demonstrated that the identified sequences were novel when compared to known double stranded DNA bacteriophage genomes.

CONCLUSIONS

In this study, we have started to address the knowledge gap on treponeme bacteriophages by characterising 38 prophage sequences in 24 treponeme genomes. Using bioinformatic approaches, we have been able to identify and compare the prophage-like elements with respect to other bacteriophages, their gene content, and their potential to be a functional and inducible bacteriophage, which in turn can help focus our attention on specific prophages to investigate further.

Novel typing scheme reveals emergence and genetic diversity of Chlamydia pecorum at the local management scale across two koala populations

To overcome shortcomings in discriminating Chlamydia pecorum strains infecting the koala (Phascolarctos cinereus) at the local level, we developed a novel genotyping scheme for this pathogen to inform koala management at a fine-scale subpopulation level. We applied this scheme to two geographically distinct koala populations in New South Wales, Australia: the Liverpool Plains and the Southern Highlands to South-west Sydney (SHSWS). Our method provides greater resolution than traditional multi-locus sequence typing, and can be used to monitor strain emergence, movement, and divergence across a range of fragmented habitats. Within the Liverpool Plains population, suspected recent introduction of a novel strain was confirmed by an absence of genetic diversity at the earliest sampling events and limited diversity at recent sampling events. Across the partially fragmented agricultural landscape of the Liverpool Plains, diversity within a widespread sequence type suggests that this degree of fragmentation may hinder but not prevent spread. In the SHSWS population, our results suggest movement of a strain from the south, where diverse strains exist, into a previously Chlamydia-free area in the north, indicating the risk of expansion towards an adjacent Chlamydia-negative koala population in South-west Sydney. In the south of the SHSWS where koala subpopulations appear segregated, we found evidence of divergent strain evolution. Our tool can be used to infer the risks of strain introduction across fragmented habitats in population management, particularly through practices such as wildlife corridor constructions and translocations.

Genome sequence of the plant-growth-promoting bacterium Bacillus velezensis EU07

Many Gram-positive spore-forming rhizobacteria of the genus Bacillus show potential as biocontrol biopesticides that promise improved sustainability and ecological safety in agriculture. Here, we present a draft-quality genome sequence for Bacillus velezensis EU07, which shows growth-promotion in tomato plants and biocontrol against Fusarium head blight. We found that the genome of EU07 is almost identical to that of the commercially used strain QST713, but identified 46 single-nucleotide differences that distinguish these strains from each other. The availability of this genome sequence will facilitate future efforts to unravel the genetic and molecular basis for EU07's beneficial properties.

Single-base resolution quantitative genome methylation analysis in the model bacterium Helicobacter pylori by enzymatic methyl sequencing (EM-Seq) reveals influence of strain, growth phase, and methyl homeostasis

BACKGROUND

Bacterial epigenetics is a rapidly expanding research field. DNA methylation by diverse bacterial methyltransferases (MTases) contributes to genomic integrity and replication, and many recent studies extended MTase function also to global transcript regulation and phenotypic variation. Helicobacter pylori is currently one of those bacterial species which possess the highest number and the most variably expressed set of DNA MTases. Next-generation sequencing technologies can directly detect DNA base methylation. However, they still have limitations in their quantitative and qualitative performance, in particular for cytosine methylation.

RESULTS

As a complementing approach, we used enzymatic methyl sequencing (EM-Seq), a technology recently established that has not yet been fully evaluated for bacteria. Thereby, we assessed quantitatively, at single-base resolution, whole genome cytosine methylation for all methylated cytosine motifs in two different H. pylori strains and isogenic MTase mutants. EM-Seq reliably detected both m5C and m4C methylation. We demonstrated that three different active cytosine MTases in H. pylori provide considerably different levels of average genome-wide single-base methylation, in contrast to isogenic mutants which completely lost specific motif methylation. We found that strain identity and changed environmental conditions, such as growth phase and interference with methyl donor homeostasis, significantly influenced quantitative global and local genome-wide methylation in H. pylori at specific motifs. We also identified significantly hyper- or hypo-methylated cytosines, partially linked to overlapping MTase target motifs. Notably, we revealed differentially methylated cytosines in genome-wide coding regions under conditions of methionine depletion, which can be linked to transcript regulation.

CONCLUSIONS

This study offers new knowledge on H. pylori global and local genome-wide methylation and establishes EM-Seq for quantitative single-site resolution analyses of bacterial cytosine methylation.

Genome Sequence Comparisons between Small and Large Colony Phenotypes of Equine Clinical Isolates of Arcanobacterium hippocoleae

Arcanobacterium hippocoleae is a Gram-positive fastidious bacterium and is occasionally isolated from the reproductive tract of apparently healthy mares (Equus caballus) or from mares with reproductive tract abnormalities. Apart from a few 16S rRNA gene-based GenBank sequences and one recent report on complete genome assembly, detailed genomic sequence and clinical experimental data are not available on the bacterium. Recently, we observed an unusual increase in the detection of the organism from samples associated with mare reproductive failures in Atlantic Canada. Two colony morphotypes (i.e., small, and large) were detected in culture media, which were identified as A. hippocoleae by MALDI-TOF mass spectrometry and 16S rRNA gene sequencing. Here, we report the whole genome sequencing and characterization of the morphotype variants. The genome length of the large phenotypes was between 2.42 and 2.43, and the small phenotype was 1.99 Mbs. The orthologous nucleotide identity between the large colony phenotypes was ~99%, and the large and small colony phenotypes was between 77.86 and 78.52%, which may warrant the classification of the two morphotypes into different species. Phylogenetic analysis based on 16S rRNA genes or concatenated housekeeping genes grouped the small and large colony variants into two different genotypic clusters. The UvrA protein, which is part of the nucleotide excision repair (NER) system, and 3-isopropoylmalate dehydratase small subunit protein expressed by the leuD gene were identified as potential virulence factors in the large and small colony morphotypes, respectively. However, detailed functional studies will be required to determine the exact roles of these and other identified hypothetical proteins in the cellular metabolism and potential pathogenicity of A. hippocoleae in mares.

Phylogenetic placement of the monotypic Baolia (Amaranthaceae s.l.) based on morphological and molecular evidence

BACKGROUND

Baolia H.W.Kung & G.L.Chu is a monotypic genus only known in Diebu County, Gansu Province, China. Its systematic position is contradictory, and its morphoanatomical characters deviate from all other Chenopodiaceae. Recent study has regarded Baolia as a sister group to Corispermoideae. We therefore sequenced and compared the chloroplast genomes of this species, and resolved its phylogenetic position based on both chloroplast genomes and marker sequences.

RESULTS

We sequenced 18 chloroplast genomes of 16 samples from two populations of Baolia bracteata and two Corispermum species. These genomes of Baolia ranged in size from 152,499 to 152,508 bp. Simple sequence repeats (SSRs) were primarily located in the LSC region of Baolia chloroplast genomes, and most of them consisted of single nucleotide A/T repeat sequences. Notably, there were differences in the types and numbers of SSRs between the two populations of B. bracteata. Our phylogenetic analysis, based on both complete chloroplast genomes from 33 species and a combination of three markers (ITS, rbcL, and matK) from 91 species, revealed that Baolia and Corispermoideae (Agriophyllum, Anthochlamys, and Corispermum) form a well-supported clade and sister to Acroglochin. According to our molecular dating results, a major divergence event between Acroglochin, Baolia, and Corispermeae occurred during the Middle Eocene, approximately 44.49 mya. Ancestral state reconstruction analysis showed that Baolia exhibited symplesiomorphies with those found in core Corispermoideae characteristics including pericarp and seed coat.

CONCLUSIONS

Comparing the chloroplast genomes of B. bracteata with those of eleven typical Chenopodioideae and Corispermoideae species, we observed a high overall similarity and a one notable noteworthy case of inversion of approximately 3,100 bp. of DNA segments only in two Atriplex and four Chenopodium species. We suggest that Corispermoideae should be considered in a broader sense, it includes Corispermeae (core Corispermoideae: Agriophyllum, Anthochlamys, and Corispermum), as well as two new monotypic tribes, Acroglochineae (Acroglochin) and Baolieae (Baolia).

Warming alters life-history traits and competition in a phage community

Host-parasite interactions are highly susceptible to changes in temperature due to mismatches in species thermal responses. In nature, parasites often exist in communities, and responses to temperature are expected to vary between host-parasite pairs. Temperature change thus has consequences for both host-parasite dynamics and parasite-parasite interactions. Here, we investigate the impact of warming (37°C, 40°C, and 42°C) on parasite life-history traits and competition using the opportunistic bacterial pathogen Pseudomonas aeruginosa (host) and a panel of three genetically diverse lytic bacteriophages (parasites). We show that phages vary in their responses to temperature. While 37°C and 40°C did not have a major effect on phage infectivity, infection by two phages was restricted at 42°C. This outcome was attributed to disruption of different phage life-history traits including host attachment and replication inside hosts. Furthermore, we show that temperature mediates competition between phages by altering their competitiveness. These results highlight phage trait variation across thermal regimes with the potential to drive community dynamics. Our results have important implications for eukaryotic viromes and the design of phage cocktail therapies.IMPORTANCEMammalian hosts often elevate their body temperatures through fevers to restrict the growth of bacterial infections. However, the extent to which fever temperatures affect the communities of phages with the ability to parasitize those bacteria remains unclear. In this study, we investigate the impact of warming across a fever temperature range (37°C, 40°C, and 42°C) on phage life-history traits and competition using a bacterium (host) and bacteriophage (parasite) system. We show that phages vary in their responses to temperature due to disruption of different phage life-history traits. Furthermore, we show that temperature can alter phage competitiveness and shape phage-phage competition outcomes. These results suggest that fever temperatures have the potential to restrict phage infectivity and drive phage community dynamics. We discuss implications for the role of temperature in shaping host-parasite interactions more widely.

Plastid phylogenomics and cytonuclear discordance in Rubioideae, Rubiaceae

In this study of evolutionary relationships in the subfamily Rubioideae (Rubiaceae), we take advantage of the off-target proportion of reads generated via previous target capture sequencing projects based on nuclear genomic data to build a plastome phylogeny and investigate cytonuclear discordance. The assembly of off-target reads resulted in a comprehensive plastome dataset and robust inference of phylogenetic relationships, where most intratribal and intertribal relationships are resolved with strong support. While the phylogenetic results were mostly in agreement with previous studies based on plastome data, novel relationships in the plastid perspective were also detected. For example, our analyses of plastome data provide strong support for the SCOUT clade and its sister relationship to the remaining members of the subfamily, which differs from previous results based on plastid data but agrees with recent results based on nuclear genomic data. However, several instances of highly supported cytonuclear discordance were identified across the Rubioideae phylogeny. Coalescent simulation analysis indicates that while ILS could, by itself, explain the majority of the discordant relationships, plastome introgression may be the better explanation in some cases. Our study further indicates that plastomes across the Rubioideae are, with few exceptions, highly conserved and mainly conform to the structure, gene content, and gene order present in the majority of the flowering plants.

Comparative genomics reveals high genetic similarity among strains of Salmonella enterica serovar Infantis isolated from multiple sources in Brazil

Background

Salmonella enterica serovar Infantis (Salmonella Infantis) is a zoonotic, ubiquitous and foodborne pathogen of worldwide distribution. Despite Brazil's relevance as a major meat exporter, few studies were conducted to characterize strains of this serovar by genomic analyses in this country. Therefore, this study aimed to assess the diversity of 80 Salmonella Infantis strains isolated from veterinary, food and human sources in Brazil between 2013 and 2018 by comparative genomic analyses. Additional genomes of non-Brazilian countries (n = 18) were included for comparison purposes in some analyses.

Methods

Analyses of whole-genome multi-locus sequence typing (wgMLST), using PGAdb-builder, and of fragmented genomes, using Gegenees, were conducted to compare the 80 Brazilian strains to the 18 non-Brazilian genomes. Pangenome analyses and calculations were performed for all Salmonella Infantis genomes analyzed. The presence of prophages was determined using PHASTER for the 80 Brazilian strains. The genome plasticity using BLAST Ring Image Generator (BRIG) and gene synteny using Mauve were evaluated for 20 selected Salmonella Infantis genomes from Brazil and ten from non-Brazilian countries. Unique orthologous protein clusters were searched in ten selected Salmonella Infantis genomes from Brazil and ten from non-Brazilian countries.

Results

wgMLST and Gegenees showed a high genomic similarity among some Brazilian Salmonella Infantis genomes, and also the correlation of some clusters with non-Brazilian genomes. Gegenees also showed an overall similarity >91% among all Salmonella Infantis genomes. Pangenome calculations revealed an open pangenome for all Salmonella Infantis subsets analyzed and a high gene content in the core genomes. Fifteen types of prophages were detected among 97.5% of the Brazilian strains. BRIG and Mauve demonstrated a high structural similarity among the Brazilian and non-Brazilian isolates. Unique orthologous protein clusters related to biological processes, molecular functions, and cellular components were detected among Brazilian and non-Brazilian genomes.

Conclusion

The results presented using different genomic approaches emphasized the significant genomic similarity among Brazilian Salmonella Infantis genomes analyzed, suggesting wide distribution of closely related genotypes among diverse sources in Brazil. The data generated contributed to novel information regarding the genomic diversity of Brazilian and non-Brazilian Salmonella Infantis in comparison. The different genetically related subtypes of Salmonella Infantis from Brazil can either occur exclusively within the country, or also in other countries, suggesting that some exportation of the Brazilian genotypes may have already occurred.

Comparative chloroplast genomes: insights into the identification and phylogeny of rapid radiation genus Rhodiola

Rhodiola L. is a genus exhibiting rapid radiation and represents a typical case for studying plastid gene adaptation in species that spread from high altitudes to low altitudes. In this study, 23 samples of 18 Rhodiola species were collected from the Qinghai-Tibetan Plateau and five scattered alpine areas, and the plastid genomes (plastomes) of these species were sequenced, annotated, and compared between high-altitude and widely distributed groups. The plastomes of Rhodiola were found to be highly conserved in terms of gene size, content, and order but highly variable in several lineage-specific features, such as codon usage bias, IR boundary shifting, and distinct repeat sequence structures binding to SSRs. Codon usage in the genes of photosystem II exhibited an obvious preference, reflecting significant environmental adaptation pressures. In this study, three repeat regions compounded with trinucleotide and mononucleotide repeats were found for the first time in R. forrestii, R. himalensis, and R. yunnanensis. High-variability regions such as ndhF, ycf1, trnH-psbA, and rpoC1-rpoB were screened, laying the foundation for the precise identification of these species. The phylogenetic analysis revealed the occurrence of cyto-nuclear discordance, likely originating from the frequent interspecific hybridization events observed within Rhodiola species during rapid radiation. Dioecious and hermaphrodite species can be broadly categorized into two subclades, probably they have different environmental adaptation strategies in response to climate change. In addition, the phylogenetic tree supported the monophyly of R. forrestii and R. yunnanensis, which compose R. Sect. Pseudorhodiola. In conclusion, plastome data enrich the genetic information available for the Rhodiola genus and may provide insight into species migration events during climate change.

RecGraph: recombination-aware alignment of sequences to variation graphs

MOTIVATION

Bacterial genomes present more variability than human genomes, which requires important adjustments in computational tools that are developed for human data. In particular, bacteria exhibit a mosaic structure due to homologous recombinations, but this fact is not sufficiently captured by standard read mappers that align against linear reference genomes. The recent introduction of pangenomics provides some insights in that context, as a pangenome graph can represent the variability within a species. However, the concept of sequence-to-graph alignment that captures the presence of recombinations has not been previously investigated.

RESULTS

In this paper, we present the extension of the notion of sequence-to-graph alignment to a variation graph that incorporates a recombination, so that the latter are explicitly represented and evaluated in an alignment. Moreover, we present a dynamic programming approach for the special case where there is at most a recombination-we implement this case as RecGraph. From a modelling point of view, a recombination corresponds to identifying a new path of the variation graph, where the new arc is composed of two halves, each extracted from an original path, possibly joined by a new arc. Our experiments show that RecGraph accurately aligns simulated recombinant bacterial sequences that have at most a recombination, providing evidence for the presence of recombination events.

AVAILABILITY AND IMPLEMENTATION

Our implementation is open source and available at https://github.com/AlgoLab/RecGraph.

Molecular analyses display the increasing diversity of Podostemaceae in China

•Four newly recorded species of Podostemaceae from southern China were identified by molecular and morphological evidence.•17 plastomes of Podostemaceae were newly sequenced and two novel polymorphic barcodes (ccsA and ndhA) detected.•Our findings reveal greater species richness (15 species from five genera) of Podostemaceae in China and supply molecular resources for research on taxonomy and phylogenomics of this enigmatic aquatic family.

Adaptations to nitrogen availability drive ecological divergence of chemosynthetic symbionts

Bacterial symbionts, with their shorter generation times and capacity for horizontal gene transfer (HGT), play a critical role in allowing marine organisms to cope with environmental change. The closure of the Isthmus of Panama created distinct environmental conditions in the Tropical Eastern Pacific (TEP) and Caribbean, offering a "natural experiment" for studying how closely related animals evolve and adapt under environmental change. However, the role of bacterial symbionts in this process is often overlooked. We sequenced the genomes of endosymbiotic bacteria in two sets of sister species of chemosymbiotic bivalves from the genera Codakia and Ctena (family Lucinidae) collected on either side of the Isthmus, to investigate how differing environmental conditions have influenced the selection of symbionts and their metabolic capabilities. The lucinid sister species hosted different Candidatus Thiodiazotropha symbionts and only those from the Caribbean had the genetic potential for nitrogen fixation, while those from the TEP did not. Interestingly, this nitrogen-fixing ability did not correspond to symbiont phylogeny, suggesting convergent evolution of nitrogen fixation potential under nutrient-poor conditions. Reconstructing the evolutionary history of the nifHDKT operon by including other lucinid symbiont genomes from around the world further revealed that the last common ancestor (LCA) of Ca. Thiodiazotropha lacked nif genes, and populations in oligotrophic habitats later re-acquired the nif operon through HGT from the Sedimenticola symbiont lineage. Our study suggests that HGT of the nif operon has facilitated niche diversification of the globally distributed Ca. Thiodiazotropha endolucinida species clade. It highlights the importance of nitrogen availability in driving the ecological diversification of chemosynthetic symbiont species and the role that bacterial symbionts may play in the adaptation of marine organisms to changing environmental conditions.

Comparative chloroplast genomes of Argentina species: genome evolution and phylogenomic implications

The genus Argentina Hill belongs to the tribe Potentilleae Sweet and contains approximately 75 species predominantly distributed in the Sino-Himalayan region and the Malesian archipelago. So far we have less knowledge on the phylogenetic relationships within Argentina owing to limited sampling of Argentina taxa or gene fragments in previous studies. Moreover, to date there is no phylogenetic study on Argentina from the perspective of comparative chloroplast (cp) genomics. Here we performed comparative genomic analyses on the cp genomes of 39 accessions representing 18 taxa of Argentina. The Argentina cp genomes presented the typical quadripartite structure, with the sizes ranging from 155 096 bp to 157 166 bp. The 39 Argentina cp genomes contained a set of 112 unique genes, comprising four ribosomal RNA (rRNA) genes, 30 transfer RNA (tRNA) genes, as well as 78 protein-coding genes (PCGs). The cp genome organization, gene content and order in Argentina were highly conserved, but some visible divergences were present in IR/SC boundary regions. Ten regions (trnH-GUG-psbA, trnG-GCC-trnfM-CAU, trnD-GUC-trnY-GUA, rpl32-trnL-UAG, atpH-atpI, rps16-trnQ-UUG, trnS-GCU-trnG-UCC, ndhF-rpl32, trnR-UCU-atpA, and accD-psaI) were identified as excellent candidate DNA markers for future studies on species identification, population genetics and phylogeny of Argentina. Our results indicated that Argentina is monophyletic. In the current sampling, the A. smithiana - A. anserina clade was sister to the remainder of Argentina. Our results corroborated the previous taxonomic treatments to transfer A. phanerophlebia and A. micropetala from the genus Sibbaldia L. to Argentina. Our results showed close relationships among A. stenophylla, A. microphylla, A. taliensis, and A. tatsienluensis, congruent with previous studies based on the morphology of these species. Twenty-six genes (rps3, rps15, rps16, rps19, rpl16, rpl20, rpl22, rpoA, rpoB, rpoC1, rpoC2, atpA, atpF, psbB, psbF, ndhA, ndhB, ndhC, ndhD, ndhF, rbcL, accD, ccsA, matK, ycf1, ycf2) were with sites under positive selection, and adaptive evolution of these genes might have played crucial roles in Argentina species adaptation to the harsh mountain environment. This study will facilitate future work on taxonomy, phylogenetics, and adaptive evolution of Argentina.