IDBA-UD: a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Diverse microbiome functions, limited temporal variation and substantial genomic conservation within sedimentary and granite rock deep underground research laboratories

BACKGROUND

Underground research laboratories (URLs) provide a window on the deep biosphere and enable investigation of potential microbial impacts on nuclear waste, CO2 and H2 stored in the subsurface. We carried out the first multi-year study of groundwater microbiomes sampled from defined intervals between 140 and 400 m below the surface of the Horonobe and Mizunami URLs, Japan.

RESULTS

We reconstructed draft genomes for > 90% of all organisms detected over a four year period. The Horonobe and Mizunami microbiomes are dissimilar, likely because the Mizunami URL is hosted in granitic rock and the Horonobe URL in sedimentary rock. Despite this, hydrogen metabolism, rubisco-based CO2 fixation, reduction of nitrogen compounds and sulfate reduction are well represented functions in microbiomes from both URLs, although methane metabolism is more prevalent at the organic- and CO2-rich Horonobe URL. High fluid flow zones and proximity to subsurface tunnels select for candidate phyla radiation bacteria in the Mizunami URL. We detected near-identical genotypes for approximately one third of all genomically defined organisms at multiple depths within the Horonobe URL. This cannot be explained by inactivity, as in situ growth was detected for some bacteria, albeit at slow rates. Given the current low hydraulic conductivity and groundwater compositional heterogeneity, ongoing inter-site strain dispersal seems unlikely. Alternatively, the Horonobe URL microbiome homogeneity may be explained by higher groundwater mobility during the last glacial period. Genotypically-defined species closely related to those detected in the URLs were identified in three other subsurface environments in the USA. Thus, dispersal rates between widely separated underground sites may be fast enough relative to mutation rates to have precluded substantial divergence in species composition. Species overlaps between subsurface locations on different continents constrain expectations regarding the scale of global subsurface biodiversity.

CONCLUSIONS

Our analyses reveal microbiome stability in the sedimentary rocks and surprising microbial community compositional and genotypic overlap over sites separated by hundreds of meters of rock, potentially explained by dispersal via slow groundwater flow or during a prior hydrological regime. Overall, microbiome and geochemical stability over the study period has important implications for underground storage applications.

Alterations in Gut Microbiota Correlate With Hematological Injuries Induced by Radiation in Beagles

Dynamics of gut microbiota and their associations with the corresponding hematological injuries postradiation remain to be elucidated. Using single whole-body exposure to 60Co-γ ray radiation at the sublethal dose of 2.5 Gy, we developed a beagle model of acute radiation syndrome (ARS) and then monitored the longitudinal changes of gut microbiome and hematology for 45 days. We found that the absolute counts of circulating lymphocytes, neutrophils, and platelets were sharply declined postradiation, accompanied by a largely shifted composition of gut microbiome that manifested as a significantly increased ratio of Firmicutes to Bacteroidetes. In irradiated beagles, alterations in hematological parameters reached a nadir on day 14, sustaining for 1 week, which were gradually returned to the normal levels thereafter. However, no structural recovery of gut microbiota was observed throughout the study. Fecal metagenomics revealed that irradiation increased the relative abundances of genus Streptococcus, species Lactobacillus animalis and Lactobacillus murinus, but decreased those of genera Prevotella and Bacteroides. Metagenomic functions prediction demonstrated that 26 altered KEGG pathways were significantly enriched on Day 14 and 35 postradiation. Furthermore, a total of 43 bacterial species were found to correlate well with hematological parameters by Spearman's analysis. Our results provide an insight into the longitudinal changes in intestinal microbiota at different clinical stages during ARS in canine. Several key microbes those tightly associated with the hematological alterations may serve as biomarkers to discriminate the different phases of host with ARS.

Detection and Molecular Diversity of Brucella melitensis in Pastoral Livestock in North-Eastern Ethiopia

Brucellosis is a neglected zoonotic disease affecting livestock and humans that remains endemic in Ethiopia. Despite its prevalence, only a few studies have identified Brucella species circulating in livestock in the country. This study aimed to determine the Brucella species responsible for infections in livestock in the Afar region of Ethiopia and characterize the isolates using whole-genome single nucleotide polymorphism (wgSNP) analysis and in silico multi-locus sequence typing (MLST). Comparisons were made between Ethiopian Brucella and regional and global isolates to determine their phylogenetic relationships. Surveys conducted in May and October-November 2022 in six villages of the Amibara district involved the collection of vaginal swabs (n = 231) and milk samples (n = 17) from 32 sheep and 199 goats kept by 143 pastoral households reporting recent abortions in the animals. Brucella melitensis was detected in three sheep and 32 goats, i.e., 15% (35/231) of animals across 20% (29/143) of households using bacterial culture and PCR-based methods (bcsp31, AMOS, and Bruce-ladder multiplex PCR). Of the 35 positive animals, B. melitensis was isolated from 24 swabs, while the remaining 11 were culture-negative and detected only by PCR. The genomic DNA of the 24 isolates was sequenced using Illumina Novaseq 6000 and assembled using the SPAdes pipeline. Nine- and 21-locus MLST identified 23 isolates as genotype ST12, while one isolate could not be typed. The wgSNP-based phylogenetic analysis revealed that the Ethiopian isolates clustered within the African clade and were closely related to isolates from Somalia. Several virulence factors responsible for adhesion, intracellular survival, and regulatory functions were detected in all isolates. No antimicrobial resistance genes associated with resistance to drugs commonly used for treating brucellosis were detected. Since B. melitensis is prevalent in sheep and goats, vaccination with the B. melitensis Rev-1 vaccine is the recommended strategy in these pastoral systems to protect animal and human health.

Temperature differentially regulates estuarine microbial N2O production along a salinity gradient

Nitrous oxide (N2O) is atmospheric trace gas that contributes to climate change and affects stratospheric and ground-level ozone concentrations. Ammonia oxidizers and denitrifiers contribute to N2O emissions in estuarine waters. However, as an important climate factor, how temperature regulates microbial N2O production in estuarine water remains unclear. Here, we have employed stable isotope labeling techniques to demonstrate that the N2O production in estuarine waters exhibited differential thermal response patterns between nearshore and offshore regions. The optimal temperatures (Topt) for N2O production rates (N2OR) were higher at nearshore than offshore sites. 15N-labeled nitrite (15NO2-) experiments revealed that at the nearshore sites dominated by ammonia-oxidizing bacteria (AOB), the thermal tolerance of 15N-N2OR increases with increasing salinity, suggesting that N2O production by AOB-driven nitrifier denitrification may be co-regulated by temperature and salinity. Metatranscriptomic and metagenomic analyses of enriched water samples revealed that the denitrification pathway of AOB is the primary source of N2O, while clade II N2O-reducers dominated N2O consumption. Temperature regulated the expression patterns of nitrite reductase (nirK) and nitrous oxide reductase (nosZ) genes from different sources, thereby influencing N2O emissions in the system. Our findings contribute to understanding the sources of N2O in estuarine waters and their response to global warming.

The last of their kind: Is the genus Scutiger (Anura: Megophryidae) a relict element of the paleo-Transhimalaya biota?

The orographic evolution of the Himalaya-Tibet Mountain system continues to be a subject of controversy, leading to considerable uncertainty regarding the environment and surface elevation of the Tibetan Plateau during the Cenozoic era. As many geoscientific (but not paleontological) studies suggest, elevations close to modern heights exist in vast areas of Tibet since at least the late Paleogene, implicating the presence of large-scale alpine environments for more than 30 million years. To explore a recently proposed alternative model that assumes a warm temperate environment across paleo-Tibet, we carried out a phylogeographic survey using genomic analyses of samples covering the range of endemic lazy toads (Scutiger) across the Himalaya-Tibet orogen. We identified two main clades, with several, geographically distinct subclades. The long temporal gap between the stem and crown age of Scutiger may suggest high extinction rates. Diversification within the crown group, depending on the calibration, occurred either from the Mid-Miocene or Late-Miocene and continued until the Holocene. The present-day Himalayan Scutiger fauna could have evolved from lineages that existed on the southern edges of the paleo-Tibetan area (the Transhimalaya = Gangdese Shan), while extant species living on the eastern edge of the Plateau originated probably from the eastern edges of northern parts of the ancestral Tibetan area (Hoh Xil, Tanggula Shan). Based on the Mid-Miocene divergence time estimation and ancestral area reconstruction, we propose that uplift-associated aridification of a warm temperate Miocene-Tibet, coupled with high extirpation rates of ancestral populations, and species range shifts along drainage systems and epigenetic transverse valleys of the rising mountains, is a plausible scenario explaining the phylogenetic structure of Scutiger. This hypothesis aligns with the fossil record but conflicts with geoscientific concepts of high elevated Tibetan Plateau since the late Paleogene. Considering a Late-Miocene/Pliocene divergence time, an alternative scenario of dispersal from SE Asia into the East, Central, and West Himalaya cannot be excluded, although essential evolutionary and biogeographic aspects remain unresolved within this model.

Adjusted bacterial cooperation in anammox community to adapt to high ammonium in wastewater treatment plant

Bacterial cooperation is very important for anammox bacteria which perform low-carbon and energy-efficient nitrogen removal, yet its variation to adapt to high NH4 +-N concentration in actual wastewater treatment plants (WWTPs) remains unclear. Here, we found wide and varied cross-feedings of anammox bacteria and symbiotic bacteria in the two series connected full-scale reactors with different NH4 +-N concentrations (297.95 ± 54.84 and 76.03 ± 34.01 mg/L) treating sludge digester liquor. The uptake of vitamin B6 as highly effective antioxidants secreted by the symbiotic bacteria was beneficial for anammox bacteria to resist the high NH4 +-N concentration and varied dissolved oxygen (DO). When NH4 +-N concentration in influent (1785.46 ± 228.5 mg/L) increased, anammox bacteria tended to reduce the amino acids supply to symbiotic bacteria to save metabolic costs. A total of 26.1% bacterial generalists switched to specialists to increase the stability and functional heterogeneity of the microbial community at high NH4 +-N conditions. V/A-type ATPase for anammox bacteria to adapt to the change of NH4 +-N was highly important to strive against cellular alkalization caused by free ammonia. This study expands the understanding of the adjusted bacterial cooperation within anammox consortia at high NH4 +-N conditions, providing new insights into bacterial adaptation to adverse environments from a sociomicrobiology perspective.

Revealing the formation mechanisms of key flavor components during the fermentation of bamboo shoots by combining flavoromics and metagenomics

Microbial metabolism plays a critical role in the flavor development of Guangxi fermented bamboo shoots (GFBS). To clarify the role of microorganisms in flavor formation and predict the metabolic pathways of key characteristic flavor compounds, this study employed metabolomics, Odor Activity Value (OAV), and Taste Activity Value (TAV) calculations, integrated with Partial Least Squares Discriminant Analysis (PLS-DA), to investigate changes in GFBS flavors-represented by volatile flavor compounds, organic acids, and free amino acids-across a 30-day fermentation period. Metagenomic datasets were used to identify taxonomic and functional changes in the microbial community. As a result, 26 characteristic flavor compounds (OAV or TAV > 1) were identified in mature GFBS, and 23 differential flavor compounds were identified at different fermentation stages using PLS-DA (VIP > 1.2). The top 10 microbial genera associated with these characteristic flavor compounds were identified, including Acinetobacter, Enterobacter, Raoultella, Enterococcus, Klebsiella, Lactococcus, Leuconostoc, Weissella, Lactiplantibacillus and Limosilactobacillus. Based on these findings, a predictive metabolic network of key flavor compounds in GFBS was constructed, providing a comprehensive understanding of the diverse metabolic roles of microorganisms during fermentation. This work lays a theoretical foundation for the standardized production and quality control of GFBS flavor.

Increased rumen Prevotella enhances BCAA synthesis, leading to synergistically increased skeletal muscle in myostatin-knockout cattle

Myostatin (MSTN) is a negative regulator of muscle growth, and its relationship with the gut microbiota is not well understood. In this study, we observed increase muscle area and branched-chain amino acids (BCAAs), an energy source of muscle, in myostatin knockout (MSTN-KO) cattle. To explore the link between increased BCAAs and rumen microbiota, we performed metagenomic sequencing, metabolome analysis of rumen fluid, and muscle transcriptomics. MSTN-KO cattle showed a significant increase in the phylum Bacteroidota (formerly Bacteroidetes), particularly the genus Prevotella (P = 3.12e-04). Within this genus, Prevotella_sp._CAG:732, Prevotella_sp._MSX73, and Prevotella_sp._MA2016 showed significant upregulation of genes related to BCAA synthesis. Functional enrichment analysis indicated enrichment of BCAA synthesis-related pathways in both rumen metagenomes and metabolomes. Additionally, muscle transcriptomics indicated enrichment in muscle fiber and amino acid metabolism, with upregulation of solute carrier family genes, enhancing BCAA transport. These findings suggest that elevated rumen Prevotella in MSTN-KO cattle, combined with MSTN deletion, synergistically improves muscle growth through enhanced BCAA synthesis and transport.

MetaCompare 2.0: differential ranking of ecological and human health resistome risks

While numerous environmental factors contribute to the spread of antibiotic resistance genes (ARGs), quantifying their relative contributions remains a fundamental challenge. Similarly, it is important to differentiate acute human health risks from environmental exposure, versus broader ecological risk of ARG evolution and spread across microbial taxa. Recent studies have proposed various methods for achieving such aims. Here, we introduce MetaCompare 2.0, which improves upon original MetaCompare pipeline by differentiating indicators of human health resistome risk (potential for human pathogens of acute resistance concern to acquire ARGs) from ecological resistome risk (overall mobility of ARGs and potential for pathogen acquisition). The updated pipeline's sensitivity was demonstrated by analyzing diverse publicly-available metagenomes from wastewater, surface water, soil, sediment, human gut, and synthetic microbial communities. MetaCompare 2.0 provided distinct rankings of the metagenomes according to both human health resistome risk and ecological resistome risk, with both scores trending higher when influenced by anthropogenic impact or other stress. We evaluated the robustness of the pipeline to sequence assembly methods, sequencing depth, contig count, and metagenomic library coverage bias. The risk scores were remarkably consistent despite variations in these technological aspects. We packaged the improved pipeline into a publicly-available web service (http://metacompare.cs.vt.edu/) that provides an easy-to-use interface for computing resistome risk scores and visualizing results.

Targeted protein evolution in the gut microbiome by diversity-generating retroelements

Diversity-generating retroelements (DGRs) accelerate evolution by rapidly diversifying variable proteins. The human gastrointestinal microbiota harbors the greatest density of DGRs known in nature, suggesting they play adaptive roles in this environment. We identified >1,100 unique DGRs among human-associated Bacteroides species and discovered a subset that diversify adhesive components of Type V pili and related proteins. We show that Bacteroides DGRs are horizontally transferred across species, that some are highly active while others are tightly controlled, and that they preferentially alter the functional characteristics of ligand-binding residues on adhesive organelles. Specific variable protein sequences are enriched when Bacteroides strains compete with other commensal bacteria in gnotobiotic mice. Analysis of >2,700 DGRs from diverse phyla in mother-infant pairs shows that Bacteroides DGRs are preferentially transferred to vaginally delivered infants where they actively diversify. Our observations provide a foundation for understanding the roles of stochastic, targeted genome plasticity in shaping host-associated microbial communities.

Microplastics exacerbate the ecological risk of antibiotic resistance genes in wetland ecosystem

Wetlands are vital components of the global ecosystem, significantly influencing the retention and dissemination of microplastics (MPs) and antibiotic resistance genes (ARGs). However, the effects of different types of MPs on the environmental dynamics of ARGs within these ecosystems remain poorly understood. This study focused on the distribution and composition of ARGs associated with two primary types of MPs-polyethylene and polypropylene-within the Poyang Lake wetland, the largest freshwater lake in China, utilizing metagenomic analysis. The findings demonstrated that the bacterial communities and ARG profiles in the plastisphere were markedly distinct from those in the surrounding water. Specifically, thirteen opportunistic pathogens and forty subtypes of ARGs, primarily related to multidrug, bacitracin, and β-lactam resistance, were identified in the plastisphere. Notably, polyethylene exhibited four times more specific ARG subtypes than polypropylene. Procrustes analysis combined with network analysis indicated a lack of strong correlation between ARG abundance and bacterial populations, suggesting potential horizontal transfer of ARGs within the microbiota of the plastisphere. Additionally, three novel and functional β-lactamase genes were identified within this environment. This investigation highlights the role of MPs as reservoirs for ARGs, facilitating their exchange and posing risks to both ecological integrity and human health, thereby underscoring the need for increased attention in future research efforts.

A core microbiome signature as an indicator of health

The gut microbiota is crucial for human health, functioning as a complex adaptive system akin to a vital organ. To identify core health-relevant gut microbes, we followed the systems biology tenet that stable relationships signify core components. By analyzing metagenomic datasets from a high-fiber dietary intervention in type 2 diabetes and 26 case-control studies across 15 diseases, we identified a set of stably correlated genome pairs within co-abundance networks perturbed by dietary interventions and diseases. These genomes formed a "two competing guilds" (TCGs) model, with one guild specialized in fiber fermentation and butyrate production and the other characterized by virulence and antibiotic resistance. Our random forest models successfully distinguished cases from controls across multiple diseases and predicted immunotherapy outcomes through the use of these genomes. Our guild-based approach, which is genome specific, database independent, and interaction focused, identifies a core microbiome signature that serves as a holistic health indicator and a potential common target for health enhancement.

Comparison of microbial diversity and carbohydrate-active enzymes in the hindgut of two wood-feeding termites, Globitermes sulphureus (Blattaria: Termitidae) and Coptotermes formosanus (Blattaria: Rhinotermitidae)

BACKGROUND

Wood-feeding termites have been employed as sources of novel and highly efficient lignocellulolytic enzymes due to their ability to degrade lignocellulose efficiently. As a higher wood-feeding termite, Globitermes sulphureus (Blattaria: Termitidae) plays a crucial role as a decomposer in regions such as Vietnam, Singapore, Myanmar, and Yunnan, China. However, the diversity of its gut microbiome and carbohydrate-active enzymes (CAZymes) remains unexplored. Here, we analyzed the diversity of hindgut microbial communities and CAZymes in a higher wood-feeding termite, G. sulphureus, and a lower wood-feeding termite, Coptotermes formosanus (Blattaria: Rhinotermitidae).

RESULTS

16S rRNA sequencing revealed that Spirochaetota, Firmicutes, and Fibrobacterota were the dominant microbiota in the hindgut of the two termite species. At the phylum level, the relative abundances of Proteobacteria and Bacteroidota were significantly greater in the hindgut of C. formosanus than in G. sulphureus. At the genus level, the relative abundances of Candidatus_Azobacteroides and Escherichia-Shigella were significantly lower in the hindgut of G. sulphureus than in C. formosanus. Metagenomic analysis revealed that glycoside hydrolases (GHs) with cellulases and hemicellulases functions were not significantly different between G. sulphureus and C. formosanus. Interestingly, the cellulases in G. sulphureus were mainly GH5_2, GH5_4, GH6, GH9, and GH45, while the hemicellulases were mainly GH11, GH8, GH10, GH11, GH26, and GH53. In C. formosanus, the cellulases were mainly GH6 and GH9, and the hemicellulases were mainly GH5_7, GH5_21, GH10, GH12, and GH53. In addition, β-glucosidase, exo-β-1,4-glucanase, and endo-β-1,4-glucanase activities did not differ significantly between the two termite species, while xylanase activity was higher in G. sulphureus than in C. formosanus. The bacteria encoding GHs in G. sulphureus were mainly Firmicutes, Fibrobacterota, and Proteobacteria, whereas Bacteroidota and Spirochaetota were the main bacteria encoding GHs in C. formosanus.

CONCLUSIONS

Our findings characterized the microbial composition and differences in the hindgut microbiota of G. sulphureus and C. formosanus. Compared to C. formosanus, G. sulphureus is enriched in genes encoding for hemicellulase and debranching enzymes. It also highlights the rich diversity of GHs in the hindgut microbiota of G. sulphureus, including the GH5 subfamily, GH6, and GH48, with the GH6 and GH48 not previously reported in other higher termites. These results strengthen the understanding of the diversity of termite gut microbiota and CAZymes.

Metagenomes and metagenome assembled genomes from anaerobic digesters at three Canadian pulp and paper mills

We present a dataset of six metagenomes and 323 metagenome assembled genomes (MAGs) describing the microbial community of anaerobic digesters at three Canadian pulp and paper mills. Our objective was to assess the coding potential of the microbial community and obtain draft genomes of key organisms in the digesters.

Comparative metagenomics of tropical reef fishes show conserved core gut functions across hosts and diets with diet-related functional gene enrichments

Fish gut microbial communities are important for the breakdown and energy harvesting of the host diet. Microbes within the fish gut are selected by environmental and evolutionary factors. To understand how fish gut microbial communities are shaped by diet, three tropical fish species (hawkfish, Paracirrhites arcatus; yellow tang, Zebrasoma flavescens; and triggerfish, Rhinecanthus aculeatus) were fed piscivorous (fish meal pellets), herbivorous (seaweed), and invertivorous (shrimp) diets, respectively. From fecal samples, a total of 43 metagenome assembled genomes (MAGs) were recovered from all fish diet treatments. Each host-diet treatment harbored distinct microbial communities based on taxonomy, with Proteobacteria, Bacteroidota, and Firmicutes being the most represented. Based on their metagenomes, MAGs from all three host-diet treatments demonstrated a baseline ability to degrade proteinaceous, fatty acid, and simple carbohydrate inputs and carry out central carbon metabolism, lactate and formate fermentation, acetogenesis, nitrate respiration, and B vitamin synthesis. The herbivorous yellow tang harbored more functionally diverse MAGs with some complex polysaccharide degradation specialists, while the piscivorous hawkfish's MAGs were more specialized for the degradation of proteins. The invertivorous triggerfish's gut MAGs lacked many carbohydrate degrading capabilities, resulting in them being more specialized and functionally uniform. Across all treatments, several MAGs were able to participate in only individual steps of the degradation of complex polysaccharides, suggestive of microbial community networks that degrade complex inputs.

Genetic potential for aerobic respiration and denitrification in globally distributed respiratory endosymbionts

The endosymbiont Candidatus Azoamicus ciliaticola was proposed to generate ATP for its eukaryotic host, an anaerobic ciliate of the Plagiopylea class, fulfilling a function analogous to mitochondria in other eukaryotic cells. The discovery of this respiratory endosymbiosis has major implications for both evolutionary history and ecology of microbial eukaryotes. However, with only a single species described, knowledge of its environmental distribution and diversity is limited. Here we report four complete, circular metagenome assembled genomes (cMAGs) representing respiratory endosymbionts inhabiting groundwater in California, Ohio, and Germany. These cMAGs form two lineages comprising a monophyletic clade within the uncharacterized gammaproteobacterial order UBA6186, enabling evolutionary analysis of their key protein complexes. Strikingly, all four cMAGs encode a cytochrome cbb3 oxidase, which indicates that these endosymbionts have the capacity for aerobic respiration. Accordingly, we detect these respiratory endosymbionts in diverse habitats worldwide, thus further expanding the ecological scope of this respiratory symbiosis.

MAGqual: a stand-alone pipeline to assess the quality of metagenome-assembled genomes

BACKGROUND

Metagenomics, the whole genome sequencing of microbial communities, has provided insight into complex ecosystems. It has facilitated the discovery of novel microorganisms, explained community interactions and found applications in various fields. Advances in high-throughput and third-generation sequencing technologies have further fuelled its popularity. Nevertheless, managing the vast data produced and addressing variable dataset quality remain ongoing challenges. Another challenge arises from the number of assembly and binning strategies used across studies. Comparing datasets and analysis tools is complex as it requires the quantitative assessment of metagenome quality. The inherent limitations of metagenomic sequencing, which often involves sequencing complex communities, mean community members are challenging to interrogate with traditional culturing methods leading to many lacking reference sequences. MIMAG standards aim to provide a method to assess metagenome quality for comparison but have not been widely adopted.

RESULTS

To address the need for simple and quick metagenome quality assignation, here we introduce the pipeline MAGqual (Metagenome-Assembled Genome qualifier) and demonstrate its effectiveness at determining metagenomic dataset quality in the context of the MIMAG standards.

CONCLUSIONS

The MAGqual pipeline offers an accessible way to evaluate metagenome quality and generate metadata on a large scale. MAGqual is built in Snakemake to ensure readability and scalability, and its open-source nature promotes accessibility, community development, and ease of updates. MAGqual is built in Snakemake, R, and Python and is available under the MIT license on GitHub at https://github.com/ac1513/MAGqual . Video Abstract.

Deciphering the functional assembly of microbial communities driven by heavy metals in the tidal soils of Hangzhou Bay

Understanding the interaction between heavy metals and soil microbiomes is essential for maintaining ecosystem health and functionality in the face of persistent human-induced challenges. This study investigated the complex relationships between heavy metal contamination and the functional characteristics of soil microbial communities in the tidal soils of Hangzhou Bay, a region experiencing substantial environmental pressure due to its proximity to densely populated and industrialized regions. The north-shore sampling site showed moderate contaminations (mg/kg) of total arsenic (16.61 ± 1.13), cadmium (0.3 ± 0.05), copper (31.28 ± 1.23), nickel (37.44 ± 2.74), lead (34.29 ± 5.99), and zinc (120.8 ± 5.96), which are 1.29-2.94 times higher than the geochemical background values in Hangzhou Bay and adjacent areas. In contrast, the south-shore sampling site showed slightly higher levels of total arsenic (13.76 ± 1.35) and cadmium (0.13 ± 0.02) than the background values. Utilizing metagenomic sequencing, we decoded microbial functional genes essential for nitrogen, phosphorus, sulfur, and methane biogeochemical cycles. Although soil available nickel content was relatively low at 1 mg/kg, it exhibited strong associations with diverse microbial genes and biogeochemical pathways. Four key genes-hxlB, glpX, opd, and phny-emerged as pivotal players in the interactions with available nickel, suggesting the adaptability of microbial metabolic responses to heavy metal. Additionally, microbial genera such as Gemmatimonas and Ilumatobacter, which harbored diverse functional genes, demonstrated potential interactions with soil nickel. These findings highlight the importance of understanding heavy metal-soil microbiome dynamics for effective environmental management strategies in the tidal soils of Hangzhou Bay, with the goal of preserving ecosystem health and functionality amidst ongoing anthropogenic challenges.

Unveiling the role of emerging metagenomics for the examination of hypersaline environments

Hypersaline ecosystems are distributed all over the globe. They are subjected to poly-extreme stresses and are inhabited by halophilic microorganisms possessing multiple adaptations. The halophiles have many biotechnological applications such as nutrient supplements, antioxidant synthesis, salt tolerant enzyme production, osmolyte synthesis, biofuel production, electricity generation etc. However, halophiles are still underexplored in terms of complex ecological interactions and functions as compared to other niches. The advent of metagenomics and the recent advancement of next-generation sequencing tools have made it feasible to investigate the microflora of an ecosystem, its interactions and functions. Both target gene and shotgun metagenomic approaches are commonly employed for the taxonomic, phylogenetic, and functional analyses of the hypersaline microbial communities. This review discusses different types of hypersaline niches, their residential microflora, and an overview of the metagenomic approaches used to investigate them. Various applications, hurdles and the recent advancements in metagenomic approaches have also been focused on here for their better understanding and utilization in the study of hypersaline microbiome.

Plant virus community structuring is shaped by habitat heterogeneity and traits for host plant resource utilisation

Host plants provide resources critical to viruses and the spatial structuring of plant communities affects the niches available for colonisation and disease emergence. However, large gaps remain in the understanding of mechanisms that govern plant-virus disease ecology across heterogeneous plant assemblages. We combine high-throughput sequencing, network, and metacommunity approaches to test whether habitat heterogeneity in plant community composition corresponded with virus resource utilisation traits of transmission mode and host range. A majority of viruses exhibited habitat specificity, with communities connected by key generalist viruses and potential host reservoirs. There was an association between habitat heterogeneity and virus community structuring, and between virus community structuring and resource utilisation traits of host range and transmission. The relationship between virus species distributions and virus trait responses to habitat heterogeneity was scale-dependent, being stronger at finer (site) than larger (habitat) spatial scales. Results indicate that habitat heterogeneity has a part in plant virus community assembly, and virus community structuring corresponds to virus trait responses that vary with the scale of observation. Distinctions in virus communities caused by plant resource compartmentalisation can be used to track trait responses of viruses to hosts important in forecasting disease emergence.

A stepwise guide for pangenome development in crop plants: an alfalfa (Medicago sativa) case study

BACKGROUND

The concept of pangenomics and the importance of structural variants is gaining recognition within the plant genomics community. Due to advancements in sequencing and computational technology, it has become feasible to sequence the entire genome of numerous individuals of a single species at a reasonable cost. Pangenomes have been constructed for many major diploid crops, including rice, maize, soybean, sorghum, pearl millet, peas, sunflower, grapes, and mustards. However, pangenomes for polyploid species are relatively scarce and are available in only few crops including wheat, cotton, rapeseed, and potatoes.

MAIN BODY

In this review, we explore the various methods used in crop pangenome development, discussing the challenges and implications of these techniques based on insights from published pangenome studies. We offer a systematic guide and discuss the tools available for constructing a pangenome and conducting downstream analyses. Alfalfa, a highly heterozygous, cross pollinated and autotetraploid forage crop species, is used as an example to discuss the concerns and challenges offered by polyploid crop species. We conducted a comparative analysis using linear and graph-based methods by constructing an alfalfa graph pangenome using three publicly available genome assemblies. To illustrate the intricacies captured by pangenome graphs for a complex crop genome, we used five different gene sequences and aligned them against the three graph-based pangenomes. The comparison of the three graph pangenome methods reveals notable variations in the genomic variation captured by each pipeline.

CONCLUSION

Pangenome resources are proving invaluable by offering insights into core and dispensable genes, novel gene discovery, and genome-wide patterns of variation. Developing user-friendly online portals for linear pangenome visualization has made these resources accessible to the broader scientific and breeding community. However, challenges remain with graph-based pangenomes including compatibility with other tools, extraction of sequence for regions of interest, and visualization of genetic variation captured in pangenome graphs. These issues necessitate further refinement of tools and pipelines to effectively address the complexities of polyploid, highly heterozygous, and cross-pollinated species.

Eubacterium siraeum suppresses fat deposition via decreasing the tyrosine-mediated PI3K/AKT signaling pathway in high-fat diet-induced obesity

BACKGROUND

Obesity in humans can lead to chronic diseases such as diabetes and cardiovascular disease. Similarly, subcutaneous fat (SCF) in pigs affects feed utilization, and excessive SCF can reduce the feed efficiency of pigs. Therefore, identifying factors that suppress fat deposition is particularly important. Numerous studies have implicated the gut microbiome in pigs' fat deposition, but research into its suppression remains scarce. The Lulai black pig (LL) is a hybrid breed derived from the Laiwu pig (LW) and the Yorkshire pig, with lower levels of SCF compared to the LW. In this study, we focused on these breeds to identify microbiota that regulate fat deposition. The key questions were: Which microbial populations reduce fat in LL pigs compared to LW pigs, and what is the underlying regulatory mechanism?

RESULTS

In this study, we identified four different microbial strains, Eubacterium siraeum, Treponema bryantii, Clostridium sp. CAG:413, and Jeotgalibaca dankookensis, prevalent in both LW and LL pigs. Blood metabolome analysis revealed 49 differential metabolites, including tanshinone IIA and royal jelly acid, known for their anti-adipogenic properties. E. siraeum was strongly correlated with these metabolites, and its genes and metabolites were enriched in pathways linked to fatty acid degradation, glycerophospholipid, and glycerolipid metabolism. In vivo mouse experiments confirmed that E. siraeum metabolites curb weight gain, reduce SCF adipocyte size, increase the number of brown adipocytes, and regulate leptin, IL-6, and insulin secretion. Finally, we found that one important pathway through which E. siraeum inhibits fat deposition is by suppressing the phosphorylation of key proteins in the PI3K/AKT signaling pathway through the reduction of tyrosine.

CONCLUSIONS

We compared LW and LL pigs using fecal metagenomics, metabolomics, and blood metabolomics, identifying E. siraeum as a strain linked to fat deposition. Oral administration experiments in mice demonstrated that E. siraeum effectively inhibits fat accumulation, primarily through the suppression of the PI3K/AKT signaling pathway, a critical regulator of lipid metabolism. These findings provide a valuable theoretical basis for improving pork quality and offer insights relevant to the study of human obesity and related chronic metabolic diseases. Video Abstract.

Revealing salt concentration for microbial balance and metabolite enrichment in secondary fortified fermented soy sauce: A multi-omics perspective

This study examined the impact of varying salt concentrations on microbiota, physicochemical properties, and metabolites in a secondary fortified fermentation process using multi-omics techniques. It aimed to determine the influence of salt stress on microbiota shifts and metabolic activities. The findings demonstrated that moderate salt reduction (MS) was found to enhance moromi's flavor and quality, while mitigating the negative effects of excessive low salt (LS). MS samples had 1.22, 1.13, and 2.92 times more amino acid nitrogen (AAN), non-volatiles, and volatiles, respectively, than high salt (HS) samples. In contrast, lactic acid and biogenic amines in LS samples were 1.56 g/100 g and 4115.11 mg/kg, respectively, decreasing to 0.15 g/100 g and 176.76 mg/kg in MS samples. Additionally, the contents of ethanol and small peptides increased in MS due to the growth of specific functional microorganisms such as Staphylococcus gallinarum, Weissella confusa, and Zygosaccharomyces rouxii, while food-borne pathogens were inhibited. Network analysis revealed that the core microbial interactions were enhanced in MS samples, promoting a balanced fermentation environment. Redundancy analysis (RDA) and correlation analyses underscored that the physicochemical properties significantly impacted bacterial community structure and the correlations between key microbes and flavor compounds. These findings provided a theoretical foundation for developing innovative reduced-salt fermentation techniques, contributing to the sustainable production of high-quality soy sauce.

Biosynthetic gene clusters from uncultivated soil bacteria of the Atacama Desert

Soil microorganisms mediate several biological processes through the secretion of natural products synthesized in specialized metabolic pathways, yet functional characterization in ecological contexts remains challenging. Using culture-independent metagenomic analyses of microbial DNA derived directly from soil samples, we examined the potential of biosynthetic gene clusters (BGCs) from six bacterial communities distributed along an altitudinal gradient of the Andes Mountains in the Atacama Desert. We mined 38 metagenome-assembled genomes (MAGs) and identified 168 BGCs. Results indicated that most predicted BGCs were classified as non-ribosomal-peptides (NRP), post-translational modified peptides (RiPP), and terpenes, which were mainly identified in genomes of species from Acidobacteriota and Proteobacteria phyla. Based on BGC composition according to types of core biosynthetic genes, six clusters of MAGs were observed, three of them with predominance for a single phylum, of which two also showed specificity to a single sampling site. Comparative analyses of accessory genes in BGCs showed associations between membrane transporters and other protein domains involved in specialized metabolism with classes of biosynthetic cores, such as resistance-nodulation-cell division (RND) multidrug efflux pumps with RiPPs and the iron-dependent transporter TonB with terpenes. Our findings increase knowledge regarding the biosynthetic potential of uncultured bacteria inhabiting pristine locations from one of the oldest and driest nonpolar deserts on Earth.IMPORTANCEMuch of what we know about specialized metabolites in the Atacama Desert, including Andean ecosystems, comes from isolated microorganisms intended for drug development and natural product discovery. To complement research on the metabolic potential of microbes in extreme environments, comparative analyses on functional annotations of biosynthetic gene clusters (BGCs) from uncultivated bacterial genomes were carried out. Results indicated that in general, BGCs encode for structurally unique metabolites and that metagenome-assembled genomes did not show an obvious relationship between the composition of their core biosynthetic potential and taxonomy or geographic distribution. Nevertheless, some members of Acidobacteriota showed a phylogenetic relationship with specific metabolic traits and a few members of Proteobacteria and Desulfobacterota exhibited niche adaptations. Our results emphasize that studying specialized metabolism in environmental samples may significantly contribute to the elucidation of structures, activities, and ecological roles of microbial molecules.

Segatella clades adopt distinct roles within a single individual's gut

Segatella is a prevalent genus within individuals' gut microbiomes worldwide, especially in non-Western populations. Although metagenomic assembly and genome isolation have shed light on its genetic diversity, the lack of available isolates from this genus has resulted in a limited understanding of how members' genetic diversity translates into phenotypic diversity. Within the confines of a single gut microbiome, we have isolated 63 strains from diverse lineages of Segatella. We performed comparative analyses that exposed differences in cellular morphologies, preferences in polysaccharide utilization, yield of short-chain fatty acids, and antibiotic resistance across isolates. We further show that exposure to Segatella isolates either evokes strong or muted transcriptional responses in human intestinal epithelial cells. Our study exposes large phenotypic differences within related Segatella isolates, extending this to host-microbe interactions.

Synergy between Arbuscular Mycorrhizal Fungi and Rhizosphere Bacterial Communities Increases the Utilization of Insoluble Phosphorus and Potassium in the Soil by Maize

Arbuscular mycorrhizal (AM) fungi can enhance plant uptake of phosphorus (P) and potassium (K), but it is not yet clear whether rhizosphere bacteria can enhance the ability of AM fungi to acquire insoluble P and K from the soil. Here, pot experiments confirmed that AM fungus-promoted insoluble P and K uptake by plants requires rhizosphere bacteria. The changes of rhizosphere bacterial communities associated with AM fungi were explored by 16S rRNA amplicon sequencing and metagenomic sequencing. Five core bacteria genera identified were involved in P and K cycles. Synthetic community (SynCom) inoculation revealed that SynCom increased soil available P and K and its coinoculation with AM fungi increased P and K concentration in the plants. This study revealed that AM fungi interact with rhizosphere bacteria and promote insoluble P and K acquisition, which provided a foundation for the application of AM fungal-bacterial biofertilizers and was beneficial for the sustainable development of agriculture.

Taxonomic variation, plastic degradation, and antibiotic resistance traits of plastisphere communities in the maturation pond of a wastewater treatment plant

Wastewater treatment facilities can filter out some plastics before they reach the open environment, yet microplastics often persist throughout these systems. As they age, microplastics in wastewater may both leach and sorb pollutants and fragment to provide an increased surface area for bacterial attachment and conjugation, possibly impacting antimicrobial resistance (AMR) traits. Despite this, little is known about the effects of persistent plastic pollution on microbial functioning. To address this knowledge gap, we deployed five different artificially weathered plastic types and a glass control into the final maturation pond of a municipal wastewater treatment plant in Ōtautahi-Christchurch, Aotearoa/New Zealand. We sampled the plastic-associated biofilms (plastisphere) at 2, 6, 26, and 52 weeks, along with the ambient pond water, at three different depths (20, 40, and 60 cm from the pond water surface). We investigated the changes in plastisphere microbial diversity and functional potential through metagenomic sequencing. Bacterial 16S ribosomal RNA genes composition did not vary among plastic types and glass controls (P = 0.997) but varied among sampling times [permutational multivariate analysis of variance (PERMANOVA), P = 0.001] and depths (PERMANOVA, P = 0.011). Overall, there was no polymer-substrate specificity evident in the total composition of genes (PERMANOVA, P = 0.67), but sampling time (PERMANOVA, P = 0.002) and depth were significant factors (PERMANOVA, P = 0.001). The plastisphere housed diverse AMR gene families, potentially influenced by biofilm-meditated conjugation. The plastisphere also harbored an increased abundance of genes associated with the biodegradation of nylon, or nylon-associated substances, including nylon oligomer-degrading enzymes and hydrolases.IMPORTANCEPlastic pollution is pervasive and ubiquitous. Occurrences of plastics causing entanglement or ingestion, the leaching of toxic additives and persistent organic pollutants from environmental plastics, and their consequences for marine macrofauna are widely reported. However, little is known about the effects of persistent plastic pollution on microbial functioning. Shotgun metagenomics sequencing provides us with the necessary tools to examine broad-scale community functioning to further investigate how plastics influence microbial communities. This study provides insight into the functional consequence of continued exposure to waste plastic by comparing the prokaryotic functional potential of biofilms on five types of plastic [linear low-density polyethylene (LLDPE), nylon-6, polyethylene terephthalate, polylactic acid, and oxygen-degradable LLDPE], glass, and ambient pond water over 12 months and at different depths (20, 40, and 60 cm) within a tertiary maturation pond of a municipal wastewater treatment plant.

Efficient alleviation granular sludge floatation in a high-rate anammox reactor by dosing folate

Although granular floatation has been recognized as a significant issue hindering the application of high-rate anammox biotechnology, limited knowledge is available about its causes and control strategies. This study proposed a novel control strategy by adding folate, and demonstrated its role in the granular floatation alleviation through long-term operation and granular characterizations. It was found that the floatation of anammox granular sludge was obviously relieved with the decreased sludge floatation potential by 67.1% after dosing with folate (8 mg/L) at a high nitrogen loading rate of 12.3 kg-N/(m3·d). Physiochemical analyses showed that the decrease of extracellular polymeric substances (EPS) content (mainly protein), the alleviation of granular surface pore plugging in conjunction with the smooth discharge of generated nitrogen gas were collectively responsible for efficient floatation control. Moreover, metagenomic analysis suggested that the synergistic interactions between anammox bacteria and their symbionts were attenuated after dosing exogenous folate. Anammox bacteria would reduce their synergistic dependence on the symbionts, and decline the supply of metabolites (e.g., amino acids and carbohydrates in EPS) to symbiotic bacteria. The declined EPS excretion contributed to the alleviation of granular floatation by dredging pores blockage, thus leading to a stable system performance. The findings not only offer insights into the role of microbial interaction in granular sludge floatation, but also provide a feasible approach for controlling the floatation issue in anammox granular-based processes.

Pangolin scales as adaptations for innate immunity against pathogens

BACKGROUND

Pangolins are the only mammals that have overlapping scales covering most of their bodies, and they play a crucial role in the ecosystem, biological research, and human health and disease. Previous studies indicated pangolin scale might provide an important mechanical defense to themselves. The origin and exact functions of this unique trait remain a mystery. Using a multi-omics analysis approach, we report a novel functional explanation for how mammalian scales can provide host-pathogen defense.

RESULTS

Our data suggest that pangolin scales have a sophisticated structure that could potentially trap pathogens. We identified numerous proteins and metabolites exhibiting antimicrobial activity, which could suggest a role for scales in pathogen defense. Notably, we found evidence suggesting the presence of exosomes derived from diverse cellular origins, including mesenchymal stem cells, immune cells, and keratinocytes. This observation suggests a complex interplay where various cell types may contribute to the release of exosomes and antimicrobial compounds at the interface between scales and viable tissue. These findings indicate that pangolin scales may serve as a multifaceted defense system, potentially contributing to innate immunity. Comparisons with human nail and hair revealed pangolin-specific proteins that were enriched in functions relating to sensing, immune responses, neutrophil degranulation, and stress responses. We demonstrated the antimicrobial activity of key pangolin scale components on pathogenic bacteria by antimicrobial assays.

CONCLUSIONS

This study identifies a potential role of pangolin scales and implicates scales, as possible determinants of pathogen defense due to their structure and contents. We indicate for the first time the presence of exosomes in pangolin scales and propose the new functions of scales and their mechanisms. This new mechanism could have implications for multiple fields, including providing interesting new research directions and important insights that can be useful for synthesizing and implementing new biomimetic antimicrobial approaches.

Slaughtering processes impact microbial communities and antimicrobial resistance genes of pig carcasses

Several steps in the abattoir can influence the presence of microbes and associated resistance genes (ARGs) on the animal carcasses used for further meat processing. We investigated how these processes influence the resistome-microbiome of groups of pigs with different on-farm antimicrobial exposure status, from the moment they entered the abattoir until the end of carcass processing. Using a targeted enrichment metagenomic approach, we identified 672 unique ARGs conferring resistance to 43 distinct AMR classes from pooled skin (N = 42) and carcass swabs (N = 63) collected sequentially before, during, and after the slaughter process and food safety interventions. We observed significant variations in the resistome and microbial profiles of pigs before and after slaughter, as well as a significant decline in ARG counts, diversity, and microbial DNA load during slaughter and carcass processing, irrespective of prior antimicrobial treatments on the farm. These results suggest that existing interventions in the abattoir are effective in reducing not only the pathogen load but also the overall bacterial burden, including ARGs on pork carcasses. Concomitant with reductions in microbial and ARG counts, we observed an increase in the relative abundance of non-drug-specific ARGs, such as those conferring resistance to metals and biocides, and in particular mercury. Using a strict colocalization procedure, we found that most mercury ARGs were associated with genomes from the Pseudomonadaceae and Enterobacteriaceae families. Collectively, these findings demonstrate that slaughter and processing practices within the abattoir can shape the microbial and ARG profiles of pork carcasses during the transition from living muscle to meat.

HyLight: Strain aware assembly of low coverage metagenomes

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

Characterization of larval gut microbiota of two endoparasitoid wasps associated with their common host, Plutella xylostella (Linnaeus) (Lepidoptera: Plutellidae)

Insect gut microbes play important roles in digestion, metabolism, development, and environmental adaptation. Parasitoid wasps are one of the most important biological control agents in pest control, while the gut microbial species compositions and the associated functions have been poorly investigated. Two endoparasitoid wasps, Cotesia vestalis and Diadromus collaris, parasitize the larval stage and pupal stage of the diamondback moth, Plutella xylostella, respectively. Using whole-genome shotgun metagenomic sequencing, we characterized the gut microbial composition, diversity, and potential functional roles associated with the two parasitoid wasp larvae. The results reveal that Proteobacteria and Firmicutes are the dominant phyla in the gut of C. vestalis and D. collaris larvae, with Rhizobium and Enterococcus being the dominant genera. The putative microbial functions associated with the two parasitoid wasps might play a virtual role in assisting in consuming the host's nutritional composition. The enriched CAZymes family genes are primarily involved in the degradation and synthesis of chitin. Despite the richness of microbial species and communities, the microbes species and the microbial community structure exhibit significant similarity between the two parasitoid wasps and between the parasitoid wasp and the host P. xylostella. Notably, the prevalence of the genus Enterococcus shared among them suggests a possible link of gut microbes between the host and their associated parasitoids. Our study offers insights into the gut microbe-based interactions between the host and parasitoid wasps for the first time, potentially paving the way for the development of an ecologically friendly biocontrol strategy against the pest P. xylostella.IMPORTANCEEndoparasitoid wasps spend the majority of their lifespan within their host and heavily rely on the host's nutrition for survival. There is limited understanding regarding the composition and physiological impacts of gut microbial communities in parasitoid wasps, particularly during the larval stage, which is directly linked to the host. Based on a thorough characterization of the gut microbe and comprehensive comparative analysis, we found the microbial species of the larval parasitoid wasp Cotesia vestalis and the pupal parasitoid wasp Diadromus collaris were similar, sharing 159 genera and 277 species, as were the microbial community structure. Certain of the dominant microbial strains of the two parasitoid wasps were similar to that of their host Plutella xylostella larvae, revealing host insect may affect the microbial community of the parasitoid wasps. The putative microbial functions associated with the parasitoid wasp larvae play an important role in dietary consumption.

Simulation of 69 microbial communities indicates sequencing depth and false positives are major drivers of bias in prokaryotic metagenome-assembled genome recovery

We hypothesize that sample species abundance, sequencing depth, and taxonomic relatedness influence the recovery of metagenome-assembled genomes (MAGs). To test this hypothesis, we assessed MAG recovery in three in silico microbial communities composed of 42 species with the same richness but different sample species abundance, sequencing depth, and taxonomic distribution profiles using three different pipelines for MAG recovery. The pipeline developed by Parks and colleagues (8K) generated the highest number of MAGs and the lowest number of true positives per community profile. The pipeline by Karst and colleagues (DT) showed the most accurate results (~ 92%), outperforming the 8K and Multi-Metagenome pipeline (MM) developed by Albertsen and collaborators. Sequencing depth influenced the accurate recovery of genomes when using the 8K and MM, even with contrasting patterns: the MM pipeline recovered more MAGs found in the original communities when employing sequencing depths up to 60 million reads, while the 8K recovered more true positives in communities sequenced above 60 million reads. DT showed the best species recovery from the same genus, even though close-related species have a low recovery rate in all pipelines. Our results highlight that more bins do not translate to the actual community composition and that sequencing depth plays a role in MAG recovery and increased community resolution. Even low MAG recovery error rates can significantly impact biological inferences. Our data indicates that the scientific community should curate their findings from MAG recovery, especially when asserting novel species or metabolic traits.

Exposure to doxycycline increases risk of carrying a broad range of enteric antimicrobial resistance determinants in an elderly cohort

OBJECTIVES

High rates of antibiotic prescription in residential aged care are likely to promote enteric carriage of antibiotic-resistant pathogens and increase the risk of antibiotic treatment failure. Despite their importance, relationships between antibiotic exposures and patterns of enteric resistance carriage in this population remain poorly understood.

METHODS

We conducted a cross-sectional metagenomic cohort analysis of stool samples from residents of five long-term aged-care facilities in South Australia. Taxonomic composition was determined, and enteric carriage of antibiotic resistance genes (ARGs) was identified and quantified against the Comprehensive Antibiotic Resistance Database. Both the detection and abundance of stool taxa and ARGs were related to antibiotic exposures up to 12 months prior. Factors associated with the abundance of ARGs of high clinical concern were identified.

RESULTS

Stool samples were provided by 164 participants (median age: 88 years, IQR 81-93; 72% female). Sixty-one percent (n = 100) of participants were prescribed antibiotics at least once in the prior 12 months (median prescriptions: 4, range: 1-52), most commonly a penicillin (n = 55, 33.5%), cephalosporin (n = 53, 32.3%), diaminopyrimidine (trimethoprim) (n = 36, 22%), or tetracycline (doxycycline) (n = 21, 12.8%). More than 1100 unique ARGs, conferring resistance to 38 antibiotic classes, were identified, including 20 ARGs of high clinical concern. Multivariate logistic regression showed doxycycline exposure to be the greatest risk factor for high ARG abundance (adjusted odds ratio [aOR]=14.8, q<0.001) and a significant contributor to inter-class selection, particularly for ARGs relating to penicillins (aOR=3.1, q=0.0004) and cephalosporins (aOR=3.4, q=0.003). High enteric ARG abundance was associated with the number of separate antibiotic exposures (aOR: 6.4, q<0.001), exposures within the prior 30 days (aOR: 4.6, q=0.008) and prior 30-100 days (aOR: 2.6, q=0.008), high duration of antibiotic exposure (aOR: 7.9, q<0.001), and exposure to 3 or more antibiotic classes (aOR: 7.4, q<0.001). Carriage of one or more ARGs of high clinical concern was identified in 99% of participants (n = 162, median: 3, IQR: 2-4), involving 11 ARGs conferring resistance to aminoglycosides, four to beta-lactams, one to glycopeptides, three to fluoroquinolones, and one to oxazolidinones. Carriage of ARGs of high clinical concern was positively associated with exposure to doxycycline (aminoglycoside, fluoroquinolone, and oxazolidinone ARGs) and trimethoprim (fluoroquinolone and beta-lactam ARGs). Analysis of doxycycline impact on microbiota composition suggested that observed resistome changes arose principally through direct ARG selection, rather than through the antibiotic depletion of sensitive bacterial populations.

CONCLUSIONS

The gut microbiome of aged care residents is a major reservoir of antibiotic resistance. As a critical antibiotic in medical practice, a comprehensive understanding of the impact of doxycycline exposure on the gut resistome is paramount for informed antibiotic use, particularly in an evolving landscape of prophylactic applications. Near-universal asymptomatic carriage of clinically critical resistance determinants is highly concerning and reinforces the urgent need for improved management of antibiotic use in long-term aged care.

Comparing microbial communities in mucilage and seawater samples: Metagenomic insights into mucilage formation in the Marmara Sea

Marine environments are subject to various naturally occurring phenomena, including marine snow and mucilage. In 2021, the rapid emergence of mucilage in the Marmara Sea raised concerns about its environmental impact. This study investigates the microbial communities in mucilage and seawater samples from the Marmara Sea using metagenomic-scale comparative analyses. The results indicate significant differences in microbial composition and diversity, with mucilage samples showing higher levels of polysaccharide biosynthesis-related enzymes. Over 50% of reads in mucilage samples remained unclassified (dark matter), highlighting unknown microbial taxa. Clean seawater was characterized by a higher presence of Euryarchaeota, Proteobacteria, and Rhodothermaeota, while Chlamydiae and Fusobacteria were dominant in mucilage. The study underscores the necessity for comprehensive metagenomic analyses to understand microbial roles in mucilage formation and persistence. Early detection of microbial shifts could serve as a warning system for mucilage outbreaks, aiding in the development of management strategies.

Disclosing the virome of Aedes, Anopheles and Culex female mosquitoes, Alto Pantanal of Mato Grosso, Brazil, 2019

Using Illumina NextSeq sequencing and bioinformatics, we identified and characterized thirty-three viral sequences of unsegmented and multipartite viral families in Aedes spp., Culex sp. and Anopheles darlingi female mosquito pools from Porto São Luiz and Pirizal, Alto Pantanal. Seventeen sequences belong to unsegmented viral families, twelve represent putative novel insect-specific viruses (ISVs) within families Chuviridae (3/33; partial genomes) and coding-complete sequences of Xinmoviridae (1/33), Rhabdoviridae (2/33) and Metaviridae (6/33); and five coding-complete sequences of already-known ISVs. Notably, two putative novel rhabdoviruses, Corixo rhabdovirus 1 and 2, were phylogenetically related to Coxipo dielmovirus, but separated from other Alpharhabdovirinae genera, sharing Anopheles spp. as host. Regarding multipartite families, sixteen segments of different putative novel viruses were identified (13 coding-complete segments) within Durnavirales (4/33), Elliovirales (1/33), Hareavirales (3/33) and Reovirales (8/33) orders. Overall, this study describes twenty-eight (28/33) putative novel ISVs and five (5/33) already described viruses using metagenomics approach.

Mitochondrial Phylogenomics of Scoliidae from China, with Evidence to Challenge the Former Placement of the Colpa Group

Scoliidae, also known as scarab hunters or flower wasps, are important in the biological control of scarabs and for pollination. Mitogenomic and phylogenetic studies are rare for this group. In this study, 10 mitochondrial genomes representing eight genera in two tribes of the family Scoliidae were determined. The general features and rearrangements of the mitochondrial genomes for 15 Scoliidae species representing all genera distributed in China were described and compared and the phylogenetic relationships among them were inferred using MrBayes and IQtree based on four data matrices. Most sequences of Scoliidae have one extra trnM gene. Species belonging to Campsomerini have lower A + T content than all Scoliini species except for Colpa tartara in this study. The AT-skew is positive in 7 out of 15 species. All 15 Scoliidae sequences have similar conserved gene arrangements with the same arrangements of PCGs and rRNA genes, except for Campsomeriella annulata. The tRNA genes have the highest frequency of rearrangement, and C. tartara is always rearranged as in its Scoliini counterparts. Our phylogenetic results support most of the relationships between genera and tribes of Scoliidae in former morphological studies. However, Colpa tartara is proved to be closer to Scoliini according to genome features, phylogenetic analyses and some morphological evidence, which challenges the former attribution of the Colpa group.

Four complete mitochondrial genomes of the subgenus Pterelachisus (Diptera, Tipulidae, Tipula) and implications for the higher phylogeny of the family Tipulidae

The complete mitochondrial genomes of Tipula (Pterelachisus) cinereocincta mesacantha Alexander, 1934, T. (P.) legalis Alexander, 1933, T. (P.) varipennis Meigen, 1818, and T. (P.) yasumatsuana Alexander, 1954 are reported, three of them being sequenced for the first time. The mitochondrial genome lengths of the four species are 15,907 bp, 15,625 bp, 15,772 bp, and 15,735 bp, respectively. All genomes exhibit a high AT base composition, with A + T content of 76.7%, 75.0%, 77.8%, and 75.4%, respectively. The newly reported mitogenomes herein show a general similarity in overall structure, gene order, base composition, and nucleotide content to those of the previously studied species within the family Tipulidae. Phylogenetic analyses were conducted to investigate the relationships within Tipulidae, using both Maximum Likelihood and Bayesian Inference approaches. The results show that the four target species of the subgenus T. (Pterelachisus) basically form a monophyletic group within Tipulidae, clustering with species of the Tipula subgenera T. (Lunatipula), T. (Vestiplex), and T. (Formotipula); however, the genus Tipula is not monophyletic. Moreover, neither the tipulid subfamily Tipulinae nor the family Limoniidae is supported to be a monophyletic group. The monophyly of the family Tipulidae, and the sister relationship between Tipulidae and Cylindrotomidae are reconfirmed. These research findings could contribute to deep insights into the systematic and evolutionary patterns of crane flies.

Multi-View Integrative Approach For Imputing Short-Chain Fatty Acids and Identifying Key factors predicting Blood SCFA

Short-chain fatty acids (SCFAs) are the main metabolites produced by bacterial fermentation of dietary fiber within gastrointestinal tract. SCFAs produced by gut microbiotas (GMs) are absorbed by host, reach bloodstream, and are distributed to different organs, thus influencing host physiology. However, due to the limited budget or the poor sensitivity of instruments, most studies on GMs have incomplete blood SCFA data, limiting our understanding of the metabolic processes within the host. To address this gap, we developed an innovative multi-task multi-view integrative approach (M2AE, Multi-task Multi-View Attentive Encoders), to impute blood SCFA levels using gut metagenomic sequencing (MGS) data, while taking into account the intricate interplay among the gut microbiome, dietary features, and host characteristics, as well as the nuanced nature of SCFA dynamics within the body. Here, each view represents a distinct type of data input (i.e., gut microbiome compositions, dietary features, or host characteristics). Our method jointly explores both view-specific representations and cross-view correlations for effective predictions of SCFAs. We applied M2AE to two in-house datasets, which both include MGS and blood SCFAs profiles, host characteristics, and dietary features from 964 subjects and 171 subjects, respectively. Results from both of two datasets demonstrated that M2AE outperforms traditional regression-based and neural-network based approaches in imputing blood SCFAs. Furthermore, a series of gut bacterial species (e.g., Bacteroides thetaiotaomicron and Clostridium asparagiforme), host characteristics (e.g., race, gender), as well as dietary features (e.g., intake of fruits, pickles) were shown to contribute greatly to imputation of blood SCFAs. These findings demonstrated that GMs, dietary features and host characteristics might contribute to the complex biological processes involved in blood SCFA productions. These might pave the way for a deeper and more nuanced comprehension of how these factors impact human health.

An animal charcoal contaminated cottage industry soil highlighted by halophilic archaea dominance and decimation of bacteria

An animal charcoal contaminated cottage industry soil in Lagos, Nigeria (ACGT) was compared in an ex post facto study with a nearby unimpacted soil (ACGC). Hydrocarbon content was higher than regulatory limits in ACGT (180.2 mg/kg) but lower in ACGC (19.28 mg/kg). Heavy metals like nickel, cadmium, chromium and lead were below detection limit in ACGC. However, all these metals, except cadmium, were detected in ACGT, but at concentrations below regulatory limits. Furthermore, copper (253.205 mg/kg) and zinc (422.630 mg/kg) were above regulatory limits in ACGT. Next generation sequencing revealed that the procaryotic community was dominated by bacteria in ACGC (62%) while in ACGT archaea dominated (76%). Dominant phyla in ACGC were Euryarchaeota (37%), Pseudomonadota (16%) and Actinomycetota (12%). In ACGT it was Euryarchaeota (76%), Bacillota (9%), Pseudomonadota (7%) and Candidatus Nanohaloarchaeota (5%). Dominant Halobacteria genera in ACGT were Halobacterium (16%), Halorientalis (16%), unranked halophilic archaeon (13%) Salarchaeum (6%) and Candidatus Nanohalobium (5%), whereas ACGC showed greater diversity dominated by bacterial genera Salimicrobium (7%) and Halomonas (3%). Heavy metals homeostasis genes, especially for copper, were fairly represented in both soils but with bacterial taxonomic affiliations. Sites like ACGT, hitherto poorly studied and understood, could be sources of novel bioresources.

Protocol for the construction and functional profiling of metagenome-assembled genomes for microbiome analyses

Constructing metagenome-assembled genomes (MAGs) from complex metagenomic samples involves a series of bioinformatics operations, each requiring deep bioinformatics knowledge. Here, we present a protocol for constructing MAGs and conducting functional profiling to address biological questions. We describe steps for system configuration, data downloads, read processing, removal of human DNA contamination, metagenomic assembly, and statistical quality assessment of the final assembly. Additionally, we detail procedures for the construction and refinement of MAGs, as well as the functional profiling of MAGs.

New evidence supports the prophage origin of RcGTA

Gene transfer agents (GTAs) are phage-like entities that package and transfer random host genome fragments between prokaryotes. RcGTA, produced by Rhodobacter capsulatus, is hypothesized to originate from a prophage ancestor. Most of the evidence supporting this hypothesis came from the finding of RcGTA-like genes in phages. More than 75% of the RcGTA genes have a phage homolog. However, only a few RcGTA homologs have been identified in a (pro)phage genome, leaving the hypothesis that GTAs evolved from prophages through gene loss with only weak evidence. We herein report the discovery of an inducible prophage (vB_MseS-P1) from a Mesorhizobium sediminum strain that contains the largest number (12) of RcGTA homologs found in a phage genome to date. We also identified three putative prophages and two prophage remnants harboring 12-14 RcGTA homologs in a Methylobacterium nodulans strain. The protein remote homology detection also revealed more RcGTA homologs from other phages than we previously thought. Moreover, the head-tail gene architecture of these newly discovered prophage-related elements closely resembles that of RcGTA. Furthermore, vB_MseS-P1 virions have structural proteins similar to RcGTA particles. Close phylogenetic relationships between certain prophage genes and RcGTA-like genes in Alphaproteobacteria further support the shared ancestry between RcGTA and prophages. Our findings provide new relatively direct evidence of the origin of RcGTA from a prophage progenitor.IMPORTANCEGTAs are important genetic elements in certain groups of bacteria and contribute to the genetic diversification, evolution, and ecological adaptation of bacteria. RcGTA, a common type of GTA, is known to package and transfer random fragments of the bacterial genome to recipient cells. However, the origin of RcGTA is still elusive. It has been hypothesized that RcGTA evolved from a prophage ancestor through gene loss. However, the few RcGTA homologs identified in a (pro)phage genome leave the hypothesis lacking direct evidence. This study uncovers the presence of a large number of RcGTA homologs in an inducible prophage and several putative prophages. The similar head-tail gene architecture and structural protein compositions of these newly discovered prophage-related elements and RcGTA further demonstrate an unprecedentedly observed close evolutionary relationship between prophages and RcGTA. Together, our findings provide more direct evidence supporting the origin of RcGTA from prophage.

Metagenomic Profiling of Internationally Sourced Sewage Influents and Effluents Yields Insight into Selecting Targets for Antibiotic Resistance Monitoring

It has been debated whether wastewater treatment plants (WWTPs) primarily act to attenuate or amplify antibiotic resistance genes (ARGs). However, ARGs are highly diverse with respect to their resistance mechanisms, mobilities, and taxonomic hosts and therefore their behavior in WWTPs should not be expected to be universally conserved. We applied metagenomic sequencing to wastewater influent and effluent samples from 12 international WWTPs to classify the behavior of specific ARGs entering and exiting WWTPs. In total, 1079 different ARGs originating from a variety of bacteria were detected. This included ARGs that could be mapped to assembled scaffolds corresponding to nine human pathogens. While the relative abundance (per 16S rRNA gene) of ARGs decreased during treatment at 11 of the 12 WWTPs sampled and absolute abundance (per mL) decreased at all 12 WWTPs, increases in relative abundance were observed for 40% of the ARGs detected at the 12th WWTP. Also, the relative abundance of mobile genetic elements (MGE) increased during treatment, but the fraction of ARGs known to be transmissible between species decreased, thus demonstrating that increased MGE prevalence may not be generally indicative of an increase in ARGs. A distinct conserved resistome was documented in both influent and effluent across samples, suggesting that well-functioning WWTPs generally attenuate influent antibiotic resistance loads. This work helps inform strategies for wastewater surveillance of antibiotic resistance, highlighting the utility of tracking ARGs as indicators of treatment performance and relative risk reduction.

Full-Length Transcriptome Profiling of the Complete Mitochondrial Genome of Sericothrips houjii (Thysanoptera: Thripidae: Sericothripinae) Featuring Extensive Gene Rearrangement and Duplicated Control Regions

The mitochondrial genome (mitogenome) of Thysanoptera has extensive gene rearrangement, and some species have repeatable control regions. To investigate the characteristics of the gene expression, transcription and post-transcriptional processes in such extensively gene-rearranged mitogenomes, we sequenced the mitogenome and mitochondrial transcriptome of Sericothrips houjii to analyze. The mitogenome was 14,965 bp in length and included two CRs contains 140 bp repeats between COIII-trnN (CR1) and trnT-trnP (CR2). Unlike the putative ancestral arrangement of insects, S. houjii exhibited only six conserved gene blocks encompassing 14 genes (trnL2-COII, trnD-trnK, ND2-trnW, ATP8-ATP6, ND5-trnH-ND4-ND4L and trnV-lrRNA). A quantitative transcription map showed the gene with the highest relative expression in the mitogenome was ND4-ND4L. Based on analyses of polycistronic transcripts, non-coding RNAs (ncRNAs) and antisense transcripts, we proposed a transcriptional model of this mitogenome. Both CRs contained the transcription initiation sites (TISs) and transcription termination sites (TTSs) of both strands, and an additional TIS for the majority strand (J-strand) was found within antisense lrRNA. The post-transcriptional cleavage processes followed the "tRNA punctuation" model. After the cleavage of transfer RNAs (tRNAs), COI and ND3 matured as bicistronic mRNA COI/ND3 due to the translocation of intervening tRNAs, and the 3' untranslated region (UTR) remained in the mRNAs for COII, COIII, CYTB and ND5. Additionally, isoform RNAs of ND2, srRNA and lrRNA were identified. In summary, the relative mitochondrial gene expression levels, transcriptional model and post-transcriptional cleavage process of S. houjii are notably different from those insects with typical mitochondrial gene arrangements. In addition, the phylogenetic tree of Thripidae including S. houjii was reconstructed. Our study provides insights into the phylogenetic status of Sericothripinae and the transcriptional and post-transcriptional regulation processes of extensively gene-rearranged insect mitogenomes.

Isolation and sequencing of Orthoflavivirus ilheusense from mosquitoes collected in the Brazilian cerrado, West-Central region, Brazil

The Orthoflavivirus ilheusense (ILHV) is an arbovirus that was first isolated in Brazil in 1944 during an epidemiologic investigation of yellow fever. Is a member of the Flaviviridae family and it belongs to the antigenic complex of the Ntaya virus group. Psorophora ferox is the primary vector of ILHV and this study presents the isolation and phylogenetic analysis of ILHV in a pool of Ps. ferox collected in the state of Goiás in 2021. Viral isolation tests were performed on Vero cells and C6/36 clones. The indirect immunofluorescence test (IFI) was used to confirm the positivity of the sample. The positive sample underwent RT-qPCR, sequencing, and phylogenetic analysis. This is the first report of ILHV circulation in this municipality and presented close relationship between this isolate and another ILHV isolate collected in the city of Belém (PA).

Genome assembly in the telomere-to-telomere era

Genome sequences largely determine the biology and encode the history of an organism, and de novo assembly - the process of reconstructing the genome sequence of an organism from sequencing reads - has been a central problem in bioinformatics for four decades. Until recently, genomes were typically assembled into fragments of a few megabases at best, but now technological advances in long-read sequencing enable the near-complete assembly of each chromosome - also known as telomere-to-telomere assembly - for many organisms. Here, we review recent progress on assembly algorithms and protocols, with a focus on how to derive near-telomere-to-telomere assemblies. We also discuss the additional developments that will be required to resolve remaining assembly gaps and to assemble non-diploid genomes.

Food safety assessment and 28-day toxicity study of the synbiotic medical food consortium SBD121

The human gut microbiome plays a crucial role in immune function. The synbiotic consortium or Defined Microbial Assemblage™ (DMA™) Medical Food product, SBD121, consisting of probiotic microbes and prebiotic fibers was designed for the clinical dietary management of rheumatoid arthritis. A 28-day repeated administration study was performed to evaluate the oral toxicity of SBD121 in male and female rats (age/weight at study start: 60 days/156-264 g) administered levels of 0, 4.96 x 1010, 2.48 x 1011, or 4.96 x 1011 colony forming units (CFU)/kg-bw. No treatment related changes were observed in ophthalmological effects, mortality, morbidity, general health and clinical observations, urinalysis, hematology, serum chemistry, absolute or relative organ weights, gross necropsy, or histopathology. A significant decrease in body weight was reported in females in the low and high-concentration groups, which corresponded in part with a significant decrease in food consumption. Results of the functional observation battery indicated front grip strength was significantly greater in the high-concentration males compared to the controls; however, this effect was not considered adverse. Based on these findings, the administration of the Medical Food SBD121 to male and female rats has a no-observable adverse effect level (NOAEL) at the highest level tested of 4.96 x 1011 CFU/kg-bw.

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.