Search NASA⌕ Search

SEARCH · Search NASA

Results for “sequencing test”

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

Anatomy of a Summertime Convective Event over the Arabian Region

This study investigates the structure and evolution of a summertime convective event that occurred on 14 July 2015 over the Arabian region. We use the WRF Model with 1-km horizontal grid spacing and test three PBL parameterizations: the Mellor–Yamada–Nakanishi–Niino (MYNN) scheme; the Asymmetrical Convective Model, version 2, (ACM2) scheme; and the quasi-normal scale-elimination (QNSE) scheme. Convection initiates near the Al Hajar Mountains of northern Oman at around 1100 local time (LT; 0700 UTC) and propagates northwestward. A nonorographic convective band along the west coast of the United Arab Emirates (UAE) develops after 1500 LT as a result of the convergence of cold pools with the sea breeze from the Arabian Gulf. The model simulation employing the QNSE scheme simulates the convection initiation and propagation well. Although the MYNN and ACM2 simulations show convective initiation near the Al Hajar Mountains, they fail to simulate the development of the convective band along the UAE west coast. The MYNN run simulates colder near-surface temperatures and a weaker sea breeze, whereas the ACM2 run simulates a stronger sea breeze but a drier lower troposphere. Sensitivity simulations using horizontal grid spacings of 9 and 3 km show that lower-resolution runs develop broader convective structures and weaker cold pools and horizontal wind divergence, affecting the development of convection along the west coast of the UAE. The 1-km run using the QNSE PBL scheme realistically captures the sequence of events that leads to the moist convection over the UAE and adjacent mountains.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A Perspective on Biomass-Derived Biofuels: From Catalyst Design Principles to Fuel Properties

The hazards to health and the environment associated with the transportation sector include smog, particulate matter, and greenhouse gas emissions. Conversion of lignocellulosic biomass into biofuels has the potential to provide significant amounts of infrastructure-compatible liquid transportation fuels that reduce those hazardous materials. However, the development of these technologies is inefficient, due to: (i) the lack of a priori fuel property consideration, (ii) poor shared vocabulary between process chemists and fuel engineers, and (iii) modern and future engines operating outside the range of traditional autoignition metrics such as octane or cetane numbers. In this perspective, we describe an approach where we follow a “fuel-property first” design methodology with a sequence of (i) identifying the desirable fuel properties for modern engines, (ii) defining molecules capable of delivering those properties, and (iii) designing catalysts and processes that can produce those molecules from a candidate feedstock in a specific conversion process. Computational techniques need to be leveraged to minimize expenses and experimental efforts on low-promise options. This concept is illustrated with current research information available for biomass conversion to fuels via catalytic fast pyrolysis and hydrotreating; outstanding challenges and research tools necessary for a successful outcome are presented.

BIOMASS FUELS↗
DOE OSTI2024

How hydrophobicity, side chains, and salt affect the dimensions of disordered proteins

Abstract Despite the generally accepted role of the hydrophobic effect as the driving force for folding, many intrinsically disordered proteins (IDPs), including those with hydrophobic content typical of foldable proteins, behave nearly as self‐avoiding random walks (SARWs) under physiological conditions. Here, we tested how temperature and ionic conditions influence the dimensions of the N‐terminal domain of pertactin (PNt), an IDP with an amino acid composition typical of folded proteins. While PNt contracts somewhat with temperature, it nevertheless remains expanded over 10–58°C, with a Flory exponent, ν , >0.50. Both low and high ionic strength also produce contraction in PNt, but this contraction is mitigated by reducing charge segregation. With 46% glycine and low hydrophobicity, the reduced form of snow flea anti‐freeze protein (red‐sfAFP) is unaffected by temperature and ionic strength and persists as a near‐SARW, ν ~ 0.54, arguing that the thermal contraction of PNt is due to stronger interactions between hydrophobic side chains. Additionally, red‐sfAFP is a proxy for the polypeptide backbone, which has been thought to collapse in water. Increasing the glycine segregation in red‐sfAFP had minimal effect on ν . Water remained a good solvent even with 21 consecutive glycine residues ( ν > 0.5), and red‐sfAFP variants lacked stable backbone hydrogen bonds according to hydrogen exchange. Similarly, changing glycine segregation has little impact on ν in other glycine‐rich proteins. These findings underscore the generality that many disordered states can be expanded and unstructured, and that the hydrophobic effect alone is insufficient to drive significant chain collapse for typical protein sequences.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Pangenomes suggest ecological-evolutionary responses to experimental soil warming

ABSTRACT Below-ground carbon transformations that contribute to healthy soils represent a natural climate change mitigation, but newly acquired traits adaptive to climate stress may alter microbial feedback mechanisms. To better define microbial evolutionary responses to long-term climate warming, we study microorganisms from an ongoing in situ soil warming experiment where, for over three decades, temperate forest soils are continuously heated at 5°C above ambient. We hypothesize that across generations of chronic warming, genomic signatures within diverse bacterial lineages reflect adaptations related to growth and carbon utilization. From our bacterial culture collection isolated from experimental heated and control plots, we sequenced genomes representing dominant taxa sensitive to warming, including lineages of Actinobacteria, Alphaproteobacteria, and Betaproteobacteria. We investigated genomic attributes and functional gene content to identify signatures of adaptation. Comparative pangenomics revealed accessory gene clusters related to central metabolism, competition, and carbon substrate degradation, with few functional annotations explicitly associated with long-term warming. Trends in functional gene patterns suggest genomes from heated plots were relatively enriched in central carbohydrate and nitrogen metabolism pathways, while genomes from control plots were relatively enriched in amino acid and fatty acid metabolism pathways. We observed that genomes from heated plots had less codon bias, suggesting potential adaptive traits related to growth or growth efficiency. Codon usage bias varied for organisms with similar 16S rrn operon copy number, suggesting that these organisms experience different selective pressures on growth efficiency. Our work suggests the emergence of lineage-specific trends as well as common ecological-evolutionary microbial responses to climate change. IMPORTANCE Anthropogenic climate change threatens soil ecosystem health in part by altering below-ground carbon cycling carried out by microbes. Microbial evolutionary responses are often overshadowed by community-level ecological responses, but adaptive responses represent potential changes in traits and functional potential that may alter ecosystem function. We predict that microbes are adapting to climate change stressors like soil warming. To test this, we analyzed the genomes of bacteria from a soil warming experiment where soil plots have been experimentally heated 5°C above ambient for over 30 years. While genomic attributes were unchanged by long-term warming, we observed trends in functional gene content related to carbon and nitrogen usage and genomic indicators of growth efficiency. These responses may represent new parameters in how soil ecosystems feedback to the climate system.

Choudoir, Mallory J. (ORCID:0000000291175150)↗
DOE OSTI2020

Programmable Gene Knockdown in Diverse Bacteria Using Mobile‐CRISPRi

Abstract Facile bacterial genome sequencing has unlocked a veritable treasure trove of novel genes awaiting functional exploration. To make the most of this opportunity requires powerful genetic tools that can target all genes in diverse bacteria. CRISPR interference (CRISPRi) is a programmable gene‐knockdown tool that uses an RNA‐protein complex comprised of a single guide RNA (sgRNA) and a catalytically inactive Cas9 nuclease (dCas9) to sterically block transcription of target genes. We previously developed a suite of modular CRISPRi systems that transfer by conjugation and integrate into the genomes of diverse bacteria, which we call Mobile‐CRISPRi. Here, we provide detailed protocols for the modification and transfer of Mobile‐CRISPRi vectors for the purpose of knocking down target genes in bacteria of interest. We further discuss strategies for optimizing Mobile‐CRISPRi knockdown, transfer, and integration. We cover the following basic protocols: sgRNA design, cloning new sgRNA spacers into Mobile‐CRISPRi vectors, Tn 7 transfer of Mobile‐CRISPRi to Gram‐negative bacteria, and ICE Bs1 transfer of Mobile‐CRISPRi to Bacillales. © 2020 The Authors. Basic Protocol 1 : sgRNA design Basic Protocol 2 : Cloning of new sgRNA spacers into Mobile‐CRISPRi vectors Basic Protocol 3 : Tn 7 transfer of Mobile‐CRISPRi to Gram‐negative bacteria Basic Protocol 4 : ICE Bs1 transfer of Mobile‐CRISPRi to Bacillales Support Protocol 1 : Quantification of CRISPRi repression using fluorescent reporters Support Protocol 2 : Testing for gene essentiality using CRISPRi spot assays on plates Support Protocol 3 : Transformation of E. coli by electroporation Support Protocol 4 : Transformation of CaCl 2 ‐competent E. coli

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Multifunctional Catalysts for the Tandem Reactions of Oxygenates

Industrially-relevant catalytic reactions rarely consist of a simple sequence of elementary steps. Moreover, kinetic coupling of multiple reactions on a catalyst surface is highly desired for process intensification and improved energy efficiency for large scale chemical transformations. The shifting landscape of hydrocarbon chemical feedstocks in the US also motivates research on the selective conversion of more complex molecules. One desirable type of catalytic reaction is the reduction of carboxylic acids that are produced from biomass feedstocks to their corresponding alcohols. The proposed research explores the fundamental importance of hydrogen spillover on a multifunctional catalyst for carboxylic acid reduction with H 2 composed of metal particles coupled to metal oxide particles. Recent work has demonstrated the excellent performance of supported tungsten oxide clusters for carboxylic acid reduction, but only after they are promoted with a late transition metal such as palladium. Elucidating the active state of the catalyst and the associated reaction mechanism for acid reduction on that active state are the overall goals of the proposed project and successful completion will enable future design of efficient multifunctional catalysts. The critically important role of the metal promoter is hypothesized to be its ability to dissociate H2 and spillover atomic H to the support. Although hydrogen spillover is a well-recognized phenomenon in catalysis, its role in both catalyst activation and catalytic turnover are still unresolved. The study combined materials synthesis, characterization, reactivity testing, and molecular simulations, to explore the effect of hydrogen chemical potential on the formation of the active catalytic sites and on the steady state catalytic reduction of carboxylic acid. Varying the hydrogen chemical potential through modification of the gas conditions, support composition, and metal loading to modulated the structure and catalytic performance of the reducible metal oxide. Dual function catalysts containing supported Pd and WO x species co-located on a non-reducible carrier (silica) and a reducible carrier (titania) were synthesized and characterized by electron microscopy, temperature-programmed reduction, and chemisorption. Spectroscopic methods such as X-ray absorption and UV-vis were also used to evaluate the catalysts, which were used in the reduction of propionic acid to aldehyde and alcohol. Quantum chemical calculations, including ab initio phase diagrams provided molecular insights into the H spillover phenomenon.

09 BIOMASS FUELS↗
DOE OSTI2026

Finding the missing pieces: filling gaps that impede the translation of omics data into models

High-throughput omics technologies such as DNA sequencing have made the sequencing and computational assembly of microbial genomes recovered from the environment relatively routine. Computational inference of the protein products encoded by these genomes, and the associated biochemical functions, should enable the accurate prediction and modeling of microbial metabolism, organismal interactions, and ecosystem processes. However, a lack of scalable, probabilistic protein annotation tools limits the full potential of modeling for understanding the metabolism and biogeochemical cycles of microbial communities. Our approach to improve inference of protein annotations and metabolic models relied on learning from and emulating expert manual curation, leveraging software engineering and data science best practices to scale up the throughput and accuracy of annotations and metabolic model construction, building software to objectively evaluate different annotation strategies, and more closely linking the protein annotation and metabolic model inference process. Outcomes of this research include several improved or new computational tools, including DRAM (Distilled and Refined Annotation of Metabolism) for annotating microbial genomes with protein function and metabolic traits, CAMPER (Curated Annotations for Microbial Polyphenol Enzymes and Reactions) for annotating key polyphenol metabolisms, EC-Bench for comprehensive and unbiased benchmarking of annotation tools, and several apps available via the DOE Systems Biology Knowledgebase (KBase) for building genome-scale metabolic models. We demonstrate that these tools allow us to scalably annotate and understand thousands of genomes for microbial communities from a variety of systems and test cases, including rivers, thawing permafrost, and gut microbiomes. All of these computational tools are available as open-source software, with most broadly and easily accessible to the scientific community via KBase apps.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Mitochondrial Genomes of the United States Distribution of Gray Fox (Urocyon cinereoargenteus) Reveal a Major Phylogeographic Break at the Great Plains Suture Zone

We examined phylogeographic structure in gray fox (Urocyon cinereoargenteus) across the United States to identify the location of secondary contact zone(s) between eastern and western lineages and investigate the possibility of additional cryptic intraspecific divergences. We generated and analyzed complete mitochondrial genome sequence data from 75 samples and partial control region mitochondrial DNA sequences from 378 samples to investigate levels of genetic diversity and structure through population- and individual-based analyses including estimates of divergence (FST and SAMOVA), median joining networks, and phylogenies. We used complete mitochondrial genomes to infer phylogenetic relationships and date divergence times of major lineages of Urocyon in the United States. Despite broad-scale sampling, we did not recover additional major lineages of Urocyon within the United States, but identified a deep east-west split (~0.8 million years) with secondary contact at the Great Plains Suture Zone and confirmed the Channel Island fox (Urocyon littoralis) is nested within U. cinereoargenteus. Genetic diversity declined at northern latitudes in the eastern United States, a pattern concordant with post-glacial recolonization and range expansion. Beyond the east-west divergence, morphologically-based subspecies did not form monophyletic groups, though unique haplotypes were often geographically limited. Gray foxes in the United States displayed a deep, cryptic divergence suggesting taxonomic revision is needed. Secondary contact at a common phylogeographic break, the Great Plains Suture Zone, where environmental variables show a sharp cline, suggests ongoing evolutionary processes may reinforce this divergence. Follow-up study with nuclear markers should investigate whether hybridization is occurring along the suture zone and characterize contemporary population structure to help identify conservation units. Comparative work on other wide-ranging carnivores in the region should test whether similar evolutionary patterns and processes are occurring.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

High-order partitioned spectral deferred correction solvers for multiphysics problems

We present an arbitrarily high-order, conditionally stable, partitioned spectral deferred correction (SDC) method for solving multiphysics problems using a sequence of pre-existing single-physics solvers. This method extends the work in [1], [2], which used implicit-explicit Runge-Kutta methods (IMEX) to build high-order, partitioned multiphysics solvers. We consider a generic multiphysics problem modeled as a system of coupled ordinary differential equations (ODEs), coupled through coupling terms that can depend on the state of each subsystem; therefore the method applies to both a semi-discretized system of partial differential equations (PDEs) or problems naturally modeled as coupled systems of ODEs. The sufficient conditions to build arbitrarily high-order partitioned SDC schemes are derived. Based on these conditions, various of partitioned SDC schemes are designed. The stability of the first-order partitioned SDC scheme is analyzed in detail on a coupled, linear model problem. We show that the scheme is conditionally stable, and under conditions on the coupling strength, the scheme can be unconditionally stable. We demonstrate the performance of the proposed partitioned solvers on several classes of multiphysics problems with moderate coupling strength. They include a stiff linear system of ODEs, advection-diffusion-reaction systems, and fluid-structure interaction problems with both incompressible and compressible flows, where we verify the design order of the SDC schemes and study various stability properties. We also directly compare the accuracy, stability, and cost of the proposed partitioned SDC solver with the partitioned IMEX method in [1], [2] on this suite of test problems. The results suggest that the high-order partitioned SDC solvers are more robust than the partitioned IMEX solvers for the numerical examples considered in this work, while the IMEX methods require fewer implicit solves.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

The eROSITA Final Equatorial-Depth Survey (eFEDS)

In 2019, the eROSITA telescope on board the Russian-German satellite Spectrum-Roentgen-Gamma (SRG) began to perform a deep all-sky X-ray survey with the aim of identifying ~100 000 clusters and groups over the course of four years. As part of its performance verification phase, a ~140 deg 2 survey, called eROSITA Final Equatorial-Depth Survey (eFEDS), was performed. With a depth typical of the all-sky survey after four years, it allows tests of tools and methods as well as improved predictions for the all-sky survey. As part of this effort, a catalog of 542 X-ray selected galaxy group and cluster candidates was compiled. Here, in this paper, we present the optical follow-up, with the aim of providing redshifts and cluster confirmation for the full sample. Furthermore, we aim to provide additional information on the dynamical state, richness, and optical center of the clusters. Finally, we aim to evaluate the impact of optical cluster confirmation on the purity and completeness of the X-ray selected sample. We used optical imaging data from the Hyper Suprime-Cam Subaru Strategic Program and from the Legacy Survey to identify optical counterparts to the X-ray detected cluster candidates. We make use of the multi-component matched filter cluster confirmation tool (MCMF), as well as of the optical cluster finder CAMIRA to derive cluster redshifts and richnesses. MCMF provided the probabilities with which an optical structure would be a chance superposition with the X-ray candidate. These probabilities were used to identify the best optical counterpart as well as to confirm an X-ray candidate as a cluster. The impact of this confirmation process on catalog purity and completeness was estimated using optical to X-ray scaling relations as well as simulations. The resulting catalog was furthermore matched with public group and cluster catalogs. Optical estimators of the cluster dynamical state were constructed based on density maps of the red-sequence galaxies at the cluster redshift. By providing redshift estimates for all 542 candidates, we construct an optically confirmed sample of 477 clusters and groups with a residual contamination of 6%. Of these, 470 (98.5%) are confirmed using MCMF, and 7 systems are added through cross-matching with spectroscopic group catalogs. Using observable-to-observable scaling and the applied confirmation threshold, we predict that 8 ± 2 real systems have been excluded with the MCMF cut required to build this low-contamination sample. This number agrees well with the 7 systems found through cross-matching that were not confirmed with MCMF. The predicted redshift and mass distribution of this catalog agree well with simulations. Thus, we expect that these 477 systems include >99% of all true clusters in the candidate list. Using an MCMF-independent method, we confirm that the catalog contamination of the confirmed subsample is 6 ± 3%. Application of the same method to the full candidate list yields 17 ± 3%, consistent with estimates coming from the fraction of confirmed systems of ~17% and with expectations from simulations of ~20%. We also present a sample of merging cluster candidates based on the derived estimators of the cluster dynamical state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Profiling the BLAST bioinformatics application for load balancing on high-performance computing clusters

Abstract Background The Basic Local Alignment Search Tool (BLAST) is a suite of commonly used algorithms for identifying matches between biological sequences. The user supplies a database file and query file of sequences for BLAST to find identical sequences between the two. The typical millions of database and query sequences make BLAST computationally challenging but also well suited for parallelization on high-performance computing clusters. The efficacy of parallelization depends on the data partitioning, where the optimal data partitioning relies on an accurate performance model. In previous studies, a BLAST job was sped up by 27 times by partitioning the database and query among thousands of processor nodes. However, the optimality of the partitioning method was not studied. Unlike BLAST performance models proposed in the literature that usually have problem size and hardware configuration as the only variables, the execution time of a BLAST job is a function of database size, query size, and hardware capability. In this work, the nucleotide BLAST application BLASTN was profiled using three methods: shell-level profiling with the Unix “time” command, code-level profiling with the built-in “profiler” module, and system-level profiling with the Unix “gprof” program. The runtimes were measured for six node types, using six different database files and 15 query files, on a heterogeneous HPC cluster with 500+ nodes. The empirical measurement data were fitted with quadratic functions to develop performance models that were used to guide the data parallelization for BLASTN jobs. Results Profiling results showed that BLASTN contains more than 34,500 different functions, but a single function, RunMTBySplitDB, takes 99.12% of the total runtime. Among its 53 child functions, five core functions were identified to make up 92.12% of the overall BLASTN runtime. Based on the performance models, static load balancing algorithms can be applied to the BLASTN input data to minimize the runtime of the longest job on an HPC cluster. Four test cases being run on homogeneous and heterogeneous clusters were tested. Experiment results showed that the runtime can be reduced by 81% on a homogeneous cluster and by 20% on a heterogeneous cluster by re-distributing the workload. Discussion Optimal data partitioning can improve BLASTN’s overall runtime 5.4-fold in comparison with dividing the database and query into the same number of fragments. The proposed methodology can be used in the other applications in the BLAST+ suite or any other application as long as source code is available.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Peatland hydrology across scale: a probabilistic framework for confronting variability, heterogeneity, and uncertainty (Final Report)

Our project aimed to guide the representations of peatland hydrology and water-carbon feedbacks within Earth System Models, by improving process understanding and developing parsimonious models that elucidate the effects of hydroclimatic variability and spatial heterogeneity. Our research questions centered on the three key drivers of peatland hydrology: (1) seasonal and interannual hydroclimatic fluctuations, (2) spatial heterogeneity of peatland microtopography, and (3) hydrological connectivity across landscape units within a peatland watershed. Our analyses and model testing used existing long-term datasets within the Marcell Experimental Forest (MEF), with new data collected whenever necessary. To date, we have demonstrated the strong influence of water table elevations on the temperature sensitivity of CH 4 emissions, using a newly developed eddy covariance dataset spanning eleven years at MEF. Specifically, higher water tables dampen the increase in CH 4 emissions in the spring as well as their decrease in the fall, resulting in hysteresis. These results imply that any hydroclimatological changes in peatlands that shift seasonal water availability from winter to summer will increase annual CH 4 emissions, even if soil temperature remains unchanged. To further investigate water table elevation change within the peatland microtopography, we installed automated water table gauges across four bog-forest transects (the “lagg”) in two peatland watersheds at MEF, which are hotspots of intense biogeochemical activity. These measurements will give us new information about the extent of lagg expansion and contraction during snowmelt and high intensity rainfall events, as well as the directionality of water and nutrient flow into the bog from the surrounding forests. Finally, we investigated the hydrological connectivity across the peatland watershed complex during early spring – a critical period where water table is elevated from snowmelt – using extended hydrological records from MEF (e.g., snow and frost depths, water table elevations, streamflow). The importance of thermal and wintertime processes within the Energy Exascale Earth System Model (E3SM) was determined through sensitivity analyses. Results from data analyses show that (i) streamflow has decreased over decades likely due to increased evapotranspiration rates (despite no detectable trends in precipitation), (ii) spring streamflow generation is controlled by a sequence of fill-and-spill pathways from the snowpack to the stream through peatland water stores, and (iii) frost depth is a key explanatory variable for the timing and magnitude of streamflow. These results suggest that frost plays an important role in connecting surface water storage in peatlands to stream outlets, and that the hydrological connectivity across the peatland watershed complex can mediate the sensitivity of hydrological responses to climate variations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Photonuclear Physics in SCALE

Radiation shielding modeling involves a wide range of applications, such as power reactors, critical experiments, portable sources, linear accelerators, and fusion systems. Particle energy ranges vary from thermal energies to fast energy spectra. MAVRIC, as the radiation shielding sequence of the SCALE code system [1] developed and maintained by Oak Ridge National Laboratory, is capable of modeling both neutron and photon interactions in matter and is widely used by the US Nuclear Regulatory Commission, Department of Energy, the National Nuclear Security Administration, and radiation shielding practitioners all around the world. Although most applications don’t involve highly energetic photons that can cause photonuclear reactions and create coupled chains of high energy, ionized electrons and bremsstrahlung radiation within a material medium, when present, these reactions can have a significant impact upon the downstream particle fluxes from incident photons above 20 MeV. For example, linear accelerators used for medical isotope production are examples of applications in which high-energetic photons are generated and used for driving or controlling other components of a complex facility. Radiation shielding analyses of such a facility would require accounting for any and all significant sources of radiation. Unfortunately, MAVRIC included no data and methods to model either photonuclear reactions or the bremsstrahlung photons until this work, which enabled accounting for both. Photonuclear reaction cross sections were imported from the Evaluated Nuclear Data File libraries, and corresponding methods were implemented in the SCALE code system; see Section 2 for how the nuclear data are imported and used. Bremsstrahlung photons were modeled by assuming thick-target approximation whereby all generated electrons and positrons due to photon interactions are assumed to be absorbed within the same material in which they were born, yielding a quick particle transport simulation method by using fast lookup tables for generating subsequent gammas; see Section 4 for details. A suite of verification models was used to test implemented methods and associated data. Tests yielded satisfactory results for the new features. Enhancements will be released as part of the next major SCALE release (v7.0), which is anticipated to be deployed in 2025, with beta releases available in the interim.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Local-scale panmixia in the lichenized fungus Xanthoria parietina contrasts with substantial genetic structure in its Trebouxia photobionts

Microsatellite markers can provide valuable information about gene flow and population history. We developed and tested new microsatellites for the nitrophilic lichenized fungus Xanthoria parietina and studied its genetic diversity and structure within the urban area of Munich, Bavaria. We compared its local genetic pattern with that of its photobiont partner Trebouxia decolorans, for which existing microsatellites were applied. For comparison, a reference site with clean air was included in the sampling. We found support for three genetic clusters in the fungus X. parietina, which occurred intermingled in collecting sites. There was a high degree of admixture within fungal populations and individuals, and analysis of molecular variance revealed a lack of population structure in the mycobiont. The Trebouxia photobiont, in contrast, exhibited structured populations which grouped into two to five genetic clusters, and individuals showed less admixture than in the mycobiont. This indicates that the two lichen partners differ in their ability to move around in the landscape. The microsatellite markers we report are polymorphic and are suitable for population genetic studies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

P finder: genomic and metagenomic annotation of RNase P RNA gene (rnpB)

Abstract Background The rnpB gene encodes for an essential catalytic RNA (RNase P). Like other essential RNAs, RNase P’s sequence is highly variable. However, unlike other essential RNAs (i.e. tRNA, 16 S, 6 S,...) its structure is also variable with at least 5 distinct structure types observed in prokaryotes. This structural variability makes it labor intensive and challenging to create and maintain covariance models for the detection of RNase P RNA in genomic and metagenomic sequences. The lack of a facile and rapid annotation algorithm has led to the rnpB gene being the most grossly under annotated essential gene in completed prokaryotic genomes with only a 24% annotation rate. Here we describe the coupling of the largest RNase P RNA database with the local alignment scoring algorithm to create the most sensitive and rapid prokaryote rnpB gene identification and annotation algorithm to date. Results Of the 2772 completed microbial genomes downloaded from GenBank only 665 genomes had an annotated rnpB gene. We applied P Finder to these genomes and were able to identify 2733 or nearly 99% of the 2772 microbial genomes examined. From these results four new rnpB genes that encode the minimal T-type P RNase P RNAs were identified computationally for the first time. In addition, only the second C-type RNase P RNA was identified in Sphaerobacter thermophilus . Of special note, no RNase P RNAs were detected in several obligate endosymbionts of sap sucking insects suggesting a novel evolutionary adaptation. Conclusions The coupling of the largest RNase P RNA database and associated structure class identification with the P Finder algorithm is both sensitive and rapid, yielding high quality results to aid researchers annotating either genomic or metagenomic data. It is the only algorithm to date that can identify challenging RNAse P classes such as C-type and the minimal T-type RNase P RNAs. P Finder is written in C# and has a user-friendly GUI that can run on multiple 64-bit windows platforms (Windows Vista/7/8/10). P Finder is free available for download at https://github.com/JChristopherEllis/P-Finder as well as a small sample RNase P RNA file for testing.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Whole metagenome sequencing and 16S rRNA gene amplicon analyses reveal the complex microbiome responsible for the success of enhanced in-situ reductive dechlorination (ERD) of a tetrachloroethene-contaminated Superfund site

The North Railroad Avenue Plume (NRAP) Superfund site in New Mexico, USA exemplifies successful chlorinated solvent bioremediation. NRAP was the result of leakage from a dry-cleaning that operated for 37 years. The presence of tetrachloroethene biodegradation byproducts, organohalide respiring genera (OHRG), and reductive dehalogenase (rdh) genes detected in groundwater samples indicated that enhanced reductive dechlorination (ERD) was the remedy of choice. This was achieved through biostimulation by mixing emulsified vegetable oil into the contaminated aquifer. This report combines metagenomic techniques with site monitoring metadata to reveal new details of ERD. DNA extracts from groundwater samples collected prior to and at four, 23 and 39 months after remedy implementation were subjected to whole metagenome sequencing (WMS) and 16S rRNA gene amplicon (16S) analyses. The response of the indigenous NRAP microbiome to ERD protocols is consistent with results obtained from microcosms, dechlorinating consortia, and observations at other contaminated sites. WMS detects three times as many phyla and six times as many genera as 16S. Both techniques reveal abundance changes in Dehalococcoides and Dehalobacter that reflect organohalide form and availability. Methane was not detected before biostimulation but appeared afterwards, corresponding to an increase in methanogenic Archaea. Assembly of WMS reads produced scaffolds containing rdh genes from Dehalococcoides, Dehalobacter, Dehalogenimonas, Desulfocarbo, and Desulfobacula. Anaerobic and aerobic cometabolic organohalide degrading microbes that increase in abundance include methanogenic Archaea, methanotrophs, Dechloromonas, and Xanthobacter, some of which contain hydrolytic dehalogenase genes. Aerobic cometabolism may be supported by oxygen gradients existing in aquifer microenvironments or by microbes that produce O 2 via microbial dismutation. The NRAP model for successful ERD is consistent with the established pathway and identifies new taxa and processes that support this syntrophic process. This project explores the potential of metagenomic tools (MGT) as the next advancement in bioremediation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Microseismic at HFTS2: A Story of Three Stimulated Wells

A microseismic data set was acquired at Hydraulic Fracturing Test Site 2 (HFTS2) in the Delaware Basin (Texas, USA) by the HFTS2 consortium as a part of data-collection program designed to help maximize the efficiency of hydraulic fracturing, understand the impacts of nearby depleted zones, and increase the production of hydrocarbons. Advanced methods of microseismic monitoring applied at HFTS2 include 1) native synchronization of five monitor arrays, 2) real-time location of microseismic events, 3) joint event hypocenter/ velocity-model inversion that locates microseismic events simultaneously with constructing azimuthally anisotropic velocity models, 4) high-resolution reservoir imaging, utilizing microseismic events as downhole sources. HFTS2 microseismic data set illustrated the influence of geologic complexity, variations in treatment designs, and the presence of parent wells on the efficiency of hydraulic fracturing. Three Wolfcamp wells, Boxwood 1H, 2H, and 4H, were stimulated sequentially, exhibiting very different patterns of microseismicity. When the Boxwood 1H well is stimulated, its toe stages quickly reacted to the presence of depleted Bitterroot section to the east of the Boxwood pad, extending microseismic events over 3,000 ft to the east. This behavior changes dramatically after the treatment passes the Bitterroot section at the toe, and the heel stages in 1H showed tight microseismic grouping around the stimulated well, indicative of treatment of virgin reservoir rock. The Boxwood 2H demonstrates this same pattern of behavior between the toe and heel stages of the well, suggesting that the depleted Bitterroot section remains easier to refracture than adjacent non-depleted rock. Microseismic activity associated with Boxwood 4H treatment is distinctly different from those in Boxwood wells 1H and 2H. The events generated during stimulation of 4H are hardly influenced by the depleted Bitterroot section and fill the space between the previously stimulated 1H and 2H wells. Altogether, microseismic patterns suggest that treatment fluids tend to find the weakest passages through geologic formations, and the relative well positions and treatment sequence have a significant impact on fracture geometry, as indicated by microseismic event locations.

58 GEOSCIENCES↗
DOE OSTI2024

Comparing Transcriptomic Points of Departure to Apical Effect Concentrations For Larval Fathead Minnow Exposed to Chemicals with Four Different Modes Of Action

It is postulated that below a transcriptomic-based point of departure, adverse effects are unlikely to occur, thereby providing a chemical concentration to use in screening level hazard assessment. The present study extends previous work describing a high-throughput fathead minnow assay that can provide full transcriptomic data after exposure to a test chemical. One-day post-hatch fathead minnows were exposed to ten concentrations of three representatives of four chemical modes of action: organophosphates, ecdysone receptor agonists, plant photosystem II inhibitors, and estrogen receptor agonists for 24 h. Concentration response modeling was performed on whole body gene expression data from each exposure, using measured chemical concentrations when available. Transcriptomic points of departure in larval fathead minnow were lower than apical effect concentrations across fish species but not always lower than toxic effect concentrations in other aquatic taxa like crustaceans and insects. The point of departure was highly dependent on measured chemical concentration which were often lower than the nominal concentration. Differentially expressed genes between chemicals within modes of action were compared and often showed statistically significant overlap. In addition, reproducibility between identical exposures using a positive control chemical (CuSO 4 ) and variability associated with the transcriptomic point of departure using in silico sampling were considered. Results extend a transcriptomic-compatible fathead minnow high-throughput assay for possible use in ecological hazard screening.

59 BASIC BIOLOGICAL SCIENCES↗