Mauve: multiple alignment of conserved genomic sequence with rearrangements

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.

Characterization of the diversity of type IV secretion system-encoding plasmids in Acinetobacter

The multi-drug resistant pathogen Acinetobacter baumannii has gained global attention as an important clinical challenge. Owing to its ability to survive on surfaces, its capacity for horizontal gene transfer, and its resistance to front-line antibiotics, A. baumannii has established itself as a successful pathogen. Bacterial conjugation is a central mechanism for pathogen evolution. The epidemic multidrug-resistant A. baumannii ACICU harbours a plasmid encoding a Type IV Secretion System (T4SS) with homology to the E. coli F-plasmid, and plasmids with homologous gene clusters have been identified in several A. baumannii sequence types. However the genetic and host strain diversity, global distribution, and functional ability of this group of plasmids is not fully understood. Using systematic analysis, we show that pACICU2 belongs to a group of almost 120 T4SS-encoding plasmids within four different species of Acinetobacter and one strain of Klebsiella pneumoniae from human and environmental origin, and globally distributed across 20 countries spanning 4 continents. Genetic diversity was observed both outside and within the T4SS-encoding cluster, and 47% of plasmids harboured resistance determinants, with two plasmids harbouring eleven. Conjugation studies with an extensively drug-resistant (XDR) strain showed that the XDR plasmid could be successfully transferred to a more divergent A. baumanii, and transconjugants exhibited the resistance phenotype of the plasmid. Collectively, this demonstrates that these T4SS-encoding plasmids are globally distributed and more widespread among Acinetobacter than previously thought, and that they represent an important potential reservoir for future clinical concern.

Loss and gain of ceftazidime-avibactam susceptibility in a non-carbapenemase-producing K1-ST23 hypervirulent Klebsiella pneumoniae

OBJECTIVES

This study aimed at revealing the underlying mechanisms of the loss and gain of ceftazidime-avibactam susceptibility in a non-carbapenemase-producing hypervirulent Klebsiella pneumoniae (hvKp).

METHODS

Here we longitudinally recovered 3 non-carbapenemase-producing K1-ST23 hvKp strains at a one-month interval (KP29105, KP29499 and KP30086) from an elderly male. Antimicrobial susceptibility testing, whole genome sequencing, transcriptomic sequencing, gene cloning, plasmid conjugation, quantitative real-time PCR (qRT-PCR), and SDS-PAGE (sodium dodecyl sulfate-polyacrylamide gel electrophoresis) were conducted.

RESULTS

Among the 3 hvKp strains, KP29105 was resistant to the third- and fourth-generation cephalosporins, KP29499 acquired resistance to both ceftazidime-avibactam and carbapenems, while KP30086 restored its susceptibility to ceftazidime-avibactam, imipenem and meropenem but retained low-level resistance to ertapenem. KP29105 and KP29499 carried plasmid-encoded genes blaCTX-M-15 and blaCTX-M-71, respectively, but KP30086 lost both. Cloning of gene blaCTX-M-71 and conjugation experiment of blaCTX-M-71-carrying plasmid showed that the transformant and transconjugant were susceptible to ceftazidime-avibactam but had a more than 8-fold increase in MICs. Supplementation with an outer membrane permeabilizer could reduce the MIC of ceftazidime-avibactam by 32 folds, indicating that porins play a key role in ceftazidime-avibactam resistance. The OmpK35 of the 3 isolates was not expressed, and the OmpK36 of KP29499 and KP30086 had a novel amino acid substitution (L359R). SDS-PAGE and qRT-PCR showed that the expression of porin OmpK36 of KP29499 and KP30086 was significantly down-regulated compared with KP29105.

CONCLUSIONS

In summary, we reported the rare ceftazidime-avibactam resistance in a non-carbapenemase-producing hvKp strain. Resistance plasmid carrying blaCTX-M-71 and mutated OmpK36 had a synergetic effect on the resistance.

Comparative genomic analysis of ESBL-selected and non-selected Escherichia coli in Australian wastewater: Elucidating differences in diversity, antimicrobial resistance, and virulence profiles

Extended-spectrum β-lactamases (ESBLs)-producing E. coli have been proposed as an indicator bacterium for antimicrobial resistance (AMR) surveillance within a OneHealth framework. However, it is important to understand the effects and potential biases ESBL-selection has on E. coli populations. Utilising whole genome sequencing, this study compared 80 ESBL-selected E. coli isolates with 201 non-selected isolates from Australian wastewater. The findings revealed significant variations between these cohorts in genetic diversity, AMR profiles, and carriage of virulence-associated genes (VAGs), plasmids, and the transmissible Locus of Stress Tolerance (tLST), a genomic island that imparts resistance to extreme heat and chlorination. The study highlights the predominance of certain sequence types (STs), particularly ST131 (75 % clade A), in ESBL-selected isolates (40 % vs 2 %) and overall the ESBL-selected isolates were largely multidrug-resistant (MDR), predominantly carrying genes for resistance to aminoglycosides, extended-spectrum β-lactams, fluoroquinolone, macrolides, sulphonamides/trimethoprim, and tetracyclines. The ESBLs identified were almost exclusively blaCTX-M genes, most commonly blaCTX-M-15 > blaCTX-M-27 > blaCTX-M-14. These were predominately carried on IncF plasmids or chromosomally (always ISEcp1 associated), in equal numbers. In contrast, 80 % of non-selected isolates carried no acquired ARGs, and none carried blaCTX-M genes. In both cohorts, extraintestinal pathogenic E. coli (ExPEC) was the dominate pathotype (35 % total) with few (4 % total) intestinal pathogenic E. coli pathotypes identified (aEPEC > ETEC > EAEC). Nevertheless, some clinically important genes were only identified in the non-selected group, namely tigecycline-resistance gene tet(X4) and AmpC ESBL blaCMY-2. Additionally, the presence of tLST, associated with higher metal resistance gene carriage (Ag, As, Cu, Hg, Ni), in a substantial portion of non-selected isolates (20 % vs 0 %), underscores environmental pressures shaping bacterial populations in wastewater ecosystems. These insights are important for developing comprehensive, less biased genomic surveillance strategies to understand and manage public health threats posed by pathogenic E. coli and AMR.

MHConstructor: a high-throughput, haplotype-informed solution to the MHC assembly challenge

The extremely high levels of genetic polymorphism within the human major histocompatibility complex (MHC) limit the usefulness of reference-based alignment methods for sequence assembly. We incorporate a short-read, de novo assembly algorithm into a workflow for novel application to the MHC. MHConstructor is a containerized pipeline designed for high-throughput, haplotype-informed, reproducible assembly of both whole genome sequencing and target capture short-read data in large, population cohorts. To-date, no other self-contained tool exists for the generation of de novo MHC assemblies from short-read data. MHConstructor facilitates wide-spread access to high-quality, alignment-free MHC sequence analysis.

An ancient ecospecies of Helicobacter pylori

Helicobacter pylori disturbs the stomach lining during long-term colonization of its human host, with sequelae including ulcers and gastric cancer1,2. Numerous H. pylori virulence factors have been identified, showing extensive geographic variation1. Here we identify a 'Hardy' ecospecies of H. pylori that shares the ancestry of 'Ubiquitous' H. pylori from the same region in most of the genome but has nearly fixed single-nucleotide polymorphism differences in 100 genes, many of which encode outer membrane proteins and host interaction factors. Most Hardy strains have a second urease, which uses iron as a cofactor rather than nickel3, and two additional copies of the vacuolating cytotoxin VacA. Hardy strains currently have a limited distribution, including in Indigenous populations in Siberia and the Americas and in lineages that have jumped from humans to other mammals. Analysis of polymorphism data implies that Hardy and Ubiquitous coexisted in the stomachs of modern humans since before we left Africa and that both were dispersed around the world by our migrations. Our results also show that highly distinct adaptive strategies can arise and be maintained stably within bacterial populations, even in the presence of continuous genetic exchange between strains.

Distribution and genetic diversity of Bombyx mori nucleopolyhedrovirus in mass-reared silkworms in Thailand

Silk is an economically and culturally important product of Thailand that is threatened by the outbreaks of the Bombyx mori nucleopolyhedrovirus (BmNPV) in silkworm rearings. BmNPV infects the larval stages of the silkworm Bombyx mori, causing fatal viral infection and preventing the larvae from reaching the pupal stage. Numerous BmNPV isolates have been described from silk producing countries from all over the world. In this study, the geographic distribution of BmNPV in Thailand was analyzed by collecting 18 samples from different regions. The BmNPV samples were analyzed by PCR and whole genome sequencing was performed for eight specimens. Their genome size ranged from 125,888 bp to 126,783 bp, comprising 138 open reading frames. Although the sequenced BmNPV genomes could be phylogenetically differentiated, a correlation between geographic and genetic distance was not observed indicating a close relationship between the BmNPV from Thailand.

First report on whole genome sequencing and comparative genomics of Salmonella enterica serovar Abortusequi isolated from Donkey in India

Salmonella enterica subspecies enterica serovar Abortusequi (S. Abortusequi) is a leading cause of abortion in equines that hinders the rapid growth of equine industry. S. Abortusequi infection in equids has re-emerged over last ten years. In the present study, S. Abortusequi was isolated and characterized from donkeys during an abortion storm in the southern peninsular region of India. Further, whole genome sequencing and phylogenomic analysis revealed that the present isolate was clustered among S. Abortusequi clade. The core genome MLST (cgMLST) analysis based on hierarchical clustering and single nucleotide polymorphism (SNP) core-genome dendrogram of the present isolate against 10 S. Abortusequi isolates revealed that the present isolate established a distinct clade compared to all previously reported isolates. A comparison of cgMLST and SNP analyses revealed the same clustering concordance between isolates. In addition, comparative genomics and phylogenetic analysis was carried out with six S. Abortusequi serovars showed a higher number of core genes than accessory genes. Further, comparative analysis of phenotype and genotype antimicrobial resistance revealed a concordance of 32% and discordance of 68% respectively.

Viromes vs. mixed community metagenomes: choice of method dictates interpretation of viral community ecology

BACKGROUND

Viruses, the majority of which are uncultivated, are among the most abundant biological entities on Earth. From altering microbial physiology to driving community dynamics, viruses are fundamental members of microbiomes. While the number of studies leveraging viral metagenomics (viromics) for studying uncultivated viruses is growing, standards for viromics research are lacking. Viromics can utilize computational discovery of viruses from total metagenomes of all community members (hereafter metagenomes) or use physical separation of virus-specific fractions (hereafter viromes). However, differences in the recovery and interpretation of viruses from metagenomes and viromes obtained from the same samples remain understudied.

RESULTS

Here, we compare viral communities from paired viromes and metagenomes obtained from 60 diverse samples across human gut, soil, freshwater, and marine ecosystems. Overall, viral communities obtained from viromes had greater species richness and total viral genome abundances than those obtained from metagenomes, although there were some exceptions. Despite this, metagenomes still contained many viral genomes not detected in viromes. We also found notable differences in the predicted lytic state of viruses detected in viromes vs metagenomes at the time of sequencing. Other forms of variation observed include genome presence/absence, genome quality, and encoded protein content between viromes and metagenomes, but the magnitude of these differences varied by environment.

CONCLUSIONS

Overall, our results show that the choice of method can lead to differing interpretations of viral community ecology. We suggest that the choice of whether to target a metagenome or virome to study viral communities should be dependent on the environmental context and ecological questions being asked. However, our overall recommendation to researchers investigating viral ecology and evolution is to pair both approaches to maximize their respective benefits. Video Abstract.

Adaptation of a methanogen to Fe0 corrosion via direct contact

Due to unique genomic adaptations, Methanococcus maripaludis Mic1c10 is highly corrosive when in direct contact with Fe0. A critical adaptation involves increased glycosylation of an extracellular [NiFe]-hydrogenase, facilitating its anchoring to cell surface proteins. Corrosive strains adapt to the constructed environment via horizontal gene transfer while retaining ancestral genes important for intraspecies competition and surface attachment. This calls for a reevaluation of how the built environment impacts methane cycling.

Mobilome impacts on physiology in the widely used non-toxic mutant Microcystis aeruginosa PCC 7806 ΔmcyB and toxic wildtype

The Microcystis mobilome is a well-known but understudied component of this bloom-forming cyanobacterium. Through genomic and transcriptomic comparisons, we found five families of transposases that altered the expression of genes in the well-studied toxigenic type-strain, Microcystis aeruginosa PCC 7086, and a non-toxigenic genetic mutant, Microcystis aeruginosa PCC 7806 ΔmcyB. Since its creation in 1997, the ΔmcyB strain has been used in comparative physiology studies against the wildtype strain by research labs throughout the world. Some differences in gene expression between what were thought to be otherwise genetically identical strains have appeared due to insertion events in both intra- and intergenic regions. In our ΔmcyB isolate, a sulfate transporter gene cluster (sbp-cysTWA) showed differential expression from the wildtype, which may have been caused by the insertion of a miniature inverted repeat transposable element (MITE) in the sulfate-binding protein gene (sbp). Differences in growth in sulfate-limited media also were also observed between the two isolates. This paper highlights how Microcystis strains continue to "evolve" in lab conditions and illustrates the importance of insertion sequences / transposable elements in shaping genomic and physiological differences between Microcystis strains thought otherwise identical. This study forces the necessity of knowing the complete genetic background of isolates in comparative physiological experiments, to facilitate the correct conclusions (and caveats) from experiments.

Genomic evidence of two-staged transmission of the early seventh cholera pandemic

The seventh cholera pandemic started in 1961 in Indonesia and spread across the world in three waves in the decades that followed. Here, we utilised genomic evidence to detail the first wave of the seventh pandemic. Genomes of 22 seventh pandemic Vibrio cholerae isolates from 1961 to 1979 were completely sequenced. Together with 152 publicly available genomes from the same period, they fell into seven phylogenetic clusters (CL1-CL7). By multilevel genome typing (MGT), all were assigned to MGT2 ST1 (Wave 1) except three isolates in CL7 which were typed as MGT2 ST2 (Wave 2). The Wave 1 seventh pandemic expanded in two stages, with Stage 1 (CL1-CL5) spread across Asia and Stage 2 (CL6 and CL7) spread to the Middle East and Africa. Three non-synonymous mutations, one each, in three regulatory genes, csrD (global regulator), acfB (chemotaxis), and luxO (quorum sensing) may have critically contributed to its pandemicity. The three MGT2 ST2 isolates in CL7 were the progenitors of Wave 2 and evolved from within Wave 1 with acquisition of a novel IncA/C plasmid. Our findings provide new insight into the evolution and transmission of the early seventh pandemic, which may aid future cholera prevention and control.

Tn4661-mediated transfer of bla CTX-M-15 from Klebsiella michiganensis to an outbreak clone of Pseudomonas aeruginosa

Carriage of CTX-M-type extended-spectrum β-lactamase (ESBL) is rare in Pseudomonas aeruginosa. During routine surveillance of an endemic ST-621 P. aeruginosa at a large hospital, isolate MRSN 100690 carrying bla CTX-M-15 was cultured from a patient (P2). This was the first detection of this ESBL in the endemic ST-621 lineage. All 1 488 bacterial isolates collected from the same facility in the 12 months prior to the incidence of 100 690 were screened for the presence of bla CTX-M-15. A set of 183 isolates was identified, in which corresponding patient metadata was evaluated for spatiotemporal overlaps with P2. The resulting three isolates, along with 100 690, were long-read sequenced using the Oxford Nanopore MinION platform to determine a potential donor of bla CTX-M-15. The screen revealed a single Klebsiella michiganensis isolate, MRSN 895358, which carried an IncA/C2 plasmid harbouring bla CTX-M-15. Notably, the patient harbouring 895358, P1, occupied the same hospital room as P2 9 months prior. Genomic alignment revealed that both isolates shared an identical 80.8 kb region containing the IncA/C2 plasmid replicon and bla CTX-M-15. This region was plasmid bound in 895 358, but chromosomally bound in 100 690 due to Tn4661-mediated transposition. ESBL bla CTX-M-15 was acquired and subsequently integrated into the chromosome of a ST-621 P. aeruginosa, likely initiated by plasmid transfer from a K. michiganensis strain.

Characterising virulence in a nontoxigenic non-O1/non-O139 Vibrio cholerae isolate imported from Vietnam

Vibrio cholerae is a major human pathogen that can cause life-threatening acute diarrhea. V. cholerae are classified according to O-antigen polysaccharide outer membrane properties, where the serotypes O1 and O139 are strains that cause pandemics and epidemics while non-O1/non-O139 usually cause mild disease. The dynamic evolution of V. cholerae involves acquisition of new virulence factors through horizontal gene transfer and formerly nontoxigenic serogroups are increasingly being reported to cause severe forms of human disease. In this study we have serotyped one isolate (ST588-CPH) of imported V. cholerae from Vietnam to Denmark and performed whole genome sequencing to identify known virulence genes and furthermore studied the pattern of virulence in closely related pathogenic strains of V. cholerae. ST558-CPH was found to be a non-O1/non-O139 strain. Initial analysis from the whole genome sequencing gave a 96,6 % match to the O139-specific wbfZ gene, but in a second analysis with a higher identification threshold, the wbfZ gene was absent. We suggest a "de novo" display of a database misannotation, which explains the conflicting results. The MLST analysis revealed that the isolate belongs to the nontoxigenic non-O1/non-O139 sequence type ST558. ST558 has recently been reported as a sequence type forming a cluster of ST's that should be monitored, as it has shown to have virulence causing moderate to severe illness. Our analysis of virulence genes identified MakA, a recently discovered toxin, which seems to be generally present in both toxigenic and nontoxigenic strains.

Deep longitudinal lower respiratory tract microbiome profiling reveals genome-resolved functional and evolutionary dynamics in critical illness

The lower respiratory tract (LRT) microbiome impacts human health, especially among critically ill patients. However, comprehensive characterizations of the LRT microbiome remain challenging due to low microbial mass and host contamination. We develop a chelex100-based low-biomass microbial-enrichment method (CMEM) that enables deep metagenomic profiling of LRT samples to recover near-complete microbial genomes. We apply the method to 453 longitudinal LRT samples from 157 intensive care unit (ICU) patients in three geographically distant hospitals. We recover 120 high-quality metagenome-assembled genomes (MAGs) and associated plasmids without culturing. We detect divergent longitudinal microbiome dynamics and hospital-specific dominant opportunistic pathogens and resistomes in pneumonia patients. Diagnosed pneumonia and the ICU stay duration were associated with the abundance of specific antibiotic-resistance genes (ARGs). Moreover, CMEM can serve as a robust tool for genome-resolved analyses. MAG-based analyses reveal strain-specific resistome and virulome among opportunistic pathogen strains. Evolutionary analyses discover increased mobilome in prevailing opportunistic pathogens, highly conserved plasmids, and new recombination hotspots associated with conjugative elements and prophages. Integrative analysis with epidemiological data reveals frequent putative inter-patient strain transmissions in ICUs. In summary, we present a genome-resolved functional, transmission, and evolutionary landscape of the LRT microbiota in critically ill patients.

Chloroplast genome evolution of Berberis (Berberidaceae): Implications for phylogeny and metabarcoding

Berberidis Radix (Sankezhen), a typical multi-origin Chinese medicinal material, originates from the dried roots of plants of the Berberis genus and is used to treat various ailments. These species have similar morphologies, potentially leading to misidentifications that can impact medicine efficacy. Therefore, developing suitable molecular markers to identify medicinal species is imperative. Furthermore, discrepancies exist in the taxonomy of the Berberis genus. In the present study, we de novo assembled the chloroplast genomes of six Berberis species (Berberis woomungensis C. Y. Wu, Berberis pruinosa Franch., Berberis thunbergii DC., Berberis chinensis Poir., Berberis wilsoniae Hemsl., and Berberis sp.) that commonly constitute Berberidis Radix and compared them with previously reported genomes. Our comparative analysis revealed similarities in genome structure, relative synonymous codon usage, amino acid frequency, repeats, and substitutions. Higher synonymous substitutions, indicative of predominant purifying selection on protein-coding genes, were observed compared to non-synonymous substitutions. However, positive selection was identified in six genes across 29 Berberis species-accD, matK, ndhD, rbcL, ycf1, and ycf2-highlighting their potential roles in adaptive responses to specific environmental conditions within the genus. Inverted repeats expansion and contraction affected the rate of mutations and were associated with the phylogenetic classification of Berberis. Our phylogenetic analysis supported the division of the Berberis complex into four genera, which corroborates previous studies involving extensive sampling. We identified the ndhD-ccsA region as the most polymorphic region and applied this region to Chinese patent medicines containing Berberidis Radix through metabarcoding. The metabarcoding analysis confirmed that five Berberis species commonly constitute Berberidis Radix in Chinese patent medicines. In conclusion, this study provides insight into the molecular evolution of the chloroplast genome and the phylogeny of the Berberis genus. In addition, metabarcoding provides insight into the species composition of Berberidis Radix in Chinese patent medicines.

Evolving Escherichia coli to use a tRNA with a non-canonical fold as an adaptor of the genetic code

All known bacterial tRNAs adopt the canonical cloverleaf 2D and L-shaped 3D structures. We aimed to explore whether alternative tRNA structures could be introduced in bacterial translation. To this end, we crafted a vitamin-based genetic system to evolve Escherichia coli toward activity of structurally non-canonical tRNAs. The system reliably couples (escape frequency <10-12) growth with the activities of a novel orthogonal histidine suppressor tRNA (HisTUAC) and of the cognate ARS (HisS) via suppression of a GTA valine codon in the mRNA of an enzyme in thiamine biosynthesis (ThiN). Suppression results in the introduction of an essential histidine and thereby confers thiamine prototrophy. We then replaced HisTUAC in the system with non-canonical suppressor tRNAs and selected for growth. A strain evolved to utilize mini HisT, a tRNA lacking the D-arm, and we identified the responsible mutation in an RNase gene (pnp) involved in tRNA degradation. This indicated that HisS, the ribosome, and EF-Tu accept mini HisT ab initio, which we confirmed genetically and through in vitro translation experiments. Our results reveal a previously unknown flexibility of the bacterial translation machinery for the accepted fold of the adaptor of the genetic code and demonstrate the power of the vitamin-based suppression system.

Comparison of organelle genomes between endangered mangrove plant Dolichandrone spathacea to terrestrial relative provides insights into its origin and adaptative evolution

Dolichandrone spathacea is a mangrove associate with high medicinal and ecological values. However, due to the dual-pressure of climate change and human activities, D. spathacea has become endangered in China. Moreover, misidentification between D. spathacea and its terrestrial relative D. cauda-felina poses further challenges to field protection and proper medicinal usage of D. spathacea. Thus, to address these problems, we sequenced and assembled mitochondrial (mt) and chloroplast (cp) genomes for both D. spathacea and D. cauda-felina. Comparative analysis revealed apparently different size and scaffold number between the two mt genomes, but a high similarity between the cp genomes. Eight regions with high sequence divergence were identified between the two cp genomes, which might be used for developing candidate DNA markers for distinguishing the two species. The splitting between D. spathacea and D. cauda-felina was inferred to occur at ~6.8 - 7.7 million years ago (Mya), which may be driven by the environment fluctuations in late Miocene. In the cp genome, 12 genes related to the expression of photosynthesis-associated proteins were detected with signatures of positive selection, which may contribute to the origin and evolutionary adaptation of Dolichandrone mangrove species. These new findings do not only enrich organelle genomic resources of Dolichandrone species, but also provide important genetic clues for improving the conservation and proper usage of endangered mangrove associate D. spathacea.

Unravelling the main genomic features of Mycoplasma equirhinis

BACKGROUND

Mycoplasma spp. are wall-less bacteria with small genomes (usually 0.5-1.5 Mb). Many Mycoplasma (M.) species are known to colonize the respiratory tract of both humans and livestock animals, where they act as primary pathogens or opportunists. M. equirhinis was described for the first time in 1975 in horses but has been poorly studied since, despite regular reports of around 14% prevalence in equine respiratory disorders. We recently showed that M. equirhinis is not a primary pathogen but could play a role in co-infections of the respiratory tract. This study was a set up to propose the first genomic characterization to better our understanding of the M. equirhinis species.

RESULTS

Four circularized genomes, two of which were generated here, were compared in terms of synteny, gene content, and specific features associated with virulence or genome plasticity. An additional 20 scaffold-level genomes were used to analyse intra-species diversity through a pangenome phylogenetic approach. The M. equirhinis species showed consistent genomic homogeneity, pointing to potential clonality of isolates despite their varied geographical origins (UK, Japan and various places in France). Three different classes of mobile genetic elements have been detected: insertion sequences related to the IS1634 family, a putative prophage related to M. arthritidis and integrative conjugative elements related to M. arginini. The core genome harbours the typical putative virulence-associated genes of mycoplasmas mainly involved in cytoadherence and immune escape.

CONCLUSION

M. equirhinis is a highly syntenic, homogeneous species with a limited repertoire of mobile genetic elements and putative virulence genes.

Comparative genomics-based insights into Xanthomonas indica, a non-pathogenic species of healthy rice microbiome with bioprotection function

Xanthomonas species are major pathogens of plants and have been studied extensively. There is increasing recognition of the importance of non-pathogenic species within the same genus. With this came the need to understand the genomic and functional diversity of non-pathogenic Xanthomonas (NPX) at the species and strain level. This study reports isolation and investigation into the genomic diversity and variation in NPX isolates, chiefly Xanthomonas indica, a newly discovered NPX species from rice. The study establishes the relationship of X. indica strains within clade I of Xanthomonads with another NPX species, X. sontii, also associated with rice seeds. Identification of highly diverse strains, open-pan genome, and systematic hyper-variation at the lipopolysaccharide biosynthetic locus when compared to pathogenic Xanthomonas indicates the acquisition of new functions for adaptation. Furthermore, comparative genomics studies established the absence of major virulence genes such as type III secretion system and effectors, which are present in the pathogens, and the presence of a known bacterial-killing type IV secretion system (X-T4SS). The diverse non-pathogenic strains of X. indica and X. sontii were found to protect rice from bacterial leaf blight pathogen, X. oryzae pv. oryzae (Xoo). The absence of phenotype of an X-T4SS mutant suggests redundancy in the genetic basis of the mechanisms involved in the bioprotection function, which may include multiple genetic loci, such as putative bacteriocin-encoding gene clusters and involvement of other factors such as nutrient and niche competition apart from induction of innate immunity through shared microbial-associated molecular patterns. The rice-NPX community and its pathogenic counterpart can be a promising model for understanding plant-microbe-microbiome interaction studies.IMPORTANCEThe Xanthomonas group of bacteria is known for its characteristic lifestyle as a phytopathogen. However, the discovery of non-pathogenic Xanthomonas (NPX) species is a major shift in understanding this group of bacteria. Multi-strain, in-depth genomic, evolutionary and functional studies on each of these NPX species are still lacking. This study on diverse non-pathogenic strains provides novel insights into genome diversity, dynamics, and evolutionary trends of NPX species from rice microbiome apart from its relationship with other relatives that form a sub-clade. Interestingly, we also uncovered that NPX species protect rice from pathogenic Xanthomonas species. The plant protection property shows their importance as a part of a healthy plant microbiome. Furthermore, finding an open pan-genome and large-scale variation at lipopolysaccharide biosynthetic locus indicates a significant role of the NPX community in host adaptation. The findings and high-quality genomic resources of NPX species and the strains will allow further systematic molecular and host-associated microbial community studies for plant health.

Genomic Analysis of Human-infecting Leptospira borgpetersenii isolates in Sri Lanka: expanded PF07598 gene family repertoire, less overall genome reduction than bovine isolates

Leptospira borgpetersenii commonly causes human leptospirosis, including severe disease. The first published analysis of L. borgpetersenii, performed on two strains of serovar Hardjo (L550 and JB197), concluded that the L. borgpetersenii genome is in the process of genome decay with functional consequences leading to a more obligately host-dependent life cycle. Yet whole genome analysis has only been carried out on few strains of L. borgpetersenii, with limited closed genomes and comprehensive analysis. Herein we report the complete, circularized genomes of seven non-Hardjo Leptospira borgpetersenii isolates from human leptospirosis patients in Sri Lanka. These isolates (all ST144) were found to be nearly identical by whole genome analysis; serotyping showed they are a novel serovar. We show that the L. borgpetersenii isolated from humans in Sri Lanka are less genomically decayed than previously reported isolates: fewer pseudogenes (N=141) and Insertion Sequence (IS) elements (N=46) compared to N=248, N=270, and N=400 pseudogenes, and N=121 and N=116 IS elements in published L. borgpetersenii Hardjo genomes (L550, JB197 and TC112). Compared to previously published L. borgpetersenii whole genome analyses showing two to three VM proteins in L. borgpetersenii isolates from cattle, rats and humans, we found that all of the human L. borgpetersenii isolates from Sri Lanka, including previously reported serovar Piyasena, have 4 encoded VM proteins, one ortholog of L. interrogans Copenhageni LIC12339 and 3 orthologs of LIC12844. Our findings of fewer pseudogenes, IS elements and expansion of the LIC12844 homologs of the PF07598 family in these human isolates suggests that this newly identified L. borgpetersenii serovar from Sri Lanka has unique pathogenicity. Comparative genome analysis and experimental studies of these L. borgpetersenii isolates will enable deeper insights into the molecular and cellular mechanisms of leptospirosis pathogenesis.

Complete chloroplast genome of a montane plant Spathoglottis aurea Lindl.: Comparative analyses and phylogenetic relationships among members of tribe collabieae

The yellow-flowered Spathoglottis aurea (tribe Collabieae; family Orchidaceae) is native to the mountainous areas of Peninsular Malaysia. The species is well known as an ornamental plant and for its role in artificial hybrid breeding. There is an interesting evolutionary relationship between S. aurea and the geographically isolated S. microchilina from Borneo that has encouraged further study of the S. aurea populations, but the genomic resource for S. aurea has not yet been reported. The present study reports the first work to characterize a chloroplast (cp) genome among the Spathoglottis genus. The complete cp genome of S. aurea was assembled from a sequence generated by the Illumina platform and analysed in comparison with other Collabieae species available in the GenBank database. The cp genome of S. aurea is 157,957 base pairs (bp) in length with guanine-cytosine (GC) content of 37.3%. The genome possessed a typical quadripartite cp genome structure with large single-copy (LSC) (86,888 bp), small single-copy (SSC) (18,125 bp) and inverted repeat (IR) (26,472 bp) sequences. A total of 134 genes were annotated, with 88 protein coding genes (PCGs), 38 transfer RNA (tRNA) genes and eight ribosomal RNA (rRNA) genes. Overall, 80 simple sequence repeats (SSR) or microsatellites were identified. Comparative analysis with other Collabieae species revealed high conservation in the cp genome arrangements with minimal difference in genome lengths. However, several mutational hotspots were also detected, with high potential to be developed as genetic markers for phylogenetic analysis. Characterization of the S. aurea cp genome revealed its conserved nature without gene loss or rearrangements when compared to other species of the Collabieae tribe. Phylogenetic analysis of Collabieae species also revealed that S. aurea has a distant evolutionary relationship to other members of the Collabieae species, despite the presence of problematic genera such as Phaius and Cephalantheropsis.

Resistance and virulence genes characteristic of a South Asia Clade (I) Candida auris strain isolated from blood in Beijing

INTRODUCTION

Candida auris is a globally disseminated invasive ascomycetous yeast, that imposes a substantial burden on healthcare systems. It has been documented to have spread to over 40 countries across six continents, necessitating in-depth comprehension through advanced techniques like Whole-Genome Sequencing.

METHOD

This study entailed the isolation and Whole-Genome Sequencing of a fluconazole-resistant C. auris strain (CA01) obtained from a patient's blood in Beijing. Genome analysis was conducted to classify the strain, and molecular docking was performed to understand the impact of mutations on drug resistance.

RESULTS

Genome analysis revealed that CA01 belongs to the South Asia Clade (I) and shares the closest genetic relationship with previously reported strains BJCA001 and BJCA002. Notably, unlike BJCA001, CA01 exhibits significant resistance to fluconazole primarily due to the A395T mutation in the ERG11 gene. Molecular docking studies demonstrated that this mutation leads to geometric changes in the active site where fluconazole binds, resulting in decreased binding affinity. Additionally, the present findings have identified several core virulence genes in C. auris, such as RBF1.

DISCUSSION

The findings from this study expand the understanding of the genetic diversity and adaptive mechanisms of C. auris within the South Asia Clade (I). The observed fluconazole resistance driven by the ERG11 mutation A395T highlights the need for heightened awareness and adaptation in clinical treatment strategies in China. This study provides critical insights into drug resistance and virulence profiles at a genetic level, which could guide future therapeutic and management strategies for C. auris infections.

Cytoplasmic genomes of Jasminum sambac reveal divergent sub-mitogenomic conformations and a large nuclear chloroplast-derived insertion

BACKGROUND

Jasminum sambac, a widely recognized ornamental plant prized for its aromatic blossoms, exhibits three flora phenotypes: single-petal ("SP"), double-petal ("DP"), and multi-petal ("MP"). The lack of detailed characterization and comparison of J. sambac mitochondrial genomes (mitogenomes) hinders the exploration of the genetic and structural diversity underlying the varying floral phenotypes in jasmine accessions.

RESULTS

Here, we de novo assembled three mitogenomes of typical phenotypes of J. sambac, "SP", "DP", and "MP-hutou" ("HT"), with PacBio reads and the "HT" chloroplast (cp) genome with Illumina reads, and verified them with read mapping and fluorescence in situ hybridization (FISH). The three mitogenomes present divergent sub-genomic conformations, with two, two, and four autonomous circular chromosomes ranging in size from 35.7 kb to 405.3 kb. Each mitogenome contained 58 unique genes. Ribosome binding sites with conserved AAGAAx/AxAAAG motifs were detected upstream of uncanonical start codons TTG, CTG and GTG. The three mitogenomes were similar in genomic content but divergent in structure. The structural variations were mainly attributed to recombination mediated by a large (~ 5 kb) forward repeat pair and several short repeats. The three jasmine cp. genomes showed a well-conserved structure, apart from a 19.9 kb inversion in "HT". We identified a 14.3 kb "HT"-specific insertion on Chr7 of the "HT" nuclear genome, consisting of two 7 kb chloroplast-derived fragments with two intact ndhH and rps15 genes, further validated by polymerase chain reaction (PCR). The well-resolved phylogeny suggests faster mitogenome evolution in J. sambac compared to other Oleaceae species and outlines the mitogenome evolutionary trajectories within Lamiales. All evidence supports that "DP" and "HT" evolved from "SP", with "HT" being the most recent derivative of "DP".

CONCLUSION

The comprehensive characterization of jasmine organelle genomes has added to our knowledge of the structural diversity and evolutionary trajectories behind varying jasmine traits, paving the way for in-depth exploration of mechanisms and targeted genetic research.

Thermophilic Hemicellulases Secreted by Microbial Consortia Selected from an Anaerobic Digester

The rise of agro-industrial activities over recent decades has exponentially increased lignocellulose biomasses (LCB) production. LCB serves as a cost-effective source for fermentable sugars and other renewable chemicals. This study explores the use of microbial consortia, particularly thermophilic consortia, for LCB deconstruction. Thermophiles produce stable enzymes that retain activity under industrial conditions, presenting a promising approach for LCB conversion. This research focused on two microbial consortia (i.e., microbiomes) that were analyzed for enzyme production using a cheap medium, i.e., a mixture of spent mushroom substrate (SMS) and digestate. The secreted xylanolytic enzymes were characterized in terms of temperature and pH optima, thermal stability, and hydrolysis products from LCB-derived polysaccharides. These enzymes showed optimal activity aligning with common biorefinery conditions and outperformed a formulated enzyme mixture in thermostability tests in the digestate. Phylogenetic and genomic analyses highlighted the genetic diversity and metabolic potential of these microbiomes. Bacillus licheniformis was identified as a key species, with two distinct strains contributing to enzyme production. The presence of specific glycoside hydrolases involved in the cellulose and hemicellulose degradation underscores these consortia's capacity for efficient LCB conversion. These findings highlight the potential of thermophilic microbiomes, isolated from an industrial environment, as a robust source of robust enzymes, paving the way for more sustainable and cost-effective bioconversion processes in biofuel and biochemical production and other biotechnological applications.

Whole Mitochondrial Genome Sequencing and Phylogenetic Tree Construction for Procypris mera (Lin 1933)

Procypris mera (Lin, 1933), also known as the Chinese ink carp, currently has a second-class protection status in China. Understanding the structure and characteristics of mitochondrial genes provides essential information for resource conservation and phylogenetic studies of P. mera. Here, we sequenced the mitochondrial genomes of three P. mera (WYL1-3) from three sites and performed phylogenetic analysis. The generated three genomes were 16,587 bp in length, comprising 13 protein-coding genes (PCGs), 22 tRNAs, two rRNAs, and two non-coding regions (control region (CR), D-loop, and light-stranded replication start OL), with a preference for codons ending in A or C. The mitochondrial genomes of WYL2 and WYL3 were identical, differing from that of WYL1 by only five single-nucleotide polymorphisms (SNPs). All mitochondrial PCGs had Ka/Ks ratios of less than one, suggesting purifying selection. Phylogenetic tree analysis based on amino acid sequences suggested that the genus Puntioplites is sister to all other genera of the subfamily Cyprinidae of China; the genus Procypris forms a monophyletic group; and the genera Carassioides, Carassius, and Cyprinus form a monophyletic group. This study contributes to our understanding of the phylogenetic relationships in subfamily Cyprininae in China and lays the foundation for resource conservation and management of P. mera.

Chloroplast genomes of Simarouba Aubl., molecular evolution and comparative analyses within Sapindales

Simarouba, a neotropical genus in the family Simaroubaceae, currently lacks comprehensive genomic data in existing databases. This study aims to fill this gap by providing genomic resources for three Simarouba species, S. amara, S. versicolor, and S. glauca. It also aims to perform comparative molecular evolutionary analyses in relation to other species within the order Sapindales. The analysis of these three Simarouba species revealed the presence of the typical quadripartite structure expected in plastomes. However, some pseudogenization events were identified in the psbC, infA, rpl22, and ycf1 genes. In particular, the CDS of the psbC gene in S. amara was reduced from 1422 bp to 584 bp due to a premature stop codon. Nucleotide diversity data pointed to gene and intergenic regions as promising candidates for species and family discrimination within the group, specifically matK, ycf1, ndhF, rpl32, petA-psbJ, and trnS-trnG. Selection signal analyses showed strong evidence for positive selection on the rpl23 gene. Phylogenetic analyses indicated that S. versicolor and S. glauca have a closer phylogenetic relationship than S. amara. We provide chloroplast genomes of three Simaruba species and use them to elucidate plastome evolution, highlight the presence of pseudogenization, and identify potential DNA barcode regions.

Genomics Insights into Mycolicibacterium Hassiacum Causing Infection in a Cat with Pyogranulomatous Dermatitis and Panniculitis

Mycolicibacterium hassiacum (homotypic synonym: Mycobacterium hassiacum) represents an ungrouped thermotolerant rapidly growing mycobacteria (RGM) species occasionally associated with infections and disease in humans. In this report, we describe a case of pyogranulomatous dermatitis and panniculitis due to M. hassiacum in an immunocompetent adult cat. To the best of our knowledge, this represents the first report of M. hassiacum infection in animals. We also report the results of the in-depth genome characterization of the isolate using a combined short- and long-read whole-genome sequencing (WGS) approach. We observed the lack of acquired-resistance genes and no evidence of mutations in housekeeping genes associated with resistance to rifampicin and isoniazid. We detected some virulence factors in our isolate, such as some associated with the interaction of mycobacteria with host cells, and the presence of multiple copies of heavy metal resistance genes (arsB, arsR, and arsL/cadL). In conclusion, M. hassiacum should be included among the RGM species associated with feline subcutaneous atypical mycobacteriosis (SAM). A reliable and fast RGM laboratory identification and characterization is important not only for an accurate etiological diagnosis but also for a correct approach to SAM treatment options.

Genomic Analyses of Methicillin-Susceptible and Methicillin-Resistant Staphylococcus pseudintermedius Strains Involved in Canine Infections: A Comprehensive Genotypic Characterization

Staphylococcus pseudintermedius is frequently associated with several bacterial infections in dogs, highlighting a One Health concern due to the zoonotic potential. Given the clinical significance of this pathogen, we performed comprehensive genomic analyses of 28 S. pseudintermedius strains isolated from canine infections throughout whole-genome sequencing using Illumina HiSeq, and compared the genetic features between S. pseudintermedius methicillin-resistant (MRSP) and methicillin-susceptible (MSSP) strains. Our analyses determined that MRSP genomes are larger than MSSP strains, with significant changes in antimicrobial resistance genes and virulent markers, suggesting differences in the pathogenicity of MRSP and MSSP strains. In addition, the pangenome analysis of S. pseudintermedius from canine and human origins identified core and accessory genomes with 1847 and 3037 genes, respectively, which indicates that most of the S. pseudintermedius genome is highly variable. Furthermore, phylogenomic analysis clearly separated MRSP from MSSP strains, despite their infection sites, showing phylogenetic differences according to methicillin susceptibility. Altogether our findings underscore the importance of studying the evolutionary dynamics of S. pseudintermedius, which is crucial for the development of effective prevention and control strategies of resistant S. pseudintermedius infections.

Genomic insights into indole-3-acetic acid catabolism in the marine algae-associated bacterium, Marinomonas sp. NFXS50

Auxins, mainly in the form of indole-3-acetic acid (IAA), regulate several aspects of plant and algal growth and development. Consequently, plant and algae-associated bacteria developed the ability to modulate IAA levels, including IAA catabolism. In this work, we present and analyse the genome sequence of the IAA-degrading and marine algae-associated bacterium, Marinomonas sp. NFXS50, analyse its IAA catabolism gene cluster and study the prevalence of IAA catabolism genes in other Marinomonas genomes. Our findings revealed the presence of homologs of the Pseudomonas iac gene cluster, implicated in IAA catabolism, in the genome of strain NFXS50; however, differences were observed in the content and organization of the Marinomonas iac gene cluster when compared to that of the model iac-containing Pseudomonas putida 1290. These variations suggest potential adaptations in the IAA catabolism pathway, possibly influenced by substrate availability and evolutionary factors. The prevalence of iac genes across several Marinomonas species underscores the significance of IAA catabolism in marine environments, potentially influencing plant/algae-bacteria interactions. This study provides novel insights into the IAA catabolism in Marinomonas, laying the groundwork for future investigations into the role of iac genes in Marinomonas physiology and the regulation of marine plant/algae-bacteria interactions.

The adaptive evolution of Quercus section Ilex using the chloroplast genomes of two threatened species

Chloroplast (cp) genome sequences have been extensively used for phylogenetic and evolutionary analyses, as many have been sequenced in recent years. Identification of Quercus is challenging because many species overlap phenotypically owing to interspecific hybridization, introgression, and incomplete lineage sorting. Therefore, we wanted to gain a better understanding of this genus at the level of the maternally inherited chloroplast genome. Here, we sequenced, assembled, and annotated the cp genomes of the threatened Quercus marlipoensis (160,995 bp) and Q. kingiana (161,167 bp), and mined these genomes for repeat sequences and codon usage bias. Comparative genomic analyses, phylogenomics, and selection pressure analysis were also performed in these two threatened species along with other species of Quercus. We found that the guanine and cytosine content of the two cp genomes were similar. All 131 annotated genes, including 86 protein-coding genes, 37 transfer RNA genes, and 8 ribosomal RNA genes, had the same order in the two species. A strong A/T bias was detected in the base composition of simple sequence repeats. Among the 59 synonymous codons, the codon usage pattern of the cp genomes in these two species was more inclined toward the A/U ending. Comparative genomic analyses indicated that the cp genomes of Quercus section Ilex are highly conserved. We detected eight highly variable regions that could be used as molecular markers for species identification. The cp genome structure was consistent and different within and among the sections of Quercus. The phylogenetic analysis showed that section Ilex was not monophyletic and was divided into two groups, which were respectively nested with section Cerris and section Cyclobalanopsis. The two threatened species sequenced in this study were grouped into the section Cyclobalanopsis. In conclusion, the analyses of cp genomes of Q. marlipoensis and Q. kingiana promote further study of the taxonomy, phylogeny and evolution of these two threatened species and Quercus.

The plastomes of Lepismium cruciforme (Vell.) Miq and Schlumbergera truncata (Haw.) Moran reveal tribe-specific rearrangements and the first loss of the trnT-GGU gene in Cactaceae

BACKGROUND

Recent studies have revealed atypical features in the plastomes of the family Cactaceae, the largest lineage of succulent species adapted to arid and semi-arid regions. Most plastomes sequenced to date are from short-globose and cylindrical cacti, while little is known about plastomes of epiphytic cacti. Published cactus plastomes reveal reduction and complete loss of IRs, loss of genes, pseudogenization, and even degeneration of tRNA structures. Aiming to contribute with new insights into the plastid evolution of Cactaceae, particularly within the tribe Rhipsalideae, we de novo assembled and analyzed the plastomes of Lepismium cruciforme and Schlumbergera truncata, two South American epiphytic cacti.

METHODS AND RESULTS

Our data reveal many gene losses in both plastomes and the first loss of functionality of the trnT-GGU gene in Cactaceae. The trnT-GGU is a pseudogene in L. cruciforme plastome and appears to be degenerating in the tribe Rhipsalideae. Although the plastome structure is conserved among the species of the tribe Rhipsalideae, with tribe-specific rearrangements, we mapped around 200 simple sequence repeats and identified nine nucleotide polymorphism hotspots, useful to improve the phylogenetic resolutions of the Rhipsalideae. Furthermore, our analysis indicated high gene divergence and rapid evolution of RNA editing sites in plastid protein-coding genes in Cactaceae.

CONCLUSIONS

Our findings show that some characteristics of the Rhipsalideae tribe are conserved, such as plastome structure with IRs containing only the ycf2 and two tRNA genes, structural degeneration of the trnT-GGU gene and ndh complex, and lastly, pseudogenization of rpl33 and rpl23 genes, both plastid translation-related genes.

Rapid Loss of Nutritional Symbionts in an Endemic Hawaiian Herbivore Radiation Is Associated with Plant Galling Habit

Insect herbivores frequently cospeciate with symbionts that enable them to survive on nutritionally unbalanced diets. While ancient symbiont gain and loss events have been pivotal for insect diversification and feeding niche specialization, evidence of recent events is scarce. We examine the recent loss of nutritional symbionts (in as little as 1 MY) in sap-feeding Pariaconus, an endemic Hawaiian insect genus that has undergone adaptive radiation, evolving various galling and free-living ecologies on a single host-plant species, Metrosideros polymorpha within the last ∼5 MY. Using 16S rRNA sequencing, we investigated the bacterial microbiomes of 19 Pariaconus species and identified distinct symbiont profiles associated with specific host-plant ecologies. Phylogenetic analyses and metagenomic reconstructions revealed significant differences in microbial diversity and functions among psyllids with different host-plant ecologies. Within a few millions of years, Pariaconus species convergently evolved the closed-gall habit twice. This shift to enclosed galls coincided with the loss of the Morganella-like symbiont that provides the essential amino acid arginine to free-living and open-gall sister species. After the Pariaconus lineage left Kauai and colonized younger islands, both open- and closed-gall species lost the Dickeya-like symbiont. This symbiont is crucial for synthesizing essential amino acids (phenylalanine, tyrosine, and lysine) as well as B vitamins in free-living species. The recurrent loss of these symbionts in galling species reinforces evidence that galls are nutrient sinks and, combined with the rapidity of the evolutionary timeline, highlights the dynamic role of insect-symbiont relationships during the diversification of feeding ecologies. We propose new Candidatus names for the novel Morganella-like and Dickeya-like symbionts.

Determining Gene Order Patterns in the Suillus and Boletales through Comparative Analysis of Their Mitogenomes

Suillus is one of the most important genera of ectomycorrhizal fungi. As a model for studying host specificity, its molecular fragments and nuclear genome have been analyzed. However, its mitochondrial genome has not yet been reported. In this study, we assembled five mitogenomes of Suillus and analyzed and compared their basic characteristics. Owing to the large number of introns as well as intergenic regions, the mitogenomic lengths of species of Suillus were greater than those of other species of Boletales. We identified two main patterns of gene order arrangement in the members of the order Boletales. The Ka/Ks values of 15 protein-coding genes were <1 for the mitochondrial genes of 39 Boletales species, indicating their conserved evolution. Phylogenetic trees, reconstructed using the mitogenomes, indicated that the genus Suillus was monophyletic. Phylogenetic results based on the internal transcribed spacer region and mitogenome were used to confirm the distribution of Suillus placidus in China. The results showed that the mitogenome was superior in distinguishing species compared with a single molecular fragment. This is the first study to investigate the mitogenome of Suillus, enriching the mitogenome information and providing basic data for the phylogeny, resource conservation, and genetic diversity of this genus.

The assembly and comparative analysis of the first complete mitogenome of Lindera aggregata

Lindera aggregata, a member belongs to the genus Lindera of Lauraceae family. Its roots and leaves have been used as traditional Chinese medicine or functional food for thousands of years. However, its mitochondrial genome has not been explored. Our aim is to sequence and assemble the mitogenome of L. aggregata to elucidate the genetic mechanism and evolutionary pathway. The results had shown that the mitogenome was extremely complex and had a unique multi-branched conformation with total size of 912,473 bp. Comprehensive analysis of protein coding genes of 7 related species showed that there were 40 common genes in their mitogenome. Interestingly, positive selection had become an important factor in the evolution of ccmB, ccmFC, rps10, rps11 and rps7 genes. Furthermore, our data highlighted the repeated trend of homologous fragment migrations between chloroplast and mitochondrial organelles, and 38 homologous fragments were identified. Phylogenetic analysis identified a tree that was basically consistent with the phylogeny of Laurales species described in the APG IV system. To sum up, this study will be helpful to the study of population genetics and evolution of Lindera species.

Comparative plastome analyses and evolutionary relationships of 25 East Asian species within the medicinal plant genus Scrophularia (Scrophulariaceae)

Backgroud

Scrophularia L., a genus of the Scrophulariaceae, is a group of important medicinal plants used for eliminating heat and detoxifying. East Asia has an abundance of potentially medicinal Scrophularia species, and it serves as a secondary diversity center of the genus. However, the genomic resources available for germplasm identification and pharmaceutical exploration of East Asian Scrophularia are insufficient, hindering its commercial and industrial development. Additionally, the interspecific relationships of most East Asian Scrophularia species remain unclear.

Methods

In this study, we sequenced the leaves of 25 East Asian species of the genus Scrophularia, assembled and annotated the complete chloroplast genomes, and subsequently performed comparative and phylogenetic analyses on these genomes.

Results and discussion

The conserved plastome length of these 25 species ranged from 151,582 bp to 153,239 bp, containing a total of 132 coding genes, including 18 duplicated genes and 114 unique genes. Through genome alignment of these 25 species, 38-53 repeated sequences and 7 shared SSRs were identified, along with regions with high nucleotide polymorphism (Pi), which could potentially serve as molecular markers for species identification. The genome structure, gene content, and arrangement showed conservation, while variations were observed in the IR boundary regions and IGS. Phylogenetic inferences based on whole plastomes or on coding sequences (CDS) only yielded congruent results. We categorized the 25 East Asian Scrophularia species into six distinct clades and further explored their interspecies relationships using morphological characteristics, such as flower color, the relative position of stamens and corolla, and plant height. This could lay a genetic basis for future resource development of Scrophularia in East Asia.

Fundamental Patterns of Structural Evolution Revealed by Chromosome-Length Genomes of Cactophilic Drosophila

A thorough understanding of adaptation and speciation requires model organisms with both a history of ecological and phenotypic study as well as a complete set of genomic resources. In particular, high-quality genome assemblies of ecological model organisms are needed to assess the evolution of genome structure and its role in adaptation and speciation. Here, we generate new genomes of cactophilic Drosophila, a crucial model clade for understanding speciation and ecological adaptation in xeric environments. We generated chromosome-level genome assemblies and complete annotations for seven populations across Drosophila mojavensis, Drosophila arizonae, and Drosophila navojoa. We use these data first to establish the most robust phylogeny for this clade to date, and to assess patterns of molecular evolution across the phylogeny, showing concordance with a priori hypotheses regarding adaptive genes in this system. We then show that structural evolution occurs at constant rate across the phylogeny, varies by chromosome, and is correlated with molecular evolution. These results advance the understanding of the D. mojavensis clade by demonstrating core evolutionary genetic patterns and integrating those patterns to generate new gene-level hypotheses regarding adaptation. Our data are presented in a new public database (cactusflybase.arizona.edu), providing one of the most in-depth resources for the analysis of inter- and intraspecific evolutionary genomic data. Furthermore, we anticipate that the patterns of structural evolution identified here will serve as a baseline for future comparative studies to identify the factors that influence the evolution of genome structure across taxa.

Phylogenomics analysis of Scutellaria (Lamiaceae) of the world

BACKGROUND

Scutellaria, a sub-cosmopolitan genus, stands as one of the Lamiaceae family's largest genera, encompassing approximately 500 species found in both temperate and tropical montane regions. Recognized for its significant medicinal properties, this genus has garnered attention as a research focus, showcasing anti-cancer, anti-inflammatory, antioxidant, and hepatoprotective qualities. Additionally, it finds application in agriculture and horticulture. Comprehending Scutellaria's taxonomy is pivotal for its effective utilization and conservation. However, the current taxonomic frameworks, primarily based on morphological characteristics, are inadequate. Despite several phylogenetic studies, the species relationships and delimitations remain ambiguous, leaving the genus without a stable and reliable classification system.

RESULTS

This study analyzed 234 complete chloroplast genomes, comprising 220 new and 14 previously published sequences across 206 species, subspecies, and varieties worldwide. Phylogenetic analysis was conducted using six data matrices through Maximum Likelihood and Bayesian Inference, resulting in a robustly supported phylogenetic framework for Scutellaria. We propose three subgenera, recommending the elevation of Section Anaspis to subgeneric rank and the merging of Sections Lupulinaria and Apeltanthus. The circumscription of Subgenus Apeltanthus and Section Perilomia needs to be reconsidered. Comparative analysis of chloroplast genomes highlighted the IR/SC boundary feature as a significant taxonomic indicator. We identified a total of 758 SSRs, 558 longer repetitive sequences, and ten highly variable regions, including trnK-rps16, trnC-petN, petN-psbM, accD-psaI, petA-psbJ, rpl32-trnL, ccsA-ndhD, rps15-ycf1, ndhF, and ycf1. These findings serve as valuable references for future research on species identification, phylogeny, and population genetics.

CONCLUSIONS

The phylogeny of Scutellaria, based on the most comprehensive sample collection to date and complete chloroplast genome analysis, has significantly enhanced our understanding of its infrageneric relationships. The extensive examination of chloroplast genome characteristics establishes a solid foundation for the future development and utilization of Scutellaria, an important medicinal plant globally.

Horizontal gene transfer from chloroplast to mitochondria of seagrasses in the yellow-Bohai seas

Seagrasses are ideal for studying plant adaptation to marine environments. In this study, the mitochondrial (mt) and chloroplast (cp) genomes of Ruppia sinensis were sequenced. The results showed an extensive gene loss in seagrasses, including a complete loss of cp-rpl19 genes in Zosteraceae, most cp-ndh genes in Hydrocharitaceae, and mt-rpl and mt-rps genes in all seagrasses, except for the mt-rpl16 gene in Phyllospadix iwatensis. Notably, most ribosomal protein genes were lost in the mt and cp genomes. The deleted cp genes were not transferred to the mt genomes through horizontal gene transfer. Additionally, a significant DNA transfer between seagrass organelles was found, with the mt genomes of Zostera containing numerous sequences from the cp genome. Rearrangement analyses revealed an unreported inversion of the cp genome in R. sinensis. Moreover, four positively selected genes (atp8, nad5, atp4, and ccmFn) and five variable regions (matR, atp4, atp8, rps7, and ccmFn) were identified.

New insights on the phylogeny, evolutionary history, and ecological adaptation mechanism in cycle-cup oaks based on chloroplast genomes

Cycle-cup oaks (Quercus section Cyclobalanopsis) are one of the principal components of forests in the tropical and subtropical climates of East and Southeast Asia. They have experienced relatively recent increases in the diversification rate, driven by changing climates and the Himalayan orogeny. However, the evolutionary history and adaptive mechanisms at the chloroplast genome level in cycle-cup oaks remain largely unknown. Therefore, we studied this problem by conducting chloroplast genomics on 50 of the ca. 90 species. Comparative genomics and other analyses showed that Quercus section Cyclobalanopsis had a highly conserved chloroplast genome structure. Highly divergent regions, such as the ndhF and ycf1 gene regions and the petN-psbM and rpoB-trnC-GCA intergenic spacer regions, provided potential molecular markers for subsequent analysis. The chloroplast phylogenomic tree indicated that Quercus section Cyclobalanopsis was not monophyletic, which mixed with the other two sections of subgenus Cerris. The reconstruction of ancestral aera inferred that Palaeotropics was the most likely ancestral range of Quercus section Cyclobalanopsis, and then dispersed to Sino-Japan and Sino-Himalaya. Positive selection analysis showed that the photosystem genes had the lowest ω values among the seven functional gene groups. And nine protein-coding genes containing sites for positive selection: ndhA, ndhD, ndhF, ndhH, rbcL, rpl32, accD, ycf1, and ycf2. This series of analyses together revealed the phylogeny, evolutionary history, and ecological adaptation mechanism of the chloroplast genome of Quercus section Cyclobalanopsis in the long river of earth history. These chloroplast genome data provide valuable information for deep insights into phylogenetic relationships and intraspecific diversity in Quercus.

A novel antigen biomarker for detection of high-level of Loa loa microfilaremia

BACKGROUND

Loiasis is a disease caused by the nematode Loa loa. Serious adverse events sometimes occur in people with heavy L. loa microfilaremia after ivermectin treatment. In regions of Central Africa where loiasis is endemic, this significantly impedes global elimination programs for lymphatic filariasis and onchocerciasis that use mass distribution of ivermectin. Improved diagnostic tests to identify individuals at increased risk of serious adverse events could facilitate efforts to eliminate lymphatic filariasis and onchocerciasis in this region.

METHODS AND FINDINGS

We previously identified the L. loa protein Ll-Bhp-1 in loiasis patient sera. Here, we further characterize Ll-Bhp-1 and report development of an antigen capture ELISA to detect this antigen. This assay detected Ll-Bhp-1 in 74 of 116 (63.8%) loiasis patient sera. Ll-Bhp-1 levels were significantly correlated with L. loa microfilarial counts, and the sensitivity of the assay was highest for samples from people with high counts, (94% and 100% in people with ≥20,000 and ≥50,000 microfilaria per milliliter of blood, respectively). The antigen was not detected in 112 sera from people with other filarial infections, or in 34 control sera from the USA.

CONCLUSIONS

This Ll-Bhp-1 antigen assay is specific for loiasis, and highly sensitive for identifying people with high L. loa microfilarial counts who are at increased risk for serious adverse events after ivermectin treatment. L. loa antigen detection has the potential to facilitate loiasis mapping efforts and programs to eliminate lymphatic filariasis and onchocerciasis in Central Africa.

Mammaliicoccus sciuri's Pan-Immune System and the Dynamics of Horizontal Gene Transfer Among Staphylococcaceae: a One-Health CRISPR Tale

Recently emancipated from the Staphylococcus genus due to genomic differences, Mammaliicoccus sciuri, previously classified as an occasional pathogen, emerges as a significant player in the landscape of resistance gene dissemination among Staphylococcaceae. Despite its classification, its role remained enigmatic. In this study, we delved into the genomic repertoire of M. sciuri to unravel its contribution to resistance and virulence gene transfer in the context of One Health. Through comprehensive analysis of publicly available genomes, we unveiled a diverse pan-immune system adept at defending against exogenous genetic elements, yet concurrently fostering horizontal gene transfer (HGT). Specifically, exploration of CRISPR-Cas systems, with spacer sequences as molecular signatures, elucidated a global dissemination pattern spanning environmental, animal, and human hosts. Notably, we identified the integration of CRISPR-Cas systems within SCCmecs (Staphylococcal Cassette Chromosome mec), harboring key genes associated with pathogenicity and resistance, especially the methicillin resistance gene mecA, suggesting a strategic adaptation to outcompete other mobile genetic elements. Our findings underscored M. sciuri's active engagement in HGT dynamics and evolutionary trajectories within Staphylococcaceae, emphasizing its central role in shaping microbial communities and highlighting the significance of understanding its implications in the One Health framework, an interdisciplinary approach that recognizes the interconnectedness of human, animal, and environmental health to address global health challenges.

Exploitation and Application of a New Genetic Sex Marker Based on Intron Insertion Variation of erc2 Gene in Oplegnathus punctatus

Spotted knifejaw (Oplegnathus punctatus), one of the most valuable mariculture species, grows with significant sexual dimorphism, with males growing significantly faster than females. O. punctatus not only has excellent growth characteristics and high food value, but also shows high economic value in aquaculture, which has become a hotspot in the field of aquaculture. The current insufficiency of sex marker identification in O. punctatus restricts the process of its unisexual breeding. Rapid identification of sex will help to study the mechanisms of sex determination and accelerate the development of sex-controlled breeding. With the completion of the sequencing of the male and female genomes of O. punctatus, the efficient and precise development of genetic sex markers has been made possible. In this study, we used genome-wide information combined with molecular biology techniques from marker sequences to further establish a rapid method for DNA insertion variant detection in the intron of O. punctatus erc2 gene, which can be used to rapidly, accurately, and efficiently identify whether DNA insertion occurs in the intron of O. punctatus erc2 gene to be detected, and to identify the sex of O. punctatus to be detected. It could also be distinguished by agarose gel electrophoresis, which would shorten the time for accurate identification and improves the detection efficiency. Homozygous comparison of male and female individuals showed that the length of the DNA fragment of the erc2 gene was 239 bp on chromosome X1 and 1173 bp on chromosome Y. It can therefore be inferred that a 934 bp insertion fragment exists on the Y chromosome. The PCR amplification results showed that two DNA fragments of 1173 bp and 239 bp could be amplified in male O. punctatus, and the 1173 bp fragment was a marker fragment specific to the variant intron erc2 gene, while only a single DNA fragment of 239 bp was amplified in female O. punctatus. It has important significance and application value in the study of neurotransmitter transmission and environmental adaptability of female and male fish based on erc2 gene, as well as the identification of male and female sex, the preparation of high male fry, and family breeding.

Understanding the biosynthesis, metabolic regulation, and anti-phytopathogen activity of 3,7-dihydroxytropolone in Pseudomonas spp

The genus Pseudomonas is a prolific source of specialized metabolites with significant biological activities, including siderophores, antibiotics, and plant hormones. These molecules play pivotal roles in environmental interactions, influencing pathogenicity, inhibiting microorganisms, responding to nutrient limitation and abiotic challenges, and regulating plant growth. These properties mean that pseudomonads are suitable candidates as biological control agents against plant pathogens. Multiple transposon-based screens have identified a Pseudomonas biosynthetic gene cluster (BGC) associated with potent antibacterial and antifungal activities, which produces 7-hydroxytropolone (7-HT). In this study, we show that this BGC also makes 3,7-dihydroxytropolone (3,7-dHT), which has strong antimicrobial activity toward Streptomyces scabies, a potato pathogen. Through metabolomics and reporter assays, we unveil the involvement of cluster-situated genes in generating phenylacetyl-coenzyme A, a key precursor for tropolone biosynthesis via the phenylacetic acid catabolon. The clustering of these phenylacetic acid genes within tropolone BGCs is unusual in other Gram-negative bacteria. Our findings support the interception of phenylacetic acid catabolism via an enoyl-CoA dehydratase encoded in the BGC, as well as highlighting an essential role for a conserved thioesterase in biosynthesis. Biochemical assays were used to show that this thioesterase functions after a dehydrogenation-epoxidation step catalyzed by a flavoprotein. We use this information to identify diverse uncharacterized BGCs that encode proteins with homology to flavoproteins and thioesterases involved in tropolone biosynthesis. This study provides insights into tropolone biosynthesis in Pseudomonas, laying the foundation for further investigations into the ecological role of tropolone production.IMPORTANCEPseudomonas bacteria produce various potent chemicals that influence interactions in nature, such as metal-binding molecules, antibiotics, or plant hormones. This ability to synthesize bioactive molecules means that Pseudomonas bacteria may be useful as biological control agents to protect plants from agricultural pathogens, as well as a source of antibiotic candidates. We have identified a plant-associated Pseudomonas strain that can produce 3,7-dihydroxytropolone, which has broad biological activity and can inhibit the growth of Streptomyces scabies, a bacterium that causes potato scab. Following the identification of this molecule, we used a combination of genetic, chemical, and biochemical experiments to identify key steps in the production of tropolones in Pseudomonas species. Understanding this biosynthetic process led to the discovery of an array of diverse pathways that we predict will produce new tropolone-like molecules. This work should also help us shed light on the natural function of antibiotics in nature.

Comparative Genomic Analysis of Mycobacterium tuberculosis Isolates Circulating in North Santander, Colombia

Tuberculosis (TB) is an important infectious disease in relation to global public health and is caused species of the Mycobacterium tuberculosis complex (MTBC). In this study, we used whole-genome sequencing (WGS) and comparative genomics to investigate the genetic diversity of M. tuberculosis (Mtb) isolates circulating in North Santander (NS), Colombia. WGS was used for the phylogenetic and lineage characterization of 18 isolates of Mtb typed with orphan genotypes from 11 municipalities of NS between 2015 and 2018. The isolates studied were included in six sublineages from L4; the most frequent were 4.1.2.1, 4.3.3, and 4.3.4.2, corresponding to a proportion of 22.2%. The genome analysis conducted allowed the identification of a set of genetic variants mainly associated with determinants of virulence and evasion of the immune system (PPE34 and PE_PGRS2); adaptation and survival (PGL/p-HBAD); stress response (sigJ and sigM); geographic variability (PPE34); and carbohydrate and lipid metabolism (aldA, rocA, and cyp144). This is the first description of the molecular epidemiology of Mtb isolates circulating in NS achieved through WGS. It was possible to perform comparative genomics analyses between Mtb isolates against the universal reference H37Rv and Colombian UT205 genome, which can help us to understand the local genetic diversity and is relevant for epidemiological studies, providing insight into TB transmission dynamics in NS.

Exploring the mitochondrial genomes and phylogenetic relationships of trans-Andean Bryconidae species (Actinopterygii: Ostariophysi: Characiformes)

Comparative mitogenomics and its evolutionary relationships within Bryconidae remains largely unexplored. To bridge this gap, this study assembled 15 mitogenomes from 11 Bryconidae species, including five newly sequenced. Salminus mitogenomes, exceeding 17,700 bp, exhibited the largest size, contrasting with a median size of 16,848 bp in the remaining species (Brycon and Chilobrycon). These mitogenomes encode 37 typical mitochondrial genes, including 13 protein-coding, 2 ribosomal RNA, and 22 transfer RNA genes, and exhibit the conserved gene arrangement found in most fish species. Phylogenetic relationships, based on the maximum-likelihood method, revealed that the trans-Andean species (found in northwestern South America) clustered into two main sister clades. One clade comprised the trans-Andean species from the Pacific slope, Brycon chagrensis and Chilobrycon deuterodon. The other clade grouped the trans-Andean species from the Magdalena-Cauca Basin Brycon moorei and Salminus affinis, with their respective cis-Andean congeners (found in eastern South America), with Brycon rubricauda as its sister clade. Since the current members of Brycon are split in three separated lineages, the systematic classification of Bryconidae requires further examination. This study provides novel insights into mitogenome characteristics and evolutionary pathways within Bryconidae, standing as crucial information for prospective phylogenetic and taxonomic studies, molecular ecology, and provides a valuable resource for environmental DNA applications.

DNA polymerase swapping in Caudoviricetes bacteriophages

BACKGROUND

Viruses with double-stranded (ds) DNA genomes in the realm Duplodnaviria share a conserved structural gene module but show a broad range of variation in their repertoires of DNA replication proteins. Some of the duplodnaviruses encode (nearly) complete replication systems whereas others lack (almost) all genes required for replication, relying on the host replication machinery. DNA polymerases (DNAPs) comprise the centerpiece of the DNA replication apparatus. The replicative DNAPs are classified into 4 unrelated or distantly related families (A-D), with the protein structures and sequences within each family being, generally, highly conserved. More than half of the duplodnaviruses encode a DNAP of family A, B or C. We showed previously that multiple pairs of closely related viruses in the order Crassvirales encode DNAPs of different families.

METHODS

Groups of phages in which DNAP swapping likely occurred were identified as subtrees of a defined depth in a comprehensive evolutionary tree of tailed bacteriophages that included phages with DNAPs of different families. The DNAP swaps were validated by constrained tree analysis that was performed on phylogenetic tree of large terminase subunits, and the phage genomes encoding swapped DNAPs were aligned using Mauve. The structures of the discovered unusual DNAPs were predicted using AlphaFold2.

RESULTS

We identified four additional groups of tailed phages in the class Caudoviricetes in which the DNAPs apparently were swapped on multiple occasions, with replacements occurring both between families A and B, or A and C, or between distinct subfamilies within the same family. The DNAP swapping always occurs "in situ", without changes in the organization of the surrounding genes. In several cases, the DNAP gene is the only region of substantial divergence between closely related phage genomes, whereas in others, the swap apparently involved neighboring genes encoding other proteins involved in phage genome replication. In addition, we identified two previously undetected, highly divergent groups of family A DNAPs that are encoded in some phage genomes along with the main DNAP implicated in genome replication.

CONCLUSIONS

Replacement of the DNAP gene by one encoding a DNAP of a different family occurred on many independent occasions during the evolution of different families of tailed phages, in some cases, resulting in very closely related phages encoding unrelated DNAPs. DNAP swapping was likely driven by selection for avoidance of host antiphage mechanisms targeting the phage DNAP that remain to be identified, and/or by selection against replicon incompatibility.

Plastomes of Nothofagus reflect a shared biogeographic history in Patagonia

Plastomes are used in phylogenetic reconstructions because of their relatively conserved nature. Nonetheless, some limitations arise, particularly at lower taxonomic levels due to reduced interspecific polymorphisms and frequent hybridization events that result in unsolved phylogenies including polytomies and reticulate evolutionary patterns. Next-generation sequencing technologies allow access to genomic data and strongly supported phylogenies, yet biased topologies may be obtained due to insufficient taxon sampling. We analyze the hypothesis that intraspecific plastome diversity reflects biogeographic history and hybridization cycles among taxa. We generated 12 new plastome sequences covering distinct latitudinal locations of all species of subgenus Nothofagus from North Patagonia. Chloroplast genomes were assembled, annotated, and searched for simple sequence repeats (SSRs). Phylogenetic reconstructions included species and sampled locations. The six Nothofagus species analyzed were of similar size and structure; only Nothofagus obliqua of subgenus Lophozonia, used as an outgroup, presented slight differences in size. We detected a variable number of SSRs in distinct species and locations. Phylogenetic analyses of plastomes confirmed that subgenus Nothofagus organizes into two monophyletic clades each consisting of individuals of different species. We detected a geographic structure within subgenus Nothofagus and found evidence of local chloroplast sharing due to past hybridization, followed by adaptive introgression and ecological divergence. These contributions enrich the comprehension of transversal evolutionary mechanisms such as chloroplast capture and its implications for phylogenetic and phylogenomic analyses.

Yushaniadezhui (Poaceae, Bambusoideae), a new bamboo species from Yunnan, China

A new bamboo species, Yushaniadezhui, from Kunming, Yunnan, China is described and illustrated in this paper. The new species used to be misidentified as Y.polytricha. Based on careful comparison of morphological features and molecular phylogeny evidence, we confirmed its identity as a new member of the genus Yushania. Yushaniadezhui resembles Y.maculata, Y.polytricha and Y.weixiensis in several aspects, such as culm height and branch complement structure. However, the glabrous culm leaf sheaths and internodes, the absence of auricles and oral setae on most foliage leaves, except the one-year-old foliage leaves, the pubescence on the adaxial surface of the one-year-old foliage leaves and its limestone habitat preference can readily distinguish this new species from its related taxa. Moreover, we emphasise that individuals from various populations and molecular markers with different inheritance patterns for phylogeny reconstruction should be included in new species discovery, especially in plant groups with complex evolutionary histories.

A novel lytic phage infecting MDR Salmonella enterica and its application as effective food biocontrol

Salmonella enterica is a foodborne pathogen associated with both typhoid and non-typhoid illness in humans and animals. This problem is further exacerbated by the emergence of antibiotic-resistant strains of Salmonella enterica. Therefore, to meet public health and safety, there is a need for an alternative strategy to tackle antibiotic-resistant bacteria. Bacteriophages or (bacterial viruses), due to their specificity, self-dosing, and antibiofilm activity, serve as a better approach to fighting against drug-resistant bacteria. In the current study, a broad-host range lytic phage phiSalP219 was isolated against multidrug-resistant Salmonella enterica serotypes Paratyphi from a pond water sample. Salmonella phage phiSalP219 was able to lyse 28/30 tested strains of Salmonella enterica. Salmonella phage phiSalP219 exhibits activity in acidic environments (pH3) and high temperatures (70°C). Electron microscopy and genome analysis revealed that phage phiSalP219 is a member of class Caudoviricetes. The genome of Salmonella phage phiSalP219 is 146Kb in size with 44.5% GC content. A total of 250 Coding Sequence (CDS) and 25 tRNAs were predicted in its genome. Predicted open reading frames (ORFs) were divided into five groups based on their annotation results: (1) nucleotide metabolism, (2) DNA replication and transcription, (3) structural proteins, (4) lysis protein, and (5) other proteins. The absence of lysogeny-related genes in their genome indicates that Salmonella phage phiSalP219 is lytic in nature. Phage phiSalP219 was also found to be microbiologically safe (due to the absence of toxin or virulence-related genes) in the control of Salmonella enterica serovar Typhimurium infections in the ready-to-eat meat and also able to eradicate biofilm formed by the same bacterium on the borosilicate glass surface.