Search NASA⌕ Search

SEARCH · Search NASA

Results for “Long-read sequencing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Integration of Yeast Episomal/Integrative Plasmid Causes Genotypic and Phenotypic Diversity and Improved Sesquiterpene Production in Metabolically Engineered Saccharomyces cerevisiae

The variability in phenotypic outcomes among biological replicates in engineered microbial factories presents a captivating mystery. Establishing the association between phenotypic variability and genetic drivers is important to solve this intricate puzzle. Here, we applied a previously developed auxin-inducible depletion of hexokinase 2 as a metabolic engineering strategy for improved nerolidol production in Saccharomyces cerevisiae, and biological replicates exhibit a dichotomy in nerolidol production of either 3.5 or 2.5 g L –1 nerolidol. Harnessing Oxford Nanopore’s long-read genomic sequencing, we reveal a potential genetic cause—the chromosome integration of a 2μ sequence-based yeast episomal plasmid, encoding the expression cassettes for nerolidol synthetic enzymes. This finding was reinforced through chromosome integration revalidation, engineering nerolidol and valencene production strains, and generating a diverse pool of yeast clones, each uniquely fingerprinted by gene copy numbers, plasmid integrations, other genomic rearrangements, protein expression levels, growth rate, and target product productivities. Τhe best clone in two strains produced 3.5 g L –1 nerolidol and ~0.96 g L –1 valencene. Comparable genotypic and phenotypic variations were also generated through the integration of a yeast integrative plasmid lacking 2μ sequences. Our work shows that multiple factors, including plasmid integration status, subchromosomal location, gene copy number, sesquiterpene synthase expression level, and genome rearrangement, together play a complicated determinant role on the productivities of sesquiterpene product. Integration of yeast episomal/integrative plasmids may be used as a versatile method for increasing the diversity and optimizing the efficiency of yeast cell factories, thereby uncovering metabolic control mechanisms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Phylogenetic analyses and reclassification of the oleaginous marine species Nannochloris sp. “desiccata" (Trebouxiophyceae, Chlorophyta), formerly Chlorella desiccata , supported by a high-quality genome assembly

Microalgae are diverse, with many gaps remaining in phylogenetic and physiological understanding. Thus, studying new microalgae species increases our broader comprehension of biological diversity, and evaluation of new candidates as algal production platforms can lead to improved productivity under a variety of cultivation conditions. Chlorella is a genus of fast-growing species often isolated from freshwater habitats and cultivated as a source of nutritional supplements. However, the use of freshwater increases competition with other freshwater needs. We identified Chlorella desiccata to be worthy of further investigation as a potential algae production strain, due to its isolation from a marine environment and its promising growth and biochemical composition properties. Long-read genomic sequencing was conducted for C. desiccata UTEX 2526, resulting in a high-quality, near chromosome level, diploid genome with an assembly length of 21.55 Mbp in only 18 contigs. We also report complete circular mitochondrial and chloroplast genomes. Phylogenomic and phylogenetic analyses using nuclear, chloroplast, 18S rRNA, and actin sequences revealed that this species clades within strains currently identified as Nannochloris (Trebouxiophyceae, Chlorophyta), leading to its reclassification as Nannochloris sp. “desiccata” UTEX 2526. The mode of cell division for this species is autosporulation, differing from the type species N. bacillaris. As has occurred across multiple microalgae genera, there are repeated examples of Nannochloris species reclassification in the literature. This high-quality genome assembly and phylogenetic analysis of the potential algal production strain Nannochloris sp. “desiccata” UTEX 2526 provides an important reference and useful tool for further studying this region of the phylogenetic tree.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Comparative physiological and genomic characterization of a novel Nitrobacter vulgaris strain from a nitrate-contaminated subsurface

Nitrite-oxidizing bacteria (NOB) represent a crucial node in the global nitrogen cycle. By catalyzing the second step of nitrification—the oxidation of nitrite to nitrate to generate energy for growth—NOB activity controls the fate of nitrite (NO 2 - ) in aerobic environments. Despite thriving in diverse environments, including soils, freshwater, marine ecosystems, subsurface habitats, and water treatment systems, organisms capable of nitrite oxidation are confined to Nitrobacter, Nitrospira, Nitrospina, Nitrotoga, and a few other specific lineages. The genus Nitrobacter, recognized for its facultative heterotrophic metabolism, is often associated with high-nitrogen environments. Here, we report the physiological characterization of a novel strain, Nitrobacter vulgaris strain MLSD-S22, isolated from a nitrate- and heavy-metal-contaminated subsurface. Growth inhibition experiments revealed that strain MLSD-S22 and the N. vulgaris type strain Z exhibited similar sensitivities to nitrite and nitrate, with nitrite being the most inhibitory. Microrespirometry demonstrated that the two N. vulgaris strains and Nitrobacter winogradskyi Nb-255 possessed higher affinities for nitrite and oxygen than previously reported for Nitrobacter, suggesting potential to compete in low-substrate environments. Long-read DNA sequencing provided a complete genome for strain MLSD-S22, revealing two plasmids and an intact nitrous oxide (N 2 O) reduction operon—an unexpected feature for Nitrobacter. While N 2 O reduction activity was not observed under the tested conditions, this discovery raises questions about the contribution of Nitrobacter NOB to the N 2 O sink. These findings broaden the physiological and genomic diversity of Nitrobacter, offering new insights into their adaptation strategies and providing a framework for future evaluation of their potential roles in nitrogen loss.

Nitrobacter↗
DOE OSTI2023

Optimizing Cell-based Antimicrobials through Pooled Genomic Libraries

DNA synthesis and assembly technologies ushered in through synthetic biology have great promise for biomanufacturing, bioremediation, and the development of living therapeutics. Unfortunately, predicting sequence to function relationships, including for biosynthetic pathways expressed in a new host organism, is difficult and often requires many iterative cycles of design, construction, and testing. We are working to develop data-driven approaches to identify the genetic determinants of growth defects and productivity for the expression of a cell-based antimicrobial. We assayed the growth, pigment production, and antimicrobial activity of a collection of over 10,000 genetic mutants of the violacein biosynthetic pathway and sequenced the genetic variation of these mutants. Through this project, we have developed an innovative codebase to automate the determination of pigmentation and antimicrobial clearing diameter for tens of thousands of genetic mutants cultivated on agar dishes. Further, we have written DNA sequence analysis code to demultiplex & provide consensus sequences from high-throughput PacBio long-read circular consensus sequencing (CCS) datasets. From this foundation, we plan to map DNA sequence to function to predict an optimal genetic design to maximize antimicrobial activity while minimizing deleterious growth effects. The workflows and algorithms developed through this project can be broadly applied to other engineered functions in microbes, uncovering sequence to function relationships for complex phenotypes where function impacts fitness.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Complete genomes of Asgard archaea reveal diverse integrated and mobile genetic elements

Asgard archaea are of great interest as the progenitors of Eukaryotes, but little is known about the mobile genetic elements (MGEs) that may shape their ongoing evolution. Here, we describe MGEs that replicate in Atabeyarchaeia, a wetland Asgard archaea lineage represented by two complete genomes. We used soil depth–resolved population metagenomic data sets to track 18 MGEs for which genome structures were defined and precise chromosome integration sites could be identified for confident host linkage. Additionally, we identified a complete 20.67 kbp circular plasmid and two family-level groups of viruses linked to Atabeyarchaeia, via CRISPR spacer targeting. Closely related 40 kbp viruses possess a hypervariable genomic region encoding combinations of specific genes for small cysteine-rich proteins structurally similar to restriction-homing endonucleases. One 10.9 kbp integrative conjugative element (ICE) integrates genomically into theAtabeyarchaeum deiterrae-1chromosome and has a 2.5 kbp circularizable element integrated within it. The 10.9 kbp ICE encodes an expressed Type IIG restriction-modification system with a sequence specificity matching an active methylation motif identified by Pacific Biosciences (PacBio) high-accuracy long-read (HiFi) metagenomic sequencing. Restriction-modification of Atabeyarchaeia differs from that of another coexisting Asgard archaea, Freyarchaeia, which has few identified MGEs but possesses diverse defense mechanisms, including DISARM and Hachiman, not found in Atabeyarchaeia. Overall, defense systems and methylation mechanisms of Asgard archaea likely modulate their interactions with MGEs, and integration/excision and copy number variation of MGEs in turn enable host genetic versatility.

Biochemistry & Molecular Biology↗
DOE OSTI2024

LevSeq: Rapid Generation of Sequence-Function Data for Directed Evolution and Machine Learning

Sequence-function data provides valuable information about the protein functional landscape but is rarely obtained during directed evolution campaigns. Here, we present Long-read every variant Sequencing (LevSeq), a pipeline that combines a dual barcoding strategy with nanopore sequencing to rapidly generate sequence-function data for entire protein-coding genes. LevSeq integrates into existing protein engineering workflows and comes with open-source software for data analysis and visualization. The pipeline facilitates data-driven protein engineering by consolidating sequence-function data to inform directed evolution and provide the requisite data for machine learning-guided protein engineering (MLPE). LevSeq enables quality control of mutagenesis libraries prior to screening, which reduces time and resource costs. Simulation studies demonstrate LevSeq’s ability to accurately detect variants under various experimental conditions. Lastly, we show LevSeq’s utility in engineering protoglobins for new-to-nature chemistry. Widespread adoption of LevSeq and sharing of the data will enhance our understanding of protein sequence-function landscapes and empower data-driven directed evolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Uncultivated Viral Populations Dominate Estuarine Viromes on the Spatiotemporal Scale

Viruses are ubiquitous and abundant in the oceans, and viral metagenomes (viromes) have been investigated extensively via several large-scale ocean sequencing projects. However, there have not been any systematic viromic studies in estuaries. Here, we investigated the viromes of the Delaware Bay and Chesapeake Bay, two Mid-Atlantic estuaries. Deep sequencing generated a total of 48,190 assembled viral sequences (>5 kb) and 26,487 viral populations (9,204 virus clusters and 17,845 singletons), including 319 circular viral contigs between 7.5 kb and 161.8 kb. Unknown viruses represented the vast majority of the dominant populations, while the composition of known viruses, such as pelagiphage and cyanophage, appeared to be relatively consistent across a wide range of salinity gradients and in different seasons. A difference between estuarine and ocean viromes was reflected by the proportions of Myoviridae, Podoviridae, Siphoviridae, Phycodnaviridae, and a few well-studied virus representatives. The difference in viral community between the Delaware Bay and Chesapeake Bay is significantly more pronounced than the difference caused by temperature or salinity, indicating strong local profiles caused by the unique ecology of each estuary. Interestingly, a viral contig similar to phages infecting Acinetobacter baumannii (“Iraqibacter”) was found to be highly abundant in the Delaware Bay but not in the Chesapeake Bay, the source of which is yet to be identified. Highly abundant viruses in both estuaries have close hits to viral sequences derived from the marine single-cell genomes or long-read single-molecule sequencing, suggesting that important viruses are still waiting to be discovered in the estuarine environment.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Luteolibacter sp. strain Populi

Luteolibacter sp. strain Populi is bacterium from the phylum Verrucomicrobiota, isolated from the rhizosphere of a black cottonwood tree, Populus trichocarpa, from the Cascade mountains in Washington. Its 6.6 Mb chromosome was completely sequenced using Oxford Nanopore long-reads and is predicted to encode 5301 proteins and 60 RNAs. The bacteria was isolated from the rhizosphere of a mature Populus trichocarpa from the Tieton riverwatershed of Washington state, USA (Lat: 46°42’9” N, Lon: 120°25 39’36” W). A rhizosphere sample (fine roots and adhering soil) was used to obtain a microbial fraction by centrifugation on Histodenz (12) and stained with 5µM Syto59 (Thermo Fisher Scientific Inc). A Cytopeia Influx cell sorter (BD, Franklin Lakes, NJ) was used to sort and array single cells (100 per plate) based on forward-side scatter and fluorescence intensity on asparagine-glucose nutrient agar (ATCC medium 184). The Luteolibacter sp. Populi genome sequence has been deposited in GenBank under the accession number CP161812. A draft genome annotated with Prokka and DRAM is available in this Narrative as Luteolibacter_sp_Prokka.240711.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Prevalence and diversity of TAL effector-like proteins in fungal endosymbiotic Mycetohabitans spp.

EndofungalMycetohabitans(formerlyBurkholderia) spp. rely on a type III secretion system to deliver mostly unidentified effector proteins when colonizing their host fungus,Rhizopus microsporus. The one known secreted effector family fromMycetohabitansconsists of homologues of transcription activator-like (TAL) effectors, which are used by plant pathogenicXanthomonasandRalstoniaspp. to activate host genes that promote disease. These ‘BurkholderiaTAL-like (Btl)’ proteins bind corresponding specific DNA sequences in a predictable manner, but their genomic target(s) and impact on transcription in the fungus are unknown. Recent phenotyping of Btl mutants of twoMycetohabitansstrains revealed that the single Btl in oneMycetohabitans endofungorumstrain enhances fungal membrane stress tolerance, while others in aMycetohabitans rhizoxinicastrain promote bacterial colonization of the fungus. The phenotypic diversity underscores the need to assess the sequence diversity and, given that sequence diversity translates to DNA targeting specificity, the functional diversity of Btl proteins. Using a dual approach to maximize capture of Btl protein sequences for our analysis, we sequenced and assembled nineMycetohabitansspp. genomes using long-read PacBio technology and also mined available short-read Illumina fungal–bacterial metagenomes. We show thatbtlgenes are present across diverseMycetohabitansstrains from Mucoromycota fungal hosts yet vary in sequences and predicted DNA binding specificity. Phylogenetic analysis revealed distinct clades of Btl proteins and suggested thatMycetohabitansmight contain more species than previously recognized. Within our data set, Btl proteins were more conserved acrossM. rhizoxinicastrains than acrossM. endofungorum, but there was also evidence of greater overall strain diversity within the latter clade. Overall, the results suggest that Btl proteins contribute to bacterial–fungal symbioses in myriad ways.

Genetics & Heredity↗
DOE OSTI2024

Genomic surveillance indicates clonal replacement of hypervirulent Klebsiella pneumoniae ST881 and ST29 lineage strains in vivo

The emergence of hypervirulent Klebsiella pneumoniae (hvKp) poses a significant public health threat, particularly regarding its carriage in the healthy population. However, the genomic epidemiological characteristics and population dynamics of hvKp within a single patient across distinct infection episodes remain largely unknown. This study aimed to investigate the clonal replacement of hvKp K2-ST881 and K54-ST29 lineage strains in a single patient experiencing multiple-site infections during two independent episodes. Two strains, designated EDhvKp-1 and EDhvKp-2, were obtained from blood and cerebrospinal fluid during the first admission, and the strain isolated from blood on the second admission was named EDhvKp-3. Whole-genome sequencing, utilizing both short-read Illumina and long-read Oxford Nanopore platforms, was conducted. In silico multilocus sequence typing (MLST), identification of antimicrobial resistance and virulence genes, and the phylogenetic relationship between our strains and other K. pneumoniae ST881 and ST29 genomes retrieved from the public database were performed. Virulence potentials were assessed through a mouse lethality assay. Our study indicated that the strains were highly susceptible to multiple antimicrobial agents. Plasmid sequence analysis confirmed that both virulence plasmids, pEDhvKp-1 (166,008 bp) and pEDhvKp-3 (210,948 bp), belonged to IncFIB type. Multiple virulence genes, including rmpA , rmpA2 , rmpC , rmpD , iroBCDN , iucABCD , and iutA , were identified. EDhvKp-1 and EDhvKp-2 showed the closest relationship to strain 502 (differing by 51 SNPs), while EDhvKp-3 exhibited 69 SNPs differences compared to strain TAKPN-1, which all recovered from Chinese patients in 2020. In the mouse infection experiment, both ST881 EDhvKp-1 and ST29 EDhvKp-3 displayed similar virulence traits, causing 90 and 100% of the mice to die within 72 h after intraperitoneal infection, respectively. Our study expands the spectrum of hvKp lineages and highlights genomic alterations associated with clonal switching between two distinct lineages of hvKP that successively replaced each other in vivo . The development of novel strategies for the surveillance, diagnosis, and treatment of high-risk hvKp is urgently needed.

Liu, Ning↗
DOE OSTI2021

Impact of short-read sequencing on the misassembly of a plant genome

Abstract Background Availability of plant genome sequences has led to significant advances. However, with few exceptions, the great majority of existing genome assemblies are derived from short read sequencing technologies with highly uneven read coverages indicative of sequencing and assembly issues that could significantly impact any downstream analysis of plant genomes. In tomato for example, 0.6% (5.1 Mb) and 9.7% (79.6 Mb) of short-read based assembly had significantly higher and lower coverage compared to background, respectively. Results To understand what the causes may be for such uneven coverage, we first established machine learning models capable of predicting genomic regions with variable coverages and found that high coverage regions tend to have higher simple sequence repeat and tandem gene densities compared to background regions. To determine if the high coverage regions were misassembled, we examined a recently available tomato long-read based assembly and found that 27.8% (1.41 Mb) of high coverage regions were potentially misassembled of duplicate sequences, compared to 1.4% in background regions. In addition, using a predictive model that can distinguish correctly and incorrectly assembled high coverage regions, we found that misassembled, high coverage regions tend to be flanked by simple sequence repeats, pseudogenes, and transposon elements. Conclusions Our study provides insights on the causes of variable coverage regions and a quantitative assessment of factors contributing to plant genome misassembly when using short reads and the generality of these causes and factors should be tested further in other species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Discovery of additional ancient genome duplications in yeasts

Whole-genome duplication (WGD) has had profound macroevolutionary impacts on diverse lineages, preceding adaptive radiations in vertebrates, teleost fish, and angiosperms. In contrast to the many known ancient WGDs in animals, and especially plants, we are aware of evidence for only four WGDs in fungi. The oldest of these occurred ∼100 million years ago (mya) and is shared by ∼60 extant Saccharomycetales species, including the baker’s yeast Saccharomyces cerevisiae. Notably, this is the only known ancient WGD event in the yeast subphylum Saccharomycotina. The dearth of ancient WGD events in fungi remains a mystery. Some studies have suggested that fungal lineages that experience chromosome and genome duplication quickly go extinct, leaving no trace in the genomic record, while others contend that the lack of known WGDs is due to an absence of data. Under the second hypothesis, additional sampling and deeper sequencing of fungal genomes should lead to the discovery of more WGD events. Coupling hundreds of recently published genomes from nearly every described Saccharomycotina species, with three additional long-read assemblies, we discovered three novel WGD events. Although the functions of retained duplicate genes originating from these events are broad, they bear similarities to the well-known WGD that occurred in the Saccharomycetales. In conclusion, our results suggest that WGD may be a more common evolutionary force in fungi than previously believed.

convergent evolution↗
DOE OSTI2022

Critical Assessment of Metagenome Interpretation: the second round of challenges

Abstract Evaluating metagenomic software is key for optimizing metagenome interpretation and focus of the Initiative for the Critical Assessment of Metagenome Interpretation (CAMI). The CAMI II challenge engaged the community to assess methods on realistic and complex datasets with long- and short-read sequences, created computationally from around 1,700 new and known genomes, as well as 600 new plasmids and viruses. Here we analyze 5,002 results by 76 program versions. Substantial improvements were seen in assembly, some due to long-read data. Related strains still were challenging for assembly and genome recovery through binning, as was assembly quality for the latter. Profilers markedly matured, with taxon profilers and binners excelling at higher bacterial ranks, but underperforming for viruses and Archaea. Clinical pathogen detection results revealed a need to improve reproducibility. Runtime and memory usage analyses identified efficient programs, including top performers with other metrics. The results identify challenges and guide researchers in selecting methods for analyses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Chromosome assembled and annotated genome sequence of Aspergillus flavus NRRL 3357

Abstract Aspergillus flavus is an opportunistic pathogen of crops, including peanuts and maize, and is the second leading cause of aspergillosis in immunocompromised patients. A. flavus is also a major producer of the mycotoxin, aflatoxin, a potent carcinogen, which results in significant crop losses annually. The A. flavus isolate NRRL 3357 was originally isolated from peanut and has been used as a model organism for understanding the regulation and production of secondary metabolites, such as aflatoxin. A draft genome of NRRL 3357 was previously constructed, enabling the development of molecular tools and for understanding population biology of this particular species. Here, we describe an updated, near complete, telomere-to-telomere assembly and re-annotation of the eight chromosomes of A. flavus NRRL 3357 genome, accomplished via long-read PacBio and Oxford Nanopore technologies combined with Illumina short-read sequencing. A total of 13,715 protein-coding genes were predicted. Using RNA-seq data, a significant improvement was achieved in predicted 5’ and 3’ untranslated regions, which were incorporated into the new gene models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A Genome-Scale Atlas Reveals Complex Interplay of Transcription and Translation in an Archaeon

The scale of post-transcriptional regulation and the implications of its interplay with other forms of regulation in environmental acclimation are underexplored for organisms of the domain Archaea. Here, we have investigated the scale of post-transcriptional regulation in the extremely halophilic archaeon Halobacterium salinarum NRC-1 by integrating the transcriptome-wide locations of transcript processing sites (TPSs) and SmAP1 binding, the genome-wide locations of antisense RNAs (asRNAs), and the consequences of RNase_2099C knockout on the differential expression of all genes. This integrated analysis has discovered that 54% of all protein-coding genes in the genome of this haloarchaeon are likely targeted by multiple mechanisms for putative post-transcriptional processing and regulation, with about 20% of genes likely being regulated by combinatorial schemes involving SmAP1, asRNAs, and RNase_2099C. Comparative analysis of mRNA levels (transcriptome sequencing [RNA-Seq]) and protein levels (sequential window acquisition of all theoretical fragment ion spectra mass spectrometry [SWATH-MS]) for 2,579 genes over four phases of batch culture growth in complex medium generated additional evidence for the conditional post-transcriptional regulation of 7% of all protein-coding genes. We demonstrate that post-transcriptional regulation may act to fine-tune specialized and rapid acclimation to stressful environments, e.g., as a switch to turn on gas vesicle biogenesis to promote vertical relocation under anoxic conditions and modulate the frequency of transposition by insertion sequence (IS) elements of the IS200/IS605, IS4, and ISH3 families. Findings from this study are provided as an atlas in a public Web resource (https://halodata.systemsbiology.net).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Genomic insights into the origin, domestication and diversification of Brassica juncea

Despite early domestication around 3000 BC, the evolutionary history of the ancient allotetraploid species Brassica juncea (L.) Czern & Coss remains uncertain. Here, we report a chromosome-scale de novo assembly of a yellow-seeded B. juncea genome by integrating long-read and short-read sequencing, optical mapping and Hi-C technologies. Nuclear and organelle phylogenies of 480 accessions worldwide supported that B. juncea is most likely a single origin in West Asia, 8,000–14,000 years ago, via natural interspecific hybridization. Subsequently, new crop types evolved through spontaneous gene mutations and introgressions along three independent routes of eastward expansion. Selective sweeps, genome-wide trait associations and tissue-specific RNA-sequencing analysis shed light on the domestication history of flowering time and seed weight, and on human selection for morphological diversification in this versatile species. Our data provide a comprehensive insight into the origin and domestication and a foundation for genomics-based breeding of B. juncea.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Reference-free structural variant detection in microbiomes via long-read co-assembly graphs

Motivation: The study of bacterial genome dynamics is vital for understanding the mechanisms underlying microbial adaptation, growth, and their impact on host phenotype. Structural variants (SVs), genomic alterations of 50 base pairs or more, play a pivotal role in driving evolutionary processes and maintaining genomic heterogeneity within bacterial populations. While SV detection in isolate genomes is relatively straightforward, metagenomes present broader challenges due to the absence of clear reference genomes and the presence of mixed strains. In response, our proposed method rhea, forgoes reference genomes and metagenome-assembled genomes (MAGs) by encompassing all metagenomic samples in a series (time or other metric) into a single co-assembly graph. The log fold change in graph coverage between successive samples is then calculated to call SVs that are thriving or declining. Results: We show rhea to outperform existing methods for SV and horizontal gene transfer (HGT) detection in two simulated mock metagenomes, particularly as the simulated reads diverge from reference genomes and an increase in strain diversity is incorporated. We additionally demonstrate use cases for rhea on series metagenomic data of environmental and fermented food microbiomes to detect specific sequence alterations between successive time and temperature samples, suggesting host advantage. Our approach leverages previous work in assembly graph structural and coverage patterns to provide versatility in studying SVs across diverse and poorly characterized microbial communities for more comprehensive insights into microbial gene flux.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Adaptive modification of antiviral defense systems in microbial community under Cr-induced stress

Background The prokaryotic antiviral defense systems are crucial for mediating prokaryote-virus interactions that influence microbiome functioning and evolutionary dynamics. Despite the prevalence and significance of prokaryotic antiviral defense systems, their responses to abiotic stress and ecological consequences remain poorly understood in soil ecosystems. We established microcosm systems with varying concentrations of hexavalent chromium (Cr(VI)) to investigate the adaptive modifications of prokaryotic antiviral defense systems under abiotic stress. Results Utilizing hybrid metagenomic assembly with long-read and short-read sequencing, we discovered that anti- viral defense systems were more diverse and prevalent in heavily polluted soils, which was corroborated by meta-analyses of public datasets from various heavy metal-contaminated sites. As the Cr(VI) concentration increased, prokaryotes with defense systems favoring prokaryote-virus mutualism gradually supplanted those with defense systems incurring high adaptive costs. Additionally, as Cr(VI) concentrations increased, enriched antiviral defense systems exhibited synchronization with microbial heavy metal resistance genes. Furthermore, the proportion of antiviral defense systems carried by mobile genetic elements (MGEs), including plasmids and viruses, increased by approximately 43% and 39%, respectively, with rising Cr concentrations. This trend is conducive to strengthening the dissemination and sharing of defense resources within microbial communities. Conclusions Overall, our study reveals the adaptive modification of prokaryotic antiviral defense systems in soil ecosystems under abiotic stress, as well as their positive contributions to establishing prokaryote-virus mutualism and the evolution of microbial heavy metal resistance. These findings advance our understanding of microbial adaptation in stressful environments and may inspire novel approaches for microbiome manipulation and bioremediation.

59 BASIC BIOLOGICAL SCIENCES↗