Pan-genome Analysis Research Articles

Two new clades recovered at high temperatures provide novel phylogenetic and genomic insights into Candidatus Accumulibacter.

Candidatus Accumulibacter, a key genus of polyphosphate-accumulating organisms, plays key roles in lab- and full-scale enhanced biological phosphorus removal (EBPR) systems. A total of 10 high-quality Ca. Accumulibacter genomes were recovered from EBPR systems operated at high temperatures, providing significantly updated phylogenetic and genomic insights into the Ca. Accumulibacter lineage. Among these genomes, clade IIF members SCELSE-3, SCELSE-4, and SCELSE-6 represent the to-date known genomes encoding a complete denitrification pathway, suggesting that Ca. Accumulibacter alone could achieve complete denitrification. Clade IIC members SSA1, SCUT-1, SCELCE-2, and SCELSE-8 lack the entire set of denitrifying genes, representing to-date known non-denitrifying Ca. Accumulibacter. A pan-genomic analysis with other Ca. Accumulibacter members suggested that all Ca. Accumulibacter likely has the potential to use dicarboxylic amino acids. Ca. Accumulibacter aalborgensis AALB and Ca. Accumulibacter affinis BAT3C720 seemed to be the only two members capable of using glucose for EBPR. A heat shock protein Hsp20 encoding gene was found exclusively in genomes recovered at high temperatures, which was absent in clades IA, IC, IG, IIA, IIB, IID, IIG, and II-I members. High transcription of this gene in clade IIC members SCUT-2 and SCUT-3 suggested its role in surviving high temperatures for Ca. Accumulibacter. Ambiguous clade identity was observed for newly recovered genomes (SCELSE-9 and SCELSE-10). Five machine learning models were developed using orthogroups as input features. Prediction results suggested that they belong to a new clade (IIK). The phylogeny of Ca. Accumulibacter was re-evaluated based on the laterally derived polyphosphokinase 2 gene, showing improved resolution in differentiating different clades.

Metagenomic characterization of a novel non-ammonia-oxidizing Thaumarchaeota from hadal sediment

BackgroundThe hadal sediment, found at an ocean depth of more than 6000 m, is geographically isolated and under extremely high hydrostatic pressure, resulting in a unique ecosystem. Thaumarchaeota are ubiquitous marine microorganisms predominantly present in hadal environments. While there have been several studies on Thaumarchaeota there, most of them have primarily focused on ammonia-oxidizing archaea (AOA). However, systematic metagenomic research specifically targeting heterotrophic non-AOA Thaumarchaeota is lacking.ResultsIn this study, we explored the metagenomes of Challenger Deep hadal sediment, focusing on the Thaumarchaeota. Functional analysis of sequence reads revealed the potential contribution of Thaumarchaeota to recalcitrant dissolved organic matter degradation. Metagenome assembly binned one new group of hadal sediment-specific and ubiquitously distributed non-AOA Thaumarchaeota, named Group-3.unk. Pathway reconstruction of this new type of Thaumarchaeota also supports heterotrophic characteristics of Group-3.unk, along with ABC transporters for the uptake of amino acids and carbohydrates and catabolic utilization of these substrates. This new clade of Thaumarchaeota also contains aerobic oxidation of carbon monoxide-related genes. Complete glyoxylate cycle is a distinctive feature of this clade in supplying intermediates of anabolic pathways. The pan-genomic and metabolic analyses of metagenome-assembled genomes belonging to Group-3.unk Thaumarchaeota have highlighted distinctions, including the dihydroxy phthalate decarboxylase gene associated with the degradation of aromatic compounds and the absence of genes related to the synthesis of some types of vitamins compared to AOA. Notably, Group-3.unk shares a common feature with deep ocean AOA, characterized by their high hydrostatic pressure resistance, potentially associated with the presence of V-type ATP and di-myo-inositol phosphate syntheses-related genes. The enrichment of organic matter in hadal sediments might be attributed to the high recruitment of sequence reads of the Group-3.unk clade of heterotrophic Thaumarchaeota in the trench sediment. Evolutionary and genetic dynamic analyses suggest that Group-3 non-AOA consists of mesophilic Thaumarchaeota organisms. These results indicate a potential role in the transition from non-AOA to AOA Thaumarchaeota and from thermophilic to mesophilic Thaumarchaeota, shedding light on recent evolutionary pathways.ConclusionsOne novel clade of heterotrophic non-AOA Thaumarchaeota was identified through metagenome analysis of sediments from Challenger Deep. Our study provides insight into the ecology and genomic characteristics of the new sub-group of heterotrophic non-AOA Thaumarchaeota, thereby extending the knowledge of the evolution of Thaumarchaeota.5ZgcybfKsG2Wo48uk2DkGDVideo

In Silico Safety Assessment of Bacillus Isolated from Polish Bee Pollen and Bee Bread as Novel Probiotic Candidates.

Bacillus species isolated from Polish bee pollen (BP) and bee bread (BB) were characterized for in silico probiotic and safety attributes. A probiogenomics approach was used, and in-depth genomic analysis was performed using a wide array of bioinformatics tools to investigate the presence of virulence and antibiotic resistance properties, mobile genetic elements, and secondary metabolites. Functional annotation and Carbohydrate-Active enZYmes (CAZYme) profiling revealed the presence of genes and a repertoire of probiotics properties promoting enzymes. The isolates BB10.1, BP20.15 (isolated from bee bread), and PY2.3 (isolated from bee pollen) genome mining revealed the presence of several genes encoding acid, heat, cold, and other stress tolerance mechanisms, adhesion proteins required to survive and colonize harsh gastrointestinal environments, enzymes involved in the metabolism of dietary molecules, antioxidant activity, and genes associated with the synthesis of vitamins. In addition, genes responsible for the production of biogenic amines (BAs) and D-/L-lactate, hemolytic activity, and other toxic compounds were also analyzed. Pan-genome analyses were performed with 180 Bacillus subtilis and 204 Bacillus velezensis genomes to mine for any novel genes present in the genomes of our isolates. Moreover, all three isolates also consisted of gene clusters encoding secondary metabolites.

Comparative genomics and characterization of a multidrug-resistant Acinetobacter baumannii VRL-M19 isolated from a crowded setting in India

A crowded vegetable market serves as a mass gathering, posing a potential risk for infection transmission. In this study, we isolated a multidrug-resistant Acinetobacter baumannii strain, VRL-M19, from the air of such a market and conducted comparative genomics and phenotypic characterization. Antimicrobial susceptibility testing, genome sequencing using Illumina HiSeq X10, and pan-genome analysis with 788 clinical isolates identified core, accessory, and unique drug-resistant determinants. Mutational analysis of drug-resistance genes, virulence factor annotation, in vitro pathogenicity assessment, subsystem analysis, Multilocus sequence typing, and whole genome phylogenetic analysis were performed. VRL-M19 exhibited multidrug resistance with 69 determinants, and analysis across 788 clinical isolates and 350 Indian isolates revealed more accessory genes (52 out of 69) in the Indian isolates. Multiple mutations were observed in drug target modification genes, and the strain was identified as a moderate biofilm-former with 55 virulence factors. Whole genome phylogenetics indicated a close relationship between VRL-M19 and clinical A. baumannii strains. In conclusion, our comprehensive study suggests that VRL-M19 is a multidrug-resistant, potential pathogen with biofilm-forming capabilities, closely associated with clinical A. baumannii strains.

Combined transcriptomic and pangenomic analyses guide metabolic amelioration to enhance tiancimycins production.

Exploration of high-yield mechanism is important for further titer improvement of valuable antibiotics, but how to achieve this goal is challenging. Tiancimycins (TNMs) are anthraquinone-fused enediynes with promising drug development potentials, but their prospective applications are limited by low titers. This work aimed to explore the intrinsic high-yield mechanism in previously obtained TNMs high-producing strain Streptomyces sp. CB03234-S for the further titer amelioration of TNMs. First, the typical ribosomal RpsL(K43N) mutation in CB03234-S was validated to be merely responsible for the streptomycin resistance but not the titer improvement of TNMs. Subsequently, the combined transcriptomic, pan-genomic and KEGG analyses revealed that the significant changes in the carbon and amino acid metabolisms could reinforce the metabolic fluxes of key CoA precursors, and thus prompted the overproduction of TNMs in CB03234-S. Moreover, fatty acid metabolism was considered to exert adverse effects on the biosynthesis of TNMs by shunting and reducing the accumulation of CoA precursors. Therefore, different combinations of relevant genes were respectively overexpressed in CB03234-S to strengthen fatty acid degradation. The resulting mutants all showed the enhanced production of TNMs. Among them, the overexpression of fadD, a key gene responsible for the first step of fatty acid degradation, achieved the highest 21.7 ± 1.1mg/L TNMs with a 63.2% titer improvement. Our studies suggested that comprehensive bioinformatic analyses are effective to explore metabolic changes and guide rational metabolic reconstitution for further titer improvement of target products. KEY POINTS: • Comprehensive bioinformatic analyses effectively reveal primary metabolic changes. • Primary metabolic changes cause precursor enrichment to enhance TNMs production. • Strengthening of fatty acid degradation further improves the titer of TNMs.

Comparative pangenome analysis of Enterococcus faecium and Enterococcus lactis provides new insights into the adaptive evolution by horizontal gene acquisitions

BackgroundEnterococcus faecium and E. lactis are phylogenetically closely related lactic acid bacteria that are ubiquitous in nature and are known to be beneficial or pathogenic. Despite their considerable industrial and clinical importance, comprehensive studies on their evolutionary relationships and genomic, metabolic, and pathogenic traits are still lacking. Therefore, we conducted comparative pangenome analyses using all available dereplicated genomes of these species.ResultsE. faecium was divided into two subclades: subclade I, comprising strains derived from humans, animals, and food, and the more recent phylogenetic subclade II, consisting exclusively of human-derived strains. In contrast, E. lactis strains, isolated from diverse sources including foods, humans, animals, and the environment, did not display distinct clustering based on their isolation sources. Despite having similar metabolic features, noticeable genomic differences were observed between E. faecium subclades I and II, as well as E. lactis. Notably, E. faecium subclade II strains exhibited significantly larger genome sizes and higher gene counts compared to both E. faecium subclade I and E. lactis strains. Furthermore, they carried a higher abundance of antibiotic resistance, virulence, bacteriocin, and mobile element genes. Phylogenetic analysis of antibiotic resistance and virulence genes suggests that E. faecium subclade II strains likely acquired these genes through horizontal gene transfer, facilitating their effective adaptation in response to antibiotic use in humans.ConclusionsOur study offers valuable insights into the adaptive evolution of E. faecium strains, enabling their survival as pathogens in the human environment through horizontal gene acquisitions.

Pangenome Analysis of Helicobacter pylori Isolates from Selected Areas of Africa Indicated Diverse Antibiotic Resistance and Virulence Genes.

The challenge facing Helicobacter pylori (H. pylori) infection management in some parts of Africa is the evolution of drug-resistant species, the lack of gold standard in diagnostic methods, and the ineffectiveness of current vaccines against the bacteria. It is being established that even though clinical consequences linked to the bacteria vary geographically, there is rather a generic approach to treatment. This situation has remained problematic in the successful fight against the bacteria in parts of Africa. As a result, this study compared the genomes of selected H. pylori isolates from selected areas of Africa and evaluated their virulence and antibiotic drug resistance, those that are highly pathogenic and are associated with specific clinical outcomes and those that are less virulent and rarely associated with clinical outcomes. 146 genomes of H. pylori isolated from selected locations of Africa were sampled, and bioinformatic tools such as Abricate, CARD RGI, MLST, Prokka, Roary, Phandango, Google Sheets, and iTOLS were used to compare the isolates and their antibiotic resistance or susceptibility. Over 20 k virulence and AMR genes were observed. About 95% of the isolates were genetically diverse, 90% of the isolates harbored shell genes, and 50% harbored cloud and core genes. Some isolates did not retain the cagA and vacA genes. Clarithromycin, metronidazole, amoxicillin, and tinidazole were resistant to most AMR genes (vacA, cagA, oip, and bab). Conclusion. This study found both virulence and AMR genes in all H. pylori strains in all the selected geographies around Africa with differing quantities. MLST, Pangenome, and ORF analyses showed disparities among the isolates. This in general could imply diversities in terms of genetics, evolution, and protein production. Therefore, generic administration of antibiotics such as clarithromycin, amoxicillin, and erythromycin as treatment methods in the African subregion could be contributing to the spread of the bacterium's antibiotic resistance.

Whole genome sequencing distinguishes skin colonizing from infection-associated Cutibacterium acnes isolates.

Cutibacterium acnes can both be a helpful colonizer of the human skin as well as the causative agent of acne and purulent infections. Until today, it is a moot point whether there are C. acnes strains exclusively devoted to be part of the skin microbiome and others, that carry special features enabling them to cause disease. So far, the search for the molecular background of such diverse behavior has led to inconsistent results. In the present study, we prospectively collected C. acnes strains from 27 infected persons and 18 healthy controls employing rigid selection criteria to ensure their role as infectious agent or colonizer. The genome sequences from these strains were obtained and carefully controlled for quality. Deduced traditional phylotyping assigned almost all superficial isolates to type IA1, while the clinical strains were evenly distributed between types IA1, IB, and II. Single locus sequence typing (SLST) showed a predominance of A1 type for the control strains, whereas 56% of the clinical isolates belonged to types A1, H1 and K8. Pangenome analysis from all the present strains and 30 published genomes indicated the presence of an open pangenome. Except for three isolates, the colonizing strains clustered in clades separate from the majority of clinical strains, while 4 clinical strains clustered with the control strains. Identical results were obtained by a single nucleotide polymorphism (SNP) analysis. However, there were no significant differences in virulence gene contents in both groups. Genome-wide association studies (GWAS) from both the pangenome and SNP data consistently showed genomic differences between both groups located in metabolic pathway and DNA repair genes. Thus, the different behavior of colonizing and infectious C. acnes strains could be due to special metabolic capacities or flexibilities rather than specific virulence traits.

WGS Analysis of Staphylococcus warneri Outbreak in a Neonatal Intensive Care Unit.

Staphylococcus warneri is an opportunistic pathogen responsible for hospital-acquired infections (HAIs). The aim of this study was to describe an outbreak caused by S. warneri infection in a neonatal intensive care unit (NICU) and provide investigation, prevention and control strategies for this outbreak. We conducted an epidemiological investigation of the NICU S. warneri outbreak, involving seven neonates, staff, and environmental screening, to identify the source of infection. WGS analyses were performed on S. warneri isolates, including species identification, core genome single-nucleotide polymorphism (cgSNP) analysis, pan-genome analysis, and genetic characterization assessment of the prevalence of specific antibiotic resistance and virulence genes. Eight S. warneri strains were isolated from this outbreak, with seven from neonates and one from environment. Six clinical cases within three days in 2021 were linked to one strain isolated from environmental samples; isolates varied by 0-69 SNPs and were confirmed to be from an outbreak through WGS. Multiple infection prevention measures were implemented, including comprehensive environmental disinfection and stringent protocols, and all affected neonates were transferred to the isolation wards. Following these interventions, no further cases of S. warneri infections were observed. Furthermore, pan-genome analysis results suggested that in human S. warneri may exhibit host specificity. The investigation has revealed that the outbreak was linked to the milk preparation workbench by the WGS. It is recommended that there be a stronger focus on environmental disinfection management in order to raise awareness, improve identification, and prevention of healthcare-associated infections that are associated with the hospital environment.

A graph-based pan-genome of Brassica oleracea provides new insights into its domestication and morphotype diversification

The domestication of Brassica oleracea has resulted in diverse morphological types with distinct patterns of organ development. Here we report a graph-based pan-genome of B. oleracea constructed from high-quality genome assemblies of different morphotypes. The pan-genome harbors over 200 structural variant hotspot regions enriched in auxin- and flowering-related genes. Population genomic analyses revealed that early domestication of B. oleracea focused on leaf or stem development. Gene flows resulting from agricultural practices and variety improvement were detected among different morphotypes. Selective-sweep and pan-genome analyses identified an auxin-responsive small auxin up-regulated RNA gene and a CLAVATA3/ESR-RELATED family gene as crucial players in leaf–stem differentiation during the early stage of B. oleracea domestication and the BoKAN1 gene as instrumental in shaping the leafy heads of cabbage and Brussels sprouts. Our pan-genome and functional analyses further revealed that variations in the BoFLC2 gene play key roles in the divergence of vernalization and flowering characteristics among different morphotypes, and variations in the first intron of BoFLC3 are involved in fine-tuning the flowering process in cauliflower. This study provides a comprehensive understanding of the pan-genome of B. oleracea and sheds light on the domestication and differential organ development of this globally important crop species.

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

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. 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. 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. 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.

Comparative Genomics and Characterization of Shigella flexneri Isolated from Urban Wastewater.

Shigella species are a group of highly transmissible Gram-negative pathogens. Increasing reports of infection with extensively drug-resistant varieties of this stomach bug has convinced the World Health Organization to prioritize Shigella for novel therapeutic interventions. We herein coupled the whole-genome sequencing of a natural isolate of Shigella flexneri with a pangenome ana-lysis to characterize pathogen genomics within this species, which will provide us with an insight into its existing genomic diversity and highlight the root causes behind the emergence of quick vaccine escape variants. The isolated novel strain of S. flexneri contained ~4,500 protein-coding genes, 57 of which imparted resistance to antibiotics. A comparative pan-genomic ana-lysis revealed genomic variability of ~64%, the shared conservation of core genes in central metabolic processes, and the enrichment of unique/accessory genes in virulence and defense mechanisms that contributed to much of the observed antimicrobial resistance (AMR). A pathway ana-lysis of the core genome mapped 22 genes to 2 antimicrobial resistance pathways, with the bulk coding for multidrug efflux pumps and two component regulatory systems that are considered to work synergistically towards the development of resistance phenotypes. The prospective evolvability of Shigella species as witnessed by the marked difference in genomic content, the strain-specific essentiality of unique/accessory genes, and the inclusion of a potent resistance mechanism within the core genome, strengthens the possibility of novel serotypes emerging in the near future and emphasizes the importance of tracking down genomic diversity in drug/vaccine design and AMR governance.

Genomic and functional diversity of the human-derived isolates of Faecalibacterium.

Faecalibacterium is one of the most abundant bacteria in the gut microbiota of healthy adults, highly regarded as a next-generation probiotic. However, the functions of Faecalibacterium genomes from cultured strains and the distribution of different species in populations may differ among different sources. We here performed an extensive analysis of pan-genomes, functions, and safety evaluation of 136 Faecalibacterium genomes collected from 10 countries. The genomes are clustered into 11 clusters, with only five of them were characterized and validly nomenclated. Over 80% of the accessory genes and unique genes of Faecalibacterium are found with unknown function, which reflects the importance of expanding the collection of Faecalibacterium strains. All the genomes have the potential to produce acetic acid and butyric acid. Nine clusters of Faecalibacterium are found significantly enriched in the healthy individuals compared with patients with type II diabetes.. This study provides a comprehensive view of genomic characteristic and functions and of culturable Faecalibacterium bacterium from human gut, and enables clinical advances in the future.

Large-scale Pan Genomic Analysis of Mycobacterium tuberculosis Reveals Key Insights Into Molecular Evolutionary Rate of Specific Processes and Functions.

Mycobacterium tuberculosis (Mtb) is the causative agent of tuberculosis (TB), an infectious disease that is a major killer worldwide. Due to selection pressure caused by the use of antibacterial drugs, Mtb is characterised by mutational events that have given rise to multi drug resistant (MDR) and extensively drug resistant (XDR) phenotypes. The rate at which mutations occur is an important factor in the study of molecular evolution, and it helps understand gene evolution. Within the same species, different protein-coding genes evolve at different rates. To estimate the rates of molecular evolution of protein-coding genes, a commonly used parameter is the ratio dN/dS, where dN is the rate of non-synonymous substitutions and dS is the rate of synonymous substitutions. Here, we determined the estimated rates of molecular evolution of select biological processes and molecular functions across 264 strains of Mtb. We also investigated the molecular evolutionary rates of core genes of Mtb by computing the dN/dS values, and estimated the pan genome of the 264 strains of Mtb. Our results show that the cellular amino acid metabolic process and the kinase activity function evolve at a significantly higher rate, while the carbohydrate metabolic process evolves at a significantly lower rate for M. tuberculosis. These high rates of evolution correlate well with Mtb physiology and pathogenicity. We further propose that the core genome of M. tuberculosis likely experiences varying rates of molecular evolution which may drive an interplay between core genome and accessory genome during M. tuberculosis evolution.

Genomic comparison between Mycobacterium bovis and Mycobacterium microti and in silico analysis of peptide-based biomarkers for serodiagnosis.

In recent years, there has been an increase in the number of reported cases of Mycobacterium microti infection in various animals, which can interfere with the ante-mortem diagnosis of animal tuberculosis caused by Mycobacterium bovis. In this study, whole genome sequencing (WGS) was used to search for protein-coding genes to distinguish M. microti from M. bovis. In addition, the population structure of the available M. microti genomic WGS datasets is described, including three novel Belgian isolates from infections in alpacas. Candidate genes were identified by examining the presence of the regions of difference and by a pan-genome analysis of the available WGS data. A total of 80 genes showed presence-absence variation between the two species, including genes encoding Proline-Glutamate (PE), Proline-Proline-Glutamate (PPE), and Polymorphic GC-Rich Sequence (PE-PGRS) proteins involved in virulence and host interaction. Filtering based on predicted subcellular localization, sequence homology and predicted antigenicity resulted in 28 proteins out of 80 that were predicted to be potential antigens. As synthetic peptides are less costly and variable than recombinant proteins, an in silico approach was performed to identify linear and discontinuous B-cell epitopes in the selected proteins. From the 28 proteins, 157 B-cell epitope-based peptides were identified that discriminated between M. bovis and M. microti species. Although confirmation by in vitro testing is still required, these candidate synthetic peptides containing B-cell epitopes could potentially be used in serological tests to differentiate cases of M. bovis from M. microti infection, thus reducing misdiagnosis in animal tuberculosis surveillance.

GastronOmics: Edibility and safety of mycelium of the oyster mushroom Pleurotusostreatus

Food production is one of the most environmentally damaging human activities. In the face of climate change, it is essential to rethink our dietary habits and explore potential alternative foods catering both towards human and planetary needs. Fungal mycelium might be an attractive alternative protein source due to its rapid growth on sustainable substrates as well as promising nutritional and organoleptic properties. The natural biodiversity of filamentous fungi is vast and represents an untapped reservoir for food innovation. However, fungi are known to produce bioactive compounds that may affect human health, both positively and negatively. To narrow the search for safe and culinarily attractive fungal species, mycelia of edible fruiting-body forming fungi provide a promising starting point. Here, we explore whether the culinary attractiveness and safety of the commonly eaten mushroom, Pleurotus ostreatus, can also be translated to its mycelium.Whole-genome sequencing and pan-genome analysis revealed a high degree of genetic variability within the genus Pleurotus, suggesting that gastronomic traits as well as food safety may differ between strains. A representative strain, P. ostreatus M2191, was further analyzed for the food safety, nutritional properties and culinary applicability of its mycelium. No regulated mycotoxins were detected in either the fruiting body nor the mycelium. Yet, P. ostreatus is known to produce four peptide toxins, Ostreatin, Ostreolysin and Pleurotoysin A/B. These were found to be lower in the mycelium compared to fruiting bodies, which are already considered safe for consumption. Instead, a number of secondary metabolites with potential health benefits were detected in the fungal mycelium. In silico analysis of the proteome suggested low allergenicity. In addition, the fruiting body and the mycelium showed similar nutritional value, which was dependent on the growth substrate.To highlight the culinary potential of mycelium, we created a dish served at the two-star restaurant the Alchemist in Copenhagen, Denmark. Sensory analysis of the mycelium dish by an untrained consumer panel indicated consumer liking and openness to fungal mycelia.Based on sustainability, safety, culinary potential, and consumer acceptance, our findings suggest that P. ostreatus mycelium has great potential for use as a novel food source.

Pan-Genomic Analysis Identifies the Chinese Strain as a New Subspecies of Xanthomonas fragariae.

Xanthomonas fragariae is classified as a quarantine pathogen by the European and Mediterranean Plant Protection Organization. It commonly induces typical angular leaf spot (ALS) symptoms in strawberry leaves. X. fragariae strains from China (YL19, SHAQP01, and YLX21) exhibit ALS symptoms in leaves and more severe symptoms of dry cavity rot in strawberry crowns. Conversely, strains from other countries do not cause severe dry cavity rot symptoms in strawberries. After employing multilocus sequence analysis (MLSA), average nucleotide identity (ANI), and amino acid identity (AAI), we determined that Chinese strains of X. fragariae are genetically distinct from other strains and can be considered a new subspecies. Subsequent analysis of 63 X. fragariae genomes published at NCBI using IPGA and EDGAR3.0 revealed the pan-genomic profile, with 1,680 shared genes present in all 63 strains, including 71 virulence-related genes. Additionally, we identified 123 genes exclusive to all the Chinese strains, encompassing 12 virulence-related genes. The qRT-PCR analysis demonstrated that the expression of XopD, XopG1, CE8, GT2, and GH121 out of 12 virulence-related genes of Chinese strains (YL19) exhibited a constant increase in the early stages (6, 24, 54, and 96 hours postinoculation [hpi]) of strawberry leaf infected by YL19. So, the presence of XopD, XopG1, CE8, GT2, and GH121 in Chinese strains may play important roles in the early infection process of Chinese strains. These findings offer novel insights into comprehending the population structure and variation in the pathogenic capacity of X. fragariae.

Deciphering the genomic character of the multidrug-resistant <i>Staphylococcus aureus</i> from Dhaka, Bangladesh

<p><italic>Staphylococcus aureus</italic> is one of the leading agents of nosocomial and community-acquired infections. In this study, we explored the genomic characterization of eight methicillin-resistant clinical isolates of <italic>S. aureus</italic> from Dhaka, Bangladesh. Notably, all strains were resistant to penicillin, cephalosporins, and monobactams, with partial susceptibility to meropenem and complete susceptibility to amikacin, vancomycin, and tigecycline antibiotics. The strains were found to have an average genome size of 2.73 Mbp and an average of 32.64% GC content. Multi-locus sequence typing analysis characterized the most predominant sequence type as ST361, which belongs to the clonal complex CC361. All isolates harbored the <italic>mecA</italic> gene, often linked to SCC<italic>mec</italic>_type IV variants. Multidrug resistance was attributed to efflux pumps NorA, NorC, SdrM, and LmrS alongside genes encoding beta-lactamase BlaZ and factors like ErmC and MepA. Additionally, virulence factors including <italic>adsA</italic>, <italic>sdrC</italic>, <italic>cap8D</italic>, <italic>harA</italic>, <italic>esaA</italic>, essC, <italic>isdB</italic>, <italic>geh</italic>, and <italic>lip</italic> were commonly identified. Furthermore, genes associated with heme uptake and clumping were present, highlighting their roles in <italic>S. aureus</italic> colonization and pathogenesis. Nine secondary metabolite biosynthetic gene clusters were found, of which six were common in all the strains. Numerous toxin-antitoxin systems were predicted, with ParE and ParB-like nuclease domains found to be the most prevalent toxin and antitoxin, respectively. Pan-genome analysis revealed 2007 core genes and 229 unique genes in the studied strains. Finally, the phylogenomic analysis showed that most Bangladeshi strains were grouped into two unique clades. This study provides a genomic and comparative insight into the multidrug resistance and pathogenicity of <italic>S. aureus</italic> strains, which will play a crucial role in the future antibiotic stewardship of Bangladesh.</p>

Graph-based pan-genome analysis reveals diversity of structural variations in native and commercial chicken

Graph-based pan-genome analysis reveals diversity of structural variations in native and commercial chicken

Whole-genome comparison using complete genomes from Campylobacter fetus strains revealed single nucleotide polymorphisms on non-genomic islands for subspecies differentiation.

Bovine Genital Campylobacteriosis (BGC), caused by Campylobacter fetus subsp. venerealis, is a sexually transmitted bacterium that significantly impacts cattle reproductive performance. However, current detection methods lack consistency and reliability due to the close genetic similarity between C. fetus subsp. venerealis and C. fetus subsp. fetus. Therefore, this study aimed to utilize complete genome analysis to distinguish genetic features between C. fetus subsp. venerealis and other subspecies, thereby enhancing BGC detection for routine screening and epidemiological studies. This study reported the complete genomes of four C. fetus subsp. fetus and five C. fetus subsp. venerealis, sequenced using long-read sequencing technologies. Comparative whole-genome analyses (n = 25) were conducted, incorporating an additional 16 complete C. fetus genomes from the NCBI database, to investigate the genomic differences between these two closely related C. fetus subspecies. Pan-genomic analyses revealed a core genome consisting of 1,561 genes and an accessory pangenome of 1,064 genes between the two C. fetus subspecies. However, no unique predicted genes were identified in either subspecies. Nonetheless, whole-genome single nucleotide polymorphisms (SNPs) analysis identified 289 SNPs unique to one or the C. fetus subspecies. After the removal of SNPs located on putative genomic islands, recombination sites, and those causing synonymous amino acid changes, the remaining 184 SNPs were functionally annotated. Candidate SNPs that were annotated with the KEGG "Peptidoglycan Biosynthesis" pathway were recruited for further analysis due to their potential association with the glycine intolerance characteristic of C. fetus subsp. venerealis and its biovar variant. Verification with 58 annotated C. fetus genomes, both complete and incomplete, from RefSeq, successfully classified these seven SNPs into two groups, aligning with their phenotypic identification as CFF (Campylobacter fetus subsp. fetus) or CFV/CFVi (Campylobacter fetus subsp. venerealis and its biovar variant). Furthermore, we demonstrated the application of mraY SNPs for detecting C. fetus subspecies using a quantitative PCR assay. Our results highlighted the high genetic stability of C. fetus subspecies. Nevertheless, Campylobacter fetus subsp. venerealis and its biovar variants encoded common SNPs in genes related to glycine intolerance, which differentiates them from C. fetus subsp. fetus. This discovery highlights the potential of employing a multiple-SNP assay for the precise differentiation of C. fetus subspecies.