Ten years of pan-genome analyses

Phylogenomic and Pangenomic Assessment of a Mediterranean Strain of Raphidiopsis raciborskii Extends Knowledge of the Global Distribution and Characteristics of a Potentially Toxigenic Cyanobacterium

Among potentially toxigenic cyanobacteria, Raphidiopsis raciborskii has attracted considerable attention due to its ability to produce massive blooms and its recent spread to temperate regions. In this work, we reported for the first time a taxonomic and functional assessment of a R. raciborskii strain isolated from the Mediterranean region, contributing to filling a gap in the global distribution and characteristics of this species. The strain LT_0923 was isolated from Lake Trasimeno, a large and shallow lake in central Italy. The phylogenomic analyses based on selected marker genes and the core genome obtained from a pangenomic analysis based on a selection of available high-quality genomes showed a strong correspondence of the Lake Trasimeno strain with the North American and, at a lower average nucleotide identity, with the South American genomes. The LT_0923 genome did not show the presence of gene clusters encoding legacy cyanotoxins or emerging toxigenic compounds. The open pangenome and the large fraction of distinct gene families identified in the cloud and partly shell genome, enriched with genes specialised in environmental-specific functions and defence mechanisms, are consistent with the development of Raphidiopsis in geographically distinct regions.

Distantly related bacteria share a rigid proteome allocation strategy with flexible enzyme kinetics

Bacteria are known to allocate their proteomes according to how fast they grow, and the allocation strategies employed strongly affect bacterial adaptation to different environments. Much of what is currently known about proteome allocation is based on extensive studies of the model organism Escherichia coli. It is not clear how much of E. coli's proteome allocation strategy is applicable to other species, particularly since different species can grow at vastly different rates even in the same growth condition. In this study, we investigate differences in nutrient-dependent proteome allocation programs adopted by several distantly related bacterial species, including Vibrio natriegens, one of the fastest-growing bacteria known. Extensive quantitative proteome characterization across conditions reveals an invariant allocation program in response to changing nutrients despite systemic, species-specific differences in enzyme kinetics. This invariant program is not organized according to the growth rate but is based on a common internal metric of nutrient quality after scaling away species-specific differences in enzyme kinetics, with the faster species behaving as if it is growing under a higher temperature. The flexibility of enzyme kinetics and the rigidity of proteome allocation programs across species defy common notions of evolvability and resource optimization. Our results suggest the existence of a blueprint of proteome allocation shared by diverse bacterial species, with implications on common underlying regulatory strategies. Further knowledge on the existence and organization of such phylogeny-transcending relations also promises to simplify the bottom-up description and understanding of bacterial behaviors in ecological communities.

Lake Malawi cichlid pangenome graph reveals extensive structural variation driven by transposable elements

Pangenome methods have the potential to uncover hitherto undiscovered sequences missing from established reference genomes, making them useful to study evolutionary and speciation processes in diverse organisms. The cichlid fishes of the East African Rift Lakes represent one of nature's most phenotypically diverse vertebrate radiations, but single-nucleotide polymorphism (SNP)-based studies have revealed little sequence difference, with 0.1%-0.25% pairwise divergence between Lake Malawi species. These were based on aligning short reads to a single linear reference genome and ignored the contribution of larger-scale structural variants (SVs). We constructed a pangenome graph that integrates six new and two existing long-read genome assemblies of Lake Malawi haplochromine cichlids. This graph intuitively represents complex and nested variation between the genomes and reveals that the SV landscape is dominated by large insertions, many exclusive to individual assemblies. The graph incorporates a substantial amount of extra sequence across seven species, the total size of which is 33.1% longer than that of a single cichlid genome. Approximately 4.73% to 9.86% of the assembly lengths are estimated as interspecies structural variation between cichlids, suggesting substantial genomic diversity underappreciated in SNP studies. Although coding regions remain highly conserved, our analysis uncovers a significant proportion of SV sequences as transposable element (TE) insertions, especially DNA, LINE, and LTR TEs. These findings underscore that the cichlid genome is shaped both by small-nucleotide mutations and large, TE-derived sequence alterations, both of which merit study to understand their interplay in cichlid evolution.

Genomic insights into the population structure, antimicrobial resistance, and virulence of Brachyspira hyodysenteriae from diverse geographical regions

Swine dysentery, caused by the anaerobic spirochete Brachyspira hyodysenteriae, leads to mucohemorrhagic diarrhea in grower-finisher pigs, impacting swine production. Knowledge regarding its genomic epidemiology is limited. We performed a whole-genome sequence analysis for 251 B. hyodysenteriae genomes from 10 countries, including 117 isolates sequenced in this study. Phylogenomic analysis based on core-genome single nucleotide polymorphisms (SNPs) revealed nine lineages, with L7 (72 isolates, 28.69%), L9 (67 isolates, 26.69%), and L2 (53 isolates, 21.12%) predominating. Geographical clustering was observed with distinct lineage distributions. Multilocus sequence typing identified 69 sequence types (STs), including 20 novel STs across 251 genomes. Association between specific lineages, STs, and geographical regions was evident, highlighting evolutionary and regional patterns. The pan-genome analysis identified 5,231 genes, categorized into core (1,648), accessory (2,619), and unique (964) components. Functional annotation linked core genes to essential cellular processes, while accessory and unique genes were enriched in genetic variability, defense mechanisms, and secondary metabolism. The pan-genome exhibited a high proportion of hypothetical genes, necessitating further functional characterization. Antimicrobial resistance (AMR) screening detected the tva(A) and lnu(C) genes associated with tiamulin and lincomycin resistance, respectively, in specific lineages and STs. Virulence factor analysis identified genes linked to hemolysin production, iron uptake, and survival in host environments in most isolates, with a subset of genes demonstrating lineage-specific associations that are further linked to pathogenic potential. This comprehensive genomic epidemiological analysis elucidates the genetic diversity, antimicrobial resistance, and virulence of B. hyodysenteriae globally, enhancing understanding of its epidemiology and guiding interventions to mitigate swine dysentery.

IMPORTANCE

Brachyspira hyodysenteriae, the primary causative agent of swine dysentery, remains a less-studied pathogen than other bacterial species that impact animal health. This study uses whole-genome sequencing and advanced phylogenomic approaches to reveal the genetic diversity and geographical distribution of B. hyodysenteriae isolates, focusing on U.S. populations. The identification of nine distinct phylogenetic lineages and associated sublineages highlights the pathogen's complex population structure and regional variation. Importantly, the study detects AMR genes, including tva(A) and lnu(C), linked to tiamulin and lincomycin resistance, that may pose significant challenges to disease management. The analysis also identifies virulence-associated genes, shedding light on molecular mechanisms underlying pathogenicity. By combining core-genome SNP phylogenies with multilocus sequence typing and accessory genome insights, this work provides a robust framework for a better understanding of B. hyodysenteriae evolution. Overall, these findings underscore the importance of genomic surveillance in informing control strategies and improving swine health worldwide.

Salty secrets of Halobacterium salinarum AD88: a new archaeal ecotype isolated from Cuatro Cienegas Basin

The Cuatro Cienegas Basin (CCB) in Mexico, represents a unique ecological habitat, characterized by extreme and fluctuating conditions, providing a window into ancient evolutionary processes. This basin, characterized by hypersalinity and phosphorus scarcity, harbors diverse microbial communities that exhibit remarkable adaptations to oligotrophic conditions. Among these, Halobacterium salinarum, a halophilic archaeon known for its polyploid genome and metabolic versatility, has been extensively studied as a model for extremophile survival. However, only a limited number of H. salinarum strains have been successfully cultured and characterized to date. Here, we report the isolation and genomic analysis of a novel Halobacterium salinarum strain, AD88, from microbial mats at the Archaean Domes site in the CCB. This strain displays unique genomic features, including smaller plasmid sizes and distinctive metabolic pathways for phosphorus and sulfur utilization. Comparative analyses with other Halobacterium strains revealed genetic innovations, such as genes involved in sulfolipid biosynthesis, enabling membrane stability in phosphorus-depleted environments, and adaptations for horizontal gene transfer, which facilitate genomic flexibility in response to environmental pressures. This study reveals that H. salinarum AD88 is the first recorded diploid strain of Halobacterium, a feature previously undocumented in this genus. Phylogenomic reconstruction positioned AD88 tightly within the Halobacterium clade, reflecting its evolutionary history within the genus. Pangenome analysis further highlighted the open nature of the Halobacterium genus, with AD88 contributing novel accessory genes linked to ecological specialization. These findings emphasize the evolutionary significance of the CCB as a natural laboratory for studying microbial adaptation and expand our understanding of archaeal genomic diversity and functional innovation under extreme conditions.

Comparative Genomic and Transcriptomic Analysis Provides New Insights into the Aflatoxin B1 Biodegradability by Kocuria rosea from Deep Sea

As a mycotoxin, aflatoxin B1 (AFB1) has attracted increasing attention due to its highly toxic effects, such as teratogen, carcinogen, and mutagen. AFB1-producing fungi have been found in both terrestrial and marine sources. Over the last two years, the marine-derived bacteria Kocuria rosea has shown the ability to degrade AFB1. However, no studies have yet explored the aflatoxin degradation potential in the genus Kocuria, and the molecular basis of aflatoxin metabolisms by strain has not yet been investigated. In this study, we first compared and analyzed the genomic features of 13 bacteria of the genus Kocuria by comparative genomic analysis and investigated the evolutionary patterns (evolutionary selection pressure analysis and gene family expansion analysis) and pan and core genomes of Kocuria strains. It was found that Kocuria genus strains showed high potential in the bioremediation of aromatic compounds and aflatoxin degradation. In addition, this study revealed 1421 differentially expressed genes and 752 upregulated genes between the aflatoxin group and the control group through transcriptome analysis and conducted functional annotation and analysis of functional enrichment. According to the expression and functional analysis of upregulated genes, the enzymes and genes (cytochrome P450, oxidoreductase, intradiol ring-cleavage dioxygenase, hydrolase, and monooxygenase) involved in the aflatoxin degradation of Kocuria rosea were screened. These findings indicate the potential of Kocuria rosea in the remediation of aflatoxin contamination and provide a foundation and direction for the further excavation of novel degradation genes, exploration of degradation mechanisms, and genomic modification.

Pangenome Analysis of the Plant Pathogen Pseudomonas syringae Reveals Unique Natural Products for Niche Adaptation

Pseudomonas syringae is a soil-dwelling bacterium that exhibits remarkable niche adaptability, and it is known for its devastating impact as a plant pathogen. This bacterium has an outstanding capability to produce a wide array of biologically active natural products. P. syringae coexists with amoebal predators and fungal strains, which drives the production of secondary metabolites for predator evasion in addition to niche adaptation. In this study, we conducted a broad pangenomic analysis of 18 taxonomically distinct P. syringae strains, leading to the identification of 231 biosynthetic gene clusters (BGCs). Among these, nonribosomal peptide synthetases (NRPSs) were particularly abundant, indicating their potential significance within this ecological context. We discovered and elucidated the structures of two novel classes of bioactive compounds, the syrilipamides and chlorosecimides. Furthermore, a bioinformatic analysis enabled the identification of an undescribed halogenase, SecA, essential for the chlorination of secimide A. We observed that syrilipamides and secimides and in particular mixtures thereof, exhibit amoebicidal activities. Additionally, secimides showed selective antifungal activity. These findings provide valuable insights into the ecological roles of P. syringae natural products and highlight their potential for biotechnological and therapeutic applications.

CropGene: a software package for the analysis of genomic and transcriptomic data of agricultural plants

Currently, the breeding of agricultural plants is increasingly based on the use of molecular biological data on genetic sequences, which makes it possible to significantly accelerate the breeding process, create new plant varieties through genomic editing. These data have a large volume, variety and require a large amount of resources, both labor and computing, to analyze the costs. Data analysis of such volume and complexity can be effective only when using modern bioinformatics methods, which include algorithms for identifying genes, predicting their function, and evaluating the effect of mutation on plant phenotype. Such an analysis has recently become impossible without the use of integrated software systems that solve problems of different levels by executing computational pipelines. The paper describes the CropGene software package developed for the comprehensive analysis of genomic and transcriptomic data of agricultural plants. CropGene includes several blocks of bioinformatic analysis, such as analysis of gene variations, assembly of genomes and transcriptomes, as well as annotation of genes and proteins. CropGene implements new methods for analyzing long non-coding RNAs, protein domains, searching and analyzing polymorphisms, and genome-wide association research. CropGene has a user-friendly interface and supports working with various types of data, which greatly simplifies its use for researchers who do not have deep knowledge in the field of bioinformatics. The paper provides examples of the use of CropGene for the analysis of agricultural organisms such as Solanum tuberosum and Zea mays. With CropGene, genetic markers have been identified that explain up to 50 % of the variability in seed color parameters; potential genes that may become promising material for producing potato varieties; more than 100 thousand new long non-coding RNAs. Orthogroups were also found, the domain structure of which shows a marked similarity with the domain architecture of characteristic secreted A2 phospholipases. Thus, CropGene is an important tool for scientists and practitioners working in the field of agrobiotechnology and plant genetics.

Genomics Research on the Road of Studying Biology and Virulence of Cereal Rust Fungi

Rust fungi are highly destructive pathogens that pose a significant threat to crop production worldwide, especially cereals. Obligate biotrophy and, in many cases, complex life cycles make rust fungi particularly challenging to study. However, recent rapid advances in sequencing technologies and genomic analysis tools have revolutionised rust fungal research. It is anticipated that the increasing availability and ongoing substantial improvements in genome assemblies will propel the field of rust biology into the post-genomic era, instigating a cascade of research endeavours encompassing multi-omics and gene discoveries. This is especially the case for many cereal rust pathogens, for which continental-scale studies of virulence have been conducted over many years and historical collections of viable isolates have been sequenced and assembled. Genomic analysis plays a crucial role in uncovering the underlying causes of the high variability of virulence and the complexity of population dynamics in rust fungi. Here, we provide an overview of progress in rust genomics, discuss the strategies employed in genomic analysis, and elucidate the strides that will drive cereal rust biology into the post-genomic era.

Metronidazole response profiles of Gardnerella species are congruent with phylogenetic and comparative genomic analyses

BACKGROUND

Bacterial vaginosis (BV) affects 20-50% of reproductive-age female patients annually, arising when opportunistic pathogens outcompete healthy vaginal flora. Many patients fail to resolve symptoms with a course of metronidazole, the current first-line treatment for BV. Our study was designed to identify genomic variation associated with metronidazole resistance among strains of Gardnerella vaginalis spp. (GV), a genus of biogenic-amine-producing bacteria closely associated with BV pathogenesis, for the development of a companion molecular diagnostic.

METHODS

Whole-genome sequencing and comparative genomic metrics, including average nucleotide identity and GC content, were performed on a diverse set of 129 GV genomes to generate data for detailed taxonomic analyses. Pangenomic analyses were employed to construct a phylogenetic tree and cluster highly related strains within genospecies. G. vaginalis spp. clinical isolates within our collection were subjected to plate-based minimum inhibitory concentration (MIC) testing of metronidazole (n = 60) and clindamycin (n = 63). DECIPHER and MAFFT were used to identify genospecies-specific primers associated with antibiotic-resistance phenotypes. PCR-based analyses with these primers were used to confirm their specificity for the relevant genospecies.

RESULTS

Eleven distinct genospecies based on standard ANI criteria were identified among the GV strains in our collection. Metronidazole MIC testing revealed six genospecies within a closely related phylogenetic clade contained only highly metronidazole-resistant strains (MIC ≥ 32 µg/mL) and suggested at least two mechanisms of metronidazole resistance within the eleven GV genospecies. All strains within the six highly metronidazole-resistant genospecies displayed susceptibility to clinically relevant clindamycin concentrations (MIC ≤ 2 µg/mL). A PCR-based molecular diagnostic assay was developed to distinguish between members of the metronidazole-resistant and mixed-response genospecies, which should be useful for determining the clade membership of various GV strains and could assist in the selection of appropriate antibiotic therapies for BV cases.

CONCLUSIONS

This study provides comparative genomic and phylogenetic evidence for eleven distinct genospecies within the genus Gardnerella vaginalis spp., and identifies genospecies-specific responses to metronidazole, the first-line treatment for BV. A companion molecular diagnostic assay was developed that is capable of identifying essentially all highly metronidazole-resistant strains that phylogenetically cluster together within the GV genospecies, which is informative for antibiotic treatment options.

Machine learning uncovers the transcriptional regulatory network for the production host Streptomyces albidoflavus

Streptomyces albidoflavus is a widely used strain for natural product discovery and production through heterologous biosynthetic gene clusters (BGCs). However, the transcriptional regulatory network (TRN) and its impact on secondary metabolism remain poorly understood. Here, we characterize the TRN using independent component analysis on 218 RNA sequencing (RNA-seq) transcriptomes across 88 unique growth conditions. We identify 78 independently modulated sets of genes (iModulons) that quantitatively describe the TRN across diverse conditions. Our analyses reveal (1) TRN adaptation to different growth conditions, (2) conserved and unique characteristics of the TRN across diverse lineages, (3) transcriptional activation of several endogenous BGCs, including surugamide, minimycin, and paulomycin, and (4) inferred functions of 40% of uncharacterized genes in the S. albidoflavus genome. These findings provide a comprehensive and quantitative understanding of the S. albidoflavus TRN, offering a knowledge base for further exploration and experimental validation.

Comparative genomics reveals genetic diversity and differential metabolic potentials of the species of Arachnia and suggests reclassification of Arachnia propionica E10012 (=NBRC_14587) as novel species

The genus Arachnia, including Arachnia propionica and Arachnia rubra, are part of the normal oral and respiratory microbiota but can act as opportunistic pathogens in humans. This study investigates the functional, phylogenomic and taxonomic characteristics of 10 completely sequenced Arachnia strains, to elucidate their evolutionary relationships and divergence patterns, focusing on genomic variability and functional diversity. Phylogenetic analyses revealed distinct patterns, with Arachnia propionica strains showing significant divergence compared to the conserved Arachnia rubra strains. Notably, E10012 (=NBRC 14587) emerged as a distinct lineage with unique adaptations, while NCTC11666 exhibited a unique phylogenetic position, suggesting subspecies-level classification. Functional analyses highlighted variability among Arachnia propionica strains, with E10012 (=NBRC 14587) showing genes linked to choline metabolism and metal resistance, and NCTC11666 enriched in carbohydrate-active enzymes like GH179. In contrast, Arachnia rubra demonstrated genomic conservation, indicative of evolutionary specialization. This study reveals that strains E10012 (=NBRC 14587) and NCTC11666 displayed unique genomic features and distinct phylogenetic positioning, suggesting their reclassification as potential novel species and subspecies respectively. This underscores the balance between genomic conservation and diversification in Arachnia, reflecting their ecological adaptability and functional roles in the oral microbiome.

Genomic configuration of Bacillus subtilis (NMB01) unveils its antiviral activity against Orthotospovirus arachinecrosis infecting tomato

Orthotospovirus arachinecrosis (groundnut bud necrosis virus, GBNV) infecting tomato is a devastating viral pathogen responsible for severe yield losses of up to 100%. Considering the significance of the plant growth-promoting bacteria to induce innate immunity, attempts were made to evaluate the antiviral efficacy of Bacillus subtilis NMB01 against GBNV in cowpea and tomato. Foliar application of B. subtilis NMB01 at 1.5% onto the leaves of cowpea and tomato followed by challenge inoculation with GBNV significantly reduced the incidence of GBNV from 80% to 90% in response to the untreated inoculated control. Hence, we had a quest to understand if any genes were contributing toward the suppression of GBNV in assay hosts. To unveil the secrecy, whole-genome sequencing of B. subtilis NMB01 was carried out. The genome sequence of NMB01 revealed the presence of secondary metabolite biosynthetic gene clusters, including non-ribosomal peptide synthetases (NRPSs) and polyketide synthases (PKSs) which also encoded bacteriocins and antimicrobial peptides. The pan-genome analysis identified 1,640 core genes, 4,885 dispensable genes, and 60 unique genes, including MAMP genes that induce host immune responses. Comparative genome and proteome analysis with other genomes of B. subtilis strains in a public domain through OrthoVenn analysis revealed the presence of 4,241 proteins, 3,695 clusters, and 655 singletons in our study isolate. Furthermore, the NMB01-treated tomato plants increased the levels of defense-related genes (MAPKK1, WRKY33, PR1, PAL, and NPR1), enhancing immune system priming against GBNV infection. These findings suggest that B. subtilis NMB01 can be used as a promising biological control agent for managing plant viral disease sustainably.

Genomic analysis of multidrug-resistant Escherichia coli isolated from dairy cows in Shihezi city, Xinjiang, China

Introduction

Dairy farming plays a vital role in agriculture and nutrition; however, the emergence of antimicrobial resistance (AMR) among bacterial pathogens poses significant risks to public health and animal welfare. Multidrug-resistant (MDR) Escherichia coli strains are of particular concern due to their potential for zoonotic transmission and resistance to multiple antibiotics. In this study, we investigated the prevalence of AMR and analyzed the genomes of two MDR E. coli isolated from dairy cows in Shihezi City.

Methods

Fecal samples were collected from dairy cows, and E. coli strains were isolated. Antibiotic susceptibility testing was conducted using the Kirby-Bauer disk diffusion method against 14 antibiotics. Two MDR isolates (E.coli_30 and E.coli_45) were selected for whole-genome sequencing and comparative genomic analysis. The Comprehensive Antibiotic Resistance Database (CARD) was used to identify AMR genes, and virulence factors were analyzed. Phylogenetic analysis was performed to determine the evolutionary relationships of the isolates, and a pangenome analysis of 50 E. coli strains was conducted to assess genetic diversity. The presence of mobile genetic elements (MGEs), including insertion sequences (IS) and transposons, was also examined.

Results

Among the E. coli isolates, 22.9% exhibited MDR, with high resistance to imipenem and ciprofloxacin, while gentamicin and tetracycline remained the most effective antibiotics. Genomic analysis revealed key AMR genes, including mphA, qnrS1, and bla CTX-M-55 (the latter found only in E.coli_45), conferring resistance to macrolides, quinolones, and beta-lactams, respectively. Virulence genes encoding type III secretion systems (TTSS) and adhesion factors were identified, indicating pathogenic potential. Phylogenetic analysis showed that E.coli_30 and E.coli_45 originated from distinct ancestral lineages. The presence of two extended-spectrum β-lactamase (ESBL) genes in E.coli_45 was noticeable, so we studied their global and national distribution using evolutionary analysis. We found that they are endemic in E. coli, Salmonella enterica, and Klebsiella pneumoniae. Pangenome analysis revealed significant genetic diversity among E. coli strains, with unique genes related to metabolism and stress response. This indicates the bacteria's adaptation to various environments. MGEs were identified as key contributors to genetic variability and adaptation.

Discussion

This study highlights the growing threat of MDR E. coli in dairy farms, emphasizing the critical role of MGEs in the spread of resistance genes. The genetic diversity observed suggests strong adaptive capabilities, justifying the need for continuous AMR surveillance in livestock. Effective monitoring and mitigation strategies are essential to prevent the dissemination of MDR bacteria, thereby protecting both animal and public health.

Developing pangenomes for large and complex plant genomes and their representation formats

BACKGROUND

The development of pangenomes has revolutionized genomic studies by capturing the complete genetic diversity within a species. Pangenome assembly integrates data from multiple individuals to construct a comprehensive genomic landscape, revealing both core and accessory genomic elements. This approach enables the identification of novel genes, structural variations, and gene presence-absence variations, providing insights into species evolution, adaptation, and trait variation. Representing pangenomes requires innovative visualization formats that effectively convey the complex genomic structures and variations.

AIM

This review delves into contemporary methodologies and recent advancements in constructing pangenomes, particularly in plant genomes. It examines the structure of pangenome representation, including format comparison, conversion, visualization techniques, and their implications for enhancing crop improvement strategies.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Earlier comparative studies have illuminated novel gene sequences, copy number variations, and presence-absence variations across diverse crop species. The concept of a pan-genome, which captures multiple genetic variations from a broad spectrum of genotypes, offers a holistic perspective of a species' genetic makeup. However, constructing a pan-genome for plants with larger genomes poses challenges, including managing vast genome sequence data and comprehending the genetic variations within the germplasm. To address these challenges, researchers have explored cost-effective alternatives to encapsulate species diversity in a single assembly known as a pangenome. This involves reducing the volume of genome sequences while focusing on genetic variations. With the growing prominence of the pan-genome concept in plant genomics, several software tools have emerged to facilitate pangenome construction. This review sheds light on developing and utilizing software tools tailored for constructing pan-genomes in plants. It also discusses representation formats suitable for downstream analyses, offering valuable insights into the genetic landscape and evolutionary dynamics of plant species. In summary, this review underscores the significance of pan-genome construction and representation formats in resolving the genetic architecture of plants, particularly those with complex genomes. It provides a comprehensive overview of recent advancements, aiding in exploring and understanding plant genetic diversity.

Population genetic analysis of clinical Mycobacterium abscessus complex strains in China

Background

To explore the genetic characteristics of the Mycobacterium abscessus complex (MABC) population in China, given its rising clinical importance among nontuberculous mycobacteria.

Methods

We conducted population genetic analyses on 360 MABC genomes from China, focusing on core genome multilocus sequence typing (cgMLST), pan-genome characterization, population genetics, and antimicrobial resistance gene profiling.

Results

Our analysis identified 273 M. abscessus subsp. abscessus (MabA) and 87 M. abscessus subsp. massiliense (MabM) isolates, uncovering 68 sequence types (STs), with ST5 being the most common. cgMLST classified 33.3% of isolates into six dominant circulating clones (DCCs) and 49.4% into 59 genomic clusters at a threshold of 25 different alleles, including 18 international clusters linking Chinese isolates with seven other countries. The MABC pan-genome is open, with MabA exhibiting greater accessory gene diversity and higher gene turnover compared to MabM. Mobile genetic elements (MGEs), such as prophages and genomic islands, were prevalent across all genomes. 139 to 151 virulence factors (VFs) were identified per genome, with distinct accessory VFs in MabA and MabM affecting immune modulation and metabolism. Resistance gene profiling revealed ubiquitous mtrA, RbpA, and bla MAB, with MabA-specific erm(41) conferring resistance to macrolides and β-lactams. Common rrs and rrl gene mutations indicated widespread resistance to aminoglycosides and macrolides, while gyrA mutations suggested emerging fluoroquinolone resistance. An acquired erm(46) gene, likely obtained via phage-mediated horizontal gene transfer, was detected in one MabA strain.

Conclusion

This study provides key genetic insights into the dynamics of MABC in China. The widespread distribution of DCCs, high genomic clustering rates, open pan-genome, and distinct resistance patterns between MabA and MabM, along with MGEs, highlight the need for targeted surveillance and tailored therapies to address emerging challenges in MABC infections.

Pangenome and pantranscriptome as the new reference for gene-family characterization: A case study of basic helix-loop-helix (bHLH) genes in barley

Genome-wide identification and comparative gene-family analyses have commonly been performed to investigate species-specific evolution linked to various traits and molecular pathways. However, most previous studies have been limited to gene screening in a single reference genome, failing to account for the gene presence/absence variations (gPAVs) in a species. Here, we propose an innovative pangenome-based approach for gene-family analyses based on orthologous gene groups (OGGs). Using the basic helix-loop-helix (bHLH) transcription factor family in barley as an example, we identified 161-176 bHLHs in 20 barley genomes, which can be classified into 201 OGGs. These 201 OGGs were further classified into 140 core, 12 softcore, 29 shell, and 20 line-specific/cloud bHLHs, revealing the complete profile of bHLH genes in barley. Using a genome-scanning approach, we overcame the genome annotation bias and identified an average of 1.5 un-annotated core bHLHs per barley genome. We found that whole-genome/segmental duplicates are predominant mechanisms contributing to the expansion of most core/softcore bHLHs, whereas dispensable bHLHs are more likely to result from small-scale duplication events. Interestingly, we noticed that the dispensable bHLHs tend to be enriched in the specific subfamilies SF13, SF27, and SF28, implying the potentially biased expansion of specific bHLHs in barley. We found that 50% of the bHLHs contain at least 1 intact transposon element (TE) within the 2-kb upstream-to-downstream region. bHLHs with copy-number variations (CNVs) have 1.48 TEs on average, significantly more than core bHLHs without CNVs (1.36), supporting a potential role of TEs in bHLH expansion. Analyses of selection pressure showed that dispensable bHLHs have experienced clear relaxation of selection compared with core bHLHs, consistent with their conservation patterns. We also integrated the pangenome data with recently available barley pantranscriptome data from 5 tissues and discovered apparent transcriptional divergence within and across bHLH subfamilies. We conclude that pangenome-based gene-family analyses can better describe the previously untapped, genuine evolutionary status of bHLHs and provide novel insights into bHLH evolution in barley. We expect that this study will inspire similar analyses in many other gene families and species.

A sweeping view of avian mycoplasmas biology drawn from comparative genomic analyses

BACKGROUND

Avian mycoplasmas are small bacteria associated with several pathogenic conditions in many wild and poultry bird species. Extensive genomic data are available for many avian mycoplasmas, yet no comparative studies focusing on this group of mycoplasmas have been undertaken so far.

RESULTS

Here, based on the comparison of forty avian mycoplasma genomes belonging to ten different species, we provide insightful information on the phylogeny, pan/core genome, energetic metabolism, and virulence of these avian pathogens. Analyses disclosed considerable inter- and intra-species genomic variabilities, with genome sizes that can vary by twice as much. Phylogenetic analysis based on concatenated orthologous genes revealed that avian mycoplasmas fell into either Hominis or Pneumoniae groups within the Mollicutes and could split into various clusters. No host co-evolution of avian mycoplasmas can be inferred from the proposed phylogenetic scheme. With 3,237 different gene clusters, the avian mycoplasma group under study proved diverse enough to have an open pan genome. However, a set of 150 gene clusters was found to be shared between all avian mycoplasmas, which is likely encoding essential functions. Comparison of energy metabolism pathways showed that avian mycoplasmas rely on various sources of energy. Superposition between phylogenetic and energy metabolism groups revealed that the glycolytic mycoplasmas belong to two distinct phylogenetic groups (Hominis and Pneumoniae), while all the arginine-utilizing mycoplasmas belong only to Hominis group. This can stand for different evolutionary strategies followed by avian mycoplasmas and further emphasizes the diversity within this group. Virulence determinants survey showed that the involved gene arsenals vary significantly within and between species, and could even be found in species often reported apathogenic. Immunoglobulin-blocking proteins were detected in almost all avian mycoplasmas. Although these systems are not exclusive to this group, they seem to present some particular features making them unique among mycoplasmas.

CONCLUSION

This comparative genomic study uncovered the significant variable nature of avian mycoplasmas, furthering our knowledge on their biological attributes and evoking new hallmarks.

A pan-genome perspective on the evolutionary dynamics of polyphyly, virulence, and antibiotic resistance in Salmonella enterica serovar Mbandaka highlights emerging threats to public health and food safety posed by cloud gene families

Salmonella enterica serovar Mbandaka, a prevalent foodborne pathogen, poses a threat to public health but remains poorly understood. We have determined the phylogenomic tree, genetic diversity, virulence, and antimicrobial resistance (AMR) profiles on a large genomic scale to elucidate the evolutionary dynamics within the Mbandaka pan-genome. The polyphyletic nature of this serovar is characterized by two distinct phylogenetic groups and inter-serovar recombination boundaries, that potentially arising from recombination events at the H2-antigen loci. The open pan-genome exhibited a flexible gene repertoire, with numerous cloud gene families involved in virulence and AMR. Extensive gene gain and loss observed at the terminal nodes of the phylogenetic tree indicate that Mbandaka individuals have undergone frequent gene turnover. The resulting changes in virulence and AMR genes potentially pose emerging threats to public health. We explored serovar conversion due to recombination of H-antigen loci, inter-serovar divergences in gene gain and loss, prophage-mediated acquisition of virulence factors, and the role of incompatibility group plasmids in acquiring resistance determinants as key molecular mechanisms driving the pathogenicity and antibiotic resistance of Mbandaka. Our work contributes to a comprehensive understanding of the complex mechanisms of pathogenesis and the ongoing evolutionary arms race with current therapeutic approaches in serovar Mbandaka.

Genomic and metabolomic analysis of Latilactobacillus sakei DCF0720 for black soybean yogurt fermentation

Lactic acid bacteria are commonly used in plant-based fermentation to reduce off-flavor and improve sensory characteristics. However, there have been few studies on Latilactobacillus sakei for plant-based yogurt fermentation and, particularly, its metabolic features at the genomic level remain unclear. This study aims to analyze the fermentation characteristics of the L. sakei DCF0720 strain and compare genetics and metabolic relations. For this, DCF0720 was used to ferment the black soybean milk and conduct the physicochemical analysis and sensory test. The genomic and metabolic analyses were performed by complete genome sequencing and 500 MHz 1H NMR, respectively. As a result, DCF0720 exhibited enhanced fermentation performance and sensory evaluations at 37 °C compared to 30 °C, which is generally recognized as the optimal growth temperature for most L. sakei strains. It also produced flavor enhancing volatile compounds such as acetoin and hydroxyacetone, possessing all three key genes for acetoin biosynthesis. DCF0720 lacks 2,3-butanediol dehydrogenase, which leads to the inhibition of acetoin production. DCF0720 possesses a complete pathway to utilize primary black soybean carbon sources such as sucrose, raffinose, and stachyose. DCF0720 also possesses genes for the GH28 family, including the key enzymes in the hydrolysis of pectin substances, which means eliminating the main soybean nonstarch polysaccharides. This study demonstrates that DCF0720 is a suitable starter for plant-based yogurt fermentation, providing a better understanding of fermentation conditions with genetic and metabolic features for black soybean yogurt. Various carbon source utilization abilities with depth metabolic pathway analysis provide that DCF0720 can be employed to develop enhanced starter cultures for black soybean yogurt and diverse plant-based yogurts.

Pangenome analysis of Paenibacillus polymyxa strains reveals the existence of multiple and functionally distinct Paenibacillus species

Paenibacillus polymyxa, a Gram-positive bacterium commonly found in soil and plant roots, plays an important role in the environment due to its nitrogen-fixing ability and is renowned for producing antibiotics like polymyxin. In this study, we present a robust framework for investigating the evolutionary and taxonomic connections of strains belonging to P. polymyxa available at the National Center for Biotechnology Information, as well as five new additional strains isolated at the University of Camerino (Italy), through pangenome analysis. These strains can produce secondary metabolites active against Staphylococcus aureus and Klebsiella pneumoniae. Employing techniques such as digital DNA-DNA hybridization (dDDH), average nucleotide identity (ANI) estimation, OrthoFinder, and ribosomal multilocus sequence typing, we consistently divided these P. polymyxa strains into four clusters, which differ significantly in terms of ANI and dDDH percentages, both considered as reference indices for separating bacterial species. Moreover, the strains of Cluster 2 were re-classified as belonging to the Paenibacillus ottowii species. By comparing the pangenomes, we identified the core genes of each cluster and analyzed them to recognize distinctive features in terms of biosynthetic/metabolic potential. The comparison of pangenomes also allowed us to pinpoint differences between clusters in terms of genetic variability and the percentage of the genome dedicated to core and accessory genes. In conclusion, the data obtained from our analyses of strains belonging to the P. polymyxa species converge toward a necessary reclassification, which will require a fundamental contribution from microbiologists in the near future.

IMPORTANCE

The development of sequencing technologies has led to an exponential increase in microbial sequencing data. Accurately identifying bacterial species remains a challenge because of extensive intra-species variability, the need for multiple identification methods, and the rapid rate of taxonomic changes. A substantial contribution to elucidating the relationships among related bacterial strains comes from comparing their genomic sequences. This comparison also allows for the identification of the "pangenome," which is the set of genes shared by all individuals of a species, as well as the set of genes that are unique to subpopulations. Here, we applied this approach to Paenibacillus polymyxa, a species studied for its potential as a biofertilizer and biocontrol agent and known as an antibiotic producer. Our work highlights the need for a more efficient classification of this bacterial species and provides a better delineation of strains with different properties.

Comprehensive pan-genome analysis of Mycobacterium marinum: insights into genomic diversity, evolution, and pathogenicity

Mycobacteria is a diverse genus that includes both innocuous environmental species and serious pathogens like Mycobacterium tuberculosis, Mycobacterium leprae, and Mycobacterium ulcerans, the causative agents of tuberculosis, leprosy, and Buruli ulcer, respectively. This study focuses on Mycobacterium marinum, a closely related species known for its larger genome and ability to infect ectothermic species and cooler human extremities. Utilizing whole-genome sequencing, we conducted a comprehensive pan-genome analysis of 100 M. marinum strains, exploring genetic diversity and its impact on pathogenesis and host specificity. Our findings highlight significant genomic diversity, with clear distinctions in core, dispensable, and unique genes among the isolates. Phylogenetic analysis revealed a broad distribution of genetic lineages, challenging previous classifications into distinct clusters. Additionally, we examined the synteny and diversity of the virulence factor CpnT, noting a wide range of C-terminal domain variations across strains, which points to potential adaptations in pathogenic mechanisms. This study enhances our understanding of M. marinum's genomic architecture and its evolutionary relationship with other mycobacterial pathogens, providing insights that could inform disease control strategies for M. tuberculosis and other mycobacteria.

Temporal changes in the mouse gut bacteriota influenced by host sex, diet, and exercise

The gut microbiota plays an important role in host physiology. However, the effects of host sex, lifestyle, and temporal influences on the bacterial community within the gut remain ill-defined. To address this gap, we evaluated 56 male and female mice over a 10-wk study to assess the effects of sex, diet, and exercise on gut community dynamics. Mice were randomly assigned to high-fat or control diet feeding and had free access to running wheels or remained sedentary throughout the study period. The fecal bacterial community was characterized by rRNA operon profiling via nanopore sequencing. Differential abundance testing indicated that ∼200 bacterial taxa were significantly influenced by sex, diet, or exercise (4.2% of total community), which also changed over time (82 taxa, 1.7% of total community). Phylogenetic analysis of taxa closely related to Dysosmobacter welbionis, and several members of the family Muribaculaceae were examined more closely and demonstrated distinct species/strain level subclustering by host sex, diet, and exercise. Collectively, these data suggest that sex and lifestyle can alter the gut bacteriota at the species/strain level that may play a role in host health. These results also highlight the need for improved characterization methods to survey microbial communities at finer taxonomic resolution.NEW & NOTEWORTHY This study demonstrates that host sex, diet, and exercise can each modulate gut bacterial community structure, which may have consequences to host physiology. Our analysis shows selection of novel strains and genera for some members of the Oscillospiraceae and Muribaculaceae by host sex, diet, and activity status. Overall, these findings provide a framework for detecting the next generation of beneficial bacteria targeting obesity and associated metabolic diseases in a sex-specific manner.

When less is more: sketching with minimizers in genomics

The exponential increase in sequencing data calls for conceptual and computational advances to extract useful biological insights. One such advance, minimizers, allows for reducing the quantity of data handled while maintaining some of its key properties. We provide a basic introduction to minimizers, cover recent methodological developments, and review the diverse applications of minimizers to analyze genomic data, including de novo genome assembly, metagenomics, read alignment, read correction, and pangenomes. We also touch on alternative data sketching techniques including universal hitting sets, syncmers, or strobemers. Minimizers and their alternatives have rapidly become indispensable tools for handling vast amounts of data.

Chromosomal gene order defines several structural classes of Staphylococcus epidermidis genomes

The original methodology for describing the pangenome of a prokaryotic species is based on modeling genomes as unordered sets of genes. More recent findings have underlined the importance of considering the ordering of genes along the genetic material as well, when making comparisons among genomes. To further investigate the benefits of gene order when describing genomes of a given species, we applied two distance metrics on a dataset of 84 genomes of Staphylococcus epidermidis. The first metric, GeLev, depends on the order of genes and is a derivative of the Levenshtein distance. The second, the Jaccard distance, depends on gene sets only. The application of these distances reveals information about the global structure of the genomes, and allows clustering of the genomes into classes. The main biological result is that, while genomes within the same class are structurally similar, genomes of different classes have an additional characteristic. Between genomes in different classes we can discover instances where a large segment of the first genome appears in reverse order in the second. This feature suggests that genome rearrangements in S. epidermidis happen on a large scale, while micro-rearrangements of single or a small number of genes are rare. Thus, this paper describes a straight-forward method to classify genomes into structural classes with the same order of genes and makes it possible to visualize reversed segments in pairs of genomes. The method can be readily applied to other species.

A gentle introduction to pangenomics

Pangenomes have emerged in response to limitations associated with traditional linear reference genomes. In contrast to a traditional reference that is (usually) assembled from a single individual, pangenomes aim to represent all of the genomic variation found in a group of organisms. The term 'pangenome' is currently used to describe multiple different types of genomic information, and limited language is available to differentiate between them. This is frustrating for researchers working in the field and confusing for researchers new to the field. Here, we provide an introduction to pangenomics relevant to both prokaryotic and eukaryotic organisms and propose a formalization of the language used to describe pangenomes (see the Glossary) to improve the specificity of discussion in the field.

Deciphering the gastrointestinal carriage of Klebsiella pneumoniae

Bacterial infections pose a significant global health threat, accounting for an estimated 7.7 million deaths. Hospital outbreaks driven by multi-drug-resistant pathogens, notably Klebsiella pneumoniae (K. pneumoniae), are of grave concern. This opportunistic pathogen causes pneumonia, urinary tract infections, and bacteremia, particularly in immunocompromised individuals. The rise of hypervirulent K. pneumoniae adds complexity, as it increasingly infects healthy individuals. Recent epidemiological data suggest that asymptomatic gastrointestinal carriage serves as a reservoir for infections in the same individual and allows for host-to-host transmission via the fecal-oral route. This review focuses on K. pneumoniae's gastrointestinal colonization, delving into epidemiological evidence, current animal models, molecular colonization mechanisms, and the protective role of the resident gut microbiota. Moreover, the review sheds light on in vivo high-throughput approaches that have been crucial for identifying K. pneumoniae factors in gut colonization. This comprehensive exploration aims to enhance our understanding of K. pneumoniae gut pathogenesis, guiding future intervention and prevention strategies.

Whole Genome Sequence Analysis of Brucella spp. from Human, Livestock, and Wildlife in South Africa

Brucellosis is an economically important zoonotic disease affecting humans, livestock, and wildlife health globally and especially in Africa. Brucella abortus and B. melitensis have been isolated from human, livestock (cattle and goat), and wildlife (sable) in South Africa (SA) but with little knowledge of the population genomic structure of this pathogen in SA. As whole genome sequencing can assist to differentiate and trace the origin of outbreaks of Brucella spp. strains, the whole genomes of retrospective isolates (n = 19) from previous studies were sequenced. Sequences were analysed using average nucleotide identity (ANI), pangenomics, and whole genome single nucleotide polymorphism (wgSNP) to trace the geographical origin of cases of brucellosis circulating in human, cattle, goats, and sable from different provinces in SA. Pangenomics analysis of B. melitensis (n = 69) and B. abortus (n = 56) was conducted with 19 strains that included B. abortus from cattle (n = 3) and B. melitensis from a human (n = 1), cattle (n = 1), goat (n = 1), Rev1 vaccine strain (n = 1), and sable (n = 12). Pangenomics analysis of B. melitensis genomes, highlighted shared genes, that include 10 hypothetical proteins and genes that encodes for acetyl-coenzyme A synthetase (acs), and acylamidase (aam) amongst the sable genomes. The wgSNP analysis confirmed the B. melitensis isolated from human was more closely related to the goat from the Western Cape Province from the same outbreak than the B. melitensis cattle sample from different cases in the Gauteng Province. The B. melitensis sable strains could be distinguished from the African lineage, constituting their own African sub-clade. The sequenced B. abortus strains clustered in the C2 lineage that is closely related to the isolates from Mozambique and Zimbabwe. This study identified genetically diverse Brucella spp. among various hosts in SA. This study expands the limited known knowledge regarding the presence of B. melitensis in livestock and humans in SA, further building a foundation for future research on the distribution of the Brucella spp. worldwide and its evolutionary background.

Unveiling a Comprehensive Multi-epitope Subunit Vaccine Strategy Against Salmonella subsp. enterica: Bridging Core, Subtractive Proteomics, and Immunoinformatics

Salmonella subsp. enterica (SE) presents a significant global health challenge in both developed and developing countries. Current SE vaccines have limitations, targeting specific strains and demonstrating moderate efficacy in adults, while also being unsuitable for young children and often unaffordable in regions with lower income levels where the disease is prevalent. To address these challenges, this study employed a computational approach integrating core proteomics, subtractive proteomics, and immunoinformatics to develop a universal SE vaccine and identify potential drug targets. Analysis of the core proteome of 185 SE strains revealed 1964 conserved proteins. Subtractive proteomics identified 9 proteins as potential vaccine candidates and 41 as novel drug targets. Using reverse vaccinology-based immunoinformatics, four multi-epitope-based subunit vaccine constructs (MESVCs) were designed, aiming to stimulate cytotoxic T lymphocyte, helper T lymphocyte, and linear B lymphocyte responses. These constructs underwent comprehensive evaluations for antigenicity, immunogenicity, toxicity, hydropathicity, and physicochemical properties. Predictive modeling, refinement, and validation were conducted to determine the secondary and tertiary structures of the SE-MESVCs, followed by docking studies with MHC-I, MHC-II, and TLR4 receptors. Molecular docking assessments showed favorable binding with all three receptors, with SE-MESVC-4 exhibiting the most promising binding energy. Molecular dynamics simulations confirmed the binding affinity and stability of SE-MESVC-4 with the TLR4/MD2 complex. Additionally, codon optimization and in silico cloning verified the efficient translation and successful expression of SE-MESVC-4 in Escherichia coli (E. coli) str. K12. Subsequent in silico immune simulation evaluated the efficacy of SE-MESVC-4 in triggering an effective immune response. These results suggest that SE-MESVC-4 may induce both humoral and cellular immune responses, making it a potential candidate for an effective SE vaccine. However, further experimental investigations are necessary to validate the immunogenicity and efficacy of SE-MESVC-4, bringing us closer to effectively combating SE infections.

Characterization of the prevalence of Salmonella in different retail chicken supply modes using genome-wide and machine-learning analyses

Salmonella is a foodborne pathogen that causes salmonellosis, of which retail chicken meat is a major source. However, the prevalence of Salmonella in different retail chicken supply modes and the threat posed to consumers remains unclear. The prevalence, serotype distribution, antibiotic resistance, and genomic characteristics of Salmonella in three supply modes of retail chicken (live poultry, frozen, and chilled) were investigated using whole-genome sequencing (WGS) and machine learning (ML). In this study, 480 retail chicken samples from live poultry, frozen, and chilled supply modes in Guangzhou from 2020 to 2021, as well as 253 Salmonella isolates (total isolation rate = 53.1 %), were collected. The prevalence of isolates in the live poultry mode (67.5 %, 81/120) was statistically higher than in the frozen (50.0 %, 120/240) and chilled (43.3 %, 52/120) (P < 0.05) modes. Serotype identification showed significant differences in the serotype distribution of Salmonella in different supply modes. S. Enteritis (46.7 %) and S. Indiana (14.2 %) were predominant in the frozen mode. S. Agona (23.5 %) and S. Saintpaul (13.6 %) were predominant in live poultry, while S. Enteritis (40.4 %) and S. Kentucky (17.3 %) were predominant in chilled mode. Antibiotic testing showed that frozen mode isolates were more resistant; the multidrug-resistant (MDR) rate of isolates in the frozen mode reached 91.8 %, significantly higher than in the chilled (86.5 %) and live (74.1 %) (P < 0.05) modes. WGS was performed on 155 top serotypes (S. Enteritidis, S. Kentucky, S. Indiana, and S. Agona). The antibiotic resistance gene analysis showed that the abundance and carrying rate of antibiotic resistance genes of Salmonella in the frozen mode (54 types, 16.1 %) were significantly higher than in other modes (live poultry: 36 types, 9.4 %, P < 0.05; chilled: 31 types, 11.6 %). The blaNDM-1 and blaNDM-9 genes encoding carbapenem resistance were found in frozen mode isolates on a complex transposon consisting of TnAS3-IS26. Virulence factors and plasmid replicons were abundant in the studied frozen mode isolates. In addition, single nucleotide polymorphism (SNP) phylogenetic tree results showed that in the frozen supply mode, the S. Enteritidis clonal clade continued to contaminate retail chicken meat and was homologous to S. Enteritidis strains found in farm chicken embryos, slaughterhouse chicken carcasses, and patients from hospitals in China (SNP 0 = 10). Notably, the pan-genome-based ML model showed that characteristic genes in frozen and live poultry isolates differed. The narZ gene was a key characteristic gene in frozen isolates, encoding nitrate reductase, relating to anaerobic bacterial growth. The ydgJ gene is a key characteristic gene in the live mode and encodes an oxidoreductase related to oxidative function in bacteria. The high prevalence of live poultry mode Salmonella and the transmission of frozen mode MDR Salmonella in this study pose serious risks to food safety and public health, emphasizing the importance of improving disinfection and cold storage measures to reduce Salmonella contamination and transmission. In conclusion, the continued surveillance of Salmonella across different supply models and the development of an epidemiological surveillance system based on WGS is necessary.

Genome and pan-genome analysis of a new exopolysaccharide-producing bacterium Pyschrobacillus sp. isolated from iron ores deposit and insights into iron uptake

Bacterial exopolysaccharides (EPS) have emerged as one of the key players in the field of heavy metal-contaminated environmental bioremediation. This study aimed to characterize and evaluate the metal biosorption potential of EPS produced by a novel Psychrobacillus strain, NEAU-3TGS, isolated from an iron ore deposit at Tamra iron mine, northern Tunisia. Genomic and pan-genomic analysis of NEAU-3TGS bacterium with nine validated published Psychrobacillus species was also performed. The results showed that the NEAU-3TGS genome (4.48 Mb) had a mean GC content of 36%, 4,243 coding sequences and 14 RNA genes. Phylogenomic analysis and calculation of nucleotide identity (ANI) values (less than 95% for new species with all strains) confirmed that NEAU-3TGS represents a potential new species. Pangenomic analysis revealed that Psychrobacillus genomic diversity represents an "open" pangenome model with 33,091 homologous genes, including 65 core, 3,738 shell, and 29,288 cloud genes. Structural EPS characterization by attenuated total reflectance-Fourier transform infrared (ATR-FTIR) spectroscopy showed uronic acid and α-1,4-glycosidic bonds as dominant components of the EPS. X-ray diffraction (XRD) analysis revealed the presence of chitin, chitosan, and calcite CaCO3 and confirmed the amorphous nature of the EPS. Heavy metal bioabsorption assessment showed that iron and lead were more adsorbed than copper and cadmium. Notably, the optimum activity was observed at 37°C, pH=7 and after 3 h contact of EPS with each metal. Genomic insights on iron acquisition and metabolism in Psychrobacillus sp. NEAU-3TGS suggested that no genes involved in siderophore biosynthesis were found, and only the gene cluster FeuABCD and trilactone hydrolase genes involved in the uptake of siderophores, iron transporter and exporter are present. Molecular modelling and docking of FeuA (protein peptidoglycan siderophore-binding protein) and siderophores ferrienterobactine [Fe+3 (ENT)]-3 and ferribacillibactine [Fe+3 (BB)]-3 ligand revealed that [Fe+3 (ENT)]-3 binds to Phe122, Lys127, Ile100, Gln314, Arg215, Arg217, and Gln252. Almost the same for [Fe+3 (ENT)]-3 in addition to Cys222 and Tyr229, but not Ile100.To the best of our knowledge, this is the first report on the characterization of EPS and the adsorption of heavy metals by Psychrobacillus species. The heavy metal removal capabilities may be advantageous for using these organisms in metal remediation.

Beyond the Human Genome Project: The Age of Complete Human Genome Sequences and Pangenome References

The Human Genome Project was an enormous accomplishment, providing a foundation for countless explorations into the genetics and genomics of the human species. Yet for many years, the human genome reference sequence remained incomplete and lacked representation of human genetic diversity. Recently, two major advances have emerged to address these shortcomings: complete gap-free human genome sequences, such as the one developed by the Telomere-to-Telomere Consortium, and high-quality pangenomes, such as the one developed by the Human Pangenome Reference Consortium. Facilitated by advances in long-read DNA sequencing and genome assembly algorithms, complete human genome sequences resolve regions that have been historically difficult to sequence, including centromeres, telomeres, and segmental duplications. In parallel, pangenomes capture the extensive genetic diversity across populations worldwide. Together, these advances usher in a new era of genomics research, enhancing the accuracy of genomic analysis, paving the path for precision medicine, and contributing to deeper insights into human biology.

Pangenome and genomic signatures linked to the dominance of the lineage-4 of Mycobacterium tuberculosis isolated from extrapulmonary tuberculosis patients in western Ethiopia

BACKGROUND

The lineage 4 (L4) of Mycobacterium tuberculosis (MTB) is not only globally prevalent but also locally dominant, surpassing other lineages, with lineage 2 (L2) following in prevalence. Despite its widespread occurrence, factors influencing the expansion of L4 and its sub-lineages remain poorly understood both at local and global levels. Therefore, this study aimed to conduct a pan-genome and identify genomic signatures linked to the elevated prevalence of L4 sublineages among extrapulmonary TB (EPTB) patients in western Ethiopia.

METHODS

A cross-sectional study was conducted at an institutional level involving confirmed cases of extrapulmonary tuberculosis (EPTB) patients from August 5, 2018, to December 30, 2019. A total of 75 MTB genomes, classified under lineage 4 (L4), were used for conducting pan-genome and genome-wide association study (GWAS) analyses. After a quality check, variants were identified using MTBseq, and genomes were de novo assembled using SPAdes. Gene prediction and annotation were performed using Prokka. The pan-genome was constructed using GET_HOMOLOGUES, and its functional analysis was carried out with the Bacterial Pan-Genome Analysis tool (BPGA). For GWAS analysis, Scoary was employed with Benjamini-Hochberg correction, with a significance threshold set at p-value ≤ 0.05.

RESULTS

The analysis revealed a total of 3,270 core genes, predominantly associated with orthologous groups (COG) functions, notably in the categories of '[R] General function prediction only' and '[I] Lipid transport and metabolism'. Conversely, functions related to '[N] Cell motility' and '[Q] Secondary metabolites biosynthesis, transport, and catabolism' were primarily linked to unique and accessory genes. The pan-genome of MTB L4 was found to be open. Furthermore, the GWAS study identified genomic signatures linked to the prevalence of sublineages L4.6.3 and L4.2.2.2.

CONCLUSIONS

Apart from host and environmental factors, the sublineage of L4 employs distinct virulence factors for successful dissemination in western Ethiopia. Given that the functions of these newly identified genes are not well understood, it is advisable to experimentally validate their roles, particularly in the successful transmission of specific L4 sublineages over others.

Average nucleotide identity-based Staphylococcus aureus strain grouping allows identification of strain-specific genes in the pangenome

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

IMPORTANCE

We analyzed the genomic diversity of Staphylococcus aureus, a globally prevalent bacterial species that causes serious infections in humans. Our goal was to build a genetic picture of the different strains of S. aureus and which genes may be associated with them. We reprocessed >84,000 genomes and subsampled to remove redundancy. We found that individual samples sharing >99.5% of their genome could be grouped into strains. We also showed that a portion of genes that are present in intermediate frequency in the species are strongly associated with some strains but completely absent from others, suggesting a role in strain specificity. This work lays the foundation for understanding individual gene histories of the S. aureus species and also outlines strategies for processing large bacterial genomic data sets.

Progress on Salt Tolerance in Brassica napus

In China, saline-alkali lands constitute 5.01% of the total land area, having a significant impact on both domestic and international food production. Rapeseed (Brassica napus L.), as one of the most important oilseed crops in China, has garnered considerable attention due to its potential adaptability to saline conditions. Breeding and improving salt-tolerant varieties is a key strategy for the effective utilization of saline lands. Hence, it is important to conduct comprehensive research into the adaptability and salt tolerance mechanisms of Brassica napus in saline environments as well as to breed novel salt-tolerant varieties. This review summarizes the molecular mechanism of salt tolerance, physiological and phenotypic indexes, research strategies for the screening of salt-tolerant germplasm resources, and genetic engineering tools for salt stress in Brassica napus. It also introduces various agronomic strategies for applying exogenous substances to alleviate salt stress and provide technological tools and research directions for future research on salt tolerance in Brassica napus.

Pan-genomics: Insight into the Functional Genome, Applications, Advancements, and Challenges

A pan-genome is a compilation of the common and unique genomes found in a given species. It incorporates the genetic information from all of the genomes sampled, producing a big and diverse set of genetic material. Pan-genomic analysis has various advantages over typical genomics research. It creates a vast and varied spectrum of genetic material by combining the genetic data from all the sampled genomes. Comparing pan-genomics analysis to conventional genomic research, there are a number of benefits. Although the most recent era of pan-genomic studies has used cutting-edge sequencing technology to shed fresh light on biological variety and improvement, the potential uses of pan-genomics in improvement have not yet been fully realized. Pan-genome research in various organisms has demonstrated that missing genetic components and the detection of significant Structural Variants (SVs) can be investigated using pan-genomic methods. Many individual-specific sequences have been linked to biological adaptability, phenotypic, and key economic attributes. This study aims to focus on how pangenome analysis uncovers genetic differences in various organisms, including human, and their effects on phenotypes, as well as how this might help us comprehend the diversity of species. The review also concentrated on potential problems and the prospects for future pangenome research.

Comparative genomic analysis of Bradyrhizobium strains with natural variability in the efficiency of nitrogen fixation, competitiveness, and adaptation to stressful edaphoclimatic conditions

Bradyrhizobium is known for fixing atmospheric nitrogen in symbiosis with agronomically important crops. This study focused on two groups of strains, each containing eight natural variants of the parental strains, Bradyrhizobium japonicum SEMIA 586 (=CNPSo 17) or Bradyrhizobium diazoefficiens SEMIA 566 (=CNPSo 10). CNPSo 17 and CNPSo 10 were used as commercial inoculants for soybean crops in Brazil at the beginning of the crop expansion in the southern region in the 1960s-1970s. Variants derived from these parental strains were obtained in the late 1980s through a strain selection program aimed at identifying elite strains adapted to a new cropping frontier in the central-western Cerrado region, with a higher capacity of biological nitrogen fixation (BNF) and competitiveness. Here, we aimed to detect genetic variations possibly related to BNF, competitiveness for nodule occupancy, and adaptation to the stressful conditions of the Brazilian Cerrado soils. High-quality genome assemblies were produced for all strains. The core genome phylogeny revealed that strains of each group are closely related, as confirmed by high average nucleotide identity values. However, variants accumulated divergences resulting from horizontal gene transfer, genomic rearrangements, and nucleotide polymorphisms. The B. japonicum group presented a larger pangenome and a higher number of nucleotide polymorphisms than the B. diazoefficiens group, possibly due to its longer adaptation time to the Cerrado soil. Interestingly, five strains of the B. japonicum group carry two plasmids. The genetic variability found in both groups is discussed considering the observed differences in their BNF capacity, competitiveness for nodule occupancy, and environmental adaptation.IMPORTANCEToday, Brazil is a global leader in the study and use of biological nitrogen fixation with soybean crops. As Brazilian soils are naturally void of soybean-compatible bradyrhizobia, strain selection programs were established, starting with foreign isolates. Selection searched for adaptation to the local edaphoclimatic conditions, higher efficiency of nitrogen fixation, and strong competitiveness for nodule occupancy. We analyzed the genomes of two parental strains of Bradyrhizobium japonicum and Bradyrhizobium diazoefficiens and eight variant strains derived from each parental strain. We detected two plasmids in five strains and several genetic differences that might be related to adaptation to the stressful conditions of the soils of the Brazilian Cerrado biome. We also detected genetic variations in specific regions that may impact symbiotic nitrogen fixation. Our analysis contributes to new insights into the evolution of Bradyrhizobium, and some of the identified differences may be applied as genetic markers to assist strain selection programs.

Are Viruses Taxonomic Units? A Protein Domain and Loop-Centric Phylogenomic Assessment

Virus taxonomy uses a Linnaean-like subsumption hierarchy to classify viruses into taxonomic units at species and higher rank levels. Virus species are considered monophyletic groups of mobile genetic elements (MGEs) often delimited by the phylogenetic analysis of aligned genomic or metagenomic sequences. Taxonomic units are assumed to be independent organizational, functional and evolutionary units that follow a 'natural history' rationale. Here, I use phylogenomic and other arguments to show that viruses are not self-standing genetically-driven systems acting as evolutionary units. Instead, they are crucial components of holobionts, which are units of biological organization that dynamically integrate the genetics, epigenetic, physiological and functional properties of their co-evolving members. Remarkably, phylogenomic analyses show that viruses share protein domains and loops with cells throughout history via massive processes of reticulate evolution, helping spread evolutionary innovations across a wider taxonomic spectrum. Thus, viruses are not merely MGEs or microbes. Instead, their genomes and proteomes conduct cellularly integrated processes akin to those cataloged by the GO Consortium. This prompts the generation of compositional hierarchies that replace the 'is-a-kind-of' by a 'is-a-part-of' logic to better describe the mereology of integrated cellular and viral makeup. My analysis demands a new paradigm that integrates virus taxonomy into a modern evolutionarily centered taxonomy of organisms.

Pan-genome insights into adaptive evolution of bacterial symbionts in mixed host-microbe symbioses represented by human gut microbiota Bacteroides cellulosilyticus

Animal hosts harbor diverse assemblages of microbial symbionts that play crucial roles in the host's lifestyle. The link between microbial symbiosis and host development remains poorly understood. In particular, little is known about the adaptive evolution of gut bacteria in host-microbe symbioses. Recently, symbiotic relationships have been categorized as open, closed, or mixed, reflecting their modes of inter-host transmission and resulting in distinct genomic features. Members of the genus Bacteroides are the most abundant human gut microbiota and possess both probiotic and pathogenic potential, providing an excellent model for studying pan-genome evolution in symbiotic systems. Here, we determined the complete genome of an novel clinical strain PL2022, which was isolated from a blood sample and performed pan-genome analyses on a representative set of Bacteroides cellulosilyticus strains to quantify the influence of the symbiotic relationship on the evolutionary dynamics. B. cellulosilyticus exhibited correlated genomic features with both open and closed symbioses, suggesting a mixed symbiosis. An open pan-genome is characterized by abundant accessory gene families, potential horizontal gene transfer (HGT), and diverse mobile genetic elements (MGEs), indicating an innovative gene pool, mainly associated with genomic islands and plasmids. However, massive parallel gene loss, weak purifying selection, and accumulation of positively selected mutations were the main drivers of genome reduction in B. cellulosilyticus. Metagenomic read recruitment analyses showed that B. cellulosilyticus members are globally distributed and active in human gut habitats, in line with predominant vertical transmission in the human gut. However, existence and/or high abundance were also detected in non-intestinal tissues, other animal hosts, and non-host environments, indicating occasional horizontal transmission to new niches, thereby creating arenas for the acquisition of novel genes. This case study of adaptive evolution under a mixed host-microbe symbiosis advances our understanding of symbiotic pan-genome evolution. Our results highlight the complexity of genetic evolution in this unusual intestinal symbiont.

Globin phylogeny, evolution and function, the newest update

Our globin census update allows us to refine our vision of globin origin, evolution, and structure to function relationship in the context of the currently accepted tree of life. The modern globin domain originates as a single domain, three-over-three α-helical folded structure before the diversification of the kingdoms of life (Bacteria, Archaea, Eukarya). Together with the diversification of prokaryotes, three monophyletic globin families (M, S, and T) emerged, most likely in Proteobacteria and Actinobacteria, displaying specific sequence and structural features, and spread by vertical and horizontal gene transfer, most probably already present in the last universal common ancestor (LUCA). Non-globin domains were added, and eventually lost again, creating multi-domain structures in key branches of M- (FHb and Adgb) and the vast majority of S globins, which with their coevolved multi-domain architectures, have predominantly "sensor" functions. Single domain T-family globins diverged into four major groups and most likely display functions related to reactive nitrogen and oxygen species (RNOS) chemistry, as well as oxygen storage/transport which drives the evolution of its major branches with their characteristic key distal residues (B10, E11, E7, and G8). M-family evolution also lead to distinctive major types (FHb and Fgb, Ngb, Adgb, GbX vertebrate Gbs), and shows the shift from high oxygen affinity controlled by TyrB10-Gln/AsnE11 likely related to RNOS chemistry in microorganisms, to a moderate oxygen affinity storage/transport function controlled by hydrophobic B10/E11-HisE7 in multicellular animals.

Leaf morphology related genes revealed by integrating Pan-transcriptome, GWAS and eQTL analyses in a Liriodendron population

Leaf plays an indispensable role in plant development and growth. Although many known genes related to leaf morphology development have been identified, elucidating the complex genetic basis of leaf morphological traits remains a challenge. Liriodendron plants are common ornamental trees due to their unique leaf shapes, while the molecular mechanism underlying Liriodendron leaf morphogenesis has remained unknown. Herein, we firstly constructed a population-level pan-transcriptome of Liriodendron from 81 accessions to explore the expression presence or absence variations (ePAVs), global expression differences at the population level, as well as differentially expressed genes (DEGs) between the Liriodendron chinense and Liriodendron tulipifera accessions. Subsequently, we integrated a genome-wide association study (GWAS), expression quantitative trait loci (eQTL), and transcriptome-wide association study (TWAS) to identify candidate genes related to leaf morphology. Through GWAS analysis, we identified 18 and 17 significant allelic loci in the leaf size and leaf shape modules, respectively. In addition, we discerned 16 candidate genes in relation to leaf morphological traits via TWAS. Further, integrating the co-localization results of GWAS and eQTL, we determined two regulatory hotspot regions, hot88 and hot758, related to leaf size and leaf shape, respectively. Finally, co-expression analysis, eQTL, and linkage mapping together demonstrated that Lchi_4g10795 regulate their own expression levels through cis-eQTL to affect the expression of downstream genes and cooperatively participate in the development of Liriodendron leaf morphology. These findings will improve our understanding of the molecular regulatory mechanism of Liriodendron leaf morphogenesis and will also accelerate molecular breeding of Liriodendron.

The Net of Life, a short story: Intricate patterns of gene flows across hundreds of extant genomes, all the way to LUCA

Over the past quarter-century, the field of evolutionary biology has been transformed by the emergence of complete genome sequences and the conceptual framework known as the 'Net of Life.' This paradigm shift challenges traditional notions of evolution as a tree-like process, emphasizing the complex, interconnected network of gene flow that may blur the boundaries between distinct lineages. In this context, gene loss, rather than horizontal gene transfer, is the primary driver of gene content, with vertical inheritance playing a principal role. The 'Net of Life' not only impacts our understanding of genome evolution but also has profound implications for classification systems, the rapid appearance of new traits, and the spread of diseases. Here, we explore the core tenets of the 'Net of Life' and its implications for genome-scale phylogenetic divergence, providing a comprehensive framework for further investigations in evolutionary biology.

The contours of evolution: In defence of Darwin's tree of life paradigm

Both the concept of a Darwinian tree of life (TOL) and the possibility of its accurate reconstruction have been much criticized. Criticisms mostly revolve around the extensive occurrence of lateral gene transfer (LGT), instances of uptake of complete organisms to become organelles (with the associated subsequent gene transfer to the nucleus), as well as the implications of more subtle aspects of the biological species concept. Here we argue that none of these criticisms are sufficient to abandon the valuable TOL concept and the biological realities it captures. Especially important is the need to conceptually distinguish between organismal trees and gene trees, which necessitates incorporating insights of widely occurring LGT into modern evolutionary theory. We demonstrate that all criticisms, while based on important new findings, do not invalidate the TOL. After considering the implications of these new insights, we find that the contours of evolution are best represented by a TOL.

Genotypic Profiling, Functional Analysis, Cholesterol-Lowering Ability, and Angiotensin I-Converting Enzyme (ACE) Inhibitory Activity of Probiotic Lactiplantibacillus plantarum K25 via Different Approaches

Due to its alleged health advantages, several uses in biotechnology and food safety, the well-known probiotic strain Lactiplantibacillus plantarum K25 has drawn interest. This in-depth investigation explores the genetic diversity, makeup, and security characteristics of the microbial genome of L. plantarum K25, providing insightful knowledge about its genotypic profile and functional characteristics. Utilizing cutting-edge bioinformatics techniques like comparative genomics, pan-genomics, and genotypic profiling was carried out to reveal the strain's multidimensional potential in various fields. The results not only add to our understanding of the genetic makeup of L. plantarum K25 but also show off its acceptability in various fields, notably in biotechnology and food safety. The explanation of evolutionary links, which highlights L. plantarum K25's aptitude as a probiotic, is one notable finding from this research. Its safety profile, which is emphasized by the absence of genes linked to antibiotic resistance, is crucial and supports its status as a promising probiotic option.

A k-mer-based bulked segregant analysis approach to map seed traits in unphased heterozygous potato genomes

High-throughput sequencing-based methods for bulked segregant analysis (BSA) allow for the rapid identification of genetic markers associated with traits of interest. BSA studies have successfully identified qualitative (binary) and quantitative trait loci (QTLs) using QTL mapping. However, most require population structures that fit the models available and a reference genome. Instead, high-throughput short-read sequencing can be combined with BSA of k-mers (BSA-k-mer) to map traits that appear refractory to standard approaches. This method can be applied to any organism and is particularly useful for species with genomes diverged from the closest sequenced genome. It is also instrumental when dealing with highly heterozygous and potentially polyploid genomes without phased haplotype assemblies and for which a single haplotype can control a trait. Finally, it is flexible in terms of population structure. Here, we apply the BSA-k-mer method for the rapid identification of candidate regions related to seed spot and seed size in diploid potato. Using a mixture of F1 and F2 individuals from a cross between 2 highly heterozygous parents, candidate sequences were identified for each trait using the BSA-k-mer approach. Using parental reads, we were able to determine the parental origin of the loci. Finally, we mapped the identified k-mers to a closely related potato genome to validate the method and determine the genomic loci underlying these sequences. The location identified for the seed spot matches with previously identified loci associated with pigmentation in potato. The loci associated with seed size are novel. Both loci are relevant in future breeding toward true seeds in potato.

Functional analysis and comparative genomics of Rahnella perminowiae S11P1 and Variovorax sp. S12S4, two plant growth-promoting rhizobacteria isolated from Crocus sativus L. (saffron) rhizosphere

BACKGROUND

Rahnella perminowiae S11P1 and Variovorax sp. S12S4 are two plant growth-promoting rhizobacteria that were previously isolated from the rhizosphere of Crocus sativus L. (saffron), and have demonstrated interesting PGP activities and promising results when used as inoculants in field trials. To further elucidate the molecular mechanisms underlying their beneficial effects on plant growth, comprehensive genome mining of S11P1 and S12S4 and comparative genomic analysis with closely related strains were conducted.

RESULTS

Functional annotation of the two strains predicted a large number of genes involved in auxin and siderophore production, nitrogen fixation, sulfur metabolism, organic acid biosynthesis, pyrroloquinoline quinone production, 1-aminocyclopropane-1-carboxylate (ACC) deaminase activity, volatile organic compounds production, and polyamine biosynthesis. In addition, numerous genes implicated in plant-bacteria interactions, such as those involved in chemotaxis and quorum sensing, were predicted. Moreover, the two strains carried genes involved in bacterial fitness under abiotic stress conditions. Comparative genomic analysis revealed an open pan-genomic structure for the two strains. COG annotation showed that higher fractions of core and accessory genes were involved in the metabolism and transport of carbohydrates and amino acids, suggesting the metabolic versatility of the two strains as effective rhizosphere colonizers. Furthermore, this study reports the first comparison of Multilocus sequence analysis (MLSA) and core-based phylogenies of the Rahnella and Variovorax genera.

CONCLUSIONS

The present study unveils the molecular mechanisms underlying plant growth promotion and biocontrol activity of S11P1 and S12S4, and provides a basis for their further biotechnological application in agriculture.

Pan-Genome Plasticity and Virulence Factors: A Natural Treasure Trove for Acinetobacter baumannii

Acinetobacter baumannii is a Gram-negative pathogen responsible for a variety of community- and hospital-acquired infections. It is recognized as a life-threatening pathogen among hospitalized individuals and, in particular, immunocompromised patients in many countries. A. baumannii, as a member of the ESKAPE group, encompasses high genomic plasticity and simultaneously is predisposed to receive and exchange the mobile genetic elements (MGEs) through horizontal genetic transfer (HGT). Indeed, A. baumannii is a treasure trove that contains a high number of virulence factors. In accordance with these unique pathogenic characteristics of A. baumannii, the authors aim to discuss the natural treasure trove of pan-genome and virulence factors pertaining to this bacterial monster and try to highlight the reasons why this bacterium is a great concern in the global public health system.

Speeding genomic island discovery through systematic design of reference database composition

BACKGROUND

Genomic islands (GIs) are mobile genetic elements that integrate site-specifically into bacterial chromosomes, bearing genes that affect phenotypes such as pathogenicity and metabolism. GIs typically occur sporadically among related bacterial strains, enabling comparative genomic approaches to GI identification. For a candidate GI in a query genome, the number of reference genomes with a precise deletion of the GI serves as a support value for the GI. Our comparative software for GI identification was slowed by our original use of large reference genome databases (DBs). Here we explore smaller species-focused DBs.

RESULTS

With increasing DB size, recovery of our reliable prophage GI calls reached a plateau, while recovery of less reliable GI calls (FPs) increased rapidly as DB sizes exceeded ~500 genomes; i.e., overlarge DBs can increase FP rates. Paradoxically, relative to prophages, FPs were both more frequently supported only by genomes outside the species and more frequently supported only by genomes inside the species; this may be due to their generally lower support values. Setting a DB size limit for our SMAll Ranked Tailored (SMART) DB design speeded runtime ~65-fold. Strictly intra-species DBs would tend to lower yields of prophages for small species (with few genomes available); simulations with large species showed that this could be partially overcome by reaching outside the species to closely related taxa, without an FP burden. Employing such taxonomic outreach in DB design generated redundancy in the DB set; as few as 2984 DBs were needed to cover all 47894 prokaryotic species.

CONCLUSIONS

Runtime decreased dramatically with SMART DB design, with only minor losses of prophages. We also describe potential utility in other comparative genomics projects.

Diversification of gene content in the Mycobacterium tuberculosis complex is determined by phylogenetic and ecological signatures

We analyzed the pan-genome and gene content modulation of the most diverse genome data set of the Mycobacterium tuberculosis complex (MTBC) gathered to date. The closed pan-genome of the MTBC was characterized by reduced accessory and strain-specific genomes, compatible with its clonal nature. However, significantly fewer gene families were shared between MTBC genomes as their phylogenetic distance increased. This effect was only observed in inter-species comparisons, not within-species, which suggests that species-specific ecological characteristics are associated with changes in gene content. Gene loss, resulting from genomic deletions and pseudogenization, was found to drive the variation in gene content. This gene erosion differed among MTBC species and lineages, even within M. tuberculosis, where L2 showed more gene loss than L4. We also show that phylogenetic proximity is not always a good proxy for gene content relatedness in the MTBC, as the gene repertoire of Mycobacterium africanum L6 deviated from its expected phylogenetic niche conservatism. Gene disruptions of virulence factors, represented by pseudogene annotations, are mostly not conserved, being poor predictors of MTBC ecotypes. Each MTBC ecotype carries its own accessory genome, likely influenced by distinct selective pressures such as host and geography. It is important to investigate how gene loss confer new adaptive traits to MTBC strains; the detected heterogeneous gene loss poses a significant challenge in elucidating genetic factors responsible for the diverse phenotypes observed in the MTBC. By detailing specific gene losses, our study serves as a resource for researchers studying the MTBC phenotypes and their immune evasion strategies.IMPORTANCEIn this study, we analyzed the gene content of different ecotypes of the Mycobacterium tuberculosis complex (MTBC), the pathogens of tuberculosis. We found that changes in their gene content are associated with their ecological features, such as host preference. Gene loss was identified as the primary driver of these changes, which can vary even among different strains of the same ecotype. Our study also revealed that the gene content relatedness of these bacteria does not always mirror their evolutionary relationships. In addition, some genes of virulence can be variably lost among strains of the same MTBC ecotype, likely helping them to evade the immune system. Overall, our study highlights the importance of understanding how gene loss can lead to new adaptations in these bacteria and how different selective pressures may influence their genetic makeup.

Whole-genome sequencing of extensively drug-resistant Salmonella enterica serovar Typhi clinical isolates from the Peshawar region of Pakistan

BACKGROUND

Typhoid fever, caused by Salmonella enterica serovar Typhi, is a significant public health concern due to the escalating of antimicrobial resistance (AMR), with limited treatment options for extensively drug-resistant (XDR) S. Typhi strains pose a serious threat to disease management and control. This study aimed to investigate the genomic characteristics, epidemiology and AMR genes of XDR S. Typhi strains from typhoid fever patients in Pakistan.

METHODOLOGY

We assessed 200 patients with enteric fever symptoms, confirming 65 S. Typhi cases through culturing and biochemical tests. Subsequent antimicrobial susceptibility testing revealed 40 cases of extensively drug-resistant (XDR) and 25 cases of multi-drug resistance (MDR). Thirteen XDR strains were selected for whole-genome sequencing, to analyze their sequence type, phylogenetics, resistance genes, pathogenicity islands, and plasmid sequences using variety of data analysis resources. Pangenome analysis was conducted for 140 XDR strains, including thirteen in-house and 127 strains reported from other regions of Pakistan, to assess their genetic diversity and functional annotation.

RESULTS

MLST analysis classified all isolates as sequence type 1 (ST-1) with 4.3.1.1. P1 genotype characterization. Prophage and Salmonella Pathogenicity Island (SPI) analysis identified intact prophages and eight SPIs involved in Salmonella's invasion and replication within host cells. Genome data analysis revealed numerous AMR genes including dfrA7, sul1, qnrS1, TEM-1, Cat1, and CTX-M-15, and SNPs associated with antibiotics resistance. IncY, IncQ1, pMAC, and pAbTS2 plasmids, conferring antimicrobial resistance, were detected in a few XDR S. Typhi strains. Phylogenetic analysis inferred a close epidemiological linkage among XDR strains from different regions of Pakistan. Pangenome was noted closed among these strains and functional annotation highlighted genes related to metabolism and pathogenesis.

CONCLUSION

This study revealed a uniform genotypic background among XDR S. Typhi strains in Pakistan, signifying a persistence transmission of a single, highly antibiotic-resistant clone. The closed pan-genome observed underscores limited genetic diversity and highlights the importance of genomic surveillance for combating drug-resistant typhoid infections.