Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sequence comparison”

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 19 records

Benchmarking of alignment-free sequence comparison methods

Background: Alignment-free (AF) sequence comparison is attracting persistent interest driven by data-intensive applications. Hence, many AF procedures have been proposed in recent years, but a lack of a clearly defined benchmarking consensus hampers their performance assessment. Results: Here, we present a community resource (http://afproject.org) to establish standards for comparing alignment-free approaches across different areas of sequence-based research. We characterize 74 AF methods available in 24 software tools for five research applications, namely, protein sequence classification, gene tree inference, regulatory element detection, genome-based phylogenetic inference, and reconstruction of species trees under horizontal gene transfer and recombination events. Conclusion: The interactive web service allows researchers to explore the performance of alignment-free tools relevant to their data types and analytical goals. It also allows method developers to assess their own algorithms and compare them with current state-of-the-art tools, accelerating the development of new, more accurate AF solutions

59 BASIC BIOLOGICAL SCIENCES↗

Comparison of Sequence Component-Based Fault Detection and Relay Coordination Algorithms in Inverter-Based Networks

Protection of inverter-based microgrids using sequence component-based relaying schemes is a promising solution. These methods offer several advantages, including lower computational requirements, compatibility with commercial relay systems, and cost-effectiveness compared to communication-based approaches. This article investigate the performance of various sequence component based schemes with the objective of identifying the algorithms that provide the best fault detection and relay coordination, solely relying on local voltages and current at relay terminals. Positive, negative and zero sequence impedance, admittance and power detection algorithms were tested on modified IEEE 13 bus test network for various shunt faults (LG, LL, LLG, LLL). Hardware-in-the-loop validation was achieved using the Typhoon real-time simulator, interfacing with a SEL 751 relay. This research demonstrates that while several algorithms are capable of detecting faults with sufficient accuracy, only a few are effective in achieving proper coordination. Validation results indicate that the negative sequence power approach provides the best performance in both fault detection and coordination.

Patel, Deepika [ORNL] (ORCID:0000000341099994)↗

Atypical Divergence of SARS-CoV-2 Orf8 from Orf7a within the Coronavirus Lineage Suggests Potential Stealthy Viral Strategies in Immune Evasion

Orf8, one of the most puzzling genes in the SARS lineage of coronaviruses, marks a unique and striking difference in genome organization between SARS-CoV-2 and SARS-CoV-1. Here, using sequence comparisons, we unequivocally reveal the distant sequence similarities between SARS-CoV-2 Orf8 with its SARS-CoV-1 counterparts and the X4-like genes of coronaviruses, including its highly divergent “paralog” gene Orf7a, whose product is a potential immune antagonist of known structure. Supervised sequence space walks unravel identity levels that drop below 10% and yet exhibit subtle conservation patterns in this novel superfamily, characterized by an immunoglobulin-like beta sandwich topology. We document the high accuracy of the sequence space walk process in detail and characterize the subgroups of the superfamily in sequence space by systematic annotation of gene and taxon groups. While SARS-CoV-1 Orf7a and Orf8 genes are most similar to bat virus sequences, their SARS-CoV-2 counterparts are closer to pangolin virus homologs, reflecting the fine structure of conservation patterns within the SARS-CoV-2 genomes. The divergence between Orf7a and Orf8 is exceptionally idiosyncratic, since Orf7a is more constrained, whereas Orf8 is subject to rampant change, a peculiar feature that may be related to hitherto-unknown viral infection strategies. Despite their common origin, the Orf7a and Orf8 protein families exhibit different modes of evolutionary trajectories within the coronavirus lineage, which might be partly attributable to their complex interactions with the mammalian host cell, reflected by a multitude of functional associations of Orf8 in SARS-CoV-2 compared to a very small number of interactions discovered for Orf7a.

60 APPLIED LIFE SCIENCES↗

Sequence information signal processor for local and global string comparisons

A sequence information signal processing integrated circuit chip designed to perform high speed calculation of a dynamic programming algorithm based upon the algorithm defined by Waterman and Smith. The signal processing chip of the present invention is designed to be a building block of a linear systolic array, the performance of which can be increased by connecting additional sequence information signal processing chips to the array. The chip provides a high speed, low cost linear array processor that can locate highly similar global sequences or segments thereof such as contiguous subsequences from two different DNA or protein sequences. The chip is implemented in a preferred embodiment using CMOS VLSI technology to provide the equivalent of about 400,000 transistors or 100,000 gates. Each chip provides 16 processing elements, and is designed to provide 16 bit, two's compliment operation for maximum score precision of between -32,768 and +32,767. It is designed to provide a comparison between sequences as long as 4,194,304 elements without external software and between sequences of unlimited numbers of elements with the aid of external software. Each sequence can be assigned different deletion and insertion weight functions. Each processor is provided with a similarity measure device which is independently variable. Thus, each processor can contribute to maximum value score calculation using a different similarity measure.

Peterson, John C.↗

Calmodulin Gene Family in Potato: Developmental and Touch-Induced Expression of the mRNA Encoding a Novel Isoform

Eight genomic clones of potato calmodulin (PCM1 to 8) were isolated and characterized. Sequence comparisons of different genes revealed that the deduced amino acid sequence of PCM1 had several unique substitutions, especially in the fourth Ca(2+)-binding area. The expression patterns of different genes were studied by northern analysis using the 3'-untranslated regions as probes. The expression of PCM1, 5, and 8 was highest in the stolon tip and it decreased during tuber development. The expression of PCM6 did not vary much in the tissues tested, except in the leaves, where the expression was lower; whereas, the expression of PCM4 was very low in all the tissues. The expression of PCM2 and PCM3 was not detected in any of the tissues tested. Among these genes, only PCM1 showed increased expression following touch stimulation. To study the regulation of PCM1, transgenic potato plants carrying the PCM1 promoter fused to the beta-glucuronidase (GUS) reporter gene were produced. GUS expression was found to be developmentally regulated and touch-responsive, indicating a positive correlation between the expression of PCM1 and GUS mRNAs. These results suggest that the 5'-flanking region of PCM1 controls developmental and touch-induced expression. X-Gluc staining patterns revealed that GUS localization is high in meristematic tissues such as the stem apex, stolon tip, and vascular regions.

Takezawa, D.↗

A General Framework to Learn Tertiary Structure for Protein Sequence Characterization

During the past five years, deep-learning algorithms have enabled ground-breaking progress towards the prediction of tertiary structure from a protein sequence. Very recently, we developed SAdLSA, a new computational algorithm for protein sequence comparison via deep-learning of protein structural alignments. SAdLSA shows significant improvement over established sequence alignment methods. In this contribution, we show that SAdLSA provides a general machine-learning framework for structurally characterizing protein sequences. By aligning a protein sequence against itself, SAdLSA generates a fold distogram for the input sequence, including challenging cases whose structural folds were not present in the training set. About 70% of the predicted distograms are statistically significant. Although at present the accuracy of the intra-sequence distogram predicted by SAdLSA self-alignment is not as good as deep-learning algorithms specifically trained for distogram prediction, it is remarkable that the prediction of single protein structures is encoded by an algorithm that learns ensembles of pairwise structural comparisons, without being explicitly trained to recognize individual structural folds. As such, SAdLSA can not only predict protein folds for individual sequences, but also detects subtle, yet significant, structural relationships between multiple protein sequences using the same deep-learning neural network. The former reduces to a special case in this general framework for protein sequence annotation.

59 BASIC BIOLOGICAL SCIENCES↗

Kinematics and spectra of planetary nebulae with O VI-sequence nuclei

Spectral features of NGC 5189 and NGC 6905 are tabulated. Fabry-Perot profiles around H alpha and O III lambda 5007 of NGC 5189, NGC 6905, NGC 246, and NGC 1535, are illustrated. The latter planetary nebula is a non-O VI-sequence, comparison object of high excitation. The kinematics of the four planetary nebulae are simply analyzed. Discussion of these data is motivated by the possibility of collisional excitation by high-speed ejecta from broad-lined O VI-sequence nuclei, and by the opportunity to make a comparison with conditions in the supernova remnant or ring nebula, G2.4 + 1.4, which contains an O VI-sequence nucleus of Population I.

Johnson, H. M.↗

FDA-ARGOS is a database with public quality-controlled reference genomes for diagnostic use and regulatory science

FDA proactively invests in tools to support innovation of emerging technologies, such as infectious disease next generation sequencing (ID-NGS). Here, we introduce FDA-ARGOS quality-controlled reference genomes as a public database for diagnostic purposes and demonstrate its utility on the example of two use cases. We provide quality control metrics for the FDA-ARGOS genomic database resource and outline the need for genome quality gap filling in the public domain. In the first use case, we show more accurate microbial identification of Enterococcus avium from metagenomic samples with FDA-ARGOS reference genomes compared to non-curated GenBank genomes. In the second use case, we demonstrate the utility of FDA-ARGOS reference genomes for Ebola virus target sequence comparison as part of a composite validation strategy for ID-NGS diagnostic tests. The use of FDA-ARGOS as an in silico target sequence comparator tool combined with representative clinical testing could reduce the burden for completing ID-NGS clinical trials

59 BASIC BIOLOGICAL SCIENCES↗

How close is close: 16S rRNA sequence identity may not be sufficient to guarantee species identity

16S rRNA (genes coding for rRNA) sequence comparisons were conducted with the following three psychrophilic strains: Bacillus globisporus W25T (T = type strain) and Bacillus psychrophilus W16AT, and W5. These strains exhibited more than 99.5% sequence identity and within experimental uncertainty could be regarded as identical. Their close taxonomic relationship was further documented by phenotypic similarities. In contrast, previously published DNA-DNA hybridization results have convincingly established that these strains do not belong to the same species if current standards are used. These results emphasize the important point that effective identity of 16S rRNA sequences is not necessarily a sufficient criterion to guarantee species identity. Thus, although 16S rRNA sequences can be used routinely to distinguish and establish relationships between genera and well-resolved species, very recently diverged species may not be recognizable.

Non-NASA Center↗

Determining divergence times with a protein clock: update and reevaluation

A recent study of the divergence times of the major groups of organisms as gauged by amino acid sequence comparison has been expanded and the data have been reanalyzed with a distance measure that corrects for both constraints on amino acid interchange and variation in substitution rate at different sites. Beyond that, the availability of complete genome sequences for several eubacteria and an archaebacterium has had a great impact on the interpretation of certain aspects of the data. Thus, the majority of the archaebacterial sequences are not consistent with currently accepted views of the Tree of Life which cluster the archaebacteria with eukaryotes. Instead, they are either outliers or mixed in with eubacterial orthologs. The simplest resolution of the problem is to postulate that many of these sequences were carried into eukaryotes by early eubacterial endosymbionts about 2 billion years ago, only very shortly after or even coincident with the divergence of eukaryotes and archaebacteria. The strong resemblances of these same enzymes among the major eubacterial groups suggest that the cyanobacteria and Gram-positive and Gram-negative eubacteria also diverged at about this same time, whereas the much greater differences between archaebacterial and eubacterial sequences indicate these two groups may have diverged between 3 and 4 billion years ago.

NASA Discipline Exobiology↗

CABO-16S—a Combined Archaea, Bacteria, Organelle 16S rRNA database framework for amplicon analysis of prokaryotes and eukaryotes in environmental samples

Abstract Identification of both prokaryotic and eukaryotic microorganisms in environmental samples is currently challenged by the need for additional sequencing to obtain separate 16S and 18S ribosomal RNA (rRNA) amplicons or the constraints imposed by “universal” primers. Organellar 16S rRNA sequences are amplified and sequenced along with prokaryote 16S rRNA and provide an alternative method to identify eukaryotic microorganisms. CABO-16S combines bacterial and archaeal sequences from the SILVA database with 16S rRNA sequences of plastids and other organelles from the PR2 database to enable identification of all 16S rRNA sequences. Comparison of CABO-16S with SILVA 138.2 results in equivalent taxonomic classification of mock communities and increased classification of diverse environmental samples. In particular, identification of phototrophic eukaryotes in shallow seagrass environments, marine waters, and lake waters was increased. The CABO-16S framework allows users to add custom sequences for further classification of underrepresented clades and can be easily updated with future releases of reference databases. Addition of sequences obtained from Sanger sequencing of methane seep sediments and curated sequences of the polyphyletic SEEP-SRB1 clade resulted in differentiation of syntrophic and non-syntrophic SEEP-SRB1 in hydrothermal vent sediments. CABO-16S highlights the benefit of combining and amending existing training sets when studying microorganisms in diverse environments.

Eitel, Eryn M. (ORCID:0009000723919297)↗

Epigenetic Footprints of CRISPR/Cas9-Mediated Genome Editing in Plants

CRISPR/Cas9 has been widely applied to various plant species accelerating the pace of plant genome editing and precision breeding in crops. Unintended effects beyond off-target nucleotide mutations are still somewhat unexplored. We investigated the degree and patterns of epigenetic changes after gene editing. We examined changes in DNA methylation in genome-edited promoters of naturally hypermethylated genes (AT1G72350 and AT1G09970) and hypomethylated genes (AT3G17320 and AT5G28770) from Arabidopsis. Transgenic plants were developed via Agrobacterium-mediated floral dip transformation. Homozygous edited lines were selected from segregated T2 plants by an in vitro digestion assay using ribonucleoprotein complex. Bisulfite sequencing comparisons were made between paired groups of edited and non-edited plants to identify changes in DNA methylation of the targeted loci. We found that directed mutagenesis via CRISPR/Cas9 resulted in no unintended morphological or epigenetic alterations. Phenotypes of wild-type, transgenic empty vector, and transgenic edited plants were similar. Epigenetic profiles revealed that methylation patterns of promoter regions flanking target sequences were identical among wild-type, transgenic empty vector, and transgenic edited plants. There was no effect of mutation type on epigenetic status. We also evaluated off-target mutagenesis effects in the edited plants. Potential off-target sites containing up to 4-bp mismatch of each target were sequenced. No off-target mutations were detected in candidate sites. Our results showed that CRISPR/Cas9 did not leave an epigenetic footprint on either the immediate gene-edited DNA and flanking DNA or introduce off-target mutations.

59 BASIC BIOLOGICAL SCIENCES↗

Genomic analysis and identification of a novel superantigen, SargEY, in Staphylococcus argenteus isolated from atopic dermatitis lesions

During surveillance of Staphylococcus aureus in lesions from patients with atopic dermatitis (AD), we isolated Staphylococcus argenteus, a species registered in 2011 as a new member of the genus Staphylococcus and previously considered a lineage of S. aureus. Genome sequence comparisons between S. argenteus isolates and representative S. aureus clinical isolates from various origins revealed that the S. argenteus genome from AD patients closely resembles that of S. aureus causing skin infections. We previously reported that 17%–22% of S. aureus isolated from skin infections produce staphylococcal enterotoxin Y (SEY), which predominantly induces T-cell proliferation via the T-cell receptor (TCR) Vα pathway. Complete genome sequencing of S. argenteus isolates revealed a gene encoding a protein similar to superantigen SEY, designated as SargEY, on its chromosome. Population structure analysis of S. argenteus revealed that these isolates are ST2250 lineage, which was the only lineage positive for the SEY-like gene among S. argenteus. Recombinant SargEY demonstrated immunological cross-reactivity with anti-SEY serum. SargEY could induce proliferation of human CD4 + and CD8 + T cells, as well as production of TNF-α and IFN-γ. SargEY showed emetic activity in a marmoset monkey model. S arg EY and SET (a phylogenetically close but uncharacterized SE) revealed their dependency on TCR Vα in inducing human T-cell proliferation. Additionally, TCR sequencing revealed other previously undescribed Vα repertoires induced by SEH. S arg EY and SEY may play roles in exacerbating the respective toxin-producing strains in AD.

59 BASIC BIOLOGICAL SCIENCES↗

Were the original eubacteria thermophiles?

Thermotoga maritima is one of the more unusual eubacteria: It is highly thermophilic, growing at temperatures higher than any other eubacterium; its cell wall appears to have a unique structure and its lipids a unique composition; and the organism is surrounded by a loose-fitting sheath of unknown function. Its phenotypic uniqueness is matched by its phylogenetic position; Thermotoga maritima represents the deepest known branching in the eubacterial line of descent, as measured by ribosomal RNA sequence comparisons. T. maritima also represents the most slowly evolving of eubacterial lineages. The fact that the two deepest branchings in the eubacterial line of descent (the other, the green non-sulfur bacteria and relatives, i.e. Chloroflexus, Thermomicrobium, etc.) are both basically thermophilic and slowly evolving, strongly suggests that all eubacteria have ultimately arisen from a thermophilic ancestor.

NASA Discipline Exobiology↗

Transcription in archaea

Using the sequences of all the known transcription-associated proteins from Bacteria and Eucarya (a total of 4,147), we have identified their homologous counterparts in the four complete archaeal genomes. Through extensive sequence comparisons, we establish the presence of 280 predicted transcription factors or transcription-associated proteins in the four archaeal genomes, of which 168 have homologs only in Bacteria, 51 have homologs only in Eucarya, and the remaining 61 have homologs in both phylogenetic domains. Although bacterial and eukaryotic transcription have very few factors in common, each exclusively shares a significantly greater number with the Archaea, especially the Bacteria. This last fact contrasts with the obvious close relationship between the archaeal and eukaryotic transcription mechanisms per se, and in particular, basic transcription initiation. We interpret these results to mean that the archaeal transcription system has retained more ancestral characteristics than have the transcription mechanisms in either of the other two domains.

NASA Discipline Exobiology↗

Calcium-stimulated autophosphorylation site of plant chimeric calcium/calmodulin-dependent protein kinase

The existence of two molecular switches regulating plant chimeric Ca(2+)/calmodulin-dependent protein kinase (CCaMK), namely the C-terminal visinin-like domain acting as Ca(2+)-sensitive molecular switch and calmodulin binding domain acting as Ca(2+)-stimulated autophosphorylation-sensitive molecular switch, has been described (Sathyanarayanan, P. V., Cremo, C. R., and Poovaiah, B. W. (2000) J. Biol. Chem. 275, 30417-30422). Here we report the identification of Ca(2+)-stimulated autophosphorylation site of CCaMK by matrix-assisted laser desorption ionization time of flight-mass spectrometry. Thr(267) was confirmed as the Ca(2+)-stimulated autophosphorylation site by post-source decay experiments and by site-directed mutagenesis. The purified T267A mutant form of CCaMK did not show Ca(2+)-stimulated autophosphorylation, autophosphorylation-dependent variable calmodulin affinity, or Ca(2+)/calmodulin stimulation of kinase activity. Sequence comparison of CCaMK from monocotyledonous plant (lily) and dicotyledonous plant (tobacco) suggests that the autophosphorylation site is conserved. This is the first identification of a phosphorylation site specifically responding to activation by second messenger system (Ca(2+) messenger system) in plants. Homology modeling of the kinase and calmodulin binding domain of CCaMK with the crystal structure of calcium/calmodulin-dependent protein kinase 1 suggests that the Ca(2+)-stimulated autophosphorylation site is located on the surface of the kinase and far from the catalytic site. Analysis of Ca(2+)-stimulated autophosphorylation with increasing concentration of CCaMK indicates the possibility that the Ca(2+)-stimulated phosphorylation occurs by an intermolecular mechanism.

Non-NASA Center↗

The green non-sulfur bacteria: a deep branching in the eubacterial line of descent

Ribosomal RNA sequence comparisons define a phylogenetic grouping, the green non-sulfur bacteria and relatives (GNS), known to contain the genera Chloroflexus, Herpetosiphon and Thermomicrobium--organisms that have little phenotypic similarity. The unit is phylogenetically deep, but entirely distinct from any other eubacterial division (phylum). It is also relatively ancient--branching from the common eubacterial stem earlier than any other group of eubacteria reported thus far. The group phenotype is predominantly thermophilic, and its thermophilic members, especially Thermomicrobium, are more slowly evolving than Herpetosiphon, a mesophile. The GNS unit appears significantly older than either the green sulfur bacteria or the cyanobacteria--making it likely that organisms such as Chloroflexus, not the cyanobacteria, generated the oldest stromatolites, which formed over three billion years ago.

Non-NASA Center↗

Insights revealed by the co-crystal structure of the Saccharomyces cerevisiae histidine phosphotransfer protein Ypd1 and the receiver domain of its downstream response regulator Ssk1

Two-component signaling systems are the primary means by which bacteria, archaea, and certain plants and fungi react to their environments. The model yeast, Saccharomyces cerevisiae, uses the Sln1 signaling pathway to respond to hyperosmotic stress. This pathway contains a hybrid histidine kinase (Sln1) that autophosphorylates and transfers a phosphoryl group to its own receiver domain (R1). The phosphoryl group is then transferred to a histidine phosphotransfer protein (Ypd1) that finally passes it to the receiver domain (R2) of a downstream response regulator (Ssk1). Under normal conditions, Ssk1 is constitutively and preferentially phosphorylated in the phosphorelay. Upon detecting hyperosmotic stress, Ssk1 rapidly dephosphorylates and activates the high-osmolarity glycerol (HOG) pathway, initiating a response. Despite their distinct physiological roles, both Sln1 and Ssk1 bind to Ypd1 at a common docking site. Co-crystal structures of response regulators in complex with their phosphorelay partners are scarce, leaving many mechanistic and structural details uncharacterized for systems like the Sln1 pathway. In this work, we present the co-crystal structure of Ypd1 and a near wild-type variant of the receiver domain of Ssk1 (Ssk1-R2-W638A) at a resolution of 2.80 Å. Our structural analyses of Ypd1-receiver domain complexes, biochemical determination of binding affinities for Ssk1-R2 variants, in silico free energy estimates, and sequence comparisons reveal distinctive electrostatic properties of the Ypd1/ Ssk1-R2-W638A complex that may provide insight into the regulation of the Sln1 pathway as a function of dynamic osmolyte concentration.

59 BASIC BIOLOGICAL SCIENCES↗