Search NASA⌕ Search

SEARCH · Search NASA

Results for “NUCLEIC ACID”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A standardized quantitative analysis strategy for stable isotope probing metagenomics

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

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Structural and functional analyses of SARS-CoV-2 Nsp3 and its specific interactions with the 5’ UTR of the viral genome

ABSTRACT Non-structural protein 3 (Nsp3) is the largest open reading frame encoded in the SARS-CoV-2 genome, essential for the formation of double-membrane vesicles (DMV) wherein viral RNA replication occurs. We conducted an extensive structure-function analysis of Nsp3 and determined the crystal structures of the ubiquitin-like 1 (Ubl1), nucleic acid binding (NAB), β-coronavirus-specific marker (βSM) domains, and a sub-region of the Y domain of this protein. We show that the Ubl1, ADP-ribose phosphatase (ADRP), human SARS Unique (HSUD), NAB, and Y domains of Nsp3 bind the 5’ UTR of the viral genome and that the Ubl1 and Y domains possess affinity for recognition of this region, suggesting high specificity. The Ubl1-Nucleocapsid (N) protein complex binds the 5’ UTR with greater affinity than the individual proteins alone. Our results suggest that multiple domains of Nsp3, particularly Ubl1 and Y, shepherd the 5’ UTR of the viral genome during translocation through the DMV membrane, priming the Ubl1 domain to load the genome onto N protein. IMPORTANCE The largest protein encoded by the SARS-CoV-2 genome is Nsp3. In infected cells, this multi-domain protein forms a pore structure in the virus-induced double-membrane vesicles (DMV). We have incomplete data on Nsp3 molecular structure, and here, we describe crystal structures for multiple domains of Nsp3. It is thought that newly replicated viral RNA transits through the DMV pore; however, we possess incomplete data on which regions of Nsp3 actually interact with RNA. Here, we present data showing that five domains of Nsp3 interact with the 5’ UTR of the SARS-CoV-2 RNA, including the Y domain for which no function has ever been discovered. These data suggest that the pore structure plays an active role in recognizing the terminal end of the genome, transiting and loading the viral RNA onto the cytoplasmic nucleocapsid protein. These data help expand our knowledge of Nsp3 structure and function and the SARS-CoV-2 replication cycle.

Microbiology↗
DOE OSTI2018

CNT-Based Label-Free Electrochemical Sensing of Native DNA with Allele Single Nucleotide Polymorphism

A label-free high-performance electrochemical sensing of single nucleotide polymorphism in native genomic DNA samples from cancer tissues, in oncogene-specific complementary DNA and in amplicons has been carried out on bundles of few-walled carboxylated carbon nanotubes, surface of which was coated with electron-dense layer of probe DNA. High selectivity and sensitivity of the novel DNA-nanosensor are ensured by the following features. Firstly, a complementary-hybridization signal is enhanced by resonance between plasmons in carbon nanotubes and dipole oscillations of hydrated nucleic acids located in near-electrode Helmholtz double layer of the impedance sensor. Secondly, a decoration of carbon nanotubes by adhered metal atoms in a high-spin state leads to Fermi level stabilization for the carbon nanotube bundles.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2016

Permafrost Meta-Omics and Climate Change

Permanently frozen soil, or permafrost, covers a large portion of the Earth's terrestrial surface and represents a unique environment for cold-adapted microorganisms. As permafrost thaws, previously protected organic matter becomes available for microbial degradation. Microbes that decompose soil carbon produce carbon dioxide and other greenhouse gases, contributing substantially to climate change. Next-generation sequencing and other -omics technologies offer opportunities to discover the mechanisms by which microbial communities regulate the loss of carbon and the emission of greenhouse gases from thawing permafrost regions. Analysis of nucleic acids and proteins taken directly from permafrost-associated soils has provided new insights into microbial communities and their functions in Arctic environments that are increasingly impacted by climate change. In this article we review current information from various molecular -omics studies on permafrost microbial ecology and explore the relevance of these insights to our current understanding of the dynamics of permafrost loss due to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Vacuum Ultraviolet Photoionization of Complex Chemical Systems

Tunable vacuum ultraviolet (VUV) radiation coupled to mass spectrometry is applied to the study of complex chemical systems in this paper. The identification of novel reactive intermediates and radicals is revealed in flame, pulsed photolysis, and pyrolysis reactors, leading to the elucidation of spectroscopy, reaction mechanisms, and kinetics. Mass-resolved threshold photoelectron photoion coincidence measurements provide unprecedented access to vibrationally resolved spectra of free radicals present in high-temperature reactors. Photoionization measurements in water clusters, nucleic acid base dimers, and their complexes with water provide signatures of proton transfer in hydrogen-bonded and π-stacked systems. Experimental and theoretical methods to track ion–molecule reactions and fragmentation pathways in intermolecular and intramolecular hydrogen-bonded systems in sugars and alcohols are described. Photoionization of laser-ablated molecules, clusters, and their reaction products inform thermodynamics and spectroscopy that are relevant to astrochemistry and catalysis. Finally, new directions in coupling VUV radiation to interrogate complex chemical systems are discussed.

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

A Grand Canonical Monte Carlo simulation program for computing ion distributions around biomolecules in hard sphere solvents

The GIBS software program is a Grand Canonical Monte Carlo (GCMC) simulation program (written in C++) that can be used for 1) computing the excess chemical potential of ions and the mean activity coefficients of salts in homogeneous electrolyte solutions; and, 2) for computing the distribution of ions around fixed macromolecules such as, nucleic acids and proteins. The solvent can be represented as neutral hard spheres or as a dielectric continuum. The ions are represented as charged hard spheres that can interact via Coulomb, hard-sphere, or Lennard-Jones potentials. In addition to hard-sphere repulsions, the ions can also be made to interact with the solvent hard spheres via short-ranged attractive square-well potentials.

Source record↗
DOE OSTI2018

Sequedex V2

Sequedex is computer software to rapidly assign phylogeny and function to short sequences of nucleic acid or protein. In addition, Sequedex comprises functions required to generate underlying data modules as well as tools for analyzing assignments in a variety of contexts.

Cohn, Judith↗
DOE OSTI2022

G4Boost: a machine learning-based tool for quadruplex identification and stability prediction

Background: G-quadruplexes (G4s), formed within guanine-rich nucleic acids, are secondary structures involved in important biological processes. Although every G4 motif has the potential to form a stable G4 structure, not every G4 motif would, and accurate energy-based methods are needed to assess their structural stability. Here, we present a decision tree-based prediction tool, G4Boost, to identify G4 motifs and predict their secondary structure folding probability and thermodynamic stability based on their sequences, nucleotide compositions, and estimated structural topologies. Results: G4Boost predicted the quadruplex folding state with an accuracy greater then 93% and an F1-score of 0.96, and the folding energy with an RMSE of 4.28 and R 2 of 0.95 only by the means of sequence intrinsic feature. G4Boost was successfully applied and validated to predict the stability of experimentally-determined G4 structures, including for plants and humans. Conclusion: G4Boost outperformed the three machine-learning based prediction tools, DeepG4, Quadron, and G4RNA Screener, in terms of both accuracy and F1-score, and can be highly useful for G4 prediction to understand gene regulation across species including plants and humans.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Carbon nanotube biocompatibility in plants is determined by their surface chemistry

Abstract Background Agriculture faces significant global challenges including climate change and an increasing food demand due to a growing population. Addressing these challenges will require the adoption of transformative innovations into biotechnology practice, such as nanotechnology. Recently, nanomaterials have emerged as unmatched tools for their use as biosensors, or as biomolecule delivery vehicles. Despite their increasingly prolific use, plant-nanomaterial interactions remain poorly characterized, drawing into question the breadth of their utility and their broader environmental compatibility. Results Herein, we characterize the response of Arabidopsis thaliana to single walled carbon nanotube (SWNT) exposure with two different surface chemistries commonly used for biosensing and nucleic acid delivery: oligonucleotide adsorbed-pristine SWNTs, and polyethyleneimine-SWNTs loaded with plasmid DNA (PEI-SWNTs), both introduced by leaf infiltration. We observed that pristine SWNTs elicit a mild stress response almost undistinguishable from the infiltration process, indicating that these nanomaterials are well-tolerated by the plant. However, PEI-SWNTs induce a much larger transcriptional reprogramming that involves stress, immunity, and senescence responses. PEI-SWNT-induced transcriptional profile is very similar to that of mutant plants displaying a constitutive immune response or treated with stress-priming agrochemicals. We selected molecular markers from our transcriptomic analysis and identified PEI as the main cause of this adverse reaction. We show that PEI-SWNT response is concentration-dependent and, when persistent over time, leads to cell death. We probed a panel of PEI variant-functionalized SWNTs across two plant species and identified biocompatible SWNT surface functionalizations. Conclusions While SWNTs themselves are well tolerated by plants, SWNTs surface-functionalized with positively charged polymers become toxic and produce cell death. We use molecular markers to identify more biocompatible SWNT formulations. Our results highlight the importance of nanoparticle surface chemistry on their biocompatibility and will facilitate the use of functionalized nanomaterials for agricultural improvement. Graphical Abstract

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

The waiting time problem in a model hominin population

Functional information is normally communicated using specific, context-dependent strings of symbolic characters. This is true within the human realm (texts and computer programs), and also within the biological realm (nucleic acids and proteins). In biology, strings of nucleotides encode much of the information within living cells. How do such information-bearing nucleotide strings arise and become established? This paper uses comprehensive numerical simulation to understand what types of nucleotide strings can realistically be established via the mutation/selection process, given a reasonable timeframe. The program Mendel’s Accountant realistically simulates the mutation/selection process, and was modified so that a starting string of nucleotides could be specified, and a corresponding target string of nucleotides could be specified. We simulated a classic pre-human hominin population of at least 10,000 individuals, with a generation time of 20 years, and with very strong selection (50 % selective elimination). Random point mutations were generated within the starting string. Whenever an instance of the target string arose, all individuals carrying the target string were assigned a specified reproductive advantage. When natural selection had successfully amplified an instance of the target string to the point of fixation, the experiment was halted, and the waiting time statistics were tabulated. Using this methodology we tested the effect of mutation rate, string length, fitness benefit, and population size on waiting time to fixation. Biologically realistic numerical simulations revealed that a population of this type required inordinately long waiting times to establish even the shortest nucleotide strings. To establish a string of two nucleotides required on average 84 million years. To establish a string of five nucleotides required on average 2 billion years. We found that waiting times were reduced by higher mutation rates, stronger fitness benefits, and larger population sizes. However, even using the most generous feasible parameters settings, the waiting time required to establish any specific nucleotide string within this type of population was consistently prohibitive. We show that the waiting time problem is a significant constraint on the macroevolution of the classic hominin population. Routine establishment of specific beneficial strings of two or more nucleotides becomes very problematic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

A comparison of methods used to unveil the genetic and metabolic pool in the built environment

BACKGROUND:A majority of indoor residential microbes originate from humans, pets, and outdoor air and are not adapted to the built environment (BE). Consequently, a large portion of the microbes identified by DNA-based methods are either dead or metabolically inactive. Although many exceptions have been noted, the ribosomal RNA fraction of the sample is more likely to represent either viable or metabolically active cells. We examined methodological variations in sample processing using a defined, mock BE microbial community to better understand the scope of technique-based vs. biological-based differences in both ribosomal transcript (rRNA) and gene (DNA) sequence community analysis. Based on in vitro tests, a protocol was adopted for the analysis of the genetic and metabolic pool (DNA vs. rRNA) of air and surface microbiomes within a residential setting. RESULTS:We observed differences in DNA/RNA co-extraction efficiency for individual microbes, but overall, a greater recovery of rRNA using FastPrep (> 50%). Samples stored with various preservation methods at - 80°C experienced a rapid decline in nucleic acid recovery starting within the first week, although post-extraction rRNA had no significant degradation when treated with RNAStable. We recommend that co-extraction samples be processed as quickly as possible after collection. The in vivo analysis revealed significant differences in the two components (genetic and metabolic pool) in terms of taxonomy, community structure, and microbial association networks. Rare taxa present in the genetic pool showed higher metabolic potential (RNA:DNA ratio), whereas commonly detected taxa of outdoor origins based on DNA sequencing, especially taxa of the Sphingomonadales order, were present in lower relative abundances in the viable community. CONCLUSIONS:Although methodological variations in sample preparations are high, large differences between the DNA and RNA fractions of the total microbial community demonstrate that direct examination of rRNA isolated from a residential BE microbiome has the potential to identify the more likely viable or active portion of the microbial community. In an environment that has primarily dead and metabolically inactive cells, we suggest that the rRNA fraction of BE samples is capable of providing a more ecologically relevant insight into the factors that drive indoor microbial community dynamics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

HT-SIP: a semi-automated stable isotope probing pipeline identifies cross-kingdom interactions in the hyphosphere of arbuscular mycorrhizal fungi

Abstract Background Linking the identity of wild microbes with their ecophysiological traits and environmental functions is a key ambition for microbial ecologists. Of many techniques that strive for this goal, Stable-isotope probing—SIP—remains among the most comprehensive for studying whole microbial communities in situ. In DNA-SIP, actively growing microorganisms that take up an isotopically heavy substrate build heavier DNA, which can be partitioned by density into multiple fractions and sequenced. However, SIP is relatively low throughput and requires significant hands-on labor. We designed and tested a semi-automated, high-throughput SIP (HT-SIP) pipeline to support well-replicated, temporally resolved amplicon and metagenomics experiments. We applied this pipeline to a soil microhabitat with significant ecological importance—the hyphosphere zone surrounding arbuscular mycorrhizal fungal (AMF) hyphae. AMF form symbiotic relationships with most plant species and play key roles in terrestrial nutrient and carbon cycling. Results Our HT-SIP pipeline for fractionation, cleanup, and nucleic acid quantification of density gradients requires one-sixth of the hands-on labor compared to manual SIP and allows 16 samples to be processed simultaneously. Automated density fractionation increased the reproducibility of SIP gradients compared to manual fractionation, and we show adding a non-ionic detergent to the gradient buffer improved SIP DNA recovery. We applied HT-SIP to 13 C-AMF hyphosphere DNA from a 13 CO 2 plant labeling study and created metagenome-assembled genomes (MAGs) using high-resolution SIP metagenomics (14 metagenomes per gradient). SIP confirmed the AMF Rhizophagus intraradices and associated MAGs were highly enriched (10–33 atom% 13 C), even though the soils’ overall enrichment was low (1.8 atom% 13 C). We assembled 212 13 C-hyphosphere MAGs; the hyphosphere taxa that assimilated the most AMF-derived 13 C were from the phyla Myxococcota, Fibrobacterota, Verrucomicrobiota, and the ammonia-oxidizing archaeon genus Nitrososphaera . Conclusions Our semi-automated HT-SIP approach decreases operator time and improves reproducibility by targeting the most labor-intensive steps of SIP—fraction collection and cleanup. We illustrate this approach in a unique and understudied soil microhabitat—generating MAGs of actively growing microbes living in the AMF hyphosphere (without plant roots). The MAGs’ phylogenetic composition and gene content suggest predation, decomposition, and ammonia oxidation may be key processes in hyphosphere nutrient cycling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Deep-ultraviolet ptychographic pocket-scope (DART): mesoscale lensless molecular imaging with label-free spectroscopic contrast

The mesoscale characterization of biological specimens has traditionally required compromises between resolution, field-of-view, depth-of-field, and molecular specificity, with most approaches relying on external labels. Here we present the Deep-ultrAviolet ptychogRaphic pockeT-scope (DART), a handheld platform that transforms label-free molecular imaging through intrinsic deep-ultraviolet spectroscopic contrast. By leveraging biomolecules’ natural absorption fingerprints and combining them with lensless ptychographic microscopy, DART resolves down to 308-nm linewidths across centimeter-scale areas while maintaining millimeter-scale depth-of-field. The system’s virtual error-bin methodology effectively eliminates artifacts from limited temporal coherence and other optical imperfections, enabling high-fidelity molecular imaging without lenses. Through differential spectroscopic imaging at deep-ultraviolet wavelengths, DART quantitatively maps nucleic acid and protein distributions with femtogram sensitivity, providing an intrinsic basis for explainable virtual staining. We demonstrate DART’s capabilities through imaging of tissue sections, cytopathology specimens, blood cells, and neural populations, revealing detailed molecular contrast without external labels. The combination of high-resolution molecular mapping and broad mesoscale imaging in a portable platform opens new possibilities from rapid clinical diagnostics, tissue analysis, to biological characterization in space exploration.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Fluorogenic aptamers resolve the flexibility of RNA junctions using orientation-dependent FRET

To further understand the transcriptome, new tools capable of measuring folding, interactions, and localization of RNA are needed. Although Förster resonance energy transfer (FRET) is an angle- and distance-dependent phenomenon, the majority of FRET measurements have been used to report distances, by assuming rotationally averaged donor–acceptor pairs. Angle-dependent FRET measurements have proven challenging for nucleic acids due to the difficulties in incorporating fluorophores rigidly into local substructures in a biocompatible manner. Fluorescence turn-on RNA aptamers are genetically encodable tags that appear to rigidly confine their cognate fluorophores, and thus have the potential to report angular-resolved FRET. In this work, we use the fluorescent aptamers Broccoli and Mango-III as donor and acceptor, respectively, to measure the angular dependence of FRET. Joining the two fluorescent aptamers by a helix of variable length allowed systematic rotation of the acceptor fluorophore relative to the donor. FRET oscillated in a sinusoidal manner as a function of helix length, consistent with simulated data generated from models of oriented fluorophores separated by an inflexible helix. Analysis of the orientation dependence of FRET allowed us to demonstrate structural rigidification of the NiCo riboswitch upon transition metal-ion binding. This application of fluorescence turn-on aptamers opens the way to improved structural interpretation of ensemble and single-molecule FRET measurements of RNA.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Structure of 3'-PO 4 /5'-OH RNA ligase RtcB in complex with a 5'-OH oligonucleotide

RtcB enzymes comprise a widely distributed family of manganese- and GTP-dependent RNA repair enzymes that join 2',3'-cyclic phosphate ends to 5'-OH ends via RtcB-(histidinyl-N)–GMP, RNA 3'-phosphate, and RNA 3 'pp 5 'G intermediates. RtcB can ligate either 5'-OH RNA or 5'-OH DNA strands in vitro. The nucleic acid contacts of RtcB are uncharted. In this work, we report a 2.7 Å crystal structure of Pyrococcus horikoshii RtcB in complex with a 6-mer 5'-OH DNA oligonucleotide HO A 1 pT 2 pG 3 pT 4 pC 5 pC 6 , which reveals enzymic contacts of Asn202 to the terminal 5'-OH nucleophile; Arg238 to the A1pT2 and T2pG3 phosphates; Arg190 and Gln194 to the T 2 pG 3 phosphate; and an Arg190 π-cation interaction with the G 3 nucleobase. The structural insights affirm functional studies of E. coli RtcB that implicated the conserved counterpart of Arg238 in engagement of the 5'-OH strand for ligation. The essential active site Cys98 that coordinates two manganese ions is oxidized to cysteine sulfonic acid in our structure, raising the prospect that RtcB activity might be sensitive to modulation during oxidative stress.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Structures of RNA phosphotransferase Tpt1 reveal distinct binding modes for an RNA 2′-PO 4 splice junction versus a 5′-PO 4 mononucleotide

Tpt1 is a widely distributed enzyme that removes an internal RNA 2′-phosphate by transfer to NAD + , via a two-step reaction in which: (i) the RNA 2′-PO 4 attacks NAD + to form an RNA-2′-phospho-(ADP-ribose) intermediate and expel nicotinamide; and (ii) the ADP-ribose O2″ attacks the RNA 2′-phosphodiester to form 2′-OH RNA and ADP-ribose-1″,2″-cyclic phosphate products. Tpt1 can also execute a single-step ADP-ribosyltransferase reaction at a 5′-monophosphate nucleic acid terminus that installs a 5′-phospho-ADP-ribose cap structure. Here we present crystal structures of Tpt1 bound to an RNA containing an internal 2′-PO 4 mark (the substrate for the canonical Tpt1 pathway) and in a complex with 5′-AMP. We find that Tpt1 has distinct binding modes, whereby the RNA 2′-PO 4 and the AMP 5′-PO 4 are engaged by the same set of active site amino acids, but the 2′-PO 4 nucleoside and the 5′-nucleoside occupy different sites on the enzyme.

Biochemistry & Molecular Biology↗
DOE OSTI2021

Data from PNA-Assisted DNAzymes to Cleave Double-Stranded DNA for Genetic Engineering with High Sequence Fidelity

DNAzymes have been widely used in many sensing and imaging applications but have rarely been used for genetic engineering since their discovery in 1994, because their substrate scope is mostly limited to single-stranded DNA or RNA, whereas genetic information is stored mostly in double-stranded DNA (dsDNA). To overcome this major limitation, we herein report peptide nucleic acid (PNA)-assisted double-stranded DNA nicking by DNAzymes (PANDA) as the first example to expand DNAzyme activity toward dsDNA. We show that PANDA is programmable in efficiently nicking or causing double strand breaks on target dsDNA, which mimics protein nucleases and can act as restriction enzymes in molecular cloning. In addition to being much smaller than protein enzymes, PANDA has a higher sequence fidelity compared with CRISPR/Cas under the condition we tested, demonstrating its potential as a novel alternative tool for genetic engineering and other biochemical applications.

Conversion↗
DOE OSTI2020

Kraken2 Metagenomic Virus Database

The Database: Kraken2 [1] database built from a classification tree containing over 700k metagenomic viruses from JGI IMG/VR [2]. (1) Wood, D. E., Lu, J., and Langmead, B. (2019). Improved metagenomic analysis with Kraken 2. Genome Biol., 20(1), 1–13. doi: 10.1186/s13059-019-1891-0 (2) Paez-Espino D, Chen I-MA, Palaniappan K, Ratner A, Chu K, Szeto E, et al. IMG/VR: a database of cultured and uncultured DNA Viruses and retroviruses. Nucleic Acids Res. 2017;45:D457-65. For Paper: Title: A k-mer based approach for virus classification in metatranscriptomic and metagenomic samples identifies viral associations in the Populus phytobiome and autism brains Abstract Background Viruses are an underrepresented taxa in the study and identification of microbiome constituents; however, they play an important role in health, microbiome regulation, and transfer of genetic material. Only a few thousand viruses have been isolated, sequenced, and assigned a taxonomy, which further limits the ability to identify and quantify viruses in the microbiome. Additionally, the vast diversity of viruses represents a challenge for classification, not only in constructing a viral taxonomy, but also in identifying similarities between a virus' genotype and its phenotype. However, the diversity of viral sequences can be leveraged to classify their sequences in metagenomic and metatranscriptomic samples. Methods To identify viruses in transcriptomic and genomic samples, we developed a dynamic programming algorithm for creating a classification tree out of 715,672 metagenome viruses. To create the classification tree, we clustered proportional similarity scores generated from the k-mer profiles of each of the metagenome viruses. We then integrated the viral classification tree with the NCBI taxonomy for use with ParaKraken, a metagenomic/transcriptomic classifier. Results To illustrate the breadth of our utility for classifying viruses with ParaKraken, we analyzed data from a plant metagenome study identifying the differences between two Populus genotypes in three different compartments and on a human metatranscriptome study identifying the differences between Autism Spectrum Disorder patients and controls in post mortem brain biopsies. In the Populus study, we identified genotype and compartment specific viral signatures, while in the Autism study we identified a significant increased abundance of eight viral sequences in Autism brain biopsies. Conclusion Viruses represent an important aspect of the microbiome. The ability to classify viruses represents the first step in being able to better understand their role in the microbiome. The viral classification method presented here allows for more complete identification of viral sequences for use in identifying associations between viruses and the host and viruses and other microbiome members. Acknowledgements and Funding This research used resources of the Oak Ridge Leadership Computing Facility at the Oak Ridge National Laboratory, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725. This research was also supported by the Plant-Microbe Interfaces Scientific Focus Area in the Genomic Science Program, the Office of Biological and Environmental Research (BER) in the U.S. Department of Energy Office of Science, and by the Department of Energy, Laboratory Directed Research and Development funding (ProjectID 8321), at the Oak Ridge National Laboratory. Oak Ridge National Laboratory is managed by UT-Battelle, LLC, for the US DOE under contract DE-AC05-00OR22725. This research used resources of the Compute and Data Environment for Science (CADES).

59 BASIC BIOLOGICAL SCIENCES↗