Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “evolutionary biology”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The phylogenetic roots of addiction: compulsive drug seeking, natural and drug-sensitive reward, and the acquisition of learned habits

Our rational faculties permit us humans to maximize the utility of our actions. We perform a fundamental type of cost-benefit analysis in which we frame a problem, assign values to the different paths, and then choose from among a set of available options, the course that promises the most favorable outcome. So why then does addiction appear to be so impervious to the associated costs, so unaffected by undesired consequences, and ultimately so resistant to cognitive oversight? The answer may likely be found in the fact that the drivers for compulsive drug seeking and drug taking are located in affective brain circuits, circuits that are structured and patterned by learning with repeated activation. These deep processes exhibit significant resistance to control by our cognitive faculties. The general consensus is that addiction arises from mechanisms that overvalue the magnitude of reward, discount the associated risks, and thereby bias individuals towards compulsive pursuit of addictive drugs. Behavioral disruption and dependence appear to arise at the intersection of a number of connected but separate phenomena: expectations for the occurrence of specific events, behaviors that seek encounters with them, the ability to notice and learn nonrandom, co-occurring conditions, the prediction and valuation of consequences, the forming of enduring memories, and the drivers of focused behavior through compulsion, habits, and acquired routines. It is important to recognize that each one of these individual faculties are present and well developed across the entire phylogenetic tree of bilateral metazoans. The goal of this special volume is dedicated to exploring the degree to which inherent elements can account for addiction and addiction-associated phenomena. In much of the literature on addiction, the underlying processes are often viewed as distinctly mammalian, arising, in part, from the strong cognitive capacities of this taxon. Some phenomena may even be regarded to exist only in primates, or even solely in humans. This supposition arises from the fact that studies are conducted almost exclusively in mammals and primates, while evolutionary antecedents of the behavior are rarely considered. A more comprehensive perspective that examines drug reward and reinforcement in a wider range of organisms demonstrates that many of the component traits are actually well developed across the greater metazoan lineage, and they may well predate the emergence of a mammalian clade by a wide margin. The collection of papers assembled here supports the notion that the capacity to associate cues and quences has not arisen in mammals de novo. Rather, the neural mechanisms for detecting contingencies and for predicting future outcomes are very deeply rooted across broad phylogenetic divisions. Our understanding of an ability to associate paired events has been enriched by work in invertebrate preparations in both classical and operant conditioning scenarios [Cook and Carew, 1986, 1989a, 1989b]. The ability to learn allows us to connect cues and behavioral actions to their associated consequences. Pavlovian conditioning enriches surrounding cues with predictive value. Outcomes with positive valence generate appetitive responses and approach to the associated cues, while those perceived as aversive bring cue avoidance and withdrawal. Humans are not the only life forms capable of such short- and long-term modulations of behavior

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Novel Toxin Biosynthetic Gene Cluster in Harmful Algal Bloom-Causing Heteroscytonema crispum : Insights into the Origins of Paralytic Shellfish Toxins

Caused by both eukaryotic dinoflagellates and prokaryotic cyanobacteria, harmful algal blooms are events of severe ecological, economic, and public health consequence, and their incidence has become more common of late. Despite coordinated research efforts to identify and characterize the genomes of harmful algal bloom-causing organisms, the genomic basis and evolutionary origins of paralytic shellfish toxins produced by harmful algal blooms remain at best incomplete. The paralytic shellfish toxin saxitoxin has an especially complex genomic architecture and enigmatic phylogenetic distribution, spanning dinoflagellates and multiple cyanobacterial genera. Using filtration and extraction techniques to target the desired cyanobacteria from nonaxenic culture, coupled with a combination of short- and long-read sequencing, we generated a reference-quality hybrid genome assembly for Heteroscytonema crispum UTEX LB 1556, a freshwater, paralytic shellfish toxin-producing cyanobacterium thought to have the largest known genome in its phylum. We report a complete, novel biosynthetic gene cluster for the paralytic shellfish toxin saxitoxin. Leveraging this biosynthetic gene cluster, we find support for the hypothesis that paralytic shellfish toxin production has appeared in divergent Cyanobacteria lineages through widespread and repeated horizontal gene transfer. This work demonstrates the utility of long-read sequencing and metagenomic assembly toward advancing our understanding of paralytic shellfish toxin biosynthetic gene cluster diversity and suggests a mechanism for the origin of paralytic shellfish toxin biosynthetic genes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Exploring the binding properties and activities of ancestral expansins

Bacterial expansins are non-lytic proteins capable of loosening cellulose networks, offering promising applications in agriculture, biotechnology, and material science. Their ability to disrupt noncovalent interactions in biopolymer matrices such as cellulose and chitin positions them as valuable tools for upgrading abundant natural materials. However, their industrial use remains limited due to their relatively low wall-loosening activity compared to plant expansins. To address this limitation, we applied Ancestral Sequence Resurrection (ASR) to reconstruct and characterize ancient variants of the Bacillus subtilis expansin BsEXLX1. ASR is a powerful evolutionary tool that enables the inference and synthesis of ancestral proteins, allowing researchers to explore functional traits that may have been lost over time. This approach not only provides insights into protein evolution but also facilitates the design of proteins with enhanced properties, such as improved substrate affinity or structural stability. In this study, we combined biochemical and biophysical assays to evaluate the activity and binding behavior of ancestral expansins. Our results reveal that ancestral variants exhibit increased cellulose affinity, reduced binding to acidic polysaccharides, and greater salt resistance. Furthermore, these traits enhance their wall-loosening activity and demonstrate the utility of ASR in engineering surface-active proteins for industrial applications, particularly in biomass processing and cellulose modification.

09 BIOMASS FUELS↗
DOE OSTI2024

Evolution of connectivity architecture in the Drosophila mushroom body

Brain evolution has primarily been studied at the macroscopic level by comparing the relative size of homologous brain centers between species. How neuronal circuits change at the cellular level over evolutionary time remains largely unanswered. Here, using a phylogenetically informed framework, we compare the olfactory circuits of three closely related Drosophila species that differ in their chemical ecology: the generalists Drosophila melanogaster and Drosophila simulans and Drosophila sechellia that specializes on ripe noni fruit. We examine a central part of the olfactory circuit that, to our knowledge, has not been investigated in these species—the connections between projection neurons and the Kenyon cells of the mushroom body—and identify species-specific connectivity patterns. We found that neurons encoding food odors connect more frequently with Kenyon cells, giving rise to species-specific biases in connectivity. These species-specific connectivity differences reflect two distinct neuronal phenotypes: in the number of projection neurons or in the number of presynaptic boutons formed by individual projection neurons. Finally, behavioral analyses suggest that such increased connectivity enhances learning performance in an associative task. Our study shows how fine-grained aspects of connectivity architecture in an associative brain center can change during evolution to reflect the chemical ecology of a species.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1992

On the early emergence of reverse transcription: theoretical basis and experimental evidence

Reverse transcriptase (RT) was first discovered as an essential catalyst in the biological cycle of retroviruses. However, in the past years evidence has accumulated showing that RTs are involved in a surprisingly large number of RNA-mediated transpositional events that include both viral and nonviral genetic entities. Although it is probable that some RT-bearing genetic elements like the different types of AIDS viruses and the mammalian LINE family have arisen in recent geological times, the possibility that reverse transcription first took place in the early Archean is supported by (1) the hypothesis that RNA preceded DNA as cellular genetic material; (2) the existence of homologous regions of the subunit tau of the E. coli DNA polymerase III with the simian immunodeficiency virus RT, the hepatitis B virus RT, and the beta' subunit of the E. coli RNA polymerase (McHenry et al. 1988); (3) the presence of several conserved motifs, including a 14-amino-acid segment that consists of an Asp-Asp pair flanked by hydrophobic amino acids, which are found in all RTs and in most cellular and viral RNA polymerases. However, whether extant RTs descend from the primitive polymerase involved in the RNA-to-DNA transition remains unproven. Substrate specificity of the AMV and HIV-1 RTs can be modified in the presence of Mn2+, a cation which allows them to add ribonucleotides to an oligo (dG) primer in a template-dependent reaction. This change in specificity is comparable to that observed under similar conditions in other nucleic acid polymerases. This experimentally induced change in RT substrate specificity may explain previous observations on the misincorporation of ribonucleotides by the Maloney murine sarcoma virus RT in the minus and plus DNA of this retrovirus (Chen and Temin 1980). Our results also suggest that HIV-infected macrophages and T-cell cells may contain mixed polynucleotides containing both ribo- and deoxyribonucleotides. The evolutionary significance of these changes in substrate specificities of nucleic acid polymerases is also discussed.

Non-NASA Center↗
DOE OSTI2026

The landscape of regulatory element evolution in a C4 perennial grass

Gene regulatory evolution is a well-known source of phenotypic diversity and adaptive evolution. Although cis-regulatory elements (CREs) play a vital role in gene expression evolution, the molecular evolution of CREs remains mostly unknown due to the difficulty in identifying and characterizing these functional elements. Comparative genomic analyses of noncoding DNA can be leveraged to identify conserved noncoding sequences (CNS), many of which may harbor functional CREs conserved by purifying selection. However, purely computational inference of CREs from putative CNS can be erroneous due to the complex genomic architecture in plants. One promising experimental approach to identify CREs is by profiling accessible chromatin regions (ACRs) that are often associated with the location of CREs. In this study, we use comparative genomics along with the profiling of ACRs to study the molecular evolution of putative functional noncoding regulatory regions in Panicoid grasses. We identified sets of CNS that varied in relationship to the degree of evolutionary divergence among the studied taxa, including identifying core-Panicoid-CNS. We augmented this analysis by profiling ACRs in Panicum hallii ecotypes using ATAC-seq. ACRs had low SNP density at the summit, harbored a high frequency of core-Panicoid-CNS, and were enriched with expression QTL. These data help to annotate the P. hallii genome for putative functional elements and suggest that a large proportion of these ACRs are evolving under purifying selection. Turnover in CNS and ACR between ecotypes of P. hallii identifies a small set of putatively divergent CREs that may underlie differences in gene regulation between genotypes from inland and coastal habitats. In summary, we profiled ACRs in Panicoid grasses and integrated this data with our putative CNS prediction framework, which provides unique insight into patterns of polymorphism and divergence in CREs in C4 perennial grasses.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Large Scale MD to Predict Epitope Regions in HIV Env [Slides]

Highly dense carbohydrates located on the surface of the HIV Env protein play a key role in immune evasion. Such evolutionary adaptation hampers any attempt to obtain a full mechanistic understanding of the role played by the glycans in protecting the virus against an effective immune response. Moreover, and due to their chemical variability, an accurate molecular understanding of the so called “glycan shield” is still limited by the lack of effective resolution of state-of-the-art experimental technics. Here, we have used extensive computational modelling in order to fill this gap, addressing the presence of a large glycan variability as observed experimentally. Based on an automated pipeline, we were able to assemble, set-up and simulate via Molecular dynamics hundreds of different glycosylated Env variants at nearly atomic resolution, recapitulating the glycosylation distributions observed experimentally. Results from these simulations were subjected to machine learning and very accurate prediction of simulation derived glycan shielding areas of each glycan as a function of static sequence features. Such predictive models of per-glycan shielding, incorporating both glycan dynamics and heterogeneity, were used to develop a novel sequence-based glycan shield mapping strategy. Parallel to these studies, we also developed an accurate machine learning approach to predict glycan heterogeneity data using sequence features and found good prediction accuracy.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Complete biosynthesis of salicylic acid from phenylalanine in plants

Salicylic acid (SA) is a pivotal phytohormone for plant responses to biotic and abiotic stresses. Plants have evolved two pathways to produce SA: the isochorismate synthase and phenylalanine ammonia lyase (PAL) pathways. Whereas the isochorismate synthase pathway has been fully identified, the PAL pathway remains incomplete. Here we report the full characterization of the PAL pathway for SA biosynthesis via functional analysis of rice (Oryza sativa) SA-DEFICIENT GENE 1 (OSD1) to OSD4. The cinnamoyl-coenzyme A (CoA) ligase OSD1 catalyses the conversion of trans-cinnamic acid to cinnamoyl-CoA, which is subsequently transformed to benzoyl-CoA via the β-oxidative pathway in peroxisomes. The resulting benzoyl-CoA is further converted to benzyl benzoate by the peroxisomal benzoyltransferase OSD2. Benzyl benzoate is subsequently hydroxylated to benzyl salicylate by the endoplasmic reticulum membrane-resident cytochrome P450 OSD3, which is ultimately hydrolysed to salicylic acid by the cytoplasmic carboxylesterase OSD4. Evolutionary analyses reveal that the PAL pathway was first assembled before the divergence of gymnosperms and has been conserved in most seed plants. Activation of the PAL pathway in rice significantly enhances salicylic acid levels and plant immunity. Completion of the PAL pathway provides critical insights into the primary salicylic acid biosynthetic pathway across plant species and offers a precise target for modulating crop immunity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

What are grana in chloroplasts of vascular plants good for?

All plants and green algae contain stacked grana thylakoid membranes in their chloroplasts, underscoring an evolutionary pressure to maintain this unique structural feature. In addition, numerous studies have demonstrated that particular lateral and vertical dimensions of grana facilitate the function, regulation and repair of the photosynthetic machinery responsible for energy conversion. In this review, we present an updated overview of our understanding concerning the structure of grana thylakoids, the forces that contribute to their formation and their architectural dynamics. After establishing the structural foundation, we continue by reviewing the implications of grana formation on light harvesting, electron transport and protein maintenance in the thylakoid membranes of vascular plants. Regarding light harvesting, we discuss the implications of grana formation on energy spillover, macromolecular crowding, lateral versus vertical excitation energy transfer, and state transition. With respect to electron transport, we summarize recent findings regarding the functional connectivity of spatially separated photosystems facilitated by grana formation through small mobile electron carriers. We also explore how the dynamic responses of grana architecture regulate electron transport. Finally, we address how grana formation contributes to the organization of protein repair and biogenesis within thylakoid membranes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

OrthoPhyl—streamlining large-scale, orthology-based phylogenomic studies of bacteria at broad evolutionary scales

Abstract There are a staggering number of publicly available bacterial genome sequences (at writing, 2.0 million assemblies in NCBI's GenBank alone), and the deposition rate continues to increase. This wealth of data begs for phylogenetic analyses to place these sequences within an evolutionary context. A phylogenetic placement not only aids in taxonomic classification but informs the evolution of novel phenotypes, targets of selection, and horizontal gene transfer. Building trees from multi-gene codon alignments is a laborious task that requires bioinformatic expertise, rigorous curation of orthologs, and heavy computation. Compounding the problem is the lack of tools that can streamline these processes for building trees from large-scale genomic data. Here we present OrthoPhyl, which takes bacterial genome assemblies and reconstructs trees from whole genome codon alignments. The analysis pipeline can analyze an arbitrarily large number of input genomes (>1200 tested here) by identifying a diversity-spanning subset of assemblies and using these genomes to build gene models to infer orthologs in the full dataset. To illustrate the versatility of OrthoPhyl, we show three use cases: E. coli/Shigella, Brucella/Ochrobactrum and the order Rickettsiales. We compare trees generated with OrthoPhyl to trees generated with kSNP3 and GToTree along with published trees using alternative methods. We show that OrthoPhyl trees are consistent with other methods while incorporating more data, allowing for greater numbers of input genomes, and more flexibility of analysis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

REC protein family expansion by the emergence of a new signaling pathway

This report presents multi-genome evidence that REC protein family expansion occurs when the emergence of new pathways gives rise to functional discordance. Specificity between residues in REC domain containing response regulators with paired histidine kinases is under negative purifying selection, constrained by the presence of other bacterial two-component systems signaling cascades that share sequence and structural identity. Presuming that the two-component systems can evolve by neutral amino acid changes (neutral drift) when purifying evolutionary constraints are relaxed, how might the REC protein family expand by amino acid changes when these constraints remain intact? Using an unsupervised machine learning approach to observe the sequence landscape of REC domains across long phylogenetic distances, we find that within-gene recombination, a subcategory of gene conversion, switched the effector domain and, consequently, the regulatory context of a duplicated response regulator from transcriptional regulation by σ54 to that by σ70. We determined that the recombined response regulator diverged from its parent by episodic diversifying selection and neutral drift. Functional experiments of the parent of recombined response regulators in a model Pseudomonas putida KT2440 model system revealed that the parent and recombined response regulators sense and respond to different carboxylic acids. Finally, a residue-switching experiment using structural predictions and functional characterization suggests that the new residues in the recombined regulator could form a new interaction interface and mediate condition-specific phosphotransfer. Overall, our study finds that genetic perturbations can create conditions of functional discordance, whereby the REC protein family can evolve by episodic diversifying selection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Independent Evidence for the Preservation of Endogenous Bone Biochemistry in a Specimen of Tyrannosaurus rex

Biomolecules preserved in deep time have potential to shed light on major evolutionary questions, driving the search for new and more rigorous methods to detect them. Despite the increasing body of evidence from a wide variety of new, high resolution/high sensitivity analytical techniques, this research is commonly met with skepticism, as the long standing dogma persists that such preservation in very deep time (>1 Ma) is unlikely. The Late Cretaceous dinosaur Tyrannosaurus rex (MOR 1125) has been shown, through multiple biochemical studies, to preserve original bone chemistry. Here, we provide additional, independent support that deep time bimolecular preservation is possible. We use synchrotron X-ray fluorescence imaging (XRF) and X-ray absorption spectroscopy (XAS) to investigate a section from the femur of this dinosaur, and demonstrate preservation of elements (S, Ca, and Zn) associated with bone remodeling and redeposition. We then compare these data to the bone of an extant dinosaur (bird), as well as a second non-avian dinosaur, Tenontosaurus tilletti (OMNH 34784) that did not preserve any sign of original biochemistry. Our data indicate that MOR 1125 bone cortices have similar bone elemental distributions to that of an extant bird, which supports preservation of original endogenous chemistry in this specimen.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

A deep dilated convolutional residual network for predicting interchain contacts of protein homodimers

Abstract Motivation Deep learning has revolutionized protein tertiary structure prediction recently. The cutting-edge deep learning methods such as AlphaFold can predict high-accuracy tertiary structures for most individual protein chains. However, the accuracy of predicting quaternary structures of protein complexes consisting of multiple chains is still relatively low due to lack of advanced deep learning methods in the field. Because interchain residue–residue contacts can be used as distance restraints to guide quaternary structure modeling, here we develop a deep dilated convolutional residual network method (DRCon) to predict interchain residue–residue contacts in homodimers from residue–residue co-evolutionary signals derived from multiple sequence alignments of monomers, intrachain residue–residue contacts of monomers extracted from true/predicted tertiary structures or predicted by deep learning, and other sequence and structural features. Results Tested on three homodimer test datasets (Homo_std dataset, DeepHomo dataset and CASP-CAPRI dataset), the precision of DRCon for top L/5 interchain contact predictions (L: length of monomer in a homodimer) is 43.46%, 47.10% and 33.50% respectively at 6 Å contact threshold, which is substantially better than DeepHomo and DNCON2_inter and similar to Glinter. Moreover, our experiments demonstrate that using predicted tertiary structure or intrachain contacts of monomers in the unbound state as input, DRCon still performs well, even though its accuracy is lower than using true tertiary structures in the bound state are used as input. Finally, our case study shows that good interchain contact predictions can be used to build high-accuracy quaternary structure models of homodimers. Availability and implementation The source code of DRCon is available at https://github.com/jianlin-cheng/DRCon. The datasets are available at https://zenodo.org/record/5998532#.YgF70vXMKsB. Supplementary information Supplementary data are available at Bioinformatics online.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Cytochrome b5 diversity in green lineages preceded the evolution of syringyl lignin biosynthesis

Abstract Lignin production marked a milestone in vascular plant evolution, and the emergence of syringyl (S) lignin is lineage specific. S-lignin biosynthesis in angiosperms, mediated by ferulate 5-hydroxylase (F5H, CYP84A1), has been considered a recent evolutionary event. F5H uniquely requires the cytochrome b5 protein CB5D as an obligatory redox partner for catalysis. However, it remains unclear how CB5D functionality originated and whether it coevolved with F5H. We reveal here the ancient evolution of CB5D-type function supporting F5H-catalyzed S-lignin biosynthesis. CB5D emerged in charophyte algae, the closest relatives of land plants, and is conserved and proliferated in embryophytes, especially in angiosperms, suggesting functional diversification of the CB5 family before terrestrialization. A sequence motif containing acidic amino residues in Helix 5 of the CB5 heme-binding domain contributes to the retention of CB5D function in land plants but not in algae. Notably, CB5s in the S-lignin–producing lycophyte Selaginella lack these residues, resulting in no CB5D-type function. An independently evolved S-lignin biosynthetic F5H (CYP788A1) in Selaginella relies on NADPH-dependent cytochrome P450 reductase as sole redox partner, distinct from angiosperms. These results suggest that angiosperm F5Hs coopted the ancient CB5D, forming a modern cytochrome P450 monooxygenase system for aromatic ring meta-hydroxylation, enabling the reemergence of S-lignin biosynthesis in angiosperms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Systematic and scalable genome-wide essentiality mapping to identify nonessential genes in phages

Phages are one of the key ecological drivers of microbial community dynamics, function, and evolution. Despite their importance in bacterial ecology and evolutionary processes, phage genes are poorly characterized, hampering their usage in a variety of biotechnological applications. Methods to characterize such genes, even those critical to the phage life cycle, are labor intensive and are generally phage specific. Here, we develop a systematic gene essentiality mapping method scalable to new phage–host combinations that facilitate the identification of nonessential genes. As a proof of concept, we use an arrayed genome-wide CRISPR interference (CRISPRi) assay to map gene essentiality landscape in the canonical coliphages λ and P1. Results from a single panel of CRISPRi probes largely recapitulate the essential gene roster determined from decades of genetic analysis for lambda and provide new insights into essential and nonessential loci in P1. We present evidence of how CRISPRi polarity can lead to false positive gene essentiality assignments and recommend caution towards interpreting CRISPRi data on gene essentiality when applied to less studied phages. Finally, we show that we can engineer phages by inserting DNA barcodes into newly identified inessential regions, which will empower processes of identification, quantification, and tracking of phages in diverse applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Physiological potential and evolutionary trajectories of syntrophic sulfate-reducing bacterial partners of anaerobic methanotrophic archaea

Sulfate-coupled anaerobic oxidation of methane (AOM) is performed by multicellular consortia of anaerobic methanotrophic archaea (ANME) in obligate syntrophic partnership with sulfate-reducing bacteria (SRB). Diverse ANME and SRB clades co-associate but the physiological basis for their adaptation and diversification is not well understood. In this work, we used comparative metagenomics and phylogenetics to investigate the metabolic adaptation among the 4 main syntrophic SRB clades (HotSeep-1, Seep-SRB2, Seep-SRB1a, and Seep-SRB1g) and identified features associated with their syntrophic lifestyle that distinguish them from their non-syntrophic evolutionary neighbors in the phylum Desulfobacterota. We show that the protein complexes involved in direct interspecies electron transfer (DIET) from ANME to the SRB outer membrane are conserved between the syntrophic lineages. In contrast, the proteins involved in electron transfer within the SRB inner membrane differ between clades, indicative of convergent evolution in the adaptation to a syntrophic lifestyle. Our analysis suggests that in most cases, this adaptation likely occurred after the acquisition of the DIET complexes in an ancestral clade and involve horizontal gene transfers within pathways for electron transfer (CbcBA) and biofilm formation (Pel). We also provide evidence for unique adaptations within syntrophic SRB clades, which vary depending on the archaeal partner. Among the most widespread syntrophic SRB, Seep-SRB1a, subclades that specifically partner ANME-2a are missing the cobalamin synthesis pathway, suggestive of nutritional dependency on its partner, while closely related Seep-SRB1a partners of ANME-2c lack nutritional auxotrophies. Our work provides insight into the features associated with DIET-based syntrophy and the adaptation of SRB towards it.

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 OSTI2023

IMG/PR: a database of plasmids from genomes and metagenomes with rich annotations and metadata

Plasmids are mobile genetic elements found in many clades of Archaea and Bacteria. They drive horizontal gene transfer, impacting ecological and evolutionary processes within microbial communities, and hold substantial importance in human health and biotechnology. To support plasmid research and provide scientists with data of an unprecedented diversity of plasmid sequences, we introduce the IMG/PR database, a new resource encompassing 699 973 plasmid sequences derived from genomes, metagenomes and metatranscriptomes. IMG/PR is the first database to provide data of plasmid that were systematically identified from diverse microbiome samples. IMG/PR plasmids are associated with rich metadata that includes geographical and ecosystem information, host taxonomy, similarity to other plasmids, functional annotation, presence of genes involved in conjugation and antibiotic resistance. The database offers diverse methods for exploring its extensive plasmid collection, enabling users to navigate plasmids through metadata-centric queries, plasmid comparisons and BLAST searches. The web interface for IMG/PR is accessible at https://img.jgi.doe.gov/pr. Plasmid metadata and sequences can be downloaded from https://genome.jgi.doe.gov/portal/IMG_PR.

59 BASIC BIOLOGICAL SCIENCES↗