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

Autotrophic biofilms sustained by deeply sourced groundwater host diverse bacteria implicated in sulfur and hydrogen metabolism

BACKGROUND

Biofilms in sulfide-rich springs present intricate microbial communities that play pivotal roles in biogeochemical cycling. We studied chemoautotrophically based biofilms that host diverse CPR bacteria and grow in sulfide-rich springs to investigate microbial controls on biogeochemical cycling.

RESULTS

Sulfide springs biofilms were investigated using bulk geochemical analysis, genome-resolved metagenomics, and scanning transmission X-ray microscopy (STXM) at room temperature and 87 K. Chemolithotrophic sulfur-oxidizing bacteria, including Thiothrix and Beggiatoa, dominate the biofilms, which also contain CPR Gracilibacteria, Absconditabacteria, Saccharibacteria, Peregrinibacteria, Berkelbacteria, Microgenomates, and Parcubacteria. STXM imaging revealed ultra-small cells near the surfaces of filamentous bacteria that may be CPR bacterial episymbionts. STXM and NEXAFS spectroscopy at carbon K and sulfur L2,3 edges show that filamentous bacteria contain protein-encapsulated spherical elemental sulfur granules, indicating that they are sulfur oxidizers, likely Thiothrix. Berkelbacteria and Moranbacteria in the same biofilm sample are predicted to have a novel electron bifurcating group 3b [NiFe]-hydrogenase, putatively a sulfhydrogenase, potentially linked to sulfur metabolism via redox cofactors. This complex could potentially contribute to symbioses, for example, with sulfur-oxidizing bacteria such as Thiothrix that is based on cryptic sulfur cycling. One Doudnabacteria genome encodes adjacent sulfur dioxygenase and rhodanese genes that may convert thiosulfate to sulfite. We find similar conserved genomic architecture associated with CPR bacteria from other sulfur-rich subsurface ecosystems.

CONCLUSIONS

Our combined metagenomic, geochemical, spectromicroscopic, and structural bioinformatics analyses of biofilms growing in sulfide-rich springs revealed consortia that contain CPR bacteria and sulfur-oxidizing Proteobacteria, including Thiothrix, and bacteria from a new family within Beggiatoales. We infer roles for CPR bacteria in sulfur and hydrogen cycling. Video Abstract.

Unraveling the phylogenomic diversity of Methanomassiliicoccales and implications for mitigating ruminant methane emissions

BACKGROUND

Methanomassiliicoccales are a recently identified order of methanogens that are diverse across global environments particularly the gastrointestinal tracts of animals; however, their metabolic capacities are defined via a limited number of cultured strains.

RESULTS

Here, we profile and analyze 243 Methanomassiliicoccales genomes assembled from cultured representatives and uncultured metagenomes recovered from various biomes, including the gastrointestinal tracts of different animal species. Our analyses reveal the presence of numerous undefined genera and genetic variability in metabolic capabilities within Methanomassiliicoccales lineages, which is essential for adaptation to their ecological niches. In particular, gastrointestinal tract Methanomassiliicoccales demonstrate the presence of co-diversified members with their hosts over evolutionary timescales and likely originated in the natural environment. We highlight the presence of diverse clades of vitamin transporter BtuC proteins that distinguish Methanomassiliicoccales from other archaeal orders and likely provide a competitive advantage in efficiently handling B12. Furthermore, genome-centric metatranscriptomic analysis of ruminants with varying methane yields reveal elevated expression of select Methanomassiliicoccales genera in low methane animals and suggest that B12 exchanges could enable them to occupy ecological niches that possibly alter the direction of H2 utilization.

CONCLUSIONS

We provide a comprehensive and updated account of divergent Methanomassiliicoccales lineages, drawing from numerous uncultured genomes obtained from various habitats. We also highlight their unique metabolic capabilities involving B12, which could serve as promising targets for mitigating ruminant methane emissions by altering H2 flow.

metaProbiotics: a tool for mining probiotic from metagenomic binning data based on a language model

Beneficial bacteria remain largely unexplored. Lacking systematic methods, understanding probiotic community traits becomes challenging, leading to various conclusions about their probiotic effects among different publications. We developed language model-based metaProbiotics to rapidly detect probiotic bins from metagenomes, demonstrating superior performance in simulated benchmark datasets. Testing on gut metagenomes from probiotic-treated individuals, it revealed the probioticity of intervention strains-derived bins and other probiotic-associated bins beyond the training data, such as a plasmid-like bin. Analyses of these bins revealed various probiotic mechanisms and bai operon as probiotic Ruminococcaceae's potential marker. In different health-disease cohorts, these bins were more common in healthy individuals, signifying their probiotic role, but relevant health predictions based on the abundance profiles of these bins faced cross-disease challenges. To better understand the heterogeneous nature of probiotics, we used metaProbiotics to construct a comprehensive probiotic genome set from global gut metagenomic data. Module analysis of this set shows that diseased individuals often lack certain probiotic gene modules, with significant variation of the missing modules across different diseases. Additionally, different gene modules on the same probiotic have heterogeneous effects on various diseases. We thus believe that gene function integrity of the probiotic community is more crucial in maintaining gut homeostasis than merely increasing specific gene abundance, and adding probiotics indiscriminately might not boost health. We expect that the innovative language model-based metaProbiotics tool will promote novel probiotic discovery using large-scale metagenomic data and facilitate systematic research on bacterial probiotic effects. The metaProbiotics program can be freely downloaded at https://github.com/zhenchengfang/metaProbiotics.

Responses of drinking water bulk and biofilm microbiota to elevated water age in bench-scale simulated distribution systems

Reductions in nonresidential water demand during the COVID-19 pandemic highlighted the importance of understanding how water age impacts drinking water quality and microbiota in piped distribution systems. Using benchtop model distribution systems, we aimed to characterize the impacts of elevated water age on microbiota in bulk water and pipe wall biofilms. Five replicate constant-flow reactors were fed with municipal chloraminated tap water for 6 months prior to building closures and 7 months after. After building closures, chloramine levels entering the reactors dropped; in the reactor bulk water and biofilms the mean cell counts and ATP concentrations increased over an order of magnitude while the detection of opportunistic pathogens remained low. Water age, and the corresponding physicochemical changes, strongly influenced microbial abundance and community composition. Differential initial microbial colonization also had a lasting influence on microbial communities in each reactor (i.e., historical contingency).

Microbial organic fertilizer prepared by co-composting of Trichoderma dregs mitigates dissemination of resistance, virulence genes, and bacterial pathogens in soil and rhizosphere

The use of manure, mycelium dregs and other waste as organic fertilizer is the main source of antibiotic resistance genes (ARGs) and pathogens in farmland. Composting of waste may effectively remove ARGs and pathogens. However, the profiles and drivers of changes in metal resistance genes (MRGs), biocide resistance genes (BRGs), and virulence genes (VGs) in soil-crop rhizosphere systems after compost application remain largely unknown. Here, we prepared two kinds of microbial organic fertilizers (MOF) by using Trichoderma dregs (TDs) and organic fertilizer mixing method (MOF1) and TDs co-composting method (MOF2). The effects of different types and doses of MOF on resistance genes, VGs and pathogens in soil-rhizosphere system and their potential mechanisms were studied. The results showed that co-composting of TDs promoted the decomposition of organic carbon and decreased the absolute abundance of ARGs and mobile genetic elements (MGEs) by 53.4-65.0%. MOF1 application significantly increased the abundance and diversity of soil ARGs, BRGs, and VGs, while low and medium doses of MOF2 significantly decreased their abundance and diversity in soil and rhizosphere. Patterns of positive co-occurrence between MGEs and VGs/MRGs/BRGs/ARGs were observed through statistical analysis and gene arrangements. ARGs/MRGs reductions in MOF2 soil were directly driven by weakened horizontal gene transfer triggered by MGEs. Furthermore, MOF2 reduced soil BRGs/VGs levels by shifting bacterial communities (e.g., reduced bacterial host) or improving soil property. Our study provided new insights into the rational use of waste to minimize the spread of resistomes and VGs in soil.

Infant gut DNA bacteriophage strain persistence during the first 3 years of life

Bacteriophages are key components of gut microbiomes, yet the phage colonization process in the infant gut remains uncertain. Here, we establish a large phage sequence database and use strain-resolved analyses to investigate DNA phage succession in infants throughout the first 3 years of life. Analysis of 819 fecal metagenomes collected from 28 full-term and 24 preterm infants and their mothers revealed that early-life phageome richness increases over time and reaches adult-like complexity by age 3. Approximately 9% of early phage colonizers, which are mostly maternally transmitted and infect Bacteroides, persist for 3 years and are more prevalent in full-term than in preterm infants. Although rare, phages with stop codon reassignment are more likely to persist than non-recoded phages and generally display an increase in in-frame reassigned stop codons over 3 years. Overall, maternal seeding, stop codon reassignment, host CRISPR-Cas locus prevalence, and diverse phage populations contribute to stable viral colonization.

Metagenomics reveals the divergence of gut microbiome composition and function in two common pika species (Ochotona curzoniae and Ochotona daurica) in China

Gut microbiome plays crucial roles in animal adaptation and evolution. However, research on adaptation and evolution of small wild high-altitude mammals from the perspective of gut microbiome is still limited. In this study, we compared differences in intestinal microbiota composition and function in Plateau pikas (Ochotona curzoniae) and Daurian pikas (O. daurica) using metagenomic sequencing. Our results showed that microbial community structure had distinct differences in different pika species. Prevotella, Methanosarcina, Rhizophagus, and Podoviridae were abundant bacteria, archaea, eukaryotes, and viruses in Plateau pikas, respectively. However, Prevotella, Methanosarcina, Ustilago, and Retroviridae were dominated in Daurian pikas. Functional pathways related to carbohydrate metabolism that refer to the utilization of pectin, hemicellulose, and debranching enzymes were abundant in Plateau pikas, while the function for degradation of chitin, lignin, and cellulose was more concentrated in Daurian pikas. Pika gut had abundant multidrug resistance genes, followed by glycopeptide and beta-lactamase resistance genes, as well as high-risk antibiotic resistance genes, such as mepA, tetM, and bacA. Escherichia coli and Klebsiella pneumoniae may be potential hosts of mepA. This research provided new insights for adaptation and evolution of wild animals from perspective of gut microbiome and broadened our understanding of high-risk antibiotic resistance genes and potential pathogens of wild animals.

An Improved Chromosome-Level Genome Assembly of the Firefly Pyrocoelia pectoralis

The endemic and endangered Chinese firefly Pyrocoelia pectoralis is a sexually dimorphic, nocturnal species. A previous attempt by this team to assemble a draft genome of P. pectoralis using PacBio and Illumina HiSeq X Ten platforms was limited in its usefulness by high redundancy and contamination. This prompted us to conduct an improved chromosome-level genome assembly of P. pectoralis. Ten chromosomes were further assembled based on Hi-C data to a 532.25 Mb final size with a 52.87 Mb scaffold N50. The total repeat lengths in the genome of P. pectoralis amount to 227.69 Mb; 42.78%. In total, 12,789 genes could be functionally annotated using at least one public database. Phylogenetic inference indicated that P. pectoralis and P. pyralis diverged ~51.41 million years ago. Gene family expansion and contraction analysis of 12 species were performed, and 546 expanded and 2660 contracted gene families were identified in P. pectoralis. We generated a high-quality draft of the P. pectoralis genome. This genome assembly should help promote research on the species' sexual dimorphism and its unique courtship behavior, which involves a combination of pheromonal and bioluminescent signals. It also can serve as a resource for accelerating genome-assisted improvements in the conservation of this species.

Niche differentiation within bacterial key-taxa in stratified surface waters of the Southern Pacific Gyre

One of the most hostile marine habitats on Earth is the surface of the South Pacific Gyre (SPG), characterized by high solar radiation, extreme nutrient depletion, and low productivity. During the SO-245 "UltraPac" cruise through the center of the ultra-oligotrophic SPG, the marine alphaproteobacterial group AEGEAN169 was detected by fluorescence in situ hybridization at relative abundances up to 6% of the total microbial community in the uppermost water layer, with two distinct populations (Candidatus Nemonibacter and Ca. Indicimonas). The high frequency of dividing cells combined with high transcript levels suggests that both clades may be highly metabolically active. Comparative metagenomic and metatranscriptomic analyses of AEGEAN169 revealed that they encoded subtle but distinct metabolic adaptions to this extreme environment in comparison to their competitors SAR11, SAR86, SAR116, and Prochlorococcus. Both AEGEAN169 clades had the highest percentage of transporters per predicted proteins (9.5% and 10.6%, respectively). In particular, the high expression of ABC transporters in combination with proteorhodopsins and the catabolic pathways detected suggest a potential scavenging lifestyle for both AEGEAN169 clades. Although both AEGEAN169 clades may share the genomic potential to utilize phosphonates as a phosphorus source, they differ in their metabolic pathways for carbon and nitrogen. Ca. Nemonibacter potentially use glycine-betaine, whereas Ca. Indicimonas may catabolize urea, creatine, and fucose. In conclusion, the different potential metabolic strategies of both clades suggest that both are well adapted to thrive resource-limited conditions and compete well with other dominant microbial clades in the uppermost layers of SPG surface waters.

Experimental evolution at ecological scales allows linking of viral genotypes to specific host strains

Viruses shape microbial community structure and activity through the control of population diversity and cell abundances. Identifying and monitoring the dynamics of specific virus-host pairs in nature is hampered by the limitations of culture-independent approaches such as metagenomics, which do not always provide strain-level resolution, and culture-based analyses, which eliminate the ecological background and in-situ interactions. Here, we have explored the interaction of a specific "autochthonous" host strain and its viruses within a natural community. Bacterium Salinibacter ruber strain M8 was spiked into its environment of isolation, a crystallizer pond from a coastal saltern, and the viral and cellular communities were monitored for one month using culture, metagenomics, and microscopy. Metagenome sequencing indicated that the M8 abundance decreased sharply after being added to the pond, likely due to forces other than viral predation. However, the presence of M8 selected for two species of a new viral genus, Phoenicisalinivirus, for which 120 strains were isolated. During this experiment, an assemblage of closely related viral genomic variants was replaced by a single population with the ability to infect M8, a scenario which was compatible with the selection of a genomic variant from the rare biosphere. Further analysis implicated a viral genomic region putatively coding for a tail fiber protein to be responsible for M8 specificity. Our results indicate that low abundance viral genotypes provide a viral seed bank that allows for a highly specialized virus-host response within a complex ecological background.

Shotgun Metagenomics Reveals Minor Micro"bee"omes Diversity Defining Differences between Larvae and Pupae Brood Combs

Bees represent not only a valuable asset in agriculture, but also serve as a model organism within contemporary microbiology. The metagenomic composition of the bee superorganism has been substantially characterized. Nevertheless, traditional cultural methods served as the approach to studying brood combs in the past. Indeed, the comb microbiome may contribute to determining larval caste differentiation and hive immunity. To further this understanding, we conducted a shotgun sequencing analysis of the brood comb microbiome. While we found certain similarities regarding species diversity, it exhibits significant differentiation from all previously described hive metagenomes. Many microbiome members maintain a relatively constant ratio, yet taxa with the highest abundance level tend to be ephemeral. More than 90% of classified metagenomes were Gammaproteobacteria, Bacilli and Actinobacteria genetic signatures. Jaccard dissimilarity between samples based on bacteria genus classifications hesitate from 0.63 to 0.77, which for shotgun sequencing indicates a high consistency in bacterial composition. Concurrently, we identified antagonistic relationships between certain bacterial clusters. The presence of genes related to antibiotic synthesis and antibiotic resistance suggests potential mechanisms underlying the stability of comb microbiomes. Differences between pupal and larval combs emerge in the total metagenome, while taxa with the highest abundance remained consistent. All this suggests that a key role in the functioning of the comb microbiome is played by minor biodiversity, the function of which remains to be established experimentally.

Metagenomic assembly is the main bottleneck in the identification of mobile genetic elements

Antimicrobial resistance genes (ARG) are commonly found on acquired mobile genetic elements (MGEs) such as plasmids or transposons. Understanding the spread of resistance genes associated with mobile elements (mARGs) across different hosts and environments requires linking ARGs to the existing mobile reservoir within bacterial communities. However, reconstructing mARGs in metagenomic data from diverse ecosystems poses computational challenges, including genome fragment reconstruction (assembly), high-throughput annotation of MGEs, and identification of their association with ARGs. Recently, several bioinformatics tools have been developed to identify assembled fragments of plasmids, phages, and insertion sequence (IS) elements in metagenomic data. These methods can help in understanding the dissemination of mARGs. To streamline the process of identifying mARGs in multiple samples, we combined these tools in an automated high-throughput open-source pipeline, MetaMobilePicker, that identifies ARGs associated with plasmids, IS elements and phages, starting from short metagenomic sequencing reads. This pipeline was used to identify these three elements on a simplified simulated metagenome dataset, comprising whole genome sequences from seven clinically relevant bacterial species containing 55 ARGs, nine plasmids and five phages. The results demonstrated moderate precision for the identification of plasmids (0.57) and phages (0.71), and moderate sensitivity of identification of IS elements (0.58) and ARGs (0.70). In this study, we aim to assess the main causes of this moderate performance of the MGE prediction tools in a comprehensive manner. We conducted a systematic benchmark, considering metagenomic read coverage, contig length cutoffs and investigating the performance of the classification algorithms. Our analysis revealed that the metagenomic assembly process is the primary bottleneck when linking ARGs to identified MGEs in short-read metagenomics sequencing experiments rather than ARGs and MGEs identification by the different tools.

An Overview of Aquatica Fu et al., a Phylogeny of Aquatic Fireflies Using Mitochondrial Genomes, a Description of Two New Species, and a New Record of Aquatic Fireflies in China (Coleoptera: Lampyridae: Luciolinae)

There are now seven species in the genus Aquatica Fu & Ballantyne, with all but one known from the characteristics of males, females, and larvae. Molecular information is combined with morphological taxonomy for the delimitation of both genus and species. The monophyly of the genus Aquatica is strongly supported in six trees, and its position as sister to the genus Nipponoluciola Ballantyne Kawashima Jusoh et Suzuki is stable across maximum likelihood and Bayesian inference results. Two new species of Aquatica Fu & Ballantyne, A. qingshen sp. nov. and A. xianning sp. nov., described from the features of males, females, and larvae, are closely related within a single clade within the genus Aquatica. Females have distinctively shaped median oviduct plates. A Chinese population identified as A. lateralis was found to be morphologically similar to the Japanese population, but genetic distances suggest that it is a distinctive species. No larvae are associated with this species. Definitions of the aquatic status of Luciolinae fireflies are expanded.

The Role of the Gut Microbiome in Hematological Cancers

Humans are in a complex symbiotic relationship with a wide range of microbial organisms, including bacteria, viruses, and fungi. The evolution and composition of the human microbiome can be an indicator of how it may affect human health and susceptibility to diseases. Microbiome alteration, termed as dysbiosis, has been linked to the pathogenesis and progression of hematological cancers. A variety of mechanisms, including epithelial barrier disruption, local chronic inflammation response trigger, antigen dis-sequestration, and molecular mimicry, have been proposed to be associated with gut microbiota. Dysbiosis may be induced or worsened by cancer therapies (such as chemotherapy and/or hematopoietic stem cell transplantation) or infection. The use of antibiotics during treatment may also promote dysbiosis, with possible long-term consequences. The aim of this review is to provide a succinct summary of the current knowledge describing the role of the microbiome in hematological cancers, as well as its influence on their therapies. Modulation of the gut microbiome, involving modifying the composition of the beneficial microorganisms in the management and treatment of hematological cancers is also discussed. Additionally discussed are the latest developments in modeling approaches and tools used for computational analyses, interpretation and better understanding of the gut microbiome data.

Hecatomb: an integrated software platform for viral metagenomics

BACKGROUND

Modern sequencing technologies offer extraordinary opportunities for virus discovery and virome analysis. Annotation of viral sequences from metagenomic data requires a complex series of steps to ensure accurate annotation of individual reads and assembled contigs. In addition, varying study designs will require project-specific statistical analyses.

FINDINGS

Here we introduce Hecatomb, a bioinformatic platform coordinating commonly used tasks required for virome analysis. Hecatomb means "a great sacrifice." In this setting, Hecatomb is "sacrificing" false-positive viral annotations using extensive quality control and tiered-database searches. Hecatomb processes metagenomic data obtained from both short- and long-read sequencing technologies, providing annotations to individual sequences and assembled contigs. Results are provided in commonly used data formats useful for downstream analysis. Here we demonstrate the functionality of Hecatomb through the reanalysis of a primate enteric and a novel coral reef virome.

CONCLUSION

Hecatomb provides an integrated platform to manage many commonly used steps for virome characterization, including rigorous quality control, host removal, and both read- and contig-based analysis. Each step is managed using the Snakemake workflow manager with dependency management using Conda. Hecatomb outputs several tables properly formatted for immediate use within popular data analysis and visualization tools, enabling effective data interpretation for a variety of study designs. Hecatomb is hosted on GitHub (github.com/shandley/hecatomb) and is available for installation from Bioconda and PyPI.

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

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

Freshwater Viral Metagenome Analyses Targeting dsDNA Viruses

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

pH selects for distinct N2O-reducing microbiomes in tropical soil microcosms

Nitrous oxide (N2O), a greenhouse gas with ozone destruction potential, is mitigated by the microbial reduction to dinitrogen catalyzed by N2O reductase (NosZ). Bacteria with NosZ activity have been studied at circumneutral pH but the microbiology of low pH N2O reduction has remained elusive. Acidic (pH < 5) tropical forest soils were collected in the Luquillo Experimental Forest in Puerto Rico, and microcosms maintained with low (0.02 mM) and high (2 mM) N2O assessed N2O reduction at pH 4.5 and 7.3. All microcosms consumed N2O, with lag times of up to 7 months observed in microcosms with 2 mM N2O. Comparative metagenome analysis revealed that Rhodocyclaceae dominated in circumneutral microcosms under both N2O feeding regimes. At pH 4.5, Peptococcaceae dominated in high-N2O, and Hyphomicrobiaceae in low-N2O microcosms. Seventeen high-quality metagenome-assembled genomes (MAGs) recovered from the N2O-reducing microcosms harbored nos operons, with all eight MAGs derived from acidic microcosms carrying the Clade II type nosZ and lacking nitrite reductase genes (nirS/K). Five of the eight MAGs recovered from pH 4.5 microcosms represent novel taxa indicating an unexplored N2O-reducing diversity exists in acidic tropical soils. A survey of pH 3.5-5.7 soil metagenome datasets revealed that nosZ genes commonly occur, suggesting broad distribution of N2O reduction potential in acidic soils.

Bioinformatic Approaches for Comparative Analysis of Viruses

The field of viral genomic studies has experienced an unprecedented increase in data volume. New strains of known viruses are constantly being added to the GenBank database and so are completely new species with little or no resemblance to our databases of sequences. In addition to this, metagenomic techniques have the potential to further increase the number and rate of sequenced genomes. Besides, it is important to consider that viruses have a set of unique features that often break down molecular biology dogmas, e.g., the flux of information from RNA to DNA in retroviruses and the use of RNA molecules as genomes. As a result, extracting meaningful information from viral genomes remains a challenge and standard methods for comparing the unknown and our databases of characterized sequences may need adaptations. Thus, several bioinformatic approaches and tools have been created to address the challenge of analyzing viral data. This chapter offers descriptions and protocols of some of the most important bioinformatic techniques for comparative analysis of viruses. The authors also provide comments and discussion on how viruses' unique features can affect standard analyses and how to overcome some of the major sources of problems. Protocols and topics emphasize online tools (which are more accessible to users) and give the real experience of what most bioinformaticians do in day-by-day work with command-line pipelines. The topics discussed include (1) clustering related genomes, (2) whole genome multiple sequence alignments for small RNA viruses, (3) protein alignment for marker genes and species affiliation, (4) variant calling and annotation, and (5) virome analyses and pathogen identification.

New mitogenomes from the genus Cricotopus (Diptera: Chironomidae, Orthocladiinae): Characterization and phylogenetic implications

Cricotopus is a large and diverse genus of non-biting midges composed of several subgenera. Complete mitogenome sequences are available for very few Cricotopus species. The subgenus Pseudocricotopus unites species with unusual morphological structures in adult male and pupal stages, however, molecular methods are needed to verify the placement of this subgenus within Cricotopus. We obtained mitogenomes of C. (Pseudocricotopus) cf. montanus and nine other Cricotopus species for phylogenetic analysis, coupled with two Rheocricotopus species and one Synorthocladius species as outgroups. The structure of the mitogenome was similar among these Cricotopus species, exhibiting A+T bias and retaining ancestral gene order. Mutation rate, estimated as Ka/Ks, varied among genes, and was highest for ATP8 and lowest for COI. The phylogenetic relationships among species of Cricotopus, Rheocricotopus and Synorthocladius was reconstructed using Bayesian inference and maximum likelihood estimation. The phylogenetic trees confirmed placement of subgenus Pseudocricotopus, represented by Cricotopus cf. montanus, within Cricotopus. Our study increases the library of chironomid mitogenomes and provides insight into the properties of their constituent genes.

Halomonas rhizosphaerae sp. nov. and Halomonas kalidii sp. nov., two novel moderate halophilic phenolic acid-degrading species isolated from saline soil

Four vanillic acid-degrading bacterial strains, named LR5S13T, LR5S20, and M4R5S39T and LN1S58, were isolated from Kalidium cuspidatum rhizosphere and bulk soils, respectively. Phylogenetic analysis based on 16S rRNA gene as well as core genome revealed that LR5S13T and LR5S20 clustered closely with each other and with Halomonas ventosae Al12T, and that the two strains shared the highest similarities (both 99.3 %) with H. ventosae Al12T, in contrast, M4R5S39T and LN1S58 clustered together and with Halomonas heilongjiangensis 9-2T, and the two strains shared the highest similarities (99.4 and 99.2 %, respectively) with H. heilongjiangensis 9-2T. The average nucleotides identities based on BLAST (ANIb) and digital DNA-DNA hybridization (dDDH) values of strains LR5S13T to LR5S20, and M4R5S39T to LN1S58, were both higher than the threshold values for delineation of a species. The ANIb and dDDH values of the four strains to their closely relatives were lower than the threshold values. All four strains take phosphatidylethanolamine, phosphatidylglycerol, and diphosphatidylglycerol as the major polar lipids, Summed Feature 8, Summed Feature 3, and C16:0 as the major fatty acids. Based on the phylogenetic and phenotypic results, the four strains should be classified as two novel Halomonas species. Therefore, Halomonas rhizosphaerae sp. nov. (type strain LR5S13T = KCTC 8016T = CGMCC 1.62049T) and Halomonas kalidii (type strain M4R5S39T = KCTC 8015T = CGMCC 1.62047T) are proposed. The geographical distribution analysis based on 16S rRNA gene revealed that the two novel species are widely distributed across the globe, specifically in highly saline habits, especially in Central and Eastern Asia.

Limits of mitochondrial genes in delimiting species within a Carbula species complex (Hemiptera: Pentatomidae)

Molecular data has become a powerful tool for species delimitation, particularly among those that present limited morphological differences; while the mitochondrial genome, with its moderate length, low cost of sequencing and fast lineage sorting, has emerged as a practical data set. Due to the limited morphological differences among the closely related species of Carbula Stål 1865, the species boundaries between Carbula abbreviata (Motschulsky, 1866), Carbula humerigera (Uhler, 1860), and Carbula putoni (Jakovlev, 1876) have remained particularly unclear. In this study, we applied two phylogenetic reconstruction methods to two data sets (mitogenome and COI) to assess the phylogeny of Carbula distributed in Asia, and five species delimitation methods to determine the boundaries between East Asian Carbula species. Our phylogenetic analyses showed Carbula to be paraphyletic; the seven known species distributed within East Asia to form a single monophyletic group, and within this, C. abbreviata, C. humerigera, C. putoni and middle-type to comprise a C. humerigera species complex. Our results show that mitogenome data alone, while effective in the differentiation of more distantly related Carbula species, is not sufficient to accurately delimit the species within this newly described complex.

Relaxed purifying selection pressure drives accelerated and dynamic gene rearrangements in thrips (Insecta: Thysanoptera) mitochondrial genomes

Insect mitochondrial genomes (mitogenome) generally present a typical gene order, which is considered as the ancestral arrangement. All sequenced mitogenomes in the Thysanoptera display high levels of gene rearrangement. Due to limited number of thrips mitogenomes sequenced, how gene rearrangement may be shaped by evolution remain unclear. Here, we analyzed 33 thrips mitogenomes, including 14 newly sequenced. These mitogenomes were diverse in organization, nucleotides substitution and gene arrangements. We found 28 highly rearranged gene orders with the breakpoints of gene rearrangements from 25 to 33. Reconstruction of the ancestors mitochondrial gene arrangements states indicated that Tubulifera have more complex pathways than Terebrantia in the gene order evolution. Molecular calibration estimated that divergence of two suborders occurred in the middle Triassic while the radiation of thrips was associated with the arose and flourish of angiosperm. Our evolutionary hypothesis testing suggests that relaxation of selection pressure enabled the early phase of Thysanoptera evolution, followed by a stronger selective pressure fixed diversification. Our analyses found gene inversion increases the nonsynonymous substitution rates and provide an evolutionary hypothesis driving the diverse gene orders.

Probiotics Bifidobacterium lactis M8 and Lactobacillus rhamnosus M9 prevent high blood pressure via modulating the gut microbiota composition and host metabolic products

IMPORTANCE

Elevated blood pressure affects 40% of the adult population, which accounts for high cardiovascular disease risk and further high mortality yearly. The global understanding of the gut microbiome for hypertension may provide important insights into the prevention. Bifidobacterium lactis M8 and Lactobacillus rhamnosus M9 originated from human breast milk, were able to decrease blood pressure, and modified metabolites in a high fructose-induced elevated blood pressure mouse model. Moreover, we found there was a close relationship between unexplored gut microbes and elevated blood pressure. Also, subsequently, the cross-link was explored among gut microbes, metabolites, and some metabolic pathways in gut microbial environment through introducing novel prediction methodology and bioinformatic analysis. It allowed us to hypothesize that probiotics can prevent elevated blood pressure via gut microbiota and related metabolism.Thus, utilization of dietary strategies (such as probiotics) to maintain the blood pressure level is of crucial importance.

Fe(II) Oxidation Shaped Functional Genes and Bacteria Involved in Denitrification and Dissimilatory Nitrate Reduction to Ammonium from Different Paddy Soils

Microbial nitrate reduction can drive Fe(II) oxidation in anoxic environments, affecting the nitrous oxide emission and ammonium availability. The nitrate-reducing Fe(II) oxidation usually causes severe cell encrustation via chemodenitrification and potentially inhibits bacterial activity due to the blocking effect of secondary minerals. However, it remains unclear how Fe(II) oxidation and subsequent cell encrustation affect the functional genes and bacteria for denitrification and dissimilatory nitrate reduction to ammonium (DNRA). Here, bacteria were enriched from different paddy soils with and without Fe(II) under nitrate-reducing conditions. Fe(II) addition decelerated nitrate reduction and increased NO2- accumulation, due to the rapid Fe(II) oxidation and cell encrustation in the periplasm and on the cell surface. The N2O accumulation was lower in the treatment with Fe(II) and nitrate than that in the treatment with nitrate only, although the proportions of N2O and NH4+ to the reduced NO3- were low (3.25% ∼ 6.51%) at the end of incubation regardless of Fe(II) addition. The dominant bacteria varied from soils under nitrate-reducing conditions, while Fe(II) addition shaped a similar microbial community, including Dechloromonas, Azospira, and Pseudomonas. Fe(II) addition increased the relative abundance of napAB, nirS, norBC, nosZ, and nirBD genes but decreased that of narG and nrfA, suggesting that Fe(II) oxidation favored denitrification in the periplasm and NO2--to-NH4+ reduction in the cytoplasm. Dechloromonas dominated the NO2--to-N2O reduction, while Thauera mediated the periplasmic nitrate reduction and cytoplasmic NO2--to-NH4+ during Fe(II) oxidation. However, Thauera showed much lower abundance than the dominant genera, resulting in slow nitrate reduction and limited NH4+ production. These findings provide new insights into the response of denitrification and DNRA bacteria to Fe(II) oxidation and cell encrustation in anoxic environments.

Tandem gene duplications contributed to high-level azole resistance in a rapidly expanding Candida tropicalis population

Invasive diseases caused by the globally distributed commensal yeast Candida tropicalis are associated with mortality rates of greater than 50%. Notable increases of azole resistance have been observed in this species, particularly within Asia-Pacific regions. Here, we carried out a genetic population study on 1571 global C. tropicalis isolates using multilocus sequence typing (MLST). In addition, whole-genome sequencing (WGS) analysis was conducted on 629 of these strains, comprising 448 clinical invasive strains obtained in this study and 181 genomes sourced from public databases. We found that MLST clade 4 is the predominant azole-resistant clone. WGS analyses demonstrated that dramatically increasing rates of azole resistance are associated with a rapid expansion of cluster AZR, a sublineage of clade 4. Cluster AZR isolates exhibited a distinct high-level azole resistance, which was induced by tandem duplications of the ERG11A395T gene allele. Ty3/gypsy-like retrotransposons were found to be highly enriched in this population. The alarming expansion of C. tropicalis cluster AZR population underscores the urgent need for strategies against growing threats of antifungal resistance.

Different microbial assemblage colonized on microplastics and clay particles in aerobic sludge treatment

In this study, a combination of property analysis and high-throughput sequencing was used to investigate the microbial colonization ability and their community structures and functions in polypropylene microplastics (PPMPs), polystyrene microplastics (PSMPs) and montmorillonite (MMT), respectively as the representatives of artificial and natural substrates in aerobic sludge treatment. After 45 d of incubation, the surface properties of substrates were altered with the increased oxygen functional groups and surface roughness, indicating microbial settlement. Moreover, MPs had different microbial structures from that of MMT, and PSMPs exhibited higher microbial diversity and abundance than PPMPs and MMT. Also, these substrates changed the inherent ecological niche in sludge. Especially, the abundance of some pathogens (e.g., Pseudomonas, Klebsiella and Flavobacterium) was increased in MPs, and the disease risk of Kyoto Encyclopedia of Genes and Genomes metabolic pathway (e.g., Infectious diseases: Bacterial, Infectious diseases: Parasitic and Immune diseases) was higher. Also, the presence of MPs inhibited the decomposition of organic matter including soluble chemical oxygen demand and protein compared to natural substrates. The findings revealed the crucial vector role of MPs for microbes and the effect on aerobic sludge treatment, highlighting the necessity of MP removal in sludge.

A novel botybirnavirus with a unique satellite dsRNA causes latent infection in Didymella theifolia isolated from tea plants

IMPORTANCE

A novel botybirnavirus, infecting the tea plant pathogen Didymella theifolia and tentatively named Didymella theifolia botybirnavirus 1 (DtBRV1), together with an additional double-stranded RNA (dsRNA), was characterized. DtBRV1 comprises two dsRNAs (1 and 2) encapsidated in isometric virions, while dsRNA3 is a satellite. The satellite represents a unique specimen since it contains a duplicated region and has high similarity to the two botybirnavirus dsRNAs, supporting the notion that it most likely originated from a deficient genomic component. The biological characteristics of DtBRV1 were further determined. With their unique molecular traits, DtBRV1 and its related dsRNA expand our understanding of virus diversity, taxonomy, and evolution.

Molecular characterization and pathogenicity evaluation of enterovirus G isolated from diarrheic piglets

IMPORTANCE

Enterovirus G is a species of positive-sense single-stranded RNA viruses associated with several mammalian diseases. The porcine enterovirus strains isolated here were chimeric viruses with the PLCP gene of porcine torovirus, which grouped together with global EV-G1 strains. The isolated EV-G strain could infect various cell types from different species, suggesting its potential cross-species infection risk. Animal experiment showed the pathogenic ability of the isolated EV-G to piglets. Additionally, the EV-Gs were widely distributed in the swine herds. Our findings suggest that EV-G may have evolved a novel mechanism for broad tropism, which has important implications for disease control and prevention.

A metagenomics method for the quantitative detection of bacterial pathogens causing hospital-associated and ventilator-associated pneumonia

IMPORTANCE

The management of ventilator-associated pneumonia and hospital-acquired pneumonia requires rapid and accurate quantitative detection of the infecting pathogen. To this end, we propose a metagenomic sequencing assay that includes the use of an internal sample processing control for the quantitative detection of 20 relevant bacterial species from bronchoalveolar lavage samples.

Comparative and phylogenetic analyses of mitochondrial genomes in Elateridae (Coleoptera: Elateroidea)

The click-beetles (Elateridae) are a species-rich beetle family that is easily recognizable. They are distributed in all zoogeographical regions with over 11,000 species. Comparative studies of the structural characteristics of mitochondrial genomes (mitogenomes), as well as phylogenetic relationships of click-beetles, can improve our understanding of mitogenomic evolution. In this study, we determined four mitogenomes from Elateridae by next-generation sequencing. The four mitogenomes were 16,005 to 16,930 bp in length with 37 typical genes and a control region (A + T-rich region). Combined with previously reported elaterid mitogenomes, all PCGs initiate with either the standard start codon of ATN or TTG. According to the nonsynonymous/synonymous mutation ratio (Ka/Ks) of all PCGs, the highest and the lowest evolutionary rates were found for atp8 and cox1, respectively. Among the control regions of the four mitogenomes, several different patterns and numbers of tandem repeats were identified, which was the primary cause of the length variation in control regions. Phylogenetic analyses were conducted based on 13 protein-coding genes and two ribosomal RNA genes from 33 species of Elateridae and two outgroups. The Bayesian inference and maximum likelihood trees had an identical topological structure. The monophyly of Cardiophorinae, Agrypninae and Elaterinae was recovered with high support in all topologies, and the Tetralobinae was placed as the earliest branch in the Elateridae. Expanding the availability of mitogenomic and genomic data from a broader range of click-beetles could provide more clarity on the disputed relationships among subfamilies within Elateridae.

Metagenomics reveals the self-recovery and risk of antibiotic resistomes during carcass decomposition of wild mammals

Animal carcass decomposition may bring serious harm to the environment, including pathogenic viruses, toxic gases and metabolites, and antibiotic resistance genes (ARGs). However, how wild mammal corpses decomposition influence and change ARGs in the environment has less explored. Through metagenomics, 16S rRNA gene sequencing, and physicochemical analysis, this study explored the succession patterns, influencing factors, and assembly process of ARGs and mobile genetic elements (MGEs) in gravesoil during long-term corpse decomposition of wild mammals. Our results indicate that the ARG and MGE communities related to wildlife corpses exhibited a pattern of differentiation first and then convergence. Different from the farmed animals, the decomposition of wild animals first reduced the diversity of ARGs and MGEs, and then recovered to a level similar to that of the control group (untreated soil). ARGs and MGEs of the gravesoil are mainly affected by deterministic processes in different stages. MGEs and bacterial community are the two most important factors affecting ARGs in gravesoil. It is worth noting that the decomposition of wild animal carcasses enriched different high-risk ARGs at different stages (bacA, mecA and floR), which have co-occurrence patterns with opportunistic pathogens (Comamonas and Acinetobacter), thereby posing a great threat to public health. These results are of great significance for wildlife corpse management and environmental and ecological safety.

Why are so many enteric pathogen infections asymptomatic? Pathogen and gut microbiome characteristics associated with diarrhea symptoms and carriage of diarrheagenic E. coli in northern Ecuador

A high proportion of enteric infections, including those caused by diarrheagenic Escherichia coli (DEC), are asymptomatic for diarrhea. The factors responsible for the development of diarrhea symptoms, or lack thereof, remain unclear. Here, we used DEC isolate genome and whole stool microbiome data from a case-control study of diarrhea in Ecuador to examine factors associated with diarrhea symptoms accompanying DEC carriage. We investigated i) pathogen abundance, ii) gut microbiome characteristics, and iii) strain-level pathogen characteristics from DEC infections with diarrhea symptoms (symptomatic infections) and without diarrhea symptoms (asymptomatic infections). We also included data from individuals with and without diarrhea who were not infected with DEC (uninfected cases and controls). i) E. coli relative abundance in the gut microbiome was highly variable, but higher on-average in individuals with symptomatic compared to asymptomatic DEC infections. Similarly, the number and relative abundances of virulence genes in the gut were higher in symptomatic than asymptomatic DEC infections. ii) Measures of microbiome diversity were similar regardless of diarrhea symptoms or DEC carriage. Proteobacterial families that have been described as pathobionts were enriched in symptomatic infections and uninfected cases, whereas potentially beneficial taxa, including the Bacteroidaceae and Bifidobacteriaceae, were more abundant in individuals without diarrhea. An analysis of high-level gene functions recovered in metagenomes revealed that genes that were differentially abundant by diarrhea and DEC infection status were more abundant in symptomatic than asymptomatic DEC infections. iii) DEC isolates from symptomatic versus asymptomatic individuals showed no significant differences in virulence or accessory gene content, and there was no phylogenetic signal associated with diarrhea symptoms. Together, these data suggest signals that distinguish symptomatic from asymptomatic DEC infections. In particular, the abundance of E. coli, the virulence gene content of the gut microbiome, and the taxa present in the gut microbiome have an apparent role.

GlcNac produced by the gut microbiome enhances host influenza resistance by modulating NK cells

Microbiota are known to modulate the host response to influenza infection, but the mechanisms remain largely unknown. Gut metabolites are the key mediators through which gut microbes play anti-influenza effect. Transferring fecal metabolites from mice with high influenza resistance into antibiotic-treated recipient mice conferred resistance to influenza infections. By comparing the metabolites of different individuals with high or low influenza resistance, we identified and validated N-acetyl-D-glucosamine (GlcNAc) and adenosine showed strong positive correlations with influenza resistance and exerted anti-influenza effects in vivo or in vitro, respectively. Especially, GlcNAc mediated the anti-influenza effect by increasing the proportion and activity of NK cells. Several gut microbes, including Clostridium sp., Phocaeicola sartorii, and Akkermansia muciniphila, were positively correlated with influenza resistance, and can upregulate the level of GlcNAc in the mouse gut by exogenous supplementation. Subsequent studies confirmed that administering a combination of the three bacteria to mice via gavage resulted in similar modulation of NK cell responses as observed with GlcNAc. This study demonstrates that gut microbe-produced GlcNAc protects the host against influenza by regulating NK cells, facilitating the elucidation of the action mechanism of gut microbes mediating host influenza resistance.

Invasive non-typhoidal Salmonella (iNTS) aminoglycoside-resistant ST313 isolates feature unique pathogenic mechanisms to reach the bloodstream

Invasive non-typhoidal Salmonella (iNTS) from the clonal type ST313 (S. Typhimurium ST313) is the primary cause of invasive salmonellosis in Africa. Recently, in Brazil, iNTS ST313 strains have been isolated from different sources, but there is a lack of understanding of the mechanisms behind how these gut bacteria can break the gut barrier and reach the patient's bloodstream. Here, we compare 13 strains of S. Typhimurium ST313, previously unreported isolates, from human blood cultures, investigating aspects of virulence and mechanisms of resistance. Initially, RNAseq analyses between ST13-blood isolate and SL1344 (ST19) prototype revealed 15 upregulated genes directly related to cellular invasion and replication, such as sopD2, sifB, and pipB. Limited information is available about S. Typhimurium ST313 pathogenesis and epidemiology, especially related to the global distribution of strains. Herein, the correlation of strains isolated from different sources in Brazil was employed to compare clinical and non-clinical isolates, a total of 22 genomes were studied by single nucleotide polymorphism (SNPs). The epidemiological analysis of 22 genomes of S. Typhimurium ST313 strains grouped them into three distinct clusters (A, B, and C) by SNP analysis, where cluster A comprised five, group B six, and group C 11. The 13 clinical blood isolates were all resistant to streptomycin, 92.3% of strains were resistant to ampicillin and 15.39% were resistant to kanamycin. The resistance genes acrA, acrB, mdtK, emrB, emrR, mdsA, and mdsB related to the production of efflux pumps were detected in all (100%) strains studied, similar to pathogenic traits investigated. In conclusion, we evidenced that S. Typhimurium ST313 strains isolated in Brazil have unique epidemiology. The elevated frequencies of virulence genes such as sseJ, sopD2, and pipB are a major concern in these Brazilian isolates, showing a higher pathogenic potential.

Hybrid-hybrid correction of errors in long reads with HERO

Although generally superior, hybrid approaches for correcting errors in third-generation sequencing (TGS) reads, using next-generation sequencing (NGS) reads, mistake haplotype-specific variants for errors in polyploid and mixed samples. We suggest HERO, as the first "hybrid-hybrid" approach, to make use of both de Bruijn graphs and overlap graphs for optimal catering to the particular strengths of NGS and TGS reads. Extensive benchmarking experiments demonstrate that HERO improves indel and mismatch error rates by on average 65% (27[Formula: see text]95%) and 20% (4[Formula: see text]61%). Using HERO prior to genome assembly significantly improves the assemblies in the majority of the relevant categories.

Rhizocompartmental microbiomes of arrow bamboo (Fargesia nitida) and their relation to soil properties in Subalpine Coniferous Forests

Arrow bamboo (Fargesia nitida) is a pioneer plant in secondary forest succession in the Sichuan Province mountains. To comprehensively investigate the microbial communities and their functional variations in different rhizocompartments (root endosphere, rhizosphere, and root zone) of arrow bamboo (Fargesia nitida), a high-throughput metagenomic study was conducted in the present study. The results showed that the abundances of the dominant bacterial phyla Proteobacteria and Actinobacteria in the bamboo root endosphere were significantly lower than those in the rhizosphere and root zones. In contrast, the dominant fungal phyla, Ascomycota and Basidiomycota, showed the opposite tendency. Lower microbial diversity, different taxonomic composition and functional profiles, and a greater abundance of genes involved in nitrogen fixation (nifB), cellulose degradation (beta-glucosidase), and cellobiose transport (cellulose 1, 4-beta-cellobiosidase) were found in the bamboo root endosphere than in the other rhizocompartments. Greater soil total carbon, total nitrogen, NH4+-N, microbial biomass carbon, and greater activities of invertase and urease were found in the bamboo root zone than in the adjacent soil (spruce root zone). In contrast, the soil microbial community and functional profiles were similar. At the phylum level, invertase was significantly related to 31 microbial taxa, and the effect of NH4+-N on the microbial community composition was greater than that of NO3--N. The soil physicochemical properties and enzyme activities were significantly correlated with microbial function. These results indicate that the root endosphere microbiomes of arrow bamboo were strongly selected by the host plant, which caused changes in the soil nutrient properties in the subalpine coniferous forest.

Infant microbiome cultivation and metagenomic analysis reveal Bifidobacterium 2'-fucosyllactose utilization can be facilitated by coexisting species

The early-life gut microbiome development has long-term health impacts and can be influenced by factors such as infant diet. Human milk oligosaccharides (HMOs), an essential component of breast milk that can only be metabolized by some beneficial gut microorganisms, ensure proper gut microbiome establishment and infant development. However, how HMOs are metabolized by gut microbiomes is not fully elucidated. Isolate studies have revealed the genetic basis for HMO metabolism, but they exclude the possibility of HMO assimilation via synergistic interactions involving multiple organisms. Here, we investigate microbiome responses to 2'-fucosyllactose (2'FL), a prevalent HMO and a common infant formula additive, by establishing individualized microbiomes using fecal samples from three infants as the inocula. Bifidobacterium breve, a prominent member of infant microbiomes, typically cannot metabolize 2'FL. Using metagenomic data, we predict that extracellular fucosidases encoded by co-existing members such as Ruminococcus gnavus initiate 2'FL breakdown, thus critical for B. breve's growth. Using both targeted co-cultures and by supplementation of R. gnavus into one microbiome, we show that R. gnavus can promote extensive growth of B. breve through the release of lactose from 2'FL. Overall, microbiome cultivation combined with genome-resolved metagenomics demonstrates that HMO utilization can vary with an individual's microbiome.

Multi-genome comparisons reveal gain-and-loss evolution of the anti-Mullerian hormone receptor type 2 gene, an old master sex determining gene, in Percidae

The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems. We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplications (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome-18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variants (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in testis than ovary. Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

Metagenomics for the microbiological diagnosis of hospital-acquired pneumonia and ventilator-associated pneumonia (HAP/VAP) in intensive care unit (ICU): a proof-of-concept study

BACKGROUND

Hospital-acquired and ventilator-associated-pneumonia (HAP/VAP) are one of the most prevalent health-care associated infections in the intensive care unit (ICU). Culture-independent methods were therefore developed to provide faster route to diagnosis and treatment. Among these, metagenomic next-generation sequencing (mNGS) has shown considerable promise.

METHODS

This proof-of-concept study describes the technical feasibility and evaluates the clinical validity of the mNGS for the detection and characterization of the etiologic agents causing hospital-acquired and ventilator-associated pneumonia. We performed a prospective study of all patients with HAP/VAP hospitalized in our intensive care unit for whom a bronchoalveolar lavage (BAL) was performed between July 2017 and November 2018. We compared BAL fluid culture and mNGS results of these patients.

RESULTS

A total of 32 BAL fluids were fully analyzed. Of these, 22 (69%) were positive by culture and all pathogens identified were also reported by mNGS. Among the culture-positive BAL samples, additional bacterial species were revealed by mNGS for 12 patients, raising the issue of their pathogenic role (colonization versus coinfection). Among BALF with culture-negative test, 5 were positive in mNGS test.

CONCLUSIONS

This study revealed concordant results for pneumonia panel pathogens between mNGS and culture-positive tests and identified additional pathogens potentially implicated in pneumonia without etiologic diagnosis by culture. mNGS has emerged as a promising methodology for infectious disease diagnoses to support conventional methods. Prospective studies with real-time mNGS are warranted to examine the impact on antimicrobial decision-making and clinical outcome.

Microbial associates of an endemic Mediterranean seagrass enhance the access of the host and the surrounding seawater to inorganic nitrogen under ocean acidification

Seagrasses are important primary producers in oceans worldwide. They live in shallow coastal waters that are experiencing carbon dioxide enrichment and ocean acidification. Posidonia oceanica, an endemic seagrass species that dominates the Mediterranean Sea, achieves high abundances in seawater with relatively low concentrations of dissolved inorganic nitrogen. Here we tested whether microbial metabolisms associated with P. oceanica and surrounding seawater enhance seagrass access to nitrogen. Using stable isotope enrichments of intact seagrass with amino acids, we showed that ammonification by free-living and seagrass-associated microbes produce ammonium that is likely used by seagrass and surrounding particulate organic matter. Metagenomic analysis of the epiphytic biofilm on the blades and rhizomes support the ubiquity of microbial ammonification genes in this system. Further, we leveraged the presence of natural carbon dioxide vents and show that the presence of P. oceanica enhanced the uptake of nitrogen by water column particulate organic matter, increasing carbon fixation by a factor of 8.6-17.4 with the greatest effect at CO2 vent sites. However, microbial ammonification was reduced at lower pH, suggesting that future ocean climate change will compromise this microbial process. Thus, the seagrass holobiont enhances water column productivity, even in the context of ocean acidification.

Gastrointestinal microbiome, resistance genes, and risk assessment of heavy metals in wild giant pandas

The gastrointestinal microbiome (GM) of giant panda (GP) plays an important role in food utilization and health and is also an essential reservoir of resistance genes. Currently, little knowledge is available on the GM, acid resistance genes (AcRGs), antibiotic resistance genes (ARGs), metal resistance genes (MRGs), and mobile genetic elements (MGEs) in wild GPs. We sampled the gastrointestinal tract of a dead GP and explored the composition and function of GM and resistance genes through cryo-scanning electron microscopy, metagenomic sequencing, and genome-resolved metagenomics. The concentration of metals in the gastrointestinal lumen, feces, bamboo, and soil was measured by inductively coupled plasma mass spectrometry. Results showed that the composition of the microbiota varied in different gastrointestinal regions. Fecal microbiota was highly associated with small intestinal and colonic microbes. The lignocellulosic cross-linked structure of bamboo was destroyed in the stomach initially and destroying degree increased from stomach to anus. Reconstruction of metagenome-assembled-genomes confirmed that core GM, e.g., Streptococcus, Clostridium, Lactococcus, Leuconostoc, and Enterococcus, carried genes encoding the lignocellulose degradation enzyme. There were no significant differences of resistance genes between gastrointestinal and fecal samples, except MGEs. Multidrug and multi-metal resistance genes were predominant in all samples, while the transposase gene tnpA was the major type of MGE. Significant correlations were observed among the abundance of GM, resistance genes, and MGEs. Gastrointestinal and fecal mercury and chromium were the main metals influencing GM and resistance genes. The content of gastrointestinal and fecal metals was significantly associated with the presence of the same metals in bamboo, which could pose a threat to the health of wild GPs. This study characterized the gastrointestinal microbiome of wild GPs, providing new evidence for the role of the gastrointestinal microbiome in degrading lignocellulose from bamboo and highlighting the urgent need to monitor metal levels in soil and bamboo.

Integrating metagenomics and metabolomics to study the response of microbiota in black soil degradation

As the largest commercial food production base and ecological security barrier, land degradation in black soil areas seriously threatens the global food supply and natural ecosystems. Therefore, determining the response of soil microbiota is crucial to restoring degraded soils. This study combined metagenomics and metabolomics to investigate the effect of different degrees of soil degradation on microbial community composition and metabolic function in black soils. It was found that alpha diversity in degraded soils (Shannon: 22.3) was higher than in nondegraded soil (ND) (Shannon: 21.8), and the degree of degradation significantly altered the structure and composition of soil microbial communities. The results of LEfSe analysis obtained 9 (ND), 7 (lightly degraded, LD), 10 (moderately degraded, MD), and 1 (severely degraded, SD) biomarkers in four samples. Bradyrhizobium, Sphingomonas, and Ramlibacter were significantly affected by soil degradation and can be considered biomarkers of ND, MD, and SD, respectively. Soil nutrient and enzyme activities decreased significantly with increasing black soil degradation, soil organic matter (SOM) content decreased from 11.12 % to 1.97 %, and Sucrase decreased from 23.53 to 6.59 mg/g/d. In addition, C was the critical driver affecting microbial community structure, contributing 61.2 % to differences in microbial community distribution, and microbial altering relative abundance which participle in the carbon cycle to respond to soil degradation. Metabolomic analyses indicated that soil degradation significantly modified the soil metabolite spectrum, and the metabolic functions of most microorganisms responding to soil degradation were adversely affected. The combined multi-omics analysis further indicated that biomarkers dominate in accumulating metabolites. These findings confirmed that due to their role in the composition and functioning of these degraded soils, these biomarkers could be employed in strategies for managing and restoring degraded black soils.

New Mitogenomes of the Green Lacewing Tribe Ankylopterygini (Neuroptera: Chrysopidae: Chrysopinae) and Phylogenetic Implications of Chrysopidae

Chrysopidae (green lacewings) are a cosmopolitan and species-rich family of Neuroptera, with remarkable significance of biological control against various agricultural and forestry pests. However, the phylogenetic position of Chrysopidae in Neuroptera and the internal relationships within the family remain equivocal among previous studies based on different types of data and sampling. Here we sequenced the mitochondrial genomes (mitogenomes) of two species of the genus Ankylopteryx in the chrysopine tribe Ankylopterygini for the first time. The characteristics of these mitogenomes were analyzed in comparison with other green lacewing mitogenomes published to date. In the phylogeny herein reconstructed based on mitogenomes, Chrysopinae were recovered as the sister group to Apochrysinae + Nothochrysinae. Within the subfamily of Chrysopinae, Nothancylini were recovered as the sister group to (Leucochrysini + Belonopterygini) + (Ankylopterygini + Chrysopini). The divergence time estimation suggested an Early Cretaceous initial divergence within the extant Chrysopidae. Within Chrysopinae, the four tribes except Nothancylini diverged around mid-Cretaceous.

Phylogenetics and biomineralization of a novel magnetotactic Gammaproteobacterium from a freshwater lake in Beijing, China

Magnetotactic bacteria (MTB) have the remarkable capability of producing intracellularly membrane-enveloped magnetic nanocrystals (i.e. magnetosomes) and swimming along geomagnetic field lines. Despite more than 50 years of research, bacterial diversity and magnetosome biomineralization within MTB are relatively less known in the Gammaproteobacteria class than other groups. This is incompatible with the status of Gammaproteobacteria as the most diverse class of gram-negative bacteria with a number of ecologically important bacteria. Here, we identify a novel MTB strain YYHR-1 affiliated with the Gammaproteobacteria class of the Pseudomonadota phylum from a freshwater lake. In YYHR-1, most magnetosome crystals are organized into a long chain aligned along the cell long axis; unusually, a few small superparamagnetic crystals are located at the side of the chain, off the main chain axis. Micromagnetic simulations indicate that magnetostatic interactions among adjacent crystals within a chain reduce the Gibbs energy to enhance chain stability. Genomic analysis suggests that duplication of magnetosome gene clusters may result in off-chain magnetosomes formation. By integrating available genomic data from Gammaproteobacteria, the phylogenetic position of MTB in this class is reassigned here. Our new findings expand knowledge about MTB diversity and magnetosome biomineralization, and deepen understanding of the phylogenetics of the Gammaproteobacteria.

A novel long-tailed myovirus represents a new T4-like cyanophage cluster

Cyanophages affect the abundance, diversity, metabolism, and evolution of picocyanobacteria in marine ecosystems. Here we report an estuarine Synechococcus phage, S-CREM2, which represents a novel viral genus and leads to the establishment of a new T4-like cyanophage clade named cluster C. S-CREM2 possesses the longest tail (~418 nm) among isolated cyanomyoviruses and encodes six tail-related proteins that are exclusively homologous to those predicted in the cluster C cyanophages. Furthermore, S-CREM2 may carry three regulatory proteins in the virion, which may play a crucial role in optimizing the host intracellular environment for viral replication at the initial stage of infection. The cluster C cyanophages lack auxiliary metabolic genes (AMGs) that are commonly found in cyanophages of the T4-like clusters A and B and encode unique AMGs like an S-type phycobilin lyase gene. A variation in the composition of tRNA and cis-regulatory RNA genes was observed between the marine and freshwater phage strains in cluster C, reflecting their different modes of coping with hosts and habitats. The cluster C cyanophages are widespread in estuarine and coastal regions and exhibit equivalent or even higher relative abundance compared to those of clusters A and B cyanophages in certain estuarine regions. The isolation of cyanophage S-CREM2 provides new insights into the phage-host interactions mediated by both newly discovered AMGs and virion-associated proteins and emphasizes the ecological significance of cluster C cyanophages in estuarine environments.

Inulin-enriched Megamonas funiformis ameliorates metabolic dysfunction-associated fatty liver disease by producing propionic acid

Accumulated evidence supports the beneficial role of inulin in alleviating metabolic dysfunction-associated fatty liver disease (MAFLD) by modulating gut microbiota. However, the underlying mechanisms are not fully understood. Here we used high-fat diet (HFD)-induced laying hen model of MAFLD to investigate the effect of inulin on ameliorating MAFLD and found that the inulin-enriched Megamonas genus was inversely correlated with hepatic steatosis-related parameters. Oral administration of a newly isolated commensal bacterium by culturomics, M. funiformis CML154, to HFD-fed hens and mice ameliorated MAFLD, changed liver gene expression profiles, and increased intestinal propionate concentration. Further evidence demonstrated that the anti-MAFLD effect of M. funiformis CML154 is attributed to propionate-mediated activation of the APN-AMPK-PPARα signaling pathway, thereby inhibiting fatty acid de novo synthesis and promoting β-oxidation. These findings establish the causal relationships among inulin, M. funiformis, and MAFLD, and suggest that M. funiformis CML154 is a probiotic candidate for preventative or therapeutic intervention of MAFLD.

Novel broad-spectrum bacteriophages against Xanthomonas oryzae and their biocontrol potential in rice bacterial diseases

Bacterial leaf blight (BLB) and bacterial leaf streak (BLS)-caused by Xanthomonas oryzae pv. oryzae (Xoo) and Xanthomonas oryzae pv. oryzicola (Xoc), respectively-are two major bacterial diseases that threaten the safe production of rice, one of the most important food crops. Bacteriophages are considered potential biocontrol agents against rice bacterial pathogens, due to their host specificity and environmental safety. It is common for BLB and BLS to occur together in fields, which highlights the need for broad-spectrum phages capable of infecting both Xoo and Xoc. In this study, two lytic broad-spectrum phages (pXoo2106 and pXoo2107) that can infect various strains of Xoo and Xoc were assessed. Both phages belong to the class Caudoviricetes and one of them to the family Autographiviridae, while the other belongs to an unclassified family. Two phages alone or combined in a phage cocktail could effectively inhibit Xoo and Xoc growth in vitro. In an in vivo biocontrol experiment, the phage cocktail reduced the total CFU and significantly eased the symptoms caused by Xoo or Xoc. Our results suggest that pXoo2106 and pXoo2107 have a broad-spectrum host range targeting different X. oryzae strains, and have strong biocontrol potential in field applications against both BLB and BLS.

Metagenomic Analysis of Bacterial Community and Isolation of Representative Strains from Vranjska Banja Hot Spring, Serbia

The hot spring Vranjska Banja is the hottest spring on the Balkan Peninsula with a water temperature of 63-95 °C and a pH value of 7.1, in situ. According to the physicochemical analysis, Vranjska Banja hot spring belongs to the bicarbonated and sulfated hyperthermal waters. The structures of microbial community of this geothermal spring are still largely unexplored. In order to determine and monitor the diversity of microbiota of the Vranjska Banja hot spring, a comprehensive culture-independent metagenomic analysis was conducted in parallel with a culture-dependent approach for the first time. Microbial profiling using amplicon sequencing analysis revealed the presence of phylogenetically novel taxa, ranging from species to phyla. Cultivation-based methods resulted in the isolation of 17 strains belonging to the genera Anoxybacillus, Bacillus, Geobacillus, and Hydrogenophillus. Whole-genome sequencing of five representative strains was then performed. The genomic characterization and OrthoANI analysis revealed that the Vranjska Banja hot spring harbors phylogenetically novel species of the genus Anoxybacillus, proving its uniqueness. Moreover, these isolates contain stress response genes that enable them to survive in the harsh conditions of the hot springs. The results of the in silico analysis show that most of the sequenced strains have the potential to produce thermostable enzymes (proteases, lipases, amylases, phytase, chitinase, and glucanase) and various antimicrobial molecules that can be of great importance for industrial, agricultural, and biotechnological applications. Finally, this study provides a basis for further research and understanding of the metabolic potential of these microorganisms.