Search

BACKGROUND: Despite being among the most abundant biological entities on earth, bacteriophage (phage) remain an understudied component of host-associated systems. One limitation to studying host-associated phage is the lack of consensus on methods for sampling phage communities. Here, we compare paired total metagenomes and viral size fraction metagenomes (viromes) as methods for investigating the dsDNA viral communities associated with the GI tract of two bee species: the European honey bee Apis mellifera and the eastern bumble bee Bombus impatiens. RESULTS: We find that viromes successfully enriched for phage, thereby increasing phage recovery, but only in honey bees. In contrast, for bumble bees, total metagenomes recovered greater phage diversity. Across both bee species, viromes better sampled low occupancy phage, while total metagenomes were biased towards sampling temperate phage. Additionally, many of the phage captured by total metagenomes were absent altogether from viromes. Comparing between bees, we show that phage communities in commercially reared bumble bees are significantly reduced in diversity compared to honey bees, likely reflecting differences in bacterial titer and diversity. In a broader context, these results highlight the complementary nature of total metagenomes and targeted viromes, especially when applied to host-associated environments. CONCLUSIONS: Overall, we suggest that studies interested in assessing total communities of host-associated phage should consider using both approaches. However, given the constraints of virome sampling, total metagenomes may serve to sample phage communities with the understanding that they will preferentially sample dominant and temperate phage. Video Abstract.

UNLABELLED: Gastrointestinal diseases are the most frequently reported clinical problems in captive common marmosets (Callithrix jacchus), often affecting the health and welfare of the animal and ultimately their use as a research subject. The microbiome has been shown to be intimately connected to diet and gastrointestinal health. Here, we use shotgun metagenomics and untargeted metabolomics in fecal samples of common marmosets collected before, during, and after a dietary transition from a biscuit to a gel diet. The overall health of marmosets, measured as weight recovery and reproductive outcome, improved after the diet transition. Moreover, each marmoset pair had significant shifts in the microbiome and metabolome after the diet transition. In general, we saw a decrease in Escherichia coli and Prevotella species and an increase in Bifidobacterium species. Untargeted metabolic profiles indicated that polyamine levels, specifically cadaverine and putrescine, were high after diet transition, suggesting either an increase in excretion or a decrease in intestinal reabsorption at the intestinal level. In conclusion, our data suggest that Bifidobacterium species could potentially be useful as probiotic supplements to the laboratory marmoset diet. Future studies with a larger sample size will be beneficial to show that this is consistent with the diet change. IMPORTANCE: Appropriate diet and health of the common marmoset in captivity are essential both for the welfare of the animal and to improve experimental outcomes. Our study shows that a gel diet compared to a biscuit diet improves the health of a marmoset colony, is linked to increases in Bifidobacterium species, and increases the removal of molecules associated with disease. The diet transition had an influence on the molecular changes at both the pair and time point group levels, but only at the pair level for the microbial changes. It appears to be more important which genes and functions present changed rather than specific microbes. Further studies are needed to identify specific components that should be considered when choosing an appropriate diet and additional supplementary foods, as well as to validate the benefits of providing probiotics. Probiotics containing Bifidobacterium species appear to be useful as probiotic supplements to the laboratory marmoset diet, but additional work is needed to validate these findings.

The gut microbiota (GM) can regulate bone mass, but its association with incident fractures is unknown. We used Cox regression models to determine whether the GM composition is associated with incident fractures in the large FINRISK 2002 cohort (n = 7043, 1092 incident fracture cases, median follow-up time 18 years) with information on GM composition and functionality from shotgun metagenome sequencing. Higher alpha diversity was associated with decreased fracture risk (hazard ratio [HR] 0.92 per standard deviation increase in Shannon index, 95% confidence interval 0.87-0.96). For beta diversity, the first principal component was associated with fracture risk (Aitchison distance, HR 0.90, 0.85-0.96). In predefined phyla analyses, we observed that the relative abundance of Proteobacteria was associated with increased fracture risk (HR 1.14, 1.07-1.20), while the relative abundance of Tenericutes was associated with decreased fracture risk (HR 0.90, 0.85-0.96). Explorative sub-analyses within the Proteobacteria phylum showed that higher relative abundance of Gammaproteobacteria was associated with increased fracture risk. Functionality analyses showed that pathways related to amino acid metabolism and lipopolysaccharide biosynthesis associated with fracture risk. The relative abundance of Proteobacteria correlated with pathways for amino acid metabolism, while the relative abundance of Tenericutes correlated with pathways for butyrate synthesis. In conclusion, the overall GM composition was associated with incident fractures. The relative abundance of Proteobacteria, especially Gammaproteobacteria, was associated with increased fracture risk, while the relative abundance of Tenericutes was associated with decreased fracture risk. Functionality analyses demonstrated that pathways known to regulate bone health may underlie these associations.

BACKGROUND: Chronic exposure to microorganisms inside homes can impact respiratory health. Few studies have used advanced sequencing methods to examine adult respiratory outcomes, especially continuous measures. We aimed to identify metagenomic profiles in house dust related to the quantitative traits of pulmonary function and airway inflammation in adults. Microbial communities, 1264 species (389 genera), in vacuumed bedroom dust from 779 homes in a US cohort were characterized by whole metagenome shotgun sequencing. We examined two overall microbial diversity measures: richness (the number of individual microbial species) and Shannon index (reflecting both richness and relative abundance). To identify specific differentially abundant genera, we applied the Lasso estimator with high-dimensional inference methods, a novel framework for analyzing microbiome data in relation to continuous traits after accounting for all taxa examined together. RESULTS: Pulmonary function measures (forced expiratory volume in one second (FEV1), forced vital capacity (FVC), and FEV1/FVC ratio) were not associated with overall dust microbial diversity. However, many individual microbial genera were differentially abundant (p-value < 0.05 controlling for all other microbial taxa examined) in relation to FEV1, FVC, or FEV1/FVC. Similarly, fractional exhaled nitric oxide (FeNO), a marker of airway inflammation, was unrelated to overall microbial diversity but associated with differential abundance for many individual genera. Several genera, including Limosilactobacillus, were associated with a pulmonary function measure and FeNO, while others, including Moraxella to FEV1/FVC and Stenotrophomonas to FeNO, were associated with a single trait. CONCLUSIONS: Using state-of-the-art metagenomic sequencing, we identified specific microorganisms in indoor dust related to pulmonary function and airway inflammation. Some were previously associated with respiratory conditions; others were novel, suggesting specific environmental microbial components contribute to various respiratory outcomes. The methods used are applicable to studying microbiome in relation to other continuous outcomes. Video Abstract.

Historically neglected by microbial ecologists, soil viruses are now thought to be critical to global biogeochemical cycles. However, our understanding of their global distribution, activities and interactions with the soil microbiome remains limited. Here we present the Global Soil Virus Atlas, a comprehensive dataset compiled from 2,953 previously sequenced soil metagenomes and composed of 616,935 uncultivated viral genomes and 38,508 unique viral operational taxonomic units. Rarefaction curves from the Global Soil Virus Atlas indicate that most soil viral diversity remains unexplored, further underscored by high spatial turnover and low rates of shared viral operational taxonomic units across samples. By examining genes associated with biogeochemical functions, we also demonstrate the viral potential to impact soil carbon and nutrient cycling. This study represents an extensive characterization of soil viral diversity and provides a foundation for developing testable hypotheses regarding the role of the virosphere in the soil microbiome and global biogeochemistry.

IntroductionEmerging evidence links changes in the gut microbiome to late-onset Alzheimer's disease (LOAD), necessitating examination of AD mouse models with consideration of the microbiome.MethodsWe used shotgun metagenomics and untargeted metabolomics to study the human amyloid beta knock-in (hAβ-KI) murine model for LOAD compared to both wild-type (WT) mice and a model for early-onset AD (3xTg-AD).ResultsEighteen-month female (but not male) hAβ-KI microbiomes were distinct from WT microbiomes, with AD genotype accounting for 18% of the variance by permutational multivariate analysis of variance (PERMANOVA). Metabolomic diversity differences were observed in females, however no individual metabolites were differentially abundant. hAβ-KI mice microbiomes were distinguishable from 3xTg-AD animals (81% accuracy by random forest modeling), with separation primarily driven by Romboutsia ilealis and Turicibacter species. Microbiomes were highly cage specific, with cage assignment accounting for more than 40% of the PERMANOVA variance between the groups.DiscussionThese findings highlight a sex-dependent variation in the microbiomes of hAβ-KI mice and underscore the importance of considering the microbiome when designing studies that use murine models for AD.HighlightsMicrobial diversity and the abundance of several species differed in human amyloid beta knock-in (hAβ-KI) females but not males. Correlations to Alzheimer's disease (AD) genotype were stronger for the microbiome than the metabolome. Microbiomes from hAβ-KI mice were distinct from 3xTg-AD mice. Cage effects accounted for most of the variance in the microbiome and metabolome.

Bile acids (BA) are among the most abundant metabolites produced by the gut microbiome. Primary BAs produced in the liver are converted by gut bacterial 7-α-dehydroxylation into secondary BAs, which can differentially regulate host health via signaling based on their varying affinity for BA receptors. Despite the importance of secondary BAs in host health, the regulation of 7-α-dehydroxylation and the role of diet in modulating this process is incompletely defined. Understanding this process could lead to dietary guidelines that beneficially shift BA metabolism. Dietary fiber regulates gut microbial composition and metabolite production. We tested the hypothesis that feeding mice a diet rich in a fermentable dietary fiber, resistant starch (RS), would alter gut bacterial BA metabolism. Male and female wild-type mice were fed a diet supplemented with RS or an isocaloric control diet (IC). Metabolic parameters were similar between groups. RS supplementation increased gut luminal deoxycholic acid (DCA) abundance. However, gut luminal cholic acid (CA) abundance, the substrate for 7-α-dehydroxylation in DCA production, was unaltered by RS. Further, RS supplementation did not change the mRNA expression of hepatic BA producing enzymes or ileal BA transporters. Metagenomic assessment of gut bacterial composition revealed no change in the relative abundance of bacteria known to perform 7-α-dehydroxylation. P. ginsenosidimutans and P. multiformis were positively correlated with gut luminal DCA abundance and increased in response to RS supplementation. These data demonstrate that RS supplementation enriches gut luminal DCA abundance without increasing the relative abundance of bacteria known to perform 7-α-dehydroxylation.

The fields of Metagenomics and Metatranscriptomics involve the examination of complete nucleotide sequences, gene identification, and analysis of potential biological functions within diverse organisms or environmental samples. Despite the vast opportunities for discovery in metagenomics, the sheer volume and complexity of sequence data often present challenges in processing analysis and visualization. This article highlights the critical role of advanced visualization tools in enabling effective exploration, querying, and analysis of these complex datasets. Emphasizing the importance of accessibility, the article categorizes various visualizers based on their intended applications and highlights their utility in empowering bioinformaticians and non-bioinformaticians to interpret and derive insights from meta-omics data effectively.

Intraspecific genomic diversity brings the potential for an unreported and diverse reservoir of cryptic antibiotic resistance genes in pathogens, as cryptic resistance can occur without major mutations and horizontal transmission. Here, we predicted the differences in the types of antibiotics and genes that induce cryptic and latent resistance between micro-diverse Escherichia coli strains. For example, we hypothesize that known resistance genes will be the culprit of latent resistance within clinical strains. We used a modified functional metagenomics method to induce expression in eight E. coli strains. We found a total of 66 individual genes conferring phenotypic resistance to 11 out of 16 antibiotics. A total of 14 known antibiotic resistance genes comprised 21% of total identified genes, whereas the majority (52 genes) were unclassified cryptic resistance genes. Between the eight strains, 1.2% of core orthologous genes were positive (conferred resistance in at least one strain). Sixty-four percent of positive orthologous genes conferred resistance to only one strain, demonstrating high intraspecific variability of latent resistance genes. Cryptic resistance genes comprised most resistance genes among laboratory and clinical strains as well as natural, semisynthetic, and synthetic antibiotics. Known antibiotic resistance genes primarily conferred resistance to multiple antibiotics from varying origins and within multiple strains. Hence, it is uncommon for E. coli to develop cross-cryptic resistance to antibiotics from multiple origins or within multiple strains. We have uncovered prospective and previously unknown resistance genes as well as antibiotics that have the potential to trigger latent antibiotic resistance in E. coli strains from varying origins.IMPORTANCEIntraspecific genomic diversity may be a driving force in the emergence of adaptive antibiotic resistance. Adaptive antibiotic resistance enables sensitive bacterial cells to acquire temporary antibiotic resistance, creating an optimal window for the development of permanent mutational resistance. In this study, we investigate cryptic resistance, an adaptive resistance mechanism, and unveil novel (cryptic) antibiotic resistance genes that confer resistance when amplified within eight E. coli strains derived from clinical and laboratory origins. We identify the potential of cryptic resistance genes to confer cross-resistance to antibiotics from varying origins and within multiple strains. We discern antibiotic characteristics that promote latent resistance in multiple strains, considering intraspecific diversity. This study may help detect novel resistance genes and functional genes that could become responsible for cryptic resistance among diverse strains and antibiotics, thus also identifying potential novel antibiotic targets and mechanisms.

Ticks are increasingly important vectors of human and agricultural diseases. While many studies have focused on tick-borne bacteria, far less is known about tick-associated viruses and their roles in public health or tick physiology. To address this, we investigated patterns of bacterial and viral communities across two field populations of western black-legged ticks (Ixodes pacificus). Through metatranscriptomic analysis of 100 individual ticks, we quantified taxon prevalence, abundance, and co-occurrence with other members of the tick microbiome. In addition to commonly found tick-associated microbes, we assembled 11 novel RNA virus genomes from Rhabdoviridae, Chuviridae, Picornaviridae, Phenuiviridae, Reoviridae, Solemovidiae, Narnaviridae and two highly divergent RNA virus genomes lacking sequence similarity to any known viral families. We experimentally verified the presence of these in I. pacificus ticks across several life stages. We also unexpectedly identified numerous virus-like transcripts that are likely encoded by tick genomic DNA, and which are distinct from known endogenous viral element-mediated immunity pathways in invertebrates. Taken together, our work reveals that I. pacificus ticks carry a greater diversity of viruses than previously appreciated, in some cases resulting in evolutionarily acquired virus-like transcripts. Our findings highlight how pervasive and intimate tick-virus interactions are, with major implications for both the fundamental biology and vectorial capacity of I. pacificus ticks.ImportanceTicks are increasingly important vectors of disease, particularly in the United States where expanding tick ranges and intrusion into previously wild areas has resulted in increasing human exposure to ticks. Emerging human pathogens have been identified in ticks at an increasing rate, and yet little is known about the full community of microbes circulating in various tick species, a crucial first step to understanding how they interact with each and their tick host, as well as their ability to cause disease in humans. We investigated the bacterial and viral communities of the Western blacklegged tick in California and found 11 previously uncharacterized viruses circulating in this population.

Studies using 16S rRNA and shotgun metagenomics typically yield different results, usually attributed to PCR amplification biases. We introduce Greengenes2, a reference tree that unifies genomic and 16S rRNA databases in a consistent, integrated resource. By inserting sequences into a whole-genome phylogeny, we show that 16S rRNA and shotgun metagenomic data generated from the same samples agree in principal coordinates space, taxonomy and phenotype effect size when analyzed with the same tree.

Microbes are essential for the functioning of all ecosystems, and as global warming and anthropogenic pollution threaten ecosystems, it is critical to understand how microbes respond to these changes. We investigated the climate response of Sphingomonas, a widespread gram-negative bacterial genus, during an 18-month microbial community reciprocal transplant experiment across a Southern California climate gradient. We hypothesized that after 18 months, the transplanted Sphingomonas clade and functional composition would correspond with site conditions and reflect the Sphingomonas composition of native communities. We extracted Sphingomonas sequences from metagenomic data across the gradient and assessed their clade and functional composition. Representatives of at least 12 major Sphingomonas clades were found at varying relative abundances along the climate gradient, and transplanted Sphingomonas clade composition shifted after 18 months. Site had a significant effect (PERMANOVA; P < 0.001) on the distribution of both Sphingomonas functional (R2 = 0.465) and clade composition (R2 = 0.400), suggesting that Sphingomonas composition depends on climate parameters. Additionally, for both Sphingomonas clade and functional composition, ordinations revealed that the transplanted communities shifted closer to the native Sphingomonas composition of the grassland site compared with the site they were transplanted into. Overall, our results indicate that climate and substrate collectively determine Sphingomonas clade and functional composition.IMPORTANCESphingomonas is the most abundant gram-negative bacterial genus in litter-degrading microbial communities of desert, grassland, shrubland, and forest ecosystems in Southern California. We aimed to determine whether Sphingomonas responds to climate change in the same way as gram-positive bacteria and whole bacterial communities in these ecosystems. Within Sphingomonas, both clade composition and functional genes shifted in response to climate and litter chemistry, supporting the idea that bacteria respond similarly to climate at different scales of genetic variation. This understanding of how microbes respond to perturbation across scales may aid in future predictions of microbial responses to climate change.

The infant gut microbiome is increasingly recognized as a reservoir of antibiotic resistance genes, yet the assembly of gut resistome in infants and its influencing factors remain largely unknown. We characterized resistome in 4132 metagenomes from 963 infants in six countries and 4285 resistance genes were observed. The inherent resistome pattern of healthy infants (N = 272) could be distinguished by two stages: a multicompound resistance phase (Months 0-7) and a tetracycline-mupirocin-β-lactam-dominant phase (Months 8-14). Microbial taxonomy explained 40.7% of the gut resistome of healthy infants, with Escherichia (25.5%) harboring the most resistance genes. In a further analysis with all available infants (N = 963), we found age was the strongest influencer on the resistome and was negatively correlated with the overall resistance during the first 3 years (p

MOTIVATION: Taxonomic classification of short reads and taxonomic profiling of metagenomic samples are well-studied yet challenging problems. The presence of species belonging to groups without close representation in a reference dataset is particularly challenging. While k-mer-based methods have performed well in terms of running time and accuracy, they tend to have reduced accuracy for such novel species. Thus, there is a growing need for methods that combine the scalability of k-mers with increased sensitivity. RESULTS: Here, we show that using locality-sensitive hashing (LSH) can increase the sensitivity of the k-mer-based search. Our method, which combines LSH with several heuristics techniques including soft lowest common ancestor labeling and voting, is more accurate than alternatives in both taxonomic classification of individual reads and abundance profiling. AVAILABILITY AND IMPLEMENTATION: CONSULT-II is implemented in C++, and the software, together with reference libraries, is publicly available on GitHub https://github.com/bo1929/CONSULT-II.

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.

Survival and growth strategies of Antarctic endolithic microbes residing in Earth's driest and coldest desert remain virtually unknown. From 109 endolithic microbiomes, 4539 metagenome-assembled genomes were generated, 49.3 % of which were novel candidate bacterial species. We present evidence that trace gas oxidation and atmospheric chemosynthesis may be the prevalent strategies supporting metabolic activity and persistence of these ecosystems at the fringe of life and the limits of habitability.

BACKGROUND: Prior to soil formation, phosphate liberated by rock weathering is often sequestered into highly insoluble lanthanide phosphate minerals. Dissolution of these minerals releases phosphate and lanthanides to the biosphere. Currently, the microorganisms involved in phosphate mineral dissolution and the role of lanthanides in microbial metabolism are poorly understood. RESULTS: Although there have been many studies of soil microbiology, very little research has investigated microbiomes of weathered rock. Here, we sampled weathered granite and associated soil to identify the zones of lanthanide phosphate mineral solubilisation and genomically define the organisms implicated in lanthanide utilisation. We reconstructed 136 genomes from 11 bacterial phyla and found that gene clusters implicated in lanthanide-based metabolism of methanol (primarily xoxF3 and xoxF5) are surprisingly common in microbial communities in moderately weathered granite. Notably, xoxF3 systems were found in Verrucomicrobia for the first time, and in Acidobacteria, Gemmatimonadetes and Alphaproteobacteria. The xoxF-containing gene clusters are shared by diverse Acidobacteria and Gemmatimonadetes, and include conserved hypothetical proteins and transporters not associated with the few well studied xoxF systems. Given that siderophore-like molecules that strongly bind lanthanides may be required to solubilise lanthanide phosphates, it is notable that candidate metallophore biosynthesis systems were most prevalent in bacteria in moderately weathered rock, especially in Acidobacteria with lanthanide-based systems. CONCLUSIONS: Phosphate mineral dissolution, putative metallophore production and lanthanide utilisation by enzymes involved in methanol oxidation linked to carbonic acid production co-occur in the zone of moderate granite weathering. In combination, these microbial processes likely accelerate the conversion of granitic rock to soil.

We present six whole community shotgun metagenomic sequencing data sets of two types of biological soil crusts sampled at the ecotone of the Mojave Desert and Colorado Desert in California. These data will help us understand the diversity and function of biocrust microbial communities, which are essential for desert ecosystems.

Contamination of genomes is an increasingly recognized problem affecting several downstream applications, from comparative evolutionary genomics to metagenomics. Here we introduce ContScout, a precise tool for eliminating foreign sequences from annotated genomes. It achieves high specificity and sensitivity on synthetic benchmark data even when the contaminant is a closely related species, outperforms competing tools, and can distinguish horizontal gene transfer from contamination. A screen of 844 eukaryotic genomes for contamination identified bacteria as the most common source, followed by fungi and plants. Furthermore, we show that contaminants in ancestral genome reconstructions lead to erroneous early origins of genes and inflate gene loss rates, leading to a false notion of complex ancestral genomes. Taken together, we offer here a tool for sensitive removal of foreign proteins, identify and remove contaminants from diverse eukaryotic genomes and evaluate their impact on phylogenomic analyses.

Multi-omics approaches have been successfully applied to investigate pregnancy and health outcomes at a molecular and genetic level in several studies. As omics technologies advance, research areas are open to study further. Here we discuss overall trends and examples of successfully using omics technologies and techniques (e.g., genomics, proteomics, metabolomics, and metagenomics) to investigate the molecular epidemiology of pregnancy. In addition, we outline omics applications and study characteristics of pregnancy for understanding fundamental biology, causal health, and physiological relationships, risk and prediction modeling, diagnostics, and correlations.

Microbial communities have evolved to colonize all ecosystems of the planet, from the deep sea to the human gut. Microbes survive by sensing, responding, and adapting to immediate environmental cues. This process is driven by signal transduction proteins such as histidine kinases, which use their sensing domains to bind or otherwise detect environmental cues and transduce signals to adjust internal processes. We hypothesized that an ecosystems unique stimuli leave a sensor fingerprint, able to identify and shed insight on ecosystem conditions. To test this, we collected 20,712 publicly available metagenomes from Host-associated, Environmental, and Engineered ecosystems across the globe. We extracted and clustered the collections nearly 18M unique sensory domains into 113,712 similar groupings with MMseqs2. We built gradient-boosted decision tree machine learning models and found we could classify the ecosystem type (accuracy: 87%) and predict the levels of different physical parameters (R2 score: 83%) using the sensor cluster abundance as features. Feature importance enables identification of the most predictive sensors to differentiate between ecosystems which can lead to mechanistic interpretations if the sensor domains are well annotated. To demonstrate this, a machine learning model was trained to predict patients disease state and used to identify domains related to oxygen sensing present in a healthy gut but missing in patients with abnormal conditions. Moreover, since 98.7% of identified sensor domains are uncharacterized, importance ranking can be used to prioritize sensors to determine what ecosystem function they may be sensing. Furthermore, these new predictive sensors can function as targets for novel sensor engineering with applications in biotechnology, ecosystem maintenance, and medicine.IMPORTANCEMicrobes infect, colonize, and proliferate due to their ability to sense and respond quickly to their surroundings. In this research, we extract the sensory proteins from a diverse range of environmental, engineered, and host-associated metagenomes. We trained machine learning classifiers using sensors as features such that it is possible to predict the ecosystem for a metagenome from its sensor profile. We use the optimized models feature importance to identify the most impactful and predictive sensors in different environments. We next use the sensor profile from human gut metagenomes to classify their disease states and explore which sensors can explain differences between diseases. The sensors most predictive of environmental labels here, most of which correspond to uncharacterized proteins, are a useful starting point for the discovery of important environment signals and the development of possible diagnostic interventions.

Western Eurasia witnessed several large-scale human migrations during the Holocene1-5. Here, to investigate the cross-continental effects of these migrations, we shotgun-sequenced 317 genomes-mainly from the Mesolithic and Neolithic periods-from across northern and western Eurasia. These were imputed alongside published data to obtain diploid genotypes from more than 1,600 ancient humans. Our analyses revealed a 'great divide' genomic boundary extending from the Black Sea to the Baltic. Mesolithic hunter-gatherers were highly genetically differentiated east and west of this zone, and the effect of the neolithization was equally disparate. Large-scale ancestry shifts occurred in the west as farming was introduced, including near-total replacement of hunter-gatherers in many areas, whereas no substantial ancestry shifts happened east of the zone during the same period. Similarly, relatedness decreased in the west from the Neolithic transition onwards, whereas, east of the Urals, relatedness remained high until around 4,000 BP, consistent with the persistence of localized groups of hunter-gatherers. The boundary dissolved when Yamnaya-related ancestry spread across western Eurasia around 5,000 BP, resulting in a second major turnover that reached most parts of Europe within a 1,000-year span. The genetic origin and fate of the Yamnaya have remained elusive, but we show that hunter-gatherers from the Middle Don region contributed ancestry to them. Yamnaya groups later admixed with individuals associated with the Globular Amphora culture before expanding into Europe. Similar turnovers occurred in western Siberia, where we report new genomic data from a 'Neolithic steppe' cline spanning the Siberian forest steppe to Lake Baikal. These prehistoric migrations had profound and lasting effects on the genetic diversity of Eurasian populations.

Functional traits influence the assembly of microbial communities, but identifying these traits in the environment has remained challenging. We studied ectomycorrhizal fungal (EMF) communities inhabiting Populus trichocarpa roots distributed across a precipitation gradient in the Pacific Northwest, USA. We profiled these communities using taxonomic (meta-barcoding) and functional (metagenomic) approaches. We hypothesized that genes involved in fungal drought-stress tolerance and fungal mediated plant water uptake would be most abundant in drier soils. We were unable to detect support for this hypothesis; instead, the abundance of genes involved in melanin synthesis, hydrophobins, aquaporins, trehalose-synthases, and other gene families exhibited no significant shifts across the gradient. Finally, we studied variation in sequence homology for certain genes, finding that fungal communities in dry soils are composed of distinct aquaporin and hydrophobin gene sequences. Altogether, our results suggest that while EMF communities exhibit significant compositional shifts across this gradient, coupled functional turnover, at least as inferred using community metagenomics is limited. Accordingly, the consequences of these distinct EMF communities on plant water uptake remain critically unknown, and future studies targeting the expression of genes involved in drought stress tolerance are required.

Although long-read sequencing has enabled obtaining high-quality and complete genomes from metagenomes, many challenges still remain to completely decompose a metagenome into its constituent prokaryotic and viral genomes. This study focuses on decomposing an estuarine metagenome to obtain a more accurate estimate of microbial diversity. To achieve this, we developed a new bead-based DNA extraction method, a novel bin refinement method, and obtained 150 Gbp of Nanopore sequencing. We estimate that there are ~500 bacterial and archaeal species in our sample and obtained 68 high-quality bins (>90% complete, 100 kbp, and all ribosomal and tRNA genes). We also obtained many contigs of picoeukaryotes, environmental DNA of larger eukaryotes such as mammals, and complete mitochondrial and chloroplast genomes and detected ~40,000 viral populations. Our analysis indicates that there are only a few strains that comprise most of the species abundances. Ocean and estuarine microbiomes play critical roles in global element cycling and ecosystem function. Despite the importance of these microbial communities, many species still have not been cultured in the lab. Environmental sequencing is the primary way the function and population dynamics of these communities can be studied. Long-read sequencing provides an avenue to overcome limitations of short-read technologies to obtain complete microbial genomes but comes with its own technical challenges, such as needed sequencing depth and obtaining high-quality DNA. We present here new sampling and bioinformatics methods to attempt decomposing an estuarine microbiome into its constituent genomes. Our results suggest there are only a few strains that comprise most of the species abundances from viruses to picoeukaryotes, and to fully decompose a metagenome of this diversity requires 1 Tbp of long-read sequencing. We anticipate that as long-read sequencing technologies continue to improve, less sequencing will be needed.

Comparative analysis of (meta)genomes necessitates aggregation, integration, and synthesis of well-annotated data using standards. The Genomic Standards Consortium (GSC) collaborates with the research community to develop and maintain the Minimum Information about any (x) Sequence (MIxS) reporting standard for genomic data. To facilitate the use of the GSC's MIxS reporting standard, we provide a description of the structure and terminology, how to navigate ontologies for required terms in MIxS, and demonstrate practical usage through a soil metagenome example.

Culture-based microbial natural product discovery strategies fail to realize the extraordinary biosynthetic potential detected across earth's microbiomes. Here we introduce Small Molecule In situ Resin Capture (SMIRC), a culture-independent method to obtain natural products directly from the environments in which they are produced. We use SMIRC to capture numerous compounds including two new carbon skeletons that were characterized using NMR and contain structural features that are, to the best of our knowledge, unprecedented among natural products. Applications across diverse marine habitats reveal biome-specific metabolomic signatures and levels of chemical diversity in concordance with sequence-based predictions. Expanded deployments, in situ cultivation, and metagenomics facilitate compound discovery, enhance yields, and link compounds to candidate producing organisms, although microbial community complexity creates challenges for the later. This compound-first approach to natural product discovery provides access to poorly explored chemical space and has implications for drug discovery and the detection of chemically mediated biotic interactions.

Giant viruses (GVs) are key players in ecosystem functioning, biogeochemistry, and eukaryotic genome evolution. GV diversity and abundance in aquatic systems can exceed that of prokaryotes, but their diversity and ecology in lakes, especially polar ones, remain poorly understood. We conducted a comprehensive survey and meta-analysis of GV diversity across 20 lakes, spanning polar to temperate regions, combining our extensive lake metagenome database from the Canadian Arctic and subarctic with publicly available datasets. Leveraging a novel GV genome identification tool, we identified 3304 GV metagenome-assembled genomes, revealing lakes as untapped GV reservoirs. Phylogenomic analysis highlighted their dispersion across all Nucleocytoviricota orders. Strong GV population endemism emerged between lakes from similar regions and biomes (Antarctic and Arctic), but a polar/temperate barrier in lacustrine GV populations and differences in their gene content could be observed. Our study establishes a robust genomic reference for future investigations into lacustrine GV ecology in fast changing polar environments.

Python for Population Genomics (PyPop) is a software package that processes genotype and allele data and performs large-scale population genetic analyses on highly polymorphic multi-locus genotype data. In particular, PyPop tests data conformity to Hardy-Weinberg equilibrium expectations, performs Ewens-Watterson tests for selection, estimates haplotype frequencies, measures linkage disequilibrium, and tests significance. Standardized means of performing these tests is key for contemporary studies of evolutionary biology and population genetics, and these tests are central to genetic studies of disease association as well. Here, we present PyPop 1.0.0, a new major release of the package, which implements new features using the more robust infrastructure of GitHub, and is distributed via the industry-standard Python Package Index. New features include implementation of the asymmetric linkage disequilibrium measures and, of particular interest to the immunogenetics research communities, support for modern nomenclature, including colon-delimited allele names, and improvements to meta-analysis features for aggregating outputs for multiple populations. Code available at: https://zenodo.org/records/10080668 and https://github.com/alexlancaster/pypop.

Cyanobacteria form diverse communities and are important primary producers in Antarctic freshwater environments, but their geographic distribution patterns in Antarctica and globally are still unresolved. There are however few genomes of cultured cyanobacteria from Antarctica available and therefore metagenome-assembled genomes (MAGs) from Antarctic cyanobacteria microbial mats provide an opportunity to explore distribution of uncultured taxa. These MAGs also allow comparison with metagenomes of cyanobacteria enriched communities from a range of habitats, geographic locations, and climates. However, most MAGs do not contain 16S rRNA gene sequences, making a 16S rRNA gene-based biogeography comparison difficult. An alternative technique is to use large-scale k-mer searching to find genomes of interest in public metagenomes. This paper presents the results of k-mer based searches for 5 Antarctic cyanobacteria MAGs from Lake Fryxell and Lake Vanda, assigned the names Phormidium pseudopriestleyi FRX01, Microcoleus sp. MP8IB2.171, Leptolyngbya sp. BulkMat.35, Pseudanabaenaceae cyanobacterium MP8IB2.15, and Leptolyngbyaceae cyanobacterium MP9P1.79 in 498,942 unassembled metagenomes from the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA). The Microcoleus sp. MP8IB2.171 MAG was found in a wide variety of environments, the P. pseudopriestleyi MAG was found in environments with challenging conditions, the Leptolyngbyaceae cyanobacterium MP9P1.79 MAG was only found in Antarctica, and the Leptolyngbya sp. BulkMat.35 and Pseudanabaenaceae cyanobacterium MP8IB2.15 MAGs were found in Antarctic and other cold environments. The findings based on metagenome matches and global comparisons suggest that these Antarctic cyanobacteria have distinct distribution patterns ranging from locally restricted to global distribution across the cold biosphere and other climatic zones.

Exploring the natural diversity of functional genes/proteins from environmental DNA in high throughput remains challenging. In this study, we developed a sequence-based functional metagenomics procedure for mining the diversity of copper (Cu) resistance gene copA in global microbiomes, by combining the metagenomic assembly technology, local BLAST, evolutionary trace analysis (ETA), chemical synthesis, and conventional functional genomics. In total, 87 metagenomes were collected from a public database and subjected to copA detection, resulting in 93,899 hits. Manual curation of 1214 hits of high confidence led to the retrieval of 517 unique CopA candidates, which were further subjected to ETA. Eventually, 175 novel copA sequences of high quality were discovered. Phylogenetic analysis showed that almost all these putative CopA proteins were distantly related to known CopA proteins, with 55 sequences from totally unknown species. Ten novel and three known copA genes were chemically synthesized for further functional genomic tests using the Cu-sensitive Escherichia coli (ΔcopA). The growth test and Cu uptake determination showed that five novel clones had positive effects on host Cu resistance and uptake. One recombinant harboring copA-like 15 (copAL15) successfully restored Cu resistance of the host with a substantially enhanced Cu uptake. Two novel copA genes were fused with the gfp gene and expressed in E. coli for microscopic observation. Imaging results showed that they were successfully expressed and their proteins were localized to the membrane. The results here greatly expand the diversity of known CopA proteins, and the sequence-based procedure developed overcomes biases in length, screening methods, and abundance of conventional functional metagenomics.

Coral reefs are both exceptionally biodiverse and threatened by climate change and other human activities. Here, we review population genomic processes in coral reef taxa and their importance for understanding responses to global change. Many taxa on coral reefs are characterized by weak genetic drift, extensive gene flow, and strong selection from complex biotic and abiotic environments, which together present a fascinating test of microevolutionary theory. Selection, gene flow, and hybridization have played and will continue to play an important role in the adaptation or extinction of coral reef taxa in the face of rapid environmental change, but research remains exceptionally limited compared to the urgent needs. Critical areas for future investigation include understanding evolutionary potential and the mechanisms of local adaptation, developing historical baselines, and building greater research capacity in the countries where most reef diversity is concentrated.

Acarospora socialis, the bright cobblestone lichen, is commonly found in southwestern North America. This charismatic yellow lichen is a species of key ecological significance as it is often a pioneer species in new environments. Despite their ecological importance virtually no research has been conducted on the genomics of A. socialis. To address this, we used long-read sequencing to generate the first high-quality draft genome of A. socialis. Lichen thallus tissue was collected from Pinkham Canyon in Joshua Tree National Park, California and deposited in the UC Riverside herbarium under accession #295874. The de novo assembly of the mycobiont partner of the lichen was generated from Pacific Biosciences HiFi long reads and Dovetail Omni-C chromatin capture data. After removing algal and bacterial contigs, the fungal genome was approximately 31.2 Mb consisting of 38 scaffolds with contig and scaffold N50 of 2.4 Mb. The BUSCO completeness score of the assembled genome was 97.5% using the Ascomycota gene set. Information on the genome of A. socialis is important for California conservation purposes given that this lichen is threatened in some places locally by wildfires due to climate change. This reference genome will be used for understanding the genetic diversity, population genomics, and comparative genomics of A. socialis species. Genomic resources for this species will support population and landscape genomics investigations, exploring the use of A. socialis as a bioindicator species for climate change, and in studies of adaptation by comparing populations that occur across aridity gradients in California.

Metagenomes encode an enormous diversity of proteins, reflecting a multiplicity of functions and activities1,2. Exploration of this vast sequence space has been limited to a comparative analysis against reference microbial genomes and protein families derived from those genomes. Here, to examine the scale of yet untapped functional diversity beyond what is currently possible through the lens of reference genomes, we develop a computational approach to generate reference-free protein families from the sequence space in metagenomes. We analyse 26,931 metagenomes and identify 1.17 billion protein sequences longer than 35 amino acids with no similarity to any sequences from 102,491 reference genomes or the Pfam database3. Using massively parallel graph-based clustering, we group these proteins into 106,198 novel sequence clusters with more than 100 members, doubling the number of protein families obtained from the reference genomes clustered using the same approach. We annotate these families on the basis of their taxonomic, habitat, geographical and gene neighbourhood distributions and, where sufficient sequence diversity is available, predict protein three-dimensional models, revealing novel structures. Overall, our results uncover an enormously diverse functional space, highlighting the importance of further exploring the microbial functional dark matter.

As coral reef ecosystems experience unprecedented change, effective monitoring of reef features supports management, conservation, and intervention efforts. Omic techniques show promise in quantifying key components of reef ecosystems including dissolved metabolites and microorganisms that may serve as invisible sensors for reef ecosystem dynamics. Dissolved metabolites are released by reef organisms and transferred among microorganisms, acting as chemical currencies and contributing to nutrient cycling and signaling on reefs. Here, we applied four omic techniques (taxonomic microbiome via amplicon sequencing, functional microbiome via shotgun metagenomics, targeted metabolomics, and untargeted metabolomics) to waters overlying Floridas Coral Reef, as well as microbiome profiling on individual coral colonies from these reefs to understand how microbes and dissolved metabolites reflect biogeographical, benthic, and nutrient properties of this 500-km barrier reef. We show that the microbial and metabolite omic approaches each differentiated reef habitats based on geographic zone. Further, seawater microbiome profiling and targeted metabolomics were significantly related to more reef habitat characteristics, such as amount of hard and soft coral, compared to metagenomic sequencing and untargeted metabolomics. Across five coral species, microbiomes were also significantly related to reef zone, followed by species and disease status, suggesting that the geographic water circulation patterns in Florida also impact the microbiomes of reef builders. A combination of differential abundance and indicator species analyses revealed metabolite and microbial signatures of specific reef zones, which demonstrates the utility of these techniques to provide new insights into reef microbial and metabolite features that reflect broader ecosystem processes.

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

Next-generation sequencing technologies have enabled many advances across diverse areas of biology, with many benefiting from increased sample size. Although the cost of running next-generation sequencing instruments has dropped substantially over time, the cost of sample preparation methods has lagged behind. To counter this, researchers have adapted library miniaturization protocols and large sample pools to maximize the number of samples that can be prepared by a certain amount of reagents and sequenced in a single run. However, due to high variability of sample quality, over and underrepresentation of samples in a sequencing run has become a major issue in high-throughput sequencing. This leads to misinterpretation of results due to increased noise, and additional time and cost rerunning underrepresented samples. To overcome this problem, we present a normalization method that uses shallow iSeq sequencing to accurately inform pooling volumes based on read distribution. This method is superior to the widely used fluorometry methods, which cannot specifically target adapter-ligated molecules that contribute to sequencing output. Our normalization method not only quantifies adapter-ligated molecules but also allows normalization of feature space; for example, we can normalize to reads of interest such as non-ribosomal reads. As a result, this normalization method improves the efficiency of high-throughput next-generation sequencing by reducing noise and producing higher average reads per sample with more even sequencing depth. IMPORTANCE High-throughput next generation sequencing (NGS) has significantly contributed to the field of genomics; however, further improvements can maximize the potential of this important tool. Uneven sequencing of samples in a multiplexed run is a common issue that leads to unexpected extra costs or low-quality data. To mitigate this problem, we introduce a normalization method based on read counts rather than library concentration. This method allows for an even distribution of features of interest across samples, improving the statistical power of data sets and preventing the financial loss associated with resequencing libraries. This method optimizes NGS, which already has huge importance across many areas of biology.

Stable isotope probing (SIP) facilitates culture-independent identification of active microbial populations within complex ecosystems through isotopic enrichment of nucleic acids. Many DNA-SIP studies rely on 16S rRNA gene sequences to identify active taxa, but connecting these sequences to specific bacterial genomes is often challenging. Here, we describe a standardized laboratory and analysis framework to quantify isotopic enrichment on a per-genome basis using shotgun metagenomics instead of 16S rRNA gene sequencing. To develop this framework, we explored various sample processing and analysis approaches using a designed microbiome where the identity of labeled genomes and their level of isotopic enrichment were experimentally controlled. With this ground truth dataset, we empirically assessed the accuracy of different analytical models for identifying active taxa and examined how sequencing depth impacts the detection of isotopically labeled genomes. We also demonstrate that using synthetic DNA internal standards to measure absolute genome abundances in SIP density fractions improves estimates of isotopic enrichment. In addition, our study illustrates the utility of internal standards to reveal anomalies in sample handling that could negatively impact SIP metagenomic analyses if left undetected. Finally, we present SIPmg, an R package to facilitate the estimation of absolute abundances and perform statistical analyses for identifying labeled genomes within SIP metagenomic data. This experimentally validated analysis framework strengthens the foundation of DNA-SIP metagenomics as a tool for accurately measuring the in situ activity of environmental microbial populations and assessing their genomic potential. IMPORTANCE Answering the questions, "who is eating what?" and "who is active?" within complex microbial communities is paramount for our ability to model, predict, and modulate microbiomes for improved human and planetary health. These questions can be pursued using stable isotope probing to track the incorporation of labeled compounds into cellular DNA during microbial growth. However, with traditional stable isotope methods, it is challenging to establish links between an active microorganism's taxonomic identity and genome composition while providing quantitative estimates of the microorganism's isotope incorporation rate. Here, we report an experimental and analytical workflow that lays the foundation for improved detection of metabolically active microorganisms and better quantitative estimates of genome-resolved isotope incorporation, which can be used to further refine ecosystem-scale models for carbon and nutrient fluxes within microbiomes.

Petabases of environmental metagenomic data are publicly available, presenting an opportunity to characterize complex environments and discover novel lineages of life. Metagenome coassembly, in which many metagenomic samples from an environment are simultaneously analyzed to infer the underlying genomes' sequences, is an essential tool for achieving this goal. We applied MetaHipMer2, a distributed metagenome assembler that runs on supercomputing clusters, to coassemble 3.4 terabases (Tbp) of metagenome data from a tropical soil in the Luquillo Experimental Forest (LEF), Puerto Rico. The resulting coassembly yielded 39 high-quality (>90% complete,

Viruses are the most ubiquitous biological entities on Earth. Even so, elucidating the impact of viruses on microbial communities and associated ecosystem processes often requires identification of unambiguous host-virus linkages-an undeniable challenge in many ecosystems. Subsurface fractured shales present a unique opportunity to first make these strong linkages via spacers in CRISPR-Cas arrays and subsequently reveal complex long-term host-virus dynamics. Here, we sampled two replicated sets of fractured shale wells for nearly 800 days, resulting in 78 metagenomes from temporal sampling of six wells in the Denver-Julesburg Basin (Colorado, USA). At the community level, there was strong evidence for CRISPR-Cas defense systems being used through time and likely in response to viral interactions. Within our host genomes, represented by 202 unique MAGs, we also saw that CRISPR-Cas systems were widely encoded. Together, spacers from host CRISPR loci facilitated 2,110 CRISPR-based viral linkages across 90 host MAGs spanning 25 phyla. We observed less redundancy in host-viral linkages and fewer spacers associated with hosts from the older, more established wells, possibly reflecting enrichment of more beneficial spacers through time. Leveraging temporal patterns of host-virus linkages across differing well ages, we report how host-virus co-existence dynamics develop and converge through time, possibly reflecting selection for viruses that can evade host CRISPR-Cas systems. Together, our findings shed light on the complexities of host-virus interactions as well as long-term dynamics of CRISPR-Cas defense among diverse microbial populations.

Orthonairovirus is a genus of viruses in the family Nairoviridae, order Bunyavirales, with a segmented circular RNA genome. They typically infect birds and mammals and are primarily transmitted by ectoparasites such as ticks. Four of nine Orthonairovirus genogroups can infect humans, with Crimean-Congo hemorrhagic fever virus infections displaying case fatality rates up to 40%. Here, we discover and describe a novel Orthonairovirus as Cencurut virus (CENV). CENV was detected in 34 of 37 Asian house shrews (Suncus murinus) sampled in Singapore and in a nymphal Amblyomma helvolum tick collected from an infected shrew. Pairwise comparison of CENV S, M, and L segments had 95.0 to 100% nucleotide and 97.5 to 100% amino acid homology within CENV genomes, suggesting a diverse viral population. Phylogenetic analysis of the individual gene segments showed that CENV is related to Erve, Lamgora, Lamusara, and Thiafora viruses, with only 49.0 to 58.2% nucleotide and 41.7 to 61.1% amino acid homology, which has previously been detected in other shrew species from France, Gabon, and Senegal respectively. The high detection frequency suggests that CENV is endemic among S. murinus populations in Singapore. The discovery of CENV, from a virus family with known zoonotic potential, underlines the importance of surveillance of synanthropic small mammals that are widely distributed across Southeast Asia.

RNA viruses are ubiquitous components of the global virosphere, yet relatively little is known about their genetic diversity or the cellular mechanisms by which they exploit the biology of their diverse eukaryotic hosts. A hallmark of (+)ssRNA (positive single-stranded RNA) viruses is the ability to remodel host endomembranes for their own replication. However, the subcellular interplay between RNA viruses and host organelles that harbor gene expression systems, such as mitochondria, is complex and poorly understood. Here we report the discovery of 763 new virus sequences belonging to the family Mitoviridae by metatranscriptomic analysis, the identification of previously uncharacterized mitovirus clades, and a putative new viral class. With this expanded understanding of the diversity of mitovirus and encoded RNA-dependent RNA polymerases (RdRps), we annotate mitovirus-specific protein motifs and identify hallmarks of mitochondrial translation, including mitochondrion-specific codons. This study expands the known diversity of mitochondrial viruses and provides additional evidence that they co-opt mitochondrial biology for their survival. IMPORTANCE Metatranscriptomic studies have rapidly expanded the cadre of known RNA viruses, yet our understanding of how these viruses navigate the cytoplasmic milieu of their hosts to survive remains poorly characterized. In this study, we identify and assemble 763 new viral sequences belonging to the Mitoviridae, a family of (+)ssRNA viruses thought to interact with and remodel host mitochondria. We exploit this genetic diversity to identify new clades of Mitoviridae, annotate clade-specific sequence motifs that distinguish the mitoviral RdRp, and reveal patterns of RdRp codon usage consistent with translation on host cell mitoribosomes. These results serve as a foundation for understanding how mitoviruses co-opt mitochondrial biology for their proliferation.

Deforestation threatens the integrity of the Amazon biome and the ecosystem services it provides, including greenhouse gas mitigation. Forest-to-pasture conversion has been shown to alter the flux of methane gas (CH4 ) in Amazonian soils, driving a switch from acting as a sink to a source of atmospheric CH4 . This study aimed to better understand this phenomenon by investigating soil microbial metagenomes, focusing on the taxonomic and functional structure of methane-cycling communities. Metagenomic data from forest and pasture soils were combined with measurements of in situ CH4 fluxes and soil edaphic factors and analysed using multivariate statistical approaches. We found a significantly higher abundance and diversity of methanogens in pasture soils. As inferred by co-occurrence networks, these microorganisms seem to be less interconnected within the soil microbiota in pasture soils. Metabolic traits were also different between land uses, with increased hydrogenotrophic and methylotrophic pathways of methanogenesis in pasture soils. Land-use change also induced shifts in taxonomic and functional traits of methanotrophs, with bacteria harbouring genes encoding the soluble form of methane monooxygenase enzyme (sMMO) depleted in pasture soils. Redundancy analysis and multimodel inference revealed that the shift in methane-cycling communities was associated with high pH, organic matter, soil porosity and micronutrients in pasture soils. These results comprehensively characterize the effect of forest-to-pasture conversion on the microbial communities driving the methane-cycling microorganisms in the Amazon rainforest, which will contribute to the efforts to preserve this important biome.

BackgroundAtrial fibrillation (AF) is an important heart rhythm disorder in aging populations. The gut microbiome composition has been previously related to cardiovascular disease risk factors. Whether the gut microbial profile is also associated with the risk of AF remains unknown.MethodsWe examined the associations of prevalent and incident AF with gut microbiota in the FINRISK 2002 study, a random population sample of 6763 individuals. We replicated our findings in an independent case-control cohort of 138 individuals in Hamburg, Germany.FindingsMultivariable-adjusted regression models revealed that prevalent AF (N = 116) was associated with nine microbial genera. Incident AF (N = 539) over a median follow-up of 15 years was associated with eight microbial genera with false discovery rate (FDR)-corrected P

BackgroundCommon seadragons (Phyllopteryx taeniolatus, Syngnathidae) are an emblem of the diverse endemic fauna of Australia's southern rocky reefs, the newly recognized "Great Southern Reef." A lack of assessments spanning this global biodiversity hotspot in its entirety is currently hampering an understanding of the factors that have contributed to its diversity. The common seadragon has a wide range across Australia's entire temperate south and includes a geogenetic break over a former land bridge, which has called its status as a single species into question. As a popular aquarium display that sells for high prices, common seadragons are also vulnerable to illegal capture.ResultsHere, we provide range-wide nuclear sequences (986 variable Ultraconserved Elements) for 198 individuals and mitochondrial genomes for 140 individuals to assess species status, identify genetic units and their diversity, and trace the source of two poached individuals. Using published data of the other two seadragon species, we found that lineages of common seadragons have diverged relatively recently (

High-throughput sequencing (HTS), more specifically RNA sequencing of plant tissues, has become an indispensable tool for plant virologists to detect and identify plant viruses. During the data analysis step, plant virologists typically compare the obtained sequences to reference virus databases. In this way, they are neglecting sequences without homologies to viruses, which usually represent the majority of sequencing reads. We hypothesized that traces of other pathogens might be detected in this unused sequence data. In the present study, our goal was to investigate whether total RNA-seq data, as generated for plant virus detection, is also suitable for the detection of other plant pathogens and pests. As proof of concept, we first analyzed RNA-seq datasets of plant materials with confirmed infections by cellular pathogens in order to check whether these non-viral pathogens could be easily detected in the data. Next, we set up a community effort to re-analyze existing Illumina RNA-seq datasets used for virus detection to check for the potential presence of non-viral pathogens or pests. In total, 101 datasets from 15 participants derived from 51 different plant species were re-analyzed, of which 37 were selected for subsequent in-depth analyses. In 29 of the 37 selected samples (78%), we found convincing traces of non-viral plant pathogens or pests. The organisms most frequently detected in this way were fungi (15/37 datasets), followed by insects (13/37) and mites (9/37). The presence of some of the detected pathogens was confirmed by independent (q)PCRs analyses. After communicating the results, 6 out of the 15 participants indicated that they were unaware of the possible presence of these pathogens in their sample(s). All participants indicated that they would broaden the scope of their bioinformatic analyses in future studies and thus check for the presence of non-viral pathogens. In conclusion, we show that it is possible to detect non-viral pathogens or pests from total RNA-seq datasets, in this case primarily fungi, insects, and mites. With this study, we hope to raise awareness among plant virologists that their data might be useful for fellow plant pathologists in other disciplines (mycology, entomology, bacteriology) as well.

With blood stream infections (BSIs) representing a major cause of mortality and morbidity worldwide, blood cultures play a crucial role in diagnosis, but their clinical application is dampened by the long turn-around time and the detection of only culturable pathogens. In this study, we developed and validated a shotgun metagenomics next-generation sequencing (mNGS) test directly from positive blood culture fluid, allowing for the identification of fastidious or slow growing microorganisms more rapidly. The test was built based on previously validated next-generation sequencing tests, which rely on several key marker genes for bacterial and fungal identification. The new test utilizes an open-source metagenomics CZ-ID platform for the initial analysis to generate the most likely candidate species, which is then used as a reference genome for downstream, confirmatory analysis. This approach is innovative because it takes advantage of an open-source software's agnostic taxonomic calling capability while still relying on the more established and previously validated marker gene-based identification scheme, increasing the confidence in the final results. The test showed high accuracy (100%, 30/30) for both bacterial and fungal microorganisms. We further demonstrated its clinical utility especially for anaerobes and mycobacteria that are either fastidious, slow growing, or unusual. Although applicable in only limited settings, the Positive Blood Culture mNGS test provides an incremental improvement in solving the unmet clinical needs for the diagnosis of challenging BSIs.

The human gut microbiome harbors substantial ecological diversity at the species level as well as at the strain level within species. In healthy hosts, species abundance fluctuations in the microbiome are thought to be stable, and these fluctuations can be described by macroecological laws. However, it is less clear how strain abundances change over time. An open question is whether individual strains behave like species themselves, exhibiting stability and following the macroecological relationships known to hold at the species level, or whether strains have different dynamics, perhaps due to the relatively close phylogenetic relatedness of cocolonizing lineages. Here, we analyze the daily dynamics of intraspecific genetic variation in the gut microbiomes of four healthy, densely longitudinally sampled hosts. First, we find that the overall genetic diversity of a large majority of species is stationary over time despite short-term fluctuations. Next, we show that fluctuations in abundances in approximately 80% of strains analyzed can be predicted with a stochastic logistic model (SLM), an ecological model of a population experiencing environmental fluctuations around a fixed carrying capacity, which has previously been shown to capture statistical properties of species abundance fluctuations. The success of this model indicates that strain abundances typically fluctuate around a fixed carrying capacity, suggesting that most strains are dynamically stable. Finally, we find that the strain abundances follow several empirical macroecological laws known to hold at the species level. Together, our results suggest that macroecological properties of the human gut microbiome, including its stability, emerge at the level of strains. IMPORTANCE To date, there has been an intense focus on the ecological dynamics of the human gut microbiome at the species level. However, there is considerable genetic diversity within species at the strain level, and these intraspecific differences can have important phenotypic effects on the host, impacting the ability to digest certain foods and metabolize drugs. Thus, to fully understand how the gut microbiome operates in times of health and sickness, its ecological dynamics may need to be quantified at the level of strains. Here, we show that a large majority of strains maintain stable abundances for periods of months to years, exhibiting fluctuations in abundance that can be well described by macroecological laws known to hold at the species level, while a smaller percentage of strains undergo rapid, directional changes in abundance. Overall, our work indicates that strains are an important unit of ecological organization in the human gut microbiome.

Bacteria catalyze the formation and destruction of soil organic matter, but the bacterial dynamics in soil that govern carbon (C) cycling are not well understood. Life history strategies explain the complex dynamics of bacterial populations and activities based on trade-offs in energy allocation to growth, resource acquisition, and survival. Such trade-offs influence the fate of soil C, but their genomic basis remains poorly characterized. We used multisubstrate metagenomic DNA stable isotope probing to link genomic features of bacteria to their C acquisition and growth dynamics. We identify several genomic features associated with patterns of bacterial C acquisition and growth, notably genomic investment in resource acquisition and regulatory flexibility. Moreover, we identify genomic trade-offs defined by numbers of transcription factors, membrane transporters, and secreted products, which match predictions from life history theory. We further show that genomic investment in resource acquisition and regulatory flexibility can predict bacterial ecological strategies in soil. IMPORTANCE Soil microbes are major players in the global carbon cycle, yet we still have little understanding of how the carbon cycle operates in soil communities. A major limitation is that carbon metabolism lacks discrete functional genes that define carbon transformations. Instead, carbon transformations are governed by anabolic processes associated with growth, resource acquisition, and survival. We use metagenomic stable isotope probing to link genome information to microbial growth and carbon assimilation dynamics as they occur in soil. From these data, we identify genomic traits that can predict bacterial ecological strategies which define bacterial interactions with soil carbon.

The power of next-generation sequencing has resulted in an explosive growth in the number of projects aiming to understand the metagenomic diversity of complex microbial environments. The interdisciplinary nature of this microbiome research community, along with the absence of reporting standards for microbiome data and samples, poses a significant challenge for follow-up studies. Commonly used names of metagenomes and metatranscriptomes in public databases currently lack the essential information necessary to accurately describe and classify the underlying samples, which makes a comparative analysis difficult to conduct and often results in misclassified sequences in data repositories. The Genomes OnLine Database (GOLD) (https:// gold.jgi.doe.gov/) at the Department of Energy Joint Genome Institute has been at the forefront of addressing this challenge by developing a standardized nomenclature system for naming microbiome samples. GOLD, currently in its twenty-fifth anniversary, continues to enrich the research community with hundreds of thousands of metagenomes and metatranscriptomes with well-curated and easy-to-understand names. Through this manuscript, we describe the overall naming process that can be easily adopted by researchers worldwide. Additionally, we propose the use of this naming system as a best practice for the scientific community to facilitate better interoperability and reusability of microbiome data.

BACKGROUNDLower respiratory tract infection (LRTI) is a leading cause of death in children worldwide. LRTI diagnosis is challenging because noninfectious respiratory illnesses appear clinically similar and because existing microbiologic tests are often falsely negative or detect incidentally carried microbes, resulting in antimicrobial overuse and adverse outcomes. Lower airway metagenomics has the potential to detect host and microbial signatures of LRTI. Whether it can be applied at scale and in a pediatric population to enable improved diagnosis and treatment remains unclear.METHODSWe used tracheal aspirate RNA-Seq to profile host gene expression and respiratory microbiota in 261 children with acute respiratory failure. We developed a gene expression classifier for LRTI by training on patients with an established diagnosis of LRTI (n = 117) or of noninfectious respiratory failure (n = 50). We then developed a classifier that integrates the host LRTI probability, abundance of respiratory viruses, and dominance in the lung microbiome of bacteria/fungi considered pathogenic by a rules-based algorithm.RESULTSThe host classifier achieved a median AUC of 0.967 by cross-validation, driven by activation markers of T cells, alveolar macrophages, and the interferon response. The integrated classifier achieved a median AUC of 0.986 and increased the confidence of patient classifications. When applied to patients with an uncertain diagnosis (n = 94), the integrated classifier indicated LRTI in 52% of cases and nominated likely causal pathogens in 98% of those.CONCLUSIONLower airway metagenomics enables accurate LRTI diagnosis and pathogen identification in a heterogeneous cohort of critically ill children through integration of host, pathogen, and microbiome features.FUNDINGSupport for this study was provided by the Eunice Kennedy Shriver National Institute of Child Health and Human Development and the National Heart, Lung, and Blood Institute (UG1HD083171, 1R01HL124103, UG1HD049983, UG01HD049934, UG1HD083170, UG1HD050096, UG1HD63108, UG1HD083116, UG1HD083166, UG1HD049981, K23HL138461, and 5R01HL155418) as well as by the Chan Zuckerberg Biohub.