Search

Amphibians represent a diverse group of tetrapods, marked by deep divergence times between their three systematic orders and families. Studying amphibian biology through the genomics lens increases our understanding of the features of this animal class and that of other terrestrial vertebrates. The need for amphibian genomic resources is more urgent than ever due to the increasing threats to this group. Amphibians are one of the most imperiled taxonomic groups, with approximately 41% of species threatened with extinction due to habitat loss, changes in land use patterns, disease, climate change, and their synergistic effects. Amphibian genomic resources have provided a better understanding of ontogenetic diversity, tissue regeneration, diverse life history and reproductive modes, anti-predator strategies, and resilience and adaptive responses. They also serve as essential models for studying broad genomic traits, such as evolutionary genome expansions and contractions, as they exhibit the widest range of genome sizes among all animal taxa and possess multiple mechanisms of genetic sex determination. Despite these features, genome sequencing of amphibians has significantly lagged behind that of other vertebrates, primarily due to the challenges of assembling their large, repeat-rich genomes and the relative lack of societal support. The emergence of long-read sequencing technologies, combined with advanced molecular and computational techniques that improve scaffolding and reduce computational workloads, is now making it possible to address some of these challenges. To promote and accelerate the production and use of amphibian genomics research through international coordination and collaboration, we launched the Amphibian Genomics Consortium (AGC, https://mvs.unimelb.edu.au/amphibian-genomics-consortium ) in early 2023. This burgeoning community already has more than 282 members from 41 countries. The AGC aims to leverage the diverse capabilities of its members to advance genomic resources for amphibians and bridge the implementation gap between biologists, bioinformaticians, and conservation practitioners. Here we evaluate the state of the field of amphibian genomics, highlight previous studies, present challenges to overcome, and call on the research and conservation communities to unite as part of the AGC to enable amphibian genomics research to leap to the next level.

While numerous computational frameworks and workflows are available for recovering prokaryote and eukaryote genomes from metagenome data, only a limited number of pipelines are designed specifically for viromics analysis. With many viromics tools developed in the last few years alone, it can be challenging for scientists with limited bioinformatics experience to easily recover, evaluate quality, annotate genes, dereplicate, assign taxonomy, and calculate relative abundance and coverage of viral genomes using state-of-the-art methods and standards. Here, we describe Modular Viromics Pipeline (MVP) v.1.0, a user-friendly pipeline written in Python and providing a simple framework to perform standard viromics analyses. MVP combines multiple tools to enable viral genome identification, characterization of genome quality, filtering, clustering, taxonomic and functional annotation, genome binning, and comprehensive summaries of results that can be used for downstream ecological analyses. Overall, MVP provides a standardized and reproducible pipeline for both extensive and robust characterization of viruses from large-scale sequencing data including metagenomes, metatranscriptomes, viromes, and isolate genomes. As a typical use case, we show how the entire MVP pipeline can be applied to a set of 20 metagenomes from wetland sediments using only 10 modules executed via command lines, leading to the identification of 11,656 viral contigs and 8,145 viral operational taxonomic units (vOTUs) displaying a clear beta-diversity pattern. Further, acting as a dynamic wrapper, MVP is designed to continuously incorporate updates and integrate new tools, ensuring its ongoing relevance in the rapidly evolving field of viromics. MVP is available at https://gitlab.com/ccoclet/mvp and as versioned packages in PyPi and Conda.IMPORTANCEThe significance of our work lies in the development of Modular Viromics Pipeline (MVP), an integrated and user-friendly pipeline tailored exclusively for viromics analyses. MVP stands out due to its modular design, which ensures easy installation, execution, and integration of new tools and databases. By combining state-of-the-art tools such as geNomad and CheckV, MVP provides high-quality viral genome recovery and taxonomy and host assignment, and functional annotation, addressing the limitations of existing pipelines. MVP's ability to handle diverse sample types, including environmental, human microbiome, and plant-associated samples, makes it a versatile tool for the broader microbiome research community. By standardizing the analysis process and providing easily interpretable results, MVP enables researchers to perform comprehensive studies of viral communities, significantly advancing our understanding of viral ecology and its impact on various ecosystems.

Fungal fermentation of food and agricultural by-products holds promise for improving food sustainability and security. However, the molecular basis of fungal waste-to-food upcycling remains poorly understood. Here we use a multi-omics approach to characterize oncom, a fermented food traditionally produced from soymilk by-products in Java, Indonesia. Metagenomic sequencing of samples from small-scale producers in Western Java indicated that the fungus Neurospora intermedia dominates oncom. Further transcriptomic, metabolomic and phylogenomic analysis revealed that oncom-derived N. intermedia utilizes pectin and cellulose degradation during fermentation and belongs to a genetically distinct subpopulation associated with human-generated by-products. Finally, we found that N. intermedia grew on diverse by-products such as fruit and vegetable pomace and plant-based milk waste, did not encode mycotoxins, and could create foods that were positively perceived by consumers outside Indonesia. These results showcase the traditional significance and future potential of fungal fermentation for creating delicious and nutritious foods from readily available by-products.

This article summarises the state of the science on the role of the gut microbiota (GM) in diabetes from a recent international expert forum organised by Diabetes, Diabetes Care, and Diabetologia, which was held at the European Association for the Study of Diabetes 2023 Annual Meeting in Hamburg, Germany. Forum participants included clinicians and basic scientists who are leading investigators in the field of the intestinal microbiome and metabolism. Their conclusions were as follows: (1) the GM may be involved in the pathophysiology of type 2 diabetes, as microbially produced metabolites associate both positively and negatively with the disease, and mechanistic links of GM functions (e.g. genes for butyrate production) with glucose metabolism have recently emerged through the use of Mendelian randomisation in humans; (2) the highly individualised nature of the GM poses a major research obstacle, and large cohorts and a deep-sequencing metagenomic approach are required for robust assessments of associations and causation; (3) because single time point sampling misses intraindividual GM dynamics, future studies with repeated measures within individuals are needed; and (4) much future research will be required to determine the applicability of this expanding knowledge to diabetes diagnosis and treatment, and novel technologies and improved computational tools will be important to achieve this goal.

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.

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.

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 Exabiome project seeks to improve the understanding of microbiomes through the development of methods for accelerating metagenomic science using exascale computing. This article gives an overview of scientific impact of the three components of the project: metagenome assembly, protein family detection, and comparative analysis of metagenomes. Exabiome developed MetaHipMer, the only metagenome assembler capable of scaling to full exascale systems. MetaHipMer has enabled ground-breaking assemblies on the Frontier supercomputer, with many scientific benefits, such as the discovery of rare species and viral genomes. To investigate protein families, Exabiome developed two exascale tools, PASTIS and HipMCL. Together, these can utilize exascale resources to understand the functional diversity of billions of dark matter proteins and novel protein families. For comparative analysis, Exabiome developed kmerprof, a tool that can be used to compare huge metagenomes for many different scientific purposes, for example, grouping human microbiomes according to body location.

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.

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 ecosystem's 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 collection's 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 patient's 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 model's 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.

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.

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.

DNA/RNA-stable isotope probing (SIP) is a powerful tool to link in situ microbial activity to sequencing data. Every SIP dataset captures distinct information about microbial community metabolism, process rates, and population dynamics, offering valuable insights for a wide range of research questions. Data reuse maximizes the information derived from the labor and resource-intensive SIP approaches. Yet, a review of publicly available SIP sequencing metadata showed that critical information necessary for reproducibility and reuse was often missing. Here, we outline the Minimum Information for any Stable Isotope Probing Sequence (MISIP) according to the Minimum Information for any (x) Sequence (MIxS) framework and include examples of MISIP reporting for common SIP experiments. Our objectives are to expand the capacity of MIxS to accommodate SIP-specific metadata and guide SIP users in metadata collection when planning and reporting an experiment. The MISIP standard requires 5 metadata fields-isotope, isotopolog, isotopolog label, labeling approach, and gradient position-and recommends several fields that represent best practices in acquiring and reporting SIP sequencing data (e.g., gradient density and nucleic acid amount). The standard is intended to be used in concert with other MIxS checklists to comprehensively describe the origin of sequence data, such as for marker genes (MISIP-MIMARKS) or metagenomes (MISIP-MIMS), in combination with metadata required by an environmental extension (e.g., soil). The adoption of the proposed data standard will improve the reuse of any sequence derived from a SIP experiment and, by extension, deepen understanding of in situ biogeochemical processes and microbial ecology.

The North Temperate Lakes Long-Term Ecological Research (NTL-LTER) program has been extensively used to improve understanding of how aquatic ecosystems respond to environmental stressors, climate fluctuations, and human activities. Here, we report on the metagenomes of samples collected between 2000 and 2019 from Lake Mendota, a freshwater eutrophic lake within the NTL-LTER site. We utilized the distributed metagenome assembler MetaHipMer to coassemble over 10 terabases (Tbp) of data from 471 individual Illumina-sequenced metagenomes. A total of 95,523,664 contigs were assembled and binned to generate 1,894 non-redundant metagenome-assembled genomes (MAGs) with ≥50% completeness and ≤10% contamination. Phylogenomic analysis revealed that the MAGs were nearly exclusively bacterial, dominated by Pseudomonadota (Proteobacteria, N = 623) and Bacteroidota (N = 321). Nine eukaryotic MAGs were identified by eukCC with six assigned to the phylum Chlorophyta. Additionally, 6,350 high-quality viral sequences were identified by geNomad with the majority classified in the phylum Uroviricota. This expansive coassembled metagenomic dataset provides an unprecedented foundation to advance understanding of microbial communities in freshwater ecosystems and explore temporal ecosystem dynamics.

Environmental DNA (eDNA) is becoming an increasingly important tool in diverse scientific fields from ecological biomonitoring to wastewater surveillance of viruses. The fundamental challenge in eDNA analyses has been the bioinformatical assignment of reads to taxonomic groups. It has long been known that full probabilistic methods for phylogenetic assignment are preferable, but unfortunately, such methods are computationally intensive and are typically inapplicable to modern next-generation sequencing data. We present a fast approximate likelihood method for phylogenetic assignment of DNA sequences. Applying the new method to several mock communities and simulated datasets, we show that it identifies more reads at both high and low taxonomic levels more accurately than other leading methods. The advantage of the method is particularly apparent in the presence of polymorphisms and/or sequencing errors and when the true species is not represented in the reference database.

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.

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.

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.

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.