Pan-genome Analysis Research Articles

Bacterial Pangenome: A Review on the Current Strategies, Tools and Applications

The genomic research has grown by the advancement of sequencing technologies, which aim to capture the "totality" of life diversity by progressing from single strain to species or even higher taxa. Out of this scientific progress and inquisitiveness, the term "pangenome" was created. The goal of a pangenome is to represent all of a species genetic variation. It divides the genome into core, accessory and unique genes. Since the introduction of pangenomics, many software programs have been actively used for pangenome analysis. This review focuses on providing a snapshot of useful tools with their use, execution and computational details. In addition to tools for pangenome, the components, approaches, mechanism of variation, applications and challenges associated with pangenome analysis are also discussed in this article. Here, we have also discussed the increasing trend of bacterial pangenome studies with few studies and their significant findings. It has been proved that pangenomics is useful for understanding the clade of bacteria that helps in developing strategies for various application based on biological similarities and differences.

Identification of genetic markers of resistance to macrolide class antibiotics in Mannheimia haemolytica isolates from a Saskatchewan feedlot.

Mannheimia haemolytica is a major contributor to bovine respiratory disease (BRD), which causes substantial economic losses to the beef industry, and there is an urgent need for rapid and accurate diagnostic tests to provide evidence for treatment decisions and support antimicrobial stewardship. Diagnostic sequencing can provide information about antimicrobial resistance genes in M. haemolytica more rapidly than conventional diagnostics. Realizing the full potential of diagnostic sequencing requires a comprehensive understanding of the genetic markers of antimicrobial resistance. We identified genetic markers of resistance in M. haemolytica to macrolide class antibiotics commonly used for control of BRD. Genome sequences were determined for 99 M. haemolytica isolates with six different susceptibility phenotypes collected over 2 years from a feedlot in Saskatchewan, Canada. Known macrolide resistance genes estT, msr(E), and mph(E) were identified in most resistant isolates within predicted integrative and conjugative elements (ICEs). ICE sequences lacking antibiotic resistance genes were detected in 10 of 47 susceptible isolates. No resistance-associated polymorphisms were detected in ribosomal RNA genes, although previously unreported mutations in the L22 and L23 ribosomal proteins were identified in 12 and 27 resistant isolates, respectively. Pangenome analysis led to the identification of 79 genes associated with resistance to gamithromycin, of which 95% (75 of 79) had no functional annotation. Most of the observed phenotypic resistance was explained by previously identified antibiotic resistance genes, although resistance to the macrolides gamithromycin and tulathromycin was not explained in 39 of 47 isolates, demonstrating the need for continued surveillance for novel determinants of macrolide resistance.IMPORTANCEBovine respiratory disease is the costliest disease of beef cattle in North America and the most common reason for injectable antibiotic use in beef cattle. Metagenomic sequencing offers the potential to make economically significant reductions in turnaround time for diagnostic information for evidence-based selection of antibiotics for use in the feedlot. The success of diagnostic sequencing depends on a comprehensive catalog of antimicrobial resistance genes and other genome features associated with reduced susceptibility. We analyzed the genome sequences of isolates of Mannheimia haemolytica, a major bovine respiratory disease pathogen, and identified both previously known and novel genes associated with reduced susceptibility to macrolide class antimicrobials. These findings reinforce the need for ongoing surveillance for markers of antimicrobial resistance to support improved diagnostics and antimicrobial stewardship.

Pangenomic analysis identifies correlations between Akkermansia species and subspecies and human health outcomes.

Aim: Akkermansia are common members of the human gastrointestinal microbiota. The prevalence of these mucophilic bacteria, especially Akkermansia muciniphila (A. muciniphila), correlates with immunological and metabolic health. The genus Akkermansia in humans includes species with significantly larger genomes than A. muciniphila, leading us to postulate that this added genetic content may influence how they impact human metabolic and immunological health. Methods: We conducted a pangenomic analysis of 234 Akkermansia complete or near-complete genomes. We also used high-resolution species and subspecies assignments to reanalyze publicly available metagenomic datasets to determine if there are relationships between Akkermansia species and A. muciniphila clades with various disease outcomes. Results: Analysis of genome-wide average nucleotide identity, 16S rRNA gene identity, conservation of core Akkermansia genes, and analysis of the fatty acid composition of representative isolates support the partitioning of the genus Akkermansia into several species. In addition, A. muciniphila sensu stricto, the most prevalent Akkermansia species in humans, should be subdivided into two subspecies. For a pediatric cohort, we observed species-specific correlations between Akkermansia abundance with baseline obesity or after various interventions. For inflammatory bowel disease cohorts, we identified a decreased abundance of Akkermansia in patients with ulcerative colitis or Crohn's disease, which was species and subspecies-dependent. In patients undergoing immune checkpoint inhibitor therapies for non-small cell lung carcinoma, we observed a significant association between one A. muciniphila subspecies and survival outcomes. Conclusion: Our findings suggest that the prevalence of specific Akkermansia species and/or subspecies can be crucial in evaluating their association with human health, particularly in different disease contexts, and is an important consideration for their use as probiotics.

Whole-Genome Sequence Analysis of Antibiotic Resistance, Virulence, and Plasmid Dynamics in Multidrug-Resistant E. coli Isolates from Imported Shrimp.

We analyzed antimicrobial resistance and virulence traits in multidrug-resistant (MDR) E. coli isolates obtained from imported shrimp using whole-genome sequences (WGSs). Antibiotic resistance profiles were determined phenotypically. WGSs identified key characteristics, including their multilocus sequence type (MLST), serotype, virulence factors, antibiotic resistance genes, and mobile elements. Most of the isolates exhibited resistance to gentamicin, streptomycin, ampicillin, chloramphenicol, nalidixic acid, ciprofloxacin, tetracycline, and trimethoprim/sulfamethoxazole. Multilocus sequence type (MLST), serotype, average nucleotide identity (ANI), and pangenome analysis showed high genomic similarity among isolates, except for EC15 and ECV01. The EC119 plasmid contained a variety of efflux pump genes, including those encoding the acid resistance transcriptional activators (gadE, gadW, and gadX), resistance-nodulation-division-type efflux pumps (mdtE and mdtF), and a metabolite, H1 symporter (MHS) family major facilitator superfamily transporter (MNZ41_23075). Virulence genes displayed diversity, particularly EC15, whose plasmids carried genes for adherence (faeA and faeC-I), invasion (ipaH and virB), and capsule (caf1A and caf1M). This comprehensive analysis illuminates antimicrobial resistance, virulence, and plasmid dynamics in E. coli from imported shrimp and has profound implications for public health, emphasizing the need for continued surveillance and research into the evolution of these important bacterial pathogens.

Morphological, Metabolomic and Genomic Evidences on Drought Stress Protective Functioning of the Endophyte Bacillus safensis Ni7.

The metabolomic and genomic characterization of an endophytic Bacillus safensis Ni7 was carried out in this study. This strain has previously been isolated from the xerophytic plant Nerium indicum L. and reported to enhance the drought tolerance in Capsicum annuum L. seedlings. The effects of drought stress on the morphology, biofilm production, and metabolite production of B. safensis Ni7 are analyzed in the current study. From the results obtained, the organism was found to have multiple strategies such as aggregation and clumping, robust biofilm production, and increased production of surfactin homologues under the drought induced condition when compared to non-stressed condition. Further the whole genome sequencing (WGS) based analysis has demonstrated B. safensis Ni7 to have a genome size of 3,671,999bp, N50 value of 3,527,239, and a mean G+C content of 41.58%. Interestingly the organism was observed to have the presence of various stress-responsive genes (13, 20U, 16U,160, 39, 17M, 18, 26, and ctc) and genes responsible for surfactin production (srfAA, srfAB, srfAC, and srfAD), biofilm production (epsD, epsE, epsF, epsG, epsH, epsI, epsK, epsL, epsM, epsN, and pel), chemotaxis (cheB_1, cheB_2, cheB_3, cheW_1, cheW_2 cheR, cheD, cheC, cheA, cheY, cheV, and cheB_4), flagella synthesis (flgG_1, flgG_2, flgG_3, flgC, and flgB) as supportive to the drought tolerance. Besides these, the genes responsible for plant growth promotion (PGP), including the genes for nitrogen (nasA, nasB, nasC, nasD, and nasE) and sulfur assimilation (cysL_1&L_2, cysI) and genes for phosphate solubilization (phoA, phoP_1& phoP_2, and phoR) could also be predicted. Along with the same, the genes for catalase, superoxide dismutase, protein homeostasis, cellular fitness, osmoprotectants production, and protein folding could also be predicted from its WGS data. Further pan-genome analysis with plant associated B. safensis strains available in the public databases revealed B. safensis Ni7 to have the presence of a total of 5391 gene clusters. Among these, 3207 genes were identified as core genes, 954 as shell genes and 1230 as cloud genes. This variation in gene content could be taken as an indication of evolution of strains of Bacillus safensis as per specific conditions and hence in the case of B. safensis Ni7 its role in habitat adaptation of plant is well expected. This diversity in endophytic bacterial genes may attribute its role to support the plant system to cope up with stress conditions. Overall, the study provides genomic evidence on Bacillus safensis Ni7 as a stress alleviating microbial partner in plants.

Phylogenetic and pangenomic analyses of members of the family Micrococcaceae related to a plant-growth-promoting rhizobacterium isolated from the rhizosphere of potato (Solanum tuberosum L.).

We report the results of taxonomic studies on members of the family Micrococcaceae that, according to the 16S rRNA, internal transcribed spacer 1 (ITS1), average nucleotide identity (ANI), and average amino acid identity (AAI) tests, are related to Kocuria rosea strain RCAM04488, a plant-growth-promoting rhizobacterium (PGPR) isolated from the rhizosphere of potato (Solanum tuberosum L.). In these studies, we used whole-genome phylogenetic tests and pangenomic analysis. According to the ANI > 95 % criterion, several known members of K. salina, K. polaris, and K. rosea (including K. rosea type strain ATCC 186T) that are related most closely to isolate RCAM04488 in the ITS1 test should be assigned to the same species with appropriate strain verification. However, these strains were isolated from strongly contrasting ecological and geographical habitats, which could not but affect their genotypes and phenotypes and which should be taken into account in evaluation of their systematic position. This contradiction was resolved by a pangenomic analysis, which showed that the strains differed strongly in the number of accessory and strain-specific genes determining their individuality and possibly their potential for adaptation to different ecological niches. Similar results were obtained in a full-scale AAI test against the UniProt database (about 250 million records), by using the AAI-profiler program and the proteome of K. rosea strain ATCC 186T as a query. According to the AAI > 65 % criterion, members of the genus Arthrobacter and several other genera belonging to the class Actinomycetes, with a very wide geographical and ecological range of sources of isolation, should be placed into the same genus as Kocuria. Within the paradigm with vertically inherited phylogenetic markers, this could be regarded as a signal for their following taxonomic reclassification. An important factor in this case may be the detailing of the gene composition of the strains and the taxonomic ratios resulting from analysis of the pangenomes of the corresponding clades.

Cofitness network connectivity determines a fuzzy essential zone in open bacterial pangenome.

Most in silico evolutionary studies commonly assumed that core genes are essential for cellular function, while accessory genes are dispensable, particularly in nutrient-rich environments. However, this assumption is seldom tested genetically within the pangenome context. In this study, we conducted a robust pangenomic Tn-seq analysis of fitness genes in a nutrient-rich medium for Sinorhizobium strains with a canonical open pangenome. To evaluate the robustness of fitness category assignment, Tn-seq data for three independent mutant libraries per strain were analyzed by three methods, which indicates that the Hidden Markov Model (HMM)-based method is most robust to variations between mutant libraries and not sensitive to data size, outperforming the Bayesian and Monte Carlo simulation-based methods. Consequently, the HMM method was used to classify the fitness category. Fitness genes, categorized as essential (ES), advantage (GA), and disadvantage (GD) genes for growth, are enriched in core genes, while nonessential genes (NE) are over-represented in accessory genes. Accessory ES/GA genes showed a lower fitness effect than core ES/GA genes. Connectivity degrees in the cofitness network decrease in the order of ES, GD, and GA/NE. In addition to accessory genes, 1599 out of 3284 core genes display differential essentiality across test strains. Within the pangenome core, both shared quasi-essential (ES and GA) and strain-dependent fitness genes are enriched in similar functional categories. Our analysis demonstrates a considerable fuzzy essential zone determined by cofitness connectivity degrees in Sinorhizobium pangenome and highlights the power of the cofitness network in understanding the genetic basis of ever-increasing prokaryotic pangenome data.

Human epicardial adipose tissue contains innate and adaptive lymphoid cells and a higher proportion of innate type 2 lymphoid cells compared to other adipose tissues

ImportanceEpicardial adipose tissue (EAT) is a biologically active organ surrounding myocardium and coronary arteries that has been associated with coronary artery disease (CAD) and atrial fibrillation. Previous work has shown that EAT exhibits beige features. ObjectiveOur objective was to determine whether the stromal vascular fraction of the human EAT contains innate or adaptive lymphoid cells compared to thoracic subcutaneous (thSAT), visceral abdominal (VAT) and subcutaneous abdominal (abSAT). ParticipantsNew pangenomic microarray analysis was performed on previous transcriptomic dataset using significance analysis of microarray and ingenuity pathway analysis (n=41) to identify specific immune signature and its link with browning genes. EAT, thSAT, VAT and abSAT samples from explanted patients with severe cardiomyopathies and multi-organ donor patients (n=17) were used for flow cytometry (FC) immunophenotyping assay. Patients were on average 55±16 years-old; 47% had hypertension and 6% CAD. Phenotypic adaptive and innate immune profiles were performed using a TBNK panel and a specific ILC1-2-3 panel including CD127, CD117, CRTH2 (CD294) and activation markers such as CD25 and CD69. ResultsTranscriptomic analysis showed a significant positive correlation between the TH2 immune pathway (IL-4, IL-5, IL-13, IL-25, IL-33) and browning genes (UCP-1, PRDM16, TMEM26, CITED1, TBX1) in EAT versus thSAT (R=0.82, P<0.0001). Regarding adaptive immune cells, a preponderance of CD8T cells, a contingent of CD4T cells, and a few B cells were observed in all ATs (P<0.0001). In innate lymphoid cells (ILCs), an increase was observed in visceral ATs (i.e. EAT; VAT 35±8ILCs/g of tissue) compared to their subcutaneous counterpart (i.e. thSAT+abSAT: 8±3 ILCs/g of AT, P=0.002), with a difference in the proportion of the 3 subtypes of ILCs (ILC1>ILC3>ILC2). In addition, we observed an increase in EAT-ILC2 compared to other ATs and almost all these EAT-ILC2 expressed CD69 and/or CD25 activation markers (99.75±0.16%; P<0.0001). We also observed more NKs in EAT and VAT (1520±71 cells/g of AT) than in SATs (562±17 cells/g of AT); P=0.01. ConclusionThis is the first study to provide a comparison between innate and adaptive lymphoid cells in human epicardial versus abdominal or thoracic adipose tissues. Further studies are ongoing to decipher whether these cells could be involved in EAT beiging. Trial RegistrationCODECOH No. DC-2021-4518 The French agency of biomedicine PFS21-005.

The mobilome of Staphylococcus aureus from wild ungulates reveals epidemiological links at the animal-human interface

Staphylococcus aureus thrives at animal-human-environment interfaces. A large-scale work from our group indicated that antimicrobial resistance (AMR) in commensal S. aureus strains from wild ungulates is associated with agricultural land cover and livestock farming, raising the hypothesis that AMR genes in wildlife strains may originate from different hosts, namely via exchange of mobile genetic elements (MGE). In this work, we generate the largest available dataset of S. aureus draft genomes from wild ungulates in Portugal and explore their mobilome, which can determine important traits such as AMR, virulence, and host specificity, to understand MGE exchange. Core genome multi-locus sequence typing based on 98 newly generated draft genomes and 101 publicly available genomes from Portugal demonstrated that the genomic relatedness of S. aureus from wild ungulates assigned to livestock-associated sequence types (ST) is greater compared to wild ungulate isolates assigned to human-associated STs. Screening of host specificity determinants disclosed the unexpected presence in wildlife of the immune evasion cluster encoded in φSa3 prophage, described as a human-specific virulence determinant. Additionally, two plasmids, pAVX and pETB, previously associated with avian species and humans, respectively, and the Tn553 transposon were detected. Both pETB and Tn553 encode penicillin resistance through blaZ. Pangenome analysis of wild ungulate isolates shows a core genome fraction of 2133 genes, with isolates assigned to ST72 and ST3224 being distinguished from the remaining by MGEs, although there is no reported role of these in adaptation to wildlife. AMR related gene clusters found in the shell genome are directly linked to resistance against penicillin, macrolides, fosfomycin, and aminoglycosides, and they represent mobile ARGs. Altogether, our findings support epidemiological interactions of human and non-human hosts at interfaces, with MGE exchange, including AMR determinants, associated with putative indirect movements of S. aureus among human and wildlife hosts that might be bridged by livestock.

Comparative Genomics of the First Resistant Candida auris Strain Isolated in Mexico: Phylogenomic and Pan-Genomic Analysis and Mutations Associated with Antifungal Resistance.

Candida auris is an emerging multidrug-resistant and opportunistic pathogenic yeast. Whole-genome sequencing analysis has defined five major clades, each from a distinct geographic region. The current study aimed to examine the genome of the C. auris 20-1498 strain, which is the first isolate of this fungus identified in Mexico. Based on whole-genome sequencing, the draft genome was found to contain 70 contigs. It had a total genome size of 12.86 Mbp, an N50 value of 1.6 Mbp, and an average guanine-cytosine (GC) content of 45.5%. Genome annotation revealed a total of 5432 genes encoding 5515 proteins. According to the genomic analysis, the C. auris 20-1498 strain belongs to clade IV (containing strains endemic to South America). Of the two genes (ERG11 and FKS1) associated with drug resistance in C. auris, a mutation was detected in K143R, a gene located in a mutation hotspot of ERG11 (lanosterol 14-α-demethylase), an antifungal drug target. The focus on whole-genome sequencing and the identification of mutations linked to the drug resistance of fungi could lead to the discovery of new therapeutic targets and new antifungal compounds.

Genome characterisation and comparative analysis of Schaalia dentiphila sp. nov. and its subspecies, S. dentiphila subsp. denticola subsp. nov., from the human oral cavity.

Schaalia species are primarily found among the oral microbiota of humans and other animals. They have been associated with various infections through their involvement in biofilm formation, modulation of host responses, and interaction with other microorganisms. In this study, two strains previously indicated as Actinomyces spp. were found to be novel members of the genus Schaalia based on their whole genome sequences. Whole-genome sequencing revealed both strains with a genome size of 2.3 Mbp and GC contents of 65.5%. Phylogenetics analysis for taxonomic placement revealed strains NCTC 9931 and C24 as distinct species within the genus Schaalia. Overall genome-relatedness indices including digital DNA-DNA hybridization (dDDH), and average nucleotide/amino acid identity (ANI/AAI) confirmed both strains as distinct species, with values below the species boundary thresholds (dDDH < 70%, and ANI and AAI < 95%) when compared to nearest type strain Schaalia odontolytica NCTC 9935T. Pangenome and orthologous analyses highlighted their differences in gene properties and biological functions compared to existing type strains. Additionally, the identification of genomic islands (GIs) and virulence-associated factors indicated their genetic diversity and potential adaptive capabilities, as well as potential implications for human health. Notably, CRISPR-Cas systems in strain NCTC 9931 underscore its adaptive immune mechanisms compared to strain C24. Based on these findings, strain NCTC 9931T (= ATCC 17982T = DSM 43331T = CIP 104728T = CCUG 18309T = NCTC 14978T = CGMCC 1.90328T) represents a novel species, for which the name Schaalia dentiphila subsp. dentiphila sp. nov. subsp. nov. is proposed, while strain C24T (= NCTC 14980T = CGMCC 1.90329T) represents a distinct novel subspecies, for which the name Schaalia dentiphila subsp. denticola. subsp. nov. is proposed. This study enriches our understanding of the genomic diversity of Schaalia species and paves the way for further investigations into their roles in oral health. This research reveals two Schaalia strains, NCTC 9931T and C24T, as novel entities with distinct genomic features. Expanding the taxonomic framework of the genus Schaalia, this study offers a critical resource for probing the metabolic intricacies and resistance patterns of these bacteria. This work stands as a cornerstone for microbial taxonomy, paving the way for significant advances in clinical diagnostics.

Putative Probiotic Ligilactobacillus salivarius Strains Isolated from the Intestines of Meat-Type Pigeon Squabs.

This study aims to screen for potential probiotic lactic acid bacteria from the intestines of meat-type pigeon squabs. Ligilactobacillus salivarius YZU37 was identified as the best comprehensive performed strain. Being acid- and bile salt-tolerant, it displayed growth-inhibition activities against Staphylococcus aureus ATCC25923, Escherichia coli ATCC25922, and Salmonella typhimurium SL1344, exhibited sensitivity to 6 commonly used antibiotics, and endowed with good cell surface hydrophobicity, auto-aggregation property, and anti-oxidant activities. Results of in vitro experiments indicated that the bacteriostatic effects of this strain were related to the production of proteinaceous substances that depend on acidic conditions. Whole-genome sequencing of L. salivarius YZU37 was performed to elucidate the genetic basis underlying its probiotic potential. Pangenome analysis of L. salivarius YZU37 and other 212 L. salivarius strains available on NCBI database revealed a pigeon-unique gene coding choloylglycine hydrolase (CGH), which had higher enzyme-substrate binding affinity than that of the common CGH shared by L. salivarius strains of other sources. Annotation of the functional genes in the genome of L. salivarius YZU37 revealed genes involved in responses to acid, bile salt, heat, cold, heavy metal, and oxidative stresses. The whole genome analysis also revealed the absence of virulence and toxin genes and the presence of 65 genes distributed under 4 CAZymes classes, 2 CRISPR-cas regions, and 3 enterolysin A clusters which may confer the acid-dependent antimicrobial potential of L. salivarius YZU37. Altogether, our results highlighted the probiotic potential of L. salivarius YZU37. Further in vivo investigations are required to elucidate its beneficial effects on pigeons.

Microbiological and molecular studies on a multidrug-resistant Pseudomonas aeruginosa from a liver transplant patient with urinary tract infection in Egypt.

Pseudomonas aeruginosa is an opportunistic pathogen responsible for complicated UTIs and exhibits high antibiotic resistance, leading to increased mortality rates, especially in cases of multidrug-resistant strains. This study aimed to investigate the antibiotic susceptibility patterns and genomic characterization of XDR strains identified in end-stage liver disease patients who underwent liver transplants. In this study, a number of 30 individuals who underwent liver transplants were registered. Ninety urine and 60 wound site swab samples were collected and processed for culturing, identification, and antimicrobial sensitivity. Extensively drug-resistant strain EMARA01 was confirmed through Sanger sequencing and was then processed for whole genome sequencing to characterize the genomic pattern. Sequencing data were processed for de novo assembly using various tools and databases, including genome annotation, serotype identification, virulence factor genes, and antimicrobial resistance gene. Pangenome analysis of randomly selected 147 reference strains and EMAR01 sequenced strain was performed using the Bacterial Pan Genome Analysis (BPGA) software. Of these total examined samples, nosocomial infection due to P. aeruginosa was detected in twelve patients' samples. AST analysis showed that P. aeruginosa strains exhibit resistance to tobramycin, erythromycin, and gentamicin, followed by piperacillin and ofloxacin, and no strains exhibit resistance to meropenem and imipenem. The CARD database identified 59 AMR genes similar to the EMAR01 strain genome and mostly belong to the family involved in the resistance-nodulation-cell division (RND) antibiotic efflux pump. Five genes; nalC, nalD, MexR, MexA, and MexB, exhibit resistance to 14 classes of antibiotics, while two AMR; CpxR, and OprM, exhibit resistance to 15 classes of drugs. Pangenome analysis revealed that the pan-genome remained open, suggesting the potential for acquiring accessory and unique genes. Notably, the genes predominantly involved in amino acid transport metabolism were identified using the KEGG database. This study provides valuable insights into the antimicrobial resistance profile, genetic features, and genomic evolution of P. aeruginosa strains causing UTIs in liver transplant patients. The findings emphasize the significance of comprehending AMR mechanisms and genetic diversity in P. aeruginosa for developing effective treatment strategies and infection control measures.

The complete catalog of antimicrobial resistance secondary active transporters in Clostridioides difficile: evolution and drug resistance perspective

Secondary active transporters shuttle substrates across eukaryotic and prokaryotic membranes, utilizing different electrochemical gradients. They are recognized as one of the antimicrobial efflux pumps among pathogens. While primary active transporters within the genome of C. difficile 630 have been completely cataloged, the systematical study of secondary active transporters remains incomplete. Here, we not only identify secondary active transporters but also disclose their evolution and role in drug resistance in C. difficile 630. Our analysis reveals that C. difficile 630 carries 147 secondary active transporters belonging to 27 (super)families. Notably, 50 (34%) of them potentially contribute to antimicrobial resistance (AMR). AMR-secondary active transporters are structurally classified into five (super)families: the p-aminobenzoyl-glutamate transporter (AbgT), drug/metabolite transporter (DMT) superfamily, major facilitator (MFS) superfamily, multidrug and toxic compound extrusion (MATE) family, and resistance-nodulation-division (RND) family. Surprisingly, complete RND genes found in C. difficile 630 are likely an evolutionary leftover from the common ancestor with the diderm. Through protein structure comparisons, we have potentially identified six novel AMR-secondary active transporters from DMT, MATE, and MFS (super)families. Pangenome analysis revealed that half of the AMR-secondary transporters are accessory genes, which indicates an important role in adaptive AMR function rather than innate physiological homeostasis. Gene expression profile firmly supports their ability to respond to a wide spectrum of antibiotics. Our findings highlight the evolution of AMR-secondary active transporters and their integral role in antibiotic responses. This marks AMR-secondary active transporters as interesting therapeutic targets to synergize with other antibiotic activity.

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

Salmonella enterica serovar Infantis (Salmonella Infantis) is a zoonotic, ubiquitous and foodborne pathogen of worldwide distribution. Despite Brazil's relevance as a major meat exporter, few studies were conducted to characterize strains of this serovar by genomic analyses in this country. Therefore, this study aimed to assess the diversity of 80 Salmonella Infantis strains isolated from veterinary, food and human sources in Brazil between 2013 and 2018 by comparative genomic analyses. Additional genomes of non-Brazilian countries (n=18) were included for comparison purposes in some analyses. Analyses of whole-genome multi-locus sequence typing (wgMLST), using PGAdb-builder, and of fragmented genomes, using Gegenees, were conducted to compare the 80 Brazilian strains to the 18 non-Brazilian genomes. Pangenome analyses and calculations were performed for all Salmonella Infantis genomes analyzed. The presence of prophages was determined using PHASTER for the 80 Brazilian strains. The genome plasticity using BLAST Ring Image Generator (BRIG) and gene synteny using Mauve were evaluated for 20 selected Salmonella Infantis genomes from Brazil and ten from non-Brazilian countries. Unique orthologous protein clusters were searched in ten selected Salmonella Infantis genomes from Brazil and ten from non-Brazilian countries. wgMLST and Gegenees showed a high genomic similarity among some Brazilian Salmonella Infantis genomes, and also the correlation of some clusters with non-Brazilian genomes. Gegenees also showed an overall similarity >91% among all Salmonella Infantis genomes. Pangenome calculations revealed an open pangenome for all Salmonella Infantis subsets analyzed and a high gene content in the core genomes. Fifteen types of prophages were detected among 97.5% of the Brazilian strains. BRIG and Mauve demonstrated a high structural similarity among the Brazilian and non-Brazilian isolates. Unique orthologous protein clusters related to biological processes, molecular functions, and cellular components were detected among Brazilian and non-Brazilian genomes. The results presented using different genomic approaches emphasized the significant genomic similarity among Brazilian Salmonella Infantis genomes analyzed, suggesting wide distribution of closely related genotypes among diverse sources in Brazil. The data generated contributed to novel information regarding the genomic diversity of Brazilian and non-Brazilian Salmonella Infantis in comparison. The different genetically related subtypes of Salmonella Infantis from Brazil can either occur exclusively within the country, or also in other countries, suggesting that some exportation of the Brazilian genotypes may have already occurred.

A de novo germline pathogenic BRCA1 variant identified following an osteosarcoma pangenomic molecular analysis.

De novo germline pathogenic variants (gPV) of the BReast CAncer 1 (BRCA1) gene are very rare. Only a few have been described up to date, usually in patients with a history of ovarian or breast cancer. Here, we report the first case of an incidental de novo BRCA1 germline pathogenic variant which was identified within the framework of the Plan France Médecine Génomique (PFMG) 2025 French national tumor sequencing program. The proband was a 29-year-old man diagnosed with metastatic osteosarcoma. Tumor whole exome sequencing identified a BRCA1 c.3756_3759del p.(Ser1253Argfs*10) pathogenic variant without loss-of-heterozygosity. A low genomic instability score and the absence of single base substitution signatures of homologous recombination deficiency suggested that the BRCA1 variant was not driver in the osteosarcoma tumorigenesis. Germline whole genome sequencing asserted the germline nature of this variant, with a 36% allele frequency, suggesting a mosaicism caused by a post-zygotic mutational event. The proband's family (parents and siblings) were not carriers of this variant confirming the de novo occurrence. Tumor sequencing programs like the French PFMG 2025 have been implemented worldwide and may help identify new gPV, including de novo variants.

Divergence within the Taxon 'Candidatus Phytoplasma asteris' Confirmed by Comparative Genome Analysis of Carrot Strains.

Phytoplasmas are linked to diseases in hundreds of economically important crops, including carrots. In carrots, phytoplasmosis is associated with leaf chlorosis and necrosis, coupled with inhibited root system development, ultimately leading to significant economic losses. During a field study conducted in Baden-Württemberg (Germany), two strains of the provisional taxon 'Candidatus Phytoplasma asteris' were identified within a carrot plot. For further analysis, strains M8 and M33 underwent shotgun sequencing, utilising single-molecule-real-time (SMRT) long-read sequencing and sequencing-by-synthesis (SBS) paired-end short-read sequencing techniques. Hybrid assemblies resulted in complete de novo assemblies of two genomes harboring circular chromosomes and two plasmids. Analyses, including average nucleotide identity and sequence comparisons of established marker genes, confirmed the phylogenetic divergence of 'Ca. P. asteris' and a different assignment of strains to the 16S rRNA subgroup I-A for M33 and I-B for M8. These groups exhibited unique features, encompassing virulence factors and genes, associated with the mobilome. In contrast, pan-genome analysis revealed a highly conserved gene set related to metabolism across these strains. This analysis of the Aster Yellows (AY) group reaffirms the perception of phytoplasmas as bacteria that have undergone extensive genome reduction during their co-evolution with the host and an increase of genome size by mobilome.

First Pangenome of Corynebacterium rouxii, a Potentially Toxigenic Species of Corynebacterium diphtheriae Complex

Corynebacterium rouxii is one of the recently described species of the Corynebacterium diphtheriae complex. As this species can potentially infect different hosts and harbor the tox gene, producing diphtheria toxin, we present its first pangenomic analysis in this work. A total of fifteen genomes deposited in online databases were included. After confirming the taxonomic position of the isolates by genomic taxonomy, the genomes were submitted to genomic plasticity, gene synteny, and pangenome prediction analyses. In addition, virulence and antimicrobial resistance genes were investigated. Finally, epidemiological data were obtained through molecular typing, clustering, and phylogenetic analysis. Our data demonstrated genetic diversity within the species with low synteny. However, the gene content is extensively conserved, and the pangenome is composed of 2606 gene families, of which 1916 are in the core genome and 80 are related to unique genes. Prophages, insertion sequences, and genomic islands were found. A type I-E CRISPR-Cas system was also detected. Besides the tox gene, determinants involved in adhesion and iron acquisition and two putative antimicrobial resistance genes were predicted. These findings provide valuable insight about this species’ pathogenicity, evolution, and diversity. In the future, our data can contribute to different areas, including vaccinology and epidemiology.

Antimicrobial resistance and population genomics of emerging multidrug-resistant Salmonella 4,[5],12:i:- in Guangdong, China.

Salmonella 4,[5],12:i:-, a monophasic variant of Salmonella Typhimurium, has emerged as a global cause of multidrug-resistant salmonellosis and has become endemic in many developing and developed countries, especially in China. Here, we have sequenced 352 clinical isolates in Guangdong, China, during 2009-2019 and performed a large-scale collection of Salmonella 4,[5],12:i:- with whole genome sequencing (WGS) data across the globe, to better understand the population structure, antimicrobial resistance (AMR) genomic characterization, and transmission routes of Salmonella 4,[5],12:i:- across Guangdong. Salmonella 4,[5],12:i:- strains showed broad genetic diversity; Guangdong isolates were found to be widely distributed among the global lineages. Of note, we identified the formation of a novel Guangdong clade (Bayesian analysis of population structure lineage 1 [BAPS1]) genetically diversified from the global isolates and likely emerged around 1990s. BAPS1 exhibits unique genomic features, including large pan-genome, decreased ciprofloxacin susceptibility due to mutation in gyrA and carriage of plasmid-mediated quinolone resistance (PMQR) genes, and the multidrug-resistant IncHI2 plasmid. Furthermore, high genetic similarity was found between strains collected from Guangdong, Europe, and North America, indicating the association with multiple introductions from overseas. These results suggested that global dissemination and local clonal expansion simultaneously occurred in Guangdong, China, and horizontally acquired resistance to first-line and last-line antimicrobials at local level, underlying emergences of extensive drug and pan-drug resistance. Our findings have increased the knowledge of global and local epidemics of Salmonella 4,[5],12:i:- in Guangdong, China, and provided a comprehensive baseline data set essential for future molecular surveillance.IMPORTANCESalmonella 4,[5],12:i:- has been regarded as the predominant pandemic serotype causing diarrheal diseases globally, while multidrug resistance (MDR) constitutes great public health concerns. This study provided a detailed and comprehensive genome-scale analysis of this important Salmonella serovar in the past decade in Guangdong, China. Our results revealed the complexity of two distinct transmission modes, namely global transmission and local expansion, circulating in Guangdong over a decade. Using phylogeography models, the origin of Salmonella 4,[5],12:i:- was predicted from two aspects, year and country, that is, Salmonella 4,[5],12:i:- emerged in 1983, and was introduced from the UK, and subsequently differentiated into the local endemic lineage circa 1991. Additionally, based on the pan-genome analysis, it was found that the gene accumulation rate in local endemic BAPS 1 lineage was higher than in other lineages, and the horizontal transmission of MDR IncHI2 plasmid associated with high resistance played a major role, which showed the potential threat to public health.

Human epicardial fat has a beige profile and contains higher type 2 innate lymphoid cells than subcutaneous fat.

Epicardial adipose tissue (EAT) is a visceral fat that has been associated with coronary artery disease and atrial fibrillation. Previous work has revealed that EAT exhibits beige features. First, a new pan-genomic microarray analysis was performed on previously collected paired human EAT and thoracic subcutaneous AT (thSAT) from the EPICAR study (n = 31) to decipher a specific immune signature and its link with browning genes. Then, adaptive (T and B cells) and innate lymphoid cell (ILC1, ILC2, and ILC3) immunophenotyping assay panels, including CD127, CD117, and prostaglandin D2 receptor 2, were performed on prospectively collected paired human multiorgan donors (n = 18; INTERFACE study). In the EPICAR study, a positive correlation between the T helper cell subtype Th2 immune pathway and browning genes was found in EAT versus thSAT (r = 0.82; p < 0.0001). In the INTERFACE study, this correlation was also observed (r = 0.31; p = 0.017), and a preponderance of CD4+T cells, CD8+T cells, and a few B cells was observed in all ATs (p < 0.0001). An increase in ILCs was observed in visceral AT (VAT) (i.e., EAT + VAT; 30 ± 5 ILCs per gram of AT) compared with subcutaneous counterparts (i.e., thSAT + abdominal SAT; 8 ± 2 ILCs per gram of AT; p = 0.001), with ILC1 being the most frequent (ILC1 > ILC3 > ILC2). Numbers of ILCs per gram of AT correlated with several Th2 or browning genes (IL-13, TNF receptor superfamily member 9 [TNFRSF9], and alkaline phosphatase, biomineralization associated [ALPL]). Interestingly, a specific increase in EAT-ILC2 compared with other ATs was observed, including a significant proportion expressing CD69 and/or CD25 activation markers (97.9% ± 1.2%; p < 0.0001). Finally, more natural killer cells were observed in EAT + VAT than in thSAT + abdominal SAT (p = 0.01). Exclusion of patients with coronary artery disease in the EPICAR and INTERFACE studies did not modify the main findings. Gene expression phenotyping confirmed specific upregulation of Th2 pathway and browning genes (IL-33 and uncoupling protein 1 [UCP-1]) in EAT. This is the first study, to our knowledge, to provide a comparison between innate and adaptive lymphoid cells in human EAT. Further studies are ongoing to decipher whether these cells could be involved in EAT beiging.