Search NASA⌕ Search

SEARCH · Search NASA

Results for “sequence analysis”

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 451 records · Page 25

Understanding the transient large amplitude oscillatory shear behavior of yield stress fluids

A full understanding of the sequence of processes exhibited by yield stress fluids under large amplitude oscillatory shearing is developed using multiple experimental and analytical approaches. A novel component rate Lissajous curve, where the rates at which strain is acquired unrecoverably and recoverably are plotted against each other, is introduced and its utility is demonstrated by application to the analytical responses of four simple viscoelastic models. Using the component rate space, yielding and unyielding are identified by changes in the way strain is acquired, from recoverably to unrecoverably and back again. The behaviors are investigated by comparing the experimental results with predictions from the elastic Bingham model that is constructed using the Oldroyd–Prager formalism and the recently proposed continuous model by Kamani, Donley, and Rogers in which yielding is enhanced by rapid acquisition of elastic strain. The physical interpretation gained from the transient large amplitude oscillatory shear (LAOS) data is compared to the results from the analytical sequence of physical processes framework and a novel time-resolved Pipkin space. The component rate figures, therefore, provide an independent test of the interpretations of the sequence of physical processes analysis that can also be applied to other LAOS analysis frameworks. Each of these methods, the component rates, the sequence of physical processes analysis, and the time-resolved Pipkin diagrams, unambigiously identifies the same material physics, showing that yield stress fluids go through a sequence of physical processes that includes elastic deformation, gradual yielding, plastic flow, and gradual unyielding.

Kamani, Krutarth M. (ORCID:0000000338975420)↗

Kinetic characterization of Vibrio cholerae ApbE: Substrate specificity and regulatory mechanisms

ApbE is a member of a novel family of flavin transferases that incorporates flavin mononucleotide (FMN) to subunits of diverse respiratory complexes, which fulfill important homeostatic functions. In this work a detailed characterization of Vibrio cholerae ApbE physiologic activity, substrate specificity and pH dependency was carried out. The data obtained show novel characteristics of the regulation and function of this family. For instance, our experiments indicate that divalent cations are essential for ApbE function, and that the selectivity depends largely on size and the coordination sphere of the cation. Our data also show that ApbE regulation by pH, ADP and potassium is an important mechanism that enhances the adaptation, survival and colonization of V. cholerae in the small intestine. Moreover, studies of the pH-dependency of the activity show that the reaction is favored under alkaline conditions, with a pKa of 8.4. These studies, together with sequence and structure analysis allowed us to identify His257, which is absolutely conserved in the family, as a candidate for the residue whose deprotonation controls the activity. Remarkably, the mutant H257G abolished the flavin transfer activity, strongly indicating that this residue plays an important role in the catalytic mechanism of ApbE.

59 BASIC BIOLOGICAL SCIENCES↗

Data Analysis of the 2020 Central Idaho Mainshock-Aftershock Sequence

In an effort to inform the Senior Seismic Hazard Analysis Committee for the Idaho National Laboratory, we provide an improved aftershock catalog related to the March 31, 2020, Mw6.5 Stanley, Idaho earthquake from picks related to a temporary network of two real-time and 15 non-telemetered seismometers within the epicentral area. From the permanent and temporary (XP) real-time network, the USGS cataloged 1,946 aftershocks between April 1, 2020 and October 31, 2020. To improve aftershock location and magnitudes, we manually picked arrival times of P and S waves from off-line stations in the XP temporary network, generated a new crustal velocity model, and independently relocated each event using the HypoDD double-difference earthquake algorithm. We created our new velocity model from existing broadband and active source seismic campaign data that were acquired near the epicentral region prior to the 2020 earthquake. We compare arrival time differences, epicentral locations and depths between aftershocks recorded with the two catalogs. We find the addition of local stations provides tighter aftershock clustering that suggests an improved aftershock locations. To detect lower magnitude events, we employed deep learning. Our method solves common problems associated with detecting many events that have a low signal-to-noise ratio. From the machine learning database, we detected more than 74,000 aftershocks. Based on the number of identified earthquakes and Gutenberg-Richter relationships derived from the USGS catalog, we estimate that we have reduced the completion magnitude for the Stanley earthquake sequence to below M1 using this machine learning approach. We located each aftershock with our new velocity model. Our new velocity model and picks suggests aftershocks occurred mostly at shallower depths than assessed in the USGS catalog. These aftershocks align along two linear trends that suggest the activation of two unnamed primary faults.

58 GEOSCIENCES↗

Reinspection of a Clinical Proteomics Tumor Analysis Consortium (CPTAC) Dataset with Cloud Computing Reveals Abundant Post-Translational Modifications and Protein Sequence Variants

The Clinical Proteomic Tumor Analysis Consortium (CPTAC) has provided some of the most in-depth analyses of the phenotypes of human tumors ever constructed. Today, the majority of proteomic data analysis is still performed using software housed on desktop computers which limits the number of sequence variants and post-translational modifications that can be considered. The original CPTAC studies limited the search for PTMs to only samples that were chemically enriched for those modified peptides. Similarly, the only sequence variants considered were those with strong evidence at the exon or transcript level. In this multi-institutional collaborative reanalysis, we utilized unbiased protein databases containing millions of human sequence variants in conjunction with hundreds of common post-translational modifications. Using these tools, we identified tens of thousands of high-confidence PTMs and sequence variants. We identified 4132 phosphorylated peptides in nonenriched samples, 93% of which were confirmed in the samples which were chemically enriched for phosphopeptides. In addition, our results also cover 90% of the high-confidence variants reported by the original proteogenomics study, without the need for sample specific next-generation sequencing. Finally, we report fivefold more somatic and germline variants that have an independent evidence at the peptide level, including mutations in ERRB2 and BCAS1. In this reanalysis of CPTAC proteomic data with cloud computing, we present an openly available and searchable web resource of the highest-coverage proteomic profiling of human tumors described to date.

60 APPLIED LIFE SCIENCES↗

Small-Magnitude Seismic Swarms in Central Utah (US): Interactions of Regional Tectonics, Local Structures and Hydrothermal Systems

Swarms in Central Utah are situated in the complex transition between the Basin and Range (BR) province and the Colorado Plateau. Transecting transverse structures, volcanic deposits, and hydrothermal systems complicate the extensional BR horst and graben structures and provide a multitude of plausible triggering mechanisms. Revisiting the catalog of the University of Utah Seismograph Stations (1981–2022), we analyze spatio-temporal patterns and characteristic features of seismic sequences. Swarms with alternating seismicity rates, bursts, and longer swarms with persistent moment release exhibit a remarkable diversity in temporal evolution. Swarm durations do not scale with cumulative seismic moment: swarms lasting less than 1 day can have similar cumulative seismic moments as month-long swarms. Here, we observe stationary swarms re-occurring for years (e.g., Mineral Mountains), as well as singular swarms in low-seismicity areas (e.g., activating a local structure). The swarms show a pronounced heterogeneity in triggering and driving mechanisms, observed in the detailed analysis of exemplary sequences (detections, relocations, moment tensors, waveform-based clustering, and repeater analysis). The 2022 Sevier Valley sequence activated a BR-related normal fault, the first resolved fault plane in the valley since 1983. The 2011 Circleville sequence is interpreted as a swarm triggered by mainshock-aftershock activity characterized by increasing magnitudes, changing rupture mechanisms, and a concentration of highly similar events in the second part of the sequence. By jointly discussing exemplary sequences and catalog statistics, we draw a comprehensive picture of swarm activity and its relation to geothermal and tectonic activity.

58 GEOSCIENCES↗

Diverse Thaumarchaeota Dominate Subsurface Ammonia-oxidizing Communities in Semi-arid Floodplains in the Western United States

Subsurface microbial communities mediate biogeochemical transformations that drive both local and ecosystem-level cycling of essential elements, including nitrogen. However, their study has been largely limited to the deep ocean, terrestrial mines, caves, and topsoils (< 30 cm). Here, we present regional insights into the microbial ecology of aerobic ammonia oxidation within the terrestrial subsurface of five semi-arid riparian sites spanning a 900-km N-S transect. We sampled sediments, profiled communities to depths of ≤10 m, and compared them to reveal trends regionally within and surrounding the Upper Colorado River Basin (CRB). The diversity and abundance of ammonia-oxidizing microbial communities were evaluated in the context of subsurface geochemistry by applying a combination of amoA (encoding ammonia monooxygenase subunit A) gene sequencing, quantitative PCR, and geochemical techniques. Analysis of 898 amoA sequences from ammonia-oxidizing archaea (AOA) and bacteria (AOB) revealed extensive ecosystem-scale diversity, including archaeal amoA sequences from four of the five major AOA lineages currently found worldwide as well as distinct AOA ecotypes associated with naturally reduced zones (NRZs) and hydrogeochemical zones (unsaturated, capillary fringe, and saturated). Overall, AOA outnumber AOB by 2- to 5000-fold over this regional scale, suggesting that AOA may play a prominent biogeochemical role in nitrification within terrestrial subsurface sediments.

54 ENVIRONMENTAL SCIENCES↗

Data for Intra- and inter-annual variability of nitrification in the rhizosphere of field-grown bioenergy sorghum

These data were collected in 2018 and 2019 at the University of Illinois Energy Farm (N 40.063607, W 88.206926). During each growing season, bulk and rhizosphere soil were collected from replicate Sorghum bicolor nitrogen use efficiency trial plots at three separate time points (approximately July 1, August 1, and September 1). We measured soil moisture, pH, soil nitrate and ammonium, potential nitrification, potential denitrification, and extracted and sequenced the V4 region of the 16S rRNA gene for microbial community analysis. All microbial sequence data is archived in the National Center for Biotechnology Information’s (NCBI) Sequence Read Archive (accession number SRP326979, project number PRJNA741261).

bioenergy↗

Calibration and Sequence Development Status for the Sample Analysis at Mars Investigation on the Mars Science Laboratory

The measurement goals of the Sample Analysis at Mars (SAM) instrument suite on the "Curiosity" Rover of the Mars Science Laboratory (MSL) include chemical and isotopic analysis of organic and inorganic volatiles for both atmospheric and solid samples [1,2]. SAM directly supports the ambitious goals of the MSL mission to provide a quantitative assessment of habitability and preservation in Gale crater by means of a range of chemical and geological measurements [3]. The SAM FM combined calibration and environmental testing took place primarily in 2010 with a limited set of tests implemented after integration into the rover in January 2011. The scope of SAM FM testing was limited both to preserve SAM consumables such as life time of its electromechanical elements and to minimize the level of terrestrial contamination in the SAM instrument. A more comprehensive calibration of a SAM-like suite of instruments will be implemented in 2012 with calibration runs planned for the SAM testbed. The SAM Testbed is nearly identical to the SAM FM and operates in a ambient pressure chamber. The SAM Instrument Suite: SAM's instruments are a Quadrupole Mass Spectrometer (QMS), a 6-column Gas Chromatograph (GC), and a 2-channel Tunable Laser Spectrometer (TLS). Gas Chromatography Mass Spectrometry is designed for identification of even trace organic compounds. The TLS [5] secures the C, H, and O isotopic composition in carbon dioxide, water, and methane. Sieved materials are delivered from the MSL sample acquisition and processing system to one of68 cups of the Sample Manipulation System (SMS). 59 of these cups are fabricated from inert quartz. After sample delivery, a cup is inserted into one of 2 ovens for evolved gas analysis (EGA ambient to >9500C) by the QMS and TLS. A portion of the gas released can be trapped and subsequently analyzed by GCMS. Nine sealed cups contain liquid solvents and chemical derivatization or thermochemolysis agents to extract and transform polar molecules such as amino acids, nucleobases, and carboxylic acids into compounds that are sufficiently volatile to transmit through the GC columns. The remaining 6 cups contain calibrants. SAM FM Calibration Overview: The SAM FM calibration in the Mars chamber employed a variety of pure gases, gas mixtures, and solid materials. Isotope calibration runs for the TLS utilized 13C enriched C02 standards and 0 enriched CH4. A variety of fluorocarbon compounds that spanned the entire mass range of the QMS as well as C3-C6 hydrocarbons were utilized for calibration of the GCMS. Solid samples consisting of a mixture of calcite, melanterite, and inert silica glass either doped or not with fluorocarbons were introduced into the SAM FM cups through the SAM inlet funnel/tube system.

Mahaffy, Paul R.↗

Statistical and linguistic features of DNA sequences

We present evidence supporting the idea that the DNA sequence in genes containing noncoding regions is correlated, and that the correlation is remarkably long range--indeed, base pairs thousands of base pairs distant are correlated. We do not find such a long-range correlation in the coding regions of the gene. We resolve the problem of the "non-stationary" feature of the sequence of base pairs by applying a new algorithm called Detrended Fluctuation Analysis (DFA). We address the claim of Voss that there is no difference in the statistical properties of coding and noncoding regions of DNA by systematically applying the DFA algorithm, as well as standard FFT analysis, to all eukaryotic DNA sequences (33 301 coding and 29 453 noncoding) in the entire GenBank database. We describe a simple model to account for the presence of long-range power-law correlations which is based upon a generalization of the classic Levy walk. Finally, we describe briefly some recent work showing that the noncoding sequences have certain statistical features in common with natural languages. Specifically, we adapt to DNA the Zipf approach to analyzing linguistic texts, and the Shannon approach to quantifying the "redundancy" of a linguistic text in terms of a measurable entropy function. We suggest that noncoding regions in plants and invertebrates may display a smaller entropy and larger redundancy than coding regions, further supporting the possibility that noncoding regions of DNA may carry biological information.

Non-NASA Center↗

Sequence Length of HIV-1 Subtype B Increases over Time: Analysis of a Cohort of Patients with Hemophilia over 30 Years

We aimed to investigate whether the sequence length of HIV-1 increases over time. We performed a longitudinal analysis of full-length coding region sequences (FLs) during an HIV-1 outbreak among patients with hemophilia and local controls infected with the Korean subclade B of HIV-1 (KSB). Genes were amplified by overlapping RT-PCR or nested PCR and subjected to direct sequencing. Overall, 141 FLs were sequentially determined over 30 years in 62 KSB-infected patients. Phylogenetic analysis indicated that within KSB, two FLs from plasma donors O and P comprised two clusters, together with 8 and 12 patients with hemophilia, respectively. Signature pattern analysis of the KSB of HIV-1 revealed 91 signature nucleotide residues (1.1%). In total, 48 and 43 signature nucleotides originated from clusters O and P, respectively. Six positions contained 100% specific nucleotide(s) in clusters O and P. In-depth FL analysis for over 30 years indicated that the KSB FL significantly increased over time before combination antiretroviral therapy (cART) and decreased with cART. This increase occurred due to the significant increase in env and nef genes, originating in the variable regions of both genes. The increase in sequence length of HIV-1 over time suggests an evolutionary direction.

59 BASIC BIOLOGICAL SCIENCES↗

MemGaze: Rapid and Effective Load-Level Memory Trace Analysis

A major challenge of memory analysis tools is combining high-resolution analysis and low overhead measurement. Currently, hardware/software-based analysis of load-level sequences incurs time slowdowns of O(100×). We present MemGaze, a tool for low-overhead, high-resolution memory analysis. MemGaze uses Intel’s Processor Tracing (PT) instruction ptwrite to collect sampled and compressed memory address traces for load-level, sequence-aware analysis of data reuse. We describe multi-resolution analysis for locations vs. operations, accesses vs. spatio-temporal reuse, and reuse (distance, rate, volume) vs. access patterns. Both trace size and resolution are controllable. We use MemGaze to elucidate the memory effects of different data structures and algorithms. For sampled traces that are ˜1% of a full one, analysis metrics have 1-25% MAPE for histograms of varying dynamic sequence lengths. With current suboptimal kernel support (PT runs continuously), MemGaze’s time overhead is typically 10–95%; 7× at worst. However, when PT runs only during samples, overhead is 10–35% on memory intensive regions and correlates with executed ptwrites.

Kilic, Ozgur O.↗

RISK-INFORMED EXTERNAL HAZARDS ANALYSIS BY EEVE TOOLKIT

Introducing risk-informed safety analyses for external events can bring sensible benefits to the nuclear power industry easing the regulatory burden introduced in the aftermath of the Fukushima accident. In the framework of the US-DOE Light-Water Reactor Sustainability program, Risk-Informed Margin Characterization project (LWRS-RISMC), research activities were conducted at INL for developing and applying a toolkit for performing risk-informed external hazards analysis. The “EEVE” (External EVEnts) toolkit, a chain of INL and state-of-the-art codes, was set-up and applied for simulating earthquake-induced transients with internal flooding in a PWR. This was done using methodology set-up for evaluating new data and determining critical areas that would benefit most from advanced analysis methods. The methodology is agnostic and it can be applied to other external hazards than just seismic and flooding (e.g., high winds, intense precipitation, etc.). The first step was the deterministic calculations for evaluating the propagation of a set of earthquakes and for determining the effects on a fire suppression system in an auxiliary building. These deterministic calculations were performed using the LS-DYNA and the OPENSEE structural mechanics codes. Probability failures and event trees of a generic 3-loops Westinghouse PWR were derived from a generic SAPHIRE model. Seismically induced Loss of Offsite Power, Station Blackout, and fire suppression piping failure events were considered. SAPHIRE code allowed a pre-screening of all possible branches of the event tree and a selection of those with the largest Core Damage Frequency increase for further analysis using dynamic PRA methods. INL's EMRALD code performed dynamic PRA calculations of the identified sequences, integrating risk analysis with on-line deterministic safety analysis results (seismically-induced internal flooding calculations by NEUTRINO code and system analysis by RELAP5-3D code). RELAP5-3D Best-Estimate calculations coupled with RAVEN statistical code were used for deriving limit surfaces as function of relevant transient parameters. This data was used by EMRALD for determining the final core status (core damage/core safe). Finally, the coupling of seismic, flooding, system thermal-hydraulic, uncertainty quantification and probabilistic risk analysis for an advanced risk-informed external hazards analysis was demonstrated. Results showed the sensible decrease in the level of conservatism achieved using the developed risk-informed methodology, which allowed also improving the evaluation of the failure modes and failure probabilities of different components.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Benchmarking Computational Tools for Calling SNPs and Indels in Complex Microbial Populations

The NASA BioNutrients missions seek to understand the suitability of microorganisms for bioproduction during space flight. One topic of interest is the stability of microbial genomes during long-term ambient storage and subsequent rehydration and growth. To address these questions, samples from 8 species were flown to ISS for 5 years of desiccated storage at ambient temperature (Stasis Packs) and 2 species were packaged along with powdered media inside a bioreactor system to allow hydration and growth in microgravity (Production Packs). For both systems, Whole Genome Sequencing (WGS) of the DNA extracted from the returned samples and paired ground controls will be conducted to identify changes in genome stability due to time, storage conditions and growth in space. Across the technical replicates, ground controls, 10 timepoints, and multiple experimental conditions, ~300 samples have been selected for initial analysis with WGS sequencing to 100x coverage. A flexible and resource efficient mutation calling pipeline is needed to process this large dataset and allow for comparisons between species. Many bioinformatics tools for calling Indels and Single Nucleotide Variants (SNVs) are designed for use with pure isolates, where true variations from the reference genome are expected to dominate the reads aligning to the location of mutation. In contrast, DNA from the Stasis Pack (SP) samples was collected directly after recovery from desiccated storage and the Production Pack (PP) samples were collected after fermentation. In this context, reads with mutations are expected to be less frequent than reads that align with the reference genome, as each sample will include multiple lines of cells. Thus, BioNutrients samples are expected to be similar to samples from cancer cell or “pooled” sequencing approaches. In preparation for the analysis of the BioNutrients samples, we have tested three mutation calling tools (GATK for Microbes, BreSeq and DiscoSNP) designed for complex samples. A challenge of validating mutation identification pipelines is a lack of “Ground Truth” datasets, especially for complex samples. To compare these three tools, we sought to identify mutations in pre-existing WGS data collected from populations of Chlamydomonas reinhardtii that were exposed to UV mutagenesis and growth in LEO as part of the Space Algae-1 mission. Here we present a summary of these tools against the analysis originally conducted using the CRISP tool. Critical metrics are compared such as runtime, the number of SNPs, the number and size of Indels, and patterns of transversion and transitions identified by each tool are reported. By sharing these benchmarking results collected in support of the BioNutrients mission, we aim to guide others seeking to identify SNVs in similarly complex microbial samples.

Biology↗

The metagenomic data life-cycle: standards and best practices

Metagenomics data analyses from independent studies can only be compared if the analysis workflows are described in a harmonised way. In this overview, we have mapped the landscape of data standards available for the description of essential steps in metagenomics: (1) material sampling, (2) material sequencing (3) data analysis and (4) data archiving & publishing. Taking examples from marine research, we summarise essential variables used to describe material sampling processes and sequencing procedures in a metagenomics experiment. These aspects of metagenomics dataset generation have been to some extent addressed by the scientific community but greater awareness and adoption is still needed. We emphasise the lack of standards relating to reporting how metagenomics datasets are analysed and how the metagenomics data analysis outputs should be archived and published. We propose best practice as a foundation for a community standard to enable reproducibility and better sharing of metagenomics datasets, leading ultimately to greater metagenomics data reuse and repurposing.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗