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

The first mitogenomic phylogenetic framework of Dorcus sensu lato (Coleoptera: Lucanidae), with an emphasis on generic taxonomy in Eastern Asia

BACKGROUND

Dorcus stag beetles in broad sense are one of the most diverse group in Lucanidae and important saproxylic insects playing a crucial role in nutrient recycling and forest biomonitoring. However, the dazzling morphological differentiations have caused numerous systematic confusion within the big genus, especially the puzzlingly generic taxonomy. So far, there is lack of molecular phylogenetic study to address the chaotic situation. In this study, we undertook mitochondrial genome sequencing of 42 representative species including 18 newly-sequenced ones from Eastern Asia and reconstructed the phylogenetic framework of stag beetles in Dorcus sensu lato for the first time.

RESULTS

The mitogenome datasets of Dorcus species have indicated the variable mitogenomic lengths ranged from 15,785 to 19,813 bp. Each mitogenome contained 13 PCGs, 2 rRNAs, 22 tRNAs, and a control region, and all PCGs were under strong purifying selection (Ka/Ks < 1). Notably, we have identified the presence of a substantial intergenic spacer (IGS) between the trnAser (UCN) and NAD1 genes, with varying lengths ranging from 129 bp (in D. hansi) to 158 bp (in D. tityus). The mitogenomic phylogenetic analysis of 42 species showed that Eastern Asia Dorcus was monophyletic, and divided into eight clades with significant genetic distance. Four of them, Clade VIII, VII, VI and I are clustered by the representative species of Serrognathus Motschulsky, Kirchnerius Schenk, Falcicornis Séguy and Dorcus s.s. respectively, which supported their fully generic positions as the previous morphological study presented. The topology also showed the remaining clades were distinctly separated from the species of Dorcus sensu lato, which implied that each of them might demonstrate independent generic status. The Linnaeus nomenclatures were suggested as Eurydorcus Didier stat. res., Eurytrachellelus Didier stat. res., Hemisodorcus Thomson stat. res. and Velutinodorcus Maes stat. res. For Clade V, IV, III and II respectively.

CONCLUSION

This study recognized the monophyly of Dorcus stag beetles and provided a framework for the molecular phylogeny of this group for the first time. The newly generated mitogenomic data serves as a valuable resource for future investigations on lucanid beetles. The generic relationship would facilitate the systematics of Dorcus stag beetles and thus be useful for exploring their evolutionary, ecological, and conservation aspects.

Whole genome sequence analysis of shallot virus X from India reveals it to be a natural recombinant with positive selection pressure

BACKGROUND

Shallots are infected by various viruses like Onion yellow dwarf virus (OYDV), Leek yellow stripe virus (LYSV), Shallot latent virus (SLV) and Shallot virus X (ShVX). In India, they have been found to be persistently infected by ShVX. ShVX also infects onion and garlic in combination with other carlaviruses and potyviruses. ShVX is a member of genus Allexivirus of family Alphaflexiviridae. ShVX has a monopartite genome, which is represented by positive sense single-stranded RNA. Globally, only six complete and 3 nearly complete genome sequences of ShV X are reported to date. This number is insufficient to measure a taxon's true molecular diversity. Moreover, the complete genome sequence of ShVX from Asia has not been reported as yet. Therefore, this study was undertaken to generate a complete genome sequence of ShVX from India.

RESULTS

Shallot virus X (ShVX) is one of the significant threats to Allium crop production. In this study, we report the first complete genome sequence of the ShVX from India through Next-generation sequencing (NGS). The complete genome of the ShVX (Accession No. OK104171), from this study comprised 8911 nucleotides. In-silico analysis of the sequence revealed variability between this isolate and isolates from other countries. The dissimilarities are spread all over the genome specifically some non-coding intergenic regions. Statistical analysis of individual genes for site-specific selection indicates a positive selection in NABP region. The presence of a recombination event was detected in coat protein region. The sequence similarity percentage and phylogenetic analysis indicate ShVX Indian isolate is a distinctly different isolate. Recombination and site-specific selection may have a function in the evolution of this isolate. This is the first detailed study of the ShVX complete genome sequence from Southeast Asia.

CONCLUSION

This study presents the first report of the entire genome sequence of an Indian isolate of ShVX along with an in-depth exploration of its evolutionary traits. The findings highlight the Indian variant as a naturally occurring recombinant, emphasizing the substantial role of recombination in the evolution of this viral species. This insight into the molecular diversity of strains within a specific geographical region holds immense significance for comprehending and forecasting potential epidemics. Consequently, the insights garnered from this research hold practical value for shaping ShVX management strategies and providing a foundation for forthcoming studies delving into its evolutionary trajectory.

Metagenomics of two aquifers with thermal anomalies in Mallorca Island, and proposal of new uncultivated taxa named following the rules of SeqCode

Groundwater offers an intriguing blend of distinctive physical and chemical conditions, constituting a challenge for microbial life. In Mallorca, the largest island of Balearic archipelago, harbours a variety of thermal anomalies (i.e., geothermal manifestation where surface aquifers exhibiting temperatures exceeding the regional average). The metagenomes of two aquifers in the centre and southern of the island showed Pseudomonadota to be the most represented phylum when using extracted 16S rRNA gene sequences. However, the microbial structures within and between aquifers were remarkably diverse but similar in their metabolic profiles as revealed by the metagenome-assembled genomes (MAGs) pointing to a prevalence of aerobic chemolithoautotrophic and heterotrophic metabolisms, especially in the Llucmajor aquifer. Also, some evidences of anaerobic lifestyles were detected, which would indicate that these environments either could suffer episodes of oxygen depletion or the anaerobes had been transported from deeper waters. We believe that the local environmental factors (temperature, external inputs or chemistry) seem to be more relevant than the connection and, eventually, transport of microbial cells within the aquifer in determining the highly divergent structures. Notably, most of the reconstructed genomes belonged to undescribed bacterial lineages and from them two high-quality MAGs could be classified as novel taxa named following the rules of the Code for Nomenclature of Prokaryotes Described from Sequence Data (SeqCode). Accordingly, we propose the new species and genus Costitxia debesea gen. nov., sp. nov., affiliated with the novel family Costitxiaceae fam. nov., order Costitxiales ord. nov. and class Costitxiia class. nov.; and the new new species and genus Lloretia debesea gen. nov. sp. nov. affiliated with the novel family Lloretiaceae fam. nov.

A metagenomic approach to unveil the association between fecal gut microbiota and short-chain fatty acids in diarrhea caused by diarrheagenic Escherichia coli in children

Diarrheagenic Escherichia coli (DEC) is the main cause of diarrhea in children under five years old. The virulence of DEC is tightly regulated by environmental signals influenced by the gut microbiota and its metabolites. Short-chain fatty acids (SCFAs) are the main metabolic product of anaerobic fermentation in the gut, but their role in DEC diarrhea has not yet been established. In this study, we determine the levels of acetate, propionate, and butyrate in stool samples from children with diarrhea caused by DEC, and we identify bacteria from the fecal gut microbiota associated with the production of SCFAs. The microbiota and SCFAs levels in stool samples obtained from 40 children with diarrhea and 43 healthy children were determined by 16S rRNA gene sequencing and HPLC, respectively. Additionally, shotgun metagenomics was used to identify metagenome-assembled genomes (MAGs) in a subgroup of samples. The results showed significantly higher levels of all SCFAs tested in diarrheal samples than in healthy controls. The abundance of Streptococcus sp., Limosilactobacillus, Blautia, Escherichia, Bacteroides, Megamonas, and Roseburia was higher in the DEC group than in healthy individuals. Functional analysis of bacteria and their main metabolic pathways made it possible to identify species MAGs that could be responsible for the detected SCFAs levels in DEC-positive diarrhea. In conclusion, based on our results and published data, we suggest that SCFAs may be important in the crosstalk between the microbiota and DEC pathogens in the gut.

Sulfur oxidation and reduction are coupled to nitrogen fixation in the roots of the salt marsh foundation plant Spartina alterniflora

Heterotrophic activity, primarily driven by sulfate-reducing prokaryotes, has traditionally been linked to nitrogen fixation in the root zone of coastal marine plants, leaving the role of chemolithoautotrophy in this process unexplored. Here, we show that sulfur oxidation coupled to nitrogen fixation is a previously overlooked process providing nitrogen to coastal marine macrophytes. In this study, we recovered 239 metagenome-assembled genomes from a salt marsh dominated by the foundation plant Spartina alterniflora, including diazotrophic sulfate-reducing and sulfur-oxidizing bacteria. Abundant sulfur-oxidizing bacteria encode and highly express genes for carbon fixation (RuBisCO), nitrogen fixation (nifHDK) and sulfur oxidation (oxidative-dsrAB), especially in roots stressed by sulfidic and reduced sediment conditions. Stressed roots exhibited the highest rates of nitrogen fixation and expression level of sulfur oxidation and sulfate reduction genes. Close relatives of marine symbionts from the Candidatus Thiodiazotropha genus contributed ~30% and ~20% of all sulfur-oxidizing dsrA and nitrogen-fixing nifK transcripts in stressed roots, respectively. Based on these findings, we propose that the symbiosis between S. alterniflora and sulfur-oxidizing bacteria is key to ecosystem functioning of coastal salt marshes.

Phylogeographic Pattern of the Assassin Bug Sycanus bifidus Inferred from Mitochondrial Genomes and Nuclear Genes

The assassin bug Sycanus bifidus has a wide distribution across southern China. This study explored its distribution and evolution by analyzing mitochondrial and nuclear ribosomal RNA genes, revealing how Pleistocene climate and geological changes shaped its phylogeography. We identified two main clades, A and B, that diverged in the Middle Pleistocene. Hainan Island's populations form a unique group within Clade A, suggesting that the Qiongzhou Strait served as a dispersal corridor during glaciation. Rising sea levels likely separated the Hainan population afterward. Ecological niche modeling showed that both populations have been viable since the last interglacial period, with demographic analyses indicating possible expansions during the Middle and Late Pleistocene, driven by favorable climates. This study highlights the significant effects of Pleistocene sea-level and climatic changes on the distribution and evolution of S. bifidus in China.

Phytoplankton-derived polysaccharides and microbial peptidoglycans are key nutrients for deep-sea microbes in the Mariana Trench

BACKGROUND

The deep sea represents the largest marine ecosystem, driving global-scale biogeochemical cycles. Microorganisms are the most abundant biological entities and play a vital role in the cycling of organic matter in such ecosystems. The primary food source for abyssal biota is the sedimentation of particulate organic polymers. However, our knowledge of the specific biopolymers available to deep-sea microbes remains largely incomplete. One crucial rate-limiting step in organic matter cycling is the depolymerization of particulate organic polymers facilitated by extracellular enzymes (EEs). Therefore, the investigation of active EEs and the microbes responsible for their production is a top priority to better understand the key nutrient sources for deep-sea microbes.

RESULTS

In this study, we conducted analyses of extracellular enzymatic activities (EEAs), metagenomics, and metatranscriptomics from seawater samples of 50-9305 m from the Mariana Trench. While a diverse array of microbial groups was identified throughout the water column, only a few exhibited high levels of transcriptional activities. Notably, microbial populations actively transcribing EE genes involved in biopolymer processing in the abyssopelagic (4700 m) and hadopelagic zones (9305 m) were primarily associated with the class Actinobacteria. These microbes actively transcribed genes coding for enzymes such as cutinase, laccase, and xyloglucanase which are capable of degrading phytoplankton polysaccharides as well as GH23 peptidoglycan lyases and M23 peptidases which have the capacity to break down peptidoglycan. Consequently, corresponding enzyme activities including glycosidases, esterase, and peptidases can be detected in the deep ocean. Furthermore, cell-specific EEAs increased at 9305 m compared to 4700 m, indicating extracellular enzymes play a more significant role in nutrient cycling in the deeper regions of the Mariana Trench.

CONCLUSIONS

Transcriptomic analyses have shed light on the predominant microbial population actively participating in organic matter cycling in the deep-sea environment of the Mariana Trench. The categories of active EEs suggest that the complex phytoplankton polysaccharides (e.g., cutin, lignin, and hemicellulose) and microbial peptidoglycans serve as the primary nutrient sources available to deep-sea microbes. The high cell-specific EEA observed in the hadal zone underscores the robust polymer-degrading capacities of hadal microbes even in the face of the challenging conditions they encounter in this extreme environment. These findings provide valuable new insights into the sources of nutrition, the key microbes, and the EEs crucial for biopolymer degradation in the deep seawater of the Mariana Trench. Video Abstract.

Evolution of myxozoan mitochondrial genomes: insights from myxobolids

BACKGROUND

Myxozoa is a class of cnidarian parasites that encompasses over 2,400 species. Phylogenetic relationships among myxozoans remain highly debated, owing to both a lack of informative morphological characters and a shortage of molecular markers. Mitochondrial (mt) genomes are a common marker in phylogeny and biogeography. However, only five complete myxozoan mt genomes have been sequenced: four belonging to two closely related genera, Enteromyxum and Kudoa, and one from the genus Myxobolus. Interestingly, while cytochrome oxidase genes could be identified in Enteromyxum and Kudoa, no such genes were found in Myxobolus squamalis, and another member of the Myxobolidae (Henneguya salminicola) was found to have lost its entire mt genome. To evaluate the utility of mt genomes to reconstruct myxozoan relationships and to understand if the loss of cytochrome oxidase genes is a characteristic of myxobolids, we sequenced the mt genome of five myxozoans (Myxobolus wulii, M. honghuensis, M. shantungensis, Thelohanellus kitauei and, Sphaeromyxa zaharoni) using Illumina and Oxford Nanopore platforms.

RESULTS

Unlike Enteromyxum, which possesses a partitioned mt genome, the five mt genomes were encoded on single circular chromosomes. An mt plasmid was found in M. wulii, as described previously in Kudoa iwatai. In all new myxozoan genomes, five protein-coding genes (cob, cox1, cox2, nad1, and nad5) and two rRNAs (rnl and rns) were recognized, but no tRNA. We found that Myxobolus and Thelohanellus species shared unidentified reading frames, supporting the view that these mt open reading frames are functional. Our phylogenetic reconstructions based on the five conserved mt genes agree with previously published trees based on the 18S rRNA gene.

CONCLUSIONS

Our results suggest that the loss of cytochrome oxidase genes is not a characteristic of all myxobolids, the ancestral myxozoan mt genome was likely encoded on a single circular chromosome, and mt plasmids exist in a few lineages. Our findings indicate that myxozoan mt sequences are poor markers for reconstructing myxozoan phylogenetic relationships because of their fast-evolutionary rates and the abundance of repeated elements, which complicates assembly.

Potential routes of plastics biotransformation involving novel plastizymes revealed by global multi-omic analysis of plastic associated microbes

Despite increasing efforts across various disciplines, the fate, transport, and impact of synthetic plastics on the environment and public health remain poorly understood. To better elucidate the microbial ecology of plastic waste and its potential for biotransformation, we conducted a large-scale analysis of all publicly available meta-omic studies investigating plastics (n = 27) in the environment. Notably, we observed low prevalence of known plastic degraders throughout most environments, except for substantial enrichment in riverine systems. This indicates rivers may be a highly promising environment for discovery of novel plastic bioremediation products. Ocean samples associated with degrading plastics showed clear differentiation from non-degrading polymers, showing enrichment of novel putative biodegrading taxa in the degraded samples. Regarding plastisphere pathogenicity, we observed significant enrichment of antimicrobial resistance genes on plastics but not of virulence factors. Additionally, we report a co-occurrence network analysis of 10 + million proteins associated with the plastisphere. This analysis revealed a localized sub-region enriched with known and putative plastizymes-these may be useful for deeper investigation of nature's ability to biodegrade man-made plastics. Finally, the combined data from our meta-analysis was used to construct a publicly available database, the Plastics Meta-omic Database (PMDB)-accessible at plasticmdb.org. These data should aid in the integrated exploration of the microbial plastisphere and facilitate research efforts investigating the fate and bioremediation potential of environmental plastic waste.

Metagenomic sequencing of a Patescibacteria-containing enrichment from Zodletone spring in Oklahoma, USA

An enrichment of sulfidic sediments from Zodletone spring was sequenced as a metagenome. Draft genomes representing Cloacimonadota, Deltabacterota, Firmicutes, and Patescibacteria were binned and annotated and will aid functional genomics and cultivation efforts.

Metagenomic insights into the impact of litter from poultry Concentrated Animal Feeding Operations (CAFOs) to adjacent soil and water microbial communities

In recent decades, human food consumption has led to an increased demand for animal-based foods, particularly chicken meat production. The state of Georgia, USA is one of the top broiler chicken producers in the United States, where animals are raised in Concentrated Animal Feeding Operations (CAFOs). Without proper management, CAFOs could negatively impact the environment and become a public health risk as a source of water and air pollution and/or by spreading antimicrobial resistance genes. In this study, we used metagenome sequencing to investigate the impact of the application of the CAFO's litter on adjacent soils and downstream creek waters in terms of microbial diversity and antimicrobial resistance profile changes. Our data indicate that while a few microbial groups increased in abundance within a short period of time after litter application, these populations subsequently decreased to levels similar to those found prior to the litter application or to below the detection limit of our metagenome sequencing effort. Microbial taxonomic composition analyses, relative abundance of Metagenome-Assembled Genomes (MAGs) and detection of Antimicrobial Resistance Genes (ARGs) allow us to conclude that this practice of litter application had a negligible effect on the microbiome or resistome profile of these soils and nearby waterways, likely due to its dilution in the field and/or outcompetition by indigenous microbes, revealing a minimal impact of these poultry facilities on the natural microbial communities.

The ecological role of microbiome at community-, taxonomic - and genome-levels in black-odorous waters

Black-odorous waters (BOWs) are heavily polluted waters where microbial information remains elusive mechanistically. Based on gene amplicon and metagenomics sequencing, a comprehensive study was conducted to investigate the microbial communities in urban and rural BOWs. The results revealed that microbial communities' assembly in urban and rural BOWs was predominantly governed by stochastic factors at the community level. At the taxonomic level, there were 62 core species (58.48%) in water and 207 core species (44.56%) in sediment across urban and rural areas. Notably, significant differences were observed in the functional genetic composition of BOWs between urban and rural areas. Specifically, rural areas exhibited an enhanced abundance of genes involved in nitrogen fixation, Fe2+ transport, and sulfate reduction. Conversely, urban areas showed higher abundances of some genes associated with carbon fixation, nitrification and denitrification. A sulfur-centered ecological model of microbial communities was constructed by integrating data from the three levels of analysis, and 14 near-complete draft genomes were generated, representing a substantial portion of the microbial community (35.04% in rural BOWs and 29.97% in urban BOWs). This research provides significant insights into the sustainable management and preservation of aquatic ecosystems affected by BOWs.

Screening potential biomarkers associated with insulin resistance in high-fat diet-fed mice by integrating metagenomics and untargeted metabolomics

Insulin resistance is the primary pathophysiological basis for metabolic syndrome and type 2 diabetes. Gut microbiota and microbiota-derived metabolites are pivotal in insulin resistance. However, identifying the specific microbes and key metabolites with causal roles is a challenging task, and the underlying mechanisms require further exploration. Here, we successfully constructed a model of insulin resistance in mice induced by a high-fat diet (HFD) and screened potential biomarkers associated with insulin resistance by integrating metagenomics and untargeted metabolomics. Our findings showed a significant increase in the abundance of 30 species of Alistipes in HFD mice compared to normal diet (ND) mice, while the abundance of Desulfovibrio and Candidatus Amulumruptor was significantly lower in HFD mice than in ND mice. Non-targeted metabolomics analysis identified 21 insulin resistance-associated metabolites, originating from the microbiota or co-metabolized by both the microbiota and the host. These metabolites were primarily enriched in aromatic amino acid metabolism (tryptophan metabolism, tyrosine metabolism, and phenylalanine metabolism) and arginine biosynthesis. Further analysis revealed a significant association between the three distinct genera and 21 differentiated metabolites in the HFD and ND mice. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of representative genomes from 12 species of the three distinct genera further revealed the functional potential in aromatic amino acid metabolism and arginine biosynthesis. This study lays the groundwork for future investigations into the mechanisms through which the gut microbiota and its metabolites impact insulin resistance.

IMPORTANCE

In this study, we aim to identify the microbes and metabolites linked to insulin resistance, some of which have not been previously reported in insulin resistance-related studies. This adds a complementary dimension to existing research. Furthermore, we establish a correlation between alterations in the gut microbiota and metabolite levels. These findings serve as a foundation for identifying the causal bacterial species and metabolites. They also offer insights that guide further exploration into the mechanisms through which these factors influence host insulin resistance.

Metagenomic evaluation of the performance of passive Moore swabs for sewage monitoring relative to composite sampling over time resolved deployments

Moore swabs have re-emerged as a versatile tool in the field of wastewater-based epidemiology during the COVID-19 pandemic and offer unique advantages for monitoring pathogens in sewer systems, especially at the neighborhood-level. However, whether Moore swabs provide comparable results to more commonly used composite samples remains to be rigorously tested including the optimal duration of Moore swab deployment. This study provides new insights into these issues by comparing the results from Moore swab samples to those of paired composite samples collected from the same sewer lines continuously over six to seventy-two hours post-deployment, during low COVID-19 prevalence periods. Our results show that Moore swabs accumulated approximately 10-fold higher PMMoV concentrations (on a basis of mL of Moore swab squeezed filtrate to mL of composite sewage) and showed comparable trends in terms of bacterial species abundance when compared to composite samples. Moore swabs also generally captured higher SARS-CoV-2 N1/N2 RNA concentrations than composite samples. Moore swabs showed comparable trends in terms of abundance dynamics of the sewage microbiome to composite samples and variable signs of saturation over time that were site and/or microbial population-specific. Based on our dual ddRT-PCR and shotgun metagenomic approach, we find that Moore swabs at our sites were optimally deployed for 6 h at a time at two sites.

Convergent photophysiology and prokaryotic assemblage structure in epilithic cyanobacterial tufts and algal turf communities

As global change spurs shifts in benthic community composition on coral reefs globally, a better understanding of the defining taxonomic and functional features that differentiate proliferating benthic taxa is needed to predict functional trajectories of reef degradation better. This is especially critical for algal groups, which feature dramatically on changing reefs. Limited attention has been given to characterizing the features that differentiate tufting epilithic cyanobacterial communities from ubiquitous turf algal assemblages. Here, we integrated an in situ assessment of photosynthetic yield with metabarcoding and shotgun metagenomic sequencing to explore photophysiology and prokaryotic assemblage structure within epilithic tufting benthic cyanobacterial communities and epilithic algal turf communities. Significant differences were not detected in the average quantum yield. However, variability in yield was significantly higher in cyanobacterial tufts. Neither prokaryotic assemblage diversity nor structure significantly differed between these functional groups. The sampled cyanobacterial tufts, predominantly built by Okeania sp., were co-dominated by members of the Proteobacteria, Firmicutes, and Bacteroidota, as were turf algal communities. Few detected ASVs were significantly differentially abundant between functional groups and consisted exclusively of taxa belonging to the phyla Proteobacteria and Firmicutes. Assessment of the distribution of recovered cyanobacterial amplicons demonstrated that alongside sample-specific cyanobacterial diversification, the dominant cyanobacterial members were conserved across tufting cyanobacterial and turf algal communities. Overall, these data suggest a convergence in taxonomic identity and mean photosynthetic potential between tufting epilithic cyanobacterial communities and algal turf communities, with numerous implications for consumer-resource dynamics on future reefs and trajectories of reef functional ecology.

Efficient Recovery of Complete Gut Viral Genomes by Combined Short- and Long-Read Sequencing

Current metagenome assembled human gut phage catalogs contained mostly fragmented genomes. Here, comprehensive gut virome detection procedure is developed involving virus-like particle (VLP) enrichment from ≈500 g feces and combined sequencing of short- and long-read. Applied to 135 samples, a Chinese Gut Virome Catalog (CHGV) is assembled consisting of 21,499 non-redundant viral operational taxonomic units (vOTUs) that are significantly longer than those obtained by short-read sequencing and contained ≈35% (7675) complete genomes, which is ≈nine times more than those in the Gut Virome Database (GVD, ≈4%, 1,443). Interestingly, the majority (≈60%, 13,356) of the CHGV vOTUs are obtained by either long-read or hybrid assemblies, with little overlap with those assembled from only the short-read data. With this dataset, vast diversity of the gut virome is elucidated, including the identification of 32% (6,962) novel vOTUs compare to public gut virome databases, dozens of phages that are more prevalent than the crAssphages and/or Gubaphages, and several viral clades that are more diverse than the two. Finally, the functional capacities are also characterized of the CHGV encoded proteins and constructed a viral-host interaction network to facilitate future research and applications.

A pan-genomic approach reveals novel Sulfurimonas clade in the ferruginous meromictic Lake Pavin

The permanently anoxic waters in meromictic lakes create suitable niches for the growth of bacteria using sulphur metabolisms like sulphur oxidation. In Lake Pavin, the anoxic water mass hosts an active cryptic sulphur cycle that interacts narrowly with iron cycling, however the metabolisms of the microorganisms involved are poorly known. Here we combined metagenomics, single-cell genomics, and pan-genomics to further expand our understanding of the bacteria and the corresponding metabolisms involved in sulphur oxidation in this ferruginous sulphide- and sulphate-poor meromictic lake. We highlighted two new species within the genus Sulfurimonas that belong to a novel clade of chemotrophic sulphur oxidisers exclusive to freshwaters. We moreover conclude that this genus holds a key-role not only in limiting sulphide accumulation in the upper part of the anoxic layer but also constraining carbon, phosphate and iron cycling.

Correlation between the gut microbiome and neurodegenerative diseases: a review of metagenomics evidence

A growing body of evidence suggests that the gut microbiota contributes to the development of neurodegenerative diseases via the microbiota-gut-brain axis. As a contributing factor, microbiota dysbiosis always occurs in pathological changes of neurodegenerative diseases, such as Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis. High-throughput sequencing technology has helped to reveal that the bidirectional communication between the central nervous system and the enteric nervous system is facilitated by the microbiota's diverse microorganisms, and for both neuroimmune and neuroendocrine systems. Here, we summarize the bioinformatics analysis and wet-biology validation for the gut metagenomics in neurodegenerative diseases, with an emphasis on multi-omics studies and the gut virome. The pathogen-associated signaling biomarkers for identifying brain disorders and potential therapeutic targets are also elucidated. Finally, we discuss the role of diet, prebiotics, probiotics, postbiotics and exercise interventions in remodeling the microbiome and reducing the symptoms of neurodegenerative diseases.

Insights into the Adolescent Cystic Fibrosis Airway Microbiome Using Shotgun Metagenomics

Cystic fibrosis (CF) is an inherited genetic disorder which manifests primarily in airway disease. Recent advances in molecular technologies have unearthed the diverse polymicrobial nature of the CF airway. Numerous studies have characterised the genus-level composition of this airway community using targeted 16S rDNA sequencing. Here, we employed whole-genome shotgun metagenomics to provide a more comprehensive understanding of the early CF airway microbiome. We collected 48 sputum samples from 11 adolescents and children with CF over a 12-month period and performed shotgun metagenomics on the Illumina NextSeq platform. We carried out functional and taxonomic analysis of the lung microbiome at the species and strain levels. Correlations between microbial diversity measures and independent demographic and clinical variables were performed. Shotgun metagenomics detected a greater diversity of bacteria than culture-based methods. A large proportion of the top 25 most-dominant species were anaerobes. Samples dominated by Staphylococcus aureus and Prevotella melaninogenica had significantly higher microbiome diversity, while no CF pathogen was associated with reduced microbial diversity. There was a diverse resistome present in all samples in this study, with 57.8% agreement between shotgun metagenomics and culture-based methods for detection of resistance. Pathogenic sequence types (STs) of S. aureus, Pseudomonas aeruginosa, Haemophilus influenzae and Stenotrophomonas maltophilia were observed to persist in young CF patients, while STs of S. aureus were both persistent and shared between patients. This study provides new insight into the temporal changes in strain level composition of the microbiome and the landscape of the resistome in young people with CF. Shotgun metagenomics could provide a very useful one-stop assay for detecting pathogens, emergence of resistance and conversion to persistent colonisation in early CF disease.

Rescue of Mycobacterium bovis DNA Obtained from Cultured Samples during Official Surveillance of Animal TB: Key Steps for Robust Whole Genome Sequence Data Generation

Epidemiological surveillance of animal tuberculosis (TB) based on whole genome sequencing (WGS) of Mycobacterium bovis has recently gained track due to its high resolution to identify infection sources, characterize the pathogen population structure, and facilitate contact tracing. However, the workflow from bacterial isolation to sequence data analysis has several technical challenges that may severely impact the power to understand the epidemiological scenario and inform outbreak response. While trying to use archived DNA from cultured samples obtained during routine official surveillance of animal TB in Portugal, we struggled against three major challenges: the low amount of M. bovis DNA obtained from routinely processed animal samples; the lack of purity of M. bovis DNA, i.e., high levels of contamination with DNA from other organisms; and the co-occurrence of more than one M. bovis strain per sample (within-host mixed infection). The loss of isolated genomes generates missed links in transmission chain reconstruction, hampering the biological and epidemiological interpretation of data as a whole. Upon identification of these challenges, we implemented an integrated solution framework based on whole genome amplification and a dedicated computational pipeline to minimize their effects and recover as many genomes as possible. With the approaches described herein, we were able to recover 62 out of 100 samples that would have otherwise been lost. Based on these results, we discuss adjustments that should be made in official and research laboratories to facilitate the sequential implementation of bacteriological culture, PCR, downstream genomics, and computational-based methods. All of this in a time frame supporting data-driven intervention.

To live or die: "Fine-tuning" adaptation revealed by systemic analyses in symbiotic bathymodiolin mussels from diverse deep-sea extreme ecosystems

Hydrothermal vents (HVs) and cold seeps (CSs) are typical deep-sea extreme ecosystems with their own geochemical characteristics to supply the unique living conditions for local communities. Once HVs or CSs stop emission, the dramatic environmental change would pose survival risks to deep-sea organisms. Up to now, limited knowledge has been available to understand the biological responses and adaptive strategy to the extreme environments and their transition from active to extinct stage, mainly due to the technical difficulties and lack of representative organisms. In this study, bathymodiolin mussels, the dominant and successful species surviving in diverse deep-sea extreme ecosystems, were collected from active and extinct HVs (Southwest Indian Ocean) or CSs (South China Sea) via two individual cruises. The transcriptomic analysis and determination of multiple biological indexes in stress defense and metabolic systems were conducted in both gills and digestive glands of mussels, together with the metagenomic analysis of symbionts in mussels. The results revealed the ecosystem- and tissue-specific transcriptional regulation in mussels, addressing the autologous adaptations in antioxidant defense, energy utilization and key compounds (i.e. sulfur) metabolism. In detail, the successful antioxidant defense contributed to conquering the oxidative stress induced during the unavoidable metabolism of xenobiotics commonly existing in the extreme ecosystems; changes in metabolic rate functioned to handle toxic matters in different surroundings; upregulated gene expression of sulfide:quinone oxidoreductase indicated an active sulfide detoxification in mussels from HVs and active stage of HVs & CSs. Coordinately, a heterologous adaptation, characterized by the functional compensation between symbionts and mussels in energy utilization, sulfur and carbon metabolism, was also evidenced by the bacterial metagenomic analysis. Taken together, a new insight was proposed that symbiotic bathymodiolin mussels would develop a "finetuning" strategy combining the autologous and heterologous regulations to fulfill the efficient and effective adaptations for successful survival.

Metagenomic assembled genomes indicated the potential application of hypersaline microbiome for plant growth promotion and stress alleviation in salinized soils

Climate change is causing unpredictable seasonal variations globally. Due to the continuously increasing earth's surface temperature, the rate of water evaporation is enhanced, conceiving a problem of soil salinization, especially in arid and semi-arid regions. The accumulation of salt degrades soil quality, impairs plant growth, and reduces agricultural yields. Salt-tolerant, plant-growth-promoting microorganisms may offer a solution, enhancing crop productivity and soil fertility in salinized areas. In the current study, genome-resolved metagenomic analysis has been performed to investigate the salt-tolerating and plant growth-promoting potential of two hypersaline ecosystems, Sambhar Lake and Drang Mine. The samples were co-assembled independently by Megahit, MetaSpades, and IDBA-UD tools. A total of 67 metagenomic assembled genomes (MAGs) were reconstructed following the binning process, including 15 from Megahit, 26 from MetaSpades, and 26 from IDBA_UD assembly tools. As compared to other assemblers, the MAGs obtained by MetaSpades were of superior quality, with a completeness range of 12.95%-96.56% and a contamination range of 0%-8.65%. The medium and high-quality MAGs from MetaSpades, upon functional annotation, revealed properties such as salt tolerance (91.3%), heavy metal tolerance (95.6%), exopolysaccharide (95.6%), and antioxidant (60.86%) biosynthesis. Several plant growth-promoting attributes, including phosphate solubilization and indole-3-acetic acid (IAA) production, were consistently identified across all obtained MAGs. Conversely, characteristics such as iron acquisition and potassium solubilization were observed in a substantial majority, specifically 91.3%, of the MAGs. The present study indicates that hypersaline microflora can be used as bio-fertilizing agents for agricultural practices in salinized areas by alleviating prevalent stresses.

IMPORTANCE

The strategic implementation of metagenomic assembled genomes (MAGs) in exploring the properties and harnessing microorganisms from ecosystems like hypersaline niches has transformative potential in agriculture. This approach promises to redefine our comprehension of microbial diversity and its ecosystem roles. Recovery and decoding of MAGs unlock genetic resources, enabling the development of new solutions for agricultural challenges. Enhanced understanding of these microbial communities can lead to more efficient nutrient cycling, pest control, and soil health maintenance. Consequently, traditional agricultural practices can be improved, resulting in increased yields, reduced environmental impacts, and heightened sustainability. MAGs offer a promising avenue for sustainable agriculture, bridging the gap between cutting-edge genomics and practical field applications.

BASALT refines binning from metagenomic data and increases resolution of genome-resolved metagenomic analysis

Metagenomic binning is an essential technique for genome-resolved characterization of uncultured microorganisms in various ecosystems but hampered by the low efficiency of binning tools in adequately recovering metagenome-assembled genomes (MAGs). Here, we introduce BASALT (Binning Across a Series of Assemblies Toolkit) for binning and refinement of short- and long-read sequencing data. BASALT employs multiple binners with multiple thresholds to produce initial bins, then utilizes neural networks to identify core sequences to remove redundant bins and refine non-redundant bins. Using the same assemblies generated from Critical Assessment of Metagenome Interpretation (CAMI) datasets, BASALT produces up to twice as many MAGs as VAMB, DASTool, or metaWRAP. Processing assemblies from a lake sediment dataset, BASALT produces ~30% more MAGs than metaWRAP, including 21 unique class-level prokaryotic lineages. Functional annotations reveal that BASALT can retrieve 47.6% more non-redundant opening-reading frames than metaWRAP. These results highlight the robust handling of metagenomic sequencing data of BASALT.

Host genotype-specific rhizosphere fungus enhances drought resistance in wheat

BACKGROUND

The severity and frequency of drought are expected to increase substantially in the coming century and dramatically reduce crop yields. Manipulation of rhizosphere microbiomes is an emerging strategy for mitigating drought stress in agroecosystems. However, little is known about the mechanisms underlying how drought-resistant plant recruitment of specific rhizosphere fungi enhances drought adaptation of drought-sensitive wheats. Here, we investigated microbial community assembly features and functional profiles of rhizosphere microbiomes related to drought-resistant and drought-sensitive wheats by amplicon and shotgun metagenome sequencing techniques. We then established evident linkages between root morphology traits and putative keystone taxa based on microbial inoculation experiments. Furthermore, root RNA sequencing and RT-qPCR were employed to explore the mechanisms how rhizosphere microbes modify plant response traits to drought stresses.

RESULTS

Our results indicated that host plant signature, plant niche compartment, and planting site jointly contribute to the variation of soil microbiome assembly and functional adaptation, with a relatively greater effect of host plant signature observed for the rhizosphere fungi community. Importantly, drought-resistant wheat (Yunhan 618) possessed more diverse bacterial and fungal taxa than that of the drought-sensitive wheat (Chinese Spring), particularly for specific fungal species. In terms of microbial interkingdom association networks, the drought-resistant variety possessed more complex microbial networks. Metagenomics analyses further suggested that the enriched rhizosphere microbiomes belonging to the drought-resistant cultivar had a higher investment in energy metabolism, particularly in carbon cycling, that shaped their distinctive drought tolerance via the mediation of drought-induced feedback functional pathways. Furthermore, we observed that host plant signature drives the differentiation in the ecological role of the cultivable fungal species Mortierella alpine (M. alpina) and Epicoccum nigrum (E. nigrum). The successful colonization of M. alpina on the root surface enhanced the resistance of wheats in response to drought stresses via activation of drought-responsive genes (e.g., CIPK9 and PP2C30). Notably, we found that lateral roots and root hairs were significantly suppressed by co-colonization of a drought-enriched fungus (M. alpina) and a drought-depleted fungus (E. nigrum).

CONCLUSIONS

Collectively, our findings revealed host genotypes profoundly influence rhizosphere microbiome assembly and functional adaptation, as well as it provides evidence that drought-resistant plant recruitment of specific rhizosphere fungi enhances drought tolerance of drought-sensitive wheats. These findings significantly underpin our understanding of the complex feedbacks between plants and microbes during drought, and lay a foundation for steering "beneficial keystone biome" to develop more resilient and productive crops under climate change. Video Abstract.

Comparative Analysis of the Mitochondrial Genomes of Chloropidae and Their Implications for the Phylogeny of the Family

Chloropidae, commonly known as grass flies, represent the most taxonomically diverse family of Diptera Carnoidea, comprising over 3000 described species worldwide. Previous phylogenetic studies of this family have predominantly relied on morphological characters, with mitochondrial genomes being reported in a few species. This study presents 11 newly sequenced mitochondrial genomes (10 Chloropidae and 1 Milichiidae) and provides the first comprehensive comparative analysis of mitochondrial genomes for Chloropidae. Apart from 37 standard genes and the control region, three conserved intergenic sequences across Diptera Cyclorrhapha were identified in all available chloropid mitochondrial genomes. Evolutionary rates within Chloropidae exhibit significant variation across subfamilies, with Chloropinae displaying higher rates than the other three subfamilies. Phylogenetic relationships based on mitochondrial genomes were inferred using maximum likelihood and Bayesian methods. The monophyly of Chloropidae and all four subfamilies is consistently strongly supported, while subfamily relationships within Chloropidae remain poorly resolved, possibly due to rapid evolution.

A novel and diverse family of filamentous DNA viruses associated with parasitic wasps

Large dsDNA viruses from the Naldaviricetes class are currently composed of four viral families infecting insects and/or crustaceans. Since the 1970s, particles described as filamentous viruses (FVs) have been observed by electronic microscopy in several species of Hymenoptera parasitoids but until recently, no genomic data was available. This study provides the first comparative morphological and genomic analysis of these FVs. We analyzed the genomes of seven FVs, six of which were newly obtained, to gain a better understanding of their evolutionary history. We show that these FVs share all genomic features of the Naldaviricetes while encoding five specific core genes that distinguish them from their closest relatives, the Hytrosaviruses. By mining public databases, we show that FVs preferentially infect Hymenoptera with parasitoid lifestyle and that these viruses have been repeatedly integrated into the genome of many insects, particularly Hymenoptera parasitoids, overall suggesting a long-standing specialization of these viruses to parasitic wasps. Finally, we propose a taxonomical revision of the class Naldaviricetes in which FVs related to the Leptopilina boulardi FV constitute a fifth family. We propose to name this new family, Filamentoviridae.

Microbiologic surveys for Baijiu fermentation are affected by experimental design

The traditional Chinese alcoholic beverage Baijiu is produced by spontaneous fermentation of grains under anaerobic conditions. While numerous studies have used metagenomic technology to investigate the microbiome of Baijiu brewing, the microbial succession mechanism of Baijiu brewing has not been fully clarified, and metagenomic strategies for microecology surveys have not been comprehensively evaluated. Using the fermentation process of strong-flavor Baijiu as a model, we compared the data for bacterial communities based on short read 16S rRNA variable regions, V3-V4, and full-length 16S regions, V1-V9, to whole metagenomic shotgun sequencing (WMS) to measure the effect of technology selection on phylogenetic and functional profiles. The results showed differences in bacterial compositions and their relation to volatiles and physicochemical variables between sequencing methods. Furthermore, the percentage of V3-V4 sequences assigned to species level was higher than the percentage of V1-V9 sequences according to SILVA taxonomy, but lower according to NCBI taxonomy (P < 0.05). In both SILVA and NCBI taxonomies, the relative abundances of bacterial communities at both the genus and family levels were more positively correlated with WMS data in the V3-V4 dataset than in the V1-V9 dataset. The WMS identified changes in abundances of multiple metabolic pathways during fermentation (P < 0.05), including "starch and sucrose metabolism," "galactose metabolism," and "fatty acid biosynthesis." Although functional predictions derived from 16S data show similar patterns to WMS, most metabolic pathway changes are uncorrelated (P > 0.05). This study provided new technical and biological insights into Baijiu production that may assist in selection of methodologies for studies of fermentation systems.

COBRA improves the completeness and contiguity of viral genomes assembled from metagenomes

Viruses are often studied using metagenome-assembled sequences, but genome incompleteness hampers comprehensive and accurate analyses. Contig Overlap Based Re-Assembly (COBRA) resolves assembly breakpoints based on the de Bruijn graph and joins contigs. Here we benchmarked COBRA using ocean and soil viral datasets. COBRA accurately joined the assembled sequences and achieved notably higher genome accuracy than binning tools. From 231 published freshwater metagenomes, we obtained 7,334 bacteriophage clusters, ~83% of which represent new phage species. Notably, ~70% of these were circular, compared with 34% before COBRA analyses. We expanded sampling of huge phages (≥200 kbp), the largest of which was curated to completion (717 kbp). Improved phage genomes from Rotsee Lake provided context for metatranscriptomic data and indicated the in situ activity of huge phages, whiB-encoding phages and cysC- and cysH-encoding phages. COBRA improves viral genome assembly contiguity and completeness, thus the accuracy and reliability of analyses of gene content, diversity and evolution.

Temporal differentiation in the adaptation of functional bacteria to low-temperature stress in partial denitrification and anammox system

Despite reliable nitrite supply through partial denitrification, the adaptation of denitrifying bacteria to low temperatures remains elusive in partial denitrification and anammox (PDA) systems. Here, temporal differentiations of the structure, activity, and relevant cold-adaptation mechanism of functional bacteria were investigated in a lab-scale PDA bioreactor at decreased temperature. Although distinct denitrifying bacteria dominated after low-temperature stress, both short- and long-term stresses exerted differential selectivity towards the species with close phylogenetic distance. Species Azonexus sp.149 showed high superiority over Azonexus sp.384 under short-term stress, and long-term stress improved the adaptation of Aquabacterium sp.93 instead of Aquabacterium sp.184. The elevated transcription of nitrite reductase genes suggested that several denitrifying bacteria (e.g., Azonexus sp.149) could compete with anammox bacteria for nitrite. Species Rivicola pingtungensis and Azonexus sp.149 could adapt through various adaptation pathways, such as the two-component system, cold shock protein (CSP), membrane alternation, and electron transport chain. By contrast, species Zoogloea sp.273 and Aquabacterium sp.93 mainly depended on the CSP and oxidative stress response. This study largely deepens our understanding of the performance deterioration in PDA systems during cold shock and provides several references for efficient adaptation to seasonal temperature fluctuation.

Resistant starch intake facilitates weight loss in humans by reshaping the gut microbiota

Emerging evidence suggests that modulation of gut microbiota by dietary fibre may offer solutions for metabolic disorders. In a randomized placebo-controlled crossover design trial (ChiCTR-TTRCC-13003333) in 37 participants with overweight or obesity, we test whether resistant starch (RS) as a dietary supplement influences obesity-related outcomes. Here, we show that RS supplementation for 8 weeks can help to achieve weight loss (mean -2.8 kg) and improve insulin resistance in individuals with excess body weight. The benefits of RS are associated with changes in gut microbiota composition. Supplementation with Bifidobacterium adolescentis, a species that is markedly associated with the alleviation of obesity in the study participants, protects male mice from diet-induced obesity. Mechanistically, the RS-induced changes in the gut microbiota alter the bile acid profile, reduce inflammation by restoring the intestinal barrier and inhibit lipid absorption. We demonstrate that RS can facilitate weight loss at least partially through B. adolescentis and that the gut microbiota is essential for the action of RS.

Diverse plasmid systems and their ecology across human gut metagenomes revealed by PlasX and MobMess

Plasmids alter microbial evolution and lifestyles by mobilizing genes that often confer fitness in changing environments across clades. Yet our ecological and evolutionary understanding of naturally occurring plasmids is far from complete. Here we developed a machine-learning model, PlasX, which identified 68,350 non-redundant plasmids across human gut metagenomes and organized them into 1,169 evolutionarily cohesive 'plasmid systems' using our sequence containment-aware network-partitioning algorithm, MobMess. Individual plasmids were often country specific, yet most plasmid systems spanned across geographically distinct human populations. Cargo genes in plasmid systems included well-known determinants of fitness, such as antibiotic resistance, but also many others including enzymes involved in the biosynthesis of essential nutrients and modification of transfer RNAs, revealing a wide repertoire of likely fitness determinants in complex environments. Our study introduces computational tools to recognize and organize plasmids, and uncovers the ecological and evolutionary patterns of diverse plasmids in naturally occurring habitats through plasmid systems.

Assembly collapsing versus heterozygosity oversizing: detection of homokaryotic and heterokaryotic Laccaria trichodermophora strains by hybrid genome assembly

Genome assembly and annotation using short-paired reads is challenging for eukaryotic organisms due to their large size, variable ploidy and large number of repetitive elements. However, the use of single-molecule long reads improves assembly quality (completeness and contiguity), but haplotype duplications still pose assembly challenges. To address the effect of read length on genome assembly quality, gene prediction and annotation, we compared genome assemblers and sequencing technologies with four strains of the ectomycorrhizal fungus Laccaria trichodermophora. By analysing the predicted repertoire of carbohydrate enzymes, we investigated the effects of assembly quality on functional inferences. Libraries were generated using three different sequencing platforms (Illumina Next-Seq, Mi-Seq and PacBio Sequel), and genomes were assembled using single and hybrid assemblies/libraries. Long reads or hybrid assemby resolved the collapsing of repeated regions, but the nuclear heterozygous versions remained unresolved. In dikaryotic fungi, each cell includes two nuclei and each nucleus has differences not only in allelic gene version but also in gene composition and synteny. These heterokaryotic cells produce fragmentation and size overestimation of the genome assembly of each nucleus. Hybrid assembly revealed a wider functional diversity of genomes. Here, several predicted oxidizing activities on glycosyl residues of oligosaccharides and several chitooligosaccharide acetylase activities would have passed unnoticed in short-read assemblies. Also, the size and fragmentation of the genome assembly, in combination with heterozygosity analysis, allowed us to distinguish homokaryotic and heterokaryotic strains isolated from L. trichodermophora fruit bodies.

Update on the proposed minimal standards for the use of genome data for the taxonomy of prokaryotes

The field of microbial taxonomy is dynamic, aiming to provide a stable and contemporary classification system for prokaryotes. Traditionally, reliance on phenotypic characteristics limited the comprehensive understanding of microbial diversity and evolution. The introduction of molecular techniques, particularly DNA sequencing and genomics, has transformed our perception of prokaryotic diversity. In the past two decades, advancements in genome sequencing have transitioned from traditional methods to a genome-based taxonomic framework, not only to define species, but also higher taxonomic ranks. As technology and databases rapidly expand, maintaining updated standards is crucial. This work seeks to revise the 2018 guidelines for applying genome sequencing data in microbial taxonomy, adapting minimal standards and recommendations to reflect technological progress during this period.

A cryptic plasmid is among the most numerous genetic elements in the human gut

Plasmids are extrachromosomal genetic elements that often encode fitness-enhancing features. However, many bacteria carry "cryptic" plasmids that do not confer clear beneficial functions. We identified one such cryptic plasmid, pBI143, which is ubiquitous across industrialized gut microbiomes and is 14 times as numerous as crAssphage, currently established as the most abundant extrachromosomal genetic element in the human gut. The majority of mutations in pBI143 accumulate in specific positions across thousands of metagenomes, indicating strong purifying selection. pBI143 is monoclonal in most individuals, likely due to the priority effect of the version first acquired, often from one's mother. pBI143 can transfer between Bacteroidales, and although it does not appear to impact bacterial host fitness in vivo, it can transiently acquire additional genetic content. We identified important practical applications of pBI143, including its use in identifying human fecal contamination and its potential as an alternative approach to track human colonic inflammatory states.

High Frankia abundance and low diversity of microbial community are associated with nodulation specificity and stability of sea buckthorn root nodule

Introduction

Actinorhizal symbioses are gaining attention due to the importance of symbiotic nitrogen fixation in sustainable agriculture. Sea buckthorn (Hippophae L.) is an important actinorhizal plant, yet research on the microbial community and nitrogen cycling in its nodules is limited. In addition, the influence of environmental differences on the microbial community of sea buckthorn nodules and whether there is a single nitrogen-fixing actinomycete species in the nodules are still unknown.

Methods

We investigated the diversity, community composition, network associations and nitrogen cycling pathways of the microbial communities in the root nodule (RN), nodule surface soil (NS), and bulk soil (BS) of Mongolian sea buckthorn distributed under three distinct ecological conditions in northern China using 16S rRNA gene and metagenomic sequencing. Combined with the data of environmental factors, the effects of environmental differences on different sample types were analyzed.

Results

The results showed that plants exerted a clear selective filtering effect on microbiota, resulting in a significant reduction in microbial community diversity and network complexity from BS to NS to RN. Proteobacteria was the most abundant phylum in the microbiomes of BS and NS. While RN was primarily dominated by Actinobacteria, with Frankia sp. EAN1pec serving as the most dominant species. Correlation analysis indicated that the host determined the microbial community composition in RN, independent of the ecological and geographical environmental changes of the sea buckthorn plantations. Nitrogen cycle pathway analyses showed that RN microbial community primarily functions in nitrogen fixation, and Frankia sp. EAN1pec was a major contributor to nitrogen fixation genes in RN.

Discussion

This study provides valuable insights into the effects of eco-geographical environment on the microbial communities of sea buckthorn RN. These findings further prove that the nodulation specificity and stability of sea buckthorn root and Frankia sp. EAN1pec may be the result of their long-term co-evolution.

Dominant nitrogen metabolisms of a warm, seasonally anoxic freshwater ecosystem revealed using genome resolved metatranscriptomics

Nitrogen (N) availability is one of the principal drivers of primary productivity across aquatic ecosystems. However, the microbial communities and emergent metabolisms that govern N cycling in tropical lakes are both distinct from and poorly understood relative to those found in temperate lakes. This latitudinal difference is largely due to the warm (>20°C) temperatures of tropical lake anoxic hypolimnions (deepest portion of a stratified water column), which result in unique anaerobic metabolisms operating without the temperature constraints found in lakes at temperate latitudes. As such, tropical hypolimnions provide a platform for exploring microbial membership and functional diversity. To better understand N metabolism in warm anoxic waters, we combined measurements of geochemistry and water column thermophysical structure with genome-resolved metatranscriptomic analyses of the water column microbiome in Lake Yojoa, Honduras. We sampled above and below the oxycline in June 2021, when the water column was stratified, and again at the same depths and locations in January 2022, when the water column was mixed. We identified 335 different lineages and significantly different microbiome membership between seasons and, when stratified, between depths. Notably, nrfA (indicative of dissimilatory nitrate reduction to ammonium) was upregulated relative to other N metabolism genes in the June hypolimnion. This work highlights the taxonomic and functional diversity of microbial communities in warm and anoxic inland waters, providing insight into the contemporary microbial ecology of tropical ecosystems as well as inland waters at higher latitudes as water columns continue to warm in the face of global change.IMPORTANCEIn aquatic ecosystems where primary productivity is limited by nitrogen (N), whether continuously, seasonally, or in concert with additional nutrient limitations, increased inorganic N availability can reshape ecosystem structure and function, potentially resulting in eutrophication and even harmful algal blooms. Whereas microbial metabolic processes such as mineralization and dissimilatory nitrate reduction to ammonium increase inorganic N availability, denitrification removes bioavailable N from the ecosystem. Therefore, understanding these key microbial mechanisms is critical to the sustainable management and environmental stewardship of inland freshwater resources. This study identifies and characterizes these crucial metabolisms in a warm, seasonally anoxic ecosystem. Results are contextualized by an ecological understanding of the study system derived from a multi-year continuous monitoring effort. This unique data set is the first of its kind in this largely understudied ecosystem (tropical lakes) and also provides insight into microbiome function and associated taxa in warm, anoxic freshwaters.

Mitogenomes Provide Insights into the Species Boundaries and Phylogenetic Relationships among Three Dolycoris Sloe Bugs (Hemiptera: Pentatomidae) from China

(1) Background: The three sloe bugs, Dolycoris baccarum, Dolycoris indicus, and Dolycoris penicillatus, are found in the Chinese mainland and are morphologically similar. The species boundaries and phylogenetic relationships of the three species remain uncertain; (2) Methods: In this study, we generated multiple mitochondrial genomes (mitogenomes) for each of the three species and conducted comparative mitogenomic analysis, species delimitation, and phylogenetic analysis based on these data; (3) Results: Mitogenomes of the three Dolycoris species are conserved in nucleotide composition, gene arrangement, and codon usage. All protein-coding genes (PCGs) were found to be under purifying selection, and the ND4 evolved at the fastest rate. Most species delimitation analyses based on the COI gene and the concatenated 13 PCGs retrieved three operational taxonomic units (OTUs), which corresponded well with the three Dolycoris species identified based on morphological characters. A clear-cut barcode gap was discovered between the interspecific and intraspecific genetic distances of the three Dolycoris species. Phylogenetic analyses strongly supported the monophyly of Dolycoris, with interspecific relationship inferred as (D. indicus + (D. baccarum + D. penicillatus)); (4) Conclusions: Our study provides the first insight into the species boundaries and phylogenetic relationships of the three Dolycoris species distributed across the Chinese mainland.

Whole-genome resequencing confirms the genetic effects of dams on an endangered fish Hemibagrus guttatus (Siluriformes: Bagridae): A case study in a tributary of the Pearl River

Dam construction in riverine ecosystems has fragmented natural aquatic habitats and has altered environmental conditions. As a result, damming has been demonstrated to threaten aquatic biodiversity by reducing species distribution ranges and hindering gene exchange, leading to the inability to adapt to environmental changes. Knowledge of the contemporary genetic diversity and genetic structure of fish populations that are separated by dams is vital to developing effective conservation strategies, particularly for endangered fish species. We chose the Lianjiang River, a tributary of the Pearl River, as a case study to assess the effects of dams on the genetic diversity and genetic structure of an endangered fish species, Hemibagrus guttatus, using whole-genome resequencing data from 63 fish samples. The results indicated low levels of genetic diversity, high levels of inbreeding and decreasing trend of effective population size in fragmented H. guttatus populations. In addition, there were significant genetic structure and genetic differentiation among populations, suggesting that the dams might have affected H. guttatus populations. Our findings may benefit management and conservation practices for this endangered species that is currently suffering from the effects of dam construction.

Novel insights into the rhizosphere and seawater microbiome of Zostera marina in diverse mariculture zones

BACKGROUND

Seagrasses offer various ecosystem services and possess high levels of primary productivity. However, the development of mariculture has affected the homeostasis of seagrass meadow ecosystems. Plant-microbiome associations are essential for seagrasses health, but little is known about the role of environmental microbiomes and how they affect seagrass in a mariculture environment. In this study, we investigated the influence of mariculture on the rhizosphere and seawater microbiome surrounding Zostera marina and focused on the bacterial, eukaryotic, and fungal components in the composition, diversity, metabolism, and responses to mariculture-related environmental factors.

RESULTS

Significant differences in the composition, richness, diversity, and internal relations of the bacterial community between the seawater and rhizosphere sediment surrounding Z. marina were observed, while differences in the eukaryotic and fungal communities were less significant. More complex bacterial and fungal co-occurrence networks were found in the seawater and rhizosphere sediment of the Saccharina japonica (SJ) and sea cucumber (SC) culture zones. The seawater in the SJ zone had higher levels of dissimilatory and assimilatory nitrate reduction, denitrification, and nitrogen fixation processes than the other three zones. The assimilatory sulfate reduction enzymes were higher in the rhizosphere sediments of the SJ zone than in the other three zones. Tetracycline, sulfonamide, and diaminopyrimidine resistance genes were enriched in the mariculture SJ and SC zones.

CONCLUSIONS

Our findings might contribute to a better understanding of the effects of mariculture on the seagrass and the meadow ecosystems and thus revealing their potential operating mechanisms. These insights may serve to raise awareness of the effects of human activities on natural ecosystems, regulation of antibiotic usage, and environmental restoration. Video Abstract.

The complete mitochondrial genome of Aphidius colemani (Hymenoptera: Braconidae: Aphidiinae)

The genome-level features are crucial genetic resources for species identification and phylogenetic analysis. Here, the complete mitochondrial genome of Aphidius colemani Viereck 1912 (Hymenoptera: Braconidae: Aphidiinae) was sequenced, determined and analyzed. The circular genome is 16,372 bp in length with an overall base composition of 38.9% for A, 46.2% for T, 6.7% for C, and 8.2% for G. The mitochondrial genome of A. colemani contained 13 protein-coding genes that initiated by the ATN codon, 22 transfer RNA genes, two ribosomal RNA genes (rRNAs), and a control region (CR). It shared the same gene arrangement patterns that occurred in two tRNA clusters of trnI-trnQ-trnM and trnW-trnC-trnY with Aphidius gifuensis. Phylogenetic analyses using Bayesian inference and Maximum-likelihood methods supported that the two species of Aphidiinae formed a clade and sister to other subfamilies of Braconidae.

Two new species of Erythroneurini (Hemiptera, Cicadellidae, Typhlocybinae) from southern China based on morphology and complete mitogenomes

Erythroneurine leafhoppers (Hemiptera, Cicadellidae, Typhlocybinae, Erythroneurini) are utilized to resolve the relationship between the four erythroneurine leafhopper (Hemiptera, Cicadellidae, Typhlocybinae, Erythroneurini): Arboridia (Arboridia) rongchangensis sp. nov., Thaia (Thaia) jiulongensis sp. nov., Mitjaevia bifurcata Luo, Song & Song, 2021 and Mitjaevia diana Luo, Song & Song, 2021, the two new species are described and illustrated. The mitochondrial gene sequences of these four species were determined to update the mitochondrial genome database of Erythroneurini. The mitochondrial genomes of four species shared high parallelism in nucleotide composition, base composition and gene order, comprising 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs) and an AT control region, which was consistent with majority of species in Cicadellidae; all genes revealed common trait of a positive AT skew and negative GC skew. The mitogenomes of four species were ultra-conservative in structure, and which isanalogous to that of others in size and A + T content. Phylogenetic trees based on the mitogenome data of these species and another 24 species were built employing the maximum likelihood and Bayesian inference methods. The results indicated that the four species belong to the tribe Erythroneurini, M. diana is the sister-group relationship of M. protuberanta + M. bifurcata. The two species Arboridia (Arboridia) rongchangensis sp. nov. and Thaia (Thaia) jiulongensis sp. nov. also have a relatively close genetic relationship with the genus Mitjaevia.

Genome-Resolved Metagenomics and Denitrifying Strain Isolation Reveal New Insights into Microbial Denitrification in the Deep Vadose Zone

The heavy use of nitrogen fertilizer in intensive agricultural areas often leads to nitrate accumulation in subsurface soil and nitrate contamination in groundwater, which poses a serious risk to public health. Denitrifying microorganisms in the subsoil convert nitrate to gaseous forms of nitrogen, thereby mitigating the leaching of nitrate into groundwater. Here, we investigated denitrifying microorganisms in the deep vadose zone of a typical intensive agricultural area in China through microcosm enrichment, genome-resolved metagenomic analysis, and denitrifying bacteria isolation. A total of 1000 metagenome-assembled genomes (MAGs) were reconstructed, resulting in 98 high-quality, dereplicated MAGs that contained denitrification genes. Among them, 32 MAGs could not be taxonomically classified at the genus or species level, indicating that a broader spectrum of taxonomic groups is involved in subsoil denitrification than previously recognized. A denitrifier isolate library was constructed by using a strategy combining high-throughput and conventional cultivation techniques. Assessment of the denitrification characteristics of both the MAGs and isolates demonstrated the dominance of truncated denitrification. Functional screening revealed the highest denitrification activity in two complete denitrifiers belonging to the genus Pseudomonas. These findings greatly expand the current knowledge of the composition and function of denitrifying microorganisms in subsoils. The constructed isolate library provided the first pool of subsoil-denitrifying microorganisms that could facilitate the development of microbe-based technologies for nitrate attenuation in groundwater.

Soil Humic Acid Stimulates Potentially Active Dissimilatory Arsenate-Reducing Bacteria in Flooded Paddy Soil as Revealed by Metagenomic Stable Isotope Probing

Dissimilatory arsenate reduction contributes a large proportion of arsenic flux from flooded paddy soil, which is closely linked to soil organic carbon input and efflux. Humic acid (HA) represents a natural ingredient in soil and is shown to enhance microbial arsenate respiration to promote arsenic mobility. However, the community and function profiles of metabolically active arsenate-respiring bacteria and their interactions with HA in paddy soil remain unclear. To probe this linkage, we performed a genome-centric comparison of potentially active arsenate-respiring bacteria in anaerobic microcosms amended with 13C-lactate and HA by combining stable-isotope probing with genome-resolved metagenomics. Indeed, HA greatly accelerated the microbial reduction of arsenate to arsenite. Enrichment of bacteria that harbor arsenate-respiring reductase genes (arrA) in HA-enriched 13C-DNA was confirmed by metagenomic binning, which are affiliated with Firmicutes (mainly Desulfitobacterium, Bacillus, Brevibacillus, and Clostridia) and Acidobacteria. Characterization of reference extracellular electron transfer (EET)-related genes in these arrA-harboring bacteria supports the presence of EET-like genes, with partial electron-transport chain genes identified. This suggests that Gram-positive Firmicutes- and Acidobacteria-related members may harbor unspecified EET-associated genes involved in metal reduction. Our findings highlight the link between soil HA and potentially active arsenate-respiring bacteria, which can be considered when using HA for arsenic removal.

Unusual rearrangements of mitogenomes in Diptera revealed by comparative analysis of 135 tachinid species (Insecta, Diptera, Tachinidae)

The Tachinidae is one of the most speciose families in Diptera, and the exclusively parasitoid species play an important role in regulating populations of many herbivorous insects in ecosystems, including many agricultural pests. To better comprehend the characteristics and evolution of the mitochondrial genome for the Tachinidae, we are adding a massive amount of new molecular data by assembling the mitogenomes for 71 genera and 135 species from all four tachinid subfamilies through next-generation sequencing, and we are presenting the most comprehensive mitogenomic phylogenetic analysis of this family so far. Extensive rearrangements observed in the mitogenome of Admontia podomyia (Exoristinae) are unique for the entire suborder Cyclorrhapha. The rearrangement pattern suggests that the process involved a tandem duplication of the complete mitogenome, followed by both random and nonrandom loss of one copy of each gene. Additionally, five minor mitogenome rearrangements are discovered and described in three subfamilies. We present the largest species-level phylogenetic hypothesis for Tachinidae to date, based on mitogenomes of 152 species of Tachinidae, representing all four subfamilies and with five non-tachinid outgroups. Our analyses support the monophyly of the Tachinidae and most tribes and genera were recovered with good support, but the higher-level phylogenetic relationships within Tachinidae were poorly resolved, indicating that mitogenome data alone are not enough to unambiguously resolve the deeper phylogenetic relationships within Tachinidae.

Gene inversion led to the emergence of brackish archaeal heterotrophs in the aftermath of the Cryogenian Snowball Earth

Land-ocean interactions greatly impact the evolution of coastal life on earth. However, the ancient geological forces and genetic mechanisms that shaped evolutionary adaptations and allowed microorganisms to inhabit coastal brackish waters remain largely unexplored. In this study, we infer the evolutionary trajectory of the ubiquitous heterotrophic archaea Poseidoniales (Marine Group II archaea) presently occurring across global aquatic habitats. Our results show that their brackish subgroups had a single origination, dated to over 600 million years ago, through the inversion of the magnesium transport gene corA that conferred osmotic-stress tolerance. The subsequent loss and gain of corA were followed by genome-wide adjustment, characterized by a general two-step mode of selection in microbial speciation. The coastal family of Poseidoniales showed a rapid increase in the evolutionary rate during and in the aftermath of the Cryogenian Snowball Earth (∼700 million years ago), possibly in response to the enhanced phosphorus supply and the rise of algae. Our study highlights the close interplay between genetic changes and ecosystem evolution that boosted microbial diversification in the Neoproterozoic continental margins, where the Cambrian explosion of animals soon followed.

Diverse hydrocarbon degradation genes, heavy metal resistome, and microbiome of a fluorene-enriched animal-charcoal polluted soil

Environmental compartments polluted with animal charcoal from the skin and hide cottage industries are rich in toxic heavy metals and diverse hydrocarbon classes, some of which are carcinogenic, mutagenic, and genotoxic, and thus require a bio-based eco-benign decommission strategies. A shotgun metagenomic approach was used to decipher the microbiome, hydrocarbon degradation genes, and heavy metal resistome of a microbial consortium (FN8) from an animal-charcoal polluted site enriched with fluorene. Structurally, the FN8 microbial consortium consists of 26 phyla, 53 classes, 119 orders, 245 families, 620 genera, and 1021 species. The dominant phylum, class, order, family, genus, and species in the consortium are Proteobacteria (51.37%), Gammaproteobacteria (39.01%), Bacillales (18.09%), Microbulbiferaceae (11.65%), Microbulbifer (12.21%), and Microbulbifer sp. A4B17 (19.65%), respectively. The microbial consortium degraded 57.56% (28.78 mg/L) and 87.14% (43.57 mg/L) of the initial fluorene concentration in 14 and 21 days. Functional annotation of the protein sequences (ORFs) of the FN8 metagenome using the KEGG GhostKOALA, KofamKOALA, NCBI's conserved domain database, and BacMet revealed the detection of hydrocarbon degradation genes for benzoate, aminobenzoate, polycyclic aromatic hydrocarbons (PAHs), chlorocyclohexane/chlorobenzene, chloroalkane/chloroalkene, toluene, xylene, styrene, naphthalene, nitrotoluene, and several others. The annotation also revealed putative genes for the transport, uptake, efflux, and regulation of heavy metals such as arsenic, cadmium, chromium, mercury, nickel, copper, zinc, and several others. Findings from this study have established that members of the FN8 consortium are well-adapted and imbued with requisite gene sets and could be a potential bioresource for on-site depuration of animal charcoal polluted sites.

Mitochondrial genomes provide insights into the Euholognatha (Insecta: Plecoptera)

BACKGROUND

Euholognatha is a monophyletic group within stoneflies comprised by a superfamily Nemouroidea and a family Scopuridae. Based on morphological data, the family-level phylogenetic relationships within Euholognatha are widely accepted, but there is still controversy among different molecular studies. To better understand the phylogeny of all six extant euholognathan families, we sequenced and analyzed seven euholognathan mitogenomes.

RESULTS

The sequence heterogeneity analysis observed a low degree of compositional heterogeneity in euholognathan mitogenomes. Meanwhile, leuctrid mitogenomes were more heterogeneous than other euholognathan families, which may affect the phylogenetic reconstruction. Phylogenetic analyses with various datasets generated three topologies. The Leuctridae was recovered as the earliest branching lineage, and the sister relationship of Capniidae and Taeniopterygidae was supported by most tree topologies and FcLM analyses. When separately excluding sparsely sampled Scopuridae or high heterogeneity leuctrid taxa, phylogenetic analyses under the same methods generated more stable and consistent tree topologies. Finally, based on the results of this study, we reconstructed the relationships within Euholognatha as: Leuctridae + (Scopuridae + ((Taeniopterygidae + Capniidae) + (Nemouridae + Notonemouridae))).

CONCLUSION

Our research shows the potential of data optimizing strategies in reconstructing phylogeny within Euholognatha and provides new insight into the phylogeny of this group.

Comparative metagenomics and microbial dynamics of jute retting environment

Jute, eco-friendly natural fiber, depends on conventional water-based microbial retting process that suffers from the production of low-quality fiber, restricting its diversified applications. The efficiency of water retting of jute depends on plant polysaccharide fermenting pectinolytic microorganisms. Understanding the phase difference in retting microbial community composition is crucial to provide knowledge on the functions of each member of microbiota for the improvement of retting and fiber quality. The retting microbiota profiling of jute was commonly performed previously using only one retting phase with culture-dependent methods which has limited coverage and accuracy. Here, for the first we have analyzed jute retting water through WGS metagenome approach in three phases (pre-retting, aerobic retting, and anaerobic retting phases) and characterized the microbial communities both culturable and non-culturable along with their dynamics with the fluctuation of oxygen availability. Our analysis revealed a total of 25.99 × 104 unknown proteins (13.75%), 16.18 × 105 annotated proteins (86.08%), and 32.68 × 102 ribosomal RNA (0.17%) in the pre-retting phase, 15.12 × 104 unknown proteins (8.53%), 16.18 × 105 annotated proteins (91.25%), and 38.62 × 102 ribosomal RNA (0.22%) in the aerobic retting phase, and 22.68 × 102 ribosomal RNA and 80.14 × 104 (99.72%) annotated protein in the anaerobic retting phase. Taxonomically, we identified 53 different phylotypes in the retting environment, with Proteobacteria being the dominant taxa comprising over 60% of the population. We have identified 915 genera from Archaea, Viruses, Bacteria, and Eukaryota in the retting habitat, with anaerobic or facultative anaerobic pectinolytic microflora being enriched in the anoxic, nutrient-rich retting niche, such as Aeromonas (7%), Bacteroides (3%), Clostridium (6%), Desulfovibrio (4%), Acinetobacter (4%), Enterobacter (1%), Prevotella (2%), Acidovorax (3%), Bacillus (1%), Burkholderia (1%), Dechloromonas (2%), Caulobacter (1%) and Pseudomonas (7%). We observed an increase in the expression of 30 different KO functional level 3 pathways in the final retting stage compared to the middle and pre-retting stages. The main functional differences among the retting phases were found to be related to nutrient absorption and bacterial colonization. These findings reveal the bacterial groups that are involved in fiber retting different phases and will facilitate to develop future phase-specific microbial consortia for the improvement of jute retting process.

MuDoGeR: Multi-Domain Genome recovery from metagenomes made easy

Several computational frameworks and workflows that recover genomes from prokaryotes, eukaryotes and viruses from metagenomes exist. Yet, it is difficult for scientists with little bioinformatics experience to evaluate quality, annotate genes, dereplicate, assign taxonomy and calculate relative abundance and coverage of genomes belonging to different domains. MuDoGeR is a user-friendly tool tailored for those familiar with Unix command-line environment that makes it easy to recover genomes of prokaryotes, eukaryotes and viruses from metagenomes, either alone or in combination. We tested MuDoGeR using 24 individual-isolated genomes and 574 metagenomes, demonstrating the applicability for a few samples and high throughput. While MuDoGeR can recover eukaryotic viral sequences, its characterization is predominantly skewed towards bacterial and archaeal viruses, reflecting the field's current state. However, acting as a dynamic wrapper, the MuDoGeR is designed to constantly incorporate updates and integrate new tools, ensuring its ongoing relevance in the rapidly evolving field. MuDoGeR is open-source software available at https://github.com/mdsufz/MuDoGeR. Additionally, MuDoGeR is also available as a Singularity container.

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.