Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation

Genetic variation and evolutionary characteristics of Echovirus 11: new variant within genotype D5 associated with neonatal death found in China

Echovirus 11 (E11) has gained attention owing to its association with severe neonatal infections. From 2018 to 2023, a surge in severe neonatal cases and fatalities linked to a novel variant of genotype D5 was documented in China, France, and Italy. However, the prevention and control of E11 variants have been hampered by limited background data on the virus circulation and genetic variance. Therefore, the present study investigated the circulating dynamics of E11 and the genetic variation and molecular evolution of genotype D5 through the collection of strains from the national acute flaccid paralysis (AFP) and hand, foot, and mouth disease (HFMD) surveillance system in China during 2000-2022 and genetic sequences published in the GenBank database. The results of this study revealed a prevalent dynamic of E11 circulation, with D5 being the predominant genotype worldwide. Further phylogenetic analysis of genotype D5 indicated that it could be subdivided into three important geographic clusters (D5-CHN1: 2014-2019, D5-CHN2: 2016-2022, and D5-EUR: 2022-2023). Additionally, variant-specific (144) amino acid mutation sites and positive-selection pressure sites (132, 262) were identified in the VP1 region. Cluster-specific recombination patterns were also identified, with CVB5, E6, and CVB4 as the major recombinant viruses. These findings provide a preliminary landscape of E11 circulation worldwide and basic scientific data for further study of the pathogenicity of E11 variants.

RES-Xre toxin-antitoxin locus knaAT maintains the stability of the virulence plasmid in Klebsiella pneumoniae

Hypervirulent Klebsiella pneumoniae isolates have been increasingly reported worldwide, especially hypervirulent drug-resistant variants owing to the acquisition of a mobilizable virulence plasmid by a carbapenem-resistant strain. This pLVPK-like mobilizable plasmid encodes various virulence factors; however, information about its genetic stability is lacking. This study aimed to investigate the type II toxin-antitoxin (TA) modules that facilitate the virulence plasmid to remain stable in K. pneumoniae. More than 3,000 TA loci in 2,000 K. pneumoniae plasmids were examined for their relationship with plasmid cargo genes. TA loci from the RES-Xre family were highly correlated with virulence plasmids of hypervirulent K. pneumoniae. Overexpression of the RES toxin KnaT, encoded by the virulence plasmid-carrying RES-Xre locus knaAT, halts the cell growth of K. pneumoniae and E. coli, whereas co-expression of the cognate Xre antitoxin KnaA neutralizes the toxicity of KnaT. knaA and knaT were co-transcribed, representing the characteristics of a type II TA module. The knaAT deletion mutation gradually lost its virulence plasmid in K. pneumoniae, whereas the stability of the plasmid in E. coli was enhanced by adding knaAT, which revealed that the knaAT operon maintained the genetic stability of the large virulence plasmid in K. pneumoniae. String tests and mouse lethality assays subsequently confirmed that a loss of the virulence plasmid resulted in reduced pathogenicity of K. pneumoniae. These findings provide important insights into the role of the RES-Xre TA pair in stabilizing virulence plasmids and disseminating virulence genes in K. pneumoniae.

High-yield synthesis of 2-O-α-D-glucosyl-D-glycerate by a bifunctional glycoside phosphorylase

Osmolytes are produced by various microorganisms as a defense mechanism to protect cells and macromolecules from damage caused by external stresses in harsh environments. Due to their useful stabilizing properties, these molecules are applied as active ingredients in a wide range of cosmetics and healthcare products. The metabolic pathways and biocatalytic syntheses of glycosidic osmolytes such as 2-O-α-D-glucosyl-D-glycerate often involve the action of a glycoside phosphorylase. Here, we report the discovery of a glucosylglycerate phosphorylase from carbohydrate-active enzyme family GH13 that is also active on sucrose, which contrasts the strict specificity of known glucosylglycerate phosphorylases that can only use α-D-glucose 1-phosphate as glycosyl donor in transglycosylation reactions. The novel enzyme can be distinguished from other phosphorylases from the same family by the presence of an atypical conserved sequence motif at specificity-determining positions in the active site. The promiscuity of the sucrose-active glucosylglycerate phosphorylase can be exploited for the high-yielding and rapid synthesis of 2-O-α-D-glucosyl-D-glycerate from sucrose and D-glycerate. KEY POINTS: • A Xylanimonas protaetiae glycoside phosphorylase can use both d-glycerate and fructose as glucosyl acceptor with high catalytic efficiency • Biocatalytic synthesis of the osmolyte 2-O-α-d-glucosyl-d-glycerate • Positions in the active site of GH13 phosphorylases act as convenient specificity fingerprints.

Molecular characterization of the meq oncogene of Marek's disease virus in vaccinated Brazilian poultry farms reveals selective pressure on prevalent strains

Marek's disease virus (MDV) has become an increasingly virulent pathogen in the poultry industry despite vaccination efforts to control it. Brazil has experienced a significant rise of Marek's disease (MD) outbreaks in recent years. Our study aimed to analyze the complete meq gene sequences to understand the molecular epidemiological basis of MD outbreaks in Brazilian vaccinated layer farms. We detected a high incidence rate of visceral MD (67.74%) and multiple circulating MDV strains. The most prevalent and geographically widespread genotype presented several clinical and molecular characteristics of a highly virulent strain and evolving under positive selective pressure. Phylogenetic and phylogeographic analysis revealed a closer relationship with strains from the USA and Japan. This study sheds light on the circulation of MDV strains capable of infecting vaccinated birds. We emphasize the urgency of adopting preventive measures to manage MDV outbreaks threatening the poultry farming industry.

The optimal metric for viral genome space

Understanding the structural similarity between genomes is pivotal in classification and phylogenetic analysis. As the number of known genomes rockets, alignment-free methods have gained considerable attention. Among these methods, the natural vector method stands out as it represents sequences as vectors using statistical moments, enabling effective clustering based on families in biological taxonomy. However, determining an optimal metric that combines different elements in natural vectors remains challenging due to the absence of a rigorous theoretical framework for weighting different k-mers and orders. In this study, we address this challenge by transforming the determination of optimal weights into an optimization problem and resolving it through gradient-based techniques. Our experimental results underscore the substantial improvement in classification accuracy achieved by employing these optimal weights, reaching an impressive 92.73% on the testing set, surpassing other alignment-free methods. On one hand, our method offers an outstanding metric for virus classification, and on the other hand, it provides valuable insights into feature integration within alignment-free methods.

Antimicrobial susceptibility and resistome of Actinobacillus pleuropneumoniae in Taiwan: a next-generation sequencing analysis

Actinobacillus pleuropneumoniae infection causes a high mortality rate in porcine animals. Antimicrobial resistance poses global threats to public health. The current study aimed to determine the antimicrobial susceptibilities and probe the resistome of A. pleuropneumoniae in Taiwan. Herein, 133 isolates were retrospectively collected; upon initial screening, 38 samples were subjected to next-generation sequencing (NGS). Over the period 2017-2022, the lowest frequencies of resistant isolates were found for ceftiofur, cephalexin, cephalothin, and enrofloxacin, while the highest frequencies of resistant isolates were found for oxytetracycline, streptomycin, doxycycline, ampicillin, amoxicillin, kanamycin, and florfenicol. Furthermore, most isolates (71.4%) showed multiple drug resistance. NGS-based resistome analysis revealed aminoglycoside- and tetracycline-related genes at the highest prevalence, followed by genes related to beta-lactam, sulfamethoxazole, florphenicol, and macrolide. A plasmid replicon (repUS47) and insertion sequences (IS10R and ISVAp11) were identified in resistant isolates. Notably, the multiple resistance roles of the insertion sequence IS10R were widely proposed in human medicine; however, this is the first time IS10R has been reported in veterinary medicine. Concordance analysis revealed a high consistency of phenotypic and genotypic susceptibility to florphenicol, tilmicosin, doxycycline, and oxytetracycline. The current study reports the antimicrobial characterization of A. pleuropneumoniae for the first time in Taiwan using NGS.

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

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

Evaluation of outbreak persistence caused by multidrug-resistant and echinocandin-resistant Candida parapsilosis using multidimensional experimental and epidemiological approaches

Candida parapsilosis is known to cause severe and persistent outbreaks in clinical settings. Patients infected with multidrug-resistant C. parapsilosis (MDR Cp) isolates were identified in a large Turkish hospital from 2017-2020. We subsequently identified three additional patients infected with MDR Cp isolates in 2022 from the same hospital and two echinocandin-resistant (ECR) isolates from a single patient in another hospital. The increasing number of MDR and ECR isolates contradicts the general principle that the severe fitness cost associated with these phenotypes could prevent their dominance in clinical settings. Here, we employed a multidimensional approach to systematically assess the fitness costs of MDR and ECR C. parapsilosis isolates. Whole-genome sequencing revealed a novel MDR genotype infecting two patients in 2022. Despite severe in vitro defects, the levels and tolerances of the biofilms of our ECR and MDR isolates were generally comparable to those of susceptible wild-type isolates. Surprisingly, the MDR and ECR isolates showed major alterations in their cell wall components, and some of the MDR isolates consistently displayed increased tolerance to the fungicidal activities of primary human neutrophils and were more immunoevasive during exposure to primary human macrophages. Our systemic infection mouse model showed that MDR and ECR C. parapsilosis isolates had comparable fungal burden in most organs relative to susceptible isolates. Overall, we observed a notable increase in the genotypic diversity and frequency of MDR isolates and identified MDR and ECR isolates potentially capable of causing persistent outbreaks in the future.

Elucidation of critical chemical moieties of metallo-β-lactamase inhibitors and prioritisation of target metallo-β-lactamases

The urgent demand for effective countermeasures against metallo-β-lactamases (MBLs) necessitates development of novel metallo-β-lactamase inhibitors (MBLIs). This study is dedicated to identifying critical chemical moieties within previously developed MBLIs, and critical MBLs should serve as the target in MBLI evaluations. Using the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA), a systematic literature analysis was conducted, and the NCBI RefSeq genome database was exploited to access the abundance profile and taxonomic distribution of MBLs and their variant types. Through the implementation of two distinct systematic approaches, we elucidated critical chemical moieties of MBLIs, providing pivotal information for rational drug design. We also prioritised MBLs and their variant types, highlighting the imperative need for comprehensive testing to ensure the potency and efficacy of the newly developed MBLIs. This approach contributes valuable information to advance the field of antimicrobial drug discovery.

Genetic characteristics and potential pathogenic agents in Campylobacter upsaliensis based on genomic analysis

Campylobacter upsaliensis was the most common Campylobacter species in pets' gastrointestinal tracts and has been isolated from patients with bacteremia, hemolytic-uremic syndrome, spontaneous abortion, and Guillain-Barré syndrome. However, the genetic characteristics and the full extent of its significance as a human pathogen remain to be fully understood. This study involved an investigation for genomic analysis of 154 strains from different sources and additional antimicrobial resistance profiles of 26 strains for this species. The genomes contained 1,558-1,971 CDS and the genome sizes were estimated to vary from 1.53 Mb to 1.86 Mb, with an average GC content of  34.71%. The entire analyzed genomes could be divided into three clades (A, B, and C) based on ANI and phylogenomic analysis. Significantly, nearly all strains in Clade B were isolated from patient samples, and the virulence-related sequences FlgD, GmhA, and CdtC might serve as determining factors for the classification of Clade B. Half of the tested isolates had MIC values over 64 μg mL-1 for nalidixic acid, gentamicin, and streptomycin. Isolates from pets in China carried more resistant elements in the genomes. This study both provided a comprehensive profile of C. upsaliensis for its genomic features and suggested some pathogenic agents for human infection with this species.

Prevalence and characteristics of ertapenem-mono-resistant isolates among carbapenem-resistant Enterobacterales in China

Carbapenem-resistant Enterobacterales (CRE), specifically those resistant to only ertapenem among carbapenems (ETP-mono-resistant), are increasingly reported, while the optimal therapy options remain uncertain. To investigate the prevalence and characteristics of ETP-mono-resistant CRE, CRE strains were systematically collected from 102 hospitals across China between 2018 and 2021. A 1:1 randomized matching study was conducted with ETP-mono-resistant strains to meropenem- and/or imipenem-resistant (MEM/IPM-resistant) strains. Antimicrobial susceptibility testing, whole-genome sequencing, carbapenem-hydrolysing activity and the expression of carbapenemase genes were determined. In total, 18.8% of CRE strains were ETP-mono-resistant, with relatively low ertapenem MIC values. ETP-mono-resistant strains exhibited enhanced susceptibility to β-lactams, β-lactam/β-lactamase inhibitor combinations, levofloxacin, fosfomycin, amikacin and polymyxin than MEM/IPM-resistant strains (P < 0.05). Phylogenetic analysis revealed high genetic diversity among ETP-mono-resistant strains. Extended-spectrum β-lactamases (ESBLs) and/or AmpC, as well as porin mutations, were identified as potential major mechanisms mediating ETP-mono-resistance, while the presence of carbapenemases was found to be the key factor distinguishing the carbapenem-resistant phenotypes between the two groups (P < 0.001). Compared with the MEM/IPM-resistant group, limited carbapenemase-producing CRE (CP-CRE) strains in the ETP-mono-resistant group showed a significantly lower prevalence of ESBLs and porin mutations, along with reduced expression of carbapenemase. Remarkably, spot assays combined with modified carbapenem inactivation method indicated that ETP-mono-resistant CP-CRE isolates grew at meropenem concentrations eightfold above their corresponding MIC values, accompanied by rapidly enhanced carbapenem-hydrolysing ability. These findings illustrate that ETP-mono-resistant CRE strains are relatively prevalent and that caution should be exercised when using meropenem alone for treatment. The detection of carbapenemase should be prioritized.

Microbial diversity and oil biodegradation potential of northern Barents Sea sediments

The Arctic, an essential ecosystem on Earth, is subject to pronounced anthropogenic pressures, most notable being the climate change and risks of crude oil pollution. As crucial elements of Arctic environments, benthic microbiomes are involved in climate-relevant biogeochemical cycles and hold the potential to remediate upcoming contamination. Yet, the Arctic benthic microbiomes are among the least explored biomes on the planet. Here we combined geochemical analyses, incubation experiments, and microbial community profiling to detail the biogeography and biodegradation potential of Arctic sedimentary microbiomes in the northern Barents Sea. The results revealed a predominance of bacterial and archaea phyla typically found in the deep marine biosphere, such as Chloroflexi, Atribacteria, and Bathyarcheaota. The topmost benthic communities were spatially structured by sedimentary organic carbon, lacking a clear distinction among geographic regions. With increasing sediment depth, the community structure exhibited stratigraphic variability that could be correlated to redox geochemistry of sediments. The benthic microbiomes harbored multiple taxa capable of oxidizing hydrocarbons using aerobic and anaerobic pathways. Incubation of surface sediments with crude oil led to proliferation of several genera from the so-called rare biosphere. These include Alkalimarinus and Halioglobus, previously unrecognized as hydrocarbon-degrading genera, both harboring the full genetic potential for aerobic alkane oxidation. These findings increase our understanding of the taxonomic inventory and functional potential of unstudied benthic microbiomes in the Arctic.

Inoculum source determines the stress resistance of electroactive functional taxa in biofilms: A metagenomic perspective

The inoculum has a crucial impact on bioreactor initialization and performance. However, there is currently a lack of guidance on selecting appropriate inocula for applications in environmental biotechnology. In this study, we applied microbial electrolysis cells (MECs) as models to investigate the differences in the functional potential of electroactive microorganisms (EAMs) within anodic biofilms developed from four different inocula (natural or artificial), using shotgun metagenomic techniques. We specifically focused on extracellular electron transfer (EET) function and stress resistance, which affect the performance and stability of MECs. Community profiling revealed that the family Geobacteraceae was the key EAM taxon in all biofilms, with Geobacter as the dominant genus. The c-type cytochrome gene imcH showed universal importance for Geobacteraceae EET and was utilized as a marker gene to evaluate the EET potential of EAMs. Additionally, stress response functional genes were used to assess the stress resistance potential of Geobacter species. Comparative analysis of imcH gene abundance revealed that EAMs with comparable overall EET potential could be enriched from artificial and natural inocula (P > 0.05). However, quantification of stress response gene copy numbers in the genomes demonstrated that EAMs originating from natural inocula possessed superior stress resistance potential (196 vs. 163). Overall, this study provides novel perspectives on the inoculum effect in bioreactors and offers theoretical guidance for selecting inoculum in environmental engineering applications.

Weak purifying selection in allelic diversity of the ADSL gene in indigenous and local chicken breeds and red junglefowl in Thailand

High levels of purine and uric acid, which are associated with health issues such as gout and cardiovascular disease, are found in the meat of fast-growing broiler chickens, which raises concerns about the quality of chicken meat and the health of the consumers who consume it. High genetic homogeneity and uniformity, particularly in genes involved in the synthesis of inosine monophosphate (IMP) and subsequent process of purine synthesis, which are associated with the meat quality, are exhibited in commercial broiler chickens owing to intensive inbreeding programs. Adenosine succinate lyase (ADSL) is a key enzyme involved in de novo purine biosynthetic pathway and its genetic polymorphisms affect IMP metabolism and purine content. In this study, we investigated the polymorphism of the ADSL gene in indigenous and local chicken breeds and red junglefowl in Thailand, using metabarcoding and genetic diversity analyses. Five alleles with 73 single nucleotide polymorphisms in exon 2, including missense and silent mutations, which may act on the synthesis efficiency of IMP and purine. Their protein structures revealed changes in amino acid composition that may affect ADSL enzyme activity. Weak purifying selection in these ADSL alleles was observed in the chicken population studied, implying that the variants have minor fitness impacts and a greater probability of fixation of beneficial mutations than strong purifying selection. A potential selective sweep was observed in Mae Hong Son chickens, whose purine content was lower than that in other breeds. This suggests a potential correlation between variations of the ADSL gene and reduced purine content and an impact of ADSL expression on the quality of chicken meat. However, further studies are required to validate its potential availability as a genetic marker for selecting useful traits that are beneficial to human health and well-being.

An entropy-based study on the mutational landscape of SARS-CoV-2 in USA: Comparing different variants and revealing co-mutational behavior of proteins

COVID-19 emergency has pushed the international scientific community to use every resource to combat the spread of the virus, to understand its biology and predict its possible evolution in terms of new variants. Since the first SARS-CoV-2 virus nucleotide and amino acid sequences were made available, information theory was used to study how viral information content was changing over time and then trace the evolution of its mutational landscape. In this work we analyzed SARS-CoV-2 sequences collected mainly in the USA in a period from March 2020 until December 2022 and computed mutation profiles of viral proteins over time through an entropy-based approach using Shannon Entropy and Hellinger distance. This representation allows an at-a-glance view of the mutational landscape of viral proteins over time and can provide new insights on the evolution of the virus from different points of view. Non-structural proteins typically showed flat mutation profiles, characterized by a very low Average mutation Entropy, while accessory and structural proteins showed mostly non uniform and high mutation profiles, often coupled with the predominance of variants. Interestingly NSP2 protein, whose function is currently still debated, falls in the same branch of NSP14 and NSP10 in the phylogenetic tree of mutations constructed through correlations of mutation profiles, suggesting a co-evolution of those proteins and a possible functional link with each other. To the best of our knowledge this is the first study based on a massive amount of data (n = 107,939,973) that analyzes from an entropy point of view the mutational landscape of SARS-CoV-2 over time and depicts a mutational temporal profile of each protein of the virus.

Molecular characterization and phylogenetic analysis of babA gene of Helicobacter pylori isolated from Indian patients with gastrointestinal diseases

INTRODUCTION

Outer membrane protein (OMP) of Helicobacter pylori (H. pylori) i.e., blood group antigen binding adhesin (babA) is responsible for the attachment of H. pylori in the gastric epithelium. Its adherence is causative for gastric pathology such as gastritis, peptic ulcer disease (PUD), or digestive tract disorders like erosive reflux disease (ERD) and (NERD) non-erosive reflux disease and together called Gastroesophageal reflux disease (GERD). BabA manifests rapid and varied selection via substitution of amino acid in its Leb-carbohydrate binding domain (CBD) which enables better binding preferences for distinct human populations and ABO blood group phenotypes. The positive evolutionary selection of the pathogenic factor of this genetically diverse bacterium has enabled it to adapt to the host gastric environment. Analyzing the association of virulent genes (cagA, vacA) and babA will help us better understand bacteria's pathogenicity.

METHOD

109 H. pylori strains from patients with distinct gastrointestinal diseases were genotyped using Polymerase Chain Reaction(PCR) for cagA, vacA, and babA followed by Sanger sequencing and phylogenetic analysis.

RESULT

In the babA + ve genotype, a statistically significant association with p = 0.04 and < 0.0001 is seen in gastritis and ERD respectively. A significant association of genotype vacAs1m2 (p = 0.0002) was seen in gastritis, vacAs1m1 (p = 0.02) in NERD, vacAs1m1 (p < 0.0001) and vacAs1m2 (p = 0.002) in ERD. This relationship helps to detect gastritis or ERD where BabA gene can be used as an independent marker for detecting their presence.

CONCLUSION

The appearance of variants within distinct disease categories is due to local genetic variation.

Whole genome analysis reveals unique traits of SARS-CoV-2 in pediatric patients

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) continues to challenge the global healthcare with emerging variants and higher infectivity as well as morbidities. This study investigated potential age-related variations through genomic characterization of the virus under common clinical settings. A cohort comprising 71 SARS-CoV-2 strains from both infected infants and accompanying adults, diagnosed via RT-PCR at a tertiary pediatric hospital and research center, underwent Illumina paired-end sequencing. The subsequent analysis involved standard genomic screening, phylogeny construction, and mutational analyses. The analyzed SARSCoV- 2 strains were compared with globally circulating variants. The overall distribution revealed 67.61 % Delta, 25.7 % Omicron, and 1 % either Kappa or Alpha variants. In 2021, Delta predominated at ∼ 94 %, with Alpha/Kappa accounting for around 5 %. However, in 2022, over 94 % of the samples were Omicron variants, signifying a substantial shift from Delta dominance. Delta variants constituted 69.5 % of infections in adults and 78.5 % in infants, while Omicron variants were responsible for 31 % of cases in infants and 18 % in adults. The Spike region harbored the majority of mutations, with T19R being the most prevalent mutation in the Delta lineage. Notably, the frequencies of this mutation varied between infants and adults. In Omicron samples, G142D emerged as the most prevalent mutation. Our dataset predominantly featured clade 21A and lineage B.1.617.2. This study underscores the differential clinical presentations and genomic characteristics of SARS-CoV-2 in pediatric patients and accompanying adults. Understanding the dynamic evolution of the SARS- CoV-2 in both pediatric and adults can help in strengthening prophylactic measures.

Prevalence of opportunistic pathogens and anti-microbial resistance in urban aquaculture ponds

Bacterial antimicrobial resistance (AMR) has emerged as a significant concern worldwide. The microbial community profile and potential AMR level in aquaculture ponds are often undervalued and attract less attention than other aquatic environments. We used amplicon and metagenomic shotgun sequencing to study microbial communities and AMR in six freshwater polyculture ponds in rural and urban areas of Bangladesh. Amplicon sequencing revealed different community structures between rural and urban ponds, with urban ponds having a higher bacterial diversity and opportunistic pathogens including Streptococcus, Staphylococcus, and Corynebacterium. Despite proteobacterial dominance, Firmicutes was the most interactive in the community network, especially in the urban ponds. Metagenomes showed that drug resistance was the most common type of AMR found, while metal resistance was only observed in urban ponds. AMR and metal resistance genes were found mainly in beta and gamma-proteobacteria in urban ponds, while AMR was found primarily in alpha-proteobacteria in rural ponds. We identified potential pathogens with a high profile of AMR and metal resistance in urban aquaculture ponds. As these ponds provide a significant source of protein for humans, our results raise significant concerns for the environmental sustainability of this food source and the dissemination of AMR into the food chain.

The proteome of Nicotiana benthamiana is shaped by extensive protein processing

Processing by proteases irreversibly regulates the fate of plant proteins and hampers the production of recombinant proteins in plants, yet only few processing events have been described in agroinfiltrated Nicotiana benthamiana, which has emerged as the main transient protein expression platform in plant science and molecular pharming. Here, we used in-gel digests and mass spectrometry to monitor the migration and topography of 5040 plant proteins within a protein gel. By plotting the peptides over the gel slices, we generated peptographs that reveal where which part of each protein was detected within the protein gel. These data uncovered that 60% of the detected proteins have proteoforms that migrate at lower than predicted molecular weights, implicating extensive proteolytic processing. This analysis confirms the proteolytic removal and degradation of autoinhibitory prodomains of most but not all proteases, and revealed differential processing within pectinemethylesterase and lipase families. This analysis also uncovered intricate processing of glycosidases and uncovered that ectodomain shedding might be common for a diverse range of receptor-like kinases. Transient expression of double-tagged candidate proteins confirmed processing events in vivo. This large proteomic dataset implicates an elaborate proteolytic machinery shaping the proteome of N. benthamiana.

Syntrichia ruralis: emerging model moss genome reveals a conserved and previously unknown regulator of desiccation in flowering plants

Water scarcity, resulting from climate change, poses a significant threat to ecosystems. Syntrichia ruralis, a dryland desiccation-tolerant moss, provides valuable insights into survival of water-limited conditions. We sequenced the genome of S. ruralis, conducted transcriptomic analyses, and performed comparative genomic and transcriptomic analyses with existing genomes and transcriptomes, including with the close relative S. caninervis. We took a genetic approach to characterize the role of an S. ruralis transcription factor, identified in transcriptomic analyses, in Arabidopsis thaliana. The genome was assembled into 12 chromosomes encompassing 21 169 protein-coding genes. Comparative analysis revealed copy number and transcript abundance differences in known desiccation-associated gene families, and highlighted genome-level variation among species that may reflect adaptation to different habitats. A significant number of abscisic acid (ABA)-responsive genes were found to be negatively regulated by a MYB transcription factor (MYB55) that was upstream of the S. ruralis ortholog of ABA-insensitive 3 (ABI3). We determined that this conserved MYB transcription factor, uncharacterized in Arabidopsis, acts as a negative regulator of an ABA-dependent stress response in Arabidopsis. The new genomic resources from this emerging model moss offer novel insights into how plants regulate their responses to water deprivation.

The global genetic diversity of planktonic foraminifera reveals the structure of cryptic speciation in plankton

The nature and extent of diversity in the plankton has fascinated scientists for over a century. Initially, the discovery of many new species in the remarkably uniform and unstructured pelagic environment appeared to challenge the concept of ecological niches. Later, it became obvious that only a fraction of plankton diversity had been formally described, because plankton assemblages are dominated by understudied eukaryotic lineages with small size that lack clearly distinguishable morphological features. The high diversity of the plankton has been confirmed by comprehensive metabarcoding surveys, but interpretation of the underlying molecular taxonomies is hindered by insufficient integration of genetic diversity with morphological taxonomy and ecological observations. Here we use planktonic foraminifera as a study model and reveal the full extent of their genetic diversity and investigate geographical and ecological patterns in their distribution. To this end, we assembled a global data set of ~7600 ribosomal DNA sequences obtained from morphologically characterised individual foraminifera, established a robust molecular taxonomic framework for the observed diversity, and used it to query a global metabarcoding data set covering ~1700 samples with ~2.48 billion reads. This allowed us to extract and assign 1 million reads, enabling characterisation of the structure of the genetic diversity of the group across ~1100 oceanic stations worldwide. Our sampling revealed the existence of, at most, 94 distinct molecular operational taxonomic units (MOTUs) at a level of divergence indicative of biological species. The genetic diversity only doubles the number of formally described species identified by morphological features. Furthermore, we observed that the allocation of genetic diversity to morphospecies is uneven. Only 16 morphospecies disguise evolutionarily significant genetic diversity, and the proportion of morphospecies that show genetic diversity increases poleward. Finally, we observe that MOTUs have a narrower geographic distribution than morphospecies and that in some cases the MOTUs belonging to the same morphospecies (cryptic species) have different environmental preferences. Overall, our analysis reveals that even in the light of global genetic sampling, planktonic foraminifera diversity is modest and finite. However, the extent and structure of the cryptic diversity reveals that genetic diversification is decoupled from morphological diversification, hinting at different mechanisms acting at different levels of divergence.

PanVA: Pangenomic Variant Analysis

Genomics researchers increasingly use multiple reference genomes to comprehensively explore genetic variants underlying differences in detectable characteristics between organisms. Pangenomes allow for an efficient data representation of multiple related genomes and their associated metadata. However, current visual analysis approaches for exploring these complex genotype-phenotype relationships are often based on single reference approaches or lack adequate support for interpreting the variants in the genomic context with heterogeneous (meta)data. This design study introduces PanVA, a visual analytics design for pangenomic variant analysis developed with the active participation of genomics researchers. The design uniquely combines tailored visual representations with interactions such as sorting, grouping, and aggregation, allowing users to navigate and explore different perspectives on complex genotype-phenotype relations. Through evaluation in the context of plants and pathogen research, we show that PanVA helps researchers explore variants in genes and generate hypotheses about their role in phenotypic variation.

Identification of pine wood nematode (Bursaphelenchus xylophilus) loading response genes in Japanese pine sawyer (Monochamus alternatus) through comparative genomics and transcriptomics

BACKGROUND

Pine wood nematode (PWN; Bursaphelenchus xylophilus) is the causative agent of pine wilt disease (PWD), which is considered the most dangerous biohazard to conifer trees globally. The transmission of PWN relies on insect vectors, particularly the Japanese pine sawyer (JPS; Monochamus alternatus). However, the molecular mechanism underlying PWN-JPS assembly remains largely unknown.

RESULTS

Here, we found that both geographical and gender could significantly affect the PCA (PWN carrying amount) of JPS; thus, JPS transcriptomes from diverse locations and genders were explored regard to PWN loading. Due to the shortage of genomes, we developed a full-length reference transcriptome for analyzing next-generation sequencing data. A comparative genomic study was performed, and 11 248 potential PWN-carrying associate genes (β) were nominated in JPS by using the reported genomes of PWN and non-PWN carrier insect species. Then, 151 differentially expressed transcripts (DETs), 28 of them overlapped with β, correlated with the PCA of JPS were nominated by RNA-Seq, and found that fatty acid β-oxidation might be the key factor that affected the PCA of JPS. Furthermore, JPS fatty acid β-oxidation rates were experimentally decreased using the inhibitor Etomoxir, leading to an increased PCA of JPS. Meanwhile, silencing MaCPT1 in JPS by RNA interference led to a decreased fatty acid β-oxidation rate and increased PCA of JPS.

CONCLUSIONS

In conclusion, MaCPT1 was able to decrease the PWN-JPS assembly formation through the fatty acid β-oxidation of JPS. These results provide new insights for exploring the impact of PWN invasion on JPS. © 2024 Society of Chemical Industry.

Bioinformatics-based identification of GH12 endoxyloglucanases in citrus-pathogenic Penicillium spp

Millions of tons of citrus peel waste are produced every year as a byproduct of the juice industry. Citrus peel is rich in pectin and xyloglucan, but while the pectin is extracted for use in the food industry, the xyloglucan is currently not valorized. To target hydrolytic degradation of citrus peel xyloglucan into oligosaccharides, we have used bioinformatics to identify three glycoside hydrolase 12 (GH12) endoxyloglucanases (EC 3.2.1.151) from the citrus fruit pathogens Penicillium italicum GL-Gan1 and Penicillium digitatum Pd1 and characterized them on xyloglucan obtained by alkaline extraction from citrus peel. The enzymes displayed pH-temperature optima of pH 4.6-5.3 and 35-37°C. PdGH12 from P. digitatum and PiGH12A from P. italicum share 84% sequence identity and displayed similar kinetics, although kcat was highest for PdGH12. In contrast, PiGH12B from P. italicum, which has the otherwise conserved Trp in subsite -4 replaced with a Tyr, displayed a 3 times higher KM and a 4 times lower kcat/KM than PiGH12A, but was the most thermostable enzyme of the three Penicillium-derived endoxyloglucanases. The benchmark enzyme AnGH12 from Aspergillus nidulans was more thermally stable and had a higher pH-temperature optimum than the enzymes from Penicillum spp. The difference in structure of the xyloglucan oligosaccharides extracted from citrus peel xyloglucan and tamarind xyloglucan by the new endoxyloglucanases was determined by LC-MS. The inclusion of citrus peel xyloglucan demonstrated that the endoxyloglucanases liberated fucosylated xyloglucan oligomers, implying that these enzymes have the potential to upgrade citrus peel residues to produce oligomers useful as intermediates or bioactive compounds.

Combined metagenomic- and culture-based approaches to investigate bacterial strain-level associations with medication-controlled mild-moderate atopic dermatitis

Background

The skin microbiome is disrupted in atopic dermatitis (AD). Existing research focuses on moderate to severe, unmedicated disease.

Objective

We sought to investigate metagenomic- and culture-based bacterial strain-level differences in mild, medicated AD and the effects these have on human keratinocytes (HKs).

Methods

Skin swabs from anterior forearms were collected from 20 pediatric participants (11 participants with AD sampled at lesional and nonlesional sites and 9 age- and sex-matched controls). Participants had primarily mild to moderate AD and maintained medication use. Samples were processed for microbial metagenomic sequencing and bacterial isolation. Isolates identified as Staphylococcus aureus were tested for enterotoxin production. HK cultures were treated with cell-free conditioned media from representative Staphylococcus species to measure barrier effects.

Results

Metagenomic sequencing identified significant differences in microbiome composition between AD and control groups. Differences were seen at the species and strain levels for Staphylococci, with S aureus found only in participants with AD and differences in Staphylococcus epidermidis strains between control and AD swabs. These strains showed differences in toxin gene presence, which was confirmed in vitro for S aureus enterotoxins. The strain from the participant with the most severe AD produced enterotoxin B levels more than 100-fold higher than the other strains (P < .001). Strains also displayed differential effects on HK metabolism and barrier function.

Conclusions

Strain-level differences in toxin genes from Staphylococcus strains may explain varying effects on HK, with S aureus and non-aureus strains negatively affecting viability and barrier function. These differences are likely important in AD pathogenesis.

Genome-wide characterization of Related to ABI3/VP1 transcription factors among U's triangle Brassica species reveals a negative role for BnaA06.RAV3L in seed size

The transcription factors Related to ABI3/VP1 (RAV) are crucial for various plant processes and stress responses. Although the U's triangle Brassica species genomes have been released, the knowledge regarding the RAV family is still limited. In this study, we identified 123 putative RAV genes across the six U's triangle Brassica species (Brassica rapa, 14; Brassica oleracea, 14; Brassica nigra, 13; Brassica carinata, 27; Brassica juncea, 28; Brassica napus, 27). Phylogenetic analysis categorized them into three groups. The RAV genes exhibited diversity in both functional and structural aspects, particularly in gene structure and cis-acting elements within their promoters. The expression analysis revealed that BnaRAV genes in Group 1/2 exhibited diverse expression patterns across various tissues, while those in Group 3 did not show expression except for BnaRAV3L-2 and BnaRAV3L-6, which were exclusively expressed in seeds. Furthermore, the seed-specific expression of BnaA06. RAV3L (BnaRAV3L-2) was confirmed through promoter-GUS staining. Subcellular localization studies demonstrated that BnaA06.RAV3L is localized to the nucleus. The overexpression of BnaA06. RAV3L in Arabidopsis led to a remarkable inhibition of seed-specific traits such as seed width, seed length, seed area, and seed weight. This study provides insights into the functional evolution of the RAV gene family in U triangle Brassica species. It establishes a foundation for uncovering the molecular mechanisms underlying the negative role of RAV3L in seed development.

Efficient activity screening of new glucuronoyl esterases using a pNP-based assay

Glucuronoyl esterases (CE15, EC 3.1.1.117) catalyze the hydrolysis of ester bonds between lignin and carbohydrates in lignocellulose. They are widespread within fungi and bacteria, and are subjects to research interest due to their potential applicability in lignocellulose processing. Identifying new and relevant glucuronoyl esterase candidates is challenging because available model substrates poorly represent the natural substrate, which leads to inefficient screening for the activity. In this study, we demonstrate how fifteen novel, fungal, putative glucuronoyl esterases from family CE15 were expressed and screened for activity towards a commercially available, colorimetric assay based on the methyl-ester of 4-O-methyl-aldotriuronic acid linked to para-nitrophenol (methyl ester-UX-β-pNP) and coupled with the activity of GH67 (α-glucuronidase) and GH43 (β-xylosidase) activity. The assay provides easy means for accurately establishing activity and determining specific activity of glucuronoyl esterases. Out of the fifteen expressed CE15 proteins, seven are active and were purified to determine their specific activity. The seven active enzymes originate from Auricularia subglabra (3 proteins), Ganoderma sinensis (2 proteins) and Neocallimastix californiae (2 proteins). Among the CE15 proteins not active towards the screening substrate (methyl ester-UX-β-pNP) were proteins originating from Schizophyllum commune, Podospora anserina, Trametes versicolor, and Coprinopsis cinerea. It is unexpected that CE15 proteins from such canonical lignocellulose degraders do not have the anticipated activity, and these observations call for deeper investigations.

Genomic and phenotypic characterization of Thermosynechococcus-like strains reveals eight species within the genus Thermosynechococcus and a novel genus Parathermosynechococcus gen. nov

Thermophilic unicellular cyanobacteria of the family Thermosynechococcaceae are essential primary producers and integral components of many microbial mats found in hot springs of Asia and North America. Historically, based on their simple morphology, these organisms, along with members of taxonomically unrelated thermophilic Thermostichaceae have been described with a generic term, "Synechococcus", used for elongated unicellular cyanobacteria. This has created significant misperception in the scientific literature regarding the taxonomic status of these essential thermophilic primary producers and their relationship with Synechococcus sensu stricto. In this manuscript, we attempted a genome-driven taxonomic reevaluation of the family Thermosynechococcaceae. Application of genomic analyses such as GTDB classification, ANI/AAI and phylogenomics support the delineation of eight species within genus Thermosynechococcus. Two subspecies were further identified within T. taiwanensis by dDDH and phylogenomics. Moreover, the results also suggest the presence of two putative new genera phylogenetically alongside genus Thermosynechococcus, a thermophilic genus Parathermosynechococcus represented by PCC 6715 and a non-thermophilic genus represented by PCC 6312. The proposed genospecies and new genera were further integrated with morphological and/or ecological information. Interestingly, the phylogeny of 16S-23S ITS achieved a better taxonomic relationship than that of 16S rRNA and supported the genome-based classification of Thermosynechococcus spp. Finally, the pan-genome analysis indicated a conserved pattern of genomic core among known members of Thermosynechococcus.

Microbial consortium involved in ferromanganese and francolite biomineralization in an anchialine environment (Zinzulùsa Cave, Castro, Italy)

Tidally-influenced subterranean settings represent natural geomicrobiological laboratories, relatively unexplored, that facilitate the investigation of new biomineralization processes. The unusual water chemistry of Zinzulùsa Cave and its oligotrophic and aphotic conditions have allowed the development of a unique ecosystem in which complex bacterial activities induce rare biomineralization processes. A diversified microbial community develops on centimeter-thick crusts that form in the submerged part of the cave. The crusts are formed of Ca-phosphate minerals, mostly carbonate-fluoroapatite (francolite), covered by a black crust, few microns in thickness, composed of ferromanganiferous oxides (hematite and vernadite). Diffuse coccoidal and filamentous bacteria and amorphous organic matter are mixed with the minerals. The micromorphologies and comparative 16S rRNA gene-based metabarcoding analyses identify a "core microbiota" also common to other natural environments characterized by FeMn and Ca-phosphate mineralization. The microbiota is characterized by nitrifying, sulfide/sulfur/thiosulfate-oxidizing and sulfate/thiosulfate/sulfur-reducing bacteria. In addition, manganese-oxidizing bacteria include the recently described "Ca. Manganitrophus noduliformans" and an abundance of bacteria belonging to the Planctomycetes-Verrucomicrobia-Chlamydiae (PVC) superphylum, as well as Haliangiales (fruiting body-forming bacteria) and Hyphomicrobiales (stalked and budding bacteria) that are known to produce extracellular polymers that trap iron and manganese oxides. 16S rRNA gene metabarcoding analysis showed the presence of bacteria able to utilize many organic P substrates, including Ramlibacter, and SEM images revealed traces of fossilized microorganisms resembling "cable bacteria", which may play a role in Ca-phosphate biomineralization. Overall, the data indicate biomineralization processes induced by microbial metabolic activities for both ferromanganiferous oxide and francolite components of these crusts.

Novel oceanic cyanobacterium isolated from Bangaram island with profound acid neutralizing ability is proposed as Leptolyngbya iicbica sp. nov. strain LK

An acid-neutralizing, filamentous, non-heterocytous, marine cyanobacterium named 'LK' has been isolated from the seashore of Bangaram Island, an atoll of Lakshadweep, India, and is described here as a novel species. LK has been characterized using morphological, ecological, and genomic features. Based on 16S rRNA, whole-genome sequencing, and marker gene-based analysis, LK has been identified as a new species. LK clustered with Leptolyngbya-like strains belonging to the LPP group but diverged from Leptolyngbya sensu stricto, indicating the polyphyletic nature of the Leptolyngbya genus. Leptolyngbya sp. SIOISBB and Halomicronema sp. CCY15110 were identified as LK's two closest phylogenetic neighbors in various phylogenetic studies. The analysis of 16S rRNA, ITS secondary structures, and genome relatedness indices such as AAI, ANI, and gANI strongly support LK as a novel species of the Leptolyngbya genus. The mechanism behind acid neutralization in LK has been delineated, attributing it to a surface phenomenon most likely due to the presence of salts of calcium, magnesium, sodium, and potassium. We name LK as Leptolyngbya iicbica strain LK which is a novel species with prominent acidic pH-neutralizing properties.

A complete biosynthetic pathway of the long-chain polyunsaturated fatty acids in an amphidromous fish, ayu sweetfish Plecoglossus altivelis (Stomiati; Osmeriformes)

The biosynthetic capability of the long-chain polyunsaturated fatty acids (LC-PUFA) in teleosts are highly diversified due to evolutionary events such as gene loss and subsequent neo- and/or sub-functionalisation of enzymes encoded by existing genes. In the present study, we have comprehensively characterised genes potentially involved in LC-PUFA biosynthesis, namely one front-end desaturase (fads2) and eight fatty acid elongases (elovl1a, elovl1b, elovl4a, elovl4b, elovl5, elovl7, elovl8a and elovl8b) from an amphidromous teleost, Ayu sweetfish, Plecoglossus altivelis. Functional analysis confirmed Fads2 with Δ6, Δ5 and Δ8 desaturase activities towards multiple PUFA substrates and several Elovl enzymes exhibited elongation capacities towards C18-20 or C18-22 PUFA substrates. Consequently, P. altivelis possesses a complete enzymatic capability to synthesise physiologically important LC-PUFA including arachidonic acid (ARA, 20:4n-6), eicosapentaenoic acid (EPA, 20:5n-3) and docosahexaenoic acid (DHA, 22:6n-3) from their C18 precursors. Interestingly, the loss of elovl2 gene in P. altivelis was corroborated by genomic and phylogenetic analyses. However, this constraint would possibly be overcome by the function of alternative Elovl enzymes, such as Elovl1b, which has not hitherto been functionally characterised in teleosts. The present study contributes novel insights into LC-PUFA biosynthesis in the relatively understudied teleost group, Osmeriformes (Stomiati), thereby enhancing our understanding of the complement of LC-PUFA biosynthetic genes within teleosts.

Comprehensive bioinformatics-based annotation and functional characterization of bovine chymosin protein revealed novel biological insights

Chymosin, an aspartic protease present in the stomachs of young ruminants like cows (bovine), causes milk coagulation and cheese production through the breakdown of κ-casein peptide bonds at the Met105-Phe106 site. Bovine chymosin is first synthesized as a pre-prochymosin that is cleaved to produce the mature chymosin protein. Despite significant strides in research, our understanding of this crucial enzyme remains incomplete. The purpose of this work was to perform in silico evolutionary and functional analysis and to gain unique insights into the structure of this protein. For this, the sequence of Bos taurus chymosin from UniProt database was subjected to various bioinformatics analyses. We found that bovine chymosin is a low molecular weight and hydrophilic protein that has homologs in other Bovidae species. Two active sites of aspartic peptidases, along with a functional domain, were identified. Gene Ontology analysis further confirmed chymosin's involvement in proteolysis and aspartic endopeptidase activity. Potential disordered residues and post-translational modification sites were also uncovered. It was revealed that the secondary structure of bovine chymosin is comprised of beta strands (44.27%), coils (43.65%), and alpha helices (12.07%). A highly optimized 3D structure was also obtained. Moreover, crucial protein-protein interactions were unveiled. Altogether, these findings provide valuable insights that could guide future research on bovine chymosin and its biological roles.

Characterization and efficacy of Salmonella phage cocktail PHA46 in the control of Salmonella Newport and Typhimurium internalized into cherry tomatoes

Non-typhoid Salmonella enterica causes salmonellosis illness, and this bacterium can contaminate food throughout the production chain, including those that are consumed as raw products. Salmonella enterica can adhere to and internalize into fresh produce such as cherry tomatoes. It has been reported that lytic bacteriophages (phages) can be used as a biocontrol agent in the agricultural field, being an alternative for the control of Salmonella in red meat, fish, lettuce, and cabbage. The aim of this study was to characterize the two phages present in the PHA46 cocktail to determine their morphology, genome, host range, and resistance to different temperatures and pHs values; and later evaluate their lytic activity to reduce the adherence to and internalization of Salmonella enterica serovars Newport and Typhimurium into cherry tomatoes. In addition, in this work, we also explored the effect of the PHA46 cocktail on the virulence of S. Newport-45 and S. Typhimurium SL1344, recovered from the interior of cherry tomatoes, on the lifespan of the animal model Caenorhabditis elegans. The nematode C. elegans, recently has been used to test the virulence of Salmonella and it is easy to maintain and work with in the laboratory. The results revealed that the morphology obtained by Transmission Electron Microscopy of two phages from the PHA46 cocktail correspond to a myovirus, the analyses of their genomes sequences did not report virulence or antimicrobial resistance genes. The PHA46 sample is specific for 33 different serovars from different Salmonella strains and shows stability at 7 °C and pH 6. Also, the PHA46 cocktail was effective in reducing the adherence of S. Newport-45 and S. Typhimurium SL1344 to cherry tomatoes, at an average of 0.9 log10, respectively. Regarding internalized bacteria, the reduction was at an average of 1.2 log10, of the serovars mentioned above. The lifespan experiments in C. elegans showed by itself, that the PHA46 cocktail was harmless to the nematode, and the virulence from both Salmonella strains grown in vitro is diminished in the presence of the PHA46 cocktail. In conclusion, these results showed that the PHA46 cocktail could be a good candidate to be used as a biocontrol agent against Salmonella enterica.

Pan-Genome-Wide Association Study reveals a key role of the salmochelin receptor IroN in the biofilm formation of Salmonella Typhimurium and its monophasic variant 4,[5],12:i:

Salmonella enterica subsp. enterica serovar Typhimurium variant 4,[5],12:i:- (so called S. 4,[5],12:i:-) has rapidly become one of the most prevalent serovars in humans in Europe, with clinical cases associated with foodborne from pork products. The mechanisms, genetic basis and biofilms relevance by which S. 4,[5],12:i:- maintains and spreads its presence in pigs remain unclear. In this study, we examined the genetic basis of biofilm production in 78 strains of S. 4,[5],12:i:- (n = 57) and S. Typhimurium (n = 21), from human gastroenteritis, food products and asymptomatic pigs. The former showed a lower Specific Biofilm Formation index (SBF) and distant phylogenetic clades, suggesting that the ability to form biofilms is not a crucial adaptation for the S. 4,[5],12:i:- emerging success in pigs. However, using a pan-Genome-Wide Association Study (pan-GWAS) we identified genetic determinants of biofilm formation, revealing 167 common orthologous groups and genes associated with the SBF. The analysis of annotated sequences highlighted specific genetic deletions in three chromosomal regions of S. 4,[5],12:i:- correlating with SBF values: i) the complete fimbrial operon stbABCDE widely recognized as the most critical factor involved in Salmonella adherence; ii) the hxlA, hlxB, and pgiA genes, which expression in S. Typhimurium is induced in the tonsils during swine infection, and iii) the entire iroA locus related to the characteristic deletion of the second-phase flagellar genomic region in S. 4,[5],12:i:-. Consequently, we further investigated the role of the iro-genes on biofilm by constructing S. Typhimurium deletion mutants in iroBCDE and iroN. While iroBCDE showed no significant impact, iroN clearly contributed to S. Typhimurium biofilm formation. In conclusion, the pan-GWAS approach allowed us to uncover complex interactions between genetic and phenotypic factors influencing biofilm formation in S. 4,[5],12:i:- and S. Typhimurium.

Wastewater surveillance together with metaviromic data revealed the unusual resurgence of infectious diseases after the first wave of the COVID-19 outbreak

How to address public health priorities after COVID-19 is becoming a critical task. To this end, we conducted wastewater surveillance for six leading pathogens, namely, SARS-CoV-2, norovirus, rotavirus, influenza A virus (IAV), enteroviruses and respiratory syncytial virus (RSV), in Nanchang city from January to April 2023. Metaviromic sequencing was conducted at the 1st, 4th, 7th, 9th, 12th and 14th weeks to reveal the dynamics of viral pathogens that were not covered by qPCR. Amplicon sequencing of the conserved region of norovirus GI and GII and the rotavirus and region encoding nonstructural protein of RSV was also conducted weekly. The results showed that after a rapid decrease in SARS-CoV-2 sewage concentrations occurred in January 2023, surges of norovirus, rotavirus, IAV and RSV started at the 6th, 7th, 8th and 11th weeks, respectively. The dynamics of the sewage concentrations of norovirus, rotavirus, IAV and RSV were consistent with the off-season resurgence of the above infectious diseases. Notably, peak sewage concentrations of norovirus GI, GII, rotavirus, IAV and RSV were found at the 6th, 3rd, 7th, 7th and 8th weeks, respectively. Astroviruses also resurge after the 7th week, as revealed by metaviromic data, suggesting that wastewater surveillance together with metaviromic data provides an essential early warning tool for revealing patterns of infectious disease resurgence.

Dental complications as a potential indicator of Redondovirus infection: a cross-sectional study

BACKGROUND

Redondoviridae is a newly discovered virus family linked to oral and respiratory conditions in people, while there is still debate about whether it is also coinfected with other respiratory viruses. This study aimed to determine the frequency of Redondovirus (ReDoV) in nasopharyngeal samples and to investigate any possible links to SARS-CoV-2 infections.

METHODS

A polymerase chain reaction (PCR) test was conducted on 731 nasopharyngeal samples from individuals referred to medical centers in Tehran, Iran, for SARS-CoV-2 testing to investigate the prevalence of ReDoV. An oral interview was performed to complete information on dental issues and the individuals' demographics, symptoms, and vaccination history.

RESULTS

The prevalence of ReDoV was 25.99%, and 15.26% had a coinfection with SARS-CoV-2. No notable correlation was found regarding ReDoVs and SARS-CoV-2 infections (p > 0.05). Women had a higher ReDoV positivity rate of 18.47% compared to men at 7.52% (p = 0.12), and there was no significant correlation between age groups and ReDoV presence. Nonetheless, a significant association was noted between ReDoVs and dental/gum issues (p < 0.0001, OR: 13.0326). A phylogenetic analysis showed that ReDoVs originated from various human-related clusters.

CONCLUSIONS

These results highlight the potential for detecting ReDoVs in nasopharyngeal samples of people with gum or dental issues. Additionally, conducting more ReDoV epidemiological research and proposing oral health as a possible marker for ReDoV infections is important.

Three novel Enterobacter cloacae bacteriophages for therapeutic use from Ghanaian natural waters

Infections caused by multidrug-resistant (MDR) bacteria are a growing global concern. Enterobacter cloacae complex (ECC) species are particularly adept at developing antibiotic resistance. Phage therapy is proposed as an alternative treatment for pathogens that no longer respond to antibiotics. Unfortunately, ECC phages are understudied when compared to phages of many other bacterial species. In this Ghanaian-Finnish study, we isolated two ECC strains from ready-to-eat food samples and three novel phages from natural waters against these strains. We sequenced the genomic DNA of the novel Enterobacter phages, fGh-Ecl01, fGh-Ecl02, and fGh-Ecl04, and assessed their therapeutic potential. All of the phages were found to be lytic, easy to propagate, and lacking any toxic, integrase, or antibiotic resistance genes and were thus considered suitable for therapy purposes. They all were found to be related to T4-type viruses: fGh-Ecl01 and fGh-Ecl04 to karamviruses and fGh-Ecl02 to agtreviruses. Testing of Finnish clinical ECC strains showed promising susceptibility to these novel phages. As many as 61.1% of the strains were susceptible to fGh-Ecl01 and fGh-Ecl04, and 7.4% were susceptible to fGh-Ecl02. Finally, we investigated the susceptibility of the newly isolated ECC strains to three antibiotics - meropenem, ciprofloxacin, and cefepime - in combination with the novel phages. The use of phages and antibiotics together had synergistic effects. When using an antibiotic-phage combination, even low concentrations of antibiotics fully inhibited the growth of bacteria.

Interactive Tree of Life (iTOL) v6: recent updates to the phylogenetic tree display and annotation tool

The Interactive Tree Of Life (https://itol.embl.de) is an online tool for the management, display, annotation and manipulation of phylogenetic and other trees. It is freely available and open to everyone. iTOL version 6 introduces a modernized and completely rewritten user interface, together with numerous new features. A new dataset type has been introduced (colored/labeled ranges), greatly upgrading the functionality of the previous simple colored range annotation function. Additional annotation options have been implemented for several existing dataset types. Dataset template files now support simple assignment of annotations to multiple tree nodes through substring matching, including full regular expression support. Node metadata handling has been greatly extended with novel display and exporting options, and it can now be edited interactively or bulk updated through annotation files. Tree labels can be displayed using multiple simultaneous font styles, with precise positioning, sizing and styling of each individual label part. Various bulk label editing functions have been implemented, simplifying large scale changes of all tree node labels. iTOL's automatic taxonomy assignment functions now support trees based on the Genome Taxonomy Database (GTDB), in addition to the NCBI taxonomy. The functionality of the optional user account pages has been expanded, simplifying the management, navigation and sharing of projects and trees. iTOL currently handles more than one and a half million trees from >130 000 individual user accounts.

Diversification of the Rho transcription termination factor in bacteria

Correct termination of transcription is essential for gene expression. In bacteria, factor-dependent termination relies on the Rho factor, that classically has three conserved domains. Some bacteria also have a functional insertion region. However, the variation in Rho structure among bacteria has not been analyzed in detail. This study determines the distribution, sequence conservation, and predicted features of Rho factors with diverse domain architectures by analyzing 2730 bacterial genomes. About half (49.8%) of the species analyzed have the typical Escherichia coli like Rho while most of the other species (39.8%) have diverse, atypical forms of Rho. Besides conservation of the main domains, we describe a duplicated RNA-binding domain present in specific species and novel variations in the bicyclomycin binding pocket. The additional regions observed in Rho proteins exhibit remarkable diversity. Commonly, however, they have exceptional amino acid compositions and are predicted to be intrinsically disordered, to undergo phase separation, or have prion-like behavior. Phase separation has recently been shown to play roles in Rho function and bacterial fitness during harsh conditions in one species and this study suggests a more widespread role. In conclusion, diverse atypical Rho factors are broadly distributed among bacteria, suggesting additional cellular roles.

Emergence of novel methicillin resistant Staphylococcus pseudintermedius lineages revealed by whole genome sequencing of isolates from companion animals and humans in Scotland

Staphylococcus pseudintermedius is an opportunistic pathogen in dogs, and infection in humans is increasingly found, often linked to contact with dogs. We conducted a retrospective genotyping and antimicrobial susceptibility testing study of 406 S. pseudintermedius isolates cultured from animals (dogs, cats and an otter) and humans across Scotland, from 2007 to 2020. Seventy-five sequence types (STs) were identified, among the 130 isolates genotyped, with 59 seen only once. We observed the emergence of two methicillin resistant Staphylococcus pseudintermedius (MRSP) clones in Scotland: ST726, a novel locally-evolving clone, and ST551, first reported in 2015 in Poland, possibly linked to animal importation to Scotland from Central Europe. While ST71 was the most frequent S. pseudintermedius strain detected, other lineages that have been replacing ST71 in other countries, in addition to ST551, were detected. Multidrug resistance (MDR) was detected in 96.4% of MRSP and 8.4% of MSSP. A single MRSP isolate was resistant to mupirocin. Continuous surveillance for the emergence and dissemination of novel MDR MRSP in animals and humans and changes in antimicrobial susceptibility in S. pseudintermedius is warranted to minimise the threat to animal and human health.

Unveiling the microevolution of antimicrobial resistance in selected Pseudomonas aeruginosa isolates from Egyptian healthcare settings: A genomic approach

The incidence of Pseudomonas aeruginosa infections in healthcare environments, particularly in low-and middle-income countries, is on the rise. The purpose of this study was to provide comprehensive genomic insights into thirteen P. aeruginosa isolates obtained from Egyptian healthcare settings. Phenotypic analysis of the antimicrobial resistance profile and biofilm formation were performed using minimum inhibitory concentration and microtiter plate assay, respectively. Whole genome sequencing was employed to identify sequence typing, resistome, virulome, and mobile genetic elements. Our findings indicate that 92.3% of the isolates were classified as extensively drug-resistant, with 53.85% of these demonstrating strong biofilm production capabilities. The predominant clone observed in the study was ST773, followed by ST235, both of which were associated with the O11 serotype. Core genome multi-locus sequence typing comparison of these clones with global isolates suggested their potential global expansion and adaptation. A significant portion of the isolates harbored Col plasmids and various MGEs, all of which were linked to antimicrobial resistance genes. Single nucleotide polymorphisms in different genes were associated with the development of antimicrobial resistance in these isolates. In conclusion, this pilot study underscores the prevalence of extensively drug-resistant P. aeruginosa isolates and emphasizes the role of horizontal gene transfer facilitated by a diverse array of mobile genetic elements within various clones. Furthermore, specific insertion sequences and mutations were found to be associated with antibiotic resistance.

Identification and characterization of the RcTCP gene family and its expression in response to abiotic stresses in castor bean

BACKGROUND

The TCP (teosinte branched1/cincinnata/proliferating cell factor) family plays a prominent role in plant development and stress responses. However, TCP family genes have thus far not been identified in castor bean, and therefore an understanding of the expression and functional aspects of castor bean TCP genes is lacking. To identify the potential biological functions of castor bean (RcTCP) TCP members, the composition of RcTCP family members, their basic physicochemical properties, subcellular localizations, interacting proteins, miRNA target sites, and gene expression patterns under stress were assessed.

RESULTS

The presence of 20 RcTCP genes on the nine chromosomes of castor bean was identified, all of which possess TCP domains. Phylogenetic analysis indicated a close relationship between RcTCP genes and Arabidopsis AtTCP genes, suggesting potential functional similarity. Subcellular localization experiments confirmed that RcTC01/02/03/10/16/18 are all localized in the nucleus. Protein interaction analysis revealed that the interaction quantity of RcTCP03/06/11 proteins is the highest, indicating a cascade response in the functional genes. Furthermore, it was found that the promoter region of RcTCP genes contains a large number of stress-responsive elements and hormone-induced elements, indicating a potential link between RcTCP genes and stress response functions. qRT-PCR showed that all RcTCP genes exhibit a distinct tissue-specific expression pattern and their expression is induced by abiotic stress (including low temperature, abscisic acid, drought, and high salt). Among them, RcTCP01/03/04/08/09/10/14/15/18/19 genes may be excellent stress-responsive genes.

CONCLUSION

We discovered that RcTCP genes play a crucial role in various activities, including growth and development, the stress response, and transcription. This study provides a basis for studying the function of RcTCP gene in castor.

The ACE1 secondary metabolite gene cluster is a pathogenicity factor of wheat blast fungus

Wheat blast caused by Pyricularia oryzae pathotype Triticum is now becoming a very serious threat to global food security. Here, we report an essential pathogenicity factor of the wheat blast fungus that is recognized and may be targeted by a rice resistance gene. Map-based cloning of Pwt2 showed that its functional allele is the ACE1 secondary metabolite gene cluster of the wheat blast fungus required for its efficient penetration of wheat cell walls. ACE1 is required for the strong aggressiveness of Triticum, Eleusine, and Lolium pathotypes on their respective hosts, but not for that of Oryza and Setaria pathotypes on rice and foxtail millet, respectively. All ACE1 alleles found in wheat blast population are recognized by a rice resistance gene, Pi33, when introduced into rice blast isolates. ACE1 mutations for evading the recognition by Pi33 do not affect the aggressiveness of the rice blast fungus on rice but inevitably impair the aggressiveness of the wheat blast fungus on wheat. These results suggest that a blast resistance gene already defeated in rice may be revived as a durable resistance gene in wheat by targeting an Achilles heel of the wheat blast fungus.

Viperin immunity evolved across the tree of life through serial innovations on a conserved scaffold

Evolutionary arms races between cells and viruses drive the rapid diversification of antiviral genes in diverse life forms. Recent discoveries have revealed the existence of immune genes that are shared between prokaryotes and eukaryotes and show molecular and mechanistic similarities in their response to viruses. However, the evolutionary dynamics underlying the conservation and adaptation of these antiviral genes remain mostly unexplored. Here, we show that viperins constitute a highly conserved family of immune genes across diverse prokaryotes and eukaryotes and identify mechanisms by which they diversified in eukaryotes. Our findings indicate that viperins are enriched in Asgard archaea and widely distributed in all major eukaryotic clades, suggesting their presence in the last eukaryotic common ancestor and their acquisition in eukaryotes from an archaeal lineage. We show that viperins maintain their immune function by producing antiviral nucleotide analogues and demonstrate that eukaryotic viperins diversified through serial innovations on the viperin gene, such as the emergence and selection of substrate specificity towards pyrimidine nucleotides, and through partnerships with genes maintained through genetic linkage, notably with nucleotide kinases. These findings unveil biochemical and genomic transitions underlying the adaptation of immune genes shared by prokaryotes and eukaryotes. Our study paves the way for further understanding of the conservation of immunity across domains of life.

VirRep: a hybrid language representation learning framework for identifying viruses from human gut metagenomes

Identifying viruses from metagenomes is a common step to explore the virus composition in the human gut. Here, we introduce VirRep, a hybrid language representation learning framework, for identifying viruses from human gut metagenomes. VirRep combines a context-aware encoder and an evolution-aware encoder to improve sequence representation by incorporating k-mer patterns and sequence homologies. Benchmarking on both simulated and real datasets with varying viral proportions demonstrates that VirRep outperforms state-of-the-art methods. When applied to fecal metagenomes from a colorectal cancer cohort, VirRep identifies 39 high-quality viral species associated with the disease, many of which cannot be detected by existing methods.

Genetic factors explaining anthocyanin pigmentation differences

BACKGROUND

Anthocyanins are important contributors to coloration across a wide phylogenetic range of plants. Biological functions of anthocyanins span from reproduction to protection against biotic and abiotic stressors. Owing to a clearly visible phenotype of mutants, the anthocyanin biosynthesis and its sophisticated regulation have been studied in numerous plant species. Genes encoding the anthocyanin biosynthesis enzymes are regulated by a transcription factor complex comprising MYB, bHLH and WD40 proteins.

RESULTS

A systematic comparison of anthocyanin-pigmented vs. non-pigmented varieties was performed within numerous plant species covering the taxonomic diversity of flowering plants. The literature was screened for cases in which genetic factors causing anthocyanin loss were reported. Additionally, transcriptomic data sets from four previous studies were reanalyzed to determine the genes possibly responsible for color variation based on their expression pattern. The contribution of different structural and regulatory genes to the intraspecific pigmentation differences was quantified. Differences concerning transcription factors are by far the most frequent explanation for pigmentation differences observed between two varieties of the same species. Among the transcription factors in the analyzed cases, MYB genes are significantly more prone to account for pigmentation differences compared to bHLH or WD40 genes. Among the structural genes, DFR genes are most often associated with anthocyanin loss.

CONCLUSIONS

These findings support previous assumptions about the susceptibility of transcriptional regulation to evolutionary changes and its importance for the evolution of novel coloration phenotypes. Our findings underline the particular significance of MYBs and their apparent prevalent role in the specificity of the MBW complex.

Susceptibility Testing of Environmental and Clinical Aspergillus sydowii Demonstrates Potent Activity of Various Antifungals

The genus Aspergillus consists of a vast number of medically and environmentally relevant species. Aspergillus species classified in series Versicolores are ubiquitous in the environment and include the opportunistic pathogen Aspergillus sydowii, which is associated with onychomycosis and superficial skin infections. Despite frequent clinical reports of A. sydowii and related series Versicolores species, antifungal susceptibility data are scarce, hampering optimal treatment choices and subsequent patient outcomes. Here, we employed antifungal susceptibility testing (AFST) based on microbroth dilution on a set of 155 series Versicolores strains using the common antifungals amphotericin B, itraconazole, voriconazole, posaconazole, isavuconazole and micafungin with the addition of luliconazole and olorofim. All strains were identified using partial calmodulin gene sequencing, with 145 being A. sydowii, seven A. creber and three A. versicolor, using the latest taxonomic insights. Overall, tested antifungals were potent against the entire strain collection. In comparison to A. fumigatus, azole and amphotericin B MICs were slightly elevated for some strains. AFST with luliconazole and olorofim, here reported for the first time, displayed the highest in vitro activity, making these antifungals interesting alternative drugs but clinical studies are warranted for future therapeutic use.

Prediction of anastomotic insufficiency based on the mucosal microbiome prior to colorectal surgery: a proof-of-principle study

Anastomotic leakage (AL) is a potentially life-threatening complication following colorectal cancer (CRC) resection. In this study, we aimed to unravel longitudinal changes in microbial structure before, during, and after surgery and to determine if microbial alterations may be predictive for risk assessment between sufficient anastomotic healing (AS) and AL prior surgery. We analysed the microbiota of 134 colon mucosal biopsies with 16S rRNA V1-V2 gene sequencing. Samples were collected from three location sites before, during, and after surgery, and patients received antibiotics after the initial collection and during surgery. The microbial structure showed dynamic surgery-related changes at different time points. Overall bacterial diversity and the abundance of some genera such as Faecalibacterium or Alistipes decreased over time, while the genera Enterococcus and Escherichia_Shigella increased. The distribution of taxa between AS and AL revealed significant differences in the abundance of genera such as Prevotella, Faecalibacterium and Phocaeicola. In addition to Phocaeicola, Ruminococcus2 and Blautia showed significant differences in abundance between preoperative sample types. ROC analysis of the predictive value of these genera for AL revealed an AUC of 0.802 (p = 0.0013). In summary, microbial composition was associated with postoperative outcomes, and the abundance of certain genera may be predictive of postoperative complications.

Historic dog Furs Unravel the Origin and Artificial Selection of Modern Nordic Lapphund and Elkhound dog Breeds

The origins and extreme morphological evolution of the modern dog breeds are poorly studied because the founder populations are extinct. Here, we analyse eight 100 to 200 years old dog fur samples obtained from traditional North Swedish clothing, to explore the origin and artificial selection of the modern Nordic Lapphund and Elkhound dog breeds. Population genomic analysis confirmed the Lapphund and Elkhound breeds to originate from the local dog population, and showed a distinct decrease in genetic diversity in agreement with intense breeding. We identified eleven genes under positive selection during the breed development. In particular, the MSRB3 gene, associated with breed-related ear morphology, was selected in all Lapphund and Elkhound breeds, and functional assays showed that a SNP mutation in the 3'UTR region suppresses its expression through miRNA regulation. Our findings demonstrate analysis of near-modern dog artifacts as an effective tool for interpreting the origin and artificial selection of the modern dog breeds.

Comparative genomics of canine parvovirus in South America: Diversification patterns in local populations

Canine parvovirus (CPV) is a significant pathogen in domestic dogs worldwide, causing a severe and often fatal disease. CPV comprises three antigenic variants (2a, 2b, and 2c) distributed unevenly among several phylogenetic groups. The present study compared genetic variability and evolutionary patterns in South American CPV populations. We collected samples from puppies suspected of CPV infection in the neighboring Argentina and Uruguay. Antigenic variants were preliminarily characterized using PCR-RFLP and partial vp2 sequencing. Samples collected in Argentina during 2008-2018 were mainly of the 2c variant. In the Uruguayan strains (2012-2019), the 2a variant wholly replaced the 2c from 2014. Full-length coding genome and vp2 sequences were compared with global strains. The 2c and 2a strains fell by phylogenetic analysis into two phylogroups (Europe I and Asia I). The 2c strains from Argentina and Uruguay clustered in the Europe I group, with strains from America, Europe, Asia, and Oceania. Europe I is widely distributed in South America in the dog population and is also being detected in the wildlife population. The 2a strains from Uruguay formed the distinct Asia I group with strains from Asia, Africa, America, and Oceania. This Asia I group is increasing its distribution in South America and worldwide. Our research reveals high genetic variability in adjacent synchronic samples and different evolutionary patterns in South American CPV. We also highlight the importance of ancestral migrations and local diversification in the evolution of global CPV strains.