metaFlye: scalable long-read metagenome assembly using repeat graphs

Ecophysiological and genomic approaches to cyanobacterial hardening for restoration

Cyanobacteria inhabit extreme environments, including drylands, providing multiple benefits to the ecosystem. Soil degradation in warm drylands is increasing due to land use intensification. Restoration methods adapted to the extreme stress in drylands are being developed, such as cyanobacteria inoculation to recover biocrusts. For this type of restoration method to be a success, it is crucial to optimize the survival of inoculated cyanobacteria in the field. One strategy is to harden them to be acclimated to stressful conditions after laboratory culturing. Here, we analyzed the genome and ecophysiological response to osmotic desiccation and UVR stresses of an Antarctic cyanobacterium, Stenomitos frigidus ULC029, which is closely related to other cyanobacteria from warm and cold dryland soils. Chlorophyll a concentrations showed that preculturing ULC029 under moderate osmotic stress improved its survival during an assay of desiccation plus rehydration under UVR. Additionally, its sequential exposure to these stress factors increased the production of exopolysaccharides, carotenoids, and scytonemin. Desiccation, but not osmotic stress, increased the concentrations of the osmoprotectants trehalose and sucrose. However, osmotic stress might induce the production of other osmoprotectants, for which the complete pathways were observed in the ULC029 genome. In total, 140 genes known to be involved in stress resistance were annotated. Here, we confirm that the sequential application of moderate osmotic stress and dehydration could improve cyanobacterial hardening for soil restoration by inducing several resistance mechanisms. We provide a high-quality genome of ULC029 and a description of the main resistance mechanisms (i.e., production of exopolysaccharides, osmoprotectants, chlorophyll, and carotenoids; DNA repair; and oxidative stress protection).

Comprehensive assessment of 11 de novo HiFi assemblers on complex eukaryotic genomes and metagenomes

Pacific Biosciences (PacBio) HiFi sequencing technology generates long reads (>10 kbp) with very high accuracy (<0.01% sequencing error). Although several de novo assembly tools are available for HiFi reads, there are no comprehensive studies on the evaluation of these assemblers. We evaluated the performance of 11 de novo HiFi assemblers on (1) real data for three eukaryotic genomes; (2) 34 synthetic data sets with different ploidy, sequencing coverage levels, heterozygosity rates, and sequencing error rates; (3) one real metagenomic data set; and (4) five synthetic metagenomic data sets with different composition abundance and heterozygosity rates. The 11 assemblers were evaluated using quality assessment tool (QUAST) and benchmarking universal single-copy ortholog (BUSCO). We also used several additional criteria, namely, completion rate, single-copy completion rate, duplicated completion rate, average proportion of largest category, average distance difference, quality value, run-time, and memory utilization. Results show that hifiasm and hifiasm-meta should be the first choice for assembling eukaryotic genomes and metagenomes with HiFi data. We performed a comprehensive benchmarking study of commonly used assemblers on complex eukaryotic genomes and metagenomes. Our study will help the research community to choose the most appropriate assembler for their data and identify possible improvements in assembly algorithms.

Nanopore adaptive sampling effectively enriches bacterial plasmids

Bacterial plasmids play a major role in the spread of antibiotic resistance genes. However, their characterization via DNA sequencing suffers from the low abundance of plasmid DNA in those samples. Although sample preparation methods can enrich the proportion of plasmid DNA before sequencing, these methods are expensive and laborious, and they might introduce a bias by enriching only for specific plasmid DNA sequences. Nanopore adaptive sampling could overcome these issues by rejecting uninteresting DNA molecules during the sequencing process. In this study, we assess the application of adaptive sampling for the enrichment of low-abundant plasmids in known bacterial isolates using two different adaptive sampling tools. We show that a significant enrichment can be achieved even on expired flow cells. By applying adaptive sampling, we also improve the quality of de novo plasmid assemblies and reduce the sequencing time. However, our experiments also highlight issues with adaptive sampling if target and non-target sequences span similar regions.

IMPORTANCE

Antimicrobial resistance causes millions of deaths every year. Mobile genetic elements like bacterial plasmids are key drivers for the dissemination of antimicrobial resistance genes. This makes the characterization of plasmids via DNA sequencing an important tool for clinical microbiologists. Since plasmids are often underrepresented in bacterial samples, plasmid sequencing can be challenging and laborious. To accelerate the sequencing process, we evaluate nanopore adaptive sampling as an in silico method for the enrichment of low-abundant plasmids. Our results show the potential of this cost-efficient method for future plasmid research but also indicate issues that arise from using reference sequences.

Refinement of the "Candidatus Accumulibacter" genus based on metagenomic analysis of biological nutrient removal (BNR) pilot-scale plants operated with reduced aeration

Members of the "Candidatus Accumulibacter" genus are widely studied as key polyphosphate-accumulating organisms (PAOs) in biological nutrient removal (BNR) facilities performing enhanced biological phosphorus removal (EBPR). This diverse lineage includes 18 "Ca. Accumulibacter" species, which have been proposed based on the phylogenetic divergence of the polyphosphate kinase 1 (ppk1) gene and genome-scale comparisons of metagenome-assembled genomes (MAGs). Phylogenetic classification based on the 16S rRNA genetic marker has been difficult to attain because most "Ca. Accumulibacter" MAGs are incomplete and often do not include the rRNA operon. Here, we investigate the "Ca. Accumulibacter" diversity in pilot-scale treatment trains performing BNR under low dissolved oxygen (DO) conditions using genome-resolved metagenomics. Using long-read sequencing, we recovered medium- and high-quality MAGs for 5 of the 18 "Ca. Accumulibacter" species, all with rRNA operons assembled, which allowed a reassessment of the 16S rRNA-based phylogeny of this genus and an analysis of phylogeny based on the 23S rRNA gene. In addition, we recovered a cluster of MAGs that based on 16S rRNA, 23S rRNA, ppk1, and genome-scale phylogenetic analyses do not belong to any of the currently recognized "Ca. Accumulibacter" species for which we propose the new species designation "Ca. Accumulibacter jenkinsii" sp. nov. Relative abundance evaluations of the genus across all pilot plant operations revealed that regardless of the operational mode, "Ca. A. necessarius" and "Ca. A. propinquus" accounted for more than 40% of the "Ca. Accumulibacter" community, whereas the newly proposed "Ca. A. jenkinsii" accounted for about 5% of the "Ca. Accumulibacter" community.IMPORTANCEOne of the main drivers of energy use and operational costs in activated sludge processes is the amount of oxygen provided to enable biological phosphorus and nitrogen removal. Wastewater treatment facilities are increasingly considering reduced aeration to decrease energy consumption, and whereas successful BNR has been demonstrated in systems with minimal aeration, an adequate understanding of the microbial communities that facilitate nutrient removal under these conditions is still lacking. In this study, we used genome-resolved metagenomics to evaluate the diversity of the "Candidatus Accumulibacter" genus in pilot-scale plants operating with minimal aeration. We identified the "Ca. Accumulibacter" species enriched under these conditions, including one novel species for which we propose "Ca. Accumulibacter jenkinsii" sp. nov. as its designation. Furthermore, the MAGs obtained for five additional "Ca. Accumulibacter" species further refine the phylogeny of the "Ca. Accumulibacter" genus and provide new insight into its diversity within unconventional biological nutrient removal systems.

Short macrocyclic peptides in sponge genomes

Sponges (Porifera) contain many peptide-specialized metabolites with potent biological activities and significant roles in shaping marine ecology. It is well established that symbiotic bacteria produce bioactive "sponge" peptides, both on the ribosome (RiPPs) and nonribosomally. Here, we demonstrate that sponges themselves also produce many bioactive macrocyclic peptides, such as phakellistatins and related proline-rich macrocyclic peptides (PRMPs). Using the Stylissa carteri sponge transcriptome, methods were developed to find sequences encoding 46 distinct RiPP-type core peptides, of which ten encoded previously identified PRMP sequences. With this basis set, the genome and transcriptome of the sponge Axinella corrugata was interrogated to find 35 PRMP precursor peptides encoding 31 unique core peptide sequences. At least 11 of these produced cyclic peptides that were present in the sponge and could be characterized by mass spectrometry, including stylissamides A-D and seven previously undescribed compounds. Precursor peptides were encoded in the A. corrugata genome, confirming their animal origin. The peptides contained signal peptide sequences and highly repetitive recognition sequence-core peptide elements with up to 25 PRMP copies in a single precursor. In comparison to sponges without PRMPs, PRMP sponges are incredibly enriched in potentially secreted polypeptides, with >23,000 individual signal peptide encoding genes found in a single transcriptome. The similarities between PRMP biosynthetic genes and neuropeptides in terms of their biosynthetic logic suggest a fundamental biology linked to circular peptides, possibly indicating a widespread and underappreciated diversity of signaling peptide post-translational modifications across the animal kingdom.

Unveiling the Microbial Realm with VEBA 2.0: A modular bioinformatics suite for end-to-end genome-resolved prokaryotic, (micro)eukaryotic, and viral multi-omics from either short- or long-read sequencing

The microbiome is a complex community of microorganisms, encompassing prokaryotic (bacterial and archaeal), eukaryotic, and viral entities. This microbial ensemble plays a pivotal role in influencing the health and productivity of diverse ecosystems while shaping the web of life. However, many software suites developed to study microbiomes analyze only the prokaryotic community and provide limited to no support for viruses and microeukaryotes. Previously, we introduced the Viral Eukaryotic Bacterial Archaeal (VEBA) open-source software suite to address this critical gap in microbiome research by extending genome-resolved analysis beyond prokaryotes to encompass the understudied realms of eukaryotes and viruses. Here we present VEBA 2.0 with key updates including a comprehensive clustered microeukaryotic protein database, rapid genome/protein-level clustering, bioprospecting, non-coding/organelle gene modeling, genome-resolved taxonomic/pathway profiling, long-read support, and containerization. We demonstrate VEBA's versatile application through the analysis of diverse case studies including marine water, Siberian permafrost, and white-tailed deer lung tissues with the latter showcasing how to identify integrated viruses. VEBA represents a crucial advancement in microbiome research, offering a powerful and accessible platform that bridges the gap between genomics and biotechnological solutions.

MCSS: microbial community simulator based on structure

De novo assembly plays a pivotal role in metagenomic analysis, and the incorporation of third-generation sequencing technology can significantly improve the integrity and accuracy of assembly results. Recently, with advancements in sequencing technology (Hi-Fi, ultra-long), several long-read-based bioinformatic tools have been developed. However, the validation of the performance and reliability of these tools is a crucial concern. To address this gap, we present MCSS (microbial community simulator based on structure), which has the capability to generate simulated microbial community and sequencing datasets based on the structure attributes of real microbiome communities. The evaluation results indicate that it can generate simulated communities that exhibit both diversity and similarity to actual community structures. Additionally, MCSS generates synthetic PacBio Hi-Fi and Oxford Nanopore Technologies (ONT) long reads for the species within the simulated community. This innovative tool provides a valuable resource for benchmarking and refining metagenomic analysis methods. Code available at: https://github.com/panlab-bio/mcss.

Oceanimonas pelagia sp. nov., a novel biosurfactant-producing and plastic-degrading potential bacterium isolated from marine coastal sediment

A marine bacterial strain, named NTOU-MSR1T, was isolated from marine sediment of northern coast of Taiwan. This bacterium was Gram-stain-negative, aerobic, and motile, with a single flagellum. Its rod-shaped cells measured approximately 0.5-0.6 µm in width and 1.8-2.0 μm in length. NTOU-MSR1T grew at temperatures ranging from 10 to 45 °C, optimally at 30 °C. The pH range for growth was 7.0-10.0, with optimal growth at pH 7.0-8.0. It tolerated NaCl concentrations up to 12%. The cell membrane predominantly contained fatty acids such C16:1ω7c, C18:1ω7c, and C16:0. The overall genome relatedness indices indicated that strain NTOU-MSR1T had an average nucleotide identity (ANI) of 87.88% and a digital DNA-DNA hybridization (dDDH) value of 35.40% compared to its closest related species, O. marisflavi 102-Na3T. These values fell below the 95% and 70% threshold for species delineation, respectively. These findings suggested that the strain NTOU-MSR1T was a new member of the Oceanimonas genus. Its genomic DNA had a G + C content of 61.0 mol%. Genomic analysis revealed genes associated with the catechol branch of β- ketoadipate pathway for degrading polycyclic aromatic hydrocarbons, resistance to heavy metal, biosynthesis of polyhydroxybutyrate and the production of glycoside hydrolases (GH19, GH23, and GH103) for chitin and glycan digestion. Additionally, NTOU-MSR1T was capable of synthesizing biosurfactants and potentially degrading plastic. The proposed name for this new species is Oceanimonas pelagia, with the type strain designated as NTOU-MSR1T (= BCRC 81403T = JCM 36023T).

High-Quality Genome Assemblies of 4 Members of the Podospora anserina Species Complex

The filamentous fungus Podospora anserina is a model organism used extensively in the study of molecular biology, senescence, prion biology, meiotic drive, mating-type chromosome evolution, and plant biomass degradation. It has recently been established that P. anserina is a member of a complex of 7 closely related species. In addition to P. anserina, high-quality genomic resources are available for 2 of these taxa. Here, we provide chromosome-level annotated assemblies of the 4 remaining species of the complex, as well as a comprehensive data set of annotated assemblies from a total of 28 Podospora genomes. We find that all 7 species have genomes of around 35 Mb arranged in 7 chromosomes that are mostly collinear and less than 2% divergent from each other at genic regions. We further attempt to resolve their phylogenetic relationships, finding significant levels of phylogenetic conflict as expected from a rapid and recent diversification.

Macrolide resistance in Mannheimia haemolytica isolates associated with bovine respiratory disease from the German national resistance monitoring program GERM-Vet 2009 to 2020

Data collected from the German national resistance monitoring program GERM-Vet showed slowly increasing prevalence of macrolide resistance among bovine respiratory disease (BRD)-associated Pasteurellacae from cattle over the last decade. The focus of this study was to analyze the genetic basis of antimicrobial resistance (AMR) and the prevalence of multidrug-resistance (MDR)-mediating integrative and conjugative elements (ICEs) in 13 German BRD-associated Mannheimia haemolytica isolates collected between 2009 and 2020 via whole-genome sequencing. Antimicrobial susceptibility testing (AST) was performed via broth microdilution according to the recommendations of the Clinical and Laboratory Standards Institute for the macrolides erythromycin, tilmicosin, tulathromycin, gamithromycin, tildipirosin, and tylosin as well as 25 other antimicrobial agents. All isolates either had elevated MICs or were resistant to at least one of the macrolides tested. Analysis of whole-genome sequences obtained by hybrid assembly of Illumina MiSeq and Oxford Nanopore MinION reads revealed the presence of seven novel Tn7406-like ICEs, designated Tn7694, and Tn7724- Tn7729. These ICEs harbored the antimicrobial resistance genes erm(T), mef (C), mph(G), floR, catA3, aad(3")(9), aph(3')-Ia, aac(3)-IIa, strA, strB, tet(Y), and sul2 in different combinations. In addition, mutational changes conferring resistance to macrolides, nalidixic acid or streptomycin, respectively, were detected among the M. haemolytica isolates. In addition, four isolates carried a 4,613-bp plasmid with the β-lactamase gene blaROB - 1. The detection of the macrolide resistance genes erm(T), mef (C), and mph(G) together with other resistance genes on MDR-mediating ICEs in bovine M. haemolytica may explain the occurrence of therapeutic failure when treating BRD with regularly used antimicrobial agents, such as phenicols, penicillins, tetracyclines, or macrolides. Finally, pathogen identification and subsequent AST is essential to ensure the efficacy of the antimicrobial agents applied to control BRD in cattle.

Xanthomonas rydalmerensis sp. nov., a non-pathogenic member of Group 1 Xanthomonas

Five bacterial isolates were isolated from Fragaria × ananassa in 1976 in Rydalmere, Australia, during routine biosecurity surveillance. Initially, the results of biochemical characterisation indicated that these isolates represented members of the genus Xanthomonas. To determine their species, further analysis was conducted using both phenotypic and genotypic approaches. Phenotypic analysis involved using MALDI-TOF MS and BIOLOG GEN III microplates, which confirmed that the isolates represented members of the genus Xanthomonas but did not allow them to be classified with respect to species. Genome relatedness indices and the results of extensive phylogenetic analysis confirmed that the isolates were members of the genus Xanthomonas and represented a novel species. On the basis the minimal presence of virulence-associated factors typically found in genomes of members of the genus Xanthomonas, we suggest that these isolates are non-pathogenic. This conclusion was supported by the results of a pathogenicity assay. On the basis of these findings, we propose the name Xanthomonas rydalmerensis, with DAR 34855T = ICMP 24941 as the type strain.

Novel functional insights into the microbiome inhabiting marine plastic debris: critical considerations to counteract the challenges of thin biofilms using multi-omics and comparative metaproteomics

BACKGROUND

Microbial functioning on marine plastic surfaces has been poorly documented, especially within cold climates where temperature likely impacts microbial activity and the presence of hydrocarbonoclastic microorganisms. To date, only two studies have used metaproteomics to unravel microbial genotype-phenotype linkages in the marine 'plastisphere', and these have revealed the dominance of photosynthetic microorganisms within warm climates. Advancing the functional representation of the marine plastisphere is vital for the development of specific databases cataloging the functional diversity of the associated microorganisms and their peptide and protein sequences, to fuel biotechnological discoveries. Here, we provide a comprehensive assessment for plastisphere metaproteomics, using multi-omics and data mining on thin plastic biofilms to provide unique insights into plastisphere metabolism. Our robust experimental design assessed DNA/protein co-extraction and cell lysis strategies, proteomics workflows, and diverse protein search databases, to resolve the active plastisphere taxa and their expressed functions from an understudied cold environment.

RESULTS

For the first time, we demonstrate the predominance and activity of hydrocarbonoclastic genera (Psychrobacter, Flavobacterium, Pseudomonas) within a primarily heterotrophic plastisphere. Correspondingly, oxidative phosphorylation, the citrate cycle, and carbohydrate metabolism were the dominant pathways expressed. Quorum sensing and toxin-associated proteins of Streptomyces were indicative of inter-community interactions. Stress response proteins expressed by Psychrobacter, Planococcus, and Pseudoalteromonas and proteins mediating xenobiotics degradation in Psychrobacter and Pseudoalteromonas suggested phenotypic adaptations to the toxic chemical microenvironment of the plastisphere. Interestingly, a targeted search strategy identified plastic biodegradation enzymes, including polyamidase, hydrolase, and depolymerase, expressed by rare taxa. The expression of virulence factors and mechanisms of antimicrobial resistance suggested pathogenic genera were active, despite representing a minor component of the plastisphere community.

CONCLUSION

Our study addresses a critical gap in understanding the functioning of the marine plastisphere, contributing new insights into the function and ecology of an emerging and important microbial niche. Our comprehensive multi-omics and comparative metaproteomics experimental design enhances biological interpretations to provide new perspectives on microorganisms of potential biotechnological significance beyond biodegradation and to improve the assessment of the risks associated with microorganisms colonizing marine plastic pollution. Video Abstract.

Genomic analysis of diverse environmental Acinetobacter isolates identifies plasmids, antibiotic resistance genes, and capsular polysaccharides shared with clinical strains

Acinetobacter baumannii, an important pathogen known for its widespread antibiotic resistance, has been the focus of extensive research within its genus, primarily involving clinical isolates. Consequently, data on environmental A. baumannii and other Acinetobacter species remain limited. Here, we utilized Illumina and Nanopore sequencing to analyze the genomes of 10 Acinetobacter isolates representing 6 different species sourced from aquatic environments in South Australia. All 10 isolates were phylogenetically distinct compared to clinical and other non-clinical Acinetobacter strains, often tens of thousands of single-nucleotide polymorphisms from their nearest neighbors. Despite the genetic divergence, we identified pdif modules (sections of mobilized DNA) carrying clinically important antimicrobial resistance genes in species other than A. baumannii, including carbapenemase oxa58, tetracycline resistance gene tet(39), and macrolide resistance genes msr(E)-mph(E). These pdif modules were located on plasmids with high sequence identity to those circulating in globally distributed A. baumannii ST1 and ST2 clones. The environmental A. baumannii isolate characterized here (SAAb472; ST350) did not possess any native plasmids; however, it could capture two clinically important plasmids (pRAY and pACICU2) with high transfer frequencies. Furthermore, A. baumannii SAAb472 possessed virulence genes and a capsular polysaccharide type analogous to clinical strains. Our findings highlight the potential for environmental Acinetobacter species to acquire and disseminate clinically important antimicrobial resistance genes, underscoring the need for further research into the ecology and evolution of this important genus.IMPORTANCEAntimicrobial resistance (AMR) is a global threat to human, animal, and environmental health. Studying AMR in environmental bacteria is crucial to understand the emergence and dissemination of resistance genes and pathogens, and to identify potential reservoirs and transmission routes. This study provides novel insights into the genomic diversity and AMR potential of environmental Acinetobacter species. By comparing the genomes of aquatic Acinetobacter isolates with clinical and non-clinical strains, we revealed that they are highly divergent yet carry pdif modules that encode resistance to antibiotics commonly used in clinical settings. We also demonstrated that an environmental A. baumannii isolate can acquire clinically relevant plasmids and carries virulence factors similar to those of hospital-associated strains. These findings suggest that environmental Acinetobacter species may serve as reservoirs and vectors of clinically important genes. Consequently, further research is warranted to comprehensively understand the ecology and evolution of this genus.

Particle-attached bacteria act as gatekeepers in the decomposition of complex phytoplankton polysaccharides

BACKGROUND

Marine microalgae (phytoplankton) mediate almost half of the worldwide photosynthetic carbon dioxide fixation and therefore play a pivotal role in global carbon cycling, most prominently during massive phytoplankton blooms. Phytoplankton biomass consists of considerable proportions of polysaccharides, substantial parts of which are rapidly remineralized by heterotrophic bacteria. We analyzed the diversity, activity, and functional potential of such polysaccharide-degrading bacteria in different size fractions during a diverse spring phytoplankton bloom at Helgoland Roads (southern North Sea) at high temporal resolution using microscopic, physicochemical, biodiversity, metagenome, and metaproteome analyses.

RESULTS

Prominent active 0.2-3 µm free-living clades comprised Aurantivirga, "Formosa", Cd. Prosiliicoccus, NS4, NS5, Amylibacter, Planktomarina, SAR11 Ia, SAR92, and SAR86, whereas BD1-7, Stappiaceae, Nitrincolaceae, Methylophagaceae, Sulfitobacter, NS9, Polaribacter, Lentimonas, CL500-3, Algibacter, and Glaciecola dominated 3-10 µm and > 10 µm particles. Particle-attached bacteria were more diverse and exhibited more dynamic adaptive shifts over time in terms of taxonomic composition and repertoires of encoded polysaccharide-targeting enzymes. In total, 305 species-level metagenome-assembled genomes were obtained, including 152 particle-attached bacteria, 100 of which were novel for the sampling site with 76 representing new species. Compared to free-living bacteria, they featured on average larger metagenome-assembled genomes with higher proportions of polysaccharide utilization loci. The latter were predicted to target a broader spectrum of polysaccharide substrates, ranging from readily soluble, simple structured storage polysaccharides (e.g., laminarin, α-glucans) to less soluble, complex structural, or secreted polysaccharides (e.g., xylans, cellulose, pectins). In particular, the potential to target poorly soluble or complex polysaccharides was more widespread among abundant and active particle-attached bacteria.

CONCLUSIONS

Particle-attached bacteria represented only 1% of all bloom-associated bacteria, yet our data suggest that many abundant active clades played a pivotal gatekeeping role in the solubilization and subsequent degradation of numerous important classes of algal glycans. The high diversity of polysaccharide niches among the most active particle-attached clades therefore is a determining factor for the proportion of algal polysaccharides that can be rapidly remineralized during generally short-lived phytoplankton bloom events. Video Abstract.

Complete genome sequences of 12 bacterial strains from the honey bee gut, resolved with long-read nanopore sequencing

We report the de novo sequencing of six bacterial strains isolated from the Western honey bee, as well as the resequencing of six strains that have existing draft genomes, to obtain complete, chromosomal-level assemblies. These strains include the bee gut symbiont genera Bartonella, Bifidobacterium, Snodgrassella, Gilliamella, Lactobacillus, and the opportunistic pathogen Serratia marcescens KZ11.

MarkerScan: Separation and assembly of cobionts sequenced alongside target species in biodiversity genomics projects

Contamination of public databases by mislabelled sequences has been highlighted for many years and the avalanche of novel sequencing data now being deposited has the potential to make databases difficult to use effectively. It is therefore crucial that sequencing projects and database curators perform pre-submission checks to remove obvious contamination and avoid propagating erroneous taxonomic relationships. However, it is important also to recognise that biological contamination of a target sample with unexpected species' DNA can also lead to the discovery of fascinating biological phenomena through the identification of environmental organisms or endosymbionts. Here, we present a novel, integrated method for detection and generation of high-quality genomes of all non-target genomes co-sequenced in eukaryotic genome sequencing projects. After performing taxonomic profiling of an assembly from the raw data, and leveraging the identity of small rRNA sequences discovered therein as markers, a targeted classification approach retrieves and assembles high-quality genomes. The genomes of these cobionts are then not only removed from the target species' genome but also available for further interrogation. Source code is available from https://github.com/CobiontID/MarkerScan. MarkerScan is written in Python and is deployed as a Docker container.

Recovering high-quality bacterial genomes from cross-contaminated cultures: a case study of marine Vibrio campbellii

BACKGROUND

Environmental monitoring of bacterial pathogens is critical for disease control in coastal marine ecosystems to maintain animal welfare and ecosystem function and to prevent significant economic losses. This requires accurate taxonomic identification of environmental bacterial pathogens, which often cannot be achieved by commonly used genetic markers (e.g., 16S rRNA gene), and an understanding of their pathogenic potential based on the information encoded in their genomes. The decreasing costs of whole genome sequencing (WGS), combined with newly developed bioinformatics tools, now make it possible to unravel the full potential of environmental pathogens, beyond traditional microbiological approaches. However, obtaining a high-quality bacterial genome, requires initial cultivation in an axenic culture, which is a bottleneck in environmental microbiology due to cross-contamination in the laboratory or isolation of non-axenic strains.

RESULTS

We applied WGS to determine the pathogenic potential of two Vibrio isolates from coastal seawater. During the analysis, we identified cross-contamination of one of the isolates and decided to use this dataset to evaluate the possibility of bioinformatic contaminant removal and recovery of bacterial genomes from a contaminated culture. Despite the contamination, using an appropriate bioinformatics workflow, we were able to obtain high quality and highly identical genomes (Average Nucleotide Identity value 99.98%) of one of the Vibrio isolates from both the axenic and the contaminated culture. Using the assembled genome, we were able to determine that this isolate belongs to a sub-lineage of Vibrio campbellii associated with several diseases in marine organisms. We also found that the genome of the isolate contains a novel Vibrio plasmid associated with bacterial defense mechanisms and horizontal gene transfer, which may offer a competitive advantage to this putative pathogen.

CONCLUSIONS

Our study shows that, using state-of-the-art bioinformatics tools and a sufficient sequencing effort, it is possible to obtain high quality genomes of the bacteria of interest and perform in-depth genomic analyses even in the case of a contaminated culture. With the new isolate and its complete genome, we are providing new insights into the genomic characteristics and functional potential of this sub-lineage of V. campbellii. The approach described here also highlights the possibility of recovering complete bacterial genomes in the case of non-axenic cultures or obligatory co-cultures.

In Silico Analysis of Shiga Toxin-Producing Escherichia coli O157:H7 Strains from Presumptive Super- and Low-Shedder Cattle

Cattle are the primary reservoir for STEC O157, with some shedding >104 CFU/g in feces, a phenomenon known as super-shedding (SS). The mechanism(s) responsible for SS are not understood but have been attributed to the environment, host, and pathogen. This study aimed to compare genetic characteristics of STEC O157 strains from cattle in the same commercial feedlot pens with SS or low-shedding (LS) status. Strains from SS (n = 35) and LS (n = 28) collected from 11 pens in three feedlots were analyzed for virulence genes, Shiga toxin-carrying bacteriophage insertion sites, and phylogenetic relationships. In silico analysis showed limited variation regarding virulence gene profiles. Stx-encoding prophage insertion sites mrlA and wrbA for stx1a and stx2a, respectively, were all occupied, but two isolates had fragments of the stx-carrying phage in mrlA and wrbA loci without stx1a and stx2a. All strains screened for lineage-specific polymorphism assay (LSPA-6) were 111111, lineage I. Of the isolates, 61 and 2 were clades 1 and 8, respectively. Phylogenetic analysis revealed that pens with more than one SS had multiple distantly related clusters of SS and LS isolates. Although virulence genes and lineage were largely similar within and across feedlots, multiple genetic origins of strains within a single feedlot pen illustrate challenges for on-farm control of STEC.

Closing the genome of unculturable cable bacteria using a combined metagenomic assembly of long and short sequencing reads

Many environmentally relevant micro-organisms cannot be cultured, and even with the latest metagenomic approaches, achieving complete genomes for specific target organisms of interest remains a challenge. Cable bacteria provide a prominent example of a microbial ecosystem engineer that is currently unculturable. They occur in low abundance in natural sediments, but due to their capability for long-distance electron transport, they exert a disproportionately large impact on the biogeochemistry of their environment. Current available genomes of marine cable bacteria are highly fragmented and incomplete, hampering the elucidation of their unique electrogenic physiology. Here, we present a metagenomic pipeline that combines Nanopore long-read and Illumina short-read shotgun sequencing. Starting from a clonal enrichment of a cable bacterium, we recovered a circular metagenome-assembled genome (5.09 Mbp in size), which represents a novel cable bacterium species with the proposed name Candidatus Electrothrix scaldis. The closed genome contains 1109 novel identified genes, including key metabolic enzymes not previously described in incomplete genomes of cable bacteria. We examined in detail the factors leading to genome closure. Foremost, native, non-amplified long reads are crucial to resolve the many repetitive regions within the genome of cable bacteria, and by analysing the whole metagenomic assembly, we found that low strain diversity is key for achieving genome closure. The insights and approaches presented here could help achieve genome closure for other keystone micro-organisms present in complex environmental samples at low abundance.

mEnrich-seq: methylation-guided enrichment sequencing of bacterial taxa of interest from microbiome

Metagenomics has enabled the comprehensive study of microbiomes. However, many applications would benefit from a method that sequences specific bacterial taxa of interest, but not most background taxa. We developed mEnrich-seq (in which 'm' stands for methylation and seq for sequencing) for enriching taxa of interest from metagenomic DNA before sequencing. The core idea is to exploit the self versus nonself differentiation by natural bacterial DNA methylation and rationally choose methylation-sensitive restriction enzymes, individually or in combination, to deplete host and background taxa while enriching targeted taxa. This idea is integrated with library preparation procedures and applied in several applications to enrich (up to 117-fold) pathogenic or beneficial bacteria from human urine and fecal samples, including species that are hard to culture or of low abundance. We assessed 4,601 bacterial strains with mapped methylomes so far and showed broad applicability of mEnrich-seq. mEnrich-seq provides microbiome researchers with a versatile and cost-effective approach for selective sequencing of diverse taxa of interest.

Genomic surveillance for antimicrobial resistance - a One Health perspective

Antimicrobial resistance (AMR) - the ability of microorganisms to adapt and survive under diverse chemical selection pressures - is influenced by complex interactions between humans, companion and food-producing animals, wildlife, insects and the environment. To understand and manage the threat posed to health (human, animal, plant and environmental) and security (food and water security and biosecurity), a multifaceted 'One Health' approach to AMR surveillance is required. Genomic technologies have enabled monitoring of the mobilization, persistence and abundance of AMR genes and mutations within and between microbial populations. Their adoption has also allowed source-tracing of AMR pathogens and modelling of AMR evolution and transmission. Here, we highlight recent advances in genomic AMR surveillance and the relative strengths of different technologies for AMR surveillance and research. We showcase recent insights derived from One Health genomic surveillance and consider the challenges to broader adoption both in developed and in lower- and middle-income countries.

Direct sequencing of insect symbionts via nanopore adaptive sampling

Insect symbionts can alter their host phenotype and their effects can range from beneficial to pathogenic. Moreover, many insects exhibit co-infections, making their study more challenging. Less than 1% of insect species have high-quality referenced genomes available and fewer still also have their symbionts sequenced. Two methods are commonly used to sequence symbionts: whole-genome sequencing to concomitantly capture the host and bacterial genomes, or isolation of the symbiont's genome before sequencing. These methods are limited when dealing with rare or poorly characterized symbionts. Long-read technology is an important tool to generate high-quality genomes as they can overcome high levels of heterozygosity, repeat content, and transposable elements that confound short-read methods. Oxford Nanopore (ONT) adaptive sampling allows a sequencing instrument to select or reject sequences in real time. We describe a method based on ONT adaptive sampling (subtractive) approach that readily permitted the sequencing of the complete genomes of mitochondria, Buchnera and its plasmids (pLeu, pTrp), and Wolbachia genomes in two aphid species, Aphis glycines and Pentalonia nigronervosa. Adaptive sampling is able to retrieve organelles such as mitochondria and symbionts that have high representation in their hosts such as Buchnera and Wolbachia, but is less successful at retrieving symbionts in low concentrations.

Host DNA depletion methods and genome-centric metagenomics of bovine hindmilk microbiome

Bovine mastitis is a multi-etiological and complex disease, resulting in serious economic consequences for dairy farmers and industry. In recent years, the microbiological evaluation of raw milk has been investigated in-depth using next-generation sequencing approaches such as metataxonomic analysis. Despite this, host DNA is a major concern in the shotgun metagenomic sequencing of microbial communities in milk samples, and it represents a big challenge. In this study, we aimed to evaluate different methods for host DNA depletion and/or microbial DNA enrichment and assess the use of PCR-based whole genome amplification in milk samples with high somatic cell count (SCC) by using short- and long-read sequencing technologies. Our results evidenced that DNA extraction performed differently in terms of host DNA removal, impacting metagenome composition and functional profiles.. Moreover, the ratio of SCC/bacteria ultimately impacts microbial DNA yield, and samples with low SCC (SCC below 100,000 cells/mL) are the most problematic. When milk samples with high SCC (SCC above 200,000 cells/mL) underwent multiple-displacement amplification (MDA), we successfully recovered high-quality metagenome-assembled genomes (MAGs), and long-read sequencing was feasible even for samples with low DNA concentration. By associating MDA and short-read sequencing, we recovered two times more MAGs than in untreated samples, and an ongoing co-infection not reported by traditional methods was detected for mastitis pathogen. Overall, this new approach will improve the detection of mastitis-associated microorganisms and make it possible to examine host-microbiome interactions in bovine mastitis.IMPORTANCENext-generation sequencing technologies have been widely used to gain new insights into the diversity of the microbial community of milk samples and dairy products for different purposes such as microbial safety, profiling of starter cultures, and host-microbiome interactions. Milk is a complex food matrix, and additionally, the presence of host nucleic acid sequences is considered a contaminant in untargeted high-throughput sequencing studies. Therefore, genomic-centric metagenomic studies of milk samples focusing on the health-disease status in dairy cattle are still scarce, which makes it difficult to evaluate the microbial ecophysiology of bovine hindmilk. This study provides an alternative method for genome-centric metagenome studies applied to hindmilk samples with high somatic cell content, which is indispensable to examining host-microbiome interactions in bovine mastitis.

Reference-free Structural Variant Detection in Microbiomes via Long-read Coassembly Graphs

Bacterial genome dynamics are vital for understanding the mechanisms underlying microbial adaptation, growth, and their broader impact on host phenotype. Structural variants (SVs), genomic alterations of 10 base pairs or more, play a pivotal role in driving evolutionary processes and maintaining genomic heterogeneity within bacterial populations. While SV detection in isolate genomes is relatively straightforward, metagenomes present broader challenges due to absence of clear reference genomes and presence of mixed strains. In response, our proposed method rhea, forgoes reference genomes and metagenome-assembled genomes (MAGs) by encompassing a single metagenome coassembly graph constructed from all samples in a series. The log fold change in graph coverage between subsequent samples is then calculated to call SVs that are thriving or declining throughout the series. We show rhea to outperform existing methods for SV and horizontal gene transfer (HGT) detection in two simulated mock metagenomes, which is particularly noticeable as the simulated reads diverge from reference genomes and an increase in strain diversity is incorporated. We additionally demonstrate use cases for rhea on series metagenomic data of environmental and fermented food microbiomes to detect specific sequence alterations between subsequent time and temperature samples, suggesting host advantage. Our innovative approach leverages raw read patterns rather than references or MAGs to include all sequencing reads in analysis, and thus provide versatility in studying SVs across diverse and poorly characterized microbial communities for more comprehensive insights into microbial genome dynamics.

Unraveling metagenomics through long-read sequencing: a comprehensive review

The study of microbial communities has undergone significant advancements, starting from the initial use of 16S rRNA sequencing to the adoption of shotgun metagenomics. However, a new era has emerged with the advent of long-read sequencing (LRS), which offers substantial improvements over its predecessor, short-read sequencing (SRS). LRS produces reads that are several kilobases long, enabling researchers to obtain more complete and contiguous genomic information, characterize structural variations, and study epigenetic modifications. The current leaders in LRS technologies are Pacific Biotechnologies (PacBio) and Oxford Nanopore Technologies (ONT), each offering a distinct set of advantages. This review covers the workflow of long-read metagenomics sequencing, including sample preparation (sample collection, sample extraction, and library preparation), sequencing, processing (quality control, assembly, and binning), and analysis (taxonomic annotation and functional annotation). Each section provides a concise outline of the key concept of the methodology, presenting the original concept as well as how it is challenged or modified in the context of LRS. Additionally, the section introduces a range of tools that are compatible with LRS and can be utilized to execute the LRS process. This review aims to present the workflow of metagenomics, highlight the transformative impact of LRS, and provide researchers with a selection of tools suitable for this task.

Meteora sporadica, a protist with incredible cell architecture, is related to Hemimastigophora

"Kingdom-level" branches are being added to the tree of eukaryotes at a rate approaching one per year, with no signs of slowing down.1,2,3,4 Some are completely new discoveries, whereas others are morphologically unusual protists that were previously described but lacked molecular data. For example, Hemimastigophora are predatory protists with two rows of flagella that were known since the 19th century but proved to represent a new deep-branching eukaryote lineage when phylogenomic analyses were conducted.2Meteora sporadica5 is a protist with a unique morphology; cells glide over substrates along a long axis of anterior and posterior projections while a pair of lateral "arms" swing back and forth, a motility system without any obvious parallels. Originally, Meteora was described by light microscopy only, from a short-term enrichment of deep-sea sediment. A small subunit ribosomal RNA (SSU rRNA) sequence was reported recently, but the phylogenetic placement of Meteora remained unresolved.6 Here, we investigated two cultivated Meteora sporadica isolates in detail. Transmission electron microscopy showed that both the anterior-posterior projections and the arms are supported by microtubules originating from a cluster of subnuclear microtubule organizing centers (MTOCs). Neither have a flagellar axoneme-like structure. Sequencing the mitochondrial genome showed this to be among the most gene-rich known, outside jakobids. Remarkably, phylogenomic analyses of 254 nuclear protein-coding genes robustly support a close relationship with Hemimastigophora. Our study suggests that Meteora and Hemimastigophora together represent a morphologically diverse "supergroup" and thus are important for resolving the tree of eukaryote life and early eukaryote evolution.

Structural and Functional Disparities within the Human Gut Virome in Terms of Genome Topology and Representative Genome Selection

Circularity confers protection to viral genomes where linearity falls short, thereby fulfilling the form follows function aphorism. However, a shift away from morphology-based classification toward the molecular and ecological classification of viruses is currently underway within the field of virology. Recent years have seen drastic changes in the International Committee on Taxonomy of Viruses' operational definitions of viruses, particularly for the tailed phages that inhabit the human gut. After the abolition of the order Caudovirales, these tailed phages are best defined as members of the class Caudoviricetes. To determine the epistemological value of genome topology in the context of the human gut virome, we designed a set of seven experiments to assay the impact of genome topology and representative viral selection on biological interpretation. Using Oxford Nanopore long reads for viral genome assembly coupled with Illumina short-read polishing, we showed that circular and linear virus genomes differ remarkably in terms of genome quality, GC skew, transfer RNA gene frequency, structural variant frequency, cross-reference functional annotation (COG, KEGG, Pfam, and TIGRfam), state-of-the-art marker-based classification, and phage-host interaction. Furthermore, the disparity profile changes during dereplication. In particular, our phage-host interaction results demonstrated that proportional abundances cannot be meaningfully compared without due regard for genome topology and dereplication threshold, which necessitates the need for standardized reporting. As a best practice guideline, we recommend that comparative studies of the human gut virome always report the ratio of circular to linear viral genomes along with the dereplication threshold so that structural and functional metrics can be placed into context when assessing biologically relevant metagenomic properties such as proportional abundance.

Complete genome of Flavobacterium pectinovorum str. ZE23VCel01 obtained through Nanopore Q20+ chemistry

This study reports the complete genome of Flavobacterium pectinovorum str. ZE23VCel01 isolated from a freshwater environment. By means of Nanopore Q20+ chemistry, the chromosome was assembled as a circular element with a length of 5.9 Mbp, a GC content of 33.58%, and a coverage of 122×.

Comprehensive Analysis of the Complete Mitochondrial Genome of Rehmannia chingii: An Autotrophic Species in the Orobanchaceae Family

Rehmannia chingii is an important medicinal plant with immense value in scientific research. However, its mitochondrial genome (mitogenome) has not yet been characterized. Herein, based on whole-genome Illumina short reads and PacBio HiFi reads, we obtained the complete mitogenome of R. chingii through a de novo assembly strategy. We carried out comparative genomic analyses and found that, in comparison with the plastid genome (plastome) showing a high degree of structural conservation, the R. chingii mitogenome structure is relatively complex, showing an intricate ring structure with 16 connections, owing to five repetitive sequences. The R. chingii mitogenome was 783,161 bp with a GC content of 44.8% and contained 77 genes, comprising 47 protein-coding genes (CDS), 27 tRNA genes, and 3 rRNA genes. We counted 579 RNA editing events in 47 CDS and 12,828 codons in all CDSs of the R. chingii mitogenome. Furthermore, 24 unique sequence transfer fragments were found between the mitogenome and plastome, comprising 8 mitogenome CDS genes and 16 plastome CDS genes, corresponding to 2.39% of the R. chingii mitogenome. Mitogenomes had shorter but more collinear regions, evidenced by a comparison of the organelles of non-parasitic R. chingii, hemiparasitic Pedicularis chinensis, and holoparasitic Aeginetia indica in the Orobanchaceae family. Moreover, from non-parasitic to holoparasitic species, the genome size in the mitogenomes of Orobanchaceae species did not decrease gradually. Instead, the smallest mitogenome was found in the hemiparasitic species P. chinensis, with a size of 225,612 bp. The findings fill the gap in the mitogenome research of the medicinal plant R. chingii, promote the progress of the organelle genome research of the Orobanchaceae family, and provide clues for molecular breeding.

The origins of the Guinness stout yeast

Beer is made via the fermentation of an aqueous extract predominantly composed of malted barley flavoured with hops. The transforming microorganism is typically a single strain of Saccharomyces cerevisiae, and for the majority of major beer brands the yeast strain is a unique component. The present yeast used to make Guinness stout brewed in Dublin, Ireland, can be traced back to 1903, but its origins are unknown. To that end, we used Illumina and Nanopore sequencing to generate whole-genome sequencing data for a total of 22 S. cerevisiae yeast strains: 16 from the Guinness collection and 6 other historical Irish brewing. The origins of the Guinness yeast were determined with a SNP-based analysis, demonstrating that the Guinness strains occupy a distinct group separate from other historical Irish brewing yeasts. Assessment of chromosome number, copy number variation and phenotypic evaluation of key brewing attributes established Guinness yeast-specific SNPs but no specific chromosomal amplifications. Our analysis also demonstrated the effects of yeast storage on phylogeny. Altogether, our results suggest that the Guinness yeast used today is related to the first deposited Guinness yeast; the 1903 Watling Laboratory Guinness yeast.

IS26 drives the dissemination of bla CTX-M genes in an Ecuadorian community

IMPORTANCE

The horizontal gene transfer events are the major contributors to the current spread of CTX-M-encoding genes, the most common extended-spectrum β-lactamase (ESBL), and many clinically crucial antimicrobial resistance (AMR) genes. This study presents evidence of the critical role of IS26 transposable element for the mobility of bla CTX-M gene among Escherichia coli isolates from children and domestic animals in the community. We suggest that the nucleotide sequences of IS26-bla CTX-M could be used to study bla CTX-M transmission between humans, domestic animals, and the environment, because understanding of the dissemination patterns of AMR genes is critical to implement effective measures to slow down the dissemination of these clinically important genes.

Aerobic denitrification as an N2O source from microbial communities

Nitrous oxide (N2O) is a potent greenhouse gas of primarily microbial origin. Oxic and anoxic emissions are commonly ascribed to autotrophic nitrification and heterotrophic denitrification, respectively. Beyond this established dichotomy, we quantitatively show that heterotrophic denitrification can significantly contribute to aerobic nitrogen turnover and N2O emissions in complex microbiomes exposed to frequent oxic/anoxic transitions. Two planktonic, nitrification-inhibited enrichment cultures were established under continuous organic carbon and nitrate feeding, and cyclic oxygen availability. Over a third of the influent organic substrate was respired with nitrate as electron acceptor at high oxygen concentrations (>6.5 mg/L). N2O accounted for up to one-quarter of the nitrate reduced under oxic conditions. The enriched microorganisms maintained a constitutive abundance of denitrifying enzymes due to the oxic/anoxic frequencies exceeding their protein turnover-a common scenario in natural and engineered ecosystems. The aerobic denitrification rates are ascribed primarily to the residual activity of anaerobically synthesised enzymes. From an ecological perspective, the selection of organisms capable of sustaining significant denitrifying activity during aeration shows their competitive advantage over other heterotrophs under varying oxygen availabilities. Ultimately, we propose that the contribution of heterotrophic denitrification to aerobic nitrogen turnover and N2O emissions is currently underestimated in dynamic environments.

Chlamydiae as symbionts of photosynthetic dinoflagellates

Chlamydiae are ubiquitous intracellular bacteria and infect a wide diversity of eukaryotes, including mammals. However, chlamydiae have never been reported to infect photosynthetic organisms. Here, we describe a novel chlamydial genus and species, Candidatus Algichlamydia australiensis, capable of infecting the photosynthetic dinoflagellate Cladocopium sp. (originally isolated from a scleractinian coral). Algichlamydia australiensis was confirmed to be intracellular by fluorescence in situ hybridization and confocal laser scanning microscopy and temporally stable at the population level by monitoring its relative abundance across four weeks of host growth. Using a combination of short- and long-read sequencing, we recovered a high-quality (completeness 91.73% and contamination 0.27%) metagenome-assembled genome of A. australiensis. Phylogenetic analyses show that this chlamydial taxon represents a new genus and species within the Simkaniaceae family. Algichlamydia australiensis possesses all the hallmark genes for chlamydiae-host interactions, including a complete type III secretion system. In addition, a type IV secretion system is encoded on a plasmid and has previously been observed for only three other chlamydial species. Twenty orthologous groups of genes are unique to A. australiensis, one of which is structurally similar to a protein known from Cyanobacteria and Archaeplastida involved in thylakoid biogenesis and maintenance, hinting at potential chlamydiae interactions with the chloroplasts of Cladocopium cells. Our study shows that chlamydiae infect dinoflagellate symbionts of cnidarians, the first photosynthetic organism reported to harbor chlamydiae, thereby expanding the breadth of chlamydial hosts and providing a new contribution to the discussion around the role of chlamydiae in the establishment of the primary plastid.

Evolutionary history of tyrosine-supplementing endosymbionts in pollen-feeding beetles

Many insects feeding on nutritionally challenging diets like plant sap, leaves, or wood engage in ancient associations with bacterial symbionts that supplement limiting nutrients or produce digestive or detoxifying enzymes. However, the distribution, function, and evolutionary dynamics of microbial symbionts in insects exploiting other plant tissues or relying on a predacious diet remain poorly understood. Here, we investigated the evolutionary history and function of the intracellular gamma-proteobacterial symbiont "Candidatus Dasytiphilus stammeri" in soft-winged flower beetles (Coleoptera, Melyridae, Dasytinae) that transition from saprophagy or carnivory to palynivory (pollen-feeding) between larval and adult stage. Reconstructing the distribution of the symbiont within the Dasytinae phylogeny unraveled not only a long-term coevolution, originating from a single acquisition event with subsequent host-symbiont codiversification, but also several independent symbiont losses. The analysis of 20 different symbiont genomes revealed that their genomes are severely eroded. However, the universally retained shikimate pathway indicates that the core metabolic contribution to their hosts is the provisioning of tyrosine for cuticle sclerotization and melanization. Despite the high degree of similarity in gene content and order across symbiont strains, the capacity to synthesize additional essential amino acids and vitamins and to recycle urea is retained in some but not all symbionts, suggesting ecological differences among host lineages. This report of tyrosine-provisioning symbionts in insects with saprophagous or carnivorous larvae and pollen-feeding adults expands our understanding of tyrosine supplementation as an important symbiont-provided benefit across a broad range of insects with diverse feeding ecologies.

WGS of a cluster of MDR Shigella sonnei utilizing Oxford Nanopore R10.4.1 long-read sequencing

OBJECTIVES

To utilize long-read nanopore sequencing (R10.4.1 flowcells) for WGS of a cluster of MDR Shigella sonnei, specifically characterizing genetic predictors of antimicrobial resistance (AMR).

METHODS

WGS was performed on S. sonnei isolates identified from stool and blood between September 2021 and October 2022. Bacterial DNA from clinical isolates was extracted on the MagNA Pure 24 and sequenced on the GridION utilizing R10.4.1 flowcells. Phenotypic antimicrobial susceptibility testing was interpreted based on CLSI breakpoints. Sequencing data were processed with BugSeq, and AMR was assessed with BugSplit and ResFinder.

RESULTS

Fifty-six isolates were sequenced, including 53 related to the cluster of cases. All cluster isolates were identified as S. sonnei by sequencing, with global genotype 3.6.1.1.2 (CipR.MSM5), MLST 152 and PopPUNK cluster 3. Core genome MLST (cgMLST, examining 2513 loci) and reference-based MLST (refMLST, examining 4091 loci) both confirmed the clonality of the isolates. Cluster isolates were resistant to ampicillin (blaTEM-1), trimethoprim/sulfamethoxazole (dfA1, dfrA17; sul1, sul2), azithromycin (ermB, mphA) and ciprofloxacin (gyrA S83L, gyrA D87G, parC S80I). No genomic predictors of resistance to carbapenems were identified.

CONCLUSIONS

WGS with R10.4.1 enabled rapid sequencing and identification of an MDR S. sonnei community cluster. Genetic predictors of AMR were concordant with phenotypic antimicrobial susceptibility testing.

Gapless genome assembly and epigenetic profiles reveal gene regulation of whole-genome triplication in lettuce

BACKGROUND

Lettuce, an important member of the Asteraceae family, is a globally cultivated cash vegetable crop. With a highly complex genome (∼2.5 Gb; 2n = 18) rich in repeat sequences, current lettuce reference genomes exhibit thousands of gaps, impeding a comprehensive understanding of the lettuce genome.

FINDINGS

Here, we present a near-complete gapless reference genome for cutting lettuce with high transformability, using long-read PacBio HiFi and Nanopore sequencing data. In comparison to stem lettuce genome, we identify 127,681 structural variations (SVs, present in 0.41 Gb of sequence), reflecting the divergence of leafy and stem lettuce. Interestingly, these SVs are related to transposons and DNA methylation states. Furthermore, we identify 4,612 whole-genome triplication genes exhibiting high expression levels associated with low DNA methylation levels and high N6-methyladenosine RNA modifications. DNA methylation changes are also associated with activation of genes involved in callus formation.

CONCLUSIONS

Our gapless lettuce genome assembly, an unprecedented achievement in the Asteraceae family, establishes a solid foundation for functional genomics, epigenomics, and crop breeding and sheds new light on understanding the complexity of gene regulation associated with the dynamics of DNA and RNA epigenetics in genome evolution.

High-quality metagenome assembly from long accurate reads with metaMDBG

We introduce metaMDBG, a metagenomics assembler for PacBio HiFi reads. MetaMDBG combines a de Bruijn graph assembly in a minimizer space with an iterative assembly over sequences of minimizers to address variations in genome coverage depth and an abundance-based filtering strategy to simplify strain complexity. For complex communities, we obtained up to twice as many high-quality circularized prokaryotic metagenome-assembled genomes as existing methods and had better recovery of viruses and plasmids.

Community-scale models of microbiomes: Articulating metabolic modelling and metagenome sequencing

Building models is essential for understanding the functions and dynamics of microbial communities. Metabolic models built on genome-scale metabolic network reconstructions (GENREs) are especially relevant as a means to decipher the complex interactions occurring among species. Model reconstruction increasingly relies on metagenomics, which permits direct characterisation of naturally occurring communities that may contain organisms that cannot be isolated or cultured. In this review, we provide an overview of the field of metabolic modelling and its increasing reliance on and synergy with metagenomics and bioinformatics. We survey the means of assigning functions and reconstructing metabolic networks from (meta-)genomes, and present the variety and mathematical fundamentals of metabolic models that foster the understanding of microbial dynamics. We emphasise the characterisation of interactions and the scaling of model construction to large communities, two important bottlenecks in the applicability of these models. We give an overview of the current state of the art in metagenome sequencing and bioinformatics analysis, focusing on the reconstruction of genomes in microbial communities. Metagenomics benefits tremendously from third-generation sequencing, and we discuss the opportunities of long-read sequencing, strain-level characterisation and eukaryotic metagenomics. We aim at providing algorithmic and mathematical support, together with tool and application resources, that permit bridging the gap between metagenomics and metabolic modelling.

High Arctic seawater and coastal soil microbiome co-occurrence and composition structure and their potential hydrocarbon biodegradation

The accelerated decline in Arctic sea-ice cover and duration is enabling the opening of Arctic marine passages and improving access to natural resources. The increasing accessibility to navigation and resource exploration and production brings risks of accidental hydrocarbon releases into Arctic waters, posing a major threat to Arctic marine ecosystems where oil may persist for many years, especially in beach sediment. The composition and response of the microbial community to oil contamination on Arctic beaches remain poorly understood. To address this, we analyzed microbial community structure and identified hydrocarbon degradation genes among the Northwest Passage intertidal beach sediments and shoreline seawater from five high Arctic beaches. Our results from 16S/18S rRNA genes, long-read metagenomes, and metagenome-assembled genomes reveal the composition and metabolic capabilities of the hydrocarbon microbial degrader community, as well as tight cross-habitat and cross-kingdom interactions dominated by lineages that are common and often dominant in the polar coastal habitat, but distinct from petroleum hydrocarbon-contaminated sites. In the polar beach sediment habitats, Granulosicoccus sp. and Cyclocasticus sp. were major potential hydrocarbon-degraders, and our metagenomes revealed a small proportion of microalgae and algal viruses possessing key hydrocarbon biodegradative genes. This research demonstrates that Arctic beach sediment and marine microbial communities possess the ability for hydrocarbon natural attenuation. The findings provide new insights into the viral and microalgal communities possessing hydrocarbon degradation genes and might represent an important contribution to the removal of hydrocarbons under harsh environmental conditions in a pristine, cold, and oil-free environment that is threatened by oil spills.

Colicins and T6SS-based competition systems enhance enterotoxigenic E. coli (ETEC) competitiveness

Diarrheal diseases are still a significant problem for humankind, causing approximately half a million deaths annually. To cause diarrhea, enteric bacterial pathogens must first colonize the gut, which is a niche occupied by the normal bacterial microbiota. Therefore, the ability of pathogenic bacteria to inhibit the growth of other bacteria can facilitate the colonization process. Although enterotoxigenic Escherichia coli (ETEC) is one of the major causative agents of diarrheal diseases, little is known about the competition systems found in and used by ETEC and how they contribute to the ability of ETEC to colonize a host. Here, we collected a set of 94 fully assembled ETEC genomes by performing whole-genome sequencing and mining the NCBI RefSeq database. Using this set, we performed a comprehensive search for delivered bacterial toxins and investigated how these toxins contribute to ETEC competitiveness in vitro. We found that type VI secretion systems (T6SS) were widespread among ETEC (n = 47). In addition, several closely related ETEC strains were found to encode Colicin Ia and T6SS (n = 8). These toxins provide ETEC competitive advantages during in vitro competition against other E. coli, suggesting that the role of T6SS as well as colicins in ETEC biology has until now been underappreciated.

The Oxford Nanopore MinION as a Versatile Technology for the Diagnosis and Characterization of Emerging Plant Viruses

The emergence of novel viral epidemics that could affect major crops represents a serious threat to global food security. The early and accurate identification of the causative viral agent is the most important step for a rapid and effective response to disease outbreaks. Over the last years, the Oxford Nanopore Technologies (ONT) MinION sequencer has been proposed as an effective diagnostic tool for the early detection and identification of emerging viruses in plants, providing many advantages compared with different high-throughput sequencing (HTS) technologies. Here, we provide a step-by-step protocol that we optimized to obtain the virome of "Lamon bean" plants (Phaseolus vulgaris L.), an agricultural product with Protected Geographical Indication (PGI) in North-East of Italy, which is frequently subjected to multiple infections caused by different RNA viruses. The conversion of viral RNA in ds-cDNA enabled the use of Genomic DNA Ligation Sequencing Kit and Native Barcoding DNA Kit, which have been originally developed for DNA sequencing. This allowed the simultaneous diagnosis of both DNA- and RNA-based pathogens, providing a more versatile alternative to the use of direct RNA and/or direct cDNA sequencing kits.

Freshwater Viral Metagenome Analyses Targeting dsDNA Viruses

Viral metagenomics is one of the most widely used approaches to study viral population genomics. With the recent development of bioinformatic tools, the number of molecular biological methods, programs, and software to analyze viral metagenome data have greatly increased. Here, we describe the basic analysis workflow along with bioinformatic tools that can be used to analyze viral metagenome data. Although this chapter assumes that the viral metagenome data are prepared from the freshwater samples and are subjected to dsDNA sequencing, the protocol can be applied and modified for other types of metagenome data collected from a variety of sources.

Andean soil-derived lignocellulolytic bacterial consortium as a source of novel taxa and putative plastic-active enzymes

An easy and straightforward way to engineer microbial environmental communities is by setting up liquid enrichment cultures containing a specific substrate as the sole source of carbon. Here, we analyzed twenty single-contig high-quality metagenome-assembled genomes (MAGs) retrieved from a microbial consortium (T6) that was selected by the dilution-to-stimulation approach using Andean soil as inoculum and lignocellulose as a selection pressure. Based on genomic metrics (e.g., average nucleotide and amino acid identities) and phylogenomic analyses, 15 out of 20 MAGs were found to represent novel bacterial species, with one of those (MAG_26) belonging to a novel genus closely related to Caenibius spp. (Sphingomonadaceae). Following the rules and requirements of the SeqCode, we propose the name Andeanibacterium colombiense gen. nov., sp. nov. for this taxon. A subsequent functional annotation of all MAGs revealed that MAG_7 (Pseudobacter hemicellulosilyticus sp. nov.) contains 20, 19 and 16 predicted genes from carbohydrate-active enzymes families GH43, GH2 and GH92, respectively. Its lignocellulolytic gene profile resembles that of MAG_2 (the most abundant member) and MAG_3858, both of which belong to the Sphingobacteriaceae family. Using a database that contains experimentally verified plastic-active enzymes (PAZymes), twenty-seven putative bacterial polyethylene terephthalate (PET)-active enzymes (i.e., alpha/beta-fold hydrolases) were detected in all MAGs. A maximum of five putative PETases were found in MAG_3858, and two PETases were found to be encoded by A. colombiense. In conclusion, we demonstrate that lignocellulose-enriched liquid cultures coupled with genome-resolved metagenomics are suitable approaches to unveil the hidden bacterial diversity and its polymer-degrading potential in Andean soil ecosystems.

First single-strain enrichments of Electrothrix cable bacteria, description of E. aestuarii sp. nov. and E. rattekaaiensis sp. nov., and proposal of a cable bacteria taxonomy following the rules of the SeqCode

Cable bacteria are electrically conductive, filamentous Desulfobulbaceae, which are morphologically, functionally, and phylogenetically distinct from the other members of this family. Cable bacteria have not been obtained in pure culture and were therefore previously described as candidate genera, Candidatus Electrothrix and Ca. Electronema; a representative of the latter is available as single-strain sediment enrichment. Here we present an improved workflow to obtain the first single-strain enrichments of Ca. Electrothrix and report their metagenome-assembled genomes (MAGs) and morphology. Based on these results and on previously published high-quality MAGs and morphological data of cable bacteria from both candidate genera, we propose to adopt the genus names Electrothrix and Electronema following the rules of the Code of Nomenclature of Prokaryotes Described from Sequence Data (SeqCode), with Electrothrix communis RBTS and Electronema aureum GSTS, respectively, as the nomenclatural types of the genera. Furthermore, based on average nucleotide identity (ANI) values < 95 % with any described species, we propose two of our three single-strain enrichment cultures as novel species of the genus Electrothrix, with the names E. aestuarii sp. nov. and E. rattekaaiensis sp. nov., according to the SeqCode.

Air monitoring by nanopore sequencing

While the air microbiome and its diversity are essential for human health and ecosystem resilience, comprehensive air microbial diversity monitoring has remained rare, so that little is known about the air microbiome's composition, distribution, or functionality. Here we show that nanopore sequencing-based metagenomics can robustly assess the air microbiome in combination with active air sampling through liquid impingement and tailored computational analysis. We provide fast and portable laboratory and computational approaches for air microbiome profiling, which we leverage to robustly assess the taxonomic composition of the core air microbiome of a controlled greenhouse environment and of a natural outdoor environment. We show that long-read sequencing can resolve species-level annotations and specific ecosystem functions through de novo metagenomic assemblies despite the low amount of fragmented DNA used as an input for nanopore sequencing. We then apply our pipeline to assess the diversity and variability of an urban air microbiome, using Barcelona, Spain, as an example; this randomized experiment gives first insights into the presence of highly stable location-specific air microbiomes within the city's boundaries, and showcases the robust microbial assessments that can be achieved through automatable, fast, and portable nanopore sequencing technology.

Putative novel hydrogen- and iron-oxidizing sheath-producing Zetaproteobacteria thrive at the Fåvne deep-sea hydrothermal vent field

IMPORTANCE

Knowledge on microbial iron oxidation is important for understanding the cycling of iron, carbon, nitrogen, nutrients, and metals. The current study yields important insights into the niche sharing, diversification, and Fe(III) oxyhydroxide morphology of Ghiorsea, an iron- and hydrogen-oxidizing Zetaproteobacteria representative belonging to Zetaproteobacteria operational taxonomic unit 9. The study proposes that Ghiorsea exhibits a more extensive morphology of Fe(III) oxyhydroxide than previously observed. Overall, the results increase our knowledge on potential drivers of Zetaproteobacteria diversity in iron microbial mats and can eventually be used to develop strategies for the cultivation of sheath-forming Zetaproteobacteria.

Phosphate-related genomic islands as drivers of environmental adaptation in the streamlined marine alphaproteobacterial HIMB59

IMPORTANCE

These results shed light on the evolutionary strategies of microbes with streamlined genomes to adapt and survive in the oligotrophic conditions that dominate the surface waters of the global ocean. At the individual level, these microbes have been subjected to evolutionary constraints that have led to a more efficient use of nutrients, removing non-essential genes named as "streamlining theory." However, at the population level, they conserve a highly diverse gene pool in flexible genomic islands resulting in polyclonal populations on the same genomic background as an evolutionary response to environmental pressures. Localization of these islands at equivalent positions in the genome facilitates horizontal transfer between clonal lineages. This high level of environmental genomic heterogeneity could explain their cosmopolitan distribution. In the case of the order HIMB59 within the class Alphaproteobacteria, two factors exert evolutionary pressure and determine this intraspecific diversity: phages and the concentration of P in the environment.

Indole Diketopiperazine Alkaloids from the Marine Sediment-Derived Fungus Aspergillus chevalieri against Pancreatic Ductal Adenocarcinoma

A new prenylated indole diketopiperazine alkaloid, rubrumline P (1), was isolated along with six more analogues and characterized from the fermentation culture of a marine sediment-derived fungus, Aspergillus chevalieri, collected at a depth of 15 m near the lighthouse in Dahab, Red Sea, Egypt. In the current study, a bioassay-guided fractionation allowed for the identification of an active fraction displaying significant cytotoxic activity against the human pancreatic adenocarcinoma cell line PANC-1 from the EtOAc extract of the investigated fungus compared to the standard paclitaxel. The structures of the isolated compounds from the active fraction were established using 1D/2D NMR spectroscopy and mass spectrometry, together with comparisons with the literature. The absolute configuration of the obtained indole diketopiperazines was established based on single-crystal X-ray diffraction analyses of rubrumline I (2) and comparisons of optical rotations and NMR data, as well as on biogenetic considerations. Genome sequencing indicated the formation of prenyltransferases, which was subsequently confirmed by the isolation of mono-, di-, tri-, and tetraprenylated compounds. Compounds rubrumline P (1) and neoechinulin D (4) confirmed preferential cytotoxic activity against PANC-1 cancer cells with IC50 values of 25.8 and 23.4 µM, respectively. Although the underlying mechanism-of-action remains elusive in this study, cell cycle analysis indicated a slight increase in the sub-G1 peak after treatment with compounds 1 and 4.

The undiscovered biosynthetic potential of the Greenland Ice Sheet microbiome

The Greenland Ice Sheet is a biome which is mainly microbially driven. Several different niches can be found within the glacial biome for those microbes able to withstand the harsh conditions, e.g., low temperatures, low nutrient conditions, high UV radiation in summer, and contrasting long and dark winters. Eukaryotic algae can form blooms during the summer on the ice surface, interacting with communities of bacteria, fungi, and viruses. Cryoconite holes and snow are also habitats with their own microbial community. Nevertheless, the microbiome of supraglacial habitats remains poorly studied, leading to a lack of representative genomes from these environments. Under-investigated extremophiles, like those living on the Greenland Ice Sheet, may provide an untapped reservoir of chemical diversity that is yet to be discovered. In this study, an inventory of the biosynthetic potential of these organisms is made, through cataloging the presence of biosynthetic gene clusters in their genomes. There were 133 high-quality metagenome-assembled genomes (MAGs) and 28 whole genomes of bacteria obtained from samples of the ice sheet surface, cryoconite, biofilm, and snow using culturing-dependent and -independent approaches. AntiSMASH and BiG-SCAPE were used to mine these genomes and subsequently analyze the resulting predicted gene clusters. Extensive sets of predicted Biosynthetic Gene Clusters (BGCs) were collected from the genome collection, with limited overlap between isolates and MAGs. Additionally, little overlap was found in the biosynthetic potential among different environments, suggesting specialization of organisms in specific habitats. The median number of BGCs per genome was significantly higher for the isolates compared to the MAGs. The most talented producers were found among Proteobacteria. We found evidence for the capacity of these microbes to produce antimicrobials, carotenoid pigments, siderophores, and osmoprotectants, indicating potential survival mechanisms to cope with extreme conditions. The majority of identified BGCs, including those in the most prevalent gene cluster families, have unknown functions, presenting a substantial potential for bioprospecting. This study underscores the diverse biosynthetic potential in Greenland Ice Sheet genomes, revealing insights into survival strategies and highlighting the need for further exploration and characterization of these untapped resources.

Hybrid RNA sequencing of broad bean wilt virus 2 from faba beans

IMPORTANCE

Globally, viral diseases impair the growth and vigor of cultivated crops such as grains, leading to a significant reduction in quality, marketability, and competitiveness. As an island nation, Australia has a distinct advantage in using its border to prevent the introduction of damaging viruses, which threaten the continental agricultural sector. However, breeding programs in Australia rely on imported seeds as new sources of genetic diversity. As such, it is critical to remain vigilant in identifying new and emerging viral pathogens, by ensuring the availability of accurate genomic diagnostic tools at the grain biosecurity border. High-throughput sequencing offers game-changing opportunities in biosecurity routine testing. Genomic results are more accurate and informative compared to traditional molecular methods or biological indexing. The present work contributes to strengthening accurate phytosanitary screening, to safeguard the Australian grains industry, and expedite germplasm release to the end users.