Metagenomic Research Articles

Isolate-anchored comparisons reveal evolutionary and functional differentiation across SAR86 marine bacteria.

SAR86 is one of the most abundant groups of bacteria in the global surface ocean. However, since its discovery over 30years ago, it has remained recalcitrant to isolation and many details regarding this group are still unknown. Here we report the cellular characteristics from the first SAR86 isolate brought into culture, Magnimaribacter mokuoloeensis strain HIMB1674, and use its closed genome in concert with over 700 environmental genomes to assess the phylogenomic and functional characteristics of this order-level lineage of marine Gammaproteobacteria. The SAR86 order Magnimaribacterales invest significant genomic resources into the capacity for $\beta$-oxidation, which is present in most genomes in high gene copy numbers. This cyclical set of reactions appears to be fed by components of cell membranes that include lipids such as phosphatidylcholine, phosphatidylethanolamine, glycolipids, and sulfolipids. In addition to the widespread capacity to degrade the side chain of steroidal compounds via $\beta$-oxidation, several SAR86 sublineages also appear able to fully degrade the steroid polycyclic ring structure as well as other aromatic, polycyclic, and heterocyclic molecules. Read recruitment from publicly available metagenomes reveals that the Magnimaribacterales compose up to 6% of the global surface ocean microbial community. Only a subset of genera drive these high relative abundances, with some more globally dominant and others restricted to specific oceanic regions. This study provides an unprecedented foundation through which to understand this highly abundant yet poorly understood lineage of marine bacteria, and charts a path to bring more representatives of this order into laboratory culture.

Acidification-based mineral weathering mechanism involves a glucose/methanol/choline oxidoreductase in Caballeronia mineralivorans PML1(12).

This work deciphers the molecular and genetic bases used by strain PML1(12) of Caballeronia mineralivorans to weather minerals. Through bioinformatics analyses, we identified a total of four GMC-FAD oxidoreductases in the genome of strain PML1(12) and a putative PQQ-dependent glucose dehydrogenase. Through a combination of physiological and geochemical analyses, we revealed that one of them (i.e., GMC3) was the enzyme responsible for the acidification-based mineral weathering mechanism used by strain PML1(12). To date, a single representative of this enzyme family has been identified in the effective mineral-weathering bacterial strain Collimonas pratensis PMB3(1). Phylogenetic analyses revealed that this new system appeared conserved in the Collimonas genus. The new findings presented in this work demonstrate that GMC oxidoreductases can have an active role in other effective MWe bacteria outside of collimonads and that Caballeronia are capable of weathering minerals using this type of enzyme. Our findings offer relevant information for different fields of research, such as environmental genomics, microbiology, chemistry, evolutionary biology, and soil sciences.

Clinical and Genomic Epidemiology of Coxsackievirus A21 and Enterovirus D68 in Homeless Shelters, King County, Washington, USA, 2019-2021.

Congregate homeless shelters are disproportionately affected by infectious disease outbreaks. We describe enterovirus epidemiology across 23 adult and family shelters in King County, Washington, USA, during October 2019-May 2021, by using repeated cross-sectional respiratory illness and environmental surveillance and viral genome sequencing. Among 3,281 participants >3 months of age, we identified coxsackievirus A21 (CVA21) in 39 adult residents (3.0% [95% CI 1.9%-4.8%] detection) across 7 shelters during October 2019-February 2020. We identified enterovirus D68 (EV-D68) in 5 adult residents in 2 shelters during October-November 2019. Of 812 environmental samples, 1 was EV-D68-positive and 5 were CVA21-positive. Other enteroviruses detected among residents, but not in environmental samples, included coxsackievirus A6/A4 in 3 children. No enteroviruses were detected during April 2020-May 2021. Phylogenetically clustered CVA21 and EV-D68 cases occurred in some shelters. Some shelters also hosted multiple CVA21 lineages.

Diversity and distribution of a prevalent Microviridae group across the global oceans

Small single-stranded DNA phages of the Microviridae family are diverse and prevalent in oceans. Our understanding of Microviridae phages that infect the ecologically important marine Roseobacter is currently limited, comprising few isolates. Here, we report six roseophages that infect Roseobacter RCA strains. Genomic and phylogenetic analyses revealed that they were new members of the previously identified subfamily Occultatumvirinae. Additionally, 232 marine uncultivated virus genomes (UViGs) affiliated to Occultatumvirinae were obtained from environmental genome datasets. Phylogenomic analysis revealed that marine Occultatumvirinae phages could be further grouped into 11 subgroups. Moreover, meta-omics based read-mapping analysis showed that Occultatumvirinae phages were globally distributed, with two low G + C subgroups showing the most prevalent distribution. Furthermore, one phage in subgroup 2 was found to be extremely ubiquitous. Overall, this study expands our understanding of the diversity and ecology of the Occultatumvirinae microviruses in the ocean and highlights their ecological impacts.

Chromosome-scale genome assembly and annotation of the tetraploid potato cultivar Diacol Capiro adapted to the Andean region.

Potato (Solanum tuberosum) is an essential crop for food security and is ranked as the third most important crop worldwide for human consumption. The Diacol Capiro cultivar holds the dominant position in Colombian cultivation, primarily catering to the food processing industry. This highly heterozygous, autotetraploid cultivar belongs to the Andigenum group and it stands out for its adaptation to a wide variety of environments spanning altitudes from 1,800 to 3,200 meters above sea level. Here, a chromosome-scale assembly, referred to as DC, is presented for this cultivar. The assembly was generated by combining circular consensus sequencing with proximity ligation Hi-C for the scaffolding and represents 2.369 Gb with 48 pseudochromosomes covering 2,091 Gb and an anchor rate of 88.26%. The reference genome metrics, including an N50 of 50.5 Mb, a BUSCO (Benchmarking Universal Single-Copy Orthologue) score of 99.38%, and an Long Terminal Repeat Assembly Index score of 13.53, collectively signal the achieved high assembly quality. A comprehensive annotation yielded a total of 154,114 genes, and the associated BUSCO score of 95.78% for the annotated sequences attests to their completeness. The number of predicted NLR (Nucleotide-Binding and Leucine-Rich-Repeat genes) was 2107 with a large representation of NBARC (for nucleotide binding domain shared by Apaf-1, certain R gene products, and CED-4) containing domains (99.85%). Further comparative analysis of the proposed annotation-based assembly with high-quality known potato genomes, showed a similar genome metrics with differences in total gene numbers related to the ploidy status. The genome assembly and annotation of DC presented in this study represent a valuable asset for comprehending potato genetics. This resource aids in targeted breeding initiatives and contributes to the creation of enhanced, resilient, and more productive potato varieties, particularly beneficial for countries in Latin America.

Abstracts From the 55th Annual Meeting of the Environmental Mutagenesis and Genomics Society, September 7 – 11, 2024 ‐ Palm Springs, CA, USA, An Oasis of DNA Discoveries

Abstracts From the 55th Annual Meeting of the Environmental Mutagenesis and Genomics Society, September 7 – 11, 2024 ‐ Palm Springs, CA, USA, An Oasis of DNA Discoveries

MOSCA 2.0: A bioinformatics framework for metagenomics, metatranscriptomics and metaproteomics data analysis and visualization.

The analysis of meta-omics data requires the utilization of several bioinformatics tools and proficiency in informatics. The integration of multiple meta-omics data is even more challenging, and the outputs of existing bioinformatics solutions are not always easy to interpret. Here, we present a meta-omics bioinformatics pipeline, Meta-Omics Software for Community Analysis (MOSCA), which aims to overcome these limitations. MOSCA was initially developed for analysing metagenomics (MG) and metatranscriptomics (MT) data. Now, it also performs MG and metaproteomics (MP) integrated analysis, and MG/MT analysis was upgraded with an additional iterative binning step, metabolic pathways mapping, and several improvements regarding functional annotation and data visualization. MOSCA handles raw sequencing data and mass spectra and performs pre-processing, assembly, annotation, binning and differential gene/protein expression analysis. MOSCA shows taxonomic and functional analysis in large tables, performs metabolic pathways mapping, generates Krona plots and shows gene/protein expression results in heatmaps, improving omics data visualization. MOSCA is easily run from a single command while also providing a web interface (MOSGUITO). Relevant features include an extensive set of customization options, allowing tailored analyses to suit specific research objectives, and the ability to restart the pipeline from intermediary checkpoints using alternative configurations. Two case studies showcased MOSCA results, giving a complete view of the anaerobic microbial communities from anaerobic digesters and insights on the role of specific microorganisms. MOSCA represents a pivotal advancement in meta-omics research, offering an intuitive, comprehensive, and versatile solution for researchers seeking to unravel the intricate tapestry of microbial communities.

The smallest angiosperm genomes may be the price for effective traps of bladderworts.

Species of the carnivorous family Lentibulariaceae exhibit the smallest genomes in flowering plants. We explored the hypothesis that their minute genomes result from the unique mitochondrial cytochrome c oxidase (COX) mutation. The mutation may boost mitochondrial efficiency, which is especially useful for suction-bladder traps of Utricularia, but also increase DNA-damaging reactive oxygen species, leading to genome shrinkage through deletion-biased DNA repair. We aimed to explore this mutation's impact on genome size, providing insights into genetic mutation roles in plant genome evolution under environmental pressures. We compiled and measured genome and mean chromosome sizes for 127 and 67 species, respectively, representing all three genera (Genlisea, Pinguicula, and Utricularia) of Lentibulariaceae. We also isolated and analyzed COX sequences to detect the mutation. Through phylogenetic regressions and Ornstein-Uhlenbeck models of trait evolution, we assessed the impact of the COX mutation on the genome and chromosome sizes across the family. Our findings reveal significant correlations between the COX mutations and smaller genome and chromosome sizes. Specifically, species carrying the ancestral COX sequence exhibited larger genomes and chromosomes than those with the mutation. This evidence supports the notion that the COX mutation contributes to genome downsizing, with statistical analyses confirming a directional evolution towards smaller genomes in species harboring these mutations. Our study confirms that the COX mutation in Lentibulariaceae is associated with genome downsizing, likely driven by increased reactive oxygen species production and subsequent DNA damage requiring deletion-biased repair mechanisms. While boosting mitochondrial energy output, this genetic mutation compromises genome integrity and may potentially affect recombination rates, illustrating a complex trade-off between evolutionary advantages and disadvantages. Our results highlight the intricate processes by which genetic mutations and environmental pressures shape genome size evolution in carnivorous plants.

Archaeology meets environmental genomics: implementing sedaDNA in the study of the human past

Sedimentary ancient DNA (sedaDNA) has become one of the standard applications in the field of paleogenomics in recent years. It has been used for paleoenvironmental reconstructions, detecting the presence of prehistoric species in the absence of macro remains and even investigating the evolutionary history of a few species. However, its application in archaeology has been limited and primarily focused on humans. This article argues that sedaDNA holds significant potential in addressing key archaeological questions concerning the origins, lifestyles, and environments of past human populations. Our aim is to facilitate the integration of sedaDNA into the standard workflows in archaeology as a transformative tool, thereby unleashing its full potential for studying the human past. Ultimately, we not only underscore the challenges inherent in the sedaDNA field but also provide a research agenda for essential enhancements needed for implementing sedaDNA into the archaeological workflow.

Environmental activity-based protein profiling for function-driven enzyme discovery from natural communities

BackgroundMicrobial communities are important drivers of global biogeochemical cycles, xenobiotic detoxification, as well as organic matter decomposition. Their major metabolic role in ecosystem functioning is ensured by a unique set of enzymes, providing a tremendous yet mostly hidden enzymatic potential. Exploring this enzymatic repertoire is therefore not only relevant for a better understanding of how microorganisms function in their natural environment, and thus for ecological research, but further turns microbial communities, in particular from extreme habitats, into a valuable resource for the discovery of novel enzymes with potential applications in biotechnology. Different strategies for their uncovering such as bioprospecting, which relies mainly on metagenomic approaches in combination with sequence-based bioinformatic analyses, have emerged; yet accurate function prediction of their proteomes and deciphering the in vivo activity of an enzyme remains challenging.ResultsHere, we present environmental activity-based protein profiling (eABPP), a multi-omics approach that extends genome-resolved metagenomics with mass spectrometry-based ABPP. This combination allows direct profiling of environmental community samples in their native habitat and the identification of active enzymes based on their function, even without sequence or structural homologies to annotated enzyme families. eABPP thus bridges the gap between environmental genomics, correct function annotation, and in vivo enzyme activity. As a showcase, we report the successful identification of active thermostable serine hydrolases from eABPP of natural microbial communities from two independent hot springs in Kamchatka, Russia.ConclusionsBy reporting enzyme activities within an ecosystem in their native state, we anticipate that eABPP will not only advance current methodological approaches to sequence homology-guided enzyme discovery from environmental ecosystems for subsequent biocatalyst development but also contributes to the ecological investigation of microbial community interactions by dissecting their underlying molecular mechanisms.

Custom scoring based on ecological topology of gut microbiota associated with cancer immunotherapy outcome

The gut microbiota influences the clinical responses of cancer patients to immunecheckpoint inhibitors (ICIs). However, there is no consensus definition of detrimental dysbiosis. Based on metagenomics (MG) sequencing of 245 non-small cell lung cancer (NSCLC) patient feces, we constructed species-level co-abundance networks that were clustered into species-interacting groups (SIGs) correlating with overall survival. Thirty-seven and forty-five MG species (MGSs) were associated with resistance (SIG1) and response (SIG2) to ICIs, respectively. When combined with the quantification of Akkermansia species, this procedure allowed a person-based calculation of a topological score (TOPOSCORE) that was validated in an additional 254 NSCLC patients and in 216 genitourinary cancer patients. Finally, this TOPOSCORE was translated into a 21-bacterial probe set-based qPCR scoring that was validated in a prospective cohort of NSCLC patients as well as in colorectal and melanoma patients. This approach could represent a dynamic diagnosis tool for intestinal dysbiosis to guide personalized microbiota-centered interventions.

Cross-continental environmental and genome-wide association study on children and adolescent anxiety and depression.

Anxiety and depression in children and adolescents warrant special attention as a public health concern given their devastating and long-term effects on development and mental health. Multiple factors, ranging from genetic vulnerabilities to environmental stressors, influence the risk for the disorders. This study aimed to understand how environmental factors and genomics affect children and adolescents anxiety and depression across three cohorts: Adolescent Brain and Cognitive Development Study (US, age of 9-10; N=11,875), Consortium on Vulnerability to Externalizing Disorders and Addictions (INDIA, age of 6-17; N=4,326) and IMAGEN (EUROPE, age of 14; N=1888). We performed data harmonization and identified the environmental impact on anxiety/depression using a linear mixed-effect model, recursive feature elimination regression, and the LASSO regression model. Subsequently, genome-wide association analyses with consideration of significant environmental factors were performed for all three cohorts by mega-analysis and meta-analysis, followed by functional annotations. The results showed that multiple environmental factors contributed to the risk of anxiety and depression during development, where early life stress and school support index had the most significant and consistent impact across all three cohorts. In both meta, and mega-analysis, SNP rs79878474 in chr11p15 emerged as a particularly promising candidate associated with anxiety and depression, despite not reaching genomic significance. Gene set analysis on the common genes mapped from top promising SNPs of both meta and mega analyses found significant enrichment in regions of chr11p15 and chr3q26, in the function of potassium channels and insulin secretion, in particular Kv3, Kir-6.2, SUR potassium channels encoded by the KCNC1, KCNJ11, and ABCCC8 genes respectively, in chr11p15. Tissue enrichment analysis showed significant enrichment in the small intestine, and a trend of enrichment in the cerebellum. Our findings provide evidences of consistent environmental impact from early life stress and school support index on anxiety and depression during development and also highlight the genetic association between mutations in potassium channels, which support the stress-depression connection via hypothalamic-pituitary-adrenal axis, along with the potential modulating role of potassium channels.

First time whole genome sequencing of Mycobacterium bovis from the environment supports transmission at the animal-environment interface

Spreading of Mycobacterium bovis causing animal tuberculosis (TB) at livestock-wildlife-environment interfaces remains a significant problem. Recently, we provided evidence of widespread environmental contamination of an endemic animal TB setting with viable and dormant M. bovis cells able to recover metabolic activity, making indirect transmission via environmental contamination plausible. We now report the first whole genome sequences of M. bovis recovered from the environment. We establish epidemiological links at the environment-animal interface by phylogenomic comparison of these M. bovis genomes with those isolated from livestock and wild ungulates from the same area. Environmental and animal genomes are highly intertwined and distribute similarly into the same M. bovis lineages, supporting several instances of environmental contamination. This study provides compelling evidence of M. bovis excretion into the environment and viability maintenance, supporting the environment as a potential source of new infection. These insights have clear implications for policy formulation, advocating environmental surveillance and an ecosystem perspective in TB control programs. Environmental ImplicationWe report the first whole genome sequences of M. bovis from the environment and establish epidemiological links at the environment-animal interface, demonstrating close phylogenomic relatedness of animal and environmental M. bovis. Definitive evidence of M. bovis excretion into the environment with viability maintenance is provided, supporting the environment as a potential source of new infection. Implications of this work include methodological innovations offering a tool to resolve indirect transmission chains and support customized biosecurity measures. Policy formulation aiming at the control of animal tuberculosis and cost mitigation should consider these findings, encouraging environmental surveillance in official eradication programmes.

Low-level PM2.5 induces the occurrence of early pulmonary injury by regulating circ_0092363

The significant correlation between particulate matter with aerodynamic diameters of ≤ 2.5 µm (PM2.5) and the high morbidity and mortality of respiratory diseases has become the consensus of the research. Epidemiological studies have clearly pointed out that there is no safe concentration of PM2.5, and mechanism studies have also shown that exposure to PM2.5 will first cause pulmonary inflammation. Therefore, the purpose of this study is to explore the mechanism of early lung injury induced by low-level PM2.5 from the perspective of epigenetics. Based on the previous results of population samples, combined with an in vitro/vivo exposure model of PM2.5, it was found that low-level PM2.5 promoted the transport of circ_0092363 from intracellular to extracellular spaces.The decreased expression of intracellular circ_0092363 resulted in reduced absorption of miR-31-5p, leading to inhibition of Rho associated coiled-coil containing protein kinase 1 (ROCK1) and the subsequent abnormal expression of tight junction proteins such as Zonula occludens protein 1 (ZO-1) and Claudin-1, ultimately inducing the occurrence of early pulmonary injury. Furthermore, this study innovatively introduced organoid technology and conducted a preliminary exploration for a study of the relationship among environmental exposure genomics, epigenetics and disease genomics in organoids. The role of circ_0092363 in early pulmonary injury induced by low-level PM2.5 was elucidated, and its value as a potential diagnostic biomarker was confirmed.

Influence of Environmental Factors and Genome Diversity on Cumulative COVID-19 Cases in the Highland Region of China: Comparative Correlational Study.

The novel coronavirus SARS-CoV-2 caused the global COVID-19 pandemic. Emerging reports support lower mortality and reduced case numbers in highland areas; however, comparative studies on the cumulative impact of environmental factors and viral genetic diversity on COVID-19 infection rates have not been performed to date. The aims of this study were to determine the difference in COVID-19 infection rates between high and low altitudes, and to explore whether the difference in the pandemic trend in the high-altitude region of China compared to that of the lowlands is influenced by environmental factors, population density, and biological mechanisms. We examined the correlation between population density and COVID-19 cases through linear regression. A zero-shot model was applied to identify possible factors correlated to COVID-19 infection. We further analyzed the correlation of meteorological and air quality factors with infection cases using the Spearman correlation coefficient. Mixed-effects multiple linear regression was applied to evaluate the associations between selected factors and COVID-19 cases adjusting for covariates. Lastly, the relationship between environmental factors and mutation frequency was evaluated using the same correlation techniques mentioned above. Among the 24,826 confirmed COVID-19 cases reported from 40 cities in China from January 23, 2020, to July 7, 2022, 98.4% (n=24,430) were found in the lowlands. Population density was positively correlated with COVID-19 cases in all regions (ρ=0.641, P=.003). In high-altitude areas, the number of COVID-19 cases was negatively associated with temperature, sunlight hours, and UV index (P=.003, P=.001, and P=.009, respectively) and was positively associated with wind speed (ρ=0.388, P<.001), whereas no correlation was found between meteorological factors and COVID-19 cases in the lowlands. After controlling for covariates, the mixed-effects model also showed positive associations of fine particulate matter (PM2.5) and carbon monoxide (CO) with COVID-19 cases (P=.002 and P<.001, respectively). Sequence variant analysis showed lower genetic diversity among nucleotides for each SARS-CoV-2 genome (P<.001) and three open reading frames (P<.001) in high altitudes compared to 300 sequences analyzed from low altitudes. Moreover, the frequencies of 44 nonsynonymous mutations and 32 synonymous mutations were significantly different between the high- and low-altitude groups (P<.001, mutation frequency>0.1). Key nonsynonymous mutations showed positive correlations with altitude, wind speed, and air pressure and showed negative correlations with temperature, UV index, and sunlight hours. By comparison with the lowlands, the number of confirmed COVID-19 cases was substantially lower in high-altitude regions of China, and the population density, temperature, sunlight hours, UV index, wind speed, PM2.5, and CO influenced the cumulative pandemic trend in the highlands. The identified influence of environmental factors on SARS-CoV-2 sequence variants adds knowledge of the impact of altitude on COVID-19 infection, offering novel suggestions for preventive intervention.

Time-shifted expression of acetoclastic and methylotrophic methanogenesis by a single Methanosarcina genomospecies predominates the methanogen dynamics in Philippine rice field soil

BackgroundThe final step in the anaerobic decomposition of biopolymers is methanogenesis. Rice field soils are a major anthropogenic source of methane, with straw commonly used as a fertilizer in rice farming. Here, we aimed to decipher the structural and functional responses of the methanogenic community to rice straw addition during an extended anoxic incubation (120 days) of Philippine paddy soil. The research combined process measurements, quantitative real-time PCR and RT-PCR of particular biomarkers (16S rRNA, mcrA), and meta-omics (environmental genomics and transcriptomics).ResultsThe analysis methods collectively revealed two major bacterial and methanogenic activity phases: early (days 7 to 21) and late (days 28 to 60) community responses, separated by a significant transient decline in microbial gene and transcript abundances and CH4 production rate. The two methanogenic activity phases corresponded to the greatest rRNA and mRNA abundances of the Methanosarcinaceae but differed in the methanogenic pathways expressed. While three genetically distinct Methanosarcina populations contributed to acetoclastic methanogenesis during the early activity phase, the late activity phase was defined by methylotrophic methanogenesis performed by a single Methanosarcina genomospecies. Closely related to Methanosarcina sp. MSH10X1, mapping of environmental transcripts onto metagenome-assembled genomes (MAGs) and population-specific reference genomes revealed this genomospecies as the key player in acetoclastic and methylotrophic methanogenesis. The anaerobic food web was driven by a complex bacterial community, with Geobacteraceae and Peptococcaceae being putative candidates for a functional interplay with Methanosarcina. Members of the Methanocellaceae were the key players in hydrogenotrophic methanogenesis, while the acetoclastic activity of Methanotrichaceae members was detectable only during the very late community response.ConclusionsThe predominant but time-shifted expression of acetoclastic and methylotrophic methanogenesis by a single Methanosarcina genomospecies represents a novel finding that expands our hitherto knowledge of the methanogenic pathways being highly expressed in paddy soils.7k5euw-hvBknz_jds4sXXfVideo

Influence of microbiota-associated metabolic reprogramming on clinical outcome in patients with melanoma from the randomized adjuvant dendritic cell-based MIND-DC trial

Tumor immunosurveillance plays a major role in melanoma, prompting the development of immunotherapy strategies. The gut microbiota composition, influencing peripheral and tumoral immune tonus, earned its credentials among predictors of survival in melanoma. The MIND-DC phase III trial (NCT02993315) randomized (2:1 ratio) 148 patients with stage IIIB/C melanoma to adjuvant treatment with autologous natural dendritic cell (nDC) or placebo (PL). Overall, 144 patients collected serum and stool samples before and after 2 bimonthly injections to perform metabolomics (MB) and metagenomics (MG) as prespecified exploratory analysis. Clinical outcomes are reported separately. Here we show that different microbes were associated with prognosis, with the health-related Faecalibacterium prausnitzii standing out as the main beneficial taxon for no recurrence at 2 years (p = 0.008 at baseline, nDC arm). Therapy coincided with major MB perturbations (acylcarnitines, carboxylic and fatty acids). Despite randomization, nDC arm exhibited MG and MB bias at baseline: relative under-representation of F. prausnitzii, and perturbations of primary biliary acids (BA). F. prausnitzii anticorrelated with BA, medium- and long-chain acylcarnitines. Combined, these MG and MB biomarkers markedly determined prognosis. Altogether, the host-microbial interaction may play a role in localized melanoma. We value systematic MG and MB profiling in randomized trials to avoid baseline differences attributed to host-microbe interactions.

Late Pleistocene stickleback environmental genomes reveal the chronology of freshwater adaptation

Directly observing the chronology and tempo of adaptation in response to ecological change is rarely possible in natural ecosystems. Sedimentary ancient DNA (sedaDNA) has been shown to be a tractable source of genome-scale data of long-dead organisms1,2,3 and to thereby potentially provide an understanding of the evolutionary histories of past populations.4,5 To date, time series of ecosystem biodiversity have been reconstructed from sedaDNA, typically using DNA metabarcoding or shotgun sequence data generated from less than 1g of sediment.6,7 Here, we maximize sequence coverage by extracting DNA from ∼50× more sediment per sample than the majority of previous studies1,2,3 to achieve genotype resolution. From a time series of Late Pleistocene sediments spanning from a marine to freshwater ecosystem, we compare adaptive genotypes reconstructed from the environmental genomes of three-spined stickleback at key time points of this transition. We find a staggered temporal dynamic in which freshwater alleles at known loci of large effect in marine-freshwater divergence of three-spined stickleback (e.g., EDA)8 were already established during the brackish phase of the formation of the isolation basin. However, marine alleles were still detected across the majority of marine-freshwater divergence-associated loci, even after the complete isolation of the lake from marine ingression. Our retrospective approach to studying adaptation from environmental genomes of three-spined sticklebacks at the end of the last glacial period complements contemporary experimental approaches9,10,11 and highlights the untapped potential for retrospective "evolve and resequence" natural experiments using sedaDNA.

Large-scale single-virus genomics uncovers hidden diversity of river water viruses and diversified gene profiles.

Environmental viruses (primarily bacteriophages) are widely recognized as playing an important role in ecosystem homeostasis through the infection of host cells. However, the majority of environmental viruses are still unknown as their mosaic structure and frequent mutations in their sequences hinder genome construction in current metagenomics. To enable the large-scale acquisition of environmental viral genomes, we developed a new single-viral genome sequencing platform with microfluidic-generated gel beads. Amplification of individual DNA viral genomes in mass-produced gel beads allows high-throughput genome sequencing compared to conventional single-virus genomics. The sequencing analysis of river water samples yielded 1431 diverse viral single-amplified genomes, whereas viral metagenomics recovered 100 viral metagenome-assembled genomes at the comparable sequence depth. The 99.5% of viral single-amplified genomes were determined novel at the species level, most of which could not be recovered by a metagenomic assembly. The large-scale acquisition of diverse viral genomes identified protein clusters commonly detected in different viral strains, allowing the gene transfer to be tracked. Moreover, comparative genomics within the same viral species revealed that the profiles of various methyltransferase subtypes were diverse, suggesting an enhanced escape from host bacterial internal defense mechanisms. Our use of gel bead-based single-virus genomics will contribute to exploring the nature of viruses by accelerating the accumulation of draft genomes of environmental DNA viruses.

Environmental Compatibility and Genome Flexibility of Klebsiella oxytoca Isolated from Eight Species of Aquatic Animals

Klebsiella oxytoca is an emerging pathogen that can cause life-threatening infectious diseases in humans. Recently, we firstly reported for the first time the presence of K. oxytoca in edible aquatic animals. In this study, we further investigated its bacterial environmental fitness and genome evolution signatures. The results revealed that K. oxytoca isolates (n = 8), originating from eight species of aquatic animals, were capable of growing under a broad spectrum of environmental conditions (pH 4.5–8.5, 0.5–6.5% NaCl), with different biofilm formation and swimming mobility profiles. The genome sequences of the K. oxytoca isolates were determined (5.84–6.02 Mb, 55.07–56.06% GC content). Strikingly, numerous putative mobile genetic elements (MGEs), particularly genomic islands (GIs, n = 105) and prophages (n = 24), were found in the K. oxytoca genomes, which provided the bacterium with specific adaptation traits, such as resistance, virulence, and material metabolism. Interestingly, the identified prophage-related clusters were derived from Burkholderia spp., Enterobacter spp., Klebsiella spp., Pseudomonas spp., and Haemophilus spp., suggesting phage transmission across Klebsiella and the other four genera. Many strain-specific (n = 10–447) genes were present in the K. oxytoca genomes, whereas the CRISPR-Cas protein-encoding gene was absent, indicating likely active horizontal gene transfer (HGT) and considerable genome variation in K. oxytoca evolution. Overall, the results of this study are the first to demonstrate the environmental compatibility and genome flexibility of K. oxytoca of aquatic animal origins.