Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sequence Alignment”

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 397 records · Page 22

UAE6 - Wind Tunnel Tests Data - UAE6 - Sequence 9 - Raw Data

Sequences 8 and 9: Downwind Sonics (F,P) and Downwind Sonics Parked (P) This test sequence used an upwind, rigid turbine with a 0° cone angle. The wind speed ranged from 5 m/s to 25 m/s. Yaw angles of 0° to 60° were achieved. The blade tip pitch was 3°. The rotor rotated at 72 RPM during Sequence 8, but it was parked during Sequence 9. Blade pressure measurements were collected. The five-hole probes were removed and the plugs were installed. Plastic tape 0.03-mm-thick was used to smooth the interface between the plugs and the blade. The teeter dampers were replaced with rigid links, and these two channels were flagged as not applicable by setting the measured values in the data file to –99999.99 Nm. The teeter link load cell was pre-tensioned to 40,000 N. During post-processing, the probe channels were set to read -99999.99. Sonic anemometers were mounted on a strut downwind of the turbine. The strut was mounted to the T-frame, which was rotated to align the anemometers aft of the 9% and 49% radius locations at hub height. Because of this configuration, the tunnel balance data are considered invalid. Sequence 9 was designed to compare the downwind sonic anemometer readings with the upwind sonic anemometers without interference from the turbine. The rotor was parked with the instrumented blade at 0° azimuth. All pressure measurements obtained in Sequence 9 are invalid because sufficient time for temperature stabilization did not occur, thus all associated data values were flagged as not applicable by setting the measured values in the data file to 0.0000 Pa. This test is further described in Appendix G.

17 WIND ENERGY↗

The Youngest Planet to Have a Spin-Orbit Alignment Measurement AU Mic b

We report measurements of the sky-projected spin–orbit angle for AU Mic b, a Neptune-size planet orbiting a very young (∼20 Myr) nearby pre-main-sequence M-dwarf star, which also hosts a bright, edge-on, debris disk. The planet was recently discovered from preliminary analysis of radial-velocity observations and confirmed to be transiting its host star from photometric data from the NASA’s TESS mission. We obtained radial-velocity measurements of AU Mic over the course of two partially observable transits and one full transit of planet b from high-resolution spectroscopic observations made with the Minerva-Australis telescope array. Only a marginal detection of the Rossiter–McLaughlin effect signal was obtained from the radial velocities, in part due to AU Mic being an extremely active star and the lack of full transit coverage plus sufficient out-of-transit baseline. As such, a precise determination of the obliquity for AU Mic b is not possible in this study and we find a sky-projected spin–orbit angle of λ=47{sub −54}{sup +26°}. This result is consistent with both the planet’s orbit being aligned or highly misaligned with the spin axis of its host star. Our measurement independently agrees with, but is far less precise than observations carried out on other instruments around the same time that measure a low-obliquity orbit for the planet. AU Mic is the youngest exoplanetary system for which the projected spin–orbit angle has been measured, making it a key data point in the study of the formation and migration of exoplanets—particularly given that the system is also host to a bright debris disk.

47 OTHER INSTRUMENTATION↗

CBCT with specification of imaging dose and CNR by anatomical volume of interest

Purpose: A novel method has been developed for volume of interest (VOI) cone-beam CT (CBCT) imaging using a 2.35 MV/Carbon target linac imaging beam line combined with dynamic multileaf collimator sequences. Methods: The authors demonstrate the concept of acquisition of multiple, separate imaging volumes, where volumes can be either completely separated or nested, and are associated with predetermined imaging dose and contrast-to-noise ratio (CNR) characteristics. Two individual MLC sequences were established in the planning system (Eclipse, Varian Medical) to collimate the beam according to a defined inner VOI (e.g., containing a target volume under image guidance) and an outer VOI (e.g., including surrounding landmarks or organs-at-risk). MLC sequences were interleaved as a function of gantry angle to produce a reconstructed CBCT image with nested VOIs. By controlling the ratio of inner-to-outer ratio of MLC segments (and thus Monitor Units) during acquisition, the relative dose and CNR in the two volumes can be controlled. Inner-to-outer ratios of 2:1 to 6:1 were examined. Results: The concept was explored using an anatomical head phantom to assess image quality. A geometric phantom was used to quantify absolute dose and CNR values for the various sequences. The authors found that the dose in the outer VOI decreased by a functional relationship dependent on the inner-to-outer sequence ratio, while the CNR varied by the square root of dose, as expected. Conclusions: In this study the authors demonstrate flexibility in VOI CBCT by tailoring the imaging dose and CNR distribution in separate volumes within the patient anatomy. This would allow for high quality imaging of a target volume for alignment purposes, with simultaneous low dose imaging of the surrounding anatomy (e.g., for coregistration)

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Genome analysis and physiological comparison of Alicycliphilus denitrificans strains BC and K601T

The genomes of the Betaproteobacteria Alicycliphilus denitrificans strains BC and K601T have been sequenced to get insight into the physiology of the two strains. Strain BC degrades benzene with chlorate as electron acceptor. The cyclohexanol-degrading denitrifying strain K601T is not able to use chlorate as electron acceptor, while strain BC cannot degrade cyclohexanol. The 16S rRNA sequences of strains BC and K601T are identical and the fatty acid methyl ester patterns of the strains are similar. Basic Local Alignment Search Tool (BLAST) analysis of predicted open reading frames of both strains showed most hits with Acidovorax sp. JS42, a bacterium that degrades nitro-aromatics. The genomes include strain-specific plasmids (pAlide201 in strain K601T and pAlide01 and pAlide02 in strain BC). Key genes of chlorate reduction in strain BC were located on a 120 kb megaplasmid (pAlide01), which was absent in strain K601T. Genes involved in cyclohexanol degradation were only found in strain K601T. Benzene and toluene are degraded via oxygenase-mediated pathways in both strains. Genes involved in the meta-cleavage pathway of catechol are present in the genomes of both strains. Strain BC also contains all genes of the ortho-cleavage pathway. The large number of mono- and dioxygenase genes in the genomes suggests that the two strains have a broader substrate range than known thus far.

Oosterkamp, Margreet J.↗

Bayesian computational approaches for gene regulation studies of bioethanol and biohydrogen production

It has recently become clear that regulatory RNAs play a major role in regulation of gene expression in bacteria. RNA secondary structures play a major role in the function of many regulatory RNAs, and structural features are often key to their interaction with other cellular components. Thus, there has been considerable interest in the prediction of the secondary structures for RNA families. A paper describing our new algorithm, RNAG, to predict consensus secondary structures for unaligned sequences using the blocked Gibbs sampler has been published[1]. This sampling algorithm iteratively samples from the conditional probability distributions: P(Structure | Alignment) and P(Alignment | Structure). Subsequent to publication of the RNAG paper we have employed the technology from RNAG in the development of an RNA motif finding algorithm. To develop and RNA motif finding algorithm, RGibbs, we capitalized on our long experience in DNA motif finding and RNA secondary structure prediction. We applied RGibbs to three data sets from the literature and compared it to existing methods: one for training and two others for tests sets. In both test sets we found RGibbs out performed existing procedures.

59 BASIC BIOLOGICAL SCIENCES↗

Comparative genomics analysis of Trichoderma reesei strains

Trichoderma reesei is a key fungus for industrial production of lignocellulolytic enzymes. The genome sequences of the Trichoderma reesei hyper-cellulolytic strain RUT-C30 and its parental strain QM6a were compared at the nucleotide level. Approximately 97 % of the 87 genomic sequence scaffolds of T. reesei QM6a (33Mb) were detected to have the correspondent nucleotide in the 182 genome sequence scaffolds of RUT-C30 (32Mb). There are 455 loci within the QM6 sequence not detected in the RUT-C30 sequence. Regions at the termini of QM6a scaffolds as well as 14 small scaffolds do not have corresponding regions in RUT-C30 genomic scaffolds. Seventy-eight protein encoding genes are included within these regions. Mutated nucleotide(s) in 2,371 positions, including short insertion/deletions (indels) were detected in the aligned regions. The predicted protein coding regions of 97 gene models contain mutations, 34 of which were not previously described. Twenty-seven out of 34 newly discovered genes were found to have mutations in the peptide amino acid sequence. This is in addition to 63 genes described in a previous study based on low coverage sequencing of RUT-C30. These newly identified proteins are involved in signal transduction, transcription, RNA processing and modification, post-translational modification according to their annotations. Similar distributions of eukaryotic orthologous group (KOG) categories between the mutated and all other proteins suggest random mutation. The roles of the mutated genes and potential regulatory regions in the observed phenotype of RUT-C30 remain to be explored in a targeted fashion.

Trichoderma, genomics, biofuels, enzymes, cellulas↗

First observation of rotational structures in Re 168

We assigned first rotational sequences to the odd-odd nucleus 168 Re. Coincidence relationships of these structures with rhenium x rays confirm the isotopic assignment, while arguments based on the γ-ray multiplicity (K-fold) distributions observed with the new bands lead to the mass assignment. Configurations for the two bands were determined through analysis of the rotational alignments of the structures and a comparison of the experimental B(M1)/B(E2) ratios with theory. Tentative spin assignments are proposed for the πh 11/2 νi 13/2 band, based on energy level systematics for other known sequences in neighboring odd-odd rhenium nuclei, as well as on systematics seen for the signature inversion feature that is well known in this region. Furthermore, the spin assignment for the πh 11/2 ν(h 9/2 /f 7/2 ) structure provides additional validation of the proposed spins and configurations for isomers in the 176 Au → 172 Ir → 168 Re α-decay chain.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Accretion disks in the IRAS 23151+5912 region

We present observations of radio continuum emission at 1.3 and 3.6 cm and H{sub 2}O masers toward the high-mass star-forming region IRAS 23151+5912 carried out with the Very Large Array-Expanded Very Large Array (in transition phase) in configuration A. We detected one continuum source at 1.3 cm in the region, but the counterpart of this radio continuum source at 3.6 cm was not detected at a level of 3σ. We also detected 13 water maser spots toward IRAS 23151+5912, which are distributed in three groups aligned along the northeast-southwest direction. Our results suggest that the 1.3 cm emission is consistent with a hypercompact H II region, probably with an embedded zero-age main-sequence star of type B2. In particular, we find that this young star is spatially associated with a maser group, which is tracing a disk-like structure of about 460 AU. We also find that the masers of the second group are probably describing a circumstellar disk of about 86 AU, whose central protostar, still undetected, should have a mass of ∼11 M {sub ☉}. We also suggest that the third water maser group is possibly associated with the SiO outflow and the undetected driving source of the system. Finally, we noted that the 1.3 cm continuum source and the three maser groups are aligned in the northeast-southwest direction, similar to the elongation of the large ionized region, which probably is the result of shock-wave induced star formation on the densest region of the medium.

79 ASTRONOMY AND ASTROPHYSICS↗

Aftershocks of the 2012 Off-Coast of Sumatra Earthquake Sequence

Aftershocks of the 2012 Off-Coast of Sumatra Earthquake Sequence exhibit a complex and diffuse spatial distribution. The first-order complexity in aftershock distribution is clear and well beyond the influence of typical earthquake location uncertainty. The sequence included rupture of multiple fault segments, spatially separated. We use surface-wave based relative centroid locations to examine whether, at the small scale, the distribution of the aftershocks was influenced by location errors. Surface-wave based relative location has delineated precise oceanic transform fault earthquake locations in multiple regions. However, the relocated aftershocks off the coast of Sumatra seldom align along simple linear trends that are compatible with the corresponding fault strikes as estimated for the GCMT catalog. The relocation of roughly 60 moderate-earthquake epicentroids suggests that the faulting involved in the 2012 earthquake aftershock sequence included strain release along many short fault segments. Statistical analysis and temporal variations of aftershocks show a typical decay of the aftershocks but a relatively low number of aftershocks, as is common for intraplate oceanic earthquakes. Coulomb stress calculations indicate that most of the moderate-magnitude aftershocks are compatible with stress changes predicted by the large-event slip models. Furthermore the patterns in the aftershocks suggest that the formation of the boundary and eventual localization of deformation between the Indian and Australian plate is a complicated process.

2012 Sumatra Earthquake↗

FOCUS: an alignment-free model to identify organisms in metagenomes using non-negative least squares

One of the major goals in metagenomics is to identify the organisms present in a microbial community from unannotated shotgun sequencing reads. Taxonomic profiling has valuable applications in biological and medical research, including disease diagnostics. Most currently available approaches do not scale well with increasing data volumes, which is important because both the number and lengths of the reads provided by sequencing platforms keep increasing. Here we introduce FOCUS, an agile composition based approach using non-negative least squares (NNLS) to report the organisms present in metagenomic samples and profile their abundances. FOCUS was tested with simulated and real metagenomes, and the results show that our approach accurately predicts the organisms present in microbial communities. FOCUS was implemented in Python. The source code and web-sever are freely available at http://edwards.sdsu.edu/FOCUS.

59 BASIC BIOLOGICAL SCIENCES↗

IDENTIFICATION OF NICOTINAMIDE MONONUCLEOTIDE DEAMIDASE OF THE BACTERIAL PYRIDINE NUCLEOTIDE CYCLE REVEALS A NOVEL BROADLY CONSERVED AMIDOHYDROLASE FAMILY

The pyridine nucleotide cycle (PNC) is a network of salvage and recycling routes maintaining homeostasis of NAD(P) cofactor pool in the cell. Nicotinamide mononucleotide (NMN) deamidase (EC 3.5.1.42), one of the key enzymes of the bacterial PNC was originally described in Enterobacteria, but the corresponding gene eluded identification for over 30 years. A genomics-based reconstruction of NAD metabolism across hundreds bacterial species suggested that NMN deamidase reaction is the only possible way of nicotinamide salvage in the marine bacterium Shewanella oneidensis. This prediction was verified via purification of native NMN deamidase from S. oneidensis followed by the identification of the respective gene, termed pncC. Enzymatic characterization of the PncC protein, as well as phenotype analysis of deletion mutants, confirmed its proposed biochemical and physiological function in S. oneidensis. Of the three PncC homologs present in E. coli, NMN deamidase activity was confirmed only for the recombinant purified product of the ygaD gene. A comparative analysis at the level of sequence and three dimensional structure, which is available for one of the PncC family member, shows no homology with any previously described amidohydrolases. Multiple alignment analysis of functional and non functional PncC homologs, together with NMN docking experiments, allowed us to tentatively identify the active site area and conserved residues therein. An observed broad phylogenomic distribution of predicted functional PncCs in bacterial kingdom is consistent with a possible role in detoxification of NMN, resulting from NAD utilization by DNA ligase.

59 BASIC BIOLOGICAL SCIENCES↗

Structure of the Core Postfusion Porcine Endogenous Retrovirus Fusion Protein

Retroviral elements from endogenous retroviruses have functions in mammalian physiology. The best-known examples are the envelope proteins that function in placenta development and immune suppression. Porcine endogenous retroviruses (PERVs) are an understudied class of endogenous retroviruses that infect cultured human cells, raising concern regarding porcine xenografts. The PERV envelope glycoprotein has also been proposed as a possible swine syncytin with a role in placental development. Despite the growing interest in PERVs, their envelope glycoproteins remain poorly characterized. Here, we successfully determined the postfusion crystal structure of the PERV core fusion ectodomain. The PERV fusion protein structure reveals a conserved class I viral fusion protein six-helix bundle. Biophysical experiments demonstrated that the thermodynamic stability of the PERV fusion protein secondary structure was the same at physiological and acidic pHs. A conserved surface analysis highlights the high degree of sequence conservation among retroviral fusogens in the chain reversal region that facilitates the large-scale conformational change required for membrane fusion. Further structural alignment of class I viral fusogens revealed a phylogenetic clustering that shows evolution into various lineages that correlate with virus mechanisms of cell entry. Our work indicates that structural dendrograms can be used to qualitatively infer insights into the fusion mechanisms of newly discovered class I viral fusogen structures.

59 BASIC BIOLOGICAL SCIENCES↗

Strand-Specific RNA-Seq Analyses of Fruiting Body Development in Coprinopsis cinerea

We report that the basidiomycete fungus Coprinopsis cinerea is an important model system for multicellular development. Fruiting bodies of C. cinerea are typical mushrooms, which can be produced synchronously on defined media in the laboratory. To investigate the transcriptome in detail during fruiting body development, high-throughput sequencing (RNA-seq) was performed using cDNA libraries strand-specifically constructed from 13 points (stages/tissues) with two biological replicates. The reads were aligned to 14,245 predicted transcripts, and counted for forward and reverse transcripts. Differentially expressed genes (DEGs) between two adjacent points and between vegetative mycelium and each point were detected by Tag Count Comparison (TCC). To validate RNA-seq data, expression levels of selected genes were compared using RPKM values in RNA-seq data and qRT-PCR data, and DEGs detected in microarray data were examined in MA plots of RNA-seq data by TCC. We discuss events deduced from GO analysis of DEGs. In addition, we uncovered both transcription factor candidates and antisense transcripts that are likely to be involved in developmental regulation for fruiting.

59 BASIC BIOLOGICAL SCIENCES↗

Recruiting Human Microbiome Shotgun Data to Site-Specific Reference Genomes

The human body consists of innumerable multifaceted environments that predispose colonization by a number of distinct microbial communities, which play fundamental roles in human health and disease. In addition to community surveys and shotgun metagenomes that seek to explore the composition and diversity of these microbiomes, there are significant efforts to sequence reference microbial genomes from many body sites of healthy adults. To illustrate the utility of reference genomes when studying more complex metagenomes, we present a reference-based analysis of sequence reads generated from 55 shotgun metagenomes, selected from 5 major body sites, including 16 sub-sites. Interestingly, between 13% and 92% (62.3% average) of these shotgun reads were aligned to a then-complete list of 2780 reference genomes, including 1583 references for the human microbiome. However, no reference genome was universally found in all body sites. For any given metagenome, the body site-specific reference genomes, derived from the same body site as the sample, accounted for an average of 58.8% of the mapped reads. While different body sites did differ in abundant genera, proximal or symmetrical body sites were found to be most similar to one another. The extent of variation observed, both between individuals sampled within the same microenvironment, or at the same site within the same individual over time, calls into question comparative studies across individuals even if sampled at the same body site. This study illustrates the high utility of reference genomes and the need for further site-specific reference microbial genome sequencing, even within the already well-sampled human microbiome.

59 BASIC BIOLOGICAL SCIENCES↗

pGraph: Efficient Parallel Construction of Large-Scale Protein Sequence Homology Graphs

Detecting sequence homology between protein sequences is a fundamental problem in computational molecular biology, with a pervasive application in nearly all analyses that aim to structurally and functionally characterize protein molecules. While detecting the homology between two protein sequences is relatively inexpensive, detecting pairwise homology for a large number of protein sequences can become computationally prohibitive for modern inputs, often requiring millions of CPU hours. Yet, there is currently no robust support to parallelize this kernel. In this paper, we identify the key characteristics that make this problemparticularly hard to parallelize, and then propose a new parallel algorithm that is suited for detecting homology on large data sets using distributed memory parallel computers. Our method, called pGraph, is a novel hybrid between the hierarchical multiple-master/worker model and producer-consumer model, and is designed to break the irregularities imposed by alignment computation and work generation. Experimental results show that pGraph achieves linear scaling on a 2,048 processor distributed memory cluster for a wide range of inputs ranging from as small as 20,000 sequences to 2,560,000 sequences. In addition to demonstrating strong scaling, we present an extensive report on the performance of the various system components and related parametric studies.

Environmental Molecular Sciences Laboratory↗

SCALE depletion capabilities for molten salt reactors and other liquid-fueled systems

Nuclear reactor systems that use fuel dissolved in a liquid have the potential for enhanced safety characteristics, improved fuel-cycle outcomes, and more efficient isotope-production configurations. In these reactor systems, the fueled liquid may simultaneously undergo irradiation, physical and chemical removal processes, and fueling. The modeling and simulation of this transmutation and decay with material additions and removals is an ongoing research area. An accurate simulation tool is critical to the reactor and fuel-cycle design, reactor deployment, and source-term characterization for these advanced reactor systems. The work described herein involved implementing, testing, and applying the capability to perform reactor physics simulations within the Oak Ridge National Laboratory-developed SCALE suite for nuclear systems analyses and design, leveraging much of its pedigree in quality-assurance and reactor-analysis capabilities. The functionalities to simulate irradiation with material feeds and removals had been added in ORIGEN, and the TRITON reactor physics sequence was extended to calculate the total removed material and track external nonirradiated mixtures to estimate separate processing or waste streams. Results from these capabilities align with analytical expectations obtained from ORIGEN for simplified test cases and with expectations for a molten salt reactor application. This implementation, available with the SCALE 6.3 release, provides for a more efficient and accurate material accountability methodology, allowing for the characterization, design, and analysis of the complete isotopic material inventory of advanced liquid-fueled systems for a variety of applications.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Identifying DQ-Domain Admittance Models of a 2.3-MVA Commercial Grid-Following Inverter Via Frequency-Domain and Time-Domain Data

Here, we present two methods to identify the admittance of a 2.3-MVA commercial grid-following inverter through frequency-domain data and time-domain data. In addition to the well-known harmonic injection method to obtain the admittance at frequency points followed by vector fitting, we adopt a step response-based method where the step response data of the converter are collected and their s-domain expressions are obtained. In turn, the s-domain admittance model of the converter is found. The two methods can be used for cross validation. Step responses and frequency-domain responses obtained from the two methods are compared. Results show that the time- domain data-based method leads to an admittance comparable with that from the frequency-domain measurement in the range less than 60 Hz. In addition, two insights are obtained from the frequency-domain admittance measurements. First, when the dq- frame is aligned to the point of common coupling (PCC) voltage, the per unit DQ-domain admittance directly reflects operating conditions. Second, existence of negative-sequence control can be detected via sequence-domain admittances.

24 POWER TRANSMISSION AND DISTRIBUTION↗