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 181 records · Page 10

Sequence information signal processor

An electronic circuit is used to compare two sequences, such as genetic sequences, to determine which alignment of the sequences produces the greatest similarity. The circuit includes a linear array of series-connected processors, each of which stores a single element from one of the sequences and compares that element with each successive element in the other sequence. For each comparison, the processor generates a scoring parameter that indicates which segment ending at those two elements produces the greatest degree of similarity between the sequences. The processor uses the scoring parameter to generate a similar scoring parameter for a comparison between the stored element and the next successive element from the other sequence. The processor also delivers the scoring parameter to the next processor in the array for use in generating a similar scoring parameter for another pair of elements. The electronic circuit determines which processor and alignment of the sequences produce the scoring parameter with the highest value.

Peterson, John C.↗

Microbial Dark Matter: Unusual intervening sequences in 16S rRNA genes of candidate phyla from the deep subsurface

The Microbial Dark Matter project has sequenced genomes from over 200 single cells from candidate phyla, greatly expanding our knowledge of the ecology, inferred metabolism, and evolution of these widely distributed, yet poorly understood lineages. The second phase of this project aims to sequence an additional 800 single cells from known as well as potentially novel candidate phyla derived from a variety of environments. In order to identify whole genome amplified single cells, screening based on phylogenetic placement of 16S rRNA gene sequences is being conducted. Briefly, derived 16S rRNA gene sequences are aligned to a custom version of the Greengenes reference database and added to a reference tree in ARB using parsimony. In multiple samples from deep subsurface habitats but not from other habitats, a large number of sequences proved difficult to align and therefore to place in the tree. Based on comparisons to reference sequences and structural alignments using SSU-ALIGN, many of these ?difficult? sequences appear to originate from candidate phyla, and contain intervening sequences (IVSs) within the 16S rRNA genes. These IVSs are short (39 - 79 nt) and do not appear to be self-splicing or to contain open reading frames. IVSs were found in the loop regions of stem-loop structures in several different taxonomic groups. Phylogenetic placement of sequences is strongly affected by IVSs; two out of three groups investigated were classified as different phyla after their removal. Based on data from samples screened in this project, IVSs appear to be more common in microbes occurring in deep subsurface habitats, although the reasons for this remain elusive.

09 BIOMASS FUELS↗

CRITICA: coding region identification tool invoking comparative analysis

Gene recognition is essential to understanding existing and future DNA sequence data. CRITICA (Coding Region Identification Tool Invoking Comparative Analysis) is a suite of programs for identifying likely protein-coding sequences in DNA by combining comparative analysis of DNA sequences with more common noncomparative methods. In the comparative component of the analysis, regions of DNA are aligned with related sequences from the DNA databases; if the translation of the aligned sequences has greater amino acid identity than expected for the observed percentage nucleotide identity, this is interpreted as evidence for coding. CRITICA also incorporates noncomparative information derived from the relative frequencies of hexanucleotides in coding frames versus other contexts (i.e., dicodon bias). The dicodon usage information is derived by iterative analysis of the data, such that CRITICA is not dependent on the existence or accuracy of coding sequence annotations in the databases. This independence makes the method particularly well suited for the analysis of novel genomes. CRITICA was tested by analyzing the available Salmonella typhimurium DNA sequences. Its predictions were compared with the DNA sequence annotations and with the predictions of GenMark. CRITICA proved to be more accurate than GenMark, and moreover, many of its predictions that would seem to be errors instead reflect problems in the sequence databases. The source code of CRITICA is freely available by anonymous FTP (rdp.life.uiuc.edu in/pub/critica) and on the World Wide Web (http:/(/)rdpwww.life.uiuc.edu).

Non-NASA Center↗

Evolution, language and analogy in functional genomics

Almost a century ago, Wittgenstein pointed out that theory in science is intricately connected to language. This connection is not a frequent topic in the genomics literature. But a case can be made that functional genomics is today hindered by the paradoxes that Wittgenstein identified. If this is true, until these paradoxes are recognized and addressed, functional genomics will continue to be limited in its ability to extrapolate information from genomic sequences.

Evolution, Molecular↗

An early ethylene up-regulated gene encoding a calmodulin-binding protein involved in plant senescence and death

35S-Labeled calmodulin (CaM) was used to screen a tobacco anther cDNA library. A positive clone (NtER1) with high homology to an early ethylene-up-regulated gene (ER66) in tomato, and an Arabidopsis homolog was isolated and characterized. Based on the helical wheel projection, a 25-mer peptide corresponding to the predicted CaM-binding region of NtER1 (amino acids 796-820) was synthesized. The gel-mobility shift assay showed that the peptide formed a stable complex with CaM only in the presence of Ca(2+). CaM binds to NtER1 with high affinity (K(d) approximately 12 nm) in a calcium-dependent manner. Tobacco flowers at different stages of development were treated with ethylene or with 1-methylcyclopropene for 2 h before treating with ethylene. Northern analysis showed that the NtER1 was rapidly induced after 15 min of exposure to ethylene. However, the 2-h 1-methylcyclopropene treatment totally blocked NtER1 expression in flowers at all stages of development, suggesting that NtER1 is an early ethylene-up-regulated gene. The senescing leaves and petals had significantly increased NtER1 induction as compared with young leaves and petals, implying that NtER1 is developmentally regulated and acts as a trigger for senescence and death. This is the first documented evidence for the involvement of Ca(2+)/CaM-mediated signaling in ethylene action.

Non-NASA Center↗

Identification of a new EF-hand superfamily member from Trypanosoma brucei

We identified several open reading frames between the regions encoding calmodulin and ubiquitin-EP52/1 in the genome of Trypanosoma brucei. One of these, EFH5, encodes a protein 192 amino acids long. The EFH5 transcript is present in poly(A)+ mRNA and is present at similar levels in the mammalian bloodstream form and the insect procyclic form. EFH5 contains four EF-hand homolog domains, two of which are inferred to bind Ca2+ ions. We expressed EFH5 as a fusion protein in Escherichia coli and demonstrated calcium-binding activity of the fusion protein using the 45Ca-overlay technique. The function of EFH5 remains unknown; however, as the fourth EF-hand homolog identified in trypanosomes, it attests to the broad range of functions assumed by calcium functioning as a second messenger. EFH5, which is most closely related to LAV1-2 from Physarum, represents a distinct subfamily among the EF-hand-containing proteins.

Non-NASA Center↗

Generate Optimized Genetic Rhythm for Enzyme Expression in Non-native systems

Most amino acids are represented by more than one codon, resulting in redundancy in the genetic code. Silent codon substitutions that do not alter the amino acid sequence still have an effect on protein expression. We have developed an algorithm, GoGREEN, to enhance the expression of foreign proteins in a host organism. GoGREEN selects codons according to frequency patterns seen in the gene of interest using the codon usage table from the host organism. GoGREEN is also designed to accommodate gaps in the sequence.This software takes for input (1) the aligned protein sequences for genes the user wishes to express, (2) the codon usage table for the host organism, (3) and the DNA sequence for the target protein found in the host organism. The program will select codons based on codon usage patterns for the target DNA sequence. The program will also select codons for “gaps” found in the aligned protein sequences using the codon usage table from the host organism.

Source record↗

Protein remote homology detection and structural alignment using deep learning

Exploiting sequence–structure–function relationships in biotechnology requires improved methods for aligning proteins that have low sequence similarity to previously annotated proteins. We develop two deep learning methods to address this gap, TM-Vec and DeepBLAST. TM-Vec allows searching for structure–structure similarities in large sequence databases. It is trained to accurately predict TM-scores as a metric of structural similarity directly from sequence pairs without the need for intermediate computation or solution of structures. Once structurally similar proteins have been identified, DeepBLAST can structurally align proteins using only sequence information by identifying structurally homologous regions between proteins. It outperforms traditional sequence alignment methods and performs similarly to structure-based alignment methods. We show the merits of TM-Vec and DeepBLAST on a variety of datasets, including better identification of remotely homologous proteins compared with state-of-the-art sequence alignment and structure prediction methods.

59 BASIC BIOLOGICAL SCIENCES↗

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↗

ModBase, a database of annotated comparative protein structure models, and associated resources

ModBase is a database of annotated comparative protein structure models. The models are calculated by ModPipe, an automated modeling pipeline that relies primarily on Modeller for fold assignment, sequence–structure alignment, model building and model assessment. ModBase currently contains 10 355 444 reliable models for domains in 2 421 920 unique protein sequences. ModBase allows users to update comparative models on demand, and request modeling of additional sequences through an interface to the ModWeb modeling server. ModBase models are available through the ModBase interface as well as the Protein Model Portal. Recently developed associated resources include the SALIGN server for multiple sequence and structure alignment, the ModEval server for predicting the accuracy of protein structure models, the PCSS server for predicting which peptides bind to a given protein and the FoXS server for calculating and fitting Small Angle X-ray Scattering profiles.

59 BASIC BIOLOGICAL SCIENCES↗

Amplicon Sequencing Reveals Microbiological Signatures in Spent Nuclear Fuel Storage Basins

Water quality is an important determinant for the structural integrity of alloy cladded fuels and assemblies during long-term wet storage. Detailed characterization of a water filled storage basin for spent nuclear reactor fuel was performed following the formation and proliferation of an amorphous white flocculent. White precipitant was sampled throughout the storage basin for chemical and spectroscopic characterization, and eDNA was extracted for pyrosequencing of bacterial rRNA gene diversity. Accordingly, spectroscopic analyses indicated the precipitant to be primarily amorphous to crystalline aluminum (oxy) hydroxides with minor associated elemental components including Fe, Si, Ti, and U. High levels of dissolved carbon were co-localized with the precipitant relative to bulk water. Bacterial densities were highly variable between sampling locations and with depth; cell numbers (log scale) ranged from 5.6 to 4.89 cells / mL. Bacterial diversity that was physically associated with the aluminum (oxy) hydroxide complexes exceeded an estimated 4,000 OTUs / amplicon library (3% cutoff) and the greatest percent majority of sequences were aligned to the families Burkholderiales (23%), Nitrospiraceae (23%), Hyphomicrobiaceae (17%), and Comamonadaceae (6%). We surmise that episodic changes in the physical and chemical properties of the basin contribute to the polymerization of aluminum (oxy) hydroxides, which in turn can chemisorb nutrients, carbon ligands and bacterial cells from the surrounding bulk aqueous phase. As such, these precipitants should establish favorable microhabitats for bacterial colonization and growth. Comparative analyses of 16S rRNA gene amplicon libraries across diverse environmental landscapes were performed and microbiological signatures unique to the spent nuclear fuel storage basin environment were revealed. These insights could spur the development of tractable bioindicators that are specific of and diagnostic for water quality at discrete locations and finer scales of resolution, marking an important contribution for improved water quality and management of spent nuclear fuel storage facilities.

aluminum↗

Mitochondrial gene arrangement of the horseshoe crab Limulus polyphemus L.: conservation of major features among arthropod classes

Numerous complete mitochondrial DNA sequences have been determined for species within two arthropod groups, insects and crustaceans, but there are none for a third, the chelicerates. Most mitochondrial gene arrangements reported for crustaceans and insect species are identical or nearly identical to that of Drosophila yakuba. Sequences across 36 of the gene boundaries in the mitochondrial DNA (mtDNA) of a representative chelicerate. Limulus polyphemus L., also reveal an arrangement like that of Drosophila yakuba. Only the position of the tRNA(LEU)(UUR) gene differs; in Limulus it is between the genes for tRNA(LEU)(CUN) and ND1. This positioning is also found in onychophorans, mollusks, and annelids, but not in insects and crustaceans, and indicates that tRNA(LEU)(CUN)-tRNA(LEU)(UUR)-ND1 was the ancestral gene arrangement for these groups, as suggested earlier. There are no differences in the relative arrangements of protein-coding and ribosomal RNA genes between Limulus and Drosophila, and none have been observed within arthropods. The high degree of similarity of mitochondrial gene arrangements within arthropods is striking, since some taxa last shared a common ancestor before the Cambrian, and contrasts with the extensive mtDNA rearrangements occasionally observed within some other metazoan phyla (e.g., mollusks and nematodes).

Non-NASA Center↗

A high-density SNP genetic map consisting of a complete set of homologous groups in autohexaploid sweetpotato (Ipomoea batatas)

Sweetpotato (Ipomoea batatas) is an autohexaploid species with 90 chromosomes (2n = 6x = 90) and a basic chromosome number of 15, and is therefore regarded as one of the most challenging species for high-density genetic map construction. Here, we used single nucleotide polymorphisms (SNPs) identified by double-digest restriction site-associated DNA sequencing based on next-generation sequencing technology to construct a map for sweetpotato. We then aligned the sequence reads onto the reference genome sequence of I. trifida, a likely diploid ancestor of sweetpotato, to detect SNPs. In addition, to simplify analysis of the complex genetic mode of autohexaploidy, we used an S1 mapping population derived from self-pollination of a single parent. As a result, 28,087 double-simplex SNPs showing a Mendelian segregation ratio in the S1 progeny could be mapped onto 96 linkage groups (LGs), covering a total distance of 33,020.4 cM. Based on the positions of the SNPs on the I. trifida genome, the LGs were classified into 15 groups, each with roughly six LGs and six small extra groups. The molecular genetic techniques used in this study are applicable to high-density mapping of other polyploid plant species, including important crops.

59 BASIC BIOLOGICAL SCIENCES↗

Two genes with similarity to bacterial response regulators are rapidly and specifically induced by cytokinin in Arabidopsis

Cytokinins are central regulators of plant growth and development, but little is known about their mode of action. By using differential display, we identified a gene, IBC6 (for induced by cytokinin), from etiolated Arabidopsis seedlings, that is induced rapidly by cytokinin. The steady state level of IBC6 mRNA was elevated within 10 min by the exogenous application of cytokinin, and this induction did not require de novo protein synthesis. IBC6 was not induced by other plant hormones or by light. A second Arabidopsis gene with a sequence highly similar to IBC6 was identified. This IBC7 gene also was induced by cytokinin, although with somewhat slower kinetics and to a lesser extent. The pattern of expression of the two genes was similar, with higher expression in leaves, rachises, and flowers and lower transcript levels in roots and siliques. Sequence analysis revealed that IBC6 and IBC7 are similar to the receiver domain of bacterial two-component response regulators. This homology, coupled with previously published work on the CKI1 histidine kinase homolog, suggests that these proteins may play a role in early cytokinin signaling.

Non-NASA Center↗

Isolation, characterization, and amino acid sequences of auracyanins, blue copper proteins from the green photosynthetic bacterium Chloroflexus aurantiacus

Three small blue copper proteins designated auracyanin A, auracyanin B-1, and auracyanin B-2 have been isolated from the thermophilic green gliding photosynthetic bacterium Chloroflexus aurantiacus. All three auracyanins are peripheral membrane proteins. Auracyanin A was described previously (Trost, J. T., McManus, J. D., Freeman, J. C., Ramakrishna, B. L., and Blankenship, R. E. (1988) Biochemistry 27, 7858-7863) and is not glycosylated. The two B forms are glycoproteins and have almost identical properties to each other, but are distinct from the A form. The sodium dodecyl sulfate-polyacrylamide gel electrophoresis apparent monomer molecular masses are 14 (A), 18 (B-2), and 22 (B-1) kDa. The amino acid sequences of the B forms are presented. All three proteins have similar absorbance, circular dichroism, and resonance Raman spectra, but the electron spin resonance signals are quite different. Laser flash photolysis kinetic analysis of the reactions of the three forms of auracyanin with lumiflavin and flavin mononucleotide semiquinones indicates that the site of electron transfer is negatively charged and has an accessibility similar to that found in other blue copper proteins. Copper analysis indicates that all three proteins contain 1 mol of copper per mol of protein. All three auracyanins exhibit a midpoint redox potential of +240 mV. Light-induced absorbance changes and electron spin resonance signals suggest that auracyanin A may play a role in photosynthetic electron transfer. Kinetic data indicate that all three proteins can donate electrons to cytochrome c-554, the electron donor to the photosynthetic reaction center.

NASA Discipline Exobiology↗

The LANL hemorrhagic fever virus database, a new platform for analyzing biothreat viruses

Hemorrhagic fever viruses (HFVs) are a diverse set of over 80 viral species, found in 10 different genera comprising five different families: arena-, bunya-, flavi-, filo- and togaviridae. All these viruses are highly variable and evolve rapidly, making them elusive targets for the immune system and for vaccine and drug design. About 55 000 HFV sequences exist in the public domain today. A central website that provides annotated sequences and analysis tools will be helpful to HFV researchers worldwide. The HFV sequence database collects and stores sequence data and provides a user-friendly search interface and a large number of sequence analysis tools, following the model of the highly regarded and widely used Los Alamos HIV database [Kuiken, C., B. Korber, and R.W. Shafer, HIV sequence databases. AIDS Rev, 2003. 5: p. 52–61]. The database uses an algorithm that aligns each sequence to a species-wide reference sequence. The NCBI RefSeq database [Sayers et al. (2011) Database resources of the National Center for Biotechnology Information. Nucleic Acids Res., 39, D38–D51.] is used for this; if a reference sequence is not available, a Blast search finds the best candidate. Using this method, sequences in each genus can be retrieved pre-aligned. The HFV website can be accessed via http://hfv.lanl.gov.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Characterization of Two Microbial Isolates from Andean Lakes in Bolivia

We are currently investigating the biological population present in the highest and least explored perennial lakes on earth in the Bolivian and Chilean Andes, including several volcanic crater lakes of more than 6000 m elevation, in combination of microbiological and molecular biological methods. Our samples were collected in saline lakes of the Laguna Blanca Laguna Verde area in the Bolivian Altiplano and in the Licancabur volcano crater (27 deg. 47 min S/67 deg. 47 min. W) in the ongoing project studying high altitude lakes. The main goal of the project is to look for analogies with Martian paleolakes. These Bolivian lakes can be described as Andean lakes following the classification of Chong. We have attempted to isolate pure cultures and phylogenetically characterize prokaryotes that grew under laboratory conditions. Sediment samples taken from the Licancabur crater lake (LC), Laguna Verde (LV), and Laguna Blanca (LB) were analyzed and cultured using enriched liquid media under both aerobic and anaerobic conditions. All cultures were incubated at room temperature (15 to 20 C) and under light exposure. For the reported isolates, 36 hours incubation were necessary for reaching optimal optical densities to consider them viable cultures. Ten serial dilutions starting from 1% inoculum were required to obtain a suitable enriched cell culture to transfer into solid media. Cultures on solid medium were necessary to verify the formation of colonies in order to isolate pure cultures. Different solid media were prepared using several combinations of both trace minerals and carbohydrates sources in order to fit their nutrient requirements. The microorganisms formed individual colonies on solid media enriched with tryptone, yeast extract and sodium chloride. Cells morphology was studied by optical and electronic microscopy. Rodshape morphologies were observed in most cases. Total bacterial genomic DNA was isolated from 50 ml late-exponential phase culture by using the CTAB miniprep protocol. The 16S rRNA genes were amplified by PCR using both Bacteria- and Archaeauniversal primer sets: 27f and 1492r, 21f and 1492r respectively. Sequences of 16S rRNA gene were determined and initially compared with reference sequences contained in the EMBL nucleotide sequence database by using the BLAST program and were subsequently aligned with 16S rRNA reference sequences in the ARB package (http://www.mikro.biologie.tu-muenchen.de). Aligned sequences were inserted within a stable phylogenetic tree by using the ARB parsimony tool. In this work we report the morphology and phylogenetic characterization of two isolates belonged to Laguna Blanca sediments.

Demergasso, C.↗