metaFlye: scalable long-read metagenome assembly using repeat graphs

Carbohydrates and carbohydrate degradation gene abundance and transcription in Atlantic waters of the Arctic

Carbohydrates are chemically and structurally diverse, represent a substantial fraction of marine organic matter and are key substrates for heterotrophic microbes. Studies on carbohydrate utilisation by marine microbes have been centred on phytoplankton blooms in temperate regions, while far less is known from high-latitude waters and during later seasonal stages. Here, we combine glycan microarrays and analytical chromatography with metagenomics and metatranscriptomics to show the spatial heterogeneity in glycan distribution and potential carbohydrate utilisation by microbes in Atlantic waters of the Arctic. The composition and abundance of monomers and glycan structures in POM varied with location and depth. Complex fucose-containing sulfated polysaccharides, known to accumulate in the ocean, were consistently detected, while the more labile β-1,3-glucan exhibited a patchy distribution. Through 'omics analysis, we identify variations in the abundance and transcription of carbohydrate degradation-related genes across samples at the community and population level. The populations contributing the most to transcription were taxonomically related to those known as primary responders and key carbohydrate degraders in temperate ecosystems, such as NS4 Marine Group and Formosa. The unique transcription profiles for these populations suggest distinct substrate utilisation potentials, with predicted glycan targets corresponding to those structurally identified in POM from the same sampling sites. By combining cutting-edge technologies and protocols, we provide insights into the carbohydrate component of the carbon cycle in the Arctic during late summer and present a high-quality dataset that will be of great value for future comparative analyses.

Dancing the Nanopore limbo - Nanopore metagenomics from small DNA quantities for bacterial genome reconstruction

BACKGROUND

While genome-resolved metagenomics has revolutionized our understanding of microbial and genetic diversity in environmental samples, assemblies of short-reads often result in incomplete and/or highly fragmented metagenome-assembled genomes (MAGs), hampering in-depth genomics. Although Nanopore sequencing has increasingly been used in microbial metagenomics as long reads greatly improve the assembly quality of MAGs, the recommended DNA quantity usually exceeds the recoverable amount of DNA of environmental samples. Here, we evaluated lower-than-recommended DNA quantities for Nanopore library preparation by determining sequencing quality, community composition, assembly quality and recovery of MAGs.

RESULTS

We generated 27 Nanopore metagenomes using the commercially available ZYMO mock community and varied the amount of input DNA from 1000 ng (the recommended minimum) down to 1 ng in eight steps. The quality of the generated reads remained stable across all input levels. The read mapping accuracy, which reflects how well the reads match a known reference genome, was consistently high across all libraries. The relative abundance of the species in the metagenomes was stable down to input levels of 50 ng. High-quality MAGs (> 95% completeness, ≤ 5% contamination) could be recovered from metagenomes down to 35 ng of input material. When combined with publicly available Illumina reads for the mock community, Nanopore reads from input quantities as low as 1 ng improved the quality of hybrid assemblies.

CONCLUSION

Our results show that the recommended DNA amount for Nanopore library preparation can be substantially reduced without any adverse effects to genome recovery and still bolster hybrid assemblies when combined with short-read data. We posit that the results presented herein will enable studies to improve genome recovery from low-biomass environments, enhancing microbiome understanding.

Genome sequence and characterisation of a freshwater photoarsenotroph, Cereibacter azotoformans strain ORIO, isolated from sediments capable of cyclic light-dark arsenic oxidation and reduction

A freshwater photosynthetic arsenite-oxidizing bacterium, Cereibacter azotoformans strain ORIO, was isolated from Owens River, CA, USA. The waters from Owens River are elevated in arsenic and serve as the headwaters to the Los Angeles Aqueduct. The complete genome sequence of strain ORIO is 4.8 Mb genome (68% G + C content) and comprises two chromosomes and six plasmids. Taxonomic analysis placed ORIO within the Cereibacter genus (formerly Rhodobacter). The ORIO genome contains arxB2 AB1 CD (encoding an arsenite oxidase), arxXSR (regulators) and several ars arsenic resistance genes all co-localised on a 136 kb plasmid, named pORIO3. Phylogenetic analysis of ArxA, the molybdenum-containing arsenite oxidase catalytic subunit, demonstrated photoarsenotrophy is likely to occur within members of the Alphaproteobacteria. ORIO is a mixotroph, oxidises arsenite to arsenate (As(V)) photoheterotrophically, and expresses arxA in cultures grown with arsenite. Further ecophysiology studies with Owens River sediment demonstrated the interconversion of arsenite and As(V) was dependent on light-dark cycling. arxA and arrA (As(V) respiratory reductase) genes were detected in the light-dark cycled sediment metagenomes suggesting syntrophic interactions among arsenotrophs. This work establishes C. azotoformans str. ORIO as a new model organism for studying photoarsenotrophy and light-dark arsenic biogeochemical cycling.

Diverse electron carriers drive syntrophic interactions in an enriched anaerobic acetate-oxidizing consortium

In many anoxic environments, syntrophic acetate oxidation (SAO) is a key pathway mediating the conversion of acetate into methane through obligate cross-feeding interactions between SAO bacteria (SAOB) and methanogenic archaea. The SAO pathway is particularly important in engineered environments such as anaerobic digestion (AD) systems operating at thermophilic temperatures and/or with high ammonia. Despite the widespread importance of SAOB to the stability of the AD process, little is known about their in situ physiologies due to typically low biomass yields and resistance to isolation. Here, we performed a long-term (300-day) continuous enrichment of a thermophilic (55 °C) SAO community from a municipal AD system using acetate as the sole carbon source. Over 80% of the enriched bioreactor metagenome belonged to a three-member consortium, including an acetate-oxidizing bacterium affiliated with DTU068 encoding for carbon dioxide, hydrogen, and formate production, along with two methanogenic archaea affiliated with Methanothermobacter_A. Stable isotope probing was coupled with metaproteogenomics to quantify carbon flux into each community member during acetate conversion and inform metabolic reconstruction and genome-scale modeling. This effort revealed that the two Methanothermobacter_A species differed in their preferred electron donors, with one possessing the ability to grow on formate and the other only consuming hydrogen. A thermodynamic analysis suggested that the presence of the formate-consuming methanogen broadened the environmental conditions where ATP production from SAO was favorable. Collectively, these results highlight how flexibility in electron partitioning during SAO likely governs community structure and fitness through thermodynamic-driven mutualism, shedding valuable insights into the metabolic underpinnings of this key functional group within methanogenic ecosystems.

Correcting modification-mediated errors in nanopore sequencing by nucleotide demodification and reference-based correction

The accuracy of Oxford Nanopore Technology (ONT) sequencing has significantly improved thanks to new flowcells, sequencing kits, and basecalling algorithms. However, novel modification types untrained in the basecalling models can seriously reduce the quality. Here we reports a set of ONT-sequenced genomes with unexpected low quality due to novel modification types. Demodification by whole-genome amplification significantly improved the quality but lost the epigenome. We also developed a reference-based method, Modpolish, for correcting modification-mediated errors while retaining the epigenome when a sufficient number of closely-related genomes is publicly available (default: top 20 genomes with at least 95% identity). Modpolish not only significantly improved the quality of in-house sequenced genomes but also public datasets sequenced by R9.4 and R10.4 (simplex). Our results suggested that novel modifications are prone to ONT systematic errors. Nevertheless, these errors are correctable by nucleotide demodification or Modpolish without prior knowledge of modifications.

Energetically exploiting lignocellulose-rich residues in anaerobic digestion technologies: from bioreactors to proteogenomics

The biogas produced through anaerobic digestion (AD) of renewable feedstocks is one of the promising alternatives to replace fossil-derived energy. Even though lignocellulosic biomass is the most abundant biomass on earth, only a small fraction is being used towards resources recovery, leaving a great potential unexploited. In this study, the combination of state-of-art genomic techniques and engineered systems were used to further advance the knowledge on biogas production from lignocellulosic-rich residues and the microbiome involved in the anaerobic digestion hereof. A long-term adapted anaerobic microbiome capable of degrading wheat straw as the sole substrate was investigated using protein stable isotope probing (protein-SIP). The results indicated that a diverse microbial community, primarily composed of Firmicutes and Methanogens, played crucial roles in cellulose degradation and methane production. Notably, Defluviitoga tunisiensis, Syntrophothermus lipocalidus, and Pelobacter carbinolicus were identified as direct metabolizers of cellulose, while Dehalobacterium assimilated labelled carbon through cross-feeding. This study provides direct evidence of primary cellulose degraders and sheds light on their genomic composition. By harnessing the potential of lignocellulosic biomass and understanding the microbial communities involved, we can promote sustainable biogas production, contributing to energy security and environmental preservation.

Diurnal cycles drive rhythmic physiology and promote survival in facultative phototrophic bacteria

Bacteria have evolved many strategies to spare energy when nutrients become scarce. One widespread such strategy is facultative phototrophy, which helps heterotrophs supplement their energy supply using light. Our knowledge of the impact that such behaviors have on bacterial fitness and physiology is, however, still limited. Here, we study how a representative of the genus Porphyrobacter, in which aerobic anoxygenic phototrophy is ancestral, responds to different light regimes under nutrient limitation. We show that bacterial survival in stationary phase relies on functional reaction centers and varies depending on the light regime. Under dark-light alternance, our bacterial model presents a diphasic life history dependent on phototrophy: during dark phases, the cells inhibit DNA replication and part of the population lyses and releases nutrients, while subsequent light phases allow for the recovery and renewed growth of the surviving cells. We correlate these cyclic variations with a pervasive pattern of rhythmic transcription which reflects global changes in diurnal metabolic activity. Finally, we demonstrate that, compared to either a phototrophy mutant or a bacteriochlorophyll a overproducer, the wild type strain is better adapted to natural environments, where regular dark-light cycles are interspersed with additional accidental dark episodes. Overall, our results highlight the importance of light-induced biological rhythms in a new model of aerobic anoxygenic phototroph representative of an ecologically important group of environmental bacteria.

Similar but different: Characterization of dddD gene-mediated DMSP metabolism among coral-associated Endozoicomonas

Endozoicomonas are often predominant bacteria and prominently important in coral health. Their role in dimethylsulfoniopropionate (DMSP) degradation has been a subject of discussion for over a decade. A previous study found that Endozoicomonas degraded DMSP through the dddD pathway. This process releases dimethyl sulfide, which is vital for corals coping with thermal stress. However, little is known about the related gene regulation and metabolic abilities of DMSP metabolism in Endozoicomonadaceae. In this study, we isolated a novel Endozoicomonas DMSP degrader and observed a distinct DMSP metabolic trend in two phylogenetically close dddD-harboring Endozoicomonas species, confirmed genetically by comparative transcriptomic profiling and visualization of the change of DMSP stable isotopes in bacterial cells using nanoscale secondary ion spectrometry. Furthermore, we found that DMSP cleavage enzymes are ubiquitous in coral Endozoicomonas with a preference for having DddD lyase. We speculate that harboring DMSP degrading genes enables Endozoicomonas to successfully colonize various coral species across the globe.

Viral Diversity in Benthic Abyssal Ecosystems: Ecological and Methodological Considerations

Viruses are the most abundant 'biological entities' in the world's oceans. However, technical and methodological constraints limit our understanding of their diversity, particularly in benthic abyssal ecosystems (>4000 m depth). To verify advantages and limitations of analyzing virome DNA subjected either to random amplification or unamplified, we applied shotgun sequencing-by-synthesis to two sample pairs obtained from benthic abyssal sites located in the North-eastern Atlantic Ocean at ca. 4700 m depth. One amplified DNA sample was also subjected to single-molecule long-read sequencing for comparative purposes. Overall, we identified 24,828 viral Operational Taxonomic Units (vOTUs), belonging to 22 viral families. Viral reads were more abundant in the amplified DNA samples (38.5-49.9%) compared to the unamplified ones (4.4-5.8%), with the latter showing a greater viral diversity and 11-16% of dsDNA viruses almost undetectable in the amplified samples. From a procedural point of view, the viromes obtained by direct sequencing (without amplification step) provided a broader overview of both ss and dsDNA viral diversity. Nevertheless, our results suggest that the contextual use of random amplification of the same sample and long-read technology can improve the assessment of viral assemblages by reducing off-target reads.

Hybrid assemblies of microbiome Blastocystis protists reveal evolutionary diversification reflecting host ecology

The most prevalent microbial eukaryote in the human gut is Blastocystis , an obligate commensal protist also common in many other vertebrates. Blastocystis is descended from free-living stramenopile ancestors; how it has adapted to thrive within humans and a wide range of hosts is unclear. Here, we cultivated six Blastocystis strains spanning the diversity of the genus and generated highly contiguous, annotated genomes with long-read DNA-seq, Hi-C, and RNA-seq. Comparative genomics between these strains and two closely related stramenopiles with different lifestyles, the lizard gut symbiont Proteromonas lacertae and the free-living marine flagellate Cafeteria burkhardae , reveal the evolutionary history of the Blastocystis genus. We find substantial gene content variability between Blastocystis strains. Blastocystis isolated from an herbivorous tortoise has many plant carbohydrate metabolizing enzymes, some horizontally acquired from bacteria, likely reflecting fermentation within the host gut. In contrast, human- isolated Blastocystis have gained many heat shock proteins, and we find numerous subtype- specific expansions of host-interfacing genes, including cell adhesion and cell surface glycan genes. In addition, we observe that human-isolated Blastocystis have substantial changes in gene structure, including shortened introns and intergenic regions, as well as genes lacking canonical termination codons. Finally, our data indicate that the common ancestor of Blastocystis lost nearly all ancestral genes for heterokont flagella morphology, including cilia proteins, microtubule motor proteins, and ion channel proteins. Together, these findings underscore the huge functional variability within the Blastocystis genus and provide candidate genes for the adaptations these lineages have undergone to thrive in the gut microbiomes of diverse vertebrates.

Diversity, Taxonomic Novelty, and Encoded Functions of Salar de Ascotán Microbiota, as Revealed by Metagenome-Assembled Genomes

Salar de Ascotán is a high-altitude arsenic-rich salt flat exposed to high ultraviolet radiation in the Atacama Desert, Chile. It hosts unique endemic flora and fauna and is an essential habitat for migratory birds, making it an important site for conservation and protection. However, there is limited information on the resident microbiota's diversity, genomic features, metabolic potential, and molecular mechanisms that enable it to thrive in this extreme environment. We used long- and short-read metagenomics to investigate the microbial communities in Ascotán's water, sediment, and soil. Bacteria predominated, mainly Pseudomonadota, Acidobacteriota, and Bacteroidota, with a remarkable diversity of archaea in the soil. Following hybrid assembly, we recovered high-quality bacterial (101) and archaeal (6) metagenome-assembled genomes (MAGs), including representatives of two putative novel families of Patescibacteria and Pseudomonadota and two novel orders from the archaeal classes Halobacteriota and Thermoplasmata. We found different metabolic capabilities across distinct lineages and a widespread presence of genes related to stress response, DNA repair, and resistance to arsenic and other metals. These results highlight the remarkable diversity and taxonomic novelty of the Salar de Ascotán microbiota and its rich functional repertoire, making it able to resist different harsh conditions. The highly complete MAGs described here could serve future studies and bioprospection efforts focused on salt flat extremophiles, and contribute to enriching databases with microbial genome data from underrepresented regions of our planet.

Functional enrichment of integrons: Facilitators of antimicrobial resistance and niche adaptation

Integrons are genetic elements, found among diverse bacteria and archaea, that capture and rearrange gene cassettes to rapidly generate genetic diversity and drive adaptation. Despite their broad taxonomic and geographic prevalence, and their role in microbial adaptation, the functions of gene cassettes remain poorly characterized. Here, using a combination of bioinformatic and experimental analyses, we examined the functional diversity of gene cassettes from different environments. We find that cassettes encode diverse antimicrobial resistance (AMR) determinants, including those conferring resistance to antibiotics currently in the developmental pipeline. Further, we find a subset of cassette functions is universally enriched relative to their broader metagenomes. These are largely involved in (a)biotic interactions, including AMR, phage defense, virulence, biodegradation, and stress tolerance. The remainder of functions are sample-specific, suggesting that they confer localised functions relevant to their microenvironment. Together, they comprise functional profiles different from bulk metagenomes, representing niche-adaptive components of the prokaryotic pangenome.

Genome sequence of a symbiotic cyanobacterium from the flowering plant Gunnera tinctoria

Metagenomic analysis of the symbiotic cyanobacteria colonies within Gunnera tinctoria stems revealed a new strain of Nostoc. Here, we report its genome sequence.

Pseudomonas kielensis str. Ze23jcel16 complete genome obtained through R10.4.1 Nanopore Flow cell chemistry

Here, we report the complete genome of Pseudomonas kielensis str. Ze23jcel16 isolated from a freshwater sample. The high-quality chromosome was obtained employing R10.4.1 Nanopore Flow cell chemistry and was assembled as a circular element at 45× coverage, a length of 5.8 Mbp, and a G+C content of 61.15%.

Mobilization of the blaKPC-14 gene among heterogenous plasmids in extensively drug-resistant hypervirulent Klebsiella pneumoniae

Introduction

Ceftazidime/avibactam (CZA) is an effective alternative for the treatment of infections caused by KPC-producing carbapenem-resistant Klebsiella pneumoniae (CRKP). However, KPC variants with CZA resistance have been observed in clinical isolates, further limiting the treatment options of clinical use.

Methods

In this study, we isolated three KPC-14-producing CRKP from two patients in intensive care units without CZA therapy. The antimicrobial susceptibility was determined using the broth microdilution method. Three CRKP were subjected to whole-genome sequencing to analyze the phylogenetic relatedness and the carriage of antimicrobial resistance genes and virulence factors. Long-read sequencing was also performed to obtain the complete sequences of the plasmids. The horizontal transfer of the blaKPC-14 gene was evaluated by conjugation experiments.

Results

Three CRKP displayed resistance or reduced susceptibility to ceftazidime/avibactam, colistin, and tigecycline. Single-nucleotide polymorphism (SNP) analysis demonstrated the close phylogenetic distance between these strains. A highly similar IncFII/IncR plasmid encoding blaKPC-14 was shared by three CRKP, with blaKPC-14 located in an NTEKPC-Ib element with the core region of ISKpn27- blaKPC-14-ISKpn6. This structure containing blaKPC-14 was also observed in another tet(A)-carrying plasmid that belonged to an unknown Inc-type in two out of three isolates. The horizontal transferability of these integrated plasmids to Escherichia coli EC600 was confirmed by the cotransmission of tet(A) and blaKPC-14 genes, but the single transfer of blaKPC-14 on the IncFII/IncR plasmid failed. Three CRKP expressed yersiniabactin and carried a hypervirulence plasmid encoding rmpA2 and aerobactin-related genes, and were thus classified as carbapenem-resistant hypervirulent K. pneumoniae (hvKP).

Discussion

In this study, we reported the evolution of a mosaic plasmid encoding the blaKPC-14 gene via mobile elements in extensively drug-resistant hvKP. The blaKPC-14 gene is prone to integrate into other conjugative plasmids via the NTEKPC-Ib element, further facilitating the spread of ceftazidime/avibactam resistance.

Role of Relebactam in the Antibiotic Resistance Acquisition in Pseudomonas aeruginosa: In Vitro Study

BACKGROUND

Pseudomonas aeruginosa shows resistance to several antibiotics and often develops such resistance during patient treatment.

OBJECTIVE

Develop an in vitro model, using clinical isolates of P. aeruginosa, to compare the ability of the imipenem and imipenem/relebactam to generate resistant mutants to imipenem and to other antibiotics. Perform a genotypic analysis to detect how the selective pressure changes their genomes.

METHODS

The antibiotics resistance was studied by microdilution assays and e-test, and the genotypic study was performed by NGS.

RESULTS

The isolates acquired resistance to imipenem in an average of 6 days, and to imipenem/relebactam in 12 days (p value = 0.004). After 30 days of exposure, 75% of the isolates reached a MIC > 64 mg/L for imipenem and 37.5% for imipenem/relebactam (p value = 0.077). The 37.5% and the 12.5% imipenem/relebactam mutants developed resistance to piperacillin/tazobactam and ceftazidime, respectively, while the 87.5% and 37.5% of the imipenem mutants showed resistance to these drugs (p value = 0.003, p value = 0.015). The main biological processes altered by the SNPs were the glycosylation pathway, transcriptional regulation, histidine kinase response, porins, and efflux pumps.

DISCUSSION

The addition of relebactam delays the generation of resistance to imipenem and limits the cross-resistance to other beta-lactams. The clinical relevance of this phenomenon, which has the limitation that it has been performed in vitro, should be evaluated by stewardship programs in clinical practice, as it could be useful in controlling multi-drug resistance in P. aeruginosa.

A novel culture-enriched metagenomic sequencing strategy effectively guarantee the microbial safety of drinking water by uncovering the low abundance pathogens

Assessing the presence of waterborne pathogens and antibiotic resistance genes (ARGs) is crucial for managing the environmental quality of drinking water sources. However, detecting low abundance pathogens in such settings is challenging. In this study, a workflow was developed to enrich for broad spectrum pathogens from drinking water samples. A mock community was used to evaluate the effectiveness of various enrichment broths in detecting low-abundance pathogens. Monthly metagenomic surveillance was conducted in a drinking water source from May to September 2021, and water samples were subjected to five enrichment procedures for 6 h to recover the majority of waterborne bacterial pathogens. Oxford Nanopore Technology (ONT) was used for metagenomic sequencing of enriched samples to obtain high-quality pathogen genomes. The results showed that selective enrichment significantly increased the proportions of targeted bacterial pathogens. Compared to direct metagenomic sequencing of untreated water samples, targeted enrichment followed by ONT sequencing significantly improved the detection of waterborne pathogens and the quality of metagenome-assembled genomes (MAGs). Eighty-six high-quality MAGs, including 70 pathogen MAGs, were obtained from ONT sequencing, while only 12 MAGs representing 10 species were obtained from direct metagenomic sequencing of untreated water samples. In addition, ONT sequencing improved the recovery of mobile genetic elements and the accuracy of phylogenetic analysis. This study highlights the urgent need for efficient methodologies to detect and manage microbial risks in drinking water sources. The developed workflow provides a cost-effective approach for environmental management of drinking water sources with microbial risks. The study also uncovered pathogens that were not detected by traditional methods, thereby advancing microbial risk management of drinking water sources.

Candidatus Kirkpatrickella diaphorinae gen. nov., sp. nov., an uncultured endosymbiont identified in a population of Diaphorina citri from Hawaii

Diaphorina citri is the hemipteran pest and vector of a devastating bacterial pathogen of citrus worldwide. In addition to the two core bacterial endosymbionts of D. citri, Candidatus Carsonella ruddii and Candidatus Profftella armatura, the genome of a novel endosymbiont and as of yet undescribed microbe was discovered in a Hawaiian D. citri population through deep sequencing of multiple D. citri populations. Found to be closely related to the genus Asaia in the family Acetobacteraceae by 16S rRNA gene sequence analysis, it forms a sister clade along with other insect-associated 16S rRNA gene sequences from uncultured bacterium found associated with Aedes koreicus and Sogatella furcifera. Multilocus sequence analysis confirmed the phylogenetic placement sister to the Asaia clade. Despite the culturable Asaia clade being the closest phylogenetic neighbour, attempts to culture this newly identified bacterial endosymbiont were unsuccessful. On the basis of these distinct genetic differences, the novel endosymbiont is proposed to be classified into a candidate genus and species 'Candidatus Kirkpatrickella diaphorinae'. The full genome was deposited in GenBank (accession number CP107052; prokaryotic 16S rRNA OP600170).

Division of labor within psyllids: metagenomics reveals an ancient dual endosymbiosis with metabolic complementarity in the genus Cacopsylla

IMPORTANCE

Heritable beneficial bacterial endosymbionts have been crucial for the evolutionary success of numerous insects by enabling the exploitation of nutritionally limited food sources. Herein, we describe a previously unknown dual endosymbiosis in the psyllid genus Cacopsylla, consisting of the primary endosymbiont "Candidatus Carsonella ruddii" and a co-occurring Enterobacteriaceae bacterium for which we propose the name "Candidatus Psyllophila symbiotica." Its localization within the bacteriome and its small genome size confirm that Psyllophila is a co-primary endosymbiont widespread within the genus Cacopsylla. Despite its highly eroded genome, Psyllophila perfectly complements the tryptophan biosynthesis pathway that is incomplete in the co-occurring Carsonella. Moreover, the genome of Psyllophila is almost as small as Carsonella's, suggesting an ancient dual endosymbiosis that has now reached a precarious stage where any additional gene loss would make the system collapse. Hence, our results shed light on the dynamic interactions of psyllids and their endosymbionts over evolutionary time.

Metagenomics uncovers dietary adaptations for chitin digestion in the gut microbiota of convergent myrmecophagous mammals

IMPORTANCE

Myrmecophagous mammals are specialized in the consumption of ants and/or termites. They do not share a direct common ancestor and evolved convergently in five distinct placental orders raising questions about the underlying adaptive mechanisms involved and the relative contribution of natural selection and phylogenetic constraints. Understanding how these species digest their prey can help answer these questions. More specifically, the role of their gut microbial symbionts in the digestion of the insect chitinous exoskeleton has not been investigated in all myrmecophagous orders. We generated 29 new gut metagenomes from nine myrmecophagous species to reconstruct more than 300 bacterial genomes in which we identified chitin-degrading enzymes. Studying the distribution of these chitinolytic bacteria among hosts revealed both shared and specific bacteria between ant-eating species. Overall, our results highlight the potential role of gut symbionts in the convergent dietary adaptation of myrmecophagous mammals and the evolutionary mechanisms shaping their gut microbiota.

A comparison of Oxford nanopore library strategies for bacterial genomics

BACKGROUND

Oxford nanopore Technologies (ONT) provides three main library preparation strategies to sequence bacterial genomes. These include tagmentation (TAG), ligation (LIG) and amplification (PCR). Despite ONT's recommendations, making an informed decision for preparation choice remains difficult without a side-by-side comparison. Here, we sequenced 12 bacterial strains to examine the overall output of these strategies, including sequencing noise, barcoding efficiency and assembly quality based on mapping to curated genomes established herein.

RESULTS

Average read length ranged closely for TAG and LIG (> 5,000 bp), while being drastically smaller for PCR (< 1,100 bp). LIG produced the largest output with 33.62 Gbp vs. 11.72 Gbp for TAG and 4.79 Gbp for PCR. PCR produced the most sequencing noise with only 22.7% of reads mappable to the curated genomes, vs. 92.9% for LIG and 87.3% for TAG. Output per channel was most homogenous in LIG and most variable in PCR, while intermediate in TAG. Artifactual tandem content was most abundant in PCR (22.5%) and least in LIG and TAG (0.9% and 2.2%). Basecalling and demultiplexing of barcoded libraries resulted in ~ 20% data loss as unclassified reads and 1.5% read leakage.

CONCLUSION

The output of LIG was best (low noise, high read numbers of long lengths), intermediate in TAG (some noise, moderate read numbers of long lengths) and less desirable in PCR (high noise, high read numbers of short lengths). Overall, users should not accept assembly results at face value without careful replicon verification, including the detection of plasmids assembled from leaked reads.

Complete mitochondrial genome and draft chloroplastic genome of Haslea ostrearia (Simonsen 1974)

The first completed, circular mitochondrial genome and the first draft, linear chloroplastic genome of the blue diatom Haslea ostrearia (Simonsen 1974, Naviculaceae, Bacillariophyceae) were assembled from Illumina and PacBio sequencing. The mitochondrial genome was composed of 38,696 bases and contained 64 genes, including 31 protein-coding genes (CDS), 2 ribosomal RNA (rRNA) genes and 23 transfer RNA (tRNA) genes. For the chloroplast, the genome was composed of 130,200 bases with 169 genes (131 CDS, 6 rRNA genes, 31 tRNA genes, and 1 transfer messenger RNA gene). Phylogenetic trees, using the maximum-likehood method and partial genes currently available for Haslea ostrearia and other diatom species, suggested the proximity of all the Haslea ostrearia strains/isolates and the possibility of using these genomes as future references.

Transcriptomics of Lactobacillus paracasei: metabolism patterns and cellular responses under high-density culture conditions

Lactobacillus paracasei has significant potential for development and application in the environmental field, particularly in addressing malodor pollution. This study aims to investigate the cellular response of L. paracasei B1 under high-density culture conditions. The selected strain has previously shown effective deodorizing and bacteriostatic abilities. Transcriptomics techniques are employed to dissect the nutrient metabolism pattern of L. paracasei B1 and its response mechanism under environmental stress. The study characterizes the functions of key differentially expressed genes during growth before and after optimizing the culture conditions. The optimization of fermentation culture conditions provides a suitable growth environment for L. paracasei B1, inducing an enhancement of its phosphotransferase system for sugar source uptake and maintaining high levels of glycolysis and pyruvate metabolism. Consequently, the strain is able to grow and multiply rapidly. Under acid stress conditions, glycolysis and pyruvate metabolism are inhibited, and L. paracasei B1 generates additional energy through aerobic respiration to meet the energy demand. The two-component system and quorum sensing play roles in the response and regulation of L. paracasei B1 to adverse environments. The strain mitigates oxygen stress damage through glutathione metabolism, cysteine and methionine metabolism, base excision repair, and purine and pyrimidine metabolism. Additionally, the strain enhances lysine synthesis, the alanine, aspartate, and glutamate metabolic pathways, and relies on the ABC transport system to accumulate amino acid-compatible solutes to counteract acid stress and osmotic stress during pH regulation. These findings establish a theoretical basis for the further development and application of L. paracasei B1 for its productive properties.

PandaGUT provides new insights into bacterial diversity, function, and resistome landscapes with implications for conservation

BACKGROUND

The gut microbiota play important roles in host adaptation and evolution, but are understudied in natural population of wild mammals. To address host adaptive evolution and improve conservation efforts of threatened mammals from a metagenomic perspective, we established a high-quality gut microbiome catalog of the giant panda (pandaGUT) to resolve the microbiome diversity, functional, and resistome landscapes using approximately 7 Tbp of long- and short-read sequencing data from 439 stool samples.

RESULTS

The pandaGUT catalog comprises 820 metagenome-assembled genomes, including 40 complete closed genomes, and 64.5% of which belong to species that have not been previously reported, greatly expanding the coverage of most prokaryotic lineages. The catalog contains 2.37 million unique genes, with 74.8% possessing complete open read frames, facilitating future mining of microbial functional potential. We identified three microbial enterotypes across wild and captive panda populations characterized by Clostridium, Pseudomonas, and Escherichia, respectively. We found that wild pandas exhibited host genetic-specific microbial structures and functions, suggesting host-gut microbiota phylosymbiosis, while the captive cohorts encoded more multi-drug resistance genes.

CONCLUSIONS

Our study provides largely untapped resources for biochemical and biotechnological applications as well as potential intervention avenues via the rational manipulation of microbial diversity and reducing antibiotic usage for future conservation management of wildlife. Video Abstract.

MetaCC allows scalable and integrative analyses of both long-read and short-read metagenomic Hi-C data

Metagenomic Hi-C (metaHi-C) can identify contig-to-contig relationships with respect to their proximity within the same physical cell. Shotgun libraries in metaHi-C experiments can be constructed by next-generation sequencing (short-read metaHi-C) or more recent third-generation sequencing (long-read metaHi-C). However, all existing metaHi-C analysis methods are developed and benchmarked on short-read metaHi-C datasets and there exists much room for improvement in terms of more scalable and stable analyses, especially for long-read metaHi-C data. Here we report MetaCC, an efficient and integrative framework for analyzing both short-read and long-read metaHi-C datasets. MetaCC outperforms existing methods on normalization and binning. In particular, the MetaCC normalization module, named NormCC, is more than 3000 times faster than the current state-of-the-art method HiCzin on a complex wastewater dataset. When applied to one sheep gut long-read metaHi-C dataset, MetaCC binning module can retrieve 709 high-quality genomes with the largest species diversity using one single sample, including an expansion of five uncultured members from the order Erysipelotrichales, and is the only binner that can recover the genome of one important species Bacteroides vulgatus. Further plasmid analyses reveal that MetaCC binning is able to capture multi-copy plasmids.

Chlorine disinfection modifies the microbiome, resistome and mobilome of hospital wastewater - A nanopore long-read metagenomic approach

The aim of the present study was to analyze changes in the microbiome, resistome, and mobilome of hospital wastewater (HWW) induced by disinfection with chlorine compounds. Changes in bacterial communities and specific antibiotic resistance genes (ARGs) in HWW were determined with the use of a nanopore long-read metagenomic approach. The main hosts of ARGs in HWW were identified, and the mobility of resistance mechanisms was analyzed. Special attention was paid to the prevalence of critical-priority pathogens in the HWW microbiome, which pose the greatest threat to human health. The results of this study indicate that chlorine disinfection of HWW can induce significant changes in the structure of the total bacterial population and antibiotic resistant bacteria (ARB) communities, and that it can modify the resistome and mobilome of HWW. Disinfection favored the selection of ARGs, decreased their prevalence in HWW, while increasing their diversity. The mobility of the HWW resistome increased after disinfection. Disinfection led to the emergence of new drug resistance mechanisms in previously sensitive bacterial taxa. In conclusion, this study demonstrated that HWW disinfected with low (sublethal) concentrations of free chlorine significantly contributes to the mobility and transfer of drug resistance mechanisms (including critical mechanisms) between bacteria (including pathogens).

Biosynthetic gene cluster synteny: Orthologous polyketide synthases in Hypogymnia physodes, Hypogymnia tubulosa, and Parmelia sulcata

Lichens are symbiotic associations consisting of a photobiont (algae or cyanobacteria) and a mycobiont (fungus), which together generate a variety of unique secondary metabolites. To access this biosynthetic potential for biotechnological applications, deeper insights into the biosynthetic pathways and corresponding gene clusters are necessary. Here, we provide a comparative view of the biosynthetic gene clusters of three lichen mycobionts derived from Hypogymnia physodes, Hypogymnia tubulosa, and Parmelia sulcata. In addition, we present a high-quality PacBio metagenome of Parmelia sulcata, from which we extracted the mycobiont bin containing 214 biosynthetic gene clusters. Most biosynthetic gene clusters in these genomes were associated with T1PKSs, followed by NRPSs and terpenes. This study focused on biosynthetic gene clusters related to polyketide synthesis. Based on ketosynthase homology, we identified nine highly syntenic clusters present in all three species. Among the four clusters belonging to nonreducing PKSs, two are putatively linked to lichen substances derived from orsellinic acid (orcinol depsides and depsidones, e.g., lecanoric acid, physodic acid, lobaric acid), one to compounds derived from methylated forms of orsellinic acid (beta orcinol depsides, e.g., atranorin), and one to melanins. Five clusters with orthologs in all three species are linked to reducing PKSs. Our study contributes to sorting and dereplicating the vast PKS diversity found in lichenized fungi. High-quality sequences of biosynthetic gene clusters of these three common species provide a foundation for further exploration into biotechnological applications and the molecular evolution of lichen substances.

Whole genome sequence of the deep-sea sponge Geodia barretti (Metazoa, Porifera, Demospongiae)

Sponges are among the earliest branching extant animals. As such, genetic data from this group are valuable for understanding the evolution of various traits and processes in other animals. However, like many marine organisms, they are notoriously difficult to sequence, and hence, genomic data are scarce. Here, we present the draft genome assembly for the North Atlantic deep-sea high microbial abundance species Geodia barretti Bowerbank 1858, from a single individual collected on the West Coast of Sweden. The nuclear genome assembly has 4,535 scaffolds, an N50 of 48,447 bp and a total length of 144 Mb; the mitochondrial genome is 17,996 bp long. BUSCO completeness was 71.5%. The genome was annotated using a combination of ab initio and evidence-based methods finding 31,884 protein-coding genes.

Comparison of de novo assembly using long-read shotgun metagenomic sequencing of viruses in fecal and serum samples from marine mammals

Introduction

Viral diseases of marine mammals are difficult to study, and this has led to a limited knowledge on emerging known and unknown viruses which are ongoing threats to animal health. Viruses are the leading cause of infectious disease-induced mass mortality events among marine mammals.

Methods

In this study, we performed viral metagenomics in stool and serum samples from California sea lions (Zalophus californianus) and bottlenose dolphins (Tursiops truncates) using long-read nanopore sequencing. Two widely used long-read de novo assemblers, Canu and Metaflye, were evaluated to assemble viral metagenomic sequencing reads from marine mammals.

Results

Both Metaflye and Canu assembled similar viral contigs of vertebrates, such as Parvoviridae, and Poxviridae. Metaflye assembled viral contigs that aligned with one viral family that was not reproduced by Canu, while Canu assembled viral contigs that aligned with seven viral families that was not reproduced by Metaflye. Only Canu assembled viral contigs from dolphin and sea lion fecal samples that matched both protein and nucleotide RefSeq viral databases using BLASTx and BLASTn for Anelloviridae, Parvoviridae and Circoviridae families. Viral contigs assembled with Canu aligned with torque teno viruses and anelloviruses from vertebrate hosts. Viruses associated with invertebrate hosts including densoviruses, Ambidensovirus, and various Circoviridae isolates were also aligned. Some of the invertebrate and vertebrate viruses reported here are known to potentially cause mortality events and/or disease in different seals, sea stars, fish, and bivalve species.

Discussion

Canu performed better by producing the most viral contigs as compared to Metaflye with assemblies aligning to both protein and nucleotide databases. This study suggests that marine mammals can be used as important sentinels to surveil marine viruses that can potentially cause diseases in vertebrate and invertebrate hosts.

Unveiling metabolic pathways involved in the extreme desiccation tolerance of an Atacama cyanobacterium

Gloeocapsopsis dulcis strain AAB1 is an extremely xerotolerant cyanobacterium isolated from the Atacama Desert (i.e., the driest and oldest desert on Earth) that holds astrobiological significance due to its ability to biosynthesize compatible solutes at ultra-low water activities. We sequenced and assembled the G. dulcis genome de novo using a combination of long- and short-read sequencing, which resulted in high-quality consensus sequences of the chromosome and two plasmids. We leveraged the G. dulcis genome to generate a genome-scale metabolic model (iGd895) to simulate growth in silico. iGd895 represents, to our knowledge, the first genome-scale metabolic reconstruction developed for an extremely xerotolerant cyanobacterium. The model's predictive capability was assessed by comparing the in silico growth rate with in vitro growth rates of G. dulcis, in addition to the synthesis of trehalose. iGd895 allowed us to explore simulations of key metabolic processes such as essential pathways for water-stress tolerance, and significant alterations to reaction flux distribution and metabolic network reorganization resulting from water limitation. Our study provides insights into the potential metabolic strategies employed by G. dulcis, emphasizing the crucial roles of compatible solutes, metabolic water, energy conservation, and the precise regulation of reaction rates in their adaptation to water stress.

Complete genome sequence of mcr-9 containing Leclercia adecarboxylata

Here, we provide the genome sequence of a Leclercia adecarboxylata isolated from a screen of an environmental bacterial isolate library for resistance to the plant flavonoid berberine. We detected the colistin resistance gene mcr-9, located on an IncFII(pECLA) plasmid.

Illuminating the oral microbiome and its host interactions: recent advancements in omics and bioinformatics technologies in the context of oral microbiome research

The oral microbiota has an enormous impact on human health, with oral dysbiosis now linked to many oral and systemic diseases. Recent advancements in sequencing, mass spectrometry, bioinformatics, computational biology, and machine learning are revolutionizing oral microbiome research, enabling analysis at an unprecedented scale and level of resolution using omics approaches. This review contains a comprehensive perspective of the current state-of-the-art tools available to perform genomics, metagenomics, phylogenomics, pangenomics, transcriptomics, proteomics, metabolomics, lipidomics, and multi-omics analysis on (all) microbiomes, and then provides examples of how the techniques have been applied to research of the oral microbiome, specifically. Key findings of these studies and remaining challenges for the field are highlighted. Although the methods discussed here are placed in the context of their contributions to oral microbiome research specifically, they are pertinent to the study of any microbiome, and the intended audience of this includes researchers would simply like to get an introduction to microbial omics and/or an update on the latest omics methods. Continued research of the oral microbiota using omics approaches is crucial and will lead to dramatic improvements in human health, longevity, and quality of life.

Masculinizer gene controls male sex determination in the codling moth, Cydia pomonella

The molecular mechanisms of sex determination in moths and butterflies (Lepidoptera) with female heterogamety (WZ/ZZ) are poorly understood, except in the silkworm Bombyx mori. However, the Masculinizer (Masc) gene that controls male development and dosage compensation in B. mori, appears to be conserved in Lepidoptera, as its masculinizing function was recently confirmed in several moth species. In this work, we investigated the role of the Masc gene in sex determination of the codling moth Cydia pomonella (Tortricidae), a globally important pest of pome fruits and walnuts. The gene structure of the C. pomonella Masc ortholog, CpMasc, is similar to B. mori Masc. However, unlike B. mori, we identified 14 splice variants of CpMasc in the available transcriptomes. Subsequent screening for sex specificity and genetic variation using publicly available data and RT-PCR revealed three male-specific splice variants. Then qPCR analysis of these variants revealed sex-biased expression showing a peak only in early male embryos. Knockdown of CpMasc by RNAi during early embryogenesis resulted in a shift from male-to female-specific splicing of the C. pomonella doublesex (Cpdsx) gene, its downstream effector, in ZZ embryos, leading to a strongly female-biased sex ratio. These data clearly demonstrate that CpMasc functions as a masculinizing gene in the sex-determining cascade of C. pomonella. Our study also showed that CpMasc transcripts are provided maternally, as they were detected in unfertilized eggs after oviposition and in mature eggs dissected from virgin females. This finding is unique, as maternal provision of mRNA has rarely been studied in Lepidoptera.

A dynamic epibiont community associated with the bone-eating polychaete genus Osedax

Osedax, the deep-sea annelid found at sunken whalefalls, is known to host Oceanospirillales bacterial endosymbionts intracellularly in specialized roots, which help it feed exclusively on vertebrate bones. Past studies, however, have also made mention of external bacteria on their trunks. During a 14-yr study, we reveal a dynamic, yet persistent, shift of Campylobacterales integrated into the epidermis of Osedax, which change over time as the whale carcass degrades on the sea floor. The Campylobacterales associated with seven species of Osedax, which comprise 67% of the bacterial community on the trunk, appear initially dominated by the genus Arcobacter (at early time points <24 mo), the Sulfurospirillum at intermediate stages (~50 mo), and the Sulfurimonas at later stages (>140 mo) of whale carcass decomposition. Metagenome analysis of the epibiont metabolic capabilities suggests potential for a transition from heterotrophy to autotrophy and differences in their capacity to metabolize oxygen, carbon, nitrogen, and sulfur. Compared to free-living relatives, the Osedax epibiont genomes were enriched in transposable elements, implicating genetic exchange on the host surface, and contained numerous secretions systems with eukaryotic-like protein (ELP) domains, suggesting a long evolutionary history with these enigmatic, yet widely distributed deep-sea worms. IMPORTANCE Symbiotic associations are widespread in nature and we can expect to find them in every type of ecological niche. In the last twenty years, the myriad of functions, interactions and species comprising microbe-host associations has fueled a surge of interest and appreciation for symbiosis. During this 14-year study, we reveal a dynamic population of bacterial epibionts, integrated into the epidermis of 7 species of a deep-sea worm group that feeds exclusively on the remains of marine mammals. The bacterial genomes provide clues of a long evolutionary history with these enigmatic worms. On the host surface, they exchange genes and appear to undergo ecological succession, as the whale carcass habitat degrades over time, similar to what is observed for some free-living communities. These, and other annelid worms are important keystone species for diverse deep-sea environments, yet the role of attached external bacteria in supporting host health has received relatively little attention.

Long-read assembled metagenomic approaches improve our understanding on metabolic potentials of microbial community in mangrove sediments

BACKGROUND

Mangrove wetlands are coastal ecosystems with important ecological features and provide habitats for diverse microorganisms with key roles in nutrient and biogeochemical cycling. However, the overall metabolic potentials and ecological roles of microbial community in mangrove sediment are remained unanswered. In current study, the microbial and metabolic profiles of prokaryotic and fungal communities in mangrove sediments were investigated using metagenomic analysis based on PacBio single-molecule real time (SMRT) and Illumina sequencing techniques.

RESULTS

Comparing to Illumina short reads, the incorporation of PacBio long reads significantly contributed to more contiguous assemblies, yielded more than doubled high-quality metagenome-assembled genomes (MAGs), and improved the novelty of the MAGs. Further metabolic reconstruction for recovered MAGs showed that prokaryotes potentially played an essential role in carbon cycling in mangrove sediment, displaying versatile metabolic potential for degrading organic carbons, fermentation, autotrophy, and carbon fixation. Mangrove fungi also functioned as a player in carbon cycling, potentially involved in the degradation of various carbohydrate and peptide substrates. Notably, a new candidate bacterial phylum named as Candidatus Cosmopoliota with a ubiquitous distribution is proposed. Genomic analysis revealed that this new phylum is capable of utilizing various types of organic substrates, anaerobic fermentation, and carbon fixation with the Wood-Ljungdahl (WL) pathway and the reverse tricarboxylic acid (rTCA) cycle.

CONCLUSIONS

The study not only highlights the advantages of HiSeq-PacBio Hybrid assembly for a more complete profiling of environmental microbiomes but also expands our understanding of the microbial diversity and potential roles of distinct microbial groups in biogeochemical cycling in mangrove sediment. Video Abstract.

Comparative analysis of organellar genomes between diploid and tetraploid Chrysanthemum indicum with its relatives

Chrysanthemum indicum, a species native to Eastern Asia is well known as one of the progenitor species of the cultivated Chrysanthemum which is grown for its ornamental and medicinal value. Previous genomic studies on Chrysanthemum have largely ignored the dynamics of plastid genome (plastome) and mitochondria genome (mitogenome) evolution when analyzing this plant lineage. In this study, we sequenced and assembled the plastomes and mitogenomes of diploid and tetraploid C. indicum as well as the morphologically divergent variety C. indicum var. aromaticum. We used published data from 27 species with both plastome and mitogenome complete sequences to explore differences in sequence evolution between the organellar genomes. The size and structure of organellar genome between diploid and tetraploid C. indicum were generally similar but the tetraploid C. indicum and C. indicum var. aromaticum were found to contain unique sequences in the mitogenomes which also contained previously undescribed open reading frames (ORFs). Across Chrysanthemum mitogenome structure varied greatly but sequences transferred from plastomes in to the mitogenomes were conserved. Finally, differences observed between mitogenome and plastome gene trees may be the result of the difference in the rate of sequence evolution between genes in these two genomes. In total the findings presented here greatly expand the resources for studying Chrysanthemum organellar genome evolution with possible applications to conservation, breeding, and gene banking in the future.

Long-Read Metagenomics of Marine Microbes Reveals Diversely Expressed Secondary Metabolites

Microbial secondary metabolites play crucial roles in microbial competition, communication, resource acquisition, antibiotic production, and a variety of other biotechnological processes. The retrieval of full-length BGC (biosynthetic gene cluster) sequences from uncultivated bacteria is difficult due to the technical constraints of short-read sequencing, making it impossible to determine BGC diversity. Using long-read sequencing and genome mining, 339 mainly full-length BGCs were recovered in this study, illuminating the wide range of BGCs from uncultivated lineages discovered in seawater from Aoshan Bay, Yellow Sea, China. Many extremely diverse BGCs were discovered in bacterial phyla such as Proteobacteria, Bacteroidota, Acidobacteriota, and Verrucomicrobiota as well as the previously uncultured archaeal phylum "Candidatus Thermoplasmatota." The data from metatranscriptomics showed that 30.1% of secondary metabolic genes were being expressed, and they also revealed the expression pattern of BGC core biosynthetic genes and tailoring enzymes. Taken together, our results demonstrate that long-read metagenomic sequencing combined with metatranscriptomic analysis provides a direct view into the functional expression of BGCs in environmental processes. IMPORTANCE Genome mining of metagenomic data has become the preferred method for the bioprospecting of novel compounds by cataloguing secondary metabolite potential. However, the accurate detection of BGCs requires unfragmented genomic assemblies, which have been technically difficult to obtain from metagenomes until recently with new long-read technologies. We used high-quality metagenome-assembled genomes generated from long-read data to determine the biosynthetic potential of microbes found in the surface water of the Yellow Sea. We recovered 339 highly diverse and mostly full-length BGCs from largely uncultured and underexplored bacterial and archaeal phyla. Additionally, we present long-read metagenomic sequencing combined with metatranscriptomic analysis as a potential method for gaining access to the largely underutilized genetic reservoir of specialized metabolite gene clusters in the majority of microbes that are not cultured. The combination of long-read metagenomic and metatranscriptomic analyses is significant because it can more accurately assess the mechanisms of microbial adaptation to the environment through BGC expression based on metatranscriptomic data.

No Isolate, No Problem: Using a Novel Insertion Sequence PCR to Link Rats to Human Shigellosis Cases in an Underserved Urban Community

During an investigation into a cluster of Shigella flexneri serotype 2a cases in an underserved community, we assessed the relatedness of human and rat S. flexneri isolates utilizing a novel PCR targeting insertion sites (IS-PCR) of mobile elements in the Shigella genome characteristic of the cluster strain. Whole-genome sequences of S. flexneri (n = 50) associated with the cluster were analyzed. De novo genome assemblies were analyzed by a Geneious V10.2.6 motif search, and two unique IS were identified in all human Shigella sequences of the local cluster. Hydrolysis probe PCR assays were designed to detect these sequences consisting of forward and reverse primers to amplify across each insertion site and a hydrolysis probe spanning the insertion site. IS-PCR was performed for three Shigella PCR-positive culture-negative rat intestine specimens from this community. Both insertion sites were detected in the de novo genome assemblies of all clinical S. flexneri isolates (n = 50). Two of the three PCR-positive culture-negative rat samples were positive for both unique ISs identified in the human S. flexneri isolates, suggesting that the rat Shigella species strains were closely related to the human strains in the cluster. The cycle threshold (Ct) values were >35, indicating that the bacterial load was very low in the rat samples. Two unique IS were identified in clinical isolates from a community S. flexneri cluster. Both IS targets were identified in PCR-positive (Shigella spp.), culture-negative rat tissue and clinical isolates from humans, indicating relatedness. IMPORTANCE This article describes a novel molecular method to show relatedness between bacterial infections, which may not be able to grow in the laboratory due to treatment with antibiotics or for bacteria requiring unique conditions to grow well. Uniquely, we applied this technique to Shigella isolates from human cases associated with a local cluster in an underserved community, as well as rat samples from the same community. We believe that this novel approach can serve as a complementary method to support outbreak/cluster investigation for Shigella spp.

Chromosome-level genome assemblies of two parasitoid biocontrol wasps reveal the parthenogenesis mechanism and an associated novel virus

BACKGROUND

Biocontrol is a key technology for the control of pest species. Microctonus parasitoid wasps (Hymenoptera: Braconidae) have been released in Aotearoa New Zealand as biocontrol agents, targeting three different pest weevil species. Despite their value as biocontrol agents, no genome assemblies are currently available for these Microctonus wasps, limiting investigations into key biological differences between the different species and strains.

METHODS AND FINDINGS

Here we present high-quality genomes for Microctonus hyperodae and Microctonus aethiopoides, assembled with short read sequencing and Hi-C scaffolding. These assemblies have total lengths of 106.7 Mb for M. hyperodae and 129.2 Mb for M. aethiopoides, with scaffold N50 values of 9 Mb and 23 Mb respectively. With these assemblies we investigated differences in reproductive mechanisms, and association with viruses between Microctonus wasps. Meiosis-specific genes are conserved in asexual Microctonus, with in-situ hybridisation validating expression of one of these genes in the ovaries of asexual Microctonus aethiopoides. This implies asexual reproduction in these Microctonus wasps involves meiosis, with the potential for sexual reproduction maintained. Investigation of viral gene content revealed candidate genes that may be involved in virus-like particle production in M. aethiopoides, as well as a novel virus infecting M. hyperodae, for which a complete genome was assembled.

CONCLUSION AND SIGNIFICANCE

These are the first published genomes for Microctonus wasps which have been deployed as biocontrol agents, in Aotearoa New Zealand. These assemblies will be valuable resources for continued investigation and monitoring of these biocontrol systems. Understanding the biology underpinning Microctonus biocontrol is crucial if we are to maintain its efficacy, or in the case of M. hyperodae to understand what may have influenced the significant decline of biocontrol efficacy. The potential for sexual reproduction in asexual Microctonus is significant given that empirical modelling suggests this asexual reproduction is likely to have contributed to biocontrol decline. Furthermore the identification of a novel virus in M. hyperodae highlights a previously unknown aspect of this biocontrol system, which may contribute to premature mortality of the host pest. These findings have potential to be exploited in future in attempt to increase the effectiveness of M. hyperodae biocontrol.

Comparative Genomics Identifies Conserved and Variable TAL Effectors in African Strains of the Cotton Pathogen Xanthomonas citri pv. malvacearum

Strains of Xanthomonas citri pv. malvacearum cause bacterial blight of cotton, a potentially serious threat to cotton production worldwide, including in sub-Saharan countries. Development of disease symptoms, such as water soaking, has been linked to the activity of a class of type 3 effectors, called transcription activator-like (TAL) effectors, which induce susceptibility genes in the host's cells. To gain further insight into the global diversity of the pathogen, to elucidate their repertoires of TAL effector genes, and to better understand the evolution of these genes in the cotton-pathogenic xanthomonads, we sequenced the genomes of three African strains of X. citri pv. malvacearum using nanopore technology. We show that the cotton-pathogenic pathovar of X. citri is a monophyletic lineage containing at least three distinct genetic subclades, which appear to be mirrored by their repertoires of TAL effectors. We observed an atypical level of TAL effector gene pseudogenization, which might be related to resistance genes that are deployed to control the disease. Our work thus contributes to a better understanding of the conservation and importance of TAL effectors in the interaction with the host plant, which can inform strategies for improving resistance against bacterial blight in cotton.

Enrichment and characterization of a nitric oxide-reducing microbial community in a continuous bioreactor

Nitric oxide (NO) is a highly reactive and climate-active molecule and a key intermediate in the microbial nitrogen cycle. Despite its role in the evolution of denitrification and aerobic respiration, high redox potential and capacity to sustain microbial growth, our understanding of NO-reducing microorganisms remains limited due to the absence of NO-reducing microbial cultures obtained directly from the environment using NO as a substrate. Here, using a continuous bioreactor and a constant supply of NO as the sole electron acceptor, we enriched and characterized a microbial community dominated by two previously unknown microorganisms that grow at nanomolar NO concentrations and survive high amounts (>6 µM) of this toxic gas, reducing it to N2 with little to non-detectable production of the greenhouse gas nitrous oxide. These results provide insight into the physiology of NO-reducing microorganisms, which have pivotal roles in the control of climate-active gases, waste removal, and evolution of nitrate and oxygen respiration.

Complete genome sequence of Staphylococcus casei strain DSM 15096

We present the first complete genome sequence of the species Staphylococcus casei . Strain DSM 15096 was sequenced with a hybrid approach using Oxford Nanopore Technologies long-read sequencing and Illumina short-read sequencing. The assembled sequences produced a 2 808 898 bp chromosomal molecule containing 2705 predicted genes, plus eight plasmids.

Characteristics and Whole-Genome Analysis of Limosilactobacillus fermentum Phage LFP02

Limosilactobacillus fermentum is a bacterium widely used in food production, medicine, and industrial fermentation. However, fermentation could fail due to phage contamination. L. fermentum bacteriophage LFP02 can be induced from L. fermentum IMAU 32579 using mitomycin C. To better understand the characteristics of this phage, its physiological and genomic characteristics were evaluated. The results showed that its optimal multiplicity of infection was 0.01, and the burst size was 148.03 ± 2.65 pfu/infective center. Compared to temperature, pH had a more obvious influence on phage viability, although its adsorption capacity was not affected by the divalent cations (Ca²⁺ and Mg²⁺) or chloramphenicol. Its genome size was 43,789 bp and the GC content was 46.06%, including 53 functional proteins. Compared to other L. fermentum phages, phage LFP02 had chromosome deletion, insertion, and inversion, which demonstrated that it was a novel phage. This study could expand the knowledge of the biological characteristics of L. fermentum bacteriophages and provide some theoretical basis for bacteriophage prevention during fermentation.

Distinct horizontal transfer mechanisms for type I and type V CRISPR-associated transposons

CRISPR-associated transposons (CASTs) co-opt CRISPR-Cas proteins and Tn7-family transposons for RNA-guided vertical and horizontal transmission. CASTs encode minimal CRISPR arrays but can't acquire new spacers. Here, we show that CASTs instead co-opt defense-associated CRISPR arrays for horizontal transmission. A bioinformatic analysis shows that all CAST sub-types co-occur with defense-associated CRISPR-Cas systems. Using an E. coli quantitative transposition assay, we show that CASTs use CRISPR RNAs (crRNAs) from these defense systems for horizontal gene transfer. A high-resolution structure of the type I-F CAST-Cascade in complex with a type III-B crRNA reveals that Cas6 recognizes direct repeats via sequence-independent π - π interactions. In addition to using heterologous CRISPR arrays, type V CASTs can also transpose via a crRNA-independent unguided mechanism, even when the S15 co-factor is over-expressed. Over-expressing S15 and the trans-activating CRISPR RNA (tracrRNA) or a single guide RNA (sgRNA) reduces, but does not abrogate, off-target integration for type V CASTs. Exploiting new spacers in defense-associated CRISPR arrays explains how CASTs horizontally transfer to new hosts. More broadly, this work will guide further efforts to engineer the activity and specificity of CASTs for gene editing applications.

Efficient High-Quality Metagenome Assembly from Long Accurate Reads using Minimizer-space de Bruijn Graphs

We introduce a novel metagenomics assembler for high-accuracy long reads. Our approach, implemented as metaMDBG, combines highly efficient de Bruijn graph assembly in minimizer space, with both a multi-k' approach for dealing with variations in genome coverage depth and an abundance-based filtering strategy for simplifying strain complexity. The resulting algorithm is more efficient than the state-of-the-art but with better assembly results. metaMDBG was 1.5 to 12 times faster than competing assemblers and requires between one-tenth and one-thirtieth of the memory across a range of data sets. We obtained up to twice as many high-quality circularised prokaryotic metagenome assembled genomes (MAGs) on the most complex communities, and a better recovery of viruses and plasmids. metaMDBG performs particularly well for abundant organisms whilst being robust to the presence of strain diversity. The result is that for the first time it is possible to efficiently reconstruct the majority of complex communities by abundance as near-complete MAGs.

Assessing labelled carbon assimilation from poly butylene adipate-co-terephthalate (PBAT) monomers during thermophilic anaerobic digestion

PBAT (poly butylene adipate-co-terephthalate) is a widely used biodegradable plastic, but the knowledge about its metabolization in anaerobic environments is very limited. In this study, the anaerobic digester sludge from a municipal wastewater treatment plant was used as inoculum to investigate the biodegradability of PBAT monomers in thermophilic conditions. The research employs a combination of ¹³C-labelled monomers and proteogenomics to track the labelled carbon and identify the microorganisms involved. A total of 122 labelled peptides of interest were identified for adipic acid (AA) and 1,4-butanedio (BD). Through the time-dependent isotopic enrichment and isotopic profile distributions, Bacteroides, Ichthyobacterium, and Methanosarcina were proven to be directly involved in the metabolization of at least one monomer. This study provides a first insight into the identity and genomic potential of microorganisms responsible for biodegradability of PBAT monomers during anaerobic digestion under thermophilic conditions.

SemiBin2: self-supervised contrastive learning leads to better MAGs for short- and long-read sequencing

MOTIVATION

Metagenomic binning methods to reconstruct metagenome-assembled genomes (MAGs) from environmental samples have been widely used in large-scale metagenomic studies. The recently proposed semi-supervised binning method, SemiBin, achieved state-of-the-art binning results in several environments. However, this required annotating contigs, a computationally costly and potentially biased process.

RESULTS

We propose SemiBin2, which uses self-supervised learning to learn feature embeddings from the contigs. In simulated and real datasets, we show that self-supervised learning achieves better results than the semi-supervised learning used in SemiBin1 and that SemiBin2 outperforms other state-of-the-art binners. Compared to SemiBin1, SemiBin2 can reconstruct 8.3-21.5% more high-quality bins and requires only 25% of the running time and 11% of peak memory usage in real short-read sequencing samples. To extend SemiBin2 to long-read data, we also propose ensemble-based DBSCAN clustering algorithm, resulting in 13.1-26.3% more high-quality genomes than the second best binner for long-read data.

AVAILABILITY AND IMPLEMENTATION

SemiBin2 is available as open source software at https://github.com/BigDataBiology/SemiBin/ and the analysis scripts used in the study can be found at https://github.com/BigDataBiology/SemiBin2_benchmark.

Best Practices in Microbial Experimental Evolution: Using Reporters and Long-Read Sequencing to Identify Copy Number Variation in Experimental Evolution

Copy number variants (CNVs), comprising gene amplifications and deletions, are a pervasive class of heritable variation. CNVs play a key role in rapid adaptation in both natural, and experimental, evolution. However, despite the advent of new DNA sequencing technologies, detection and quantification of CNVs in heterogeneous populations has remained challenging. Here, we summarize recent advances in the use of CNV reporters that provide a facile means of quantifying de novo CNVs at a specific locus in the genome, and nanopore sequencing, for resolving the often complex structures of CNVs. We provide guidance for the engineering and analysis of CNV reporters and practical guidelines for single-cell analysis of CNVs using flow cytometry. We summarize recent advances in nanopore sequencing, discuss the utility of this technology, and provide guidance for the bioinformatic analysis of these data to define the molecular structure of CNVs. The combination of reporter systems for tracking and isolating CNV lineages and long-read DNA sequencing for characterizing CNV structures enables unprecedented resolution of the mechanisms by which CNVs are generated and their evolutionary dynamics.

Translocation of Oral Microbiota into the Pancreatic Ductal Adenocarcinoma Tumor Microenvironment

Oral dysbiosis has long been associated with pancreatic ductal adenocarcinoma (PDAC). In this work, we explore the relationship between the oral and tumor microbiomes of patients diagnosed with PDAC. Salivary and tumor microbiomes were analyzed using a variety of sequencing methods, resulting in a high prevalence and relative abundance of oral bacteria, particularly Veillonella and Streptococcus, within tumor tissue. The most prevalent and abundant taxon found within both saliva and tumor tissue samples, Veillonella atypica, was cultured from patient saliva, sequenced and annotated, identifying genes that potentially contribute to tumorigenesis. High sequence similarity was observed between sequences recovered from patient matched saliva and tumor tissue, indicating that the taxa found in PDAC tumors may derive from the mouth. These findings may have clinical implications in the care and treatment of patients diagnosed with PDAC.

Distinct genomic routes underlie transitions to specialised symbiotic lifestyles in deep-sea annelid worms

Bacterial symbioses allow annelids to colonise extreme ecological niches, such as hydrothermal vents and whale falls. Yet, the genetic principles sustaining these symbioses remain unclear. Here, we show that different genomic adaptations underpin the symbioses of phylogenetically related annelids with distinct nutritional strategies. Genome compaction and extensive gene losses distinguish the heterotrophic symbiosis of the bone-eating worm Osedax frankpressi from the chemoautotrophic symbiosis of deep-sea Vestimentifera. Osedax's endosymbionts complement many of the host's metabolic deficiencies, including the loss of pathways to recycle nitrogen and synthesise some amino acids. Osedax's endosymbionts possess the glyoxylate cycle, which could allow more efficient catabolism of bone-derived nutrients and the production of carbohydrates from fatty acids. Unlike in most Vestimentifera, innate immunity genes are reduced in O. frankpressi, which, however, has an expansion of matrix metalloproteases to digest collagen. Our study supports that distinct nutritional interactions influence host genome evolution differently in highly specialised symbioses.