hybridSPAdes: an algorithm for hybrid assembly of short and long reads

Characterization of a new lytic bacteriophage (vB_KpnM_KP1) targeting Klebsiella pneumoniae strains associated with nosocomial infections

A new bacteriophage, vB_KpnM_KP1, was identified and characterized, exhibiting a strong lytic effect on Klebsiella pneumoniae. Host range analysis revealed its effectiveness against 77.4% of clinical strains, achieving complete lysis of those associated with urinary tract infections (UTIs). Phage stability tests demonstrated that vB_KpnM_KP1 remained stable at neutral pH and across all tested temperatures. However, inactivation was observed at high ethanol concentrations and extreme pH levels. Transmission electron microscopy (TEM) analysis identified vB_KpnM_KP1 as a Myo-type phage with an icosahedral head and a contractile tail. Moreover, genome annotation of vB_KpnM_KP1 revealed a linear DNA genome of 174,802 bp, containing 307 open reading frames. Functional predictions suggest the presence of genes involved in DNA replication, transcription, morphogenesis, and cell lysis. Phylogenetic analysis classified vB_KpnM_KP1 within the Slopekvirus genus of the Straboviridae family, showing high sequence identity with phages that infect Enterobacter, Escherichia and Klebsiella species. These findings highlight the potential of phage vB_KpnM_KP1 as an alternative treatment for multidrug-resistant K. pneumoniae infections, particularly in UTIs, while offering valuable insights into its stability and genetic composition.

Cardinium disrupts Wolbachia-host dynamics in the domestic mite Tyrophagus putrescentiae: evidence from manipulative experiments

We investigated the tripartite interactions between two intracellular bacterial symbionts, Cardinium and Wolbachia in Tyrophagus putrescentiae. Cultures of Tyrophagus putrescentiae are typically single-infected by one intracellular symbiont. However, co-infection can be experimentally induced by mixing single-infected cultures, resulting in 10% of mite individuals being double-infected (Cardinium + Wolbachia) and a corresponding reduction in host fitness. Here, we assembled the genomes of Cardinium and Wolbachia and analyzed their gene expression in parental single-infected and mixed mite cultures using population-level samples (ranging from 7,500 to 10,000 mites). Wolbachia interacts more extensively with its mite host than Cardinium in single-infected cultures. However, in mixed cultures, (i) Wolbachia exhibited reduced regulation of the host compared with Cardinium; (ii) the gene expression profile of Cardinium shifted, increasing its interactions with the host, whereas the gene expression profile of Wolbachia remained unchanged; and (iii) Wolbachia genes exhibited a loss of interactions with mite gene expression, as indicated by reduced correlations (for example with host MAPK, endocytosis, and calcium signaling pathways). The experiments show that at the mite population level, symbiont infection disrupts gene expression interaction between the two symbionts and their host in different ways. Wolbachia was more influenced by Cardinium gene expression than vice versa. Cardinium can inhibit the growth of Wolbachia by disrupting its interaction with the host, leading to a loss of Wolbachia's influence on mite immune and regulatory pathways. The reasons for responses are due to co-infection or the reduced frequency of Wolbachia single-infected individuals due to the analyses of population-level samples.

IMPORTANCE

We found that Cardinium disrupts the interaction between Wolbachia and mite host. In Wolbachia single-infected cultures, strong correlations exist between symbiont and host gene expressions. Interestingly, although Cardinium can also interact with the host, this interaction appears weaker compared with Wolbachia in single-infected cultures. These results suggest that both symbionts affect mite host gene expression, particularly in immune and regulatory pathways. In mixed samples, Cardinium appears to outcompete Wolbachia by disrupting its host interaction. It indicates competition between these two intracellular symbionts in mite populations. Wolbachia belongs to a mite-specific supergroup Q, distinct from the more commonly studied Wolbachia supergroups. As these mite-specific bacteria exhibit pathogen-blocking effects, our findings may have relevance for other systems, such as ticks and tick-borne diseases. The study sheds light on intracellular symbiont interaction within a novel mite-symbiont model.

Contaminated drinking water facilitates Escherichia coli strain-sharing within households in urban informal settlements

Identifying bacterial transmission pathways is crucial to inform strategies that limit the spread of pathogenic and antibiotic-resistant bacteria. Here we assessed Escherichia coli strain-sharing and overlap of antibiotic resistance genes (ARGs) across humans, poultry, canines, soil, and drinking water within and between households in urban informal settlements in Nairobi, Kenya. We collected 321 samples from 50 households with half having access to chlorinated water. We performed Pooling Isolated Colonies-seq, which sequences pools of up to five E. coli colonies per sample to capture strain diversity. Pooling Isolated Colonies-seq captured 1,516 colonies and identified 154 strain-sharing events, overcoming limitations of single-isolate sequencing and metagenomics. Within households, strain-sharing rates and resistome similarities across sample types were strongly correlated, suggesting clonal transmission of ARGs. E. coli isolated from the environment carried clinically relevant ARGs. Strain-sharing was rare between animals and humans but frequent between humans and drinking water. E. coli-contaminated stored drinking water was associated with higher human-human strain-sharing within households. These results suggest that contaminated drinking water facilitates human to human strain-sharing, and water treatment can disrupt transmission.

Genomic and Phenotypic Characterization of a Novel Virulent Strain of Cyvirus cyprinidallo2 Originating from an Outbreak in The Netherlands

Cyvirus cyprinidallo2 (CyHV-2) is the causative agent of herpesviral hematopoietic necrosis in several economically important farmed freshwater fish species of the genus Carassius. Despite several CyHV-2 strains being isolated and fully sequenced, there is a lack of detailed characterization and consistent information on strains that exhibit high virulence in adult goldfish through viral challenge by immersion, particularly in the context of European strains and host populations. Strains that can cause highly virulent disease via this inoculation route are much more compatible with experimental designs that are representative of natural infection; thus, their utilization provides greater biological relevance. Consequently, in this study, we isolated three novel strains of CyHV-2 (designated NL-1, NL-2, and NL-3), originating from outbreaks in The Netherlands. Full-length genome sequencing and phylogenetic analyses revealed that these newly isolated strains are distinct from known strains and from each other. Significant differences were observed between the strains, in terms of in vitro growth kinetics, with NL-2 exhibiting stable passaging and superior fitness in vitro. Importantly, the challenge of adult Shubunkin goldfish with the NL-2 strain via immersion (2000 PFU/mL) induced an average mortality of ~40%, while parallel experiments with the CyHV-2 reference strain ST-J1 resulted in no mortality. Taken together, this study represents the characterization of a new CyHV-2 in vivo infection model, much more compatible with experimental designs that are required to be representative of natural infection. This model will be extremely useful in many aspects of CyHV-2 research in the future. Importantly, the genetic and phenotypic characterization performed in this study generates hypotheses on the potential roles of CyHV-2 genes in adaptation of the virus in vitro or in vivo.

A novel Erwiniaceae gut symbiont modulates gene expression of the intracellular bacterium Cardinium in the stored product mite Tyrophagus putrescentiae

We examined host and bacterial gene expression profiles in the stored product mite Tyrophagus putrescentiae co-infected with Wolbachia (wTPut) and Cardinium (cTPut) while varying the presence of the Erwiniaceae symbiont (SLS). SLS, a novel symbiont in the family Erwiniaceae, with a genome size of 1.7 Mb, is found in 16% of mite species in infected cultures. In addition, SLS was detected in mite feces but not in their eggs. Although Wolbachia expression remained unchanged, the presence or absence of SLS significantly affected Cardinium expression. It indicated that the effect of Wolbachia on SLS was neutral. In SLS-positive samples, Cardinium exhibited 29 upregulated and 48 downregulated genes compared to SLS-negative samples. Furthermore, Cardinium gene expression strongly correlated with mite KEGG gene expression in SLS-positive samples. Positive Spearman's correlations between Cardinium gene expression and mite KEGG immune and regulatory pathways were doubled in SLS-positive compared to SLS-negative samples. The diversity of expressed genes in the mite host decreased in the presence of SLS. Cardinium had more interacting genes to mite host in SLS-positive samples than without SLS. Transposases are the most affected Cardinium genes, showing upregulation in the presence of SLS. Correlation analyses revealed interactions between Cardinium and SLS via mite immune and regulatory pathways, including lysosome, ubiquitin-mediated proteolysis, PIK3_Akt, and cGMP-PKG. The results showed that Cardinium indirectly affects the gut symbionts of mites.IMPORTANCEThis study introduces a new model to analyze interactions between intracellular bacterial symbionts, gut bacterial symbionts, and their mite hosts. Using gene expression correlations, we investigated how the intracellular Cardinium responds to the novel Erwiniaceae gut symbiont in the mold mite Tyrophagus putrescentiae. The data showed that both mite and Cardinium gene expression are different in the samples with and without Erwiniaceae symbionts. In the presence of Erwiniaceae symbionts, Cardinium increased the interaction with the mite host in terms of changes in gene expression. The mite immune and regulatory pathway gene expression is differently correlated to Cardinium genes in relation to Erwiniaceae symbionts. As a well-known producer of allergens, T. putrescentiae physiology and thus its allergen production are influenced by both symbionts, potentially affecting the release of allergens into human environments.

Novel and diverse features identified in the genomes of bacteria isolated from a hydrothermal vent plume

Hydrothermal vent plumes (HVPs), formed by high-temperature vent emissions, are rich in compounds that support chemosynthesis and serve as reservoirs of microbial diversity and genetic innovation. Through turbulence, mixing, and interaction with subsea currents, vent communities are thought to disperse across ocean basins. In this study, we focused on the plume of the Moytirra hydrothermal vent field, a relatively unexplored site, to investigate its microbial inhabitants. We cultured bacteria from the Moytirra HVP using 11 different media types and performed complete genome sequencing on 12 isolates. Our analyses revealed four putatively novel species from the Thalassobaculum, Sulfitobacter, Idiomarina, and Christiangramia genera. Comparative genomics identified unique genomic islands containing biosynthetic gene clusters, including a novel Non-Ribosomal Peptide Synthetase/Polyketide Synthase cluster, toxin-antitoxin systems, and evidence of horizontal gene transfer facilitated by prophages. These findings underscore the potential of HVPs as a source of novel microbial species and biotechnologically relevant genes, contributing to our understanding of the biodiversity and genetic complexity of these extreme environments.IMPORTANCEHydrothermal vents are dynamic environments that offer unique nutrients for chemosynthetic organisms to drive biology in the deep-sea. The dynamics of these ecosystems are thought to drive genomic innovation in resident populations. Hydrothermal vent plumes (HVPs) mix with surrounding water, carrying local microbiota with them and dispersing for hundreds of kilometers. This study isolated bacteria from a HVP to capture a genomic snapshot of the microbial community, revealing four putatively novel species of bacteria within three taxonomic classes. The addition of these genomes to public databases provides valuable insights into the genomic function, architecture, and novel biosynthetic gene clusters of bacteria found in these extreme environments.

Bacteroides maternus sp. nov., a novel species isolated from human faeces

A novel bacterial strain, MSB163, was isolated from the stool sample of a healthy mother, 4 weeks after giving birth via vaginal delivery. Taxonomic identification tools revealed that MSB163 belongs to the genus Bacteroides, but it is distinct from any currently known species. The closest related species is Bacteroides cellulosilyticus strain BFG- 250, with an average nucleotide identity (fastANI) of 94.51%. The genome length of MSB163 is 6,440,948 bp and the GC content 42.95%. Two plasmids were identified in the whole genome sequence. MSB163 is a Gram-negative, rod-shaped, non-motile anaerobic bacterium. The optimum growth conditions were at 37 °C, pH 7 and 0% (w/v) NaCl. The respiratory quinones were the menaquinones MK- 10 and MK- 11 and C15:0 ANTEISO was the major fatty acid. The predominant polar lipids were phosphatidylethanolamine, diphosphatidylglycerol and phospholipid. According to the taxonomic results and physiological analysis, strain MSB163 represents a novel species of the genus Bacteroides, for which we propose the name Bacteroides maternus, since the type strain was isolated from the stool sample of a mother. B. maternus type strain (MSB163) sequencing can be accessed under the biosample ID SAMN3953129 on NCBI. The strain was deposited on BCCM/LMG Bacteria Collection under the accession number LMG 33,374 and Leibniz Institut DSMZ GMBH under the accession number DSM 117,047.

Kastor: a reference-based comparative approach for assessment and correction of gene-fragmenting errors in long-read assemblies of small genomes

Long read sequencing technologies provide an efficient approach to generating highly contiguous and informative assemblies. However, higher relative error rates can introduce frameshifts and premature stop codons that pseudogenize genes, hindering downstream analyses. We developed a software tool that detects gene-fragmenting errors in draft assemblies of small genomes through comparison with a curated set of reference genome sequences and raw read information. In our presented example, detected errors represent less than 0.05% of the genome, but when corrected reduced the rate of pseudogenes from 23.3 to 5.6% in example long read assemblies, comparable to the rate of pseudogenes in short read assemblies. We demonstrate that this software can detect assembly errors in long read assemblies generated from small genomes and correct them to de-fragment genes.

Genomes of Prochlorococcus, Synechococcus, bacteria, and viruses recovered from marine picocyanobacteria cultures based on Illumina and Qitan nanopore sequencing

Prochlorococcus and Synechococcus are key contributors to marine primary production and play essential roles in global biogeochemical cycles. Despite the ecological importance of these two picocyanobacterial genera, current genomic datasets still lack comprehensive representation of under-sampled ocean regions, associated bacteria and viruses. To address this gap, we used a combination of second- and third-generation sequencing technologies to assemble comprehensive genomic data from 105 Picocyanobacterial enrichment cultures isolated from the Indian Ocean, the South China Sea, and the western Pacific Ocean. This dataset includes 55 Prochlorococcus and 50 Synechococcus genomes with high completeness (>98%) and low contamination (<2%), along with 308 non-redundant associated bacterial genomes derived from 1,457 medium- and high-quality non-cyanobacteria metagenome-assembled genomes (MAGs, completeness ≥50% and contamination ≤10%). Additionally, 2,113 non-redundant viral operational taxonomic units (vOTUs) were derived from a total of 7632 qualified viral contigs. This dataset provides a valuable resource for improving our understanding of the complex interactions among Prochlorococcus, Synechococcus, and their associated bacteria and viruses in marine ecosystems, offering a foundation to study their ecological roles and evolutionary dynamics.

Gut microbiome evolution from infancy to 8 years of age

The human gut microbiome is most dynamic in early life. Although sweeping changes in taxonomic architecture are well described, it remains unknown how, and to what extent, individual strains colonize and persist and how selective pressures define their genomic architecture. In this study, we combined shotgun sequencing of 1,203 stool samples from 26 mothers and their twins (52 infants), sampled from childbirth to 8 years after birth, with culture-enhanced, deep short-read and long-read stool sequencing from a subset of 10 twins (20 infants) to define transmission, persistence and evolutionary trajectories of gut species from infancy to middle childhood. We constructed 3,995 strain-resolved metagenome-assembled genomes across 399 taxa, and we found that 27.4% persist within individuals. We identified 726 strains shared within families, with Bacteroidales, Oscillospiraceae and Lachnospiraceae, but not Bifidobacteriaceae, vertically transferred. Lastly, we identified weaning as a critical inflection point that accelerates bacterial mutation rates and separates functional profiles of genes accruing mutations.

Three reference genomes for freshwater diatom ecology and evolution

Diatoms are an important component of marine and freshwater ecosystems. Although the majority of described diatom species live in freshwater systems, genome sequencing efforts have focused primarily on marine species. Genomic resources for freshwater species have the potential to improve our understanding of diatom ecology and evolution, particularly in the context of major environmental shifts. We used long- and short-read sequencing platforms to assemble reference genomes for three freshwater diatom species, all in the order Thalalassiosirales, which are abundant in the plankton of oceans, lakes, reservoirs, and rivers worldwide. We targeted three species that cover the breadth of phylogenetic diversity in the cyclostephanoid clade of Thalassiosirales: Cyclostephanos tholiformis (JALLPB020000000), Discostella pseudostelligera (JALLBG020000000), and Praestephanos triporus (JALLAZ020000000). The reference genome for D. pseudostelligera was considerably smaller (39 Mb) than those of both P. triporus (73 Mb) and C. tholiformis (177 Mb). Long-read sequencing allowed for the assembly of scaffold-level genomes, including regions rich in repetitive DNA. Compared to short-read assemblies, long-read assemblies increased the contig N50 length as much as 37-fold and reduced the number of contigs by more than 88%. Transcriptome-guided annotation of the protein-coding genes identified between 10,000 and 12,000 genes. This work provides further demonstration of the value of long-read sequencing and provides novel genomic resources for understanding the ecology and evolution of freshwater diatoms.

Classification and prediction of variants associated with hearing loss using sequence information in the vicinity of mutation sites

Hearing impairment is a major global health problem, affecting more than 5 % of the world's population at various ages, from neonates to the elderly. Among the common genetic variations in humans, single nucleotide variations and small insertions or deletions predominate. The study of hearing loss resulting from these variations is proving invaluable in the analysis and diagnosis of hearing disorders. The identification of pathogenic mutations is frequently a lengthy and laborious process. Existing computational prediction tools have been developed primarily for common diseases and genome-wide analyses, with less focus on deafness. This study proposes a novel approach that focuses on the regions surrounding mutation sites. Mutation sites associated with deafness and their flanking regions of different lengths were extracted from relevant databases and combined into seven distinct segments of different lengths. The information-theoretic features of these segments were computed. Five machine learning algorithms were then used for training, resulting in the construction of a model capable of classifying and predicting deafness-related mutations. For fragments encompassing the 250 bp regions upstream and downstream of the mutations, the average AUC of the five classifiers on the independent test set is 0.89 and the average ACC is 0.85, indicating that the model has a high recognition rate of the pathogenic deafness mutation site. An ensemble approach was also applied to predict variants of uncertain significance (VUS) that may be associated with deafness. These variants were then scored and ranked to assess their likelihood of contributing to the condition.

Whole-Genome Sequencing of Linezolid-Resistant and Linezolid-Intermediate-Susceptibility Enterococcus faecalis Clinical Isolates in a Mexican Tertiary Care University Hospital

Linezolid-non-susceptible Enterococcus faecalis (LNSEf) has emerged as a critical clinical concern worldwide, yet data from Latin American settings remain scarce. This study aimed to investigate the molecular epidemiology and mechanisms underlying LNSEf in a Mexican tertiary care university hospital, focusing on clinical correlates and clonal relationships. A total of 392 non-duplicated E. faecalis isolates were collected over 12 months, of which 24 with minimum inhibitory concentrations ≥4 µg/mL underwent whole-genome sequencing to identify specific resistance determinants (optrA, cfrA, 23S rRNA mutations) and to perform multilocus sequence typing (MLST) and phylogenetic analyses. Of the 392 isolates, 6.12% showed linezolid non-susceptibility, predominantly linked to plasmid- or chromosomally encoded optrA; only two isolates carried cfrA. No mutations were detected in 23S rRNA domain V or ribosomal proteins L3/L4. Clinically, LNSEf strains were associated with immunosuppression, previous surgical interventions, and prolonged hospital stays. Although most LNSEf isolates retained susceptibility to ampicillin, vancomycin, and daptomycin, they exhibited high rates of resistance to other antibiotic classes, particularly aminoglycosides and fluoroquinolones. These findings underscore the emergence of LNSEf in this region, highlighting the need for robust genomic surveillance, strict infection control, and judicious antimicrobial stewardship to curb further dissemination.

Tracing the Origin and Evolution of the Fungal Mycophenolic Acid Biosynthesis Pathway

Like bacteria and plants, fungi produce a remarkable diversity of small molecules with potent activities for human health known as natural products or secondary metabolites. One such example is mycophenolic acid, a powerful immunosuppressant drug that is administered daily to millions of transplant recipients worldwide. Production of mycophenolic acid is restricted to a very limited number of filamentous fungi, and little is known about its biosynthetic modalities. It is therefore a particular challenge to improve our knowledge of the biosynthesis of this valuable natural compound, as this would contribute to a better understanding of the specialized metabolism of fungi and could also lead to the identification of new fungal producers for the supply of immunosuppressants. Here, we were interested in deciphering the origin and evolution of the fungal mycophenolic acid biosynthetic pathway. Large-scale analyses of fungal genomic resources led us to identify several new species that harbor a gene cluster for mycophenolic acid biosynthesis. Phylogenomic analysis suggests that the mycophenolic acid biosynthetic gene cluster originated early in a common ancestor of the fungal family Aspergillaceae but was repeatedly lost and it is now present in a narrow but diverse set of filamentous fungi. Moreover, a comparison of the inosine 5'-monophosphate dehydrogenase protein sequences that are the target of the mycophenolic acid drug as well as analysis of mycophenolic acid production and susceptibility suggest that all mycophenolic acid fungal producers are resistant to this toxic compound, but that this resistance is likely to be based on different molecular mechanisms. Our study provides new insight into the evolution of the biosynthesis of the important secondary metabolite mycophenolic acid in fungi.

Genomic Insights and Antimicrobial Potential of Newly Streptomyces cavourensis Isolated from a Ramsar Wetland Ecosystem

The growing threat of antimicrobial resistance underscores the urgent need to identify new bioactive compounds. In this study, a Streptomyces strain, ACT158, was isolated from a Ramsar wetland ecosystem and found to exhibit broad-spectrum effects against Gram-positive and Gram-negative bacteria, as well as fungal pathogens. The active strain was characterized as S. cavourensis according to its morphology, phylogenetic analysis, average nucleotide identity (ANI), and digital DNA-DNA hybridization (dDDH). Whole-genome sequencing (WGS) and annotation revealed a genome size of 6.86 Mb with 5122 coding sequences linked to carbohydrate metabolism, secondary metabolite biosynthesis, and stress responses. Genome mining through antiSMASH revealed 32 biosynthetic gene clusters (BGCs), including those encoding polyketides, nonribosomal peptides, and terpenes, many of which showed low similarity to known clusters. Comparative genomic analysis, showing high genomic synteny with closely related strains. Unique genomic features of ACT158 included additional BGCs and distinct genes associated with biosynthesis pathways and stress adaptation. These findings highlight the strain's potential as a rich source of bioactive compounds and provide insights into its genomic basis for antimicrobial production and its ecological and biotechnological significance.

Whole-genome automated assembly pipeline for Chlamydia trachomatis strains from reference, in vitro and clinical samples using the integrated CtGAP pipeline

Whole genome sequencing (WGS) is pivotal for the molecular characterization of Chlamydia trachomatis (Ct)-the leading bacterial cause of sexually transmitted infections and infectious blindness worldwide. Ct WGS can inform epidemiologic, public health and outbreak investigations of these human-restricted pathogens. However, challenges persist in generating high-quality genomes for downstream analyses given its obligate intracellular nature and difficulty with in vitro propagation. No single tool exists for the entirety of Ct genome assembly, necessitating the adaptation of multiple programs with varying success. Compounding this issue is the absence of reliable Ct reference strain genomes. We, therefore, developed CtGAP-Chlamydia trachomatisGenome Assembly Pipeline-as an integrated 'one-stop-shop' pipeline for assembly and characterization of Ct genome sequencing data from various sources including isolates, in vitro samples, clinical swabs and urine. CtGAP, written in Snakemake, enables read quality statistics output, adapter and quality trimming, host read removal, de novo and reference-guided assembly, contig scaffolding, selective ompA, multi-locus-sequence and plasmid typing, phylogenetic tree construction, and recombinant genome identification. Twenty Ct reference genomes were also generated. Successfully validated on a diverse collection of 363 samples containing Ct, CtGAP represents a novel pipeline requiring minimal bioinformatics expertise with easy adaptation for use with other bacterial species.

A genomic and phenotypic investigation of pigeon-adaptive Salmonella

Salmonella, a significant threat to public safety, inflicts substantial economic losses on the poultry industry. The unique "parental feeding" breeding model of pigeon farms, against the "all-in & all-out" biosecurity strategy, makes them susceptible to Salmonella infections and subsequent outbreaks of pigeon paratyphoid. This study initially studied three pigeon paratyphoid outbreak incidents in Henan, China, in which 53 strains of pigeon-origin Salmonella Typhimurium (STM) were identified. Whole-genome sequencing (WGS) and antimicrobial-resistant profile analysis revealed that the three outbreaks were caused by distinct STM clones (ST128-DT2, ST19-DT99). Global phylogenetic analysis suggested that the United States is a possible origin, indicating a risk of intercontinental transmission via pigeon eggs. Further bacterial virulence and invasion assays, including in vitro and in vivo assays, revealed that pigeon-host-adaptive STM, compared to broad-host-range STM, carried fewer resistance genes, exhibited higher invasion indices and pseudogene levels, displayed a non-rdar (red dry and rough) phenotype, and had strong biofilm formation capability. Additionally, they showed reduced virulence and invasiveness in mice but a pigeon-adaptive feature in cogent models. The collective results support the host adaptation for pigeons among DT2 and DT99 phage-type isolates.

Optimized production of concanamycins using a rational metabolic engineering strategy

Plecomacrolides, such as concanamycins and bafilomycins, are potent and specific inhibitors of vacuolar-type ATPase. Concanamycins are 18-membered macrolides with promising therapeutic potential against multiple diseases, including viral infection, osteoporosis, and cancer. Due to the complexity of their total synthesis, the production of concanamycins is only achieved through microbial fermentation. However, the low titers of concanamycin A and its analogs in the native producing strains are a significant bottleneck for scale-up, robust structure-activity relationship studies, and drug development. To address this challenge, we designed a library of engineered Streptomyces strains for the overproduction of concanamycin A-C by combining the overexpression of target regulatory genes with the optimization of fermentation media. Integration of two endogenous regulators from the concanamycin biosynthetic gene cluster (cms) and one heterologous regulatory gene from the bafilomycin biosynthetic gene cluster significantly increased production of concanamycin A and its less abundant analog concanamycin B in Streptomyces eitanensis. The highest titers reported to date were observed in the engineered S. eitanensis DHS10676, which produced over 900 mg/L of concanamycin A and 300 mg/L of concanamycin B. Heterologous overexpression of the identified target regulatory genes across a panel of Streptomyces spp. harboring a putative concanamycin biosynthetic gene cluster confirmed its identity, and significantly improved concanamycin A production in all tested strains. Strain engineering, optimization of fermentation, and extraction purification protocols enabled swift access to these structurally complex plecomacrolides for semi-synthetic medicinal chemistry-based approaches. Together, this work established a platform for robust overproduction of concanamycin analogs across species.

Methodological aspects of investigating the resistome in pig farm environments

A typical One Health issue, antimicrobial resistance (AMR) development and its spread among people, animals, and the environment attracts significant research attention. The animal sector is one of the major contributors to the development and dissemination of AMR and accounts for more than 50 % of global antibiotics usage. The use of antibiotics exerts a selective pressure for resistant bacteria in the exposed microbiome, but many questions about the epidemiology of AMR in farm environments remain unanswered. This is connected to several methodological challenges and limitations, such as inconsistent sampling methods, complexity of farm environment samples and the lack of standardized protocols for sample collection, processing and bioinformatical analysis. In this project, we combined metagenomics and bioinformatics to optimise the methodology for reproducible research on the resistome in complex samples from the indoor farm environment. The work included optimizing sample collection, transportation, and storage, as well as DNA extraction, sequencing, and bioinformatic analysis, such as metagenome assembly and antibiotic resistance gene (ARG) detection. Our studies suggest that the current most optimal and cost-effective pipeline for ARG search should be based on Illumina sequencing of sock sample material at high depth (at least 25 M 250 bp PE for AMR gene families and 43 M for gene variants). We present a computational analysis utilizing MEGAHIT assembly to balance the identification of bacteria carrying ARGs with the potential loss of diversity and abundance of resistance genes. Our findings indicate that searching against multiple ARG databases is essential for detecting the highest diversity of ARGs.

Emerging Mobile Colistin Resistance Gene Mcr-1 and Mcr-10 in Enterobacteriaceae Isolates From Urban Sewage in China

Purpose

This study aimed to investigate the epidemiology and dissemination of mcr-positive Enterobacteriaceae in urban sewage in Yangzhou, China.

Methods

A total of 366 sewage samples were collected from the Yangzhou Wastewater Treatment Plant in Jiangsu Province. Colistin-resistant Enterobacteriaceae was identified through PCR targeting mcr-1 to mcr-10 genes. The isolates underwent antimicrobial susceptibility testing, and whole-genome sequencing was performed to analyze their genomic features. Additionally, conjugation experiments were conducted to assess the transferability of mcr-positive plasmids.

Results

Three mcr-positive Enterobacteriaceae isolates were identified, representing an isolation rate of 0.82%. These included one mcr-1-positive Escherichia coli (ST167) and two mcr-10-positive Klebsiella pneumoniae complex strains with novel sequence types ST6801 and ST6825. The mcr-1 gene was located on an IncI2 plasmid (pYZ22WS208_3) and successfully transferred to recipient strains. In contrast, the mcr-10 gene was carried on IncF plasmids (pYZ22WS067_1 and pYZ22WS223_1) but was not transferable in this study. Phylogenetic analysis revealed that the mcr-1-positive E. coli strain clustered within Clade II, alongside strains from various countries and sources. Phylogenomic analysis of mcr-10-positive isolates showed their sporadic distribution across 13 countries, with associations to diverse hosts and environments, indicating potential for widespread transmission.

Conclusion

This study demonstrates the presence of mcr-1 and mcr-10-positive Enterobacteriaceae in wastewater, emphasizing the importance of wastewater surveillance for tracking antibiotic resistance. The horizontal transfer of mcr-1 and potential spread of mcr-10 across various hosts underscore the need for ongoing monitoring and preventive measures.

An Inversion Polymorphism Under Balancing Selection, Involving Giant Mobile Elements, in an Invasive Fungal Pathogen

Recombination suppression can evolve in sex or mating-type chromosomes, or in autosomal supergenes, with different haplotypes being maintained by balancing selection. In the invasive chestnut blight fungus Cryphonectria parasitica, a genomic region was suggested to lack recombination and to be partially physically linked to the mating-type (MAT) locus based on segregation analyses. Using hundreds of available C. parasitica genomes and generating new high-quality genome assemblies, we show that a ca. 1.2 Mb genomic region proximal to the mating-type locus lacks recombination, with the segregation of two highly differentiated haplotypes in balanced proportions in invasive populations. High-quality genome assemblies further revealed an inversion in one of the haplotypes in the invaded range. The two haplotypes were estimated to have diverged 1.5 million years ago, and each harboured specific genes, some of which likely belonging to Starships. These are large transposable elements, mobilized by tyrosine recombinases, able to move accessory genes, and involved in adaptation in multiple fungi. The MAT-proximal region carried genes upregulated under virus infection or vegetative incompatibility reaction. In the native range, the MAT-proximal region also appeared to have a different evolutionary history than the rest of the genome. In all continents, the MAT-Proximal region was enriched in nonsynonymous substitutions, in gene presence/absence polymorphism, in tyrosine recombinases and in transposable elements. This study thus sheds light on a case of a large nonrecombining region partially linked to a mating compatibility locus, likely maintained by balancing selection on differentiated haplotypes, possibly involved in adaptation in a devastating tree pathogen.

Genomic, Transcriptomic and Suspect/Non-Target Screening Analyses Reveal the Role of CYP450s in the Degradation of Imazalil and Delineate Its Transformation Pathway by Cladosporium herbarum

Imazalil (IMZ), a major surface water contaminant characterised by high environmental recalcitrance and toxicity, is used in fruit-packaging plants to control fungal infestations during storage. This leads to the production of wastewaters which should be treated on site before their environmental release. We previously isolated a Cladosporium herbarum strain, the first microorganism that could degrade IMZ. Here we describe the genetic network utilised by the fungus to degrade IMZ and its detailed transformation. Genomic and transcriptomic analysis of C. herbarum pointed to the involvement of strongly upregulated CYP450s in IMZ degradation, as further verified by cessation of its biodegradation by CYP450 inhibitors. LC-QTOF-HRMS analysis and suspect/non-target screening identified nine transformation products (TPs) of IMZ. IMZ biotransformation mainly proceeded through O-dealkylation, while other less important paths, most probably controlled by upregulated oxidases, were operative involving successive hydroxylation reactions. These lead to the formation of TPs like IMZ_313 and IMZ_331, with the former being further transformed through imidazole ring scission to IMZ_288, a TP reported for the first time. We provide first evidence for the transformation mechanism of IMZ by C. herbarum and the potential genes/enzymes involved, paving the way for the use of C. herbarum in the biodepuration of agro-industrial effluents.

Whole genome sequence data of Trichoderma yunnanense strain TM10, a plant growth-promoting fungus and biocontrol agent

Trichoderma yunnanense strain TM10 was isolated from rhizosphere soil of rice plants cultivated under system of rice intensification (SRI) practises in West Java, Indonesia. It exhibits significant potential as a plant growth promoter and biocontrol agent in rice plants. Although this strain has shown promise in promoting plant growth and suppressing phytopathogens under in vitro and in planta conditions, there is still a lack of genomic data to elucidate the molecular mechanisms underlying its plant growth-promoting and biocontrol capabilities. This study reports the whole genome sequence of T. yunnanense strain TM10. The genome of the fungus was sequenced using the MGI DNBSEQ-G400 high-throughput sequencing platform. The assembled genome of T. yunnanense strain TM10 was approximately 36 Mbp in length, comprising 385 contigs with a GC content of 48 % and a sequencing coverage of 43.8×. This genomic data provides a foundation for harnessing the plant growth-promoting and biocontrol potential of this strain. The complete genome sequence has been deposited at the National Center for Biotechnology Information (NCBI) under Bioproject accession number PRJNA1181959, BioSample ID SAMN44575400, and genome accession number JBIYZQ000000000. These data are valuable for further research into the biotechnological potential of this strain and for exploring the molecular mechanisms underlying its plant growth-promoting and biocontrol activities.

Determining the incidence, risk factors and biological drivers of irritable bowel syndrome (IBS) as part of the constellation of postacute sequelae of SARS-CoV-2 infection (PASC) outcomes in the Arizona CoVHORT-GI: a longitudinal cohort study

INTRODUCTION

Postacute sequelae of SARS-CoV-2 infection (PASC) are extensive. Also known as long COVID, primary outcomes reported are neurologic, cardiac and respiratory in nature. However, several studies have also reported an increase in gastrointestinal (GI) symptoms and syndromes following COVID-19. This study of PASC will include extensive analyses of GI symptoms, determine if people with pre-existing irritable bowel syndrome (IBS) are at higher risk of developing PASC generally or PASC-GI, and which biomarkers are impacted and to what degree. This R01 study is being funded by the National Institute of Diabetes and Digestive and Kidney Diseases (1R01DK135483-01) from 2023 to 2028.

METHODS AND ANALYSES

This study combines a longitudinal epidemiologic cohort study and in-depth, novel biologic analyses. In collaboration with a pre-existing study, the Arizona CoVID-19 Cohort (CoVHORT)-GI will recruit participants based on the history of COVID infection(s), new or ongoing GI symptoms 3-6 months postinfection, and pre-existing or incident IBS diagnosis to represent five study groups for comparison and analyses. A subset (n=1000) of those recruited will submit both stool and blood samples. Both samples will undergo a novel method to quantitate humoral and mucosal immune responses to host-derived faecal communities in conjunction with magnetic bead-based separation and high-depth shotgun microbial sequencing. Stool samples will also undergo traditional microbiome analyses (diversity and abundance) and faecal calprotectin assays. Additional serum analyses will aim to determine if a proteomics-based signature exists that differentiates a unique biomarker compositional signature discriminating PASC-GI versus no PASC. All laboratory data will be linked with in-depth epidemiologic data on demographics, symptoms and chronic conditions.

ETHICS AND DISSEMINATION

This study involves human participants and was approved by the University of Arizona Institutional Review Board (IRB (#00002332) and has been deemed minimal risk. Participants gave informed consent to participate in the study before taking part. All publications from the study will be shared back to participants along with alternative lay summaries and webinars to communicate key findings. The data management plan has been published and is publicly available online, including protocols for data requests.

The plastomes of Cypripedium (Orchidaceae: Cypripedioideae) exhibit atypical GC content and genome size based on different sequencing strategies

Most of the sequenced plastomes of photosynthetic angiosperms exhibit conservation in size, gene content, gene order, and GC content. In contrast, the sequenced plastomes of Cypripedium are distinguished by genome size expansion, AT-biased base composition, structural variation, and a low substitution rate. Additionally, the impact of sequencing methods is seldom addressed in prior studies, and the species represented in these studies are underrepresented. These atypical plastome features render the genus an ideal candidate for investigating plastome evolution. Besides, the backbone relationships within the genus remain poorly resolved. In this study, we sequenced twelve Cypripedium plastomes using three distinct sequencing strategies and obtained an additional 27 sequences from GenBank for comparative analysis. We classified the plastomes of the genus into two types: one resembling those of most other angiosperms, and the other characterized by inverted repeat (IR) expansion and small single copy (SSC) contraction. The plastomes within this genus exhibit significant size variations (∼72 kb), variations in GC content, and structural differences at the genus level. Furthermore, our comparative analysis revealed that the choice of sequencing strategy significantly impacts the assembly results. The uncovered regions in samples sequenced with short-read technology are predominantly AT-rich, suggesting that short-read sequencing may lead to assembly errors in plastomes with AT-rich regions and long repeats. Additionally, we have reconstructed the phylogeny of the genus using plastome-level data. However, the phylogenetic relationships within the genus remain partially solved. This study provides new insights into the evolution of plastomes, particularly those with AT-rich base compositions and genomes containing long repeat regions.

Rumen DNA virome and its relationship with feed efficiency in dairy cows

BACKGROUND

The rumen harbors a diverse virome that interacts with other microorganisms, playing pivotal roles in modulating metabolic processes within the rumen environment. However, the characterization of rumen viruses remains incomplete, and their association with production traits, such as feed efficiency (FE), has not been documented. In this study, rumen fluid from 30 Chinese Holstein dairy cows was analyzed using next-generation sequencing (NGS) and High-Fidelity (HiFi) sequencing to elucidate the rumen DNA virome profile and uncover potential viral mechanisms influencing FE.

RESULTS

Integrated NGS and HiFi sequencing enhanced the length, completeness, and resolution of viral operational taxonomic units (vOTUs) compared to NGS. A total of 6,922 vOTUs were identified, including 4,716 lytic and 1,961 temperate vOTUs. At the family level, lytic viruses were predominantly from Siphoviridae (30.35%) and Schitoviridae (23.93%), while temperate viruses were primarily Siphoviridae (67.21%). The study annotated 2,382 auxiliary metabolic genes (AMGs), comprising 1,752 lytic virus-associated AMGs across 51 functional categories and 589 temperate virus-associated AMGs across 29 categories. Additionally, 2,232 vOTU-host metagenome-assembled genome (hMAG) linkages were predicted, with Firmicutes_A (33.60%) and Bacteroidota (33.24%) being the most prevalent host phyla. Significant differences in viral populations were observed between high and low FE groups across multiple taxonomic levels (P < 0.05). Two pathways were proposed to explain how rumen viruses might modulate FE: (1) Lytic viruses could lyse beneficial host bacteria linked to favorable cattle phenotypes, such as vOTU1836 targeting Ruminococcaceae, resulting in diminished organic acid production and consequently lower FE; (2) AMG-mediated host metabolism modulation, exemplified by GT2 carried by vOTU0897, which may enhance Lachnospiraceae fermentation capacity, increasing organic acid production and thereby improving FE.

CONCLUSIONS

This study constructed a comprehensive rumen DNA virome profile for Holstein dairy cows, elucidating the structural and functional complexity of rumen viruses, the roles of AMGs, and vOTU-hMAG linkages. The integration of these data offers novel insights into the mechanisms by which rumen viruses may regulate nutrient utilization, potentially influencing FE in dairy cows. Video Abstract.

Navigating the Complex Terrain of Methane Synthesis: Multienzyme Control Points and Data-Driven Strategies

Anaerobic digestion is a crucial process in wastewater treatment, renowned for its sustainable biogas production capabilities and the simultaneous reduction of environmental pollution. However, dysregulation of vital biological processes and pathways can lead to reduced efficiency and suboptimal biogas output, which can be seen through low counts per million of sequences related to three critical control points for methane synthesis. Namely, tetrahydromethanopterin S-methyltransferase (MTR), methyl-coenzyme reductase M (MCR), and CoB/CoM heterodisulfide oxidoreductase (HDR) are the last reactions that must occur. This study leveraged sequencing data from NCBI's Sequence Read Archive (SRA) database whose sample origins were of two model full-scale biodigesters. By assembling a genome BBBAS3_2 (86.4% complete), this study was able to align, assemble, and follow expression patterns related to KEGG pathways and sample conditions. This study detected and estimated expression patterns of conserved (in methanogenic archae) alleles for electron cycling by the heterodisulfide reductase complex, methylation and demethylation by methyltransferases, and oxidation of active coenzymes A and B by methyl-coenzyme reductase to surveil counts of sequences at critical control points. Utilizing a streamlined cloud bioinformatics approach with Omicsbox and Kbase, this study bridges the division between data-intensive sequencing and innovative solutions for methane optimization.

Complementary insights into gut viral genomes: a comparative benchmark of short- and long-read metagenomes using diverse assemblers and binners

BACKGROUND

Metagenome-assembled viral genomes have significantly advanced the discovery and characterization of the human gut virome. However, we lack a comparative assessment of assembly tools on the efficacy of viral genome identification, particularly across next-generation sequencing (NGS) and third-generation sequencing (TGS) data.

RESULTS

We evaluated the efficiency of NGS, TGS, and hybrid assemblers for viral genome discovery using 95 viral-like particle (VLP)-enriched fecal samples sequenced on both Illumina and PacBio platforms. MEGAHIT, metaFlye, and hybridSPAdes emerged as the optimal choices for NGS, TGS, and hybrid datasets, respectively. Notably, these assemblers recovered distinct viral genomes, demonstrating a remarkable degree of complementarity. By combining individual assembler results, we expanded the total number of nonredundant high-quality viral genomes by 4.83 ~ 21.7-fold compared to individual assemblers. Among them, viral genomes from NGS and TGS data have the least overlap, indicating the impact of data type on viral genome recovery. We also evaluated four binning methods, finding that CONCOCT incorporated more unrelated contigs into the same bins, while MetaBAT2, AVAMB, and vRhyme balanced inclusiveness and taxonomic consistency within bins.

CONCLUSIONS

Our findings highlight the challenges in metagenome-driven viral discovery, underscoring tool limitations. We advocate for combined use of multiple assemblers and sequencing technologies when feasible and highlight the urgent need for specialized tools tailored to gut virome assembly. This study contributes essential insights for advancing viral genome research in the context of gut metagenomics. Video Abstract.

Sequencing-based analysis of microbiomes

Microbiomes occupy a range of niches and, in addition to having diverse compositions, they have varied functional roles that have an impact on agriculture, environmental sciences, and human health and disease. The study of microbiomes has been facilitated by recent technological and analytical advances, such as cheaper and higher-throughput DNA and RNA sequencing, improved long-read sequencing and innovative computational analysis methods. These advances are providing a deeper understanding of microbiomes at the genomic, transcriptional and translational level, generating insights into their function and composition at resolutions beyond the species level.

Comparative genomics analysis of Salmonella Enteritidis isolated from clinical cases associated with chicken

Salmonella Enteritidis is a major foodborne pathogen, and the emergence of multidrug-resistant (MDR) S. Enteritidis poses a serious public health challenge. In this study, we report the genomic characterization of five S. Enteritidis isolates from clinical. These isolates exhibited resistance to seven classes of antimicrobials with four of the five characterized as MDR. Isolate 33 A exhibited resistance to colistin and polymyxin B, while no associated antimicrobial resistance genes (ARGs) were identified in its genome. Isolate 21 A and 44 A were extended-spectrum beta-lactamases-producing (ESBLs). Whole genome sequencing analysis revealed the presence of multiple mobile genetic elements (MGEs), including plasmids, prophages, and genomic islands, which may have facilitated the acquisition and dissemination of ARGs. Notably, several ARGs, including blaCTX-M-55, blaTEM-141, blaTEM-1B, aph(3')-IIa, aph(3'')-Ib, aph(6)-Id, tet(A), floR, fosA3, and sul2, were identified on plasmids. In addition, chromosomal point mutations in gyrA (D87G and D87Y) and acrB (F28L and L40P) were also observed in each isolate. Multiple virulence genes associated with the type III secretion system were identified on Salmonella pathogenicity islands (SPIs) SPI-1 and SPI-2. Phylogenetic analysis revealed that the five isolates, along with a clinical and chicken origin isolates in the database, clustered together, suggesting a probable common source of infection. Our findings highlight the intricate genetic mechanisms behind MDR in S. Enteritidis, emphasizing the ongoing necessity for surveillance and appropriate antimicrobial usage. This contributes to our understanding of S. Enteritidis transmission within the food chain.

Metagenome-Assembled Genomes of Pig Fecal Samples in Nine European Countries: Insights into Antibiotic Resistance Genes and Viruses

Gut microbiota plays a crucial role in the health and productivity of pigs. However, the spread of antibiotic resistance genes (ARGs) and viruses within the pig intestinal microbiota poses significant threats to animal and public health. This study utilized 181 pig samples from nine European countries and employed metagenomic assembly methods to investigate the dynamics and distribution of ARGs and viruses within the pig intestinal microbiota, aiming to observing their associations with potential bacterial hosts. We identified 4605 metagenome-assembled genomes (MAGs), corresponding to 19 bacterial phyla, 97 families, 309 genera, and a total of 449 species. Additionally, 44 MAGs were classified as archaea. Analysis of ARGs revealed 276 ARG types across 21 ARG classes, with Glycopeptide being the most abundant ARG class, followed by the class of Multidrug. Treponema D sp016293915 was identified as a primary potential bacterial host for Glycopeptide. Aligning nucleotide sequences with a viral database, we identified 1044 viruses. Among the viral genome families, Peduoviridae and Intestiviridae were the most prevalent, with CAG-914 sp000437895 being the most common potential host species for both. These findings highlight the importance of MAGs in enhancing our understanding of the gut microbiome, revealing microbial diversity, antibiotic resistance, and virus-bacteria interactions. The data analysis for the article was based on the public dataset PRJEB22062 in the European Nucleotide Archive.

Unlocking the growth-promoting and antagonistic power: A comprehensive whole genome study on Bacillus velezensis strains

Bacillus species are extensively documented as plant growth-promoting rhizobacteria, contributing significantly to the enhancement of soil fertility, nutrient recycling, and the control of phytopathogens. Utilizing them as biocontrol agents represents an environmentally friendly strategy, particularly within the rhizospheric community. This study presents the comprehensive genome sequences of three B. velezensis strains (LGMB12, LGMB319, and LGMB426) which were previously isolated from root samples of maize (Zea mays L.), along with a type strain FZB42. The research assesses the capability of the three strains for antagonizing fungi, specifically Fusarium graminearum, Fusarium verticillioides, Colletotrichum graminicola, and Stenocarpella sp. In paired cultures involving maize fungi, treatments containing bacteria B. velezensis exhibited statistically significant differences compared to both negative and positive treatments in terms of antagonism. Furthermore, genome mining techniques were employed to explore their inherent antagonistic potential. The assembly revealed that strains LGMB12, LGMB319, LGMB426, and FZB42 exhibit genome sizes of 4,187,541 bp, 4,244,954 bp, 3,976,537 bp, and 3,990,518 respectively. Their respective G + C content stands at 46.42 %, 46.50 %, 46.51 %, and 46.38 %. Moreover, the genomes present multiple gene clusters responsible for the synthesis of secondary metabolites and carbohydrate-active enzymes (CAZymes). These clusters highlight a diverse array of antibacterial and antifungal properties, complemented by numerous plant growth-promoting genes. These results highlight the potential of B. velezensis LGMB12, LGMB319, and LGMB426 strains as biocontrol and plant growth promotion agents, being promising candidates for further studies in agricultural production, including field trials.

Two distinct Dehalobacter metagenome-assembled genomes from anaerobic chloroform and dichloromethane degrading consortia

Here we present two metagenomes and two Dehalobacter metagenome-assembled genomes from subcultures of an anaerobic chloroform and dichloromethane degrading microbial community used for bioremediation. Our objective was to assemble and curate the genome(s) of Dehalobacter, key biodegraders in the culture, through repeated sequencing and joint assembly with previous datasets.

Simple, reference-independent assessment to empirically guide correction and polishing of hybrid microbial community metagenomic assembly

Hybrid metagenomic assembly of microbial communities, leveraging both long- and short-read sequencing technologies, is becoming an increasingly accessible approach, yet its widespread application faces several challenges. High-quality references may not be available for assembly accuracy comparisons common for benchmarking, and certain aspects of hybrid assembly may benefit from dataset-dependent, empiric guidance rather than the application of a uniform approach. In this study, several simple, reference-free characteristics-particularly coding gene content and read recruitment profiles-were hypothesized to be reliable indicators of assembly quality improvement during iterative error-fixing processes. These characteristics were compared to reference-dependent genome- and gene-centric analyses common for microbial community metagenomic studies. Two laboratory-scale bioreactors were sequenced with short- and long-read platforms, and assembled with commonly used software packages. Following long read assembly, long read correction and short read polishing were iterated up to ten times to resolve errors. These iterative processes were shown to have a substantial effect on gene- and genome-centric community compositions. Simple, reference-free assembly characteristics, specifically changes in gene fragmentation and short read recruitment, were robustly correlated with advanced analyses common in published comparative studies, and therefore are suitable proxies for hybrid metagenome assembly quality to simplify the identification of the optimal number of correction and polishing iterations. As hybrid metagenomic sequencing approaches will likely remain relevant due to the low added cost of short-read sequencing for differential coverage binning or the ability to access lower abundance community members, it is imperative that users are equipped to estimate assembly quality prior to downstream analyses.

Pangenomic insights into Dehalobacter evolution and acquisition of functional genes for bioremediation

Dehalobacter is a genus of organohalide-respiring bacteria that is recognized for its fastidious growth using reductive dehalogenases (RDases). In the SC05 culture, however, a Dehalobacter population also mineralizes dichloromethane (DCM) produced by chloroform dechlorination using the mec cassette, just downstream of its active RDase. A closed genome of this DCM-mineralizing lineage has previously evaded assembly. Here, we present the genomes of two novel Dehalobacter strains, each of which was assembled from the metagenome of a distinct subculture from SC05. A pangenomic analysis of the Dehalobacter genus, including RDase synteny and phylogenomics, reveals at least five species of Dehalobacter based on average nucleotide identity, RDase and core gene synteny, as well as differential functional genes. An integration hotspot is also pinpointed in the Dehalobacter genome, in which many recombinase islands have accumulated. This nested recombinase island encodes the active RDase and mec cassette in both SC05 Dehalobacter genomes, indicating the transfer of key functional genes between species of Dehalobacter. Horizontal gene transfer between these two novel Dehalobacter strains has implications for the evolutionary history within the SC05 subcultures and of the Dehalobacter genus as a whole, especially regarding adaptation to anthropogenic chemicals.

Haplotype-aware sequence alignment to pangenome graphs

Modern pangenome graphs are built using haplotype-resolved genome assemblies. When mapping reads to a pangenome graph, prioritizing alignments that are consistent with the known haplotypes improves genotyping accuracy. However, the existing rigorous formulations for colinear chaining and alignment problems do not consider the haplotype paths in a pangenome graph. This often leads to spurious read alignments to those paths that are unlikely recombinations of the known haplotypes. In this paper, we develop novel formulations and algorithms for sequence-to-graph alignment and chaining problems. Inspired by the genotype imputation models, we assume that a query sequence is an imperfect mosaic of reference haplotypes. Accordingly, we introduce a recombination penalty in the scoring functions for each haplotype switch. First, we solve haplotype-aware sequence-to-graph alignment in [Formula: see text] time, where Q is the query sequence, E is the set of edges, and H is the set of haplotypes represented in the graph. To complement our solution, we prove that an algorithm significantly faster than [Formula: see text] is impossible under the strong exponential time hypothesis (SETH). Second, we propose a haplotype-aware chaining algorithm that runs in [Formula: see text] time after graph preprocessing, where N is the count of input anchors. We then establish that a chaining algorithm significantly faster than [Formula: see text] is impossible under SETH. As a proof-of-concept, we implemented our chaining algorithm in the Minichain aligner. By aligning sequences sampled from the human major histocompatibility complex (MHC) to a pangenome graph of 60 MHC haplotypes, we demonstrate that our algorithm achieves better consistency with ground-truth recombinations compared with a haplotype-agnostic algorithm.

Metagenome from ACT-3/CF: an anaerobic chloroform-degrading microbial community

Here, we present the metagenome of an anaerobic chloroform respiring mixed microbial community used for bioremediation. Our objective was to obtain draft genomes of key microorganisms in the culture.

HyLight: Strain aware assembly of low coverage metagenomes

Different strains of identical species can vary substantially in terms of their spectrum of biomedically relevant phenotypes. Reconstructing the genomes of microbial communities at the level of their strains poses significant challenges, because sequencing errors can obscure strain-specific variants. Next-generation sequencing (NGS) reads are too short to resolve complex genomic regions. Third-generation sequencing (TGS) reads, although longer, are prone to higher error rates or substantially more expensive. Limiting TGS coverage to reduce costs compromises the accuracy of the assemblies. This explains why prior approaches agree on losses in strain awareness, accuracy, tendentially excessive costs, or combinations thereof. We introduce HyLight, a metagenome assembly approach that addresses these challenges by implementing the complementary strengths of TGS and NGS data. HyLight employs strain-resolved overlap graphs (OG) to accurately reconstruct individual strains within microbial communities. Our experiments demonstrate that HyLight produces strain-aware and contiguous assemblies at minimal error content, while significantly reducing costs because utilizing low-coverage TGS data. HyLight achieves an average improvement of 19.05% in preserving strain identity and demonstrates near-complete strain awareness across diverse datasets. In summary, HyLight offers considerable advances in metagenome assembly, insofar as it delivers significantly enhanced strain awareness, contiguity, and accuracy without the typical compromises observed in existing approaches.

Comparative genomics unravels a rich set of biosynthetic gene clusters with distinct evolutionary trajectories across fungal species (Termitomyces) farmed by termites

The use of compounds produced by hosts or symbionts for defence against antagonists has been identified in many organisms, including in fungus-farming termites (Macrotermitinae). The obligate mutualistic fungus Termitomyces plays a pivotal role in plant biomass decomposition and as the primary food source for these termites. Despite the isolation of various specialized metabolites from different Termitomyces species, our grasp of their natural product repertoire remains incomplete. To address this knowledge gap, we conducted a comprehensive analysis of 39 Termitomyces genomes, representing 21 species associated with members of five termite host genera. We identified 754 biosynthetic gene clusters (BGCs) coding for specialized metabolites and categorized 660 BGCs into 61 biosynthetic gene cluster families (GCFs) spanning five compound classes. Seven GCFs were shared by all 21 Termitomyces species and 21 GCFs were present in all genomes of subsets of species. Evolutionary constraint analyses on the 25 most abundant GCFs revealed distinctive evolutionary histories, signifying that millions of years of termite-fungus symbiosis have influenced diverse biosynthetic pathways. This study unveils a wealth of non-random and largely undiscovered chemical potential within Termitomyces and contributes to our understanding of the intricate evolutionary trajectories of biosynthetic gene clusters in the context of long-standing symbiosis.

Comparative Analysis of the Mitochondrial Genome of Eggplant (Solanum melongena L.) to Identify Cytoplasmic Male Sterility Candidate Genes

Cytoplasmic male sterility (CMS) is important for commercial hybrid seed production. However, it is still not used in eggplant (Solanum melongena L.), and corresponding regulatory genes and mechanisms of action have not been reported. We report CMS line 327A, which was derived from the hybridization between cultivated and wild eggplants. By looking at different stages of anther development under a microscope, we saw that the 327A anther's tapetum layer vacuolized during meiosis, which caused abortion. To investigate the 327A CMS regulatory genes, the mitochondrial genomes of 327A and its maintainer line 327B were assembled de novo. It was found that 15 unique ORFs (Open Reading Frame) were identified in 327A. RT-PCR and RT-QPCAR tests confirmed that orf312a and orf172a, 327A-specific ORFs with a transmembrane domain, were strongly expressed in sterile anthers of 327A. In addition, orf312a has a chimeric structure with the ribosomal protein subunit rpl16. Therefore, orf312a and orf172a can be considered strong candidate genes for CMS. Concurrently, we analyzed the characteristics of CMS to develop a functional molecular marker, CMS312, targeting a future theoretical basis for eggplant CMS three-line molecular breeding.

Characterizing a stable five-species microbial community for use in experimental evolution and ecology

Model microbial communities are regularly used to test ecological and evolutionary theory as they are easy to manipulate and have fast generation times, allowing for large-scale, high-throughput experiments. A key assumption for most model microbial communities is that they stably coexist, but this is rarely tested experimentally. Here we report the (dis)assembly of a five-species microbial community from a metacommunity of soil microbes that can be used for future experiments. Using reciprocal invasion-from-rare experiments we show that all species can coexist and we demonstrate that the community is stable for a long time (~600 generations). Crucially for future work, we show that each species can be identified by their plate morphologies, even after >1 year in co-culture. We characterise pairwise species interactions and produce high-quality reference genomes for each species. This stable five-species community can be used to test key questions in microbial ecology and evolution.

Occurrence and characteristics of blaOXA-181-carrying Klebsiella aerogenes from swine in China

OBJECTIVES

Klebsiella aerogenes is a largely understudied opportunistic pathogen that can cause sepsis and lead to high mortality rates. In this study, we reported the occurrence of carbapenem-resistant blaOXA-181-carrying Klebsiella aerogenes from swine in China and elucidate their genomic characteristics.

METHODS

A total of 126 samples, including 109 swine fecal swabs, 14 environmental samples, and three feed samples were collected from a pig farm in China. The samples were enriched with LB broth culture and then inoculated into MacConkey agar plates for bacterial isolation. After PCR detection of carbapenemases genes, the blaOXA-181-carrying isolates were subjected to antimicrobial susceptibility testing, and whole-genome sequence analysis.

RESULTS

Four Klebsiella aerogenes isolates carrying the blaOXA-181 gene were obtained from swine faecal samples. All the 4 strains were belonged to ST438. The blaOXA-181 genes were located in IncX3-ColKP3 hybrid plasmids with the core genetic structure of IS26-ΔIS3000-ΔISEcp1-blaOXA-181-ΔlysR-ΔereA-ΔrepA-ISKpn19-tinR-qnrS1-ΔIS2-IS26, which suggests the potential for horizontal transfer and further dissemination of this resistance gene among Enterobacteriaceae and other sources.

CONCLUSIONS

This study represents the first instance of OXA-181-producing K. aerogenes being identified from swine faeces in China. It is crucial to maintain continuous monitoring and ongoing attention to the detection of K. aerogenes carrying blaOXA-181 and other resistance genes in pigs.

Trans2express - de novo transcriptome assembly pipeline optimized for gene expression analysis

BACKGROUND

As genomes of many eukaryotic species, especially plants, are large and complex, their de novo sequencing and assembly is still a difficult task despite progress in sequencing technologies. An alternative to genome assembly is the assembly of transcriptome, the set of RNA products of the expressed genes. While a bunch of de novo transcriptome assemblers exists, the challenges of transcriptomes (the existence of isoforms, the uneven expression levels across genes) complicates the generation of high-quality assemblies suitable for downstream analyses.

RESULTS

We developed Trans2express - a web-based tool and a pipeline of de novo hybrid transcriptome assembly and postprocessing based on rnaSPAdes with a set of subsequent filtrations. The pipeline was tested on Arabidopsis thaliana cDNA sequencing data obtained using Illumina and Oxford Nanopore Technologies platforms and three non-model plant species. The comparison of structural characteristics of the transcriptome assembly with reference Arabidopsis genome revealed the high quality of assembled transcriptome with 86.1% of Arabidopsis expressed genes assembled as a single contig. We tested the applicability of the transcriptome assembly for gene expression analysis. For both Arabidopsis and non-model species the results showed high congruence of gene expression levels and sets of differentially expressed genes between analyses based on genome and based on the transcriptome assembly.

CONCLUSIONS

We present Trans2express - a protocol for de novo hybrid transcriptome assembly aimed at recovering of a single transcript per gene. We expect this protocol to promote the characterization of transcriptomes and gene expression analysis in non-model plants and web-based tool to be of use to a wide range of plant biologists.

Trojan Horse virus delivering CRISPR-AsCas12f1 controls plant bacterial wilt caused by Ralstonia solanacearum

Plant bacterial wilt caused by Ralstonia solanacearum results in huge losses. Accordingly, developing an effective control method for this disease is urgently required. Filamentous phages, which do not lyse host bacteria and exert minimal burden, offer a potential biocontrol solution. A filamentous phage RSCq that infects R. solanacearum was isolated in this study through genome mining. We constructed engineered filamentous phages based on RSCq by employing our proposed approach with wide applicability to non-model phages, enabling the exogenous genes delivery into bacterial cells. CRISPR-AsCas12f1 is a miniature class 2 type V-F CRISPR-Cas system. A CRISPR-AsCas12f1-based gene editing system that targets the key virulence regulator gene hrpB was developed, generating the engineered phage RSCqCRISPR-Cas. Similar to the Greek soldiers in the Trojan Horse, our findings demonstrated that the engineered phage-delivered CRISPR-Cas system could disarm the key "weapon," hrpB, of R. solanacearum, in medium and plants. Remarkably, pretreatment with RSCqCRISPR-Cas significantly controlled tobacco bacterial wilt, highlighting the potential of engineered filamentous phages as promising biocontrol agents against plant bacterial diseases.IMPORTANCEBacterial disease, one of the major plant diseases, causes huge food and economic losses. Phage therapy, an environmentally friendly control strategy, has been frequently reported in plant bacterial disease control. However, host specificity, sensitivity to ultraviolet light and certain conditions, and bacterial resistance to phage impede the widespread application of phage therapy in crop production. Filamentous phages, which do not lyse host bacteria and exert minimal burden, offer a potential solution to overcome the limitations of lytic phage biocontrol. This study developed a genetic engineering approach with wide applicability to non-model filamentous phages and proved the application possibility of engineered phage-based gene delivery in plant bacterial disease biocontrol for the first.

Making a Pathogen? Evaluating the Impact of Protist Predation on the Evolution of Virulence in Serratia marcescens

Opportunistic pathogens are environmental microbes that are generally harmless and only occasionally cause disease. Unlike obligate pathogens, the growth and survival of opportunistic pathogens do not rely on host infection or transmission. Their versatile lifestyles make it challenging to decipher how and why virulence has evolved in opportunistic pathogens. The coincidental evolution hypothesis postulates that virulence results from exaptation or pleiotropy, i.e. traits evolved for adaptation to living in one environment that have a different function in another. In particular, adaptation to avoid or survive protist predation has been suggested to contribute to the evolution of bacterial virulence (the training ground hypothesis). Here, we used experimental evolution to determine how the selective pressure imposed by a protist predator impacts the virulence and fitness of a ubiquitous environmental opportunistic bacterial pathogen that has acquired multidrug resistance: Serratia marcescens. To this aim, we evolved S. marcescens in the presence or absence of generalist protist predator, Tetrahymena thermophila. After 60 d of evolution, we evaluated genotypic and phenotypic changes by comparing evolved S. marcescens with the ancestral strain. Whole-genome shotgun sequencing of the entire evolved populations and individual isolates revealed numerous cases of parallel evolution, many more than statistically expected by chance, in genes associated with virulence. Our phenotypic assays suggested that evolution in the presence of a predator maintained virulence, whereas evolution in the absence of a predator resulted in attenuated virulence. We also found a significant correlation between virulence, biofilm formation, growth, and grazing resistance. Overall, our results provide evidence that bacterial virulence and virulence-related traits are maintained by selective pressures imposed by protist predation.

From soil to sequence: filling the critical gap in genome-resolved metagenomics is essential to the future of soil microbial ecology

Soil microbiomes are heterogeneous, complex microbial communities. Metagenomic analysis is generating vast amounts of data, creating immense challenges in sequence assembly and analysis. Although advances in technology have resulted in the ability to easily collect large amounts of sequence data, soil samples containing thousands of unique taxa are often poorly characterized. These challenges reduce the usefulness of genome-resolved metagenomic (GRM) analysis seen in other fields of microbiology, such as the creation of high quality metagenomic assembled genomes and the adoption of genome scale modeling approaches. The absence of these resources restricts the scale of future research, limiting hypothesis generation and the predictive modeling of microbial communities. Creating publicly available databases of soil MAGs, similar to databases produced for other microbiomes, has the potential to transform scientific insights about soil microbiomes without requiring the computational resources and domain expertise for assembly and binning.

Vitamin B12 as a source of variability in isotope effects for chloroform biotransformation by Dehalobacter

Carbon and chlorine isotope effects for biotransformation of chloroform by different microbes show significant variability. Reductive dehalogenases (RDase) enzymes contain different cobamides, affecting substrate preferences, growth yields, and dechlorination rates and extent. We investigate the role of cobamide type on carbon and chlorine isotopic signals observed during reductive dechlorination of chloroform by the RDase CfrA. Microcosm experiments with two subcultures of a Dehalobacter-containing culture expressing CfrA-one with exogenous cobamide (Vitamin B12, B12+) and one without (to drive native cobamide production)-resulted in a markedly smaller carbon isotope enrichment factor (εC, bulk) for B12- (-22.1 ± 1.9‰) compared to B12+ (-26.8 ± 3.2‰). Both cultures exhibited significant chlorine isotope fractionation, and although a lower εCl, bulk was observed for B12- (-6.17 ± 0.72‰) compared to B12+ (-6.86 ± 0.77‰) cultures, these values are not statistically different. Importantly, dual-isotope plots produced identical slopes of ΛCl/C (ΛCl/C, B12+ = 3.41 ± 0.15, ΛCl/C, B12- = 3.39 ± 0.15), suggesting the same reaction mechanism is involved in both experiments, independent of the lower cobamide bases. A nonisotopically fractionating masking effect may explain the smaller fractionations observed for the B12- containing culture.

Aspergillus brasiliensis E_15.1: A Novel Thermophilic Endophyte from a Volcanic Crater Unveiled through Comprehensive Genome-Wide, Phenotypic Analysis, and Plant Growth-Promoting Trails

Thermophilic fungi have been seldom studied despite the fact that they can contribute to understanding ecological mechanisms of adaptation in diverse environments and have attractive toolboxes with a wide range of biotechnological applications. This work describes for the first time an endophytic and thermophilic strain of Aspergillus brasiliensis that was isolated in the crater of the active volcano "El Chichonal" in Mexico. This strain was capable of surviving in soil with a temperature of 60 °C and a pH of neutral acidity, which preluded a high thermostability and a potential in industrial application. The complete genome of A. brasiliensis E_15.1 was sequenced and assembled in 37 Mb of genomic DNA. We performed a comprehensive phylogenomic analysis for the precise taxonomic identification of this species as a novel strain of Aspergillus brasiliensis. Likewise, the predicted coding sequences were classified according to various functions including Carbohydrate-Active Enzymes (CAZymes), biosynthetic gene clusters of secondary metabolites (BGCs), and metabolic pathways associated with plant growth promotion. A. brasiliensis E_15.1 was found to degrade chitin, chitooligosaccharides, xylan, and cellulose. The genes to biosynthesize clavaric acid (a triterpene with antitumor activity) were found, thus probably having antitumor activity. In addition to the genomic analysis, a set of enzymatic assays confirmed the thermostability of extracellular xylanases and cellulases of A. brasiliensis E_15.1. The enzymatic repertoire of A. brasiliensis E_15.1 suggests that A. brasiliensis E_15.1 has a high potential for industrial application due to its thermostability and can promote plant growth at high temperatures. Finally, this strain constitutes an interesting source of terpenoids with pharmacological activity.

Metagenomic sequences from anaerobic chloroform and dichloromethane degrading microbial communities

Here, we present metagenomes from two cultures derived from an anaerobic microbial consortium used for bioremediation. One culture dechlorinates chloroform to dichloromethane, which is further mineralized to CO2. A second subculture was amended with only dichloromethane. We sought draft genomes of key microorganisms to identify metabolic potential in these consortia.