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

Enteric Populations of Escherichia coli are Likely to be Resistant to Phages Due to O Antigen Expression

There is a surfeit of bioinformatic data showing that bacteriophages abound in the enteric microbiomes of humans. What is the contribution of these viruses in shaping the bacterial strain and species composition of the gut microbiome and how are these phages maintained over time? To address these questions, we performed experiments with Escherichia coli and phages isolated from four fecal microbiota transplantation (FMT) doses as representative samples of non-dysbiotic enteric microbiota and develop and analyze the properties of a mathematical model of the population and evolutionary dynamics of bacteria and phage. Our models predict and experiments confirm that due to production of the O antigen, E. coli in the enteric microbiome are likely to be resistant to infection with co-occurring phages. Furthermore, our modeling suggests that the phages can be maintained in the population due to the high rates of host transition between resistant and sensitive states, which we call leaky resistance. Based on our observations and model predictions, we postulate that the phages found in the human gut are likely to play little role in shaping the composition of E. coli at the strain level in the enteric microbiome in healthy individuals. How general this is for other species of bacteria in the enteric flora is not yet clear, although O antigen expression is common across many taxa.

Crape myrtle LiGAoxs displaying activities of gibberellin oxidases respond to branching architecture

In the realm of ornamental horticulture, crape myrtle (Lagerstroemia indica) stands out for its aesthetic appeal, attributed largely to its vibrant flowers and distinctive branching architecture. This study embarked on a comprehensive exploration of the gibberellin oxidase (GAox) gene family in crape myrtle, illuminating its pivotal role in regulating GA levels, a key determinant of plant developmental processes. We identified and characterized 36 LiGAox genes, subdivided into GA2ox, GA3ox, GA20ox, and GAox-like subgroups, through genomic analyses. These genes' evolutionary trajectories were delineated, revealing significant gene expansions attributed to segmental duplication events. Functional analyses highlighted the divergent expression patterns of LiGAox genes across different crape myrtle varieties, associating them with variations in flower color and branching architecture. Enzymatic activity assays on selected LiGA2ox enzymes exhibited pronounced GA2 oxidase activity, suggesting a potential regulatory role in GA biosynthesis. Our findings offered a novel insight into the molecular underpinnings of GA-mediated growth and development in L. indica, providing a foundational framework for future genetic enhancements aimed at optimizing ornamental traits.

N-glycosylation on hemagglutinin head reveals inter-branch antigenic variability of avian influenza virus H5-subtypes

H5-subtype avian influenza virus (AIV) is globally prevalent and undergoes frequent antigenic drift, necessitating regular updates to vaccines. One of the many influencing elements that cause incompatibility between vaccinations and epidemic strains is the dynamic alteration of glycosylation sites. However, the biological significance of N-glycosylation in the viral evolution and antigenic changes is unclear. Here, we performed a systematic analysis of glycosylation sites on the HA1 subunit of H5N1, providing insights into the changes of primary glycosylation sites, including 140 N, 156 N, and 170 N within the antigenic epitopes of HA1 protein. Multiple recombinant viruses were then generated based on HA genes of historical vaccine strains and deactivated for immunizing SPF chickens. Inactivated recombinant strains showed relatively closer antigenicity compared to which has identical N-glycosylation patterns. The N-glycosylation modification discrepancy highlights the inter-branch antigenic diversity of H5-subtype viruses in avian influenza and serves as a vital foundation for improving vaccination tactics.

A Multilocus Sequence Typing Scheme for Rapid Identification of Xanthomonas citri Based on Whole-Genome Sequencing Data

Xanthomonas citri is a plant-pathogenic bacterium associated with a diverse range of host plant species. It has undergone substantial reclassification and currently consists of 14 different subspecies or pathovars that are responsible for a wide range of plant diseases. Whole-genome sequencing (WGS) provides a cutting-edge advantage over other diagnostic techniques in epidemiological and evolutionary studies of X. citri because it has a higher discriminatory power and is replicable across laboratories. WGS also allows for the improvement of multilocus sequence typing (MLST) schemes. In this study, we used genome sequences of Xanthomonas isolates from the NCBI RefSeq database to develop a seven-gene MLST scheme that yielded 19 sequence types (STs) that correlated with phylogenetic clades of X. citri subspecies or pathovars. Using this MLST scheme, we examined 2,911 Xanthomonas species assemblies from NCBI GenBank and identified 15 novel STs from 37 isolates that were misclassified in NCBI. In total, we identified 545 X. citri assemblies from GenBank with 95% average nucleotide identity to the X. citri type strain, and all were classified as one of the 34 STs. All MLST classifications correlated with a phylogenetic position inferred from alignments using 92 conserved genes. We observed several instances where strains from different pathovars formed closely related monophyletic clades and shared the same ST, indicating that further investigation of the validity of these pathovars is required. Our MLST scheme described here is a robust tool for rapid classification of X. citri pathovars using WGS and a powerful method for further comprehensive taxonomic revision of X. citri pathovars.

Diversity and structure of the deep-sea sponge microbiome in the equatorial Atlantic Ocean

Sponges (phylum Porifera) harbour specific microbial communities that drive the ecology and evolution of the host. Understanding the structure and dynamics of these communities is emerging as a primary focus in marine microbial ecology research. Much of the work to date has focused on sponges from warm and shallow coastal waters, while sponges from the deep ocean remain less well studied. Here, we present a metataxonomic analysis of the microbial consortia associated with 23 individual deep-sea sponges. We identify a high abundance of archaea relative to bacteria across these communities, with certain sponge microbiomes comprising more than 90 % archaea. Specifically, the archaeal family Nitrosopumilaceae is prolific, comprising over 99 % of all archaeal reads. Our analysis revealed that sponge microbial communities reflect the host sponge phylogeny, indicating a key role for host taxonomy in defining microbiome composition. Our work confirms the contribution of both evolutionary and environmental processes to the composition of microbial communities in deep-sea sponges.

Delineating the evolutionary pathway to multidrug-resistant outbreaks of a Mycobacterium tuberculosis L4.1.2.1/Haarlem sublineage

OBJECTIVES

We sought to capture the evolutionary itinerary of the Mycobacterium tuberculosis L4.1.2.1/Haarlem sublineage in northern Tunisia, where it caused a major multidrug-resistant (MDR) tuberculosis outbreak in a context strictly negative for HIV infection.

METHODS

We combined whole genome sequencing and Bayesian approaches using a representative collection of drug-susceptible and drug-resistant L4.1.2.1/Haarlem clinical strains (n = 121) recovered from the outbreak region over 16 years.

RESULTS

In the absence of drug resistance, the L4.1.2.1/Haarlem sublineage showed a propensity for rapid transmission as witnessed by the high clustering (44.6%) and recent transmission rates (25%), as well as the reduced mean distance between genome pairs. The entire pool of L4.1.2.1/Haarlem MDR strains was found to be linked to either the aforementioned major outbreak (68 individuals, 2001-2016) or to a minor, newly uncovered outbreak (six cases, 2001-2011). Strikingly, the two outbreaks descended independently from a common ancestor that can be dated back to 1886.

CONCLUSIONS

Our data point to the intrinsic propensity for rapid transmission of the M. tuberculosis L4.1.2.1/Haarlem sublineage in northern Tunisia, linking the overall MDR tuberculosis epidemic to a single ancestor. These findings bring out the important role of the bacillus' genetic background in the emergence of successful MDR M. tuberculosis clones.

Associations among MHC genes, latitude, and avian malaria infections in the rufous-collared sparrow (Zonotrichia capensis)

The major histocompatibility complex (MHC) is a genetic region in jawed vertebrates that contains key genes involved in the immune response. Associations between the MHC and avian malaria infections in wild birds have been observed and mainly explored in the Northern Hemisphere, while a general lack of information remains in the Southern Hemisphere. Here, we investigated the associations between the MHC genes and infections with Plasmodium and Haemoproteus blood parasites along a latitudinal gradient in South America. We sampled 93 rufous-collared sparrows (Zonotrichia capensis) individuals from four countries, Colombia, Ecuador, Peru, and Chile, and estimated MHC-I and MHC-II allele diversity. We detected between 1-4 (MHC-I) and 1-6 (MHC-II) amino acidic alleles per individual, with signs of positive selection. We obtained generalized additive mixed models to explore the associations between MHC-I and MHC-II diversity and latitude. We also explored the relationship between infection status and latitude/biome. We found a non-linear association between the MHC-II amino acidic allele diversity and latitude. Individuals from north Chile presented a lower MHC genetic diversity than those from other locations. We also found an association between deserts and xeric shrublands and a lower prevalence of Haemoproteus parasites. Our results support a lower MHC genetic in arid or semi-arid habitats in the region with the lower prevalence of Haemoproteus parasites.

Emergence of novel hypervirulent Acinetobacter baumannii strain and herpes simplex type 1 virus in a case of community-acquired pneumonia in China

BACKGROUND

A. baumannii is an important and common clinical pathogen, especially in the intensive care unit (ICU). This study aimed to characterize one hypervirulent A. baumannii strain in a patient with community-acquired pneumonia and herpes simplex type 1 virus infection.

METHODS

Minimum inhibitory concentrations (MICs) were determined using the Kirby-Bauer (K-B) and broth microdilution methods. Galleria mellonella infection model experiment was conducted. Whole-genome sequencing (WGS) was performed using the Illumina and Nanopore platforms. The resistance and virulence determinants were identified using the ABRicate program with ResFinder and the VFDB database. The capsular polysaccharide locus (K locus) and lipooligosaccharide outer core locus (OC locus) were identified using Kleborate with Kaptive. Phylogenetic analyses were conducted using the BacWGSTdb server.

RESULTS

A. baumannii XH2146 strain belongs to ST10Pas and ST447Oxf. The strain was resistant to cefazolin, ciprofloxacin, and trimethoprim/sulfamethoxazole (TMP-SMX). Bautype and Kaptive analyses showed that XH2146 contains OCL2 and KL49. WGS analysis revealed that the strain harbored blaADC-76, blaOXA-68, ant(3'')-IIa, tet(B), and sul2. Notably, tet(B) and sul2, both were located within a 114,700-bp plasmid (designated pXH2146-1). Virulence assay revealed A. baumannii XH2146 possessed higher virulence than A. baumannii AB5075 at 12 h. Comparative genomic analysis showed that A. baumannii ST447 strains were mainly isolated from the USA and exhibited a relatively close genetic relationship. Importantly, 11 strains were observed to carry blaOXA-58; blaOXA-23 was identified in 11 isolates and three ST447 A. baumannii strains harbored blaNDM-1.

CONCLUSIONS

Early detection of community-acquired hypervirulent Acinetobacter baumannii strains is recommended to prevent their extensive spread in hospitals.

Genomic and functional insights of a mucin foraging Rhodopirellula halodulae sp. nov

Nine novel strains were obtained from various algal and seagrass samples. The analysis of the 16S rRNA gene-based phylogenetic tree revealed monophyletic placement of all novel strains within the Rhodopirellula genus. The type strain was identified as JC737T, which shared 99.1 % 16S rRNA gene sequence identity with Rhodopirellula baltica SH1T, while strain JC740 was designated as an additional strain. The genome sizes of strains JC737T and JC740 were 6.6 and 6.7 Mb, respectively, and the G + C content was 56.2 %. The strains cladded distinctly in the phylogenomic tree, and the ANI and dDDH values of the strain JC737T were 75.8-76.1 % and 20.8-21.3 %, respectively, in comparison to other Rhodopirellula members. The strain demonstrated a versatile degradation capability, exhibiting a diverse array of complex polysaccharides, including mucin which had not been previously identified within the members of the phylum Planctomycetota. The phylogenomic, pan-genomic, morphological, physiological, and genomic characterization of the strain lead to the proposal to describe the strain as Rhodopirellula halodulae sp. nov.

Zinc mediates control of nitrogen fixation via transcription factor filamentation

Plants adapt to fluctuating environmental conditions by adjusting their metabolism and gene expression to maintain fitness1. In legumes, nitrogen homeostasis is maintained by balancing nitrogen acquired from soil resources with nitrogen fixation by symbiotic bacteria in root nodules2-8. Here we show that zinc, an essential plant micronutrient, acts as an intracellular second messenger that connects environmental changes to transcription factor control of metabolic activity in root nodules. We identify a transcriptional regulator, FIXATION UNDER NITRATE (FUN), which acts as a sensor, with zinc controlling the transition between an inactive filamentous megastructure and an active transcriptional regulator. Lower zinc concentrations in the nodule, which we show occur in response to higher levels of soil nitrate, dissociates the filament and activates FUN. FUN then directly targets multiple pathways to initiate breakdown of the nodule. The zinc-dependent filamentation mechanism thus establishes a concentration readout to adapt nodule function to the environmental nitrogen conditions. In a wider perspective, these results have implications for understanding the roles of metal ions in integration of environmental signals with plant development and optimizing delivery of fixed nitrogen in legume crops.

Varroa destructor shapes the unique viral landscapes of the honey bee populations of the Azores archipelago

The worldwide dispersal of the ectoparasitic mite Varroa destructor from its Asian origins has fundamentally transformed the relationship of the honey bee (Apis mellifera) with several of its viruses, via changes in transmission and/or host immunosuppression. The extent to which honey bee-virus relationships change after Varroa invasion is poorly understood for most viruses, in part because there are few places in the world with several geographically close but completely isolated honey bee populations that either have, or have not, been exposed long-term to Varroa, allowing for separate ecological, epidemiological, and adaptive relationships to develop between honey bees and their viruses, in relation to the mite's presence or absence. The Azores is one such place, as it contains islands with and without the mite. Here, we combined qPCR with meta-amplicon deep sequencing to uncover the relationship between Varroa presence, and the prevalence, load, diversity, and phylogeographic structure of eight honey bee viruses screened across the archipelago. Four viruses were not detected on any island (ABPV-Acute bee paralysis virus, KBV-Kashmir bee virus, IAPV-Israeli acute bee paralysis virus, BeeMLV-Bee macula-like virus); one (SBV-Sacbrood virus) was detected only on mite-infested islands; one (CBPV-Chronic bee paralysis virus) occurred on some islands, and two (BQCV-Black queen cell virus, LSV-Lake Sinai virus,) were present on every single island. This multi-virus screening builds upon a parallel survey of Deformed wing virus (DWV) strains that uncovered a remarkably heterogeneous viral landscape featuring Varroa-infested islands dominated by DWV-A and -B, Varroa-free islands naïve to DWV, and a refuge of the rare DWV-C dominating the easternmost Varroa-free islands. While all four detected viruses investigated here were affected by Varroa for one or two parameters (usually prevalence and/or the Richness component of ASV diversity), the strongest effect was observed for the multi-strain LSV. Varroa unambiguously led to elevated prevalence, load, and diversity (Richness and Shannon Index) of LSV, with these results largely shaped by LSV-2, a major LSV strain. Unprecedented insights into the mite-virus relationship were further gained from implementing a phylogeographic approach. In addition to enabling the identification of a novel LSV strain that dominated the unique viral landscape of the easternmost islands, this approach, in combination with the recovered diversity patterns, strongly suggests that Varroa is driving the evolutionary change of LSV in the Azores. This study greatly advances the current understanding of the effect of Varroa on the epidemiology and adaptive evolution of these less-studied viruses, whose relationship with Varroa has thus far been poorly defined.

TnpB homologues exapted from transposons are RNA-guided transcription factors

Transposon-encoded tnpB and iscB genes encode RNA-guided DNA nucleases that promote their own selfish spread through targeted DNA cleavage and homologous recombination1-4. These widespread gene families were repeatedly domesticated over evolutionary timescales, leading to the emergence of diverse CRISPR-associated nucleases including Cas9 and Cas12 (refs. 5,6). We set out to test the hypothesis that TnpB nucleases may have also been repurposed for novel, unexpected functions other than CRISPR-Cas adaptive immunity. Here, using phylogenetics, structural predictions, comparative genomics and functional assays, we uncover multiple independent genesis events of programmable transcription factors, which we name TnpB-like nuclease-dead repressors (TldRs). These proteins use naturally occurring guide RNAs to specifically target conserved promoter regions of the genome, leading to potent gene repression in a mechanism akin to CRISPR interference technologies invented by humans7. Focusing on a TldR clade found broadly in Enterobacteriaceae, we discover that bacteriophages exploit the combined action of TldR and an adjacently encoded phage gene to alter the expression and composition of the host flagellar assembly, a transformation with the potential to impact motility8, phage susceptibility9, and host immunity10. Collectively, this work showcases the diverse molecular innovations that were enabled through repeated exaptation of transposon-encoded genes, and reveals the evolutionary trajectory of diverse RNA-guided transcription factors.

Phylogenomic analyses of the genus Actinoplanes: description of four novel genera

The genus Actinoplanes comprises 57 species (in February 2024) that are important components of ecosystems and are widely used in biotechnology, especially pharmaceuticals. Phylogenetic analysis of the family Micromonosporaceae (based on the 16S rRNA gene sequence) allowed us to group members of different genera into separate clades; however, the genus Actinoplanes was divided into three separate clades. Such phylogenetic heterogeneity could be due to the limitations of 16S rRNA gene analysis. In response to this heterogeneity, genomic phylogeny was performed. Phylogenomic reconstruction based on 324 single-copy orthologous genes allowed us to divide the genus Actinoplanes first into four clades and then, based on average nucleotide identity analysis, into five clades. Finally, chemotaxonomic analysis of each clade confirmed each clade's distinctiveness and the necessity to reclassify the genus Actinoplanes. The obtained data allowed us to divide the genus Actinoplanes into five genera: Actinoplanes, Paractinoplanes, Winogradskya, Symbioplanes and Amorphoplanes.

Five centuries of genome evolution and multi-host adaptation of Campylobacter jejuni in Brazil

Consumption of raw, undercooked or contaminated animal food products is a frequent cause of Campylobacter jejuni infection. Brazil is the world's third largest producer and a major exporter of chicken meat, yet population-level genomic investigations of C. jejuni in the country remain scarce. Analysis of 221 C. jejuni genomes from Brazil shows that the overall core and accessory genomic features of C. jejuni are influenced by the identity of the human or animal source. Of the 60 sequence types detected, ST353 is the most prevalent and consists of samples from chicken and human sources. Notably, we identified the presence of diverse bla genes from the OXA-61 and OXA-184 families that confer beta-lactam resistance as well as the operon cmeABCR related to multidrug efflux pump, which contributes to resistance against tetracyclines, macrolides and quinolones. Based on limited data, we estimated the most recent common ancestor of ST353 to the late 1500s, coinciding with the time the Portuguese first arrived in Brazil and introduced domesticated chickens into the country. We identified at least two instances of ancestral chicken-to-human infections in ST353. The evolution of C. jejuni in Brazil was driven by the confluence of clinically relevant genetic elements, multi-host adaptation and clonal population growth that coincided with major socio-economic changes in poultry farming.

Development and validation of YARN: A novel SE-400 MPS kit for East Asian paternal lineage analysis

Y-chromosomal short tandem repeat polymorphisms (Y-STRs) and Y-chromosomal single nucleotide polymorphisms (Y-SNPs) are valuable genetic markers used in paternal lineage identification and population genetics. Currently, there is a lack of an effective panel that integrates Y-STRs and Y-SNPs for studying paternal lineages, particularly in East Asian populations. Hence, we developed a novel Y-chromosomal targeted panel called YARN (Y-chromosome Ancestry and Region Network) based on multiplex PCR and a single-end 400 massive parallel sequencing (MPS) strategy, consisting of 44 patrilineage Y-STRs and 260 evolutionary Y-SNPs. A total of 386 reactions were validated for the effectiveness and applicability of YARN according to SWGDAM validation guidelines, including sensitivity (with a minimum input gDNA of 0.125 ng), mixture identification (ranging from 1:1-1:10), PCR inhibitor testing (using substances such as 50 μM hematin, 100 μM hemoglobin, 100 μM humic acid, and 2.5 mM indigo dye), species specificity (successfully distinguishing humans from other animals), repeatability study (achieved 100% accuracy), and concordance study (with 99.91% accuracy for 1121 Y-STR alleles). Furthermore, we conducted a pilot study using YARN in a cohort of 484 Han Chinese males from Huaiji County, Zhaoqing City, Guangdong, China (GDZQHJ cohort). In this cohort, we identified 52 different Y-haplogroups and 73 different surnames. We found weak to moderate correlations between the Y-haplogroups, Chinese surnames, and geographical locations of the GDZQHJ cohort (with λ values ranging from 0.050 to 0.340). However, when we combined two different categories into a new independent variable, we observed stronger correlations (with λ values ranging from 0.617 to 0.754). Overall, the YARN panel, which combines Y-STR and Y-SNP genetic markers, meets forensic DNA quality assurance guidelines and holds potential for East Asian geographical origin inference and paternal lineage analysis.

Single-cell transcriptomics across 2,534 microbial species reveals functional heterogeneity in the rumen microbiome

Deciphering the activity of individual microbes within complex communities and environments remains a challenge. Here we describe the development of microbiome single-cell transcriptomics using droplet-based single-cell RNA sequencing and pangenome-based computational analysis to characterize the functional heterogeneity of the rumen microbiome. We generated a microbial genome database (the Bovine Gastro Microbial Genome Map) as a functional reference map for the construction of a single-cell transcriptomic atlas of the rumen microbiome. The atlas includes 174,531 microbial cells and 2,534 species, of which 172 are core active species grouped into 12 functional clusters. We detected single-cell-level functional roles, including a key role for Basfia succiniciproducens in the carbohydrate metabolic niche of the rumen microbiome. Furthermore, we explored functional heterogeneity and reveal metabolic niche trajectories driven by biofilm formation pathway genes within B. succiniciproducens. Our results provide a resource for studying the rumen microbiome and illustrate the diverse functions of individual microbial cells that drive their ecological niche stability or adaptation within the ecosystem.

Ancient Plasmodium genomes shed light on the history of human malaria

Malaria-causing protozoa of the genus Plasmodium have exerted one of the strongest selective pressures on the human genome, and resistance alleles provide biomolecular footprints that outline the historical reach of these species1. Nevertheless, debate persists over when and how malaria parasites emerged as human pathogens and spread around the globe1,2. To address these questions, we generated high-coverage ancient mitochondrial and nuclear genome-wide data from P. falciparum, P. vivax and P. malariae from 16 countries spanning around 5,500 years of human history. We identified P. vivax and P. falciparum across geographically disparate regions of Eurasia from as early as the fourth and first millennia BCE, respectively; for P. vivax, this evidence pre-dates textual references by several millennia3. Genomic analysis supports distinct disease histories for P. falciparum and P. vivax in the Americas: similarities between now-eliminated European and peri-contact South American strains indicate that European colonizers were the source of American P. vivax, whereas the trans-Atlantic slave trade probably introduced P. falciparum into the Americas. Our data underscore the role of cross-cultural contacts in the dissemination of malaria, laying the biomolecular foundation for future palaeo-epidemiological research into the impact of Plasmodium parasites on human history. Finally, our unexpected discovery of P. falciparum in the high-altitude Himalayas provides a rare case study in which individual mobility can be inferred from infection status, adding to our knowledge of cross-cultural connectivity in the region nearly three millennia ago.

All-RNA-mediated targeted gene integration in mammalian cells with rationally engineered R2 retrotransposons

All-RNA-mediated targeted gene integration methods, rendering reduced immunogenicity, effective deliverability with non-viral vehicles, and a low risk of random mutagenesis, are urgently needed for next-generation gene addition technologies. Naturally occurring R2 retrotransposons hold promise in this context due to their site-specific integration profile. Here, we systematically analyzed the biodiversity of R2 elements and screened several R2 orthologs capable of full-length gene insertion in mammalian cells. Robust R2 system gene integration efficiency was attained using combined donor RNA and protein engineering. Importantly, the all-RNA-delivered engineered R2 system showed effective integration activity, with efficiency over 60% in mouse embryos. Unbiased high-throughput sequencing demonstrated that the engineered R2 system exhibited high on-target integration specificity (99%). In conclusion, our study provides engineered R2 tools for applications based on hit-and-run targeted DNA integration and insights for further optimization of retrotransposon systems.

Phage anti-CRISPR control by an RNA- and DNA-binding helix-turn-helix protein

In all organisms, regulation of gene expression must be adjusted to meet cellular requirements and frequently involves helix-turn-helix (HTH) domain proteins1. For instance, in the arms race between bacteria and bacteriophages, rapid expression of phage anti-CRISPR (acr) genes upon infection enables evasion from CRISPR-Cas defence; transcription is then repressed by an HTH-domain-containing anti-CRISPR-associated (Aca) protein, probably to reduce fitness costs from excessive expression2-5. However, how a single HTH regulator adjusts anti-CRISPR production to cope with increasing phage genome copies and accumulating acr mRNA is unknown. Here we show that the HTH domain of the regulator Aca2, in addition to repressing Acr synthesis transcriptionally through DNA binding, inhibits translation of mRNAs by binding conserved RNA stem-loops and blocking ribosome access. The cryo-electron microscopy structure of the approximately 40 kDa Aca2-RNA complex demonstrates how the versatile HTH domain specifically discriminates RNA from DNA binding sites. These combined regulatory modes are widespread in the Aca2 family and facilitate CRISPR-Cas inhibition in the face of rapid phage DNA replication without toxic acr overexpression. Given the ubiquity of HTH-domain-containing proteins, it is anticipated that many more of them elicit regulatory control by dual DNA and RNA binding.

Meta-analysis reveals the effects of microbial inoculants on the biomass and diversity of soil microbial communities

Microbial inoculation involves transplanting microorganisms from their natural habitat to new plants or soils to improve plant performance, and it is being increasingly used in agriculture and ecological restoration. However, microbial inoculants can invade and alter the composition of native microbial communities; thus, a comprehensive analysis is urgently needed to understand the overall impact of microbial inoculants on the biomass, diversity, structure and network complexity of native communities. Here we provide a meta-analysis of 335 studies revealing a positive effect of microbial inoculants on soil microbial biomass. This positive effect was weakened by environmental stress and enhanced by the use of fertilizers and native inoculants. Although microbial inoculants did not alter microbial diversity, they induced major changes in the structure and bacterial composition of soil microbial communities, reducing the complexity of bacterial networks and increasing network stability. Finally, higher initial levels of soil nutrients amplified the positive impact of microbial inoculants on fungal biomass, actinobacterial biomass, microbial biomass carbon and microbial biomass nitrogen. Together, our results highlight the positive effects of microbial inoculants on soil microbial biomass, emphasizing the benefits of native inoculants and the important regulatory roles of soil nutrient levels and environmental stress.

Phylogenomic analyses of multidrug resistant Corynebacterium striatum strains isolated from patients in a tertiary care hospital in the UK

Corynebacterium striatum is an emerging nosocomial pathogen. This is the first report showing the presence of three distinct multidrug resistant lineages of C. striatum among patients in a UK hospital. The presence of ErmX, Tet(W), Bla and AmpC proteins, and mutations in gyrA gene are associated with the resistance to clindamycin, doxycycline, penicillin and moxifloxacin, respectively. These strains are equipped with several corynebacterial virulence genes including two SpaDEF-type and a novel pilus gene cluster, which needs further molecular characterisation. This study highlights a need of developing an active surveillance strategy for routine monitoring and preventing potential cross-transmission among susceptible patients.

Phylogenomic and genomic analysis reveals unique and shared genetic signatures of Mycobacterium kansasii complex species

Species belonging to the Mycobacterium kansasii complex (MKC) are frequently isolated from humans and the environment and can cause serious diseases. The most common MKC infections are caused by the species M. kansasii (sensu stricto), leading to tuberculosis-like disease. However, a broad spectrum of virulence, antimicrobial resistance and pathogenicity of these non-tuberculous mycobacteria (NTM) are observed across the MKC. Many genomic aspects of the MKC that relate to these broad phenotypes are not well elucidated. Here, we performed genomic analyses from a collection of 665 MKC strains, isolated from environmental, animal and human sources. We inferred the MKC pangenome, mobilome, resistome, virulome and defence systems and show that the MKC species harbours unique and shared genomic signatures. High frequency of presence of prophages and different types of defence systems were observed. We found that the M. kansasii species splits into four lineages, of which three are lowly represented and mainly in Brazil, while one lineage is dominant and globally spread. Moreover, we show that four sub-lineages of this most distributed M. kansasii lineage emerged during the twentieth century. Further analysis of the M. kansasii genomes revealed almost 300 regions of difference contributing to genomic diversity, as well as fixed mutations that may explain the M. kansasii's increased virulence and drug resistance.

Lacrimispora brassicae sp. nov. isolated from fermented cabbage, and proposal of Clostridium indicum Gundawar et al. 2019 and Clostridium methoxybenzovorans Mechichi et al. 1999 as heterotypic synonyms of Lacrimispora amygdalina (Parshina et al. 2003) Haas and Blanchard 2020 and Lacrimispora indolis (McClung and McCoy 1957) Haas and Blanchard 2020, respectively

A Gram-stain-negative, endospore-forming, rod-shaped, indole-producing bacterial strain, designated YZC6T, was isolated from fermented cabbage. Strain YZC6T grew at 10-37  °C, pH 5.5-8.5, and with up to 2  % (w/v) NaCl. The major cellular fatty acids were C16 : 0 and C18 : 1 cis 11 dimethyl acetal. Phylogenetic analysis of the 16S rRNA gene revealed that strain YZC6T belonged to the genus Lacrimispora and was closely related to Lacrimispora aerotolerans DSM 5434T (98.3  % sequence similarity), Lacrimispora saccharolytica WM1T (98.1  %), and Lacrimispora algidixylanolytica SPL73T (98.1  %). The average nucleotide identity based on blast (below 87.8  %) and digital DNA-DNA hybridization (below 36.1 %) values between the novel isolate and its corresponding relatives showed that strain YZC6T could be readily distinguished from its closely related species. Based on genotypic, phenotypic, and chemotaxonomic data, a novel Lacrimispora species, Lacrimispora brassicae sp. nov., was proposed, with YZC6T as the type strain (=MAFF 212518T=JCM 32810T=DSM 112100T). This study also proposed Clostridium indicum Gundawar et al. 2019 as a later heterotypic synonym of Lacrimispora amygdalina (Parshina et al. 2003) Haas and Blanchard 2020 and Clostridium methoxybenzovorans Mechichi et al. 1999 as a later heterotypic synonym of Lacrimispora indolis (McClung and McCpy 1957) Haas and Blanchard 2020.

Exotic medicinal plants from the Silk Road promote the diversification of traditional Chinese medicines

BACKGROUND

As a political, economic, and cultural exchange channel between ancient China and countries in Asia, Europe, and Africa, the Silk Road has promoted political, trade, and cultural exchanges between China and foreign countries in Chinese history and also promoted the development of traditional Chinese medicine.

METHODS

This article summarizes the introduction of medicinal materials from the Han to Qing Dynasties, spanning approximately 2000 years.

RESULTS

A total of 235 types of medicinal plant materials were imported. An analysis of 178 medicinal herbs of known origin, belonging to 72 families revealed their effectiveness in treating 20 diseases. The maximum number of medicinal herbs used to treat gastrointestinal and digestive disorders (GAS) was 122. The applications and origin of exotic medicinal materials, including draconis sanguis and olibanum have changed during the development of the Silk Road. Imported medicinal materials are affected by five factors, including local demand, adaptability, cultural exchange, scarcity, and medical theory. Five modes for introducing medicinal materials include the onshore Silk Road, the maritime Silk Road, diplomatic envoys and gifts, overseas Chinese, cultural exchange, and medical integration. The application of exotic medicinal materials expands the resources and application fields of traditional Chinese medicine, enriching the theory of traditional Chinese medicine.

CONCLUSION

Traditional Chinese medicinal compounds introduced to China through the ancient Silk Road not only promoted their integration into foreign medicine but also had long-lasting impacts to date and over a wide range, thereby considerably affecting the pharmaceutical and general healthcare industries.

Pseudomonas aeruginosa breaches respiratory epithelia through goblet cell invasion in a microtissue model

Pseudomonas aeruginosa, a leading cause of severe hospital-acquired pneumonia, causes infections with up to 50% mortality rates in mechanically ventilated patients. Despite some knowledge of virulence factors involved, it remains unclear how P. aeruginosa disseminates on mucosal surfaces and invades the tissue barrier. Using infection of human respiratory epithelium organoids, here we observed that P. aeruginosa colonization of apical surfaces is promoted by cyclic di-GMP-dependent asymmetric division. Infection with mutant strains revealed that Type 6 Secretion System activities promote preferential invasion of goblet cells. Type 3 Secretion System activity by intracellular bacteria induced goblet cell death and expulsion, leading to epithelial rupture which increased bacterial translocation and dissemination to the basolateral epithelium. These findings show that under physiological conditions, P. aeruginosa uses coordinated activity of a specific combination of virulence factors and behaviours to invade goblet cells and breach the epithelial barrier from within, revealing mechanistic insight into lung infection dynamics.

Microbial Ecology of an Arctic Travertine Geothermal Spring: Implications for Biosignature Preservation and Astrobiology

Jotun springs in Svalbard, Norway, is a rare warm environment in the Arctic that actively forms travertine. In this study, we assessed the microbial ecology of Jotun's active (aquatic) spring and dry spring transects. We evaluated the microbial preservation potential and mode, as well as the astrobiological relevance of the travertines to marginal carbonates mapped at Jezero Crater on Mars (the Mars 2020 landing site). Our results revealed that microbial communities exhibited spatial dynamics controlled by temperature, fluid availability, and geochemistry. Amorphous carbonates and silica precipitated within biofilm and on the surface of filamentous microorganisms. The water discharged at the source is warm, with near neutral pH, and undersaturated in silica. Hence, silicification possibly occurred through cooling, dehydration, and partially by a microbial presence or activities that promote silica precipitation. CO2 degassing and possible microbial contributions induced calcite precipitation and travertine formation. Jotun revealed that warm systems that are not very productive in carbonate formation may still produce significant carbonate buildups and provide settings favorable for fossilization through silicification and calcification. Our findings suggest that the potential for amorphous silica precipitation may be essential for Jezero Crater's marginal carbonates because it significantly increases the preservation potential of putative martian organisms.

Characterization of antimicrobial resistance profiles in Escherichia coli isolated from captive mammals in Ecuador

BACKGROUND

This study focuses on the AMR profiles in E. coli isolated from captive mammals at EcoZoo San Martín, Baños de Agua Santa, Ecuador, highlighting the role of wildlife as reservoirs of resistant bacteria.

AIMS

The aim of this research is to investigate the antimicrobial resistance profiles of E. coli strains isolated from various species of captive mammals, emphasizing the potential zoonotic risks and the necessity for integrated AMR management strategies.

MATERIALS & METHODS

A total of 189 fecal samples were collected from 70 mammals across 27 species. These samples were screened for E. coli, resulting in 90 identified strains. The resistance profiles of these strains to 16 antibiotics, including 10 β-lactams and 6 non-β-lactams, were determined using the disk diffusion method. Additionally, the presence of Extended-Spectrum Beta-Lactamase (ESBL) genes and other resistance genes was analyzed using PCR.

RESULTS

Significant resistance was observed, with 52.22% of isolates resistant to ampicillin, 42.22% to ceftriaxone and cefuroxime, and 27.78% identified as ESBL-producing E. coli. Multiresistance (resistance to more than three antibiotic groups) was found in 35.56% of isolates. Carnivorous and omnivorous animals, particularly those with prior antibiotic treatments, were more likely to harbor resistant strains.

DISCUSSION

These findings underscore the role of captive mammals as indicators of environmental AMR. The high prevalence of resistant E. coli in these animals suggests that zoos could be significant reservoirs for the spread of antibiotic-resistant bacteria. The results align with other studies showing that diet and antibiotic treatment history influence resistance profiles.

CONCLUSION

The study highlights the need for an integrated approach involving veterinary care, habitat management, and public awareness to prevent captive wildlife from becoming reservoirs of antibiotic-resistant bacteria. Improved waste management practices and responsible antibiotic use are crucial to mitigate the risks of AMR in zoo environments and reduce zoonotic threats.

Genomic analysis of canine pneumoviruses and canine respiratory coronavirus from New Zealand

AIMS

To isolate canine respiratory coronavirus (CRCoV) and canine pneumovirus (CnPnV) in cell culture and to compare partial genomic sequences of CRCoV and CnPnV from New Zealand with those from other countries.

METHODS

Oropharyngeal swab samples from dogs affected by canine infectious respiratory disease syndrome that were positive for CnPnV (n = 15) or CRCoV (n = 1) by virus-specific reverse transcriptase quantitative PCR (RT-qPCR) in a previous study comprised the starting material. Virus isolation was performed in HRT-18 cells for CRCoV and RAW 264.7 and Vero cells for CnPnV. The entire sequence of CnPnV G protein (1,266 nucleotides) and most (8,063/9,707 nucleotides) of the 3' region of CRCoV that codes for 10 structural and accessory proteins were amplified and sequenced. The sequences were analysed and compared with other sequences available in GenBank using standard molecular tools including phylogenetic analysis.

RESULTS

Virus isolation was unsuccessful for both CRCoV and CnPnV. Pneumovirus G protein was amplified from 3/15 (20%) samples that were positive for CnPnV RNA by RT-qPCR. Two of these (NZ-048 and NZ-049) were 100% identical to each other, and 90.9% identical to the third one (NZ-007). Based on phylogenetic analysis of the G protein gene, CnPnV NZ-048 and NZ-049 clustered with sequences from the USA, Thailand and Italy in group A, and CnPnV NZ-007 clustered with sequences from the USA in group B. The characteristics of the predicted genes (length, position) and their putative protein products (size, predicted structure, presence of N- and O-glycosylation sites) of the New Zealand CRCoV sequence were consistent with those reported previously, except for the region located between open reading frame (ORF)3 (coding for S protein) and ORF6 (coding for E protein). The New Zealand virus was predicted to encode 5.9 kDa, 27 kDa and 12.7 kDa proteins, which differed from the putative coding capacity of this region reported for CRCoV from other countries.

CONCLUSIONS

This report represents the first characterisation of partial genomic sequences of CRCoV and CnPnV from New Zealand. Our results suggest that the population of CnPnV circulating in New Zealand is not homogeneous, and that the viruses from two clades described overseas are also present here. Limited conclusions can be made based on only one CRCoV sequence, but the putative differences in the coding capacity of New Zealand CRCoV support the previously reported variability of this region. The reasons for such variability and its biological implications need to be further elucidated.

Propionate outperforms conventional acetate as electron donors for highly-sensitive electrochemical active biofilm sensors in water biotoxicity early-warning

With the ability to generate in situ real-time electric signals, electrochemically active biofilm (EAB) sensors have attracted wide attention as a promising water biotoxicity early-warning device. Organic matters serving as the electron donors potentially affect the electric signal's output and the sensitivity of the EAB sensor. To explore the influence of organic matters on EAB sensor's performance, this study tested six different organic matters during the sensor's inoculation. Besides the acetate, a conventional and widely used organic matter, propionate and lactate were also found capable of starting up the sensor. Moreover, the propionate-fed (PF) sensor delivered the highest sensitivity, which are respectively 1.4 times and 2.8 times of acetate-fed (AF) sensor and lactate-fed (LF) sensor. Further analysis revealed that EAB of PF sensor had more vulnerable intracellular metabolism than the others, which manifested as the most severe energy metabolic suppression and reactive oxygen species attack. Regarding the microbial function, a two-component system that was deemed as an environment awareness system was found in the EAB of PF, which also contributed to its high sensitivity. Finally, PF sensor was tested in real water environment to deliver early-warning signals.

Molecular Characterization and Potential Host-switching of Swine Farm associated Clostridioides difficile ST11

OBJECTIVE

To conduct molecular prevalence and genetic polymorphism analysis of 24 Swine Farm associated C. difficile ST11 strains, in addition to other representative sequenced ST strains.

METHODS

The collected C. difficile strains underwent whole genome sequencing and bioinformatic analysis using the illumina NovaSeq platform, SPAdes, Prokka, MOB-suite, and FastTree. Virulence and antibiotic resistance genes were identified through NCBI Pathogen Database. Cytotoxicity tests were conducted on HT-29 cells and Vero cells to verify the function of toxin A and toxin B.

RESULTS

The most prevalent resistance genes in ST11 were found to be against β-lactamases, aminoglycosides, and tetracycline. A C. difficile isolate (strain 27) with tcdA deletion and high antibiotic resistance genes was far apart from other swine farm associated ST11 isolates in the phylogenetic branch. The remarkable genetic similarity between animal and human C. difficile strains suggests potential transmission of ST11 strains between animals and humans. The plasmid replicon sequences repUS43 were identified in all ST11 strains except one variant (strain 27), and 91.67% (22/24) of these were assessed by MOB-typer as having mobilizable plasmids.

CONCLUSION

Swine farm associated C. difficile ST11 carried fewer virulence genes than ST11 strains collected from NCBI database. It is critical to monitor the evolution of C. difficile strains to understand their changing characteristics, host-switching, and develop effective control and prevention strategies.

Time-calibrated phylogenetic and chromosomal mobilome analyses of Staphylococcus aureus CC398 reveal geographical and host-related evolution

An international collection of Staphylococcus aureus of clonal complex (CC) 398 from diverse hosts spanning all continents and a 30 year-period is studied based on whole-genome sequencing (WGS) data. The collection consists of publicly available genomic data from 2994 strains and 134 recently sequenced Swiss methicillin-resistant S. aureus (MRSA) CC398 strains. A time-calibrated phylogeny reveals the presence of distinct phylogroups present in Asia, North and South America and Europe. European MRSA diverged from methicillin-susceptible S. aureus (MSSA) at the beginning of the 1950s. Two major European phylogroups (EP4 and EP5), which diverged approximately 1974, are the main drivers of MRSA CC398 spread in Europe. Within EP5, an emergent MRSA lineage spreading among the European horse population (EP5-Leq) diverged approximately 1996 from the pig lineage (EP5-Lpg), and also contains human-related strains. EP5-Leq is characterized by staphylococcal cassette chromosome mec (SCCmec) IVa and spa type t011 (CC398-IVa-t011), and EP5-Lpg by CC398-SCCmecVc-t011. The lineage-specific antibiotic resistance and virulence gene patterns are mostly mediated by the acquisition of mobile genetic elements like SCCmec, S. aureus Genomic Islands (SaGIs), prophages and transposons. Different combinations of virulence factors are present on S. aureus pathogenicity islands (SaPIs), and novel antimicrobial resistance gene containing elements are associated with certain lineages expanding in Europe. This WGS-based analysis reveals the actual evolutionary trajectory and epidemiological trend of the international MRSA CC398 population considering host, temporal, geographical and molecular factors. It provides a baseline for global WGS-based One-Health studies of adaptive evolution of MRSA CC398 as well as for local outbreak investigations.

Self-assembling viral histones are evolutionary intermediates between archaeal and eukaryotic nucleosomes

Nucleosomes are DNA-protein complexes composed of histone proteins that form the basis of eukaryotic chromatin. The nucleosome was a key innovation during eukaryotic evolution, but its origin from histone homologues in Archaea remains unclear. Viral histone repeats, consisting of multiple histone paralogues within a single protein, may reflect an intermediate state. Here we examine the diversity of histones encoded by Nucleocytoviricota viruses. We identified 258 histones from 168 viral metagenomes with variable domain configurations including histone singlets, doublets, triplets and quadruplets, the latter comprising the four core histones arranged in series. Viral histone repeats branch phylogenetically between Archaea and eukaryotes and display intermediate functions in Escherichia coli, self-assembling into eukaryotic-like nucleosomes that stack into archaeal-like oligomers capable of impacting genomic activity and condensing DNA. Histone linkage also facilitates nucleosome formation, promoting eukaryotic histone assembly in E. coli. These data support the hypothesis that viral histone repeats originated in stem-eukaryotes and that nucleosome evolution proceeded through histone repeat intermediates.

Complete chloroplast genomes of the hemiparasitic genus Cymbaria: Insights into comparative analysis, development of molecular markers, and phylogenetic relationships

The hemiparasitic tribe Cymbarieae (Orobanchaceae) plays a crucial role in elucidating the initial stage of the transition from autotrophism to heterotrophism. However, the complete chloroplast genome of the type genus Cymbaria has yet to be reported. In addition, the traditional Mongolian medicine Cymbaria daurica is frequently subjected to adulteration or substitution because of the minor morphological differences with Cymbaria mongolica. In this study, the complete chloroplast genomes of the two Cymbaria species were assembled and annotated, and those of other published 52 Orobanchaceae species were retrieved for comparative analyses. We found that the Cymbaria chloroplast genomes are characterized by pseudogenization or loss of stress-relevant genes (ndh) and a unique rbcL-matK inversion. Unlike the high variability observed in holoparasites, Cymbaria and other hemiparasites exhibit high similarity to autotrophs in genome size, guanine-cytosine (GC) content, and intact genes. Notably, four pairs of specific DNA barcodes were developed and validated to distinguish the medicinal herb from its adulterants. Phylogenetic analyses revealed that the genus Cymbaria and the Schwalbea-Siphonostegia clade are grouped into the tribe Cymbarieae, which forms a sister clade to the remaining Orobanchaceae parasitic lineages. Moreover, the diversification of monophyletic Cymbaria occurred during the late Miocene (6.72 Mya) in the Mongol-Chinese steppe region. Our findings provide valuable genetic resources for studying the phylogeny of Orobanchaceae and plant parasitism, and genetic tools to validate the authenticity of the traditional Mongolian medicine "Xinba.".

Genetic and epigenetic signatures for improved breeding of cultivated blueberry

Blueberry belongs to the Vaccinium genus and is a highly popular fruit crop with significant economic importance. It was not until the early twentieth century that they began to be domesticated through extensive interspecific hybridization. Here, we collected 220 Vaccinium accessions from various geographical locations, including 154 from the United States, 14 from China, eight from Australia, and 29 from Europe and other countries, comprising 164 Vaccinium corymbosum, 15 Vaccinium ashei, 10 lowbush blueberries, seven half-high blueberries, and others. We present the whole-genome variation map of 220 accessions and reconstructed the hundred-year molecular history of interspecific hybridization of blueberry. We focused on the two major blueberry subgroups, the northern highbush blueberry (NHB) and southern highbush blueberry (SHB) and identified candidate genes that contribute to their distinct traits in climate adaptability and fruit quality. Our analysis unveiled the role of gene introgression from Vaccinium darrowii and V. ashei into SHB in driving the differentiation between SHB and NHB, potentially facilitating SHB's adaptation to subtropical environments. Assisted by genome-wide association studies, our analysis suggested VcTBL44 as a pivotal gene regulator governing fruit firmness in SHB. Additionally, we conducted whole-genome bisulfite sequencing on nine NHB and 12 SHB cultivars, and characterized regions that are differentially methylated between the two subgroups. In particular, we discovered that the β-alanine metabolic pathway genes were enriched for DNA methylation changes. Our study provides high-quality genetic and epigenetic variation maps for blueberry, which offer valuable insights and resources for future blueberry breeding.

Unraveling the potential of uninvestigated thermoalkaliphilic lipases by molecular docking and molecular dynamic simulation: an in silico characterization study

Thermoalkaliphilic lipase enzymes are mostly favored for use in the detergent industry. While there has been considerable research on Geobacillus lipases, a significant portion of these enzymes remains unexplored or undocumented in the scientific literature. This work performed in silico phylogeny, sequence alignment, structural and enzyme-substrate interaction analyses of the five thermoalkaliphilic lipases belonging to different Geobacillus species (Geobacillus stearothermophilus lipase = GsLip, Geobacillus sp. B4113_201601 lipase = Gb4Lip, Geobacillus kaustophilus HTA426 lipase = GkLip, Geobacillus sp. SP22 lipase = GspLip, Geobacillus sp. NTU 03 lipase = GntLip). For this purpose, unreviewed enzyme sequences of five Geobacillus thermoalkaliphilic lipases were analyzed at sequence and phylogeny levels. 3D homology enzyme models were built, validated, and investigated by different bioinformatics tools. The ligand interactions screening using seven para-nitrophenyl (pNP) esters and enzyme-ligand interactions were analyzed on Gb4Lip:pNP-C12 and BTL2:pNP-C12 by MD simulation. Biophysicochemical characteristic analysis showed that Gb4Lip had a theoretical T m value of above 65 ºC, and a higher aliphatic index indicating greater thermal stability. Sequence alignment showed a hydrophilic threonine in the α6 helix of Gb4Lip, indicating high enzymatic activity. A normalized temperature factor B (B'-factor) analysis showed that the lid domains of five lipases significantly possessed lower B'-factor values, compared to G. thermocatenulatus lipase 2 (BTL2), indicating that they had higher rigidity. Molecular docking results indicated that the five lipases had the highest binding affinity toward pNP-C12. The RMSF investigation revealed that the thermostability of Gb4Lip is influenced by specific molecular elements: D202-S203 within the αB region of the lid domain, and E274-Q275 within the b3 strand, as well as W278 in the b3-b4 loop, and H282 in the b4 strand of the Ca2+-binding region. MD simulation analysis showed that catalytic residue S114 and at least one oxyanion hole residue (F17 and/or Q114) in Gb4Lip frequently formed hydrogen bonds with the pNP-C12 ligand at 343 K and 348 K throughout the simulation process, indicating that Gb4Lip might catalyze relatively long-chain ligand pNP-C12 with high performance. In conclusion, Gb4Lip might be a more suitable candidate as the detergent additive. In addition, this investigation can offer valuable perspectives on Family I.5 lipases such as Gb4Lip for future exploration in the field of protein engineering.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04023-5.

Genomic study of European Clostridioides difficile ribotype 002/sequence type 8

Clostridioides difficile has significant clinical importance as a leading cause of healthcare-associated infections, with symptoms ranging from mild diarrhoea to severe colitis, and possible life-threatening complications. C. difficile ribotype (RT) 002, mainly associated with MLST sequence type (ST) 8, is one of the most common RTs found in humans. This study aimed at investigating the genetic characteristics of 537 C. difficile genomes of ST8/RT002. To this end, we sequenced 298 C. difficile strains representing a new European genome collection, with strains from Germany, Denmark, France and Portugal. These sequences were analysed against a global dataset consisting of 1,437 ST8 genomes available through Enterobase. Our results showed close genetic relatedness among the studied ST8 genomes, a diverse array of antimicrobial resistance (AMR) genes and the presence of multiple mobile elements. Notably, the pangenome analysis revealed an open genomic structure. ST8 shows relatively low overall variation. Thus, clonal isolates were found across different One Health sectors (humans, animals, environment and food), time periods, and geographical locations, suggesting the lineage's stability and a universal environmental source. Importantly, this stability did not hinder the acquisition of AMR genes, emphasizing the adaptability of this bacterium to different selective pressures. Although only 2.4 % (41/1,735) of the studied genomes originated from non-human sources, such as animals, food, or the environment, we identified 9 cross-sectoral core genome multilocus sequence typing (cgMLST) clusters. Our study highlights the importance of ST8 as a prominent lineage of C. difficile with critical implications in the context of One Health. In addition, these findings strongly support the need for continued surveillance and investigation of non-human samples to gain a more comprehensive understanding of the epidemiology of C. difficile.

Genome-wide identification of the LEA gene family in Panax ginseng: Evidence for the role of PgLEA2-50 in plant abiotic stress response

Ginseng frequently encounters environmental stress during its growth and development. Late Embryogenesis Abundant (LEA) proteins play a crucial role in combating adversity stress, particularly against abiotic challenges In this study, 107 LEA genes from ginseng, spanning eight subfamilies, were identified, demonstrating significant evolutionary conservation, with the LEA2 subfamily being most prominent. Gene duplication events, primarily segmental duplications, have played a major role in the expansion of the LEA gene family, which has undergone strong purifying selection. PgLEAs were unevenly distributed across 22 chromosomes, with each subfamily featuring unique structural domains and conserved motifs. PgLEAs were expressed in various tissues, exhibiting distinct variations in abundance and tissue specificity. Numerous regulatory cis-elements, related to abiotic stress and hormones, were identified in the promoter region. Additionally, PgLEAs were regulated by a diverse array of abiotic stress-related transcription factors. A total of 35 PgLEAs were differentially expressed following treatments with ABA, GA, and IAA. Twenty-three PgLEAs showed significant but varied responses to drought, extreme temperatures, and salinity stress. The transformation of tobacco with the key gene PgLEA2-50 enhanced osmoregulation and antioxidant levels in transgenic lines, improving their resistance to abiotic stress. This study offers insights into functional gene analysis, focusing on LEA proteins, and establishes a foundational framework for research on ginseng's resilience to abiotic stress.

Recombination and repeat-induced point mutation landscapes reveal trade-offs between the sexual and asexual cycles of Magnaporthe oryzae

The fungal disease caused by Magnaporthe oryzae is one of the most devastating diseases that endanger many crops worldwide. Evidence shows that sexual reproduction can be advantageous for fungal diseases as hybridization facilitates host-jumping. However, the pervasive clonal lineages of M. oryzae observed in natural fields contradict this expectation. A better understanding of the roles of recombination and the fungi-specific repeat-induced point mutation (RIP) in shaping its evolutionary trajectory is essential to bridge this knowledge gap. Here we systematically investigate the RIP and recombination landscapes in M. oryzae using a whole genome sequencing data from 252 population samples and 92 cross progenies. Our data reveal that the RIP can robustly capture the population history of M. oryzae, and we provide accurate estimations of the recombination and RIP rates across different M. oryzae clades. Significantly, our results highlight a parent-of-origin bias in both recombination and RIP rates, tightly associating with their sexual potential and variations of effector proteins. This bias suggests a critical trade-off between generating novel allelic combinations in the sexual cycle to facilitate host-jumping and stimulating transposon-associated diversification of effectors in the asexual cycle to facilitate host coevolution. These findings provide unique insights into understanding the evolution of blast fungus.

Characterization of drug resistance genes in Indian Plasmodium falciparum and Plasmodium vivax field isolates

One of the major challenges for malaria control and elimination is the spread and emergence of antimalarial drug resistance. Mutations in Plasmodium falciparum (Pf) and Plasmodium vivax (Pv) field isolates for five drug resistance genes viz. crt, mdr1, dhps, dhfr and kelch known to confer resistance to choloroquine (CQ), sulfadoxine-pyrimethamine (SP) and artemisinin (ART) and its derivatives were analyzed. A total of 342 symptomatic isolates of P. falciparum (Pf) and P. vivax (Pv) from 1993 to 2014 were retrieved from malaria parasite repository at National Institute of Malaria Research (NIMR). Sample DNA was extracted from dried blood spots and various targeted single nucleotide polymorphisms (SNPs) associated with antimalarial drug resistance were analysed for these isolates. 72S (67.7%) and 76T (83.8%) mutations along with SVMNT haplotype (67.7%) predominated the study population for Pfcrt. The most prevalent SNPs were 108N (73.2%) and 437G (24.8%) and the most prevalent haplotypes were ACNRNI (51.9%) and SAKAA (74.5%) in Pfdhfr and Pfdhps respectively. Only two mutations in Pfmdr1, 86Y (26.31%) and 184F (56.26%), were seen frequently in our study population. No mutations associated with Pfk13 were observed. For Pv, all the studied isolates showed two Pvdhps mutations, 383G and 553G, and two Pfdhfr mutations, 58R and 117N. Similarly, three mutations, viz. 958M, 908L and 1076L were found in Pvmdr1. No variations were observed in Pvcrt-o and Pvk12 genes. Overall, our study demonstrates an increase in mutations associated with SP resistance in both Pf and Pv, however, no single nucleotide polymorphisms (SNPs) associated with ART resistance have been observed for either species. Various SNPs associated with CQ resistance were seen in Pf; whereas only Pvmdr1 associated resistant SNPs were observed in Pv. Therefore, molecular characterization of drug resistance genes is essential for timely monitoring and prevention of malaria by identifying the circulating drug resistant parasites in the country.

Whole-genome sequencing, multilocus sequence typing, and resistance mechanism of the carbapenem-resistant Pseudomonas aeruginosa in China

Pseudomonas aeruginosa is a significant pathogen responsible for severe multisite infections with high morbidity and mortality rates. This study analyzed carbapenem-resistant Pseudomonas aeruginosa (CRPA) at a tertiary hospital in Shandong, China, using whole-genome sequencing (WGS). The objective was to explore the mechanisms and molecular characteristics of carbapenem resistance. A retrospective analysis of 91 isolates from January 2022 to March 2023 was performed, which included strain identification and antimicrobial susceptibility testing. WGS was utilized to determine the genome sequences of these CRPA strains, and the species were precisely identified using average nucleotide identification (ANI), with further analysis on multilocus sequence typing and strain relatedness. Some strains were found to carry the ampD and oprD genes, while only a few harbored carbapenemase genes or related genes. Notably, all strains possessed the mexA, mexE, and mexX genes. The major lineage identified was ST244, followed by ST235. The study revealed a diverse array of carbapenem resistance mechanisms among hospital isolates, differing from previous studies in mainland China. It highlighted that carbapenem resistance is not due to a single mechanism but rather a combination of enzyme-mediated resistance, AmpC overexpression, OprD dysfunction, and efflux pump overexpression. This research provides valuable insights into the evolutionary mechanisms and molecular features of CRPA resistance in this region, aiding in the national prevention and control of CRPA, and offering references for targeting and developing new drugs.

Functional identification of specialized diterpene synthases from Chamaecyparis obtusa and C. obtusa var. formosana to illustrate the putative evolution of diterpene synthases in Cupressaceae

Chamaecyparis obtusa and C. obtusa var. formosana of the Cupressaceae family are well known for their fragrance and excellent physical properties. To investigate the biosynthesis of unique diterpenoid compounds, diterpene synthase genes for specialized metabolite synthesis were cloned from C. obtusa and C. obtusa var. formosana. Using an Escherichia coli co-expression system, eight diterpene synthases (diTPSs) were characterized. CoCPS and CovfCPS are class II monofunctional (+)-copalyl diphosphate synthases [(+)-CPSs]. Class I monofunctional CoLS and CovfLS convert (+)-copalyl diphosphate [(+)-CPP] to levopimaradiene, CoBRS, CovfBRS1, and CovfBRS3 convert (+)-CPP to (-)-beyerene, and CovfSDS converts (+)-CPP to (-)-sandaracopimaradiene. These enzymes are all monofunctional diterpene syntheses in Cupressaceae family of gymnosperm, and differ from those in Pinaceae. The discovery of the enzyme responsible for the biosynthesis of tetracyclic diterpene (-)-beyerene was characterized for the first time. Diterpene synthases with different catalytic functions exist in closely related species within the Cupressaceae family, indicating that this group of monofunctional diterpene synthases is particularly prone to the evolution of new functions and development of species-specific specialized diterpenoid constituents.

Evidence of Lactobacillus strains shared between the female urinary and vaginal microbiota

Lactobacillus species are common inhabitants of the 'healthy' female urinary and vaginal communities, often associated with a lack of symptoms in both anatomical sites. Given identification by prior studies of similar bacterial species in both communities, it has been hypothesized that the two microbiotas are in fact connected. Here, we carried out whole-genome sequencing of 49 Lactobacillus strains, including 16 paired urogenital samples from the same participant. These strains represent five different Lactobacillus species: L. crispatus, L. gasseri, L. iners, L. jensenii, and L. paragasseri. Average nucleotide identity (ANI), alignment, single-nucleotide polymorphism (SNP), and CRISPR comparisons between strains from the same participant were performed. We conducted simulations of genome assemblies and ANI comparisons and present a statistical method to distinguish between unrelated, related, and identical strains. We found that 50 % of the paired samples have identical strains, evidence that the urinary and vaginal communities are connected. Additionally, we found evidence of strains sharing a common ancestor. These results establish that microbial sharing between the urinary tract and vagina is not limited to uropathogens. Knowledge that these two anatomical sites can share lactobacilli in females can inform future clinical approaches.

Genome-Wide Identification of the WRKY Gene Family in Four Cotton Varieties and the Positive Role of GhWRKY31 in Response to Salt and Drought Stress

The WRKY gene family is ubiquitously distributed in plants, serving crucial functions in stress responses. Nevertheless, the structural organization and evolutionary dynamics of WRKY genes in cotton have not been fully elucidated. In this study, a total of 112, 119, 217, and 222 WRKY genes were identified in Gossypium arboreum, Gossypium raimondii, Gossypium hirsutum, and Gossypium barbadense, respectively. These 670 WRKY genes were categorized into seven distinct subgroups and unequally distributed across chromosomes. Examination of conserved motifs, domains, cis-acting elements, and gene architecture collectively highlighted the evolutionary conservation and divergence within the WRKY gene family in cotton. Analysis of synteny and collinearity further confirmed instances of expansion, duplication, and loss events among WRKY genes during cotton evolution. Furthermore, GhWRKY31 transgenic Arabidopsis exhibited heightened germination rates and longer root lengths under drought and salt stress. Silencing GhWRKY31 in cotton led to reduced levels of ABA, proline, POD, and SOD, along with downregulated expression of stress-responsive genes. Yeast one-hybrid and molecular docking assays confirmed the binding capacity of GhWRKY31 to the W box of GhABF1, GhDREB2, and GhRD29. The findings collectively offer a systematic and comprehensive insight into the evolutionary patterns of cotton WRKYs, proposing a suitable regulatory framework for developing cotton cultivars with enhanced resilience to drought and salinity stress.

Public health aspects of Vibrio spp. related to the consumption of seafood in the EU

Vibrio parahaemolyticus, Vibrio vulnificus and non-O1/non-O139 Vibrio cholerae are the Vibrio spp. of highest relevance for public health in the EU through seafood consumption. Infection with V. parahaemolyticus is associated with the haemolysins thermostable direct haemolysin (TDH) and TDH-related haemolysin (TRH) and mainly leads to acute gastroenteritis. V. vulnificus infections can lead to sepsis and death in susceptible individuals. V. cholerae non-O1/non-O139 can cause mild gastroenteritis or lead to severe infections, including sepsis, in susceptible individuals. The pooled prevalence estimate in seafood is 19.6% (95% CI 13.7-27.4), 6.1% (95% CI 3.0-11.8) and 4.1% (95% CI 2.4-6.9) for V. parahaemolyticus, V. vulnificus and non-choleragenic V. cholerae, respectively. Approximately one out of five V. parahaemolyticus-positive samples contain pathogenic strains. A large spectrum of antimicrobial resistances, some of which are intrinsic, has been found in vibrios isolated from seafood or food-borne infections in Europe. Genes conferring resistance to medically important antimicrobials and associated with mobile genetic elements are increasingly detected in vibrios. Temperature and salinity are the most relevant drivers for Vibrio abundance in the aquatic environment. It is anticipated that the occurrence and levels of the relevant Vibrio spp. in seafood will increase in response to coastal warming and extreme weather events, especially in low-salinity/brackish waters. While some measures, like high-pressure processing, irradiation or depuration reduce the levels of Vibrio spp. in seafood, maintaining the cold chain is important to prevent their growth. Available risk assessments addressed V. parahaemolyticus in various types of seafood and V. vulnificus in raw oysters and octopus. A quantitative microbiological risk assessment relevant in an EU context would be V. parahaemolyticus in bivalve molluscs (oysters), evaluating the effect of mitigations, especially in a climate change scenario. Knowledge gaps related to Vibrio spp. in seafood and aquatic environments are identified and future research needs are prioritised.

Single-fly genome assemblies fill major phylogenomic gaps across the Drosophilidae Tree of Life

Long-read sequencing is driving rapid progress in genome assembly across all major groups of life, including species of the family Drosophilidae, a longtime model system for genetics, genomics, and evolution. We previously developed a cost-effective hybrid Oxford Nanopore (ONT) long-read and Illumina short-read sequencing approach and used it to assemble 101 drosophilid genomes from laboratory cultures, greatly increasing the number of genome assemblies for this taxonomic group. The next major challenge is to address the laboratory culture bias in taxon sampling by sequencing genomes of species that cannot easily be reared in the lab. Here, we build upon our previous methods to perform amplification-free ONT sequencing of single wild flies obtained either directly from the field or from ethanol-preserved specimens in museum collections, greatly improving the representation of lesser studied drosophilid taxa in whole-genome data. Using Illumina Novaseq X Plus and ONT P2 sequencers with R10.4.1 chemistry, we set a new benchmark for inexpensive hybrid genome assembly at US $150 per genome while assembling genomes from as little as 35 ng of genomic DNA from a single fly. We present 183 new genome assemblies for 179 species as a resource for drosophilid systematics, phylogenetics, and comparative genomics. Of these genomes, 62 are from pooled lab strains and 121 from single adult flies. Despite the sample limitations of working with small insects, most single-fly diploid assemblies are comparable in contiguity (>1 Mb contig N50), completeness (>98% complete dipteran BUSCOs), and accuracy (>QV40 genome-wide with ONT R10.4.1) to assemblies from inbred lines. We present a well-resolved multi-locus phylogeny for 360 drosophilid and 4 outgroup species encompassing all publicly available (as of August 2023) genomes for this group. Finally, we present a Progressive Cactus whole-genome, reference-free alignment built from a subset of 298 suitably high-quality drosophilid genomes. The new assemblies and alignment, along with updated laboratory protocols and computational pipelines, are released as an open resource and as a tool for studying evolution at the scale of an entire insect family.

Diversity, abundance, and expression of proteorhodopsin genes in the northern South China Sea

Proteorhodopsins have been suggested as an important strategy among phototrophs to capture solar energy in marine environments. The goals of this study was to investigate the diversity of proteorhodopsin genes and to explore their abundance, distribution, and expression in the coastal surface waters of the northern South China Sea, one of the largest marginal seas of the western North Pacific Ocean. Using 21 metagenomes, we recovered proteorhodopsin genes from a wide range of prokaryotic taxa, and chlorophyll a contributed significantly to the community composition of proteorhodopsin-containing microbes. Most proteorhodopsin sequences were predicted to encode green light-absorbing proton pumps and green light-absorbing proteorhodopsin genes were more abundant than blue-absorbing ones. The variations in the conserved residues involved in ion pumping and several uncharacterized proteorhodopsins were observed. The gene abundance pattern of proteorhodopsin types were significantly influenced by the levels of total organic carbon and soluble reactive phosphorus. Gene expression analysis confirmed the importance of proteorhodopsin-based phototrophy and revealed different expressional patterns among major phyla. In tandem, we screened 2295 metagenome-assembled genomes to describe the taxonomic distribution of proteorhodopsins. Bacteroidota are the key lineages encoding proteorhodopsins, but proteorhodopsins were predicated from members of Proteobacteria, Marinisomatota, Myxococcota, Verrucomicrobiota and Thermoplasmatota. Our study expanded the diversity of proteorhodopsins and improve our understanding on the significance of proteorhodopsin-mediated phototrophy in the marine ecosystem.

The origins of haplotype 58 (H58) Salmonella enterica serovar Typhi

Antimicrobial resistance (AMR) poses a serious threat to the clinical management of typhoid fever. AMR in Salmonella Typhi (S. Typhi) is commonly associated with the H58 lineage, a lineage that arose comparatively recently before becoming globally disseminated. To better understand when and how H58 emerged and became dominant, we performed detailed phylogenetic analyses on contemporary genome sequences from S. Typhi isolated in the period spanning the emergence. Our dataset, which contains the earliest described H58 S. Typhi organism, indicates that ancestral H58 organisms were already multi-drug resistant (MDR). These organisms emerged spontaneously in India in 1987 and became radially distributed throughout South Asia and then globally in the ensuing years. These early organisms were associated with a single long branch, possessing mutations associated with increased bile tolerance, suggesting that the first H58 organism was generated during chronic carriage. The subsequent use of fluoroquinolones led to several independent mutations in gyrA. The ability of H58 to acquire and maintain AMR genes continues to pose a threat, as extensively drug-resistant (XDR; MDR plus resistance to ciprofloxacin and third generation cephalosporins) variants, have emerged recently in this lineage. Understanding where and how H58 S. Typhi originated and became successful is key to understand how AMR drives successful lineages of bacterial pathogens. Additionally, these data can inform optimal targeting of typhoid conjugate vaccines (TCVs) for reducing the potential for emergence and the impact of new drug-resistant variants. Emphasis should also be placed upon the prospective identification and treatment of chronic carriers to prevent the emergence of new drug resistant variants with the ability to spread efficiently.

Structural characterization of macro domain-containing Thoeris antiphage defense systems

Thoeris defense systems protect bacteria from infection by phages via abortive infection. In these systems, ThsB proteins serve as sensors of infection and generate signaling nucleotides that activate ThsA effectors. Silent information regulator and SMF/DprA-LOG (SIR2-SLOG) containing ThsA effectors are activated by cyclic ADP-ribose (ADPR) isomers 2'cADPR and 3'cADPR, triggering abortive infection via nicotinamide adenine dinucleotide (NAD+) depletion. Here, we characterize Thoeris systems with transmembrane and macro domain (TM-macro)-containing ThsA effectors. We demonstrate that ThsA macro domains bind ADPR and imidazole adenine dinucleotide (IAD), but not 2'cADPR or 3'cADPR. Combining crystallography, in silico predictions, and site-directed mutagenesis, we show that ThsA macro domains form nucleotide-induced higher-order oligomers, enabling TM domain clustering. We demonstrate that ThsB can produce both ADPR and IAD, and we identify a ThsA TM-macro-specific ThsB subfamily with an active site resembling deoxy-nucleotide and deoxy-nucleoside processing enzymes. Collectively, our study demonstrates that Thoeris systems with SIR2-SLOG and TM-macro ThsA effectors trigger abortive infection via distinct mechanisms.

Brainwide silencing of prion protein by AAV-mediated delivery of an engineered compact epigenetic editor

Prion disease is caused by misfolding of the prion protein (PrP) into pathogenic self-propagating conformations, leading to rapid-onset dementia and death. However, elimination of endogenous PrP halts prion disease progression. In this study, we describe Coupled Histone tail for Autoinhibition Release of Methyltransferase (CHARM), a compact, enzyme-free epigenetic editor capable of silencing transcription through programmable DNA methylation. Using a histone H3 tail-Dnmt3l fusion, CHARM recruits and activates endogenous DNA methyltransferases, thereby reducing transgene size and cytotoxicity. When delivered to the mouse brain by systemic injection of adeno-associated virus (AAV), Prnp-targeted CHARM ablates PrP expression across the brain. Furthermore, we have temporally limited editor expression by implementing a kinetically tuned self-silencing approach. CHARM potentially represents a broadly applicable strategy to suppress pathogenic proteins, including those implicated in other neurodegenerative diseases.

A novel Enterococcus faecium phage EF-M80: unveiling the effects of hydrogel-encapsulated phage on wound infection healing

Background

Enterococcus faecium is one of the members of ESKAPE pathogens. Due to its resistance to antimicrobial agents, treating this bacterium has become challenging. The development of innovative approaches to combat antibiotic resistance is necessary. Phage therapy has emerged as a promising method for curing antibiotic-resistant bacteria.

Methods

In this study, E. faecium phages were isolated from wastewater. Phage properties were characterized through in vitro assays (e.g. morphological studies, and physicochemical properties). In addition, whole genome sequencing was performed. A hydrogel-based encapsulated phage was obtained and its structure characteristics were evaluated. Wound healing activity of the hydrogel-based phage was assessed in a wound mice model.

Results

The purified phage showed remarkable properties including broad host range, tolerance to high temperature and pH and biofilm degradation feature as a stable and reliable therapeutic agent. Whole genome sequencing revealed that the genome of the EF-M80 phage had a length of 40,434 bp and harbored 65 open reading frames (ORFs) with a GC content of 34.9% (GenBank accession number is OR767211). Hydrogel-based encapsulated phage represented an optimized structure. Phage-loaded hydrogel-treated mice showed that the counting of neutrophils, fibroblasts, blood vessels, hair follicles and percentage of collagen growth were in favor of the wound healing process in the mice model.

Conclusion

These findings collectively suggest the promising capability of this phage-based therapeutic strategy for the treatment of infections associated with the antibiotic-resistant E. faecium. In the near future, we hope to expect the presence of bacteriophages in the list of antibacterial compounds used in the clinical settings.