Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 505 records · Page 28

A NASTRAN DMAP alter for linear buckling analysis under dynamic loading

A unique modification to the NASTRAN solution sequence for transient analysis with direct time integration (COSMIC NASTRAN rigid format 9) was developed and incorporated into a DMAP alter. This DMAP alter calculates the buckling stability of a dynamically loaded structure, and is used to predict the onset of structural buckling under stress wave loading conditions. The modified solution sequence incorporates the linear buckling analysis capability (rigid format 5) of NASTRAN into the existing Transient solution rigid format in such a way as to provide a time dependent eigensolution which is used to assess the buckling stability of the structure as it responds to the impulsive load. As a demonstration of the validity of this modified solution procedure, the dynamic buckling of a prismatic bar subjected to an impulsive longitudinal compression is analyzed and compared to the known theoretical solution. In addition, a dynamic buckling analysis is performed for the analytically less tractable problem of the localized dynamic buckling of an initially flawed composite laminate under transverse impact loading. The addition of this DMAP alter to the transient solution sequence in NASTRAN facilitates the prediction of both time and mode of buckling.

Aiello, Robert A.↗

Scaling features of noncoding DNA

We review 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, and utilize this fact to build a Coding Sequence Finder Algorithm, which uses statistical ideas to locate the coding regions of an unknown DNA sequence. Finally, we describe briefly some recent work adapting 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, and reporting that noncoding regions in eukaryotes display a larger redundancy than coding regions. Specifically, we consider the possibility that this result is solely a consequence of nucleotide concentration differences as first noted by Bonhoeffer and his collaborators. We find that cytosine-guanine (CG) concentration does have a strong "background" effect on redundancy. However, we find that for the purine-pyrimidine binary mapping rule, which is not affected by the difference in CG concentration, the Shannon redundancy for the set of analyzed sequences is larger for noncoding regions compared to coding regions.

Non-NASA Center↗

The “Most Wanted” Taxa from the Human Microbiome for Whole Genome Sequencing

The goal of the Human Microbiome Project (HMP) is to generate a comprehensive catalog of human-associated microorganisms including reference genomes representing the most common species. Toward this goal, the HMP has characterized the microbial communities at 18 body habitats in a cohort of over 200 healthy volunteers using 16S rRNA gene (16S) sequencing and has generated nearly 1,000 reference genomes from human-associated microorganisms. To determine how well current reference genome collections capture the diversity observed among the healthy microbiome and to guide isolation and future sequencing of microbiome members, we compared the HMP’s 16S data sets to several reference 16S collections to create a ‘most wanted’ list of taxa for sequencing. Our analysis revealed that the diversity of commonly occurring taxa within the HMP cohort microbiome is relatively modest, few novel taxa are represented by these OTUs and many common taxa among HMP volunteers recur across different populations of healthy humans. Taken together, these results suggest that it should be possible to perform whole-genome sequencing on a large fraction of the human microbiome, including the ‘most wanted’, and that these sequences should serve to support microbiome studies across multiple cohorts. Also, in stark contrast to other taxa, the ‘most wanted’ organisms are poorly represented among culture collections suggesting that novel culture- and single-cell-based methods will be required to isolate these organisms for sequencing.

59 BASIC BIOLOGICAL SCIENCES↗

Mechanism of H 2 S Oxidation by the Dissimilatory Perchlorate-Reducing Microorganism Azospira suillum PS

The genetic and biochemical basis of perchlorate-dependent H 2 S oxidation (PSOX) was investigated in the dissimilatory perchlorate-reducing microorganism (DPRM) Azospira suillum PS (PS). Previously, it was shown that all known DPRMs innately oxidize H 2 S, producing elemental sulfur (S o ). Although the process involving PSOX is thermodynamically favorable (ΔG ° ' = -206 kJ ∙ mol -1 H 2 S), the underlying biochemical and genetic mechanisms are currently unknown. Interestingly, H 2 S is preferentially utilized over physiological electron donors such as lactate or acetate although no growth benefit is obtained from the metabolism. Here, we determined that PSOX is due to a combination of enzymatic and abiotic interactions involving reactive intermediates of perchlorate respiration. Using various approaches, including barcode analysis by sequencing (Bar-seq), transcriptome sequencing (RNA-seq), and proteomics, along with targeted mutagenesis and biochemical characterization, we identified all facets of PSOX in PS. In support of our proposed model, deletion of identified upregulated PS genes traditionally known to be involved in sulfur redox cycling (e.g., Sox, sulfide:quinone reductase [SQR]) showed no defect in PSOX activity. Proteomic analysis revealed differential abundances of a variety of stress response metal efflux pumps and divalent heavy-metal transporter proteins, suggesting a general toxicity response. Furthermore, in vitro biochemical studies demonstrated direct PSOX mediated by purified perchlorate reductase (PcrAB) in the absence of other electron transfer proteins. The results of these studies support a model in which H 2 S oxidation is mediated by electron transport chain short-circuiting in the periplasmic space where the PcrAB directly oxidizes H 2 S to S o . The biogenically formed reactive intermediates (ClO 2 - and O 2 ) subsequently react with additional H 2 S, producing polysulfide and S o as end products.

59 BASIC BIOLOGICAL SCIENCES↗

Direct profiling of environmental microbial populations by thermal dissociation analysis of native rRNAs hybridized to oligonucleotide microarrays

Oligonucleotide microarrays were used to profile directly extracted rRNA from environmental microbial populations without PCR amplification. In our initial inspection of two distinct estuarine study sites, the hybridization patterns were reproducible and varied between estuarine sediments of differing salinities. The determination of a thermal dissociation curve (i.e., melting profile) for each probe-target duplex provided information on hybridization specificity, which is essential for confirming adequate discrimination between target and nontarget sequences.

Non-NASA Center↗

Timing and source of subtype-C HIV-1 superinfection in the newly infected partner of Zambian couples with disparate viruses

Background: HIV-1 superinfection occurs at varying frequencies in different at risk populations. Though seroincidence is decreased, in the negative partner of HIV-discordant couples after joint testing and counseling in the Zambia Emory HIV Research Project (ZEHRP) cohort, the annual infection rate remains relatively high at 7-8%. Based on sequencing within the gp41 region of each partner’s virus, 24% of new infections between 2004 and 2008 were the result of transmission from a non-spousal partner. Since these seroconvertors and their spouses have disparate epidemiologically-unlinked viruses, there is a risk of superinfection within the marriage. We have, therefore, investigated the incidence and viral origin of superinfection in these couples. Results: Superinfection was detected by heteroduplex mobility assay (HMA), degenerate base counting of the gp41 sequence, or by phylogenetic analysis of the longitudinal sequences. It was confirmed by full-length env single genome amplification and phylogenetic analysis. In 22 couples (44 individuals), followed for up to five years, three of the newly infected (initially HIV uninfected) partners became superinfected. In each case superinfection occurred during the first 12 months following initial infection of the negative partner, and in each case the superinfecting virus was derived from a non-spousal partner. In addition, one probable case of intra-couple HIV-1 superinfection was observed in a chronically infected partner at the time of his seroconverting spouse’s initial viremia. Extensive recombination within the env gene was observed following superinfection. Conclusions: In this subtype-C discordant couple cohort, superinfection, during the first year after HIV-1 infection of the previously negative partner, occurred at a rate similar to primary infection (13.6% [95% CI 5.2-34.8] vs 7.8% [7.1-8.6]). While limited intra-couple superinfection may in part reflect continued condom usage within couples, this and our lack of detecting newly superinfected individuals after one year of primary infection raise the possibility that immunological resistance to intra-subtype superinfection may develop over time in subtype C infected individuals.

59 BASIC BIOLOGICAL SCIENCES↗

Dynamic stress analysis of smooth and notched fiber composite flexural specimens

A detailed analysis of the dynamic stress field in smooth and notched fiber composite (Charpy-type) specimens is reported in this paper. The analysis is performed with the aid of the direct transient response analysis solution sequence of MSC/NASTRAN. Three unidirectional composites were chosen for the study. They are S-Glass/Epoxy, Kevlar/Epoxy and T-300/Epoxy composite systems. The specimens are subjected to an impact load which is modeled as a triangular impulse with a maximum of 2000 lb and a duration of 1 ms. The results are compared with those of static analysis of the specimens subjected to a peak load of 2000 lb. For the geometry and type of materials studied, the static analysis results gave close conservative estimates for the dynamic stresses. Another interesting inference from the study is that the impact induced effects are felt by S-Glass/Epoxy specimens sooner than Kevlar/Epoxy or T-300/Epoxy specimens.

Murthy, P. L. N.↗

Dynamic stress analysis of smooth and notched fiber composite flexural specimens

A detailed analysis of the dynamic stress field in smooth and notched fiber composite (Charpy-type) specimens is reported in this paper. The analysis is performed with the aid of the direct transient response analysis solution sequence of MSC/NASTRAN. Three unidirectional composites were chosen for the study. They are S-Glass/Epoxy, Kevlar/Epoxy and T-300/Epoxy composite systems. The specimens are subjected to an impact load which is modeled as a triangular impulse with a maximum of 2000 lb and a duration of 1 ms. The results are compared with those of static analysis of the specimens subjected to a peak load of 2000 lb. For the geometry and type of materials studied, the static analysis results gave close conservative estimates for the dynamic stresses. Another interesting inference from the study is that the impact induced effects are felt by S-Glass/Epoxy specimens sooner than Kevlar/Epoxy or T-300/Epoxy specimens.

Murthy, P. L. N.↗

Selective Whole-Genome Amplification as a Tool to Enrich Specimens with Low Treponema pallidum Genomic DNA Copies for Whole-Genome Sequencing

Downstream next-generation sequencing (NGS) of the syphilis spirochete Treponema pallidum subspecies pallidum (T. pallidum) is hindered by low bacterial loads and the overwhelming presence of background metagenomic DNA in clinical specimens. In this study, we investigated selective whole-genome amplification (SWGA) utilizing multiple displacement amplification (MDA) in conjunction with custom oligonucleotides with an increased specificity for the T. pallidum genome and the capture and removal of 5'-C-phosphate-G-3' (CpG) methylated host DNA using the NEBNext Microbiome DNA enrichment kit followed by MDA with the REPLI-g single cell kit as enrichment methods to improve the yields of T. pallidum DNA in isolates and lesion specimens from syphilis patients. Sequencing was performed using the Illumina MiSeq v2 500 cycle or NovaSeq 6000 SP platform. These two enrichment methods led to 93 to 98% genome coverage at 5 reads/site in 5 clinical specimens from the United States and rabbit-propagated isolates, containing >14 T. pallidum genomic copies/μL of sample for SWGA and >129 genomic copies/μL for CpG methylation capture with MDA. Variant analysis using sequencing data derived from SWGA-enriched specimens showed that all 5 clinical strains had the A2058G mutation associated with azithromycin resistance. SWGA is a robust method that allows direct whole-genome sequencing (WGS) of specimens containing very low numbers of T. pallidum, which has been challenging until now.

59 BASIC BIOLOGICAL SCIENCES↗

Multiple Lines of Evidences Reveal Mechanisms Underpinning Mercury Resistance and Volatilization by Stenotrophomonas sp. MA5 Isolated from the Savannah River Site (SRS), USA

A largely understudied microbially mediated mercury (Hg) bioremediative pathway includes the volatilization of Hg 2+ to Hg 0 . Therefore, studies on Hg resistant bacteria (HgR), isolated from historically long-term contaminated environments, can serve as models to understand mechanisms underpinning Hg cycling. Towards this end, a mercury resistant bacterial strain, identified as Stenotrophomonas sp., strain MA5, was isolated from Mill Branch on the Savannah River Site (SRS); an Hg-impacted ecosystem. Minimum inhibitory concentration (MIC) analysis showed Hg resistance of up to 20 µg/mL by MA5 with 95% of cells retaining viability. Microcosm studies showed that the strain depleted more than 90% of spiked Hg 2+ within the first 24 h of growth and the detection of volatilized mercury indicated that the strain was able to reduce Hg 2+ to Hg 0 . To understand molecular mechanisms of Hg volatilization, a draft whole genome sequence was obtained, annotated and analyzed, which revealed the presence of a transposon-derived mer operon (merRTPADE) in MA5, known to transport and reduce Hg 2+ into Hg 0 . Based on the whole genome sequence of strain MA5, qRT-PCR assays were designed on merRTPADE, we found a ~40-fold higher transcription of mer T, P, A, D and E when cells were exposed to 5 µg/mL Hg 2+ . Interestingly, strain MA5 increased cellular size as a function of increasing Hg concentrations, which is likely an evolutionary response mechanism to cope with Hg stress. Moreover, metal contaminated environments are shown to co-select for antibiotic resistance. When MA5 was screened for antibiotic resistance, broad resistance against penicillin, streptomycin, tetracycline, ampicillin, rifampicin, and erythromycin was found; this correlated with the presence of multiple gene determinants for antibiotic resistance within the whole genome sequence of MA5. Overall, this study provides an in-depth understanding of the underpinnings of Stenotrophomonas-mercury interactions that facilitate cellular survival in a contaminated soil habitat.

59 BASIC BIOLOGICAL SCIENCES↗

Deletion analysis of the itaconic acid biosynthesis gene cluster components in Aspergillus pseudoterreus ATCC32359

The filamentous fungus Aspergillus terreus has been successfully used for industrial production of itaconic acid (IA) for many years. The IA biosynthesis pathway has recently been characterized at a molecular genetic level as an IA gene cluster (IA cluster) by a clone-based transcriptomic approach. The cluster consists of four genes, including genes for cis-aconitic acid decarboxylase (cadA), a predicted transcription factor (tf), mitochondrial organic acid transporter (mttA), and MFS (Major Facilitator Superfamily) type transporter (mfsA). In this research, we performed EST (Expressed Sequence Tag) analysis and systematic gene deletions to investigate the role of those genes in IA biosynthesis in A. pseudoterreus ATCC32359. EST analysis showed a similar expression pattern among those four genes that was distinct from neighbouring genes and further confirmed that they belong to the same biosynthesis cluster. Systematic gene deletion analysis demonstrated that tf, cadA, mttA, and mfsA genes in the cluster are essential for IA production. Deletion of any of them will either completely abolish the IA production or dramatically decrease amount of IA produced. More importantly, a significant amount of aconitic acid was detected in the cadA deletion strain but not in the other deletion strains. Therefore, the cadA deletion strain is a novel microorganism host for production of aconitic acid and other value-added chemicals from lignocellulosic biomass.

Deng, Shuang↗

ClassyFlu: Classification of Influenza A Viruses with Discriminatively Trained Profile-HMMs

Accurate and rapid characterization of influenza A virus (IAV) hemagglutinin (HA) and neuraminidase (NA) sequences with respect to subtype and clade is at the basis of extended diagnostic services and implicit to molecular epidemiologic studies. ClassyFlu is a new tool and web service for the classification of IAV sequences of the HA and NA gene into subtypes and phylogenetic clades using discriminatively trained profile hidden Markov models (HMMs), one for each subtype or clade. ClassyFlu merely requires as input unaligned, full-length or partial HA or NA DNA sequences. It enables rapid and highly accurate assignment of HA sequences to subtypes H1–H17 but particularly focusses on the finer grained assignment of sequences of highly pathogenic avian influenza viruses of subtype H5N1 according to the cladistics proposed by the H5N1 Evolution Working Group. NA sequences are classified into subtypes N1–N10. ClassyFlu was compared to semiautomatic classification approaches using BLAST and phylogenetics and additionally for H5 sequences to the new ‘‘Highly Pathogenic H5N1 Clade Classification Tool’’ (IRD-CT) proposed by the Influenza Research Database. Our results show that both web tools (ClassyFlu and IRD-CT), although based on different methods, are nearly equivalent in performance and both are more accurate and faster than semiautomatic classification. A retraining of ClassyFlu to altered cladistics as well as an extension of ClassyFlu to other IAV genome segments or fragments thereof is undemanding. This is exemplified by unambiguous assignment to a distinct cluster within subtype H7 of sequences of H7N9 viruses which emerged in China early in 2013 and caused more than 130 human infections.

59 BASIC BIOLOGICAL SCIENCES↗

IMG-ABC: An Atlas of Biosynthetic Gene Clusters to Fuel the Discovery of Novel Secondary Metabolites

In the discovery of secondary metabolites (SMs), large-scale analysis of sequence data is a promising exploration path that remains largely underutilized due to the lack of relevant computational resources. We present IMG-ABC (https://img.jgi.doe.gov/abc/) -- An Atlas of Biosynthetic gene Clusters within the Integrated Microbial Genomes (IMG) system1. IMG-ABC is a rich repository of both validated and predicted biosynthetic clusters (BCs) in cultured isolates, single-cells and metagenomes linked with the SM chemicals they produce and enhanced with focused analysis tools within IMG. The underlying scalable framework enables traversal of phylogenetic dark matter and chemical structure space -- serving as a doorway to a new era in the discovery of novel molecules.

99 GENERAL AND MISCELLANEOUS↗

Metatranscriptomic insights on gene expression and regulatory controls in Candidatus Accumulibacter phosphatis

Previous studies on enhanced biological phosphorus removal (EBPR) have focused on reconstructing genomic blueprints for the model polyphosphate-accumulating organism Candidatus Accumulibacter phosphatis. Here, a time series metatranscriptome generated from enrichment cultures of Accumulibacter was used to gain insight into anerobic/aerobic metabolism and regulatory mechanisms within an EBPR cycle. Co-expressed gene clusters were identified displaying ecologically relevant trends consistent with batch cycle phases. Transcripts displaying increased abundance during anerobic acetate contact were functionally enriched in energy production and conversion, including upregulation of both cytoplasmic and membrane-bound hydrogenases demonstrating the importance of transcriptional regulation to manage energy and electron flux during anerobic acetate contact. We hypothesized and demonstrated hydrogen production after anerobic acetate contact, a previously unknown strategy for Accumulibacter to maintain redox balance. Genes involved in anerobic glycine utilization were identified and phosphorus release after anerobic glycine contact demonstrated, suggesting that Accumulibacter routes diverse carbon sources to acetyl-CoA formation via previously unrecognized pathways. A comparative genomics analysis of sequences upstream of co-expressed genes identified two statistically significant putative regulatory motifs. One palindromic motif was identified upstream of genes involved in PHA synthesis and acetate activation and is hypothesized to be a phaR binding site, hence representing a hypothetical PHA modulon. A second motif was identified ~35 base pairs (bp) upstream of a large and diverse array of genes and hence may represent a sigma factor binding site. As a result, this analysis provides a basis and framework for further investigations into Accumulibacter metabolism and the reconstruction of regulatory networks in uncultured organisms.

54 ENVIRONMENTAL SCIENCES↗

Combined theoretical and experimental study of the valence, Rydberg and ionic states of fluorobenzene

New photoelectron spectra (PES) and ultra violet (UV) and vacuum UV (VUV) absorption spectra of fluorobenzene recorded at higher resolution than previously, have been combined with mass-resolved (2 + 1) and (3 + 1) resonance enhanced multiphoton ionization (REMPI) spectra; this has led to the identification of numerous Rydberg states. The PES have been compared with earlier mass-analyzed threshold ionization and photoinduced Rydberg ionization (PIRI) spectra to give an overall picture of the ionic state sequence. The analysis of these spectra using both equations of motion with coupled cluster singles and doubles (EOM-CCSD) configuration interaction and time dependent density functional theory (TDDFT) calculations have been combined with vibrational analysis of both the hot and cold bands of the spectra, in considerable detail. The results extend several earlier studies on the vibronic coupling leading to conical intersections between the X{sup 2}B{sub 1} and A{sup 2}A{sub 2} states, and a further trio (B, C, and D) of states. The conical intersection of the X and A states has been explicitly identified, and its structure and energetics evaluated. The energy sequence of the last group is only acceptable to the present study if given as B{sup 2}B{sub 2}<C{sup 2}B{sub 1}<D{sup 2}A{sub 1}, a conclusion which is in agreement with most previous EOM-CCSD and other calculations. However, this symmetry ordering of the B and C states forces reconsideration of the nature of the PIRI spectrum. The coupling between these two states is induced by the a{sub 2} modes, ν{sub 12} and ν{sub 14} and we propose that the 14{sup 1} band is observed in the B{sup 2}B{sub 2} band in the PES for the first time, because of the improved resolution. This same assignment is given to the lowest energy band in the PIRI spectrum which was previously assigned as the origin band and further conclude that the entire PIRI spectrum is induced by ν{sub 12} and ν{sub 14}. The relative intensities of the various Rydberg state peaks in the VUV absorption and REMPI spectra of fluorobenzene are very similar to those observed in the equivalent spectra of benzene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A deep learning approach to real-time HIV outbreak detection using genetic data

Pathogen genomic sequence data are increasingly made available for epidemiological monitoring. A main interest is to identify and assess the potential of infectious disease outbreaks. While popular methods to analyze sequence data often involve phylogenetic tree inference, they are vulnerable to errors from recombination and impose a high computational cost, making it difficult to obtain real-time results when the number of sequences is in or above the thousands. Here, we propose an alternative strategy to outbreak detection using genomic data based on deep learning methods developed for image classification. The key idea is to use a pairwise genetic distance matrix calculated from viral sequences as an image, and develop convolutional neutral network (CNN) models to classify areas of the images that show signatures of active outbreak, leading to identification of subsets of sequences taken from an active outbreak. We showed that our method is efficient in finding HIV-1 outbreaks with R0 ≥ 2.5, and overall a specificity exceeding 98% and sensitivity better than 92%. We validated our approach using data from HIV-1 CRF01 in Europe, containing both endemic sequences and a well-known dual outbreak in intravenous drug users. Our model accurately identified known outbreak sequences in the background of slower spreading HIV. Importantly, we detected both outbreaks early on, before they were over, implying that had this method been applied in real-time as data became available, one would have been able to intervene and possibly prevent the extent of these outbreaks. This approach is scalable to processing hundreds of thousands of sequences, making it useful for current and future real-time epidemiological investigations, including public health monitoring using large databases and especially for rapid outbreak identification.

59 BASIC BIOLOGICAL SCIENCES↗

Multiple melanocortin receptors are expressed in bone cells

Melanocortin receptors belong to the seven transmembrane domain, G-protein coupled family of receptors. There are five members of this receptor family labeled MC1R-MC5R. These receptors are activated by fragments derived from a larger molecule, proopiomelanocortin (POMC) and include ACTH, alpha beta and gamma-MSH and beta-endorphin. Because of in vitro and in vivo data suggesting direct effects of these POMC molecules on bone and bone turnover, we examined bone and bone derived cells for the presence of the various members of the melanocortin receptor family. We report that the five known melanocortin receptors are expressed to varying degrees in osteoblast-like and osteoclastic cells. POMC fragments increased proliferation and expression of a variety of genes in osteoblastic cells. Furthermore, POMC mRNA was detected in osteoclastic cells. These data demonstrate that POMC-derived peptide hormones acting through high affinity melanocortin receptors have specific effects on bone cells. Thus, in addition to the indirect effects of POMC-derived hormones on bone turnover through their modulation of steroid hormone secretion, POMC fragments may have direct and specific effects on bone cell subpopulations.

NASA Discipline Musculoskeletal↗